- Open Access
miR-638 regulates gene expression networks associated with emphysematous lung destruction
Genome Medicine volume 5, Article number: 114 (2013)
Chronic obstructive pulmonary disease (COPD) is a heterogeneous disease characterized by varying degrees of emphysematous lung destruction and small airway disease, each with distinct effects on clinical outcomes. There is little known about how microRNAs contribute specifically to the emphysema phenotype. We examined how genome-wide microRNA expression is altered with regional emphysema severity and how these microRNAs regulate disease-associated gene expression networks.
We profiled microRNAs in different regions of the lung with varying degrees of emphysema from 6 smokers with COPD and 2 controls (8 regions × 8 lungs = 64 samples). Regional emphysema severity was quantified by mean linear intercept. Whole genome microRNA and gene expression data were integrated in the same samples to build co-expression networks. Candidate microRNAs were perturbed in human lung fibroblasts in order to validate these networks.
The expression levels of 63 microRNAs (P < 0.05) were altered with regional emphysema. A subset, including miR-638, miR-30c, and miR-181d, had expression levels that were associated with those of their predicted mRNA targets. Genes correlated with these microRNAs were enriched in pathways associated with emphysema pathophysiology (for example, oxidative stress and accelerated aging). Inhibition of miR-638 expression in lung fibroblasts led to modulation of these same emphysema-related pathways. Gene targets of miR-638 in these pathways were amongst those negatively correlated with miR-638 expression in emphysema.
Our findings demonstrate that microRNAs are altered with regional emphysema severity and modulate disease-associated gene expression networks. Furthermore, miR-638 may regulate gene expression pathways related to the oxidative stress response and aging in emphysematous lung tissue and lung fibroblasts.
Chronic obstructive pulmonary disease (COPD) is a growing epidemic and currently the third leading cause of death in the US . Mechanisms leading to the development of COPD are poorly understood, and there are currently no effective therapies that modify long-term lung function decline in patients with the disease.
The heterogeneity of COPD (that is, patients affected by varying degrees of emphysema and small airway disease) may account for some of the difficulty in identifying effective therapies . Although COPD is diagnosed as airflow obstruction that is not fully reversible, measurements of airflow obstruction are weakly correlated with emphysema severity and provide no information on regional differences in emphysema within the lung [3, 4]. Emphysema is independently associated with poor clinical outcomes and there is evidence of an emphysema-predominant COPD phenotype in which airflow limitation is less important . Understanding the molecular dysregulation leading to emphysema may lead to a better understanding of this phenotype and more targeted therapeutics.
Prior studies of global gene and microRNA expression in COPD were based on airflow limitation [6–13]. Thus, they were unable to capture the effects of emphysema independent of small airway disease. Furthermore, as these studies were all based on a case–control design, they were unable to account for differences in gene expression associated with varying disease severity within an individual’s lung. We have developed an approach to quantify emphysema severity within different regions of the lung from patients with severe COPD, using micro-computed tomography (micro-CT) data to evaluate mean linear intercept (Lm), a measure of alveolar destruction . Using these Lm measurements, we recently identified a signature of gene expression associated with increasing regional emphysema severity within individual lungs . We found that many of these genes are involved in biological pathways such as inflammation, extracellular matrix (ECM) remodeling, and tissue repair. Furthermore, we were able to predict novel therapeutics using our emphysema gene signature and the Connectivity Map, showing that the compound tripeptide Gly-His-Lys (GHK) reversed the emphysema-associated gene expression changes . However, the molecular mechanisms responsible for the disease-associated changes in gene expression within the lung are not well understood.
MicroRNAs are 19 to 25 nucleotide non-coding RNAs that regulate gene expression through inhibition of mRNA translation or induction of mRNA degradation . It is thought that they fine tune complex regulatory networks by inhibiting multiple genes at once [18, 19]. They have repeatedly been implicated in lung disease and are altered with airflow obstruction; however, their role in emphysema is unclear [12, 13, 20]. MicroRNAs are beginning to be studied in clinical trials as potential biomarkers and therapeutics; thus, understanding how they contribute to disease pathology could lead to significant clinical impact .
In this study, we assessed whether microRNA expression changes as regional emphysema severity increases within an individual’s lung and whether these microRNA changes regulate gene expression alterations associated with disease. We performed microRNA profiling in the same lung tissue samples in which we analyzed mRNA expression and evaluated the relationship between microRNA and mRNA expression changes and emphysema severity. A candidate microRNA, miR-638, was further studied in human lung fibroblasts derived from patients with severe COPD to better understand its role in regulating gene expression patterns associated with emphysema pathogenesis in vivo.
A schematic of the study design is shown in Additional file 1: Figure S1.
Sample acquisition and preparation
A detailed description of sample preparation and acquisition has been described previously [14, 22]. Briefly, lungs were obtained from six subjects undergoing lung transplantation for severe COPD and two donors without COPD. Lungs were cut into 2 cm slices along the axial plane and a cluster of 1 cm cores was taken from each slice. One core from each slice was evaluated by micro-CT while RNA and microRNA were extracted from an adjacent core. Emphysema severity was calculated for each cluster by mean linear intercept (Lm), an estimate of alveolar size, using the micro-CT images. Approximately 1,000 contiguous images were obtained from each scanned core and Lm was measured at 20 regularly spaced intervals in each scan. Lm was calculated by placing a grid of parallel lines over each image and determining the sum of the length of all grid lines divided by the number of intercepts between the alveolar septae and grid lines. Thus, a high Lm corresponds to more severe emphysema as there will be fewer alveolar septae for the grid lines to intercept, and thus a lower denominator.
This study was approved by the institutional review boards at University of British Columbia, Boston University Medical Campus, University of Groningen (for fibroblast work), and University of Pennsylvania (where lungs were procured). The study 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.
RNA isolation and microarray processing
RNA was extracted from tissue cores and fractionated into high molecular weight (mRNA-containing) and low molecular weight (LMW; microRNA-containing) portions using the miRNeasy mini kit (QIAGEN, Valencia, CA, USA). RNA fraction integrity was assessed using the Agilent 2100 Bioanalyzer and purity was measured using the NanoDrop spectrophotometer.
LMW RNA (<200 ng ng) was processed and hybridized onto NCode version 3 microRNA microarrays (Invitrogen, Carlsbad, CA, USA) as previously described . These arrays contained probes for 1,053 microRNAs from 6 species, including 467 human microRNAs, printed in triplicate. Raw data were extracted using GenePix Pro 4.0 (Axon/Molecular Devices, Sunnyvale, CA, USA).
The mRNA data were generated using Human Exon 1.0 ST microarrays as described in our previously published gene-expression study .
Raw microRNA data underwent preprocessing and filtering as previously described . Briefly, the data were quantile normalized, and log transformed (limma package, Bioconductor, R) [24, 25]. A microRNA probe was called 'present' if its signal intensity was two standard deviations above the average of the background probes in the sample. MicroRNA probes were filtered out if they were not 'present' in at least 80% of the samples. If two or three replicate probes were called present, the median of the present replicates was taken to obtain a single expression value within each sample. If 0 or 1 probes were called present, the microRNA was filtered out. Samples were then filtered based on two quality metrics: (1) principal component analysis, whereby samples were excluded if they were significant outliers based on the first two principle components; and (2) present/absent filtering, whereby samples were excluded if less than 50% of all human microRNA probesets were 'absent' in the sample. Three samples were excluded by at least one of these two quality metrics.
Transcript-level gene expression values were generated for our previously published gene expression paper via the robust multichip average (RMA) algorithm (affy package, Bioconductor, R) [15, 26]. Gene symbols were retrieved using the Entrez Gene Custom chip definition file . Raw data for microRNA and gene expression are publicly available at the Gene Expression Omnibus (GEO) under the accessions GSE49881 and GSE27597, respectively.
MicroRNA (and gene) expression values were related to emphysema severity using linear mixed-effects models that controlled for differences between patients and position in the lung from which the sample was taken (nlme package, Bioconductor, R) . Two linear mixed effects models were used:
where miRNA i is the response variable representing the log2 expression of microRNA i. βPatient*Patient represents the random effect of 'patient', which controls for samples originating from the same lung. βSlice*Slice is a fixed effect that controls for the region in the lung the sample came from (slice number from apex to base). β0 is the intercept and ϵ i is the error for miRNA i . In Equation 2 there is a second fixed effect, βLm*Lm, which controls for emphysema severity as defined by the natural log of Lm. miRNA i was associated with emphysema severity if the model in Equation 2 (with the emphysema term) fit better than the model in Equation 1 using a likelihood ratio test. Adjustment for multiple comparisons was done by applying the Benjamini-Hochberg false discovery rate (FDR) .
The correlation analysis between microRNAs and mRNAs was done using the same equations with genes ('Genei') substituted as the response variable for microRNAs ('miRNA i ') and 'βmiRNA*miRNA' substituted as the fixed effect 'βLm*Lm' in Equation 2. Thus, this analysis evaluated the correlation between each gene and each microRNA after correction for patient and region of lung. Correlated microRNA-gene pairs (FDR <0.25) were filtered to include only microRNA-predicted target gene pairs as determined by at least one of five target prediction algorithms (Targetscan, Pictar, Tarbase, miRBase, and microRNA.org) [30–34]. The lenient FDR cutoff of 0.25 for significance of microRNA-target gene pairs (corresponding to P = 0.002) was used given the hypothesis-generating nature of this correlation analysis. Furthermore, genes were then filtered to include only predicted targets, and targets were validated through in vitro experiments.
Real time PCR validation
Differential expression of eight microRNAs was evaluated with quantitative RT-PCR (qRT-PCR). Three to four samples from two of the eight patients were used to evaluate each microRNA. qRT-PCR was performed using the Taq-Man Small RNA Assay (Applied Biosystems) with 10 ng LMW RNA per manufacturer’s protocol for miR-181d, miR-30c, miR-150, miR18a-3p, miR-211, miR-296-5p, miR-483-3p and miR-638. PCR was carried out using the StepOnePlus Real-Time PCR system with 40 cycles of amplification and data acquisition. Samples were run in triplicate and normalized to U6. Analysis was done using the comparative CT method and expression levels were compared to microarray by correlation analysis.
Gene set enrichment analysis (GSEA) was used to determine whether pathways that may be biologically important in COPD were associated with microRNAs in our data set . Canonical pathway gene sets were obtained from the Molecular Signature Database (MSigDB) . For each differentially expressed microRNA, a ranked gene list was generated using the t-statistics from the linear mixed effect models that correlated mRNA and microRNA expression.
Isolation, culture, and microRNA transfection in primary lung fibroblasts
Primary lung fibroblasts were cultured from peripheral lung tissue of a subject with stage IV COPD undergoing lung transplant surgery. Cells were isolated and characterized using our explant technique as described previously . Fibroblasts were cultured in complete culture medium (Ham’s F12, 10% fetal bovine serum, penicillin, streptomycin and glutamin (all from Lonza, Verviers, Belgium)) and used for experiments after five passages.
Transfection experiments (n = 3) with miR-638 inhibitors and controls were performed using the HiPerFect transfection reagent (QIAGEN). In short, 2.5 × 105 fibroblasts per well were seeded in six-wells plates. MicroRNA inhibitors in a final concentration of 50 nM (QIAGEN) for mir-638 were added dropwise to the cells. The Inhibitor negative control from QIAGEN was used as the negative control. The transfection reagent was removed after 24 hours and cells were washed with Ham’s culture medium without fetal bovine serum and harvested using trypsin and used for RNA isolation. RNA was isolated using the miRNeasy kit.
RNA was hybridized to Affymetrix Human Gene 1.0 ST arrays. Differential expression between inhibitor and control experiments was determined by empirical Bayes-moderated t-tests. Differences in overall expression of miR-638 predicted targets between inhibitor and control experiments was determined by two-tailed Kolmogorov-Smirnov testing (limma, R) .
Pathway analysis in fibroblasts and comparison of fibroblast and emphysema datasets
GSEA analysis was done to determine biological pathways enriched with miR-638 inhibition. Genes were ranked by fold change difference between control and inhibitor and canonical pathway gene sets were run against this ranked list. Our fibroblast dataset was compared to our emphysema dataset by using a gene set created from the 100 most upregulated miR-638 predicted targets (by fold change) between miR-638 inhibitor experiments and control experiments. This gene set was run against a list of genes ranked by their correlation with miR-638.
Eight sample pairs from different lung slices were obtained from six subjects requiring lung transplant for COPD and two organ donors (64 samples total). Detailed characteristics of subjects and lung specimens were described previously and are shown in Additional file 2. Briefly, the six subjects with COPD had a forced expiratory volume in 1 second (FEV1) <25% predicted (severe disease). Subjects with COPD had higher and more variable mean linear intercept (Lm) values for their tissue samples than donors except for one subject who had small airways disease without emphysema. All subjects were former smokers except one donor. Four samples were excluded from further analysis: one as an outlier by Lm, the others as outliers by microRNA microarray quality metrics.
MicroRNA expression is altered with regional emphysema severity
After microarray preprocessing, 397 of the initial 467 human microRNAs were detected above background and used for further analysis. Sixty-three microRNAs were correlated with Lm and thus with regional emphysema severity (P-value < 0.05, 20 expected by chance at this threshold; Figure 1; Additional file 3). In contrast, only 18 microRNAs were correlated with position in the lung the sample came from (for example, apical or basal, P < 0.05), and only three of these 18 microRNAs correlated with Lm (Additional file 4). The microRNAs upregulated with increasing emphysema severity included broadly conserved and poorly conserved microRNAs, with the most upregulated including miR-520e and miR-302d from the miR-93 family, miR-92a, miR-638, miR-211, and miR-150. The downregulated microRNAs included clusters from several broadly conserved families. These include Let-7c, Let-7d, Let-7e, and Let-7f from the Let-7/miR-98 family, miR-181c and miR-181d from the miR-181/4262 family, and miR-30a-3p, miR-30c, miR-30e-5p, and miR-30e-3p from the miR-30/384-5p family. The association with regional emphysema severity was validated using qRT-PCR. Expression values derived from microarray and RT-PCR data were significantly correlated for four of the eight miRNAs (Pearson; P ≤ 0.05; Additional file 1: Figure S2) and just below statistical significance but trending in the correct direction for a fifth microRNA, miR-181d (P = 0.06). The expression of these five microRNAs and two of the other three all trended in the expected direction (same as in the array data) with increasing Lm.
Gene expression is correlated with microRNA expression in the emphysematous lung
To identify which microRNAs may regulate gene expression associated with emphysema pathogenesis, we generated a microRNA-mRNA correlation network by integrating the microRNA expression data with global mRNA expression data generated from the same lung specimens. We found that 1,079 microRNA-mRNA pairs were significantly correlated (371 negative correlations, 708 positive correlations) using linear mixed effect models (FDR <0.25) and were predicted to interact by at least one of five sequence-based target prediction algorithms.
Of the 63 differentially expressed microRNAs, 51 had at least one positively or negatively correlated predicted target at FDR <0.25 (Figure 2; Additional file 5). The most highly connected microRNAs in the network included miR-638, miR-18a-3p, miR-483-3p, miR-181d, and miR-30c, which had greater than 50 positively or negatively correlated predicted targets, suggesting that these microRNAs may be important regulators of gene expression associated with emphysema severity.
MicroRNAs regulate pathways involved in COPD pathogenesis
To gain insight into functional significance of the microRNA-mRNA interactions, we next used GSEA to determine the association between each microRNA and 880 canonical pathway gene sets. Despite the microRNA-mRNA target analysis showing both positive and negative correlations, the significantly enriched canonical pathways were almost all inversely correlated with upregulated microRNAs, suggesting an inhibitory role for microRNAs with increasing disease severity. Unsupervised hierarchical clustering of the pathways significantly enriched amongst upregulated microRNAs identified several distinct clusters of downregulated functional categories associated with expression of multiple upregulated microRNAs, including miR-638 (Figure 3). Pathway gene sets most downregulated in emphysema, including those involved in cell surface signaling and ECM maintenance (for example, transforming growth factor β, integrin, and wnt-signaling), were inversely correlated with the expression of multiple upregulated microRNAs. We also identified additional pathways important to emphysema pathology that were inversely correlated with upregulated microRNAs, including pathways in tissue and cell repair (for example, DNA repair and synthesis, cell cycle maintenance and traverse, RNA processing, and proteasome maintenance/ubiquitination), and mitochondrial biogenesis (for example, oxidative phosphorylation, electron transport chain).
Inhibition of miR-638 in lung fibroblasts leads to overexpression of its predicted targets
We next examined the specific effects of miR-638 on the expression of its predicted targets in primary human lung fibroblasts obtained from a subject with severe COPD, as fibroblasts play a key role in the aberrant tissue repair and ECM maintenance observed in emphysema pathogenesis [37, 38]. miR-638 was chosen as it had many correlated targets, was anti-correlated with multiple pathways involved in ECM dysregulation and tissue repair by GSEA, and has been shown previously to be expressed in normal human lung fibroblasts . We knocked down miR-638 in primary human lung fibroblasts and measured global gene expression using Affymetrix Human Gene 1.0 ST arrays. We observed a significant enrichment of miR-638 predicted targets amongst upregulated genes in miR-638 knockdown experiments compared to controls (P = 2.6 × 10-10, Kolmogorov-Smirnov test; Figure 4A). miR-638 inhibition caused modest changes in the expression of many target genes (1.15 to 1.4 fold change for significant genes in our experiments), similar to previously published studies on microRNA transfection [40, 41].
miR-638 targets are involved in oxidative stress-related pathways and are anti-correlated with miR-638 expression in emphysematous lung tissue
To evaluate the functional significance of miR-638 expression in emphysema, we identified the canonical pathways affected by miR-638 knock-down in fibroblasts, and studied the overlap with genes and pathways correlated with miR-638 in emphysematous lung tissue. GSEA analysis revealed enrichment of many of the same biological pathways among the genes increased in fibroblasts after miR-638 inhibition as those whose expression was anti-correlated with miR-638 expression in emphysema (FDR <0.25; Figure 4C; Additional file 6: Table S5). miR-638 inhibition in fibroblasts led to upregulation of pathway gene sets involved in DNA repair, cell-cycle maintenance and traverse, RNA transcription, mitochondrial biogenesis, telomere maintenance, and processing of damaged proteins. Sixty-nine (56%) of these pathways were also anti-correlated with miR-638 expression in emphysematous lung tissue, and previously implicated in oxidative stress and accelerated lung aging responses.
The 100 most upregulated predicted targets by fold change with miR-638 inhibition in fibroblasts were enriched amongst predicted targets anti-correlated with miR-638 expression in emphysema (P < 0.001, Kolmogorov-Smirnov test; Figure 4B). These predicted targets included genes involved in cell proliferation (ADAM15, ARHGDIA, COMMD1, DHCR7, HDAC5, MAD2L2, PFKL, YPEL3), autophagy and protein degradation (ATG9A, GANAB), mitochondrial functioning (DHCR7, ERAL1, SLC25A1, STARD3, TOMM40), DNA damage response (APBB1, KRT7, MAD2L2), oxidative stress response (CARHSP1), and ECM remodeling (LTBP4). These results demonstrate that by modulating miR-638 in lung fibroblasts obtained from a patient with COPD, we can reproduce emphysema-related patterns of gene expression observed in lung tissue and further suggest a role for miR-638 in emphysema pathogenesis.
In this study, we have shown that regional emphysema severity is associated with alterations in microRNA expression and that a subset of these microRNA changes are correlated with changes in their predicted gene target expression. While 51 of the 63 altered microRNAs correlated with at least one of their predicted targets, a much smaller subset had many correlated targets. Five of these microRNAs (miR-638, miR-181d, miR-18a-3p, miR-30c, and miR-483-3p) were correlated with ≥50 of their predicted targets, suggesting that these microRNAs may play a key role in gene regulation in emphysema.
Unique strengths of our study include our ability to relate microRNA changes to the emphysematous component of COPD, to account for regional differences in emphysema within an individual patient’s lung, and to relate microRNA changes to mRNA changes within these same samples. Other studies of microRNA expression in COPD have used airflow obstruction alone as the COPD phenotype to associate with molecular alterations [12, 13]. However, emphysema correlates poorly with severity of airflow obstruction and contributes independently to outcomes in COPD [2, 4, 42]. In addition, emphysematous destruction and small airway disease occur to different degrees throughout the lungs of patients with COPD, and airflow obstruction measurements cannot capture these differences. We isolated multiple samples from regions within an individual patient’s lung varying by degree of emphysema and by location (for example, apical, basal), and used Lm to quantify the degree of emphysema severity within each sample. By using this study design, we were able to specifically identify those microRNAs whose expression is related to emphysema severity/progression within the lung, controlling both for microRNA expression differences between individuals and for differences in lung location the sample came from (for example, apical, basal). Furthermore, the unique study design allowed us to identify microRNA-mRNA networks that underlie the pathogenesis of emphysematous destruction within an individual’s lung.
An analysis of microRNA-target interactions revealed numerous significant negative and positive interactions suggesting that microRNAs are an important part of the regulatory network in emphysema. Somewhat surprisingly, 66% of significant correlations were in the positive direction. Although this seems counterintuitive, given that microRNAs inhibit gene expression, positive correlations are noted throughout the literature, including in a study of human airway microRNA and gene expression [19, 43–45]. Positive correlations are attributed to microRNAs acting together in a network with other regulatory processes to influence gene expression (Additional file 1: Figure S3). Both microRNAs and the mRNAs they directly inhibit may be under the control of common upstream regulators. This can lead to positive correlation between the expression of the microRNA and mRNA pairs even though the microRNAs are functioning to decrease mRNA expression levels in this context (a coexpressed or 'incoherent' circuit). Negative microRNA-target gene correlations suggest upregulated microRNAs act with other processes to lower target gene expression (an expected or 'coherent' circuit).
Our pathway analysis suggests microRNA-gene interactions may contribute to emphysema pathogenesis by influencing pathways involved in ECM maintenance, the oxidative stress response and lung aging. The pathway gene sets most downregulated with increasing emphysema severity, including those involved in ECM maintenance, were anti-correlated with upregulated microRNA expression. In addition, pathways involved in tissue repair processes (for example, DNA repair, cell cycle traverse, and lysosome maintenance) and mitochondrial functioning, which are dysregulated in accelerating lung aging processes and critical to emphysema pathogenesis, were also significantly enriched amongst genes anti-correlated with upregulated microRNAs. Multiple microRNAs have previously been found to regulate tissue repair processes and to act cooperatively or redundantly to control these functions . Our findings support this and suggest that dysregulation of multiple microRNAs may contribute to emphysema pathogenesis.
We studied the function of miR-638, a microRNA whose expression increases with increasing emphysema severity, in regulating gene expression in human lung fibroblasts. Peripheral lung fibroblasts were studied as they have been implicated in COPD as key regulators of accelerated lung aging and dysregulated ECM repair. Our analysis suggests that the microRNAs changing with emphysema severity may be regulators of these processes. Fibroblasts are the main producers of the ECM and have been shown by us and others to be deficient in collagen contraction, a model property of the repair process, in patients with COPD [37, 38]. Accelerated aging processes (for example, cellular senescence) are upregulated in fibroblasts in COPD and have been implicated as mediators of fibroblast functional insufficiency, and thus their decreased contribution to ECM repair, in COPD . miR-638 was chosen as it has many correlated targets and was previously found to be associated with cellular senescence and DNA damage [39, 48], suggesting it contributes to accelerated lung aging as well. It is expressed in human lung fibroblasts and its expression increases as a function of cellular senescence . miR-638 overexpression in bronchial epithelial cells also increases DNA damage and decreases DNA repair, both of which contribute to cellular senescence and tissue repair dysregulation, after exposure to benzo(a)pyrene, a polycyclic aromatic hydrocarbon and component of cigarette smoke .
We found that miR-638 knock-down in COPD fibroblasts preferentially leads to over-expression of its predicted targets and that overexpressed targets are enriched amongst gene targets anti-correlated with miR-638 expression in emphysematous lung tissue. Thus, although miR-638 is positively correlated with two-thirds of its predicted targets in emphysema, those targets that are most affected with modulation of this microRNA in COPD fibroblasts are anti-correlated with miR-638 expression in emphysema. While the positive correlations imply that miR-638 is part of a regulatory network in which upstream regulators lead to microRNA-mRNA target coexpression (incoherent regulatory circuit), we show evidence that miR-638 is an important regulator in the network as well, as it is anti-correlated with those targets it directly inhibits most in COPD fibroblasts. These findings clearly suggest that miR-638 contributes to the gene expression changes related to emphysema severity.
miR-638 inhibition in fibroblasts led to upregulation of pathways that were also dysregulated in emphysematous lung tissue and involved in responses consistent with accelerated lung aging and the oxidative stress response. Furthermore, several miR-638 gene targets involved in these processes were dysregulated in vitro with miR-638 inhibition and anti-correlated with miR-638 expression in lung tissue with increasing emphysema severity. Accelerated aging, which can occur independently of and earlier than chronological aging, is a collection of molecular and cellular alterations that has been implicated in several diseases (for example, emphysema, Parkinson’s disease, diabetes, coronary artery disease). Some of these processes include: (1) cellular senescence, telomere attrition, persistent DNA damage response, and cell cycle arrest, leading to decreased cell proliferation and stimulation of T-helper type 1-associated inflammation; (2) mitochondrial dysfunction leading to increased production of free radicals and oxidative stress; and (3) decreased removal of damaged proteins leading to further cellular senescence and increased tissue damage [49–52]. Chronic exposure to cigarette smoke leads to chronic oxidative stress that has been implicated in turning on these aging processes prematurely. This is hypothesized to lead to the tissue destruction and inability to replace cells that we observe in emphysema. Given our pathway and microRNA-target interaction analyses, we propose that miR-638 contributes to accelerated lung aging and oxidative stress responses in emphysema (Figure 5).
While most studies of microRNAs in the lung show regulation of just one or a few genes, a strength of our study is that we are able to show that microRNAs can alter the expression of multiple genes leading to pathology. Amongst these genes is HDAC5, which encodes histone deacetylase 5, and has been shown to decrease with increasing COPD severity. Through epigenetic modification, this decrease in HDAC5 has been hypothesized to contribute to cell cycle arrest and increased inflammation . Another of these genes, ADAM15, which encodes one of the disintegrin and metalloproteinase domain-containing proteins, has been shown to be protective against emphysema as ADAM15−/− mice develop more severe emphysema in response to cigarette smoke . Many of the other gene targets have never been studied in COPD, and thus represent an area for future investigation as possible important regulators of disease.
While only a relatively small number of patients were profiled in this study, we obtained multiple samples per patient, allowing us to examine the microRNA alterations associated with regional differences in emphysema severity. Sequence-based computational algorithms for predicting microRNA-mRNA interactions generally have poor sensitivity and specificity. However, by including predictions from multiple algorithms and by intersecting this with microRNA-gene coexpression networks generated in subjects with severe emphysema we were able to identify a microRNA capable of modulating its targets in emphysema-related pathways in vitro.
Using our analytical approach we have predicted which microRNAs are most important to emphysema pathogenesis and how they are interacting with target genes to contribute to disease. Additional studies are needed to investigate how the microRNAs specifically contribute to the emphysema phenotype. In addition, given the vast public health impact of COPD, the scarcity of treatment options, and the recent advances in the development of microRNA-directed therapeutics, future work may focus on these microRNAs as potential treatment targets in emphysema. We have shown that the downstream inhibition of multiple genes by at least one of these microRNAs, miR-638, likely functions to affect emphysema pathogenesis. Thus, targeting microRNAs in emphysema may have a significant impact by altering the expression of many genes relevant to disease pathogenesis.
This study demonstrates that microRNA expression is altered with regional emphysema severity and that these microRNAs can regulate the gene expression changes associated with emphysema. Our study design allowed us to examine emphysematous destruction independent of airflow limitation and to control for heterogeneity of emphysema severity between and within patients, which limited previous analyses of microRNA expression in COPD. We identified several microRNAs that may be involved in emphysema pathogenesis and appear to modulate disease-specific gene expression networks. We demonstrated that miR-638, which was upregulated in emphysema, likely participates in COPD pathology by responding to oxidative stress, contributing to an accelerated lung aging response and leading to an impaired ability of the damaged lung to replace the damaged ECM. MicroRNAs identified as important to emphysema pathogenesis are intriguing candidates for further study in biomarker and therapeutic target development in emphysema.
chronic obstructive pulmonary disease
false discovery rate
gene set enrichment analysis
mean linear intercept
low molecular weight
Molecular Signature Database
quantitative real time polymerase chain reaction.
Mannino DM, Buist AS: Global burden of COPD: risk factors, prevalence, and future trends. Lancet. 2007, 370: 765-773. 10.1016/S0140-6736(07)61380-4.
Friedlander AL, Lynch D, Dyar LA, Bowler RP: Phenotypes of chronic obstructive pulmonary disease. COPD. 2007, 4: 355-384. 10.1080/15412550701629663.
Celli BR, MacNee W, ATS/ERS Task Force: Standards for the diagnosis and treatment of patients with COPD: a summary of the ATS/ERS position paper. Eur Respir J. 2004, 23: 932-946. 10.1183/09031936.04.00014304.
Gelb AF, Hogg JC, Muller NL, Schein MJ, Kuei J, Tashkin DP, Epstein JD, Kollin J, Green RH, Zamel N, Elliott WM, Hadjiaghai L: Contribution of emphysema and small airways in COPD. Chest. 1996, 109: 353-359. 10.1378/chest.109.2.353.
Makita H, Nasuhara Y, Nagai K, Ito Y, Hasegawa M, Betsuyaku T, Onodera Y, Hizawa N, Nishimura M, Hokkaido COPD, Cohort Study Group: Characterization of phenotypes based on severity of emphysema in chronic obstructive pulmonary disease. Thorax. 2007, 62: 932-937. 10.1136/thx.2006.072777.
Francis SM, Larsen JE, Pavey SJ, Bowman RV, Hayward NK, Fong KM, Yang IA: Expression profiling identifies genes involved in emphysema severity. Respir Res. 2009, 10: 81-10.1186/1465-9921-10-81.
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.
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.
Ning W, Lee J, Kaminski N, Feghali-Bostwick CA, Watkins SC, Pilewski JM, Peters DG, Hogg JC, Choi AM: Comprehensive analysis of gene expression on GOLD-2 versus GOLD-0 smokers reveals novel genes important in the pathogenesis of COPD. Proc Am Thorac Soc. 2006, 3: 466-10.1513/pats.200603-031MS.
Spira A, Beane J, Shah V, Liu G, Schembri F, Yang X, Palma J, Brody JS: Effects of cigarette smoke on the human airway epithelial cell transcriptome. Proc Natl Acad Sci U S A. 2004, 101: 10143-10148. 10.1073/pnas.0401422101.
Wang IM, 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.
Ezzie ME, Crawford M, Cho JH, Orellana R, Zhang S, Gelinas R, Batte K, Yu L, Nuovo G, Galas D, Diaz P, Wang K, Nana-Sinkam SP: Gene expression networks in COPD: MicroRNA and mRNA regulation. Thorax. 2012, 67: 122-131. 10.1136/thoraxjnl-2011-200089.
Pottelberge GR, Mestdagh P, Bracke KR, Thas O, Durme YM, Joos GF, Vandesompele J, Brusselle GG: MicroRNA expression in induced sputum of smokers and patients with chronic obstructive pulmonary disease. Am J Respir Crit Care Med. 2011, 183: 898-906. 10.1164/rccm.201002-0304OC.
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.
Campbell JD, McDonough JE, Zeskind JE, Hackett TL, Pechkovsky DV, Brandsma CA, Suzuki M, Gosselink JV, Liu G, Alekseyev YO, Xiao J, Zhang X, Hayashi S, Cooper JD, Timens W, Postma DS, Knight DA, Marc LE, James HC, Avrum S: A gene expression signature of emphysema-related lung destruction and its reversal by the tripeptide GHK. Genome Med. 2012, 4: 67-
Lamb J, Crawford ED, Peck D, Modell JW, Blat IC, Wrobel MJ, Lerner J, Brunet JP, Subramanian A, Ross KN, Reich M, Hieronymus H, Wei G, Armstrong SA, Haggarty SJ, Clemons PA, Wei R, Carr SA, 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.
Bartel DP: MicroRNAs: target recognition and regulatory functions. Cell. 2009, 136: 215-233. 10.1016/j.cell.2009.01.002.
Satoh J: Molecular network analysis of human microRNA targetome: from cancers to Alzheimer's disease. BioData Min. 2012, 5: 17-0381-5-17
Tsang J, Zhu J, van Oudenaarden A: MicroRNA-mediated feedback and feedforward loops are recurrent network motifs in mammals. Mol Cell. 2007, 26: 753-767. 10.1016/j.molcel.2007.05.018.
Yang IA, Relan V, Wright CM, Davidson MR, Sriram KB, Savarimuthu Francis SM, Clarke BE, Duhig EE, Bowman RV, Fong KM: Common pathogenic mechanisms and pathways in the development of COPD and lung cancer. Expert Opin Ther Targets. 2011, 15: 439-456.
Nana-Sinkam SP, Croce CM: Clinical applications of microRNAs in cancer. Clin Pharmacol Ther. 2013, 93: 98-104. 10.1038/clpt.2012.192.
Choong CK, Haddad FJ, Martinez C, Hu DZ, Pierce JA, 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.
Schembri F, Sridhar S, Perdomo C, Gustafson AM, Zhang X, Ergun A, Lu J, Liu G, Zhang X, Bowers J, Vaziri C, Ott K, Sensinger K, Collins JJ, Brody JS, Getts R, Lenburg ME, Spira A: MicroRNAs as modulators of smoking-induced gene expression changes in human airway epithelium. Proc Natl Acad Sci U S A. 2009, 106: 2319-2324. 10.1073/pnas.0806383106.
Smyth GK: Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004, 3: Article3-
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 JY, Zhang J: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5: R80-10.1186/gb-2004-5-10-r80.
Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, Speed TP: Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003, 4: 249-264. 10.1093/biostatistics/4.2.249.
Dai M, Wang P, Boyd AD, Kostov G, Athey B, Jones EG, Bunney WE, Myers RM, Speed TP, Akil H, Watson SJ, Meng F: Evolving gene/transcript definitions significantly alter the interpretation of GeneChip data. Nucleic Acids Res. 2005, 33: e175-10.1093/nar/gni179.
Pinheiro J, Bates D, DebRoy S, Sarkar D: Nlme: linear and nonlinear mixed effects models. R package version 3. 2008, Vienna, Austria: R Foundation for Statistical Computing, 3:1–104
Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Statist Soc B. 1995, 57: 289-300.
Lewis BP, Burge CB, Bartel DP: Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell. 2005, 120: 15-20. 10.1016/j.cell.2004.12.035.
Rajewsky N, Socci ND: Computational identification of microRNA targets. Dev Biol. 2004, 267: 529-535. 10.1016/j.ydbio.2003.12.003.
Betel D, Wilson M, Gabow A, Marks DS, Sander C: The microRNA.org resource: targets and expression. Nucleic Acids Res. 2008, 36: D149-D153.
Sethupathy P, Corda B, Hatzigeorgiou AG: TarBase: a comprehensive database of experimentally supported animal microRNA targets. RNA. 2006, 12: 192-197.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, 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 U S A. 2005, 102: 15545-15550. 10.1073/pnas.0506580102.
Ito K, Ito M, Elliott WM, Cosio B, Caramori G, Kon OM, Barczyk A, Hayashi S, Adcock IM, Hogg JC, Barnes PJ: Decreased histone deacetylase activity in chronic obstructive pulmonary disease. N Engl J Med. 2005, 352: 1967-1976. 10.1056/NEJMoa041892.
The Molecular Signature Database. [http://www.broadinstitute.org/gsea/msigdb/index.jsp]
Noordhoek JA, Postma DS, Chong LL, Chong LL, Menkema L, Kauffman HF, Timens W, van Straaten JF, 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.
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.
Maes OC, Sarojini H, Wang E: Stepwise up-regulation of microRNA expression levels from replicating to reversible and irreversible growth arrest states in WI-38 human fibroblasts. J Cell Physiol. 2009, 221: 109-119. 10.1002/jcp.21834.
Farh KK, Grimson A, Jan C, Lewis BP, Johnston WK, Lim LP, Burge CB, Bartel DP: The widespread impact of mammalian MicroRNAs on mRNA repression and evolution. Science. 2005, 310: 1817-1821. 10.1126/science.1121158.
Brennecke J, Stark A, Russell RB, Cohen SM: Principles of microRNA-target recognition. PLoS Biol. 2005, 3: e85-10.1371/journal.pbio.0030085.
Han MK, Agusti A, Calverley PM, Celli BR, Criner G, Curtis JL, Fabbri LM, Goldin JG, Jones PW, Macnee W, Make BJ, Rabe KF, Rennard SI, Sciurba FC, Silverman EK, Vestbo J, Washko GR, Wouters EF, Martinez FJ: Chronic obstructive pulmonary disease phenotypes: the future of COPD. Am J Respir Crit Care Med. 2010, 182: 598-604. 10.1164/rccm.200912-1843CC.
Solberg OD, Ostrin EJ, Love MI, Peng JC, Bhakta NR, Hou L, Nguyen C, Solon M, Nguyen C, Barczak AJ, Zlock LT, Blagev DP, Finkbeiner WE, Ansel KM, Arron JR, Erle DJ, Woodruff PG: Airway epithelial miRNA expression is altered in asthma. Am J Respir Crit Care Med. 2012, 186: 965-974. 10.1164/rccm.201201-0027OC.
El Baroudi M, Cora D, Bosia C, Osella M, Caselle M: A curated database of miRNA mediated feed-forward loops involving MYC as master regulator. PLoS One. 2011, 6: e14742-10.1371/journal.pone.0014742.
Osella M, Bosia C, Cora D, Caselle M: The role of incoherent microRNA-mediated feedforward loops in noise buffering. PLoS Comput Biol. 2011, 7: e1001101-10.1371/journal.pcbi.1001101.
Su WL, Kleinhanz RR, Schadt EE: Characterizing the role of miRNAs within gene regulatory networks using integrative genomics techniques. Mol Syst Biol. 2011, 7: 490-
Lee J, Sandford A, Man P, Sin DD: Is the aging process accelerated in chronic obstructive pulmonary disease?. Curr Opin Pulm Med. 2011, 17: 90-97. 10.1097/MCP.0b013e328341cead.
Li D, Wang Q, Liu C, Duan H, Zeng X, Zhang B, Li X, Zhao J, Tang S, Li Z, Xing X, Yang P, Chen L, Zeng J, Zhu X, Zhang S, Zhang Z, Ma L, He Z, Wang E, Xiao Y, Zheng Y, Chen W: Aberrant expression of miR-638 contributes to benzo(a)pyrene-induced human cell transformation. Toxicol Sci. 2012, 125: 382-391. 10.1093/toxsci/kfr299.
Karrasch S, Holz O, Jorres RA: Aging and induced senescence as factors in the pathogenesis of lung emphysema. Respir Med. 2008, 102: 1215-1230. 10.1016/j.rmed.2008.04.013.
Wu JJ, Quijano C, Chen E, Liu H, Cao L, Fergusson MM, Rovira II, Gutkind S, Daniels MP, Komatsu M, Finkel T: Mitochondrial dysfunction and oxidative stress mediate the physiological impairment induced by the disruption of autophagy. Aging (Albany NY). 2009, 1: 425-437.
Finkel T, Holbrook NJ: Oxidants, oxidative stress and the biology of ageing. Nature. 2000, 408: 239-247. 10.1038/35041687.
Provinciali M, Cardelli M, Marchegiani F: Inflammation, chronic obstructive pulmonary disease and aging. Curr Opin Pulm Med. 2011, 17: S3-10.
Nakajima T, McManus D, Tesfaigzi Y, McDonald JD, Owen CA: ADAM 15 protects the lung from cigarette smoke-induced lung inflammation and pulmonary emphysema in mice [abstract]. Am J Respir Crit Care Med Conference Abstract Issue. 2011, 183: A6426-
We thank JE McDonough, JV Gosselink and M Suzuki for help with tissue sampling and LM analyses; G. Lui, Y Alekseyev, X Zhang and J Xiao for array processing; JD Cooper for providing lung tissue; MR Jonker and W Kooistra for fibroblast experiments; and PG Woodruff for his guidance regarding microRNA-mRNA interactions.
This work was funded by the National Heart, Lung and Blood Institute (R01HL095388 to AS and MEL), National Center for Advancing Translational Science (UO1 TR000157 to AS and MEL), British Columbia Lung Association, Canadian Institute for Health Research (to JH and DAK), and the Dutch Asthma Foundation and European Respiratory Society (International Research Fellowships to CAB) and the ‘Stichting Astma Bestrijding’(to CAB), National Science Foundation (Integrative Graduate Education and Research Traineeship to JDC).
The authors declare that they have no competing interests.
Stephanie A Christenson and Corry-Anke Brandsma contributed as co-first authors. AS, SAC, CAB, JDC, DAK, DSP, ML, and WT conceived and designed the experiments. SC analyzed the array data with the assistance of AS, JC, and ML. CB and DVP performed the in vitro experiments. JCH, CAB, DSP, and WT contributed materials. SAC, CAB, JDC, MEL, JCH, DSP, WT and AS wrote and/or revised the manuscript. All authors read and approved the final version for publication.
Stephanie A Christenson, Corry-Anke Brandsma contributed equally to this work.
Electronic supplementary material
Additional file 2: Table S1: Subject demographics for the eight participants, six requiring transplantation for COPD and two organ donors. Those with COPD had a forced expiratory volume in 1 second / forced vital capacity (FEV1/FVC) <70% predicted and FEV1 <25% predicted. One COPD subject had small airway disease only (no emphysema). Some subjects carried other diagnoses: one COPD subject had von Willebrand disease, one COPD subject and one donor had hypertension, and one COPD subject had α1 antitrypsin disease. All subjects were former smokers with the exception of one donor who was a never-smoker. (XLS 32 KB)
Additional file 6: Table S5: Pathway gene sets anti-correlated with miR-638 in emphysema and upregulated with miR-638 knock-down in fibroblasts. Normalized enrichment score (NES), P-value, and FDR value for GSEA of all overlapping gene sets in both datasets at FDR <0.25. Gene sets are sorted by function, all associated with the oxidative stress and abnormal lung aging response: protein processing, mitochondrial dysfunction, DNA damage response, carbohydrate metabolism, diseases associated with accelerated aging, cell cycle traverse, and RNA transcription and processing. (XLS 17 KB)
About this article
Cite this article
Christenson, S.A., Brandsma, C., Campbell, J.D. et al. miR-638 regulates gene expression networks associated with emphysematous lung destruction. Genome Med 5, 114 (2013). https://doi.org/10.1186/gm519
- Chronic Obstructive Pulmonary Disease
- microRNA Expression
- Human Lung Fibroblast
- Emphysema Severity
- Small Airway Disease