Integrative analyses of transcriptome sequencing identify novel functional lncRNAs in esophageal squamous cell carcinoma

The transcriptional landscape and clinico-biological characterization of human endogenous retroviruses in esophageal squamous cell carcinoma

Human endogenous retroviruses (HERVs) are emerging as critical elements in host genomic regulation. Aberrant HERV transcription has been implicated in developmental and tissue-specific aging and pathological processes. In this study, we presented a comprehensive locus-specific characterization of the HERV expression landscape in esophageal squamous cell carcinoma (ESCC). We demonstrated the transcriptional diversity among patients and identified 12 clinically relevant HERVs in the SCH cohort, which were experimentally validated by Real-Time Quantitative Polymerase Chain Reaction (RT-qPCR) in the CAMS cohort. ESCC patients were stratified into three HERV-based subtypes (HERVhigh, HERVmedian and HERVlow) with distinct clinical and biological characteristics. The HERVhigh subtype was associated with worse survival, increased CD4+ T cells infiltration and decreased metabolic activity, whereas the HERVlow subtype was characterized by abundant CD8+ T cells, increased metabolic activity, and better survival. The HERV-based tumor subtyping was further robustly validated by RNA sequencing and RT-qPCR in two additional external cohorts. Our findings demonstrate the clinical significance of HERVs for tumor subtyping and prognosis, provide insights into the functional role of HERVs and a valuable resource for developing novel biomarkers and therapeutic targets in ESCC.

An integrated analysis of the anticarcinogenic role of forkhead box protein 1 in oesophageal squamous cell carcinoma

Forkhead box protein 1 (FOXP1) serves as a tumour promoter or suppressor depending on different cancers, but its effect in oesophageal squamous cell carcinoma has not been fully elucidated. This study investigated the role of FOXP1 in oesophageal squamous cell carcinoma through bioinformatics analysis and experimental verification. We determined through public databases that FOXP1 expresses low in oesophageal squamous cell carcinoma compared with normal tissues, while high expression of FOXP1 indicates a better prognosis. We identified potential target genes regulated by FOXP1, and explored the potential biological processes and signalling pathways involved in FOXP1 in oesophageal squamous cell carcinoma through GO and KEGG enrichment, gene co-expression analysis, and protein interaction network construction. We also analysed the correlation between FOXP1 and tumour immune infiltration levels. We further validated the inhibitory effect of FOXP1 on the proliferation of oesophageal squamous cell carcinoma cells through CCK-8, colony formation and subcutaneous tumour formation assays. This study revealed the anticarcinogenic effect of FOXP1 in oesophageal squamous cell carcinoma, which may serve as a novel biological target for the treatment of tumour.

RSPO2-associated mitochondrial metabolism defines molecular subtypes with distinct clinical and immune features in esophageal cancer

BACKGROUND

Esophageal cancer is a highly aggressive malignancy with limited treatment options and poor prognosis. The identification of novel molecular subtypes and therapeutic targets is crucial for improving clinical outcomes.

METHOD

In this study, we investigated the role of R-spondin 2 (RSPO2) in esophageal cancer and its association with mitochondrial metabolism. Using bioinformatics analysis of publicly available datasets, we identified a panel of RSPO2-related mitochondrial metabolism genes and their expression patterns in esophageal cancer. Based on these genes, we stratified esophageal cancer patients into distinct molecular subtypes with different survival rates, immune cell infiltration profiles, and drug sensitivities.

RESULTS

Our findings suggest that RSPO2-related mitochondrial metabolism genes may serve as potential therapeutic targets and prognostic markers for esophageal cancer. These genes play an important role in the prognosis, immune cell infiltration and drug sensitivity of esophageal cancer.

CONCLUSION

The identified molecular subtypes provide valuable insights into the underlying molecular mechanisms of esophageal cancer and could guide personalized treatment strategies in the future.

Comprehensive landscape and future perspectives of non-coding RNAs in esophageal squamous cell carcinoma, a bibliometric analysis from 2008 to 2023

Objectives: Summarize the progress and hot topic evolution of non-coding RNAs (ncRNAs) research in esophageal squamous cell carcinoma (ESCC) in recent years and predict future research directions. Methods: Relevant articles from the Web of Science until 31 October 2023 were obtained. Bibliometric analysis of included articles was performed using software (VOSviewer, CiteSpace, and Bibliometrix). The volume and citation of publications, as well as the country, institution, author, journal, keywords of the articles were used as variables to analyze the research trends and hot spot evolution. Results: 1,118 literature from 2008 to 2023 were retrieved from database, with 25 countries/regions, 793 institutions, 5,426 authors, 261 journals involved. Global cooperation was centered on China, Japan, and the United States. Zhengzhou University, an institution from China, had the highest publication. The most prolific author was Guo Wei, and the most prolific journal was Oncology Letters. Analysis of keywords revealed that the research in this field revolved around the role of ncRNAs in the occurrence, development, diagnosis, treatment, and prognosis of ESCC, mainly including micro RNAs, long non-coding RNAs, and then circular RNAs. Conclusion: Overall, research on ncRNAs in ESCC remains strong. Previous research has mainly focused on the basic research, with a focus on the mechanism of ncRNAs in the occurrence, development, diagnosis, treatment, and prognosis of ESCC. Combining current research with emerging disciplines to further explore its mechanisms of action or shifting the focus of research from preclinical research to clinical research based on diagnosis, treatment, and prognosis, will be the main breakthrough in this field in the future.

Validation of serum cystatin SN detection for diagnosis and poor prognosis of esophageal squamous cell carcinoma

Background

The identification of effective tumor markers is of paramount importance for the early diagnosis, treatment, and prognosis of esophageal squamous cell carcinoma (ESCC). The present study endeavors to identify efficacious serological markers that can differentiate patients with early-stage ESCC from those with benign esophageal lesions and healthy controls (HC). Cystatin-SN (CST1), an active cysteine protease inhibitor belonging to the Cystatin (CST) superfamily, is implicated in the pathogenesis of inflammation and tumorigenesis. The objective of this investigation is to assess the diagnostic, therapeutic, and prognostic potential of serum CST1 in ESCC.

Methods

In our prior RNA sequencing and screening endeavors, we have identified ten genes that are up-regulated in relation to esophageal cancer. Subsequently, we have verified the gene CST1 from the transcriptome data of the The Cancer Genome Atlas Program (TCGA) and Gene Expression Profiling Interactive Analysis (GEPIA) database. Following this, we conducted an enzyme-linked immunosorbent assay (ELISA) to ascertain the expression levels of CST1 in serum samples from clinical cohorts.

Results

The study revealed a significant elevation in serum CST1 levels among patients with early-stage esophageal squamous cell carcinoma (ESCC) (7.41 ± 4.32 ng/ml) compared to those with esophageal benign lesions (4.67 ± 2.43 ng/ml) (p < 0.0001) and healthy controls (4.87 ± 2.77 ng/ml) (p < 0.0001). The diagnostic sensitivity of CST1 for ESCC was 75.68% (specificity 70.83%, AUC 0.775). Combination of CST1 and SCC-Ag exhibited the AUC up to 0.819. Additionally, serum CST1 levels exhibited a significant decrease at 1-2 weeks post-surgery (4.49 ± 3.31 ng/ml) compared to pre-surgery levels (7.68 ± 3.71 ng/ml) (p<0.0001). Survival analysis demonstrated a strong association between high (844/415-1543 d) or low (1490/645-1710 d) serum CST1 levels at diagnosis and overall survival time (p < 0.001). Furthermore, multivariate regression analysis confirmed CST1 (p=0.024, HR=2.023, 95%CI 1.099-3.725) as an independent prognostic factor.

Conclusion

Serum CST1 has the potential to function as a diagnostic indicator for distinguishing early-stage esophageal squamous cell carcinoma (ESCC) from individuals with benign esophageal lesions and healthy individuals. Additionally, it could serve as a prognostic predictor and therapeutic efficacy indicator for patients with ESCC.

Comparative transcriptome characterization of esophageal squamous cell carcinoma and adenocarcinoma

Background

Esophageal cancers are primarily categorized as esophageal squamous cell carcinoma (ESCC) and esophageal adenocarcinoma (EAC). While various (epi) genomic alterations associated with tumor development in ESCC and EAC have been documented, a comprehensive comparison of the transcriptomes in these two cancer subtypes remains lacking.

Methods

We collected 551 gene expression profiles from publicly available sources, including normal, ESCC, and EAC tissues or cell lines. Subsequently, we conducted a systematic analysis to compare the transcriptomes of these samples at various levels, including gene expression, promoter activity, alternative splicing (AS), alternative polyadenylation (APA), and gene fusion.

Results

Seven distinct cluster gene expression patterns were identified among the differentially expressed genes in normal, ESCC, and EAC tissues. These patterns were enriched in the PI3K-Akt signaling pathway and the activation of extracellular matrix organization and exhibited repression of epidermal development. Notably, we observed additional genes or unique expression levels enriched in these shared pathways and biological processes related to tumor development and immune activation. In addition to the differentially expressed genes, there was an enrichment of lncRNA co-expression networks and downregulation of promoter activity associated with the repression of epidermal development in both ESCC and EAC. This indicates a common feature between these two cancer subtypes. Furthermore, differential AS and APA patterns in ESCC and EAC appear to partially affect the expression of host genes associated with bacterial or viral infections in these subtypes. No gene fusions were observed between ESCC and EAC, thus highlighting the distinct molecular mechanisms underlying these two cancer subtypes.

Conclusions

We conducted a comprehensive comparison of ESCC and EAC transcriptomes and uncovered shared and distinct transcriptomic signatures at multiple levels. These findings suggest that ESCC and EAC may exhibit common and unique mechanisms involved in tumorigenesis.

Molecular Characterization of Esophageal Squamous Cell Carcinoma Using Quantitative Proteomics

Esophageal squamous cell carcinoma (ESCC) is a heterogeneous cancer associated with a poor prognosis in advanced stages. In India, it is the sixth most common cause of cancer-related mortality. In this study, we employed high-resolution mass spectrometry-based quantitative proteomics to characterize the differential protein expression pattern associated with ESCC. We identified several differentially expressed proteins including PDPN, TOP2A, POSTN and MMP2 that were overexpressed in ESCC. In addition, we identified downregulation of esophagus tissue-enriched proteins such as SLURP1, PADI1, CSTA, small proline-rich proteins such as SPRR3, SPRR2A, SPRR1A, KRT4, and KRT13, involved in squamous cell differentiation. We identified several overexpressed proteins mapped to the 3q24-29 chromosomal region, aligning with CNV alterations in this region reported in several published studies. Among these, we identified overexpression of SOX2, TP63, IGF2BP2 and RNF13 that are encoded by genes in the 3q26 region. Functional enrichment analysis revealed proteins involved in cell cycle pathways, DNA replication, spliceosome, and DNA repair pathways. We identified the overexpression of multiple proteins that play a major role in alleviating ER stress, including SYVN1 and SEL1L. The SYVN1/SEL1L complex is an essential part of the ER quality control machinery clearing misfolded proteins from the ER. SYVN1 is an E3 ubiquitin ligase that ubiquitinates ER-resident proteins. Interestingly, there are also other non-canonical substrates of SYVN1 which are known to play a crucial role in tumor progression. Thus, SYVN1 could be a potential therapeutic target in ESCC.

NIMA-related kinase-6 (NEK6) as an executable target in cancer

Cancer is a disease that develops when cells begin to divide uncontrollably and spreads to other parts of the body. Proliferation and invasion of cancerous cells are generally known to be influenced by cell cycle-related proteins in human malignancies. Therefore, in this review, we have emphasized on the serine/threonine kinase named NEK6. NEK6 is been deliberated to play a critical role in mitosis progression that includes mitotic spindle formation, metaphase to anaphase transition, and centrosome separation. Moreover, it has a mechanistic role in DNA repair and can cause apoptosis when inhibited. Past studies have connected NEK6 protein expression to cancer cell senescence. Besides, there are reports relating NEK6 to a range of malignancies including breast, lung, ovarian, prostate, kidney, liver, and others. Given its significance, this review attempts to describe the structural and functional aspects of NEK6 in various cellular processes, as well as how it is linked to different forms of cancer. Lastly, we have accentuated, on some of the plausible inhibitors that have been explored against NEK6 overexpression.

A preliminary analysis of prognostic genes in advanced laryngeal squamous cell carcinoma patients with postoperative radiotherapy

Advanced laryngeal squamous cell carcinoma (LSCC) has a high mortality rate, and the prognosis is poor. However, the underlying molecular biological mechanisms bringing about the development and progression of advanced LSCC are not entirely clarified. This study aimed to find out the potential biomarkers to predict the prognosis in advanced LSCC patients who had undergone postoperative radiotherapy alone. The next-generation sequencing of RNA was performed to detect the mRNAs expression profiling in 10 advanced LSCC samples, comprised of 5 samples from LSCC patients with favorable outcome and 5 samples from paired patients with poor outcome. Then bioinformatics analysis including Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway were used to find out functional core genes that were significantly different between the two groups. 1630 differentially expressed genes (DEGs) were confirmed to have significant differences between the two groups. 53 GO terms and 19 pathways which were closely related to the DEGs were identified. Finally, 52 intersection DEGs which were both related to the top three GO terms and pathways were identified. The expression of several core genes was confirmed with RT-qPCR in tissues from another 75 patients. RT-qPCR confirmed that the genes of c-JUN, LYN, PIK3R2, and TNFAIP3 were significantly differentially expressed between the two groups, which was in accordance with the RNA sequencing data. The DEGs identified above may be potential prognostic markers for advanced LSCC patients with postoperative radiotherapy, and may provide essential guidance for following-up.

Systematic identification of aberrant non-coding RNAs and their mediated modules in rotator cuff tears

Background: Rotator cuff tears (RCT) is the most common cause of shoulder dysfunction, however, its molecular mechanisms remain unclear. Non-coding RNAs(ncRNAs), such as long ncRNA (lncRNA), microRNA (miRNA) and circular RNA (circRNA), are involved in a variety of diseases, but little is known about their roles in RCT. Therefore, the purpose of this study is to identify dysregulated ncRNAs and understand how they influence RCT.

Methods: We performed RNA sequencing and miRNA sequencing on five pairs of torn supraspinatus muscles and matched unharmed subscapularis muscles to identify RNAs dysregulated in RCT patients. To better comprehend the fundamental biological processes, we carried out enrichment analysis of these dysregulated mRNAs or the co-expressed genes of dysregulated ncRNAs. According to the competing endogenous RNA (ceRNA) theory, we finally established ceRNA networks to explore the relationship among dysregulated RNAs in RCT.

Results: A total of 151 mRNAs, 38 miRNAs, 20 lncRNAs and 90 circRNAs were differentially expressed between torn supraspinatus muscles and matched unharmed subscapularis muscles, respectively. We found that these dysregulated mRNAs, the target mRNAs of these dysregulated miRNAs or the co-expressed mRNAs of these dysregulated ncRNAs were enriched in muscle structure development, actin-mediated cell contraction and actin binding. Then we constructed and analyzed the ceRNA network and found that the largest module in the ceRNA network was associated with vasculature development. Based on the topological properties of the largest module, we identified several important ncRNAs including hsa_circ_0000722, hsa-miR-129-5p and hsa-miR-30c-5p, whose interacting mRNAs related to muscle diseases, fat and inflammation.

Conclusion: This study presented a systematic dissection of the expression profile of mRNAs and ncRNAs in RCT patients and revealed some important ncRNAs which may contribute to the development of RCT. Such results could provide new insights for further research on RCT.

Emerging Roles of lncRNAs Regulating RNA-Mediated Type-I Interferon Signaling Pathway

The type-I interferon (IFN-I) signaling pathway plays pivot roles in defending against pathogen invasion. Exogenous ssRNA and dsRNA could be immunogenic. RNA-mediated IFN signaling is extensively studied in the field. The incorrect functioning of this pathway leads to either autoimmune diseases or suffering from microorganism invasion. From the discrimination of “self” and “non-self” molecules by receptors to the fine-tune modulations in downstream cascades, all steps are under the surveillance featured by complex feedbacks and regulators. Studies in recent years highlighted the emerging roles of long noncoding RNAs (lncRNAs) as a reservoir for signaling regulation. LncRNAs bind to targets through the structure and sequence, and thus the mechanisms of action can be complex and specific. Here, we summarized lncRNAs modulating the RNA-activated IFN-I signaling pathway according to the event order during the signaling. We hope this review help understand how lncRNAs are participating in the regulation of IFN-I signaling.

Long non-coding RNAs in gastrointestinal cancers: implications for protein phosphorylation

Phosphorylation of proteins is one of the most extensively investigated post-translational protein modifications. Threonine, serine and tyrosine in proteins are the most commonly phosphorylated amino acids. Dysregulated cancer-related signaling pathways due to aberrant phosphorylation status of the key protein(s) in these pathways exist in most malignancies. Intensive studies in the recent decade have implicated long non-coding RNAs (lncRNAs) in the precise regulation of protein phosphorylation in cancers. In this review, we systematically delve into recent advance that underlines the multidimensional role of lncRNAs in modulating protein phosphorylation, regulating cancerous signaling and impacting prognosis of gastrointestinal (GI) cancers including hepatocellular carcinoma, colorectal cancer, gastric cancer, esophageal cancer, and pancreatic cancer. LncRNAs regulate protein phosphorylation via directly binding to the target protein(s), interacting with the partner protein(s) of the target protein(s) or lncRNAs-encoded small peptides. Although there are still extensive studies on disclosing the intricate interactions between lncRNAs and proteins and their impacts on protein phosphorylation, we believe that targeting lncRNAs controlling phosphorylation of key protein(s) in cancerous signaling pathways might provide novel paths for precision therapeutics of GI cancers in the future.

CRNDE/ETS1/GPR17 Facilitates the Proliferation, Migration, and Invasion of Glioma

BACKGROUND

Numerous lncRNAs were found as regulatory factors for occurrence and progression of various tumors, but there is still less research on the role of lncRNAs in malignant progression of glioma.

METHODS

Bioinformatics analysis analyzed differential genes (DEGs) in the TCGA database. MTT, flow cytometry, and Transwell assays were performed to test the proliferation, apoptosis, migration, and invasion of cells. qRT-PCR and western blot were conducted to detect RNA and protein expressions of each gene, respectively. CHIP assay verified the binding relationship between genes. FISH assayed subcellar location of CRNDE, and xenograft in nude mice was performed for in vivo verification.

RESULTS

CRNDE was upregulated in glioma cells, and overexpression of CRNDE facilitated malignant progression of glioma cells. CRNDE regulated occurrence and development of glioma through the CRNDE-ETS1-GPR17 axis. ETS1 was proved to target promoter region of GPR17. Overexpression of CRNDE promoted the binding between ETS1 and the promoter region of GPR17, thus, promoting the transcription of GPR17, while silencing of GPR17 inhibited promotion of CRNDE on proliferation, migration, and invasion of glioma cells.

CONCLUSIONS

These results demonstrated that CRNDE regulated GPR17 expression by binding ETS1, a transcription factor, thereby affecting glioma development. The results also indicated that CRNDE could serve as a possible therapeutic target and prognostic biomarker for glioma.

FoXA2 promotes esophageal squamous cell carcinoma progression by ZEB2 activation

Background

It has been reported that Forkhead transcription family member (FOXA2) regulates esophageal squamous cell carcinoma (ESCC) progression. However, the specific mechanism, by which FOXA2 promotes ESCC malignant progression, remains unclear.

Materials and methods

QRT-PCR and western blotting were applied to measure FOXA2 expression in ESCC tissues, while CCK-8 assay and Transwell assays were used to investigate the effect of FOXA2 on ESCC. Luciferase reporter assay, followed by fast chromatin immunoprecipitation (ChIP) assay, was used to study the relationship between FOXA2 and ZEB2.

Results

FOXA2 was significantly increased in ESCC tissues, when compared to normal tissues. Moreover, high expression of FOXA2 was also found in ESCC cells. Knockdown of FOXA2 inhibited ESCC cell proliferation, invasion, and migration. Mechanically, FOXA2 was verified to regulate ZEB2 expression at transcription level. Moreover, ZEB2 reversed the inhibitory effect of FOXA2 on ESCC proliferation, invasion, and migration. The relationship between ZEB2 and FOXA2 in ESCC tissues was negative.

Conclusions

These results indicate that FOXA2 plays a critical role in ESCC progression and may become a potential candidate target for ESCC treatment.

Interplay and cooperation between SREBF1 and master transcription factors regulate lipid metabolism and tumor-promoting pathways in squamous cancer

Squamous cell carcinomas (SCCs) comprise one of the most common histologic types of human cancer. Transcriptional dysregulation of SCC cells is orchestrated by tumor protein p63 (TP63), a master transcription factor (TF) and a well-researched SCC-specific oncogene. In the present study, both Gene Set Enrichment Analysis (GSEA) of SCC patient samples and in vitro loss-of-function assays establish fatty-acid metabolism as a key pathway downstream of TP63. Further studies identify sterol regulatory element binding transcription factor 1 (SREBF1) as a central mediator linking TP63 with fatty-acid metabolism, which regulates the biosynthesis of fatty-acids, sphingolipids (SL), and glycerophospholipids (GPL), as revealed by liquid chromatography tandem mass spectrometry (LC-MS/MS)-based lipidomics. Moreover, a feedback co-regulatory loop consisting of SREBF1/TP63/Kruppel like factor 5 (KLF5) is identified, which promotes overexpression of all three TFs in SCCs. Downstream of SREBF1, a non-canonical, SCC-specific function is elucidated: SREBF1 cooperates with TP63/KLF5 to regulate hundreds of cis-regulatory elements across the SCC epigenome, which converge on activating cancer-promoting pathways. Indeed, SREBF1 is essential for SCC viability and migration, and its overexpression is associated with poor survival in SCC patients. Taken together, these data shed light on mechanisms of transcriptional dysregulation in cancer, identify specific epigenetic regulators of lipid metabolism, and uncover SREBF1 as a potential therapeutic target and prognostic marker in SCC.

Long non-coding RNA OIP5-AS1 promotes the progression of esophageal cancer by regulating miR-30a/VOPP1 expression

Long non-coding RNAs (lncRNAs) serve an important role in the development of esophageal cancer (EC), which is the eighth most common type of cancer worldwide. lncRNA opa-interacting protein 5 antisense transcript 1 (OIP5-AS1) is associated with human malignancy. However, the biological roles of OIP5-AS1 in the development of EC remain unclear. In the present study, transfection was conducted, and reverse transcription-quantitative PCR and western blot analysis were used for the detection of mRNA and protein expression, respectively. Furthermore, dual-luciferase reporter and RNA immunoprecipitation assays were used to study the interaction between miRNA and lncRNA or genes. The results revealed that OIP5-AS1 expression in EC tissues and cultured EC cells was upregulated, microRNA-30a (miR-30a) expression was downregulated. OIP5-AS1-knockdown suppressed the proliferation, migration and invasion of EC9706 and EC109 cells. miR-30a was confirmed to interact with OIP5-AS1, and miR-30a-mimics transfection ameliorated the effects of OIP5-AS1 in EC cells. Vesicular overexpressed in cancer prosurvival protein 1 (VOPP1) was verified as the direct target of miR-30a. VOPP1 expression was positively correlated with OIP5-AS1 expression in EC cells. Overexpression of VOPP1 ameliorated the negative effects of OIP5-AS1-knockdown on EC9706 and EC109 cells. In conclusion, OIP5-AS1 promoted the proliferation, migration and invasion of EC cells by increasing VOPP1 expression by sponging miR-30a.

linc01305 promotes metastasis and proliferation of esophageal squamous cell carcinoma through interacting with IGF2BP2 and IGF2BP3 to stabilize HTR3A mRNA

Evidence shows that long noncoding RNAs (lncRNAs) modulate mRNAs of multiple genes by post-transcriptional regulation. However, in esophageal squamous cell carcinoma, lncRNAs involvement in post-transcriptional regulation of mRNAs have been rarely reported. In this study, we investigated a novel mechanism of linc01305 promoting metastasis and proliferation of ESCC. The results for real-time quantitative reverse transcription PCR (qRT-PCR) and fluorescence in situ hybridization showed that linc01305 was highly expressed and predominantly located in cytoplasm of human esophageal cancer cells. Transwell and colony formation assays confirmed that linc01305 promoted migration and proliferation of esophageal cancer cells. RNA-seq, linc01305 pulldown, mass spectrometry, RNA immunoprecipitation and mRNA stability assays demonstrated that linc01305 stabilized mRNA of target gene HTR3A through interacting with IGF2BP2 and IGF2BP3. Taken together, our data unveils a novel mechanism in which cytoplasmic linc01305 stabilizes HTR3A mRNA through interacting with IGF2BP2 and IGF2BP3 and thereby promotes metastasis and proliferation of ESCC.

Upregulated long non‑coding RNA LincIN promotes tumor progression via the regulation of nuclear factor 90/microRNA‑7/HOXB13 in esophageal squamous cell carcinoma

Long non-coding RNA LincIN has been reported to be overexpressed and to be involved in the metastasis of breast cancer. However, the expression and role of LincIN in esophageal squamous cell carcinoma (ESCC) remain unsolved. In the present study, LincIN expression was examined in ESCC by RT-qPCR, and the roles of LincIN in ESCC were determined using cell growth, migration and invasion assays. In addition, the effects of LincIN on nuclear factor 90 (NF90) and microRNA/miR (miR)-7 were examined by RNA immunoprecipitation assay, RT-qPCR, dual-luciferase reporter assay and western blot analysis. The results revealed that LincIN expression was significantly increased in ESCC tissues and cell lines. The increased expression of LincIN was positively associated with invasion depth, lymph node metastasis, TNM stage and a poor prognosis. Functional assays revealed that the overexpression of LincIN promoted ESCC cell growth, migration and invasion. Mechanistic analysis revealed that LincIN physically bound to NF90, enhanced the binding between NF90 and primary miR-7 (pri-miR-7), and further enhanced the inhibitory effects of NF90 on miR-7 biogenesis. Therefore, LincIN downregulated miR-7 expression in ESCC. The expression of miR-7 inversely correlated with that of LincIN in ESCC tissues. By downregulating miR-7, LincIN increased the expression of HOXB13, a target of miR-7. The overexpression of miR-7 or the depletion of HOXB13 both attenuated the tumor-promoting roles of LincIN in ESCC cell growth, migration and invasion. On the whole, the findings of the present study suggest that LincIN is overexpressed and plays an oncogenic role in ESCC via the regulation of the NF90/miR-7/HOXB13 axis. Thus, LincIN may prove to be a promising prognostic biomarker and therapeutic target for ESCC.

Involvement of long non-coding RNAs in the progression of esophageal cancer

Esophageal cancer (EC) is one of the most common malignant tumors of the digestive system with high incidence and mortality rate worldwide. Therefore, exploring the pathogenesis of EC and searching for new targeted therapies are the current research hotspot for EC treatment. Long non‐coding RNAs (lncRNAs) are endogenous RNAs with more than 200 nucleotides, but without protein‐coding function. In recent years, lncRNAs have gradually become the focuses in the field of non‐coding RNA. Some lncRNAs have been proved to be closely related to the pathogenesis of EC. Many lncRNAs are abnormally expressed in EC and participate in many biological processes including cell proliferation, apoptosis, and metastasis by inhibiting or promoting target gene expression. LncRNAs can also regulate the progression of EC through epithelial‐mesenchymal transformation (EMT), which is closely related to the occurrence, development, and prognosis of EC. In this article, we review and discuss the involvement of lncRNAs in the progression of EC.

Alterations of RNA splicing patterns in esophagus squamous cell carcinoma

Alternative splicing (AS) is an important biological process for regulating the expression of various isoforms from a single gene and thus to promote proteome diversity. In this study, RNA-seq data from 15 pairs of matched esophageal squamous cell carcinoma (ESCC) and normal tissue samples as well as two cell lines were analyzed. AS events with significant differences were identified between ESCC and matched normal tissues, which were re-annotated to find protein coding genes or non-coding RNAs. A total of 45,439 AS events were found. Of these, 6019 (13.25%) significant differentially AS events were identified. Exon skipping (SE) events occupied the largest proportion of abnormal splicing events. Fifteen differential splicing events with the same trends of ΔΨ values in ESCC tissues, as well in the two cell lines were found. Four pathways and 20 biological processes related to pro-metastasis cell junction and migration were significantly enriched for the differentially spliced genes. The upregulated splicing factor SF3B4, which regulates 92 gene splicing events, could be a potential prognostic factor of ESCC. Differentially spliced genes, including HNRNPC, VCL, ZNF207, KIAA1217, TPM1 and CALD1 are shown with a sashimi plot. These results suggest that cell junction- and migration-related biological processes are influenced by AS abnormalities, and aberrant splicing events can be affected by splicing factor expression changes. The involved splicing factor SF3B4 was found to be a survival-related gene in ESCC and is presumed to regulate AS in multiple cancers. In summary, we identified significant differentially expressed AS events which may be related to the development of ESCC.

Systematic Proteome and Lysine Succinylome Analysis Reveals Enhanced Cell Migration by Hyposuccinylation in Esophageal Squamous Cell Carcinoma

Esophageal squamous cell carcinoma (ESCC) is an aggressive malignancy with poor therapeutic outcomes. However, the alterations in proteins and posttranslational modifications (PTMs) leading to the pathogenesis of ESCC remain unclear. Here, we provide the comprehensive characterization of the proteome, phosphorylome, lysine acetylome, and succinylome for ESCC and matched control cells using quantitative proteomic approach. We identify abnormal protein and PTM pathways, including significantly downregulated lysine succinylation sites in cancer cells. Focusing on hyposuccinylation, we reveal that this altered PTM was enriched on enzymes of metabolic pathways inextricably linked with cancer metabolism. Importantly, ESCC malignant behaviors such as cell migration are inhibited once the level of succinylation was restored in vitro or in vivo. This effect was further verified by mutations to disrupt succinylation sites in candidate proteins. Meanwhile, we found that succinylation has a negative regulatory effect on histone methylation to promote cancer migration. Finally, hyposuccinylation is confirmed in primary ESCC specimens. Our findings together demonstrate that lysine succinylation may alter ESCC metabolism and migration, providing new insights into the functional significance of PTM in cancer biology.

Silencing FAM135B enhances radiosensitivity of esophageal carcinoma cell

FAM135B (family with sequence similarity 135, member B) is related to the progression of esophageal squamous cell carcinoma (ESCC). However, the role played by the gene in radiosensitivity remains unknown. Herein, we examined the relationship between FAM135B and radiosensitivity. According to the results, FAM135B is highly expressed in ESCC cells, and ESCC cells with high levels of FAM135B are resistant to irradiation. Silencing FAM135B inhibits colony formation capability and cell cycle protein expression (pP53, CDK1), promotes cell cycle arrest at the G2/M phase following irradiation. Moreover, transcriptome sequencing analysis demonstrates that FAM135B regulates downstream PI3K/Akt/mTOR signaling pathway, and western blot verifies the result. One of the mechanisms of increasing radiosensitivity by silencing FAM135B expression in ESCC cells may be achieved by regulating the PI3K/Akt/mTOR signaling pathway. Silencing FAM135B shows synergy with PI3K/Akt/mTOR pathway inhibitor (rapamycin) in increasing radiosensitivity, regulating the expression of cell cycle protein and inducing apoptosis of ESCC cells. The results indicate that FAM135B could be a potential treatment target for ESCC in management of radiosensitivity.

Multi-Omics Analysis to Characterize Cigarette Smoke Induced Molecular Alterations in Esophageal Cells

Though smoking remains one of the established risk factors of esophageal squamous cell carcinoma, there is limited data on molecular alterations associated with cigarette smoke exposure in esophageal cells. To investigate molecular alterations associated with chronic exposure to cigarette smoke, non-neoplastic human esophageal epithelial cells were treated with cigarette smoke condensate (CSC) for up to 8 months. Chronic treatment with CSC increased cell proliferation and invasive ability of non-neoplastic esophageal cells. Whole exome sequence analysis of CSC treated cells revealed several mutations and copy number variations. This included loss of high mobility group nucleosomal binding domain 2 (HMGN2) and a missense variant in mediator complex subunit 1 (MED1). Both these genes play an important role in DNA repair. Global proteomic and phosphoproteomic profiling of CSC treated cells lead to the identification of 38 differentially expressed and 171 differentially phosphorylated proteins. Bioinformatics analysis of differentially expressed proteins and phosphoproteins revealed that most of these proteins are associated with DNA damage response pathway. Proteomics data revealed decreased expression of HMGN2 and hypophosphorylation of MED1. Exogenous expression of HMGN2 and MED1 lead to decreased proliferative and invasive ability of smoke exposed cells. Immunohistochemical labeling of HMGN2 in primary ESCC tumor tissue sections (from smokers) showed no detectable expression while strong to moderate staining of HMGN2 was observed in normal esophageal tissues. Our data suggests that cigarette smoke perturbs expression of proteins associated with DNA damage response pathways which might play a vital role in development of ESCC.

Characterization of super-enhancer-associated functional lncRNAs acting as ceRNAs in ESCC

Long noncoding RNAs (lncRNAs) have important regulatory roles in cancer biology. Although some lncRNAs have well-characterized functions, the vast majority of this class of molecules remains functionally uncharacterized. To systematically pinpoint functional lncRNAs, a computational approach was proposed for identification of lncRNA-mediated competing endogenous RNAs (ceRNAs) through combining global and local regulatory direction consistency of expression. Using esophageal squamous cell carcinoma (ESCC) as model, we further identified many known and novel functional lncRNAs acting as ceRNAs (ce-lncRNAs). We found that most of them significantly regulated the expression of cancer-related hallmark genes. These ce-lncRNAs were significantly regulated by enhancers, especially super-enhancers (SEs). Landscape analyses for lncRNAs further identified SE-associated functional ce-lncRNAs in ESCC, such as HOTAIR, XIST, SNHG5, and LINC00094. THZ1, a specific CDK7 inhibitor, can result in global transcriptional downregulation of SE-associated ce-lncRNAs. We further demonstrate that a SE-associated ce-lncRNA, LINC00094 can be activated by transcription factors TCF3 and KLF5 through binding to SE regions and promoted ESCC cancer cell growth. THZ1 downregulated expression of LINC00094 through inhibiting TCF3 and KLF5. Our data demonstrated the important roles of SE-associated ce-lncRNAs in ESCC oncogenesis and might serve as targets for ESCC diagnosis and therapy.

Proteomic variations of esophageal squamous cell carcinoma revealed by combining RNA-seq proteogenomics and G-PTM search strategy

BACKGROUND

Cancer that arises from epithelial cells of the esophagus is called esophagus squamous cell carcinoma (ESCC) and is mostly observed in developing nations. Evaluation of cancer genomes and its regulation into proteins plays a predominant role in understanding the cancer progressions. Mass-spectrometry-based proteomics is a consequential tool to estimate proteomic variation and posttranslational modifications (PTMs) from standard protein databases. Post-translational modifications play a crucial role in protein folding and PTMs can be accounted for as a biological signal to interpret the structural changes and transition order of proteins. Functional validation of cancer-related mutations can explain the effects of mutations on genes and the identification of Oncogenes and tumor suppressor genes. Therefore, we present a study on protein variations to interpret the structural changes and transition order of proteins in ESCC carcinogenesis.

METHODOLOGY

We are using a bottom-up proteomics approach with Galaxy-P framework and RNA sequence data analysis to generate the sample-specific databases containing details of RNA splicing and variant peptides. Once the database generated with information on variable modification, only the curated PTMs at specific positions are considered to perform spectral matching. Proteogenomics mapping was performed to identify protein variations in ESCC.

RESULTS

RNA-sequence proteogenomics with G-PTM (Global Post-Translational Modification) searching strategy has revealed proteomic events including several peptides that contain single amino acid variations, novel splice junction peptides and posttranslationally modified peptides. Proteogenomic mapping exhibited the splice junction peptides mapped predominantly for Malic enzyme exon type (ME-3) and MCM7 protein-coding genes that promote cancer progression, found to be exhibited in ESCC samples. Approximately 25 ± types of PTM modifications were recorded, and Protein Phosphorylation was largely noted.

CONCLUSION

ESCC cancer prognosis at the molecular level enables a better understanding of cancer carcinogenesis and protein modifications can be used as potential biomarkers.

Long non-coding RNAs in ovarian granulosa cells

Granulosa cells (GCs) are somatic cells surrounding oocytes within follicles and are essential for folliculogenesis. Pathological changes in GCs are found in several ovarian disorders. Recent reports have indicated that long non-coding RNAs (lncRNAs), which modulate gene expression via multiple mechanisms, are key regulators of the normal development of GCs, follicles, and ovaries. In addition, accumulating evidence has suggested that lncRNAs can be utilized as biomarkers for the diagnosis and prognosis of GC-related diseases, such as polycystic ovary syndrome (PCOS) and premature ovarian insufficiency (POI). Therefore, lncRNAs not only play a role in GCs that are involved in normal folliculogenesis, but they may also be considered as potential candidate biomarkers and therapeutic targets in GCs under pathological conditions. In the future, a detailed investigation of the in vivo delivery or targeting of lncRNAs and large-cohort-validation of the clinical applicability of lncRNAs is required.

Identification of lncRNA-associated differential subnetworks in oesophageal squamous cell carcinoma by differential co-expression analysis

Differential expression analysis has led to the identification of important biomarkers in oesophageal squamous cell carcinoma (ESCC). Despite enormous contributions, it has not harnessed the full potential of gene expression data, such as interactions among genes. Differential co-expression analysis has emerged as an effective tool that complements differential expression analysis to provide better insight of dysregulated mechanisms and indicate key driver genes. Here, we analysed the differential co-expression of lncRNAs and protein-coding genes (PCGs) between normal oesophageal tissue and ESCC tissues, and constructed a lncRNA-PCG differential co-expression network (DCN). DCN was characterized as a scale-free, small-world network with modular organization. Focusing on lncRNAs, a total of 107 differential lncRNA-PCG subnetworks were identified from the DCN by integrating both differential expression and differential co-expression. These differential subnetworks provide a valuable source for revealing lncRNA functions and the associated dysfunctional regulatory networks in ESCC. Their consistent discrimination suggests that they may have important roles in ESCC and could serve as robust subnetwork biomarkers. In addition, two tumour suppressor genes (AL121899.1 and ELMO2), identified in the core modules, were validated by functional experiments. The proposed method can be easily used to investigate differential subnetworks of other molecules in other cancers.

RNA sequencing of exosomes revealed differentially expressed long noncoding RNAs in early-stage esophageal squamous cell carcinoma and benign esophagitis

Aim: To explore the roles of exosomal long noncoding RNAs (lncRNAs) in early-stage esophageal squamous cell carcinoma (ESCC) and benign esophagitis. Materials & methods: Exosomal lncRNAs were analyzed using RNA-seq and validated by quantitative real-time PCR, loss-of-function, co-culture and RNA pulldown assays. Results: Exosomal lncRNAs displayed tighter tissue-specificity, higher expression level and lower splicing efficiency than that of mRNAs. A total of 152 exosomal lncRNAs were differentially expressed between ESCC and controls. A total of 124 exosomal lncRNAs were dysregulated between ESCC and esophagitis. Knockdown of 13 ESCC-associated lncRNAs modified proliferation, migration, and apoptosis of ESCC cells. A novel lncRNA RP5-1092A11.2 was highly expressed in ESCC-derived exosomes, ESCC cells and tumor tissues. Exosomes released from RP5-1092A11.2-knockdown cells inhibited ESCC cell proliferation. Conclusion: Dysregulated exosomal lncRNAs were functionally associated with different disease status in esophagus.

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.

Identification of key genes by integrating DNA methylation and next-generation transcriptome sequencing for esophageal squamous cell carcinoma

Aberrant DNA methylation leads to abnormal gene expression, making it a significant regulator in the progression of cancer and leading to the requirement for integration of gene expression with DNA methylation. Here, we identified 120 genes demonstrating an inverse correlation between DNA methylation and mRNA expression in esophageal squamous cell carcinoma (ESCC). Sixteen key genes, such as SIX4, CRABP2, and EHD3, were obtained by filtering 10 datasets and verified in paired ESCC samples by qRT-PCR. 5-Aza-dC as a DNA methyltransferase (DNMT) inhibitor could recover their expression and inhibit clonal growth of cancer cells in seven ESCC cell lines. Furthermore, 11 of the 16 genes were correlated with OS (overall survival) and DFS (disease-free survival) in 125 ESCC patients. ChIP-Seq data and WGBS data showed that DNA methylation and H3K27ac histone modification of these key genes displayed inverse trends, suggesting that there was collaboration between DNA methylation and histone modification in ESCC. Our findings illustrate that the integrated multi-omics data (transcriptome and epigenomics) can accurately obtain potential prognostic biomarkers, which may provide important insight for the effective treatment of cancers.

Identification of key genes in invasive clinically non-functioning pituitary adenoma by integrating analysis of DNA methylation and mRNA expression profiles

BACKGROUND

Tumor surrounding the internal carotid artery or invading to the cavernous sinus is an important characteristic of invasive pituitary adenoma, and a pivotal factor of tumor residue and regrowth. Without specific changes in serum hormone related to the adenohypophyseal cell of origin, clinically non-functioning pituitary adenoma is more likely to be diagnosed at invasive stages compared with functioning pituitary adenoma. The underlying mechanism of tumor invasion remains unknown. In this study, we aimed to identify key genes in tumor invasion by integrating analyses of DNA methylation and gene expression profiles.

METHOD

Genome-wide DNA methylation and mRNA microarray analysis were performed for tumor samples from 68 patients at the Beijing Tiantan Hospital. Differentially expressed genes and methylated probes were identified based on an invasive vs non-invasive grouping. Differentially methylated probes in the promoter region of targeted genes were assessed. Pearson correlation analysis was used to identify genes with a strong association between DNA methylation status and expression levels. Pyrosequencing and RT-PCR were used to validate the methylation status and expression levels of candidate genes, respectively.

RESULTS

A total of 8842 differentially methylated probes, located on 4582 genes, and 661 differentially expressed genes were identified. Both promoter methylation and expression alterations were observed for 115 genes with 58 genes showing a negative correlation between DNA methylation status and expression level. Nineteen genes that exhibited notably negative correlations between DNA methylation and gene expression levels, are involved in various gene ontologies and pathways, or played an important role in different diseases, were regarded as candidate genes. We found an increased methylation with a decreased expression of PHYHD1, LTBR, C22orf42, PRR5, ANKDD1A, RAB13, CAMKV, KIFC3, WNT4 and STAT6, and a decreased methylation with an increased expression of MYBPHL. The methylation status and expression levels of these genes were validated by pyrosequencing and RT-PCR.

CONCLUSIONS

The DNA methylation and expression levels of PHYHD1, LTBR, MYBPHL, C22orf42, PRR5, ANKDD1A, RAB13, CAMKV, KIFC3, WNT4 and STAT6 are associated with tumor invasion, and these genes may become the potential genes for targeted therapy.

LncRNA625 inhibits STAT1-mediated transactivation potential in esophageal cancer cells

Although Signal transducer and activator of transcription 1 (STAT1)-mediated transactivation potential is inhibited in cancer cells, the mechanism is poorly understood. In the present study, we implicated long non-coding RNA lncRNA625 in the inhibition of STAT1 activity. LncRNA625 knockdown up-regulated STAT1-mediated transcription and resulted in an increase of STAT1-mediated expression of IFITM2. Conversely, lncRNA625 upregulation inhibited STAT1 reporter activity. Mechanistically, lncRNA625 inhibited STAT1 binding to the promoter of IFITM2 in both untreated cells and following interferon-gamma (IFN-γ) stimulation. LncRNA625 interacted with the DNA-binding (DB) domain of STAT1 and promoted STAT1 interaction with T-cell protein tyrosine phosphatase TC45 to dephosphorylate pSTAT1. Taken together, the results show that lncRNA625 inhibits STAT1-mediated transactivation potential by causing formation of STAT1-TC45 complexes, resulting in STAT1 dephosphorylation.

LncRNA PICSAR promotes cell proliferation, migration and invasion of fibroblast-like synoviocytes by sponging miRNA-4701-5p in rheumatoid arthritis

Background

Long non-coding RNAs (lncRNAs) have drawn increasing attention because they play a pivotal role in various types of autoimmune diseases, including rheumatoid arthritis (RA). Fibroblast-like synoviocytes (FLSs), a prominent component of hyperplastic synovial pannus tissue, are the primary effector cells in RA synovial hyperplasia and invasion which can lead to joint destruction. In this study, we investigated whether lncRNAs could act as competing endogenous RNAs to regulate the pathological behaviors of RA-FLSs.

Methods

LncRNA microarray was conducted to establish lncRNA expression profiles in FLSs isolated from RA patients and healthy controls (HCs). Differentially expressed lncRNAs were verified by quantitative real-time PCR (qRT-PCR) on RA-FLSs and synovial fluid. The functional role of lncRNA PICSAR downregulation was evaluated in RA-FLSs. We conducted molecular biological analysis to predict miRNAs which have a potential binding site for PICSAR and further refined the results by qRT-PCR. Luciferase reporter assay was adopted to validate the interaction of lncRNA PICSAR and miR-4701-5p. Western Blot and qPCR were used to identify the target gene and protein. The functional role of miR-4701-5p upregulation was examined in RA-FLSs.

Findings

We identified a long intergenic non-protein-coding RNA162 (LINC00162), also known as lncRNA PICSAR (p38 inhibited cutaneous squamous cell carcinoma associated lincRNA), has significantly higher expression in RA-FLSs and RA synovial fluid. The cell proliferation, migration, invasion and proinflammatory cytokines production of RA-FLSs showed significant alterations after the lncRNA PICSAR suppression. Mechanistically, lncRNA PICSAR functioned through sponging miR-4701-5p in RA-FLSs.

Interpretation

Our results reveal PICSAR may exert an essential role in promoting synovial invasion and joint destruction by sponging miR-4701-5p in RA and that lncRNA PICSAR may act as a biomarker of RA.

HERES, a lncRNA that regulates canonical and noncanonical Wnt signaling pathways via interaction with EZH2

Aberrant lncRNA expression is responsible for cancer progression and metastasis, positioning lncRNAs not only as biomarkers but also as promising therapeutic targets for curing cancer. A number of lncRNAs have been reported in ESCC but their mechanistic roles largely remain unknown. Wnt signaling pathways are often dysregulated in ESCC; however, the role of lncRNAs in such dysregulation was also undetermined. We found 6 lncRNAs that are significantly dysregulated and correlated with outcomes in ESCC patients. The most upregulated lncRNA, HERES, promotes cancer progression and epigenetically regulates canonical and noncanonical Wnt signaling pathways simultaneously through interaction with EZH2. These results show that HERES represents an early diagnostic and therapeutic target for squamous-cell-type cancers caused by defects in Wnt signaling pathways.

Wnt signaling through both canonical and noncanonical pathways plays a core role in development. Dysregulation of these pathways often causes cancer development and progression. Although the pathways independently contribute to the core processes, a regulatory molecule that commonly activates both of them has not yet been reported. Here, we describe a long noncoding RNA (lncRNA), HERES, that epigenetically regulates both canonical and noncanonical Wnt signaling pathways in esophageal squamous cell carcinoma (ESCC). For this study, we performed RNA-seq analysis on Korean ESCC patients and validated these results on a larger ESCC cohort to identify lncRNAs commonly dysregulated in ESCCs. Six of the dysregulated lncRNAs were significantly associated with the clinical outcomes of ESCC patients and defined 4 ESCC subclasses with different prognoses. HERES reduction repressed cell proliferation, migration, invasion, and colony formation in ESCC cell lines and tumor growth in xenograft models. HERES appears to be a transacting factor that regulates CACNA2D3, SFRP2, and CXXC4 simultaneously to activate Wnt signaling pathways through an interaction with EZH2 via its G-quadruple structure-like motif. Our results suggest that HERES holds substantial potential as a therapeutic target for ESCC and probably other cancers caused by defects in Wnt signaling pathways.

RAC1 inhibition reverses cisplatin resistance in esophageal squamous cell carcinoma and induces downregulation of glycolytic enzymes

Development of chemoresistance remains a major challenge in treating esophageal squamous cell carcinoma (ESCC) patients despite treatment advances. However, the role of RAC1 in chemoresistance of ESCC and the underlying mechanisms remain largely unknown. In this study, we found that higher levels of RAC1 expression were associated with poorer prognosis in ESCC patients. Enhanced RAC1 expression increased cell proliferation, migration, and chemoresistance in vitro. Combination therapy using RAC1 inhibitor EHop-016 and cisplatin significantly promoted cell viability inhibition, G2/M phase cycle arrest, and apoptosis when compared to each monotherapy. Mechanistically, glycolysis was significantly downregulated in the RAC1 inhibitor monotherapy group and the combination group via inhibiting AKT/FOXO3a signaling when compared to the control group. Moreover, the silencing of RAC1 inhibited AKT/FOXO3a signaling and cell glycolysis while the upregulation of RAC1 produced an opposite effect. In murine xenograft models, the tumor volume and the expression of glycolytic enzymes were significantly reduced in combination therapy when compared to each monotherapy group. Overall, our study demonstrates that targeting RAC1 with an inhibitor overcomes cisplatin resistance in ESCC by suppressing glycolytic enzymes, which provides a promising strategy for treatment of ESCC in clinical practice.

The interaction of lncRNA EZR-AS1 with SMYD3 maintains overexpression of EZR in ESCC cells

EZR, a member of the ezrin-radixin-moesin (ERM) family, is involved in multiple aspects of cell migration and cancer. SMYD3, a histone H3–lysine 4 (H3–K4)-specific methyltransferase, regulates EZR gene transcription, but the molecular mechanisms of epigenetic regulation remain ill-defined. Here, we show that antisense lncRNA EZR-AS1 was positively correlated with EZR expression in both human esophageal squamous cell carcinoma (ESCC) tissues and cell lines. Both in vivo and in vitro studies revealed that EZR-AS1 promoted cell migration through up-regulation of EZR expression. Mechanistically, antisense lncRNA EZR-AS1 formed a complex with RNA polymerase II to activate the transcription of EZR. Moreover, EZR-AS1 could recruit SMYD3 to a binding site, present in a GC-rich region downstream of the EZR promoter, causing the binding of SMYD3 and local enrichment of H3K4me3. Finally, the interaction of EZR-AS1 with SMYD3 further enhanced EZR transcription and expression. Our findings suggest that antisense lncRNA EZR-AS1, as a member of an RNA polymerase complex and through enhanced SMYD3-dependent H3K4 methylation, plays an important role in enhancing transcription of the EZR gene to promote the mobility and invasiveness of human cancer cells.

Lapachol is a novel ribosomal protein S6 kinase 2 inhibitor that suppresses growth and induces intrinsic apoptosis in esophageal squamous cell carcinoma cells

Lapachol is a 1,4-naphthoquinone that is isolated from the Bignoniaceae family. It has been reported to exert anti-inflammatory, antibacterial, and anticancer activities. However, the anticancer activity of lapachol and its molecular mechanisms against esophageal squamous cell carcinoma (ESCC) cells have not been fully investigated. Herein, we report that lapachol is a novel ribosomal protein S6 kinase 2 (RSK2) inhibitor that suppresses growth and induces intrinsic apoptosis in ESCC cells. We found that lapachol strongly attenuates downstream signaling molecules of RSK2 in ESCC cells and also directly inhibits RSK2 activity in vitro. The RSK protein is highly activated in ESCC cells and knockdown of RSK2 significantly suppresses anchorage-dependent and anchorage-independent growth of ESCC cells. Additionally, lapachol inhibits anchorage-dependent and anchorage-independent growth of ESCC cells, and the inhibition of cell growth by lapachol is dependent on the expression of RSK2. We also found that lapachol induces mitochondria-mediated cellular apoptosis by activating caspases-3, -7, and PARP, inducing the expression of cytochrome c and BAX by inhibiting downstream molecules of RSK2. Overall, lapachol is a potent RSK2 inhibitor that might be used for chemotherapy against ESCC.

Screening key lncRNAs for human rectal adenocarcinoma based on lncRNA-mRNA functional synergistic network

Background

Rectal adenocarcinoma (READ) is one of the deadliest malignancies, and the molecular mechanisms underlying the initiation and development of READ remain largely unknown. In this study, we aimed to find key long noncoding RNAs (lncRNAs) and mRNAs in READ by RNA sequencing.

Methods

RNA sequencing was performed to identify differentially expressed mRNAs (DEmRNAs) and lncRNAs (DElncRNAs) between READ and normal tissue. READ‐specific protein‐protein interaction (PPI), DElncRNA‐DEmRNA coexpression, and DElncRNA‐nearby DEmRNA interaction networks were constructed. DEmRNAs and DEmRNAs coexpressed with DElncRNAs were functionally annotated.

Results

A total of 2113 DEmRNAs and 150 DElncRNAs between READ and normal tissue were identified. The PPI network identified several hub proteins, including CDK1, AURKB, CDC6, FOXQ1, NUF2, and TOP2A. The DElncRNA‐DEmRNA coexpression and DElncRNA‐nearby DEmRNA interaction networks identified some hub lncRNAs, including CCAT1, LOC105374879, GAS5, and B3GALT5‐AS1. The colorectal cancer pathway, the intestinal immune network for IgA production and the p53 signaling pathway were three pathways significantly enriched in DEmRNAs and DEmRNAs coexpressed with DElncRNAs. MSH6 coexpressed with two DElncRNAs (LOC105374879 and CASC15) and BCL2 coexpressed with B3GALT5‐AS1 were significantly enriched in the colorectal cancer signaling pathway. TNFRSF17 coexpressed with B3GALT5‐AS1 was enriched in the intestinal immune network for IgA production. CCNB2 coexpressed with LOC105374879 was enriched in the p53 signaling pathway.

Conclusion

A total of four DEmRNAs (MSH6, BCL2, TNFRSF17, and CCNB2) and three DElncRNAs (LOC105374879, CASC15, and B3GALT5‐AS1) may be involved in the pathogenesis of READ; this data may contribute to understanding the mechanisms of READ and the development of therapeutic strategies for READ.

Bioinformatics identification of lncRNA biomarkers associated with the progression of esophageal squamous cell carcinoma

The poor outcome of patients with esophageal squamous cell carcinoma (ESCC) highlights the importance of the identification of novel effective prognostic biomarkers. Long non-coding RNAs (lncRNAs) serve regulatory roles in various types of cancer. The aim of the present study was to investigate the lncRNA expression profile in ESCC and to identify lncRNAs associated with the prognosis of ESCC by performing comprehensive bioinformatics analyses. The RNA-sequencing (Seq) expression dataset GSE53625 generated from ESCC samples was used as a training dataset. Additional RNA-Seq datasets relative to ESCC samples were downloaded from The Cancer Genome Atlas and used as a validation dataset. Data were screened using the limma package, and differentially expressed lncRNAs between early- and late-stage ESCC were identified. A random forest algorithm was used to select the optimal lncRNA biomarkers, which were then analyzed using the support vector machine (SVM) algorithm with R software. The identified lncRNA biomarkers were examined in the validation dataset by bidirectional hierarchical clustering and using an SVM classifier. Subsequently, univariate and multivariate Cox regression analyses were performed to analyze the potential ability lncRNAs to predict the survival rate of patients with ESCC. By examining the training group, 259 deregulated lncRNAs between early- and advanced-stage ESCC were identified. Further bioinformatics analyses identified a nine-lncRNA signature, including AC098973, AL133493, RP11-51M24, RP11-317N8, RP11-834C11, RP11-69C17, LINC00471, LINC01193 and RP1-124C. This nine-lncRNA signature was used to predict the tumor stage and patient survival rate with high reliability and accuracy in the training and validation datasets. Furthermore, these nine lncRNA biomarkers were primarily involved in regulating the cell cycle and DNA replication, and these processes were previously identified to be associated with the progression of ESCC. The identified nine-lncRNA signature was identified to be associated with the tumor stage, and could be used as predictor of the survival rate of patients with ESCC.

Clinical significance of LOXL4 expression and features of LOXL4-associated protein-protein interaction network in esophageal squamous cell carcinoma

Lysyl oxidase-like 4 (LOXL4), a member of the LOX family proteins, catalyzes oxidative deamination of lysine residues in collagen and elastin, which are responsible for maintaining extracellular matrix homeostasis. In this study, the mRNA expression of LOXL4 in seven esophageal squamous cell carcinoma (ESCC) cell lines and 15 ESCC pairs of clinical samples were examined. Furthermore, LOXL4 protein levels in the ESCC cell lines were determined using western blotting. With the use of immunofluorescence, LOXL4 was observed to be localized primarily in the cytoplasm, but was also present in the nucleus. In addition, the results indicated that the upregulated expression of LOXL4 was associated with poor survival in patients with ESCC even following curative resection (P = 0.010). Similar Kaplan-Meier estimator curves for proteins that interact with LOXL4, SUV39H1 (P = 0.014) and COL2A1 (P = 0.011), were plotted. The analyses based on the protein-protein interaction network depicted the expression of LOXL4 and its associated proteins as well as their functions, suggesting that LOXL4 and its associated proteins may serve a significant role in the development and progression of ESCC. In conclusion, the results of the present study suggest that LOXL4 is a potential biomarker for patients with ESCC, as well as SUV39H1 and COL2A1, and high expression levels of these genes are associated with poor prognosis in patients with ESCC.

Co-expression network modeling identifies key long non-coding RNA and mRNA modules in altering molecular phenotype to develop stress-induced depression in rats

Long non-coding RNAs (lncRNAs) have recently emerged as one of the critical epigenetic controllers, which participate in several biological functions by regulating gene transcription, mRNA splicing, protein interaction, etc. In a previous study, we reported that lncRNAs may play a role in developing depression pathophysiology. In the present study, we have examined how lncRNAs are co-expressed with gene transcripts and whether specific lncRNA/mRNA modules are associated with stress vulnerability or resiliency to develop depression. Differential regulation of lncRNAs and coding RNAs were determined in hippocampi of three group of rats comprising learned helplessness (LH, depression vulnerable), non-learned helplessness (NLH, depression resilient), and tested controls (TC) using a single-microarray-based platform. Weighted gene co-expression network analysis (WGCNA) was conducted to correlate the expression status of protein-coding transcripts with lncRNAs. The associated co-expression modules, hub genes, and biological functions were analyzed. We found signature co-expression networks as well as modules that underlie normal as well as aberrant response to stress. We also identified specific hub and driver genes associated with vulnerability and resilience to develop depression. Altogether, our study provides evidence that lncRNA associated complex trait-specific networks may play a crucial role in developing depression.

Long Read Single-Molecule Real-Time Sequencing Elucidates Transcriptome-Wide Heterogeneity and Complexity in Esophageal Squamous Cells

Esophageal squamous cell carcinoma is a leading cause of cancer death. Mapping the transcriptional landscapes such as isoforms, fusion transcripts, as well as long noncoding RNAs have played a central role to understand the regulating mechanism during malignant processes. However, canonical methods such as short-read RNA-seq are difficult to define the entire polyadenylated RNA molecules. Here, we combined single-molecule real-time sequencing with RNA-seq to generate high-quality long reads and to survey the transcriptional program in esophageal squamous cells. Compared with the recent annotations of human transcriptome (Ensembl 38 release 91), single-molecule real-time data identified many unannotated transcripts, novel isoforms of known genes and an expanding repository of long intergenic noncoding RNAs (lincRNAs). By integrating with annotation of lincRNA catalog, 1,521 esophageal-cancer-specific lincRNAs were defined from single-molecule real-time reads. Kyoto Encyclopedia of Genes and Genomes enrichment analysis indicated that these lincRNAs and their target genes are involved in a variety of cancer signaling pathways. Isoform usage analysis revealed the shifted alternative splicing patterns, which can be recaptured from clinical samples or supported by previous studies. Utilizing vigorous searching criteria, we also detected multiple transcript fusions, which are not documented in current gene fusion database or readily identified from RNA-seq reads. Two novel fusion transcripts were verified based on real-time PCR and Sanger sequencing. Overall, our long-read single-molecule sequencing largely expands current understanding of full-length transcriptome in esophageal cells and provides novel insights on the transcriptional diversity during oncogenic transformation.

lncRNA profile study reveals the mRNAs and lncRNAs associated with docetaxel resistance in breast cancer cells

Resistance to adjuvant systemic treatment, including taxanes (docetaxel and paclitaxel) is a major clinical problem for breast cancer patients. lncRNAs (long non-coding RNAs) are non-coding transcripts, which have recently emerged as important players in a variety of biological processes, including cancer development and chemotherapy resistance. However, the contribution of lncRNAs to docetaxel resistance in breast cancer and the relationship between lncRNAs and taxane-resistance genes are still unclear. Here, we performed comprehensive RNA sequencing and analyses on two docetaxel-resistant breast cancer cell lines (MCF7-RES and MDA-RES) and their docetaxel-sensitive parental cell lines. We identified protein coding genes and pathways that may contribute to docetaxel resistance. More importantly, we identified lncRNAs that were consistently up-regulated or down-regulated in both the MCF7-RES and MDA-RES cells. The co-expression network and location analyses pinpointed four overexpressed lncRNAs located within or near the ABCB1 (ATP-binding cassette subfamily B member 1) locus, which might up-regulate the expression of ABCB1. We also identified the lncRNA EPB41L4A-AS2 (EPB41L4A Antisense RNA 2) as a potential biomarker for docetaxel sensitivity. These findings have improved our understanding of the mechanisms underlying docetaxel resistance in breast cancer and have provided potential biomarkers to predict the response to docetaxel in breast cancer patients.

Screening and clinical significance of tumor markers in head and neck squamous cell carcinoma through bioinformatics analysis

In order to identify potential diagnostic and prognostic biomarkers, and treatment targets for head and neck squamous cell carcinoma (HNSCC), the present study obtained the gene expression profiles in HNSCC through public data mining, and core genes were identified using a series of bioinformatics analysis methods and databases. A total of nine hub genes (SPP1, ITGA6, TMPRSS11D, MMP1, LAMC2, FAT1, ACTA1, SERPINE1 and CEACAM1) were identified to be significantly correlated with HNSCC. Furthermore, overall survival analysis demonstrated that the expression values of hub genes were associated with overall survival in HNSCC. Furthermore, certain of the identified genes, including, TMPRSS11D, ACTA1 and CEACAM1, have not been thoroughly investigated in HNSCC previously. Taken together, the nine hub genes obtained by screening in the present study may serve as potential tumor markers and important prognostic indicators for HNSCC.

Gossypetin is a novel MKK3 and MKK6 inhibitor that suppresses esophageal cancer growth in vitro and in vivo

Gossypetin as a hexahydroxylated flavonoid found in many flowers and Hibiscus. It exerts various pharmacological activities, including antioxidant, antibacterial and anticancer activities. However, the anticancer capacity of gossypetin has not been fully elucidated. In this study, gossypetin was found to inhibit anchorage-dependent and -independent growth of esophageal cancer cells. To identify the molecular target(s) of gossypetin, various signaling protein kinases were screened and results indicate that gossypetin strongly attenuates the MKK3/6-p38 signaling pathway by directly inhibiting MKK3 and MKK6 protein kinase activity in vitro. Mechanistic investigations showed that arginine-61 in MKK6 is critical for binding with gossypetin. Additionally, the inhibition of cell growth by gossypetin is dependent on the expression of MKK3 and MKK6. Gossypetin caused G2 phase cell cycle arrest and induced intrinsic apoptosis by activating caspases 3 and 7 and increasing the expression of BAX and cytochrome c. Notably, gossypetin suppressed patient-derived esophageal xenograft tumor growth in an in vivo mouse model. Our findings suggest that gossypetin is an MKK3 and MKK6 inhibitor that could be useful for preventing or treating esophageal cancer.

Mutations, Differential Gene Expression, and Chimeric Transcripts in Esophageal Squamous Cell Carcinoma Show High Heterogeneity

Esophageal squamous cell carcinoma (ESCC) is a frequent and lethal neoplasia. As recent advances in targeted therapy have not improved ESCC prognosis, characterization of molecular alterations associated to this tumor is of foremost relevance. In this study, we analyze, for the first time, the complete genomic profile of ESCC by RNA-seq. TP53 was the most frequently mutated gene in the investigation and validation sets (78.6% and 67.4%, respectively). Differential expression analysis between tumor and nontumor adjacent mucosa showed 6698 differentially expressed genes, most of which were overexpressed (74%). Enrichment analysis identified overrepresentation of Wnt pathway, with overexpressed activators and underexpressed inactivators, suggesting activation of canonical and noncanonical Wnt signaling pathways. Higher WNT7B expression was associated with poor prognosis. Twenty-one gene fusions were identified in 50% of tumors, none of which involving the same genes in different patients; 71% of fusions involved syntenic genes. Comparisons with TCGA data showed co-amplification of seven gene pairs involved in fusions in the present study (~33%), suggesting that these rearrangements might have been driven by chromoanagenesis. In conclusion, genomic alterations in ESCC are highly heterogeneous, impacting negatively in target therapy development.

LncRNA co-expression network analysis reveals novel biomarkers for pancreatic cancer

High mortality and low survival rates for pancreatic ductal adenocarcinoma (PDAC) mainly result from the delay in diagnosis and treatment. Therefore, there is an urgent need to identify early PDAC biomarkers and new therapeutic targets. In this study, we applied a commonly used systems biology approach, the weighted gene co-expression network analysis (WGCNA), on lncRNA expression data. Eleven lncRNAs, namely A2M-AS1, DLEU2, LINC01133, LINC00675, MIR155HG, SLC25A25-AS1, LINC01857, LOC642852 (LINC00205), ITGB2-AS1, TSPOAP1-AS1 and PSMB8-AS1 have been identified and validated on an independent PDAC expression dataset. Furthermore, we characterized them by functional and pathway enrichment analysis and identified which lncRNAs showed differential expression, differential promoter methylation levels and copy number alterations between normal and PDAC samples. Finally, we also performed a survival analysis and identified A2M-AS1, LINC01133, LINC00205 and TSPOAP1-AS1 as prognostic biomarkers for PDAC. Interestingly, although only a few cancer-associated lncRNAs have been functionally characterized, LINC00675 and LINC01133 lncRNAs have already been demonstrated to be involved in PDAC development and progression. Therefore, our results provide new potential diagnostic/prognostic biomarkers and therapeutic targets for PDAC that deserve to be further investigated. Moreover, these lncRNAs may improve the understanding about molecular pathogenesis of PDAC.

SLC52A3 expression is activated by NF-κB p65/Rel-B and serves as a prognostic biomarker in esophageal cancer

The human riboflavin transporter-3 (encoded by SLC52A3) plays a prominent role in riboflavin absorption. Interestingly, abnormal expression patterns of SLC52A3 in multiple types of human cancers have been recently noted. However, the molecular mechanisms underlying its dysregulation remain unclear. In this study, we find that SLC52A3 has two transcript variants that differ in the transcriptional start site, and encode different proteins: SLC52A3a and SLC52A3b. Importantly, aberrant expressions of SLC52A3 are associated with stepwise development of esophageal squamous cell carcinoma (ESCC) as well as the survival rates of ESCC patients. Functionally, SLC52A3a, but not SLC52A3b, strongly promotes the proliferation and colony formation of ESCC cells. Furthermore, SLC52A3 5'-flanking regions contain NF-κB p65/Rel-B-binding sites, which are crucial for mediating SLC52A3 transcriptional activity in ESCC cells. Chromatin immunoprecipitation and electrophoretic mobility shift assay reveal that p65/Rel-B bind to 5'-flanking regions of SLC52A3. Accordingly, NF-κB signaling upregulates SLC52A3 transcription upon TNFα stimulation. Taken together, these results elucidate the mechanisms underlying SLC52A3 overexpression in ESCC. More importantly, our findings identify SLC52A3 as both a predictive and prognostic biomarker for this deadly cancer.