EpiFactors: a comprehensive database of human epigenetic factors and complexes

Key chromatin regulator-related genes associated with the risk of coronary artery disease regulate the expression of HCFC1, RNF8, TNP1 and SET

Chromatin regulators are indispensable upstream epigenetic regulators.The emergence and progression of atherosclerosis has been demonstrated to be influenced by smooth muscle-related chromatin regulators, such as ZEB2 and MAFF. However, specific chromatin regulators and their possible roles have not been clarified. Information was gathered from 51 patients diagnosed with coronary artery disease (CAD) and 50 individuals in good health from the GEO database. 440 genes were identified as having differential expression across the two datasets, and these genes were linked to cellular reactions. Enrichment of pathways related to histone modification and transcriptional regulatory factors was observed in GO and KEGG analyses. Four machine learning models (RF, SVM, GLM, and XGB) were developed using the expression profiles of 440 chromatin-associated genes in the CAD cohort to pinpoint genes with significant diagnostic potential. After evaluating residuals, root mean square errors, receiver operating characteristic curves, and immune-infiltration, four key genes (HCFC1, RNF8, TNP1, and SET) were identified. Gene expression in different blood vessel levels in atherosclerotic plaques in a mouse model of coronary artery disease showed significant variations. The gene expression levels in macrophages aligned with clinical data from the GEO database as expected. This discovery is crucial for future analysis and the prediction of drug and miRNA targets. In conclusion, we found that the four hub genes are important in the mechanism of CAD. These findings provide new ideas for the study of potential epigenetic predictive markers and therapeutic targets to be used in determining a treatment strategy for CAD.

Annelid adult cell type diversity and their pluripotent cellular origins

Many annelids can regenerate missing body parts or reproduce asexually, generating all cell types in adult stages. However, the putative adult stem cell populations involved in these processes, and the diversity of cell types generated by them, are still unknown. To address this, we recover 75,218 single cell transcriptomes of the highly regenerative and asexually-reproducing annelid Pristina leidyi. Our results uncover a rich cell type diversity including annelid specific types as well as novel types. Moreover, we characterise transcription factors and gene networks that are expressed specifically in these populations. Finally, we uncover a broadly abundant cluster of putative stem cells with a pluripotent signature. This population expresses well-known stem cell markers such as vasa, piwi and nanos homologues, but also shows heterogeneous expression of differentiated cell markers and their transcription factors. We find conserved expression of pluripotency regulators, including multiple chromatin remodelling and epigenetic factors, in piwi+ cells. Finally, lineage reconstruction analyses reveal computational differentiation trajectories from piwi+ cells to diverse adult types. Our data reveal the cell type diversity of adult annelids by single cell transcriptomics and suggest that a piwi+ cell population with a pluripotent stem cell signature is associated with adult cell type differentiation.

Polyploidy Promotes Hypertranscription, Apoptosis Resistance, and Ciliogenesis in Cancer Cells and Mesenchymal Stem Cells of Various Origins: Comparative Transcriptome In Silico Study

Mesenchymal stem cells (MSC) attract an increasing amount of attention due to their unique therapeutic properties. Yet, MSC can undergo undesirable genetic and epigenetic changes during their propagation in vitro. In this study, we investigated whether polyploidy can compromise MSC oncological safety and therapeutic properties. For this purpose, we compared the impact of polyploidy on the transcriptome of cancer cells and MSC of various origins (bone marrow, placenta, and heart). First, we identified genes that are consistently ploidy-induced or ploidy-repressed through all comparisons. Then, we selected the master regulators using the protein interaction enrichment analysis (PIEA). The obtained ploidy-related gene signatures were verified using the data gained from polyploid and diploid populations of early cardiomyocytes (CARD) originating from iPSC. The multistep bioinformatic analysis applied to the cancer cells, MSC, and CARD indicated that polyploidy plays a pivotal role in driving the cell into hypertranscription. It was evident from the upregulation of gene modules implicated in housekeeping functions, stemness, unicellularity, DNA repair, and chromatin opening by means of histone acetylation operating via DNA damage associated with the NUA4/TIP60 complex. These features were complemented by the activation of the pathways implicated in centrosome maintenance and ciliogenesis and by the impairment of the pathways related to apoptosis, the circadian clock, and immunity. Overall, our findings suggest that, although polyploidy does not induce oncologic transformation of MSC, it might compromise their therapeutic properties because of global epigenetic changes and alterations in fundamental biological processes. The obtained results can contribute to the development and implementation of approaches enhancing the therapeutic properties of MSC by removing polyploid cells from the cell population.

Rapid profiling of transcription factor-cofactor interaction networks reveals principles of epigenetic regulation

Transcription factor (TF)-cofactor (COF) interactions define dynamic, cell-specific networks that govern gene expression; however, these networks are understudied due to a lack of methods for high-throughput profiling of DNA-bound TF-COF complexes. Here we describe the Cofactor Recruitment (CoRec) method for rapid profiling of cell-specific TF-COF complexes. We define a lysine acetyltransferase (KAT)-TF network in resting and stimulated T cells. We find promiscuous recruitment of KATs for many TFs and that 35% of KAT-TF interactions are condition specific. KAT-TF interactions identify NF-κB as a primary regulator of acutely induced H3K27ac. Finally, we find that heterotypic clustering of CBP/P300-recruiting TFs is a strong predictor of total promoter H3K27ac. Our data supports clustering of TF sites that broadly recruit KATs as a mechanism for widespread co-occurring histone acetylation marks. CoRec can be readily applied to different cell systems and provides a powerful approach to define TF-COF networks impacting chromatin state and gene regulation.

The omics era: a nexus of untapped potential for Mendelian chromatinopathies

The OMICs cascade describes the hierarchical flow of information through biological systems. The epigenome sits at the apex of the cascade, thereby regulating the RNA and protein expression of the human genome and governs cellular identity and function. Genes that regulate the epigenome, termed epigenes, orchestrate complex biological signaling programs that drive human development. The broad expression patterns of epigenes during human development mean that pathogenic germline mutations in epigenes can lead to clinically significant multi-system malformations, developmental delay, intellectual disabilities, and stem cell dysfunction. In this review, we refer to germline developmental disorders caused by epigene mutation as "chromatinopathies". We curated the largest number of human chromatinopathies to date and our expanded approach more than doubled the number of established chromatinopathies to 179 disorders caused by 148 epigenes. Our study revealed that 20.6% (148/720) of epigenes cause at least one chromatinopathy. In this review, we highlight key examples in which OMICs approaches have been applied to chromatinopathy patient biospecimens to identify underlying disease pathogenesis. The rapidly evolving OMICs technologies that couple molecular biology with high-throughput sequencing or proteomics allow us to dissect out the causal mechanisms driving temporal-, cellular-, and tissue-specific expression. Using the full repertoire of data generated by the OMICs cascade to study chromatinopathies will provide invaluable insight into the developmental impact of these epigenes and point toward future precision targets for these rare disorders.

Identification and validation of the clinical prediction model and biomarkers based on chromatin regulators in colon cancer by integrated analysis of bulk- and single-cell RNA sequencing data

Background

Chromatin regulators (CRs) are implicated in the development of cancer, but a comprehensive investigation of their role in colon adenocarcinoma (COAD) is inadequate. The purpose of this study is to find CRs that can provide recommendations for clinical diagnosis and treatment, and to explore the reasons why they serve as critical CRs.

Methods

We obtained data from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. Weighted Gene Co-Expression Network Analysis (WGCNA) screened tumor-associated CRs. LASSO-Cox regression was used to construct the model and to screen key CRs together with support vector machine (SVM), the univariate Cox regression. We used single-cell data to explore the expression of CRs in cells and their communication. Immune infiltration, immune checkpoints, mutation, methylation, and drug sensitivity analyses were performed. Gene expression was verified by quantitative real-time reverse transcription-polymerase chain reaction (qRT-PCR). Pan-cancer analysis was used to explore the importance of hub CRs.

Results

We finally obtained 32 tumor-associated CRs. The prognostic model was constructed based on RCOR2, PPARGC1A, PKM, RAC3, PHF19, MYBBP1A, ORC1, and EYA2 by the LASSO-Cox regression. Single-cell data revealed that the model was immune-related. Combined with immune infiltration analysis, immune checkpoint analysis, and tumor immune dysfunction and exclusion (TIDE) analysis, the low-score risk group had more immune cell infiltration and better immune response. Mutation and methylation analysis showed that multiple CRs may be mutated and methylated in colon cancer. Drug sensitivity analysis revealed that the low-risk group may be more sensitive to several drugs and PKM was associated with multiple drugs. Combined with machine learning, PKM is perhaps the most critical gene in CRs. Pan-cancer analysis showed that PKM plays a role in the prognosis of cancers.

Conclusions

We developed a prognostic model for COAD based on CRs. Increased expression of the core gene PKM is linked with a poor prognosis in several malignancies.

MLL1 regulates cytokine-driven cell migration and metastasis

Cell migration is a critical contributor to metastasis. Cytokine production and its role in cancer cell migration have been traditionally associated with immune cells. We find that the histone methyltransferase Mixed-Lineage Leukemia 1 (MLL1) controls 3D cell migration via cytokines, IL-6, IL-8, and TGF-β1, secreted by the cancer cells themselves. MLL1, with its scaffold protein Menin, controls actin filament assembly via the IL-6/8/pSTAT3/Arp3 axis and myosin contractility via the TGF-β1/Gli2/ROCK1/2/pMLC2 axis, which together regulate dynamic protrusion generation and 3D cell migration. MLL1 also regulates cell proliferation via mitosis-based and cell cycle-related pathways. Mice bearing orthotopic MLL1-depleted tumors exhibit decreased lung metastatic burden and longer survival. MLL1 depletion leads to lower metastatic burden even when controlling for the difference in primary tumor growth rates. Combining MLL1-Menin inhibitor with paclitaxel abrogates tumor growth and metastasis, including preexistent metastasis. These results establish MLL1 as a potent regulator of cell migration and highlight the potential of targeting MLL1 in patients with metastatic disease.

STRN3 promotes tumour growth in hepatocellular carcinoma by inhibiting the hippo pathway

HCC is a globally high-incidence malignant tumour, and its pathogenesis is still unclear. Recently, STRN3 has been found to be elevated in various tumours, but its expression and biological functions in HCC have not been studied. In the study, clinical correlation analysis was performed on 371 liver cancer patients from TCGA database and liver cancer tissues and normal tissues from the GEO database. qRT-PCR and western blotting were used to detect relevant proteins in cells, and CCK8 and colony formation experiments were performed to analyse cell proliferation ability. Transwell and wound healing experiments were performed to detect cell invasion ability, and flow cytometry was used to detect cell apoptosis. Single-cell sequencing data and multiple immunofluorescence were analysed for the expression abundance and distribution of certain proteins. Immunohistochemistry was used to assess the expression of STRN3 in patients' tumour and adjacent non-cancerous tissues. The results indicated STRN3 was highly expressed in liver tumour tissues and was closely associated with poor prognosis. Knockdown of STRN3 could significantly inhibit cell proliferation and migration ability. At the same time, we found that STRN3 could inhibit the Hippo pathway and promote the entry of YAP protein into the nucleus. Our study first found that STRN3 could promote tumour growth by inhibiting the Hippo pathway. The study of STRN3 can promote the understanding and treatment of the occurrence and development of HCC.

Decoding chromatin states by proteomic profiling of nucleosome readers

DNA and histone modifications combine into characteristic patterns that demarcate functional regions of the genome1,2. While many 'readers' of individual modifications have been described3-5, how chromatin states comprising composite modification signatures, histone variants and internucleosomal linker DNA are interpreted is a major open question. Here we use a multidimensional proteomics strategy to systematically examine the interaction of around 2,000 nuclear proteins with over 80 modified dinucleosomes representing promoter, enhancer and heterochromatin states. By deconvoluting complex nucleosome-binding profiles into networks of co-regulated proteins and distinct nucleosomal features driving protein recruitment or exclusion, we show comprehensively how chromatin states are decoded by chromatin readers. We find highly distinctive binding responses to different features, many factors that recognize multiple features, and that nucleosomal modifications and linker DNA operate largely independently in regulating protein binding to chromatin. Our online resource, the Modification Atlas of Regulation by Chromatin States (MARCS), provides in-depth analysis tools to engage with our results and advance the discovery of fundamental principles of genome regulation by chromatin states.

A high throughput CRISPR perturbation screen identifies epigenetic regulators impacting primordial germ cell development

Certain environmental factors can impact fertility and reproductive parameters such as the number and quality of sperm and eggs. One possible mechanism is the perturbation of epigenetic landscapes in the germline. To explore this possibility, we conducted a CRISPRi screen of epigenetic-related genes to identify those that specifically perturb the differentiation of embryonic stem cells (ESCs) into primordial germ cell-like cells (PGCLCs), exploiting a highly scalable cytokine-free platform. Of the 701 genes screened, inhibition of 53 decreased the efficiency of PGCLC formation. NCOR2, a transcriptional repressor that acts via recruitment of Class I and Class IIa histone deacetylases (HDACs) to gene targets, was particularly potent in suppressing PGCLC differentiation. Consistent with evidence that histone deacetylation is crucial for germline differentiation, we found that the HDAC inhibitors (HDACi) valproic acid (VPA; an anti-convulsant) and sodium butyrate (SB; a widely-used dietary supplement) also suppressed ESC>PGCLC differentiation. Furthermore, exposure of developing mouse embryos to SB or VPA caused hypospermatogenesis. Transcriptome analyses of HDACi-treated, differentiating ESC>PGCLC cultures revealed suppression of germline-associated pathways and enhancement of somatic pathways. This work demonstrates the feasibility of conducting large-scale functional screens of genes, chemicals, or other agents that may impact germline development.

Construction and validation of a chromatin regulator-related gene signature for prognostic and therapeutic significance of clear cell renal cell carcinoma

Background

Epigenetic alterations driven by chromatin regulators (CRs) are well-recognized cancer hallmarks. Growing evidence suggests that the imbalance of CRs may lead to the occurrence of various diseases including tumors. However, the role and prognostic value of CRs in clear cell renal cell carcinoma (ccRCC) remain undefined.

Methods

Consensus clustering analysis was used to identify different subtypes. Univariate and multivariate Cox regression analysis were performed to identify prognosis-related CRs and constructed a risk model. Transcriptome sequencing was used to verify gene expression levels. Kaplan-Meier survival analysis was used to compare overall survival (OS) between high- and low-risk groups. The area under the curve (AUC) value of the receiver operating characteristic (ROC) curve was used to evaluate the performance of the model. The ESTIMATE algorithm and single-sample gene set enrichment analysis (ssGSEA) were executed to evaluate the immune characteristics of samples. Correlation analysis was used to assess the relationship between risk score and immune checkpoint genes, the relationship between expression levels of CRs and immune cell infiltration and drug therapeutic response. Finally, we also compared differences in drug sensitivity between low- and high-risk groups.

Results

We identified three CRs-related subtypes with different characteristics. A prognostic model was built with four CRs and can precisely predict the OS of patients in different risk groups. The model has good stability and applicability and was further verified in the internal and external dataset. The transcriptomic levels of the four CRs were also validated, and the risk score was an independent prognostic factor for ccRCC. Obvious differences in the immune microenvironment and the expression levels of immune checkpoints were observed in low- and high-risk group. Higher immune activity and immune cell infiltration were found in the high-risk group. Besides, the expression levels of CRs were associated with drug therapeutic response. Patients with high-risk score may be more sensitive to gemcitabine, vinblastine, paclitaxel, axitinib, sunitinib, and temsirolimus.

Conclusions

CRs were strongly associated with the occurrence and development of ccRCC. Targeting CRs may become a new therapeutic strategy for ccRCC. Besides, CRs-related gene signature can predict the prognosis and therapeutic significance of ccRCC, which provides an important reference for clinical decision-making.

A CIC-related-epigenetic factors-based model associated with prediction, the tumor microenvironment and drug sensitivity in osteosarcoma

Osteosarcoma is generally considered a cold tumor and is characterized by epigenetic alterations. Although tumor cells are surrounded by many immune cells such as macrophages, T cells may be suppressed, be inactivated, or not be presented due to various mechanisms, which usually results in poor prognosis and insensitivity to immunotherapy. Immunotherapy is considered a promising anti-cancer therapy in osteosarcoma but requires more research, but osteosarcoma does not currently respond well to this therapy. The cancer immunity cycle (CIC) is essential for anti-tumor immunity, and is epigenetically regulated. Therefore, it is possible to modulate the immune microenvironment of osteosarcoma by targeting epigenetic factors. In this study, we explored the correlation between epigenetic modulation and CIC in osteosarcoma through bioinformatic methods. Based on the RNA data from TARGET and GSE21257 cohorts, we identified epigenetic related subtypes by NMF clustering and constructed a clinical prognostic model by the LASSO algorithm. ESTIMATE, Cibersort, and xCell algorithms were applied to analyze the tumor microenvironment. Based on eight epigenetic biomarkers (SFMBT2, SP140, CBX5, HMGN2, SMARCA4, PSIP1, ACTR6, and CHD2), two subtypes were identified, and they are mainly distinguished by immune response and cell cycle regulation. After excluding ACTR6 by LASSO regression, the prognostic model was established and it exhibited good predictive efficacy. The risk score showed a strong correlation with the tumor microenvironment, drug sensitivity and many immune checkpoints. In summary, our study sheds a new light on the CIC-related epigenetic modulation mechanism of osteosarcoma and helps search for potential drugs for osteosarcoma treatment.

A modular dCas9-based recruitment platform for combinatorial epigenome editing

Targeted epigenome editing tools allow precise manipulation and investigation of genome modifications, however they often display high context dependency and variable efficacy between target genes and cell types. While systems that simultaneously recruit multiple distinct 'effector' chromatin regulators can improve efficacy, they generally lack control over effector composition and spatial organisation. To overcome this we have created a modular combinatorial epigenome editing platform, called SSSavi. This system is an interchangeable and reconfigurable docking platform fused to dCas9 that enables simultaneous recruitment of up to four different effectors, allowing precise control of effector composition and spatial ordering. We demonstrate the activity and specificity of the SSSavi system and, by testing it against existing multi-effector targeting systems, demonstrate its comparable efficacy. Furthermore, we demonstrate the importance of the spatial ordering of the recruited effectors for effective transcriptional regulation. Together, the SSSavi system enables exploration of combinatorial effector co-recruitment to enhance manipulation of chromatin contexts previously resistant to targeted editing.

Development and validation of epigenetic modification-related signals for the diagnosis and prognosis of colorectal cancer

BACKGROUND

Colorectal cancer (CRC) is one of the world's most common malignancies. Epigenetics is the study of heritable changes in characteristics beyond the DNA sequence. Epigenetic information is essential for maintaining specific expression patterns of genes and the normal development of individuals, and disorders of epigenetic modifications may alter the expression of oncogenes and tumor suppressor genes and affect the development of cancer. This study elucidates the relationship between epigenetics and the prognosis of CRC patients by developing a predictive model to explore the potential value of epigenetics in the treatment of CRC.

METHODS

Gene expression data of CRC patients' tumor tissue and controls were downloaded from GEO database. Combined with the 720 epigenetic-related genes (ERGs) downloaded from EpiFactors database, prognosis-related epigenetic genes were selected by univariate cox and LASSO analyses. The Kaplan-Meier and ROC curve were used to analyze the accuracy of the model. Data of 238 CRC samples with survival data downloaded from the GSE17538 were used for validation. Finally, the risk model is combined with the clinical characteristics of CRC patients to perform univariate and multivariate cox regression analysis to obtain independent risk factors and draw nomogram. Then we evaluated the accuracy of its prediction by calibration curves.

RESULTS

A total of 2906 differentially expressed genes (DEGs) were identified between CRC and control samples. After overlapping DEGs with 720 ERGs, 56 epigenetic-related DEGs (DEERGs) were identified. Combining univariate and LASSO regression analysis, the 8 epigenetic-related genes-based risk score model of CRC was established. The ROC curves and survival difference of high and low risk groups revealed the good performance of the risk score model based on prognostic biomarkers in both training and validation sets. A nomogram with good performance to predict the survival of CRC patients were established based on age, NM stage and risk score. The calibration curves showed that the prognostic model had good predictive performance.

CONCLUSION

In this study, an epigenetically relevant 8-gene signature was constructed that can effectively predict the prognosis of CRC patients and provide potential directions for targeted therapies for CRC.

Temporal and Spatial Signatures of Scylla paramamosain Transcriptome Reveal Mechanistic Insights into Endogenous Ovarian Maturation under Risk of Starvation

Variability in food availability leads to condition-dependent investments in reproduction. This study is aimed at understanding the metabolic response and regulatory mechanism of female Scylla paramamosain in response to starvation in a temporal- and tissue-specific manner. The mud crabs were starved for 7 (control), 14, 28, and 40 days for histological and biochemical analysis in the hepatopancreas, ovary, and serum, as well as for RNA sequencing on the hepatopancreas and ovary. We further highlighted candidate gene modules highly linked to physiological traits. Collectively, our observations suggested that starvation triggered endogenous ovarian maturation at the expense of hepatopancreas mass, with both metabolic adjustments to optimize energy and fatty acid supply from hepatopancreas to ovary in the early phase, followed by the activation of autophagy-related pathways in both organs over prolonged starvation. These specific adaptive responses might be considered efficient strategies to stimulate ovarian maturation of Scylla paramamosain under fasting stress, which improves the nutritional value of female mud crabs and other economically important crustaceans.

Identification and validation of chromatin regulator-related signatures as a novel prognostic model for low-grade gliomas using translational bioinformatics

AIMS

The purpose of this study is to explore the potential biological role and prognostic significance of chromatin regulators (CRs) in low-grade gliomas (LGGs).

MAIN METHODS

CRs were obtained from the FACER database. Transcription profiles of LGG patients were collected from the TCGA and CGGA databases. Differentially expressed CRs (DECRs) between LGGs and normal controls were identified using DESeq2. The consensus clustering algorithm was employed to distinguish subtypes of LGGs based on prognosis-related DECRs. The differences in clinical and molecular characteristics between different subtypes were explored. R packages, GSVA, ssGSEA, and ESTIMATE were utilized to elucidate the tumor microenvironment and activated pathways in different subtypes. Subsequently, a CRs-related signature was developed using LASSO Cox regression. Its performance was evaluated by Kaplan-Meier curve and ROC curve analyses. In vitro experiments were performed to explore the function of JADE3 in LGGs, which predominantly expressed in glioma cells.

KEY FINDINGS

We identified 43 DECRs and two CRs-related subtypes of LGGs. The subtype characterized by shorter survival displayed significant enrichment for pathways associated with DNA damage response and repair, along with heightened immune cell infiltration. Furthermore, the CRs-based signature exhibited excellent prognostic performance in both the TCGA and CGGA databases. Knockdown of JADE3 significantly increased the invasion, migration, and proliferation abilities of Hs683.

SIGNIFICANCE

Our study reveals the aberrant expression and prognostic value of CRs in LGGs. It emphasizes the potential regulatory role of CRs in the microenvironment and DNA damage repair in LGGs. JADE3 could be a possible therapeutic target for LGGs.

Chromatin regulator SMARCAL1 modulates cellular lipid metabolism

Biallelic mutations of the chromatin regulator SMARCAL1 cause Schimke Immunoosseous Dysplasia (SIOD), characterized by severe growth defects and premature mortality. Atherosclerosis and hyperlipidemia are common among SIOD patients, yet their onset and progression are poorly understood. Using an integrative approach involving proteomics, mouse models, and population genetics, we investigated SMARCAL1's role. We found that SmarcAL1 interacts with angiopoietin-like 3 (Angptl3), a key regulator of lipoprotein metabolism. In vitro and in vivo analyses demonstrate SmarcAL1's vital role in maintaining cellular lipid homeostasis. The observed translocation of SmarcAL1 to cytoplasmic peroxisomes suggests a potential regulatory role in lipid metabolism through gene expression. SmarcAL1 gene inactivation reduces the expression of key genes in cellular lipid catabolism. Population genetics investigations highlight significant associations between SMARCAL1 genetic variations and body mass index, along with lipid-related traits. This study underscores SMARCAL1's pivotal role in cellular lipid metabolism, likely contributing to the observed lipid phenotypes in SIOD patients.

A comprehensive atlas of epigenetic regulators reveals tissue-specific epigenetic regulation patterns

Epigenetic machinery contributes to gene regulation in eukaryotic species. However, the machinery including more than 600 epigenetic regulator (ER) genes responsible for reading, writing, and erasing histone modifications and DNA modifications remains largely uncharacterized across species. We compile a comprehensive list of ERs based on an evolutionary analysis across 23 species, which is the most comprehensive ER list in various species until recently. We further perform comparative transcriptomic analyses across different tissues in humans, mice, as well as other amniote species. We observe a consistent tissue-of-origin expression specificity pattern of duplicated ER genes across species and suggest links between expression specificity and ER gene evolution as well as ER function. Additional analyses further suggest that ER duplication can generate tissue-specific ER genes with the same epigenetic substrates, which may be closely related to their regulatory specificity in tissue development. Our work can serve as a foundation to better comprehend the tissue-specific expression patterns of ER genes from an evolutionary perspective and also the functional implications of ERs in tissue-specific epigenetic regulation.

Pan-cancer landscape of epigenetic factor expression predicts tumor outcome

Oncogenic pathways that drive cancer progression reflect both genetic changes and epigenetic regulation. Here we stratified primary tumors from each of 24 TCGA adult cancer types based on the gene expression patterns of epigenetic factors (epifactors). The tumors for five cancer types (ACC, KIRC, LGG, LIHC, and LUAD) separated into two robust clusters that were better than grade or epithelial-to-mesenchymal transition in predicting clinical outcomes. The majority of epifactors that drove the clustering were also individually prognostic. A pan-cancer machine learning model deploying epifactor expression data for these five cancer types successfully separated the patients into poor and better outcome groups. Single-cell analysis of adult and pediatric tumors revealed that expression patterns associated with poor or worse outcomes were present in individual cells within tumors. Our study provides an epigenetic map of cancer types and lays a foundation for discovering pan-cancer targetable epifactors.

MORC3 represses the HCMV major immediate early promoter in myeloid cells in the absence of PML nuclear bodies

Human cytomegalovirus (HCMV) can undergo either a latent or a lytic infection in cells of the myeloid lineage. Whilst the molecular mechanisms which determine the outcome of infection are far from clear, it is well established that a key factor is the differential regulation of the major immediate early promoter (MIEP) responsible for driving lytic immediate early gene expression. Using a myelomonocytic cell line stably transduced with a GFP reporter under the control of the MIEP, which recapitulates MIEP regulation in the context of virus infection, we have used an unbiased CRISPR-Cas9 sub-genomic, epigenetic library screen to identify novel cellular factors involved in MIEP repression during establishment and maintenance of latency in myeloid cells. One such cellular factor identified was MORC3. Consistent with MORC3 being a robust repressor of the MIEP, we show that THP1 cells devoid of MORC3 fail to establish latency. We also show that MORC3 is induced during latent infection, recruited to the MIEP and forms MORC3 nuclear bodies (MORC3-NBs) which, interestingly, co-localize with viral genomes. Finally, we show that the latency-associated functions of MORC3 are regulated by the deSUMOylase activity of the viral latency-associated LUNA protein likely to prevent untimely HCMV reactivation.

Seven chromatin regulators as immune cell infiltration characteristics, potential diagnostic biomarkers and drugs prediction in hepatocellular carcinoma

Treatment is challenging due to the heterogeneity of hepatocellular carcinoma (HCC). Chromatin regulators (CRs) are important in epigenetics and are closely associated with HCC. We obtained HCC-related expression data and relevant clinical data from The Cancer Genome Atlas (TCGA) databases. Then, we crossed the differentially expressed genes (DEGs), immune-related genes and CRs to obtain immune-related chromatin regulators differentially expressed genes (IRCR DEGs). Least absolute shrinkage and selection operator (LASSO) Cox regression analysis was performed to select the prognostic gene and construct a risk model for predicting prognosis in HCC, followed by a correlation analysis of risk scores with clinical characteristics. Finally, we also carried out immune microenvironment analysis and drug sensitivity analysis, the correlation between risk score and clinical characteristics was analyzed. In addition, we carried out immune microenvironment analysis and drug sensitivity analysis. Functional analysis suggested that IRCR DEGs was mainly enriched in chromatin-related biological processes. We identified and validated PPARGC1A, DUSP1, APOBEC3A, AIRE, HDAC11, HMGB2 and APOBEC3B as prognostic biomarkers for the risk model construction. The model was also related to immune cell infiltration, and the expression of CD48, CTLA4, HHLA2, TNFSF9 and TNFSF15 was higher in high-risk group. HCC patients in the high-risk group were more sensitive to Axitinib, Docetaxel, Erlotinib, and Metformin. In this study, we construct a prognostic model of immune-associated chromatin regulators, which provides new ideas and research directions for the accurate treatment of HCC.

A neuron-specific microexon ablates the novel DNA-binding function of a histone H3K4me0 reader PHF21A

How cell-type-specific chromatin landscapes emerge and progress during metazoan ontogenesis remains an important question. Transcription factors are expressed in a cell-type-specific manner and recruit chromatin-regulatory machinery to specific genomic loci. In contrast, chromatin-regulatory proteins are expressed broadly and are assumed to exert the same intrinsic function across cell types. However, human genetics studies have revealed an unexpected vulnerability of neurodevelopment to chromatin factor mutations with unknown mechanisms. Here, we report that 14 chromatin regulators undergo evolutionary-conserved neuron-specific splicing events involving microexons. Of the 14 chromatin regulators, two are integral components of a histone H3K4 demethylase complex; the catalytic subunit LSD1 and an H3K4me0-reader protein PHF21A adopt neuron-specific forms. We found that canonical PHF21A (PHF21A-c) binds to DNA by AT-hook motif, and the neuronal counterpart PHF21A-n lacks this DNA-binding function yet maintains H3K4me0 recognition intact. In-vitro reconstitution of the canonical and neuronal PHF21A-LSD1 complexes identified the neuronal complex as a hypomorphic H3K4 demethylating machinery with reduced nucleosome engagement. Furthermore, an autism-associated PHF21A missense mutation, 1285 G>A, at the last nucleotide of the common exon immediately upstream of the neuronal microexon led to impaired splicing of PHF21A -n. Thus, ubiquitous chromatin regulatory complexes exert unique intrinsic functions in neurons via alternative splicing of their subunits and potentially contribute to faithful human brain development.

Epigenetic target identification strategy based on multi-feature learning

The identification of potential epigenetic targets for a known bioactive compound is essential and promising as more and more epigenetic drugs are used in cancer clinical treatment and the availability of chemogenomic data related to epigenetics increases. In this study, we introduce a novel epigenetic target identification strategy (ETI-Strategy) that integrates a multi-task graph convolutional neural network prior model and a protein-ligand interaction classification discriminating model using large-scale bioactivity data for a panel of 55 epigenetic targets. Our approach utilizes machine learning techniques to achieve an AUC value of 0.934 for the prior model and 0.830 for the discriminating model, outperforming inverse docking in predicting protein-ligand interactions. When comparing with other open-source target identification tools, it was found that only our tool was able to accurately predict all the targets corresponding to each compound. This further demonstrates the ability of our strategy to take full advantage of molecular-level information as well as protein-level information in molecular activity prediction. Our work highlights the contribution of machine learning in the identification of potential epigenetic targets and offers a novel approach for epigenetic drug discovery and development.Communicated by Ramaswamy H. Sarma.

The chromatin remodeling factors EP300 and TRRAP are novel SMYD3 interactors involved in the emerging 'nonmutational epigenetic reprogramming' cancer hallmark

SMDY3 is a histone-lysine N-methyltransferase involved in several oncogenic processes and is believed to play a major role in various cancer hallmarks. Recently, we identified ATM, BRCA2, CHK2, MTOR, BLM, MET, AMPK, and p130 as direct SMYD3 interactors by taking advantage of a library of rare tripeptides, which we first tested for their in vitro binding affinity to SMYD3 and then used as in silico probes to systematically search the human proteome. Here, we used this innovative approach to identify further SMYD3-interacting proteins involved in crucial cancer pathways and found that the chromatin remodeling factors EP300 and TRRAP interact directly with SMYD3, thus linking SMYD3 to the emerging 'nonmutational epigenetic reprogramming' cancer hallmark. Of note, we validated these interactions in gastrointestinal cancer cell lines, including HCT-116 cells, which harbor a C-terminal truncating mutation in EP300, suggesting that EP300 binds to SMYD3 via its N-terminal region. While additional studies are required to ascertain the functional mechanisms underlying these interactions and their significance, the identification of two novel SMYD3 interactors involved in epigenetic cancer hallmark pathways adds important pieces to the puzzle of how SMYD3 exerts its oncogenic role.

Methamphetamine-induced region-specific transcriptomic and epigenetic changes in the brain of male rats

Psychostimulant methamphetamine (METH) is neurotoxic to the brain and, therefore, its misuse leads to neurological and psychiatric disorders. The gene regulatory network (GRN) response to neurotoxic METH binge remains unclear in most brain regions. Here we examined the effects of binge METH on the GRN in the nucleus accumbens, dentate gyrus, Ammon's horn, and subventricular zone in male rats. At 24 h after METH, ~16% of genes displayed altered expression and over a quarter of previously open chromatin regions - parts of the genome where genes are typically active - showed shifts in their accessibility. Intriguingly, most changes were unique to each area studied, and independent regulation between transcriptome and chromatin accessibility was observed. Unexpectedly, METH differentially impacted gene activity and chromatin accessibility within the dentate gyrus and Ammon's horn. Around 70% of the affected chromatin-accessible regions in the rat brain have conserved DNA sequences in the human genome. These regions frequently act as enhancers, ramping up the activity of nearby genes, and contain mutations linked to various neurological conditions. By sketching out the gene regulatory networks associated with binge METH in specific brain regions, our study offers fresh insights into how METH can trigger profound, region-specific molecular shifts.

Expanding adverse outcome pathways towards one health models for nanosafety

The ever-growing production of nano-enabled products has generated the need for dedicated risk assessment strategies that ensure safety for humans and the environment. Transdisciplinary approaches are needed to support the development of new technologies while respecting environmental limits, as also highlighted by the EU Green Deal Chemicals Strategy for Sustainability and its safe and sustainable by design (SSbD) framework. The One Health concept offers a holistic multiscale approach for the assessment of nanosafety. However, toxicology is not yet capable of explaining the interaction between chemicals and biological systems at the multiscale level and in the context of the One Health framework. Furthermore, there is a disconnect between chemical safety assessment, epidemiology, and other fields of biology that, if unified, would enable the adoption of the One Health model. The development of mechanistic toxicology and the generation of omics data has provided important biological knowledge of the response of individual biological systems to nanomaterials (NMs). On the other hand, epigenetic data have the potential to inform on interspecies mechanisms of adaptation. These data types, however, need to be linked to concepts that support their intuitive interpretation. Adverse Outcome Pathways (AOPs) represent an evolving framework to anchor existing knowledge to chemical risk assessment. In this perspective, we discuss the possibility of integrating multi-level toxicogenomics data, including toxicoepigenetic insights, into the AOP framework. We anticipate that this new direction of toxicogenomics can support the development of One Health models applicable to groups of chemicals and to multiple species in the tree of life.

An effective prognostic model in colon adenocarcinoma composed of cuproptosis-related epigenetic regulators

Background: Colorectal adenocarcinoma (COAD) is a common malignant tumor with little effective prognostic markers. Cuproptosis is a newly discovered mode of cell death that may be related to epigenetic regulators. This study aimed to explore the association between epigenetic regulators and cuproptosis, and to establish a prognostic prediction model for COAD based on epigenetic regulators associated with cuproptosis (EACs). Methods: RNA sequencing data and clinical data of 524 COAD patients were obtained from the TCGA-COAD database, cuproptosis-related genes were from the FerrDb database, and epigenetic-related genes were from databases such as GO and EpiFactors. LASSO regression analysis and other methods were used to screen out epigenetic regulators associated with cuproptosis and prognosis. The risk score of each patient was calculated and the patients were divided into high-risk group and low-risk group. Next, the survival difference, functional enrichment analyses, tumor mutation burden, chemotherapy drug sensitivity and other indicators between the two groups were compared and analyzed. Results: We found 716 epigenetic regulators closely related to cuproptosis, among which 35 genes were related to prognosis of COAD. We further screened out 7 EACs from the 35 EACs to construct a prognostic prediction model. We calculated the risk score of each patient based on these 7 genes, and divided the patients into high-risk group and low-risk group. We found that the overall survival rate and progression-free survival rate of the high-risk group were significantly lower than those of the low-risk group. This model showed good predictive ability in the training set, test set and overall data set. We also constructed a prognostic prediction model based on risk score and other clinical features, and drew the corresponding Nomogram. In addition, we found significant differences between the high-risk group and the low-risk group in tumor mutation burden, chemotherapy drug sensitivity and other clinical aspects. Conclusion: We established an effective predictive prediction model for COAD based on EACs, revealing the association between epigenetic regulators and cuproptosis in COAD. We hope that this model can not only facilitate the treatment decision of COAD patients, but also promote the research progress in the field of cuproptosis.

Short heat shock has a long-term effect on mesenchymal stem cells' transcriptome

The adverse effects of heat stress (HS) on physiological systems are well documented, yet the underlying molecular mechanisms behind it remain poorly understood. To address this knowledge gap, we conducted a comprehensive investigation into the impact of HS on mesenchymal stem cells (MSCs), focusing on their morphology, phenotype, proliferative capacity, and fate determination. Our in-depth analysis of the MSCs' transcriptome revealed a significant influence of HS on the transcriptional landscape. Notably, even after a short period of stress, we observed a persistent alteration in cell identity, potentially mediated by the activation of bivalent genes. Furthermore, by comparing the differentially expressed genes following short HS with their transcriptional state after recovery, we identified the transient upregulation of MLL and other histone modifiers, providing a potential mechanistic explanation for the stable activation of bivalent genes. This could be used to predict and modify the long-term effect of HS on cell identity.

Integrated multiomic analysis and high-throughput screening reveal potential gene targets and synergetic drug combinations for osteosarcoma therapy

Although great advances have been made over the past decades, therapeutics for osteosarcoma are quite limited. We performed long-read RNA sequencing and tandem mass tag (TMT)-based quantitative proteome on osteosarcoma and the adjacent normal tissues, next-generation sequencing (NGS) on paired osteosarcoma samples before and after neoadjuvant chemotherapy (NACT), and high-throughput drug combination screen on osteosarcoma cell lines. Single-cell RNA sequencing data were analyzed to reveal the heterogeneity of potential therapeutic target genes. Additionally, we clarified the synergistic mechanisms of doxorubicin (DOX) and HDACs inhibitors for osteosarcoma treatment. Consequently, we identified 2535 osteosarcoma-specific genes and several alternative splicing (AS) events with osteosarcoma specificity and/or patient heterogeneity. Hundreds of potential therapeutic targets were identified among them, which showed the core regulatory roles in osteosarcoma. We also identified 215 inhibitory drugs and 236 synergistic drug combinations for osteosarcoma treatment. More interestingly, the multiomic analysis pointed out the pivotal role of HDAC1 and TOP2A in osteosarcoma. HDAC inhibitors synergized with DOX to suppress osteosarcoma both in vitro and in vivo. Mechanistically, HDAC inhibitors synergized with DOX by downregulating SP1 to transcriptionally modulate TOP2A expression. This study provided a comprehensive view of molecular features, therapeutic targets, and synergistic drug combinations for osteosarcoma.

Landscape of mSWI/SNF chromatin remodeling complex perturbations in neurodevelopmental disorders

DNA sequencing-based studies of neurodevelopmental disorders (NDDs) have identified a wide range of genetic determinants. However, a comprehensive analysis of these data, in aggregate, has not to date been performed. Here, we find that genes encoding the mammalian SWI/SNF (mSWI/SNF or BAF) family of ATP-dependent chromatin remodeling protein complexes harbor the greatest number of de novo missense and protein-truncating variants among nuclear protein complexes. Non-truncating NDD-associated protein variants predominantly disrupt the cBAF subcomplex and cluster in four key structural regions associated with high disease severity, including mSWI/SNF-nucleosome interfaces, the ATPase-core ARID-armadillo repeat (ARM) module insertion site, the Arp module and DNA-binding domains. Although over 70% of the residues perturbed in NDDs overlap with those mutated in cancer, ~60% of amino acid changes are NDD-specific. These findings provide a foundation to functionally group variants and link complex aberrancies to phenotypic severity, serving as a resource for the chromatin, clinical genetics and neurodevelopment communities.

Restricting epigenetic activity promotes the reprogramming of transformed cells to pluripotency in a line-specific manner

Somatic cell reprogramming and oncogenic transformation share surprisingly similar features, yet transformed cells are resistant to reprogramming. Epigenetic barriers must block transformed cells from reprogramming, but the nature of those barriers is unclear. In this study, we generated a systematic panel of transformed mouse embryonic fibroblasts (MEFs) using oncogenic transgenes and discovered transformed cell lines compatible with reprogramming when transfected with Oct4/Sox2/Klf4/Myc. By comparing the reprogramming-capable and incapable transformed lines we identified multiple stages of failure in the reprogramming process. Some transformed lines failed at an early stage, whilst other lines seemed to progress through a conventional reprogramming process. Finally, we show that MEK inhibition overcomes one critical reprogramming barrier by indirectly suppressing a hyperacetylated active epigenetic state. This study reveals that diverse epigenetic barriers underly resistance to reprogramming of transformed cells.

Cooperative regulation of Zhx1 and hnRNPA1 drives the cardiac progenitor-specific transcriptional activation during cardiomyocyte differentiation

The zinc finger proteins (ZNFs) mediated transcriptional regulation is critical for cell fate transition. However, it is still unclear how the ZNFs realize their specific regulatory roles in the stage-specific determination of cardiomyocyte differentiation. Here, we reported that the zinc fingers and homeoboxes 1 (Zhx1) protein, transiently expressed during the cell fate transition from mesoderm to cardiac progenitors, was indispensable for the proper cardiomyocyte differentiation of mouse and human embryonic stem cells. Moreover, Zhx1 majorly promoted the specification of cardiac progenitors via interacting with hnRNPA1 and co-activated the transcription of a wide range of genes. In-depth mechanistic studies showed that Zhx1 was bound with hnRNPA1 by the amino acid residues (Thr111-His120) of the second Znf domain, thus participating in the formation of cardiac progenitors. Together, our study highlights the unrevealed interaction of Zhx1/hnRNPA1 for activating gene transcription during cardiac progenitor specification and also provides new evidence for the specificity of cell fate determination in cardiomyocyte differentiation.

Oxygen Level Modifies the Expression of Genes Involved in the Epigenetic Regulation of Multipotent Stromal Cells In Vitro

Changes in the transcriptional activity of genes involved in the epigenetic regulation of adipose tissue multipotent mesenchymal stromal cells were analyzed in vitro at different O2 levels. DNA microarray study showed that the most pronounced changes in gene expression, including genes responsible for the epigenetic regulation of mesenchymal stromal cells, occurred at 3% O2. A lower number of genes changed the expression at 1% O2, and a minimum response was observed at 5% O2 in comparison with standard culturing conditions (20% O2). The greatest number of differentially expressed genes were genes responsible for the regulation of histones; the genes encoding products that regulate chromatin, DNA, and RNA constituted a lower part. Thus, the degree of hypoxia can modify the response of multipotent mesenchymal stromal cells at the level of epigenetic regulators.

Tumor subtypes and signature model construction based on chromatin regulators for better prediction of prognosis in uveal melanoma

Background: Uveal Melanoma (UM) is the most prevalent primary intraocular malignancy in adults. This study assessed the importance of chromatin regulators (CRs) in UM and developed a model to predict UM prognosis. Methods: Gene expression data and clinical information for UM were obtained from public databases. Samples were typed according to the gene expression of CRs associated with UM prognosis. The prognostic key genes were further screened by the protein interaction network, and the risk model was to predict UM prognosis using the least absolute shrinkage and selection operator (LASSO) regression analysis and performed a test of the risk mode. In addition, we performed gene set variation analysis, tumor microenvironment, and tumor immune analysis between subtypes and risk groups to explore the mechanisms influencing the development of UM. Results: We constructed a signature model consisting of three CRs (RUVBL1, SIRT3, and SMARCD3), which was shown to be accurate, and valid for predicting prognostic outcomes in UM. Higher immune cell infiltration in poor prognostic subtypes and risk groups. The Tumor immune analysis and Tumor Immune Dysfunction and Exclusion (TIDE) score provided a basis for clinical immunotherapy in UM. Conclusion: The risk model has prognostic value for UM survival and provides new insights into the treatment of UM.

A polyphenol-rich green Mediterranean diet enhances epigenetic regulatory potential: the DIRECT PLUS randomized controlled trial

BACKGROUND

The capacity of a polyphenol-enriched diet to modulate the epigenome in vivo is partly unknown. Given the beneficial metabolic effects of a Mediterranean (MED) diet enriched in polyphenols and reduced in red/processed meat (green-MED), as previously been proven by the 18-month DIRECT PLUS randomized controlled trial, we analyzed the effects of the green-MED diet on methylome and transcriptome levels to highlight molecular mechanisms underlying the observed metabolic improvements.

METHODS

Our study included 260 participants (baseline BMI = 31.2 kg/m², age = 5 years) of the DIRECT PLUS trial, initially randomized to one of the intervention arms: A. healthy dietary guidelines (HDG), B. MED (440 mg polyphenols additionally provided by walnuts), C. green-MED (1240 mg polyphenols additionally provided by walnuts, green tea, and Mankai: green duckweed shake). Blood methylome and transcriptome of all study subjects were analyzed at baseline and after completing the 18-month intervention using Illumina EPIC and RNA sequencing technologies.

RESULTS

A total of 1573 differentially methylated regions (DMRs; false discovery rate (FDR) < 5 %) were found in the green-MED compared to the MED (177) and HDG (377) diet participants. This corresponded to 1753 differentially expressed genes (DEGs; FDR < 5 %) in the green-MED intervention compared to MED (7) and HDG (738). Consistently, the highest number (6 %) of epigenetic modulating genes was transcriptionally changed in subjects participating in the green-MED intervention. Weighted cluster network analysis relating transcriptional and phenotype changes among participants subjected to the green-MED intervention identified candidate genes associated with serum-folic acid change (all P < 1 × 10^-3) and highlighted one module including the KIR3DS1 locus, being negatively associated with the polyphenol changes (e.g. P < 1 × 10^-4), but positively associated with the MRI-assessed superficial subcutaneous adipose area-, weight- and waist circumference- 18-month change (all P < 0.05). Among others, this module included the DMR gene Cystathionine Beta-Synthase, playing a major role in homocysteine reduction.

CONCLUSIONS

The green-MED high polyphenol diet, rich in green tea and Mankai, renders a high capacity to regulate an individual's epigenome. Our findings suggest epigenetic key drivers such as folate and green diet marker to mediate this capacity and indicate a direct effect of dietary polyphenols on the one‑carbon metabolism.

Epigenetic regulation of human-specific gene expression in the prefrontal cortex

BACKGROUND

Changes in gene expression levels during brain development are thought to have played an important role in the evolution of human cognition. With the advent of high-throughput sequencing technologies, changes in brain developmental expression patterns, as well as human-specific brain gene expression, have been characterized. However, interpreting the origin of evolutionarily advanced cognition in human brains requires a deeper understanding of the regulation of gene expression, including the epigenomic context, along the primate genome. Here, we used chromatin immunoprecipitation sequencing (ChIP-seq) to measure the genome-wide profiles of histone H3 lysine 4 trimethylation (H3K4me3) and histone H3 lysine 27 acetylation (H3K27ac), both of which are associated with transcriptional activation in the prefrontal cortex of humans, chimpanzees, and rhesus macaques.

RESULTS

We found a discrete functional association, in which ^H3K4me3HP gain was significantly associated with myelination assembly and signaling transmission, while ^H3K4me3HP loss played a vital role in synaptic activity. Moreover, ^H3K27acHP gain was enriched in interneuron and oligodendrocyte markers, and ^H3K27acHP loss was enriched in CA1 pyramidal neuron markers. Using strand-specific RNA sequencing (ssRNA-seq), we first demonstrated that approximately 7 and 2% of human-specific expressed genes were epigenetically marked by ^H3K4me3HP and ^H3K27acHP, respectively, providing robust support for causal involvement of histones in gene expression. We also revealed the co-activation role of epigenetic modification and transcription factors in human-specific transcriptome evolution. Mechanistically, histone-modifying enzymes at least partially contribute to an epigenetic disturbance among primates, especially for the H3K27ac epigenomic marker. In line with this, peaks enriched in the macaque lineage were found to be driven by upregulated acetyl enzymes.

CONCLUSIONS

Our results comprehensively elucidated a causal species-specific gene-histone-enzyme landscape in the prefrontal cortex and highlighted the regulatory interaction that drove transcriptional activation.

CTCF and Its Multi-Partner Network for Chromatin Regulation

Architectural proteins are essential epigenetic regulators that play a critical role in organizing chromatin and controlling gene expression. CTCF (CCCTC-binding factor) is a key architectural protein responsible for maintaining the intricate 3D structure of chromatin. Because of its multivalent properties and plasticity to bind various sequences, CTCF is similar to a Swiss knife for genome organization. Despite the importance of this protein, its mechanisms of action are not fully elucidated. It has been hypothesized that its versatility is achieved through interaction with multiple partners, forming a complex network that regulates chromatin folding within the nucleus. In this review, we delve into CTCF's interactions with other molecules involved in epigenetic processes, particularly histone and DNA demethylases, as well as several long non-coding RNAs (lncRNAs) that are able to recruit CTCF. Our review highlights the importance of CTCF partners to shed light on chromatin regulation and pave the way for future exploration of the mechanisms that enable the finely-tuned role of CTCF as a master regulator of chromatin.

Comprehensive analyses reveal the prognosis and biological function roles of chromatin regulators in lung adenocarcinoma

The present study explored the prognosis and biological function roles of chromatin regulators (CRs) in patients with lung adenocarcinoma (LUAD). Using transcriptome profile and clinical follow-up data of LUAD dataset, we explored the molecular classification, developed, and validated a CR prognostic model, built an individual risk scoring system in LUAD, and compared the clinical and molecular characteristics between different subtypes and risk stratifications. We investigated the chemotherapy sensitivity and predicted potential immunotherapy response. Lastly, we collected the clinical samples and validated the prognosis and potential function role of NAPS2. Our study indicated that LUAD patients could be classified into two subtypes that had obviously different clinical background and molecular features. We constructed a prognostic model with eight CR genes, which was well validated in several other population cohort. We built high- and low-risk stratifications for LUAD patients. Patients from high-risk group were totally different from low-risk groups in clinical, biological function, gene mutation, microenvironment, and immune infiltration levels. We idented several potential molecular compounds for high-risk group treatment. We predicted that high-risk group may have poor immunotherapy response. We finally found that Neuronal PAS Domain Protein 2 (NPAS2) involved in the progression of LUAD via regulating cell adhesion. Our study indicated that CR involved in the progression of LUAD and affect their prognosis. Different therapeutic strategies should be developed for different molecular subtypes and risk stratifications. Our comprehensive analyses uncover specific determinants of CRs in LUAD and provides implications for investigating disease-associated CRs.

Epigenomic signatures reveal mechanistic clues and predictive markers for autism spectrum disorder

Autism spectrum disorder (ASD) comprises a heterogeneous group of neurodevelopmental outcomes in children with a commonality in deficits in social communication and language combined with repetitive behaviors and interests. The etiology of ASD is heterogeneous, as several hundred genes have been implicated as well as multiple in utero environmental exposures. Over the past two decades, epigenetic investigations, including DNA methylation, have emerged as a novel way to capture the complex interface of multivariate ASD etiologies. More recently, epigenome-wide association studies using human brain and surrogate accessible tissues have revealed some convergent genes that are epigenetically altered in ASD, many of which overlap with known genetic risk factors. Unlike transcriptomes, epigenomic signatures defined by DNA methylation from surrogate tissues such as placenta and cord blood can reflect past differences in fetal brain gene transcription, transcription factor binding, and chromatin. For example, the discovery of NHIP (neuronal hypoxia inducible, placenta associated) through an epigenome-wide association in placenta, identified a common genetic risk for ASD that was modified by prenatal vitamin use. While epigenomic signatures are distinct between different genetic syndromic causes of ASD, bivalent chromatin and some convergent gene pathways are consistently epigenetically altered in both syndromic and idiopathic ASD, as well as some environmental exposures. Together, these epigenomic signatures hold promising clues towards improved early prediction and prevention of ASD as well genes and gene pathways to target for pharmacological interventions. Future advancements in single cell and multi-omic technologies, machine learning, as well as non-invasive screening of epigenomic signatures during pregnancy or newborn periods are expected to continue to impact the translatability of the recent discoveries in epigenomics to precision public health.

Annelid adult cell type diversity and their pluripotent cellular origins

Annelids are a broadly distributed, highly diverse, economically and environmentally important group of animals. Most species can regenerate missing body parts, and many are able to reproduce asexually. Therefore, many annelids can generate all adult cell types in adult stages. However, the putative adult stem cell populations involved in these processes, as well as the diversity of adult cell types generated by them, are still unknown. Here, we recover 75,218 single cell transcriptomes of Pristina leidyi, a highly regenerative and asexually-reproducing freshwater annelid. We characterise all major annelid adult cell types, and validate many of our observations by HCR in situ hybridisation. Our results uncover complex patterns of regionally expressed genes in the annelid gut, as well as neuronal, muscle and epidermal specific genes. We also characterise annelid-specific cell types such as the chaetal sacs and globin+ cells, and novel cell types of enigmatic affinity, including a vigilin+ cell type, a lumbrokinase+ cell type, and a diverse set of metabolic cells. Moreover, we characterise transcription factors and gene networks that are expressed specifically in these populations. Finally, we uncover a broadly abundant cluster of putative stem cells with a pluripotent signature. This population expresses well-known stem cell markers such as vasa, piwi and nanos homologues, but also shows heterogeneous expression of differentiated cell markers and their transcription factors. In these piwi+ cells, we also find conserved expression of pluripotency regulators, including multiple chromatin remodelling and epigenetic factors. Finally, lineage reconstruction analyses reveal the existence of differentiation trajectories from piwi+ cells to diverse adult types. Our data reveal the cell type diversity of adult annelids for the first time and serve as a resource for studying annelid cell types and their evolution. On the other hand, our characterisation of a piwi+ cell population with a pluripotent stem cell signature will serve as a platform for the study of annelid stem cells and their role in regeneration.

Epigenetic and Tumor Microenvironment for Prognosis of Patients with Gastric Cancer

BACKGROUND

Epigenetics studies heritable or inheritable mechanisms that regulate gene expression rather than altering the DNA sequence. However, no research has investigated the link between TME-related genes (TRGs) and epigenetic-related genes (ERGs) in GC.

METHODS

A complete review of genomic data was performed to investigate the relationship between the epigenesis tumor microenvironment (TME) and machine learning algorithms in GC.

RESULTS

Firstly, TME-related differential expression of genes (DEGs) performed non-negative matrix factorization (NMF) clustering analysis and determined two clusters (C1 and C2). Then, Kaplan-Meier curves for overall survival (OS) and progression-free survival (PFS) rates suggested that cluster C1 predicted a poorer prognosis. The Cox-LASSO regression analysis identified eight hub genes (SRMS, MET, OLFML2B, KIF24, CLDN9, RNF43, NETO2, and PRSS21) to build the TRG prognostic model and nine hub genes (TMPO, SLC25A15, SCRG1, ISL1, SOD3, GAD1, LOXL4, AKR1C2, and MAGEA3) to build the ERG prognostic model. Additionally, the signature's area under curve (AUC) values, survival rates, C-index scores, and mean squared error (RMS) curves were evaluated against those of previously published signatures, which revealed that the signature identified in this study performed comparably. Meanwhile, based on the IMvigor210 cohort, a statistically significant difference in OS between immunotherapy and risk scores was observed. It was followed by LASSO regression analysis which identified 17 key DEGs and a support vector machine (SVM) model identified 40 significant DEGs, and based on the Venn diagram, eight co-expression genes (ENPP6, VMP1, LY6E, SHISA6, TMEM158, SYT4, IL11, and KLK8) were discovered.

CONCLUSION

The study identified some hub genes that could be useful in predicting prognosis and management in GC.

Large-scale mapping and mutagenesis of human transcriptional effector domains

Human gene expression is regulated by more than 2,000 transcription factors and chromatin regulators1,2. Effector domains within these proteins can activate or repress transcription. However, for many of these regulators we do not know what type of effector domains they contain, their location in the protein, their activation and repression strengths, and the sequences that are necessary for their functions. Here, we systematically measure the effector activity of more than 100,000 protein fragments tiling across most chromatin regulators and transcription factors in human cells (2,047 proteins). By testing the effect they have when recruited at reporter genes, we annotate 374 activation domains and 715 repression domains, roughly 80% of which are new and have not been previously annotated3-5. Rational mutagenesis and deletion scans across all the effector domains reveal aromatic and/or leucine residues interspersed with acidic, proline, serine and/or glutamine residues are necessary for activation domain activity. Furthermore, most repression domain sequences contain sites for small ubiquitin-like modifier (SUMO)ylation, short interaction motifs for recruiting corepressors or are structured binding domains for recruiting other repressive proteins. We discover bifunctional domains that can both activate and repress, some of which dynamically split a cell population into high- and low-expression subpopulations. Our systematic annotation and characterization of effector domains provide a rich resource for understanding the function of human transcription factors and chromatin regulators, engineering compact tools for controlling gene expression and refining predictive models of effector domain function.

Comprehensive characterization of the embryonic factor LEUTX

The paired-like homeobox transcription factor LEUTX is expressed in human preimplantation embryos between the 4- and 8-cell stages, and then silenced in somatic tissues. To characterize the function of LEUTX, we performed a multiomic characterization of LEUTX using two proteomics methods and three genome-wide sequencing approaches. Our results show that LEUTX stably interacts with the EP300 and CBP histone acetyltransferases through its 9 amino acid transactivation domain (9aaTAD), as mutation of this domain abolishes the interactions. LEUTX targets genomic cis-regulatory sequences that overlap with repetitive elements, and through these elements it is suggested to regulate the expression of its downstream genes. We find LEUTX to be a transcriptional activator, upregulating several genes linked to preimplantation development as well as 8-cell-like markers, such as DPPA3 and ZNF280A. Our results support a role for LEUTX in preimplantation development as an enhancer binding protein and as a potent transcriptional activator.

A neurodevelopmental epigenetic programme mediated by SMARCD3-DAB1-Reelin signalling is hijacked to promote medulloblastoma metastasis

How abnormal neurodevelopment relates to the tumour aggressiveness of medulloblastoma (MB), the most common type of embryonal tumour, remains elusive. Here we uncover a neurodevelopmental epigenomic programme that is hijacked to induce MB metastatic dissemination. Unsupervised analyses of integrated publicly available datasets with our newly generated data reveal that SMARCD3 (also known as BAF60C) regulates Disabled 1 (DAB1)-mediated Reelin signalling in Purkinje cell migration and MB metastasis by orchestrating cis-regulatory elements at the DAB1 locus. We further identify that a core set of transcription factors, enhancer of zeste homologue 2 (EZH2) and nuclear factor I X (NFIX), coordinates with the cis-regulatory elements at the SMARCD3 locus to form a chromatin hub to control SMARCD3 expression in the developing cerebellum and in metastatic MB. Increased SMARCD3 expression activates Reelin-DAB1-mediated Src kinase signalling, which results in a MB response to Src inhibition. These data deepen our understanding of how neurodevelopmental programming influences disease progression and provide a potential therapeutic option for patients with MB.

Phthalates impact on the epigenetic factors contributed specifically by the father at fertilization

BACKGROUND

Preconception exposure to phthalates such as the anti-androgenic dibutyl-phthalate (DBP) impacts both male and female reproduction, yet how this occurs largely remains unknown. Previously we defined a series of RNAs expressly provided by sperm at fertilization and separately, and in parallel, those that responded to high DBP exposure. Utilizing both populations of RNAs, we now begin to unravel the impact of high-DBP exposure on those RNAs specifically delivered by the father.

RESULTS

Enrichment of RNAs altered by DBP exposure within the Molecular Signature Database highlighted cellular stress, cell cycle, apoptosis, DNA damage response, and gene regulation pathways. Overlap within each of these five pathways identified those RNAs that were specifically (≥ fivefold enriched) or primarily (≥ twofold enriched) provided as part of the paternal contribution compared to the oocyte at fertilization. Key RNAs consistently altered by DBP, including CAMTA2 and PSME4, were delivered by sperm reflective of these pathways. The majority (64/103) of overlapping enriched gene sets were related to gene regulation. Many of these RNAs (45 RNAs) corresponded to key interconnected CRREWs (Chromatin remodeler cofactors, RNA interactors, Readers, Erasers, and Writers). Modeling suggests that CUL2, PHF10, and SMARCC1 may coordinate and mechanistically modulate the phthalate response.

CONCLUSIONS

Mediated through a CRREW regulatory network, the cell responded to exposure presenting stressed-induced changes in the cell cycle-DNA damage-apoptosis. Interestingly, the majority of these DBP-responsive epigenetic mediators' direct acetylation or deacetylation, impacting the sperm's cargo delivered at fertilization and that of the embryo.

Defining a Correlative Transcriptional Signature Associated with Bulk Histone H3 Acetylation Levels in Adult Glioblastomas

Glioblastoma (GB) is the most prevalent primary brain cancer and the most aggressive form of glioma because of its poor prognosis and high recurrence. To confirm the importance of epigenetics in glioma, we explored The Cancer Gene Atlas (TCGA) database and we found that several histone/DNA modifications and chromatin remodeling factors were affected at transcriptional and genetic levels in GB compared to lower-grade gliomas. We associated these alterations in our own cohort of study with a significant reduction in the bulk levels of acetylated lysines 9 and 14 of histone H3 in high-grade compared to low-grade tumors. Within GB, we performed an RNA-seq analysis between samples exhibiting the lowest and highest levels of acetylated H3 in the cohort; these results are in general concordance with the transcriptional changes obtained after histone deacetylase (HDAC) inhibition of GB-derived cultures that affected relevant genes in glioma biology and treatment (e.g., A2ML1, CD83, SLC17A7, TNFSF18). Overall, we identified a transcriptional signature linked to histone acetylation that was potentially associated with good prognosis, i.e., high overall survival and low rate of somatic mutations in epigenetically related genes in GB. Our study identifies lysine acetylation as a key defective histone modification in adult high-grade glioma, and offers novel insights regarding the use of HDAC inhibitors in therapy.

CRdb: a comprehensive resource for deciphering chromatin regulators in human

Chromatin regulators (CRs) regulate epigenetic patterns on a partial or global scale, playing a critical role in affecting multi-target gene expression. As chromatin immunoprecipitation sequencing (ChIP-seq) data associated with CRs are rapidly accumulating, a comprehensive resource of CRs needs to be built urgently for collecting, integrating, and processing these data, which can provide abundant annotated information on CR upstream and downstream regulatory analyses as well as CR-related analysis functions. This study established an integrative CR resource, named CRdb (http://cr.liclab.net/crdb/), with the aim of curating a large number of available resources for CRs and providing extensive annotations and analyses of CRs to help biological researchers clarify the regulation mechanism and function of CRs. The CRdb database comprised a total of 647 CRs and 2,591 ChIP-seq samples from more than 300 human tissues and cell types. These samples have been manually curated from NCBI GEO/SRA and ENCODE. Importantly, CRdb provided the abundant and detailed genetic annotations in CR-binding regions based on ChIP-seq. Furthermore, CRdb supported various functional annotations and upstream regulatory information on CRs. In particular, it embedded four types of CR regulatory analyses: CR gene set enrichment, CR-binding genomic region annotation, CR-TF co-occupancy analysis, and CR regulatory axis analysis. CRdb is a useful and powerful resource that can help in exploring the potential functions of CRs and their regulatory mechanism in diseases and biological processes.

Databases and prospects of dynamic gene regulation in eukaryotes: A mini review

In eukaryotes, dynamic regulation enables DNA polymerases to catalyze a variety of RNA products in spatial and temporal patterns. Dynamic gene expression is regulated by transcription factors (TFs) and epigenetics (DNA methylation and histone modification). The applications of biochemical technology and high-throughput sequencing enhance the understanding of mechanisms of these regulations and affected genomic regions. To provide a searchable platform for retrieving such metadata, numerous databases have been developed based on the integration of genome-wide maps (e.g., ChIP-seq, whole-genome bisulfite sequencing, RNA-seq, ATAC-seq, DNase-seq, and MNase-seq data) and functionally genomic annotation. In this mini review, we summarize the main functions of TF-related databases and outline the prevalent approaches used in inferring epigenetic regulations, their associated genes, and functions. We review the literature on crosstalk between TF and epigenetic regulation and the properties of non-coding RNA regulation, which are challenging topics that promise to pave the way for advances in database development.

EpiFactors 2022: expansion and enhancement of a curated database of human epigenetic factors and complexes

We present an update of EpiFactors, a manually curated database providing information about epigenetic regulators, their complexes, targets, and products which is openly accessible at http://epifactors.autosome.org. An updated version of the EpiFactors contains information on 902 proteins, including 101 histones and protamines, and, as a main update, a newly curated collection of 124 lncRNAs involved in epigenetic regulation. The amount of publications concerning the role of lncRNA in epigenetics is rapidly growing. Yet, the resource that compiles, integrates, organizes, and presents curated information on lncRNAs in epigenetics is missing. EpiFactors fills this gap and provides data on epigenetic regulators in an accessible and user-friendly form. For 820 of the genes in EpiFactors, we include expression estimates across multiple cell types assessed by CAGE-Seq in the FANTOM5 project. In addition, the updated EpiFactors contains information on 73 protein complexes involved in epigenetic regulation. Our resource is practical for a wide range of users, including biologists, bioinformaticians and molecular/systems biologists.

A novel chromatin regulator-related immune checkpoint related gene prognostic signature and potential candidate drugs for endometrial cancer patients

Background

Endometrial cancer (EC) is the most common gynecologic malignancy in developed countries and its prevalence is increasing. As an emerging therapy with a promising efficacy, immunotherapy has been extensively applied in the treatment of solid tumors. In addition, chromatin regulators (CRs), as essential upstream regulators of epigenetics, play a significant role in tumorigenesis and cancer development.

Methods

CRs and immune checkpoint-related genes (ICRGs) were obtained from the previous top research. The Genome Cancer Atlas (TCGA) was utilized to acquire the mRNA expression and clinical information of patients with EC. Correlation analysis was utilized for screen CRs-related ICRGs (CRRICRGs). By Cox regression and least absolute shrinkage and selection operator (LASSO) analysis, prognosis related CRRICRGs were screened out and risk model was constructed. The Kaplan–Meier curve was used to estimate the prognosis between high- and low-risk group. By comparing the IC50 value, the drugs sensitivity difference was explored. We obtained small molecule drugs for the treatment of UCEC patients based on CAMP dataset.

Results

We successfully constructed a 9 CRRICRs-based prognostic signature for patients with UCEC and found the riskscore was an independent prognostic factor. The results of functional analysis suggested that CRRICRGs may be involved in immune processes associated with cancer. Immune characteristics analysis provided further evidence that the CRRICRGs-based model was correlated with immune cells infiltration and immune checkpoint. Eight small molecule drugs that may be effective for the treatment of UCEC patients were screened. Effective drugs identified by drug sensitivity profiling in high- and low-risk groups.

Conclusion

In summary, our study provided novel insights into the function of CRRICRGs in UCEC. We also developed a reliable prognostic panel for the survival of patients with UCEC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s41065-022-00253-w.