Transcriptomic analysis reveals key lncRNAs associated with ribosomal biogenesis and epidermis differentiation in head and neck squamous cell carcinoma

Predicting the prognosis, immune response, and immunotherapy in head and neck squamous cell carcinoma using a novel risk model based on anoikis-related lncRNAs

BACKGROUND

Head and neck squamous cell carcinoma (HNSCC) is an extremely heterogeneous and metastatic disease. Anoikis, which is a specific type of programmed apoptosis, is involved in tumor metastasis, tissue homeostasis, and development. Herein, we constructed an anoikis-related long non-coding RNA (lncRNA) signature to predict the prognosis, immune responses, and therapeutic effects in HNSCC patients.

METHODS

A total of 501 HNSCC samples were acquired from the TCGA database and randomly classified into the training and validation groups (1:1 ratio). Thereafter, the results derived from the training set were analyzed with the LASSO regression analysis, and a novel anoikis-related lncRNA risk model was constructed. Time-dependent ROC curves and Kaplan-Meier analysis were carried out to assess the diagnostic value and survival outcomes. A nomogram was utilized to predict the prognostic accuracy. Furthermore, we studied the tumor microenvironment, tumor mutation burden, enrichment pathways, and the response to chemotherapy and immunotherapy.

RESULTS

Seven anoikis-related lncRNAs (AC015878.1, CYTOR, EMSLR, LINC01503, LINC02084, RAB11B-AS1, Z97200.1) were screened to design a novel risk model, which was recognized as the independent prognostic factor for HNSCC patients. The findings implied that low-risk patients showed significantly longer OS, PFS, and DSS compared to those high-risk patients. The two groups that were classified using the risk model showed significant differences in their immune landscape. The risk model also predicted that low-risk HNSCC patients could attain a better response to immunotherapy, while high-risk patients would be more sensitive to gemcitabine, docetaxel, and cisplatin.

CONCLUSIONS

We constructed a novel risk model that could be employed for effectively predicting patient prognosis with a good independent prognostic value for HNSCC patients. Furthermore, this model could be used for designing new immunotherapeutic and chemotherapeutic strategies, and it helps clinicians establish personalized and detailed strategies for HNSCC patients.

Identification of competing endogenous RNA network in laryngeal squamous cell carcinoma

OBJECTIVE

This study was conducted to investigate key long noncoding RNAs (lncRNAs) involved in competitive endogenous RNA (ceRNA) network associated with laryngeal squamous cell carcinoma (LSCC).

MATERIALS AND METHODS

Three mRNA datasets, two miRNA datasets, and one lncRNA dataset of LSCC were downloaded from GEO database. Following the identification of differentially expressed mRNAs (DEmRNAs), (microRNAs) miRNAs (DEmiRNAs), and lncRNAs (DElncRNAs) in LSCC compared with adjacent tissues, functional enrichment of DEmRNAs was performed. Then, construction of the ceRNA (DElncRNA-DEmiRNA-DEmRNA) regulatory network and functional analyses of all DEmRNAs in ceRNA regulatory network were conducted. Quantitative real-time polymerase chain reactions (qRT-PCR) were used to detect the expression levels of selected DEmRNAs, DEmiRNAs, and DElncRNAs.

RESULTS

A total of 3449 DEmRNAs, 40 DEmiRNAs, and 100 DElncRNAs were identified in LSCC. The ceRNA networks, which contained 132 DElncRNA-DEmiRNA pairs and 287 DEmiRNA-DEmRNA pairs, involving 44 lncRNAs, 3 miRNAs, and 271 mRNAs, were obtained. DEmRNAs in ceRNA regulatory networks were significantly enriched in pathways in cancer, prostate cancer, and aldosterone-regulated sodium reabsorption. Except for HCG22 and hsa-miR-1246, expressions of the others in the qRT-PCR results played the same pattern with that in our integrated analysis, generally.

CONCLUSIONS

We concluded that HCG22/EGOT-hsa-miR-1275-FAM107A and HCG22/EGOT-hsa-miR-1246-Glycerol-3-phosphate dehydrogenase 1 like interaction pairs may play a central role in LSCC.

Long non-coding RNA LINC00152 in cancer: Roles, mechanisms, and chemotherapy and radiotherapy resistance

Long non-coding RNA LINC00152 (cytoskeleton regulator, or LINC00152) is an 828-bp lncRNA located on chromosome 2p11.2. LINC00152 was originally discovered during research on hepatocarcinogenesis and has since been regarded as a crucial oncogene that regulates gene expression in many cancer types. LINC00152 is aberrantly expressed in various cancers, including gastric, breast, ovarian, colorectal, hepatocellular, and lung cancer, and glioma. Several studies have indicated that LINC00152 is correlated with cell proliferation, apoptosis, migration, invasion, cell cycle, epithelial-mesenchymal transition (EMT), chemotherapy and radiotherapy resistance, and tumor growth and metastasis. High LINC00152 expression in most tumors is significantly associated with poor patient prognosis. Mechanistic analysis has demonstrated that LINC00152 can serve as a competing endogenous RNA (ceRNA) by sponging miRNA, regulating the abundance of the protein encoded by a particular gene, or modulating gene expression at the epigenetic level. LINC00152 can serve as a diagnostic or prognostic biomarker, as well as a therapeutic target for most cancer types. In the present review, we discuss the roles and mechanisms of LINC00152 in human cancer, focusing on its functions in chemotherapy and radiotherapy resistance.

7SK Acts as an Anti-tumor Factor in Tongue Squamous Cell Carcinoma

Increasing evidence has shown the mechanistic insights about non-coding RNA 7SK in controlling the transcription. However, the biological function and mechanism of 7SK in cancer are largely unclear. Here, we show that 7SK is down-regulated in human tongue squamous carcinoma (TSCC) and acts as a TSCC suppressor through multiple cell-based assays including a migration assay and a xenograft mouse model. The expression level of 7SK was negatively correlated with the size of tumors in the 73 in-house collected TSCC patients. Through combined analysis of 7SK knockdown of RNA-Seq and available published 7SK ChIRP-seq data, we identified 27 of 7SK-regulated genes that were involved in tumor regulation and whose upstream regulatory regions were bound by 7SK. Motif analysis showed that the regulatory sequences of these genes were enriched for transcription factors FOXJ3 and THRA, suggesting a potential involvement of FOXJ3 and THRA in 7SK-regulated genes. Interestingly, the augmented level of FOXJ3 in TSCC patients and previous reports on THRA in other cancers have suggested that these two factors may promote TSCC progression. In support of this idea, we found that 21 out of 27 aforementioned 7SK-associated genes were regulated by FOXJ3 and THRA, and 12 of them were oppositely regulated by 7SK and FOXJ3/THRA. We also found that FOXJ3 and THRA dramatically promoted migration in SCC15 cells. Collectively, we identified 7SK as an antitumor factor and suggested a potential involvement of FOXJ3 and THRA in 7SK-mediated TSCC progression.

Integrative Analysis of DNA Methylation Data and Transcriptome Data Identified a DNA Methylation-Dysregulated Four-LncRNA Signature for Predicting Prognosis in Head and Neck Squamous Cell Carcinoma

Increasing evidence has demonstrated the crosstalk between DNA epigenetic alterations and aberrant expression of long non-coding RNAs (lncRNAs) during carcinogenesis. However, epigenetically dysregulated lncRNAs and their functional and clinical roles in Head and Neck Squamous Cell Carcinoma (HNSCC) are still not explored. In this study, we performed an integrative analysis of DNA methylation data and transcriptome data and identified a DNA methylation-dysregulated four-lncRNA signature (DNAMeFourLncSig) from 596 DNA methylation-dysregulated lncRNAs using a machine-learning-based feature selection method, which classified the patients of the discovery cohort into two risk groups with significantly different survival including overall survival, disease-specific survival, and progression-free survival. Then the DNAMeFourLncSig was implemented to another two HNSCC patient cohorts and showed similar prognostic values in both. Results from multivariable Cox regression analysis revealed that the DNAMeFourLncSig might be an independent prognostic factor. Furthermore, the DNAMeFourLncSig was substantially correlated with the complete response rate of chemotherapy and may predict chemotherapy response. Functional in silico analysis found that DNAMeFourLncSig-related mRNAs were mainly enriched in cell differentiation, tissue development and immune-related pathways. Overall, our study will improve our understanding of underlying transcriptional and epigenetic mechanisms in HNSCC carcinogenesis and provided a new potential biomarker for the prognosis of patients with HNSCC.

Assessment of multiple pathways involved in the inhibitory effect of HCG22 on oral squamous cell carcinoma progression

LncRNAs have been proposed to be associated with the tumorigenesis and progression of oral squamous cell carcinoma (OSCC). LncRNA HLA complex group 22 (HCG22) was reported to be lowly expressed and associated with poor prognosis in head and neck squamous cell carcinoma (HNSCC). However, the biological role and related mechanism of HCG22 in OSCC have not been characterized. HCG22 expression in OSCC cells was detected by qRT-PCR. Cell proliferation, invasion, and apoptosis were evaluated by Bromodeoxyuridine (BrdU) proliferation assay, Transwell invasion assay, and flow cytometry analysis, respectively. The protein levels of proliferating cell nuclear antigen (PCNA), E-cadherin, Vimentin, Bcl-2, Bax, protein kinase B (Akt), phosphorylated Akt (p-Akt), mammalian target of rapamycin (mTOR), phosphorylated mTOR (p-mTOR), and β-catenin were detected by western blot. Cell growth evaluation was performed using in vitro colony formation assay and in vivo tumor xenograft assay. We found that HCG22 was weakly expressed in OSCC cells. HCG22 overexpression inhibited cell proliferation and invasion and induced apoptosis in OSCC cells. The levels of PCNA, Vimentin, and Bcl-2 were decreased and E-cadherin and Bax expression was elevated in OSCC cells after HCG22 overexpression. Additionally, HCG22 overexpression inhibited the Akt, mTOR and Wnt/β-catenin pathways. Activation of Akt, mTOR, and Wnt/β-catenin pathways attenuated the anti-tumor property of HCG22 in OSCC cells. Furthermore, HCG22 overexpression inhibited the growth of OSCC cells in vitro and in vivo. In conclusion, HCG22 exerted anti-tumor property in OSCC by inhibiting the Akt, mTOR, and Wnt/β-catenin pathways.

MYCT1 Inhibits the Adhesion and Migration of Laryngeal Cancer Cells Potentially Through Repressing Collagen VI

MYCT1, a target of c-Myc, inhibits laryngeal cancer cell migration, but the underlying mechanism remains unclear. In the study, we detected differentially expressed genes (DEGs) from laryngeal cancer cells transfected by MYCT1 using RNA-seq (GSE123275). DEGs from head and neck squamous cell carcinoma (HNSCC) were first screened by comparison of transcription data from the Gene Expression Omnibus (GSE6631) and the Cancer Genome Atlas (TCGA) datasets using weighted gene co-expression network analysis (WGCNA). GO and KEGG pathway analysis explained the functions of the DEGs. The DEGs overlapped between GSE6631and TCGA datasets were then compared with ours to find the key DEGs downstream of MYCT1 related to the adhesion and migration of laryngeal cancer cells. qRT-PCR and Western blot were applied to validate gene expression at mRNA and protein levels, respectively. Finally, the cell adhesion, migration, and wound healing assays were to check cell adhesion and migration abilities, respectively. As results, 39 overlapping genes were enriched in the GSE6631 and TCGA datasets, and most of them revealed adhesion function. Thirteen of 39 genes including COL6 members COL6A1, COL6A2, and COL6A3 were overlapped in GSE6631, TCGA, and GSE123275 datasets. Similar to our RNA-seq results, we confirmed that COL6 is a target of MYCT1 in laryngeal cancer cells. We also found that MYCT1 inhibited the adhesion and migration of laryngeal cancer cells via COL6. These indicate that COL6 is a potential target of MYCT1 and participates the adhesion and migration of laryngeal cancer cells, which provides an important clue for further study on how MYCT1 regulating COL6 in laryngeal cancer progression.

Identification of PTPRR and JAG1 as key genes in castration-resistant prostate cancer by integrated bioinformatics methods. J Zhejiang Univ Sci B 2021;21:246-255. [PMID: 32133801 DOI: 10.1631/jzus.b1900329]

To identify novel genes in castration-resistant prostate cancer (CRPC), we downloaded three microarray datasets containing CRPC and primary prostate cancer in Gene Expression Omnibus (GEO). R packages affy and limma were performed to identify differentially expressed genes (DEGs) between primary prostate cancer and CRPC. After that, we performed functional enrichment analysis including gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) pathway. In addition, protein-protein interaction (PPI) analysis was used to search for hub genes. Finally, to validate the significance of these genes, we performed survival analysis. As a result, we identified 53 upregulated genes and 58 downregulated genes that changed in at least two datasets. Functional enrichment analysis showed significant changes in the positive regulation of osteoblast differentiation pathway and aldosterone-regulated sodium reabsorption pathway. PPI network identified hub genes like cortactin-binding protein 2 (CTTNBP2), Rho family guanosine triphosphatase (GTPase) 3 (RND3), protein tyrosine phosphatase receptor-type R (PTPRR), Jagged1 (JAG1), and lumican (LUM). Based on PPI network analysis and functional enrichment analysis, we identified two genes (PTPRR and JAG1) as key genes. Further survival analysis indicated a relationship between high expression of the two genes and poor prognosis of prostate cancer. In conclusion, PTPRR and JAG1 are key genes in the CRPC, which may serve as promising biomarkers of diagnosis and prognosis of CRPC.

RNA sequencing: new technologies and applications in cancer research

Over the past few decades, RNA sequencing has significantly progressed, becoming a paramount approach for transcriptome profiling. The revolution from bulk RNA sequencing to single-molecular, single-cell and spatial transcriptome approaches has enabled increasingly accurate, individual cell resolution incorporated with spatial information. Cancer, a major malignant and heterogeneous lethal disease, remains an enormous challenge in medical research and clinical treatment. As a vital tool, RNA sequencing has been utilized in many aspects of cancer research and therapy, including biomarker discovery and characterization of cancer heterogeneity and evolution, drug resistance, cancer immune microenvironment and immunotherapy, cancer neoantigens and so on. In this review, the latest studies on RNA sequencing technology and their applications in cancer are summarized, and future challenges and opportunities for RNA sequencing technology in cancer applications are discussed.

Multi-Omics Integration Reveals a Competitive Endogenous RNAs Network for the Identification of Progression Biomarkers and the Stratification of Patients Diagnosed With Nephroblastoma

Specific types of nephroblastoma (Wilms' tumor, WT) are known to associate with poor overall survival. Emerging experimental evidence has demonstrated that competitive endogenous RNA (ceRNA) networks have important roles in regulating cancer occurrence, but the roles of ceRNAs in regulating the WT progression and the patient outcomes remain unclear. Using the multi-omics data of 132 WT patients collected from TARGET database, an integration analysis pipeline was performed to construct a highly reliable ceRNA network. As results, a total of 147 nodes (116 mRNAs, 15 miRNAs, and 16 lncRNAs) were identified and used to explore the underlying mechanism for WT progression. WGCNA analysis further identified several prognostic molecules, including hsa-mir-93, LINC00087 and RP5-1086K13, that significantly associated with the overall survival rate. And, enrichment analysis verified the participation of these molecules in tumor-related pathways, such as those controlling autophagy and cadherin-mediated adhesion. Importantly, the WT patients were classified into three categories according to the ceRNA network, which significantly correlated with the overall survival. In conclusion, the ceRNA network could be a promising tool to further validate the prognostic biomarkers and categories of patients diagnosed with WT.

LncRNA CYTOR attenuates sepsis-induced myocardial injury via regulating miR-24/XIAP

This work aims to investigate the function and mechanism of long non-coding RNA (lncRNA) cytoskeleton regulator RNA (CYTOR) in myocardial injury induced by sepsis. The sepsis-induced myocardial injury model in mice was established by intraperitoneal injection of LPS (10 mg/kg) in vivo, and cardiomyocyte H9c2 was treated with LPS to mimic sepsis in vitro. CYTOR expression and miR-24 expression were detected by qRT-PCR. After up-regulation or down-regulation of CYTOR and miR-24 expression in the H9c2 cells, and the viability of the cells was detected via MTT assay, and cell apoptosis was detected by TUNEL assay. Western blot was applied to determine the expression level of caspase 3, Bax and X-chromosome-linked inhibitor of apoptosis (XIAP). Interaction between CYTOR and miR-24 was determined by bioinformatics analysis, RT-PCR and dual luciferase reporter assay. Interaction between miR-24 and XIAP was determined through bioinformatics analysis, RT-PCR, western blot and dual luciferase reporter assay. CYTOR was markedly down-regulated. CYTOR interacted with miR-24, and negatively regulated its expression level. Over-expression of CYTOR or transfection of miR-24 inhibitors significantly promoted viability and inhibited apoptosis of H9c2 cells, while the knockdown of CYTOR and transfection of miR-24 mimics had opposite effects. CYTOR suppressed the expression level of apoptosis-related proteins, but miR-24 increased them. miR-24 directly targeted the 3'UTR of XIAP, and suppressed it, and XIAP was modulated indirectly by CYTOR. Down-regulation of CYTOR aggravates sepsis-induced cardiac injury via regulating miR-24 and XIAP.

The cox-filter method identifies respective subtype-specific lncRNA prognostic signatures for two human cancers

BACKGROUND

The most common histological subtypes of esophageal cancer are squamous cell carcinoma (ESCC) and adenocarcinoma (EAC). It has been demonstrated that non-marginal differences in gene expression and somatic alternation exist between these two subtypes; consequently, biomarkers that have prognostic values for them are expected to be distinct. In contrast, laryngeal squamous cell cancer (LSCC) has a better prognosis than hypopharyngeal squamous cell carcinoma (HSCC). Likewise, subtype-specific prognostic signatures may exist for LSCC and HSCC. Long non-coding RNAs (lncRNAs) hold promise for identifying prognostic signatures for a variety of cancers including esophageal cancer and head and neck squamous cell carcinoma (HNSCC).

METHODS

In this study, we applied a novel feature selection method capable of identifying specific prognostic signatures uniquely for each subtype - the Cox-filter method - to The Cancer Genome Atlas esophageal cancer and HSNCC RNA-Seq data, with the objectives of constructing subtype-specific prognostic lncRNA expression signatures for esophageal cancer and HNSCC.

RESULTS

By incorporating biological relevancy information, the lncRNA lists identified by the Cox-filter method were further refined. The resulting signatures include genes that are highly related to cancer, such as H19 and NEAT1, which possess perfect prognostic values for esophageal cancer and HNSCC, respectively.

CONCLUSIONS

The Cox-filter method is indeed a handy tool to identify subtype-specific prognostic lncRNA signatures. We anticipate the method will gain wider applications.

Characterization of long non-coding RNA and messenger RNA profiles in laryngeal cancer by weighted gene co-expression network analysis

Laryngeal cancer (LC) is a malignant tumor in the head and neck region. It was recently elucidated that long non-coding RNAs (lncRNAs) participate in the pathogenesis of LC. However, the detailed mechanism of lncRNA in LC and whether long non-coding RNAs serve as effective biomarkers remains unclear. Ribonucleic acid (RNA) sequence data of LC and 11 patient clinical traits were extracted from The Cancer Genome Atlas (TCGA) database and analyzed by weighted gene co-expression network analysis (WGCNA). A total of 9 co-expression modules were identified. The co-expression Pink module significantly correlated with four clinical traits, including history of smoking, lymph node count, tumor status, and the success of follow-up treatment. Based on the co-expression Pink module, lncRNA-microRNA (miRNA)-messenger RNA (mRNA) and lncRNA-RNA binding protein-mRNA networks were constructed. We found that 8 lncRNAs significantly impacted overall survival (OS) in LC patients. These identified lncRNA and hub gene biomarkers were also validated in multiple LC cells in vitro via qPCR. Taken together, this study provided the framework of co-expression gene modules of LC and identified some important biomarkers in LC development and disease progression.

Role of cytoplasmic lncRNAs in regulating cancer signaling pathways

Cancer remains a serious healthcare problem despite significant improvements in early detection and treatment approaches in the past few decades. Novel biomarkers for diagnosis and therapeutic strategies are urgently needed. In recent years, long noncoding RNAs (lncRNAs) have been reported to be aberrantly expressed in tumors and show crosstalk with key cancer-related signaling pathways. In this review, we summarized the current progress of research on cytoplasmic lncRNAs and their roles in regulating cancer signaling and tumor progression, further characterization of which may lead to effective approaches for cancer prevention and therapy.

The challenge of determining the impact of FUT3 tumor-associated polymorphism rs2306969 (-6951 C> T) in invasive breast cancer cells

FUT3 gene is responsible for encode an homonymous α1,3/4-fucosyltransferase involved in the synthesis of sialyl-Lewis antigens. FUT3-fucosylated glycoconjugates play key roles in pathways involved in tumor biology and metastasis, such as cellular ligation to E-selectins, TGF-β-induced epithelial-mesenchymal transition, NK cell-mediated tumor cytotoxicity and apoptosis. Tumor-associated FUT3 promoter polymorphism rs2306969 (-6951 C> T, position related to the gene's translation start site) has been linked to breast, ovarian and intestinal gastric cancer. Although non-coding polymorphisms accounts for the majority of variations founded in breast cancer, their functional roles are still poorly understood. This study aimed to investigate the impact of different alleles for this variation in FUT3 expression of invasive breast tumors. A luciferase reporter assay was performed using two breast tumor cell lines to evaluate respectively the impact of FUT3 rs2306969 (-6951 CC) and (-6951 TT) on protein expression. Gene and protein expressions were also measured in twenty-nine fresh biopsies of invasive breast tumors. Rs2306969 did not significantly influence FUT3 expression in both used systems. However, this study is defiant since the biological role of this polymorphism in breast cancer and other tumor types could be linked to cis/trans modulation of other genes, respond to different environmental stimuli or impact gene expression only in association with other variations. Rs2306969 did not modulate FUT3 expression in breast tumors under non-stimulated conditions. Nevertheless, our study contributes to the notably challenging task that is to understand how non-coding polymorphisms can drive the overall risk in cancer development.