The E2F6 protein functions as an Rb-independent repressor of gene transcription. We have previously provided evidence suggesting a role for E2F6 in repression of E2F-responsive genes at S phase. Here, we have identified BRG1, the ATPase subunit of the SWI/SNF chromatin-remodeling complex, as an E2F6 interacting protein. Immunoprecipitation experiments demonstrate that BRG1 binds specifically to E2F6 and E2F4 but not the activator E2Fs. E2F6 was also able to interact with BAF155, a BRG1-associated factor, in the SWI/SNF complex. Chromatin immunoprecipitation assays demonstrate the binding of BRG1 coincident with E2F6 on G1/S gene promoters during S phase. Collectively, our studies suggest that E2F6 may recruit BRG1 in transcriptional regulation of genes important for G1/S phase transition of the cell cycle.

In mouse development, long-term silencing by CpG island DNA methylation is specifically targeted to germline genes; however, the molecular mechanisms of this specificity remain unclear. Here, we demonstrate that the transcription factor E2F6, a member of the polycomb repressive complex 1.6 (PRC1.6), is critical to target and initiate epigenetic silencing at germline genes in early embryogenesis. Genome-wide, E2F6 binds preferentially to CpG islands in embryonic cells. E2F6 cooperates with MGA to silence a subgroup of germline genes in mouse embryonic stem cells and in embryos, a function that critically depends on the E2F6 marked box domain. Inactivation of E2f6 leads to a failure to deposit CpG island DNA methylation at these genes during implantation. Furthermore, E2F6 is required to initiate epigenetic silencing in early embryonic cells but becomes dispensable for the maintenance in differentiated cells. Our findings elucidate the mechanisms of epigenetic targeting of germline genes and provide a paradigm for how transient repression signals by DNA-binding factors in early embryonic cells are translated into long-term epigenetic silencing during mouse development.

E2F6 is a member of the E2F family of transcription factors involved in regulation of a wide variety of genes through both activation and repression. E2F6 has been reported as overexpressed in breast cancers but whether or not this is important for tumor development is unclear. We first checked E2F6 expression in tumor cDNAs and the protein level in a range of breast cancer cell lines. RNA interference-mediated depletion was then used to assess the importance of E2F6 expression in cell lines with regard to cell cycle profile using fluorescence-activated cell sorting and a cell survival assay using (3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT). The overexpression of E2F6 was confirmed in breast tumor cDNA samples and breast cancer cell lines. Depletion of E2F6 in the breast cancer cells reduced cell viability in MCF-7, T-47D, and MDA-MB-231 cells. There was little effect in the nontumor breast cell line MCF-10A. The deleterious effect on cancer cells was greater during replication stress, leading to an increase in the proportion of breast cancer cells with sub-G1 DNA content. These results suggest that E2F6 might be essential for the survival of breast cancer cells experiencing replication stress, and therefore it could be a target for combined therapy.

Glioma is the most common primary tumour of the central nervous system and is considered one of the greatest challenges for neurosurgery. Mounting evidence has shown that lncRNAs participate in various biological processes of tumours, including glioma. This study aimed to reveal the role and relevant mechanism of COX10-AS1 in glioma. The expression of COX10-AS1, miR-641 and E2F6 was measured by qRT-PCR and/or western blot. Clone formation assays, EdU assays, Transwell assays and tumour xenograft experiments were performed to evaluate the effects of COX10-AS1, miR-641 and E2F6 on glioma proliferation, migration and invasion. Luciferase reporter assays, RNA pull-down assays and ChIP assays were conducted to analyse the relationship among COX10-AS1, miR-641 and E2F6. We demonstrated that COX10-AS1 was upregulated in glioma tissues and cell lines, which was related to the grade of glioma and patient survival. Next, through functional assays, we found that COX10-AS1 influenced the proliferation, migration and invasion of glioma cell lines. Then, with the help of bioinformatics analysis, we confirmed that COX10-AS1 regulated glioma progress by acting as a sponge of miR-641 to regulate E2F6. Moreover, further study indicated that E2F6 could promote COX10-AS1 expression by binding to its promoter region. Taken together, the data indicated that COX10-AS1 acts as an oncogene in combination with COX10-AS1/miR-641/E2F6 in glioma, which may be beneficial to the diagnosis and treatment of glioma.

Ovarian cancer is the most lethal cancer of the female reproductive system. In that regard, several epidemiological studies suggest that long-term exposure to estrogen could increase ovarian cancer risk, although its precise role remains controversial. To decipher a mechanism for this, we previously generated a mathematical model of how estrogen-mediated upregulation of the transcription factor, E2F6, upregulates the ovarian cancer stem/initiating cell marker, c-Kit, by epigenetic silencing the tumor suppressor miR-193a, and a competing endogenous (ceRNA) mechanism. In this study, we tested that previous mathematical model, showing that estrogen treatment of immortalized ovarian surface epithelial cells upregulated both E2F6 and c-KIT, but downregulated miR-193a. Luciferase assays further confirmed that microRNA-193a targets both E2F6 and c-Kit. Interestingly, ChIP-PCR and bisulphite pyrosequencing showed that E2F6 also epigenetically suppresses miR-193a, through recruitment of EZH2, and by a complex ceRNA mechanism in ovarian cancer cell lines. Importantly, cell line and animal experiments both confirmed that E2F6 promotes ovarian cancer stemness, whereas E2F6 or EZH2 depletion derepressed miR-193a, which opposes cancer stemness, by alleviating DNA methylation and repressive chromatin. Finally, 118 ovarian cancer patients with miR-193a promoter hypermethylation had poorer survival than those without hypermethylation. These results suggest that an estrogen-mediated E2F6 ceRNA network epigenetically and competitively inhibits microRNA-193a activity, promoting ovarian cancer stemness and tumorigenesis.

BACKGROUND Pancreatic adenocarcinoma (PAAD) is one of the deadliest types of cancer. In the early stages, patients often have atypical symptoms, making diagnosis difficult. The prognosis of diagnosed patients is very poor and treating PAAD is challenging. Therefore, determining reliable risk factors related to PAAD development is critical for improving patient prognosis. E2F family transcription factors (TFs) are essential regulators of DNA synthesis and cell cycle progression in eukaryotic cells, and they have been identified as prognostic biomarkers associated with multiple cancer types. However, further research is necessary to establish the prognostic relevance of these TFs in PAAD patients. MATERIAL AND METHODS We assessed PAAD patient transcriptional and outcome data using the TIMER, ONCOMINE, STRING, GEPIA, cBioPortal, Kaplan-Meier Plotter, GSCALite, and starBase databases. RESULTS PAAD tumor tissues exhibited increased expression of E2F1/3/5/7/8 relative to that in normal tissues, while the expression of E2F2/3/6/8 was associated with a more advanced tumor stage. Survival analyses indicated that PAAD patients expressing higher levels of E2F1/2/3/7/8 exhibited shorter overall survival (OS) and disease-free survival (DFS) than patients expressing lower levels of these TFs. In addition, E2F4 and E2F6 overexpression was associated with poorer DFS and OS, respectively. We also found that the expression of E2Fs was significantly correlated with immune infiltrates, including CD8+ T cells, CD4+ T cells, B cells, dendritic cells, neutrophils, and macrophages. CONCLUSIONS Our study may provide new insights into the optimal choice of immunotherapy and promising novel targets for therapeutic intervention in PAAD patients.

Long noncoding RNA LAMTOR5 antisense RNA 1 (LAMTOR5-AS1) has been certified as a risk predictor and diagnostic biomarker of prostate cancer. However, the expression and exact roles of LAMTOR5-AS1 in non-small cell lung cancer (NSCLC) remain unclear. Thus, we measured LAMTOR5-AS1 expression in NSCLC and gauged its clinical value. The detailed roles and downstream working mechanism of LAMTOR5-AS1 in NSCLC were comprehensively unraveled. qRT-PCR was applied to measure gene expression. Functionally, utilizing small interfering RNA, LAMTOR5-AS1 was ablated, and the functional alterations were addressed by means of different experiments. The targeting activities between LAMTOR5-AS1 and microRNA-506-3p (miR-506-3p) and between miR-506-3p and E2F transcription factor 6 (E2F6) were confirmed by RNA immunoprecipitation and luciferase reporter assays. LAMTOR5-AS1 overexpression in NSCLC was verified in TCGA datasets and our own cohort and manifested an evident relationship with poor prognosis. Interference with LAMTOR5-AS1 led to repression of the proliferation, cloning, and metastasis abilities of NSCLC cells in vitro. We further confirmed an obvious increase in LAMTOR5-AS1-silenced NSCLC cell apoptosis. Furthermore, the absence of LAMTOR5-AS1 restricted tumor growth in vivo. Mechanistically, LAMTOR5-AS1 sponged miR-506-3p in NSCLC cells. Furthermore, E2F6, a downstream target of miR-506-3p, was under the control of LAMTOR5-AS1, which was realized by decoying miR-506-3p. Rescue experiments showed that miR-506-3p suppression or E2F6 reintroduction was capable of remitting LAMTOR5-AS1 deficiency-triggered anticarcinogenic actions in NSCLC. Our study confirmed the exact roles of LAMTOR5-AS1 for the first time and revealed that LAMTOR5-AS1 knockdown disrupts the malignancy of NSCLC by targeting the miR-506-3p/E2F6 axis. Targeting the LAMTOR5-AS1/miR-506-3p/E2F6 pathway may be instrumental for managing patients with NSCLC.

Long noncoding RNA SLC9A3 antisense RNA 1 (SLC9A3‑AS1) plays a central role in lung cancer; yet, its functions in nasopharyngeal carcinoma (NPC) have not been elucidated. The present study revealed the roles of SLC9A3‑AS1 in NPC and dissected the mechanisms downstream of SLC9A3‑AS1. SLC9A3‑AS1 levels in NPC were assessed by applying RT‑qPCR. The modulatory role of SLC9A3‑AS1 interference on NPC cells was examined using numerous functional experiments. High expression of SLC9A3‑AS1 was observed in NPC samples. Patients with NPC with a high level of SLC9A3‑AS1 experienced a shorter overall survival than those with a low SLC9A3‑AS1 level. Loss of SLC9A3‑AS1 reduced NPC cell proliferation, colony formation, migration, and invasion but induced cell apoptosis in vitro. Animal experiments further revealed that the depletion of SLC9A3‑AS1 hindered NPC tumour growth in vivo. As a competitive endogenous RNA, SLC9A3‑AS1 sponged microRNA‑486‑5p (miR‑486‑5p), consequently upregulating E2F transcription factor 6 (E2F6). Finally, the effects of SLC9A3‑AS1 silencing on NPC cells were reversed by inhibiting miR‑486‑5p or overexpressing E2F6. In summary, SLC9A3‑AS1 exerted carcinogenic effects on NPC cells by adjusting the miR‑486‑5p/E2F6 axis. Accordingly, the newly identified SLC9A3‑AS1/miR‑486‑5p/E2F6 pathway may offer attractive therapeutic targets for future development.

Dysregulation of long noncoding RNA (lncRNA) is known to be involved in numerous human diseases, including lung cancer. However, the precise biological functions of most lncRNA remain to be elucidated. Here, we identified a novel up-regulated lncRNA, LINC01436 (RefSeq: NR_110419.1), in non-small cell lung cancer (NSCLC). High expression of LINC01436 was significantly associated with poor overall survival. Notably, LINC01436 expression was transcriptionally repressed by E2F6 under normoxia, and the inhibitory effect was relieved in a hypoxic microenvironment. Gain- and loss-of-function studies revealed that LINC01436 acted as a proto-oncogene by promoting lung cancer cell growth, migration and invasion in vitro. Xenograft tumor assays in nude mice confirmed that LINC01436 promoted tumor growth and metastasis in vivo. Mechanistically, LINC01436 exerted biological functions by acting as a microRNA (miR)-30a-3p sponge to regulate the expression of its target gene EPAS1. Our findings characterize LINC01436 as a new hypoxia-sensitive lncRNA with oncogenic function in NSCLC, suggesting that LINC01436 may be a potential biomarker for prognosis and a potential target for treatment.

Loss of expression of surface antigens represents a significant problem for cancer immunotherapy. Microphthalmia-associated transcription factor (MITF-M) regulates melanocyte fate by driving expression of many differentiation genes, whose protein products can be recognized by cytolytic T lymphocytes. We previously reported that interleukin-1ß (IL-1ß) can downregulate MITF-M levels. Here we show that downregulation of MITF-M expression by IL-1ß was paralleled by an upregulation of miR-155 expression in four melanoma lines. We confirmed that miR-155 was able to target endogenous MITF-M in melanoma cells and demonstrated a role for miR-155 in the IL-1ß-induced repression of MITF-M by using an antagomiR. Notably, we also observed a strong negative correlation between MITF-M and miR-155 levels in a mouse model of melanoma. Taken together, our results indicate that MITF-M downregulation by inflammatory stimuli might be partly due to miR-155 upregulation. This could represent a novel mechanism of melanoma immune escape in an inflammatory microenvironment.

Transcription factor TFII-I is a multifunctional protein implicated in the regulation of cell cycle and stress-response genes. Previous studies have shown that a subset of TFII-I associated genomic sites contained DNA-binding motifs for E2F family transcription factors. We analyzed the co-association of TFII-I and E2Fs in more detail using bioinformatics, chromatin immunoprecipitation, and co-immunoprecipitation experiments. The data show that TFII-I interacts with E2F transcription factors. Furthermore, TFII-I, E2F4, and E2F6 interact with DNA-regulatory elements of several genes implicated in the regulation of the cell cycle, including DNMT1, HDAC1, CDKN1C, and CDC27. Inhibition of TFII-I expression led to a decrease in gene expression and in the association of E2F4 and E2F6 with these gene loci in human erythroleukemia K562 cells. Finally, TFII-I deficiency reduced the proliferation of K562 cells and increased the sensitivity toward doxorubicin toxicity. The results uncover novel interactions between TFII-I and E2Fs and suggest that TFII-I mediates E2F function at specific cell cycle genes.

E2Fs, as a family of pivotal transcription factors, have been implicated in multiple biological functions in human cancer; however, the expression and prognostic significance of E2Fs in breast cancer remains unknown. In the present study, the mRNA expression patterns of E2Fs in breast cancer were investigated with Oncomine and The Cancer Genome Atlas data. Prognostic values of E2Fs for patients with breast cancer were determined using the Kaplan-Meier plotter database. The results strongly indicated that E2F1, E2F2, E2F3, E2F5, E2F7 and E2F8 were overexpressed in patients with breast cancer, whereas E2F4 and E2F6 exhibited no expression difference between patients with cancer and healthy controls. In survival analyses, elevated E2F1, E2F3, E2F5, E2F7 and E2F8 expression levels were significantly associated with lower overall survival, relapse-free survival (RFS), distant metastasis-free survival (DMFS) or post-progression survival for patients with breast cancer. Furthermore, high expression of E2F4 indicated improved RFS but reduced DMFS. Subgroup analyses based on four clinicopathological factors further revealed that E2Fs were associated with the prognosis of patients with breast cancer in an estrogen receptor-, progesterone receptor-, human epidermal growth factor 2- and lymph node status-specific manner. These data indicated that E2Fs may serve as promising biomarkers and therapeutic targets for breast cancer.

Ctf18-replication factor C complex including Dscc1 (DNA replication and sister chromatid cohesion 1) is implicated in sister chromatid cohesion, DNA replication, and genome stability in S. cerevisiae and C. elegans. We previously performed gene expression profiling in primary colorectal cancer cells in order to identify novel molecular targets for the treatment of colorectal cancer. A feature of the cancer-associated transcriptional signature revealed from this effort is the elevated expression of the proto-oncogene DSCC1. Here, we have interrogated the molecular basis for deviant expression of human DSCC1 in colorectal cancer and its ability to promote survival of cancer cells. Quantitative PCR and immunohistochemical analyses corroborated that the expression level of DSCC1 is elevated in 60-70% of colorectal tumors compared to their matched noncancerous colonic mucosa. An in silico evaluation of the presumptive DSCC1 promoter region for consensus DNA transcriptional regulatory elements revealed a potential role for the E2F family of DNA-binding proteins in controlling DSCC1 expression. RNAi-mediated reduction of E2F1 reduced expression of DSCC1 in colorectal cancer cells. Gain- and loss-of-function experiments demonstrated that DSCC1 is involved in the viability of cancer cells in response to genotoxic stimuli. We reveal that E2F-dependent expression of DSCC1 confers anti-apoptotic properties in colorectal cancer cells, and that its suppression may be a useful option for the treatment of colorectal cancer.

We identify the cell cycle-regulated mRNA transcripts genome-wide in the osteosarcoma-derived U2OS cell line. This results in 2140 transcripts mapping to 1871 unique cell cycle-regulated genes that show periodic oscillations across multiple synchronous cell cycles. We identify genomic loci bound by the G2/M transcription factor FOXM1 by chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq) and associate these with cell cycle-regulated genes. FOXM1 is bound to cell cycle-regulated genes with peak expression in both S phase and G2/M phases. We show that ChIP-seq genomic loci are responsive to FOXM1 using a real-time luciferase assay in live cells, showing that FOXM1 strongly activates promoters of G2/M phase genes and weakly activates those induced in S phase. Analysis of ChIP-seq data from a panel of cell cycle transcription factors (E2F1, E2F4, E2F6, and GABPA) from the Encyclopedia of DNA Elements and ChIP-seq data for the DREAM complex finds that a set of core cell cycle genes regulated in both U2OS and HeLa cells are bound by multiple cell cycle transcription factors. These data identify the cell cycle-regulated genes in a second cancer-derived cell line and provide a comprehensive picture of the transcriptional regulatory systems controlling periodic gene expression in the human cell division cycle.

Eukaryotic translation initiation factor 4E binding protein 1 (EIF4EBP1) encodes the 4EBP1 protein, a negative regulator of mRNA translation and a substrate of the mechanistic target of rapamycin (mTOR), whose function and relevance in cancer is still under debate. Here, we analyzed EIF4EBP1 expression in different glioma patient cohorts and investigated its mode of transcriptional regulation in glioblastoma cells. We verified that EIF4EBP1 mRNA is overexpressed in malignant gliomas, including isocitrate dehydrogenase (IDH)-wildtype glioblastomas, relative to non-neoplastic brain tissue in multiple publically available datasets. Our analyses revealed that EIF4EBP1 overexpression in malignant gliomas is neither due to gene amplification nor to altered DNA methylation, but rather results from aberrant transcriptional activation by distinct transcription factors. We found seven transcription factor candidates co-expressed with EIF4EBP1 in gliomas and bound to the EIF4EBP1 promoter, as revealed by chromatin immunoprecipitation (ChIP)-sequencing data. We investigated the ability of these candidates to activate the EIF4EBP1 promoter using luciferase reporter assays, which supported four transcription factors as candidate EIF4EBP1 regulators, namely MYBL2, ETS1, HIF-1A, and E2F6. Finally, by employing transient knock-down experiments to repress either of these transcription factors, we identified MYBL2 and ETS1 as the relevant transcriptional drivers of enhanced EIF4EBP1 expression in malignant glioma cells. Taken together, our findings confirm enhanced expression of EIF4EBP1 in malignant gliomas relative to non-neoplastic brain tissue and characterize the underlying molecular pathomechanisms.

Liver hepatocellular carcinoma (HCC) remains a leading cause of cancer-related death. Poor understanding of the mechanisms underlying HCC prevents early detection and leads to high mortality. We developed a random forest model that incorporates copy-number variation, DNA methylation, transcription factor, and microRNA binding information as features to predict gene expression in HCC. Our model achieved a highly significant correlation between predicted and measured expression of held-out genes. Furthermore, we identified potential regulators of gene expression in HCC. Many of these regulators have been previously found to be associated with cancer and are differentially expressed in HCC. We also evaluated our predicted target sets for these regulators by making comparison with experimental results. Lastly, we found that the transcription factor E2F6, one of the candidate regulators inferred by our model, is predictive of survival rate in HCC. Results of this study will provide directions for future prospective studies in HCC.

In eukaryotic genomes, about 10% of genes are arranged in a head-to-head (H2H) orientation, and the distance between the transcription start sites of each gene pair is closer than 1 kb. Two genes in an H2H pair are prone to co-express and co-function. There have been many studies on bidirectional promoters. However, the mechanism by which H2H genes are regulated at the transcriptional level still needs further clarification, especially with regard to the co-regulation of H2H pairs. In this study, we first used the Hi-C data of chromatin linkages to identify spatially interacting H2H pairs, and then integrated ChIP-seq data to compare H2H gene pairs with and without evidence of spatial interactions in terms of their binding transcription factors (TFs). Using ChIP-seq and DNase-seq data, histones and DNase associated with H2H pairs were identified. Furthermore, we looked into the connections between H2H genes in a human co-expression network.

Dysregulation of the normal gene expression program is the cause of a broad range of diseases, including cancer. Detecting the specific perturbed regulators that have an effect on the generation and the development of the disease is crucial for understanding the disease mechanism and for taking decisions on efficient preventive and curative therapies. Moreover, detecting such perturbations at the patient level is even more important from the perspective of personalized medicine. We applied the Transcription Factor Target Enrichment Analysis, a method that detects the activity of transcription factors based on the quantification of the collective transcriptional activation of their targets, to a large collection of 5607 cancer samples covering eleven cancer types. We produced for the first time a comprehensive catalogue of altered transcription factor activities in cancer, a considerable number of them significantly associated to patient's survival. Moreover, we described several interesting TFs whose activity do not change substantially in the cancer with respect to the normal tissue but ultimately play an important role in patient prognostic determination, which suggest they might be promising therapeutic targets. An additional advantage of this method is that it allows obtaining personalized TF activity estimations for individual patients.

We present a network framework for analyzing multi-level regulation in higher eukaryotes based on systematic integration of various high-throughput datasets. The network, namely the integrated regulatory network, consists of three major types of regulation: TF→gene, TF→miRNA and miRNA→gene. We identified the target genes and target miRNAs for a set of TFs based on the ChIP-Seq binding profiles, the predicted targets of miRNAs using annotated 3'UTR sequences and conservation information. Making use of the system-wide RNA-Seq profiles, we classified transcription factors into positive and negative regulators and assigned a sign for each regulatory interaction. Other types of edges such as protein-protein interactions and potential intra-regulations between miRNAs based on the embedding of miRNAs in their host genes were further incorporated. We examined the topological structures of the network, including its hierarchical organization and motif enrichment. We found that transcription factors downstream of the hierarchy distinguish themselves by expressing more uniformly at various tissues, have more interacting partners, and are more likely to be essential. We found an over-representation of notable network motifs, including a FFL in which a miRNA cost-effectively shuts down a transcription factor and its target. We used data of C. elegans from the modENCODE project as a primary model to illustrate our framework, but further verified the results using other two data sets. As more and more genome-wide ChIP-Seq and RNA-Seq data becomes available in the near future, our methods of data integration have various potential applications.

Lung cancer accounts for the highest number of cancer-related deaths worldwide. Early diagnosis significantly increases the disease-free survival rate and a large amount of effort has been expended in screening trials and the development of early molecular diagnostics. However, a gold standard diagnostic strategy is not yet available. Here, based on miRNA expression profile in lung cancer and using a novel in silico reverse-transcriptomics approach, followed by analysis of the interactome; we have identified potential transcription factor (TF) markers that would facilitate diagnosis of subtype specific lung cancer. A subset of seven TF markers has been used in a microarray screen and was then validated by blood-based qPCR using stage-II and IV non-small cell lung carcinomas (NSCLC). Our results suggest that overexpression of HMGA1, E2F6, IRF1, and TFDP1 and downregulation or no expression of SUV39H1, RBL1, and HNRPD in blood is suitable for diagnosis of lung adenocarcinoma and squamous cell carcinoma sub-types of NSCLC. Here, E2F6 was, for the first time, found to be upregulated in NSCLC blood samples. The miRNA-TF-miRNA interaction based molecular mechanisms of these seven markers in NSCLC revealed that HMGA1 and TFDP1 play vital roles in lung cancer tumorigenesis. The strategy developed in this work is applicable to any other cancer or disease and can assist in the identification of potential biomarkers.