Orit Rozenblatt-Rosen

Single-cell RNA-seq highlights intratumoral heterogeneity in primary glioblastoma

Cancer at single-cell resolution Single-cell sequencing can illuminate the genetic properties of brain cancers and reveal heterogeneity within a tumor. Patel et al. examined the genome sequence of single cells isolated from brain glioblastomas. The findings revealed shared chromosomal changes but also extensive transcription variation, including genes related to signaling, which represent potential therapeutic targets. The authors suggest that the variation in tumor cells reflects neural development and that such variation among cancer cells may prove to have clinical significance. Science, this issue p. 1396 Screening individual cancer cells within a brain tumor may help to guide treatment and predict prognosis. Human cancers are complex ecosystems composed of cells with distinct phenotypes, genotypes, and epigenetic states, but current models do not adequately reflect tumor composition in patients. We used single-cell RNA sequencing (RNA-seq) to profile 430 cells from five primary glioblastomas, which we found to be inherently variable in their expression of diverse transcriptional programs related to oncogenic signaling, proliferation, complement/immune response, and hypoxia. We also observed a continuum of stemness-related expression states that enabled us to identify putative regulators of stemness in vivo. Finally, we show that established glioblastoma subtype classifiers are variably expressed across individual cells within a tumor and demonstrate the potential prognostic implications of such intratumoral heterogeneity. Thus, we reveal previously unappreciated heterogeneity in diverse regulatory programs central to glioblastoma biology, prognosis, and therapy.

SOX2 is an amplified lineage-survival oncogene in lung and esophageal squamous cell carcinomas

Lineage-survival oncogenes are activated by somatic DNA alterations in cancers arising from the cell lineages in which these genes play a role in normal development. Here we show that a peak of genomic amplification on chromosome 3q26.33 found in squamous cell carcinomas (SCCs) of the lung and esophagus contains the transcription factor gene SOX2, which is mutated in hereditary human esophageal malformations, is necessary for normal esophageal squamous development, promotes differentiation and proliferation of basal tracheal cells and cooperates in induction of pluripotent stem cells. SOX2 expression is required for proliferation and anchorage-independent growth of lung and esophageal cell lines, as shown by RNA interference experiments. Furthermore, ectopic expression of SOX2 here cooperated with FOXE1 or FGFR2 to transform immortalized tracheobronchial epithelial cells. SOX2-driven tumors show expression of markers of both squamous differentiation and pluripotency. These characteristics identify SOX2 as a lineage-survival oncogene in lung and esophageal SCC.

Short interfering RNAs can induce unexpected and divergent changes in the levels of untargeted proteins in mammalian cells

RNA interference (RNAi) mediated by short interfering RNAs (siRNAs) is a widely used method to analyze gene function. To use RNAi knockdown accurately to infer gene function, it is essential to determine the specificity of siRNA-mediated RNAi. We have assessed the specificity of 10 different siRNAs corresponding to the MEN1 gene by examining the expression of two additional genes, TP53 (p53) and CDKN1A (p21), which are considered functionally unrelated to menin but are sensitive markers of cell state. MEN1 RNA and corresponding protein levels were all reduced after siRNA transfection of HeLa cells, although the degree of inhibition mediated by individual siRNAs varied. Unexpectedly, we observed dramatic and significant changes in protein levels of p53 and p21 that were unrelated to silencing of the target gene. The modulations in p53 and p21 levels were not abolished on titration of the siRNAs, and similar results were obtained in three other cell lines; in none of the cell lines tested did we see an effect on the protein levels of actin. These data suggest that siRNAs can induce nonspecific effects on protein levels that are siRNA sequence dependent but that these effects may be difficult to detect until genes central to a pivotal cellular response, such as p53 and p21, are studied. We find no evidence that activation of the double-stranded RNA-triggered IFN-associated antiviral pathways accounts for these effects, but we speculate that partial complementary sequence matches to off-target genes may result in a micro-RNA-like inhibition of translation.

Menin associates with a trithorax family histone methyltransferase complex and with the Hoxc8 locus

The cellular function of the menin tumor suppressor protein, product of the MEN1 gene mutated in familial multiple endocrine neoplasia type 1, has not been defined. We now show that menin is associated with a histone methyltransferase complex containing two trithorax family proteins, MLL2 and Ash2L, and other homologs of the yeast Set1 assembly. This menin-associated complex methylates histone H3 on lysine 4. A subset of tumor-derived menin mutants lacks the associated histone methyltransferase activity. In addition, menin is associated with RNA polymerase II whose large subunit carboxyl-terminal domain is phosphorylated on Ser 5. Men1 knockout embryos and cells show decreased expression of the homeobox genes Hoxc6 and Hoxc8. Chromatin immunoprecipitation experiments reveal that menin is bound to the Hoxc8 locus. These results suggest that menin activates the transcription of differentiation-regulating genes by covalent histone modification, and that this activity is related to tumor suppression by MEN1.

Dissecting the multicellular ecosystem of metastatic melanoma by single-cell RNA-seq

Single-cell expression profiles of melanoma Tumors harbor multiple cell types that are thought to play a role in the development of resistance to drug treatments. Tirosh et al. used single-cell sequencing to investigate the distribution of these differing genetic profiles within melanomas. Many cells harbored heterogeneous genetic programs that reflected two different states of genetic expression, one of which was linked to resistance development. Following drug treatment, the resistance-linked expression state was found at a much higher level. Furthermore, the environment of the melanoma cells affected their gene expression programs. Science, this issue p. 189 Melanoma cells show transcriptional heterogeneity. To explore the distinct genotypic and phenotypic states of melanoma tumors, we applied single-cell RNA sequencing (RNA-seq) to 4645 single cells isolated from 19 patients, profiling malignant, immune, stromal, and endothelial cells. Malignant cells within the same tumor displayed transcriptional heterogeneity associated with the cell cycle, spatial context, and a drug-resistance program. In particular, all tumors harbored malignant cells from two distinct transcriptional cell states, such that tumors characterized by high levels of the MITF transcription factor also contained cells with low MITF and elevated levels of the AXL kinase. Single-cell analyses suggested distinct tumor microenvironmental patterns, including cell-to-cell interactions. Analysis of tumor-infiltrating T cells revealed exhaustion programs, their connection to T cell activation and clonal expansion, and their variability across patients. Overall, we begin to unravel the cellular ecosystem of tumors and how single-cell genomics offers insights with implications for both targeted and immune therapies.

The Menin Tumor Suppressor Protein Is an Essential Oncogenic Cofactor for MLL-Associated Leukemogenesis

The Mixed-Lineage Leukemia (MLL) protein is a histone methyltransferase that is mutated in clinically and biologically distinctive subsets of acute leukemia. MLL normally associates with a cohort of highly conserved cofactors to form a macromolecular complex that includes menin, a product of the MEN1 tumor suppressor gene, which is mutated in heritable and sporadic endocrine tumors. We demonstrate here that oncogenic MLL fusion proteins retain an ability to stably associate with menin through a high-affinity, amino-terminal, conserved binding motif and that this interaction is required for the initiation of MLL-mediated leukemogenesis. Furthermore, menin is essential for maintenance of MLL-associated but not other oncogene induced myeloid transformation. Acute genetic ablation of menin reverses aberrant Hox gene expression mediated by MLL-menin promoter-associated complexes, and specifically abrogates the differentiation arrest and oncogenic properties of MLL-transformed leukemic blasts. These results demonstrate that a human oncoprotein is critically dependent on direct physical interaction with a tumor suppressor protein for its oncogenic activity, validate a potential target for molecular therapy, and suggest central roles for menin in altered epigenetic functions underlying the pathogenesis of hematopoietic cancers.

Interpreting cancer genomes using systematic host network perturbations by tumour virus proteins

Genotypic differences greatly influence susceptibility and resistance to disease. Understanding genotype–phenotype relationships requires that phenotypes be viewed as manifestations of network properties, rather than simply as the result of individual genomic variations. Genome sequencing efforts have identified numerous germline mutations, and large numbers of somatic genomic alterations, associated with a predisposition to cancer. However, it remains difficult to distinguish background, or ‘passenger’, cancer mutations from causal, or ‘driver’, mutations in these data sets. Human viruses intrinsically depend on their host cell during the course of infection and can elicit pathological phenotypes similar to those arising from mutations. Here we test the hypothesis that genomic variations and tumour viruses may cause cancer through related mechanisms, by systematically examining host interactome and transcriptome network perturbations caused by DNA tumour virus proteins. The resulting integrated viral perturbation data reflects rewiring of the host cell networks, and highlights pathways, such as Notch signalling and apoptosis, that go awry in cancer. We show that systematic analyses of host targets of viral proteins can identify cancer genes with a success rate on a par with their identification through functional genomics and large-scale cataloguing of tumour mutations. Together, these complementary approaches increase the specificity of cancer gene identification. Combining systems-level studies of pathogen-encoded gene products with genomic approaches will facilitate the prioritization of cancer-causing driver genes to advance the understanding of the genetic basis of human cancer.

The parafibromin tumor suppressor protein is part of a human Paf1 complex

ABSTRACT Parafibromin, the product of the HRPT2 (hyperparathyroidism-jaw tumor syndrome 2) tumor suppressor gene, is the human homologue of yeast Cdc73, part of the yeast RNA polymerase II/Paf1 complex known to be important for histone modification and connections to posttranscriptional events. By purifying cellular parafibromin and characterizing its associated proteins, we have identified a human counterpart to the yeast Paf1 complex including homologs of Leo1, Paf1, and Ctr9. Like the yeast complex, the parafibromin complex associates with the nonphosphorylated and Ser2 and Ser5 phosphorylated forms of the RNA polymerase II large subunit. Immunofluorescence experiments show that parafibromin is a nuclear protein. In addition, cotransfection data suggest that parafibromin can interact with a histone methyltransferase complex that methylates histone H3 on lysine 4. Some mutant forms of parafibromin lack association with hPaf1 complex members and with the histone methyltransferase complex, suggesting that disruption of these complexes may correlate with the oncogenic process.

Single-cell RNA-seq reveals new types of human blood dendritic cells, monocytes, and progenitors

Whats in a drop of blood? Blood contains many types of cells, including many immune system components. Immune cells used to be characterized by marker-based assays, but now classification relies on the genes that cells express. Villani et al. used deep sequencing at the single-cell level and unbiased clustering to define six dendritic cell and four monocyte populations. This refined analysis has identified, among others, a previously unknown dendritic cell population that potently activates T cells. Further cell culture revealed possible differentiation progenitors within the different cell populations. Science, this issue p. eaah4573 Discovery of additional immune cell subtypes will help identify functions and immune monitoring during disease. INTRODUCTION Dendritic cells (DCs) and monocytes consist of multiple specialized subtypes that play a central role in pathogen sensing, phagocytosis, and antigen presentation. However, their identities and interrelationships are not fully understood, as these populations have historically been defined by a combination of morphology, physical properties, localization, functions, developmental origins, and expression of a restricted set of surface markers. RATIONALE To overcome this inherently biased strategy for cell identification, we performed single-cell RNA sequencing of ~2400 cells isolated from healthy blood donors and enriched for HLA-DR+ lineage− cells. This single-cell profiling strategy and unbiased genomic classification, together with follow-up profiling and functional and phenotypic characterization of prospectively isolated subsets, led us to identify and validate six DC subtypes and four monocyte subtypes, and thus revise the taxonomy of these cells. RESULTS Our study reveals: 1) A new DC subset, representing 2 to 3% of the DC populations across all 10 donors tested, characterized by the expression of AXL, SIGLEC1, and SIGLEC6 antigens, named AS DCs. The AS DC population further divides into two populations captured in the traditionally defined plasmacytoid DC (pDC) and CD1C+ conventional DC (cDC) gates. This split is further reflected through AS DC gene expression signatures spanning a spectrum between cDC-like and pDC-like gene sets. Although AS DCs share properties with pDCs, they more potently activate T cells. This discovery led us to reclassify pDCs as the originally described “natural interferon-producing cells (IPCs)” with weaker T cell proliferation induction ability. 2) A new subdivision within the CD1C+ DC subset: one defined by a major histocompatibility complex class II–like gene set and one by a CD14+ monocyte–like prominent gene set. These CD1C+ DC subsets, which can be enriched by combining CD1C with CD32B, CD36, and CD163 antigens, can both potently induce T cell proliferation. 3) The existence of a circulating and dividing cDC progenitor giving rise to CD1C+ and CLEC9A+ DCs through in vitro differentiation assays. This blood precursor is defined by the expression of CD100+CD34int and observed at a frequency of ~0.02% of the LIN–HLA-DR+ fraction. 4) Two additional monocyte populations: one expressing classical monocyte genes and cytotoxic genes, and the other with unknown functions. 5) Evidence for a relationship between blastic plasmacytoid DC neoplasia (BPDCN) cells and healthy DCs. CONCLUSION Our revised taxonomy will enable more accurate functional and developmental analyses as well as immune monitoring in health and disease. The discovery of AS DCs within the traditionally defined pDC population explains many of the cDC properties previously assigned to pDCs, highlighting the need to revisit the definition of pDCs. Furthermore, the discovery of blood cDC progenitors represents a new therapeutic target readily accessible in the bloodstream for manipulation, as well as a new source for better in vitro DC generation. Although the current results focus on DCs and monocytes, a similar strategy can be applied to build a comprehensive human immune cell atlas. Establishing a human blood monocyte and dendritic cell atlas. We isolated ~2400 cells enriched from the healthy human blood lineage− HLA-DR+ compartment and subjected them to single-cell RNA sequencing. This strategy, together with follow-up profiling and functional and phenotypic characterization, led us to update the original cell classification to include six DCs, four monocyte subtypes, and one conventional DC progenitor. Dendritic cells (DCs) and monocytes play a central role in pathogen sensing, phagocytosis, and antigen presentation and consist of multiple specialized subtypes. However, their identities and interrelationships are not fully understood. Using unbiased single-cell RNA sequencing (RNA-seq) of ~2400 cells, we identified six human DCs and four monocyte subtypes in human blood. Our study reveals a new DC subset that shares properties with plasmacytoid DCs (pDCs) but potently activates T cells, thus redefining pDCs; a new subdivision within the CD1C+ subset of DCs; the relationship between blastic plasmacytoid DC neoplasia cells and healthy DCs; and circulating progenitor of conventional DCs (cDCs). Our revised taxonomy will enable more accurate functional and developmental analyses as well as immune monitoring in health and disease.

Single-cell RNA-seq supports a developmental hierarchy in human oligodendroglioma

Although human tumours are shaped by the genetic evolution of cancer cells, evidence also suggests that they display hierarchies related to developmental pathways and epigenetic programs in which cancer stem cells (CSCs) can drive tumour growth and give rise to differentiated progeny. Yet, unbiased evidence for CSCs in solid human malignancies remains elusive. Here we profile 4,347 single cells from six IDH1 or IDH2 mutant human oligodendrogliomas by RNA sequencing (RNA-seq) and reconstruct their developmental programs from genome-wide expression signatures. We infer that most cancer cells are differentiated along two specialized glial programs, whereas a rare subpopulation of cells is undifferentiated and associated with a neural stem cell expression program. Cells with expression signatures for proliferation are highly enriched in this rare subpopulation, consistent with a model in which CSCs are primarily responsible for fuelling the growth of oligodendroglioma in humans. Analysis of copy number variation (CNV) shows that distinct CNV sub-clones within tumours display similar cellular hierarchies, suggesting that the architecture of oligodendroglioma is primarily dictated by developmental programs. Subclonal point mutation analysis supports a similar model, although a full phylogenetic tree would be required to definitively determine the effect of genetic evolution on the inferred hierarchies. Our single-cell analyses provide insight into the cellular architecture of oligodendrogliomas at single-cell resolution and support the cancer stem cell model, with substantial implications for disease management.