A chromatin code for limb segment identity in axolotl limb regeneration.

The salamander limb correctly regenerates missing limb segments because connective tissue cells have segment-specific identities, termed "positional information". How positional information is molecularly encoded at the chromatin level has been unknown. Here, we performed genome-wide chromatin profiling in mature and regenerating axolotl limb connective tissue cells. We find segment-specific levels of histone H3K27me3 as the major positional mark, especially at limb homeoprotein gene loci but not their upstream regulators, constituting an intrinsic segment information code. During regeneration, regeneration-specific regulatory elements became active prior to the re-appearance of developmental regulatory elements. In the hand, the permissive chromatin state of the homeoprotein gene HoxA13 engages with the regeneration program bypassing the upper limb program. Comparison of regeneration regulatory elements with those found in other regenerative animals identified a core shared set of transcription factors, supporting an ancient, conserved regeneration program.

Cellular senescence promotes progenitor cell expansion during axolotl limb regeneration.

Axolotl limb regeneration is accompanied by the transient induction of cellular senescence within the blastema, the structure that nucleates regeneration. The precise role of this blastemal senescent cell (bSC) population, however, remains unknown. Here, through a combination of gain- and loss-of-function assays, we elucidate the functions and molecular features of cellular senescence in vivo. We demonstrate that cellular senescence plays a positive role during axolotl regeneration by creating a pro-proliferative niche that supports progenitor cell expansion and blastema outgrowth. Senescent cells impact their microenvironment via Wnt pathway modulation. Further, we identify a link between Wnt signaling and senescence induction and propose that bSC-derived Wnt signals facilitate the proliferation of neighboring cells in part by preventing their induction into senescence. This work defines the roles of cellular senescence in the regeneration of complex structures.

Conserved enhancers control notochord expression of vertebrate Brachyury.

The cell type-specific expression of key transcription factors is central to development and disease. Brachyury/T/TBXT is a major transcription factor for gastrulation, tailbud patterning, and notochord formation; however, how its expression is controlled in the mammalian notochord has remained elusive. Here, we identify the complement of notochord-specific enhancers in the mammalian Brachyury/T/TBXT gene. Using transgenic assays in zebrafish, axolotl, and mouse, we discover three conserved Brachyury-controlling notochord enhancers, T3, C, and I, in human, mouse, and marsupial genomes. Acting as Brachyury-responsive, auto-regulatory shadow enhancers, in cis deletion of all three enhancers in mouse abolishes Brachyury/T/Tbxt expression selectively in the notochord, causing specific trunk and neural tube defects without gastrulation or tailbud defects. The three Brachyury-driving notochord enhancers are conserved beyond mammals in the brachyury/tbxtb loci of fishes, dating their origin to the last common ancestor of jawed vertebrates. Our data define the vertebrate enhancers for Brachyury/T/TBXTB notochord expression through an auto-regulatory mechanism that conveys robustness and adaptability as ancient basis for axis development.

Conserved enhancer logic controls the notochord expression of vertebrate Brachyury.

The cell type-specific expression of key transcription factors is central to development. Brachyury/T/TBXT is a major transcription factor for gastrulation, tailbud patterning, and notochord formation; however, how its expression is controlled in the mammalian notochord has remained elusive. Here, we identify the complement of notochord-specific enhancers in the mammalian Brachyury/T/TBXT gene. Using transgenic assays in zebrafish, axolotl, and mouse, we discover three Brachyury-controlling notochord enhancers T3, C, and I in human, mouse, and marsupial genomes. Acting as Brachyury-responsive, auto-regulatory shadow enhancers, deletion of all three enhancers in mouse abolishes Brachyury/T expression selectively in the notochord, causing specific trunk and neural tube defects without gastrulation or tailbud defects. Sequence and functional conservation of Brachyury-driving notochord enhancers with the brachyury/tbxtb loci from diverse lineages of fishes dates their origin to the last common ancestor of jawed vertebrates. Our data define the enhancers for Brachyury/T/TBXTB notochord expression as ancient mechanism in axis development.

Tig1 regulates proximo-distal identity during salamander limb regeneration.

Salamander limb regeneration is an accurate process which gives rise exclusively to the missing structures, irrespective of the amputation level. This suggests that cells in the stump have an awareness of their spatial location, a property termed positional identity. Little is known about how positional identity is encoded, in salamanders or other biological systems. Through single-cell RNAseq analysis, we identified Tig1/Rarres1 as a potential determinant of proximal identity. Tig1 encodes a conserved cell surface molecule, is regulated by retinoic acid and exhibits a graded expression along the proximo-distal axis of the limb. Its overexpression leads to regeneration defects in the distal elements and elicits proximal displacement of blastema cells, while its neutralisation blocks proximo-distal cell surface interactions. Critically, Tig1 reprogrammes distal cells to a proximal identity, upregulating Prod1 and inhibiting Hoxa13 and distal transcriptional networks. Thus, Tig1 is a central cell surface determinant of proximal identity in the salamander limb.

Postembryonic development and aging of the appendicular skeleton in Ambystoma mexicanum.

BACKGROUND: The axolotl is a key model to study appendicular regeneration. The limb complexity resembles that of humans in structure and tissue components; however, axolotl limbs develop postembryonically. In this work, we evaluated the postembryonic development of the appendicular skeleton and its changes with aging. RESULTS: The juvenile limb skeleton is formed mostly by Sox9/Col1a2 cartilage cells. Ossification of the appendicular skeleton starts when animals reach a length of 10 cm, and cartilage cells are replaced by a primary ossification center, consisting of cortical bone and an adipocyte-filled marrow cavity. Vascularization is associated with the ossification center and the marrow cavity formation. We identified the contribution of Col1a2-descendants to bone and adipocytes. Moreover, ossification progresses with age toward the epiphyses of long bones. Axolotls are neotenic salamanders, and still ossification remains responsive to l-thyroxine, increasing the rate of bone formation. CONCLUSIONS: In axolotls, bone maturation is a continuous process that extends throughout their life. Ossification of the appendicular bones is slow and continues until the complete element is ossified. The cellular components of the appendicular skeleton change accordingly during ossification, creating a heterogenous landscape in each element. The continuous maturation of the bone is accompanied by a continuous body growth.

A conserved regulatory program initiates lateral plate mesoderm emergence across chordates.

Cardiovascular lineages develop together with kidney, smooth muscle, and limb connective tissue progenitors from the lateral plate mesoderm (LPM). How the LPM initially emerges and how its downstream fates are molecularly interconnected remain unknown. Here, we isolate a pan-LPM enhancer in the zebrafish-specific draculin (drl) gene that provides specific LPM reporter activity from early gastrulation. In toto live imaging and lineage tracing of drl-based reporters captures the dynamic LPM emergence as lineage-restricted mesendoderm field. The drl pan-LPM enhancer responds to the transcription factors EomesoderminA, FoxH1, and MixL1 that combined with Smad activity drive LPM emergence. We uncover specific activity of zebrafish-derived drl reporters in LPM-corresponding territories of several chordates including chicken, axolotl, lamprey, Ciona, and amphioxus, revealing a universal upstream LPM program. Altogether, our work provides a mechanistic framework for LPM emergence as defined progenitor field, possibly representing an ancient mesodermal cell state that predates the primordial vertebrate embryo.

Publisher Correction: Application and optimization of CRISPR-Cas9-mediated genome engineering in axolotl (Ambystoma mexicanum).

In the version of this protocol originally published, the recipe for CAS9 buffer was incorrectly identified as a recipe for sodium acetate solution, and vice versa. These errors have been corrected in the PDF and HTML versions of the paper.

Application and optimization of CRISPR-Cas9-mediated genome engineering in axolotl (Ambystoma mexicanum).

Genomic manipulation is essential to the use of model organisms to understand development, regeneration and adult physiology. The axolotl (Ambystoma mexicanum), a type of salamander, exhibits an unparalleled regenerative capability in a spectrum of complex tissues and organs, and therefore serves as a powerful animal model for dissecting mechanisms of regeneration. We describe here an optimized stepwise protocol to create genetically modified axolotls using the CRISPR-Cas9 system. The protocol, which takes 7-8 weeks to complete, describes generation of targeted gene knockouts and knock-ins and includes site-specific integration of large targeting constructs. The direct use of purified CAS9-NLS (CAS9 containing a C-terminal nuclear localization signal) protein allows the prompt formation of guide RNA (gRNA)-CAS9-NLS ribonucleoprotein (RNP) complexes, which accelerates the creation of double-strand breaks (DSBs) at targeted genomic loci in single-cell-stage axolotl eggs. With this protocol, a substantial number of F0 individuals harboring a homozygous-type frameshift mutation can be obtained, allowing phenotype analysis in this generation. In the presence of targeting constructs, insertions of exogenous genes into targeted axolotl genomic loci can be achieved at efficiencies of up to 15% in a non-homologous end joining (NHEJ) manner. Our protocol bypasses the long generation time of axolotls and allows direct functional analysis in F0 genetically manipulated axolotls. This protocol can be potentially applied to other animal models, especially to organisms with a well-characterized transcriptome but lacking a well-characterized genome.

Single-cell analysis uncovers convergence of cell identities during axolotl limb regeneration.

Amputation of the axolotl forelimb results in the formation of a blastema, a transient tissue where progenitor cells accumulate prior to limb regeneration. However, the molecular understanding of blastema formation had previously been hampered by the inability to identify and isolate blastema precursor cells in the adult tissue. We have used a combination of Cre-loxP reporter lineage tracking and single-cell messenger RNA sequencing (scRNA-seq) to molecularly track mature connective tissue (CT) cell heterogeneity and its transition to a limb blastema state. We have uncovered a multiphasic molecular program where CT cell types found in the uninjured adult limb revert to a relatively homogenous progenitor state that recapitulates an embryonic limb bud-like phenotype including multipotency within the CT lineage. Together, our data illuminate molecular and cellular reprogramming during complex organ regeneration in a vertebrate.

Author Correction: The axolotl genome and the evolution of key tissue formation regulators.

In the originally published version of this Article, the sequenced axolotl strain (the homozygous white mutant) was denoted as 'D/D' rather than 'd/d' in Fig. 1a and the accompanying legend, the main text and the Methods section. The original Article has been corrected online.

The axolotl genome and the evolution of key tissue formation regulators.

Salamanders serve as important tetrapod models for developmental, regeneration and evolutionary studies. An extensive molecular toolkit makes the Mexican axolotl (Ambystoma mexicanum) a key representative salamander for molecular investigations. Here we report the sequencing and assembly of the 32-gigabase-pair axolotl genome using an approach that combined long-read sequencing, optical mapping and development of a new genome assembler (MARVEL). We observed a size expansion of introns and intergenic regions, largely attributable to multiplication of long terminal repeat retroelements. We provide evidence that intron size in developmental genes is under constraint and that species-restricted genes may contribute to limb regeneration. The axolotl genome assembly does not contain the essential developmental gene Pax3. However, mutation of the axolotl Pax3 paralogue Pax7 resulted in an axolotl phenotype that was similar to those seen in Pax3-/- and Pax7-/- mutant mice. The axolotl genome provides a rich biological resource for developmental and evolutionary studies.

Efficient gene knockin in axolotl and its use to test the role of satellite cells in limb regeneration.

Salamanders exhibit extensive regenerative capacities and serve as a unique model in regeneration research. However, due to the lack of targeted gene knockin approaches, it has been difficult to label and manipulate some of the cell populations that are crucial for understanding the mechanisms underlying regeneration. Here we have established highly efficient gene knockin approaches in the axolotl (Ambystoma mexicanum) based on the CRISPR/Cas9 technology. Using a homology-independent method, we successfully inserted both the Cherry reporter gene and a larger membrane-tagged Cherry-ERT2-Cre-ERT2 (∼5-kb) cassette into axolotl Sox2 and Pax7 genomic loci. Depending on the size of the DNA fragments for integration, 5-15% of the F0 transgenic axolotl are positive for the transgene. Using these techniques, we have labeled and traced the PAX7-positive satellite cells as a major source contributing to myogenesis during axolotl limb regeneration. Our work brings a key genetic tool to molecular and cellular studies of axolotl regeneration.

Tissue- and time-directed electroporation of CAS9 protein-gRNA complexes in vivo yields efficient multigene knockout for studying gene function in regeneration.

A rapid method for temporally and spatially controlled CRISPR-mediated gene knockout in vertebrates will be an important tool to screen for genes involved in complex biological phenomena like regeneration. Here we show that in vivo injection of CAS9 protein-guide RNA (gRNA) complexes into the spinal cord lumen of the axolotl and subsequent electroporation leads to comprehensive knockout of Sox2 gene expression in SOX2+ neural stem cells with corresponding functional phenotypes from the gene knockout. This is particularly surprising considering the known prevalence of RNase activity in cerebral spinal fluid, which apparently the CAS9 protein protects against. The penetrance/efficiency of gene knockout in the protein-based system is far higher than corresponding electroporation of plasmid-based CRISPR systems. We further show that simultaneous delivery of CAS9-gRNA complexes directed against Sox2 and GFP yields efficient knockout of both genes in GFP-reporter animals. Finally, we show that this method can also be applied to other tissues such as skin and limb mesenchyme. This efficient delivery method opens up the possibility for rapid in vivo genetic screens during axolotl regeneration and can in principle be applied to other vertebrate tissue systems.

Germline transgenic methods for tracking cells and testing gene function during regeneration in the axolotl.

The salamander is the only tetrapod that regenerates complex body structures throughout life. Deciphering the underlying molecular processes of regeneration is fundamental for regenerative medicine and developmental biology, but the model organism had limited tools for molecular analysis. We describe a comprehensive set of germline transgenic strains in the laboratory-bred salamander Ambystoma mexicanum (axolotl) that open up the cellular and molecular genetic dissection of regeneration. We demonstrate tissue-dependent control of gene expression in nerve, Schwann cells, oligodendrocytes, muscle, epidermis, and cartilage. Furthermore, we demonstrate the use of tamoxifen-induced Cre/loxP-mediated recombination to indelibly mark different cell types. Finally, we inducibly overexpress the cell-cycle inhibitor p16 (INK4a) , which negatively regulates spinal cord regeneration. These tissue-specific germline axolotl lines and tightly inducible Cre drivers and LoxP reporter lines render this classical regeneration model molecularly accessible.

Comparative transcriptional profiling of the axolotl limb identifies a tripartite regeneration-specific gene program.

Understanding how the limb blastema is established after the initial wound healing response is an important aspect of regeneration research. Here we performed parallel expression profile time courses of healing lateral wounds versus amputated limbs in axolotl. This comparison between wound healing and regeneration allowed us to identify amputation-specific genes. By clustering the expression profiles of these samples, we could detect three distinguishable phases of gene expression - early wound healing followed by a transition-phase leading to establishment of the limb development program, which correspond to the three phases of limb regeneration that had been defined by morphological criteria. By focusing on the transition-phase, we identified 93 strictly amputation-associated genes many of which are implicated in oxidative-stress response, chromatin modification, epithelial development or limb development. We further classified the genes based on whether they were or were not significantly expressed in the developing limb bud. The specific localization of 53 selected candidates within the blastema was investigated by in situ hybridization. In summary, we identified a set of genes that are expressed specifically during regeneration and are therefore, likely candidates for the regulation of blastema formation.

Connective tissue cells, but not muscle cells, are involved in establishing the proximo-distal outcome of limb regeneration in the axolotl.

During salamander limb regeneration, only the structures distal to the amputation plane are regenerated, a property known as the rule of distal transformation. Multiple cell types are involved in limb regeneration; therefore, determining which cell types participate in distal transformation is important for understanding how the proximo-distal outcome of regeneration is achieved. We show that connective tissue-derived blastema cells obey the rule of distal transformation. They also have nuclear MEIS, which can act as an upper arm identity regulator, only upon upper arm amputation. By contrast, myogenic cells do not obey the rule of distal transformation and display nuclear MEIS upon amputation at any proximo-distal level. These results indicate that connective tissue cells, but not myogenic cells, are involved in establishing the proximo-distal outcome of regeneration and are likely to guide muscle patterning. Moreover, we show that, similarly to limb development, muscle patterning in regeneration is influenced by ß-catenin signalling.

Regeneration and reprogramming.

Recent reprogramming studies indicate that mammalian, somatic cells have the potential to achieve pluripotent states and undergo cell type switching. Such cellular traits are observed under natural conditions in animals that regenerate complex organs. A number of invertebrates display the amazing trait of whole body regeneration. Underlying this trait is the maintenance of pluripotent cells in somatic tissue, and molecular studies indicate the use of common players associated with pluripotency and germ cell properties between these invertebrates and mammalian pluripotent cells. In regenerative vertebrates, heart regeneration, lens regeneration, and retinal regeneration provide good examples of dedifferentiation and transdifferentiation. The molecular factors associated with these phenomena are discussed.