Adrenergic receptor signaling induced by Klf15, a regulator of regeneration enhancer, promotes kidney reconstruction.

Despite the recent discovery of tissue regeneration enhancers in highly regenerative animals, upstream and downstream genetic programs connected by these enhancers still remain unclear. Here, we performed a genome-wide analysis of enhancers and associated genes in regenerating nephric tubules of Xenopus laevis. Putative enhancers were identified using assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) and H3K27ac chromatin immunoprecipitation sequencing (ChIP-seq) analyses. Their target genes were predicted based on their proximity to enhancers on genomic DNA and consistency of their transcriptome profiles to ATAC-seq/ChIP-seq profiles of the enhancers. Motif enrichment analysis identified the central role of Krüppel-like factors (Klf) in the enhancer. Klf15, a member of the Klf family, directly binds enhancers and stimulates expression of regenerative genes, including adrenoreceptor alpha 1A (adra1a), whereas inhibition of Klf15 activity results in failure of nephric tubule regeneration. Moreover, pharmacological inhibition of Adra1a-signaling suppresses nephric tubule regeneration, while its activation promotes nephric tubule regeneration and restores organ size. These results indicate that Klf15-dependent adrenergic receptor signaling through regeneration enhancers plays a central role in the genetic network for kidney regeneration.

Retinoic acid control of pax8 during renal specification of Xenopus pronephros involves hox and meis3.

Development of the Xenopus pronephros relies on renal precursors grouped at neurula stage into a specific region of dorso-lateral mesoderm called the kidney field. Formation of the kidney field at early neurula stage is dependent on retinoic (RA) signaling acting upstream of renal master transcriptional regulators such as pax8 or lhx1. Although lhx1 might be a direct target of RA-mediated transcriptional activation in the kidney field, how RA controls the emergence of the kidney field remains poorly understood. In order to better understand RA control of renal specification of the kidney field, we have performed a transcriptomic profiling of genes affected by RA disruption in lateral mesoderm explants isolated prior to the emergence of the kidney field and cultured at different time points until early neurula stage. Besides genes directly involved in pronephric development (pax8, lhx1, osr2, mecom), hox (hoxa1, a3, b3, b4, c5 and d1) and the hox co-factor meis3 appear as a prominent group of genes encoding transcription factors (TFs) downstream of RA. Supporting the idea of a role of meis3 in the kidney field, we have observed that meis3 depletion results in a severe inhibition of pax8 expression in the kidney field. Meis3 depletion only marginally affects expression of lhx1 and aldh1a2 suggesting that meis3 principally acts upstream of pax8. Further arguing for a role of meis3 and hox in the control of pax8, expression of a combination of meis3, hoxb4 and pbx1 in animal caps induces pax8 expression, but not that of lhx1. The same combination of TFs is also able to transactivate a previously identified pax8 enhancer, Pax8-CNS1. Mutagenesis of potential PBX-Hox binding motifs present in Pax8-CNS1 further allows to identify two of them that are necessary for transactivation. Finally, we have tested deletions of regulatory sequences in reporter assays with a previously characterized transgene encompassing 36.5 â€‹kb of the X. tropicalis pax8 gene that allows expression of a truncated pax8-GFP fusion protein recapitulating endogenous pax8 expression. This transgene includes three conserved pax8 enhancers, Pax8-CNS1, Pax8-CNS2 and Pax8-CNS3. Deletion of Pax8-CNS1 alone does not affect reporter expression, but deletion of a 3.5 â€‹kb region encompassing Pax8-CNS1 and Pax8-CNS2 results in a severe inhibition of reporter expression both in the otic placode and kidney field domains.

The shh limb enhancer is activated in patterned limb regeneration but not in hypomorphic limb regeneration in Xenopus laevis.

Xenopus young tadpoles regenerate a limb with the anteroposterior (AP) pattern, but metamorphosed froglets regenerate a hypomorphic limb after amputation. The key gene for AP patterning, shh, is expressed in a regenerating limb of the tadpole but not in that of the froglet. Genomic DNA in the shh limb-specific enhancer, MFCS1 (ZRS), is hypermethylated in froglets but hypomethylated in tadpoles: shh expression may be controlled by epigenetic regulation of MFCS1. Is MFCS1 specifically activated for regenerating the AP-patterned limb? We generated transgenic Xenopus laevis lines that visualize the MFCS1 enhancer activity with a GFP reporter. The transgenic tadpoles showed GFP expression in hoxd13-and shh-expressing domains of developing and regenerating limbs, whereas the froglets showed no GFP expression in the regenerating limbs despite having hoxd13 expression. Genome sequence analysis and co-transfection assays using cultured cells revealed that Hoxd13 can activate Xenopus MFCS1. These results suggest that MFCS1 activation correlates with regeneration of AP-patterned limbs and that re-activation of epigenetically inactivated MFCS1 would be crucial to confer the ability to non-regenerative animals for regenerating a properly patterned limb.

Development of a heat-stable alkaline phosphatase reporter system for cis-regulatory analysis and its application to 3D digital imaging of Xenopus embryonic tissues.

Xenopus is one of the essential model systems for studying vertebrate development. However, one drawback of this system is that, because of the opacity of Xenopus embryos, 3D imaging analysis is limited to surface structures, explant cultures, and post-embryonic tadpoles. To develop a technique for 3D tissue/organ imaging in whole Xenopus embryos, we identified optimal conditions for using placental alkaline phosphatase (PLAP) as a transgenic reporter and applied it to the correlative light microscopy and block-face imaging (CoMBI) method for visualization of PLAP-expressing tissues/organs. In embryos whose endogenous alkaline phosphatase activities were heat-inactivated, PLAP staining visualized various tissue-specific enhancer/promoter activities in a manner consistent with green fluorescent protein (GFP) fluorescence. Furthermore, PLAP staining appeared to be more sensitive than GFP fluorescence as a reporter, and the resulting expression patterns were not mosaic, in striking contrast to the mosaic staining pattern of ß-galactosidase expressed from the lacZ gene that was introduced by the same transgenesis method. Owing to efficient penetration of alkaline phosphatase substrates, PLAP activity was detected in deep tissues, such as the developing brain, spinal cord, heart, and somites, by whole-mount staining. The stained embryos were analyzed by the CoMBI method, resulting in the digital reconstruction of 3D images of the PLAP-expressing tissues. These results demonstrate the efficacy of the PLAP reporter system for detecting enhancer/promoter activities driving deep tissue expression and its combination with the CoMBI method as a powerful approach for 3D digital imaging analysis of specific tissue/organ structures in Xenopus embryos.

X-ray micro-computed tomography of Xenopus tadpole reveals changes in brain ventricular morphology during telencephalon regeneration.

Xenopus tadpoles serve as an exceptional model organism for studying post-embryonic development in vertebrates. During post-embryonic development, large-scale changes in tissue morphology, including organ regeneration and metamorphosis, occur at the organ level. However, understanding these processes in a three-dimensional manner remains challenging. In this study, the use of X-ray micro-computed tomography (microCT) for the three-dimensional observation of the soft tissues of Xenopus tadpoles was explored. The findings revealed that major organs, such as the brain, heart, and kidneys, could be visualized with high contrast by phosphotungstic acid staining following fixation with Bouin's solution. Then, the changes in brain shape during telencephalon regeneration were analyzed as the first example of utilizing microCT to study organ regeneration in Xenopus tadpoles, and it was found that the size of the amputated telencephalon recovered to >80% of its original length within approximately 1 week. It was also observed that the ventricles tended to shrink after amputation and maintained this state for at least 3 days. This shrinkage was transient, as the ventricles expanded to exceed their original size within the following week. Temporary shrinkage and expansion of the ventricles, which were also observed in transgenic or fluorescent dye-injected tadpoles with telencephalon amputation, may be significant in tissue homeostasis in response to massive brain injury and subsequent repair and regeneration. This established method will improve experimental analyses in developmental biology and medical science using Xenopus tadpoles.

Phenotype-genotype relationships in Xenopus sox9 crispants provide insights into campomelic dysplasia and vertebrate jaw evolution.

Since CRISPR-based genome editing technology works effectively in the diploid frog Xenopus tropicalis, a growing number of studies have successfully modeled human genetic diseases in this species. However, most of their targets were limited to non-syndromic diseases that exhibit abnormalities in a small fraction of tissues or organs in the body. This is likely because of the complexity of interpreting the phenotypic variations resulting from somatic mosaic mutations generated in the founder animals (crispants). In this study, we attempted to model the syndromic disease campomelic dysplasia (CD) by generating sox9 crispants in X. tropicalis. The resulting crispants failed to form neural crest cells at neurula stages and exhibited various combinations of jaw, gill, ear, heart, and gut defects at tadpole stages, recapitulating part of the syndromic phenotype of CD patients. Genotyping of the crispants with a variety of allelic series of mutations suggested that the heart and gut defects depend primarily on frame-shift mutations expected to be null, whereas the jaw, gill, and ear defects could be induced not only by such mutations but also by in-frame deletion mutations expected to delete part of the jawed vertebrate-specific domain from the encoded Sox9 protein. These results demonstrate that Xenopus crispants are useful for investigating the phenotype-genotype relationships behind syndromic diseases and examining the tissue-specific role of each functional domain within a single protein, providing novel insights into vertebrate jaw evolution.

Cryo-injury procedure-induced cardiac regeneration shows unique gene expression profiles in the newt Pleurodeles waltl.

BACKGROUND: Cardiac regeneration in the adult mouse is not substantial. Some vertebrates, such as newts and zebrafish, regenerate the heart throughout their lives. To understand how regenerative abilities differ among animal species, comparative research has been conducted in animals like mouse, zebrafish, and newt. For those purposes, cryo-injury is suitable as an experimental model for the pathological condition of human myocardial infarction. In fact, cryo-injury procedures are common in mouse and zebrafish. RESULTS: In the present study, we induced cryo-damage on the ventricle in Iberian ribbed newts using a liquid nitrogen-chilled probe. We observed that the injured area recovered within 8 weeks, with remodeling of scar tissue and proliferation of cardiomyocytes. We investigated the subsequent recovery of cryo-injured and amputated tissues by comparative analysis of the gene expression profiles following these two procedures. CONCLUSIONS: Notably, we established a cryo-injury procedure for the newt and confirmed that regeneration of the cryo-damaged myocardial tissue is achieved by changes in gene expression that are milder than those observed in the amputation model. Our results suggest that the cryo-injury method is suitable for comparing the process of cardiac regeneration in the newt with that in other animal models.

Optimization of CRISPR/Cas9-mediated gene disruption in Xenopus laevis using a phenotypic image analysis technique.

The CRISPR/Cas9 method has become popular for gene disruption experiments in Xenopus laevis. However, the experimental conditions that influence the efficiency of CRISPR/Cas9 remain unclear. To that end, we developed an image analysis technique for the semi-quantitative evaluation of the pigment phenotype resulting from the disruption of tyrosinase genes in X. laevis using a CRISPR/Cas9 approach, and then examined the effects of varying five experimental parameters (timing of the CRISPR reagent injection into developing embryos; amount of Cas9 mRNA in the injection reagent; total injection volume per embryo; number of injection sites per embryo; and the culture temperature of the injected embryos) on the gene disruption efficiency. The results of this systematic analysis suggest that the highest possible efficiency of target gene disruption can be achieved by injecting a total of 20 nL of the CRISPR reagent containing 1500 pg of Cas9 mRNA or 4 ng of Cas9 protein into two separate locations (10 nL each) of one-cell stage embryos cultured at 22°C. This study also highlights the importance of balancing the experimental parameters for increasing gene disruption efficiency and provides valuable insights into the optimal conditions for applying the CRISPR/Cas9 system to new experimental organisms.

Genome evolution in the allotetraploid frog Xenopus laevis.

To explore the origins and consequences of tetraploidy in the African clawed frog, we sequenced the Xenopus laevis genome and compared it to the related diploid X. tropicalis genome. We characterize the allotetraploid origin of X. laevis by partitioning its genome into two homoeologous subgenomes, marked by distinct families of 'fossil' transposable elements. On the basis of the activity of these elements and the age of hundreds of unitary pseudogenes, we estimate that the two diploid progenitor species diverged around 34 million years ago (Ma) and combined to form an allotetraploid around 17-18 Ma. More than 56% of all genes were retained in two homoeologous copies. Protein function, gene expression, and the amount of conserved flanking sequence all correlate with retention rates. The subgenomes have evolved asymmetrically, with one chromosome set more often preserving the ancestral state and the other experiencing more gene loss, deletion, rearrangement, and reduced gene expression.

Asymmetrically reduced expression of hand1 homeologs involving a single nucleotide substitution in a cis-regulatory element.

During vertebrate evolution, whole genome duplications resulted in a number of duplicated genes, some of which eventually changed their expression patterns and/or levels via alteration of cis-regulatory sequences. However, the initial process involved in such cis-regulatory changes remains unclear. Therefore, we investigated this process by analyzing the duplicated hand1 genes of Xenopus laevis (hand1.L and hand1.S), which were generated by allotetraploidization 17-18 million years ago, and compared these with their single ortholog in the ancestral-type diploid species X. tropicalis. A dN/dS analysis indicated that hand1.L and hand1.S are still under purifying selection, and thus, their products appear to retain ancestral functional properties. RNA-seq and in situ hybridization analyses revealed that hand1.L and hand1.S have similar expression patterns to each other and to X. tropicalis hand1, but the hand1.S expression level was much lower than the hand1.L expression level in the primordial heart. A comparative sequence analysis, luciferase reporter analysis, ChIP-PCR analysis, and transgenic reporter analysis showed that a single nucleotide substitution in the hand1.S promoter was responsible for the reduced expression in the heart. These findings demonstrated that a small change in the promoter sequence can trigger diversification of duplicated gene expression prior to diversification of their encoded protein functions in a young duplicated genome.

Co-accumulation of cis-regulatory and coding mutations during the pseudogenization of the Xenopus laevis homoeologs six6.L and six6.S.

Common models for the evolution of duplicated genes after genome duplication are subfunctionalization, neofunctionalization, and pseudogenization. Although the crucial roles of cis-regulatory mutations in subfunctionalization are well-documented, their involvement in pseudogenization and/or neofunctionalization remains unclear. We addressed this issue by investigating the evolution of duplicated homeobox genes, six6.L and six6.S, in the allotetraploid frog Xenopus laevis. Based on a comparative expression analysis, we observed similar eye-specific expression patterns for the two loci and their single ortholog in the ancestral-type diploid species Xenopus tropicalis. However, we detected lower levels of six6.S expression than six6.L expression. The six6.S enhancer sequence was more highly diverged from the orthologous enhancer of X. tropicalis than the six6.L enhancer, and showed weaker activity in a transgenic reporter assay. Based on a phylogenetic analysis of the protein sequences, we observed greater divergence between X. tropicalis Six6 and Six6.S than between X. tropicalis Six6 and Six6.L, and the observed mutations were reminiscent of a microphthalmia mutation in human SIX6. Misexpression experiments showed that six6.S has weaker eye-enlarging activity than six6.L, and targeted disruption of six6.L reduced the eye size more significantly than that of six6.S. These results suggest that enhancer attenuation stimulates the accumulation of hypomorphic coding mutations, or vice versa, in one duplicated gene copy and facilitates pseudogenization. We also underscore the value of the allotetraploid genome of X. laevis as a resource for studying latent pathogenic mutations.

Conservatism and variability of gene expression profiles among homeologous transcription factors in Xenopus laevis.

Xenopus laevis has an allotetraploid genome of 3.1Gb, in contrast to the diploid genome of a closely related species, Xenopus tropicalis. Here, we identified 412 genes (189 homeolog pairs, one homeologous gene cluster pair, and 28 singletons) encoding transcription factors (TFs) in the X. laevis genome by comparing them with their orthologs from X. tropicalis. Those genes include the homeobox gene family (Mix/Bix, Lhx, Nkx, Paired, POU, and Vent), Sox, Fox, Pax, Dmrt, Hes, GATA, T-box, and some clock genes. Most homeolog pairs for TFs are retained in two X. laevis subgenomes, named L and S, at higher than average rates (87.1% vs 60.2%). Among the 28 singletons, 82.1% were deleted from chromosomes of the S subgenome, a rate similar to the genome-wide average (82.1% vs 74.6%). Interestingly, nkx2-1, nkx2-8, and pax9, which reside consecutively in a postulated functional gene cluster, were deleted from the S chromosome, suggesting cluster-level gene regulation. Transcriptome correlation analysis demonstrated that TF homeolog pairs tend to have more conservative developmental expression profiles than most other types of genes. In some cases, however, either of the homeologs may show strongly different spatio-temporal expression patterns, suggesting neofunctionalization, subfunctionalization, or nonfunctionalization after allotetraploidization. Analyses of otx1 suggests that homeologs with much lower expression levels have undergone greater amino acid sequence diversification. Our comprehensive study implies that TF homeologs are highly conservative after allotetraploidization, possibly because the DNA sequences that they bind were also duplicated, but in some cases, they differed in expression levels or became singletons due to dosage-sensitive regulation of their target genes.

Epigenetic modification maintains intrinsic limb-cell identity in Xenopus limb bud regeneration.

Many amphibians can regenerate limbs, even in adulthood. If a limb is amputated, the stump generates a blastema that makes a complete, new limb in a process similar to developmental morphogenesis. The blastema is thought to inherit its limb-patterning properties from cells in the stump, and it retains the information despite changes in morphology, gene expression, and differentiation states required by limb regeneration. We hypothesized that these cellular properties are maintained as epigenetic memory through histone modifications. To test this hypothesis, we analyzed genome-wide histone modifications in Xenopus limb bud regeneration. The trimethylation of histone H3 at lysine 4 (H3K4me3) is closely related to an open chromatin structure that allows transcription factors access to genes, whereas the trimethylation of histone H3 at lysine 27 (H3K27me3) is related to a closed chromatin state that blocks the access of transcription factors. We compared these two modification profiles by high-throughput sequencing of samples prepared from the intact limb bud and the regenerative blastema by chromatin immunoprecipitation. For many developmental genes, histone modifications at the transcription start site were the same in the limb bud and the blastema, were stable during regeneration, and corresponded well to limb properties. These results support our hypothesis that histone modifications function as a heritable cellular memory to maintain limb cell properties, despite dynamic changes in gene expression during limb bud regeneration in Xenopus.

Transcriptional regulators in the Hippo signaling pathway control organ growth in Xenopus tadpole tail regeneration.

The size and shape of tissues are tightly controlled by synchronized processes among cells and tissues to produce an integrated organ. The Hippo signaling pathway controls both cell proliferation and apoptosis by dual signal-transduction states regulated through a repressive kinase cascade. Yap1 and Tead, transcriptional regulators that act downstream of the Hippo signaling kinase cascade, have essential roles in regulating cell proliferation. In amphibian limb or tail regeneration, the local tissue outgrowth terminates when the correct size is reached, suggesting that organ size is strictly controlled during epimorphic organ-level regeneration. We recently demonstrated that Yap1 is required for the regeneration of Xenopus tadpole limb buds (Hayashi et al., 2014, Dev. Biol. 388, 57-67), but the molecular link between the Hippo pathway and organ size control in vertebrate epimorphic regeneration is not fully understood. To examine the requirement of Hippo pathway transcriptional regulators in epimorphic regeneration, including organ size control, we inhibited these regulators during Xenopus tadpole tail regeneration by overexpressing a dominant-negative form of Yap (dnYap) or Tead4 (dnTead4) under a heat-shock promoter in transgenic animal lines. Each inhibition resulted in regeneration defects accompanied by reduced cell mitosis and increased apoptosis. Single-cell gene manipulation experiments indicated that Tead4 cell-autonomously regulates the survival of neural progenitor cells in the regenerating tail. In amphibians, amputation at the proximal level of the tail (deep amputation) results in faster regeneration than that at the distal level (shallow amputation), to restore the original-sized tail with similar timing. However, dnTead4 overexpression abolished the position-dependent differential growth rate of tail regeneration. These results suggest that the transcriptional regulators in the Hippo pathway, Tead4 and Yap1, are required for general vertebrate epimorphic regeneration as well as for organ size control in appendage regeneration. In regenerative medicine, these findings should contribute to the development of three-dimensional organs with the correct size for a patient's body.

Dynamic in vivo binding of transcription factors to cis-regulatory modules of cer and gsc in the stepwise formation of the Spemann-Mangold organizer.

How multiple developmental cues are integrated on cis-regulatory modules (CRMs) for cell fate decisions remains uncertain. The Spemann-Mangold organizer in Xenopus embryos expresses the transcription factors Lim1/Lhx1, Otx2, Mix1, Siamois (Sia) and VegT. Reporter analyses using sperm nuclear transplantation and DNA injection showed that cerberus (cer) and goosecoid (gsc) are activated by the aforementioned transcription factors through CRMs conserved between X. laevis and X. tropicalis. ChIP-qPCR analysis for the five transcription factors revealed that cer and gsc CRMs are initially bound by both Sia and VegT at the late blastula stage, and subsequently bound by all five factors at the gastrula stage. At the neurula stage, only binding of Lim1 and Otx2 to the gsc CRM, among others, persists, which corresponds to their co-expression in the prechordal plate. Based on these data, together with detailed expression pattern analysis, we propose a new model of stepwise formation of the organizer, in which (1) maternal VegT and Wnt-induced Sia first bind to CRMs at the blastula stage; then (2) Nodal-inducible Lim1, Otx2, Mix1 and zygotic VegT are bound to CRMs in the dorsal endodermal and mesodermal regions where all these genes are co-expressed; and (3) these two regions are combined at the gastrula stage to form the organizer. Thus, the in vivo dynamics of multiple transcription factors highlight their roles in the initiation and maintenance of gene expression, and also reveal the stepwise integration of maternal, Nodal and Wnt signaling on CRMs of organizer genes to generate the organizer.

Six1 is a key regulator of the developmental and evolutionary architecture of sensory neurons in craniates.

BACKGROUND: Various senses and sensory nerve architectures of animals have evolved during adaptation to exploit diverse environments. In craniates, the trunk sensory system has evolved from simple mechanosensory neurons inside the spinal cord (intramedullary), called Rohon-Beard (RB) cells, to multimodal sensory neurons of dorsal root ganglia (DRG) outside the spinal cord (extramedullary). The fish and amphibian trunk sensory systems switch from RB cells to DRG during development, while amniotes rely exclusively on the DRG system. The mechanisms underlying the ontogenic switching and its link to phylogenetic transition remain unknown. RESULTS: In Xenopus, Six1 overexpression promoted precocious apoptosis of RB cells and emergence of extramedullary sensory neurons, whereas Six1 knockdown delayed the reduction in RB cell number. Genetic ablation of Six1 and Six4 in mice led to the appearance of intramedullary sensory neuron-like cells as a result of medial migration of neural crest cells into the spinal cord and production of immature DRG neurons and fused DRG. Restoration of SIX1 expression in the neural crest-linage partially rescued the phenotype, indicating the cell autonomous requirements of SIX1 for normal extramedullary sensory neurogenesis. Mouse Six1 enhancer that mediates the expression in DRG neurons activated transcription in Xenopus RB cells earlier than endogenous six1 expression, suggesting earlier onset of mouse SIX1 expression than Xenopus during sensory development. CONCLUSIONS: The results indicated the critical role of Six1 in transition of RB cells to DRG neurons during Xenopus development and establishment of exclusive DRG system of mice. The study provided evidence that early appearance of SIX1 expression, which correlated with mouse Six1 enhancer, is essential for the formation of DRG-dominant system in mice, suggesting that heterochronic changes in Six1 enhancer sequence play an important role in alteration of trunk sensory architecture and contribute to the evolution of the trunk sensory system.