Open Access
Open Peer Review

This article has Open Peer Review reports available.

How does Open Peer Review work?

Whole blood gene expression in infants with respiratory syncytial virus bronchiolitis

  • Hans-Olav Fjaerli1Email author,
  • Geir Bukholm2,
  • Anne Krog2,
  • Camilla Skjaeret2,
  • Marit Holden3 and
  • Britt Nakstad1
BMC Infectious Diseases20066:175

DOI: 10.1186/1471-2334-6-175

Received: 12 May 2006

Accepted: 13 December 2006

Published: 13 December 2006



Respiratory syncytial virus (RSV) is a major cause of viral bronchiolitis in infants worldwide, and environmental, viral and host factors are all of importance for disease susceptibility and severity. To study the systemic host response to this disease we used the microarray technology to measure mRNA gene expression levels in whole blood of five male infants hospitalised with acute RSV, subtype B, bronchiolitis versus five one year old male controls exposed to RSV during infancy without bronchiolitis. The gene expression levels were further evaluated in a new experiment using quantitative real-time polymerase chain reaction (QRT-PCR) both in the five infants selected for microarray and in 13 other infants hospitalised with the same disease.


Among the 30 genes most differentially expressed by microarray nearly 50% were involved in immunological processes. We found the highly upregulated interferon, alpha-inducible protein 27 (IFI27) and the highly downregulated gene Charcot-Leyden crystal protein (CLC) to be the two most differentially expressed genes in the microarray study. When performing QRT-PCR on these genes IFI27 was upregulated in all but one infant, and CLC was downregulated in all 18 infants, and similar to that given by microarray.


The gene IFI27 is upregulated and the gene CLC is downregulated in whole blood of infants hospitalised with RSV, subtype B, bronchiolitis and is not reported before. More studies are needed to elucidate the specificity of these gene expressions in association with host response to this virus in bronchiolitis of moderate severity.


Why some infants develop bronchiolitis when exposed to respiratory syncytial virus (RSV) is poorly understood, and several aspects of environmental and host immunity have been extensively studied [1, 2]. The role of the virus itself in modulating host immune response has also been investigated. Recent studies have revealed an important role of certain RSV surface proteins to counteract the host interferon (IFN) response, a response known to be important for viral clearance [3, 4].

Microarray can detect the simultaneous expressions and interactions of thousands of genes [5]. This technology was used as a hypothesis generating tool to identify the most differentially expressed genes in whole blood of five male infants hospitalised with RSV, subtype B, bronchiolitis versus five one year old male controls exposed to RSV during infancy without being hospitalized and/or treated for acute bronchiolitis. To further validate our results quantitative real-time polymerase chain reaction (QRT-PCR) was performed on the most differentially expressed genes given by the microarray experiment both in the infants selected for microarray and in a new population of 13 more infants admitted to our hospital with the same disease [6]. Four of the five one year old male controls were used as a pooled exogenous control in the QRT-PCR study.

The infants were diagnosed with bronchiolitis if having symptoms from the lower airways characterized by wheezing, dyspnea, respiratory distress, poor feeding, tachypnea and fine crepitations upon auscultation when examined by the doctor on call in the emergency room [7]. The doctor also conducted a structured clinical interview with the parents and NPA was taken for viral analysis with multiplex RT-PCR [8]. Median duration of symptoms prior to hospitalisation was 4 (range 2–7) days, median age at admission was 3 (range 1–8) months and median duration of hospitalisation was 3 (range 1–8) days. Severity of illness was comparable in all 18 infants, except for one infant who had a prolonged hospitalisation period of 8 days. However, all infants received the same routine treatment including inhalations with nebulised racemic epinephrine and none of the infants needed artificial ventilation. More information on the selection of cases and controls is described elsewhere [See Additional file 1].

Microarray results

The Feature Extraction software identified 19663 unique gene sequences. After normalization and filtering of control spots and spots with signal intensities below 300 in either channel a total of 15761 genes were available for further evaluation. When performing statistical analysis of the data using the BAMarray software we found 439 genes to be significantly differentially expressed between cases and controls. Furthermore, by also analyzing the same data with the Rosetta Luminator software, a list of 1544 genes with P value ≤ 0.01 was identified. 54 of these genes had mean fold changes of at least ± 2.50. Finally, by comparing the list of 439 genes found by BAMarray with the list of 54 genes identified with the Rosetta Luminator software, a total of 30 genes occurred in both lists and were selected for final presentation (tab 1). More methodological details are described elsewhere [See Additional file 1].
Table 1

Significant genes with mean fold changes1 of at least ± 2.50 as given by microarray2

GeneBank identifier

Gene symbol

Gene name

Biological process

Fold change



Interferon, alpha-inducible protein 27

Immune response





Defence response




Interferon-induced protein 44

Immune response




Membrane-spanning 4-domains, subfamily A, member 4

Signal transduction




Matrix metalloproteinase 9

Collagen catabolism




Hemoglobin, zeta

Oxygen transport




DNA polymerase-transactivated protein 6





Hemoglobin, delta

Oxygen transport




Hemoglobin, epsilon 1

Oxygen transport




2'-5'-oligoadenylate synthetase 3, 100kDa

Immune response




Fc fragment of IgG, high affinity Ia, receptor (CD64)

Immune response




Syntaxin binding protein 2

Intracellular protein transport




Epithelial stromal interaction 1 (breast)





Eukaryotic translation initiation factor 2-alpha kinase 2

Immune response




Erythroid associated factor

Hemoglobin metabolism




Macrophage receptor with collagenous structure

Response to pathogenic bacteria




2,3-bisphosphoglycerate mutase

Carbohydrate metabolism




NAD(P)H dehydrogenase, quinone 2

Electron transport




Interferon-induced protein 44-like

Immune response




Interferon, alpha-inducible protein (clone IFI-15K)

Immune response




UDP-N-acetyl-alpha-D-galactosamine:polypeptide N-acetylgalactosaminyltransferase 14






Immune response




Lutheran blood group (Auberger b antigen included)

Cell adhesion




G protein-coupled receptor 78

G-protein coupled receptor protein signaling pathway




Dopamine receptor D4

G-protein coupled receptor protein signaling pathway




Proline-rich protein BstNI subfamily 3

G-protein coupled receptor protein signaling pathway




Homo sapiens major histocompatibility complex, class II, DR beta 4

Immune response




Hemoglobin alpha 1

Oxygen transport




Major histocompatibility complex, class II, DQ alpha 1

Immune response




Charcot-Leyden crystal protein

Antimicrobial humoral response


1Gene expressions given as mean of the individual fold changes (signal intensity in a case minus signal intensity in the corresponding control). 2Gene expressions as given by microarray in whole blood of five male infants hospitalised with respiratory syncytial virus, subtype B, bronchiolitis versus five one year old male controls exposed to RSV during infancy but not hospitalised and/or treated for bronchiolitis

22 of these 30 genes were upregulated while 8 genes were downregulated, and 13 of the 30 genes were involved in several aspects of immune response. Among these were the upregulated interferon inducible genes IFI27, IFI44, EIF2AK2, IFI44L, OAS3 and G1P2 and the downregulated genes HLA-DRB4, HLA-DQA1 and CLC (tab 1).

QRT-PCR results

A QRT-PCR study with TaqMan Low Density Array (TLDA) cards was performed on 23 of the most differentially expressed genes given by microarray to compare the two methods before the study was applied on a new patient population. Six genes failed amplification and were excluded from further analyses. More methodological details are described elsewhere [See Additional file 1].

Five of the 17 genes available for further evaluation (DNAPTP6, HP, IFI27, MS4A4A and STXBP2) were upregulated and one gene (CLC) was downregulated in all five male infants when compared to the pooled exogenous control. Seven genes (BPGM, G1P2, HBD, IFI44L, MARCO, MMP9 and NQO2) were upregulated in all but one infant while four genes (EPSTI1, ERAF, HBE1 and IFI44) were upregulated in all but two infants (fig 1). When comparing the gene expression results of the TLDA study and the microarray experiment for each of the five infants the expression profiles were quite similar for the genes CLC, DNAPTP6, IFI27 and MS4A4A. The remaining genes showed more variations (fig 1).
Figure 1

Evaluation of genes differentially expressed by microarray (blue bars) with QRT-PCR (red bars). The microarray experiment was conducted as a case-control study with whole blood gene expressions given as ratio of signal intensity in each of five male infants hospitalised with respiratory syncytial virus, subtype B, bronchiolitis versus signal intensity in a corresponding one year old male control exposed to the virus during infancy but not hospitalised and/or treated for bronchiolitis. The quantitative real-time polymerase chain reaction (QRT-PCR) study used TaqMan Low Density Array (TLDA) cards with gene expression given as mean relative quantification (RQ) from triplets of each gene and analyzed in the same five male infants as selected for microarray. A pooled sample from four of the five one year old male controls was used as exogenous control and beta-glucuronidase (GUSB) was used as endogenous control in the QRT-PCR study. The following genes (gene symbol and name) were evaluated; BPGM: 2,3-bisphosphoglycerate mutase, CLC: Charcot-Leyden crystal protein, DNAPTP6: DNA polymerase-transactivated protein 6, EPSTI1: Epithelial stromal interaction 1 (breast), ERAF: Erythroid associated factor, G1P2: Interferon, alpha-inducible protein (clone IFI-15K), HBD: Hemoglobin, delta, HBE1: Hemoglobin, epsilon 1, HP: Haptoglobin, IFI27: Interferon, alpha-inducible protein 27, IFI44: Interferon-induced protein 44, IFI44L: Interferon-induced protein 44-like, MARCO: Macrophage receptor with collagenous structure, MMP9: Matrix metalloproteinase 9, MS4A4A: Membrane-spanning 4-domains, subfamily A, member 4, NQO2: NAD(P)H dehydrogenase, quinone 2, STXBP2: Syntaxin binding protein 2

To evaluate the impact of the 17 genes expressed by microarray a new cohort was tested using QRT-PCR. In the new cohort of 13 infants hospitalised with the same disease 11 of the 17 genes were significantly differentially expressed in accordance with the results from the microarray study (p < 0.01; independent samples t-test). In order to see if some of the genes associated with immunological processes could also be differentially expressed in the new cohort, as well as in the five infants selected for microarray, we evaluated the expressions of the six immune response genes available for studies with TLDA cards. The gene IFI27 was upregulated in 17 infants and slightly downregulated in one infant while the gene IFI44 was downregulated in 15 infants and upregulated in three infants. The genes IFI44L and MARCO were downregulated in 14 infants and upregulated in four infants, and the gene G1P2 was downregulated and upregulated in nine infants, respectively. Finally, the gene CLC was downregulated in all 18 infants (fig 2).
Figure 2

QRT-PCR study of six whole blood immune response genes as given by microarray. A quantitative real-time polymerase chain reaction (QRT-PCR) study using TaqMan Low Density Array (TLDA) cards was performed with gene expressions given as mean relative quantification (RQ) from triplets of each gene. Analyzed in 18 infants hospitalised with respiratory syncytial virus, subtype B, bronchiolitis versus a pooled sample from four one year old male children exposed to RSV during infancy but not treated and/or hospitalised for bronchiolitis during infancy as exogenous control and with beta-glucuronidase (GUSB) as endogenous control. The following genes (gene symbol and name) were studied; CLC: Charcot-Leyden crystal protein, G1P2: Interferon, alpha-inducible protein (clone IFI-15K), IFI27: Interferon, alpha-inducible protein 27, IFI44: Interferon-induced protein 44, IFI44L: Interferon-induced protein 44-like, MARCO:Macrophage receptor with collagenous structure


The aim of the present study was to use microarray as a hypothesis generating tool to identify the most differentially expressed genes in whole blood of infants hospitalised with RSV, subtype B, bronchiolitis, and to further evaluate some of these genes with QRT-PCR. We found the interferon, alpha-inducible protein 27 (IFI27) and Charcot-Leyden crystal protein (CLC) to be the most differentially expressed genes in this study and not to be reported before in conjunction with acute viral bronchiolitis.

All 18 hospitalised infants in our study suffered from acute bronchiolitis of moderate severity. They were hospitalised because of severe respiratory difficulties, but were not in need of assisted ventilation including treatment with continuous positive airway pressure (CPAP) or mechanical ventilation. Our population is therefore representative for the majority of cases hospitalised during infancy worldwide [9]. The children in the control group were all healthy one year old children at the time of clinical examination and blood collection. They were representing the majority of infants undergoing RSV infection with only mild symptoms. This group was chosen to exclude children from the control group that could have a possible genomic predisposition for a more abnormal immune response, characterizing the infants with overt bronchiolitis. We would also be interested in differentiating between the gene expression profile of mild cases and severe cases of RSV infection, but this will require a new study and a design specific to support this question. The gene expression profile found in the present study is the result of RSV infection. It would in further studies be interesting to examine if other respiratory viral infections causing bronchiolitis give a similar or different expression pattern.

One limitation of the present study is the discrepancy in age between the cases and controls. It is well known that the immunological responses mature very rapidly from birth to age one year and the differences observed may therefore simply be due to age rather than to any infection related host response [10]. In order to confirm our gene expression results new studies, which also include a control group being more closely age matched, are very much needed. This study is also rather small, consisting of only 18 infants, and warrants us to be cautious to the interpretations of the results given.

The immune response to RSV infection, as with any other infection is comprised of an innate response, and subsequently by activation of humoral and cellular specific immunity. The cellular immune response seems to be important in controlling the infection once initiated and for clearance of the virus, while the humoral immune response is important in providing protection against subsequent infection [11]. The interferons (IFNs) are a diverse family of cytokines consisting in humans mainly of IFN-alpha and IFN-beta (type I) and IFN-gamma (type II). Both types play an essential role in host immunity by inhibiting the replication and spread of viral, bacterial, and parasitic pathogens. Several studies have shown the importance of IFN-gamma as regulators of host response to RSV infection [12, 13]. However, the sensitivity of RSV to the antiviral activity of INF-alpha is known to be low, at least when compared with that of other viruses [14]. IFNs mediate their effects via transcription of interferon stimulated genes (ISGs) by binding to cell surface receptors activating members of the JAK-STAT pathway [15]. Studies have shown an important role of the RSV surface protein NS1 to antagonize the type I IFN-mediated antiviral response [16]. Also, recent in vitro studies have shown that knocking out the NS1 gene results in upregulation of several ISGs, increased viral clearance and modulation of the host response towards a Th1 phenotype [3, 4, 17]. The gene IFI27 belong to a group of small ISGs [18, 19]. Its biological functions have yet to be revealed, but studies have found this gene to be highly upregulated during cancer development in several tissues and in children with untreated juvenile dermatomyositis [20]. It has also been shown to have a direct antiviral effect against certain viruses [21]. In our study IFI27 was highly upregulated in all but one infant even after evaluation with QRT-PCR. However, the specific biological role(s) of this gene, if any, in early-life RSV bronchiolitis need more studies.

Major basic protein (MBP), eosinophil peroxidase (EPO), eosinophil cationic protein (ECP), eosinophil-derived neurotoxin (EDN) and Charcot-Leyden crystal (CLC) protein are major secretory effector proteins of eosinophils [22]. Eosinophilia of blood and tissue is classically associated with parasitic and especially invasive helminthic infections and is not typical for either bacterial or viral infections. However, several groups have shown evidence of eosinophil degranulation in the lung parenchyma during RSV infection, and both ECP and EDN are shown to have antiviral effects against RSV in vitro [23]. Furthermore, in vivo studies of infants with RSV bronchiolitis have found blood levels of ECP to be higher during convalescence than during acute disease and therefore possibly play a role in the development of the long-standing inflammatory reaction seen in the airways after RSV bronchiolitis, and in which eosinophils play an important role [24, 25]. In our study we found the gene coding for the eosinophil-derived phospholipase CLC to be downregulated in peripheral blood of all infants even after evaluation with QRT-PCR. Studies have shown that eosinophil-derived phospholipases may contribute to surfactant dysfunction in the asthmatic lung [26]. However, the possible importance of CLC downregulation in peripheral blood of infants with RSV bronchiolitis need more studies.

The present study shows the importance of further evaluation of microarray data with other methods such as QRT-PCR and protein expression analyses. Microarray and QRT-PCR are suitable tools for gene expression studies related to acute RSV bronchiolitis in infants and can increase our knowledge of the host response to this disease. The present study identified the genes IFI27 and CLC to be highly differentially expressed in whole blood of infants hospitalised with RSV, subtype B, bronchiolitis of moderate severity. However, more studies are needed to identify the specificity of these two genes in relation to this disease.


Authors’ Affiliations

University of Oslo, Faculty Division Akershus University Hospital, Department of Paediatrics, Akershus University Hospital
Institute of Clinical Epidemiology and Molecular Biology, Akershus University Hospital
Norwegian Computing Center


  1. Ogra PL: Respiratory syncytial virus: the virus, the disease and the immune response. Paediatr Respir Rev. 2004, 5 Suppl A: S119-S126. 10.1016/S1526-0542(04)90023-1.View ArticlePubMedGoogle Scholar
  2. Black CP: Systematic review of the biology and medical management of respiratory syncytial virus infection. Respir Care. 2003, 48: 209-231.PubMedGoogle Scholar
  3. Bitko V, Musiyenko A, Shulyayeva O, Barik S: Inhibition of respiratory viruses by nasally administered siRNA. Nat Med. 2005, 11: 50-55. 10.1038/nm1164.View ArticlePubMedGoogle Scholar
  4. Zhang W, Yang H, Kong X, Mohapatra S, San Juan-Vergara H, Hellermann G, Behera S, Singam R, Lockey RF, Mohapatra SS: Inhibition of respiratory syncytial virus infection with intranasal siRNA nanoparticles targeting the viral NS1 gene. Nat Med. 2005, 11: 56-62. 10.1038/nm1174.View ArticlePubMedGoogle Scholar
  5. Lockhart DJ, Winzeler EA: Genomics, gene expression and DNA arrays. Nature. 2000, 405: 827-836. 10.1038/35015701.View ArticlePubMedGoogle Scholar
  6. Kubista M, Andrade JM, Bengtsson M, Forootan A, Jonak J, Lind K, Sindelka R, Sjoback R, Sjogreen B, Strombom L, Stahlberg A, Zoric N: The real-time polymerase chain reaction. Mol Aspects Med. 2006, 27: 95-125. 10.1016/j.mam.2005.12.007.View ArticlePubMedGoogle Scholar
  7. Ruuskanen O, Ogra PL: Respiratory syncytial virus. Curr Probl Pediatr. 1993, 23: 50-79. 10.1016/0045-9380(93)90003-U.PubMedGoogle Scholar
  8. Kehl SC, Henrickson KJ, Hua W, Fan J: Evaluation of the Hexaplex assay for detection of respiratory viruses in children. J Clin Microbiol. 2001, 39: 1696-1701. 10.1128/JCM.39.5.1696-1701.2001.View ArticlePubMedPubMed CentralGoogle Scholar
  9. Wang EE, Law BJ, Stephens D: Pediatric Investigators Collaborative Network on Infections in Canada (PICNIC) prospective study of risk factors and outcomes in patients hospitalized with respiratory syncytial viral lower respiratory tract infection. J Pediatr. 1995, 126: 212-219. 10.1016/S0022-3476(95)70547-3.View ArticlePubMedGoogle Scholar
  10. Bont L, Kimpen JL: Immunological mechanisms of severe respiratory syncytial virus bronchiolitis. Intensive Care Med. 2002, 28: 616-621. 10.1007/s00134-002-1256-z.View ArticlePubMedGoogle Scholar
  11. Hacking D, Hull J: Respiratory syncytial virus--viral biology and the host response. J Infect. 2002, 45: 18-24. 10.1053/jinf.2002.1015.View ArticlePubMedGoogle Scholar
  12. Aberle JH, Aberle SW, Rebhandl W, Pracher E, Kundi M, Popow-Kraupp T: Decreased interferon-gamma response in respiratory syncytial virus compared to other respiratory viral infections in infants. Clin Exp Immunol. 2004, 137: 146-150. 10.1111/j.1365-2249.2004.02504.x.View ArticlePubMedPubMed CentralGoogle Scholar
  13. Schauer U, Hoffjan S, Rothoeft T, Bartz H, Konig S, Fuchs E, Bittscheidt J, Kochling A, Stephan V: Severe respiratory syncytial virus infections and reduced interferon-gamma generation in vitro. Clin Exp Immunol. 2004, 138: 102-109. 10.1111/j.1365-2249.2004.02582.x.View ArticlePubMedPubMed CentralGoogle Scholar
  14. Guerrero-Plata A, Baron S, Poast JS, Adegboyega PA, Casola A, Garofalo RP: Activity and regulation of alpha interferon in respiratory syncytial virus and human metapneumovirus experimental infections. J Virol. 2005, 79: 10190-10199. 10.1128/JVI.79.16.10190-10199.2005.View ArticlePubMedPubMed CentralGoogle Scholar
  15. Haque SJ, Williams BR: Signal transduction in the interferon system. Semin Oncol. 1998, 25: 14-22.PubMedGoogle Scholar
  16. Bossert B, Conzelmann KK: Respiratory syncytial virus (RSV) nonstructural (NS) proteins as host range determinants: a chimeric bovine RSV with NS genes from human RSV is attenuated in interferon-competent bovine cells. J Virol. 2002, 76: 4287-4293. 10.1128/JVI.76.9.4287-4293.2002.View ArticlePubMedPubMed CentralGoogle Scholar
  17. Conzelmann KK: Transcriptional activation of alpha/beta interferon genes: interference by nonsegmented negative-strand RNA viruses. J Virol. 2005, 79: 5241-5248. 10.1128/JVI.79.9.5241-5248.2005.View ArticlePubMedPubMed CentralGoogle Scholar
  18. Martensen PM, Sogaard TM, Gjermandsen IM, Buttenschon HN, Rossing AB, Bonnevie-Nielsen V, Rosada C, Simonsen JL, Justesen J: The interferon alpha induced protein ISG12 is localized to the nuclear membrane. Eur J Biochem. 2001, 268: 5947-5954. 10.1046/j.0014-2956.2001.02545.x.View ArticlePubMedGoogle Scholar
  19. Parker N, Porter AC: Identification of a novel gene family that includes the interferon-inducible human genes 6-16 and ISG12. BMC Genomics. 2004, 5: 8-10.1186/1471-2164-5-8.View ArticlePubMedPubMed CentralGoogle Scholar
  20. Martensen PM, Justesen J: Small ISGs coming forward. J Interferon Cytokine Res. 2004, 24: 1-19. 10.1089/107999004772719864.View ArticlePubMedGoogle Scholar
  21. Labrada L, Liang XH, Zheng W, Johnston C, Levine B: Age-dependent resistance to lethal alphavirus encephalitis in mice: analysis of gene expression in the central nervous system and identification of a novel interferon-inducible protective gene, mouse ISG12. J Virol. 2002, 76: 11688-11703. 10.1128/JVI.76.22.11688-11703.2002.View ArticlePubMedPubMed CentralGoogle Scholar
  22. Egesten A, Calafat J, Janssen H, Knol EF, Malm J, Persson T: Granules of human eosinophilic leucocytes and their mobilization. Clin Exp Allergy. 2001, 31: 1173-1188. 10.1046/j.1365-2222.2001.01138.x.View ArticlePubMedGoogle Scholar
  23. Domachowske JB, Rosenberg HF: Respiratory syncytial virus infection: immune response, immunopathogenesis, and treatment. Clin Microbiol Rev. 1999, 12: 298-309.PubMedPubMed CentralGoogle Scholar
  24. Smyth RL, Fletcher JN, Thomas HM, Hart CA: Immunological responses to respiratory syncytial virus infection in infancy. Arch Dis Child. 1997, 76: 210-214.View ArticlePubMedPubMed CentralGoogle Scholar
  25. Priftis KN, Papadopoulou A, Liatsis E, Katsikas D, Nicolaidou P, Kanariou M: Serum eosinophil cationic protein and CD23 in acute RSV bronchiolitis. Med Sci Monit. 2005, 11: CR493-CR497.PubMedGoogle Scholar
  26. Ackerman SJ, Kwatia MA, Doyle CB, Enhorning G: Hydrolysis of surfactant phospholipids catalyzed by phospholipase A2 and eosinophil lysophospholipases causes surfactant dysfunction: a mechanism for small airway closure in asthma. Chest. 2003, 123: 355S-10.1378/chest.123.3_suppl.355S.View ArticlePubMedGoogle Scholar
  27. Pre-publication history

    1. The pre-publication history for this paper can be accessed here:


© Fjaerli et al; licensee BioMed Central Ltd. 2006

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 (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.