- Research article
- Open Access
- Open Peer Review
Gene expression analysis in asthma using a targeted multiplex array
BMC Pulmonary Medicinevolume 17, Article number: 189 (2017)
Gene expression changes in the structural cells of the airways are thought to play a role in the development of asthma and airway hyperresponsiveness. This includes changes to smooth muscle contractile machinery and epithelial barrier integrity genes. We used a targeted gene expression arrays to identify changes in the expression and co-expression of genes important in asthma pathology.
RNA was isolated from the airways of donor lungs from 12 patients with asthma (8 fatal) and 12 non-asthmatics controls and analyzed using a multiplexed, hypothesis-directed platform to detect differences in gene expression. Genes were grouped according to their role in airway dysfunction: airway smooth muscle contraction, cytoskeleton structure and regulation, epithelial barrier function, innate and adaptive immunity, fibrosis and remodeling, and epigenetics.
Differential gene expression and gene co-expression analyses were used to identify disease associated changes in the airways of asthmatics. There was significantly decreased abundance of integrin beta 6 and Ras-Related C3 Botulinum Toxin Substrate 1 (RAC1) in the airways of asthmatics, genes which are known to play an important role in barrier function. Significantly elevated levels of Collagen Type 1 Alpha 1 (COL1A1) and COL3A1 which have been shown to modulate cell proliferation and inflammation, were found in asthmatic airways. Additionally, we identified patterns of differentially co-expressed genes related to pathways involved in virus recognition and regulation of interferon production. 7 of 8 pairs of differentially co-expressed genes were found to contain CCCTC-binding factor (CTCF) motifs in their upstream promoters.
Changes in the abundance of genes involved in cell-cell and cell-matrix interactions could play an important role in regulating inflammation and remodeling in asthma. Additionally, our results suggest that alterations to the binding site of the transcriptional regulator CTCF could drive changes in gene expression in asthmatic airways. Several asthma susceptibility loci are known to contain CTCF motifs and so understanding the role of this transcription factor may expand our understanding of asthma pathophysiology and therapeutic options.
Asthma is a chronic inflammatory disease of the airways, characterized by symptoms of breathlessness, wheezing and cough, associated with variable airflow limitation and airway hyperresponsiveness (AHR). Asthma is also characterized by airway remodeling which includes goblet cell metaplasia, epithelial damage, subepithelial fibrosis, basement membrane thickening, and increased airway smooth muscle (ASM) mass . The pathogenesis of asthma is believed to involve an interaction between the innate and adaptive immune systems , and phenotypic changes within the epithelial-mesenchymal trophic unit . Genetic and genomic analyses have been used to discover the molecular mechanisms underlying these phenotypic changes. Large genome-wide association studies (GWAS) have reproducibly identified single nucleotide polymorphisms (SNPs) in or near genes predominantly expressed in the airway epithelium and immune cells as susceptibility factors for asthma. The specific genes include gasdermin B (GSDMB) , interleukin 33 (IL33), and thymic stromal lymphopoietin (TSLP) . Another candidate gene ADAM33, expressed in ASM, has been shown to be associated with asthma in linkage analysis . Genetic variants associated with susceptibility for asthma may exert their effect by altering gene expression levels; indeed many of the SNPs associated with asthma and AHR have been shown to be expression quantitative trait loci (eQTL) in lung tissue, epithelial and blood cells, and altered protein expression of some of these genes has been found in cells and tissue from asthmatic individuals .
In this hypothesis driven study we used a reproducible multiplexed technology (Nanostring®) to quantify the expression of 334 genes potentially involved in phenotypic changes in asthmatic airways. This technology is highly sensitive, reproducible, and is suitable for archived tissue specimens as it is insensitive to RNA degradation . We hypothesized that changes in the expression patterns of genes involved in ASM contraction, the cytoskeleton, epithelial barrier function, innate/adaptive immunity, fibrosis and remodeling, and epigenetics would be present in the airway tissue of asthmatics compared to non-asthmatics. Our results suggest that alterations in the expression of genes involved in cell-cell and cell-matrix interactions may contribute to the pathogenesis of asthma, particularly severe asthma. The identification of altered gene co-expression networks may identify changes in transcriptional regulation that could be pathogenic and missed with commonly used analyses for differential expression.
Additional description of methods is provided in the online supplement.
Subject selection and RNA isolation
Human lungs were donated with consent from the IIAM and used with approval from the University of British Columbia and St. Paul’s Hospital ethics committee. Diagnosis of asthma was determined through patient medical history and asthma medication usage as determined by family interview. Non-asthma donor deaths were primarily due to head trauma while 8 of the 12 donors with asthma died during exacerbations of their asthma. The other four donors with asthma died due to other, accidental causes (eg. head trauma). Subject demographics can be seen in Table 1 with full subject characteristics found in Additional file 1: Table S1. After surgical removal the lungs were flushed with Custodiol HTK solution (Odyssey Pharmaceuticals, East Hanover, NJ, USA) and transported on ice by plane. The time between harvesting and arrival at the University of British Columbia was 15–20 h. Tissues from the lungs have been used in previous studies [9,10,11]. Inflated frozen lungs were processed into tissue cores for sectioning on a crytostat. A total of twenty 10 μm thick sections per core were cut and stored at −80 °C until RNA was isolated. Sections 1, 5, 10, 15 and 20 were stained with hematoxylin and eosin (H&E) for morphometric measurements. For the remaining 15 sections, airways and a small amount of surrounding parenchyma were macroscopically dissected using a scalpel (Fisher Scientific® No. 11) for RNA isolation. Samples were only used if the airways seen on the first section were continuous for the 20 sequential sections. Large vessels were avoided. Sample airway is seen in Fig. 1. RNA was isolated using the Qiagen® RNeasy Mini Kit according to manufacturers protocol.
Characterization of airway dimensions
Measurements of ASM, epithelial, collagen, and total wall area in addition to basement membrane thickness were measured to quantify the degree of airway remodeling. The five sections stained with H&E (Fig. 2) from each core were digitally scanned and airway wall compartments were quantified using the Aperio® system (Leica Biosystems, Germany). ASM area, epithelial area and total airway wall area were quantified using a point counting method, where a grid of 4000 points was overlaid onto each airway of interest and the points falling on the area of interest were counted (Image Pro Plus ®, Media Cybernetics, Maryland). The measurements for all airways were normalized to the internal perimeter (Pi) of the airway and averaged across all the sections for each subject. The Pi was measured by tracing along the luminal side of the epithelium. For the measurement of basement membrane thickness and collagen area, two additional sections were cut from each core (same as used for RNA isolation and other measurements) to stain with Masson’s trichrome (basement membrane) and Picrosirius red (collagen). To quantify the basement membrane thickness, a random series of line segments was placed over each image and the thickness of the basement membrane was measured at the points where any line segment crossed the basement membrane. All thickness measurements were made perpendicular to the epithelium. A minimum of 40 measurements were made for each airway. For the measurement of collagen content, Picrosirius red stained slides were visualized under polarized light where the Picrosirius red stain shows birefringence. The collagen appears red on a black background and the amount of collagen was quantified using color segmentation in Image Pro Plus®. All measurements were carried out in a blinded manner.
Gene expression analysis
Expression of mRNA for the 334 candidate genes and 12 housekeeping genes was measured with the Nanostring® system using a custom codeset panel. The most stable housekeeping genes were selected by measuring 12 common housekeeping genes. By comparing the % Coefficient of variation (%CV) across the 12 housekeeping genes we were able to determine that the 5 most stable genes for data normalization were: RNA Polymerase R2A (POLR2A), TATA box binding protein (TBP), Ribosomal Protein L19 (RPL19), Guanine nucleotide-binding protein subunit beta-2-like 1 (GNB2L1), and β-Glucuronidase (GUSB). These 5 genes were also selected because they spanned a range of counts from low (average 198 counts for GUSB) to high (average 22,486 counts for RPL19). Data were normalized in accordance with Nanostring® guidelines. See Additional file 1: Table S2 in the online supplement for a complete list of genes in the panel.
Candidate gene selection
We selected the candidate genes based on a priori hypotheses and grouped them based on their function and potential role in the pathogenesis of asthma and/or AHR. These are: ASM contraction and relaxation, structure and regulation of the cytoskeleton, epithelial barrier function, innate and adaptive immunity, fibrosis and remodeling and epigenetics. The rationale for the choice of groups of genes is provided below and the list of the genes by category is in Additional file 1: Table S2
The contraction and relaxation of ASM and its regulation; MYH11, MLCK, SM-22 and actin have been previously examined in the context of asthma . Many of the genes in the contractile machinery group were taken from the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway for vascular smooth muscle contraction. Additionally, work by Sieck et al.  led to the selection of CD38 and other calcium handling proteins that have been shown to be expressed in cultured ASM cells. Genes within this pathway are involved in either contraction or relaxation of smooth muscle and so have the potential to play a role in AHR in asthma.
The structure and regulation of the cytoskeleton; genes within the cytoskeletal group were selected based on previous work by Gunst et al. . Genes in this group have been shown to be important in transmitting ASM force to the external environment at adherens junctions, in maintaining the actin filament lattice or or regulating ASM stiffness independent of force generation .
Epithelial barrier function; a number of observations suggest that the airway epithelium is disrupted in asthma and that this may in part result from altered repair mechanisms . Disrupted features include detachment of columnar ciliated cells, the presence of epithelial cell aggregates (Creola bodies) in sputum, decreased expression of epithelial cell-cell junction proteins (E-cadherin, ZO-1, protocadherin-1) and increased expression of epithelial repair markers (TGF-β, EGFR and CD44), mucins, and altered expression of repair-associated fucosylated glycoproteins . A defective epithelial barrier may have important consequences in asthma as it is thought to lead to increased accessibility of allergens to immune and structural cells within the mucosal and submucosal spaces.
Innate/Adaptive immunity; innate immune receptors and related mediators have been implicated in the pathogenesis of asthma (e.g. IL-33, TSLP, and ST2) . Furthermore, it has been demonstrated that antiviral immune responses are compromised in airway epithelial cells from asthmatics . Intrinsic differences in innate immune responses in airway epithelial cells may therefore contribute to disease development and exacerbations in response to environmental exposures including allergens, viruses, and air pollution . We therefore determined the expression patterns of all Toll-like receptors (TLRs), Nod-like receptors (NLRs), and Rig-like receptors (RLRs) and related mediators to provide a comprehensive screen of these candidate genes in asthmatic and non-asthmatic airway wall samples.
Fibrosis and remodeling; the role of the myocardin pathway in proliferation of ASM cells  led to the selection of genes in this pathway. In addition, we interrogated members of the Notch family as these genes are integral to airway development and differentiation  and may play important roles in asthma. ECM proteins are altered in the airways of asthmatics and as such, we included a number of genes that code for ECM components thought to be involved in the remodeling of asthmatic airways. In vivo work has shown a role for matrix metalloproteinases (MMPs) in the development of airway inflammation and hyperresponsiveness .
Epigenetics; as a first line of contact with the external environment, the airway epithelium is an attractive target for epigenetic research. Alterations in DNA methylation and histone modifications have been reported in the airway epithelium of asthmatic subjects , however many of the mediators involved have not been studied. We targeted the histone acetyltransferases KAT2A, CREBBP, and EP300 as they are responsible for acetylating lysine 18 on histone 3 which is up regulated in asthmatic epithelial cells . We also focused on AURKA, PRMT5, SUV39H1, and HDAC10 which have been identified to be potentially involved in the pathogenesis of asthma based on preliminary data from an array analysis of epigenetic modifying enzymes (58). A number of genes chosen for this study were previously found to be differentially methylated in preliminary (PTK7, BCL3, DNMT3b, and PTPRO)  and final (CRIP1, STAT5A, FGFR1, S100A2, ITGA2, EGR4, EID1, and IGSF4C) analyses of DNA methylation in asthmatic airway epithelial cells .
Additional genes were added to the list given their reproducible association in asthma GWAS and observation that the SNP’s in these genes act as eQTLs.
Data and statistical analyses
Final Nanostring results were filtered to keep only genes that had an average count of at least 30. Normalized mRNA expression values were compared between asthmatic and non-asthmatic subjects using a linear model with a negative binomial distribution controlling for age, sex, and inhaled corticosteroid use. Differential gene expression data are presented as volcano plots as well as in a summary table showing the top differentially regulated genes. The level of expression of each transcript is not completely independent since there was strong co-expression between the 344 genes. To account for this, we employed the Matrix Spectral Decomposition analysis of Nyholt and Li et al. [27, 28] to identify the effective number of independent genes. This lead to the multiple comparison correction shown in Table 2 that is based on an effective n of 31. Adjusted p-values will be listed as p.adj and genes with a nominal unadjusted p < 0.05 will be listed as p.unadj.
In addition to the analysis of differential gene expression, we performed an analysis of differential co-expression in our data set using the analysis package CoXpress for R. Differential co-expression analysis identifies pairs of genes that are differentially co-expressed i.e. have opposite correlation patterns in cases vs. controls or show correlations in one condition only. Genes that were differentially co-expressed were entered into WebGestalt for pathway enrichment analysis and are presented in table form. Network analyst was used to understand the protein-protein interaction network [PMID: 25,950,236} of all nominally significant genes (p.unadj < 0.05). Data are plotted using GraphPad version 5.04 (La Jolla California USA).
The total number of airways analyzed was 52 in asthmatic and 53 in non-asthmatic subjects; with an average of 4.3 airways per subject (p > 0.05, asthmatic vs. non-asthmatic). The average internal perimeter (Pi) of asthmatic subjects was 5.1 ± 1.5 mm (geometric mean 4.0 ± 1.7 mm) and in non-asthmatic subjects was 5.3 ± 1.5 mm (geometric mean 4.9 ± 1.3 mm) (p > 0.05 for both arithmetic and geometric mean). Airway wall area per unit length of Pi was significantly greater in asthmatics (0.22 ± 0.024 mm2/mm) than in non-asthmatics (0.13 ± 0.019 mm2/mm, p < 0.01). There was also an increase in the ASM area per unit Pi (0.018 ± 0.0024 mm2/mm vs. 0.011 ± 0.0015 mm2/mm, p < 0.05) and in basement membrane thickness (6.9 ± 0.81 mm vs. 3.9 ± 0.73 mm, p < 0.01) in asthmatics versus non-asthmatics. There was no significant difference in the area of epithelium or collagen per unit Pi between donor groups (p > 0.05). Side by side comparison of asthmatic and non-asthmatic airway can be seen in Fig. 2.
Differential gene expression analysis
Gene expression changes in all genes are summarized in Fig. 3 with the candidate gene hypothesis categories plotted in Fig. 4. In total there were 51 genes differentially expressed based on a threshold p-value of p < 0.05 and three genes that were significant after p-value correction (Table 2, p.adj). In brief, there were three genes whose significance reached the adjusted p-value cutoff, Collagen Type 1 Alpha 1 (COL1A1), COL3A1, and integrin beta 6 (ITGB6). The gene for COL1A1 was the most significantly up-regulated gene (1.83-fold increase, p.adj = 0.01) and integrin beta 6 (ITGB6) was the most significantly down-regulated gene (1.29-fold decrease, p.adj = 0.002) (Additional file 1: Table S2). COL1A1 expression was positively associated with the amount of collagen in the airway in asthmatics and non-asthmatics (Collagen/Pi, R2 = 0.2221 and 0.2182 for asthma and non-asthma respectively, p < 0.01) and the thickness of the basement membrane in both groups combined (R2 = 0.1313, p = 0.01). COL3A1 expression was positively associated with ASM/Pi (R2 = 0.2089, p = 0.02) and Collagen/Pi (R2 = 0.2083, p = 0.03) in asthma. ITGB6 expression was negatively associated with basement membrane thickness in asthma (R2 = 0.1801, p = 0.04) and negatively associated with Epithelial area/Pi in both groups combined (R2 = 0.1158, p = 0.02). There was an association between ITGB6 expression and Collagen/Pi in asthma that did not quite reach significance. (R2 = 0.144, p = 0.06). In each hypothesis group there were a number of differentially expressed genes that did not reach significance after p-value adjustment, these included: contractile apparatus structure – Smoothelin (1.40 fold increase, p.unadj = 0.01); regulation of contraction – CD38 (1.66 fold decrease, p.unadj =0.003); cytoskeletal structure – ITGB6 (see above); cytoskeletal regulation – RAC1 (1.60 fold decrease, p.unadj =0.002); epigenetic regulation – PRMT5 (1.52 fold decrease, p.unadj =0.03); epithelial function – LAMC2 (1.59 fold decrease, p.unadj =0.006); fibrosis and remodeling – COL1A1 (see above); innate and adaptive immunity – PTGFR (2.18 fold decrease, p.unadj =0.004). The counts, p-values, and adjusted p-values for all significant genes can be seen in Additional file 1: Table S2. Of the 12 genes identified in GWAS or linkage analysis, only ADAM33 (1.56-fold increase, p.unadj =0.0057) was up-regulated. There were no pathways significantly enriched in the differentially up or down-regulated genes.
Using Network Analyst, we identified a minimum protein-protein interaction network and key nodes from our nominally differentially expressed genes (Fig. 5). Green nodes indicate down-regulated genes, red nodes indicate up-regulated nodes (both relative to non-asthmatics), and grey nodes indicate first order interactions. This network highlights a number of key nodes in our data set including: Mitogen-Activated Protein Kinase 1 (MAPK1, degree = 24), c-FOS (degree = 22, and Calmodulin 3 (CALM3, degree = 21). Using this network, we were able to identify pathways significantly associated with up and down-regulated nodes (Table 3). These included key pathways in collagen degradation and remodeling and alterations to cell-cell communication.
Differential co-expression analysis
We identified groups of genes that were differentially co-expressed between asthmatics and non-asthmatics. In this analysis, genes are clustered together based on how their expression values correlate with each other. The analysis was performed twice with the comparison group being the non-asthmatics or asthmatics in the different analyses. Clusters of genes that are significantly co-expressed in one condition and not in the other are said to be differentially co-expressed. Figure 5 shows an example; expression of genes in cluster 11 changes from subject to subject in asthmatics (left) and non-asthmatics (right). Each line represents one gene in the group. Genes in cluster 11 follow a similar pattern of expression in asthmatics (p < 0.001) but not in non-asthmatics (p = 0.32). The rest of these figures can be seen in Additional file 1: Fig. S3 and S4.
In non-asthmatics, there were 3 groups of genes that were differentially co-expressed compared to asthmatics. These clusters (10, 35, and 53) had an average co-expression correlation coefficient of 0.772 (p < 0.0001) in non-asthmatics and 0.223 (p > 0.05) in asthmatics. The genes in these groupings are summarized in Additional file 1: Table S3. We performed pathways analysis on each group of genes with Webgestalt. These results are shown in Table 4. There were no pathways that were significantly enriched for genes differentially co-expressed in non-asthmatics. Further analysis identified specific gene pairs from cluster 10 that were both significantly positively correlated in non-asthmatics and significantly negatively correlated in asthmatics (Table 5). These included: chitinase 3-like 1(CHI3L1) and GSDMB (R = 0.760 and R = −0.707); CHI3L1 and histone deacetylase 10 (HDAC10) (R = 0.731 and R = −0.676); HDAC10 and thymocyte antigen 1 (THY1 or CD90) (R = 0.727 and R = −0.596); and indoleamine 2,3-dioxygenase 1 (IDO1) and nuclear factor of activated T-cells, cytoplasmic 2 (NFATC2) (R = 0.659 and R = −0.604), in non-asthmatics and asthmatics respectively. All showed significant but opposite direction of correlation (p < 0.05 Fig. 6).
In asthmatic samples there were 6 clusters of genes that were found to be differentially co-expressed. These clusters had an average correlation coefficient of 0.728 (p < 0.0001) in asthmatics and 0.169 (p > 0.05) in non-asthmatics. The genes in these 6 clusters are summarized in Additional file 1: Table S3. Each of the clusters was also analyzed with Webgestalt (Table 4). In brief, asthmatic co-expressed genes were significantly enriched in pathways for cytoplasmic virus pattern recognition signaling (p = 3.0 × 10−4), positive, and negative regulation of type 1 interferon production (p = 3.9 × 10−3, and p = 2.0 × 10−3 respectively). Within the cluster of genes that were differentially co-expressed, there were 4 pairs of genes whose expression were positively and significantly correlated in asthmatics but negatively and significantly correlated in non-asthmatics (Table 5). In brief: laminin B2 (LAMB2) and matrix metallopeptidase 9 (MMP9) (R = 0.727 and R = −0.686); BPI fold containing family A (BPFIA1) and DEAD box polypeptide 58 (DDX58 or RIG-1, retinoic acid inducible gene 1 protein) (R = 0.814 and R = −0.646); bone morphogenetic protein 2 (BMP2) and NFATC2 (R = 0.707 and R = −0.660); and G-protein alpha stimulating (GNAS) and NFATC2 (R = 0.665 and R = −0.798), in asthmatics and non-asthmatic respectively. All were significantly correlated but in opposite directions (p < 0.05, Fig. 6).
All eight pairs of differentially co-expressed genes (Table 5) were entered into the ENCODE ChIP-SEQ significance tool. All but one of the genes (IDO1) from the pairs was regulated by the transcriptional repressor CCCTC-binding factor or CTCF. Since only 47% of the candidate genes have a CTCF binding site this represents a significant enrichment (p = 0.0035, Chi-squared test).
In this study the Nanostring® platform was used to quantify the airway expression of candidate genes hypothesized to be important in the pathophysiology of asthma. mRNA was obtained from the airways of asthmatic (mainly fatal asthmatics) and non-asthmatic donor lungs and measures of airway remodeling were made on the sampled airways. Candidate genes were grouped into the following categories: ASM contraction, the cytoskeleton, epithelial barrier function, innate/adaptive immunity, fibrosis/remodeling, and epigenetics. 51 genes (15%) were nominally differentially expressed (p.unadj <0.05) in asthmatic airway tissue and included many genes important in cell-cell and cell-matrix interactions (COL1A1, COL3A1, ITGB6, LAMC2, RAC1). Of the 51 genes differentially expressed based on a nominal p-value, only 3 were significant following multiple comparison correction (ITGB6, COL1A1, COL3A1).
Cell-cell junctions are altered in the airway epithelium of asthmatics  and this may result in greater permeability of the epithelial layer and ultimately hypersensitivity of the ASM to agonist challenge . The most significant differentially expressed gene in the data set was ITGB6. ITGB6 rapidly accumulates following injury to the epithelial layer and is considered to be important for normal wound healing . In the mouse loss of ITGB6 causes an increase in the number of B-cells and T-cells around the airways, an increase in IL-4 production, and airway hyperresponsivenss in naïve mice without allergen challenge . Additionally, influenza virus has been shown to interact with ITGB6 to cause epithelial cell death and collagen deposition in a TGF-ß dependent manner . Another gene with lower abundance in asthmatics that did not meet the adjusted p-value cut off was RAC1. RAC1 has been shown to be important in the formation of tight junctions in an EGFR dependent manner , specifically by regulating tight junction protein 1 (or zona occluden 1)  which was also lower in abundance in the asthmatic samples. Changes in these genes could also affect epithelial-mesenchymal transition (EMT) , cytoskeletal stability, actin filament assembly/disassembly, cell stiffness and/or cell migration . In addition, a protein-protein interaction network highlighted that down-regulated genes were enriched in pathways for cell-cell communication and integrin cell surface interactions.
The most significantly up-regulated gene in asthmatic samples was COL1A1 which codes for the alpha chain in type 1 collagen. Collagen 1 is the major type of collagen in basement membrane and a major protein found in remodeled airways. Collagen 1 is important in airway remodeling , in particular thickening of the subepithelial space which is associated with worsening of asthma symptoms . Collagen type 3 was also significantly elevated in asthmatic subjects and is also significantly elevated in the basement membrane of asthmatic subjects . Type 1 and 3 collagen have both been associated with worsening lung function in a horse model of asthma . Beyond the ability of collagens to affect distensibility of the airways, collagen 1 has been shown to stimulate ASM to produce MMP1  and proliferate, in conjunction with FAK . Collagen 1 and 3 expression has also been shown to be unaffected by corticosteroid usage in severe asthmatics  and collagen 1 and 3 contribute to the loss of the anti-mitotic effect of corticosteroids . Enrichment of pathways involved in collagen remodeling was seen in our network analysis and highlights the importance of understanding fibrosis as it relates to fiber production, degradation, and organization and how this impacts normal cell function.
ADAM33 has been described in candidate gene studies to be associated with asthma  and was significantly elevated in asthmatic airways. ADAM33 has been implicated in smooth muscle development, cell-cell connections, and cell proliferation and differentiation . Over-expression of ADAM33 in asthmatic airways has been described  and may be an important determinant of disease progression. There is evidence that ADAM33 can stimulate angiogenesis ex vivo and in vivo and by this mechanisms may contribute to airway remodeling . Furthermore, ADAM33 family member TACE/ADAM17 can mediate release of TNF-α and fracktalkine (or CX3CL1) from the cell membrane [48, 49] and ADAM9 may mediate the release of growth factor HB-EGF . If ADAM33 has a similar capacity for cytokine or growth factor cleavage this could make it a major contributor to airway remodeling in asthma.
Other genes in the 5 most differentially up or down regulated genes that were nominaly significant include: Cyclic ADP Ribose (CD38), Interleukin 13 Receptor Alpha 1 (IL13RA1), Prostaglandin F Receptor (PTGFR), Heat Shock Protein Beta 1 (HSPB1), and Interferon Induced with Helicase C Domain 1 (IFIH1). CD38 is a protein that generates the second messenger cADPR to cause calcium release. Recent work has explored the role of CD38 in asthma and has suggested that increased CD38 expression causes hypercontractility in ASM cells from asthmatics  although in our samples we saw no CD38 staining in the ASM layer (Additional file 1: Fig. S1 and S2). Additionally, CD38 deficient mice have reduced AHR following ovalbumin challenge . Decreased IL13RA1 expression is surprising in the context of asthma but this could be due to a compensatory response to continued eosinophilia and IL13/IL4 exposure in the lung . Prostaglandin’s can be both pro and anti inflammatory but there is little research on the role of prostaglandin F in the context of asthma. HSPB1 (or heat shock protein 27) is a chaperone protein that has been implicated in cellular differentiation, apoptosis, and smooth muscle contraction . Up-regulation of the gene could contribute to ASM hypercontractility in asthma but further work investigating the phosphorylation state and activity of HSP27 in asthma is needed to answer this question. IFIH1, also known as MDA5, is a DEAD box double stranded (ds) RNA helicases that can detect intracellular viral dsRNA and lead to the production of interferons . MDA5 and TLR3 signaling have been shown to be deficient in bronchial epithelial cells from asthmatic subjects  and this could be responsible for the defective epithelial release of interferon I and III in response to rhino virus infection [57, 58].
Co-expression of genes does not imply interaction between their proteins but instead may suggest similarities in their regulation by transcription factors or epigenetic mechanisms . Co-expression analyses can reveal changes in the regulation of gene expression  and have been used to identify epigenetic changes that affect gene co-expression in cancer . In our study, genes that were differentially co-expressed between asthmatics and non-asthmatics were significantly enriched for pathways involved in virus recognition and regulation of interferon production (Table 4). The genes enriched in these pathways were from cluster 16 and were RIG-I (DDX58), RIG-1-like receptor 3 (DHX58), and interferon induced with helicase C domain 1 (IFIH1). This finding suggests that a central molecular mechanism may regulate diverse antiviral immune molecules in response to viral infections that may trigger asthma exacerbations and/or pathogenesis . IFIH1, as discussed earlier, was also differentially expressed.
One of the most intriguing results was the identification of a single transcriptional repressor, CTCF, that controls the expression of all but one of the differentially co-expressed pairs of genes. CTCF influences gene expression through chromatin modifications  resulting in insulation of the target regions . CTCF is an architectural protein that mediates inter- and intra-chromosomal interactions at distant genomic sites, and regulates three-dimensional genome architecture . There are examples of CTCF silencing one gene while activating another . Specific to asthma, differential expression at the ZPBP2/GSDMB/ORMDL3 locus was identified resulting from allele-specific chromatin remodeling mediated by CTCF . A SNP in ZPBP2 created a CTCF binding site resulting in increased expression of ZPBP2 but diminished expression of GSDMB and ORMDL3 . Additionally, CTCF is highly sensitive to DNA methylation at CTCF binding sites ; changes to the methylome can have direct effects on the regulation of CTCF target genes. CTCF could play a crucial role in controlling the many gene expression changes observed in the airways of asthmatics and is worthy of more intense research that is beyond the scope of this paper.
There are several limitations of this study. Firstly, the use of whole airway RNA rather than RNA from specific cell types precludes us from conclusively identifying the site of gene expression. Secondly, the majority of the asthmatic patients were fatal asthmatics and experienced hypoxia and treatment with steroids during their fatal attack which can affect mRNA expression in tissues taken for research purposes. We addressed this by controlling for steroid use in our analysis of differential gene expression. Additionally, the use of more severe asthmatics may mean that these results are not generalizable to asthmatics as a whole. But considering severe asthmatic populations have the most hospital visits and are most at risk for exacerbations, we believe our results provide significant insight into the genes that are altered in fatal disease. Finally, the relatively small sample size limits our ability to detect differences in gene expression less than ~1.5 fold on average, although this also means that the changes we see are likely to be real. Procurement of donor lungs is time consuming and costly so increasing the number of patients for this study was not feasible, however future studies in asthmatic biopsies or cell culture experiments could confirm these results with the ability for much larger sample size. A limited number of donors with non-fatal asthma (n = 4) means we were unable to test for differences between these two groups of donors (fatal vs. non-fatal).
This study identifies changes in the expression and co-expression of genes thought to be important in asthma and AHR. Specifically, we identified altered abundance of genes involved in cell-cell and cell-matrix connections as well as those involved in the immune response and cell homeostasis. We also identified changes in the co-expression of genes involved in virus recognition and interferon production. The transcription factor CTCF could be an important regulator of the asthmatic phenotype and warrants further investigation. Future work should focus on elucidating the potential mechanisms behind altered CTCF binding as it relates to asthma pathophysiology.
Percent coefficient of variation
ADAM metallopeptidase domain 17
ADAM metallopeptidase domain 33
ADAM metallopeptidase domain 9
Airway smooth muscle
Aurora kinase A
B-Cell CLL/lymphoma 3
Bone morphogenetic protein 2
BPI fold containing family A member 1
Cyclic ADP Ribose
Cluster of differentiation 38
Cluster of differentiation 44
Collagen Type 1 Alpha 1
Collagen Type 3 Alpha 1
CREB binding protein
Cysteine rich protein 1
DExD/H-box helicase 58
DNA methyltransferase 3 Beta
Epidermal growth factor receptor
Early growth response 4
EP300 interacting inhibitor of differentiation 1
E1A binding protein P300
Expression quantitative trait loci
Focal adhesion kinase
Fibroblast growth factor receptor 1
Adenylate cyclase-stimulating G alpha protein
Guanine nucleotide-binding protein subunit Beta-2-Like 1
Genome wide association study
Hemoatoxylin and eosin
Heparin-binding EGF-like growth factor
Histone deacetylase 10
Heat shock protein 27
Indoleamine 2,3-dioxygenase 1
Interferon induced with helicase C domain 1
Immunoglobulin superfamily member 4C
International Institute of the Advancement of Medicine
Integrin subunit Alpha 2
Integrin subunit Beta 6
Lysine acetyltransferase 2A
Kyoto encyclopedia of genes and genomes
Laminin subunit Beta 2
Laminin subunit Gamma 2
Melanoma differentiation-associated protein 5
Myosin light chain kinase
Matrix metallopeptidase 1
Matrix metallopeptidase 9
Myosin heavy chain 11
Nuclear factor of activated T-cells, cytoplasmic 2
ORMDL sphingolipid biosynthesis regulator 3
Internal perimeter of the airway
RNA polymerase R2A
Protein arginine methyltransferase 5
Prostaglandin F receptor
Protein tyrosine kinase 7
Protein tyrosine phosphatase, receptor type O
Ras-related C3 botulinum toxin substrate 1
Retinoic acid inducible gene 1
Ribosomal protein L19
S100 calcium binding protein A2
Standard error of the mean
Smooth muscle protein 22-alpha/transgelin
Single nucleotide polymorphism
Suppression of tumorigenicity 2
Signal transducer and activator of transcription 5A
Suppressor of variegation 3–9 homolog 1
Tumor necrosis factor-α-converting enzyme
TATA box binding protein
Transforming growth factor Beta 1
Thymocyte antigen 1
- TNF- α:
Tumor necoris factor alpha
Thymic stromal lymphopoietin
Zona pellucida binding protein 2
James AL, Elliot JG, Jones RL, Carroll ML, Mauad T, Bai TR, et al. Airway smooth muscle hypertrophy and hyperplasia in asthma. Am J Respir Crit Care Med. 2012;185:1058–64.
Barnes PJ. Inhaled Glucocorticoids for asthma. N Engl J Med. 1995;332:868–75.
Holgate ST, Holloway J, Wilson S, Bucchieri F, Puddicombe S, Davies DE. Epithelial–Mesenchymal communication in the pathogenesis of chronic asthma. Proc Am Thorac Soc. 2004;1:93–8.
Tulah AS, Holloway JW, Sayers I. Defining the contribution of SNPs identified in asthma GWAS to clinical variables in asthmatic children. BMC Med Genet. 2013;14:100.
Torgerson DG, Ampleford EJ, Chiu GY, Gauderman WJ, Gignoux CR, Graves PE, et al. Meta-analysis of genome-wide association studies of asthma in ethnically diverse north American populations. Nat Genet. 2011;43:887–92.
van Eerdewegh P, Little RD, Dupuis J, del Mastro RG, Falls K, Simon J, et al. Association of the ADAM33 gene with asthma and bronchial hyperresponsiveness. Nature. 2002;418:426.
Hao K, Bosse Y, Nickle DC, Pare PD, Postma DS, Laviolette M, et al. Lung eQTLs to help reveal the molecular underpinnings of asthma. PLoS Genet. 2012;8:e1003029.
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:2587–93.
Chin LY, Bosse Y, Pascoe CD, Hackett TL, Seow CY, Pare PD. Mechanical properties of asthmatic airway smooth muscle. Eur Respir J. 2012;40:45–54.
Hackett T-L, Warner SM, Stefanowicz D, Shaheen F, Pechkovsky DV, Murray LA, et al. Induction of epithelial-mesenchymal transition in primary airway epithelial cells from patients with asthma by transforming growth factor-beta1. Am J Respir Crit Care Med. 2009;180:122–33.
Hirota JA, Hirota SA, Warner SM, Stefanowicz D, Shaheen F, Beck PL, et al. The airway epithelium nucleotide-binding domain and leucine-rich repeat protein 3 inflammasome is activated by urban particulate matter. J Allergy Clin Immunol. 2012;129:1116–25. e6
Woodruff PG, Dolganov GM, Ferrando RE, Donnelly S, Hays SR, Solberg OD, et al. Hyperplasia of smooth muscle in mild to moderate asthma without changes in cell size or gene expression. Am J Respir Crit Care Med. 2004;169:1001–6.
Sieck GC, White TA, Thompson MA, Pabelick CM, Wylam ME, Prakash YS. Regulation of store-operated Ca2+ entry by CD38 in human airway smooth muscle. Am J Physiol - Lung Cell Mol Physiol. 2008;294:L378–85.
Gunst SJ, Zhang W. Actin cytoskeletal dynamics in smooth muscle: a new paradigm for the regulation of smooth muscle contraction. Am J Physiol - Cell Physiol. 2008;295:C576–87.
Lan B, Wang L, Zhang J, Pascoe CD, Norris BA, Liu JC-Y, et al. Rho-kinase mediated cytoskeletal stiffness in skinned smooth muscle. J Appl Physiol Bethesda Md. 1985;2013 https://doi.org/10.1152/japplphysiol.00654.2013.
Allahverdian S, Harada N, Singhera GK, Knight DA, Dorscheid DR. Secretion of IL-13 by airway epithelial cells enhances epithelial repair via HB-EGF. Am J Respir Cell Mol Biol. 2008;38:153–60.
Dorscheid DR, Wojcik KR, Yule K, White SR. Role of cell surface glycosylation in mediating repair of human airway epithelial cell monolayers. Am J Physiol Lung Cell Mol Physiol. 2001;281:L982–92.
Moffatt MF, Gut IG, Demenais F, Strachan DP, Bouzigon E, Heath S, et al. A large-scale, consortium-based genomewide association study of asthma. N Engl J Med. 2010;363:1211–21.
Xu J, Singhera GK, Dorscheid DR. Expression of surfactant protein D in airways of asthmatics and interleukin-13 modulation of surfactant protein D in human models of airway epithelium. Respir Res. 2015;16:26.
Hirota JA, Gold MJ, Hiebert PR, Parkinson LG, Wee T, Smith D, et al. The nucleotide-binding domain, leucine-rich repeat protein 3 inflammasome/IL-1 receptor I axis mediates innate, but not adaptive, immune responses after exposure to particulate matter under 10 μm. Am J Respir Cell Mol Biol. 2015;52:96–105.
Wang D-Z, Olson EN. Control of smooth muscle development by the myocardin family of transcriptional coactivators. Curr Opin Genet Dev. 2004;14:558.
Tsao P-N, Vasconcelos M, Izvolsky KI, Qian J, Lu J, Cardoso WV. Notch signaling controls the balance of ciliated and secretory cell fates in developing airways. Dev Camb Engl. 2009;136:2297–307.
Kumagai K, Ohno I, Okada S, Ohkawara Y, Suzuki K, Shinya T, et al. Inhibition of matrix metalloproteinases prevents allergen-induced airway inflammation in a murine model of asthma. J Immunol Baltim Md 1950. 1999;162:4212–9.
Stefanowicz D, Lee JY, Lee K, Shaheen F, Koo H-K, Booth S, et al. Elevated H3K18 acetylation in airway epithelial cells of asthmatic subjects. Respir Res. 2015;16:95.
Stefanowicz D. Profiling the epigenetic landscape of the airway epithelium in asthma. 2014. https://circle.ubc.ca/handle/2429/47071. Accessed 2 Oct 2015.
Stefanowicz D, Hackett T-L, Garmaroudi FS, Günther OP, Neumann S, Sutanto EN, et al. DNA methylation profiles of airway epithelial cells and PBMCs from healthy, atopic and asthmatic children. PLoS One. 2012;7:e44213.
Nyholt DR. A simple correction for multiple testing for single-nucleotide polymorphisms in linkage disequilibrium with each other. Am J Hum Genet. 2004;74:765–9.
Li J, Ji L. Adjusting multiple testing in multilocus analyses using the eigenvalues of a correlation matrix. Heredity. 2005;95:221–7.
de Boer WI, Sharma HS, Baelemans SMI, Hoogsteden HC, Lambrecht BN, Braunstahl GJ. Altered expression of epithelial junctional proteins in atopic asthma: possible role in inflammation. Can J Physiol Pharmacol. 2008;86:105–12.
Barnes PJ, Cuss FM, Palmer JB. The effect of airway epithelium on smooth muscle contractility in bovine trachea. Br J Pharmacol. 1985;86:685–91.
Huang XZ, Wu JF, Cass D, Erle DJ, Corry D, Young SG, et al. Inactivation of the integrin beta 6 subunit gene reveals a role of epithelial integrins in regulating inflammation in the lung and skin. J Cell Biol. 1996;133:921–8.
Jolly L, Stavrou A, Vanderstoken G, Meliopoulos VA, Habgood A, Tatler AL, et al. Influenza promotes collagen deposition via αvβ6 Integrin-mediated transforming growth factor β activation. J Biol Chem. 2014;289:35246–63.
Terakado M, Gon Y, Sekiyama A, Takeshita I, Kozu Y, Matsumoto K, et al. The Rac1/JNK pathway is critical for EGFR-dependent barrier formation in human airway epithelial cells. Am J Physiol Lung Cell Mol Physiol. 2011;300:L56–63.
Jou T-S, Schneeberger EE, Nelson WJ. Structural and functional regulation of tight junctions by RhoA and Rac1 small GTPases. J Cell Biol. 1998;142:101–15.
Hackett T-L. Epithelial-mesenchymal transition in the pathophysiology of airway remodelling in asthma. Curr Opin Allergy Clin Immunol. 2012;12:53–9.
Parsons JT, Horwitz AR, Schwartz MA. Cell adhesion: integrating cytoskeletal dynamics and cellular tension. Nat Rev Mol Cell Biol. 2010;11:633–43.
Liesker JJ, Hacken NHT, Zeinstra-Smith M, Rutgers SR, Postma DS, Timens W. Reticular basement membrane in asthma and COPD: similar thickness, yet different composition. Int J Chron Obstruct Pulmon Dis. 2009;4:127–35.
Chetta A, Foresi A, Del Donno M, Bertorelli G, Pesci A, Olivieri D. Airways remodeling is a distinctive feature of asthma and is related to severity of disease. Chest. 1997;111:852–7.
Wilson JW, Li X. The measurement of reticular basement membrane and submucosal collagen in the asthmatic airway. Clin Exp Allergy. 1997;27:363–71.
Setlakwe EL, Lemos KR, Lavoie-Lamoureux A, Duguay J-D, Lavoie J-P. Airway collagen and elastic fiber content correlates with lung function in equine heaves. Am J Physiol - Lung Cell Mol Physiol. 2014;307:L252–60.
Rogers NK, Clements D, Dongre A, Harrison TW, Shaw D, Johnson SR. Extra-cellular matrix proteins induce matrix Metalloproteinase-1 (MMP-1) activity and increase airway smooth muscle contraction in asthma. PLoS One. 2014;9:e90565.
Dekkers BGJ, Spanjer AIR, van der Schuyt RD, Kuik WJ, Zaagsma J, Meurs H. Focal adhesion Kinase regulates collagen I–induced airway smooth muscle phenotype switching. J Pharmacol Exp Ther. 2013;346:86–95.
Chakir J, Shannon J, Molet S, Fukakusa M, Elias J, Laviolette M, et al. Airway remodeling-associated mediators in moderate to severe asthma: effect of steroids on TGF-β, IL-11, IL-17, and type I and type III collagen expression. J Allergy Clin Immunol. 2003;111:1293–8.
Bonacci JV, Schuliga M, Harris T, Stewart AG. Collagen impairs glucocorticoid actions in airway smooth muscle through integrin signalling. Br J Pharmacol. 2006;149:365–73.
Shapiro SD, Owen CA. ADAM-33 surfaces as an asthma gene. N Engl J Med. 2002;347:936–8.
Foley SC, Mogas AK, Olivenstein R, Fiset PO, Chakir J, Bourbeau J, et al. Increased expression of ADAM33 and ADAM8 with disease progression in asthma. J Allergy Clin Immunol. 2007;119:863–71.
Puxeddu I, Pang YY, Harvey A, Haitchi HM, Nicholas B, Yoshisue H, et al. The soluble form of a disintegrin and metalloprotease 33 promotes angiogenesis: implications for airway remodeling in asthma. J Allergy Clin Immunol. 2008;121:1400–6. 1406.e1–4
Black RA. Tumor necrosis factor-α converting enzyme. Int J Biochem Cell Biol. 2002;34:1–5.
Garton KJ, Gough PJ, Blobel CP, Murphy G, Greaves DR, Dempsey PJ, et al. Tumor necrosis factor-α-converting enzyme (ADAM17) mediates the cleavage and shedding of Fractalkine (CX3CL1). J Biol Chem. 2001;276:37993–8001.
Izumi Y, Hirata M, Hasuwa H, Iwamoto R, Umata T, Miyado K, et al. A metalloprotease-disintegrin, MDC9/meltrin-gamma/ADAM9 and PKCdelta are involved in TPA-induced ectodomain shedding of membrane-anchored heparin-binding EGF-like growth factor. EMBO J. 1998;17:7260–72.
Jude JA, Dileepan M, Panettieri RA, Walseth TF, Kannan MS. Altered CD38/cyclic ADP-ribose Signaling contributes to the asthmatic phenotype. J Allergy. 2012;2012:289468.
Deshpande DA, White TA, Guedes AGP, Milla C, Walseth TF, Lund FE, et al. Altered airway responsiveness in CD38-deficient mice. Am J Respir Cell Mol Biol. 2005;32:149–56.
Myrtek D, Knoll M, Matthiesen T, Krause S, Lohrmann J, Schillinger D, et al. Expression of interleukin-13 receptor alpha 1-subunit on peripheral blood eosinophils is regulated by cytokines. Immunology. 2004;112:597–604.
Salinthone S, Tyagi M, Gerthoffer WT. Small heat shock proteins in smooth muscle. Pharmacol Ther. 2008;119:44–54.
Reikine S, Nguyen JB, Modis Y. Pattern recognition and Signaling mechanisms of RIG-I and MDA5. Front Immunol. 2014;5 https://doi.org/10.3389/fimmu.2014.00342.
Parsons KS, Hsu AC, Wark PAB. TLR3 and MDA5 signalling, although not expression, is impaired in asthmatic epithelial cells in response to rhinovirus infection. Clin Exp Allergy J Br Soc Allergy Clin Immunol. 2014;44:91–101.
Wark PAB, Johnston SL, Bucchieri F, Powell R, Puddicombe S, Laza-Stanca V, et al. Asthmatic bronchial epithelial cells have a deficient innate immune response to infection with rhinovirus. J Exp Med. 2005;201:937–47.
Contoli M, Message SD, Laza-Stanca V, Edwards MR, Wark PAB, Bartlett NW, et al. Role of deficient type III interferon-lambda production in asthma exacerbations. Nat Med. 2006;12:1023–6.
Hsu C-L, Juan H-F, Huang H-C. Functional analysis and characterization of differential Coexpression networks. Sci Rep. 2015;5:13295.
de la Fuente A. From “differential expression” to “differential networking” – identification of dysfunctional regulatory networks in diseases. Trends Genet. 2010;26:326–33.
Ii M, Yamamoto H, Taniguchi H, Adachi Y, Nakazawa M, Ohashi H, et al. Co-expression of laminin β3 and γ2 chains and epigenetic inactivation of laminin α3 chain in gastric cancer. Int J Oncol. 2011;39:593–9.
Porritt RA, Hertzog PJ. Dynamic control of type I IFN signalling by an integrated network of negative regulators. Trends Immunol. 2015;36:150–60.
Ong C-T, Corces VG. CTCF: an architectural protein bridging genome topology and function. Nat Rev Genet. 2014;15:234–46.
Mishiro T, Ishihara K, Hino S, Tsutsumi S, Aburatani H, Shirahige K, et al. Architectural roles of multiple chromatin insulators at the human apolipoprotein gene cluster. EMBO J. 2009;28:1234–45.
Verlaan DJ, Berlivet S, Hunninghake GM, Madore A-M, Larivière M, Moussette S, et al. Allele-specific chromatin Remodeling in the ZPBP2/GSDMB/ORMDL3 locus associated with the risk of asthma and autoimmune disease. Am J Hum Genet. 2009;85:377–93.
Lu Wang PhD for her assistance in collection of histological images.
Research was funded by the Canadian Institute for Health Research (CIHR) and no role in the design, collection, analysis, or interpretation of data and results, and did not contribute to the writing of the manuscript.
Availability of data and materials
Complete gene expression data is available in Additional file 1: Table S2 in the online supplement. Data is also available from corresponding author upon request.
Ethics approval and consent to participate
Human lungs were donated with consent from the IIAM and used with approval from the University of British Columbia and St. Paul’s Hospital ethics committee.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Methods – Methods describing selection of house keeping genes and immunohistochemical staining procedure. Supplementary Tables – Tables containing clinical demographics for subjects, average counts, fold change, and p-value for all genes studied, and all differentially co-expressed genes. Supplementary Figures and Legends – Figures showing sample immunohistochemical staining for proteins of significantly altered genes, co-expression plots. (DOCX 35 kb)