Histone chaperone HIRA, promyelocytic leukemia protein, and p62/SQSTM1 coordinate to regulate inflammation during cell senescence.

Cellular senescence, a stress-induced stable proliferation arrest associated with an inflammatory senescence-associated secretory phenotype (SASP), is a cause of aging. In senescent cells, cytoplasmic chromatin fragments (CCFs) activate SASP via the anti-viral cGAS/STING pathway. Promyelocytic leukemia (PML) protein organizes PML nuclear bodies (NBs), which are also involved in senescence and anti-viral immunity. The HIRA histone H3.3 chaperone localizes to PML NBs in senescent cells. Here, we show that HIRA and PML are essential for SASP expression, tightly linked to HIRA's localization to PML NBs. Inactivation of HIRA does not directly block expression of nuclear factor κB (NF-κB) target genes. Instead, an H3.3-independent HIRA function activates SASP through a CCF-cGAS-STING-TBK1-NF-κB pathway. HIRA physically interacts with p62/SQSTM1, an autophagy regulator and negative SASP regulator. HIRA and p62 co-localize in PML NBs, linked to their antagonistic regulation of SASP, with PML NBs controlling their spatial configuration. These results outline a role for HIRA and PML in the regulation of SASP.

Hypomorphic and dominant-negative impact of truncated SOX9 dysregulates Hedgehog-Wnt signaling, causing campomelia.

Haploinsufficiency for SOX9, the master chondrogenesis transcription factor, can underlie campomelic dysplasia (CD), an autosomal dominant skeletal malformation syndrome, because heterozygous Sox9 null mice recapitulate the bent limb (campomelia) and some other phenotypes associated with CD. However, in vitro cell assays suggest haploinsufficiency may not apply for certain mutations, notably those that truncate the protein, but in these cases in vivo evidence is lacking and underlying mechanisms are unknown. Here, using conditional mouse mutants, we compared the impact of a heterozygous Sox9 null mutation (Sox9+/-) with the Sox9+/Y440X CD mutation that truncates the C-terminal transactivation domain but spares the DNA-binding domain. While some Sox9+/Y440X mice survived, all Sox9+/- mice died perinatally. However, the skeletal defects were more severe and IHH signaling in developing limb cartilage was significantly enhanced in Sox9+/Y440X compared with Sox9+/-. Activating Sox9Y440X specifically in the chondrocyte-osteoblast lineage caused milder campomelia, and revealed cell- and noncell autonomous mechanisms acting on chondrocyte differentiation and osteogenesis in the perichondrium. Transcriptome analyses of developing Sox9+/Y440X limbs revealed dysregulated expression of genes for the extracellular matrix, as well as changes consistent with aberrant WNT and HH signaling. SOX9Y440X failed to interact with ß-catenin and was unable to suppress transactivation of Ihh in cell-based assays. We propose enhanced HH signaling in the adjacent perichondrium induces asymmetrically localized excessive perichondrial osteogenesis resulting in campomelia. Our study implicates combined haploinsufficiency/hypomorphic and dominant-negative actions of SOX9Y440X, cell-autonomous and noncell autonomous mechanisms, and dysregulated WNT and HH signaling, as the cause of human campomelia.

Accurate identification of structural variations from cancer samples.

Structural variations (SVs) are commonly found in cancer genomes. They can cause gene amplification, deletion and fusion, among other functional consequences. With an average read length of hundreds of kilobases, nano-channel-based optical DNA mapping is powerful in detecting large SVs. However, existing SV calling methods are not tailored for cancer samples, which have special properties such as mixed cell types and sub-clones. Here we propose the Cancer Optical Mapping for detecting Structural Variations (COMSV) method that is specifically designed for cancer samples. It shows high sensitivity and specificity in benchmark comparisons. Applying to cancer cell lines and patient samples, COMSV identifies hundreds of novel SVs per sample.

Extensive in vivo metabolite-protein interactions revealed by large-scale systematic analyses.

Natural small compounds comprise most cellular molecules and bind proteins as substrates, products, cofactors, and ligands. However, a large-scale investigation of in vivo protein-small metabolite interactions has not been performed. We developed a mass spectrometry assay for the large-scale identification of in vivo protein-hydrophobic small metabolite interactions in yeast and analyzed compounds that bind ergosterol biosynthetic proteins and protein kinases. Many of these proteins bind small metabolites; a few interactions were previously known, but the vast majority are new. Importantly, many key regulatory proteins such as protein kinases bind metabolites. Ergosterol was found to bind many proteins and may function as a general regulator. It is required for the activity of Ypk1, a mammalian AKT/SGK kinase homolog. Our study defines potential key regulatory steps in lipid biosynthetic pathways and suggests that small metabolites may play a more general role as regulators of protein activity and function than previously appreciated.

Quantifying full-length circular RNAs in cancer.

Circular RNAs (circRNAs) are abundantly expressed in cancer. Their resistance to exonucleases enables them to have potentially stable interactions with different types of biomolecules. Alternative splicing can create different circRNA isoforms that have different sequences and unequal interaction potentials. The study of circRNA function thus requires knowledge of complete circRNA sequences. Here we describe psirc, a method that can identify full-length circRNA isoforms and quantify their expression levels from RNA sequencing data. We confirm the effectiveness and computational efficiency of psirc using both simulated and actual experimental data. Applying psirc on transcriptome profiles from nasopharyngeal carcinoma and normal nasopharynx samples, we discover and validate circRNA isoforms differentially expressed between the two groups. Compared with the assumed circular isoforms derived from linear transcript annotations, some of the alternatively spliced circular isoforms have 100 times higher expression and contain substantially fewer microRNA response elements, showing the importance of quantifying full-length circRNA isoforms.

Pipeline Therapies for Gout.

PURPOSE OF REVIEW: Despite effective available treatments, gout management is often unsuccessful in getting patients to target serum urate goal and in managing flares in the setting of comorbidities. Studies addressing future treatment options for short- and long-term management are reviewed. RECENT FINDINGS: URAT-1 blocking agents have been helpful but have had limitations related to effects on renal function, lack of efficacy with renal impairment, and potential to increase renal stones. Dotinurad may function in the setting of decreased renal function. Arhalofenate has anti-URAT-1 activity and may also blunt gout flares. A new xanthine oxidase inhibitor (XOI), tigulixostat, is under study. New uricase treatments manufactured in combination with agents that can reduce immunogenicity may make uricase treatment simpler. A unique strategy of inhibiting gut uricase may offer the benefits of avoiding systemic absorption. For gout flares, IL-1ß inhibitor studies in progress include different dosing schedules. Dapansutrile, an oral agent under investigation, inhibits activation of the NLRP3 inflammasome and may be an effective anti-inflammatory. New treatments for gout that are under study may work in the setting of comorbidities, simplify management, utilize new mechanisms, or have reduced side effects.

Targeting PPAR-gamma counteracts tumour adaptation to immune-checkpoint blockade in hepatocellular carcinoma.

OBJECTIVE: Therapy-induced tumour microenvironment (TME) remodelling poses a major hurdle for cancer cure. As the majority of patients with hepatocellular carcinoma (HCC) exhibits primary or acquired resistance to antiprogrammed cell death (ligand)-1 (anti-PD-[L]1) therapies, we aimed to investigate the mechanisms underlying tumour adaptation to immune-checkpoint targeting. DESIGN: Two immunotherapy-resistant HCC models were generated by serial orthotopic implantation of HCC cells through anti-PD-L1-treated syngeneic, immunocompetent mice and interrogated by single-cell RNA sequencing (scRNA-seq), genomic and immune profiling. Key signalling pathway was investigated by lentiviral-mediated knockdown and pharmacological inhibition, and further verified by scRNA-seq analysis of HCC tumour biopsies from a phase II trial of pembrolizumab (NCT03419481). RESULTS: Anti-PD-L1-resistant tumours grew >10-fold larger than parental tumours in immunocompetent but not immunocompromised mice without overt genetic changes, which were accompanied by intratumoral accumulation of myeloid-derived suppressor cells (MDSC), cytotoxic to exhausted CD8+ T cell conversion and exclusion. Mechanistically, tumour cell-intrinsic upregulation of peroxisome proliferator-activated receptor-gamma (PPARγ) transcriptionally activated vascular endothelial growth factor-A (VEGF-A) production to drive MDSC expansion and CD8+ T cell dysfunction. A selective PPARγ antagonist triggered an immune suppressive-to-stimulatory TME conversion and resensitised tumours to anti-PD-L1 therapy in orthotopic and spontaneous HCC models. Importantly, 40% (6/15) of patients with HCC resistant to pembrolizumab exhibited tumorous PPARγ induction. Moreover, higher baseline PPARγ expression was associated with poorer survival of anti-PD-(L)1-treated patients in multiple cancer types. CONCLUSION: We uncover an adaptive transcriptional programme by which tumour cells evade immune-checkpoint targeting via PPARγ/VEGF-A-mediated TME immunosuppression, thus providing a strategy for counteracting immunotherapeutic resistance in HCC.

Whole-genome analysis of noncoding genetic variations identifies multiscale regulatory element perturbations associated with Hirschsprung disease.

It is widely recognized that noncoding genetic variants play important roles in many human diseases, but there are multiple challenges that hinder the identification of functional disease-associated noncoding variants. The number of noncoding variants can be many times that of coding variants; many of them are not functional but in linkage disequilibrium with the functional ones; different variants can have epistatic effects; different variants can affect the same genes or pathways in different individuals; and some variants are related to each other not by affecting the same gene but by affecting the binding of the same upstream regulator. To overcome these difficulties, we propose a novel analysis framework that considers convergent impacts of different genetic variants on protein binding, which provides multiscale information about disease-associated perturbations of regulatory elements, genes, and pathways. Applying it to our whole-genome sequencing data of 918 short-segment Hirschsprung disease patients and matched controls, we identify various novel genes not detected by standard single-variant and region-based tests, functionally centering on neural crest migration and development. Our framework also identifies upstream regulators whose binding is influenced by the noncoding variants. Using human neural crest cells, we confirm cell stage-specific regulatory roles of three top novel regulatory elements on our list, respectively in the RET, RASGEF1A, and PIK3C2B loci. In the PIK3C2B regulatory element, we further show that a noncoding variant found only in the patients affects the binding of the gliogenesis regulator NFIA, with a corresponding up-regulation of multiple genes in the same topologically associating domain.

Supervised enhancer prediction with epigenetic pattern recognition and targeted validation.

Enhancers are important non-coding elements, but they have traditionally been hard to characterize experimentally. The development of massively parallel assays allows the characterization of large numbers of enhancers for the first time. Here, we developed a framework using Drosophila STARR-seq to create shape-matching filters based on meta-profiles of epigenetic features. We integrated these features with supervised machine-learning algorithms to predict enhancers. We further demonstrated that our model could be transferred to predict enhancers in mammals. We comprehensively validated the predictions using a combination of in vivo and in vitro approaches, involving transgenic assays in mice and transduction-based reporter assays in human cell lines (153 enhancers in total). The results confirmed that our model can accurately predict enhancers in different species without re-parameterization. Finally, we examined the transcription factor binding patterns at predicted enhancers versus promoters. We demonstrated that these patterns enable the construction of a secondary model that effectively distinguishes enhancers and promoters.

Deep learning identifies and quantifies recombination hotspot determinants.

MOTIVATION: Recombination is one of the essential genetic processes for sexually reproducing organisms, which can happen more frequently in some regions, called recombination hotspots. Although several factors, such as PRDM9 binding motifs, are known to be related to the hotspots, their contributions to the recombination hotspots have not been quantified, and other determinants are yet to be elucidated. Here, we propose a computational method, RHSNet, based on deep learning and signal processing, to identify and quantify the hotspot determinants in a purely data-driven manner, utilizing datasets from various studies, populations, sexes and species. RESULTS: RHSNet can significantly outperform other sequence-based methods on multiple datasets across different species, sexes and studies. In addition to being able to identify hotspot regions and the well-known determinants accurately, more importantly, RHSNet can quantify the determinants that contribute significantly to the recombination hotspot formation in the relation between PRDM9 binding motif, histone modification and GC content. Further cross-sex, cross-population and cross-species studies suggest that the proposed method has the generalization power and potential to identify and quantify the evolutionary determinant motifs. AVAILABILITY AND IMPLEMENTATION: https://github.com/frankchen121212/RHSNet. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Measuring the impact of steroid therapy on health-related quality of life in patients with rheumatic diseases: international development of a glucocorticoid treatment-specific patient-reported outcome measure.

OBJECTIVES: Glucocorticoids (GCs) ('steroids') are used to treat rheumatic diseases but adverse effects are common. We aimed to explore the impact of GC therapy on health-related quality of life (HRQoL), to inform the development of a treatment-specific patient-reported outcome measure (PROM) for use in clinical trials and practice. METHODS: Semi-structured qualitative interviews were conducted with patients from the UK, USA and Australia, treated for a rheumatic condition with GCs in the last 2 years. Purposive sampling was used to select participants with a range of demographic and disease features. An initial conceptual framework informed interview prompts and cues. Interviews elicited GC-related physical and psychological symptoms and salient aspects of HRQoL in relation to GC therapy. Interview data were analysed inductively to develop initial individual themes and domains. Candidate questionnaire items were developed and refined. RESULTS: Sixty semi-structured qualitative interviews were conducted (UK n = 34, USA n = 10, Australia n = 16). The mean age was 58 years; 39/60 were female; and 18 rheumatic diseases were represented. Some 126 individual themes were identified and organized into six domains: physical symptoms; psychological symptoms; psychological impact of steroids; impact of steroids on participation; impact of steroids on relationships; and benefits of steroids. Candidate questionnaire items were tested and refined by piloting with patient research partners, iterative rounds of cognitive interviews and linguistic translatability assessment, informing a draft questionnaire. CONCLUSION: We describe an international qualitative study to develop candidate items for a treatment-specific PROM for patients with rheumatic diseases. A future survey will enable the validation of a final version of the PROM.

2023 American College of Rheumatology and American Association of Hip and Knee Surgeons Clinical Practice Guideline for the Optimal Timing of Elective Hip or Knee Arthroplasty for Patients With Symptomatic Moderate-to-Severe Osteoarthritis or Advanced Symptomatic Osteonecrosis With Secondary Arthritis for Whom Nonoperative Therapy Is Ineffective.

OBJECTIVE: To develop evidence-based consensus recommendations for the optimal timing of hip and knee arthroplasty to improve patient-important outcomes including, but not limited to, pain, function, infection, hospitalization, and death at 1 year for patients with symptomatic and radiographic moderate-to-severe osteoarthritis or advanced symptomatic osteonecrosis with secondary arthritis of the hip or knee who have previously attempted nonoperative therapy, and for whom nonoperative therapy was ineffective, and who have chosen to undergo elective hip or knee arthroplasty (collectively referred to as TJA). METHODS: We developed 13 clinically relevant population, intervention, comparator, outcomes (PICO) questions. After a systematic literature review, the Grading of Recommendations Assessment, Development and Evaluation (GRADE) approach was used to rate the quality of evidence (high, moderate, low, or very low), and evidence tables were created. A Voting Panel, including 13 physicians and patients, discussed the PICO questions until consensus was achieved on the direction (for/against) and strength (strong/conditional) of the recommendations. RESULTS: The panel conditionally recommended against delaying TJA to pursue additional nonoperative treatment including physical therapy, nonsteroidal antiinflammatory drugs, ambulatory aids, and intraarticular injections. It conditionally recommended delaying TJA for nicotine reduction or cessation. The panel conditionally recommended delay for better glycemic control for patients who have diabetes mellitus, although no specific measure or level was identified. There was consensus that obesity by itself was not a reason for delay, but that weight loss should be strongly encouraged, and the increase in operative risk should be discussed. The panel conditionally recommended against delay in patients who have severe deformity or bone loss, or in patients who have a neuropathic joint. Evidence for all recommendations was graded as low or very low quality. CONCLUSION: This guideline provides evidence-based recommendations regarding the optimal timing of TJA in patients who have symptomatic and radiographic moderate-to-severe osteoarthritis or advanced symptomatic osteonecrosis with secondary arthritis for whom nonoperative therapy was ineffective to improve patient-important outcomes, including pain, function, infection, hospitalization, and death at 1 year. We acknowledge that the evidence is of low quality primarily due to indirectness and hope future research will allow for further refinement of the recommendations.

Accurate reconstruction of viral genomes in human cells from short reads using iterative refinement.

BACKGROUND: After an infection, human cells may contain viral genomes in the form of episomes or integrated DNA. Comparing the genomic sequences of different strains of a virus in human cells can often provide useful insights into its behaviour, activity and pathology, and may help develop methods for disease prevention and treatment. To support such comparative analyses, the viral genomes need to be accurately reconstructed from a large number of samples. Previous efforts either rely on customized experimental protocols or require high similarity between the sequenced genomes and a reference, both of which limit the general applicability of these approaches. In this study, we propose a pipeline, named ASPIRE, for reconstructing viral genomes accurately from short reads data of human samples, which are increasingly available from genome projects and personal genomics. ASPIRE contains a basic part that involves de novo assembly, tiling and gap filling, and additional components for iterative refinement, sequence corrections and wrapping. RESULTS: Evaluated by the alignment quality of sequencing reads to the reconstructed genomes, these additional components improve the assembly quality in general, and in some particular samples quite substantially, especially when the sequenced genome is significantly different from the reference. We use ASPIRE to reconstruct the genomes of Epstein Barr Virus (EBV) from the whole-genome sequencing data of 61 nasopharyngeal carcinoma (NPC) samples and provide these sequences as a resource for EBV research. CONCLUSIONS: ASPIRE improves the quality of the reconstructed EBV genomes in published studies and outperforms TRACESPipe in some samples considered.

New guidelines for DNA methylome studies regarding 5-hydroxymethylcytosine for understanding transcriptional regulation.

Many DNA methylome profiling methods cannot distinguish between 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC). Because 5mC typically acts as a repressive mark whereas 5hmC is an intermediate form during active demethylation, the inability to separate their signals could lead to incorrect interpretation of the data. Is the extra information contained in 5hmC signals worth the additional experimental and computational costs? Here we combine whole-genome bisulfite sequencing (WGBS) and oxidative WGBS (oxWGBS) data in various human tissues to investigate the quantitative relationships between gene expression and the two forms of DNA methylation at promoters, transcript bodies, and immediate downstream regions. We find that 5mC and 5hmC signals correlate with gene expression in the same direction in most samples. Considering both types of signals increases the accuracy of expression levels inferred from methylation data by a median of 18.2% as compared to having only WGBS data, showing that the two forms of methylation provide complementary information about gene expression. Differential analysis between matched tumor and normal pairs is particularly affected by the superposition of 5mC and 5hmC signals in WGBS data, with at least 25%-40% of the differentially methylated regions (DMRs) identified from 5mC signals not detected from WGBS data. Our results also confirm a previous finding that methylation signals at transcript bodies are more indicative of gene expression levels than promoter methylation signals. Overall, our study provides data for evaluating the cost-effectiveness of some experimental and analysis options in the study of DNA methylation in normal and cancer samples.

Shaping the nebulous enhancer in the era of high-throughput assays and genome editing.

Since the 1st discovery of transcriptional enhancers in 1981, their textbook definition has remained largely unchanged in the past 37 years. With the emergence of high-throughput assays and genome editing, which are switching the paradigm from bottom-up discovery and testing of individual enhancers to top-down profiling of enhancer activities genome-wide, it has become increasingly evidenced that this classical definition has left substantial gray areas in different aspects. Here we survey a representative set of recent research articles and report the definitions of enhancers they have adopted. The results reveal that a wide spectrum of definitions is used usually without the definition stated explicitly, which could lead to difficulties in data interpretation and downstream analyses. Based on these findings, we discuss the practical implications and suggestions for future studies.

Identification of a major determinant for serine-threonine kinase phosphoacceptor specificity.

Eukaryotic protein kinases are generally classified as being either tyrosine or serine-threonine specific. Though not evident from inspection of their primary sequences, many serine-threonine kinases display a significant preference for serine or threonine as the phosphoacceptor residue. Here we show that a residue located in the kinase activation segment, which we term the "DFG+1" residue, acts as a major determinant for serine-threonine phosphorylation site specificity. Mutation of this residue was sufficient to switch the phosphorylation site preference for multiple kinases, including the serine-specific kinase PAK4 and the threonine-specific kinase MST4. Kinetic analysis of peptide substrate phosphorylation and crystal structures of PAK4-peptide complexes suggested that phosphoacceptor residue preference is not mediated by stronger binding of the favored substrate. Rather, favored kinase-phosphoacceptor combinations likely promote a conformation optimal for catalysis. Understanding the rules governing kinase phosphoacceptor preference allows kinases to be classified as serine or threonine specific based on their sequence.

Racial, ethnic, and healthcare disparities in rheumatoid arthritis.

PURPOSE OF REVIEW: This review highlights the available data describing racial and ethnic health disparities among patients with rheumatoid arthritis in the United States from an epidemiological, disease activity, and wider socioeconomic standpoint. RECENT FINDINGS: Despite centralized government initiatives to include more underrepresentative minority populations into research, many of the studies that examined rheumatoid arthritis still fail to include sizeable cohorts of races or ethnic groups other than whites. Evidence is slowly mounting that individual, provider, and system-level barriers exist and contribute to unequal care that leads to poorer outcomes amongst patients with rheumatoid arthritis. As rheumatoid arthritis is a progressive disease, early treatment is crucial to delay functional decline - a narrow window for many minority patients who are disproportionality affected by disability. SUMMARY: To combat the inequality that exists amongst rheumatoid arthritis patients we must focus on why discrepancies exist on every level, system, physician, patient, and illness. Further research is needed to tease the complex interplay between race, social economic status, medical access, and outcomes to explain the disparities found in rheumatoid arthritis.

A general near-exact k-mer counting method with low memory consumption enables de novo assembly of 106× human sequence data in 2.7 hours.

MOTIVATION: In de novo sequence assembly, a standard pre-processing step is k-mer counting, which computes the number of occurrences of every length-k sub-sequence in the sequencing reads. Sequencing errors can produce many k-mers that do not appear in the genome, leading to the need for an excessive amount of memory during counting. This issue is particularly serious when the genome to be assembled is large, the sequencing depth is high, or when the memory available is limited. RESULTS: Here, we propose a fast near-exact k-mer counting method, CQF-deNoise, which has a module for dynamically removing noisy false k-mers. It automatically determines the suitable time and number of rounds of noise removal according to a user-specified wrong removal rate. We tested CQF-deNoise comprehensively using data generated from a diverse set of genomes with various data properties, and found that the memory consumed was almost constant regardless of the sequencing errors while the noise removal procedure had minimal effects on counting accuracy. Compared with four state-of-the-art k-mer counting methods, CQF-deNoise consistently performed the best in terms of memory usage, consuming 49-76% less memory than the second best method. When counting the k-mers from a human dataset with around 60× coverage, the peak memory usage of CQF-deNoise was only 10.9 GB (gigabytes) for k = 28 and 21.5 GB for k = 55. De novo assembly of 106× human sequencing data using CQF-deNoise for k-mer counting required only 2.7 h and 90 GB peak memory. AVAILABILITY AND IMPLEMENTATION: The source codes of CQF-deNoise and SH-assembly are available at https://github.com/Christina-hshi/CQF-deNoise.git and https://github.com/Christina-hshi/SH-assembly.git, respectively, both under the BSD 3-Clause license.

Asymptomatic hyperuricemia: is it really asymptomatic?

PURPOSE OF REVIEW: Hyperuricemia is highly prevalent, affecting approximately 38 million individuals in the United States. However, the significance of asymptomatic hyperuricemia - hyperuricemia in the absence of gout - continues to be debated. RECENT FINDINGS: Asymptomatic hyperuricemia results in monosodium urate crystal deposition in tissues, which may promote chronic inflammation. Intracellularly, hyperuricemia inhibits the master regulator adenosine monophosphate (AMP)-associated protein kinase and may condition innate immune responses through durable epigenetic modifications. At the population level, asymptomatic hyperuricemia is associated with multiple comorbidities, including hypertension, chronic kidney disease, coronary artery disease, and diabetes; limitations of these studies include that most are retrospective and some do not rigorously distinguish between asymptomatic hyperuricemia and gout. Treatment studies suggest that urate lowering may reduce the risk of incidence or progression of some of these comorbidities; unfortunately, many of these treatment studies are small or flawed, and not all study results are consistent. SUMMARY: Accumulating evidence suggests that asymptomatic hyperuricemia contributes to the comorbidities with which it associates and that proper asymptomatic hyperuricemia treatment may reduce future risk. Additional prospective trials are needed to definitely establish causality and support decision-making as to whether, and which patients with asymptomatic hyperuricemia would warrant urate-lowering treatment.

Loss of tumor suppressor IGFBP4 drives epigenetic reprogramming in hepatic carcinogenesis.

Genomic sequencing of hepatocellular carcinoma (HCC) uncovers a paucity of actionable mutations, underscoring the necessity to exploit epigenetic vulnerabilities for therapeutics. In HCC, EZH2-mediated H3K27me3 represents a major oncogenic chromatin modification, but how it modulates the therapeutic vulnerability of signaling pathways remains unknown. Here, we show EZH2 acts antagonistically to AKT signaling in maintaining H3K27 methylome through epigenetic silencing of IGFBP4. ChIP-seq revealed enrichment of Ezh2/H3K27me3 at silenced loci in HBx-transgenic mouse-derived HCCs, including Igfbp4 whose down-regulation significantly correlated with EZH2 overexpression and poor survivals of HCC patients. Functional characterizations demonstrated potent growth- and invasion-suppressive functions of IGFBP4, which was associated with transcriptomic alterations leading to deregulation of multiple signaling pathways. Mechanistically, IGFBP4 stimulated AKT/EZH2 phosphorylation to abrogate H3K27me3-mediated silencing, forming a reciprocal feedback loop that suppressed core transcription factor networks (FOXA1/HNF1A/HNF4A/KLF9/NR1H4) for normal liver homeostasis. Consequently, the in vivo tumorigenicity of IGFBP4-silenced HCC cells was vulnerable to pharmacological inhibition of EZH2, but not AKT. Our study unveils chromatin regulation of a novel liver tumor suppressor IGFBP4, which constitutes an AKT-EZH2 reciprocal loop in driving H3K27me3-mediated epigenetic reprogramming. Defining the aberrant chromatin landscape of HCC sheds light into the mechanistic basis of effective EZH2-targeted inhibition.