A β-catenin-driven switch in TCF/LEF transcription factor binding to DNA target sites promotes commitment of mammalian nephron progenitor cells

  1. Qiuyu Guo
  2. Albert Kim
  3. Bin Li
  4. Andrew Ransick
  5. Helena Bugacov
  6. Xi Chen
  7. Nils Lindström
  8. Aaron Brown
  9. Leif Oxburgh
  10. Bing Ren
  11. Andrew P McMahon  Is a corresponding author
  1. Department of Stem Cell Biology and Regenerative Medicine, Eli and Edythe Broad-CIRM Center for Regenerative Medicine and Stem Cell Research, Keck School of Medicine of the University of Southern California, United States
  2. The Rogosin Institute, United States
  3. Center for Molecular Medicine, Maine Medical Center Research Institute, United States
  4. Ludwig Institute for Cancer Research, Department of Cellular and Molecular Medicine, Institute of Genomic Medicine, Moores Cancer Center, University of California San Diego, United States
5 figures, 1 table and 7 additional files

Figures

Figure 1 with 1 supplement
Nephron progenitor expansion medium (NPEM) supplemented with differential levels of CHIR99021 (CHIR) models nephron progenitor cell maintenance or differentiation in a plate.

(A) Immunofluorescence (IF) staining showing expression level of Six2, non-phospho (NP) β-catenin and Jag1 in nephron progenitor cell (NPC) cultured in NPEM supplemented with various CHIR dosages. (B) Relative intensity of IF signals from individual cells in experiment associated with A. (C) Heatmap/hierarchical cluster of expression levels of NPC signature, self-renewal, and differentiation marker genes. (D) Top five enriched Gene Ontology (GO) terms of indicated differentially expressed gene lists, analyzed by DAVID. Link to high-definition figure: https://www.dropbox.com/s/76chyrs07m0toen/Fig%201.pdf?dl=0.

Figure 1—figure supplement 1
Supplementary RNA-Seq data analysis.

(A) Overview of experiment design and data available (gray). (B) Bar plots show RT-qPCR measurement of relative expression of the indicated genes, as verification of results in Figure 1C. (C) Hierarchical cluster of R-square values between transcriptome-wide TPM of the indicated pair of replicate RNA-Seq data sets. (D) Top five enriched Gene Ontology (GO) terms of genes differentially expressed between low CHIR99021 (CHIR) condition and uncultured nephron progenitor cell (NPC). (E) Top five enriched GO terms of genes differentially expressed between low CHIR and no CHIR conditions. (F) Venn diagram shows overlap of β-catenin target genes with all genes associated to the GO term ‘cell cycle’ that are highly expressed in low CHIR versus no CHIR condition. (G) Heatmap and hierarchical cluster showing log2 relative TPM (TPM divided by mean across samples) to reflect change of gene expression (at isoform level) of Karner et al., 2011) class II genes (differentially expressed low CHIR > no CHIR or high CHIR > low CHIR) in our data set. The highlighted gene is supported by genetic evidence of regulation by β-catenin in the current data. (H) Immunofluorescence (IF) staining of SIX2 in unpurified NPC-free cortex (NFC) cells and NPC freshly purified from NFC; bar plots showing percentage of SIX2+ cells in the preps.

Figure 2 with 1 supplement
High dosage of CHIR99021 (CHIR) triggered change of nephron progenitor cell (NPC) epigenome.

(A, B) Genome browser view of RNA-Seq, ATAC-Seq, as well as Six2, H3K27ac, H3K27me3, and Ser5P ChIP-Seq data near Six2 (A) and Wnt4 (B) in low CHIR (left) and high CHIR (right) conditions. Black arrow indicates Six2DE and Wnt4DE, respectively. (C) Display of differentially accessible regions (DARs) generated by the indicated comparisons. (Left) Heatmaps showing log2 normalized read counts of top 500 most significant DARs; (middle) top five most significant Gene Ontology (GO) terms associated with the DARs; (right) top five most enriched motifs discovered de novo in the DARs; * indicates less well-conserved motifs for the factor. NFC: NPC-free cortex.

Figure 2—figure supplement 1
Supplementary ATAC-Seq data analysis.

(A) Hierarchical cluster of R-square values between normalized ATAC-Seq reads within merged peaks from all samples. (B) Histograms of distances from differentially accessible (DA) chromatin regions to transcriptional start site (TSS) implicate a predominant enhancer feature. CHIR: CHIR99021; NPC: nephron progenitor cell; NFC: NPC-free cortex.

Figure 3 with 2 supplements
Differential expression of TCF/LEF family transcription factors in nephron progenitor cell (NPC) in response to distinct level of CHIR99021 (CHIR).

(A) Bar plots showing RNA-Seq measured expression levels of TCF/LEF family factors in NPC cultured in nephron progenitor expansion medium (NPEM) culture supplemented with various concentrations of CHIR. (B) Immunofluorescence (IF) staining of TCF/LEF family factors in NPEM cultured with conditions indicated. (C) Relative intensity of IF signals from individual cells in experiment associated with B. Link to high-definition figure: https://www.dropbox.com/s/pvhu1ffhoujt39p/Fig%203.pdf?dl=0. NFC: NPC-free cortex.

Figure 3—figure supplement 1
Expression of TCF/LEF factors measured by immunoblots.

(A) Immunoblots of TCF/LEF family factors in nephron progenitor cell cultured in nephron progenitor expansion medium with the indicated treatment; red arrows indicate target band. (B) Quantification of target band intensity normalized to histone H3 intensity. CHIR: CHIR99021.

Figure 3—figure supplement 2
Supplementary evidence for differential expression of TCF/LEF factors.

(A) Bar plots show RT-qPCR measurement of relative expression of TCF/LEF family factors, as verification for results in Figure 3A. (B) Bar plots show expression of individual transcripts of TCF/LEF factors in our RNA-Seq data. (C) tSNE plot displaying unbiased cluster of nephron lineage cells profiled by single-cell RNA-Seq. (D) Feature plots displaying distribution of self-renewal (red) and differentiation (green) marker genes transcripts on the tSNE plot. (E) Feature plots showing distribution of TCF/LEF factors transcripts on the tSNE plot. (F) Dot plots showing accumulated expression level of marker genes as well as TCF/LEF factor in selected clusters of cells. NPC: nephron progenitor cell; NFC: NPC-free cortex; CHIR: CHIR99021.

Figure 4 with 3 supplements
Increased CHIR99021 (CHIR) dosage induces a switch of TCF/LEF factors binding to the genome.

(A, B) Genome browser view of Wnt4 enhancer locus showing ChIP-Seq signal of TCF/LEF factors in low and high CHIR conditions. (C) Histograms showing binding intensity of TCF/LEF factors and chromatin markers on the three sets of TCF/LEF binding sites assigned by overlap of low CHIR Tcf7l1 and high CHIR Lef1 binding sites. (D) Result from de novo motif discovery of the three sets of TCF/LEF binding sites described in A. (E) Top Gene Ontology (GO) terms associated with the corresponding sets of TCF/LEF binding sites shown in B. (F) Percentage of TCF/LEF target genes belonging to different sets in differential expressed genes specific to low CHIR or high CHIR conditions as described in Figure 1. (G) Venn diagram showing overlap of set 2 and 3 target genes assigned by GREAT (McLean et al., 2010). (H) Histograms showing quantification of reads from the indicated data sets in ±2 kb of β-catenin binding sites in high CHIR condition.

Figure 4—figure supplement 1
Supplementary ChIP-Seq data analysis.

(A) Numbers of and overlap between binding sites of the same factors between different conditions. (B) Numbers of and overlap between binding sites of different factors in the same conditions. (C) Most enriched motifs by de novo discovery (Homer) from binding sites of the data sets indicated. (D) Hierarchical clustering of normalized read counts of ChIP-Seq data sets on merged binding sites. NPC: nephron progenitor cell; CHIR: CHIR99021.

Figure 4—figure supplement 2
Analysis of Tcf7l1 binding in low CHIR99021 (CHIR).

(A) Overlap of direct and indirect Tcf7l1 binding in low CHIR with Lef1 binding in high CHIR (left) and enrichment of ChIP-Seq reads on the two types of binding events. (B) Distribution of ATAC-Seq, H3K27ac ChIP-Seq, and Ser5P ChIP-Seq signals on Tcf7l1 binding sites in low CHIR with or without TCF/LEF motifs. (C) Gene Ontology (GO) analysis by GREAT on Tcf7l1 binding sites in low CHIR with or without TCF/LEF motifs. (D) Overlap of predicted target genes between Tcf7l1 binding sites in low CHIR with or without TCF/LEF motifs. (E) Histograms showing intensities of ATAC-Seq or ChIP-Seq signals around set 1, 2, and 3 TCF/LEF binding sites in low CHIR and high CHIR conditions.

Figure 4—figure supplement 3
β-Catenin binding sites near (A) Pla2g7 and (B) Tafa5, two β-catenin target genes reported in Karner et al., 2011.

The ChIP-qPCR target sites were marked as red bins on the top track. The highlighted region is where our data is consistent with Karner’s.

Figure 5 with 1 supplement
β-Catenin activates gene expression through both stable and de novo enhancer–promoter loops.

(A) β-Catenin binding sites that overlap with an anchor of chromatin loop in high CHIR99021 (CHIR), the proportion that connects to a transcriptional start site (TSS) (gray in the pie chart) and segregation between two types of loops defined in B. (B) Overlap of chromatin loops between low CHIR and high CHIR conditions. Examples (C, E) and schematics (D, F) of β-catenin utilizing low/high CHIR-shared enhancer–promoter loops to activate Wnt4 (C, D) or high CHIR-specific loops to activate Lef1 (E, F). Black arrow at the bottom indicates the β-catenin binding sites involved in the loops.

Figure 5—figure supplement 1
Supplementary Hi-C data analysis.

(A) Reproducible loops between replicates. (B) Bar plots showing the percentage of loops in the corresponding category (high CHIR99021 [CHIR] only, low CHIR only or share) that are linked to differentially expressed genes (DEG) that are high in low CHIR or high CHIR.

Tables

Key resources table
Reagent type
(species) or resource
DesignationSource or referenceIdentifiersAdditional information
AntibodyRabbit polyclonal anti-Six2ProteinTech11562-1-AP; RRID:AB_2189084ChIP (1:200)
AntibodyMouse monoclonal IgG1 anti-Six2AbnovaH00010736-M01; RRID:AB_436993IF (1:1000)
AntibodyRabbit polyclonal anti-β-cateninThermoFisher71–2700; RRID:AB_2533982ChIP (1:200)
AntibodyRabbit monoclonal anti-non-phopho-β-cateninCell Signaling Technology8814; RRID:AB_11127203ChIP (1:200), IF (1:500) and western blot (1:1000)
AntibodyMouse monoclonal anti-Tcf7l1Santa Cruzsc-166411; RRID:AB_2302942IF (1:1000), western blot (1:1000)
AntibodyRabbit polyclonal anti-Tcf7l1Thermo ScientificPA5-40327; RRID:AB_2577173ChIP (1:200)
AntibodyRabbit monoclonal anti-Tcf7l2Cell Signaling Technology2569; RRID:AB_2199816ChIP (1:200), IF (1:500), western blot (1:1000)
AntibodyRabbit monoclonal anti-Tcf7Cell Signaling Technology2203; RRID:AB_2199302ChIP (1:200), western blot (1:1000)
AntibodyRat monoclonal IgG2b anti-Tcf7R and D SystemsMAB8224IF (1:100)
AntibodyRabbit monoclonal anti-Lef1Cell Signaling Technology2230; RRID:AB_823558ChIP (1:200), IF (1:500), western blot (1:1000)
AntibodyRabbit monoclonal anti-histone H3Abcamab1791; RRID:AB_302613Western blot (1:1000)
AntibodyRabbit monoclonal anti-histone H3K27acAbcamab4729; RRID:AB_2118291ChIP (1:200)
AntibodyMouse monoclonal anti-histone H3K27me2me3Active Motif39536; RRID:AB_2793247ChIP (1:200)
AntibodyMouse monoclonal anti-Ser5P-RNAPIIMillipore05–623; RRID:AB_309852ChIP (1:200)
AntibodyGoat polyclonal anti-Jag1R and D SystemsAF599; RRID:AB_2128257IF (1:50)
AntibodyHRP-conjugated goat anti-mouse IgG1Thermo ScientificA10551Western blot (1:1000)
AntibodyHRP-conjugated goat anti-rabbit IgGCell Signaling Technology7074Western blot (1:1000)
AntibodyAlexa647 mouse anti-rat IgG2bAbcamab172335IF (1:500)
AntibodyAlexa555 donkey anti-rabbit IgGAbcamab150074IF (1:500)
AntibodyAlexa488 goat anti-mouse IgG1Thermo ScientificA10551IF (1:500)
Commercial assay or kitSuperScript IV VILO Master Mix with ezDNase EnzymeThermo Fisher Scientific11766050
Commercial assay or kitLuna Universal qPCR Master MixNew England BiolabM3003
Commercial assay or kitRNeasy micro kitQiagen74004
Commercial assay or kitKAPA Stranded mRNA-Seq KitKapa BiosystemsKK8420
Commercial assay or kitSimpleChIP chromatin IP buffersCell Signaling Technology14231
Commercial assay or kitProtein A/G agarose beadsThermo Fisher Scientific20423
Commercial assay or kitminElute reaction cleanup kitQiagen28204
Commercial assay or kitSEA blockThermo Fisher Scientific107452659
Commercial assay or kitThruplex DNA library prep kitClontechR400523
Software, algorithmPartek Flow platformPartek (https://www.partek.com/partek-flow/)
Software, algorithmDESeq2Love et al., 2014RRID:SCR_015687
Software, algorithmDAVIDHuang et al., 2009 (http://david.abcc.ncifcrf.gov/)RRID:SCR_001881
Software, algorithmHomerHeinz et al., 2010 (http://homer.ucsd.edu/)RRID:SCR_010881
Software, algorithmQuESTValouev et al., 2008
Software, algorithmGREATMcLean et al., 2010

Additional files

Supplementary file 1

Differentially expressed genes (DEG) between conditions.

https://cdn.elifesciences.org/articles/64444/elife-64444-supp1-v3.xlsx
Supplementary file 2

Differentially accessible regions (DAR) between conditions.

https://cdn.elifesciences.org/articles/64444/elife-64444-supp2-v3.xlsx
Supplementary file 3

Tcf7l1 and Lef1 peak coordinates by category as in Figure 4.

https://cdn.elifesciences.org/articles/64444/elife-64444-supp3-v3.xlsx
Supplementary file 4

Hi-C loops to transcriptional start site (TSS).

https://cdn.elifesciences.org/articles/64444/elife-64444-supp4-v3.xlsx
Supplementary file 5

Hi-C loops that are anchored on β-catenin binding sites in high CHIR99021 (CHIR) condition.

https://cdn.elifesciences.org/articles/64444/elife-64444-supp5-v3.xlsx
Supplementary file 6

Expression of Wnt target genes published in Karner et al., 2011.

https://cdn.elifesciences.org/articles/64444/elife-64444-supp6-v3.xlsx
Transparent reporting form
https://cdn.elifesciences.org/articles/64444/elife-64444-transrepform-v3.docx

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Qiuyu Guo
  2. Albert Kim
  3. Bin Li
  4. Andrew Ransick
  5. Helena Bugacov
  6. Xi Chen
  7. Nils Lindström
  8. Aaron Brown
  9. Leif Oxburgh
  10. Bing Ren
  11. Andrew P McMahon
(2021)
A β-catenin-driven switch in TCF/LEF transcription factor binding to DNA target sites promotes commitment of mammalian nephron progenitor cells
eLife 10:e64444.
https://doi.org/10.7554/eLife.64444