Interaction of E2F7 transcription factor with E2F1 and C-terminal-binding protein (CtBP) provides a mechanism for E2F7-dependent transcription repression

Atlas of regulated target genes of transcription factors (ART-TF) in human ES cells

Background

Transcription factors (TFs) play central roles in maintaining “stemness” of embryonic stem (ES) cells and their differentiation into several hundreds of adult cell types. The regulatory competence of TFs is routinely assessed by detecting target genes to which they bind. However, these data do not indicate which target genes are activated, repressed, or not affected by the change of TF abundance. There is a lack of large-scale studies that compare the genome binding of TFs with the expression change of target genes after manipulation of each TF.

Results

In this paper we associated human TFs with their target genes by two criteria: binding to genes, evaluated from published ChIP-seq data (n = 1868); and change of target gene expression shortly after induction of each TF in human ES cells. Lists of direction- and strength-specific regulated target genes are generated for 311 TFs (out of 351 TFs tested) with expected proportion of false positives less than or equal to 0.30, including 63 new TFs not present in four existing databases of target genes. Our lists of direction-specific targets for 152 TFs (80.0%) are larger that in the TRRUST database. In average, 30.9% of genes that respond greater than or equal to twofold to the induction of TFs are regulated targets. Regulated target genes indicate that the majority of TFs are either strong activators or strong repressors, whereas sets of genes that responded greater than or equal to twofold to the induction of TFs did not show strong asymmetry in the direction of expression change. The majority of human TFs (82.1%) regulated their target genes primarily via binding to enhancers. Repression of target genes is more often mediated by promoter-binding than activation of target genes. Enhancer-promoter loops are more abundant among strong activator and repressor TFs.

Conclusions

We developed an atlas of regulated targets of TFs (ART-TF) in human ES cells by combining data on TF binding with data on gene expression change after manipulation of individual TFs. Sets of regulated gene targets were identified with a controlled rate of false positives. This approach contributes to the understanding of biological functions of TFs and organization of gene regulatory networks. This atlas should be a valuable resource for ES cell-based regenerative medicine studies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04924-3.

Circular RNAs in Alzheimer's Disease: A New Perspective of Diagnostic and Therapeutic Targets

BACKGROUND

Circular RNAs (circRNAs), as covalently closed single-stranded noncoding RNA molecules, have been recently identified to involve in several biological processes, principally through targeting microRNAs. Among various neurodegenerative diseases (NDs), accumulating evidence has proposed key roles for circRNAs in the pathogenesis of Alzheimer's disease (AD); although the exact relationship between these RNA molecules and AD progression is not clear, they have been believed to mostly act as miRNA sponges or gene transcription modulators through correlating with multiple proteins, involved in the accumulation of Amyloid β (Aβ) peptides, as well as tau protein, as AD's pathological hallmark. More interestingly, circRNAs have also been reported to play diagnostic and therapeutic roles during AD progression.

OBJECTIVE

Literature review indicated that circRNAs could essentially contribute to the onset and development of AD. Thus, in the current review, the circRNAs' biogenesis and functions are addressed at first, and then the interplay between particular circRNAs and AD is comprehensively discussed. Eventually, the diagnostic and therapeutic significance of these noncoding RNAs is highlighted in brief.

RESULTS

A large number of circRNAs are expressed in the brain. Thereby, these RNA molecules are noticed as potential regulators of neural functions in healthy circumstances, as well as neurological disorders. Moreover, circRNAs have also been reported to have potential diagnostic and therapeutic capacities in relation to AD, the most prevalent ND.

CONCLUSION

CircRNAs have been shown to act as sponges for miRNAs, thereby regulating the function of related miRNAs, including oxidative stress, reduction of neuroinflammation, and the formation and metabolism of Aβ, all of which developed in AD. CircRNAs have also been proposed as biomarkers that have potential diagnostic capacities in AD. Despite these characteristics, the use of circRNAs as therapeutic targets and promising diagnostic biomarkers will require further investigation and characterization of the function of these RNA molecules in AD.

Altered Transcriptional Regulation of Glycolysis in Circulating CD8+ T Cells of Rheumatoid Arthritis Patients

Peripheral T lymphocytes of rheumatoid arthritis (RA) patients show pathological changes in their metabolic pathways, especially glycolysis. These changes may drive the increased proliferation and tissue invasiveness of RA T cells. In order to study the transcriptional regulation underlying these alterations, we analysed publicly available RNA sequencing data from circulating T lymphocyte subsets of healthy individuals, untreated RA patients, and patients undergoing treatment for RA. Differential co-expression networks were created using sample-wise edge weights from an analysis called “linear interpolation to obtain network estimates for single sample” (lionessR), and annotated using the Gene Transcription Regulation Database (GTRD). Genes with high centrality scores were identified. CD8+ effector memory cells (Tem) and CD8+CD45RA+ effector memory cells (Temra) showed large changes in the transcriptional regulation of glycolysis in untreated RA. PFKFB3 and GAPDH were differentially regulated and had high centrality scores in CD8+ Tem cells. PFKFB3 downregulation may be due to HIF1A post transcriptional inhibition. Tocilizumab treatment partially reversed the RA-associated differential expression of several metabolic and regulatory genes. MYC was upregulated and had high centrality scores in RA CD8+ Temra cells; however, its glycolysis targets were unaltered. The upregulation of the PI3K-AKT and mTOR pathways may explain MYC upregulation.

Differential Expression of Mitosis and Cell Cycle Regulatory Genes during Recovery from an Acute Respiratory Virus Infection

Acute respiratory virus infections can have profound and long-term effects on lung function that persist even after the acute responses have fully resolved. In this study, we examined gene expression by RNA sequencing in the lung tissue of wild-type BALB/c mice that were recovering from a sublethal infection with the pneumonia virus of mice (PVM), a natural rodent pathogen of the same virus family and genus as the human respiratory syncytial virus. We compared these responses to gene expression in PVM-infected mice treated with Lactobacillus plantarum, an immunobiotic agent that limits inflammation and averts the negative clinical sequelae typically observed in response to acute infection with this pathogen. Our findings revealed prominent differential expression of inflammation-associated genes as well as numerous genes and gene families implicated in mitosis and cell-cycle regulation, including cyclins, cyclin-dependent kinases, cell division cycle genes, E2F transcription factors, kinesins, centromere proteins, and aurora kinases, among others. Of particular note was the differential expression of the cell division cycle gene Cdc20b, which was previously identified as critical for the ex vivo differentiation of multi-ciliated cells. Collectively, these findings provided us with substantial insight into post-viral repair processes and broadened our understanding of the mechanisms underlying Lactobacillus-mediated protection.

Efficient induction and sustenance of pluripotent stem cells from bovine somatic cells

Although derivation of naïve bovine embryonic stem cells is unachieved, the possibility for generation of bovine induced pluripotent stem cells (biPSCs) has been generally reported. However, attempts to sustain biPSCs by promoting self-renewal have not been successful. Methods established for maintaining murine and human induced pluripotent stem cells (iPSCs) do not support self-renewal of iPSCs for any bovid species. In this study, we examined methods to enhance complete reprogramming and concurrently investigated signaling relevant to pluripotency of the bovine blastocyst inner cell mass (ICM). First, we identified that forced expression of SV40 large T antigen together with the reprogramming genes (OCT4, SOX2, KLF4 and MYC) substantially enhanced the reprogramming efficacy of bovine fibroblasts to biPSCs. Second, we uncovered that TGFβ signaling is actively perturbed in the ICM. Inhibition of ALK4/5/7 to block TGFβ/activin/nodal signaling together with GSK3β and MEK1/2 supported robust in vitro self-renewal of naïve biPSCs with unvarying colony morphology, steady expansion, expected pluripotency gene expression and committed differentiation plasticity. Core similarities between biPSCs and stem cells of the 16-cell-stage bovine embryo indicated a stable ground state of pluripotency; this allowed us to reliably gain predictive understanding of signaling in bovine pluripotency using systems biology approaches. Beyond defining a high-fidelity platform for advancing biPSC-based biotechnologies that have not been previously practicable, these findings also represent a significant step towards understanding corollaries and divergent aspects of bovine pluripotency.

This article has an associated First Person interview with the joint first authors of the paper.

Summary: Pluripotency reprogramming by overcoming the stable epigenome of bovine cells, and uncovering precise early embryo self-renewal mechanisms enables sustenance and expansion of authentic induced pluripotent stem cells in vitro.

Identification of E2F transcription factor 7 as a novel potential biomarker for oral squamous cell carcinoma

BACKGROUND

As a tumor-accelerating transcriptional factor, E2F transcription factor 7 (E2F7) was up-regulated in many forms of cancers. Nevertheless, little has been reported about the impacts of E2F7 on oral squamous cell carcinoma (OSCC). Here, we aimed to probe whether E2F7 had influences on OSCC and its potential mechanism.

METHODS

The expression of E2F7 in OSCC tissues was analyzed using the data acquired from TCGA and ONCOMINE databases. E2F7 prognostic value in OSCC patients was analyzed utilizing TCGA database. The expression of E2F7 in OSCC cell lines was detected by qRT-PCR. Gain-and loss-function of E2F7 assays in TCA-83 and CAL27 cells were performed respectively to inquire the function of E2F7. Western blotting was applied to test the alternations of EMT-related markers.

RESULTS

In OSCC tissues, E2F7 was highly expressed. Besides, high expression of E2F7 predicted worse prognosis in OSCC patients. Moreover, E2F7 was over-expressed in TCA-83, HSC-4 and CAL27 (all OSCC cell lines) cells relative to that in HNOK (a normal cell line) cells. Gain-and loss-function assays displayed that deficiency of E2F7 suppresses CAL27 cell growth, migration, invasion and E2F7 high-expression resulted in inverse outcomes in TCA-83 cells. Finally, we found that silencing of E2F7 facilitated E-cadherin protein expression level and reduced N-cadherin, Vimentin and Snail protein levels in CAL27 cells, whilst E2F7 high-expression exhibited the opposite effects in TCA-83 cells.

CONCLUSIONS

These outcomes indicated that E2F7 performs a carcinogenic role in OSCC, which provides a theoretical basis for the therapeutic strategies of OSCC.

A FOXO-dependent replication checkpoint restricts proliferation of damaged cells

DNA replication is challenged by numerous exogenous and endogenous factors that can interfere with the progression of replication forks. Substantial accumulation of single-stranded DNA during DNA replication activates the DNA replication stress checkpoint response that slows progression from S/G2 to M phase to protect genomic integrity. Whether and how mild replication stress restricts proliferation remains controversial. Here, we identify a cell cycle exit mechanism that prevents S/G2 phase arrested cells from undergoing mitosis after exposure to mild replication stress through premature activation of the anaphase promoting complex/cyclosome (APC/C^CDH1). We find that replication stress causes a gradual decrease of the levels of the APC/C^CDH1 inhibitor EMI1/FBXO5 through Forkhead box O (FOXO)-mediated inhibition of its transcription factor E2F1. By doing so, FOXOs limit the time during which the replication stress checkpoint is reversible and thereby play an important role in maintaining genomic stability.

An integrated multi-omics approach identifies epigenetic alterations associated with Alzheimer's disease

Protein aggregation is the hallmark of neurodegeneration but the molecular mechanisms underlying late-onset Alzheimer’s disease (AD) remain unclear. Here we integrated transcriptomic, proteomic and epigenomic analyses of post-mortem human brains to identify molecular pathways involved in AD. RNA-seq analysis revealed upregulation of transcription- and chromatin-related genes, including the histone acetyltransferases for H3K27ac and H3K9ac. An unbiased proteomic screening singled out H3K27ac and H3K9ac as main enrichments specific to AD. In turn, epigenomic profiling revealed gains of H3K27ac and H3K9ac linked to transcription, chromatin, and disease pathways in AD. Increasing genome-wide H3K27ac and H3K9ac in a fly model of AD exacerbated amyloid-β42-driven neurodegeneration. Together, these findings suggest that AD involves a reconfiguration of the epigenome, where H3K27ac and H3K9ac impact disease pathways by dysregulating transcription- and chromatin-gene feedback loops. The identification of this process highlights potential epigenetic strategies for early-stage disease treatment.

Cell cycle oscillators underlying orderly proteolysis of E2F8

E2F8 is a transcriptional repressor that antagonizes E2F1 at the crossroads of the cell cycle, apoptosis, and cancer. Previously, we discovered that E2F8 is a direct target of the APC/C ubiquitin ligase. Nevertheless, it remains unknown how E2F8 is dynamically controlled throughout the entirety of the cell cycle. Here, using newly developed human cell-free systems that recapitulate distinct inter-mitotic and G1 phases and a continuous transition from prometaphase to G1, we reveal an interlocking dephosphorylation switch coordinating E2F8 degradation with mitotic exit and the activation of APC/C^Cdh1. Further, we uncover differential proteolysis rates for E2F8 at different points within G1 phase, accounting for its accumulation in late G1 while APC/C^Cdh1 is still active. Finally, we demonstrate that the F-box protein Cyclin F regulates E2F8 in G2-phase. Altogether, our data define E2F8 regulation throughout the cell cycle, illuminating an extensive coordination between phosphorylation, ubiquitination and transcription in mammalian cell cycle.

Circular RNA PRKCI promotes glioma cell progression by inhibiting microRNA-545

We here tested expression and potential functions of circular RNA PRKCI (circPRKCI) in human glioma. Our results show that circPRKCI is upregulated in human glioma tissues and glioma cells, correlating with downregulation of its potential target, microRNA-545 (miR-545). In A172 and primary human glioma cells, shRNA-mediated silencing of circPRKCI inhibited cancer cell growth, survival, proliferation, and migration. Conversely, ectopic circPRKCI overexpression promoted A172 cell progression. miR-545 is the primary target of circPRKCI in glioma cells. Forced overexpression of miR-545 mimicked circPRKCI shRNA-induced actions, inhibiting glioma cell survival and proliferation. In contrast, miR-545 inhibition, by a lentiviral antagomiR-545 construct, reversed circPRKCI shRNA-induced anti-A172 cell activity. Importantly, neither circPRKCI shRNA nor circPRKCI overexpression was effective in miR-545-knockout (Cas9 method) A172 cells. Importantly, the subcutaneous and orthotopic A172 xenograft growth was significantly inhibited by circPRKCI silencing. Collectively, circPRKCI promotes human glioma cell progression possibly by inhibiting miR-545. Targeting circPRKCI-miR-545 cascade could efficiently inhibit human glioma cells.

Zinc finger proteins orchestrate active gene silencing during embryonic stem cell differentiation

Transcription factors and chromatin remodeling proteins control the transcriptional variability for ESC lineage commitment. During ESC differentiation, chromatin modifiers are recruited to the regulatory regions by transcription factors, thereby activating the lineage-specific genes or silencing the transcription of active ESC genes. However, the underlying mechanisms that link transcription factors to exit from pluripotency are yet to be identified. In this study, we show that the Ctbp2-interacting zinc finger proteins, Zfp217 and Zfp516, function as linkers for the chromatin regulators during ESC differentiation. CRISPR-Cas9-mediated knock-outs of both Zfp217 and Zfp516 in ESCs prevent the exit from pluripotency. Both zinc finger proteins regulate the Ctbp2-mediated recruitment of the NuRD complex and polycomb repressive complex 2 (PRC2) to active ESC genes, subsequently switching the H3K27ac to H3K27me3 during ESC differentiation for active gene silencing. We therefore suggest that some zinc finger proteins orchestrate to control the concise epigenetic states on active ESC genes during differentiation, resulting in natural lineage commitment.

Association of Single-Nucleotide Polymorphism REX1 rs6815391, OCT4 rs13409 or rs3130932, and CTBP2 rs3740535 with Primary Lung Cancer Susceptibility: A Case-Control Study in a Chinese Population

The purpose of the current study is to explore the contribution of single-nucleotide polymorphisms (SNPs) of REX1 rs6815391, OCT4 rs13409 or rs3130932, and CTBP2 rs3740535 to the risk of lung cancer. A questionnaire survey was used to obtain basic information of the included subjects. A case control study was performed in 1121 patients and 1121 controls. All subjects were subjected to blood sampling for genomic DNA extraction and genotyping of the cancer stem cell-associated gene SNPs, including REX1 rs6815391, OCT4 rs13409 or rs3130932, and CTBP2 rs3740535 by real-time PCR. The association with the risk of primary lung cancer and interaction with environmental factors were assessed using unconditional logistic regression for the odds ratios and corresponding 95% confidence intervals. The genotype frequency distribution of OCT4 rs13409 loci was statistically significant, but there was no significant difference in the rest of the loci between lung cancer patients and healthy controls. The OCT4 gene was also related with lung cancer susceptibility in the genetic model after adjusting for lung cancer-related factors. Despite the presence of the dominant or recessive model, the four loci polymorphisms were associated with pollution near the place of residence, house type, worse ventilation situation, smoking, passive smoking, cooking oil fumes (COF), and family history of cancer, which increased the risk of lung cancer. Nonmarried status, 18.5≤BMI, COF, smoking, passive smoking, family history of cancer, and history of lung disease were independent risk factors of lung cancer susceptibility. Additionally, college degree or above, no pollution near the place of residence, protective genotype 1 or 2, and well ventilation can reduce the occurrence of lung cancer. There is an interaction between the four loci and environmental factors, and OCT4 rs13409 is a risk factor of primary lung cancer.

E2F1 and E2F7 differentially regulate KPNA2 to promote the development of gallbladder cancer

Karyopherin alpha 2 (KPNA2) is a nuclear import factor that is elevated in multiple cancers. However, its molecular regulation at the transcriptional levels is poorly understood. Here we found that KPNA2 was significantly upregulated in gallbladder cancer (GBC), and the increased levels were correlated with short survival of patients. Gene knocking down of KPNA2 inhibited tumor cell proliferation and migration in vitro as well as xenografted tumor development in vivo. A typical transcription factor E2F1 associated with its DNA-binding partner DP1 bond to the promoter region of KPNA2 and induced KPNA2 expression. In contrast, an atypical transcription factor E2F7 competed against DP1 and blocked E2F1-induced KPNA2 gene activation. Mutation of the dimerization residues of E2F7 or DNA-binding domain of E2F1 abolished the suppressive effects of E2F7 on KPNA2 gene expression. In addition, KPNA2 mediated nuclear localization of E2F1 and E2F7, where they in turn controlled KPNA2 expression. Taken together, our data provided mechanistic insights into divergently transcriptional regulation of KPNA2, thus pointing to KPNA2 as a potential target for cancer therapy.

Comprehensive molecular biomarker identification in breast cancer brain metastases

Background

Breast cancer brain metastases (BCBM) develop in about 20–30% of breast cancer (BC) patients. BCBM are associated with dismal prognosis not at least due to lack of valuable molecular therapeutic targets. The aim of the study was to identify new molecular biomarkers and targets in BCBM by using complementary state-of-the-art techniques.

Methods

We compared array expression profiles of three BCBM with 16 non-brain metastatic BC and 16 primary brain tumors (prBT) using a false discovery rate (FDR) p < 0.05 and fold change (FC) > 2. Biofunctional analysis was conducted on the differentially expressed probe sets. High-density arrays were employed to detect copy number variations (CNVs) and whole exome sequencing (WES) with paired-end reads of 150 bp was utilized to detect gene mutations in the three BCBM.

Results

The top 370 probe sets that were differentially expressed between BCBM and both BC and prBT were in the majority comparably overexpressed in BCBM and included, e.g. the coding genes BCL3, BNIP3, BNIP3P1, BRIP1, CASP14, CDC25A, DMBT1, IDH2, E2F1, MYCN, RAD51, RAD54L, and VDR. A number of small nucleolar RNAs (snoRNAs) were comparably overexpressed in BCBM and included SNORA1, SNORA2A, SNORA9, SNORA10, SNORA22, SNORA24, SNORA30, SNORA37, SNORA38, SNORA52, SNORA71A, SNORA71B, SNORA71C, SNORD13P2, SNORD15A, SNORD34, SNORD35A, SNORD41, SNORD53, and SCARNA22. The top canonical pathway was entitled, role of BRCA1 in DNA damage response. Network analysis revealed key nodes as Akt, ERK1/2, NFkB, and Ras in a predicted activation stage. Downregulated genes in a data set that was shared between BCBM and prBT comprised, e.g. BC cell line invasion markers JUN, MMP3, TFF1, and HAS2. Important cancer genes affected by CNVs included TP53, BRCA1, BRCA2, ERBB2, IDH1, and IDH2. WES detected numerous mutations, some of which affecting BC associated genes as CDH1, HEPACAM, and LOXHD1.

Conclusions

Using complementary molecular genetic techniques, this study identified shared and unshared molecular events in three highly aberrant BCBM emphasizing the challenge to detect new molecular biomarkers and targets with translational implications. Among new findings with the capacity to gain clinical relevance is the detection of overexpressed snoRNAs known to regulate some critical cellular functions as ribosome biogenesis.

Electronic supplementary material

The online version of this article (10.1186/s12967-017-1370-x) contains supplementary material, which is available to authorized users.

Pattern-based Search of Epigenomic Data Using GeNemo

Compared with the robust text-based search tools for genomic or RNA sequencing data, current methodologies for pattern-based searches of epigenomic and other functional genomic data are very limited. GeNemo is the first online search tool that accomplishes this goal. Users input their functional genomic data in the Browser Extensible Data (BED), Peaks, and bigWig formats, and may search for data in any of the three formats. Users may specify which types of datasets to search against, choosing from a variety of online datasets, with the Encyclopedia of DNA Elements (ENCODE) representing different epigenomic marks, transcriptional factor binding sites, and chromatin hypersensitivities or accessibilities in specific cell types, and developmental stages or species (mouse or human). GeNemo returns a list of genomic regions with matching patterns to the input data, which may be viewed in the browser as well as downloaded in the BED file format. The upgraded GeNemo has improved graphical display, has more robust interface, and is no longer prone to errors due to changes in the University of California, Santa Cruz (UCSC) genome browser. Troubleshooting steps for common problems are discussed. As the amount of functional genomic data is expanding exponentially, there is a critical need to develop and refine new bioinformatic tools such as GeNemo for data analyses and interpretation.

Identification of a Low-Frequency Missense Variant in E2F Transcription Factor 7 Associated with Colorectal Cancer Risk In A Chinese Population

Background: Transcription factors regulate gene expression and play important role in tumor genesis. Especially, the E2F transcription factor family controls the cell cycle and regulate many tumor suppressors. Missense variants in E2F family genes, which change the amino acid sequence, may alter the capacity for DNA binding or the protein structure, leading to a functional alteration. Material and Methods: We here searched for missense variants in E2F transcription family genes (E2F1~E2F8) and identified two (rs2075995 for E2F2 and rs3829295 for E2F7) with minor allele frequencies >0.01 in Chinese Han Beijing population from the 1000 genome project. We genotyped these two variants in 1,055 colorectal cancer (CRC) patients and 1,936 healthy controls using Taqman genotyping assays and assessed associations between SNPs and risk of CRC using logistic regression adjusted for gender and age. Results: We found rs3829295 at E2F7 to be significantly associated with risk of CRC. Compared with TT genotype carriers, CT and CT+CC genotype carriers had lower risks of CRC with ORs of 0.61 (95% CI: 0.44-0.85, P=0.003) and 0.61 (95% CI: 0.44-0.84, P=0.003), respectively. When stratified by gender and age, significant associations were observed in males (OR= 0.56, 95% CI: 0.38-0.83, P=0.004) for rs3829295, but not females (OR= 0.73, 95% CI: 0.43-1.22, P=0.232). Conclusion: Through a systematic assessment of variants in the E2F transcription factor family, we identified a lowfrequent missense variant in E2F7 significantly associated with CRC risk, indicating that E2F7 may play an important role in development of this tumor type.

Information-dependent enrichment analysis reveals time-dependent transcriptional regulation of the estrogen pathway of toxicity

The twenty-first century vision for toxicology involves a transition away from high-dose animal studies to in vitro and computational models (NRC in Toxicity testing in the 21st century: a vision and a strategy, The National Academies Press, Washington, DC, 2007). This transition requires mapping pathways of toxicity by understanding how in vitro systems respond to chemical perturbation. Uncovering transcription factors/signaling networks responsible for gene expression patterns is essential for defining pathways of toxicity, and ultimately, for determining the chemical modes of action through which a toxicant acts. Traditionally, transcription factor identification is achieved via chromatin immunoprecipitation studies and summarized by calculating which transcription factors are statistically associated with up- and downregulated genes. These lists are commonly determined via statistical or fold-change cutoffs, a procedure that is sensitive to statistical power and may not be as useful for determining transcription factor associations. To move away from an arbitrary statistical or fold-change-based cutoff, we developed, in the context of the Mapping the Human Toxome project, an enrichment paradigm called information-dependent enrichment analysis (IDEA) to guide identification of the transcription factor network. We used a test case of activation in MCF-7 cells by 17β estradiol (E2). Using this new approach, we established a time course for transcriptional and functional responses to E2. ERα and ERβ were associated with short-term transcriptional changes in response to E2. Sustained exposure led to recruitment of additional transcription factors and alteration of cell cycle machinery. TFAP2C and SOX2 were the transcription factors most highly correlated with dose. E2F7, E2F1, and Foxm1, which are involved in cell proliferation, were enriched only at 24 h. IDEA should be useful for identifying candidate pathways of toxicity. IDEA outperforms gene set enrichment analysis (GSEA) and provides similar results to weighted gene correlation network analysis, a platform that helps to identify genes not annotated to pathways.

E2F7 overexpression leads to tamoxifen resistance in breast cancer cells by competing with E2F1 at miR-15a/16 promoter

About 50-70% of breast cancers are estrogen receptor α (ERα) positive and most of them are sensitive to endocrine therapy including tamoxifen. However, one third of these patients will eventually develop resistance and relapse. We found that the expression of miR-15a and miR-16 were significantly decreased in tamoxifen resistant ER positive breast cancer cell lines. Exogenous expression of miR-15a/16 mimics re-sensitized resistant cells to tamoxifen by inhibiting Cyclin E1 and B cell lymphoma-2 (Bcl-2) to induce cell growth arrest and apoptosis respectively. Further, we identified that a repressive member of E2F family, E2F7, was responsible for the suppression of miR-15a/16 cluster by competing with E2F1 for E2F binding site at the promoter of their host gene DLEU2. Moreover, high expression of E2F7 is correlated with high risk of relapse and poor prognosis in breast cancer patients receiving tamoxifen treatment. Together, our results suggest that overexpression of E2F7 represses miR-15a/16 and then increases Cyclin E1 and Bcl-2 that result in tamoxifen resistance. E2F7 may be a valuable prognostic marker and a therapeutic target of tamoxifen resistance in breast cancer.

CDKN2A-independent role of BMI1 in promoting growth and survival of Ph+ acute lymphoblastic leukemia

BMI1 is a key component of the PRC1 complex (polycomb repressive complex-1) required for maintenance of normal and cancer stem cells. Its aberrant expression is detected in chronic myeloid leukemia and Ph+ acute lymphoblastic leukemia (ALL), but no data exist on BMI1 requirement in ALL cells.

We show here that BMI1 expression is important for proliferation and survival of Ph+ ALL cells and for leukemogenesis of Ph+ cells in vivo.

Levels of BIM, interferon-α (IFNα)-regulated genes, and E2F7 were upregulated in BMI1-silenced cells, suggesting that repressing their expression is important for BMI1 biological effects. Consistent with this hypothesis, we found that: i) downregulation of BIM or E2F7 abrogated apoptosis or rescued, in part, the reduced proliferation and colony formation of BMI1 silenced BV173 cells; ii) BIM/E2F7-double silencing further enhanced colony formation and in vivo leukemogenesis of BMI1-silenced cells; iii) overexpression of BIM and E2F7 mimicked the effect of BMI1 silencing in BV173 and SUP-B15 cells and iv) treatment with IFNα suppressed proliferation and colony formation of Ph+ ALL cells.

These studies indicate that the growth-promoting effects of BMI1 in Ph+ ALL cells depend on suppression of multiple pathways and support the use of IFNα in the therapy of Ph+ ALL.

E2F Transcription Factors Control the Roller Coaster Ride of Cell Cycle Gene Expression

Initially, the E2F transcription factor was discovered as a factor able to bind the adenovirus E2 promoter and activate viral genes. Afterwards it was shown that E2F also binds to promoters of nonviral genes such as C-MYC and DHFR, which were already known at that time to be important for cell growth and DNA metabolism, respectively. These findings provided the first clues that the E2F transcription factor might be an important regulator of the cell cycle. Since this initial discovery in 1987, several additional E2F family members have been identified, and more than 100 targets genes have been shown to be directly regulated by E2Fs, the majority of these are important for controlling the cell cycle. The progression of a cell through the cell cycle is accompanied with the increased expression of a specific set of genes during one phase of the cell cycle and the decrease of the same set of genes during a later phase of the cell cycle. This roller coaster ride, or oscillation, of gene expression is essential for the proper progression through the cell cycle to allow accurate DNA replication and cell division. The E2F transcription factors have been shown to be critical for the temporal expression of the oscillating cell cycle genes. This review will focus on how the oscillation of E2Fs and their targets is regulated by transcriptional, post-transcriptional and post-translational mechanism in mammals, yeast, flies, and worms. Furthermore, we will discuss the functional impact of E2Fs on the cell cycle progression and outline the consequences when E2F expression is disturbed.

Genome-wide analysis reveals NRP1 as a direct HIF1α-E2F7 target in the regulation of motorneuron guidance in vivo

In this study, we explored the existence of a transcriptional network co-regulated by E2F7 and HIF1α, as we show that expression of E2F7, like HIF1α, is induced in hypoxia, and because of the previously reported ability of E2F7 to interact with HIF1α. Our genome-wide analysis uncovers a transcriptional network that is directly controlled by HIF1α and E2F7, and demonstrates both stimulatory and repressive functions of the HIF1α -E2F7 complex. Among this network we reveal Neuropilin 1 (NRP1) as a HIF1α-E2F7 repressed gene. By performing in vitro and in vivo reporter assays we demonstrate that the HIF1α-E2F7 mediated NRP1 repression depends on a 41 base pairs ‘E2F-binding site hub’, providing a molecular mechanism for a previously unanticipated role for HIF1α in transcriptional repression. To explore the biological significance of this regulation we performed in situ hybridizations and observed enhanced nrp1a expression in spinal motorneurons (MN) of zebrafish embryos, upon morpholino-inhibition of e2f7/8 or hif1α. Consistent with the chemo-repellent role of nrp1a, morpholino-inhibition of e2f7/8 or hif1α caused MN truncations, which was rescued in TALEN-induced nrp1a^hu10012 mutants, and phenocopied in e2f7/8 mutant zebrafish. Therefore, we conclude that repression of NRP1 by the HIF1α-E2F7 complex regulates MN axon guidance in vivo.

Structural insights into the DNA-binding specificity of E2F family transcription factors

The mammalian cell cycle is controlled by the E2F family of transcription factors. Typical E2Fs bind to DNA as heterodimers with the related dimerization partner (DP) proteins, whereas the atypical E2Fs, E2F7 and E2F8 contain two DNA-binding domains (DBDs) and act as repressors. To understand the mechanism of repression, we have resolved the structure of E2F8 in complex with DNA at atomic resolution. We find that the first and second DBDs of E2F8 resemble the DBDs of typical E2F and DP proteins, respectively. Using molecular dynamics simulations, biochemical affinity measurements and chromatin immunoprecipitation, we further show that both atypical and typical E2Fs bind to similar DNA sequences in vitro and in vivo. Our results represent the first crystal structure of an E2F protein with two DBDs, and reveal the mechanism by which atypical E2Fs can repress canonical E2F target genes and exert their negative influence on cell cycle progression.

Unbiased proteomic profiling strategy for discovery of bacterial effector proteins reveals that Salmonella protein PheA is a host cell cycle regulator

Salmonella utilizes a type III secretion system to inject bacterial effector proteins into the host cell cytosol. Once in the cytosol, these effectors hijack various biochemical pathways to regulate virulence. Despite the importance of effector proteins, especially for understanding host-pathogen interactions, a potentially large number of effectors are yet to be identified. Here, we demonstrate that unbiased chemical proteomic profiling using off-the-shelf fluorescent probes leads to the discovery of a host cell cycle regulator encoded in the Salmonella genome. Our profiling combined with bioinformatic analysis implicates 29 Salmonella as potential effectors. We follow up on the top candidate, chorismate mutase-P/prehenate dehydratase, PheA, and present evidence that PheA is an effector that mimics E2F7 transcription factor of the host cell and promotes G1/S cell cycle arrest. This validates our strategy and opens opportunities for effector identification in the future.

Downregulation of E2F1 during ER stress is required to induce apoptosis

The endoplasmic reticulum (ER) has recently emerged as an alternative target to induce cell death in tumours, because prolonged ER stress results in the induction of apoptosis even in chemoresistant transformed cells. Here, we show that the DNA-damage-responsive pro-apoptotic factor E2F1 is unexpectedly downregulated during the ER stress-mediated apoptotic programme. E2F1 decline is a late event during the ER response and is mediated by the two unfolded protein response (UPR) sensors ATF6 and IRE1 (also known as ERN1). Whereas ATF6 directly interacts with the E2F1 promoter, IRE1 requires the involvement of the known E2F1 modulator E2F7, through the activation of its main target Xbp-1. Importantly, inhibition of the E2F1 decrease prevents ER-stress-induced apoptosis, whereas E2F1 knockdown efficiently sensitises cells to ER stress-dependent apoptosis, leading to the upregulation of two main factors in the UPR pro-apoptotic execution phase, Puma and Noxa (also known as BBC3 and PMAIP1, respectively). Our results point to a novel key role of E2F1 in the cell survival/death decision under ER stress, and unveil E2F1 inactivation as a valuable novel potential therapeutic strategy to increase the response of tumour cells to ER stress-based anticancer treatments.

Transcriptome-wide signatures of tumor stage in kidney renal clear cell carcinoma: connecting copy number variation, methylation and transcription factor activity

BACKGROUND

Comparative analysis of expression profiles between early and late stage cancers can help to understand cancer progression and metastasis mechanisms and to predict the clinical aggressiveness of cancer. The observed stage-dependent expression changes can be explained by genetic and epigenetic alterations as well as transcription dysregulation. Unlike genetic and epigenetic alterations, however, activity changes of transcription factors, generally occurring at the post-transcriptional or post-translational level, are hard to detect and quantify.

METHODS

Here we developed a statistical framework to infer the activity changes of transcription factors by simultaneously taking into account the contributions of genetic and epigenetic alterations to mRNA expression variations.

RESULTS

Applied to kidney renal clear cell carcinoma (KIRC), the model underscored the role of methylation as a significant contributor to stage-dependent expression alterations and identified key transcription factors as potential drivers of cancer progression.

CONCLUSIONS

Integrating copy number, methylation, and transcription factor activity signatures to explain stage-dependent expression alterations presented a precise and comprehensive view on the underlying mechanisms during KIRC progression.

Influence of GSH synthesis inhibition on temporal distribution of NAD+/NADH during vascular endothelial cells proliferation

Pathological conditions states such as stroke, diabetes mellitus, hypertension, dyslipidemia are associated with increased levels of free radicals that alter normal function of the vascular endothelium and perturb vascular homeostasis. The redox couples reduced glutathione (GSH)/oxidized glutathione (GSSG), NADH/NAD+, and NADPH/NADP+ play major functions in the intracellular redox balance. Any decrease in tissue or systemic GSH levels under the aforementioned pathologies would enhance oxidative damage to the vascular endothelium. Beside their role as coenzyme that participate in cellular metabolism, pyridine nucleotides serve also as substrate for enzymes involved in DNA repair and longevity. There is scant data on NAD+/NADH kinetics and distribution during human cells proliferation. Here, we determined the influence of cellular GSH status on the early dynamics of nuclear-to-cytosol (N-to-C) NAD+ and nuclear NADH kinetics (6 h interval) over 72 h of endothelial cell proliferation. The IHEC cell line was used as a surrogate for human brain micro vascular endothelial cells. Inhibition of GSH synthesis by buthionine sulfoximine (BSO) and sustained low cellular GSH significantly increased nuclear NADH levels (p<0.01), which correlated with lower nuclear GSH and prolonged cell cycle S-phase. When BSO was removed the pattern of nuclear NAD+ resembled that of control group, but nuclear NADH concentrations remained elevated, as in GSH deficient cells (p<0.01). The coincidence of high nuclear NADH and lower nuclear NAD+ with S-phase prolongation are suggestive of CtBP and NAD+-dependent DNA repair enzyme activation under conditions of decreased cellular GSH. These results provide important insights into GSH control of vascular endothelial growth and restitution, key processes in the restoration of the endothelium adjacent to the post-injury lesion site.

CtBP2 proteome: Role of CtBP in E2F7-mediated repression and cell proliferation

C-terminal binding protein (CtBP) family transcriptional corepressors include CtBP1 and CtBP2. While CtBP1 and CtBP2 share significant amino acid sequence homology, CtBP2 possesses a unique N-terminal domain that is modified by acetylation and contributes to exclusive nuclear localization. Although CtBP1 and CtBP2 are functionally redundant for certain activities during vertebrate development, they also perform unique functions. Previous studies have identified several CtBP1-interacting proteins that included other transcriptional corepressors, DNA-binding repressors and histone modifying enzymatic components such as the histone deacetylases and the histone demethylase LSD-1. Here, we carried out an unbiased proteomic analysis of CtBP2-associated proteins and discovered the association of several components of the CtBP1 proteome as well as novel interactions. The CtBP2 proteome contained components of the NuRD complex and the E2F family member E2F7. E2F7 interacted with the hydrophobic cleft region of CtBP1 and CtBP2 through a prototypical CtBP binding motif, PIDLS. E2F7 repressed E2F1 transcription, inhibited cell proliferation in a CtBP-dependent fashion. Our study identified CtBP as a corepressor of E2F7 and as a regulator of DNA damage response.