Peroxisomes and cancer: The role of a metabolic specialist in a disease of aberrant metabolism

Proteomic insight towards key modulating proteins regulated by the aryl hydrocarbon receptor involved in ovarian carcinogenesis and chemoresistance

Gynecological malignancies pose a severe threat to female lives. Ovarian cancer (OC), the most lethal gynecological malignancy, is clinically presented with chemoresistance and a higher relapse rate. Several studies have highly correlated the incidence of OC to exposure to environmental pollutants, such as 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), a process mainly mediated through activating the aryl hydrocarbon receptor (AhR). We have previously reported that exposure of OC cells to TCDD, an AhR activator, significantly modulated the expression of several genes that play roles in stemness and chemoresistance. However, the effect of AhR activation on the whole OC cell proteome aiming at identifying novel druggable targets for both prevention and treatment intervention purposes remains unrevealed. For this purpose, we conducted a comparative proteomic analysis of OC cells A2780 untreated/treated with TCDD for 24 h using a mass spectrometry-based label-free shotgun proteomics approach. The most significantly dysregulated proteins were validated by Western blot analysis. Our results showed that upon AhR activation by TCDD, out of 2598 proteins identified, 795 proteins were upregulated, and 611 were downregulated. STRING interaction analysis and KEGG-Reactome pathway analysis approaches identified several significantly dysregulated proteins that were categorized to be involved in chemoresistance, cancer progression, invasion and metastasis, apoptosis, survival, and prognosis in OC. Importantly, selected dysregulated genes identified by the proteomic study were validated at the protein expression levels by Western blot analysis. In conclusion, this study provides a better understanding of the the cross-talk between AhR and several other molecular signaling pathways and the role and involvement of AhR in ovarian carcinogenesis and chemoresistance. Moreover, the study suggests that AhR is a potential therapeutic target for OC prevention and maintenance. SIGNIFICANCE: To our knowledge, this is the first study that investigates the role and involvement of AhR and its regulated genes in OC by performing a comparative proteomic analysis to identify the critical proteins with a modulated expression upon AhR activation. We found AhR activation to play a tumor-promoting and chemoresistance-inducing role in the pathogenesis of OC. The results of our study help to devise novel therapeutics for better management and prevention and open the doors to finding novel biomarkers for the early detection and prognosis of OC.

Identification and Validation of a PEX5-Dependent Signature for Prognostic Prediction in Glioma

Gliomas, the most prevalent and lethal form of brain cancer, are known to exhibit metabolic alterations that facilitate tumor growth, invasion, and resistance to therapies. Peroxisomes, essential organelles responsible for fatty acid oxidation and reactive oxygen species (ROS) homeostasis, rely on the receptor PEX5 for the import of metabolic enzymes into their matrix. However, the prognostic significance of peroxisomal enzymes for glioma patients remains unclear. In this study, we elucidate that PEX5 is indispensable for the cell growth, migration, and invasion of glioma cells. We establish a robust prognosis model based on the expression of peroxisomal enzymes, whose localization relies on PEX5. This PEX5-dependent signature not only serves as a robust prognosis model capable of accurately predicting outcomes for glioma patients, but also effectively distinguishes several clinicopathological features, including the grade, isocitrate dehydrogenase (IDH) mutation, and 1p19q codeletion status. Furthermore, we developed a nomogram that integrates the prognostic model with other clinicopathological factors, demonstrating highly accurate performance in estimating patient survival. Patients classified into the high-risk group based on our prognostic model exhibited an immunosuppressive microenvironment. Finally, our validation reveals that the elevated expression of GSTK1, an antioxidant enzyme within the signature, promotes the cell growth and migration of glioma cells, with this effect dependent on the peroxisomal targeting signal recognized by PEX5. These findings identify the PEX5-dependent signature as a promising prognostic tool for gliomas.

A novel peroxisome-related gene signature predicts clinical prognosis and is associated with immune microenvironment in low-grade glioma

Low-grade glioma (LGG), a common primary tumor, mainly originates from astrocytes and oligodendrocytes. Increasing evidence has shown that peroxisomes function in the regulation of tumorigenesis and development of cancer. However, the prognostic value of peroxisome-related genes (PRGs) in LGG has not been reported. Therefore, it is necessary to construct a prognostic risk model for LGG patients based on the expression profiles of peroxisome-related genes. Our study mainly concentrated on developing a peroxisome-related gene signature for overall survival (OS) prediction in LGG patients. First, according to these peroxisome-related genes, all LGG patients from The Cancer Genome Atlas (TCGA) database could be divided into two subtypes. Univariate Cox regression analysis was used to find prognostic peroxisome-related genes in TCGA_LGG dataset, and least absolute shrinkage and selection operator Cox regression analysis was employed to establish a 14-gene signature. The risk score based on the signature was positively associated with unfavorable prognosis. Then, multivariate Cox regression incorporating additional clinical characteristics showed that the 14-gene signature was an independent predictor of LGG. Time-dependent ROC curves revealed good performance of this prognostic signature in LGG patients. The performance about predicting OS of LGG was validated using the GSE107850 dataset derived from the Gene Expression Omnibus (GEO) database. Furethermore, we constructed a nomogram model based on the gene signature and age, which showed a better prognostic power. Gene ontology (GO) and Kyoto Encylopedia of Genes and Genomes (KEGG) analyses showed that neuroactive ligand-receptor interaction and phagosome were enriched and that the immune status was decreased in the high-risk group. Finally, cell counting kit-8 (CCK8) were used to detect cell proliferation of U251 and A172 cells. Inhibition of ATAD1 (ATPase family AAA domain-containing 1) and ACBD5 (Acyl-CoA binding-domain-containing-5) expression led to significant inhibition of U251 and A172 cell proliferation. Flow cytometry detection showed that ATAD1 and ACBD5 could induce apoptosis of U251 and A172 cells. Therefore, through bioinformatics methods and cell experiments, our study developed a new peroxisome-related gene signature that migh t help improve personalized OS prediction in LGG patients.

High-content image screening to identify chemical modulators for peroxisome and ferroptosis

BACKGROUND

The peroxisome is a dynamic organelle with variety in number, size, shape, and activity in different cell types and physiological states. Recent studies have implicated peroxisomal homeostasis in ferroptosis susceptibility. Here, we developed a U-2OS cell line with a fluorescent peroxisomal tag and screened a target-selective chemical library through high-content imaging analysis.

METHODS

U-2OS cells stably expressing the mOrange2-Peroxisomes2 tag were generated to screen a target-selective inhibitor library. The nuclear DNA was counterstained with Hoechst 33342 for cell cycle analysis. Cellular images were recorded and quantitatively analyzed through a high-content imaging platform. The effect of selected compounds on ferroptosis induction was analyzed in combination with ferroptosis inducers (RSL3 and erastin). Flow cytometry analysis was conducted to assess the level of reactive oxygen species (ROS) and cell death events.

RESULTS

Through the quantification of DNA content and peroxisomal signals in single cells, we demonstrated that peroxisomal abundance was closely linked with cell cycle progression and that peroxisomal biogenesis mainly occurred in the G1/S phase. We further identified compounds that positively and negatively regulated peroxisomal abundance without significantly affecting the cell cycle distribution. Some compounds promoted peroxisomal signals by inducing oxidative stress, while others regulated peroxisomal abundance independent of redox status. Importantly, compounds with peroxisome-enhancing activity potentiated ferroptosis induction.

CONCLUSIONS

Our findings pinpoint novel cellular targets that might be involved in peroxisome homeostasis and indicate that compounds promoting peroxisomal abundance could be jointly applied with ferroptosis inducers to potentiate anticancer effect.

Peroxisomal hydrogen peroxide signaling: A new chapter in intracellular communication research

Hydrogen peroxide (H2O2), a natural metabolite commonly found in aerobic organisms, plays a crucial role in numerous cellular signaling processes. One of the key organelles involved in the cell's metabolism of H2O2 is the peroxisome. In this review, we first provide a concise overview of the current understanding of H2O2 as a molecular messenger in thiol redox signaling, along with the role of peroxisomes as guardians and modulators of cellular H2O2 balance. Next, we direct our focus toward the recently identified primary protein targets of H2O2 originating from peroxisomes, emphasizing their importance in unraveling the complex interplay between peroxisomal H2O2 and cell signaling. We specifically focus on three areas: signaling through peroxiredoxin redox relay complexes, calcium signaling, and phospho-signaling. Finally, we highlight key research directions that warrant further investigation to enhance our comprehension of the molecular and biochemical mechanisms linking alterations in peroxisomal H2O2 metabolism with disease.

Establishment and validation of a novel peroxisome-related gene prognostic risk model in kidney clear cell carcinoma

BACKGROUND

Kidney clear cell carcinoma (KIRC) is the most common subtype of renal cell carcinoma. Peroxisomes play a role in the regulation of tumorigenesis and cancer progression, yet the prognostic significance of peroxisome-related genes (PRGs) remains rarely studied. The study aimed to establish a novel prognostic risk model and identify potential biomarkers in KIRC.

METHODS

The significant prognostic PRGs were screened through differential and Cox regression analyses, and LASSO Cox regression analysis was performed to establish a prognostic risk model in the training cohort, which was validated internally in the testing and entire cohorts, and further assessed in the GSE22541 cohort. Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed to explore the function and pathway differences between the high-risk and low-risk groups. The relationship between risk score and immune cell infiltration levels was evaluated in the CIBERSORT, ESTIMATE and TIMER databases. Finally, potential biomarkers were identified and validated from model genes, using immunohistochemistry.

RESULTS

Fourteen significant prognostic PRGs were identified using multiple analyses, and 9 genes (ABCD1, ACAD11, ACAT1, AGXT, DAO, EPHX2, FNDC5, HAO1, and HNGCLL1) were obtained to establish a prognostic model via LASSO Cox regression analysis. Combining the risk score with clinical factors to construct a nomogram, which provided support for personalized treatment protocols for KIRC patients. GO and KEGG analyses highlighted associations with substance metabolism, transport, and the PPAR signaling pathways. Tumor immune infiltration indicated immune suppression in the high-risk group, accompanied by higher tumor purity and the expression of 9 model genes was positively correlated with the level of immune cell infiltration. ACAT1 has superior prognostic capabilities in predicting the outcomes of KIRC patients.

CONCLUSIONS

The peroxisome-related prognostic risk model could better predict prognosis in KIRC patients.

Ether-lipids and cellular signaling: A differential role of alkyl- and alkenyl-ether-lipids?

Ether-lipids (EL) are specific lipids bearing a characteristic sn-1 ether bond. Depending on the ether or vinyl-ether nature of this bond, they are present as alkyl- or alkenyl-EL, respectively. Among EL, alkenyl-EL, also referred as plasmalogens in the literature, attract most of the scientific interest as they are the predominant EL species in eukaryotic cells, thus less is known about alkyl-EL. EL have been implicated in various signaling pathways and alterations in their quantity are frequently observed in pathologies such as neurodegenerative and cardiovascular diseases or cancer. However, it remains unknown whether both alkyl- and alkenyl-EL play the same roles in these processes. This review summarizes the roles and mechanisms of action of EL in cellular signaling and tries to discriminate between alkyl- and alkenyl-EL. We also focus on the involvement of EL-mediated alterations of cellular signaling in diseases and discuss the potential interest for EL in therapy.

Peroxisome disruption alters lipid metabolism and potentiates antitumor response with MAPK-targeted therapy in melanoma

Melanomas reprogram their metabolism to rapidly adapt to therapy-induced stress conditions, allowing them to persist and ultimately develop resistance. We report that a subpopulation of melanoma cells tolerate MAPK pathway inhibitors (MAPKis) through a concerted metabolic reprogramming mediated by peroxisomes and UDP-glucose ceramide glycosyltransferase (UGCG). Compromising peroxisome biogenesis, by repressing PEX3 expression, potentiated the proapoptotic effects of MAPKis via an induction of ceramides, an effect limited by UGCG-mediated ceramide metabolism. Cotargeting PEX3 and UGCG selectively eliminated a subset of metabolically active, drug-tolerant CD36+ melanoma persister cells, thereby sensitizing melanoma to MAPKis and delaying resistance. Increased levels of peroxisomal genes and UGCG were found in patient-derived MAPKi-relapsed melanomas, and simultaneously inhibiting PEX3 and UGCG restored MAPKi sensitivity in multiple models of therapy resistance. Finally, combination therapy consisting of a newly identified inhibitor of the PEX3-PEX19 interaction, a UGCG inhibitor, and MAPKis demonstrated potent antitumor activity in preclinical melanoma models, thus representing a promising approach for melanoma treatment.

GPR37 expression as a prognostic marker in gliomas: a bioinformatics-based analysis

BACKGROUND

Gliomas are the most frequently diagnosed primary brain tumors, and are associated with multiple molecular aberrations during their development and progression. GPR37 is an orphan G protein-coupled receptor (GPCR) that is implicated in different physiological pathways in the brain, and has been linked to various malignancies. The aim of this study was to explore the relationship between GPR37 gene expression and the clinicopathological factors, patient prognosis, tumor-infiltrating immune cell signature GSEA and methylation levels in glioma.

METHODS

We explored the diagnostic value, clinical relevance, and molecular function of GPR37 in glioma using TCGA, STRING, cBioPortal, Tumor Immunity Estimation Resource (TIMER) database and MethSurv databases. Besides, the "ssGSEA" algorithm was conducted to estimate immune cells infiltration abundance, with 'ggplot2' package visualizing the results. Immunohistochemical staining of clinical samples were used to verify the speculations of bioinformatics analysis.

RESULTS

GPR37 expression was significantly higher in the glioma tissues compared to the normal brain tissues, and was linked to poor prognosis. Functional annotation of GPR37 showed enrichment of ether lipid metabolism, fat digestion and absorption, and histidine metabolism. In addition, GSEA showed that GPR37 was positively correlated to the positive regulation of macrophage derived foam cell differentiation, negative regulation of T cell receptor signaling pathway, neuroactive ligand receptor interaction, calcium signaling pathway, and negatively associated with immunoglobulin complex, immunoglobulin complex circulating, ribosome and spliceosome mediated by circulating immunoglobulin etc. TIMER2.0 and ssGSEA showed that GPR37 expression was significantly associated with the infiltration of T cells, CD8 T cell, eosinophils, macrophages, neutrophils, NK CD56dim cells, NK cells, plasmacytoid DCs (pDCs), T helper cells and T effector memory (Tem) cells. In addition, high GPR37 expression was positively correlated with increased infiltration of M2 macrophages, which in turn was associated with poor prognosis. Furthermore, GPR37 was positively correlated with various immune checkpoints (ICPs). Finally, hypomethylation of the GPR37 promoter was associated with its high expression levels and poor prognosis in glioma.

CONCLUSION

GPR37 had diagnostic and prognostic value in glioma. The possible biological mechanisms of GPR37 provide novel insights into the clinical diagnosis and treatment of glioma.

Structure of the peroxisomal Pex1/Pex6 ATPase complex bound to a substrate

The double-ring AAA+ ATPase Pex1/Pex6 is required for peroxisomal receptor recycling and is essential for peroxisome formation. Pex1/Pex6 mutations cause severe peroxisome associated developmental disorders. Despite its pathophysiological importance, mechanistic details of the heterohexamer are not yet available. Here, we report cryoEM structures of Pex1/Pex6 from Saccharomyces cerevisiae, with an endogenous protein substrate trapped in the central pore of the catalytically active second ring (D2). Pairs of Pex1/Pex6(D2) subdomains engage the substrate via a staircase of pore-1 loops with distinct properties. The first ring (D1) is catalytically inactive but undergoes significant conformational changes resulting in alternate widening and narrowing of its pore. These events are fueled by ATP hydrolysis in the D2 ring and disengagement of a "twin-seam" Pex1/Pex6(D2) heterodimer from the staircase. Mechanical forces are propagated in a unique manner along Pex1/Pex6 interfaces that are not available in homo-oligomeric AAA-ATPases. Our structural analysis reveals the mechanisms of how Pex1 and Pex6 coordinate to achieve substrate translocation.

Preexisting Immunity Drives the Response to Neoadjuvant Chemotherapy in Esophageal Adenocarcinoma

Current treatment for patients with locally advanced esophageal adenocarcinoma (EAC) is neoadjuvant chemotherapy (nCT), alone or combined with radiotherapy, before surgery. However, fewer than 30% of treated patients show a pathologic complete response to nCT, which correlates with increased 5-year survival compared with nonresponders. Understanding the mechanisms of response to nCT is pivotal to better stratify patients and inform more efficacious therapies. Here, we investigated the immune mechanisms involved in nCT response by multidimensional profiling of pretreatment tumor biopsies and blood from 68 patients with EAC (34 prospectively and 34 retrospectively collected), comparing complete responders versus nonresponders to nCT. At the tumor level, complete response to nCT was associated with molecular signatures of immune response and proliferation, increased putative antitumor tissue-resident memory CD39+ CD103+ CD8+ T cells, and reduced immunosuppressive T regulatory cells (Treg) and M2-like macrophages. Systemically, complete responders showed higher frequencies of immunostimulatory CD14+ CD11c+ HLA-DRhigh cells, and reduced programmed cell death ligand 1-positive (PD-L1+) monocytic myeloid-derived suppressor cells, along with high plasma GM-CSF (proinflammatory) and low IL4, CXCL10, C3a, and C5a (suppressive). Plasma proinflammatory and suppressive cytokines correlated directly and inversely, respectively, with the frequency of tumor-infiltrating CD39+ CD103+ CD8+ T cells. These results suggest that preexisting immunity in baseline tumor drives the clinical activity of nCT in locally advanced EAC. Furthermore, it may be possible to stratify patients based on predictive immune signatures, enabling tailored neoadjuvant and/or adjuvant regimens.

SIGNIFICANCE

Multidimensional profiling of pretreatment esophageal adenocarcinoma shows patient response to nCT is correlated with active preexisting immunity and indicates molecular pathways of resistance that may be targeted to improve clinical outcomes.

Marine Plasmalogens: A Gift from the Sea with Benefits for Age-Associated Diseases

Aging increases oxidative and inflammatory stress caused by a reduction in metabolism and clearance, thus leading to the development of age-associated diseases. The quality of our daily diet and exercise is important for the prevention of these diseases. Marine resources contain various valuable nutrients, and unique glycerophospholipid plasmalogens are found abundantly in some marine invertebrates, including ascidians. One of the major classes, the ethanolamine class (PlsEtn), exists in a high ratio to phospholipids in the brain and blood, while decreased levels have been reported in patients with age-associated diseases, including Alzheimer's disease. Animal studies have shown that the administration of marine PlsEtn prepared from marine invertebrates improved PlsEtn levels in the body and alleviated inflammation. Animal and human studies have reported that marine PlsEtn ameliorates cognitive impairment. In this review, we highlight the biological significance, relationships with age-associated diseases, food functions, and healthcare materials of plasmalogens based on recent knowledge and discuss the contribution of marine plasmalogens to health maintenance in aging.

Evaluating the strength of evidence for genes implicated in peroxisomal disorders using the ClinGen clinical validity framework and providing updates to the peroxisomal disease nomenclature

Peroxisomal disorders are heterogeneous in nature, with phenotypic overlap that is indistinguishable without molecular testing. Newborn screening and gene sequencing for a panel of genes implicated in peroxisomal diseases are critical tools for the early and accurate detection of these disorders. It is therefore essential to evaluate the clinical validity of the genes included in sequencing panels for peroxisomal disorders. The Peroxisomal Gene Curation Expert Panel (GCEP) assessed genes frequently included on clinical peroxisomal testing panels using the Clinical Genome Resource (ClinGen) gene-disease validity curation framework and classified gene-disease relationships as Definitive, Strong, Moderate, Limited, Disputed, Refuted, or No Known Disease Relationship. Subsequent to gene curation, the GCEP made recommendations to update the disease nomenclature and ontology in the Monarch Disease Ontology (Mondo) database. Thirty-six genes were assessed for the strength of evidence supporting their role in peroxisomal disease, leading to 36 gene-disease relationships, after two genes were removed for their lack of a role in peroxisomal disease and two genes were curated for two different disease entities each. Of these, 23 were classified as Definitive (64%), one as Strong (3%), eight as Moderate (23%), two as Limited (5%), and two as No known disease relationship (5%). No contradictory evidence was found to classify any relationships as Disputed or Refuted. The gene-disease relationship curations are publicly available on the ClinGen website (https://clinicalgenome.org/affiliation/40049/). The changes to peroxisomal disease nomenclature are displayed on the Mondo website (http://purl.obolibrary.org/obo/MONDO_0019053). The Peroxisomal GCEP-curated gene-disease relationships will inform clinical and laboratory diagnostics and enhance molecular testing and reporting. As new data will emerge, the gene-disease classifications asserted by the Peroxisomal GCEP will be re-evaluated periodically.

Treatment Strategies for KRAS-Mutated Non-Small-Cell Lung Cancer

Activating mutations in KRAS are highly prevalent in solid tumours and are frequently found in 35% of lung, 45% of colorectal, and up to 90% of pancreatic cancers. Mutated KRAS is a prognostic factor for disease-free survival (DFS) and overall survival (OS) in NSCLC and is associated with a more aggressive clinical phenotype, highlighting the need for KRAS-targeted therapy. Once considered undruggable due to its smooth shallow surface, a breakthrough showed that the activated G12C-mutated KRAS isozyme can be directly inhibited via a newly identified switch II pocket. This discovery led to the development of a new class of selective small-molecule inhibitors against the KRAS G12C isoform. Sotorasib and adagrasib are approved in locally advanced or metastatic NSCLC patients who have received at least one prior systemic therapy. Currently, there are at least twelve KRAS G12C inhibitors being tested in clinical trials, either as a single agent or in combination. In this study, KRAS mutation prevalence, subtypes, rates of occurrence in treatment-resistant invasive mucinous adenocarcinomas (IMAs), and novel drug delivery options are reviewed. Additionally, the current status of KRAS inhibitors, multiple resistance mechanisms that limit efficacy, and their use in combination treatment strategies and novel multitargeted approaches in NSCLC are discussed.

Pre-meiotic deletion of PEX5 causes spermatogenesis failure and infertility in mice

Peroxisomes are involved in the regulation of various pathological processes. Peroxisomal biogenesis factor 5 (PEX5), which plays an essential role in peroxisomal biogenesis, is critical for reactive oxygen species (ROS) accumulation. However, its underlying functions in spermatogenesis have not yet been identified. Pex5 was deleted by crossing Stra8-Cre mice with Pex5^flox/flox mice before the onset of meiosis. The morphology of testes and epididymides, spermatogenesis function, and fertility in both wild type (WT) and Pex5^-/- mice were analysed by haematoxylin and eosin (HE) and immunofluorescent staining. Mechanism of PEX5 affecting peroxisomes and spermatogenesis were validated by Western blot and transmission electron microscopy (TEM). Transcriptome RNA sequencing (RNA-seq) was used to profile the dysregulated genes in testes from WT and Pex5^-/- mice on postnatal day (P) 35. The adult Pex5 knockout male mice were completely sterile with no mature sperm production. Loss of Pex5 in spermatocytes resulted in multinucleated giant cell formation, meiotic arrest, abnormal tubulin expression, and deformed acrosome formation. Furthermore, Pex5 deletion led to delayed DNA double-strand break repair and improper crossover at the pachytene stage. Impaired peroxisome function in Pex5 knockout mice induced ROS redundancy, which in turn led to an increase in germ cell apoptosis and a decline in autophagy. Pex5 regulates ROS during meiosis and is essential for spermatogenesis and male fertility in mice.

Good things come to those who bait: the peroxisomal docking complex

Peroxisomal integrity and function are highly dependent on its membrane and soluble (matrix) components. Matrix enzymes are imported post-translationally in a folded or even oligomeric state, via a still mysterious protein translocation mechanism. They are guided to peroxisomes via the Peroxisomal Targeting Signal (PTS) sequences which are recognized by specific cytosolic receptors, Pex5, Pex7 and Pex9. Subsequently, cargo-loaded receptors bind to the docking complex in an initial step, followed by channel formation, cargo-release, receptor-recycling and -quality control. The docking complexes of different species share Pex14 as their core component but differ in composition and oligomeric state of Pex14. Here we review and highlight the latest insights on the structure and function of the peroxisomal docking complex. We summarize differences between yeast and mammals and then we integrate this knowledge into our current understanding of the import machinery.

Circular RNA Sec61 subunit alpha isoform 1 by competitive absorption of microRNA-513a-5p mediates peroxisomal biogenesis factor 5 expression and promotes the malignant phenotype of non-small cell lung cancer

Circular RNAs (circRNAs) are functional RNAs in the development and metabolism of non-small cell lung cancer (NSCLC). Therein, this paper particularly elucidated the circRNA SEC61 subunit alpha isoform 1 (circSEC61A1) in NSCLC has not been fully elucidated. Clinical analysis of circSEC61A1 expression was performed on specimens collected from 51 patients with primary NSCLC, together with patients' survival. Cell experiments were performed after interfering with circSEC61A1, microRNA (miR)-513a-5p, and peroxisomal biogenesis factor 5 (PEX5) expression, respectively, and cell malignant phenotypes and aerobic glycolysis were evaluated, as well as epithelial-to-mesenchymal transition (EMT)-related markers and Wnt/β-catenin pathway. Xenografts experiments studied the performance of circSEC61A1 in vivo. The downstream molecules of circSEC61A1 were searched. Our data demonstrated that circSEC61A1 was upregulated in NSCLC patients, showing an association with poorer survival outcomes. In cell experiments, circSEC61A1 overexpression promoted NSCLC malignant phenotypes, glycolysis, EMT, and Wnt/β-catenin pathway activation, whereas circSEC61A1 underexpression did the opposite. Knockdown of circSEC61A1 limited tumor growth and metastasis. Furthermore, circSEC61A1 could regulate PEX5 expression through competitive absorption of miR-513a-5p. Generally, circSEC61A1 is a potential biomarker for NSCLC, and circSEC61A1 serves tumor-promoting action in the progression of NSCLC.

Organelle-targeted therapies: a comprehensive review on system design for enabling precision oncology

Cancer is a major threat to human health. Among various treatment methods, precision therapy has received significant attention since the inception, due to its ability to efficiently inhibit tumor growth, while curtailing common shortcomings from conventional cancer treatment, leading towards enhanced survival rates. Particularly, organelle-targeted strategies enable precise accumulation of therapeutic agents in organelles, locally triggering organelle-mediated cell death signals which can greatly reduce the therapeutic threshold dosage and minimize side-effects. In this review, we comprehensively discuss history and recent advances in targeted therapies on organelles, specifically including nucleus, mitochondria, lysosomes and endoplasmic reticulum, while focusing on organelle structures, organelle-mediated cell death signal pathways, and design guidelines of organelle-targeted nanomedicines based on intervention mechanisms. Furthermore, a perspective on future research and clinical opportunities and potential challenges in precision oncology is presented. Through demonstrating recent developments in organelle-targeted therapies, we believe this article can further stimulate broader interests in multidisciplinary research and technology development for enabling advanced organelle-targeted nanomedicines and their corresponding clinic translations.

Inhibition of BRD4 Promotes Pexophagy by Increasing ROS and ATM Activation

Although autophagy regulates the quality and quantity of cellular compartments, the regulatory mechanisms underlying peroxisomal autophagy (pexophagy) remain largely unknown. In this study, we identified several BRD4 inhibitors, including molibresib, a novel pexophagy inducer, via chemical library screening. Treatment with molibresib promotes loss of peroxisomes selectively, but not mitochondria, ER, or Golgi apparatus in HeLa cells. Consistently, depletion of BRD4 expression also induced pexophagy in RPE cells. In addition, the inhibition of BRD4 by molibresib increased autophagic degradation of peroxisome ATG7-dependency. We further found that molibresib produced reactive oxygen species (ROS), which potentiates ATM activation. Inhibition of ROS or ATM suppressed the loss of peroxisomes in molibresib-treated cells. Taken together, our data suggest that inhibition of BRD4 promotes pexophagy by increasing ROS and ATM activation.

Classification of stomach adenocarcinoma based on fatty acid metabolism-related genes frofiling

Background: Fatty acid metabolism (FAM)-related genes play a key role in the development of stomach adenocarcinoma (STAD). Although immunotherapy has led to a paradigm shift in STAD treatment, the overall response rate of immunotherapy for STAD is low due to heterogeneity of the tumor immune microenvironment (TIME). How FAM-related genes affect TIME in STAD remains unclear.

Methods: The univariate Cox regression analysis was performed to screen prognostic FAM-related genes using transcriptomic profiles of the Cancer Genome Atlas (TCGA)-STAD cohort. Next, the consensus clustering analysis was performed to divide the STAD cohort into two groups based on the 13 identified prognostic genes. Then, gene set enrichment analysis (GSEA) was carried out to identify enriched pathways in the two groups. Furthermore, we developed a prognostic signature model based on 7 selected prognostic genes, which was validated to be capable in predicting the overall survival (OS) of STAD patients using the univariate Cox regression, least absolute shrinkage and selection operator (LASSO) regression, and multivariate Cox regression analyses. Finally, the “Estimation of STromal and Immune cells in MAlignant Tumours using Expression data” (ESTIMATE) algorithm was used to evaluate the stromal, immune, and ESTIMATE scores, and tumor purity of each STAD sample.

Results: A total of 13 FAM-related genes were identified to be significantly associated with OS in STAD patients. Two molecular subtypes, which we named Group 1 and Group 2, were identified based on these FAM-related prognostic genes using the consensus clustering analysis. We showed that Group 2 was significantly correlated with poor prognosis and displayed higher programmed cell death ligand 1 (PD-L1) expressions and distinct immune cell infiltration patterns. Furthermore, using GSEA, we showed that apoptosis and HCM signaling pathways were significantly enriched in Group 2. We constructed a prognostic signature model using 7 selected FAM-related prognostic genes, which was proven to be effective for prediction of STAD (HR = 1.717, 95% CI = 1.105–1.240, p < 0.001). After classifying the patients into the high- and low-risk groups based on our model, we found that patients in the high-risk group tend to have more advanced T stages and higher tumor grades, as well as higher immune scores. We also found that the risk scores were positively correlated with the infiltration of certain immune cells, including resting dendritic cells (DCs), and M2 macrophages. We also demonstrated that elevated expression of gamma-glutamyltransferase 5 (GGT5) is significantly associated with worse OS and disease-free survival (DFS), more advanced T stage and higher tumor grade, and increased immune cell infiltration, suggesting that STAD patients with high GGT5 expression in the tumor tissues might have a better response to immunotherapy.

Conclusion: FAM-related genes play critical roles in STAD prognosis by shaping the TIME. These genes can regulate the infiltration of various immune cells and thus are potential therapeutic targets worthy of further investigation. Furthermore, GGT5 was a promising marker for predicting immunotherapeutic response in STAD patients.

Altered expression of ACOX2 in non-small cell lung cancer

Peroxisomes are organelles that play essential roles in many metabolic processes, but also play roles in innate immunity, signal transduction, aging and cancer. One of the main functions of peroxisomes is the processing of very-long chain fatty acids into metabolites that can be directed to the mitochondria. One key family of enzymes in this process are the peroxisomal acyl-CoA oxidases (ACOX1, ACOX2 and ACOX3), the expression of which has been shown to be dysregulated in some cancers. Very little is however known about the expression of this family of oxidases in non-small cell lung cancer (NSCLC). ACOX2 has however been suggested to be elevated at the mRNA level in over 10% of NSCLC, and in the present study using both standard and bioinformatics approaches we show that expression of ACOX2 is significantly altered in NSCLC. ACOX2 mRNA expression is linked to a number of mutated genes, and associations between ACOX2 expression and tumour mutational burden and immune cell infiltration were explored. Links between ACOX2 expression and candidate therapies for oncogenic driver mutations such as KRAS were also identified. Furthermore, levels of acyl-CoA oxidases and other associated peroxisomal genes were explored to identify further links between the peroxisomal pathway and NSCLC. The results of this biomarker driven study suggest that ACOX2 may have potential clinical utility in the diagnosis, prognosis and stratification of patients into various therapeutically targetable options.

Identifying cancer tissue-of-origin by a novel machine learning method based on expression quantitative trait loci

Cancer of unknown primary (CUP) refers to cancer with primary lesion unidentifiable by regular pathological and clinical diagnostic methods. This kind of cancer is extremely difficult to treat, and patients with CUP usually have a very short survival time. Recent studies have suggested that cancer treatment targeting primary lesion will significantly improve the survival of CUP patients. Thus, it is critical to develop accurate yet fast methods to infer the tissue-of-origin (TOO) of CUP. In the past years, there are a few computational methods to infer TOO based on single omics data like gene expression, methylation, somatic mutation, and so on. However, the metastasis of tumor involves the interaction of multiple levels of biological molecules. In this study, we developed a novel computational method to predict TOO of CUP patients by explicitly integrating expression quantitative trait loci (eQTL) into an XGBoost classification model. We trained our model with The Cancer Genome Atlas (TCGA) data involving over 7,000 samples across 20 types of solid tumors. In the 10-fold cross-validation, the prediction accuracy of the model with eQTL was over 0.96, better than that without eQTL. In addition, we also tested our model in an independent data downloaded from Gene Expression Omnibus (GEO) consisting of 87 samples across 4 cancer types. The model also achieved an f1-score of 0.7-1 depending on different cancer types. In summary, eQTL was an important information in inferring cancer TOO and the model might be applied in clinical routine test for CUP patients in the future.

A New Prognostic Indicator of Immune Microenvironment and Therapeutic Response in Lung Adenocarcinoma Based on Peroxisome-Related Genes

Lung adenocarcinoma (LUAD) has been the major cause of tumor-associated mortality in recent years and exhibits a poor outcome. New data revealed that peroxisomes have a function in the regulation of the development and progression of several tumors. However, the prognostic values of peroxisome-related genes (PRGs) were rarely reported. Genomic sequence, mutation, and clinical data of 535 LUAD tissues were obtained from TCGA data sets. Within the TCGA cohort, a multigene signature was constructed with the assistance of the LASSO Cox regression model. Three GEO data sets, including GSE3141, GSE31210, and GSE72094, were obtained as validation cohorts. ROC assays, Kaplan-Meier methods, and multivariate assays were applied to examine the prognostic capacities of the novel signature. Gene Set Enrichment Analysis (GSEA) was performed to further understand the underlying molecular mechanisms. In this study, we identified 47 differentially expressed peroxisome-related genes (PRGs), including 25 increased and 22 decreased PRGs. A prognostic model of six PRGs was established. The univariate and multivariate Cox analyses both showed that the p value of risk score was less than 0.05. In LUAD patients, the strong connection between the risk score and overall survival was further verified in three other GEO data sets. TMB and cancer stem cell infiltration were shown to be significantly higher in the high-risk group in comparison to the low-risk group. The TIDE score of the group with the low risk was considerably greater than that of the group with the high risk. Several drugs, targeting PRG-related genes, were available for the treatments of LUAD. Overall, we developed a novel peroxisome-related prognostic signature for LUAD patients. This signature could successfully indicate LUAD patients' chances of survival as well as their immune system's responsiveness to treatments. In addition, it has the potential to serve as immunotherapeutic targets for LUAD patients.

Genetic variants in DDO and PEX5L in peroxisome-related pathways predict non-small cell lung cancer survival

Peroxisomes play a role in lipid metabolism and regulation of reactive oxygen species, but its role in development and progression of non-small cell lung cancer (NSCLC) is not well understood. Here, we investigated the associations between 9708 single-nucleotide polymorphisms (SNPs) in 113 genes in the peroxisome-related pathways and survival of NSCLC patients from the Prostate, Lung, Colorectal, and Ovarian Cancer Screening Trial (PLCO) and the Harvard Lung Cancer Susceptibility (HLCS) study. In 1185 NSCLC patients from the PLCO trial, we found that 213 SNPs were significantly associated with NSCLC overall survival (OS) (p ≤ 0.05, Bayesian false discovery probability [BFDP] ≤ 0.80), of which eight SNPs were validated in the HLCS data set. In a multivariate Cox proportional hazards regression model, two independent SNPs (rs9384742 DDO and rs9825224 PEX5L) were significantly associated with NSCLC survival (hazards ratios [HR] of 1.17 with 95% CI [confidence interval] of 1.06-1.28 and 0.86 with 95% CI of 0.77-0.96, respectively). Patients with one or two protective genotypes had a significantly higher OS (HR: 0.787 [95% CI: 0.620-0.998] and 0.691 [95% CI: 0.543-0.879], respectively). Further expression quantitative trait loci analysis using whole blood and lung tissue showed that the minor allele of rs9384742 DDO was significantly associated with decreased messenger RNA (mRNA) expression levels and that DDO expression was also decreased in NSCLC tumor tissue. Additionally, high PEX5L expression levels were significantly associated with lower survival of NSCLC. Our data suggest that variants in these peroxisome-related genes may influence gene regulation and are potential predictors of NSCLC OS, once validated by additional studies.

Beyond rare disorders: A new era for peroxisomal pathophysiology

Metabolism is emerging as a central influencer of multiple disease states in humans. Peroxisomes are central metabolic organelles whose decreased function gives rise to severe peroxisomal diseases. Recently, it is becoming clear that, beyond such rare inborn errors, the deterioration of peroxisomal functions contributes to multiple and prevalent diseases such as cancer, viral infection, diabetes, and neurodegeneration. Despite the clear importance of peroxisomes in common pathophysiological processes, research on the mechanisms underlying their contributions is still sparse. Here, we highlight the timeliness of focusing on peroxisomes in current research on central, abundant, and society-impacting human pathologies. As peroxisomes are now coming into the spotlight, it is clear that intensive research into these important organelles will enable a better understanding of their contribution to human health, serving as the basis to develop new diagnostic and therapeutic approaches to prevent and treat human diseases.

Mice with a deficiency in Peroxisomal Membrane Protein 4 (PXMP4) display mild changes in hepatic lipid metabolism

Peroxisomes play an important role in the metabolism of a variety of biomolecules, including lipids and bile acids. Peroxisomal Membrane Protein 4 (PXMP4) is a ubiquitously expressed peroxisomal membrane protein that is transcriptionally regulated by peroxisome proliferator-activated receptor α (PPARα), but its function is still unknown. To investigate the physiological function of PXMP4, we generated a Pxmp4 knockout (Pxmp4^-/-) mouse model using CRISPR/Cas9-mediated gene editing. Peroxisome function was studied under standard chow-fed conditions and after stimulation of peroxisomal activity using the PPARα ligand fenofibrate or by using phytol, a metabolite of chlorophyll that undergoes peroxisomal oxidation. Pxmp4^-/- mice were viable, fertile, and displayed no changes in peroxisome numbers or morphology under standard conditions. Also, no differences were observed in the plasma levels of products from major peroxisomal pathways, including very long-chain fatty acids (VLCFAs), bile acids (BAs), and BA intermediates di- and trihydroxycholestanoic acid. Although elevated levels of the phytol metabolites phytanic and pristanic acid in Pxmp4^-/- mice pointed towards an impairment in peroxisomal α-oxidation capacity, treatment of Pxmp4^-/- mice with a phytol-enriched diet did not further increase phytanic/pristanic acid levels. Finally, lipidomic analysis revealed that loss of Pxmp4 decreased hepatic levels of the alkyldiacylglycerol class of neutral ether lipids, particularly those containing polyunsaturated fatty acids. Together, our data show that while PXMP4 is not critical for overall peroxisome function under the conditions tested, it may have a role in the metabolism of (ether)lipids.

rs9390123 and rs9399451 influence the DNA repair capacity of lung cancer by regulating PEX3 and PHACTR2‑AS1 expression instead of PHACTR2

Lung cancer is a common cancer type, and has the highest mortality rate in the world. A genome‑wide association study suggests that the genetic marker rs9390123 is significantly associated with DNA repair capacity (DRC) in lung cancer. Analysis of the data derived from the 1000 Genomes Project indicated that there is another single nucleotide polymorphism (SNP), rs9399451, in strong linkage disequilibrium with rs9390123 in Caucasian individuals, thus suggesting that this SNP could be associated with DRC. However, the causal SNP and mechanism of DRC remain unclear. In the present study, dual luciferase assay results indicated that both SNPs are functional in lung cells. Through chromosome conformation capture, an enhancer containing the two functional SNPs was observed to bind the promoter of peroxisomal biogenesis factor 3 and phosphatase and actin regulator 2 antisense RNA 1 (PHACTR2‑AS1). Knockdown of PHACTR2‑AS1 could significantly influence lung cell proliferation, colony formation, migration and wound healing, which verified that PHACTR2‑AS1 is a novel oncogene for lung cancer. Through chromatin immunoprecipitation, the transcription factor POU class 2 homeobox 1 (POU2F1) was identified to bind to the surrounding segments of these two SNPs, and their interaction was investigated. The present study identified the mechanism via which rs9390123 and rs9399451 could influence DRC.

Proteomic Analysis of Mouse Kidney Tissue Associates Peroxisomal Dysfunction with Early Diabetic Kidney Disease

Background: The absence of efficient inhibitors for diabetic kidney disease (DKD) progression reflects the gaps in our understanding of DKD molecular pathogenesis. Methods: A comprehensive proteomic analysis was performed on the glomeruli and kidney cortex of diabetic mice with the subsequent validation of findings in human biopsies and omics datasets, aiming to better understand the underlying molecular biology of early DKD development and progression. Results: LC–MS/MS was employed to analyze the kidney proteome of 2 DKD models: Ins2Akita (early and late DKD) and db/db mice (late DKD). The abundance of detected proteins was defined. Pathway analysis of differentially expressed proteins in the early and late DKD versus the respective controls predicted dysregulation in DKD hallmarks (peroxisomal lipid metabolism and β-oxidation), supporting the functional relevance of the findings. Comparing the observed protein changes in early and late DKD, the consistent upregulation of 21 and downregulation of 18 proteins was detected. Among these were downregulated peroxisomal and upregulated mitochondrial proteins. Tissue sections from 16 DKD patients were analyzed by IHC confirming our results. Conclusion: Our study shows an extensive differential expression of peroxisomal proteins in the early stages of DKD that persists regardless of the disease severity, providing new perspectives and potential markers of diabetic kidney dysfunction.

The Peroxisome-Autophagy Redox Connection: A Double-Edged Sword?

Peroxisomes harbor numerous enzymes that can produce or degrade hydrogen peroxide (H2O2). Depending on its local concentration and environment, this oxidant can function as a redox signaling molecule or cause stochastic oxidative damage. Currently, it is well-accepted that dysfunctional peroxisomes are selectively removed by the autophagy-lysosome pathway. This process, known as "pexophagy," may serve a protective role in curbing peroxisome-derived oxidative stress. Peroxisomes also have the intrinsic ability to mediate and modulate H2O2-driven processes, including (selective) autophagy. However, the molecular mechanisms underlying these phenomena are multifaceted and have only recently begun to receive the attention they deserve. This review provides a comprehensive overview of what is known about the bidirectional relationship between peroxisomal H2O2 metabolism and (selective) autophagy. After introducing the general concepts of (selective) autophagy, we critically examine the emerging roles of H2O2 as one of the key modulators of the lysosome-dependent catabolic program. In addition, we explore possible relationships among peroxisome functioning, cellular H2O2 levels, and autophagic signaling in health and disease. Finally, we highlight the most important challenges that need to be tackled to understand how alterations in peroxisomal H2O2 metabolism contribute to autophagy-related disorders.

Antioxidant and cytoprotective effects of avocado oil and extract (Persea americana Mill) against rotenone using monkey kidney epithelial cells (Vero)

Oxidative stress is known to be involved in development of numerous diseases including cardiovascular, respiratory, renal, kidney and cancer. Thus, investigations that mimic oxidative stress in vitro may play an important role to find new strategies to control oxidative stress and subsequent consequences are important. Rotenone, widely used as a pesticide has been used as a model to simulate oxidative stress. However, this chemical was found to produce several diseases. Therefore, the aim of this study was to investigate the antioxidant and cytoprotective effect of avocado (Persea americana Mill) extract and oil in monkey kidney epithelial cells (VERO) exposed to rotenone. VERO cells were exposed to IC₅₀ of rotenone in conjunction with different concentrations of avocado extract and oil (ranging from 1 to 1000 µg/ml), for 24 hr. Subsequently, cell viability and oxidative metabolism were assessed. Data demonstrated that avocado extract and oil in the presence of rotenone increased cellular viability at all tested concentrations compared to cells exposed only to rotenone. In addition, extract and avocado oil exhibited antioxidant action as evidenced by decreased levels of reactive oxygen species (ROS), superoxide ion, and lipid peroxidation, generated by rotenone. Further, avocado extract and oil appeared to be safe, since these compounds did not affect cell viability and or generate oxidative stress. Therefore, avocado appears to display a promising antioxidant potential by decreasing oxidative stress.

Downregulation of ACE2 is associated with advanced pathological features and poor prognosis in clear cell renal cell carcinoma

Aims: The aim of this study was to explore the alteration in ACE2 expression and correlation between ACE2 expression and immune infiltration in clear cell renal cell carcinoma (ccRCC). Methods: The authors first analyzed the expression profiles and prognostic value of ACE2 in ccRCC patients using The Cancer Genome Atlas public database. The authors used ESTIMATE and CIBERSORT algorithms to analyze the correlation between ACE2 expression and tumor microenvironment in ccRCC samples. Results: ACE2 was correlated with sex, distant metastasis, clinical stage, tumor T stage and histological grade. Moreover, downregulation of ACE2 was correlated with unfavorable prognosis. In addition, ACE2 expression was associated with different immune cell subtypes. Conclusion: The authors' analyses suggest that ACE2 plays an important role in the development and progression of ccRCC and may serve as a potential prognostic biomarker in ccRCC patients.

Potential lncRNA Biomarkers for HBV-Related Hepatocellular Carcinoma Diagnosis Revealed by Analysis on Coexpression Network

Background

Increasing evidence demonstrated that long noncoding RNA (lncRNA) could affect inflammatory tumor immune microenvironment by modulating gene expression and could be used as a biomarker for HBC-related hepatocellular carcinoma (HCC) but still needs further research. The aim of the present study was to determine an lncRNA signature for the diagnosis of HBV-related HCC.

Methods

HBV-related HCC expression profiles (GSE55092, GSE19665, and GSE84402) were abstracted from the GEO (Gene Expression Omnibus) data resource, and R package limma and RobustRankAggreg were employed to identify common differentially expressed genes (DEGs). Using machine learning, optimal diagnostic lncRNA molecular markers for HBV-related HCC were identified. The expression of candidate lncRNAs was cross-validated in GSE121248, and an ROC (receiver operating characteristic) curve of lncRNA biomarkers was carried out. Additionally, a coexpression network and functional annotation was built, after which a PPI (protein-protein interaction) network along with module analysis were conducted with the Cytoscape open source software.

Result

A total of 38 DElncRNAs and 543 DEmRNAs were identified with a fold change larger than 2.0 and a P value < 0.05. By machine learning, AL356056.2, AL445524.1, TRIM52-AS1, AC093642.1, EHMT2-AS1, AC003991.1, AC008040.1, LINC00844, and LINC01018 were screened out as optional diagnostic lncRNA biosignatures for HBV-related HCC. The AUC (areas under the curve) of the SVM (support vector machine) model and random forest model were 0.957 and 0.904, respectively, and the specificity and sensitivity were 95.7 and 100% and 94.3 and 86.5%, respectively. The results of functional enrichment analysis showed that the integrated coexpressed DEmRNAs shared common cascades in the p53 signaling pathway, retinol metabolism, PI3K-Akt signaling cascade, and chemical carcinogenesis. The integrated DEmRNA PPI network complex was found to be comprised of 87 nodes, and two vital modules with a high degree were selected with the MCODE app.

Conclusion

The present study identified nine potential diagnostic biomarkers for HBV-related HCC, all of which could potentially modulated gene expression related to inflammatory conditions in the tumor immune microenvironment. The functional annotation of the target DEmRNAs yielded novel evidence in evaluating the precise functions of lncRNA in HBV-related HCC.

Differential Expression of Peroxisomal Proteins in Distinct Types of Parotid Gland Tumors

Salivary gland cancers are rare but aggressive tumors that have poor prognosis and lack effective cure. Of those, parotid tumors constitute the majority. Functioning as metabolic machinery contributing to cellular redox balance, peroxisomes have emerged as crucial players in tumorigenesis. Studies on murine and human cells have examined the role of peroxisomes in carcinogenesis with conflicting results. These studies either examined the consequences of altered peroxisomal proliferators or compared their expression in healthy and neoplastic tissues. None, however, examined such differences exclusively in human parotid tissue or extended comparison to peroxisomal proteins and their associated gene expressions. Therefore, we examined differences in peroxisomal dynamics in parotid tumors of different morphologies. Using immunofluorescence and quantitative PCR, we compared the expression levels of key peroxisomal enzymes and proliferators in healthy and neoplastic parotid tissue samples. Three parotid tumor subtypes were examined: pleomorphic adenoma, mucoepidermoid carcinoma and acinic cell carcinoma. We observed higher expression of peroxisomal matrix proteins in neoplastic samples with exceptional down regulation of certain enzymes; however, the degree of expression varied between tumor subtypes. Our findings confirm previous experimental results on other organ tissues and suggest peroxisomes as possible therapeutic targets or markers in all or certain subtypes of parotid neoplasms.

Silencing PEX26 as an unconventional mode to kill drug-resistant cancer cells and forestall drug resistance

Promoting the macroautophagy/autophagy-mediated degradation of specific proteins and organelles can potentially be utilized to induce apoptosis in cancer cells or sensitize tumor cells to therapy. To examine this concept, we enriched for autophagosomes from histone deacetylase inhibitor (HDACi)-sensitive U937 lymphoma cells and isogenic HDACi-resistant cells. Mass spectrometry on autophagosome-enriched fractions revealed that HDACi-resistant cells undergo elevated pexophagy, or autophagy of the peroxisome, an organelle that supports tumor growth. To disturb peroxisome homeostasis, we enhanced pexophagy in HDACi-resistant cells via genetic silencing of peroxisome exportomer complex components (PEX1, PEX6, or PEX26). This consequently sensitized resistant cells to HDACi-mediated apoptosis, which was rescued by inhibiting ATM/ataxia-telangiectasia mutated (ATM serine/threonine kinase), a mediator of pexophagy. We subsequently engineered melanoma cells to stably repress PEX26 using CRISPR interference (CRISPRi). Melanoma cells with repressed PEX26 expression showed evidence of both increased pexophagy and peroxisomal matrix protein import defects versus single guide scrambled (sgSCR) controls. In vivo studies showed that sgPEX26 melanoma xenografts recurred less compared to sgSCR xenografts, following the development of resistance to mitogen-activated protein kinase (MAPK)-targeted therapy. Finally, prognostic analysis of publicly available datasets showed that low expression levels of PEX26, PEX6 and MTOR, were significantly associated with prolonged patient survival in lymphoma, lung cancer and melanoma cohorts. Our work highlighted that drugs designed to disrupt peroxisome homeostasis may serve as unconventional therapies to combat therapy resistance in cancer.Abbreviations: ABCD3/PMP70: ATP binding cassette subfamily D member 3; ACOX1: acyl-CoA oxidase 1; AP: autophagosome; COX: cytochrome c oxidase; CQ: chloroquine; CRISPRi: clustered regularly interspaced short palindromic repeats interference; DLBCL: diffuse large B-cell lymphoma; GO: gene ontology; dCas9: Cas9 endonuclease dead, or dead Cas9; HDACi: histone deacetylase inhibitors; IHC: Immunohistochemistry; LAMP2: lysosomal associated membrane protein 2; LCFAs: long-chain fatty acids; LFQ-MS: label-free quantitation mass spectrometry; LPC: lysophoshatidylcholine; MAP1LC3B/LC3B: microtubule associated protein 1 light chain 3 beta; MTOR: mechanistic target of rapamycin kinase; PBD: peroxisome biogenesis disorders; PTS1: peroxisomal targeting signal 1; ROS: reactive oxygen species; sgRNA: single guide RNA; VLCFAs: very-long chain fatty acids; Vor: vorinostat; WO: wash-off.

PI5P4Ks drive metabolic homeostasis through peroxisome-mitochondria interplay

PI5P4Ks are a class of phosphoinositide kinases that phosphorylate PI-5-P to PI-4,5-P₂. Distinct localization of phosphoinositides is fundamental for a multitude of cellular functions. Here, we identify a role for peroxisomal PI-4,5-P₂ generated by the PI5P4Ks in maintaining energy balance. We demonstrate that PI-4,5-P₂ regulates peroxisomal fatty acid oxidation by mediating trafficking of lipid droplets to peroxisomes, which is essential for sustaining mitochondrial metabolism. Using fluorescent-tagged lipids and metabolite tracing, we show that loss of the PI5P4Ks significantly impairs lipid uptake and β-oxidation in the mitochondria. Further, loss of PI5P4Ks results in dramatic alterations in mitochondrial structural and functional integrity, which under nutrient deprivation is further exacerbated, causing cell death. Notably, inhibition of the PI5P4Ks in cancer cells and mouse tumor models leads to decreased cell viability and tumor growth, respectively. Together, these studies reveal an unexplored role for PI5P4Ks in preserving metabolic homeostasis, which is necessary for tumorigenesis.

Mechanisms and Functions of Pexophagy in Mammalian Cells

Peroxisomes play essential roles in diverse cellular metabolism functions, and their dynamic homeostasis is maintained through the coordination of peroxisome biogenesis and turnover. Pexophagy, selective autophagic degradation of peroxisomes, is a major mechanism for removing damaged and/or superfluous peroxisomes. Dysregulation of pexophagy impairs the physiological functions of peroxisomes and contributes to the progression of many human diseases. However, the mechanisms and functions of pexophagy in mammalian cells remain largely unknown compared to those in yeast. This review focuses on mammalian pexophagy and aims to advance the understanding of the roles of pexophagy in human health and diseases. Increasing evidence shows that ubiquitination can serve as a signal for pexophagy, and ubiquitin-binding receptors, substrates, and E3 ligases/deubiquitinases involved in pexophagy have been described. Alternatively, pexophagy can be achieved in a ubiquitin-independent manner. We discuss the mechanisms of these ubiquitin-dependent and ubiquitin-independent pexophagy pathways and summarize several inducible conditions currently used to study pexophagy. We highlight several roles of pexophagy in human health and how its dysregulation may contribute to diseases.

Perspectives on Organelle Interaction, Protein Dysregulation, and Cancer Disease

In recent decades, compelling evidence has emerged showing that organelles are not static structures but rather form a highly dynamic cellular network and exchange information through membrane contact sites. Although high-throughput techniques facilitate identification of novel contact sites (e.g., organelle-organelle and organelle-vesicle interactions), little is known about their impact on cellular physiology. Moreover, even less is known about how the dysregulation of these structures impacts on cellular function and therefore, disease. Particularly, cancer cells display altered signaling pathways involving several cell organelles; however, the relevance of interorganelle communication in oncogenesis and/or cancer progression remains largely unknown. This review will focus on organelle contacts relevant to cancer pathogenesis. We will highlight specific proteins and protein families residing in these organelle-interfaces that are known to be involved in cancer-related processes. First, we will review the relevance of endoplasmic reticulum (ER)-mitochondria interactions. This section will focus on mitochondria-associated membranes (MAMs) and particularly the tethering proteins at the ER-mitochondria interphase, as well as their role in cancer disease progression. Subsequently, the role of Ca²⁺ at the ER-mitochondria interphase in cancer disease progression will be discussed. Members of the Bcl-2 protein family, key regulators of cell death, also modulate Ca²⁺ transport pathways at the ER-mitochondria interphase. Furthermore, we will review the role of ER-mitochondria communication in the regulation of proteostasis, focusing on the ER stress sensor PERK (PRKR-like ER kinase), which exerts dual roles in cancer. Second, we will review the relevance of ER and mitochondria interactions with other organelles. This section will focus on peroxisome and lysosome organelle interactions and their impact on cancer disease progression. In this context, the peroxisome biogenesis factor (PEX) gene family has been linked to cancer. Moreover, the autophagy-lysosome system is emerging as a driving force in the progression of numerous human cancers. Thus, we will summarize our current understanding of the role of each of these organelles and their communication, highlighting how alterations in organelle interfaces participate in cancer development and progression. A better understanding of specific organelle communication sites and their relevant proteins may help to identify potential pharmacological targets for novel therapies in cancer control.

Prognostic role of PHYH for overall survival (OS) in clear cell renal cell carcinoma (ccRCC)

This study attempts to evaluate the prognostic role of PHYH for overall survival (OS) in clear cell renal cell carcinoma (ccRCC) by means of publicly available data from The Cancer Genome Atlas (TCGA). Clinical pathologic features and PHYH expression were downloaded from the TCGA database and relationships between them were analyzed by univariate and multivariate Cox regression analyses. Gene Set Enrichment Analysis (GSEA) and gene–gene interactions were also performed between tissues with different PHYH expression levels. PHYH expression levels were significantly lower in patient with ccRCC compared with normal tissues (p = 1.156e−19). Kaplan–Meier survival analysis showed that high expression of PHYH had a better prognosis than low expression (p = 9e−05). Moreover, PHYH expression was also significantly associated with high grade (G2-4, p = 0.025), high stage (StageIII & IV, p = 5.604e−05), and high level of stage_T (T3-4, p = 4.373e−05). Univariate and multivariate Cox regression analyses indicated that PHYH could be acted as an independent prognostic factor (p < 0.05). Nomogram including clinical pathologic features and PHYH expression were also provided. GSEA revealed that butanoate metabolism, histidine metabolism, propanoate metabolism, pyruvate metabolism, tryptophan metabolism, PPAR signalling pathway, and renin–angiotensin system were differentially enriched in PHYH high-expression phenotype. ICGC database was utilized to verify the expression level and survival benefit of PHYH (both p < 0.05). We suspect that elevated PHYH expression may be served as a potential prognostic molecular marker of better survival in ccRCC. Besides, alpha-oxidation was closely regulated by PHYH, and PPAR signalling, pyruvate metabolism, butanoate metabolism, and RAS might be the key pathways regulated by PHYH in CCRC.

The diversity and breadth of cancer cell fatty acid metabolism

Tumor cellular metabolism exhibits distinguishing features that collectively enhance biomass synthesis while maintaining redox balance and cellular homeostasis. These attributes reflect the complex interactions between cell-intrinsic factors such as genomic-transcriptomic regulation and cell-extrinsic influences, including growth factor and nutrient availability. Alongside glucose and amino acid metabolism, fatty acid metabolism supports tumorigenesis and disease progression through a range of processes including membrane biosynthesis, energy storage and production, and generation of signaling intermediates. Here, we highlight the complexity of cellular fatty acid metabolism in cancer, the various inputs and outputs of the intracellular free fatty acid pool, and the numerous ways that these pathways influence disease behavior.

Peroxisomes form intralumenal vesicles with roles in fatty acid catabolism and protein compartmentalization in Arabidopsis

Peroxisomes are vital organelles that compartmentalize critical metabolic reactions, such as the breakdown of fats, in eukaryotic cells. Although peroxisomes typically are considered to consist of a single membrane enclosing a protein lumen, more complex peroxisomal membrane structure has occasionally been observed in yeast, mammals, and plants. However, technical challenges have limited the recognition and understanding of this complexity. Here we exploit the unusually large size of Arabidopsis peroxisomes to demonstrate that peroxisomes have extensive internal membranes. These internal vesicles accumulate over time, use ESCRT (endosomal sorting complexes required for transport) machinery for formation, and appear to derive from the outer peroxisomal membrane. Moreover, these vesicles can harbor distinct proteins and do not form normally when fatty acid β-oxidation, a core function of peroxisomes, is impaired. Our findings suggest a mechanism for lipid mobilization that circumvents challenges in processing insoluble metabolites. This revision of the classical view of peroxisomes as single-membrane organelles has implications for all aspects of peroxisome biogenesis and function and may help address fundamental questions in peroxisome evolution.

Peroxisomes are organelles compartmentalising metabolic reactions such as the breakdown of fats, and are commonly thought of as single membrane-bound compartments. Here the authors show that Arabidopsis peroxisomes contain extensive internal vesicles that form from the bounding membrane in an ESCRT-dependent process.

Peroxisome Metabolism in Cancer

Peroxisomes are metabolic organelles involved in lipid metabolism and cellular redox balance. Peroxisomal function is central to fatty acid oxidation, ether phospholipid synthesis, bile acid synthesis, and reactive oxygen species homeostasis. Human disorders caused by genetic mutations in peroxisome genes have led to extensive studies on peroxisome biology. Peroxisomal defects are linked to metabolic dysregulation in diverse human diseases, such as neurodegeneration and age-related disorders, revealing the significance of peroxisome metabolism in human health. Cancer is a disease with metabolic aberrations. Despite the critical role of peroxisomes in cell metabolism, the functional effects of peroxisomes in cancer are not as well recognized as those of other metabolic organelles, such as mitochondria. In addition, the significance of peroxisomes in cancer is less appreciated than it is in degenerative diseases. In this review, I summarize the metabolic pathways in peroxisomes and the dysregulation of peroxisome metabolism in cancer. In addition, I discuss the potential of inactivating peroxisomes to target cancer metabolism, which may pave the way for more effective cancer treatment.

Targeting the Redox Landscape in Cancer Therapy

Reactive oxygen species (ROS) are produced predominantly by the mitochondrial electron transport chain and by NADPH oxidases in peroxisomes and in the endoplasmic reticulum. The antioxidative defense counters overproduction of ROS with detoxifying enzymes and molecular scavengers, for instance, superoxide dismutase and glutathione, in order to restore redox homeostasis. Mutations in the redox landscape can induce carcinogenesis, whereas increased ROS production can perpetuate cancer development. Moreover, cancer cells can increase production of antioxidants, leading to resistance against chemo- or radiotherapy. Research has been developing pharmaceuticals to target the redox landscape in cancer. For instance, inhibition of key players in the redox landscape aims to modulate ROS production in order to prevent tumor development or to sensitize cancer cells in radiotherapy. Besides the redox landscape of a single cell, alternative strategies take aim at the multi-cellular level. Extracellular vesicles, such as exosomes, are crucial for the development of the hypoxic tumor microenvironment, and hence are explored as target and as drug delivery systems in cancer therapy. This review summarizes the current pharmaceutical and experimental interventions of the cancer redox landscape.

Computational Proteome-Wide Study for the Prediction of Escherichia coli Protein Targeting in Host Cell Organelles and Their Implication in Development of Colon Cancer

Enterohemorrhagic Escherichia coli infection is associated with gastrointestinal disorders, including diarrhea and colorectal cancer. Although evidences have established the involvement of E. coli in the growth of colon cancer, the molecular mechanisms of carcinogenesis of cancer growth and development are not well understood. We analyzed E. coli protein targeting in host cell organelles and the implication in colon cancer using in silico approaches. Our results indicated that many E. coli proteins targeted the endoplasmic reticulum (ER), ER membranes, Golgi apparatus, Golgi apparatus membranes, peroxisomes, nucleus, nuclear membrane, mitochondria, and mitochondrial membrane of host cells. These targeted proteins in ER, Golgi apparatus, peroxisomes, nucleus, and mitochondria may alter the normal functioning of various pathways including DNA repair, apoptosis, replication, transcription, and protein folding in E. coli-infected host cells. The results of the current in silico study provide insights into E. coli pathogenesis and may aid in designing new preventive and therapeutic strategies.

PEX5, a novel target of microRNA-31-5p, increases radioresistance in hepatocellular carcinoma by activating Wnt/β-catenin signaling and homologous recombination

Rationale: Hepatocellular carcinoma (HCC) is the second leading cause of cancer-related death worldwide, with high recurrence and metastasis rates. Although radiation is an effective treatment for tumors, it is often limited by intrinsic radioresistance in HCC. The contributions of dysregulated microRNAs, including miR-31-5p, to HCC progression have been recently reported. However, the role of miR-31-5p in the radiation response of HCC is unknown. In this study, we aimed to investigate the impact of miR-31-5p on HCC radiosensitivity.

Methods: miR-31-5p expression in HCC tissues, paired adjacent tissues, and HCC cell lines was measured using quantitative real-time polymerase chain reaction and in situ hybridization. Bioinformatic analyses, gain- and loss-of-function experiments, and luciferase reporter assays were performed to validate peroxisomal biogenesis factor 5 (PEX5) as a direct target of miR-31-5p. The biofunctions of PEX5 and miR-31-5p in HCC were determined by Transwell, wound-healing, and Cell Counting Kit-8 (CCK8) assays. A colony formation assay was used to evaluate the radiosensitivity of HCC cells. The interaction among PEX5, β-catenin, Rac1, and JNK-2 was confirmed by coimmunoprecipitation. A xenograft tumor model was established to validate the effects of miR-31-5p and PEX5 on HCC progression and radiosensitivity in vivo.

Results: Low expression of miR-31-5p in HCC specimens, as observed in this study, predicted a poor clinical outcome. However, the expression pattern of PEX5, as a direct target of miR-31-5p, was opposite that of miR-31-5p, and high PEX5 expression indicated poor prognosis in HCC patients. Ectopic expression of PEX5 increased the proliferation, migration, and invasion abilities and enhanced the radioresistance of HCC cells in vitro and in vivo; however, these phenotypes were inhibited by miR-31-5p. Mechanistically, PEX5 stabilized cytoplasmic β-catenin and facilitated β-catenin nuclear translocation to activate Wnt/β-catenin signaling. Moreover, upon radiation exposure, PEX5 reduced excessive reactive oxygen species (ROS) accumulation and activated the homologous recombination (HR) pathway, which protected HCC cells from radiation-induced damage.

Conclusions: Our findings demonstrated a novel role for PEX5 as a miR-31-5p target and a mediator of the Wnt/β-catenin signaling and HR pathways, providing new insights into studying HCC radiation responses and implicating PEX5 and miR-31-5p as potential therapeutic targets in HCC.

Transcriptomic profiling of peroxisome-related genes reveals a novel prognostic signature in hepatocellular carcinoma

Emerging evidence suggests that peroxisomes play a role in the regulation of tumorigenesis and cancer progression. However, the prognostic value of peroxisome-related genes has been rarely investigated. This study aimed to establish a peroxisome-related gene signature for overall survival (OS) prediction in patients with hepatocellular carcinoma (HCC). First, univariate Cox regression analysis was employed to identify prognostic peroxisome-related genes in The Cancer Genome Atlas liver cancer cohort, and least absolute shrinkage and selection operator Cox regression analysis was used to construct a 10-gene signature. The risk score based on the signature was positively correlated with poor prognosis (HR = 4.501, 95% CI = 3.021–6.705, P = 1.39e−13). Second, multivariate Cox regression incorporating additional characteristics revealed that the signature was an independent predictor. Time-dependent ROC curves demonstrated good performance of the signature in predicting the OS of HCC patients. The prognostic performance was validated using International Cancer Genome Consortium HCC cohort data. Gene set enrichment analysis revealed that the signature-related alterations in biological processes mainly involved peroxisomal functions. Finally, we developed a nomogram model based on the gene signature and TNM stage, which showed a superior prognostic power (C-index = 0.702). Thus, our study revealed a novel peroxisome-related gene signature that may help improve personalized OS prediction in HCC patients.

HSD17B4, ACAA1, and PXMP4 in Peroxisome Pathway Are Down-Regulated and Have Clinical Significance in Non-small Cell Lung Cancer

To explore the potential functions and clinical significances of peroxisomes during lung cancer development and progression, we investigated the expressional profiles of peroxisome pathway genes and their correlations with clinical features in non-small cell lung cancer (NSCLC). The RNA-seq data of NSCLC including lung squamous carcinoma (LUSC) and lung adenocarcinoma (LUAD) patients with their clinical information were downloaded from The Cancer Genome Atlas (TCGA). Gene expression comparisons between tumor and normal samples were performed with edgeR package in R software and the results of the 83 peroxisome pathway genes were extracted. Through Venn diagram analysis, 38 common differentially expressed peroxisome pathway genes (C-DEPGs) in NSCLC were identified. Principal components analysis (PCA) was performed and the 38 C-DEPGs could discriminate NSCLC tumors from the non-tumor controls well. Through Kaplan-Meier survival and Cox regression analyses, 11 of the C-DEPGs were shown to have prognostic effects on NSCLC overall survival (OS) and were considered as key C-DEPGs (K-DEPGs). Through Oncomine, Human Protein Atlas (HPA) and the Clinical Proteomic Tumor Analysis Consortium (CPTAC), three K-DEPGs (HSD17B4, ACAA1, and PXMP4) were confirmed to be down-regulated in NSCLC at both mRNA and protein level. Their dy-regulation mechanisms were revealed through their correlations with their copy number variations and methylation status. Their potential functions in NSCLC were explored through their NSCLC-specific co-expression network analysis, their correlations with immune infiltrations, immunomodulator gene expressions, MKI67 expression and their associations with anti-cancer drug sensitivity. Our findings suggested that HSD17B4, ACAA1, and PXMP4 might be new markers for NSCLC diagnosis and prognosis and might provide new clues for NSCLC treatment.

Peroxisomicine A1, a potential antineoplastic agent, causes micropexophagy in addition to macropexophagy

Peroxisomicine A1 (PA1) is a potential antineoplastic agent with high and selective toxicity toward peroxisomes of tumor cells. Pexophagy is a selective autophagy process that degrades damaged peroxisomes; this process has been studied mainly in methylotrophic yeasts. There are two main modes of pexophagy in yeast: macropexophagy and micropexophagy. Previous studies showed that peroxisomes damaged by a prolonged exposition to PA1 are eliminated by macropexophagy. In this work, Candida boidinii was grown in methanol-containing media, and PA1 was added to the cultures at 2 µg/mL after they reached the mid-exponential growth phase. Samples were taken at 5, 10, 15, 20, and 25 min after the addition of PA1 and processed for ultrastructural analysis. Typical morphological characteristics of micropexophagy were observed: the direct engulfment of peroxisomes by the vacuolar membrane and the presence of the micropexophagic membrane apparatus (MIPA), which mediates the fusion between the opposing tips of the vacuole to complete sequestration of peroxisomes from the cytosol. In conclusion, here we report that, in addition to macropexophagy, peroxisomes damaged by PA1 can be eliminated by micropexophagy. This information is useful to deepen the knowledge of the mechanism of action of PA1 and of that of pexophagy per se.