Interferon regulatory factor 7 regulates airway epithelial cell responses to human rhinovirus infection
BMC Genomics volume 17, Article number: 76 (2016)
Human rhinoviruses (HRV) cause the majority of colds and trigger exacerbations of chronic lower airway diseases. Airway epithelial cells are the primary site for HRV infection and replication, and the initiation of host inflammatory responses. At present, the molecular mechanisms that underpin HRV responses in airway epithelial cells are incompletely understood. The aim of this study was to employ microarray profiling, upstream regulator analysis, and siRNA mediated gene silencing to further our understanding of the role of interferon regulatory factor 7 (IRF7) in this response.
Primary human bronchial epithelial cells (HBE) where transfected with siRNA that targets IRF7 or a non-silencing control (all-star control) using Lipofectamine. The cells were allowed to recover, and then cultured in the presence or absence of HRV-16 for 24 h. Global patterns of gene expression were profiled on microarrays. A subset of genes identified in the microarray study were validated at the mRNA and/or protein level using real time RT-qPCR, ELISA, and western blots.
Hundreds of genes were upregulated in HBE during HRV infection. Pathways analysis demonstrated that these genes were mainly involved in type I and II interferon signaling, RIG-I/MDA5 signaling, antigen processing and presentation, and apoptosis. Upstream regulator analysis of these data suggested that IRF7 was a major molecular driver of this response. Knockdown of IRF7 reduced the HRV-driven upregulation of genes involved in antiviral responses (interferon signaling, Toll-like receptor signaling, NOD-like receptor signaling, RIG-I/MDA5 signaling), and increased the expression of genes that promote inflammation (e.g. CXCL5, IL-33, IL1RL1) and the response to oxidative stress. However, the majority of genes that were perturbed by HRV in HBE cells including those that are known to be regulated by IRF7 were insensitive to IRF7 knockdown. Upstream regulator analysis of the part of the response that was insensitive to IRF7 knockdown suggested it was driven by NF-κB, STAT1, STAT3, and IRF1.
Our findings demonstrate that IRF7 regulates the expression of genes involved in antiviral immunity, inflammation, and the response to oxidative stress during HRV infections in HBE cells, and also suggests that other transcription factors play a major role in this response.
Human rhinovirus is a positive sense, single-strand RNA virus from the Picornavirus family. It causes the majority of colds, and it also triggers exacerbations of asthma, chronic obstructive pulmonary disease, and cystic fibrosis . More than 150 strains of HRV have been identified, and these have been classified into three species (A, B, C) based on genome sequencing data . Airway epithelial cells are the primary site for rhinovirus infection and replication. Depending on the serotype of HRV, cell entry is mediated by binding to the receptors ICAM-1 (major group), low-density-lipoprotein receptor family members (minor group), or cadherin-related family member 3 (CDHR3 for HRV-C) .
Host responses to HRV are thought to be initiated when viral proteins and nucleic acids trigger pathogen recognition receptors (PRR) of the innate immune system (TLR2, TLR3, TLR7, TLR8, MDA5, RIG-I) [4–6]. PRR signaling triggers intracellular signaling cascades that converge on the IRF and NF-κB family of transcription factors, which in turn activate interferon-induced antiviral and proinflammatory gene expression programs respectively. In epithelial cells, interferon programs are mainly activated via TLR3, MDA5, and RIG-I signaling [4–6]. IRF7 plays a central role in the activation of innate antiviral responses by controlling the transcription of type I and III interferon genes . Type I interferon signaling in turn upregulates IRF7 expression, thereby providing positive feedback amplification of the antiviral response . The interferon system plays a crucial role in antiviral defence by upregulating a set of effector molecules (e.g. Mx1, OAS, PKR, ISG15), which induce a robust antiviral state in infected and neighbouring cells preventing the spread of infection . The importance of IRF7 in antiviral immunity was clearly demonstrated in IRF7 deficient mice, which are highly susceptible to viral infections . Recent data in humans indicates that null mutations in IRF7 are associated with susceptibility to severe, life-threatening influenza .
Although it is well established that IRF7 is a master regulator of the antiviral response, the function (or redundancy thereof) of IRF7 can vary markedly depending on the cell type and the virus . We reported that expression of IRF7 was upregulated in nasal epithelial scrapings from human adults after experimental infection with HRV-16 . Employing a computational analysis of gene expression profiling data derived from nasal wash samples, we also showed that IRF7 was a major hub, connecting interferon-related gene networks during picornavirus-induced asthma exacerbations in children . Lewis et al. reported that expression of IRF7 was upregulated in nasal aspirates from asthmatic children during natural colds . Notably, all of the above studies focused on expression levels of IRF7; this is the first study to perform gene silencing experiments in HBE to elucidate the role of IRF7 in the regulation of HRV responses.
The following reagents were purchased from the indicated suppliers: Ham’s F-12 medium, Eagle’s minimal essential medium, Hank’s balanced salt solution (HBSS), penicillin-streptomycin-amphotericin B, L-glutamine, TRIzol reagent, sodium pyruvate, nonessential amino acids, gentamicin, fetal bovine serum (FBS), dNTPs, oligo(dT), random hexamers, and Superscript III from Invitrogen Life Technologies (Burlington, Ontario, Canada). Bronchial epithelial cell growth medium (BEGM) was from BioWhittaker (Walkersville, MD). Gene expression kits for GAPDH, IL-6, CXCL5, and HERC5 were from Applied Biosystems (Foster City, CA); antibody pairs for CCL5, CXCL5, and CXCL10, as well as recombinant standards for each assay were from R&D Systems (Minneapolis, MN); TaqMan master mix and complete protease inhibitor cocktail tablets were from Roche Diagnostics (Laval, Quebec, Canada); western blot antibody for IRF7 was from BioLegend (San Diego, CA), while the viperin antibody was from Abcam (Toronto, ON). All other chemicals were purchased from Sigma-Aldrich (St. Louis, MO).
Primary epithelial cell culture
Normal human lungs not used for transplantation were obtained from a tissue retrieval service (International Institute for the Advancement of Medicine, Edison, NJ). Primary HBE cells from a total of 9 donors (4 male; age range 13–63 years) were used for the current studies. All donors died from head trauma or stroke, and none of the donors had any inflammatory lung disease. Written informed consent was obtained from responsible next of kin for tissues to be used for either transplant or for research. Approval to use recovered organs for these studies was obtained from the Conjoint Health Research Ethics Board of the University of Calgary. Primary HBE cells were obtained by protease digestion of dissected airways as previously described . HBE cells were grown on six-well culture plates in BEGM. Hydrocortisone was removed from the medium 24 h prior to stimulation and all stimulations were performed in this hydrocortisone-free medium.
WI-38 cells were purchased from the American Type Culture Collection (Manassas, VA). HRV-16 viral stocks were propagated in WI-38 cells and purified by centrifugation through sucrose to remove ribosomes and soluble factors of WI-38 origin as previously described . HBE cells were infected with 105.5 50 % tissue culture-infective dose (TCID50) U/ml, which equates (approximately) to a multiplicity of infection of 1.
RNA extraction and real-time RT-qPCR
Total cellular RNA was extracted from cultured epithelial cells with TRIzol (1 ml per 10 cm2) as previously described . Yield of RNA was quantified based upon absorbance at 260 nm. RNA was DNase treated with DNase I (Ambion, Austin, TX) and then requantified. One μg of RNA was reverse transcribed with both Oligo dT and random hexamers with Superscript III, followed by PCR amplification in the presence of specific forward and reverse primers and fluorescently labeled probes for each gene of interest. Gene expression for CXCL10, CCL5, Viperin, IL-6, CXCL5, and HERC5 were quantified using the Applied Biosystems Model 7900 Sequence Detector (Foster City, CA). Analysis of the housekeeping gene GAPDH was performed on each sample using a primer and probe kit obtained from Applied Biosystems. Fold induction of each gene was calculated using the formula 2-ΔΔCT as previously described . Primers and probes for CXCL10, CCL5 and viperin were as follows: CXCL10 forward primer 5′-GAAATTATTCCTGCAAGCCAATTT-3′; reverse primer 5′-TCACCCTTCTTTTTCATTGTAGCA-3′; probe 5′-FAM-TCCACGTGTTGAGATCA-MGB-3′, CCL5 forward primer 5′-TCTGCGCTCCTGCATCTG-3′; reverse primer 5′-AGTGGGCGGGCAATGTAG-3′; probe 5′-FAM-ATTCCTCGGACACCACACCCTGCTG-MGB-3′, and viperin forward primer 5′-CCTGCTTGGTGCCTGAATCT-3′; reverse primer 5′- GCGCATATATTCATCCAGAATAAGG-3′; probe 5′-FAM-ACCAGAAGATGAAAGACT-MGB-3′. Gene expression kits were purchased from Applied Biosystems for IL-6 (Hs00985639_m1), CXCL5 (Hs01099660_g1), and HERC5 (Hs00180943_m1).
CXCL10, CCL5 and CXCL5 protein levels were assayed by enzyme-linked immunosorbent assay using matched antibody pairs (R&D Systems, Minneapolis, MN).
siRNA knockdown of IRF7
Subconfluent HBE cells were transfected with either 10 nM of a siRNA targeted to IRF7 or to a nontargeting siRNA, all-star control (Qiagen, Toronto, ON) for 24 h at 37 °C using Lipofectamine RNAiMAX (Life Technologies) in BEGM without antibiotics. After transfection, medium was changed to BEGM without hydrocortisone and cells were allowed to recover for 24 h. Cells were then infected with HRV-16, and supernatants, whole-cell lysates and total RNA were collected 24 h post infection. The specific forward siRNA sequences used were as follows: IRF7 siRNA #1, 5′-CCCGAGCTGCACGTTCCTATA-3′; IRF7 siRNA #2, 5′-CTGGAAGCACTTCGCGCGCAA-3′.
After stimulation, cells were washed with HBSS and then scraped in lysis buffer (1 % Triton X-100, 1X mini complete, 1 mM PMSF, 2 mM sodium orthovanadate, 20 mM sodium pyrophosphate, and 50 mM sodium fluoride). Lysates were incubated on ice for 10 min, frozen at −80 °C, thawed, sonicated, and centrifuged at 10,000 × g for 5 min at 4 °C. Protein concentrations for the Triton X-100-soluble lysates were quantified using a DC Protein Assay (Bio-Rad, Montreal, Quebec, Canada). Equivalent amounts of protein were separated by SDS-PAGE, and proteins were then transferred to a nitrocellulose membrane. Membranes were blocked with 5 % skim milk for 1 h and probed with either 1 μg/mL of the IRF7 antibody or a 1:1000 dilution of the viperin antibody overnight at 4 °C. Membranes were washed and then incubated for 1 h with a 1:10,000 dilution of horseradish peroxidase-conjugated anti-mouse Ig antibody (Jackson ImmunoResearch Laboratories, Inc., West Grove, PA). Proteins were visualized with a Pierce ECL substrate reagent (Thermo Scientific, Rockford, IL). Equal loading was further determined probing with an antibody against GAPDH.
Gene expression profiling
Total RNA samples were shipped on dry ice to Expression Analysis Inc. (Durham, NC). Purified RNA was quantified using a Nanodrop spectrophotometer and quality was determined using the Agilent (Santa Clara, CA) 2100 Lab-on-a-Chip System. The RNA samples were processed using the Enzo Single Round RNA Amplification and Biotin Labeling System (Enzo, New York, NY), and hybridized to Primeview microarrays (Affymetrix, Santa Clara CA). The raw microarray data is available at the Gene Expression Omnibus repository under accession GSE70190.
The microarray data was analyzed in R software for statistical computing (https://www.r-project.org/). The quality of the microarray data was assessed with ArrayQualityMetrics . The data was preprocessed employing the RMA algorithm . A custom chip description file (primeviewhsentrezg, version 19) was used to map probe sets to genes . The data was filtered to differentiate between relevant gene expression signals and noise . Differentially expressed genes were identified using Linear models for microarray data (LIMMA), with False Discovery Rate (FDR) control for multiple testing . Pathways analysis was performed using InnateDb. Upstream regulators were identified using Ingenuity Systems software . The IRF7 gene network was reconstructed using experimentally supported findings from the Ingenuity Systems Knowledgebase .
For assessments of mRNA and protein expression for individual genes of interest, normally distributed data were analyzed using one-way analysis of variance (ANOVA) with student Newman-Keuls post hoc analysis. Analysis of paired comparisons were performed using t-tests. Data that were not normally distributed were analyzed using Kruskal-Wallis ANOVA followed by Wilcoxon matched-pairs signed-rank test. Statistical hypotheses were two-sided and values of p < 0.05 were considered significant.
Rhinovirus-induced gene expression patterns in human bronchial epithelial (HBE) cells
Primary HBE cells from 5 donors were transfected with siRNA that targets IRF7 or a non-silencing siRNA control (all-star control). The HBE cells were allowed to recover for 24 h, and then they were exposed to HRV or medium control for another 24 h. As illustrated in Fig. 1, HRV induced the expression of IRF7 protein in HBE cells that were transfected with control reagents (all-star control siRNA, lipid reagent alone, medium only), but not in HBE cells that were transfected with siRNA reagents that target IRF7 (siRNA #1, siRNA #2). Of note, IRF7 protein was not induced by the all-star control siRNA in HBE cells in the absence of HRV exposure.
Next, global patterns of gene expression were profiled on microarrays. Initially, we analyzed the data pertaining to the response to HRV in HBE cells transfected with the all-star control. As illustrated in Fig. 2, 437 genes were upregulated in the response and 139 genes were downregulated (False Discovery Rate (FDR) < 0.05, Additional file 1). The biological function of the differentially expressed genes was interrogated using InnateDb. The data showed the upregulated genes were mainly involved in type I and II interferon signaling, RIG-I/MDA5 signaling, antigen processing and presentation, and apoptosis (Additional file 2). The downregulated genes were not strongly enriched for immune-related pathways, therefore we utilized the Gene Ontology (GO) database to interrogate their function. This analysis showed that the downregulated genes were enriched with the GO biological process term epithelial development (GO0060429; ALDH3A2, CTGF, DKK1, EMP1, ERRFI1, INSR, JAG1, KRT15, PTHLH, SCEL, SOX21; adjusted p-value = 0.0079).
IRF7 is a putative molecular driver of interferon-induced gene networks
Upstream regulator analysis was employed to infer the molecular drivers that give rise to HRV-induced gene expression patterns in HBE cells (transfected with the all-star control) . In this analysis, the differentially expressed genes from Fig. 2 were tested for enrichment of downstream target genes of known transcriptional regulators and other signal transduction pathways. As illustrated in Table 1, the most significant candidate drivers were IFNG, IFNA2, IRF7, TNF, and IFNL1. Notably, IRF7 was the only driver amongst these candidates that was also upregulated in the response (Additional file 1). To illustrate the role of IRF7 in the regulation of HRV responses, we utilized experimentally supported findings from the Ingenuity Systems database to reconstruct the IRF7 network . The resulting network is shown in Fig. 3.
Knockdown of IRF7 perturbs HBE responses to rhinovirus
We next examined the impact of IRF7 knockdown on the HRV response. Gene expression patterns were compared between HRV-induced HBE cells transfected with IRF7-siRNA #1 versus the all-star control. The data showed that knockdown of IRF7 resulted in decreased expression of 447 genes, and increased expression of 181 genes (Additional file 3, Fig. 4). The downregulated genes were mainly involved in type I and II interferon signaling, Toll-like receptor signaling, NOD-like receptor signaling, and RIG-I/MDA5 signaling (Additional file 4). As a supplementary analysis, we also examined the function of the downregulated genes after removal of known IRF7 target genes. These residual genes were associated with NOD-like receptor signaling pathways, trafficking of connexons, and metabolic pathways (Additional file 5). The upregulated genes were involved in the response to oxidative stress (GO:0006979; ADA, ARL6IP5, CYGB, FOS, MMP14, NET1, NQO1, PLK3, SNCA, TRPC6, UCP2; p-value = 0.0004) and inflammation (GO:0050729; ADORA2B, C3, IL1RL1, IL33, SERPINE1, TNFSF11; p-value = 0.0005), although it is noteworthy that these findings were not significant after adjustment of the p-values for multiple testing.
We next wanted to investigate the part of the HRV response that was insensitive to IRF7 knockdown. In HBE cells that were transfected with the all-star control siRNA, 576 genes were differentially expressed after HRV exposure, but only 177 of these genes (31 %) were sensitive to IRF7 knockdown (Fig. 5). Moreover, the IRF7 gene network depicted in Fig. 3 contained 86 genes, and only 41 of these (47 %) were sensitive to IRF7 knockdown (Fig. 4, Additional file 6). We analyzed these HRV responsive genes that were insensitive to IRF7 knockdown using upstream regulator analysis. The data showed the highest ranking candidate transcription factors driving this response were NF-κB, STAT1, STAT3, and IRF1 (Additional file 7).
Validation of the microarray findings
To confirm the findings from the microarray experiment, a subset of genes were selected for validation at the mRNA and/or protein levels. As illustrated in Fig. 6, Reverse Transcription quantitative PCR (RT-qPCR) analysis demonstrated that knockdown of IRF7 in HBE cells followed by exposure to HRV resulted in reduced expression of CCL5, CXCL10, HERC5, IL-6, RSAD2 (viperin), and increased expression of CXCL5, and this was consistent with the microarray data (Additional file 3). Expression levels of CCL5, CXCL5, CXCL10, and viperin were also examined at the protein level (Fig. 7). These data confirmed that the mRNA expression changes observed upon siRNA knockdown were associated with parallel changes at the protein level, although these changes were more pronounced with siRNA #2 than with siRNA #1.
IRF7 is a master regulator of type I interferon gene expression and antiviral immunity, but the role of IRF7 in HBE responses to HRV has not been previously investigated. Employing microarray profiling we showed that hundreds of genes were upregulated in HBE upon HRV infection. Upstream regulator analysis  suggested that IRF7 was a major molecular driver of this response. We then showed that knockdown of IRF7 in HBE altered HRV responses, by decreasing the upregulation of genes involved in type I and II interferon signaling and innate immunity (toll-like receptor signaling, NOD-like receptor signaling, RIG-I/MDA5 signaling), and increasing the expression of genes that promote inflammation (C3, CXCL5, IL-33, IL1RL1) and the response to oxidative stress. We validated a subset of the differentially regulated genes (CCL5, CXCL5, CXCL10, HERC5, IL-6, RSAD2/viperin) at the mRNA and/or protein levels, thus confirming the findings from the microarray experiment. The data also showed that only around 30 % of the HRV responsive genes and half of known IRF7 target genes were sensitive to IRF7 knockdown, and these genes exhibited modest fold changes in our experiments (less than 2 fold), suggesting that other transcription factors are likely to be important in the regulation of interferon gene expression programs in HBE cells. Upstream regulator analysis of the part of the HRV response that was insensitive to IRF7 knockdown suggested roles for NF-κB, STAT1, STAT3, and IRF1 in the regulation of the response. Overall, our findings demonstrate that IRF7 regulates the expression of genes involved in antiviral immunity, inflammation, and the response to oxidative stress in HBE cells during HRV infections, and also suggest that additional transcription factors play a major role in the regulation of this response.
We showed that the upregulation of CCL5, CXCL10, HERC5, IL-6, and RSAD2 by HRV was reduced by IRF7 knockdown in HBE cells. Regulation of these genes by IRF7 has been previously reported in other contexts [26–29]. CCL5 (RANTES) is a potent chemoattractant for T cells, monocytes, and eosinophils. CCL5 was upregulated in nasal wash samples from wheezing, HRV-infected infants in comparison to their nonwheezing counterparts . CXCL10 is chemotactic for T cells and natural killer cells, and was strongly upregulated in airway epithelial cells during HRV infections . CXCL10 promotes airways inflammation and airways hyperresponsiveness in mouse models of asthma . HERC5 is an E3 protein ligase that mediates ISGylation of protein targets. It catalyzes the conjugation of ISG15 onto IRF3, which inhibits the ubiquitination and degradation of IRF3 by preventing IRF3 from interacting with Pin1. This results in the sustained activation of IRF3 and the antiviral response . IL-6 promotes the acute phase response, stimulates T cell responses and antibody production, and it is increased during symptomatic rhinovirus infections . IRF7 was required for IL-6 production by poly (I-C) stimulated THP-1 monocytes . In contrast, IRF7 was dispensable for IL-6 production in plasmacytoid and myeloid dendritic cells stimulated by TLR7 or TLR9 ligands . Viperin (RSAD2) is an antiviral protein that can inhibit influenza budding from the plasma membrane by interfering with lipid rafts . Knockdown of viperin in airway epithelial cells increased HRV replication . It is not possible to determine from our experiments if IRF7 regulates HRV-induced airway inflammatory responses in vivo, however, Girkin et al. showed that intranasal delivery of siRNA that targets IRF7 reduced the infiltration of neutrophils and macrophages in the bronchoalveolar lavage fluid of non-allergic mice after HRV-1B infection . Taken together, these studies suggest that IRF7 regulates HRV-induced gene network patterns and ensuing cellular immune responses in the airways.
In contrast to genes involved in antiviral defense, knockdown of IRF7 resulted in increased expression of CXCL5, IL-33, and IL1RL1. CXCL5 (ENA-78) is important for neutrophil homeostasis and recruitment to the lung during infections . IL-33 is a nuclear-associated alarmin expressed in the airway epithelium. It signals through the ST2 receptor (IL1RL1), which is highly expressed on T-helper type 2 (Th2) lymphocytes and group 2 innate lymphoid cells (ILC2s) . IL-33 is upregulated in the airways during HRV-induced asthma exacerbations . Moreover, supernatants from HRV-infected epithelial cells stimulate the production of Th2 cytokines by Th2 and ILC2 cells in an IL-33 dependant manner . IRF7 binds to the IL-33 promoter and regulates IL-33 transcription in human monocyte THP-1 cell lines . Given that IL-33 is a potent driver of asthmatic ILC2 and Th2 responses, and we found that knockdown of IRF7 results in upregulation of IL-33, this suggests the hypothesis that promoting IRF7 responses during asthma exacerbations may have a dual beneficial effect of increasing antiviral immunity and supressing IL-33 responses.
Knockdown of IRF7 decreased expression of genes involved in Nod-like receptor signaling, and increased expression of genes involved in the response to oxidative stress (e.g. cytoglobin (CYGB), NAD(P)H:quinone acceptor oxidoreductase 1 (NQO1), alpha-synuclein (SNCA), short transient receptor potential channel 6 (TRPC6)) [41–44]. In this context it is noteworthy that the generation of reactive oxygen species during HRV infection compromises epithelial barrier function, and this can lead to secondary bacterial infections [45, 46]. HRV infection stimulates the epithelial production of reactive oxygen species by triggering the dsRNA receptor Nod-like receptor X-1 (NLRX-1) . Cross-talk between oxidative stress pathways and IRF7 signaling was reported in the context of respiratory syncytial virus infection . The role of IRF7 in the regulation of Nod-like receptor signaling and oxidative stress pathways during HRV infections merits further investigation in follow-up studies.
Upstream regulator analysis demonstrated that IRF7 was the most significant candidate molecular driver that was upregulated in the HRV response. However, our siRNA studies showed that only about 30 % of the HRV responsive genes and half of known IRF7 target genes were sensitive to IRF7 knockdown. Whilst variable and incomplete knockdown and the potential for non-specific activation of antiviral pathways is a well-known limitation of siRNA technology, an alternative explanation of the data is that other transcription factors may be playing a role. Employing upstream regulator analysis we showed that NF-κB, STAT1, STAT3, and IRF1 were candidate drivers of the part of the HRV response that was insensitive to IRF7 knockdown. In this context it is noteworthy we have reported that HRV-induced upregulation of CXCL10 in bronchial epithelial cells was dependent on IRF1 . Wang et al. have also reported that knockdown of IRF3 abrogated HRV-induced type I and III interferon responses in BEAS-2B epithelial cells . Taken together, these observations suggest that interferon responses in HRV-induced HBE cells can be variably mediated by interactions between IRF1, IRF3, and IRF7.
This study has limitations that should be acknowledged. First, the gene expression studies were conducted at a single time point post infection (24 h), so it was not possible to determine the role of IRF7 in the kinetics of the response. Second, siRNA-mediated gene silencing was used to investigate the role of IRF7 in the regulation of gene expression patterns, and from this data it was not possible to determine if IRF7 was playing a direct or indirect regulatory role. Third, prior knowledge was used to infer the molecular drivers of the response, and this analysis cannot reveal novel molecular drivers. This latter limitation can be addressed using unbiased data-driven algorithms to infer the molecular drivers of the response . Although we acknowledge that the use of epithelial cells grown in submersion culture, as opposed to differentiated air-liquid interface cultures, may be considered a limitation of the current study, we believe that the data obtained are still highly relevant. This is based on the observation that the profile of HRV-upregulated genes in these cultures, in both the present and in an earlier study , aligns closely with those observed in nasal epithelial scrapings from individuals experimentally infected with HRV in vivo .
Our study demonstrates that IRF7 regulates the expression of antiviral, inflammatory and oxidative stress response pathways during HRV infections in HBE cells, and suggests that additional transcription factors play an important role in the response.
Availability of supporting data
The microarray data are available at the Gene Expression Omnibus repository (accession GSE70190).
adenosine A2b receptor
aldehyde dehydrogenase 3 family member A2
analysis of variance
ADP-ribosylation factor-like 6 interacting protein 5
bronchial epithelial cell growth medium
complement component 3
chemokine (C-C motif) ligand 5, also known as regulated on activation, normal T cell expressed and secreted (RANTES)
Cadherin-related family member 3
connective tissue growth factor
chemokine (C-X-C motif) ligand 10
chemokine (C-X-C motif) ligand 5
dickkopf WNT signaling pathway inhibitor 1
enzyme-linked immunosorbent assay
epithelial membrane protein 1
ERBB receptor feedback inhibitor 1
fetal bovine serum
False Discovery Rate
FBJ murine osteosarcoma viral oncogene homolog
human bronchial epithelial cells
Hank’s balanced salt solution
HECT and RLD domain containing E3 ubiquitin protein ligase 5
human rhinovirus type 16
human rhinovirus type C
interferon, lambda 1
interleukin 1 receptor-like 1
group 2 innate lymphoid cells
interferon regulatory factor 7
interferon stimulated gene that encodes interferon-induced 17 kDa protein
jagged 1, KRT15, keratin 15 type I
linear models for microarray data
melanoma differentiation associated gene 5
matrix metallopeptidase 14
MX dynamin-like GTPase 1
neuroepithelial cell transforming 1
nuclear factor kappa-light-chain-enhancer of activated B cells
receptor nucleotide-binding oligomerization domain receptor
NAD(P)H dehydrogenase, quinone 1
2′-5′-oligoadenylate synthetase 1
protein kinase RNA-activated
polo-like kinase 3
pathogen recognition receptor
parathyroid hormone-like hormone
regulated on activation, normal T cell expressed and secreted
retinoic acid-inducible gene 1
reverse transcription quantitative polymerase chain reaction
sodium dodecyl sulfate polyacrylamide gel electrophoresis
serpin peptidase inhibitor clade E
small interfering RNA
sex determining region Y box 21
T-helper type 2
tumor necrosis factor
tumor necrosis factor (ligand) superfamily member 11, also known as receptor activator of nuclear factor-κB ligand (RANKL)
transient receptor potential cation channel, subfamily C, member 6
uncoupling protein 2
Leigh R, Proud D. Virus-induced modulation of lower airway diseases: pathogenesis and pharmacologic approaches to treatment. Pharmacol Ther. 2015;148:185–98.
Gern JE. The ABCs of rhinoviruses, wheezing, and asthma. J Virol. 2010;84(15):7418–26.
Bochkov YA, Watters K, Ashraf S, Griggs TF, Devries MK, Jackson DJ, et al. Cadherin-related family member 3, a childhood asthma susceptibility gene product, mediates rhinovirus C binding and replication. Proc Natl Acad Sci U S A. 2015.
Wang Q, Miller DJ, Bowman ER, Nagarkar DR, Schneider D, Zhao Y, et al. MDA5 and TLR3 initiate pro-inflammatory signaling pathways leading to rhinovirus-induced airways inflammation and hyperresponsiveness. PLoS Pathog. 2011;7(5), e1002070.
Wang Q, Nagarkar DR, Bowman ER, Schneider D, Gosangi B, Lei J, et al. Role of double-stranded RNA pattern recognition receptors in rhinovirus-induced airway epithelial cell responses. J Immunol. 2009;183(11):6989–97.
Slater L, Bartlett NW, Haas JJ, Zhu J, Message SD, Walton RP, et al. Co-ordinated role of TLR3, RIG-I and MDA5 in the innate response to rhinovirus in bronchial epithelium. PLoS Pathog. 2010;6(11), e1001178.
Ning S, Pagano JS, Barber GN. IRF7: activation, regulation, modification and function. Genes Immun. 2011;12(6):399–414.
Marie I, Durbin JE, Levy DE. Differential viral induction of distinct interferon-alpha genes by positive feedback through interferon regulatory factor-7. EMBO J. 1998;17(22):6660–9.
Schoggins JW, Rice CM. Interferon-stimulated genes and their antiviral effector functions. Curr Opin Virol. 2011;1(6):519–25.
Honda K, Yanai H, Negishi H, Asagiri M, Sato M, Mizutani T, et al. IRF-7 is the master regulator of type-I interferon-dependent immune responses. Nature. 2005;434(7034):772–7.
Ciancanelli MJ, Huang SX, Luthra P, Garner H, Itan Y, Volpi S, et al. Infectious disease. Life-threatening influenza and impaired interferon amplification in human IRF7 deficiency. Science. 2015;348(6233):448–53.
Daffis S, Suthar MS, Szretter KJ, Gale Jr M, Diamond MS. Induction of IFN-beta and the innate antiviral response in myeloid cells occurs through an IPS-1-dependent signal that does not require IRF-3 and IRF-7. PLoS Pathog. 2009;5(10), e1000607.
Proud D, Turner RB, Winther B, Wiehler S, Tiesman JP, Reichling TD, et al. Gene expression profiles during in vivo human rhinovirus infection: insights into the host response. Am J Respir Crit Care Med. 2008;178(9):962–8.
Bosco A, Ehteshami S, Panyala S, Martinez FD. Interferon regulatory factor 7 is a major hub connecting interferon-mediated responses in virus-induced asthma exacerbations in vivo. J Allergy Clin Immunol. 2012;129(1):88–94.
Lewis TC, Henderson TA, Carpenter AR, Ramirez IA, McHenry CL, Goldsmith AM, et al. Nasal cytokine responses to natural colds in asthmatic children. Clin Exp Allergy. 2012;42(12):1734–44.
Churchill L, Chilton FH, Resau JH, Bascom R, Hubbard WC, Proud D. Cyclooxygenase metabolism of endogenous arachidonic acid by cultured human tracheal epithelial cells. Am Rev Respir Dis. 1989;140:449–59.
Sanders SP, Siekierski ES, Porter JD, Richards SM, Proud D. Nitric oxide inhibits rhinovirus-induced cytokine production and viral replication in a human respiratory epithelial cell line. J Virol. 1998;72:934–42.
Chomczynski P, Sacchi N. Single-step method of RNA isolation by acid guanidinium thiocyanate-phenol-chloroform extraction. Anal Biochem. 1987;162:156–9.
Sanders SP, Siekierski ES, Richards SM, Porter JD, Imani F, Proud D. Rhinovirus infection induces expression of type 2 nitric oxide synthase in human respiratory epithelial cells in vitro and in vivo. J Allergy Clin Immunol. 2001;107:235–43.
Kauffmann A, Gentleman R, Huber W. arrayQualityMetrics--a bioconductor package for quality assessment of microarray data. Bioinformatics. 2009;25(3):415–6.
Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, et al. Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003;4(2):249–64.
Dai M, Wang P, Boyd AD, Kostov G, Athey B, Jones EG, et al. Evolving gene/transcript definitions significantly alter the interpretation of GeneChip data. Nucleic Acids Res. 2005;33(20), e175.
Lu J, Kerns RT, Peddada SD, Bushel PR. Principal component analysis-based filtering improves detection for Affymetrix gene expression arrays. Nucleic Acids Res. 2011;39(13), e86.
Smyth GK. Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004;3:Article3.
Kramer A, Green J, Pollard Jr J, Tugendreich S. Causal analysis approaches in Ingenuity Pathway Analysis. Bioinformatics. 2014;30(4):523–30.
Lin R, Genin P, Mamane Y, Hiscott J. Selective DNA binding and association with the CREB binding protein coactivator contribute to differential activation of alpha/beta interferon genes by interferon regulatory factors 3 and 7. Mol Cell Biol. 2000;20(17):6342–53.
Goubau D, Romieu-Mourez R, Solis M, Hernandez E, Mesplede T, Lin R, et al. Transcriptional re-programming of primary macrophages reveals distinct apoptotic and anti-tumoral functions of IRF-3 and IRF-7. Eur J Immunol. 2009;39(2):527–40.
Barnes BJ, Richards J, Mancl M, Hanash S, Beretta L, Pitha PM. Global and distinct targets of IRF-5 and IRF-7 during innate response to viral infection. J Biol Chem. 2004;279(43):45194–207.
Sweeney SE. Targeting interferon regulatory factors to inhibit activation of the type I IFN response: implications for treatment of autoimmune disorders. Cell Immunol. 2011;271(2):342–9.
Pacifico L, Iacobini M, Viola F, Werner B, Mancuso G, Chiesa C. Chemokine concentrations in nasal washings of infants with rhinovirus illnesses. Clin Infect Dis. 2000;31(3):834–8.
Spurrell JC, Wiehler S, Zaheer RS, Sanders SP, Proud D. Human airway epithelial cells produce IP-10 (CXCL10) in vitro and in vivo upon rhinovirus infection. Am J Physiol Lung Cell Mol Physiol. 2005;289(1):L85–95.
Medoff BD, Sauty A, Tager AM, Maclean JA, Smith RN, Mathew A, et al. IFN-gamma-inducible protein 10 (CXCL10) contributes to airway hyperreactivity and airway inflammation in a mouse model of asthma. J Immunol. 2002;168(10):5278–86.
Shi HX, Yang K, Liu X, Liu XY, Wei B, Shan YF, et al. Positive regulation of interferon regulatory factor 3 activation by Herc5 via ISG15 modification. Mol Cell Biol. 2010;30(10):2424–36.
Zhu Z, Tang W, Ray A, Wu Y, Einarsson O, Landry ML, et al. Rhinovirus stimulation of interleukin-6 in vivo and in vitro. Evidence for nuclear factor kappa B-dependent transcriptional activation. J Clin Invest. 1996;97(2):421–30.
Wang X, Hinson ER, Cresswell P. The interferon-inducible protein viperin inhibits influenza virus release by perturbing lipid rafts. Cell Host Microbe. 2007;2(2):96–105.
Girkin J, Hatchwell L, Foster P, Johnston SL, Bartlett N, Collison A, et al. CCL7 and IRF-7 mediate hallmark inflammatory and IFN responses following rhinovirus 1B infection. J Immunol. 2015;194(10):4924–30.
Mei J, Liu Y, Dai N, Hoffmann C, Hudock KM, Zhang P, et al. Cxcr2 and Cxcl5 regulate the IL-17/G-CSF axis and neutrophil homeostasis in mice. J Clin Invest. 2012;122(3):974–86.
Molofsky AB, Savage AK, Locksley RM. Interleukin-33 in tissue homeostasis, injury, and inflammation. Immunity. 2015;42(6):1005–19.
Jackson DJ, Makrinioti H, Rana BM, Shamji BW, Trujillo-Torralbo MB, Footitt J, et al. IL-33-dependent type 2 inflammation during rhinovirus-induced asthma exacerbations in vivo. Am J Respir Crit Care Med. 2014;190(12):1373–82.
Sun L, Zhu Z, Cheng N, Yan Q, Ye RD. Serum amyloid A induces interleukin-33 expression through an IRF7-dependent pathway. Eur J Immunol. 2014;44(7):2153–64.
Dinkova-Kostova AT, Talalay P. NAD(P)H:quinone acceptor oxidoreductase 1 (NQO1), a multifunctional antioxidant enzyme and exceptionally versatile cytoprotector. Arch Biochem Biophys. 2010;501(1):116–23.
le Thuy TT, Matsumoto Y, Thuy TT, Hai H, Suoh M, Urahara Y, et al. Cytoglobin deficiency promotes liver cancer development from hepatosteatosis through activation of the oxidative stress pathway. Am J Pathol. 2015;185(4):1045–60.
Hashimoto M, Hsu LJ, Rockenstein E, Takenouchi T, Mallory M, Masliah E. alpha-Synuclein protects against oxidative stress via inactivation of the c-Jun N-terminal kinase stress-signaling pathway in neuronal cells. J Biol Chem. 2002;277(13):11465–72.
Ding Y, Winters A, Ding M, Graham S, Akopova I, Muallem S, et al. Reactive oxygen species-mediated TRPC6 protein activation in vascular myocytes, a mechanism for vasoconstrictor-regulated vascular tone. J Biol Chem. 2011;286(36):31799–809.
Comstock AT, Ganesan S, Chattoraj A, Faris AN, Margolis BL, Hershenson MB, et al. Rhinovirus-induced barrier dysfunction in polarized airway epithelial cells is mediated by NADPH oxidase 1. J Virol. 2011;85(13):6795–808.
Sajjan U, Wang Q, Zhao Y, Gruenert DC, Hershenson MB. Rhinovirus disrupts the barrier function of polarized airway epithelial cells. Am J Respir Crit Care Med. 2008;178(12):1271–81.
Unger BL, Ganesan S, Comstock AT, Faris AN, Hershenson MB, Sajjan US. Nod-like receptor X-1 is required for rhinovirus-induced barrier dysfunction in airway epithelial cells. J Virol. 2014;88(7):3705–18.
Fang L, Choudhary S, Tian B, Boldogh I, Yang C, Ivanciuc T, et al. Ataxia telangiectasia mutated kinase mediates NF-kappaB serine 276 phosphorylation and interferon expression via the IRF7-RIG-I amplification loop in paramyxovirus infection. J Virol. 2015;89(5):2628–42.
Zaheer RS, Proud D. Human rhinovirus-induced epithelial production of CXCL10 is dependent upon IFN regulatory factor-1. Am J Respir Cell Mol Biol. 2010;43(4):413–21.
Huynh-Thu VA, Irrthum A, Wehenkel L, Geurts P. Inferring regulatory networks from expression data using tree-based methods. PLoS One. 2010;5(9):e12776.
Proud D, Hudy MH, Wiehler S, Zaheer RS, Amin MA, Pelikan JB, et al. Cigarette smoke modulates expression of human rhinovirus-induced airway epithelial host defense genes. PLoS One. 2012;7, e40762.
This work was supported by the Department of Health, Government of Western Australia, and grant number 43923 from the Canadian Institutes of Health Research. AB is supported by a BrightSpark Foundation McCusker Fellowship. D. Proud holds a Canada Research Chair in Inflammatory Airway Disease.
The authors declare that they have no competing interests.
AB participated in the design and conception of the study, analyzed the gene expression profiling experiments and drafted the manuscript. SW carried out experiments, analyzed data, and helped to draft the manuscript. DP conceived and designed the study, analyzed data, and helped to draft the manuscript. All authors read and approved the final manuscript.
This excel file contains the list of differentially expressed genes that were induced by HRV in HBE cells transfected with the all-star control siRNA. (XLSX 71 kb)
This excel file contains the list of biological pathways that were upregulated by HRV in HBE cells transfected with the all-star control siRNA. (XLSX 23 kb)
This excel file contains the list of differentially expressed genes identified in HRV-induced HBE cells transfected with IRF7-siRNA #1 versus the all-star control siRNA. (XLSX 77 kb)
This excel file contains the list of biological pathways that were decreased in HRV-induced HBE cells transfected with IRF7-siRNA#1 versus the all-star control siRNA. (XLSX 13 kb)
This excel file contains the results of a pathways analysis of remaining genes (after removal of known IRF7 target genes) that were decreased in HRV-induced HBE cells transfected with IRF7-siRNA#1 versus the all-star control siRNA. (XLSX 11 kb)
This excel file contains the results of a differential gene expression analysis of IRF7 gene networks in HBE cells transfected with IRF7-siRNA#1 versus the all-star control siRNA. (XLSX 18 kb)
This excel file contains the results from an Upstream regulator analysis of the HRV responsive genes in HBE cells that were insensitive to IRF7 knockdown. (XLSX 11 kb)
About this article
Cite this article
Bosco, A., Wiehler, S. & Proud, D. Interferon regulatory factor 7 regulates airway epithelial cell responses to human rhinovirus infection. BMC Genomics 17, 76 (2016). https://doi.org/10.1186/s12864-016-2405-z