Cancer is a disease of subverted regulatory pathways. In this paper, we reconstruct the regulatory network around E2F, a family of transcription factors whose deregulation has been associated to cancer progression, chemoresistance, invasiveness, and metastasis. We integrate gene expression profiles of cancer cell lines from two E2F1-driven highly aggressive bladder and breast tumors, and use network analysis methods to identify the tumor type-specific core of the network. By combining logic-based network modeling, in vitro experimentation, and gene expression profiles from patient cohorts displaying tumor aggressiveness, we identify and experimentally validate distinctive, tumor type-specific signatures of receptor proteins associated to epithelial–mesenchymal transition in bladder and breast cancer. Our integrative network-based methodology, exemplified in the case of E2F1-induced aggressive tumors, has the potential to support the design of cohort- as well as tumor type-specific treatments and ultimately, to fight metastasis and therapy resistance.

Malignant melanoma is highly lethal due to its aggressive invasive properties and metastatic dissemination. The transcription factor E2F1 is crucial for melanoma progression through poorly understood mechanisms. Here, we show that the miR-224/miR-452 cluster is significantly increased in advanced melanoma and invasive/ metastatic cell lines that express high levels of E2F1. miR-224/miR-452 expression is directly activated by E2F1 through transactivation of the GABRE gene. Ectopic expression of miR-224/miR-452 in less aggressive cells induces EMT and cytoskeletal rearrangements and enhances migration/invasion. Conversely, miR-224/miR-452 depletion in metastatic cells induces the reversal of EMT, inhibition of motility, loss of the invasive phenotype and an absence of lung metastases in mice. We identify the metastasis suppressor TXNIP as new target of miR-224/miR-452 that induces feedback inhibition of E2F1 and show that miR-224/452-mediated downregulation of TXNIP is essential for E2F1-induced EMT and invasion. The E2F1-miR-224/452-TXNIP axis constitutes a molecular signature that predicts patient survival and may help to set novel therapies.

No abstract

Long non-coding RNAs (lncRNAs) have emerged as integral components of E2F1-regulated gene regulatory networks (GRNs), but their implication in advanced or treatment-refractory malignancy is unknown. Methods: We combined high-throughput transcriptomic approaches with bioinformatics and structure modeling to search for lncRNAs that participate in E2F1-activated prometastatic GRNs and their phenotypic targets in the highly-relevant case of E2F1-driven aggressive bladder cancer (BC). RNA immunoprecipitation was performed to verify RNA-protein interactions. Functional analyses including qRT-PCR, immunoblotting, luciferase assays and measurement of extracellular fluxes were conducted to validate expression and target gene regulation. Results: We identified E2F1-responsive lncRNA-SLC16A1-AS1 and its associated neighboring protein-coding gene, SLC16A1/MCT1, which both promote cancer invasiveness. Mechanistically, upon E2F1-mediated co-transactivation of the gene pair, SLC16A1-AS1 associates with E2F1 in a structure-dependent manner and forms an RNA-protein complex that enhances SLC16A1/MCT1 expression through binding to a composite SLC16A1-AS1:E2F1-responsive promoter element. Moreover, SLC16A1-AS1 increases aerobic glycolysis and mitochondrial respiration and fuels ATP production by fatty acid β-oxidation. These metabolic changes are accompanied by alterations in the expression of the SLC16A1-AS1:E2F1-responsive gene PPARA, a key mediator of fatty acid β-oxidation. Conclusions: Our results unveil a new gene regulatory program by which E2F1-induced lncRNA-SLC16A1-AS1 forms a complex with its transcription factor that promotes cancer metabolic reprogramming towards the acquisition of a hybrid oxidative phosphorylation/glycolysis cell phenotype favoring BC invasiveness.

High rates of lethal outcome in tumour metastasis are associated with the acquisition of invasiveness and chemoresistance. Several clinical studies indicate that E2F1 overexpression across high-grade tumours culminates in unfavourable prognosis and chemoresistance in patients. Thus, fine-tuning the expression of E2F1 could be a promising approach for treating patients showing chemoresistance.Methods: We integrated bioinformatics, structural and kinetic modelling, and experiments to study cooperative regulation of E2F1 by microRNA (miRNA) pairs in the context of anticancer chemotherapy resistance.Results: We showed that an enhanced E2F1 repression efficiency can be achieved in chemoresistant tumour cells through two cooperating miRNAs. Sequence and structural information were used to identify potential miRNA pairs that can form tertiary structures with E2F1 mRNA. We then employed molecular dynamics simulations to show that among the identified triplexes, miR-205-5p and miR-342-3p can form the most stable triplex with E2F1 mRNA. A mathematical model simulating the E2F1 regulation by the cooperative miRNAs predicted enhanced E2F1 repression, a feature that was verified by in vitro experiments. Finally, we integrated this cooperative miRNA regulation into a more comprehensive network to account for E2F1-related chemoresistance in tumour cells. The network model simulations and experimental data indicate the ability of enhanced expression of both miR-205-5p and miR-342-3p to decrease tumour chemoresistance by cooperatively repressing E2F1.Conclusions: Our results suggest that pairs of cooperating miRNAs could be used as potential RNA therapeutics to reduce E2F1-related chemoresistance.