Combinatorial bZIP dimers display complex DNA-binding specificity landscapes

A large-scale cancer-specific protein-DNA interaction network

Cancer development and progression are generally associated with gene dysregulation, often resulting from changes in the transcription factor (TF) sequence or expression. Identifying key TFs involved in cancer gene regulation provides a framework for potential new therapeutics. This study presents a large-scale cancer gene TF-DNA interaction network, as well as an extensive promoter clone resource for future studies. Highly connected TFs bind to promoters of genes associated with either good or poor cancer prognosis, suggesting that strategies aimed at shifting gene expression balance between these two prognostic groups may be inherently complex. However, we identified potential for oncogene-targeted therapeutics, with half of the tested oncogenes being potentially repressed by influencing specific activators or bifunctional TFs. Finally, we investigate the role of intrinsically disordered regions within the key cancer-related TF ESR1 in DNA binding and transcriptional activity, and found that these regions can have complex trade-offs in TF function. Altogether, our study broadens our knowledge of the TFs involved in cancer gene regulation and provides a valuable resource for future studies and therapeutics.

Reversal of high-glucose-induced transcriptional and epigenetic memories through NRF2 pathway activation

Diabetes complications such as nephropathy, retinopathy, or cardiovascular disease arise from vascular dysfunction. In this context, it has been observed that past hyperglycemic events can induce long-lasting alterations, a phenomenon termed "metabolic memory." In this study, we evaluated the genome-wide gene expression and chromatin accessibility alterations caused by transient high-glucose exposure in human endothelial cells (ECs) in vitro. We found that cells exposed to high glucose exhibited substantial gene expression changes in pathways known to be impaired in diabetes, many of which persist after glucose normalization. Chromatin accessibility analysis also revealed that transient hyperglycemia induces persistent alterations, mainly in non-promoter regions identified as enhancers with neighboring genes showing lasting alterations. Notably, activation of the NRF2 pathway through NRF2 overexpression or supplementation with the plant-derived compound sulforaphane, effectively reverses the glucose-induced transcriptional and chromatin accessibility memories in ECs. These findings underscore the enduring impact of transient hyperglycemia on ECs' transcriptomic and chromatin accessibility profiles, emphasizing the potential utility of pharmacological NRF2 pathway activation in mitigating and reversing the high-glucose-induced transcriptional and epigenetic alterations.

Dimerization of hub protein DYNLL1 and bZIP transcription factor CREB3L1 enhances transcriptional activation of CREB3L1 target genes like arginine vasopressin

bZIP transcription factors can function as homodimers or heterodimers through interactions with their disordered coiled-coil domain. Such dimer assemblies are known to influence DNA-binding specificity and/or the recruitment of binding partners, which can cause a functional switch of a transcription factor from being an activator to a repressor. We recently identified the genomic targets of a bZIP transcription factor called CREB3L1 in rat hypothalamic supraoptic nucleus by ChIP-seq. The objective of this study was to investigate the CREB3L1 protein-to-protein interactome of which little is known. For this approach, we created and screened a rat supraoptic nucleus yeast two-hybrid prey library with the bZIP region of rat CREB3L1 as the bait. Our yeast two-hybrid approach captured five putative CREB3L1 interacting prey proteins in the supraoptic nucleus. One interactor was selected by bioinformatic analyses for more detailed investigation by co-immunoprecipitation, immunofluorescent cellular localisation, and reporter assays in vitro. Here we identify dimerisation hub protein Dynein Light Chain LC8-Type 1 as a CREB3L1 interacting protein that in vitro enhances CREB3L1 activation of target genes.

The basic leucine zipper domain (bZIP) transcription factor BbYap1 promotes evasion of host humoral immunity and regulates lipid homeostasis contributing to fungal virulence in Beauveria bassiana

Basic leucine zipper domain transcription factors (TFs), of which yeast activator protein (Yap) is a significant class, are crucial for the development of sclerotia, the stress response, vegetative growth, and spore adhesion. Nevertheless, nothing is known about how Yap TFs contribute to the pathogenicity of entomopathogenic fungus. In this work, Beauveria bassiana was used to identify and knock out the yeast gene BbYap1, which is similar to Yap. The BbYap1 gene deletion has an impact on lipid homeostasis of B. bassiana; oleic acid, for example, dropped by 95.69%. The BbYap1 mutant exhibited much less virulence and vegetative development in comparison to the wild strain, while demonstrating a greater sensitivity to chemical stress. It is noteworthy that the physiological abnormalities brought on by BbYap1 deletion were largely repaired by the addition of exogenous oleic acid, as seen by the notable increase in insect survival in the blood cavity injection group. Following infection with the BbYap1 mutant, the host exhibits a considerable down-regulation of the expression of β-1,3-glucan recognition protein, gallerimycin, gloverin, and moricin-like protein genes. Likewise, the introduction of exogenous oleic acid markedly increased the host's expression of the aforementioned genes. In summary, BbYap1 regulates cellular enzyme lipid homeostasis and fungal virulence by eluding host humoral defense, which contributes to fungal chemical stress and vegetative development.

IMPORTANCE

Entomopathogenic fungi (EPF) offer an effective and eco-friendly alternative to curb insect populations in biocontrol strategy. When EPF enter the hemolymph of their host, they encounter a variety of stress reactions, such as immunological and oxidative stress. Basic leucine zipper domain transcription factors, of which yeast activator protein (Yap) is a significant class, have diverse biological functions related to metabolism, development, reproduction, conidiation, stress responses, and pathogenicity. This study demonstrates that BbYap1 of Beauveria bassiana regulates cellular enzyme lipid homeostasis and fungal virulence by eluding host humoral defense, which contributes to fungal chemical stress and vegetative development. These findings offer fresh perspectives for comprehending molecular roles of YAP in EPF.

Arabidopsis AGB1 participates in salinity response through bZIP17-mediated unfolded protein response

BACKGROUND

Plant heterotrimeric G proteins respond to various environmental stresses, including high salinity. It is known that Gβ subunit AGB1 functions in maintaining local and systemic Na+/K+ homeostasis to accommodate ionic toxicity under salt stress. However, whether AGB1 contributes to regulating gene expression for seedling's survival under high salinity remains unclear.

RESULTS

We showed that AGB1-Venus localized to nuclei when facing excessive salt, and the induction of a set of bZIP17-dependent salt stress-responsive genes was reduced in the agb1 mutant. We confirmed both genetic and physical interactions of AGB1 and bZIP17 in plant salinity response by comparing salt responses in the single and double mutants of agb1 and bzip17 and by BiFC assay, respectively. In addition, we show that AGB1 depletion decreases nuclei-localization of transgenic mRFP-bZIP17 under salt stress, as shown in s1p s2p double mutant in the Agrobacteria-mediated transient mRFP-bZIP17 expression in young seedlings.

CONCLUSIONS

Our results indicate that AGB1 functions in S1P and/or S2P-mediated proteolytic processing of bZIP17 under salt stress to regulate the induction of salinity-responsive gene expression.

A Knockout of the Photoreceptor PtAureo1a Results in Altered Diel Expression of Diatom Clock Components

Plants and algae use light not only for driving photosynthesis but also to sense environmental cues and to adjust their circadian clocks via photoreceptors. Aureochromes are blue-light-dependent photoreceptors that also function as transcription factors, possessing both a LOV and a bZIP domain. Aureochromes so far have only been detected in Stramenopile algae, which include the diatoms. Four paralogues of aureochromes have been identified in the pennate model diatom Phaeodactylum tricornutum: PtAureo1a, 1b, 1c, and 2. While it was shown recently that diatoms have a diel rhythm, the molecular mechanisms and components regulating it are still largely unknown. Diel gene expression analyses of wild-type P. tricornutum, a PtAureo1a knockout strain, and the respective PtAureo1 complemented line revealed that all four aureochromes have a different diel regulation and that PtAureo1a has a strong co-regulatory influence on its own transcription, as well as on that of other genes encoding different blue-light photoreceptors (CPF1, 2 and 4), proteins involved in photoprotection (Lhcx1), and specific bHLH transcription factors (RITMO1). Some of these genes completely lost their circadian expression in the PtAureo1a KO mutant. Our results suggest a major involvement of aureochromes in the molecular clock of diatoms.

DeepReg: a deep learning hybrid model for predicting transcription factors in eukaryotic and prokaryotic genomes

Deep learning models (DLMs) have gained importance in predicting, detecting, translating, and classifying a diversity of inputs. In bioinformatics, DLMs have been used to predict protein structures, transcription factor-binding sites, and promoters. In this work, we propose a hybrid model to identify transcription factors (TFs) among prokaryotic and eukaryotic protein sequences, named Deep Regulation (DeepReg) model. Two architectures were used in the DL model: a convolutional neural network (CNN), and a bidirectional long-short-term memory (BiLSTM). DeepReg reached a precision of 0.99, a recall of 0.97, and an F1-score of 0.98. The quality of our predictions, the bias-variance trade-off approach, and the characterization of new TF predictions were evaluated and compared against those produced by DeepTFactor, as well as against experimental data from three model organisms. Predictions based on our DLM tended to exhibit less variance and bias than those from DeepTFactor, thus increasing reliability and decreasing overfitting.

AIRE relies on Z-DNA to flag gene targets for thymic T cell tolerization

AIRE is an unconventional transcription factor that enhances the expression of thousands of genes in medullary thymic epithelial cells and promotes clonal deletion or phenotypic diversion of self-reactive T cells1-4. The biological logic of AIRE's target specificity remains largely unclear as, in contrast to many transcription factors, it does not bind to a particular DNA sequence motif. Here we implemented two orthogonal approaches to investigate AIRE's cis-regulatory mechanisms: construction of a convolutional neural network and leveraging natural genetic variation through analysis of F1 hybrid mice5. Both approaches nominated Z-DNA and NFE2-MAF as putative positive influences on AIRE's target choices. Genome-wide mapping studies revealed that Z-DNA-forming and NFE2L2-binding motifs were positively associated with the inherent ability of a gene's promoter to generate DNA double-stranded breaks, and promoters showing strong double-stranded break generation were more likely to enter a poised state with accessible chromatin and already-assembled transcriptional machinery. Consequently, AIRE preferentially targets genes with poised promoters. We propose a model in which Z-DNA anchors the AIRE-mediated transcriptional program by enhancing double-stranded break generation and promoter poising. Beyond resolving a long-standing mechanistic conundrum, these findings suggest routes for manipulating T cell tolerance.

Identification of promoter targets by Aureochrome 1a in the diatom Phaeodactylum tricornutum

Aureochromes (AUREOs) are unique blue light receptors and transcription factors found only in stramenopile algae. While each of the four AUREOs identified in the diatom Phaeodactylum tricornutum may have a specific function, PtAUREO1a has been shown to have a strong impact on overall gene regulation, when light changes from red to blue light conditions. Despite its significance, the molecular mechanism of PtAUREO1a is largely unexplored. To comprehend the overall process of gene regulation by PtAUREO1a, we conducted a series of in vitro and in vivo experiments, including pull-down assays, yeast one-hybrid experiments, and phenotypical characterization using recombinant PtAUREOs and diatom mutant lines expressing a modified PtAureo1a gene. We describe the distinct light absorption properties of four PtAUREOs and the formation of all combinations of their potential dimers. We demonstrate the capability of PtAUREO1a and 1b to activate the genes, diatom-specific cyclin 2, PtAureo1a, and PtAureo1c under both light and dark conditions. Using mutant lines expressing a modified PtAUREO1a protein with a considerably reduced light absorption, we found novel evidence that PtAUREO1a regulates the expression of PtLHCF15, which is essential for red light acclimation. Based on current knowledge, we present a working model of PtAUREO1a gene regulation properties.

ATF4 Signaling in HIV-1 Infection: Viral Subversion of a Stress Response Transcription Factor

Cellular integrated stress response (ISR), the mitochondrial unfolded protein response (UPRmt), and IFN signaling are associated with viral infections. Activating transcription factor 4 (ATF4) plays a pivotal role in these pathways and controls the expression of many genes involved in redox processes, amino acid metabolism, protein misfolding, autophagy, and apoptosis. The precise role of ATF4 during viral infection is unclear and depends on cell hosts, viral agents, and models. Furthermore, ATF4 signaling can be hijacked by pathogens to favor viral infection and replication. In this review, we summarize the ATF4-mediated signaling pathways in response to viral infections, focusing on human immunodeficiency virus 1 (HIV-1). We examine the consequences of ATF4 activation for HIV-1 replication and reactivation. The role of ATF4 in autophagy and apoptosis is explored as in the context of HIV-1 infection programmed cell deaths contribute to the depletion of CD4 T cells. Furthermore, ATF4 can also participate in the establishment of innate and adaptive immunity that is essential for the host to control viral infections. We finally discuss the putative role of the ATF4 paralogue, named ATF5, in HIV-1 infection. This review underlines the role of ATF4 at the crossroads of multiple processes reflecting host-pathogen interactions.

Co-translational Assembly Pathways of Nuclear Multiprotein Complexes Involved in the Regulation of Gene Transcription

Most factors that regulate gene transcription in eukaryotic cells are multimeric, often large, protein complexes. The understanding of the biogenesis pathways of such large and heterogeneous protein assemblies, as well as the dimerization partner choice among transcription factors, is crucial to interpret and control gene expression programs and consequent cell fate decisions. Co-translational assembly (Co-TA) is thought to play key roles in the biogenesis of protein complexes by directing complex formation during protein synthesis. In this review we discuss the principles of Co-TA with a special focus for the assembly of transcription regulatory complexes. We outline the expected molecular advantages of establishing co-translational interactions, pointing at the available, or missing, evidence for each of them. We hypothesize different molecular mechanisms based on Co-TA to explain the allocation "dilemma" of paralog proteins and subunits shared by different transcription complexes. By taking as a paradigm the different assembly pathways employed by three related transcription regulatory complexes (TFIID, SAGA and ATAC), we discuss alternative Co-TA strategies for nuclear multiprotein complexes and the widespread - yet specific - use of Co-TA for the formation of nuclear complexes involved in gene transcription. Ultimately, we outlined a series of open questions which demand well-defined lines of research to investigate the principles of gene regulation that rely on the coordinated assembly of protein complexes.

bZIP transcription factor FabR: Redox-dependent mechanism controlling docosahexaenoic acid biosynthesis and H2O2 stress response in Schizochytrium sp

Schizochytrium sp. is an important industrial strain for commercial production of docosahexaenoic acid (DHA), which plays essential physiological roles in infant development and human health. The regulatory network for DHA biosynthesis and lipid accumulation in Schizochytrium remains poorly understood. FabR (fatty acid biosynthesis repressor), a basic leucine zipper (bZIP) transcription factor, was transcriptionally downregulated under low-nitrogen condition. Deletion of fabR gene (mutant ΔfabR) increased production of total lipids and DHA by 30.1% and 46.5%, respectively. ΔfabR displayed H2O2 stress resistance higher than that of parental strain or complementation strain CfabR. FabR bound specifically to 7-bp pseudo-palindromic sequence 5'-ATTSAAT-3' in upstream regions and repressed transcription of fatty acid biosynthesis genes (acl, fas, pfa) and antioxidant defense genes (cat, sod1, sod2, gpx). DNA binding activity of FabR was regulated in a redox-dependent manner. Under oxidative condition, FabR forms intermolecular disulfide bonds between two Cys46 residues of dimers; its DNA binding activity is thereby lost, and the transcription of its target genes is enhanced through derepression. Our findings clarify the redox-dependent mechanism that modulates FabR activity governing lipid and DHA biosynthesis and H2O2 stress response in Schizochytrium.

Integrated Stress Response (ISR) Pathway: Unraveling Its Role in Cellular Senescence

Cellular senescence is a complex process characterized by irreversible cell cycle arrest. Senescent cells accumulate with age, promoting disease development, yet the absence of specific markers hampers the development of selective anti-senescence drugs. The integrated stress response (ISR), an evolutionarily highly conserved signaling network activated in response to stress, globally downregulates protein translation while initiating the translation of specific protein sets including transcription factors. We propose that ISR signaling plays a central role in controlling senescence, given that senescence is considered a form of cellular stress. Exploring the intricate relationship between the ISR pathway and cellular senescence, we emphasize its potential as a regulatory mechanism in senescence and cellular metabolism. The ISR emerges as a master regulator of cellular metabolism during stress, activating autophagy and the mitochondrial unfolded protein response, crucial for maintaining mitochondrial quality and efficiency. Our review comprehensively examines ISR molecular mechanisms, focusing on ATF4-interacting partners, ISR modulators, and their impact on senescence-related conditions. By shedding light on the intricate relationship between ISR and cellular senescence, we aim to inspire future research directions and advance the development of targeted anti-senescence therapies based on ISR modulation.

Principles of Computation by Competitive Protein Dimerization Networks

Many biological signaling pathways employ proteins that competitively dimerize in diverse combinations. These dimerization networks can perform biochemical computations, in which the concentrations of monomers (inputs) determine the concentrations of dimers (outputs). Despite their prevalence, little is known about the range of input-output computations that dimerization networks can perform (their "expressivity") and how it depends on network size and connectivity. Using a systematic computational approach, we demonstrate that even small dimerization networks (3-6 monomers) are expressive, performing diverse multi-input computations. Further, dimerization networks are versatile, performing different computations when their protein components are expressed at different levels, such as in different cell types. Remarkably, individual networks with random interaction affinities, when large enough (≥8 proteins), can perform nearly all (~90%) potential one-input network computations merely by tuning their monomer expression levels. Thus, even the simple process of competitive dimerization provides a powerful architecture for multi-input, cell-type-specific signal processing.

Paired yeast one-hybrid assays to detect DNA-binding cooperativity and antagonism across transcription factors

Cooperativity and antagonism between transcription factors (TFs) can drastically modify their binding to regulatory DNA elements. While mapping these relationships between TFs is important for understanding their context-specific functions, existing approaches either rely on DNA binding motif predictions, interrogate one TF at a time, or study individual TFs in parallel. Here, we introduce paired yeast one-hybrid (pY1H) assays to detect cooperativity and antagonism across hundreds of TF-pairs at DNA regions of interest. We provide evidence that a wide variety of TFs are subject to modulation by other TFs in a DNA region-specific manner. We also demonstrate that TF-TF relationships are often affected by alternative isoform usage and identify cooperativity and antagonism between human TFs and viral proteins from human papillomaviruses, Epstein-Barr virus, and other viruses. Altogether, pY1H assays provide a broadly applicable framework to study how different functional relationships affect protein occupancy at regulatory DNA regions.

Genome-wide census of ATF4 binding sites and functional profiling of trait-associated genetic variants overlapping ATF4 binding motifs

Activating Transcription Factor 4 (ATF4) is an important regulator of gene expression in stress responses and developmental processes in many cell types. Here, we catalogued ATF4 binding sites in the human genome and identified overlaps with trait-associated genetic variants. We probed these genetic variants for allelic regulatory activity using a massively parallel reporter assay (MPRA) in HepG2 hepatoma cells exposed to tunicamycin to induce endoplasmic reticulum stress and ATF4 upregulation. The results revealed that in the majority of cases, the MPRA allelic activity of these SNPs was in agreement with the nucleotide preference seen in the ATF4 binding motif from ChIP-Seq. Luciferase and electrophoretic mobility shift assays in additional cellular models further confirmed ATF4-dependent regulatory effects for the SNPs rs532446 (GADD45A intronic; linked to hematological parameters), rs7011846 (LPL upstream; myocardial infarction), rs2718215 (diastolic blood pressure), rs281758 (psychiatric disorders) and rs6491544 (educational attainment). CRISPR-Cas9 disruption and/or deletion of the regulatory elements harboring rs532446 and rs7011846 led to the downregulation of GADD45A and LPL, respectively. Thus, these SNPs could represent examples of GWAS genetic variants that affect gene expression by altering ATF4-mediated transcriptional activation.

Structural basis for specific DNA sequence motif recognition by the TFAP2 transcription factors

The TFAP2 family regulates gene expression during differentiation, development, and organogenesis, and includes five homologs in humans. They all possess a highly conserved DNA binding domain (DBD) followed by a helix-span-helix (HSH) domain. The DBD-HSH tandem domain specifically binds to a GCC(N3)GGC consensus sequence, but the precise recognition mechanisms remain unclear. Here, we found that TFAP2 preferred binding to the GCC(N3)GGC sequence, and the pseudo-palindromic GCC and GGC motifs and the length of the central spacer between the two motifs determined their binding specificity. Structural studies revealed that the two flat amphipathic α-helical HSH domains of TFAP2A stacked with each other to form a dimer via hydrophobic interactions, while the stabilized loops from both DBD domains inserted into two neighboring major grooves of the DNA duplex to form base-specific interactions. This specific DNA binding mechanism controlled the length of the central spacer and determined the DNA sequence specificity of TFAP2. Mutations of the TFAP2 proteins are implicated in various diseases. We illustrated that reduction or disruption of the DNA binding ability of the TFAP2 proteins is the primary cause of TFAP2 mutation-associated diseases. Thus, our findings also offer valuable insights into the pathogenesis of disease-associated mutations in TFAP2 proteins.

Tomatidine targets ATF4-dependent signaling and induces ferroptosis to limit pancreatic cancer progression

Pancreatic ductal adenocarcinoma (PDAC) is an aggressive cancer with high metastasis and therapeutic resistance. Activating transcription factor 4 (ATF4), a master regulator of cellular stress, is exploited by cancer cells to survive. Prior research and data reported provide evidence that high ATF4 expression correlates with worse overall survival in PDAC. Tomatidine, a natural steroidal alkaloid, is associated with inhibition of ATF4 signaling in multiple diseases. Here, we discovered that in vitro and in vivo tomatidine treatment of PDAC cells inhibits tumor growth. Tomatidine inhibited nuclear translocation of ATF4 and reduced the transcriptional binding of ATF4 with downstream promoters. Tomatidine enhanced gemcitabine chemosensitivity in 3D ECM-hydrogels and in vivo. Tomatidine treatment was associated with induction of ferroptosis signaling validated by increased lipid peroxidation, mitochondrial biogenesis, and decreased GPX4 expression in PDAC cells. This study highlights a possible therapeutic approach utilizing a plant-derived metabolite, tomatidine, to target ATF4 activity in PDAC.

Prediction of cooperative homeodomain DNA binding sites from high-throughput-SELEX data

Homeodomain proteins constitute one of the largest families of metazoan transcription factors. Genetic studies have demonstrated that homeodomain proteins regulate many developmental processes. Yet, biochemical data reveal that most bind highly similar DNA sequences. Defining how homeodomain proteins achieve DNA binding specificity has therefore been a long-standing goal. Here, we developed a novel computational approach to predict cooperative dimeric binding of homeodomain proteins using High-Throughput (HT) SELEX data. Importantly, we found that 15 of 88 homeodomain factors form cooperative homodimer complexes on DNA sites with precise spacing requirements. Approximately one third of the paired-like homeodomain proteins cooperatively bind palindromic sequences spaced 3 bp apart, whereas other homeodomain proteins cooperatively bind sites with distinct orientation and spacing requirements. Combining structural models of a paired-like factor with our cooperativity predictions identified key amino acid differences that help differentiate between cooperative and non-cooperative factors. Finally, we confirmed predicted cooperative dimer sites in vivo using available genomic data for a subset of factors. These findings demonstrate how HT-SELEX data can be computationally mined to predict cooperativity. In addition, the binding site spacing requirements of select homeodomain proteins provide a mechanism by which seemingly similar AT-rich DNA sequences can preferentially recruit specific homeodomain factors.

The computational capabilities of many-to-many protein interaction networks

Many biological circuits comprise sets of protein variants that interact with one another in a many-to-many, or promiscuous, fashion. These architectures can provide powerful computational capabilities that are especially critical in multicellular organisms. Understanding the principles of biochemical computations in these circuits could allow more precise control of cellular behaviors. However, these systems are inherently difficult to analyze, due to their large number of interacting molecular components, partial redundancies, and cell context dependence. Here, we discuss recent experimental and theoretical advances that are beginning to reveal how promiscuous circuits compute, what roles those computations play in natural biological contexts, and how promiscuous architectures can be applied for the design of synthetic multicellular behaviors.

A genome-wide association analysis of loss of ambulation in dystrophinopathy patients suggests multiple candidate modifiers of disease severity

The major determinant of disease severity in Duchenne muscular dystrophy (DMD) or milder Becker muscular dystrophy (BMD) is whether the dystrophin gene (DMD) mutation truncates the mRNA reading frame or allows expression of a partially functional protein. However, even in the complete absence of dystrophin, variability in disease severity is observed, and candidate gene studies have implicated several genes as modifiers. Here we present the largest genome-wide search to date for loci influencing severity in N = 419 DMD patients. Availability of subjects for such studies is quite limited, leading to modest sample sizes, which present a challenge for GWAS design. We have therefore taken special steps to minimize heterogeneity within our dataset at the DMD locus itself, taking a novel approach to mutation classification to effectively exclude the possibility of residual dystrophin expression, and utilized statistical methods that are well adapted to smaller sample sizes, including the use of a novel linear regression-like residual for time to ambulatory loss and the application of evidential statistics for the GWAS approach. Finally, we applied an unbiased in silico pipeline, utilizing functional genomic datasets to explore the potential impact of the best supported SNPs. In all, we obtained eight SNPs (out of 1,385,356 total) with posterior probability of trait-marker association (PPLD) ≥ 0.4, representing six distinct loci. Our analysis prioritized likely non-coding SNP regulatory effects on six genes (ETAA1, PARD6G, GALNTL6, MAN1A1, ADAMTS19, and NCALD), each with plausibility as a DMD modifier. These results support both recurrent and potentially new pathways for intervention in the dystrophinopathies.

Double DAP-seq uncovered synergistic DNA binding of interacting bZIP transcription factors

Many eukaryotic transcription factors (TF) form homodimer or heterodimer complexes to regulate gene expression. Dimerization of BASIC LEUCINE ZIPPER (bZIP) TFs are critical for their functions, but the molecular mechanism underlying the DNA binding and functional specificity of homo- versus heterodimers remains elusive. To address this gap, we present the double DNA Affinity Purification-sequencing (dDAP-seq) technique that maps heterodimer binding sites on endogenous genomic DNA. Using dDAP-seq we profile twenty pairs of C/S1 bZIP heterodimers and S1 homodimers in Arabidopsis and show that heterodimerization significantly expands the DNA binding preferences of these TFs. Analysis of dDAP-seq binding sites reveals the function of bZIP9 in abscisic acid response and the role of bZIP53 heterodimer-specific binding in seed maturation. The C/S1 heterodimers show distinct preferences for the ACGT elements recognized by plant bZIPs and motifs resembling the yeast GCN4 cis-elements. This study demonstrates the potential of dDAP-seq in deciphering the DNA binding specificities of interacting TFs that are key for combinatorial gene regulation.

Activating Transcription Factor 3 Stimulates Follicle-Stimulating Hormone-β Expression In Vitro But Is Dispensable for Follicle-Stimulating Hormone Production in Murine Gonadotropes In Vivo

Follicle-stimulating hormone (FSH), a dimeric glycoprotein produced by pituitary gonadotrope cells, regulates spermatogenesis in males and ovarian follicle growth in females. Hypothalamic gonadotropin-releasing hormone (GnRH) stimulates FSHβ subunit gene (Fshb) transcription, though the underlying mechanisms are poorly understood. To address this gap in knowledge, we examined changes in pituitary gene expression in GnRH-deficient mice (hpg) treated with a regimen of exogenous GnRH that increases pituitary Fshb but not luteinizing hormone β (Lhb) messenger RNA levels. Activating transcription factor 3 (Atf3) was among the most upregulated genes. Activating transcription factor 3 (ATF3) can heterodimerize with members of the activator protein 1 family to regulate gene transcription. Co-expression of ATF3 with JunB stimulated murine Fshb, but not Lhb, promoter-reporter activity in homologous LβT2b cells. ATF3 also synergized with a constitutively active activin type I receptor to increase endogenous Fshb expression in these cells. Nevertheless, FSH production was intact in gonadotrope-specific Atf3 knockout [conditional knockout (cKO)] mice. Ovarian follicle development, ovulation, and litter sizes were equivalent between cKOs and controls. Testis weights and sperm counts did not differ between genotypes. Following gonadectomy, increases in LH secretion were enhanced in cKO animals. Though FSH levels did not differ between genotypes, post-gonadectomy increases in pituitary Fshb and gonadotropin α subunit expression were more pronounced in cKO than control mice. These data indicate that ATF3 can selectively stimulate Fshb expression in vitro but is not required for FSH production in vivo.

CREB3L2-ATF4 heterodimerization defines a transcriptional hub of Alzheimer's disease gene expression linked to neuropathology

Gene expression is changed by disease, but how these molecular responses arise and contribute to pathophysiology remains less understood. We discover that β-amyloid, a trigger of Alzheimer's disease (AD), promotes the formation of pathological CREB3L2-ATF4 transcription factor heterodimers in neurons. Through a multilevel approach based on AD datasets and a novel chemogenetic method that resolves the genomic binding profile of dimeric transcription factors (ChIPmera), we find that CREB3L2-ATF4 activates a transcription network that interacts with roughly half of the genes differentially expressed in AD, including subsets associated with β-amyloid and tau neuropathologies. CREB3L2-ATF4 activation drives tau hyperphosphorylation and secretion in neurons, in addition to misregulating the retromer, an endosomal complex linked to AD pathogenesis. We further provide evidence for increased heterodimer signaling in AD brain and identify dovitinib as a candidate molecule for normalizing β-amyloid-mediated transcriptional responses. The findings overall reveal differential transcription factor dimerization as a mechanism linking disease stimuli to the development of pathogenic cellular states.

Disease-associated non-coding variants alter NKX2-5 DNA-binding affinity

Genome-wide association studies (GWAS) have mapped over 90 % of disease- or trait-associated variants within the non-coding genome, like cis-regulatory elements (CREs). Non-coding single nucleotide polymorphisms (SNPs) are genomic variants that can change how DNA-binding regulatory proteins, like transcription factors (TFs), interact with the genome and regulate gene expression. NKX2-5 is a TF essential for proper heart development, and mutations affecting its function have been associated with congenital heart diseases (CHDs). However, establishing a causal mechanism between non-coding genomic variants and human disease remains challenging. To address this challenge, we identified 8475 SNPs predicted to alter NKX2-5 DNA-binding using a position weight matrix (PWM)-based predictive model. Five variants were prioritized for in vitro validation; four of them are associated with traits and diseases that impact cardiovascular health. The impact of these variants on NKX2-5 binding was evaluated with electrophoretic mobility shift assay (EMSA) using purified recombinant NKX2-5 homeodomain. Binding curves were constructed to determine changes in binding between variant and reference alleles. Variants rs7350789, rs7719885, rs747334, and rs3892630 increased binding affinity, whereas rs61216514 decreased binding by NKX2-5 when compared to the reference genome. Our findings suggest that differential TF-DNA binding affinity can be key in establishing a causal mechanism of pathogenic variants.

Regulatory networks driving expression of genes critical for glioblastoma are controlled by the transcription factor c-Jun and the pre-existing epigenetic modifications

BACKGROUND

Glioblastoma (GBM, WHO grade IV) is an aggressive, primary brain tumor. Despite extensive tumor resection followed by radio- and chemotherapy, life expectancy of GBM patients did not improve over decades. Several studies reported transcription deregulation in GBMs, but regulatory mechanisms driving overexpression of GBM-specific genes remain largely unknown. Transcription in open chromatin regions is directed by transcription factors (TFs) that bind to specific motifs, recruit co-activators/repressors and the transcriptional machinery. Identification of GBM-related TFs-gene regulatory networks may reveal new and targetable mechanisms of gliomagenesis.

RESULTS

We predicted TFs-regulated networks in GBMs in silico and intersected them with putative TF binding sites identified in the accessible chromatin in human glioma cells and GBM patient samples. The Cancer Genome Atlas and Glioma Atlas datasets (DNA methylation, H3K27 acetylation, transcriptomic profiles) were explored to elucidate TFs-gene regulatory networks and effects of the epigenetic background. In contrast to the majority of tumors, c-Jun expression was higher in GBMs than in normal brain and c-Jun binding sites were found in multiple genes overexpressed in GBMs, including VIM, FOSL2 or UPP1. Binding of c-Jun to the VIM gene promoter was stronger in GBM-derived cells than in cells derived from benign glioma as evidenced by gel shift and supershift assays. Regulatory regions of the majority of c-Jun targets have distinct DNA methylation patterns in GBMs as compared to benign gliomas, suggesting the contribution of DNA methylation to the c-Jun-dependent gene expression.

CONCLUSIONS

GBM-specific TFs-gene networks identified in GBMs differ from regulatory pathways attributed to benign brain tumors and imply a decisive role of c-Jun in controlling genes that drive glioma growth and invasion as well as a modulatory role of DNA methylation.

A stay of execution: ATF4 regulation and potential outcomes for the integrated stress response

ATF4 is a cellular stress induced bZIP transcription factor that is a hallmark effector of the integrated stress response. The integrated stress response is triggered by phosphorylation of the alpha subunit of the eukaryotic initiation factor 2 complex that can be carried out by the cellular stress responsive kinases; GCN2, PERK, PKR, and HRI. eIF2α phosphorylation downregulates mRNA translation initiation en masse, however ATF4 translation is upregulated. The integrated stress response can output two contradicting outcomes in cells; pro-survival or apoptosis. The mechanism for choice between these outcomes is unknown, however combinations of ATF4 heterodimerisation partners and post-translational modifications have been linked to this regulation. This semi-systematic review article covers ATF4 target genes, heterodimerisation partners and post-translational modifications. Together, this review aims to be a useful resource to elucidate the mechanisms controlling the effects of the integrated stress response. Additional putative roles of the ATF4 protein in cell division and synaptic plasticity are outlined.

MiRNA Differences Related to Treatment-Resistant Schizophrenia

Schizophrenia (SZ) is a serious mental disorder that is typically treated with antipsychotic medication. Treatment-resistant schizophrenia (TRS) is the condition where symptoms remain after pharmacological intervention, resulting in long-lasting functional and social impairments. As the identification and treatment of a TRS patient requires previous failed treatments, early mechanisms of detection are needed in order to quicken the access to effective therapy, as well as improve treatment adherence. In this study, we aim to find a microRNA (miRNA) signature for TRS, as well as to shed some light on the molecular pathways potentially involved in this severe condition. To do this, we compared the blood miRNAs of schizophrenia patients that respond to medication and TRS patients, thus obtaining a 16-miRNA TRS profile. Then, we assessed the ability of this signature to separate responders and TRS patients using hierarchical clustering, observing that most of them are grouped correctly (~70% accuracy). We also conducted a network, pathway analysis, and bibliography search to spot molecular pathways potentially altered in TRS. We found that the response to stress seems to be a key factor in TRS and that proteins p53, SIRT1, MDM2, and TRIM28 could be the potential mediators of such responses. Finally, we suggest a molecular pathway potentially regulated by the miRNAs of the TRS profile.

C/EBP-Family Redundancy Determines Patient Survival and Lymph Node Involvement in PDAC

Pancreatic ductal adenocarcinoma (PDAC) is a dismal disease with a poor clinical prognosis and unsatisfactory treatment options. We previously found that the transcription factor CCAAT/Enhancer-Binding Protein Delta (C/EBPδ) is lowly expressed in PDAC compared to healthy pancreas duct cells, and that patient survival and lymph node involvement in PDAC is correlated with the expression of C/EBPδ in primary tumor cells. C/EBPδ shares a homologous DNA-binding sequence with other C/EBP-proteins, leading to the presumption that other C/EBP-family members might act redundantly and compensate for the loss of C/EBPδ. This implies that patient stratification could be improved when expression levels of multiple C/EBP-family members are considered simultaneously. In this study, we assessed whether the quantification of C/EBPβ or C/EBPγ in addition to that of C/EBPδ might improve the prediction of patient survival and lymph node involvement using a cohort of 68 resectable PDAC patients. Using Kaplan-Meier analyses of patient groups with different C/EBP-expression levels, we found that both C/EBPβ and C/EBPγ can partially compensate for low C/EBPδ and improve patient survival. Further, we uncovered C/EBPβ as a novel predictor of a decreased likelihood of lymph node involvement in PDAC, and found that C/EBPβ and C/EBPδ can compensate for the lack of each other in order to reduce the risk of lymph node involvement. C/EBPγ, on the other hand, appears to promote lymph node involvement in the absence of C/EBPδ. Altogether, our results show that the redundancy of C/EBP-family members might have a profound influence on clinical prognoses and that the expression of both C/EPBβ and C/EBPγ should be taken into account when dichotomizing patients according to C/EBPδ expression.

Microglial activating transcription factor 3 upregulation: An indirect target to attenuate inflammation in the nervous system

Activating Transcription Factor 3 (ATF3) is upregulated in reaction to several cellular stressors found in a wide range of pathological conditions to coordinate a transcriptional response. ATF3 was first implicated in the transcriptional reaction to axotomy when its massive upregulation was measured in sensory and motor neuron cell bodies following peripheral nerve injury. It has since been shown to be critical for successful axon regeneration in the peripheral nervous system and a promising target to mitigate regenerative failure in the central nervous system. However, much of the research to date has focused on ATF3's function in neurons, leaving the expression, function, and therapeutic potential of ATF3 in glia largely unexplored. In the immunology literature ATF3 is seen as a master regulator of the innate immune system. Specifically, in macrophages following pathogen or damage associated molecular pattern receptor activation and subsequent cytokine release, ATF3 upregulation abrogates the inflammatory response. Importantly, ATF3 upregulation is not exclusively due to cellular stress exposure but has been achieved by the administration of several small molecules. In the central nervous system, microglia represent the resident macrophage population and are therefore of immediate interest with respect to ATF3 induction. It is our perspective that the potential of inducing ATF3 expression to dampen inflammatory microglial phenotype represents an unexplored therapeutic target and may have synergistic benefits when paired with concomitant neuronal ATF3 upregulation. This would be of particular benefit in pathologies that involve both detrimental inflammation and neuronal damage including spinal cord injury, multiple sclerosis, and stroke.

Structural basis of transcription regulation by CNC family transcription factor, Nrf2

Several basic leucine zipper (bZIP) transcription factors have accessory motifs in their DNA-binding domains, such as the CNC motif of CNC family or the EHR motif of small Maf (sMaf) proteins. CNC family proteins heterodimerize with sMaf proteins to recognize CNC-sMaf binding DNA elements (CsMBEs) in competition with sMaf homodimers, but the functional role of the CNC motif remains elusive. In this study, we report the crystal structures of Nrf2/NFE2L2, a CNC family protein regulating anti-stress transcriptional responses, in a complex with MafG and CsMBE. The CNC motif restricts the conformations of crucial Arg residues in the basic region, which form extensive contact with the DNA backbone phosphates. Accordingly, the Nrf2-MafG heterodimer has approximately a 200-fold stronger affinity for CsMBE than canonical bZIP proteins, such as AP-1 proteins. The high DNA affinity of the CNC-sMaf heterodimer may allow it to compete with the sMaf homodimer on target genes without being perturbed by other low-affinity bZIP proteins with similar sequence specificity.

The F-bZIP-regulated Zn deficiency response in land plants

This review describes zinc sensing and transcriptional regulation of the zinc deficiency response in Arabidopsis, and discusses how their evolutionary conservation in land plants facilitates translational approaches for improving the Zn nutritional value of crop species. Zinc is an essential micronutrient for all living organisms due to its presence in a large number of proteins, as a structural or catalytic cofactor. In plants, zinc homeostasis mechanisms comprise uptake from soil, transport and distribution throughout the plant to provide adequate cellular zinc availability. Here, I discuss the transcriptional regulation of the response to zinc deficiency and the zinc sensing mechanisms in Arabidopsis, and their evolutionary conservation in land plants. The Arabidopsis F-group basic region leucine-zipper (F-bZIP) transcription factors bZIP19 and bZIP23 function simultaneously as sensors of intracellular zinc status, by direct binding of zinc ions, and as the central regulators of the zinc deficiency response, with their target genes including zinc transporters from the ZRT/IRT-like Protein (ZIP) family and nicotianamine synthase enzymes that produce the zinc ligand nicotianamine. I note that this relatively simple mechanism of zinc sensing and regulation, together with the evolutionary conservation of F-bZIP transcription factors across land plants, offer important research opportunities. One of them is to use the F-bZIP-regulated zinc deficiency response as a tractable module for evolutionary and comparative functional studies. Another research opportunity is translational research in crop plants, modulating F-bZIP activity as a molecular switch to enhance zinc accumulation. This should become a useful plant-based solution to alleviate effects of zinc deficiency in soils, which impact crop production and crop zinc content, with consequences for human nutrition globally.

TaFDL2-1A confers drought stress tolerance by promoting ABA biosynthesis, ABA responses, and ROS scavenging in transgenic wheat

Plants have developed various protective mechanisms to survive drought stress. Previously, it was shown that a wheat bZIP transcription factor gene TaFD-Like2-1A (TaFDL2-1A) can confer drought tolerance in Arabidopsis. However, the biological functions related to drought stress tolerance of TaFDL2-1A in wheat (Triticum aestivum L.) remain unclear. In the present study, overexpression of TaFDL2-1A in the wheat cultivar Fielder improved drought resistance and conferred abscisic acid (ABA) hypersensitivity. Further analysis showed that overexpression of TaFDL2-1A increased the hypersensitivity of stomata to drought stress and endogenous ABA content under drought conditions. Genetic analysis and transcriptional regulation analysis indicated that TaFDL2-1A binds directly to the promoter fragments of TaRAB21s and TaNCED2s via ACGT core cis-elements, thereby activating their expression, leading to enhanced ABA responses and endogenous ABA accumulation. In addition, our results demonstrate that overexpression of TaFDL2-1A results in higher SOD and GPX activities in wheat under drought conditions by promoting the expression of TaSOD1 and TaGPx1-D, indicating enhanced reactive oxygen species (ROS) scavenging. These results imply that TaFDL2-1A positively regulates ABA biosynthesis, ABA responses, and ROS scavenging to improve drought stress tolerance in transgenic wheat. Our findings improve our understanding of the mechanisms that allow the wheat bZIP transcription factor to improve drought resistance and provide a useful reference gene for breeding programs to enhance drought resistance.

Relative importance of composition structures and biologically meaningful logics in bipartite Boolean models of gene regulation

Boolean networks have been widely used to model gene networks. However, such models are coarse-grained to an extent that they abstract away molecular specificities of gene regulation. Alternatively, bipartite Boolean network models of gene regulation explicitly distinguish genes from transcription factors (TFs). In such bipartite models, multiple TFs may simultaneously contribute to gene regulation by forming heteromeric complexes, thus giving rise to composition structures. Since bipartite Boolean models are relatively recent, an empirical investigation of their biological plausibility is lacking. Here, we estimate the prevalence of composition structures arising through heteromeric complexes. Moreover, we present an additional mechanism where composition structures may arise as a result of multiple TFs binding to cis-regulatory regions and provide empirical support for this mechanism. Next, we compare the restriction in BFs imposed by composition structures and by biologically meaningful properties. We find that though composition structures can severely restrict the number of Boolean functions (BFs) driving a gene, the two types of minimally complex BFs, namely nested canalyzing functions (NCFs) and read-once functions (RoFs), are comparatively more restrictive. Finally, we find that composition structures are highly enriched in real networks, but this enrichment most likely comes from NCFs and RoFs.

Control of ribosomal RNA synthesis by hematopoietic transcription factors

Ribosomal RNAs (rRNAs) are the most abundant cellular RNAs, and their synthesis from rDNA repeats by RNA polymerase I accounts for the bulk of all transcription. Despite substantial variation in rRNA transcription rates across cell types, little is known about cell-type-specific factors that bind rDNA and regulate rRNA transcription to meet tissue-specific needs. Using hematopoiesis as a model system, we mapped about 2,200 ChIP-seq datasets for 250 transcription factors (TFs) and chromatin proteins to human and mouse rDNA and identified robust binding of multiple TF families to canonical TF motifs on rDNA. Using a 47S-FISH-Flow assay developed for nascent rRNA quantification, we demonstrated that targeted degradation of C/EBP alpha (CEBPA), a critical hematopoietic TF with conserved rDNA binding, caused rapid reduction in rRNA transcription due to reduced RNA Pol I occupancy. Our work identifies numerous potential rRNA regulators and provides a template for dissection of TF roles in rRNA transcription.

Flexibility of flanking DNA is a key determinant of transcription factor affinity for the core motif

Selective gene regulation is mediated by recognition of specific DNA sequences by transcription factors (TFs). The extremely challenging task of searching out specific cognate DNA binding sites among several million putative sites within the eukaryotic genome is achieved by complex molecular recognition mechanisms. Elements of this recognition code include the core binding sequence, the flanking sequence context, and the shape and conformational flexibility of the composite binding site. To unravel the extent to which DNA flexibility modulates TF binding, in this study, we employed experimentally guided molecular dynamics simulations of ternary complex of closely related Hox heterodimers Exd-Ubx and Exd-Scr with DNA. Results demonstrate that flexibility signatures embedded in the flanking sequences impact TF binding at the cognate binding site. A DNA sequence has intrinsic shape and flexibility features. While shape features are localized, our analyses reveal that flexibility features of the flanking sequences percolate several basepairs and allosterically modulate TF binding at the core. We also show that lack of flexibility in the motif context can render the cognate site resistant to protein-induced shape changes and subsequently lower TF binding affinity. Overall, this study suggests that flexibility-guided DNA shape, and not merely the static shape, is a key unexplored component of the complex DNA-TF recognition code.

Prediction of protein-ligand binding affinity from sequencing data with interpretable machine learning

Protein-ligand interactions are increasingly profiled at high throughput using affinity selection and massively parallel sequencing. However, these assays do not provide the biophysical parameters that most rigorously quantify molecular interactions. Here we describe a flexible machine learning method, called ProBound, that accurately defines sequence recognition in terms of equilibrium binding constants or kinetic rates. This is achieved using a multi-layered maximum-likelihood framework that models both the molecular interactions and the data generation process. We show that ProBound quantifies transcription factor (TF) behavior with models that predict binding affinity over a range exceeding that of previous resources; captures the impact of DNA modifications and conformational flexibility of multi-TF complexes; and infers specificity directly from in vivo data such as ChIP-seq without peak calling. When coupled with an assay called KD-seq, it determines the absolute affinity of protein-ligand interactions. We also apply ProBound to profile the kinetics of kinase-substrate interactions. ProBound opens new avenues for decoding biological networks and rationally engineering protein-ligand interactions.

BindVAE: Dirichlet variational autoencoders for de novo motif discovery from accessible chromatin

We present a novel unsupervised deep learning approach called BindVAE, based on Dirichlet variational autoencoders, for jointly decoding multiple TF binding signals from open chromatin regions. BindVAE can disentangle an input DNA sequence into distinct latent factors that encode cell-type specific in vivo binding signals for individual TFs, composite patterns for TFs involved in cooperative binding, and genomic context surrounding the binding sites. On the task of retrieving the motifs of expressed TFs in a given cell type, BindVAE is competitive with existing motif discovery approaches.

AP-1 transcription factor network explains diverse patterns of cellular plasticity in melanoma cells

Cellular plasticity associated with fluctuations in transcriptional programs allows individual cells in a tumor to adopt heterogeneous differentiation states and switch phenotype during their adaptive responses to therapies. Despite increasing knowledge of such transcriptional programs, the molecular basis of cellular plasticity remains poorly understood. Here, we combine multiplexed transcriptional and protein measurements at population and single-cell levels with multivariate statistical modeling to show that the state of AP-1 transcription factor network plays a unifying role in explaining diverse patterns of plasticity in melanoma. We find that a regulated balance among AP-1 factors cJUN, JUND, FRA2, FRA1, and cFOS determines the intrinsic diversity of differentiation states and adaptive responses to MAPK inhibitors in melanoma cells. Perturbing this balance through genetic depletion of specific AP-1 proteins, or by MAPK inhibitors, shifts cellular heterogeneity in a predictable fashion. Thus, AP-1 may serve as a critical node for manipulating cellular plasticity with potential therapeutic implications.

The protein phosphatase 2C clade A TaPP2CA interact with calcium-dependent protein kinases, TaCDPK5/TaCDPK9-1, that phosphorylate TabZIP60 transcription factor from wheat (Triticum aestivum L.)

Previously we have found that TabZIP60 from the ABF/AREB (ABRE-binding factor/ABA-responsive element-binding protein) subfamily of bZIP transcription factor (TF) was involved in salt stress response. However, the regulatory mechanism of TabZIP60 is unknown. In the present study, we identified two calcium-dependent protein kinase (CDPK) genes, TaCDPK5/TaCDPK9-1, which were clustered into group Ⅰ and were induced by salt, abscisic acid (ABA), and polyethylene glycol (PEG) treatments. RT-qPCR results showed that the expression level of salt-induced TabZIP60 was drastically inhibited by Ca²⁺ channel blocker LaCl₃. TaCDPK5/TaCDPK9-1 were involved in interaction with TabZIP60 protein in vivo and in vitro. And TaCDPK5/TaCDPK9-1 could autophosphorylate and phosphorylate TabZIP60 protein in a Ca²⁺-dependent way. Mutational analysis indicated that Serine-110 of TabZIP60 was essential for TaCDPK5/TaCDPK9-1-TabZIP60 interaction and was the phosphorylation site of TaCDPK5/TaCDPK9-1 kinases. Yeast two-hybrid assay results showed the interactions between TaCDPK5/TaCDPK9-1 and wheat protein phosphatase 2 C clade A TaPP2CA116/ TaPP2CA121 separately. These findings demonstrate that the phosphorylation status of TabZIP60 controlled by TaPP2CA116/ TaPP2CA121 and TaCDPK5/TaCDPK9-1 might play a crucial role in wheat during salt stress.

TGA transcription factors-Structural characteristics as basis for functional variability

TGA transcription factors are essential regulators of various cellular processes, their activity connected to different hormonal pathways, interacting proteins and regulatory elements. Belonging to the basic region leucine zipper (bZIP) family, TGAs operate by binding to their target DNA sequence as dimers through a conserved bZIP domain. Despite sharing the core DNA-binding sequence, the TGA paralogues exert somewhat different DNA-binding preferences. Sequence variability of their N- and C-terminal protein parts indicates their importance in defining TGA functional specificity through interactions with diverse proteins, affecting their DNA-binding properties. In this review, we provide a short and concise summary on plant TGA transcription factors from a structural point of view, including the relation of their structural characteristics to their functional roles in transcription regulation.

The Time Sequence of Gene Expression Changes after Spinal Cord Injury

Gene expression changes following spinal cord injury (SCI) are time-dependent, and an accurate understanding of these changes can be crucial in determining time-based treatment options in a clinical setting. We performed RNA sequencing of the contused spinal cord of rats at five different time points from the very acute to chronic stages (1 hour, 1 day, 1 week, 1 month, and 3 months) following SCI. We identified differentially expressed genes (DEGs) and Gene Ontology (GO) terms at each time point, and 14,257 genes were commonly expressed at all time points. The biological process of the inflammatory response was increased at 1 hour and 1 day, and the cellular component of the integral component of the synaptic membrane was increased at 1 day. DEGs associated with cell activation and the innate immune response were highly enriched at 1 week and 1 month, respectively. A total of 2841 DEGs were differentially expressed at any of the five time points, and 18 genes (17 upregulated and 1 downregulated) showed common expression differences at all time points. We found that interleukin signaling, neutrophil degranulation, eukaryotic translation, collagen degradation, LGI–ADAM interactions, GABA receptor, and L1CAM-ankyrin interactions were prominent after SCI depending on the time post injury. We also performed gene–drug network analysis and found several potential antagonists and agonists which can be used to treat SCI. We expect to discover effective treatments in the clinical field through further studies revealing the efficacy and safety of potential drugs.

An Approach to Derive Functional Peptide Inhibitors of Transcription Factor Activity

We report the development of a high-throughput, intracellular "transcription block survival" (TBS) screening platform to derive functional transcription factor antagonists. TBS is demonstrated using the oncogenic transcriptional regulator cJun, with the development of antagonists that bind cJun and prevent both dimerization and, more importantly, DNA binding remaining a primary challenge. In TBS, cognate TRE sites are introduced into the coding region of the essential gene, dihydrofolate reductase (DHFR). Introduction of cJun leads to TRE binding, preventing DHFR expression by directly blocking RNA polymerase gene transcription to abrogate cell proliferation. Peptide library screening identified a sequence that both binds cJun and antagonizes function by preventing DNA binding, as demonstrated by restored cell viability and subsequent in vitro hit validation. TBS is an entirely tag-free genotype-to-phenotype approach, selecting desirable attributes such as high solubility, target specificity, and low toxicity within a complex cellular environment. TBS facilitates rapid library screening to accelerate the identification of therapeutically valuable sequences.

Activating Transcription Factor 3 (ATF3) is a Highly Conserved Pro-regenerative Transcription Factor in the Vertebrate Nervous System

The vertebrate nervous system exhibits dramatic variability in regenerative capacity across species and neuronal populations. For example, while the mammalian central nervous system (CNS) is limited in its regenerative capacity, the CNS of many other vertebrates readily regenerates after injury, as does the peripheral nervous system (PNS) of mammals. Comparing molecular responses across species and tissues can therefore provide valuable insights into both conserved and distinct mechanisms of successful regeneration. One gene that is emerging as a conserved pro-regenerative factor across vertebrates is activating transcription factor 3 (ATF3), which has long been associated with tissue trauma. A growing number of studies indicate that ATF3 may actively promote neuronal axon regrowth and regeneration in species ranging from lampreys to mammals. Here, we review data on the structural and functional conservation of ATF3 protein across species. Comparing RNA expression data across species that exhibit different abilities to regenerate their nervous system following traumatic nerve injury reveals that ATF3 is consistently induced in neurons within the first few days after injury. Genetic deletion or knockdown of ATF3 expression has been shown in mouse and zebrafish, respectively, to reduce axon regeneration, while inducing ATF3 promotes axon sprouting, regrowth, or regeneration. Thus, we propose that ATF3 may be an evolutionarily conserved regulator of neuronal regeneration. Identifying downstream effectors of ATF3 will be a critical next step in understanding the molecular basis of vertebrate CNS regeneration.

Induction of heparanase 2 (Hpa2) expression by stress is mediated by ATF3

Activity of heparanase, endoglycosidase that cleaves heparan sulfate side chains in heparan sulfate proteoglycans, is highly implicated in tumor progression and metastasis. Heparanase inhibitors are therefore being evaluated clinically as anti-cancer therapeutics. Heparanase 2 (Hpa2) is a close homolog of heparanase that lacks HS-degrading activity and functions as an endogenous inhibitor of heparanase. As a result, Hpa2 appears to attenuate tumor growth but mechanisms that regulate Hpa2 expression and determine the ratio between heparanase and Hpa2 are largely unknown. We have recently reported that the expression of Hpa2 is induced by endoplasmic reticulum (ER) and proteotoxic stresses, but the mechanism(s) underlying Hpa2 gene regulation was obscure. Here we expand the notion that Hpa2 is regulated by conditions of stress. We report that while ER and hypoxia, each alone, resulted in a 3-7 fold increase in Hpa2 expression, combining ER stress and hypoxia resulted in a noticeable, over 40-fold increase in Hpa2 expression. A prominent induction of Hpa2 expression was also quantified in cells exposed to heat shock, proteotoxic stress, lysosomal stress, and chemotherapy (cisplatin), strongly implying that Hpa2 is regulated by conditions of stress. Furthermore, analyses of the Hpa2 gene promoter led to the identification of activating-transcription-factor 3 (ATF3) as a transcription factor that mediates Hpa2 induction by stress, thus revealing, for the first time, a molecular mechanism that underlies Hpa2 gene regulation. Induction of Hpa2 and ATF3 by conditions of stress that often accompany the rapid expansion of tumors is likely translated to improved survival of cancer patients.

TRIB1 regulates LDL metabolism through CEBPα-mediated effects on the LDL receptor in hepatocytes

Genetic variants near the TRIB1 gene are highly significantly associated with plasma lipid traits and coronary artery disease. While TRIB1 is likely causal of these associations, the molecular mechanisms are not well understood. Here we sought to investigate how TRIB1 influences low density lipoprotein cholesterol (LDL-C) levels in mice. Hepatocyte-specific deletion of Trib1 (Trib1Δhep) in mice increased plasma cholesterol and apoB and slowed the catabolism of LDL-apoB due to decreased levels of LDL receptor (LDLR) mRNA and protein. Simultaneous deletion of the transcription factor CCAAT/enhancer-binding protein alpha (CEBPα) with TRIB1 eliminated the effects of TRIB1 on hepatic LDLR regulation and LDL catabolism. Using RNA-seq, we found that activating transcription factor 3 (Atf3) was highly upregulated in the livers of Trib1Δhep but not Trib1Δhep CebpaΔhep mice. ATF3 has been shown to directly bind to the CEBPα protein, and to repress the expression of LDLR by binding its promoter. Blunting the increase of ATF3 in Trib1Δhep mice reduced the levels of plasma cholesterol and partially attenuated the effects on LDLR. Based on these data, we conclude that deletion of Trib1 leads to a posttranslational increase in CEBPα, which increases ATF3 levels, thereby contributing to the downregulation of LDLR and increased plasma LDL-C.

The bZIP Transcription Factor ZIP-11 Is Required for the Innate Immune Regulation in Caenorhabditis elegans

Innate immunity is the first line of host defense against pathogen infection in metazoans. However, the molecular mechanisms of the complex immune regulatory network are not fully understood. Based on a transcriptome profiling of the nematode Caenorhabditis elegans, we found that a bZIP transcription factor ZIP-11 was up-regulated upon Pseudomonas aeruginosa PA14 infection. The tissue specific RNAi knock-down and rescue data revealed that ZIP-11 acts in intestine to promote host resistance against P. aeruginosa PA14 infection. We further showed that intestinal ZIP-11 regulates innate immune response through constituting a feedback loop with the conserved PMK-1/p38 mitogen-activated protein signaling pathway. Intriguingly, ZIP-11 interacts with a CCAAT/enhancer-binding protein, CEBP-2, to mediate the transcriptional response to P. aeruginosa PA14 infection independently of PMK-1/p38 pathway. In addition, human homolog ATF4 can functionally substitute for ZIP-11 in innate immune regulation of C. elegans. Our findings indicate that the ZIP-11/ATF4 genetic program activates local innate immune response through conserved PMK-1/p38 and CEBP-2/C/EBPγ immune signals in C. elegans, raising the possibility that a similar process may occur in other organisms.

The Jun-dependent axon regeneration gene program: Jun promotes regeneration over plasticity

The regeneration-associated gene (RAG) expression program is activated in injured peripheral neurons after axotomy and enables long-distance axon re-growth. Over 1000 genes are regulated, and many transcription factors are upregulated or activated as part of this response. However, a detailed picture of how RAG expression is regulated is lacking. In particular, the transcriptional targets and specific functions of the various transcription factors are unclear. Jun was the first-regeneration-associated transcription factor identified and the first shown to be functionally important. Here we fully define the role of Jun in the RAG expression program in regenerating facial motor neurons. At 1, 4 and 14 days after axotomy, Jun upregulates 11, 23 and 44% of the RAG program, respectively. Jun functions relevant to regeneration include cytoskeleton production, metabolic functions and cell activation, and the downregulation of neurotransmission machinery. In silico analysis of promoter regions of Jun targets identifies stronger over-representation of AP1-like sites than CRE-like sites, although CRE sites were also over-represented in regions flanking AP1 sites. Strikingly, in motor neurons lacking Jun, an alternative SRF-dependent gene expression program is initiated after axotomy. The promoters of these newly expressed genes exhibit over-representation of CRE sites in regions near to SRF target sites. This alternative gene expression program includes plasticity-associated transcription factors and leads to an aberrant early increase in synapse density on motor neurons. Jun thus has the important function in the early phase after axotomy of pushing the injured neuron away from a plasticity response and towards a regenerative phenotype.

ATF3 promotes the serine synthesis pathway and tumor growth under dietary serine restriction

The serine synthesis pathway (SSP) involving metabolic enzymes phosphoglycerate dehydrogenase (PHGDH), phosphoserine aminotransferase 1 (PSAT1), and phosphoserine phosphatase (PSPH) drives intracellular serine biosynthesis and is indispensable for cancer cells to grow in serine-limiting environments. However, how SSP is regulated is not well understood. Here, we report that activating transcription factor 3 (ATF3) is crucial for transcriptional activation of SSP upon serine deprivation. ATF3 is rapidly induced by serine deprivation via a mechanism dependent on ATF4, which in turn binds to ATF4 and increases the stability of this master regulator of SSP. ATF3 also binds to the enhancers/promoters of PHGDH, PSAT1, and PSPH and recruits p300 to promote expression of these SSP genes. As a result, loss of ATF3 expression impairs serine biosynthesis and the growth of cancer cells in the serine-deprived medium or in mice fed with a serine/glycine-free diet. Interestingly, ATF3 expression positively correlates with PHGDH expression in a subset of TCGA cancer samples.

Activation of the serine synthesis pathway is important for cancer cell growth, but how this pathway is regulated is not well understood. Li et al. report that ATF3 is an important regulator of this pathway and can promote serine biosynthesis and tumor growth under serine-limiting conditions.

Sequence-specific dynamics of DNA response elements and their flanking sites regulate the recognition by AP-1 transcription factors

Activator proteins 1 (AP-1) comprise one of the largest families of eukaryotic basic leucine zipper transcription factors. Despite advances in the characterization of AP-1 DNA-binding sites, our ability to predict new binding sites and explain how the proteins achieve different gene expression levels remains limited. Here we address the role of sequence-specific DNA flexibility for stability and specific binding of AP-1 factors, using microsecond-long molecular dynamics simulations. As a model system, we employ yeast AP-1 factor Yap1 binding to three different response elements from two genetic environments. Our data show that Yap1 actively exploits the sequence-specific flexibility of DNA within the response element to form stable protein–DNA complexes. The stability also depends on the four to six flanking nucleotides, adjacent to the response elements. The flanking sequences modulate the conformational adaptability of the response element, making it more shape-efficient to form specific contacts with the protein. Bioinformatics analysis of differential expression of the studied genes supports our conclusions: the stability of Yap1–DNA complexes, modulated by the flanking environment, influences the gene expression levels. Our results provide new insights into mechanisms of protein–DNA recognition and the biological regulation of gene expression levels in eukaryotes.