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.

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.

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.

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.

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.

Cells keep a memory of their tissue origin during axolotl limb regeneration.

During limb regeneration adult tissue is converted into a zone of undifferentiated progenitors called the blastema that reforms the diverse tissues of the limb. Previous experiments have led to wide acceptance that limb tissues dedifferentiate to form pluripotent cells. Here we have reexamined this question using an integrated GFP transgene to track the major limb tissues during limb regeneration in the salamander Ambystoma mexicanum (the axolotl). Surprisingly, we find that each tissue produces progenitor cells with restricted potential. Therefore, the blastema is a heterogeneous collection of restricted progenitor cells. On the basis of these findings, we further demonstrate that positional identity is a cell-type-specific property of blastema cells, in which cartilage-derived blastema cells harbour positional identity but Schwann-derived cells do not. Our results show that the complex phenomenon of limb regeneration can be achieved without complete dedifferentiation to a pluripotent state, a conclusion with important implications for regenerative medicine.

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.

Efficient regeneration by activation of neurogenesis in homeostatically quiescent regions of the adult vertebrate brain.

In contrast to mammals, salamanders and teleost fishes can efficiently repair the adult brain. It has been hypothesised that constitutively active neurogenic niches are a prerequisite for extensive neuronal regeneration capacity. Here, we show that the highly regenerative salamander, the red spotted newt, displays an unexpectedly similar distribution of active germinal niches with mammals under normal physiological conditions. Proliferation zones in the adult newt brain are restricted to the forebrain, whereas all other regions are essentially quiescent. However, ablation of midbrain dopamine neurons in newts induced ependymoglia cells in the normally quiescent midbrain to proliferate and to undertake full dopamine neuron regeneration. Using oligonucleotide microarrays, we have catalogued a set of differentially expressed genes in these activated ependymoglia cells. This strategy identified hedgehog signalling as a key component of adult dopamine neuron regeneration. These data show that brain regeneration can occur by activation of neurogenesis in quiescent brain regions.

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.

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.

Gene expression profile of the regeneration epithelium during axolotl limb regeneration.

Urodele amphibians are unique among adult vertebrates in their ability to regenerate missing limbs. The process of limb regeneration requires several key tissues including a regeneration-competent wound epidermis called the regeneration epithelium (RE). We used microarray analysis to profile gene expression of the RE in the axolotl, a Mexican salamander. A list of 125 genes and expressed sequence tags (ESTs) showed a ≥1.5-fold expression in the RE than in a wound epidermis covering a lateral cuff wound. A subset of the RE ESTs and genes were further characterized for expression level changes over the time-course of regeneration. This study provides the first large scale identification of specific gene expression in the RE.

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.

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.

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.

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.

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.

Identification and characterisation of the posteriorly-expressed Xenopus neurotrophin receptor homolog genes fullback and fullback-like.

We have identified fullback and fullback-like, two Xenopus laevis neurotrophin receptor homolog (NRH1) genes. The sequences of Fullback and Fullback-like are very similar to that of the neurotrophin receptor p75NTR, in both their extracellular and their intracellular domains. As their names imply, fullback and fullback-like are expressed in essentially identical patterns in the posterior of the embryo from the early gastrula stage onward. At tailbud and tadpole stages transcripts are also present in dorsal somites and the head, in addition to the growing tailbud. This expression pattern differs from that of p75NTR, suggesting that fullback and fullback-like have different functions from p75NTR during early development.