Wide-ranging functions of E2F4 in transcriptional activation and repression revealed by genome-wide analysis

DeepCG: A cell graph model for predicting prognosis in lung adenocarcinoma

Lung cancer is the first leading cause of cancer-related death in the United States, with lung adenocarcinoma as the major subtype accounting for 40% of all cases. To improve patient survival, image-based prognostic models were developed due to the ready availability of pathological images at diagnosis. However, the application of these models is hampered by two main challenges: the lack of publicly available image datasets with high-quality survival information and the poor interpretability of conventional convolutional neural network models. Here, we integrated matched transcriptomic and H&E staining data from TCGA (The Cancer Genome Atlas) to develop an image-based prognostic model, termed Deep-learning based Cell Graph (DeepCG) model. Instead of survival data, we used a gene signature to predict patient prognostic risks, which was then used as labels for training DeepCG. Importantly, by employing graph structures to capture cell patterns, DeepCG can provide cell-level interpretation, which was more biologically relevant than previous region-level insights. We validated the prognostic values of DeepCG in independent datasets and demonstrated its ability to identify prognostically informative cells in images.

The pRb/RBL2-E2F1/4-GCN5 axis regulates cancer stem cell formation and G0 phase entry/exit by paracrine mechanisms

The lethality, chemoresistance and metastatic characteristics of cancers are associated with phenotypically plastic cancer stem cells (CSCs). How the non-cell autonomous signalling pathways and cell-autonomous transcriptional machinery orchestrate the stem cell-like characteristics of CSCs is still poorly understood. Here we use a quantitative proteomic approach for identifying secreted proteins of CSCs in pancreatic cancer. We uncover that the cell-autonomous E2F1/4-pRb/RBL2 axis balances non-cell-autonomous signalling in healthy ductal cells but becomes deregulated upon KRAS mutation. E2F1 and E2F4 induce whereas pRb/RBL2 reduce WNT ligand expression (e.g. WNT7A, WNT7B, WNT10A, WNT4) thereby regulating self-renewal, chemoresistance and invasiveness of CSCs in both PDAC and breast cancer, and fibroblast proliferation. Screening for epigenetic enzymes identifies GCN5 as a regulator of CSCs that deposits H3K9ac onto WNT promoters and enhancers. Collectively, paracrine signalling pathways are controlled by the E2F-GCN5-RB axis in diverse cancers and this could be a therapeutic target for eliminating CSCs.

Analysis of E2F1 single-nucleotide polymorphisms reveals deleterious non-synonymous substitutions that disrupt E2F1-RB protein interaction in cancer

Cancer is a medical condition that is caused by the abnormal growth and division of cells, leading to the formation of tumors. The E2F1 and RB pathways are critical in regulating cell cycle, and their dysregulation can contribute to the development of cancer. In this study, we analyzed experimentally reported SNPs in E2F1 and assessed their effects on the binding affinity with RB. Out of 46, nine mutations were predicted as deleterious, and further analysis revealed four highly destabilizing mutations (L206W, R232C, I254T, A267T) that significantly altered the protein structure. Molecular docking of wild-type and mutant E2F1 with RB revealed a docking score of -242 kcal/mol for wild-type, while the mutant complexes had scores ranging from -217 to -220 kcal/mol. Molecular simulation analysis revealed variations in the dynamics features of both mutant and wild-type complexes due to the acquired mutations. Furthermore, the total binding free energy for the wild-type E2F1-RB complex was -64.89 kcal/mol, while those of the L206W, R232C, I254T, and A267T E2F1-RB mutants were -45.90 kcal/mol, -53.52 kcal/mol, -55.67 kcal/mol, and -61.22 kcal/mol, respectively. Our study is the first to extensively analyze E2F1 gene mutations and identifies candidate mutations for further validation and potential targeting for cancer therapeutics.

Dynamic and distinct histone modifications facilitate human trophoblast lineage differentiation

The placenta serves as an essential organ for fetal growth throughout pregnancy. Histone modification is a crucial regulatory mechanism involved in numerous biological processes and development. Nevertheless, there remains a significant gap in our understanding regarding the epigenetic regulations that influence trophoblast lineage differentiation, a fundamental aspect of placental development. Here, through comprehensive mapping of H3K4me3, H3K27me3, H3K9me3, and H3K27ac loci during the differentiation of trophoblast stem cells (TSCs) into syncytiotrophoblasts (STs) and extravillous trophoblasts (EVTs), we reveal dynamic reconfiguration in H3K4me3 and H3K27ac patterns that establish an epigenetic landscape conducive to proper trophoblast lineage differentiation. We observe that broad H3K4me3 domains are associated with trophoblast lineage-specific gene expression. Unlike embryonic stem cells, TSCs lack robust bivalent domains. Notably, the repression of ST- and EVT-active genes in TSCs is primarily attributed to the weak H3K4me3 signal rather than bivalent domains. We also unveil the inactivation of TSC enhancers precedes the activation of ST enhancers during ST formation. Our results provide a comprehensive global map of diverse histone modifications, elucidating the dynamic histone modifications during trophoblast lineage differentiation.

HPV38 impairs UV-induced transcriptional activation of the IL-18 pro-inflammatory cytokine

IMPORTANCE

Here, we demonstrate that the direct binding of p53 on the IL-18 promoter region regulates its gene expression. However, the presence of E6 and E7 from human papillomavirus type 38 impairs this mechanism via a new inhibitory complex formed by DNA methyltransferase 1 (DNMT1)/PKR/ΔNp73α, which binds to the region formerly occupied by p53 in primary keratinocytes.

E2F8 exerts cancer-promoting effects by transcriptionally activating RRM2 and E2F8 knockdown synergizes with WEE1 inhibition in suppressing lung adenocarcinoma

Ribonucleotide reductase (RR) is a rate-limiting enzyme that facilitates DNA replication and repair by reducing nucleotide diphosphates (NDPs) to deoxyribonucleotide diphosphates (dNDPs) and is thereby crucial for cell proliferation and cancer development. The E2F family of transcription factors includes key regulators of gene expression involved in cell cycle control. In this study, E2F8 expression was significantly increased in most cancer tissues of lung adenocarcinoma (LUAD) patients and was correlated with the expression of RRM2 through database and clinical samples analysis. The protein expression of E2F8 and RRM2 were positively correlated with tumor-node-metastasis (TNM) pathological stage, and high expression of E2F8 and RRM2 predicted a low 5-year overall survival rate in LUAD patients. Overexpression and knockdown experiments showed that E2F8 was essential for LUAD cell proliferation, DNA synthesis, and cell cycle progression, which were RRM2-dependent. Reporter gene, ChIP-qPCR, and DNA pulldown-Western blot assays indicated that E2F8 activated the transcription of the RRM2 gene by directly binding with the RRM2 promoter in LUAD cells. Previous studies indicated that inhibition of WEE1 kinase can suppress the phosphorylation of CDK1/2 and promote the degradation of RRM2. We further showed here that the combination of E2F8 knockdown with MK-1775, an inhibitor of WEE1 being evaluated in clinical trials, synergistically suppressed proliferation and promoted apoptosis of LUAD cells in vitro and in vivo. Thus, this study reveals a novel role of E2F8 as a proto-oncogenic transcription activator by activating RRM2 expression in LUAD, and targeting both the transcription and degradation mechanisms of RRM2 could produce a synergistic inhibitory effect for LUAD treatment in addition to conventional inhibition of RR enzyme activity.

MicroRNAs as the critical regulators of tumor angiogenesis in liver cancer

Liver cancer is one of the most common malignancies in human digestive system. Despite the recent therapeutic methods, there is a high rate of mortality among liver cancer patients. Late diagnosis in the advanced tumor stages can be one of the main reasons for the poor prognosis in these patients. Therefore, investigating the molecular mechanisms of liver cancer can be helpful for the early stage tumor detection and treatment. Vascular expansion in liver tumors can be one of the important reasons for poor prognosis and aggressiveness. Therefore, anti-angiogenic drugs are widely used in liver cancer patients. MicroRNAs (miRNAs) have key roles in the regulation of angiogenesis in liver tumors. Due to the high stability of miRNAs in body fluids, these factors are widely used as the non-invasive diagnostic and prognostic markers in cancer patients. Regarding, the importance of angiogenesis during liver tumor growth and invasion, in the present review, we discussed the role of miRNAs in regulation of angiogenesis in these tumors. It has been reported that miRNAs mainly exert an anti-angiogenic function by regulation of tumor microenvironment, transcription factors, and signaling pathways in liver tumors. This review can be an effective step to suggest the miRNAs for the non-invasive early detection of malignant and invasive liver tumors.

Interfering with Dusp2 alleviates high glucose-induced vascular endothelial cell dysfunction by promoting p38 MAPK pathway activation

BACKGROUND

Hyperglycemia-induced vascular endothelial cell dysfunction is a major factor contributing to diabetic lower extremity ischemia. We intend to investigate the role of Dusp2 in hyperglycemia-induced vascular endothelial cell dysfunction and related mechanisms.

METHODS

The human umbilical vein endothelial cells (HUVECs) were treated with high glucose (HG) as the cell model. Streptozotocin injection was performed to induce diabetes and femoral artery ligation was to induce hind limb ischemia in mice. The levels of Dusp2, p-p38 MAPK, E2F4, and p38 MAPK were evaluated by Western blot or quantitative real-time PCR. The laser Doppler perfusion imaging was conducted to measure blood flow recovery. The cell counting kit-8, transwell, and tube formation assay were performed to evaluate cell proliferation, migration, and angiogenesis, respectively. CD31 immunohistochemical staining was carried out to detect the capillary density of gastrocnemius. The dual-luciferase reporter gene assay and Chromatin immunoprecipitation assay were executed to explore the interaction between E2F4 and Dusp2.

RESULTS

Dusp2 was highly expressed in HG-induced HUVECs and diabetic lower extremity ischemia model mice. Interference with Dusp2 promoted cell proliferation, migration, and angiogenesis, as well as alleviated mouse diabetic hindlimb ischemia. Dusp2 knockdown up-regulated p-p38 MAPK levels. We verified the binding between E2F4 and Dusp2. Overexpressing E2F4 suppressed Dusp2 levels and promoted cell proliferation, migration, and angiogenesis, co-overexpression of Dusp2 reversed the results.

CONCLUSIONS

Overexpressing E2F4 promotes endothelial cell proliferation, migration, and angiogenesis by inhibiting Dusp2 expression and activating p38 MAPK to alleviate vascular endothelial cell dysfunction under HG stimulation.

The E2F4/p130 Repressor Complex Cooperates with Oncogenic ΔNp73α To Inhibit Gene Expression in Human Papillomavirus 38 E6/E7-Transformed Keratinocytes and in Cancer Cells

Tumor suppressor p53 and its related proteins, p63 and p73, can be synthesized as multiple isoforms lacking part of the N- or C-terminal regions. Specifically, high expression of the ΔNp73α isoform is notoriously associated with various human malignancies characterized by poor prognosis. This isoform is also accumulated by oncogenic viruses, such as Epstein-Barr virus (EBV), as well as genus beta human papillomaviruses (HPV) that appear to be involved in carcinogenesis. To gain additional insight into ΔNp73α mechanisms, we have performed proteomics analyses using human keratinocytes transformed by the E6 and E7 proteins of the beta-HPV type 38 virus as an experimental model (38HK). We find that ΔNp73α associates with the E2F4/p130 repressor complex through a direct interaction with E2F4. This interaction is favored by the N-terminal truncation of p73 characteristic of ΔNp73 isoforms. Moreover, it is independent of the C-terminal splicing status, suggesting that it could represent a general feature of ΔNp73 isoforms (α, β, γ, δ, ε, ζ, θ, η, and η1). We show that the ΔNp73α-E2F4/p130 complex inhibits the expression of specific genes, including genes encoding for negative regulators of proliferation, both in 38HK and in HPV-negative cancer-derived cell lines. Such genes are not inhibited by E2F4/p130 in primary keratinocytes lacking ΔNp73α, indicating that the interaction with ΔNp73α rewires the E2F4 transcriptional program. In conclusion, we have identified and characterized a novel transcriptional regulatory complex with potential implications in oncogenesis. IMPORTANCE The TP53 gene is mutated in about 50% of human cancers. In contrast, the TP63 and TP73 genes are rarely mutated but rather expressed as ΔNp63 and ΔNp73 isoforms in a wide range of malignancies, where they act as p53 antagonists. Accumulation of ΔNp63 and ΔNp73, which is associated with chemoresistance, can result from infection by oncogenic viruses such as EBV or HPV. Our study focuses on the highly carcinogenic ΔNp73α isoform and uses a viral model of cellular transformation. We unveil a physical interaction between ΔNp73α and the E2F4/p130 complex involved in cell cycle control, which rewires the E2F4/p130 transcriptional program. Our work shows that ΔNp73 isoforms can establish interactions with proteins that do not bind to the TAp73α tumor suppressor. This situation is analogous to the gain-of-function interactions of p53 mutants supporting cellular proliferation.

E2F2 is upregulated by the ERK pathway and regulates decidualization via MCM4

Cell cycle modulation is an important event during decidualization. E2F2 is a transcription regulator that plays a vital role in cell cycle regulation. However, the biological role of E2F2 in decidualization has not yet been identified. In this study, estrogen (E2) and progestin (P4)-induced in vitro and in vivo decidualization models were applied. Our data showed that the expression levels of E2F2 and its downstream target MCM4 were downregulated in the uterus tissues of E2P4-treated mice compared with control mice. In hESCs, exposure to E2P4 resulted in a significant decrease in E2F2 and MCM4 expression. E2P4 treatment reduced hESC proliferation and ectopic expression of E2F2 or MCM4 elevated the viability of E2P4-treated hESCs. In addition, ectopic expression of E2F2 or MCM4 restored the expression of G1 phase-associated proteins. The ERK pathway was inactivated in E2P4-treated hESCs. Treatment with ERK agonist Ro 67-7476 restored the expression of E2F2, MCM4, and G1 phase-associated proteins that were inhibited by E2P4. Moreover, Ro 67-7476 retracted the levels of IGFBP1 and PRL that were induced by E2P4. Collectively, our results indicate that E2F2 is regulated by ERK signaling and contributes to decidualization via regulation of MCM4. Therefore, E2F2/MCM4 cascade may serve as promising targets for alleviating decidualization dysfunction.

GEMC1 and MCIDAS interactions with SWI/SNF complexes regulate the multiciliated cell-specific transcriptional program

Multiciliated cells (MCCs) project dozens to hundreds of motile cilia from their apical surface to promote the movement of fluids or gametes in the mammalian brain, airway or reproductive organs. Differentiation of MCCs requires the sequential action of the Geminin family transcriptional activators, GEMC1 and MCIDAS, that both interact with E2F4/5-DP1. How these factors activate transcription and the extent to which they play redundant functions remains poorly understood. Here, we demonstrate that the transcriptional targets and proximal proteomes of GEMC1 and MCIDAS are highly similar. However, we identified distinct interactions with SWI/SNF subcomplexes; GEMC1 interacts primarily with the ARID1A containing BAF complex while MCIDAS interacts primarily with BRD9 containing ncBAF complexes. Treatment with a BRD9 inhibitor impaired MCIDAS-mediated activation of several target genes and compromised the MCC differentiation program in multiple cell based models. Our data suggest that the differential engagement of distinct SWI/SNF subcomplexes by GEMC1 and MCIDAS is required for MCC-specific transcriptional regulation and mediated by their distinct C-terminal domains.

E2F4DN Transgenic Mice: A Tool for the Evaluation of E2F4 as a Therapeutic Target in Neuropathology and Brain Aging

E2F4 was initially described as a transcription factor with a key function in the regulation of cell quiescence. Nevertheless, a number of recent studies have established that E2F4 can also play a relevant role in cell and tissue homeostasis, as well as tissue regeneration. For these non-canonical functions, E2F4 can also act in the cytoplasm, where it is able to interact with many homeostatic and synaptic regulators. Since E2F4 is expressed in the nervous system, it may fulfill a crucial role in brain function and homeostasis, being a promising multifactorial target for neurodegenerative diseases and brain aging. The regulation of E2F4 is complex, as it can be chemically modified through acetylation, from which we present evidence in the brain, as well as methylation, and phosphorylation. The phosphorylation of E2F4 within a conserved threonine motif induces cell cycle re-entry in neurons, while a dominant negative form of E2F4 (E2F4DN), in which the conserved threonines have been substituted by alanines, has been shown to act as a multifactorial therapeutic agent for Alzheimer’s disease (AD). We generated transgenic mice neuronally expressing E2F4DN. We have recently shown using this mouse strain that expression of E2F4DN in 5xFAD mice, a known murine model of AD, improved cognitive function, reduced neuronal tetraploidization, and induced a transcriptional program consistent with modulation of amyloid-β (Aβ) peptide proteostasis and brain homeostasis recovery. 5xFAD/E2F4DN mice also showed reduced microgliosis and astrogliosis in both the cerebral cortex and hippocampus at 3-6 months of age. Here, we analyzed the immune response in 1 year-old 5xFAD/E2F4DN mice, concluding that reduced microgliosis and astrogliosis is maintained at this late stage. In addition, the expression of E2F4DN also reduced age-associated microgliosis in wild-type mice, thus stressing its role as a brain homeostatic agent. We conclude that E2F4DN transgenic mice represent a promising tool for the evaluation of E2F4 as a therapeutic target in neuropathology and brain aging.

Genetic Associations and Differential mRNA Expression Levels of Host Genes Suggest a Viral Trigger for Endemic Pemphigus Foliaceus

The long search for the environmental trigger of the endemic pemphigus foliaceus (EPF, fogo selvagem) has not yet resulted in any tangible findings. Here, we searched for genetic associations and the differential expression of host genes involved in early viral infections and innate antiviral defense. Genetic variants could alter the structure, expression sites, or levels of the gene products, impacting their functions. By analyzing 3063 variants of 166 candidate genes in 227 EPF patients and 194 controls, we found 12 variants within 11 genes associated with differential susceptibility (p < 0.005) to EPF. The products of genes TRIM5, TPCN2, EIF4E, EIF4E3, NUP37, NUP50, NUP88, TPR, USP15, IRF8, and JAK1 are involved in different mechanisms of viral control, for example, the regulation of viral entry into the host cell or recognition of viral nucleic acids and proteins. Only two of nine variants were also associated in an independent German cohort of sporadic PF (75 patients, 150 controls), aligning with our hypothesis that antiviral host genes play a major role in EPF due to a specific virus−human interaction in the endemic region. Moreover, CCL5, P4HB, and APOBEC3G mRNA levels were increased (p < 0.001) in CD4+ T lymphocytes of EPF patients. Because there is limited or no evidence that these genes are involved in autoimmunity, their crucial role in antiviral responses and the associations that we observed support the hypothesis of a viral trigger for EPF, presumably a still unnoticed flavivirus. This work opens new frontiers in searching for the trigger of EPF, with the potential to advance translational research that aims for disease prevention and treatment.

KLF15 cistromes reveal a hepatocyte pathway governing plasma corticosteroid transport and systemic inflammation

Circulating corticosteroids orchestrate stress adaptation, including inhibition of inflammation. While pathways governing corticosteroid biosynthesis and intracellular signaling are well understood, less is known about mechanisms controlling plasma corticosteroid transport. Here, we show that hepatocyte KLF15 (Kruppel-like factor 15) controls plasma corticosteroid transport and inflammatory responses through direct transcriptional activation of Serpina6, which encodes corticosteroid-binding globulin (CBG). Klf15-deficient mice have profoundly low CBG, reduced plasma corticosteroid binding capacity, and heightened mortality during inflammatory stress. These defects are completely rescued by reconstituting CBG, supporting that KLF15 works primarily through CBG to control plasma corticosterone homeostasis. To understand transcriptional mechanisms, we generated the first KLF15 cistromes using newly engineered Klf153xFLAG mice. Unexpectedly, liver KLF15 is predominantly promoter enriched, including Serpina6, where it binds a palindromic GC-rich motif, opens chromatin, and transactivates genes with minimal associated direct gene repression. Overall, we provide critical mechanistic insight into KLF15 function and identify a hepatocyte-intrinsic transcriptional module that potently regulates systemic corticosteroid transport and inflammation.

A Mutant Variant of E2F4 Triggers Multifactorial Therapeutic Effects in 5xFAD Mice

Alzheimer’s disease (AD) has a complex etiology, which requires a multifactorial approach for an efficient treatment. We have focused on E2 factor 4 (E2F4), a transcription factor that regulates cell quiescence and tissue homeostasis, controls gene networks affected in AD, and is upregulated in the brains of Alzheimer’s patients and of APP^swe/PS1^dE9 and 5xFAD transgenic mice. E2F4 contains an evolutionarily conserved Thr-motif that, when phosphorylated, modulates its activity, thus constituting a potential target for intervention. In this study, we generated a knock-in mouse strain with neuronal expression of a mouse E2F4 variant lacking this Thr-motif (E2F4DN), which was mated with 5xFAD mice. Here, we show that neuronal expression of E2F4DN in 5xFAD mice potentiates a transcriptional program consistent with the attenuation of the immune response and brain homeostasis. This correlates with reduced microgliosis and astrogliosis, modulation of amyloid-β peptide proteostasis, and blocking of neuronal tetraploidization. Moreover, E2F4DN prevents cognitive impairment and body weight loss, a known somatic alteration associated with AD. We also show that our finding is significant for AD, since E2F4 is expressed in cortical neurons from Alzheimer patients in association with Thr-specific phosphorylation, as evidenced by an anti-E2F4/anti-phosphoThr proximity ligation assay. We propose E2F4DN-based gene therapy as a promising multifactorial approach against AD.

ZMIZ2 promotes the development of triple-receptor negative breast cancer

Background

Triple-receptor negative breast cancer (TNBC) is an aggressive breast tumor subtype that generally has a poor prognosis. This study aimed to investigate the role and regulatory mechanisms of Zinc finger MIZ-type containing 2 (ZMIZ2) in relation to TNBC.

Methods

Based on data from The Cancer Genome Atlas (TCGA), the expression of ZMIZ2 in different subtypes and its correlation with androgen receptor (AR) were analyzed, and a regulatory mechanism network was constructed. The expression and prognostic value of ZMIZ2 in clinical TNBC tissue samples were also investigated. Furthermore, in vitro studies were conducted to investigate the effects of ZMIZ2 knockdown on the malignant behaviors of TNBC cells and target gene expression.

Results

Based on TCGA data, ZMIZ2 was found to be significantly upregulated in TNBC tissues and its expression was negatively correlated with AR expression. Key relationships, such as the ZMIZ2-CCL5, ZMIZ2/AR-MCM3, ZMIZ2/AR-E2F4, and the ZMIZ2/AR-DHX38 were identified, which were enriched in NOD-like receptor signaling pathway/toll-like receptor signaling pathway, DNA replication, cell cycle, and spliceosome, respectively. Moreover, ZMIZ2 was upregulated in clinical breast cancer tissues and its high expression was correlated with the poor prognosis of TNBC patients. Furthermore, ZMIZ2 expression was increased in breast cancer cells, and a knockdown of ZMIZ2 inhibited MDA-MB-231 cell proliferation, migration, and invasion, induced cell cycle arrest in the G1 phase, and promoted cell apoptosis. Furthermore, ZMIZ2 knockdown inhibited the mRNA and protein expression of CCL5, MCM3, E2F4, and DHX38.

Conclusion

Our findings reveal that ZMIZ2 is upregulated in TNBC tissues and is associated with its poor prognosis. ZMIZ2 may promote TNBC progression by promoting the expression of its target genes and affecting the corresponding pathways. Consequently, ZMIZ2 may serve as a promising target for future TNBC treatments.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12935-021-02393-x.

Retinoblastoma-E2F Transcription Factor Interplay Is Essential for Testicular Development and Male Fertility

The retinoblastoma (RB) protein family members (pRB, p107 and p130) are key regulators of cell cycle progression, but also play crucial roles in apoptosis, and stem cell self-renewal and differentiation. RB proteins exert their effects through binding to E2F transcription factors, which are essential developmental and physiological regulators of tissue and organ homeostasis. According to the canonical view, phosphorylation of RB results in release of E2Fs and induction of genes needed for progress of the cell cycle. However, there are eight members in the E2F transcription factor family with both activator (E2F1-3a) and repressor (E2F3b-E2F8) roles, highlighting the functional diversity of RB-E2F pathway. In this review article we summarize the data showing that RB-E2F interaction is a key cell-autonomous mechanism responsible for establishment and maintenance of lifelong male fertility. We also review the expression pattern of RB proteins and E2F transcription factors in the testis and male germ cells. The available evidence supports that RB and E2F family members are widely and dynamically expressed in the testis, and they are known to have versatile roles during spermatogenesis. Knowledge of the function and significance of RB-E2F interplay for testicular development and spermatogenesis comes primarily from gene knock-out (KO) studies. Several studies conducted in Sertoli cell-specific pRB-KO mice have demonstrated that pRB-mediated inhibition of E2F3 is essential for Sertoli cell functional maturation and cell cycle exit, highlighting that RB-E2F interaction in Sertoli cells is paramount to male fertility. Similarly, ablation of either pRB or E2F1 in the germline results in progressive testicular atrophy due to germline stem cell (GSC) depletion, emphasizing the importance of proper RB-E2F interplay for germline maintenance and lifelong sperm production. In summary, while balanced RB-E2F interplay is essential for cell-autonomous maintenance of GSCs and, the pRB-E2F3 system in Sertoli cells is critical for providing GSC niche thus laying the basis for spermatogenesis.

E2F4-Based Gene Therapy Mitigates the Phenotype of the Alzheimer's Disease Mouse Model 5xFAD

After decades of unfruitful work, no effective therapies are available for Alzheimer's disease (AD), likely due to its complex etiology that requires a multifactorial therapeutic approach. We have recently shown using transgenic mice that E2 factor 4 (E2F4), a transcription factor that regulates cell quiescence and tissue homeostasis, and controls gene networks affected in AD, represents a good candidate for a multifactorial targeting of AD. Here we show that the expression of a dominant negative form of human E2F4 (hE2F4DN), unable to become phosphorylated in a Thr-conserved motif known to modulate E2F4 activity, is an effective and safe AD multifactorial therapeutic agent. Neuronal expression of hE2F4DN in homozygous 5xFAD (h5xFAD) mice after systemic administration of an AAV.PHP.B-hSyn1.hE2F4DN vector reduced the production and accumulation of Aβ in the hippocampus, attenuated reactive astrocytosis and microgliosis, abolished neuronal tetraploidization, and prevented cognitive impairment evaluated by Y-maze and Morris water maze, without triggering side effects. This treatment also reversed other alterations observed in h5xFAD mice such as paw-clasping behavior and body weight loss. Our results indicate that E2F4DN-based gene therapy is a promising therapeutic approach against AD.

Hydroxyurea and Caffeine Impact pRb-like Protein-Dependent Chromatin Architecture Profiles in Interphase Cells of Vicia faba

The survival of cells depends on their ability to replicate correctly genetic material. Cells exposed to replication stress can experience a number of problems that may lead to deregulated proliferation, the development of cancer, and/or programmed cell death. In this article, we have induced prolonged replication arrest via hydroxyurea (HU) treatment and also premature chromosome condensation (PCC) by co-treatment with HU and caffeine (CF) in the root meristem cells of Vicia faba. We have analyzed the changes in the activities of retinoblastoma-like protein (RbS807/811ph). Results obtained from the immunocytochemical detection of RbS807/811ph allowed us to distinguish five unique activity profiles of pRb. We have also performed detailed 3D modeling using Blender 2.9.1., based on the original data and some final conclusions. 3D models helped us to visualize better the events occurring within the nuclei and acted as a high-resolution aid for presenting the results. We have found that, despite the decrease in pRb activity, its activity profiles were mostly intact and clearly recognizable, with some local alterations that may correspond to the increased demand in transcriptional activity. Our findings suggest that Vicia faba’s ability to withstand harsh environments may come from its well-developed and highly effective response to replication stress.

Dynamic effects of genetic variation on gene expression revealed following hypoxic stress in cardiomyocytes

One life-threatening outcome of cardiovascular disease is myocardial infarction, where cardiomyocytes are deprived of oxygen. To study inter-individual differences in response to hypoxia, we established an in vitro model of induced pluripotent stem cell-derived cardiomyocytes from 15 individuals. We measured gene expression levels, chromatin accessibility, and methylation levels in four culturing conditions that correspond to normoxia, hypoxia, and short- or long-term re-oxygenation. We characterized thousands of gene regulatory changes as the cells transition between conditions. Using available genotypes, we identified 1,573 genes with a cis expression quantitative locus (eQTL) in at least one condition, as well as 367 dynamic eQTLs, which are classified as eQTLs in at least one, but not in all conditions. A subset of genes with dynamic eQTLs is associated with complex traits and disease. Our data demonstrate how dynamic genetic effects on gene expression, which are likely relevant for disease, can be uncovered under stress.

Gene co-expression in the interactome: moving from correlation toward causation via an integrated approach to disease module discovery

In this study, we integrate the outcomes of co-expression network analysis with the human interactome network to predict novel putative disease genes and modules. We first apply the SWItch Miner (SWIM) methodology, which predicts important (switch) genes within the co-expression network that regulate disease state transitions, then map them to the human protein-protein interaction network (PPI, or interactome) to predict novel disease-disease relationships (i.e., a SWIM-informed diseasome). Although the relevance of switch genes to an observed phenotype has been recently assessed, their performance at the system or network level constitutes a new, potentially fascinating territory yet to be explored. Quantifying the interplay between switch genes and human diseases in the interactome network, we found that switch genes associated with specific disorders are closer to each other than to other nodes in the network, and tend to form localized connected subnetworks. These subnetworks overlap between similar diseases and are situated in different neighborhoods for pathologically distinct phenotypes, consistent with the well-known topological proximity property of disease genes. These findings allow us to demonstrate how SWIM-based correlation network analysis can serve as a useful tool for efficient screening of potentially new disease gene associations. When integrated with an interactome-based network analysis, it not only identifies novel candidate disease genes, but also may offer testable hypotheses by which to elucidate the molecular underpinnings of human disease and reveal commonalities between seemingly unrelated diseases.

LncRNA MYLK-AS1 facilitates tumor progression and angiogenesis by targeting miR-424-5p/E2F7 axis and activating VEGFR-2 signaling pathway in hepatocellular carcinoma

BACKGROUND

Long non-coding RNAs (lncRNAs) are crucial in the invasion, angiogenesis, progression, and metastasis of hepatocellular carcinoma (HCC). The lncRNA MYLK-AS1 promotes the growth and invasion of HCC through the EGFR/HER2-ERK1/2 signaling pathway. However, the clinical significance of MYLK-AS1 in HCC still needs to be further determined.

METHODS

Bioinformatic analysis was performed to determine the potential relationship among MYLK-AS1, miRNAs and mRNAs. A total of 156 samples of normal liver and paired HCC tissues from HCC patients were used to evaluate MYLK-AS1 expression by qRT-PCR. Human HCC cell lines were used to evaluate the colony formation, cell proliferation, migration, invasion, cell cycle and apoptosis after transfection of lentiviral short-hairpin RNAs (shRNAs) targeting MYLK-AS1 or MYLK-AS1 vectors. The competitive endogenous RNA (ceRNA) mechanism was clarified using fluorescence in situ hybridization (FISH), Western blotting, qPCR, RNA binding protein immunoprecipitation (RIP), and dual luciferase reporter analysis.

RESULTS

MYLK-AS1 up-regulation was detected in the HCC tumor tissues and cell lines associated with the enhancement of the angiogenesis and tumor progression. The down-regulation of MYLK-AS1 reversed the effects on angiogenesis, proliferation, invasion and metastasis in the HCC cells and in vivo. MYLK-AS1 acted as ceRNA, capable of regulating the angiogenesis in HCC, while the microRNA miR-424-5p was the direct target of MYLK-AS1. Promoting the angiogenesis and the tumor proliferation, the complex MYLK-AS1/miR-424-5p activated the VEGFR-2 signaling through E2F7, whereas the specific targeting of E2F transcription factor 7 (E2F7) by miR-424-5p, was indicated by the mechanism studies.

CONCLUSIONS

MYLK-AS1 and E2F7 are closely related to some malignant clinicopathological features and prognosis of HCC, thus the MYLK-AS1/ miR-424-5p/E2F7 signaling pathway might represent a promising treatment strategy to combat HCC.

E2F7-EZH2 axis regulates PTEN/AKT/mTOR signalling and glioblastoma progression

BACKGROUND

E2F transcription factors are considered to be important drivers of tumour growth. E2F7 is an atypical E2F factor, and its role in glioblastoma remains undefined.

METHODS

E2F7 expression was examined in patients by IHC and qRT-PCR. The overall survival probability was determined by statistical analyses. MTT assay, colony formation, cell-cycle assay, cell metastasis and the in vivo model were employed to determine the functional role of E2F7 in glioblastoma. Chromatin immunoprecipitation, luciferase assay and western blot were used to explore the underlying mechanisms.

RESULTS

E2F7 was found to be up-regulated in glioblastoma patients, and high E2F7 expression was associated with poor overall survival in glioblastoma patients. Functional studies showed that E2F7 promoted cell proliferation, cell-cycle progression, cell metastasis and tumorigenicity abilities in vitro and in vivo. E2F7 promoted the transcription of EZH2 by binding to its promoter and increased H3K27me3 level. EZH2 recruited H3K27me3 to the promoter of PTEN and inhibited PTEN expression, and then activated the AKT/mTOR signalling pathway. In addition, restored expression of EZH2 recovered the abilities of cell proliferation and metastasis in E2F7-silencing cells.

CONCLUSION

Collectively, our findings indicate that E2F7 promotes cell proliferation, cell metastasis and tumorigenesis via EZH2-mediated PTEN/AKT/mTOR pathway in glioblastoma.

E2F8 Induces Cell Proliferation and Invasion through the Epithelial-Mesenchymal Transition and Notch Signaling Pathways in Ovarian Cancer

Background: Despite the recent research implicating E2F8 (E2F Transcription Factor 8) in cancer, the role of E2F8 in the progression of ovarian cancer has remained unclear. Hence, we explored the bio-functional effects of E2F8 knockdown on ovarian cancer cell lines in vitro and in vivo. Methods: The expression of E2F8 was compared between ovarian cancer and noncancer tissues, and its association with the progression-free survival of ovarian cancer patients was analyzed. To demonstrate the function of E2F8 in cell proliferation, migration, and invasion, we employed RNA interference to suppress E2F8 expression in ovarian cancer cell lines. Finally, the effect of E2F8 knockdown was investigated in a xenograft mouse model of ovarian cancer. Results: Ovarian cancer tissue exhibited significantly higher E2F8 expression compared to that of normal ovarian tissue. Clinical data showed that E2F8 was a significant predictor of progression-free survival. Moreover, the prognosis of the ovarian cancer patients with high E2F8 expression was poorer than that of the patients with low E2F8 expression. In vitro experiments using E2F8-knockdown ovarian cancer cell lines demonstrated that E2F8 knockdown inhibited cell proliferation, migration, and tumor invasion. Additionally, E2F8 was a potent inducer and modulator of the expression of epithelial–mesenchymal transition and Notch signaling pathway-related markers. We confirmed the function of E2F8 in vivo, signifying that E2F8 knockdown was significantly correlated with reduced tumor size and weight. Conclusions: Our findings indicate that E2F8 is highly correlated with ovarian cancer progression. Hence, E2F8 can be utilized as a prognostic marker and therapeutic target against ovarian malignancy.

Human papillomavirus type 38 alters wild-type p53 activity to promote cell proliferation via the downregulation of integrin alpha 1 expression

Tumor suppressors can exert pro-proliferation functions in specific contexts. In the beta human papillomavirus type 38 (HPV38) experimental model, the viral proteins E6 and E7 promote accumulation of a wild-type (WT) p53 form in human keratinocytes (HKs), promoting cellular proliferation. Inactivation of p53 by different means strongly decreases the proliferation of HPV38 E6/E7 HKs. This p53 form is phosphorylated at S392 by the double-stranded RNA-dependent protein kinase PKR, which is highly activated by HPV38. PKR-mediated S392 p53 phosphorylation promotes the formation of a p53/DNMT1 complex, which inhibits expression of integrin alpha 1 (ITGA1), a repressor of epidermal growth factor receptor (EGFR) signaling. Ectopic expression of ITGA1 in HPV38 E6/E7 HKs promotes EGFR degradation, inhibition of cellular proliferation, and cellular death. Itga1 expression was also inhibited in the skin of HPV38 transgenic mice that have an elevated susceptibility to UV-induced skin carcinogenesis. In summary, these findings reveal the existence of a specific WT p53 form that displays pro-proliferation properties.

Beyond What Your Retina Can See: Similarities of Retinoblastoma Function between Plants and Animals, from Developmental Processes to Epigenetic Regulation

The Retinoblastoma protein (pRb) is a key cell cycle regulator conserved in a wide variety of organisms. Experimental analysis of pRb's functions in animals and plants has revealed that this protein participates in cell proliferation and differentiation processes. In addition, pRb in animals and its orthologs in plants (RBR), are part of highly conserved protein complexes which suggest the possibility that analogies exist not only between functions carried out by pRb orthologs themselves, but also in the structure and roles of the protein networks where these proteins are involved. Here, we present examples of pRb/RBR participation in cell cycle control, cell differentiation, and in the regulation of epigenetic changes and chromatin remodeling machinery, highlighting the similarities that exist between the composition of such networks in plants and animals.

Long non-coding RNA LINC00337 induces autophagy and chemoresistance to cisplatin in esophageal squamous cell carcinoma cells via upregulation of TPX2 by recruiting E2F4

Esophageal cancer represents the eighth most frequently occurring cancer, as well as the sixth most widespread cause of cancer-related deaths. In recent years, accumulating evidence has implicated long non-coding RNAs in the progression of esophageal squamous cell carcinoma (ESCC). The aim of the present study was to investigate the potential involvement and underlying mechanisms of LINC00337 in ESCC. Expression patterns of LINC00337 and targeting protein for Xenopus kinesin-like protein 2 (TPX2) in ESCC tissues and cells were detected using RT-qPCR and immunohistochemical staining. Next, the interactions among LINC00337, E2F4, and TPX2 were identified using chromatin immunoprecipitation, dual-luciferase reporter, and RNA-binding protein immunoprecipitation assays, suggesting that LINC00337 could recruit E2F4 to enhance the transcription of TPX2. Thereafter, the regulatory roles of LINC00337 and TPX2 in ESCC were analyzed by altering the expression of LINC00337 or TPX2 in ESCC cells following treatment with cisplatin (DDP). The levels of autophagy-related proteins Beclin1 and LC3II/LC3I, viability, autophagy, apoptosis, and chemoresistance of ESCC cells to DDP were measured following transfection treatment with different plasmids. Additionally, the role of the LINC00337/E2F4/TPX2 axis was assessed in vivo by measuring tumor formation in nude mice. The results demonstrated that LINC00337 upregulated TPX2, consequently leading to elevated levels of Beclin1 and LC3II/LC3I, promoted cell viability and autophagy, while inhibiting apoptosis and chemosensitivity to DDP in ESCC. In sum, the current study evidenced that the overexpression of LINC00337 could potentially enhance ESCC cell autophagy and chemoresistance to DDP via the upregulation of TPX2 by recruiting E2F4. Thus, LINC00337 may serve as a potential candidate for the treatment of ESCC.

Motor neurons from ALS patients with mutations in C9ORF72 and SOD1 exhibit distinct transcriptional landscapes

Amyotrophic lateral sclerosis (ALS) is a progressive motor neuron disease that culminates in paralysis and death. Here, we present our analyses of publicly available multiOMIC data sets generated using motor neurons from ALS patients and control cohorts. Functional annotation of differentially expressed genes in induced pluripotent stem cell (iPSC)-derived motor neurons generated from patients with mutations in C9ORF72 (C9-ALS) suggests elevated expression of genes that pertain to extracellular matrix (ECM) and cell adhesion, inflammation and TGFβ targets. On the other end of the continuum, we detected diminished expression of genes repressed by quiescence-promoting E2F4/DREAM complex. Proteins whose abundance was significantly altered in C9-ALS neurons faithfully recapitulated the transcriptional aberrations. Importantly, patterns of gene expression in spinal motor neurons dissected from C9-ALS or sporadic ALS patients were highly concordant with each other and with the C9-ALS iPSC neurons. In contrast, motor neurons from patients with mutations in SOD1 exhibited dramatically different signatures. Elevated expression of gene sets such as ECM and cell adhesion genes occurs in C9 and sporadic ALS but not SOD1-ALS. These analyses indicate that despite the similarities in outward manifestations, transcriptional and proteomic signatures in ALS motor neurons can vary significantly depending on the identity of the causal mutations.

Leptin promotes methionine adenosyltransferase 2A expression in hepatic stellate cells by the downregulation of E2F-4 via the β-catenin pathway

Most obese patients develop hyperleptinaemia. Leptin, mainly produced by adipocytes, demonstrates a promotional role in liver fibrosis. Hepatic stellate cell (HSC) activation, a key step in liver fibrogenesis, requires global reprogramming of gene expression. The remodeling of DNA methylation is a mechanism of the epigenetic regulation of gene expression. The biosynthesis of S-adenosylmethionine, a principle biological methyl donor, is catalyzed by methionine adenosyltransferase (MAT) such as MATⅡ which has been shown to promote HSC activation in vitro. This study was mainly aimed to determine the effect of leptin on MAT2A expression (the catalytic subunit of MATⅡ) in HSCs. Results showed that MAT2A knockdown reduced leptin-induced HSC activation and liver fibrosis in the leptin-deficient mouse model. Leptin promoted MAT2A expression in HSCs and increased MAT2A promoter activity. The axis of the β-catenin pathway/E2F-4 mediated the effect of leptin on MAT2A expression. Leptin-induced β-catenin signaling reduced E2F-4 expression and thus abated E2F-4 binding to MAT2A promoter at a site around -2779 bp, leading to an increase in the MAT2A promoter activity. These data might shed more light on the mechanisms responsible for liver fibrogenesis in obese patients with hyperleptinaemia.

Integrative prediction of gene expression with chromatin accessibility and conformation data

BACKGROUND

Enhancers play a fundamental role in orchestrating cell state and development. Although several methods have been developed to identify enhancers, linking them to their target genes is still an open problem. Several theories have been proposed on the functional mechanisms of enhancers, which triggered the development of various methods to infer promoter-enhancer interactions (PEIs). The advancement of high-throughput techniques describing the three-dimensional organization of the chromatin, paved the way to pinpoint long-range PEIs. Here we investigated whether including PEIs in computational models for the prediction of gene expression improves performance and interpretability.

RESULTS

We have extended our [Formula: see text] framework to include DNA contacts deduced from chromatin conformation capture experiments and compared various methods to determine PEIs using predictive modelling of gene expression from chromatin accessibility data and predicted transcription factor (TF) motif data. We designed a novel machine learning approach that allows the prioritization of TFs binding to distal loop and promoter regions with respect to their importance for gene expression regulation. Our analysis revealed a set of core TFs that are part of enhancer-promoter loops involving YY1 in different cell lines.

CONCLUSION

We present a novel approach that can be used to prioritize TFs involved in distal and promoter-proximal regulatory events by integrating chromatin accessibility, conformation, and gene expression data. We show that the integration of chromatin conformation data can improve gene expression prediction and aids model interpretability.

PIR promotes tumorigenesis of breast cancer by upregulating cell cycle activator E2F1

Pirin (PIR) protein belongs to the superfamily of cupin and is highly conserved between eukaryotic and prokaryotic organisms. It has been reported that PIR is upregulated in various tumors and involved in tumorigenesis. However, its biological functions particularly in promoting tumorigenesis are, to date, poorly characterized. Here we report that knockdown of PIR in MCF7 and MDA-MB-231 cell lines causes a dramatic decrease in cell proliferation and xenograft tumor growth in mice. Mechanistically, the cell cycle activator E2F1 and its target genes cdk4, cdk6, cycE, cycD and DDR1 are remarkably downregulated in PIR depleted cells, leading to G1/S phase arrest. Luciferase reporter assay and chromatin immunoprecipitation assay indicate that PIR can activate E2F1 transcription by binding to its promoter region. Consistent with the observation in PIR knockdown cells, PIR inhibitors markedly inhibit the proliferation of both cell lines. Furthermore, knockdown of PIR significantly decreases the abilities of MCF7 cells for mobility and invasion in vitro and their metastasis in mice, which may be attributed to the decrease of DDR1. In conclusion, PIR stimulates tumorigenesis and progression by activating E2F1 and its target genes. Our finding thus suggests PIR as a potential druggable target for the therapy of cancers with high expression level of PIR.

E2F6 is essential for cell viability in breast cancer cells during replication stress

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.

ZC3H18 specifically binds and activates the BRCA1 promoter to facilitate homologous recombination in ovarian cancer

Reduced BRCA1 expression causes homologous recombination (HR) repair defects in high-grade serous ovarian cancers (HGSOCs). Here, we demonstrate that BRCA1 is transcriptionally activated by a previously unknown function of ZC3H18. We show that ZC3H18 is a DNA-binding protein that interacts with an E2F site in the BRCA1 promoter where it facilitates recruitment of E2F4 to an adjacent E2F site to promote BRCA1 transcription. Consistent with ZC3H18 role in activating BRCA1 expression, ZC3H18 depletion induces BRCA1 promoter methylation, reduces BRCA1 expression, disrupts HR, and sensitizes cells to DNA crosslinkers and poly(ADP-ribose) polymerase inhibitors. Moreover, in patient-derived xenografts and primary HGSOC tumors, ZC3H18 and E2F4 mRNA levels are positively correlated with BRCA1 mRNA levels, further supporting ZC3H18 role in regulating BRCA1. Given that ZC3H18 lies within 16q24.2, a region with frequent copy number loss in HGSOC, these findings suggest that ZC3H18 copy number losses could contribute to HR defects in HGSOC.

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.

Integrative characterization of G-Quadruplexes in the three-dimensional chromatin structure

DNA molecules are highly compacted in the eukaryotic nucleus where distal regulatory elements reach their targets through three-dimensional chromosomal interactions. G-quadruplexes, stable four-stranded non-canonical DNA structures, can change local chromatin organization through the exclusion of nucleosomes. However, the relationship between G-quadruplexes and higher-order genome organization remains unknown. Here, we found that G-quadruplexes are significantly enriched at boundaries of topological associated domains (TADs). Architectural protein occupancy, which plays critical roles in the formation of TADs, was highly correlated with the content of G-quadruplexes at TAD boundaries. Moreover, adjacent boundaries containing G-quadruplexes frequently interacted with each other because of the high enrichment of architectural protein binding sites. Similar to CCCTC-binding factor (CTCF) binding sites, G-quadruplexes also showed strong insulation ability in the separation of adjacent regions. Additionally, the insulation ability of CTCF binding sites and TAD boundaries was significantly reinforced by G-quadruplexes. Furthermore, G-quadruplex motifs on different strands were associated with the orientation of CTCF binding sites. These findings suggest a potential role for G-quadruplexes in loop extrusion. The enrichment of transcription factor binding sites (TFBSs) around regulatory elements containing G-quadruplexes led to frequent interactions between regulatory elements containing G-quadruplexes. Intriguingly, more than 99% of G-quadruplexes overlapped with TFBSs. The binding sites of CTCF and cohesin proteins were preferentially located surrounding G-quadruplexes. Accordingly, we proposed a new mechanism of long-distance gene regulation in which G-quadruplexes are involved in distal interactions between enhancers and promoters.

Healthy diet is associated with gene expression in blood: the Framingham Heart Study

BACKGROUND

Genes in metabolic and nutrient signaling pathways play important roles in lifespan in model organisms and human longevity.

OBJECTIVE

The aim of this study was to examine the relation of a quantitative measure of healthy diet to gene expression in a community-based cohort.

METHODS

We used the 2015 Dietary Guidelines for Americans Adherence Index (DGAI) score to quantify key dietary recommendations of an overall healthy diet. Our current analyses included 2220 Offspring participants (mean age 66 ± 9 y, 55.4% women) and 2941 Third-Generation participants (mean age 46 ± 9 y, 54.5% women) from the Framingham Heart Study. Gene expression was profiled in blood through the use of the Affymetrix Human Exon 1.0 ST Array. We conducted a transcriptome-wide association study of DGAI adjusting for age, sex, smoking, cell counts, and technical covariates. We also constructed a combined gene score from genes significantly associated with DGAI.

RESULTS

The DGAI was significantly associated with the expression of 19 genes (false discovery rate <0.05). The most significant gene, ARRDC3, is a member of the arrestin family of proteins, and evidence in animal models and human data suggests that this gene is a regulator of obesity and energy expenditure. The DGAI gene score was associated with body mass index (P = 1.4 × 10-50), fasting glucose concentration (P = 2.5 × 10-11), type 2 diabetes (P = 1.1 × 10-5), and metabolic syndrome (P = 1.8 × 10-32).

CONCLUSIONS

Healthier diet was associated with genes involved in metabolic function. Further work is needed to replicate our findings and investigate the relation of a healthy diet to altered gene regulation.

E2F4 regulates transcriptional activation in mouse embryonic stem cells independently of the RB family

E2F transcription factors are central regulators of cell division and cell fate decisions. E2F4 often represents the predominant E2F activity in cells. E2F4 is a transcriptional repressor implicated in cell cycle arrest and whose repressive activity depends on its interaction with members of the RB family. Here we show that E2F4 is important for the proliferation and the survival of mouse embryonic stem cells. In these cells, E2F4 acts in part as a transcriptional activator that promotes the expression of cell cycle genes. This role for E2F4 is independent of the RB family. Furthermore, E2F4 functionally interacts with chromatin regulators associated with gene activation and we observed decreased histone acetylation at the promoters of cell cycle genes and E2F targets upon loss of E2F4 in RB family-mutant cells. Taken together, our findings uncover a non-canonical role for E2F4 that provide insights into the biology of rapidly dividing cells.

E2F transcription factors are regulators of cell division and cell fate decisions. Here the authors show that E2F4 is important for proliferation and survival of mouse ESCs, independent of the RB family, and that E2F4 interacts with chromatin regulators associated with gene activation.

The Toxoplasma effector TEEGR promotes parasite persistence by modulating NF-κB signalling via EZH2

The protozoan parasite Toxoplasma gondii has co-evolved with its homeothermic hosts (humans included) strategies that drive its quasi-asymptomatic persistence in hosts, hence optimizing the chance of transmission to new hosts. Persistence, which starts with a small subset of parasites that escape host immune killing and colonize the so-called immune privileged tissues where they differentiate into a low replicating stage, is driven by the interleukin 12 (IL-12)-interferon-γ (IFN-γ) axis. Recent characterization of a family of Toxoplasma effectors that are delivered into the host cell, in which they rewire the host cell gene expression, has allowed the identification of regulators of the IL-12-IFN-γ axis, including repressors. We now report on the dense granule-resident effector, called TEEGR (Toxoplasma E2F4-associated EZH2-inducing gene regulator) that counteracts the nuclear factor-κB (NF-κB) signalling pathway. Once exported into the host cell, TEEGR ends up in the nucleus where it not only complexes with the E2F3 and E2F4 host transcription factors to induce gene expression, but also promotes shaping of a non-permissive chromatin through its capacity to switch on EZH2. Remarkably, EZH2 fosters the epigenetic silencing of a subset of NF-κB-regulated cytokines, thereby strongly contributing to the host immune equilibrium that influences the host immune response and promotes parasite persistence in mice.

TREX1 suppression imparts cancer-stem-cell-like characteristics to CD133- osteosarcoma cells through the activation of E2F4 signaling

There is ongoing debate whether cancer stem cells (CSCs) could arise from the transformation of non-CSCs under specific conditions. In the present study, the role of the three prime repair exonuclease 1 (TREX1) in regulating CSC generation form human osteosarcoma cells was investigated. High, intermediate and low levels of TREX1 expression were respectively observed in low-grade, high-grade and metastatic human osteosarcoma samples, while the opposite tendency was observed for E2F4, a transcription factor associated with G2 arrest. Luciferase assay proved that TREX1 had a negative impact on the activity of E2F4 promoter. TREX1 was highly expressed in CD133^- HOS cells (non-CSC osteosarcoma cells) compared to CD133⁺ ones; whereas TREX1 knockdown endowed the CD133^- non-CSCs with CSC-like characteristics in vitro relying on E2F4 activation, as demonstrated by enlarged proportion of the subset expressing CSC markers in flow cytometry analysis, enhanced self-renewal ability in osteosphere formation assay, increased metastasis capacity in migration and invasion assays, together with improved chemoresistance to cisplatin. Furthermore, TREX1 knockdown and subsequent E2F4 activation could promote the tumorigenicity of CD133^- non-CSCs in vivo. With respect to underlying mechanisms, it was found that in CD133^- HOS cells, TREX1 suppression would allow the activation of β-catenin signaling in the dependence of E2F4, thus possibly leading to the up-regulation of the transcription factor OCT4. These findings suggested that TREX1 was probably a negative regulator of CSC formation and hence worth to be further studied for developing new treatments in cancer therapies targeting CSCs.

Identification of human silencers by correlating cross-tissue epigenetic profiles and gene expression

Compared to enhancers, silencers are notably difficult to identify and validate experimentally. In search for human silencers, we utilized H3K27me3-DNase I hypersensitive site (DHS) peaks with tissue specificity negatively correlated with the expression of nearby genes across 25 diverse cell lines. These regions are predicted to be silencers since they are physically linked, using Hi-C loops, or associated, using expression quantitative trait loci (eQTL) results, with a decrease in gene expression much more frequently than general H3K27me3-DHSs. Also, these regions are enriched for the binding sites of transcriptional repressors (such as CTCF, MECOM, SMAD4, and SNAI3) and depleted of the binding sites of transcriptional activators. Using sequence signatures of these regions, we constructed a computational model and predicted approximately 10,000 additional silencers per cell line and demonstrated that the majority of genes linked to these silencers are expressed at a decreased level. Furthermore, single nucleotide polymorphisms (SNPs) in predicted silencers are significantly associated with disease phenotypes. Finally, our results show that silencers commonly interact with enhancers to affect the transcriptional dynamics of tissue-specific genes and to facilitate fine-tuning of transcription in the human genome.

Impact of dynamin 2 on adenovirus nuclear entry

The large GTPase dynamin 2 controls both endosomal fission and microtubule acetylation. Here we report that dynamin 2 alters microtubules and regulates the trafficking of human adenovirus type 37. Dynamin 2 knockdown by siRNA in infected cells resulted in accumulation of acetylated tubulin, repositioning of microtubule organizing centers (MTOCs) closer to cell nuclei, increased virus in the cytosol (with a compensatory decrease in endosomal virus), reduced proinflammatory cytokine induction, and increased binding of virus to the nucleoporin, Nup358. These events led to increased viral DNA nuclear entry and viral replication. Overexpression of dynamin 2 generated opposite effects. Therefore, dynamin 2 inhibits adenovirus replication and promotes innate immune responses by the infected cell. MTOC transposition in dynamin 2 knockdown promotes a closer association with nuclear pore complexes to facilitate viral DNA delivery. Dynamin 2 plays a key role in adenoviral trafficking and influences host responses to infection.

Dichotomy in redundant enhancers points to presence of initiators of gene regulation

Background

The regulatory landscape of a gene locus often consists of several functionally redundant enhancers establishing phenotypic robustness and evolutionary stability of its regulatory program. However, it is unclear what mechanisms are employed by redundant enhancers to cooperatively orchestrate gene expression.

Results

By comparing redundant enhancers to single enhancers (enhancers present in a single copy in a gene locus), we observed that the DNA sequence encryption differs between these two classes of enhancers, suggesting a difference in their regulatory mechanisms. Initiator enhancers, which are a subset of redundant enhancers and show similar sequence encryption to single enhancers, differ from the rest of redundant enhancers in their sequence encryption, evolutionary conservation and proximity to target genes. Genes hosting initiator enhancers in their loci feature elevated levels of expression. Initiator enhancers show a high level of 3D chromatin contacts with both transcription start sites and regular enhancers, suggesting their roles as primary activators and intermediate catalysts of gene expression, through which the regulatory signals of redundant enhancers are propagated to the target genes. In addition, GWAS and eQTLs variants are significantly enriched in initiator enhancers compared to redundant enhancers, suggesting a key functional role these sequences play in gene regulation.

Conclusions

The specific characteristics and widespread abundance of initiator enhancers advocate for a possible universal hierarchical mechanism of tissue-specific gene regulation involving multiple redundant enhancers acting through initiator enhancers.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5335-0) contains supplementary material, which is available to authorized users.

Activation of E2F-dependent transcription by the mouse cytomegalovirus M117 protein affects the viral host range

Cytomegaloviruses (CMVs) have a highly restricted host range as they replicate only in cells of their own or closely related species. To date, the molecular mechanisms underlying the CMV host restriction remain poorly understood. However, it has been shown that mouse cytomegalovirus (MCMV) can be adapted to human cells and that adaptation goes along with adaptive mutations in several viral genes. In this study, we identify MCMV M117 as a novel host range determinant. Mutations in this gene enable the virus to cross the species barrier and replicate in human RPE-1 cells. We show that the M117 protein is expressed with early kinetics, localizes to viral replication compartments, and contributes to the inhibition of cellular DNA synthesis. Mechanistically, M117 interacts with members of the E2F transcription factor family and induces E2F target gene expression in murine and human cells. While the N-terminal part of M117 mediates E2F interaction, the C-terminal part mediates self-interaction. Both parts are required for the activation of E2F-dependent transcription. We further show that M117 is dispensable for viral replication in cultured mouse fibroblasts and endothelial cells, but is required for colonization of mouse salivary glands in vivo. Conversely, inactivation of M117 or pharmacological inhibition of E2F facilitates MCMV replication in human RPE-1 cells, whereas replacement of M117 by adenovirus E4orf6/7, a known E2F activator, prevents it. These results indicate that E2F activation is detrimental for MCMV replication in human cells. In summary, this study identifies MCMV M117 as a novel E2F activator that functions as a host range determinant by precluding MCMV replication in human cells.

An altered gene expression profile in tyramine-exposed intestinal cell cultures supports the genotoxicity of this biogenic amine at dietary concentrations

Tyramine, histamine and putrescine are the most commonly detected and most abundant biogenic amines (BA) in food. The consumption of food with high concentrations of these BA is discouraged by the main food safety agencies, but legal limits have only been set for histamine. The present work reports a transcriptomic investigation of the oncogenic potential of the above-mentioned BA, as assessed in the HT29 human intestinal epithelial cell line. Tyramine had a greater effect on the expression of genes involved in tumorigenesis than did histamine or putrescine. Since some of the genes that showed altered expression in tyramine-exposed cells are involved in DNA damage and repair, the effect of this BA on the expression of other genes involved in the DNA damage response was investigated. The results suggest that tyramine might be genotoxic for intestinal cells at concentrations easily found in BA-rich food. Moreover, a role in promoting intestinal cancer cannot be excluded.

Identification of a functional antioxidant response element at the HIF1A locus

Reactive oxygen species (ROS), which are a byproduct of oxidative metabolism, serve as signaling molecules in a number of physiological settings. However, if their levels are not tightly maintained, excess ROS lead to potentially cytotoxic oxidative stress. Accordingly, several transcriptional regulatory networks have evolved to include components that are highly ROS-responsive. Depending on the context, these regulatory networks can leverage ROS to respond to nutrient conditions, metabolism, or other physiological signals, or to respond to oxidative stress. However, ROS signaling is complex, so regulatory interactions between various ROS-responsive transcription factors are still being mapped out. Here we show that the transcription factor NRF2, a key regulator of the adaptive response to oxidative stress, directly regulates expression of HIF1A, which encodes HIF1α, a key transcriptional regulator of the adaptive response to hypoxia. We used an integrative genomics approach to identify HIF1A as a ROS-responsive transcript and we found an NRF2-bound antioxidant response element (ARE) approximately 30 kilobases upstream of HIF1A. This ARE sequence is deeply conserved, and we verified that it is directly bound and activated by NRF2. In addition, we found that HIF1A is upregulated in breast and bladder tumors with high NRF2 activity. Taken together, our results demonstrate that NRF2 targets a functional ARE at the HIF1A locus, and reveal a direct regulatory connection between two important oxygen responsive transcription factors.

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.

Functional haplotypes of ARID4A affect promoter activity and semen quality of bulls

The AT-rich interaction domain 4 A (ARID4A) has an important role in regulating Sertoli cell function and male fertility. Its molecular mechanisms, however, remain largely unknown. In this study, two single nucleotide polymorphisms (SNPs) (g.53 G > T, ss 1966531596, and g.826 G > A, rs 210809648) were identified in the promoter region of ARID4A in 215 Chinese Holstein bulls using polymerase chain reaction (PCR)-restriction fragment length polymorphism and created restriction site-PCR. Results revealed that bulls with g.53 G > T-GG and g.826 G > A-G G genotype exhibited higher sperm deformity rate than those with g.53 G > T-TT and g.826 G > A-AA genotype (P <  0.01). Furthermore, three haplotypes (H1 (GG), H3 (TG), H4 (TA)) and six haplotype combinations (H1H1, H1H3, H1H4, H3H3, H3H4, H4H4) were obtained. The bulls with H4H4 exhibited lower sperm deformity rate than those with H1H1 and H1H3 (P <  0.05). In addition, results of bioinformatics analysis revealed that ARID4A has two promoters and that two SNPs of ARID4A are located in transcription factor binding sites. Compared with g.53 G > T-G and g.826 G > A-G allele, there was a greater fluorescence intensity in g.53 G > T-T and g.826 G > A-A allele by transient transfection in MLTC-1 cells and the luciferase report assay. qRT-PCR indicated the ARID4A expression was greater in bull spermatozoa with H4H4 haplotype combination than those with H1H1 haplotype combination (P < 0.05). Results of the present study indicate that g.53 G > T and g.826 G > A are functional mutations that are involved in regulation of ARID4A gene expression by affecting promoter activity and thus semen quality of Chinese Holstein bulls.

JQ1 affects BRD2-dependent and independent transcription regulation without disrupting H4-hyperacetylated chromatin states

The bromodomain and extra-terminal domain (BET) proteins are promising drug targets for cancer and immune diseases. However, BET inhibition effects have been studied more in the context of bromodomain-containing protein 4 (BRD4) than BRD2, and the BET protein association to histone H4-hyperacetylated chromatin is not understood at the genome-wide level. Here, we report transcription start site (TSS)-resolution integrative analyses of ChIP-seq and transcriptome profiles in human non-small cell lung cancer (NSCLC) cell line H23. We show that di-acetylation at K5 and K8 of histone H4 (H4K5acK8ac) co-localizes with H3K27ac and BRD2 in the majority of active enhancers and promoters, where BRD2 has a stronger association with H4K5acK8ac than H3K27ac. Although BET inhibition by JQ1 led to complete reduction of BRD2 binding to chromatin, only local changes of H4K5acK8ac levels were observed, suggesting that recruitment of BRD2 does not influence global histone H4 hyperacetylation levels. This finding supports a model in which recruitment of BET proteins via histone H4 hyperacetylation is predominant over hyperacetylation of histone H4 by BET protein-associated acetyltransferases. In addition, we found that a remarkable number of BRD2-bound genes, including MYC and its downstream target genes, were transcriptionally upregulated upon JQ1 treatment. Using BRD2-enriched sites and transcriptional activity analysis, we identified candidate transcription factors potentially involved in the JQ1 response in BRD2-dependent and -independent manner.

Potential biomarkers and therapeutic targets in cervical cancer: Insights from the meta-analysis of transcriptomics data within network biomedicine perspective

The malignant neoplasm of the cervix, cervical cancer, has effects on the reproductive tract. Although infection with oncogenic human papillomavirus is essential for cervical cancer development, it alone is insufficient to explain the development of cervical cancer. Therefore, other risk factors such as host genetic factors should be identified, and their importance in cervical cancer induction should be determined. Although gene expression profiling studies in the last decade have made significant molecular findings about cervical cancer, adequate screening and effective treatment strategies have yet to be achieved. In the current study, meta-analysis was performed on cervical cancer-associated transcriptome data and reporter biomolecules were identified at RNA (mRNA, miRNA), protein (receptor, transcription factor, etc.), and metabolite levels by the integration of gene expression profiles with genome-scale biomolecular networks. This approach revealed already-known biomarkers, tumor suppressors and oncogenes in cervical cancer as well as various receptors (e.g. ephrin receptors EPHA4, EPHA5, and EPHB2; endothelin receptors EDNRA and EDNRB; nuclear receptors NCOA3, NR2C1, and NR2C2), miRNAs (e.g., miR-192-5p, miR-193b-3p, and miR-215-5p), transcription factors (particularly E2F4, ETS1, and CUTL1), other proteins (e.g., KAT2B, PARP1, CDK1, GSK3B, WNK1, and CRYAB), and metabolites (particularly, arachidonic acids) as novel biomarker candidates and potential therapeutic targets. The differential expression profiles of all reporter biomolecules were cross-validated in independent RNA-Seq and miRNA-Seq datasets, and the prognostic power of several reporter biomolecules, including KAT2B, PCNA, CD86, miR-192-5p and miR-215-5p was also demonstrated. In this study, we reported valuable data for further experimental and clinical efforts, because the proposed biomolecules have significant potential as systems biomarkers for screening or therapeutic purposes in cervical carcinoma.