Chimeric antigen receptor (CAR)-T immunotherapy has yielded impressive results in several B cell malignancies, establishing itself as a powerful means to redirect the natural properties of T lymphocytes. In this strategy, the T cell genome is modified by the integration of lentiviral vectors encoding CAR that direct tumor cell killing. However, this therapeutic approach is often limited by the extent of CAR-T cell expansion in vivo. A major outstanding question is whether or not CAR-T integration itself enhances the proliferative competence of individual T cells by rewiring their regulatory landscape. To address this question, it is critical to define the identity of an individual CAR-T cell and simultaneously chart where the CAR-T vector integrates into the genome. Here, we report the development of a method called EpiVIA (https://github.com/VahediLab/epiVIA) for the joint profiling of the chromatin accessibility and lentiviral integration site analysis at the population and single-cell levels. We validate our technique in clonal cells with previously defined integration sites and further demonstrate the ability to measure lentiviral integration sites and chromatin accessibility of host and viral genomes at the single-cell resolution in CAR-T cells. We anticipate that EpiVIA will enable the single-cell deconstruction of gene regulation during CAR-T therapy, leading to the discovery of cellular factors associated with durable treatment.

Genetics is a major determinant of susceptibility to autoimmune disorders. Here, we examined whether genome organization provides resilience or susceptibility to sequence variations, and how this would contribute to the molecular etiology of an autoimmune disease. We generated high-resolution maps of linear and 3D genome organization in thymocytes of NOD mice, a model of type 1 diabetes (T1D), and the diabetes-resistant C57BL/6 mice. Multi-enhancer interactions formed at genomic regions harboring genes with prominent roles in T cell development in both strains. However, diabetes risk-conferring loci coalesced enhancers and promoters in NOD, but not C57BL/6 thymocytes. 3D genome mapping of NODxC57BL/6 F1 thymocytes revealed that genomic misfolding in NOD mice is mediated in cis. Moreover, immune cells infiltrating the pancreas of humans with T1D exhibited increased expression of genes located on misfolded loci in mice. Thus, genetic variation leads to altered 3D chromatin architecture and associated changes in gene expression that may underlie autoimmune pathology.

Mutations in MECP2 cause Rett syndrome (RTT), an X-linked neurological disorder characterized by regressive loss of neurodevelopmental milestones and acquired psychomotor deficits. However, the cellular heterogeneity of the brain impedes an understanding of how MECP2 mutations contribute to RTT. Here we developed a Cre-inducible method for cell-type-specific biotin tagging of MeCP2 in mice. Combining this approach with an allelic series of knock-in mice carrying frequent RTT-associated mutations (encoding T158M and R106W) enabled the selective profiling of RTT-associated nuclear transcriptomes in excitatory and inhibitory cortical neurons. We found that most gene-expression changes were largely specific to each RTT-associated mutation and cell type. Lowly expressed cell-type-enriched genes were preferentially disrupted by MeCP2 mutations, with upregulated and downregulated genes reflecting distinct functional categories. Subcellular RNA analysis in MeCP2-mutant neurons further revealed reductions in the nascent transcription of long genes and uncovered widespread post-transcriptional compensation at the cellular level. Finally, we overcame X-linked cellular mosaicism in female RTT models and identified distinct gene-expression changes between neighboring wild-type and mutant neurons, providing contextual insights into RTT etiology that support personalized therapeutic interventions.

Copy number variations (CNVs) in the Neurexin 1 (NRXN1) gene, which encodes a presynaptic protein involved in neurotransmitter release, are some of the most frequently observed single-gene variants associated with autism spectrum disorder (ASD). To address the functional contribution of NRXN1 CNVs to behavioral phenotypes relevant to ASD, we carried out systematic behavioral phenotyping of an allelic series of Nrxn1 mouse models: one carrying promoter and exon 1 deletion abolishing Nrxn1α transcription, one carrying exon 9 deletion disrupting Nrxn1α protein translation, and one carrying an intronic deletion with no observable effect on Nrxn1α expression. We found that homozygous loss of Nrxn1α resulted in enhanced aggression in males, reduced affiliative social behaviors in females, and significantly altered circadian activities in both sexes. Heterozygous or homozygous loss of Nrxn1α affected the preference for social novelty in male mice, and notably, enhanced repetitive motor skills and motor coordination in both sexes. In contrast, mice bearing an intronic deletion of Nrxn1 did not display alterations in any of the behaviors assessed. These findings demonstrate the importance of Nrxn1α gene dosage in regulating social, circadian, and motor functions, and the variables of sex and genomic positioning of CNVs in the expression of autism-related phenotypes. Importantly, mice with heterozygous loss of Nrxn1, as found in numerous autistic individuals, show an elevated propensity to manifest autism-related phenotypes, supporting the use of models with this genomic architecture to study ASD etiology and assess additional genetic variants associated with autism.

Type 1 diabetes (T1D) is an autoimmune disease in which immune cells destroy insulin-producing beta cells. The aetiology of this complex disease is dependent on the interplay of multiple heterogeneous cell types in the pancreatic environment. Here, we provide a single-cell atlas of pancreatic islets of 24 T1D, autoantibody-positive and nondiabetic organ donors across multiple quantitative modalities including ~80,000 cells using single-cell transcriptomics, ~7,000,000 cells using cytometry by time of flight and ~1,000,000 cells using in situ imaging mass cytometry. We develop an advanced integrative analytical strategy to assess pancreatic islets and identify canonical cell types. We show that a subset of exocrine ductal cells acquires a signature of tolerogenic dendritic cells in an apparent attempt at immune suppression in T1D donors. Our multimodal analyses delineate cell types and processes that may contribute to T1D immunopathogenesis and provide an integrative procedure for exploration and discovery of human pancreatic function.

Histone post-translational modifications (PTMs) contribute to chromatin accessibility due to their chemical properties and their ability to recruit enzymes responsible for DNA readout and chromatin remodeling. To date, more than 400 different histone PTMs and thousands of combinations of PTMs have been identified, the vast majority with still unknown biological function. Identification and quantification of histone PTMs has become routine in mass spectrometry (MS) but, since raising antibodies for each PTM in a study can be prohibitive, lots of potential is lost from MS datasets when uncharacterized PTMs are found to be significantly regulated. We developed an assay that uses metabolic labeling and MS to associate chromatin accessibility with histone PTMs and their combinations. The labeling is achieved by spiking in the cell media a 5x concentration of stable isotope labeled arginine and allow cells to grow for at least one cell cycle. We quantified the labeling incorporation of about 200 histone peptides with a proteomics workflow, and we confirmed that peptides carrying PTMs with extensively characterized roles in active transcription or gene silencing were in highly or poorly labeled forms, respectively. Data were further validated using next-generation sequencing to assess the transcription rate of chromatin regions modified with five selected PTMs. Furthermore, we quantified the labeling rate of peptides carrying co-existing PTMs, proving that this method is suitable for combinatorial PTMs. We focus on the abundant bivalent mark H3K27me3K36me2, showing that H3K27me3 dominantly represses histone swapping rate even in the presence of the more permissive PTM H3K36me2. Together, we envision this method will help to generate hypotheses regarding histone PTM functions and, potentially, elucidate the role of combinatorial histone codes.

Genetic variants associated with autism spectrum disorders (ASDs) are enriched in genes encoding synaptic proteins and chromatin regulators. Although the role of synaptic proteins in ASDs is widely studied, the mechanism by which chromatin regulators contribute to ASD risk remains poorly understood. Upon profiling and analyzing the transcriptional and epigenomic features of genes expressed in the cortex, we uncovered a unique set of long genes that contain broad enhancer-like chromatin domains (BELDs) spanning across their entire gene bodies. Analyses of these BELD genes show that they are highly transcribed with frequent RNA polymerase II (Pol II) initiation and low Pol II pausing, and they exhibit frequent chromatin-chromatin interactions within their gene bodies. These BELD features are conserved from rodents to humans, are enriched in genes involved in synaptic function, and appear post-natally concomitant with synapse development. Importantly, we find that BELD genes are highly implicated in neurodevelopmental disorders, particularly ASDs, and that their expression is preferentially down-regulated in individuals with idiopathic autism. Finally, we find that the transcription of BELD genes is particularly sensitive to alternations in ASD-associated chromatin regulators. These findings suggest that the epigenomic regulation of BELD genes is important for post-natal cortical development and lend support to a model by which mutations in chromatin regulators causally contribute to ASDs by preferentially impairing BELD gene transcription.

Multi-enhancer hubs are spatial clusters of enhancers present across numerous developmental programs. Here, we studied the functional relevance of these three-dimensional structures in T cell biology. Mathematical modeling identified a highly connected multi-enhancer hub at the Ets1 locus, comprising a noncoding regulatory element that was a hotspot for sequence variation associated with allergic disease in humans. Deletion of this regulatory element in mice revealed that the multi-enhancer connectivity was dispensable for T cell development but required for CD4+ T helper 1 (Th1) differentiation. These mice were protected from Th1-mediated colitis but exhibited overt allergic responses. Mechanistically, the multi-enhancer hub controlled the dosage of Ets1 that was required for CTCF recruitment and assembly of Th1-specific genome topology. Our findings establish a paradigm wherein multi-enhancer hubs control cellular competence to respond to an inductive cue through quantitative control of gene dosage and provide insight into how sequence variation within noncoding elements at the Ets1 locus predisposes individuals to allergic responses.

Amyotrophic lateral sclerosis (ALS) is a devastating neurodegenerative disease whose causes are still poorly understood. To identify additional genetic risk factors, we assessed the role of de novo mutations in ALS by sequencing the exomes of 47 ALS patients and both of their unaffected parents (n = 141 exomes). We found that amino acid-altering de novo mutations were enriched in genes encoding chromatin regulators, including the neuronal chromatin remodeling complex (nBAF) component SS18L1 (also known as CREST). CREST mutations inhibited activity-dependent neurite outgrowth in primary neurons, and CREST associated with the ALS protein FUS. These findings expand our understanding of the ALS genetic landscape and provide a resource for future studies into the pathogenic mechanisms contributing to sporadic ALS.

Chromatin loops enable transcription-factor-bound distal enhancers to interact with their target promoters to regulate transcriptional programs. Although developmental transcription factors such as active forms of Notch can directly stimulate transcription by activating enhancers, the effect of their oncogenic subversion on the 3D organization of cancer genomes is largely undetermined. By mapping chromatin looping genome-wide in Notch-dependent triple-negative breast cancer and B cell lymphoma, we show that beyond the well-characterized role of Notch as an activator of distal enhancers, Notch regulates its direct target genes by instructing enhancer repositioning. Moreover, a large fraction of Notch-instructed regulatory loops form highly interacting enhancer and promoter spatial clusters termed "3D cliques." Loss- and gain-of-function experiments show that Notch preferentially targets hyperconnected 3D cliques that regulate the expression of crucial proto-oncogenes. Our observations suggest that oncogenic hijacking of developmental transcription factors can dysregulate transcription through widespread effects on the spatial organization of cancer genomes.

Isolated reports of new-onset diabetes in individuals with COVID-19 have led to the hypothesis that SARS-CoV-2 is directly cytotoxic to pancreatic islet β cells. This would require binding and entry of SARS-CoV-2 into β cells via co-expression of its canonical cell entry factors, angiotensin-converting enzyme 2 (ACE2) and transmembrane serine protease 2 (TMPRSS2); however, their expression in human pancreas has not been clearly defined. We analyzed six transcriptional datasets of primary human islet cells and found that ACE2 and TMPRSS2 were not co-expressed in single β cells. In pancreatic sections, ACE2 and TMPRSS2 protein was not detected in β cells from donors with and without diabetes. Instead, ACE2 protein was expressed in islet and exocrine tissue microvasculature and in a subset of pancreatic ducts, whereas TMPRSS2 protein was restricted to ductal cells. These findings reduce the likelihood that SARS-CoV-2 directly infects β cells in vivo through ACE2 and TMPRSS2.