Identification of common microRNA between COPD and non-small cell lung cancer through pathway enrichment analysis
BMC Genomic Data volume 22, Article number: 41 (2021)
Different factors have been introduced which influence the pathogenesis of chronic obstructive pulmonary disease (COPD) and non-small cell lung cancer (NSCLC). COPD as an independent factor is involved in the development of lung cancer. Moreover, there are certain resemblances between NSCLC and COPD, such as growth factors, activation of intracellular pathways, as well as epigenetic factors. One of the best approaches to understand the possible shared pathogenesis routes between COPD and NSCLC is to study the biological pathways that are activated. MicroRNAs (miRNAs) are critical biomolecules that implicate the regulation of several biological and cellular processes. As such, the main goal of this study was to use a systems biology approach to discover common dysregulated miRNAs between COPD and NSCLC, one that targets most genes within common enriched pathways.
To reconstruct the miRNA-pathways for each disease, we used the microarray miRNA expression data. Then, we employed “miRNA set enrichment analysis” (MiRSEA) to identify the most significant joint miRNAs between COPD and NSCLC based on the enrichment scores. Overall, our study revealed the involvement of the targets of miRNAs (such as has-miR-15b, hsa-miR-106a, has-miR-17, has-miR-103, and has-miR-107) in the most important common biological pathways.
According to the promising results of the pathway analysis, the identified miRNAs can be utilized as the new potential signatures for therapy through understanding the molecular mechanisms of both diseases.
Chronic obstructive pulmonary disease (COPD) is a lung-related disease specified by the continuous respiratory symptoms and boosted inflammatory response owing to harmful gases and particles [1, 2]. On the one hand, COPD raises oxidative stress leading to DNA damage, chronic exposure, repression of the DNA repair mechanisms, and cellular proliferation ; on the other hand, lung cancer as the fifth cancer leading to global mortality is usually classified into two main histologic types: non-small cell lung cancer (NSCLC) and small cell lung cancer (SCLC) . Moreover, the mutations in oncogenes can also lead to lung cancer, and as a result, cell proliferation and forming a tumor [4, 5]. Furthermore, cell proliferation and unsuppressed cell growth are the known characteristics of cancer progression in which several genes and proteins are involved, especially, the kinases and kinase receptors .
The rate of lung cancer in patients with COPD is nearly five times more than that of smokers without COPD ; besides, the overlap between COPD and lung cancer can be due to joint genetic susceptibility as well as the smoking-related processes . COPD is recognized as an autonomous risk factor for lung cancer, particularly for NSCLC as the most prevalent lung cancer type . Both COPD and NSCLC are mostly caused by cigarette smoking  through inducing inflammation and oxidative stress in the lung . Some common processes would contribute to the development of COPD and lung cancer in patients, such as abnormal immunity, cell proliferation, apoptosis, and chromatin modifications .
MicroRNAs are a category of functional non-coding RNAs containing 20 ~ 24 nucleotides, what negatively governs mRNA stability and/or suppresses mRNA translation through binding to the 3′ untranslated region [12, 13]. The role of miRNAs in a wide range of cellular process, including proliferation, cell cycle, differentiation, apoptosis, and metastasis has been reported . MiRNAs are involved in the initiation and development of disparate cancer types, while they are dysregulated in many cancers. Moreover, the alteration in the mRNA expression levels is also correlated with several diseases such as cardiovascular diseases, immunity- or inflammation-related diseases, and COPD [14, 15]. Furthermore, miRNAs function as oncogenes or tumor suppressors through regulating their target genes. It should be noted that miRNAs have great potential to be used as therapeutic targets, therefore, the determination and visualization of their positions in the regulatory pathways will be helpful in the development of novel medications . If a miRNA is related to the physiological process, it certainly regulates a gene or multiple genes in a corresponding pathway.
There are many common pathways that are activated in COPD and NSCLC . Athyros and colleagues, for one, found the impairment of several steps in the reverse cholesterol transport pathway via systematic inflammation in COPD . Moreover, Aldehydes identified the elevation of histone 3 phosphorylation in cigarette smokers via the activation of proliferative pathways, including the phosphatidylinositol-3 kinase (PI3K) / protein kinase B (PKB/Akt)  and MAPK pathway [20, 21]. The KEGG database is a series of biological pathways wherein many genes, proteins and other products are involved; however, the information about miRNAs is not mentioned in them. Cong Pian et al.  designed a new pathway database with the aid of KEGG plus miRNAs and integrated the human miRNA-target interactions with KEGG pathways using the hypergeometric test. Furthermore, C. Brinkrolf et al.  introduced a platform called VANESA for reconstructing, visualizing, and analyzing biological networks, to predict human miRNAs that may be co-expressed with genes involved in the KEGG pathway.
The aim of this study is to identify the most significant miRNAs as the new biomarkers which are common between COPD and NSCLC via analyzing the shared pathways between both diseases. To this aim, we considered two miRNA datasets related to COPD and NSCLC and normalized each dataset; then, we enriched both datasets to detect those pathways that contained more target genes for each miRNA list.
Thus, we detected those miRNAs that targeted more genes within the shared pathways and had more metabolic and genetic impact on the enriched pathways; then, we introduced the common pathways with the common miRNAs between COPD and NSCLC; and finally, we analyzed the enriched miRNA-pathway sets by identifying the number of target genes for each miRNAs that contributed in a specific pathway. To have an overall view, the workflow of the different steps is visualized in Fig. 1.
This study presents common miRNA biomarkers between COPD and NSCLC of pre-processed datasets via miRNA-pathway set enrichment analysis and highlights those pathways with more target genes of miRNAs associated with COPD and NSCLC. As such, it specifies the most significant miRNAs or core miRNAs using analyzed pathways. In addition, it assesses the most significant pathways by affecting the core miRNAs on their targets as the components in the pathways. In the meanwhile and as a final step, this study has performed a literature-based search to study the identified miRNA biomarkers on the common pathways.
To construct the expression matrices for all samples, we removed zero values from both datasets. Eventually, the total number of miRNAs after normalization was equal to 1308 and 1145 miRNAs for COPD and NSCLC, respectively; which were considered for further analysis.
The workflow of steps performed in this study. This scheme shows that after collecting miRNA expression profiles, pre-processing was individually performed for each dataset, and then, the enrichment miRNA-pathways were utilized to discover dysregulated pathways though miRNA sets. Those common miRNAs that had the most effects on the enriched pathways on the basis of enrichment scores were selected, and the target genes were extracted from target prediction databases for common miRNAs between COPD and NSCLC. At the end, the pathways analysis was performed.
Enrichment analysis and identification of dysregulated pathways
The results of miRNA set enrichment analysis revealed the pathways regulated by each miRNA in each disease. We identified 149 significant enriched pathways in COPD (1 up-regulated and 148 down-regulated pathways) and 146 significant enriched pathways in NSCLC (72 up-regulated and 74 down-regulated pathways). Among all enriched pathways, similar pathways were found between down-regulated pathways in COPD and up-regulated pathways in NSCLC. In Tables 1 and 2, we only demonstrated the top 10 significant enriched pathways for COPD and NSCLC, respectively. In these tables, size of pathways based on the number of contributed features (SIZE), pathways’ enrichment scores before and after running enrichment peak (ES and NES), percentage of miRNA list before running enrichment peak (Mir%), and the enrichment signal strength are represented in the columns. The full list of common pathways between both diseases is shown in Table S2 and S3 for COPD and NSCLC, respectively.
By comparing the enrichment results, we selected 7 common dysregulated pathways with different regulations in COPD and NSCLC, including non-small cell lung cancer, cell cycle, P53 signaling pathway, VEGF signaling pathway, TGF beta signaling pathway, WNT signaling pathway, and ERBB signaling pathway.
Common miRNAs between COPD and NSCLC
Table S1 shows the enriched pathways and common core miRNAs between COPD and NSCLC in each pathway in such a way that these miRNAs were at least common between the two pathways. For better recognition, in Table S1, each miRNA is highlighted with a color scale from Green to Yellow to show well the degree of the replication of the miRNAs in all enriched pathways.
Moreover, to detect significant miRNAs among all pathways, the average enrichment scores of each miRNA for all enriched pathways as well as the mean score of core miRNAs within each pathway were calculated and shown in Table 4. The zero value in each cell means that the miRNA was not found in that pathway.
A network of common pathways is also shown in Fig. 2, in which each node in this network represents the pathway and each edge between two nodes indicates that there are common miRNAs between two pathways.
Next, since we aimed to clarify the significant miRNAs in common pathways between COPD and NSCLC, we selected the most significant enriched pathways based on two factors: the calculated scores (Table 4) and the number of core miRNAs (Fig. 2). Given pathways, including cell cycle, P53 signaling, non-small cell lung cancer, VEGF signaling, ERBB signaling, WNT signaling, and TGF beta signaling pathways in KEGG had the mean enrichment scores: − 0.0628, − 0.0554, − 0.04016, − 0.0306, − 0.0208, − 0.0201, and − 0.0079, respectively. The results showed that the average number of NES in all pathways for COPD that have more pathways than NSCLC were almost equal to − 2, this means that the NES lower than − 2 could be meaningful in biology for COPD. But for NSCLC, the changes of NES were almost stable (0 ∙ 8 ≤ NES ≤ 1 ∙ 1); thus, we considered those pathways with the average of NES less than − 2 for COPD and found the common significant pathways between both diseases. Moreover, the number of core miRNAs in each pathway, as the second factor, was determined to be equal to 15, 15, 11, 10, 9, 7, and 5 for p53 signaling, cell cycle, non-small cell lung cancer, ERBB signaling, WNT signaling, VEGF signaling, and TGF beta signaling pathways, respectively. Finally, the pathways comprising the highest average enrichment scores along with high number of common core miRNAs were selected. Therefore, three pathways including cell cycle, non-small cell lung cancer, and p53 signaling pathways were detected as the most significant pathways. As a further note, the number of shared miRNAs between p53 signaling and cell cycle pathways was 12, the common miRNA numbers between cell cycle and NSCLC pathways was 9, and the number of shared miRNAs between p53 signaling and NSCLC pathways was 8.
Significant miRNAs in the selected enriched pathways
In Fig. 3, the results for the differential expression level of miRNA set (miRNA-pathway) along with weighted miRNA correlation and ranked miRNA list based on miRScores in non-small cell lung cancer, cell cycle, and p53 signaling pathways are visualized.
Through the selected enriched pathways, we determined core miRNAs which were shared among these pathways and also had higher average of enrichment scores than the other miRNAs. In Fig. 4, the most significant common core miRNAs among three selected pathways are depicted; furthermore, the overlapping core miRNAs of cell cycle, p53 signaling, and non-small cell lung cancer pathways are shown as well.
As a result, the has-miR-15b, hsa-miR-106a, has-miR-17, has-miR-103 and has-miR-107 were selected as the most significant miRNAs with 0.094, 0.083, 0.07, 0.037, and 0.022, respectively, to be the mean of their enrichment scores among all pathways. Afterward, the target genes of the most significant miRNAs in each pathway were identified and the targets mapped to selected pathways. The converted pathways are demonstrated in Fig. 5a, b and c.
Until now, diverse functions have been introduced for miRNAs including inflammation, development of airway epithelial cells, stress responses, and formation of pulmonary surfactant; also, miRNAs have critical roles in the disease’s progression . Moreover, COPD as a progressive and incompletely reversible disease leads to chronic inflammation due to considerable dysregulation of the immune system. Some studies have been conducted to reveal the common pathogenesis mechanisms between COPD and NSCLC [24,25,26,27]. Herein, we have investigated common actuated pathways in both diseases and the role of miRNAs in them. Although we have investigated the most significant common biological pathways between both diseases, our main goal has been to identify the common miRNAs through a pathway enrichment analysis method, what has not been employed before.
To this aim, we first identified DEMs within each pathway as miRNA-pathway; and next, we considered the miRNAs which had more target genes in that pathway. In addition, to avoid the false-positive miRNA-target interactions, we considered the experimentally validated miRNA-target interaction using miRNA-target prediction databases. To select the most significant enriched pathways, a threshold was considered based on the number of DEMs that enriched the pathways. After pathway enrichment, we analyzed three pathways, including non-small cell lung cancer, cell cycle, and p53 signaling pathways as the remarkable enriched pathways based on average ES and their common involved miRNAs, including has-miR-15b, hsa-miR-106a, has-miR-17, has-miR-103 and has-miR-107. These pathways were down-regulated in COPD and up-regulated in NSCLC. The aforementioned miRNAs have been demonstrated to be associated with COPD and lung cancer [28,29,30,31,32,33,34,35,36,37]. The point is that though miRNAs’ roles are known in COPD and NSCLC, understanding their functions in common significant pathways may shed light on the pathogenic mechanism of both diseases and develop new therapeutic targets.
MiR-107 is a tumor suppressor that targets the epidermal growth factor receptor (EGFR). The deregulation of EGFR has been observed in multiple types of cancers, including NSCLC, while frequent EGFR protein overexpression was observed in NSCLC and COPD [38, 39]. Not only EGFR facilitates proliferative signaling through downstream signaling pathways, i.e. PI3K/AKT/mTOR and RAS/ERK, EGFR signaling pathway is one of the activated pathways in lung cancer and employs downstream RAS or ERK pathways to direct proliferative signaling for lung cancer cells . Moreover, the down-regulation of miR-107 may cause cell cycle and proliferation due to the up-regulation of CDK6 (also targeted by miR-103) and CDK8 and metastasis and tumor growth because of the up-regulation of BDNF as well as indirect regulating of the P13K/AKT signaling pathway [40, 41]. The specific binding of BDNF to tropomyosin-related receptor kinase B (TrkB) results in the activation of several downstream pathways such as RAS/ERK, PLC/PKC, PI3K/AKT, JAK/STAT, and AMPK/ACC . However, miR-107 was up-regulated in COPD, which might cause down-regulation of its targets leading to down-regulation of downstream pathways of cell cycle and non-small cell lung cancer pathways. Furthermore, miR-107 has a substantial role in the regulation of echinoderm microtubule-associated protein-like 4 (EML4) - anaplastic lymphoma kinase (ALK) fusion - which may result in the constitutive ALK activation, thus facilitating invasion, cell proliferation, and inhibition of apoptosis .
MiR-106a, as a member of the identified common miRNA families in this study, is an oncogenic miRNA which targets transcription factor (TF) of FOXO3 thus modulating the expression of genes implicated in apoptosis, cell cycle arrest, and autophagy. Furthermore, FOXO3a as a target gene of miR-103, can boost metastasis downstream of PI3K/AKT prohibition in collaboration with the WNT/β-catenin pathway in colon cancer. However, the final results of transcriptional activation or inactivation of FOXO are changeable depending on the context wherein that they occur, as it is a tumor suppressor within the context of PI3K/AKT signaling in neuroblastoma . Another TF targeted by miR-106a is E2F2 belonging to the E2F family, which controls the cell cycle as well as the function of tumor suppressor proteins. On the other hand, the AKT3 protein is known as a critical regulator of the PI3K-AKT-mTOR pathway and miR-106a may decrease the activities of the PI3K/AKT pathway by suppressing the transcription function of E2F2 as well as down-regulation of AKT3.
D. Yang et al.  stated that the expression level of miR-103, which we extracted as a common effective miRNA between both diseases mentioned, decreases in NSCLC and COPD tissues, while it inversely correlates with tumor stage and size. Moreover, miR-103 can prohibit cell proliferation, reduce tumor volume and weight, and increase apoptosis in NSCLC, while targeting MAP2K2 which is a member of MAPK signaling cascade and also a RAS downstream signaling pathway regulating cell proliferation, differentiation, and survival by ERK1/2 activation . On the other hand, MAP2K2 as the downstream of proto-oncogene BRAF is also targeted by miR-17 which is up-regulated in NSCLC and COPD. MAP2K2 has a key role in cell proliferation and cell cycle regulation, so several compounds have been developed to inhibit it in various diseases .
MiR-17 also targets TGFA that encodes transforming growth factor-α (TGF-α), a member of the epidermal growth factor family, what causes the activation of a signaling pathway for cell proliferation, differentiation, and development; this is while the down-regulation of TGF-α due to overexpression of miR-17 may inhibit tumorigenesis and disease progression . Furthermore, in , L. Chen et al. defined E2F3 as a transcriptional activator that is a target gene of miR-17 capable of increasing the cellular proliferation via boosting the G1/S transition; besides, the down-regulation of E2F3 might be due to the molecular mechanism employed by up-regulation of miR-17 as a tumor suppressor.
While MiR-15, the other common core miRNA between COPD and NSCLC, may also function as a tumor-suppressor through the down-regulation of PIK3R3, a gene can serve as a second messenger in growth signaling pathways and can induce cell cycle arrest in the G1-G0 phase and act as a tumorigenesis miRNA in NSCLC [49,50,51]. However, T. Yang et al.  reported that miR-15b may suppress cell proliferation and induce apoptosis to inhibit cell survival, that is it can inhibit cell proliferation and invasion through down-regulation of ATK3.
S. Lim and P. Kaldis  investigated the regulation of the cell cycle pathway in the cell growth and stated that cyclin dependent kinases (CDKs), which are the target genes of miR-17, miR-107, and miR-103, are the key regulatory enzymes that regulate the progression of cells. In addition, cyclin-CDK inhibitor (CKI) family members have involved various functions including DNA damage repair, transcription, metabolism, epigenetic regulation, proteolytic degradation, stem cell self-renewal, spermatogenesis, and neuronal functions . Moreover, EM. Gordon et al.  reported that the transcription factors E2Fs and their regulator Rb are the targets of CDKs, while the E2F proteins during the G1 phase of the cell cycle are activated by phosphorylation of Rb by CDK4/cyclin D (cycD) and Cdk2/cyclin E (cycE) complexes. As such, the transcriptional regulator E2F was found to be a crucial transcriptional regulator in cell cycle . According to our enrichment analysis results, E2F is the target of miR-15 and miR-106a, and it seems the dysregulation of CDKs by some regulator miRNAs could affect cell proliferation or division. Moreover, P53 and its transcriptional targets has a critical function in both G1 and G2 checkpoints . Furthermore, the transcription of various miRNAs is activated by p53 resulting in the repression of genes regulating DNA repair, apoptosis, and cell-cycle progression. Moreover, the regulation of p53 by miR-17, miR-103a, and miR-103b in the cell cycle pathway may be effective in apoptosis. Therefore, the effects of these miRNAs on cell cycle can provide new perspectives on the treatment.
This study has investigated the interactions between core miRNAs and target genes within enriched pathways through miRNA-pathway enrichment analysis in an attempt to identify miRNA biomarkers that are common between COPD and NSCLC. In addition to miRNA-target genes in each pathway, there exist other miRNAs that might have interactions with the miRNA-target genes in the pathways; thus, the study of gene co-expression networks between our identified target genes and genes within the pathways might be a very informative site in this study and can be considered as a research field for future work.
In conclusion, this study showed that the identified miRNAs including miR-106a, miR-17, miR-17, miR-15b, miR-107, and miR-103 were dysregulated in COPD and NSCLC. Moreover, the identification of molecular functions and investigating common pathways regulated by the discovered miRNAs and their target genes may help to look into and understand the possible mechanisms at work in such diseases. As such, it also may have the potential to be used for diagnostic as well as therapeutic goals regarding both diseases.
MiRNA expression profiles
In this project, we considered two miRNA expression profiles related to COPD and NSCLC. Both datasets were downloaded from the gene expression omnibus database (https://www.ncbi.nlm.nih.gov/geo/). The first dataset with the accession number GSE38974 contained 70 miRNAs that are differentially expressed between COPD subjects and smokers without COPD. The used platform is a miRCURY LNA microarray (GPL7723) for the miRNA dataset. Accordingly, the differential miRNA expression set between control and COPD samples were identified using a statistical test with the false positive discovery rate (pFDR) lower than 0.05 and with at least 1.5 fold changes between the groups. Finally, the dataset contained the miRNA expression profiles of 27 samples of which 19 samples were related to the COPD patients and 8 samples were related to normal subjects.
The second dataset with the accession number GSE36681 contained miRNA expression profiling of 56 pairs of fresh-frozen (FF) and 47 pairs of formalin-fixed, paraffin-embedded (FFPE) samples taken from never smokers’ lungs. The most differentially expressed miRNAs were evaluated by Cox analysis and Log-Rank test. Then, the functional significance of the most significant miRNAs was then measured via detecting the candidate targets experimentally. The Illumina human microRNA expression beadchip was used as the platform for this experiment with GPL8179, what included 206 samples of which again the same numbers of samples, as stated above, were related to NSCLC patients and normal tissues. As such, the GEOquery R package  was used for downloading the expression data.
Normalization and pre-processing of miRNA expression profiles
All expression data were quantile-normalized and log2-transformed in R using EdgeR package . Afterward, the samples were checked to exclude the ones containing the missing data or zero variances. Two expression matrices related to COPD and NSCLC cases were reconstructed.
miRNA set enrichment analysis
In the next step, both datasets were enriched using the MiRSEA package in R . This package was utilized to pathway enrichment analysis of differential expressed miRNAs (DEMs) using KEGG pathways. It is to be mentioned that MiRSEA determines the miRNAs regulating pathways and calculates miRNA-pathway weights based on the hypergeometric test (eq. (1)).
In this equation, Wij denotes the weight of association between miRNA i and pathway j, and p is measured as eq. 2.
Where m denotes the number of genes in the whole genome; t is the number of genes involved in pathway j; n is the number of targets of miRNA i; r denotes the number of overlaps between targets of miRNA i and genes in pathway j.
MiRSEA determines DEMs between the two phenotypes considering FDR < 1 and it thus carries out enrichment analysis by comparing DEMs with the miRNAs list in various pathways. Following this, it combines the differential expression levels of the miRNAs and the miRNA-pathway weight (Wij), and defines a miRScore for each enriched miRNA-pathway as eq. (3).
Where Wi is the weight of miRNA i with a given pathway and DEi is the differential expression level of miRNA i.
Thus, a miRNA in a pathway with miRScore greater than zero indicates that the miRNA would probably regulate the pathway in the specific phenotype. MiRSEA ranks miRNAs in the profile and forms miRNA list according to the decreasing miRScore. We selected core miRNAs with the highest miRScores in each pathway for COPD and NSCLC, each separately, as these core miRNAs may have key functions in their pathways through regulating their target genes.
Discovering dysregulated pathways and common miRNAs between COPD and NSCLC
To identify dysregulated KEGG pathways, we first sorted all enriched pathways based on miRNA Enrichment Score (miRES) that shows the extent of overrepresentation of pathways toward top or bottom of the ranked miRNA list. For COPD and NSCLC, we mapped both ranked list of dysregulated pathways and selected miRNAs which either were common or had higher miRESs. Then, we identified canonical pathways associated with a specific phenotype. We discovered the regulated pathways by miRNA set that were common between COPD and NSCLC and also had the differential expression of miRNAs among the two phenotypes. Finally, we selected the most significant pathways and related miRNAs (miRNA-pathways) based on ES.
For each pathway in both diseases, we evaluated the miRNA-pathways to determine which miRNAs regulated the pathway with more targets. We determined the correlated miRNAs within each pathway with a differential weighted score (dw-score) based on eq. (1) for each disease, separately. Among these miRNAs, we specified core miRNAs at and before the point where miRSEA is acquired (miRSEA(p) < 0) and then selected common miRNAs between COPD and NSCLC. After identifying common core miRNAs between the two diseases, we found those miRNAs that were common among all selected pathways, and created two lists of miRNAs (COPD and NSCLC cases) for each enriched pathway. We then combined both miRNA lists related to the diseases for each common enriched pathway in order to preserve only joint core miRNAs in each list. Finally, we calculated the mean of enrichment scores for each pathway and reconstructed a list of common miRNAs among all enriched pathways. In order to select the most significant common miRNAs among all enriched pathways, we selected the miRNAs based on the highest of average enrichment score found in all enriched pathways.
Predicting miRNA targets and analyzing significant common pathways
The target genes of the selected miRNAs were identified by MiRSEA through four target genes prediction databases, i.e. miRWalk, TarBase, miRTarBase, and miR2Disease. To better understand the regulation mechanisms of these common miRNAs within the enriched pathways, we mapped these targets into three numbers of the most common pathways. In addition, to visualize miRNAs and their target genes in a pathway, we used WikiPathway  and Pathvisio  aiming to map and analyze the miRNAs within pathways.
Availability of data and materials
The datasets used and/or analyzed during the current study are available at GEO database: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE38974 and https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE36681. Programming language: R. Other requirements: R environment. R. Packages: GEOquery, EdgR, and MiRSEA. Tested on R version 3.6.1.
Anaplastic lymphoma kinase
Adenosine monophosphate-activated protein kinase
Brain-derived neurotrophic factor
Chronic obstructive pulmonary disease
Differential expressed miRNA
Epidermal growth factor receptor
Echinoderm microtubule-associated protein-like 4
Extracellular signal-regulated protein kinase
False discovery rate
Fork head box O3
Fragments per kilo base of transcript per million mapped reads
Gene expression omnibus
Dual specificity mitogen-activated protein kinase kinase 2
Mitogen-activated protein kinase
Percentage of miRNA list before running enrichment peak
miRNA Enrichment Score
MiRNA set enrichment analysis
Non-small cell lung cancer
Phosphatidylinositol 3-kinase regulatory subunit gamma
Protein kinase B
Pearson correlation coefficient
Small cell lung cancer
Transforming growth factor
Tropomyosin-related receptor kinase B
Vascular endothelial growth factor
Vogelmeier CF, Criner GJ, Martinez FJ, Anzueto A, Barnes PJ, Bourbeau J, et al. Global strategy for the diagnosis, management, and prevention of chronic obstructive lung disease 2017 report. Am J Respir Crit Care Med. 2017;195(5):557–82. https://doi.org/10.1164/rccm.201701-0218PP.
Butler SJ, Ellerton L, Goldstein RS, Brooks D. Prevalence of lung cancer in chronic obstructive pulmonary disease: a systematic review. Respir Med X. 2019;1:100003. https://doi.org/10.1016/j.yrmex.2019.100003.
Durham AL, Adcock IM. The relationship between COPD and lung cancer. Lung Cancer. 2015;90(2):121–7. https://doi.org/10.1016/j.lungcan.2015.08.017.
Bade BC, Dela Cruz CS. Lung Cancer 2020. Clin Chest Med. 2020;41(1):1–24. https://doi.org/10.1016/j.ccm.2019.10.001.
Park HY, Kang D, Shin SH, Yoo K-H, Rhee CK, Suh GY, et al. Chronic obstructive pulmonary disease and lung cancer incidence in never smokers: a cohort study. Thorax. 2020;75(6):506–9. https://doi.org/10.1136/thoraxjnl-2019-213732.
Van Roosbroeck K, Calin GA. Cancer hallmarks and MicroRNAs: the therapeutic connection. Adv Cancer Res. 2017;(135):119-49. https://doi.org/10.1016/bs.acr.2017.06.002.
Chang JT, Anic GM, Rostron BL, Tanwar M, Chang CM. Cigarette smoking reduction and health risks: a systematic review and Meta-analysis. Nicotine Tob Res. 2020;23(4):635–42. https://doi.org/10.1093/ntr/ntaa156.
Parris BA, O’Farrell HE, Fong KM, Yang IA. Chronic obstructive pulmonary disease (COPD) and lung cancer: common pathways for pathogenesis. J Thorac Dis. 2019;11:S2155–72. https://doi.org/10.21037/jtd.2019.10.54.
Reck M, Heigener DF, Mok T, Soria JC, Rabe KF. Management of non-small-cell lung cancer: recent developments. Lancet. 2013;382(9893):709–19. https://doi.org/10.1016/S0140-6736(13)61502-0.
Yao H, Rahman I. Current concepts on the role of inflammation in COPD and lung cancer. Curr Opin Pharmacol. 2009;9(4):375–83. https://doi.org/10.1016/j.coph.2009.06.009.
Tessema M, Tassew DD, Yingling CM, Do K, Picchi MA, Wu G, et al. Identification of novel epigenetic abnormalities as sputum biomarkers for lung cancer risk among smokers and COPD patients. Lung Cancer. 2020;146:189–96. https://doi.org/10.1016/j.lungcan.2020.05.017.
Bartel DP. MicroRNAs: Genomics, Biogenesis, Mechanism, and Function. Cell. 2004;116(2):281-97. https://doi.org/10.1016/S0092-8674(04)00045-5.
Ali Syeda Z, Langden SSS, Munkhzul C, Lee M, Song SJ. Regulatory mechanism of MicroRNA expression in Cancer. Int J Mol Sci. 2020;21(5):1723. https://doi.org/10.3390/ijms21051723.
Okubo M, Tahara T, Shibata T, Yamashita H, Nakamura M, Yoshioka D, et al. Association study of common genetic variants in pre-microRNAs in patients with ulcerative colitis. J Clin Immunol. 2011;31(1):69–73. https://doi.org/10.1007/s10875-010-9461-y.
Liao J, Shen J, Leng Q, Qin M, Zhan M, Jiang F. MicroRNA-based biomarkers for diagnosis of non-small cell lung cancer (NSCLC). Thorac Cancer. 2020;11(3):762–8. https://doi.org/10.1111/1759-7714.13337.
Hamzeiy H, Suluyayla R, Brinkrolf C, Janowski SJ, Hofestaedt R, Allmer J. Visualization and analysis of MicroRNAs within KEGG pathways using VANESA. J Integr Bioinform. 2017;14:9. https://doi.org/10.1515/jib-2016-0004.
Yang IA, Relan V, Wright CM, Davidson MR, Sriram KB, Savarimuthu Francis SM, et al. Common pathogenic mechanisms and pathways in the development of COPD and lung cancer. Expert Opin Ther Targets. 2011;15(4):439–56. https://doi.org/10.1517/14728222.2011.555400.
Chan SMH, Selemidis S, Bozinovski S, Vlahos R. Pathobiological mechanisms underlying metabolic syndrome (MetS) in chronic obstructive pulmonary disease (COPD): clinical significance and therapeutic strategies. Pharmacol Ther. 2019;198:160–88. https://doi.org/10.1016/j.pharmthera.2019.02.013.
Ibuki Y, Toyooka T, Zhao X, Yoshida I. Cigarette sidestream smoke induces histone H3 phosphorylation via JNK and PI3K/Akt pathways, leading to the expression of proto-oncogenes. Carcinogenesis. 2014;35(6):1228–37. https://doi.org/10.1093/carcin/bgt492.
Yun CL, Zierath JR. AMP-activated protein kinase signaling in metabolic regulation. J Clin Invest. 2006;100:328–41. https://doi.org/10.1161/01.
Pian C, Zhang G, Gao L, Fan X, Li F. miR+pathway: the integration and visualization of miRNA and KEGG pathways. Brief Bioinform. 2020;21(2):699–708. https://doi.org/10.1093/bib/bby128.
Brinkrolf C, Janowski SJ, Kormeier B, Lewinski M, Hippe K, Borck D, et al. VANESA - a software application for the visualization and analysis of networks in system biology applications. J Integr Bioinform. 2014;11(2):239. https://doi.org/10.2390/biecoll-jib-2014-239.
Akbas F, Coskunpinar E, Aynaci E, Müsteri Oltulu Y, Yildiz P. Analysis of Serum Micro-RNAs as potential biomarker in chronic obstructive pulmonary disease. Exp Lung Res. 2012;38(6):286-94. https://doi.org/10.3109/01902148.2012.689088.
Powell HA, Iyen-Omofoman B, Baldwin DR, Hubbard RB, Tata LJ. Chronic obstructive pulmonary disease and risk of lung cancer: the importance of smoking and timing of diagnosis. J Thorac Oncol. 2013;8(1):6–11. https://doi.org/10.1097/JTO.0b013e318274a7dc.
Gao Y-H, Guan W-J, Liu Q, Wang H-Q, Zhu Y-N, Chen R-C, et al. Impact of COPD and emphysema on survival of patients with lung cancer: a meta-analysis of observational studies. Respirology. 2016;21(2):269–79. https://doi.org/10.1111/resp.12661.
Sanchez-Salcedo P, Zulueta JJ. Lung cancer in chronic obstructive pulmonary disease patients, it is not just the cigarette smoke. Curr Opin Pulm Med. 2016;22(4):344–9. https://doi.org/10.1097/MCP.0000000000000283.
Barreiro E, Bustamante V, Curull V, Gea J, López-Campos JL, Muñoz X. Relationships between chronic obstructive pulmonary disease and lung cancer: biological insights. J Thor Dis. 2016;8(10):E1122–35. https://doi.org/10.21037/jtd.2016.09.54.
Sato T, Baskoro H, Rennard SI, Seyama K, Takahashi K. MicroRNAs as therapeutic targets in lung disease: prospects and challenges. Chronic Obstr Pulm Dis. 2016;3(1):382–8. https://doi.org/10.15326/jcopdf.3.1.2015.0160.
Wu K, Tsai Y, Lien C, Kuo P. The Roles of MicroRNA in Lung Cancer; 2019. p. 1–25.
Pietrusińska M, Pająk A, Górski P, Kuna P, Szemraj J, Goździńska-Nielepkowicz A, et al. Preliminary studies: differences in microRNA expression in asthma and chronic obstructive pulmonary disease. Adv Dermatol Allergol. 2016;4(4):276–80. https://doi.org/10.5114/ada.2016.61603.
Wang R, Xu J, Liu H, Zhao Z. Peripheral leukocyte microRNAs as novel biomarkers for COPD. Int J COPD. 2017;Volume 12:1101–12. https://doi.org/10.2147/COPD.S130416.
Soeda S, Ohyashiki JH, Ohtsuki K, Umezu T, Setoguchi Y, Ohyashiki K. Clinical relevance of plasma miR-106b levels in patients with chronic obstructive pulmonary disease. Int J Mol Med. 2013;31(3):533–9. https://doi.org/10.3892/ijmm.2013.1251.
Ezzie ME, Crawford M, Cho J-H, Orellana R, Zhang S, Gelinas R, et al. Gene expression networks in COPD: microRNA and mRNA regulation. Thorax. 2012;67(2):122–31. https://doi.org/10.1136/thoraxjnl-2011-200089.
Xie X, Liu HT, Mei J, Ding FB, Xiao HB, Hu FQ, et al. miR-106a promotes growth and metastasis of non-small cell lung cancer by targeting PTEN. Int J Clin Exp Pathol. 2015;8(4):3827-34. https://doi.org/10.3892/ol.2018.9697.
Yanaihara N, Caplen N, Bowman E, Seike M, Kumamoto K, Yi M, et al. Unique microRNA molecular profiles in lung cancer diagnosis and prognosis. Cancer Cell. 2006;9(3):189-98. https://doi.org/10.1016/j.ccr.2006.01.025.
Eapen MS, McAlinden KD, Myers S, Lu W, Sohal SS. microRNAs are key regulators in Chronic Lung Disease: Exploring the vital link between disease progression and Lung Cancer. J Clin Med. 2019;8(11):1986. https://doi.org/10.3390/jcm8111986.
Denkçeken T, Pala E. Investigation of key miRNAs and potential mechanisms in non-small cell lung cancer development from chronic obstructive pulmonary disease. Gen Physiol Biophys. 2020;39(01):69–77. https://doi.org/10.4149/gpb_2019042.
Brambilla E, Gazdar A. Pathogenesis of lung cancer signalling pathways: roadmap for therapies. Eur Respir J. 2009;33(6):1485–97. https://doi.org/10.1183/09031936.00014009.
Suzuki M, Wada H, Yoshino M, Tian L, Shigematsu H, Suzuki H, et al. Molecular characterization of chronic obstructive pulmonary disease-related non-small cell lung cancer through aberrant methylation and alterations of EGFR signaling. Ann Surg Oncol. 2010;17(3):878-88. https://doi.org/10.1245/s10434-009-0739-3.
Xia H, Li Y, Lv X. MicroRNA-107 inhibits tumor growth and metastasis by targeting the BDNF-mediated PI3K/AKT pathway in human non-small lung cancer. Int J Oncol. 2016;49(4):1325–33. https://doi.org/10.3892/ijo.2016.3628.
Wang P, Liu X, Shao Y, Wang H, Liang C, Han B, et al. MicroRNA-107-5p suppresses non-small cell lung cancer by directly targeting oncogene epidermal growth factor receptor. Oncotarget. 2017;8:57012–23. https://doi.org/10.18632/oncotarget.18505.
Takahashi Y, Forrest ARR, Maeno E, Hashimoto T, Daub CO, Yasuda J. MiR-107 and MiR-185 can induce cell cycle arrest in human non small cell lung cancer cell lines. PLoS One. 2009;4(8):e6677. https://doi.org/10.1371/journal.pone.0006677.
Santo EE, Stroeken P, Sluis PV, Koster J, Versteeg R, Westerhout EM. FOXO3a is a major target of inactivation by PI3K/AKT signaling in aggressive neuroblastoma. Cancer Res. 2013;73(7):2189–98. https://doi.org/10.1158/0008-5472.CAN-12-3767.
Yang D, Wang JJ, Li JS, Xu QY. MiR-103 functions as a tumor suppressor by directly targeting programmed cell death 10 in NSCLC. Oncol Res. 2018;26(4):519–28. https://doi.org/10.3727/096504017X15000757094686.
Scheffler M, Holzem A, Kron A, Nogova L, Ihle MA, von Levetzow C, et al. Co-occurrence of targetable mutations in non-small cell lung cancer (NSCLC) patients harboring MAP2K1 mutations. Lung Cancer. 2020;144:40–8. https://doi.org/10.1016/j.lungcan.2020.04.020.
Li W, Cui K, Prochownik EV, Li Y. The deubiquitinase USP21 stabilizes MEK2 to promote tumor growth. Cell Death Dis. 2018;9(5):482-94. https://doi.org/10.1038/s41419-018-0523-z.
Chen S, Sun KX, Liu BL, Zong ZH, Zhao Y. MicroRNA-505 functions as a tumor suppressor in endometrial cancer by targeting TGF-aα. Mol Cancer. 2016;15(1):11. https://doi.org/10.1186/s12943-016-0496-4.
Chen L, Kong G, Zhang C, Dong H, Yang C, Song G, et al. MicroRNA-432 functions as a tumor suppressor gene through targeting E2F3 and AXL in lung adenocarcinoma. Oncotarget. 2016;7(15):20041–53. https://doi.org/10.18632/oncotarget.7884.
Bandi N, Vassella E. miR-34a and miR-15a/16 are co-regulated in non-small cell lung cancer and control cell cycle progression in a synergistic and Rb-dependent manner. Mol Cancer. 2011;10:55. https://doi.org/10.1186/1476-4598-10-55.
Alderman C, Sehlaoui A, Xiao Z, Yang Y. MicroRNA-15a inhibits the growth and invasiveness of malignant melanoma and directly targets on CDCA4 gene. Tumor Biol. 2016;37(10):13941–50. https://doi.org/10.1007/s13277-016-5271-z.
Li B, Xie Z, Li B. miR-152 functions as a tumor suppressor in colorectal cancer by targeting PIK3R3. Tumor Biol. 2016;37(8):10075-84. https://doi.org/10.1007/s13277-016-4888-2.
Yang T, Thakur A, Chen T, Yang L, Lei G, Liang Y, et al. MicroRNA-15a induces cell apoptosis and inhibits metastasis by targeting BCL2L2 in non-small cell lung cancer. Tumor Biol. 2015;36(6):4357-65. https://doi.org/10.1007/s13277-015-3075-1.
Yamakuchi M, Lotterman CD, Bao C, Hruban RH, Karim B, Mendell JT, et al. P53-induced microRNA-107 inhibits HIF-1 and tumor angiogenesis. Proc Natl Acad Sci U S A. 2010;107(14):6334–9. https://doi.org/10.1073/pnas.0911082107.
Lim S, Kaldis P. Cdks, cyclins and CKIs: roles beyond cell cycle regulation. Development. 2013;140(15):3079–93. https://doi.org/10.1242/dev.091744.
Gordon E, Ravicz J, Liu S, Chawla S, Hall F. Cell cycle checkpoint control: the cyclin G1/Mdm2/p53 axis emerges as a strategic target for broad-spectrum cancer gene therapy - a review of molecular mechanisms for oncologists. Mol Clin Oncol. 2018. https://doi.org/10.3892/mco.2018.1657.
Sharma P, Alsharif S, Bursch K, Parvathaneni S, Anastasakis DG, Chahine J, et al. Keratin 19 regulates cell cycle pathway and sensitivity of breast cancer cells to CDK inhibitors. Sci Rep. 2019;9(1):14650. https://doi.org/10.1038/s41598-019-51195-9.
Sean D, Meltzer PS. GEOquery: a bridge between the gene expression omnibus (GEO) and BioConductor. Bioinformatics. 2007;23(14):1846–7. https://doi.org/10.1093/bioinformatics/btm254.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2009;26(1):139–40. https://doi.org/10.1093/bioinformatics/btp616.
Han J, Liu S, Zhang Y, Xu Y, Jiang Y, Zhang C, et al. MiRSEA: Discovering the pathways regulated by dysfunctional MicroRNAs. Oncotarget. 2016;7(34):55012-25. https://doi.org/10.18632/oncotarget.10839.
Pico AR, Kelder T, van Iersel MP, Hanspers K, Conklin BR, Evelo C. WikiPathways: pathway editing for the people. PLoS Biol. 2008;6(7):e184. https://doi.org/10.1371/journal.pbio.0060184.
Kutmon M, van Iersel MP, Bohler A, Kelder T, Nunes N, Pico AR, et al. PathVisio 3: an extendable pathway analysis toolbox. PLoS Comput Biol. 2015;11(2):e1004085. https://doi.org/10.1371/journal.pcbi.1004085.
There were no sources of funding for the research.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Common core miRNAs among all enriched pathways. In addition, all miRNAs are depicted with color scales from Green for more replicated miRNAs to Yellow for less replicated miRNAs. For example, hsa-miR-107 is common between five pathways: cell cycle, ERBB signaling, p53 signaling, VEGF signaling, and non-small cell lung cancer pathways, thus is highlighted with dark green, or hsa-miR-203 is shared between two pathways: cell cycle and non -small cell lung cancer pathways which is specified with light yellow. Table S2. Down-regulated enriched pathways in COPD. Also, the size of pathways based on the number of contributed features (SIZE), pathways’ enrichment scores before and after running enrichment peak (ES and NES), percentage of miRNA list before running enrichment peak (Mir%), and enrichment signal strength are represented in the columns. Moreover, the strength of NESs for all pathways is depicted by color-scaled column, which means that the red NES is more meaningful pathway in biology than the green one. Table S3. Up-regulated enriched pathways in NSCLC. Also, the size of pathways based on the number of contributed features (SIZE), pathways’ enrichment scores before and after running enrichment peak (ES and NES), percentage of miRNA list before running enrichment peak (Mir%), and enrichment signal strength are represented in the columns. Moreover, the strength of NESs for all pathways is depicted by color-scaled column, which means that the red NES is more meaningful pathway in biology than the green one.
About this article
Cite this article
Fathinavid, A., Ghobadi, M.Z., Najafi, A. et al. Identification of common microRNA between COPD and non-small cell lung cancer through pathway enrichment analysis. BMC Genom Data 22, 41 (2021). https://doi.org/10.1186/s12863-021-00986-z