Dynamic regulation of transcriptional states by chromatin and transcription factors

Pioneer factors: Emerging rules of engagement for transcription factors on chromatinized DNA

Pioneering transcription factors (TFs) can drive cell fate changes by binding their DNA motifs in a repressive chromatin environment. Recent structures illustrate emerging rules for nucleosome engagement: TFs distort the nucleosomal DNA to gain access or employ alternative DNA-binding modes with smaller footprints, they preferentially access solvent-exposed motifs near the entry/exit sites, and frequently interact with histones. The extent of TF-histone interactions, in turn, depends on the motif location on the nucleosome, the type of DNA-binding fold, and adjacent domains present. TF-histone interactions can phase TF motifs relative to nucleosomes, and we discuss how these complex and surprisingly diverse interactions between nucleosomes and TFs contribute to function.

Analysis of a detailed multi-stage model of stochastic gene expression using queueing theory and model reduction

We introduce a biologically detailed, stochastic model of gene expression describing the multiple rate-limiting steps of transcription, nuclear pre-mRNA processing, nuclear mRNA export, cytoplasmic mRNA degradation and translation of mRNA into protein. The processes in sub-cellular compartments are described by an arbitrary number of processing stages, thus accounting for a significantly finer molecular description of gene expression than conventional models such as the telegraph, two-stage and three-stage models of gene expression. We use two distinct tools, queueing theory and model reduction using the slow-scale linear-noise approximation, to derive exact or approximate analytic expressions for the moments or distributions of nuclear mRNA, cytoplasmic mRNA and protein fluctuations, as well as lower bounds for their Fano factors in steady-state conditions. We use these to study the phase diagram of the stochastic model; in particular we derive parametric conditions determining three types of transitions in the properties of mRNA fluctuations: from sub-Poissonian to super-Poissonian noise, from high noise in the nucleus to high noise in the cytoplasm, and from a monotonic increase to a monotonic decrease of the Fano factor with the number of processing stages. In contrast, protein fluctuations are always super-Poissonian and show weak dependence on the number of mRNA processing stages. Our results delineate the region of parameter space where conventional models give qualitatively incorrect results and provide insight into how the number of processing stages, e.g. the number of rate-limiting steps in initiation, splicing and mRNA degradation, shape stochastic gene expression by modulation of molecular memory.

Transcriptional synergy in human aortic endothelial cells is vulnerable to combination p300/CBP and BET bromodomain inhibition

Combinatorial signaling by proinflammatory cytokines synergizes to exacerbate toxicity to cells and tissue injury during acute infections. To explore synergism at the gene-regulatory level, we investigated the dynamics of transcription and chromatin signaling in response to dual cytokines by integrating nascent RNA imaging mass spectrometry, RNA sequencing, amplification-independent mRNA quantification, assay for transposase-accessible chromatin using sequencing (ATAC-seq), and transcription factor profiling. Costimulation with interferon-gamma (IFNγ) and tumor necrosis factor alpha (TNFα) synergistically induced a small subset of genes, including the chemokines CXCL9, -10, and -11. Gene induction coincided with increased chromatin accessibility at non-coding regions enriched for p65 and STAT1 binding sites. To discover coactivator dependencies, we conducted a targeted chemogenomic screen of transcriptional inhibitors followed by modeling of inhibitor dose-response curves. These results identified high efficacy of either p300/CREB-binding protein (CBP) or bromodomain and extra-terminal (BET) bromodomain inhibitors to disrupt induction of synergy genes. Combination p300/CBP and BET bromodomain inhibition at half-maximal inhibitory concentrations (subIC50) synergistically abrogated IFNγ/TNFα-induced chemokine gene and protein levels.

Epigenetic aging in major depressive disorder: Clocks, mechanisms and therapeutic perspectives

Depression, a chronic mental disorder characterized by persistent sadness, loss of interest, and difficulty in daily tasks, impacts millions globally with varying treatment options. Antidepressants, despite their long half-life and minimal effectiveness, leave half of patients undertreated, highlighting the need for new therapies to enhance well-being. Epigenetics, which studies genetic changes in gene expression or cellular phenotype without altering the underlying Deoxyribonucleic Acid (DNA) sequence, is explored in this article. This article delves into the intricate relationship between epigenetic mechanisms and depression, shedding light on how environmental stressors, early-life adversity, and genetic predispositions shape gene expression patterns associated with depression. We have also discussed Histone Deacetylase (HDAC) inhibitors, which enhance cognitive function and mood regulation in depression. Non-coding RNAs, (ncRNAs) such as Long Non-Coding RNAs (lncRNAs) and micro RNA (miRNAs), are highlighted as potential biomarkers for detecting and monitoring major depressive disorder (MDD). This article also emphasizes the reversible nature of epigenetic modifications and their influence on neuronal growth processes, underscoring the dynamic interplay between genetics, environment, and epigenetics in depression development. It explores the therapeutic potential of targeting epigenetic pathways in treating clinical depression. Additionally, it examines clinical findings related to epigenetic clocks and their role in studying depression and biological aging.

Oncogenes and tumor suppressor genes: functions and roles in cancers

Cancer, being the most formidable ailment, has had a profound impact on the human health. The disease is primarily associated with genetic mutations that impact oncogenes and tumor suppressor genes (TSGs). Recently, growing evidence have shown that X-linked TSGs have specific role in cancer progression and metastasis as well. Interestingly, our genome harbors around substantial portion of genes that function as tumor suppressors, and the X chromosome alone harbors a considerable number of TSGs. The scenario becomes even more compelling as X-linked TSGs are adaptive to key epigenetic processes such as X chromosome inactivation. Therefore, delineating the new paradigm related to X-linked TSGs, for instance, their crosstalk with autosome and involvement in cancer initiation, progression, and metastasis becomes utmost importance. Considering this, herein, we present a comprehensive discussion of X-linked TSG dysregulation in various cancers as a consequence of genetic variations and epigenetic alterations. In addition, the dynamic role of X-linked TSGs in sex chromosome-autosome crosstalk in cancer genome remodeling is being explored thoroughly. Besides, the functional roles of ncRNAs, role of X-linked TSG in immunomodulation and in gender-based cancer disparities has also been highlighted. Overall, the focal idea of the present article is to recapitulate the findings on X-linked TSG regulation in the cancer landscape and to redefine their role toward improving cancer treatment strategies.

Genome-scale chromatin binding dynamics of RNA Polymerase II general transcription machinery components

A great deal of work has revealed, in structural detail, the components of the preinitiation complex (PIC) machinery required for initiation of mRNA gene transcription by RNA polymerase II (Pol II). However, less-well understood are the in vivo PIC assembly pathways and their kinetics, an understanding of which is vital for determining how rates of in vivo RNA synthesis are established. We used competition ChIP in budding yeast to obtain genome-scale estimates of the residence times for five general transcription factors (GTFs): TBP, TFIIA, TFIIB, TFIIE and TFIIF. While many GTF-chromatin interactions were short-lived ( < 1 min), there were numerous interactions with residence times in the range of several minutes. Sets of genes with a shared function also shared similar patterns of GTF kinetic behavior. TFIIE, a GTF that enters the PIC late in the assembly process, had residence times correlated with RNA synthesis rates. The datasets and results reported here provide kinetic information for most of the Pol II-driven genes in this organism, offering a rich resource for exploring the mechanistic relationships between PIC assembly, gene regulation, and transcription.

Emerging Insights into the Role of BDNF on Health and Disease in Periphery

Brain-derived neurotrophic factor (BDNF) is a growth factor that promotes the survival and growth of developing neurons. It also enhances circuit formation to synaptic transmission for mature neurons in the brain. However, reduced BDNF expression and single nucleotide polymorphisms (SNP) are reported to be associated with functional deficit and disease development in the brain, suggesting that BDNF is a crucial molecule for brain health. Interestingly, BDNF is also expressed in the hypothalamus in appetite and energy metabolism. Previous reports demonstrated that BDNF knockout mice exhibited overeating and obesity phenotypes remarkably. Therefore, we could raise a hypothesis that the loss of function of BDNF may be associated with metabolic syndrome and peripheral diseases. In this review, we describe our recent finding that BDNF knockout mice develop metabolic dysfunction-associated steatohepatitis and recent reports demonstrating the role of one of the BDNF receptors, TrkB-T1, in some peripheral organ functions and diseases, and would provide an insight into the role of BDNF beyond the brain.

Real-time single-molecule imaging of transcriptional regulatory networks in living cells

Gene regulatory networks drive the specific transcriptional programmes responsible for the diversification of cell types during the development of multicellular organisms. Although our knowledge of the genes involved in these dynamic networks has expanded rapidly, our understanding of how transcription is spatiotemporally regulated at the molecular level over a wide range of timescales in the small volume of the nucleus remains limited. Over the past few decades, advances in the field of single-molecule fluorescence imaging have enabled real-time behaviours of individual transcriptional components to be measured in living cells and organisms. These efforts are now shedding light on the dynamic mechanisms of transcription, revealing not only the temporal rules but also the spatial coordination of underlying molecular interactions during various biological events.

Synthetic enhancers including TFREs improve transgene expression in CHO cells

The human cytomegalovirus major immediate early gene (CMV) promoter is currently the most preferred promoter for recombinant therapeutic proteins (RTPs) production in CHO cells. To enhance the production of RTPs, five synthetic enhancers including multiple transcription factor regulatory elements (TFREs) were evaluated to enhance recombinant protein level in transient and stably transfected CHO cells. Compared with the control, four elements can enhance the report genes expression under both two transfected states. Further, the function of these four enhancers on human serum albumin (HSA) were investigated. We found that the transient expression can increase by up to 1.5 times, and the stably expression can maximum increase by up to 2.14 times. The enhancement of transgene expression was caused by the boost of their corresponding mRNA levels. Transcriptomics analysis was performed and found that transcriptional activation and cell cycle regulation genes were involved. In conclusion, optimization of enhancers in the CMV promoter could increase the production yield of transgene in transfected CHO cells, which has significance for developing high-yield CHO cell expression system.

KIT mutations and expression: current knowledge and new insights for overcoming IM resistance in GIST

Gastrointestinal stromal tumor (GIST) is the most common sarcoma located in gastrointestinal tract and derived from the interstitial cell of Cajal (ICC) lineage. Both ICC and GIST cells highly rely on KIT signal pathway. Clinically, about 80-90% of treatment-naive GIST patients harbor primary KIT mutations, and special KIT-targeted TKI, imatinib (IM) showing dramatic efficacy but resistance invariably occur, 90% of them was due to the second resistance mutations emerging within the KIT gene. Although there are multiple variants of KIT mutant which did not show complete uniform biologic characteristics, most of them have high KIT expression level. Notably, the high expression level of KIT gene is not correlated to its gene amplification. Recently, accumulating evidences strongly indicated that the gene coding, epigenetic regulation, and pre- or post- protein translation of KIT mutants in GIST were quite different from that of wild type (WT) KIT. In this review, we elucidate the biologic mechanism of KIT variants and update the underlying mechanism of the expression of KIT gene, which are exclusively regulated in GIST, providing a promising yet evidence-based therapeutic landscape and possible target for the conquer of IM resistance. Video Abstract.

FOXA1 forms biomolecular condensates that unpack condensed chromatin to function as a pioneer factor

Eukaryotic DNA is packaged into chromatin in the nucleus, restricting the binding of transcription factors (TFs) to their target DNA sites. FOXA1 functions as a pioneer TF to bind condensed chromatin and initiate the opening of local chromatin for gene expression. However, the principles of FOXA1 recruitment and how it subsequently unpacks the condensed chromatin remain elusive. Here, we revealed that FOXA1 intrinsically forms submicron-sized condensates through its N- and C-terminal intrinsically disordered regions (IDRs). Notably, both IDRs enable FOXA1 to dissolve the condensed chromatin. In addition, the DNA-binding capacity of FOXA1 contributes to its ability to both form condensates and dissolve condensed chromatin. Further genome-wide investigation showed that IDRs enable FOXA1 to bind and unpack the condensed chromatin to regulate the proliferation and migration of breast cancer cells. This work provides a principle of how pioneer TFs function to initiate competent chromatin states using their IDRs.

Peak Scores Significantly Depend on the Relationships between Contextual Signals in ChIP-Seq Peaks

Chromatin immunoprecipitation followed by massively parallel DNA sequencing (ChIP-seq) is a central genome-wide method for in vivo analyses of DNA-protein interactions in various cellular conditions. Numerous studies have demonstrated the complex contextual organization of ChIP-seq peak sequences and the presence of binding sites for transcription factors in them. We assessed the dependence of the ChIP-seq peak score on the presence of different contextual signals in the peak sequences by analyzing these sequences from several ChIP-seq experiments using our fully enumerative GPU-based de novo motif discovery method, Argo_CUDA. Analysis revealed sets of significant IUPAC motifs corresponding to the binding sites of the target and partner transcription factors. For these ChIP-seq experiments, multiple regression models were constructed, demonstrating a significant dependence of the peak scores on the presence in the peak sequences of not only highly significant target motifs but also less significant motifs corresponding to the binding sites of the partner transcription factors. A significant correlation was shown between the presence of the target motifs FOXA2 and the partner motifs HNF4G, which found experimental confirmation in the scientific literature, demonstrating the important contribution of the partner transcription factors to the binding of the target transcription factor to DNA and, consequently, their important contribution to the peak score.

Distinct guard cell-specific remodeling of chromatin accessibility during abscisic acid- and CO2-dependent stomatal regulation

In plants, epidermal guard cells integrate and respond to numerous environmental signals to control stomatal pore apertures, thereby regulating gas exchange. Chromatin structure controls transcription factor (TF) access to the genome, but whether large-scale chromatin remodeling occurs in guard cells during stomatal movements, and in response to the hormone abscisic acid (ABA) in general, remains unknown. Here, we isolate guard cell nuclei from Arabidopsis thaliana plants to examine whether the physiological signals, ABA and CO2 (carbon dioxide), regulate guard cell chromatin during stomatal movements. Our cell type-specific analyses uncover patterns of chromatin accessibility specific to guard cells and define cis-regulatory sequences supporting guard cell-specific gene expression. We find that ABA triggers extensive and dynamic chromatin remodeling in guard cells, roots, and mesophyll cells with clear patterns of cell type specificity. DNA motif analyses uncover binding sites for distinct TFs enriched in ABA-induced and ABA-repressed chromatin. We identify the Abscisic Acid Response Element (ABRE) Binding Factor (ABF) bZIP-type TFs that are required for ABA-triggered chromatin opening in guard cells and roots and implicate the inhibition of a clade of bHLH-type TFs in controlling ABA-repressed chromatin. Moreover, we demonstrate that ABA and CO2 induce distinct programs of chromatin remodeling, whereby elevated atmospheric CO2 had only minimal impact on chromatin dynamics. We provide insight into the control of guard cell chromatin dynamics and propose that ABA-induced chromatin remodeling primes the genome for abiotic stress resistance.

Dynamic chromatin remodeling in cycling human endometrium at single-cell level

Estrogen-dependent proliferation followed by progesterone-dependent differentiation of the endometrium culminates in a short implantation window. We performed single-cell assay for transposase-accessible chromatin with sequencing on endometrial samples obtained across the menstrual cycle to investigate the regulation of temporal gene networks that control embryo implantation. We identify uniquely accessible chromatin regions in all major cellular constituents of the endometrium, delineate temporal patterns of coordinated chromatin remodeling in epithelial and stromal cells, and gain mechanistic insights into the emergence of a receptive state through integrated analysis of enriched transcription factor (TF) binding sites in dynamic chromatin regions, chromatin immunoprecipitation sequencing analyses, and gene expression data. We demonstrate that the implantation window coincides with pervasive cooption of transposable elements (TEs) into the regulatory chromatin landscape of decidualizing cells and expression of TE-derived transcripts in a spatially defined manner. Our data constitute a comprehensive map of the chromatin changes that control TF activities in a cycling endometrium at cellular resolution.

PRO-IP-seq tracks molecular modifications of engaged Pol II complexes at nucleotide resolution

RNA Polymerase II (Pol II) is a multi-subunit complex that undergoes covalent modifications as transcription proceeds through genes and enhancers. Rate-limiting steps of transcription control Pol II recruitment, site and degree of initiation, pausing duration, productive elongation, nascent transcript processing, transcription termination, and Pol II recycling. Here, we develop Precision Run-On coupled to Immuno-Precipitation sequencing (PRO-IP-seq), which double-selects nascent RNAs and transcription complexes, and track phosphorylation of Pol II C-terminal domain (CTD) at nucleotide-resolution. We uncover precise positional control of Pol II CTD phosphorylation as transcription proceeds from the initiating nucleotide (+1 nt), through early (+18 to +30 nt) and late (+31 to +60 nt) promoter-proximal pause, and into productive elongation. Pol II CTD is predominantly unphosphorylated from initiation until the early pause-region, whereas serine-2- and serine-5-phosphorylations are preferentially deposited in the later pause-region. Upon pause-release, serine-7-phosphorylation rapidly increases and dominates over the region where Pol II assembles elongation factors and accelerates to its full elongational speed. Interestingly, tracking CTD modifications upon heat-induced transcriptional reprogramming demonstrates that Pol II with phosphorylated CTD remains paused on thousands of heat-repressed genes. These results uncover dynamic Pol II regulation at rate-limiting steps of transcription and provide a nucleotide-resolution technique for tracking composition of engaged transcription complexes.

Distinct guard cell specific remodeling of chromatin accessibility during abscisic acid and CO2 dependent stomatal regulation

In plants, epidermal guard cells integrate and respond to numerous environmental signals to control stomatal pore apertures thereby regulating gas exchange. Chromatin structure controls transcription factor access to the genome, but whether large-scale chromatin remodeling occurs in guard cells during stomatal movements, and in response to the hormone abscisic acid (ABA) in general, remain unknown. Here we isolate guard cell nuclei from Arabidopsis thaliana plants to examine whether the physiological signals, ABA and CO2, regulate guard cell chromatin during stomatal movements. Our cell type specific analyses uncover patterns of chromatin accessibility specific to guard cells and define novel cis-regulatory sequences supporting guard cell specific gene expression. We find that ABA triggers extensive and dynamic chromatin remodeling in guard cells, roots, and mesophyll cells with clear patterns of cell-type specificity. DNA motif analyses uncover binding sites for distinct transcription factors enriched in ABA-induced and ABA-repressed chromatin. We identify the ABF/AREB bZIP-type transcription factors that are required for ABA-triggered chromatin opening in guard cells and implicate the inhibition of a set of bHLH-type transcription factors in controlling ABA-repressed chromatin. Moreover, we demonstrate that ABA and CO2 induce distinct programs of chromatin remodeling. We provide insight into the control of guard cell chromatin dynamics and propose that ABA-induced chromatin remodeling primes the genome for abiotic stress resistance.

Transcription Factor Dynamics: One Molecule at a Time

Cells must tightly regulate their gene expression programs and yet rapidly respond to acute biochemical and biophysical cues within their environment. This information is transmitted to the nucleus through various signaling cascades, culminating in the activation or repression of target genes. Transcription factors (TFs) are key mediators of these signals, binding to specific regulatory elements within chromatin. While live-cell imaging has conclusively proven that TF-chromatin interactions are highly dynamic, how such transient interactions can have long-term impacts on developmental trajectories and disease progression is still largely unclear. In this review, we summarize our current understanding of the dynamic nature of TF functions, starting with a historical overview of early live-cell experiments. We highlight key factors that govern TF dynamics and how TF dynamics, in turn, affect downstream transcriptional bursting. Finally, we conclude with open challenges and emerging technologies that will further our understanding of transcriptional regulation.

Epigenetic regulation in major depression and other stress-related disorders: molecular mechanisms, clinical relevance and therapeutic potential

Major depressive disorder (MDD) is a chronic, generally episodic and debilitating disease that affects an estimated 300 million people worldwide, but its pathogenesis is poorly understood. The heritability estimate of MDD is 30-40%, suggesting that genetics alone do not account for most of the risk of major depression. Another factor known to associate with MDD involves environmental stressors such as childhood adversity and recent life stress. Recent studies have emerged to show that the biological impact of environmental factors in MDD and other stress-related disorders is mediated by a variety of epigenetic modifications. These epigenetic modification alterations contribute to abnormal neuroendocrine responses, neuroplasticity impairment, neurotransmission and neuroglia dysfunction, which are involved in the pathophysiology of MDD. Furthermore, epigenetic marks have been associated with the diagnosis and treatment of MDD. The evaluation of epigenetic modifications holds promise for further understanding of the heterogeneous etiology and complex phenotypes of MDD, and may identify new therapeutic targets. Here, we review preclinical and clinical epigenetic findings, including DNA methylation, histone modification, noncoding RNA, RNA modification, and chromatin remodeling factor in MDD. In addition, we elaborate on the contribution of these epigenetic mechanisms to the pathological trait variability in depression and discuss how such mechanisms can be exploited for therapeutic purposes.

Glucocorticoid stimulation induces regionalized gene responses within topologically associating domains

Transcription-factor binding to cis-regulatory regions regulates the gene expression program of a cell, but occupancy is often a poor predictor of the gene response. Here, we show that glucocorticoid stimulation led to the reorganization of transcriptional coregulators MED1 and BRD4 within topologically associating domains (TADs), resulting in active or repressive gene environments. Indeed, we observed a bias toward the activation or repression of a TAD when their activities were defined by the number of regions gaining and losing MED1 and BRD4 following dexamethasone (Dex) stimulation. Variations in Dex-responsive genes at the RNA levels were consistent with the redistribution of MED1 and BRD4 at the associated cis-regulatory regions. Interestingly, Dex-responsive genes without the differential recruitment of MED1 and BRD4 or binding by the glucocorticoid receptor were found within TADs, which gained or lost MED1 and BRD4, suggesting a role of the surrounding environment in gene regulation. However, the amplitude of the response of Dex-regulated genes was higher when the differential recruitment of the glucocorticoid receptor and transcriptional coregulators was observed, reaffirming the role of transcription factor-driven gene regulation and attributing a lesser role to the TAD environment. These results support a model where a signal-induced transcription factor induces a regionalized effect throughout the TAD, redefining the notion of direct and indirect effects of transcription factors on target genes.

Genome-scale chromatin interaction dynamic measurements for key components of the RNA Pol II general transcription machinery

Background

A great deal of work has revealed in structural detail the components of the machinery responsible for mRNA gene transcription initiation. These include the general transcription factors (GTFs), which assemble at promoters along with RNA Polymerase II (Pol II) to form a preinitiation complex (PIC) aided by the activities of cofactors and site-specific transcription factors (TFs). However, less well understood are the in vivo PIC assembly pathways and their kinetics, an understanding of which is vital for determining on a mechanistic level how rates of in vivo RNA synthesis are established and how cofactors and TFs impact them.

Results

We used competition ChIP to obtain genome-scale estimates of the residence times for five GTFs: TBP, TFIIA, TFIIB, TFIIE and TFIIF in budding yeast. While many GTF-chromatin interactions were short-lived (< 1 min), there were numerous interactions with residence times in the several minutes range. Sets of genes with a shared function also shared similar patterns of GTF kinetic behavior. TFIIE, a GTF that enters the PIC late in the assembly process, had residence times correlated with RNA synthesis rates.

Conclusions

The datasets and results reported here provide kinetic information for most of the Pol II-driven genes in this organism and therefore offer a rich resource for exploring the mechanistic relationships between PIC assembly, gene regulation, and transcription. The relationships between gene function and GTF dynamics suggest that shared sets of TFs tune PIC assembly kinetics to ensure appropriate levels of expression.

Mega-Enhancer Bodies Organize Neuronal Long Genes in the Cerebellum

Dynamic regulation of gene expression plays a key role in establishing the diverse neuronal cell types in the brain. Recent findings in genome biology suggest that three-dimensional (3D) genome organization has important, but mechanistically poorly understood functions in gene transcription. Beyond local genomic interactions between promoters and enhancers, we find that cerebellar granule neurons undergoing differentiation in vivo exhibit striking increases in long-distance genomic interactions between transcriptionally active genomic loci, which are separated by tens of megabases within a chromosome or located on different chromosomes. Among these interactions, we identify a nuclear subcompartment enriched for near-megabase long enhancers and their associated neuronal long genes encoding synaptic or signaling proteins. Neuronal long genes are differentially recruited to this enhancer-dense subcompartment to help shape the transcriptional identities of granule neuron subtypes in the cerebellum. SPRITE analyses of higher-order genomic interactions, together with IGM-based 3D genome modeling and imaging approaches, reveal that the enhancer-dense subcompartment forms prominent nuclear structures, which we term mega-enhancer bodies. These novel nuclear bodies reside in the nuclear periphery, away from other transcriptionally active structures, including nuclear speckles located in the nuclear interior. Together, our findings define additional layers of higher-order 3D genome organization closely linked to neuronal maturation and identity in the brain.

3D chromatin reprogramming primes human memory TH2 cells for rapid recall and pathogenic dysfunction

Memory T cells provide long-lasting defense responses through their ability to rapidly reactivate, but how they efficiently "recall" an inflammatory transcriptional program remains unclear. Here, we show that human CD4⁺ memory T helper 2 (T_H2) cells carry a chromatin landscape synergistically reprogrammed at both one-dimensional (1D) and 3D levels to accommodate recall responses, which is absent in naive T cells. In memory T_H2 cells, recall genes were epigenetically primed through the maintenance of transcription-permissive chromatin at distal (super)enhancers organized in long-range 3D chromatin hubs. Precise transcriptional control of key recall genes occurred inside dedicated topologically associating domains ("memory TADs"), in which activation-associated promoter-enhancer interactions were preformed and exploited by AP-1 transcription factors to promote rapid transcriptional induction. Resting memory T_H2 cells from patients with asthma showed premature activation of primed recall circuits, linking aberrant transcriptional control of recall responses to chronic inflammation. Together, our results implicate stable multiscale reprogramming of chromatin organization as a key mechanism underlying immunological memory and dysfunction in T cells.

Dynamic switching of transcriptional regulators between two distinct low-mobility chromatin states

How chromatin dynamics relate to transcriptional activity remains poorly understood. Using single-molecule tracking, coupled with machine learning, we show that histone H2B and multiple chromatin-bound transcriptional regulators display two distinct low-mobility states. Ligand activation results in a marked increase in the propensity of steroid receptors to bind in the lowest-mobility state. Mutational analysis revealed that interactions with chromatin in the lowest-mobility state require an intact DNA binding domain and oligomerization domains. These states are not spatially separated as previously believed, but individual H2B and bound-TF molecules can dynamically switch between them on time scales of seconds. Single bound-TF molecules with different mobilities exhibit different dwell time distributions, suggesting that the mobility of TFs is intimately coupled with their binding dynamics. Together, our results identify two unique and distinct low-mobility states that appear to represent common pathways for transcription activation in mammalian cells.

Chromatin accessibility dynamics insight into crosstalk between regulatory landscapes in poplar responses to multiple treatments

Perennial trees develop and coordinate endogenous response signaling pathways, including their crosstalk and convergence, to cope with various environmental stresses which occur simultaneously in most cases. These processes are involved in gene transcriptional regulations that depend on dynamic interactions between regulatory proteins and corresponding chromatin regions, but the mechanisms remain poorly understood in trees. In this study, we detected chromatin regulatory landscapes of poplar under abscisic acid, methyl jasmonate, salicylic acid and sodium chloride (NaCl) treatment, through integrating ATAC-seq and RNA-seq data. Our results showed that the degree of chromatin accessibility for a given gene is closely related to its expression level. However, unlike the gene expression that shows treatment-specific response patterns, changes in chromatin accessibility exhibit high similarities under these treatments. We further proposed and experimentally validated that a homologous gene copy of RESPONSIVE TO DESICCATION 26 mediates the crosstalk between jasmonic acid and NaCl signaling pathways by directly regulating the stress-responsive genes and that circadian clock-related transcription factors like REVEILLE8 play a central role in response of poplar to these treatments. Overall, our study provides a chromatin insight into the molecular mechanism of transcription regulatory networks in response to different environmental stresses and raises the key roles of the circadian clock of poplar to adapt to adverse environments.

Polymer folding through active processes recreates features of genome organization

From proteins to chromosomes, polymers fold into specific conformations that control their biological function. Polymer folding has long been studied with equilibrium thermodynamics, yet intracellular organization and regulation involve energy-consuming, active processes. Signatures of activity have been measured in the context of chromatin motion, which shows spatial correlations and enhanced subdiffusion only in the presence of adenosine triphosphate. Moreover, chromatin motion varies with genomic coordinate, pointing toward a heterogeneous pattern of active processes along the sequence. How do such patterns of activity affect the conformation of a polymer such as chromatin? We address this question by combining analytical theory and simulations to study a polymer subjected to sequence-dependent correlated active forces. Our analysis shows that a local increase in activity (larger active forces) can cause the polymer backbone to bend and expand, while less active segments straighten out and condense. Our simulations further predict that modest activity differences can drive compartmentalization of the polymer consistent with the patterns observed in chromosome conformation capture experiments. Moreover, segments of the polymer that show correlated active (sub)diffusion attract each other through effective long-ranged harmonic interactions, whereas anticorrelations lead to effective repulsions. Thus, our theory offers nonequilibrium mechanisms for forming genomic compartments, which cannot be distinguished from affinity-based folding using structural data alone. As a first step toward exploring whether active mechanisms contribute to shaping genome conformations, we discuss a data-driven approach.

Inferring transcriptional bursting kinetics from single-cell snapshot data using a generalized telegraph model

Gene expression has inherent stochasticity resulting from transcription's burst manners. Single-cell snapshot data can be exploited to rigorously infer transcriptional burst kinetics, using mathematical models as blueprints. The classical telegraph model (CTM) has been widely used to explain transcriptional bursting with Markovian assumptions. However, growing evidence suggests that the gene-state dwell times are generally non-exponential, as gene-state switching is a multi-step process in organisms. Therefore, interpretable non-Markovian mathematical models and efficient statistical inference methods are urgently required in investigating transcriptional burst kinetics. We develop an interpretable and tractable model, the generalized telegraph model (GTM), to characterize transcriptional bursting that allows arbitrary dwell-time distributions, rather than exponential distributions, to be incorporated into the ON and OFF switching process. Based on the GTM, we propose an inference method for transcriptional bursting kinetics using an approximate Bayesian computation framework. This method demonstrates an efficient and scalable estimation of burst frequency and burst size on synthetic data. Further, the application of inference to genome-wide data from mouse embryonic fibroblasts reveals that GTM would estimate lower burst frequency and higher burst size than those estimated by CTM. In conclusion, the GTM and the corresponding inference method are effective tools to infer dynamic transcriptional bursting from static single-cell snapshot data.

Chromatin Accessibility and Transcriptional Landscape during Inhibition of Salmonella enterica by Lactobacillus reuteri in IPEC-J2 Cells

Lactobacillus reuteri is a probiotic with bacteriostatic effects, which can effectively inhibit the activity of pathogens. However, the molecular mechanism underlying the inhibition of pathogens by L. reuteri in intestinal cells remains unclear. Using the porcine intestinal cell line IPEC-J2 as a model, we combined RNA-seq and ATAC-seq methods to delineate the porcine genome-wide changes in biological processes and chromatin accessibility in IPEC-J2 cells stimulated by Salmonella enterica BNCC186354, as well as L. reuteri ATCC 53608. Overall, we found that many porcine transcripts were altered after S. enterica BNCC186354 treatment, while L. reuteri ATCC 53608 treatment partially restored this alteration, such as salmonella infection and PI3K/AKT and MAPK pathways. Combined analysis of these two datasets revealed that 26 genes with similar trends overlapped between gene expression and chromatin accessibility. In addition, we identified potential host functional transcription factors (TFs), such as GATA1, TAL1, TBP, RUNX1, Gmeb1, Gfi1b, RARA, and RXRG, in IPEC-J2 cells that might play a critical role and are targeted by L. reuteri ATCC 53608. Moreover, we verified that PI3K/AKT, MAPK, and apoptosis pathways are potentially regulated by S. enterica BNCC186354 but restored by L. reuteri ATCC 53608. The PI3K/AKT pathway was activated by L. reuteri ATCC 53608, thereby potentially inhibiting S. enterica BNCC186354 infection. In conclusion, our data provide new insights into the expression pattern of functional genes and the epigenetic alterations in IPEC-J2 cells underlying the bacteriostatic action of L. reuteri ATCC 53608.

Manipulation of the nucleoscaffold potentiates cellular reprogramming kinetics

Somatic cell fate is an outcome set by the activities of specific transcription factors and the chromatin landscape and is maintained by gene silencing of alternate cell fates through physical interactions with the nuclear scaffold. Here, we evaluate the role of the nuclear scaffold as a guardian of cell fate in human fibroblasts by comparing the effects of transient loss (knockdown) and mutation (progeria) of functional Lamin A/C, a core component of the nuclear scaffold. We observed that Lamin A/C deficiency or mutation disrupts nuclear morphology, heterochromatin levels, and increases access to DNA in lamina-associated domains. Changes in Lamin A/C were also found to impact the mechanical properties of the nucleus when measured by a microfluidic cellular squeezing device. We also show that transient loss of Lamin A/C accelerates the kinetics of cellular reprogramming to pluripotency through opening of previously silenced heterochromatin domains while genetic mutation of Lamin A/C into progerin induces a senescent phenotype that inhibits the induction of reprogramming genes. Our results highlight the physical role of the nuclear scaffold in safeguarding cellular fate.

Dynamic modeling of the cellular senescence gene regulatory network

Cellular senescence is a cell fate that prominently impacts physiological and pathophysiological processes. Diverse cellular stresses induce it, and dramatic gene expression changes accompany it. However, determining the interactions comprising the gene regulatory network (GRN) governing senescence remains challenging. Recent advances in signal processing techniques provide opportunities to reconstruct GRNs. Here, we describe a GRN for senescence integrating time-series transcriptome and transcription factor depletion datasets. Specifically, we infer a set of differential equations using the "Sparse Identification of Nonlinear Dynamics" (SINDy) algorithm, discriminate genes with potential hidden regulators, validate the inferred GRN for time-points not included in the training data, and comprehensively benchmark our approach. Our work is a proof of concept for a data-driven GRN reconstruction method, consolidating an iterative, powerful mathematical platform for senescence modeling that can be used to test hypotheses in silico and has the potential for future discoveries of clinical impact.

Multifaceted Transcriptional Network of Estrogen-Related Receptor Alpha in Health and Disease

Estrogen-related receptors (ERRα, β and γ in mammals) are orphan members of the nuclear receptor superfamily acting as transcription factors. ERRs are expressed in several cell types and they display various functions in normal and pathological contexts. Amongst others, they are notably involved in bone homeostasis, energy metabolism and cancer progression. In contrast to other nuclear receptors, the activities of the ERRs are apparently not controlled by a natural ligand but they rely on other means such as the availability of transcriptional co-regulators. Here we focus on ERRα and review the variety of co-regulators that have been identified by various means for this receptor and their reported target genes. ERRα cooperates with distinct co-regulators to control the expression of distinct sets of target genes. This exemplifies the combinatorial specificity of transcriptional regulation that induces discrete cellular phenotypes depending on the selected coregulator. We finally propose an integrated view of the ERRα transcriptional network.

Phosphorylation of ATF2 promotes odontoblastic differentiation via intrinsic HAT activity

Mouse dental papilla cells (mDPCs) are cranial neural crest-derived dental mesenchymal cells that give rise to dentin-secreting odontoblasts after the bell stage during odontogenesis. The odontoblastic differentiation of mDPCs is spatiotemporally regulated by transcription factors (TFs). Our previous work revealed that chromatin accessibility was correlated with the occupation of the basic leucine zipper (bZIP) TF family during odontoblastic differentiation. However, the detailed mechanism by which TFs regulate the initiation of odontoblastic differentiation remains elusive. Here, we report that phosphorylation of ATF2 (p-ATF2) is particularly increased during odontoblastic differentiation in vivo and in vitro. ATAC-seq and p-ATF2 CUT&Tag experiments further demonstrate a high correlation between p-ATF2 localization and increased chromatin accessibility of regions near mineralization-related genes. Knockdown of Atf2 inhibits the odontoblastic differentiation of mDPCs, while overexpression of p-ATF2 promotes odontoblastic differentiation. ATAC-seq after overexpression of p-ATF2 reveals that p-ATF2 increases the chromatin accessibility of regions adjacent to genes associated with matrix mineralization. Furthermore, we find that p-ATF2 physically interacts with and promotes H2BK12 acetylation. Taken together, our findings reveal a mechanism that p-ATF2 promotes odontoblastic differentiation at initiation via remodeling chromatin accessibility and emphasize the role of the phosphoswitch model of TFs in cell fate transitions.

A chromosome-scale genome assembly of Dasypyrum villosum provides insights into its application as a broad-spectrum disease resistance resource for wheat improvement

Dasypyrum villosum is one of the most valuable gene resources in wheat improvement, especially for disease resistance. The mining of favorable genes from D. villosum is frustrated by the lack of a whole genome sequence. In this study, we generated a doubled-haploid line, 91C43^DH, using microspore culture and obtained a 4.05-GB high-quality, chromosome-scale genome assembly for D. villosum. The assembly contains39 727 high-confidence genes, and 85.31% of the sequences are repetitive. Two reciprocal translocation events were detected, and 7VS-4VL is a unique translocation in D. villosum. The prolamin seed storage protein-coding genes were found to be duplicated; in particular, the genes encoding low-molecular-weight glutenin at the Glu-V3 locus were significantly expanded. RNA sequencing (RNA-seq) analysis indicated that, after Blumeria graminearum f.sp tritici (Bgt) inoculation, there were more upregulated genes involved in the pattern-triggered immunity and effector-triggered immunity defense pathways in D. villosum than in Triticum urartu. MNase hypersensitive sequencing (MH-seq) identified two Bgt-inducible MH sites (MHSs), one in the promoter and one in the 3' terminal region of the powdery mildew resistance (Pm) gene NLR1-V. Each site had two subpeaks and they were termed MHS1 (MHS1.1/1.2) and MHS2 (MHS2.1/2.2). Bgt-inducible MHS2.2 was uniquely present in D. villosum, and MHS1.1 was more inducible in D. villosum than in wheat, suggesting that MHSs may be critical for regulation of NLR1-V expression and plant defense. In summary, this study provides a valuable genome resource for functional genomics studies and wheat-D. villosum introgression breeding. The identified regulatory mechanisms may also be exploited to develop new strategies for enhancing Pm resistance by optimizing gene expression in wheat.

PRO-IP-seq Tracks Molecular Modifications of Engaged Pol II Complexes at Nucleotide Resolution

RNA Polymerase II (Pol II) is a multi-subunit complex that undergoes covalent modifications as transcription proceeds through genes and enhancers. Rate-limiting steps of transcription control Pol II recruitment, site and degree of initiation, pausing duration, productive elongation, nascent transcript processing, transcription termination, and Pol II recycling. Here, we developed Precision Run-On coupled to Immuno-Precipitation sequencing (PRO-IP-seq) and tracked phosphorylation of Pol II C-terminal domain (CTD) at nucleotide-resolution. We uncovered precise positional control of Pol II CTD phosphorylation as transcription proceeds from the initiating nucleotide, through early and late promoter-proximal pause, and into productive elongation. Pol II CTD was predominantly unphosphorylated in the early pause-region, whereas serine-2- and serine-5-phosphorylations occurred preferentially in the later pause-region. Serine-7-phosphorylation dominated after the pause-release in a region where Pol II accelerates to its full elongational speed. Interestingly, tracking transcription upon heat-induced reprogramming demonstrated that Pol II with phosphorylated CTD remains paused on heat-repressed genes.

Low HER2 expression in normal breast epithelium enables dedifferentiation and malignant transformation via chromatin opening

Overexpression of the HER2 protein in breast cancer patients is a predictor of poor prognosis and resistance to therapies. We used an inducible breast cancer transformation system that allows investigation of early molecular changes. HER2 overexpression to similar levels as those observed in a subtype of HER2-positive breast cancer patients induced transformation of MCF10A cells and resulted in gross morphological changes, increased anchorage-independent growth of cells, and altered the transcriptional programme of genes associated with oncogenic transformation. Global phosphoproteomic analysis during HER2 induction predominantly detected an increase in protein phosphorylation. Intriguingly, this correlated with chromatin opening, as measured by ATAC-seq on acini isolated from 3D cell culture. HER2 overexpression resulted in opening of many distal regulatory regions and promoted reprogramming-associated heterogeneity. We found that a subset of cells acquired a dedifferentiated breast stem-like phenotype, making them likely candidates for malignant transformation. Our data show that this population of cells, which counterintuitively enriches for relatively low HER2 protein abundance and increased chromatin accessibility, possesses transformational drive, resulting in increased anchorage-independent growth in vitro compared to cells not displaying a stem-like phenotype.

Noncoding Mutations in a Thyroid Hormone Receptor Gene That Impair Cone Photoreceptor Function

The function of a hormone receptor requires mechanisms to control precisely where, when, and at what level the receptor gene is expressed. An intriguing case concerns the selective induction of thyroid hormone receptor β2 (TRβ2), encoded by Thrb, in the pituitary and also in cone photoreceptors, in which it critically regulates expression of the opsin photopigments that mediate color vision. Here, we investigate the physiological significance of a candidate enhancer for induction of TRβ2 by mutagenesis of a conserved intron region in its natural context in the endogenous Thrb gene in mice. Mutation of e-box sites for bHLH (basic-helix-loop-helix) transcription factors preferentially impairs TRβ2 expression in cones whereas mutation of nearby sequences preferentially impairs expression in pituitary. A deletion encompassing all sites impairs expression in both tissues, indicating bifunctional activity. In cones, the e-box mutations disrupt chromatin acetylation, blunt the developmental induction of TRβ2, and ultimately impair cone opsin expression and sensitivity to longer wavelengths of light. These results demonstrate the necessity of studying an enhancer in its natural chromosomal context for defining biological relevance and reveal surprisingly critical nuances of level and timing of enhancer function. Our findings illustrate the influence of noncoding sequences over thyroid hormone functions.

Epigenetic Regulation of Profibrotic Macrophages in Systemic Sclerosis-Associated Interstitial Lung Disease

OBJECTIVE

Systemic sclerosis-associated interstitial lung disease (SSc-ILD) is the leading cause of death in patients with SSc with unclear pathogenesis and limited treatment options. Evidence strongly supports an important role for profibrotic secreted phosphoprotein 1 (SPP1)-expressing macrophages in SSc-ILD. This study was undertaken to define the transcriptome and chromatin structural changes of SPP1 SSc-ILD macrophages in order to better understand their role in promoting fibrosis and to identify transcription factors associated with open chromatin driving their altered phenotype.

METHODS

We performed single-cell RNA sequencing (scRNA-Seq) on 11 explanted SSc-ILD and healthy control lung samples, as well as single-cell assay for transposase-accessible chromatin sequencing on 5 lung samples to define altered chromatin accessibility of SPP1 macrophages. We predicted transcription factors regulating SPP1 macrophages using single-cell regulatory network inference and clustering (SCENIC) and determined transcription factor binding sites associated with global alterations in SPP1 chromatin accessibility using Signac/Seurat.

RESULTS

We identified distinct macrophage subpopulations using scRNA-Seq analysis in healthy and SSc-ILD lungs and assessed gene expression changes during the change of healthy control macrophages into SPP1 macrophages. Analysis of open chromatin validated SCENIC predictions, indicating that microphthalmia-associated transcription factor, transcription factor EB, activating transcription factor 6, sterol regulatory element binding transcription factor 1, basic helix-loop-helix family member E40, Kruppel-like factor 6, ETS variant transcription factor 5, and/or members of the activator protein 1 family of transcription factors regulate SPP1 macrophage differentiation.

CONCLUSION

Our findings shed light on the underlying changes in chromatin structure and transcription factor regulation of profibrotic SPP1 macrophages in SSc-ILD. Similar alterations in SPP1 macrophages may underpin fibrosis in other organs involved in SSc and point to novel targets for the treatment of SSc-ILD, specifically targeting profibrotic macrophages.

PPARγ lipodystrophy mutants reveal intermolecular interactions required for enhancer activation

Peroxisome proliferator-activated receptor γ (PPARγ) is the master regulator of adipocyte differentiation, and mutations that interfere with its function cause lipodystrophy. PPARγ is a highly modular protein, and structural studies indicate that PPARγ domains engage in several intra- and inter-molecular interactions. How these interactions modulate PPARγ's ability to activate target genes in a cellular context is currently poorly understood. Here we take advantage of two previously uncharacterized lipodystrophy mutations, R212Q and E379K, that are predicted to interfere with the interaction of the hinge of PPARγ with DNA and with the interaction of PPARγ ligand binding domain (LBD) with the DNA-binding domain (DBD) of the retinoid X receptor, respectively. Using biochemical and genome-wide approaches we show that these mutations impair PPARγ function on an overlapping subset of target enhancers. The hinge region-DNA interaction appears mostly important for binding and remodelling of target enhancers in inaccessible chromatin, whereas the PPARγ-LBD:RXR-DBD interface stabilizes the PPARγ:RXR:DNA ternary complex. Our data demonstrate how in-depth analyses of lipodystrophy mutants can unravel molecular mechanisms of PPARγ function.

Chromatin accessibility-based characterisation of brain gene regulatory networks in three distinct honey bee polyphenisms

The honey bee genome has the capacity to produce three phenotypically distinct organisms (two diploid female castes: queen and worker, and a haploid male drone). Previous studies have implicated metabolic flux acting via epigenetic regulation in directing nutrition-driven phenotypic plasticity in the honey bee. However, the cis-acting DNA regulatory elements that establish tissue and polyphenism -specific epigenomes and gene expression programmes, remain unclear. Using a high resolution multiomic approach including assay for transposase-accessible chromatin by sequencing (ATAC-seq), RNA-seq and ChIP-seq, we produce the first genome-wide maps of the regulatory landscape across all three adult honey bee phenotypes identifying > 5000 regulatory regions in queen, 7500 in worker and 6500 in drone, with the vast majority of these sites located within intronic regions. These regions are defined by positive enrichment of H3K27ac and depletion of H3K4me3 and show a positive correlation with gene expression. Using ATAC-seq footprinting we determine queen, worker and drone -specific transcription factor occupancy and uncover novel phenotype-specific regulatory networks identifying two key nuclear receptors that have previously been implicated in caste-determination and adult behavioural maturation in honey bees; ecdysone receptor and ultraspiracle. Collectively, this study provides novel insights into key gene regulatory networks that are associated with these distinct polyphenisms in the honey bee.

Global changes in chromatin accessibility and transcription in growth hormone-secreting pituitary adenoma

PURPOSE

Growth hormone-secreting pituitary adenoma (GHPA) is an insidious disease with persistent hypersecretion of growth hormone and insulin-like growth factor 1, causing increased morbidity and mortality. Previous studies have investigated the transcription of GHPA. However, the gene regulatory landscape has not been fully characterized. The objective of our study was to unravel the changes in chromatin accessibility and transcription in GHPA.

METHODS

Six patients diagnosed with GHPA in the Department of Neurosurgery at Huashan Hospital were enrolled in our study. Primary pituitary adenoma tissues and adjacent normal pituitary specimens with no morphologic abnormalities from these six patients were obtained at surgery. RNA sequencing (RNA-seq) and assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) were applied to investigate the underlying relationship between gene expression and chromatin accessibility changes in GHPA.

RESULTS

Totally, 1528 differential expression genes (DEGs) were identified by transcriptomics analyses, including 725 up-regulated and 803 down-regulated. Further, we obtained 64 significantly DEGs including 10 DEGs were elevated and 54 DEGs were negligibly expressed in tumors tissues. The up-regulated DEGs were mainly involved in terms related to synapse formation, nervous system development and secretory pathway. In parallel, 3916 increased and 2895 decreased chromatin-accessible regions were mapped by ATAC-seq. Additionally, the chromatin accessible changes were frequently located adjacent to transcription factor CTCF and Rfx2 binding site.

CONCLUSIONS

Our results are the first to demonstrate the landscape of chromatin accessibility in GHPA, which may contribute to illustrate the underlying transcriptional regulation mechanism of this disease.

Integrated multi-omics approach revealed cellular senescence landscape

Cellular senescence is a complex multifactorial biological phenomenon that plays essential roles in aging, and aging-related diseases. During this process, the senescent cells undergo gene expression altering and chromatin structure remodeling. However, studies on the epigenetic landscape of senescence using integrated multi-omics approaches are limited. In this research, we performed ATAC-seq, RNA-seq and ChIP-seq on different senescent types to reveal the landscape of senescence and identify the prime regulatory elements. We also obtained 34 key genes and deduced that NAT1, PBX1 and RRM2, which interacted with each other, could be the potential markers of aging and aging-related diseases. In summary, our work provides the landscape to study accessibility dynamics and transcriptional regulations in cellular senescence. The application of this technique in different types of senescence allows us to identify the regulatory elements responsible for the substantial regulation of transcription, providing the insights into molecular mechanisms of senescence.

Three-dimensional genome organization in immune cell fate and function

Immune cell development and activation demand the precise and coordinated control of transcriptional programmes. Three-dimensional (3D) organization of the genome has emerged as an important regulator of chromatin state, transcriptional activity and cell identity by facilitating or impeding long-range genomic interactions among regulatory elements and genes. Chromatin folding thus enables cell type-specific and stimulus-specific transcriptional responses to extracellular signals, which are essential for the control of immune cell fate, for inflammatory responses and for generating a diverse repertoire of antigen receptor specificities. Here, we review recent findings connecting 3D genome organization to the control of immune cell differentiation and function, and discuss how alterations in genome folding may lead to immune dysfunction and malignancy.

3D chromatin maps of the human pancreas reveal lineage-specific regulatory architecture of T2D risk

Three-dimensional (3D) chromatin organization maps help dissect cell-type-specific gene regulatory programs. Furthermore, 3D chromatin maps contribute to elucidating the pathogenesis of complex genetic diseases by connecting distal regulatory regions and genetic risk variants to their respective target genes. To understand the cell-type-specific regulatory architecture of diabetes risk, we generated transcriptomic and 3D epigenomic profiles of human pancreatic acinar, alpha, and beta cells using single-cell RNA-seq, single-cell ATAC-seq, and high-resolution Hi-C of sorted cells. Comparisons of these profiles revealed differential A/B (open/closed) chromatin compartmentalization, chromatin looping, and transcriptional factor-mediated control of cell-type-specific gene regulatory programs. We identified a total of 4,750 putative causal-variant-to-target-gene pairs at 194 type 2 diabetes GWAS signals using pancreatic 3D chromatin maps. We found that the connections between candidate causal variants and their putative target effector genes are cell-type stratified and emphasize previously underappreciated roles for alpha and acinar cells in diabetes pathogenesis.

Optimization of synthetic molecular reporters for a mesenchymal glioblastoma transcriptional program by integer programing

MOTIVATION

A recent approach to perform genetic tracing of complex biological problems involves the generation of synthetic deoxyribonucleic acid (DNA) probes that specifically mark cells with a phenotype of interest. These synthetic locus control regions (sLCRs), in turn, drive the expression of a reporter gene, such as fluorescent protein. To build functional and specific sLCRs, it is critical to accurately select multiple bona fide cis-regulatory elements from the target cell phenotype cistrome. This selection occurs by maximizing the number and diversity of transcription factors (TFs) within the sLCR, yet the size of the final sLCR should remain limited.

RESULTS

In this work, we discuss how optimization, in particular integer programing, can be used to systematically address the construction of a specific sLCR and optimize pre-defined properties of the sLCR. Our presented instance of a linear optimization problem maximizes the activation potential of the sLCR such that its size is limited to a pre-defined length and a minimum number of all TFs deemed sufficiently characteristic for the phenotype of interest is covered. We generated an sLCR to trace the mesenchymal glioblastoma program in patients by solving our corresponding linear program with the software optimizer Gurobi. Considering the binding strength of transcription factor binding sites (TFBSs) with their TFs as a proxy for activation potential, the optimized sLCR scores similarly to an sLCR experimentally validated in vivo, and is smaller in size while having the same coverage of TFBSs.

AVAILABILITY AND IMPLEMENTATION

We provide a Python implementation of the presented framework in the Supplementary Material with which an optimal selection of cis-regulatory elements can be calculated once the target set of TFs and their binding strength with their TFBSs is known.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

DNAffinity: a machine-learning approach to predict DNA binding affinities of transcription factors

We present a physics-based machine learning approach to predict in vitro transcription factor binding affinities from structural and mechanical DNA properties directly derived from atomistic molecular dynamics simulations. The method is able to predict affinities obtained with techniques as different as uPBM, gcPBM and HT-SELEX with an excellent performance, much better than existing algorithms. Due to its nature, the method can be extended to epigenetic variants, mismatches, mutations, or any non-coding nucleobases. When complemented with chromatin structure information, our in vitro trained method provides also good estimates of in vivo binding sites in yeast.

Transcription factor-nucleosome dynamics from plasma cfDNA identifies ER-driven states in breast cancer

Genome-wide binding profiles of estrogen receptor (ER) and FOXA1 reflect cancer state in ER+ breast cancer. However, routine profiling of tumor transcription factor (TF) binding is impractical in the clinic. Here, we show that plasma cell-free DNA (cfDNA) contains high-resolution ER and FOXA1 tumor binding profiles for breast cancer. Enrichment of TF footprints in plasma reflects the binding strength of the TF in originating tissue. We defined pure in vivo tumor TF signatures in plasma using ER+ breast cancer xenografts, which can distinguish xenografts with distinct ER states. Furthermore, state-specific ER-binding signatures can partition human breast tumors into groups with significantly different ER expression and mortality. Last, TF footprints in human plasma samples can identify the presence of ER+ breast cancer. Thus, plasma TF footprints enable minimally invasive mapping of the regulatory landscape of breast cancer in humans and open vast possibilities for clinical applications across multiple tumor types.

Transcriptome Analysis and Intraspecific Variation in Spanish Fir (Abies pinsapo Boiss.)

Spanish fir (Abies pinsapo Boiss.) is an endemic, endangered tree that has been scarcely investigated at the molecular level. In this work, the transcriptome of Spanish fir was assembled, providing a large catalog of expressed genes (22,769), within which a high proportion were full-length transcripts (12,545). This resource is valuable for functional genomics studies and genome annotation in this relict conifer species. Two intraspecific variations of A. pinsapo can be found within its largest population at the Sierra de las Nieves National Park: one with standard green needles and another with bluish-green needles. To elucidate the causes of both phenotypes, we studied different physiological and molecular markers and transcriptome profiles in the needles. "Green" trees showed higher electron transport efficiency and enhanced levels of chlorophyll, protein, and total nitrogen in the needles. In contrast, needles from "bluish" trees exhibited higher contents of carotenoids and cellulose. These results agreed with the differential transcriptomic profiles, suggesting an imbalance in the nitrogen status of "bluish" trees. Additionally, gene expression analyses suggested that these differences could be associated with different epigenomic profiles. Taken together, the reported data provide new transcriptome resources and a better understanding of the natural variation in this tree species, which can help improve guidelines for its conservation and the implementation of adaptive management strategies under climatic change.

NICD3 regulates the expression of MUC5AC and MUC2 by recruiting SMARCA4 and is involved in the differentiation of mucinous colorectal adenocarcinoma

Adenocarcinoma is the most prevalent histological subtype of colorectal cancer (CRC), with mucinous colorectal adenocarcinoma (MCA) being a unique form. Although the mucinous subtype is known to elicit a worse response to chemotherapy and immunotherapy than the nonmucinous subtype, its pathogenesis remains poorly understood. Neurogenic locus notch homolog protein 3 (NOTCH3), a member of the NOTCH subfamilies, is highly expressed in CRC. In the past three decades, many studies have been performed evaluating the biological role of NOTCH3 in CRC. However, the precise activities of NOTCH3 in MCA, as well as the mechanisms involved in its transcriptional control, are yet to be elucidated. Our finding showed that the critical transcriptional regulatory factor transcription activator BRG1 (SMARCA4) directly binds to the intracellular domain of NOTCH3 to control transcriptional regulation. Moreover, RNA‐sequencing results indicated a common targeting effect on the transcriptional activity of mucin‐5AC (MUC5AC) and mucin‐2 (MUC2) in CRC cells by NOTCH3 and SMARCA4. Furthermore, NOTCH3 was found to control the expressions of MUC5AC and MUC2 in a SMARCA4‐dependent manner. MUC5AC and MUC2, which encode two secreted mucins, are located on chromosome 11p15.5, and are linked to the development of MCA. This finding suggests that the interaction between NOTCH3 and SMARCA4 may be involved in MCA differentiation by jointly targeting MUC5AC and MUC2. Patients with MCA are often treated in accordance with CRC guidelines. Determining the relationship between NOTCH3 and SMARCA4 by demonstrating their interactions in the pathophysiology of MCA could provide novel therapeutic targets and help identify potential prognostic markers for MCA.

Generating specificity in genome regulation through transcription factor sensitivity to chromatin

Cell type-specific gene expression relies on transcription factors (TFs) binding DNA sequence motifs embedded in chromatin. Understanding how motifs are accessed in chromatin is crucial to comprehend differential transcriptional responses and the phenotypic impact of sequence variation. Chromatin obstacles to TF binding range from DNA methylation to restriction of DNA access by nucleosomes depending on their position, composition and modification. In vivo and in vitro approaches now enable the study of TF binding in chromatin at unprecedented resolution. Emerging insights suggest that TFs vary in their ability to navigate chromatin states. However, it remains challenging to link binding and transcriptional outcomes to molecular characteristics of TFs or the local chromatin substrate. Here, we discuss our current understanding of how TFs access DNA in chromatin and novel techniques and directions towards a better understanding of this critical step in genome regulation.

MOZ and Menin-MLL Complexes Are Complementary Regulators of Chromatin Association and Transcriptional Output in Gastrointestinal Stromal Tumor

Gastrointestinal stromal tumor (GIST) is commonly characterized by activating mutations in the receptor tyrosine kinase KIT. Tyrosine kinase inhibitors are the only approved therapy for GIST, and complementary treatment strategies are urgently needed. As GIST lacks oncogene amplification and relies upon an established network of transcription factors, we hypothesized that unique chromatin-modifying enzymes are essential in orchestrating the GIST epigenome. We identified through genome-scale CRISPR screening that MOZ and Menin-MLL chromatin regulatory complexes are cooperative and unique dependencies in GIST. These complexes were enriched at GIST-relevant genes and regulated their transcription. Inhibition of MOZ and Menin-MLL complexes decreased GIST cell proliferation by disrupting interactions with transcriptional/chromatin regulators, such as DOT1L. MOZ and Menin inhibition caused significant reductions in tumor burden in vivo, with superior effects observed with combined Menin and KIT inhibition. These results define unique chromatin regulatory dependencies in GIST and identify potential therapeutic strategies for clinical application.

SIGNIFICANCE

Although many malignancies rely on oncogene amplification, GIST instead depends upon epigenetic regulation of KIT and other essential genes. Utilizing genome-scale CRISPR dependency screens, we identified complementary chromatin-modifying complexes essential to GIST and characterize the consequences of their disruption, elucidating a novel therapeutic approach to this disease. This article is highlighted in the In This Issue feature, p. 1599.