Respiratory Mucosal Proteome Quantification in Human Influenza Infections

Abstract

Respiratory influenza virus infections represent a serious threat to human health. Underlying medical conditions and genetic make-up predispose some influenza patients to more severe forms of disease. To date, only a few studies have been performed in patients to correlate a selected group of cytokines and chemokines with influenza infection. Therefore, we evaluated the potential of a novel multiplex micro-proteomics technology, SOMAscan, to quantify proteins in the respiratory mucosa of influenza A and B infected individuals. The analysis included but was not limited to quantification of cytokines and chemokines detected in previous studies. SOMAscan quantified more than 1,000 secreted proteins in small nasal wash volumes from infected and healthy individuals. Our results illustrate the utility of micro-proteomic technology for analysis of proteins in small volumes of respiratory mucosal samples. Furthermore, when we compared nasal wash samples from influenza-infected patients with viral load ≥ 2⁸ and increased IL-6 and CXCL10 to healthy controls, we identified 162 differentially-expressed proteins between the two groups. This number greatly exceeds the number of DEPs identified in previous studies in human influenza patients. Most of the identified proteins were associated with the host immune response to infection, and changes in protein levels of 151 of the DEPs were significantly correlated with viral load. Most important, SOMAscan identified differentially expressed proteins heretofore not associated with respiratory influenza infection in humans. Our study is the first report for the use of SOMAscan to screen nasal secretions. It establishes a precedent for micro-proteomic quantification of proteins that reflect ongoing response to respiratory infection.

About the Authors:

Tony Marion

Affiliation: University of Tennessee Health Science Center, Department of Microbiology, Immunology and Biochemistry, Memphis, United States of America

Husni Elbahesh

Affiliation: University of Tennessee Health Science Center, Department of Microbiology, Immunology and Biochemistry, Memphis, United States of America

Paul G. Thomas

Affiliation: Department of Immunology, St. Jude Children’s Research Hospital, Memphis, United States of America

John P. DeVincenzo

Affiliations University of Tennessee Health Science Center, Department of Microbiology, Immunology and Biochemistry, Memphis, United States of America, University of Tennessee Health Science Center, Department of Pediatrics, Memphis, United States of America, Children’s Foundation Research Center at Le Bonheur Children’s Hospital, Memphis, United States of America

Richard Webby

Affiliation: Department of Infectious Diseases, St Jude Children’s Research Hospital, Memphis, United States of America

Klaus Schughart

* E-mail: [email protected]

Affiliations University of Tennessee Health Science Center, Department of Microbiology, Immunology and Biochemistry, Memphis, United States of America, Department of Infection Genetics, Helmholtz Centre for Infection Research, Braunschweig, Germany, University of Veterinary Medicine Hannover, Hannover, Germany

Introduction

Each year, about 500 million people are infected by the influenza A virus (IAV) worldwide, of which about 500,000 die [1]. In recent history, the emergence of new influenza subtypes has caused several pandemics [2–4]. The most severe pandemic in 1918 caused about 30–50 million deaths worldwide [5], and a new variant of a seasonal H1N1 virus, pH1N1, caused a world-wide pandemic in 2009 [6–8]. Avian viruses can also directly infect humans. In particular, two subtypes, H5N1 and H7N9, may cause severe disease with lethal outcome [9–13]. Adverse health conditions, such as obesity and diabetes, and genetic make-up predispose influenza patients to more severe forms of the disease [14–19]. Cytokines and chemokines released in a cytokine storm in response to influenza infection contribute to disease severity [20]. Unraveling the pathogenesis of influenza in humans so as to identify potential targets for human therapeutics and predictors of disease severity necessitates the evaluation of the main site of viral replication, the mucosal tissues of the respiratory tract. The majority of the disease pathogenesis caused by influenza occurs after viral replication has already started to decline [21], thus adding to the impetus to develop host-response-targeted therapies in addition to continuing evaluation of better antiviral therapeutics. Additionally, host-response-based diagnostics may improve identification of patients at highest risk of disease progression. Quantitative mucosal biomarker identification is important for such work to proceed rationally.

Most multiplex assays for disease protein biomarkers in inflammation and infection have been limited to detection of chemokines and cytokines expected to play a role in disease pathogenesis and for which prepared kits are readily available. Hence those assays have inherent bias. Here, we have attempted a new approach to biomarker identification in influenza infected patients using an aptamer-dependent, micro-proteomic approach (SOMAscan^®). SOMAscan is a recently developed technology that can simultaneously quantify more than 1,000 human proteins in small volumes of complex biological fluids [22]. We used the SOMAscan version 1.2k that had a multiplex library of 1,129 SOMAmers (Slow Off-rate Modified nucleic acid based Aptamers) that each quantifies a single soluble protein [23]. SOMAscan transforms the number of each protein-bound SOMAmer into a quantitative measure of protein concentration [24]. SOMAscan is highly sensitive with a threshold of detection of 30 femtomolar, <1 pg/ml and 10⁸-fold dynamic range for quantification of proteome changes in mice [25] and humans [22, 24, 26, 27].

Here, we identified several differentially expressed proteins in mucosal secretions that heretofore have not been associated with respiratory influenza infection in humans. Our results indicate that SOMAscan is well suited for biomarker discovery in respiratory infections.

Results

SOMAscan detects a broad spectrum of proteins in mucosal secretions

We first evaluated whether an unbiased proteomic screen like SOMAscan would identify proteins in mucosal samples other than those detected in standard cytokine and chemokine assays [28, 29] and already known to be relevant for inflammation and immune response to respiratory influenza infection. For this, we analyzed 24 nasal wash samples from a previously described influenza cohort [28] (S1 Table). The samples included 18 from patients with identified influenza infection, 11 patients infected with IAV H3 strain, two infected with IAV pH1 subtypes, and five infected with influenza B. The healthy control group consisted of samples from six individuals who were not infected with influenza, who showed no clinical signs of infection, such as fever, cough, or sore throat [28] and who appeared to be otherwise healthy. Viral loads in infected patients, determined by quantitative reverse transcription PCR (qRT-PCR) [28], ranged from 0.72 to 5.1x10⁶ copies/ml [28].

The nasal washes were analyzed with the 1.2k version of SOMAscan with 1,129 SOMAmers that simultaneously quantified 1,030 different human proteins [24]. Quantitative expression signals were log₂-transformed and quantile-normalized and then examined by principle component analysis (PCA) to visualize variation and grouping of samples. PCA is a mathematical transformation that reduces variation in a large data set to a few dimensions that project differences determined by the strongest variables. The first two dimensions that represent the highest variations of a PCA analysis can be visualized in a 2D plot. In this way, samples that are most disparate with respect to protein expression levels are found more distantly located from each other in the 2D PCA plot, and samples with similar protein expression levels are closer to each other in the plot. As shown in Fig 1, most infected patients (red circles) segregated together indicating that their protein expression levels were similar. Most influenza virus-positive patients segregated to the right of PCA1 = 0 suggesting that these patients had similar changes in their mucosal proteome compared to others in the cohort. The healthy controls (grey circles) segregated to the left except for ID_3058 that segregated with the leftmost group of virus-infected patients. Three patients who were diagnosed as virus-positive (red circles), based on qRT-PCR [28], segregated with healthy controls (grey circles). Patient ID_3045 segregated independently of all other patients. The imperfect PCA distribution of infected versus healthy patients is in part due to the small sample size but also reflects the large uncontrollable heterogeneity which is an intrinsic property of human cohorts. Here, many confounding factors, such as adverse health conditions, time after infection, genetics, and general environment, may influence mucosal protein expression levels and are difficult or impossible to control.

[Figure omitted. See PDF.]

Fig 1. PCA analysis of normalized protein expression values.

Principle component analysis (PCA) was performed with quantile normalized log₂–transformed protein expression values from nasal washes for all 24 samples. The first two principal components are shown representing 24% and 17%, respectively, of the total variation. Healthy controls are labeled gray and IAV-positive samples (in which influenza A or B was detected by PCR) are labeled red. In addition, sample identities (e.g. ID_4043) are shown. Horizontal and vertical axis represent principle component 1 and 2, respectively.

https://doi.org/10.1371/journal.pone.0153674.g001

There was good correlation between measurements in SOMAscan for four analytes that were previously quantified by Luminex in samples from the same patients [28]. Relative protein expression levels for IL-6, CXCL10/IP10, CCL7/MCP3, and CXCL8/IL-8 were highest among proteins measured by SOMAscan, and those levels correlated well with measurements of the same proteins by Luminex (S1 Fig). The majority of proteins quantified by SOMAscan were not measured in the previous Luminex study. The correlation between Luminex and SOMAscan validates our results and verifies SOMAscan's potential to identify proteins relevant to pathophysiology, inflammation, and immune response in mucosal secretions after influenza infection.

Differentially expressed proteins in influenza virus infected patients

A comparison of all 18 samples from infected patients with the six samples from healthy individuals did not identify statistically significant differentially expressed proteins (DEPs). This result is consistent with the incomplete segregation of positive and negative samples in the PCA (Fig 1) and previous chemokine and cytokine measurements in samples from the entire patient cohort from which the samples analyzed herein were obtained [28]. The low viral load in some of the infected patients (S1 Table) may account for this result. It is well known from animal experiments that immune responses and pathologies depend largely on viral loads (reviewed in [30]). Furthermore, the distribution of viral loads among all of the infected patients in the entire patient cohort [28] was bimodal (S2 Fig). Based upon the bimodal distribution between patients with higher versus lower viral RNA levels in the entire influenza infected patient cohort, we selected a subset of patients with a minimum viral load of 2⁸ (subset A, S2 Table) to identify proteins that were differentially expressed between infected patients (12 samples) and healthy controls (6 samples). Fig 2 shows the PCA for this subset. Comparison of the SOMAscan results between virally-infected patients and healthy controls from subset A by LIMMA revealed a total of 23 DEPs, 5 increased and 18 decreased (log₂-fold change ≥ 1 and adjusted p < 0.01) (S3 Table).

[Figure omitted. See PDF.]

Fig 2. PCA analysis of normalized protein expression values of subset A.

Principle component analysis (PCA) was performed with normalized log₂–transformed protein expression values from nasal washes for 18 samples from the subset A that was selected based viral load being higher than > 2⁸. The first two principal components are shown representing 42% and 16%, respectively, of the total variation. Healthy controls are labelled gray and IAV-positive samples are labelled red. In addition, sample identities (e.g., ID_4002) are shown. Horizontal and vertical axis represent principle component 1 and 2, respectively.

https://doi.org/10.1371/journal.pone.0153674.g002

In subset A, sample ID_4043 segregated separately in the PCA (Fig 2) from the other influenza virus-positive samples and had low levels of CXCL10 and IL6 (data not shown). CXCL10 and IL6 are known to be highly elevated after influenza infection [29], and the previous cytokine-chemokine analysis [28] indicated that plasma CXCL10 was correlated with viral loads. Therefore, we analyzed a second subset (subset B) which was identical to subset A except that the outlier sample ID_4043 was deleted. Thus, subset B (S2 Table) consisted of six healthy controls (same as in subset A) and 11 influenza-infected patients (12 minus 1 from subset A). The corresponding PCA showed good segregation between uninfected control and influenza virus-positive groups (S3 Fig). Accordingly, analyses of IL6 and CXCL10 expression levels showed significant differences between the two groups (Kruskal-Wallis, p < 0.01 and p < 0.05, respectively) (S4 Fig), even though samples from two healthy controls had relatively elevated CXCL10 compared to the other healthy controls. Further analysis of subset B using LIMMA revealed a total of 162 DEPs, 63 increased and 99 decreased in samples from virally infected patients compared to uninfected individuals (log₂-fold change ≥ 1 and adjusted p < 0.01) (S4 Table). Functional analysis of these DEPs using the Reactome pathway database showed enrichment for pathways that are activated during the host immune response (Fig 3). Furthermore, 151 of the 162 (up- or down-regulated) DEPs were also highly correlated with viral load. 56 of the 151 DEPs were positively and 95 were negatively correlated with viral load (0.6 < ρ <—0.6 and FDR < 0.01) (Fig 4, S5 Table).

[Figure omitted. See PDF.]

Fig 3. Functional analysis of DEPs from subset B.

Pathway enrichment analysis for normalized log₂–transformed expression values of proteins that were differentially expressed in IAV-positive versus IAV-negative patients from subset B is presented. Pathway terms that were enriched for the DEPs from subset B are indicated on the y-axis, the number of DEPs in the respective pathway category is indicated on the x-axis. The p-value for the probability that the observed distribution of expression occurred by chance is represented by colors of bars. The cut-off value for the pathway p-values was chosen at 0.05.

https://doi.org/10.1371/journal.pone.0153674.g003

[Figure omitted. See PDF.]

Fig 4. Scatter plot of RPS6KA5 expression levels and viral load.

Scatter plot of normalized log₂–transformed relative protein concentrations of RPS6KA5 which was the DEP that most strongly correlated with viral titers and log₂-transformed virus titers are presented for patients from subset B. Red line: linear regression model.

https://doi.org/10.1371/journal.pone.0153674.g004

In conclusion, the SOMAscan analysis revealed the expected large heterogeneity among the 24 samples, similar to selective cytokine-chemokine multiplex assays, and reproduced results obtained previously with Luminex assays. In addition SOMAscan identified many additional DEPs that could not be detected in the previous study of the influenza cohort from which our samples were derived. Thus, SOMAscan is well suited for the detection of proteome differences in patients’ responses to influenza virus infection.

Discussion

In this study, we demonstrate the utility of an unbiased, micro-proteomic analysis from small volumes of nasal washes to quantify proteome differences between influenza virus-infected patients and healthy individuals. SOMAscan identified more than 160 differentially expressed proteins (DEPs) from human respiratory influenza patients with a minimum viral load of 2⁸. Not surprising, most of the DEPs are associated with inflammation and/or immune response pathways. Among the 160 DEPs, 151 were correlated with viral loads in the respective subpopulation of influenza virus-infected patients

Among the DEPs with increased expression detected in our study, many are associated with the recruitment and activation of memory T cells. CXCL11 [31], CXCL10 [32, 33], IL-16 [34, 35], and GNLY (granulysisn) [36], are all chemokines chemotactic for activated and memory CD4 T cells, especially Th1 [33]. Another DEP that was correlated with viral load is amphiregulin. Amphiregulin is an epithelial and fibroblast cell mitogen and growth factor that, chronically, is also associated with inflammation [37, 38]. Recent results from lung influenza infection of mice point to a critical role for amphiregulin in limiting inflammation or otherwise protecting lung epithelium from damage after influenza virus infection [39–41]. Both innate lymphocytes (ILC) [41] and regulatory T cells [39] contributed to local amphiregulin in infected mouse lungs. It will be interesting to know what cells produce amphiregulin and whether it may be a biomarker to predict disease severity in human influenza patients.

The three DEPs with highest correlation to viral load in subset B were RPS6KA5 (MSK1, mitogen and stress activated kinase 1), KNG1 (kininogen), MAP2K1 (MEK-1, mitogen-activated protein kinase kinase-1), ANXA1 (annexin A1), and STIP1 (stress induced phosphoprotein 1). All of these proteins are associated with the initiation or resolution of inflammation and/or stress response and were strongly correlated with viral load. RPS6KA5 (MSK1) had the highest correlation with viral load in subset B (rho = 0.91, p < 0.01). MSK1 has been associated with activation [42–44] and inhibition of inflammation [44, 45] in lung epithelium. MSK1 is downstream of ERK1/2 and p38^MAPK and can phosphorylate NFκB p65 and is a key regulator of lung airway inflammation. MEK-1 is upstream of ERK1/2 and MSK-1 in the non-canonical pathway for p65/RelA phosphorylation by MSK1 leading to NFκB activated transcription of inflammatory cytokines and chemokines [46] including those detected as DEPs described above. MSK1 phosphorylation of p65/RelA is also critical for inflammatory cytokine production in an experimental model for respiratory syncytial virus (RSV) infection [47]. Herein is the first reported association of this key regulator of inflammation with respiratory influenza infection. Kininogen is a co-factor in the plasma kinin-kallikrein, contact system for blood coagulation and the generation of inflammatory mediators [48]. High molecular weight kininogen (HMWK) is precurosor for bradykinin, which is an important mediator of endothelial permeability. Hantavirus infected endothelial cells induce increased cleavage of HMWK with increased bradykinin and vascular permeability [49]. Increased nasal secretions have been correlated with increased glandular bradykinin in experimental infection of adult humans with IAV H1N1 [50]. Annexin A1 functions as a resolvin to reverse or resolve inflammation after it is initiated [51].

Most previous studies of respiratory influenza infection in humans have limited analyses of soluble proteins in serum or nasal washes to cytokines and chemokines predicted to be relevant to inflammation and immune response (reviewed in [52–54]). For example, the previous study of the large patient cohort from which our samples were derived [28] investigated changes in 41 cytokines and chemokines using Luminex. Distinct changes for IFNA2, CCL7 (MCP-3), IL6, and VEGFA (VEGF) in nasal washes and CCL7 (MCP-3), IL10, IL6, and CXCL10 (IP-10) in peripheral blood of human patients were significantly correlated with disease severity [28] but not viral load. IL6 was also identified as a DEP in our study; however, it should be noted that we used a more stringent p = 0.01. When p = 0.05, CCL7, IL6 and CXCL10 are also DEPs in our study. Another study compared plasma cytokines, metalloproteases, and complement activation in seven patients infected with H1N1 and found 16 proteins that were differentially expressed in the blood of influenza infected patients [29]. Five of the 16 proteins, CXCL10, IL6, IL17D, CSF2, PDGFB, from that study were also identified as DEPs in our study (adjusted p = 0.05). Moreover, elevated levels of CXCL10 (IP-10), IL6, IL17, and IL2 in sera of 16 hospital-admitted patients infected with H7N9 were previously reported, and IL6 and CXCL10 correlated with the 13 patients that developed severe disease [55]. A more recent study identified a correlation between serum levels of cytokines MIF, SCF, MCP-1, HGF, and SCGF-β and disease severity after H7N9 respiratory influenza infection [56]. We also found IL17D, IL6, and CXCL10 to be elevated DEPs in influenza virus-infected patients in our study.

A number of global proteomic analyses have been performed with influenza-infected cells in vitro [57–64] and mouse lungs ex vivo [65]. A recent global proteomic analysis of plasma proteins associated with influenza infection used two-dimensional gel electrophoresis and liquid chromatography mass spectrometry (LC-MS) [66]. Only FAM157A, leucine-rich alpha 2 glycoprotein, serum amyloid A protein, and dual oxidase 1 were statistically significantly different in plasma of healthy versus influenza-infected individuals by LC-MS.

In addition to up-regulated proteins, our study also identified many proteins with reduced expression in influenza virus-infected patients versus healthy controls. Of the proteins with reduced expression identified in our study, the five most-reduced DEPs in samples from influenza-infected patients that were also correlated with higher viral loads were CTSD, KLK7, MFGE8, MAPK9 and CD27, respectively. Cathepsin D (CTSD) is a lysosomal aspartic protease that has been shown to be highly expressed in alveolar resident macrophages and Kupffer cells following IAV infections in mice [67, 68]. KLK7 is a chymotrypsin-like serine protease that cleaves proteins at tyrosine, phenylalanine or leucine [69]. It is expressed in the epithelium of the upper and lower respiratory tract (nose, paranasal sinuses, larynx, trachea, bronchial tree) and in their submucosal glands but not the alveolar epithelium [70]. MFGE8 (lactadherin) is a cell adhesion protein that has been suggested to connect smooth muscle cells and elastic fibers of arteries [71]. Importantly, MFGE8 has been shown to promote opsonization of apoptotic cells and their subsequent phagocytosis through its phosphatidylserine (PS) and RGD integrin binding domains, respectively [72]. MAPK9 (JNK2) has been reported to play a role in p53 regulation and T cell differentiation and activation [73–75]. MAPKs and JNKs have been shown to play a role in the immune response to avian and human IAV infections [76–78]. CD27 is a member of the TNF-receptor family and serves as a T-cell costimulatory molecule [79]. In humans, CD27+ B cells exhibit a memory cell phenotype [80]. In T cells, CD27 is required for efficient CD4+ and CD8+ responses to influenza virus infection. This was only true for generating and maintaining antigen-specific T cell immunity, not effector cell expansion [81]. Of note, these reduced DEPs are involved in lung homeostasis and immunity and their reduced expression negatively correlates with viral load, suggesting that they could potentially be part of a biomarker or signature of infection severity. Importantly, none of these proteins are represented in other standard cytokine-chemokine multiplex panels typically used to analyze infected patient samples.

SOMAscan is capable of measuring a much larger number of cytokines and chemokines than are evaluable by standard techniques (Luminex and ELISA), especially with the limited sample volumes inherent in human translational studies. Additionally, other proteins, not known to contribute to inflammation, nor associated with inflammation or immunity are assessable. Therefore, SOMAscan allowed the identification of potentially novel pathways in the context of viral infections that may not be discovered with a pre-selected or biased multiplex assay such as Luminex and ELISA. While these proteins have reported functions, future studies could identify novel functions during viral infection that would provide additional insights into the complex host-influenza virus interactions.

The SOMAscan results were generally correlated well with Luminex and other assays that identified DEPs in mucosal fluids from influenza virus-infected patients. For example, IL-6, CXCL10/IP10, CCL7/MCP3, and CXCL8/IL-8 were among proteins present at high levels in nasal washes when measured by SOMAscan. Those proteins had similarly high relative levels in the same samples when measured by Luminex [28], but SOMAscan measured 25-fold more proteins. Among the proteins quantified by SOMAscan but not Luminex or other assay methods, many were highly correlated to influenza viral loads. As noted above, all of those proteins are integral to pathways important in inflammation and stress response. DEPs that segregated according to clinical disease severity were not identified. This result is not surprising since clinical disease severity was not correlated with viral loads in the patient samples assayed. Possibly serum protein levels measured by SOMAscan would have yielded correlation with disease severity since many of the clinical symptoms of disease are systemic. Also a larger patient sampling with more samples per disease severity group may yield significant correlations among proteins measured by SOMAscan and disease severity among influenza virus-infected patients. Finally, absolute correlation between different assay systems will always be compromised by differences in the method for collection of results, e.g., light absorbance, fluorescence intensity, etc., and differences in epitope recognition and binding for any given protein. The latter will be particularly relevant for antibody-dependent versus aptamer-dependent ligand binding. As noted in a recent commentary, SOMAscan will be no better or worse for detecting and quantifying a known protein.

The primary benefit provided by SOMAscan is in its discovery potential for identifying associations between proteins among a large catalog of proteins and a phenotype of interest [82]. It is within this context that we think the results from our study can move understanding of respiratory pathophysiology associated with influenza infection forward.

SOMAscan thus represents a very promising alternative methodology to existing technologies for the discovery of biomarkers for disease severity in respiratory infections including influenza. Furthermore, SOMAscan may also have great potential to analyze proteome changes within the mucosa in the context of other bacterial or viral respiratory infections or lung diseases in humans that are presently understudied because of limitations with other multiplex assay technologies.

Materials and Methods

Ethics statement

This study was conducted in compliance with Department of Health and Human Services regulations in 45 CFR46 and the Declaration of Helsinki. The Institutional Review Boards of St. Jude Children’s Research Hospital and the University of Tennessee Health Science Center/Le Bonheur Children’s Hospital approved the study. Nasal lavage samples were acquired from the St, Jude Children's Research Hospital Human Influenza Research Specimen Repository [28]. The samples were originally acquired after receiving written, informed consent from participants or the participants' parents/guardians. Use of the patient samples in the present study was approved by the UTHSC IRB (13-02691-XM).

Patient samples

Subject inclusion criteria, definition of respiratory influenza infection, criteria for categorizing disease severity, and methods for sample collection and storage have been described [28]. Twenty-four samples from the repository were provided for the present study based solely upon the criteria of detection of influenza viral genome and disease severity after infection. Six samples per group were selected from each of four groups: i) uninfected healthy controls: individuals with no clinical signs or symptoms of respiratory infection and otherwise healthy; ii) mild: no emergency room (ER) visit or hospitalization; iii) moderate: ER visit but no hospitalization; and iv) severe: hospitalization. All of the assayed specimens were collected on the day of initial diagnosis (study day 0 in [28]). The uninfected, healthy controls were household contacts of infected patients who did not become infected with influenza. Patient demographics and viral loads determined by qRT-PCR are presented in S1 Table.

SOMAscan analysis

The SOMAscan assay including calibration of the sample matrix for nasal lavage samples was performed by SomaLogic, Boulder, CO. Sixty μl of a 1:2 dilution of each nasal lavage supernatant was used for the SOMAscan assay.

SOMAmers are modified nucleic acid aptamers, each with both unique protein binding characteristics and unique identifying primary nucleic acid sequence that can be detected and quantified by DNA microarray. Each SOMAmer has been validated for its unique specificity, upper and lower limits of detection, and intra- and inter-assay variability. SOMAscan uses a 1,129 SOMAmer library to quantify 1,030 individual proteins or protein complexes that are detected by the library as relative fluorescence units of each SOMAmer bound to its respective protein [23]. Each SOMAmer has a unique specificity for a defined epitope. A given SOMAmer may thus bind to an isolated polypeptide or to a heterodimer of the same protein with another protein. Another SOMAmer may bind to the same heterodimer but neither polypeptide alone. Hence two different SOMAmers may be listed as detecting the same protein, but in fact bind to individual epitopes. Thus, the number of proteins (1,030) that can be detected is less than the total number of SOMAmers (1,129) used in the analysis. Details of SOMAmer validation, the list of proteins quantified by the current SOMAmer library, and SOMAscan assay procedure are at http://www.somalogic.com/Technology/SOMAscan-basic-info.aspx.

Bioinformatics

Data were analyzed using the R software package [83]. Pre-processing steps included quantile normalization and log₂ transformation of the signal intensity (relative light units) measures provided by the SOMAscan platform. Principal component analysis (PCA) analysis were performed using the affycoretools package [84]. Identification of differentially expressed probe sets (DEPs) was performed with the LIMMA package [85] using BH correction for multiple testing [86] using the indicated thresholds for fold-change and adjusted p-values (see results). Heatmaps were performed using the R software package [83]. GO, KEGG, Reactome enrichment analysis and was performed with the R package ReactomePA provided in the R package cluster Profiler [87].

Supporting Information

[Figure omitted. See PDF.]

S1 Fig. Correlation of protein quantification by SOMAscan and Luminex.

Correlation analysis of log₂ relative fluorescence units (RFU) from SOMAscan and protein concentration (log₂ of pg/ml) from Luminex study described by [28] are presented for IL6, CXCL10/IP-10, CCL7/MCP3, and CXCL8/IL8. Labels at top and left refer to individual chemokines and technology used for detection (SO: SOMAscan, LU: Luminex). The diagonal represents the histogram of the measured protein level determined by SOMAscan (SOMA) or Luminex (LUMI), the bottom left below the histogram diagonal represents scatter plots of pairs of measurements, the top right above the histogram diagonal presents correlation coefficients and p-values from a linear regression analysis of the respective pair-wise analyses. n = 24. **: p <0.01, ***: p < 0.001.

https://doi.org/10.1371/journal.pone.0153674.s001

(PDF)

S2 Fig. Distribution of viral loads.

Upper panel: Histogram of log₂-transformed viral loads showing the distribution of viral loads as measured by quantitative real-time reverse-transcription polymerase chain reaction at day 0 from all patients of the cohort [28] (n = 139 samples). The threshold of 256 (log₂ = 8) is shown as a stippled line. Lower panel: Histogram of log₂-transformed viral loads showing the distribution of viral loads as measured by quantitative real-time reverse-transcription polymerase chain reaction of the patients selected for the SOMAscan study (n = 24).

https://doi.org/10.1371/journal.pone.0153674.s002

(PDF)

S3 Fig. PCA analysis of normalized protein expression values for subset B.

Principle component analysis (PCA) was performed with normalized log₂–transformed protein expression values from nasal washes for 17 samples from the subset that were selected based on high vs. low viral loads and omitted outlier sample ID_4043 (subset B). The first two principal components are shown that represent 46% and 16%, respectively, of the total variation. Healthy control samples are labeled gray, and influenza virus-positive samples are labeled red. In addition, sample identities (e.g., ID_4002) are shown. Horizontal and vertical axis represent principle component 1 and 2, respectively.

https://doi.org/10.1371/journal.pone.0153674.s003

(PDF)

S4 Fig. Expression values for CXCL10 and IL6 from subset B.

Normalized log₂–transformed protein concentration of CXCL10 (A) and IL6 (B) for 17 samples from subset B that were selected based on high vs. low viral load. Bars represent mean expression values per group +/- 1 SEM. "negative": samples that were diagnosed IAV-negative; "positive": samples that were diagnosed IAV-positive. IL6 and CXCL10 expression levels showed significant differences between the two groups (Kruskal-Wallis, p < 0.01 and p < 0.05, respectively).

https://doi.org/10.1371/journal.pone.0153674.s004

(PDF)

S1 Table. Demographic data, type of influenza virus identified, log2 viral load, and relative severity of symptoms and clinical disease of respective patients.

https://doi.org/10.1371/journal.pone.0153674.s005

(DOCX)

S2 Table. List of patients included into subsets A and B.

https://doi.org/10.1371/journal.pone.0153674.s006

(DOCX)

S3 Table. List of proteins that were differentially expressed between uninfected control and influenza virus-positive patients from subset A.

The patient identifiers are those from S1 Table. The numerical values for each patient column are the respective normalized log2 protein levels for each patient for each protein identified in column A and B. Log2 Fold Change (FC), Average (average respective protein level for all patients), t-statistic, p, adjusted p, and B were all calculated with LIMMA.

https://doi.org/10.1371/journal.pone.0153674.s007

(TXT)

S4 Table. List of proteins that were differentially expressed between uninfected control and influenza virus-positive patients from subset B.

Column labels are the same as in S2 Table.

https://doi.org/10.1371/journal.pone.0153674.s008

(TXT)

S5 Table. List of DEPS from subset B that were correlated with viral load. cor is Spearman's correlation coefficient (ρ).

https://doi.org/10.1371/journal.pone.0153674.s009

(TXT)

Acknowledgments

The authors wish to acknowledge the study teams at St. Jude Children’s Research Hospital including Trushar Jeevan, and at LeBonheur Children’s Hospital including Lisa Harrison, Shelly Stephens, Ryan Murphy, Magan Stevens, and Lauren Lazar, Philippa Roddam, and the technology support provided by SomaLogic. This work was supported by start-up grants from UTHSC awarded to KS and HE and intramural grants from the Helmholtz-Association (Program Infection and Immunity) awarded to KS. The FLU09 influenza cohort was supported by HHS contracts HHSN266200700005C and HHSN272201400006C, the NIAID-funded St. Jude Center of Excellence for Influenza Research and Surveillance awarded to RW, supports RW and PT. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Author Contributions

Conceived and designed the experiments: TM KS. Performed the experiments: TM KS. Analyzed the data: TM HE JD KS. Contributed reagents/materials/analysis tools: PT JD RW. Wrote the paper: TM HE JD KS.

Citation: Marion T, Elbahesh H, Thomas PG, DeVincenzo JP, Webby R, Schughart K (2016) Respiratory Mucosal Proteome Quantification in Human Influenza Infections. PLoS ONE 11(4): e0153674. https://doi.org/10.1371/journal.pone.0153674

References

1. Fauci AS. Seasonal and pandemic influenza preparedness: science and countermeasures. J Infect Dis. 2006;194 Suppl 2:S73–6. pmid:17163392.

2. Klenk HD, Garten W, Matrosovich M. Molecular mechanisms of interspecies transmission and pathogenicity of influenza viruses: Lessons from the 2009 pandemic. Bioessays. 2011;33(3):180–8. Epub 2011/02/15. pmid:21319184.

3. Kilbourne ED. Influenza pandemics of the 20th century. Emerg Infect Dis. 2006;12(1):9–14. New York Medical College, Valhalla, New York, USA. [email protected] pmid:16494710

4. Russell CJ, Webster RG. The genesis of a pandemic influenza virus. Cell. 2005;123(3):368–71. Department of Infectious Diseases, St.Jude Children's Hospital, Memphis, Tennessee 38105, USA. pmid:16269328

5. Johnson NP, Mueller J. Updating the accounts: global mortality of the 1918–1920 "Spanish" influenza pandemic. Bull Hist Med. 2002;76(1):105–15. pmid:11875246

6. Garten RJ, Davis CT, Russell CA, Shu B, Lindstrom S, Balish A, et al. Antigenic and genetic characteristics of swine-origin 2009 A(H1N1) influenza viruses circulating in humans. Science. 2009;325(5937):197–201. pmid:19465683.

7. Itoh Y, Shinya K, Kiso M, Watanabe T, Sakoda Y, Hatta M, et al. In vitro and in vivo characterization of new swine-origin H1N1 influenza viruses. Nature. 2009;460(7258):1021–5. Epub 2009/08/13. pmid:19672242; PubMed Central PMCID: PMC2748827.

8. Neumann G, Noda T, Kawaoka Y. Emergence and pandemic potential of swine-origin H1N1 influenza virus. Nature. 2009;459(7249):931–9. Epub 2009/06/16. pmid:19525932; PubMed Central PMCID: PMC2873852.

9. Gambotto A, Barratt-Boyes SM, de Jong MD, Neumann G, Kawaoka Y. Human infection with highly pathogenic H5N1 influenza virus. Lancet. 2008;371(9622):1464–75. Department of Surgery, Division of Infectious Diseases, School of Medicine, University of Pittsburgh, Pittsburgh, PA 15219, USA. [email protected] pmid:18440429

10. Krug RM. Virology. Clues to the virulence of H5N1 viruses in humans. Science. 2006;311(5767):1562–3. Institute for Cellular and Molecular Biology, University of Texas at Austin, Austin, TX 78712, USA. [email protected] pmid:16543447

11. Gao R, Cao B, Hu Y, Feng Z, Wang D, Hu W, et al. Human Infection with a Novel Avian-Origin Influenza A (H7N9) Virus. N Engl J Med. 2013. Epub 2013/04/13. pmid:23577628.

12. Kageyama T, Fujisaki S, Takashita E, Xu H, Yamada S, Uchida Y, et al. Genetic analysis of novel avian A(H7N9) influenza viruses isolated from patients in China, February to April 2013. Euro Surveill. 2013;18(15):20453. Epub 2013/04/19. pmid:23594575.

13. Watanabe T, Kiso M, Fukuyama S, Nakajima N, Imai M, Yamada S, et al. Characterization of H7N9 influenza A viruses isolated from humans. Nature. 2013;501:551–5. Epub 2013/07/12. pmid:23842494.

14. Mancuso P. Obesity and respiratory infections: does excess adiposity weigh down host defense? Pulm Pharmacol Ther. 2013;26(4):412–9. Epub 2012/05/29. pmid:22634305; PubMed Central PMCID: PMC3439591.

15. Karlsson EA, Marcelin G, Webby RJ, Schultz-Cherry S. Review on the impact of pregnancy and obesity on influenza virus infection. Influenza Other Respir Viruses. 2012;6(6):449–60. Epub 2012/02/18. pmid:22335790.

16. Albright FS, Orlando P, Pavia AT, Jackson GG, Cannon Albright LA. Evidence for a heritable predisposition to death due to influenza. J Infect Dis. 2008;197(1):18–24. Department of Pharmacotherapy, University of Utah College of Pharmacy, Salt Lake City, USA. pmid:18171280

17. Gottfredsson M, Halldorsson BV, Jonsson S, Kristjansson M, Kristjansson K, Kristinsson KG, et al. Lessons from the past: familial aggregation analysis of fatal pandemic influenza (Spanish flu) in Iceland in 1918. Proc Natl Acad Sci U S A. 2008;105(4):1303–8. Landspitali University Hospital, 108 Reykjavik, Iceland. [email protected] pmid:18216264

18. Horby P, Sudoyo H, Viprakasit V, Fox A, Thai PQ, Yu H, et al. What is the evidence of a role for host genetics in susceptibility to influenza A/H5N1? Epidemiol Infect. 2010;138(11):1–9. Epub 2010/03/20. S0950268810000518 [pii] pmid:20236573.

19. O'Brien KB, Vogel P, Duan S, Govorkova EA, Webby RJ, McCullers JA, et al. Impaired wound healing predisposes obese mice to severe influenza virus infection. J Infect Dis. 2012;205(2):252–61. Epub 2011/12/08. pmid:22147799; PubMed Central PMCID: PMC3244366.

20. Tisoncik JR, Korth MJ, Simmons CP, Farrar J, Martin TR, Katze MG. Into the eye of the cytokine storm. Microbiol Mol Biol Rev. 2012;76(1):16–32. Epub 2012/03/07. pmid:22390970; PubMed Central PMCID: PMC3294426.

21. Bagga B, Woods CW, Veldman TH, Gilbert A, Mann A, Balaratnam G, et al. Comparing influenza and RSV viral and disease dynamics in experimentally infected adults predicts clinical effectiveness of RSV antivirals. Antivir Ther. 2013;18(6):785–91. Epub 2013/05/30. pmid:23714753.

22. Kraemer S, Vaught JD, Bock C, Gold L, Katilius E, Keeney TR, et al. From SOMAmer-based biomarker discovery to diagnostic and clinical applications: a SOMAmer-based, streamlined multiplex proteomic assay. PLoS ONE. 2011;6(10):e26332. Epub 2011/10/25. pmid:22022604; PubMed Central PMCID: PMC3195687.

23. Rohloff JC, Gelinas AD, Jarvis TC, Ochsner UA, Schneider DJ, Gold L, et al. Nucleic Acid Ligands With Protein-like Side Chains: Modified Aptamers and Their Use as Diagnostic and Therapeutic Agents. Molecular therapy Nucleic acids. 2014;3:e201. Epub 2014/10/08. pmid:25291143; PubMed Central PMCID: PMCPmc4217074.

24. Gold L, Ayers D, Bertino J, Bock C, Bock A, Brody EN, et al. Aptamer-based multiplexed proteomic technology for biomarker discovery. PLoS ONE. 2010;5(12):e15004. Epub 2010/12/18. pmid:21165148; PubMed Central PMCID: PMC3000457.

25. Loffredo FS, Steinhauser ML, Jay SM, Gannon J, Pancoast JR, Yalamanchi P, et al. Growth differentiation factor 11 is a circulating factor that reverses age-related cardiac hypertrophy. Cell. 2013;153(4):828–39. Epub 2013/05/15. pmid:23663781; PubMed Central PMCID: PMC3677132.

26. De Groote MA, Nahid P, Jarlsberg L, Johnson JL, Weiner M, Muzanyi G, et al. Elucidating novel serum biomarkers associated with pulmonary tuberculosis treatment. PLoS ONE. 2013;8(4):e61002. Epub 2013/05/03. pmid:23637781; PubMed Central PMCID: PMC3630118.

27. Mehan MR, Ayers D, Thirstrup D, Xiong W, Ostroff RM, Brody EN, et al. Protein signature of lung cancer tissues. PLoS ONE. 2012;7(4):e35157. Epub 2012/04/18. pmid:22509397; PubMed Central PMCID: PMC3324437.

28. Oshansky CM, Gartland AJ, Wong SS, Jeevan T, Wang D, Roddam PL, et al. Mucosal Immune Responses Predict Clinical Outcomes during Influenza Infection Independently of Age and Viral Load. Am J Respir Crit Care Med. 2014;189:449–62. Epub 2013/12/07. pmid:24308446.

29. Berdal JE, Mollnes TE, Waehre T, Olstad OK, Halvorsen B, Ueland T, et al. Excessive innate immune response and mutant D222G/N in severe A (H1N1) pandemic influenza. J Infect. 2011;63(4):308–16. Epub 2011/07/26. pmid:21781987.

30. Kollmus H, Wilk E, Schughart K. Systems biology and systems genetics-novel innovative approaches to study host-pathogen interactions during influenza infection. Curr Opin Virol. 2014;6C:47–54. Epub 2014/04/29. pmid:24769047.

31. Cole KE, Strick CA, Paradis TJ, Ogborne KT, Loetscher M, Gladue RP, et al. Interferon-inducible T cell alpha chemoattractant (I-TAC): a novel non-ELR CXC chemokine with potent activity on activated T cells through selective high affinity binding to CXCR3. J Exp Med. 1998;187(12):2009–21. Epub 1998/06/24. pmid:9625760; PubMed Central PMCID: PMCPmc2212354.

32. Qin S, Rottman JB, Myers P, Kassam N, Weinblatt M, Loetscher M, et al. The chemokine receptors CXCR3 and CCR5 mark subsets of T cells associated with certain inflammatory reactions. J Clin Invest. 1998;101(4):746–54. Epub 1998/03/21. pmid:9466968; PubMed Central PMCID: PMCPmc508621.

33. Sallusto F, Lenig D, Mackay CR, Lanzavecchia A. Flexible programs of chemokine receptor expression on human polarized T helper 1 and 2 lymphocytes. J Exp Med. 1998;187(6):875–83. Epub 1998/04/04. pmid:9500790; PubMed Central PMCID: PMCPmc2212187.

34. Cruikshank WW, Long A, Tarpy RE, Kornfeld H, Carroll MP, Teran L, et al. Early identification of interleukin-16 (lymphocyte chemoattractant factor) and macrophage inflammatory protein 1 alpha (MIP1 alpha) in bronchoalveolar lavage fluid of antigen-challenged asthmatics. Am J Respir Cell Mol Biol. 1995;13(6):738–47. Epub 1995/12/01. pmid:7576712.

35. Cruikshank WW, Center DM, Nisar N, Wu M, Natke B, Theodore AC, et al. Molecular and functional analysis of a lymphocyte chemoattractant factor: association of biologic function with CD4 expression. Proc Natl Acad Sci U S A. 1994;91(11):5109–13. Epub 1994/05/24. pmid:7910967; PubMed Central PMCID: PMCPmc43941.

36. Deng A, Chen S, Li Q, Lyu SC, Clayberger C, Krensky AM. Granulysin, a cytolytic molecule, is also a chemoattractant and proinflammatory activator. J Immunol. 2005;174(9):5243–8. Epub 2005/04/22. pmid:15843520.

37. Shoyab M, Plowman GD, McDonald VL, Bradley JG, Todaro GJ. Structure and function of human amphiregulin: a member of the epidermal growth factor family. Science. 1989;243(4894 Pt 1):1074–6. Epub 1989/02/24. pmid:2466334.

38. Berasain C, Avila MA. Amphiregulin. Seminars in cell & developmental biology. 2014;28:31–41. Epub 2014/01/28. pmid:24463227.

39. Arpaia N, Green JA, Moltedo B, Arvey A, Hemmers S, Yuan S, et al. A Distinct Function of Regulatory T Cells in Tissue Protection. Cell. 2015;162(5):1078–89. Epub 2015/09/01. pmid:26317471; PubMed Central PMCID: PMCPmc4603556.

40. Jamieson AM, Pasman L, Yu S, Gamradt P, Homer RJ, Decker T, et al. Role of tissue protection in lethal respiratory viral-bacterial coinfection. Science. 2013;340(6137):1230–4. Epub 2013/04/27. pmid:23618765; PubMed Central PMCID: PMCPmc3933032.

41. Monticelli LA, Sonnenberg GF, Abt MC, Alenghat T, Ziegler CG, Doering TA, et al. Innate lymphoid cells promote lung-tissue homeostasis after infection with influenza virus. Nat Immunol. 2011;12(11):1045–54. Epub 2011/09/29. pmid:21946417; PubMed Central PMCID: PMCPmc3320042.

42. Fujita J, Kawaguchi M, Kokubu F, Ohara G, Ota K, Huang SK, et al. Interleukin-33 induces interleukin-17F in bronchial epithelial cells. Allergy. 2012;67(6):744–50. Epub 2012/05/01. pmid:22540331.

43. Lin CH, Shih CH, Chen BC. Thrombin-induced IL-8/CXCL8 release is mediated by CK2, MSK1, and NF-kappaB pathways in human lung epithelial cells. European journal of pharmacology. 2015;767:135–43. Epub 2015/10/16. pmid:26463037.

44. Moens U, Kostenko S, Sveinbjornsson B. The Role of Mitogen-Activated Protein Kinase-Activated Protein Kinases (MAPKAPKs) in Inflammation. Genes. 2013;4(2):101–33. Epub 2013/01/01. pmid:24705157; PubMed Central PMCID: PMCPmc3899974.

45. Vettorazzi S, Bode C, Dejager L, Frappart L, Shelest E, Klassen C, et al. Glucocorticoids limit acute lung inflammation in concert with inflammatory stimuli by induction of SphK1. Nat Commun. 2015;6:7796. Epub 2015/07/18. pmid:26183376; PubMed Central PMCID: PMCPmc4518295.

46. Joo JH, Jetten AM. NF-kappaB-dependent transcriptional activation in lung carcinoma cells by farnesol involves p65/RelA(Ser276) phosphorylation via the MEK-MSK1 signaling pathway. J Biol Chem. 2008;283(24):16391–9. Epub 2008/04/22. pmid:18424438; PubMed Central PMCID: PMCPmc2423266.

47. Brasier AR, Tian B, Jamaluddin M, Kalita MK, Garofalo RP, Lu M. RelA Ser276 phosphorylation-coupled Lys310 acetylation controls transcriptional elongation of inflammatory cytokines in respiratory syncytial virus infection. J Virol. 2011;85(22):11752–69. Epub 2011/09/09. pmid:21900162; PubMed Central PMCID: PMCPmc3209292.

48. Wu Y. Contact pathway of coagulation and inflammation. Thrombosis journal. 2015;13:17. Epub 2015/05/08. pmid:25949215; PubMed Central PMCID: PMCPmc4421925.

49. Govorkova EA, Baranovich T, Marathe BM, Yang L, Taylor MA, Webster RG, et al. Sialic Acid-Binding Protein Sp2CBMTD Protects Mice Against Lethal Challenge with the Emerging A(H7N9) Influenza Virus. Antimicrob Agents Chemother. 2014;59:1495–504. Epub 2014/12/24. pmid:25534734.

50. Doyle WJ, Skoner DP, White M, Hayden F, Kaplan AP, Kaliner MA, et al. Pattern of nasal secretions during experimental influenza virus infection. Rhinology. 1996;34(1):2–8. Epub 1996/03/01. pmid:8739859.

51. Headland SE, Norling LV. The resolution of inflammation: Principles and challenges. Semin Immunol. 2015;27(3):149–60. Epub 2015/04/26. pmid:25911383.

52. Oshansky CM, Thomas PG. The human side of influenza. J Leukoc Biol. 2012;92(1):83–96. Epub 2012/03/01. pmid:22362872; PubMed Central PMCID: PMC3382310.

53. Iwasaki A, Pillai PS. Innate immunity to influenza virus infection. Nat Rev Immunol. 2014;14(5):315–28. Epub 2014/04/26. pmid:24762827.

54. Dutta A, Huang CT, Chen TC, Lin CY, Chiu CH, Lin YC, et al. IL-10 inhibits neuraminidase-activated TGF-beta and facilitates Th1 phenotype during early phase of infection. Nat Commun. 2015;6:6374. Epub 2015/03/03. pmid:25728041.

55. Chi Y, Zhu Y, Wen T, Cui L, Ge Y, Jiao Y, et al. Cytokine and chemokine levels in patients infected with the novel avian influenza A (H7N9) virus in China. J Infect Dis. 2013;208(12):1962–7. Epub 2013/08/31. pmid:23990573.

56. Crowley JJ, Zhabotynsky V, Sun W, Huang S, Pakatci IK, Kim Y, et al. Analyses of allele-specific gene expression in highly divergent mouse crosses identifies pervasive allelic imbalance. Nat Genet. 2015;47(4):353–60. Epub 2015/03/03. pmid:25730764; PubMed Central PMCID: PMCPmc4380817.

57. Liu C, Zhang A, Guo J, Yang J, Zhou H, Chen H, et al. Identification of human host proteins contributing to H5N1 influenza virus propagation by membrane proteomics. J Proteome Res. 2012;11(11):5396–405. pmid:22985253.

58. Liu L, Zhou J, Wang Y, Mason RJ, Funk CJ, Du Y. Proteome alterations in primary human alveolar macrophages in response to influenza A virus infection. J Proteome Res. 2012;11(8):4091–101. pmid:22709384; PubMed Central PMCID: PMC3412919.

59. Kroeker AL, Ezzati P, Halayko AJ, Coombs KM. Response of primary human airway epithelial cells to influenza infection: a quantitative proteomic study. J Proteome Res. 2012;11(8):4132–46. Epub 2012/06/15. pmid:22694362; PubMed Central PMCID: PMC3411195.

60. Cheung CY, Chan EY, Krasnoselsky A, Purdy D, Navare AT, Bryan JT, et al. H5N1 virus causes significant perturbations in host proteome very early in influenza virus-infected primary human monocyte-derived macrophages. J Infect Dis. 2012;206(5):640–5. Epub 2012/07/24. pmid:22822004; PubMed Central PMCID: PMC3491741.

61. Lietzen N, Ohman T, Rintahaka J, Julkunen I, Aittokallio T, Matikainen S, et al. Quantitative subcellular proteome and secretome profiling of influenza A virus-infected human primary macrophages. PLoS Pathog. 2011;7(5):e1001340. Epub 2011/05/19. pmid:21589892; PubMed Central PMCID: PMC3093355.

62. Coombs KM, Berard A, Xu W, Krokhin O, Meng X, Cortens JP, et al. Quantitative proteomic analyses of influenza virus-infected cultured human lung cells. J Virol. 2010;84(20):10888–906. Epub 2010/08/13. pmid:20702633; PubMed Central PMCID: PMC2950599.

63. Vester D, Rapp E, Gade D, Genzel Y, Reichl U. Quantitative analysis of cellular proteome alterations in human influenza A virus-infected mammalian cell lines. Proteomics. 2009;9(12):3316–27. Epub 2009/06/09. pmid:19504497.

64. Liu N, Song W, Wang P, Lee K, Chan W, Chen H, et al. Proteomics analysis of differential expression of cellular proteins in response to avian H9N2 virus infection in human cells. Proteomics. 2008;8(9):1851–8. Epub 2008/04/10. pmid:18398875.

65. Qi W, Tian J, Su S, Huang L, Li H, Liao M. Identification of Potential Virulence Determinants associated H9N2 Avian Influenza Virus PB2 E627K Mutation by Comparative Proteomics. Proteomics. 2015;15:1512–24. Epub 2015/02/03. pmid:25641917.

66. Choi S, Lim JY, Kim Y, Song MJ, Jung WW, Seo JB, et al. Plasma proteomic analysis of patients infected with H1N1 influenza virus. Proteomics. 2014;14(16):1933–42. Epub 2014/06/04. pmid:24888898.

67. Kovner AV, Potapova OV, Shkurupy VA. Morphofunctional changes in the lung vascular endotheliocytes in mice with influenza A/H5N1 A/Goose/Krasnoozerskoye/627/05. Bull Exp Biol Med. 2013;154(4):476–9. pmid:23486585.

68. Potapova OV, Sharkova TV, Shkurupiy VA, Shestopalov AM. Immunomorphologic manifestations in mice liver infected with influenza A/H5N1, A/goose/Krasnoozerskoye/627/05 strain. Clin Dev Immunol. 2013;2013:342686. pmid:24454472; PubMed Central PMCID: PMCPMC3886489.

69. Skytt A, Stromqvist M, Egelrud T. Primary substrate specificity of recombinant human stratum corneum chymotryptic enzyme. Biochem Biophys Res Commun. 1995;211(2):586–9. pmid:7794273.

70. Petraki CD, Papanastasiou PA, Karavana VN, Diamandis EP. Cellular distribution of human tissue kallikreins: immunohistochemical localization. Biol Chem. 2006;387(6):653–63. pmid:16800726.

71. Larsson A, Peng S, Persson H, Rosenbloom J, Abrams WR, Wassberg E, et al. Lactadherin binds to elastin—a starting point for medin amyloid formation? Amyloid. 2006;13(2):78–85. pmid:16911961.

72. Hanayama R, Tanaka M, Miwa K, Shinohara A, Iwamatsu A, Nagata S. Identification of a factor that links apoptotic cells to phagocytes. Nature. 2002;417(6885):182–7. pmid:12000961.

73. Hu MC, Qiu WR, Wang YP. JNK1, JNK2 and JNK3 are p53 N-terminal serine 34 kinases. Oncogene. 1997;15(19):2277–87. pmid:9393873.

74. Spurlock CF 3rd, Tossberg JT, Fuchs HA, Olsen NJ, Aune TM. Methotrexate increases expression of cell cycle checkpoint genes via JNK activation. Arthritis Rheum. 2012;64(6):1780–9. pmid:22183962; PubMed Central PMCID: PMCPMC3310965.

75. Yang DD, Conze D, Whitmarsh AJ, Barrett T, Davis RJ, Rincon M, et al. Differentiation of CD4+ T cells to Th1 cells requires MAP kinase JNK2. Immunity. 1998;9(4):575–85. pmid:9806643.

76. Hui KP, Lee SM, Cheung CY, Ng IH, Poon LL, Guan Y, et al. Induction of proinflammatory cytokines in primary human macrophages by influenza A virus (H5N1) is selectively regulated by IFN regulatory factor 3 and p38 MAPK. J Immunol. 2009;182(2):1088–98. pmid:19124752.

77. Nacken W, Ehrhardt C, Ludwig S. Small molecule inhibitors of the c-Jun N-terminal kinase (JNK) possess antiviral activity against highly pathogenic avian and human pandemic influenza A viruses. Biol Chem. 2012;393(6):525–34. Epub 2012/05/26. pmid:22628315.

78. Conze D, Krahl T, Kennedy N, Weiss L, Lumsden J, Hess P, et al. c-Jun NH(2)-terminal kinase (JNK)1 and JNK2 have distinct roles in CD8(+) T cell activation. J Exp Med. 2002;195(7):811–23. pmid:11927626; PubMed Central PMCID: PMCPMC2193724.

79. Gravestein LA, Nieland JD, Kruisbeek AM, Borst J. Novel mAbs reveal potent co-stimulatory activity of murine CD27. Int Immunol. 1995;7(4):551–7. pmid:7547681.

80. Klein U, Rajewsky K, Kuppers R. Human immunoglobulin (Ig)M+IgD+ peripheral blood B cells expressing the CD27 cell surface antigen carry somatically mutated variable region genes: CD27 as a general marker for somatically mutated (memory) B cells. J Exp Med. 1998;188(9):1679–89. pmid:9802980; PubMed Central PMCID: PMCPMC2212515.

81. Hendriks J, Gravestein LA, Tesselaar K, van Lier RA, Schumacher TN, Borst J. CD27 is required for generation and long-term maintenance of T cell immunity. Nat Immunol. 2000;1(5):433–40. pmid:11062504.

82. Powell K. New platform for cataloging hundreds of proteins gets test drive. Nat Med. 2014;20(10):1082–3. Epub 2014/10/09. pmid:25295926.

83. R_Core_Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. 2013. Available: http://www.R-project.org/.

84. MacDonald J. affycoretools: Functions useful for those doing repetitive analyses with Affymetrix GeneChips. R package version 1340. 2008.

85. Smyth GK. Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004;3:Article3. Epub 2006/05/02. pmid:16646809.

86. Benjamini Y, Hochberg Y. Controlling the false discovery rate: A practical and powerful apporach to multiple testing. Journal of the Royal Statistical Society. 1995;57:289–300.

87. Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–7. Epub 2012/03/30. pmid:22455463; PubMed Central PMCID: PMC3339379.