key: cord-0985907-ap0zwzzl authors: Herberg, Jethro A.; Kaforou, Myrsini; Gormley, Stuart; Sumner, Edward R.; Patel, Sanjay; Jones, Kelsey D. J.; Paulus, Stéphane; Fink, Colin; Martinon-Torres, Federico; Montana, Giovanni; Wright, Victoria J.; Levin, Michael title: Transcriptomic Profiling in Childhood H1N1/09 Influenza Reveals Reduced Expression of Protein Synthesis Genes date: 2013-07-29 journal: The Journal of Infectious Diseases DOI: 10.1093/infdis/jit348 sha: d9d57e6e5eaf077c1223b5295b461b7f2f1cfc82 doc_id: 985907 cord_uid: ap0zwzzl We compared the blood RNA transcriptome of children hospitalized with influenza A H1N1/09, respiratory syncytial virus (RSV) or bacterial infection, and healthy controls. Compared to controls, H1N1/09 patients showed increased expression of inflammatory pathway genes and reduced expression of adaptive immune pathway genes. This was validated on an independent cohort. The most significant function distinguishing H1N1/09 patients from controls was protein synthesis, with reduced gene expression. Reduced expression of protein synthesis genes also characterized the H1N1/09 expression profile compared to children with RSV and bacterial infection, suggesting that this is a key component of the pathophysiological response in children hospitalized with H1N1/09 infection. Between July 2009 and June 2010, we recruited 165 febrile children (<17 years) with respiratory infection admitted to St Mary's Hospital London. Whole blood was collected in PAXgene tubes, together with routine clinical samples. Patients were recruited early in their admission, before diagnostic studies were completed, and were assigned to diagnostic categories once results were available (Supplementary Figure 1) . The study was approved by St Mary's Research Ethics Committee (REC 09/H0712/58). Written, informed consent was obtained. Viral diagnostics were undertaken on nasopharyngeal aspirates using immunofluorescence (respiratory syncytial virus [RSV] , adenovirus, parainfluenza virus, influenza A + B) and nested polymerase chain reaction (PCR; RSV, coronavirus, adenovirus, parainfluenza 1-4, influenza A + B, bocavirus, metapneumovirus, rhinovirus). Bacterial cultures included blood, cerebrospinal fluid (CSF), urine, and tissue sites. Pneumococcal antigen was measured in blood and urine, and bacterial DNA was detected by meningococcal and pneumococcal PCR. PAXgene tubes were extracted using PAXgene Blood RNA extraction kits (Qiagen). After quantification and quality control, biotin-labeled cRNA was prepared using Illumina TotalPrep RNA Amplification kits (Applied Biosystems). 750 ng labeled cRNA was hybridized to Illumina HumanHT-12v3 BeadChips, and the microarrays scanned. Quality control parameters were assessed using Genome Studio software, visual inspection of the microarray images and principal component analysis (PCA). Two children in the H1N1/09, RSV, and bacterial cohorts had oral or intravenous steroids prior to blood sampling but were not outliers using PCA or unsupervised clustering algorithms (data not shown). Expression data were analyzed using 'R' Language and Environment for Statistical Computing 2.12.1 and GeneSpringGX 11.5 software (Agilent), as described (Supplementary Methods). The data have been deposited in NCBI's Gene Expression Omnibus and are accessible through Series accession number GSE42026 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE42026). We analyzed 800 network-associated molecules using IPA (Ingenuity Systems, www.ingenuity.com) to identify the most significant biological pathways for differentially expressed genes in each group. Differences between cohorts in age, proportion of white cell type, and days from illness onset to presentation were calculated using Kruskall-Wallis tests, with Dunn post-test comparisons of each cohort (Graphpad Prism). Differences in gender, severity of illness, and deaths were calculated using χ 2 tests of independence (Supplementary Table 1 ). Our consecutively recruited cohort included 25 patients with H1N1/09 infection (6 coinfected with other pathogens), 34 with RSV infection (11 coinfected) , and 18 with bacterial infection (Supplementary Figure 1) . Thirty-three healthy control children were recruited from the outpatients department. Patient demographic and clinical data are summarized in Supplementary Table 1 . When the admission blood RNA expression profile was compared in H1N1/09 patients and controls, we found 1267 significantly differentially expressed (SDE) transcripts (P < .001) (Supplementary Data 1 and 2A) that separated patients and controls on unsupervised clustering analysis ( Figure 1A ). Using Support Vector Machine (SVM) with leave-one-out validation, strong class prediction sensitivity was achieved (96%; Supplementary Data 1). When patients with RSV or bacterial infection were compared to controls, we identified 1172 and 1869 SDE transcripts (P < .001) respectively (Supplementary Data 1 and 2B and 2C). SVM modeling of RNA signatures had sensitivity of 95% and 98% for distinguishing RSV-and bacterial-infected patients respectively (Supplementary Data 1). As the RSV patients were significantly younger than the control or H1N1/09 patients, we compared the RNA expression signature of RSV patients to a subgroup of the youngest controls (n = 10, median age 9 months). The 2411 SDE genes identified (corrected P < .05) were concordant with those identified using all controls: 772 of the top 800 SDE transcripts between RSV and youngest controls were represented, with the same direction of expression, in both data sets. The biological pathways distinguishing H1N1/09, RSV, and bacterial infections from controls were investigated with IPA, which assigns SDE transcripts to known biological functions and pathways. Biological functions of SDE were enriched for infectious disease, respiratory disease and inflammatory response genes in all 3 groups (Supplementary Data 3). In the H1N1/09 vs controls data set, the functional category "protein synthesis" was significantly over-represented (P = 7.93 × 10 −26 ). Furthermore, the 3 most significant canonical pathways were related to initiation of protein translation: eukaryotic initiation factor 2 (eIF2 -P = 3.16 × 10 −25 ), eukaryotic initiation factor 4 (eIF4 -P = 3.24 × 10 −88 ), and mammalian target of rapamycin (mTOR -P = 1.91 × 10 −6 ) (Supplementary Data 4). Transcripts with increased abundance were enriched for immune response pathways and functions, including signaling from pattern recognition receptors (top increased genes OAS3, C1QB, C1QC), interferon signaling (IFI35, IFIT1, IFT3), and antigen presentation through MHC class I (TAP1, MR1, HLA-B). Transcripts with reduced abundance were enriched for pathways including antigen presentation through MHC class II (most decreased: HLA-DQA, HLA-DPB, HLA-DOA), T-cell signaling and protein synthesis -including those involved in translation, elongation, and ribosomal transport (Supplementary Data 4). There was reduced expression in 63 of 69 SDE transcripts in the eIF2 pathway, P < .0001. Protein synthesis-related transcripts with increased expression included genes with known inhibitory functions in protein synthesis, including EIF2AK2, IFIT1/ISG56, and IFIT2/ISG54 (fold change [FC] of 3.9, 10.0, and 4.0, respectively, in H1N1/09 relative to controls) [5, 6] . To validate our findings, we used a publically available adult data set for H1N1/09 and controls [4] . There was complete segregation of H1N1/09 patients and controls when SDE transcripts from our data were used for clustering (Supplementary Figure 2A) . When using the SVM prediction model trained on our H1N1/09 vs control data, and applying this to the independent data set, there was 100% sensitivity and specificity of disease and control status. Of 69 transcripts in the top canonical pathway (EIF2 signaling), 32 were SDE in the validation data set (31 with the same direction of regulation; Supplementary Figure 2B ). To identify unique aspects of the transcription profile of H1N1/09 infection, we compared expression in H1N1/09, RSV, and bacterial infection. Transcripts that were SDE between H1N1/09 and RSV (n = 601 of 7295 transcripts, corrected P < .01) and between H1N1/09 and bacterial infection (n = 734 of 6808 transcripts, corrected P < .01) were used for unsupervised clustering and for classification, H1N1/09 patients were distinguished from RSV ( Figure 1B ) and bacterial patients ( Figure 1C ) with accuracy of 93% and 92% for H1N1/09 vs RSV and H1N1/09 vs bacteria, respectively (Supplementary Data 1 and 2D and 2E). "Protein synthesis" was the most significant functional category in both the H1N1/09 vs RSV, and H1N1/09 vs bacterial infection data sets (corrected P = 2.73 × 10 −13 and 1.6 × 10 −11 , respectively). The most significant canonical pathways were related to protein synthesis: eIF2 pathway (P = 2.0 × 10 −20 and 9.0 × 10 −10 ), eIF4 (P = 4.7 × 10 −8 and 0.006), and mTOR (P = 8.1 × 10 −7 and 0.016) for H1N1/09 vs RSV, and H1N1/09 vs bacterial infection, respectively (Figure 2 and Supplementary Data 3 and 4) . Marked differences were observed between RNA transcriptional profiles in H1N1/09-infected patients compared to RSV and gram-positive bacterial infection or controls. Patients with H1N1/09 showed increased expression in well-established pathways of the innate antiviral immune response, including pattern recognition receptor, interferon signaling, and myeloid cell activation pathways. Many interferon-induced genes with anti-viral roles showed increased differential expression, including IFI27 (P = 5.8 × 10 −11 , FC 72.3), and genes implicated in repression of protein synthesis including EIF2AK2, and also IFIT1 and IFIT2, whose anti-viral roles may be mediated through effects on translation [7, 8] . In vitro and murine studies have implicated IFITM3 in restricting influenza replication [9] . Our finding of increased IFITM3 expression (P = 3 × 10 −3 , FC 6.25) supports its role in H1N1/09 infection. Differences in whole blood RNA expression can reflect altered proportions of cell types and/or altered transcript abundance. Lymphocyte proportion differed significantly only between H1N1/09 and controls (fold change 2.1), whereas the fold change of the top 800 SDE genes used for Ingenuity varied between −6.2 and 72.2. Therefore, differences in lymphocyte proportion do not explain the observed differences, although a role for lymphocyte subset populations is not excluded. H1N1/09 patients showed reduced expression in adaptive immune response pathways, including T-cell activation through NFAT, B cell activation, and MHC class II antigen presentation. Reduced expression of adaptive T-and B-cell pathways may represent a host negative regulatory feedback mechanism [3] , or viral-induced immune subversion [10] , analogous to viral inhibition of antigen presentation by inhibition of the MHC class II transactivator. The 3 most significant pathways are involved in protein synthesis and translation initiation (eIF2, eIF4, and mTOR), and within these, most genes showed reduced expression. In contrast, inhibitors of translation including EIF2AK2 (which encodes PKR) were upregulated. PKR phosphorylates eukaryotic initiation factor 2 alpha (eIF2a) to cause the arrest of translation initiation [11] , suggesting that protein synthesis in H1N1/09 patients may be impaired by repressed transcription of genes required for the translation apparatus, and by increased transcription of translation inhibitors. Reduced expression of protein synthesis genes was significant between H1N1/09 and controls (P = 7.9 × 10 −26 ) and also between H1N1/09 and RSV or bacterial infection (P = 2.73 × 10 −13 and 1.6 × 10 −11 , respectively; Figure 2 ). In an influenza A challenge study of adult volunteers, differential expression of translation initiation pathway genes was identified in symptomatic but not subclinical infections [12] , supporting a role for altered protein synthesis as a key component of the host response to H1N1/09 infection. Increased expression of protein synthesis-inhibiting genes with reduced expression of translation initiation genes is likely to impair protein synthesis. Influenza virus has been observed to shut-off protein synthesis and decrease mRNA levels in infected cells [13] , but the phenomenon we have observed in peripheral blood suggests a widespread alteration affecting cells distant from the site of infection. We hypothesize that suppression of protein synthesis may be a distinctive feature of H1N1/09 infection and may play a role in the prolonged debility, fatigue, and delayed recovery seen in severely affected patients. Many of the fatal H1N1/09 cases had acute respiratory distress syndrome requiring prolonged ventilation [1] , and late pulmonary hemorrhage or persistent respiratory failure, with death weeks after the acute illness. Suppressed protein synthesis might result in poor tissue recovery and may provide an explanation for the prolonged illness and poor outcome in some patients. In pigs, sepsis impairs protein synthesis through inhibition of mTOR and eIF2 pathways [14] . In culture, H5N1 influenza virus impairs mTOR signaling leading to autophagic cell death Figure 2 . Top canonical pathways differing between H1N1/09 and controls, RSV and bacterial infection. Each bar is colored in proportion to the number of SDE H1N1/09 transcripts increased (red) or decreased (green) in abundance relative to the comparator cohort. The total bar length is proportional to P value. Colored blocks next to each pathway are coded according to biological function. Protein synthesis pathways (yellow) were the most significant in all 3 comparisons, with predominant decreased expression in H1N1/09 patients relative to the comparator group. Innate immune pathway transcripts (brown) were increased in H1N1/09 patients, whereas adaptive immune transcripts (pink) were reduced relative to controls. [15] . Further studies are required to establish whether the changes we observed in peripheral blood RNA expression are associated with inhibition of protein synthesis and to confirm the role of this phenomenon in severe illness. The changes in expression of genes involved in protein synthesis we have observed may reflect a response to viral components entering the bloodstream, or an indirect effect mediated by host inflammatory molecules. There may be advantages to both virus and host in inhibiting protein synthesis, as immune subversion may occur through impaired leukocyte function, and limitation of production of viral proteins may favor recovery. Limitations in this study include the younger age of the RSV cohort and ethnic variation in the population studied, which was not matched between groups. Furthermore changes in RNA expression during disease may result from altered cell type proportions, as well as upregulation or downregulation of genes. The association of influenza with bacterial infection is well known, and occult bacterial infection may have contributed to disease in some patients. Our comparative study of RNA expression in children with H1N1/09 infection provides information on host response to the pandemic strain. The finding that altered expression of genes involved in protein synthesis is a key feature distinguishing H1N1/09 from other common childhood infections offers a new avenue for investigating the immunopathogenesis of influenza. Comparison of pandemic and seasonal influenza reveals higher mortality and increased prevalence of shock in children with severe H1N1/09 infection Gene expression patterns in blood leukocytes discriminate patients with acute infections Host adaptive immunity deficiency in severe pandemic influenza A distinct influenza infection signature in the blood transcriptome of patients with severe community-acquired pneumonia A new pathway of translational regulation mediated by eukaryotic initiation factor 3 Distinct induction patterns and functions of two closely related interferon-inducible human genes, ISG54 and ISG56 ISG56/IFIT1 is primarily responsible for IFN-induced changes in patterns of PIV5 transcription and protein synthesis 2'-O methylation of the viral mRNA cap evades host restriction by IFIT family members IFITM3 restricts the morbidity and mortality associated with influenza Pandemic (H1N1) 2009 influenza virus induces weaker host immune responses in vitro: a possible mechanism of high transmissibility PERK and PKR: old kinases learn new tricks Temporal dynamics of host molecular responses differentiate symptomatic and asymptomatic influenza A infection Inhibition of host protein synthesis and degradation of cellular mRNAs during infection by influenza and herpes simplex virus Sepsis and development impede muscle protein synthesis in neonatal pigs by different ribosomal mechanisms Avian influenza A virus H5N1 causes autophagy-mediated cell death through suppression of mTOR signaling Acknowledgements. We thank the clinicians and patients for taking part, Dr Jayne Dennis who gave invaluable microarray advice, and Mary Cross who gave clinical research support.Financial support. This work was supported by the Imperial College Comprehensive Biomedical Research Centre and the Wellcome Trust Centre for Respiratory Infection at Imperial College (DMPED P26077 to FMT) and the Fondo de Investigación Sanitaria del plan nacional de I+D+I and 'fondos FEDER (F. I. S.; PI070069/PI1000540 to FMT). FMT represents the GENDRES research group (www.gendres.org).Potential conflicts of interest. All authors: No reported conflicts.All authors have submitted the ICMJE Form for Disclosure of Potential Conflicts of Interest. Conflicts that the editors consider relevant to the content of the manuscript have been disclosed.