NF-Y and SP transcription factors — New insights in a long-standing liaison

Liver kinase B1 (LKB1) regulates the epigenetic landscape of mouse pancreatic beta cells

Liver kinase B1 (LKB1/STK11) is an important regulator of pancreatic β-cell identity and function. Elimination of Lkb1 from the β-cell results in improved glucose-stimulated insulin secretion and is accompanied by profound changes in gene expression, including the upregulation of several neuronal genes. The mechanisms through which LKB1 controls gene expression are, at present, poorly understood. Here, we explore the impact of β cell- selective deletion of Lkb1 on chromatin accessibility in mouse pancreatic islets. To characterize the role of LKB1 in the regulation of gene expression at the transcriptional level, we combine these data with a map of islet active transcription start sites and histone marks. We demonstrate that LKB1 elimination from β-cells results in widespread changes in chromatin accessibility, correlating with changes in transcript levels. Changes occurred in hundreds of promoter and enhancer regions, many of which were close to neuronal genes. We reveal that dysregulated enhancers are enriched in binding motifs for transcription factors important for β-cell identity, such as FOXA, MAFA or RFX6 and we identify microRNAs (miRNAs) that are regulated by LKB1 at the transcriptional level. Overall, our study provides important new insights into the epigenetic mechanisms by which LKB1 regulates β-cell identity and function.

COL1A1 proximal promoter topology regulates its transcriptional response to transforming growth factor β

OBJECTIVE

The COL1A1 proximal promoter contains two GC-rich regions and two inverted CCAAT boxes. The transcription factors Sp1 and CBF bind to the GC sequence at -122 to -115 bp and the inverted CCAAT box at -101 to -96 bp, respectively, and stimulate COL1A1 transcriptional activity.

METHODS

To further define the regulatory mechanisms controlling COL1A1 expression by Sp1 and CBF, we introduced 2, 4, 6, or 8 thymidine nucleotides (T-tracts) at position -111 bp of the COL1A1 gene promoter to increase the physical distance between these two binding sites and examined in vitro the transcriptional activities of the resulting constructs and their response to TGF-β1.`.

RESULTS

Insertion of 2 or 4 nucleotides decreased COL1A1 promoter activity by up to 70%. Furthermore, the expected increase in COL1A1 transcription in response to TGF-β1 was abolished. Computer modeling of the modified DNA structure indicated that increasing the physical distance between the Sp1 and CBF binding sites introduces a rotational change in the DNA topology that disrupts the alignment of Sp1 and CBF binding sites and likely alters protein-protein interactions among these transcription factors or their associated co-activators.

CONCLUSION

The topology of the COL1A1 proximal promoter is crucial in determining the transcriptional activity of the gene and its response to the stimulatory effects of TGF-β1.

Identification of putative promoter elements for epsilon glutathione s-transferases genes associated with resistance to DDT in the malaria vector mosquito anopheles arabiensis

The purpose of this study was to identify the putative regulatory elements in the promoter region of An. arabiensis strains which differed in susceptibility to DDT and compare with those identified in its sibling An. gambaie. Basal expression level of Epsilon class GSTs (Glutathione S - transferases) GSTe1 gene was 0.512 - 0.658 (95% CI) and GSTe2 0.672 - 1.204 (95% CI) in adults of DDT resistant KGB compared to 0.031 - 0.04 (95% CI) and 0.148 - 0.199 (95% CI) respectively in susceptible MAT strains of An. arabiensis. Induced mean expression of GSTe2 in larvae exposed to DDT for one hour was 0.901 - 1.172 (95% CI) in KGB and 0.475 - 0.724 (95% CI) in MAT strain. In present work, strain specific primers were used to amplify and sequenced the promoter regions of GSTe1 and GSTe2 in the KGB, MAT and field specimens. Computational analysis revealed presence of classical arthropod initiator sequence TCAGT and putative core promoter elements, GC, CAAT, TATA boxes. A typical TATA box was identified at 35 bp upstream Transcription Start Site (TSS) in GSTe1 but was absent in GSTe2. Several binding sites for regulatory elements downstream and multiple polymorphic sites were identified between strains. The role of these regulatory elements in transcription of these genes has not been determined. However, on comparison the 2 bp adenosine indel (insertion/deletion) which was essential in driving the promoter activity in An. gambiae was identified only DDT resistant KGB strain.

The transcription factor NF-YA is crucial for neural progenitor maintenance during brain development

In contrast to stage-specific transcription factors, the role of ubiquitous transcription factors in neuronal development remains a matter of scrutiny. Here, we demonstrated that a ubiquitous factor NF-Y is essential for neural progenitor maintenance during brain morphogenesis. Deletion of the NF-YA subunit in neural progenitors by using nestin-cre transgene in mice resulted in significant abnormalities in brain morphology, including a thinner cerebral cortex and loss of striatum during embryogenesis. Detailed analyses revealed a progressive decline in multiple neural progenitors in the cerebral cortex and ganglionic eminences, accompanied by induced apoptotic cell death and reduced cell proliferation. In neural progenitors, the NF-YA short isoform lacking exon 3 is dominant and co-expressed with cell cycle genes. ChIP-seq analysis from the cortex during early corticogenesis revealed preferential binding of NF-Y to the cell cycle genes, some of which were confirmed to be downregulated following NF-YA deletion. Notably, the NF-YA short isoform disappears and is replaced by its long isoform during neuronal differentiation. Forced expression of the NF-YA long isoform in neural progenitors resulted in a significant decline in neuronal count, possibly due to the suppression of cell proliferation. Collectively, we elucidated a critical role of the NF-YA short isoform in maintaining neural progenitors, possibly by regulating cell proliferation and apoptosis. Moreover, we identified an isoform switch in NF-YA within the neuronal lineage in vivo, which may explain the stage-specific role of NF-Y during neuronal development.

Transcription and chromatin regulation by TAF4b during cellular quiescence of developing prospermatogonia

Prospermatogonia (ProSpg) link the embryonic development of male primordial germ cells to the healthy establishment of postnatal spermatogonia and spermatogonial stem cells. While these spermatogenic precursor cells undergo the characteristic transitions of cycling and quiescence, the transcriptional events underlying these developmental hallmarks remain unknown. Here, we investigated the expression and function of TBP-associated factor 4b (Taf4b) in the timely development of quiescent mouse ProSpg using an integration of gene expression profiling and chromatin mapping. We find that Taf4b mRNA expression is elevated during the transition of mitotic-to-quiescent ProSpg and Taf4b-deficient ProSpg are delayed in their entry into quiescence. Gene ontology, protein network analysis, and chromatin mapping demonstrate that TAF4b is a direct and indirect regulator of chromatin and cell cycle-related gene expression programs during ProSpg quiescence. Further validation of these cell cycle mRNA changes due to the loss of TAF4b was accomplished via immunostaining for proliferating cell nuclear antigen (PCNA). Together, these data indicate that TAF4b is a key transcriptional regulator of the chromatin and quiescent state of the developing mammalian spermatogenic precursor lineage.

NF-YA isoforms with alternative splicing of exon-5 in Aves

NF-YA, the regulatory subunit of the trimeric CCAAT-binding transcription factor NF-Y, is present in vertebrates in two major alternative spliced isoforms: NF-YAl and NF-YAs, differing for the presence of exon-3. NF-YAx, a third isoform without exon-3/-5, was reported only in human neuronal cells and tumors. These events affect the Trans-Activation Domain. We provide here evidence for the expression of NF-YAx and for the existence of a new isoform, NF-YAg, skipping only exon-5. These isoforms are abundant in Aves, but not in reptiles, and are the prevalent transcripts in the initial phases of embryo development in chicken. Finally, we analyzed NF-YAg and NF-YAx amino acid sequence using AlphaFold: absence of exon-5 denotes a global reduction of β-stranded elements, while removal of the disordered exon-3 sequence has limited effects on TAD architecture. These data identify an expanded program of NF-YA isoforms within the TAD in Aves, implying a role during early development.

Sp3 is essential for normal lung morphogenesis and cell cycle progression during mouse embryonic development

Members of the Sp family of transcription factors regulate gene expression via binding GC boxes within promoter regions. Unlike Sp1, which stimulates transcription, the closely related Sp3 can either repress or activate gene expression and is required for perinatal survival in mice. Here, we use RNA-seq and cellular phenotyping to show how Sp3 regulates murine fetal cell differentiation and proliferation. Homozygous Sp3-/- mice were smaller than wild-type and Sp+/- littermates, died soon after birth and had abnormal lung morphogenesis. RNA-seq of Sp3-/- fetal lung mesenchymal cells identified alterations in extracellular matrix production, developmental signaling pathways and myofibroblast/lipofibroblast differentiation. The lungs of Sp3-/- mice contained multiple structural defects, with abnormal endothelial cell morphology, lack of elastic fiber formation, and accumulation of lipid droplets within mesenchymal lipofibroblasts. Sp3-/- cells and mice also displayed cell cycle arrest, with accumulation in G0/G1 and reduced expression of numerous cell cycle regulators including Ccne1. These data detail the global impact of Sp3 on in vivo mouse gene expression and development.

Crucial Abiotic Stress Regulatory Network of NF-Y Transcription Factor in Plants

Nuclear Factor-Y (NF-Y), composed of three subunits NF-YA, NF-YB and NF-YC, exists in most of the eukaryotes and is relatively conservative in evolution. As compared to animals and fungi, the number of NF-Y subunits has significantly expanded in higher plants. The NF-Y complex regulates the expression of target genes by directly binding the promoter CCAAT box or by physical interaction and mediating the binding of a transcriptional activator or inhibitor. NF-Y plays an important role at various stages of plant growth and development, especially in response to stress, which attracted many researchers to explore. Herein, we have reviewed the structural characteristics and mechanism of function of NF-Y subunits, summarized the latest research on NF-Y involved in the response to abiotic stresses, including drought, salt, nutrient and temperature, and elaborated the critical role of NF-Y in these different abiotic stresses. Based on the summary above, we have prospected the potential research on NF-Y in response to plant abiotic stresses and discussed the difficulties that may be faced in order to provide a reference for the in-depth analysis of the function of NF-Y transcription factors and an in-depth study of plant responses to abiotic stress.

Loss of the Ash2l subunit of histone H3K4 methyltransferase complexes reduces chromatin accessibility at promoters

Changes in gene expression programs are intimately linked to cell fate decisions. Post-translational modifications of core histones contribute to control gene expression. Methylation of lysine 4 of histone H3 (H3K4) correlates with active promoters and gene transcription. This modification is catalyzed by KMT2 methyltransferases, which require interaction with 4 core subunits, WDR5, RBBP5, ASH2L and DPY30, for catalytic activity. Ash2l is necessary for organismal development and for tissue homeostasis. In mouse embryo fibroblasts (MEFs), Ash2l loss results in gene repression, provoking a senescence phenotype. We now find that upon knockout of Ash2l both H3K4 mono- and tri-methylation (H3K4me1 and me3, respectively) were deregulated. In particular, loss of H3K4me3 at promoters correlated with gene repression, especially at CpG island promoters. Ash2l loss resulted in increased loading of histone H3 and reduced chromatin accessibility at promoters, accompanied by an increase of repressing and a decrease of activating histone marks. Moreover, we observed altered binding of CTCF upon Ash2l loss. Lost and gained binding was noticed at promoter-associated and intergenic sites, respectively. Thus, Ash2l loss and reduction of H3K4me3 correlate with altered chromatin accessibility and transcription factor binding. These findings contribute to a more detailed understanding of mechanistic consequences of H3K4me3 loss and associated repression of gene transcription and thus of the observed cellular consequences.

SP1 and NFY Regulate the Expression of PNPT1, a Gene Encoding a Mitochondrial Protein Involved in Cancer

The Polyribonucleotide nucleotidyltransferase 1 gene (PNPT1) encodes polynucleotide phosphorylase (PNPase), a 3′-5′ exoribonuclease involved in mitochondrial RNA degradation and surveillance and RNA import into the mitochondrion. Here, we have characterized the PNPT1 promoter by in silico analysis, luciferase reporter assays, electrophoretic mobility shift assays (EMSA), chromatin immunoprecipitation (ChIP), siRNA-based mRNA silencing and RT-qPCR. We show that the Specificity protein 1 (SP1) transcription factor and Nuclear transcription factor Y (NFY) bind the PNPT1 promoter, and have a relevant role regulating the promoter activity, PNPT1 expression, and mitochondrial activity. We also found in Kaplan–Meier survival curves that a high expression of either PNPase, SP1 or NFY subunit A (NFYA) is associated with a poor prognosis in liver cancer. In summary, our results show the relevance of SP1 and NFY in PNPT1 expression, and point to SP1/NFY and PNPase as possible targets in anti-cancer therapy.

WhichTF is functionally important in your open chromatin data?

We present WhichTF, a computational method to identify functionally important transcription factors (TFs) from chromatin accessibility measurements. To rank TFs, WhichTF applies an ontology-guided functional approach to compute novel enrichment by integrating accessibility measurements, high-confidence pre-computed conservation-aware TF binding sites, and putative gene-regulatory models. Comparison with prior sheer abundance-based methods reveals the unique ability of WhichTF to identify context-specific TFs with functional relevance, including NF-κB family members in lymphocytes and GATA factors in cardiac cells. To distinguish the transcriptional regulatory landscape in closely related samples, we apply differential analysis and demonstrate its utility in lymphocyte, mesoderm developmental, and disease cells. We find suggestive, under-characterized TFs, such as RUNX3 in mesoderm development and GLI1 in systemic lupus erythematosus. We also find TFs known for stress response, suggesting routine experimental caveats that warrant careful consideration. WhichTF yields biological insight into known and novel molecular mechanisms of TF-mediated transcriptional regulation in diverse contexts, including human and mouse cell types, cell fate trajectories, and disease-associated cells.

Transcription factors (TFs), a class of DNA binding proteins, regulate tissue- and cell-type-specific expression of genes. Identifying the critical TFs in a given cellular context leads to investigating molecular regulatory mechanisms in development, differentiation, and disease. Because there are more than 1,500 human TFs, experimental measurements of genome-wide occupancy across all TFs have been challenging. While computational approaches play pivotal roles, most existing methods rely on statistical enrichment, focusing either on sequence motif similarity recognized by TFs or the similarity of the genomic region of interest with the previously characterized TF occupancy profile. Here we propose WhichTF as an alternative, incorporating curated biomedical knowledge from ontology and integrating it with the high-confidence prediction of conserved TF binding sites in user-provided genomic regions of interest. We develop a new WhichTF score to rank TFs and demonstrate its applicability across human and mouse cell types, cellular differentiation trajectories, and disease-associated cells.

Phylogeny of NF-YA trans-activation splicing isoforms in vertebrate evolution

NF-Y is a trimeric pioneer Transcription Factor (TF) whose target sequence -the CCAAT box- is present in ~25% of mammalian promoters. We reconstruct the phylogenetic history of the regulatory NF-YA subunit in vertebrates. We find that in addition to the remarkable conservation of the subunits-interaction and DNA-binding parts, the Transcriptional Activation Domain (TAD) is also conserved (>90% identity among bony vertebrates). We infer the phylogeny of the alternatively spliced exon-3 and partial splicing events of exon-7 -7N and 7C- revealing independent clade-specific losses of these regions. These isoforms shape the TAD. Absence of exon-3 in basal deuterostomes, cartilaginous fishes and hagfish, but not in lampreys, suggests that the "short" isoform is primordial, with emergence of exon-3 in chordates. Exon 7N was present in the vertebrate common ancestor, while 7C is a molecular innovation of teleost fishes. RNA-seq analysis in several species confirms expression of all these isoforms. We identify 3 blocks of amino acids in the TAD shared across deuterostomes, yet structural predictions and sequence analyses suggest an evolutionary drive for maintenance of an Intrinsically Disordered Region -IDR- within the TAD. Overall, these data help reconstruct the logic for alternative splicing of this essential eukaryotic TF.

The functional insight into the genetics of cardiovascular disease: results from the post-GWAS study

Cardiovascular diseases (CVDs), the leading cause of death worldwide, generally refer to a range of pathological conditions with the involvement of the heart and the blood vessels. A sizable fraction of the susceptibility loci is known, but the underlying mechanisms have been established only for a small proportion. Therefore, there is an increasing need to explore the functional relevance of trait-associated variants and, moreover, to search for novel risk genetic variation. We have reported the bioinformatic approach allowing effective identification of functional non-coding variants by integrated analysis of genome-wide data. Here, the analysis of 1361 previously identified regulatory SNPs (rSNPs) was performed to provide new insights into cardiovascular risk. We found 773,471 coding co-segregating markers for input rSNPs using the 1000 Genomes Project. The intersection of GWAS-derived SNPs with a relevance to cardiovascular traits with these markers was analyzed within a window of 10 Kbp. The effects on the transcription factor (TF) binding sites were explored by DeFine models. Functional pathway enrichment and protein– protein interaction (PPI) network analyses were performed on the targets and the extended genes by STRING and DAVID. Eighteen rSNPs were functionally linked to cardiovascular risk. A significant impact on binding sites of thirteen TFs including those involved in blood cells formation, hematopoiesis, macrophage function, inflammation, and vasoconstriction was found in K562 cells. 21 rSNP gene targets and 5 partners predicted by PPI were enriched for spliceosome and endocytosis KEGG pathways, endosome sorting complex and mRNA splicing REACTOME pathways. Related Gene Ontology terms included mRNA splicing and processing, endosome transport and protein catabolic processes. Together, the findings provide further insight into the biological basis of CVDs and highlight the importance of the precise regulation of splicing and alternative splicing.

Different SP1 binding dynamics at individual genomic loci in human cells

Using a tamoxifen-inducible time-course ChIP-sequencing (ChIP-seq) approach, we show that the ubiquitous transcription factor SP1 has different binding dynamics at its target sites in the human genome. SP1 very rapidly reaches maximal binding levels at some sites, but binding kinetics at other sites is biphasic, with rapid half-maximal binding followed by a considerably slower increase to maximal binding. While ∼70% of SP1 binding sites are located at promoter regions, loci with slow SP1 binding kinetics are enriched in enhancer and Polycomb-repressed regions. Unexpectedly, SP1 sites with fast binding kinetics tend to have higher quality and more copies of the SP1 sequence motif. Different cobinding factors associate near SP1 binding sites depending on their binding kinetics and on their location at promoters or enhancers. For example, NFY and FOS are preferentially associated near promoter-bound SP1 sites with fast binding kinetics, whereas DNA motifs of ETS and homeodomain proteins are preferentially observed at sites with slow binding kinetics. At promoters but not enhancers, proteins involved in sumoylation and PML bodies associate more strongly with slow SP1 binding sites than with the fast binding sites. The speed of SP1 binding is not associated with nucleosome occupancy, and it is not necessarily coupled to higher transcriptional activity. These results with SP1 are in contrast to those of human TBP, indicating that there is no common mechanism affecting transcription factor binding kinetics. The biphasic kinetics at some SP1 target sites suggest the existence of distinct chromatin states at these loci in different cells within the overall population.

Insulin Directly Regulates the Circadian Clock in Adipose Tissue

Adipose tissue (AT) is a key metabolic organ which functions are rhythmically regulated by an endogenous circadian clock. Feeding is a "zeitgeber" aligning the clock in AT with the external time, but mechanisms of this regulation remain largely unclear. We tested the hypothesis that postprandial changes of the hormone insulin directly entrain circadian clocks in AT and investigated a transcriptional-dependent mechanism of this regulation. We analyzed gene expression in subcutaneous AT (SAT) of obese subjects collected before and after the hyperinsulinemic-euglycemic clamp or control saline infusion (SC). The expressions of core clock genes PER2, PER3, and NR1D1 in SAT were differentially changed upon insulin and saline infusion, suggesting insulin-dependent clock regulation. In human stem cell-derived adipocytes, mouse 3T3-L1 cells, and AT explants from mPer2^Luc knockin mice, insulin induced a transient increase of the Per2 mRNA and protein expression, leading to the phase shift of circadian oscillations, with similar effects for Per1 Insulin effects were dependent on the region between -64 and -43 in the Per2 promoter but not on CRE and E-box elements. Our results demonstrate that insulin directly regulates circadian clocks in AT and isolated adipocytes, thus representing a primary mechanism of feeding-induced AT clock entrainment.

Activation of γ-globin gene expression by GATA1 and NF-Y in hereditary persistence of fetal hemoglobin

Hereditary persistence of fetal hemoglobin (HPFH) ameliorates β-hemoglobinopathies by inhibiting the developmental switch from γ-globin (HBG1/HBG2) to β-globin (HBB) gene expression. Some forms of HPFH are associated with γ-globin promoter variants that either disrupt binding motifs for transcriptional repressors or create new motifs for transcriptional activators. How these variants sustain γ-globin gene expression postnatally remains undefined. We mapped γ-globin promoter sequences functionally in erythroid cells harboring different HPFH variants. Those that disrupt a BCL11A repressor binding element induce γ-globin expression by facilitating the recruitment of transcription factors NF-Y to a nearby proximal CCAAT box and GATA1 to an upstream motif. The proximal CCAAT element becomes dispensable for HPFH variants that generate new binding motifs for activators NF-Y or KLF1, but GATA1 recruitment remains essential. Our findings define distinct mechanisms through which transcription factors and their cis-regulatory elements activate γ-globin expression in different forms of HPFH, some of which are being recreated by therapeutic genome editing.

Mechanisms involved in selecting and maintaining neuroblastoma cancer stem cell populations, and perspectives for therapeutic targeting

Pediatric neuroblastomas (NBs) are heterogeneous, aggressive, therapy-resistant embryonal tumours that originate from cells of neural crest (NC) origin and in particular neuroblasts committed to the sympathoadrenal progenitor cell lineage. Therapeutic resistance, post-therapeutic relapse and subsequent metastatic NB progression are driven primarily by cancer stem cell (CSC)-like subpopulations, which through their self-renewing capacity, intermittent and slow cell cycles, drug-resistant and reversibly adaptive plastic phenotypes, represent the most important obstacle to improving therapeutic outcomes in unfavourable NBs. In this review, dedicated to NB CSCs and the prospects for their therapeutic eradication, we initiate with brief descriptions of the unique transient vertebrate embryonic NC structure and salient molecular protagonists involved NC induction, specification, epithelial to mesenchymal transition and migratory behaviour, in order to familiarise the reader with the embryonic cellular and molecular origins and background to NB. We follow this by introducing NB and the potential NC-derived stem/progenitor cell origins of NBs, before providing a comprehensive review of the salient molecules, signalling pathways, mechanisms, tumour microenvironmental and therapeutic conditions involved in promoting, selecting and maintaining NB CSC subpopulations, and that underpin their therapy-resistant, self-renewing metastatic behaviour. Finally, we review potential therapeutic strategies and future prospects for targeting and eradication of these bastions of NB therapeutic resistance, post-therapeutic relapse and metastatic progression.

The cis-regulatory effects of modern human-specific variants

The Neanderthal and Denisovan genomes enabled the discovery of sequences that differ between modern and archaic humans, the majority of which are noncoding. However, our understanding of the regulatory consequences of these differences remains limited, in part due to the decay of regulatory marks in ancient samples. Here, we used a massively parallel reporter assay in embryonic stem cells, neural progenitor cells, and bone osteoblasts to investigate the regulatory effects of the 14,042 single-nucleotide modern human-specific variants. Overall, 1791 (13%) of sequences containing these variants showed active regulatory activity, and 407 (23%) of these drove differential expression between human groups. Differentially active sequences were associated with divergent transcription factor binding motifs, and with genes enriched for vocal tract and brain anatomy and function. This work provides insight into the regulatory function of variants that emerged along the modern human lineage and the recent evolution of human gene expression.

HAT cofactor TRRAP modulates microtubule dynamics via SP1 signaling to prevent neurodegeneration

Brain homeostasis is regulated by the viability and functionality of neurons. HAT (histone acetyltransferase) and HDAC (histone deacetylase) inhibitors have been applied to treat neurological deficits in humans; yet, the epigenetic regulation in neurodegeneration remains elusive. Mutations of HAT cofactor TRRAP (transformation/transcription domain-associated protein) cause human neuropathies, including psychosis, intellectual disability, autism, and epilepsy, with unknown mechanism. Here we show that Trrap deletion in Purkinje neurons results in neurodegeneration of old mice. Integrated transcriptomics, epigenomics, and proteomics reveal that TRRAP via SP1 conducts a conserved transcriptomic program. TRRAP is required for SP1 binding at the promoter proximity of target genes, especially microtubule dynamics. The ectopic expression of Stathmin3/4 ameliorates defects of TRRAP-deficient neurons, indicating that the microtubule dynamics is particularly vulnerable to the action of SP1 activity. This study unravels a network linking three well-known, but up-to-date unconnected, signaling pathways, namely TRRAP, HAT, and SP1 with microtubule dynamics, in 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.

Expression of the zinc finger transcription factor Sp6-9 in the velvet worm Euperipatoides kanangrensis suggests a conserved role in appendage development in Panarthropoda

The Sp-family genes encode important transcription factors in animal development. Here we investigate the embryonic expression patterns of the complete set of Sp-genes in the velvet worm Euperipatoides kanangrensis (Onychophora), with a special focus on the Sp6-9 ortholog. In arthropods, Sp6-9, the ortholog of the Drosophila melanogaster D-Sp1 gene plays a conserved role in appendage development. Our data show that the expression of Sp6-9 during the development of the velvet worm is conserved, suggesting that the key function of the Sp6-9 gene dates back to at least the last common ancestor of arthropods and onychophorans and thus likely the last common ancestor of Panarthropoda.

Discovery, characterization and potential roles of a novel NF-YAx splice variant in human neuroblastoma

Background

Identification of novel cancer-associated splice variants is of potential diagnostic, prognostic and therapeutic importance. NF-Y transcription factor is comprised of NF-YA, NF-YB and NF-YC subunits, binds inverted CCAAT-boxes in ≈70% of gene promoters, regulates > 1000 cancer-associated genes and proteins involved in proliferation, staminality, differentiation, apoptosis, metabolism and is subject to component alternative splicing. RT-PCR evaluation of alternative NF-YA splicing in primary human neuroblastomas (NBs), led to discovery of a novel NF-YAx splice variant, also expressed during mouse embryo development and induced by doxorubicin in NB cells. Here, we report the discovery and characterisation of NF-YAx and discus its potential roles in NB.

Methods

NF-YAx cDNA was RT-PCR-cloned from a stage 3 NB (provided by the Italian Association of Haematology and Paediatric Oncology, Genova, IT), sequenced and expressed as a protein using standard methods and compared to known fully-spliced NF-YAl and exon B-skipped NF-YAs isoforms in: EMSAs for capacity to form NF-Y complexes; by co-transfection, co-immunoprecipitation and Western blotting for capacity to bind Sp1; by IF for localisation; in AO/EtBr cell-death and colony formation assays for relative cytotoxicity, and by siRNA knockdown, use of inhibitors and Western blotting for potential mechanisms of action. Stable SH-SY5Y transfectants of all three NF-YA isoforms were also propagated and compared by RT-PCR and Western blotting for differences in cell-death and stem cell (SC)-associated gene expression, in cell-death assays for sensitivity to doxorubicin and in in vitro proliferation, substrate-independent growth and in vivo tumour xenograft assays for differences in growth and tumourigenic capacity.

Results

NF-YAx was characterized as a novel variant with NF-YA exons B, D and partial F skipping, detected in 20% of NF-YA positive NBs, was the exclusive isoform in a stage 3 NB, expressed in mouse stage E11.5–14 embryos and induced by doxorubicin in SH-SY5Y NB cells. The NF-YAx protein exhibited nuclear localisation, competed with other isoforms in CCAAT box-binding NF-Y complexes but, in contrast to other isoforms, did not bind Sp1. NF-YAx expression in neural-related progenitor and NB cells repressed Bmi1 expression, induced KIF1Bβ expression and promoted KIF1Bβ-dependent necroptosis but in NB cells also selected tumourigenic, doxorubicin-resistant, CSC-like sub-populations, resistant to NF-YAx cytotoxicity.

Conclusions

The discovery of NF-YAx in NBs, its expression in mouse embryos and induction by doxorubicin in NB cells, unveils a novel NF-YA splice mechanism and variant, regulated by and involved in development, genotoxic-stress and NB. NF-YAx substitution of other isoforms in NF-Y complexes and loss of capacity to bind Sp1, characterises this novel isoform as a functional modifier of NF-Y and its promotion of KIF1Bβ-dependent neural-lineage progenitor and NB cell necroptosis, association with doxorubicin-induced necroptosis and expression in mouse embryos coinciding with KIF1Bβ-dependent sympathetic neuroblast-culling, confirm a cytotoxic function and potential role in suppressing NB initiation. On the other hand, the in vitro selection of CSC-like NB subpopulations resistant to NF-YAx cytotoxicity not only helps to explain high-level exclusive NF-YAx expression in a stage 3 NB but also supports a role for NF-YAx in disease progression and identifies a potential doxorubicin-inducible mechanism for post-therapeutic relapse.

Overexpression and alternative splicing of NF-YA in breast cancer

NF-Y is a CCAAT-binding trimeric transcription factor, whose regulome, interactome and oncogenic potential point to direct involvement in cellular transformation. Yet little is known about the levels of NF-Y subunits in tumors. We focused on breast carcinomas, and analyzed RNA-Seq datasets of TCGA and 54 BRCA cell lines at gene and isoforms level. We partitioned all tumors in the four major subclasses. NF-YA, but not histone-fold subunits NF-YB/NF-YC, is globally overexpressed, correlating with the proliferative Ki67 marker and a common set of 840 genes, with cell-cycle, metabolism GO terms. Their promoters are enriched in NF-Y, GC-rich and E2F sites. Surprisingly, there is an isoform switch, with the “short” isoform -NF-YAs- becoming predominant in tumors. E2F genes are also overexpressed in BRCA, but no switch in isoforms is observed. In Basal-like Claudin^low cell lines and tumors, expression of NF-YAl -long- isoform is high, together with 11 typical EMT markers and low levels of basal Keratins. Analysis of Progression-Free-Intervals indicates that tumors with unbalance of NF-YA isoforms ratios have worst clinical outcomes. The data suggest that NF-YA overexpression increases CCAAT-dependent, pro-growth genes in BRCA. NF-YAs is associated with a proliferative signature, but high levels of NF-YAl signal loss of epithelial features, EMT and acquisition of a more aggressive behavior in a subset of Claudin^low Basal-like tumors.

Curcumin ameliorates PRMT5-MEP50 arginine methyltransferase expression by decreasing the Sp1 and NF-YA transcription factors in the A549 and MCF-7 cells

The protein arginine methyltransferase 5 (PRMT5) and its catalytic partner methylosome protein MEP50 (WDR77) catalyse the mono- and symmetric di-methylation of selective arginines in various histones and non-histone target proteins. It has emerged as a crucial epigenetic regulator in cell proliferation and differentiation; which also reported to be overexpressed in many forms of cancers in humans. In this study, we aimed to assess the modulations in the expression of this enzyme upon exposure to the well-studied natural compound from the spice turmeric, curcumin. We exposed the lung and breast cancer cell lines (A549 and MCF-7) to curcumin (2 and 20 μM) and observed a highly significant inhibitory effect on the expression of both PRMT5 and MEP50. The level of symmetrical dimethylarginine (SDMA) in multiple proteins, and more specifically, the H4R3me2s mark (which predominates in GC-rich motifs in nucleosomal DNA) was also diminished significantly. We also found that curcumin significantly reduced the level and enrichment of the transcription factors Sp1 and NF-YA which shares their binding sites within the GC-rich region of the PRMT5 proximal promoter. Furthermore, the involvement of both PKC-p38-ERK-cFos and AKT-mTOR signalling was observed in reducing the Sp1 and NF-YA expression by curcumin. Therefore, we propose curcumin decreased the expression of PRMT5 in these cells by affecting at least these two transcription factors. Altogether, we report a new molecular target of curcumin and further elucidation of this proposed mechanism through which curcumin affects the PRMT5-MEP50 methyltransferase expression might be explored for its therapeutic application.

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.

The phosphorylatable Ser320 of NF-YA is involved in DNA binding of the NF-Y trimer

Nuclear factor Y (NF-Y) is a transcription factor trimer binding to the functionally important CCAAT box, present in promoters of growth-promoting and cell cycle-regulated genes. The regulatory nuclear factor YA (NF-YA) subunit confers sequence-specificity to the histone-like nuclear factor YB/YC dimer. NF-YA harbors 2 serines-Ser320 and Ser326-shown to be phosphorylated by cyclin-dependent kinase 2. High-throughput proteomics data indicate that they are phosphorylated in vivo. Specifically, Ser320 makes structural contacts with the DNA phosphate backbone; Ser320-P is the major NF-YA phosphorylation isoform following overexpression in HeLa cells, increasing upon mitotic arrest. EMSA with recombinant Ala and Glu mutants confirm a role of Ser320, but not Ser326, in stabilization of DNA binding. Transactivation assays of the CCAAT-dependent MDR1 and RHOB promoters show loss in transcription function for Ser320Glu and Ser320Ala NF-YA mutants. Phylogenetic analysis of NF-YA proteins indicates that Ser320 is indeed evolutionarily conserved. We conclude that phosphorylation of this residue belongs to the core mechanisms of DNA-binding control, possibly driven by the necessity to unfasten binding of or to evict NF-Y from CCAAT sites under specific conditions of growth regulation.-Bernardini, A., Lorenzo, M., Nardini, M., Mantovani, R., Gnesutta, N. The phosphorylatable Ser320 of NF-YA is involved in DNA binding of the NF-Y trimer.

Transcription factor Sp2 potentiates binding of the TALE homeoproteins Pbx1:Prep1 and the histone-fold domain protein Nf-y to composite genomic sites

Different transcription factors operate together at promoters and enhancers to regulate gene expression. Transcription factors either bind directly to their target DNA or are tethered to it by other proteins. The transcription factor Sp2 serves as a paradigm for indirect genomic binding. It does not require its DNA-binding domain for genomic DNA binding and occupies target promoters independently of whether they contain a cognate DNA-binding motif. Hence, Sp2 is strikingly different from its closely related paralogs Sp1 and Sp3, but how Sp2 recognizes its targets is unknown. Here, we sought to gain more detailed insights into the genomic targeting mechanism of Sp2. ChIP-exo sequencing in mouse embryonic fibroblasts revealed genomic binding of Sp2 to a composite motif where a recognition sequence for TALE homeoproteins and a recognition sequence for the trimeric histone-fold domain protein nuclear transcription factor Y (Nf-y) are separated by 11 bp. We identified a complex consisting of the TALE homeobox protein Prep1, its partner PBX homeobox 1 (Pbx1), and Nf-y as the major partners in Sp2-promoter interactions. We found that the Pbx1:Prep1 complex together with Nf-y recruits Sp2 to co-occupied regulatory elements. In turn, Sp2 potentiates binding of Pbx1:Prep1 and Nf-y. We also found that the Sp-box, a short sequence motif close to the Sp2 N terminus, is crucial for Sp2's cofactor function. Our findings reveal a mechanism by which the DNA binding-independent activity of Sp2 potentiates genomic loading of Pbx1:Prep1 and Nf-y to composite motifs present in many promoters of highly expressed genes.

Enhancement of transgene expression by nuclear transcription factor Y and CCCTC-binding factor

If a transgene is effectively delivered to a cell, its expression may still be limited by epigenetic mechanisms that silence the transgene. Indeed, once the transgene reaches the nucleus, it may be bound by histone proteins and condensed into heterochromatin or associated with repressor proteins that block transcription. In this study, we sought to enhance transgene expression by adding binding motifs for several different epigenetic enzymes either upstream or downstream of two promoters (CMV and EF1α). Screening these plasmids revealed that luciferase expression was enhanced 10-fold (10.4 ± 5.8) by the addition of a CCAAT box just upstream of the EF1α promoter to recruit nuclear transcription factor Y (NF-Y), while inserting a CCCTC-binding factor (CTCF) motif downstream of the EF1α promoter enhanced expression at least 14-fold (14.03 ± 6.54). ChIP assays confirmed that NF-Y and CTCF bound to the motifs that were added to each plasmid, but the presence of NF-Y and CTCF did not significantly affect the levels of histone acetylation (H3K9ac) or methylation (H3K9me3). Overall, these results show that transgene expression from the EF1α promoter can be significantly increased with motifs that recruit NF-Y or CTCF. © 2018 American Institute of Chemical Engineers Biotechnol. Prog., 34:1581-1588, 2018.

Virtual Sorting Has a Distinctive Advantage in Identification of Anticorrelated Genes and Further Negative Regulators of Immune Cell Subpopulations

Immune cells are highly plastic in both gene expression and cell phenotype. We have established a method of gene expressional plasticity and virtual sorting to evaluate immune cell subpopulations and their characteristic genes in human CD4⁺ T cells. In this study, we continued to investigate the informatics mechanism on the effectiveness of virtual sorting. We found that virtual sorting had an overall positive correlation to the Pearson correlation in the identification of positively correlated genes. However, owing to nonlinear biological anticorrelation, virtual sorting showed a distinctive advantage for anticorrelated genes, suggesting an important role in the identification of negative regulators. In addition, based on virtual sorting results, we identified two basic gene sets among highly plastic genes, i.e., highly plastic cell cycle-associated molecules and highly plastic immune and defense response-associated molecules. Genes within each set tended to be positively connected, but genes between two sets were often anticorrelated. Further analysis revealed preferential transcription factor binding motifs existed between highly plastic cell cycle-associated molecules and highly plastic immune and defense response-associated molecules. Our results strongly suggested predetermined regulation, which was called an immune cell internal phenotype, should exist and could be mined by virtual sorting analysis. This provided efficient functional clues to study immune cell phenotypes and their regulation. Moreover, the current substantial virtual sorting results in both CD4⁺ T cells and B cells provide a useful resource for big-data-driven experimental studies and knowledge discoveries.

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.