FOXM1/lncRNA TYMSOS/miR-214-3p–Mediated High Expression of NCAPG Correlates With Poor Prognosis and Cell Proliferation in Non–Small Cell Lung Carcinoma

Introduction

Lung cancer is the major cause of cancer-related death worldwide, according to cancer statistics in 2020, the morbidity rate of lung cancer ranks second, whereas the death rate of lung cancer ranks the top (Siegel et al., 2020). Lung cancer is composed of non–small cell lung carcinoma (NSCLC) and small cell lung carcinoma, NSCLC includes lung adenocarcinoma (LUAD), lung squamous cell carcinoma, and large-cell lung carcinoma, the NSCLC accounts for approximately 85% of all cases (Molina et al., 2008; Fan et al., 2021). The high lung cancer mortality is owing to cancer metastatic with the time of diagnosis, indicating that ameliorating patient’s survival will require more effective systemic therapies (Morgensztern et al., 2010). However, there is much progress achieved in cancer treatment methods, ranging from targeted therapies to immunotherapy (Herbst et al., 2018); the survival period of lung cancer patients was very short (Xu et al., 2020a; Liang et al., 2021). Therefore, finding a new therapeutic target or biomarker is essential for effectively preventing the development of lung cancer.

Recently, studies have started to characterize the regulatory effects that condensin complexes may have on diverse cancers, including colon adenocarcinoma and liver hepatocellular carcinoma. These condensin complexes play a different role in the progression of cancer, as well as drug resistance (Wu et al., 2019). Studies have indicated that nonstructural maintenance of chromosomes condensin I complex subunit G (NCAPG) was essential for chromatin separation (Eberlein et al., 2009). Its elevated expression has been reported in several cancers, including hepatocellular carcinoma (HCC) (Gong et al., 2019), gastric cancer (Sun et al., 2020), ovarian cancer (Xu et al., 2020b), and cardia adenocarcinoma (Zhang et al., 2021). It has been shown that elevated NCAPG level promotes cancer cell migration by activating Wnt/β-catenin signaling (Zhang et al., 2021). A recent study found that NCAPG through activating the PI3K/AKT pathway promotes HCC cell proliferation and antiapoptosis (Gong et al., 2019). Furthermore, up-regulation of NCAPG expression significantly facilitated the progression of BRCA. It has been shown that NCAPG regulated the stability of STAT3 confers trastuzumab resistance in HER2-positive breast cancer (Jiang et al., 2020). In gastric cancer, Song et al. reported that miR-193b-3p via targets the NCAPG and reduces the expression of NCAPG, resulting in inhibition of the cell proliferation of gastric cancer (Song et al., 2018). Moreover, up-regulated NCAPG expression is also linked to the cancer pathological stage of castration-resistant prostate cancer; knockdown of NCAPG inhibited cancer cell aggressiveness (Arai et al., 2018). However, a comprehensive analysis of the function and regulation mechanism of NCAPG in LUAD is still insufficient. Thus, our study utilized various bioinformatics approaches and functional studies to examine whether NCAPG is involved in NSCLC progression, immune infiltration, and its plausible molecular regulation.

In this study, we compared the expression of NCAPG between NSCLC tissues and normal samples and investigated the correlation between NCAPG expression and clinical parameters of NSCLC. In addition, we explored the prognostic value and clinical significance of NCAPG in NSCLC. Meanwhile, the correlation between NCAPG expression and immune infiltration was analyzed to explore the potential mechanisms involved in NCAPG modulation in the progression of NSCLC. Finally, the biological role of NCAPG was identified in NSCLC. In summary, we demonstrated the potential role of NCAPG in regulating tumor progression and its potential application in the diagnosis and prognostic evaluation in NSCLC.

Materials and Methods

Data Collection

The Cancer Genome Atlas (TCGA)–LUAD cohort data and corresponding clinical information of 535 LUAD patients were downloaded from the TCGA website (https://portal.gdc.cancer.gov/repository). LUAD patients were classified into low- and high-NCAPG expression groups according to the median NCAPG expression value. NCAPG expression data and clinical information from datasets GSE10072 and GSE18842 were downloaded from the GEO database and validated for survival analyses. The gene expression profiles were normalized using the scale method provided in the “limma” R package. Data analysis was performed with the R (version 3.6.3) and ggplot2 (3.3.3) packages. The expression data were normalized to transcripts per kilobase million values before further analysis. Besides, the ROC curve was used to evaluate the diagnostic value of NCAPG L using the pROC R package and ggplot2 R package.

Construction and Evaluation of the Nomogram

The RMS R package was performed to generate a nomogram. C-index and calibration curve were performed by the Hmisc R package. In our study, C-index was performed to determine the discrimination of nomogram and used a bootstrap method with 1,000 resamples to calculate C-index.

Gene Expression Profiling Interactive Analysis

Gene Expression Profiling Interactive Analysis (GEPIA) (http://gepia.cancer-pku.cn/index.html) is a user-friendly web portal for gene expression analysis based on TCGA and GTEx data (Tang et al., 2017). In the current study, GEPIA was used to analyze the expression and prognosis values of NCAPG in pan-cancer. Furthermore, we used GEPIA analysis of the correlation between NCAPG expression and pathological stage.

Kaplan–Meier Plotter Database Analysis

We used KM Plotter (http://kmplot.com), an online database that contains gene expression data and survival information of 3,452 clinical lung cancer patients, to analyze the prognostic value of NCAPG in lung cancer (Győrffy, 2021). The patient samples were separated into two groups by median expression (high expression and low expression) to analyze the overall survival (OS) with hazard ratios (HRs) with 95% confidence intervals and log-rank p values.

Immune Infiltration Analysis

TIMER (https://cistrome.shinyapps.io/timer/) web server is a comprehensive resource for systematical analysis of immune infiltrates across diverse cancer types. In this study, TIMER was used to examine the correlation between somatic copy number alterations of NCAPG and the infiltration level of B cells, CD4⁺ T cells, CD8⁺ T cells, neutrophils, macrophages, and dendritic cells.

Function and Pathway Analysis by Gene Set Enrichment Analysis

In the present research, we utilized the linkedomics database (http://www.linkedomics.org/login.php) to obtain the coexpression genes of NCAPG in LGG. The gene set KEGG v6.2. symbols.gmt,” which served as a reference gene set, was downloaded from the Molecular Signatures Database (MSigDB) (http://software.broadinstitute.org/gsea/msigdb). We used the Gene Set Enrichment Analysis (GSEA) software and clusterProfiler package to perform the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis signaling pathway of NCAPG in lung cancer (Subramanian et al., 2005; Yu et al., 2012; Vasaikar et al., 2018).

The NCAPG-Interacting Genes and Protein Analysis

The GeneMANIA database (http://www.genemania.org) was applied to construct the NCAPG interaction network (Zuberi et al., 2013). The STRING online database (https://string-db.org/) was used to construct a protein–protein interaction (PPI) network of NCAPG (Szklarczyk et al., 2019).

Prediction of Long Noncoding RNA and c8eRNA Network Construction

We used the starBase (http://starbase.sysu.edu.cn/) to predict the potential upstream microRNAs (miRNAs) of NCAPG and examine the expression, prognosis, and correlation between the miRNA-214-3p and long noncoding RNA (lncRNA); we used the starBase to predict the binding sites among the miRNA, mRNA, and lncRNA (Li et al., 2014). The lncLocator (www.csbio.sjtu.edu.cn/bioinf/lncLocator.), a subcellular localization predictor for lncRNAs based on a stacked ensemble classifier, and CPC2 (http://cpc2.cbi.pku.edu.cn), a fast and accurate coding potential calculator based on sequence intrinsic features (Kang et al., 2017; Cao et al., 2018). In this study, lncLocator and CPC2 were used to explore the subcellular localization and the protein-coding ability of thymidylate synthetase opposite strand (TYMSOS).

Cell Culture, RNA Isolation, and Real-Time Polymerase Chain Reaction

Cell culture, RNA isolation, and real-time polymerase chain reaction (PCR) assay were performed as published (Shi et al., 2019). The BEAS-2B cell line was purchased from the cell bank of Kunming Institute of Zoology and cultured in BEGM media (Lonza, CC-3170). HEK-293T was obtained from ATCC. Lung cancer cell lines, including A549, H1299, and H1975, were purchased from Cobioer (China) with STR documents; A549, H1299, and H1975 cells were all cultured in RPMI 1640 medium (Corning) supplemented with 10% fetal bovine serum (FBS) and 1% penicillin/streptomycin. HEK-293T cells were cultured in a Dulbecco modified eagle medium (Corning). The detail information of primer used in this study are as follows: NCAPG-F: GAG​GCT​GCT​GTC​GAT​TAA​GGA, NCAPG-R: AAC​TGT​CTT​ATC​ATC​CAT​CGT​GC, TYMSOS-F: ATG​ACG​CCC​GCC​TCG​GGG​GCC, TYMSOS-R: TCA​GGA​AGG​ACG​ACC​GCA​CGG​GCA​CC, FOXM1-F: CGT​CGG​CCA​CTG​ATT​CTC​AAA, FOXM1-R: GGC​AGG​GGA​TCT​CTT​AGG​TTC, miRNA-214-3p-F TTT​TTA​CTA​CTA​TGG​CGG​GTG​ATA​AAA​CGT​GTA, miRNA-214-3p-R: GCA​AGC​TGT​AAT​CGA​CGG​GAA​GAG​CAT​GCC​CAT​CC, ACTIN-F: CTTCGCGGGCGACGAT, and ACTIN-R: CCA​TAG​GAA​TCC​TTC​TGA​CC. The expression quantification was obtained with the 2^−ΔΔCt method.

Constructs, Transfection, and Infection

The FOXM1 and NCAPG shRNAs and control scrambled shRNA were transfected into HEK-293 T cells with the psPAX2/pMD2 G plasmids (Addgene) using calcium phosphate. After transfection, the cell supernatants were harvested and used to infect the A549 cell, and the stably lentiviral infected cells were selected with puromycin. Two shRNAs for knockdown of TYMSOS are as follows (sh-TYMSOS#1: 5ʹ-TAAGAAGATCTAATGCATCCT-3ʹand sh-TYMSOS#2: 5ʹ-AAC​ACT​TTA​TTA​TCA​CAT​CAG-3ʹ) were used to silence TYMSOS. Two shRNAs for knockdown of FOXM1 are as follows: sh-FOXM1 #1: 5ʹ-TCC​TGG​AGG​CTC​ACG​CCC​CCA-3ʹ and sh-FOXM1#2:5ʹ-GCCACAGGTTTCTGGCCTTGC-3ʹ, used to silence FOXM1.

Cell Proliferation, Colony Formation, and Tumorsphere Formation Assays

Cell proliferation assay was performed as previously described (Xu et al., 2020a). Indicated tumor cells were plated onto 12-well plates; the cell numbers were subsequently counted each day using an automatic cell analyzer countstar (Shanghai Ruiyu Biotech Co.). For colony formation assay, a total of 500 indicated cells/well were plated onto 6-well plates and cultured for 2 weeks at 37°C. The medium was changed every 3 days. Two weeks later, indicated cells were fixed with 4% paraformaldehyde for 30 min at room temperature and subsequently stained with 0.1% crystal violet for 30 min at room temperature.

Cell Migration Assay

The cell migration assay was performed as documented (Shi et al., 2019). For wound healing assay, to produce a wound, the monolayer cells in 6-well plate were scraped in a straight line with pipette tips. Plate was then washed with warm PBS to remove detached cells. Photographs of the scratch were taken at indicated time points using Nikon inverted microscope. For Transwell assay, 2.5 × 10⁴ cells in 100 µL serum-free medium were plated in 24-well plate chamber insert (Corning Life Sciences, catalog no. 3422), with medium containing 10% FBS at the bottom of the insert. Cells were incubated for 24 h and then fixed with 4% paraformaldehyde for 20 min. After washing, cells were stained with 0.5% crystal violet blue. The positively stained cells were examined under the light microscope.

Nuclear and Cytoplasmic Extraction

The nuclear and cytoplasmic fractions of RNA were extracted with a PARIS™ kit (Invitrogen, Thermo Fisher Scientific, Waltham, MA, USA). The cell slide was placed at the bottom of the 6-well plate, and an appropriate number of cells (3 × 10⁵/well) were cultivated. Before the experiment, the cell confluence reached 60%–70%.

Luciferase Activity Assays

The sequences of TYMSOS, NCAPG–3′ UTR, and their miR-214-3p binding sites mutant versions were synthesized and added to luciferase reporter vector pGL3 (Promega, Shanghai, China), named TYMSOS-WT, TYMSOS-Mut, NCAPG-WT, and NCAPG-Mut, respectively. The experimental steps follow the manufacturer’s protocols of Dual Luciferase Assay Kit (Promega).

Immunohistochemistry Staining

Immunohistochemistry (IHC) staining was performed as documented (Kang et al., 2017). Briefly, cancer tissues were collected; primary antibody overnight incubation and second antibody incubation were performed. Finally, develop using the instrument. The detailed information of primers used in this study is as follows: (NCAPG, catalog no. 24563-1-AP, 1:200; Proteinch, Shanghai, China). Lung cancer tissue and normal lung tissues were obtained from advanced-stage lung cancer patients from the Third Affiliated Hospital of Kunming Medical University (Yunnan Tumor Hospital), as well as ethics approval and consent to participate. This study was approved by the ethics committee of the Third Affiliated Hospital of Kunming Medical University (Yunnan Tumor Hospital), and informed consent was obtained from all patients.

Statistical Analysis

The results of Kaplan–Meier plots, PrognoScan, and GEPIA are displayed with HR and P or Cox p values from a log-rank test. The Wilcoxon signed-rank test and Wilcoxon rank-sum test were performed to investigate the expression of NCAPG in paired and nonpaired samples, respectively. The Wilcoxon signed-rank test was used to analyze the relations between NCAPG expression and clinical features. The univariate and multivariate analyses using the Cox regression model were carried out to evaluate death risk, including gender, age, number of pack-years smoked, T stage, N stage, M stage, pathological stage, and NCAPG expression. All the tests were two-sided, and p < 0.05 was considered statistically significant. p < 0.05 (*), p < 0.01 (**), and p < 0.001 (***) were considered significant.

Results

Expression Level, Clinical-Pathological Features, and Prognosis of NCAPG in Pan-Cancer

To examine the expression of NCAPG in various cancers, we applies the TIMER database to analyze the expression of NCAPG in human cancer. The results revealed that NCAPG was significantly up-regulated in 21 cancer types than normal samples, including bladder urothelial carcinoma, breast invasive carcinoma, cervical squamous cell carcinoma and endocervical adenocarcinoma, cholangiocarcinoma, colon adenocarcinoma, esophageal carcinoma, glioblastoma multiforme, head and neck squamous cell carcinoma, kidney chromophobe, kidney renal clear cell carcinoma, kidney renal papillary cell carcinoma, liver hepatocellular carcinoma, LUAD, lung squamous cell carcinoma, pheochromocytoma and paraganglioma, prostate adenocarcinoma, rectum adenocarcinoma, skin cutaneous melanoma, stomach adenocarcinoma, thyroid carcinoma, and uterine corpus endometrial carcinoma were prominently ascended when matched to normal tissues (Figure 1A). However, no significant difference of NCAPG in lymphoid neoplasm diffuse large B-cell lymphoma, brain lower-grade glioma, ovarian serous cystadenocarcinoma, sarcoma, uterine carcinosarcoma, and uveal melanoma was observed. Next, we utilized the Oncomine database to examine the differences in expression of NCAPG in various cancers. The analysis result has shown that the expression of NCAPG elevated in 17 tumors compared with the normal tissues (Figure 1B). Besides, we showed that the expression of NCAPG was decreased in acute myeloid leukemia (LAML). Furthermore, we check the level of NCAPG in TCGA pan-cancers using the GEPIA tool. The result showed that NCAPG expression in ACC, BLCA, BRCA, CESC, COAD, DLBC, esophageal carcinoma, glioblastoma multiforme, HNSC, LIHC, LUAD, LUSC, OV, pancreatic adenocarcinoma, rectum adenocarcinoma, SKCM, STAD, THYM, UCEC, and UCS was significantly elevated in cancer than corresponding normal controls. And in LAML, NCAPG was decreased (Figure 1C). Finally, we used the CCLE database to analyze the expression of NCAPG in various tumor cells; the results indicated that NCAPG was elevated in LUAD cells (Figure 1D). Taken together, these results confirm that NCAPG plays a crucial role in carcinogenesis of the human cancer.

FIGURE 1

FIGURE 1. Expression analysis for NCAPG in human cancers (A) NCAPG expression of different tumor types in pan-cancer examined by the TIMER database. (B) Increased expression of NCAPG in pan-cancer examined by the Oncomine database. (C) NCAPG expression in pan-cancer examined by TCGA and GTEx data. (D) The expression of NCAPG in different tumor cells lines was examined by the CCLE database. *p < 0.05, **p < 0.01, ***p < 0.001.

Given that NCAPG was highly expressed in diverse cancers, therefore, we next examined the expression of NCAPG whether correlated with the clinical-pathological features in NSCLC; we applied the GEPIA approach to analyze NCAPG expression and the pathological stages across various cancer types of TCGA. The results showed that high expression of NCAPG was significantly correlated with the pathological stages of human cancers, including ACC, BRCA, KIRC, KIRP, LIHC, LUAD, LUSC, OV, SKCM, TGCT, and THCA (Supplementary Figure S1).

Finally, survival time analysis of NCAPG in pan-cancer was conducted. Four prognostic indices, consisting of OS, disease-free survival (DFS), progression-free survival (PFS), and disease-specific survival (DSS) were contained. Regarding OS, up-regulation of NCAPG expression correlated with unfavorable prognosis in ACC, KIRC, KIRP, LGG, LIHC, LUAD, MESO, PAAD, PCPG, and SARC, but THYM patients elevated NCAPG expression exhibit good prognosis (Supplementary Figure S2). For DFS, about these varies cancer, up-regulation of NCAPG shown poor prognosis in ACC, KICH, KIRC, KIRP, LGG, LIHC, PAAD, PCPG, PRAD, SARC, THCA, and UVM (Supplementary Figure S3). Concerning PFS, in human cancers, up-regulation of NCAPG expression shows a bad prognosis in ACC, KICH, KIRC, KIRP, LGG, LIHC, MESO, PAAD, PCPG, PRAD, and UVM (Supplementary Figure S4A). Concerning DSS, elevated expression of NCAPG implies poor prognosis in CESC, KIRP, LIHC, MESO, SARC, and THCA (Supplementary Figure S4B). Through analysis of the OS, RFS, PFS, and DSS of NCAPG in human cancer, NCAPG may be used for an unfavorable prognostic biomarker in patients with ACC, KIRC, KIRP, LGG, LIHC, LUAD, MESO, PAAD, PCPG, and SARC. The above results demonstrated that NCAPG plays an important role in diverse cancer progression.

Correlation Between NCAPG Expression and Tumor Mutational Burden, DNA Microsatellite Instability

Studies have demonstrated that tumor mutational burden (TMB), which represents the number of somatic mutations that occur on the coding region of the tumor cell genome, may be a biomarker of sensitivity to immune checkpoint inhibitors (Campbell et al., 2017; Chalmers et al., 2017; Hu et al., 2021). Microsatellite instability (MSI) refers to the random changes in microsatellites in tumor cells. These MSI changes have been linked to increased sensitivity to checkpoint inhibitor drugs (Cortes-Ciriano et al., 2017). We figured out the TMB and MSI of every sample in human cancer and then examined the relevance between NCAPG expression and TMB, and MSI in various human cancers. The analysis result indicated that the elevated expression of NCAPG correlated positively with TMB in ACC, UCEC, TGCT, STAD, SKCM, SARC, READ, PRAD, PAAD, OV, MESO, LUSC, LUAD, LGG, KIRC, KICH, HNSC, COAD, BRCA, and BLCA but notable negative within THYM (Supplementary Figure S5A). Furthermore, up-regulated expressions of NCAPG positively correlated with MSI in COAD, ESCA, LIHC, LUSC, SARC, STAD, UCEC, and UCS but negative in DLBC, LAML, and READ (Supplementary Figure S5B).

The Expression of NCAPG in Immune and Molecular Subtypes of Human Cancers

Next, we use the TISIDB website to examine the expression of NCAPG in different immune and molecular subtypes in various tumors (Qiu et al., 2021). Immune subtype was including wound healing, IFN-gamma dominant, inflammatory, lymphocyte depleted, immunologically quit, TGF-b dominant. The results indicated that NCAPG expression was related to diverse immune subtypes in BLCA, BRCA, CESC, KIRC, KIRP, LIHC, LUAD, LUSC, MESO, OV, PRAD, SARC, SKCM, STAD, UVM, and UCEC (Supplementary Figure S6). For example, in LUAD, NCAPG had high expression in the wound healing group, IFN-gamma dominant group, and low expression in lymphocyte depleted types. For molecular subtypes of cancers, a positive connection with NCAPG expression existed in ACC, BRCA, COAD, ESCA, HNSC, KIRP, LGG, LIHC, LUSC, OV, SKCM, and UCEC (Supplementary Figure S7). Above all, these findings confirm that NCAPG has distinct patterns of expression in different cancer tissues.

NCAPG Was Highly Expressed in NSCLC

To elucidate the biological functions of NCAPG in NSCLC, we first examine the expression of NCAPG in NSCLC by using TCGA and GEO database, the results showed that the expression of NCAPG, but not of NCAPG2, was significantly up-regulated in NSCLC (Figures 2A, B), Notably, NCAPG expression was positively associated with poor clinical-pathological features, including the pathologic stage, TNM stage, primary therapy outcome, gender, and smoking status (Figures 2C, D).

FIGURE 2

FIGURE 2. NCAPG was a high expression in NSCLC. (A–D) The expression of NCAPG in NSCLC was examined by TCGA and GEO databases. (C,D) NCAPG expression is associated with clinicopathological characteristics in NSCLC. (E) Receiver operating characteristic analysis (ROC) of NCAPG in NSCLC. (*p < 0.05, **p < 0.01, ***p < 0.001).

We used logistics analysis to examine the correlation between NCAPG expression and clinicopathological characteristics. We found that NCAPG expression was positively correlated with the clinical stage including the TNM stage, pathologic stage, primary therapy outcome, gender, and smoker (Table 1). ROC curve analysis showed that the AUC for NCAPG in LUAD and LUSC patients was 0.963 and 0.995, respectively (Figure 2E). These findings indicated that the elevated expression level of NCAPG was correlated with poor pathobiology features of NSCLC patients. Given that NCAPG was significantly associated with poor pathobiology features in NSCLC patients, therefore, we further analyzed the prognosis value of NCAPG in NSCLC, the results showed that high expression of NCAPG was correlated with poor prognosis in NSCLC, this result was verified by GEO data (Figures 3A–D). The prognosis analysis for different subgroups showed that high expression of NCAPG was associated with poor prognosis (Figure 3E).

TABLE 1

TABLE 1. Logistic regression analysis of correlation between clinical–pathological characteristics and NCAPG expression in NSCLC patients.

FIGURE 3

FIGURE 3. The prognosis of NCAPG in NSCLC. (A–C) The prognosis of NCAPG in NSCLC is examined by the TCGA database. (D) The prognosis of NCAPG in NSCLC is examined by the GEO database. (E) Distinct clinical outcomes based on NCAPG expression in NSCLC patient subgroups. *p < 0.05, **p < 0.01, ***p < 0.001.

All the statistically significant prognostic factors in each multivariate Cox regression analysis were utilized to construct a prognostic nomogram, the calibration curve was drawn to examine the efficiency of the nomogram. The clinical pathologic stage, TNM stage, and NCAPG expression, were included in the nomogram to predict OS, DSS, and progression-free interval (PFI); the c-index was OS: C-index: 0.709 (0.684–0.733), DSS: C-index: 0.710 (0.679–0.741), PFI: C-index: 0.645 (0.620–0.669) (Figures 4A–F). The calibration curves all presented desirable predictions of the three nomograms for the 1-, 3-, and 5-year clinical outcomes.

FIGURE 4

FIGURE 4. Nomogram for prediction of prognosis in NSCLC patients. Construction and validation of nomograms based on NCAPG expression. Shown are the nomograms constructed to establish NCAPG expression-based risk scoring models for 1-, 3-, and 5-year overall survival (A), disease-specific survival (B), and progression-free interval (C). Calibration plots validate the efficiency of nomograms for overall survival (D), disease-specific survival (E), and progression-free interval (F). OS, overall survival; FI, progression-free interval; DSS, disease-specific survival. *p < 0.05, **p < 0.01, ***p < 0.001.

Analysis of the Function of NCAPG in LUAD

After determining the prognostic value of NCAPG in LUAD, we next explored the functions of NCAPG in NSCLC progression. The function module of linked-omics was used to examine the NCAPG coexpression gene in TCGA-NSCLC cohort (Figure 5A); the genes that were most positively related to NCAPG are exhibited in the heat map (Figures 5B). KEGG pathway enrichment analysis showed that NCAPG-coexpressed genes mainly were involved in primary immunodeficiency, hematopoietic cell lineage, T-cell receptor signaling pathway, cell adhesion molecules, cell differentiation of Th1 and Th2, cell differentiation of Th17, natural killer cell–mediated cytotoxicity, and Fc epsilon RI signaling pathway (Figure 5C). These conclusions intensively indicated that NCAPG is involved in the regulation of the immune response in NSCLC.

FIGURE 5

FIGURE 5. GO and KEGG enrichment analysis for NCAPG. (A) Volcano plot of genes differentially expressed in correlation with NCAPG. (B) Heat maps of genes positively correlated with NCAPG (top 50) (C) Analysis of the KEGG pathway of NCAPG in NSCLC by linkedomics. (D–F) GO analysis for NCAPG in NSCLC. (G) The gene interaction meshwork of NCAPG was constructed by using GeneMania. (H) Utilized the STRING to construction the protein interaction meshwork of NCAPG. *p < 0.05, **p < 0.01, ***p < 0.001.

For the GO term, NCAPG-coexpressed genes are mainly involved in the molecular function of microtubule-binding, tubulin binding, catalytic activity, acting on DNA, microtubule motor activity ATPase activity, motor activity, histone deacetylase binding, DNA helicase activity, exonuclease activity, and kinase regulator activity (Figure 5D). NCAPG-coexpressed genes are mainly involved in the biology process of mitotic nuclear division, chromosome segregation, nuclear division, microtubule cytoskeleton organization involved in mitosis, DNA replication, organelle fission, mitotic sister chromatid segregation, sister chromatid segregation, mitotic spindle organization, and spindle organization (Figure 5E). NCAPG-coexpressed genes are mainly involved in the cellular component of the chromosomal region, spindle, chromosome, centromeric region, kinetochore, microtubule, condensed chromosome kinetochore, mitotic spindle, condensed chromosome, centromeric region, kinesin complex, and intercellular bridge (Figure 5F). Furthermore, we used the Gene Mania database to analyze the gene interaction meshwork of NCAPG and the altered adjacent to genes. The results suggest the 20 constantly changing genes were closely correlated with NCAPG, including MKI87, NCAPH, NCAPD2, CDK1, HSF2, CCNB1, and AURKB (Figure 5G). Finally, we also used the STRING database to construct the PPI network of NCAPG (Figure 5H). We further analyzed the functions of this interaction gene with NCAPG; the result was suggestive of these genes mainly linked to the cell cycle and cell proliferation. Collectively, these findings emphasize the crucial role of NCAPG in regulating the progression of NSCLC.

Analysis NCAPG-Related Signaling Pathways in LUAD by GSEA

To examine the potential mechanisms influenced by NCAPG in LUAD, GSEA was conducted. The analysis result showed that NCAPG mainly participates in apoptosis, epithelial–mesenchymal transition, G2/M checkpoint, IL-2 STAT5 signal pathway, IL-6–JAK–STAT3 signal pathway, and MTORC1signaling pathway in LUAD (Figures 6A–H). Similarly, relating to the KEGG pathways, the enriched result proves that the NCAPG is mainly involved in B-cell receptor signal pathway and natural killer cell–mediated cytotoxicity signaling pathway (Figures 6I–P). The aforementioned findings intensively were indicative of NCAPG participation in lung cancer progression and immune response.

FIGURE 6

FIGURE 6. The signaling pathway for NCAPG enrichment by GSEA software. (A–P) The signaling pathway for NCAPG enrichment by GSEA software.

NCAPG Expression Is Correlated With Immune Infiltration in NSCLC

Tumor-infiltrating lymphocyte plays crucial roles in the prognosis of NSCLC (Fan et al., 2021; Zeng et al., 2021). First, we used TIMER online analysis tool to examine the relationship between the immune infiltration and the somatic copy number alterations of NCAPG in LUAD. Particularly, NCAPG CNV has memorable correlations with infiltrating levels of B cells, CD4⁺ T cells, macrophages, neutrophils, and dendritic cells (Figure 7A). Next, we further explored the correlation between NCAPG expression and the infiltration levels of different immune cells in TCGA-LUAD patients. NCAPG expression was weakly correlated with macrophages (r = 0.22, p = 2.91e-07), moderately with CD8⁺ T cells (r = 0.38, p = 3.43e-19) and B cells (r = 0.40, p = 1.1e-20), and intensively with neutrophils cells (r = 0.61, p = 3.65e-13), myeloid dendritic cells (r = 0.66, p = 1.12e-64), and CD4⁺ T cells (r = 0.67, p = 5.22e-69) (Figures 7B–G).

FIGURE 7

FIGURE 7. The correlation between NCAPG expression and immune infiltration in NSCLC. (A) The NCAPG gene copy numbers affect the infiltration level of different immune cells in LUAD. (B–G) NCAPG positively linked to the infiltration of different immune cells examined by TIMER databases. *p < 0.05, **p < 0.01, ***p < 0.001.

To further examine the functions of NCAPG in tumor immune, we used the TIMER database to explore the correlation between the NCAPG expression and immune cell–related markers in LUAD. The result indicated that NCAPG was markedly positively linked to different immune cells, containing B cells, CD8⁺ T cells, CD8⁺ T cells, M1 macrophage biomarkers, M2 macrophage biomarkers, neutrophil and dendritic cell–related markers in LUAD (Table 2).

TABLE 2

TABLE 2. The correlation between NCAPG and gene markers of immune cells examined by TIMER database.

Immune checkpoints play a significant role in cancer immune evasion (Qu et al., 2021). In consideration of the crucial roles of NCAPG in LUAD, the correlation of NCAPG with immune checkpoints related genes was assessed. NCAPG expression was significantly positively linked to CD274 (r = 0.566, p = 4.79e-45), CTLA4 (r = 0.716, p = 4.9e-82), PDCD1 (r = 0.659, p = 1.84e-65), TIGIT (r = 0.738, p = 1.06e-89), PDCD1LG2 (r = 0.596, p = 7.58e-51), HAVCR2 (r = 0.578, p = 3.21e-47), LAG3 (r = 0.551, p = 3.39e-42), and SIGLEC15 (r = 0.327, p = 2.74e-14) in LUAD (Supplementary Figures S8A–S8H). In addition, we also found significant correlations of NCAPG with immune infiltration and immune inhibitors across human heterogeneous cancers (Supplementary Figures S8I–S8J). These results indicated that NCAPG might play significant roles in the regulation of immune infiltration of LUAD.

Analysis of the Prognosis of the Immune Infiltration and NCAPG Expression in LUAD

NCAPG expression is markedly related to immune infiltration and unfavorable prognosis in LUAD (Fan et al., 2021). Next, we examined whether NCAPG expression affects the prognosis of LUAD through immune infiltrate. We used the Kaplan–Meier plotter analysis of the prognosis of NCAPG in LUAD in different related subgroups of immune cells. The analysis results suggest elevated NCAPG expression and reduced level of B cells, CD4⁺ memory T cells, CD8⁺ memory T cells, macrophages, and natural killer T cells, which had an unfavorable prognosis observed in the LUAD patient (Figures 8A–L). The prognosis of LUAD was of no relevance in other immune cell groups. These findings suggest that NCAPG may affect the survival state of LUAD patients, partially owing to immune infiltrate.

FIGURE 8

FIGURE 8. The prognosis analysis based on the expression of NCAPG in different immune cell subgroups in NSCLC. (A–L) Correlation between NCAPG expression and OS in different immune cell infiltration groups in LUAD patients. *p < 0.05, **p < 0.01, ***p < 0.001.

Depletion of NCAPG Inhibits NSCLC Progression

To examine the protein expression of NCAPG in NSCLC, we utilized an IHC assay to examine the expression of NCAPG in lung cancer. The results verify that the elevated expression of NCAPG was observed in lung cancer than the normal tissues (Figures 9A, B). Next, we examined the mRNA expression of NCAPG in BEAS2B and NSCLC cell lines. The result showed that NCAPG expressions are significantly elevated in all NSCLC cells, highest expression NCAPG observe in A549 (Figure 9C). Owing to NCAPG being highly expressed in NSCLC cells lines, we speculated that NCAPG might play a crucial role in the progression of NSCLC. We then selected the A549 cell for further studies. We inhibited the NCAPG mRNA expression using two independent lentiviral shRNAs targeting different regions of NCAPG. The efficiency and specificity of NCAPG knockdown in A549 cells were verified by quantitative reverse transcriptase (qRT)–PCR (Figure 9D). Growth curve and colony formation assays were carried out to evaluate cell proliferation, and the results showed that NCAPG silencing notably inhibited the proliferation of A549 cells (Figures 9E, F). Transwell and wound healing experiments showed that NCAPG knockdown reduced the migration of A549 cells (Figures 9G, H). Tumorsphere assay showed that NCAPG knockdown reduced the self-renewal of GSCs (Figures 9I). Collectively, these findings indicate NCAPG exerts an oncogenic role in NSCLC cells.

FIGURE 9

FIGURE 9. NCAPG silencing on cell on proliferation, migration, and self-renewal ability of NSCLC Cell. (A) IHC analysis of NCAPG in lung cancer tissues and normal tissues. (B) IHC scores statistical analysis of NCAPG in lung cancer tissues and normal tissues. (C) The expression of NCAPG in NSCLC cell lines was examined by qRT-PCR assay. (D) The establishment of NCAPG knockdown in A549 cell lines and verified by real-time qPCR (E,F) Growth curve and clone formation assays were utilized to detect silencing of NCAPG on the growth of A549 cell. (G,H) Transwell and wound healing assays were used to detect silencing of NCAPG on the migration of A549 cells. (I) A tumorsphere experiment was used to detect silencing of NCAPG on the self-renewal of A549 cells. *p < 0.05, **p < 0.01, ***p < 0.001.

Identification of miR-214-3p as the Direct Binding miRNA of NCAPG in LUAD

To uncover whether NCAPG was regulated by miRNAs, we predicted the potential miRNAs that bind with NCAPG by using the RNA22 (Loher and Rigoutsos, 2012), miRmap (Goossens et al., 2019), DIANA-microT web server v5.0 (Paraskevopoulou et al., 2013), miRanda (Rost and McGregor, 2012), and TargetScan (Agarwal et al., 2015). Eventually, only four miRNAs, including miRNA-15a, miRNA-147a, miRNA-375, and miRNA-214-3p, were commonly predicted. Next, we assessed their expression levels and prognostic values in NSCLC using TCGA-LUAD data (Figures 10A, B). According to the ceRNA theory, the miRNA could bind with 3′ UTR of mRNA and inhibit the expression of mRNA; thus, the expression of miRNA and mRNA should have a negative correlation pattern. We performed the expression correlation between NCAPG and miRNA-214-3p, the result demonstrated that NCAPG was significantly negatively correlated with miRNA-214-3p in LUAD (Figure 10C). In addition, Among the four miRNAs, only miR-214-3p was significantly decreased in LUAD than a control group. Survival analysis demonstrated that the elevated miR-214-3p expression has a favorable prognosis in LUAD. By integration of the expression analysis, correlation analysis, and survival analysis, we supposed that miR-214-3p was the most potential miRNA that may be potentially bound to NCAPG, the potential binding site between the miRNA-214-3p and NCAPG 3′ UTR were predicted base-pairing with miRNA-214-3p by starBase (Figure 10D). Next, we analyzed the expressions of miRNA-214-3p in whole blood of human lung cancer and lung cancer by using the GEO database (Patnaik et al., 2017); these data also revealed miRNA-214-3p expressions are down-regulated in the whole blood of human lung cancer and lung cancer (Figures 10E,F). Finally, we also analyzed the expressions of miRNA-214-3p in normal human bronchial epithelial cells and NSCLC cell lines. Consistent with our findings in TCGA, we found that miRNA-214-3p expressions are down-regulated in all NSCLC cell lines (Figure 10G). To investigate the effects of miRNA-214-3p on the expression of NCAPG, we conducted overexpression analysis. Overexpressing miRNA-214-3p significantly reduced the mRNA expressions of NCAPG in A549 cells (Figure 10H). The luciferase assay showed that transfection of miRNA-214-3p mimics significantly the reduced relative luciferase activity of NCAPG-3UTR-WT–treated lung cancer cells, but did not affect that of NCAPG-3UTR-MUT–treated lung cancer cells (Figure 10I). Collectively, these data imply that miRNA-214-3p might target the NCAPG and plays pivotal roles in the regulation of NCAPG expression in NSCLC.

FIGURE 10

FIGURE 10. Identification of upstream miRNAs of NCAPG in LUAD. (A) The expression of miRNA-15a, miRNA-147a, miRNA-375, and miRNA-214-3p in TCGA-LUAD. (B) The prognostic values of miRNA-15a, miRNA-147a, miRNA-375, and miRNA-214-3p in TCGA-LUAD. (C) The correlation between expression of miRNA-214-3p and NCAPG in LUAD. (D) The potential binding sites of miR-214-3p and NCAPG. (E) The expression of miRNA-214-3p in NSCLC was examined by the GEO database. (F) The expression of miRNA-214-3p in NSCLC determined by GEO database (G) The expression of miRNA-214-3p in NSCLC cell lines. (H) The expression of NCAPG in A549 cell after overexpression of miRNA-214-3p examined by the qRT-PCR assay. (I) Relative luciferase activities of wild type (WT) and mutated (MUT) NCAPG–3′ UTR reporter plasmid in A549 cells cotransfected with miRNA-214-3p mimics. *p < 0.05, **p < 0.01, ***p < 0.001.

Identification of TYMSOS as Upstream lncRNA of miRNA-214-3p in LUAD

Diverse regulatory mechanisms of gene expression have been well established in the NSCLC progression. Emerging evidence has demonstrated that lncRNAs play an indispensable role in regulating gene expression. Next, the upstream lncRNAs of miRNA-214-3p were predicted using the starBase database and DIANA-LncBase v3 (Karagkouni et al., 2020). There are four lncRNAs were predicted, including LINC01228, RUSC1-AS1, SNHG3, and TYMSOS. Identically, the expression levels of and prognostic values of lncRNAs in LUAD were determined by using TCGA-LUAD data. The analysis results confirmed that LINC01228, RUSC1-AS1, and SNHG3 were down-regulated in NSCLC, on the contrary, TYMSOS was significantly up-regulation in NSCLC and high expression of TYMSOS was correlated with unfavorable prognosis in NSCLC (Figures 11A,B). In the light of the competitive endogenous RNAs theory, lncRNAs able to up-regulation the mRNA level via reducing the expression of miRNAs. Therefore, the elevated expression of lncRNA should decrease the level of miRNA, and positive correlation between mRNA. Among all the four lncRNAs, only TYMSOS was negatively correlated with miRNA-214-3p in LUAD (Figure 11C), and TYMSOS was positively correlated with NCAPG in LUAD (Figure 11D). The target sites in the miRNA-214-3p and TYMSOS were predicted to pair with miRNA-214-3p by TarBase (Figure 11E). Moreover, the TYMSOS was markedly elevated in LUAD matched to normal controls by analysis of the GEO databases (Mezheyeuski et al., 2018) (Figure 11F). Finally, to further confirm the interaction between TYMSOS and miRNA-214-3p, the luciferase assay was performed, the luciferase assay indicated that transfection of miRNA-214-3p mimics significantly decreased the relative luciferase activity of TYMSOS-WT–treated lung cancer cells, but did not affect that of TYMSOS-MUT–treated lung cancer cells (Figure 11G). Taken together, the data suggest that TYMSOS acted as a sponge for miRNA-214-3p in NSCLC cells.

FIGURE 11

FIGURE 11. Identification of upstream lncRNAs of miRNA-214-3p in LUAD. (A) The expression of LINC01228, RUSC1-AS1, SNHG3, and TYMSOS in TCGA-LUAD. (B) The prognostic values of LINC01228, RUSC1-AS1, SNHG3, and TYMSOS in TCGA-LUAD. (C) The correlation between expression of miRNA-214-3p and TYMSOS in LUAD. (D) The correlation between expression of NCAPG and TYMSOS in LUAD. (E) The potential binding sites of miR-214-3p and TYMSOS. (F) The expression of TYMSOS in GSE81089 and control normal samples determined by GEO database. (G) Relative luciferase activities of wild-type (WT) and mutated (MUT) THTMSOS reporter plasmid in A549 cells cotransfected with miRNA-214-3p mimics. (H) The expression of FOXM1 in LUAD was examined by using the GEPIA database. (I) The pathological stage analysis for FOXM1 in LUAD was determined by the GEPIA database. (J) The overall survival analysis for FOXM1 in LUAD was examined by using the GEPIA database. (K) The correlation between expression of FOXM1 and TYMSOS in TCGA-LUAD. (L) The correlation between expression of FOXM1 and NCAPG in TCGA-LUAD. (M) TYMSOS expression in A549 cell after knockdown FOXM1 was examined by real-time PCR. *p < 0.05, **p < 0.01, ***p < 0.001.

To explore the up-regulation mechanism of TYMSOS in LUAD by using the JASPAR, PROMO, ConTra v3, and UCSC databases. We identified transcription factors that may regulate the expression of TYMSOS. Data indicated that FOXM1 was the most qualified one. We examined the mRNA expression of FOXM1 by using the GEPIA database. FOXM1expression was substantially elevated in LUAD (Figure 11H), the expression of FOXM1 was significantly elevated in LUAD patients than normal controls based on tumor grade (Figure 11I). Up-regulation of FOXM1expression is associated with the poor prognosis in lung cancer patients based on the GEPIA database (Figure 11J). We examined the expression relationship between FOXM1 and TYMSOS in TCGA-LUAD. The expression of FOXM1 was markedly linked to the expression of TYMSOS and NCAPG in LUAD (Figures 11K, L). Next, the knockdown plasmids of FOXM1 were constructed and respectively transfected into the A549 cell. Then, qPCR assays revealed that knockdown of FOXM1 can markedly reduce the expression of TYMSOS levels (Figure 11M). Collectively, these findings indicated that FOXM1 may be a potential transcriptional activator of TYMSOS.

Depletion of TYMSOS Inhibits Cell Proliferation, Migration, and Self-Renewal Abilities of NSCLC Cells

The underlying mechanisms of TYMSOS were determined by cellular localization. The subcellular localization of TYMSOS used the lncLocator and by performing a nuclear/cytoplasmic separate experiment. The results indicated that TYMSOS was mainly located in the cytoplasm of A549 cells (Figures 12A,B). We also analyzed the coding potential of TYMSOS by performing the coding potential calculator, the results showed that TYMSOS does not possess the protein-coding ability (Figure 12C). We first examined the expression of TYMSOS in NSCLC cells lines; the analysis results indicated that TYMSOS was significantly up-regulated in NSCLC cells (Figure 12D). To explore the biological function of TYMSOS in lung cancer cells, the TYMSOS knockdown plasmid was transfected into the A549 cell. The results indicated that the expression of NCAPG and miRNA-214-3p were decreased and increased, respectively, after the knockdown of TYMSOS in A549 cells (Figure 12E). The growth curve and colony formation experiment showed that TYMSOS knockdown significantly inhibits cell proliferation (Figures 12F,G). To estimate the significance of TYMSOS in NSCLC progression, migration and self-renewal of A549 cells were investigated after the knockdown of TYMSOS. The tumorsphere assays indicated that the self-renewal ability of A549 cells was significantly suppressed by down-regulation of TYMSOS (Figure 12H). The Transwell assays indicated that the migratory ability of A549 cells was significantly suppressed by the down-regulation of TYMSOS (Figure 12I). These findings indicated that TYMSOS functioned as an oncogenic lncRNA in promoting cell growth, metastasis, and self-renewal of NSCLC cells. These data suggest that the TYMSOS/miRNA-214-3p/NCAPG axis plays a significant role in the progression of lung cancer.

FIGURE 12

FIGURE 12. TYMSOS silencing on cell on proliferation, migration, and self-renewal ability of NSCLC cell. (A) The cellular localization for TYMSOS was predicted by using lncLocator. (B) The cellular localization for TYMSOS was examined by nuclear and cytoplasmic extraction experiments. (C) The coding ability for TYMSOS was predicted by using Coding Potential Calculator. (D) The expression of TYMSOS in NSCLC cells lines was examined by qRT-PCR assay. (E) The expression of miRNA-214-3p and NCAPG in A549 cells after knockdown of TYMSOS was examined by qRT-PCR assay. (F,G) The growth curve and clone formation assays were used to detect the silencing of TYMSOS on the growth of A549 cells. (H) The tumorsphere assay was utilized to detect silencing of TYMSOS on the self-renewal of the A549 cell. (I) The Transwell assay was utilized to detect silencing of TYMSOS on the migration of A549 cells. *p < 0.05, **p < 0.01, ***p < 0.001.

Discussion

According to the global cancer statistics in 2020, lung cancer results in the highest morbidity rates and brought huge economic burden worldwide (Sung et al., 2021). Despite there being many different treatments, including immunotherapy, more than half of the patients are diagnosed at an advanced stage, and the survival rate of patients is very low (Hirsch et al., 2017). The initiation and progression of lung cancer are an exceedingly complicated process that involves genetic mutations, tumor microenvironment, and the dysregulation of epigenetic pathways (Govindan et al., 2012; Duruisseaux and Esteller, 2018; Shi et al., 2019). Concurrently, epigenetic changes in lung cancer such as histone modifications (Chen et al., 2018), DNA methylation (Hulbert et al., 2017) and noncoding RNAs (ncRNAs) (Li et al., 2019) have been extensively studied. The RNAs mainly include snRNAs, snoRNAs, piRNAs, tRNAs, and ribosomal RNAs, miRNAs, lncRNAs, circular RNAs (Esteller, 2011; Hombach and Kretz, 2016; Li et al., 2018). Accumulating evidence has revealed that numerous ncRNAs are dysregulated and involved in lung cancer initiation and progression (Chen et al., 2016; Wang et al., 2018; Iqbal et al., 2019). Many studies demonstrate that miRNAs and lncRNAs are involved in the pathogenesis of lung cancer. For example, lncRNA H19 the elevated cell growth and metastasis ability by reducing the level of miR-200 (Zhao et al., 2019). LncRNA TUC338 promotes invasion of lung cancer by activating the MAPK pathway (Zhang et al., 2018).

NCAPG plays key roles in mitotic chromosome condensation and cancer progression (Sun et al., 2020). It has also been reported that NCAPG serves as a prognostic biomarker in HCC (Xiao et al., 2020). Combining bioinformatics analyses and biological function validation, our study provides evidence that NCAPG was elevated in lung cancer, and its high expression was correlated with poor clinical-pathological features, including the pathologic stage, TNM stage, primary therapy outcome, gender, and smoking status. Finally, OS analyses showed that the up-regulated expression of NCAPG was associated with a poorer survival state in NSCLC patients. Aberrant expression of NCAPG has been elucidated in various human cancers. For example, NCAPG exerted oncogenic functions in hepatocellular carcinoma via regulating PI3K/AKT signaling (Gong et al., 2019). Furthermore, NCAPG has been reported to inhibit cardia adenocarcinoma apoptosis and promote epithelial–mesenchymal transition through the Wnt/β-catenin signaling pathway (Zhang et al., 2021). All these studies indicated that NCAPG might exert a pivotal function in human cancers. In our study, we observed increased NCAPG in NSCLC, which was also associated with a poor prognosis. Similarly, high expression of NCAPG is associated with poor prognosis in ovarian cancer (Xu et al., 2020b).

For the function of NCAPG in NSCLC, we performed the GSEA and KEGG pathway enrichment analyses in TCGA-LUAD; the results indicated that NCAPG participates in the immune system in LUAD. It has been well recognized that the tumor microenvironment is able to affect the productiveness of immunotherapy and involved in the process of cancer (Fridman et al., 2017; Jarosz-Biej et al., 2019; McLaughlin et al., 2020). It has been demonstrated that immune modulators play crucial roles in tumor cell evasion of immune surveillance (Xue et al., 2017). Studies have confirmed that the increased immune checkpoints such as PD-L1, CTLA-4, TIM-3, and LAG3 in glioma help tumor immune evasion, leading to T-cell dysfunction (Xue et al., 2017; Woroniecka et al., 2018). Here, we take the lead to review that up-regulation of NCAPG was positively correlated with the immune cell infiltration in lung cancer. Moreover, we also found that the NCAPG expression was not only significantly positively associated with various immune cell markers but also positively correlated with the expression of immune checkpoints–related genes. More significantly, we showed that NCAPG expression could affect the OS state of LUAD patients through immune infiltration. These results suggest that NCAPG could as an effective immunotherapy target in LUAD. However, the crucial functions of NCAPG in the tumor-immune regulation still need further research.

A growing number of studies have reported that lncRNA plays an important role in gene expression regulation (Zhang et al., 2020). To figure out the upstream regulatory microRNA about NCAPG, we used different databases to predict potential microRNAs that may be linked to NCAPG. After performing comprehensive analysis and functional experiments, miRNA-214-3p was chosen as miRNA that may be binding to the NCAPG. Furthermore, previous findings as well indicated that miRNA-214-3p played inhibitory roles in tumor cell proliferation (Li et al., 2020). These results indicate that miRNA-214-3p could act as an upstream regulatory miRNA of NCAPG in LUAD. Through conducting lncRNA molecular characteristics, expression, survival, and correlation analysis, one of the most potential elevated lncRNA, TYMSOS, was identified. It has been reported to function as oncogenes in gastric cancer (Gu et al., 2021). For instance, TYMSOS through binding the miR-4739 and resulting in the elevated expression of ZNF703 gives rise to the progression of gastric cancer (Gu et al., 2021). In this study, we found that the TYMSOS and NCAPG expression was markedly elevated in LUAD matched to the control group, and the elevated expression of TYMSOS and NCAPG was associated with the poor outcome by survival analysis. On the contrary, miR-214-3p expression was reduced in LUAD, and the elevated miR-214-3p expression indicated better OS prognosis in LUAD. In addition, we also identified FOXM1 as a potential transcriptional activator of TYMSOS. We also performed functional assay to verify the function of TYMSOS and NCAPG in NSCLC; the results verified that knockdown of TYMSOS and NCAPG significantly inhibited the cell proliferation, migration, and self-renewal abilities of NSCLC cells. We also confirm the interaction between TYMSOS, miRNA-214-3p, and NCAPG. Although various molecular mechanisms can lead to increased gene expression, such as the DNA methylation, regulation of transcription factors, and gene copy number amplification, lncRNA/miRNA axis is one of the main regulatory mechanisms of gene expression.

We revealed that NCAPG was up-regulated in different cancers and positively correlated to poor prognosis in NSCLC. The upstream complex molecular regulation mechanism of NCAPG was first revealed by us, that is, FOXM1/TYMSOS/miRNA-214-3p (Figure 13). Furthermore, our work demonstrates that NCAPG was further implicated in the alteration of the tumor microenvironment. TYMSOS and NCAPG knockdown inhibited cell proliferation, cell migration, and cancer stem cell self-renewal in A549 cells. These findings suggest that NCAPG may play an oncogene role in the progression of lung cancer and represent an effective target of immunotherapy in the NSCLC.

FIGURE 13

FIGURE 13. LncRNA-TYMSOS/miRNA-214-3p/NCAPG axis facilitated NSCLC progression.

To the best of our knowledge, this is the first study to explore the correlation between NCAPG and NSCLC. However, there are some limitations to our research. First, our study was based on expression data extracted from TCGA but may be more convincing if supported by a prospective clinical study. Furthermore, the biological functions of NCAPG need to be further explored in in vivo experiments. In the future, we will pay more attention to the function of NCAPG in tumor progression and tumor microenvironment regulation of NSCLC. Furthermore, we will perform more in vivo and in vitro experiments to explore the function and the potential molecular mechanisms of NCAPG in tumor progression and tumor microenvironment regulation of NSCLC.

Conclusion

Our study found that NCAPG expression was increased in NSCLC, which was also associated with a poor prognosis. Furthermore, NCAPG might be involved in the progression of NSCLC by regulating the function of immune-infiltrating cells and immune response–related signaling pathways. Herein, we revealed the biological functions of NCAPG in NSCLC and offered a potential strategy for the diagnosis and treatment of NSCLC.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.

Author Contributions

YY, XJ, and LT designed this work and performed related function assay. JW, and DZ analyzed data. WC and LD supervised and wrote the manuscript. All authors have read and approved the final version of the manuscript.

Funding

This work was supported by National Nature Science Foundation of China (No. 82160508) and Yunnan Applied Basic Research Projects (YNWRMY-2019-067, 2019FE001) and Yunnan Province Specialized Training Grant for High-Level Healthcare Professionals (D-201614).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The handling editor IH is currently organizing a Research Topic with one of the authors WC.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

We would like to thank the Core Technology Facility, Kunming Institute of Zoology, Chinese Academy of Sciences for providing us with bioinformatics analysis. And we are grateful to Guolan Ma for her technical support.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2021.785767/full#supplementary-material

Abbreviations

ACC, adrenocortical carcinoma; BLCA, bladder urothelial carcinoma; BRCA, breast invasive carcinoma; CESC, cervical squamous cell carcinoma and endocervical adenocarcinoma; COAD, colon adenocarcinoma; CSC, cancer stem cell; ESCA, esophageal carcinoma; GBM, glioblastoma multiforme; HNSC, head and neck squamous cell carcinoma; KICH, kidney chromophobe; KIRC, kidney renal clear cell carcinoma; KIRP, kidney renal papillary cell carcinoma; LAML, acute myeloid leukemia; LGG, brain lower grade glioma; LIHC, liver hepatocellular carcinoma; LUAD, lung adenocarcinoma; LUSC, lung squamous cell carcinoma; MESO, mesothelioma; MSI, microsatellite instability; OV, ovarian serous cystadenocarcinoma; PAAD, pancreatic adenocarcinoma; PCPG, pheochromocytoma and paraganglioma; PRAD, prostate adenocarcinoma; READ, colon adenocarcinoma; READ, rectum adenocarcinoma; STAD, stomach adenocarcinoma; TGCT, testicular germ cell tumors; THCA, thyroid carcinoma; TMB, tumor mutational burden; UCEC, uterine corpus endometrial carcinoma; UCS, uterine carcinosarcoma.

References

Agarwal, V., Bell, G. W., Nam, J. W., and Bartel, D. P. (2015). Predicting Effective microRNA Target Sites in Mammalian mRNAs. Elife 4, e05005. doi:10.7554/elife.05005

CrossRef Full Text | Google Scholar

Arai, T., Okato, A., Yamada, Y., Sugawara, S., Kurozumi, A., Kojima, S., et al. (2018). Regulation ofNCAPGbymiR-99a-3p(passenger Strand) Inhibits Cancer Cell Aggressiveness and Is Involved in CRPC. Cancer Med. 7 (5), 1988–2002. doi:10.1002/cam4.1455

PubMed Abstract | CrossRef Full Text | Google Scholar

Campbell, B. B., Light, N., Fabrizio, D., Zatzman, M., Fuligni, F., de Borja, R., et al. (2017). Comprehensive Analysis of Hypermutation in Human Cancer. Cell 171 (5), 1042–1056. e10. doi:10.1016/j.cell.2017.09.048

PubMed Abstract | CrossRef Full Text | Google Scholar

Cao, Z., Pan, X., Yang, Y., Huang, Y., and Shen, H.-B. (2018). The lncLocator: a Subcellular Localization Predictor for Long Non-coding RNAs Based on a Stacked Ensemble Classifier. Bioinformatics 34 (13), 2185–2194. doi:10.1093/bioinformatics/bty085

PubMed Abstract | CrossRef Full Text | Google Scholar

Chalmers, Z. R., Connelly, C. F., Fabrizio, D., Gay, L., Ali, S. M., Ennis, R., et al. (2017). Analysis of 100,000 Human Cancer Genomes Reveals the Landscape of Tumor Mutational burden. Genome Med. 9 (1), 34. doi:10.1186/s13073-017-0424-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Y., Liu, X., Li, Y., Quan, C., Zheng, L., and Huang, K. (2018). Lung Cancer Therapy Targeting Histone Methylation: Opportunities and Challenges. Comput. Struct. Biotechnol. J. 16, 211–223. doi:10.1016/j.csbj.2018.06.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Z., Li, J.-L., Lin, S., Cao, C., Gimbrone, N. T., Yang, R., et al. (2016). cAMP/CREB-Regulated LINC00473 marks LKB1-Inactivated Lung Cancer and Mediates Tumor Growth. J. Clin. Invest. 126 (6), 2267–2279. doi:10.1172/jci85250

PubMed Abstract | CrossRef Full Text | Google Scholar

Cortes-Ciriano, I., Lee, S., Park, W.-Y., Kim, T.-M., and Park, P. J. (2017). A Molecular Portrait of Microsatellite Instability across Multiple Cancers. Nat. Commun. 8, 15180. doi:10.1038/ncomms15180

PubMed Abstract | CrossRef Full Text | Google Scholar

Duruisseaux, M., and Esteller, M. (2018). Lung Cancer Epigenetics: From Knowledge to Applications. Semin. Cancer Biol. 51, 116–128. doi:10.1016/j.semcancer.2017.09.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Eberlein, A., Takasuga, A., Setoguchi, K., Pfuhl, R., Flisikowski, K., Fries, R., et al. (2009). Dissection of Genetic Factors Modulating Fetal Growth in Cattle Indicates a Substantial Role of the Non-SMC Condensin I Complex, Subunit G (NCAPG) Gene. Genetics 183 (3), 951–964. doi:10.1534/genetics.109.106476

PubMed Abstract | CrossRef Full Text | Google Scholar

Esteller, M. (2011). Non-coding RNAs in Human Disease. Nat. Rev. Genet. 12 (12), 861–874. doi:10.1038/nrg3074

PubMed Abstract | CrossRef Full Text | Google Scholar

Fan, Y., Liu, B., Chen, F., Song, Z., Han, B., Meng, Y., et al. (2021). Hepcidin Upregulation in Lung Cancer: A Potential Therapeutic Target Associated with Immune Infiltration. Front. Immunol. 12, 612144. doi:10.3389/fimmu.2021.612144

PubMed Abstract | CrossRef Full Text | Google Scholar

Fridman, W. H., Zitvogel, L., Sautès–Fridman, C., and Kroemer, G. (2017). The Immune Contexture in Cancer Prognosis and Treatment. Nat. Rev. Clin. Oncol. 14 (12), 717–734. doi:10.1038/nrclinonc.2017.101

PubMed Abstract | CrossRef Full Text | Google Scholar

Gong, C., Ai, J., Fan, Y., Gao, J., Liu, W., Feng, Q., et al. (2019). NCAPG Promotes the Proliferation of Hepatocellular Carcinoma through PI3K/AKT Signaling. Onco Targets Ther. 12, 8537–8552. doi:10.2147/ott.s217916

PubMed Abstract | CrossRef Full Text | Google Scholar

Goossens, E. A. C., de Vries, M. R., Simons, K. H., Putter, H., Quax, P. H. A., and Nossent, A. Y. (2019). miRMap: Profiling 14q32 microRNA Expression and DNA Methylation throughout the Human Vasculature. Front. Cardiovasc. Med. 6, 113. doi:10.3389/fcvm.2019.00113

PubMed Abstract | CrossRef Full Text | Google Scholar

Govindan, R., Ding, L., Griffith, M., Subramanian, J., Dees, N. D., Kanchi, K. L., et al. (2012). Genomic Landscape of Non-small Cell Lung Cancer in Smokers and Never-Smokers. Cell 150 (6), 1121–1134. doi:10.1016/j.cell.2012.08.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Gu, Y., Wan, C., Zhou, G., Zhu, J., Shi, Z., and Zhuang, Z. (2021). TYMSOS Drives the Proliferation, Migration, and Invasion of Gastric Cancer Cells by Regulating ZNF703 via Sponging miR-4739. Cell Biol Int 45 (8), 1710–1719. doi:10.1002/cbin.11610

PubMed Abstract | CrossRef Full Text | Google Scholar

Győrffy, B. (2021). Survival Analysis across the Entire Transcriptome Identifies Biomarkers with the Highest Prognostic Power in Breast Cancer. Comput. Struct. Biotechnol. J. 19, 4101–4109. doi:10.1016/j.csbj.2021.07.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Herbst, R. S., Morgensztern, D., and Boshoff, C. (2018). The Biology and Management of Non-small Cell Lung Cancer. Nature 553 (7689), 446–454. doi:10.1038/nature25183

PubMed Abstract | CrossRef Full Text | Google Scholar

Hirsch, F. R., Scagliotti, G. V., Mulshine, J. L., Kwon, R., Curran, W. J., Wu, Y.-L., et al. (2017). Lung Cancer: Current Therapies and New Targeted Treatments. The Lancet 389 (10066), 299–311. doi:10.1016/s0140-6736(16)30958-8

CrossRef Full Text | Google Scholar

Hombach, S., and Kretz, M. (2016). Non-coding RNAs: Classification, Biology and Functioning. Adv. Exp. Med. Biol. 937, 3–17. doi:10.1007/978-3-319-42059-2_1

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, J., Qiu, D., Yu, A., Hu, J., Deng, H., Li, H., et al. (2021). YTHDF1 Is a Potential Pan-Cancer Biomarker for Prognosis and Immunotherapy. Front. Oncol. 11, 607224. doi:10.3389/fonc.2021.607224

PubMed Abstract | CrossRef Full Text | Google Scholar

Hulbert, A., Jusue-Torres, I., Stark, A., Chen, C., Rodgers, K., Lee, B., et al. (2017). Early Detection of Lung Cancer Using DNA Promoter Hypermethylation in Plasma and Sputum. Clin. Cancer Res. 23 (8), 1998–2005. doi:10.1158/1078-0432.ccr-16-1371

PubMed Abstract | CrossRef Full Text | Google Scholar

Iqbal, M. A., Arora, S., Prakasam, G., Calin, G. A., and Syed, M. A. (2019). MicroRNA in Lung Cancer: Role, Mechanisms, Pathways and Therapeutic Relevance. Mol. Aspects Med. 70, 3–20. doi:10.1016/j.mam.2018.07.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Jarosz-Biej, M., Smolarczyk, R., Cichoń, T., and Kułach, N. (2019). Tumor Microenvironment as A "Game Changer" in Cancer Radiotherapy. Int. J. Mol. Sci. 20 (13), 3212. doi:10.3390/ijms20133212

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, L., Ren, L., Chen, H., Pan, J., Zhang, Z., Kuang, X., et al. (2020). NCAPG Confers Trastuzumab Resistance via Activating SRC/STAT3 Signaling Pathway in HER2-Positive Breast Cancer. Cell Death Dis 11 (7), 547. doi:10.1038/s41419-020-02753-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Kang, Y.-J., Yang, D.-C., Kong, L., Hou, M., Meng, Y.-Q., Wei, L., et al. (2017). CPC2: a Fast and Accurate Coding Potential Calculator Based on Sequence Intrinsic Features. Nucleic Acids Res. 45 (W1), W12–w16. doi:10.1093/nar/gkx428

PubMed Abstract | CrossRef Full Text | Google Scholar

Karagkouni, D., Paraskevopoulou, M. D., Tastsoglou, S., Skoufos, G., Karavangeli, A., Pierros, V., et al. (2020). DIANA-LncBase V3: Indexing Experimentally Supported miRNA Targets on Non-coding Transcripts. Nucleic Acids Res. 48 (D1), D101–d110. doi:10.1093/nar/gkz1036

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, J. H., Liu, S., Zhou, H., Qu, L. H., and Yang, J. H. (2014). starBase v2.0: Decoding miRNA-ceRNA, miRNA-ncRNA and Protein-RNA Interaction Networks from Large-Scale CLIP-Seq Data. Nucleic Acids Res. 42, D92–D97. Database issue. doi:10.1093/nar/gkt1248

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, S., Sun, Y., Tian, T., Qin, X., Lin, S., Zhang, T., et al. (2020). MicroRNA-214-3p Modified Tetrahedral Framework Nucleic Acids Target Survivin to Induce Tumour Cell Apoptosis. Cell Prolif 53 (1), e12708. doi:10.1111/cpr.12708

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, X., Yang, L., and Chen, L.-L. (2018). The Biogenesis, Functions, and Challenges of Circular RNAs. Mol. Cel 71 (3), 428–442. doi:10.1016/j.molcel.2018.06.034

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Y., Yin, Z., Fan, J., Zhang, S., and Yang, W. (2019). The Roles of Exosomal miRNAs and lncRNAs in Lung Diseases. Sig Transduct Target. Ther. 4, 47. doi:10.1038/s41392-019-0080-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Liang, R., Li, X., Li, W., Zhu, X., and Li, C. (2021). DNA Methylation in Lung Cancer Patients: Opening a "window of Life" under Precision Medicine. Biomed. Pharmacother. 144, 112202. doi:10.1016/j.biopha.2021.112202

PubMed Abstract | CrossRef Full Text | Google Scholar

Loher, P., and Rigoutsos, I. (2012). Interactive Exploration of RNA22 microRNA Target Predictions. Bioinformatics 28 (24), 3322–3323. doi:10.1093/bioinformatics/bts615

PubMed Abstract | CrossRef Full Text | Google Scholar

McLaughlin, M., Patin, E. C., Pedersen, M., Wilkins, A., Dillon, M. T., Melcher, A. A., et al. (2020). Inflammatory Microenvironment Remodelling by Tumour Cells after Radiotherapy. Nat. Rev. Cancer 20 (4), 203–217. doi:10.1038/s41568-020-0246-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Mezheyeuski, A., Bergsland, C. H., Backman, M., Djureinovic, D., Sjöblom, T., Bruun, J., et al. (2018). Multispectral Imaging for Quantitative and Compartment-specific Immune Infiltrates Reveals Distinct Immune Profiles that Classify Lung Cancer Patients. J. Pathol. 244 (4), 421–431. doi:10.1002/path.5026

CrossRef Full Text | Google Scholar

Molina, J. R., Yang, P., Cassivi, S. D., Schild, S. E., and Adjei, A. A. (2008). Non-small Cell Lung Cancer: Epidemiology, Risk Factors, Treatment, and Survivorship. Mayo Clinic Proc. 83 (5), 584–594. doi:10.4065/83.5.584

CrossRef Full Text | Google Scholar

Morgensztern, D., Ng, S. H., Gao, F., and Govindan, R. (2010). Trends in Stage Distribution for Patients with Non-small Cell Lung Cancer: a National Cancer Database Survey. J. Thorac. Oncol. : official Publ. Int. Assoc. Study Lung Cancer 5 (1), 29–33. doi:10.1097/JTO.0b013e3181c5920c

CrossRef Full Text | Google Scholar

Paraskevopoulou, M. D., Georgakilas, G., Kostoulas, N., Vlachos, I. S., Vergoulis, T., Reczko, M., et al. (2013). DIANA-microT Web Server v5.0: Service Integration into miRNA Functional Analysis Workflows. Nucleic Acids Res. 41, W169–W173. Web Server issue. doi:10.1093/nar/gkt393

PubMed Abstract | CrossRef Full Text | Google Scholar

Patnaik, S. K., Kannisto, E. D., Mallick, R., Vachani, A., and Yendamuri, S. (2017). Whole Blood microRNA Expression May Not Be Useful for Screening Non-small Cell Lung Cancer. PLoS One 12 (7), e0181926. doi:10.1371/journal.pone.0181926

PubMed Abstract | CrossRef Full Text | Google Scholar

Qiu, M.-j., Wang, Q.-s., Li, Q.-t., Zhu, L.-s., Li, Y.-n., Yang, S.-l., et al. (2021). KIF18B Is a Prognostic Biomarker and Correlates with Immune Infiltrates in Pan-Cancer. Front. Mol. Biosci. 8, 559800. doi:10.3389/fmolb.2021.559800

PubMed Abstract | CrossRef Full Text | Google Scholar

Qu, S., Liu, J., and Wang, H. (2021). EVA1B to Evaluate the Tumor Immune Microenvironment and Clinical Prognosis in Glioma. Front. Immunol. 12, 648416. doi:10.3389/fimmu.2021.648416

PubMed Abstract | CrossRef Full Text | Google Scholar

Rost, G. C., and McGregor, K. K. (2012). Miranda Rights Comprehension in Young Adults with Specific Language Impairment. Am. J. Speech Lang. Pathol. 21 (2), 101–108. doi:10.1044/1058-0360(2011/10-0094)

CrossRef Full Text | Google Scholar

Shi, Y., Fan, S., Wu, M., Zuo, Z., Li, X., Jiang, L., et al. (2019). YTHDF1 Links Hypoxia Adaptation and Non-small Cell Lung Cancer Progression. Nat. Commun. 10 (1), 4892. doi:10.1038/s41467-019-12801-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Siegel, R. L., Miller, K. D., and Jemal, A. (2020). Cancer Statistics, 2020. CA A. Cancer J. Clin. 70 (1), 7–30. doi:10.3322/caac.21590

CrossRef Full Text | Google Scholar

Song, B., Du, J., Song, D.-f., Ren, J.-c., and Feng, Y. (2018). Dysregulation of NCAPG, KNL1, miR-148a-3p, miR-193b-3p, and miR-1179 May Contribute to the Progression of Gastric Cancer. Biol. Res. 51 (1), 44. doi:10.1186/s40659-018-0192-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Subramanian, A., Tamayo, P., Mootha, V. K., Mukherjee, S., Ebert, B. L., Gillette, M. A., et al. (2005). Gene Set Enrichment Analysis: a Knowledge-Based Approach for Interpreting Genome-wide Expression Profiles. Proc. Natl. Acad. Sci. 102 (43), 15545–15550. doi:10.1073/pnas.0506580102

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, D.-P., Lin, C.-C., Hung, S.-T., Kuang, Y.-Y., Hseu, Y.-C., Fang, C.-L., et al. (2020). Aberrant Expression of NCAPG Is Associated with Prognosis and Progression of Gastric Cancer. Cancer Manag. Res. 12, 7837–7846. doi:10.2147/cmar.s248318

PubMed Abstract | CrossRef Full Text | Google Scholar

Sung, H., Ferlay, J., Siegel, R. L., Laversanne, M., Soerjomataram, I., Jemal, A., et al. (2021). Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA A. Cancer J. Clin. 71 (3), 209–249. doi:10.3322/caac.21660

CrossRef Full Text | Google Scholar

Szklarczyk, D., Gable, A. L., Lyon, D., Junge, A., Wyder, S., Huerta-Cepas, J., et al. (2019). STRING V11: Protein-Protein Association Networks with Increased Coverage, Supporting Functional Discovery in Genome-wide Experimental Datasets. Nucleic Acids Res. 47 (D1), D607–d613. doi:10.1093/nar/gky1131

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, Z., Li, C., Kang, B., Gao, G., Li, C., and Zhang, Z. (2017). GEPIA: a Web Server for Cancer and normal Gene Expression Profiling and Interactive Analyses. Nucleic Acids Res. 45 (W1), W98–w102. doi:10.1093/nar/gkx247

PubMed Abstract | CrossRef Full Text | Google Scholar

Vasaikar, S. V., Straub, P., Wang, J., and Zhang, B. (2018). LinkedOmics: Analyzing Multi-Omics Data within and across 32 Cancer Types. Nucleic Acids Res. 46 (D1), D956–d963. doi:10.1093/nar/gkx1090

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, L., Tong, X., Zhou, Z., Wang, S., Lei, Z., Zhang, T., et al. (2018). Circular RNA Hsa_circ_0008305 (circPTK2) Inhibits TGF-β-Induced Epithelial-Mesenchymal Transition and Metastasis by Controlling TIF1γ in Non-small Cell Lung Cancer. Mol. Cancer 17 (1), 140. doi:10.1186/s12943-018-0889-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Woroniecka, K., Chongsathidkiet, P., Rhodin, K., Kemeny, H., Dechant, C., Farber, S. H., et al. (2018). T-cell Exhaustion Signatures Vary with Tumor Type and Are Severe in Glioblastoma. Clin. Cancer Res. 24 (17), 4175–4186. doi:10.1158/1078-0432.ccr-17-1846

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, F., Huang, W., Yang, L., and Xu, F. (2019). MicroRNA-101-3p Regulates Gastric Cancer Cell Proliferation, Invasion and Apoptosis by Targeting PIM 1 Expression. Cel Mol Biol (Noisy-le-grand) 65 (7), 118–122. doi:10.14715/cmb/2019.65.7.20

CrossRef Full Text | Google Scholar

Xiao, C., Gong, J., Jie, Y., Cao, J., Chen, Z., Li, R., et al. (2020). NCAPG Is a Promising Therapeutic Target across Different Tumor Types. Front. Pharmacol. 11, 387. doi:10.3389/fphar.2020.00387

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, P., Jiang, L., Yang, Y., Wu, M., Liu, B., Shi, Y., et al. (2020). PAQR4 Promotes Chemoresistance in Non-small Cell Lung Cancer through Inhibiting Nrf2 Protein Degradation. Theranostics 10 (8), 3767–3778. doi:10.7150/thno.43142

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, T., Dong, M., Wang, Z., Li, H., and Li, X. (2020). Elevated mRNA Expression Levels of NCAPG Are Associated with Poor Prognosis in Ovarian Cancer. Cmar 12, 5773–5786. doi:10.2147/cmar.s253349

PubMed Abstract | CrossRef Full Text | Google Scholar

Xue, S., Hu, M., Iyer, V., and Yu, J. (2017). Blocking the PD-1/pd-L1 Pathway in Glioma: a Potential New Treatment Strategy. J. Hematol. Oncol. 10 (1), 81. doi:10.1186/s13045-017-0455-6

CrossRef Full Text | Google Scholar

Yu, G., Wang, L.-G., Han, Y., and He, Q.-Y. (2012). clusterProfiler: an R Package for Comparing Biological Themes Among Gene Clusters. OMICS: A J. Integr. Biol. 16 (5), 284–287. doi:10.1089/omi.2011.0118

CrossRef Full Text | Google Scholar

Zeng, Y., Shi, Y., Xu, L., Zeng, Y., Cui, X., Wang, Y., et al. (2021). Prognostic Value and Related Regulatory Networks of MRPL15 in Non-small-cell Lung Cancer. Front. Oncol. 11, 656172. doi:10.3389/fonc.2021.656172

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, X., Xie, K., Zhou, H., Wu, Y., Li, C., Liu, Y., et al. (2020). Role of Non-coding RNAs and RNA Modifiers in Cancer Therapy Resistance. Mol. Cancer 19 (1), 47. doi:10.1186/s12943-020-01171-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, X., Zhu, M., Wang, H., Song, Z., Zhan, D., Cao, W., et al. (2021). Overexpression of NCAPG Inhibits Cardia Adenocarcinoma Apoptosis and Promotes Epithelial-Mesenchymal Transition through the Wnt/β-Catenin Signaling Pathway. Gene 766, 145163. doi:10.1016/j.gene.2020.145163

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y. X., Yuan, J., Gao, Z. M., and Zhang, Z. G. (2018). LncRNA TUC338 Promotes Invasion of Lung Cancer by Activating MAPK Pathway. Eur. Rev. Med. Pharmacol. Sci. 22 (2), 443–449. doi:10.26355/eurrev_201801_14193

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, Y., Feng, C., Li, Y., Ma, Y., and Cai, R. (2019). LncRNA H19 Promotes Lung Cancer Proliferation and Metastasis by Inhibiting miR-200a Function. Mol. Cel Biochem 460 (1-2), 1–8. doi:10.1007/s11010-019-03564-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Zuberi, K., Franz, M., Rodriguez, H., Montojo, J., Lopes, C. T., Bader, G. D., et al. (2013). GeneMANIA Prediction Server 2013 Update. Nucleic Acids Res. 41, W115–W122. Web Server issue. doi:10.1093/nar/gkt533

PubMed Abstract | CrossRef Full Text | Google Scholar