3.1. Characterization of the Newt nZRS Enhancer
We isolated the newt (
C. pyrrhogaster) nZRS fragment (819 bp) from genomic DNA, and a single PCR product was observed (
Figure 1a). The nZRS fragment was sequenced (
Figure 1b), and is available at Genbank with accession number PP691624. Next, we confirmed that the newt nZRS enhancer sequence contained DNA binding sites to drive reporter expression. Therefore, the nZRS sequence was characterized with previous mouse conserved DNA-protein binding sites, conserved uncharacterized sites, and urodele uncharacterized sites (
Figure 1b). MEME Suite was used to produce aligned conserved sequence blocks with other known vertebrate ZRS sequence regions (
Tables S1-S8). Similar to the mouse ZRS, nZRS contains ETS (erythroblast transformation-specific) conserved sites E0, E1 (containing a snake-specific deletion site), E3, and E4 (
Figure 1b;
Tables S1, S3, S7, and S8), but devoid of the E2 site (
Figure S3) [
20,
37]. EST sites are known key regulators of Shh positional expression in the limb [
20,
37]. The nZRS contains two ETV2 sites (ACTTCCTT and AAGGAAGT). ETV2, a master regulator of hematoendothelial lineages, was recently shown to be an upstream regulator of Shh expression, chromatin regulation, and transcriptional activation during limb bud development [
38]. Four Hox sites, WMS (Werner mesomelic syndrome, a point mutation site), and HAND2/Ebox [
39,
40] were also identified in nZRS (
Figure 1b;
Table S6). HAND2 is a known regulator of Shh in the limb bud mesenchyme [
39]. Two Pitx1 sites (GGATTA and TGGTGCGC) that modulate hindlimb expression [
41] were identified in the nZRS sequence. Two vertebrate uncharacterized conserved regions were also identified, an Ebox site and a 32-bp box (
Figure 1b; Tables S3, and S5). Unique uncharacterized regions of nZRS were identified using Pro-Coffee alignment, providing a heat map of regions with low-high homology (
Figure S3). Beyond development, adult newts regulate Shh expression during limb regeneration [
19]. Therefore, we examined potential unique sites that may explain Shh regulation beyond metamorphosis. We identified the following urodele sites: two copies of the TTTTCTTTTTG sequence found in newt ZRS (
C. pyrrhogaster,
C. orientalis, and
P. waltl), an upstream sequence containing an Ebox site (CTACCTT
CATATG), and TCCGAAAAGCCGCGAAGCAACAGAGAGCG (
Figure 1b). Since the ZRS region is sensitive to even single point mutations associated with human disease [
40], it would be interesting to examine if these urodele sites have any functional binding sites, a subject for future research.
3.2. Functional Analysis of the nZRS Enhancer
For functional analysis of the nZRS sequence, we developed the plasmid pnZRS-hmp-mcherry01 (
Figure 1c-d), and used it for F0 transgenic newts (
Table S9). Reporter expression was detected at larval stage 38 in the limb bud (
Figure 2a, b; n=6). To confirm the expression of Shh protein, we developed a newt anti-Shh rabbit polyclonal antibody, and tested for specificity in known Shh-expressing tissue, the retina (
Figure S1). We detected Shh protein localization in newt retinal ganglion cells, as was previously described in fish [
27] (
Figure S1c). In the WT limb bud stage 38.5, Shh protein expression coincided with nZRS reporter expression (
Figure 2c, d). We further validated the nZRS reporter expression at other development stages and in the hindlimb. The newt nZRS reporter was detected in the hindlimb bud and digit patterning (
Figure 3), suggesting that this nZRS reporter possessed spatial and temporal regulation during limb development. Unlike development in the mouse limb, where the fore limb and hind limbs develop almost in synchrony, in urodele amphibians the fore limb with four digits develops first followed by the hind limbs with five digits. During development, newt larval digits are patterned with pre-axial dominance [
2,
22]. Here, the nZRS reporter was observed during forelimb digit formation at stage 52, with reporter expression localized with higher intensity in the fourth digit along the posterior-ventral limb region (
Figure 3b). At stage 52, during hindlimb bud development, a ZPA-like cluster of mCherry cells polarized in the posterior (
Figure 3c). At stage 57, the fifth digit of the hindlimb also expressed mCherry at the posterior-ventral limb region (
Figure 3c). Reporter mCherry+ cells could be tracked above the elbow/knee region (
Figure 3b, c). Transgenic nZRS reporter newts having the phenotype described in
Figure 2 and
Figure 3 were further screened by genomic PCR (
Figure 3d), and to ensure that the full functional transgene region (1.58 kbp) was inserted. Individuals without the functional transgene region, were considered negative individuals and later separated.
3.3. nZRS Reporter Activity Following Metamorphosed Limb Regeneration
Transgenic nZRS reporter newts (described in
Figure 2 and
Figure 3, with positive functional transgene insertions) were reared beyond metamorphosis. When newts were one year old (with a body length 5.5-6 cm), terrestrial metamorphosed limb regeneration was examined, n = 9 (
Figure 4). No mCherry reporter was detected in the intact limb, prior to limb amputations (
Figure 4a). Shh protein expression in the intact limbs of adult WT newts was also not detected at mid stylopod or zeugopod (
Figure S2a, b). Next, we amputated the limb of nZRS transgenic newts at the mid stylopod (humerus) and tracked the mCherry reporter during regeneration. In transgenic nZRS reporter newts, a signal was observed as early as 2 h following amputation, and not detected at 14 days post amputation (dpa) (
Figure 4b). At 35 dpa, nZRS reporter showed mCherry+ cells clustering in the blastema, displaying signs of early ZPA organization (
Figure 4b). Note that the adult WT blastema displayed Shh protein expression at 30 dpa (
Figure S2c). By 40 dpa, the nZRS reporter was localized in the posterior blastema (
Figure 4b).
To rule out autofluorescence, we used non-injected embryos (control, WT) from the same clutch, reared in parallel with the transgenic newts to the same age (
Table S9). Non-injection newts of the same age did not show any reporter signal (
Figure 4c). From 40-49 dpa, there was a 90° posterior tilt between the regenerating blastema and the base of the limb. At the late blastema to early palette stage, at 49 dpa, just before the regeneration of digits, the blastemal ZPA became polarized in the posterior region in preparation for digit regeneration (
Figure 4b, d, e). At 58 dpa, digit regeneration and cartilage condensation were observed, where digits III and IV displayed a reporter signal (
Figure 4f). However, a decreasing signal at 65 dpa was observed in digit IV while a weak signal was detected in digits I-III (
Figure 4f). By 133 dpa, the nZRS reporter signal was not detected, and pigmentation had increased throughout the regenerating digits.
3.4. nZRS Reporter Overlaps with Shh Protein Expression
To identify the nZRS reporter mCherry+ and Shh+ cells observed in the blastemal ZPA (
Figure 4e), we collected limb samples at 65 dpa (n = 3 limbs) (
Figure 4f) and prepared sections for IHC. This day provided immature to mature regenerating tissue along the proximal to distal axis (
Figure 5a). Here, we screened tissue sections to confirm if our nZRS reporter newts had overlapping Shh protein expression in cells.
Sections were stained with anti-Shh rabbit polyclonal and anti-RFP mouse monoclonal antibodies (
Figure 5). At 65 dpa, in transgenic nZRS reporter newts, skeletal condensation and bone maturation were observed in the humerus, ulna, radius, and digit patterning (
Figure 4f, 5a). We examined our nZRS reporter for the mCherry signal in regions distal (
Figure 5c-e) and proximal (
Figure 5f) to the amputation plane at the mid humerus (
Figure 5a). Here, with nZRS reporter newts, we identified and co-localized mCherry+ and Shh+ cells in the dermal region of the regenerating limb (
Figure 5c-f). In the distal region, in a few cells expressing Shh+/mCherry+, the cytoplasm had protruding structures, n= 1/3 (
Figure 5e). mCherry protein expression was identified in the cytoplasm, giving the appearance of a webbed nucleus (
Figure 5e;
Figure 4e). Shh protein expression was co-localized in the cytoplasm and cell membrane (
Figure 5e). Previously, newts expressing mCherry and acetylated tubulin in Schwann cells (SCs) gave a similar webbed nuclear appearance with a flattened morphology [
8]. In the adult WT at 14 and 30 dpa, Shh expression was detected in myelinating Schwann cells (mSCs) of dermal nerves adjacent to glands (
Figure S2c-i). Shh+/mCherry+ cells of the nZRS reporter were commonly observed in group formation along the regenerating dermis (
Figure 5c-f). In proximal regions of the limb where tissue had regenerated, weak expression of mCherry+/Shh- or mCherry+/Shh+ cells was detected, and these cells had a flattened SC-like morphology (
Figure 5f), similar to a previous study [
8]. We reasoned that these distal Shh+/mCherry+ cells were probably immature SCs (iSCs) or dedifferentiated blastemal cells (
Figure 4e) derived from mSCs observed in adult limb regeneration at 14dpa, 30 dpa (
Figure S2) and proximal 65 dpa (
Figure 5f). We therefore applied the nZRS reporter with other known SC markers that could provide evidence for the involvement of SCs in our system (
Figure 6).
3.5. Dermal Regenerating Schwann Cells and Glands Express Shh
We revisited which mSCs were actually contributing to the blastema. Here we co-labeled 30 dpa blastema’s of WT limbs with anti-Shh and anti-Acetylated tubulin, n= 3 adult blastema. In the newt, we previously showed that mSCs express acetylated tubulin [
8]. Here we found that Shh+/acetylated tubulin+ SCs of the dermis, but not acetylated tubulin+ SCs, were participating in the blastema (
Figure 6a-e;
Figure S4). At 30 dpa, Shh+ SCs cells appear to be dedifferentiating, and change morphology (
Figure 6d, e;
Figure S2, S4). We used our nZRS reporter that aligned with Shh protein expression (
Figure 5) to examine another marker of SCs, neural adhesion molecule 1 (NCAM) [
42]. NCAM, also known as CD56, is commonly expressed in immature, sensory, and non-myelinating SCs [
42,
43,
44]. NCAM has also been reported in mesenchymal cells of the newt blastema [
45]. In nZRS reporter newts, iSCs expressed mCherry+/NCAM+, but maturing mCherry+ mSCs with a flattened morphology did not express NCAM at 65 dpa (
Figure 6f-l). At 133 dpa, mSCs with residual mCherry reporter protein expressed acetylated tubulin (
Figure 6m-p), but not Shh, as observed at 65 dpa (
Figure 6f). Dermal glands in the regenerating limb expressed Shh protein (
Figure S2c-i) and the nZRS mCherry reporter (
Figure 6f-p).