Clustering of DNA sequences in human promoters

Induction of the human CDC45 gene promoter activity by natural compound trans‑resveratrol

GGAA motifs in the human TP53 and HELB gene promoters play a part in responding to trans‑resveratrol (Rsv) in HeLa S3 cells. This sequence is also present in the 5'‑upstream region of the human CDC45 gene, which encodes a component of CMG DNA helicase protein complex. The cells were treated with Rsv (20 µM), then transcripts and the translated protein were analyzed by quantitative RT‑PCR and western blotting, respectively. The results showed that the CDC45 gene and protein expression levels were induced after the treatment. To examine whether they were due to the activation of transcription, a 5'‑upstream 556‑bp of the CDC45 gene was cloned and inserted into a multi‑cloning site of the Luciferase (Luc) expression vector. In the present study, various deletion/point mutation‑introduced Luc expression plasmids were constructed and they were used for the transient transfection assay. The results showed that the GGAA motif, which is included in a putative RELB protein recognizing sequence, plays a part in the promoter activity with response to Rsv in HeLa S3 cells.

Lineage-determining transcription factor-driven promoters regulate cell type-specific macrophage gene expression

Mammalian promoters consist of multifarious elements, which make them unique and support the selection of the proper transcript variants required under diverse conditions in distinct cell types. However, their direct DNA-transcription factor (TF) interactions are mostly unidentified. Murine bone marrow-derived macrophages (BMDMs) are a widely used model for studying gene expression regulation. Thus, this model serves as a rich source of various next-generation sequencing data sets, including a large number of TF cistromes. By processing and integrating the available cistromic, epigenomic and transcriptomic data from BMDMs, we characterized the macrophage-specific direct DNA-TF interactions, with a particular emphasis on those specific for promoters. Whilst active promoters are enriched for certain types of typically methylatable elements, more than half of them contain non-methylatable and prototypically promoter-distal elements. In addition, circa 14% of promoters-including that of Csf1r-are composed exclusively of 'distal' elements that provide cell type-specific gene regulation by specialized TFs. Similar to CG-rich promoters, these also contain methylatable CG sites that are demethylated in a significant portion and show high polymerase activity. We conclude that this unusual class of promoters regulates cell type-specific gene expression in macrophages, and such a mechanism might exist in other cell types too.

Therapeutic Potential of Clusterin Inhibition in Human Cancer

Clusterin (CLU) protein is involved in various pathophysiological processes including carcinogenesis and tumor progression. In recent years, the role of the secretory isoform has been demonstrated in tumor cells, where it inhibits apoptosis and favors the acquisition of resistance to conventional treatments used to treat cancer. To determine the possible therapeutic potential of inhibiting this protein, numerous studies have been carried out in this field. In this article, we present the existing knowledge to date on the inhibition of this protein in different types of cancer and analyze the importance it could have in the development of new therapies targeted against this disease.

Regulatory effect of NK homeobox 1 (NKX2.1) on melanocortin 4 receptor (Mc4r) promoter in Mandarin fish

The melanocortin 4 receptor (MC4R) is a G protein-coupled transporter that mediates the regulation of thyroid hormones and leptin on energy balance and food intake. However, the mechanisms of transcriptional regulation of Mc4r by thyroid hormone and leptin in fish have been rarely reported. The messenger RNA expression of Mc4r gene was significantly higher in brain than those in other tissues of mandarin fish. We analyzed the structure and function of a 2029 bp sequence of Mc4r promoter. Meanwhile, overexpression of NKX2.1 and incubation with leptin significantly increased Mc4r promoter activity, but triiodothyronine showed the opposite effect. In addition, mutations in the NKX2.1 binding site abolished not only the activation of Mc4r promoter activity by leptin but also the inhibitory effect of thyroid hormones on Mc4r promoter activity. In summary, these results suggested that thyroid hormones and leptin might regulate the transcriptional expression of Mc4r through NKX2.1.

Transcription factor IID parks and drives preinitiation complexes at sharp or broad promoters

Core promoters are sites where transcriptional regulatory inputs of a gene are integrated to direct the assembly of the preinitiation complex (PIC) and RNA polymerase II (Pol II) transcription output. Until now, core promoter functions have been investigated by distinct methods, including Pol II transcription initiation site mappings and structural characterization of PICs on distinct promoters. Here, we bring together these previously unconnected observations and hypothesize how, on metazoan TATA promoters, the precisely structured building up of transcription factor (TF) IID-based PICs results in sharp transcription start site (TSS) selection; or, in contrast, how the less strictly controlled positioning of the TATA-less promoter DNA relative to TFIID-core PIC components results in alternative broad TSS selections by Pol II.

Widespread perturbation of ETS factor binding sites in cancer

Although >90% of somatic mutations reside in non-coding regions, few have been reported as cancer drivers. To predict driver non-coding variants (NCVs), we present a transcription factor (TF)-aware burden test based on a model of coherent TF function in promoters. We apply this test to NCVs from the Pan-Cancer Analysis of Whole Genomes cohort and predict 2555 driver NCVs in the promoters of 813 genes across 20 cancer types. These genes are enriched in cancer-related gene ontologies, essential genes, and genes associated with cancer prognosis. We find that 765 candidate driver NCVs alter transcriptional activity, 510 lead to differential binding of TF-cofactor regulatory complexes, and that they primarily impact the binding of ETS factors. Finally, we show that different NCVs within a promoter often affect transcriptional activity through shared mechanisms. Our integrated computational and experimental approach shows that cancer NCVs are widespread and that ETS factors are commonly disrupted.

Gall-specific promoter, an alternative to the constitutive CaMV35S promoter, drives host-derived RNA interference targeting Mi-msp2 gene to confer effective nematode resistance

One of the major obligate plant parasites causing massive economic crop losses belongs to the class of root-knot nematodes (RKNs). Targeting of major nematode parasitism genes via Host Delivered-RNAi (HD-RNAi) to confer silencing is established as one of the most effective approaches to curb nematode infection. Utilizing nematode-responsive root-specific (NRRS) promoters to design a dsRNA molecule targeting approach to hamper nematode parasitism. Here, a previously validated peroxidase gall specific promoter, pAt2g18140, from Arabidopsis was employed to express the dsRNA construct of the nematode effector gene Mi-msp2 from Meloidogyne incognita. Arabidopsis RNAi lines of CaMV35S::Mi-msp2-RNAi and pAt2g18140::Mi-msp2-RNAi were compared with control plants to assess the decrease in plant nematode infection. When subjected to infection, the maximum reductions in the numbers of galls, females and egg masses in the CaMV35S::Mi-msp2-RNAi lines were 61%, 66% and 95%, respectively, whereas for the pAt2g18140::Mi-msp2-RNAi lines, they were 63%, 68% and 100%, respectively. The reduction in transcript level ranged from 79%-82% for CaMV35S::Mi-msp2-RNAi and 72%-79% for the pAt2g18140::Mi-msp2-RNAi lines. Additionally, a reduction in female size and a subsequent reduction in next-generation fecundity demonstrate the efficacy and potential of the gall specific promoter pAt2g18140 for utilization in the development of HD-RNAi constructs against RKN, as an excellent alternative to the CaMV35S promoter.

Transposable Elements: Distribution, Polymorphism, and Climate Adaptation in Populus

Transposable elements (TEs) are a class of mobile genetic elements that make effects on shaping rapid phenotypic traits of adaptive significance. TE insertions are usually related to transcription changes of nearby genes, and thus may be subjected to purifying selection. Based on the available genome resources of Populus, we found that the composition of Helitron DNA family were highly variable and could directly influence the transcription of nearby gene expression, which are involving in stress-responsive, programmed cell death, and apoptosis pathway. Next, we indicated TEs are highly enriched in Populus trichocarpa compared with three other congeneric poplar species, especially located at untranslated regions (3'UTRs and 5'UTRs) and Helitron transposons, particularly 24-nt siRNA-targeted, are significantly associated with reduced gene expression. Additionally, we scanned a representative resequenced Populus tomentosa population, and identified 9,680 polymorphic TEs loci. More importantly, we identified a Helitron transposon located at the 3'UTR, which could reduce WRKY18 expression level. Our results highlight the importance of TE insertion events, which could regulate gene expression and drive adaptive phenotypic variation in Populus.

BANP opens chromatin and activates CpG-island-regulated genes

The majority of gene transcripts generated by RNA polymerase II in mammalian genomes initiate at CpG island (CGI) promoters^1,2, yet our understanding of their regulation remains limited. This is in part due to the incomplete information that we have on transcription factors, their DNA-binding motifs and which genomic binding sites are functional in any given cell type^3-5. In addition, there are orphan motifs without known binders, such as the CGCG element, which is associated with highly expressed genes across human tissues and enriched near the transcription start site of a subset of CGI promoters^6-8. Here we combine single-molecule footprinting with interaction proteomics to identify BTG3-associated nuclear protein (BANP) as the transcription factor that binds this element in the mouse and human genome. We show that BANP is a strong CGI activator that controls essential metabolic genes in pluripotent stem and terminally differentiated neuronal cells. BANP binding is repelled by DNA methylation of its motif in vitro and in vivo, which epigenetically restricts most binding to CGIs and accounts for differential binding at aberrantly methylated CGI promoters in cancer cells. Upon binding to an unmethylated motif, BANP opens chromatin and phases nucleosomes. These findings establish BANP as a critical activator of a set of essential genes and suggest a model in which the activity of CGI promoters relies on methylation-sensitive transcription factors that are capable of chromatin opening.

Characterization of the human TARDBP gene promoter

The expression of TDP-43, the main component of neuronal intracellular inclusions across a broad spectrum of ALS and FTD disorders, is developmentally regulated and studies in vivo have shown that TDP-43 overexpression can be toxic, even before observation of pathological aggregates. Starting from these observations, the regulation of its expression at transcriptional level might represent a further key element for the pathogenesis of neurodegenerative diseases. Therefore, we have characterized the human TARDBP promoter, in order to study the transcriptional mechanisms of expression. Mapping of cis-acting elements by luciferase assays in different cell outlined that the activity of the promoter seems to be higher in SH-SY5Y, Neuro2A, and HeLa than in HEK293. In addition, we tested effects of two SNPs found in the promoter region of ALS patients and observed no significant effect on transcription levels in all tested cell lines. Lastly, while TDP-43 overexpression did not affect significantly the activity of its promoter (suggesting that TDP-43 does not influence its own transcription), the presence of the 5'UTR sequence and of intron-1 splicing seem to impact positively on TDP-43 expression without affecting transcript stability. In conclusion, we have identified the region spanning nucleotides 451-230 upstream from the transcription start site as the minimal region with a significant transcription activity. These results lay an important foundation for exploring the regulation of the TARDBP gene transcription by exogenous and endogenous stimuli and the implication of transcriptional mechanisms in the pathogenesis of TDP-43 proteinopathies.

Dysregulation of PGC-1α-Dependent Transcriptional Programs in Neurological and Developmental Disorders: Therapeutic Challenges and Opportunities

Substantial evidence indicates that mitochondrial impairment contributes to neuronal dysfunction and vulnerability in disease states, leading investigators to propose that the enhancement of mitochondrial function should be considered a strategy for neuroprotection. However, multiple attempts to improve mitochondrial function have failed to impact disease progression, suggesting that the biology underlying the normal regulation of mitochondrial pathways in neurons, and its dysfunction in disease, is more complex than initially thought. Here, we present the proteins and associated pathways involved in the transcriptional regulation of nuclear-encoded genes for mitochondrial function, with a focus on the transcriptional coactivator peroxisome proliferator-activated receptor gamma coactivator-1alpha (PGC-1α). We highlight PGC-1α’s roles in neuronal and non-neuronal cell types and discuss evidence for the dysregulation of PGC-1α-dependent pathways in Huntington’s Disease, Parkinson’s Disease, and developmental disorders, emphasizing the relationship between disease-specific cellular vulnerability and cell-type-specific patterns of PGC-1α expression. Finally, we discuss the challenges inherent to therapeutic targeting of PGC-1α-related transcriptional programs, considering the roles for neuron-enriched transcriptional coactivators in co-regulating mitochondrial and synaptic genes. This information will provide novel insights into the unique aspects of transcriptional regulation of mitochondrial function in neurons and the opportunities for therapeutic targeting of transcriptional pathways for neuroprotection.

Dysregulation of PGC-1α-Dependent Transcriptional Programs in Neurological and Developmental Disorders: Therapeutic Challenges and Opportunities

Substantial evidence indicates that mitochondrial impairment contributes to neuronal dysfunction and vulnerability in disease states, leading investigators to propose that the enhancement of mitochondrial function should be considered a strategy for neuroprotection. However, multiple attempts to improve mitochondrial function have failed to impact disease progression, suggesting that the biology underlying the normal regulation of mitochondrial pathways in neurons, and its dysfunction in disease, is more complex than initially thought. Here, we present the proteins and associated pathways involved in the transcriptional regulation of nuclear-encoded genes for mitochondrial function, with a focus on the transcriptional coactivator peroxisome proliferator-activated receptor gamma coactivator-1alpha (PGC-1α). We highlight PGC-1α's roles in neuronal and non-neuronal cell types and discuss evidence for the dysregulation of PGC-1α-dependent pathways in Huntington's Disease, Parkinson's Disease, and developmental disorders, emphasizing the relationship between disease-specific cellular vulnerability and cell-type-specific patterns of PGC-1α expression. Finally, we discuss the challenges inherent to therapeutic targeting of PGC-1α-related transcriptional programs, considering the roles for neuron-enriched transcriptional coactivators in co-regulating mitochondrial and synaptic genes. This information will provide novel insights into the unique aspects of transcriptional regulation of mitochondrial function in neurons and the opportunities for therapeutic targeting of transcriptional pathways for neuroprotection.

On the NF-Y regulome as in ENCODE (2019)

NF-Y is a trimeric Transcription Factor -TF- which binds with high selectivity to the conserved CCAAT element. Individual ChIP-seq analysis as well as ENCODE have progressively identified locations shared by other TFs. Here, we have analyzed data introduced by ENCODE over the last five years in K562, HeLa-S3 and GM12878, including several chromatin features, as well RNA-seq profiling of HeLa cells after NF-Y inactivation. We double the number of sequence-specific TFs and co-factors reported. We catalogue them in 4 classes based on co-association criteria, infer target genes categorizations, identify positional bias of binding sites and gene expression changes. Larger and novel co-associations emerge, specifically concerning subunits of repressive complexes as well as RNA-binding proteins. On the one hand, these data better define NF-Y association with single members of major classes of TFs, on the other, they suggest that it might have a wider role in the control of mRNA production.

The ongoing ENCODE consortium represents a useful compendium of locations of TFs, chromatin marks, gene expression data. In previous reports, we identified modules of CCAAT-binding NF-Y with individual TFs. Here, we analyzed all 363 factors currently present: 68 with enrichment of CCAAT in their locations, 38 with overlap of peaks. New sequence-specific TFs, co-activators and co-repressors are reported. Co-association patterns correspond to specific targeted genes categorizations and gene expression changes, as assessed by RNA-seq after NF-Y inactivation. These data widen and better define a coherent model of synergy of NF-Y with selected groups of TFs and co-factors.

A reference catalog of DNA palindromes in the human genome and their variations in 1000 Genomes

A palindrome in DNA is like a palindrome in language, but when read backwards, it is a complement of the forward sequence; effectively, the two halves of a sequence complement each other from its midpoint like in a double strand of DNA. Palindromes are distributed throughout the human genome and play significant roles in gene expression and regulation. Palindromic mutations are linked to many human diseases, such as neuronal disorders, mental retardation, and various cancers. In this work, we computed and analyzed the palindromic sequences in the human genome and studied their conservation in personal genomes using 1000 Genomes data. We found that ~30% of the palindromes exhibit variation, some of which are caused by rare variants. The analysis of disease/trait-associated single-nucleotide polymorphisms in palindromic regions showed that disease-associated risk variants are 14 times more likely to be present in palindromic regions than in other regions. The catalog of palindromes in the reference genome and 1000 Genomes is being made available here with details on their variations in each individual genome to serve as a resource for future and retrospective whole-genome studies identifying statistically significant palindrome variations associated with diseases or traits and their roles in disease mechanisms.

NRF-1 and HIF-1α contribute to modulation of human VDAC1 gene promoter during starvation and hypoxia in HeLa cells

VDAC (Voltage Dependent Anion Channel) is a family of pore forming protein located in the outer mitochondrial membrane. Its channel property ensures metabolites exchange between mitochondria and the rest of the cell resulting in metabolism and bioenergetics regulation, and in cell death and life switch. VDAC1 is the best characterized and most abundant isoform, and is involved in many pathologies, as cancer or neurodegenerative diseases. However, little information is available about its gene expression regulation in normal and/or pathological conditions. In this work, we explored VDAC1 gene expression regulation in normal conditions and in the contest of some metabolic and energetic mitochondrial dysfunction and cell stress as example. The core of the putative promoter region was characterized in terms of transcription factors responsive elements both by bioinformatic studies and promoter activity experiments. In particular, we found an abundant presence of NRF-1 sites, together with other transcription factors binding sites involved in cell growth, proliferation, development, and we studied their prevalence in gene activity. Furthermore, upon depletion of nutrients or controlled hypoxia, as detected in various pathologies, we found that VDAC1 transcripts levels were significantly increased in a time related manner. VDAC1 promoter activity was also validated by gene reporter assays. According to PCR real-time experiments, it was confirmed that VDAC1 promoter activity is further stimulated when cells are exposed to stress. A bioinformatic survey suggested HIF-1α, besides NRF-1, as a most active TFBS. Their validation was obtained by TFBS mutagenesis and TF overexpression experiments. In conclusion, we experimentally demonstrated the involvement of both NRF-1 and HIF-1α in the regulation of VDAC1 promoter activation at basal level and in some peculiar cell stress conditions.

Bioinformatic analysis reveals MIR502 as a potential tumour suppressor in ovarian cancer

BACKGROUND

Ovarian cancer (OC) is a major cause of death among women due to the lack of early screening methods and its complex pathological progression. Increasing evidence has indicated that microRNAs regulate gene expression in tumours by interacting with mRNAs. Although the research regarding OC and microRNAs is extensive, the vital role of MIR502 in OC remains unclear.

METHODS

We integrated two microRNA expression arrays from GEO to identify differentially expressed genes. The Kaplan-Meier method was used to screen for miRNAs that had an influence on survival outcome. Upstream regulators of MIR502 were predicted by JASPAR and verified by ChIP-seq data. The LinkedOmics database was used to study genes that were correlated with MIR502. Gene Set Enrichment Analysis (GSEA) was conducted for functional annotation with GO and KEGG pathway enrichment analyses by using the open access WebGestalt tool. We constructed a PPI network by using STRING to further explore the core proteins.

RESULTS

We found that the expression level of MIR502 was significantly downregulated in OC, which was related to poor overall survival. NRF1, as an upstream regulator of MIR502, was predicted by JASPAR and verified by ChIP-seq data. In addition, anti-apoptosis and pro-proliferation genes in the Hippo signalling pathway, including CCND1, MYC, FGF1 and GLI2, were negatively regulated by MIR502, as shown in the GO and KEGG pathway enrichment results. The PPI network further demonstrated that CCND1 and MYCN were at core positions in the development of ovarian cancer.

CONCLUSIONS

MIR502, which is regulated by NRF1, acts as a tumour suppressor gene to accelerate apoptosis and suppress proliferation by targeting the Hippo signalling pathway in ovarian cancer.

Transcriptional regulatory network controlling the ontogeny of hematopoietic stem cells

In this study from Gao et al., the authors performed RNA-seq and histone mark ChIP-seq to define the transcriptomes and epigenomes of cells representing key developmental stages of HSC ontogeny in mice. Using a novel computational algorithm, target inference via physical connection (TIPC), they constructed developmental stage-specific transcriptional regulatory networks by linking enhancers and predicted bound transcription factors to their target promoters, thus providing a useful resource for uncovering regulators of HSC formation.

Hematopoietic stem cell (HSC) ontogeny is accompanied by dynamic changes in gene regulatory networks. We performed RNA-seq and histone mark ChIP-seq to define the transcriptomes and epigenomes of cells representing key developmental stages of HSC ontogeny in mice. The five populations analyzed were embryonic day 10.5 (E10.5) endothelium and hemogenic endothelium from the major arteries, an enriched population of prehematopoietic stem cells (pre-HSCs), fetal liver HSCs, and adult bone marrow HSCs. Using epigenetic signatures, we identified enhancers for each developmental stage. Only 12% of enhancers are primed, and 78% are active, suggesting the vast majority of enhancers are established de novo without prior priming in earlier stages. We constructed developmental stage-specific transcriptional regulatory networks by linking enhancers and predicted bound transcription factors to their target promoters using a novel computational algorithm, target inference via physical connection (TIPC). TIPC predicted known transcriptional regulators for the endothelial-to-hematopoietic transition, validating our overall approach, and identified putative novel transcription factors, including the broadly expressed transcription factors SP3 and MAZ. Finally, we validated a role for SP3 and MAZ in the formation of hemogenic endothelium. Our data and computational analyses provide a useful resource for uncovering regulators of HSC formation.

The landscape of transcription initiation across latent and lytic KSHV genomes

Precise promoter annotation is required for understanding the mechanistic basis of transcription initiation. In the context of complex genomes, such as herpesviruses where there is extensive genic overlap, identification of transcription start sites (TSSs) is particularly problematic and cannot be comprehensively accessed by standard RNA sequencing approaches. Kaposi's sarcoma-associated herpesvirus (KSHV) is an oncogenic gammaherpesvirus and the etiological agent of Kaposi's sarcoma and the B cell lymphoma primary effusion lymphoma (PEL). Here, we leverage RNA annotation and mapping of promoters for analysis of gene expression (RAMPAGE) and define KSHV TSSs transcriptome-wide and at nucleotide resolution in two widely used models of KSHV infection, namely iSLK.219 cells and the PEL cell line TREx-BCBL1-RTA. By mapping TSSs over a 96 h time course of reactivation we confirm 48 of 50 previously identified TSSs. Moreover, we identify over 100 novel transcription start site clusters (TSCs) in each cell line. Our analyses identified cell-type specific differences in TSC positions as well as promoter strength, and defined motifs within viral core promoters. Collectively, by defining TSSs at high resolution we have greatly expanded the transcriptional landscape of the KSHV genome and identified transcriptional control mechanisms at play during KSHV lytic reactivation.

Regulation of differential proton-coupled folate transporter gene expression in human tumors: transactivation by KLF15 with NRF-1 and the role of Sp1

Tumors can be therapeutically targeted with novel antifolates (e.g. AGF94) that are selectively transported by the human proton-coupled folate transporter (hPCFT). Studies were performed to determine the transcription regulation of hPCFT in tumors and identify possible mechanisms that contribute to the highly disparate levels of hPCFT in HepG2 versus HT1080 tumor cells. Transfection of hPCFT-null HT1080 cells with hPCFT restored transport and sensitivity to AGF94 Progressive deletions of the hPCFT promoter construct (-2005 to +96) and reporter gene assays in HepG2 and HT1080 cells confirmed differences in hPCFT transactivation and localized a minimal promoter to between positions -50 and +96. The minimal promoter included KLF15, GC-Box and NRF-1 cis-binding elements whose functional importance was confirmed by promoter deletions and mutations of core consensus sequences and reporter gene assays. In HepG2 cells, NRF-1, KLF15 and Sp1 transcripts were increased over HT1080 cells by ∼5.1-, ∼44-, and ∼2.4-fold, respectively. In Drosophila SL2 cells, transfection with KLF15 and NRF-1 synergistically activated the hPCFT promoter; Sp1 was modestly activating or inhibitory. Chromatin immunoprecipitation and electrophoretic mobility shift assay (EMSA) and supershifts confirmed differential binding of KLF15, Sp1, and NRF-1 to the hPCFT promoter in HepG2 and HT1080 cells that paralleled hPCFT levels. Treatment of HT1080 nuclear extracts (NE) with protein kinase A increased Sp1 binding to its consensus sequence by EMSA, suggesting a role for Sp1 phosphorylation in regulating hPCFT transcription. A better understanding of determinants of hPCFT transcriptional control may identify new therapeutic strategies for cancer by modulating hPCFT levels in combination with hPCFT-targeted antifolates.

Bioinformatic and functional analysis of promoter region of human SLC25A13 gene

The human SLC25A13 gene encodes the liver type aspartate/glutamate carrier isoform 2 (AGC2, commonly named as citrin), which plays a key role in the main NADH-shuttle of human hepatocyte. Biallelic SLC25A13 mutations result in Citrin deficiency (CD). In order to identify the important regulatory region of SLC25A13 gene and elucidate the way how potential promoter mutations affect the citrin expression, we performed promoter deletion analysis and established the reporter constructs of luciferase gene-carrying SLC25A13 promoter containing several mutations located in putative transcription factor-binding sites. The luciferase activities of all promoter constructs were measured using a Dual-Luciferase Reporter Assay System. Bioinformatic analysis showed that the promoter of SLC25A13 gene lacks TATA box and obviously typical initiator element, but contains a CCAAT box and two GC box. Promoter deletion analysis confirmed the region from -221 to -1 upstream ATG was essential for SLC25A13 to maintain the promoter activity. We utilized dual-luciferase reporter system as function analytical model to tentatively assess the effect of artificially constructed promoter mutations on citrin expression, and our analysis revealed that mutated putative CCAAT box and GC box could significantly affect the citrin expression. Our study confirmed the important SLC25A13 promoter regions that influenced citrin expression in HL7702 cells, and constructed a function analytical model. This work may be useful to further identify the pathogenic mutations leading to CD in the promoter region.

Structure of human TFIID and mechanism of TBP loading onto promoter DNA

The general transcription factor IID (TFIID) is a critical component of the eukaryotic transcription preinitiation complex (PIC) and is responsible for recognizing the core promoter DNA and initiating PIC assembly. We used cryo-electron microscopy, chemical cross-linking mass spectrometry, and biochemical reconstitution to determine the complete molecular architecture of TFIID and define the conformational landscape of TFIID in the process of TATA box-binding protein (TBP) loading onto promoter DNA. Our structural analysis revealed five structural states of TFIID in the presence of TFIIA and promoter DNA, showing that the initial binding of TFIID to the downstream promoter positions the upstream DNA and facilitates scanning of TBP for a TATA box and the subsequent engagement of the promoter. Our findings provide a mechanistic model for the specific loading of TBP by TFIID onto the promoter.

Classifying human promoters by occupancy patterns identifies recurring sequence elements, combinatorial binding, and spatial interactions

Background

Characterizing recurring sequence patterns in human promoters has been a challenging undertaking even nowadays where a near-complete overview of promoters exists. However, with the more recent availability of genomic location (ChIP-seq) data, one can approach that question through the identification of characteristic patterns of transcription factor occupancy and histone modifications.

Results

Based on the ENCODE annotation and integration of sequence motifs as well as three-dimensional chromatin data, we have undertaken a re-analysis of occupancy and sequence patterns in human promoters. We identify clear groups of CAAT-box and E-box sequence motif containing promoters, as well as a group of promoters whose interaction with an enhancer appears to be mediated by CCCTC-binding factor (CTCF) binding on the promoter. We also extend our analysis to inactive promoters, showing that only a surprisingly small number of inactive promoters is repressed by the polycomb complex. We also identify combinatorial patterns of transcription factor interactions indicated by the ChIP-seq signals.

Conclusion

Our analysis defines subgroups of promoters characterized by stereotypic patterns of transcription factor occupancy, and combinations of specific sequence patterns which are required for their binding. This grouping provides new hypotheses concerning the assembly and dynamics of transcription factor complexes at their respective promoter groups, as well as questions on the evolutionary origin of these groups.

Electronic supplementary material

The online version of this article (10.1186/s12915-018-0585-5) contains supplementary material, which is available to authorized users.

Control of human gene expression: High abundance of divergent transcription in genes containing both INR and BRE elements in the core promoter

Background

DNA sequence elements in the core promoter can play a central role in regulation of gene expression. Core elements (e.g. INR and TATA box) are located within ~50bp of the transcription start site and both upstream and downstream elements are known. Although all can affect the level of gene expression, their mechanism of action has yet to be fully defined. The studies described here are focused on two core promoter elements, INR and BRE, in the human genome. The locations of the two elements were determined in a large number of human promoters and the results were interpreted in terms of overall promoter function.

Results

A total of 13,406 promoters were collected from the reference version of the human genome and found to contain 62,891 INR sequences and 32,290 BRE. An INR sequence was found in the core region of 1231 (9.2%) promoters and a BRE in 2592 (19.3%); 121 promoters (0.9%) have both INR and BRE elements. Counts support the view that most human promoters lack an INR or BRE element in the core promoter. Further analysis was carried out with the aligned aggregate of promoters from each chromosome. The results showed distinct INR distributions in separate chromosome groups indicating a degree of chromosome specificity to the way core promoter elements are deployed in the genome. The rare promoters with both INR and BRE elements were found to be enriched among the genes with divergent transcription. Enrichment raises the possibility that core promoter elements can have a function in chromosome organization as well as in initiation of transcription.

Genome-wide analysis of the human malaria parasite Plasmodium falciparum transcription factor PfNF-YB shows interaction with a CCAAT motif

Little is known about transcription factor regulation during the Plasmodium falciparum intraerythrocytic cycle. In order to elucidate the role of the P. falciparum (Pf)NF-YB transcription factor we searched for target genes in the entire genome. PfNF-YB mRNA is highly expressed in late trophozoite and schizont stages relative to the ring stage. In order to determine the candidate genes bound by PfNF-YB a ChIP-on-chip assay was carried out and 297 genes were identified. Ninety nine percent of PfNF-YB binding was to putative promoter regions of protein coding genes of which only 16% comprise proteins of known function. Interestingly, our data reveal that PfNF-YB binding is not exclusively to a canonical CCAAT box motif. PfNF-YB binds to genes coding for proteins implicated in a range of different biological functions, such as replication protein A large subunit (DNA replication), hypoxanthine phosphoribosyltransferase (nucleic acid metabolism) and multidrug resistance protein 2 (intracellular transport).

HTLV-1 bZIP factor suppresses TDP1 expression through inhibition of NRF-1 in adult T-cell leukemia

Adult T-cell leukemia (ATL) is an aggressive T-cell malignancy caused by human T-cell leukemia virus type 1 (HTLV-1). We recently reported that abacavir, an anti-HIV-1 drug, potently and selectively kills ATL cells. This effect was attributed to the reduced expression of tyrosyl-DNA-phosphodiesterase 1 (TDP1), a DNA repair enzyme, in ATL cells. However, the molecular mechanism underlying the downregulation of TDP1 in ATL cells remains elusive. Here we identified the core promoter of the TDP1 gene, which contains a conserved nuclear respiratory factor 1 (NRF-1) binding site. Overexpression of NRF-1 increased TDP1-promoter activity, whereas the introduction of dominant-negative NRF-1 repressed such activity. Overexpression of NRF-1 also upregulated endogenous TDP-1 expression, while introduction of shNRF-1 suppressed TDP1 in Jurkat T cells, making them susceptible to abacavir. These results indicate that NRF-1 is a positive transcriptional regulator of TDP1-gene expression. Importantly, we revealed that HTLV-1 bZIP factor (HBZ) protein which is expressed in all ATL cases physically interacts with NRF-1 and inhibits the DNA-binding ability of NRF-1. Taken together, HBZ suppresses TDP1 expression by inhibiting NRF-1 function in ATL cells. The HBZ/NRF-1/TDP1 axis provides new therapeutic targets against ATL and might explain genomic instability leading to the pathogenesis of ATL.

Genome-wide analysis of cis-regulatory element structure and discovery of motif-driven gene co-expression networks in grapevine

Coordinated transcriptional and metabolic reprogramming ensures a plant’s continued growth and survival under adverse environmental conditions. Transcription factors (TFs) act to modulate gene expression through complex cis-regulatory element (CRE) interactions. Genome-wide analysis of known plant CREs was performed for all currently predicted protein-coding gene promoters in grapevine (Vitis vinifera L.). Many CREs such as abscisic acid (ABA)-responsive, drought-responsive, auxin-responsive, and evening elements, exhibit bona fide CRE properties such as strong position bias towards the transcription start site (TSS) and over-representation when compared with random promoters. Genes containing these CREs are enriched in a large repertoire of plant biological pathways. Large-scale transcriptome analyses also show that these CREs are highly implicated in grapevine development and stress response. Numerous CRE-driven modules in condition-specific gene co-expression networks (GCNs) were identified and many of these modules were highly enriched for plant biological functions. Several modules corroborate known roles of CREs in drought response, pathogen defense, cell wall metabolism, and fruit ripening, whereas others reveal novel functions in plants. Comparisons with Arabidopsis suggest a general conservation in promoter architecture, gene expression dynamics, and GCN structure across species. Systems analyses of CREs provide insights into the grapevine cis-regulatory code and establish a foundation for future genomic studies in grapevine.

Role of Y Box Protein-1 in cancer: As potential biomarker and novel therapeutic target

The Y-box binding protein (YB-1) is known to be a multifunctional transcription and translation factor during expression of several proteins. It is a vital oncoprotein that regulates cancer cell progression and proliferation. YB-1 is over-expressed in various human cancers such as breast cancer, colon cancer, lung cancer, gastric cancer, oesophageal cancer and glioblastoma. Nuclear expression of YB-1 is found to be associated with multidrug resistance and cancer cell progression. YB-1 is reported to regulate many cellular signalling pathways in different types of cancer proliferation. Knowledge about nuclear localization and nuclear level expression of YB-1 in different cancers has been correlated with prospective prognosis of cancer. This review discusses the prospects of YB-1 as a potential biomarker as well as therapeutic target in lieu of their role during cancer progression and multidrug resistance.

Krüppel-like factors compete for promoters and enhancers to fine-tune transcription

Krüppel-like factors (KLFs) are a family of 17 transcription factors characterized by a conserved DNA-binding domain of three zinc fingers and a variable N-terminal domain responsible for recruiting cofactors. KLFs have diverse functions in stem cell biology, embryo patterning, and tissue homoeostasis. KLF1 and related family members function as transcriptional activators via recruitment of co-activators such as EP300, whereas KLF3 and related members act as transcriptional repressors via recruitment of C-terminal Binding Proteins. KLF1 directly activates the Klf3 gene via an erythroid-specific promoter. Herein, we show KLF1 and KLF3 bind common as well as unique sites within the erythroid cell genome by ChIP-seq. We show KLF3 can displace KLF1 from key erythroid gene promoters and enhancers in vivo. Using 4sU RNA labelling and RNA-seq, we show this competition results in reciprocal transcriptional outputs for >50 important genes. Furthermore, Klf3-/- mice displayed exaggerated recovery from anemic stress and persistent cell cycling consistent with a role for KLF3 in dampening KLF1-driven proliferation. We suggest this study provides a paradigm for how KLFs work in incoherent feed-forward loops or networks to fine-tune transcription and thereby control diverse biological processes such as cell proliferation.

Towards a mechanistic understanding of core promoter recognition from cryo-EM studies of human TFIID

TFIID is a critical component of the eukaryotic transcription pre-initiation complex (PIC) required for the recruitment of RNA Pol II to the start site of protein-coding genes. Within the PIC, TFIID's role is to recognize and bind core promoter sequences and recruit the rest of the PIC components. Due to its size and its conformational complexity, TFIID poses a serious challenge for structural characterization. The small amounts of purified TFIID that can be obtained by present methods of purification from endogenous sources has limited structural studies to cryo-EM visualization, which requires very small amounts of sample. Previous cryo-EM studies have shed light on how the extreme conformational flexibility of TFIID is involved in core promoter DNA binding. Recent progress in cryo-EM methodology has facilitated a parallel progress in the study of human TFIID, leading to an improvement in resolution and the identification of the structural elements in the complex directly involved in DNA interaction. While many questions remain unanswered, the present structural knowledge of human TFIID suggests a mechanism for the sequential engagement with different core promoter sequences and how it could be influenced by regulatory factors.

Transcribing β-cell mitochondria in health and disease

Background

The recent genome-wide association studies (GWAS) of Type 2 Diabetes (T2D) have identified the pancreatic β-cell as the culprit in the pathogenesis of the disease. Mitochondrial metabolism plays a crucial role in the processes controlling release of insulin and β-cell mass. This notion implies that mechanisms controlling mitochondrial function have the potential to play a decisive pathogenetic role in T2D.

Scope of the review

This article reviews studies demonstrating that there is indeed mitochondrial dysfunction in islets in T2D, and that GWAS have identified a variant in the gene encoding transcription factor B1 mitochondrial (TFB1M), predisposing to T2D due to mitochondrial dysfunction and impaired insulin secretion. Mechanistic studies of the nature of this pathogenetic link, as well as of other mitochondrial transcription factors, are described.

Major conclusions

Based on this, it is argued that transcription and translation in mitochondria are critical processes determining mitochondrial function in β-cells in health and disease.

Identification of Arabidopsis genic and non-genic promoters by paired-end sequencing of TSS tags

Information about transcription start sites (TSSs) provides baseline data for the analysis of promoter architecture. In this paper we used paired- and single-end deep sequencing to analyze Arabidopsis TSS tags from several libraries prepared from roots, shoots, flowers and etiolated seedlings. The clustering of approximately 33 million mapped TSS tags led to the identification of 324 461 promoters that covered 79.7% (21 672/27 206) of protein-coding genes in the Arabidopsis genome. In addition we identified intragenic, antisense and orphan promoters that were not associated with any gene models. Of these, intragenic promoters exhibited unique characteristics regarding dinucleotide sequences at TSSs and core promoter element composition, suggesting that these promoters use different mechanisms of transcriptional initiation. An analysis of base composition with regard to promoter position revealed a low GC content throughout the promoter region and several local strand biases that were evident for TATA-type promoters, but not for Coreless-type promoters. Most observed strand biases coincided with strand biases of single nucleotide polymorphism rate. Our analysis also revealed that transcription of a gene is supported by an average of 2.7 genic promoters, among which one specific promoter, designated as a top promoter, substantially determines the expression level of the gene.

The garlic NF-YC gene, AsNF-YC8, positively regulates non-ionic hyperosmotic stress tolerance in tobacco

To investigate the relationship between nuclear factor Y (NF-Y) and stress tolerance in garlic, we cloned a NF-Y family gene AsNF-YC8 from garlic, which was largely upregulated at dehydrate stage. Expression pattern analyses in garlic revealed that AsNF-YC8 is induced through abscisic acid (ABA) and abiotic stresses, such as NaCl and PEG. Compared with wild-type plants, the overexpressing-AsNF-YC8 transgenic tobacco plants showed higher seed germination rates, longer root length and better plant growth under salt and drought stresses. Under drought stress, the transgenic plants maintained higher relative water content (RWC), net photosynthesis, lower levels of malondialdehyde (MDA), and less ion leakage (IL) than wild-type control plants. These results indicate the high tolerance of the transgenic plants to drought stress compared to the WT. The transgenic tobacco lines accumulated less reactive oxygen species (ROS) and exhibited higher antioxidative enzyme activities compared with wild-type (WT) plants under drought stress, which suggested that the overexpression of AsNF-YC8 improves the antioxidant defense system by regulating the activities of these antioxidant enzymes, which in turn protect transgenic lines against drought stress. These results suggest that AsNF-YC8 plays an important role in tolerance to drought and salt stresses.

High constitutive activity of a broad panel of housekeeping and tissue-specific cis-regulatory elements depends on a subset of ETS proteins

Here, Curina et al. found that a subgroup of ETS family transcription factors (the ELF proteins) shows a strong bias in its genomic distribution by binding in very close proximity (<60 nt) to the transcription start sites of housekeeping genes. They show that a limited number of highly active transcription factors can equip cis-regulatory elements with disparate functional roles and cell type specificity with the ability to efficiently promote transcription.

Enhancers and promoters that control the transcriptional output of terminally differentiated cells include cell type-specific and broadly active housekeeping elements. Whether the high constitutive activity of these two groups of cis-regulatory elements relies on entirely distinct or instead also on shared regulators is unknown. By dissecting the cis-regulatory repertoire of macrophages, we found that the ELF subfamily of ETS proteins selectively bound within 60 base pairs (bp) from the transcription start sites of highly active housekeeping genes. ELFs also bound constitutively active, but not poised, macrophage-specific enhancers and promoters. The role of ELFs in promoting high-level constitutive transcription was suggested by multiple evidence: ELF sites enabled robust transcriptional activation by endogenous and minimal synthetic promoters, ELF recruitment was stabilized by the transcriptional machinery, and ELF proteins mediated recruitment of transcriptional and chromatin regulators to core promoters. These data suggest that the co-optation of a limited number of highly active transcription factors represents a broadly adopted strategy to equip both cell type-specific and housekeeping cis-regulatory elements with the ability to efficiently promote transcription.

DNA structural features of eukaryotic TATA-containing and TATA-less promoters

Eukaryotic genes can be broadly classified as TATA‐containing and TATA‐less based on the presence of TATA box in their promoters. Experiments on both classes of genes have revealed a disparity in the regulation of gene expression and cellular functions between the two classes. In this study, we report characteristic differences in promoter sequences and associated structural properties of the two categories of genes in six different eukaryotes. We have analyzed three structural features, DNA duplex stability, bendability, and curvature along with the distribution of A‐tracts, G‐quadruplex motifs, and CpG islands. The structural feature analyses reveal that while the two classes of gene promoters are distinctly different from each other, the properties are also distinguishable across the six organisms.

Epigenomics of human CD8 T cell differentiation and aging

The efficacy of the adaptive immune response declines dramatically with age, but the cell-intrinsic mechanisms driving immune aging in humans remain poorly understood. Immune aging is characterized by a loss of self-renewing naïve cells and the accumulation of differentiated but dysfunctional cells within the CD8 T cell compartment. Using ATAC-seq, we inferred the transcription factor binding activities correlated with naive and central and effector memory CD8 T cell states in young adults. Integrating our results with RNA-seq, we identified transcription networks associated with CD8 T cell differentiation, with prominent roles implicated for BATF, ETS1, Eomes, and Sp1. Extending our analysis to aged humans, we found that the differences between the memory and naive subsets were largely preserved across age, but that naive and central memory cells from older individuals exhibited a shift toward more differentiated patterns of chromatin openness. Additionally, aged naive cells displayed a loss in chromatin accessibility at gene promoters, largely associated with a decrease in NRF1 binding. This shift was implicated in a marked drop-off in the ability of the aged naive cells to transcribe respiratory chain genes, which may explain the reduced capacity of oxidative phosphorylation in older naïve cells. Our findings identify BATF- and NRF1-driven gene regulation as potential targets for delaying CD8 T cell aging and restoring function.

A pilot study on the prevalence of DNA palindromes in breast cancer genomes

Background

DNA palindromes are a unique pattern of repeat sequences that are present in the human genome. It consists of a sequence of nucleotides in which the second half is the complement of the first half but appearing in reverse order. These palindromic sequences may have a significant role in DNA replication, transcription and gene regulation processes. They occur frequently in human cancers by clustering at specific locations of the genome that undergo gene amplification and tumorigenesis. Moreover, some studies showed that palindromes are clustered in amplified regions of breast cancer genomes especially in chromosomes (chr) 8 and 11. With the large number of personal genomes and cancer genomes becoming available, it is now possible to study their association to diseases using computational methods. Here, we conducted a pilot study on chromosomes 8 and 11 of cancer genomes to identify computationally the differentially occurring palindromes.

Methods

We processed 69 breast cancer genomes from The Cancer Genome Atlas including serum-normal and tumor genomes, and 1000 Genomes to serve as control group. The Biological Language Modelling Toolkit (BLMT) computes palindromes in whole genomes. We developed a computational pipeline integrating BLMT to compute and compare prevalence of palindromes in personal genomes.

Results

We carried out a pilot study on chr 8 and chr 11 taking into account single nucleotide polymorphisms, insertions and deletions. Of all the palindromes that showed any variation in cancer genomes, 38% of what were near breast cancer genes happened to be the most differentiated palindromes in tumor (i.e. they ranked among the top 25% by our heuristic measure).

Conclusions

These results will shed light on the prevalence of palindromes in oncogenes and the mutations that are present in the palindromic regions that could contribute to genomic rearrangements, and breast cancer progression.

Most of the tight positional conservation of transcription factor binding sites near the transcription start site reflects their co-localization within regulatory modules

Background

Transcription factors (TFs) form complexes that bind regulatory modules (RMs) within DNA, to control specific sets of genes. Some transcription factor binding sites (TFBSs) near the transcription start site (TSS) display tight positional preferences relative to the TSS. Furthermore, near the TSS, RMs can co-localize TFBSs with each other and the TSS. The proportion of TFBS positional preferences due to TFBS co-localization within RMs is unknown, however. ChIP experiments confirm co-localization of some TFBSs genome-wide, including near the TSS, but they typically examine only a few TFs at a time, using non-physiological conditions that can vary from lab to lab. In contrast, sequence analysis can examine many TFs uniformly and methodically, broadly surveying the co-localization of TFBSs with tight positional preferences relative to the TSS.

Results

Our statistics found 43 significant sets of human motifs in the JASPAR TF Database with positional preferences relative to the TSS, with 38 preferences tight (±5 bp). Each set of motifs corresponded to a gene group of 135 to 3304 genes, with 42/43 (98%) gene groups independently validated by DAVID, a gene ontology database, with FDR < 0.05. Motifs corresponding to two TFBSs in a RM should co-occur more than by chance alone, enriching the intersection of the gene groups corresponding to the two TFs. Thus, a gene-group intersection systematically enriched beyond chance alone provides evidence that the two TFs participate in an RM. Of the 903 = 43*42/2 intersections of the 43 significant gene groups, we found 768/903 (85%) pairs of gene groups with significantly enriched intersections, with 564/768 (73%) intersections independently validated by DAVID with FDR < 0.05. A user-friendly web site at http://go.usa.gov/3kjsH permits biologists to explore the interaction network of our TFBSs to identify candidate subunit RMs.

Conclusions

Gene duplication and convergent evolution within a genome provide obvious biological mechanisms for replicating an RM near the TSS that binds a particular TF subunit. Of all intersections of our 43 significant gene groups, 85% were significantly enriched, with 73% of the significant enrichments independently validated by gene ontology. The co-localization of TFBSs within RMs therefore likely explains much of the tight TFBS positional preferences near the TSS.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-1354-5) contains supplementary material, which is available to authorized users.

The multifaceted roles of NUCLEAR FACTOR-Y in Arabidopsis thaliana development and stress responses

NUCLEAR FACTOR-Y (NF-Y) is a heterotrimeric transcription factor (TF) consisting of evolutionarily distinct NF-YA, NF-YB and NF-YC subunits. The functional NF-Y heterotrimer binds to CCAAT elements in eukaryotic gene promoters and influences their expression. The genome of the model organism Arabidopsis thaliana encodes 10 distinct NF-YA, NF-YB, and NF-YC proteins, allowing for enormous combinatorial and functional diversity. Two decades of research have elucidated the importance of NF-Ys in plant growth, development and stress responses; however, the molecular mechanisms of action remain largely unexplored. Intriguingly, recent evidence suggests that NF-Ys are frequently associated with other groups of TFs, expanding the potential NF-Y combinatorial complexity. Further, information regarding the regulation of individual NF-Y subunits at the transcriptional and post-transcriptional level is beginning to emerge. In this review, we will identify developing trends within the NF-Y field and discuss recent progress towards a better understanding of NF-Y function, molecular action, and regulation in the context of Arabidopsis. This article is part of a Special Issue entitled: Nuclear Factor Y in Development and Disease, edited by Prof. Roberto Mantovani.

Genome wide features, distribution and correlations of NF-Y binding sites

NF-Y is a trimeric transcription factor that binds on DNA the CCAAT-box motif. In this article we reviewed and complemented with additional bioinformatic analysis existing data on genome-wide NF-Y binding characterization in human, reaching the following main conclusions: (1) about half of NF-Y binding sites are located at promoters, about 60-80 base pairs from transcription start sites; NF-Y binding to distal genomic regions takes place at inactive chromatin loci and/or DNA repetitive elements more often than active enhancers; (2) on almost half of its binding sites, regardless of their genomic localization (promoters or distal regions), NF-Y finds on DNA more than one CCAAT-box, and most of those multiple CCAAT binding loci present precise spacing and organization of the elements composing them; (3) there exists a well defined class of transcription factors that show genome-wide co-localization with NF-Y. Some of them lack their canonical binding site in binding regions overlapping with NF-Y, hence hinting at NF-Y mediated recruitment, while others show a precise positioning on DNA of their binding sites with respect to the CCAAT box bound by NF-Y. This article is part of a Special Issue entitled: Nuclear Factor Y in Development and Disease, edited by Prof. Roberto Mantovani.

Promoter Architecture and Sex-Specific Gene Expression in Daphnia pulex

Large-scale transcription start site (TSS) profiling produces a high-resolution, quantitative picture of transcription initiation and core promoter locations within a genome. However, application of TSS profiling to date has largely been restricted to a small set of prominent model systems. We sought to characterize the cis-regulatory landscape of the water flea Daphnia pulex, an emerging model arthropod that reproduces both asexually (via parthenogenesis) and sexually (via meiosis). We performed Cap Analysis of Gene Expression (CAGE) with RNA isolated from D. pulex within three developmental states: sexual females, asexual females, and males. Identified TSSs were utilized to generate a "Daphnia Promoter Atlas," i.e., a catalog of active promoters across the surveyed states. Analysis of the distribution of promoters revealed evidence for widespread alternative promoter usage in D. pulex, in addition to a prominent fraction of compactly-arranged promoters in divergent orientations. We carried out de novo motif discovery using CAGE-defined TSSs and identified eight candidate core promoter motifs; this collection includes canonical promoter elements (e.g., TATA and Initiator) in addition to others lacking obvious orthologs. A comparison of promoter activities found evidence for considerable state-specific differential gene expression between states. Our work represents the first global definition of transcription initiation and promoter architecture in crustaceans. The Daphnia Promoter Atlas presented here provides a valuable resource for comparative study of cis-regulatory regions in metazoans, as well as for investigations into the circuitries that underpin meiosis and parthenogenesis.

TRIB3 increases cell resistance to arsenite toxicity by limiting the expression of the glutathione-degrading enzyme CHAC1

Arsenic, a metalloid with cytotoxic and carcinogenic effects related to the disruption of glutathione homeostasis, induces the expression of ATF4, a central transcription factor in the cellular stress response. However, the interplay between factors downstream of ATF4 is incompletely understood. In this article, we investigate the role of Tribbles homolog 3 (TRIB3), a regulatory member of the ATF4 pathway, in determining cell sensitivity to arsenite. Our results show that arsenite potently upregulates Trib3 mRNA and protein in an ATF4-dependent manner in mouse embryonic fibroblasts. Trib3-deficient cells display increased susceptibility to arsenite-induced cell death, which is rescued by re-expressing TRIB3. In cells lacking TRIB3, arsenite stress leads to markedly elevated mRNA and protein levels of Chac1, a gene that encodes a glutathione-degrading enzyme and is not previously known to be repressed by TRIB3. Analysis of the Chac1 promoter identified two regulatory elements that additively mediate the induction of Chac1 by arsenite and ATF4, as well as the robust suppression of Chac1 by TRIB3. Crucially, Chac1 silencing enhances glutathione levels and eliminates the increased susceptibility of Trib3-deficient cells to arsenite stress. Moreover, Trib3-deficient cells demonstrate an increased rate of glutathione consumption, which is abolished by Chac1 knockdown. Taken together, these data indicate that excessive Chac1 expression is detrimental to arsenite-treated cell survival and that TRIB3 is critical for restraining the pro-death potential of Chac1 during arsenite stress, representing a novel mechanism of cell viability regulation that occurs within the ATF4 pathway.

A high definition look at the NF-Y regulome reveals genome-wide associations with selected transcription factors

NF-Y is a trimeric transcription factor (TF), binding the CCAAT box element, for which several results suggest a pioneering role in activation of transcription. In this work, we integrated 380 ENCODE ChIP-Seq experiments for 154 TFs and cofactors with sequence analysis, protein–protein interactions and RNA profiling data, in order to identify genome-wide regulatory modules resulting from the co-association of NF-Y with other TFs. We identified three main degrees of co-association with NF-Y for sequence-specific TFs. In the most relevant one, we found TFs having a significant overlap with NF-Y in their DNA binding loci, some with a precise spacing of binding sites with respect to the CCAAT box, others (FOS, Sp1/2, RFX5, IRF3, PBX3) mostly lacking their canonical binding site and bound to arrays of well spaced CCAAT boxes. As expected, NF-Y binding also correlates with RNA Pol II General TFs and with subunits of complexes involved in the control of H3K4 methylations. Co-association patterns are confirmed by protein–protein interactions, and correspond to specific functional categorizations and expression level changes of target genes following NF-Y inactivation. These data define genome-wide rules for the organization of NF-Y-centered regulatory modules, supporting a model of distinct categorization and synergy with well defined sets of TFs.

C/EBPβ (CEBPB) protein binding to the C/EBP|CRE DNA 8-mer TTGC|GTCA is inhibited by 5hmC and enhanced by 5mC, 5fC, and 5caC in the CG dinucleotide

During mammalian development, some methylated cytosines (5mC) in CG dinucleotides are iteratively oxidized by TET dioxygenases to 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC), and 5-carboxylcytosine (5caC). The effect of these cytosine oxidative products on the sequence-specific DNA binding of transcription factors is being actively investigated. Here, we used the electrophoretic mobility shift assay (EMSA) to examine C/EBPα and C/EBPβ homodimers binding to all 25 chemical forms of a CG dinucleotide for two DNA sequences: the canonical C/EBP 8-mer TTGC|GCAA and the chimeric C/EBP|CRE 8-mer TTGC|GTCA. 5hmC in the CG dinucleotide in the C/EBP|CRE motif 8-mer TGAC|GCAA inhibits binding of C/EBPβ but not C/EBPα. Binding was increased by 5mC, 5fC and 5caC. Circular dichroism monitored thermal denaturations for C/EBPβ bound to the C/EBP|CRE motif confirmed the EMSA. The structural differences between C/EBPα and C/EBPβ that may account for the difference in binding 5hmC in the 8-mer TGAC|GCAA are explored.

Foxtail Millet NF-Y Families: Genome-Wide Survey and Evolution Analyses Identified Two Functional Genes Important in Abiotic Stresses

It was reported that Nuclear Factor Y (NF-Y) genes were involved in abiotic stress in plants. Foxtail millet (Setaria italica), an elite stress tolerant crop, provided an impetus for the investigation of the NF-Y families in abiotic responses. In the present study, a total of 39 NF-Y genes were identified in foxtail millet. Synteny analyses suggested that foxtail millet NF-Y genes had experienced rapid expansion and strong purifying selection during the process of plant evolution. De novo transcriptome assembly of foxtail millet revealed 11 drought up-regulated NF-Y genes. SiNF-YA1 and SiNF-YB8 were highly activated in leaves and/or roots by drought and salt stresses. Abscisic acid (ABA) and H2O2 played positive roles in the induction of SiNF-YA1 and SiNF-YB8 under stress treatments. Transient luciferase (LUC) expression assays revealed that SiNF-YA1 and SiNF-YB8 could activate the LUC gene driven by the tobacco (Nicotiana tobacam) NtERD10, NtLEA5, NtCAT, NtSOD, or NtPOD promoter under normal or stress conditions. Overexpression of SiNF-YA1 enhanced drought and salt tolerance by activating stress-related genes NtERD10 and NtCAT1 and by maintaining relatively stable relative water content (RWC) and contents of chlorophyll, superoxide dismutase (SOD), peroxidase (POD), catalase (CAT) and malondialdehyde (MDA) in transgenic lines under stresses. SiNF-YB8 regulated expression of NtSOD, NtPOD, NtLEA5, and NtERD10 and conferred relatively high RWC and chlorophyll contents and low MDA content, resulting in drought and osmotic tolerance in transgenic lines under stresses. Therefore, SiNF-YA1 and SiNF-YB8 could activate stress-related genes and improve physiological traits, resulting in tolerance to abiotic stresses in plants. All these results will facilitate functional characterization of foxtail millet NF-Ys in future studies.

Nucleosomes are enriched at the boundaries of hypomethylated regions (HMRs) in mouse dermal fibroblasts and keratinocytes

Background

The interplay between epigenetic modifications and chromatin structure are integral to our understanding of genome function. Methylation of cytosine (5mC) at CG dinucleotides, traditionally associated with transcriptional repression, is the most highly studied chemical modification of DNA, occurring at over 70% of all CG dinucleotides in the genome. Hypomethylated regions (HMRs) often occur in CG islands (CGIs), however, they also occur outside of CGIs and function as cell-type specific enhancers. During the process of differentiation, reorganization of chromatin and nucleosome arrangement at regulatory regions is thought to occur in order for the establishment of cell-type specific transcriptional programs. However, the specifics regarding the organization of nucleosomes at HMRs and the potential mechanisms regulating nucleosome occupancy in these regions are unknown. Here, we have investigated nucleosome organization around hypomethylated regions (HMRs) identified in two mouse primary cells.

Results

Microccocal nuclease (MNase) digested mononucleosomes from primary cultures of new-born female mouse dermal fibroblasts and keratinocytes were mapped and compared to the HMRs obtained from single base-pair resolution methylomes. In both cell types, we find that nucleosomes are enriched at HMR boundaries. In contrast to the nucleosomes found at boundaries of HMRs in CGIs, HMRs outside of CGIs are calculated to be preferentially bound by nucleosomes, with phased nucleosomes propagating into the methylated region. Nucleosomes are enriched at the tissue-specific HMRs (TS-HMR) boundaries in both cell types suggesting that nucleosome organization surrounding HMR boundaries is independent of methylation status. In addition, we find potential transcription factor (TF) binding sites (E-box motifs) enriched in non-CGI TS-HMR boundaries.

Conclusions

Our results show that intrinsic nucleosome occupancy score (INOS) positively correlate with the nucleosome organization surrounding non-CGI TS-HMRs, suggesting that DNA sequence plays a role in the establishment of HMRs in the genome. Since nucleosomes impact all processes involving the genome, our results provide a link between epigenetic modifications, chromatin structure, and regulatory function.

Electronic supplementary material

The online version of this article (doi:10.1186/1756-8935-7-34) contains supplementary material, which is available to authorized users.

Determination and inference of eukaryotic transcription factor sequence specificity

Transcription factor (TF) DNA sequence preferences direct their regulatory activity, but are currently known for only ∼1% of eukaryotic TFs. Broadly sampling DNA-binding domain (DBD) types from multiple eukaryotic clades, we determined DNA sequence preferences for >1,000 TFs encompassing 54 different DBD classes from 131 diverse eukaryotes. We find that closely related DBDs almost always have very similar DNA sequence preferences, enabling inference of motifs for ∼34% of the ∼170,000 known or predicted eukaryotic TFs. Sequences matching both measured and inferred motifs are enriched in chromatin immunoprecipitation sequencing (ChIP-seq) peaks and upstream of transcription start sites in diverse eukaryotic lineages. SNPs defining expression quantitative trait loci in Arabidopsis promoters are also enriched for predicted TF binding sites. Importantly, our motif "library" can be used to identify specific TFs whose binding may be altered by human disease risk alleles. These data present a powerful resource for mapping transcriptional networks across eukaryotes.

GA binding protein augments autophagy via transcriptional activation of BECN1-PIK3C3 complex genes

Macroautophagy is a vesicular catabolic trafficking pathway that is thought to protect cells from diverse stressors and to promote longevity. Recent studies have revealed that transcription factors play important roles in the regulation of autophagy. In this study, we have identified GA binding protein (GABP) as a transcriptional regulator of the combinatorial expression of BECN1-PIK3C3 complex genes involved in autophagosome initiation. We performed bioinformatics analyses that demonstrated highly conserved putative GABP sites in genes that encode BECN1/Beclin 1, several BECN1 interacting proteins, and downstream autophagy proteins including the ATG12-ATG5-ATG16L1 complex. We demonstrate that GABP binds to the promoter regions of BECN1-PIK3C3 complex genes and activates their transcriptional activities. Knockdown of GABP reduced BECN1-PIK3C3 complex transcripts, BECN1-PIK3C3 complex protein levels and autophagy in cultured cells. Conversely, overexpression of GABP increased autophagy. Nutrient starvation increased GABP-dependent transcriptional activity of BECN1-PIK3C3 complex gene promoters and increased the recruitment of GABP to the BECN1 promoter. Our data reveal a novel function of GABP in the regulation of autophagy via transcriptional activation of the BECN1-PIK3C3 complex.

Leveraging cross-species transcription factor binding site patterns: from diabetes risk loci to disease mechanisms

Genome-wide association studies have revealed numerous risk loci associated with diverse diseases. However, identification of disease-causing variants within association loci remains a major challenge. Divergence in gene expression due to cis-regulatory variants in noncoding regions is central to disease susceptibility. We show that integrative computational analysis of phylogenetic conservation with a complexity assessment of co-occurring transcription factor binding sites (TFBS) can identify cis-regulatory variants and elucidate their mechanistic role in disease. Analysis of established type 2 diabetes risk loci revealed a striking clustering of distinct homeobox TFBS. We identified the PRRX1 homeobox factor as a repressor of PPARG2 expression in adipose cells and demonstrate its adverse effect on lipid metabolism and systemic insulin sensitivity, dependent on the rs4684847 risk allele that triggers PRRX1 binding. Thus, cross-species conservation analysis at the level of co-occurring TFBS provides a valuable contribution to the translation of genetic association signals to disease-related molecular mechanisms.