{'fraction_non_alphanumeric': 0.03979373749024663, 'fraction_numerical': 0.03192319435492079, 'mean_word_length': 4.976885644768856, 'pattern_counts': {'":': 0, '<': 4, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 2, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
LIM homeodomain transcription factor Isl1 affects urethral epithelium differentiation and apoptosis via Shh

Urethral hypoplasia, including failure of urethral tube closure, is one of the common phenotypes observed in hereditary human disorders, the mechanism of which remains unclear. The present study was thus designed to study the expression, functions, and related mechanisms of the LIM homeobox transcription factor Isl1 throughout mouse urethral development. Results showed that Isl1 was highly expressed in urethral epithelial cells and mesenchymal cells of the genital tubercle (GT). Functional studies were carried out by utilizing the tamoxifen-inducible Isl1-knockout mouse model. Histological and morphological results indicated that Isl1 deletion caused urethral hypoplasia and inhibited maturation of the complex urethral epithelium. In addition, we show that Isl1-deleted mice failed to maintain the progenitor cell population required for renewal of urethral epithelium during tubular morphogenesis and exhibited significantly increased cell death within the urethra. Dual-Luciferase reporter assays and yeast one-hybrid assays showed that ISL1 was essential for normal urethral development by directly targeting the Shh gene. Collectively, results presented here demonstrated that Isl1 plays a crucial role in mouse urethral development, thus increasing our potential for understanding the mechanistic basis of hereditary urethral hypoplasia.


Introduction
External genital development of mammals is regulated by a sophisticated progression of budding and fusion events. In mice, paired genital swellings arise lateral to the cloacal membrane at around embryonic day (E) 10.5, then these swellings merge to form a single genital tubercle (GT) at E11.5, and finally give rise to the prepuce by E13.5 1 . Up to E15.5, male and female GTs are morphologically indistinguishable, and their development is presumably controlled by the same genetic program called androgen-independent GT patterning 2-4 . After E16.5, urethra in males canalize in the presence of androgen signaling, whereas they remain as epithelial cords in females 5,6 . Three germ layers participate in external genital development: endodermally derived urethral plate epithelium forms the entire urethral tube, mesoderm forms the stromal tissue of the genital, and a jacket of ectodermal epithelium forms the skin [7][8][9] . Epithelial differentiation is an important process for reproductive organ development. Urethral defects lead to severe developmental disease; a prime example being the most common urethra defect in human, hypospadias, a failure of urethral tube closure often accompanied by agenesis of the ventral aspect of the genitals, which affects ∼1 in 250 live births 10,11 .
Genetic mutation is one of the crucial causes of urethral defects. Up to now, several factors have been identified that participate in urethral development. Sonic hedgehog (Shh) mutants result in hypospadias associated with increased cell death and failure to maintain Wnt-Fgf8 signaling in endoderm 2,7,[12][13][14][15] . The Fgf8-expressing distal urethral epithelium has been shown to have a growthpromoting function 16,17 . Fgfr2 is required for outgrowth and patterning of the GT, and regulates cell number by controlling the length of specific cell cycle phases 10,18,19 . Hoxa13 mutant hypospadias occur as a result of combined loss of Fgf8 and Bmp7 expression in the urethral epithelium 16 . Genes including Lef1 20,21 , Msx1 1 , bcatenin 22 , and Noggin 6 also have been shown to be involved in urethral development. However, a comprehensive network of gene regulation during urethral development remains unclear.
The LIM homeobox gene Isl1 plays critical roles in multiple tissues in mouse embryonic development, including the nerve 23,24 , stomach 25,26 , limb 27,28 , and heart [29][30][31][32] , by functioning in cell proliferation, apoptosis 33 , and differentiation 34,35 . ISL1 is also involved in hormone biosynthesis and secretion in endocrine tissues such as pancreatic islets 36,37 , pituitary 38 , and pineal glands 39,40 . Isl1 expression has been detected in genital mesenchyme, and Isl1 knockout in mesenchyme results in urogenital malformations, although urethral development in these Isl1 mutants was unaffected 41 . Additionally, genome-wide association studies have revealed that Isl1 is a major susceptibility gene for human congenital anomalies of the kidney and urinary tract (CAKUT) and bladder exstrophy-epispadias complex (BEEC) [42][43][44] , implying potential roles for Isl1 in urethra development. The present study was thus designed to examine a potential role for the transcription factor Isl1 in urethral development.

Mice maintenance and treatment
Adult (6-to 8-week-old) male and female C57BL/6 mice were used for this study. The age of embryos was determined by the appearance of the vaginal plug, which was taken to be E0.5.
Generations of Isl1 MCM/Del and Isl1 F/+ mice were described previously 45 . Briefly, we used a "floxed" Isl1 allele (Isl1F) in which LoxP sites were inserted into the introns flanking exon 4 of the Isl1 locus, and a tamoxifeninducible knock in Isl1-MERCreMER allele. Isl1 F/F mice were mated with Isl1 MCM/+ mice to generate litters with equal numbers of Isl1 MCM/F -inducible knockouts and Isl1 F/+ controls. To induce excision in Isl1 MCM/F embryos, pregnant females were administered an oral gavage of 75 mg/kg body weight of tamoxifen (T5648; Sigma, St. Louis, MO, USA) in corn oil (10 mg/ml) at E9.5 for three consecutive days just before Isl1 expression sharply increased. Embryos were harvested from pregnant mice obtained by timed matings at the desired stages of development and genotyped by common PCR. Sex of embryos were identified by common PCR before E15.5 46 and later morphological analysis of gonads. For studies involving embryos, only males were presented (except for specially marked parts). All animal studies were approved by Ethics Committee of China Agricultural University and performed in accordance with the guidelines and regulatory standards of the Institutional Animal Care and Use of Animals for Scientific Purposes.

Common PCR and qPCR
GTs were dissected from stage-matched embryos and were pooled according to genotype for each litter collected. Total RNA was extracted using Trizol reagent (9109; TaKaRa, Dalian, China). RNA quantity and purity were determined using a NanoDrop (ND-2000, USA). One microgram of high-quality RNA (260/280 ratios slightly higher than 2.0 and 260/230 ratios higher than 1.7 for each pooled sample was reverse transcribed into complementary DNA (cDNA) using M-MLV (M170A; Promega, USA). Quantitative PCR (qPCR) amplication was performed four times (DRR420A; Takara, Dalian, China) in the ABI 7500 system (Applied Biosystems, Foster City, USA). For normalization purposes, an identical set of reactions were prepared for glyceraldehyde 3phosphate dehydrogenase (GAPDH).
Genomic DNA was isolated from tail or GT following the HotSHOT method 47 and genotyping was performed using standard PCR methods with specific primers 45 . The relative location of primers used to identify the wild-type (WT), floxed (Flox), and rearranged alleles are shown as solid arrows in Fig. 2b; primer pair F1/R1 amplifies a 406 bp WT and 502 bp Flox bands. In Isl1 MCM/Del mice, the expected PCR products using the primer pairs MCM-F and MCM-R were 289 bp. Primer pair F1/R2 amplifies a 730-bp rearranged floxed allele in Isl1 MCM/Del genital after Cre activation; there was no product under the same PCR condition before Cre activation. The common PCR (Bio-Rad Laboratories) was performed using the following protocol: 95°C for 5 min; 95°C for 30 s, 60°C for 30 s, 72°C for 45 s (35 cycles); 72°C for 5 min; 4°C holding. PCR products were observed on a 1.5% agarose gel. PCR primers designed for this study were listed in Supplementary Tables S1 and S2.

Western blot
Briefly, total proteins were extracted in RIPA (radio immunoprecipitation assay) buffer (9806; Cell Signaling, Danvers, MA, USA) containing 1 mM phenylmethylsulfonyl fluoride (8553S; Cell Signaling, Danvers, MA, USA) according to the manufacturer's protocol. The BCA Protein Assay Kit (HX18651; Hoaxing, China) was used to measure protein concentration. Electrophoresis was performed with 30 μg total proteins separated by 12% sodium dodecyl sulfate-polyacrylamide gel electrophoresis and transferred to polyvinylidene difluoride membranes (IPVH00010, Millipore, USA). The membrane was blocked with 5% (w/v) nonfat dry milk in 0.05 M Tris-buffered saline and 0.1% Tween-20 (TBST, pH 7.4) for 1 h and incubated with anti-ISL1(1:300; AF1837; R&D, Minneapolis, MN, USA) antibody and internal control GAPDH antibody (1:10,000; Ambion, USA) overnight at 4°C. The secondary antibody, horseradish peroxidase-conjugated donkey anti-goat IgG (1:10,000; ab205723; Abcam, Cambridge, MA, USA), was diluted 1:5000 in TBST. The membranes were visualized using the SuperSignal West Pico Kit (Thermo Scientific, Waltham, MA, USA) substrate at room temperature. We used the ImageJ Software to assay the relative intensity of each blot. The intensity values of each group were normalized to GAPDH (internal control) in the same group.

Hematoxylin and eosin staining
The embryonic specimens were fixed overnight in 4% paraformaldehyde/phosphate-buffered saline (PBS), dehydrated in ethanol, and embedded in paraffin. Sections of 5 μm were cut from wax-embedded embryos and floated onto slides coated with 3-triethoxysilylpropylamine (440140, Sigma). The slides were dried at 37°C overnight, dewaxed through xylene, and then rehydrated through decreasing concentrations of ethanol. Next, they were stained with hematoxylin, counterstained with eosin, dehydrated, and equilibrated with xylene. Sections were photographed under bright-field microscope photograph system (Leica Microsystems, Buffalo Grove, IL, USA).

Immunofluorescence and immunohistochemistry
The sections were deparaffinized, rehydrated, and subjected to microwave antigen retrieval with 0.01 M sodium citrate buffer (pH 6.0). The sections were then blocked with 10% normal donkey or goat serum in PBS at room temperature for 1 h, followed by the incubation with anti-ISL1 antibody (1:50; 40.2D6; Developmental Studies Hybridoma Bank, Iowa City, IA, USA), antibromodeoxyuridine (BrdU) (1:300; G3G4; Developmental Studies Hybridoma Bank) overnight at 4°C. After rinsing thoroughly with PBS, the sections were incubated with the secondary antibodies, which were cy3conjuncted goat anti-mouse IgG (1:50; 115-165-003; Jackson ImmunoResearch, West Grove, PA, USA) for 2 h at room temperature. The sections were then rinsed with PBS and stained with 4′,6-diamidino-2-phenylindole (10236276001; Roche Applied Science, Basel, Switzerland) for 10 min. Finally, 20 μl Vectashield mounting medium (H-1000; Vector Laboratories, Burlingame, CA, USA) was applied to each slide, and a coverslip was sealed in place. Leica Microsystems was used for imaging immunofluorescent sections. IgG was used as the negative control.

Detection of proliferation and apoptotic cells
To detect the proliferation of urethral cells, pregnant mice were injected intraperitoneally with BrdU (100 mg/ kg body weight) 4 h before sacrifice. Embryos were removed and processed for immunofluorescence staining, following the procedure described above. Terminal deoxynucleotidyl transferase dUTP nick-end labeling (TUNEL) assay for the detection of apoptotic cells was performed with the In Situ Apoptosis Detection Kit (Roche Applied Science) according to the manufacturer's instructions. Four embryos of each group were included for the analysis, and six sections of each embryo were examined for the proliferation and apoptosis.

Transient transfection and Dual-Luciferase reporter assays
The Shh, Fgfr2, and Fgf8 promoter fragment (0 to −2000 bp) was cloned from mouse genomic DNA and inserted into the pGL3.0 vector (E1751, Promega). Isl1 was amplified using the primers containing the NheI-XhoI restriction sites and inserted into the pcDNA3.1 vector. An empty luciferase reporter vector was used as a control. The primer sequence is in Supplementary  Table S3.
The human embryonic kidney 293FT cells were maintained in Dulbecco's modified Eagle's medium supplemented with 10% fetal bovine serum and 1% penicillin-streptomycin at 37°C with 5% CO 2 . Cells were transfected with Isl1 expression vector, Fgf8, Fgfr2, Shh luciferase reporter vector, and pTK-Ranilla vector using Lipofectamine 2000 (Invitrogen) according to the manufacturer's instruction. Twenty four hours post transfection, luciferase activity was measured by using Dual-Luciferase Reporter Assay Kit (E1910; Promega) on a Modulus II Microplate Multimode Reader (Turner Biosystems, Sunnyvale, CA, USA). The values were normalized against Renilla luciferase activity. At least four independent experiments were performed.

Yeast one-hybrid assay
The Fgf8, Fgfr2, Shh promoter fragment (0 to −2000 bp) was inserted into the pLacZi vector and Isl1 fragments was inserted into the pB42AD vector as described previously. The primer sequence is in Supplementary Table S4.
Plasmids for pB42AD DNA-binding domain fusions (BD) were co-transformed with the pLacZi:LacZ reporter gene plasmid (AD) into the yeast strain EGY48 using standard transformation techniques. Transformants were grown on proper dropout plates containing X-gal (5-bromo-4chloro-3-indolyl-β-D-galactopyranoside) for color development. Yeast transformation and liquid assay were conducted as described in the Yeast Protocols Handbook (Clontech). At least four independent experiments were performed.

Statistics
The results were expressed as means ± SEM of at least four independent experiments. The differences among groups were determined using a Student's t test or a oneway analysis of variance. P < 0.05 was considered to be statistically significant.

Isl1 expression in the developing mouse urethra
First, we examined Isl1 messenger RNA (mRNA) levels in developing male mouse urethra from E12.5 to adult by using qPCR. Results showed that Isl1 mRNA levels were highest at E12.5, the earliest stage examined, decreasing by 40% at E13.5 and E15.5, and by 62.5% at E18.5 relative to levels observed at E12.5, followed by no significant changes up until adult stages (Fig. 1a). In addition, we investigated ISL1 expression by IHC, and results showed that ISL1 was weakly stained in the outgrowth at E11.5 (Fig. 1d, d1). At E12.5 and E15.5, ISL1 was strongly stained in most urethral epithelium (Fig. 1e, f, e1, f1) and around mesenchymal cells, but by E18.5, only a few ISL1positive cells were observed within the urethral epithelium (Fig. 1g, g1). This Isl1 expression pattern in during urethral development suggests that Isl1 might be involved in regulating urethral development.
Isl1 expression is effectively ablated in Isl1 MCM/Del mice To identify whether Isl1 was involved in regulating urethral development, the tamoxifen-inducible Cre/loxP system was utilized to knockout Isl1 in mouse as in our previous report 25 . Briefly, Isl1-MERCreMER males 45 were used to mate with Isl1 loxP/loxP female mice 48 to produce embryos carrying Isl1 mER-Cre-mER; Isl1 loxP/+ (Fig. 2a), where the Isl1 promoter drives the expression of Cre. We then injected tamoxifen into pregnant dams at E9.5 to generate Isl1-inducible knockout embryos, hereafter referred to as Isl1 MCM/Del . Isl1 F/+ sibling embryos served as controls. Genotype of control and mutant embryos were detected by PCR (Fig. 2c). Western blot results showed that ISL1 protein levels in knockout mice decreased more than 90% compared with control embryos Fig. 1 Isl1 expression pattern during GT urethral development. a Isl1 mRNA relative expression levels at different developmental stages in male GTs. Twenty-one-day mice pituitaries served as a positive control. Each sample contained at least four independent replications. Data were presented as means ± SEM. b, c Abridged general view of d-g1 slice position. Red line represents the urethra. d-g1 ISL1 localization in GTs at E11.5, E12.5, E15.5, and E18.5. Scale bars: d-g, 50, μm; d1-g1, 10 μm. GT, genital tubercle; Mes, genital mesenchyme; Ure, urethra; bl, bladder (Fig. 2d). IHC results demonstrated that ISL1 staining was reduced in urethral epithelium of Isl1 MCM/Del mice compared to controls (Fig. 2e-e3). These data demonstrated that Isl1 was effectively down-regulated in urethral epithelium of Isl1 MCM/Del embryos.

Isl1 deletion disrupts urethral internalization
To investigate effects of Isl1 on GT urethral development, morphological and histological examinations were performed at E15.5 and E18.5 for Isl1 MCM/Del embryos and control littermates. Isl1 deletion resulted in hypoplasia of the GT in male and female embryos (Fig. 3), while GTs had developed with a urethral seam along the ventral midline in control male embryos at E15.5 (Fig. 3a). The overall size of mutant GTs was smaller than control GTs. In addition, the proximal urethral opening of mutants was both wider mediolaterally and extended farther proximodistally compared with controls (Fig. 3b). By E18.5, in Isl1 MCM/Del male embryos, the fusion of the GT prepuce failed, which led to an ectopic opening (Fig. 3e, f).
Further histological analyses revealed that at E15.5, the proximal urethral plate of control embryos had separated medially to form a urethral tube (Fig. 3a2), whereas the distal end remained a bilaminar epithelial plate (Fig. 3a1). In Isl1 MCM/Del mutants at the same stage, urethral tube epithelium was not detected proximally (Fig. 3b2), but the urethral remnant and a shallow groove on the ventral surface of the genital distal were observed (Fig. 3b1). By E18.5, Isl1 deletion resulted in the absence of internal urethra of male GTs (Fig. 3, compare e1 with f1, Fig. 3, compare e2 with f2).
Female Isl1 MCM/Del mutants had hypoplasia similar to that observed in males, but effects on prepuce fusion were not as severe as those observed in male mutants either at E15.5 (Fig. 3c, d) or E18.5 (Fig. 3g, h). Isl1 deletion resulted in an open ventral urethra in females (Fig. 3, compare c2 with d2 at E15.5 and compare g2 with h2 at E18.5). Collectively, these data demonstrated that Isl1 was crucial for urethral formation of GTs.

Isl1 affects epithelial stratification and urethral cell differentiation
Formation of the GT epithelium initiates as early as E12.5 49 . In male embryos, histological examination showed that, in control embryos, basal cells were columnar and elongated perpendicular to the basement membrane, intermediate cells were smaller, and apical cells were rounded or squamous. In contrast, urethral cells in male Isl1 MCM/Del embryos were more homogenous, all cells being rounded, with no obvious orientation within the epithelium (Fig. 4a, b). Keratin 14 (K14), a differentiation marker, was not detected either in male controls or mutants (Fig. 4c, d). By E15.5, the GT epithelium in male controls was stratified in 4-6 cell layers, whereas in Isl1 MCM/Del embryos, GT epithelium was significantly thinner, being composed of 2-4 cell layers (Fig. 4e, f). In addition, in controls, K14 was evenly expressed in all urethral epithelial cells, but in Isl1 MCM/Del embryos, K14 was expressed only in urethral epithelium close to the lumen (Fig. 4g, h). Female embryos exhibited

Isl1 knockout enhances apoptosis of urethral epithelial cells
To investigate the causes by which Isl1 ablation resulted in reductions in urethral epithelial layers, we examined urethral epithelial cell apoptosis. Results showed that the proportion of TUNEL-positive cells within epithelial cells drastically increased in Isl1 MCM/Del male embryos (28.3%) compared with those in control embryos (6.3%) at E12.5 (Fig. 5a, b, j). However, no TUNEL-positive cells were detected in urethral epithelium of controls or Isl1 MCM/Del mutants at E15.5 (Fig. 5e, f, j). In addition, we analyzed mRNA levels of the anti-apoptosis gene Bcl-2 and the proapoptosis gene Bax by qPCR analysis of RNA isolated from whole GT tissue. Results showed that Bax/Bcl-2 ratios rose significantly both in E12.5 and E15.5 Isl1 MCM/Del male embryos GTs relative to controls (Fig. 5l, m). Increased apoptosis indicated by the qPCR results at E15.5 were likely owing to inclusion of urethral mesenchymal cells, as we observed increased apoptosis in mesenchymal cells in mutants relative to controls at this stage (Fig. 5e, f).
To monitor cell proliferation, we performed BrdU labeling and immunostaining studies. In E12.5 male embryos, the proportion of BrdU-positive cells within total urethral epithelium was not significantly different between controls and mutants (49.1% versus 50.3%, Fig. 5c, d, k). The overall proportion of BrdU-positive cells significantly decreased in both controls and mutants in going from E12.5 to E15.5 (11.3% and 15.7%, respectively, Fig. 5g, h, k), consistent with a decreased proliferation rate at the later stage 10 .
Furthermore, we analyzed effects of Isl1 expression on cell cycle genes, including Cyclin A1, Cyclin D2, Cyclin E1, Dp1, and E2f1 in GTs of the Isl1 MCM/Del and control male embryos at E12.5 and E15.5. Results showed mRNA levels of these genes exhibited no significant differences between Isl1 MCM/Del and control male embryos (Fig. 5n, o). Similar results were observed for female embryos (data not shown). These results indicated that Isl1 ablation enhanced apoptosis of urethral epithelial cells at E12.5, but did not have significant effects on cell proliferation.

ISL1 maintains urethral epithelium by targeting Shh
To find potential direct downstream targets of ISL1 affecting differentiation and apoptosis of urethral epithelial cells in mutants, we examined expression of genes required for development of urethral epithelium, including Lef1, Bmp7, Msx1, Fgf8, Shh, Fgfr2, and Noggin, by qPCR analyses of mRNA from mutants and controls. Results showed that levels of Shh mRNA were decreased in GTs of Isl1 MCM/Del embryos relative to controls at E12.5 and E15.5, whereas Fgfr2 and Fgf8 mRNA levels were significantly increased. Isl1 ablation did not have significant effects on mRNA levels of Lef1, Bmp7, Msx1 either at E12.5 or E15.5. At E15.5, Noggin mRNA levels were significantly lower in Isl1 MCM/Del embryos when compared with controls (Fig. 6a, b). These data suggested that Shh, Fgf8, and Fgfr2 were potential direct downstream targets of ISL1.
To investigate whether Shh, Fgf8, and/or Fgfr2 were direct target genes of ISL1 regulating urethral epithelium development, we performed Dual-Luciferase reporter and  Expression of urethral epithelium candidate genes were compared by qPCR in normal and knockout mice at E12.5 and E15.5 in male GTs. c Dual-Luciferase reporter assays. Only the Shh promoter exhibited a significant increase in relative luciferase activity in response to ISL1. d Yeast one-hybrid assays. BD-ISL1 activated the LacZ reporter gene driven by the Shh promoter fragments (AD-Shh), but not Fgf8 (AD-Fgf8) and Fgfr2 (AD-Fgfr2) promoter fragments. e, f SHH protein levels were examined by Western blot. g SHH expressed in the urethral epithelium at E12.5. The white dotted line marks the boundary of urethral epithelium. Red staining is SHH, and DAPI nuclear counterstaining (DNA) is blue. Scale bars: 20 μm. Results were presented as means ± SEM. At least four independent experiments for qPCR and Dual-Luciferase reporter assays were used. *P < 0.05, **P < 0.01, and ***P < 0.001 yeast one-hybrid assays. Promoter regions of Shh, Fgf8, and Fgfr2 were cloned, selecting 1-2000 bp upstream of the 5′-untranslated regions as core promoter regions. Bioinformatics analyses revealed potential target sites for ISL1 (ATTA/TAAT) in promoters of all three genes. Results of Dual-Luciferase reporter assays demonstrated that ISL1 overexpression enhanced activity of the Shh promoter-Luciferase reporter~3-fold, while Fgf8 and Fgfr2 promoter-Luciferase reporter activities exhibited no significant differences (Fig. 6c). Similarly, yeast one-hybrid assays showed that the ISL1 only activated the Shh promoter (Fig. 6d). In addition, Western blot results showed that SHH protein levels decreased by 40% in Isl1 MCM/Del embryos (Fig. 6e, f) relative to controls. Additionally, immunofluorescence analyses showed that SHH was expressed in urethral epithelial cells at E12.5, and ISL1 and SHH proteins were co-expressed within cells of the urethral epithelium. SHH protein decreased after Isl1 depletion, which was consistent with qPCR and Western blot data (Fig. 6g). Together, these data supported Shh as a direct target of ISL1, indicating that Isl1 affects urethral epithelium development via Shh.

Discussion
Expression and function of the LIM homeobox protein ISL1 has been extensively studied in the nervous system 23,26,50 , digestive system 24,25 , and endocrine glands 39,40 . Results of the present study demonstrated that Isl1 was expressed in urethral epithelial cells throughout urethral development, and that ablation of Isl1 just prior to GT development in Isl1 MCM/Del mutants led to GT morphological abnormalities and urethral hypoplasia.
Results indicated that Isl1 was highly expressed in urethral epithelial cells and affected urethral formation throughout GT ontogeny. However, there is a report that conditional knockout mutants of Isl1 using Tbx4-Cre showed GT morphological abnormalities, but no significant effect on urethral development 41 . Tbx4-Cre is expressed in mesenchymal cells beginning in early embryogenesis 51 ; thus, Isl1 was specifically knocked out in mesenchymal cells. In the study presented here, we used Isl1-MERCreMER 45 induced with tamoxifen at E9.5 as the development of the urinary tract occurs after E10.5 1 and Isl1-null embryos die between E7.5 and E9.5 45 . In our model, Isl1 was knocked out both in mesenchymal and urethral epithelial cells from E9.5. Our results showed that Isl1 deficiency resulted in GT morphological abnormalities and complete absence of the male urethra at E18.5. In support, our results are in agreement with reports that Isl1 is expressed in mesenchymal cells derived from mesoderm, which subsequently forms the prepuce 41,42 .Thus, ISL1 is crucial for several aspects of urethral formation, as in the nervous system 23,24,50 , digestive system 25,26 , and endocrine glands systems 39,40 .
In addition, results of this study revealed that ISL1 affects urethra formation by influencing cell survival and differentiation of the epithelial cells. Histological examination revealed reduced layers of stratified urethral epithelial cells in mutants at E12.5 and E15.5, with cells closer to the lumen being more elongated in Isl1-knockout mice. In addition, the proportion of TUNEL-positive cells within epithelial cells was over four times higher in Isl1 MCM/Del male embryos than in the control in the early stage of urethra formation. However, Ching et al. 41 reported that conditional knockout of Isl1 in GT mesenchymal cells resulted in reduced apoptosis and had no significant effects on cell proliferation in GT 51 .
Shh plays important roles in regulating formation and development of multiple organs 52-57 , including urethral development 2,58 . Shh is specifically expressed in urethral epithelial cells, and Shh knockout leads to abnormal cell cycle 59 of the urinary tract, increased cell death, and results in urethral hypoplasia 2,58 , similar to phenotypes observed with the Isl1 MCM/Del mice model presented here. These observations suggest that Isl1 may affect GT development by regulating expression of Shh. In support, our results showed that both Isl1 and Shh were expressed in all urethral epithelial cells. In addition, results of Dual-Luciferase reporter and yeast one-hybrid assays demonstrated that the Shh promoter contained consensus binding sites for ISL1, and that ISL1 overexpression enhanced activity of the Shh promoter. These results are consistent with ISL1 influencing urethral epithelium development by regulation of Shh expression.
Urethral hypoplasia is a common birth defect in humans, yet its etiology and pattern of onset are largely unknown. There are many factors that cause urethral hypoplasia, including environment, diet, and heredity factors. Several studies have detected significant changes in Isl1 in human urethral hypoplasia samples compared to the normal sample by sequencing 42,60 . Draaken et al. 42 reported that Isl1 DNA was significantly decreased in BEEC patients, suggesting that Isl1 plays a role in BEEC. Our laboratory studies showed that Isl1 affected stomach pyloric sphincter development. We investigated the role of Isl1 in the urethral sphincter, which may lead to BEEC. In addition, we tested the urethral sphincter and bladder of Isl1 MCM/Del mutants, and no significant abnormalities were found until E18.5 (data not shown). It may be that Isl1 MCM/Del mutants are not suitable for BEEC studies, but Isl1 MCM/Del mutants can be used as an animal model to study diseases related to urethral hypoplasia.
Lastly, we presented the novel mechanism of Isl1 regulating urethral development (Fig. 7). In brief, Isl1 affects urethral epithelium fate by directly targeting Shh promoter, subsequently affects commitment of progenitor cell, and enhance cell death. These findings are important for our understanding of diseases resulting from abnormalities of urethral development.