
==== Front
iScience
iScience
iScience
2589-0042
Elsevier

S2589-0042(24)01950-3
10.1016/j.isci.2024.110725
110725
Article
Targeting ALDOA to modulate tumorigenesis and energy metabolism in retinoblastoma
Wang Yinghao 12
Tang Junjie 12
Liu Yaoming 12
Zhang Zhihui 1
Zhang Hongwei 1
Ma Yujun 1
Wang Xinyue 1
Ai Siming 1
Mao Yuxiang 1
Zhang Ping 1
Chen Shuxia 1
Li Jinmiao 1
Gao Yang 1
Cheng Chao 1
Li Cheng 1
Su Shicai 1
Lu Rong lurong@gzzoc.com
13∗
1 State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center, Sun Yat-Sen University, Guangdong Provincial Key Laboratory of Ophthalmology and Visual Science, Guangzhou 510060, China
∗ Corresponding author lurong@gzzoc.com
2 These authors contributed equally

3 Lead contact

15 8 2024
20 9 2024
15 8 2024
27 9 11072529 11 2023
18 5 2024
9 8 2024
© 2024 The Author(s)
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).
Summary

This study aims to elucidate the pivotal role of aldolase A (ALDOA) in retinoblastoma (RB) and evaluate the potential of the ALDOA inhibitor itaconate in impeding RB progression. Utilizing single-cell RNA sequencing, ALDOA consistently exhibits overexpression across diverse cell types, particularly in cone precursor cells, retinoma-like cells, and retinoblastoma-like cells. This heightened expression is validated in RB tissues and cell lines. ALDOA knockdown significantly diminishes RB cell viability, impedes colony formation, and induces notable metabolic alterations. RNA-seq analysis identifies SUSD2, ARHGAP27, and CLK2 as downstream genes associated with ALDOA. The application of itaconate demonstrates efficacy in inhibiting RB cell proliferation, validated through in vitro and in vivo models. This study emphasizes ALDOA as a promising target for innovative RB therapies, with potential implications for altering tumor energy metabolism.

Graphical abstract

Highlights

• ALDOA overexpression in various retinoblastoma cell types and tissues

• ALDOA knockdown reduces retinoblastoma cell viability and alters metabolism

• Itaconate inhibits retinoblastoma proliferation in both in vitro and in vivo models

• Targeting ALDOA holds potential for advancing retinoblastoma therapies

Molecular biology; Cancer; Transcriptomics

Subject areas

Molecular biology
Cancer
Transcriptomics
Published: August 15, 2024
==== Body
pmcIntroduction

Retinoblastoma (RB) is the most prevalent ocular malignancy in infants and young children. It constitutes 2–4% of all pediatric tumors, wielding a profound impact on both the visual acuity and life expectancy of affected individuals.1 The incidence of RB typically ranges from 1 in 20,000 to 1 in 15,000, yielding approximately 8,000–9,000 new cases among children annually.2 Among the RB patients who experience metastasis, an alarming 80% eventually face a fatal outcome, with the invasive extraocular extension emerging as pivotal contributors to the metastatic process.3 Recent years have witnessed notable progress in the diagnosis and management of RB, leading to substantial improvements in both survival and rates of eye preservation. However, it is crucial to note that the advanced RB remain refractory to current therapeutic modalities, resulting in a less favorable prognosis.

The phenomenon, known as aerobic glycolysis or the Warburg effect, has evolved beyond a mere adaptation to the environment, now recognized as an integral component of malignant phenotypes and tumor metabolic reprogramming.4 Enhanced glycolysis in tumors not only fuels their proliferation but also provides intermediates for synthesizing critical biomolecules like nucleotides, lipids, and proteins.5 Consequently, targeting tumor glycolysis has arisen as a promising strategy in cancer therapy, as tumor cells adjust their energy metabolism to meet the demands of growth, biosynthesis, and redox reactions. Aldolase is a pivotal enzyme within the glycolysis pathway, facilitating the breakdown of 1,6-diphosphate-D-fructose into 3-phosphate-D-glyceraldehyde, and subsequently interacting with α-dihydroxyacetone phosphate.6 The aldolase family comprises three members: ALDOA, ALDOB, and ALDOC. Notably, ALDOA has been observed to be overexpressed in a variety of malignant tumors, including gastric cancer, hepatocellular carcinoma, colorectal cancer, and cervical cancer.7,8,9,10

Our recent research into the metabolic microenvironment of RB has revealed an upregulation of the glycolytic pathway.11 In this study, we aim to investigate the expression of key glycolytic enzymes in RB at the single-cell level. By integrating single-cell RNA sequencing (scRNA-seq) datasets from RB and normal retinal tissues, we compared gene expression profiles and confirmed the upregulation of ALDOA in RB. Subsequent in vitro and in vivo experiments demonstrated the effectiveness of targeted ALDOA inhibition in suppressing RB growth.

Results

Elevated ALDOA expression in retinoblastoma

GSEA of DEGs between retina and RB revealed significant enrichment in gene sets related to amide biosynthesis, glycolysis, catalytic activity/metabolism, and ATP synthesis (Figure 1A). This suggests a significant alteration in energy metabolism in RB, highlighting the potential crucial role of the glycolytic pathway in RB. Expression patterns of key genes within the glycolytic pathway was showed in heatmap (Figure 1B). Through dimensionality reduction and subsequent clustering analysis utilizing Uniform Manifold Approximation and Projection (UMAP), we visualized the expression profiles of Aldolase A (ALDOA) in RB and normal retina (Figure 2A). The analysis unveiled a high expression level of ALDOA in RB (Figure 2B), particularly elevated in extraocular RB compared to intraocular cases (Figure 2C). Branched expression analysis modeling unveiled a clear upregulation of ALDOA expression in retinoma-like cells as they progressed through pseudotime. However, this unique expression profile was not evident in other cell types (Figure S1). To validate ALDOA expression in RB, we examined its levels in human RB tissue samples and RB cell lines (WERI-RB1, Y79). The results demonstrated upregulated ALDOA mRNA and protein levels in RB cells compared to ARPE-19 (Figures 3A and 3B). Additionally, we analyzed ALDOA mRNA expression using the GSE125903 dataset12 from the Gene Expression Omnibus (GEO) database, distinguishing between invasive and non-invasive RB cases. The analysis revealed higher ALDOA expression in invasive RB cases compared to non-invasive ones (Figure 3C). Proteins extracted from tissues of 5 intraocular and 5 extraocular RB patients revealed higher ALDOA protein expression in extraocular RB, as detected by western blotting (Figure 3D). IHC results further confirmed significantly higher ALDOA expression in extraocular RB compared to intraocular RB and normal retina (Figure 3E). Clinicopathological analysis revealed a trend toward higher ALDOA expression in cases with high-risk features like scleral and optic nerve invasion. Notably, ALDOA expression was significantly lower in younger patients (<2 years old) (Figure S2).Figure 1 Enrichment analysis and glycolytic pathway expression in retinoblastoma

(A) Gene set enrichment analysis (GSEA) of differentially expressed genes (DEGs) between normal retina and retinoblastoma (RB) samples.

(B) Heatmap representation of the expression patterns of key genes within the glycolytic pathway, highlighting their differential expression in RB compared to normal retina.

Figure 2 Visualization of ALDOA expression profiles in retinoblastoma and normal retina

(A) Visualization of the expression profiles of Aldolase A (ALDOA) in retinoblastoma (RB) and normal retina using uniform manifold approximation and projection (UMAP).

(B) Comparison of ALDOA expression between RB and normal retina.

(C) Comparison of ALDOA expression between extraocular RB and intraocular case.

(D) ALDOA expression across various cell types.

Figure 3 ALDOA overexpression in various retinoblastoma cell types and tissues

(A and B) Quantitative PCR (A) and western blot assays (B) revealed significantly higher messenger RNA (mRNA) and protein levels of ALDOA in retinoblastoma (RB) cells compared to normal cells.

(C) ALDOA mRNA expression levels in non-invasive and invasive retinoblastoma cases from the GSE125903 dataset.

(D) Western blot analysis of ALDOA protein expression in tissues obtained from 5 intraocular and 5 extraocular RB patients, demonstrating elevated ALDOA protein levels in extraocular RB.

(E) Immunohistochemistry (IHC) staining and a statistical analysis of ALDOA expression in normal retina, intraocular RB tissue, and extraocular RB tissue. R, retina; RB, retinoblastoma; S, sclera; VC, vitreous cavity; ON, optic nerve. Data are presented as means ± standard deviations. Statistical significance was determined using Student’s t test, ∗p < 0.05, ∗∗p < 0.001, ∗∗∗p < 0.001, ∗∗∗∗p < 0.0001.

Disruption of retinoblastoma cell growth and energy metabolism by ALDOA downregulation

Given the observed ALDOA overexpression in RB, we investigated the functional consequences of ALDOA downregulation in RB cell lines. The qRT-PCR and western blotting results showed a decrease in ALDOA mRNA and protein expression levels after transfection (Figures 4A and 4B). The CCK-8 assay indicated that upon ALDOA knockdown in Y79 and WERI-RB1 cells, noticeable alterations in cell viability were observed from the second day, with more pronounced differences becoming evident by the third day when compared to the si-NC groups (Figure 4C). Additionally, the EdU assay results showed a decrease in the proliferation ability of Y79 and WERI-RB1 cells following ALDOA knockdown (Figure S3). The colony formation assays revealed a significant decrease in RB cells colony formation abilities following the knockdown of ALDOA (Figure 4D). Lactate assay results demonstrated a marked reduction in extracellular lactate level in RB cells following ALDOA knockdown (Figure 4E), accompanied by elevated levels of reactive oxygen species and decreased ATP levels (Figures 4F and 4G).Figure 4 Functional implications of ALDOA downregulation in retinoblastoma cell lines

(A and B) Quantitative reverse transcription-polymerase chain reaction (qRT-PCR) (A) and western blotting analysis (B) revealed a substantial reduction in ALDOA mRNA and protein expression levels following transfection of RB cell lines.

(C) The CCK-8 assay demonstrated inhibited cell proliferation in response to ALDOA downregulation.

(D) Colony formation experiments indicated a reduced tumorigenic potential of RB cell lines after ALDOA knockdown.

(E) Extracellular lactate analysis exhibited a notable decrease in lactate levels in RB cells treated with siALDOA for 48 h compared to the negative control.

(F) Levels of reactive oxygen species (ROS) in RB cells treated with siALDOA for 48 h were increased compared to the negative control.

(G) Intracellular ATP levels in RB cells treated with siALDOA for 48 h were reduced in comparison to the negative control. Data are presented as means ± standard deviations. Statistical significance was determined using Student’s t test, with ∗p < 0.05, ∗∗p < 0.001, and ∗∗∗p < 0.001 denoting significance.

Signaling alteration in response to ALDOA downregulation in retinoblastoma

Next, we performed transcriptome sequencing on both the si-ALDOA and si-NC WERI-RB1 to explore the downstream alterations resulting from ALDOA downregulation. Our transcriptome analysis revealed significant changes in gene expression patterns, highlighting 162 upregulated genes including PHRF1, SUSD2, IGFN1, UNC79, and 136 downregulated genes including PTP4A3, ARHGAP27, MYT1, CLK2 (Figures 5A and 5B). Further pathway enrichment analysis of these differentially expressed genes (DEGs) indicated potential shifts in biological functions in RB following ALDOA downregulation. Significantly, the upregulated DEGs were notably enriched in pathways related to “Organic acid transport”, “Amino acid transport”, and “Carboxylic acid transport” (Figure 5C). Conversely, the downregulated DEGs were predominantly associated with pathways involving “DNA-binding transcription activator activity” and “histone methyltransferase activity” (Figure 5D). Gene set enrichment analysis (GSEA) further revealed that after ALDOA knockdown, upregulated DEGs were enriched in pathways such as glycosphingolipid biosynthesis, while downregulated DEGs were enriched in metabolic pathways like oxidative phosphorylation, glutathione metabolism, and fatty acid metabolism (Figure S4). To validate the expression levels of the identified DEGs, we conducted qRT-PCR on selected genes displaying significant expression differences. The results demonstrated a substantial upregulation of SUSD2, ARHGAP27, and CLK2 mRNA levels in RB cell lines compared to the ARPE-19 (Figure 5E). Furthermore, in RB cell lines where ALDOA was knocked down, these genes exhibited altered mRNA expression patterns, with increased expression of SUSD2 and decreased expression of ARHGAP27 and CLK2 (Figure 5F). Western blot results showed that after ALDOA knockdown, the protein level of SUSD2 increased while the levels of CLK2 and ARHGAP27 decreased in WERI-RB1 cells. However, in Y79 cell lines, an increase in SUSD2 protein level was not detected after ALDOA knockdown, and the trends for the other two proteins were also decreasing (Figure 5G). Furthermore, to individually knock down SUSD2, ARHGAP27, and CLK2, three pairs of siRNAs were designed for each gene (Figure 5H). CCK8 experiments revealed that knocking down ARHGAP27 led to a decrease in the proliferation capacity of RB cells (Figure 5I).Figure 5 Effects of ALDOA downregulation on WERI-RB1 cell line

(A) RNA-seq analysis was employed to examine alterations in the mRNA expression profile of WERI-RB1 cells following ALDOA downregulation.

(B) Heatmap analysis illustrating differentially expressed genes (DEGs) in WERI-RB1 cells subjected to ALDOA knockdown.

(C and D) Pathway enrichment analysis focused on the upregulated DEGs (C) and downregulated DEGs (D). Upregulated DEGs were significantly associated with pathways related to “Organic acid transport”, “Amino acid transport”, and “Carboxylic acid transport”, while downregulated DEGs were primarily linked to “DNA-binding transcription activator activity” and “Histone methyltransferase activity”.

(E) Quantitative PCR analysis of selected DEGs (SUSD2, ARHGAP27, and CLK2).

(F) Quantitative PCR analysis of SUSD2, ARHGAP27, and CLK2 in RB cell lines with ALDOA knockdown.

(G) Western blot analysis showing changes in SUSD2, ARHGAP27, and CLK2 protein levels following ALDOA knockdown.

(H) Western blot analysis of SUSD2, ARHGAP27, and CLK2 protein knockdown.

(I) Proliferation changes in WERI-RB1 cells following knockdown of SUSD2, ARHGAP27, and CLK2. Statistical significance was assessed using Student’s t test, with ∗p < 0.05 and ∗∗p < 0.001 indicating significance.

Itaconate suppresses retinoblastoma cell proliferation and inhibits WERI-RB1 xenograft growth

Itaconate is an enzymatic inhibitor of ALDOA, which suppresses the catalytic activity of ALDOA in the glycolytic pathway through covalent modification of Cys73 and Cys339, without altering the protein expression level of ALDOA.13 Hence, we utilized itaconate as a means to inhibit the enzymatic activity of ALDOA, while keeping its protein expression level unaltered (Figure 6A). Results from the CCK8 assay demonstrate a reduction in RB cell viability upon the addition of itaconate (Figure 6B). Moreover, the introduction of itaconate led to alterations in energy metabolism within RB cells. Mass spectrometry analysis revealed alterations in cellular metabolite levels upon treatment with itaconate in RB cell lines. Specifically, there was a downward trend in the levels of metabolites involved in glycolysis pathways such as 3-Phospho-D-glycerate, Beta-D-Fructose 6 phosphate, and D-Glucose-6-phosphate, although not statistically significant (Figure 6F). However, metabolites associated with oxidative phosphorylation pathways including L-Malate acid, Fumarate, and Oxaloacetate exhibited a significant decrease (Figure 6G). To further investigate, an RB murine model was established by injecting WERI-RB1 cells into the eyes of nude mice. 14 days post-grafting, the animals were randomized into groups and received intraocular administrations of either PBS, itaconate, or melphalan, with enucleation performed on day 40 (Figures 7A and 7B). Analysis of eye weight revealed that both the itaconate and melphalan treatment groups had lower eye weights compared to the control group, although these differences were not statistically significant (Figure 7C). External examination of the eyes and H&E staining showed significant corneal neovascularization and extraocular invasion in the PBS group, with extensive tumor cell infiltration in the vitreous body and anterior chamber. In contrast, the itaconate treatment group exhibited tumor cells primarily located anterior to the retina. The melphalan-treated mice had fewer intraocular tumor cells but showed signs of ophthalmatrophia and retinal structural abnormalities (Figure 7D). Histological analysis revealed ophthalmatrophia in 40% of the melphalan-treated eyes and 20% of the itaconate-treated eyes, while extraocular invasion was observed in 30% of the PBS group and 10% of the itaconate group (Figure 7E). TUNEL staining showed a higher incidence of TUNEL-positive cells in the tumors of the itaconate group compared to PBS controls (Figure 7F). Additionally, the itaconate group had a lower proportion of TUNEL-positive cells in the retina compared to the melphalan-treated group (Figure 7G). These observations indicate that itaconate is effective in targeting intraocular RB cells similarly to melphalan, while maintaining a more favorable safety profile for retinal tissue.Figure 6 Itaconate-mediated inhibition of ALDOA and its impact on retinoblastoma in vitro and in vivo

(A) Western blotting analysis Itaconate serves as an enzymatic inhibitor of ALDOA while preserving ALDOA’s protein expression level.

(B) CCK8 assay reveal a reduction in RB cell viability following the addition of itaconate.

(C, D, and E) Introduction of itaconate induces alterations in energy metabolism within RB cells, leading to reduced lactate production (C), increased generation of reactive oxygen species (ROS) (D), and decreased ATP generation (E).

(F) Metabolite changes in the glycolytic pathway in RB cells treated with itaconate.

(G) Metabolite changes in the oxidative phosphorylation pathway in RB cells treated with itaconate. Statistical significance was assessed using Student’s t test, with ∗p < 0.05 and ∗∗p < 0.001 indicating significance.

Figure 7 Itaconate treatment in an RB murine model

(A) Study schematic illustrating the experimental design.

(B) Gross appearance of eyeballs 40 days post-xenograft.

(C) Eye weight analysis (n = 10).

(D) Representative views of enucleated eyes and H&E staining from each treatment group.

(E) Proportion of eyes with ophthalmatrophia and extraocular RB in the different treatment groups.

(F) Representative TUNEL images and quantitative analysis of intratumoral cells comparing the itaconate and melphalan groups.

(G) Representative TUNEL images and quantitative analysis of retinal cells comparing the itaconate and PBS groups. Data are presented as means ± SD. Statistical significance was determined by unpaired, two-tailed Student’s t test, with ∗p < 0.05 and ∗∗p < 0.01 indicating significance.

Discussion

In this study, we integrated the scRNA-seq datasets of normal retinal tissue and RB tissue, uncovering a significant upregulation of the glycolytic enzyme ALDOA in tumor tissue. Additionally, GSEA identified enrichment of gene sets containing ALDOA related to RB, particularly those associated with ATP generation, protein synthesis, and glycolysis. Next, we confirmed the elevated expression of ALDOA in both cell lines and patient tissues, and noted a significant increase in extraocular RB. Through the application of siRNA interference and the ALDOA inhibitor itaconate, we demonstrated that ALDOA inhibition alter tumor energy and metabolism and efficiently diminished the viability of RB cells and the tumor progression in the RB mouse model. Furthermore, by using transcriptome sequencing, we uncovered possible downstream targets under the control of ALDOA, such as SUSD2, ARHGAP27, and CLK2.

The mutation or loss of RB1 is a fundamental etiological factor in the development of RB. RB1 encodes the RB protein, which directs cellular metabolism reprogramming through E2F-induced transcriptional mechanisms, involving processes such as glucose oxidation phosphorylation, fatty acid oxidation, and amino acid synthesis.14 Nonetheless, the exact role of metabolic reprogramming in tumor progression and the specific metabolic enzymes implicated in the development of RB remain incompletely understood. ALDOA, a member of the aldolase enzyme family, has been the subject of extensive research, shedding light on its intricate interactions with hypoxia-inducible factor-1 (HIF-1) and its crucial role in modulating the glycolytic pathway and the process of epithelial-mesenchymal transition (EMT).15 It is worth noting that ALDOA expression levels have been thoroughly examined across various malignancies, revealing a direct link between ALDOA expression and the aggressiveness and invasiveness of tumors. Furthermore, comprehensive survival analyses consistently emphasize a strong positive correlation between increased ALDOA expression and an adverse prognosis.16 In line with previous findings, our research results validate a noticeable upregulation of ALDOA in RB tissue. Importantly, in cases of extraocular RB, ALDOA demonstrates higher expression levels compared to intraocular cases, highlighting the pivotal role of ALDOA in promoting the initiation and progression of RB. Recent research has increasingly focused on RB invasiveness, identifying targets such as UBE2C, SOX4, and the MCM family.17,18,19 However, additional investigation is needed to clarify the specific role of ALDOA in the extraocular invasion of RB. Liquid biopsies, especially aqueous humor biopsies, show promise for predicting RB prognosis and assessing chemotherapy efficacy.20,21,22,23 Investigations utilizing aqueous humor biopsy have shown potential in aiding RB prognosis and assessing chemotherapy efficacy.24,25 Thus, further exploration of tumor-related glycolytic enzymes in aqueous humor or blood samples may provide additional support for the diagnostic significance of ALDOA in extraocular RB.

In the course of RB development, distinct cell types play varying roles. Our analysis of scRNA-seq dataset reveals a pronounced upregulation of ALDOA in specific cell types, particularly cone cell precursors, RB-like cells, and retinoma like cells. Cone precursor like cells represent an early-stage cell population with the inherent potential for differentiation into cones or malignant cells.26 RB-like cells embody a highly proliferative malignant cell population characterized by a reduction in photoreceptor characteristics.27 Retinoma like cells serve as transitional cell stages between premalignant cone precursors and fully developed tumor cells.28 The specific high expression of ALDOA within these three cell types suggests its potential oncogenic role in RB. Apart from tumor cells, the tumor microenvironment comprises a complex network of cells, stroma, blood vessels, and immune cells surrounding the tumor.29,30 Within the tumor microenvironment, the expression levels of ALDOA may be regulated, and ALDOA can influence the metabolism and survival status of tumor cells. Specifically, ALDOA may influence the function and status of other cells within the tumor microenvironment by affecting the energy metabolism and signaling pathways of tumor cells.31,32 However, in this study, the impact of ALDOA on other cells within the RB microenvironment, such as immune cells, stromal cells, and endothelial cells, has not yet been investigated. Additionally, ALDOA may modulate factors such as oxygen concentration, pH balance, and nutrient supply within the tumor microenvironment, affecting the growth and chemoresistance of tumor cells.33,34,35 Moreover, previous studies have implicated the POU2F1-ALDOA axis in promoting tumor cell proliferation and chemoresistance.36 In this study, RNA sequencing revealed a significant alteration in POU2F1 expression following ALDOA knockdown, suggesting that the POU2F1-ALDOA axis may play a critical regulatory role in RB.

Itaconate has the capability to inhibit the catalytic activity of ALDOA in the glycolytic pathway without altering the expression level of the ALDOA protein. Initially discovered as a metabolite activated in macrophages, itaconate plays a critical role in the negative regulation of inflammatory responses by suppressing the production of inflammatory factors associated with macrophage activation.37 This regulatory process primarily functions through the alkylation of KEAP1, leading to the activation of NRF2.38 In the research of tumor, there has been a hypothesis that itaconate might indirectly influence tumor progression by its impact on macrophages.39 However, it remains uncertain whether itaconate directly regulates tumor progression. In this study, we found that although the intravitreal injection of itaconate did not result in a reduction in the rate of intraocular tumorigenesis, there was a notable decrease in the size of intraocular tumors, as well as a reduction in structural damage to the eyeball following the treatment. These observations suggest the potential of itaconate as an adjunctive therapy for RB. One of the primary challenges in translating preclinical findings into clinical applications is navigating the complex regulatory landscape. For itaconate to be considered for clinical use in treating RB or other tumors, detailed pharmacokinetic studies must be conducted. Additionally, thorough investigation of potential side effects is necessary to ensure patient safety.

Genes such as PHRF1, SUSD2, IGFN1, and UNC79 were identified as the top upregulated genes in ALDOA knockdown samples. SUSD2, situated on chromosome 22, encodes a transmembrane protein that plays a crucial role in cell-cell interactions and cell-matrix adhesion. Studies have indicated that SUSD2 acts as a tumor suppressor in RB by regulating the occurrence, development, and angiogenesis of RB, making it a potential target for RB treatment.40 Conversely, genes including PTP4A3, ARHGAP27, MYT1, and CLK2 were identified as the top downregulated genes in ALDOA knockdown samples. ARHGAP27 functions as a Rho GTPase-activating protein and is associated with pediatric chronic myeloid leukemia as well as ovarian cancer.41,42 CLK2, a phosphorylation-involved splicing protein kinase, exerts a promotional role in various cancers.43,44 These differentially expressed genes might play important roles in RB progression, part of which were validated in this study, while the intrinsic interaction network in RB still requires further exploration.

In conclusion, our study illuminates the important role of ALDOA in driving tumor growth in RB. These findings not only offer insights into the metabolic dynamics of RB but also present potential target candidates for the diagnosis and therapy of RB patients.

Limitations of the study

Our study has certain limitations that warrant discussion. First, the relatively small number of clinical samples used may constrain the analysis, potentially limiting the generalizability of our findings. Expanding the sample size and diversifying patient demographics would strengthen the robustness of our conclusions. Second, incorporating primary cells into future investigations could better capture the heterogeneity of RB, thus enhancing the validity of our results. Additionally, further exploration into patient outcomes and correlation with relevant clinical parameters would provide deeper insights into the clinical significance of our findings. Lastly, while our study successfully identifies downstream targets, a more comprehensive investigation is needed to elucidate the mechanistic intricacies of how ALDOA influences these targets and the broader signaling pathways involved.

Resource availability

Lead contact

Further information and requests for resources should be directed to and will be fulfilled by the lead contact, Rong Lu (lurong@gzzoc.com).

Materials availability

This study did not generate new unique reagents.

Data and code availability

• The raw sequence data reported in this paper have been deposited in the Genome Sequence Archive in National Genomics Data Center, China National Center for Bioinformation/Beijing Institute of Genomics, Chinese Academy of Sciences (GSA-Human: HRA008175) that are publicly accessible at https://ngdc.cncb.ac.cn/gsa-human. Original western blot images have been provided in the supplementary material.

• This paper does not report original code.

• Any additional information required to reanalyze the data reported in this paper is available from the lead contact upon request.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (82203175 ), the Natural Science Foundation of Guangdong Province, China (2022A1515011609 ), and the Science and Technology Program of Guangzhou, China (2024A03J0266 , 2024A03J0267 ).

Author contributions

Conceptualization, Y.W. and R.L.; Methodology, Y.W., J.T., Y.M., Z.Z., S.C., and Y.L.; Software, H.Z. and Y.L.; Investigation, J.L., Y.G., C.L., C.C., S.S., and R.L.; Sample collection: S.A., Y.M., P.Z., and R.L.; Visualization, J.T., H.Z., and Y.L.; Writing – original draft, J.T. and Y.L; Writing – review and editing, Y.W., X.W., J.L., and R.L.; Funding acquisition, Y.L. and R.L.; Supervision, R.L.

Declaration of interests

The authors declare no conflict of interest.

STAR★Methods

Key resources table

REAGENT or RESOURCE	SOURCE	IDENTIFIER	
Antibodies	
	
Recombinant Anti-SUSD2 antibody	Abcam	Cat# ab182147, RRID: AB_3082988	
Anti-Aldolase antibody	Abcam	Cat# ab252953, RRID: AB_3099487	
CLK2 antibody	Proteintech	Cat# 11092-1-AP, RRID: AB_2082254	
ARHGAP27 antibody	Proteintech	Cat# 22077-1-AP, RRID: AB_2878985	
Beta Tubulin antibody	Proteintech	Cat# 66240-1-Ig, RRID: AB_2881629	
Beta Actin Antibody	Proteintech	Cat# 66009-1-lg, RRID: AB_2782959	
	
Chemicals, peptides, and recombinant proteins	
	
RPMI 1640	Corning	#10-040-CV	
DMEM/F-12	Gibco	#11320033	
Fetal Bovine Serum	Gibco	#10099141C	
Penicillin-streptomycin	Gibco	#15070063	
HiPerFect Reagent	Qiagen	# 301704	
Melphalan	Selleck	#S8266	
Itaconic acid	Selleck	#S3095	
	
Critical commercial assays	
	
EZ-press RNA Purification Kit	EZBioscience	#B0004D	
4× Reverse Transcription Master Mix	EZBioscience	#A001GQ	
TUNEL Apoptosis Assay Kit	Beyotime	#C1088	
Cell Counting Kit-8	Vazyme	#A311-01	
BeyoClick™ EdU Cell Proliferation Kit with Alexa Fluor 488	Beyotime	# C0071S	
Lactate Colorimetric Assay Kit	Abbkine	#KTB1100	
ADP/ATP Ratio Assay Kit	Abnova	#KA1673	
	
Deposited data	
	
Raw data	Liu et al.18	GSE249995	
Raw data	Wu et al.27	GSE168434	
Raw data	Rajasekaran et al.12	GSE125903	
Raw data	Lukowski et al.45	E-MTAB-7316	
Raw data	This study	HRA008175	
	
Experimental models: Cell lines	
	
WERI-RB1	American Type Culture Collection	# HTB-169	
Y79	American Type Culture Collection	# HTB-18	
ARPE-19	American Type Culture Collection	# CRL-2302	
	
Experimental models: Organisms/strains	
	
Mouse: BALB/C nude	BesTest Bio-Tech	N/A	
	
Oligonucleotides	
	
Primers for Q-RT PCR, see Table S2	This paper	N/A	
	
Software and algorithms	
	
GraphPad Prism	GraphPad	V 8.0.2	
R (v 4.0.3)	R	r-project.org	
DESeq2	Love et al.46	https://bioconductor.org/packages/release/bioc/html/DESeq2.html	
GSEA	UC San Diego and Broad Institute	http://gsea-msigdb.org	

Experimental model and study participant details

Clinical samples

Experiments involving human retinoblastoma tissues were ethically approved by the Committee of Ethics at Zhongshan Ophthalmic Center in Guangzhou, China, ensuring compliance with the guidelines set forth in the Declaration of Helsinki (Ethical Approval Number: 2023KYPJ308). All patients included in the study provided informed consent. Samples from retinoblastoma and uveal melanoma patients were obtained from individuals undergoing surgery at the ocular tumor division of Zhongshan Ophthalmic Center. A section of uveal melanoma patient tissue was utilized as the control representing normal retina in the study. Extraocular RB refers to a tumor that extends beyond the eye, involving structures such as the orbit, optic nerve, meningeal spaces around the optic nerve, and full-thickness invasion of the sclera.47 It may also invade adjacent adipose tissue, extraocular muscles, bones, conjunctiva, or eyelids. Comprehensive clinical data for the patients involved were documented in Table S1. All patients included in this study were Chinese, of Asian ethnicity, ensuring demographic homogeneity and relevance to the specific population under study.

scRNA-seq data processing and integration

We integrated scRNA-seq datasets from multiple sources for our analysis. Specifically, scRNA-seq data from 4 RB samples (2 intraocular and 2 extraocular RB patients) from our previous cohort were accessed via the GEO database under accession ID GSE249995.18 Additionally, scRNA-seq data from 7 RB samples (4 intraocular and 3 extraocular RB patients) were obtained from the GEO database with accession number GSE168434.27 For comparative analysis, we included publicly available scRNA-seq datasets from 5 normal retina samples, which were sourced from ArrayExpress under the accession number E-MTAB-7316.45 The procedures for unsupervised clustering analysis, dimensionality reduction, and cell type assignment were meticulously followed in accordance with established protocols.18 The SeuratR package (version v4.0.3) was employed for scRNA-seq data integration.48 Differentially expressed genes were defined with a significance threshold of p < 0.05. Gene Set Enrichment Analysis (GSEA) was performed using the OmicStudio cloud platform (https://www.omicstudio.cn)49 to understand the biological significance of the differentially expressed genes between RB and normal retina samples. The analysis utilized gene sets from the Molecular Signatures Database (MSigDB) and provided enrichment scores and p-values to determine the significance of each gene set. Graphical representations were generated to visualize the results, highlighting key pathways and biological processes altered in retinoblastoma. The Monocle3 software (version 1.2.9) was utilized to infer the pseudo-time trajectory.50 DEGs between different states were determined using the fit_models function, employing a Q value threshold of <0.001. For cell type assignment, we gathered known gene markers from previous studies representing various cell types.51,52,53 Using a semi-supervised machine learning method, we mapped the expression patterns of these markers across all cell clusters. Cells showing distinct expression levels of marker genes were assigned corresponding cell types, with manual adjustments made based on marker gene expression profiles. Visualization of the trajectory was conducted through UMAP graphs, while dynamic expression heatmaps were generated using the Heatmap function within the ComplexHeatmap package.

Cell lines and reagents

Two RB cell lines, WERI-RB1 and Y79, were employed, along with a non-tumor human retinal pigment epithelial cell line known as ARPE-19 (American Type Culture Collection, Manassas, VA, USA). Y79 and WERI-RB1 cells were cultured in RPMI 1640 medium (Corning, USA). ARPE-19 were cultured in Dulbecco’s Modified Eagle Medium/Nutrient Mixture F-12 (DMEM/F-12) (Gibco, Thermo Fisher Scientific, USA). These cultures were all supplemented with 10% fetal bovine serum (Gibco, Thermo Fisher Scientific, USA) and 1% penicillin-streptomycin (Invitrogen, USA). All cell cultures were maintained at 37°C humidified incubator with 5% CO2. Itaconic acid was purchased from Selleck (S3095, Houston, TX, USA), and stock solutions were diluted in DMSO at an initial concentration of 1 mM. Melphalan was purchased from Selleck (S8266, Houston, TX, USA).

Tumor xenograft model

All animal experiments adhered to the ARVO Statement for the Use of Animals in Ophthalmic and Vision Research and were approved and supervised by the Institutional Animal Care and Use Committee at the Zhongshan Ophthalmic Center. Female BALB/C nude mice were procured from Zhuhai BesTest Bio-Tech Co (Zhuhai, China). At four weeks of age, the mice were randomly allocated into three groups: PBS (N = 10 eyes), Itaconate (N = 10 eyes), and melphalan (N = 10 eyes). Intravitreal injections of 1 × 10ˆ5 freshly prepared WERI-RB1 cells were administered using a 33g Hamilton needle, with neomycin ointment applied to the eyes afterward. Mice were monitored every three days for survival and ocular changes. On day 14, mice were examined under a microscope to confirm intraocular tumor presence prior to allocation. Each mouse was then randomly assigned an identification number to ensure unbiased group assignment before receiving intravitreal injections of PBS, 1 mM itaconate, or 1mM melphalan. On day 40, mice were euthanized, and their eyeballs were enucleated, fixed, and paraffin-embedded for subsequent H&E staining. TUNEL staining was performed using the One-Step TUNEL Apoptosis Assay Kit as per the manufacturer’s instructions. The percentage of apoptotic cells per tumor-bearing eye was calculated by averaging the number of TUNEL-positive nuclei divided by the total number of nuclei in the field.

Method details

Western blotting

Cells and RB tissues were collected and then lysed in RIPA buffer (K1020, APE×Bio, USA), supplemented with a protease inhibitor. The protein lysates obtained were assessed using a standard Western blotting protocol to detect the expression levels of ALDOA (1:1000, Abcam, ab252953), SUSD2 (1:1000, Abcam, ab182147), CLK2 (1:500, Proteintech, 11092-1-AP), ARHGAP27 (1:500, Proteintech, 22077-1-AP), β-actin (1:10000, Proteintech, 66009-1-lg) and β-tubulin (1:10000, Proteintech, 66240-1-Ig). Following this, the bound antibodies were detected by incubating the membranes with HRP-conjugated Affinipure Goat Anti-Rabbit IgG(H + L) (1:10000, ProteinTech, SA00001-2), and subsequently, the membranes underwent three rounds of 10-min washing in TBS-0.1% Tween 20. Finally, the bands were visualized using enhanced chemiluminescence (ECL; Tanon, China) and the Tanon 5200 MultiImage System (Tanon, China).

RNA extraction and quantitative real-time qPCR

RNA was extracted using the EZBioscience EZ-press RNA Purification Kit (B0004D, EZBioscience, USA). For cDNA synthesis, 1000 ng of total RNA was utilized along with the EZBioscience 4× Reverse Transcription Master Mix (A001GQ). The resulting reverse transcription product was then combined with other components in a 10 μL reaction for PCR amplification. GAPDH was chosen as the reference gene, and the relative expression levels of the target genes were calculated using the 2-ΔΔCT method. Each experiment was performed in triplicate. Primer sequences for each gene were provided in Table S2.

siRNA transfection of cell line

The small interfering RNA (siRNA) target sequences were synthesized by Shanghai GimaGen Corporation. For transfection, Y79 and WERI-RB1 cells were cultured in 1640 culture medium. On the day of transfection, 100 μL of the same medium was added to each well of a 24-well plate containing 2 × 10ˆ5 cells. 750 ng of siRNA, resulting in a final concentration of 100 nM, was diluted in 100 μL of serum-free 1640 culture medium, and 6 μL of HiPerFect Reagent (301704, Qiagen, Netherlands) was added. After incubation at room temperature for 5–10 min, the transfection complex was added to the cells, followed by co-culture for 6 h. Subsequently, 400 μL of complete 1640 culture medium was added, and cells were cultured for 24–48 h before collection for further experiments. Sequences of siRNA were provided in Table S3.

Cell viability assay

Cell viability of Y79 and WERI-RB1 cells was assessed using the Cell Counting Kit-8 (CCK-8) assay (A311-01, Vazyme, China). In brief, cells were seeded in 96-well plates at a density of 5000 cells per well and subsequently treated with siRNA and itaconic acid. Next, 10 μL of the CCK-8 reagent was introduced into each well, and the plates were incubated for an additional 4 h. The absorbance of the samples was then measured at a wavelength of 450 nm. Cell viability inhibition was calculated using GraphPad Prism software (version 8.0; GraphPad Software, La Jolla, CA, USA), and all experiments were performed in triplicate.

EdU assay

Cells were cultured in 6-well plates and allowed to recover overnight before treatment. A 2X EdU working solution was prepared by diluting EdU (10 mM) 1:500 in cell culture medium to achieve a final concentration of 20 μM. This was added in equal volume to the wells to obtain a final EdU concentration of 10 μM. Cells were incubated for 2 h, then fixed with 1 mL of 4% paraformaldehyde for 15 min at room temperature. Cells were washed three times with 1 mL of wash solution (3–5 min each), permeabilized with 1 mL of PBS containing 0.3% Triton X-100 for 10–15 min, and washed again. The Click Additive Solution was prepared as per the manufacturer’s instructions, and the Click reaction mixture was added (0.5 mL per well) and incubated at room temperature in the dark for 30 min. Cells were washed three more times, and nuclear staining was performed. Fluorescence microscopy was used to observe the cells, with Azide 488 having a maximum excitation wavelength of 495 nm and an emission wavelength of 519 nm. Three biological replicates were conducted.

Cell colony formation assay

Following transfection, Y79 and WERI-RB1 cells were plated onto poly-D-lysine-coated coverslips within 6-well plates, with a cell density of 1×105 cells per well. This culture was maintained for a 24-h period, with regular medium replacements every 3 days. Subsequently, after a 7-day incubation, the cells underwent staining with 0.1% crystal violet dye for a duration of 1 h. Colonies were identified as cell clusters consisting of more than 50 cells and were quantified. Each experiment was repeated three times.

Measurement of extracellular lactate level

The extracellular lactate content was quantified using the Lactate Colorimetric Assay Kit (KTB1100, Abbkine, China). Following cell lysis with lactate extract buffer, the intracellular lactate levels were measured in accordance with the manufacturer’s instructions. Specifically, a total of 5 × 106 cultured cells were washed with PBS and the resulting supernatant (0.8 mL) was collected in a sterile tube post-centrifugation. Standard curves were constructed per the manufacturer’s protocol, and the samples were incubated with the reaction mixture for 30 min. After centrifugation to remove insoluble materials, the supernatants were dissolved in ethanol and the absorbance was measured at a wavelength of 570 nm. Calculation of lactate levels was based on the standard curve.

ADP/ATP ratio assay

The ADP/ATP ratio was determined using the ADP/ATP Ratio Assay Kit (KA1673, Abnova) following the manufacturer’s instructions. Briefly, 1 × 104 treated cells were plated in white opaque 96-well plates. 90 μL of ATP reagent was added to each well, and luminescence (RLU A) was measured after a 1-min incubation using a luminometer. Ten minutes later, luminescence (RLU B) was measured again. Subsequently, 5 μL of ADP reagent was added to each well, and luminescence (RLU C) was quantified after another 1-min incubation. The ADP/ATP ratio was calculated using the formula: ADP/ATP ratio = ((RLU C) - (RLU B))/(RLU A). All experiments were performed in triplicate to ensure accuracy.

LC‒MS analysis

Liquid chromatography‒mass spectrometry (LC‒MS) (Novogene, China) was utilized for mass spectrometry analysis. Cell samples were processed, and supernatants were collected after centrifugation. Following drying and reconstitution, samples were analyzed on a Waters UPLC BEH Amide column with ESI source conditions set at 600°C. Multiple reaction monitoring (MRM) was employed for scanning, and data processing was performed using SCIEX OS software. Metabolite structures were identified through mass spectra analysis and database retrieval, followed by standardization and statistical analysis for data interpretation.

RNA sequencing

Total RNA was recovered from cell cultures using Trizol Reagent (Thermo Fisher Scientific-Invitrogen, Carlsbad, CA), and frozen at −80°C prior to use for RNA sequencing. Sequencing libraries were generated using NEBNext Ultra RNA Library Prep Kit for Illumina (NEB, USA, Catalog #: E7530L) following manufacturer’s recommendations and index codes were added to attribute sequences to each sample. Library quality was assessed on the Agilent 5400 system (Agilent, USA) and quantified by QPCR (1.5 nM). The qualified libraries were pooled and sequenced on the Illumina NovaSeq 6000 using 2 × 150 paired-end reads. Fastp (v0.23.4) was used to trim adapter sequences and low-quality bases in raw FASTQ files.54 The survived reads were mapped onto GRCh38. p13 through STAR (v2.7.11b) and then counted based on GENCODE v39 GTF.55,56 Finally, RSEM (v1.3.3) was utilized to estimate gene expression levels from alignment files.57 The RSEM outputs summarized expected count, transcripts per million (TPM), and fragments per kilobase per million mapped fragments (FPKM). For downstream analysis, genes with zero TPM value were excluded. Identification of differentially expressed genes was performed using R package DESeq2 (v1.39.8).46 The identified DEGs were sorted based on their P-values for further analysis. Annotation, visualization, and integrated discovery were performed using the DAVID.58 This analytical platform was utilized for Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, enabling a comprehensive exploration of the molecular functions of the DEGs. GO analysis was conducted within the biological process (BP) category, with statistical significance defined at p < 0.05. Additionally, Gene Set Enrichment Analysis (GSEA) was performed to further elucidate the functional significance of the DEGs.59 The raw sequence data reported in this paper have been deposited in the Genome Sequence Archive60 in National Genomics Data Center,61 China National Center for Bioinformation/Beijing Institute of Genomics, Chinese Academy of Sciences (GSA-Human: HRA008175) that are publicly accessible at https://ngdc.cncb.ac.cn/gsa-human.

Histological and immunohistochemistry

Four-micrometer-thick paraffin sections obtained from both human and mouse tissues underwent standard Hematoxylin and Eosin (H&E) staining and immunohistochemistry (IHC) procedures following optimized protocols. Tissue specimens for IHC analysis were collected from a cohort comprising 15 extraocular retinoblastoma (RB) patients and 15 intraocular RB patients who had received treatment at the Zhongshan Ophthalmic Center. The assessment of IHC staining in RB samples was conducted by two authors (J.T. and Y.W.), under the supervision of pathologists (P.Z. and S.C.). ALDOA expression levels were evaluated using a 1:100 dilution of anti-ALDOA antibody (Abcam, ab252953), and quantification was performed by determining the proportion of positively stained tumor cells.

Quantification and statistical analysis

The data are presented as means ± standard deviation (SD). Data analysis and presentation were conducted using R software (R Foundation for Statistical Computing, Vienna, Austria, http://www.r-project.org) and GraphPad Prism Software v 8.0.2 (GraphPad, Inc., La Jolla, CA, USA). Statistical analyses were performed utilizing unpaired two-tailed Student’s t test, one-way ANOVA, or Fisher’s exact test as appropriate. Significance levels are indicated in the respective figures.

Supplemental information

Document S1. Figures S1–S4, Tables S1–S3, and Data S1–S11

Table S4. The differentially expressed genes (DEGs) between normal retina and retinoblastoma (RB), related to Figure 1

Supplemental information can be found online at https://doi.org/10.1016/j.isci.2024.110725.
==== Refs
References

1 Kivelä T.T. Hadjistilianou T. Neonatal Retinoblastoma Asia. Pac. J. Oncol. Nurs. 4 2017 197 204 10.4103/apjon.apjon_18_17 28695165
2 Fabian I.D. Onadim Z. Karaa E. Duncan C. Chowdhury T. Scheimberg I. Ohnuma S.I. Reddy M.A. Sagoo M.S. The management of retinoblastoma Oncogene 37 2018 1551 1560 10.1038/s41388-017-0050-x 29321667
3 Zhou C. Wen X. Ding Y. Ding J. Jin M. Liu Z. Wang S. Han M. Yuan H. Xiao Y. Eye-Preserving Therapies for Advanced Retinoblastoma: A Multicenter Cohort of 1678 Patients in China Ophthalmology 129 2022 209 219 10.1016/j.ophtha.2021.09.002 34536465
4 Vander Heiden M.G. Cantley L.C. Thompson C.B. Understanding the Warburg effect: the metabolic requirements of cell proliferation Science 324 2009 1029 1033 10.1126/science.1160809 19460998
5 Levine A.J. Puzio-Kuter A.M. The control of the metabolic switch in cancers by oncogenes and tumor suppressor genes Science 330 2010 1340 1344 10.1126/science.1193494 21127244
6 Chang Y.C. Yang Y.C. Tien C.P. Yang C.J. Hsiao M. Roles of Aldolase Family Genes in Human Cancers and Diseases Trends Endocrinol. Metab. 29 2018 549 559 10.1016/j.tem.2018.05.003 29907340
7 Jiang Z. Wang X. Li J. Yang H. Lin X. Aldolase A as a prognostic factor and mediator of progression via inducing epithelial-mesenchymal transition in gastric cancer J. Cell Mol. Med. 22 2018 4377 4386 10.1111/jcmm.13732 29992789
8 Li X. Jiang F. Ge Z. Chen B. Yu J. Xin M. Wang J. An L. Wei J. Wu L. Fructose-Bisphosphate Aldolase A Regulates Hypoxic Adaptation in Hepatocellular Carcinoma and Involved with Tumor Malignancy Dig. Dis. Sci. 64 2019 3215 3227 10.1007/s10620-019-05642-2 31041640
9 Ye F. Chen Y. Xia L. Lian J. Yang S. Aldolase A overexpression is associated with poor prognosis and promotes tumor progression by the epithelial-mesenchymal transition in colon cancer Biochem. Biophys. Res. Commun. 497 2018 639 645 10.1016/j.bbrc.2018.02.123 29453983
10 Saito Y. Takasawa A. Takasawa K. Aoyama T. Akimoto T. Ota M. Magara K. Murata M. Hirohashi Y. Hasegawa T. Aldolase A promotes epithelial-mesenchymal transition to increase malignant potentials of cervical adenocarcinoma Cancer Sci. 111 2020 3071 3081 10.1111/cas.14524 32530543
11 Tang J. Liu Y. Wang Y. Zhang Z. Nie J. Wang X. Ai S. Li J. Gao Y. Li C. Deciphering metabolic heterogeneity in retinoblastoma unravels the role of monocarboxylate transporter 1 in tumor progression Biomark. Res. 12 2024 48 10.1186/s40364-024-00596-8 38730450
12 Rajasekaran S. Nagarajha Selvan L.D. Dotts K. Kumar R. Rishi P. Khetan V. Bisht M. Sivaraman K. Krishnakumar S. Sahoo D. Non-coding and Coding Transcriptional Profiles Are Significantly Altered in Pediatric Retinoblastoma Tumors Front. Oncol. 9 2019 221 10.3389/fonc.2019.00221 31058073
13 Qin W. Qin K. Zhang Y. Jia W. Chen Y. Cheng B. Peng L. Chen N. Liu Y. Zhou W. S-glycosylation-based cysteine profiling reveals regulation of glycolysis by itaconate Nat. Chem. Biol. 15 2019 983 991 10.1038/s41589-019-0323-5 31332308
14 Yao Y. Gu X. Xu X. Ge S. Jia R. Novel insights into RB1 mutation Cancer Lett. 547 2022 215870 10.1016/j.canlet.2022.215870
15 Chang Y.C. Chan Y.C. Chang W.M. Lin Y.F. Yang C.J. Su C.Y. Huang M.S. Wu A.T.H. Hsiao M. Feedback regulation of ALDOA activates the HIF-1α/MMP9 axis to promote lung cancer progression Cancer Lett. 403 2017 28 36 10.1016/j.canlet.2017.06.001 28610954
16 Li X. Yu C. Luo Y. Lin J. Wang F. Sun X. Gao Y. Tan W. Xia Q. Kong X. Aldolase A Enhances Intrahepatic Cholangiocarcinoma Proliferation and Invasion through Promoting Glycolysis Int. J. Biol. Sci. 17 2021 1782 1794 10.7150/ijbs.59068 33994862
17 Yang J. Li Y. Han Y. Feng Y. Zhou M. Zong C. He X. Jia R. Xu X. Fan J. Single-cell transcriptome profiling reveals intratumoural heterogeneity and malignant progression in retinoblastoma Cell Death Dis. 12 2021 1100 10.1038/s41419-021-04390-4 34815392
18 Liu Y. Hu W. Xie Y. Tang J. Ma H. Li J. Nie J. Wang Y. Gao Y. Cheng C. Single-cell transcriptomics enable the characterization of local extension in retinoblastoma Commun. Biol. 7 2024 11 10.1038/s42003-023-05732-y 38172218
19 Tang J. Liu Y. Zhang Z. Ren Y. Ma Y. Wang Y. Li J. Gao Y. Li C. Cheng C. Heterogeneous Expression Patterns of the Minichromosome Maintenance Complex Members in Retinoblastoma Unveil Its Clinical Significance Invest. Ophthalmol. Vis. Sci. 65 2024 31 10.1167/iovs.65.1.31
20 Berry J.L. Xu L. Kooi I. Murphree A.L. Prabakar R.K. Reid M. Stachelek K. Le B.H.A. Welter L. Reiser B.J. Genomic cfDNA Analysis of Aqueous Humor in Retinoblastoma Predicts Eye Salvage: The Surrogate Tumor Biopsy for Retinoblastoma Mol. Cancer Res. 16 2018 1701 1712 10.1158/1541-7786.Mcr-18-0369 30061186
21 Ghose N. Kaliki S. Liquid biopsy in Retinoblastoma: A review Semin. Ophthalmol. 37 2022 813 819 10.1080/08820538.2022.2078165 35604935
22 Schmidt M.J. Prabakar R.K. Pike S. Yellapantula V. Peng C.C. Kuhn P. Hicks J. Xu L. Berry J.L. Simultaneous Copy Number Alteration and Single-Nucleotide Variation Analysis in Matched Aqueous Humor and Tumor Samples in Children with Retinoblastoma Int. J. Mol. Sci. 24 2023 8606 10.3390/ijms24108606
23 Ju S. Chen C. Zhang J. Xu L. Zhang X. Li Z. Chen Y. Zhou J. Ji F. Wang L. Detection of circulating tumor cells: opportunities and challenges Biomark. Res. 10 2022 58 10.1186/s40364-022-00403-2 35962400
24 Li H.T. Xu L. Weisenberger D.J. Li M. Zhou W. Peng C.C. Stachelek K. Cobrinik D. Liang G. Berry J.L. Characterizing DNA methylation signatures of retinoblastoma using aqueous humor liquid biopsy Nat. Commun. 13 2022 5523 10.1038/s41467-022-33248-2 36130950
25 Luo Y. Xu M. Yang L. Yao Y. Berry J.L. Xu L. Wen X. He X. Han M. Fan X. Correlating somatic copy number alteration in aqueous humour cfDNA with chemotherapy history, eye salvage and pathological features in retinoblastoma Br. J. Ophthalmol. 108 2024 449 456 10.1136/bjo-2022-322866 36931696
26 Liu H. Zhang Y. Zhang Y.Y. Li Y.P. Hua Z.Q. Zhang C.J. Wu K.C. Yu F. Zhang Y. Su J. Jin Z.B. Human embryonic stem cell-derived organoid retinoblastoma reveals a cancerous origin Proc. Natl. Acad. Sci. USA 117 2020 33628 33638 10.1073/pnas.2011780117 33318192
27 Wu C. Yang J. Xiao W. Jiang Z. Chen S. Guo D. Zhang P. Liu C. Yang H. Xie Z. Single-cell characterization of malignant phenotypes and microenvironment alteration in retinoblastoma Cell Death Dis. 13 2022 438 10.1038/s41419-022-04904-8 35523772
28 Dimaras H. Khetan V. Halliday W. Orlic M. Prigoda N.L. Piovesan B. Marrano P. Corson T.W. Eagle R.C. Jr. Squire J.A. Gallie B.L. Loss of RB1 induces non-proliferative retinoma: increasing genomic instability correlates with progression to retinoblastoma Hum. Mol. Genet. 17 2008 1363 1372 10.1093/hmg/ddn024 18211953
29 Qian Y. Yin Y. Zheng X. Liu Z. Wang X. Metabolic regulation of tumor-associated macrophage heterogeneity: insights into the tumor microenvironment and immunotherapeutic opportunities Biomark. Res. 12 2024 1 10.1186/s40364-023-00549-7 38185636
30 Vitale I. Manic G. Coussens L.M. Kroemer G. Galluzzi L. Macrophages and Metabolism in the Tumor Microenvironment Cell Metab. 30 2019 36 50 10.1016/j.cmet.2019.06.001 31269428
31 Tian W. Zhou J. Chen M. Qiu L. Li Y. Zhang W. Guo R. Lei N. Chang L. Bioinformatics analysis of the role of aldolase A in tumor prognosis and immunity Sci. Rep. 12 2022 11632 10.1038/s41598-022-15866-4
32 Lu Y. Zhang Y. Wang X. Zhang H. Zhu Y. Zhang J. Sha H. Zou R. Gan Y. Sui Y. Aldolase A Promotes Colorectal Cancer Progression through Targeting COPS6 and Regulating MAPK Signaling Pathway Dis. Markers 2023 2023 1702125 10.1155/2023/1702125
33 Zhao Y. Li Z.X. Zhu Y.J. Fu J. Zhao X.F. Zhang Y.N. Wang S. Wu J.M. Wang K.T. Wu R. Single-Cell Transcriptome Analysis Uncovers Intratumoral Heterogeneity and Underlying Mechanisms for Drug Resistance in Hepatobiliary Tumor Organoids Adv. Sci. 8 2021 e2003897 10.1002/advs.202003897
34 Chen L. Wu Z. Guo J. Wang X. Zhao Z. Liang H. Zhang R. Deng J. Initial clinical and experimental analyses of ALDOA in gastric cancer, as a novel prognostic biomarker and potential therapeutic target Clin. Exp. Med. 23 2023 2443 2456 10.1007/s10238-022-00952-8 36422738
35 Wang C. Xu J. Yuan D. Bai Y. Pan Y. Zhang J. Shao C. Exosomes carrying ALDOA and ALDH3A1 from irradiated lung cancer cells enhance migration and invasion of recipients by accelerating glycolysis Mol. Cell. Biochem. 469 2020 77 87 10.1007/s11010-020-03729-3 32297178
36 Lin J. Xia L. Oyang L. Liang J. Tan S. Wu N. Yi P. Pan Q. Rao S. Han Y. The POU2F1-ALDOA axis promotes the proliferation and chemoresistance of colon cancer cells by enhancing glycolysis and the pentose phosphate pathway activity Oncogene 41 2022 1024 1039 10.1038/s41388-021-02148-y 34997215
37 Lampropoulou V. Sergushichev A. Bambouskova M. Nair S. Vincent E.E. Loginicheva E. Cervantes-Barragan L. Ma X. Huang S.C.C. Griss T. Itaconate Links Inhibition of Succinate Dehydrogenase with Macrophage Metabolic Remodeling and Regulation of Inflammation Cell Metab. 24 2016 158 166 10.1016/j.cmet.2016.06.004 27374498
38 Mills E.L. Ryan D.G. Prag H.A. Dikovskaya D. Menon D. Zaslona Z. Jedrychowski M.P. Costa A.S.H. Higgins M. Hams E. Itaconate is an anti-inflammatory metabolite that activates Nrf2 via alkylation of KEAP1 Nature 556 2018 113 117 10.1038/nature25986 29590092
39 Weiss J.M. Davies L.C. Karwan M. Ileva L. Ozaki M.K. Cheng R.Y. Ridnour L.A. Annunziata C.M. Wink D.A. McVicar D.W. Itaconic acid mediates crosstalk between macrophage metabolism and peritoneal tumors J. Clin. Invest. 128 2018 3794 3805 10.1172/jci99169 29920191
40 Watson A.P. Evans R.L. Egland K.A. Multiple functions of sushi domain containing 2 (SUSD2) in breast tumorigenesis Mol. Cancer Res. 11 2013 74 85 10.1158/1541-7786.Mcr-12-0501-t 23131994
41 Youn M. Smith S.M. Lee A.G. Chae H.D. Spiteri E. Erdmann J. Galperin I. Jones L.M. Donato M. Abidi P. Comparison of the Transcriptomic Signatures in Pediatric and Adult CML Cancers 13 2021 6263 10.3390/cancers13246263
42 Yang Y. Wu L. Shu X. Lu Y. Shu X.O. Cai Q. Beeghly-Fadiel A. Li B. Ye F. Berchuck A. Genetic Data from Nearly 63,000 Women of European Descent Predicts DNA Methylation Biomarkers and Epithelial Ovarian Cancer Risk Cancer Res. 79 2019 505 517 10.1158/0008-5472.Can-18-2726 30559148
43 Yoshida T. Kim J.H. Carver K. Su Y. Weremowicz S. Mulvey L. Yamamoto S. Brennan C. Mei S. Long H. CLK2 Is an Oncogenic Kinase and Splicing Regulator in Breast Cancer Cancer Res. 75 2015 1516 1526 10.1158/0008-5472.Can-14-2443 25670169
44 Park S.Y. Piao Y. Thomas C. Fuller G.N. de Groot J.F. Cdc2-like kinase 2 is a key regulator of the cell cycle via FOXO3a/p27 in glioblastoma Oncotarget 7 2016 26793 26805 10.18632/oncotarget.8471 27050366
45 Lukowski S.W. Lo C.Y. Sharov A.A. Nguyen Q. Fang L. Hung S.S. Zhu L. Zhang T. Grünert U. Nguyen T. A single-cell transcriptome atlas of the adult human retina Embo J. 38 2019 e100811 10.15252/embj.2018100811
46 Love M.I. Huber W. Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2 Genome Biol. 15 2014 550 10.1186/s13059-014-0550-8 25516281
47 Fabian I.D. Reddy A. Sagoo M.S. Classification and staging of retinoblastoma Community Eye Health 31 2018 11 13 29915461
48 Hao Y. Hao S. Andersen-Nissen E. Mauck W.M. 3rd Zheng S. Butler A. Lee M.J. Wilk A.J. Darby C. Zager M. Integrated analysis of multimodal single-cell data Cell 184 2021 3573 3587.e29 10.1016/j.cell.2021.04.048 34062119
49 Lyu F. Han F. Ge C. Mao W. Chen L. Hu H. Chen G. Lang Q. Fang C. OmicStudio: A composable bioinformatics cloud platform with real-time feedback that can generate high-quality graphs for publication iMeta 2 2023 e85 38868333
50 Trapnell C. Cacchiarelli D. Grimsby J. Pokharel P. Li S. Morse M. Lennon N.J. Livak K.J. Mikkelsen T.S. Rinn J.L. The dynamics and regulators of cell fate decisions are revealed by pseudotemporal ordering of single cells Nat. Biotechnol. 32 2014 381 386 10.1038/nbt.2859 24658644
51 Cowan C.S. Renner M. De Gennaro M. Gross-Scherf B. Goldblum D. Hou Y. Munz M. Rodrigues T.M. Krol J. Szikra T. Cell Types of the Human Retina and Its Organoids at Single-Cell Resolution Cell 182 2020 1623 1640.e34 10.1016/j.cell.2020.08.013 32946783
52 Menon M. Mohammadi S. Davila-Velderrain J. Goods B.A. Cadwell T.D. Xing Y. Stemmer-Rachamimov A. Shalek A.K. Love J.C. Kellis M. Hafler B.P. Single-cell transcriptomic atlas of the human retina identifies cell types associated with age-related macular degeneration Nat. Commun. 10 2019 4902 10.1038/s41467-019-12780-8 31653841
53 Collin J. Queen R. Zerti D. Steel D.H. Bowen C. Parulekar M. Lako M. Dissecting the Transcriptional and Chromatin Accessibility Heterogeneity of Proliferating Cone Precursors in Human Retinoblastoma Tumors by Single Cell Sequencing-Opening Pathways to New Therapeutic Strategies? Invest. Ophthalmol. Vis. Sci. 62 2021 18 10.1167/iovs.62.6.18
54 Chen S. Ultrafast one-pass FASTQ data preprocessing, quality control, and deduplication using fastp Imeta 2 2023 e107 10.1002/imt2.107
55 Dobin A. Davis C.A. Schlesinger F. Drenkow J. Zaleski C. Jha S. Batut P. Chaisson M. Gingeras T.R. STAR: ultrafast universal RNA-seq aligner Bioinformatics 29 2013 15 21 10.1093/bioinformatics/bts635 23104886
56 Harrow J. Frankish A. Gonzalez J.M. Tapanari E. Diekhans M. Kokocinski F. Aken B.L. Barrell D. Zadissa A. Searle S. GENCODE: the reference human genome annotation for The ENCODE Project Genome Res. 22 2012 1760 1774 10.1101/gr.135350.111 22955987
57 Li B. Dewey C.N. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome BMC Bioinf. 12 2011 323 10.1186/1471-2105-12-323
58 Dennis G. Jr. Sherman B.T. Hosack D.A. Yang J. Gao W. Lane H.C. Lempicki R.A. DAVID: Database for Annotation, Visualization, and Integrated Discovery Genome Biol. 4 2003 P3 12734009
59 Subramanian A. Tamayo P. Mootha V.K. Mukherjee S. Ebert B.L. Gillette M.A. Paulovich A. Pomeroy S.L. Golub T.R. Lander E.S. Mesirov J.P. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles Proc. Natl. Acad. Sci. USA 102 2005 15545 15550 10.1073/pnas.0506580102 16199517
60 Chen T. Chen X. Zhang S. Zhu J. Tang B. Wang A. Dong L. Zhang Z. Yu C. Sun Y. The Genome Sequence Archive Family: Toward Explosive Data Growth and Diverse Data Types Dev. Reprod. Biol. 19 2021 578 583 10.1016/j.gpb.2021.08.001
61 CNCB-NGDC Members and Partners Database Resources of the National Genomics Data Center, China National Center for Bioinformation in 2022 Nucleic Acids Res. 50 2021 D27 D38 10.1093/nar/gkab951
