The time-resolved genomic impact of Wnt/β-catenin signaling

Divergent WNT signaling and drug sensitivity profiles within hepatoblastoma tumors and organoids

Hepatoblastoma, the most prevalent pediatric liver cancer, almost always carries a WNT-activating CTNNB1 mutation, yet exhibits notable molecular heterogeneity. To characterize this heterogeneity and identify novel targeted therapies, we perform comprehensive analysis of hepatoblastomas and tumor-derived organoids using single-cell RNA-seq/ATAC-seq, spatial transcriptomics, and high-throughput drug profiling. We identify two distinct tumor epithelial signatures: hepatic 'fetal' and WNT-high 'embryonal', displaying divergent WNT signaling patterns. The fetal group is enriched for liver-specific WNT targets, while the embryonal group is enriched in canonical WNT target genes. Gene regulatory network analysis reveals enrichment of regulons related to hepatic functions such as bile acid, lipid and xenobiotic metabolism in the fetal subtype but not in the embryonal subtype. In addition, the dichotomous expression pattern of the transcription factors HNF4A and LEF1 allows for a clear distinction between the fetal and embryonal tumor cells. We also perform high-throughput drug screening using patient-derived tumor organoids and identify sensitivity to HDAC inhibitors. Intriguingly, embryonal and fetal tumor organoids are sensitive to FGFR and EGFR inhibitors, respectively, indicating a dependency on EGF/FGF signaling in hepatoblastoma tumorigenesis. In summary, our data uncover the molecular and drug sensitivity landscapes of hepatoblastoma and pave the way for the development of targeted therapies.

Deleterious ZNRF3 germline variants cause neurodevelopmental disorders with mirror brain phenotypes via domain-specific effects on Wnt/β-catenin signaling

Zinc and RING finger 3 (ZNRF3) is a negative-feedback regulator of Wnt/β-catenin signaling, which plays an important role in human brain development. Although somatically frequently mutated in cancer, germline variants in ZNRF3 have not been established as causative for neurodevelopmental disorders (NDDs). We identified 12 individuals with ZNRF3 variants and various phenotypes via GeneMatcher/Decipher and evaluated genotype-phenotype correlation. We performed structural modeling and representative deleterious and control variants were assessed using in vitro transcriptional reporter assays with and without Wnt-ligand Wnt3a and/or Wnt-potentiator R-spondin (RSPO). Eight individuals harbored de novo missense variants and presented with NDD. We found missense variants associated with macrocephalic NDD to cluster in the RING ligase domain. Structural modeling predicted disruption of the ubiquitin ligase function likely compromising Wnt receptor turnover. Accordingly, the functional assays showed enhanced Wnt/β-catenin signaling for these variants in a dominant negative manner. Contrarily, an individual with microcephalic NDD harbored a missense variant in the RSPO-binding domain predicted to disrupt binding affinity to RSPO and showed attenuated Wnt/β-catenin signaling in the same assays. Additionally, four individuals harbored de novo truncating or de novo or inherited large in-frame deletion variants with non-NDD phenotypes, including heart, adrenal, or nephrotic problems. In contrast to NDD-associated missense variants, the effects on Wnt/β-catenin signaling were comparable between the truncating variant and the empty vector and between benign variants and the wild type. In summary, we provide evidence for mirror brain size phenotypes caused by distinct pathomechanisms in Wnt/β-catenin signaling through protein domain-specific deleterious ZNRF3 germline missense variants.

Wnt target gene activation requires β-catenin separation into biomolecular condensates

The Wnt/β-catenin signaling pathway plays numerous essential roles in animal development and tissue/stem cell maintenance. The activation of genes regulated by Wnt/β-catenin signaling requires the nuclear accumulation of β-catenin, a transcriptional co-activator. β-catenin is recruited to many Wnt-regulated enhancers through direct binding to T-cell factor/lymphoid enhancer factor (TCF/LEF) family transcription factors. β-catenin has previously been reported to form phase-separated biomolecular condensates (BMCs), which was implicated as a component of β-catenin's mechanism of action. This function required aromatic amino acid residues in the intrinsically disordered regions (IDRs) at the N- and C-termini of the protein. In this report, we further explore a role for β-catenin BMCs in Wnt target gene regulation. We find that β-catenin BMCs are miscible with LEF1 BMCs in vitro and in cultured cells. We characterized a panel of β-catenin mutants with different combinations of aromatic residue mutations in human cell culture and Drosophila melanogaster. Our data support a model in which aromatic residues across both IDRs contribute to BMC formation and signaling activity. Although different Wnt targets have different sensitivities to loss of β-catenin's aromatic residues, the activation of every target examined was compromised by aromatic substitution. These mutants are not defective in nuclear import or co-immunoprecipitation with several β-catenin binding partners. In addition, residues in the N-terminal IDR with no previously known role in signaling are clearly required for the activation of various Wnt readouts. Consistent with this, deletion of the N-terminal IDR results in a loss of signaling activity, which can be rescued by the addition of heterologous IDRs enriched in aromatic residues. Overall, our work supports a model in which the ability of β-catenin to form biomolecular condensates in the nucleus is tightly linked to its function as a transcriptional co-regulator.

Modulating β-catenin homeostasis for cancer therapy

β-Catenin is a well-established driver of many cancers; however, there are challenges in developing agents targeting β-catenin for clinical use. Recent progress has indicated that most of the pathological changes in β-catenin may be commonly caused by loss of protein homeostasis. Modulation of β-catenin homeostasis, especially by hyperactivation of β-catenin, potentially leads to robust antitumor outcomes. Here, we comprehensively dissect the protein homeostasis of β-catenin in terms of time, compartmentalization, supramolecular assemblies, and dynamics, with emphasis on changes in β-catenin homeostasis upon oncogenic mutations. We propose that altered β-catenin homeostasis could be deleterious for β-catenin-dependent cancers and that modulation of β-catenin homeostasis offers a novel avenue for targeting β-catenin for cancer therapy.

Exhaustive identification of genome-wide binding events of transcriptional regulators

Genome-wide binding assays aspire to map the complete binding pattern of gene regulators. Common practice relies on replication-duplicates or triplicates-and high stringency statistics to favor false negatives over false positives. Here we show that duplicates and triplicates of CUT&RUN are not sufficient to discover the entire activity of transcriptional regulators. We introduce ICEBERG (Increased Capture of Enrichment By Exhaustive Replicate aGgregation), a pipeline that harnesses large numbers of CUT&RUN replicates to discover the full set of binding events and chart the line between false positives and false negatives. We employed ICEBERG to map the full set of H3K4me3-marked regions, the targets of the co-factor β-catenin, and those of the transcription factor TBX3, in human colorectal cancer cells. The ICEBERG datasets allow benchmarking of individual replicates, comparing the performance of peak calling and replication approaches, and expose the arbitrary nature of strategies to identify reproducible peaks. Instead of a static view of genomic targets, ICEBERG establishes a spectrum of detection probabilities across the genome for a given factor, underlying the intrinsic dynamicity of its mechanism of action, and permitting to distinguish frequent from rare regulation events. Finally, ICEBERG discovered instances, undetectable with other approaches, that underlie novel mechanisms of colorectal cancer progression.

Transcriptional responses of cancer cells to heat shock-inducing stimuli involve amplification of robust HSF1 binding

Responses of cells to stimuli are increasingly discovered to involve the binding of sequence-specific transcription factors outside of known target genes. We wanted to determine to what extent the genome-wide binding and function of a transcription factor are shaped by the cell type versus the stimulus. To do so, we induced the Heat Shock Response pathway in two different cancer cell lines with two different stimuli and related the binding of its master regulator HSF1 to nascent RNA and chromatin accessibility. Here, we show that HSF1 binding patterns retain their identity between basal conditions and under different magnitudes of activation, so that common HSF1 binding is globally associated with distinct transcription outcomes. HSF1-induced increase in DNA accessibility was modest in scale, but occurred predominantly at remote genomic sites. Apart from regulating transcription at existing elements including promoters and enhancers, HSF1 binding amplified during responses to stimuli may engage inactive chromatin.

β-catenin repositions over time

How β-catenin, the nuclear activator of the Wnt pathway, affects the chromatin environment of its targets is unknown. Over a time course of stimulation, β-catenin repositions itself around the genome in a cell-type-specific manner, eliciting transient chromatin changes in differentiated cells and progressive shaping of undifferentiated cells.