Allergy, Asthma & Clinical Immunology
Whole blood vs PBMC: compartmental differences in gene expression profiling exemplified in asthma
Daniel He1,2,
Chen Xi Yang1,2,
Basak Sahin1,
Amrit Singh1,2,
Casey P. Shannon2,
John-Paul Oliveria3,4,
Gail M. Gauvreau3 &
Scott J. Tebbutt ORCID: orcid.org/0000-0002-7908-15811,2,5
Allergy, Asthma & Clinical Immunology volume 15, Article number: 67 (2019) Cite this article
Blood has proven to be a useful resource for molecular analysis in numerous biomedical studies, with peripheral blood mononuclear cells (PBMCs) and whole blood being the major specimen types. However, comparative analyses between these two major compartments (PBMCs and whole blood) are few and far between. In this study, we compared gene expression profiles of PBMCs and whole blood samples obtained from research subjects with or without mild allergic asthma.
Whole blood (PAXgene) and PBMC samples were obtained from 5 mild allergic asthmatics and 5 healthy controls. RNA from both sample types was measured for expression of 730 immune-related genes using the NanoString nCounter platform.
We identified 64 uniquely expressed transcripts in whole blood that reflected a variety of innate, humoral, and adaptive immune processes, and 13 uniquely expressed transcripts in PBMCs which were representative of T-cell and monocyte-mediated processes. Furthermore, analysis of mild allergic asthmatics versus non-asthmatics revealed 47 differentially expressed transcripts in whole blood compared to 1 differentially expressed transcript in PBMCs (FDR < 0.25). Finally, through simultaneous measurement of PBMC proteins on the nCounter assay, we identified CD28 and OX40 (TNFRSF4), both of which are critical co-stimulatory molecules during T-cell activation, as significantly upregulated in asthmatics.
Whole blood RNA preserved in PAXgene tubes is excellent for producing gene expression data with minimal variability and good sensitivity, suggesting its utility in multi-centre studies requiring measurement of blood gene expression.
Gene expression profiling is a fundamental feature of systems biology studies, which are ever-increasing in use to unravel the complexity of biological systems and their disease states. In particular, measurement of gene expression ('transcriptomics') in blood is an effective tool for the discovery of biomarkers of disease, and remains a key component of studies requiring large numbers of samples due to ease of collection and storage. Although the advent of single-cell RNA sequencing has made it possible to measure gene expression of specific blood cell types, there remain financial and technical barriers to its widespread use. Hence, the most common methods of blood RNA profiling either investigate RNA extracted from all blood cell types (i.e. whole blood) or only peripheral blood mononuclear cells (PBMCs). While PBMCs (lymphocytes and monocytes) play a significant role in the immune system, they remain a subset of all immune cells and do not include other cell types such as eosinophils and neutrophils. Thus, consideration of the compartment being investigated is required when carrying out blood gene expression studies, particularly with respect to the disease or other condition of interest.
Asthma is a chronic disease of the lungs that is characterized by frequent dyspnea, wheezing, coughing, and chest tightness due to inflammation of the airways. Asthma affects an estimated 334 million people worldwide and its prevalence continues to rise [1]. This is accompanied by a high economic burden, with recent estimates suggesting an annual cost (both from medical costs and loss of productivity) of €72.2 billion in the European Union [2] and $56 billion in the United States [3]. Asthma is a heterogeneous disorder with numerous immunological pathways contributing to pathogenesis. Historically, allergic asthma has been described as a T helper type 2 (TH2) cell-mediated inflammatory response, characterized by airway and blood eosinophilia. Further investigations have identified a subset of patients with a predominantly neutrophilic asthma characterized by involvement of T helper type 17 (TH17) cells. While asthma can be controlled with corticosteroid treatment, some patients are steroid-resistant and require specific targeted therapies to manage their symptoms. This is more common in those with TH17-predominant asthma [4].
The use of transcriptomics has been a key method in uncovering the phenotypes of asthma. Microarray analysis of airway epithelial brushings revealed two distinct clusters of TH2-high and TH2-low asthmatic patients, the latter of which displayed no significant difference in expression of TH2 cytokines (IL-4, IL-5, IL-13) when compared to healthy controls [5]. Although there is a systemic component to the pathophysiology of asthma, differentially expressed genes (DEGs) found in the blood of asthmatic individuals may not translate to those found in lung. Results from a recent U-BIOPRED study showed that DEGs associated with eosinophilic TH2 asthma in the lung were not found in whole blood; in fact, only one of the transcripts was detectable (SERPINB2) and not significantly different between asthmatics and healthy controls [6]. Similarly, investigations of PBMCs from asthmatic individuals also displayed a unique gene expression profile differing from that found in asthmatic lung as well as whole blood [7, 8]. Thus, in order to determine blood biomarkers of asthma, it is evident that the merits and limitations of investigating either blood compartment require thorough evaluation. In this study, we sought to compare the gene expression profiles of whole blood and PBMC compartments in a cohort of mild allergic asthmatic individuals and healthy controls.
Whole blood (in PAXgene RNA tubes, PreAnalytiX, Hombrechtikon, Switzerland) and PBMCs (Ficoll preparation) were obtained from 8 mild allergic asthmatics defined as skin test positive to common aeroallergens and having FEV1 (forced expiratory volume in 1 s) ≤ 80% of predicted and methacholine PC20 (provocative concentration causing 20% fall in FEV1) ≤ 16 mg/ml, and 8 healthy controls (skin test negative, FEV1 ˃ 80% of predicted, methacholine PC20 ˃ 16 mg/ml) at McMaster University (Hamilton, ON, Canada) (Table 1). Prior to methacholine challenge, subjects were instructed to avoid the use of asthma medication prior to testing (for further details, see Ref. [9]). RNA was extracted from whole blood samples using PAXgene Blood miRNA Kit (PreAnalytiX, Hombrechtikon, Switzerland), then quantified (100 ng) for expression of 730 transcripts using the NanoString nCounter panCancer Immune Profiling panel (NanoString Technologies, Seattle, WA). PBMCs were lysed with RLT buffer (Qiagen, Hilden Germany) and analyzed for gene transcript and protein expression using the nCounter panCancer Immune Profiling plus Vantage 3D Protein Immune Cell Profiling panel (NanoString Technologies, Seattle, WA). In order to directly compare subject-specific pairs of whole blood and PBMCs we excluded data from 6 subjects (3 allergic asthmatic, 3 healthy controls) due to poor data quality of the PBMC sample as determined by NanoString quality control standard operating procedures (SOPs). Data were normalized using total sum scaling, and transcripts were filtered for low abundance using the geometric mean of negative controls (no target RNA). Normalization of protein data was also performed using total sum scaling, with an additional normalization against IgG negative controls. Differential expression testing between whole blood and PBMC fractions, as well as asthmatic and control groups, was performed using Linear Models for Microarrays and RNA-Seq (LIMMA, version 3.38.3) in the R statistical environment (version 3.5.2). Where appropriate, subjects were included into our linear modeling as random effects. False positives introduced by multiple testing were controlled for using the Benjamini–Hochberg false discovery rate, FDR (LIMMA default parameter). Gene expression data has been deposited in the Gene Expression Omnibus under accession number GSE132006, and normalized protein values can be found in Additional file 1.
Table 1 Characteristics of subjects analysed for both whole blood and PBMC gene expression
To examine overall gene expression, analysis of all measured transcripts on the NanoString Immune Profiling Panel (730 genes in total) revealed 704 differentially expressed genes (FDR < 0.25) between whole blood and PBMC compartments, of which only 6 genes (EGR1, IL32, FOS, CCL3L1, IFNL1, and EGR2) had increased levels of expression in PBMC samples compared to whole blood samples (Table 2; Fig. 1a, b). Pearson correlation analysis revealed an average correlation of 0.97 between whole blood samples and 0.90 between PBMC samples. Between whole blood and PBMC samples taken from the same subject, we calculated an average correlation of 0.87 (Table 3, Fig. 1c).
Table 2 Top 6 upregulated differentially expressed genes between PBMC and whole blood ordered by log2 fold change (log2FC)
Comparison of gene expression profiles of PBMC and whole blood samples. Top 6 differentially expressed genes with increased expression (FDR < 0.25) in a PBMC samples compared to whole blood samples and b whole blood samples compared to PBMC samples. Y-axis values reflect log2-transformed gene expression. c Correlation matrix of gene expression measured in PBMC and whole blood samples
Table 3 Pearson correlation of RNA expression values between whole blood and PBMCs within subjects
In order to examine compartment-specific pathways, we filtered for low abundant transcripts (all measured counts of a gene transcript falling below a background threshold set by the geometric mean of negative controls) and found 64 transcripts detected above background in whole blood but not PBMC samples, and 13 transcripts detected above background in PBMC but not whole blood samples (Fig. 2a). Notably, GATA3 gene expression was detectable only in the PBMC samples along with IL4, both of which are mediators of TH2-immunity, in addition to other unique PBMC transcripts such as CXCL2, CXCL3, EGR1, and EGR2. We next analyzed tissue-specific genes using NanoString-curated transcript annotations. Unsurprisingly, the top pathways associated with unique transcripts in whole blood reflected a variety of immune processes such as the adaptive, innate, and humoral immune responses (Fig. 2b). For the PBMC-specific genes, the annotated pathways reflected T-lymphocyte and monocyte/macrophage-specific processes (Fig. 2c).
Detectable transcripts in whole blood and PBMC after removing lowly abundant transcripts. a Uniquely detected transcripts in whole blood versus PBMC. Intersection plots (generating using R package 'UpSetR' version 1.3.3) showing pathways associated with b unique whole blood transcripts and c unique PBMC transcripts. On the left are the number of uniquely detected whole blood (WB) or PBMC genes identified in each pathway. The matrix identifies gene intersections and combinations of intersections between pathways, with the top bar plot representing the number of genes within each association. Both intersection plots have been truncated to show the top 8 pathways by number of intersecting genes and a maximum of 10 between-pathway intersections
In addition to being uniquely detected in PBMC samples, EGR2 was the only differentially expressed transcript (FDR < 0.25) when comparing asthmatics to healthy controls (see blue point in Fig. 3a, Additional file 2). Of note, GATA3 was differentially expressed at a nominal p-value of 0.0028, but did not survive FDR cut-off. Compared to controls, we identified 41 downregulated and 6 upregulated differentially expressed transcripts (FDR < 0.25) in asthmatic whole blood. Pathway analysis of the downregulated transcripts reflected general immune processes in addition to specific functions relating to NK cells and TH1-mediated immunity (Fig. 3b).
Differentially expressed genes in allergic asthma. a Volcano plot of differentially expressed genes (FDR < 0.25) in allergic asthmatics from whole blood and PBMC. b UpSet plot of downregulated genes in allergic asthma whole blood. Top 10 pathways by number of intersecting genes are shown, with a maximum of 20 between-pathway intersections
The NanoString nCounter Vantage 3D Immune Cell Profiling assay allows for simultaneous measurement of RNA and protein expression. Of the 30 proteins on the panel, 10 were identified as positively correlated (p < 0.05, Pearson) with RNA expression levels measured in PBMCs (Additional file 3, Fig. 4a). NCR1 (natural cytotoxicity triggering receptor 1) protein expression was found to be significantly negatively correlated (p < 0.05, Pearson) with whole blood gene expression of NCR1 (Additional file 3, Fig. 4a). Differential expression analysis of measured proteins in all PBMC samples (8 asthmatic, 7 controls) revealed upregulated OX40 (TNFRSF4) and CD28 in the asthmatic cohort (Table 4, Fig. 4b); neither of these were differentially expressed in PBMC-extracted RNA, but the gene encoding for the OX40 ligand (TNFSF4) was upregulated in asthmatic whole blood-extracted RNA (Fig. 3a).
Analysis of protein expression in PBMCs. a Significant correlations (p < 0.05) of gene and protein expression in PBMCs (blue) and whole blood (orange), as measured using Pearson correlation. b Scatterplots of the two differentially expressed proteins in the allergic asthmatic cohort and their corresponding RNA scatterplots
Table 4 Differentially expressed proteins between allergic asthmatics (n = 8) and healthy controls (n = 7)
Unraveling the complexity of human disease has been enhanced with the rapid development of high-throughput molecular profiling. The examination of immune cells within blood samples is a commonly used method to elucidate disease mechanisms [10]. Here, we compared two types of blood sampling commonly used in immunological studies: PBMCs and whole blood.
Our results show that RNA extracted from whole blood stored in PAXgene tubes yielded higher counts of the 730 measured transcripts compared to RNA from PBMC lysates. Only 6 transcripts were expressed in greater quantities in PBMC lysates compared to whole blood, whereas whole blood had 704 transcripts with higher expression. A reason for this may be the lack of neutrophils, basophils, and eosinophils in PBMC samples. Basophils and eosinophils are only a small subset of all immune cells (0–2% and 1–7%, respectively). Neutrophils make up a majority of circulating nucleated blood cells (45–75%), but a relatively low RNA content relative to normal cells [11]. While this may be partially mitigated through overrepresentation of genes important in immune functions within our 730 measured transcripts (almost half of the 64 uniquely expressed transcripts we identified in whole blood are expressed by neutrophils), the strong correlation observed between PBMC and whole blood samples from the same subject (mean 0.87) suggests that granulocytic involvement is unlikely to explain the higher signal obtained from whole blood. We also cannot rule out technical differences in sample processing. Namely, there may have been disparities in the quality of total RNA extracted from the PAXgene RNA preservation reagent compared to the RLT buffer used for PBMC lysates. Within sample types, whole blood samples (r = 0.97) appear to be more robust than PBMC samples (r = 0.90), though both are acceptable as high positive correlations [12].
The data generated from our study present a contrast to previous studies comparing gene expression profiles obtained from PBMCs and whole blood. Feezor et al., Debey et al., Palmer et al., Bondar et al., Min et al. [13,14,15,16,17] found that RNA extracted from PBMCs had a higher abundance of gene expression compared to whole blood (PAXgene) as measured through microarray technology. This is likely due to the abundance of globin transcripts present in whole blood (80–90%), which causes decreased sensitivity in detection of other transcripts [14]. Hemoglobin was identified in a majority of these studies as being highly expressed in the whole blood gene expression profiles, and four of the five studies acknowledged globin as a potential cause of low signal-to-noise ratios in whole blood. Our targeted measurement approach to gene expression (730 genes), though more limited in scope, mitigates these effects as globin transcripts were not measured on the NanoString panCancer Immune Profiling Panel. Despite this, correlations between expression profiles of PBMC and whole blood were found to be between 0.78 and 0.91, which are in line with our results [13, 16].
Interestingly, another group used a similar approach to the aforementioned studies, yet found that whole blood (PAXgene) transcripts produced a greater signal compared to PBMC transcripts despite globin (HBB) being highly expressed [18]. A reason for this may be the use of the Affymetrix Human Exon (HuEx) 1.0 ST microarray, which measures gene expression through exon probing and is unique when compared to the aforementioned studies which used the U133A (Affymetrix, Santa Clara, CA), Lymphochip, HumanHT-12 v4 (Illumina Inc, San Diego, CA), and Sentrix Human-6 v2 (Illumina) microarrays. While the HuEx technology is able to discern alternatively spliced transcripts and thus gene isoforms, it has worse detection rates compared to conventional microarrays and has poorer reproducibility for genes with fewer exons [19]. With this in mind, it is possible that the effects of globin were dampened as the globin genes HBA1, HBA2, HBB, and HBD each contain 3 exons, which is low when compared to the mean number of exons per gene (8.8) [20].
The upregulated transcripts in PBMC samples (EGR1, IL32, FOS, CCL3L1, IFNL1, and EGR2) have distinct roles in lymphocytes and monocytes. FOS, which encodes for the proto-oncogene c-Fos, acts as a suppressor of the immune system [21]. IL-32 is primarily expressed in lymphocytes, monocytes, and natural killer cells, and induces production of other pro-inflammatory cytokines such as IL-8 and TNF-α [22]. IFNL1 encodes for interferon lambda (IFN-λ, also known as IL-29), which is produced in response to viral infections by myeloid cells [23]. CCL3L1 is produced by PBMCs as a monocyte chemoattractant [24].
The early growth response transcription factors EGR1 and EGR2 play opposing roles in immunity; while EGR1 positively regulates T and B lymphocyte activation upon antigenic stimulation, EGR2 is involved in T cell anergy and is involved in the FasL-mediated apoptotic pathway [25]. In our analysis, we found that both EGR1 and EGR2 were only detected in PBMC samples, and EGR2 was differentially expressed in our asthmatic cohort. EGR2 has been shown to be a crucial mediator of allergic asthma, as its expression is necessary in mast cells to direct CD4+ T cell migration to inflamed lung [26]. Though EGR2 negatively regulates T-cell activation, it is required for TH1, TH2, TH9, TH17, and cytotoxic T lymphocyte (CTL) differentiation and cytokine production [27]. TH2-mediated immunity was also represented in the uniquely expressed PBMC transcripts, as the GATA3 transcription factor and IL-4 cytokine were detected in PBMC but not whole blood samples, and both are crucial in the differentiation of naïve helper T cells towards a TH2 orientation.
Many of the statistically significant transcripts in the asthma whole blood cohort were related to cytotoxicity and natural killer (NK) cells. Differentiation of cytotoxic T-lymphocytes (CTLs) requires TBX21, RUNX3, EOMES, and STAT4 [28], all of which are associated with TH1-mediated immunity and were significantly downregulated in addition to the NK surface markers KIR2DS1, KIR3DS1, KLRC1, KLRD1, KLRG1, KLRK1. Transcripts encoding for cytolytic proteins such as PRF1, GZMB, and GZMM were all found to be significantly downregulated in the asthmatic cohort. Although CTLs and NK cells are both mononuclear cells, the differential expression of these transcripts was not detected in PBMCs. CD9 was also upregulated, and it has been used as a surface marker of eosinophils in asthma [29]. Overall, the differentially expressed transcripts found in whole blood describe a downregulation of TH1 differentiation and cell-mediated immunity. For a list of gene/protein abbreviations and their official HUGO names, please see Additional file 4.
The squared correlation coefficient of RNA transcript expression to protein expression has been reported to be ~ 0.4, meaning that mRNA levels account for roughly 40% of the variation in protein concentrations [30]. In our 30 measured proteins extracted from PBMCs, we calculated a squared coefficient of correlation of 0.12 between gene transcripts and proteins. However, we identified 12 proteins with strong correlations to their RNA abundance, which is consistent with a gene-specific level of protein expression [31].
The analysis of differential gene and protein expression identified an interesting insight into the pathophysiology of asthma. The two differentially expressed proteins in our asthmatic cohort were CD28 and OX40 (TNFRSF4), both of which are critical co-stimulatory molecules during T-cell activation. CD28 knockout mice do not develop an inflammatory response during allergen challenge [32], and OX40 is essential for optimizing CD4 and CD8 T-cell responses upon activation [33]. In our gene expression data, the ligand for OX40 (OX40L, TNFSF4) was also significantly upregulated in the whole blood samples of our asthmatic cohort. The OX40/OX40L axis is a central player in promoting TH2 polarization of naïve T-cells in the lymph node, with knockout mouse models of either molecule showing a marked reduction in airway hyperresponsiveness, eosinophilia, and pulmonary inflammation when challenged with ovalbumin [34, 35]. OX40 is expressed on T-cells, but can also be found on NK cells, NKT cells, and neutrophils [36], whereas OX40L is primarily expressed on antigen presenting cells such as B cells, macrophages, and dendritic cells but can be induced on other cell types such as mast cells [37], basophils [38], NK cells [39], and neutrophils [40]. Given the predominant role of TH2 cells in the mechanisms of asthma, it is no surprise that its activator OX40 is found to be upregulated in PBMCs of asthmatics and that widespread OX40L gene expression from a variety of cell types can be detected in whole blood, suggesting a consistently upregulated phenotype of T-cell activation.
In summary, we show here that PAXgene-preserved whole blood is excellent for producing gene expression data with minimal variability and good sensitivity compared to PBMC samples. Despite the small sample size of our study, our paired analysis of whole blood and PBMC samples increased our power to detect differences. Studies requiring blood RNA samples across multiple sites should consider RNA from whole blood instead of PBMC due to its ease of processing (obtained directly from blood draw vs Ficoll isolation) and storage. Although we were limited by the number of transcripts examined, our gene panel did not contain globin genes and thus our data was not hindered by the suppressive effects of globin on the measurement of other transcripts. An added benefit of utilizing the NanoString platform is its precision and robustness; previous studies examining its reproducibility have found that it is comparable to quantitative PCR (qPCR) and exceeds that of RNA-seq and microarrays [41,42,43]. For the first time, we have shown gene expression profiling differences between blood compartments when comparing asthmatic subjects with non-asthmatic healthy controls. Furthermore, in whole blood transcript and PBMC protein expression, we identified the upregulation of the OX40/OX40L axis in the peripheral blood of asthmatic individuals, which has previously been identified as upregulated in serum of pediatric asthmatic patients [44, 45] and investigated as a potential therapeutic target in asthma [46].
The gene expression dataset supporting the conclusions of this article is available in the Gene Expression Omnibus (GEO), GSE132006 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE132006), while normalized protein data is included within the article (Additional file 1).
differentially expressed gene
FDR:
false discovery rate
FEV1 :
forced expiratory volume, 1 s
PBMC:
peripheral blood mononuclear cell
PC20 :
provocative concentration for 20% fall in FEV1
quantitative polymerase chain reaction
U-BIOPRED:
Unbiased BIOmarkers in PREDiction of respiratory disease outcomes
Anandan C, Nurmatov U, van Schayck OC, Sheikh A. Is the prevalence of asthma declining? Systematic review of epidemiological studies. Allergy. 2010;65(2):152–67.
CAS PubMed Article PubMed Central Google Scholar
Gibson GJ, Loddenkemper R, Sibille Y, Lundback B. The European lung white book: respiratory health and disease in Europe. Sheffield: European Respiratory Society; 2013.
Barnett SB, Nurmagambetov TA. Costs of asthma in the United States: 2002–2007. J Allergy Clin Immunol. 2011;127(1):145–52.
Lambrecht BN, Hammad H. The immunology of asthma. Nat Immunol. 2015;16(1):45–56.
Woodruff PG, Modrek B, Choy DF, Jia G, Abbas AR, Ellwanger A, et al. T-helper type 2-driven inflammation defines major subphenotypes of asthma. Am J Respir Crit Care Med. 2009;180(5):388–95.
Bigler J, Boedigheimer M, Schofield JPR, Skipp PJ, Corfield J, Rowe A, et al. A severe asthma disease signature from gene expression profiling of peripheral blood from U-BIOPRED cohorts. Am J Respir Crit Care Med. 2017;195(10):1311–20.
Shin SW, Oh TJ, Park SM, Park JS, Jang AS, Park SW, et al. Asthma-predictive genetic markers in gene expression profiling of peripheral blood mononuclear cells. Allergy Asthma Immunol Res. 2011;3(4):265–72.
Hakonarson H, Bjornsdottir US, Halapi E, Bradfield J, Zink F, Mouy M, et al. Profiling of genes expressed in peripheral blood mononuclear cells predicts glucocorticoid sensitivity in asthma patients. Proc Natl Acad Sci USA. 2005;102(41):14789–94.
Diamant Z, Gauvreau GM, Cockcroft DW, Boulet LP, Sterk PJ, de Jongh FH, et al. Inhaled allergen bronchoprovocation tests. J Allergy Clin Immunol. 2013;132(5):1045–1055.e6.
Singh A, Yamamoto M, Ruan J, Choi JY, Gauvreau GM, Olek S, et al. Th17/Treg ratio derived using DNA methylation analysis is associated with the late phase asthmatic response. Allergy Asthma Clin Immunol. 2014;10(1):32.
PubMed PubMed Central Article CAS Google Scholar
Subrahmanyam YV, Yamaga S, Prashar Y, Lee HH, Hoe NP, Kluger Y, et al. RNA expression patterns change dramatically in human neutrophils exposed to bacteria. Blood. 2001;97(8):2457–68.
Hinkle DE, Wiersma W, Jurs SG. Applied statistics for behavioural sciences. Boston: Houghton Mifflin; 2003.
Bondar G, Cadeiras M, Wisniewski N, Maque J, Chittoor J, Chang E, et al. Comparison of whole blood and peripheral blood mononuclear cell gene expression for evaluation of the perioperative inflammatory response in patients with advanced heart failure. PLoS ONE. 2014;9(12):e115097.
Debey S, Schoenbeck U, Hellmich M, Gathof BS, Pillai R, Zander T, et al. Comparison of different isolation techniques prior gene expression profiling of blood derived cells: impact on physiological responses, on overall expression and the role of different cell types. Pharmacogenom J. 2004;4(3):193–207.
Feezor RJ, Baker HV, Mindrinos M, Hayden D, Tannahill CL, Brownstein BH, et al. Whole blood and leukocyte RNA isolation for gene expression analyses. Physiol Genomics. 2004;19(3):247–54.
Min JL, Barrett A, Watts T, Pettersson FH, Lockstone HE, Lindgren CM, et al. Variability of gene expression profiles in human blood and lymphoblastoid cell lines. BMC Genomics. 2010;11:96.
Palmer C, Diehn M, Alizadeh AA, Brown PO. Cell-type specific gene expression profiles of leukocytes in human peripheral blood. BMC Genomics. 2006;7:115.
Joehanes R, Johnson AD, Barb JJ, Raghavachari N, Liu P, Woodhouse KA, et al. Gene expression analysis of whole blood, peripheral blood mononuclear cells, and lymphoblastoid cell lines from the Framingham Heart Study. Physiol Genomics. 2012;44(1):59–75.
Robinson MD, Speed TP. A comparison of Affymetrix gene expression arrays. BMC Bioinform. 2007;8:449.
Sakharkar MK, Chow VT, Kangueane P. Distributions of exons and introns in the human genome. Silico Biol. 2004;4(4):387–93.
Ray N, Kuwahara M, Takada Y, Maruyama K, Kawaguchi T, Tsubone H, et al. c-Fos suppresses systemic inflammatory response to endotoxin. Int Immunol. 2006;18(5):671–7.
Kim SH, Han SY, Azam T, Yoon DY, Dinarello CA. Interleukin-32: a cytokine and inducer of TNFalpha. Immunity. 2005;22(1):131–42.
Lauterbach H, Bathke B, Gilles S, Traidl-Hoffmann C, Luber CA, Fejer G, et al. Mouse CD8alpha+ DCs and human BDCA3+ DCs are major producers of IFN-lambda in response to poly IC. J Exp Med. 2010;207(12):2703–17.
Proost P, Menten P, Struyf S, Schutyser E, De Meester I, Van Damme J. Cleavage by CD26/dipeptidyl peptidase IV converts the chemokine LD78beta into a most efficient monocyte attractant and CCR1 agonist. Blood. 2000;96(5):1674–80.
Gomez-Martin D, Diaz-Zamudio M, Galindo-Campos M, Alcocer-Varela J. Early growth response transcription factors and the modulation of immune response: implications towards autoimmunity. Autoimmun Rev. 2010;9(6):454–8.
Wu Z, Macneil AJ, Junkins R, Li B, Berman JN, Lin TJ. Mast cell FcepsilonRI-induced early growth response 2 regulates CC chemokine ligand 1-dependent CD4+ T cell migration. J Immunol (Baltimore, Md: 1950). 2013;190(9):4500–7.
Du N, Kwon H, Li P, West EE, Oh J, Liao W, et al. EGR2 is critical for peripheral naive T-cell differentiation and the T-cell response to influenza. Proc Natl Acad Sci USA. 2014;111(46):16484–9.
Cruz-Guilloty F, Pipkin ME, Djuretic IM, Levanon D, Lotem J, Lichtenheld MG, et al. Runx3 and T-box proteins cooperate to establish the transcriptional program of effector CTLs. J Exp Med. 2009;206(1):51–9.
Fernvik E, Halldéen G, Hed J, Lundahl J. Intracellular and surface distribution of CD9 in human eosinophils. Apmis. 1995;103(7–8):699–706.
Vogel C, Marcotte EM. Insights into the regulation of protein abundance from proteomic and transcriptomic analyses. Nat Rev Genet. 2012;13(4):227–32.
Schwanhausser B, Busse D, Li N, Dittmar G, Schuchhardt J, Wolf J, et al. Global quantification of mammalian gene expression control. Nature. 2011;473(7347):337–42.
PubMed Article CAS PubMed Central Google Scholar
Chapoval SP, David CS. CD28 costimulation is critical for experimental allergic asthma in HLA-DQ8 transgenic mice. Clin Immunol (Orlando, Fla). 2003;106(2):83–94.
Redmond WL, Ruby CE, Weinberg AD. The role of OX40-mediated co-stimulation in T-cell activation and survival. Crit Rev Immunol. 2009;29(3):187–201.
Jember AG, Zuberi R, Liu FT, Croft M. Development of allergic inflammation in a murine model of asthma is dependent on the costimulatory receptor OX40. J Exp Med. 2001;193(3):387–92.
Hoshino A, Tanaka Y, Akiba H, Asakura Y, Mita Y, Sakurai T, et al. Critical role for OX40 ligand in the development of pathogenic Th2 cells in a murine model of asthma. Eur J Immunol. 2003;33(4):861–9.
Croft M, So T, Duan W, Soroosh P. The significance of OX40 and OX40L to T-cell biology and immune disease. Immunol Rev. 2009;229(1):173–91.
Kashiwakura J, Yokoi H, Saito H, Okayama Y. T cell proliferation by direct cross-talk between OX40 ligand on human mast cells and OX40 on human T cells: comparison of gene expression profiles between human tonsillar and lung-cultured mast cells. J Immunol. 2004;173(8):5247–57.
Di C, Lin X, Zhang Y, Zhong W, Yuan Y, Zhou T, et al. Basophil-associated OX40 ligand participates in the initiation of Th2 responses during airway inflammation. J Biol Chem. 2015;290(20):12523–36.
Zingoni A, Sornasse T, Cocks BG, Tanaka Y, Santoni A, Lanier LL. Cross-talk between activated human NK cells and CD4+ T cells via OX40-OX40 ligand interactions. J Immunol (Baltimore, Md: 1950). 2004;173(6):3716–24.
Karulf M, Kelly A, Weinberg AD, Gold JA. OX40 ligand regulates inflammation and mortality in the innate immune response to sepsis. J Immunol (Baltimore, Md: 1950). 2010;185(8):4856–62.
Veldman-Jones MH, Brant R, Rooney C, Geh C, Emery H, Harbron CG, et al. Evaluating robustness and sensitivity of the NanoString Technologies nCounter platform to enable multiplexed gene expression analysis of clinical samples. Cancer Res. 2015;75(13):2587–93.
Malkov VA, Serikawa KA, Balantac N, Watters J, Geiss G, Mashadi-Hossein A, et al. Multiplexed measurements of gene signatures in different analytes using the Nanostring nCounter Assay System. BMC Res Notes. 2009;2:80.
Prokopec SD, Watson JD, Waggott DM, Smith AB, Wu AH, Okey AB, et al. Systematic evaluation of medium-throughput mRNA abundance platforms. RNA. 2013;19(1):51–62.
Ma SL, Zhang L. Elevated serum OX40L is a biomarker for identifying corticosteroid resistance in pediatric asthmatic patients. BMC Pulm Med. 2019;19(1):66.
Ezzat MH, Imam SS, Shaheen KY, Elbrhami EM. Serum OX40 ligand levels in asthmatic children: a potential biomarker of severity and persistence. Allergy Asthma Proc. 2011;32(4):313–8.
Gauvreau GM, Boulet LP, Cockcroft DW, FitzGerald JM, Mayers I, Carlsten C, et al. OX40L blockade and allergen-induced airway responses in subjects with mild asthma. Clin Exp Allergy. 2014;44(1):29–37.
The authors thank the research participants for volunteering to take part in these studies, as well as Christopher DeHeer (NanoString) for technical assistance.
Support for this project was provided by AllerGen NCE and the BC Lung Association, as well as an unrestricted grant-in-aid from NanoString for the Vantage 3D assays.
Centre for Heart Lung Innovation, University of British Columbia, Room 166, 1081 Burrard Street, Vancouver, BC, V6Z1Y6, Canada
Daniel He, Chen Xi Yang, Basak Sahin, Amrit Singh & Scott J. Tebbutt
Prevention of Organ Failure (PROOF) Centre of Excellence, Vancouver, BC, V6Z2K5, Canada
Daniel He, Chen Xi Yang, Amrit Singh, Casey P. Shannon & Scott J. Tebbutt
Department of Medicine, McMaster University, Hamilton, ON, L8N3Z5, Canada
John-Paul Oliveria & Gail M. Gauvreau
Department of Pathology, Stanford University, Palo Alto, CA, 94043, USA
John-Paul Oliveria
Department of Medicine (Division of Respiratory Medicine), University of British Columbia, Vancouver, BC, V6Z1Y6, Canada
Scott J. Tebbutt
Daniel He
Chen Xi Yang
Basak Sahin
Amrit Singh
Casey P. Shannon
Gail M. Gauvreau
J-PO and GMG recruited and collected participant samples; CXY, BS, and SJT designed and performed experiments; DH and CXY analyzed the data; DH, AS, CPS, J-PO, GMG, and SJT discussed results; and DH and SJT wrote the manuscript. All authors read and approved the final manuscript.
Correspondence to Scott J. Tebbutt.
All study subjects provided written, informed consent prior to sample collection, and the study was approved by the research ethics boards of McMaster University and the University of British Columbia.
Normalized protein data obtained from healthy controls (HC) and allergic asthmatics (AA). Raw protein values were normalized using total sum scaling and IgG negative controls.
Differentially expressed genes between allergic asthmatics and healthy controls in PBMC and whole blood samples. Comparison of gene expression between asthmatics and healthy controls using linear modeling (FDR cutoff set at 0.25).
Spearman correlation of RNA expression values measured in whole blood and PBMC with measured protein expression in PBMCs. Correlation analysis of RNA expression values measured in whole blood and PBMC samples with protein expression values in PBMCs. Bold indicates significant correlation (p < 0.05).
List of official HUGO gene/protein names. Table of gene/protein abbreviations with their corresponding official Human Gene Organization (HUGO) names.
He, D., Yang, C.X., Sahin, B. et al. Whole blood vs PBMC: compartmental differences in gene expression profiling exemplified in asthma. Allergy Asthma Clin Immunol 15, 67 (2019). https://doi.org/10.1186/s13223-019-0382-x