Heterogeneity in retinoblastoma: a tale of molecules and models

Recent progress in retinoblastoma: Pathogenesis, presentation, diagnosis and management

Retinoblastoma, the primary ocular malignancy in pediatric patients, poses a substantial threat to mortality without prompt and effective management. The prognosis for survival and preservation of visual acuity hinges upon the disease severity at the time of initial diagnosis. Notably, retinoblastoma has played a crucial role in unraveling the genetic foundations of oncogenesis. The process of tumorigenesis commonly begins with the occurrence of biallelic mutation in the RB1 tumor suppressor gene, which is then followed by a cascade of genetic and epigenetic alterations that correspond to the clinical stage and pathological features of the tumor. The RB1 gene, recognized as a tumor suppressor, encodes the retinoblastoma protein, which plays a vital role in governing cellular replication through interactions with E2F transcription factors and chromatin remodeling proteins. The diagnosis and treatment of retinoblastoma necessitate consideration of numerous factors, including disease staging, germline mutation status, family psychosocial factors, and the resources available within the institution. This review has systematically compiled and categorized the latest developments in the diagnosis and treatment of retinoblastoma which enhanced the quality of care for this pediatric malignancy.

Retinoblastoma vulnerability to combined de novo and salvage pyrimidine ribonucleotide synthesis pharmacologic blockage

Retinoblastoma is an eye cancer that commonly affects young children. Despite significant advances, current treatments cause side effects even when administered locally, and patients may still have to undergo enucleation. This is particularly disheartening in cases of bilateral retinoblastoma. Hence, there is an urgent need for novel therapeutic strategies. Inhibitors of the enzyme dihydroorotate dehydrogenase (DHODH), which is involved in the de novo pyrimidine ribonucleotide synthesis pathway, have proven to be effective in preclinical trials against several cancers including pediatric cancers. Here we tested whether blocking pyrimidine ribonucleotide synthesis promotes retinoblastoma cell death. Cultured retinoblastoma cell lines were treated with small molecule inhibitors of DHODH alone or in combination with inhibitors of nucleoside uptake to also block the salvage pathway for pyrimidine ribonucleotide formation. On their own, DHODH inhibitors had a moderate killing effect. However, the combination with nucleoside uptake inhibitors greatly enhanced the effect of DHODH inhibition. In addition, we observed that pyrimidine ribonucleotide synthesis blockage can cause cell death in a p53 mutant retinoblastoma cell line derived from a patient with metastasis. Explaining these results, the analysis of a published patient cohort revealed that loss of chr16q22.2 (containing the DHODH gene) is amongst the most frequent alterations in retinoblastoma and that these tumors often show gains in chromosome regions expressing pyrimidine ribonucleotide salvage factors. Furthermore, these genome alterations associate with malignancy. These results indicate that targeting pyrimidine ribonucleotide synthesis may be an effective therapeutic strategy to consider as a treatment for retinoblastoma.

Serum Metabolomics of Retinoblastoma: Assessing the Differential Serum Metabolic Signatures of Unilateral and Bilateral Patients

Retinoblastoma (Rb) is the most common pediatric eye cancer. To identify the biomarkers for early diagnosis and monitoring the progression of Rb in patients, mapping of the alterations in their metabolic profiles is essential. The present study aims at exploring the metabolic disparity in serum from Rb patients and controls using NMR-based metabolomics. A total of 72 metabolites, including carbohydrates, amino acids, and organic acids, were quantified in serum samples from 24 Rb patients and 26 controls. Distinct clusters of Rb patients and controls were obtained using the partial least-squares discriminant analysis (PLS-DA) model. Further, univariate and multivariate analyses of unilateral and bilateral Rb patients with respect to their age-matched controls depicted their distinct metabolic fingerprints. Metabolites including 2-phosphoglycerate, 4-aminobutyrate, proline, O-phosphocholine, O-phosphoethanolamine, and Sn-glycero-3-phosphocholine (Sn-GPC) showed significant perturbation in both unilateral and bilateral Rb patients. However, metabolic differences among the bilateral Rb cases were more pronounced than those in unilateral Rb cases with respect to controls. In addition to major discriminatory metabolites for Rb, unilateral and bilateral Rb cases showed specific metabolic changes, which might be the result of their differential genetic/somatic mutational backgrounds. This further suggests that the aberrant metabolic perturbation in bilateral patients signifies the severity of the disease in Rb patients. The present study demonstrated that identified serum metabolites have potential to serve as a noninvasive method for detection of Rb, discriminate bilateral from unilateral Rb patients, and aid in better understanding of the RB tumor biology.

Transcriptome analysis reveals molecularly distinct subtypes in retinoblastoma

Retinoblastoma is the most frequent intraocular malignancy in children. Little is known on the molecular basis underlying the biological and clinical behavior of this cancer. Here, using gene expression profiles, we demonstrate the existence of two major retinoblastoma subtypes that can be divided into six subgroups. Subtype 1 has higher expression of cone related genes and higher percentage of RB1 germline mutation. By contrast, subtype 2 tumors harbor more genes with ganglion/neuronal features. The dedifferentiation in subtype 2 is associated with stemness features including low immune infiltration. Gene Otology analysis demonstrates that immune response regulations and visual related pathways are the key molecular difference between subtypes. Subtype 1b has the highest risk of invasiveness across all subtypes. The recognition of these molecular subtypes shed a light on the important biological and clinical perspectives for retinoblastomas.

Retinal pigment epithelium exhibits gene expression and phagocytic activity alterations when exposed to retinoblastoma chemotherapeutics

Retinoblastoma (Rb) is a rare malignant disorder affecting the developing retina of children under the age of five. Chemotherapeutic agents used for treating Rb have been associated with defects of the retinal pigment epithelium (RPE), such as hyperplasia, gliosis, and mottling. Herein, we have developed two pluripotent stem cell (PSC)-RPE models to assess the cytotoxicity of known Rb chemotherapeutics such as Melphalan, Topotecan and TW-37. Our findings demonstrate that these drugs alter the RPE by decreasing the monolayer barrier's trans-epithelial resistance and affecting the cells' phagocytic activity. Transcriptional analyses demonstrate an altered expression of genes involved in melanin and retinol processing, tight junction and apical-basal polarity pathways in both models. When applied within the clinical range, none of the drug treatments caused significant cytotoxic effects, changes to the apical-basal polarity, tight junction network or cell cycle. Together, our results demonstrate that although the most commonly used Rb chemotherapeutic drugs do not cause cytotoxicity in RPE, their application in vitro leads to compromised phagocytosis and strength of the barrier function, in addition to changes in gene expression that could alter the visual cycle in vivo. Our data demonstrate that widely used Rb chemotherapeutic drugs can have a deleterious impact on RPE cells and thus great care has to be exercised with regard to their delivery so the adjacent healthy RPE is not damaged during the course of tumor eradication.

Identification of novel genes by targeted exome sequencing in Retinoblastoma

BACKGROUND

Retinoblastoma (RB) is initiated by mutation in both alleles of RB1 gene. However, few cases may occur even in the absence of RB1 mutation suggesting the role of genes other than RB1.

METHODOLOGY

The current study was planned to utilize targeted exome sequencing in Indian RB patients affected with unilateral non-familial RB. 75 unilateral RB patients below 5 years of age were enrolled. Genomic DNA was extracted from blood and tumor tissue. From peripheral blood DNA, all coding and exon/intron regions were amplified using PCR and direct sequencing. Cases which did not harbor pathogenic variants in peripheral blood DNA were further screened for mutations in their tumor tissue DNA using targeted exome sequencing. Three pathogenicity prediction tools (Mutation Taster, SIFT, and PolyPhen-2) were used to determine the pathogenicity of non-synonymous variations. An in-house bioinformatics pipeline was devised for the mutation screening by targeted exome sequencing. Protein modeling studies were also done to predict the effect of the mutations on the protein structure and function.

RESULTS

Using the mentioned approach, we found two novel variants (g.69673_69674insT and g.48373314C>A) in RB1 gene in peripheral blood DNA. We also found novel variants in eight genes (RB1, ACAD11, GPR151, KCNA1, OTOR, SOX30, ARL11, and MYCT1) that may be associated with RB pathogenesis.

CONCLUSION

The present study expands our current knowledge regarding the genomic landscape of RB and also highlights the importance of NGS technologies to detect genes and novel variants that may play an important role in cancer initiation, progression, and prognosis.

Development and characterization of DNA aptamer against Retinoblastoma by Cell-SELEX

Retinoblastoma (RB) is the most common paediatric intraocular tumour. The management of RB has improved the survival and vision with recent advances in the treatment. Improved therapeutic approaches focussing on targeting tumours and minimizing the treatment-associated side effects are being developed. In this study, we generated a ssDNA aptamer against RB by cell-SELEX and high-throughput sequencing using Weri-RB1 cell line as the target, and Muller glial cell line Mio-M1 as the control. Three aptamers were selected based on the number of repetitions in NGS and phylogenetic relationship and evaluated by flow cytometry to assess their binding affinity and selectivity. The dissociation constant, Kd values of three selected aptamers were found to be in the nanomolar range. Aptamer VRF-CSRB-01 with the best binding affinity and a Kd value of 49.41 ± 7.87 nM was further characterized. The proteinase and temperature treatment indicated that VRF-CSRB-01 targets surface proteins, and has a good binding affinity and excellent selectivity under physiological conditions. The aptamer VRF-CSRB-01 was stable over 72 h in serum and 96 h in cerebral spinal fluid and vitreous. With the high affinity, specificity, stability and specific recognition of clinical RB tumours, VRF-CSRB-01 aptamer holds potential for application in diagnosis and targeting RB.

Single-cell characterization of malignant phenotypes and microenvironment alteration in retinoblastoma

Retinoblastoma (RB) is the most common primary intraocular malignancy of childhood. It is known that the tumor microenvironment (TME) regulates tumorigenesis and metastasis. However, how the malignant progression in RB is determined by the heterogeneity of tumor cells and TME remains uncharacterized. Here, we conducted integrative single-cell transcriptome and whole-exome sequencing analysis of RB patients with detailed pathological and clinical measurements. By single-cell transcriptomic sequencing, we profiled around 70,000 cells from tumor samples of seven RB patients. We identified that the major cell types in RB were cone precursor-like (CP-like) and MKI67+ cone precursor (MKI67+ CP) cells. By integrating copy number variation (CNV) analysis, we found that RB samples had large clonal heterogeneity, where the malignant MKI67+ CP cells had significantly larger copy number changes. Enrichment analysis revealed that the conversion of CP-like to MKI67+ CP resulted in the loss of photoreceptor function and increased cell proliferation ability. The TME in RB was composed of tumor-associated macrophages (TAMs), astrocyte-like, and cancer-associated fibroblasts (CAFs). Particularly, during the invasion process, TAMs created an immunosuppressive environment, in which the proportion of TAMs decreased, M1-type macrophage was lost, and the TAMs-related immune functions were depressed. Finally, we identified that TAMs regulated tumor cells through GRN and MIF signaling pathways, while TAMs self-regulated through inhibition of CCL and GALECTIN signaling pathways during the invasion process. Altogether, our study creates a detailed transcriptomic map of RB with single-cell characterization of malignant phenotypes and provides novel molecular insights into the occurrence and progression of RB.

The Oncogene MYCN Modulates Glycolytic and Invasive Genes to Enhance Cell Viability and Migration in Human Retinoblastoma

Retinoblastoma is usually initiated by biallelic RB1 gene inactivation. In addition, MYCN copy number alterations also contribute to RB pathogenesis. However, MYCN expression, its role in disease progression and correlation with RB histological risk factors are not well understood. We studied the expression of MYCN in enucleated RB patient specimens by immunohistochemistry. MYCN is overexpressed in RB compared to control retina. Our microarray gene expression analysis followed by qRT-PCR validation revealed that genes involved in glucose metabolism and migration are significantly downregulated in MYCN knockdown cells. Further, targeting MYCN in RB cells using small molecule compounds or shRNAs led to decreased cell survival and migration, increased apoptosis and cell cycle arrest, suggesting that MYCN inhibition can be a potential therapeutic strategy. We also noted that MYCN inhibition results in reduction in glucose uptake, lactate production, ROS levels and gelatinolytic activity of active-MMP9, explaining a possible mechanism of MYCN in RB. Taking clues from our findings, we tested a combination treatment of RB cells with carboplatin and MYCN inhibitors to find enhanced therapeutic efficacy compared to single drug treatment. Thus, MYCN inhibition can be a potential therapeutic strategy in combination with existing chemotherapy drugs to restrict tumor cell growth in RB.

Giant Y79 retinoblastoma cells contain functionally active T-type calcium channels

Retinoblastoma is the most common malignant intraocular tumor in children. Y79 human retinoblastoma cells are in vitro models of retinal tumors used for drug screening. Undifferentiated Y79 cells originate from a primitive multi-potential neuroectodermal cell and express neuronal and glial properties. However, the nature of cellular heterogeneity in Y79 cells is unclear because functional methods to characterize neurons or glial cells have not been employed to Y79 cells. Here, we perform patch-clamp recordings to characterize electrophysiological properties in retinoblastoma cells. We identified a population of large-sized Y79 cells (i.e., giant cells, ~ 40-µm diameter), hyperpolarized resting membrane potential (-54 mV), and low input resistance (~ 600 MΩ), indicating electrically mature cells. We also found that giant Y79 cells contain increased density of T-type calcium channels. Finally, we found that T-type calcium channels are active only in giant cells suggesting that cancer treatments aimed to prevent calcium influx in retinoblastomas should be tested in giant cells.

Pre-clinical evidence that salinomycin is active against retinoblastoma via inducing mitochondrial dysfunction, oxidative damage and AMPK activation

The poor outcomes in retinoblastoma necessitate new treatments. Salinomycin is an attractive candidate, and has demonstrated selective anti-cancer properties in different cancer types. This work addressed the efficacy of salinomycin in retinoblastoma models and probe the associated mechanisms. Cellular functional assays were conducted to determine the effects salinomycin in vitro. Xenograft retinoblastoma mouse model was established to investigate the efficacy of salinomycin in vivo. Biochemical assays were conducted to analyze the mechanism of salinomycin's action focusing on mitochondrial functions, energy reduction-related signaling pathways. Salinomycin has positive effects towards retinoblastoma cells regardless of heterogeneity through suppressing growth and inducing apoptosis. Salinomycin also specifically inhibits cells displaying stemness and highly invasive phenotypes. Using retinoblastoma xenograft mouse model, we show that salinomycin at non-toxic dose effectively inhibits growth and induces apoptosis. Mechanistic studies show that salinomycin inhibits mitochondrial respiration via specifically suppressing complex I and II activities, reduces mitochondrial membrane potential and decreases energy reduction, followed by induction of oxidative stress and damage, AMPK activation and mTOR inhibition. Our study highlights that adding salinomycin to the existing treatment armamentarium for retinoblastoma is beneficial.

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?

Purpose

Retinoblastoma (Rb) is a malignant neoplasm arising during retinal development from mutations in the RB1 gene. Loss or inactivation of both copies of RB1 results in initiation of retinoblastoma tumors; however, additional genetic changes are needed for the continued growth and spread of the tumor. Ex vivo research has shown that in humans, retinoblastoma may initiate from RB1-depleted cone precursors. Notwithstanding, it has not been possible to assess the full spectrum of clonal types within the tumor itself in vivo and the molecular changes occurring at the cells of origin, enabling their malignant conversion. To overcome these challenges, we have performed the first single cell (sc) RNA- and ATAC-Seq analyses of primary tumor tissues, enabling us to dissect the transcriptional and chromatin accessibility heterogeneity of proliferating cone precursors in human Rb tumors.

Methods

Two Rb tumors each characterized by two pathogenic RB1 mutations were dissociated to single cells and subjected to scRNA-Seq and scATAC-Seq using the 10× Genomics platform. In addition, nine human embryonic and fetal retina samples were dissociated to single cells and subjected to scRNA- and ATAC-Seq analyses. The scRNA- and ATAC-Seq data were embedded using Uniform Manifold Approximation and Projection and clustered with Seurat graph-based clustering. Integrated scATAC-Seq analysis of Rb tumors and human embryonic/fetal retina samples was performed to identify Rb cone enriched subclusters. Pseudo time analysis of proliferating cones in the Rb samples was performed with Monocle. Ingenuity Pathway Analysis was used to identify the signaling pathway and upstream regulators in the Rb cone-enriched subclusters.

Results

Our single cell analyses revealed the predominant presence of cone precursors at different stages of the cell cycle in the Rb tumors and among those identified the G2/M subset as the cell type of origin. scATAC-Seq analysis identified two Rb enriched cone subclusters, each characterized by activation of different upstream regulators and signaling pathways, enabling proliferating cone precursors to escape cell cycle arrest and/or apoptosis.

Conclusions

Our study provides evidence of Rb tumor heterogeneity and defines molecular pathways that can be targeted to define new treatment strategies.

Chromosome 6p Amplification in Aqueous Humor Cell-Free DNA Is a Prognostic Biomarker for Retinoblastoma Ocular Survival

Aqueous humor contains tumor-derived cell-free DNA (cfDNA) and can serve as a liquid biopsy for retinoblastoma. We previously associated somatic copy-number alteration (SCNA) 6p gain with a 10-fold increased risk of enucleation. Here we provide a 2-year update to further explore 6p gain as a prognostic biomarker for ocular survival. Patients diagnosed with retinoblastoma from December 2014 to July 2019 from whom aqueous humor was sampled were included. cfDNA was extracted and shallow whole-genome sequencing performed to identify highly recurrent retinoblastoma SCNAs (gain of 1q, 2p, 6p, loss of 13q, 16q). 116 aqueous humor samples from 50 eyes of 46 patients were included: 27 eyes were salvaged, 23 were enucleated. Highly recurrent retinoblastoma SCNAs were found in 66% eyes. 6p gain was the most prevalent SCNA (50% eyes). It was particularly more prevalent in enucleated eyes (73.9%) than in salvaged eyes (29.6%; P = 0.004). 6p gain in aqueous humor cfDNA portended nearly 10-fold increased odds of enucleation (OR = 9.87; 95% confidence interval = 1.75-55.65; P = 0.009). In the enucleated eyes, 6p gain was associated with aggressive histopathologic features, including necrosis, higher degrees of anaplasia, and focal invasion of ocular structures. With extended follow-up and nearly double the aqueous humor samples, we continue to demonstrate 6p gain as a potential prognostic biomarker for retinoblastoma. IMPLICATIONS: Aqueous humor is a high-yield source of tumor-derived DNA in retinoblastoma eyes. Detection of 6p gain in the aqueous humor allows for targeted, patient-centered therapies based on this molecular prognostic marker. Prospective, multicenter studies with aqueous humor sampled from all eyes at diagnosis are warranted to validate these findings.

Gene expression profiling of tumor stroma interactions in retinoblastoma

We aimed to identify the critical molecular pathways altered upon tumor stroma interactions in retinoblastoma (RB). In vitro 2 D cocultures of RB tumor cells (Weri-Rb-1 and NCC-RbC-51) with primary bone marrow stromal cells (BMSC) was established. Global gene expression patterns in coculture samples were assessed using Affymetrix Prime view human gene chip microarray and followed with bioinformatics analyses. Key upregulated genes from Weri-Rb-1 + BMSC and NCC-RbC-51 + BMSC coculture were validated using qRT-PCR to ascertain their role in RB progression. Whole genome microarray experiments identified significant (P ≤ 0.05, 1.1 log 2 FC) transcriptome level changes induced upon coculture of RB cells with BMSC. A total of 1155 genes were downregulated and 1083 upregulated in Weri-Rb-1 + BMSC coculture. Similarly, 1865 genes showed downregulation and 1644 genes were upregulation in NCC-RbC-51 + BMSC coculture. The upregulated genes were significantly associated with pathways of focal adhesion, PI3K-Akt signalling, ECM-receptor interaction, JAK-STAT, TGF-β signalling thus contributing to RB progression. Validation of key genes by qRT-PCR revealed significant overexpression of IL8, IL6, MYC and SMAD3 in the case of Weri-Rb-1 + BMSC coculture and IL6 in the case of NCC-RbC-51 + BMSC coculture. The microarray expression study on in vitro RB coculture models revealed the pathways that could be involved in the progression of RB. The gene signature obtained in a stimulated model when a growing tumor interacts with its microenvironment may provide new horizons for potential targeted therapy in RB.

Crocin inhibits proliferation and induces apoptosis through suppressing MYCN expression in retinoblastoma

The pathogenetic mechanisms of retinoblastoma are still not yet fully elucidated, putting limits to efficacious treatment. Crocin is the main component of saffron, which exhibits significant antitumorigenic properties. The aim of this paper is to investigate the effect of crocin on retinoblastoma. The effects of crocin on the proliferation of human retinoblastoma cells were determined by the 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide assay, cell number assay, and colony formation assay. Cell apoptosis induced by crocin was measured by flow cytometry analysis. Cleaved poly(ADP-ribose) polymerase and cleaved caspase-3 were tested by western blot analysis. The expression levels of MYCN were assessed by western blot and quantitative polymerase chain reaction and the stability of MYCN messenger RNA was determined by in vitro RNA degradation assays. We found that crocin significantly inhibited the cell proliferation and clonogenicity and induced cell apoptosis in Y79 and WERI-RB-1 cells. In addition, crocin treatment significantly reduced the expression and the stability of MYCN. Besides, overexpression of MYCN rescued the inhibitory effect of crocin in Y79 cells. Our findings suggest that crocin exhibits antitumorigenic effects in human retinoblastoma cell lines through a MYCN-dependent manner, which may provide guidance to logical therapeutic designs in prevention and treatment of retinoblastoma.