Contextual Refinement of Regulatory Targets Reveals Effects on Breast Cancer Prognosis of the Regulome

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.

Regulator combinations identify systemic sclerosis patients with more severe disease

Systemic sclerosis (SSc) is a heterogeneous autoimmune disorder that results in skin fibrosis, autoantibody production, and internal organ dysfunction. We previously identified 4 “intrinsic” subsets of SSc based upon skin gene expression that are found across organ systems. Gene expression regulators that underlie the SSc-intrinsic subsets, or are associated with clinical covariates, have not been systematically characterized. Here, we present a computational framework to calculate the activity scores of gene expression regulators and identify their associations with SSc clinical outcomes. We found that regulator activity scores can reproduce the intrinsic molecular subsets, with distinct sets of regulators identified for inflammatory, fibroproliferative, limited, and normal-like samples. Regulators most highly correlated with modified Rodnan skin score (MRSS) also varied by intrinsic subset. We identified subgroups of patients with fibroproliferative and inflammatory SSc with more severe pathophenotypes, such as higher MRSS and increased likelihood of interstitial lung disease (ILD). Using an independent cohort, we show that the group with more severe ILD was more likely to show forced vital capacity decline over a period of 36–54 months. Our results demonstrate an association among the activation of regulators, gene expression subsets, and clinical variables that can identify patients with SSc with more severe disease.

An association between the activation of regulators, gene expression subsets, and clinical variables identifies systemic sclerosis patients with more severe disease.

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-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.

Exploration of prognosis-related microRNA and transcription factor co-regulatory networks across cancer types

The study of cancer prognosis serves as an important part of cancer research. Large-scale cancer studies have identified numerous genes and microRNAs (miRNAs) associated with prognosis. These informative genes and miRNAs represent potential biomarkers to predict survival and to elucidate the molecular mechanism of tumour progression. MiRNAs and transcription factors (TFs) can work cooperatively as essential mediators of gene expression, and their dysregulation affects cancer prognosis. A panoramic view of cancer prognosis at the system level, considering the co-regulation roles of miRNA and TF, remains elusive. Here, we establish 12 prognosis-related miRNA-TF co-regulatory networks. The characteristics of prognostic target genes and their regulators in the network are depicted. Although the target genes and co-regulatory patterns exhibit cancer-specific properties, some miRNAs and TFs are highly conserved across cancers. We illustrate and interpret the roles of these conserved regulators by building a model associated with cancer hallmarks, functional enrichment analysis, network community detection, and exhaustive literature research. The elaborated system-level prognostic miRNA-TF co-regulation landscape, including the highlighted roles of conserved regulators, provides a novel and powerful insights into further biological and medical discoveries.