The simultaneous sequencing of a single cell's genome and transcriptome offers a powerful means to dissect genetic variation and its effect on gene expression. Here we describe G&T-seq, a method for separating and sequencing genomic DNA and full-length mRNA from single cells. By applying G&T-seq to over 220 single cells from mice and humans, we discovered cellular properties that could not be inferred from DNA or RNA sequencing alone.

Most human pre-implantation embryos are mosaics of euploid and aneuploid cells. To determine the fate of aneuploid cells and the developmental potential of mosaic embryos, here we generate a mouse model of chromosome mosaicism. By treating embryos with a spindle assembly checkpoint inhibitor during the four- to eight-cell division, we efficiently generate aneuploid cells, resulting in embryo death during peri-implantation development. Live-embryo imaging and single-cell tracking in chimeric embryos, containing aneuploid and euploid cells, reveal that the fate of aneuploid cells depends on lineage: aneuploid cells in the fetal lineage are eliminated by apoptosis, whereas those in the placental lineage show severe proliferative defects. Overall, the proportion of aneuploid cells is progressively depleted from the blastocyst stage onwards. Finally, we show that mosaic embryos have full developmental potential, provided they contain sufficient euploid cells, a finding of significance for the assessment of embryo vitality in the clinic.

A major rate-limiting step in developing more effective immunotherapies for GBM is our inadequate understanding of the cellular complexity and the molecular heterogeneity of immune infiltrates in gliomas. Here, we report an integrated analysis of 201,986 human glioma, immune, and other stromal cells at the single cell level. In doing so, we discover extensive spatial and molecular heterogeneity in immune infiltrates. We identify molecular signatures for nine distinct myeloid cell subtypes, of which five are independent prognostic indicators of glioma patient survival. Furthermore, we identify S100A4 as a regulator of immune suppressive T and myeloid cells in GBM and demonstrate that deleting S100a4 in non-cancer cells is sufficient to reprogram the immune landscape and significantly improve survival. This study provides insights into spatial, molecular, and functional heterogeneity of glioma and glioma-associated immune cells and demonstrates the utility of this dataset for discovering therapeutic targets for this poorly immunogenic cancer.

TET2 somatic mutations occur in ~10% of DLBCLs but are of unknown significance. Herein we show that TET2 is required for the humoral immune response and is a DLBCL tumor suppressor. TET2 loss of function disrupts transit of B-cells through germinal centers (GC), causing GC hyperplasia, impaired class switch recombination, blockade of plasma cell differentiation and a pre-neoplastic phenotype. TET2 loss was linked to focal loss of enhancer hydroxymethylation and transcriptional repression of genes that mediate GC exit such as PRDM1. Notably, these enhancers and genes are also repressed in CREBBP-mutant DLBCLs. Accordingly, TET2 mutation in patients yields a CREBBP-mutant gene expression signature, CREBBP and TET2 mutations are generally mutually exclusive, and hydroxymethylation loss caused by TET2 deficiency impairs enhancer H3K27 acetylation. Hence TET2 plays a critical role in the GC reaction and its loss of function results in lymphomagenesis through failure to activate genes linked to GC exit signals.

Parallel sequencing of a single cell's genome and transcriptome provides a powerful tool for dissecting genetic variation and its relationship with gene expression. Here we present a detailed protocol for G&T-seq, a method for separation and parallel sequencing of genomic DNA and full-length polyA(+) mRNA from single cells. We provide step-by-step instructions for the isolation and lysis of single cells; the physical separation of polyA(+) mRNA from genomic DNA using a modified oligo-dT bead capture and the respective whole-transcriptome and whole-genome amplifications; and library preparation and sequence analyses of these amplification products. The method allows the detection of thousands of transcripts in parallel with the genetic variants captured by the DNA-seq data from the same single cell. G&T-seq differs from other currently available methods for parallel DNA and RNA sequencing from single cells, as it involves physical separation of the DNA and RNA and does not require bespoke microfluidics platforms. The process can be implemented manually or through automation. When performed manually, paired genome and transcriptome sequencing libraries from eight single cells can be produced in ∼3 d by researchers experienced in molecular laboratory work. For users with experience in the programming and operation of liquid-handling robots, paired DNA and RNA libraries from 96 single cells can be produced in the same time frame. Sequence analysis and integration of single-cell G&T-seq DNA and RNA data requires a high level of bioinformatics expertise and familiarity with a wide range of informatics tools.