Quantitative control of Ets1 dosage by a multi-enhancer hub promotes Th1 cell differentiation and protects from allergic inflammation.

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.

Allelic contribution of Nrxn1α to autism-relevant behavioral phenotypes in mice.

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.

Single-cell multi-omics analysis of human pancreatic islets reveals novel cellular states in type 1 diabetes.

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.

BET bromodomain protein inhibition reverses chimeric antigen receptor extinction and reinvigorates exhausted T cells in chronic lymphocytic leukemia.

Chimeric antigen receptor (CAR) T cells have induced remarkable antitumor responses in B cell malignancies. Some patients do not respond because of T cell deficiencies that hamper the expansion, persistence, and effector function of these cells. We used longitudinal immune profiling to identify phenotypic and pharmacodynamic changes in CD19-directed CAR T cells in patients with chronic lymphocytic leukemia (CLL). CAR expression maintenance was also investigated because this can affect response durability. CAR T cell failure was accompanied by preexisting T cell-intrinsic defects or dysfunction acquired after infusion. In a small subset of patients, CAR silencing was observed coincident with leukemia relapse. Using a small molecule inhibitor, we demonstrated that the bromodomain and extra-terminal (BET) family of chromatin adapters plays a role in downregulating CAR expression. BET protein blockade also ameliorated CAR T cell exhaustion as manifested by inhibitory receptor reduction, enhanced metabolic fitness, increased proliferative capacity, and enriched transcriptomic signatures of T cell reinvigoration. BET inhibition decreased levels of the TET2 methylcytosine dioxygenase, and forced expression of the TET2 catalytic domain eliminated the potency-enhancing effects of BET protein targeting in CAR T cells, providing a mechanism linking BET proteins and T cell dysfunction. Thus, modulating BET epigenetic readers may improve the efficacy of cell-based immunotherapies.

Exome-wide evaluation of rare coding variants using electronic health records identifies new gene-phenotype associations.

The clinical impact of rare loss-of-function variants has yet to be determined for most genes. Integration of DNA sequencing data with electronic health records (EHRs) could enhance our understanding of the contribution of rare genetic variation to human disease1. By leveraging 10,900 whole-exome sequences linked to EHR data in the Penn Medicine Biobank, we addressed the association of the cumulative effects of rare predicted loss-of-function variants for each individual gene on human disease on an exome-wide scale, as assessed using a set of diverse EHR phenotypes. After discovering 97 genes with exome-by-phenome-wide significant phenotype associations (P < 10-6), we replicated 26 of these in the Penn Medicine Biobank, as well as in three other medical biobanks and the population-based UK Biobank. Of these 26 genes, five had associations that have been previously reported and represented positive controls, whereas 21 had phenotype associations not previously reported, among which were genes implicated in glaucoma, aortic ectasia, diabetes mellitus, muscular dystrophy and hearing loss. These findings show the value of aggregating rare predicted loss-of-function variants into 'gene burdens' for identifying new gene-disease associations using EHR phenotypes in a medical biobank. We suggest that application of this approach to even larger numbers of individuals will provide the statistical power required to uncover unexplored relationships between rare genetic variation and disease phenotypes.

SARS-CoV-2 Cell Entry Factors ACE2 and TMPRSS2 Are Expressed in the Microvasculature and Ducts of Human Pancreas but Are Not Enriched in ß Cells.

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.

SARS-CoV-2 Cell Entry Factors ACE2 and TMPRSS2 are Expressed in the Pancreas but are Not Enriched in Islet Endocrine Cells.

Reports of new-onset diabetes and diabetic ketoacidosis in individuals with COVID-19 have led to the hypothesis that SARS-CoV-2, the virus that causes COVID-19, is directly cytotoxic to pancreatic islet ß cells. This would require binding and entry of SARS-CoV-2 into host ß cells via cell surface co-expression of ACE2 and TMPRSS2, the putative receptor and effector protease, respectively. To define ACE2 and TMPRSS2 expression in the human pancreas, we examined six transcriptional datasets from primary human islet cells and assessed protein expression by immunofluorescence in pancreata from donors with and without diabetes. ACE2 and TMPRSS2 transcripts were low or undetectable in pancreatic islet endocrine cells as determined by bulk or single cell RNA sequencing, and neither protein was detected in α or ß cells from these donors. Instead, ACE2 protein was expressed in the islet and exocrine tissue microvasculature and also found in a subset of pancreatic ducts, whereas TMPRSS2 protein was restricted to ductal cells. The absence of significant ACE2 and TMPRSS2 co-expression in islet endocrine cells reduces the likelihood that SARS-CoV-2 directly infects pancreatic islet ß cells through these cell entry proteins.

TooManyCells identifies and visualizes relationships of single-cell clades.

Identifying and visualizing transcriptionally similar cells is instrumental for accurate exploration of the cellular diversity revealed by single-cell transcriptomics. However, widely used clustering and visualization algorithms produce a fixed number of cell clusters. A fixed clustering 'resolution' hampers our ability to identify and visualize echelons of cell states. We developed TooManyCells, a suite of graph-based algorithms for efficient and unbiased identification and visualization of cell clades. TooManyCells introduces a visualization model built on a concept intentionally orthogonal to dimensionality-reduction methods. TooManyCells is also equipped with an efficient matrix-free divisive hierarchical spectral clustering different from prevalent single-resolution clustering methods. TooManyCells enables multiresolution and multifaceted exploration of single-cell clades. An advantage of this paradigm is the immediate detection of rare and common populations that outperforms popular clustering and visualization algorithms, as demonstrated using existing single-cell transcriptomic data sets and new data modeling drug-resistance acquisition in leukemic T cells.

Joint profiling of chromatin accessibility and CAR-T integration site analysis at population and single-cell levels.

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.

Genetic Variation in Type 1 Diabetes Reconfigures the 3D Chromatin Organization of T Cells and Alters Gene Expression.

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.

A mass spectrometry-based assay using metabolic labeling to rapidly monitor chromatin accessibility of modified histone proteins.

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.

Oncogenic Notch Promotes Long-Range Regulatory Interactions within Hyperconnected 3D Cliques.

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.

A mass spectrometry-based assay using metabolic labeling to rapidly monitor chromatin accessibility of modified histone proteins

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.

A mass spectrometry-based assay using metabolic labeling to rapidly monitor chromatin accessibility of modified histone proteins

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.

Long genes linked to autism spectrum disorders harbor broad enhancer-like chromatin domains.

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.

Biotin tagging of MeCP2 in mice reveals contextual insights into the Rett syndrome transcriptome.

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.

GR Extinguishes Inflamed Chromatin, and NF-κB Evacuates.

In this issue of Immunity,Oh et al. (2017) provide insight into the molecular effects of glucocorticoid receptor (GR) activation at a clinically relevant time point, after an inflammatory stimulus. They report that GR activation causes a global reduction in NF-κB binding, as well as time-dependent transcriptional effects.

The Crucial Role of DNA Methylation and MeCP2 in Neuronal Function.

A neuron is unique in its ability to dynamically modify its transcriptional output in response to synaptic activity while maintaining a core gene expression program that preserves cellular identity throughout a lifetime that is longer than almost every other cell type in the body. A contributing factor to the immense adaptability of a neuron is its unique epigenetic landscape that elicits locus-specific alterations in chromatin architecture, which in turn influences gene expression. One such epigenetic modification that is sensitive to changes in synaptic activity, as well as essential for maintaining cellular identity, is DNA methylation. The focus of this article is on the importance of DNA methylation in neuronal function, summarizing recent studies on critical players in the establishment of (the "writing"), the modification or erasure of (the "editing"), and the mediation of (the "reading") DNA methylation in neurodevelopment and neuroplasticity. One "reader" of DNA methylation in particular, methyl-CpG-binding protein 2 (MeCP2), is highlighted, given its undisputed importance in neuronal function.

Distinct cellular and molecular environments support aging-related DNA methylation changes in the substantia nigra.

AIM: We aimed to couple brain region-specific changes in global DNA methylation over aging to underlying cellular and molecular environments. MATERIALS & METHODS: We measured two major forms of DNA methylation and analyzed Dnmt, Tet and metabolite levels in the striatum and substantia nigra (SN) over aging in healthy male mice. RESULTS: The ratio of 5-hydroxymethylcytosine to 5-methylcytosine increases over aging in the SN, and 5-hydroxymethylcytosine increases preferentially in dopaminergic neurons. Additionally, this age-dependent alteration in methylation correlates with a reduction in the ratio of α-ketoglutarate to succinate in the SN. CONCLUSION: Distinct cellular and molecular environments correlate with aging-associated methylation changes in the SN, implicating this epigenetic mechanism in the susceptibility of this brain region to age-related cell loss.

Locus- and cell type-specific epigenetic switching during cellular differentiation in mammals.

BACKGROUND: Epigenomic reconfiguration, including changes in DNA methylation and histone modifications, is crucial for the differentiation of embryonic stem cells (ESCs) into somatic cells. However, the extent to which the epigenome is reconfigured and the interplay between components of the epigenome during cellular differentiation remain poorly defined. METHODS: We systematically analyzed and compared DNA methylation, various histone modification, and transcriptome profiles in ESCs with those of two distinct types of somatic cells from human and mouse. RESULTS: We found that global DNA methylation levels are lower in somatic cells compared to ESCs in both species. We also found that 80% of regions with histone modification occupancy differ between human ESCs and the two human somatic cell types. Approximately 70% of the reconfigurations in DNA methylation and histone modifications are locus- and cell type-specific. Intriguingly, the loss of DNA methylation is accompanied by the gain of different histone modifications in a locus- and cell type-specific manner. Further examination of transcriptional changes associated with epigenetic reconfiguration at promoter regions revealed an epigenetic switching for gene regulation-a transition from stable gene silencing mediated by DNA methylation in ESCs to flexible gene repression facilitated by repressive histone modifications in somatic cells. CONCLUSIONS: Our findings demonstrate that the epigenome is reconfigured in a locus- and cell type-specific manner and epigenetic switching is common during cellular differentiation in both human and mouse.