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

Overexpression of the PtrNF-YA6 gene inhibits secondary cell wall thickening in poplar

The NF-Y gene family in plants plays a crucial role in numerous biological processes, encompassing hormone response, stress response, as well as growth and development. In this study, we first used bioinformatics techniques to identify members of the NF-YA family that may function in wood formation. We then used molecular biology techniques to investigate the role and molecular mechanism of PtrNF-YA6 in secondary cell wall (SCW) formation in Populus trichocarpa. We found that PtrNF-YA6 protein was localized in the nucleus and had no transcriptional activating activity. Overexpression of PtrNF-YA6 had an inhibitory effect on plant growth and development and significantly suppressed hemicellulose synthesis and SCW thickening in transgenic plants. Yeast one-hybrid and ChIP-PCR assays revealed that PtrNF-YA6 directly regulated the expression of hemicellulose synthesis genes (PtrGT47A-1, PtrGT8C, PtrGT8F, PtrGT43B, PtrGT47C, PtrGT8A and PtrGT8B). In conclusion, PtrNF-YA6 can inhibit plant hemicellulose synthesis and SCW thickening by regulating the expression of downstream SCW formation-related target genes.

The pancancer overexpressed NFYC Antisense 1 controls cell cycle mitotic progression through in cis and in trans modes of action

Antisense RNAs (asRNAs) represent an underappreciated yet crucial layer of gene expression regulation. Generally thought to modulate their sense genes in cis through sequence complementarity or their act of transcription, asRNAs can also regulate different molecular targets in trans, in the nucleus or in the cytoplasm. Here, we performed an in-depth molecular characterization of NFYC Antisense 1 (NFYC-AS1), the asRNA transcribed head-to-head to NFYC subunit of the proliferation-associated NF-Y transcription factor. Our results show that NFYC-AS1 is a prevalently nuclear asRNA peaking early in the cell cycle. Comparative genomics suggests a narrow phylogenetic distribution, with a probable origin in the common ancestor of mammalian lineages. NFYC-AS1 is overexpressed pancancer, preferentially in association with RB1 mutations. Knockdown of NFYC-AS1 by antisense oligonucleotides impairs cell growth in lung squamous cell carcinoma and small cell lung cancer cells, a phenotype recapitulated by CRISPR/Cas9-deletion of its transcription start site. Surprisingly, expression of the sense gene is affected only when endogenous transcription of NFYC-AS1 is manipulated. This suggests that regulation of cell proliferation is at least in part independent of the in cis transcription-mediated effect on NFYC and is possibly exerted by RNA-dependent in trans effects converging on the regulation of G2/M cell cycle phase genes. Accordingly, NFYC-AS1-depleted cells are stuck in mitosis, indicating defects in mitotic progression. Overall, NFYC-AS1 emerged as a cell cycle-regulating asRNA with dual action, holding therapeutic potential in different cancer types, including the very aggressive RB1-mutated tumors.

Expression and function of NF-Y subunits in cancer

NF-Y is a Transcription Factor (TF) targeting the CCAAT box regulatory element. It consists of the NF-YB/NF-YC heterodimer, each containing an Histone Fold Domain (HFD), and the sequence-specific subunit NF-YA. NF-YA expression is associated with cell proliferation and absent in some post-mitotic cells. The review summarizes recent findings impacting on cancer development. The logic of the NF-Y regulome points to pro-growth, oncogenic genes in the cell-cycle, metabolism and transcriptional regulation routes. NF-YA is involved in growth/differentiation decisions upon cell-cycle re-entry after mitosis and it is widely overexpressed in tumors, the HFD subunits in some tumor types or subtypes. Overexpression of NF-Y -mostly NF-YA- is oncogenic and decreases sensitivity to anti-neoplastic drugs. The specific roles of NF-YA and NF-YC isoforms generated by alternative splicing -AS- are discussed, including the prognostic value of their levels, although the specific molecular mechanisms of activity are still to be deciphered.

Natural medicinal compounds target signal transduction pathways to overcome ABC drug efflux transporter-mediated multidrug resistance in cancer

ATP-binding cassette (ABC) transporters such as ABCB1, ABCG2, and ABCC1 are the major players in drug efflux-mediated multidrug resistance (MDR), which severely affects the efficacy of chemotherapy. Several synthetic compounds block the drug transport by ABC transporters; however, they exhibit a narrow therapeutic window, and produce side effects in non-target normal tissues. Conversely, the downregulation of the expression of ABC drug transporters seems to be a promising strategy to reverse MDR in cancer cells. Several signaling pathways, such as NF-κB, STAT3, Gli, NICD, YAP/TAZ, and Nrf2 upregulate the expression of ABC drug transporters in drug-resistant cancers. Recently, natural medicinal compounds have gained importance to overcome the ABC drug-efflux pump-mediated MDR in cancer. These compounds target transcription factors and the associated signal transduction pathways, thereby downregulating the expression of ABC transporters in drug-resistant cancer cells. Several potent natural compounds have been identified as lead candidates to synergistically enhance chemotherapeutic efficacy, and a few of them are already in clinical trials. Therefore, modulation of signal transduction pathways using natural medicinal compounds for the reversal of ABC drug transporter-mediated MDR in cancer is a novel approach for improving the efficiency of the existing chemotherapeutics. In this review, we discuss the modulatory role of natural medicinal compounds on cellular signaling pathways that regulate the expression of ABC transporters in drug-resistant cancer cells.

Myo-differentiation reporter screen reveals NF-Y as an activator of PAX3-FOXO1 in rhabdomyosarcoma

Recurrent chromosomal rearrangements found in rhabdomyosarcoma (RMS) produce the PAX3-FOXO1 fusion protein, which is an oncogenic driver and a dependency in this disease. One important function of PAX3-FOXO1 is to arrest myogenic differentiation, which is linked to the ability of RMS cells to gain an unlimited proliferation potential. Here, we developed a phenotypic screening strategy for identifying factors that collaborate with PAX3-FOXO1 to block myo-differentiation in RMS. Unlike most genes evaluated in our screen, we found that loss of any of the three subunits of the Nuclear Factor Y (NF-Y) complex leads to a myo-differentiation phenotype that resembles the effect of inactivating PAX3-FOXO1. While the transcriptomes of NF-Y- and PAX3-FOXO1-deficient RMS cells bear remarkable similarity to one another, we found that these two transcription factors occupy nonoverlapping sites along the genome: NF-Y preferentially occupies promoters, whereas PAX3-FOXO1 primarily binds to distal enhancers. By integrating multiple functional approaches, we map the PAX3 promoter as the point of intersection between these two regulators. We show that NF-Y occupies CCAAT motifs present upstream of PAX3 to function as a transcriptional activator of PAX3-FOXO1 expression in RMS. These findings reveal a critical upstream role of NF-Y in the oncogenic PAX3-FOXO1 pathway, highlighting how a broadly essential transcription factor can perform tumor-specific roles in governing cellular state.

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.

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 USR domain of USF1 mediates NF-Y interactions and cooperative DNA binding

The trimeric CCAAT-binding NF-Y is a "pioneer" Transcription Factor -TF- known to cooperate with neighboring TFs to regulate gene expression. Genome-wide analyses detected a precise stereo-alignment -10/12 bp- of CCAAT with E-box elements and corresponding colocalization of NF-Y with basic-Helix-Loop-Helix (bHLH) TFs. We dissected here NF-Y interactions with USF1 and MAX. USF1, but not MAX, cooperates in DNA binding with NF-Y. NF-Y and USF1 synergize to activate target promoters. Reconstruction of complexes by structural means shows independent DNA binding of MAX, whereas USF1 has extended contacts with NF-Y, involving the USR, a USF-specific amino acid sequence stretch required for trans-activation. The USR is an intrinsically disordered domain and adopts different conformations based on E-box-CCAAT distances. Deletion of the USR abolishes cooperative DNA binding with NF-Y. Our data indicate that the functionality of certain unstructured domains involves adapting to small variation in stereo-alignments of the multimeric TFs sites.

Characterization of genes associated with TGA7 during the floral transition

BACKGROUND

The TGACG-binding (TGA) family has 10 members that play vital roles in Arabidopsis thaliana defense responses and development. However, their involvement in controlling flowering time remains largely unknown and requires further investigation.

RESULTS

To study the role of TGA7 during floral transition, we first investigated the tga7 mutant, which displayed a delayed-flowering phenotype under both long-day and short-day conditions. We then performed a flowering genetic pathway analysis and found that both autonomous and thermosensory pathways may affect TGA7 expression. Furthermore, to reveal the differential gene expression profiles between wild-type (WT) and tga7, cDNA libraries were generated for WT and tga7 mutant seedlings at 9 days after germination. For each library, deep-sequencing produced approximately 6.67 Gb of high-quality sequences, with the majority (84.55 %) of mRNAs being between 500 and 3,000 nt. In total, 325 differentially expressed genes were identified between WT and tga7 mutant seedlings. Among them, four genes were associated with flowering time control. The differential expression of these four flowering-related genes was further validated by qRT-PCR.

CONCLUSIONS

Among these four differentially expressed genes associated with flowering time control, FLC and MAF5 may be mainly responsible for the delayed-flowering phenotype in tga7, as TGA7 expression was regulated by autonomous pathway genes. These results provide a framework for further studying the role of TGA7 in promoting flowering.

Structural Variations of the 3D Genome Architecture in Cervical Cancer Development

Human papillomavirus (HPV) integration is the major contributor to cervical cancer (CC) development by inducing structural variations (SVs) in the human genome. SVs are directly associated with the three-dimensional (3D) genome structure leading to cancer development. The detection of SVs is not a trivial task, and several genome-wide techniques have greatly helped in the identification of SVs in the cancerous genome. However, in cervical cancer, precise prediction of SVs mainly translocations and their effects on 3D-genome and gene expression still need to be explored. Here, we have used high-throughput chromosome conformation capture (Hi-C) data of cervical cancer to detect the SVs, especially the translocations, and validated it through whole-genome sequencing (WGS) data. We found that the cervical cancer 3D-genome architecture rearranges itself as compared to that in the normal tissue, and 24% of the total genome switches their A/B compartments. Moreover, translocation detection from Hi-C data showed the presence of high-resolution t(4;7) (q13.1; q31.32) and t(1;16) (q21.2; q22.1) translocations, which disrupted the expression of the genes located at and nearby positions. Enrichment analysis suggested that the disrupted genes were mainly involved in controlling cervical cancer-related pathways. In summary, we detect the novel SVs through Hi-C data and unfold the association among genome-reorganization, translocations, and gene expression regulation. The results help understand the underlying pathogenicity mechanism of SVs in cervical cancer development and identify the targeted therapeutics against cervical cancer.

NF-Y Subunits Overexpression in HNSCC

Simple Summary

Cancer cells have altered gene expression profiles. This is ultimately elicited by altered structure, expression or binding of transcription factors to regulatory regions of genomes. The CCAAT-binding trimer is a pioneer transcription factor involved in the activation of “cancer” genes. We and others have shown that the regulatory NF-YA subunit is overexpressed in epithelial cancers. Here, we examined large datasets of bulk gene expression profiles, as well as single-cell data, in head and neck squamous cell carcinomas by bioinformatic methods. We partitioned tumors according to molecular subtypes, mutations and positivity for HPV. We came to the conclusion that high levels of the histone-like subunits and the “short” NF-YAs isoform are protective in HPV-positive tumors. On the other hand, high levels of the “long” NF-YAl were found in the recently identified aggressive and metastasis-prone cell population undergoing partial epithelial to mesenchymal transition, p-EMT.

Abstract

NF-Y is the CCAAT-binding trimer formed by the histone fold domain (HFD), NF-YB/NF-YC and NF-YA. The CCAAT box is generally prevalent in promoters of “cancer” genes. We reported the overexpression of NF-YA in BRCA, LUAD and LUSC, and of all subunits in HCC. Altered splicing of NF-YA was found in breast and lung cancer. We analyzed RNA-seq datasets of TCGA and cell lines of head and neck squamous cell carcinomas (HNSCC). We partitioned all TCGA data into four subtypes, deconvoluted single-cell RNA-seq of tumors and derived survival curves. The CCAAT box was enriched in the promoters of overexpressed genes. The “short” NF-YAs was overexpressed in all subtypes and the “long” NF-YAl in Mesenchymal. The HFD subunits are overexpressed, except Basal (NF-YB) and Atypical (NF-YC); NF-YAl is increased in p53 mutated tumors. In HPV-positive tumors, high levels of NF-YAs, p16 and ΔNp63 correlate with better prognosis. Deconvolution of single cell RNA-seq (scRNA-seq) found a correlation of NF-YAl with Cancer Associated Fibroblasts (CAFs) and p-EMT cells, a population endowed with metastatic potential. We conclude that overexpression of HFD subunits and NF-YAs is protective in HPV-positive tumors; expression of NF-YAl is largely confined to mutp53 tumors and malignant p-EMT cells.

Live cell dynamics of the NF-Y transcription factor

Transcription factors (TFs) are core players in the control of gene expression, evolutionarily selected to recognise a subset of specific DNA sequences and nucleate the recruitment of the transcriptional machinery. How TFs assemble and move in the nucleus to locate and bind their DNA targets and cause a transcriptional response, remains mostly unclear. NF-Y is a highly conserved, heterotrimeric TF with important roles in both housekeeping and lineage-specific gene expression, functioning as a promoter organiser. Despite a large number of biochemical, structural and genomic studies of NF-Y, there is a lack of experiments in single living cells; therefore, basic assumptions of NF-Y biology remain unproven in vivo. Here we employ a series of dynamic fluorescence microscopy methods (FLIM-FRET, NB, RICS and FRAP) to study NF-Y dynamics and complex formation in live cells. Specifically, we provide quantitative measurement of NF-Y subunit association and diffusion kinetics in the nucleus that collectively suggest NF-Y to move and bind chromatin as a trimeric complex in vivo.

Transcription factor competition at the γ-globin promoters controls hemoglobin switching

BCL11A, the major regulator of fetal hemoglobin (HbF, α₂γ₂) level, represses γ-globin expression through direct promoter binding in adult erythroid cells in a switch to adult hemoglobin (HbA, α₂β₂). To uncover how BCL11A initiates repression, we used CRISPR-Cas9, dCas9, dCas9-KRAB and dCas9-VP64 screens to dissect the γ-globin promoters and identified an activator element near the BCL11A-binding site. Using CUT&RUN and base editing, we demonstrate that a proximal CCAAT box is occupied by the activator NF-Y. BCL11A competes with NF-Y binding through steric hindrance to initiate repression. Occupancy of NF-Y is rapidly established following BCL11A depletion, and precedes γ-globin derepression and locus control region (LCR)-globin loop formation. Our findings reveal that the switch from fetal to adult globin gene expression within the >50-kb β-globin gene cluster is initiated by competition between a stage-selective repressor and a ubiquitous activating factor within a remarkably discrete region of the γ-globin promoters.

Deep neural networks identify sequence context features predictive of transcription factor binding

Transcription factors (TFs) bind DNA by recognizing specific sequence motifs, typically of length 6-12bp. A motif can occur many thousands of times in the human genome, but only a subset of those sites are actually bound. Here we present a machine learning framework leveraging existing convolutional neural network architectures and model interpretation techniques to identify and interpret sequence context features most important for predicting whether a particular motif instance will be bound. We apply our framework to predict binding at motifs for 38 TFs in a lymphoblastoid cell line, score the importance of context sequences at base-pair resolution, and characterize context features most predictive of binding. We find that the choice of training data heavily influences classification accuracy and the relative importance of features such as open chromatin. Overall, our framework enables novel insights into features predictive of TF binding and is likely to inform future deep learning applications to interpret non-coding genetic variants.

Synergism of Proneurogenic miRNAs Provides a More Effective Strategy to Target Glioma Stem Cells

Simple Summary

miRNAs function as critical regulators of gene expression and have been defined as contributors of cancer phenotypes by acting as oncogenes or tumor suppressors. Based on these findings, miRNA-based therapies have been explored in the treatment of many different malignancies. The use of single miRNAs has faced some challenges and showed limited success. miRNAs cooperate to regulate distinct biological processes and pathways and, therefore, combination of related miRNAs could amplify the repression of oncogenic factors and the effect on cancer relevant pathways. We established that the combination of tumor suppressor miRNAs miR-124, miR-128, and miR-137 is much more effective than single miRNAs in disrupting proliferation and survival of glioma stem cells and neuroblastoma lines and promoting differentiation and response to radiation. Subsequent genomic analyses showed that other combinations of tumor suppressor miRNAs could be equally effective, and its use could provide new routes to target in special cancer-initiating cell populations.

Abstract

Tumor suppressor microRNAs (miRNAs) have been explored as agents to target cancer stem cells. Most strategies use a single miRNA mimic and present many disadvantages, such as the amount of reagent required and the diluted effect on target genes. miRNAs work in a cooperative fashion to regulate distinct biological processes and pathways. Therefore, we propose that miRNA combinations could provide more efficient ways to target cancer stem cells. We have previously shown that miR-124, miR-128, and miR-137 function synergistically to regulate neurogenesis. We used a combination of these three miRNAs to treat glioma stem cells and showed that this treatment was much more effective than single miRNAs in disrupting cell proliferation and survival and promoting differentiation and response to radiation. Transcriptomic analyses indicated that transcription regulation, angiogenesis, metabolism, and neuronal differentiation are among the main biological processes affected by transfection of this miRNA combination. In conclusion, we demonstrated the value of using combinations of neurogenic miRNAs to disrupt cancer phenotypes and glioma stem cell growth. The synergistic effect of these three miRNA amplified the repression of oncogenic factors and the effect on cancer relevant pathways. Future therapeutic approaches would benefit from utilizing miRNA combinations, especially when targeting cancer-initiating cell populations.

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.

NF-Y Overexpression in Liver Hepatocellular Carcinoma (HCC)

NF-Y is a pioneer trimeric transcription factor formed by the Histone Fold Domain (HFD) NF-YB/NF-YC subunits and NF-YA. Three subunits are required for DNA binding. CCAAT-specificity resides in NF-YA and transactivation resides in Q-rich domains of NF-YA and NF-YC. They are involved in alternative splicing (AS). We recently showed that NF-YA is overexpressed in breast and lung carcinomas. We report here on the overexpression of all subunits in the liver hepatocellular carcinoma (HCC) TCGA database, specifically the short NF-YAs and NF-YC2 (37 kDa) isoforms. This is observed at all tumor stages, in viral-infected samples and independently from the inflammatory status. Up-regulation of NF-YAs and NF-YC, but not NF-YB, is associated to tumors with mutant p53. We used a deep-learning-based method (DeepCC) to extend the partitioning of the three molecular clusters to all HCC TCGA tumors. In iCluster3, CCAAT is a primary matrix found in promoters of up-regulated genes, and cell-cycle pathways are enriched. Finally, clinical data indicate that, globally, only NF-YAs, but not HFD subunits, correlate with the worst prognosis; in iCluster1 patients, however, all subunits correlate. The data show a difference with other epithelial cancers, in that global overexpression of the three subunits is reported and clinically relevant in a subset of patients; yet, they further reinstate the regulatory role of the sequence-specific subunit.

Influenza A virus directly modulates mouse eosinophil responses

Allergic asthma and influenza are common respiratory diseases with a high probability of co-occurrence. During the 2009 influenza pandemic, hospitalized patients with influenza experienced lower morbidity if asthma was an underlying condition. We have previously demonstrated that acute allergic asthma protects mice from severe influenza and have implicated eosinophils in the airways of mice with allergic asthma as participants in the antiviral response. However, very little is known about how eosinophils respond to direct exposure to influenza A virus (IAV) or the microenvironment in which the viral burden is high. We hypothesized that eosinophils would dynamically respond to the presence of IAV through phenotypic, transcriptomic, and physiologic changes. Using our mouse model of acute fungal asthma and influenza, we showed that eosinophils in lymphoid tissues were responsive to IAV infection in the lungs and altered surface expression of various markers necessary for cell activation in a niche-specific manner. Siglec-F expression was altered in a subset of eosinophils after virus exposure, and those expressing high Siglec-F were more active (IL-5Rα^hi CD62L^lo ). While eosinophils exposed to IAV decreased their overall transcriptional activity and mitochondrial oxygen consumption, transcription of genes encoding viral recognition proteins, Ddx58 (RIG-I), Tlr3, and Ifih1 (MDA5), were up-regulated. CD8⁺ T cells from IAV-infected mice expanded in response to IAV PB1 peptide-pulsed eosinophils, and CpG methylation in the Tbx21 promoter was reduced in these T cells. These data offer insight into how eosinophils respond to IAV and help elucidate alternative mechanisms by which they regulate antiviral immune responses during IAV infection.

Integrating Peak Colocalization and Motif Enrichment Analysis for the Discovery of Genome-Wide Regulatory Modules and Transcription Factor Recruitment Rules

Chromatin immunoprecipitation followed by next-generation sequencing (ChIP-Seq) has opened new avenues of research in the genome-wide characterization of regulatory DNA-protein interactions at the genetic and epigenetic level. As a consequence, it has become the de facto standard for studies on the regulation of transcription, and literally thousands of data sets for transcription factors and cofactors in different conditions and species are now available to the scientific community. However, while pipelines and best practices have been established for the analysis of a single experiment, there is still no consensus on the best way to perform an integrated analysis of multiple datasets in the same condition, in order to identify the most relevant and widespread regulatory modules composed by different transcription factors and cofactors. We present here a computational pipeline for this task, that integrates peak summit colocalization, a novel statistical framework for the evaluation of its significance, and motif enrichment analysis. We show examples of its application to ENCODE data, that led to the identification of relevant regulatory modules composed of different factors, as well as the organization on DNA of the binding motifs responsible for their recruitment.

NF-YA Overexpression in Lung Cancer: LUAD

The trimeric transcription factor (TF) NF-Y regulates the CCAAT box, a DNA element enriched in promoters of genes overexpressed in many types of cancer. The regulatory NF-YA is present in two major isoforms, NF-YAl ("long") and NF-YAs ("short"). There is growing indication that NF-YA levels are increased in tumors. Here, we report interrogation of RNA-Seq TCGA (The Cancer Genome Atlas)-all 576 samples-and GEO (Gene Expression Ominibus) datasets of lung adenocarcinoma (LUAD). NF-YAs is overexpressed in the three subtypes, proliferative, inflammatory, and TRU (terminal respiratory unit). CCAAT is enriched in promoters of tumor differently expressed genes (DEG) and in the proliferative/inflammatory intersection, matching with KEGG (Kyoto Encyclopedia of Genes and Genomes) terms cell-cycle and signaling. Increasing levels of NF-YAs are observed from low to high CpG island methylator phenotypes (CIMP). We identified 166 genes overexpressed in LUAD cell lines with low NF-YAs/NF-YAl ratios: applying this centroid to TCGA samples faithfully predicted tumors' isoform ratio. This signature lacks CCAAT in promoters. Finally, progression-free intervals and hazard ratios concurred with the worst prognosis of patients with either a low or high NF-YAs/NF-YAl ratio. In conclusion, global overexpression of NF-YAs is documented in LUAD and is associated with aggressive tumor behavior; however, a similar prognosis is recorded in tumors with high levels of NF-YAl and overexpressed CCAAT-less genes.

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.

Genomic Characterization of the Zinc Transcriptional Regulatory Element Reveals Potential Functional Roles of ZNF658

The zinc transcriptional regulatory element (ZTRE) is a newly reported binding motif for human zinc finger protein ZNF658, which alters gene expression in response to cellular zinc. The ZTRE has two nucleotide components-the palindromic flanking pairs and the bridging "N" bases between these flanks that range in number from 0 to 100. There are 12 pairs of ZTRE flanks (designated A-L). Three thousand five hundred twenty-five genes contain one or more ZTREs - 1000 to + 200 bp from their transcriptional start site (TSS). ZTRE-E is observed at a greater frequency, and ZTRE containing 25 bridging bases are less frequent, within - 200 bp from the TSS. The genes with ZTREs in this range are enriched in processes that may compensate zinc deficiency, while other genes with ZTREs outside this range are enriched in transcriptional activation processes. The division of ZTREs into two groups may imply a dual role of ZNF658, similar to the homologous yeast protein Zap1, via binding to low or high affinity sequences dependent upon cellular zinc. The KLF/Sp1-family binding motif is prevalent within the ZTRE "N" bridging bases, suggesting ZNF658 may compete with Sp1-like transactivators to suppress transcription.

NF-YA Overexpression in Lung Cancer: LUSC

The CCAAT box is recognized by the trimeric transcription factor NF-Y, whose NF-YA subunit is present in two major splicing isoforms, NF-YAl (“long”) and NF-YAs (“short”). Little is known about the expression levels of NF-Y subunits in tumors, and nothing in lung cancer. By interrogating RNA-seq TCGA and GEO datasets, we found that, unlike NF-YB/NF-YC, NF-YAs is overexpressed in lung squamous cell carcinomas (LUSC). The ratio of the two isoforms changes from normal to cancer cells, with NF-YAs becoming predominant in the latter. NF-YA increased expression correlates with common proliferation markers. We partitioned all 501 TCGA LUSC tumors in the four molecular cohorts and verified that NF-YAs is similarly overexpressed. We analyzed global and subtype-specific RNA-seq data and found that CCAAT is the most abundant DNA matrix in promoters of genes overexpressed in all subtypes. Enriched Gene Ontology terms are cell-cycle and signaling. Survival curves indicate a worse clinical outcome for patients with increasing global amounts of NF-YA; same with hazard ratios with very high and, surprisingly, very low NF-YAs/NF-YAl ratios. We then analyzed gene expression in this latter cohort and identified a different, pro-migration signature devoid of CCAAT. We conclude that overexpression of the NF-Y regulatory subunit in LUSC has the scope of increasing CCAAT-dependent, proliferative (NF-YAs^high) or CCAAT-less, pro-migration (NF-YAl^high) genes. The data further reinstate the importance of analysis of single isoforms of TFs involved in tumor development.

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.

Coordinated Transcriptional Regulation of Cytochrome P450 3As by Nuclear Transcription Factor Y and Specificity Protein 1

The cytochrome P450 3A subfamily plays vital roles in the metabolism of endogenous chemicals and xenobiotics. Understanding the basal expression of CYP3A in humans and pigs is crucial for drug evaluation. In this study, we demonstrated that the basal transcriptional regulation of CYP3A genes in hepatocytes is evolutionarily conserved between humans and pigs. The basal expression of CYP3A genes is transactivated by two cis-acting elements, the CCAAT and GC boxes, located a constant distance apart in the proximal promoter region of six CYP3A genes. Mutation analysis of these two cis-acting elements suggested that they play important roles in mediating basal expression, but to different extents because of the nucleotide variations in the elements. Two transcription factors, nuclear transcription factor Y (NF-Y) and specificity protein 1 (Sp1), directly bind to these cis-acting elements in CYP3A proximal promoters in HepG2 cells and porcine hepatocytes. Furthermore, changing the distance between the NF-Y and Sp1 binding sites resulted in decreases in the promoter activity of CYP3A genes. Conclusively, our results show that human and porcine CYP3A genes are regulated by NF-Y and Sp1 in a coordinated manner, and that the distance between these two cis-acting elements is crucial for constitutive CYP3A expression.

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.

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.

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.

SWIM: a computational tool to unveiling crucial nodes in complex biological networks

SWItchMiner (SWIM) is a wizard-like software implementation of a procedure, previously described, able to extract information contained in complex networks. Specifically, SWIM allows unearthing the existence of a new class of hubs, called “fight-club hubs”, characterized by a marked negative correlation with their first nearest neighbors. Among them, a special subset of genes, called “switch genes”, appears to be characterized by an unusual pattern of intra- and inter-module connections that confers them a crucial topological role, interestingly mirrored by the evidence of their clinic-biological relevance. Here, we applied SWIM to a large panel of cancer datasets from The Cancer Genome Atlas, in order to highlight switch genes that could be critically associated with the drastic changes in the physiological state of cells or tissues induced by the cancer development. We discovered that switch genes are found in all cancers we studied and they encompass protein coding genes and non-coding RNAs, recovering many known key cancer players but also many new potential biomarkers not yet characterized in cancer context. Furthermore, SWIM is amenable to detect switch genes in different organisms and cell conditions, with the potential to uncover important players in biologically relevant scenarios, including but not limited to human cancer.

NF-Y in cancer: Impact on cell transformation of a gene essential for proliferation

NF-Y is a ubiquitous heterotrimeric transcription factor with a binding affinity for the CCAAT consensus motif, one of the most common cis-acting element in the promoter and enhancer regions of eukaryote genes in direct (CCAAT) or reverse (ATTGG) orientation. NF-Y consists of three subunits, NF-YA, the regulatory subunit of the trimer, NF-YB, and NF-YC, all required for CCAAT binding. Growing evidence in cells and animal models support the notion that NF-Y, driving transcription of a plethora of cell cycle regulatory genes, is a key player in the regulation of proliferation. Proper control of cellular growth is critical for cancer prevention and uncontrolled proliferation is a hallmark of cancer cells. Indeed, during cell transformation aberrant molecular pathways disrupt mechanisms controlling proliferation and many growth regulatory genes are altered in tumors. Here, we review bioinformatics, molecular and functional evidence indicating the involvement of the cell cycle regulator NF-Y in cancer-associated pathways. This article is part of a Special Issue entitled: Nuclear Factor Y in Development and Disease, edited by Prof. Roberto Mantovani.

Lotus japonicus NF-YA1 Plays an Essential Role During Nodule Differentiation and Targets Members of the SHI/STY Gene Family

Legume plants engage in intimate relationships with rhizobial bacteria to form nitrogen-fixing nodules, root-derived organs that accommodate the microsymbiont. Members of the Nuclear Factor Y (NF-Y) gene family, which have undergone significant expansion and functional diversification during plant evolution, are essential for this symbiotic liaison. Acting in a partially redundant manner, NF-Y proteins were shown, previously, to regulate bacterial infection, including selection of a superior rhizobial strain, and to mediate nodule structure formation. However, the exact mechanism by which these transcriptional factors exert their symbiotic functions has remained elusive. By carrying out detailed functional analyses of Lotus japonicus mutants, we demonstrate that LjNF-YA1 becomes indispensable downstream from the initial cortical cell divisions but prior to nodule differentiation, including cell enlargement and vascular bundle formation. Three affiliates of the SHORT INTERNODES/STYLISH transcription factor gene family, called STY1, STY2, and STY3, are demonstrated to be among likely direct targets of LjNF-YA1, and our results point to their involvement in nodule formation.

NF-Y and the immune response: Dissecting the complex regulation of MHC genes

Nuclear Factor Y (NF-Y) was first described as one of the CCAAT binding factors. Although CCAAT motifs were found to be present in various genes, NF-Y attracted a lot of interest early on, due to its role in Major Histocompatibility Complex (MHC) gene regulation. MHC genes are crucial in immune response and show peculiar expression patterns. Among other conserved elements on MHC promoters, an NF-Y binding CCAAT box was found to contribute to MHC transcriptional regulation. NF-Y along with other DNA binding factors assembles in a stereospecific manner to form a multiprotein scaffold, the MHC enhanceosome, which is necessary but not sufficient to drive transcription. Transcriptional activation is achieved by the recruitment of yet another factor, the class II transcriptional activator (CIITA). In this review, we briefly discuss basic findings on MHCII transcription regulation and we highlight NF-Y different modes of function in MHCII gene activation. 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.

Structural determinants for NF-Y/DNA interaction at the CCAAT box

The recently determined crystal structures of the sequence-specific transcription factor NF-Y have illuminated the structural mechanism underlying transcription at the CCAAT box. NF-Y is a trimeric protein complex composed by the NF-YA, NF-YB, and NF-YC subunits. NF-YB and NF-YC contain a histone-like domain and assemble on a head-to-tail fashion to form a dimer, which provides the structural scaffold for the DNA sugar-phosphate backbone binding (mimicking the nucleosome H2A/H2B-DNA assembly) and for the interaction with NF-YA. The NF-YA subunit hosts two structurally extended α-helices; one is involved in NF-YB/NF-YC binding and the other inserts deeply into the DNA minor groove, providing exquisite sequence-specificity for recognition and binding of the CCAAT box. The analysis of these structural data is expected to serve as a powerful guide for future experiments aimed at understanding the role of post-translational modification at NF-Y regulation sites and to unravel the three-dimensional architecture of higher order complexes formed between NF-Y and other transcription factors that act synergistically for transcription activation. Moreover, these structures represent an excellent starting point to challenge the formation of a stable hybrid nucleosome between NF-Y and core histone proteins, and to rationalize the fine molecular details associated with the wide combinatorial association of plant NF-Y subunits. This article is part of a Special Issue entitled: Nuclear Factor Y in Development and Disease, edited by Prof. Roberto Mantovani.

NF-Y Binding Site Architecture Defines a C-Fos Targeted Promoter Class

ChIP-seq experiments detect the chromatin occupancy of known transcription factors in a genome-wide fashion. The comparisons of several species-specific ChIP-seq libraries done for different transcription factors have revealed a complex combinatorial and context-specific co-localization behavior for the identified binding regions. In this study we have investigated human derived ChIP-seq data to identify common cis-regulatory principles for the human transcription factor c-Fos. We found that in four different cell lines, c-Fos targeted proximal and distal genomic intervals show prevalences for either AP-1 motifs or CCAAT boxes as known binding motifs for the transcription factor NF-Y, and thereby act in a mutually exclusive manner. For proximal regions of co-localized c-Fos and NF-YB binding, we gathered evidence that a characteristic configuration of repeating CCAAT motifs may be responsible for attracting c-Fos, probably provided by a nearby AP-1 bound enhancer. Our results suggest a novel regulatory function of NF-Y in gene-proximal regions. Specific CCAAT dimer repeats bound by the transcription factor NF-Y define this novel cis-regulatory module. Based on this behavior we propose a new enhancer promoter interaction model based on AP-1 motif defined enhancers which interact with CCAAT-box characterized promoter regions.

The Arabidopsis thaliana Nuclear Factor Y Transcription Factors

Nuclear factor Y (NF-Y) is an evolutionarily conserved trimeric transcription factor complex present in nearly all eukaryotes. The heterotrimeric NF-Y complex consists of three subunits, NF-YA, NF-YB, and NF-YC, and binds to the CCAAT box in the promoter regions of its target genes to regulate their expression. Yeast and mammal genomes generally have single genes with multiple splicing isoforms that encode each NF-Y subunit. By contrast, plant genomes generally have multi-gene families encoding each subunit and these genes are differentially expressed in various tissues or stages. Therefore, different subunit combinations can lead to a wide variety of NF-Y complexes in various tissues, stages, and growth conditions, indicating the potentially diverse functions of this complex in plants. Indeed, many recent studies have proved that the NF-Y complex plays multiple essential roles in plant growth, development, and stress responses. In this review, we highlight recent progress on NF-Y in Arabidopsis thaliana, including NF-Y protein structure, heterotrimeric complex formation, and the molecular mechanism by which NF-Y regulates downstream target gene expression. We then focus on its biological functions and underlying molecular mechanisms. Finally, possible directions for future research on NF-Y are also presented.