Skip to main content
Download PDF

Identification of key genes in chronic intermittent hypoxia-induced lung cancer progression based on transcriptome sequencing

BMC Cancer volume 24, Article number: 41 (2024) Cite this article

Abstract

Background

Obstructive sleep apnea (OSA) is associated with increased risk of lung cancer mortality. Nevertheless, little is known about the underlying molecular mechanisms. This research aimed to investigate differentially expressed genes (DEGs) and explore their function in Lewis lung carcinoma (LLC)-bearing mice exposed to chronic intermittent hypoxia (CIH) by transcriptome sequencing.

Methods

Lung cancer tissues in LLC-bearing mice exposed to CIH or normoxia were subjected for transcriptome sequencing to examine DEGs. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses were employed to explore the function of DEGs. To evaluate the prognostic value of DEGs, the Kaplan–Meier survival analysis in combination with Cox proportional hazard model were applied based on The Cancer Genome Atlas.

Results

A total of 388 genes with 207 up-regulated and 181 down-regulated genes were differentially expressed between the CIH and normoxia control groups. Bioinformatics analysis revealed that the DEGs were related to various signaling pathways such as chemokine signaling pathway, IL-17 signaling pathway, TGF-β signaling pathway, transcriptional misregulation in cancer, natural killer cell mediated cytotoxicity, PPAR signaling pathway. In addition, the DEGs including APOL1, ETFB, KLK8, PPP1R3G, PRL, SPTA1, PLA2G3, PCP4L1, NINJ2, MIR186, and KLRG1 were proven to be significantly correlated with poorer overall survival in lung adenocarcinoma.

Conclusions

CIH caused a significant change of gene expression profiling in LLC-bearing mice. The DEGs were found to be involved in various physiological and pathological processes and correlated with poorer prognosis in lung cancer.

Peer Review reports

Introduction

Obstructive sleep apnea (OSA) is a highly prevalent sleep disorder characterized by intermittent partial or complete collapse, causing chronic intermittent hypoxia (CIH), sleep fragmentation, and increased inspiratory efforts. These can induce a complex series of pathophysiological changes, leading to the damage of multi-organ and multi-system. OSA has been recognized as an independent risk factor for incident cardiovascular disease [1], cognitive impairment [2], and metabolic disease [3]. Accumulating evidence shows that OSA is also independently associated with increased cancer incidence and mortality [4, 5].

Lung cancer is reported to be the most common cause of cancer mortality worldwide, with an estimated 1.8 million deaths [6]. The associations between OSA and occurrence and mortality of lung cancer were supported by more and more studies [7,8,9,10,11]. For instance, a cross-sectional study of 302 subjects demonstrated that sleep apnea and its related nocturnal hypoxia were independently related with an increased prevalence of lung cancer [7]. In addition, Huang et al. [8] found that severe OSA was an risk factor of cancer mortality in stage III-IV lung cancer patients. However, the mechanisms whereby OSA results in increased risks of lung cancer incidence and mortality remain unclear.

In recent decades, RNA sequencing has become an established and powerful method for analyzing differential gene expression and differential splicing of mRNAs. RNA sequencing has been used to uncover molecular mechanisms and explore potential diagnostic and therapeutic targets for many kinds of diseases. Model of CIH has generally been used so as to mimic human OSA [12]. The present study was carried out to identify differentially expressed genes (DEGs) and explore their function in Lewis lung carcinoma (LLC)-bearing mice exposed to CIH by transcriptome sequencing. In addition, the connections between different expression levels of DEGs and prognosis in lung adenocarcinoma were also investigated.

Methods

Animal

Seven-week-old male C57BL/6 mice were obtained from Chinese Academy of Science Laboratory Animals Center located in Shanghai. Mice were allocated to either the normoxia control (NC) (n = 12) or the CIH group (n = 12). All mice were housed in standard cages with a 12-h-day/12-h-night cycle and kept in an animal room. They were free access to water and food. This protocol was approved by the IACUC and IBC Committee in Zhongshan Hospital, Xiamen University (approved number: 2017–015) and performed according to the Guide for the Care and Use of Laboratory Animal. This study is reported in accordance with ARRIVE guidelines.

Intermittent hypoxia (IH) exposure

Mice in CIH group were exposed to IH for 5 consecutive weeks. IH exposure was conducted in the light-on period, from 9 AM to 5 PM. The protocol of IH has been reported in our previous study [13]. The one cycle of IH consisted of the 50 s of nitrogen, lasting for 10 s, oxygen for 10 s, and subsequent compressed air for 50 s. The range of oxygen concentration was between 21% ± 1% to 6% ± 1%.

Lung cancer cell culture and tumor induction

LLC cells were obtained from CoBioer Biosciences Co., Ltd. China. The cells were cultured in high-glucose Dulbecco's Modification of Eagle's Medium in combination with 10% fetal bovine serum (GIBCO, USA). The right flank of each mouse was injected with the LLC cells at density of 1 × 106 LLC/100 μL PBS after 7 days of IH treatment. The tumor volume was recorded every 5 days after a tumor became palpable. The tumor width (W) and length (L) were collected and used for calculating the tumor volume (mm3) (V = W2 × L/2).

Tissue preparation

All the mice were euthanized with intraperitoneal injection of pentobarbital (150 mg/kg) after 5 weeks of the IH treatment. The tumors were excised, weighted and subsequently stored in RNA locker. 6 lung cancer tissue samples were randomly selected from each group and then subjected for sequencing.

RNA isolation, library preparation, and transcriptome sequencing

The total RNA was extracted from the lung cancer tissue using TRIzol® Reagent (Magen). Oligo(dT) magnetic beads was used for mRNA purification. Then the fragmentation was carried out with divalent cations. Subsequently, the double stranded cDNAs were synthesized. PCR amplification was performed with adaptor-ligated cDNA. Furthermore, we purified PCR products. Finally, the library preparations were sequenced on an Illumina Novaseq 6000 platform at Shanghai Applied Protein Technology Co., Ltd. Raw sequence data was deposited in the in the NCBI SRA database repository, accession number PRJNA948556.

Data processing

Initially, we processed raw data of fastq format through in-house perl scripts. In order to obtain mapped reads, the clean reads were separately aligned to reference genome with orientation mode using HISAT2 software (http://daehwankimlab.github.io/hisat2/). We used FeatureCounts (http://subread.sourceforge.net/) to count the reads numbers mapped to each gene. Then fragments per kilobase of transcript per million mapped fragments (FPKM) of each gene was calculated based on reads count mapped to the gene and the length of this gene. Differential expression analysis was performed using the DESeq2 (http://bioconductor.org/packages/release/bioc/html/DESeq2.html). Genes with p value < 0.05 and | log2FC |> 1 were judged to be DEGs.

Enrichment analysis

The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis can explain the functional enrichment of DEGs and clarify the differences between samples at the gene function level. TopGO package (http://www.bioconductor.org/packages/release/bioc/html/topGO.html) for GO function enrichment and KEGG pathway enrichment analysis (http://www.kegg.jp/kegg/kegg1.html) were used in this study. GO or KEGG function was considered to be significantly enriched when p value < 0.05.

qRT-PCR

Six genes were selected to validate the sequencing data with a real-time PCR technology. Real-time PCR reactions were conducted on an ABI 7500 thermocycler (Applied Biosystems, USA). The 2-ΔΔCt method was utilized to analyze the relative gene expressions. The primer sequences are presented in Supplementary Table S1. Data were analyzed using GraphPad Prism 5.0. Data were expressed as mean ± SD. Unpaired Student’s t-test was used for comparisons of qRT-PCR data between the two groups. Statistical significance was determined as p-value < 0.05.

Prognosis analysis

In order to evaluate the prognostic value of DEGs, the Kaplan–Meier survival analysis combined with Cox proportional hazard model were applied based on The Cancer Genome Atlas (TCGA) using R language packages (survival and survminer). R software (version 4.2.1) was used for all statistical analyses.

Results

IH promoted lung cancer growth

The final tumor weight and volume were much higher in the CIH group when compared with NC group. While CIH group mice gained less body weight. The detailed data have been published in our previous research article [13].

Summary of transcriptomic profiles

The raw reads of each library ranged from 40, 131, 160 to 51, 498, 348. The GC content and Q30 values for all libraries were beyond 49.5% and 91.93%, respectively (Supplementary Table S2). The ratio of clean reads to raw reads was over 98%. An average 94.59% of the reads, including multiple mapped reads and unique mapped reads, were mapped to the reference genome. The results suggested that the data were reliable and eligible for the subsequent analysis. The ratio of FPKM value less than 1 is high, which revealed that the majority of genes were expressed at low levels (Supplementary Figure 1A). The values of distribution were similar among the 12 samples (Supplementary Figure 1B). A principal component analysis (PCA) was conducted to check the clustering pattern of the samples. The PC1 and PC2 were 47.23% and 12%, respectively. Although not too far from one another, the distance between CIH and NC was apparent and sufficient for the analysis (Supplementary Figure 1C).

Identification of DEGs

The results presented that there were 388 DEGs, including 207 up-regulated and 181 down-regulated genes, between the CIH group and NC group (Fig. 1A). Then we conducted hierarchical clustering analysis and revealed distinct clustering of CIH samples from NC samples (Fig. 1B). In addition, volcano plots were constructed using p value and fold change value to visualize differential expression between NC and CIH group (Fig. 1C). The top 20 up-regulated and down-regulated DEGs were summarized in Tables 1 and 2, respectively.

Fig. 1
figure 1

DEGs between CIH group and NC group. A Bar graph for the number of up- and down-regulated DEGs. B Heat map generated by hierarchical clustering of DEGs in CIH group and NC group, red color representing increased expressed genes and blue color for decreased expressed genes. C Volcano plot of DEGs between CIH group and NC group, red dots representing the up-regulated DEGs, blue dots representing the down-regulated DEGs

Full size image
Table 1 List of top 20 up-regulated DEGs in CIH group as compared to controls
Full size table
Table 2 List of top 20 down-regulated DEGs in CIH group as compared to controls
Full size table

GO enrichment and KEGG pathway analyses

To better reveal the function of these DEGs, GO function enrichment and KEGG pathway enrichment analyses were conducted. The top 10 GO terms including biological process, molecular function, and cell component for up-regulated and down-regulated DEGs were displayed in Fig. 2A and B, respectively. We found that these DEGs were involved in various molecular functions and biological processes. The up-regulated DEGs were mainly enriched in pathways such as IL-17 signaling pathway, TGF-β signaling pathway, chemokine signaling pathway, transcriptional misregulation in cancer (Fig. 3A). The main KEGG pathways for down-regulated DEGs were natural killer cell mediated cytotoxicity, ABC transporters, and PPAR signaling pathway, among others (Fig. 3B).

Fig. 2
figure 2

GO analyses of DEGs. A GO terms for the up-regulated DEGs. B GO terms for the down-regulated DEGs. Y-axis represents GO terms, X-axis represents the p value

Full size image
Fig. 3
figure 3

KEGG pathway analyses of DEGs. A Top 20 pathways associated with up-regulated DEGs (B) Top 20 pathways related to down-regulated DEGs. The size of the dot represents the number of genes, the color of the dot represents risk factor, Y-axis shows pathway name, X-axis shows the p value

Full size image

Prognosis analysis

The prognostic value of the DEGs in lung adenocarcinoma patients was examined using Kaplan–Meier survival curves based on TCGA datasets. Univariate Cox regression analysis suggested that the high expression levels of up-regulated DEGs including APOL1, ETFB, KLK8, PPP1R3G, PRL, and SPTA1 were significantly correlated with poorer overall survival (OS) in lung adenocarcinoma. The low expression levels of down-regulated DEGs PLA2G3, PCP4L1, NINJ2, MIR186, and KLRG1 were also statistically correlated with a shorter OS in lung adenocarcinoma (Fig. 4).

Fig. 4
figure 4

The overall survivals of lung adenocarcinoma patients with high or low gene expression. A-F The up-regulated DEGs. (G-K) The low-regulated DEGs

Full size image

Validation by qRT-PCR

To verify the accuracy of the RNA sequencing data, 6 prognosis-related genes including 3 up-regulated and 3 down-regulated genes were selected for validation experiments. The results indicated that the expression levels of ETFB, KLK8, and SPTA1 were remarkably up-regulated, and the expression levels of KLRG1, PCP4L1 and NINJ2 were significantly down-regulated (Fig. 5A). This was consistent with sequencing data, confirming the accuracy of sequencing (Fig. 5B).

Fig. 5
figure 5

DEGs validated by qRT-PCR. A The relative expression levels of 6 DGEs by qRT-PCR analysis (n = 3). B The expression pattern of DEGs in both qRT-PCR and RNA-seq. *p < 0.05; **p < 0.01

Full size image

Discussion

This study established models of CIH and lung cancer in mice to explore the potential mechanisms whereby OSA promoted lung cancer progression. 388 genes were differentially expressed in lung cancer between the CIH-treated LLC-bearing mice and NC group. Further bioinformatics analysis found that these DEGs were related to various pathways involving chemokine signaling pathway, IL-17 signaling pathway, TGF-β signaling pathway, transcriptional misregulation in cancer, natural killer cell mediated cytotoxicity, PPAR signaling pathway. In addition, we identified 11 DEGs, which were associated with unfavorable prognosis in lung adenocarcinoma patients.

Growing evidence suggests a close relationship between OSA and cancer incidence and mortality [4, 5]. A community-based cohort including 400 residents with 20 years follow-up reported that moderate-to-severe OSA was independently correlated with an increased risk of cancer incidence and mortality [14]. Nieto et al. [15] demonstrated a dose–response relationship between OSA and cancer mortality in a community-based sample. A multicenter study with 5427 patients showed that OSA severity was associated with increased cancer mortality after follow-up of 4.5 years [16]. As for lung cancer type, a prospective study followed up 65,330 women and suggested that no independent association existed between OSA and overall cancer risk. However, significant associations were noted for the smoking-related cancers (e.g., lung cancer) [11]. A meta-analysis revealed that OSA patients had an nearly 30% higher risk of lung cancer compared with non-OSA subjects after pooling data of four observational studies with 4,885,518 patients [9]. Liu et al. [10] showed that the recurrence rate, metastasis rate and mortality increased in lung cancer patients with OSA during the one-year follow-up period. Li et al. [17] also demonstrated an association between OSA and increased risk for mortality in lung cancer patients. Furthermore, they reveal some molecular pathways involved in sustained hypoxia-induced lung cancer progression some based on Gene Expression Omnibus data base. However, the model of sustained hypoxia did not reflect the pathophysiology of OSA.

These findings have been further supported by animal experiments. In a mouse model of sleep apnea, Almendros and colleagues [18] found that CIH induced a remarkable increase in melanoma lung metastasis when compared to NC group. Another study also demonstrated that CIH promoted lung epithelial TC1 cell tumor growth and invasion toward adjacent tissues in a mouse model of OSA [19]. In addition, our previous animal study showed that CIH accelerated lung cancer development and enhanced the vascular endothelial growth factor expression. The mean SUVmax values assessed by micro-PET–CT were considerably higher in the CIH group than the NC group [20]. Collectively, these results indicated that CIH was an important component of tumor malignant properties, facilitating cancer growth, invasion and migration.

The contributing mechanisms underlying this association are not yet well understood. To gain a better understanding of the mechanisms, Illumina high-throughput technology was utilized to sequence the transcriptome of LLC-bearing mice exposed to CIH. The present study revealed that the up-regulated DEGs were markedly related to chemokine signaling pathway, IL-17 signaling pathway, TGF-β signaling pathway, transcriptional misregulation in cancer. Considerable evidence linked IL-17 with lung cancer. Numasaki et al. [21] demonstrated that IL-17 obviously increased angiogenesis and promoted the growth of lung cancer transplanted in severe combined immunodeficient mice. Another study reported that intraperitoneal injection of IL-17 resulted in remarkably larger tumors in a LLC mouse model when compared with the control group [22]. Furthermore, It was reported that TGF-β could promote tumor progression through suppression of immune surveillance, angiopoiesis, and promotion of epithelial to mesenchymal transition [23, 24]. An increasing body of evidence suggested that the TGF-β pathway activation contributed to poor prognosis in lung cancer patients [25, 26]. The down-regulated DEGs were enriched in natural killer cell mediated cytotoxicity and PPAR signaling pathway. In fact, the killing mediated by natural killer cells and cytotoxic T lymphocytes represents a crucial mechanism in the immune defense against cancers. And the impairment of natural killer cell mediated cytotoxicity facilitates the growth of lung cancer. A clinical study found that severe OSA had considerably fewer invariant natural killer T cells (iNKT) compared to mild-moderate OSA or no OSA patients and 12 months of continuous positive airway pressure therapy increased the frequency of iNKT cells. Furthermore, they found that hypoxia resulted in impaired cytotoxicity [27].

The present study also revealed several DEGs related to poorer prognosis in lung adenocarcinoma patients. Planque et al. [28] analyzed KLK8 mRNAs in 60 NSCLC tissues and in paired unaffected tissues by PCR and found that KLK8-T4 alternative splice variant, alone or in combination was independent marker of poor prognosis in lung cancer. Another study demonstrated that MIR186 could inhibit lung cancer progression through targeting SIRT6 [29]. In addition, a study showed that the deficiency of PRMT5 enhanced Klrg1+ terminal CD8+ T cell development and eliminated antitumor activity [30]. So this indicated that the DEGs identified in this study could be promising candidates for future research in investigating the association between OSA and increased cancer mortality. The identified DEGs collectively suggest a molecular landscape favoring tumor growth, immune evasion, and inflammatory responses in the context of CIH-induced lung cancer. The dysregulation of pathways related to angiogenesis, immune modulation, and transcriptional control may contribute to the observed poorer prognosis in lung adenocarcinoma.

The current study has some strengths. Firstly, high-throughput sequencing was utilized to investigate the mechanisms linking OSA with progression of lung cancer. Secondly, the findings of the key signaling pathways and prognosis related genes provided a new direction for future study on this issue. Thirdly, the lung cancer prognosis related genes were validated in both animal model and human study, which made the result more reliable. Our study also has several limitations that warrant mention. Firstly, the CIH animal models only simulated one of the major features of OSA and lacked of other features (e.g., hypercapnia and sympathetic hyperactivity). Secondly, although the DEGs and their potential function were identified and analyzed, the functional and mechanistic study on cell line were not conducted. Finally, only LLC-bearing mouse model was established. Whether these results can also be successfully applied to other lung cancer type needed further confirmation. Recognizing the heterogeneity of lung cancer, future studies should incorporate different lung cancer subtypes, such as squamous cell carcinoma and small cell carcinoma.

Conclusions

CIH caused a significant change of gene expression profiling in LLC-bearing mice. The dysregulated genes were involved in chemokine signaling pathway, IL-17 signaling pathway, TGF-β signaling pathway, transcriptional misregulation in cancer, natural killer cell mediated cytotoxicity, PPAR signaling pathway. Meanwhile, 11 DEGs were identified to be associated with poorer prognosis in lung adenocarcinoma. The collective impact of these DEGs on key pathways suggested a multi-faceted contribution to CIH-induced lung cancer progression, providing potential insights into the mechanisms underlying the observed poorer prognosis of lung adenocarcinoma in OSA patients.

Availability of data and materials

The data presented in the study are deposited in the NCBI SRA database repository, accession number PRJNA948556.

References

  1. Deng M, Huang YT, Xu JQ, Ke X, Dong YF, Cheng XS. Association Between Intermittent Hypoxia and Left Ventricular Remodeling in Patients With Obstructive Sleep Apnea-Hypopnea Syndrome. Front Physiol. 2020;11: 608347.

    Article  PubMed  Google Scholar 

  2. Gosselin N, Baril AA, Osorio RS, Kaminska M, Carrier J. Obstructive Sleep Apnea and the Risk of Cognitive Decline in Older Adults. Am J Respir Crit Care Med. 2019;199(2):142–8.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Gabryelska A, Karuga FF, Szmyd B, Bialasiewicz P. HIF-1alpha as a Mediator of Insulin Resistance, T2DM, and Its Complications: Potential Links With Obstructive Sleep Apnea. Front Physiol. 2020;11:1035.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Martinez-Garcia MA, Campos-Rodriguez F, Barbe F. Cancer and OSA: Current Evidence From Human Studies. Chest. 2016;150(2):451–63.

    Article  PubMed  Google Scholar 

  5. Gozal D, Farre R, Nieto FJ. Obstructive sleep apnea and cancer: Epidemiologic links and theoretical biological constructs. Sleep Med Rev. 2016;27:43–55.

    Article  PubMed  Google Scholar 

  6. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, Bray F. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J Clin. 2021;71(3):209–49.

    Article  PubMed  Google Scholar 

  7. Seijo LM, Perez-Warnisher MT, Giraldo-Cadavid LF, Oliveros H, Cabezas E, Troncoso MF, Gomez T, Melchor R, Pinillos EJ, El Hachem A, et al. Obstructive sleep apnea and nocturnal hypoxemia are associated with an increased risk of lung cancer. Sleep Med. 2019;63:41–5.

    Article  PubMed  Google Scholar 

  8. Huang HY, Shih-Wei L, Chuang LP, Wang CL, Sun MH, Li HY, Chang CJ, Chang SC, Yang CT, Chen NH. Severe Obstructive Sleep Apnea Associated With Higher Risk of Mortality in Stage III and IV Lung Cancer. J Clin Sleep Med. 2020;16(7):1091–8.

  9. Cheong AJY, Tan BKJ, Teo YH, Tan NKW, Yap DWT, Sia CH, Ong TH, Leow LC, See A, Toh ST. Obstructive Sleep Apnea and Lung Cancer: A Systematic Review and Meta-Analysis. Ann Am Thorac Soc. 2022;19(3):469–75.

    Article  PubMed  Google Scholar 

  10. Liu W, Luo M, Fang YY, Wei S, Zhou L, Liu K. Relationship between Occurrence and Progression of Lung Cancer and Nocturnal Intermittent Hypoxia, Apnea and Daytime Sleepiness. Curr Med Sci. 2019;39(4):568–75.

    Article  PubMed  Google Scholar 

  11. Huang T, Lin BM, Stampfer MJ, Schernhammer ES, Saxena R, Tworoger SS, Redline S. Associations of Self-reported Obstructive Sleep Apnea with Total and Site-specific Cancer Risk in Older Women: A Prospective Study. Sleep. 2021;44(3):zsaa198.

  12. Barnes LA, Mesarwi OA, Sanchez-Azofra A. The Cardiovascular and Metabolic Effects of Chronic Hypoxia in Animal Models: A Mini-Review. Front Physiol. 2022;13: 873522.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Zhang XB, Lin XL, Wu XY, Zeng YM, Chen XY, Luo X, Zeng HQ. Differential expression of microRNAs in xenografted Lewis lung carcinomas subjected to intermittent hypoxia: a next-generation sequence analysis. Transl Cancer Res. 2020;9(7):4354–65.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Marshall NS, Wong KK, Cullen SR, Knuiman MW, Grunstein RR. Sleep apnea and 20-year follow-up for all-cause mortality, stroke, and cancer incidence and mortality in the Busselton Health Study cohort. J Clin Sleep Med. 2014;10(4):355–62.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Nieto FJ, Peppard PE, Young T, Finn L, Hla KM, Farre R. Sleep-disordered breathing and cancer mortality: results from the Wisconsin Sleep Cohort Study. Am J Respir Crit Care Med. 2012;186(2):190–4.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Martinez-Garcia MA, Campos-Rodriguez F, Duran-Cantolla J, de la Pena M, Masdeu MJ, Gonzalez M, Del Campo F, Serra PC, Valero-Sanchez I, Ferrer MJ, et al. Obstructive sleep apnea is associated with cancer mortality in younger patients. Sleep Med. 2014;15(7):742–8.

    Article  PubMed  Google Scholar 

  17. Li L, Lu J, Xue W, Wang L, Zhai Y, Fan Z, Wu G, Fan F, Li J, Zhang C, et al. Target of obstructive sleep apnea syndrome merge lung cancer: based on big data platform. Oncotarget. 2017;8(13):21567–78.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Almendros I, Montserrat JM, Torres M, Dalmases M, Cabanas ML, Campos-Rodriguez F, Navajas D, Farre R. Intermittent hypoxia increases melanoma metastasis to the lung in a mouse model of sleep apnea. Respir Physiol Neurobiol. 2013;186(3):303–7.

    Article  PubMed  Google Scholar 

  19. Almendros I, Wang Y, Becker L, Lennon FE, Zheng J, Coats BR, Schoenfelt KS, Carreras A, Hakim F, Zhang SX, et al. Intermittent hypoxia-induced changes in tumor-associated macrophages and tumor malignancy in a mouse model of sleep apnea. Am J Respir Crit Care Med. 2014;189(5):593–601.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Zhang XB, Yang YY, Zeng Y, Zeng HQ, Fu BB, Ko CY, Luo X, Du YP, Chen LD, Lai YT, et al. Anti-tumor effect of endostatin in a sleep-apnea mouse model with tumor. Clin Transl Oncol. 2019;21(5):572–81.

  21. Numasaki M, Watanabe M, Suzuki T, Takahashi H, Nakamura A, McAllister F, Hishinuma T, Goto J, Lotze MT, Kolls JK, et al. IL-17 enhances the net angiogenic activity and in vivo growth of human non-small cell lung cancer in SCID mice through promoting CXCR-2-dependent angiogenesis. J Immunol. 2005;175(9):6177–89.

    Article  PubMed  Google Scholar 

  22. Wei L, Wang H, Yang F, Ding Q, Zhao J. Interleukin-17 potently increases non-small cell lung cancer growth. Mol Med Rep. 2016;13(2):1673–80.

    Article  PubMed  Google Scholar 

  23. Miyazono K, Ehata S, Koinuma D. Tumor-promoting functions of transforming growth factor-beta in progression of cancer. Upsala J Med Sci. 2012;117(2):143–52.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Heldin CH, Vanlandewijck M, Moustakas A. Regulation of EMT by TGFbeta in cancer. FEBS Lett. 2012;586(14):1959–70.

    Article  PubMed  Google Scholar 

  25. Kong F, Jirtle RL, Huang DH, Clough RW, Anscher MS. Plasma transforming growth factor-beta1 level before radiotherapy correlates with long term outcome of patients with lung carcinoma. Cancer. 1999;86(9):1712–9.

    Article  PubMed  Google Scholar 

  26. Zhao L, Ji W, Zhang L, Ou G, Feng Q, Zhou Z, Lei M, Yang W, Wang L. Changes of circulating transforming growth factor-beta1 level during radiation therapy are correlated with the prognosis of locally advanced non-small cell lung cancer. J Thoracic Oncol. 2010;5(4):521–5.

    Article  Google Scholar 

  27. Gaoatswe G, Kent BD, Corrigan MA, Nolan G, Hogan AE, McNicholas WT, O’Shea D. Invariant Natural Killer T Cell Deficiency and Functional Impairment in Sleep Apnea: Links to Cancer Comorbidity. Sleep. 2015;38(10):1629–34.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Planque C, Choi YH, Guyetant S, Heuze-Vourc’h N, Briollais L, Courty Y. Alternative splicing variant of kallikrein-related peptidase 8 as an independent predictor of unfavorable prognosis in lung cancer. Clin Chem. 2010;56(6):987–97.

    Article  PubMed  Google Scholar 

  29. Ruan L, Chen J, Ruan L, Yang T, Wang P. MicroRNA-186 suppresses lung cancer progression by targeting SIRT6. Cancer Biomarkers. 2018;21(2):415–23.

    Article  PubMed  Google Scholar 

  30. Zheng Y, Chen Z, Zhou B, Chen S, Han L, Chen N, Ma Y, Xie G, Yang J, Nie H, et al. PRMT5 Deficiency Enforces the Transcriptional and Epigenetic Programs of Klrg1(+)CD8(+) Terminal Effector T Cells and Promotes Cancer Development. J Immunol. 2022;208(2):501–13.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by the National Natural Science Foundation of China (82170103), the Natural Science Foundation of Fujian Province (2020J011292), the Young people training project from Fujian Province Health Bureau (2020GGB057), and Xiamen Medical and Health Guidance Project (3502Z20214ZD1043).

Author information

Author notes

  1. Li-Da Chen and Li Lin contributed equally to this work.

Authors and Affiliations

  1. Department of Respiratory and Critical Care Medicine, Zhangzhou Affiliated Hospital of Fujian Medical University, Zhangzhou, Fujian Province, China

    Li-Da Chen, Li Lin, Wei-Liang Zhang, Huang-Yu Li & Jia-Min Luo

  2. Department of Emergency Medicine, Zhangzhou Affiliated Hospital of Fujian Medical University, Zhangzhou, Fujian Province, China

    Ji-Zhi Chen

  3. Ningde Food and Drug Inspection Testing Center, Ningde, Fujian Province, China

    Yang Song

  4. The School of Clinical Medicine, Fujian Medical University, No. 1, Xuefu North Road, University New District, Fuzhou, Fujian Province, 350122, People’s Republic of China

    Xiao-Bin Zhang

  5. Department of Pulmonary and Critical Care Medicine, Zhongshan Hospital, Xiamen University, Xiamen, Fujian Province, China

    Xiao-Bin Zhang

Authors
  1. Li-Da Chen
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Li Lin
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Ji-Zhi Chen
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Yang Song
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Wei-Liang Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Huang-Yu Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Jia-Min Luo
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Xiao-Bin Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

L.-D.C., L.L., X.-B.Z. and J.-Z. C. performed most of the experiments with assistance; Y.S., W.-L.Z. and H.-Y.L contributed to the data analyses; L.-D.C., L.L. and X.-B.Z. wrote the manuscript; J.-M.L. edited the manuscript. All authors have read and agreed to the published version of the manuscript.

Corresponding author

Correspondence to Xiao-Bin Zhang.

Ethics declarations

Ethics approval and consent to participate

The animal study was as approved by the IACUC and IBC Committee in Zhongshan Hospital, Xiamen University (approved number: 2017–015). This study is reported in accordance with ARRIVE guidelines.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

 Supplementary Figure 1. Summary of transcriptomic profiles. (A) The percentage of different FPKM ranges in 12 samples. (B) Boxplot of FPKM distribution among the two groups. (C) Clustering pattern of the samples examined by PCA.

Additional file 2:

 Supplementary Table S1. Primers used for qRT-PCR.

Additional file 3:

 Supplementary Table S2. Quality control.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Chen, LD., Lin, L., Chen, JZ. et al. Identification of key genes in chronic intermittent hypoxia-induced lung cancer progression based on transcriptome sequencing. BMC Cancer 24, 41 (2024). https://doi.org/10.1186/s12885-023-11785-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12885-023-11785-3

Keywords

BMC Cancer

ISSN: 1471-2407

Contact us