Abstract
Retinoblastomas form in response to biallelic RB1 mutations or MYCN amplification and progress to more aggressive and therapy-resistant phenotypes through accumulation of secondary genomic changes. Progression-related changes include recurrent somatic copy number alterations and typically non-recurrent nucleotide variants, including synonymous and non-coding variants, whose significance has been unclear. To determine if nucleotide variants recurrently affect specific biological processes, we identified altered genes and over-represented variant gene ontologies in 168 exome or whole-genome-sequenced retinoblastomas and 12 tumor-matched cell lines. In addition to RB1 mutations, MYCN amplification, and established retinoblastoma somatic copy number alterations, the analyses revealed enrichment of variant genes related to diverse biological processes including histone monoubiquitination, mRNA processing (P) body assembly, and mitotic sister chromatid segregation and cytokinesis. Importantly, non-coding and synonymous variants increased the enrichment significance of each over-represented biological process term. To assess the effects of such mutations, we examined the consequences of a 3' UTR variant of PCGF3 (a BCOR-binding component of Polycomb repressive complex I), dual 3' UTR variants of CDC14B (a regulator of sister chromatid segregation), and a synonymous variant of DYNC1H1 (a regulator of P-body assembly). One PCGF3 and one of two CDC14B 3' UTR variants impaired gene expression whereas a base-edited DYNC1H1 synonymous variant altered protease sensitivity and stability. Retinoblastoma cell lines retained only ~50% of variants detected in tumors and enriched for new variants affecting p53 signaling. These findings reveal potentially important differences in retinoblastoma cell lines and tumors and implicate synonymous and non-coding variants, along with non-synonymous variants, in retinoblastoma oncogenesis.