Gene regulatory network construction identified NFYA as a diffuse subtype-specific prognostic factor in gastric cancer

The enhancer RNA ADCY10P1 is associated with the progression of ovarian cancer

BACKGROUND

Emerging evidence identifies enhancer RNAs (eRNAs) as a class of regulatory ncRNAs that can contribute to the transcription of target genes. In this study, we used an integrated data analysis method to identify the important role of eRNAs in ovarian cancer (OC).

METHODS

Gene expression profiles and clinical information from The Cancer Genome Atlas (TCGA) database were used for this study. Based on expression analysis using GEPIA2 gene and Kaplan-Meier survival was performed to ensure the significance of the selected enhancer RNA ADCY10P1 in OC. Next, we explored the correlation and clinical significance between ADCY10P1 and target gene NFYA. Furthermore, we evaluated the effects of overexpression of ADCY10P1 on the proliferation, migration, invasion and epithelial-mesenchymal transformation (EMT) of OC cell lines. We also investigated the biological function enrichment score of ADCY10P1 and verified it with OC cell lines. Finally, external validation was conducted, and the prognostic value of the ADCY10P1 in different tumors was demonstrated.

RESULTS

We selected the eRNA ADCY10P1 associated with OC prognosis, with NFYA as its predicted target gene. Low ADCY10P1 expression was found to be associated with poor overall survival, high histological grade, and advanced stage of OC. Additionally, overexpression of ADCY10P1 inhibited the proliferation, migration, invasion and EMT phenotype of OC cell lines. Furthermore, ADCY10P1 was observed to inhibit glycolysis and fatty acid metabolism, thereby affecting OC progression. Meanwhile, OC tissue samples were externally validated. In addition, the pan-cancer analysis revealed that ADCY10P1 had prognostic value in other cancers.

CONCLUSIONS

This study showed that ADCY10P1 plays a key role in OC progression and may facilitate prognosis prediction.

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.

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.

Gene networks and transcriptional regulators associated with liver cancer development and progression

BACKGROUND

Treatment options for hepatocellular carcinoma (HCC) are limited, and overall survival is poor. Despite the high frequency of this malignoma, its basic disease mechanisms are poorly understood. Therefore, the aim of this study was to use different methodological approaches and combine the results to improve our knowledge on the development and progression of HCC.

METHODS

Twenty-three HCC samples were characterized by histological, morphometric and cytogenetic analyses, as well as comparative genomic hybridization (aCGH) and genome-wide gene expression followed by a bioinformatic search for potential transcriptional regulators and master regulatory molecules of gene networks.

RESULTS

Histological evaluation revealed low, intermediate and high-grade HCCs, and gene expression analysis split them into two main sets: GE1-HCC and GE2-HCC, with a low and high proliferation gene expression signature, respectively. Array-based comparative genomic hybridization demonstrated a high level of chromosomal instability, with recurrent chromosomal gains of 1q, 6p, 7q, 8q, 11q, 17q, 19p/q and 20q in both HCC groups and losses of 1p, 4q, 6q, 13q and 18q characteristic for GE2-HCC. Gene expression and bioinformatics analyses revealed that different genes and gene regulatory networks underlie the distinct biological features observed in GE1-HCC and GE2-HCC. Besides previously reported dysregulated genes, the current study identified new candidate genes with a putative role in liver cancer, e.g. C1orf35, PAFAH1B3, ZNF219 and others.

CONCLUSION

Analysis of our findings, in accordance with the available published data, argues in favour of the notion that the activated E2F1 signalling pathway, which can be responsible for both inappropriate cell proliferation and initial chromosomal instability, plays a pivotal role in HCC development and progression. A dedifferentiation switch that manifests in exaggerated gene expression changes might be due to turning on transcriptional co-regulators with broad impact on gene expression, e.g. POU2F1 (OCT1) and NFY, as a response to accumulating cell stress during malignant development. Our findings point towards the necessity of different approaches for the treatment of HCC forms with low and high proliferation signatures and provide new candidates for developing appropriate HCC therapies.

Molecular Network Profiling in Intestinal- and Diffuse-Type Gastric Cancer

Epithelial-mesenchymal transition (EMT) plays an important role in the acquisition of cancer stem cell (CSC) feature and drug resistance, which are the main hallmarks of cancer malignancy. Although previous findings have shown that several signaling pathways are activated in cancer progression, the precise mechanism of signaling pathways in EMT and CSCs are not fully understood. In this study, we focused on the intestinal and diffuse-type gastric cancer (GC) and analyzed the gene expression of public RNAseq data to understand the molecular pathway regulation in different subtypes of gastric cancer. Network pathway analysis was performed by Ingenuity Pathway Analysis (IPA). A total of 2815 probe set IDs were significantly different between intestinal- and diffuse-type GC data in cBioPortal Cancer Genomics. Our analysis uncovered 10 genes including male-specific lethal 3 homolog (Drosophila) pseudogene 1 (MSL3P1), CDC28 protein kinase regulatory subunit 1B (CKS1B), DEAD-box helicase 27 (DDX27), golgi to ER traffic protein 4 (GET4), chromosome segregation 1 like (CSE1L), translocase of outer mitochondrial membrane 34 (TOMM34), YTH N6-methyladenosine RNA binding protein 1 (YTHDF1), ribonucleic acid export 1 (RAE1), par-6 family cell polarity regulator beta (PARD6B), and MRG domain binding protein (MRGBP), which have differences in gene expression between intestinal- and diffuse-type GC. A total of 463 direct relationships with three molecules (MYC, NTRK1, UBE2M) were found in the biomarker-filtered network generated by network pathway analysis. The networks and features in intestinal- and diffuse-type GC have been investigated and profiled in bioinformatics. Our results revealed the signaling pathway networks in intestinal- and diffuse-type GC, bringing new light for the elucidation of drug resistance mechanisms in CSCs.

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.

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.