Uphyloplot2: visualizing phylogenetic trees from single-cell RNA-seq data.

BACKGROUND: Recent advances in single cell sequencing technologies allow for greater resolution in assessing tumor clonality using chromosome copy number variations (CNVs). While single cell DNA sequencing technologies are ideal to identify tumor sub-clones, they remain expensive and in contrast to single cell RNA-seq (scRNA-seq) methods are more limited in the data they generate. However, CNV data can be inferred from scRNA-seq and bulk RNA-seq, for which several tools have been developed, including inferCNV, CaSpER, and HoneyBADGER. Inferences regarding tumor clonality from CNV data (and other sources) are frequently visualized using phylogenetic plots, which previously required time-consuming and error-prone, manual analysis. RESULTS: Here, we present Uphyloplot2, a python script that generates phylogenetic plots directly from inferred RNA-seq data, or any Newick formatted dendrogram file. The tool is publicly available at https://github.com/harbourlab/UPhyloplot2/ . CONCLUSIONS: Uphyloplot2 is an easy-to-use tool to generate phylogenetic plots to depict tumor clonality from scRNA-seq data and other sources.

Gain of function of ASXL1 truncating protein in the pathogenesis of myeloid malignancies.

Additional Sex Combs-Like 1 (ASXL1) is mutated at a high frequency in all forms of myeloid malignancies associated with poor prognosis. We generated a Vav1 promoter-driven Flag-Asxl1Y588X transgenic mouse model, Asxl1Y588X Tg, to express a truncated FLAG-ASXL1aa1-587 protein in the hematopoietic system. The Asxl1Y588X Tg mice had an enlarged hematopoietic stem cell (HSC) pool, shortened survival, and predisposition to a spectrum of myeloid malignancies, thereby recapitulating the characteristics of myeloid malignancy patients with ASXL1 mutations. ATAC- and RNA-sequencing analyses revealed that the ASXL1aa1-587 truncating protein expression results in more open chromatin in cKit+ cells compared with wild-type cells, accompanied by dysregulated expression of genes critical for HSC self-renewal and differentiation. Liquid chromatography-tandem mass spectrometry and coimmunoprecipitation experiments showed that ASXL1aa1-587 acquired an interaction with BRD4. An epigenetic drug screening demonstrated a hypersensitivity of Asxl1Y588X Tg bone marrow cells to BET bromodomain inhibitors. This study demonstrates that ASXL1aa1-587 plays a gain-of-function role in promoting myeloid malignancies. Our model provides a powerful platform to test therapeutic approaches of targeting the ASXL1 truncation mutations in myeloid malignancies.

PRAME induces genomic instability in uveal melanoma.

PRAME is a CUL2 ubiquitin ligase subunit that is normally expressed in the testis but becomes aberrantly overexpressed in many cancer types in association with aneuploidy and metastasis. Here, we show that PRAME is expressed predominantly in spermatogonia around the time of meiotic crossing-over in coordination with genes mediating DNA double strand break repair. Expression of PRAME in somatic cells upregulates pathways involved in meiosis, chromosome segregation and DNA repair, and it leads to increased DNA double strand breaks, telomere dysfunction and aneuploidy in neoplastic and non-neoplastic cells. This effect is mediated at least in part by ubiquitination of SMC1A and altered cohesin function. PRAME expression renders cells susceptible to inhibition of PARP1/2, suggesting increased dependence on alternative base excision repair pathways. These findings reveal a distinct oncogenic function of PRAME that can be targeted therapeutically in cancer.

PRAME induces genomic instability in uveal melanoma.

PRAME is a CUL2 ubiquitin ligase subunit that is normally expressed in the testis but becomes aberrantly overexpressed in many cancer types in association with aneuploidy and metastasis. Here, we show that PRAME is expressed predominantly in spermatogonia around the time of meiotic crossing-over in coordination with genes mediating DNA double strand break repair. Expression of PRAME in somatic cells upregulates pathways involved in meiosis, chromosome segregation and DNA repair, and it leads to increased DNA double strand breaks, telomere dysfunction and aneuploidy in neoplastic and non-neoplastic cells. This effect is mediated at least in part by ubiquitination of SMC1A and altered cohesin function. PRAME expression renders cells susceptible to inhibition of PARP1/2, suggesting increased dependence on alternative base excision repair pathways. These findings reveal a distinct oncogenic function of PRAME than can be targeted therapeutically in cancer.

HDAC11 activity contributes to MEK inhibitor escape in uveal melanoma.

We previously demonstrated that pan-HDAC inhibitors could limit escape from MEK inhibitor (MEKi) therapy in uveal melanoma (UM) through suppression of AKT and YAP/TAZ signaling. Here, we focused on the role of specific HDACs in therapy adaptation. Class 2 UM displayed higher expression of HDACs 1, 2, and 3 than Class 1, whereas HDACs 6, 8, and 11 were uniformly expressed. Treatment of UM cells with MEKi led to modulation of multiple HDACs, with the strongest increases observed in HDAC11. RNA-seq analysis showed MEKi to decrease the expression of multiple HDAC11 target genes. Silencing of HDAC11 significantly reduced protein deacetylation, enhanced the apoptotic response to MEKi and reduced growth in long-term colony formation assays across multiple UM cell lines. Knockdown of HDAC11 led to decreased expression of TAZ in some UM cell lines, accompanied by decreased YAP/TAZ transcriptional activity and reduced expression of multiple YAP/TAZ target genes. Further studies showed this decrease in TAZ expression to be associated with increased LKB1 activation and modulation of glycolysis. In an in vivo model of uveal melanoma, silencing of HDAC11 limited the escape to MEKi therapy, an effect associated with reduced levels of Ki67 staining and increased cleaved caspase-3. We have demonstrated a novel role for adaptive HDAC11 activity in UM cells, that in some cases modulates YAP/TAZ signaling leading to MEKi escape.

RB1 loss triggers dependence on ESRRG in retinoblastoma.

Retinoblastoma (Rb) is a deadly childhood eye cancer that is classically initiated by inactivation of the RB1 tumor suppressor. Clinical management continues to rely on nonspecific chemotherapeutic agents that are associated with treatment resistance and toxicity. Here, we analyzed 103 whole exomes, 20 whole transcriptomes, 5 single-cell transcriptomes, and 4 whole genomes from primary Rb tumors to identify previously unknown Rb dependencies. Several recurrent genomic aberrations implicate estrogen-related receptor gamma (ESRRG) in Rb pathogenesis. RB1 directly interacts with and inhibits ESRRG, and RB1 loss uncouples ESRRG from negative regulation. ESRRG regulates genes involved in retinogenesis and oxygen metabolism in Rb cells. ESRRG is preferentially expressed in hypoxic Rb cells in vivo. Depletion or inhibition of ESRRG causes marked Rb cell death, which is exacerbated in hypoxia. These findings reveal a previously unidentified dependency of Rb cells on ESRRG, and they implicate ESRRG as a potential therapeutic vulnerability in Rb.

PieParty: visualizing cells from scRNA-seq data as pie charts.

Single-cell RNA sequencing (scRNA-seq) has been a transformative technology in many research fields. Dimensional reduction techniques such as UMAP and tSNE are used to visualize scRNA-seq data in two or three dimensions for cells to be clustered in biologically meaningful ways. Subsequently, gene expression is frequently mapped onto these plots to show the distribution of gene expression across the plots, for instance to distinguish cell types. However, plotting each cell with only a single color leads to repetitive and unintuitive representations. Here, we present PieParty, which allows scRNA-seq data to be plotted such that every cell is represented as a pie chart, and every slice in the pie charts corresponds to the gene expression of a single gene. This allows for the simultaneous visualization of the expression of multiple genes and gene networks. The resulting figures are information dense, space efficient, and highly intuitive. PieParty is publicly available on GitHub at https://github.com/harbourlab/PieParty.

Multiregional genetic evolution of metastatic uveal melanoma.

Uveal melanoma (UM) is the most common primary intraocular malignancy in adults and leads to deadly metastases for which there is no approved treatment. Genetic events driving early tumor development are well-described, but those occurring later during metastatic progression remain poorly understood. We performed multiregional genomic sequencing on 22 tumors collected from two patients with widely metastatic UM who underwent rapid autopsy. We observed multiple seeding events from the primary tumors, metastasis-to-metastasis seeding, polyclonal seeding, and late driver variants in ATM, KRAS, and other genes previously unreported in UM. These findings reveal previously unrecognized temporal and anatomic complexity in the genetic evolution of metastatic uveal melanoma, and they highlight the distinction between early and late phases of UM genetic evolution with implications for novel therapeutic approaches.

Analytical Validation and Performance of a 7-Gene Next-Generation Sequencing Panel in Uveal Melanoma.

INTRODUCTION: Gene expression profiling (GEP) is widely used for prognostication in patients with uveal melanoma (UM). Because biopsy tissue is limited, it is critical to obtain as much genomic information as possible from each sample. Combined application of both GEP and next-generation sequencing (NGS) allows for analysis of RNA and DNA from a single biopsy sample, offers additional prognostic information, and can potentially inform therapy selection. This study evaluated the analytical performance of a targeted custom NGS panel for mutational profiling of 7 genes commonly mutated in UM. METHODS: One hundred five primary UM tumors were analyzed, including 37 formalin-fixed paraffin-embedded (FFPE) and 68 fine-needle aspiration biopsy specimens. Sequencing was performed on the Ion GeneStudio S5 platform to an average read depth of >500X per region of interest. RESULTS: The 7-gene panel achieved a positive percent agreement of 100% for detection of both single-nucleotide variants and insertions/deletions, with a technical positive predictive value of 98.8% and 100%, respectively. Intra-assay and inter-assay concordance studies confirmed the assay's reproducibility and repeatability. DISCUSSION/CONCLUSION: The 7-gene panel is a robust, highly accurate NGS test that can be successfully performed, along with GEP, from a single small-gauge needle biopsy sample or FFPE specimen.

On the in vivo origin of human nasal mesenchymal stem cell cultures.

OBJECTIVES: Mesenchymal stem cells (MSCs), classically expanded in culture from bone marrow, are of broad interest to the regenerative medicine community. Human nasal turbinate mesenchymal-like stem cell cultures have also been described, defined by an in vitro phenotype similar to bone marrow MSCs. Nonetheless, the identity in vivo of the cells that give rise to nasal MSC-like cultures remains unclear, and these cells are often suggested to be related to olfactory lineages. Here, we sought to define the in vivo phenotype of human nasal MSC-like cells. METHODS: Human turbinate tissue samples were used for RNA and immunohistochemical analysis. We also analyzed a recently published single cell RNA-sequencing dataset from adult human olfactory and respiratory mucosa samples from our lab, to focus on cell populations expressing MSC markers. Immunochemistry was performed to stain turbinate sections and nasal MSC cultures for selected markers. RESULTS: While there is no single MSC-specific gene, we identified a human nasal mucosal cell population in vivo that uniquely expressed transcripts characteristic of typical MSC cultures, including ENG (CD105), NES, and CD34, and lacked expression of other transcripts associated with surface epithelia. The expression of transcription factors such as SOX17, EBF1, and FOXP1 suggests cells in the MSC-like cluster maintain an ability to direct cell fate, consistent with the behavior of nasal MSC-like cells in vitro. SOX17 was found to be uniformly expressed by nasal MSC cultures, consistent with the in vivo data. Immunohistochemistry of human nasal tissue samples indicated that ENG, CD34, and SOX17 expression localized selectively to cells surrounding blood vessels in the lamina propria. CONCLUSION: Our findings provide evidence that the in vivo origin of nasal MSC-like cultures is likely a vascular or pericyte population, rather than cells related to the olfactory neuronal lineage. LEVEL OF EVIDENCE: NA.

A novel cardiomyogenic role for Isl1+ neural crest cells in the inflow tract.

The degree to which populations of cardiac progenitors (CPCs) persist in the postnatal heart remains a controversial issue in cardiobiology. To address this question, we conducted a spatiotemporally resolved analysis of CPC deployment dynamics, tracking cells expressing the pan-CPC gene Isl1 Most CPCs undergo programmed silencing during early cardiogenesis through proteasome-mediated and PRC2 (Polycomb group repressive complex 2)-mediated Isl1 repression, selectively in the outflow tract. A notable exception is a domain of cardiac neural crest cells (CNCs) in the inflow tract. These "dorsal CNCs" are regulated through a Wnt/ß-catenin/Isl1 feedback loop and generate a limited number of trabecular cardiomyocytes that undergo multiple clonal divisions during compaction, to eventually produce ~10% of the biventricular myocardium. After birth, CNCs continue to generate cardiomyocytes that, however, exhibit diminished clonal amplification dynamics. Thus, although the postnatal heart sustains cardiomyocyte-producing CNCs, their regenerative potential is likely diminished by the loss of trabeculation-like proliferative properties.

Single-cell analysis of olfactory neurogenesis and differentiation in adult humans.

The presence of active neurogenic niches in adult humans is controversial. We focused attention to the human olfactory neuroepithelium, an extracranial site supplying input to the olfactory bulbs of the brain. Using single-cell RNA sequencing analyzing 28,726 cells, we identified neural stem cell and neural progenitor cell pools and neurons. Additionally, we detailed the expression of 140 olfactory receptors. These data from the olfactory neuroepithelium niche provide evidence that neuron production may continue for decades in humans.

Single-cell analysis reveals new evolutionary complexity in uveal melanoma.

Uveal melanoma (UM) is a highly metastatic cancer that, in contrast to cutaneous melanoma, is largely unresponsive to checkpoint immunotherapy. Here, we interrogate the tumor microenvironment at single-cell resolution using scRNA-seq of 59,915 tumor and non-neoplastic cells from 8 primary and 3 metastatic samples. Tumor cells reveal novel subclonal genomic complexity and transcriptional states. Tumor-infiltrating immune cells comprise a previously unrecognized diversity of cell types, including CD8+ T cells predominantly expressing the checkpoint marker LAG3, rather than PD1 or CTLA4. V(D)J analysis shows clonally expanded T cells, indicating that they are capable of mounting an immune response. An indolent liver metastasis from a class 1B UM is infiltrated with clonally expanded plasma cells, indicative of antibody-mediated immunity. This complex ecosystem of tumor and immune cells provides new insights into UM biology, and LAG3 is identified as a potential candidate for immune checkpoint blockade in patients with high risk UM.

Intraocular Metastasis in Unilateral Multifocal Uveal Melanoma Without Melanocytosis or Germline BAP1 Mutations.

Importance: There has been speculation on the pathogenesis of unilateral multifocal uveal melanoma, but there remains no convincing explanation. Genetic analysis suggests that unilateral multifocal uveal melanoma may represent intraocular metastasis with increased risk of systemic metastasis. Objective: To evaluate the pathogenesis of unilateral multifocal uveal melanoma. Design, Setting, and Participants: This clinical case series was conducted in tertiary academic ocular oncology referral centers and included patients with unilateral multifocal uveal melanoma. Main Outcomes and Measures: Gene expression and mutation profiling of tumor samples. Results: Four patients (all male; age range, 54-77 years) who were diagnosed with uveal melanoma were treated with plaque brachytherapy, and subsequently developed a second discrete uveal melanoma in the same eye were included. None demonstrated ocular or oculodermal melanocytosis. All 8 tumors available for analysis exhibited class 2 gene expression profiles. In all 4 cases, the initial and subsequent tumors were available for targeted DNA sequencing and identical driver mutations were present in both tumors. Data were collected from September 2015 to August 2018. Conclusions and Relevance: Unilateral multifocal uveal melanoma in the absence of ocular melanocytosis appears to occur preferentially in tumors with the class 2 gene expression profile and a BRCA1-associated protein 1 gene (BAP1) mutation. The presence of identical BAP1 mutations in multiple tumors in the same eye in the absence of a germline BAP1 mutation suggests intraocular metastasis rather than independent primary tumors. These findings indicate that the first site of metastasis can be within the eye itself and suggest that patients with unilateral multifocal uveal melanoma may be at increased risk of systemic metastasis.

BAP1 regulates epigenetic switch from pluripotency to differentiation in developmental lineages giving rise to BAP1-mutant cancers.

The BAP1 tumor suppressor is mutated in many human cancers such as uveal melanoma, leading to poor patient outcome. It remains unclear how BAP1 functions in normal biology or how its loss promotes cancer progression. Here, we show that Bap1 is critical for commitment to ectoderm, mesoderm, and neural crest lineages during Xenopus laevis development. Bap1 loss causes transcriptional silencing and failure of H3K27ac to accumulate at promoters of key genes regulating pluripotency-to-commitment transition, similar to findings in uveal melanoma. The Bap1-deficient phenotype can be rescued with human BAP1, by pharmacologic inhibition of histone deacetylase (HDAC) activity or by specific knockdown of Hdac4. Similarly, BAP1-deficient uveal melanoma cells are preferentially vulnerable to HDAC4 depletion. These findings show that Bap1 regulates lineage commitment through H3K27ac-mediated transcriptional activation, at least in part, by modulation of Hdac4, and they provide insights into how BAP1 loss promotes cancer progression.

Genomic evolution of uveal melanoma arising in ocular melanocytosis.

Ocular melanocytosis is the most important predisposing condition for the eye cancer uveal melanoma (UM). Here, we present a patient who developed UM arising within ocular melanocytosis who was treated with enucleation (eye removal), which provided an invaluable opportunity to interrogate both the UM and adjacent uveal tissue containing the melanocytosis using whole-exome and deep-targeted sequencing. This analysis revealed a clonal PLCB4 mutation in the melanocytosis, confirming that this is indeed a neoplastic condition and explaining why it predisposes to UM. This mutation was present in 100% of analyzed UM cells, indicating that a PLCB4-mutant cell gave rise to the UM. The earliest aberrations specific to the tumor were loss of Chromosomes 1p, 3, and 9p, which were present in virtually all tumor cells. A mutation in BAP1 arose later on the other copy of Chromosome 3 in a tumor subclone, followed by a gain of Chromosome 8q. These findings provide a mechanistic explanation for the well-known clinical association between ocular melanocytosis and UM by showing that this predisposing condition introduces the first "hit" and thereby increases the stochastic likelihood of acquiring further aberrations leading to UM.

HDAC Inhibition Enhances the In Vivo Efficacy of MEK Inhibitor Therapy in Uveal Melanoma.

PURPOSE: The clinical use of MEK inhibitors in uveal melanoma is limited by the rapid acquisition of resistance. This study has used multiomics approaches and drug screens to identify the pan-HDAC inhibitor panobinostat as an effective strategy to limit MEK inhibitor resistance.Experimental Design: Mass spectrometry-based proteomics and RNA-Seq were used to identify the signaling pathways involved in the escape of uveal melanoma cells from MEK inhibitor therapy. Mechanistic studies were performed to evaluate the escape pathways identified, and the efficacy of the MEK-HDAC inhibitor combination was demonstrated in multiple in vivo models of uveal melanoma. RESULTS: We identified a number of putative escape pathways that were upregulated following MEK inhibition, including the PI3K/AKT pathway, ROR1/2, and IGF-1R signaling. MEK inhibition was also associated with increased GPCR expression, particularly the endothelin B receptor, and this contributed to therapeutic escape through ET-3-mediated YAP signaling. A screen of 289 clinical grade compounds identified HDAC inhibitors as potential candidates that suppressed the adaptive YAP and AKT signaling that followed MEK inhibition. In vivo, the MEK-HDAC inhibitor combination outperformed either agent alone, leading to a long-term decrease of tumor growth in both subcutaneous and liver metastasis models and the suppression of adaptive PI3K/AKT and YAP signaling. CONCLUSIONS: Together, our studies have identified GPCR-mediated YAP activation and RTK-driven AKT signaling as key pathways involved in the escape of uveal melanoma cells from MEK inhibition. We further demonstrate that HDAC inhibition is a promising combination partner for MEK inhibitors in advanced uveal melanoma.

Punctuated evolution of canonical genomic aberrations in uveal melanoma.

Cancer is thought to arise through the accumulation of genomic aberrations evolving under Darwinian selection. However, it remains unclear when the aberrations associated with metastasis emerge during tumor evolution. Uveal melanoma (UM) is the most common primary eye cancer and frequently leads to metastatic death, which is strongly linked to BAP1 mutations. Accordingly, UM is ideally suited for studying the clonal evolution of metastatic competence. Here we analyze sequencing data from 151 primary UM samples using a customized bioinformatic pipeline, to improve detection of BAP1 mutations and infer the clonal relationships among genomic aberrations. Strikingly, we find BAP1 mutations and other canonical genomic aberrations usually arise in an early punctuated burst, followed by neutral evolution extending to the time of clinical detection. This implies that the metastatic proclivity of UM is "set in stone" early in tumor evolution and may explain why advances in primary treatment have not improved survival.

Drug and disease signature integration identifies synergistic combinations in glioblastoma.

Glioblastoma (GBM) is the most common primary adult brain tumor. Despite extensive efforts, the median survival for GBM patients is approximately 14 months. GBM therapy could benefit greatly from patient-specific targeted therapies that maximize treatment efficacy. Here we report a platform termed SynergySeq to identify drug combinations for the treatment of GBM by integrating information from The Cancer Genome Atlas (TCGA) and the Library of Integrated Network-Based Cellular Signatures (LINCS). We identify differentially expressed genes in GBM samples and devise a consensus gene expression signature for each compound using LINCS L1000 transcriptional profiling data. The SynergySeq platform computes disease discordance and drug concordance to identify combinations of FDA-approved drugs that induce a synergistic response in GBM. Collectively, our studies demonstrate that combining disease-specific gene expression signatures with LINCS small molecule perturbagen-response signatures can identify preclinical combinations for GBM, which can potentially be tested in humans.

Epigenetic reprogramming and aberrant expression of PRAME are associated with increased metastatic risk in Class 1 and Class 2 uveal melanomas.

BACKGROUND: We previously identified PRAME as a biomarker for metastatic risk in Class 1 uveal melanomas. In this study, we sought to define a threshold value for positive PRAME expression (PRAME+) in a large dataset, identify factors associated with PRAME expression, evaluate the prognostic value of PRAME in Class 2 uveal melanomas, and determine whether PRAME expression is associated with aberrant hypomethylation of the PRAME promoter. RESULTS: Among 678 samples analyzed by qPCR, 498 (73.5%) were PRAME- and 180 (26.5%) were PRAME+. Class 1 tumors were more likely to be PRAME-, whereas Class 2 tumors were more likely to be PRAME+ (P < 0.0001). PRAME expression was associated with shorter time to metastasis and melanoma specific mortality in Class 2 tumors (P = 0.01 and P = 0.02, respectively). In Class 1 tumors, PRAME expression was directly associated with SF3B1 mutations (P < 0.0001) and inversely associated with EIF1AX mutations (P = 0.004). PRAME expression was strongly associated with hypomethylation at 12 CpG sites near the PRAME promoter. MATERIALS AND METHODS: Analyses included PRAME mRNA expression, Class 1 versus Class 2 status, chromosomal copy number, mutation status of BAP1, EIF1AX, GNA11, GNAQ and SF3B1, and genomic DNA methylation status. Analyses were performed on 555 de-identified samples from Castle Biosciences, 123 samples from our center, and 80 samples from the TCGA. CONCLUSIONS: PRAME is aberrantly hypomethylated and activated in Class 1 and Class 2 uveal melanomas and is associated with increased metastatic risk in both classes. Since PRAME has been successfully targeted for immunotherapy, it may prove to be a companion prognostic biomarker.