Dynamic Foxp3-chromatin interaction controls tunable Treg cell function.

Nuclear factor Foxp3 determines regulatory T (Treg) cell fate and function via mechanisms that remain unclear. Here, we investigate the nature of Foxp3-mediated gene regulation in suppressing autoimmunity and antitumor immune response. Contrasting with previous models, we find that Foxp3-chromatin binding is regulated by Treg activation states, tumor microenvironment, and antigen and cytokine stimulations. Proteomics studies uncover dynamic proteins within Foxp3 proximity upon TCR or IL-2 receptor signaling in vitro, reflecting intricate interactions among Foxp3, signal transducers, and chromatin. Pharmacological inhibition and genetic knockdown experiments indicate that NFAT and AP-1 protein Batf are required for enhanced Foxp3-chromatin binding in activated Treg cells and tumor-infiltrating Treg cells to modulate target gene expression. Furthermore, mutations at the Foxp3 DNA-binding domain destabilize the Foxp3-chromatin association. These representative settings delineate context-dependent Foxp3-chromatin interaction, suggesting that Foxp3 associates with chromatin by hijacking DNA-binding proteins resulting from Treg activation or differentiation, which is stabilized by direct Foxp3-DNA binding, to dynamically regulate Treg cell function according to immunological contexts.

Antitumor progenitor exhausted CD8+ T cells are sustained by TCR engagement.

The durability of an antitumor immune response is mediated in part by the persistence of progenitor exhausted CD8+ T cells (Tpex). Tpex serve as a resource for replenishing effector T cells and preserve their quantity through self-renewal. However, it is unknown how T cell receptor (TCR) engagement affects the self-renewal capacity of Tpex in settings of continued antigen exposure. Here we use a Lewis lung carcinoma model that elicits either optimal or attenuated TCR signaling in CD8+ T cells to show that formation of Tpex in tumor-draining lymph nodes and their intratumoral persistence is dependent on optimal TCR engagement. Notably, attenuated TCR stimulation accelerates the terminal differentiation of optimally primed Tpex. This TCR-reinforced Tpex development and self-renewal is coupled to proximal positioning to dendritic cells and epigenetic imprinting involving increased chromatin accessibility at Egr2 and Tcf1 target loci. Collectively, this study highlights the critical function of TCR engagement in sustaining Tpex during tumor progression.

Spatial transcriptomics reveals the distinct organization of mouse prefrontal cortex and neuronal subtypes regulating chronic pain.

The prefrontal cortex (PFC) is a complex brain region that regulates diverse functions ranging from cognition, emotion and executive action to even pain processing. To decode the cellular and circuit organization of such diverse functions, we employed spatially resolved single-cell transcriptome profiling of the adult mouse PFC. Results revealed that PFC has distinct cell-type composition and gene-expression patterns relative to neighboring cortical areas-with neuronal excitability-regulating genes differently expressed. These cellular and molecular features are further segregated within PFC subregions, alluding to the subregion-specificity of several PFC functions. PFC projects to major subcortical targets through combinations of neuronal subtypes, which emerge in a target-intrinsic fashion. Finally, based on these features, we identified distinct cell types and circuits in PFC underlying chronic pain, an escalating healthcare challenge with limited molecular understanding. Collectively, this comprehensive map will facilitate decoding of discrete molecular, cellular and circuit mechanisms underlying specific PFC functions in health and disease.

SNIP1 and PRC2 coordinate cell fates of neural progenitors during brain development.

Stem cell survival versus death is a developmentally programmed process essential for morphogenesis, sizing, and quality control of genome integrity and cell fates. Cell death is pervasive during development, but its programming is little known. Here, we report that Smad nuclear interacting protein 1 (SNIP1) promotes neural progenitor cell survival and neurogenesis and is, therefore, integral to brain development. The SNIP1-depleted brain exhibits dysplasia with robust induction of caspase 9-dependent apoptosis. Mechanistically, SNIP1 regulates target genes that promote cell survival and neurogenesis, and its activities are influenced by TGFß and NFκB signaling pathways. Further, SNIP1 facilitates the genomic occupancy of Polycomb complex PRC2 and instructs H3K27me3 turnover at target genes. Depletion of PRC2 is sufficient to reduce apoptosis and brain dysplasia and to partially restore genetic programs in the SNIP1-depleted brain in vivo. These findings suggest a loci-specific regulation of PRC2 and H3K27 marks to toggle cell survival and death in the developing brain.

Phosphorylation of 53BP1 by ATM enforce neurodevelopmental programs in cortical organoids.

53BP1 is a well-established DNA damage repair factor recently shown to regulate gene expression and critically influence tumor suppression and neural development. For gene regulation, how 53BP1 is regulated remains unclear. Here, we showed that 53BP1-serine 25 phosphorylation by ATM is required for neural progenitor cell proliferation and neuronal differentiation in cortical organoids. 53BP1-serine 25 phosphorylation dynamics controls 53BP1 target genes for neuronal differentiation and function, cellular response to stress, and apoptosis. Beyond 53BP1, ATM is required for phosphorylation of factors in neuronal differentiation, cytoskeleton, p53 regulation, and ATM, BNDF, and WNT signaling pathways for cortical organoid differentiation. Overall, our data suggest that 53BP1 and ATM control key genetic programs required for human cortical development.

Interrogating bromodomain inhibitor resistance in KMT2A-rearranged leukemia through combinatorial CRISPR screens.

Bromo- and extra-terminal domain inhibitors (BETi) have exhibited therapeutic activities in many cancers. However, the mechanisms controlling BETi response and resistance are not well understood. We conducted genome-wide loss-of-function CRISPR screens using BETi-treated KMT2A-rearranged (KMT2A-r) cell lines. We revealed that Speckle-type POZ protein (SPOP) gene (Speckle Type BTB/POZ Protein) deficiency caused significant BETi resistance, which was further validated in cell lines and xenograft models. Proteomics analysis and a kinase-vulnerability CRISPR screen indicated that cells treated with BETi are sensitive to GSK3 perturbation. Pharmaceutical inhibition of GSK3 reversed the BETi-resistance phenotype. Based on this observation, a combination therapy regimen inhibiting both BET and GSK3 was developed to impede KMT2A-r leukemia progression in patient-derived xenografts in vivo. Our results revealed molecular mechanisms underlying BETi resistance and a promising combination treatment regimen of ABBV-744 and CHIR-98014 by utilizing unique ex vivo and in vivo KMT2A-r PDX models.

Auxin-inducible degron 2 system deciphers functions of CTCF domains in transcriptional regulation.

BACKGROUND: CTCF is a well-established chromatin architectural protein that also plays various roles in transcriptional regulation. While CTCF biology has been extensively studied, how the domains of CTCF function to regulate transcription remains unknown. Additionally, the original auxin-inducible degron 1 (AID1) system has limitations in investigating the function of CTCF. RESULTS: We employ an improved auxin-inducible degron technology, AID2, to facilitate the study of acute depletion of CTCF while overcoming the limitations of the previous AID system. As previously observed through the AID1 system and steady-state RNA analysis, the new AID2 system combined with SLAM-seq confirms that CTCF depletion leads to modest nascent and steady-state transcript changes. A CTCF domain sgRNA library screening identifies the zinc finger (ZF) domain as the region within CTCF with the most functional relevance, including ZFs 1 and 10. Removal of ZFs 1 and 10 reveals genomic regions that independently require these ZFs for DNA binding and transcriptional regulation. Notably, loci regulated by either ZF1 or ZF10 exhibit unique CTCF binding motifs specific to each ZF. CONCLUSIONS: By extensively comparing the AID1 and AID2 systems for CTCF degradation in SEM cells, we confirm that AID2 degradation is superior for achieving miniAID-tagged protein degradation without the limitations of the AID1 system. The model we create that combines AID2 depletion of CTCF with exogenous overexpression of CTCF mutants allows us to demonstrate how peripheral ZFs intricately orchestrate transcriptional regulation in a cellular context for the first time.

eccDNAs are apoptotic products with high innate immunostimulatory activity.

Extrachromosomal circular DNA elements (eccDNAs) have been described in the literature for several decades, and are known for their broad existence across different species1,2. However, their biogenesis and functions are largely unknown. By developing a new circular DNA enrichment method, here we purified and sequenced full-length eccDNAs with Nanopore sequencing. We found that eccDNAs map across the entire genome in a close to random manner, suggesting a biogenesis mechanism of random ligation of genomic DNA fragments. Consistent with this idea, we found that apoptosis inducers can increase eccDNA generation, which is dependent on apoptotic DNA fragmentation followed by ligation by DNA ligase 3. Importantly, we demonstrated that eccDNAs can function as potent innate immunostimulants in a manner that is independent of eccDNA sequence but dependent on eccDNA circularity and the cytosolic DNA sensor Sting. Collectively, our study not only revealed the origin, biogenesis and immunostimulant function of eccDNAs but also uncovered their sensing pathway and potential clinical implications in immune response.

Acute depletion of CTCF rewires genome-wide chromatin accessibility.

BACKGROUND: The transcription factor CTCF appears indispensable in defining topologically associated domain boundaries and maintaining chromatin loop structures within these domains, supported by numerous functional studies. However, acute depletion of CTCF globally reduces chromatin interactions but does not significantly alter transcription. RESULTS: Here, we systematically integrate multi-omics data including ATAC-seq, RNA-seq, WGBS, Hi-C, Cut&Run, and CRISPR-Cas9 survival dropout screens, and time-solved deep proteomic and phosphoproteomic analyses in cells carrying auxin-induced degron at endogenous CTCF locus. Acute CTCF protein degradation markedly rewires genome-wide chromatin accessibility. Increased accessible chromatin regions are frequently located adjacent to CTCF-binding sites at promoter regions and insulator sites associated with enhanced transcription of nearby genes. In addition, we use CTCF-associated multi-omics data to establish a combinatorial data analysis pipeline to discover CTCF co-regulatory partners. We successfully identify 40 candidates, including multiple established partners. Interestingly, many CTCF co-regulators that have alterations of their respective downstream gene expression do not show changes of their own expression levels across the multi-omics measurements upon acute CTCF loss, highlighting the strength of our system to discover hidden co-regulatory partners associated with CTCF-mediated transcription. CONCLUSIONS: This study highlights that CTCF loss rewires genome-wide chromatin accessibility, which plays a critical role in transcriptional regulation.

CovidExpress: an interactive portal for intuitive investigation on SARS-CoV-2 related transcriptomes.

Infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in humans could cause coronavirus disease 2019 (COVID-19). Since its first discovery in Dec 2019, SARS-CoV-2 has become a global pandemic and caused 3.3 million direct/indirect deaths (2021 May). Amongst the scientific community's response to COVID-19, data sharing has emerged as an essential aspect of the combat against SARS-CoV-2. Despite the ever-growing studies about SARS-CoV-2 and COVID-19, to date, only a few databases were curated to enable access to gene expression data. Furthermore, these databases curated only a small set of data and do not provide easy access for investigators without computational skills to perform analyses. To fill this gap and advance open-access to the growing gene expression data on this deadly virus, we collected about 1,500 human bulk RNA-seq datasets from publicly available resources, developed a database and visualization tool, named CovidExpress (https://stjudecab.github.io/covidexpress). This open access database will allow research investigators to examine the gene expression in various tissues, cell lines, and their response to SARS-CoV-2 under different experimental conditions, accelerating the understanding of the etiology of this disease to inform the drug and vaccine development. Our integrative analysis of this big dataset highlights a set of commonly regulated genes in SARS-CoV-2 infected lung and Rhinovirus infected nasal tissues, including OASL that were under-studied in COVID-19 related reports. Our results also suggested a potential FURIN positive feedback loop that might explain the evolutional advantage of SARS-CoV-2.

Distinct dynamics and functions of H2AK119ub1 and H3K27me3 in mouse preimplantation embryos.

Polycomb repressive complexes 1 and 2 (PRC1/2) maintain transcriptional silencing of developmental genes largely by catalyzing the formation of mono-ubiquitinated histone H2A at lysine 119 (H2AK119ub1) and trimethylated histone H3 at lysine 27 (H3K27me3), respectively. How Polycomb domains are reprogrammed during mammalian preimplantation development remains largely unclear. Here we show that, although H2AK119ub1 and H3K27me3 are highly colocalized in gametes, they undergo differential reprogramming dynamics following fertilization. H3K27me3 maintains thousands of maternally biased domains until the blastocyst stage, whereas maternally biased H2AK119ub1 distribution in zygotes is largely equalized at the two-cell stage. Notably, while maternal PRC2 depletion has a limited effect on global H2AK119ub1 in early embryos, it disrupts allelic H2AK119ub1 at H3K27me3 imprinting loci including Xist. By contrast, acute H2AK119ub1 depletion in zygotes does not affect H3K27me3 imprinting maintenance, at least by the four-cell stage. Importantly, loss of H2AK119ub1, but not H3K27me3, causes premature activation of developmental genes during zygotic genome activation (ZGA) and subsequent embryonic arrest. Thus, our study reveals distinct dynamics and functions of H3K27me3 and H2AK119ub1 in mouse preimplantation embryos.

A transcriptional roadmap for 2C-like-to-pluripotent state transition.

In mouse embryonic stem cell (ESC), a small cell population displays totipotent features by expressing a set of genes that are transiently active in 2-cell-stage embryos. These 2-cell-like (2C-like) cells spontaneously transit back into the pluripotent state. We previously dissected the transcriptional dynamics of the transition from pluripotency to the totipotent 2C-like state and identified factors that modulate the process. However, how 2C-like cells transit back into the pluripotent state remains largely unknown. In this study, we analyzed the transcriptional dynamics from the 2C-like state to pluripotent ESCs and identified an intermediate state. The intermediate state characterized by two-wave step up-regulation of pluripotent genes is different from the one observed during the 2C-like entry transition. Nonsense-mediated Dux mRNA decay plays an important role in the 2C-like state exit. Thus, our study not only provides a transcriptional roadmap for 2C-like-to-pluripotent state transition but also reveals a key molecular event driving the transition.

Cell type-specific transcriptional programs in mouse prefrontal cortex during adolescence and addiction.

Coordinated activity-induced transcriptional changes across multiple neuron subtypes of the prefrontal cortex (PFC) play a pivotal role in encoding and regulating major cognitive behaviors. Yet, the specific transcriptional programs in each neuron subtype remain unknown. Using single-cell RNA sequencing (scRNA-seq), here we comprehensively classify all unique cell subtypes in the PFC. We analyze transcriptional dynamics of each cell subtype under a naturally adaptive and an induced condition. Adaptive changes during adolescence (between P21 and P60), a highly dynamic phase of postnatal neuroplasticity, profoundly impacted transcription in each neuron subtype, including cell type-specific regulation of genes implicated in major neuropsychiatric disorders. On the other hand, an induced plasticity evoked by chronic cocaine addiction resulted in progressive transcriptional changes in multiple neuron subtypes and became most pronounced upon prolonged drug withdrawal. Our findings lay a foundation for understanding cell type-specific postnatal transcriptional dynamics under normal PFC function and in neuropsychiatric disease states.

Myc and Dnmt1 impede the pluripotent to totipotent state transition in embryonic stem cells.

Totipotency refers to the ability of a cell to generate all of the cell types of an organism. Unlike pluripotency, the establishment of totipotency is poorly understood. In mouse embryonic stem cells, Dux drives a small percentage of cells into a totipotent state by expressing 2-cell-embryo-specific transcripts. To understand how this transition takes place, we performed single-cell RNA-seq, which revealed a two-step transcriptional reprogramming process characterized by downregulation of pluripotent genes in the first step and upregulation of the 2-cell-embryo-specific elements in the second step. To identify factors controlling the transition, we performed a CRISPR-Cas9-mediated screen, which revealed Myc and Dnmt1 as two factors preventing the transition. Mechanistic studies demonstrate that Myc prevents downregulation of pluripotent genes in the first step, while Dnmt1 impedes 2-cell-embryo-specific gene activation in the second step. Collectively, the findings of our study reveal insights into the establishment and regulation of the totipotent state in mouse embryonic stem cells.

Reprogramming of Chromatin Accessibility in Somatic Cell Nuclear Transfer Is DNA Replication Independent.

Mammalian oocytes have the ability to reset the transcriptional program of differentiated somatic cells into that of totipotent embryos through somatic cell nuclear transfer (SCNT). However, the mechanisms underlying SCNT-mediated reprogramming are largely unknown. To understand the mechanisms governing chromatin reprogramming during SCNT, we profiled DNase I hypersensitive sites (DHSs) in donor cumulus cells and one-cell stage SCNT embryos. To our surprise, the chromatin accessibility landscape of the donor cells is drastically changed to recapitulate that of the in vitro fertilization (IVF)-derived zygotes within 12 hr. Interestingly, this DHS reprogramming takes place even in the presence of a DNA replication inhibitor, suggesting that SCNT-mediated DHS reprogramming is independent of DNA replication. Thus, this study not only reveals the rapid and drastic nature of the changes in chromatin accessibility through SCNT but also establishes a DNA replication-independent model for studying cellular reprogramming.

FIND: difFerential chromatin INteractions Detection using a spatial Poisson process.

Polymer-based simulations and experimental studies indicate the existence of a spatial dependency between the adjacent DNA fibers involved in the formation of chromatin loops. However, the existing strategies for detecting differential chromatin interactions assume that the interacting segments are spatially independent from the other segments nearby. To resolve this issue, we developed a new computational method, FIND, which considers the local spatial dependency between interacting loci. FIND uses a spatial Poisson process to detect differential chromatin interactions that show a significant difference in their interaction frequency and the interaction frequency of their neighbors. Simulation and biological data analysis show that FIND outperforms the widely used count-based methods and has a better signal-to-noise ratio.

Alterations of specific chromatin conformation affect ATRA-induced leukemia cell differentiation.

Chromatin conformation plays a key role in regulating gene expression and controlling cell differentiation. However, the whole-genome chromatin conformation changes that occur during leukemia cell differentiation are poorly understood. Here, we characterized the changes in chromatin conformation, histone states, chromatin accessibility, and gene expression using an all-trans retinoic acid (ATRA)-induced HL-60 cell differentiation model. The results showed that the boundaries of topological associated domains (TADs) were stable during differentiation; however, the chromatin conformations within several specific TADs were obviously changed. By combining H3K4me3, H3K27ac, and Hi-C signals, we annotated the differential gene-regulatory chromatin interactions upon ATRA induction. The gains and losses of the gene-regulatory chromatin interactions are significantly correlated with gene expression and chromatin accessibility. Finally, we found that the loss of GATA2 expression and DNA binding are crucial for the differentiation process, and changes in the chromatin structure around the GATA2 regulate its expression upon ATRA induction. This study provided both statistical insights and experimental details regarding the relationship between chromatin conformation changes and transcription regulation during leukemia cell differentiation, and the results suggested that the chromatin conformation is a new type of potential drug target for cancer therapy.

BL-Hi-C is an efficient and sensitive approach for capturing structural and regulatory chromatin interactions.

In human cells, DNA is hierarchically organized and assembled with histones and DNA-binding proteins in three dimensions. Chromatin interactions play important roles in genome architecture and gene regulation, including robustness in the developmental stages and flexibility during the cell cycle. Here we propose in situ Hi-C method named Bridge Linker-Hi-C (BL-Hi-C) for capturing structural and regulatory chromatin interactions by restriction enzyme targeting and two-step proximity ligation. This method improves the sensitivity and specificity of active chromatin loop detection and can reveal the regulatory enhancer-promoter architecture better than conventional methods at a lower sequencing depth and with a simpler protocol. We demonstrate its utility with two well-studied developmental loci: the beta-globin and HOXC cluster regions.

3CPET: finding co-factor complexes from ChIA-PET data using a hierarchical Dirichlet process.

Various efforts have been made to elucidate the cooperating proteins involved in maintaining chromatin interactions; however, many are still unknown. Here, we present 3CPET, a tool based on a non-parametric Bayesian approach, to infer the set of the most probable protein complexes involved in maintaining chromatin interactions and the regions that they may control, making it a valuable downstream analysis tool in chromatin conformation studies. 3CPET does so by combining data from ChIA-PET, transcription factor binding sites, and protein interactions. 3CPET results show biologically significant and accurate predictions when validated against experimental and simulation data.