Spectrum of Mesenchymal-Epithelial Transition Aberrations and Potential Clinical Implications: Insights From Integrative Pancancer Analysis

A network pharmacology-based investigation of emodin against pancreatic adenocarcinoma

Pancreatic adenocarcinoma (PAAD) is one of the most common malignancies worldwide with an increasing incidence and poor outcome due to the lack of effective diagnostic and treatment methods. Emerging evidence implicates that emodin displays extensive spectrum anticancer properties. Differential expression genes in PAAD patients were analyzed by Gene Expression Profiling Interactive Analysis (GEPIA) website, and the targets of emodin were obtained via Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform. Subsequently, enrichment analyses were performed using R software. A protein-protein interaction (PPI) network was constructed by STRING database and Cytoscape software was used to identify the hub genes. Prognostic value and immune infiltration landscapes were explored through Kaplan-Meier plotter (KM plotter) website and the Single-Sample Gene Set Enrichment Analysis package of R. Finally, molecular docking was used to computationally verify the interaction of ligand and receptor proteins. A total of 9191 genes were significantly differentially expressed in PAAD patients and 34 potential targets of emodin were obtained. Intersections of the 2 groups were considered as potential targets of emodin against PAAD. Functional enrichment analyses illustrated that these potential targets were linked to numerous pathological processes. Hub genes identified through PPI networks were correlated with poor prognosis and infiltration level of different immune cells in PAAD patients. Perhaps emodin interacted with the key molecules and regulate the activity of them. We revealed the inherent mechanism of emodin against PAAD with the aid of network pharmacology, which provided reliable evidence and a novel guideline for clinical treatment.

Prognostic Outcome of Mesenchymal Transition Biomarkers in Correlation with EGFR Expression in Epithelial Ovarian Carcinoma Patients

BACKGROUND

CD44 is an epithelial-mesenchymal transition (EMT) surface receptor that regulates the interactivity between the cells and the extracellular matrix, thereby promoting cell migration. The epidermal growth factor receptor (EGFR) family is a trans-membrane kinase-related protein. It regulates cell adhesion proteins, which may promote cell proliferation and invasiveness. Mesenchymal epithelial transition (MET) is another EMT receptor that stimulates cell proliferation, invasion, survival, and angiogenesis. This study aimed to evaluate the prognostic impact of CD44, EGFR expressions, and MET gene amplification in epithelial ovarian cancer (EOC).

METHODS

This is a retrospective cohort study, including 85 cases of EOC. CD44 and EGFR expressions were evaluated in both epithelial and stromal cells by immunohistochemistry. Tumor cells also underwent a cytogenetic analysis using fluorescent in situ hybridization (FISH) to detect MET gene amplification.

RESULTS

High CD44 expression in tumors was significantly associated with serous subtypes (P=0.001), peritoneal deposits (P=0.002), and advanced stage (P=0.002). EGFR high tumor expression demonstrated a significant association with lymph node metastasis (P=0.038) and the advanced stage of EOC (P=0.016). Increased copy number of the MET gene was significantly associated with partial therapy response (P=0.030).  CD44 and EGFR tumor high expression was associated with poor overall survival (OS). In addition, MET gene gain in tumors was associated with a shorter OS (P=0.000).

CONCLUSION

EMT biomarkers (CD44 and MET) and EGFR expression in EOC are independent prognostic factors for OS. MET gene increase copy number was detected in cases of serous neoplasm and associated with poor survival and minimal therapy response.

Diverse Chromobox Family Members: Potential Prognostic Biomarkers and Therapeutic Targets in Head and Neck Squamous Cell Carcinoma

Background

The Chromobox (CBX) family members were involved in a variety of physiological and oncological processes through the regulation of the epigenetic modification of chromatin. However, the comprehensive analysis of the CBX family in head and neck squamous cell carcinoma (HNSC) is lacking.

Methods

In this work, we used multiple online databases and tools to investigate the roles of CBX family in aspects of gene expression, prognostic evaluation, genetic alteration, immune micro-environment of tumor, and status of methylation.

Results

The mRNA expression levels of CBX1, CBX3, and CBX5 were aberrantly increased in patients with HNSC, while CBX7 was aberrantly decreased. Higher expression of CBX7 was significantly associated with longer OS. Within the 5–11% of genetic alteration rate of CBXs, CBX3 ranked the highest and CBX5/7 ranked the lowest. SPRR1B, S100A7, CASP14, CDSN, LCE3D were the top 5 neighbor genes with the strongest association with CBXs in HNSC patients. Signaling pathways such as epidermal cell differentiation, cornification, and peptide cross-linking were demonstrated to have a strong association with CBX genes. The profiles of immune cell infiltration had high similarity for the group of HNSC patients stratified by expression of CBXs. The methylation levels of CBX1 and CBX5 significantly decreased, while that of CBX7 significantly increased in HNSC samples when compared with normal tissue.

Conclusion

In conclusion, the CBX family showed its valuation for further investigation in HNSC. Our research highlighted that CBX7 had the potential to be a novel diagnostic and prognostic biomarker for patients with HNSC.

Expression profile and prognostic values of Chromobox family members in human glioblastoma

Glioblastoma (GBM) is the most common and deadly malignant primary brain tumor. Chromobox (CBX) family proteins are essential components of the epigenetic regulatory complex and are involved in the occurrence and development of various cancers. However, the roles of CBX members in GBM is little known. In this analysis, we synthesized several mainstream bioinformatics databases to comprehensively explore the expression profiles, prognostic implications, genetic alterations, immune infiltration, and potential biological functions of the CBXs in GBM, and cell experiments were also conducted to investigate the role of CBX8 in GBM. We found that the elevated mRNA expression of CBX2/3/5/8 and reduced mRNA expression of CBX6/7 were found in GBM. The protein levels of CBX2/3/5/8 were elevated in GBM tissues, whereas the protein levels of CBX6/7 showed no significant difference. The upregulated expression of CBX2/3/8 was found to be both correlated with the tumor grade and recurrent status. The overexpression of CBX3/8 and underexpression of CBX6 mRNA were associated with the poor prognosis. These findings suggested that CBX3 and CBX8 might be useful diagnostic and prognostic biomarkers in GBM. Further cell experiment results supported that CBX8 promoted the proliferation of glioma cells. Moreover, a high genetic alteration rate of CBXs (37%) was found in GBM and to varying degrees. The expression of CBXs was significantly related to the immune cells infiltration. CBX7 methylation level was significantly increased in GBM tissues. Our results may provide novel ideas to find potential prognostic markers and new therapeutic targets among CBX family members in glioblastoma.

Integrated analysis of expression, prognostic value and immune infiltration of GSDMs in hepatocellular carcinoma

Six Gasdermins (GSDM) family members participate in various biological processes especially pyroptosis, as well as in the initiation and development of many types of cancer. However, the systematic analysis of the GSDM family in hepatocellular carcinoma (HCC) is lacking. In this study, several bioinformatics databases were recruited to analyze the roles of the GSDMs in differential expression, prognostic correlation, functional enrichment exploration, immune modulation, genetic alterations, and methylated modification in patients with HCC. Consequently, the mRNA expression of all the six GSDMs was accordantly increased in HCC, while only the protein expressions of GSDMB, GSDMD, and GSDME were apparently increased in HCC tissue. The expression of all the GSDMs (except GSDMA) was significantly higher in tumor stage 1-3 subgroups, compared with that in normal subgroups. Higher GSDME expression was significantly associated with shorter overall survival (OS) and disease specific survival (DSS) in patients with HCC. GSDMD had the highest genetic alteration rate among the GSDMs. The three signal pathways which were most likely related to GSDMs-associated molecules were the cell adhesion, growth regulation, and hormone metabolic process. The majority of GSDMs members were positively correlated with the infiltration of B cells, neutrophils, and dendritic cells, however negatively correlated with macrophage. All of the six GSDM members showed remarkably decreased methylation levels in HCC tissues. In conclusion, the GSDM family (especially GSDME) had the potential to become essential biomarkers to better improve the diagnosis and prognosis of HCC, as well as provided insight for the development of therapeutic targets.

Integrative pan-cancer analysis of MEK1 aberrations and the potential clinical implications

Alterations of mitogen-activated protein kinase kinase 1 (MEK1) are commonly associated with tumorigenesis, and MEK1 is thought to be a suitable targeted therapy for various cancers. However, abnormal MEK1 alterations and their relevant clinical implications are unknown. Our research comprehensively analyzed the MEK1 alteration spectrum and provided novel insight for targeted therapies. There were 7694 samples covering 32 types of cancer from The Cancer Genome Atlas (TCGA) database. They were used to conduct an integrative analysis of MEK1 expression, alterations, functional impacts and clinical significance. There was a dramatic difference in the alteration frequency and distribution and clinical implications in 32 types of cancer from the TCGA. Skin cutaneous melanoma (SKCM) has the most alterations and has therapeutic targets located in the protein kinase domain, and the growing expression of SKCM is positively related to patient prognosis. MEK1 expression in lung adenocarcinoma (LUAD), kidney renal papillary cell carcinoma (KIRP), esophageal carcinoma (ESCA) and liver hepatocellular carcinoma (LIHC) is decreased, which is associated with better prognosis, while MEK1 expression in thymoma (THYM), stomach adenocarcinoma (STAD), kidney renal clear cell carcinoma (KIRC), testicular germ cell tumors (TGCTs) and head and neck squamous cell carcinoma (HNSC) is increased, which is associated with better prognosis. Mesothelioma (MESO) has the second highest alterations but has no therapy targets. This study provided a great and detailed interpretation of MEK1 expression, alterations and clinical implications in 32 types of cancer and reminded us to fill the gap in MEK1 research from a new perspective.

Comprehensive Analysis of YTH Domain Family in Lung Adenocarcinoma: Expression Profile, Association with Prognostic Value, and Immune Infiltration

Background

All YTH domain family members are m⁶A reader proteins accounting for the methylation modulation involved in the process of tumorgenesis and tumor progression. However, the expression profiles and roles of the YTH domain family in lung adenocarcinoma (LUAD) remain to be further illustrated.

Methods

GEPIA2 and TNMplot databases were used to generate the expression profiles of the YTH family. Kaplan-Meier plotter database was employed to analysis the prognostic value of the YTH family. Coexpression profiles and genetic alterations analysis of the YTH family were undertaken using the cBioPortal database. YTH family protein-associated protein-protein interaction (PPI) network was identified by using STRING. Functional enrichment analysis was performed with the help of the WebGestalt database. The correlation analysis between the YTH family and immune cell infiltration in LUAD was administrated by using the TIMER2.0 database.

Results

mRNA expression of YTHDC1 and YTHDC2 was significantly lower in LUAD, whereas YTHDF1, YTHDF2, and YTHDF3 with apparently higher expression. YTHDF2 expression was observed to be the highest in the nonsmoker subgroup, and its expression gradually decreased with the increased severity of smoking habit. LUAD patients with low expression of YTHDC2, YTHDF1, and YTHDF2 were correlated with a better overall survival (OS) time. The YTHDF1 genetic alteration rate was 26%, which was the highest in the YTH family. The major cancer-associated functions of YTH family pointed in the direction of immunomodulation, especially antigen processing and presentation. Most of the YTH family members were significantly correlated with the infiltration of CD4+ T cells, CD8+ T cells, macrophages, and neutrophils, indicating the deep involvement of the YTH domain family in the immune cell infiltration in LUAD.

Conclusion

The molecular and expression profiles of the YTH family were dysregulated in LUAD. YTH family members (especially YTHDC2) were promising biomarkers and potential therapeutic targets that may bring benefit for the patients with LUAD.

Insights of fibroblast growth factor receptor 3 aberrations in pan-cancer and their roles in potential clinical treatment

Fibroblast growth factor receptor 3 (FGFR3) alters frequently across various cancer types and is a common therapeutic target in bladder urothelial carcinoma (BLCA) with FGFR3 variants. Although emerging evidence supports the role of FGFR3 in individual cancer types, no pan-cancer analysis is available. In this work, we used the open comprehensive datasets, covering a total of 10,953 patients with 10,967 samples across 32 TCGA cancer types, to identify the full alteration spectrum of FGFR3. FGFR3 abnormal expression, methylation patterns, alteration frequency, mutation location distribution, functional impact, and prognostic implications differed greatly from cancer to cancer. The overall alteration frequency of FGFR3 was relatively low in all cancers. Targetable mutations were mainly detected in BLCA, and S249C, Y373C, G370C, and R248C were hotspot mutations that could be targeted by an FDA approved erdafitinib. Genetic fusions were mainly observed in glioma, followed by BLCA. FGFR3-TACC3 was the most common fusion type which was proposed as novel therapeutic targets in glioma and was targetable with erdafitinib in BLCA. Lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC) were two lung cancer subtypes, FGFR3 fusion and hotspot mutation like S249C were observed more commonly in LUSC but not in LUAD. DNA methylation was correlated with the expression of FGFR3 and its downstream genes in some tumors. FGFG3 abnormal expression and alterations exhibited clinical correlations with patient prognosis in several tumors. This work exhibited the full alteration spectrum of FGFR3 and indicated several new clues for their application as potential therapeutic targets and prognostic indicators.