PDIA6 contributes to aerobic glycolysis and cancer progression in oral squamous cell carcinoma

The SOX2/PDIA6 axis mediates aerobic glycolysis to promote stemness in non-small cell lung cancer cells

Non-small cell lung cancer (NSCLC) is an aggressive and rapidly expanding lung cancer. Abnormal upregulation or knockdown of PDIA6 expression can predict poor prognosis in various cancers. This study aimed to investigate the biological function of PDIA6 in NSCLC. SOX2 and PDIA6 expression in NSCLC tissues and regulatory relationship between them were analyzed using bioinformatics. GSEA was performed on the enrichment pathway of PDIA6. qRT-PCR was utilized to examine expression of SOX2 and PDIA6 in NSCLC tissues and cells, and dual-luciferase reporter assay and ChIP experiments were performed to validate their regulatory relationship. CCK-8 experiment was conducted to assess cell viability, western blot was to examine levels of stem cell markers and proteins related to aerobic glycolysis pathway in cells. Cell sphere formation assay was used to evaluate efficiency of cell sphere formation. Reagent kits were used to measure glycolysis levels and glycolysis products. High expression of PDIA6 in NSCLC was linked to aerobic glycolysis. Knockdown of PDIA6 reduced cell viability, expression of stem cell surface markers, and cell sphere formation efficiency in NSCLC. Overexpression of PDIA6 could enhance cell viability and promote aerobic glycolysis, but the addition of 2-DG could reverse this result. Bioinformatics predicted the existence of upstream transcription factor SOX2 for PDIA6, and SOX2 was significantly upregulated in NSCLC, and they had a binding relationship. Further experiments revealed that PDIA6 overexpression restored repressive effect of knocking down SOX2 on aerobic glycolysis and cell stemness. This work revealed that the SOX2/PDIA6 axis mediated aerobic glycolysis to promote NSCLC cell stemness, providing new therapeutic strategies for NSCLC.

Population-enriched innate immune variants may identify candidate gene targets at the intersection of cancer and cardio-metabolic disease

Both cancer and cardio-metabolic disease disparities exist among specific populations in the US. For example, African Americans experience the highest rates of breast and prostate cancer mortality and the highest incidence of obesity. Native and Hispanic Americans experience the highest rates of liver cancer mortality. At the same time, Pacific Islanders have the highest death rate attributed to type 2 diabetes (T2D), and Asian Americans experience the highest incidence of non-alcoholic fatty liver disease (NAFLD) and cancers induced by infectious agents. Notably, the pathologic progression of both cancer and cardio-metabolic diseases involves innate immunity and mechanisms of inflammation. Innate immunity in individuals is established through genetic inheritance and external stimuli to respond to environmental threats and stresses such as pathogen exposure. Further, individual genomes contain characteristic genetic markers associated with one or more geographic ancestries (ethnic groups), including protective innate immune genetic programming optimized for survival in their corresponding ancestral environment(s). This perspective explores evidence related to our working hypothesis that genetic variations in innate immune genes, particularly those that are commonly found but unevenly distributed between populations, are associated with disparities between populations in both cancer and cardio-metabolic diseases. Identifying conventional and unconventional innate immune genes that fit this profile may provide critical insights into the underlying mechanisms that connect these two families of complex diseases and offer novel targets for precision-based treatment of cancer and/or cardio-metabolic disease.

CLSPN actives Wnt/β-catenin signaling to facilitate glycolysis and cell proliferation in oral squamous cell carcinoma

OBJECTIVE

Oral squamous cell carcinoma (OSCC) is a common malignancy with a poor prognosis. This study aimed to determine the influence and underlying mechanisms of CLSPN on OSCC.

METHODS

CLSPN expression was tested using quantitative real-time polymerase chain reaction, immunohistochemistry, and western blotting. Flow cytometry, cell counting kit, and colony formation assays were performed to determine OSCC cell apoptosis, viability, and proliferation, respectively. In OSCC cells, the extracellular acidification rate (ECAR), oxygen consumption rate (OCR), glucose uptake, and lactate production were determined using the corresponding kits. Changes in the protein levels of HK2, PKM2, LDHA, Wnt3a, and β-catenin were assessed using western blotting.

RESULTS

CLSPN expression was increased in OSCC tissues. Overexpression of CLSPN in HSC-2 cells promoted cell proliferation, increased the levels of ECAR, glucose uptake, and lactate production, and increased the protein levels of HK2, PKM2, LDHA, Wnt3a, and β-catenin, but inhibited OCR levels and apoptosis. The knockdown of CLSPN in CAL27 cells resulted in the opposite results. Moreover, the effects of CLSPN overexpression on glycolysis and OSCC cell proliferation were reversed by Wnt3a knockdown. In vivo, knockdown of CLSPN restrained tumor growth, glycolysis, and the activation of Wnt/β-catenin signaling.

CONCLUSION

CLSPN promoted glycolysis and OSCC cell proliferation, and reduced apoptosis, which was achieved by the activation of Wnt/β-catenin signaling pathway.

Comprehensive intratumoral heterogeneity landscaping of liver hepatocellular carcinoma and discerning of APLP2 in cancer progression

INTRODUCTION

As the sixth most common type of cancer worldwide, liver hepatocellular carcinoma (LIHC) emerges as grave public health danger owing to its chemotherapy-resistant feature. Disulfidoptosis is a newly discovered programmed cell death process affecting the normal actin cytoskeleton structure.

METHODS

Single-cell RNA (scRNA)-seq data were procured from GSE149614 and GSE202642 datasets. We utilized uniform manifold approximation and projection and clustering algorithm Louvian for dimensionality reduction and FindAllMarkers function for determining the differentially expressed genes (DEGs). Monocle2 and SCENIC were utilized to perform pseudo-time series and transcription factor analysis for selected subgroups. A series of in vitro experiments, including colony formation assay (CFA), flow cytometry targeting apoptosis and cell cycle, was applied to investigate how APLP2 regulated the LIHC progression. Two cell lines of LIHC cells, HepG2, and Huh7, were used for si-APLP2 transfection.

RESULTS

Tumor heterogeneity landscape of LIHC was depicted by detailed subgroup analysis. We found T and B cells were enriched with POU2F1 and HES1 activity. Inflammatory cancer-associated fibroblasts interacted with the cancer cells, uniquely through COL1A1/SDC1, COL1A2/SDC1 and LUM/ITGB1 pathways. The transformation from normal hepatocytes to malignant cells was displayed by cell trajectory analysis. State4, which was determined as malignant cells, was enriched in PI3K, hypoxia, and Epidermal growth factor receptor pathway, and enriched with Nuclear Receptor Subfamily 2 Group F Member 1 transcription factor activity. We observed an intense communication from the cancer cells to endothelial cells, mainly through the Vitronectin (VTN) to Kinase Insert Domain Receptor (KDR) pathway. A prognostic model targeting LIHC was constructed based on the disulfidoptosis-based DEGs, namely APLP2, PDIA6, YBX1, SPP1, whose accuracy was validated in multiple cohorts. Knockdown of APLP2 significantly increased the apoptosis and delayed cell cycle progression of LIHC cell line.

CONCLUSION

A prognostic model targeting LIHC was constructed based on the disulfidoptosis-related DEGs, which displayed high stability and accuracy in multiple cohorts. APLP2 played an active role in the carcinogenesis of LIHC by regulating the apoptosis and cell cycle.

Inhibition of UFM1 expression suppresses cancer progression and is linked to the dismal prognosis and immune infiltration in oral squamous cell carcinoma

BACKGROUND

Ubiquitin fold modifier 1 (UFM1) overexpression is associated with cancer cell proliferation, migration and invasion. However, the roles and pathways of UFM1 in oral squamous cell carcinoma (OSCC) has remained undefined.

METHODS

The expression of UFM1 and the relationship between UFM1 expression and prognosis were investigated using data of OSCC patients from The Cancer Genome Atlas (TCGA) database. The UFM1 co-expressed genes, and the association between the UFM1 expression and immune cells and ubiquitination were explored. The effects of UFM1 expression on the growth and migration of OSCC cells were investigated by siRNA interference, Cell Counting Kit-8 (CCK-8), Transwell, Western blotting, and wound healing experiments.

RESULTS

UFM1 was highly expressed in OSCC. UFM1 overexpression was associated with short overall survival, disease-specific survival, and progression-free interval, and was an adverse factor for prognosis in OSCC. UFM1-related nomograms were significantly associated with poor prognosis in OSCC patients. Decreased UFM1 expression could inhibit the proliferation, migration, and invasion of OSCC cells. UFM1 was associated with the immune cells (such as the Th17 cells, T helper cells, and cytotoxic cells) and ubiquitination.

CONCLUSION

Elevated UFM1 expression was associated with poor prognosis, ubiquitination and immune infiltration in OSCC, and inhibition of UFM1 expression delayed OSCC progression, showing that UFM1 could be a biomarker for prognosis and treating OSCC patients.

Hypoxia-induced tumor exosomes promote angiogenesis through miR-1825/TSC2/mTOR axis in oral squamous cell carcinoma

BACKGROUND

Oral squamous cell carcinoma (OSCC) is characterized by enhanced angiogenesis resulting in poor prognosis despite improvements in diagnostic/therapeutic techniques. Here, we aimed at investigating potential roles of miR-1825 enclosed in OSCC-derived exosomes on angiogenesis under hypoxic conditions.

METHODS

Effects of miR-1825 mimic/inhibitor as well as hypoxia-induced tumor derived exosomes on human umbilical vein endothelial cells (HUVECs) were evaluated using cell viability, migration/invasion, tube formation, and spheroid-based 3D angiogenesis assays.

RESULTS

Hypoxic conditions caused significant increase in miR-1825 levels in OSCC cells and hiTDEs. miR-1825 alone and within hiTDEs promoted endothelial cell viability, migration, invasion, and angiogenic potential, which is reversed via inhibition of miR-1825 expression. miR-1825 within hiTDEs altered the angiogenesis potential of HUVEC cells via deregulation of TSC2/mTOR axis.

CONCLUSIONS

We showed that hypoxia led to OSCC-derived exosome mediated transfer of miR-1825 to HUVECs and enhanced angiogenesis in OSCC in vitro.

Immunogenic landscape and risk score prediction based on unfolded protein response (UPR)-related molecular subtypes in hepatocellular carcinoma

Background

Hepatocellular carcinoma (HCC) is the most common type of cancer and causes a significant number of cancer-related deaths worldwide. The molecular mechanisms underlying the development of HCC are complex, and the heterogeneity of HCC has led to a lack of effective prognostic indicators and drug targets for clinical treatment of HCC. Previous studies have indicated that the unfolded protein response (UPR), a fundamental pathway for maintaining endoplasmic reticulum homeostasis, is involved in the formation of malignant characteristics such as tumor cell invasiveness and treatment resistance. The aims of our study are to identify new prognostic indicators and provide drug treatment targets for HCC in clinical treatment based on UPR-related genes (URGs).

Methods

Gene expression profiles and clinical information were downloaded from the TCGA, ICGC and GEO databases. Consensus cluster analysis was performed to classify the molecular subtypes of URGs in HCC patients. Univariate Cox regression and machine learning LASSO algorithm were used to establish a risk prognosis model. Kaplan-Meier and ROC analyses were used to evaluate the clinical prognosis of URGs. TIMER and XCell algorithms were applied to analyze the relationships between URGs and immune cell infiltration. Real time-PCR was performed to analyze the effect of sorafenib on the expression levels of four URGs.

Results

Most URGs were upregulated in HCC samples. According to the expression pattern of URGs, HCC patients were divided into two independent clusters. Cluster 1 had a higher expression level, worse prognosis, and higher expression of immunosuppressive factors than cluster 2. Patients in cluster 1 were more prone to immune escape during immunotherapy, and were more sensitive to chemotherapeutic drugs. Four key UPR genes (ATF4, GOSR2, PDIA6 and SRPRB) were established in the prognostic model and HCC patients with high risk score had a worse clinical prognosis. Additionally, patients with high expression of four URGs are more sensitive to sorafenib. Moreover, ATF4 was upregulated, while GOSR2, PDIA6 and SRPRB were downregulated in sorafenib-treated HCC cells.

Conclusion

The UPR-related prognostic signature containing four URGs exhibits high potential application value and performs well in the evaluation of effects of chemotherapy/immunotherapy and clinical prognosis.

Quercetin acts via the G3BP1/YWHAZ axis to inhibit glycolysis and proliferation in oral squamous cell carcinoma

There is increasing evidence that the GTPase-activating protein SH3 domain-binding protein 1 (G3BP1) plays important roles in the formation of various tumors. However, the biological functions and mechanism of G3BP1 in promoting the progression of oral squamous cell carcinoma remain largely unknown. The impacts of quercetin on glycolysis and proliferation of the CAL27 oral squamous cell carcinoma line were investigated, and the mediating role of the G3BP1/YWHAZ pathway was explored. CAL27 cells stably over- or underexpressing G3BP1 were treated with quercetin, and then cell proliferation was assayed together with the expression of proteins involved in glucose uptake, glycolysis, and lactate production, as well as the activity of hexokinase, pyruvate kinase, and lactate dehydrogenase. CAL27 cells expressed G3BP1 and YWHAZ at significantly higher levels than normal oral squamous cells. CAL27 cells showed the highest expression of both proteins among the three carcinoma lines (TSCCA, SCC15, 42 CAL27). Overexpressing G3BP1 in CAL27 cells markedly induced glucose uptake, glycolysis, cell proliferation, and YWHAZ expression. Knocking down G3BP1 or YWHAZ exerted the opposite effects, which were similar to the effects of inhibiting glycolysis. Quercetin repressed glucose uptake, glycolysis, cell proliferation, and G3BP1/YWHAZ signaling in a dose-dependent way, and these effects were antagonized by G3BP1 overexpression. Quercetin can inhibit glycolysis and cell proliferation of oral squamous cell carcinoma, apparently by inhibiting the G3BP1/YWHAZ axis.

Glycolysis-Related Gene Analyses Indicate That DEPDC1 Promotes the Malignant Progression of Oral Squamous Cell Carcinoma via the WNT/β-Catenin Signaling Pathway

Increasing evidence suggests that aerobic glycolysis is related to the progression of oral squamous cell carcinoma (OSCC). Hence, we focused on glycolysis-related gene sets to screen for potential therapeutic targets for OSCC. The expression profiles of OSCC samples and normal controls were obtained from The Cancer Genome Atlas (TCGA). Then, the differentially expressed gene sets were selected from the official GSEA website following extraction of the differentially expressed core genes (DECGs). Subsequently, we tried to build a risk model on the basis of DECGs to predict the prognosis of OSCC patients via Cox regression analysis. Furthermore, crucial glycolysis-related genes were selected to explore their biological roles in OSCC. Two active glycolysis-related pathways were acquired and 66 DECGs were identified. Univariate Cox regression analysis showed that six genes, including HMMR, STC2, DDIT4, DEPDC1, SLC16A3, and AURKA, might be potential prognostic factors. Subsequently, a risk formula consisting of DEPDC1, DDIT4, and SLC16A3 was established on basis of the six molecules. Furthermore, DEPDC1 was proven to be related to advanced stage cancer and lymph node metastasis. Moreover, functional experiments suggested that DEPDC1 promoted the aerobic glycolysis, migration, and invasion of OSCC via the WNT/β-catenin pathway. The risk score according to glycolysis-related gene expression might be an independent prognostic factor in OSCC. In addition, DEPDC1 was identified as playing a carcinogenic role in OSCC progression, suggesting that DEPDC1 might be a novel biomarker and therapeutic target for OSCC.

A Pan-Cancer Analysis of Heat-Shock Protein 90 Beta1(HSP90B1) in Human Tumours

Background: HSP90B1, a member of the heat-shock protein 90 family, plays a vital role as a molecular chaperone for oncogenes and stimulates tumour growth. However, its role in various cancers remains unexplored. Methods: Using the cancer genome atlas, gene expression omnibus the Human Protein Atlas databases and various other bioinformatic tools, this study investigated the involvement of HSP90B1 in 33 different tumour types. Results: The over-expression of HSP90B1 generally predicted poor overall survival and disease-free survival for patients with tumours, such as adrenocortical carcinoma, bladder urothelial carcinoma, kidney renal papillary cell carcinoma, and lung adenocarcinoma. In this study, HSP90B1 was highly expressed in the majority of tumours. A comparison was made between the phosphorylation of HSP90B1 in normal and primary tumour tissues, and putative functional mechanisms in HSP90B1-mediated oncogenesis were investigated. Additionally, the mutation burden of HSP90B1 in cancer was evaluated along with the survival rate of patients with cancer patients. Conclusion: This first pan-cancer investigation reveals the oncogenic functions of HSP90B1 in various cancers.

A signature of immune-related gene pairs (IRGPs) for risk stratification and prognosis of oral cancer patients

Background

With low response to present immunotherapy, it is imperative to identify new immune-related biomarkers for more effective immunotherapies for oral cancer.

Methods

RNA profiles for 390 oral cancer patients and 32 normal samples were downloaded from The Cancer Genome Atlas (TCGA) database and differentially expressed genes (DEGs) were analyzed. Immune genesets from ImmPort repository were overlapped with DEGs. After implementing univariate Cox analysis and the least absolute shrinkage and selection operator (LASSO) Cox regression analysis, key immune-related gene pairs (IRGPs) among the overlapped DEGs for predicting the survival risk were obtained. Then, the cutoff of risk score was calculated by the receiver operating characteristic (ROC) curve to stratify oral cancer patients into high and low-risk groups. Multivariate Cox analysis was used to analyze independent prognostic indicators for oral cancer. Besides, infiltration of immune cells, functional annotation, and mutation analysis of IRGPs were conducted. Biological functions correlated with IRGPs were enriched by Gene Set Enrichment Analysis (GSEA) method.

Results

We identified 698 differentially expressed genes (DEGs) in response to oral cancer. 17 IRGPs among the DEGs were identified and integrated into a risk score model. Patients in the high-risk group have a significantly worse prognosis than those in the low-risk group in both training (P<0.001) and test (P=0.019) cohorts. Meanwhile, the IRGP model was identified as an independent prognostic factor for oral cancer. Different infiltration patterns of immune cells were found between the high- and low-risk groups that more types of T and B cells were enriched in the low-risk group. More immune-related signaling pathways were highly enriched in the low-risk group and Tenascin C (TNC) was the most frequently mutated gene. We have developed a novel 17-IRGPs signature for risk stratification and prognostic prediction of oral cancer.

Conclusion

Our study provides a foundation for improved immunotherapy and prognosis and is beneficial to the individualized management of oral cancer patients.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12957-022-02630-1.

Phase I Trial of Cetuximab, Radiotherapy, and Ipilimumab in Locally Advanced Head and Neck Cancer

PURPOSE

Concurrent radiotherapy with cetuximab, an anti-EGFR mAb, is a standard treatment for locally advanced head and neck squamous carcinoma (HNSCC). Cytotoxic T lymphocyte antigen-4-positive (CTLA-4+) regulatory T cells (Treg) dampen cellular immunity and correlate negatively with clinical outcomes. This phase I study added ipilimumab, an anti-CTLA-4 mAb, to cetuximab-radiotherapy.

PATIENTS AND METHODS

A (3 + 3) design was used to establish the recommended phase II dose (RP2D) of ipilimumab, added at week 5 for four, every-3-week doses to fixed, standard cetuximab-radiotherapy. Eligible subjects had stage III to IVb, high-risk [human papillomavirus-negative (HPV-)] or intermediate-risk HPV-positive (HPV+)] HNSCC. Dose-limiting toxicity (DLT) was defined as any grade 4 adverse event (AE) except in-field radiation dermatitis or immune-related (ir) AE requiring ≥2 weeks of systemic steroids. Baseline tumor and serial blood specimens were collected for immune correlatives.

RESULTS

From July 2013 to May 2016, 18 patients enrolled. Two of 6 in cohort 1 (ipilimumab 3 mg/kg) experienced grade 3 dermatologic DLTs, triggering deescalation of ipilimumab to 1 mg/kg. Dose level -1 was expanded to N  =  12 without DLT. irAE included: grade 1, 2, and 3 dermatitis (2, 1, and 3 cases), grade 4 colitis (1), and grade 1 hyperthyroidism (1). Three-year disease-free survival (DFS) and overall survival were 72% [90% confidence interval (CI), 57-92] and 72% (90% CI, 56-92). High expression of coinhibitory receptors PD1/LAG3/CD39 on baseline tumor-infiltrating Treg was associated with worse DFS (HR = 5.6; 95% CI, 0.83-37.8; P = 0.08).

CONCLUSIONS

The RP2D for ipilimumab plus standard cetuximab-radiotherapy is 1 mg/kg in weeks 5, 8, 11, and 14. The regimen is tolerable and yields acceptable survival without cytotoxic chemotherapy.

Protein disulfide isomerase family 6 promotes the imatinib-resistance of renal cell carcinoma by regulation of Wnt3a-Frizzled1 axis

Imatinib is a nontoxic tyrosine kinase inhibitor, used in the treatment of advanced renal cell carcinoma. However, some patients with renal cell carcinoma develop resistance to imatinib. Protein disulfide isomerase family 6 (PDIA6) was involved in the chemo-resistance of lung adenocarcinoma. In this study, the effect of PDIA6 on imatinib-resistance of renal cell carcinoma was investigated. First, PDIA6 was found to be up-regulated in the imatinib-resistant renal cell carcinoma tissues and cells. Functional assays showed that knockdown of PDIA6 sensitized imatinib-resistant renal cell carcinoma cells to imatinib through decreasing the half-maximal inhibitory concentration (IC50) of imatinib-resistant renal cell carcinoma cells. Secondly, cell proliferation of imatinib-resistant renal cell carcinoma cells was suppressed by PDIA6 silencing, and the apoptosis was promoted with reduced Bcl-2, enhanced Bax and cleaved caspase-3. Moreover, the interference of PDIA6 increased phosphorylation of H2A histone family member X (γH2AX), while decreased Rad51 and phosphorylated DNA-dependent protein kinase (DNA-PK) (p-DNA-PK) in imatinib-resistant renal cell carcinoma cells. Lastly, protein expression levels of Wnt3a and Frizzled1 (FZD1) in imatinib-resistant renal cell carcinoma cells were down-regulated by silencing of PDIA6. Over-expression of FZD1 attenuated PDIA6 silencing-induced increase in cell apoptosis and decrease in cell proliferation in imatinib-resistant renal cell carcinoma cells. In conclusion, knockdown of PDIA6 sensitized imatinib-resistant renal cell carcinoma cells into imatinib through inactivation of Wnt3a-FZD1 axis.

Identification of a novel glycolysis-related signature to predict the prognosis of patients with breast cancer

Background

Breast cancer (BC) has a high incidence and mortality rate in females. Its conventional clinical characteristics are far from accurate for the prediction of individual outcomes. Therefore, we aimed to develop a novel signature to predict the survival of patients with BC.

Methods

We analyzed the data of a training cohort from the Cancer Genome Atlas (TCGA) database and a validation cohort from the Gene Expression Omnibus (GEO) database. After the applications of Gene Set Enrichment Analysis (GSEA) and Cox regression analyses, a glycolysis-related signature for predicting the survival of patients with BC was developed; the signature contained AK3, CACNA1H, IL13RA1, NUP43, PGK1, and SDC1. Furthermore, on the basis of expression levels of the six-gene signature, we constructed a risk score formula to classify the patients into high- and low-risk groups. The receiver operating characteristic (ROC) curve and the Kaplan-Meier curve were used to assess the predicted capacity of the model. Later, a nomogram was developed to predict the outcomes of patients with risk score and clinical features over a period of 1, 3, and 5 years. We further used Human Protein Atlas (HPA) database to validate the expressions of the six biomarkers in tumor and sample tissues, which were taken as control.

Results

We constructed a six-gene signature to predict the outcomes of patients with BC. The patients in the high-risk group showed poor prognosis than those in the low-risk group. The area under the curve (AUC) values were 0.719 and 0.702, showing that the prediction performance of the signature is acceptable. Additionally, Cox regression analysis revealed that these biomarkers could independently predict the prognosis of BC patients with BC without being affected by clinical factors. The expression levels of all six biomarkers in BC tissues were higher than that in normal tissues; however, AK3 was an exception.

Conclusion

We developed a six-gene signature to predict the prognosis of patients with BC. Our signature has been proved to have the ability to make an accurate prediction and might be useful in expanding the hypothesis in clinical research.

Comprehensive analysis of angiogenesis subtype of squamous cell carcinoma

Background

Squamous cell carcinoma (SCC) is a disease with distinct management complexities as it displays a remarkably heterogeneous molecular subtype. However, the landscape of angiogenesis for SCC is not fully investigated.

Method and materials

The angiogenesis-related subtypes of SCC were established by using the ConsensusClusterPlus package based on angiogenesis-related genes and TCGA data. We analyzed the alteration of genes and miRNAs as well as pathways associated with angiogenesis subtypes. Next, the regulation network, the correlation with genomic characteristics, immune microenvironment, and clinical features of the angiogenesis subtypes were further investigated. Finally, the prognostic impact of the angiogenesis-related subtypes for SCC was also analyzed.

Results

A total of 1368 SCC samples were included in this study. Two angiogenesis subtypes were then identified based on the one hundred and sixty-three angiogenesis-related genes with subtype1 (angiogenesis subtype) of 951 SCC patients and subtype2 (non-angiogenesis subtype) of 417 SCC. GSEA revealed that angiogenesis and epithelial-mesenchymal transition, inflammatory response, and hypoxia were enriched in the angiogenesis subtype. Eight of the 15 immune checkpoints (ADORA2A, BTLA, CD276, CYBB, HAVCR2, SIGLEC7, SIGLEC9, and VTCN1) were significantly upregulated while C10orf54 were significantly downregulated in the angiogenesis subtype. The survival analysis revealed that the patients in the angiogenesis subtype have poorer survival outcomes than those in the non-angiogenesis subtype (P = 0.017 for disease-free interval and P = 0.00013 for overall survival).

Conclusion

Our analysis revealed a novel angiogenesis subtype classification in SCC and provides new insights into a hallmark of SCC progression.

Long non-coding RNA PCAT6 regulates bladder cancer progression via the microRNA-143-3p/PDIA6 axis

Although long non-coding (lnc)RNAs have been reported to be involved in the pathological development of bladder cancer, the functions of lncRNA prostate cancer-associated transcript 6 (PCAT6) and its underlying mechanism of action in bladder cancer remain unknown. The present study aimed to investigate the effect of PCAT6 in bladder cancer progression and explore its potential application as a novel treatment target. The expression of PCAT6 and miR-143-3p in bladder cancer tissues, adjacent normal tissues and cell lines was measured using reverse transcription-quantitative PCR. Fluorescence in situ hybridization assay was used to detect the subcellular localization of PCAT6. MTT, EdU, Transwell and wound healing assays were conducted to assess the biological function of PCAT6 on cell proliferation, migration and invasion. Putative binding sites between miR-143-3p and PCAT6 or PDIA6 were predicted using starBase, Lncbase and TargetScan analyzes. Dual-luciferase reporter assay was also used to confirm the potential binding between PCAT6 and miR-143-3p. RNA immunoprecipitation assay was performed to verify the possible interaction between PCAT6 and miR-143-3p. Western blotting was used to measure the expression of PDIA6. The results demonstrated that the expression levels of PCAT6 were upregulated in bladder cancer tissues relative to those in adjacent normal bladder tissues. Knockdown of PCAT6 served a role in suppressing the proliferation, migration and invasion of T24T and EJ bladder cancer cells. PCAT6 knockdown contributed to a reduction of PDIA6 expression at the mRNA and protein levels compared with that in negative control-transfected cells, whilst the miR-143-3p inhibitor partially mitigated this reduction effect. In addition, rescue experiments revealed that the miR-143-3p inhibitors reversed the effects of PCAT6 silencing on the malignant phenotypes of bladder cancer. Collectively, the results of the present study demonstrated that PCAT6 may serve an oncogenic role in bladder cancer via the miR-143-3p/PDIA6 axis. These results may provide a potential therapeutic target for the treatment of bladder cancer.