Unraveling the Complexity of Liver Disease One Cell at a Time

The human liver is a complex organ made up of multiple specialized cell types that carry out key physiological functions. An incomplete understanding of liver biology limits our ability to develop therapeutics to prevent chronic liver diseases, liver cancers, and death as a result of organ failure. Recently, single-cell modalities have expanded our understanding of the cellular phenotypic heterogeneity and intercellular cross-talk in liver health and disease. This review summarizes these findings and looks forward to highlighting new avenues for the application of single-cell genomics to unravel unknown pathogenic pathways and disease mechanisms for the development of new therapeutics targeting liver pathology. As these technologies mature, their integration into clinical data analysis will aid in patient stratification and in developing treatment plans for patients suffering from liver disease.

Uncovering the Contribution of Moderate-Penetrance Susceptibility Genes to Breast Cancer by Whole-Exome Sequencing and Targeted Enrichment Sequencing of Candidate Genes in Women of European Ancestry

Rare variants in at least 10 genes, including BRCA1, BRCA2, PALB2, ATM, and CHEK2, are associated with increased risk of breast cancer; however, these variants, in combination with common variants identified through genome-wide association studies, explain only a fraction of the familial aggregation of the disease. To identify further susceptibility genes, we performed a two-stage whole-exome sequencing study. In the discovery stage, samples from 1528 breast cancer cases enriched for breast cancer susceptibility and 3733 geographically matched unaffected controls were sequenced. Using five different filtering and gene prioritization strategies, 198 genes were selected for further validation. These genes, and a panel of 32 known or suspected breast cancer susceptibility genes, were assessed in a validation set of 6211 cases and 6019 controls for their association with risk of breast cancer overall, and by estrogen receptor (ER) disease subtypes, using gene burden tests applied to loss-of-function and rare missense variants. Twenty genes showed nominal evidence of association (p-value < 0.05) with either overall or subtype-specific breast cancer. Our study had the statistical power to detect susceptibility genes with effect sizes similar to ATM, CHEK2, and PALB2, however, it was underpowered to identify genes in which susceptibility variants are rarer or confer smaller effect sizes. Larger sample sizes would be required in order to identify such genes.

Comprehensive multi-cohort transcriptional meta-analysis of muscle diseases identifies a signature of disease severity

Muscle diseases share common pathological features suggesting common underlying mechanisms. We hypothesized there is a common set of genes dysregulated across muscle diseases compared to healthy muscle and that these genes correlate with severity of muscle disease. We performed meta-analysis of transcriptional profiles of muscle biopsies from human muscle diseases and healthy controls. Studies obtained from public microarray repositories fulfilling quality criteria were divided into six categories: (i) immobility, (ii) inflammatory myopathies, (iii) intensive care unit (ICU) acquired weakness (ICUAW), (iv) congenital muscle diseases, (v) chronic systemic diseases, (vi) motor neuron disease. Patient cohorts were separated in discovery and validation cohorts retaining roughly equal proportions of samples for the disease categories. To remove bias towards a specific muscle disease category we repeated the meta-analysis five times by removing data sets corresponding to one muscle disease class at a time in a "leave-one-disease-out" analysis. We used 636 muscle tissue samples from 30 independent cohorts to identify a 52 gene signature (36 up-regulated and 16 down-regulated genes). We validated the discriminatory power of this signature in 657 muscle biopsies from 12 additional patient cohorts encompassing five categories of muscle diseases with an area under the receiver operating characteristic curve of 0.91, 83% sensitivity, and 85.3% specificity. The expression score of the gene signature inversely correlated with quadriceps muscle mass (r = -0.50, p-value = 0.011) in ICUAW and shoulder abduction strength (r = -0.77, p-value = 0.014) in amyotrophic lateral sclerosis (ALS). The signature also positively correlated with histologic assessment of muscle atrophy in ALS (r = 0.88, p-value = 1.62 × 10^-3) and fibrosis in muscular dystrophy (Jonckheere trend test p-value = 4.45 × 10^-9). Our results identify a conserved transcriptional signature associated with clinical and histologic muscle disease severity. Several genes in this conserved signature have not been previously associated with muscle disease severity.

The metabolic enzyme hexokinase 2 localizes to the nucleus in AML and normal haematopoietic stem and progenitor cells to maintain stemness

Mitochondrial metabolites regulate leukaemic and normal stem cells by affecting epigenetic marks. How mitochondrial enzymes localize to the nucleus to control stem cell function is less understood. We discovered that the mitochondrial metabolic enzyme hexokinase 2 (HK2) localizes to the nucleus in leukaemic and normal haematopoietic stem cells. Overexpression of nuclear HK2 increases leukaemic stem cell properties and decreases differentiation, whereas selective nuclear HK2 knockdown promotes differentiation and decreases stem cell function. Nuclear HK2 localization is phosphorylation-dependent, requires active import and export, and regulates differentiation independently of its enzymatic activity. HK2 interacts with nuclear proteins regulating chromatin openness, increasing chromatin accessibilities at leukaemic stem cell-positive signature and DNA-repair sites. Nuclear HK2 overexpression decreases double-strand breaks and confers chemoresistance, which may contribute to the mechanism by which leukaemic stem cells resist DNA-damaging agents. Thus, we describe a non-canonical mechanism by which mitochondrial enzymes influence stem cell function independently of their metabolic function.

Thomas, Egan et al. report that hexokinase 2 localizes to the nucleus of leukaemic and normal haematopoietic cells to maintain stemness by interacting with nuclear proteins and modulating chromatin accessibility independently of its kinase activity.

Identification of the global miR-130a targetome reveals a role for TBL1XR1 in hematopoietic stem cell self-renewal and t(8;21) AML

Gene expression profiling and proteome analysis of normal and malignant hematopoietic stem cells (HSCs) point to shared core stemness properties. However, discordance between mRNA and protein signatures highlights an important role for post-transcriptional regulation by microRNAs (miRNAs) in governing this critical nexus. Here, we identify miR-130a as a regulator of HSC self-renewal and differentiation. Enforced expression of miR-130a impairs B lymphoid differentiation and expands long-term HSCs. Integration of protein mass spectrometry and chimeric AGO2 crosslinking and immunoprecipitation (CLIP) identifies TBL1XR1 as a primary miR-130a target, whose loss of function phenocopies miR-130a overexpression. Moreover, we report that miR-130a is highly expressed in t(8;21) acute myeloid leukemia (AML), where it is critical for maintaining the oncogenic molecular program mediated by the AML1-ETO complex. Our study establishes that identification of the comprehensive miRNA targetome within primary cells enables discovery of genes and molecular networks underpinning stemness properties of normal and leukemic cells.

IPO11 regulates the nuclear import of BZW1/2 and is necessary for AML cells and stem cells

AML cells are arranged in a hierarchy with stem/progenitor cells giving rise to more differentiated bulk cells. Despite the importance of stem/progenitors in the pathogenesis of AML, the determinants of the AML stem/progenitor state are not fully understood. Through a comparison of genes that are significant for growth and viability of AML cells by way of a CRISPR screen, with genes that are differentially expressed in leukemia stem cells (LSC), we identified importin 11 (IPO11) as a novel target in AML. Importin 11 (IPO11) is a member of the importin β family of proteins that mediate transport of proteins across the nuclear membrane. In AML, knockdown of IPO11 decreased growth, reduced engraftment potential of LSC, and induced differentiation. Mechanistically, we identified the transcription factors BZW1 and BZW2 as novel cargo of IPO11. We further show that BZW1/2 mediate a transcriptional signature that promotes stemness and survival of LSC. Thus, we demonstrate for the first time how specific cytoplasmic-nuclear regulation supports stem-like transcriptional signature in relapsed AML.

Hypophosphorylated pRb knock-in mice exhibit hallmarks of aging and vitamin C-preventable diabetes

Despite extensive analysis of pRB phosphorylation in vitro, how this modification influences development and homeostasis in vivo is unclear. Here, we show that homozygous Rb^∆K4 and Rb^∆K7 knock‐in mice, in which either four or all seven phosphorylation sites in the C‐terminal region of pRb, respectively, have been abolished by Ser/Thr‐to‐Ala substitutions, undergo normal embryogenesis and early development, notwithstanding suppressed phosphorylation of additional upstream sites. Whereas Rb^∆K4 mice exhibit telomere attrition but no other abnormalities, Rb^∆K7 mice are smaller and display additional hallmarks of premature aging including infertility, kyphosis, and diabetes, indicating an accumulative effect of blocking pRb phosphorylation. Diabetes in Rb^∆K7 mice is insulin‐sensitive and associated with failure of quiescent pancreatic β‐cells to re‐enter the cell cycle in response to mitogens, resulting in induction of DNA damage response (DDR), senescence‐associated secretory phenotype (SASP), and reduced pancreatic islet mass and circulating insulin level. Pre‐treatment with the epigenetic regulator vitamin C reduces DDR, increases cell cycle re‐entry, improves islet morphology, and attenuates diabetes. These results have direct implications for cell cycle regulation, CDK‐inhibitor therapeutics, diabetes, and longevity.

Temporal profiling of therapy resistance in human medulloblastoma identifies novel targetable drivers of recurrence

Medulloblastoma (MB) remains a leading cause of cancer-related mortality among children. The paucity of MB samples collected at relapse has hindered the functional understanding of molecular mechanisms driving therapy failure. New models capable of accurately recapitulating tumor progression in response to conventional therapeutic interventions are urgently needed. In this study, we developed a therapy-adapted PDX MB model that has a distinct advantage of generating human MB recurrence. The comparative gene expression analysis of MB cells collected throughout therapy led to identification of genes specifically up-regulated after therapy, including one previously undescribed in the setting of brain tumors, bactericidal/permeability-increasing fold-containing family B member 4 (BPIFB4). Subsequent functional validation resulted in a markedly diminished in vitro proliferation, self-renewal, and longevity of MB cells, translating into extended survival and reduced tumor burden in vivo. Targeting endothelial nitric oxide synthase, a downstream substrate of BPIFB4, impeded growth of several patient-derived MB lines at low nanomolar concentrations.

Author-sourced capture of pathway knowledge in computable form using Biofactoid

Making the knowledge contained in scientific papers machine-readable and formally computable would allow researchers to take full advantage of this information by enabling integration with other knowledge sources to support data analysis and interpretation. Here we describe Biofactoid, a web-based platform that allows scientists to specify networks of interactions between genes, their products, and chemical compounds, and then translates this information into a representation suitable for computational analysis, search and discovery. We also report the results of a pilot study to encourage the wide adoption of Biofactoid by the scientific community.

Single-Cell, Single-Nucleus, and Spatial RNA Sequencing of the Human Liver Identifies Cholangiocyte and Mesenchymal Heterogeneity

The critical functions of the human liver are coordinated through the interactions of hepatic parenchymal and non-parenchymal cells. Recent advances in single-cell transcriptional approaches have enabled an examination of the human liver with unprecedented resolution. However, dissociation-related cell perturbation can limit the ability to fully capture the human liver's parenchymal cell fraction, which limits the ability to comprehensively profile this organ. Here, we report the transcriptional landscape of 73,295 cells from the human liver using matched single-cell RNA sequencing (scRNA-seq) and single-nucleus RNA sequencing (snRNA-seq). The addition of snRNA-seq enabled the characterization of interzonal hepatocytes at a single-cell resolution, revealed the presence of rare subtypes of liver mesenchymal cells, and facilitated the detection of cholangiocyte progenitors that had only been observed during in vitro differentiation experiments. However, T and B lymphocytes and natural killer cells were only distinguishable using scRNA-seq, highlighting the importance of applying both technologies to obtain a complete map of tissue-resident cell types. We validated the distinct spatial distribution of the hepatocyte, cholangiocyte, and mesenchymal cell populations by an independent spatial transcriptomics data set and immunohistochemistry. Conclusion: Our study provides a systematic comparison of the transcriptomes captured by scRNA-seq and snRNA-seq and delivers a high-resolution map of the parenchymal cell populations in the healthy human liver.

A microfluidic platform enables comprehensive gene expression profiling of mouse retinal stem cells

Loss of photoreceptors due to retinal degeneration is a major cause of untreatable visual impairment and blindness. Cell replacement therapy, using retinal stem cell (RSC)-derived photoreceptors, holds promise for reconstituting damaged cell populations in the retina. One major obstacle preventing translation to the clinic is the lack of validated markers or strategies to prospectively identify these rare cells in the retina and subsequently enrich them. Here, we introduce a microfluidic platform that combines nickel micromagnets, herringbone structures, and a design enabling varying flow velocities among three compartments to facilitate a highly efficient enrichment of RSCs. In addition, we developed an affinity enrichment strategy based on cell-surface markers that was utilized to isolate RSCs from the adult ciliary epithelium. We showed that targeting a panel of three cell surface markers simultaneously facilitates the enrichment of RSCs to 1 : 3 relative to unsorted cells. Combining the microfluidic platform with single-cell whole-transcriptome profiling, we successfully identified four differentially expressed cell surface markers that can be targeted simultaneously to yield an unprecedented 1 : 2 enrichment of RSCs relative to unsorted cells. We also identified transcription factors (TFs) that play functional roles in maintenance, quiescence, and proliferation of RSCs. This level of analysis for the first time identified a spectrum of molecular and functional properties of RSCs.

A clinically applicable integrative molecular classification of meningiomas

Meningiomas are the most common primary intracranial tumour in adults1. Patients with symptoms are generally treated with surgery as there are no effective medical therapies. The World Health Organization histopathological grade of the tumour and the extent of resection at surgery (Simpson grade) are associated with the recurrence of disease; however, they do not accurately reflect the clinical behaviour of all meningiomas2. Molecular classifications of meningioma that reliably reflect tumour behaviour and inform on therapies are required. Here we introduce four consensus molecular groups of meningioma by combining DNA somatic copy-number aberrations, DNA somatic point mutations, DNA methylation and messenger RNA abundance in a unified analysis. These molecular groups more accurately predicted clinical outcomes compared with existing classification schemes. Each molecular group showed distinctive and prototypical biology (immunogenic, benign NF2 wild-type, hypermetabolic and proliferative) that informed therapeutic options. Proteogenomic characterization reinforced the robustness of the newly defined molecular groups and uncovered highly abundant and group-specific protein targets that we validated using immunohistochemistry. Single-cell RNA sequencing revealed inter-individual variations in meningioma as well as variations in intrinsic expression programs in neoplastic cells that mirrored the biology of the molecular groups identified.

Nicotinamide phosphoribosyltransferase inhibitors selectively induce apoptosis of AML stem cells by disrupting lipid homeostasis

Current treatments for acute myeloid leukemia (AML) are often ineffective in eliminating leukemic stem cells (LSCs), which perpetuate the disease. Here, we performed a metabolic drug screen to identify LSC-specific vulnerabilities and found that nicotinamide phosphoribosyltransferase (NAMPT) inhibitors selectively killed LSCs, while sparing normal hematopoietic stem and progenitor cells. Treatment with KPT-9274, a NAMPT inhibitor, suppressed the conversion of saturated fatty acids to monounsaturated fatty acids, a reaction catalyzed by the stearoyl-CoA desaturase (SCD) enzyme, resulting in apoptosis of AML cells. Transcriptomic analysis of LSCs treated with KPT-9274 revealed an upregulation of sterol regulatory-element binding protein (SREBP)-regulated genes, including SCD, which conferred partial protection against NAMPT inhibitors. Inhibition of SREBP signaling with dipyridamole enhanced the cytotoxicity of KPT-9274 on LSCs in vivo. Our work demonstrates that altered lipid homeostasis plays a key role in NAMPT inhibitor-induced apoptosis and identifies NAMPT inhibition as a therapeutic strategy for targeting LSCs in AML.

Control of tissue development and cell diversity by cell cycle-dependent transcriptional filtering

Cell cycle duration changes dramatically during development, starting out fast to generate cells quickly and slowing down over time as the organism matures. The cell cycle can also act as a transcriptional filter to control the expression of long gene transcripts, which are partially transcribed in short cycles. Using mathematical simulations of cell proliferation, we identify an emergent property that this filter can act as a tuning knob to control gene transcript expression, cell diversity, and the number and proportion of different cell types in a tissue. Our predictions are supported by comparison to single-cell RNA-seq data captured over embryonic development. Additionally, evolutionary genome analysis shows that fast-developing organisms have a narrow genomic distribution of gene lengths while slower developers have an expanded number of long genes. Our results support the idea that cell cycle dynamics may be important across multicellular animals for controlling gene transcript expression and cell fate.

The mitochondrial peptidase, neurolysin, regulates respiratory chain supercomplex formation and is necessary for AML viability

Neurolysin (NLN) is a zinc metallopeptidase whose mitochondrial function is unclear. We found that NLN was overexpressed in almost half of patients with acute myeloid leukemia (AML), and inhibition of NLN was selectively cytotoxic to AML cells and stem cells while sparing normal hematopoietic cells. Mechanistically, NLN interacted with the mitochondrial respiratory chain. Genetic and chemical inhibition of NLN impaired oxidative metabolism and disrupted the formation of respiratory chain supercomplexes (RCS). Furthermore, NLN interacted with the known RCS regulator, LETM1, and inhibition of NLN disrupted LETM1 complex formation. RCS were increased in patients with AML and positively correlated with NLN expression. These findings demonstrate that inhibiting RCS formation selectively targets AML cells and stem cells and highlights the therapeutic potential of pharmacologically targeting NLN in AML.

Surveilling and Tracking COVID-19 Patients Using a Portable Quantum Dot Smartphone Device

The ability to rapidly diagnose, track, and disseminate information for SARS-CoV-2 is critical to minimize its spread. Here, we engineered a portable smartphone-based quantum barcode serological assay device for real-time surveillance of patients infected with SARS-CoV-2. Our device achieved a clinical sensitivity of 90% and specificity of 100% for SARS-CoV-2, as compared to 34% and 100%, respectively, for lateral flow assays in a head-to-head comparison. The lateral flow assay misdiagnosed ∼2 out of 3 SARS-CoV-2 positive patients. Our quantum dot barcode device has ∼3 times greater clinical sensitivity because it is ∼140 times more analytically sensitive than lateral flow assays. Our device can diagnose SARS-CoV-2 at different sampling dates and infectious severity. We developed a databasing app to provide instantaneous results to inform patients, physicians, and public health agencies. This assay and device enable real-time surveillance of SARS-CoV-2 seroprevalence and potential immunity.

scNetViz: from single cells to networks using Cytoscape

Single-cell RNA-sequencing (scRNA-seq) has revolutionized molecular biology and medicine by enabling high-throughput studies of cellular heterogeneity in diverse tissues. Applying network biology approaches to scRNA-seq data can provide useful insights into genes driving heterogeneous cell-type compositions of tissues. Here, we present scNetViz-  a Cytoscape app to aid biological interpretation of cell clusters in scRNA-seq data using network analysis. scNetViz calculates the differential expression of each gene across clusters and then creates a cluster-specific gene functional interaction network between the significantly differentially expressed genes for further analysis, such as pathway enrichment analysis. To automate a complete data analysis workflow, scNetViz integrates parts of the Scanpy software, which is a popular Python package for scRNA-seq data analysis, with Cytoscape apps such as stringApp, cyPlot, and enhancedGraphics. We describe our implementation of methods for accessing data from public single cell atlas projects, differential expression analysis, visualization, and automation. scNetViz enables users to analyze data from public atlases or their own experiments, which we illustrate with two use cases. Analysis can be performed via the Cytoscape GUI or CyREST programming interface using R (RCy3) or Python (py4cytoscape).

Single-Cell Profiling Shows Murine Forebrain Neural Stem Cells Reacquire a Developmental State when Activated for Adult Neurogenesis

The transitions from developing to adult quiescent and activated neural stem cells (NSCs) are not well understood. Here, we use single-cell transcriptional profiling and lineage tracing to characterize these transitions in the murine forebrain. We show that the two forebrain NSC parental populations, embryonic cortex and ganglionic eminence radial precursors (RPs), are highly similar even though they make glutamatergic versus gabaergic neurons. Both RP populations progress linearly to transition from a highly active embryonic to a dormant adult stem cell state that still shares many similarities with embryonic RPs. When adult NSCs of either embryonic origin become reactivated to make gabaergic neurons, they acquire a developing ganglionic eminence RP-like identity. Thus, transitions from embryonic RPs to adult NSCs and back to neuronal progenitors do not involve fundamental changes in cell identity, but rather reflect conversions between activated and dormant NSC states that may be determined by the niche environment.

Generation of mature compact ventricular cardiomyocytes from human pluripotent stem cells

Compact cardiomyocytes that make up the ventricular wall of the adult heart represent an important therapeutic target population for modeling and treating cardiovascular diseases. Here, we established a differentiation strategy that promotes the specification, proliferation and maturation of compact ventricular cardiomyocytes from human pluripotent stem cells (hPSCs). The cardiomyocytes generated under these conditions display the ability to use fatty acids as an energy source, a high mitochondrial mass, well-defined sarcomere structures and enhanced contraction force. These ventricular cells undergo metabolic changes indicative of those associated with heart failure when challenged in vitro with pathological stimuli and were found to generate grafts consisting of more mature cells than those derived from immature cardiomyocytes following transplantation into infarcted rat hearts. hPSC-derived atrial cardiomyocytes also responded to the maturation cues identified in this study, indicating that the approach is broadly applicable to different subtypes of the heart. Collectively, these findings highlight the power of recapitulating key aspects of embryonic and postnatal development for generating therapeutically relevant cell types from hPSCs.

Cardiomyocytes of heart ventricles consist of subpopulations of trabecular and compact subtypes. Here the authors describe the generation of structurally, metabolically and functionally mature compact ventricular cardiomyocytes as well as mature atrial cardiomyocytes from human pluripotent stem cells.

Tutorial: guidelines for annotating single-cell transcriptomic maps using automated and manual methods

Single-cell transcriptomics can profile thousands of cells in a single experiment and identify novel cell types, states and dynamics in a wide variety of tissues and organisms. Standard experimental protocols and analysis workflows have been developed to create single-cell transcriptomic maps from tissues. This tutorial focuses on how to interpret these data to identify cell types, states and other biologically relevant patterns with the objective of creating an annotated map of cells. We recommend a three-step workflow including automatic cell annotation (wherever possible), manual cell annotation and verification. Frequently encountered challenges are discussed, as well as strategies to address them. Guiding principles and specific recommendations for software tools and resources that can be used for each step are covered, and an R notebook is included to help run the recommended workflow. Basic familiarity with computer software is assumed, and basic knowledge of programming (e.g., in the R language) is recommended.

SBGN Bricks Ontology as a tool to describe recurring concepts in molecular networks

A comprehensible representation of a molecular network is key to communicating and understanding scientific results in systems biology. The Systems Biology Graphical Notation (SBGN) has emerged as the main standard to represent such networks graphically. It has been implemented by different software tools, and is now largely used to communicate maps in scientific publications. However, learning the standard, and using it to build large maps, can be tedious. Moreover, SBGN maps are not grounded on a formal semantic layer and therefore do not enable formal analysis. Here, we introduce a new set of patterns representing recurring concepts encountered in molecular networks, called SBGN bricks. The bricks are structured in a new ontology, the Bricks Ontology (BKO), to define clear semantics for each of the biological concepts they represent. We show the usefulness of the bricks and BKO for both the template-based construction and the semantic annotation of molecular networks. The SBGN bricks and BKO can be freely explored and downloaded at sbgnbricks.org.

Distinct fibroblast functional states drive clinical outcomes in ovarian cancer and are regulated by TCF21

Recent studies indicate that cancer-associated fibroblasts (CAFs) are phenotypically and functionally heterogeneous. However, little is known about CAF subtypes, the roles they play in cancer progression, and molecular mediators of the CAF “state.” Here, we identify a novel cell surface pan-CAF marker, CD49e, and demonstrate that two distinct CAF states, distinguished by expression of fibroblast activation protein (FAP), coexist within the CD49e⁺ CAF compartment in high-grade serous ovarian cancers. We show for the first time that CAF state influences patient outcomes and that this is mediated by the ability of FAP-high, but not FAP-low, CAFs to aggressively promote proliferation, invasion and therapy resistance of cancer cells. Overexpression of the FAP-low–specific transcription factor TCF21 in FAP-high CAFs decreases their ability to promote invasion, chemoresistance, and in vivo tumor growth, indicating that it acts as a master regulator of the CAF state. Understanding CAF states in more detail could lead to better patient stratification and novel therapeutic strategies.

PRMT5 inhibition disrupts splicing and stemness in glioblastoma

Glioblastoma (GBM) is a deadly cancer in which cancer stem cells (CSCs) sustain tumor growth and contribute to therapeutic resistance. Protein arginine methyltransferase 5 (PRMT5) has recently emerged as a promising target in GBM. Using two orthogonal-acting inhibitors of PRMT5 (GSK591 or LLY-283), we show that pharmacological inhibition of PRMT5 suppresses the growth of a cohort of 46 patient-derived GBM stem cell cultures, with the proneural subtype showing greater sensitivity. We show that PRMT5 inhibition causes widespread disruption of splicing across the transcriptome, particularly affecting cell cycle gene products. We identify a GBM splicing signature that correlates with the degree of response to PRMT5 inhibition. Importantly, we demonstrate that LLY-283 is brain-penetrant and significantly prolongs the survival of mice with orthotopic patient-derived xenografts. Collectively, our findings provide a rationale for the clinical development of brain penetrant PRMT5 inhibitors as treatment for GBM.

Gradient of Developmental and Injury Response transcriptional states defines functional vulnerabilities underpinning glioblastoma heterogeneity

Glioblastomas harbor diverse cell populations, including rare glioblastoma stem cells (GSCs) that drive tumorigenesis. To characterize functional diversity within this population, we performed single-cell RNA sequencing on >69,000 GSCs cultured from the tumors of 26 patients. We observed a high degree of inter- and intra-GSC transcriptional heterogeneity that could not be fully explained by DNA somatic alterations. Instead, we found that GSCs mapped along a transcriptional gradient spanning two cellular states reminiscent of normal neural development and inflammatory wound response. Genome-wide CRISPR-Cas9 dropout screens independently recapitulated this observation, with each state characterized by unique essential genes. Further single-cell RNA sequencing of >56,000 malignant cells from primary tumors found that the majority organize along an orthogonal astrocyte maturation gradient yet retain expression of founder GSC transcriptional programs. We propose that glioblastomas grow out of a fundamental GSC-based neural wound response transcriptional program, which is a promising target for new therapy development.

Sphingosine-1-phosphate receptor 3 potentiates inflammatory programs in normal and leukemia stem cells to promote differentiation

Acute myeloid leukemia (AML) is a caricature of normal hematopoiesis, driven from leukemia stem cells (LSC) that share some hematopoietic stem cell (HSC) programs including responsiveness to inflammatory signaling. Although inflammation dysregulates mature myeloid cells and influences stemness programs and lineage determination in HSC by activating stress myelopoiesis, such roles in LSC are poorly understood. Here, we show that S1PR3, a receptor for the bioactive lipid sphingosine-1-phosphate, is a central regulator which drives myeloid differentiation and activates inflammatory programs in both HSC and LSC. S1PR3-mediated inflammatory signatures varied in a continuum from primitive to mature myeloid states across AML patient cohorts, each with distinct phenotypic and clinical properties. S1PR3 was high in LSC and blasts of mature myeloid samples with linkages to chemosensitivity, while S1PR3 activation in primitive samples promoted LSC differentiation leading to eradication. Our studies open new avenues for therapeutic target identification specific for each AML subset.