A gene expression signature of emphysema-related lung destruction and its reversal by the tripeptide GHK
- Joshua D Campbell1, 2,
- John E McDonough3,
- Julie E Zeskind1, 2,
- Tillie L Hackett3,
- Dmitri V Pechkovsky3,
- Corry-Anke Brandsma4,
- Masaru Suzuki3,
- John V Gosselink3,
- Gang Liu1,
- Yuriy O Alekseyev5,
- Ji Xiao1,
- Xiaohui Zhang1,
- Shizu Hayashi3,
- Joel D Cooper6,
- Wim Timens4,
- Dirkje S Postma7,
- Darryl A Knight3,
- Marc E Lenburg1, 2Email author,
- James C Hogg†3 and
- Avrum Spira1, 2Email author
© Campbell et al.; licensee BioMed Central Ltd. 2012
Received: 5 March 2012
Accepted: 16 August 2012
Published: 31 August 2012
Chronic obstructive pulmonary disease (COPD) is a heterogeneous disease consisting of emphysema, small airway obstruction, and/or chronic bronchitis that results in significant loss of lung function over time.
In order to gain insights into the molecular pathways underlying progression of emphysema and explore computational strategies for identifying COPD therapeutics, we profiled gene expression in lung tissue samples obtained from regions within the same lung with varying amounts of emphysematous destruction from smokers with COPD (8 regions × 8 lungs = 64 samples). Regional emphysema severity was quantified in each tissue sample using the mean linear intercept (Lm) between alveolar walls from micro-CT scans.
We identified 127 genes whose expression levels were significantly associated with regional emphysema severity while controlling for gene expression differences between individuals. Genes increasing in expression with increasing emphysematous destruction included those involved in inflammation, such as the B-cell receptor signaling pathway, while genes decreasing in expression were enriched in tissue repair processes, including the transforming growth factor beta (TGFβ) pathway, actin organization, and integrin signaling. We found concordant differential expression of these emphysema severity-associated genes in four cross-sectional studies of COPD. Using the Connectivity Map, we identified GHK as a compound that can reverse the gene-expression signature associated with emphysematous destruction and induce expression patterns consistent with TGFβ pathway activation. Treatment of human fibroblasts with GHK recapitulated TGFβ-induced gene-expression patterns, led to the organization of the actin cytoskeleton, and elevated the expression of integrin β1. Furthermore, addition of GHK or TGFβ restored collagen I contraction and remodeling by fibroblasts derived from COPD lungs compared to fibroblasts from former smokers without COPD.
These results demonstrate that gene-expression changes associated with regional emphysema severity within an individual's lung can provide insights into emphysema pathogenesis and identify novel therapeutic opportunities for this deadly disease. They also suggest the need for additional studies to examine the mechanisms by which TGFβ and GHK each reverse the gene-expression signature of emphysematous destruction and the effects of this reversal on disease progression.
Chronic obstructive pulmonary disease (COPD) is a significant public health problem worldwide and the third leading cause of death in the United States . It is characterized by irreversible airflow limitation due to obstruction in the small conducting airways and emphysematous destruction of the gas exchanging tissue of the lung. Tobacco smoke is a significant risk factor for COPD and at least 25% of smokers develop this disease . Current theories concerning disease pathogenesis include an imbalance between protease and anti-protease activity, induced apoptosis of alveolar cells through deregulation of pathways involved in oxidative stress, chronic inflammation, and aberrant tissue remodeling that lead to the destruction of the extracellular matrix (ECM) in the lung [3, 4]. Lung repair and regeneration are potential processes to target with novel therapeutics in COPD as abnormal tissue repair by the epithelial-mesenchymal trophic unit can result in either fibrosis or destruction of the ECM . However, the molecular mechanisms responsible for the pathogenesis of COPD remain poorly understood.
Several groups have profiled gene expression in patients with and without COPD or between patients with varying levels of airflow obstruction in order to understand differences in gene expression related to COPD [6–11]. While these studies have provided an initial look into the COPD transcriptome, their results primarily relied on the use of lung function tests to define the presence or degree of COPD. Lung function phenotypes can neither distinguish between obstruction in the small airways and emphysematous destruction of the lung parenchyma nor provide information about regional differences in disease severity. Recently, McDonough et al. used micro-CT scans to quantify the degree of emphysema in different regions of lungs from patients with severe COPD by measuring the mean linear intercept (Lm), a morphological measurement of alveolar destruction . In order to gain insights into biological pathways associated with increasing emphysema severity within a patient and explore computational strategies for identifying COPD therapeutics, we obtained paired samples from eight regions at regular intervals between the apex and base of each explanted lung from six patients with severe COPD (Global Initiative for Chronic Obstructive Lung Disease (GOLD) stage IV) and two donor lungs. The degree of emphysematous destruction was quantified in one tissue sample from each region by Lm, while gene expression was profiled in the adjacent tissue sample from the same region.
We identified a number of genes whose expression is associated with increasing emphysematous destruction and found that pathways enriched among these genes were involved in the immune response and tissue remodeling. Using the Connectivity Map (CMap) , we found that the tripeptide Gly-His-Lys (GHK) was able to reverse the aberrant patterns of gene expression associated with increasing emphysema severity and induce patterns of gene expression consistent with transforming growth factor beta (TGFβ) pathway activation. Furthermore, we showed that by treating distal lung fibroblasts from COPD patients with GHK, we can restore normal contractile function through re-organization of the actin cytoskeleton and up-regulation of integrin-β1. These data further support the potential of GHK as a therapeutic in the treatment of emphysema.
Materials and methods
Sample acquisition and processing
Single lungs (n = 6) were removed from patients treated for severe COPD by double lung transplantation at the University of Pennsylvania. Donor lungs (n = 2) for which no suitable recipient was identified were released for research use from the Gift of Life Organ Procurement Organization in Philadelphia. This study was approved by the institutional review boards and conforms to the Helsinki Declaration. Written informed consent for use of these specimens and the relevant clinical and radiological data required for this research were obtained from each patient prior to surgery and from the next of kin of the persons whose donated lung was released for research. Each lung was removed from the thorax, cooled to 1.6°C, and transported to the laboratory where the bronchial stump was cannulated . The lung was then inflated using a compressed air source attached to an underwater seal to slowly increase transpulmonary pressure (PL) from 0 to 30 cmH2O. The specimen was then held at a transpulmonary pressure of 10 cmH2 Owhile frozen by liquid nitrogen vapor (-130°C). The frozen specimen had a multidetector CT scan followed by being cut into 2-cm thick slices in the same plane as the CT scan. Tissue samples were collected using a sharpened steel cylinder (cork bore diameter of 14 mm). One sample from a cluster of four core samples of lung obtained from each site was processed for micro-CT . A companion core from the same cluster was used for the gene profiling and validation studies reported here. The representative nature of these samples with respect to the entire lung was established by comparing the densities of the sampled sites with the frequency distribution of the densities in the entire lung on multidetector CT as reported in McDonough et al..
Measurement of mean linear intercept
The severity of emphysema within each core was estimated by measuring Lm. A micro-CT scan of each core provided approximately 1,000 contiguous 16-μm thick images. Lm was measured at 20 regularly spaced intervals of each of the micro-CT scans using a previously validated grid of test lines projected onto the image and a custom macro linked to specialized software (ImagePro Plus; MediaCybernetics (Rockville, MD, USA). The number of intercepts between these lines and tissue was counted. Lm was calculated as the total length of the test lines divided by the number of cross-overs with tissue (equal to the number of intercepts divided by 2).
Microarray sample processing
High molecular weight (mRNA-containing fraction) RNA was isolated from tissue cores using the miRNeasy Mini Kit (Qiagen, Valencia, CA, USA). RNA integrity was assessed using an Agilent 2100 Bioanalyzer and RNA purity was assessed using a NanoDrop spectrophotometer. RNA (1 μg) was processed and hybridized onto the Human Exon 1.0 ST array (Affymetrix Inc., Santa Clara, CA, USA) according to the manufacturer's protocol as previously described . Expression Console Version 1.1 (Affymetrix Inc.) was used to generate transcript-level gene expression estimates for the 'core' exon probesets via the robust multichip average (RMA) algorithm. Gene symbols of transcript IDs were retrieved using DAVID . These gene expression data are available through the Gene Expression Omnibus (GEO) under the accession GSE27597.
Microarray data analysis
Gene ij is the log2 expression value for sample i in patient j for a single gene. Slice is a fixed effect controlling for the position within the lung from which the sample core was obtained. The random term ε ij represents the random error, which was assumed to be normally distributed, α j represents the random effect for patient, and β0 represents the intercept. Model 2 contains an additional fixed effect term for emphysema severity measured by the natural log of Lm. A gene's expression profile was considered associated with Lm if model 2 fit better than model 1 as determined by a significant P-value from a likelihood ratio test between the two models after applying a false discovery rate (FDR) correction. In the immunohistochemistry experiments, these linear models were also used to examine the relationship between Lm and the volume fraction of tissue with positive staining by substituting volume fraction (Vv) for gene expression as a dependent variable. All statistical analyses were conducted using R statistical software v2.9.2 and the nlme package in Bioconductor v2.4 .
Functional enrichment analysis
Functional enrichment analysis was performed using DAVID 2008 or Gene Set Enrichment Analysis (GSEA) v2.0.7 [16, 18]. For DAVID, functional enrichment was examined among Gene Ontology categories, and KEGG and BIOCARTA pathways. All genes in the species Homo sapiens were used as a reference set. For GSEA, genes were ranked by the t-statistic of the βLm coefficient in the linear mixed-effects model and then analyzed for the enrichment of canonical pathways and Gene Ontology term gene sets obtained from MSigDB v2.5.
Connecting to other gene-expression datasets
Using GSEA, sets of genes reported to change with COPD-related phenotypes or with TGFβ treatment in other gene-expression studies were examined in a ranked list of genes ordered from most induced in severe emphysema to most repressed in severe emphysema by the t-statistic of the βLm coefficient in the linear mixed-effects model. Conversely, sets of genes we identified as significantly positively or negatively associated with Lm were examined within gene lists ranked by the degree of differential expression as determined by re-analyzing previously published COPD-or TGFβ-related microarray studies. See Additional file 1 for a description of the data normalization procedures and statistical analyses used to generate gene sets and/or ranked gene lists for each of the previously published gene-expression datasets.
In order to find compounds that reverse gene-expression patterns associated with emphysema severity, we generated separate signatures for each COPD or TGFβ gene-expression dataset examined in this study. Signatures were generated by identifying the 50 genes most up-regulated and the 50 genes most down-regulated with respect to a COPD or TGFβ-related phenotype. Each signature was queried against the CMap using the algorithm described by Lamb et al.. See Additional file 1 for a description of the statistical analysis used to generate each query signature for each phenotype within each dataset. The list of all CMap query signatures used in this analysis include: 1) genes that change in expression as a function of regional emphysema severity in this study; genes that change in expression with 2) forced expiratory volume in 1 second (FEV1), 3) FEV1/forced vital capacity (FVC), or 4) between cases versus controls in Bhattacharya et al.; genes that change in expression between 5) controls versus emphysema patients or between 6) controls versus α1-antitrypsin disease in Golpon et al.; genes that change in expression with 7) FEV1 or 8) diffusing capacity of carbon monoxide (DLCO) in Spira et al.; genes that change in expression with 9) FEV1, 10) FEV1/FVC, 11) DLCO, 12) non-smokers versus GOLD2, or 13) non-smokers versus GOLD3 in Wang et al.; and genes that change in expression with TGFβ treatment from 14) Qin et al., 15) Classen et al., 16) Renzoni et al., 17) Koinuma et al., and 18) Malizia et al.. For comparison of the CMap data to our in vitro studies of the effects of GHK in primary lung fibroblasts, raw data for GHK-treated and control samples were downloaded from the CMap website and normalized using MAS5.0 with the Affymetrix CDF. Genes were ranked by a paired t-test between treatment and controls of different batches and compared to gene sets of GHK and TGFβ treatment using GSEA.
Isolation and culture of lung fibroblasts
Lung tissue of former smokers (defined as quitting smoking for at least one year before surgery) with normal lung function or GOLD stage IV COPD was obtained from patients undergoing surgery for resection for pulmonary carcinoma or lung transplantation. Fibroblast cultures were established from parenchymal lung tissue by an explant technique as previously described . Isolated cells were characterized as fibroblasts by morphological appearance and expression pattern of specific proteins as described previously [24, 25]. Fibroblast cultures were stored into liquid nitrogen until use.
Fibroblast cultures at passage 3 were cultured in eight-well chamber slides (Gibco, Burlington, ON, Canada) in growth medium (DMEM, 10% fetal bovine serum (FBS), penicillin, and streptavidin from Invitrogen, Burlington, ON, Canada). After reaching 70% confluence, fibroblasts were cultured for 24 h in 1% FBS DMEM and then incubated with either TGFβ1 10ng/ml (Peprotech, Dollard des Ormeaux, Quebec, Canada), GHK 10 nM (Sigma, Markham, Ontario, Canada) or control media (1% FBS DMEM, penicillin, streptavidin) for a further 48 h. After stimulations, chamber slides were fixed with 4% paraformaldehyde for 20 minutes, blocked in 10% goat serum in phosphate-buffered saline (PBS) with 0.1% saponin for 1 h and then stained with integrin-β1 antibody (M-106, Santa Cruz Biotechnology, Santa Cruz, CA, USA) in 0.1% saponin in PBS for 2 h at room temperature. Following washing in PBS with 0.1% saponin and 0.1% Tween 20, secondary antibody conjugated with goat anti-Mouse IgG Alexa Fluor 488 and Phalloidin conjugated with Alexa Fluor 594 were incubated for 2 h at room temperature. Following final washes, cultures were incubated with DAPI 1 ng/ml and then coverslipped with cytoseal. Confocal images were acquired with a Leica AOBS SP2 laser scanning confocal microscope (Leica, Heidelberg, Germany). The images were overlaid and the contrast enhancements were performed on the images using Volocity software™ (Improvisions Inc., Boston, MA, USA) as previously described .
Collagen gel contraction assays
Fibroblast cultures at passage 3 were cultured in six-well tissue culture plates (Gibco, Canada) in growth medium (DMEM, 10% FBS, penicillin, and streptavidin from Invitrogen, Canada). After reaching 70% confluence, fibroblasts were cultured for 24 h in 1% FBS DMEM and then incubated with TGFβ1 10 ng/ml (Peprotec, Canada), GHK 10 nM (Sigma, Canada) or growth media control for a further 48 h. Prior to the end of the treatment time point, a 12-well tissue culture plate was incubated with 1% bovine serum albumin in DMEM for 2 h. The medium was removed and then 500 µl of 0.4 mg/ml type I collagen (BD Biosciences, Mississauga, ON, Canada) was added and allowed to polymerize for 8 h at 37°C. The treated fibroblasts were then trypsinized and seeded at 2 × 105 cells/500 µl of 1% FBS DMEM, penicillin, streptavidin in duplicate on the collagen gels and cultured for an additional 24 h at 37°C in 5% CO2. The gels were imaged before and after and the extent of gel contraction measured using Image Pro Software.
Multi-photon and second harmonic generation microscopy
The collagen gels were fixed with 4% paraformaldehyde for 20 minutes and washed in PBS with 0.1% saponin and 0.1% Tween 20 before being incubated with phalloidin conjugated with Alexa Fluor 594 for 1 h at room temp. Gels were then mounted on to a glass slide using Secure-seal™ imaging spacers (size 20 mm; Sigma) and aqueous mounting media. The gels were then imaged using second harmonic generation microscopy to determine fibrilar collagen as previously described . For each cell volume, Z-section images were compiled and the three-dimensional image restoration was performed using Volocity software (Improvisions, Inc.). A noise-removal filter with a kernel size of 3 × 3 was applied to these three-dimensional images.
Subject demographics for lung tissue samples
Lm mean ± SD (μm)
Lm range (μm)
716 ± 164
414 ± 82
724 ± 252
1,822 ± 1270
1,352 ± 599
1,097 ± 441
384 ± 47
289 ± 41
Pathways associated with regional emphysema severity
In order to demonstrate that the 127 gene signature is related to regional emphysema severity within individuals and not to differences between donors and COPD patients or to differences in levels of emphysema between COPD patients, we repeated the same statistical analysis while only including the five COPD patients with emphysema and standardizing the Lm measurements within each patient core to a mean of zero and a standard deviation of one (Z-score). Using GSEA, the sets of up-and down-regulated genes in the 127-gene signature identified in the previous analysis with all eight patients and unscaled Lm measurements were concordantly enriched among genes differentially expressed when only the five emphysema patients were analyzed with Z-scored Lm measurements, indicating that this gene signature is associated with regional emphysema severity (FDR <0.001, GSEA; see Additional file 6 for the enrichment plot).
Validation of pathways up-regulated in regions of severe emphysema
Validation of pathways down-regulated in regions of severe emphysema
To further validate these findings, we cultured human lung fibroblasts with and without TGFβ1 and found that the set of genes most induced by TGFβ1 were enriched among genes that decrease in expression with increasing regional emphysema severity (FDR <0.05, GSEA; see Additional file 8 for the GSEA enrichment plot and Additional file 1 for the fibroblast culture methods). Immunostaining of lung tissue from the same regions on which we performed gene expression analysis localized SMAD2, a down-stream signal transducer of TGFβ, to the alveolar and airway walls while members of the bone morphogenetic protein (BMP) pathway, including SMAD6 and SMAD1, were primarily seen in vascular endothelial cells (Additional file 9).
Relationship to expression profiles in other COPD studies
In order to show that the gene expression signature of regional emphysema severity is present in larger cohorts of patients with earlier stages of disease, we used GSEA to examine the relationship between genes associated with regional emphysema severity in this dataset and genes associated with COPD phenotypes in other cross-sectional studies [6–11]. The genes that decreased in expression with increasing emphysema severity were significantly enriched among genes down-regulated as a function of COPD-related phenotypes in four of the five previously published datasets that we examined (FDR <0.05, GSEA; Additional file 10). In addition, genes that increased in expression with increasing emphysema severity were enriched amongst genes up-regulated as a function of COPD-related phenotypes in three of the five datasets (FDR <0.05, GSEA). As an example, the enrichment of genes associated with regional emphysema severity among genes differentially expressed with the presence of emphysema in the dataset from Golpon et al. is shown in Figure 4c,d. Conversely, sets of genes reported to be differentially expressed with COPD in four of the six other cross-sectional studies were enriched among the genes changing in expression with increasing regional emphysema severity (FDR <0.05, GSEA; Additional file 10). Examples of genes validated by quantitative RT-PCR in this study and concordantly differentially expressed in other studies are shown in Additional file 11. Many of these datasets, such as Bhattacharya et al. and Wang et al., contained larger numbers of patients with a variety of stages of disease (for example, GOLD stage 0 through GOLD stage IV). The enrichment of genes associated with regional emphysema severity with COPD-related phenotypes in these other datasets suggests that the biological processes associated with increasing emphysema severity within a patient with severe COPD also vary in individuals with earlier stages of disease.
Prediction of novel therapeutics for emphysema
In order to identify compounds that might reverse the gene-expression pattern associated with progression of emphysema, we utilized the CMap , a compendium of microarray experiments that measure the effect of therapeutic compounds on gene expression in cancer cell lines. Signatures of genes that 1) change in expression with regional emphysema severity in this dataset, 2) change in expression with lung function measures in other datasets [6, 7, 9, 10], or 3) change in expression with TGFβ treatment in other datasets [19–23] were each used as separate queries into the CMap data. We found that gene expression changes resulting from treatment with the tripeptide GHK, a compound thought to accelerate wound healing [32, 33], were negatively correlated with expression patterns associated with increasing regional emphysema severity (P = 0.006) and the COPD-related expression patterns observed in Bhattacharya et al. and Golpon et al. (P < 0.05). In addition, the gene expression effects of GHK are similar to the effects of TGFβ treatment observed by Malizia et al. (P = 0.004).
Reversal of COPD-related phenotypes in fibroblasts by GHK
The goal of this study was to identify gene expression changes associated with regional emphysema severity in order to elucidate biological processes underlying the progression of emphysema and to identify potential COPD therapeutics. By measuring gene expression from regions of varying emphysema severity within the same lung and by using a morphologic measurement of airspace size (Lm), which reflects the degree of alveolar destruction, we were able to identify gene expression changes associated specifically with the emphysematous component of COPD.
Interestingly, there was significant enrichment between genes differentially expressed in COPD or associated with worsening lung function in other datasets and those we found to be associated with regional emphysema severity. Importantly, this similarity supports the notion that regional differences in emphysema severity reflect the processes that occur with general COPD pathogenesis and progression and are not only present in patients with end-stage disease. Overall, these observations suggest a similarity in the gene expression alterations that accompany airflow obstruction, gas exchange abnormalities, and alveolar destruction measured by Lm.
A common characteristic in the pathology of COPD is progressive lymphocyte infiltration of the small airways and alveolar walls . In addition, the formation of tertiary lymphoid organs within this infiltration suggests the presence of an adaptive immune response to persistent foreign or autoimmune antigens [37, 38]. The present study extends these observations by showing that the expression patterns of several components of the B-cell receptor signaling pathway have increased expression in regions of severe emphysema. Igα (CD79A) and Igβ(CD79B) are proteins that associate with the B-cell receptor and transmit its signal upon stimulation. Immunohistochemistry showed a significant relationship between the volume fraction of the airway wall and alveolar tissue positively stained for CD79A and an increase in Lm. This relationship supports an increased number of B cells in both airway wall and alveolar tissues and is consistent with the induction of CD79A during tissue destruction associated with the increase in Lm.
The TGFβ signaling pathway is involved in a variety of cellular processes, including immune response, extracellular matrix remodeling, angiogenesis, and cell differentiation. This pathway has also been implicated in a variety of diseases such as cancer and fibrosis . It has been hypothesized that the TGFβ pathway could play a role in COPD pathogenesis, but its role is not completely understood . Togo et al. found that fibroblasts isolated from COPD patients exhibited reduced chemotaxis, reduced nuclear to cytoplasmic ratios of phosphorylated SMAD3, and decreased α-smooth muscle actin production compared to controls when treated with TGFβ. Decreased mRNA expression or protein levels for TGFβ1, TGFBR1 , SMAD3 , SMAD6 , and SMAD7 [41, 43] have been reported in more advanced stages of COPD or fibroblasts from COPD patients. In both alveolar and bronchiolar epithelium of emphysematous lungs, a decrease in phosphorylated SMAD2 has been shown by immunohistochemistry . In normal human lung parenchyma, repair processes in response to mechanical injury are associated with increased TGFβ signaling, while a decrease in expression has been observed for TGFβ-related genes with worsening lung function in patients with COPD [25, 45]. Furthermore, association studies have identified both promoter and coding region polymorphisms in the TGFβ1 gene that associate with increased risk for COPD [46–48]. In the present study, we identified several components of the TGFβ and BMP pathways that have decreasing expression with increasing emphysema severity. In the BMP pathway, ACVRL1 and ENG are receptors involved in the phosphorylation of SMAD1 and are expressed in the mature lung vasculature. The changing expression of SMAD6 and SMAD1, their localization predominantly to vascular endothelial cells, and the roles of ACVRL1 and ENG in angiogenesis support the hypothesis of aberrant tissue remodeling in the lung vasculature during emphysema pathogenesis. In the TGFβ pathway, TGFBR2 is a receptor involved in the phosphorylation of SMAD2/3 and is important for many tissue remodeling processes, including wound repair. Moreover, genes found to be induced by TGFβ in diverse studies were down-regulated in regions of severe emphysema. The localization of SMAD2 to alveolar and airway tissue and the decreased TGFβ pathway activity seen with increasing emphysema severity support the hypothesis that a decrease in TGFβ pathway activity also contributes to emphysema pathogenesis.
As COPD remains a major public health concern due to lack of effective therapeutic strategies, we sought to use computational methods to identify compounds that might modulate molecular processes associated with emphysema pathogenesis. The CMap is a large compendium of microarray experiments that measures the effect of over 1,000 compounds on gene expression in several cell lines . By querying a gene expression signature of disease pathogenesis against the CMap dataset, one can find compounds that elicit a pattern of gene expression that is the opposite to the disease-related gene expression profile. This can lead to the hypothesis that such compounds, since they reverse the disease-related gene expression pattern, are potential therapeutics for that disease. This approach has been recently successful in the therapeutic repositioning of the antiulcer drug cimetidine to lung adenocarcinoma and the anticonvulsant drug topiramate to inflammatory bowel disease [49, 50]. In these studies, signatures for each disease were derived using several publicly available gene-expression datasets and queried in the CMap. Candidate compounds or drugs that could significantly reverse the disease-related signatures of gene expression were further validated in vitro, showing that this computational method is a viable approach for identifying novel therapeutics.
Using the CMap dataset, we identified a relationship between the gene expression changes induced by the tripeptide GHK and those that are repressed with increasing emphysema severity. Intriguingly, we further found that GHK-treatment induced a pattern of gene expression similar to that resulting from TGFβ pathway activation. We replicated both of these findings in human lung fibroblasts, which are the major interstitial cells that maintain tissue structural integrity by sculpting the connective tissue. GHK-Cu is a natural tripeptide that, in human plasma, can be found at a concentration of 200 ng/ml at the age of 20 years but drops to around 80 ng/ml by the age of 60 years . Characterization GHK-Cu in skin wound repair models suggests that it induces wound contraction, cell proliferation, angiogenesis, and increased expression of antioxidant enzymes and integrins [34, 51]. Direct evidence for the ability of GHK-Cu to promote wound healing comes from experimental rat models where GHK treatment causes an acceleration of healing and a concentration-dependent increase of connective tissue and other ECM components [32, 33]. These effects are consistent with the gene expression alterations induced by GHK and TGFβ treatment. Moreover, we confirmed these similarities by demonstrating that GHK and TGFβ induced significantly higher expression and re-organization of actin and integrin-β1 in distal lung fibroblasts.
We further assessed the ability of GHK and TGFβ to induce tissue contraction. As in previous studies , we demonstrated that distal lung fibroblasts derived from COPD patients have intrinsic defects in collagen I contraction compared to fibroblasts derived from former smokers without COPD. When fibroblasts from COPD lungs were treated with GHK or TGFβcontraction and remodeling of collagen gels was induced to levels comparable to fibroblasts from former smokers without COPD. We further demonstrated that the collagen contraction induced in COPD fibroblasts by GHK involves the organization of collagen I gels into collagen fibrils using multi-photon microscopy. Taken together, these data further support the hypothesis in which a wound-healing-like process is diminished as a function of emphysema progression and further suggest that this process is related to the TGFβ pathway.
While the number of subjects in this study for genomic analysis was small, the analysis of eight specimens per lung representing different degrees of emphysema from each individual allowed us to detect gene expression changes specifically associated with regional emphysema severity. We further demonstrated that these genes are concordantly differentially expressed in previous cross-sectional studies involving larger numbers of individuals with varying degrees of airflow limitation. These results validate the gene expression differences associated with regional emphysema severity in independent cohorts from different clinical settings and support the hypothesis that the genes whose expression is associated with regional emphysema severity reflect the activity of true disease-associated processes. As demonstrated by our micro-CT data, COPD is a heterogeneous disease within the lung . Further studies will be required to assess whether COPD-associated differences in ECM remodeling by distal fibroblasts in vitro is associated with the regional disease severity in the tissue from which the fibroblasts are derived.
This study has provided insights into molecular processes associated with emphysematous destruction of the lung and revealed mechanisms that contribute to the pathogenesis of COPD. Whole genome gene-expression analysis supports the role of the immune response in regional emphysema and elucidates additional pathways involved in the process of emphysematous destruction. The suggestion that progressive emphysematous destruction is associated with down-regulation of genes involved in or downstream of tissue remodeling and wound repair pathways supports a role for defects in ECM homeostasis and angiogenesis in the emphysematous destruction that occurs with chronic inflammation in COPD. We propose that these processes could be linked through decreased TGFβ pathway activation. These data are supported by our identification of GHK as a compound with the potential to mimic TGFβ pathway activity and induce collagen contraction, an important functional component of wound repair.
bone morphogenetic protein
chronic obstructive pulmonary disease
diffusing capacity of carbon monoxide
Dulbecco's modified Eagle's medium
fetal bovine serum
false discovery rate
forced expiratory volume in 1 second
forced vital capacity
Global Initiative for Chronic Obstructive Lung Disease
Gene Set Enrichment Analysis
mean linear intercept
transforming growth factor
We thank A Wright, D Horng, and P Sanchez for supporting studies in this manuscript; W Elliott for help with the immunohistochemistry; F Shaheen for help with cell culture; K Steiling for reviewing this manuscript; T Abraham at the UBC James Hogg Research Centre Imaging Cellular Imaging and Biophysics Core facility for help with imaging the collagen gels. This work was funded by the National Heart, Lung and Blood Institute (R01HL095388 to AS and MEL), National Center for Advancing Translational Science (UL1 TR000157 to AS and MEL), National Science Foundation (Integrative Graduate Education and Research Traineeship to JDC and JEZ), British Columbia Lung Association, Canadian Institute for Health Research, Parker B Francis Foundation (Senior Fellowship to TLH) and the Dutch Asthma Foundation and European Respiratory Society (International Research Fellowship to CAB).
- Miniño AM, Xu J, Kochanek KD, Statistics V: Deaths: Preliminary Data for 2008. National Vital Statistics Reports Volume 59. 2010, NVSS, [http://www.cdc.gov/nchs/data/nvsr/nvsr59/nvsr59_02.pdf]
- Lange P, Scharling H, Fabricius P, Vestbo J: Developing COPD: a 25 year follow up study of the general population. Thorax. 2006, 61: 935-939. 10.1136/thx.2006.062802.PubMedPubMed CentralView Article
- Park JW, Ryter SW, Choi AMK: Functional significance of apoptosis in chronic obstructive pulmonary disease. COPD. 2007, 4: 347-353. 10.1080/15412550701603775.PubMedView Article
- Postma DS, Timens W: Remodeling in asthma and chronic obstructive pulmonary disease. Proc Am Thorac Soc. 2006, 3: 434-439. 10.1513/pats.200601-006AW.PubMedView Article
- Rennard SI, Bailey KL: Chronic obstructive pulmonary disease exacerbations: accurate and easy measurement promises much. Am J Respir Crit Care Med. 2012, 185: 1139-1141. 10.1164/rccm.201202-0227ED.PubMedView Article
- Golpon HA, Coldren CD, Zamora MR, Cosgrove GP, Moore MD, Tuder RM, Geraci MW, Voelkel NF: Emphysema lung tissue gene expression profiling. Am J Respir Cell Mol Biol. 2004, 31: 595-600. 10.1165/rcmb.2004-0008OC.PubMedView Article
- Bhattacharya S, Srisuma S, Demeo DL, Shapiro SD, Bueno R, Silverman EK, Reilly JJ, Mariani TJ: Molecular biomarkers for quantitative and discrete COPD phenotypes. Am J Respir Cell Mol Biol. 2009, 40: 359-367. 10.1165/rcmb.2008-0114OC.PubMedPubMed CentralView Article
- Ning W, Li C-J, Kaminski N, Feghali-Bostwick C, Alber SM, Di YP, Otterbein SL, Song R, Hayashi S, Zhou Z, Pinsky DJ, Watkins SC, Pilewski JM, Sciurba FC, Peters DG, Hogg JC, Choi AMK: Comprehensive gene expression profiles reveal pathways related to the pathogenesis of chronic obstructive pulmonary disease. Proc Natl Acad Sci USA. 2004, 101: 14895-14900. 10.1073/pnas.0401168101.PubMedPubMed CentralView Article
- Spira A, Beane J, Pinto-Plata V, Kadar A, Liu G, Shah V, Celli B, Brody JS: Gene expression profiling of human lung tissue from smokers with severe emphysema. Am J Respir Cell Mol Biol. 2004, 31: 601-610. 10.1165/rcmb.2004-0273OC.PubMedView Article
- Wang I-M, Stepaniants S, Boie Y, Mortimer JR, Kennedy B, Elliott M, Hayashi S, Loy L, Coulter S, Cervino S, Harris J, Thornton M, Raubertas R, Roberts C, Hogg JC, Crackower M, O'Neill G, Paré PD: Gene expression profiling in patients with chronic obstructive pulmonary disease and lung cancer. Am J Respir Crit Care Med. 2008, 177: 402-411. 10.1164/rccm.200703-390OC.PubMedView Article
- Francis SMS, Larsen JE, Pavey SJ, Bowman RV, Hayward NK, Fong KM, Yang I: Expression profiling identifies genes involved in emphysema severity. Respir Res. 2009, 10: 81-10.1186/1465-9921-10-81.PubMedPubMed CentralView Article
- McDonough JE, Yuan R, Suzuki M, Seyednejad N, Elliott WM, Sanchez PG, Wright AC, Gefter WB, Litzky L, Coxson HO, Paré PD, Sin DD, Pierce RA, Woods JC, McWilliams AM, Mayo JR, Lam SC, Cooper JD, Hogg JC: Small-airway obstruction and emphysema in chronic obstructive pulmonary disease. N Engl J Med. 2011, 365: 1567-1575. 10.1056/NEJMoa1106955.PubMedPubMed CentralView Article
- Lamb J, Crawford ED, Peck D, Modell JW, Blat IC, Wrobel MJ, Lerner J, Brunet J-P, Subramanian A, Ross KN, Reich M, Hieronymus H, Wei G, Armstrong S, Haggarty SJ, Clemons P, Wei R, Carr S, Lander ES, Golub TR: The Connectivity Map: using gene-expression signatures to connect small molecules, genes, and disease. Science. 2006, 313: 1929-1935. 10.1126/science.1132939.PubMedView Article
- Choong CK, Haddad FJ, Martinez C, Hu DZ, Pierce J, Meyers BF, Patterson GA, Cooper JD: A simple, reproducible, and inexpensive technique in the preparation of explanted emphysematous lungs for ex vivo studies. J Thorac Cardiovasc Surg. 2005, 130: 922-923.PubMed
- Zhang X, Liu G, Lenburg ME, Spira A: Comparison of smoking-induced gene expression on Affymetrix Exon and 3'-based expression arrays. Genome Inform. 2007, 18: 247-257.PubMed
- Dennis G, Sherman BT, Hosack D, Yang J, Gao W, Lane HC, Lempicki R: DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003, 4: P3-10.1186/gb-2003-4-5-p3.PubMedView Article
- Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, Hornik K, Hothorn T, Huber W, Iacus S, Irizarry R, Leisch F, Li C, Maechler M, Rossini AJ, Sawitzki G, Smith C, Smyth G, Tierney L, Yang JYH, Zhang J: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5: R80-10.1186/gb-2004-5-10-r80.PubMedPubMed CentralView Article
- Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette M, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP: Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005, 102: 15545-15550. 10.1073/pnas.0506580102.PubMedPubMed CentralView Article
- Qin H, Chan MWY, Liyanarachchi S, Balch C, Potter D, Souriraj IJ, Cheng ASL, Agosto-Perez FJ, Nikonova EV, Yan PS, Lin H-J, Nephew KP, Saltz JH, Showe LC, Huang THM, Davuluri RV: An integrative ChIP-chip and gene expression profiling to model SMAD regulatory modules. BMC Syst Biol. 2009, 3: 73-10.1186/1752-0509-3-73.PubMedPubMed CentralView Article
- Classen S, Zander T, Eggle D, Chemnitz JM, Brors B, Büchmann I, Popov A, Beyer M, Eils R, Debey S, Schultze JL: Human resting CD4+ T cells are constitutively inhibited by TGF beta under steady-state conditions. J Immunol. 2007, 178: 6931-6940.PubMedView Article
- Renzoni E, Abraham DJ, Howat S, Shi-Wen X, Sestini P, Bou-Gharios G, Wells AU, Veeraraghavan S, Nicholson AG, Denton CP, Leask A, Pearson JD, Black CM, Welsh KI, du Bois RM: Gene expression profiling reveals novel TGFbeta targets in adult lung fibroblasts. Respir Res. 2004, 5: 24-10.1186/1465-9921-5-24.PubMedPubMed CentralView Article
- Koinuma D, Tsutsumi S, Kamimura N, Taniguchi H, Miyazawa K, Sunamura M, Imamura T, Miyazono K, Aburatani H: Chromatin immunoprecipitation on microarray analysis of Smad2/3 binding sites reveals roles of ETS1 and TFAP2A in transforming growth factor beta signaling. Mol Cell Biol. 2009, 29: 172-186. 10.1128/MCB.01038-08.PubMedPubMed CentralView Article
- Malizia AP, Keating DT, Smith SM, Walls D, Doran PP, Egan JJ: Alveolar epithelial cell injury with Epstein-Barr virus upregulates TGFbeta1 expression. Am J Physiol Lung Cell Mol Physiol. 2008, 295: L451-460. 10.1152/ajplung.00376.2007.PubMedView Article
- Noordhoek JA, Postma DS, Chong LL, Menkema L, Kauffman HF, Timens W, van Straaten JFM, van der Geld YM: Different modulation of decorin production by lung fibroblasts from patients with mild and severe emphysema. COPD. 2005, 2: 17-25. 10.1081/COPD-200050678.PubMedView Article
- Pechkovsky DV, Hackett TL, An SS, Shaheen F, Murray L, Knight D: Human lung parenchyma but not proximal bronchi produces fibroblasts with enhanced TGF-beta signaling and alpha-SMA expression. Am J Respir Cell Mol Biol. 2010, 43: 641-651. 10.1165/rcmb.2009-0318OC.PubMedView Article
- Hackett T-L, Warner SM, Stefanowicz D, Shaheen F, Pechkovsky DV, Murray L, Argentieri R, Kicic A, Stick SM, Bai TR, Knight D: 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-133. 10.1164/rccm.200811-1730OC.PubMedView Article
- Abraham T, Carthy J, McManus B: Collagen matrix remodeling in 3-dimensional cellular space resolved using second harmonic generation and multiphoton excitation fluorescence. J Struct Biol. 2010, 169: 36-44. 10.1016/j.jsb.2009.07.023.PubMedView Article
- Bignon J, Khoury F, Even P, Andre J, Brouet G: Morphometric study in chronic obstructive bronchopulmonary disease. Pathologic, clinical, and physiologic correlations. Am Rev Respir Dis. 1969, 99: 669-695.PubMed
- Faith JJ, Hayete B, Thaden JT, Mogno I, Wierzbowski J, Cottarel G, Kasif S, Collins JJ, Gardner TS: Large-scale mapping and validation of Escherichia coli transcriptional regulation from a compendium of expression profiles. PLoS Biol. 2007, 5: e8-10.1371/journal.pbio.0050008.PubMedPubMed CentralView Article
- Chambers RC, Leoni P, Kaminski N, Laurent GJ, Heller RA: Global expression profiling of fibroblast responses to transforming growth factor-beta1 reveals the induction of inhibitor of differentiation-1 and provides evidence of smooth muscle cell phenotypic switching. Am J Pathol. 2003, 162: 533-546. 10.1016/S0002-9440(10)63847-3.PubMedPubMed CentralView Article
- Verrecchia F, Chu ML, Mauviel A: Identification of novel TGF-beta/Smad gene targets in dermal fibroblasts using a combined cDNA microarray/promoter transactivation approach. J Biol Chem. 2001, 276: 17058-17062. 10.1074/jbc.M100754200.PubMedView Article
- Siméon A, Wegrowski Y, Bontemps Y, Maquart FX: Expression of glycosaminoglycans and small proteoglycans in wounds: modulation by the tripeptide-copper complex glycyl-L-histidyl-L-lysine-Cu(2+). J Invest Dermatol. 2000, 115: 962-968. 10.1046/j.1523-1747.2000.00166.x.PubMedView Article
- Maquart FX, Bellon G, Chaqour B, Wegrowski J, Patt LM, Trachy RE, Monboisse JC, Chastang F, Birembaut P, Gillery P: In vivo stimulation of connective tissue accumulation by the tripeptide-copper complex glycyl-L-histidyl-L-lysine-Cu2+ in rat experimental wounds. J Clin Invest. 1993, 92: 2368-2376. 10.1172/JCI116842.PubMedPubMed CentralView Article
- Pickart L: The human tri-peptide GHK and tissue remodeling. J Biomater Sci Polym Ed. 2008, 19: 969-988. 10.1163/156856208784909435.PubMedView Article
- Togo S, Holz O, Liu X, Sugiura H, Kamio K, Wang X, Kawasaki S, Ahn Y, Fredriksson K, Skold CM, Mueller KC, Branscheid D, Welker L, Watz H, Magnussen H, Rennard SI: Lung fibroblast repair functions in patients with chronic obstructive pulmonary disease are altered by multiple mechanisms. Am J Respir Crit Care Med. 2008, 178: 248-260. 10.1164/rccm.200706-929OC.PubMedPubMed CentralView Article
- Hogg JC, Chu F, Utokaparch S, Woods R, Elliott WM, Buzatu L, Cherniack RM, Rogers RM, Sciurba FC, Coxson HO, Paré PD: The nature of small-airway obstruction in chronic obstructive pulmonary disease. N Engl J Med. 2004, 350: 2645-2653. 10.1056/NEJMoa032158.PubMedView Article
- Brusselle GG, Demoor T, Bracke KR, Brandsma C-a, Timens W: Lymphoid follicles in (very) severe COPD: beneficial or harmful?. Eur Respir J. 2009, 34: 219-230. 10.1183/09031936.00150208.PubMedView Article
- van der Strate BW, Postma DS, Brandsma C-A, Melgert BN, Luinge M, Geerlings M, Hylkema MN, van den Berg A, Timens W, Kerstjens HaM: Cigarette smoke-induced emphysema: A role for the B cell?. Am J Respir Crit Care Med. 2006, 173: 751-758. 10.1164/rccm.200504-594OC.PubMedView Article
- Blobe GC, Schiemann WP, Lodish HF: Role of transforming growth factor beta in human disease. N Engl J Med. 2000, 342: 1350-1358. 10.1056/NEJM200005043421807.PubMedView Article
- Morty RE, Königshoff M, Eickelberg O: Transforming growth factor-beta signaling across ages: from distorted lung development to chronic obstructive pulmonary disease. Proc Am Thorac Soc. 2009, 6: 607-613. 10.1513/pats.200908-087RM.PubMedView Article
- Zandvoort A, Postma DS, Jonker MR, Noordhoek J, Vos JTWM, van der Geld YM, Timens W: Altered expression of the Smad signalling pathway: implications for COPD pathogenesis. Eur Respir J. 2006, 28: 533-541. 10.1183/09031936.06.00078405.PubMedView Article
- Zandvoort A, Postma DS, Jonker MR, Noordhoek J, Vos JTWM, Timens W: Smad gene expression in pulmonary fibroblasts: indications for defective ECM repair in COPD. Respir Res. 2008, 9: 83-10.1186/1465-9921-9-83.PubMedPubMed CentralView Article
- Springer J, Scholz FR, Peiser C, Groneberg D, Fischer A: SMAD-signaling in chronic obstructive pulmonary disease: transcriptional down-regulation of inhibitory SMAD 6 and 7 by cigarette smoke. Biol Chem. 2004, 385: 649-653.PubMedView Article
- Leppäranta O, Myllärniemi M, Salmenkivi K, Kinnula VL, Keski-Oja J, Koli K: Reduced phosphorylation of the TGF-Beta signal transducer Smad2 in emphysematous human lung. COPD. 2009, 6: 234-241. 10.1080/15412550903049173.PubMedView Article
- Gosselink JV, Hayashi S, Elliott WM, Xing L, Chan B, Yang L, Wright C, Sin D, Paré PD, Pierce J, Pierce R, Patterson A, Cooper J, Hogg JC: Differential expression of tissue repair genes in the pathogenesis of chronic obstructive pulmonary disease. Am J Respir Crit Care Med. 2010, 181: 1329-1335. 10.1164/rccm.200812-1902OC.PubMedPubMed CentralView Article
- Smolonska J, Wijmenga C, Postma DS, Boezen HM: Meta-analyses on suspected chronic obstructive pulmonary disease genes: a summary of 20 years' research. Am J Respir Crit Care Med. 2009, 180: 618-631. 10.1164/rccm.200905-0722OC.PubMedView Article
- van Diemen CC, Postma DS, Vonk JM, Bruinenberg M, Nolte IM, Boezen HM: Decorin and TGF-beta1 polymorphisms and development of COPD in a general population. Respir Res. 2006, 7: 89-10.1186/1465-9921-7-89.PubMedPubMed CentralView Article
- van Diemen CC, Postma DS, Aulchenko YS, Snijders PJLM, Oostra B, van Duijn CM, Boezen HM: Novel strategy to identify genetic risk factors for COPD severity: a genetic isolate. Eur Respir J. 2010, 35: 768-775. 10.1183/09031936.00054408.PubMedView Article
- Dudley JT, Sirota M, Shenoy M, Pai RK, Roedder S, Chiang AP, Morgan A, Sarwal MM, Pasricha PJ, Butte AJ: Computational repositioning of the anticonvulsant topiramate for inflammatory bowel disease. Sci Transl Med. 2011, 3: 96ra76-10.1126/scitranslmed.3002648.PubMedPubMed Central
- Sirota M, Dudley JT, Kim J, Chiang AP, Morgan A, Sweet-Cordero A, Sage J, Butte AJ: Discovery and preclinical validation of drug indications using compendia of public gene expression data. Sci Transl Med. 2011, 3: 96ra77-10.1126/scitranslmed.3001318.PubMedPubMed Central
- Kang Y-A, Choi H-R, Na J-I, Huh C-H, Kim M-J, Youn S-W, Kim K-H, Park K-C: Copper-GHK increases integrin expression and p63 positivity by keratinocytes. Arch Dermatol Res. 2009, 301: 301-306. 10.1007/s00403-009-0942-x.PubMedView Article
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.