PRC2, Chromatin Regulation, and Human Disease: Insights From Molecular Structure and Function

Biochemical properties of chromatin domains define genome compartmentalization

Chromatin three-dimensional (3D) organization inside the cell nucleus determines the separation of euchromatin and heterochromatin domains. Their segregation results in the definition of active and inactive chromatin compartments, whereby the local concentration of associated proteins, RNA and DNA results in the formation of distinct subnuclear structures. Thus, chromatin domains spatially confined in a specific 3D nuclear compartment are expected to share similar epigenetic features and biochemical properties, in terms of accessibility and solubility. Based on this rationale, we developed the 4f-SAMMY-seq to map euchromatin and heterochromatin based on their accessibility and solubility, starting from as little as 10 000 cells. Adopting a tailored bioinformatic data analysis approach we reconstruct also their 3D segregation in active and inactive chromatin compartments and sub-compartments, thus recapitulating the characteristic properties of distinct chromatin states. A key novelty of the new method is the capability to map both the linear segmentation of open and closed chromatin domains, as well as their compartmentalization in one single experiment.

Nonlinear DNA methylation trajectories in aging male mice

Although DNA methylation data yields highly accurate age predictors, little is known about the dynamics of this quintessential epigenomic biomarker during lifespan. To narrow the gap, we investigate the methylation trajectories of male mouse colon at five different time points of aging. Our study indicates the existence of sudden hypermethylation events at specific stages of life. Precisely, we identify two epigenomic switches during early-to-midlife (3-9 months) and mid-to-late-life (15-24 months) transitions, separating the rodents' life into three stages. These nonlinear methylation dynamics predominantly affect genes associated with the nervous system and enrich in bivalently marked chromatin regions. Based on groups of nonlinearly modified loci, we construct a clock-like classifier STageR (STage of aging estimatoR) that accurately predicts murine epigenetic stage. We demonstrate the universality of our clock in an independent mouse cohort and with publicly available datasets.

Early detection of malignant and pre-malignant peripheral nerve tumors using cell-free DNA fragmentomics

Early detection of neurofibromatosis type 1 (NF1) associated peripheral nerve sheath tumors (PNST) informs clinical decision-making, potentially averting deadly outcomes. Here, we describe a cell-free DNA (cfDNA) fragmentomic approach which distinguishes non-malignant, pre-malignant and malignant forms of NF1 PNST. Using plasma samples from a novel cohort of 101 NF1 patients and 21 healthy controls, we validated that our previous cfDNA copy number alteration (CNA)-based approach identifies malignant peripheral nerve sheath tumor (MPNST) but cannot distinguish among benign and premalignant states. We therefore investigated the ability of fragment-based cfDNA features to differentiate NF1-associated tumors including binned genome-wide fragment length ratios, end motif analysis, and non-negative matrix factorization deconvolution of fragment lengths. Fragmentomic methods were able to differentiate pre-malignant states including atypical neurofibromas (AN). Fragmentomics also adjudicated AN cases suspicious for MPNST, correctly diagnosing samples noninvasively, which could have informed clinical management. Overall, this study pioneers the early detection of malignant and premalignant peripheral nerve sheath tumors in NF1 patients using plasma cfDNA fragmentomics. In addition to screening applications, this novel approach distinguishes atypical neurofibromas from benign plexiform neurofibromas and malignant peripheral nerve sheath tumors, enabling more precise clinical diagnosis and management.

MAL expression downregulation through suppressive H3K27me3 marks at the promoter in HPV16-related cervical cancers is prognostically relevant and manifested by the interplay of novel MAL antisense long noncoding RNA AC103563.8, E7 oncoprotein and EZH2

BACKGROUND

MAL (T-lymphocyte maturation-associated protein) is highly downregulated in most cancers, including cervical cancer (CaCx), attributable to promoter hypermethylation. Long noncoding RNA genes (lncGs) play pivotal roles in CaCx pathogenesis, by interacting with human papillomavirus (HPV)-encoded oncoproteins, and epigenetically regulating coding gene expression. Hence, we attempted to decipher the impact and underlying mechanisms of MAL downregulation in HPV16-related CaCx pathogenesis, by interrogating the interactive roles of MAL antisense lncRNA AC103563.8, E7 oncoprotein and PRC2 complex protein, EZH2.

RESULTS

Employing strand-specific RNA-sequencing, we confirmed the downregulated expression of MAL in association with poor overall survival of CaCx patients bearing HPV16, along with its antisense long noncoding RNA (lncRNA) AC103563.8. The strength of positive correlation between MAL and AC103563.8 was significantly high among patients compared to normal individuals. While downregulated expression of MAL was significantly associated with poor overall survival of CaCx patients bearing HPV16, AC103563.8 did not reveal any such association. We confirmed the enrichment of chromatin suppressive mark, H3K27me3 at MAL promoter, using ChIP-qPCR in HPV16-positive SiHa cells. Subsequent E7 knockdown in such cells significantly increased MAL expression, concomitant with decreased EZH2 expression and H3K27me3 marks at MAL promoter. In silico analysis revealed that both E7 and EZH2 bear the potential of interacting with AC103563.8, at the same binding domain. RNA immunoprecipitation with anti-EZH2 and anti-E7 antibodies, respectively, and subsequent quantitative PCR analysis in E7-silenced and unperturbed SiHa cells confirmed the interaction of AC103563.8 with EZH2 and E7, respectively. Apparently, AC103563.8 seems to preclude EZH2 and bind with E7, failing to block EZH2 function in patients. Thereby, enhanced EZH2 expression in the presence of E7 could potentially inactivate the MAL promoter through H3K27me3 marks, corroborating our previous results of MAL expression downregulation in patients.

CONCLUSION

AC103563.8-E7-EZH2 axis, therefore, appears to crucially regulate the expression of MAL, through chromatin inactivation in HPV16-CaCx pathogenesis, warranting therapeutic strategy development.

Integrated single-cell multiomics uncovers foundational regulatory mechanisms of lens development and pathology

Ocular lens development entails epithelial to fiber cell differentiation, defects in which cause congenital cataracts. We report the first single-cell multiomic atlas of lens development, leveraging snRNA-seq, snATAC-seq and CUT&RUN-seq to discover previously unreported mechanisms of cell fate determination and cataract-linked regulatory networks. A comprehensive profile of cis- and trans-regulatory interactions, including for the cataract-linked transcription factor MAF, is established across a temporal trajectory of fiber cell differentiation. Furthermore, we identify an epigenetic paradigm of cellular differentiation, defined by progressive loss of the H3K27 methylation writer Polycomb repressive complex 2 (PRC2). PRC2 localizes to heterochromatin domains across master-regulator transcription factor gene bodies, suggesting it safeguards epithelial cell fate. Moreover, we demonstrate that FGF hyper-stimulation in vivo leads to MAF network activation and the emergence of novel lens cell states. Collectively, these data depict a comprehensive portrait of lens fiber cell differentiation, while defining regulatory effectors of cell identity and cataract formation.

Estrogen promotes the proliferation and migration of endometrial cancer cells by upregulating the expression of lncRNA HOTAIR

OBJECTIVE

Estrogen (E2) is the main contributor to the progression of endometrial cancer (EC). The long noncoding RNA HOX antisense intergenic RNA (HOTAIR) is emerging as a new regulator in several cancer types. This study aimed to investigate the role of HOTAIR in EC development and identify the underlying molecular mechanisms.

METHODS

HOTAIR expression levels in human EC tissues and the corresponding adjacent tissues and human EC Ishikawa cells were determined by quantitative PCR. Ishikawa cells were treated with E2 or estrogen receptor (ER) inhibitor ICI182780, transfected with siHOTAIR oligo, or infected with lentivirus expressing shHOTAIR/shNC, alone or in combinations. The protein expression of polycomb repressive complex 2 (PRC2) was evaluated by western blotting, and cell migration was measured by transwell assays. A xenograft tumorigenic model was established by inoculating control or stable shHOTAIR-infected Ishikawa cells into nude mice and implanting 17β-estradiol release pellets.

RESULTS

HOTAIR expression was significantly elevated in human EC tissues. E2 exposure markedly increased HOTAIR levels in Ishikawa cells. Notably, E2 increased the protein expression of PRC2 and promoted EC cell migration, which were dependent on HOTAIR expression, as HOTAIR knockdown abolished these effects of E2. Similarly, E2 promoted the in vivo proliferation of grafted Ishikawa cells via upregulated HOTAIR expression in nude mice.

CONCLUSIONS

Human EC tissues highly express HOTAIR, and E2-induced EC progression depends on HOTAIR expression. This work suggests that the E2-HOTAIR axis is a potential therapeutic target in EC therapy.

Tau protein modulates an epigenetic mechanism of cellular senescence in human SH-SY5Y neuroblastoma cells

Introduction: Progressive Tau deposition in neurofibrillary tangles and neuropil threads is the hallmark of tauopathies, a disorder group that includes Alzheimer's disease. Since Tau is a microtubule-associated protein, a prevalent concept to explain the pathogenesis of tauopathies is that abnormal Tau modification contributes to dissociation from microtubules, assembly into multimeric β-sheets, proteotoxicity, neuronal dysfunction and cell loss. Tau also localizes in the cell nucleus and evidence supports an emerging function of Tau in DNA stability and epigenetic modulation. Methods: To better characterize the possible role of Tau in regulation of chromatin compaction and subsequent gene expression, we performed a bioinformatics analysis of transcriptome data obtained from Tau-depleted human neuroblastoma cells. Results: Among the transcripts deregulated in a Tau-dependent manner, we found an enrichment of target genes for the polycomb repressive complex 2. We further describe decreased cellular amounts of the core components of the polycomb repressive complex 2 and lower histone 3 trimethylation in Tau deficient cells. Among the de-repressed polycomb repressive complex 2 target gene products, IGFBP3 protein was found to be linked to increased senescence induction in Tau-deficient cells. Discussion: Our findings propose a mechanism for Tau-dependent epigenetic modulation of cell senescence, a key event in pathologic aging.

Activation of AKT induces EZH2-mediated β-catenin trimethylation in colorectal cancer

Colorectal cancer (CRC) develops in part through the deregulation of different signaling pathways, including activation of the WNT/β-catenin and PI3K/AKT pathways. Additionally, the lysine methyltransferase enhancer of zeste homologue 2 (EZH2) is commonly overexpressed in CRC. EZH2 canonically represses gene transcription by trimethylating lysine 27 of histone H3, but also has non-histone substrates. Here, we demonstrated that in CRC, active AKT phosphorylated EZH2 on serine 21. Phosphorylation of EZH2 by AKT induced EZH2 to interact with and methylate β-catenin at lysine 49, which increased β-catenin's binding to the chromatin. Additionally, EZH2-mediated β-catenin trimethylation induced β-catenin to interact with TCF1 and RNA polymerase II and resulted in dramatic gains in genomic regions with β-catenin occupancy. EZH2 catalytic inhibition decreased stemness but increased migratory phenotypes of CRC cells with active AKT. Overall, we demonstrated that EZH2 modulates AKT-induced changes in gene expression through the AKT/EZH2/β-catenin axis in CRC.

Integrative Multi-Omics Analysis of Oncogenic EZH2 Mutants: From Epigenetic Reprogramming to Molecular Signatures

Somatic heterozygous mutations in the active site of the enhancer of zeste homolog 2 (EZH2) are prevalent in diffuse large B-cell lymphoma (DLBCL) and acute myeloid leukemia (AML). The methyltransferase activity of EZH2 towards lysine 27 on histone H3 (H3K27) and non-histone proteins is dysregulated by the presence of gain-of-function (GOF) and loss-of-function (LOF) mutations altering chromatin compaction, protein complex recruitment, and transcriptional regulation. In this study, a comprehensive multi-omics approach was carried out to characterize the effects of differential H3K27me3 deposition driven by EZH2 mutations. Three stable isogenic mutants (EZH2Y641F, EZH2A677G, and EZH2H689A/F667I) were examined using EpiProfile, H3K27me3 CUT&Tag, ATAC-Seq, transcriptomics, label-free proteomics, and untargeted metabolomics. A discrete set of genes and downstream targets were identified for the EZH2 GOF and LOF mutants that impacted pathways involved in cellular proliferation, differentiation, and migration. Disruption of protein networks and metabolic signatures able to sustain aberrant cell behavior was observed in response to EZH2 mutations. This systems biology-based analysis sheds light on EZH2-mediated cell transformative processes, from the epigenetic to the phenotypic level. These studies provide novel insights into aberrant EZH2 function along with targets that can be explored for improved diagnostics/treatment in hematologic malignancies with mutated EZH2.

Integrated single-cell multiomics uncovers foundational regulatory mechanisms of lens development and pathology

Ocular lens development entails epithelial to fiber cell differentiation, defects in which cause congenital cataract. We report the first single-cell multiomic atlas of lens development, leveraging snRNA-seq, snATAC-seq, and CUT&RUN-seq to discover novel mechanisms of cell fate determination and cataract-linked regulatory networks. A comprehensive profile of cis- and trans-regulatory interactions, including for the cataract-linked transcription factor MAF, is established across a temporal trajectory of fiber cell differentiation. Further, we divulge a conserved epigenetic paradigm of cellular differentiation, defined by progressive loss of H3K27 methylation writer Polycomb repressive complex 2 (PRC2). PRC2 localizes to heterochromatin domains across master-regulator transcription factor gene bodies, suggesting it safeguards epithelial cell fate. Moreover, we demonstrate that FGF hyper-stimulation in vivo leads to MAF network activation and the emergence of novel lens cell states. Collectively, these data depict a comprehensive portrait of lens fiber cell differentiation, while defining regulatory effectors of cell identity and cataract formation.

lncRNAs-EZH2 interaction as promising therapeutic target in cutaneous melanoma

Melanoma is the most lethal skin cancer with increasing incidence worldwide. Despite a great improvement of diagnostics and treatment of melanoma patients, this disease is still a serious clinical problem. Therefore, novel druggable targets are in focus of research. EZH2 is a component of the PRC2 protein complex that mediates epigenetic silencing of target genes. Several mutations activating EZH2 have been identified in melanoma, which contributes to aberrant gene silencing during tumor progression. Emerging evidence indicates that long non-coding RNAs (lncRNAs) are molecular "address codes" for EZH2 silencing specificity, and targeting lncRNAs-EZH2 interaction may slow down the progression of many solid cancers, including melanoma. This review summarizes current knowledge regarding the involvement of lncRNAs in EZH2-mediated gene silencing in melanoma. The possibility of blocking lncRNAs-EZH2 interaction in melanoma as a novel therapeutic option and plausible controversies and drawbacks of this approach are also briefly discussed.

Imagawa-Matsumoto syndrome: SUZ12-related overgrowth disorder

The SUZ12 gene encodes a subunit of polycomb repressive complex 2 (PRC2) that is essential for development by silencing the expression of multiple genes. Germline heterozygous variants in SUZ12 have been found in Imagawa-Matsumoto syndrome (IMMAS) characterized by overgrowth and multiple dysmorphic features. Similarly, both EZH2 and EED also encode a subunit of PRC2 each and their pathogenic variants cause Weaver syndrome and Cohen-Gibson syndrome, respectively. Clinical manifestations of these syndromes significantly overlap, although their different prevalence rates have recently been noted: generalized overgrowth, intellectual disability, scoliosis, and excessive loose skin appear to be less prevalent in IMMAS than in the other two syndromes. We could not determine any apparent genotype-phenotype correlation in IMMAS. The phenotype of neurofibromatosis type 1 arising from NF1 deletion was also shown to be modified by the deletion of SUZ12, 560 kb away. This review deepens our understanding of the clinical and genetic characteristics of IMMAS together with other overgrowth syndromes related to PRC2. We also report on a novel IMMAS patient carrying a splicing variant (c.1023+1G>C) in SUZ12. This patient had a milder phenotype than other previously reported IMMAS cases, with no macrocephaly or overgrowth phenotypes, highlighting the clinical variation in IMMAS.

Transcriptomic Signatures of Single-Suture Craniosynostosis Phenotypes

Craniosynostosis is a birth defect where calvarial sutures close prematurely, as part of a genetic syndrome or independently, with unknown cause. This study aimed to identify differences in gene expression in primary calvarial cell lines derived from patients with four phenotypes of single-suture craniosynostosis, compared to controls. Calvarial bone samples (N = 388 cases/85 controls) were collected from clinical sites during reconstructive skull surgery. Primary cell lines were then derived from the tissue and used for RNA sequencing. Linear models were fit to estimate covariate adjusted associations between gene expression and four phenotypes of single-suture craniosynostosis (lambdoid, metopic, sagittal, and coronal), compared to controls. Sex-stratified analysis was also performed for each phenotype. Differentially expressed genes (DEGs) included 72 genes associated with coronal, 90 genes associated with sagittal, 103 genes associated with metopic, and 33 genes associated with lambdoid craniosynostosis. The sex-stratified analysis revealed more DEGs in males (98) than females (4). There were 16 DEGs that were homeobox (HOX) genes. Three TFs (SUZ12, EZH2, AR) significantly regulated expression of DEGs in one or more phenotypes. Pathway analysis identified four KEGG pathways associated with at least one phenotype of craniosynostosis. Together, this work suggests unique molecular mechanisms related to craniosynostosis phenotype and fetal sex.

EZH2-mediated H3K27me3 is a predictive biomarker and therapeutic target in uveal melanoma

Although gene mutations and aberrant chromosomes are associated with the pathogenesis and prognosis of uveal melanoma (UM), potential therapeutic targets still need to be explored. We aim to determine the predictive value and potential therapeutic target of EZH2 in uveal melanoma. Eighty-five uveal melanoma samples were recruited in our study, including 19 metastatic and 66 nonmetastatic samples. qRT-PCR, immunohistochemistry staining, and western blotting were applied to detect the expression of EZH2 and H3K27me3. We found that EZH2 (41/85, 48.24%) and H3K27me3 (49/85, 57.65%) were overexpressed in uveal melanoma. The expression of EZH2 was not significantly associated with metastasis. High H3K27me3 expression was correlated with poor patient prognosis. UNC 1999, an EZH2 inhibitor, can downregulate H3K27me3 expression and has the most potency to inhibit OMM1 cell growth by the cell cycle and ferroptosis pathway. These results indicate that H3K27me3 can be a biomarker predicting a poor prognosis of UM. EZH2 is the potential therapeutic target for UM.