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.

ChIPprimersDB: a public repository of verified qPCR primers for chromatin immunoprecipitation (ChIP).

Chromatin immunoprecipitation (ChIP) has ushered in a new era of scientific discovery by allowing new insights into DNA-protein interactions. ChIP is used to quantify enriched genomic regions using qPCR, and more recently is combined with next generation sequencing (ChIP-seq) to obtain a genome wide profile of protein binding sites. Nevertheless, ChIP-qPCR remains an integral component of this technology for quality control purposes, before the library preparation and sequencing steps. In addition, ChIP-qPCR remains more time- and cost-effective for many focused projects in which the DNA regions of interest are already known. However, the DNA oligonucleotide primers needed for ChIP-qPCR are more challenging to design than for other qPCR projects. Here, we present the first public repository for ChIP oligonucleotides that have been verified to perform well in ChIP-qPCR experiments. ChIPprimersDB was developed by manual screening of publications to ensure primer quality and provide additional specific information on the ChIP experiments where the primers have been used. In addition to the primer sequences, the database includes information about the antibody, cells and tissues used in the experiment, information on the experimental design, and a direct link to the original publication. The database is linked at https://umiamihealth.org/bascom-palmer-eye-institute/research/clinical-and-laboratory-research/ocular-oncology-laboratory/chip-primers and hosted at https://www.chipprimers.com/.

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.

Impact of the Usher syndrome on olfaction.

Usher syndrome is a genetically and clinically heterogeneous disease in humans, characterized by sensorineural hearing loss, retinitis pigmentosa and vestibular dysfunction. This disease is caused by mutations in genes encoding proteins that form complex networks in different cellular compartments. Currently, it remains unclear whether the Usher proteins also form networks within the olfactory epithelium (OE). Here, we describe Usher gene expression at the mRNA and protein level in the OE of mice and showed interactions between these proteins and olfactory signaling proteins. Additionally, we analyzed the odor sensitivity of different Usher syndrome mouse models using electro-olfactogram recordings and monitored significant changes in the odor detection capabilities in mice expressing mutant Usher proteins. Furthermore, we observed changes in the expression of signaling proteins that might compensate for the Usher protein deficiency. In summary, this study provides novel insights into the presence and purpose of the Usher proteins in olfactory signal transduction.

Isolation, culture optimization and functional characterization of stem cell neurospheres from mouse neonatal olfactory bulb and epithelium.

The olfactory epithelium contains basal cells with stem cell characteristics, which have the capacity to differentiate throughout life into olfactory receptor neurons (ORNs). Here we investigate the in vitro characteristics of stem cells taken from the olfactory bulb (OB) and the olfactory epithelium (OE) of neonatal TIS21 knock-in mice. The major aim of the study was the generation of olfactory neurospheres (ONS) derived from OB and OE of neonatal mice as a tool to further analyze the elementary processes of ORN development. Our data showed that the presence of epidermal growth factor (EGF) and fibroblast growth factor (FGF) leads to a significant increase in number of ONS derived from OB but not from OE. The differentiation of ONSs led to the formation of different neuronal cell types, in particular to bipolar-shaped cells as well as putative pyramidal-neurons, astrocytes and oligodendrocytes. Immunohistochemical staining confirmed the presence of astrocytes and neurons in both types of ONSs. In order to investigate the functionality of the neurons we performed calcium imaging and patch-clamp experiments. Calcium imaging experiments revealed that the application of high potassium concentration provokes calcium transients. No excitable properties, neither sodium currents nor action potentials, were observed for the bipolar-shaped cells derived from OB and OE neurospheres, which means that these types of cells morphologically defined as putative neuronal cells, were not physiologically active. Interestingly, patch-clamp recordings performed in the pyramidal-shaped cells of OB neurospheres showed sodium and potassium currents as well as action potentials. Our study will help to establish further models in the field of olfactology.

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.

The Inflammasome-Dependent Dysfunction and Death of Retinal Ganglion Cells after Repetitive Intraocular Pressure Spikes.

The dysfunction and selective loss of retinal ganglion cells (RGCs) is a known cause of vision loss in glaucoma and other neuropathies, where ocular hypertension (OHT) is the major risk factor. We investigated the impact of transient non-ischemic OHT spikes (spOHT) on RGC function and viability in vivo to identify cellular pathways linking low-grade repetitive mechanical stress to RGC pathology. We found that repetitive spOHT had an unexpectedly high impact on intraocular homeostasis and RGC viability, while exposure to steady OHT (stOHT) of a similar intensity and duration failed to induce pathology. The repetitive spOHT induced the rapid activation of the inflammasome, marked by the upregulation of NLRP1, NLRP3, AIM2, caspases -1, -3/7, -8, and Gasdermin D (GSDMD), and the release of interleukin-1ß (IL-1ß) and other cytokines into the vitreous. Similar effects were also detected after 5 weeks of exposure to chronic OHT in an induced glaucoma model. The onset of these immune responses in both spOHT and glaucoma models preceded a 50% deficit in pattern electroretinogram (PERG) amplitude and a significant loss of RGCs 7 days post-injury. The inactivation of inflammasome complexes in Nlrp1-/-, Casp1-/-, and GsdmD-/- knockout animals significantly suppressed the spOHT-induced inflammatory response and protected RGCs. Our results demonstrate that mechanical stress produced by acute repetitive spOHT or chronic OHT is mechanistically linked to inflammasome activation, which leads to RGC dysfunction and death.

Co-Targeting FASN and mTOR Suppresses Uveal Melanoma Growth.

Uveal melanoma (UM) displays a high frequency of metastasis; however, effective therapies for metastatic UM are limited. Identifying unique metabolic features of UM may provide a potential targeting strategy. A lipid metabolism protein expression signature was induced in a normal choroidal melanocyte (NCM) line transduced with GNAQ (Q209L), a driver in UM growth and development. Consistently, UM cells expressed elevated levels of fatty acid synthase (FASN) compared to NCMs. FASN upregulation was associated with increased mammalian target of rapamycin (mTOR) activation and sterol regulatory element-binding protein 1 (SREBP1) levels. FASN and mTOR inhibitors alone significantly reduced UM cell growth. Concurrent inhibition of FASN and mTOR further reduced UM cell growth by promoting cell cycle arrest and inhibiting glucose utilization, TCA cycle metabolism, and de novo fatty acid biosynthesis. Our findings indicate that FASN is important for UM cell growth and co-inhibition of FASN and mTOR signaling may be considered for treatment of UM.

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.

Early Mechanisms of Chemoresistance in Retinoblastoma.

Retinoblastoma is the most common eye cancer in children and is fatal if left untreated. Over the past three decades, chemotherapy has become the mainstay of eye-sparing treatment. Nevertheless, chemoresistance continues to represent a major challenge leading to ocular and systemic toxicity, vision loss, and treatment failure. Unfortunately, the mechanisms leading to chemoresistance remain incompletely understood. Here, we engineered low-passage human retinoblastoma cells to study the early molecular mechanisms leading to resistance to carboplatin, one of the most widely used agents for treating retinoblastoma. Using single-cell next-generation RNA sequencing (scRNA-seq) and single-cell barcoding technologies, we found that carboplatin induced rapid transcriptomic reprogramming associated with the upregulation of PI3K-AKT pathway targets, including ABC transporters and metabolic regulators. Several of these targets are amenable to pharmacologic inhibition, which may reduce the emergence of chemoresistance. We provide evidence to support this hypothesis using a third-generation inhibitor of the ABCB1 transporter.

BAP1 Loss Promotes Suppressive Tumor Immune Microenvironment via Upregulation of PROS1 in Class 2 Uveal Melanomas.

Uveal melanoma (UM) is the most common primary cancer of the eye and is associated with a high rate of metastatic death. UM can be stratified into two main classes based on metastatic risk, with class 1 UM having a low metastatic risk and class 2 UM having a high metastatic risk. Class 2 UM have a distinctive genomic, transcriptomic, histopathologic, and clinical phenotype characterized by biallelic inactivation of the BAP1 tumor-suppressor gene, an immune-suppressive microenvironment enriched for M2-polarized macrophages, and poor response to checkpoint-inhibitor immunotherapy. To identify potential mechanistic links between BAP1 loss and immune suppression in class 2 UM, we performed an integrated analysis of UM samples, as well as genetically engineered UM cell lines and uveal melanocytes (UMC). Using RNA sequencing (RNA-seq), we found that the most highly upregulated gene associated with BAP1 loss across these datasets was PROS1, which encodes a ligand that triggers phosphorylation and activation of the immunosuppressive macrophage receptor MERTK. The inverse association between BAP1 and PROS1 in class 2 UM was confirmed by single-cell RNA-seq, which also revealed that MERTK was upregulated in CD163+ macrophages in class 2 UM. Using ChIP-seq, BAP1 knockdown in UM cells resulted in an accumulation of H3K27ac at the PROS1 locus, suggesting epigenetic regulation of PROS1 by BAP1. Phosphorylation of MERTK in RAW 264.7 monocyte-macrophage cells was increased upon coculture with BAP1-/- UMCs, and this phosphorylation was blocked by depletion of PROS1 in the UMCs. These findings were corroborated by multicolor immunohistochemistry, where class 2/BAP1-mutant UMs demonstrated increased PROS1 expression in tumor cells and increased MERTK phosphorylation in CD163+ macrophages compared with class 1/BAP1-wildtype UMs. Taken together, these findings provide a mechanistic link between BAP1 loss and the suppression of the tumor immune microenvironment in class 2 UMs, and they implicate the PROS1-MERTK pathway as a potential target for immunotherapy in UM.

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.

Molecular evolution of a chordate specific family of G protein-coupled receptors.

BACKGROUND: Chordate evolution is a history of innovations that is marked by physical and behavioral specializations, which led to the development of a variety of forms from a single ancestral group. Among other important characteristics, vertebrates obtained a well developed brain, anterior sensory structures, a closed circulatory system and gills or lungs as blood oxygenation systems. The duplication of pre-existing genes had profound evolutionary implications for the developmental complexity in vertebrates, since mutations modifying the function of a duplicated protein can lead to novel functions, improving the evolutionary success. RESULTS: We analyzed here the evolution of the GPRC5 family of G protein-coupled receptors by comprehensive similarity searches and found that the receptors are only present in chordates and that the size of the receptor family expanded, likely due to genome duplication events in the early history of vertebrate evolution. We propose that a single GPRC5 receptor coding gene originated in a stem chordate ancestor and gave rise by duplication events to a gene family comprising three receptor types (GPRC5A-C) in vertebrates, and a fourth homologue present only in mammals (GPRC5D). Additional duplications of GPRC5B and GPRC5C sequences occurred in teleost fishes. The finding that the expression patterns of the receptors are evolutionarily conserved indicates an important biological function of these receptors. Moreover, we found that expression of GPRC5B is regulated by vitamin A in vivo, confirming previous findings that linked receptor expression to retinoic acid levels in tumor cell lines and strengthening the link between the receptor expression and the development of a complex nervous system in chordates, known to be dependent on retinoic acid signaling. CONCLUSIONS: GPRC5 receptors, a class of G protein-coupled receptors with unique sequence characteristics, may represent a molecular novelty that helped non-chordates to become chordates.

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.

Dual Screen for Efficacy and Toxicity Identifies HDAC Inhibitor with Distinctive Activity Spectrum for BAP1-Mutant Uveal Melanoma.

Drug screens leading to successful targeted therapies in cancer have been mainly based on cell viability assays identifying inhibitors of dominantly acting oncogenes. In contrast, there has been little success in discovering targeted therapies that reverse the effects of inactivating mutations in tumor-suppressor genes. BAP1 is one such tumor suppressor that is frequently inactivated in a variety of cancers, including uveal melanoma, renal cell carcinoma, and mesothelioma. Because BAP1 is an epigenetic transcriptional regulator of developmental genes, we designed a two-phase drug screen involving a cell-based rescue screen of transcriptional repression caused by BAP1 loss, followed by an in vivo screen of lead compounds for rescue of a BAP1-deficient phenotype with minimal toxicity in Xenopus embryos. The first screen identified 9 compounds, 8 of which were HDAC inhibitors. The second screen eliminated all except one compound due to inefficacy or toxicity. The resulting lead compound, quisinostat, has a distinctive activity spectrum, including high potency against HDAC4, which was recently shown to be a key target of BAP1. Quisinostat was further validated in a mouse model and found to prevent the growth of BAP1-mutant uveal melanomas. This innovative strategy demonstrates the potential for identifying therapeutic compounds that target tumor-suppressor mutations in cancer. IMPLICATIONS: Few drugs have been identified that target mutations in tumor suppressors. Using a novel 2-step screening approach, strategy, we identified quisinostat as a candidate for therapy in BAP1-mutant uveal melanoma. HDAC4 is implicated as a key target in uveal melanoma and perhaps other BAP1-mutant cancers.

Establishment and Characterization of a Novel Human Ocular Adnexal Sebaceous Carcinoma Cell Line.

Purpose: Sebaceous carcinoma (SC) is a malignant eyelid tumor of the ocular adnexa that is primarily treated via surgical excision. Few therapies exist in advanced cases, and medical therapy is limited because of our incomplete understanding of SC biology. Herein, we describe a technique to culture human ocular adnexal SC for use as an in vitro model. Methods: Human ocular adnexal SC tumor cells were isolated from a patient undergoing orbital exenteration surgery and named Bascom Palmer 50 (BP50). They were cultured in Dulbecco's modified Eagle medium/nutrient mixture F-12 supplemented with 10% fetal bovine serum and antibiotics and were maintained at 37°C in humidified 5% CO2. The cells were characterized by immunohistochemistry, exome sequencing, and short tandem repeats analysis. In vitro drug screening against mitomycin-C (MMC) was performed using a cell viability assay. Results: BP50 grew past 40 passages with a doubling time of 52.3 hours. Immunocytochemical staining revealed expression of SC-associated markers adipophilin, epithelial membrane antigen, p53, and androgen receptor. Whole exome sequencing showed a significant carryover in somatic mutations between the tumor tissue and corresponding cell line, revealing genetic markers consistent with SC. MMC affected cell viability in a dose-dependent manner. Conclusions: BP50 displays characteristics of ocular adnexal SC and therefore may facilitate improved understanding of SC biology and the high throughput assessment of novel therapeutic compounds and new drug combinatorial approaches targeted for this disease. Translational Relevance: Drug screening with MMC against these cells shows in vitro evidence to support its continued clinical use in SC.

Tmem16b is specifically expressed in the cilia of olfactory sensory neurons.

Calcium-activated chloride channels (CaCCs) are involved in many physiological processes, including sensory signal transduction, but only little is known to date about their structure and function. We performed a proteome analysis of the olfactory epithelium (OE) membrane proteome and identified so far uncharacterized membrane proteins as candidate channels. One of the most abundant membrane proteins in olfactory sensory neurons (OSNs) was Tmem16b, a member of a recently identified family of CaCCs. In addition to former studies performed on Tmem16b, we show here that Tmem16b expression is highly specific for the OE, in contrast to the closely related Tmem16a, which shows a broad expression pattern in secretory epithelial cells. Native Tmem16b is localized in the cilia of the OSNs, which is in agreement with previous electrophysiological recordings.

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.

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.