E2F - at the crossroads of life and death

The evolution and functional characterization of transcription factors E2Fs in lamprey, Lethenteron reissneri

The E2 promoter binding factors (E2Fs) are a group of transcriptional regulators that govern the cell cycle and play crucial roles in various cellular physiological processes, including proliferation and embryonic development. In this study, we identified four homologous genes-Lr-E2F3, Lr-E2F4, Lr-E2F5, and Lr-E2F8-from the lamprey (Lethenteron reissneri) genome database. Phylogenetic tree analysis was conducted to elucidate the evolutionary relationships within the E2F family across different species. Furthermore, analyses of motifs, domains, gene structures, and 3D structures reinforced the conservation of the E2F family. Notably, synteny analysis revealed that the neighboring genes of the Lr-E2Fs exhibited greater diversity compared to those in jawed vertebrates. Activity assays indicated that Lr-E2Fs may be involved in lamprey innate immunity mediated by NF-кB. Additionally, morphological observations of embryos microinjected with Cas9/sgRNA demonstrated that E2F-deficient lamprey embryos displayed embryonic lethality, suggesting that Lr-E2Fs play a significant role in lamprey embryonic development. In summary, our research not only provides new insights into the evolution of Lr-E2Fs but also offers valuable clues regarding their functional roles.

E2F1 Promotes the Occurrence of Head and Neck Squamous Cell Carcinoma and Serves as a Prognostic Biomarker

Head and neck squamous cell carcinoma (HNSCC) is a common malignant tumor occurring in various sites such as the oral cavity, pharynx, larynx, and nasal cavity. This study aimed to explore the biological functions and prognostic value of E2F transcription factor 1 (E2F1) in HNSCC. Transcriptome and single-cell sequencing (scRNA-seq) data of HNSCC patients were analyzed using data from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO). All samples were divided into high and low expression groups based on the expression levels of E2F1. A risk model was constructed based on Lasso-Cox regression, and the differences between the two groups in terms of prognosis were explored. The scRNA-seq data of HNSCC samples were analyzed using the Seurat package to identify cell types. AUCell was used to score different types of cells, and subsequently, the interaction pathways between the high-scoring cell population and other cell populations were explored using the CellChat package. The expression level of E2F1 in tumor tissues was higher than that in normal tissues, which was confirmed by in vitro experiments. Analysis of transcriptome data from TCGA revealed significant differences in overall survival (OS) between the high and low expression groups. Prognostic genes were selected based on DEGs between the two groups, and a risk model was constructed. Subsequently, a nomogram model was constructed based on clinical factors and risk scores, which exhibited good predictive performance. The expression landscape of prognostic genes in different cell types was explored using scRNA-seq data of HNSCC samples. Dendritic cell populations were identified as high-scoring cell populations, and the pathways of interaction between this cell population and other cell populations were explored. We identified E2F1 as an independent prognostic factor closely associated with the prognosis and immune response of HNSCC.

Revisiting the female germline cell development

The formation of the female germline is the fundamental process in most flowering plants' sexual reproduction. In Arabidopsis, only one somatic cell obtains the female germline fate, and this process is regulated by different pathways. Megaspore mother cell (MMC) is the first female germline, and understanding MMC development is essential for comprehending the complex mechanisms of plant reproduction processes. Recently, more advanced technologies such as whole-mount single-molecule fluorescence in situ hybridization (smFISH), laser-assisted microdissection (LCM), chromatin immunoprecipitation/sequencing, and CRISPR gene editing have provided opportunities to reveal the mechanism of female germline development at different stages. Single-cell transcriptome/spatial transcriptomics analysis helps to investigate complex cellular systems at the single-cell level, reflecting the biological complexity of different cell types. In this review, we highlight recent progress that facilitates the development of the female germline to explore the roles of crucial gene regulatory networks, epigenetic pathways, cell-cycle regulators, and phytohormones in this process. This review discusses three key phases in female germline development and provides the possibility of distinct pathways restricting germline development in the future.

Long noncoding RNA DREAMer bridges the DREAM complex and E2f1 to regulate endoreplication in Drosophila

Rb/E2f and DREAM complexes play vital roles in regulating cell cycle progression. To date, how they coordinate their functions to regulate cell cycle-dependent gene expression is not clear. Here, we identified a long noncoding RNA (lncRNA), which we named DREAMer, that bridges the interaction between E2f1 and the dREAM complex to regulate endoreplication specifically in Drosophila salivary gland. We show that E2f1 directly stimulates DREAMer expression, whereas DREAMer mediates the repression of e2f1 transcription by modulating the recruitment of the dREAM complex to the e2f1 promoter via a direct interaction with the dREAM component E2f2. The depletion of DREAMer impairs dREAM binding, leading to derepression of e2f1 transcription, which ultimately increases E2f1 activity and promotes the endoreplication. Furthermore, the transcriptomic analysis revealed profound changes in cell cycle-related gene expression in DREAMerKO salivary glands. Together, our findings reveal an lncRNA-mediated link between the dREAM complex and E2f1, which regulates endoreplication during development.

Quiescence-Origin Senescence: A New Paradigm in Cellular Aging

Cellular senescence, traditionally viewed as a consequence of proliferating and growing cells overwhelmed by extensive stresses and damage, has long been recognized as a critical cellular aging mechanism. Recent research, however, has revealed a novel pathway termed "quiescence-origin senescence", where cells directly transition into senescence from the quiescent state, bypassing cell proliferation and growth. This opinion paper presents a framework conceptualizing a continuum between quiescence and senescence with quiescence deepening as a precursor to senescence entry. We explore the triggers and controllers of this process and discuss its biological implications. Given that the majority of cells in the human body are dormant rather than proliferative, understanding quiescence-origin senescence has significant implications for tissue homeostasis, aging, cancer, and various disease processes. The new paradigm in exploring this previously overlooked senescent cell population may reshape our intervention strategies for age-related diseases and tissue regeneration.

Single-cell RNA-seq reveals the transcriptional program underlying tumor progression and metastasis in neuroblastoma

Neuroblastoma (NB) is one of the most common childhood malignancies. Sixty percent of patients present with widely disseminated clinical signs at diagnosis and exhibit poor outcomes. However, the molecular mechanisms triggering NB metastasis remain largely uncharacterized. In this study, we generated a transcriptomic atlas of 15 447 NB cells from eight NB samples, including paired samples of primary tumors and bone marrow metastases. We used time-resolved analysis to chart the evolutionary trajectory of NB cells from the primary tumor to the metastases in the same patient and identified a common 'starter' subpopulation that initiates tumor development and metastasis. The 'starter' population exhibited high expression levels of multiple cell cycle-related genes, indicating the important role of cell cycle upregulation in NB tumor progression. In addition, our evolutionary trajectory analysis demonstrated the involvement of partial epithelial-to-mesenchymal transition (p-EMT) along the metastatic route from the primary site to the bone marrow. Our study provides insights into the program driving NB metastasis and presents a signature of metastasis-initiating cells as an independent prognostic indicator and potential therapeutic target to inhibit the initiation of NB metastasis.

Genome-wide characterization and expression profiling of E2F/DP gene family members in response to abiotic stress in tomato (Solanum lycopersicum L.)

BACKGROUND

E2F/DP (Eukaryotic 2 transcription factor/dimerization partner) family proteins play an essential function in the cell cycle development of higher organisms. E2F/DP family genes have been reported only in a few plant species. However, comprehensive genome-wide characterization analysis of the E2F/DP gene family of Solanum lycopersicum has not been reported so far.

RESULTS

This study identified eight nonredundant SlE2F/DP genes that were classified into seven groups in the phylogenetic analysis. All eight genes had a single E2F-TDP domain and few genes had additional domains. Two segmental duplication gene pairs were observed within tomato, in addition to cis-regulatory elements, miRNA target sites and phosphorylation sites which play an important role in plant development and stress response in tomato. To explore the three-dimensional (3D) models and gene ontology (GO) annotations of SlE2F/DP proteins, we pointed to their putative transporter activity and their interaction with several putative ligands. The localization of SlE2F/DP-GFP fused proteins in the nucleus and endoplasmic reticulum suggested that they may act in other biological functions. Expression studies revealed the differential expression pattern of most of the SlE2F/DP genes in various organs. Moreover, the expression of E2F/DP genes against abiotic stress, particularly SlE2F/DP2 and/or SlE2F/DP7, was upregulated in response to heat, salt, cold and ABA treatment. Furthermore, the co-expression analysis of SlE2F/DP genes with multiple metabolic pathways was co-expressed with defence genes, transcription factors and so on, suggested their crucial role in various biological processes.

CONCLUSIONS

Overall, our findings provide a way to understand the structure and function of SlE2F/DP genes; it might be helpful to improve fruit development and tolerance against abiotic stress through marker-assisted selection or transgenic approaches.

The association of E2F1 and E2F2 single nucleotide polymorphisms with laryngeal squamous cell carcinoma pathomorphological features

BACKGROUND

Laryngeal squamous cell carcinoma (LSCC) is one of the most common types of cancer in the upper respiratory tract. It is well-known that it has a high mortality rate and poor prognosis in advanced stages. There are well-known risk factors for LSCC, though new specific and prognostic blood-based markers for LSCC development and prognosis are essential. The current study aimed to evaluate the impact of four different single nucleotide polymorphisms (SNPs), E2F1 (rs3213183 and rs3213180) and E2F2 (rs2075993 and rs3820028), on LSCC development, morphological features, and patient 5-year survival rate.

METHODS

A total of 200 LSCC patients and 200 controls were included in this study; both groups were matched by age and sex. In the present study, we analyzed four single nucleotide polymorphisms (SNPs) in the genes E2F1 (rs3213183 and rs3213180) and E2F2 (rs2075993 and rs3820028) and evaluated their associations with the risk of LSCC development, its clinical and morphological manifestation, and patients 5-year survival rate. Genotyping was carried out using RT-PCR.

RESULTS

None of the analyzed SNPs showed a direct association with LSCC development. E2F2 rs2075993 G allele carriers (OR = 4.589, 95% CI 1.050-20.051, p = 0.043) and rs3820028 A allele carriers (OR = 4.750, 95% CI 1.088-20.736, p = 0.038) had a statistically significantly higher risk for poor differentiated or undifferentiated LSCC than non-carriers. E2F1 rs3213180 GC heterozygotes were found to have a 3.7-fold increased risk for lymph node involvement (OR = 3.710, 95% CI 1.452-9.479, p = 0.006). There was no statistically significant association between investigated SNPs and patient 5-year survival rate.

CONCLUSIONS

The present study indicates that E2F2 rs2075993 and rs3820028 impact LSCC differentiation, whereas E2F1 rs3213180 - on lymph node involvement.

Upregulation of CCNB2 and a novel lncRNAs-related risk model predict prognosis in clear cell renal cell carcinoma

BACKGROUND

Clear cell renal cell carcinoma (ccRCC) is the main type of renal cell carcinoma. Cyclin B2 (CCNB2) is a subtype of B-type cyclin that is associated with the prognosis of several cancers. This study aimed to identify the relationship between CCNB2 and progression of ccRCC and construct a novel lncRNAs-related model to predict prognosis of ccRCC patients.

METHODS

The data were obtained from public databases. We identified CCNB2 in ccRCC using Kaplan-Meier survival analysis, univariate and multivariate Cox regression, and Gene Ontology analysis. External validation was then performed. The risk model was constructed based on prognostic lncRNAs by the LASSO algorithm and multivariate Cox regression. Receiver operating characteristics (ROC) curves were used to evaluate the model. Consensus clustering analysis was performed to re-stratify the patients. Finally, we analyzed the tumor-immune microenvironment and performed screening of potential drugs.

RESULTS

CCNB2 associated with late clinicopathological parameters and poor prognosis in ccRCC and was an independent predictor for disease-free survival. In addition, CCNB2 shared the same expression pattern with known suppressive immune checkpoints. A risk model dependent on the expression of three prognostic CCNB2-related lncRNAs (SNHG17, VPS9D1-AS1, and ZMIZ1-AS1) was constructed. The risk signature was an independent predictor of ccRCC. The area under the ROC (AUC) curve for overall survival at 1-, 3-, 5-, and 8-year was 0.704, 0.702, 0.741, and 0.763. The high-risk group and cluster 2 had stronger immunogenicity and were more sensitive to immunotherapy.

CONCLUSION

CCNB2 could be an important biomarker for predicting prognosis in ccRCC patients. Furthermore, we developed a novel lncRNAs-related risk model and identified two CCNB2-related molecular clusters. The risk model performed well in predicting overall survival and immunological microenvironment of ccRCC.

E2F7 serves as a potential prognostic biomarker for lung adenocarcinoma

E2F transcription factors (E2Fs) are a family of transcription factors critical regulators of the cell cycle, apoptosis, and differentiation, thus influencing tumorigenesis. However, the specific roles of E2Fs in lung adenocarcinoma (LUAD) remain unclear. Data from The Cancer Genome Atlas (TCGA) were used. R version. 4.0.3 and multiple databases (TIMER, cBioportal, gene expression profile interaction analysis [GEPIA], LinkedOmics, and CancerSEA) were utilized to investigate mRNA expression, mutational analysis, prognosis, clinical correlations, co-expressed gene, pathway and network, and single-cell analyses. Immunohistochemistry (IHC) confirmed that E2F transcription factor 7 (E2F7) correlated with LUAD. Among the E2Fs, E2F7 was identified by constructing a prognostic model most significantly associated with overall survival (OS) in LUAD patients. The univariate and multivariate Cox regression analyses showed that E2F7, p-T stage, and p-TNM stage were closely related to OS and progression-free survival (PFS) (P < .05) in LUAD. E2F 7/8 were also identified as significantly associated with tumor stage in the GEPIA database. Compared with paracancerous tissues, E2F7 was up-regulated in LUAD by IHC, and E2F7 might be positively correlated with larger tumors and higher TNM stages. E2F7 may primarily regulate DNA repair, damage, and cell cycle processes and thus affect LUAD tumorigenesis, invasion, and metastasis by LinkedOmics and CancerSEA. E2F7 serves as a potential prognostic biomarker for LUAD.

B4GALNT1 promotes carcinogenesis by regulating epithelial-mesenchymal transition in hepatocellular carcinoma based on pan-cancer analysis

BACKGROUND

β-1,4-N-Acetyl galactosaminyltransferase1 (B4GALNT1) has been reported to play important roles in tumor progression and metastasis. Herein, we investigate the oncogenic roles of B4GALNT1 for hepatocellular carcinoma (HCC) based on immune microenvironment.

METHODS

The Cancer Genome Atlas database and Genotype-Tissue Expression, cBioportal, ImmuCellAI, TIMER2 and other databases were searched to analyze the expression and clinical applications of B4GALNT1 in liver cancer patients. Kaplan-Meier survival analysis, Cox regression analysis, Kyoto Encyclopedia of Genes and Genomes and Gene Ontology enrichment analysis were utilized. Moreover, western blot assay, immune histochemistry staining, Cell Counting Kit-8 (CCK-8) assay, invasion and migration assay were performed to evaluate the function of B4GALNT1 in HCC.

RESULTS

B4GALNT1 is overexpressed in 14 tumors, and the mRNA expression levels of B4GALNT1 were remarkably elevated in most tumor types, including HCC. In addition, B4GALNT1 expression was an independent prognostic factor, and a low expression of B4GALNT1 showed a better overall survival and disease-specific recurrence-free survival in patients with HCC. Gene set variation analysis indicated that B4GALNT1 presented a positive correlation with the epithelial-mesenchymal transition (EMT) pathway in HCC. B4GALNT1 expression was closely associated with immune activation genes in the HCC microenvironment. Moreover, B4GALNT1 expression was higher in HCC tissue than that in surrounding tissues. B4GALNT1 promoted the expression of apoptosis-related or EMT-related proteins and then decreased the expression of Bcl-2 and Bcl-xl in HCC cells, suggesting that B4GALNT1 knockdown significantly inhibited the proliferation and invasion of HCC cells.

CONCLUSIONS

B4GALNT1 may promote HCC development through regulating the EMT pathway, which suggests that B4GALNT1 may serve as a promising predictive biomarker and a potential therapeutic target for HCC.

Histone methyltransferase activity affects metabolism in human cells independently of transcriptional regulation

The N-terminal tails of eukaryotic histones are frequently posttranslationally modified. The role of these modifications in transcriptional regulation is well-documented. However, the extent to which the enzymatic processes of histone posttranslational modification might affect metabolic regulation is less clear. Here, we investigated how histone methylation might affect metabolism using metabolomics, proteomics, and RNA-seq data from cancer cell lines, primary tumour samples and healthy tissue samples. In cancer, the expression of histone methyltransferases (HMTs) was inversely correlated to the activity of NNMT, an enzyme previously characterised as a methyl sink that disposes of excess methyl groups carried by the universal methyl donor S-adenosyl methionine (SAM or AdoMet). In healthy tissues, histone methylation was inversely correlated to the levels of an alternative methyl sink, PEMT. These associations affected the levels of multiple histone marks on chromatin genome-wide but had no detectable impact on transcriptional regulation. We show that HMTs with a variety of different associations to transcription are co-regulated by the Retinoblastoma (Rb) tumour suppressor in human cells. Rb-mutant cancers show increased total HMT activity and down-regulation of NNMT. Together, our results suggest that the total activity of HMTs affects SAM metabolism, independent of transcriptional regulation.

Integrative bioinformatics analysis of WDHD1: a potential biomarker for pan-cancer prognosis, diagnosis, and immunotherapy

BACKGROUND

Although WD repeat and high-mobility group box DNA binding protein 1 (WDHD1) played an essential role in DNA replication, chromosome stability, and DNA damage repair, the panoramic picture of WDHD1 in human tumors remains unclear. Hence, this study aims to comprehensively characterize WDHD1 across 33 human cancers.

METHODS

Based on publicly available databases such as TCGA, GTEx, and HPA, we used a bioinformatics approach to systematically explore the genomic features and biological functions of WDHD1 in pan-cancer.

RESULTS

WDHD1 mRNA levels were significantly increased in more than 20 types of tumor tissues. Elevated WDHD1 expression was associated with significantly shorter overall survival (OS) in 10 tumors. Furthermore, in uterine corpus endometrial carcinoma (UCEC) and liver hepatocellular carcinoma (LIHC), WDHD1 expression was significantly associated with higher histological grades and pathological stages. In addition, WDHD1 had a high diagnostic value among 16 tumors (area under the ROC curve [AUC] > 0.9). Functional enrichment analyses suggested that WDHD1 probably participated in many oncogenic pathways such as E2F and MYC targets (false discovery rate [FDR] < 0.05), and it was involved in the processes of DNA replication and DNA damage repair (p.adjust < 0.05). WDHD1 expression also correlated with the half-maximal inhibitory concentrations (IC50) of rapamycin (4 out of 10 cancers) and paclitaxel (10 out of 10 cancers). Overall, WDHD1 was negatively associated with immune cell infiltration and might promote tumor immune escape. Our analysis of genomic alterations suggested that WDHD1 was altered in 1.5% of pan-cancer cohorts and the "mutation" was the predominant type of alteration. Finally, through correlation analysis, we found that WDHD1 might be closely associated with tumor heterogeneity, tumor stemness, mismatch repair (MMR), and RNA methylation modification, which were all processes associated with the tumor progression.

CONCLUSIONS

Our pan-cancer analysis of WDHD1 provides valuable insights into the genomic characterization and biological functions of WDHD1 in human cancers and offers some theoretical support for the future use of WDHD1-targeted therapies, immunotherapies, and chemotherapeutic combinations for the management of tumors.

PP2A and its adapter protein IER5 induce the DNA-binding ability and target gene expression of E2F1 via dephosphorylation at serine 375

The transcription factor E2F1 participates in cell cycle control through transcriptional activation of genes that promote S-phase entry. E2F1 is also linked to the expression of proapoptotic genes, and the loss of E2F1 activity facilitates tumor progression by reducing cellular apoptosis. Phosphorylation controlled by protein kinases and phosphatases is the major posttranslational modification and regulates the cellular levels and transactivator function of E2F1. Here, we characterize the regulatory roles of serine-375 (S375), one of the major phosphorylation sites of E2F1. Cyclin-dependent kinases such as CDK8 phosphorylate at S375 of E2F1, which is dephosphorylated by protein phosphatase 2A (PP2A) containing the B55 regulatory subunit. The PP2A adapter protein IER5 binds to both PP2A/B55 and E2F1 and assists dephosphorylation at S375 by PP2A. S375-dephosphorylated E2F1 exhibits higher DNA-binding affinity than the phosphorylated form. Although the promoter regions of proapoptotic genes are less occupied by E2F1 in cells, an increase in S375-dephosphorylated E2F1 induces preferential binding of E2F1 to the proapoptotic gene promoters and their expression. Our data identify PP2A/B55-IER5 as a critical regulator of E2F1 and suggest that the phosphorylation state of E2F1 is an important determinant for the expression of proapoptotic genes.

HBx promotes hepatocellular carcinoma progression by repressing the transcription level of miR-187-5p

HBV-associated hepatitis B virus x protein (HBx) plays multiple roles in the development of hepatocellular carcinoma. In our prior study, we discovered that miR-187-5p expression was inhibited by HBx. To investigate the underlying molecular mechanism of HBx-mediated miR-187-5p downregulation in hepatocellular carcinoma cells, effects of HBx and miR-187-5p on hepatoma carcinoma cell were observed, as well as their interactions. Through in vitro and in vivo experiments, we demonstrated that overexpression of miR-187-5p inhibited proliferation, migration, and invasion. Simultaneously, we observed a dysregulation in the expression of miR-187-5p in liver cancer cell lines, which may be attributed to transcriptional inhibition through the E2F1/FoxP3 axis. Additionally, we noted that HBx protein is capable of enhancing the expression of E2F1, a transcription factor that promotes the expression of FoxP3. In conclusion, our results suggest that the inhibitory effect of HBx on miR-187-5p is mediated through the E2F1/FoxP3 axis. As shown in this work, HBx promotes hepatoma carcinoma cell proliferation, migration, and invasion through the E2F1/FoxP3/miR-187 axis. It provides a theoretical basis for finding therapeutic targets that will help clinic treatment for HCC.

CDK4-E2F3 signals enhance oxidative skeletal muscle fiber numbers and function to affect myogenesis and metabolism

Understanding how skeletal muscle fiber proportions are regulated is vital to understanding muscle function. Oxidative and glycolytic skeletal muscle fibers differ in their contractile ability, mitochondrial activity, and metabolic properties. Fiber-type proportions vary in normal physiology and disease states, although the underlying mechanisms are unclear. In human skeletal muscle, we observed that markers of oxidative fibers and mitochondria correlated positively with expression levels of PPARGC1A and CDK4 and negatively with expression levels of CDKN2A, a locus significantly associated with type 2 diabetes. Mice expressing a constitutively active Cdk4 that cannot bind its inhibitor p16INK4a, a product of the CDKN2A locus, were protected from obesity and diabetes. Their muscles exhibited increased oxidative fibers, improved mitochondrial properties, and enhanced glucose uptake. In contrast, loss of Cdk4 or skeletal muscle-specific deletion of Cdk4's target, E2F3, depleted oxidative myofibers, deteriorated mitochondrial function, and reduced exercise capacity, while increasing diabetes susceptibility. E2F3 activated the mitochondrial sensor PPARGC1A in a Cdk4-dependent manner. CDK4, E2F3, and PPARGC1A levels correlated positively with exercise and fitness and negatively with adiposity, insulin resistance, and lipid accumulation in human and rodent muscle. All together, these findings provide mechanistic insight into regulation of skeletal muscle fiber-specification that is of relevance to metabolic and muscular diseases.

Regeneration of Pancreatic β-Cells for Diabetes Therapeutics by Natural DYRK1A Inhibitors

The pathogenesis of diabetes mellitus is characterized by insulin resistance and islet β-cell dysfunction. Up to now, the focus of diabetes treatment has been to control blood glucose to prevent diabetic complications. There is an urgent need to develop a therapeutic approach to restore the mass and function of β-cells. Although exogenous islet cell transplantation has been used to help patients control blood glucose, it is costly and has very narrow application scenario. So far, small molecules have been reported to stimulate β-cell proliferation and expand β-cell mass, increasing insulin secretion. Dual-specificity tyrosine-regulated kinase 1A (DYRK1A) inhibitors can induce human β-cell proliferation in vitro and in vivo, and show great potential in the field of diabetes therapeutics. From this perspective, we elaborated on the mechanism by which DYRK1A inhibitors regulate the proliferation of pancreatic β-cells, and summarized several effective natural DYRK1A inhibitors, hoping to provide clues for subsequent structural optimization and drug development in the future.

MNX1-AS1 Promotes Phase Separation of IGF2BP1 to Drive c-Myc-Mediated Cell-Cycle Progression and Proliferation in Lung Cancer

c-Myc and E2F1 play critical roles in many human cancers. As long noncoding RNAs (lncRNA) are known to regulate various tumorigenic processes, elucidation of mechanisms of cross-talk between lncRNAs and c-Myc/E2F1-related signaling pathways could provide important insights into cancer biology. In this study, we used integrated bioinformatic analyses and found that the lncRNA MNX1-AS1 is upregulated in non-small cell lung cancer (NSCLC) via copy-number gain and c-Myc-mediated transcriptional activation. High levels of MNX1-AS1 were associated with poor clinical outcomes in patients with lung cancer. MNX1-AS1 promoted cell proliferation and colony formation in vitro and tumor growth in vivo. MNX1-AS1 bound and drove phase separation of IGF2BP1, which increased the interaction of IGF2BP1 with the 3'-UTR (untranslated region) of c-Myc and E2F1 mRNA to promote their stability. The c-Myc/MNX1-AS1/IGF2BP1 positive feedback loop accelerated cell-cycle progression and promoted continuous proliferation of lung cancer cells. In a lung cancer patient-derived xenograft model, inhibition of MNX1-AS1 suppressed cancer cell proliferation and tumor growth. These findings offer new insights into the regulation and function of c-Myc and E2F1 signaling in NSCLC tumorigenesis and suggest that the MNX1-AS1/IGF2BP1 axis may serve as a potential biomarker and therapeutic target in NSCLC.

SIGNIFICANCE

MNX1-AS1 drives phase separation of IGF2BP1 to increase c-Myc and E2F1 signaling and to activate cell-cycle progression to promote proliferation in NSCLC.

Network analysis reveals that the tumor suppressor lncRNA GAS5 acts as a double-edged sword in response to DNA damage in gastric cancer

The lncRNA GAS5 acts as a tumor suppressor and is downregulated in gastric cancer (GC). In contrast, E2F1, an important transcription factor and tumor promoter, directly inhibits miR-34c expression in GC cell lines. Furthermore, in the corresponding GC cell lines, lncRNA GAS5 directly targets E2F1. However, lncRNA GAS5 and miR-34c remain to be studied in conjunction with GC. Here, we present a dynamic Boolean network to classify gene regulation between these two non-coding RNAs (ncRNAs) in GC. This is the first study to show that lncRNA GAS5 can positively regulate miR-34c in GC through a previously unknown molecular pathway coupling lncRNA/miRNA. We compared our network to several in-vivo/in-vitro experiments and obtained an excellent agreement. We revealed that lncRNA GAS5 regulates miR-34c by targeting E2F1. Additionally, we found that lncRNA GAS5, independently of p53, inhibits GC proliferation through the ATM/p38 MAPK signaling pathway. Accordingly, our results support that E2F1 is an engaging target of drug development in tumor growth and aggressive proliferation of GC, and favorable results can be achieved through tumor suppressor lncRNA GAS5/miR-34c axis in GC. Thus, our findings unlock a new avenue for GC treatment in response to DNA damage by these ncRNAs.

Identification and Expression Analysis of the Solanum tuberosum StATG8 Family Associated with the WRKY Transcription Factor

Autophagy is an evolutionarily well-conserved cellular catabolic pathway in eukaryotic cells and plays an important role in cellular processes. Autophagy is regulated by autophagy-associated (ATG) proteins. Among these ATG proteins, the ubiquitin-like protein ATG8/LC3 is essential for autophagosome formation and function. In this study, the potato StATG8 family showed clade I and clade II with significantly different sequences. Expression of the StATG8 family was also increased in senescence. Interestingly, the expression of the StATG8 and other core StATG genes decreased in potato tubers as the tubers matured. The StATG8 family also responded to a variety of stresses such as heat, wounding, salicylic acid, and salt stress. We also found that some Arabidopsis WRKY transcription factors interacted with the StATG8 protein in planta. Based on group II-a WRKY, StATG8-WRKY interaction is independent of the ATG8 interacting motif (AIM) or LC3 interacting region (LIR) motif. This study showed that the StATG8 family had diverse functions in tuber maturation and multiple stress responses in potatoes. Additionally, StATG8 may have an unrelated autophagy function in the nucleus with the WRKY transcription factor.

Comprehensive analysis identifies as a critical prognostic prediction gene in breast cancer

BACKGROUND

Aurora kinases (AURKs) family plays a vital role not only in cell division but also in tumorigenesis. However, there are still rare systematic analyses of the diverse expression patterns and prognostic value of the AURKs family in breast cancer (BC). Systematic bioinformatics analysis was conducted to explore the biological role, prognostic value, and immunologic function of AURKs family in BC.

METHODS

The expression, prognostic value, and clinical functions of AURKs family in BC were evaluated with several bioinformatics web portals: ONCOMINE Gene Expression Profiling Interactive Analysis, Kaplan-Meier plotter, cBioPortal, Metascape, GeneMANIA, and LinkedOmics; and the result was verified using human tissues.

RESULTS

The expression of AURKA and AURKB were upregulated in BC in subgroup analyses based on tumor stage (all P   <  0.05). BC patients with high AURKA and AURKB expression had a worse overall survival, relapse-free survival, and distant metastasis-free survival (all P   <  0.05). Verification experiment revealed that AURKA and AURKB were upregulated in BC ( P  < 0.05). AURKA and AURKB were specifically associated with several tumor-associated kinases (polo-like kinase 1 and cyclin-dependent kinase 1), miRNAs (miR-507 and miR-381), and E2F transcription factor 1. Moreover, AURKA and AURKB were correlated with immune cell infiltration. Functional enrichment analysis revealed that AURKA and AURKB were involved in the cell cycle signaling pathway, platinum drug resistance signaling pathway, ErbB signaling pathway, Hippo signaling pathway, and nucleotide-binding and oligomerization domain-like receptor signaling pathway.

CONCLUSIONS

Aurora kinases AURKA and AURKB could be employed as novel prognostic biomarkers or promising therapeutic targets for BC.

Genetic determinants of lung cancer: Understanding the oncogenic potential of somatic missense mutations

BACKGROUND

Treatment of lung cancer is getting more personalized nowadays and medical practitioners are moving away from conventional histology-driven empirical treatments, platinum-based chemotherapy, and other invasive surgical resections and have started adopting alternate therapies in which therapeutic targets are patient's molecular oncogenic drivers.

AIM

The aim of the current study is to extract meaningful information from the online somatic mutation data (retrieved from cBioPortal) of 16 most significantly mutated oncogenes in non-small-cell lung cancer (NSCLC), namely EGFR, NRAS, KRAS, HER2 (ERBB2), RET, MET, ROS1, FGFR1, BRAF, AKT1, MEK1 (MAP2K1), PIK3CA, PTEN, DDR2, LKB1 (STK11) and ALK, for improving our understanding of the pathobiology of the lung cancer that can aid decision-making on critical clinical and therapeutic considerations.

METHODS

Using an integrated approach comprising 4 steps, the oncogenic potential of 661 missense non-synonymous single nucleotide polymorphisms (nsSNPs) in 16 genes was ascertained using 2059 NSCLC (1575 lung adenocarcinomas, 484 lung squamous cell carcinomas) patients' online mutation data. The steps used comprise sequence/structure homology-based prediction, scoring of conservation of mutated residues and positions, prediction of resulting molecular and functional consequences using machine-learning and structure-guided approach.

RESULTS

Out of a total of 661 nsSNPs analyzed, a set of 29 nsSNPs has been identified as conserved high confidence mutations in 10 of 16 genes relevant to the under study. Out of 29 conserved high confidence nsSNPs, 4 nsSNPs (EGFR N1094Y, BRAF M620I, DDR2 R307L, ALK P1350T) have been found to be putative novel rare genetic markers for NSCLC.

CONCLUSIONS

The current study, the first of its kind, has provided a list of deleterious non-synonymous somatic mutations in a selected pool of oncogenes that can be considered as a promising target for future drug design and therapy for patients with lung adenocarcinomas and squamous cell carcinomas.

The miR-532-E2F1 feedback loop contributes to gastric cancer progression

Gastric cancer (GC) ranks fourth in incidence and mortality worldwide, ascertaining the pathogenesis of GC is crucial for its treatment. E2F1, which regulates the transcription of genes encoding proteins involved in DNA repair, DNA replication, mitosis and survival of cancer patients, functions as a key regulator in GC progression. However, the underneath mechanism of these processes is not fully elucidated. Here, TCGA database analysis, microarray immunohistochemical technique and western blot showed that E2F1 was highly upregulated in clinical GC tissues and correlated with tumor malignancy. In vitro and in vivo assays confirmed the oncogenic function of E2F1. MiR-532 was decreased and negatively correlated with E2F1 in GC tissues. MiR-532 directly targeted and inhibited E2F1 expression, leading to the decrease of ASK1 and elevation of TXNIP, and affected proliferation, cell cycle, apoptosis and DNA damage in vitro and tumor growth in vivo. Moreover, E2F1 serves as a transcriptional repressor to suppress miR-532 expression and a double-negative feedback loop was formed between them. This study demonstrates the significant roles of the E2F1-miR-532 double-negative feedback loop in GC progression and may represent a potential target for GC therapy.

Construction of Bone Metastasis-Specific Regulation Network Based on Prognostic Stemness-Related Signatures in Prostate Cancer

Background

We planned to uncover the cancer stemness-related genes (SRGs) in prostate cancer (PCa) and its underlying mechanism in PCa metastasis.

Methods

We acquired the RNA-seq data of 406 patients with PCa from the TCGA database. Based on the mRNA stemness index (mRNAsi) calculated by one-class logistic regression (OCLR) algorithm, SRGs in PCa were extracted by WGCNA. Univariate and multivariate regression analyses were applied to uncover OS-associated SRGs. Gene Set Variation Analysis (GSVA), Gene Set Enrichment Analysis (GSEA), and Pearson's correlation analysis were performed to discover the possible mechanism of PCa metastasis. The significantly correlated transcription factors of OS-associated SRGs were also identified by Pearson's correlation analysis. ChIP-seq was applied to validate the binding relationship of TFs and OS-associated SRGs and spatial transcriptome and single-cell sequencing were performed to uncover the location of key biomarkers expression. Lastly, we explored the specific inhibitors for SRGs using CMap algorithm.

Results

We identified 538 differentially expressed genes (DEGs) between non-metastatic and metastatic PCa. Furthermore, OS-associated SRGs were identified. The Pearson correlation analysis revealed that FOXM1 was significantly correlated with NEIL3 (correlation efficient =0.89, p < 0.001) and identified hallmark_E2F_targets as the potential pathway mechanism of NEIL3 promoting PCa metastasis (correlation efficient =0.58, p < 0.001). Single-cell sequencing results indicated that FOXM1 regulating NEIL3 may get involved in the antiandrogen resistance of PCa. Rottlerin was discovered to be a potential target drug for PCa.

Conclusion

We constructed a regulatory network based on SRGs associated with PCa metastasis and explored possible mechanism.

Targeted Re-Sequencing of the 2p21 Locus Identifies Non-Syndromic Cleft Lip Only Novel Susceptibility Gene ZFP36L2

rs7590268 present on the 2p21 locus was identified to be associated with non-syndromic cleft lip with or without cleft palate (NSCL/P) in several populations, including the Chinese Han population, indicating that 2p21 was a susceptibility locus for NSCL/P. However, previous studies have only identified common single-nucleotide polymorphism (SNP) within the THADA gene, neglecting the rare variants and other genes in 2p21; thus, this study was designed to investigate additional variants and novel susceptibility genes in 2p21. A total of 159 NSCL/P patients and 542 controls were recruited in the discovery phase, whereas 1830 NSCL/P patients and 2,436 controls were recruited in the replication phase. After targeted region sequencing, we performed association and burden analyses for the common and rare variants, respectively. Furthermore, RNA-seq, proliferation assay and cell cycle analysis were performed to clarify the possible function of the candidate gene ZFP36L2. Association analysis showed that four SNPs were specifically associated with non-syndromic cleft lip only (NSCLO) and two SNPs were associated with both NSCLO and NSCL/P. Burden analysis indicated that ZFP36L2 was associated with NSCLO (p = .0489, OR = 2.41, 95% CI: 0.98–5.90). Moreover, SNPs in the ZFP36L2 targeted gene JUP were also associated with NSCLO. ZFP36L2 also inhibited cell proliferation and induced G2 phase arrest in the GMSM-K cell line. Therefore, we proposed that ZFP36L2 is a novel susceptibility gene of NSCLO in the 2p21 locus, which could lead to NSCLO by modulating cell proliferation and cycle.

Systematic Elucidation of the Aneuploidy Landscape and Identification of Aneuploidy Driver Genes in Prostate Cancer

Aneuploidy is widely identified as a remarkable feature of malignancy genomes. Increasing evidences suggested aneuploidy was involved in the progression and metastasis of prostate cancer (PCa). Nevertheless, no comprehensive analysis was conducted in PCa about the effects of aneuploidy on different omics and, especially, about the driver genes of aneuploidy. Here, we validated the association of aneuploidy with the progression and prognosis of PCa and performed a systematic analysis in mutation profile, methylation profile, and gene expression profile, which detailed the molecular process aneuploidy implicated. By multi-omics analysis, we managed to identify 11 potential aneuploidy driver genes (GSTM2, HAAO, C2orf88, CYP27A1, FAXDC2, HFE, C8orf88, GSTP1, EFS, HIF3A, and WFDC2), all of which were related to the development and metastasis of PCa. Meanwhile, we also found aneuploidy and its driver genes were correlated with the immune microenvironment of PCa. Our findings could shed light on the tumorigenesis of PCa and provide a better understanding of the development and metastasis of PCa; additionally, the driver genes could be promising and actionable therapeutic targets pointing to aneuploidy.

Exploring the Prognostic Value, Immune Implication and Biological Function of H2AFY Gene in Hepatocellular Carcinoma

Background

Hepatocellular carcinoma (HCC) is an extremely malignant cancer with poor survival. H2AFY gene encodes for a variant of H2A histone, and it has been found to be dysregulated in various tumors. However, the clinical value, biological functions and correlations with immune infiltration of H2AFY in HCC remain unclear.

Methods

We analyzed the expression and clinical significance of H2AFY in HCC using multiple databases, including Oncomine, HCCDB, TCGA, ICGC, and so on. The genetic alterations of H2AFY were analyzed by cBioPortal and COSMIC databases. Co-expression networks of H2AFY and its regulators were investigated by LinkedOmics. The correlations between H2AFY and tumor immune infiltration were explored using TIMER, TISIDB databases, and CIBERSORT method. Finally, H2AFY was knocked down with shRNA lentiviruses in HCC cell lines for functional assays in vitro.

Results

H2AFY expression was upregulated in the HCC tissues and cells. Kaplan-Meier and Cox regression analyses revealed that high H2AFY expression was an independent prognostic factor for poor survival in HCC patients. Functional network analysis indicated that H2AFY and its co-expressed genes regulates cell cycle, mitosis, spliceosome and chromatin assembly through pathways involving many cancer-related kinases and E2F family. Furthermore, we observed significant correlations between H2AFY expression and immune infiltration in HCC. H2AFY knockdown suppressed the cell proliferation and migration, promoted cycle arrest, and apoptosis of HCC cells in vitro.

Conclusion

Our study revealed that H2AFY is a potential biomarker for unfavorable prognosis and correlates with immune infiltration in HCC.

CHPF promotes gastric cancer tumorigenesis through the activation of E2F1

Chondroitin polymerizing factor (CHPF) is an important glycosyltransferase involved in the biosynthesis of chondroitin sulfate. However, the relationship between CHPF and gastric cancer has not been fully investigated. CHPF expression in gastric cancer tissues was detected by immunohistochemistry and correlated with gastric cancer patient prognosis. Cultured gastric cancer cells and human gastric epithelial cell line GES1 were used to investigate the effects of shCHPF and shE2F1 on the development and progression of gastric cancer by MTT, western blotting, flow cytometry analysis of cell apoptosis, colony formation, transwell and gastric cancer xenograft mouse models, in vitro and in vivo. In gastric cancer tissues, CHPF was found to be significantly upregulated, and its expression correlated with tumor infiltration and advanced tumor stage and shorter patient survival in gastric cancer. CHPF may promote gastric cancer development by regulating cell proliferation, colony formation, cell apoptosis and cell migration, while knockdown induced the opposite effects. Moreover, the results from in vivo experiments demonstrated that tumor growth was suppressed by CHPF knockdown. Additionally, E2F1 was identified as a potential downstream target of CHPF in the regulation of gastric cancer, and its knockdown decreased the CHPF-induced promotion of gastric cancer. Mechanistic study revealed that CHPF may regulate E2F1 through affecting UBE2T-mediated E2F1 ubiquitination. This study showed, for the first time, that CHPF is a potential prognostic indicator and tumor promoter in gastric cancer whose function is likely carried out through the regulation of E2F1.

Phytochemicals targeting NF-κB signaling: Potential anti-cancer interventions

Nuclear factor κB (NF-κB) is a ubiquitous regulator of the signalome and is indispensable for various biological cell functions. NF-κB consists of five transcription factors that execute both cytoplasmic and nuclear signaling processes in cells. NF-κB is the only signaling molecule that governs both pro- and anti-apoptotic, and pro- and anti-inflammatory responses. This is due to the canonical and non-canonical components of the NF-κB signaling pathway. Together, these pathways orchestrate cancer-related inflammation, hyperplasia, neoplasia, and metastasis. Non-canonical NF-κB pathways are particularly involved in the chemoresistance of cancer cells. In view of its pivotal role in cancer progression, NF-κB represents a potentially significant therapeutic target for modifying tumor cell behavior. Several phytochemicals are known to modulate NF-κB pathways through the stabilization of its inhibitor, IκB, by inhibiting phosphorylation and ubiquitination thereof. Several natural pharmacophores are known to inhibit the nuclear translocation of NF-κB and associated pro-inflammatory responses and cell survival pathways. In view of this and the high degree of specificity exhibited by various phytochemicals for the NF-κB component, we herein present an in-depth overview of these phytochemicals and discuss their mode of interaction with the NF-κB signaling pathways for controlling the fate of tumor cells for cancer-directed interventions.

•

NF-κB plays a pivotal role in the maintenance of homeostasis and various inflammation-mediated pathologies.

•

NF-κB is involved in cancer development and progression by modulating growth signaling and apoptosis pathways.

•

Phytochemicals modulating NF-κB activity should be exploited to design anticancer drugs with minimal side effects.

•

Use of these phytochemicals in adjunctive chemotherapy may enhance the chemosensitivity of existing chemotherapeutic drugs.

Suppression of HSF1 activity by wildtype p53 creates a driving force for p53 loss-of-heterozygosity

The vast majority of human tumors with p53 mutations undergo loss of the remaining wildtype p53 allele (loss-of-heterozygosity, p53LOH). p53LOH has watershed significance in promoting tumor progression. However, driving forces for p53LOH are poorly understood. Here we identify the repressive WTp53-HSF1 axis as one driver of p53LOH. We find that the WTp53 allele in AOM/DSS chemically-induced colorectal tumors (CRC) of p53R248Q/+ mice retains partial activity and represses heat-shock factor 1 (HSF1), the master regulator of the proteotoxic stress response (HSR) that is ubiquitously activated in cancer. HSR is critical for stabilizing oncogenic proteins including mutp53. WTp53-retaining CRC tumors, tumor-derived organoids and human CRC cells all suppress the tumor-promoting HSF1 program. Mechanistically, retained WTp53 activates CDKN1A/p21, causing cell cycle inhibition and suppression of E2F target MLK3. MLK3 links cell cycle with the MAPK stress pathway to activate the HSR response. In p53R248Q/+ tumors WTp53 activation by constitutive stress represses MLK3, thereby weakening the MAPK-HSF1 response necessary for tumor survival. This creates selection pressure for p53LOH which eliminates the repressive WTp53-MAPK-HSF1 axis and unleashes tumor-promoting HSF1 functions, inducing mutp53 stabilization enabling invasion.

PRMT5 Selective Inhibitor Enhances Therapeutic Efficacy of Cisplatin in Lung Cancer Cells

As a therapeutic approach, epigenetic modifiers have the potential to enhance the efficacy of chemotherapeutic agents. Protein arginine methyltransferase 5 (PRMT5), highly expressed in lung adenocarcinoma, was identified to be involved in tumorigenesis. In the current study, we examined the potential antineoplastic activity of PRMT5 inhibitor, arginine methyltransferase inhibitor 1 (AMI-1), and cisplatin on lung adenocarcinoma. Bioinformatic analyses identified apoptosis, DNA damage, and cell cycle progression as the main PRMT5-associated functional pathways, and survival analysis linked the increased PRMT5 gene expression to worse overall survival in lung adenocarcinoma. Combined AMI-1 and cisplatin treatment significantly reduced cell viability and induced apoptosis. Cell cycle arrest in A549 and DMS 53 cells was evident after AMI-1, and was reinforced after combination treatment. Western blot analysis showed a reduction in demethylation histone 4, a PRMT5- downstream target, after treatment with AMI-1 alone or in combination with cisplatin. While the combination approach tackled lung cancer cell survival, it exhibited cytoprotective abilities on HBEpC (normal epithelial cells). The survival of normal bronchial epithelial cells was not affected by using AMI-1. This study highlights evidence of novel selective antitumor activity of AMI-1 in combination with cisplatin in lung adenocarcinoma cells.

A novel E2F1-regulated lncRNA, LAPAS1, is required for S phase progression and cell proliferation

The transcription factor E2F1 induces both proliferation and apoptosis and is a critical downstream target of the tumor suppressor RB.

Long non-coding RNAs (lncRNAs) are major regulators of many cellular processes, including cell cycle progression and cell proliferation. However, the mode of action as well as the transcriptional regulation of most lncRNAs are only beginning to be understood.

Here, we report that a novel human lncRNA, LAPAS1, is an E2F1- regulated lncRNA that affects S phase progression. Inhibition of LAPAS1 expression increases percentage of S phase cells, and its silencing in synchronized cells delays their progression through S phase. In agreement with its suggested role in cell cycle progression, prolonged inhibition of LAPAS1 attenuates proliferation of human cancer cells.

Our data demonstrate that LAPAS1 predominantly functions in trans to repress expression of Sphingolipid Transporter 2 (SPNS2). Importantly, knockdown of SPNS2 rescues the effect of LAPAS1 silencing on cell cycle and cell proliferation.

Notably, low levels of LAPAS1 are associated with increased survival of kidney cancer patients.

Summarily, we identify LAPAS1 as a novel E2F-regulated lncRNA that has a potential role in human cancer and regulates cell-cycle progression and cell proliferation, at least in part, via regulation of SPNS2.

Expression and prognosis analyses of CASP1 in acute myeloid leukemia

Caspase1 (CASP1) is a gene that encodes multiple proteins related to cell death. Nevertheless, the function of CASP1 in the pathogenesis of AML is still unclear. In the present study, a detailed analysis of cancer versus normal samples was performed to explore the relationship between CASP1 and leukemia. We used sequencing data from multiple cancer gene databases to analyze the gene expression and regulatory network of CASP1 in leukemia. We discovered that mRNA expression levels of CASP1 are increased in leukemia cell lines, especially in acute myelocytic leukemia (AML). Then, we verified the mRNA expression of CASP1 in AML clinical samples and observed significantly higher expression of CASP1 in relapsed AML patients. High CASP1 expression was associated with poor prognosis and CASP1 inhibition could impair the proliferation of AML cells. Related functional network identification suggests that CASP1 regulates apoptosis, immune and inflammatory response via pathways involving LYN, LCK, and the E2F family. These findings suggest that CASP1 probably contributes to the pathogenesis, and identify CASP1 as a factor for predicting the prognosis and as a therapeutic target of AML patients.

Systematic expression analysis of the CELSR family reveals the importance of CELSR3 in human lung adenocarcinoma

Cadherin EGF LAG seven‐pass G‐type receptors (CELSRs) are involved in the progression of various types of cancer. CELSR3, a crucial signalling molecule in the WNT/PCP pathway, is believed to be associated with tumorigenesis and metastasis. However, its role in lung adenocarcinoma (LUAD) remains unclear. In this paper, we analysed the expression of CELSR family members using the Oncomine, GEPIA and UALCAN databases. We used a Kaplan‐Meier plotter to assess the effect of CELSRs on tumour prognosis. Next, gene ontology (GO), KEGG pathway, miRNA target, kinase target and transcription factor‐target enrichment were analysed by GSEA. Simultaneously, we conducted functional assays including cell viability, colony formation and transwell assays, to determine the oncogenic role of CELSR3 in LUAD. Finally, we used the TIMER and TISIDB databases to analyse the correlation between CELSR3 and immune infiltration and the potential chemokine receptor axis causing immune cell expression. High expression of CELSR3 is in LUAD predicts poor prognosis and early progression of the tumour. KEGG and GO enrichment analysis revealed the functional relationship between CELSR3 and cell adhesion, the cell cycle, and DNA replication. Down‐regulation of CELSR3 suppressed cell proliferation to a significant extent, in addition to inhibiting invasion and migration in LUAD cells. Finally, CELSR3 expression was significantly correlated with the infiltration level of CD8+T cells through the CCL17/CCR4 axis in LUAD. These results indicate that CELSR3 can serve as a prognostic biomarker for determining prognosis and immune infiltration in LUAD.

Genome-wide identification and functional characterization of the PheE2F/DP gene family in Moso bamboo

BACKGROUND

E2F/DP proteins have been shown to regulate genes implicated in cell cycle control and DNA repair. However, to date, research into the potential role of the Moso bamboo E2F/DP family has been limited.

RESULTS

Here, we identified 23 E2F/DPs in the Moso bamboo genome, including nine E2F genes, six DP genes, eight DEL genes and one gene with a partial E2F domain. An estimation of the divergence time of the paralogous gene pairs suggested that the E2F/DP family expansion primarily occurred through a whole-genome duplication event. A regulatory element and coexpression network analysis indicated that E2F/DP regulated the expression of cell cycle-related genes. A yeast two-hybrid assay and expression analysis based on transcriptome data and in situ hybridization indicated that the PheE2F-PheDP complex played important roles in winter Moso bamboo shoot growth. The qRT-PCR results showed that the PheE2F/DPs exhibited diverse expression patterns in response to drought and salt treatment and diurnal cycles.

CONCLUSION

Our findings provide novel insights into the Moso bamboo E2F/DP family and partial experimental evidence for further functional verification of the PheE2F/DPs.

Potential role of glucosamine-phosphate N-acetyltransferase 1 in the development of lung adenocarcinoma

Glucosamine-phosphate N-acetyltransferase 1 (GNPNAT1) is a key enzyme associated with glucose metabolism and uridine diphosphate-N-acetylglucosamine biosynthesis. Abnormal GNPNAT1 expression might be associated with carcinogenesis. We analyzed multiple lung adenocarcinoma (LUAD) gene expression databases and verified GNPNAT1 higher expression in LUAD tumor tissues than in normal tissues. Moreover, we analyzed the survival relationship between LUAD patients’ clinical status and GNPNAT1 expression, and found higher GNPNAT1 expression in LUAD patients with unfavorable prognosis. We built GNPNAT1 gene co-expression networks and further annotated the co-expressed genes’ Gene Ontology (GO) terms, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, and various associated regulatory factors. These co-expression genes’ functional networks mainly participate in chromosome segregation, RNA metabolic process, and RNA transport. We analyzed GNPNAT1 genetic alterations and co-occurrence networks, and the functional networks of these genes showed that GNPNAT1 participates in multiple steps of cell cycle transition and in the development of some cancers. We assessed the correlation between GNPNAT1 expression and cancer immune infiltrates and showed that GNPNAT1 expression is correlated with several immune cells, chemokines, and immunomodulators in LUAD. We found that GNPNAT1 correlates with LUAD development and prognosis, laying a foundation for further research, especially in immunotherapy.

Exploration of Prognostic Biomarkers among Replication Factor C Family in the Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is one of the common cancers with a high incidence and mortality. The human replication factor C (RFC) family contains 5 subunits that play an important role in DNA replication and DNA damage repair. RFCs are abnormally expressed in a variety of cancers; some of them are differentially expressed in HCC tissues and related to tumor growth. However, the expression, prognostic value, and effect targets of the whole RFC family in HCC are still unclear. To address these issues, we performed a multidimensional analysis of RFCs in HCC patients by Oncomine, UALCAN, GEPIA, Human protein atlas, Kaplan-Meier plotter, cBioPortal, GeneMANIA, String, and LinkedOmics. mRNA expression of RFCs was significantly increased in HCC tissues. There was a significant correlation between the expression of RFC2/3/4/5 and tumor stage of HCC patients. Besides, high mRNA expression of RFC2/4 was associated with worse overall survival (OS). Moreover, genetic alterations of RFCs were associated with worse OS in HCC patients. We found that genes co-expressed with RFC2/4 were mainly involved in biological processes, such as chromosome segregation, mitotic cell cycle phase transition, and telomere organization and they activated the cell cycle and spliceosome pathways. The gene set is mainly enriched in cancer-related kinases AURKA, ATR, CDK1, PLK1, and CHEK1. E2F family members were the key transcription factors for RFCs. Our results suggest that differentially expressed RFC2 and RFC4 are potential prognostic biomarkers in HCC and may act on E2F transcription factors and some kinase targets to dysregulate the cell cycle pathway. These efforts may provide new research directions for prognostic biomarkers and therapeutic targets in HCC.

p73: From the p53 shadow to a major pharmacological target in anticancer therapy

p73, along with p53 and p63, belongs to the p53 family of transcription factors. Besides the p53-like tumor suppressive activities, p73 has unique roles, namely in neuronal development and differentiation. In addition, the TP73 gene is rarely mutated in tumors. This makes p73 a highly appealing therapeutic target, particularly towards cancers with a null or disrupted p53 pathway. Distinct isoforms are transcribed from the TP73 locus either with (TAp73) and without (ΔNp73) the N-terminal transactivation domain. Conversely to TA tumor suppressors, ΔN proteins exhibit oncogenic properties by inhibiting p53 and TA protein functions. As such, p73 isoforms compose a puzzled and challenging regulatory pathway. This state-of-the-art review affords an update overview on p73 structure, biological functions and pharmacological regulation. Importantly, it addresses the relevance of p73 isoforms in carcinogenesis, highlighting their potential as drug targets in anticancer therapy. A critical discussion of major pharmacological approaches to promote p73 tumor suppressive activities, with relevant survival outcomes for cancer patients, is also provided.

EZH2 reduction is an essential mechanoresponse for the maintenance of super-enhancer polarization against compressive stress in human periodontal ligament stem cells

Despite the ubiquitous mechanical cues at both spatial and temporal dimensions, cell identities and functions are largely immune to the everchanging mechanical stimuli. To understand the molecular basis of this epigenetic stability, we interrogated compressive force-elicited transcriptomic changes in mesenchymal stem cells purified from human periodontal ligament (PDLSCs), and identified H3K27me3 and E2F signatures populated within upregulated and weakly downregulated genes, respectively. Consistently, expressions of several E2F family transcription factors and EZH2, as core methyltransferase for H3K27me3, decreased in response to mechanical stress, which were attributed to force-induced redistribution of RB from nucleoplasm to lamina. Importantly, although epigenomic analysis on H3K27me3 landscape only demonstrated correlating changes at one group of mechanoresponsive genes, we observed a genome-wide destabilization of super-enhancers along with aberrant EZH2 retention. These super-enhancers were tightly bounded by H3K27me3 domain on one side and exhibited attenuating H3K27ac deposition and flattening H3K27ac peaks along with compensated EZH2 expression after force exposure, analogous to increased H3K27ac entropy or decreased H3K27ac polarization. Interference of force-induced EZH2 reduction could drive actin filaments dependent spatial overlap between EZH2 and super-enhancers and functionally compromise the multipotency of PDLSC following mechanical stress. These findings together unveil a specific contribution of EZH2 reduction for the maintenance of super-enhancer stability and cell identity in mechanoresponse.

Long noncoding RNA ERICD interacts with ARID3A via E2F1 and regulates migration and proliferation of osteosarcoma cells

Long noncoding RNA (lncRNA) dysregulation is known to be taking part in majority of cancers, including osteosarcoma. In one of our previous studies, we showed that lncRNA MEG3 is being regulated by microRNA-664a (miR-664a) suppresses the migratory potential of osteosarcoma cells (U-2OS). We now report a novel lncRNA, namely, ERICD, which is linked to the transcription factor AT-rich interaction domain 3A (ARID3A) in U-2OS cells. We show that ARID3A binds to ERICD and indirectly interacts with each other via the E2F transcription factor 1 (E2F1). Furthermore, small interfering RNA (siRNA)-mediated knockdown of ERICD inhibited cell migration, formation of colonies, and proliferation in U-2OS cells. Overexpression of ARID3A inhibited cell migration, colony formation, and proliferation, whereas siRNA-mediated knockdown of ARID3A promoted cell migration, colony formation, and proliferation. Our findings indicate that ARID3A and lncRNA ERICD have plausible tumor suppressive and oncogenic functions, respectively, in osteosarcoma. Our data demonstrate the converse interaction between ARID3A and lncRNA ERICD that target DNA-binding proteins and dysregulation of their expression through E2F1 augments osteosarcoma progression. The cell rescue experiment also indicated E2F1 to be involved in the regulation of ARID3A and ERICD.

E2F transcription factor 8 promotes proliferation and radioresistance in glioblastoma

INTRODUCTION

Glioblastoma (GBM) is the most lethal brain tumor that has a median survival of less than 2 years. Tumor heterogeneity and high therapeutic resistance are hallmarks of GBM. Transcription factors (TFs) play a critical role in tumor progression by regulating the transcriptional events associated with tumor cells transition into more malignant cellular phenotypes. The E2 F transcription factor 8 (E2 F8) is a recently identified TF in the E2 F family. Studies have revealed that E2 F8 is involved in tumorigenesis of multiple cancer types; however, the oncogenic role of E2 F8 in GBM was rarely studied and the prognostic value of E2 F8 has not been explored.

AIMS

In this study, we investigated the expression profile, prognosis correlation and oncogenic role of E2 F8 to explore its potential use as a GBM therapeutic target.

RESULTS

E2 F8 was significantly enriched in GBM compared with normal brain tissues and low-grade glioma. E2 F8 high expression was strongly associated with worse outcome of GBM patients. E2 F8 silencing significantly attenuated the proliferation of tumor cells in vitro and tumorigenicity in vivo, while its overexpression promoted the proliferation of GBM tumor cells. Bioinformatics analysis revealed that E2 F8 was tightly linked to multiple oncogenic processes in GBM, including aggressive cell cycle, DNA repair, STAT3, TGFRβ and WNT pathways. E2 F8 high expression correlated with the expression of a variety of well-known oncogenes in GBM. E2 F8 was identified as a crucial transcriptional regulator of CHEK1 via its directly binding CHEK1 promoter area. Finally, E2 F8 conferred significant radioresistance to GBM tumor cells in vitro and in vivo.

CONCLUSION

E2 F8 is highly expressed in GBM and associated with worse outcome in GBM patients. It promotes tumorigenesis and radioresistance of GBM tumor cells and has oncogenic roles via its involvement in multiple oncogenic processes and pathways such as the regulation of CHEK1 transcriptional activity.

Identification of nsSNPs of transcription factor E2F1 predisposing individuals to lung cancer and head and neck cancer

E2Fs transcription factors family is involved in the G1/S transition and DNA replication and their deregulated expression have been reported in various human cancers. Studies have shown that the genetic variants of E2F1 family members play an important role in head and neck carcinogenesis. In this study, we predicted six highly deleterious nsSNPs (C227F, R252H, V295D, C298Y, R56W, and Y59C) of E2F1 gene through in silico analyses. The latter was based on protein structure, function, and amino acid conservation. Molecular dynamics studies showed a deviation of the structures of the mutant proteins from the global protein parameters. Further, a case-control study that included total 535 samples (305 cancer patients and 230 controls) was conducted to find the association of the predicted SNPs with the susceptibility to lung cancer (LC) and head and neck cancer (HNC). The genotyping was done applying in-house artificial-RFLP method. Statistical analysis showed that the mutant alleles/genotypes of rs3213172 (R252H) were found to increase ∼ 2-5 fold risk of LC and HNC in all the genetic models. These results suggest that the rs3213172C/T polymorphism of the E2F1 gene could be used as an effective biomarker for genetic susceptibility to LC and HNC in our population.

Expression and gene regulation network of RBM8A in hepatocellular carcinoma based on data mining

RNA binding motif protein 8A (RBM8A) is an RNA binding protein in a core component of the exon junction complex. Abnormal RBM8A expression is associated with carcinogenesis. We used sequencing data from the Cancer Genome Atlas database and Gene Expression Omnibus, analyzed RBM8A expression and gene regulation networks in hepatocellular carcinoma (HCC). Expression was analyzed using Oncomine^TM and Gene Expression Profiling Interactive Analysis tools, while RBM8A alterations and related functional networks were identified using cBioPortal. LinkedOmics was used to identify differential gene expression with RBM8A and to analyze Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathways. Gene enrichment analysis examined target networks of kinases, miRNAs and transcription factors. We found that RBM8A is overexpressed and the RBM8A gene often amplified in HCC. Expression of this gene is linked to functional networks involving the ribosome and RNA metabolic signaling pathways. Functional network analysis suggested that RBM8A regulates the spliceosome, ribosome, DNA replication and cell cycle signaling via pathways involving several cancer-related kinases, miRNAs and E2F Transcription Factor 1. Our results demonstrate that data mining efficiently reveals information about RBM8A expression and potential regulatory networks in HCC, laying a foundation for further study of the role of RBM8A in carcinogenesis.

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.

MiR-1285-5p/TMEM194A axis affects cell proliferation in breast cancer

The onset of breast cancer among young patients is a major issue in cancer etiology. Our previous study has shown that poor prognosis in young women with breast cancer is associated with lower expression of the microRNA miR‐1285‐5p. In this study, we showed that the expression of miR‐1285‐5p is lower in tumor tissues than in normal tissues. Accumulating evidence suggests that miR‐1285‐5p plays critical roles in various types of cancers. However, the functional role of miR‐1285‐5p in breast cancer remains to be elucidated. Here, we showed the tumor‐suppressive role of miR‐1285‐5p and detailed its mechanism of action in breast cancer. Overexpression of miR‐1285‐5p significantly inhibited cell proliferation in breast cancer cells regardless of the tumor subtype. Among the target genes of miR‐1285‐5p, we found that transmembrane protein 194A (TMEM194A) was directly regulated by miR‐1285‐5p. Notably, separation of centrosomes from the nuclear envelope was observed upon knockdown of TMEM194A or overexpression of miR‐1285‐5p. In conclusion, our findings show that miR‐1285‐5p is a tumor suppressor via TMEM194A inhibition in breast cancer.

Identification of E2F8 as a Transcriptional Regulator of Gluconeogenesis in Primary Mouse Hepatocytes

The dysregulation of hepatic gluconeogenesis is a major factor in the pathogenesis of type 2 diabetes mellitus (T2DM). Hepatic gluconeogenesis is known to be tightly regulated at the transcription/expression level. The aim of this study was to evaluate the role of the E2F8 transcription factor in glucose metabolism. Here, we found that hepatic expression levels of E2F8 were increased in db/db and high-fat-diet-induced obese mice. Adenovirus-mediated overexpression of E2F8 in primary mouse hepatocytes upregulated expression of gluconeogenic genes, including those for PGC-1α, PEPCK, and G6Pase, subsequently increasing cellular glucose output. We demonstrated that E2F8 overexpression impairs insulin sensitivity in vitro. Furthermore, knockdown of E2F8 expression increased insulin sensitivity in primary hepatocytes. In summary, these findings indicated that E2F8 is involved in gluconeogenesis and insulin resistance and may represent a new therapeutic target in T2DM prevention.

The conflicting role of E2F1 in prostate cancer: A matter of cell context or interpretational flexibility?

The transcription factor E2F1 plays a crucial role in mediating multiple cancer hallmark capabilities that regulate cell cycle, survival, apoptosis, metabolism, and metastasis. Aberrant activation of E2F1 is closely associated with a poor clinical outcome in various human cancers. However, E2F1 has conflictingly been reported to exert tumor suppressive activity, raising a question as to the nature of its substantive role in the control of cell fate. In this review, we summarize deregulated E2F1 activity and its role in prostate cancer. We highlight the recent advances in understanding the molecular mechanism by which E2F1 regulates the development and progression of prostate cancer, providing insight into how cell context or data interpretation shapes the role of E2F1 in prostate cancer. This review will aid in translating biomedical knowledge into therapeutic strategies for prostate cancer.