Resolving medulloblastoma cellular architecture by single-cell genomics.
2019; 572 (7767): 74–79
LC3-Associated Phagocytosis in Myeloid Cells Promotes Tumor Immune Tolerance.
2018; 175 (2): 429–41.e16
Medulloblastoma is a malignant childhood cerebellar tumour type that comprises distinct molecular subgroups. Whereas genomic characteristics of these subgroups are well defined, the extent to which cellular diversity underlies their divergent biology and clinical behaviour remains largely unexplored. Here we used single-cell transcriptomics to investigate intra- and intertumoral heterogeneity in 25 medulloblastomas spanning all molecular subgroups. WNT, SHH and Group 3 tumours comprised subgroup-specific undifferentiated and differentiated neuronal-like malignant populations, whereas Group 4 tumours consisted exclusively of differentiated neuronal-like neoplastic cells. SHH tumours closely resembled granule neurons of varying differentiation states that correlated with patient age. Group 3 and Group 4 tumours exhibited a developmental trajectory from primitive progenitor-like to more mature neuronal-like cells, the relative proportions of which distinguished these subgroups. Cross-species transcriptomics defined distinct glutamatergic populations as putative cells-of-origin for SHH and Group 4 subtypes. Collectively, these data provide insights into the cellular and developmental states underlying subtype-specific medulloblastoma biology.
View details for DOI 10.1038/s41586-019-1434-6
View details for PubMedID 31341285
Pan-cancer genome and transcriptome analyses of 1,699 paediatric leukaemias and solid tumours.
2018; 555 (7696): 371–76
Targeting autophagy in cancer cells and in the tumor microenvironment are current goals of cancer therapy. However, components of canonical autophagy play roles in other biological processes, adding complexity to this goal. One such alternative function of autophagy proteins is LC3-associated phagocytosis (LAP), which functions in phagosome maturation and subsequent signaling events. Here, we show that impairment of LAP in the myeloid compartment, rather than canonical autophagy, induces control of tumor growth by tumor-associated macrophages (TAM) upon phagocytosis of dying tumor cells. Single-cell RNA sequencing (RNA-seq) analysis revealed that defects in LAP induce pro-inflammatory gene expression and trigger STING-mediated type I interferon responses in TAM. We found that the anti-tumor effects of LAP impairment require tumor-infiltrating T cells, dependent upon STING and the type I interferon response. Therefore, autophagy proteins in the myeloid cells of the tumor microenvironment contribute to immune suppression of T lymphocytes by effecting LAP.
View details for DOI 10.1016/j.cell.2018.08.061
View details for PubMedID 30245008
View details for PubMedCentralID PMC6201245
A Single-Cell Transcriptional Atlas of the Developing Murine Cerebellum.
Current biology : CB
2018; 28 (18): 2910–20.e2
Analysis of molecular aberrations across multiple cancer types, known as pan-cancer analysis, identifies commonalities and differences in key biological processes that are dysregulated in cancer cells from diverse lineages. Pan-cancer analyses have been performed for adult but not paediatric cancers, which commonly occur in developing mesodermic rather than adult epithelial tissues. Here we present a pan-cancer study of somatic alterations, including single nucleotide variants, small insertions or deletions, structural variations, copy number alterations, gene fusions and internal tandem duplications in 1,699 paediatric leukaemias and solid tumours across six histotypes, with whole-genome, whole-exome and transcriptome sequencing data processed under a uniform analytical framework. We report 142 driver genes in paediatric cancers, of which only 45% match those found in adult pan-cancer studies; copy number alterations and structural variants constituted the majority (62%) of events. Eleven genome-wide mutational signatures were identified, including one attributed to ultraviolet-light exposure in eight aneuploid leukaemias. Transcription of the mutant allele was detectable for 34% of protein-coding mutations, and 20% exhibited allele-specific expression. These data provide a comprehensive genomic architecture for paediatric cancers and emphasize the need for paediatric cancer-specific development of precision therapies.
View details for DOI 10.1038/nature25795
View details for PubMedID 29489755
View details for PubMedCentralID PMC5854542
Dissecting the clonal origins of childhood acute lymphoblastic leukemia by single-cell genomics.
Proceedings of the National Academy of Sciences of the United States of America
2014; 111 (50): 17947-17952
The cerebellum develops from a restricted number of cell types that precisely organize to form the circuitry that controls sensory-motor coordination and some higher-order cognitive processes. To acquire an enhanced understanding of the molecular processes that mediate cerebellar development, we performed single-cell RNA-sequencing of 39,245 murine cerebellar cells at twelve critical developmental time points. Using recognized lineage markers, we confirmed that the single-cell data accurately recapitulate cerebellar development. We then followed distinct populations from emergence through migration and differentiation, and determined the associated transcriptional cascades. After identifying key lineage commitment decisions, focused analyses uncovered waves of transcription factor expression at those branching points. Finally, we created Cell Seek, a flexible online interface that facilitates exploration of the dataset. Our study provides a transcriptional summarization of cerebellar development at single-cell resolution that will serve as a valuable resource for future investigations of cerebellar development, neurobiology, and disease.
View details for DOI 10.1016/j.cub.2018.07.062
View details for PubMedID 30220501
Circular RNAs Are the Predominant Transcript Isoform from Hundreds of Human Genes in Diverse Cell Types
2012; 7 (2)
Many cancers have substantial genomic heterogeneity within a given tumor, and to fully understand that diversity requires the ability to perform single cell analysis. We performed targeted sequencing of a panel of single nucleotide variants (SNVs), deletions, and IgH sequences in 1,479 single tumor cells from six acute lymphoblastic leukemia (ALL) patients. By accurately segregating groups of cooccurring mutations into distinct clonal populations, we identified codominant clones in the majority of patients. Evaluation of intraclonal mutation patterns identified clone-specific punctuated cytosine mutagenesis events, showed that most structural variants are acquired before SNVs, determined that KRAS mutations occur late in disease development but are not sufficient for clonal dominance, and identified clones within the same patient that are arrested at varied stages in B-cell development. Taken together, these data order the sequence of genetic events that underlie childhood ALL and provide a framework for understanding the development of the disease at single-cell resolution.
View details for DOI 10.1073/pnas.1420822111
View details for PubMedID 25425670
Single-cell RNA sequencing reveals the impact of chromosomal instability on glioblastoma cancer stem cells.
BMC medical genomics
2019; 12 (1): 79
Most human pre-mRNAs are spliced into linear molecules that retain the exon order defined by the genomic sequence. By deep sequencing of RNA from a variety of normal and malignant human cells, we found RNA transcripts from many human genes in which the exons were arranged in a non-canonical order. Statistical estimates and biochemical assays provided strong evidence that a substantial fraction of the spliced transcripts from hundreds of genes are circular RNAs. Our results suggest that a non-canonical mode of RNA splicing, resulting in a circular RNA isoform, is a general feature of the gene expression program in human cells.
View details for DOI 10.1371/journal.pone.0030733
View details for Web of Science ID 000301977500016
View details for PubMedID 22319583
View details for PubMedCentralID PMC3270023
Sequencing the Genomes of Single Cells.
Methods in molecular biology (Clifton, N.J.)
2019; 1979: 227–34
Intra-tumor heterogeneity stems from genetic, epigenetic, functional, and environmental differences among tumor cells. A major source of genetic heterogeneity comes from DNA sequence differences and/or whole chromosome and focal copy number variations (CNVs). Whole chromosome CNVs are caused by chromosomal instability (CIN) that is defined by a persistently high rate of chromosome mis-segregation. Accordingly, CIN causes constantly changing karyotypes that result in extensive cell-to-cell genetic heterogeneity. How the genetic heterogeneity caused by CIN influences gene expression in individual cells remains unknown.We performed single-cell RNA sequencing on a chromosomally unstable glioblastoma cancer stem cell (CSC) line and a control normal, diploid neural stem cell (NSC) line to investigate the impact of CNV due to CIN on gene expression. From the gene expression data, we computationally inferred large-scale CNVs in single cells. Also, we performed copy number adjusted differential gene expression analysis between NSCs and glioblastoma CSCs to identify copy number dependent and independent differentially expressed genes.Here, we demonstrate that gene expression across large genomic regions scales proportionally to whole chromosome copy number in chromosomally unstable CSCs. Also, we show that the differential expression of most genes between normal NSCs and glioblastoma CSCs is largely accounted for by copy number alterations. However, we identify 269 genes whose differential expression in glioblastoma CSCs relative to normal NSCs is independent of copy number. Moreover, a gene signature derived from the subset of genes that are differential expressed independent of copy number in glioblastoma CSCs correlates with tumor grade and is prognostic for patient survival.These results demonstrate that CIN is directly responsible for gene expression changes and contributes to both genetic and transcriptional heterogeneity among glioblastoma CSCs. These results also demonstrate that the expression of some genes is buffered against changes in copy number, thus preserving some consistency in gene expression levels from cell-to-cell despite the continuous change in karyotype driven by CIN. Importantly, a gene signature derived from the subset of genes whose expression is buffered against copy number alterations correlates with tumor grade and is prognostic for patient survival that could facilitate patient diagnosis and treatment.
View details for DOI 10.1186/s12920-019-0532-5
View details for PubMedID 31151460
View details for PubMedCentralID PMC6545015
High-resolution transcriptional dissection of in vivo Atoh1-mediated hair cell conversion in mature cochleae identifies Isl1 as a co-reprogramming factor.
2018; 14 (7): e1007552
Single-cell genome sequencing can detect low-frequency genetic alterations present in complex tissues. However, the experimental procedures are technically challenging. This includes dissociation of the tissue, isolation of single cells, whole-genome amplification, sequencing library preparation, and an optional target enrichment. Here we describe how to perform each of these processes to obtain high-quality single-cell genome sequencing data.
View details for DOI 10.1007/978-1-4939-9240-9_14
View details for PubMedID 31028641
Murine hematopoietic stem cell activity is derived from pre-circulation embryos but not yolk sacs.
2018; 9 (1): 5405
In vivo direct conversion of differentiated cells holds promise for regenerative medicine; however, improving the conversion efficiency and producing functional target cells remain challenging. Ectopic Atoh1 expression in non-sensory supporting cells (SCs) in mouse cochleae induces their partial conversion to hair cells (HCs) at low efficiency. Here, we performed single-cell RNA sequencing of whole mouse sensory epithelia harvested at multiple time points after conditional overexpression of Atoh1. Pseudotemporal ordering revealed that converted HCs (cHCs) are present along a conversion continuum that correlates with both endogenous and exogenous Atoh1 expression. Bulk sequencing of isolated cell populations and single-cell qPCR confirmed 51 transcription factors, including Isl1, are differentially expressed among cHCs, SCs and HCs. In transgenic mice, co-overexpression of Atoh1 and Isl1 enhanced the HC conversion efficiency. Together, our study shows how high-resolution transcriptional profiling of direct cell conversion can identify co-reprogramming factors required for efficient conversion.
View details for DOI 10.1371/journal.pgen.1007552
View details for PubMedID 30063705
View details for PubMedCentralID PMC6086484
Early somatic mosaicism is a rare cause of long-QT syndrome
PROCEEDINGS OF THE NATIONAL ACADEMY OF SCIENCES OF THE UNITED STATES OF AMERICA
2016; 113 (41): 11555-11560
The embryonic site of definitive hematopoietic stem cell (dHSC) origination has been debated for decades. Although an intra-embryonic origin is well supported, the yolk sac (YS) contribution to adult hematopoiesis remains controversial. The same developmental origin makes it difficult to identify specific markers that discern between an intraembryonic versus YS-origin using a lineage trace approach. Additionally, the highly migratory nature of blood cells and the inability of pre-circulatory embryonic cells (i.e., 5-7 somite pairs (sp)) to robustly engraft in transplantation, even after culture, has precluded scientists from properly answering these questions. Here we report robust, multi-lineage and serially transplantable dHSC activity from cultured 2-7sp murine embryonic explants (Em-Ex). dHSC are undetectable in 2-7sp YS explants. Additionally, the engraftment from Em-Ex is confined to an emerging CD31+CD45+c-Kit+CD41- population. In sum, our work supports a model in which the embryo, not the YS, is the major source of lifelong definitive hematopoiesis.
View details for DOI 10.1038/s41467-018-07769-8
View details for PubMedID 30573729
View details for PubMedCentralID PMC6302089
Single-cell genome sequencing: current state of the science
NATURE REVIEWS GENETICS
2016; 17 (3): 175-188
Somatic mosaicism, the occurrence and propagation of genetic variation in cell lineages after fertilization, is increasingly recognized to play a causal role in a variety of human diseases. We investigated the case of life-threatening arrhythmia in a 10-day-old infant with long QT syndrome (LQTS). Rapid genome sequencing suggested a variant in the sodium channel NaV1.5 encoded by SCN5A, NM_000335:c.5284G > T predicting p.(V1762L), but read depth was insufficient to be diagnostic. Exome sequencing of the trio confirmed read ratios inconsistent with Mendelian inheritance only in the proband. Genotyping of single circulating leukocytes demonstrated the mutation in the genomes of 8% of patient cells, and RNA sequencing of cardiac tissue from the infant confirmed the expression of the mutant allele at mosaic ratios. Heterologous expression of the mutant channel revealed significantly delayed sodium current with a dominant negative effect. To investigate the mechanism by which mosaicism might cause arrhythmia, we built a finite element simulation model incorporating Purkinje fiber activation. This model confirmed the pathogenic consequences of cardiac cellular mosaicism and, under the presenting conditions of this case, recapitulated 2:1 AV block and arrhythmia. To investigate the extent to which mosaicism might explain undiagnosed arrhythmia, we studied 7,500 affected probands undergoing commercial gene-panel testing. Four individuals with pathogenic variants arising from early somatic mutation events were found. Here we establish cardiac mosaicism as a causal mechanism for LQTS and present methods by which the general phenomenon, likely to be relevant for all genetic diseases, can be detected through single-cell analysis and next-generation sequencing.
View details for DOI 10.1073/pnas.1607187113
View details for PubMedID 27681629
Dynamic ASXL1 Exon Skipping and Alternative Circular Splicing in Single Human Cells.
2016; 11 (10)
The field of single-cell genomics is advancing rapidly and is generating many new insights into complex biological systems, ranging from the diversity of microbial ecosystems to the genomics of human cancer. In this Review, we provide an overview of the current state of the field of single-cell genome sequencing. First, we focus on the technical challenges of making measurements that start from a single molecule of DNA, and then explore how some of these recent methodological advancements have enabled the discovery of unexpected new biology. Areas highlighted include the application of single-cell genomics to interrogate microbial dark matter and to evaluate the pathogenic roles of genetic mosaicism in multicellular organisms, with a focus on cancer. We then attempt to predict advances we expect to see in the next few years.
View details for DOI 10.1038/nrg.2015.16
View details for Web of Science ID 000371739200012
View details for PubMedID 26806412
Noninvasive in vivo monitoring of tissue-specific global gene expression in humans.
Proceedings of the National Academy of Sciences of the United States of America
2014; 111 (20): 7361-7366
Circular RNAs comprise a poorly understood new class of noncoding RNA. In this study, we used a combination of targeted deletion, high-resolution splicing detection, and single-cell sequencing to deeply probe ASXL1 circular splicing. We found that efficient circular splicing required the canonical transcriptional start site and inverted AluSx elements. Sequencing-based interrogation of isoforms after ASXL1 overexpression identified promiscuous linear splicing between all exons, with the two most abundant non-canonical linear products skipping the exons that produced the circular isoforms. Single-cell sequencing revealed a strong preference for either the linear or circular ASXL1 isoforms in each cell, and found the predominant exon skipping product is frequently co-expressed with its reciprocal circular isoform. Finally, absolute quantification of ASXL1 isoforms confirmed our findings and suggests that standard methods overestimate circRNA abundance. Taken together, these data reveal a dynamic new view of circRNA genesis, providing additional framework for studying their roles in cellular biology.
View details for DOI 10.1371/journal.pone.0164085
View details for PubMedID 27736885
View details for PubMedCentralID PMC5063410
A quantitative comparison of single-cell whole genome amplification methods.
2014; 9 (8)
Circulating cell-free RNA in the blood provides a potential window into the health, phenotype, and developmental programs of a variety of human organs. We used high-throughput methods of RNA analysis such as microarrays and next-generation sequencing to characterize the global landscape circulating RNA in a cohort of human subjects. By focusing on genes whose expression is highly specific to certain tissues, we were able to identify the relative contributions of these tissues to circulating RNA and to monitor changes in tissue development and health. As one application of this approach, we performed a longitudinal study on pregnant women and analyzed their combined cell-free RNA transcriptomes across all three trimesters of pregnancy and after delivery. In addition to the analysis of mRNA, we observed and characterized noncoding species such as long noncoding RNA and circular RNA transcripts whose presence had not been previously observed in human plasma. We demonstrate that it is possible to track specific longitudinal phenotypic changes in both the mother and the fetus and that it is possible to directly measure transcripts from a variety of fetal tissues in the maternal blood sample. We also studied the role of neuron-specific transcripts in the blood of healthy adults and those suffering from the neurodegenerative disorder Alzheimer's disease and showed that disease specific neural transcripts are present at increased levels in the blood of affected individuals. Characterization of the cell-free transcriptome in its entirety may thus provide broad insights into human health and development without the need for invasive tissue sampling.
View details for DOI 10.1073/pnas.1405528111
View details for PubMedID 24799715
Massive evolution of the immunoglobulin heavy chain locus in children with B precursor acute lymphoblastic leukemia
2012; 120 (22): 4407-4417
Single-cell sequencing is emerging as an important tool for studies of genomic heterogeneity. Whole genome amplification (WGA) is a key step in single-cell sequencing workflows and a multitude of methods have been introduced. Here, we compare three state-of-the-art methods on both bulk and single-cell samples of E. coli DNA: Multiple Displacement Amplification (MDA), Multiple Annealing and Looping Based Amplification Cycles (MALBAC), and the PicoPLEX single-cell WGA kit (NEB-WGA). We considered the effects of reaction gain on coverage uniformity, error rates and the level of background contamination. We compared the suitability of the different WGA methods for the detection of copy-number variations, for the detection of single-nucleotide polymorphisms and for de-novo genome assembly. No single method performed best across all criteria and significant differences in characteristics were observed; the choice of which amplifier to use will depend strongly on the details of the type of question being asked in any given experiment.
View details for DOI 10.1371/journal.pone.0105585
View details for PubMedID 25136831
View details for PubMedCentralID PMC4138190
Towards molecular medicine - A case for a biological periodic table
AMERICAN JOURNAL OF PHARMACOGENOMICS
2005; 5 (4): 207-211
The ability to distinguish clonal B-cell populations based on the sequence of their rearranged immunoglobulin heavy chain (IgH) locus is an important tool for diagnosing B-cell neoplasms and monitoring treatment response. Leukemic precursor B cells may continue to undergo recombination of the IgH gene after malignant transformation; however, the magnitude of evolution at the IgH locus is currently unknown. We used next-generation sequencing to characterize the repertoire of IgH sequences in diagnostic samples of 51 children with B precursor acute lymphoblastic leukemia (B-ALL). We identified clonal IgH rearrangements in 43 of 51 (84%) cases and found that the number of evolved IgH sequences per patient ranged dramatically from 0 to 4024. We demonstrate that the evolved IgH sequences are not the result of amplification artifacts and are unique to leukemic precursor B cells. In addition, the evolution often follows an allelic exclusion pattern, where only 1 of 2 rearranged IgH loci exhibit ongoing recombination. Thus, precursor B-cell leukemias maintain evolution at the IgH locus at levels that were previously underappreciated. This finding sheds light on the mechanisms associated with leukemic clonal evolution and may fundamentally change approaches for monitoring minimal residual disease burden.
View details for DOI 10.1182/blood-2012-05-429811
View details for Web of Science ID 000313111300023
View details for PubMedID 22932801
View details for PubMedCentralID PMC3507147
The recently amplified pace of development in the technologies to study both normal and aberrant cellular physiology has allowed for a transition from the traditional reductionist approaches to global interrogations of human biology. This transformation has created the anticipation that we will soon more effectively treat or contain most types of diseases through a 'systems-based' approach to understanding and correcting the underlying etiology of these processes. However, to accomplish these goals, we must first have a more comprehensive understanding of all the elements involved in human cellular physiology, as well as why and how they interact. With the vast number of biological components that have and are being discovered, creating methods with modern computational techniques to better organize biological elements is the next requisite step in this process. This article aims to articulate the importance of the organization of chemical elements into a periodic table had on the conversion of chemistry into a quantitative, translatable science, as well as how we can apply the lessons learned in that transition to the current transformation taking place in biology.
View details for Web of Science ID 000231842000001
View details for PubMedID 16078857