A Single-Cell Atlas of In Vivo Mammalian Chromatin Accessibility.

We applied a combinatorial indexing assay, sci-ATAC-seq, to profile genome-wide chromatin accessibility in ∼100,000 single cells from 13 adult mouse tissues. We identify 85 distinct patterns of chromatin accessibility, most of which can be assigned to cell types, and ∼400,000 differentially accessible elements. We use these data to link regulatory elements to their target genes, to define the transcription factor grammar specifying each cell type, and to discover in vivo correlates of heterogeneity in accessibility within cell types. We develop a technique for mapping single cell gene expression data to single-cell chromatin accessibility data, facilitating the comparison of atlases. By intersecting mouse chromatin accessibility with human genome-wide association summary statistics, we identify cell-type-specific enrichments of the heritability signal for hundreds of complex traits. These data define the in vivo landscape of the regulatory genome for common mammalian cell types at single-cell resolution.

A time-resolved, multi-symbol molecular recorder via sequential genome editing.

DNA is naturally well suited to serve as a digital medium for in vivo molecular recording. However, contemporary DNA-based memory devices are constrained in terms of the number of distinct 'symbols' that can be concurrently recorded and/or by a failure to capture the order in which events occur1. Here we describe DNA Typewriter, a general system for in vivo molecular recording that overcomes these and other limitations. For DNA Typewriter, the blank recording medium ('DNA Tape') consists of a tandem array of partial CRISPR-Cas9 target sites, with all but the first site truncated at their 5' ends and therefore inactive. Short insertional edits serve as symbols that record the identity of the prime editing guide RNA2 mediating the edit while also shifting the position of the 'type guide' by one unit along the DNA Tape, that is, sequential genome editing. In this proof of concept of DNA Typewriter, we demonstrate recording and decoding of thousands of symbols, complex event histories and short text messages; evaluate the performance of dozens of orthogonal tapes; and construct 'long tape' potentially capable of recording as many as 20 serial events. Finally, we leverage DNA Typewriter in conjunction with single-cell RNA-seq to reconstruct a monophyletic lineage of 3,257 cells and find that the Poisson-like accumulation of sequential edits to multicopy DNA tape can be maintained across at least 20 generations and 25 days of in vitro clonal expansion.

Multiplex profiling of developmental cis-regulatory elements with quantitative single-cell expression reporters.

The inability to scalably and precisely measure the activity of developmental cis-regulatory elements (CREs) in multicellular systems is a bottleneck in genomics. Here we develop a dual RNA cassette that decouples the detection and quantification tasks inherent to multiplex single-cell reporter assays. The resulting measurement of reporter expression is accurate over multiple orders of magnitude, with a precision approaching the limit set by Poisson counting noise. Together with RNA barcode stabilization via circularization, these scalable single-cell quantitative expression reporters provide high-contrast readouts, analogous to classic in situ assays but entirely from sequencing. Screening >200 regions of accessible chromatin in a multicellular in vitro model of early mammalian development, we identify 13 (8 previously uncharacterized) autonomous and cell-type-specific developmental CREs. We further demonstrate that chimeric CRE pairs generate cognate two-cell-type activity profiles and assess gain- and loss-of-function multicellular expression phenotypes from CRE variants with perturbed transcription factor binding sites. Single-cell quantitative expression reporters can be applied in developmental and multicellular systems to quantitatively characterize native, perturbed and synthetic CREs at scale, with high sensitivity and at single-cell resolution.

Editing TINF2 as a potential therapeutic approach to restore telomere length in dyskeratosis congenita.

Mutations in the TINF2 gene, encoding the shelterin protein TIN2, cause telomere shortening and the inherited bone marrow (BM) failure syndrome dyskeratosis congenita (DC). A lack of suitable model systems limits the mechanistic understanding of telomere shortening in the stem cells and thus hinders the development of treatment options for BM failure. Here, we endogenously introduced TIN2-DC mutations in human embryonic stem cells (hESCs) and human hematopoietic stem and progenitor cells (HSPCs) to dissect the disease mechanism and identify a gene-editing strategy that rescued the disease phenotypes. The hESCs with the T284R disease mutation exhibited the short telomere phenotype observed in DC patients. Yet, telomeres in mutant hESCs did not trigger DNA damage responses at telomeres or show exacerbated telomere shortening when differentiated into telomerase-negative cells. Disruption of the mutant TINF2 allele by introducing a frameshift mutation in exon 2 restored telomere length in stem cells and the replicative potential of differentiated cells. Similarly, we introduced TIN2-DC disease variants in human HSPCs to assess the changes in telomere length and proliferative capacity. Lastly, we showed that editing at exon 2 of TINF2 that restored telomere length in hESCs could be generated in TINF2-DC patient HSPCs. Our study demonstrates a simple genetic intervention that rescues the TIN2-DC disease phenotype in stem cells and provides a versatile platform to assess the efficacy of potential therapeutic approaches in vivo.

Genetic and molecular identification of three human TPP1 functions in telomerase action: recruitment, activation, and homeostasis set point regulation.

Telomere length homeostasis is essential for the long-term survival of stem cells, and its set point determines the proliferative capacity of differentiated cell lineages by restricting the reservoir of telomeric repeats. Knockdown and overexpression studies in human tumor cells showed that the shelterin subunit TPP1 recruits telomerase to telomeres through a region termed the TEL patch. However, these studies do not resolve whether the TPP1 TEL patch is the only mechanism for telomerase recruitment and whether telomerase regulation studied in tumor cells is representative of nontransformed cells such as stem cells. Using genome engineering of human embryonic stem cells, which have physiological telomere length homeostasis, we establish that the TPP1 TEL patch is genetically essential for telomere elongation and thus long-term cell viability. Furthermore, genetic bypass, protein fusion, and intragenic complementation assays define two distinct additional mechanisms of TPP1 involvement in telomerase action at telomeres. We demonstrate that TPP1 provides an essential step of telomerase activation as well as feedback regulation of telomerase by telomere length, which is necessary to determine the appropriate telomere length set point in human embryonic stem cells. These studies reveal and resolve multiple TPP1 roles in telomere elongation and stem cell telomere length homeostasis.

Retinoic acid induces human gastruloids with posterior embryo-like structures.

Gastruloids are a powerful in vitro model of early human development. However, although elongated and composed of all three germ layers, human gastruloids do not morphologically resemble post-implantation human embryos. Here we show that an early pulse of retinoic acid (RA), together with later Matrigel, robustly induces human gastruloids with posterior embryo-like morphological structures, including a neural tube flanked by segmented somites and diverse cell types, including neural crest, neural progenitors, renal progenitors and myocytes. Through in silico staging based on single-cell RNA sequencing, we find that human RA-gastruloids progress further than other human or mouse embryo models, aligning to E9.5 mouse and CS11 cynomolgus monkey embryos. We leverage chemical and genetic perturbations of RA-gastruloids to confirm that WNT and BMP signalling regulate somite formation and neural tube length in the human context, while transcription factors TBX6 and PAX3 underpin presomitic mesoderm and neural crest, respectively. Looking forward, RA-gastruloids are a robust, scalable model for decoding early human embryogenesis.

Induction and in silico staging of human gastruloids with neural tube, segmented somites & advanced cell types.

Embryonic organoids are emerging as powerful models for studying early mammalian development. For example, stem cell-derived 'gastruloids' form elongating structures containing all three germ layers1-4. However, although elongated, human gastruloids do not morphologically resemble post-implantation embryos. Here we show that a specific, discontinuous regimen of retinoic acid (RA) robustly induces human gastruloids with embryo-like morphological structures, including a neural tube and segmented somites. Single cell RNA-seq (sc-RNA-seq) further reveals that these human 'RA-gastruloids' contain more advanced cell types than conventional gastruloids, including neural crest cells, renal progenitor cells, skeletal muscle cells, and, rarely, neural progenitor cells. We apply a new approach to computationally stage human RA-gastruloids relative to somite-resolved mouse embryos, early human embryos and other gastruloid models, and find that the developmental stage of human RA-gastruloids is comparable to that of E9.5 mouse embryos, although some cell types show greater or lesser progression. We chemically perturb WNT and BMP signaling in human RA-gastruloids and find that these signaling pathways regulate somite patterning and neural tube length, respectively, while genetic perturbation of the transcription factors PAX3 and TBX6 markedly compromises the formation of neural crest and somites/renal cells, respectively. Human RA-gastruloids complement other embryonic organoids in serving as a simple, robust and screenable model for decoding early human embryogenesis.

Multiplex single-cell chemical genomics reveals the kinase dependence of the response to targeted therapy.

Chemical genetic screens are a powerful tool for exploring how cancer cells' response to drugs is shaped by their mutations, yet they lack a molecular view of the contribution of individual genes to the response to exposure. Here, we present sci-Plex-Gene-by-Environment (sci-Plex-GxE), a platform for combined single-cell genetic and chemical screening at scale. We highlight the advantages of large-scale, unbiased screening by defining the contribution of each of 522 human kinases to the response of glioblastoma to different drugs designed to abrogate signaling from the receptor tyrosine kinase pathway. In total, we probed 14,121 gene-by-environment combinations across 1,052,205 single-cell transcriptomes. We identify an expression signature characteristic of compensatory adaptive signaling regulated in a MEK/MAPK-dependent manner. Further analyses aimed at preventing adaptation revealed promising combination therapies, including dual MEK and CDC7/CDK9 or nuclear factor κB (NF-κB) inhibitors, as potent means of preventing transcriptional adaptation of glioblastoma to targeted therapy.

A non-genetic switch triggers alternative telomere lengthening and cellular immortalization in ATRX deficient cells.

Alternative Lengthening of Telomeres (ALT) is an aberrant DNA recombination pathway which grants replicative immortality to approximately 10% of all cancers. Despite this high prevalence of ALT in cancer, the mechanism and genetics by which cells activate this pathway remain incompletely understood. A major challenge in dissecting the events that initiate ALT is the extremely low frequency of ALT induction in human cell systems. Guided by the genetic lesions that have been associated with ALT from cancer sequencing studies, we genetically engineered primary human pluripotent stem cells to deterministically induce ALT upon differentiation. Using this genetically defined system, we demonstrate that disruption of the p53 and Rb pathways in combination with ATRX loss-of-function is sufficient to induce all hallmarks of ALT and results in functional immortalization in a cell type-specific manner. We further demonstrate that ALT can be induced in the presence of telomerase, is neither dependent on telomere shortening nor crisis, but is rather driven by continuous telomere instability triggered by the induction of differentiation in ATRX-deficient stem cells.

Multiplex single-cell chemical genomics reveals the kinase dependence of the response to targeted therapy.

Chemical genetic screens are a powerful tool for exploring how cancer cells' response to drugs is shaped by their mutations, yet they lack a molecular view of the contribution of individual genes to the response to exposure. Here, we present sci-Plex-Gene-by-Environment (sci-Plex-GxE), a platform for combined single-cell genetic and chemical screening at scale. We highlight the advantages of large-scale, unbiased screening by defining the contribution of each of 522 human kinases to the response of glioblastoma to different drugs designed to abrogate signaling from the receptor tyrosine kinase pathway. In total, we probed 14,121 gene-by-environment combinations across 1,052,205 single-cell transcriptomes. We identify an expression signature characteristic of compensatory adaptive signaling regulated in a MEK/MAPK-dependent manner. Further analyses aimed at preventing adaptation revealed promising combination therapies, including dual MEK and CDC7/CDK9 or NF-kB inhibitors, as potent means of preventing transcriptional adaptation of glioblastoma to targeted therapy.

Multiplex, single-cell CRISPRa screening for cell type specific regulatory elements.

CRISPR-based gene activation (CRISPRa) is a promising therapeutic approach for gene therapy, upregulating gene expression by targeting promoters or enhancers in a tissue/cell-type specific manner. Here, we describe an experimental framework that combines highly multiplexed perturbations with single-cell RNA sequencing (sc-RNA-seq) to identify cell-type-specific, CRISPRa-responsive cis- regulatory elements and the gene(s) they regulate. Random combinations of many gRNAs are introduced to each of many cells, which are then profiled and partitioned into test and control groups to test for effect(s) of CRISPRa perturbations of both enhancers and promoters on the expression of neighboring genes. Applying this method to candidate cis- regulatory elements in both K562 cells and iPSC-derived excitatory neurons, we identify gRNAs capable of specifically and potently upregulating target genes, including autism spectrum disorder (ASD) and neurodevelopmental disorder (NDD) risk genes. A consistent pattern is that the responsiveness of individual enhancers to CRISPRa is restricted by cell type, implying a dependency on either chromatin landscape and/or additional trans- acting factors for successful gene activation. The approach outlined here may facilitate large-scale screens for gRNAs that activate therapeutically relevant genes in a cell type-specific manner.

Telomere length set point regulation in human pluripotent stem cells critically depends on the shelterin protein TPP1.

Telomere maintenance is essential for the long-term proliferation of human pluripotent stem cells, while their telomere length set point determines the proliferative capacity of their differentiated progeny. The shelterin protein TPP1 is required for telomere stability and elongation, but its role in establishing a telomere length set point remains elusive. Here, we characterize the contribution of the shorter isoform of TPP1 (TPP1S) and the amino acid L104 outside the TEL patch, TPP1's telomerase interaction domain, to telomere length control. We demonstrate that cells deficient for TPP1S (TPP1S knockout [KO]), as well as the complete TPP1 KO cell lines, undergo telomere shortening. However, TPP1S KO cells are able to stabilize short telomeres, while TPP1 KO cells die. We compare these phenotypes with those of TPP1L104A/L104A mutant cells, which have short and stable telomeres similar to the TPP1S KO. In contrast to TPP1S KO cells, TPP1L104A/L104A cells respond to increased telomerase levels and maintain protected telomeres. However, TPP1L104A/L104A shows altered sensitivity to expression changes of shelterin proteins suggesting the mutation causes a defect in telomere length feedback regulation. Together this highlights TPP1L104A/L104A as the first shelterin mutant engineered at the endogenous locus of human stem cells with an altered telomere length set point.

Minimized human telomerase maintains telomeres and resolves endogenous roles of H/ACA proteins, TCAB1, and Cajal bodies.

We dissected the importance of human telomerase biogenesis and trafficking pathways for telomere maintenance. Biological stability of human telomerase RNA (hTR) relies on H/ACA proteins, but other eukaryotes use other RNP assembly pathways. To investigate additional rationale for human telomerase assembly as H/ACA RNP, we developed a minimized cellular hTR. Remarkably, with only binding sites for telomerase reverse transcriptase (TERT), minimized hTR assembled biologically active enzyme. TERT overexpression was required for cellular interaction with minimized hTR, indicating that H/ACA RNP assembly enhances endogenous hTR-TERT interaction. Telomere maintenance by minimized telomerase was unaffected by the elimination of the telomerase holoenzyme Cajal body chaperone TCAB1 or the Cajal body scaffold protein Coilin. Surprisingly, wild-type hTR also maintained and elongated telomeres in TCAB1 or Coilin knockout cells, with distinct changes in telomerase action. Overall, we elucidate trafficking requirements for telomerase biogenesis and function and expand mechanisms by which altered telomere maintenance engenders human disease.

Human intestinal tissue with adult stem cell properties derived from pluripotent stem cells.

Genetically engineered human pluripotent stem cells (hPSCs) have been proposed as a source for transplantation therapies and are rapidly becoming valuable tools for human disease modeling. However, many applications are limited due to the lack of robust differentiation paradigms that allow for the isolation of defined functional tissues. Here, using an endogenous LGR5-GFP reporter, we derived adult stem cells from hPSCs that gave rise to functional human intestinal tissue comprising all major cell types of the intestine. Histological and functional analyses revealed that such human organoid cultures could be derived with high purity and with a composition and morphology similar to those of cultures obtained from human biopsies. Importantly, hPSC-derived organoids responded to the canonical signaling pathways that control self-renewal and differentiation in the adult human intestinal stem cell compartment. This adult stem cell system provides a platform for studying human intestinal disease in vitro using genetically engineered hPSCs.

Restoring systemic GDF11 levels reverses age-related dysfunction in mouse skeletal muscle.

Parabiosis experiments indicate that impaired regeneration in aged mice is reversible by exposure to a young circulation, suggesting that young blood contains humoral "rejuvenating" factors that can restore regenerative function. Here, we demonstrate that the circulating protein growth differentiation factor 11 (GDF11) is a rejuvenating factor for skeletal muscle. Supplementation of systemic GDF11 levels, which normally decline with age, by heterochronic parabiosis or systemic delivery of recombinant protein, reversed functional impairments and restored genomic integrity in aged muscle stem cells (satellite cells). Increased GDF11 levels in aged mice also improved muscle structural and functional features and increased strength and endurance exercise capacity. These data indicate that GDF11 systemically regulates muscle aging and may be therapeutically useful for reversing age-related skeletal muscle and stem cell dysfunction.