Next Article in Journal
A Galactooligosaccharide Product Decreases the Rotavirus Infection in Suckling Rats
Previous Article in Journal
Features of ZED1227: The First-In-Class Tissue Transglutaminase Inhibitor Undergoing Clinical Evaluation for the Treatment of Celiac Disease
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cells
Volume 11
Issue 10
10.3390/cells11101668
cells-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Integrated Analysis of Cancer Tissue and Vitreous Humor from Retinoblastoma Eyes Reveals Unique Tumor-Specific Metabolic and Cellular Pathways in Advanced and Non-Advanced Tumors

by
Vishnu Suresh Babu
1,2,
Ashwin Mallipatna
3,
Deepak SA
4,
Gagan Dudeja
3,
Ramaraj Kannan
1,
Rohit Shetty
1,
Archana Padmanabhan Nair
1,
Seetharamanjaneyulu Gundimeda
4,
Shyam S. Chaurasia
5,
Navin Kumar Verma
6,7,
Rajamani Lakshminarayanan
7,
Stephane Heymans
2,8,*,
Veluchamy A. Barathi
7,9,10,
Nilanjan Guha
4,* and
Arkasubhra Ghosh
1,*
1
GROW Research Laboratory, Narayana Nethralaya Foundation, Bangalore 560099, India
2
Department of Cardiology, Cardiovascular Research Institute Maastricht (CARIM), Maastricht University, 6229 ER Maastricht, The Netherlands
3
Retinoblastoma Service, Narayana Nethralaya, Bangalore 560099, India
4
Agilent Technologies India Pvt Ltd., Bangalore 560048, India
5
Ophthalmology and Visual Sciences, Medical College of Wisconsin, Milwaukee, WI 53226, USA
6
Lee Kong Chian School of Medicine, Nanyang Technological University Singapore, Singapore 308232, Singapore
7
Singapore Eye Research Institute, Singapore 169856, Singapore
8
Centre for Molecular and Vascular Biology, Department of Cardiovascular Sciences, KU Leuven, Herestraat 49, Bus 911, 3000 Leuven, Belgium
9
The Ophthalmology and Visual Sciences ACP, Duke-NUS Medical School, Singapore 169857, Singapore
10
Department of Ophthalmology, Yong Loo Lin School of Medicine, National University of Singapore, Singapore 117597, Singapore
*
Authors to whom correspondence should be addressed.
Cells 2022, 11(10), 1668; https://doi.org/10.3390/cells11101668
Submission received: 24 April 2022 / Revised: 10 May 2022 / Accepted: 11 May 2022 / Published: 18 May 2022
(This article belongs to the Topic Cell Signaling Pathways)
Browse Figures
Review Reports Versions Notes

Abstract

:
Retinoblastoma (Rb) is a pediatric intraocular malignancy that is proposed to originate from maturing cone cell precursors in the developing retina. The molecular mechanisms underlying the biological and clinical behaviors are important to understand in order to improve the management of advanced-stage tumors. While the genetic causes of Rb are known, an integrated understanding of the gene expression and metabolic processes in tumors of human eyes is deficient. By integrating transcriptomic profiling from tumor tissues and metabolomics from tumorous eye vitreous humor samples (with healthy, age-matched pediatric retinae and vitreous samples as controls), we uncover unique functional associations between genes and metabolites. We found distinct gene expression patterns between clinically advanced and non-advanced Rb. Global metabolomic analysis of the vitreous humor of the same Rb eyes revealed distinctly altered metabolites, indicating how tumor metabolism has diverged from healthy pediatric retina. Several key enzymes that are related to cellular energy production, such as hexokinase 1, were found to be reduced in a manner corresponding to altered metabolites; notably, a reduction in pyruvate levels. Similarly, E2F2 was the most significantly elevated E2F family member in our cohort that is part of the cell cycle regulatory circuit. Ectopic expression of the wild-type RB1 gene in the Rb-null Y79 and WERI-Rb1 cells rescued hexokinase 1 expression, while E2F2 levels were repressed. In an additional set of Rb tumor samples and pediatric healthy controls, we further validated differences in the expression of HK1 and E2F2. Through an integrated omics analysis of the transcriptomics and metabolomics of Rb, we uncovered a significantly altered tumor-specific metabolic circuit that reduces its dependence on glycolytic pathways and is governed by Rb1 and HK1.

1. Introduction

Retinoblastoma (Rb) is the most common pediatric intraocular tumor, characterized by the presence of retinal lesions with vascularization, vitreous seeding, retinal detachment, and calcification [1]. The major clinical goal in treating children with Rb involves preserving life and vision. It is caused by the loss of the functional retinoblastoma protein (pRb), affecting approximately 8000 children annually [2]. Since pRb is a tumor suppressor protein, bi-allelic mutations in RB1 and the amplification of oncogenes such as MYCN relieve cell-cycle checkpoints, thereby driving neoplastic transformation in precursor retinal cells [3]. RB1 loss causes various alterations to cellular functions through transcriptional and epigenetic mechanisms [4], a greater understanding of which can provide novel treatment strategies. Although significant advances have been made to understand the genetics [5] and pathogenesis of retinoblastoma [6], it is still the most aggressive and fatal pediatric cancer in developing countries [7,8]. Rb is curable with global salvation strategies [9], and is universally fatal if left untreated. Most Rb tumors are diagnosed at advanced stages, which exhibit resistance to conventional therapy and possess a metastatic risk [10]. Therefore, understanding the molecular mechanism of the Rb tumor subtypes, particularly the molecular milieu in advanced tumors to develop targeted therapies, remains the most important prerequisite for better treatment outcomes [11].
Although there have been numerous reports on specific molecular aspects of retinoblastoma [12], concerted efforts to identify molecular pathways and factors that are conserved across the genetic, transcriptomic, and metabolomic profiles in the same tumor tissues are lacking. Such a unified investigation strategy is well suited to reveal better correlations between genetic, gene expression, and metabolic profiles that can be mechanistically linked to disease progression [13]. In this study, we analyzed the transcriptomic and metabolomic profiles of Rb tumors and vitreous humor to identify markers that are associated with the different stages of the cancer. Our goal was to identify, in human tumor samples, specific markers and enriched pathways that are specific to the cancer stage, which can be used to better understand the disease and to tailor treatment modalities more accurately, an unmet clinical need.

2. Materials and Methods

2.1. Clinical Samples

The study was conducted by the Declaration of Helsinki principles under a protocol approved by the institutional ethics committee of Narayana Nethralaya (EC Ref. no: C/2013/03/02). Informed written consent was received from all parents before inclusion in the study. Histology confirmed Rb tumors (n = 9) comprising Group E and Group D of the age range 0.2–4 years, and pediatric controls (n = 2) of the age range 0.2–0.3 years, were used for the microarray and metabolomics study. The details of clinical samples, including age, gender, laterality, tumor viability, and clinical and histopathology details, are mentioned in Table 1. For the immunohistochemistry validations, we have used additional Rb subjects (n = 25) comprising Group E and Group D of the age range 0.2–4 years. Clinical and histopathology details are mentioned in Supplementary Table S1.

2.2. Tumor mRNA Profiling

Total RNA was isolated from 9 Rb tumors and 2 control pediatric retinal samples, using the Agilent Absolutely RNA miRNA kit (cat#400814, Agilent Technologies, Santa Clara, CA, USA) according to the manufacturer’s instructions. For the cell culture microarray, total RNA was isolated from RB1-null Y79 and RB1-overexpressed Y79 cells, using the Agilent Absolutely RNA miRNA kit. Twenty-five nanograms of RNA was taken from each Rb tumor and each control pediatric retina, and the cells were labeled with Cy3 dye using an Agilent Low Input Quick Amp Labeling Kit (p/n 5190-2305, Agilent Technologies, Santa Clara, CA, USA). Gene expression microarray analysis was performed using the Agilent SurePrint G3 Human GE 8 × 60 K V2 Microarray (Agilent Technologies, Santa Clara, CA, USA) and an Agilent SureScan Microarray scanner. The gene expression data were extracted using Agilent Feature Extraction Software (11.5.1.1) and analyzed using Agilent GeneSpring GX 13.1. The analysis was carried out using a t-test unpaired statistical method with the Benjamini–Hochberg FDR method.

2.3. Metabolomics

The samples from 9 patients and 2 controls were extracted using methanol:ethanol (1:1 v/v). The extracted samples were subjected to LC/QTOF-MS analysis, using an Agilent 1290 Infinity LC System coupled to an Agilent 6550 Accurate mass QTOF LC-MS system (Agilent Technologies, Santa Clara, CA, USA). The data was acquired using electrospray ionization in positive and negative ion modes using a modified polar reverse-phase C18 column and an HILIC column. Molecular features were detected using Agilent MassHunter Profinder (v. B.06.00), and were searched and confirmed by matching them against the Agilent METLIN MS/MS library. Agilent Mass Profiler Professional software (MPP) was used for the statistical comparison of the LC/MS data from the Rb and control samples.

2.4. Multi-Omic Data Analysis

The metabolomics and gene microarray results were combined and analyzed using a pathway-centric approach. Transcriptomics and metabolomics data were co-visualized in a pathway context using the Multi-Omics Analysis tool of GeneSpring GX 13.1. This enabled the identification of the differential pathways and entities across both gene expression and metabolomics. Information related to Agilent Technologies products described in this manuscript is for research use only, and not for use in diagnostic procedures.

2.5. Pathway Enrichment Analysis and Regulatory Network Analysis

The KEGG functional enrichment of the tumor microarray, metabolomics, and the cell culture microarray was carried out using the GeneSpring GX 13.1, NetworkAnalyst 3.0, and MetaboAnalyst 5.0 packages; and the pathways with p < 0.05 and fold-change > 2.0 were significantly enriched. KEGG pathway enrichment analysis using GeneSpring GX 13.1 was also conducted on the up- and downregulated genes and metabolites that were involved in relationship pairs. Based on the interaction information, the construction of gene interactions was performed using NetworkAnalyst and Cytoscape 3.2. The entire metabolome mapping and gene-metabolite interactions network was constructed using MetaboAnalyst 5.0.

2.6. Gene Expression Analysis

Total RNA that was extracted for microarray from the clinical subjects was also used for RT-PCR validation. RT-PCR was performed with Agilent Brilliant III Ultra-Fast RT-PCR reagent (cat#600884, Agilent Technologies, Santa Clara, CA, USA), using Agilent AriaMX real-time PCR instruments. The relative mRNA expression levels were quantified using the C(t) method. For in vitro assays, total RNA was isolated from cells using the Trizol reagent (Invitrogen, Waltham, MA, USA) according to the manufacturer’s protocol. A total of 1 µg of RNA was reverse transcribed using the Bio-Rad iScript cDNA synthesis kit (cat#1708890, Bio-Rad, Hercules, CA, USA), and quantitative real-time PCR was performed using the Kappa Sybr Fast qPCR kit (cat#KK4601, Kapa Biosystems Pty (Ltd.), Cape Town, South Africa) using the Bio-Rad CFX96 system. Relative mRNA expression levels were quantified using the C(t) method. The results were normalized to housekeeping human β-actin. The details of the primers used are described in Table S3.

2.7. Cell Lines

WERI-Rb1 and Y79 cells were obtained from the American Type Culture Collection (ATCC, Manassas, VA, USA). Cells were cultured in RPMI 1640 medium (Gibco, Grand Island, NY, USA, cat#11875093) supplemented with 10% FBS and 1% Pen Strep (Penicillin–Streptomycin), and maintained at 37 °C in a humidified atmosphere of 5% CO2, with intermittent shaking in an upright T25 flask.

2.8. Histopathology & Light Microscopy

Paraffin-embedded specimens of Rb tumor (n = 9) and control retina (n = 2) were used. A total of 4 µm paraffin sections were dewaxed at 60 °C and rehydrated in decreasing concentrations of ethanol. Slides were stained with hematoxylin and eosin according to standard procedures. Brightfield images were captured using an Olympus CKX53 microscope and Carl Zeiss Axioplan 2 microscope (Carl Zeiss, Oberkochen, Germany).

2.9. Immunohistochemistry

For IHC, 4 µm sections of Rb tumor (n = 25) and pediatric retina (n = 2) were deparaffinized and rehydrated, and were subjected to heat-induced epitope retrieval using citrate buffer for 20 min at 100 °C. After an endogenous peroxidase block, tissues were incubated overnight at 4 °C with antibodies for Ki67 (1:1000; cat#ab16667, Abcam, Cambridge, UK), Rb (1:500; cat#9309, Cell Signaling Technology, Danvers, MA, USA), phospho-Rb (1:500; cat#8516, Cell Signaling Technology, Danvers, MA, USA), E2F1(1:1000; cat#sc251, Santa Cruz Biotechnologies, Dallas, TX, USA), E2F2 (1:500, cat#209662, Abcam, Cambridge, UK), and HK1 (1:500, cat#ab55144, Abcam, Cambridge, UK). The signals were developed using the DAB substrate (Dako Envision System, Agilent Technologies, Santa Clara, CA, USA). Images were analyzed and captured at Brightfield using an Olympus CKX53 microscope. Following the assessment of the staining using manual scoring, the scores were tabulated as ranks over a range of 0–3, with 0 indicating less than 10% positive cells and 3 indicating more than 80% positive cells.

2.10. Lentiviral Plasmid and Vector

We constructed a lentiviral plasmid expressing the RB1 gene in the pCL20 backbone. The cloning strategy was designed using ApE software (Version 8.5.2.0). The pCL20c vector was digested with the 5′ SmaI site and the 3′ KpnI site to incorporate the RB1 gene. Forward and reverse primers with compatible restriction sites were designed to separate the RB1 gene via PCR amplification from pSG5L HA-RB (Addgene, Cambridge, MA, USA). The details of the cloning primers used are:
RB1 Forward: 5′-CCTACGACGTGCCCGACTACG-3′.
RB1 Reverse: 3′-AGTGACCGGTTCATTTCTCTTCCTTGTTTGAG-5′.
Post-ligation, the product was transformed using DH5α and the colonies were screened using restriction mapping and confirmed using sequencing.

2.11. Western Blotting

For Western blot analysis, cells were lysed in RIPA buffer (20 mM Tris pH 8.0, 0.1% SDS, 150 mM NaCl, 0.08% sodium deoxycholate, and 1% NP40 supplemented with 1 tablet of protease inhibitor (Complete ultra mini-tablet, Roche, Indianapolis, IN, USA) and phosphatase inhibitor (PhosStop tablet, Roche, Indianapolis, IN, USA)). A total of 20 µg of total protein was loaded per lane, and protein was separated using SDS-PAGE. The separated proteins on the gel were transferred onto a PVDF membrane and were probed for specific antibodies against Rb (cat#9390; Cell signaling, Danvers, MA, USA), phospho-Rb (cat#8516, Cell signaling, Danvers, MA, USA), E2F2 (ab209662; Abcam, Cambridge, UK), HK1(cat#2024; Cell signaling, Danvers, MA, USA), and GAPDH (cat#5174; Cell signaling, Danvers, MA, USA) at 1:1000 dilution in 5% BSA in 1× TBST, overnight at 4 °C. After 4 washes with 1× TBST for 10 min, membranes were incubated with HRP-conjugated anti-mouse (cat#7076; Cell signaling, Danvers, MA, USA) or anti-rabbit antibodies (cat#7074; Cell signaling, Danvers, MA, USA) at 1:2000 dilution for 2 h. Images were visualized using the Image Quant LAS 500 system (GE Healthcare Life Sciences, Piscataway, NJ, USA).

2.12. Statistical Analysis

A statistical analysis was performed using GraphPad Prism 8. Data are presented as the mean ± S.D. unless indicated otherwise, and p < 0.05 was considered as being statistically significant. For all representative images, results were reproduced at least three times in independent experiments. For all quantitative data, the statistical test used is indicated in the legends. A statistical ‘decision tree’ is shown in Supplementary Figure S4. Heatmaps of the Z-transformed gene expression level of the mRNA microarray were created using Python 3.7 Seaborne 0.9.0.

3. Results

3.1. Transcriptomic Profiling of Retinoblastoma Tumors Identifies Distinct Expression Profiles in Rb Subtypes

To unravel the molecular networks in Rb tumors, we performed total mRNA profiling using gene expression microarrays in enucleated tumor tissues from nine retinoblastoma patients and two age-matched pediatric retinae as the control (Table 1), as the discovery cohort. We identified distinct dysregulated gene clusters, implicating gross differences between Rb tumors and the controls (Figure 1A). The top upregulated genes identified in the microarray comprised RB1 pathway-related transcription factors [14] such as E2F1 (p < 0.05, FC = 29.2), E2F2 (p < 0.05, FC = 540.2), and key cell cycle checkpoint genes [15] such as CCNB2 (p < 0.05, FC = 319.4), CCNE2 (p < 0.05, FC = 46.81), CDK1 (p < 0.05, FC = 39.81), CDKN2A (p < 0.05, FC = 29.28), and CHEK2 (p < 0.05, FC = 26.9). We also identified the significant upregulation of immune system-related genes such as CD86 (p < 0.05, FC = 493.5) and CD19 (p < 0.05, FC = 13.27) [16], and epigenetic regulators such as SYK [17] (p < 0.05, FC = 14.9) and PRDM1 [18] (p < 0.05, FC = 15.89) in Rb tumors. Strikingly, we also identified the mitochondrial TCA-related FAHD1 gene [19] (p < 0.05, FC = 2.38) to be significantly upregulated in Rb subjects. The top downregulated genes comprised photoreceptor-related genes such as RHO [20] (p < 0.05, FC = −1058.9), NRL [21] (p < 0.05, FC = −104.00), PDE6D [22] (p < 0.05, FC = −7.178), CRABP1 [23] (p < 0.05, FC = −140.60), and glycolytic factors such as HK1 [24] (p < 0.05, FC = −17.46), SLC2A1 [25] (p < 0.05, FC = −8.16), and FOXO3 [26] (p < 0.05, FC = −5.08). The methyltransferase MGMT [27] (p < 0.05, FC = −4.8) and RB1 (p < 0.05, FC = −4.36) were significantly downregulated in Rb tumors. We applied a standard QC procedure to the dataset considering the sample quality, hybridization quality, signal comparability, and array correlation. The boxplot of raw intensities (Figure S1A) and the density histogram of the log intensity distribution (Figure 1B) of each array before normalization provide an overview of the dataset quality. Normalization was able to sufficiently remove smaller discrepancies, leading to comparable distributions between all arrays (Figure 1C and Figure S1B). The findings prompted us to elucidate the microarray analysis in the Rb subtypes. We segregated the Rb cohort into advanced Rb (defined as AJCC Stage [28] cT3 or IIRC [29]) and non-advanced Rb (defined as AJCC Stage cT2 or IIRC Group D), based on their clinical and histopathological information (advanced Rb, n = 5; non-advanced Rb, n = 4). Using the clustering analysis functionality of GeneSpring GX 13.1 (Agilent Technologies, Santa Clara, CA, USA), we found that the advanced Rb subjects (P1–P5) strongly correlated with each other, compared to the non-advanced (P6–P9) and the controls (Figure S1C). However, advanced Rb subjects exhibited a weak correlation with the controls in the microarray (Figure 1D, rs = 0.42), which is suggestive of their greater de-differentiation state, while non-advanced Rb subjects showed a better correlation with the controls (Figure 1E, rs = 0.9). We identified distinct clusters of gene expression profiles between the advanced Rb, non-advanced Rb, and the controls (Figure 1F). In the advanced Rb tumors, we identified 6089 genes using microarray, of which 1027 gene sets were unique to the clinical pathology-defined high-risk tumor (p < 0.05, FC > 2). We identified 2633 genes that were unique to the non-advanced Rb subtype out of 7695 genes detected in the microarray (p < 0.05, FC > 2), compared to age-matched healthy retinae (Figure 1G).

3.2. Rb Tumor Clinical Subtypes Demonstrate Altered Molecular Pathways Unique to Their Stage

We compared the transcriptomic profiles of two Rb subtypes and identified a set of significant genes belonging to the top enriched pathways in the Rb subtypes (p < 0.005). Among the previously identified genes in Rb [2], several were found to be differentially expressed in advanced and non-advanced Rb (E2F’s, CDK’s, MYCN, and SYK). We further arranged differentially expressed significant genes in advanced and non-advanced Rb into four gene groups showing their interaction networks (Figure 2A). Gene group 1 genes were significantly upregulated in advanced Rb tumors, compared to non-advanced Rb and controls (CDK1, CDKN2A, CCNB2, CCNE2, PTTG1, and DDB1). We performed KEGG enrichment analysis on gene group 1 using the gene sets from the gene ontology biological processes (GOBP) and the MSigDB hallmarks (HALLMARK), and revealed the gene interaction network of group 1 genes (Figure 2B) and their functional overlapping pathways (p < 0.005). Group 1 genes were associated with pathways such as cell cycle regulation, mitotic spindle checkpoint, DNA damage and response, and cytokinesis (Figure 2C). Gene group 2 displayed a significant upregulation of E2F transcription activators (E2F1, E2F2, E2F3, and E2F6) [30] and tumor-specific transcription factors (MYCN [31], RUNX1 [32], and GABPB1 [33]) in advanced Rb. Strikingly, we found a significant downregulation of transcription factors such as FOXO3 and repressor E2Fs (E2F7 and E2F8) [30] in the advanced Rb cohort compared to non-advanced tumors and controls. The epigenetic regulators show a high degree of expression for SYK, PRDM1, and TK1 in advanced Rb, while relatively lower expression of these factors was detected in non-advanced Rb. Notably, advanced Rb displayed a low degree of expression of BRD4, DNMT3A, and MGMT, while non-advanced Rb displayed a high degree of expression of these epigenetic factors. Immune-related genes such as CD86, CD19, and CD36 were significantly upregulated in advanced Rb while displaying a low degree of expression of CD81 and CD163. KEGG enrichment analysis revealed gene interaction networks of group 2 genes (Figure 2D) and their overlapping pathways (p < 0.005). Gene group 2 shows an association with pathways including cellular senescence, B-cell receptor signaling, pathways in cancer, transcriptional misregulation, JAK/STAT signaling, E2F targets, and MYC targets (Figure 2E). Gene group 3 shows differential expression of metabolic genes related to glycolysis, Krebs cycle, and fatty acid metabolism. Notably, glycolytic genes displayed a significantly low degree of expression in advanced and non-advanced subjects compared to the controls (HK1, HK2, HK3, GCKR, SLC2A1, G6PC, ALDOC1, and ENO1). However, Krebs cycle-related genes were significantly altered across all Rb subjects compared to the controls. Advanced Rb subjects showed a significantly high degree of expression for fatty acid metabolism genes such as MYCLD and ACSL1, compared to non-advanced Rb tumors. Notably, the PPAR pathway-related genes (PPARGC1A, PPARA, and PPARD) displayed a low degree of expression in advanced and non-advanced Rb tumors compared to the controls. KEGG enrichment analysis revealed gene interaction networks of group 2 genes (Figure 2F) and their overlapping pathways (p < 0.005). Group 3 genes were associated with pathways such as the integration of energy metabolism, glutamine biosynthesis, the regulation of lipid metabolism, pyrimidine salvage, ETC, and glycolysis (Figure 2G). Gene group 4 shows the enrichment of phototransduction specific genes that were significantly low in advanced and non-advanced Rb tumors. Network analysis shows the interaction of RHO with CRX and NRL, and the functional pathways regulated by the Group 4 genes (Figure S1D). Thus, the analysis revealed several previously unreported gene expression clusters and the associations that are unique to the Rb tumor stage.

3.3. Quantitative Validation of Microarray-Identified Targets Confirm Distinct Transcriptomic Profiles in Rb Subtypes

We performed RT-PCR to further evaluate the expression patterns of microarray-identified targets in advanced and non-advanced Rb groups. We found a significant degree of upregulation of cell cycle checkpoint genes such as CDK1, CCNB2, CCNE2, CDKN2A, CDKN3, PTTG1, CHEK2, and RPA4 (p ≤ 0.05) in advanced Rb tumors compared to non-advanced Rb tumors (Figure 3A). Microarray-identified epigenetic targets including PRDM1, SYK, and TK1 were also found to be upregulated in an RT-PCR analysis of advanced Rb subjects compared to non-advanced tumors (Figure 3B). Key photoreceptor genes such as RLBP1, RDH12, CRABP1, NRL, SOX8, RHO, PAX6, and SAG were significantly low in advanced Rb tumors, compared to non-advanced Rb tumors in RT-PCR analysis (p ≤ 0.05), further confirming the gene expression pattern identified in the tumor microarray (Figure 3C). RT-PCR analysis further confirmed the microarray expression patterns of glycolytic genes such as HK1 and G6PC, which shows significant downregulation in advanced Rb compared to non-advanced Rb (p ≤ 0.05) (Figure 3D). RT-PCR analysis confirmed the microarray expression of the E2F transcription activator E2F2 to be significantly high in advanced Rb tumors, compared to non-advanced Rb tumors. The E2F transcription repressor E2F7 was significantly low in advanced Rb tumors, compared to non-advanced Rb (p ≤ 0.05), further confirming the microarray-identified expression pattern (Figure 3E). RT-PCR log2-normalized values of the validated set of genes associated with the above-mentioned pathways showed a strong correlation with log2-normalized values of the same set of genes identified using the microarray analysis (r = 0.61, p < 0.001) (Figure 3F).

3.4. Differentially Accumulated Metabolites Reveal the Enrichment of Key Metabolic Pathways in Rb Vitreous Humor

The transcriptomic profile of the Rb tumor prompted us to elucidate the differential metabolites in the vitreous humor samples of the same Rb subjects and controls using metabolomics analysis. The principal component analysis of the samples shows a marked difference between Rb subjects and controls (Figure 4A). Metabolomics analysis was conducted using the LC-MS positive and negative (HILIC and C18 columns) method for a broader coverage of metabolites. A total of 1190 metabolites were detected in LC-MS analysis and the data were normalized to reduce or eliminate the effects of total sample amount variation on the quantification of individual metabolites (Figure 4B and Figure S2A). We performed a pairwise analysis of the samples using the clustering analysis functionality of GeneSpring. The correlation analysis followed by clustering showed a strong relationship between advanced Rb subjects (P1–P5) compared to non-advanced Rb (P6–P9) and controls (Figure S2B). Statistical analysis of the metabolomics data with stringent filters (p < 0.05, FC > 2) revealed 350 differentially expressed metabolites in the Rb vitreous humor (Figure 4C). Hierarchical clustering analysis performed on differentially accumulated metabolites identified components of lipid metabolism (sphingosine and ceramide) and fatty acid metabolism (okadaic acid, methyl palmitate, decanoyl-CoA, phytanic acid, and hexacosanoic acid) to be upregulated in advanced Rb compared to the non-advanced group. Strikingly, we found a low expression of soraphen A, an ACC inhibitor [34], in the advanced group, further highlighting the role of long-chain fatty acid synthesis in high-risk Rb tumors. A subset of the amino acids, namely creatine, taurine, and isoleucine–histidine–lysine displayed varied degrees of regulation in Rb subtypes. However, key metabolites of glycolysis (pyruvic acid, lactobionic acid, and galactosamine) were significantly downregulated in advanced Rb (Figure 4D). KEGG pathway analysis showed the enrichment of metabolites associated with linoleic acid metabolism, taurine, and hypo-taurine metabolism, and tyrosine metabolism, etc., as the top targets in Rb vitreous. Notably, metabolites that were associated with the citrate cycle and fatty acid metabolism showed a higher degree of enrichment than glycolysis and pyruvate metabolism in Rb vitreous (Figure 4E). This analysis highlights the unique details of significantly altered tumor metabolism in Rb eyes compared to healthy ones.

3.5. Integrated Transcriptomic and Metabolome Analysis Reveals Gene–Metabolite Interaction Networks Associated with Retinoblastoma

To gain insights into the overlapping pathways in Rb, we performed an integrated network analysis of transcriptomics and metabolomics to identify unique functional relationships between gene expression and metabolism. Using the KEGG global metabolic network, we mapped differentially accumulated metabolites detected in Rb vitreous samples (Figure 5A). Further, the enrichment analysis module of MetaboAnalyst was used, which further confirmed the significant alteration of fatty acid metabolism (yellow edges), amino acid metabolism (violet edges), and the pentose phosphate pathway (green edges) in the global metabolome map of Rb (Figure 5A). The metabolome map also highlighted the relatively low activity of glycolysis (blue edges) and the low abundance of the pyruvate metabolite. We constructed a metabolite–gene interaction network by integrating the differentially expressed transcriptomic targets and the associated metabolites (Figure 5B). The interactive map highlights the low abundance of the pyruvic acid metabolite (p < 0.05, FC < 2) and the significant downregulation of glycolytic genes (p < 0.05, FC < 2). Similarly, a low abundance of valproic acid (p < 0.05, FC < 2) is linked to a significant upregulation of HDAC2 (Histone Deacetylase 2) gene expression (p < 0.05, FC > 2) in the interactome map. The HDAC2 gene plays a redundant role in lipid accumulation [35], further confirming the high fatty acid metabolic state in Rb. A relatively high expression of FN1 (fibronectin1) and BAAT (Bile Acid-CoA:Amino Acid N-Acyltransferase) was also highlighted with the low abundance of the taurine metabolite in the interactome map (Figure 5B). FN1 overexpression can trigger ER stress and facilitate lipid accumulation [36], while BAAT transcription in conjugation with glycine or taurine regulates cholesterol and phospholipid synthesis [37], favoring cancer growth. Highly abundant metabolites such as salicylic acid, creatine, and myristic acid (p < 0.05, FC > 2) were found to interact with downregulated metabolism-associated genes such as HIF1A (Hypoxia inducible factor 1 subunit alpha), PRKAA2 (Protein Kinase AMP-Activated Catalytic Subunit Alpha 2), and GUCA1A (Guanylate Cyclase Activator 1A) (p < 0.05, FC < 2). A low expression of HIF1A reduces glycolytic metabolism and enhances mitochondrial oxidative phosphorylation [38,39], further providing evidence for low glycolysis in Rb. PRKAA2 encodes the AMPKα protein, which functions as a metabolic sensor in many diseases, including cancer [40]. The GUCA1A gene encodes for the GCAP1 protein, which regulates N-terminal myristoylation and lipid modifications in photoreceptors [41]. A combined pathway analysis of transcriptomics and metabolomics identified glycolysis, AMPKα signaling, HIF-1 signaling, and fatty acid biosynthesis to be altered in Rb tumors (Figure S3A).

3.6. Validation of Molecular Signatures That Specify Clinical and Histopathological Grades of Retinoblastoma

To validate the multi-omics-identified targets in Rb, we first performed immunohistochemical staining on the same set of advanced and non-advanced Rb tumor tissues, using pediatric retina as the controls to evaluate the classical molecular signs. H and E staining confirmed densely packed tumor cells with little cytoplasm partly arranged in perivascular cuffs in advanced Rb (Figure 6A). In contrast, the non-advanced Rb tissues contained tumor cells surrounding cystic spaces with relatively more cytoplasm and smaller nuclei. Healthy pediatric retina shows clear and visible retinal layers in H and E staining. Advanced Rb subjects showed strong and consistent immunostaining positivity for the proliferation marker Ki67 and the transcription factor E2F1, while their signals were relatively weaker, though not significant in non-advanced Rb, and undetectable in the pediatric retina (Figure 6A,C). The total Rb expression, observed in the healthy retina, was not detectable in tumors using IHC (Figure 6B,C). We further profiled the microarray expression of Ki67, E2F1, RB1, RBL1, and RBL2 in Rb subjects based on their laterality and severity. Ki67 expression was higher in bilateral and unilateral Rb tumors, and this consistent expression trend was evident in advanced and non-advanced Rb subjects (Figure 6D). E2F1 expression was higher in the bilateral group compared to the unilateral (p < 0.005), and in advanced Rb compared to non-advanced Rb (p < 0.05), providing clues as to Rb-loss mediated transcriptional effects. RB1 expression was significantly low in unilateral Rb compared to bilateral, advanced, and non-advanced Rb groups (Figure 6F), while its pocket protein family gene RBL2 (p130) [42] was significantly low in the bilateral and unilateral Rb groups, and in the advanced and non-advanced Rb compared to the pediatric controls (Figure 6H). However, RBL1 (p107) was upregulated in Rb subjects, confirming its cell cycle compensation effects [43] compared to the controls (Figure 6G).

3.7. Validation of Multi-Omics Findings in Rb Subjects and the In Vitro Model

The combined analysis of transcriptomics and metabolomics revealed glycolysis and pathways in cancer to be significantly altered in Rb (p < 0.001). We corroborated these findings by validating the expression of their key candidate markers: HK1, which regulates glycolysis [24], and the E2F2 transcription factor associated with key cancer pathways [44]. We performed an immunohistochemical analysis on an additional cohort of advanced (n = 15) and non-advanced Rb tissues (n = 10), and detected strong and consistent immunostaining positivity for E2F2 (IHC score= 2–3), while the signals were undetectable in control pediatric retinae (n = 2, IHC score = 0). HK1 showed weak signals in advanced (n = 15) and non-advanced Rb (n = 10) tumor areas (IHC score = 0–1), while strongly staining the intact retina portions within non-advanced Rb tissues, thereby serving as an internal control. Notably, HK1 exhibited strong and consistent signals in the photoreceptor layer of healthy pediatric controls (n = 2; IHC score= 3) (Figure 7A,B). RT-PCR validations further confirmed the expression pattern of E2F2 and HK1 in advanced and non-advanced Rb tumors (Figure 7C). To elucidate how RB1 influences HK1 and E2F2 expression in the context of cancer, we overexpressed the wild-type RB1 in Rb null WERI-Rb1 retinoblastoma cells in vitro. The ectopic expression of RB1 induced HK1 protein and transcript in WERI-Rb1 cells (Figure 7D,E,G, p < 0.001), while reducing E2F2 transcript and protein levels (Figure 7D, F, p < 0.001). We further performed total mRNA profiling using RB1-null Y79 and RB1-overexpressed Y79 microarrays. We applied the standard QC procedure, and the data was normalized, leading to comparable distributions between all arrays (Figure S3B). Statistical analysis revealed a distinct cluster of expression profiles between RB1-null and RB1-overexpressed Y79 (Figure S3C). We identified 65 differentially expressed genes between RB1-null and RB1-overexpressed Y79 (p < 0.05, FC > 2). The E2F transcription factor genes (E2F1, E2F2) and cell cycle genes (CDK1, CCNB2) were significantly upregulated (p < 0.05, FC > 2), and the glycolytic genes were significantly downregulated (HK1, ENO1, and G6PC3) (p < 0.05, FC < 2) (Figure 7H). The transcriptomic analysis also revealed a high degree of expression for the aldolase isoenzymes ALDOA, ALDOB, and ALDOC in the Rb-null condition, while RB1 complementation showed a significant downregulation of aldolase isoenzymes (Figure 7I). KEGG enrichment pathway analysis revealed sphingolipid metabolism, branched-chain amino-acid metabolism, apoptosis, and the DNA repair pathway to be significantly altered between the RB1 null and RB1-complemented Y79 cells, indicating that the cell line microarray had significant similarities with the transcriptomic landscape of the Rb subjects (Figure 2).

4. Discussion

This integrated omics analysis of the Rb tumor transcriptome and the vitreous humor metabolome within the human Rb tumor subtypes using the pediatric retina as controls was performed to enrich our understanding of altered gene function with changes in tumor metabolism. Although there have recently been a few reports on alterations in the Rb tumor-associated metabolites [45,46], we found significantly more functional associations between the metabolic targets and the enriched pathways owing to the coordinated samples used for analysis and the use of age-matched pediatric controls. Although we have focused on intraocular advanced and non-advanced retinoblastoma tumors, our findings can be extended to other cancer systems with persistent RB1 mutations.
Rb is a complex disease with predominant genetic and epigenetic events [47], and it is important to focus on improving our understanding of the regulatory mechanisms that promote tumor growth and intra-tumoral heterogeneity. To accomplish this, we classified our Rb cohort based on the AJCC or IIRC guidelines as advanced and non-advanced Rb, and used age-matched pediatric controls. The different clinical and pathological features of the two subtypes identified the relevance of this classification. The limitation of our study was the relatively small sample size due to the prerequisite of different sample types that were collected at the same time from each eye in the investigation. However, the major strength of the study, in comparison with the previous publications, is the ability to integrate different omics data sets from tumors and vitreous humors from the same patients. Further, we validated the protein levels in FFPE sections of the same tumors, as well as an additional independent cohort of 25 Rb tissues. These data sets allowed us to discover new insights into the interactions between the data sets, and to reveal new insights into the disease mechanism that was lacking in previously reported single omics studies [48,49,50].
Using transcriptomic profiling, we found there was a high degree of expression of cell cycle regulators and DNA damage and response checkpoint genes in both Rb subtypes, a consequence of pRb loss in the tumors [51,52]. Moreover, we report on the significant enrichment of E2F pathway genes in advanced Rb compared to non-advanced Rb, highlighting major flaws in transcriptional events due to RB1 inactivation in high-risk tumors. In contrast to previous findings on E2F1 [53], we report on the significantly high expression of E2F2 in advanced Rb tumors, along with the low expression of the E2F repressors E2F7 and E2F8, compared to non-advanced Rb in our study, which is also supported by a recent study [54]. Various studies have highlighted MYCN expression in Rb tumors [55], while we show the presence of the differential expression of MYCN in Rb tumor subtypes. Furthermore, we report on the elevation of tumor-specific transcription factors such as RUNX1 and GABPB1, and the reduced expression of FOXO3 that is low in the advanced group. While these transcriptional alterations are expected from the literature on other high-risk pediatric tumors [56,57,58,59], their distinct transcriptional regulation in advanced Rb tumors is a new observation.
Previous studies have highlighted the role of epigenetic modifications and SYK in targeting Rb tumorigenesis [4]. In agreement with these findings, we report on the high expression of SYK in advanced Rb tumors compared to non-advanced Rb. In addition, we also discovered the differential expression of PRDM1 and TK1 in Rb subtypes, which are known epigenetic regulators in other cancers [60,61]. Our study also highlights the low expression of known canonical epigenetic readers such as BRD4, DNMT3A, and MGMT in advanced Rb tumors, highlighting a diverse epigenetic landscape in advanced tumors, extending our insights into Rb pathogenesis. These observations give support for the further investigation of broad epigenetic modulators for adjunctive treatments in advanced Rb. We further report on the differential expression of several immune cell surface receptors within the tumor tissues. B-cell maturation genes such as CD19 and its binding partner CD81 were significantly altered in the Rb subtypes, which have roles in regulating the immune response and receptor signaling [62]. We also observed high CD86 levels in the advanced tumor, which are supported by the presence of high CD86 B cell signatures in advanced stages of other cancers [63]. However, a high degree of CD80 expression, along with CD86 and increased transcript numbers of T-cell counter receptors (CD28 and CTLA4) (data not shown) in tumors also point out the costimulatory roles of CD86/CD80, which are known to induce T cell anergy [64]. The differential expression of these immune markers highlights a unique immunophenotype of the Rb tumors that requires further exploration to understand their roles in tumor growth and its progression of severity. We also observed a differential loss of visual cycle-associated genes among the clinical subgroups and compared to controls, aspects of which are supported by certain studies [65,66]. Both the Rb subtypes expressed the cone-specific markers RXR γ, CRX, OTX1, and TRβ2 (data not shown), compared to the pediatric retina, supporting the cone photoreceptor origin of cancer [67].
We found that advanced Rb was associated with a low expression of glycolytic genes and altered Krebs cycle with a higher propensity for fatty acid metabolism. Our findings are also supported by a recent study that highlighted the role of lipid metabolism in Rb tumors [45]. Furthermore, our global metabolomic profiling of Rb vitreous humor identified a variety of differentially accumulated metabolites that are distinct from the pediatric controls. Importantly, we observed a significantly low abundance of pyruvate in Rb vitreous humor, which directly correlated with the loss of glycolytic gene expression, particularly HK1. Pathway analysis further confirmed our findings, highlighting the low enrichment of glycolysis and pyruvate metabolism pathways compared to high fatty acid and lipid metabolism in advanced Rb. Therefore, the metabolomics findings in vitreous humor corroborated with the transcriptomics profile of the Rb tumor and revealed a significant distinction in metabolite abundance between the Rb subtypes. The data integration of transcriptomic and metabolomic profiles of Rb revealed altered pathways, including cellular metabolism and transcriptional machinery, and we identified HK1 and E2F2 as being representative genes that regulate distinct transcriptional networks in Rb tumors. Through immunohistochemistry, we show the distribution of HK1 and E2F2 proteins in Rb tissue subtypes. The RB1 complex binds to E2F2 promoters and represses its function during the development and differentiation programs [68], and we assume that the high expression of E2F2 in tumors is a direct consequence of RB1 inactivation. However, HK1, which controls glycolytic flux [69], has been found to be upregulated in other tumors [70] and is associated with oncogenes such as KRAS [71] mutants. In this aspect, we find Rb tumors to behave in a metabolically opposite manner to other types of solid tumors, thus providing a potentially unique therapeutic target. In addition, recent reports highlight the glycolytic control exerted by pRb proteins [72], and we hypothesize that the low enrichment of glycolysis in tumors is a direct consequence of the loss of pRb proteins. Therefore, our observation in Rb tumors is critical in terms of delineating the unique metabolic phenotypes of tumors of diverse origins. Using in vitro models, we show that RB1 complementation enhanced the glycolytic HK1 protein and its transcript, while reducing E2F2 levels, further confirming our findings in human tumors. The transcriptomic profiling of RB1-complemented Y79 shows differentially regulated cell cycle and metabolic pathways that corroborated with findings in Rb tumors.

5. Conclusions

In conclusion, the comprehensive transcriptomic profiling of the Rb tumor subtypes revealed differentially expressed genes belonging to cell cycle, metabolism, epigenetics, immune regulation, and phototransduction pathways. The global metabolomic analysis of Rb vitreous humor revealed metabolites that are associated with critical signaling pathways that are committed to the cancer state, which are otherwise not apparent. The integration of multi-omics data identified co-expression network modules within Rb clinical subtypes that may be crucial for prognostication and future drug development efforts.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cells11101668/s1, Figure S1: Identification of distinct molecular signatures in advanced and non-advanced Rb tumors from primary enucleated eyes; Figure S2: Differentially accumulated metabolites revealed enrichment of key metabolic pathways in Rb vitreous humor; Figure S3: Pathways associated with Rb and in vitro models; Figure S4: Statistical decision tree. Table S1: Clinical and histopathological details of samples used for validation; Table S2: RT-PCR primer details.

Author Contributions

Conceptualization, A.M., N.G. and A.G.; methodology, D.S., S.G., V.S.B., N.G. and A.G.; investigation, V.S.B., N.G. and A.G.; resources, R.S., A.M., G.D., S.S.C. and S.H.; data curation, V.S.B., D.S., S.G., N.G. and A.G.; writing—original draft preparation, V.S.B., R.K., D.S., A.P.N. and A.G.; writing—review and editing, V.S.B., A.M., R.L., N.K.V., V.A.B., N.G., S.H. and AG.; supervision, V.A.B., N.G., S.H. and A.G.; funding acquisition, N.G., R.S. and A.G. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the Narayana Nethralaya Foundation, Bangalore, India. The funders had no role in the study design, data collection, and analysis.

Institutional Review Board Statement

This cross-sectional study was approved by the Narayana Nethralaya institutional ethics committee (EC Ref. no: C/2013/03/02; 5 April 2013). Subject recruitment and sample collection procedures were conducted in accordance with the tenets of the Declaration of Helsinki.

Informed Consent Statement

Informed written consents were received from all legal guardians of the subjects before inclusion in the study.

Data Availability Statement

Data will be available on request.

Acknowledgments

The authors thank G. Kumarmanickavel and Swaminathan Sethu for their expertise and assistance throughout all aspects of our study. We gratefully acknowledge the contribution of the Retinoblastoma clinic at NN for their help and guidance during the study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lohmann, D.R.; Gallie, B.L. GeneReviews®; Adam, M.P., Ardinger, H.H., Pagon, R.A., Wallace, S.E., Bean, L.J.H., Grip, K.W., Mirzaa, G.M., Amemiya, A., Eds.; University of Washington: Seattle, WA, USA, 1993. [Google Scholar]
  2. Dimaras, H.; Corson, T.W.; Cobrinik, D.; White, A.; Zhao, J.; Munier, F.L.; Abramson, D.L.; Shields, C.L.; Chantada, G.L.; Njuguna, F.; et al. Retinoblastoma. Nat. Rev. Dis. Prim. 2015, 1, 15021. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Xu, X.L.; Singh, H.P.; Wang, L.; Qi, D.L.; Poulos, B.K.; Abramson, D.H.; Jhanwar, S.C.; Cobrinik, D. Rb suppresses human cone-precursor-derived retinoblastoma tumours. Nature 2014, 514, 385–388. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Zhang, J.; Benavente, C.A.; McEvoy, J.; Flores-Otero, J.; Ding, L.; Chen, X.; Ulyanov, A.; Wu, G.; Wilson, M.; Wang, J.; et al. A novel retinoblastoma therapy from genomic and epigenetic analyses. Nature. 2012, 481, 329–334. [Google Scholar] [CrossRef] [PubMed]
  5. Gupta, H.; Malaichamy, S.; Mallipatna, A.; Murugan, S.; Jeyabalan, N.; Suresh Babu, V.; Ghosh, A.; Ghosh, A.; Santhosh, S.; Seshagiri, S.; et al. Retinoblastoma genetics screening and clinical management. BMC Med. Genom. 2021, 14, 188. [Google Scholar] [CrossRef]
  6. Singh, L.; Kashyap, S. Update on pathology of retinoblastoma. Int. J. Ophthalmol. 2018, 11, 2011–2016. [Google Scholar]
  7. Kaliki, S.; Patel, A.; Iram, S.; Ramappa, G.; Mohamed, A.; Palkonda, V.A.R. Retinoblastoma in India: Clinical presentation and outcome in 1457 Patients (2074 Eyes). Retina 2019, 39, 379–391. [Google Scholar] [CrossRef]
  8. Wong, E.S.; Choy, R.W.; Zhang, Y.; Chu, W.K.; Chen, L.J.; Pang, C.P.; Yam, J.C. Global retinoblastoma survival and globe preservation: A systematic review and meta-analysis of associations with socioeconomic and health-care factors. Lancet Glob. Health 2022, 10, e380–e389. [Google Scholar] [CrossRef]
  9. Fabian, I.D.; Onadim, Z.; Karaa, E.; Duncan, C.; Chowdhury, T.; Scheimberg, I.; Ohnuma, I.; Reddy, M.A.; Sagoo, M.S. The management of retinoblastoma. Oncogene 2018, 37, 1551–1560. [Google Scholar] [CrossRef]
  10. Berry, J.L.; Kogachi, K.; Aziz, H.A.; McGovern, K.; Zolfaghari, E.; Murphree, A.L.; Jubran, R.; Kim, J.W. Risk of metastasis and orbital recurrence in advanced retinoblastoma eyes treated with systemic chemoreduction versus primary enucleation. Pediatr. Blood Cancer 2017, 64, e26270. [Google Scholar] [CrossRef]
  11. Canturk, S.; Qaddoumi, I.; Khetan, V.; Ma, Z.; Furmanchuk, A.; Antoneli, C.B.; Sultan, I.; Kebudi, R.; Sharma, T.; Galindo, C.R.; et al. Survival of retinoblastoma in less-developed countries impact of socioeconomic and health-related indicators. Br. J. Ophthalmol. 2010, 94, 1432–1436. [Google Scholar]
  12. Zeng, Y.; He, T.; Liu, J.; Li, Z.; Xie, F.; Chen, C.; Xing, Y. Bioinformatics analysis of multi-omics data identifying molecular biomarker candidates and epigenetically regulatory targets associated with retinoblastoma. Medicine 2020, 99, e23314. [Google Scholar] [CrossRef] [PubMed]
  13. Zhao, Y.; Gao, Y.; Xu, X.; Zhou, J.; Wang, H. Multi-omics analysis of genomics, epigenomics and transcriptomics for molecular subtypes and core genes for lung adenocarcinoma. BMC Cancer 2021, 21, 257. [Google Scholar] [CrossRef] [PubMed]
  14. Swiss, V.A.; Casaccia, P. Cell-context specific role of the E2F/Rb pathway in development and disease. Glia 2010, 58, 377–390. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Barnum, K.J.; O’Connell, M.J. Cell cycle regulation by checkpoints. Methods Mol. Biol. 2014, 1170, 29–40. [Google Scholar]
  16. Kin, N.W.; Sanders, V.M. CD86 regulates IgG1 production via a CD19-dependent mechanism. J. Immunol. 2007, 179, 1516–1523. [Google Scholar] [CrossRef] [Green Version]
  17. McEvoy, J.D.; Dyer, M.A. Genetic and Epigenetic Discoveries in Human Retinoblastoma. Crit. Rev. Oncog. 2015, 20, 217–225. [Google Scholar]
  18. Liu, C.; Banister, C.E.; Weige, C.C.; Altomare, D.; Richardson, J.H.; Contreras, C.M.; Buckhaults, P.J. PRDM1 silences stem cell-related genes and inhibits proliferation of human colon tumor organoids. Proc. Natl. Acad. Sci. USA 2018, 115, E5066–E5075. [Google Scholar] [CrossRef] [Green Version]
  19. Etemad, S.; Petit, M.; Weiss, A.K.H.; Schrattenholz, A.; Baraldo, G.; Jansen-Durr, P. Oxaloacetate decarboxylase FAHD1—A new regulator of mitochondrial function and senescence. Mech. Ageing Dev. 2019, 177, 22–29. [Google Scholar] [CrossRef]
  20. Fulton, A.B.; Dodge, J.; Hansen, R.M.; Williams, T.P. The rhodopsin content of human eyes. Investig. Ophthalmol. Vis. Sci. 1999, 40, 1878–1883. [Google Scholar]
  21. Mears, A.J.; Kondo, M.; Swain, P.K.; Takada, Y.; Bush, R.A.; Saunders, T.L.; Sieving, P.A.; Swaroop, A. Nrl is required for rod photoreceptor development. Nat. Genet. 2001, 29, 447–452. [Google Scholar] [CrossRef]
  22. Cote, R.H. Characteristics of photoreceptor PDE (PDE6): Similarities and differences to PDE5. Int. J. Impot. Res. 2004, 16, S28–S33. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Napoli, J.L. Cellular retinoid binding-proteins, CRBP, CRABP, FABP5: Effects on retinoid metabolism, function and related diseases. Pharmacol. Ther. 2017, 173, 19–33. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Tseng, P.L.; Chen, C.W.; Hu, K.H.; Cheng, H.C.; Lin, Y.H.; Tsai, W.H.; Cheng, T.J.; Wu, W.H.; Yeh, C.W.; Lin, C.C.; et al. The decrease of glycolytic enzyme hexokinase 1 accelerates tumor malignancy via deregulating energy metabolism but sensitizes cancer cells to 2-deoxyglucose inhibition. Oncotarget 2018, 9, 18949–18969. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Lu, Y.-Y.; Wu, C.-H.; Hong, C.-H.; Chang, K.-L.; Lee, C.-H. GLUT-1 Enhances Glycolysis, Oxidative Stress, and Fibroblast Proliferation in Keloid. Life 2021, 11, 505. [Google Scholar] [CrossRef]
  26. Khatri, S.; Yepiskoposyan, H.; Gallo, C.A.; Tandon, P.; Plas, D.R. FOXO3a regulates glycolysis via transcriptional control of tumor suppressor TSC1. J. Biol. Chem. 2010, 285, 15960–15965. [Google Scholar]
  27. Christmann, M.; Verbeek, B.; Roos, W.P.; Kaina, B. O(6)-Methylguanine-DNA methyltransferase (MGMT) in normal tissues and tumors: Enzyme activity, promoter methylation and immunohistochemistry. Biochim. Biophys. Acta 2011, 1816, 179–190. [Google Scholar]
  28. Mallipatna, A.; Gallie, B.; Chévez-Barrios, P.; Rouic, L.; Chantada, G.; Doz, F.; Brisse, H.J.; Munier, F.L.; Albert, D.M.; Mora, J.C.; et al. AJCC Staging for Retinoblastoma, 8th ed.; Springer: Berlin, Germany, 2017; pp. 819–831. [Google Scholar]
  29. Linn Murphree, A. Intraocular retinoblastoma: The case for a new group classification. Ophthalmol. Clin. North Am. 2005, 18, 41–53. [Google Scholar] [CrossRef]
  30. Cuitino, M.C.; Pecot, T.; Sun, D.; Kladney, R.; Okano-Uchida, T.; Shinde, N.; Saeed, R.; Perez Castro, A.J.; Webb, A.; Liu, T.; et al. Two Distinct E2F Transcriptional Modules Drive Cell Cycles and Differentiation. Cell Rep. 2019, 27, 3547–3560.e5. [Google Scholar]
  31. Liu, Z.; Chen, S.S.; Clarke, S.; Veschi, V.; Thiele, C.J. Targeting MYCN in Pediatric and Adult Cancers. Front. Oncol. 2020, 10, 623679. [Google Scholar] [CrossRef]
  32. Hong, D.; Fritz, A.J.; Gordon, J.A.; Tye, C.E.; Boyd, J.R.; Tracy, K.M.; Frietze, S.E.; Carr, F.E.; Nickerson, J.A.; van Wijnen, A.J.; et al. RUNX1-dependent mechanisms in biological control and dysregulation in cancer. J. Cell Physiol. 2019, 234, 8597–8609. [Google Scholar]
  33. Bell, R.J.; Rube, H.T.; Kreig, A.; Mancini, A.; Fouse, S.D.; Nagarajan, R.P.; Choi, S.; Hong, C.; He, D.; Pekmezci, M.; et al. Cancer. The transcription factor GABP selectively binds and activates the mutant TERT promoter in cancer. Science 2015, 348, 1036–1039. [Google Scholar] [PubMed] [Green Version]
  34. Jump, D.B.; Torres-Gonzalez, M.; Olson, L.K. Soraphen A, an inhibitor of acetyl CoA carboxylase activity, interferes with fatty acid elongation. Biochem. Pharmacol. 2011, 81, 649–660. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Haberland, M.; Carrer, M.; Mokalled, M.H.; Montgomery, R.L.; Olson, E.N. Redundant control of adipogenesis by histone deacetylases 1 and 2. J. Biol. Chem. 2010, 285, 14663–14670. [Google Scholar] [PubMed] [Green Version]
  36. Du, H.; Wang, Y.; Zhang, Z.; Yang, J.; Zhang, J.; Zhang, Y. Fibronectin Overexpression Modulates Formation of Macrophage Foam Cells by Activating SREBP2 Involved in Endoplasmic Reticulum Stress. Cell Physiol. Biochem. 2015, 36, 1821–1834. [Google Scholar] [CrossRef]
  37. Chiang, J.Y. Bile acid metabolism and signaling. Compr. Physiol. 2013, 3, 1191–1212. [Google Scholar] [PubMed] [Green Version]
  38. Del Rey, M.J.; Valin, A.; Usategui, A.; Garcia-Herrero, C.M.; Sanchez-Arago, M.; Cuezva, J.M.; Galindo, M.; Bravo, B.; Canete, J.D.; Blanco, F.J.; et al. Hif-1alpha knockdown reduces glycolytic metabolism and induces cell death of human synovial fibroblasts under normoxic conditions. Sci. Rep. 2017, 7, 3644. [Google Scholar] [CrossRef] [Green Version]
  39. Miska, J.; Lee-Chang, C.; Rashidi, A.; Muroski, M.E.; Chang, A.L.; Lopez-Rosas, A.; Zhang, P.; Panek, W.K.; Cordero, A.; Han, Y.; et al. HIF-1alpha is a metabolic switch between glycolytic-driven migration and oxidative phosphorylation-driven immunosuppression of Tregs in glioblastoma. Cell Rep. 2019, 27, 226–237.e4. [Google Scholar]
  40. Herzig, S.; Shaw, R.J. AMPK: Guardian of metabolism and mitochondrial homeostasis. Nat. Rev. Mol. Cell Biol. 2018, 19, 121–135. [Google Scholar] [CrossRef] [Green Version]
  41. Hwang, J.Y.; Koch, K.W. Calcium- and myristoyl-dependent properties of guanylate cyclase-activating protein-1 and protein-2. Biochemistry 2002, 41, 13021–13028. [Google Scholar] [CrossRef]
  42. Giacinti, C.; Giordano, A. RB and cell cycle progression. Oncogene 2006, 25, 5220–5227. [Google Scholar]
  43. Sage, J.; Miller, A.L.; Perez-Mancera, P.A.; Wysocki, J.M.; Jacks, T. Acute mutation of retinoblastoma gene function is sufficient for cell cycle re-entry. Nature 2003, 424, 223–228. [Google Scholar] [CrossRef] [PubMed]
  44. Zeng, Z.; Cao, Z.; Tang, Y. Increased E2F2 predicts poor prognosis in patients with HCC based on TCGA data. BMC Cancer 2020, 20, 1037. [Google Scholar] [CrossRef] [PubMed]
  45. Sahoo, S.; Ravi Kumar, R.K.; Nicolay, B.; Mohite, O.; Sivaraman, K.; Khetan, V.; Rishi, P.; Ganesan, S.; Subramanyam, K.; Raman, K.; et al. Metabolite systems profiling identifies exploitable weaknesses in retinoblastoma. FEBS Lett. 2019, 593, 23–41. [Google Scholar] [CrossRef] [PubMed]
  46. Kohe, S.; Brundler, M.A.; Jenkinson, H.; Parulekar, M.; Wilson, M.; Peet, A.C.; McConville, C.M. Metabolite profiling in retinoblastoma identifies novel clinicopathological subgroups. Br. J. Cancer 2015, 113, 1216–1224. [Google Scholar] [CrossRef] [Green Version]
  47. Benavente, C.A.; Dyer, M.A. Genetics and epigenetics of human retinoblastoma. Annu. Rev. Pathol. 2015, 10, 547–562. [Google Scholar] [CrossRef] [Green Version]
  48. Elchuri, S.V.; Rajasekaran, S.; Miles, W.O. RNA-Sequencing of Primary Retinoblastoma Tumors Provides New Insights and Challenges Into Tumor Development. Front. Genet. 2018, 9, 170. [Google Scholar] [CrossRef] [Green Version]
  49. Chakraborty, S.; Khare, S.; Dorairaj, S.K.; Prabhakaran, V.C.; Prakash, D.R.; Kumar, A. Identification of genes associated with tumorigenesis of retinoblastoma by microarray analysis. Genomics 2007, 90, 344–353. [Google Scholar] [CrossRef]
  50. Galardi, A.; Colletti, M.; Lavarello, C.; Di Paolo, V.; Mascio, P.; Russo, I.; Cozza, R.; Romanzo, A.; Valente, P.; De Vito, R.; et al. Proteomic Profiling of Retinoblastoma-Derived Exosomes Reveals Potential Biomarkers of Vitreous Seeding. Cancers 2020, 12, 1555. [Google Scholar] [CrossRef]
  51. Kent, L.N.; Leone, G. The broken cycle: E2F dysfunction in cancer. Nat. Rev. Cancer 2019, 19, 326–338. [Google Scholar] [CrossRef]
  52. Shamma, A.; Takegami, Y.; Miki, T.; Kitajima, S.; Noda, M.; Obara, T.; Okamoto, T.; Takahashi, C. Rb Regulates DNA damage response and cellular senescence through E2F-dependent suppression of N-ras isoprenylation. Cancer Cell 2009, 15, 255–269. [Google Scholar] [CrossRef] [Green Version]
  53. Sangwan, M.; McCurdy, S.R.; Livne-Bar, I.; Ahmad, M.; Wrana, J.L.; Chen, D.; Bremner, R. Established and new mouse models reveal E2f1 and Cdk2 dependency of retinoblastoma, and expose effective strategies to block tumor initiation. Oncogene 2012, 31, 5019–5028. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Fan, X.; Wang, Y.; Jiang, T.; Liu, T.; Jin, Y.; Du, K.; Niu, Y.; Zhang, C.; Liu, Z.; Lei, Y.; et al. B-Myb accelerates colorectal cancer progression through reciprocal feed-forward transactivation of E2F2. Oncogene 2021, 40, 5613–5625. [Google Scholar] [CrossRef] [PubMed]
  55. Rushlow, D.E.; Mol, B.M.; Kennett, J.Y.; Yee, S.; Pajovic, S.; Theriault, B.L.; Progoda Lee, N.L.; Spencer, C.; Dimaras, H.; Corson, T.W.; et al. Characterisation of retinoblastomas without RB1 mutations: Genomic, gene expression, and clinical studies. Lancet Oncol. 2013, 14, 327–334. [Google Scholar] [CrossRef] [Green Version]
  56. Huang, M.; Weiss, W.A. Neuroblastoma and MYCN. Cold Spring Harb. Perspect. Med. 2013, 3, a014415. [Google Scholar] [CrossRef] [PubMed]
  57. Utnes, P.; Lokke, C.; Flaegstad, T.; Einvik, C. Clinically Relevant Biomarker Discovery in High-Risk Recurrent Neuroblastoma. Cancer Inform. 2019, 18, 1176935119832910. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Yang, B.; Dai, J.X.; Pan, Y.B.; Ma, Y.B.; Chu, S.H. Examining the biomarkers and molecular mechanisms of medulloblastoma based on bioinformatics analysis. Oncol. Lett. 2019, 18, 433–441. [Google Scholar] [CrossRef] [PubMed]
  59. Rupp, M.; Hagenbuchner, J.; Rass, B.; Fiegl, H.; Kiechl-Kohlendorfer, U.; Obexer, P.; Ausserlechner, M.J. FOXO3-mediated chemo-protection in high-stage neuroblastoma depends on wild-type TP53 and SESN3. Oncogene 2017, 36, 6190–6203. [Google Scholar] [CrossRef] [Green Version]
  60. Romero-Garcia, R.; Gomez-Jaramillo, L.; Mateos, R.M.; Jimenez-Gomez, G.; Pedreno-Horrillo, N.; Foncubierta, E.; Rodriguez-Gutierrez, J.F.; Garzon, S.; Lopez, F.M.; Rodriguez, C.; et al. Differential epigenetic regulation between the alternative promoters, PRDM1alpha and PRDM1beta, of the tumour suppressor gene PRDM1 in human multiple myeloma cells. Sci. Rep. 2020, 10, 15899. [Google Scholar] [CrossRef]
  61. Malvi, P.; Janostiak, R.; Nagarajan, A.; Cai, G.; Wajapeyee, N. Loss of thymidine kinase 1 inhibits lung cancer growth and metastatic attributes by reducing GDF15 expression. PLoS Genet. 2019, 15, e1008439. [Google Scholar] [CrossRef] [Green Version]
  62. Susa, K.J.; Seegar, T.C.; Blacklow, S.C.; Kruse, A.C. A dynamic interaction between CD19 and the tetraspanin CD81 controls B cell co-receptor trafficking. eLife 2020, 9, e52337. [Google Scholar] [CrossRef]
  63. Schlosser, H.A.; Thelen, M.; Lechner, A.; Wennhold, K.; Garcia-Marquez, M.A.; Rothschild, S.I.; Staib, E.; Zander, T.; Beutner, D.; Gathof, B.; et al. B cells in esophago-gastric adenocarcinoma are highly differentiated, organize in tertiary lymphoid structures and produce tumor-specific antibodies. Oncoimmunology 2019, 8, e1512458. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Halliday, N.; Williams, C.; Kennedy, A.; Waters, E.; Pesenacker, A.M.; Soskic, B.; Hinze, C.; Hou, T.Z.; Rowshanravan, B.; Janman, D.; et al. CD86 Is a Selective CD28 Ligand Supporting FoxP3+ Regulatory T Cell Homeostasis in the Presence of High Levels of CTLA-4. Front. Immunol. 2020, 11, 600000. [Google Scholar] [CrossRef] [PubMed]
  65. Liu, J.; Ottaviani, D.; Sefta, M.; Desbrousses, C.; Chapeaublanc, E.; Aschero, R. A high-risk retinoblastoma subtype with stemness features, dedifferentiated cone states and neuronal/ganglion cell gene expression. Nat. Commun. 2021, 12, 5578. [Google Scholar] [CrossRef] [PubMed]
  66. Kapatai, G.; Brundler, M.A.; Jenkinson, H.; Kearns, P.; Parulekar, M.; Peet, A.C.; McConville, C.M. Gene expression profiling identifies different sub-types of retinoblastoma. Br. J. Cancer 2013, 109, 512–525. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Xu, X.L.; Fang, Y.; Lee, T.C.; Forrest, D.; Gregory-Evans, C.; Almeida, D.; Liu, A.; Jhanwar, S.C.; Abramson, D.H.; Cobrinik, D. Retinoblastoma has properties of a cone precursor tumor and depends upon cone-specific MDM2 signaling. Cell 2009, 137, 1018–1031. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Sun, H.; Chang, Y.; Schweers, B.; Dyer, M.A.; Zhang, X.; Hayward, S.W.; Goodrich, D.W. An E2F binding-deficient Rb1 protein partially rescues developmental defects associated with Rb1 nullizygosity. Mol. Cell Biol. 2006, 26, 1527–1537. [Google Scholar] [CrossRef] [Green Version]
  69. Tanner, L.B.; Goglia, A.G.; Wei, M.H.; Sehgal, T.; Parsons, L.R.; Park, J.O.; White, E.; Toettcher, J.E.; Rabinowitz, J.D. Four Key Steps Control Glycolytic Flux in Mammalian Cells. Cell Syst. 2018, 7, 49–62.e8. [Google Scholar] [CrossRef]
  70. He, X.; Lin, X.; Cai, M.; Zheng, X.; Lian, L.; Fan, D.; Wu, X.; Lan, P.; Wang, J. Overexpression of Hexokinase 1 as a poor prognosticator in human colorectal cancer. Tumor Biol. 2016, 37, 3887–3895. [Google Scholar] [CrossRef]
  71. Amendola, C.R.; Mahaffey, J.P.; Parker, S.J.; Ahearn, I.M.; Chen, W.-C.; Zhou, M.; Court, H.; Shi, J.; Mendoza, S.L.; Morten, M.J.; et al. KRAS4A directly regulates hexokinase 1. Nature 2019, 576, 482–486. [Google Scholar] [CrossRef]
  72. Clem, B.F.; Chesney, J. Molecular pathways: Regulation of metabolism by RB. Clin. Cancer Res. 2012, 18, 6096–6100. [Google Scholar] [CrossRef] [Green Version]
Cells 11 01668 g001 550
Figure 1. Identification of distinct transcriptional signatures in advanced and non-advanced Rb tumors from primary enucleated eyes. (A) Volcano plot of differentially expressed genes. Density curve plot showing the overall signal distribution of all probe sets on the microarray. (B) Before normalization. (C) After normalization. Spearman’s correlation analysis of genes detected in microarray platform across: (D) Control vs. Advanced Rb (r = 0.42, p < 0.001); (E) Control vs. Non-advanced Rb (r = 0.90, p < 0.001. (F) Principal component analysis on the microarray results shows distinct clusters, implicating gross differences between the advanced Rb, non-advanced Rb, and the controls. (G) Venn diagram showing differential genes identified in Group E (advanced) and Group D (non-advanced) Rb tumor microarrays compared to pediatric retina.
Figure 1. Identification of distinct transcriptional signatures in advanced and non-advanced Rb tumors from primary enucleated eyes. (A) Volcano plot of differentially expressed genes. Density curve plot showing the overall signal distribution of all probe sets on the microarray. (B) Before normalization. (C) After normalization. Spearman’s correlation analysis of genes detected in microarray platform across: (D) Control vs. Advanced Rb (r = 0.42, p < 0.001); (E) Control vs. Non-advanced Rb (r = 0.90, p < 0.001. (F) Principal component analysis on the microarray results shows distinct clusters, implicating gross differences between the advanced Rb, non-advanced Rb, and the controls. (G) Venn diagram showing differential genes identified in Group E (advanced) and Group D (non-advanced) Rb tumor microarrays compared to pediatric retina.
Cells 11 01668 g001
Cells 11 01668 g002 550
Figure 2. Distinct genes and biological processes are altered in advanced and non-advanced retinoblastoma. (A) Heatmap showing significantly differentially expressed genes identified four main gene groups in advanced (n = 5) and non-advanced Rb (n = 4), compared to pediatric retina controls (n = 2). Genes were grouped based on the pathways reported in GOBP. Results are represented as a network of enriched genes sets (nodes) connected to their overlapping pathways (edges). The node size is proportional to the total number of genes in the gene set of interest. (B) Gene group 1 represents an enriched cluster of genes in the cell cycle, and the DNA damage and response pathway. (C) List of functional pathways regulated by the gene sets in Group 1. (D) Gene group 2 represents an enriched group of genes belonging to transcriptional regulation, immune system regulation, and epigenetic factors. (E) List of functional pathways regulated by the gene sets in group 2. (F) Gene group 3 represents an enriched set of genes belonging to cellular metabolic processes such as glycolysis, Krebs cycle, and fatty acid metabolism. (G) List of functional pathways regulated by the gene sets in group 3.
Figure 2. Distinct genes and biological processes are altered in advanced and non-advanced retinoblastoma. (A) Heatmap showing significantly differentially expressed genes identified four main gene groups in advanced (n = 5) and non-advanced Rb (n = 4), compared to pediatric retina controls (n = 2). Genes were grouped based on the pathways reported in GOBP. Results are represented as a network of enriched genes sets (nodes) connected to their overlapping pathways (edges). The node size is proportional to the total number of genes in the gene set of interest. (B) Gene group 1 represents an enriched cluster of genes in the cell cycle, and the DNA damage and response pathway. (C) List of functional pathways regulated by the gene sets in Group 1. (D) Gene group 2 represents an enriched group of genes belonging to transcriptional regulation, immune system regulation, and epigenetic factors. (E) List of functional pathways regulated by the gene sets in group 2. (F) Gene group 3 represents an enriched set of genes belonging to cellular metabolic processes such as glycolysis, Krebs cycle, and fatty acid metabolism. (G) List of functional pathways regulated by the gene sets in group 3.
Cells 11 01668 g002
Cells 11 01668 g003 550
Figure 3. Validation of selected genes for mRNA expression in retinoblastoma subtypes. RT-PCR validation of microarray-identified targets (A) Cell cycle checkpoint genes. (B) Epigenetic regulators. (C) Photoreceptor-specific genes. (D) Energy and growth regulator genes. (E) Transcription factor genes. (F) Pearson correlation plot between microarray and RT-PCR expression of genes involved in cell cycle checkpoint, epigenetic regulation, photoreceptor-specific, energy regulation, and transcription factors. Values represent mean ± SEM. Unpaired two-sided Student’s t-test was used for statistical analysis. * p < 0.05, ns = no significance.
Figure 3. Validation of selected genes for mRNA expression in retinoblastoma subtypes. RT-PCR validation of microarray-identified targets (A) Cell cycle checkpoint genes. (B) Epigenetic regulators. (C) Photoreceptor-specific genes. (D) Energy and growth regulator genes. (E) Transcription factor genes. (F) Pearson correlation plot between microarray and RT-PCR expression of genes involved in cell cycle checkpoint, epigenetic regulation, photoreceptor-specific, energy regulation, and transcription factors. Values represent mean ± SEM. Unpaired two-sided Student’s t-test was used for statistical analysis. * p < 0.05, ns = no significance.
Cells 11 01668 g003
Cells 11 01668 g004 550
Figure 4. Differentially accumulated metabolites revealed the enrichment of key metabolic pathways in Rb vitreous humor. (A) Principal component analysis of metabolite profiling in the vitreous humor of Rb and pediatric controls. (B) Density plot showing the distribution of metabolite intensity before and after normalization. (C) Volcano plot of differentially accumulated metabolites. p < 0.05 (in red), FC > 2: significantly increased; p < 0.05, FC < 2 (in blue): significantly reduced relative to the control. (D) Hierarchical clustering heatmap of significant metabolites in the vitreous humor of Rb and controls using Pearson correlation as the metric distance. The color code indicates the metabolite’s abundance. (E) KEGG pathway enrichment analysis of differentially accumulated metabolites.
Figure 4. Differentially accumulated metabolites revealed the enrichment of key metabolic pathways in Rb vitreous humor. (A) Principal component analysis of metabolite profiling in the vitreous humor of Rb and pediatric controls. (B) Density plot showing the distribution of metabolite intensity before and after normalization. (C) Volcano plot of differentially accumulated metabolites. p < 0.05 (in red), FC > 2: significantly increased; p < 0.05, FC < 2 (in blue): significantly reduced relative to the control. (D) Hierarchical clustering heatmap of significant metabolites in the vitreous humor of Rb and controls using Pearson correlation as the metric distance. The color code indicates the metabolite’s abundance. (E) KEGG pathway enrichment analysis of differentially accumulated metabolites.
Cells 11 01668 g004
Cells 11 01668 g005 550
Figure 5. Integrated transcriptomic and metabolome analysis provides unique insights into pathways that are associated with retinoblastoma. (A) Differentially accumulated metabolites identified in Rb vitreous were mapped onto the KEGG global metabolome map associated with Rb. The nodes in the figure represent metabolic compounds. Edges are enzymatic transformations. Highlighted nodes represent the differentially accumulated metabolites identified in Rb vitreous. p < 0.05, FC > 2; upregulated (in red), p < 0.05, FC < 2; downregulated (in green). Edges highlighted in violet represent the glycolysis/gluconeogenesis map in Rb. Edges highlighted in yellow represent the fatty acid metabolism-related enzyme map in Rb. Edges in pink represent the amino-acid metabolism map in Rb. Edges in green represent the pentose phosphate pathway enzymes in Rb. (B) Integration of transcriptomic and metabolomic interaction network. Circular nodes indicate differentially regulated genes identified from the microarray; (in green as downregulated set: p < 0.05, FC < 2; in red as upregulated set: p < 0.05, FC > 2. Square blue nodes indicate differentially accumulated metabolites from metabolomics analysis (p < 0.05). Blue square nodes indicate upregulated metabolites (FC > 2), and yellow square nodes indicate downregulated metabolites (FC < 2).
Figure 5. Integrated transcriptomic and metabolome analysis provides unique insights into pathways that are associated with retinoblastoma. (A) Differentially accumulated metabolites identified in Rb vitreous were mapped onto the KEGG global metabolome map associated with Rb. The nodes in the figure represent metabolic compounds. Edges are enzymatic transformations. Highlighted nodes represent the differentially accumulated metabolites identified in Rb vitreous. p < 0.05, FC > 2; upregulated (in red), p < 0.05, FC < 2; downregulated (in green). Edges highlighted in violet represent the glycolysis/gluconeogenesis map in Rb. Edges highlighted in yellow represent the fatty acid metabolism-related enzyme map in Rb. Edges in pink represent the amino-acid metabolism map in Rb. Edges in green represent the pentose phosphate pathway enzymes in Rb. (B) Integration of transcriptomic and metabolomic interaction network. Circular nodes indicate differentially regulated genes identified from the microarray; (in green as downregulated set: p < 0.05, FC < 2; in red as upregulated set: p < 0.05, FC > 2. Square blue nodes indicate differentially accumulated metabolites from metabolomics analysis (p < 0.05). Blue square nodes indicate upregulated metabolites (FC > 2), and yellow square nodes indicate downregulated metabolites (FC < 2).
Cells 11 01668 g005
Cells 11 01668 g006 550
Figure 6. Molecular signatures specific to the clinical and histopathological grades of retinoblastoma in primary enucleated eyes. (A) H and E and immunohistochemistry profile showing the expression of Ki67, Rb, and E2F1 on 4 Rb tissues. High levels of expression for Ki67 and E2F1 were observed in the advanced Rb group, while expression of both proteins were lower in the non-advanced Rb group. Scalebar = 50 µm. (B) H and E and immunohistochemistry profile showing expression of Ki67, E2F1 on 2 pediatric retina tissues. Rb expression was high in pediatric control tissues while no E2F1 expression was observed. Scalebar = 50 µm. (C) IHC scores of Ki67, Rb, and E2F1 in 9 Rb tissues and 2 pediatric retina tissues. (D) Gene expression of Ki67 in bilateral, unilateral, Group E, and Group D identified in microarray compared to control. Ki67 expression was high in the Rb tumor microarray. (E) Gene expression of E2F1 in bilateral, unilateral, Group E, and Group D identified in microarray compared to control. E2F1 expression was significantly higher in Bilateral and Group E Rb subjects compared to unilateral and Group D Rb subjects. (F) Gene expression of RB1 in bilateral, unilateral, Group E, and Group D identified in microarray compared to control. RB1 expression was found to be significantly low in both Group E and Group D compared to pediatric controls. RB1 expression was low in the unilateral and bilateral groups compared to pediatric control. (G) Microarray expression of RBL1 (p107) in Group E and Group D, unilateral and bilateral Rb subjects. RBL1 expression was found to be significantly high in both Group E and Group D compared to pediatric controls. RBL1 expression was high in the unilateral and bilateral group compared to pediatric control. (H) Microarray expression of RBL2 (p130) in Group E, Group D, unilateral, and bilateral Rb subjects. RBL2 expression was found to be significantly low in both Group E and Group D compared to pediatric controls. RBL2 expression was low in the unilateral and bilateral groups compared to pediatric control. Values represent mean ± SEM. Unpaired two-sided Student’s t-test was used for statistical analysis. * p < 0.05, ** p < 0.01, *** p < 0.001, ns = no significance.
Figure 6. Molecular signatures specific to the clinical and histopathological grades of retinoblastoma in primary enucleated eyes. (A) H and E and immunohistochemistry profile showing the expression of Ki67, Rb, and E2F1 on 4 Rb tissues. High levels of expression for Ki67 and E2F1 were observed in the advanced Rb group, while expression of both proteins were lower in the non-advanced Rb group. Scalebar = 50 µm. (B) H and E and immunohistochemistry profile showing expression of Ki67, E2F1 on 2 pediatric retina tissues. Rb expression was high in pediatric control tissues while no E2F1 expression was observed. Scalebar = 50 µm. (C) IHC scores of Ki67, Rb, and E2F1 in 9 Rb tissues and 2 pediatric retina tissues. (D) Gene expression of Ki67 in bilateral, unilateral, Group E, and Group D identified in microarray compared to control. Ki67 expression was high in the Rb tumor microarray. (E) Gene expression of E2F1 in bilateral, unilateral, Group E, and Group D identified in microarray compared to control. E2F1 expression was significantly higher in Bilateral and Group E Rb subjects compared to unilateral and Group D Rb subjects. (F) Gene expression of RB1 in bilateral, unilateral, Group E, and Group D identified in microarray compared to control. RB1 expression was found to be significantly low in both Group E and Group D compared to pediatric controls. RB1 expression was low in the unilateral and bilateral groups compared to pediatric control. (G) Microarray expression of RBL1 (p107) in Group E and Group D, unilateral and bilateral Rb subjects. RBL1 expression was found to be significantly high in both Group E and Group D compared to pediatric controls. RBL1 expression was high in the unilateral and bilateral group compared to pediatric control. (H) Microarray expression of RBL2 (p130) in Group E, Group D, unilateral, and bilateral Rb subjects. RBL2 expression was found to be significantly low in both Group E and Group D compared to pediatric controls. RBL2 expression was low in the unilateral and bilateral groups compared to pediatric control. Values represent mean ± SEM. Unpaired two-sided Student’s t-test was used for statistical analysis. * p < 0.05, ** p < 0.01, *** p < 0.001, ns = no significance.
Cells 11 01668 g006
Cells 11 01668 g007 550
Figure 7. Differentially expressed targets confirm multi-omics findings in enucleated Rb tissue samples and in vitro. (A) Immunohistochemistry profile showing the expression of phospho-Rb, E2F2, and HK1 on 4 Rb tissues and 2 pediatric retinae. Scalebar = 200 µm. (B) IHC scores of phospho-Rb, E2F2, and HK1 in 25 Rb tissues and 2 pediatric retina tissues. (C) Gene expression of E2F2 and HK1 in Rb subtypes relative to the pediatric retina. (D) Protein expression of phospho-Rb, Rb, E2F2, and HK1 in the control and RB1-complemented WERI-Rb1 cells. Relative gene expression of (E) RB1, (F) E2F2, and (G) HK1 in control and RB1-complemented WERI-Rb1 cells, which illustrate their differences at the molecular level. (H) Volcano plot of differentially regulated genes in RB1-null and RB1-complemented Y79 cells. (I) Heatmap showing significantly differentially expressed genes identified in the RB1-null and RB1-complemented Y79 cells. Values represent mean ± SEM. Unpaired two-sided Student’s t-test was used for statistical analysis. * p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001, ## p = 0.005.
Figure 7. Differentially expressed targets confirm multi-omics findings in enucleated Rb tissue samples and in vitro. (A) Immunohistochemistry profile showing the expression of phospho-Rb, E2F2, and HK1 on 4 Rb tissues and 2 pediatric retinae. Scalebar = 200 µm. (B) IHC scores of phospho-Rb, E2F2, and HK1 in 25 Rb tissues and 2 pediatric retina tissues. (C) Gene expression of E2F2 and HK1 in Rb subtypes relative to the pediatric retina. (D) Protein expression of phospho-Rb, Rb, E2F2, and HK1 in the control and RB1-complemented WERI-Rb1 cells. Relative gene expression of (E) RB1, (F) E2F2, and (G) HK1 in control and RB1-complemented WERI-Rb1 cells, which illustrate their differences at the molecular level. (H) Volcano plot of differentially regulated genes in RB1-null and RB1-complemented Y79 cells. (I) Heatmap showing significantly differentially expressed genes identified in the RB1-null and RB1-complemented Y79 cells. Values represent mean ± SEM. Unpaired two-sided Student’s t-test was used for statistical analysis. * p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001, ## p = 0.005.
Cells 11 01668 g007
Table
Table 1. Clinical and histopathological details of samples.
Table 1. Clinical and histopathological details of samples.
IDSexLateralityAge at PresentationClinical RiskIIRC GroupAJCC Staging
P1MBilateral15 monthsAdvancedGroup EcT3b
P2FUnilateral20 monthsAdvancedGroup EcT3b
P3MUnilateral24 monthsAdvancedGroup EcT3a
P4FBilateral4 monthsAdvancedGroup EcT3b
P5MBilateral30 monthsAdvancedGroup EcT3b
P6FBilateral21 monthsNon-advancedGroup DcT2b
P7FUnilateral28 monthsNon-advancedGroup DcT2b
P8MUnilateral20 monthsNon-advanced Group DcT2b
P9MUnilateral21 monthsNon-advancedGroup DcT2a
Control 1FNA3 monthsCardiac Arrest (no ocular complications)
Control 2FNA2 monthsMultiple organ dysfunction (no ocular complications)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Babu, V.S.; Mallipatna, A.; SA, D.; Dudeja, G.; Kannan, R.; Shetty, R.; Nair, A.P.; Gundimeda, S.; Chaurasia, S.S.; Verma, N.K.; et al. Integrated Analysis of Cancer Tissue and Vitreous Humor from Retinoblastoma Eyes Reveals Unique Tumor-Specific Metabolic and Cellular Pathways in Advanced and Non-Advanced Tumors. Cells 2022, 11, 1668. https://doi.org/10.3390/cells11101668

AMA Style

Babu VS, Mallipatna A, SA D, Dudeja G, Kannan R, Shetty R, Nair AP, Gundimeda S, Chaurasia SS, Verma NK, et al. Integrated Analysis of Cancer Tissue and Vitreous Humor from Retinoblastoma Eyes Reveals Unique Tumor-Specific Metabolic and Cellular Pathways in Advanced and Non-Advanced Tumors. Cells. 2022; 11(10):1668. https://doi.org/10.3390/cells11101668

Chicago/Turabian Style

Babu, Vishnu Suresh, Ashwin Mallipatna, Deepak SA, Gagan Dudeja, Ramaraj Kannan, Rohit Shetty, Archana Padmanabhan Nair, Seetharamanjaneyulu Gundimeda, Shyam S. Chaurasia, Navin Kumar Verma, and et al. 2022. "Integrated Analysis of Cancer Tissue and Vitreous Humor from Retinoblastoma Eyes Reveals Unique Tumor-Specific Metabolic and Cellular Pathways in Advanced and Non-Advanced Tumors" Cells 11, no. 10: 1668. https://doi.org/10.3390/cells11101668

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop