Pioneer transcription factors: establishing competence for gene expression

Transcriptional Regulation of Mouse Mast Cell Differentiation and the Role of Human Lung Mast Cells in Airway Inflammation

Mast cells (MCs) play a critical role in allergic inflammation, anaphylaxis, and chronic inflammatory diseases such as asthma, COPD, and osteoarthritis. Dysregulated MC activation can lead to MC activation syndrome (MACS), which is observed in patients with long COVID. MCs express the high-affinity receptor for IgE and, upon activation, release mediators and cytokines that trigger anaphylactic shock and promote allergic inflammation. They also interact with epithelial and nerve cells, which are crucial in forming a complex network of cell-cell and gene-gene interactions driving chronic inflammation that can confer resistance to treatment. In this review, in the context of the literature, we focus on experiments conducted in our laboratory investigating how transcription factors and enhancers regulate genes critical in mouse MC differentiation and function related to human lung inflammation.

The Link Between TGF-β Signaling and Chromatin Dynamics Takes an Unexpected Path

The TGF-β/BMP signaling superfamily plays critical roles in regulating a wide range of physiological and pathological processes. A recent study on Tgfbr1 function during mouse embryonic development found that, in the absence of this gene, the primordium of the external genitalia changes its fate, forming an additional set of hindlimbs and resulting in the so-called "six-legged" mice. Analysis of these embryos demonstrated that Tgfbr1 induces global remodeling of chromatin accessibility, establishing distinct regulatory elements that other morphogenetic factors use to drive either hindlimb or genital development. While global changes in chromatin structure have been previously associated with TGF-β signaling, the alterations seen in six-legged embryos differ significantly from previously described mechanisms. In this review, I discuss these distinctions and underscore the importance of studying biological processes within their natural tissue context, where essential tissue interactions are preserved.

SWI/SNF complex-mediated ZNF410 cooperative binding maintains chromatin accessibility and enhancer activity

The clustering of multiple transcription factor binding sites (TFBSs) for the same TF has proved to be a pervasive feature of cis-regulatory elements in the eukaryotic genome. However, the contribution of binding sites within the homotypic clusters of TFBSs (HCTs) to TF binding and target gene expression remains to be understood. Here, we characterize the CHD4 enhancers that harbor unique functional ZNF410 HCTs genome wide. We uncover that ZNF410 controls chromatin accessibility and activity of the CHD4 enhancer regions. We demonstrate that ZNF410 binds to the HCTs in a collaborative fashion, further conferring transcriptional activation. In particular, three ZNF410 motifs (sub-HCTs) located at 3' end of the distal enhancer act as "switch motifs" to control chromatin accessibility and enhancer activity. Mechanistically, the SWI/SNF complex is selectively required to mediate cooperative ZNF410 binding for CHD4 expression. Together, our findings expose a complex functional hierarchy of homotypic clustered motifs, which cooperate to fine-tune target gene expression.

Tumors and their microenvironments: Learning from pediatric brain pathologies

Early clues to tumors and their microenvironments come from embryonic development. Here we review the literature and consider whether the embryonic brain and its pathologies can serve as a better model. Among embryonic organs, the brain is the most heterogenous and complex, with multiple lineages leading to wide spectrum of cell states and types. Its dysregulation promotes neurodevelopmental brain pathologies and pediatric tumors. Embryonic brain pathologies point to the crucial importance of spatial heterogeneity over time, akin to the tumor microenvironment. Tumors dedifferentiate through genetic mutations and epigenetic modulations; embryonic brains differentiate through epigenetic modulations. Our innovative review proposes learning developmental brain pathologies to target tumor evolution-and vice versa. We describe ways through which tumor pharmacology can learn from embryonic brains and their pathologies, and how learning tumor, and its microenvironment, can benefit targeting neurodevelopmental pathologies. Examples include pediatric low-grade versus high-grade brain tumors as in rhabdomyosarcomas and gliomas.

Mutant CEBPA promotes tolerance to inflammatory stress through deficient AP-1 activation

The CEBPA transcription factor is frequently mutated in acute myeloid leukemia (AML). Mutations in the CEBPA gene, which are typically biallelic, result in the production of a shorter isoform known as p30. Both the canonical 42-kDa isoform (p42) and the AML-associated p30 isoform bind chromatin and activate transcription, but the specific transcriptional programs controlled by each protein and how they are linked to a selective advantage in AML is not well understood. Here, we show that cells expressing the AML-associated p30 have reduced baseline inflammatory gene expression and display altered dynamics of transcriptional induction in response to LPS, consequently impacting cytokine secretion. This confers p30-expressing cells an increased resistance to the adverse effects of prolonged exposure to inflammatory signals. Mechanistically, we show that these differences primarily arise from the differential regulation of AP-1 family proteins. In addition, we find that the impaired function of the AP-1 member ATF4 in p30-expressing cells alters their response to ER stress. Collectively, these findings uncover a link between mutant CEBPA, inflammation and the stress response, potentially revealing a vulnerability in AML.

Persistent epigenetic memory of SARS-CoV-2 mRNA vaccination in monocyte-derived macrophages

Immune memory plays a critical role in the development of durable antimicrobial immune responses. How precisely mRNA vaccines train innate immune cells to shape protective host defense mechanisms remains unknown. Here we show that SARS-CoV-2 mRNA vaccination significantly establishes histone H3 lysine 27 acetylation (H3K27ac) at promoters of human monocyte-derived macrophages, suggesting epigenetic memory. However, we found that two consecutive vaccinations were required for the persistence of H3K27ac, which matched with pro-inflammatory innate immune-associated transcriptional changes and antigen-mediated cytokine secretion. H3K27ac at promoter regions were preserved for six months and a single mRNA booster vaccine potently restored their levels and release of macrophage-derived cytokines. Interestingly, we found that H3K27ac at promoters is enriched for G-quadruplex DNA secondary structure-forming sequences in macrophage-derived nucleosome-depleted regions, linking epigenetic memory to nucleic acid structure. Collectively, these findings reveal that mRNA vaccines induce a highly dynamic and persistent training of innate immune cells enabling a sustained pro-inflammatory immune response.

Competing dynamic gene regulatory networks involved in fibroblast reprogramming to hematopoietic progenitor cells

Direct reprogramming of somatic cells offers a potentially safer therapeutic approach to generate patient-specific hematopoietic cells. However, this strategy is limited by stochasticity of reprogramming. Investigating the gene regulatory networks involved during reprogramming would help generate functional cells in adequate numbers. To address this, we developed an inducible system to reprogram fibroblasts to hematopoietic progenitor cells by ectopically expressing the two transcription factors SCL and LMO2. Transcriptome and epigenome analysis at different stages of reprogramming revealed uniform silencing of fibroblast genes and upregulation of the hemogenic endothelial program. Integrated analysis suggested that the transcription factors FLI1, GATA1/2, and KLF14 are direct targets of SCL/LMO2, which subsequently induce the hematopoietic program. Single-cell RNA sequencing revealed conflicting and competing fate decisions at intermediate stages of reprogramming. Inhibiting signaling pathways associated with competing neuronal fate enhanced reprogramming efficiency. In conclusion, this study identifies early/intermediate reprogramming events and associated pathways that could be targeted to improve reprogramming efficiency.

Gonadotropin-releasing hormone regulates transcription of the inhibin B co-receptor, TGFBR3L, via early growth response one

Follicle-stimulating hormone (FSH), a product of pituitary gonadotrope cells, regulates gonadal function and fertility. FSH production is stimulated by gonadotropin-releasing hormone (GnRH) and activin-class ligands of the TGFβ family. Inhibin A and B are TGFβ proteins that suppress FSH synthesis by competitively binding activin type II receptors in concert with the co-receptors betaglycan (TGFBR3) and TGFBR3L. Betaglycan mediates the actions of both inhibins and is broadly expressed. In contrast, TGFBR3L is inhibin B-specific and selectively expressed in gonadotropes. This cell-restricted expression is driven, in part, by steroidogenic factor 1 (SF-1, NR5A1), which stimulates Tgfbr3l/TGFBR3L transcription via two conserved promoter elements. Tgfbr3l expression is lost in mice lacking SF-1 in gonadotropes. However, SF-1 alone is unlikely to fully explain gonadotrope-restricted Tgfbr3l/TGFBR3L expression. Here, we report that GnRH induces binding of the transcription factor, early growth response 1 (EGR1), to the murine Tgfbr3l and human TGFBR3L promoters at a conserved cis-element between the two SF-1 binding sites. In homologous LβT2 cells, GnRH stimulation of Tgfbr3l/TGFBR3L promoter-reporters depends on EGR1 binding to this cis-element. In heterologous cells, over-expressed EGR1 independently and synergistically with SF-1 activates Tgfbr3l/TGFBR3L promoter-reporter activities. In vivo, Tgfbr3l mRNA expression is reduced in the pituitaries of: 1) GnRH-deficient mice, 2) wild-type mice treated with a GnRH receptor antagonist, and 3) gonadotrope-specific Egr1 knockout mice. Gonadectomy, which increases GnRH pulse frequency, enhances Tgfbr3l expression in control but not gonadotrope-specific Egr1 knockouts. Collectively, these data indicate that GnRH stimulates Tgfbr3l/TGFBR3L transcription via EGR1, which acts with SF-1 through conserved promoter elements.

The AaFoxA factor regulates female reproduction through chromatin remodeling in the mosquito vector Aedes aegypti

Female mosquitoes are vectors of many devastating human diseases because they require blood feeding to initiate reproduction. Thus, elucidation of molecular mechanisms managing female mosquito reproduction is essential. Although the regulation of gene expression during the mosquito gonadotrophic cycle has been studied in detail, how this process is controlled at the chromatin level remains unclear. Chromatin must be accessible for transcription factors (TFs) governing gene expression. A specialized class of TFs, called pioneer factors (PFs), binds and remodels closed chromatin, permitting other TFs to bind DNA and activate the gene expression. Here, we identified a homolog of the vertebrate PF FoxA in the mosquito Aedes aegypti and used the CRISPR-Cas9 system to generate mosquitoes deficient in AaFoxA. We found that ovary development was severely retarded in mutant females. Multiomics and molecular biology analyses have shown that AaFoxA increased histone acetylation and decreased methylation of H3K27 by controlling the chromatin accessibility of histone modification enzymes and chromatin remodelers. AaFoxA is bound to the loci of chromatin remodelers, changing their chromatin accessibility and modulating their temporal expression patterns. AaFoxA increased the accessibility of the ecdysone receptor (EcR) and E74 loci, indicating the important role of AaFoxA in the hormonal regulation of mosquito reproductive events. Further, the CUT&RUN and ATAC-seq analyses revealed that AaFoxA temporarily bound closed chromatin, making it differentially accessible during the mosquito gonadotrophic cycle. Hence, this study demonstrates that AaFoxA modulates chromatin dynamics throughout female mosquito reproduction.

A novel function of the Wnt antagonist secreted frizzled-related protein 4 as a transcriptional regulator of Dickkopf-1, another Wnt antagonist, in glioblastoma cell line U87MG

Wnt/β-catenin pathway dysregulation is associated with glioblastoma multiforme (GBM) pathogenesis and Wnt antagonists are downregulated in GBM. Wnt antagonist secreted frizzled-related protein 4 (sFRP4) has a tissue-specific, anti-metastatic and anti-stemness property. Our lab previously reported that gene silencing of sFRP4 in GBM cell line U87MG increases expression of another Wnt antagonist, Dickkopf-1 (Dkk1) and sFRP4 has a DNA binding ability. These findings in accordance with the nuclear localization of sFRP4 led to our present hypothesis that sFRP4 presumably negatively regulates Dkk1 and it probably interacts with the promoter region of Dkk1. Methylation-specific PCR (MSP), chromatin accessibility real-time PCR (ChART-PCR) assay, chromatin immunoprecipitation (ChIP), and quantitative DNA-protein interaction enzyme-linked immunosorbent assay (qDPI-ELISA) were carried out to test our hypothesis. We demonstrated that sFRP4 overexpression does not alter the methylation status of the Dkk1 promoter region. sFRP4 overexpression inhibits DNA-transcription factor interaction and enables chromatin accessibility to DNase I. Pertinently, sFRP4 has strong putative binding sites in the Dkk1 promoter region and its overexpression disrupts its interaction with the Dkk1 promoter. Interestingly, sFRP4 has the strongest affinity towards the -282 to +118 bp region. Downregulation of Dkk1 by overexpressed sFRP4 occurs by inhibition of the direct interaction of sFRP4 with the promoter region of Dkk1 as observed with low concentrations of sFRP4. We report for the first time a novel function of the Wnt antagonist sFRP4 acting as a transcription factor for another Wnt antagonist Dkk1, throwing open a new vista in the complex interplay between different antagonists of the Wnt pathway.

Structural insights into the recognition of native nucleosomes by pioneer transcription factors

Pioneer transcription factors possess the unique ability to bind to nucleosomal DNA and locally open closed chromatin, enabling the binding of additional chromatin-associated factors. These factors are pivotal in determining cell fate. Structural studies of pioneer transcription factors interacting with nucleosomes have predominantly relied on model systems incorporating canonical DNA motifs within synthetic, strongly positioned DNA. However, recent advances have revealed structures of several pioneer transcription factors bound to their native nucleosome targets at gene enhancers involved in cell reprogramming. These findings offer fresh insights into how pioneer transcription factors recognize and disrupt compact chromatin. In this review, we summarize these recent discoveries and explore their broader implications.

OLIG2 mediates a rare targetable stem cell fate transition in sonic hedgehog medulloblastoma

Functional cellular heterogeneity in tumours often underlies incomplete response to therapy and relapse. Previously, we demonstrated that the growth of the paediatric brain malignancy, sonic hedgehog subgroup medulloblastoma, is rooted in a dysregulated developmental hierarchy, the apex of which is defined by characteristically quiescent SOX2+ stem-like cells. Integrating gene expression and chromatin accessibility patterns in distinct cellular compartments, we identify the transcription factor Olig2 as regulating the stem cell fate transition from quiescence to activation, driving the generation of downstream neoplastic progenitors. Inactivation of Olig2 blocks stem cell activation and tumour output. Targeting this rare OLIG2-driven proliferative programme with a small molecule inhibitor, CT-179, dramatically attenuates early tumour formation and tumour regrowth post-therapy, and significantly increases median survival in vivo. We demonstrate that targeting transition from quiescence to proliferation at the level of the tumorigenic cell could be a pivotal medulloblastoma treatment strategy.

TaDL interacts with TaB3 and TaNF-YB1 to synergistically regulate the starch synthesis and grain quality in bread wheat

Starch biosynthesis is a critical factor in wheat (Triticum aestivum L.) quality and yield. However, the full scope of its regulation is not fully understood. Here we report that TaDL interacts with TaB3 and TaNF-YB1 to synergistically regulate starch biosynthesis and quality in wheat. Genome-edited tadl mutant lines had smaller and lighter grains with lower total starch and amylose contents compared to wild type (WT). Correspondingly, the transcript levels of starch biosynthesis-related genes, including TaSUS1, TaSUS2, TaAGPL2, TaSBEIIa, TaGBSSII, and TaSWEET2a, were markedly lower at 15 d after flowering (DAF) in tadl mutants. TaDL physically interacted with TaB3 and TaNF-YB1 and activated the transcription of TaSUS2 and TaAGPL2 through direct binding to their promoter regions. A null mutant of TaB3 also affected grain filling, with phenotypes similar to those of tadl mutants, whereas overexpression of TaNF-YB1 promoted grain filling. Our study demonstrated that TaDL plays an essential role in starch biosynthesis and identified an elite allele (TaDL-BI) associated with starch content, providing insights into the underlying molecular mechanism of wheat grain filling, which may be useful in breeding of high-yielding wheat and quality improvement.

The role of chromatin-related epigenetic modulations in CAKUT

Congenital anomalies of the kidney and urinary tract (CAKUT) represent a major health burden in humans. Phenotypes range from renal hypoplasia or renal agenesis, cystic renal dysplasia, duplicated or horseshoe kidneys to obstruction of the ureteropelvic junction, megaureters, duplicated ureters, urethral valves or bladder malformations. Over the past decade, next-generation sequencing has identified numerous causative genes; however, the genetic basis of most cases remains unexplained. It is assumed that environmental factors have a significant impact on the phenotype, but, overall, the pathogenesis has remained poorly understood. Interestingly however, CAKUT is a common phenotypic feature in two human syndromes, Kabuki and Koolen-de Vries syndrome, caused by dysfunction of genes encoding for KMT2D and KANSL1, both members of protein complexes playing an important role in histone modifications. In this chapter, we discuss current knowledge regarding epigenetic modulation in renal development and a putatively under-recognized role of epigenetics in CAKUT.

Transcription factor-dependent regulatory networks of sexual reproduction in Fusarium graminearum

Transcription factors (TFs) involved in sexual reproduction in filamentous fungi have been characterized. However, we have little understanding of how these TFs synergize within regulatory networks resulting in sexual development. We investigated 13 TFs in Fusarium graminearum, whose knockouts exhibited abortive or arrested phenotypes during sexual development to elucidate the transcriptional regulatory cascade underlying the development of the sexual fruiting bodies. A Bayesian network of the TFs was inferred based on transcriptomic data from key stages of sexual development. We evaluated in silico knockout impacts to the networks of the developmental phenotypes among the TFs and guided knockout transcriptomics experiments to properly assess regulatory roles of genes with same developmental phenotypes. Additional transcriptome data were collected for the TF knockouts guided by the stage at which their phenotypes appeared and by the cognate in silico prediction. Global TF networks revealed that TFs within the mating-type locus (MAT genes) trigger a transcriptional cascade involving TFs that affected early stages of sexual development. Notably, PNA1, whose knockout mutants produced exceptionally small protoperithecia, was shown to be an upstream activator for MAT genes and several TFs essential for ascospore production. In addition, knockout mutants of SUB1 produced excessive numbers of protoperithecia, wherein MAT genes and pheromone-related genes exhibited dysregulated expression. We conclude that PNA1 and SUB1 play central and suppressive roles in initiating sexual reproduction, respectively. This comprehensive investigation contributes to our understanding of the transcriptional framework governing the multicellular body plan during sexual reproduction in F. graminearum.IMPORTANCEUnderstanding transcriptional regulation of sexual development is crucial to the elucidation of the complex reproductive biology in Fusarium graminearum. We performed gene knockouts on 13 transcription factors (TFs), demonstrating knockout phenotypes affecting distinct stages of sexual development. Using transcriptomic data across stages of sexual development, we inferred a Bayesian network of these TFs that guided experiments to assess the robustness of gene interactions using a systems biology approach. We discovered that the mating-type locus (MAT genes) initiates a transcriptional cascade, with PNA1 identified as an upstream activator essential for early sexual development and ascospore production. Conversely, SUB1 was found to play a suppressive role, with knockout mutants exhibiting excessive protoperithecia due to abnormally high expression of MAT and pheromone-related genes. These findings highlight the central roles of PNA1 and SUB1 in regulating other gene activity related to sexual reproduction, contributing to a deeper understanding of the mechanisms of the multiple TFs that regulate sexual development.

ChromBPNet: bias factorized, base-resolution deep learning models of chromatin accessibility reveal cis-regulatory sequence syntax, transcription factor footprints and regulatory variants

Despite extensive mapping of cis-regulatory elements (cREs) across cellular contexts with chromatin accessibility assays, the sequence syntax and genetic variants that regulate transcription factor (TF) binding and chromatin accessibility at context-specific cREs remain elusive. We introduce ChromBPNet, a deep learning DNA sequence model of base-resolution accessibility profiles that detects, learns and deconvolves assay-specific enzyme biases from regulatory sequence determinants of accessibility, enabling robust discovery of compact TF motif lexicons, cooperative motif syntax and precision footprints across assays and sequencing depths. Extensive benchmarks show that ChromBPNet, despite its lightweight design, is competitive with much larger contemporary models at predicting variant effects on chromatin accessibility, pioneer TF binding and reporter activity across assays, cell contexts and ancestry, while providing interpretation of disrupted regulatory syntax. ChromBPNet also helps prioritize and interpret regulatory variants that influence complex traits and rare diseases, thereby providing a powerful lens to decode regulatory DNA and genetic variation.

Chromatin enables precise and scalable gene regulation with factors of limited specificity

Biophysical constraints limit the specificity with which transcription factors (TFs) can target regulatory DNA. While individual nontarget binding events may be low affinity, the sheer number of such interactions could present a challenge for gene regulation by degrading its precision or possibly leading to an erroneous induction state. Chromatin can prevent nontarget binding by rendering DNA physically inaccessible to TFs, at the cost of energy-consuming remodeling orchestrated by pioneer factors (PFs). Under what conditions and by how much can chromatin reduce regulatory errors on a global scale? We use a theoretical approach to compare two scenarios for gene regulation: one that relies on TF binding to free DNA alone and one that uses a combination of TFs and chromatin-regulating PFs to achieve desired gene expression patterns. We find, first, that chromatin effectively silences groups of genes that should be simultaneously OFF, thereby allowing more accurate graded control of expression for the remaining ON genes. Second, chromatin buffers the deleterious consequences of nontarget binding as the number of OFF genes grows, permitting a substantial expansion in regulatory complexity. Third, chromatin-based regulation productively co-opts nontarget TF binding for ON genes in order to establish a "leaky" baseline expression level, which targeted activator or repressor binding subsequently up- or down-modulates. Thus, on a global scale, using chromatin simultaneously alleviates pressure for high specificity of regulatory interactions and enables an increase in genome size with minimal impact on global expression error.

Asymmetry of Motif Conservation Within Their Homotypic Pairs Distinguishes DNA-Binding Domains of Target Transcription Factors in ChIP-Seq Data

Transcription factors (TFs) are the main regulators of eukaryotic gene expression. The cooperative binding of at least two TFs to genomic DNA is a major mechanism of transcription regulation. Massive analysis of the co-occurrence of overrepresented pairs of motifs for different target TFs studied in ChIP-seq experiments can clarify the mechanisms of TF cooperation. We categorized the target TFs from M. musculus ChIP-seq and A. thaliana ChIP-seq/DAP-seq experiments according to the structure of their DNA-binding domains (DBDs) into classes. We studied homotypic pairs of motifs, using the same recognition model for each motif. Asymmetric and symmetric pairs consist of motifs of remote and close recognition scores. We found that asymmetric pairs of motifs predominate for all TF classes. TFs from the murine/plant 'Basic helix-loop-helix (bHLH)', 'Basic leucine zipper (bZIP)', and 'Tryptophan cluster' classes and murine 'p53 domain' and 'Rel homology region' classes showed the highest enrichment of asymmetric homotypic pairs of motifs. Pioneer TFs, despite their DBD types, have a higher significance of asymmetry within homotypic pairs of motifs compared to other TFs. Asymmetry within homotypic CEs is a promising new feature decrypting the mechanisms of gene transcription regulation.

Dynamic 3D genome reorganization during senescence: defining cell states through chromatin

Cellular senescence, a cell state characterized by growth arrest and insensitivity to growth stimulatory hormones, is accompanied by a massive change in chromatin organization. Senescence can be induced by a range of physiological signals and pathological stresses and was originally thought to be an irreversible state, implicated in normal development, wound healing, tumor suppression and aging. Recently cellular senescence was shown to be reversible in some cases, with exit being triggered by the modulation of the cell's transcriptional program by the four Yamanaka factors, the suppression of p53 or H3K9me3, PDK1, and/or depletion of AP-1. Coincident with senescence reversal are changes in chromatin organization, most notably the loss of senescence-associated heterochromatin foci (SAHF) found in oncogene-induced senescence. In addition to fixed-cell imaging, chromatin conformation capture and multi-omics have been used to examine chromatin reorganization at different spatial resolutions during senescence. They identify determinants of SAHF formation and other key features that differentiate distinct types of senescence. Not surprisingly, multiple factors, including the time of induction, the type of stress experienced, and the type of cell involved, influence the global reorganization of chromatin in senescence. Here we discuss how changes in the three-dimensional organization of the genome contribute to the regulation of transcription at different stages of senescence. In particular, the distinct contributions of heterochromatin- and lamina-mediated interactions, changes in gene expression, and other cellular control mechanisms are discussed. We propose that high-resolution temporal and spatial analyses of the chromatin landscape during senescence will identify early markers of the different senescence states to help guide clinical diagnosis.

Profiling of snoRNAs in Exosomes Secreted from Cells Infected with Influenza A Virus

Small nucleolar RNAs (snoRNAs) are non-coding RNAs (ncRNAs) that regulate many cellular processes. Changes in the profiles of cellular ncRNAs and those secreted in exosomes are observed during viral infection. In our study, we analysed differences in expression profiles of snoRNAs isolated from exosomes of influenza (IAV)-infected and non-infected MDCK cells using high-throughput sequencing. The analysis revealed 133 significantly differentially regulated snoRNAs (131 upregulated and 2 downregulated), including 93 SNORD, 38 SNORA, and 2 SCARNA. The most upregulated was SNORD58 (log2FoldChange = 9.61), while the only downregulated snoRNAs were SNORD3 (log2FC = -2.98) and SNORA74 (log2FC = -2.67). Several snoRNAs previously described as involved in viral infections were upregulated, including SNORD27, SNORD28, SNORD29, SNORD58, and SNORD44. In total, 533 interactors of dysregulated snoRNAs were identified using the RNAinter database with an assigned confidence score ≥ 0.25. The main groups of predicted interactors were transcription factors (TFs, 169 interactors) and RNA-binding proteins (RBPs, 130 interactors). Among the most important were pioneer TFs such as POU5F1, SOX2, CEBPB, and MYC, while in the RBP category, notable interactors included Polr2a, TNRC6A, IGF2BP3, and FMRP. Our results suggest that snoRNAs are involved in pro-viral activity, although follow-up studies including experimental validation would be beneficial.

Phospho-regulation of ASCL1-mediated chromatin opening during cellular reprogramming

The proneural transcription factor ASCL1 regulates neurogenesis and drives somatic cell reprogramming into neurons. However, not all cell types can be reprogrammed by ASCL1, raising the questions of what provides competence and how we can overcome barriers to enable directed differentiation. Here, we investigate how levels of ASCL1 and its phosphorylation modulate its activity over progressive lineage restriction of mouse embryonic stem cells. We find that inhibition of ASCL1 phosphorylation enhances reprogramming of both mesodermal and neuroectodermal cells, while pluripotent cells remain refractory to ASCL1-directed neuronal differentiation. By performing RNA-seq and ATAC-seq in neuroectoderm, we find that un(der)phosphorylated ASCL1 causes increased chromatin accessibility at sites proximal to neuronal genes, accompanied by their increased expression. Combined analysis of protein stability and proneural function of phosphomutant and phosphomimetic ASCL1 reveals that protein stability plays only a marginal role in regulating activity, while changes in amino acid charge cannot fully explain enhanced activity of the serine-proline mutant variants of ASCL1. Our work provides new insights into proneural factor activity and regulation, and suggests ways to optimize reprogramming protocols in cancer and regenerative medicine.

Runx3, Brn3a and Isl1 interplay orchestrates the transcriptional program in the early stages of proprioceptive neuron development

BACKGROUND

The development and diversification of sensory proprioceptive neurons, which reside in the dorsal root ganglia (DRG) and express the tropomyosin receptor kinase C (TrkC), depend on the transcription factor (TF) Runx3. Runx3-deficient mice develop severe limb ataxia due to TrkC neuron cell death. Two additional TFs Pou4f1 (also called Brn3a) and Isl1 also play an important role in sensory neuron development. Thus, we aimed to unravel the chromatin state of early-developing TrkC neurons and decipher the Runx3 high-confidence target genes (HCT) and the possible cooperation between Runx3, Brn3a and Isl1 in the regulation of these genes.

METHODS

Runx3 expression is driven by the gene proximal P2 promoter. Transcriptome analysis was conducted by RNA-seq on RNA isolated from heterozygous (P2+/-) vs. homozygous (P2-/-) TrkC neurons and differentially expressed genes (DEGs) were determined. Genome-wide occupancy of Runx3, Brn3a, Isl1 and histone H3 acetylated on lysine 27 (H3K27Ac) was determined using CUT&RUN. The landscape of Transposase-accessible chromatin was analyzed via ATAC-seq.

FINDINGS

The intersection of Runx3 genomic occupancy-associated genes and DEG data discovered 244 Runx3 HCT. Brn3a and Isl1 were found to bind to numerous genomic loci, some of which overlapped with Runx3. Most genomic regions bound by each of these three TFs or co-bound by them resided in distantly located enhancer regions rather than in gene promoters. In activated and suppressed neuronal Runx3 HCT, Runx3 cooperated mainly with Brn3a to regulate expression through distantly located enhancers. Interestingly, suppression of non-neuronal immune genes was mainly managed via Runx3 without Brn3a. The distribution of ATAC and H3K27Ac marked regions in Runx3 peaks containing at least one RUNX binding site (Runx3_RBS) revealed that while most promoter regions were marked by ATAC, a prominent fraction of intron/intergenic regions occupied by Runx3, Brn3a or Isl1 were unmarked by ATAC and/or H3K27Ac.

CONCLUSIONS

These analyses shed new light on the interplay of Runx3, Brn3a, Isl1, and open chromatin regions in regulating the Runx3 HCT in the early developmental stages of TrkC neurons.

Transcriptional regulation of Cry2Ab toxin receptor ABCA2 gene in insects involves GATAe and splicing of a 5' UTR intron

Bacillus thuringiensis (Bt) produces Cry toxins that are used to control insect pests worldwide. However, evolution of insect resistance threatens the sustainable application of these toxins. In some cases, Cry toxin resistance has been linked to mutations affecting toxin receptors expression. Previous work identified HaGATAe transcriptional factor (TF) to be involved in the expression of multiple Cry1 receptor genes. Also, it was reported that 5´untranslated region (UTR) could be involved in regulation of gene expression in eukaryotic cells. The ABCA2 protein functions as Cry2A toxin receptor in multiple lepidopteran species. Here, we investigated regulation of HaABCA2 expression in Helicoverpa armigera and in different insect cell lines. Transient expression of HaABCA2 gene resulted in susceptibility to Cry2Ab in Sf9 cells. Transient expression of HaGATAe transcriptional factor in Sf9 cells enhanced the expression of multiple larval midgut proteins including SfABCA2, increasing the susceptibility to activated Cry2Ab. The silencing of HaGATAe expression in H. armigera larvae by RNAi, resulted in lower expression of HaABCA2 which correlated with reduced susceptibility to Cry2Ab. The GATAe-binding site in the promoter of HaABCA2 gene was identified by systematic truncations, site directed mutagenesis and DNA Pull-down analysis. In addition, 5' RACE analysis revealed that HaABCA2 transcripts in larval midgut cells had at least three different 5' UTRs. Here we also show that the retention of an intron in one of these 5' UTRs significantly inhibited the HaABCA2 expression. A short sequence after the start codon of translation of HaABCA2 was identified to be required for the intron removal. These findings provide new insight for mechanism of Cry2Ab resistance in H. armigera.

Single-nucleus chromatin accessibility and transcriptomic map of breast tissues of women of diverse genetic ancestry

Single-nucleus analysis allows robust cell-type classification and helps to establish relationships between chromatin accessibility and cell-type-specific gene expression. Here, using samples from 92 women of several genetic ancestries, we developed a comprehensive chromatin accessibility and gene expression atlas of the breast tissue. Integrated analysis revealed ten distinct cell types, including three major epithelial subtypes (luminal hormone sensing, luminal adaptive secretory precursor (LASP) and basal-myoepithelial), two endothelial and adipocyte subtypes, fibroblasts, T cells, and macrophages. In addition to the known cell identity genes FOXA1 (luminal hormone sensing), EHF and ELF5 (LASP), TP63 and KRT14 (basal-myoepithelial), epithelial subtypes displayed several uncharacterized markers and inferred gene regulatory networks. By integrating breast epithelial cell gene expression signatures with spatial transcriptomics, we identified gene expression and signaling differences between lobular and ductal epithelial cells and age-associated changes in signaling networks. LASP cells and fibroblasts showed genetic ancestry-dependent variability. An estrogen receptor-positive subpopulation of LASP cells with alveolar progenitor cell state was enriched in women of Indigenous American ancestry. Fibroblasts from breast tissues of women of African and European ancestry clustered differently, with accompanying gene expression differences. Collectively, these data provide a vital resource for further exploring genetic ancestry-dependent variability in healthy breast biology.

Opportunities and challenges of single-cell and spatially resolved genomics methods for neuroscience discovery

Over the past decade, single-cell genomics technologies have allowed scalable profiling of cell-type-specific features, which has substantially increased our ability to study cellular diversity and transcriptional programs in heterogeneous tissues. Yet our understanding of mechanisms of gene regulation or the rules that govern interactions between cell types is still limited. The advent of new computational pipelines and technologies, such as single-cell epigenomics and spatially resolved transcriptomics, has created opportunities to explore two new axes of biological variation: cell-intrinsic regulation of cell states and expression programs and interactions between cells. Here, we summarize the most promising and robust technologies in these areas, discuss their strengths and limitations and discuss key computational approaches for analysis of these complex datasets. We highlight how data sharing and integration, documentation, visualization and benchmarking of results contribute to transparency, reproducibility, collaboration and democratization in neuroscience, and discuss needs and opportunities for future technology development and analysis.

Nanoscale imaging of DNA-RNA identifies transcriptional plasticity at heterochromatin

The three-dimensional structure of DNA is a biophysical determinant of transcription. The density of chromatin condensation is one determinant of transcriptional output. Chromatin condensation is generally viewed as enforcing transcriptional suppression, and therefore, transcriptional output should be inversely proportional to DNA compaction. We coupled stable isotope tracers with multi-isotope imaging mass spectrometry to quantify and image nanovolumetric relationships between DNA density and newly made RNA within individual nuclei. Proliferative cell lines and cycling cells in the murine small intestine unexpectedly demonstrated no consistent relationship between DNA density and newly made RNA, even though localized examples of this phenomenon were detected at nuclear-cytoplasmic transitions. In contrast, non-dividing hepatocytes demonstrated global reduction in newly made RNA and an inverse relationship between DNA density and transcription, driven by DNA condensates at the nuclear periphery devoid of newly made RNA. Collectively, these data support an evolving model of transcriptional plasticity that extends at least to a subset of chromatin at the extreme of condensation as expected of heterochromatin.

Foxa deficiency restricts hepatitis B virus biosynthesis through epigenic silencing

In the hepatis B virus (HBV) transgenic mouse model of chronic infection, the forkhead box protein A/hepatocyte nuclear factor 3 (Foxa/HNF3) family of pioneer transcription factors are required to support postnatal viral demethylation and subsequent HBV transcription and replication. Liver-specific Foxa-deficient mice with hepatic expression of only Foxa3 do not support HBV replication but display biliary epithelial hyperplasia with bridging fibrosis. However, liver-specific Foxa-deficient mice with hepatic expression of only Foxa1 or Foxa2 also successfully restrict viral transcription and replication but display only minimal alterations in liver physiology. These observations suggest that the level of Foxa activity, rather than the combination of specific Foxa genes, is a key determinant of HBV biosynthesis. Together, these findings suggest that targeting Foxa activity could lead to HBV DNA methylation and transcriptional inactivation, resulting in the resolution of chronic HBV infections that are responsible for approximately one million deaths annually worldwide.

IMPORTANCE

The current absence of curative therapies capable of resolving chronic hepatis B virus (HBV) infection is a major clinical problem associated with considerable morbidity and mortality. The small viral genome limits molecular targets for drug development, suggesting that the identification of cellular factors essential for HBV biosynthesis may represent alternative targets for therapeutic intervention. Genetic Foxa deficiency in the neonatal liver of HBV transgenic mice leads to the transcriptional silencing of viral DNA by CpG methylation without affecting viability or displaying an obvious phenotype. Therefore, limiting liver Foxa activity therapeutically may lead to the methylation of viral covalently closed circular DNA (cccDNA), resulting in its transcriptional silencing and ultimately the resolution of chronic HBV infection.

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.

Opening and changing: mammalian SWI/SNF complexes in organ development and carcinogenesis

The switch/sucrose non-fermentable (SWI/SNF) subfamily are evolutionarily conserved, ATP-dependent chromatin-remodelling complexes that alter nucleosome position and regulate a spectrum of nuclear processes, including gene expression, DNA replication, DNA damage repair, genome stability and tumour suppression. These complexes, through their ATP-dependent chromatin remodelling, contribute to the dynamic regulation of genetic information and the maintenance of cellular processes essential for normal cellular function and overall genomic integrity. Mutations in SWI/SNF subunits are detected in 25% of human malignancies, indicating that efficient functioning of this complex is required to prevent tumourigenesis in diverse tissues. During development, SWI/SNF subunits help establish and maintain gene expression patterns essential for proper cellular identity and function, including maintenance of lineage-specific enhancers. Moreover, specific molecular signatures associated with SWI/SNF mutations, including disruption of SWI/SNF activity at enhancers, evasion of G0 cell cycle arrest, induction of cellular plasticity through pro-oncogene activation and Polycomb group (PcG) complex antagonism, are linked to the initiation and progression of carcinogenesis. Here, we review the molecular insights into the aetiology of human malignancies driven by disruption of the SWI/SNF complex and correlate these mechanisms to their developmental functions. Finally, we discuss the therapeutic potential of targeting SWI/SNF subunits in cancer.

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.

A molecular switch from tumor suppressor to oncogene in ER+ve breast cancer: Role of androgen receptor, JAK-STAT, and lineage plasticity

Cancers develop resistance to inhibitors of oncogenes mainly due to target-centric mechanisms such as mutations and splicing. While inhibitors or antagonists force targets to unnatural conformation contributing to protein instability and resistance, activating tumor suppressors may maintain the protein in an agonistic conformation to elicit sustainable growth inhibition. Due to the lack of tumor suppressor agonists, this hypothesis and the mechanisms underlying resistance are not understood. In estrogen receptor (ER)-positive breast cancer (BC), androgen receptor (AR) is a druggable tumor suppressor offering a promising avenue for this investigation. Spatial genomics suggests that the molecular portrait of AR-expressing BC cells in tumor microenvironment corresponds to better overall patient survival, clinically confirming AR's role as a tumor suppressor. Ligand activation of AR in ER-positive BC xenografts reprograms cistromes, inhibits oncogenic pathways, and promotes cellular elasticity toward a more differentiated state. Sustained AR activation results in cistrome rearrangement toward transcription factor PROP paired-like homeobox 1, transformation of AR into oncogene, and activation of the Janus kinase/signal transducer (JAK/STAT) pathway, all culminating in lineage plasticity to an aggressive resistant subtype. While the molecular profile of AR agonist-sensitive tumors corresponds to better patient survival, the profile represented in the resistant phenotype corresponds to shorter survival. Inhibition of activated oncogenes in resistant tumors reduces growth and resensitizes them to AR agonists. These findings indicate that persistent activation of a context-dependent tumor suppressor may lead to resistance through lineage plasticity-driven tumor metamorphosis. Our work provides a framework to explore the above phenomenon across multiple cancer types and underscores the importance of factoring sensitization of tumor suppressor targets while developing agonist-like drugs.

p53motifDB: integration of genomic information and tumor suppressor p53 binding motifs

The tumor suppressor gene TP53 encodes the DNA binding transcription factor p53 and is one of the most commonly mutated genes in human cancer. Tumor suppressor activity requires binding of p53 to its DNA response elements and subsequent transcriptional activation of a diverse set of target genes. Despite decades of close study, the logic underlying p53 interactions with its numerous potential genomic binding sites and target genes is not yet fully understood. Here, we present a database of DNA and chromatin-based information focused on putative p53 binding sites in the human genome to allow users to generate and test new hypotheses related to p53 activity in the genome. Users can query genomic locations based on experimentally observed p53 binding, regulatory element activity, genetic variation, evolutionary conservation, chromatin modification state, and chromatin structure. We present multiple use cases demonstrating the utility of this database for generating novel biological hypotheses, such as chromatin-based determinants of p53 binding and potential cell type-specific p53 activity. All database information is also available as a precompiled sqlite database for use in local analysis or as a Shiny web application.

Dynamic phosphorylation of FOXA1 by Aurora B guides post-mitotic gene reactivation

FOXA1 serves as a crucial pioneer transcription factor during developmental processes and plays a pivotal role as a mitotic bookmarking factor to perpetuate gene expression profiles and maintain cellular identity. During mitosis, the majority of FOXA1 dissociates from specific DNA binding sites and redistributes to non-specific binding sites; however, the regulatory mechanisms governing molecular dynamics and activity of FOXA1 remain elusive. Here, we show that mitotic kinase Aurora B specifies the different DNA binding modes of FOXA1 and guides FOXA1 biomolecular condensation in mitosis. Mechanistically, Aurora B kinase phosphorylates FOXA1 at Serine 221 (S221) to liberate the specific, but not the non-specific, DNA binding. Interestingly, the phosphorylation of S221 attenuates the FOXA1 condensation that requires specific DNA binding. Importantly, perturbation of the dynamic phosphorylation impairs accurate gene reactivation and cell proliferation, suggesting that reversible mitotic protein phosphorylation emerges as a fundamental mechanism for the spatiotemporal control of mitotic bookmarking.

Deciphering lineage-relevant gene regulatory networks during endoderm formation by InPheRNo-ChIP

Deciphering the underlying gene regulatory networks (GRNs) that govern early human embryogenesis is critical for understanding developmental mechanisms yet remains challenging due to limited sample availability and the inherent complexity of the biological processes involved. To address this, we developed InPheRNo-ChIP, a computational framework that integrates multimodal data, including RNA-seq, transcription factor (TF)-specific ChIP-seq, and phenotypic labels, to reconstruct phenotype-relevant GRNs associated with endoderm development. The core of this method is a probabilistic graphical model that models the simultaneous effect of TFs on their putative target genes to influence a particular phenotypic outcome. Unlike the majority of existing GRN inference methods that are agnostic to the phenotypic outcomes, InPheRNo-ChIP directly incorporates phenotypic information during GRN inference, enabling the distinction between lineage-specific and general regulatory interactions. We integrated data from three experimental studies and applied InPheRNo-ChIP to infer the GRN governing the differentiation of human embryonic stem cells into definitive endoderm. Benchmarking against a scRNA-seq CRISPRi study demonstrated InPheRNo-ChIP's ability to identify regulatory interactions involving endoderm markers FOXA2, SMAD2, and SOX17, outperforming other methods. This highlights the importance of incorporating the phenotypic context during network inference. Furthermore, an ablation study confirms the synergistic contribution of ChIP-seq, RNA-seq, and phenotypic data, highlighting the value of multimodal integration for accurate phenotype-relevant GRN reconstruction.

Genome access is transcription factor-specific and defined by nucleosome position

Mammalian gene expression is controlled by transcription factors (TFs) that engage sequence motifs in a chromatinized genome, where nucleosomes can restrict DNA access. Yet, how nucleosomes affect individual TFs remains unclear. Here, we measure the ability of over one hundred TF motifs to recruit TFs in a defined chromosomal locus in mouse embryonic stem cells. This identifies a set sufficient to enable the binding of TFs with diverse tissue specificities, functions, and DNA-binding domains. These chromatin-competent factors are further classified when challenged to engage motifs within a highly phased nucleosome. The pluripotency factors OCT4-SOX2 preferentially engage non-nucleosomal and entry-exit motifs, but not nucleosome-internal sites, a preference that also guides binding genome wide. By contrast, factors such as BANP, REST, or CTCF engage throughout, causing nucleosomal displacement. This supports that TFs vary widely in their sensitivity to nucleosomes and that genome access is TF specific and influenced by nucleosome position in the cell.

Enrichment of Cis-Acting Regulatory Elements in Differentially Methylated Regions Following Lipopolysaccharide Treatment of Bovine Endometrial Epithelial Cells

Endometritis is an inflammatory disease that negatively influences fertility and is common in milk-producing cows. An in vitro model for bovine endometrial inflammation was used to identify enrichment of cis-acting regulatory elements in differentially methylated regions (DMRs) in the genome of in vitro-cultured primary bovine endometrial epithelial cells (bEECs) before and after treatment with lipopolysaccharide (LPS) from E. coli, a key player in the development of endometritis. The enriched regulatory elements contain binding sites for transcription factors with established roles in inflammation and hypoxia including NFKB and Hif-1α. We further showed co-localization of certain enriched cis-acting regulatory motifs including ARNT, Hif-1α, and NRF1. Our results show an intriguing interplay between increased mRNA levels in LPS-treated bEECs of the mRNAs encoding the key transcription factors such as AHR, EGR2, and STAT1, whose binding sites were enriched in the DMRs. Our results demonstrate an extraordinary cis-regulatory complexity in these DMRs having binding sites for both inflammatory and hypoxia-dependent transcription factors. Obtained data using this in vitro model for bacterial-induced endometrial inflammation have provided valuable information regarding key transcription factors relevant for clinical endometritis in both cattle and humans.

Structural insights into the cooperative nucleosome recognition and chromatin opening by FOXA1 and GATA4

Mouse FOXA1 and GATA4 are prototypes of pioneer factors, initiating liver cell development by binding to the N1 nucleosome in the enhancer of the ALB1 gene. Using cryoelectron microscopy (cryo-EM), we determined the structures of the free N1 nucleosome and its complexes with FOXA1 and GATA4, both individually and in combination. We found that the DNA-binding domains of FOXA1 and GATA4 mainly recognize the linker DNA and an internal site in the nucleosome, respectively, whereas their intrinsically disordered regions interact with the acidic patch on histone H2A-H2B. FOXA1 efficiently enhances GATA4 binding by repositioning the N1 nucleosome. In vivo DNA editing and bioinformatics analyses suggest that the co-binding mode of FOXA1 and GATA4 plays important roles in regulating genes involved in liver cell functions. Our results reveal the mechanism whereby FOXA1 and GATA4 cooperatively bind to the nucleosome through nucleosome repositioning, opening chromatin by bending linker DNA and obstructing nucleosome packing.

Systematic dissection of sequence features affecting binding specificity of a pioneer factor reveals binding synergy between FOXA1 and AP-1

Despite the unique ability of pioneer 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 chromatin immunoprecipitation with integrated synthetic oligonucleotides (ChIP-ISO) to systematically dissect sequence features affecting the binding specificity of a classic PF, FOXA1, in human A549 cells. Combining ChIP-ISO with in vitro and neural network analyses, we find that (1) FOXA1 binding is strongly affected by co-binding transcription factors (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.

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.

Temporally resolved early bone morphogenetic protein-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 bone morphogenetic protein (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 hr 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.

How Transcription Factor Clusters Shape the Transcriptional Landscape

In eukaryotic cells, gene transcription typically occurs in discrete periods of promoter activity, interspersed with intervals of inactivity. This pattern deviates from simple stochastic events and warrants a closer examination of the molecular interactions that activate the promoter. Recent studies have identified transcription factor (TF) clusters as key precursors to transcriptional bursting. Often, these TF clusters form at chromatin segments that are physically distant from the promoter, making changes in chromatin conformation crucial for promoter-TF cluster interactions. In this review, I explore the formation and constituents of TF clusters, examining how the dynamic interplay between chromatin architecture and TF clustering influences transcriptional bursting. Additionally, I discuss techniques for visualizing TF clusters and provide an outlook on understanding the remaining gaps in this field.

Data enhancement in the age of spatial biology

Unveiling the intricate interplay of cells in their native environment lies at the heart of understanding fundamental biological processes and unraveling disease mechanisms, particularly in complex diseases like cancer. Spatial transcriptomics (ST) offers a revolutionary lens into the spatial organization of gene expression within tissues, empowering researchers to study both cell heterogeneity and microenvironments in health and disease. However, current ST technologies often face limitations in either resolution or the number of genes profiled simultaneously. Integrating ST data with complementary sources, such as single-cell transcriptomics and detailed tissue staining images, presents a powerful solution to overcome these limitations. This review delves into the computational approaches driving the integration of spatial transcriptomics with other data types. By illuminating the key challenges and outlining the current algorithmic solutions, we aim to highlight the immense potential of these methods to revolutionize our understanding of cancer biology.

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.

Structural dynamics in chromatin unraveling by pioneer transcription factors

Pioneer transcription factors are proteins with a dual function. First, they regulate transcription by binding to nucleosome-free DNA regulatory elements. Second, they bind to DNA while wrapped around histone proteins in the chromatin and mediate chromatin opening. The molecular mechanisms that connect the two functions are yet to be discovered. In recent years, pioneer factors received increased attention mainly because of their crucial role in promoting cell fate transitions that could be used for regenerative therapies. For example, the three factors required to induce pluripotency in somatic cells, Oct4, Sox2, and Klf4 were classified as pioneer factors and studied extensively. With this increased attention, several structures of complexes between pioneer factors and chromatin structural units (nucleosomes) have been resolved experimentally. Furthermore, experimental and computational approaches have been designed to study two unresolved, key scientific questions: First, do pioneer factors induce directly local opening of nucleosomes and chromatin fibers upon binding? And second, how do the unstructured tails of the histones impact the structural dynamics involved in such conformational transitions? Here we review the current knowledge about transcription factor-induced nucleosome dynamics and the role of the histone tails in this process. We discuss what is needed to bridge the gap between the static views obtained from the experimental structures and the key structural dynamic events in chromatin opening. Finally, we propose that integrating nuclear magnetic resonance spectroscopy with molecular dynamics simulations is a powerful approach to studying pioneer factor-mediated dynamics of nucleosomes and perhaps small chromatin fibers using native DNA sequences.

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.