RFX5 promotes the progression of hepatocellular carcinoma through transcriptional activation of KDM4A

Targeting histone methylation and demethylation for non-alcoholic fatty liver disease

Non-alcoholic fatty liver disease (NAFLD) is the leading chronic liver disease worldwide, facing increasing challenges in terms of prevention and treatment. The methylation of lysine and arginine residues on histone proteins is dynamically controlled by histone methyltransferases (HMTs) and histone demethylases (HDMs), regulating chromatin structure and gene transcription. Mutations, genetic translocations, and altered gene expression involving HMTs and HDMs are frequently observed in NAFLD. HMTs and HDMs are receiving increasing attention in regulating NALFD. Targeting specific HMTs and HDMs for drug development is becoming a new strategy for treating NAFLD. This review provides a comprehensive summary of the regulatory mechanism of histone methylation/demethylation in NAFLD. Additionally, we discuss the potential applications of HMTs and HDMs inhibitors in preventing NAFLD, which may provide a scientific basis for the treatment of NAFLD.

The emerging roles of lysine-specific demethylase 4A in cancer: Implications in tumorigenesis and therapeutic opportunities

Lysine-specific demethylase 4 A (KDM4A, also named JMJD2A, KIA0677, or JHDM3A) is a demethylase that can remove methyl groups from histones H3K9me2/3, H3K36me2/3, and H1.4K26me2/me3. Accumulating evidence suggests that KDM4A is not only involved in body homeostasis (such as cell proliferation, migration and differentiation, and tissue development) but also associated with multiple human diseases, especially cancers. Recently, an increasing number of studies have shown that pharmacological inhibition of KDM4A significantly attenuates tumor progression in vitro and in vivo in a range of solid tumors and acute myeloid leukemia. Although there are several reviews on the roles of the KDM4 subfamily in cancer development and therapy, all of them only briefly introduce the roles of KDM4A in cancer without systematically summarizing the specific mechanisms of KDM4A in various physiological and pathological processes, especially in tumorigenesis, which greatly limits advances in the understanding of the roles of KDM4A in a variety of cancers, discovering targeted selective KDM4A inhibitors, and exploring the adaptive profiles of KDM4A antagonists. Herein, we present the structure and functions of KDM4A, simply outline the functions of KDM4A in homeostasis and non-cancer diseases, summarize the role of KDM4A and its distinct target genes in the development of a variety of cancers, systematically classify KDM4A inhibitors, summarize the difficulties encountered in the research of KDM4A and the discovery of related drugs, and provide the corresponding solutions, which would contribute to understanding the recent research trends on KDM4A and advancing the progression of KDM4A as a drug target in cancer therapy.

Identification and validation of STAT4 as a prognostic biomarker in acute myeloid leukemia

Acute myelogenous leukemia (AML) is a common malignancy and is supposed to have the ability to escape host immune surveillance. The present study aimed to identify key genes in AML that may affect tumor immunity and to provide prognosis biomarkers of AML. The Cancer Genome Atlas (TCGA) dataset was screened for transcription factors (TFs) involved in immunity and influencing survival, combining Gene Expression Omnibus (GEO) data to validate the impact on patient survival. A prognostic signature was established using four transcription factors, and these genes play an important role in the immune system, with higher regulatory T cell (Treg) scores in high-risk patients compared with the low-risk group. Analysis of individual genes showed that STAT4 and Treg are closely related, which may be due to STAT4 transcribing related genes that affect immunity. STAT4 expression was positively correlated with the proportion of abnormal cells and promoted AML recurrence as verified by AML clinical patient samples. In addition, silencing of STAT4 significantly slowed down the proliferation capacity of HL60 cells. In conclusion, these findings suggest that STAT4 may be a potential biomarker for AML prognosis. As a key gene affecting the prognosis of AML patients, STAT4 has the potential to be a candidate diagnostic and prognostic biomarker for AML.

Whole-genome resequencing reveals the uniqueness of Subei yak

Subei yak is an essential local yak in the Gansu Province, which genetic resource has recently been discovered. It is a meat-milk dual-purpose variety with high fecundity and relatively stable population genetic structure. However, its population genetic structure and genetic diversity are yet to be reported. Therefore, this study aimed to identify molecular markers of Subei yak genome by whole-genome resequencing, and to analyze the population structure and genetic diversity of Subei yak. This study screened 12,079,496 single nucleotide polymorphism (SNP) molecular markers in the 20 Subei yaks genome using whole-genome resequencing technology. Of these SNPs, 32.09% were located in the intronic region of the genome. Principal component analysis, phylogenetic analysis, and population structure analysis revealed that the Subei yak belonged to an independent group in the domestic yak population. A selective clearance analysis was carried out on Subei yak and other domestic yaks, and the genes under positive selection were annotated. The functional enrichment analysis showed that Subei yak possessed prominent selection characteristics in terms of external environment perception, hypoxia adaptation, and muscle development. Furthermore, Subei yak showed excellent muscle fat deposition and meat quality traits. Thus, this study will serve as a reference for discovering population structure, genetic evolution, and other unique traits of Subei yak and for expanding the genetic variation catalog of yaks.

Histone Demethylation Profiles in Nonalcoholic Fatty Liver Disease and Prognostic Values in Hepatocellular Carcinoma: A Bioinformatic Analysis

Nonalcoholic fatty liver disease (NAFLD) is the most common chronic liver disease with multifactorial pathogenesis; histone demethylases (HDMs) are emerging as attractive targets. We identified HDM genes (including KDM5C, KDM6B, KDM8, KDM4A, and JMJD7) that were differentially expressed in NAFLD and normal samples by exploring gene expression profiling datasets. There was no significant difference in the expression of genes related to histone demethylation between mild and advanced NAFLD. In vitro and in vivo studies indicated that KDM6B and JMJD7 were upregulated at the mRNA level in NAFLD. We explored the expression levels and prognostic values of the identified HDM genes in hepatocellular carcinoma (HCC). KDM5C and KDM4A were upregulated in HCC compared to normal tissue, while KDM8 showed downregulation. The abnormal expression levels of these HDMs could provide prognostic values. Furthermore, KDM5C and KDM4A were associated with immune cell infiltration in HCC. HDMs were associated with cellular and metabolic processes and may be involved in the regulation of gene expression. Differentially expressed HDM genes identified in NAFLD may provide value to understanding pathogenesis and in the development of epigenetic therapeutic targets. However, on the basis of the inconsistent results of in vitro studies, future in vivo experiments combined with transcriptomic analysis are needed for further validation.

Engineering human JMJD2A tudor domains for an improved understanding of histone peptide recognition

JMJD2A is a histone lysine demethylase which recognizes and demethylates H3K9me3 and H3K36me3 residues and is overexpressed in various cancers. It utilizes a tandem tudor domain to facilitate its own recruitment to histone sites, recognizing various di- and tri-methyl lysine residues with moderate affinity. In this study, we successfully engineered the tudor domain of JMJD2A to specifically bind to H4K20me3 with a 20-fold increase of affinity and improved selectivity. To reveal the molecular basis, we performed molecular dynamics and free energy decomposition analysis on the human JMJD2A tandem tudor domains bound to H4K20me2, H4K20me3, and H3K23me3 peptides to uncover the residues and conformational changes important for the enhanced binding affinity and selectivity toward H4K20me2/3. These analyses revealed new insights into understanding chromatin reader domains recognizing histone modifications and improving binding affinity and selectivity of these domains. Furthermore, we showed that the tight binding of JMJD2A to H4K20me2/3 is not sufficient to improve the efficiency of CRISPR-CAS9 mediated homology directed repair (HDR), suggesting a complicated relationship between JMJD2A and the DNA damage response beyond binding affinity toward the H4K20me2/3 mark.

Integrative network-based analysis on multiple Gene Expression Omnibus datasets identifies novel immune molecular markers implicated in non-alcoholic steatohepatitis

INTRODUCTION

Non-alcoholic steatohepatitis (NASH), an advanced subtype of non-alcoholic fatty liver disease (NAFLD), has becoming the most important aetiology for end-stage liver disease, such as cirrhosis and hepatocellular carcinoma. This study were designed to explore novel genes associated with NASH.

METHODS

Here, five independent Gene Expression Omnibus (GEO) datasets were combined into a single cohort and analyzed using network biology approaches.

RESULTS

11 modules identified by weighted gene co-expression network analysis (WGCNA) showed significant association with the status of NASH. Further characterization of four gene modules of interest demonstrated that molecular pathology of NASH involves the upregulation of hub genes related to immune response, cholesterol and lipid metabolic process, extracellular matrix organization, and the downregulation of hub genes related to cellular amino acid catabolic, respectively. After DEGs enrichment analysis and module preservation analysis, the Turquoise module associated with immune response displayed a remarkably correlation with NASH status. Hub genes with high degree of connectivity in the module, including CD53, LCP1, LAPTM5, NCKAP1L, C3AR1, PLEK, FCER1G, HLA-DRA and SRGN were further verified in clinical samples and mouse model of NASH. Moreover, single-cell RNA-seq analysis showed that those key genes were expressed by distinct immune cells such as microphages, natural killer, dendritic, T and B cells. Finally, the potential transcription factors of Turquoise module were characterized, including NFKB1, STAT3, RFX5, ILF3, ELF1, SPI1, ETS1 and CEBPA, the expression of which increased with NASH progression.

DISCUSSION

In conclusion, our integrative analysis will contribute to the understanding of NASH and may enable the development of potential biomarkers for NASH therapy.

Cantharidin Inhibits Proliferation of Liver Cancer by Inducing DNA Damage via KDM4A-Dependent Histone H3K36 Demethylation

Objective

To investigate the effect of cantharidin on DNA damage in hepatocellular carcinoma cells and its possible mechanism.

Methods

Cell proliferation assay and terminal deoxynucleotidyl transferase-mediated dUTP nick end labeling (TUNEL) assay were used to analyze the effects of cantharidin on cell proliferation and apoptosis of hepatocellular carcinoma cells. The expression levels of DNA damage markers H2AX and P21 were analyzed by qRT-PCR. The expression of KDM4A and H3K36me3 was observed by western blot. The expression of KDM4A was regulated by siRNA or plasmid transfection. The effect of KDM4A on DNA damage induced by cantharidin in liver cancer was observed after overexpression and addiction of KDM4A.

Results

Cantharidin can significantly inhibit the growth of hepatocellular carcinoma cells and induce apoptosis of hepatocellular carcinoma cells. Cantharidin enhances the chemotherapy sensitivity of liver cancer by targeting the upregulation of KDM4A and the regulation of DNA damage induced by H3K36me3. Overexpression of KDM4A enhances DNA damage induced by cantharidin in HCC. KDM4A silencing attenuated the damage of cantharidin to the DNA of HCC cells.

Conclusion

Cantharidin can inhibit the growth and promote apoptosis of hepatocellular carcinoma cells. Meanwhile, cantharidin can induce DNA damage in HCC cells. Mechanism studies have shown that cantharidin induces DNA damage through the demethylation of KDM4A-dependent histone H3K36.

Expression of IER3 in hepatocellular carcinoma: clinicopathology, prognosis, and potential regulatory pathways

Background

Immediate early response 3 (IER3) is correlated to the prognosis of several cancers, but the precise mechanisms underlying the regulation by IER3 of the occurrence and development of hepatocellular carcinoma (HCC) remain unknown.

Methods

The expression level of IER3 was examined by using in-house immunohistochemistry (IHC), public gene chip, and public RNA-sequencing (RNA-seq). The standardized mean difference (SMD) was calculated to compare the expression levels of IER3 between HCC patients and controls. The summary receiver operating characteristics (sROC) was plotted to comprehensively understand the discriminatory capability of IER3 between HCC and non-HCC group. The Kaplan-Meier curves and the combined hazard ratios (HRs) were used to determine the prognostic value of IER3 in HCC. Moreover, differentially expressed genes (DEGs) and co-expression genes (CEGs) were used to explored the molecular mechanisms of IER3 underlying HCC. hTFtarget was used to predict the transcription factors (TFs) of IER3. The binding site of TFs and the IER3 promoter region was forecasted using the JASPAR website. The relevant ChIP-seq data were used to determine whether TF peaks were present in the IER3 transcription initiation.

Results

A significantly increased expression of IER3 protein was found in HCC tissue relative to non-HCC tissue as detected by IHC (p < 0.001). Compared to 1,263 cases of non-HCC tissues, IER3 in 1483 cases of HCC tissues was upregulated (SMD = 0.42, 95% confidence interval [CI] [0.09-0.76]). The sROC showed that IER3 had a certain ability at differentiating HCC tissues (area under the curve (AUC) = 0.65, 95% CI [0.61-0.69]). Comprehensive analysis of the effect of IER3 on the prognosis of patients with HCC demonstrated that higher IER3 expression was associated with poor prognosis in HCC (HRs = 1.30, 95% CI [1.03-1.64]). Pathway enrichment analysis revealed that IER3-related genes were mostly enriched in the PI3K-Akt signaling pathway, cancer-related signaling pathways, the p53 signaling pathway, and other signaling pathways. Regulatory factor X5 (RFX5) was identified as a possible regulator of IER3-related TF.

Conclusion

IER3 may be a potential prognostic marker for HCC. The molecular mechanisms of IER3 in HCC warrant further study.

The transcription factor RFX5 positively regulates expression of MHCIa in the red-spotted grouper (Epinephelus akaara)

Regulatory factor X 5 (RFX 5) is a member of the RFX family, and it forms the transcription factor complex RFX with RFXANK/B and RFXAP. The RFX complex can activate MHC expression by binding to the MHC promoter. However, the regulate mechanism of RFX in fish species is not been fully elucidated. In this study, we investigated the transcriptional regulation of Epinephelus akaara RFX5 (EaRFX5) on EaMHCI, and its effect on immune pathways. The genomic sequence of EaRFX5 was 35,774 bp and consisted of ten exons and nine introns. The length of EaRFX5 ORF sequence is 2,160 bp, encoding 719 amino acids. By qRT-PCR, EaRFX5 was detected constitutively expressed in twelve selected tissues, showing a wide range of expression. EaRFX5 expression parttern in response to poly (I:C), LPS, Zymosan A, SGIV, and NNV challenges showed that EaRFX5 plays a differentiated immunomodulatory role in response to various stimuli in different tissues, and EaRFX5 was most significantly upregulated in the kidney after challenge with SGIV. Subcellular localization assays showed that EaRFX5 is a typical nuclear protein. Based on the in vitro overexpression experiments, EaRFX5 appeared to promote the expression of EaMHCIa gene, interferon signalling pathway and inflammatory cytokine. Luciferase reporter assay showed that the -267 bp to +82 bp region of EaMHCIa promoter was the core region where EaRFX5 modulated. Additionally, point mutations and electrophoretic mobility shift assays indicating M3 is the EaRFX5 binding sites in the EaMHCIa promoter. These results contribute to elucidating the function of EaRFX5 in fish immune response, and provide the first evidence of positive regulation of MHCIa expression by RFX5 in fish.

PHOSPHO1 Serves as a Key Metabolism-Related Biomarker in the Tumorigenesis of Diffuse Large B-cell Lymphoma

OBJECTIVE

Diffuse large B-cell lymphoma (DLBCL) is an aggressive type of non-Hodgkin lymphoma. Due to its genetic heterogeneity and abnormal metabolism, many DLBCL patients have a poor prognosis. This study investigated the key metabolism-related genes and potential mechanisms.

METHODS

Differentially expressed genes, differentially expressed transcription factors (TFs), and differentially expressed metabolism-related genes (DEMRGs) of glucose and lipid metabolic processes were identified using the edgeR package. Key DEMRGs were screened by Lasso regression, and a prediction model was constructed. The cell type identification by estimating relative subsets of RNA transcripts algorithm was utilized to assess the fraction of immune cells, and Gene Set Enrichment Analysis was used to determine immune-related pathways. A regulatory network was constructed with significant co-expression interactions among TFs, DEMRGs, immune cells/pathways, and hallmark pathways.

RESULTS

A total of 1551 DEMRGs were identified. A prognostic model with a high applicability (area under the curve=0.921) was constructed with 13 DEMRGs. Tumorigenesis of DLBCL was highly related to the neutrophil count. Four DEMRGs (PRXL2AB, CCN1, DECR2 and PHOSPHO1) with 32 TF-DEMRG, 36 DEMRG-pathway, 14 DEMRG-immune-cell, 9 DEMRG-immune-gene-set, and 67 DEMRG-protein-chip interactions were used to construct the regulatory network.

CONCLUSION

We provided a prognostic prediction model based on 13 DEMRGs for DLBCL. We found that phosphatase, orphan 1 (PHOSPHO1) is positively regulated by regulatory factor X5 (RFX5) and mediates MYC proto-oncogene (MYC) targeting the V2 pathway and neutrophils.

Analysis of m6A-Related lncRNAs for Prognosis Value and Response to Immune Checkpoint Inhibitors Therapy in Hepatocellular Carcinoma

Introduction

N6-methyladenosine (m6A) modification and long non-coding RNAs (lncRNAs) play pivotal roles in the progression of hepatocellular carcinoma (HCC). However, how their interaction is involved in the prognostic value of HCC and immune checkpoint inhibitors (ICIs) therapy remains unclear.

Methods

The RNA sequencing and clinical data of HCC patients were collected from TCGA database. The prognostic m6A-related lncRNAs were screened out with Pearson correlation test, univariate Cox analysis and least absolute shrinkage and selection operator (LASSO) Cox regression. Patients with HCC were classified into 2 subtypes by consensus clustering. Survival analyses were performed to assess the prognostic value of different clusters and risk models. Potential tumor correlated biological pathways correlated with different clusters were explored through gene set enrichment analysis. We also identified the relationship of the risk model and clusters with response to immune checkpoint inhibitors (ICIs) therapy and tumor microenvironment (TME). Furthermore, the prognostic value of the 9 m6A-related lncRNAs was validated in the external cohort. Finally, the role of SNHG4 was explored by silencing and overexpression of SNHG4 through conducting proliferation, migration and invasion experiments.

Results

Patients from 2 clusters and different risk groups based on m6A-related lncRNAs had significantly different clinicopathological characteristics and overall survival outcomes. Tumor-correlated biological pathways were found to be correlated with Cluster 2 through GSEA. Moreover, we found that patients from different clusters and risk groups expressed higher levels of immune checkpoint genes and had distinct TME and different responses for ICIs therapy. Prognostic value of this risk model was further confirmed in the external cohort. Finally, consistent with the discovery, SNHG4 played an oncogenic role in vitro.

Conclusion

Our study demonstrated that the 9 m6A-related lncRNA signature may serve as a novel predictor in the prognosis of HCC and optimize (ICIs) therapy. SNHG4 plays an oncogenic role in HCC.

tRNA-Derived Fragment tRF-17-79MP9PP Attenuates Cell Invasion and Migration via THBS1/TGF-β1/Smad3 Axis in Breast Cancer

tRNA derivatives have been identified as a new kind of potential biomarker for cancer. Previous studies have identified that there were 30 differentially expressed tRNAs derivatives in breast cancer tissue with the high-throughput sequencing technique. This study aimed to investigate the possible biological function and mechanism of tRNA derivatives in breast cancer cells. One such tRF, a 5'-tRF fragment of tRF-17-79MP9PP (tRF-17) was screened in this study, which is processed from the mature tRNA-Val-AAC and tRNA-Val-CAC. tRF-17 with significantly low expression in breast cancer tissues and serum. The level of tRF-17 differentiated breast cancer from healthy controls with sensitivity of 70.4% and specificity of 68.4%. Overexpression of tRF-17 suppressed cells malignant activity. THBS1 (Thrombospondin-1) as a downstream target of tRF-17, and reduction of THBS1 expression also partially recovered the effects of tRF-17 inhibition on breast cancer cell viability, invasion and migration. Besides, THBS1, TGF-β1, Smad3, p-Smad3 and epithelial-to-mesenchymal transition related genes N-cadherin, MMP3, MMP9 were markedly down-regulated in tRF-17 overexpressing cells. Moreover, tRF-17 attenuated the THBS1-mediated TGF-β1/Smad3 signaling pathway in breast cancer cells. In general, the tRF-17/THBS1/TGF-β1/smad3 axis elucidates the molecular mechanism of breast cancer cells invasion and migration and could lead to a potential therapeutic target for breast cancer.

Identification of a 14-Gene Prognostic Signature for Diffuse Large B Cell Lymphoma (DLBCL)

Although immunotherapy is a potential strategy to resist cancers, due to the inadequate acknowledge, this treatment is not always effective for diffuse large B cell lymphoma (DLBCL) patients. Based on the current situation, it is critical to systematically investigate the immune pattern. According to the result of univariate and multivariate cox proportional hazards, LASSO regression and Kaplan-Meier survival analysis on immune-related genes (IRGs), a prognostic signature, containing 14 IRGs (AQP9, LMBR1L, FGF20, TANK, CRP, ORM1, JAK1, BACH2, MTCP1, IFITM1, TNFSF10, FGF12, RFX5, and LAP3), was built. This model was validated by external data, and performed well. DLBCL patients were divided into low- and high-risk groups, according to risk scores from risk formula. The results of CIBERSORT showed that different immune status and infiltration pattern were observed in these two groups. Gene set enrichment analysis (GSEA) indicated 12 signaling pathways were significantly enriched in the high-risk group, such as natural killer cell-mediated cytotoxicity, toll-like receptor signaling pathway, and so on. In summary, 14 clinically significant IRGs were screened to build a risk score formula. This formula was an accurate tool to provide a certain basis for the treatment of DLBCL patients.