AtlasXplore: a web platform for visualizing and sharing spatial epigenome data.

MOTIVATION: In recent years, a growing number of spatial epigenome datasets have been generated, presenting rich opportunities for studying the regulation mechanisms in solid tissue sections. However, visual exploration of these datasets requires extensive computational processing of raw data, presenting a challenge for researchers without advanced computational skills to fully explore and analyze such datasets. RESULTS: Here, we introduce AtlasXplore, a web-based platform that enables scientists to interactively navigate a growing collection of spatial epigenome data using an expanding set of tools. AVAILABILITY AND IMPLEMENTATION: https://web.atlasxomics.com.

Single-cell RNA-seq reveals developmental plasticity with coexisting oncogenic states and immune evasion programs in ETP-ALL.

Lineage plasticity and stemness have been invoked as causes of therapy resistance in cancer, because these flexible states allow cancer cells to dedifferentiate and alter their dependencies. We investigated such resistance mechanisms in relapsed/refractory early T-cell progenitor acute lymphoblastic leukemia (ETP-ALL) carrying activating NOTCH1 mutations via full-length single-cell RNA sequencing (scRNA-seq) of malignant and microenvironmental cells. We identified 2 highly distinct stem-like states that critically differed with regard to cell cycle and oncogenic signaling. Fast-cycling stem-like leukemia cells demonstrated Notch activation and were effectively eliminated in patients by Notch inhibition, whereas slow-cycling stem-like cells were Notch independent and rather relied on PI3K signaling, likely explaining the poor efficacy of Notch inhibition in this disease. Remarkably, we found that both stem-like states could differentiate into a more mature leukemia state with prominent immunomodulatory functions, including high expression of the LGALS9 checkpoint molecule. These cells promoted an immunosuppressive leukemia ecosystem with clonal accumulation of dysfunctional CD8+ T cells that expressed HAVCR2, the cognate receptor for LGALS9. Our study identified complex interactions between signaling programs, cellular plasticity, and immune programs that characterize ETP-ALL, illustrating the multidimensionality of tumor heterogeneity. In this scenario, combination therapies targeting diverse oncogenic states and the immune ecosystem seem most promising to successfully eliminate tumor cells that escape treatment through coexisting transcriptional programs.

GPR6 Structural Insights: Homology Model Construction and Docking Studies.

GPR6 is an orphan G protein-coupled receptor that has been associated with the cannabinoid family because of its recognition of a sub-set of cannabinoid ligands. The high abundance of GPR6 in the central nervous system, along with high constitutive activity and a link to several neurodegenerative diseases make GPR6 a promising biological target. In fact, diverse research groups have demonstrated that GPR6 represents a possible target for the treatment of neurodegenerative disorders such as Parkinson's disease, Alzheimer's disease, and Huntington's disease. Several patents have claimed the use of a wide range of pyrazine derivatives as GPR6 inverse agonists for the treatment of Parkinson's disease symptoms and other dyskinesia syndromes. However, the full pharmacological importance of GPR6 has not yet been fully explored due to the lack of high potency, readily available ligands targeting GPR6. The long-term goal of the present study is to develop such ligands. In this paper, we describe our initial steps towards this goal. A human GPR6 homology model was constructed using a suite of computational techniques. This model permitted the identification of unique GPR6 structural features and the exploration of the GPR6 binding crevice. A subset of patented pyrazine analogs were docked in the resultant GPR6 inactive state model to validate the model, rationalize the structure-activity relationships from the reported patents and identify the key residues in the binding crevice for ligand recognition. We will take this structural knowledge into the next phase of GPR6 project, in which scaffold hopping will be used to design new GPR6 ligands.

Towards A Molecular Understanding of The Cannabinoid Related Orphan Receptor GPR18: A Focus on Its Constitutive Activity.

The orphan G-protein coupled receptor (GPCR), GPR18, has been recently proposed as a potential member of the cannabinoid family as it recognizes several endogenous, phytogenic, and synthetic cannabinoids. Potential therapeutic applications for GPR18 include intraocular pressure, metabolic disorders, and cancer. GPR18 has been reported to have high constitutive activity, i.e., activation/signaling occurs in the absence of an agonist. This activity can be reduced significantly by the A3.39N mutation. At the intracellular (IC) ends of (transmembrane helices) TMH3 and TMH6 in GPCRs, typically, a pair of oppositely charged amino acids form a salt bridge called the "ionic lock". Breaking of this salt bridge creates an IC opening for coupling with G protein. The GPR18 "ionic lock" residues (R3.50/S6.33) can form only a hydrogen bond. In this paper, we test the hypothesis that the high constitutive activity of GPR18 is due to the weakness of its "ionic lock" and that the A3.39N mutation strengthens this lock. To this end, we report molecular dynamics simulations of wild-type (WT) GPR18 and the A3.39N mutant in fully hydrated (POPC) phophatidylcholine lipid bilayers. Results suggest that in the A3.39N mutant, TMH6 rotates and brings R3.50 and S6.33 closer together, thus strengthening the GPR18 "ionic lock".

Regulatory Programs of B-cell Activation and Germinal Center Reaction Allow B-ALL Escape from CD19 CAR T-cell Therapy.

Chimeric antigen receptor (CAR) T-cell therapy has led to tremendous successes in the treatment of B-cell malignancies. However, a large fraction of treated patients relapse, often with disease expressing reduced levels of the target antigen. Here, we report that exposing CD19+ B-cell acute lymphoblastic leukemia (B-ALL) cells to CD19 CAR T cells reduced CD19 expression within hours. Initially, CD19 CAR T cells caused clustering of CD19 at the T cell-leukemia cell interface followed by CD19 internalization and decreased CD19 surface expression on the B-ALL cells. CD19 expression was then repressed by transcriptional rewiring. Using single-cell RNA sequencing and single-cell assay for transposase-accessible chromatin using sequencing, we demonstrated that a subset of refractory CD19low cells sustained decreased CD19 expression through transcriptional programs of physiologic B-cell activation and germinal center reaction. Inhibiting B-cell activation programs with the Bruton's tyrosine kinase inhibitor ibrutinib increased the cytotoxicity of CD19 CAR T cells without affecting CAR T-cell viability. These results demonstrate transcriptional plasticity as an underlying mechanism of escape from CAR T cells and highlight the importance of combining CAR T-cell therapy with targeted therapies that aim to overcome this plasticity. See related Spotlight by Zhao and Melenhorst, p. 1040.

Dynamic transcriptional reprogramming leads to immunotherapeutic vulnerabilities in myeloma.

While there is extensive evidence for genetic variation as a basis for treatment resistance, other sources of variation result from cellular plasticity. Using multiple myeloma as an example of an incurable lymphoid malignancy, we show how cancer cells modulate lineage restriction, adapt their enhancer usage and employ cell-intrinsic diversity for survival and treatment escape. By using single-cell transcriptome and chromatin accessibility profiling, we show that distinct transcriptional states co-exist in individual cancer cells and that differential transcriptional regulon usage and enhancer rewiring underlie these alternative transcriptional states. We demonstrate that exposure to standard treatment further promotes transcriptional reprogramming and differential enhancer recruitment while simultaneously reducing developmental potential. Importantly, treatment generates a distinct complement of actionable immunotherapy targets, such as CXCR4, which can be exploited to overcome treatment resistance. Our studies therefore delineate how to transform the cellular plasticity that underlies drug resistance into immuno-oncologic therapeutic opportunities.

Single-Cell Profiling Reveals Metabolic Reprogramming as a Resistance Mechanism in BRAF-Mutated Multiple Myeloma.

PURPOSE: Although remarkably effective in some patients, precision medicine typically induces only transient responses despite initial absence of resistance-conferring mutations. Using BRAF-mutated myeloma as a model for resistance to precision medicine we investigated if BRAF-mutated cancer cells have the ability to ensure their survival by rapidly adapting to BRAF inhibitor treatment. EXPERIMENTAL DESIGN: Full-length single-cell RNA (scRNA) sequencing (scRNA-seq) was conducted on 3 patients with BRAF-mutated myeloma and 1 healthy donor. We sequenced 1,495 cells before, after 1 week, and at clinical relapse to BRAF/MEK inhibitor treatment. We developed an in vitro model of dabrafenib resistance using genetically homogeneous single-cell clones from two cell lines with established BRAF mutations (U266, DP6). Transcriptional and epigenetic adaptation in resistant cells were defined by RNA-seq and H3K27ac chromatin immunoprecipitation sequencing (ChIP-seq). Mitochondrial metabolism was characterized by metabolic flux analysis. RESULTS: Profiling by scRNA-seq revealed rapid cellular state changes in response to BRAF/MEK inhibition in patients with myeloma and cell lines. Transcriptional adaptation preceded detectable outgrowth of genetically discernible drug-resistant clones and was associated with widespread enhancer remodeling. As a dominant vulnerability, dependency on oxidative phosphorylation (OxPhos) was induced. In treated individuals, OxPhos was activated at the time of relapse and showed inverse correlation to MAPK activation. Metabolic flux analysis confirmed OxPhos as a preferential energetic resource of drug-persistent myeloma cells. CONCLUSIONS: This study demonstrates that cancer cells have the ability to rapidly adapt to precision treatments through transcriptional state changes, epigenetic adaptation, and metabolic rewiring, thus facilitating the development of refractory disease while simultaneously exposing novel vulnerabilities.

Therapeutic Exploitation of GPR18: Beyond the Cannabinoids?

GPR18 is a G-protein-coupled receptor that belongs to the orphan class A family. Even though it shares low sequence homology with the cannabinoid receptors CB1R and CB2R, a growing body of research suggests its relationship with the endocannabinoid system, not only because it is able to recognize cannabinoid ligands but also because of its expression and ability to heteromerize with CBRs. In this review, we aim to analyze the biological relevance, reported modulators, and structural features of GPR18. In order to guide future drug design in this field, highlights from molecular modeling of GPR18 will be provided.