A β-catenin-driven switch in TCF/LEF transcription factor binding to DNA target sites promotes commitment of mammalian nephron progenitor cells
- 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
- The Rogosin Institute, United States
- Center for Molecular Medicine, Maine Medical Center Research Institute, United States
- 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
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 | Designation | Source or reference | Identifiers | Additional information |
|---|---|---|---|---|
| Antibody | Rabbit polyclonal anti-Six2 | ProteinTech | 11562-1-AP; RRID:AB_2189084 | ChIP (1:200) |
| Antibody | Mouse monoclonal IgG1 anti-Six2 | Abnova | H00010736-M01; RRID:AB_436993 | IF (1:1000) |
| Antibody | Rabbit polyclonal anti-β-catenin | ThermoFisher | 71–2700; RRID:AB_2533982 | ChIP (1:200) |
| Antibody | Rabbit monoclonal anti-non-phopho-β-catenin | Cell Signaling Technology | 8814; RRID:AB_11127203 | ChIP (1:200), IF (1:500) and western blot (1:1000) |
| Antibody | Mouse monoclonal anti-Tcf7l1 | Santa Cruz | sc-166411; RRID:AB_2302942 | IF (1:1000), western blot (1:1000) |
| Antibody | Rabbit polyclonal anti-Tcf7l1 | Thermo Scientific | PA5-40327; RRID:AB_2577173 | ChIP (1:200) |
| Antibody | Rabbit monoclonal anti-Tcf7l2 | Cell Signaling Technology | 2569; RRID:AB_2199816 | ChIP (1:200), IF (1:500), western blot (1:1000) |
| Antibody | Rabbit monoclonal anti-Tcf7 | Cell Signaling Technology | 2203; RRID:AB_2199302 | ChIP (1:200), western blot (1:1000) |
| Antibody | Rat monoclonal IgG2b anti-Tcf7 | R and D Systems | MAB8224 | IF (1:100) |
| Antibody | Rabbit monoclonal anti-Lef1 | Cell Signaling Technology | 2230; RRID:AB_823558 | ChIP (1:200), IF (1:500), western blot (1:1000) |
| Antibody | Rabbit monoclonal anti-histone H3 | Abcam | ab1791; RRID:AB_302613 | Western blot (1:1000) |
| Antibody | Rabbit monoclonal anti-histone H3K27ac | Abcam | ab4729; RRID:AB_2118291 | ChIP (1:200) |
| Antibody | Mouse monoclonal anti-histone H3K27me2me3 | Active Motif | 39536; RRID:AB_2793247 | ChIP (1:200) |
| Antibody | Mouse monoclonal anti-Ser5P-RNAPII | Millipore | 05–623; RRID:AB_309852 | ChIP (1:200) |
| Antibody | Goat polyclonal anti-Jag1 | R and D Systems | AF599; RRID:AB_2128257 | IF (1:50) |
| Antibody | HRP-conjugated goat anti-mouse IgG1 | Thermo Scientific | A10551 | Western blot (1:1000) |
| Antibody | HRP-conjugated goat anti-rabbit IgG | Cell Signaling Technology | 7074 | Western blot (1:1000) |
| Antibody | Alexa647 mouse anti-rat IgG2b | Abcam | ab172335 | IF (1:500) |
| Antibody | Alexa555 donkey anti-rabbit IgG | Abcam | ab150074 | IF (1:500) |
| Antibody | Alexa488 goat anti-mouse IgG1 | Thermo Scientific | A10551 | IF (1:500) |
| Commercial assay or kit | SuperScript IV VILO Master Mix with ezDNase Enzyme | Thermo Fisher Scientific | 11766050 | |
| Commercial assay or kit | Luna Universal qPCR Master Mix | New England Biolab | M3003 | |
| Commercial assay or kit | RNeasy micro kit | Qiagen | 74004 | |
| Commercial assay or kit | KAPA Stranded mRNA-Seq Kit | Kapa Biosystems | KK8420 | |
| Commercial assay or kit | SimpleChIP chromatin IP buffers | Cell Signaling Technology | 14231 | |
| Commercial assay or kit | Protein A/G agarose beads | Thermo Fisher Scientific | 20423 | |
| Commercial assay or kit | minElute reaction cleanup kit | Qiagen | 28204 | |
| Commercial assay or kit | SEA block | Thermo Fisher Scientific | 107452659 | |
| Commercial assay or kit | Thruplex DNA library prep kit | Clontech | R400523 | |
| Software, algorithm | Partek Flow platform | Partek (https://www.partek.com/partek-flow/) | ||
| Software, algorithm | DESeq2 | Love et al., 2014 | RRID:SCR_015687 | |
| Software, algorithm | DAVID | Huang et al., 2009 (http://david.abcc.ncifcrf.gov/) | RRID:SCR_001881 | |
| Software, algorithm | Homer | Heinz et al., 2010 (http://homer.ucsd.edu/) | RRID:SCR_010881 | |
| Software, algorithm | QuEST | Valouev et al., 2008 | ||
| Software, algorithm | GREAT | McLean et al., 2010 |
Additional files
-
Supplementary file 1
Differentially expressed genes (DEG) between conditions.
- https://cdn.elifesciences.org/articles/64444/elife-64444-supp1-v2.xlsx
-
Supplementary file 2
Differentially accessible regions (DAR) between conditions.
- https://cdn.elifesciences.org/articles/64444/elife-64444-supp2-v2.xlsx
-
Supplementary file 3
Tcf7l1 and Lef1 peak coordinates by category as in Figure 4.
- https://cdn.elifesciences.org/articles/64444/elife-64444-supp3-v2.xlsx
-
Supplementary file 4
Hi-C loops to transcriptional start site (TSS).
- https://cdn.elifesciences.org/articles/64444/elife-64444-supp4-v2.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-v2.xlsx
-
Supplementary file 6
Expression of Wnt target genes published in Karner et al., 2011).
- https://cdn.elifesciences.org/articles/64444/elife-64444-supp6-v2.xlsx
-
Transparent reporting form
- https://cdn.elifesciences.org/articles/64444/elife-64444-transrepform-v2.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)
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
