Genomic binding of NF-Y in mouse and human cells.

NF-Y is a Transcription Factor that regulates transcription through binding to the CCAAT-box. To understand its strategy, we analyzed 16 ChIP-seq datasets from human and mouse cells. Shared loci, mostly located in promoters of expressed genes of cell cycle, metabolism and gene expression pathways, are associated with histone marks of active chromatin and specific modules of TFs. Other peaks are in enhancers and Transposable Elements -TE- of retroviral origin in human and mouse. We evaluated the relationship with USF1, a common synergistic partner in promoters and MLT1 TEs, upon NF-YB inactivation: USF1 binding decreases in promoters, modestly in MLT1, suggesting a pioneering role of NF-Y in formers, not in the latters. These data define a common set of NF-Y functional targets across different mammalian cell types, suggesting a pioneering role in promoters with respect to TEs.

The inhibition of inner mitochondrial fusion in hepatocytes reduces non-alcoholic fatty liver and improves metabolic profile during obesity by modulating bile acid conjugation.

AIMS: Mitochondria are plastic organelles that continuously undergo biogenesis, fusion, fission, and mitophagy to control cellular energy metabolism, calcium homeostasis, hormones, sterols, and bile acids (BAs) synthesis. Here, we evaluated how the impairment of mitochondrial fusion in hepatocytes affects diet-induced liver steatosis and obesity. METHODS AND RESULTS: Male mice selectively lacking the key protein involved in inner mitochondrial fusion, optic atrophy 1 (OPA1) (OPA1ΔHep) were fed a high fat diet (HFD) for 20 weeks. OPA1ΔHep mice were protected from the development of hepatic steatosis and obesity because of reduced lipid absorption; a profile which was accompanied by increased respiratory exchange ratio in vivo, suggesting a preference for carbohydrates in OPA1ΔHep compared to controls. At the molecular level, this phenotype emerged as a consequence of poor mitochondria-peroxisome- endoplasmic reticulum (ER) tethering in OPA1 deficient hepatocytes, which impaired BAs conjugation and release in the bile, thus impacting lipid absorption from the diet. Concordantly, the liver of subjects with non-alcoholic fatty liver disease (NAFLD) presented an increased expression of OPA1 and of the network of proteins involved in mitochondrial function when compared with controls. CONCLUSION: Patients with NAFLD present increased expression of proteins involved in mitochondrial fusion in the liver. The selective deficency of OPA1 in hepatocytes protects mice from HFD-induced metabolic dysfunction by reducing BAs secretion and dietary lipids absorption as a consequence of reduced liver mitochondria-peroxisome-ER tethering.

The NF-Y splicing signature controls hybrid EMT and ECM-related pathways to promote aggressiveness of colon cancer.

Aberrant splicing events are associated with colorectal cancer (CRC) and provide new opportunities for tumor diagnosis and treatment. The expression of the splice variants of NF-YA, the DNA binding subunit of the transcription factor NF-Y, is deregulated in multiple cancer types compared to healthy tissues. NF-YAs and NF-YAl isoforms differ in the transactivation domain, which may result in distinct transcriptional programs. In this study, we demonstrated that the NF-YAl transcript is higher in aggressive mesenchymal CRCs and predicts shorter patients' survival. In 2D and 3D conditions, CRC cells overexpressing NF-YAl (NF-YAlhigh) exhibit reduced cell proliferation, rapid single cell amoeboid-like migration, and form irregular spheroids with poor cell-to-cell adhesion. Compared to NF-YAshigh, NF-YAlhigh cells show changes in the transcription of genes involved in epithelial-mesenchymal transition, extracellular matrix and cell adhesion. NF-YAl and NF-YAs bind similarly to the promoter of the E-cadherin gene, but oppositely regulate its transcription. The increased metastatic potential of NF-YAlhigh cells in vivo was confirmed in zebrafish xenografts. These results suggest that the NF-YAl splice variant could be a new CRC prognostic factor and that splice-switching strategies may reduce metastatic CRC progression.

NF-YAl drives EMT in Claudinlow tumours.

NF-Y is a trimeric transcription factor whose binding site -the CCAAT box- is enriched in cancer-promoting genes. The regulatory subunit, the sequence-specificity conferring NF-YA, comes in two major isoforms, NF-YA long (NF-YAl) and short (NF-YAs). Extensive expression analysis in epithelial cancers determined two features: widespread overexpression and changes in NF-YAl/NF-YAs ratios (NF-YAr) in tumours with EMT features. We performed wet and in silico experiments to explore the role of the isoforms in breast -BRCA- and gastric -STAD- cancers. We generated clones of two Claudinlow BRCA lines SUM159PT and BT549 ablated of exon-3, thus shifting expression from NF-YAl to NF-YAs. Edited clones show normal growth but reduced migratory capacities in vitro and ability to metastatize in vivo. Using TCGA, including upon deconvolution of scRNA-seq data, we formalize the clinical importance of high NF-YAr, associated to EMT genes and cell populations. We derive a novel, prognostic 158 genes signature common to BRCA and STAD Claudinlow tumours. Finally, we identify splicing factors associated to high NF-YAr, validating RBFOX2 as promoting expression of NF-YAl. These data bring three relevant results: (i) the definition and clinical implications of NF-YAr and the 158 genes signature in Claudinlow tumours; (ii) genetic evidence of 28 amino acids in NF-YAl with EMT-promoting capacity; (iii) the definition of selected splicing factors associated to NF-YA isoforms.

NF-Y subunits overexpression in gastric adenocarcinomas (STAD).

NF-Y is a pioneer transcription factor-TF-formed by the Histone-like NF-YB/NF-YC subunits and the regulatory NF-YA. It binds to the CCAAT box, an element enriched in promoters of genes overexpressed in many types of cancer. NF-YA is present in two major isoforms-NF-YAs and NF-YAl-due to alternative splicing, overexpressed in epithelial tumors. Here we analyzed NF-Y expression in stomach adenocarcinomas (STAD). We completed the partitioning of all TCGA tumor samples (450) according to molecular subtypes proposed by TCGA and ACRG, using the deep learning tool DeepCC. We analyzed differentially expressed genes-DEG-for enriched pathways and TFs binding sites in promoters. CCAAT is the predominant element only in the core group of genes upregulated in all subtypes, with cell-cycle gene signatures. NF-Y subunits are overexpressed, particularly NF-YA. NF-YAs is predominant in CIN, MSI and EBV TCGA subtypes, NF-YAl is higher in GS and in the ACRG EMT subtypes. Moreover, NF-YAlhigh tumors correlate with a discrete Claudinlow cohort. Elevated NF-YB levels are protective in MSS;TP53+ patients, whereas high NF-YAl/NF-YAs ratios correlate with worse prognosis. We conclude that NF-Y isoforms are associated to clinically relevant features of gastric cancer.

Alternative splicing of NF-YA promotes prostate cancer aggressiveness and represents a new molecular marker for clinical stratification of patients.

BACKGROUND: Approaches based on expression signatures of prostate cancer (PCa) have been proposed to predict patient outcomes and response to treatments. The transcription factor NF-Y participates to the progression from benign epithelium to both localized and metastatic PCa and is associated with aggressive transcriptional profile. The gene encoding for NF-YA, the DNA-binding subunit of NF-Y, produces two alternatively spliced transcripts, NF-YAs and NF-YAl. Bioinformatic analyses pointed at NF-YA splicing as a key transcriptional signature to discriminate between different tumor molecular subtypes. In this study, we aimed to determine the pathophysiological role of NF-YA splice variants in PCa and their association with aggressive subtypes. METHODS: Data on the expression of NF-YA isoforms were extracted from the TCGA (The Cancer Genome Atlas) database of tumor prostate tissues and validated in prostate cell lines. Lentiviral transduction and CRISPR-Cas9 technology allowed the modulation of the expression of NF-YA splice variants in PCa cells. We characterized 3D cell cultures through in vitro assays and RNA-seq profilings. We used the rank-rank hypergeometric overlap approach to identify concordant/discordant gene expression signatures of NF-YAs/NF-YAl-overexpressing cells and human PCa patients. We performed in vivo studies in SHO-SCID mice to determine pathological and molecular phenotypes of NF-YAs/NF-YAl xenograft tumors. RESULTS: NF-YA depletion affects the tumorigenic potential of PCa cells in vitro and in vivo. Elevated NF-YAs levels are associated to aggressive PCa specimens, defined by Gleason Score and TNM classification. NF-YAl overexpression increases cell motility, while NF-YAs enhances cell proliferation in PCa 3D spheroids and xenograft tumors. The transcriptome of NF-YAs-spheroids has an extensive overlap with localized and metastatic human PCa signatures. According to PCa PAM50 classification, NF-YAs transcript levels are higher in LumB, characterized by poor prognosis compared to LumA and basal subtypes. A significant decrease in NF-YAs/NF-YAl ratio distinguishes PCa circulating tumor cells from cancer cells in metastatic sites, consistently with pro-migratory function of NF-YAl. Stratification of patients based on NF-YAs expression is predictive of clinical outcome. CONCLUSIONS: Altogether, our results indicate that the modulation of NF-YA isoforms affects prostate pathophysiological processes and contributes to cancer-relevant phenotype, in vitro and in vivo. Evaluation of NF-YA splicing may represent a new molecular strategy for risk assessment of PCa patients.

The transcription factor NF-Y participates to stem cell fate decision and regeneration in adult skeletal muscle.

The transcription factor NF-Y promotes cell proliferation and its activity often declines during differentiation through the regulation of NF-YA, the DNA binding subunit of the complex. In stem cell compartments, the shorter NF-YA splice variant is abundantly expressed and sustains their expansion. Here, we report that satellite cells, the stem cell population of adult skeletal muscle necessary for its growth and regeneration, express uniquely the longer NF-YA isoform, majorly associated with cell differentiation. Through the generation of a conditional knock out mouse model that selectively deletes the NF-YA gene in satellite cells, we demonstrate that NF-YA expression is fundamental to preserve the pool of muscle stem cells and ensures robust regenerative response to muscle injury. In vivo and ex vivo, satellite cells that survive to NF-YA loss exit the quiescence and are rapidly committed to early differentiation, despite delayed in the progression towards later states. In vitro results demonstrate that NF-YA-depleted muscle stem cells accumulate DNA damage and cannot properly differentiate. These data highlight a new scenario in stem cell biology for NF-Y activity, which is required for efficient myogenic differentiation.

NF-Y Subunits Overexpression in HNSCC.

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.

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.

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.

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.

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-YAshigh) or CCAAT-less, pro-migration (NF-YAlhigh) genes. The data further reinstate the importance of analysis of single isoforms of TFs involved in tumor development.