- Open Access
Haploinsufficiency of Hedgehog interacting protein causes increased emphysema induced by cigarette smoke through network rewiring
- Taotao Lao1,
- Kimberly Glass1, 2, 3,
- Weiliang Qiu1,
- Francesca Polverino4, 5, 6,
- Kushagra Gupta4,
- Jarrett Morrow1,
- John Dominic Mancini1,
- Linh Vuong1,
- Mark A Perrella4, 7,
- Craig P Hersh1, 4,
- Caroline A Owen4, 5,
- John Quackenbush1, 2, 3,
- Guo-Cheng Yuan2, 3,
- Edwin K Silverman1, 4Email author and
- Xiaobo Zhou1, 4Email author
© Lao et al.; licensee BioMed Central. 2015
- Received: 16 December 2014
- Accepted: 23 January 2015
- Published: 14 February 2015
The HHIP gene, encoding Hedgehog interacting protein, has been implicated in chronic obstructive pulmonary disease (COPD) by genome-wide association studies (GWAS), and our subsequent studies identified a functional upstream genetic variant that decreased HHIP transcription. However, little is known about how HHIP contributes to COPD pathogenesis.
We exposed Hhip haploinsufficient mice (Hhip +/- ) to cigarette smoke (CS) for 6 months to model the biological consequences caused by CS in human COPD risk-allele carriers at the HHIP locus. Gene expression profiling in murine lungs was performed followed by an integrative network inference analysis, PANDA (Passing Attributes between Networks for Data Assimilation) analysis.
We detected more severe airspace enlargement in Hhip +/- mice vs. wild-type littermates (Hhip +/+ ) exposed to CS. Gene expression profiling in murine lungs suggested enhanced lymphocyte activation pathways in CS-exposed Hhip +/- vs. Hhip +/+ mice, which was supported by increased numbers of lymphoid aggregates and enhanced activation of CD8+ T cells after CS-exposure in the lungs of Hhip +/- mice compared to Hhip +/+ mice. Mechanistically, results from PANDA network analysis suggested a rewired and dampened Klf4 signaling network in Hhip +/- mice after CS exposure.
In summary, HHIP haploinsufficiency exaggerated CS-induced airspace enlargement, which models CS-induced emphysema in human smokers carrying COPD risk alleles at the HHIP locus. Network modeling suggested rewired lymphocyte activation signaling circuits in the HHIP haploinsufficiency state.
- Chronic Obstructive Pulmonary Disease
- Cigarette Smoke
- Cigarette Smoke Exposure
- Hedgehog Pathway
- Lymphoid Aggregate
Chronic obstructive pulmonary disease (COPD), the third leading cause of death in the US , is a complex disease strongly influenced by cigarette smoke (CS) and genetic predisposition [2,3]. COPD is characterized by emphysematous destruction of the alveoli and thickening of the small airway walls in response to chronic exposure to cigarette smoke (CS). While the cause of emphysematous destruction of the alveoli is likely due to a combination of the cytotoxic and pro-inflammatory activities of CS, the pathogenic mechanisms underlying the disease remain inadequately defined .
Recently, progress in genome-wide association studies (GWAS) has provided compelling evidence for several COPD susceptibility loci, including an intergenic region on chromosome 4q31 [5,6]. Subsequent work by our group confirmed hedgehog interacting protein (HHIP) as the causative gene for this locus and identified a potential functional variant . HHIP is a negative regulator of the Hedgehog pathway that is important for morphogenesis of the lungs and other organs [8-10]. HHIP competitively binds to all three ligands of the Hedgehog pathway: Sonic Hedgehog (Shh), Indian Hedgehog (Ihh), and Desert Hedgehog (Dhh). Recent reports have suggested potential roles for the Hedgehog pathway in bleomycin-induced pulmonary fibrosis  and CS-induced lung cancer . However, the mechanisms by which HHIP influences COPD susceptibility remain to be determined. We have now extended our previous work to an in vivo murine model by demonstrating that Hhip haploinsufficiency exaggerated CS-induced airspace enlargement and rewired lymphocyte activation signaling pathways under CS exposure.
Mice and cigarette smoke exposure in mice
Hhip +/- mice in a mixed genetic background (Jackson Laboratory) were backcrossed for 12 generations to C57BL/6J wild-type mice that were purchased from Jackson Laboratory. All mice were housed in the animal facility at Harvard Medical School with a 12 h light/12 h dark cycle. This study was performed in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health. The protocol was approved by the Harvard Medical Area (HMA) Standing Committee on Animals (Protocol Number: 04833). Every effort was made to minimize suffering.
Female Hhip +/- and littermate mice (approximately 10 weeks old) were exposed to mixed main-stream and side-stream CS from 3R4F Kentucky Research cigarettes for 5 days/week in Teague TE 10z chambers (Total Suspended Particulates approximately 100 to 200 mg/m3 and CO levels approximately 6 ppm). As a control, mice were exposed to filtered air for the same duration. At the end of the 6-month exposure period, both respiratory mechanics and airspace enlargement were measured in mice exposed to CS or air. All mice were euthanized by CO2 narcosis and cervical dislocation before removing lungs for other studies.
Measurement of lung mechanics
To measure respiratory mechanics, mice were anesthetized with a cocktail of 100 mg/kg body weight ketamine, 10 mg/kg xylazine, and 3 mg/kg acepromazine. A tracheostomy was performed, and an 18-gauge cannula was inserted and secured in the trachea using sutures. Mice were then connected via the cannula to a digitally controlled mechanical ventilator (Flexi Vent device; Scireq Inc., Montreal, QC, Canada). Ventilator settings were f = 150/min, FiO2 = 0.21, tidal volume = 10 mL/kg body weight, and positive end-expiratory pressure (PEEP) = 3 cm H2O. Murine lungs were inflated to total lung capacity (TLC; 25 cm H2O) three times for a volume history. Quasi-lung compliance and tissue elastance at PEEP of 3 cm H2O were then measured.
Morphometric analysis of airspace size in mice
After completion of the respiratory mechanics measurements, mice were euthanized; lungs were inflated with PBS to a constant pressure of 25 cm H2O and then fixed with 10% formaldehyde for 48 h. The mean alveolar chord length was measured after slides were stained with Gill’s solution as described previously . At least 15 images per mouse lung were randomly captured for analysis using methods described previously . Images were processed and analyzed for mean alveolar chord length using Scion image software.
Gene expression profiling
Six mice from each of four groups with different genotypes (Hhip +/+ or Hhip +/- ) and treatments (air or CS) were randomly chosen for gene expression profiling. RNA was extracted from murine lung tissues using Allprep kit (Qiagen, Valencia, CA, USA); 750 ng of total RNA for each sample was hybridized onto a Sentrix MouseRef-8 v2.0 Expression BeadChip Array (Illumina, San Diego, CA, USA), based on the manufacturer’s protocol. The BeadChip Array was scanned using an Illumina BeadArray Reader and quality control analysis was performed using Illumina GenomeStudio v3.1.3 software, Gene Expression Module version 1.1.1, and the statistical package R v2.15.2. Each of the samples and all probes passed quality control, and the data were background corrected, log2 transformed, and quantile normalized. Only probes with gene symbol and chromosome number annotations were retained, providing 22,121 of the 25,697 probes for analysis. A linear regression model was fit to the data for each probe to detect differential expression between genotypes (Hhip +/+ /Hhip +/- ) and treatments (air/CS), and having the processing information available for each array. The chip barcode was included as a covariate in the model to correct for batch effects. Clustering and MDS plots provided evidence of acceptable batch correction. In order to provide insight into probes demonstrating differential expression in Hhip +/+ and Hhip +/- mice in response to air/CS treatment, the interaction term genotype*treatment was included in the linear model. The genotype-by-treatment interaction was evaluated using the linear model expression ~ treatment + genotype + treatment*genotype + batch. Lastly, to stabilize the standard error of the estimated regression coefficient, we applied an empirical Bayes shrinkage method to obtain a moderated t-test statistic and its P value (R package limma v3.14.3). To adjust for multiple testing and to control false discovery rate (FDR), we corrected the P values using the Benjamin and Hochberg (1995) FDR correction method. Genes with P adjusted <0.05 were defined as significantly different gene expression. The network-based functional prediction was performed on the GeneMANIA webserver [15,16] using the default list of networks and weighting.
The microarray data have been deposited in NCBI’s Gene Expression Omnibus  and are accessible through GEO Series accession number GSE65124.
Detection of gene expression by real-time RT-PCR
Real-time PCR with gene-specific TaqMan primers/probes was performed as described previously . Relative gene expression was calculated based on the standard 2-ΔΔCT method, using GAPDH (glyceraldehyde 3-phosphate dehydrogenase) as a reference gene. Gene expression levels were analyzed by both two-way ANOVA analysis and student two-sample unpaired t test.
Quantification of lymphoid aggregates in murine lungs exposed chronically to cigarette smoke
The number of lymphoid aggregates was determined by counting total number of lymphoid aggregates around airways with 200 to 1,000 μm diameter in each murine lung section that was stained with hematoxylin and eosin. The correlation between the number of lymphoid aggregates per airway and airspace size represented by mean airspace chord length (MCL) was determined by Pearson’s correlation coefficient using commercial statistical software STATA (version 12.1). A scatter plot analysis of MCL vs. lymphoid aggregates per airway was performed in the four groups of mice.
Immunohistochemistry staining of T cells and B cells in lung sections
We performed triple immunofluorescence staining of lung sections from Hhip+/- mice exposed to CS for 6 months. Briefly, slides were deparaffinized, and antigen retrieval was achieved by heating slides in 10 mM sodium citrate and 2 mM citric acid buffer (pH 6.0). Lung sections were incubated for 1 h at 37°C with goat anti-CD4 IgG (Abcam, Cambridge, MA, USA; diluted 1:50), then for 1 h at 37°C with rat anti-CD8 IgG (Novus Biologicals, Littleton, CO, USA; diluted 1:100), and then 1 h at 37°C with rabbit anti-CD20 IgG (Abcam, Cambridge, MA, USA; diluted 1:100) followed by incubation with second antibodies including Alexa 647-conjugated rabbit anti-goat IgG (diluted 1:100), Alexa 488-conjugated goat anti-rat F(ab')2 (diluted 1:100), and Alexa 546-conjugated goat anti-rabbit IgG (diluted 1:100). Images of the stained lung sections were analyzed with a confocal microscope (Leica Microsystems, Buffalo Grove, IL, USA).
Quantification of CD8+ T cells in lung parenchyma
The number of CD8+ T cells was measured in lung sections from Hhip +/+ and Hhip +/- mice exposed to air (n = 4 mice/group) or CS (n = 5 Hhip +/+ and n = 9 Hhip +/- mice) for 6 months. After antigen retrieval in 10 mM sodium citrate and 2 mM citric acid buffer (pH 6.0), lung sections were incubated with a rat IgG to murine CD8 (Novus Biologicals, Littleton, CO, USA; diluted 1:100) or non-immune murine IgG as control for 1 h at 37°C. After secondary antibody incubation and counterstaining with DAPI, the number of CD8 positive cells was then normalized by unit of alveolar area that was quantified using MetaMorph software (Molecular Devices, Wetzlar, Germany) and compared between groups using two-way ANOVA analysis and two-sided unpaired t-test.
Measurement of cytokines by enzyme-linked immunosorbent assay (ELISA)
Interleukin 2 (IL-2) and Interferon gamma (IFN-γ) levels in murine lungs were measured using mouse Quantikine ELISA Kit and mouse IFN-γ Duoset (R&D, Minneapolis, MN, USA) following the manufacturer’s protocol and then normalized to total lung protein amount measured using Bradford Assay (Biorad, Hercules, CA, USA).
Leukocyte counting in bronchoalveolar lavage samples after acute smoke exposure
Hhip +/+ and Hhip +/- mice were exposed to CS or air for 2 months and then bronchoalveolar lavage was performed using 3 x 1 mL aliquots of PBS. Total leukocyte counts were performed using a hemocytometer.
Flow cytometric analysis to quantify activation of T lymphocytes
Immediately after CS exposure and lung removal, lung tissues were digested with collagenase and DNase (Sigma, St. Louis, MO, USA) to prepare single-cell suspensions by passing the dissociated tissue through a 70-μm cell strainer (BD Falcon). The suspensions of cells were stained in PBS containing 2% FBS with the following antibodies: FITC-conjugated anti-CD4, PE-Cy7-conjugated anti-CD3, APC-Cy7-conjugated anti-CD45, PE-conjugated anti-CD69, and APC-conjugated anti-CD8 antibodies (BD Bioscience). The cells were gated on FSC (forward scatter) and SSC (side scatter). The cells were then gated on SSC versus CD45 followed by SSC versus CD3 in order to gate on all T cells. The cells were further gated to obtain the percentage of CD4+ T and CD8+ T cells. Stained cells were examined on a FACS Canto flow cytometer (BD Bioscience, San Jose, CA, USA) and analyzed using FlowJo software (TreeStar, Ashland, OR, USA). The percentages of activated CD4+ T cells and activated CD8+ T cells were determined based on double positive cells (CD69+ CD4+ or CD69+ CD8+ cells) among total CD4+ or CD8+ T cells.
We used the ‘PANDA’ method  to construct four genome-wide regulatory network models based on Hhip genotype (Hhip +/- or Hhip +/+ ) and treatment (CS or Air). To begin, we collapsed expression values for probe sets that map to the same gene symbol by selecting, for each experiment, the signal from the probe set with the most significant detection P value, resulting in expression information for 15,124 unique genes. We then obtained the position-weight-matrix for 374 unique vertebrate TFs with known binding site motifs, 130 recorded in the Jaspar database , and 290 from the uniprobe database  (46 TFs have motifs in both databases). Defining a ‘true hit’ as a log-odds ratio greater than 7.5, we used HOMER  to map each of these motifs to putative regulatory regions, defined as (-750, +250) base-pairs around the transcription start site, for 15,124 genes with expression data in the samples. Because transcriptional regulation involves assembly of protein complexes and allows for combinatorial regulatory processes, we also collected information regarding physical protein interaction data between transcription factors estimated from mouse-2-hybrid analysis . Finally, we used PANDA  to integrate information from transcription factor binding motifs, protein-protein interaction data, and gene expression for each experimental condition, and constructed directed networks for each experimental condition. To reduce potential background noise due to microarray sample size (6 mice/group), we applied a jack-knife  method to repeatedly sample four arrays out of six from each experimental condition (Hhip +/+ -Air or CS and Hhip +/- Air or CS), building 15 networks per group.
After reconstructing the networks, we specifically investigated the differential targeting of HHIP in the predicted networks. To begin, we identified 32 transcription factors that have a putative target in the promoter region of HHIP. We then identified the 32 ‘edges’ that extend between these transcription factors and HHIP in our reconstructed network models. To evaluate whether any of these 32 edges are significantly different between the sets of networks associated with each of the experimental conditions, we performed an unpaired two-sample t-test, comparing the distribution of edge weights in each pair of conditions.
We next identified subnetworks of edges that are most distinct between each pair of our reconstructed network models. For each edge connecting a transcription factor to its target gene, PANDA assigns a Z-score weight reflecting the confidence level of the inferred regulatory relationship. We contrasted both the average and the distribution (using an unpaired two sample t-test) of this weight value between sets of networks (those reconstructed for each experimental condition using the jack-knife procedure, see above). Using this information, we identified subnetworks of edges that differ between each pair of networks by selecting high confidence edges, defined as those with an absolute difference in average weight greater than 3 and a P value from the t-test of less than 1×10-3. We used DAVID  to test for pathway enrichment in genes targeted in each of these identified subnetworks. The P values for KEGG Pathways enriched at a Benjamini-Hochberg significance of less than 0.01 in at least one of these subnetworks were illustrated in a heat map. At this confidence level the selected pathways all had an overlap of greater than 10 genes with their respective differentially targeted gene sets.
We also characterized the differential targeting patterns of transcription factors in these pairs of subnetworks by calculating, for each transcription factor, a fold-enrichment for its number of outgoing edges in one subnetwork compared to the other, as well as an associated P value. To calculate the P value we performed Fisher’s exact test and determined the significance of the overlap in the transcription factor’s outgoing edges with the edges in its associated ‘enriched’ subnetwork, using edges in either subnetwork as a background. A list of core TFs for each of the subnetworks was then generated by filtering TFs with total edges > =3 and absolute log2 fold change greater than 1 or -log10 P value greater than 3 in either subnetwork. Finally, we took the union of these core TFs and determined their average expression across samples in each of the four experimental groups. We also used an unpaired t-test to evaluate their pairwise differential expression between the groups. Note that several of these core TFs identified through the network analysis (which utilizes motif information) did not have probes on the array, including ELK1, SOX11, and ZFP740, and these are excluded from expression analysis.
Statistical analysis methods
For MCL analysis, we applied a linear mixed effects model for mice that were exposed to CS, in which the probit-transformed lung damage (quantified alveolar chord length for each image from each mouse) is the outcome variable, the intercept is the random effect to account for the dependence among the longitudinal measurements of the lung damage measurement, and mice type is the fixed effect (wild-type mice as the reference group). We used two-way ANOVA analysis first to determine whether there are overall differences among all four groups of mice followed by unpaired t tests to determine whether there are significant differences between each pair of groups. For gene expression of Gstp1, Cstw, Fyb, Lat2, and Klf4, we first evaluated if there is an interaction effect between genotype and treatment on gene expression using two-way ANOVA analysis. We then compared the two conditions of a factor (genotype or CS treatment) for a given level of the other factor (CS treatment or genotype) using un-paired two sample two-sided t tests. However, for statistical analysis on CD8+ T cells quantification, to handle the non-normality due to outliers we applied Kruskal-Wallis test, which is a non-parametric test without assumption of normality, and detected significant differences among groups (P <0.001). We then applied Wilcoxon rank sum test for pair-wise comparisons. All other data were also analyzed by both two-way ANOVA analysis and student two-sample unpaired t tests.
Hhip +/- mice demonstrated more severe airspace enlargement induced by cigarette smoke compared with wild-type littermate control mice
For lung mechanics, CS-exposed Hhip +/- mice had increased quasi-static lung compliance (Figure 1C) and correspondingly decreased tissue elastance (Figure 1D) compared to CS-exposed Hhip +/+ mice, indicating that Hhip +/- mice exposed to CS have more compliant lungs and a greater loss of elastic recoil than Hhip +/+ mice exposed to CS.
Assessments of the Hedgehog pathway in murine lungs
Given that HHIP is known to inhibit the Hedgehog pathway, we assessed the expression of major components in the Hedgehog pathway in the chronic CS exposure murine model. Interestingly, Gli1 and Patch1 showed increased expression in air-exposed Hhip +/- mice vs. air-exposed Hhip +/+ mice while Gli1 and Gli2 showed increased expression in CS-exposed Hhip +/- mice vs. CS-exposed Hhip +/+ mice. No significant differences in Gli3 were observed among these four groups of mice (Additional file 1: Figure S2). Hence, Hhip +/- mice may have relatively higher activity of the Hedgehog pathway at baseline independent of CS treatment.
Gene expression profiling in murine lungs after cigarette smoke exposure
In order to identify molecular targets that contributed to increased airspace size in Hhip +/- mice exposed to CS, we performed gene expression profiling in lung tissues from randomly chosen Hhip +/- and Hhip +/+ mice exposed to either air or CS for 6 months (n = 6 mice/group). Data were preprocessed using the LIMMA package in R/Bioconductor, as we have done previously .
We then validated expression changes of 16 genes that demonstrated significant gene-by-CS treatment interactions and also were related to lymphocyte activation pathways by real-time RT-PCR (Additional file 2: Table S4). One gene, Gstp1 (Glutathione S-Transferase P1), which showed differential expression between both Hhip +/- -CS vs. Hhip +/+ -CS (green in Figure 2A) and Hhip +/+ -CS vs. Hhip +/+ -air (yellow in Figure 2A), demonstrated diminished induction by CS in Hhip +/- vs. Hhip +/+ mice (Figure 2B). Moreover, several genes that were identified as having significant genotype-by-treatment interactions and were related to lymphocyte activation showed increased expression in the lungs of Hhip +/- -CS compared to Hhip +/+ -CS mice (Figure 2C) and significant genotype-by-treatment interaction as suggested by two-way ANOVA analysis (Figure 2C, P <0.01). For example, Fyb (Fyn binding protein), also called ADAP (adhesion and degranulation-promoting adapter protein), promotes T cell conjugation to antigen-laden antigen presenting cells (APC) and enhances T cell survival . Linker for activation of T cells family member 2 (Lat2) showed significantly increased expression in the lungs of Hhip +/- -CS compared to Hhip +/+ -CS mice, in contrast to slightly reduced expression in Hhip +/+ mice after CS exposure. Ctsw (Cathepsin W), encoding a papain-like cysteine proteinase that is specifically expressed in natural killer cells (NK) and cytotoxic T cells, showed increased expression in CS-exposed Hhip +/- mice compared to Hhip +/+ mice.
HHIP haploinsufficiency led to activation of CD8+ T cells
Furthermore, a T lymphocyte-related cytokine, IL-2 , was also significantly increased in lung tissues of Hhip +/- mice compared to Hhip +/+ mice, despite the lack of CS-induced increase in IL-2 levels in Hhip +/+ mice (Additional file 1: Figure S5B, two-way ANOVA analysis, P <0.05). Interferon gamma (IFN-γ) also showed a trend toward increased levels in the lungs of Hhip +/- mice compared to Hhip +/+ mice after CS exposure (Additional file 1: Figure S5C). These changes in lung cytokine levels suggested a genotype-dependent response to CS in Hhip +/- mice compared to Hhip +/+ mice.
In addition to chronic CS exposure, we also compared lung inflammatory responses in bronchoalvoelar lavage (BAL) samples from Hhip +/- versus Hhip +/+ mice exposed to CS for 2 months. As expected, total cells in BAL samples increased in Hhip +/+ and Hhip +/- mice exposed to 2 months of CS (Figure 3E) but total BAL cell counts were comparable in CS-exposed Hhip +/- and CS-exposed Hhip +/+ mice. More importantly, the percentage of activated cytotoxic T cells, represented by positive staining for both CD8 and CD69, significantly increased in enzymatically digested lung samples from CS-exposed Hhip +/- mice compared with those in CS-exposed Hhip +/+ mice as detected by flow cytometry (Figure 3F). In contrast, the percentage of activated CD4+ T cells increased in the lungs after CS exposure in both Hhip +/+ and Hhip +/- mice, and CS-exposed Hhip +/- mice had a similar percentage of activated CD4+ T cells in their lungs as CS-exposed Hhip +/+ mice (Additional file 1: Figure S5D), indicating that HHIP haploinsufficiency had less effect on the CS-induced activation of CD4+ T cells in murine lungs.
Application of network models in mice exposed chronically to cigarette smoke
To interrogate potential regulatory mechanisms that differed in Hhip +/- vs. Hhip +/+ mice and infer molecular driving factors that may contribute to the murine phenotypes described above, we constructed network models from microarray gene expression data using the PANDA (Passing Attributes between Networks for Data Assimilation) method . PANDA integrates multiple sources of information, including protein-protein interaction, gene expression, and sequence motif data, to construct genome-wide, condition-specific regulatory networks. PANDA starts with a ‘prior’ map of potential transcription factor (TF)-gene interactions, often estimated by scanning the promoter regions of genes for the binding sites of TFs with known position-weight matrices. It then integrates this information with gene expression and protein-protein interaction data, using a message-passing approach to iteratively update the given initial map by evaluating (1) the similarity between gene co-expression and TF-targeting patterns as well as (2) the similarity between TF-TF (protein) interactions and input gene regulatory patterns. The result is a regulatory network estimate consistent with all sources of input data.
In PANDA networks, HHIP is predicted to be connected/regulated by 32 TFs with a putative binding site in the promoter of HHIP (Additional file 1: Figure S6). We then assessed significantly differentially targeted genes in each identified subnetwork by KEGG pathway analysis. The T cell receptor signaling pathways and chemokine signaling pathways were identified to be differentially regulated when comparing Hhip +/- vs. Hhip +/+ mice either exposed to air or CS (indicated by blue stars in Figure 4C).
Recent genome-wide association studies (GWAS) of COPD have identified several genomic regions that are clearly associated with COPD susceptibility; however, the biological mechanisms by which these genetic loci influence COPD remain undefined. Within one of the most well established COPD GWAS loci upstream from the HHIP gene on chromosome 4q31, we have previously identified a functional variant that exerts long-range regulation of the HHIP gene . We have now extended those studies by exposing Hhip +/- mice to chronic CS to model the biological impact of HHIP haploinsufficiency in CS-induced emphysema. We demonstrated that HHIP haploinsufficiency exaggerated CS-induced airspace enlargement. Gene expression profiling was used to investigate the molecular basis for the observed lung morphological and functional changes in Hhip haploinsufficient mice with chronic smoke exposure. Using PANDA, a recently developed network modeling approach, we characterized transcriptional regulatory networks that contribute to the CS-induced alterations in the Hhip +/- murine lung. Among 15 identified TFs that control differential regulation in Hhip +/- -CS mice compared with Hhip +/+ -CS mice, the top highly-expressed TFs in murine lungs include Klf4, which also showed reduced expression in lungs of CS-exposed Hhip +/- mice compared to Hhip +/+ mice.
Functional interpretation of GWAS candidate regions/genes in complex diseases has been challenging in human genetics . Given the early post-natal lethality in homozygous Hhip knockout mice and the moderate effects of the previously identified human functional variant on expression of HHIP , we evaluated a heterozygous murine model to dissect the function of this gene in a disease-relevant murine model. It is not unprecedented to observe significant phenotypic changes in heterozygous mice with deficiency of genes that are crucial for morphogenesis, such as BECN1 , TSC2 , and NPM1  - especially when tissue injury is induced in mice.
Increased lymphoid aggregates in human COPD lung samples were reported 10 years ago to be associated with progression of COPD . We also observed similar lymphoid aggregates formed around airways in CS-exposed mice, and the number of these lymphoid aggregates significantly increased in Hhip +/- mice exposed to CS. The presence of lymphoid aggregates was highly correlated with airspace enlargement in Hhip+/- mice exposed to CS (Figure 3C), suggesting dysregulated immune responses in Hhip +/- mice associated with the severity of airspace enlargement. Furthermore, subsequent FACS analysis demonstrated increased CD8+ T cell activation in Hhip +/- mice exposed to 2 months of CS. Thus, CD8+ T cells may represent the major type of activated lymphocytes in Hhip +/- mice exposed to CS, which is consistent with findings in human COPD lungs . However, more in-depth investigations are needed to assess the involvement of B cell signaling pathways, as suggested by the pathway analysis based on gene-by-treatment interaction genes (Additional file 2: Table S2B).
In addition to airspace enlargement, we also observed increased lung compliance and reduced tissue elastance in Hhip+/- mice exposed to CS (Figure 1B). However, we observed no significant changes of lung compliance and tissue elastance in wild type mice exposed to 6 months of chronic CS, similar to a previous study that failed to detect any significant differences in lung compliance between C57BL/6 wild-type mice exposed to air and smoke for 6 months . The most likely reason for our inability to detect a CS-induced change in lung compliance in wild-type mice is that the emphysema that develops in this murine strain is relatively mild in severity. However, significant increases in lung compliance that we observed in CS-exposed Hhip +/- mice when compared with CS-exposed Hhip +/+ mice are consistent with the increase in distal airspace size that we detected in CS-exposed Hhip +/- mice. Thus, lung function changes in CS-exposed Hhip +/- mice strongly supported that haploinsufficiency of Hhip increased the severity of CS-induced emphysema that develops in C57BL/6 strain mice.
Recent studies in murine models [36-38] and human samples  strongly supported pathological roles of activation of T cells in CS-induced emphysema development. Although similar contributions of CD4+ T helper or CD8+ T cytotoxic cells to airspace enlargement were observed in various CS-induced emphysema murine models [36-38], we observed significantly enhanced CD8+ T cell activation in Hhip +/- mice that can produce deleterious effects leading to emphysema progression in both humans [27,40] and mice  due to recruitment of NK cells. In addition, activated CD8+ T cells in CS-exposed lungs release cytotoxic proteins, including perforin and granzyme, that promote cellular protein degradation and cleavage, leading to apoptotic alveolar cell death in CS-exposed lungs [40-42], thus contributing to emphysematous destruction in lungs.
Our network analysis, based on the PANDA method, provided novel gene regulatory insights into Hhip haploinsufficiency-imposed signaling rewiring in a systematic way. Among TFs that uniquely regulate the pulmonary response of Hhip +/- mice to CS, Klf4 not only showed differential connectivity but also was significantly reduced in expression in Hhip +/- murine lungs. Klf4 is important in maintaining T lymphocyte homeostasis as: (1) reduced expression of Klf4 is required for lineage commitment of T lymphocytes ; and (2) Klf4 directly regulates thymocyte proliferation . More interestingly, the proliferation and homing of CD8+ T cells is controlled through a transcriptional regulation axis of Elf4-Klf4 and Elf4-Klf2 , which provides one plausible mechanism by which HHIP haploinsufficiency modulates lymphocyte activation through Klf4 signaling. This altered Klf4 signaling may be caused by either reduced expression of Klf4 in Hhip +/- mice (Additional file 1: Figure S7) or differential targets of Klf4 in Hhip +/- mice compared with Hhip +/+ mice in response to chronic CS exposure (Figure 5A). Additionally, a recent report showed that activation of the Hedgehog pathway in CD8+ T cells controlled killing function of cytotoxic T cells by polarizing the cytoskeleton to facilitate delivery of cytotoxic granules to the targets of cytotoxic T cells . This also suggests that the Hedgehog pathway may contribute to the enhanced CD8+ T cell function in Hhip +/- mice through a similar mechanism .
Together, our results indicate that after chronic CS exposure, Hhip haploinsufficient mice developed increased emphysema that was related to increased CD8+ T lymphocyte activation and differential Klf4 signaling in Hhip +/- mice. In the future, additional comprehensive immunological experiments are needed to determine whether and how HHIP regulates lymphocyte proliferation, differentiation, and/or cytotoxic killing through Klf4 signaling and/or the Hedgehog pathway. Furthermore, whether Klf4 signaling and the Hedgehog pathway cooperate to modulate activation of naïve CD8+ T cells needs additional investigation.
There are several limitations to our study. First, the HHIP gene locus has also been associated with lung function in general population samples [46,47]. However, in our current system, no significant differences in lung morphology (Figure 1) between Air-exposed Hhip +/- -mice and Hhip +/+ -mice were observed. This may be linked to the relatively modest reduction in HHIP expression in Hhip +/- mice which may have been insufficient to result in significant lung function changes at the single time point studied here (8 months of age). Additional investigations are needed to explore how HHIP regulates lung function in non-smokers and whether there are common mechanisms shared between HHIP locus-regulated lung function and susceptibility to CS-induced emphysema. Second, transcription factors in the Hedgehog pathway, Gli1, Gli2, and Gli3, were not included in the murine motif database when we built PANDA networks. As a larger number of TFs are included in the existing motif database, our future network modeling studies using the PANDA method may reveal additional insights into the activities of Hhip in regulating CS-induced lung effects in mice. Third, we also attempted to compare gene expression changes caused by Hhip haploinsufficiency in murine lungs (Hhip +/- -air vs. Hhip +/+ -air) with our previous microarray analysis in an airway epithelial cell line (Beas-2B cells), treated with either non-targeting shRNA controls or HHIP shRNAs . No overlapping genes were identified (data not shown). This could be due to sample type differences (whole murine lung vs. a single human cell line), differential reductions in HHIP expression (approximately 33% vs. >70% reductions at mRNA level, respectively), or timing of sample collection (total lung lysates from mice at 8 months of age vs. cultured human cell lines). Lastly, we did not observe significantly increased lung levels of IL-2 and IFN-γ in wild type C57BL/6 mice exposed to 6 months of CS. Our results are consistent with a previous report showing that the C57BL/6 murine strain was a mildly susceptible strain with moderate emphysema development without any increases in lung IFN-γ levels after 6 months of CS exposure . However, under the same CS exposure condition, we observed a trend toward increased IFN-γ in Hhip+/- mice, suggesting the possibility of genotype-specific CD8+ T activation after CS exposure.
In summary, combined with our previous finding that a COPD risk allele was associated with reduced enhancer activity for HHIP , our murine model has: (1) successfully demonstrated that modest reductions in Hhip gene expression, similar to smokers carrying HHIP COPD GWAS risk-alleles, led to increased susceptibility to CS-induced emphysema in Hhip +/- mice; (2) shown that HHIP modulates lymphocyte activation especially in CD8+ T cells; (3) characterized networks with genes and TFs that are differentially regulated in murine lungs responsive to either CS treatment or haploinsufficiency of Hhip; and (4) revealed roles of Hhip in maintaining lung homeostasis that is disrupted in the haploinsufficiency state in the adult murine lung in addition to its known roles in embryonic lung development (identified using a loss-of-function approach) . Identification of strategies to correct alterations in network rewiring in Hhip +/- mice could eventually facilitate the development of novel therapies for human COPD.
We sincerely thank Dr. Nandini Krishnamoorthy for her generous help with the flow cytometry data analysis on lung samples from mice exposed to 2 months of cigarette smoke.
This work was supported by U.S. National Institutes of Health (NIH) grants R01 HL075478 and P01 HL105339 (EKS), R01HL111759 (EKS, JQ, GCY), R21HL120794 (EKS and XZ), R01 AI111475-01 and R21 HL111 835 (CO), the BWH-LRRI Research Consortium (CO), and the Flight Attendants Medical Research Institute grant CIA123046 (CO).
- Arialdi M, Miniño MJX, Kenneth D, Kochanek MA. Division of vital statistics. Deaths: preliminary data for 2008. Natl Vital Stat Rep. 2010;59:1–72.Google Scholar
- McCloskey SC, Patel BD, Hinchliffe SJ, Reid ED, Wareham NJ, Lomas DA. Siblings of patients with severe chronic obstructive pulmonary disease have a significant risk of airflow obstruction. Am J Respir Crit Care Med. 2001;164:1419–24.View ArticlePubMedGoogle Scholar
- Silverman EK, Chapman HA, Drazen JM, Weiss ST, Rosner B, Campbell EJ, et al. Genetic epidemiology of severe, early-onset chronic obstructive pulmonary disease, Risk to relatives for airflow obstruction and chronic bronchitis. Am J Respir Crit Care Med. 1998;157:1770–8.View ArticlePubMedGoogle Scholar
- Salazar LM, Herrera AM. Fibrotic response of tissue remodeling in COPD. Lung. 2011;189:101–9.View ArticlePubMedGoogle Scholar
- Pillai SG, Ge D, Zhu G, Kong X, Shianna KV, Need AC, et al. A genome-wide association study in chronic obstructive pulmonary disease (COPD): identification of two major susceptibility loci. PLoS Genet. 2009;5:e1000421.View ArticlePubMed CentralPubMedGoogle Scholar
- Sovio U, Bennett AJ, Millwood IY, Molitor J, O’Reilly PF, Timpson NJ, et al. Genetic determinants of height growth assessed longitudinally from infancy to adulthood in the northern Finland birth cohort 1966. PLoS Genet. 2009;5:e1000409.View ArticlePubMed CentralPubMedGoogle Scholar
- Zhou X, Baron RM, Hardin M, Cho MH, Zielinski J, Hawrylkiewicz I, et al. Identification of a chronic obstructive pulmonary disease genetic determinant that regulates HHIP. Hum Mol Genet. 2012;21:1325–35.View ArticlePubMed CentralPubMedGoogle Scholar
- Chuang PT, McMahon AP. Vertebrate Hedgehog signalling modulated by induction of a Hedgehog-binding protein. Nature. 1999;397:617–21.View ArticlePubMedGoogle Scholar
- Kawahira H, Ma NH, Tzanakakis ES, McMahon AP, Chuang PT, Hebrok M. Combined activities of hedgehog signaling inhibitors regulate pancreas development. Development. 2003;130:4871–9.View ArticlePubMedGoogle Scholar
- Chuang PT, Kawcak T, McMahon AP. Feedback control of mammalian Hedgehog signaling by the Hedgehog-binding protein, Hip1, modulates Fgf signaling during branching morphogenesis of the lung. Genes Dev. 2003;17:342–7.View ArticlePubMed CentralPubMedGoogle Scholar
- Liu L, Kugler MC, Loomis CA, Samdani R, Zhao Z, Chen GJ, et al. Hedgehog signaling in neonatal and adult lung. Am J Respir Cell Mol Biol. 2013;48:703–10.View ArticlePubMed CentralPubMedGoogle Scholar
- Lemjabbar-Alaoui H, Dasari V, Sidhu SS, Mengistab A, Finkbeiner W, Gallup M, et al. Wnt and Hedgehog are critical mediators of cigarette smoke-induced lung cancer. PLoS One. 2006;1:e93.View ArticlePubMed CentralPubMedGoogle Scholar
- Moghadaszadeh B, Rider BE, Lawlor MW, Childers MK, Grange RW, Gupta K, et al. Selenoprotein N deficiency in mice is associated with abnormal lung development. FASEB J. 2013;27:1585–99.View ArticlePubMed CentralPubMedGoogle Scholar
- An CH, Wang XM, Lam HC, Ifedigbo E, Washko GR, Ryter SW, et al. TLR4 Deficiency Promotes Autophagy during Cigarette Smoke-Induced Pulmonary Emphysema. Am J Physiol Lung Cell Mol Physiol. 2012;303:L748–57Google Scholar
- Mostafavi S, Ray D, Warde-Farley D, Grouios C, Morris Q. GeneMANIA: a real-time multiple association network integration algorithm for predicting gene function. Genome Biol. 2008;9 Suppl 1:S4.Google Scholar
- GeneMANIA. [http://www.genemania.org/]
- Edgar R, Domrachev M, Lash AE. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002;30:207–10.View ArticlePubMed CentralPubMedGoogle Scholar
- Zhou X, Qiu W, Sathirapongsasuti JF, Cho MH, Mancini JD, Lao T, et al. Gene expression analysis uncovers novel hedgehog interacting protein (HHIP) effects in human bronchial epithelial cells. Genomics. 2013;101:263–72.View ArticlePubMed CentralPubMedGoogle Scholar
- Glass K, Huttenhower C, Quackenbush J, Yuan GC. Passing Messages Between Networks to Refine Predicted Interactions. PLoS One. 2013;8:e64832.View ArticlePubMed CentralPubMedGoogle Scholar
- Sandelin A, Alkema W, Engstrom P, Wasserman WW, Lenhard B. JASPAR: an open-access database for eukaryotic transcription factor binding profiles. Nucleic Acids Res. 2004;32:D91–4.View ArticlePubMed CentralPubMedGoogle Scholar
- Robasky K, Bulyk ML. UniPROBE, update 2011: expanded content and search tools in the online database of protein-binding microarray data on protein-DNA interactions. Nucleic Acids Res. 2011;39:D124–8.View ArticlePubMed CentralPubMedGoogle Scholar
- Heinz S, Benner C, Spann N, Bertolino E, Lin YC, Laslo P, et al. Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities. Mol Cell. 2010;38:576–89.View ArticlePubMed CentralPubMedGoogle Scholar
- Ravasi T, Suzuki H, Cannistraci CV, Katayama S, Bajic VB, Tan K, et al. An atlas of combinatorial transcriptional regulation in mouse and man. Cell. 2010;140:744–52.View ArticlePubMedGoogle Scholar
- Wu CFJ. Jackknife, bootstrap and other resampling methods in regression analysis. Annals Stat. 1986;14:1261–95.View ArticleGoogle Scholar
- da Huang W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57.View ArticleGoogle Scholar
- Mueller KL, Thomas MS, Burbach BJ, Peterson EJ, Shimizu Y. Adhesion and degranulation-promoting adapter protein (ADAP) positively regulates T cell sensitivity to antigen and T cell survival. J Immunol. 2007;179:3559–69.View ArticlePubMedGoogle Scholar
- Hogg JC, Chu F, Utokaparch S, Woods R, Elliott WM, Buzatu L, et al. The nature of small-airway obstruction in chronic obstructive pulmonary disease. N Engl J Med. 2004;350:2645–53.View ArticlePubMedGoogle Scholar
- Maeno T, Houghton AM, Quintero PA, Grumelli S, Owen CA, Shapiro SD. CD8+ T Cells are required for inflammation and destruction in cigarette smoke-induced emphysema in mice. J Immunol. 2007;178:8090–6.View ArticlePubMedGoogle Scholar
- D’Armiento JM, Scharf SM, Roth MD, Connett JE, Ghio A, Sternberg D, et al. Eosinophil and T cell markers predict functional decline in COPD patients. Respir Res. 2009;10:113.View ArticlePubMed CentralPubMedGoogle Scholar
- An J, Golech S, Klaewsongkram J, Zhang Y, Subedi K, Huston GE, et al. Kruppel-like factor 4 (KLF4) directly regulates proliferation in thymocyte development and IL-17 expression during Th17 differentiation. FASEB J. 2011;25:3634–45.View ArticlePubMed CentralPubMedGoogle Scholar
- Freedman ML, Monteiro AN, Gayther SA, Coetzee GA, Risch A, Plass C, et al. Principles for the post-GWAS functional characterization of cancer risk loci. Nat Genet. 2011;43:513–8.View ArticlePubMed CentralPubMedGoogle Scholar
- Qu X, Yu J, Bhagat G, Furuya N, Hibshoosh H, Troxel A, et al. Promotion of tumorigenesis by heterozygous disruption of the beclin 1 autophagy gene. J Clin Invest. 2003;112:1809–20.View ArticlePubMed CentralPubMedGoogle Scholar
- Ehninger D, Han S, Shilyansky C, Zhou Y, Li W, Kwiatkowski DJ, et al. Reversal of learning deficits in a Tsc2+/− mouse model of tuberous sclerosis. Nat Med. 2008;14:843–8.View ArticlePubMed CentralPubMedGoogle Scholar
- Grisendi S, Bernardi R, Rossi M, Cheng K, Khandker L, Manova K, et al. Role of nucleophosmin in embryonic development and tumorigenesis. Nature. 2005;437:147–53.View ArticlePubMedGoogle Scholar
- Guerassimov A, Hoshino Y, Takubo Y, Turcotte A, Yamamoto M, Ghezzo H, et al. The development of emphysema in cigarette smoke-exposed mice is strain dependent. Am J Respir Crit Care Med. 2004;170:974–80.View ArticlePubMedGoogle Scholar
- Podolin PL, Foley JP, Carpenter DC, Bolognese BJ, Logan GA, Long 3rd E, et al. T cell depletion protects against alveolar destruction due to chronic cigarette smoke exposure in mice. Am J Physiol Lung Cell Mol Physiol. 2013;304:L312–23.View ArticlePubMedGoogle Scholar
- Eppert BL, Wortham BW, Flury JL, Borchers MT. Functional characterization of T cell populations in a mouse model of chronic obstructive pulmonary disease. J Immunol. 2013;190:1331–40.View ArticlePubMed CentralPubMedGoogle Scholar
- Motz GT, Eppert BL, Wesselkamper SC, Flury JL, Borchers MT. Chronic cigarette smoke exposure generates pathogenic T cells capable of driving COPD-like disease in Rag2−/− mice. Am J Respir Crit Care Med. 2010;181:1223–33.View ArticlePubMed CentralPubMedGoogle Scholar
- Wang J, Urbanowicz RA, Tighe PJ, Todd I, Corne JM, Fairclough LC. Differential activation of killer cells in the circulation and the lung: a study of current smoking status and chronic obstructive pulmonary disease (COPD). PLoS One. 2013;8:e58556.View ArticlePubMed CentralPubMedGoogle Scholar
- Chrysofakis G, Tzanakis N, Kyriakoy D, Tsoumakidou M, Tsiligianni I, Klimathianaki M, et al. Perforin expression and cytotoxic activity of sputum CD8+ lymphocytes in patients with COPD. Chest. 2004;125:71–6.View ArticlePubMedGoogle Scholar
- Freeman CM, Han MK, Martinez FJ, Murray S, Liu LX, Chensue SW, et al. Cytotoxic potential of lung CD8(+) T cells increases with chronic obstructive pulmonary disease severity and with in vitro stimulation by IL-18 or IL-15. J Immunol. 2010;184:6504–13.View ArticlePubMed CentralPubMedGoogle Scholar
- Shiratsuchi N, Asai K, Kanazawa H, Kyoh S, Tochino Y, Kodama T, et al. Measurement of soluble perforin, a marker of CD8+ T lymphocyte activation in epithelial lining fluid. Respir Med. 2011;105:1885–90.View ArticlePubMedGoogle Scholar
- Wen X, Liu H, Xiao G, Liu X. Downregulation of the transcription factor KLF4 is required for the lineage commitment of T cells. Cell Res. 2011;21:1701–10.View ArticlePubMed CentralPubMedGoogle Scholar
- Yamada T, Park CS, Mamonkin M, Lacorazza HD. Transcription factor ELF4 controls the proliferation and homing of CD8+ T cells via the Kruppel-like factors KLF4 and KLF2. Nat Immunol. 2009;10:618–26.View ArticlePubMed CentralPubMedGoogle Scholar
- de la Roche M, Ritter AT, Angus KL, Dinsmore C, Earnshaw CH, Reiter JF, et al. Hedgehog signaling controls T cell killing at the immunological synapse. Science. 2013;342:1247–50.View ArticlePubMed CentralPubMedGoogle Scholar
- Hancock DB, Eijgelsheim M, Wilk JB, Gharib SA, Loehr LR, Marciante KD, et al. Meta-analyses of genome-wide association studies identify multiple loci associated with pulmonary function. Nat Genet. 2010;42:45–52.View ArticlePubMed CentralPubMedGoogle Scholar
- Repapi E, Sayers I, Wain LV, Burton PR, Johnson T, Obeidat M, et al. Genome-wide association study identifies five loci associated with lung function. Nat Genet. 2010;42:36–44.View ArticlePubMed CentralPubMedGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.