Engineered cell entry links receptor biology with single-cell genomics.

Cells communicate with each other via receptor-ligand interactions. Here, we describe lentiviral-mediated cell entry by engineered receptor-ligand interaction (ENTER) to display ligand proteins, deliver payloads, and record receptor specificity. We optimize ENTER to decode interactions between T cell receptor (TCR)-MHC peptides, antibody-antigen, and other receptor-ligand pairs. A viral presentation strategy allows ENTER to capture interactions between B cell receptor and any antigen. We engineer ENTER to deliver genetic payloads to antigen-specific T or B cells to selectively modulate cellular behavior in mixed populations. Single-cell readout of ENTER by RNA sequencing (ENTER-seq) enables multiplexed enumeration of antigen specificities, TCR clonality, cell type, and states of individual T cells. ENTER-seq of CMV-seropositive patient blood samples reveals the viral epitopes that drive effector memory T cell differentiation and inter-clonal vs. intra-clonal phenotypic diversity targeting the same epitope. ENTER technology enables systematic discovery of receptor specificity, linkage to cell fates, and antigen-specific cargo delivery.

B cell-specific XIST complex enforces X-inactivation and restrains atypical B cells.

The long non-coding RNA (lncRNA) XIST establishes X chromosome inactivation (XCI) in female cells in early development and thereafter is thought to be largely dispensable. Here, we show XIST is continually required in adult human B cells to silence a subset of X-linked immune genes such as TLR7. XIST-dependent genes lack promoter DNA methylation and require continual XIST-dependent histone deacetylation. XIST RNA-directed proteomics and CRISPRi screen reveal distinctive somatic cell-type-specific XIST complexes and identify TRIM28 that mediates Pol II pausing at promoters of X-linked genes in B cells. Single-cell transcriptome data of female patients with either systemic lupus erythematosus or COVID-19 infection revealed XIST dysregulation, reflected by escape of XIST-dependent genes, in CD11c+ atypical memory B cells (ABCs). XIST inactivation with TLR7 agonism suffices to promote isotype-switched ABCs. These results indicate cell-type-specific diversification and function for lncRNA-protein complexes and suggest expanded roles for XIST in sex-differences in biology and medicine.

Discovery and functional interrogation of SARS-CoV-2 RNA-host protein interactions.

SARS-CoV-2 is the cause of a pandemic with growing global mortality. Using comprehensive identification of RNA-binding proteins by mass spectrometry (ChIRP-MS), we identified 309 host proteins that bind the SARS-CoV-2 RNA during active infection. Integration of this data with ChIRP-MS data from three other RNA viruses defined viral specificity of RNA-host protein interactions. Targeted CRISPR screens revealed that the majority of functional RNA-binding proteins protect the host from virus-induced cell death, and comparative CRISPR screens across seven RNA viruses revealed shared and SARS-specific antiviral factors. Finally, by combining the RNA-centric approach and functional CRISPR screens, we demonstrated a physical and functional connection between SARS-CoV-2 and mitochondria, highlighting this organelle as a general platform for antiviral activity. Altogether, these data provide a comprehensive catalog of functional SARS-CoV-2 RNA-host protein interactions, which may inform studies to understand the host-virus interface and nominate host pathways that could be targeted for therapeutic benefit.

Genome-wide programmable transcriptional memory by CRISPR-based epigenome editing.

A general approach for heritably altering gene expression has the potential to enable many discovery and therapeutic efforts. Here, we present CRISPRoff-a programmable epigenetic memory writer consisting of a single dead Cas9 fusion protein that establishes DNA methylation and repressive histone modifications. Transient CRISPRoff expression initiates highly specific DNA methylation and gene repression that is maintained through cell division and differentiation of stem cells to neurons. Pairing CRISPRoff with genome-wide screens and analysis of chromatin marks establishes rules for heritable gene silencing. We identify single guide RNAs (sgRNAs) capable of silencing the large majority of genes including those lacking canonical CpG islands (CGIs) and reveal a wide targeting window extending beyond annotated CGIs. The broad ability of CRISPRoff to initiate heritable gene silencing even outside of CGIs expands the canonical model of methylation-based silencing and enables diverse applications including genome-wide screens, multiplexed cell engineering, enhancer silencing, and mechanistic exploration of epigenetic inheritance.

Divergent clonal differentiation trajectories of T cell exhaustion.

Chronic antigen exposure during viral infection or cancer promotes an exhausted T cell (Tex) state with reduced effector function. However, whether all antigen-specific T cell clones follow the same Tex differentiation trajectory remains unclear. Here, we generate a single-cell multiomic atlas of T cell exhaustion in murine chronic viral infection that redefines Tex phenotypic diversity, including two late-stage Tex subsets with either a terminal exhaustion (Texterm) or a killer cell lectin-like receptor-expressing cytotoxic (TexKLR) phenotype. We use paired single-cell RNA and T cell receptor sequencing to uncover clonal differentiation trajectories of Texterm-biased, TexKLR-biased or divergent clones that acquire both phenotypes. We show that high T cell receptor signaling avidity correlates with Texterm, whereas low avidity correlates with effector-like TexKLR fate. Finally, we identify similar clonal differentiation trajectories in human tumor-infiltrating lymphocytes. These findings reveal clonal heterogeneity in the T cell response to chronic antigen that influences Tex fates and persistence.

Differential usage of transcriptional repressor Zeb2 enhancers distinguishes adult and embryonic hematopoiesis.

The transcriptional repressor ZEB2 regulates development of many cell fates among somatic, neural, and hematopoietic lineages, but the basis for its requirement in these diverse lineages is unclear. Here, we identified a 400-basepair (bp) region located 165 kilobases (kb) upstream of the Zeb2 transcriptional start site (TSS) that binds the E proteins at several E-box motifs and was active in hematopoietic lineages. Germline deletion of this 400-bp region (Zeb2Δ-165mice) specifically prevented Zeb2 expression in hematopoietic stem cell (HSC)-derived lineages. Zeb2Δ-165 mice lacked development of plasmacytoid dendritic cells (pDCs), monocytes, and B cells. All macrophages in Zeb2Δ-165 mice were exclusively of embryonic origin. Using single-cell chromatin profiling, we identified a second Zeb2 enhancer located at +164-kb that was selectively active in embryonically derived lineages, but not HSC-derived ones. Thus, Zeb2 expression in adult, but not embryonic, hematopoiesis is selectively controlled by the -165-kb Zeb2 enhancer.

ecDNA hubs drive cooperative intermolecular oncogene expression.

Extrachromosomal DNA (ecDNA) is prevalent in human cancers and mediates high expression of oncogenes through gene amplification and altered gene regulation1. Gene induction typically involves cis-regulatory elements that contact and activate genes on the same chromosome2,3. Here we show that ecDNA hubs-clusters of around 10-100 ecDNAs within the nucleus-enable intermolecular enhancer-gene interactions to promote oncogene overexpression. ecDNAs that encode multiple distinct oncogenes form hubs in diverse cancer cell types and primary tumours. Each ecDNA is more likely to transcribe the oncogene when spatially clustered with additional ecDNAs. ecDNA hubs are tethered by the bromodomain and extraterminal domain (BET) protein BRD4 in a MYC-amplified colorectal cancer cell line. The BET inhibitor JQ1 disperses ecDNA hubs and preferentially inhibits ecDNA-derived-oncogene transcription. The BRD4-bound PVT1 promoter is ectopically fused to MYC and duplicated in ecDNA, receiving promiscuous enhancer input to drive potent expression of MYC. Furthermore, the PVT1 promoter on an exogenous episome suffices to mediate gene activation in trans by ecDNA hubs in a JQ1-sensitive manner. Systematic silencing of ecDNA enhancers by CRISPR interference reveals intermolecular enhancer-gene activation among multiple oncogene loci that are amplified on distinct ecDNAs. Thus, protein-tethered ecDNA hubs enable intermolecular transcriptional regulation and may serve as units of oncogene function and cooperative evolution and as potential targets for cancer therapy.

A fluorogenic array for temporally unlimited single-molecule tracking.

We describe three optical tags, ArrayG, ArrayD and ArrayG/N, for intracellular tracking of single molecules over milliseconds to hours. ArrayG is a fluorogenic tag composed of a green fluorescent protein-nanobody array and monomeric wild-type green fluorescent protein binders that are initially dim but brighten ~26-fold on binding with the array. By balancing the rates of binder production, photobleaching and stochastic binder exchange, we achieve temporally unlimited tracking of single molecules. High-speed tracking of ArrayG-tagged kinesins and integrins for thousands of frames reveals novel dynamical features. Tracking of single histones at 0.5 Hz for >1 hour with the import competent ArrayG/N tag shows that chromosomal loci behave as Rouse polymers with visco-elastic memory and exhibit a non-Gaussian displacement distribution. ArrayD, based on a dihydrofolate reductase nanobody array and dihydrofolate reductase-fluorophore binder, enables dual-color imaging. The arrays combine brightness, fluorogenicity, fluorescence replenishment and extended fluorophore choice, opening new avenues for tracking single molecules in living cells.

Rapid disorganization of mechanically interacting systems of mammary acini.

Cells and multicellular structures can mechanically align and concentrate fibers in their ECM environment and can sense and respond to mechanical cues by differentiating, branching, or disorganizing. Here we show that mammary acini with compromised structural integrity can interconnect by forming long collagen lines. These collagen lines then coordinate and accelerate transition to an invasive phenotype. Interacting acini begin to disorganize within 12.5 ± 4.7 h in a spatially coordinated manner, whereas acini that do not interact mechanically with other acini disorganize more slowly (in 21.8 ± 4.1 h) and to a lesser extent (P < 0.0001). When the directed mechanical connections between acini were cut with a laser, the acini reverted to a slowly disorganizing phenotype. When acini were fully mechanically isolated from other acini and also from the bulk gel by box-cuts with a side length <900 µm, transition to an invasive phenotype was blocked in 20 of 20 experiments, regardless of waiting time. Thus, pairs or groups of mammary acini can interact mechanically over long distances through the collagen matrix, and these directed mechanical interactions facilitate transition to an invasive phenotype.

Cadherin-dependent mechanotransduction depends on ligand identity but not affinity.

This study investigates the relationship between classical cadherin binding affinities and mechanotransduction through cadherin-mediated adhesions. The mechanical properties of cadherin-dependent intercellular junctions are generally attributed to differences in the binding affinities of classical cadherin subtypes that contribute to cohesive energies between cells. However, cell mechanics and mechanotransduction may also regulate intercellular contacts. We used micropipette measurements to quantify the two-dimensional affinities of cadherins at the cell surface, and two complementary mechanical measurements to assess ligand-dependent mechanotransduction through cadherin adhesions. At the cell surface, the classical cadherins investigated in this study form both homophilic and heterophilic bonds with two-dimensional affinities that differ by less than threefold. In contrast, mechanotransduction through cadherin adhesions is strongly ligand dependent such that homophilic, but not heterophilic ligation mediates mechanotransduction, independent of the cadherin binding affinity. These findings suggest that ligand-selective mechanotransduction may supersede differences in cadherin binding affinities in regulating intercellular contacts.

Cadherin point mutations alter cell sorting and modulate GTPase signaling.

This study investigated the impact of cadherin binding differences on both cell sorting and GTPase activation. The use of N-terminal domain point mutants of Xenopus C-cadherin enabled us to quantify binding differences and determine their effects on cadherin-dependent functions without any potential complications arising as a result of differences in cytodomain interactions. Dynamic cell-cell binding measurements carried out with the micropipette manipulation technique quantified the impact of these mutations on the two-dimensional binding affinities and dissociation rates of cadherins in the native context of the cell membrane. Pairwise binding affinities were compared with in vitro cell-sorting specificity and ligation-dependent GTPase signaling. Two-dimensional affinity differences greater than five-fold correlated with cadherin-dependent in vitro cell segregation, but smaller differences failed to induce cell sorting. Comparison of the binding affinities with GTPase signaling amplitudes further demonstrated that differential binding also proportionally modulates intracellular signaling. These results show that differential cadherin affinities have broader functional consequences than merely controlling cell-cell cohesion.

Bidirectional epigenetic editing reveals hierarchies in gene regulation.

CRISPR perturbation methods are limited in their ability to study non-coding elements and genetic interactions. In this study, we developed a system for bidirectional epigenetic editing, called CRISPRai, in which we apply activating (CRISPRa) and repressive (CRISPRi) perturbations to two loci simultaneously in the same cell. We developed CRISPRai Perturb-seq by coupling dual perturbation gRNA detection with single-cell RNA sequencing, enabling study of pooled perturbations in a mixed single-cell population. We applied this platform to study the genetic interaction between two hematopoietic lineage transcription factors, SPI1 and GATA1, and discovered novel characteristics of their co-regulation on downstream target genes, including differences in SPI1 and GATA1 occupancy at genes that are regulated through different modes. We also studied the regulatory landscape of IL2 (interleukin-2) in Jurkat T cells, primary T cells and chimeric antigen receptor (CAR) T cells and elucidated mechanisms of enhancer-mediated IL2 gene regulation. CRISPRai facilitates investigation of context-specific genetic interactions, provides new insights into gene regulation and will enable exploration of non-coding disease-associated variants.

Spatiotemporal co-dependency between macrophages and exhausted CD8+ T cells in cancer.

T cell exhaustion is a major impediment to antitumor immunity. However, it remains elusive how other immune cells in the tumor microenvironment (TME) contribute to this dysfunctional state. Here, we show that the biology of tumor-associated macrophages (TAMs) and exhausted T cells (Tex) in the TME is extensively linked. We demonstrate that in vivo depletion of TAMs reduces exhaustion programs in tumor-infiltrating CD8+ T cells and reinvigorates their effector potential. Reciprocally, transcriptional and epigenetic profiling reveals that Tex express factors that actively recruit monocytes to the TME and shape their differentiation. Using lattice light sheet microscopy, we show that TAM and CD8+ T cells engage in unique, long-lasting, antigen-specific synaptic interactions that fail to activate T cells but prime them for exhaustion, which is then accelerated in hypoxic conditions. Spatially resolved sequencing supports a spatiotemporal self-enforcing positive feedback circuit that is aligned to protect rather than destroy a tumor.

Inducible lncRNA transgenic mice reveal continual role of HOTAIR in promoting breast cancer metastasis.

HOTAIR is a 2.2-kb long noncoding RNA (lncRNA) whose dysregulation has been linked to oncogenesis, defects in pattern formation during early development, and irregularities during the process of epithelial-to-mesenchymal transition (EMT). However, the oncogenic transformation determined by HOTAIR in vivo and its impact on chromatin dynamics are incompletely understood. Here, we generate a transgenic mouse model with doxycycline-inducible expression of human HOTAIR in the context of the MMTV-PyMT breast cancer-prone background to systematically interrogate the cellular mechanisms by which human HOTAIR lncRNA acts to promote breast cancer progression. We show that sustained high levels of HOTAIR over time increased breast metastatic capacity and invasiveness in breast cancer cells, promoting migration and subsequent metastasis to the lung. Subsequent withdrawal of HOTAIR overexpression reverted the metastatic phenotype, indicating oncogenic lncRNA addiction. Furthermore, HOTAIR overexpression altered both the cellular transcriptome and chromatin accessibility landscape of multiple metastasis-associated genes and promoted EMT. These alterations are abrogated within several cell cycles after HOTAIR expression is reverted to basal levels, indicating an erasable lncRNA-associated epigenetic memory. These results suggest that a continual role for HOTAIR in programming a metastatic gene regulatory program. Targeting HOTAIR lncRNA may potentially serve as a therapeutic strategy to ameliorate breast cancer progression.

Genome-wide CRISPR screens of T cell exhaustion identify chromatin remodeling factors that limit T cell persistence.

T cell exhaustion limits antitumor immunity, but the molecular determinants of this process remain poorly understood. Using a chronic stimulation assay, we performed genome-wide CRISPR-Cas9 screens to systematically discover regulators of T cell exhaustion, which identified an enrichment of epigenetic factors. In vivo CRISPR screens in murine and human tumor models demonstrated that perturbation of the INO80 and BAF chromatin remodeling complexes improved T cell persistence in tumors. In vivo Perturb-seq revealed distinct transcriptional roles of each complex and that depletion of canonical BAF complex members, including Arid1a, resulted in the maintenance of an effector program and downregulation of exhaustion-related genes in tumor-infiltrating T cells. Finally, Arid1a depletion limited the acquisition of exhaustion-associated chromatin accessibility and led to improved antitumor immunity. In summary, we provide an atlas of the genetic regulators of T cell exhaustion and demonstrate that modulation of epigenetic state can improve T cell responses in cancer immunotherapy.

The evolutionary dynamics of extrachromosomal DNA in human cancers.

Oncogene amplification on extrachromosomal DNA (ecDNA) is a common event, driving aggressive tumor growth, drug resistance and shorter survival. Currently, the impact of nonchromosomal oncogene inheritance-random identity by descent-is poorly understood. Also unclear is the impact of ecDNA on somatic variation and selection. Here integrating theoretical models of random segregation, unbiased image analysis, CRISPR-based ecDNA tagging with live-cell imaging and CRISPR-C, we demonstrate that random ecDNA inheritance results in extensive intratumoral ecDNA copy number heterogeneity and rapid adaptation to metabolic stress and targeted treatment. Observed ecDNAs benefit host cell survival or growth and can change within a single cell cycle. ecDNA inheritance can predict, a priori, some of the aggressive features of ecDNA-containing cancers. These properties are facilitated by the ability of ecDNA to rapidly adapt genomes in a way that is not possible through chromosomal oncogene amplification. These results show how the nonchromosomal random inheritance pattern of ecDNA contributes to poor outcomes for patients with cancer.

B cell-specific XIST complex enforces X-inactivation and restrains atypical B cells.

The long noncoding RNA (lncRNA) XIST establishes X chromosome inactivation (XCI) in female cells in early development and thereafter is thought to be largely dispensable. Here we show XIST is continually required in adult human B cells to silence a subset of X-linked immune genes such as TLR7 . XIST-dependent genes lack promoter DNA methylation and require continual XIST-dependent histone deacetylation. XIST RNA-directed proteomics and CRISPRi screen reveal distinctive somatic cell-specific XIST complexes, and identify TRIM28 that mediates Pol II pausing at promoters of X-linked genes in B cells. XIST dysregylation, reflected by escape of XIST-dependent genes, occurs in CD11c+ atypical memory B cells across single-cell transcriptome data in patients with female-biased autoimmunity and COVID-19 infection. XIST inactivation with TLR7 agonism suffices to promote isotype-switched atypical B cells. These results suggest cell-type-specific diversification of lncRNA-protein complexes increase lncRNA functionalities, and expand roles for XIST in sex-differences in biology and medicine. HIGHLIGHTS: XIST prevents escape of genes with DNA hypomethylated promoters in B cells.XIST maintains X-inactivation through continuous deacetylation of H3K27ac.XIST ChIRP-MS and allelic CRISPRi screen reveal a B cell-specific XIST cofactor TRIM28.XIST loss and TLR7 stimulation promotes CD11c+ atypical B cell formation.

A human mutation in STAT3 promotes type 1 diabetes through a defect in CD8+ T cell tolerance.

Naturally occurring cases of monogenic type 1 diabetes (T1D) help establish direct mechanisms driving this complex autoimmune disease. A recently identified de novo germline gain-of-function (GOF) mutation in the transcriptional regulator STAT3 was found to cause neonatal T1D. We engineered a novel knock-in mouse incorporating this highly diabetogenic human STAT3 mutation (K392R) and found that these mice recapitulated the human autoimmune diabetes phenotype. Paired single-cell TCR and RNA sequencing revealed that STAT3-GOF drives proliferation and clonal expansion of effector CD8+ cells that resist terminal exhaustion. Single-cell ATAC-seq showed that these effector T cells are epigenetically distinct and have differential chromatin architecture induced by STAT3-GOF. Analysis of islet TCR clonotypes revealed a CD8+ cell reacting against known antigen IGRP, and STAT3-GOF in an IGRP-reactive TCR transgenic model demonstrated that STAT3-GOF intrinsic to CD8+ cells is sufficient to accelerate diabetes onset. Altogether, these findings reveal a diabetogenic CD8+ T cell response that is restrained in the presence of normal STAT3 activity and drives diabetes pathogenesis.

Allosteric cross talk between cadherin extracellular domains.

Atomic force microscopy and surface force apparatus measurements determined the functional impact of the cadherin point mutation W2A and domain deletion mutations on C-cadherin binding signatures. Direct comparison of results obtained using both experimental approaches demonstrates that C-cadherin ectodomains form multiple independent bonds that require different structural regions. The results presented reveal significant interdomain cross talk. They further demonstrate that the mutation W2A not only abolishes adhesion between N-terminal domains, but allosterically modulates other binding states that require functional domains distal to the N-terminal binding site. Such allosteric effects may play a prominent role in modulating adhesion by Type I classic cadherins, cadherin oligomerization at junctional contacts, and propagation of binding information to the cytoplasmic region.

Concerted localization-resets precede YAP-dependent transcription.

Yes-associated protein 1 (YAP) is a transcriptional regulator with critical roles in mechanotransduction, organ size control, and regeneration. Here, using advanced tools for real-time visualization of native YAP and target gene transcription dynamics, we show that a cycle of fast exodus of nuclear YAP to the cytoplasm followed by fast reentry to the nucleus ("localization-resets") activates YAP target genes. These "resets" are induced by calcium signaling, modulation of actomyosin contractility, or mitosis. Using nascent-transcription reporter knock-ins of YAP target genes, we show a strict association between these resets and downstream transcription. Oncogenically-transformed cell lines lack localization-resets and instead show dramatically elevated rates of nucleocytoplasmic shuttling of YAP, suggesting an escape from compartmentalization-based control. The single-cell localization and transcription traces suggest that YAP activity is not a simple linear function of nuclear enrichment and point to a model of transcriptional activation based on nucleocytoplasmic exchange properties of YAP.