Signaling pathway network alterations in human ovarian cancers identified with quantitative mitochondrial proteomics

Clinically relevant immune subtypes based on alternative splicing landscape of immune-related genes for lung cancer advanced PPPM approach

Background

Alternative splicing (AS) occurs in the process of gene post-transcriptional process, which is very important for the correct synthesis and function of protein. The change of AS pattern may lead to the change of expression level or function of lung cancer-related genes, and then affect the occurrence and development of lung cancers. The specific AS pattern might be used as a biomarker for early warning and prognostic assessment of a cancer in the framework of predictive, preventive, and personalized medicine (PPPM; 3PM). AS events of immune-related genes (IRGs) were closely associated with tumor progression and immunotherapy. We hypothesize that IRG-AS events are significantly different in lung adenocarcinomas (LUADs) vs. controls or in lung squamous cell carcinomas (LUSCs) vs. controls. IRG-AS alteration profiling was identified to construct IRG-differentially expressed AS (IRG-DEAS) signature models. Study on the selective AS events of specific IRGs in lung cancer patients might be of great significance for further exploring the pathogenesis of lung cancer, realizing early detection and effective monitoring of lung cancer, finding new therapeutic targets, overcoming drug resistance, and developing more effective therapeutic strategies, and better used for the prediction, diagnosis, prevention, and personalized medicine of lung cancer.

Methods

The transcriptomic, clinical, and AS data of LUADs and LUSCs were downloaded from TCGA and its SpliceSeq databases. IRG-DEAS events were identified in LUAD and LUSC, followed by their functional characteristics, and overall survival (OS) analyses. OS-related IRG-DEAS prognostic models were constructed for LUAD and LUSC with Lasso regression, which were used to classify LUADs and LUSCs into low- and high-risk score groups. Furthermore, the immune cell distribution, immune-related scores, drug sensitivity, mutation status, and GSEA/GSVA status were analyzed between low- and high-risk score groups. Also, low- and high-immunity clusters and AS factor (SF)-OS-related-AS co-expression network and verification of cell function of CELF6 were analyzed in LUAD and LUSC.

Results

Comprehensive analysis of transcriptomic, clinical, and AS data of LUADs and LUSCs identified IRG-AS events in LUAD (n = 1607) and LUSC (n = 1656), including OS-related IRG-AS events in LUAD (n = 127) and LUSC (n = 105). A total of 66 IRG-DEAS events in LUAD and 89 IRG-DEAS events in LUSC were identified compared to controls. The overlapping analysis between IRG-DEASs and OS-related IRG-AS events revealed 14 OS-related IRG-DEAS events for LUAD and 16 OS-related IRG-DEAS events for LUSC, which were used to identify and optimize a 12-OS-related-IRG-DEAS signature prognostic model for LUAD and an 11-OS-related-IRG-DEAS signature prognostic model for LUSC. These two prognostic models effectively divided LUAD or LUSC samples into low- and high-risk score groups that were closely associated with OS, clinical characteristics, and tumor immune microenvironment, with significant gene sets and pathways enriched in the two groups. Moreover, weighted gene co-expression network (WGCNA) and nonnegative matrix factorization method (NMF) analyses identified four OS-relevant subtypes of LUAD and six OS-relevant subtypes of LUSC, and ssGSEA identified five immunity-relevant subtypes of LUAD and five immunity-relevant subtypes of LUSC. Interestingly, splicing factors-OS-related-AS network revealed hub molecule CELF6 was significantly related to the malignant phenotype in lung cancer cells.

Conclusions

This study established two reliable IRG-DEAS signature prognostic models and constructed interesting splicing factor-splicing event networks in LUAD and LUSC, which can be used to construct clinically relevant immune subtypes, patient stratification, prognostic prediction, and personalized medical services in the PPPM practice.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13167-024-00366-4.

Mitochondrial dysfunction route as a possible biomarker and therapy target for human cancer

Mitochondria are vital organelles found within living cells and have signalling, biosynthetic, and bioenergetic functions. Mitochondria play a crucial role in metabolic reprogramming, which is a characteristic of cancer cells and allows them to assure a steady supply of proteins, nucleotides, and lipids to enable rapid proliferation and development. Their dysregulated activities have been associated with the growth and metastasis of different kinds of human cancer, particularly ovarian carcinoma. In this review, we briefly demonstrated the modified mitochondrial function in cancer, including mutations in mtDNA, reactive oxygen species production, dynamics, apoptosis of cells, autophagy, and calcium excess to maintain cancer genesis, progression, and metastasis. Furthermore, the mitochondrial dysfunction pathway for some genomic, proteomic, and metabolomics modifications in ovarian cancer has been studied. Additionally, ovarian cancer has been linked to targeted therapies and biomarkers found through various alteration processes underlying mitochondrial dysfunction, notably targeting reactive oxygen species, metabolites, rewind metabolic pathways, and chemo-resistant ovarian carcinoma cells.

The Impact of Mitochondria in Ovarian Cancer Cell Metabolism, Proliferation, and Metastasis

Mitochondrial dysfunctions are significantly implicated in cancer initiation, progression, and metastasis, which have been shown for several cancers including ovarian cancer.An increase in mitochondrial dysfunction is also associated with drug resistance along with cancer progression, which in part is related to its specific microenvironment that is characterized by ascites, low glucose levels, and hypoxia that causes ovarian cancer cells to switch to mitochondrial respiration to enable their survival. Peritoneal ascitic fluid accumulation is a specific feature of ovarian cancer, and it is a major cause of its metastatic spread that also presents challenges for effective treatment. Among the treatment difficulties for ovarian cancer is the mutation rate and frequency of mtDNA in ovarian cancer tissue that can affect the efficiency of chemotherapeutic drugs. The varied and multiple mutations of different types enable metabolic reprogramming, cancer cell proliferation, and drug resistance.New specific information on mechanisms underlying several of the mitochondrial dysfunctions has led to proposing various mitochondrial determinants as targets for ovarian cancer therapy, which include targeting specific mitochondrial proteins and phosphoproteins as well as reactive oxygen species (ROS) that accumulate abnormally in cancer cells. Because of the genetically and histologically heterogeneous nature of the disease, combination therapy approaches will be necessary to combat the disease and achieve progress in effective treatment of ovarian cancer. This chapter will address (1) mitochondrial vulnerabilities underlying dysfunction and disease; (2) mitochondrial dysfunction in ovarian cancer; (3) present treatment difficulties for ovarian cancer and new potential treatment strategies to target ovarian cancer mitochondrial metabolism; and (4) biobehavioral factors influencing ovarian cancer development.

Ovarian Cancer: A Landscape of Mitochondria with Emphasis on Mitochondrial Dynamics

Ovarian cancer (OC) represents the main cause of death from gynecological malignancies in western countries. Altered cellular and mitochondrial metabolism are considered hallmarks in cancer disease. Several mitochondrial aspects have been found altered in OC, such as the oxidative phosphorylation system, oxidative stress and mitochondrial dynamics. Mitochondrial dynamics includes cristae remodeling, fusion, and fission processes forming a dynamic mitochondrial network. Alteration of mitochondrial dynamics is associated with metabolic change in tumour development and, in particular, the mitochondrial shaping proteins appear also to be responsible for the chemosensitivity and/or chemoresistance in OC. In this review a focus on the mitochondrial dynamics in OC cells is presented.

Identification of HIBCH as a Fatty Acid Metabolism-Related Biomarker in Aortic Valve Calcification Using Bioinformatics

Objective. To identify fatty acid metabolism-related biomarkers of aortic valve calcification (AVC) using bioinformatics and to research the role of immune cell infiltration for AVC. Methods. The AVC dataset was retrieved from the Gene Expression Omnibus database. R package is used for differential expression genes analysis and weighted gene coexpression analysis. The differentially coexpressed genes were identified by the Venn diagram, followed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of differentially coexpressed genes. Functions closely related to AVC were identified by GO and KEGG enrichment analyses of differentially coexpressed genes. Genes related to fatty acid metabolism were retrieved from the Molecular Signatures Database (MSigDB) database. After removing duplicate genes, least absolute shrinkage and selection operator (LASSO) regression analysis, support vector machine recursive feature elimination (SVM-RFE), and random forest were applied to recognize biomarkers related to fatty acid metabolism in AVC. The CIBERSORT tool was used to analyze infiltration of immune cells in normal and AVC samples. Correlations between biomarkers and immune cells were calculated. Finally, HIBCH-related pathway was predicted by single-gene gene set enrichment analysis (GSEA). Results. 2416 differentially expressed genes and one coexpression module were identified. A total of 1473 differentially coexpressed genes were acquired. GO and KEGG enrichment analyses demonstrated that differentially coexpressed genes were closely related to fatty acid metabolism. LASSO regression analysis, SVM-REF, and random forest revealed that 3-hydroxyisobutyryl-CoA hydrolase (HIBCH) was a biomarker of fatty acid metabolism-related genes in AVC. Significant high levels of memory B cells were found in AVC than normal samples, while activated natural killer (NK) cells were significantly low in AVC than normal samples. A significantly positive relevance was observed between HIBCH and activated NK cells, regulatory T cells, monocytes, naïve B cells, activated dendritic cells, resting memory CD4 T cells, resting NK cells, and CD8 T cells. A significantly negative relevance was observed between HIBCH and activated memory CD4 T cells, memory B cells, neutrophils, gamma delta T cells, M0 macrophages, and plasma cells. The single-gene GSEA results suggest that HIBCH may work through the inhibition of multiple immune-related pathways. Conclusion. HIBCH is closely relevant to immune cell infiltration in AVC and could be applied as a diagnostic marker for AVC.

Molecular Analysis of Short- versus Long-Term Survivors of High-Grade Serous Ovarian Carcinoma

Despite having similar histologic features, patients with high-grade serous ovarian carcinoma (HGSC) often experience highly variable outcomes. The underlying determinants for long-term survival (LTS, ≥10 years) versus short-term survival (STS, <3 years) are largely unknown. The present study sought to identify molecular predictors of LTS for women with HGSC. A cohort of 24 frozen HGSC samples was collected (12 LTS and 12 STS) and analyzed at DNA, RNA, and protein levels. OVCAR5 and OVCAR8 cell lines were used for in vitro validation studies. For in vivo studies, we injected OVCAR8 cells into the peritoneal cavity of female athymic nude mice. From RNAseq analysis, 11 genes were found to be differentially expressed between the STS and LTS groups (fold change > 2; false discovery rate < 0.01). In the subsequent validation cohort, transmembrane protein 62 (TMEM62) was found to be related to LTS. CIBERSORT analysis showed that T cells (follicular helper) were found at higher levels in tumors from LTS than STS groups. In vitro data using OVCAR5 and OVCAR8 cells showed decreased proliferation with TMEM62 overexpression and positive correlation with a longevity-regulating pathway (KEGG HSA04213) at the RNA level. In vivo analysis using the OVCAR8-TMEM62-TetON model showed decreased tumor burden in mice with high- vs. low-expressing TMEM62 tumors. Our results demonstrate that restoring TMEM62 may be a novel approach for treatment of HGSC. These findings may have implications for biomarker and intervention strategies to help improve patient outcomes

KAZN as a diagnostic marker in ovarian cancer: a comprehensive analysis based on microarray, mRNA-sequencing, and methylation data

Background

Ovarian cancer (OC) is among the deadliest malignancies in women and the lack of appropriate markers for early diagnosis leads to poor prognosis in most cases. Previous studies have shown that KAZN is involved in multiple biological processes during development, such as cell proliferation, differentiation, and apoptosis, so defects or aberrant expression of KAZN might cause queer cell behaviors such as malignancy. Here we evaluated the KAZN expression and methylation levels for possible use as an early diagnosis marker for OC.

Methods

We used data from Gene Expression Omnibus (GEO) microarrays, The Cancer Genome Atlas (TCGA), and Clinical Proteomic Tumor Analysis Consortium (CPTAC) to investigate the correlations between KAZN expression and clinical characteristics of OC by comparing methylation levels of normal and OC samples. The relationships among differentially methylated sites in the KAZN gene, corresponding KAZN mRNA expression levels and prognosis were analyzed.

Results

KAZN was up-regulated in ovarian epithelial tumors and the expression of KAZN was correlated with the patients’ survival time. KAZN CpG site cg17657618 was positively correlated with the expression of mRNA and the methylation levels were significantly differential between the group of stage “I and II” and the group of stage “III and IV”. This study also presents a new method to classify tumor and normal tissue in OC using DNA methylation pattern in the KAZN gene body region.

Conclusions

KAZN was involved in ovarian cancer pathogenesis. Our results demonstrate a new direction for ovarian cancer research and provide a potential diagnostic biomarker as well as a novel therapeutic target for clinical application.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-022-09747-2.

The Value of H2BC12 for Predicting Poor Survival Outcomes in Patients With WHO Grade II and III Gliomas

Purpose: Glioma is a common primary malignant brain tumor. Grade II (GII) gliomas are prone to develop into anaplastic grade III (GIII) gliomas, which indicate a higher malignancy and poorer survival outcome. This study aimed to satisfy the increasing demand for novel sensitive biomarkers and potential therapeutic targets in the treatment of GII and GIII gliomas.

Methods: A TCGA dataset was used to investigate the expression of H2BC12 mRNA in GII and GIII gliomas and its relation to clinical pathologic characteristics. Glioma tissues were collected to verify results from the TCGA dataset, and H2BC12 mRNA was detected by RT-qPCR. ROC analysis was employed to evaluate the classification power for GII and GIII. The significance of H2BC12 mRNA GII and GIII gliomas was also investigated. In addition, H2BC12 expression-related pathways were enriched by gene set enrichment analysis (GSEA). DNA methylation level and mutation of H2BC12 were analyzed by the UALCAN and CBioPortal databases, respectively.

Results: Based on the sample data from multiple databases and RT-qPCR, higher expression of H2BC12 mRNA was found in GII and GIII glioma tissue compared to normal tissue, which was consistent with a trend with our clinical specimen. H2BC12 mRNA had a better power in distinguishing between GII and GIII and yielded an AUC of 0.706 with a sensitivity of 76.9% and specificity of 81.8%. Meanwhile, high H2BC12 levels were associated with IDH status, 1p/19q codeletion, primary therapy outcome, and the histological type of gliomas. Moreover, the overall survival (OS), disease-specific survival (DSS), and progress-free interval (PFI) of GII glioma patients with higher levels of H2BC12 were shorter than those of patients with lower levels as well as GIII patients. In the multivariate analysis, a high H2BC12 level was an independent predictor for poor survival outcomes of gliomas. The Wnt or PI3K-AKT signaling pathways, DNA repair, cellular senescence, and DNA double-strand break repair were differentially activated in phenotypes that were positively associated with H2BC12. H2BC12 DNA methylation was high in TP53 nonmutant patients, and no H2BC12 mutation was observed in gliomas patients.

Conclusion: H2BC12 is a promising biomarker for the diagnosis and prognosis of patients with WHO grade II and III gliomas.

Mitochondrial Dysfunction Pathway Alterations Offer Potential Biomarkers and Therapeutic Targets for Ovarian Cancer

The mitochondrion is a very versatile organelle that participates in some important cancer-associated biological processes, including energy metabolism, oxidative stress, mitochondrial DNA (mtDNA) mutation, cell apoptosis, mitochondria-nuclear communication, dynamics, autophagy, calcium overload, immunity, and drug resistance in ovarian cancer. Multiomics studies have found that mitochondrial dysfunction, oxidative stress, and apoptosis signaling pathways act in human ovarian cancer, which demonstrates that mitochondria play critical roles in ovarian cancer. Many molecular targeted drugs have been developed against mitochondrial dysfunction pathways in ovarian cancer, including olive leaf extract, nilotinib, salinomycin, Sambucus nigra agglutinin, tigecycline, and eupatilin. This review article focuses on the underlying biological roles of mitochondrial dysfunction in ovarian cancer progression based on omics data, potential molecular relationship between mitochondrial dysfunction and oxidative stress, and future perspectives of promising biomarkers and therapeutic targets based on the mitochondrial dysfunction pathway for ovarian cancer.

Comparative microsomal proteomics of a model lung cancer cell line NCI-H23 reveals distinct differences between molecular profiles of 3D and 2D cultured cells

Lung cancer is the leading cause of cancer-related deaths in the USA and worldwide. Yet, about 95% of new drug candidates validated in preclinical phase eventually fail in clinical trials. Such a high attrition rate is attributed mostly to the inability of conventional two-dimensionally (2D) cultured cancer cells to mimic native three-dimensional (3D) growth of malignant cells in human tumors. To ascertain phenotypical differences between these two distinct culture conditions, we carried out a comparative proteomic analysis of a membrane fraction obtained from 3D- and 2D-cultured NSCLC model cell line NCI-H23. This analysis revealed a map of 1,166 (24%) protein species regulated in culture dependent manner, including differential regulation of a subset of cell surface-based CD molecules. We confirmed exclusive expression of CD99, CD146 and CD239 in 3D culture. Furthermore, label-free quantitation, targeting KRas proteoform-specific peptides, revealed upregulation of both wild type and monoallelic KRas4BG12C mutant at the surface of 3D cultured cells. In order to reduce the high attrition rate of new drug candidates, the results of this study strongly suggests exploiting base-line molecular profiling of a large number of patient-derived NSCLC cell lines grown in 2D and 3D culture, prior to actual drug candidate testing.

HSP60 Regulates Lipid Metabolism in Human Ovarian Cancer

Accumulating evidence demonstrates that cancer is an oxidative stress-related disease, and oxidative stress is closely linked with heat shock proteins (HSPs). Lipid oxidative stress is derived from lipid metabolism dysregulation that is closely associated with the development and progression of malignancies. This study sought to investigate regulatory roles of HSPs in fatty acid metabolism abnormality in ovarian cancer. Pathway network analysis of 5115 mitochondrial expressed proteins in ovarian cancer revealed various lipid metabolism pathway alterations, including fatty acid degradation, fatty acid metabolism, butanoate metabolism, and propanoate metabolism. HSP60 regulated the expressions of lipid metabolism proteins in these lipid metabolism pathways, including ADH5, ECHS1, EHHADH, HIBCH, SREBP1, ACC1, and ALDH2. Further, interfering HSP60 expression inhibited migration, proliferation, and cell cycle and induced apoptosis of ovarian cancer cells in vitro. In addition, mitochondrial phosphoproteomics and immunoprecipitation-western blot experiments identified and confirmed that phosphorylation occurred at residue Ser70 in protein HSP60, which might regulate protein folding of ALDH2 and ACADS in ovarian cancers. These findings clearly demonstrated that lipid metabolism abnormality occurred in oxidative stress-related ovarian cancer and that HSP60 and its phosphorylation might regulate this lipid metabolism abnormality in ovarian cancer. It opens a novel vision in the lipid metabolism reprogramming in human ovarian cancer.

The Anterior GRadient (AGR) family proteins in epithelial ovarian cancer

Epithelial ovarian cancer (EOC) is the most common gynecologic disorder. Even with the recent progresses made towards the use of new therapeutics, it still represents the most lethal gynecologic malignancy in women from developed countries.

The discovery of the anterior gradient proteins AGR2 and AGR3, which are highly related members belonging to the protein disulfide isomerase (PDI) family, attracted researchers’ attention due to their putative involvement in adenocarcinoma development. This review compiles the current knowledge on the role of the AGR family and the expression of its members in EOC and discusses the potential clinical relevance of AGR2 and AGR3 for EOC diagnosis, prognosis, and therapeutics.

A better understanding of the role of the AGR family may thus provide new handling avenues for EOC patients.

Protective Effects of Flavonoids Against Mitochondriopathies and Associated Pathologies: Focus on the Predictive Approach and Personalized Prevention

Multi-factorial mitochondrial damage exhibits a “vicious circle” that leads to a progression of mitochondrial dysfunction and multi-organ adverse effects. Mitochondrial impairments (mitochondriopathies) are associated with severe pathologies including but not restricted to cancers, cardiovascular diseases, and neurodegeneration. However, the type and level of cascading pathologies are highly individual. Consequently, patient stratification, risk assessment, and mitigating measures are instrumental for cost-effective individualized protection. Therefore, the paradigm shift from reactive to predictive, preventive, and personalized medicine (3PM) is unavoidable in advanced healthcare. Flavonoids demonstrate evident antioxidant and scavenging activity are of great therapeutic utility against mitochondrial damage and cascading pathologies. In the context of 3PM, this review focuses on preclinical and clinical research data evaluating the efficacy of flavonoids as a potent protector against mitochondriopathies and associated pathologies.

Development of a Prognostic Five-Gene Signature for Diffuse Lower-Grade Glioma Patients

Background: Diffuse lower-grade gliomas (LGGs) are infiltrative and heterogeneous neoplasms. Gene signature including multiple protein-coding genes (PCGs) is widely used as a tumor marker. This study aimed to construct a multi-PCG signature to predict survival for LGG patients.

Methods: LGG data including PCG expression profiles and clinical information were downloaded from The Cancer Genome Atlas (TCGA) and the Chinese Glioma Genome Atlas (CGGA). Survival analysis, receiver operating characteristic (ROC) analysis, and random survival forest algorithm (RSFVH) were used to identify the prognostic PCG signature.

Results: From the training (n = 524) and test (n = 431) datasets, a five-PCG signature which can classify LGG patients into low- or high-risk group with a significantly different overall survival (log rank P < 0.001) was screened out and validated. In terms of prognosis predictive performance, the five-PCG signature is stronger than other clinical variables and IDH mutation status. Moreover, the five-PCG signature could further divide radiotherapy patients into two different risk groups. GO and KEGG analysis found that PCGs in the prognostic five-PCG signature were mainly enriched in cell cycle, apoptosis, DNA replication pathways.

Conclusions: The new five-PCG signature is a reliable prognostic marker for LGG patients and has a good prospect in clinical application.

The mitochondrial landscape of ovarian cancer: emerging insights

Ovarian cancer (OC) is known to be the most lethal cancer in women worldwide, and its etiology is poorly understood. Recent studies show that mitochondrial DNA (mtDNA) content as well as mtDNA and nuclear genes encoding mitochondrial proteins influence OC risk. This review presents an overview of role of mitochondrial genetics in influencing OC development and discusses the contribution of mitochondrial proteome in OC development, progression and therapy. A role of mitochondrial genetics in racial disparity is also highlighted. In-depth understanding of role of mitochondria in OC will help develop strategies toward prevention and treatment and improving overall survival in women with OC.

Modelling the Functions of Polo-Like Kinases in Mice and Their Applications as Cancer Targets with a Special Focus on Ovarian Cancer

Polo-like kinases (PLKs) belong to a five-membered family of highly conserved serine/threonine kinases (PLK1-5) that play differentiated and essential roles as key mitotic kinases and cell cycle regulators and with this in proliferation and cellular growth. Besides, evidence is accumulating for complex and vital non-mitotic functions of PLKs. Dysregulation of PLKs is widely associated with tumorigenesis and by this, PLKs have gained increasing significance as attractive targets in cancer with diagnostic, prognostic and therapeutic potential. PLK1 has proved to have strong clinical relevance as it was found to be over-expressed in different cancer types and linked to poor patient prognosis. Targeting the diverse functions of PLKs (tumor suppressor, oncogenic) are currently at the center of numerous investigations in particular with the inhibition of PLK1 and PLK4, respectively in multiple cancer trials. Functions of PLKs and the effects of their inhibition have been extensively studied in cancer cell culture models but information is rare on how these drugs affect benign tissues and organs. As a step further towards clinical application as cancer targets, mouse models therefore play a central role. Modelling PLK function in animal models, e.g., by gene disruption or by treatment with small molecule PLK inhibitors offers promising possibilities to unveil the biological significance of PLKs in cancer maintenance and progression and give important information on PLKs' applicability as cancer targets. In this review we aim at summarizing the approaches of modelling PLK function in mice so far with a special glimpse on the significance of PLKs in ovarian cancer and of orthotopic cancer models used in this fatal malignancy.

SILAC quantitative proteomics analysis of ivermectin-related proteomic profiling and molecular network alterations in human ovarian cancer cells

The antiparasitic agent ivermectin offers more promises to treat a diverse range of diseases. However, a comprehensive proteomic analysis of ivermectin-treated ovarian cancer (OC) cells has not been performed. This study sought to identify ivermectin-related proteomic profiling and molecular network alterations in human OC cells. Stable isotope labeling with amino acids in cell culture (SILAC)-based quantitative proteomics was used to study the human OC TOV-21G cells. After TOV-21G cells underwent 10 passages in SILAC-labeled growth media, TOV-21G cells were treated with 10 ml of 20 μmol/L ivermectin in cell growing medium for 24 h. The SILAC-labeled proteins were digested with trypsin; tryptic peptides were identified with mass spectrometry (MS). Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis was used to mine signaling pathway alterations with ivermectin-related proteins in TOV-21G cells. Gene ontology (GO) analysis was used to explore biological functions of ivermectin-related proteins, including biological processes (BPs), cellular components (CCs), and molecular functions (MFs). The protein-protein interaction network was analyzed with molecular complex detection (MCODE) to identify hub modules. In total, 4,447 proteins were identified in ivermectin-treated TOV-21G cells. KEGG analysis revealed 89 statistically significant signaling pathways. Interestingly, the clustering analysis of these pathways showed that ivermectin was involved in various cancer pathogenesis processes, including modulation of replication, RNA metabolism, and translational machinery. GO analysis revealed 69 statistically significant CCs, 87 MFs, and 62 BPs. Furthermore, MCODE analysis identified five hub modules, including 147 hub molecules. Those hub modules involved ribosomal proteins, RNA-binding proteins, cell-cycle progression-related proteins, proteasome subunits, and minichromosome maintenance proteins. These findings demonstrate that SILAC quantitative proteomics is an effective method to analyze ivermectin-treated cells, provide the first ivermectin-related proteomic profiling and molecular network alterations in human OC cells, and provide deeper insights into molecular mechanisms and functions of ivermectin to inhibit OC cells.

Quantitative proteomics revealed energy metabolism pathway alterations in human epithelial ovarian carcinoma and their regulation by the antiparasite drug ivermectin: data interpretation in the context of 3P medicine

Objective

Energy metabolism abnormality is the hallmark in epithelial ovarian carcinoma (EOC). This study aimed to investigate energy metabolism pathway alterations and their regulation by the antiparasite drug ivermectin in EOC for the discovery of energy metabolism pathway-based molecular biomarker pattern and therapeutic targets in the context of predictive, preventive, and personalized medicine (PPPM) in EOC.

Methods

iTRAQ-based quantitative proteomics was used to identify mitochondrial differentially expressed proteins (mtDEPs) between human EOC and control mitochondrial samples isolated from 8 EOC and 11 control ovary tissues from gynecologic surgery of Chinese patients, respectively. Stable isotope labeling with amino acids in cell culture (SILAC)-based quantitative proteomics was used to analyze the protein expressions of energy metabolic pathways in EOC cells treated with and without ivermectin. Cell proliferation, cell cycle, apoptosis, and important molecules in energy metabolism pathway were examined before and after ivermectin treatment of different EOC cells.

Results

In total, 1198 mtDEPs were identified, and various mtDEPs were related to energy metabolism changes in EOC, with an interesting result that EOC tissues had enhanced abilities in oxidative phosphorylation (OXPHOS), Kreb's cycle, and aerobic glycolysis, for ATP generation, with experiment-confirmed upregulations of UQCRH in OXPHOS; IDH2, CS, and OGDHL in Kreb's cycle; and PKM2 in glycolysis pathways. Importantly, PDHB that links glycolysis with Kreb's cycle was upregulated in EOC. SILAC-based quantitative proteomics found that the protein expression levels of energy metabolic pathways were regulated by ivermectin in EOC cells. Furthermore, ivermectin demonstrated its strong abilities to inhibit proliferation and cell cycle and promote apoptosis in EOC cells, through molecular networks to target PFKP in glycolysis; IDH2 and IDH3B in Kreb's cycle; ND2, ND5, CYTB, and UQCRH in OXPHOS; and MCT1 and MCT4 in lactate shuttle to inhibit EOC growth.

Conclusions

Our findings revealed that the Warburg and reverse Warburg effects coexisted in human ovarian cancer tissues, provided the first multiomics-based molecular alteration spectrum of ovarian cancer energy metabolism pathways (aerobic glycolysis, Kreb's cycle, oxidative phosphorylation, and lactate shuttle), and demonstrated that the antiparasite drug ivermectin effectively regulated these changed molecules in energy metabolism pathways and had strong capability to inhibit cell proliferation and cell cycle progression and promote cell apoptosis in ovarian cancer cells. The observed molecular changes in energy metabolism pathways bring benefits for an in-depth understanding of the molecular mechanisms of energy metabolism heterogeneity and the discovery of effective biomarkers for individualized patient stratification and predictive/prognostic assessment and therapeutic targets/drugs for personalized therapy of ovarian cancer patients.

Uncovering Mitochondrial Determinants of Racial Disparities in Ovarian Cancer

Ovarian cancer (OC) incidence and mortality rates differ between racial groups. Mitochondrial genetic factors are now emerging as determinants of racial disparities in OC. A comprehensive understanding of the role of mitochondria in OC health disparities will help in developing novel therapeutic strategies targeting mitochondria to reduce or eliminate racial health disparities.

RSPO3 is a marker candidate for predicting tumor aggressiveness in ovarian cancer

Background

Ovarian cancer, a highly aggressive and heterogeneous gynecological malignancy that has long been difficult for physicians to identify and treat, requires more effective and precise molecular targets. R-spondin 3 (RSPO3) is a secreted protein that plays a tumorigenic role in several human cancers. However, the functional contribution and prognostic role of RSPO3 in ovarian cancer remain unclear.

Methods

RSPO3 expression in ovarian cancer tissues was assessed using western blotting, quantitative real-time polymerase chain reaction (qRT-PCR), and immunohistochemistry, and its relationships to clinicopathological parameters were investigated using the data of 179 ovarian cancer patients. RSPO3’s biological function was evaluated using Cell Counting Kit-8, colony formation, wound healing, and Matrigel transwell assay in RSPO3-knockdown and RSPO3-overexpression ovarian cancer cell lines SKOV3 and OVCAR3. The possible biological processes associated with RSPO3 were identified using functional enrichment analysis based on the transcriptome sequencing data from The Cancer Genome Atlas (TCGA) ovarian cancer cohort and our experimental cells, and further verified using western blotting and immunofluorescence in the ovarian cancer cell model.

Results

The RSPO3 mRNA and protein levels were both upregulated in ovarian cancer tissues. High RSPO3 expression was correlated with lymphovascular space invasion (LVSI), lymph node metastasis, distant metastasis, and advanced tumor stage. Survival analysis showed that RSPO3 is an independent prognostic marker in ovarian cancer. Moreover, in vitro RSPO3 knockdown significantly inhibited the invasion ability of ovarian cancer cells, while overexpression significantly promoted it. Using transcriptome sequencing and pathway validation experiments, we demonstrated for the first time that RSPO3 promotes ovarian cancer invasiveness through activation of the PI3K/AKT pathway and modulation of epithelial-mesenchymal transition (EMT), while the common Wnt/β-catenin signaling pathway was not involved.

Conclusions

RSPO3 plays a definite oncogenic role and promotes tumor aggressiveness in ovarian cancer, which may serve as a potential prognostic marker and therapeutic target for this disease.

Rare genetic variants suggest dysregulation of signaling pathways in low- and high-risk patients developing severe ovarian hyperstimulation syndrome

PURPOSE

To investigate if rare gene variants in women with severe ovarian hyperstimulation syndrome (OHSS) provide clues to the mechanisms involved in the syndrome.

METHODS

Among participants in a prospective randomized study (Toftager et al. 2016), six women with predicted low and six women with predicted high risk of OHSS developing severe OHSS (grades 4 and 5, Golan classification) were selected. In the same cohort, six plus six matched controls developing no signs of OHSS (Golan grade 0) were selected. Whole-exome sequencing was performed. Analysis using a predefined in silico OHSS gene panel, variant filtering, and pathway analyses was done.

RESULTS

We found no convincing monogenetic association with the development of OHSS using the in silico gene panel. Pathway analysis of OHSS variant lists showed substantial overlap in highly enriched top pathways (p value range p < 0.0001 and p > 9.8E-17) between the low- and high-risk group developing severe OHSS, i.e., "the integrin-linked kinase (ILK) signaling pathway" and the "axonal guidance signaling pathway," both being connected to vasoactive endothelial growth factor (VEGF) and endothelial function.

CONCLUSION

Rare variants in OHSS cases with two distinct risk profiles enrich the same signaling pathways linked to VEGF and endothelial function. Clarification of the mechanism as well as potentially defining genetic predisposition of the high vascular permeability is important for future targeted treatment and prevention of OHSS; the potential roles of ILK signaling and the axonal guidance signaling need to be validated by functional studies.

MASS SPECTROMETRY-BASED MITOCHONDRIAL PROTEOMICS IN HUMAN OVARIAN CANCERS

The prominent characteristics of mitochondria are highly dynamic and regulatory, which have crucial roles in cell metabolism, biosynthetic, senescence, apoptosis, and signaling pathways. Mitochondrial dysfunction might lead to multiple serious diseases, including cancer. Therefore, identification of mitochondrial proteins in cancer could provide a global view of tumorigenesis and progression. Mass spectrometry-based quantitative mitochondrial proteomics fulfils this task by enabling systems-wide, accurate, and quantitative analysis of mitochondrial protein abundance, and mitochondrial protein posttranslational modifications (PTMs). Multiple quantitative proteomics techniques, including isotope-coded affinity tag, stable isotope labeling with amino acids in cell culture, isobaric tags for relative and absolute quantification, tandem mass tags, and label-free quantification, in combination with different PTM-peptide enrichment methods such as TiO₂ enrichment of tryptic phosphopeptides and antibody enrichment of other PTM-peptides, increase flexibility for researchers to study mitochondrial proteomes. This article reviews isolation and purification of mitochondria, quantitative mitochondrial proteomics, quantitative mitochondrial phosphoproteomics, mitochondrial protein-involved signaling pathway networks, mitochondrial phosphoprotein-involved signaling pathway networks, integration of mitochondrial proteomic and phosphoproteomic data with whole tissue proteomic and transcriptomic data and clinical information in ovarian cancers (OC) to in-depth understand its molecular mechanisms, and discover effective mitochondrial biomarkers and therapeutic targets for predictive, preventive, and personalized treatment of OC. This proof-of-principle model about OC mitochondrial proteomics is easily implementable to other cancer types. © 2020 John Wiley & Sons Ltd. Mass Spec Rev.

Liquid Biopsy is Instrumental for 3PM Dimensional Solutions in Cancer Management

One in every four deaths is due to cancer in Europe. In view of its increasing incidence, cancer became the leading cause of death and disease burden in Denmark, France, the Netherlands, and the UK. Without essential improvements in cancer prevention, an additional 775,000 cases of annual incidence have been prognosed until 2040. Between 1995 and 2018, the direct costs of cancer doubled from EUR 52 billion to EUR 103 billion in Europe, and per capita health spending on cancer increased by 86% from EUR 105 to EUR 195 in general, whereby Austria, Germany, Switzerland, Benelux, and France spend the most on cancer care compared to other European countries. In view of the consequent severe socio-economic burden on society, the paradigm change from a reactive to a predictive, preventive, and personalized medical approach in the overall cancer management is essential. Concepts of predictive, preventive, and personalized medicine (3PM) demonstrate a great potential to revise the above presented trends and to implement cost-effective healthcare that benefits the patient and society as a whole. At any stage, application of early and predictive diagnostics, targeted prevention, and personalization of medical services are basic pillars making 3PM particularly attractive for the patients as well as ethical and cost-effective healthcare. Optimal 3PM approach requires novel instruments such as well-designed liquid biopsy application. This review article highlights current achievements and details liquid biopsy approaches specifically in cancer management. 3PM-relevant expert recommendations are provided.

Altered Organelle Calcium Transport in Ovarian Physiology and Cancer

Calcium levels have a huge impact on the physiology of the female reproductive system, in particular, of the ovaries. Cytosolic calcium levels are influenced by regulatory proteins (i.e., ion channels and pumps) localized in the plasmalemma and/or in the endomembranes of membrane-bound organelles. Imbalances between plasma membrane and organelle-based mechanisms for calcium regulation in different ovarian cell subtypes are contributing to ovarian pathologies, including ovarian cancer. In this review, we focused our attention on altered calcium transport and its role as a contributor to tumor progression in ovarian cancer. The most important proteins described as contributing to ovarian cancer progression are inositol trisphosphate receptors, ryanodine receptors, transient receptor potential channels, calcium ATPases, hormone receptors, G-protein-coupled receptors, and/or mitochondrial calcium uniporters. The involvement of mitochondrial and/or endoplasmic reticulum calcium imbalance in the development of resistance to chemotherapeutic drugs in ovarian cancer is also discussed, since Ca²⁺ channels and/or pumps are nowadays regarded as potential therapeutic targets and are even correlated with prognosis.

Identification of pathology-specific regulators of m6A RNA modification to optimize lung cancer management in the context of predictive, preventive, and personalized medicine

Relevance

Lung cancer is the most common malignant tumor with high morbidity (11.6% of the total diagnosed cancer cases) and mortality (18.4% of the total cancer deaths), and its 5-year survival rate is very low (20%). Clarification of any molecular events and the discovery of effective biomarkers will offer increasing promise for lung canner management. N⁶-methyladenosine (m⁶A) modification is one of the important RNA modifications that are closely associated with lung cancer, and are tightly regulated by m⁶A regulators. Elucidation of pathology-specific m⁶A regulators will directly contribute to lung cancer medical services in the context of predictive, preventive, and personalized medicine (PPPM).

Purpose

To investigate pathology-specific regulators of m⁶A RNA modifications in lung cancer and further inspect the m⁶A regulator gene signature as useful tools for PPPM in lung cancers.

Methods

The gene expression data of 19 m⁶A regulators (m⁶A-methyltransferases-ZC3H13, KIAA1429, RBM15/15B, WTAP, and METTL3/14; demethylases-FTO and ALKBH5; and m⁶A-binding proteins-HNRNPC, YTHDF1/2/3, YTHDC1/2, IGF2BP1/2/3, and HNRNPA2B1) and clinical data of 1013 lung cancer patients [511 lung adenocarcinoma (LUAD) and 502 lung squamous carcinoma (LUSC)] and 109 controls (Con) were obtained from the TCGA database. Quantitative real-time PCR (qRT-PCR) was used to verify m⁶A regulators in lung cancer cell lines. Protein-protein interaction (PPI), gene co-expression, survival analysis, and heatmap were used to analyze these m⁶A regulators in this set of lung cancer clinical data. Lasso regression was used to optimize the pathology-specific m⁶A regulator gene signature. Gene set enrichment analysis (GSEA) was used to reveal the functional characteristics of m⁶A regulators.

Results

Those 19 m⁶A regulator profiling was significantly differentially expressed in lung cancer tissues relative to control tissues, which was also verified in lung cancer cell lines. Those m⁶A regulators interacted mutually, and those regulator-based sample clusters were correlated with clinical traits, including survival status, gender, tobacco smoking history, primary disease, and pathologic stage. Further, lasso regression based on the 19 m⁶A regulators optimized and identified a three-m⁶A-regulator signature (KIAA1429, METTL3, and IGF2BP1) as independent prognostic factor, which classified 1013 lung cancer patients into high-risk and low-risk groups according to median value (0.84) of the lasso regression risk scores. This three-m⁶A-regulator signature profiling was significantly related to lung cancer overall survival, cancer status, and the above-described clinical traits. Further, GSEA revealed that KIAA1429, METTL3, and IGF2BP1 were significantly related to multiple biological behaviors, including proliferation, apoptosis, metastasis, energy metabolism, drug resistance, and recurrence, and that KIAA1429 and IGF2BP1 had potential target genes, including E2F3, WTAP, CCND1, CDK4, EGR2, YBX1, and TLX, which were associated with cancers.

Conclusion

This study provided the first view of the pathology-specific regulators of m⁶A RNA modification in lung cancers and identified the three-m⁶A-regulator signature (KIAA1429, METTL3, and IGF2BP1) as an independent prognostic model to classify lung cancers into high- and low-risk groups for patient stratification, prognostic assessment, and personalized treatment toward PPPM in lung cancers.

RAS/MAPK Pathway Driver Alterations Are Significantly Associated With Oncogenic KIT Mutations in Germ-cell Tumors

OBJECTIVE

To report the mutational profile and clinical outcomes of a cohort of patients with KIT-mutant seminomas and nonseminomatous germ-cell tumors (SGCT/NSGCTs).

PATIENTS AND METHODS

Retrospective cohort study of all patients with KIT-mutant GCTs sequenced at Memorial Sloan Kettering between March 2014 and March 2020. Tumors were assessed with MSK-IMPACT, a DNA next-generation sequencing assay for targeted sequencing of up to 468 key cancer genes.

RESULTS

Among 568 patients with GCTs, 8.1% had somatic KIT mutations, including 28 seminomas and 18 mixed/NSGCTs. Exons 17 (67.3%), 11 (22.4%), and 13 (6.1%) were most commonly affected. KIT-mutant cases were enriched for oncogenic RAS/MAPK pathway alterations compared to KIT-wildtype cases (34.8% vs 19.2%, P = .02). Among KIT-mutant cases, concurrent mutations were noted in KRAS (21.7%), RRAS2 (11.8%), CBL (6.5%), NRAS (4.3%), MAP2K1 (2.2%), and RAC1 (2.2%). Mutations in KRAS, RRAS2, and NRAS were mutually exclusive. In all, 73.9% of patients developed metastases and 95.7% received chemotherapy. No patients received KIT-directed tyrosine kinase inhibitors (TKIs). Classification as a NSGCT rather than a SGCT was associated with an increased risk of death (hazard ratio 9.1, 95% confidence interval 1.1-78.4, P = .04) while the presence of a concurrent RAS/MAPK pathway alteration was not (hazard ratio 0.8, 95% confidence interval 0.1-4.3, P = .76).

CONCLUSION

Mitogenic driver alterations can co-occur with activating KIT mutations, which may explain the lack of efficacy of KIT-directed TKIs in prior trials. Novel KIT-directed TKIs that target exon 17 mutations may benefit chemotherapy-refractory patients with KIT-mutant GCTs without RAS/MAPK alterations. Dual MEK/KIT inhibitor therapy in KIT-mutant GCTs with concurrent RAS/MAPK alterations could also be a plausible therapeutic strategy.

The Novel Zinc Finger Protein 587B Gene, ZNF587B, Regulates Cell Proliferation and Metastasis in Ovarian Cancer Cells in vivo and in vitro

Background

The zinc finger protein 587B (ZNF587B) is a novel cisplatin-sensitive gene that was identified in our previous research by using a genome-scale CRISPR-Cas9 knockout library in ovarian cancer (OC) cell lines. ZNF587B belongs to the C2H2-type zinc finger protein (ZFP) family. Many ZFP protein could inhibit tumor development and malignancy. However, the function of ZNF587B remains unknown.

Methods

Quantitative PCR (qPCR) was utilized to compare ZNF587B mRNA expression levels in OC and normal ovarian cell lines. The small interfering RNA (siRNA) and full-length ZNF587B eukaryotic expression plasmid were constructed and transfected into OC cells later. Colony formation, 5-ethynyl-2′-deoxyuridine (EdU) assay, transwell assay, and xenograft experiment were conducted to evaluate the effect of ZNF587B on OC cells.

Results

ZNF587B was downregulated by approximately 43% and 17% in the OC cell lines SKOV3 and A2780, respectively, compared with that in the normal ovarian cell line IOSE80. Overexpression of ZNF587B reduced cell proliferation, colony formation, migration, and invasion, which could be reversed by knockdown of ZNF587B via siRNA. Xenograft experiments also confirmed that ZNF587B could suppress tumor growth. Survival data of OC patients in the SurvExpress database showed that with respect to overall survival, low-risk patients grouped by the prognostic index had a higher expression of ZNF587B and a better prognosis than high-risk group (HR = 1.77, 95% CI: 0.55–0.70, p = 0.023). Moreover, overexpression of ZNF587B promoted OC cells apoptosis when pretreated with cisplatin.

Conclusion

ZNF587B is a novel potential tumor suppressor of OC and may be a therapeutic target for OC.

Rho GTPases in Gynecologic Cancers: In-Depth Analysis toward the Paradigm Change from Reactive to Predictive, Preventive, and Personalized Medical Approach Benefiting the Patient and Healthcare

Rho guanosine triphospatases (GTPases) resemble a conserved family of GTP-binding proteins regulating actin cytoskeleton dynamics and several signaling pathways central for the cell. Rho GTPases create a so-called Ras-superfamily of GTPases subdivided into subgroups comprising at least 20 members. Rho GTPases play a key regulatory role in gene expression, cell cycle control and proliferation, epithelial cell polarity, cell migration, survival, and apoptosis, among others. They also have tissue-related functions including angiogenesis being involved in inflammatory and wound healing processes. Contextually, any abnormality in the Rho GTPase function may result in severe consequences at molecular, cellular, and tissue levels. Rho GTPases also play a key role in tumorigenesis and metastatic disease. Corresponding mechanisms include a number of targets such as kinases and scaffold/adaptor-like proteins initiating GTPases-related signaling cascades. The accumulated evidence demonstrates the oncogenic relevance of Rho GTPases for several solid malignancies including breast, liver, bladder, melanoma, testicular, lung, central nervous system (CNS), head and neck, cervical, and ovarian cancers. Furthermore, Rho GTPases play a crucial role in the development of radio- and chemoresistance e.g. under cisplatin-based cancer treatment. This article provides an in-depth overview on the role of Rho GTPases in gynecological cancers, highlights relevant signaling pathways and pathomechanisms, and sheds light on their involvement in tumor progression, metastatic spread, and radio/chemo resistance. In addition, insights into a spectrum of novel biomarkers and innovative approaches based on the paradigm shift from reactive to predictive, preventive, and personalized medicine are provided.

High Levels of HIST1H2BK in Low-Grade Glioma Predicts Poor Prognosis: A Study Using CGGA and TCGA Data

A number of biomarkers have been identified for various cancers. However, biomarkers associated with glioma remain largely to be explored. In the current study, we investigated the relationship between the expression and prognostic value of the HIST1H2BK gene in glioma. Sequential data filtering (survival analysis, independent prognostic analysis, ROC curve analysis, and clinical correlation analysis) was performed, which resulted in identification of the association between the HIST1H2BK gene and glioma. Then, the HIST1H2BK gene was analyzed using bioinformatics (Kaplan–Meier survival analysis, univariate Cox analysis, multivariate Cox analysis, and ROC curve analysis). The results showed that low expression of HIST1H2BK was associated with better prognosis, and high expression of HIST1H2BK was associated with poor prognosis. In addition, HIST1H2BK was an independent prognostic indicator for patients with glioma. We also evaluated the association between HIST1H2BK and clinical characteristics. Furthermore, gene set enrichment analysis (GSEA) and analysis of immune infiltration were performed. The results showed that HIST1H2BK was associated with intensity of immune infiltration in glioma. Finally, co-expression analysis was performed. The results showed that HIST1H2BK was positively correlated with HIST1H2AG, HIST2H2AA4, HIST1H2BJ, HIST2H2BE, and HIST1H2AC, and negatively correlated with PDZD4, CRY2, GABBR1, rp5-1119a7.17, and KCNJ11. This study showed that upregulation of HIST1H2BK in low-grade glioma (LGG) tissue was an indicator of poor prognosis. Moreover, this study demonstrated that HIST1H2BK may be a promising biomarker for the treatment of LGG.

Calcium Channels as Novel Therapeutic Targets for Ovarian Cancer Stem Cells

Drug resistance in epithelial ovarian cancer (EOC) is reportedly attributed to the existence of cancer stem cells (CSC), because in most cancers, CSCs still remain after chemotherapy. To overcome this limitation, novel therapeutic strategies are required to prevent cancer recurrence and chemotherapy-resistant cancers by targeting cancer stem cells (CSCs). We screened an FDA-approved compound library and found four voltage-gated calcium channel blockers (manidipine, lacidipine, benidipine, and lomerizine) that target ovarian CSCs. Four calcium channel blockers (CCBs) decreased sphere formation, viability, and proliferation, and induced apoptosis in ovarian CSCs. CCBs destroyed stemness and inhibited the AKT and ERK signaling pathway in ovarian CSCs. Among calcium channel subunit genes, three L- and T-type calcium channel genes were overexpressed in ovarian CSCs, and downregulation of calcium channel genes reduced the stem-cell-like properties of ovarian CSCs. Expressions of these three genes are negatively correlated with the survival rate of patient groups. In combination therapy with cisplatin, synergistic effect was shown in inhibiting the viability and proliferation of ovarian CSCs. Moreover, combinatorial usage of manidipine and paclitaxel showed enhanced effect in ovarian CSCs xenograft mouse models. Our results suggested that four CCBs may be potential therapeutic drugs for preventing ovarian cancer recurrence.

Label-free quantitative identification of abnormally ubiquitinated proteins as useful biomarkers for human lung squamous cell carcinomas

BACKGROUND

Ubiquitination is an important molecular event in lung squamous cell carcinoma (LSCC), which currently is mainly studied in nonsmall cell lung carcinoma cell models but lacking of ubiquitination studies on LSCC tissues. Here, we presented the ubiquitinated protein profiles of LSCC tissues to explore ubiquitination-involved molecular network alterations and identify abnormally ubiquitinated proteins as useful biomarkers for predictive, preventive, and personalized medicine (PPPM) in LSCC.

METHODS

Anti-ubiquitin antibody-based enrichment coupled with LC-MS/MS was used to identify differentially ubiquitinated proteins (DUPs) between LSCC and control tissues, followed by integrative omics analyses to identify abnormally ubiquitinated protein biomarkers for LSCC.

RESULTS

Totally, 400 DUPs with 654 ubiquitination sites were identified,, and motifs A-X (1/2/3)-K* were prone to be ubiquitinated in LSCC tissues. Those DUPs were involved in multiple molecular network systems, including the ubiquitin-proteasome system (UPS), cell metabolism, cell adhesion, and signal transduction. Totally, 44 hub molecules were revealed by protein-protein interaction network analysis, followed by survival analysis in TCGA database (494 LSCC patients and 20,530 genes) to obtain 18 prognosis-related mRNAs, of which the highly expressed mRNAs VIM and IGF1R were correlated with poorer prognosis, while the highly expressed mRNA ABCC1 was correlated with better prognosis. VIM-encoded protein vimentin and ABCC1-encoded protein MRP1 were increased in LSCC, which were all associated with poor prognosis. Proteasome-inhibited experiments demonstrated that vimentin and MRP1 were degraded through UPS. Quantitative ubiquitinomics found ubiquitination level was decreased in vimentin and increased in MRP1 in LSCC. These findings showed that the increased vimentin in LSCC might be derived from its decreased ubiquitination level and that the increased MRP1 in LSCC might be derived from its protein synthesis > degradation. GSEA and co-expression gene analyses revealed that VIM and MRP1 were involved in multiple crucial biological processes and pathways. Further, TRIM2 and NEDD4L were predicted as E3 ligases to regulate ubiquitination of vimentin and MRP1, respectively.

CONCLUSION

These findings revealed ubiquitinomic variations and molecular network alterations in LSCC, which is in combination with multiomics analysis to identify ubiquitination-related biomarkers for in-depth insight into the molecular mechanism and therapeutic targets and for prediction, diagnosis, and prognostic assessment of LSCC.