Pioneer transcription factors: establishing competence for gene expression

The importance of DNA sequence for nucleosome positioning in transcriptional regulation

Nucleosome positioning is a key factor for transcriptional regulation. Nucleosomes regulate the dynamic accessibility of chromatin and interact with the transcription machinery at every stage. Influences to steer nucleosome positioning are diverse, and the according importance of the DNA sequence in contrast to active chromatin remodeling has been the subject of long discussion. In this study, we evaluate the functional role of DNA sequence for all major elements along the process of transcription. We developed a random forest classifier based on local DNA structure that assesses the sequence-intrinsic support for nucleosome positioning. On this basis, we created a simple data resource that we applied genome-wide to the human genome. In our comprehensive analysis, we found a special role of DNA in mediating the competition of nucleosomes with cis-regulatory elements, in enabling steady transcription, for positioning of stable nucleosomes in exons, and for repelling nucleosomes during transcription termination. In contrast, we relate these findings to concurrent processes that generate strongly positioned nucleosomes in vivo that are not mediated by sequence, such as energy-dependent remodeling of chromatin.

Pioneer Transcription Factors: The First Domino in Zygotic Genome Activation

Zygotic genome activation (ZGA) is a pivotal event in mammalian embryogenesis, marking the transition from maternal to zygotic control of development. During the ZGA process that is characterized by the intricate cascade of gene expression, who tipped the first domino in a meticulously arranged sequence is a subject of paramount interest. Recently, Dux, Obox and Nr5a2 were identified as pioneer transcription factors that reside at the top of transcriptional hierarchy. Through co-option of retrotransposon elements as hubs for transcriptional activation, these pioneer transcription factors rewire the gene regulatory network, thus initiating ZGA. In this review, we provide a snapshot of the mechanisms underlying the functions of these pioneer transcription factors. We propose that ZGA is the starting point where the embryo's own genome begins to influence development trajectory, therefore in-depth dissecting the functions of pioneer transcription factors during ZGA will form a cornerstone of our understanding for early embryonic development, which will pave the way for advancing our grasp of mammalian developmental biology and optimizing in vitro production (IVP) techniques.

The C. elegans Myc-family of transcription factors coordinate a dynamic adaptive response to dietary restriction

Dietary restriction (DR), the process of decreasing overall food consumption over an extended period of time, has been shown to increase longevity across evolutionarily diverse species and delay the onset of age-associated diseases in humans. In Caenorhabditis elegans, the Myc-family transcription factors (TFs) MXL-2 (Mlx) and MML-1 (MondoA/ChREBP), which function as obligate heterodimers, and PHA-4 (orthologous to FOXA) are both necessary for the full physiological benefits of DR. However, the adaptive transcriptional response to DR and the role of MML-1::MXL-2 and PHA-4 remains elusive. We identified the transcriptional signature of C. elegans DR, using the eat-2 genetic model, and demonstrate broad changes in metabolic gene expression in eat-2 DR animals, which requires both mxl-2 and pha-4. While the requirement for these factors in DR gene expression overlaps, we found many of the DR genes exhibit an opposing change in relative gene expression in eat-2;mxl-2 animals compared to wild-type, which was not observed in eat-2 animals with pha-4 loss. Surprisingly, we discovered more than 2000 genes synthetically dysregulated in eat-2;mxl-2, out of which the promoters of down-regulated genes were substantially enriched for PQM-1 and ELT-1/3 GATA TF binding motifs. We further show functional deficiencies of the mxl-2 loss in DR outside of lifespan, as eat-2;mxl-2 animals exhibit substantially smaller brood sizes and lay a proportion of dead eggs, indicating that MML-1::MXL-2 has a role in maintaining the balance between resource allocation to the soma and to reproduction under conditions of chronic food scarcity. While eat-2 animals do not show a significantly different metabolic rate compared to wild-type, we also find that loss of mxl-2 in DR does not affect the rate of oxygen consumption in young animals. The gene expression signature of eat-2 mutant animals is consistent with optimization of energy utilization and resource allocation, rather than induction of canonical gene expression changes associated with acute metabolic stress, such as induction of autophagy after TORC1 inhibition. Consistently, eat-2 animals are not substantially resistant to stress, providing further support to the idea that chronic DR may benefit healthspan and lifespan through efficient use of limited resources rather than broad upregulation of stress responses, and also indicates that MML-1::MXL-2 and PHA-4 may have distinct roles in promotion of benefits in response to different pro-longevity stimuli.

Structure-based learning to predict and model protein-DNA interactions and transcription-factor co-operativity in cis-regulatory elements

Transcription factor (TF) binding is a key component of genomic regulation. There are numerous high-throughput experimental methods to characterize TF-DNA binding specificities. Their application, however, is both laborious and expensive, which makes profiling all TFs challenging. For instance, the binding preferences of ∼25% human TFs remain unknown; they neither have been determined experimentally nor inferred computationally. We introduce a structure-based learning approach to predict the binding preferences of TFs and the automated modelling of TF regulatory complexes. We show the advantage of using our approach over the classical nearest-neighbor prediction in the limits of remote homology. Starting from a TF sequence or structure, we predict binding preferences in the form of motifs that are then used to scan a DNA sequence for occurrences. The best matches are either profiled with a binding score or collected for their subsequent modeling into a higher-order regulatory complex with DNA. Co-operativity is modelled by: (i) the co-localization of TFs and (ii) the structural modeling of protein-protein interactions between TFs and with co-factors. We have applied our approach to automatically model the interferon-β enhanceosome and the pioneering complexes of OCT4, SOX2 (or SOX11) and KLF4 with a nucleosome, which are compared with the experimentally known structures.

Pink adipose tissue: A paradigm of adipose tissue plasticity

Adipose tissue is highly plastic, as illustrated mainly by the transdifferentiation of white adipocytes into beige adipocytes, depending on environmental conditions. However, during gestation and lactation in rodent, there is an amazing phenomenon of transformation of subcutaneous adipose tissue into mammary glandular tissue, known as pink adipose tissue, capable of synthesizing and secreting milk. Recent work using transgenic lineage-tracing experiments, mainly carried out in Saverio Cinti's team, has demonstrated very convincingly that this process does indeed correspond to a transdifferentiation of white adipocytes into mammary alveolar cells (pink adipocytes) during gestation and lactation. This phenomenon is reversible, since during the post-lactation phase, pink adipocytes revert to the white adipocyte phenotype. The molecular mechanisms underlying this reversible transdifferentiation remain poorly understood.

Adipo-Epithelial Transdifferentiation in In Vitro Models of the Mammary Gland

Subcutaneous adipocytes are crucial for mammary gland epithelial development during pregnancy. Our and others' previous data have suggested that adipo-epithelial transdifferentiation could play a key role in the mammary gland alveolar development. In this study, we tested whether adipo-epithelial transdifferentiation occurs in vitro. Data show that, under appropriate co-culture conditions with mammary epithelial organoids (MEOs), mature adipocytes lose their phenotype and acquire an epithelial one. Interestingly, even in the absence of MEOs, extracellular matrix and diffusible growth factors are able to promote adipo-epithelial transdifferentiation. Gene and protein expression studies indicate that transdifferentiating adipocytes exhibit some characteristics of milk-secreting alveolar glands, including significantly higher expression of milk proteins such as whey acidic protein and β-casein. Similar data were also obtained in cultured human multipotent adipose-derived stem cell adipocytes. A miRNA sequencing experiment on the supernatant highlighted mir200c, which has a well-established role in the mesenchymal-epithelial transition, as a potential player in this phenomenon. Collectively, our data show that adipo-epithelial transdifferentiation can be reproduced in in vitro models where this phenomenon can be investigated at the molecular level.

Bioinformatics analysis for constructing a cellular senescence-related age-related macular degeneration diagnostic model and identifying relevant disease subtypes to guide treatment

Age-related macular degeneration (AMD) is a condition causing progressive central vision loss. Growing evidence suggests a link between cellular senescence and AMD. However, the exact mechanism by which cellular senescence leads to AMD remains unclear. Employing machine learning, we established an AMD diagnostic model. Through unsupervised clustering, two distinct AMD subtypes were identified. GO, KEGG, and GSVA analyses explored the diverse biological functions associated with the two subtypes. By WGCNA, we constructed a coexpression network of differential genes between the subtypes, revealing the regulatory role of hub genes at the level of transcription factors and miRNAs. We identified 5 genes associated with inflammation for the construction of the AMD diagnostic model. Additionally, we observed that the level of cellular senescence and pathways related to programmed cell death (PCD), such as ferroptosis, necroptosis, and pyroptosis, exhibited higher expression levels in subtype B than A. Immune microenvironments also differed between the subtypes, indicating potentially distinct pathogenic mechanisms and therapeutic targets. In summary, by leveraging cellular senescence-associated gene expression, we developed an AMD diagnostic model. Furthermore, we identified two subtypes with varying expression patterns of senescence genes, revealing their differential roles in programmed cell death, disease progression, and immune microenvironments within AMD.

STR mutations on chromosome 15q cause thyrotropin resistance by activating a primate-specific enhancer of MIR7-2/MIR1179

Thyrotropin (TSH) is the master regulator of thyroid gland growth and function. Resistance to TSH (RTSH) describes conditions with reduced sensitivity to TSH. Dominantly inherited RTSH has been linked to a locus on chromosome 15q, but its genetic basis has remained elusive. Here we show that non-coding mutations in a (TTTG)4 short tandem repeat (STR) underlie dominantly inherited RTSH in all 82 affected participants from 12 unrelated families. The STR is contained in a primate-specific Alu retrotransposon with thyroid-specific cis-regulatory chromatin features. Fiber-seq and RNA-seq studies revealed that the mutant STR activates a thyroid-specific enhancer cluster, leading to haplotype-specific upregulation of the bicistronic MIR7-2/MIR1179 locus 35 kb downstream and overexpression of its microRNA products in the participants' thyrocytes. An imbalance in signaling pathways targeted by these micro-RNAs provides a working model for this cause of RTSH. This finding broadens our current knowledge of genetic defects altering pituitary-thyroid feedback regulation.

Contribution of the transcription factor SfGATAe to Bt Cry toxin resistance in Spodoptera frugiperda through reduction of ABCC2 expression

Insect resistance evolution poses a significant threat to the advantages of biopesticides and transgenic crops utilizing insecticidal Cry-toxins from Bacillus thuringiensis (Bt). However, there is limited research on the relationship between transcriptional regulation of specific toxin receptors in lepidopteran insects and their resistance to Bt toxins. Here, we report the positive regulatory role of the SfGATAe transcription factor on the expression of the ABCC2 gene in Spodoptera frugiperda. DNA regions in the SfABCC2 promoter that are vital for regulation by SfGATAe, utilizing DAP-seq technology and promoter deletion mapping. Through yeast one-hybrid assays, DNA pull-down experiments, and site-directed mutagenesis, we confirmed that the transcription factor SfGATAe regulates the core control site PBS2 in the ABCC2 target gene. Tissue-specific expression analysis has revealed that SfGATAe is involved in the regulation and expression of midgut cells in the fall armyworm. Silencing SfGATAe in fall armyworm larvae resulted in reduced expression of SfABCC2 and decreased sensitivity to Cry1Ac toxin. Overall, this study elucidated the regulatory mechanism of the transcription factor SfGATAe on the expression of the toxin receptor gene SfABCC2 and this transcriptional control mechanism impacts the resistance of the fall armyworm to Bt toxins.

Nucleosome-bound NR5A2 structure reveals pioneer factor mechanism by DNA minor groove anchor competition

Gene expression during natural and induced reprogramming is controlled by pioneer transcription factors that initiate transcription from closed chromatin. Nr5a2 is a key pioneer factor that regulates zygotic genome activation in totipotent embryos, pluripotency in embryonic stem cells and metabolism in adult tissues, but the mechanism of its pioneer activity remains poorly understood. Here, we present a cryo-electron microscopy structure of human NR5A2 bound to a nucleosome. The structure shows that the conserved carboxy-terminal extension (CTE) loop of the NR5A2 DNA-binding domain competes with a DNA minor groove anchor of the nucleosome and releases entry-exit site DNA. Mutational analysis showed that NR5A2 D159 of the CTE is dispensable for DNA binding but required for stable nucleosome association and persistent DNA 'unwrapping'. These findings suggest that NR5A2 belongs to an emerging class of pioneer factors that can use DNA minor groove anchor competition to destabilize nucleosomes and facilitate gene expression during reprogramming.

Individual transcription factors modulate both the micromovement of chromatin and its long-range structure

The control of eukaryotic gene expression is intimately connected to highly dynamic chromatin structures. Gene regulation relies on activator and repressor transcription factors (TFs) that induce local chromatin opening and closing. However, it is unclear how nucleus-wide chromatin organization responds dynamically to the activity of specific TFs. Here, we examined how two TFs with opposite effects on local chromatin accessibility modulate chromatin dynamics nucleus-wide. We combine high-resolution diffusion mapping and dense flow reconstruction and correlation in living cells to obtain an imaging-based, nanometer-scale analysis of local diffusion processes and long-range coordinated movements of both chromatin and TFs. We show that the expression of either an individual transcriptional activator (CDX2) or repressor (SIX6) with large numbers of binding sites increases chromatin mobility nucleus-wide, yet they induce opposite coherent chromatin motions at the micron scale. Hi-C analysis of higher-order chromatin structures shows that induction of the pioneer factor CDX2 leads both to changes in local chromatin interactions and the distribution of A and B compartments, thus relating the micromovement of chromatin with changes in compartmental structures. Given that inhibition of transcription initiation and elongation by RNA Pol II has a partial impact on the global chromatin dynamics induced by CDX2, we suggest that CDX2 overexpression alters chromatin structure dynamics both dependently and independently of transcription. Our biophysical analysis shows that sequence-specific TFs can influence chromatin structure on multiple architectural levels, arguing that local chromatin changes brought by TFs alter long-range chromatin mobility and its organization.

Induced pluripotent stem cells (iPSCs): molecular mechanisms of induction and applications

The induced pluripotent stem cell (iPSC) technology has transformed in vitro research and holds great promise to advance regenerative medicine. iPSCs have the capacity for an almost unlimited expansion, are amenable to genetic engineering, and can be differentiated into most somatic cell types. iPSCs have been widely applied to model human development and diseases, perform drug screening, and develop cell therapies. In this review, we outline key developments in the iPSC field and highlight the immense versatility of the iPSC technology for in vitro modeling and therapeutic applications. We begin by discussing the pivotal discoveries that revealed the potential of a somatic cell nucleus for reprogramming and led to successful generation of iPSCs. We consider the molecular mechanisms and dynamics of somatic cell reprogramming as well as the numerous methods available to induce pluripotency. Subsequently, we discuss various iPSC-based cellular models, from mono-cultures of a single cell type to complex three-dimensional organoids, and how these models can be applied to elucidate the mechanisms of human development and diseases. We use examples of neurological disorders, coronavirus disease 2019 (COVID-19), and cancer to highlight the diversity of disease-specific phenotypes that can be modeled using iPSC-derived cells. We also consider how iPSC-derived cellular models can be used in high-throughput drug screening and drug toxicity studies. Finally, we discuss the process of developing autologous and allogeneic iPSC-based cell therapies and their potential to alleviate human diseases.

ARID1A orchestrates SWI/SNF-mediated sequential binding of transcription factors with ARID1A loss driving pre-memory B cell fate and lymphomagenesis

ARID1A, a subunit of the canonical BAF nucleosome remodeling complex, is commonly mutated in lymphomas. We show that ARID1A orchestrates B cell fate during the germinal center (GC) response, facilitating cooperative and sequential binding of PU.1 and NF-kB at crucial genes for cytokine and CD40 signaling. The absence of ARID1A tilts GC cell fate toward immature IgM+CD80-PD-L2- memory B cells, known for their potential to re-enter new GCs. When combined with BCL2 oncogene, ARID1A haploinsufficiency hastens the progression of aggressive follicular lymphomas (FLs) in mice. Patients with FL with ARID1A-inactivating mutations preferentially display an immature memory B cell-like state with increased transformation risk to aggressive disease. These observations offer mechanistic understanding into the emergence of both indolent and aggressive ARID1A-mutant lymphomas through the formation of immature memory-like clonal precursors. Lastly, we demonstrate that ARID1A mutation induces synthetic lethality to SMARCA2/4 inhibition, paving the way for potential precision therapy for high-risk patients.

The role of histone acetylation in transcriptional regulation and seed development

Histone acetylation is highly conserved across eukaryotes and has been linked to gene activation since its discovery nearly 60 years ago. Over the past decades, histone acetylation has been evidenced to play crucial roles in plant development and response to various environmental cues. Emerging data indicate that histone acetylation is one of the defining features of "open chromatin," while the role of histone acetylation in transcription remains controversial. In this review, we briefly describe the discovery of histone acetylation, the mechanism of histone acetylation regulating transcription in yeast and mammals, and summarize the research progress of plant histone acetylation. Furthermore, we also emphasize the effect of histone acetylation on seed development and its potential use in plant breeding. A comprehensive knowledge of histone acetylation might provide new and more flexible research perspectives to enhance crop yield and stress resistance.

DNA controls the dimerization of the human FoxP1 forkhead domain

Transcription factors (TFs) regulate gene expression by binding to specific DNA sequences and gating access to genes. Even when the binding of TFs and their cofactors to DNA is reversible, indicating a reversible control of gene expression, there is little knowledge about the molecular effect DNA has on TFs. Using single-molecule multiparameter fluorescence spectroscopy, molecular dynamics simulations, and biochemical assays, we find that the monomeric form of the forkhead (FKH) domain of the human FoxP1 behaves as a disordered protein and increases its folded population when it dimerizes. Notably, DNA binding promotes a disordered FKH dimer bound to DNA, negatively controlling the stability of the dimeric FoxP1:DNA complex. The DNA-mediated reversible regulation on FKH dimers suggests that FoxP1-dependent gene suppression is unstable, and it must require the presence of other dimerization domains or cofactors to revert the negative impact exerted by the DNA.

A Foxf1-Wnt-Nr2f1 cascade promotes atrial cardiomyocyte differentiation in zebrafish

Nr2f transcription factors (TFs) are conserved regulators of vertebrate atrial cardiomyocyte (AC) differentiation. However, little is known about the mechanisms directing Nr2f expression in ACs. Here, we identified a conserved enhancer 3' to the nr2f1a locus, which we call 3'reg1-nr2f1a (3'reg1), that can promote Nr2f1a expression in ACs. Sequence analysis of the enhancer identified putative Lef/Tcf and Foxf TF binding sites. Mutation of the Lef/Tcf sites within the 3'reg1 reporter, knockdown of Tcf7l1a, and manipulation of canonical Wnt signaling support that Tcf7l1a is derepressed via Wnt signaling to activate the transgenic enhancer and promote AC differentiation. Similarly, mutation of the Foxf binding sites in the 3'reg1 reporter, coupled with gain- and loss-of-function analysis supported that Foxf1 promotes expression of the enhancer and AC differentiation. Functionally, we find that Wnt signaling acts downstream of Foxf1 to promote expression of the 3'reg1 reporter within ACs and, importantly, both Foxf1 and Wnt signaling require Nr2f1a to promote a surplus of differentiated ACs. CRISPR-mediated deletion of the endogenous 3'reg1 abrogates the ability of Foxf1 and Wnt signaling to produce surplus ACs in zebrafish embryos. Together, our data support that downstream members of a conserved regulatory network involving Wnt signaling and Foxf1 function on a nr2f1a enhancer to promote AC differentiation in the zebrafish heart.

Transcriptional regulation of the alternative sex hormone-binding globulin promoter by KLF4

In most mammals the major site of sex hormone-binding globulin (SHBG) synthesis is the liver wherefrom it is secreted into the bloodstream and is the primary determinant of sex steroid access to target tissues. The minor site of SHBG synthesis is the testis and in lower mammals testicular SHBG has long been known to be synthesized and secreted by Sertoli cells. However, human testicular SHBG is expressed in developing germ cells from an upstream alternative promoter (altP-SHBG). Transcripts arising from this region comprise an alternative first exon (1A) with the resultant protein confined to the acrosomal compartment of the mature spermatozoa. I have dissected the regulatory components of the alternative SHBG promoter and identified motifs that are required for optimal transcriptional activity from this region. Transcriptional activity is driven by two CACCC elements that appear to be functionally redundant. The transcription factor KLF4 interacts with promoter the region spanning these elements in vivo. Knockdown of Klf4 results in decreased altP-SHBG activity, while Klf4 overexpression relieves the effects of knockdown. Based on their shared patterns of expression in vivo, I conclude that KLF4 is a transcriptional regulator of SHBG in male germ cells.

Temporally resolved early BMP-driven transcriptional cascade during human amnion specification

Amniogenesis, a process critical for continuation of healthy pregnancy, is triggered in a collection of pluripotent epiblast cells as the human embryo implants. Previous studies have established that BMP signaling is a major driver of this lineage specifying process, but the downstream BMP-dependent transcriptional networks that lead to successful amniogenesis remain to be identified. This is, in part, due to the current lack of a robust and reproducible model system that enables mechanistic investigations exclusively into amniogenesis. Here, we developed an improved model of early amnion specification, using a human pluripotent stem cell-based platform in which the activation of BMP signaling is controlled and synchronous. Uniform amniogenesis is seen within 48 hours after BMP activation, and the resulting cells share transcriptomic characteristics with amnion cells of a gastrulating human embryo. Using detailed time-course transcriptomic analyses, we established a previously uncharacterized BMP-dependent amniotic transcriptional cascade, and identified markers that represent five distinct stages of amnion fate specification; the expression of selected markers was validated in early post-implantation macaque embryos. Moreover, a cohort of factors that could potentially control specific stages of amniogenesis was identified, including the transcription factor TFAP2A. Functionally, we determined that, once amniogenesis is triggered by the BMP pathway, TFAP2A controls the progression of amniogenesis. This work presents a temporally resolved transcriptomic resource for several previously uncharacterized amniogenesis states and demonstrates a critical intermediate role for TFAP2A during amnion fate specification.

Distinct regulatory networks control toxin gene expression in elapid and viperid snakes

BACKGROUND

Venom systems are ideal models to study genetic regulatory mechanisms that underpin evolutionary novelty. Snake venom glands are thought to share a common origin, but there are major distinctions between venom toxins from the medically significant snake families Elapidae and Viperidae, and toxin gene regulatory investigations in elapid snakes have been limited. Here, we used high-throughput RNA-sequencing to profile gene expression and microRNAs between active (milked) and resting (unmilked) venom glands in an elapid (Eastern Brown Snake, Pseudonaja textilis), in addition to comparative genomics, to identify cis- and trans-acting regulation of venom production in an elapid in comparison to viperids (Crotalus viridis and C. tigris).

RESULTS

Although there is conservation in high-level mechanistic pathways regulating venom production (unfolded protein response, Notch signaling and cholesterol homeostasis), there are differences in the regulation of histone methylation enzymes, transcription factors, and microRNAs in venom glands from these two snake families. Histone methyltransferases and transcription factor (TF) specificity protein 1 (Sp1) were highly upregulated in the milked elapid venom gland in comparison to the viperids, whereas nuclear factor I (NFI) TFs were upregulated after viperid venom milking. Sp1 and NFI cis-regulatory elements were common to toxin gene promoter regions, but many unique elements were also present between elapid and viperid toxins. The presence of Sp1 binding sites across multiple elapid toxin gene promoter regions that have been experimentally determined to regulate expression, in addition to upregulation of Sp1 after venom milking, suggests this transcription factor is involved in elapid toxin expression. microRNA profiles were distinctive between milked and unmilked venom glands for both snake families, and microRNAs were predicted to target a diversity of toxin transcripts in the elapid P. textilis venom gland, but only snake venom metalloproteinase transcripts in the viperid C. viridis venom gland. These results suggest differences in toxin gene posttranscriptional regulation between the elapid P. textilis and viperid C. viridis.

CONCLUSIONS

Our comparative transcriptomic and genomic analyses between toxin genes and isoforms in elapid and viperid snakes suggests independent toxin regulation between these two snake families, demonstrating multiple different regulatory mechanisms underpin a venomous phenotype.

Cell-type-directed design of synthetic enhancers

Transcriptional enhancers act as docking stations for combinations of transcription factors and thereby regulate spatiotemporal activation of their target genes1. It has been a long-standing goal in the field to decode the regulatory logic of an enhancer and to understand the details of how spatiotemporal gene expression is encoded in an enhancer sequence. Here we show that deep learning models2-6, can be used to efficiently design synthetic, cell-type-specific enhancers, starting from random sequences, and that this optimization process allows detailed tracing of enhancer features at single-nucleotide resolution. We evaluate the function of fully synthetic enhancers to specifically target Kenyon cells or glial cells in the fruit fly brain using transgenic animals. We further exploit enhancer design to create 'dual-code' enhancers that target two cell types and minimal enhancers smaller than 50 base pairs that are fully functional. By examining the state space searches towards local optima, we characterize enhancer codes through the strength, combination and arrangement of transcription factor activator and transcription factor repressor motifs. Finally, we apply the same strategies to successfully design human enhancers, which adhere to enhancer rules similar to those of Drosophila enhancers. Enhancer design guided by deep learning leads to better understanding of how enhancers work and shows that their code can be exploited to manipulate cell states.

FXR controls insulin content by regulating Foxa2-mediated insulin transcription

Farnesoid X receptor (FXR) is a nuclear ligand-activated receptor of bile acids that plays a role in the modulation of insulin content. However, the underlying molecular mechanisms remain unclear. Forkhead box a2 (Foxa2) is an important nuclear transcription factor in pancreatic β-cells and is involved in β-cell function. We aimed to explore the signaling mechanism downstream of FXR to regulate insulin content and underscore its association with Foxa2 and insulin gene (Ins) transcription. All experiments were conducted on FXR transgenic mice, INS-1 823/13 cells, and diabetic Goto-Kakizaki (GK) rats undergoing sham or Roux-en-Y gastric bypass (RYGB) surgery. Islets from FXR knockout mice and INS-1823/13 cells with FXR knockdown exhibited substantially lower insulin levels than that of controls. This was accompanied by decreased Foxa2 expression and Ins transcription. Conversely, FXR overexpression increased insulin content, concomitant with enhanced Foxa2 expression and Ins transcription in INS-1 823/13 cells. Moreover, FXR knockdown reduced FXR recruitment and H3K27 trimethylation in the Foxa2 promoter. Importantly, Foxa2 overexpression abrogated the adverse effects of FXR knockdown on Ins transcription and insulin content in INS-1 823/13 cells. Notably, RYGB surgery led to improved insulin content in diabetic GK rats, which was accompanied by upregulated FXR and Foxa2 expression and Ins transcription. Collectively, these data suggest that Foxa2 serves as the target gene of FXR in β-cells and mediates FXR-enhanced Ins transcription. Additionally, the upregulated FXR/Foxa2 signaling cascade could contribute to the enhanced insulin content in diabetic GK rats after RYGB.

HNF4A guides the MLL4 complex to establish and maintain H3K4me1 at gene regulatory elements

Hepatocyte nuclear factor 4A (HNF4A/NR2a1), a transcriptional regulator of hepatocyte identity, controls genes that are crucial for liver functions, primarily through binding to enhancers. In mammalian cells, active and primed enhancers are marked by monomethylation of histone 3 (H3) at lysine 4 (K4) (H3K4me1) in a cell type-specific manner. How this modification is established and maintained at enhancers in connection with transcription factors (TFs) remains unknown. Using analysis of genome-wide histone modifications, TF binding, chromatin accessibility and gene expression, we show that HNF4A is essential for an active chromatin state. Using HNF4A loss and gain of function experiments in vivo and in cell lines in vitro, we show that HNF4A affects H3K4me1, H3K27ac and chromatin accessibility, highlighting its contribution to the establishment and maintenance of a transcriptionally permissive epigenetic state. Mechanistically, HNF4A interacts with the mixed-lineage leukaemia 4 (MLL4) complex facilitating recruitment to HNF4A-bound regions. Our findings indicate that HNF4A enriches H3K4me1, H3K27ac and establishes chromatin opening at transcriptional regulatory regions.

BORIS/CTCFL epigenetically reprograms clustered CTCF binding sites into alternative transcriptional start sites

BACKGROUND

Pervasive usage of alternative promoters leads to the deregulation of gene expression in carcinogenesis and may drive the emergence of new genes in spermatogenesis. However, little is known regarding the mechanisms underpinning the activation of alternative promoters.

RESULTS

Here we describe how alternative cancer-testis-specific transcription is activated. We show that intergenic and intronic CTCF binding sites, which are transcriptionally inert in normal somatic cells, could be epigenetically reprogrammed into active de novo promoters in germ and cancer cells. BORIS/CTCFL, the testis-specific paralog of the ubiquitously expressed CTCF, triggers the epigenetic reprogramming of CTCF sites into units of active transcription. BORIS binding initiates the recruitment of the chromatin remodeling factor, SRCAP, followed by the replacement of H2A histone with H2A.Z, resulting in a more relaxed chromatin state in the nucleosomes flanking the CTCF binding sites. The relaxation of chromatin around CTCF binding sites facilitates the recruitment of multiple additional transcription factors, thereby activating transcription from a given binding site. We demonstrate that the epigenetically reprogrammed CTCF binding sites can drive the expression of cancer-testis genes, long noncoding RNAs, retro-pseudogenes, and dormant transposable elements.

CONCLUSIONS

Thus, BORIS functions as a transcription factor that epigenetically reprograms clustered CTCF binding sites into transcriptional start sites, promoting transcription from alternative promoters in both germ cells and cancer cells.

Why Are Nucleosome Breathing Dynamics Asymmetric?

In eukaryotic cells, DNA is bound to nucleosomes, but DNA segments occasionally unbind in the process known as nucleosome breathing. Although DNA can unwrap simultaneously from both ends of the nucleosome (symmetric breathing), experiments indicate that DNA prefers to dissociate from only one end (asymmetric breathing). However, the molecular origin of the asymmetry is not understood. We developed a new theoretical approach that gives microscopic explanations of asymmetric breathing. It is based on a stochastic description that leads to a comprehensive evaluation of dynamics by using effective free-energy landscapes. It is shown that asymmetric breathing follows the kinetically preferred pathways. In addition, it is also found that asymmetric breathing leads to a faster target search by transcription factors. Theoretical predictions, supported by computer simulations, agree with experiments. It is proposed that nature utilizes the symmetry of nucleosome breathing to achieve a better dynamic accessibility of chromatin for more efficient genetic regulation.

The Neural Progenitor Cell-Associated Transcription Factor FoxG1 Regulates Cardiac Epicardial Cell Proliferation

The epicardium is a layer of mesothelial cells that covers the surface of the heart. During development, epicardial cells undergo epithelial-to-mesenchymal transition (EMT) to form multipotent precursors that migrate into the heart and contribute to the coronary vasculature by differentiating into adventitial fibroblasts, smooth muscle cells, and endothelial cells. Epicardial cells also provide paracrine signals to cardiac myocytes that are required for appropriate heart growth. In adult hearts, a similar process of epicardial cell EMT, migration, and differentiation occurs after myocardial infarction (MI, heart attack). Pathological cardiac hypertrophy is associated with fibrosis, negative remodeling, and reduced cardiac function. In contrast, aerobic exercises such as swimming and running promote physiological (i.e., beneficial) hypertrophy, which is associated with angiogenesis and improved cardiac function. As epicardial cell function(s) during physiological hypertrophy are poorly understood, we analyzed and compared the native epicardial cells isolated directly from the hearts of running-exercised mice and age-matched, nonrunning littermates. To obtain epicardial cells, we enzymatically digested the surfaces of whole hearts and performed magnetic-activated cell sorting (MACS) with antibodies against CD104 (integrin β4). By cDNA microarray assays, we identified genes with increased transcription in epicardial cells after running exercise; these included FoxG1, a transcription factor that controls neural progenitor cell proliferation during brain development and Snord116, a small noncoding RNA that coordinates expression of genes with epigenetic, circadian, and metabolic functions. In cultured epicardial cells, shRNA-mediated FoxG1 knockdown significantly decreased cell proliferation, as well as Snord116 expression. Our results demonstrate that FoxG1 regulates epicardial proliferation, and suggest it may affect cardiac remodeling.

It Takes Two to Tango: The Interplay between Prostate Cancer and Its Microenvironment from an Epigenetic Perspective

Prostate cancer is the second most common cancer in men worldwide and is associated with high morbidity and mortality. Consequently, there is an urgent unmet need for novel treatment avenues. In addition to somatic genetic alterations, deviations in the epigenetic landscape of cancer cells and their tumor microenvironment (TME) are critical drivers of prostate cancer initiation and progression. Unlike genomic mutations, epigenetic modifications are potentially reversible. Therefore, the inhibition of aberrant epigenetic modifications represents an attractive and exciting novel treatment strategy for castration-resistant prostate cancer patients. Moreover, drugs targeting the epigenome also exhibit synergistic interactions with conventional therapeutics by directly enhancing their anti-tumorigenic properties by "priming" the tumor and tumor microenvironment to increase drug sensitivity. This review summarizes the major epigenetic alterations in prostate cancer and its TME, and their involvement in prostate tumorigenesis, and discusses the impact of epigenome-targeted therapies.

MethMotif.Org 2024: a database integrating context-specific transcription factor-binding motifs with DNA methylation patterns

MethMotif (https://methmotif.org) is a publicly available database that provides a comprehensive repository of transcription factor (TF)-binding profiles, enriched with DNA methylation patterns. In this release, we have enhanced the platform, expanding our initial collection to over 700 position weight matrices (PWM), all of which include DNA methylation profiles. One of the key advancements in this release is the segregation of TF-binding motifs based on their cofactors and DNA methylation status. We have previously demonstrated that gene ontology (GO) enriched terms associated with TF target genes may differ based on their association with alternative cofactors and DNA methylation status. MethMotif provides precomputed GO annotations for each human TF of interest, as well as for TF-co-TF complexes, enabling a comprehensive analysis of TF functions in the context of their co-factors. Additionally, MethMotif has been updated to encompass data for two new species, Mus musculus and Arabidopsis thaliana, widening its applicability to a broader community. MethMotif stands out as the first and only TF-binding motifs database to incorporate context-specific PWM coupled with epigenetic information, thereby enlightening context-specific TF functions. This enhancement allows the community to explore and gain deeper insights into the regulatory mechanisms governing transcriptional processes.

Systems Biology Approach to Analyze Microarray Datasets for Identification of Disease-Causing Genes: Case Study of Oral Squamous Cell Carcinoma

The discovery of potential disease-causing genes can aid medical progress. The post-genomic era has made this a more difficult task. Modern high-throughput methods have not solved the problem of identifying disease genes. Conventional methods cannot be used to investigate many rare or lethal diseases. Monitoring gene expression values in different samples using microarray technology is one of the best and most accurate ways to identify disease-causing genes. One of the most recent advances in experimental molecular biology is microarrays, which allow researchers to simultaneously monitor the expression levels of thousands of genes. Statistical analysis of microarray data might aid gene discovery by revealing pathways related to the target gene and facilitating identification of candidate genes. Systems biology, an interdisciplinary approach, has emerged as a crucial analytic tool with the potential to reveal previously unidentified causes and consequences of human illness. Genetic, environmental, immunological, or neurological factors have been implicated in the developing complex disorders like cancer. Because of this, it is important to approach the study of such disease from a novel perspective. The system biology approach allows us to rapidly identify disease-causing genes and assess their viability as therapeutic targets. This chapter demonstrates systems biology approaches to identify candidate genes using public database. Oral squamous cell carcinoma (OSCC) is used as a model disease to show how systems biology can be used successfully to identify and prioritize disease genes.

Is euchromatin really open in the cell?

Genomic DNA is wrapped around a core histone octamer and forms a nucleosome. In higher eukaryotic cells, strings of nucleosomes are irregularly folded as chromatin domains that act as functional genome units. According to a typical textbook model, chromatin can be categorized into two types, euchromatin and heterochromatin, based on its degree of compaction. Euchromatin is open, while heterochromatin is closed and condensed. However, is euchromatin really open in the cell? New evidence from genomics and advanced imaging studies has revealed that euchromatin consists of condensed liquid-like domains. Condensed chromatin seems to be the default chromatin state in higher eukaryotic cells. We discuss this novel view of euchromatin in the cell and how the revealed organization is relevant to genome functions.

FOXO transcription factors as mediators of stress adaptation

The forkhead box protein O (FOXO, consisting of FOXO1, FOXO3, FOXO4 and FOXO6) transcription factors are the mammalian orthologues of Caenorhabditis elegans DAF-16, which gained notoriety for its capability to double lifespan in the absence of daf-2 (the gene encoding the worm insulin receptor homologue). Since then, research has provided many mechanistic details on FOXO regulation and FOXO activity. Furthermore, conditional knockout experiments have provided a wealth of data as to how FOXOs control development and homeostasis at the organ and organism levels. The lifespan-extending capabilities of DAF-16/FOXO are highly correlated with their ability to induce stress response pathways. Exogenous and endogenous stress, such as cellular redox stress, are considered the main drivers of the functional decline that characterizes ageing. Functional decline often manifests as disease, and decrease in FOXO activity indeed negatively impacts on major age-related diseases such as cancer and diabetes. In this context, the main function of FOXOs is considered to preserve cellular and organismal homeostasis, through regulation of stress response pathways. Paradoxically, the same FOXO-mediated responses can also aid the survival of dysfunctional cells once these eventually emerge. This general property to control stress responses may underlie the complex and less-evident roles of FOXOs in human lifespan as opposed to model organisms such as C. elegans.

Functional Genome Analyses Reveal the Molecular Basis of Oil Accumulation in Developing Seeds of Castor Beans

Castor (Ricinus communis L.) seeds produce abundant ricinoleic acid during seed maturation, which is important for plant development and human demands. Ricinoleic acid, as a unique hydroxy fatty acid (HFA), possesses a distinct bond structure that could be used as a substitute for fossil fuels. Here, we identified all homologous genes related to glycolysis, hydroxy fatty acid biosynthesis, and triacylglycerol (TAG) accumulation in castor seeds. Furthermore, we investigated their expression patterns globally during five seed development stages. We characterized a total of 66 genes involved in the glycolysis pathway, with the majority exhibiting higher expression levels during the early stage of castor bean seed development. This metabolic process provided abundant acetyl-CoA for fatty acid (FA) biosynthesis. Subsequently, we identified 82 genes involved in the processes of de novo FA biosynthesis and TAG assembly, with the majority exhibiting high expression levels during the middle or late stages. In addition, we examined the expression patterns of the transcription factors involved in carbohydrate and oil metabolism. For instance, RcMYB73 and RcERF72 exhibited high expression levels during the early stage, whereas RcWRI1, RcABI3, and RcbZIP67 showed relatively higher expression levels during the middle and late stages, indicating their crucial roles in seed development and oil accumulation. Our study suggests that the high HFA production in castor seeds is attributed to the interaction of multiple genes from sugar transportation to lipid droplet packaging. Therefore, this research comprehensively characterizes all the genes related to glycolysis, fatty acid biosynthesis, and triacylglycerol (TAG) accumulation in the castor and provides novel insight into exploring the genetic mechanisms underlying seed oil accumulation in the endosperm of castor beans.

Inhibition of the lysine demethylase LSD1 modulates the balance between inflammatory and antiviral responses against coronaviruses

Innate immune responses to coronavirus infections are highly cell specific. Tissue-resident macrophages, which are infected by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in patients but are inconsistently infected in vitro, exert critical but conflicting effects by secreting both antiviral type I interferons (IFNs) and tissue-damaging inflammatory cytokines. Steroids, the only class of host-targeting drugs approved for the treatment of coronavirus disease 2019 (COVID-19), indiscriminately suppress both responses, possibly impairing viral clearance. Here, we established in vitro cell culture systems that enabled us to separately investigate the cell-intrinsic and cell-extrinsic proinflammatory and antiviral activities of mouse macrophages infected with the prototypical murine coronavirus MHV-A59. We showed that the nuclear factor κB-dependent inflammatory response to viral infection was selectively inhibited by loss of the lysine demethylase LSD1, which was previously implicated in innate immune responses to cancer, with negligible effects on the antiviral IFN response. LSD1 ablation also enhanced an IFN-independent antiviral response, blocking viral egress through the lysosomal pathway. The macrophage-intrinsic antiviral and anti-inflammatory activity of Lsd1 inhibition was confirmed in vitro and in a humanized mouse model of SARS-CoV-2 infection. These results suggest that LSD1 controls innate immune responses against coronaviruses at multiple levels and provide a mechanistic rationale for potentially repurposing LSD1 inhibitors for COVID-19 treatment.

Characterizing Different Modes of Interplay Between Rap1 and H3 Using Inducible H3-depletion Yeast

Histones and transcription factors (TFs) are two important DNA-binding proteins that interact, compete, and together regulate transcriptional processes in response to diverse internal and external stimuli. Condition-specific depletion of histones in Saccharomyces cerevisiae using a galactose-inducible H3 promoter provides a suitable framework for examining transcriptional alteration resulting from reduced nucleosome content. However, the effect on DNA binding activities of TFs is yet to be fully explored. In this work, we combine ChIP-seq of H3 with RNA-seq to elucidate the genome-scale relationships between H3 occupancy patterns and transcriptional dynamics before and after global H3 depletion. ChIP-seq of Rap1 is also conducted in the H3-depletion and control treatments, to investigate the interplay between this master regulator TF and nucleosomal H3, and to explore the impact on diverse transcriptional responses of different groups of target genes and functions. Ultimately, we propose a working model and testable hypotheses regarding the impact of global and local H3 depletion on transcriptional changes. We also demonstrate different potential modes of interaction between Rap1 and H3, which sheds light on the potential multifunctional regulatory capabilities of Rap1 and potentially other pioneer factors.

Multifaceted SOX2-chromatin interaction underpins pluripotency progression in early embryos

Pioneer transcription factors (TFs), such as OCT4 and SOX2, play crucial roles in pluripotency regulation. However, the master TF-governed pluripotency regulatory circuitry was largely inferred from cultured cells. In this work, we investigated SOX2 binding from embryonic day 3.5 (E3.5) to E7.5 in the mouse. In E3.5 inner cell mass (ICM), SOX2 regulates the ICM-trophectoderm program but is dispensable for opening global enhancers. Instead, SOX2 occupies preaccessible enhancers in part opened by early-stage expressing TFs TFAP2C and NR5A2. SOX2 then widely redistributes when cells adopt naive and formative pluripotency by opening enhancers or poising them for rapid future activation. Hence, multifaceted pioneer TF-enhancer interaction underpins pluripotency progression in embryos, including a distinctive state in E3.5 ICM that bridges totipotency and pluripotency.

Old school, new rules: floral meristem development revealed by 3D gene expression atlases and high-resolution transcription factor-chromatin dynamics

The intricate morphology of the flower is primarily established within floral meristems in which floral organs will be defined and from where the developing flower will emerge. Floral meristem development involves multiscale-level regulation, including lineage and positional mechanisms for establishing cell-type identity, and transcriptional regulation mediated by changes in the chromatin environment. However, many key aspects of floral meristem development remain to be determined, such as: 1) the exact role of cellular location in connecting transcriptional inputs to morphological outcomes, and 2) the precise interactions between transcription factors and chromatin regulators underlying the transcriptional networks that regulate the transition from cell proliferation to differentiation during floral meristem development. Here, we highlight recent studies addressing these points through newly developed spatial reconstruction techniques and high-resolution transcription factor-chromatin environment interactions in the model plant Arabidopsis thaliana. Specifically, we feature studies that reconstructed 3D gene expression atlases of the floral meristem. We also discuss how the precise timing of floral meristem specification, floral organ patterning, and floral meristem termination is determined through temporally defined epigenetic dynamics for fine-tuning of gene expression. These studies offer fresh insights into the well-established principles of floral meristem development and outline the potential for further advances in this field in an age of integrated, powerful, multiscale resolution approaches.

Transcription factor condensates and signaling driven transcription

Transcription Factor (TF) condensates are a heterogenous mix of RNA, DNA, and multiple co-factor proteins capable of modulating the transcriptional response of the cell. The dynamic nature and the spatial location of TF-condensates in the 3D nuclear space is believed to provide a fast response, which is on the same pace as the signaling cascade and yet ever-so-specific in the crowded environment of the nucleus. However, the current understanding of how TF-condensates can achieve these feet so quickly and efficiently is still unclear. In this review, we draw parallels with other protein condensates and share our speculations on how the nucleus uses these TF-condensates to achieve high transcriptional specificity and fidelity. We discuss the various constituents of TF-condensates, their properties, and the known and unknown functions of TF-condensates with a particular focus on steroid signaling-induced transcriptional programs.

The role of pioneer transcription factors in the induction of direct cellular reprogramming

Regenerative medicine is a highly advanced medical field that aims to restore tissues and organs lost due to diseases and injury using a person's own cells or those of others. Direct cellular reprogramming is a promising technology that can directly induce cell-fate conversion from terminally differentiated cells to other cell types and is expected to play a pivotal role in applications in regenerative medicine. The induction of direct cellular reprogramming requires one or more master transcription factors with the potential to reconstitute cell type-specific transcription factor networks. The set of master transcription factors may contain unique transcription factors called pioneer factors that can open compacted chromatin structures and drive the transcriptional activation of target genes. Therefore, pioneer factors may play a central role in direct cellular reprogramming. However, our understanding of the molecular mechanisms by which pioneer factors induce cell-fate conversion is still limited. This review briefly summarizes the outcomes of recent findings and discusses future perspectives, focusing on the role of pioneer factors in direct cellular reprogramming.

ΔNp63α facilitates proliferation and migration, and modulates the chromatin landscape in intrahepatic cholangiocarcinoma cells

p63 plays a crucial role in epithelia-originating tumours; however, its role in intrahepatic cholangiocarcinoma (iCCA) has not been completely explored. Our study revealed the oncogenic properties of p63 in iCCA and identified the major expressed isoform as ΔNp63α. We collected iCCA clinical data from The Cancer Genome Atlas database and analyzed p63 expression in iCCA tissue samples. We further established genetically modified iCCA cell lines in which p63 was overexpressed or knocked down to study the protein function/function of p63 in iCCA. We found that cells overexpressing p63, but not p63 knockdown counterparts, displayed increased proliferation, migration, and invasion. Transcriptome analysis showed that p63 altered the iCCA transcriptome, particularly by affecting cell adhesion-related genes. Moreover, chromatin accessibility decreased at p63 target sites when p63 binding was lost and increased when p63 binding was gained. The majority of the p63 bound sites were located in the distal intergenic regions and showed strong enhancer marks; however, active histone modifications around the Transcription Start Site changed as p63 expression changed. We also detected an interaction between p63 and the chromatin structural protein YY1. Taken together, our results suggest an oncogenic role for p63 in iCCA.

The C. elegans Myc-family of transcription factors coordinate a dynamic adaptive response to dietary restriction

Dietary restriction (DR), the process of decreasing overall food consumption over an extended period of time, has been shown to increase longevity across evolutionarily diverse species and delay the onset of age-associated diseases in humans. In Caenorhabditis elegans, the Myc-family transcription factors (TFs) MXL-2 (Mlx) and MML-1 (MondoA/ChREBP), which function as obligate heterodimers, and PHA-4 (orthologous to forkhead box transcription factor A) are both necessary for the full physiological benefits of DR. However, the adaptive transcriptional response to DR and the role of MML-1::MXL-2 and PHA-4 remains elusive. We identified the transcriptional signature of C. elegans DR, using the eat-2 genetic model, and demonstrate broad changes in metabolic gene expression in eat-2 DR animals, which requires both mxl-2 and pha-4. While the requirement for these factors in DR gene expression overlaps, we found many of the DR genes exhibit an opposing change in relative gene expression in eat-2;mxl-2 animals compared to wild-type, which was not observed in eat-2 animals with pha-4 loss. We further show functional deficiencies of the mxl-2 loss in DR outside of lifespan, as eat-2;mxl-2 animals exhibit substantially smaller brood sizes and lay a proportion of dead eggs, indicating that MML-1::MXL-2 has a role in maintaining the balance between resource allocation to the soma and to reproduction under conditions of chronic food scarcity. While eat-2 animals do not show a significantly different metabolic rate compared to wild-type, we also find that loss of mxl-2 in DR does not affect the rate of oxygen consumption in young animals. The gene expression signature of eat-2 mutant animals is consistent with optimization of energy utilization and resource allocation, rather than induction of canonical gene expression changes associated with acute metabolic stress -such as induction of autophagy after TORC1 inhibition. Consistently, eat-2 animals are not substantially resistant to stress, providing further support to the idea that chronic DR may benefit healthspan and lifespan through efficient use of limited resources rather than broad upregulation of stress responses, and also indicates that MML-1::MXL-2 and PHA-4 may have different roles in promotion of benefits in response to different pro-longevity stimuli.

Epigenome editing in cancer: Advances and challenges for potential therapeutic options

Cancers are diseases caused by genetic and non-genetic environmental factors. Epigenetic alterations, some attributed to non-genetic factors, can lead to cancer development. Epigenetic changes can occur in tumor suppressors or oncogenes, or they may contribute to global cell state changes, making cells abnormal. Recent advances in gene editing technology show potential for cancer treatment. Herein, we will discuss our current knowledge of epigenetic alterations occurring in cancer and epigenetic editing technologies that can be applied to developing therapeutic options.

Human pluripotent stem cell fate trajectories toward lung and hepatocyte progenitors

In this study, we interrogate molecular mechanisms underlying the specification of lung progenitors from human pluripotent stem cells (hPSCs). We employ single-cell RNA-sequencing with high temporal precision, alongside an optimized differentiation protocol, to elucidate the transcriptional hierarchy of lung specification to chart the associated single-cell trajectories. Our findings indicate that Sonic hedgehog, TGF-β, and Notch activation are essential within an ISL1/NKX2-1 trajectory, leading to the emergence of lung progenitors during the foregut endoderm phase. Additionally, the induction of HHEX delineates an alternate trajectory at the early definitive endoderm stage, preceding the lung pathway and giving rise to a significant hepatoblast population. Intriguingly, neither KDR+ nor mesendoderm progenitors manifest as intermediate stages in the lung and hepatic lineage development. Our multistep model offers insights into lung organogenesis and provides a foundation for in-depth study of early human lung development and modeling using hPSCs.

Systematic Dissection of Sequence Features Affecting the Binding Specificity of a Pioneer Factor Reveals Binding Synergy Between FOXA1 and AP-1

Despite the unique ability of pioneer transcription factors (PFs) to target nucleosomal sites in closed chromatin, they only bind a small fraction of their genomic motifs. The underlying mechanism of this selectivity is not well understood. Here, we design a high-throughput assay called ChIP-ISO to systematically dissect sequence features affecting the binding specificity of a classic PF, FOXA1. Combining ChIP-ISO with in vitro and neural network analyses, we find that 1) FOXA1 binding is strongly affected by co-binding TFs AP-1 and CEBPB, 2) FOXA1 and AP-1 show binding cooperativity in vitro, 3) FOXA1's binding is determined more by local sequences than chromatin context, including eu-/heterochromatin, and 4) AP-1 is partially responsible for differential binding of FOXA1 in different cell types. Our study presents a framework for elucidating genetic rules underlying PF binding specificity and reveals a mechanism for context-specific regulation of its binding.

FOXO1 is a master regulator of CAR T memory programming

Poor CAR T persistence limits CAR T cell therapies for B cell malignancies and solid tumors1,2. The expression of memory-associated genes such as TCF7 (protein name TCF1) is linked to response and long-term persistence in patients3-7, thereby implicating memory programs in therapeutic efficacy. Here, we demonstrate that the pioneer transcription factor, FOXO1, is responsible for promoting memory programs and restraining exhaustion in human CAR T cells. Pharmacologic inhibition or gene editing of endogenous FOXO1 in human CAR T cells diminished the expression of memory-associated genes, promoted an exhaustion-like phenotype, and impaired antitumor activity in vitro and in vivo. FOXO1 overexpression induced a gene expression program consistent with T cell memory and increased chromatin accessibility at FOXO1 binding motifs. FOXO1-overexpressing cells retained function, memory potential, and metabolic fitness during settings of chronic stimulation and exhibited enhanced persistence and antitumor activity in vivo. In contrast, TCF1 overexpression failed to enforce canonical memory programs or enhance CAR T cell potency. Importantly, endogenous FOXO1 activity correlated with CAR T and TIL responses in patients, underscoring its clinical relevance in cancer immunotherapy. Our results demonstrate that memory reprogramming through FOXO1 can enhance the persistence and potency of human CAR T cells and highlights the utility of pioneer factors, which bind condensed chromatin and induce local epigenetic remodeling, for optimizing therapeutic T cell states.

International Union of Basic and Clinical Pharmacology CXIII: Nuclear Receptor Superfamily-Update 2023

The NR superfamily comprises 48 transcription factors in humans that control a plethora of gene network programs involved in a wide range of physiologic processes. This review will summarize and discuss recent progress in NR biology and drug development derived from integrating various approaches, including biophysical techniques, structural studies, and translational investigation. We also highlight how defective NR signaling results in various diseases and disorders and how NRs can be targeted for therapeutic intervention via modulation via binding to synthetic lipophilic ligands. Furthermore, we also review recent studies that improved our understanding of NR structure and signaling. SIGNIFICANCE STATEMENT: Nuclear receptors (NRs) are ligand-regulated transcription factors that are critical regulators of myriad physiological processes. NRs serve as receptors for an array of drugs, and in this review, we provide an update on recent research into the roles of these drug targets.

Gene regulatory network inference in the era of single-cell multi-omics

The interplay between chromatin, transcription factors and genes generates complex regulatory circuits that can be represented as gene regulatory networks (GRNs). The study of GRNs is useful to understand how cellular identity is established, maintained and disrupted in disease. GRNs can be inferred from experimental data - historically, bulk omics data - and/or from the literature. The advent of single-cell multi-omics technologies has led to the development of novel computational methods that leverage genomic, transcriptomic and chromatin accessibility information to infer GRNs at an unprecedented resolution. Here, we review the key principles of inferring GRNs that encompass transcription factor-gene interactions from transcriptomics and chromatin accessibility data. We focus on the comparison and classification of methods that use single-cell multimodal data. We highlight challenges in GRN inference, in particular with respect to benchmarking, and potential further developments using additional data modalities.

Joint sequence & chromatin neural networks characterize the differential abilities of Forkhead transcription factors to engage inaccessible chromatin

The DNA-binding activities of transcription factors (TFs) are influenced by both intrinsic sequence preferences and extrinsic interactions with cell-specific chromatin landscapes and other regulatory proteins. Disentangling the roles of these binding determinants remains challenging. For example, the FoxA subfamily of Forkhead domain (Fox) TFs are known pioneer factors that can bind to relatively inaccessible sites during development. Yet FoxA TF binding also varies across cell types, pointing to a combination of intrinsic and extrinsic forces guiding their binding. While other Forkhead domain TFs are often assumed to have pioneering abilities, how sequence and chromatin features influence the binding of related Fox TFs has not been systematically characterized. Here, we present a principled approach to compare the relative contributions of intrinsic DNA sequence preference and cell-specific chromatin environments to a TF's DNA-binding activities. We apply our approach to investigate how a selection of Fox TFs (FoxA1, FoxC1, FoxG1, FoxL2, and FoxP3) vary in their binding specificity. We over-express the selected Fox TFs in mouse embryonic stem cells, which offer a platform to contrast each TF's binding activity within the same preexisting chromatin background. By applying a convolutional neural network to interpret the Fox TF binding patterns, we evaluate how sequence and preexisting chromatin features jointly contribute to induced TF binding. We demonstrate that Fox TFs bind different DNA targets, and drive differential gene expression patterns, even when induced in identical chromatin settings. Despite the association between Forkhead domains and pioneering activities, the selected Fox TFs display a wide range of affinities for preexiting chromatin states. Using sequence and chromatin feature attribution techniques to interpret the neural network predictions, we show that differential sequence preferences combined with differential abilities to engage relatively inaccessible chromatin together explain Fox TF binding patterns at individual sites and genome-wide.

Pioneer activity distinguishes activating from non-activating SOX2 binding sites

Genome-wide transcriptional activity involves the binding of many transcription factors (TFs) to thousands of sites in the genome. Pioneer TFs are a class of TFs that maintain open chromatin and allow non-pioneer TFs access to their target sites. Determining which TF binding sites directly drive transcription remains a challenge. Here, we use acute protein depletion of the pioneer TF SOX2 to establish its functionality in maintaining chromatin accessibility. We show that thousands of accessible sites are lost within an hour of protein depletion, indicating rapid turnover of these sites in the absence of the pioneer factor. To understand the relationship with transcription, we performed nascent transcription analysis and found that open chromatin sites that are maintained by SOX2 are highly predictive of gene expression, in contrast to all other SOX2 binding sites. We use CRISPR-Cas9 genome editing in the Klf2 locus to functionally validate a predicted regulatory element. We conclude that the regulatory activity of SOX2 is exerted mainly at sites where it maintains accessibility and that other binding sites are largely dispensable for gene regulation.

In Silico Verification of Predicted Potential Promoter Sequences in the Rice (Oryza sativa) Genome

The exact identification of promoter sequences remains a serious problem in computational biology, as the promoter prediction algorithms under development continue to produce false-positive results. Therefore, to fully assess the validity of predicted sequences, it is necessary to perform a comprehensive test of their properties, such as the presence of downstream transcribed DNA regions behind them, or chromatin accessibility for transcription factor binding. In this paper, we examined the promoter sequences of chromosome 1 of the rice Oryza sativa genome from the Database of Potential Promoter Sequences predicted using a mathematical algorithm based on the derivation and calculation of statistically significant promoter classes. In this paper TATA motifs and cis-regulatory elements were identified in the predicted promoter sequences. We also verified the presence of potential transcription start sites near the predicted promoters by analyzing CAGE-seq data. We searched for unannotated transcripts behind the predicted sequences by de novo assembling transcripts from RNA-seq data. We also examined chromatin accessibility in the region of the predicted promoters by analyzing ATAC-seq data. As a result of this work, we identified the predicted sequences that are most likely to be promoters for further experimental validation in an in vivo or in vitro system.

Understanding fundamental principles of enhancer biology at a model locus: Analysing the structure and function of an enhancer cluster at the α-globin locus

Despite ever-increasing accumulation of genomic data, the fundamental question of how individual genes are switched on during development, lineage-specification and differentiation is not fully answered. It is widely accepted that this involves the interaction between at least three fundamental regulatory elements: enhancers, promoters and insulators. Enhancers contain transcription factor binding sites which are bound by transcription factors (TFs) and co-factors expressed during cell fate decisions and maintain imposed patterns of activation, at least in part, via their epigenetic modification. This information is transferred from enhancers to their cognate promoters often by coming into close physical proximity to form a 'transcriptional hub' containing a high concentration of TFs and co-factors. The mechanisms underlying these stages of transcriptional activation are not fully explained. This review focuses on how enhancers and promoters are activated during differentiation and how multiple enhancers work together to regulate gene expression. We illustrate the currently understood principles of how mammalian enhancers work and how they may be perturbed in enhanceropathies using expression of the α-globin gene cluster during erythropoiesis, as a model.