# Hnf1b renal expression directed by a distal enhancer responsive to Pax8

Laura Goea, Isabelle Buisson, Valerie Bello, Alexis Eschstruth, Melanie Paces-Fessy, Ronan Le Bouffant, Albert Chesneau, Silvia Cereghini, Jean-François Riou, Muriel Umbhauer

► To cite this version:

Laura Goea, Isabelle Buisson, Valerie Bello, Alexis Eschstruth, Melanie Paces-Fessy, et al. Hnf1b renal expression directed by a distal enhancer responsive to Pax8. Scientific Reports, 2022, 12 (1), pp.19921. 10.1038/s41598-022-21171-x. hal-03813666

HAL Id: hal-03813666

https://cnrs.hal.science/hal-03813666

Submitted on 22 Nov 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for updates

OPEN

# Hnf1b renal expression directed by a distal enhancer responsive to Pax8

L. Goea¹,³,⁴, I. Buisson¹,⁴, V. Bello¹, A. Eschstruth¹, M. Paces-Fessy¹, R. Le Bouffant¹, A. Chesneau²,³, S. Cereghini¹, J. F. Riou¹ & M. Umbhauer¹

Xenopus provides a simple and efficient model system to study nephrogenesis and explore the mechanisms causing renal developmental defects in human. Hnf1b (hepatocyte nuclear factor 1 homeobox b), a gene whose mutations are the most commonly identified genetic cause of developmental kidney disease, is required for the acquisition of a proximo-intermediate nephron segment in Xenopus as well as in mouse. Genetic networks involved in Hnf1b expression during kidney development remain poorly understood. We decided to explore the transcriptional regulation of Hnf1b in the developing Xenopus pronephros and mammalian renal cells. Using phylogenetic footprinting, we identified an evolutionary conserved sequence (CNS1) located several kilobases (kb) upstream the Hnf1b transcription start and harboring epigenomic marks characteristics of a distal enhancer in embryonic and adult renal cells in mammals. By means of functional expression assays in Xenopus and mammalian renal cell lines we showed that CNS1 displays enhancer activity in renal tissue. Using CRISPR/cas9 editing in Xenopus tropicalis, we demonstrated the in vivo functional relevance of CNS1 in driving hnf1b expression in the pronephros. We further showed the importance of Pax8-CNS1 interaction for CNS1 enhancer activity allowing us to conclude that Hnf1b is a direct target of Pax8. Our work identified for the first time a Hnf1b renal specific enhancer and may open important perspectives into the diagnosis for congenital kidney anomalies in human, as well as modeling HNF1B-related diseases.

The hepatocyte nuclear factor 1 homeobox B (Hnf1b) is a transcription factor (TF) involved in the development and homeostasis of many organs in vertebrates, including liver, pancreas, hindbrain, genital tract and kidney. In the mammalian embryonic kidney, Hnf1b is expressed in tissues leading both to the collecting system and the nephrons. Gene targeting studies have shown that Hnf1b plays a crucial role at several stages of kidney development, including ureteric bud branching and early nephrogenesis as well as nephron patterning and epithelial tubular organization of the collecting ducts¹⁻⁵. The pivotal role of Hnf1b in renal tubule fate is underscored by the finding that combined expression of the transcription factors Emx2, Hnf1b, Hnf4a and Pax8 is sufficient to reprogram fibroblasts into renal tubular epithelial cells⁶. Altogether, these studies led to the identification of many Hnf1b target genes and uncovered Hnf1b coordinated transcriptional circuits involved in nephrogenesis and epithelial differentiation¹⁻⁵,⁷⁻⁹. By contrast, regulation of Hnf1b expression in the different embryonic renal compartments remains essentially unknown. Retinoic acid signaling has been shown to directly regulate Hnf1b expression in the mouse developing neural tube through a conserved enhancer element but this enhancer is not active in the developing mouse kidney (Ref.¹⁰ and S. Cereghini, unpublished).

Heterozygous mutations in HNF1B are one of the most commonly identified genetic cause of developmental kidney disease in human. They are associated with several disorders such as the Renal Cysts and Diabetes Syndrome (RCAD, OMIM #137920), a multi-organ disease, characterized by kidney abnormalities and early diabetes¹¹. Molecular genetic defects comprise whole-gene deletions and intragenic mutations, including missense, frameshift, nonsense, and small insertions or deletions¹². Mutations are inherited in an autosomal dominant pattern, although up to 50% of mutations occur de novo¹³,¹⁴. The phenotype of HNF1B mutant carriers is highly variable within families¹⁴. Consistent with these observations, heterozygous mutants in a novel mouse

¹Laboratoire de Biologie du Développement, CNRS, Institut de Biologie Paris Seine, IBPS, UMR7622, Sorbonne Université, 75005 Paris, France. ²Paris-Saclay Institute of Neuroscience, CNRS, CERTO-Retina France, Université Paris-Saclay, 91405 Orsay, France. ³Present address: R&D, Research & Development, Pharmaceuticals - Precision Nephrology, Bayer AG, 8, Building 0520, 8.025, 42096 Wuppertal, Germany. ⁴These authors contributed equally: L. Goea and I. Buisson. email: muriel.umbhauer@sorbonne-universite.fr
**Figure 1.** Identification of a conserved non-coding sequence as candidate for *Hnf1b* enhancer regulated by Pax8. (a) Localization of two conserved non-coding sequence (CNS1 and CNS2, red bars) in *M. musculus* and *X. tropicalis*. CNS1 is localized in the intergenic region between *Hnf1b* and the upstream gene *Heatr6*. CNS2 is localized in the 4th intron of *Hnf1b*. Arrows indicate the transcription initiation of *Hnf1b*. Protein coding exons are represented by blue bars, exonic sequences corresponding to 5' and 3' UTRs are in yellow. (b) Nucleotide sequence comparison of CNS1 in different species. The table indicates the % of conserved nucleotides. (c) Nucleotide sequence alignment of CNS1 in *X. tropicalis*, *H. sapiens* and *M. musculus*. Conserved nucleotides are in red. Underlined are conserved putative binding sites of TFs known to be expressed in the developing kidney according to (https://sckidney flatironinstitute.org)⁷⁹, GUDMAP (https://www.gudmap.org)⁸⁰,⁸¹ or Xenbase (https://www.xenbase.org/entry/). The blue box highlights the conserved Pax8 putative binding site. Sequences targeted by the gRNAa and gRNAb used for CRISPR editing in Fig. 6 are shown in green. (d) Pax8 binding site sequence identified in CNS1 and Pax8 binding consensus sequence as defined by Ref.³¹ shown as sequence logo. For each position, the size of the characters represents the relative frequency of the corresponding position. The lower sequence is the Pax8 binding sequence in CNS1. Nucleotides mutated in the Pax8-BS mut oligonucleotide used for Electrophoretic Mobility Shift Assay (EMSA), and in plasmids CNS1 mut-Luc and CNS1 mut-eGFP are labelled in red. (e) EMSA. Biotin-labelled oligonucleotide probes were incubated with nuclear extracts from MDCK as indicated: Pax8-BS, Pax8-BS mut and Control TPO (Pax8 binding site identified in the thyroperoxidase promoter, used as a positive control). An arrow marks the Pax8 specific complexes. Specificity of the Pax8 complex was demonstrated by competition with 10- and 100-fold molar excesses of unlabelled Pax8-BS oligonucleotide (compare lanes 5 and 7) as well as incubation with 1 µg of Pax8 monoclonal antibodies (Pax8 AB) leading to a decrease in the Pax8 complex signal (compare lanes 5 with 9). The figure is representative of 3 independent experiments. Original uncropped autoradiography of the whole membrane is presented in Supplementary Fig. S5.

<figure><img src="image_3.png" /></figure>

model carrying a splice donor site mutation identified in humans, exhibited a range of renal phenotypes similar to those described in human and, with some variability, decreased HNF1B protein levels from the normal allele⁸. Thus, any mechanism that modulates HNF1B expression levels may therefore impact on the variable pathogenic phenotypes in *HNF1B* heterozygous mutant carriers. Identification of the mechanisms underlying *HNF1B* expression is therefore of great importance for the understanding of renal anomalies associated with *HNF1B* mutations. It will potentially uncover novel mutations affecting HNF1B protein levels because they are associated with *HNF1B* regulatory sequences or TFs interacting with these sequences.

*Xenopus* provides a simple and efficient model system to study nephrogenesis and explore the mechanisms that cause renal developmental defects in human¹⁵,¹⁶. The *Xenopus* tadpole contains only one functional large nephron, the pronephros, displaying a structural and functional organization very similar to the mammalian nephron of the metanephros¹⁷,¹⁸. Many developmental genes that govern pronephros formation have similar functions during mammalian kidney development. Pax8 has been shown to be a major actor of renal specification and epithelialization in *Xenopus*; its depletion leads to a complete absence of pronephric tubule¹⁹⁻²¹. *Pax2* is expressed in the developing pronephros at later stages than *pax8*, once renal cells have been specified, and is required for the expression of several terminal differentiation markers of the pronephric tubule¹⁹. The expression of *hnf1b* begins at the late neurula stage in the kidney field and is maintained throughout the entire pronephros, with the highest levels in the proximal region of the pronephric tubule anlage at tailbud stages. Overexpression of an *Hnf1b*-dominant negative construct in *Xenopus* embryos leads to a strong reduction of the proximal and intermediate pronephric segments and downregulation of notch ligand *delta-like 1* (*dll1*) and iroquois (*irx1/2/3*) gene expression³. These observations are in agreement with the consequences of specific inactivation of *Hnf1b* in mouse nephron progenitors as well as in *Hnf1b*-deficient zebrafish embryos, pointing *hnf1b* as a major regulator of nephron tubular segmentation in vertebrates³,⁵,²².

Some lines of evidence in *Xenopus* and zebrafish indicate that *pax8* and/or *pax2* are regulators of *hnf1b* expression. *Pax2a/8*-deficient zebrafish embryos fail to activate expression of *hnf1ba* in the intermediate mesoderm. In *Xenopus*, *pax8* loss of function leads to inhibition of *hnf1b* expression in the kidney field as soon as the late neurula stage while expression of *lhx1* and *osr2* is not affected. In addition, *hnf1b* is significantly up-regulated in blastula animal caps in response to *pax8* ectopic expression even in absence of protein synthesis, suggesting that *hnf1b* might be a direct target of *pax8*¹⁹. Together, these results suggest a model where *pax2* and/or *pax8* act upstream of *hnf1b*, which in turn results in the initiation of segment-specific gene expression programs²².

In the present study, we have explored the transcriptional regulation of *Hnf1b* in the developing *Xenopus* pronephros and mammalian renal cells. Using phylogenetic footprinting, we searched for evolutionary conserved non-coding sequences as potential candidates of *Hnf1b* enhancers. We identified a 306-pb sequence (CNS1) located several kb upstream the *Hnf1b* transcription start site. CNS1 sequence is conserved in *Xenopus*, mouse and human and bears epigenomic marks characteristics of a distal enhancer in embryonic and adult renal cells in mammals. Enhancer activity of CNS1 was established by means of reporter expression assays in *Xenopus* embryos and mammalian renal cell lines. We demonstrated that CNS1 is able to drive reporter gene expression in MDCK (Madin-Darby Canine Kidney) and IMCD3 (inner medullary collecting duct) cell lines as well as in the pronephros in *Xenopus* embryos. In vivo functional relevance of CNS1 in endogenous *hnf1b* pronephric expression was proven using CRISPR/cas9 in *Xenopus*. Inspection of conserved TF binding sites in CNS1 revealed the presence of a Pax8 binding site. We further demonstrated the importance of Pax8-interaction for CNS1 enhancer activity allowing us to conclude that *Hnf1b* is a direct target of Pax8.
a **M.musculus**

<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /></figure>

X.tropicalis

<figure><img src="image_6.png" /></figure>

b

<table><thead><tr><th></th><th>H. sapiens</th><th>M. musculus</th><th>M. domestica</th><th>G.gallus</th><th>X. tropicalis</th><th>X. laevis.S</th><th>X. laevis.L</th></tr></thead><tbody><tr><td>H. sapiens</td><td>100.00</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>M. musculus</td><td>96.73</td><td>100.00</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>M. domestica</td><td>86.93</td><td>85.62</td><td>100.00</td><td></td><td></td><td></td><td></td></tr><tr><td>G. gallus</td><td>63.95</td><td>63.61</td><td>64.86</td><td>100.00</td><td></td><td></td><td></td></tr><tr><td>X. tropicalis</td><td>76.57</td><td>75.58</td><td>77.70</td><td>61.02</td><td>100.00</td><td></td><td></td></tr><tr><td>X. laevis.S</td><td>74.17</td><td>73.18</td><td>75.99</td><td>60.20</td><td>95.72</td><td>100.00</td><td></td></tr><tr><td>X. laevis.L</td><td>74.59</td><td>73.93</td><td>77.70</td><td>61.36</td><td>95.41</td><td>92.76</td><td>100.00</td></tr></tbody></table>

c

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>

Tcf4

<figure><img src="image_12.png" /></figure>

RORa

Myb

Left1

Tcf4

Tcf4

Cebp

Hnf6

Hmga

gRNAb-CCCTTTTCTCCAGTGGCCTGATA

X. tropicalis TCCCCCTCCCCTTTTCTC 305

H. sapiens CCCTCCCTGTCTCTTTCTC 306

M. musculus CCCTCCCTGTCTCTTTCTC 306

d

<figure><img src="image_13.png" /></figure>

e

<figure><img src="image_14.png" /></figure>
## Results

### Identification of a conserved non-coding sequence as candidate for *Hnf1b* enhancer regulated by *Pax8*.
In order to identify regulatory sequences that activate *Hnf1b* gene expression in the developing kidney, we searched for evolutionary conserved non-coding sequences (CNSs) as candidates for enhancers. Using in silico methods that allow detailed inspection of sequence conservation profiles across several genomes<sup>23,24</sup>, we compared the genomic sequence of a 157 kb segment encompassing the mouse (*Mus musculus*) *Hnf1b* gene with the orthologous intervals in human (*Homo sapiens*), rat (*Rattus norvegicus*), opossum (*Monodelphis domestica*), chicken (*Gallus gallus*), frog (*Xenopus tropicalis*) and zebrafish (*Danio rerio*). This analysis identified one CNS (called CNS1) conserved in all species, except zebrafish, located 33.4 kb upstream the *Hnf1b* transcription start site in mouse and a second one (CNS2) located in the 4th intron of *Hnf1b* conserved in all species examined (Fig. 1a,b). A survey of the literature revealed that CNS2 is identical to a previously characterized *Hnf1b* enhancer whose activity is restricted to *Hnf1b* neural expression in transgenic mice<sup>10</sup>. We therefore focused our analysis on the 306-bp genomic region CNS1.

Active regulatory sequences show epigenomic properties such as characteristic patterns of histone modification and accessible chromatin. CNS1 was queried against public databases for Deoxyribonuclease I (DNase I) hypersensitive site sequencing (DNase-seq), Assay for Transposase-Accessible Chromatin sequencing (ATAC-seq) and chromatin-immunoprecipitation sequencing (ChIP-seq) for histones marks or the transcriptional cofactor p300 from human, mouse (ENCODE project<sup>25,26</sup>) and *Xenopus*<sup>27</sup> embryonic and adult tissues or cell lines. We found that CNS1 correlates with a DNase I hypersensitive site in several human primary cells from adult kidney tubule known to express HNF1B at a high level but not in Caco-2 or A549 cell lines derived from colon carcinomas and lung adenocarcinomas respectively, which also express HNF1B although at a lower level (https://www.proteinatlas.org) (Fig. 2a, Fig. S1a). No DNase I hypersensitive signal was detected at the CNS1 location in spinal cord or neural progenitors (Fig. 2a, Fig. S1a). In the mouse embryo, *Hnf1b* is expressed during kidney development but also in the developing liver, pancreas and lung<sup>28,29</sup>. We looked for CNS1 epigenetic status in these embryonic tissues. We observed that CNS1 showed characteristics of an active enhancer in E15.5 embryonic kidney: it overlaps with an ATAC-seq signal which is bordered by H3K4me1, H3K4me2 and H3K27ac signals. Interestingly, no such signal was observed in lung or liver tissue (Fig. 2b,c, Fig. S1b). In *Xenopus*, we also found CNS1 associated with active enhancer marks (p300, H3K4me1) in ChIP-seq data from whole neurula and tailbud stages embryos (stages 16 and 30, respectively) when *hnf1b* is expressed in the developing pronephros (Fig. 2b,d, Fig. S1c).

Altogether, these data suggest that CNS1 is an enhancer whose activity is potentially specific to renal tissue. Notably, CNS2 is not localized in a DNase I hypersensitive site in any of the analyzed cell lines and does not correlate with H3K27ac enrichment in mouse embryonic kidney (Fig. S1a,b).

We then searched for conserved TF putative binding motifs in CNS1 using R-Vista 2.0<sup>30</sup>. We identified 23 putative TF binding sites conserved between frog, mouse and human, including Pax8 whose knockdown in *Xenopus* leads to inhibition of *hnf1b* expression in the kidney field<sup>19</sup>. Eight additional ones were conserved solely between frog and human (Tables S1, S2). The sequence of the Pax8 putative binding site in CNS1 shows high similarity to known Pax8 binding sites<sup>31</sup> and is strictly conserved in mouse and human CNS1 (Fig. 1c,d). The actual binding of Pax8 to the identified sequence was tested by electrophoretic mobility shift assays (EMSAs). As a positive control we used a known Pax8-binding site derived from the thyroperoxidase promoter, a target gene of PAX8<sup>32</sup> (Control TPO, Fig. 1e). When incubated with nuclear extracts from MDCK cells containing abundant levels of Pax8 protein<sup>33</sup>, biotin-labelled Pax8-BS and control-TPO oligonucleotides formed a similar protein-DNA complex, as judged by their identical electrophoretic mobility (Fig. 1e). When excess unlabelled wild-type oligonucleotides were added as competitors, they weakened or abolished completely the shifted bands (Fig. 1e, compare lanes 2 and 3 for control TPO and lanes 5, 6 and 7 for Pax8-BS). No complex was detected with biotinylated oligonucleotides containing mutations known to abolish Pax8 binding<sup>34</sup>. The presence of Pax8 antibodies in the reaction did not induce a supershift, but decreased the signal intensity of the electrophoretic complexes, suggesting that antibody interaction disrupts Pax8 interaction with both control-TPO and Pax8-BS oligonucleotides (Fig. 1e, compare lanes 2 and 4 for control TPO and lanes 5 and 9 for Pax8-BS), further confirming the specificity of Pax8 binding. Together, these data show that Pax8 is able to specifically bind in vitro to the in silico identified Pax8 binding site of CNS1.

### CNS1 is specifically activated by ectopic expression of *pax8* and *pax2* and functions as an enhancer in mammalian renal cells.
We subsequently tested the ability of different versions of exogenous *pax8* to transactivate CNS1 in HEK293 cells that do not express Pax8 protein (www.proteinatlas.org/ENSG00000125618-PAX8/cell) using luciferase assays. The CNS1 sequence was cloned into the luciferase pGL4.23[*luc2/minP*] vector and the resulting construct CNS1-Luc was used to transfect HEK293 cells together with a control vector carrying *Renilla* luciferase gene used to normalize for transfection efficiency. We analysed the effect of expressing wild-type *pax8*, *pax8ΔO* (a mutated version of *pax8* with deletion of the octapeptide motif mediating interaction with *Groucho/Grg4* co-repressors responsible for transcriptional repressor activity<sup>35,36</sup>), and *pax8VP16* (a chimeric protein corresponding to *pax8* fused to the potent transcriptional activation domain of the herpes simplex virus VP16 and conferring strong transcriptional activity to *pax8*<sup>19</sup>). HEK293 cells were transfected with CNS1-Luc vector alone or associated with one of the *pax8* expression vectors (*pax8*, *pax8VP16*, *pax8ΔO*). Expression of either *pax8*, *pax8VP16* or *pax8ΔO* led all to an increase in relative luciferase activity. As expected, this increase was much higher for *pax8ΔO* and *pax8VP16* (~40-fold) than for wild-type *pax8* (~five-fold), and further suggested that *Groucho/Grg4* repression is also at work in HEK293 cells (Fig. 3a,b). Mutation of the *Pax8* binding site in CNS1 almost abolished this activation in every instance (Fig. 3a,b). Our results therefore show that CNS1 transactivation in HEK293 cells is dependent on the expression of exogenous *pax8*, and the
<figure><img src="image_15.png" /><figcaption>**Figure 2.** Chromatin signature of CNS1 region in *Homo sapiens*, *Mus musculus* and *Xenopus tropicalis*. (a) Snapshots of representative DNase-seq tracks from ENCODE database in human (https://www.encodeproject.org)²⁵,²⁶. The overall profiles spanning a larger region including CNS2 and the description of the biological materials are shown in Fig. S1a. Open chromatin state is clearly detected at CNS1 level in embryonic kidney tissue and renal cells. In contrast, CNS1 is not accessible to DNase in non-renal cells expressing HNF1B (Caco-2 and A549 cells) as well as in spinal cord or neural progenitors where HNF1B is not expressed. (b) Characteristics of enrichment patterns for histone modifications and p300. (c) Snapshots of representative ATAC-seq and histone modifications tracks as indicated from ENCODE database (https://www.encodeproject.org)²⁵,²⁶ in kidney, lung and liver of E15.5 mouse embryo. The overall profiles spanning a larger region including CNS2 is shown in Fig. S2b. Although Hnf1b is expressed in developing kidney, liver and lung, CNS1 shows characteristics of an active enhancer only in the developing kidney. (d) Snapshots of representative p300 and histone modifications tracks as indicated from the epigenome reference maps in Ref.²⁷ in *Xenopus tropicalis* embryos at stages 10.5, 16 and 30. The overall profiles spanning a larger region including CNS2 is shown in Fig. S1c. CNS1 is associated with active enhancer marks (p300, H3K4me1) in embryos at neurula and tailbud stages (stage 16 and 30, respectively) when *hnf1b* is expressed in the developing pronephros but not at early gastrula (stage 10.5). In (a,c,d), the blue box indicates the 306-bp CNS1 location.</figcaption></figure>
<figure><img src="image_16.png" /><figcaption>b</figcaption></figure>

<figure><img src="image_17.png" /><figcaption>c</figcaption></figure>

<figure><img src="image_18.png" /><figcaption>d</figcaption></figure>

<figure><img src="image_19.png" /></figure>
►**Figure 3.** CNS1 is specifically activated by ectopic expression of pax8 and pax2 and functions as an enhancer in mammalian renal cells. (a) Reporter constructs used in luciferase assay experiments. The backbone plasmid pGL4.23[*luc2miniP*] was used either unmodified (Empty-Luc) or modified by insertion of the wild type CNS1 (CNS1-Luc) or CNS1 mutated in the Pax8 binding site (CNS1 mut-Luc) upstream of the minimal promoter. (b) Exogenous pax8 and pax2 are able to activate CNS1 enhancer activation in HEK293 cells. Expression vectors for pax8, pax2, pax8ΔO or pax8VP16 were co-transfected with CNS1-Luc or CNS1mut-Luc, and a control Renilla luciferase vector for normalization. Normalized luciferase activity was determined as the Firefly over Renilla luciferase activity. Fold change is expressed as the ratio of normalized values on the normalized value for CNS1-Luc vector alone. Expression of pax8, pax8ΔO, pax8VP16 and pax2, all result in an increase of reporter activation. In every case, this increase is abolished when the mutant construct, CNS1mut-Luc, is used instead of CNS1-Luc. Combined results from three (pax8, pax8ΔO and pax8VP16) or four (pax2) independent experiments. Statistical significance was determined using One-Way ANOVA followed by Tuckey's multiple comparison test. P-values are annotated on the figure (c) Pax8 immunofluorescence staining in MDCK and IMCD3 cells and counterstaining of nuclei with DAPI. Pax8 is detected in MDCK and IMCD3 cell nuclei. Scale bar: 50 µm. (d) Analysis of CNS1 activity in luciferase reporter assay in MDCK and IMCD3 cell lines. Histograms show the relative luciferase activities in MDCK and IMCD3 cells transfected with Empty-Luc, CNS1-Luc or CNS1 mut-Luc together with a control Renilla luciferase vector for normalization. Fold change is expressed as the ratio of normalized values on the normalized value for Empty-Luc. Combined results from four independent experiments. Statistical significance was determined using One-Way ANOVA followed by Tuckey's multiple comparison test. P-values are annotated on the figure.

integrity of a sequence that binds pax8 in vitro, strongly suggesting that this process results from the interaction of exogenous pax8 with this sequence.

The vertebrates *Pax* genes are classified into four sub-groups and *Pax8* belongs to sub-group 2 together with *Pax2* and *Pax5* exhibiting similar recognition sequences. Since *Pax2*, but not *Pax5*, is expressed in the developing kidney³⁷ (GUDMAP https://www.gudmap.org), we evaluated the ability of pax2 to increase transcriptional activity of the CNS1-Luc reporter construct in HEK293 cells. *Pax2* expression induced luciferase activity from the CNS1-Luc construct and was unable to activate transcription from CNS1 mut-Luc (Fig. 3a,b). These results indicate that both *pax2* and *pax8* are able to activate CNS1 in HEK293 cells when ectopically expressed.

We then investigated the functionality of the identified CNS1 in the kidney derived cell lines MDCK and IMCD3, which both express *Pax8* (Fig. 3c) as well as *Hnf1b*³³,³⁸. Luciferase activity of CNS1-Luc transfected cells was significantly increased compared to cells transfected with the empty luciferase vector pGL4.23[*luc2/minP*] (Empty-Luc) in both cell lines (Fig. 3d) although to a greater extent in IMCD3 in comparison to MDCK cells. Mutating *Pax8* binding site led to a significant decrease of reporter activity in transfected MDCK cells but not in IMCD3 (Fig. 3d).

Thus, CNS1 enhancer is active in mammalian renal cell lines and its dependency upon the *Pax8* binding site varies between cell lines.

**CNS1 functions as a nephric enhancer in vivo in Xenopus embryos and its activity is partly dependent on the transcription factor pax8.** CNS1 enhancer activity was then investigated in vivo in *Xenopus laevis* by means of F0 transgenesis³⁹. A plasmid containing CNS1 upstream of a β-globin basal promoter and eGFP coding sequence was used to generate transgenic embryos through the *I-SceI* meganuclease method⁴⁰. Reporter expression was examined at different developmental stages by in situ hybridization for eGFP mRNA. The control reporter construct without CNS1 (β-globin-eGFP), did not lead to any significant GFP expression at all stages examined (Fig. 4a,c,d). In contrast, CNS1 drove strong and reproducible expression in the developing pronephros starting at the mid tailbud stage 25 (Fig. 4a–d). Transverse sections showed expression in the developing tubule anlage but not in the splanchnic mesoderm from which the glomus forms (Fig. 4c). GFP reporter expression was therefore similar to the endogenous *hnf1b* expression in the developing pronephros except that endogenous *hnf1b* mRNA is detected earlier, as soon as the late neurula stage 19 in the kidney field (Fig. S2). GFP expression was additionally observed in posterior somites in CNS1-eGFP embryos, although somites do not express endogenous *hnf1b* (Fig. S3). Using the REMI transgenesis method that allows the production of fully transgenic-non-mosaic embryos with high efficiency⁴¹, we confirmed that CNS1 was able to drive reporter expression in the pronephros, as shown by GFP fluorescence in the pronephros in 32.6% of the injected embryos (Fig. 4e).

Mutation of the *Pax8* binding site greatly diminished the frequency of pronephric GFP expression at stage 28 when compared to the wild-type CNS1-eGFP construct (Fig. 4c,d). However, GFP expression at stage 35 was not affected, although *pax8* is known to be expressed in the pronephric tubule at this stage³⁷ (Fig. 4c,d). These results suggest that additional inputs are present in *Xenopus* embryos at stage 35 that can overcome the mutation of the *Pax8* binding site for CNS1 activity.

*Pax2* and *pax8* are both expressed in the developing pronephros at tailbud stage 28, but, in contrast to *pax8*, *pax2* is not required for *hnf1b* expression at tailbud stage¹⁹ (Fig. 5a). Although unlikely, it is still possible that in stage 28 transgenic embryos *pax2* may be involved in CNS1-driven reporter activation instead of *pax8*. We therefore monitored the effect of *pax2* depletion upon reporter gene activity in CNS1-eGFP transgenic embryos. After injecting CNS1-eGFP construct at the 1-cell stage for transgenesis, a *pax2* antisense morpholinos (Mo *pax2*) was injected at the 2-cell stage into the two blastomeres, and reporter expression analysed at stage 28 (Fig. 5b). For each experiment, the efficiency of Mo *pax2* was confirmed by verifying inhibition of *clcnkb1* expression¹⁹ in a separate group of embryos cultured until stage 35 (not shown). Comparison of the percentage of embryos
<figure><img src="image_20.png" /><figcaption>b</figcaption></figure>

<figure><img src="image_21.png" /></figure>

Percentage of transgenic embryos
expressing GFP in kidney (%)

<figure><img src="image_22.png" /><figcaption>d</figcaption></figure>

<figure><img src="image_23.png" /><figcaption>e</figcaption></figure>

<figure><img src="image_24.png" /></figure>

**Figure 4.** Analysis of CNS1 activity in transgenesis assay in *Xenopus laevis*. (a) Reporter constructs used in transgenesis. Wild type CNS1 or CNS1 mutated in the Pax8 binding site (CNS1 mut) was cloned upstream of the β-globin basal promoter and the eGFP reporter gene. (b) Percentage of F0 CNS1-eGFP transgenic embryos generated by the I-SceI meganuclease method showing GFP expression detected by in situ hybridization in the developing pronephros at stages 20, 22, 25 and 28. Results corresponding to three independent experiments. (c) GFP expression detected by in situ hybridization in representative F0 transgenic embryos at stages 28 and 35 generated with reporter constructs as indicated. Transverse cryostat section (16 µm) of CNS1-eGFP transgenic embryos showing GFP reporter gene expression at the level of the proximal part of the pronephros are shown on the right. At stage 28, the dotted line delineates the mesoderm-endoderm frontier. sp splanchnic mesoderm. (d) Percentage of F0 transgenic embryos generated with the indicated reporter constructs expressing GFP in pronephros at stages 28 and 35. Results from three independent experiments. Statistical significance was determined using Fisher's exact test. P-values are annotated on the figure. (e) Representative stage 40 F0 transgenic tadpole obtained by the REMI method with the CNS1-eGFP construct. Arrowhead indicates the fluorescent pronephros. In (b,d,e), n indicates the total number of analysed embryos.
<figure><img src="image_25.png" /></figure>

hnf1b in situ hybridization

<table><thead><tr><th></th><th>Uninjected side</th><th>Injected side</th></tr></thead><tbody><tr><th>cMo</th><td>0.4 mm</td><td></td></tr><tr><th>Mo pax2</th><td></td><td></td></tr><tr><th>Mo pax8</th><td></td><td></td></tr></tbody></table>

<figure><img src="image_26.png" /></figure>

**Figure 5.** Endogenous pax8, but not pax2, is required for *hnf1b* expression and CNS1 activity in the pronephros at tailbud stage. (a) Pax2 depletion does not affect *hnf1b* expression at stage 28. 4-cell stage embryos were injected into the two left blastomeres with either a control Mo (cMo), a mix of two morpholinos, Mo Pax8 and Mo Pax8.B, (named Mo pax8 for simplicity) or Mo Pax2.2 (named Mo pax2). They were further cultured until tailbud stage 28 and processed for *hnf1b* whole-mount in situ hybridization. *Hnf1b* pronephric expression is strongly inhibited in 80% (n=24) of pax8 depletion. In contrast, pax2 depletion has no significant effect on *hnf1b* pronephric expression (n=60), as embryos injected with cMo (n=71) (three independent experiments). (b) Pax2 depletion does not affect CNS1-driven GFP expression in the pronephros. F0 CNS1-eGFP transgenic embryos were generated by the I-SceI meganuclease method at the 1-cell stage. At the 2-cell stage, same mixes of Mo as in (a) were injected into every blastomeres. GFP expression was monitored by whole-mount in situ hybridization in embryos cultured until the tailbud stage 28. The histogram indicates the percentage of CNS1-eGFP transgenic embryos showing GFP expression in controls, or upon depletion of pax2 or pax8. In contrast to pax8 depletion, pax2 depletion has no significant effect on CNS1 activity. Statistical significance was determined using Fisher's exact test. P-values are annotated on the figure. Combined results from three independent experiments. n indicates the total number of analysed embryos.

showing GFP mRNAs expression in the developing pronephros of embryos injected with Mo pax2, or with a control morpholino (cMo), did not show any difference (Fig. 5b). These observations show that endogenous pax2 is not required for CNS1 activity in the pronephros at tailbud stage 28. In contrast, when morpholinos targeting pax8 (Mo pax8) were used, the vast majority of injected embryos did not show any GFP signal (Fig. 5b). It is unclear however if this results from the absence of pax8 protein binding to CNS1, or because early development of the pronephric tubule is severely affected by pax8 depletion as previously shown¹⁹.

Altogether, these observations argue for a role of CNS1 in pax8 regulation of *hnf1b* at tailbud stage but also indicate that this enhancer is likely to be regulated by additional not yet identified inputs at later developmental stages.

**CNS1 is required for endogenous *hnf1b* expression in the developing pronephros.** In order to evaluate the contribution of CNS1 in its normal genomic environment to the regulation of endogenous *hnf1b* expression, we attempted to mutate CNS1 using the available CRISPR/Cas9 genome editing approach in *Xenopus tropicalis*. We designed two gRNA (gRNAa and b) targeting two sequences encompassing a 180 bp region of CNS1 that includes the pax8 binding sequence (Fig. 1c). TIDE analysis indicated that gRNAa and gRNAb efficiently target CNS1 in F0 embryos (Fig. 6a). We co-injected gRNAa and gRNAb with Cas9 protein in *Xenopus tropicalis* embryos with the objective to delete the CNS1 DNA sequence comprised between sgRNAa and sgRNAb targeted sites. Electrophoresis analysis of PCR products amplified with primers surrounding CNS1 revealed a lower band corresponding to a truncated sequence in 3 out of 5 embryos, indicating that the deletion of the targeted sequence occurred with high frequency (Fig. 6b). Embryos were injected with gRNAa and
a

Chromatogram - single embryo

uninjected

<figure><img src="image_27.png" /></figure>

TIDE analysis - Average of 4 embryos

gRNAa CRISPR : Total efficiency = 51.3%

<figure><img src="image_28.png" /></figure>

Chromatogram - single embryo

uninjected

<figure><img src="image_29.png" /></figure>

TIDE analysis - Average of 4 embryos

gRNAb CRISPR : Total efficiency = 71.7%

<figure><img src="image_30.png" /></figure>

b

<figure><img src="image_31.png" /></figure>

c

<figure><img src="image_32.png" /></figure>

<figure><img src="image_33.png" /></figure>

d

<figure><img src="image_34.png" /></figure>
►**Figure 6.** CNS1 CRISPR editing inhibits endogenous *hnf1b* pronephric expression in *Xenopus tropicalis* embryos. (a) GRNAa and gRNAAb efficiently edit *Xenopus tropicalis* embryo DNA in two different locations in CNS1. Embryos were injected at the 1-cell stage with Cas9 protein and gRNAa or gRNAAb and cultured until neurula stage. The target sequences of gRNAa and gRNAAb are depicted in Fig. 1 (see also “Methods” section). Top panels: chromatogram showing CRISPR editing by gRNAa and gRNAAb in single embryos. Lower panels: TIDE analysis of sequence trace degradation at the expected Cas9 site of DNA cleavage. Percentage of CNS1 sequence containing insertions and deletions (Indels) represented as the mean from four embryos. The use of gRNAa and gRNAAb resulted in 51.3% and 71.7% editing efficiency, respectively. (b) The use of gRNAa and gRNAAb together efficiently leads to nucleotide sequence deletion into CNS1. Embryos were injected at the 1-cell stage with Cas9 protein, together with gRNAa, gRNAAb or both. At tailbud stage, their DNA were analysed by PCR using primers allowing amplification of a DNA fragment containing the entire CNS1. Electrophoresis of amplified fragments from single embryos. A band corresponding to the wild type sequence is observed in all embryos. An additional lower band is present in 3 out of 5 embryos co-injected with gRNAa and gRNAAb, indicating that a deletion occurred (c) *Hnf1b* pronephric expression in gRNAa, gRNAAb and gRNAa-gRNAAb CRISPR embryos analyzed by in situ hybridization at tailbud stage 28. Embryos were classified into three groups according to their *hnf1b* pronephric expression (normal, slightly diminished, strongly diminished). The histogram shows the percentage of embryos for each group. Average values from four independent experiments. Statistical significance was determined using Fisher's exact test (comparison to the uninjected control embryos). P-values are annotated on the figure. n indicates the total number of analysed embryos. (d) *Lhx1* pronephric expression in gRNAa, gRNAAb and gRNAa-gRNAAb CRISPR embryos analyzed by in situ hybridization at tailbud stage 28. No significant effect on *lhx1* expression is observed. n indicates the total number of embryos analysed from three independent experiments.

gRNAAb alone or in combination, reared to tailbud stage 28 and analyzed for *hnf1b* pronephric expression by in situ hybridization. About 50% of the embryos co-injected with gRNAa and gRNAAb displayed reduced *hnf1b* expression in the pronephros in comparison to control uninjected embryos, with a strong reduction for half of them. Injection of gRNAa or gRNAAb alone led to *hnf1b* expression defect in a much smaller percentage (Fig. 6c). When *lhx1* pronephric expression was analyzed, we did not see any effect of the gRNAs alone or in combination, highlighting the specific regulatory effect of CNS1 on *hnf1b* expression (Fig. 6d). We also analyzed the effect of editing the *solute carrier family 45 member 2 (slc45a2)* gene which is not expected to result in a kidney phenotype$^{42,43}$. Although we observed a slight decrease of *hnf1b* pronephric expression in some embryos, we did not see any strong effect comparable to those observed for CNS1 editing (Fig. S4). Together, these results indicate that CNS1 is required in vivo for the proper *hnf1b* pronephric expression at tailbud stage.

Discussion

Enhancers confer tissue-specific gene expression required for proper development and homeostasis. In this study,
we have identified an evolutionary conserved sequence (CNS1) located several kb upstream of the *Hnf1b* tran-
scription start site (TSS) showing enhancer activity in renal tissue. We further demonstrated that CNS1 binds and
is responsive to the TF Pax8, and is critical in vivo for *hnf1b* expression in the *Xenopus* developing pronephros.

Numerous methods have been developed with success to identify putative enhancers across entire genomes,
based on evolutionary DNA sequence conservation, chromatin accessibility, histone modifications and chromatin
interactions⁴⁴⁻⁴⁶. Whatever the method, finding out which putative enhancers are functional, which patterns of
expression they control, and how they encode these expression patterns is always a challenge. Our work high-
lights and confirms the power of the *Xenopus* model for comparative genomics-based approach coupled with
transgenesis assay²³. Although *hnf1b* is expressed in other embryonic tissues such as the neural tissue and the
endoderm, CNS1 drives expression in *hnf1b* pronephric territory only, indicating it is a tissue-specific enhancer.
During neurulation, we did not observe any GFP reporter expression in the kidney field suggesting that CNS1
is dispensable for the initiation of *hnf1b* expression in the developing kidney but is rather critical for its main-
tenance. In order to have a more comprehensive understanding of the regulatory sequences governing *hnf1b*
spatiotemporal expression, other approaches not depending to the strict sequence conservation should be used. It
is indeed acknowledged that although genes involved in vertebrate development are unusually enriched for highly
conserved non-coding sequence elements, enhancer function can be conserved without discernible sequence
conservation⁴⁷. The number, location, and type of TF binding sites can change within an enhancer over evolution
without loss of functionality⁴⁸. Notably, we could not find in the zebrafish genome any sequence exhibiting high
similarity with CNS1. It may be because the *Xenopus* genome is much closer to that of mammalian genomes
than the zebrafish genome, where many duplications and rearrangements followed by subfunctionalization of
cis regulatory elements occurred after they had diverged from other vertebrates⁴⁹,⁵⁰.

Connecting enhancers with their target genes is a major challenge. In contrast to promoters that reside in the
first 1–2 kb upstream of the TSS of a gene, enhancers can be found much far away from the genes they influence.
In *Xenopus tropicalis*, CNS1 is located 16.6 kb upstream of the *hnf1b* TSS and 23.6 kb downstream the end of
the *hnf1b* flanking gene *heatr6* (HEAT repeat containing 6 protein). Synteny is conserved in mouse and CNS1
is located 33.4 kb upstream the *Mus musculus* *Hnf1b* TSS, and 33.4 kb downstream the end of the *Heatr6* gene
raising the question of potential role of CNS1 in *Heatr6* regulation. *Heatr6* is known to be significantly expressed
in the ectoderm and neurectoderm in *Xenopus* embryos, but neither in the developing pronephros (Xenbase,
https://www.xenbase.org), nor in the developing mouse kidney (GUDMAP, https://www.gudmap.org). *Heatr6* is
therefore unlikely to be under the influence of CNS1. In contrast, our data clearly argue in favor of *hnf1b* being a
target gene regulated by CNS1: 1/CNS1 drives reporter expression in a subset of *hnf1b* expression territories 2/
CNS1 is responsive to pax8, which is known to control *hnf1b* expression in the *Xenopus pronephros*¹⁹ 3/CRISPR/
cas9 which allowed CNS1 editing within its native genomic context leads to a decrease of endogenous *hnf1b* expression. This last observation is striking since several lines of evidence suggest that there is ample enhancer redundancy among developmentally regulated genes. Enhancers regulating the same gene often display overlapping spatiotemporal activity, a mechanism that contribute to fine-tune both the levels and patterns of gene expression but also safeguard against genetic and environmental perturbations therefore conferring robustness to developmental processes⁵¹⁻⁵³. Beside pronephros expression, we found that CNS1 can promote expression of a reporter gene in domains devoid of endogenous *hnf1b* expression, such as somites (Fig. S3). This may indicate the existence of silencer elements that downregulate *hnf1b* expression in these territories. Alternatively, this ectopic expression could be a consequence of the disruption of the CNS1 genomic architecture in our assay since it is known that the genomic context plays an important role for the proper function of enhancer elements.

Enhancers bind multiple TFs to activate specific gene transcription. In the present study, we addressed the role of the AATTCAGGCAATTATCTGCAT sequence identified as a Pax8 binding site present in CNS1. Our observations allowed us to conclude that this binding site is essential for CNS1 to activate transcription in specific cellular contexts. In renal cell lines, mutation of the Pax8 binding site leads to a small decrease of CNS1 activity in MDCK cells but does not have any significant effect in IMCD3 cells. In the *Xenopus* embryo, the Pax8 binding site is required for CNS1 activity in the developing pronephros at tailbud stage when patterning of the anlage occurs. Later on, when the segmented tubule is formed, the Pax8 binding site is not required anymore for reporter gene expression, indicating that other TFs are also at work and can compensate for its mutation. Interestingly, we noticed that CNS1 harbors LEF/TCF (T-cell factor/lymphoid enhancer-binding factor) binding sites suggesting that the Wnt/beta-catenin pathway could be an input for CNS1 activity. In agreement with this hypothesis, Wnt signaling components are expressed at several stages of pronephros development (including Wnt9a, Frizzle8 and Frizzled 7⁵⁴), canonical Wnt signaling activity is detected in the pronephros of the tadpole in a *Xenopus* Wnt-responsive transgenic line⁵⁵ and Wnt/beta-catenin signaling is required for pronephros development⁵⁶. Bioinformatic analysis also identified a conserved putative binding site for RBPJk (Recombination Signal-Binding Protein 1for J-Kappa), a well-known component of the Notch pathway. Several notch ligands and notch effectors are expressed in the developing pronephros and proper notch activity is necessary for the establishment of proximal cells and differentiation of intermediate and distal tubules⁵⁷⁻⁶⁰. Human and *Xenopus* CNS1 also harbor a putative GATA3 binding site. Although GATA3 expression in the pronephros is restricted to the distal and collecting tubules at tadpole stage, it may play a role in controlling CNS1 activity in these regions. To our knowledge, the effect of GATA3 loss of function on pronephros development has not been studied yet.

Because of their high degree of homology, both within the paired domain and in their COOH-terminal transactivation domain, Pax-2 and Pax-8 proteins share similar DNA binding motifs⁶¹. Indeed, we have found that exogenous pax2 is able to activate CNS1 transcriptional activity as pax8 does, while this activation is lost when mutations are introduced in the AATTCAGGCAATTATCTGCAT consensus sequence. However, our results are not in favor of a predominant role of pax2 in CNS1 pronephric activity in the *Xenopus* embryo. Pax2 knockdown neither affects reporter gene expression driven by CNS1 nor *hnf1b* endogenous expression in the pronephros at tailbud stage. Nevertheless, we cannot rule out the possibility that pax2 may contribute to pronephric CNS1 activity in vivo; but that its absence may be compensated by pax8. Pax2 and Pax8 functional redundancy is well documented in different tissues including kidney⁶². During mouse kidney development, both genes act in a cooperative and dosage dependent manner as revealed by the correlation between the number of Pax2/8 alleles and the severity of renal defects⁶³⁻⁶⁵.

Our previous work has highlighted a role for pax8 in *hnf1* expression at the early steps of pronephros development¹⁹. To the best of our knowledge, among the studies aiming at the identification of Pax8 and Pax2 target genes in different tissues, *Hnf1b* was not identified as direct target gene. With this work we demonstrate for the first time that *hnf1b* is directly regulated by pax8 in the developing pronephros. Whether this regulation is functionally conserved during kidney development in other vertebrates remains to be proved. In zebrafish, *pax2a/8*-deficient embryos fail to activate expression of *hnf1ba* in the intermediate mesoderm²². In mouse, *Hnf1b* is co-expressed with *Pax2* or *Pax8* in several renal compartments of the developing kidney. *Hnf1b* has been shown to cooperate with *Pax2* to control common pathways governing both kidney morphogenesis and ureter differentiation⁶⁶ but a control of *Hnf1b* expression by either *Pax8* or *Pax2* in the developing mouse kidney has not yet been documented. *Hnf1b* expression appears normal in *Pax2*-deficient embryos⁶⁶ but since *Pax8* is expressed in these embryos⁶², a definitive conclusion on *Hnf1b* regulation by *Pax* genes cannot be drawn. The finding that in mouse embryonic kidney but not liver or lung, CNS1 is located in open chromatin and associated with active enhancer histone marks strongly supports the hypothesis of a functional role for CNS1 in the *Hnf1b* transcriptional regulation in the developing mouse kidney. *Pax2* and *Pax8* have been shown to recruit activating complexes that imprint distinct patterns of histone methylation resulting in gene activation⁶⁷,⁶⁸. *Pax2* and/or *Pax8* could therefore impart epigenetic changes at the *Hnf1b* locus leading to *Hnf1b* transcriptional regulation.

Variation in gene expression is a significant determinant of human disease and enhancer-related dysregulation of gene expression has been recognized as one of the main drivers in the pathogenesis of many diseases⁶⁹. There are several described Mendelian disorders for which specific mutations in enhancers disrupt enhancer–gene regulation and are unequivocally causal⁷⁰. To date, all described *HNF1B* mutations in humans are located in the HNF1B coding region or splice sites. The identification of *HNF1B* regulatory sequences and their binding TFs could contribute to the understanding of *HNF1B*-associated diseases and may open important perspectives into the diagnosis for understanding congenital anomalies of kidney and urinary tract in human.
**Methods**

**Identification of conserved non-coding regions (CNSs).** We used the dynamic graphical interface ECR Browser (http://ecrbrowser.dcode.org/), a tool for visualizing and accessing data from comparisons of multiple vertebrate genomes²⁴. With minimum similarity of 70% and 100-base pair length as set parameters, sequences were filtered by conservation analysis across different vertebrate species (Mus musculus, Homo sapiens, Rattus norvegicus, Monodelphis domestica, Gallus gallus, Xenopus tropicalis, Danio rerio) using a mouse 157 kb genomic sequence encompassing *Hnf1* gene as a reference genome. R-Vista 2.0 (https://rvista.dcode.org) was used to identify conserved putative TFs binding sites in CNS1³⁰. P300 and histones ChIP-seq, DNase-seq and ATAC-seq peaks were displayed using the ENCODE databases (https://www.encodeproject.org)²⁵,²⁶ and database from the Veenstra's laboratory (http://veenstra.science.ru.nl/trackhubx.htm)²⁷ for *Xenopus*. General analyses and visualization were done using UCSC genome browser (https://genome.ucsc.edu/) and Ensembl genome browser (http://www.ensembl.org/).

**Electrophoretic mobility shift assays.** Electrophoretic mobility shift assay (EMSA) analysis was performed using the Light Shift® Chemiluminescent EMSA Kit (Thermo Scientific). Briefly, nuclear extracts were prepared from MDCK cells, while DNA-probes (Pax8-BS, Control TPO, Pax8-BS mut) were biotin-labelled at 5' End (Eurofins). Twenty µG nuclear extracts were incubated with the different biotin-labeled probes on ice for 20 min in a 20 µL mixture. For competition EMSA, a 10- or 100-fold molar excess of the cold probes was added to the EMSA binding reaction. For supershift assays, 1 µg of anti-Pax8 antibody (Abcam ab53490), was incubated with the nuclear extracts for additional 10 min at 4 °C before adding the probe. The reaction mixture was run on 6% polyacrylamide gel in 0.5× Tris–Borate-EDTA buffer (pH 8.5) at 300 V. The gel was electro-blotted to positively charged nylon membrane using a UV-light cross-linker instrument equipped with 254 nm bulbs at 120 mJ/cm². Biotin-labeled DNA was detected using streptavidin-horseradish conjugate and the chemiluminescent substrate contained in LightShift Chemiluminescent EMSA Kit.

**Cell culture, transfection and luciferase reporter assay.** HEK293, MDCK, and IMCD3 cells were cultured in Dulbecco's modified Eagle's medium (DMEM) supplemented with 10% fetal calf serum. The reporter constructs used for luciferase assay were generated as the following: CNS1 was amplified by PCR from *Xenopus laevis* genomic DNA using primers 5'-CCGCTCGAGCAGAGCAGACAGGGTCTGTA-3' and 5'-CCCAAGCTTTGACCGTCAGTTTCATGACT-3' and inserted into pGEM*-T Easy vectors (PROMEGA). Then, CNS1 was subcloned into the pGL4.23[luc2]miniP (PROMEGA) using XhoI and HindIII to generate CNS1-Luc. Comitantly, we introduced several point mutations variants in the Pax8 binding site using the QuikChange® II XL Site-Directed Mutagenesis Kit (Agilent Technologies) to generate the CNS1 mut-Luc vector. Primers used are 5'-CCCCTAACATACGTTGAATTAAAAATTAGGGTCCTAACCTGCATGGCTTCCCTCTGATTAAG-3' and 5'-CTTAATCAGAGGGAAGCCATGCAGGTTAGGACCTAATTTTTAATTCACGTATGTTAGGGG-3'. The CNS1 and the mutated CNS1 nucleotide sequences were confirmed by sequencing. For Dual Luciferase Reporter Assay analyses, cells were plated into 12-well plates at a density of 2.5 × 10⁵ cells/well. Transfection mixtures containing Opti-MEM® Reduced Serum Media (ThermoFisher), 500 ng of reporter vector (pGL4.23[luc2]miniP called Empty-Luc, CNS1-Luc or CNS1mut-Luc), 5 ng of pCMV-renilla luciferase normalization vector and, when appropriate, 50 pg or 100 pg of the expression vector (pCS2pax2, pCS2pax8VP16, pCS2pax8ΔO, Caroll and Vize, 1999, Buisson et al.¹⁹) were prepared. Cells were transfected using 1.5 µL of X-tremeGENE HP DNA Transfection Reagent (Roche) according to the manufacturer's protocol. Transfected cells were washed twice in PBS, followed by the addition of 200 µL 1× passive lysis buffer (Promega, Madison, WI, USA). All values are shown as the mean ± S.E.M. All experiments were repeated at least three times. Statistical significance was determined using One-Way ANOVA followed by Tuckey's multiple comparison test.

**Immunofluorescence.** Cells on a glass slide were fixed with 3% formaldehyde at room temperature for 15 min and permeabilized with 0.2%. Triton X-100 in PBS for 5 min. Cells were then treated with 3% bovine serum albumin in PBS and incubated overnight at 4 °C with the primary mouse antibody anti-Pax8 (Abcam ab53490, diluted 1/100). After washing unbound antibodies, cells were incubated with the secondary antibodies Alexa Fluor 488 anti-mouse (1/2000) along with DAPI for 1 h. Slides were washed and mounted in Mowiol mounting medium (Sigma-Aldrich). Pictures were taken using Zeiss Apotome Wide Field Microscope.

**Xenopus embryos, morpholino microinjection, in situ hybridization.** *Xenopus laevis* and *Xenopus tropicalis* were purchased from the CNRS *Xenopus* breeding Center (Rennes, France). Embryos were obtained as described¹⁹,⁷¹, and were raised in Modified Barth's solution (MBS). Stages were according to the normal table of *Xenopus laevis*⁷² and as described for *Xenopus tropicalis*⁷³. All experimental protocols and the handling of *Xenopus* were done in accordance with the European Community Directive 2010/63/UE and with the ARRIVE guidelines. They were approved by the French National Ethics Committee for Science and Health report on "Ethical Principles for Animal Experimentation"—C2EA-05 Charles Darwin—under Agreement No 27754-2021050515059177. Microinjection of morpholinos was performed according to Buisson et al.¹⁹. Pax2 knock-down was obtained by injection of MoPax2.2; Pax8 knockdown was obtained by co-injection of MoPax8 and MoPax8.B¹⁹. Whole mount in situ hybridization to detect *GFP*, *hnf1b* and *lhx1* mRNAs were performed as in Ref.¹⁹.

**Transgenic reporter assay in *Xenopus laevis*.** Transgenic *Xenopus laevis* embryos were generated using the I-SceI meganuclease method⁴⁰ or the restriction-enzyme-mediated integration (REMI) method⁴¹. Constructs
used for transgenesis were obtained as follows. CNS1 or mutated CNS1 sequences were retrieved from CNS1-Luc and CNS1 mut-Luc plasmids using XhoI and Hind III, and were inserted upstream of a human β-globin basal promoter (–37 to +12) followed by an eGFP coding sequence into the beta globin-eGFP reporter vector described in Ogino et al.²³. Manipulated embryos were cultured until the desired stage, and normally developing embryos were subjected to whole-mount in situ hybridization for GFP mRNA detection.

**CRISPR/Cas9 editing in Xenopus tropicalis.** gRNAs were designed and selected for their high predicted specificity and efficiency using CRISPOR online tool⁷⁴ (http://crispor.tefor.net/). gRNAa (ACTGTGCTCAGC TTAATCAGCGG) and gRNAb (TATCAGGCCACTGAGAAAAGGGG) were designed to target CNS1. The *slc45a2* gene was targeted with the two gRNAs GGAGCCTCCCAGGAGATCTACGG and ATCCCACGCTGA GCAGAACAGGG.

They were purchased from Integrated DNA Technologies, as Alt-R crRNA and tracrRNAv (IDT, Coralville, IA, USA) and dissolved in duplex buffer (IDT) at 100 µM each. cr:tracrRNA duplexes were obtained by mixing equal amount of crRNA and tracrRNA, heating at 95 °C for five minutes and letting cool down to room temperature. GRNA:Cas9 RNP complex was obtained by incubating 1 µL of 30 µM Cas9 protein (kindly provided by TACGENE, Paris, France) with 2 µL cr:tracrRNA duplex in a final volume of 10 µL of 20 mM Hepes-NaOH, 150 mM KCl, pH 7.5 for 10 min at 28 °C.

**Xenopus tropicalis** one-cell stage embryos were injected with 2 nL of gRNA:Cas9 RNP complex solution and were cultured to the desired stage. For co-injection, both complexes were mixed equally. Single embryo genomic DNA was obtained by digesting for 1 h at 55 °C in 100 µL lysis buffer (100 mM Tris-HCl pH 7.5, 1 mM EDTA, 250 mM NaCl, 0.2% SDS, 0.1 µg/µL Proteinase K), precipitating with 1 volume of isopropanol and resuspended in 100 µL PCR-grade water. The region surrounding the gRNA binding site was amplified by PCR using 5'-CCT TACGCCTTACTGAAGAGTGC-3' as forward primer and 5'-CAGCCAGGGAACTATGTGCAA-3' as reverse primer and PCR products were analyzed by 1.5% agarose TAE gel electrophoresis and send for sequencing. TIDE was used to determine total efficiency and insertion and deletion frequencies in the amplified gene region⁷⁵.

**Data availability**

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request. Public databases and bioinformatic tools used for our studies are the following: ECR Browser (http://ecrbrowser.dcode.org/)²⁴ for comparisons of multiple vertebrate genomes; R-Vista 2.0 (https://rvista.dcode.org)³⁰ to identify conserved putative TFs binding sites; ENCODE databases (https://www.encodeproject.org)²⁵,²⁶ to display histones ChIP-seq, DNase-seq and ATAC-seq peaks in mammals and the database from the Veenstra's laboratory for Xenopus (http://veenstra.science.ru.nl/trackhubx.htm)²⁷. Accessible chromatin and ChIP-seq data are visualized with the UCSC genome browser at (https://genome.ucsc.edu/)²⁷ and Ensembl genome browser (http://www.ensembl.org/)²⁷. Sequence alignment was performed with the short sequence aligner software ClustalW (http://www.ebi.ac.uk/Tools/msa/clustalw2/)²⁸. Gene expression data are available on (https://sckidney flatironinstitute.org)²⁹, GUDMAP (https://www.gudmap.org)⁸⁰,⁸¹ or Xenbase (https://www.xenbase.org).

Received: 24 May 2022; Accepted: 23 September 2022

Published online: 19 November 2022

**References**

1. Desgrange, A. et al. HNF1B controls epithelial organization and cell polarity during uretic bud branching and collecting duct morphogenesis. *Development* **144**, 4704–4719 (2017).

2. Gresh, L. et al. A transcriptional network in polycystic kidney disease. *EMBO J.* **23**, 1657–1668 (2004).

3. Heliot, C. et al. HNF1B controls proximal-intermediate nephron segment identity in vertebrates by regulating notch signalling components and Irx1/2. *Development* **140**, 873–885 (2013).

4. Lokmane, L., Heliot, C., Garcia-Villalba, P., Fabre, M. & Cereghini, S. vHNF1 functions in distinct regulatory circuits to control uretic bud branching and early nephrogenesis. *Development* **137**, 347–357 (2010).

5. Massa, F. et al. Hepatocyte nuclear factor 1β controls nephron tubular development. *Development* **140**, 886–896 (2013).

6. Kaminski, M. M. et al. Direct reprogramming of fibroblasts into renal tubular epithelial cells by defined transcription factors. *Nat. Cell Biol.* **18**, 1269–1280 (2016).

7. Chan, S. C. et al. Mechanism of fibrosis in HNF1B-related autosomal dominant tubulointerstitial kidney disease. *J. Am. Soc. Nephrol.* **29**, 2493–2509 (2018).

8. Niborski, L. L. et al. Hnf1b haploinsufficiency differentially affects developmental target genes in a new renal cysts and diabetes mouse model. *Dis. Model Mech.* **14**, 047498 (2021).

9. Verdeguer, F. et al. A mitotic transcriptional switch in polycystic kidney disease. *Nat. Med.* **16**, 106–110 (2010).

10. Pouilhe, M., Gilardi-Hebenstreit, P., Desmarquet-Trin Dinh, C. & Charnay, P. Direct regulation of vHnf1 by retinoic acid signaling and MAF-related factors in the neural tube. *Dev. Biol.* **309**, 344–357 (2007).

11. Clissold, R. L., Hamilton, A. J., Hattersley, A. T., Ellard, S. & Bingham, C. HNF1B-associated renal and extra-renal disease-an expanding clinical spectrum. *Nat. Rev. Nephrol.* **11**, 102–112 (2015).

12. Bellanné-Chantelot, C. et al. Large genomic rearrangements in the hepatocyte nuclear factor-1beta (TCF2) gene are the most frequent cause of maturity-onset diabetes of the young type 5. *Diabetes* **54**, 3126–3132 (2005).

13. Edghill, E. L., Bingham, C., Ellard, S. & Hattersley, A. T. Mutations in hepatocyte nuclear factor-1beta and their related phenotypes. *J. Med. Genet.* **43**, 84–90 (2006).

14. Heidet, L. et al. Spectrum of HNF1B mutations in a large cohort of patients who harbor renal diseases. *Clin. J. Am. Soc. Nephrol.* **5**, 1079–1090 (2010).

15. Blackburn, A. T. M. & Miller, R. K. Modeling congenital kidney diseases in *Xenopus laevis*. *Dis. Model Mech.* **12**, 038604 (2019).

16. Getwan, M. & Lienkamp, S. S. Toolbox in a tadpole: Xenopus for kidney research. *Cell Tissue Res.* **369**, 143–157 (2017).

17. Raciti, D. et al. Organization of the pronephric kidney revealed by large-scale gene expression mapping. *Genome Biol.* **9**, R84 (2008).
18. Zhou, X. & Vize, P. D. Proximo-distal specialization of epithelial transport processes within the Xenopus pronephric kidney tubules. *Dev. Biol.* **271**, 322–338 (2004).

19. Buisson, I., Le Bouffant, R., Futel, M., Riou, J.-F. & Umbhauer, M. Pax8 and Pax2 are specifically required at different steps of Xenopus pronephros development. *Dev. Biol.* **397**, 175–190 (2015).

20. Carroll, T. J. & Vize, P. D. Synergism between Pax-8 and lim-1 in embryonic kidney development. *Dev. Biol.* **214**, 46–59 (1999).

21. del Viso, F., Bhattacharya, D., Kong, Y., Gilchrist, M. J. & Khokha, M. K. Exon capture and bulk segregant analysis: Rapid discovery of causative mutations using high-throughput sequencing. *BMC Genomics* **13**, 649 (2012).

22. Naylor, R. W., Przepierski, A., Ren, Q., Yu, J. & Davidson, A. J. HNF1β is essential for nephron segmentation during nephrogenesis. *J. Am. Soc. Nephrol.* **24**, 77–87 (2013).

23. Ogino, H., Ochi, H., Uchiyama, C., Louie, S. & Grainger, R. M. Comparative genomics-based identification and analysis of cis-regulatory elements. *Methods Mol. Biol.* **917**, 245–263 (2012).

24. Ovcharenko, I., Nobrega, M. A., Loots, G. G. & Stubbs, L. ECR Browser: A tool for visualizing and accessing data from comparisons of multiple vertebrate genomes. *Nucleic Acids Res.* **32**, W280–W286 (2004).

25. ENCODE Project Consortium. An integrated encyclopedia of DNA elements in the human genome. *Nature* **489**, 57–74 (2012).

26. Luo, Y. *et al.* New developments on the encyclopedia of DNA elements (ENCODE) data portal. *Nucleic Acids Res.* **48**, D882–D889 (2020).

27. Hontelez, S. *et al.* Embryonic transcription is controlled by maternally defined chromatin state. *Nat. Commun.* **6**, 10148 (2015).

28. Barbacci, E. *et al.* Variant hepatocyte nuclear factor 1 is required for visceral endoderm specification. *Development* **126**, 4795–4805 (1999).

29. Coffinier, C., Barra, J., Babinet, C. & Yaniv, M. Expression of the vHNF1/HNF1beta homeoprotein gene during mouse organogenesis. *Mech. Dev.* **89**, 211–213 (1999).

30. Loots, G. G. & Ovcharenko, I. rVISTA 2.0: Evolutionary analysis of transcription factor binding sites. *Nucleic Acids Res.* **32**, W217–W221 (2004).

31. Ruiz-Llorente, S. *et al.* Genome-wide analysis of Pax8 binding provides new insights into thyroid functions. *BMC Genomics* **13**, 147 (2012).

32. Zannini, M., Francis-Lang, H., Plachov, D. & Di Lauro, R. Pax-8, a paired domain-containing protein, binds to a sequence overlapping the recognition site of a homeodomain and activates transcription from two thyroid-specific promoters. *Mol. Cell Biol.* **12**, 4230–4241 (1992).

33. Khundmiri, S. J., Chen, L., Lederer, E. D., Yang, C.-R. & Knepper, M. A. Transcriptomes of major proximal tubule cell culture models. *J. Am. Soc. Nephrol.* **32**, 86–97 (2021).

34. Esposito, C., Miccadei, S., Saiardi, A. & Civitareale, D. PAX 8 activates the enhancer of the human thyroperoxidase gene. *Biochem. J.* **331**, 37–40 (1998).

35. Eberhard, D., Jiménez, G., Heavey, B. & Busslinger, M. Transcriptional repression by Pax5 (BSAP) through interaction with corepressors of the Groucho family. *EMBO J.* **19**, 2292–2303 (2000).

36. Lechner, M. S. & Dressler, G. R. Mapping of Pax-2 transcription activation domains. *J. Biol. Chem.* **271**, 21088–21093 (1996).

37. Heller, N. & Brändli, A. W. Xenopus Pax-2/5/8 orthologues: Novel insights into Pax gene evolution and identification of Pax-8 as the earliest marker for otic and pronephric cell lineages. *Dev. Genet.* **24**, 208–219 (1999).

38. Ma, Z. *et al.* Mutations of HNF-1beta inhibit epithelial morphogenesis through dysregulation of SOCS-3. *Proc. Natl. Acad. Sci. U.S.A.* **104**, 20386–20391 (2007).

39. Ogino, H. & Ochi, H. Resources and transgenesis techniques for functional genomics in Xenopus. *Dev. Growth Differ.* **51**, 387–401 (2009).

40. Pan, F. C., Chen, Y., Loeber, J., Henningfeld, K. & Pieler, T. I-Scl meganuclease-mediated transgenesis in Xenopus. *Dev. Dyn.* **235**, 247–252 (2006).

41. Kroll, K. L. & Amaya, E. Transgenic Xenopus embryos from sperm nuclear transplantations reveal FGF signaling requirements during gastrulation. *Development* **122**, 3173–3183 (1996).

42. Shigeta, M. *et al.* Rapid and efficient analysis of gene function using CRISPR-Cas9 in *Xenopus tropicalis* founders. *Genes Cells* **21**, 755–771 (2016).

43. DeLay, B. D. *et al.* Tissue-specific gene inactivation in *Xenopus laevis*: Knockout of lhx1 in the kidney with CRISPR/Cas9. *Genetics* **208**, 673–686 (2018).

44. Allende, M. L., Manzanares, M., Tena, J. J., Feijóo, C. G. & Gómez-Skarmeta, J. L. Cracking the genome's second code: Enhancer detection by combined phylogenetic footprinting and transgenic fish and frog embryos. *Methods* **39**, 212–219 (2006).

45. Fu, S. *et al.* Differential analysis of chromatin accessibility and histone modifications for predicting mouse developmental enhancers. *Nucleic Acids Res.* **46**, 11184–11201 (2018).

46. Gao, T. & Qian, J. EnhancerAtlas 2.0: An updated resource with enhancer annotation in 586 tissue/cell types across nine species. *Nucleic Acids Res.* **48**, D58–D64 (2020).

47. Meireles-Filho, A. C. A. & Stark, A. Comparative genomics of gene regulation-conservation and divergence of cis-regulatory information. *Curr. Opin. Genet. Dev.* **19**, 565–570 (2009).

48. Wong, E. S. *et al.* Deep conservation of the enhancer regulatory code in animals. *Science* **370**, 8137 (2020).

49. Cañestro, C., Yokoi, H. & Postlethwait, J. H. Evolutionary developmental biology and genomics. *Nat. Rev. Genet.* **8**, 932–942 (2007).

50. Hellsten, U. *et al.* The genome of the Western clawed frog *Xenopus tropicalis*. *Science* **328**, 633–636 (2010).

51. Hong, J.-W., Hendrix, D. A. & Levine, M. S. Shadow enhancers as a source of evolutionary novelty. *Science* **321**, 1314 (2008).

52. Kvon, E. Z., Waymack, R., Gad, M. & Wunderlich, Z. Enhancer redundancy in development and disease. *Nat. Rev. Genet.* **22**, 324–336 (2021).

53. Osterwalder, M. *et al.* Enhancer redundancy provides phenotypic robustness in mammalian development. *Nature* **554**, 239–243 (2018).

54. Zhang, B., Tran, U. & Wessely, O. Expression of Wnt signaling components during Xenopus pronephros development. *PLoS ONE* **6**, e26533 (2011).

55. Denayer, T., Van Roy, F. & Vleminkx, K. In vivo tracing of canonical Wnt signaling in Xenopus tadpoles by means of an inducible transgenic reporter tool. *FEBS Lett.* **580**, 393–398 (2006).

56. Lyons, J. P. *et al.* Requirement of Wnt/beta-catenin signaling in pronephric kidney development. *Mech. Dev.* **126**, 142–159 (2009).

57. Katada, T. & Sakurai, H. Proper notch activity is necessary for the establishment of proximal cells and differentiation of intermediate, distal, and connecting tubule in *Xenopus pronephros* development. *Dev. Dyn.* **245**, 472–482 (2016).

58. Naylor, R. W. & Jones, E. A. Notch activates Wnt-4 signalling to control medio-lateral patterning of the pronephros. *Development* **136**, 3585–3595 (2009).

59. Taelman, V., Van Campenhout, C., Sölter, M., Pieler, T. & Bellefroid, E. J. The notch-effector HRT1 gene plays a role in glomerular development and patterning of the *Xenopus pronephros* anlagen. *Development* **133**, 2961–2971 (2006).

60. McLaughlin, K. A., Rones, M. S. & Mercola, M. Notch regulates cell fate in the developing pronephros. *Dev. Biol.* **227**, 567–580 (2000).

61. Phelps, D. E. & Dressler, G. R. Identification of novel Pax-2 binding sites by chromatin precipitation. *J. Biol. Chem.* **271**, 7978–7985 (1996).
62. Bouchard, M., Souabni, A., Mandler, M., Neubüser, A. & Busslinger, M. Nephric lineage specification by Pax2 and Pax8. *Genes Dev.* **16**, 2958–2970 (2002).

63. Marcotte, M., Sharma, R. & Bouchard, M. Gene regulatory network of renal primordium development. *Pediatr. Nephrol.* **29**, 637–644 (2014).

64. Narlis, M., Grote, D., Gaitan, Y., Boualia, S. K. & Bouchard, M. Pax2 and pax8 regulate branching morphogenesis and nephron differentiation in the developing kidney. *J. Am. Soc. Nephrol.* **18**, 1121–1129 (2007).

65. Sharma, R., Sanchez-Ferras, O. & Bouchard, M. Pax genes in renal development, disease and regeneration. *Semin. Cell Dev. Biol.* **44**, 97–106 (2015).

66. Paces-Fessy, M., Fabre, M., Lesaulnier, C. & Cereghini, S. Hnf1b and Pax2 cooperate to control different pathways in kidney and ureter morphogenesis. *Hum. Mol. Genet.* **21**, 3143–3155 (2012).

67. Abraham, S. et al. The Groucho-associated phosphatase PPM1B displaces Pax transactivation domain interacting protein (PTIP) to switch the transcription factor Pax2 from a transcriptional activator to a repressor. *J. Biol. Chem.* **290**, 7185–7194 (2015).

68. Patel, S. R., Bhumbra, S. S., Paknikar, R. S. & Dressler, G. R. Epigenetic mechanisms of Groucho/Grg/TLE mediated transcriptional repression. *Mol. Cell* **45**, 185–195 (2012).

69. Chatterjee, S. & Ahituv, N. Gene regulatory elements, major drivers of human disease. *Annu. Rev. Genomics Hum. Genet.* **18**, 45–63 (2017).

70. Karnuta, J. M. & Scacheri, P. C. Enhancers: Bridging the gap between gene control and human disease. *Hum Mol. Genet.* **27**, R219–R227 (2018).

71. Jafkins, A., Abu-Daya, A., Noble, A., Zimmerman, L. B. & Guille, M. Husbandry of *Xenopus tropicalis*. *Methods Mol. Biol.* **917**, 17–31 (2012).

72. Nieuwkoop, P. D. & Faber, J. *Normal Table of Xenopus lævis* (Daudin) 2nd edn. (North-Holland Publishing Company, 1967).

73. Khokha, M. K. et al. Techniques and probes for the study of *Xenopus tropicalis* development. *Dev. Dyn.* **225**, 499–510 (2002).

74. Concordet, J.-P. & Haeussler, M. CRISPOR: Intuitive guide selection for CRISPR/Cas9 genome editing experiments and screens. *Nucleic Acids Res.* **46**, W242–W245 (2018).

75. Brinkman, E. K., Chen, T., Amendola, M. & van Steensel, B. Easy quantitative assessment of genome editing by sequence trace decomposition. *Nucleic Acids Res.* **42**, e168 (2014).

76. Kent, W. J. et al. The human genome browser at UCSC. *Genome Res.* **12**, 996–1006 (2002).

77. Howe, K. L. et al. Ensembl 2021. *Nucleic Acids Res.* **49**, D884–D891 (2021).

78. Madeira, F. et al. Search and sequence analysis tools services from EMBL-EBI in 2022. *Nucleic Acids Res.* https://doi.org/10.1093/nar/gkac240 (2022).

79. Lindström, N. O. et al. Spatial transcriptional mapping of the human nephrogenic program. *Dev. Cell* **5**, 2381–2398 (2021).

80. Harding, S. D. et al. The GUDMAP database—An online resource for genitourinary research. *Development* **138**, 2845–2853 (2011).

81. McMahon, A. P. et al. GUDMAP: The genitourinary developmental molecular anatomy project. *J. Am. Soc. Nephrol.* **19**, 667–671 (2008).

## Acknowledgements

The authors thank H. Ogino for kindly providing the beta globin-eGFP reporter vector and J.P Concordet at TACGENE (Paris, France) for providing Cas9 protein. They acknowledge the ENCODE consortium and the ENCODE laboratories of Dr. Bing Ren (USCD) and Dr. John Stamatoyannopoulos (UW) for epigenomic data on mouse fetal organs and human cells, respectively. They also acknowledge the laboratory of Dr. G. Veenstra for for ChIP-sequencong data on *Xenopus* embryos and Dr. I. Ovcharenko for help with ECR browser. European Union's Horizon 2020 Research and Innovation Program supported the research presented in this manuscript, under the Marie Skłodowska Curie Grant Agreement No. 642937 (RENALTRACT; MSCA-ITN-2014-642937) and Grants from CNRS and Sorbonne Université. Laura Goea was recipient of a PhD candidate contract from ITN RENALTRACT MSCA-ITN-2014-642937. They thank S. Authier and E. Manzoni for excellent technical assistance in the maintenance of the *Xenopus* animal facility.

## Author contributions

M.U. and J-F.R. conceived and designed the study. M.U. supervised the studies and wrote the paper with S.C. Most of the experiments were performed by L.G. and I.B. except CRISPR/cas9 experiments which were performed by A.E., morpholino experiments performed by V.B. and R.E.M.I. transgenesis by A.C. Results were analyzed by L.G., I.B., M.P.-F. and R.L.B.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-022-21171-x.

**Correspondence** and requests for materials should be addressed to M.U.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
<figure><img src="image_35.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or
format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the
Creative Commons licence, and indicate if changes were made. The images or other third party material in this
article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the
material. If material is not included in the article's Creative Commons licence and your intended use is not
permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from
the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2022