Immunology. Cooperative transcription factor complexes in control

Hold out the genome: a roadmap to solving the cis-regulatory code

Gene expression is regulated by transcription factors that work together to read cis-regulatory DNA sequences. The 'cis-regulatory code' - how cells interpret DNA sequences to determine when, where and how much genes should be expressed - has proven to be exceedingly complex. Recently, advances in the scale and resolution of functional genomics assays and machine learning have enabled substantial progress towards deciphering this code. However, the cis-regulatory code will probably never be solved if models are trained only on genomic sequences; regions of homology can easily lead to overestimation of predictive performance, and our genome is too short and has insufficient sequence diversity to learn all relevant parameters. Fortunately, randomly synthesized DNA sequences enable testing a far larger sequence space than exists in our genomes, and designed DNA sequences enable targeted queries to maximally improve the models. As the same biochemical principles are used to interpret DNA regardless of its source, models trained on these synthetic data can predict genomic activity, often better than genome-trained models. Here we provide an outlook on the field, and propose a roadmap towards solving the cis-regulatory code by a combination of machine learning and massively parallel assays using synthetic DNA.

TFSyntax: a database of transcription factors binding syntax in mammalian genomes

In mammals, transcriptional factors (TFs) drive gene expression by binding to regulatory elements in a cooperative manner. Deciphering the rules of such cooperation is crucial to obtain a full understanding of cellular homeostasis and development. Although this is a long-standing topic, there is no comprehensive database for biologists to access the syntax of TF binding sites. Here we present TFSyntax (https://tfsyntax.zhaopage.com), a database focusing on the arrangement of TF binding sites. TFSyntax maps the binding motif of 1299 human TFs and 890 mouse TFs across 382 cells and tissues, representing the most comprehensive TF binding map to date. In addition to location, TFSyntax defines motif positional preference, density and colocalization within accessible elements. Powered by a series of functional modules based on web interface, users can freely search, browse, analyze, and download data of interest. With comprehensive characterization of TF binding syntax across distinct tissues and cell types, TFSyntax represents a valuable resource and platform for studying the mechanism of transcriptional regulation and exploring how regulatory DNA variants cause disease.

CTCF DNA binding domain undergoes dynamic and selective protein–protein interactions

CTCF is a predominant insulator protein required for three-dimensional chromatin organization. However, the roles of its insulation of enhancers in a 3D nuclear organization have not been fully explained. Here, we found that the CTCF DNA-binding domain (DBD) forms dynamic self-interacting clusters. Strikingly, CTCF DBD clusters were found to incorporate other insulator proteins but are not coenriched with transcriptional activators in the nucleus. This property is not observed in other domains of CTCF or the DBDs of other transcription factors. Moreover, endogenous CTCF shows a phenotype consistent with the DBD by forming small protein clusters and interacting with CTCF motif arrays that have fewer transcriptional activators bound. Our results reveal an interesting phenomenon in which CTCF DBD interacts with insulator proteins and selectively localizes to nuclear positions with lower concentrations of transcriptional activators, providing insights into the insulation function of CTCF.

•

The CTCF DNA-binding domain forms protein clusters in vivo and in vitro

•

CTCF DBD clusters colocalize with insulator proteins but not with activators

•

Arginine residues of CTCF DBD are frequently mutated in cancers

•

Multiple transcription factor DBDs form protein clusters

c-Rel-dependent monocytes are potent immune suppressor cells in cancer

Myeloid-derived suppressor cells (MDSCs) are a heterogeneous population of leukocytes that are important for tumorigenesis and tumor immunotherapy. They comprise up to 10% of leukocytes in the blood of tumor patients and their depletion may be required for successful tumor immunotherapy. However, the identity of MDSCs remains obscure, primarily due to their heterogeneity and lack of a known lineage-specific transcription factor specifying their differentiation. Using single-cell transcriptomics and gene knockout approaches, we now describe a subset of murine and human myeloid suppressor cells, named rel-dependent monocytes (rMos), which are programmed by the transcription factor c-Rel of the NF-κB family. Unlike MDSCs described previously, the c-Rel-dependent monocytes expressed a high amount of the proinflammatory cytokine IL-1β together with a low level of suppressive molecule arginase 1. Both in vitro and in tumor-bearing mice, these c-Rel⁺ IL-1β^hi Arg1^- monocytes promoted tumor growth by potently suppressing T cell function and showed a strong migratory phenotype, all of which were impaired by c-Rel deficiency or inhibition. Mechanistic studies revealed that c-Rel controlled the expression of monocyte signature genes through a unique transcriptional complex called the c-Rel enhanceosome, and IL-1β-CCL2 crosstalk between tumor cells and the rel-dependent monocytes maintained the suppressive tumor microenvironment. Thus, c-Rel specifies the development of a suppressive monocyte population and could be selectively targeted for treating cancer.

Transcriptome-wide identification and expression profiling of the ERF gene family suggest roles as transcriptional activators and repressors of fruit ripening in durian

The involvement of the phytohormone ethylene as the main trigger of climacteric fruit ripening is well documented. However, our knowledge regarding the role of ethylene response factor (ERF) transcription factor in the transcriptional regulation of ethylene biosynthesis during fruit ripening remains limited. Here, comprehensive transcriptome analysis and expression profiling revealed 63 ERFs in durian pulps, termed DzERF1-DzERF63, of which 34 exhibited ripening-associated expression patterns at three stages (unripe, midripe, and ripe) during fruit ripening. Hierarchical clustering analysis classified 34 ripening-associated DzERFs into three distinct clades, among which, clade I consisted of downregulated DzERFs and clade III included those upregulated during ripening. Phylogenetic analysis predicted the functions of some DzERFs based on orthologs of previously characterized ERFs. Among downregulated DzERFs, DzERF6 functional prediction revealed its role as a negative regulator of ripening via ethylene biosynthetic gene repression, whereas among upregulated genes, DzERF9 was predicted to positively regulate ethylene biosynthesis. Correlation network analysis of 34 ripening-associated DzERFs with potential target genes revealed a strong negative correlation between DzERF6 and ethylene biosynthetic genes and a strong positive correlation between DzERF9 and ethylene biosynthesis. DzERF6 and DzERF9 showed differential expression patterns in association with different ripening treatments (natural, ethylene-induced, and 1-methylcyclopropene-delayed ripening). DzERF6 was downregulated, whereas DzERF9 was upregulated, during ripening and after ethylene treatment. The auxin-repressed and auxin-induced expression of DzERF6 and DzERF9, respectively, confirmed its dose-dependent responsiveness to exogenous auxin. We suggest ethylene- and auxin-mediated roles of DzERF6 and DzERF9 during fruit ripening, possibly through transcriptional regulation of ethylene biosynthetic genes.

A fruit ripening-associated transcription factor CsMADS5 positively regulates carotenoid biosynthesis in citrus

Carotenoids in citrus contribute to the quality of the fruit, but the mechanism of its transcriptional regulation is fairly unknown. Here, we characterized a citrus FRUITFULL sub-clade MADS gene, CsMADS5, that was ripening-inducible and acted as a nucleus-localized trans-activator. Transient overexpression of CsMADS5 in citrus induced fruit coloration and enhanced carotenoid concentrations. The expression of carotenogenic genes including phytoene synthase (PSY), phytoene desaturase (PDS), and lycopene β-cyclase 1 (LCYb1) was increased in the peels of fruits overexpressing CsMADS5. Similar results were observed from stable overexpression of CsMADS5 in tomato fruits and citrus calli, even though the effect of CsMADS5 on carotenoid metabolism in transgenic citrus calli was limited. Further biochemical analyses demonstrated that CsMADS5 activated the transcription of PSY, PDS, and LCYb1 by directly binding to their promoters. We concluded that CsMADS5 positively regulates carotenoid biosynthesis in fruits by directly activating the transcription of carotenogenic genes. Moreover, CsMADS5 physically interacted with a positive regulator CsMADS6, indicating that CsMADS5 may form an enhancer complex with CsMADS6 to synergistically promote carotenoid accumulation. These findings expand our understanding of the complex transcriptional regulatory hierarchy of carotenoid biosynthesis during fruit ripening.

The Multifaceted Output of c-Jun Biological Activity: Focus at the Junction of CD8 T Cell Activation and Exhaustion

c-Jun is a major component of the dimeric transcription factor activator protein-1 (AP-1), a paradigm for transcriptional response to extracellular signaling, whose components are basic-Leucine Zipper (bZIP) transcription factors of the Jun, Fos, activating transcription factor (ATF), ATF-like (BATF) and Jun dimerization protein 2 (JDP2) gene families. Extracellular signals regulate c-Jun/AP-1 activity at multiple levels, including transcriptional and posttranscriptional regulation of c-Jun expression and transactivity, in turn, establishing the magnitude and the duration of c-Jun/AP-1 activation. Another important level of c-Jun/AP-1 regulation is due to the capability of Jun family members to bind DNA as a heterodimer with every other member of the AP-1 family, and to interact with other classes of transcription factors, thereby acquiring the potential to integrate diverse extrinsic and intrinsic signals into combinatorial regulation of gene expression. Here, we review how these features of c-Jun/AP-1 regulation underlie the multifaceted output of c-Jun biological activity, eliciting quite distinct cellular responses, such as neoplastic transformation, differentiation and apoptosis, in different cell types. In particular, we focus on the current understanding of the role of c-Jun/AP-1 in the response of CD8 T cells to acute infection and cancer. We highlight the transcriptional and epigenetic regulatory mechanisms through which c-Jun/AP-1 participates in the productive immune response of CD8 T cells, and how its downregulation may contribute to the dysfunctional state of tumor infiltrating CD8 T cells. Additionally, we discuss recent insights pointing at c-Jun as a suitable target for immunotherapy-based combination approaches to reinvigorate anti-tumor immune functions.

BPC transcription factors and a Polycomb Group protein confine the expression of the ovule identity gene SEEDSTICK in Arabidopsis

The BASIC PENTACYSTEINE (BPC) GAGA (C-box) binding proteins belong to a small plant transcription factor family. We previously reported that class I BPCs bind directly to C-boxes in the SEEDSTICK (STK) promoter and the mutagenesis of these cis-elements affects STK expression in the flower. The MADS-domain factor SHORT VEGETATIVE PHASE (SVP) is another key regulator of STK. Direct binding of SVP to CArG-boxes in the STK promoter are required to repress its expression during the first stages of flower development. Here we show that class II BPCs directly interact with SVP and that MADS-domain binding sites in the STK promoter region are important for the correct spatial and temporal expression of this homeotic gene. Furthermore, we show that class I and class II BPCs act redundantly to repress STK expression in the flower, most likely by recruiting TERMINAL FLOWER 2/LIKE HETEROCHROMATIN PROTEIN 1 (TFL2/LHP1) and mediating the establishment and the maintenance of H3K27me3 repressive marks on DNA. We investigate the role of LHP1 in the regulation of STK expression. In addition to providing a better understanding of the role of BPC transcription factors in the regulation of STK expression, our results suggest the existence of a more general regulatory complex composed of BPCs, MADS-domain factors and Polycomb Repressive Complexes that co-operate to regulate gene expression in reproductive tissues. We believe that our data along with the molecular model described here could provide significant insights for a more comprehensive understanding of gene regulation in plants.

Auxin Response Factor 2A Is Part of the Regulatory Network Mediating Fruit Ripening Through Auxin-Ethylene Crosstalk in Durian

Fruit ripening is a highly coordinated developmental process driven by a complex hormonal network. Ethylene is the main regulator of climacteric fruit ripening. However, a putative role of other key phytohormones in this process cannot be excluded. We previously observed an increasing level of auxin during the post-harvest ripening of the durian fruit, which occurred concomitantly with the rise in the climacteric ethylene biosynthesis. Herein, we connect the key auxin signaling component, auxin response factors (ARFs), with the regulatory network that controls fruit ripening in durian through the identification and functional characterization of a candidate ripening-associated ARF. Our transcriptome-wide analysis identified 15 ARF members in durian (DzARFs), out of which 12 were expressed in the fruit pulp. Most of these DzARFs showed a differential expression, but DzARF2A had a marked ripening-associated expression pattern during post-harvest ripening in Monthong, a commercial durian cultivar from Thailand. Phylogenetic analysis of DzARF2A based on its tomato orthologue predicted a role in ripening through the regulation of ethylene biosynthesis. Transient expression of DzARF2A in Nicotiana benthamiana leaves significantly upregulated the expression levels of ethylene biosynthetic genes, pointing to a ripening-associated role of DzARF2A through the transcriptional regulation of ethylene biosynthesis. Dual-luciferase reporter assay determined that DzARF2A trans-activates durian ethylene biosynthetic genes. We previously reported significantly higher auxin level during post-harvest ripening in a fast-ripening cultivar (Chanee) compared to a slow-ripening one (Monthong). DzARF2A expression was significantly higher during post-harvest ripening in the fast-ripening cultivars (Chanee and Phuangmanee) compared to that of the slow-ripening ones (Monthong and Kanyao). Thus, higher auxin level could upregulate the expression of DzARF2A during ripening of a fast-ripening cultivar. The auxin-induced expression of DzARF2A confirmed its responsiveness to exogenous auxin treatment in a dose-dependent manner, suggesting an auxin-mediated role of DzARF2A in fruit ripening. We suggest that high DzARF2A expression would activate ARF2A-mediated transcription of ethylene biosynthetic genes, leading to increased climacteric ethylene biosynthesis (auxin-ethylene crosstalk) and faster ripening. Hence, we demonstrated DzARF2A as a new component of the regulatory network possibly mediating durian fruit ripening through transcriptional regulation of ethylene biosynthetic genes.

GABPα and CREB1 Binding to Double Nucleotide Polymorphisms of Their Consensus Motifs and Cooperative Binding to the Composite ETS ⇔ CRE Motif ()

Previously, cooperative binding of the bZIP domain of CREB1 and the ETS domain of GABPα was observed for the composite DNA ETS ⇔ CRE motif (A ⁰ C ¹ C ² G ³ G ⁴ A ⁵ A ⁶ G ⁷ T ⁸ G ⁹ A ¹⁰ C ¹¹ G ¹² T ¹³ C ¹⁴ A ¹⁵ ). Single nucleotide polymorphisms (SNPs) at the beginning and end of the ETS motif (ACCGGAAGT) increased cooperative binding. Here, we use an Agilent microarray of 60-mers containing all double nucleotide polymorphisms (DNPs) of the ETS ⇔ CRE motif to explore GABPα and CREB1 binding to their individual motifs and their cooperative binding. For GABPα, all DNPs were bound as if each SNP acted independently. In contrast, CREB1 binding to some DNPs was stronger or weaker than expected, depending on the locations of each SNP. CREB1 binding to DNPs where both SNPs were in the same half site, T ⁸ G ⁹ A ¹⁰ or T ¹³ C ¹⁴ A ¹⁵ , was greater than expected, indicating that an additional SNP cannot destroy binding as much as expected, suggesting that an individual SNP is enough to abolish sequence-specific DNA binding of a single bZIP monomer. If a DNP contains SNPs in each half site, binding is weaker than expected. Similar results were observed for additional ETS and bZIP family members. Cooperative binding between GABPα and CREB1 to the ETS ⇔ CRE motif was weaker than expected except for DNPs containing A ⁷ and SNPs at the beginning of the ETS motif.

Single cell RNA-Seq reveals pre-cDCs fate determined by transcription factor combinatorial dose

Background

Classic dendritic cells (cDCs) play a central role in the immune system by processing and presenting antigens to activate T cells, and consist of two major subsets: CD141⁺ cDC (cDC1) and CD1c⁺ cDC (cDC2). A population of migratory precursor cells, the pre-cDCs, is the immediate precursors to both cDC subsets. Previous studies showed that there were two pre-committed pre-cDC subpopulations. However, the key molecular drivers of pre-commitment in human pre-cDCs were not investigated.

Results

To identify the key molecular drivers for pre-commitment in human pre-cDCs, we performed single cell RNA sequencing (RNA-Seq) of two cDC subsets and pre-cDCs, and bulk RNA-Seq of pre-cDCs and cDCs from human peripheral blood. We found that pre-DC subpopulations cannot be separated by either variable genes within pre-cDCs or differentially expressed genes between cDC1 and cDC2. In contrast, they were separated by 16 transcription factors that are themselves differentially expressed or have regulated targets enriched in the differentially expressed genes between bulk cDC1 and cDC2, with one subpopulation close to cDC1 and the other close to cDC2. More importantly, these two pre-cDC sub-populations are correlated with ratio of IRF8 to IRF4 expression level more than their individual expression level. We also verified these findings using three recently published datasets.

Conclusions

In this study, we demonstrate that single cell transcriptome profiling can reveal pre-cDCs differentiation map, and our results suggest the concept that combinatorial dose of transcription factors determines cell differentiation fate.

Electronic supplementary material

The online version of this article (10.1186/s12860-019-0199-y) contains supplementary material, which is available to authorized users.

Chromatin Accessibility and Interactions in the Transcriptional Regulation of T Cells

During T cell differentiation and activation, specific stimuli, and a network of transcription factors (TFs) are involved in orchestrating chromatin accessibility, establishing enhancer-promoter interactions, and regulating gene expression. Over the past few years, there have been new insights into how chromatin interactions coordinate differentiation during T cell development and how regulatory elements are programmed to allow T cells to differentially respond to distinct stimuli. In this review, we discuss recent advances related to the roles of TFs in establishing the regulatory chromatin landscapes that orchestrate T cell development and differentiation. In particular, we focus on the role of TFs (e.g., TCF-1, BCL11B, PU.1, STAT3, STAT5, AP-1, and IRF4) in mediating chromatin accessibility and interactions and in regulating gene expression in T cells, including gene expression that is dependent on IL-2 and IL-21. Furthermore, we discuss the state of knowledge on enhancer-promoter interactions and how autoimmune disease risk variants can be linked to molecular functions of putative target genes.

Quantitative profiling of BATF family proteins/JUNB/IRF hetero-trimers using Spec-seq

Background

BATF family transcription factors (BATF, BATF2 and BATF3) form hetero-trimers with JUNB and either IRF4 or IRF8 to regulate cell fate in T cells and dendritic cells in vivo. While each combination of the hetero-trimer has a distinct role, some degree of cross-compensation was observed. The basis for the differential actions of IRF4 and IRF8 with BATF factors and JUNB is still unknown. We propose that the differences in function between these hetero-trimers may be caused by differences in their DNA binding preferences. While all three BATF family transcription factors have similar binding preferences when binding as a hetero-dimer with JUNB, the cooperative binding of IRF4 or IRF8 to the hetero-dimer/DNA complex could change the preferences. We used Spec-seq, which allows for the efficient and accurate determination of relative affinity to a large collection of sequences in parallel, to find differences between cooperative DNA binding of IRF4, IRF8 and BATF family members.

Results

We found that without IRF binding, all three hetero-dimer pairs exhibit nearly the same binding preferences to both expected wildtype binding sites TRE (TGA(C/G)TCA) and CRE (TGACGTCA). IRF4 and IRF8 show the very similar DNA binding preferences when binding with any of the three hetero-dimers. No major change of binding preferences was found in the half-sites between different hetero-trimers. IRF proteins bind with substantially lower affinity with either a single nucleotide spacer between IRF and BATF binding site or with an alternative mode of binding in the opposite orientation. In addition, the preference to CRE binding site was reduced with either IRF binding in all BATF–JUNB combinations.

Conclusions

The specificities of BATF, BATF2 and BATF3 are all very similar as are their interactions with IRF4 and IRF8. IRF proteins binding adjacent to BATF sites increases affinity substantially compared to sequences with spacings between the sites, indicating cooperative binding through protein–protein interactions. The preference for the type of BATF binding site, TRE or CRE, is also altered when IRF proteins bind. These in vitro preferences aid in the understanding of in vivo binding activities.

Electronic supplementary material

The online version of this article (10.1186/s12867-018-0106-7) contains supplementary material, which is available to authorized users.

Transcriptional orchestration of mitochondrial homeostasis in a cellular model of PGC-1-related coactivator-dependent thyroid tumor

The PGC-1 (Peroxisome proliferator-activated receptor Gamma Coactivator-1) family of coactivators (PGC-1α, PGC-1β, and PRC) plays a central role in the transcriptional control of mitochondrial biogenesis and oxidative phosphorylation (OXPHOS) processes. These coactivators integrate mitochondrial energy production into cell metabolism using complementary pathways. The XTC.UC1 cell line is a mitochondria-rich model of thyroid tumors whose biogenesis is almost exclusively dependent on PRC. Here we aim to propose an integrative view of the cellular pathways regulated by PRC through integration of cDNA and miRNA microarray data and chromatin immunoprecipitation results obtained from XTC.UC1 cells invalidated for PRC. This study showes that PRC induces a complex network of cellular functions interacting with at least one to five of the studied transcription factors (Estrogen Related Receptor alpha, ERR1; Nuclear-Respiratory Factors, NRF1 and NRF2; cAMP Response Element Binding, CREB; and Ying Yang, YY1). Our data confirm that ERR1 is a key partner of PRC in the regulation of mitochondrial functions and suggest a potential role of this complex in RNA processing. PRC is also involved in transcriptional regulatory complexes targeting 12 miRNAs, five of which are involved in the control of the OXPHOS process. Our findings demonstrate that the PRC coactivator can act in complex with several transcription factors and regulate miRNA expression to control the fine regulation of main metabolic functions in the cell. Therefore, in PGC-1α/β-associated pathologies, PRC, as a metabolic sensor, may ensure mitochondrial homeostasis.

Litchi Fruit LcNAC1 is a Target of LcMYC2 and Regulator of Fruit Senescence Through its Interaction with LcWRKY1

Senescence is a key factor resulting in deterioration of non-climacteric fruit. NAC transcription factors are important regulators in plant development and abiotic stress responses, yet little information regarding the role of NACs in regulating non-climacteric fruit senescence is available. In this study, we cloned 13 NAC genes from litchi (Litchi chinensis) fruit, and analyzed subcellular localization and expression profiles of these genes during post-harvest natural and low-temperature-delayed senescence. Of the 13 NAC genes, expression of LcNAC1 was up-regulated in the pericarp and pulp as senescence progressed, and was significantly higher in senescence-delayed fruit than that in naturally senescent fruit. LcNAC1 was induced by exogenous ABA and hydrogen peroxide. Yeast one-hybrid analysis and transient dual-luciferase reporter assay showed that LcNAC1 was positively regulated by the LcMYC2 transcription factor. LcNAC1 activated the expression of LcAOX1a, a gene associated with reactive oxygen species regulation and energy metabolism, whereas LcWRKY1 repressed LcAOX1a expression. In addition, LcNAC1 interacted with LcWRKY1 in vitro and in vivo. These results indicated that LcNAC1 and LcWRKY1 form a complex to regulate the expression of LcAOX1a antagonistically. Taken together, the results reveal a hierarchical and co-ordinated regulatory network in senescence of harvested litchi fruit.

Genomic feature extraction and comparison based on global alignment of ChIP-sequencing data

Enhanced accuracy and high-throughput capability in capturing genetic activities lead ChIP-sequencing technology to be applied prevalently in diverse study for tackling DNA-protein interaction problems. Till now, such questions as deciding suitable ChIP-seq arguments and comparing sample quality still haunt biologists. We propose the methods for answering such questions as deciding optimal argument pairs in global alignment of ChIP sequencing data; then we employ a modern signal processing approach to extract inherent genomic features from the global alignments of transcriptional binding activities; together with pairwise comparison from intra- and inter-sample perspectives; thus we can further determine alignment quality and decide the optimal candidate for multi-source heterogeneous high-throughput sequences. The work provides a practical approach to quantitatively compare the alignment quality for heterogeneous sequencing data, especially in determining the efficiency of transcriptional binding from replicate samples, thus it helps to exploit the potentiality of ChIP-seq for deep comprehension of inherent biological meanings from the high-throughput genomic sequences.

Abscisic-acid-dependent basic leucine zipper (bZIP) transcription factors in plant abiotic stress

One of the major causes of significant crop loss throughout the world is the myriad of environmental stresses including drought, salinity, cold, heavy metal toxicity, and ultraviolet-B (UV-B) rays. Plants as sessile organisms have evolved various effective mechanism which enable them to withstand this plethora of stresses. Most of such regulatory mechanisms usually follow the abscisic-acid (ABA)-dependent pathway. In this review, we have primarily focussed on the basic leucine zipper (bZIP) transcription factors (TFs) activated by the ABA-mediated signalosome. Upon perception of ABA by specialized receptors, the signal is transduced via various groups of Ser/Thr kinases, which phosphorylate the bZIP TFs. Following such post-translational modification of TFs, they are activated so that they bind to specific cis-acting sequences called abscisic-acid-responsive elements (ABREs) or GC-rich coupling elements (CE), thereby influencing the expression of their target downstream genes. Several in silico techniques have been adopted so far to predict the structural features, recognize the regulatory modification sites, undergo phylogenetic analyses, and facilitate genome-wide survey of TF under multiple stresses. Current investigations on the epigenetic regulation that controls greater accessibility of the inducible regions of DNA of the target gene to the bZIP TFs exclusively under stress situations, along with the evolved stress memory responses via genomic imprinting mechanism, have been highlighted. The potentiality of overexpression of bZIP TFs, either in a homologous or in a heterologous background, in generating transgenic plants tolerant to various abiotic stressors have also been addressed by various groups. The present review will provide a coherent documentation on the functional characterization and regulation of bZIP TFs under multiple environmental stresses, with the major goal of generating multiple-stress-tolerant plant cultivars in near future.

COPAR: A ChIP-Seq Optimal Peak Analyzer

Sequencing data quality and peak alignment efficiency of ChIP-sequencing profiles are directly related to the reliability and reproducibility of NGS experiments. Till now, there is no tool specifically designed for optimal peak alignment estimation and quality-related genomic feature extraction for ChIP-sequencing profiles. We developed open-sourced COPAR, a user-friendly package, to statistically investigate, quantify, and visualize the optimal peak alignment and inherent genomic features using ChIP-seq data from NGS experiments. It provides a versatile perspective for biologists to perform quality-check for high-throughput experiments and optimize their experiment design. The package COPAR can process mapped ChIP-seq read file in BED format and output statistically sound results for multiple high-throughput experiments. Together with three public ChIP-seq data sets verified with the developed package, we have deposited COPAR on GitHub under a GNU GPL license.

The banana fruit Dof transcription factor MaDof23 acts as a repressor and interacts with MaERF9 in regulating ripening-related genes

The DNA binding with one finger (Dof) proteins, a family of plant-specific transcription factors, are involved in a variety of plant biological processes. However, little information is available on their involvement in fruit ripening. We have characterized 25 MaDof genes from banana fruit (Musa acuminata), designated as MaDof1-MaDof25 Gene expression analysis in fruit subjected to different ripening conditions revealed that MaDofs were differentially expressed during different stages of ripening. MaDof10, 23, 24, and 25 were ethylene-inducible and nuclear-localized, and their transcript levels increased during fruit ripening. Moreover, yeast two-hybrid and bimolecular fluorescence complementation analyses demonstrated a physical interaction between MaDof23 and MaERF9, a potential regulator of fruit ripening reported in a previous study. We determined that MaDof23 is a transcriptional repressor, whereas MaERF9 is a transcriptional activator. We suggest that they might act antagonistically in regulating 10 ripening-related genes, including MaEXP1/2/3/5, MaXET7, MaPG1, MaPME3, MaPL2, MaCAT, and MaPDC, which are associated with cell wall degradation and aroma formation. Taken together, our findings provide new insight into the transcriptional regulation network controlling banana fruit ripening.

Methylated Cytosines Mutate to Transcription Factor Binding Sites that Drive Tetrapod Evolution

In mammals, the cytosine in CG dinucleotides is typically methylated producing 5-methylcytosine (5mC), a chemically less stable form of cytosine that can spontaneously deaminate to thymidine resulting in a T•G mismatched base pair. Unlike other eukaryotes that efficiently repair this mismatched base pair back to C•G, in mammals, 5mCG deamination is mutagenic, sometimes producing TG dinucleotides, explaining the depletion of CG dinucleotides in mammalian genomes. It was suggested that new TG dinucleotides generate genetic diversity that may be critical for evolutionary change. We tested this conjecture by examining the DNA sequence properties of regulatory sequences identified by DNase I hypersensitive sites (DHSs) in human and mouse genomes. We hypothesized that the new TG dinucleotides generate transcription factor binding sites (TFBS) that become tissue-specific DHSs (TS-DHSs). We find that 8-mers containing the CG dinucleotide are enriched in DHSs in both species. However, 8-mers containing a TG and no CG dinucleotide are preferentially enriched in TS-DHSs when compared with 8-mers with neither a TG nor a CG dinucleotide. The most enriched 8-mer with a TG and no CG dinucleotide in tissue-specific regulatory regions in both genomes is the AP-1 motif (TGA^C/_GTCAN), and we find evidence that TG dinucleotides in the AP-1 motif arose from CG dinucleotides. Additional TS-DHS-enriched TFBS containing the TG/CA dinucleotide are the E-Box motif (GCAGCTGC), the NF-1 motif (GGCA—TGCC), and the GR (glucocorticoid receptor) motif (G-ACA—TGT-C). Our results support the suggestion that cytosine methylation is mutagenic in tetrapods producing TG dinucleotides that create TFBS that drive evolution.

GABPα Binding to Overlapping ETS and CRE DNA Motifs Is Enhanced by CREB1: Custom DNA Microarrays

To achieve proper spatiotemporal control of gene expression, transcription factors cooperatively assemble onto specific DNA sequences. The ETS domain protein monomer of GABPα and the B-ZIP domain protein dimer of CREB1 cooperatively bind DNA only when the ETS (^C/_GCGGAAGT) and CRE (GTGACGTCAC) motifs overlap precisely, producing the ETS↔CRE motif (^C/_GCGGAAGTGACGTCAC). We designed a Protein Binding Microarray (PBM) with 60-bp DNAs containing four identical sectors, each with 177,440 features that explore the cooperative interactions between GABPα and CREB1 upon binding the ETS↔CRE motif. The DNA sequences include all 15-mers of the form ^C/_GCGGA—–CG—, the ETS↔CRE motif, and all single nucleotide polymorphisms (SNPs), and occurrences in the human and mouse genomes. CREB1 enhanced GABPα binding to the canonical ETS↔CRE motif CCGGAAGT two-fold, and up to 23-fold for several SNPs at the beginning and end of the ETS motif, which is suggestive of two separate and distinct allosteric mechanisms of cooperative binding. We show that the ETS-CRE array data can be used to identify regions likely cooperatively bound by GABPα and CREB1 in vivo, and demonstrate their ability to identify human genetic variants that might inhibit cooperative binding.

A widespread role of the motif environment in transcription factor binding across diverse protein families

Transcriptional regulation requires the binding of transcription factors (TFs) to short sequence-specific DNA motifs, usually located at the gene regulatory regions. Interestingly, based on a vast amount of data accumulated from genomic assays, it has been shown that only a small fraction of all potential binding sites containing the consensus motif of a given TF actually bind the protein. Recent in vitro binding assays, which exclude the effects of the cellular environment, also demonstrate selective TF binding. An intriguing conjecture is that the surroundings of cognate binding sites have unique characteristics that distinguish them from other sequences containing a similar motif that are not bound by the TF. To test this hypothesis, we conducted a comprehensive analysis of the sequence and DNA shape features surrounding the core-binding sites of 239 and 56 TFs extracted from in vitro HT-SELEX binding assays and in vivo ChIP-seq data, respectively. Comparing the nucleotide content of the regions around the TF-bound sites to the counterpart unbound regions containing the same consensus motifs revealed significant differences that extend far beyond the core-binding site. Specifically, the environment of the bound motifs demonstrated unique sequence compositions, DNA shape features, and overall high similarity to the core-binding motif. Notably, the regions around the binding sites of TFs that belong to the same TF families exhibited similar features, with high agreement between the in vitro and in vivo data sets. We propose that these unique features assist in guiding TFs to their cognate binding sites.

Circadian clocks: the tissue is the issue

The circadian clock uses a widely expressed pair of clock activators to drive tissue-specific rhythms in target gene expression. A new study sheds light on this tissue specificity by showing that binding of clock activators and tissue-specific transcription factors to closely associated target sites enables cooperative activation of target genes in different tissues.

Homeostatic control of regulatory T cell diversity

Regulatory T (TReg) cells constitute an essential counterbalance to adaptive immune responses. Failure to maintain appropriate TReg cell numbers or function leads to autoimmune, malignant and immunodeficient conditions. Dynamic homeostatic processes preserve the number of forkhead box P3-expressing (FOXP3(+)) TReg cells within a healthy range, with high rates of cell division being offset by apoptosis under steady-state conditions. Recent studies have shown that TReg cells become specialized for different environmental contexts, tailoring their functions and homeostatic properties to a wide range of tissues and immune conditions. In this Review, we describe new insights into the molecular controls that maintain the steady-state homeostasis of TReg cells and the cues that drive TReg cell adaptation to inflammation and/or different locations. We highlight how differing local milieu might drive context-specific TReg cell function and restoration of immune homeostasis, and how dysregulation of these processes can precipitate disease.

The Arabidopsis NAC transcription factor ANAC096 cooperates with bZIP-type transcription factors in dehydration and osmotic stress responses

Multiple transcription factors (TFs) play essential roles in plants under abiotic stress, but how these multiple TFs cooperate in abiotic stress responses remains largely unknown. In this study, we provide evidence that the NAC (for NAM, ATAF1/2, and CUC2) TF ANAC096 cooperates with the bZIP-type TFs ABRE binding factor and ABRE binding protein (ABF/AREB) to help plants survive under dehydration and osmotic stress conditions. ANAC096 directly interacts with ABF2 and ABF4, but not with ABF3, both in vitro and in vivo. ANAC096 and ABF2 synergistically activate RD29A transcription. Our genome-wide gene expression analysis revealed that a major proportion of abscisic acid (ABA)-responsive genes are under the transcriptional regulation of ANAC096. We found that the Arabidopsis thaliana anac096 mutant is hyposensitive to exogenous ABA and shows impaired ABA-induced stomatal closure and increased water loss under dehydration stress conditions. Furthermore, we found the anac096 abf2 abf4 triple mutant is much more sensitive to dehydration and osmotic stresses than the anac096 single mutant or the abf2 abf4 double mutant. Based on these results, we propose that ANAC096 is involved in a synergistic relationship with a subset of ABFs for the transcriptional activation of ABA-inducible genes in response to dehydration and osmotic stresses.

Regulators of chromatin state and transcription in CD4 T-cell polarization

Mature naive CD4 T-cells possess the potential for an array of highly specialized functions, from inflammatory to potently suppressive. This potential is encoded in regulatory DNA elements and is fulfilled through modification of chromatin and selective activation by the collaborative function of diverse transcription factors in response to environmental cues. The mechanisms and strategies employed by transcription factors for the programming of CD4 T-cell subsets will be discussed. In particular, the focus will be on co-operative activity of environmental response factors in the initial activation of regulatory DNA elements and chromatin alteration, and the subsequent role of 'master regulator' transcription factors in defining the fidelity and environmental responsiveness of different CD4 T-cell subsets.

Contribution of nucleosome binding preferences and co-occurring DNA sequences to transcription factor binding

BACKGROUND

Chromatin plays a critical role in regulating transcription factors (TFs) binding to their canonical transcription factor binding sites (TFBS). Recent studies in vertebrates show that many TFs preferentially bind to genomic regions that are well bound by nucleosomes in vitro. Co-occurring secondary motifs sometimes correlated with functional TFBS.

RESULTS

We used a logistic regression to evaluate how well the propensity for nucleosome binding and co-occurrence of a secondary motif identify which canonical motifs are bound in vivo. We used ChIP-seq data for three transcription factors binding to their canonical motifs: c-Jun binding the AP-1 motif (TGA(C)/(G)TCA), GR (glucocorticoid receptor) binding the GR motif (G-ACA---(T)/(C)GT-C), and Hoxa2 (homeobox a2) binding the Pbx (Pre-B-cell leukemia homeobox) motif (TGATTGAT). For all canonical TFBS in the mouse genome, we calculated intrinsic nucleosome occupancy scores (INOS) for its surrounding 150-bps DNA and examined the relationship with in vivo TF binding. In mouse mammary 3134 cells, c-Jun and GR proteins preferentially bound regions calculated to be well-bound by nucleosomes in vitro with the canonical AP-1 and GR motifs themselves contributing to the high INOS. Functional GR motifs are enriched for AP-1 motifs if they are within a nucleosome-sized 150-bps region. GR and Hoxa2 also bind motifs with low INOS, perhaps indicating a different mechanism of action.

CONCLUSION

Our analysis quantified the contribution of INOS and co-occurring sequence to the identification of functional canonical motifs in the genome. This analysis revealed an inherent competition between some TFs and nucleosomes for binding canonical TFBS. GR and c-Jun cooperate if they are within 150-bps. Binding of Hoxa2 and a fraction of GR to motifs with low INOS values suggesting they are not in competition with nucleosomes and may function using different mechanisms.