Mechanism-guided engineering of a minimal biological particle for genome editing.

The widespread application of genome editing to treat or even cure disease requires the delivery of genome editors into the nucleus of target cells. Enveloped Delivery Vehicles (EDVs) are engineered virally-derived particles capable of packaging and delivering CRISPR-Cas9 ribonucleoproteins (RNPs). However, the presence of lentiviral genome encapsulation and replication components in EDVs has obscured the underlying delivery mechanism and precluded particle optimization. Here we show that Cas9 RNP nuclear delivery is independent of the native lentiviral capsid structure. Instead, EDV-mediated genome editing activity corresponds directly to the number of nuclear localization sequences on the Cas9 enzyme. EDV structural analysis using cryo-electron tomography and small molecule inhibitors guided the removal of ∼80% of viral residues, creating a minimal EDV (miniEDV) that retains full RNP delivery capability. MiniEDVs are 25% smaller yet package equivalent amounts of Cas9 RNPs relative to the original EDVs, and demonstrated increased editing in cell lines and therapeutically-relevant primary human T cells. These results show that virally-derived particles can be streamlined to create efficacious genome editing delivery vehicles that could simplify production and manufacturing. SIGNIFICANCE STATEMENT: Our results highlight the importance of understanding how virally-derived particles function to eliminate unnecessary viral proteins and create more efficacious and easier-to-produce delivery vehicles for therapeutic genome editing.

Transcript-specific enrichment enables profiling rare cell states via scRNA-seq.

Single-cell genomics technologies have accelerated our understanding of cell-state heterogeneity in diverse contexts. Although single-cell RNA sequencing (scRNA-seq) identifies many rare populations of interest that express specific marker transcript combinations, traditional flow sorting limits our ability to enrich these populations for further profiling, including requiring cell surface markers with high-fidelity antibodies. Additionally, many single-cell studies require the isolation of nuclei from tissue, eliminating the ability to enrich learned rare cell states based on extranuclear protein markers. To address these limitations, we describe Programmable Enrichment via RNA Flow-FISH by sequencing (PERFF-seq), a scalable assay that enables scRNA-seq profiling of subpopulations from complex cellular mixtures defined by the presence or absence of specific RNA transcripts. Across immune populations (n = 141,227 cells) and fresh-frozen and formalin-fixed paraffin-embedded brain tissue (n = 29,522 nuclei), we demonstrate the sorting logic that can be used to enrich for cell populations via RNA-based cytometry followed by high-throughput scRNA-seq. Our approach provides a rational, programmable method for studying rare populations identified by one or more marker transcripts.

Latent human herpesvirus 6 is reactivated in CAR T cells.

Cell therapies have yielded durable clinical benefits for patients with cancer, but the risks associated with the development of therapies from manipulated human cells are understudied. For example, we lack a comprehensive understanding of the mechanisms of toxicities observed in patients receiving T cell therapies, including recent reports of encephalitis caused by reactivation of human herpesvirus 6 (HHV-6)1. Here, through petabase-scale viral genomics mining, we examine the landscape of human latent viral reactivation and demonstrate that HHV-6B can become reactivated in cultures of human CD4+ T cells. Using single-cell sequencing, we identify a rare population of HHV-6 'super-expressors' (about 1 in 300-10,000 cells) that possess high viral transcriptional activity, among research-grade allogeneic chimeric antigen receptor (CAR) T cells. By analysing single-cell sequencing data from patients receiving cell therapy products that are approved by the US Food and Drug Administration2 or are in clinical studies3-5, we identify the presence of HHV-6-super-expressor CAR T cells in patients in vivo. Together, the findings of our study demonstrate the utility of comprehensive genomics analyses in implicating cell therapy products as a potential source contributing to the lytic HHV-6 infection that has been reported in clinical trials1,6-8 and may influence the design and production of autologous and allogeneic cell therapies.

Spatiotemporal transcriptomic maps of whole mouse embryos at the onset of organogenesis.

Spatiotemporal orchestration of gene expression is required for proper embryonic development. The use of single-cell technologies has begun to provide improved resolution of early regulatory dynamics, including detailed molecular definitions of most cell states during mouse embryogenesis. Here we used Slide-seq to build spatial transcriptomic maps of complete embryonic day (E) 8.5 and E9.0, and partial E9.5 embryos. To support their utility, we developed sc3D, a tool for reconstructing and exploring three-dimensional 'virtual embryos', which enables the quantitative investigation of regionalized gene expression patterns. Our measurements along the main embryonic axes of the developing neural tube revealed several previously unannotated genes with distinct spatial patterns. We also characterized the conflicting transcriptional identity of 'ectopic' neural tubes that emerge in Tbx6 mutant embryos. Taken together, we present an experimental and computational framework for the spatiotemporal investigation of whole embryonic structures and mutant phenotypes.

High-resolution Slide-seqV2 spatial transcriptomics enables discovery of disease-specific cell neighborhoods and pathways.

High-resolution spatial transcriptomics enables mapping of RNA expression directly from intact tissue sections; however, its utility for the elucidation of disease processes and therapeutically actionable pathways remains unexplored. We applied Slide-seqV2 to mouse and human kidneys, in healthy and distinct disease paradigms. First, we established the feasibility of Slide-seqV2 in tissue from nine distinct human kidneys, which revealed a cell neighborhood centered around a population of LYVE1+ macrophages. Second, in a mouse model of diabetic kidney disease, we detected changes in the cellular organization of the spatially restricted kidney filter and blood-flow-regulating apparatus. Third, in a mouse model of a toxic proteinopathy, we identified previously unknown, disease-specific cell neighborhoods centered around macrophages. In a spatially restricted subpopulation of epithelial cells, we discovered perturbations in 77 genes associated with the unfolded protein response. Our studies illustrate and experimentally validate the utility of Slide-seqV2 for the discovery of disease-specific cell neighborhoods.

Molecular logic of cellular diversification in the mouse cerebral cortex.

The mammalian cerebral cortex has an unparalleled diversity of cell types, which are generated during development through a series of temporally orchestrated events that are under tight evolutionary constraint and are critical for proper cortical assembly and function1,2. However, the molecular logic that governs the establishment and organization of cortical cell types remains unknown, largely due to the large number of cell classes that undergo dynamic cell-state transitions over extended developmental timelines. Here we generate a comprehensive atlas of the developing mouse neocortex, using single-cell RNA sequencing and single-cell assay for transposase-accessible chromatin using sequencing. We sampled the neocortex every day throughout embryonic corticogenesis and at early postnatal ages, and complemented the sequencing data with a spatial transcriptomics time course. We computationally reconstruct developmental trajectories across the diversity of cortical cell classes, and infer their spatial organization and the gene regulatory programs that accompany their lineage bifurcation decisions and differentiation trajectories. Finally, we demonstrate how this developmental map pinpoints the origin of lineage-specific developmental abnormalities that are linked to aberrant corticogenesis in mutant mice. The data provide a global picture of the regulatory mechanisms that govern cellular diversification in the neocortex.

Highly sensitive spatial transcriptomics at near-cellular resolution with Slide-seqV2.

Measurement of the location of molecules in tissues is essential for understanding tissue formation and function. Previously, we developed Slide-seq, a technology that enables transcriptome-wide detection of RNAs with a spatial resolution of 10 µm. Here we report Slide-seqV2, which combines improvements in library generation, bead synthesis and array indexing to reach an RNA capture efficiency ~50% that of single-cell RNA-seq data (~10-fold greater than Slide-seq), approaching the detection efficiency of droplet-based single-cell RNA-seq techniques. First, we leverage the detection efficiency of Slide-seqV2 to identify dendritically localized mRNAs in neurons of the mouse hippocampus. Second, we integrate the spatial information of Slide-seqV2 data with single-cell trajectory analysis tools to characterize the spatiotemporal development of the mouse neocortex, identifying underlying genetic programs that were poorly sampled with Slide-seq. The combination of near-cellular resolution and high transcript detection efficiency makes Slide-seqV2 useful across many experimental contexts.

Slide-seq: A scalable technology for measuring genome-wide expression at high spatial resolution.

Spatial positions of cells in tissues strongly influence function, yet a high-throughput, genome-wide readout of gene expression with cellular resolution is lacking. We developed Slide-seq, a method for transferring RNA from tissue sections onto a surface covered in DNA-barcoded beads with known positions, allowing the locations of the RNA to be inferred by sequencing. Using Slide-seq, we localized cell types identified by single-cell RNA sequencing datasets within the cerebellum and hippocampus, characterized spatial gene expression patterns in the Purkinje layer of mouse cerebellum, and defined the temporal evolution of cell type-specific responses in a mouse model of traumatic brain injury. These studies highlight how Slide-seq provides a scalable method for obtaining spatially resolved gene expression data at resolutions comparable to the sizes of individual cells.

DAX1/NR0B1 was expressed during mammalian gonadal development and gametogenesis before it was recruited to the eutherian X chromosome.

The nuclear receptor subfamily 0, group B, member 1 (NR0B1) gene is an orphan nuclear receptor that is X-linked in eutherian mammals and plays a critical role in the establishment and function of the hypothalamic-pituitary-adrenal-gonadal axis. Duplication or overexpression of NR0B1 in eutherian males causes male to female sex reversal, and mutation and deletions of NR0B1 cause testicular defects. Thus, gene dosage is critical for the function of NR0B1 in normal gonadogenesis. However, NR0B1 is autosomal in all noneutherian vertebrates, including marsupials and monotreme mammals, and two active copies of the gene are compatible with both male and female gonadal development. In the current study, we examined the evolution and expression of autosomal NR0B1 during gonadal development in a marsupial (the tammar wallaby) as compared to the role of its X-linked orthologues in a eutherian (the mouse). We show that NR0B1 underwent rapid evolutionary change when it relocated from its autosomal position in the nonmammalian vertebrates, monotremes, and marsupials to an X-linked location in eutherian mammals. Despite the acquisition of a novel genomic location and a unique N-terminal domain, NR0B1 protein distribution was remarkably similar between mice and marsupials both throughout gonadal development and during gamete formation. A conserved accumulation of NR0B1 protein was observed in developing oocytes, where its function appears to be critical in the early embryo, prior to zygotic genome activation. Together these findings suggest that NR0B1 had a conserved role in gonadogenesis that existed long before it moved to the X chromosome and despite undergoing significant evolutionary change.