BAZ2A (TIP5) is involved in epigenetic alterations in prostate cancer and its overexpression predicts disease recurrence

Structure- and Property-Based Optimization of Efficient Pan-Bromodomain and Extra Terminal Inhibitors to Identify Oral and Intravenous Candidate I-BET787

The bromodomain and extra terminal (BET) family of bromodomain-containing proteins are important epigenetic regulators that elicit their effect through binding histone tail N-acetyl lysine (KAc) post-translational modifications. Recognition of such markers has been implicated in a range of oncology and immune diseases and, as such, small-molecule inhibition of the BET family bromodomain-KAc protein-protein interaction has received significant interest as a therapeutic strategy, with several potential medicines under clinical evaluation. This work describes the structure- and property-based optimization of a ligand and lipophilic efficient pan-BET bromodomain inhibitor series to deliver candidate I-BET787 (70) that demonstrates efficacy in a mouse model of inflammation and suitable properties for both oral and intravenous (IV) administration. This focused two-phase explore-exploit medicinal chemistry effort delivered the candidate molecule in 3 months with less than 100 final compounds synthesized.

An RNA-dependent and phase-separated active subnuclear compartment safeguards repressive chromatin domains

The nucleus is composed of functionally distinct membraneless compartments that undergo phase separation (PS). However, whether different subnuclear compartments are connected remains elusive. We identified a type of nuclear body with PS features composed of BAZ2A that associates with active chromatin. BAZ2A bodies depend on RNA transcription and BAZ2A non-disordered RNA-binding TAM domain. Although BAZ2A and H3K27me3 occupancies anticorrelate in the linear genome, in the nuclear space, BAZ2A bodies contact H3K27me3 bodies. BAZ2A-body disruption promotes BAZ2A invasion into H3K27me3 domains, causing H3K27me3-body loss and gene upregulation. Weak BAZ2A-RNA interactions, such as with nascent transcripts, promote BAZ2A bodies, whereas the strong binder long non-coding RNA (lncRNA) Malat1 impairs them while mediating BAZ2A association to chromatin at nuclear speckles. In addition to unraveling a direct connection between nuclear active and repressive compartments through PS mechanisms, the results also showed that the strength of RNA-protein interactions regulates this process, contributing to nuclear organization and the regulation of chromatin and gene expression.

The Transcriptional Landscape of Immune-Response 3'-UTR Alternative Polyadenylation in Melanoma

The prognosis of patients with malignant melanoma has been improved in recent decades due to advancements in immunotherapy. However, a considerable proportion of patients are refractory to treatment, particularly at advanced stages. This underscores the necessity of developing a new strategy to improve it. Alternative polyadenylation (APA), as a marker of crucial posttranscriptional regulation, has emerged as a major new type of epigenetic marker involved in tumorigenesis. However, the potential roles of APA in shaping the tumor microenvironment (TME) are largely unexplored. Herein, we collected two cohorts comprising melanoma patients who received immune checkpoint inhibitor (ICI) immunotherapy to quantify transcriptome-wide discrepancies in APA. We observed a global change in 3'-UTRs between responders and non-responders, which might involve DNA damage response, angiogenesis, PI3K-AKT signaling pathways, etc. Ten putative master APA regulatory factors for those APA events were detected via a network analysis. Notably, we established an immune response-related APA scoring system (IRAPAss), which exhibited a great performance of predicting immunotherapy response in multiple cohorts. Furthermore, we examined the correlation of APA with TME at the single-cell level using four single-cell immune profiles of tumor-infiltrating lymphocytes (TILs), which revealed an overall discrepancy in 3'-UTR length across diverse T cell populations, probably contributing to immunoregulation in melanoma. In conclusion, our study provides a transcriptional landscape of APA implicated in immunoregulation, which might lay the foundation for developing a new strategy for improving immunotherapy response for melanoma patients by targeting APA.

Pan-cancer and multi-omics analyses revealed the diagnostic and prognostic value of BAZ2A in liver cancer

BAZ2A, an epigenetic regulatory factor that affects ribosomal RNA transcription, has been shown to be highly expressed in several cancers and promotes tumor cell migration. This study explored the expression and mechanism of BAZ2A in tumorigenesis at the pan-cancer level. The Cancer Genome Atlas, Gene Expression Omnibus databases and TIMER2.0, cBioPortal and other tools were used to analyze the level of expression of BAZ2A in various tumor tissues and to examine the relationship between BAZ2A and survival, prognosis, mutation and immune invasion. In vitro experiments were performed to assess the function of BAZ2A in cancer cells. Using combined transcriptome and proteome analysis, we examined the possible mechanism of BAZ2A in tumors. BAZ2A exhibited high expression levels in multiple tumor tissues and displayed a significant association with cancer patient prognosis. The main type of BAZ2A genetic variation in cancer is gene mutation. Downregulation of BAZ2A inhibited proliferation, migration, and invasion and promoted apoptosis in LM6 liver cancer cell. The mechanism of BAZ2A in cancer development may involve lipid metabolism. These results help expand our understanding of BAZ2A in tumorigenesis and development and suggest BAZ2A may serve as a prognostic and diagnostic factor in several cancers.

Analyzing omics data based on sample network

Identifying valuable features from complex omics data is of great significance for disease diagnosis study. This paper proposes a new feature selection algorithm based on sample network (FS-SN) to mine important information from omics data. The sample network is constructed according to the sample neighbor relationship at the molecular (feature) expression level, and the distinguishing ability of the feature is evaluated based on the topology of the sample network. The sample network established on a feature with a strong discriminating ability tends to have many edges between the same group samples and few edges between the different group samples. At the same time, FS-SN removes redundant features according to the gravitational interaction between features. To show the validation of FS-SN, it was compared on ten public datasets with ERGS, mRMR, ReliefF, ATSD-DN, and INDEED which are efficient in omics data analysis. Experimental results show that FS-SN performed better than the compared methods in accuracy, sensitivity and specificity in most cases. Hence, FS-SN making use of the topology of the sample network is effective for analyzing omics data, it can identify key features that reflect the occurrence and development of diseases, and reveal the underlying biological mechanism.

Biochemical and cellular insights into the Baz2B protein, a non-catalytic subunit of the chromatin remodeling complex

Baz2B is a regulatory subunit of the ATP-dependent chromatin remodeling complexes BRF1 and BRF5, which control access to DNA during DNA-templated processes. Baz2B has been implicated in several diseases and also in unhealthy ageing, however limited information is available on the domains and cellular roles of Baz2B. To gain more insight into the Baz2B function, we biochemically characterized the TAM (Tip5/ARBP/MBD) domain with the auxiliary AT-hook motifs and the bromodomain (BRD). We observed alterations in histone code recognition in bromodomains carrying cancer-associated point mutations, suggesting their potential involvement in disease. Furthermore, the depletion of Baz2B in the Hap1 cell line resulted in altered cell morphology, reduced colony formation and perturbed transcriptional profiles. Despite that, super-resolution microscopy images revealed no changes in the overall chromatin structure in the absence of Baz2B. These findings provide insights into the biological function of Baz2B.

Nasopharyngeal carcinoma derived exosomes regulate the proliferation and migration of nasopharyngeal carcinoma cells by mediating the miR-99a-5p BAZ2A axis

OBJECTIVES

Nasopharyngeal Carcinoma (NPC) is a common malignant tumor of nasopharyngeal mucosal epithelium in clinical practice. Radiotherapy and chemotherapy are the main treatment methods at present, but the therapeutic effect is still unsatisfactory. Studies have shown that exosomes and microRNAs (miRNAs) play an important role in the development of cancer. Therefore, this study aimed to investigate the effects of NPC derived exosomes on NPC and their molecular mechanisms.

METHODS

Serum was collected from healthy subjects, Epstein-Barr Virus (EBV) infected patients and NPC patients (n = 9 group) and exosomes were extracted separately. High-throughput sequencing of exosomes was performed to screen differentially expressed miRNAs. The function of the screened miRNA was identified by treating NPC cells with exosomes. The target gene of miRNA was identified using the dual-luciferase assay. Real-Time quantitative Polymerase Chain Reaction (RT-qPCR) was used to determine the levels of miR-99a-5p and Bromodomain Adjacent Tozinc finger domain protein 2A (BAZ2A). Cell Counting Kit-8 assay, flow cytometry, and wound healing assay were utilized to detect cell viability, cell cycle and apoptosis, and migration ability. The protein levels were evaluated by Western blot.

RESULTS

MiR-99a-5p was identified as the most significant differentially expressed miRNA in exosomes (p < 0.05). The proliferation and migration of NPC cells were extremely facilitated by exosomes, accompanied by the suppressed apoptosis, upregulated BAZ2A, Monocyte Chemotactic Protein-1 (MCP1), and Vascular Endothelial Growth Factor A (VEGFA), and downregulation of Interleukin (IL)-1β and Nuclear Transcription Factor-κB (NF-κB) (p < 0.05). BAZ2A was a target gene of miR-99a-5p. Furthermore, the regulatory effect of exosomes on the proliferation, migration, and apoptosis was significantly abolished by overexpression of miR-99a-5p or downregulation of BAZ2A (p < 0.05).

CONCLUSION

NPC derived exosomes facilitated the proliferation and migration of NPC through regulating the miR-99a-5p/BAZ2A axis.

Cancer metastasis under the magnifying glass of epigenetics and epitranscriptomics

Most of the cancer-associated mortality and morbidity can be attributed to metastasis. The role of epigenetic and epitranscriptomic alterations in cancer origin and progression has been extensively demonstrated during the last years. Both regulations share similar mechanisms driven by DNA or RNA modifiers, namely writers, readers, and erasers; enzymes responsible of respectively introducing, recognizing, or removing the epigenetic or epitranscriptomic modifications. Epigenetic regulation is achieved by DNA methylation, histone modifications, non-coding RNAs, chromatin accessibility, and enhancer reprogramming. In parallel, regulation at RNA level, named epitranscriptomic, is driven by a wide diversity of chemical modifications in mostly all RNA molecules. These two-layer regulatory mechanisms are finely controlled in normal tissue, and dysregulations are associated with every hallmark of human cancer. In this review, we provide an overview of the current state of knowledge regarding epigenetic and epitranscriptomic alterations governing tumor metastasis, and compare pathways regulated at DNA or RNA levels to shed light on a possible epi-crosstalk in cancer metastasis. A deeper understanding on these mechanisms could have important clinical implications for the prevention of advanced malignancies and the management of the disseminated diseases. Additionally, as these epi-alterations can potentially be reversed by small molecules or inhibitors against epi-modifiers, novel therapeutic alternatives could be envisioned.

Reevaluation of bromodomain ligands targeting BAZ2A

BAZ2A promotes migration and invasion in prostate cancer. Two chemical probes, the specific BAZ2-ICR, and the BAZ2/BRD9 cross-reactive GSK2801, interfere with the recognition of acetylated lysines in histones by the bromodomains of BAZ2A and of its BAZ2B paralog. The two chemical probes were tested in prostate cancer cell lines with opposite androgen susceptibility. BAZ2-ICR and GSK2801 showed different cellular efficacies in accordance with their unequal selectivity profiles. Concurrent inhibition of BAZ2 and BRD9 did not reproduce the effects observed with GSK2801, indicating possible off-targets for this chemical probe. On the other hand, the single BAZ2 inhibition by BAZ2-ICR did not phenocopy genetic ablation, demonstrating that bromodomain interference is not sufficient to strongly affect BAZ2A functionality and suggesting a PROTAC-based chemical ablation as an alternative optimization strategy and a possible therapeutic approach. In this context, we also present the crystallographic structures of BAZ2A in complex with the above chemical probes. Binding poses of TP-238 and GSK4027, chemical probes for the bromodomain subfamily I, and two ligands of the CBP/EP300 bromodomains identify additional headgroups for the development of BAZ2A ligands.

BAZ2A-RNA mediated association with TOP2A and KDM1A represses genes implicated in prostate cancer

BAZ2A represses rRNA genes (rDNA) that are transcribed by RNA polymerase I. In prostate cancer (PCa), BAZ2A function goes beyond this role because it represses genes frequently silenced in metastatic disease. However, the mechanisms of this BAZ2A-mediated repression remain elusive. Here, we show that BAZ2A represses genes through its RNA-binding TAM domain using mechanisms differing from rDNA silencing. Although the TAM domain mediates BAZ2A recruitment to rDNA, in PCa, this is not required for BAZ2A association with target genes. Instead, the BAZ2A-TAM domain in association with RNA mediates the interaction with topoisomerase 2A (TOP2A) and histone demethylase KDM1A, whose expression positively correlates with BAZ2A levels in localized and metastatic PCa. TOP2A and KDM1A pharmacological inhibition up-regulate BAZ2A-repressed genes that are regulated by inactive enhancers bound by BAZ2A, whereas rRNA genes are not affected. Our findings showed a novel RNA-based mechanism of gene regulation in PCa. Furthermore, we determined that RNA-mediated interactions between BAZ2A and TOP2A and KDM1A repress genes critical to PCa and may prove to be useful to stratify prostate cancer risk and treatment in patients.

A new prediction model of hepatocellular carcinoma based on N7-methylguanosine modification

PURPOSE

Hepatocellular carcinoma (HCC) is a kind of primary liver cancer. It is a common malignant tumor of digestive system that is difficult to predict the prognosis of patients. As an important epigenetic modification, N7 methyl guanosine (m7G) is indispensable in gene regulation. This regulation may affect the development and occurrence of cancer. However, the prognosis of long non coding RNAs (lncRNAs) in HCC is limited, especially how m7G-related lncRNAs regulate the development of HCC has not been reported.

METHODS

The Cancer Genome Atlas (TCGA) provides us with the expression data and corresponding clinical information of HCC patients we need. We used a series of statistical methods to screen four kinds of m7G-related lncRNAs related to HCC prognosis and through a series of verifications, the results were in line with our expectations. Finally, we also explored the IC50 difference and correlation analysis of various common chemotherapy drugs.

RESULT

Our study identified four differentially expressed m7g-related lncRNAs associated with HCC prognosis. Survival curve analysis showed that high risk lncRNAs would lead to poor prognosis of HCC patients. M7G signature's AUC was 0.789, which shows that the prognosis model we studied has certain significance in predicting the prognosis of HCC patients. Moreover, our study found that different risk groups have different immune and tumor related pathways through gene set enrichment analysis. In addition, many immune cell functions are significantly different among different risk groups, such as T cell functions, including coordination of type I INF response and coordination of type II INF response. The expression of PDCD1, HHLA2, CTLA-4 and many other immune checkpoints in different risk groups is also different. Additionally, we analyzed the differences of IC50 and risk correlation of 15 chemotherapeutic drugs among different risk groups.

CONCLUSION

A novel lncRNAs associated with m7G predicts the prognosis of HCC.

GSK2801 Reverses Paclitaxel Resistance in Anaplastic Thyroid Cancer Cell Lines through MYCN Downregulation

Anaplastic thyroid cancer (ATC) is a very rare, but extremely aggressive form of thyroid malignancy, responsible for the highest mortality rate registered for thyroid cancer. Treatment with taxanes (such as paclitaxel) is an important approach in counteracting ATC or slowing its progression in tumors without known genetic aberrations or those which are unresponsive to other treatments. Unfortunately, resistance often develops and, for this reason, new therapies that overcome taxane resistance are needed. In this study, effects of inhibition of several bromodomain proteins in paclitaxel-resistant ATC cell lines were investigated. GSK2801, a specific inhibitor of BAZ2A, BAZ2B and BRD9, was effective in resensitizing cells to paclitaxel. In fact, when used in combination with paclitaxel, it was able to reduce cell viability, block the ability to form colonies in an anchor-independent manner, and strongly decrease cell motility. After RNA-seq following treatment with GSK2801, we focused our attention on MYCN. Based on the hypothesis that MYCN was a major downstream player in the biological effects of GSK2801, we tested a specific inhibitor, VPC-70619, which showed effective biological effects when used in association with paclitaxel. This suggests that the functional deficiency of MYCN determines a partial resensitization of the cells examined and, ultimately, that a substantial part of the effect of GSK2801 results from inhibition of MYCN expression.

Identification of molecular subtypes and a prognostic signature based on chromatin regulators related genes in prostate cancer

The clinical and molecular phenotypes of prostate cancer (PCa) exhibit substantial heterogeneity, ranging from indolent to metastatic disease. In this study, we aimed to identify PCa subtypes and construct a gene signature that can predict the recurrence-free survival (RFS) of PCa patients based on chromatin regulators genes (CRGs). Strikingly, we identified two heterogeneous subtypes with distinct clinical and molecular characteristics. Furthermore, by performing differential analysis between the two CRGs subtypes, we successfully constructed a gene signature to predict PCa prognosis. The signature, comprising four genes (MXD3, SSTR1, AMH and PPFIA2), was utilized to classify PCa patients into two risk groups; the high-risk group was characterized by poor prognosis and more aggressive clinical features. Moreover, we investigated the immune profile, mutation landscape and molecular pathways in each of the groups. Additionally, drug-susceptibility testing was performed to explore sensitive drugs for high-risk patients. Furthermore, we found that MXD3 downregulation suppressed the proliferation of PCa cell lines in vitro. Overall, our results highlight the signature based on CRGs as a powerful tool for predicting RFS of PCa patients, as well as an indicator for personalized treatment of those patients.

Genetic Polymorphisms Associated with Prothrombin Time and Activated Partial Thromboplastin Time in Chinese Healthy Population

(1) Background: The purpose of this study was to evaluate the effect of gene polymorphisms on prothrombin time (PT) and activated partial thromboplastin time (APTT) in a healthy Chinese population. (2) Methods: A total of 403 healthy volunteers from a series of novel oral anticoagulants (NOACs) bioequivalence trials in China were included. Coagulation tests for PT and APTT were performed in the central lab at Peking University First Hospital. Whole-exome sequencing (WES) and genome-wide association analysis were performed. (3) Results: In the correlation analysis of PT, 105 SNPs from 84 genes reached the genome-wide significance threshold (p < 1 × 10−5). Zinc Finger Protein 594 (ZNF594) rs184838268 (p = 4.50 × 10−19) was most significantly related to PT, and Actinin Alpha 1 (ACTN1) was found to interact most with other candidate genes. Significant associations with previously reported candidate genes Aurora Kinase B (AURKB), Complement C5(C5), Clock Circadian Regulator (CLOCK), and Histone Deacetylase 9(HDAC9) were detected in our dataset (p < 1 × 10−5). PiggyBac Transposable Element Derived 2(PGBD2) rs75935520 (p = 4.49 × 10−6), Bromodomain Adjacent To Zinc Finger Domain 2A(BAZ2A) rs199970765 (p = 5.69 × 10−6) and Protogenin (PRTG) rs80064850 (p = 8.69 × 10−6) were significantly correlated with APTT (p < 1 × 10−5). The heritability values of PT and APTT were 0.83 and 0.64, respectively; (4) Conclusion: The PT and APTT of healthy populations are affected by genetic polymorphisms. ZNF594 and ACTN1 variants could be novel genetic markers of PT, while PRTG polymorphisms might be associated with APTT levels. The findings could be attributed to ethnic differences, and need further investigation.

Role of 5-methylcytosine in determining the prognosis, tumor microenvironment, and applicability of precision medicine in patients with hepatocellular carcinoma

Background: 5-methylcytosine has a profound impact on the development and progression of hepatocellular carcinoma. The aim of this study was to investigate the usefulness of 5-methylcytosine in determining the prognosis, tumor microenvironment, and applicability of precision medicine in hepatocellular carcinoma.

Methods: We collected data of seven hepatocellular carcinoma cohorts (The Cancer Genome Atlas, International Cancer Genome Consortium, GSE14520, GSE6764, GSE9843, GSE63898, GSE76427). An unsupervised clustering method was used to identify novel subtypes of hepatocellular carcinoma based on the expression 5-methylcytosine gene signatures. The 5-methylcytosine score was determined using the least absolute shrinkage and selection operator method based on the differential expression of genes in the identified subtypes. Subsequently, we investigated the association between 5-methylcytosine-based clusters (according to the 5-methylcytosine score) and clinical outcomes, immunophenotypes, classical molecular subtypes, and therapeutic opportunities in hepatocellular carcinoma. Finally, we examined the sensitivity of patients with high 5-methylcytosine score to drugs.

Results: We identified two hepatocellular carcinoma-specific, 5-methylcytosine-based subtypes (clusters 1 and 2). Cluster 1 exhibited significantly higher 5-methylcytosine scores versus cluster 2. The 5-methylcytosine-based subtypes accurately predicted classical molecular subtypes, immunophenotypes, prognosis, and therapeutic opportunities for patients with hepatocellular carcinoma. Cluster 1 (high 5-methylcytosine score) was characterized by lower anticancer immunity and worse prognosis versus cluster 2 (low 5-methylcytosine score). Moreover, cluster 1 (high 5-methylcytosine score) exhibited low sensitivity to cancer immunotherapy, but high sensitivity to radiotherapy and targeted therapy with lenvatinib.

Conclusion: The novel 5-methylcytosine-based subtypes (according to the 5-methylcytosine score) may reflect the prognosis, tumor microenvironment, and applicability of precision medicine in patients with hepatocellular carcinoma.

Characterization of chromatin regulators in hepatocellular carcinoma to guide clinical therapy

Background: Hepatocellular carcinoma (HCC) is notorious for its high mortality and incidence. Accumulating evidence confirms that chromatin regulators (CRs) have a significant impact on cancer. Therefore, exploring the mode of action and prognostic value of CRs is imminent for the treatment of hepatocellular carcinoma.

Method: Transcriptome and clinical data for this study have been downloaded from TCGA (https://portal.gdc.cancer.gov/) and ICGC (https://dcc.icgc.org/). Univariate analysis was used to screen CRs with prognostic value, and our prognostic risk score signature was developed using least absolute shrinkage along with selection operator (lasso) Cox regression analysis. The CRs-based prognostic model was constructed in the TCGA dataset, and low-risk HCC patients had a better prognosis, which was finally validated in the ICGC dataset. We used the receiver operating characteristic curve to identify the accuracy of the prediction model and establish a line chart to prove the clinical effectiveness of the model. We also discussed the differences in drug sensitivity via CellMiner database, tumor immune microenvironment via ssGSEA algorithm, and clinical characteristics among different risk groups.

Results: A prognostic model consisting of seven CRs was constructed and verified in HCC patients. Furthermore, we found that this risk score prognostic signature could independently predict the prognosis of HCC patients. Functional enrichment analysis revealed that CRs are mainly associated with cancer-related signaling pathways and metabolic pathways. In addition, immune cell abundance correlates with risk score levels

Conclusion: In brief, we systematically explored the mode of action of CRs in HCC patients and established a reliable prognostic prediction model.

Molecular characterization, clinical relevance and immune feature of m7G regulator genes across 33 cancer types

Over 170 RNA modifications have been identified after transcriptions, involving in regulation of RNA splicing, processing, translation and decay. Growing evidence has unmasked the crucial role of N6-methyladenosine (m6A) in cancer development and progression, while, as a relative newly found RNA modification, N7-methylguanosine (m7G) is also certified to participate in tumorigenesis via different catalytic machinery from that of m6A. However, system analysis on m7G RNA modification-related regulator genes is lack. In this study, we first investigated the genetic alteration of m7G related regulator genes in 33 cancers, and found mRNA expression levels of most regulator genes were positively correlated with copy number variation (CNV) and negatively correlated with methylation in most cancers. We built a m7G RNA modification model based on the enrichment of the regulator gene scores to evaluate the m7G modification levels in 33 cancers, and investigated the connections of m7G scores to clinical outcomes. Furthermore, we paid close attention to the role of m7G in immunology due to the widely used immune checkpoint blockade therapy. Our results showed the higher m7G scores related to immunosuppression of tumor cells. Further confirmation with phase 3 clinical data with application of anti-PDL1/PDL indicated the impact of m7G modification level on immunotherapy effect. Relevance of m7G regulator genes and drug sensitivity was also evaluated to provide a better treatment choice when treating cancers. In summary, our study uncovered the profile of m7G RNA modification through various cancers, and figured out the connection of m7G modification levels with therapeutical outcomes, providing potential better options of cancer treatment.

dN/dS-H, a New Test to Distinguish Different Selection Modes in Protein Evolution and Cancer Evolution

One of the most popular measures in the analysis of protein sequence evolution is the ratio of nonsynonymous distance (d_N) to synonymous distance (d_S). Under the assumption that synonymous substitutions in the coding region are selectively neutral, the d_N/d_S ratio can be used to statistically detect the adaptive evolution (or purifying selection) if d_N/d_S > 1 (or d_N/d_S < 1) significantly. However, due to strong structural constraints and/or variable functional constraints imposed on amino acid sites, most encoding genes in most species have demonstrated d_N/d_S < 1. Consequently, the statistical power for testing d_N/d_S = 1 may be insufficient to distinguish between different selection modes. In this paper, we propose a more powerful test, called d_N/d_S-H, in which a new parameter H, a relative measure of rate variation among sites, was introduced. Given the condition of strong purifying selections at some sites, the d_N/d_S-H model predicts d_N/d_S = 1-H for neutral evolution, d_N/d_S < 1-H for nearly neutral selection, and d_N/d_S > 1-H for adaptive evolution. The potential of this new method for resolving the neutral-adaptive debates is illustrated by the protein sequence evolution in vertebrates, Drosophila and yeasts, as well as somatic cancer evolution (specialized as the C_N/C_S-H test).

Identification of a BAZ2A-Bromodomain Hit Compound by Fragment Growing

BAZ2A is an epigenetic regulator affecting transcription of ribosomal RNA. It is overexpressed in aggressive and recurrent prostate cancer, promoting cellular migration. Its bromodomain is characterized by a shallow and difficult-to-drug pocket. Here, we describe a structure-based fragment-growing campaign for the identification of ligands of the BAZ2A bromodomain. By combining docking, competition binding assays, and protein crystallography, we have extensively explored the interactions of the ligands with the rim of the binding pocket, and in particular ionic interactions with the side chain of Glu1820, which is unique to BAZ2A. We present 23 high-resolution crystal structures of the holo BAZ2A bromodomain and analyze common bromodomain/ligand motifs and favorable intraligand interactions. Binding of some of the compounds is enantiospecific, with affinity in the low micromolar range. The most potent ligand has an equilibrium dissociation constant of 7 μM and a good selectivity over the paralog BAZ2B bromodomain.

A Simple Evolutionary Model of Genetic Robustness After Gene Duplication

When a dispensable gene is duplicated (referred to the ancestral dispensability denoted by O⁺), genetic buffering and duplicate compensation together maintain the duplicate redundancy, whereas duplicate compensation is the only mechanism when an essential gene is duplicated (referred to the ancestral essentiality denoted by O^-). To investigate these evolutionary scenarios of genetic robustness, I formulated a simple mixture model for analyzing duplicate pairs with one of the following states: double dispensable (DD), semi-dispensable (one dispensable one essential, DE), or double essential (EE). This model was applied to the yeast duplicate pairs from a whole-genome duplication (WGD) occurred about 100 million years ago (mya), and the mouse duplicate pairs from a WGD occurred about more than 500 mya. Both case studies revealed that the proportion of essentiality for those duplicates with ancestral essentiality [P_E(O^-)] was much higher than that for those with ancestral dispensability [P_E(O⁺)]. While it was negligible in the yeast duplicate pairs, P_E(O⁺) (about 20%) was shown statistically significant in the mouse duplicate pairs. These findings, together, support the hypothesis that both sub-functionalization and neo-functionalization may play some roles after gene duplication, though the former may be much faster than the later.

A Systematic Analysis Reveals the Prognostic and Immunological Role of Reptin/RUVBL2 in Human Tumors

Reptin/RUVBL2 is involved in the remodeling of chromatin, DNA damage repair, and regulation of the cell cycle, all of which help to play essential roles in cancer. However, relevant pan-cancer analysis of Reptin is lacking. This study first investigated the potential oncogenic roles of Reptin and revealed a relationship between Reptin with clinicopathological characteristics and immune infiltration based on big data. Here, we showed that Reptin is overexpressed in many cancers. A significant association exists between the expression of Reptin and the prognosis of cancer cases. Reptin had a meaningful interaction with the immune infiltration of CD4⁺ Th1 cells and immune modulator genes in multiple cancer types. And negative correlation exists between Reptin and cancer-associated fibroblasts in BRCA, PRAD, TGCT, and THYM. A significant negative association exists between Reptin and regulatory T cells in TGCT and THCA. Moreover, Reptin is significantly associated with genomic heterogeneity, DNA mismatch repair genes, methyltransferase, and RNA modification genes in specific cancer types. Spliceosome, Hippo signaling pathway, DNA replication pathway, and acetyltransferase activity-associated functions were observed in the effect of Reptin on the tumor. This systematic analysis highlights Reptin as a vital cancer regulator among numerous genes and proved its potential prognosticator value and therapeutic target role for specific tumor types.

Metabolism-Associated DNA Methylation Signature Stratifies Lower-Grade Glioma Patients and Predicts Response to Immunotherapy

Metabolism and DNA methylation (DNAm) are closely linked. The value of the metabolism-DNAm interplay in stratifying glioma patients has not been explored. In the present study, we aimed to stratify lower-grade glioma (LGG) patients based on the DNAm associated with metabolic reprogramming. Four data sets of LGGs from three databases (TCGA/CGGA/GEO) were used in this study. By screening the Kendall’s correlation of DNAm with 87 metabolic processes from KEGG, we identified 391 CpGs with a strong correlation with metabolism. Based on these metabolism-associated CpGs, we performed consensus clustering and identified three distinct subgroups of LGGs. These three subgroups were characterized by distinct molecular features and clinical outcomes. We also constructed a subgroup-related, quantifiable CpG signature with strong prognostic power to stratify LGGs. It also serves as a potential biomarker to predict the response to immunotherapy. Overall, our findings provide new perspectives for the stratification of LGGs and for understanding the mechanisms driving malignancy.

Uncovering N4-Acetylcytidine-Related mRNA Modification Pattern and Landscape of Stemness and Immunity in Hepatocellular Carcinoma

N4-acetylcytidine (ac4C) is an ancient and conserved RNA modification. Previously, ac4C mRNA modification has been reported promoting proliferation and metastasis of tumor cells. However, it remains unclear whether and how ac4C-related mRNA modification patterns influencing the prognosis of hepatocellular carcinoma (HCC) patients. Hereby, we constructed an ac4Cscore model and classified patients into two groups and investigated the potential intrinsic and extrinsic characteristics of tumor. The ac4Cscore model, including COL15A1, G6PD and TP53I3, represented ac4C-related mRNA modification patterns in HCC. According to ac4Cscore, patients were stratified to high and low groups with distinct prognosis. Patients subject to high group was related to advanced tumor stage, higher TP53 mutation rate, higher tumor stemness, more activated pathways in DNA-repair system, lower stromal score, higher immune score and higher infiltrating of T cells regulatory. While patients attributed to low group were correlated with abundance of T cells CD4 memory, less aggressive immune subtype and durable therapy benefit. We also found ac4Cscore as a novel marker to predict patients’ prognosis with anti-PD1 immunotherapy and/or mTOR inhibitor treatment. Our study for the first time showed the association between ac4C-related mRNA modification patterns and tumor intrinsic and extrinsic characteristics, thus influencing the prognosis of patients.

Construction of a Novel Immune-Related mRNA Signature to Predict the Prognosis and Immune Characteristics of Human Colorectal Cancer

Background: Anti-cancer immunotherapeutic approaches have gained significant efficacy in multiple cancer types. However, not all patients with colorectal cancer (CRC) could benefit from immunotherapy due to tumor heterogeneity. The purpose of this study was to construct an immune-related signature for predicting the immune characteristics and prognosis of CRC.

Methods: RNA-sequencing data and corresponding clinical information of patients with CRC were obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO), and immune-related genes (IRGs) were downloaded from the Immunology Database and Analysis Portal (ImmPort). Then, we utilized univariate, lasso regression, and multivariate cox regression to identify prognostic IRGs and develop the immune-related signature. Subsequently, a nomogram was established based on the signature and other prognostic factors, and its predictive capacity was assessed by receiver operating characteristic (ROC) and decision curve analysis (DCA). Finally, associations between the signature and the immune characteristics of CRC were assessed.

Results: In total, 472 samples downloaded from TCGA were divided into the training cohort (236 samples) and internal validation cohort (236 samples), and the GEO cohort was downloaded as an external validation cohort (122 samples). A total of 476 differently expressed IRGs were identified, 17 of which were significantly correlated to the prognosis of CRC patients. Finally, 10 IRGs were filtered out to construct the risk score signature, and patients were divided into low- and high-risk groups according to the median of risk scores in the training cohort. The high-risk score was significantly correlated with unfavorable survival outcomes and aggressive clinicopathological characteristics in CRC patients, and the results were further confirmed in the internal validation cohort, entire TCGA cohort, and external validation cohort. Immune infiltration analysis revealed that patients in the low-risk group infiltrated with high tumor-infiltrating immune cell (TIIC) abundances compared to the high-risk group. Moreover, we also found that the immune checkpoint biomarkers were significantly overexpressed in the low-risk group.

Conclusion: The prognostic signature established by IRGs showed a promising clinical value for predicting the prognosis and immune characteristics of human CRC, which contribute to individualized treatment decisions.

Generation and Analysis of Pyroptosis-Based and Immune-Based Signatures for Kidney Renal Clear Cell Carcinoma Patients, and Cell Experiment

Background: Pyroptosis is a programmed cell death caused by inflammasomes, which is closely related to immune responses and tumor progression. The present study aimed to construct dual prognostic indices based on pyroptosis-associated and immune-associated genes and to investigate the impact of the biological signatures of these genes on Kidney Renal Clear Cell Carcinoma (KIRC).

Materials and Methods: All the KIRC samples from the Cancer Genome Atlas (TCGA) were randomly and equally divided into the training and testing datasets. Cox and Least Absolute Shrinkage and Selection Operator (LASSO) regression analysis were used to screen crucial pyroptosis-associated genes (PAGs), and a pyroptosis-associated genes prognostic index (PAGsPI) was constructed. Immune-associated genes (IAGs) related to PAGs were identified, and then screened through Cox and LASSO regression analyses, and an immune-associated genes prognostic index (IAGsPI) was developed. These two prognostic indices were verified by using the testing and the Gene Expression Omnibus (GEO) datasets and an independent cohort. The patients’ response to immunotherapy was analyzed. A nomogram was constructed and calibrated. qRT-PCR was used to detect the expression of PAGs and IAGs in the tumor tissues and normal tissues. Functional experiment was carried out.

Results: 86 PAGs and 1,774 differentially expressed genes (DEGs) were obtained. After intersecting PAGs with DEGs, 22 differentially expressed PAGs (DEPAGs) were included in Cox and LASSO regression analyses, identifying 5 crucial PAGs. The PAGsPI was generated. Patients in the high-PAGsPI group had a poor prognosis. 82 differentially expressed IAGs (DEIAGs) were highly correlated with DEPAGs. 7 key IAGs were screened out, and an IAGsPI was generated. Patients in the high-IAGsPI group had a poor prognosis. PAGsPI and IAGsPI were verified to be robust and reliable. The results revealed patients in low-PAGsPI group and high-IAGsPI group may be more sensitive to immunotherapy. The calibrated nomogram was proved to be reliable. An independent cohort study also proved that PAGsPI and IAGsPI performed well in prognosis prediction. We found that the expression of AIM2 may affect proliferation of KIRC cells.

Conclusion: PAGsPI and IAGsPI could be regarded as potential biomarkers for predicting the prognosis of patients with KIRC.

In vivo CRISPR screening identifies BAZ2 chromatin remodelers as druggable regulators of mammalian liver regeneration

Identifying new pathways that regulate mammalian regeneration is challenging due to the paucity of in vivo screening approaches. We employed pooled CRISPR knockout and activation screening in the regenerating liver to evaluate 165 chromatin regulatory proteins. Both screens identified the imitation-SWI chromatin remodeling components Baz2a and Baz2b, not previously implicated in regeneration. In vivo sgRNA, siRNA, and knockout strategies against either paralog confirmed increased regeneration. Distinct BAZ2-specific bromodomain inhibitors, GSK2801 and BAZ2-ICR, resulted in accelerated liver healing after diverse injuries. Inhibitor-treated mice also exhibited improved healing in an inflammatory bowel disease model, suggesting multi-tissue applicability. Transcriptomics on regenerating livers showed increases in ribosomal and cell cycle mRNAs. Surprisingly, CRISPRa screening to define mechanisms showed that overproducing Rpl10a or Rpl24 was sufficient to drive regeneration, whereas Rpl24 haploinsufficiency was rate limiting for BAZ2 inhibition-mediated regeneration. The discovery of regenerative roles for imitation-SWI components provides immediate strategies to enhance tissue repair.

Diagnostic, Therapeutic, and Prognostic Value of the m6A Writer Complex in Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) has poor prognosis and is usually diagnosed only at an advanced stage. Identification of novel biomarkers is critical to early diagnosis and better prognosis for HCC patients. N6-methyladenosine (m⁶A) RNA methylation regulators play important roles in the development of many tumors. However, the m⁶A writer complex, a key executor of m⁶A methylation modification, has not been independently investigated, and its specific bioinformatics analysis has not yet been performed in HCC. In this study, we used multiple public databases to evaluate the diagnostic, therapeutic, and prognostic value of the m⁶A writers in HCC. The results showed that expression levels of METTL3, VIRMA and CBLL1 were significantly increased, while expression levels of METTL14 and ZC3H13 were significantly decreased in HCC, which was closely related to clinicopathological factors, such as tumor stage and prognosis. Bioinformatics further explored the possible underlying mechanisms by which the m⁶A writer complex are involved in activation of tumor-promoting pathways and/or inhibition of tumor-suppressing pathways, including apoptosis, cell cycle, DNA damage response and EMT. Furthermore, we showed that the m⁶A writer complex is correlated with immune cell infiltration and immunoregulator expression in HCC. In conclusion, the m⁶A writer complex may represent a promising biomarker and target that can guide targeted therapy or immunotherapy for HCC patients.

Crosstalk Between Four Types of RNA Modification Writers Characterizes the Tumor Immune Microenvironment Infiltration Patterns in Skin Cutaneous Melanoma

The “writers” of four types of adenosine (A)-related RNA modifications (N6-methyladenosine, N1-methyladenosine, alternative polyadenylation, as well as A-to-inosine RNA editing) are closely related to the tumorigenesis and progression of many cancer types, including skin cutaneous melanoma (SKCM). However, the potential roles of the crosstalk between these RNA modification “writers” in the tumor microenvironment (TME) remain unclear. The RNA modification patterns were identified using an unsupervised clustering method. Subsequently, based on differentially expressed genes responsible for the aforementioned RNA modification patterns, an RNA modification “writer” scoring model (W_Score) was constructed to quantify the RNA modification-associated subtypes in individual patients. Moreover, a correlation analysis for W_Score and the TME characteristics, clinical features, molecular subtypes, drug sensitivities, immune responses, and prognosis was performed. We identified three RNA modification patterns, corresponding to distinct tumor immune microenvironment characteristics and survival outcomes. Based on the W_Score score, which was extracted from the RNA modification-related signature genes, patients with SKCM were divided into high- and low-W_Score groups. The low-W_Score group was characterized by better survival outcomes and strengthened immunocyte infiltration. Further analysis showed that the low-W_Score group was positively associated with higher tumor mutation burden and PD-L1 expression. Of note, two immunotherapy cohorts demonstrated that patients with low W_Score exhibited long-term clinical benefits and an enhanced immune response. This study is the first to systematically analyze four types of A-related RNA modifications in SKCM, revealing that these “writers” essentially contribute to TME complexity and diversity. We quantitatively evaluated the RNA modification patterns in individual tumors, which could aid in developing personalized immunotherapy strategies for patients.

LINC00885 promotes cervical cancer progression through sponging miR-3150b-3p and upregulating BAZ2A

Background

Cervical cancer (CC) is one of the most common malignancies affecting female worldwide. Long non-coding RNAs (lncRNAs) are increasingly indicated as crucial participants and promising therapeutic targets in human cancers. The main objective of this study was to explore the functions and mechanism of LINC00885 in CC.

Methods

RT-qPCR and western blot were used to detect RNA and protein levels. Functional and mechanism assays were respectively done for the analysis of cell behaviors and molecular interplays.

Results

Long intergenic non-coding RNA 885 (LINC00885) was discovered to be upregulated in CC tissues and cell lines through bioinformatics analysis and RT-qPCR. Overexpression of LINC00885 promoted proliferation and inhibited apoptosis, whereas its silence exerted opposite effects. The cytoplasmic localization of LINC00885 was ascertained and furthermore, LINC00885 competitively bound with miR-3150b-3p to upregulate BAZ2A expression in CC cells. Rescue assays confirmed that LINC00885 regulated CC proliferation and apoptosis through miR-3150b-3p/BAZ2A axis. Finally, we confirmed that LINC00885 aggravated tumor growth through animal experiments.

Conclusions

LINC00885 exerted oncogenic function in CC via regulating miR-3150b-3p/BAZ2A axis. These findings suggested LINC00885 might serve as a potential promising therapeutic target for CC patients.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13062-021-00314-6.

Identification and validation of a signature involving voltage-gated chloride ion channel genes for prediction of prostate cancer recurrence

Voltage-gated chloride ion channels (CLCs) are transmembrane proteins that maintain chloride ion homeostasis in various cells. Accumulating studies indicated CLCs were related to cell growth, proliferation, and cell cycle. Nevertheless, the role of CLCs in prostate cancer (PCa) has not been systematically profiled. The purpose of this study was to investigate the expression profiles and biofunctions of CLCs genes, and construct a novel risk signature to predict biochemical recurrence (BCR) of PCa patients. We identified five differentially expressed CLCs genes in our cohort and then constructed a signature composed of CLCN2 and CLCN6 through Lasso-Cox regression analysis in the training cohort from the Cancer Genome Atlas (TCGA). The testing and entire cohorts from TCGA and the GSE21034 from the Gene Expression Omnibus (GEO) were used as internal and independent external validation datasets. This signature could divide PCa patients into the high and low risk groups with different prognoses, was apparently correlated with clinical features, and was an independent excellent prognostic indicator. Enrichment analysis indicated our signature was primarily concentrated in cellular process and metabolic process. The expression patterns of CLCN2 and CLCN6 were detected in our own cohort based immunohistochemistry staining, and we found CLCN2 and CLCN6 were highly expressed in PCa tissues compared with benign tissues and positively associated with higher Gleason score and shorter BCR-free time. Functional experiments revealed that CLCN2 and CLCN6 downregulation inhibited cell proliferation, colony formation, invasion, and migration, but prolonged cell cycle and promoted apoptosis. Furthermore, Seahorse assay showed that silencing CLCN2 or CLCN6 exerted potential inhibitory effects on energy metabolism in PCa. Collectively, our signature could provide a novel and robust strategy for the prognostic evaluation and improve treatment decision making for PCa patients.

Long Non-coding RNAs Gabarapl2 and Chrnb2 Positively Regulate Inflammatory Signaling in a Mouse Model of Dry Eye

Purpose: To elucidate the expression profile and the potential role of long non-coding ribonucleic acids (RNAs) (lncRNAs) in a dry eye disease (DED) model.

Methods: A DED model was established in C57BL/6J mice with 0.2% benzalkonium chloride (BAC) twice a day for 14 days. The differentially expressed lncRNAs were detected by RNA-seq technology (Gene Expression Omnibus, GEO GSE186450) and the aberrantly expressed lncRNAs were further verified by RT-qPCR. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were conducted to predicate the related candidate genes and potential pathological pathways. Cells from a human corneal epithelial cell line (HCECs) were cultured under hyperosmolarity. The regulation of inflammatory factors by silencing potential targeted lncRNAs was verified in vitro in HCECs.

Results: In our study, a significant increase in corneal fluorescence staining and a reduction in tear production were observed in DED mice at all follow-ups compared with the controls, and the differences were increasing over time. In total, 2,649 upregulated and 704 downregulated lncRNAs were identified in DED mice. We selected six aberrantly expressed and most abundant lncRNAs and performed RT-qPCR using the samples for RNA-seq. Chrnb2, Gabarapl2, and Usp31 were thereby confirmed as the most significantly altered lncRNAs. Pathway analysis revealed that the neuroactive ligand–receptor interaction signaling pathway was the most enriched, followed by the calcium signaling pathway and cytokine–cytokine receptor interaction. Following treatment of Gabarapl2 siRNA and Chrnb2 siRNA, tumor necrosis factor-α (TNF-α), interleukin (IL)-1β, and IL-6 were significantly downregulated in the HCECs.

Conclusion: Our study suggests that Chrnb2 and Gabarapl2 may be involved in the inflammation response by regulating TNF-α, IL-1β, and IL-6 in DED. These candidate lncRNAs may be both potential biomarkers and therapeutic targets for DED.

Genome-Wide Analysis of H3K27me3 in Porcine Embryonic Muscle Development

The trimethylation of histone H3 lysine 27 (H3K27me3) is one of the most important chromatin modifications, which is generally presented as a repressive mark in various biological processes. However, the dynamic and global-scale distribution of H3K27me3 during porcine embryonic muscle development remains unclear. Here, our study provided a comprehensive genome-wide view of H3K27me3 and analyzed the matching transcriptome in the skeletal muscles on days 33, 65, and 90 post-coitus from Duroc fetuses. Transcriptome analysis identified 4,124 differentially expressed genes (DEGs) and revealed the key transcriptional properties in three stages. We found that the global H3K27me3 levels continually increased during embryonic development, and the H3K27me3 level was negatively correlated with gene expression. The loss of H3K27me3 in the promoter was associated with the transcriptional activation of 856 DEGs in various processes, including skeletal muscle development, calcium signaling, and multiple metabolic pathways. We also identified for the first time that H3K27me3 could enrich in the promoter of genes, such as DES, MYL1, TNNC1, and KLF5, to negatively regulate gene expression in porcine satellite cells (PSCs). The loss of H3K27me3 could promote muscle cell differentiation. Taken together, this study provided the first genome-wide landscape of H3K27me3 in porcine embryonic muscle development. It revealed the complex and broad function of H3K27me3 in the regulation of embryonic muscle development from skeletal muscle morphogenesis to myofiber maturation.

The emerging role of ISWI chromatin remodeling complexes in cancer

Disordered chromatin remodeling regulation has emerged as an essential driving factor for cancers. Imitation switch (ISWI) family are evolutionarily conserved ATP-dependent chromatin remodeling complexes, which are essential for cellular survival and function through multiple genetic and epigenetic mechanisms. Omics sequencing and a growing number of basic and clinical studies found that ISWI family members displayed widespread gene expression and genetic status abnormalities in human cancer. Their aberrant expression is closely linked to patient outcome and drug response. Functional or componential alteration in ISWI-containing complexes is critical for tumor initiation and development. Furthermore, ISWI-non-coding RNA regulatory networks and some non-coding RNAs derived from exons of ISWI member genes play important roles in tumor progression. Therefore, unveiling the transcriptional regulation mechanism underlying ISWI family sparked a booming interest in finding ISWI-based therapies in cancer. This review aims at describing the current state-of-the-art in the role of ISWI subunits and complexes in tumorigenesis, tumor progression, immunity and drug response, and presenting deep insight into the physiological and pathological implications of the ISWI transcription machinery in cancers.

Structural basis of the TAM domain of BAZ2A in binding to DNA or RNA independent of methylation status

Bromodomain adjacent to zinc finger domain protein 2A (BAZ2A) (also called transcription termination factor-1 interacting protein 5), a key component of the nucleolar remodeling complex, recruits the nucleolar remodeling complex to ribosomal RNA genes, leading to their transcriptional repression. In addition to its tandem plant homeodomain-bromodomain that is involved in binding to acetylated histone H4, BAZ2A also contains a methyl-CpG-binding domain (MBD)-like Tip5/ARBP/MBD (TAM) domain that shares sequence homology with the MBD. In contrast with the methyl-CpG-binding ability of the canonical MBD, the BAZ2A TAM domain has been shown to bind to promoter-associated RNAs of ribosomal RNA genes and promoter DNAs of other genes independent of DNA methylation. Nevertheless, how the TAM domain binds to RNA/DNA mechanistically remains elusive. Here, we characterized the DNA-/RNA-binding basis of the BAZ2A TAM domain by EMSAs, isothermal titration calorimetry binding assays, mutagenesis analysis, and X-ray crystallography. Our results showed that the TAM domain of BAZ2A selectively binds to dsDNA and dsRNA and that it binds to the backbone of dsDNA in a sequence nonspecific manner, which is distinct from the base-specific binding of the canonical MBD. Thus, our results explain why the TAM domain of BAZ2A does not specifically bind to mCG or TG dsDNA like the canonical MBD and also provide insights for further biological study of BAZ2A acting as a transcription factor in the future.

Non-cytochrome P450 enzymes involved in the oxidative metabolism of xenobiotics: Focus on the regulation of gene expression and enzyme activity

Oxidative metabolism is one of the major biotransformation reactions that regulates the exposure of xenobiotics and their metabolites in the circulatory system and local tissues and organs, and influences their efficacy and toxicity. Although cytochrome (CY)P450s play critical roles in the oxidative reaction, extensive CYP450-independent oxidative metabolism also occurs in some xenobiotics, such as aldehyde oxidase, xanthine oxidoreductase, flavin-containing monooxygenase, monoamine oxidase, alcohol dehydrogenase, or aldehyde dehydrogenase-dependent oxidative metabolism. Drugs form a large portion of xenobiotics and are the primary target of this review. The common reaction mechanisms and roles of non-CYP450 enzymes in metabolism, factors affecting the expression and activity of non-CYP450 enzymes in terms of inhibition, induction, regulation, and species differences in pharmaceutical research and development have been summarized. These non-CYP450 enzymes are detoxifying enzymes, although sometimes they mediate severe toxicity. Synthetic or natural chemicals serve as inhibitors for these non-CYP450 enzymes. However, pharmacokinetic-based drug interactions through these inhibitors have rarely been reported in vivo. Although multiple mechanisms participate in the basal expression and regulation of non-CYP450 enzymes, only a limited number of inducers upregulate their expression. Therefore, these enzymes are considered non-inducible or less inducible. Overall, this review focuses on the potential xenobiotic factors that contribute to variations in gene expression levels and the activities of non-CYP450 enzymes.

Alpha-ketoglutarate decorated iron oxide-gold core-shell nanoparticles for active mitochondrial targeting and radiosensitization enhancement in hepatocellular carcinoma

The ability of some tumours to impart radioresistance serves as a barrier in the cancer therapeutics. Mitochondrial metabolism significantly persuades this cancer cell survival, incursion and plays a crucial role in conferring radioresistance. It would be of great importance to target the active mitochondria to overcome this resistance and achieve tumoricidal efficacy. The current report investigates the improved radiosensitization effect (under Gamma irradiation) in hepatocellular carcinoma through active mitochondrial targeting of alpha-ketoglutarate decorated iron oxide-gold core-shell nanoparticles (GNP). The loading of a chemotherapeutic drug N-(4-hydroxyphenyl)retinamide in GNP allows adjuvant chemotherapy, which further sensitizes cancerous cells for radiotherapy. The GNP shows a drug loading efficiency of 8.5 wt% with a sustained drug release kinetics. The X-Ray diffraction (XRD) pattern and High-Resolution Transmission Electron microscopy (HRTEM) indicates the synthesis of core iron oxide nanoparticles with indications of a thin layer of gold shell on the surface with 1:7 ratios of Fe: Au. The GNP application significantly reduced per cent cell viability in Hepatocellular carcinoma cells through improved radiosensitization at 5 Gy gamma radiation dose. The molecular mechanism revealed a sharp increment in reactive oxygen species (ROS) generation and DNA fragmentation. The mitochondrial targeting probes confirm the presence of GNP in the mitochondria, which could be the possible reason for such improved cellular damage. In addition to the active mitochondrial targeting, the currently fabricated nanoparticles work as a potent Magnetic Resonance Imaging (MRI)/Computed Tomography (CT) contrast agent. This multifunctional therapeutic potential makes GNP as one of the most promising theragnostic molecules in cancer therapeutics.

Identification of a Costimulatory Molecule-Related Signature for Predicting Prognostic Risk in Prostate Cancer

Costimulatory molecules have been proven to enhance antitumor immune responses, but their roles in prostate cancer (PCa) remain unexplored. In this study, we aimed to explore the gene expression profiles of costimulatory molecule genes in PCa and construct a prognostic signature to improve treatment decision making and clinical outcomes. Five prognosis-related costimulatory molecule genes (RELT, TNFRSF25, EDA2R, TNFSF18, and TNFSF10) were identified, and a prognostic signature was constructed based on these five genes. This signature was an independent prognostic factor according to multivariate Cox regression analysis; it could stratify PCa patients into two subgroups with different prognoses and was highly associated with clinical features. The prognostic significance of the signature was well validated in four different independent external datasets. Moreover, patients identified as high risk based on our prognostic signature exhibited a high mutation frequency, a high level of immune cell infiltration and an immunosuppressive microenvironment. Therefore, our signature could provide clinicians with prognosis predictions and help guide treatment for PCa patients.

BAZ2A-mediated repression via H3K14ac-marked enhancers promotes prostate cancer stem cells

Prostate cancer (PCa) is one of the most prevalent cancers in men. Cancer stem cells are thought to be associated with PCa relapse. Here, we show that BAZ2A is required for PCa cells with a cancer stem‐like state. BAZ2A genomic occupancy in PCa cells coincides with H3K14ac‐enriched chromatin regions. This association is mediated by BAZ2A‐bromodomain (BAZ2A‐BRD) that specifically binds H3K14ac. BAZ2A associates with inactive enhancers marked by H3K14ac and repressing transcription of genes frequently silenced in aggressive and poorly differentiated PCa. BAZ2A‐mediated repression is also linked to EP300 that acetylates H3K14ac. BAZ2A‐BRD mutations or treatment with inhibitors abrogating BAZ2A‐BRD/H3K14ac interaction impair PCa stem cells. Furthermore, pharmacological inactivation of BAZ2A‐BRD impairs Pten‐loss oncogenic transformation of prostate organoids. Our findings indicate a role of BAZ2A‐BRD in PCa stem cell features and suggest potential epigenetic‐reader therapeutic strategies to target BAZ2A in aggressive PCa.

Functional Roles of Bromodomain Proteins in Cancer

Histone acetylation is generally associated with an open chromatin configuration that facilitates many cellular processes including gene transcription, DNA repair, and DNA replication. Aberrant levels of histone lysine acetylation are associated with the development of cancer. Bromodomains represent a family of structurally well-characterized effector domains that recognize acetylated lysines in chromatin. As part of their fundamental reader activity, bromodomain-containing proteins play versatile roles in epigenetic regulation, and additional functional modules are often present in the same protein, or through the assembly of larger enzymatic complexes. Dysregulated gene expression, chromosomal translocations, and/or mutations in bromodomain-containing proteins have been correlated with poor patient outcomes in cancer. Thus, bromodomains have emerged as a highly tractable class of epigenetic targets due to their well-defined structural domains, and the increasing ease of designing or screening for molecules that modulate the reading process. Recent developments in pharmacological agents that target specific bromodomains has helped to understand the diverse mechanisms that bromodomains play with their interaction partners in a variety of chromatin processes, and provide the promise of applying bromodomain inhibitors into the clinical field of cancer treatment. In this review, we explore the expression and protein interactome profiles of bromodomain-containing proteins and discuss them in terms of functional groups. Furthermore, we highlight our current understanding of the roles of bromodomain-containing proteins in cancer, as well as emerging strategies to specifically target bromodomains, including combination therapies using bromodomain inhibitors alongside traditional therapeutic approaches designed to re-program tumorigenesis and metastasis.

Identification and Validation of Immune-Related Methylation Clusters for Predicting Immune Activity and Prognosis in Breast Cancer

The role of DNA methylation of breast cancer-infiltrating immune cells has not been fully explored. We conducted a cohort-based retrospective study analyzing the genome-wide immune-related DNA methylation of 1057 breast cancer patients from the TCGA cohort and GSE72308 cohort. Based on patients' overall survival (OS), a prognostic risk score system using 18 immune-related methylation genes (IRMGs) was established and further validated in an independent cohort. Kaplan-Meier analysis showed a clear separation of OS between the low- and high-risk groups. Patients in the low-risk group had a higher immune score and stromal score compared with the high-risk group. Moreover, the characteristics based on 18-IRMGs signature were related to the tumor immune microenvironment and affected the abundance of tumor-infiltrating immune cells. Consistently, the 18-IRMGs signatures showed similar influences on immune modulation and survival in another external validation cohort (GSE72308). In conclusion, the proposed 18-IRMGs signature could be a potential marker for breast cancer prognostication.

Binding Selectivity of Inhibitors toward Bromodomains BAZ2A and BAZ2B Uncovered by Multiple Short Molecular Dynamics Simulations and MM-GBSA Calculations

Two Bromodomain-Containing proteins BAZ2A and BAZ2B are responsible for remodeling chromatin and regulating noncoding RNAs. As for our current studies, integration of multiple short molecular dynamics simulations (MSMDSs) with molecular mechanics generalized Born surface area (MM-GBSA) method is adopted for insights into binding selectivity of three small molecules D8Q, D9T and UO1 to BAZ2A against BAZ2B. The calculations of MM-GBSA unveil that selectivity of inhibitors toward BAZ2A and BAZ2B highly depends on the enthalpy changes and the details uncover that D8Q has better selectivity toward BAZ2A than BAZ2B, D9T more favorably bind to BAZ2B than BAZ2A, and UO1 does not show obvious selectivity toward these two proteins. The analysis of interaction network between residues and inhibitors indicates that seven residues are mainly responsible for the selectivity of D8Q, six residues for D9T and four residues provide significant contributions to associations of UO1 with two proteins. Moreover the analysis of interaction network not only reveals warm spots of inhibitor bindings to BAZ2A and BAZ2B but also unveils that common residue pairs, including (W1816, W1887), (P1817, P1888), (F1818, F1889), (V1822, V1893), (N1823, N1894),(L1826, L1897), (V1827, V1898), (F1872, F1943), (N1873, N1944) and (V1879, I1950) belonging to (BAZ2A, BAZ2B), induce mainly binding differences of inhibitors to BAZ2A and BAZ2B. Hence, insights from our current studies offer useful dynamics information relating with conformational alterations and structure-affinity relationship at atomistic levels for novel therapeutic strategies toward BAZ2A and BAZ2B.

Epigenetic reprogramming during prostate cancer progression: A perspective from development

Conrad Waddington's theory of epigenetic landscape epitomize the process of cell fate and cellular decision-making during development. Wherein the epigenetic code maintains patterns of gene expression in pluripotent and differentiated cellular states during embryonic development and differentiation. Over the years disruption or reprogramming of the epigenetic landscape has been extensively studied in the course of cancer progression. Cellular dedifferentiation being a key hallmark of cancer allow us to take cues from the biological processes involved during development. Here, we discuss the role of epigenetic landscape and its modifiers in cell-fate determination, differentiation and prostate cancer progression. Lately, the emergence of RNA-modifications has also furthered our understanding of epigenetics in cancer. The overview of the epigenetic code regulating androgen signalling, and progression to aggressive neuroendocrine stage of PCa reinforces its gene regulatory functions during the development of prostate gland as well as cancer progression. Additionally, we also highlight the clinical implications of cancer cell epigenome, and discuss the recent advancements in the therapeutic strategies targeting the advanced stage disease.

Identification of Key MicroRNAs and Mechanisms in Prostate Cancer Evolution Based on Biomarker Prioritization Model and Carcinogenic Survey

Background: Prostate cancer (PCa) is occurred with increasing incidence and heterogeneous pathogenesis. Although clinical strategies are accumulated for PCa prevention, there is still a lack of sensitive biomarkers for the holistic management in PCa occurrence and progression. Based on systems biology and artificial intelligence, translational informatics provides new perspectives for PCa biomarker prioritization and carcinogenic survey.

Methods: In this study, gene expression and miRNA-mRNA association data were integrated to construct conditional networks specific to PCa occurrence and progression, respectively. Based on network modeling, hub miRNAs with significantly strong single-line regulatory power were topologically identified and those shared by the condition-specific network systems were chosen as candidate biomarkers for computational validation and functional enrichment analysis.

Results: Nine miRNAs, i.e., hsa-miR-1-3p, hsa-miR-125b-5p, hsa-miR-145-5p, hsa-miR-182-5p, hsa-miR-198, hsa-miR-22-3p, hsa-miR-24-3p, hsa-miR-34a-5p, and hsa-miR-499a-5p, were prioritized as key players for PCa management. Most of these miRNAs achieved high AUC values (AUC > 0.70) in differentiating different prostate samples. Among them, seven of the miRNAs have been previously reported as PCa biomarkers, which indicated the performance of the proposed model. The remaining hsa-miR-22-3p and hsa-miR-499a-5p could serve as novel candidates for PCa predicting and monitoring. In particular, key miRNA-mRNA regulations were extracted for pathogenetic understanding. Here hsa-miR-145-5p was selected as the case and hsa-miR-145-5p/NDRG2/AR and hsa-miR-145-5p/KLF5/AR axis were found to be putative mechanisms during PCa evolution. In addition, Wnt signaling, prostate cancer, microRNAs in cancer etc. were significantly enriched by the identified miRNAs-mRNAs, demonstrating the functional role of the identified miRNAs in PCa genesis.

Conclusion: Biomarker miRNAs together with the associated miRNA-mRNA relations were computationally identified and analyzed for PCa management and carcinogenic deciphering. Further experimental and clinical validations using low-throughput techniques and human samples are expected for future translational studies.

Quantitative Proteomic Analysis of Porcine Intestinal Epithelial Cells Infected with Porcine Deltacoronavirus Using iTRAQ-Coupled LC-MS/MS

Porcine deltacoronavirus (PDCoV) is an emergent enteropathogenic coronavirus associated with swine diarrhea. Porcine small intestinal epithelial cells (IPEC) are the primary target cells of PDCoV infection in vivo. Here, isobaric tags for relative and absolute quantification (iTRAQ) labeling coupled to liquid chromatography-tandem mass spectrometry (LC-MS/MS) was used to quantitatively identify differentially expressed proteins (DEPs) in PDCoV-infected IPEC-J2 cells. A total of 78 DEPs, including 23 upregulated and 55 downregulated proteins, were identified at 24 h postinfection. The data are available via ProteomeXchange with identifier PXD019975. To ensure reliability of the proteomics data, two randomly selected DEPs, the downregulated anaphase-promoting complex subunit 7 (ANAPC7) and upregulated interferon-induced protein with tetratricopeptide repeats 1 (IFIT1), were verified by real-time PCR and Western blot, and the results of which indicate that the proteomics data were reliable and valid. Bioinformatics analyses, including GO, COG, KEGG, and STRING, further demonstrated that a majority of the DEPs are involved in numerous crucial biological processes and signaling pathways, such as immune system, digestive system, signal transduction, RIG-I-like receptor, mTOR, PI3K-AKT, autophagy, and cell cycle signaling pathways. Altogether, this is the first study on proteomes of PDCoV-infected host cells, which shall provide valuable clues for further investigation of PDCoV pathogenesis.

Integration of anatomy ontology data with protein-protein interaction networks improves the candidate gene prediction accuracy for anatomical entities

Background

Identification of genes responsible for anatomical entities is a major requirement in many fields including developmental biology, medicine, and agriculture. Current wet lab techniques used for this purpose, such as gene knockout, are high in resource and time consumption. Protein–protein interaction (PPI) networks are frequently used to predict disease genes for humans and gene candidates for molecular functions, but they are rarely used to predict genes for anatomical entities. Moreover, PPI networks suffer from network quality issues, which can be a limitation for their usage in predicting candidate genes. Therefore, we developed an integrative framework to improve the candidate gene prediction accuracy for anatomical entities by combining existing experimental knowledge about gene-anatomical entity relationships with PPI networks using anatomy ontology annotations. We hypothesized that this integration improves the quality of the PPI networks by reducing the number of false positive and false negative interactions and is better optimized to predict candidate genes for anatomical entities. We used existing Uberon anatomical entity annotations for zebrafish and mouse genes to construct gene networks by calculating semantic similarity between the genes. These anatomy-based gene networks were semantic networks, as they were constructed based on the anatomy ontology annotations that were obtained from the experimental data in the literature. We integrated these anatomy-based gene networks with mouse and zebrafish PPI networks retrieved from the STRING database and compared the performance of their network-based candidate gene predictions.

Results

According to evaluations of candidate gene prediction performance tested under four different semantic similarity calculation methods (Lin, Resnik, Schlicker, and Wang), the integrated networks, which were semantically improved PPI networks, showed better performances by having higher area under the curve values for receiver operating characteristic and precision-recall curves than PPI networks for both zebrafish and mouse.

Conclusion

Integration of existing experimental knowledge about gene-anatomical entity relationships with PPI networks via anatomy ontology improved the candidate gene prediction accuracy and optimized them for predicting candidate genes for anatomical entities.

Targeting epigenetic protein-protein interactions with small-molecule inhibitors

Epigenetic protein-protein interactions (PPIs) play essential roles in regulating gene expression, and their dysregulations have been implicated in many diseases. These PPIs are comprised of reader domains recognizing post-translational modifications on histone proteins, and of scaffolding proteins that maintain integrities of epigenetic complexes. Targeting PPIs have become focuses for development of small-molecule inhibitors and anticancer therapeutics. Here we summarize efforts to develop small-molecule inhibitors targeting common epigenetic PPI domains. Potent small molecules have been reported for many domains, yet small domains that recognize methylated lysine side chains on histones are challenging in inhibitor development. We posit that the development of potent inhibitors for difficult-to-prosecute epigenetic PPIs may be achieved by interdisciplinary approaches and extensive explorations of chemical space.