Background: Microarray analysis and 454 cDNA sequencing were used to investigate a centuries-old problem in regenerative biology: the basis of nerve-dependent limb regeneration in salamanders. Innervated (NR) and denervated (DL) forelimbs of Mexican axolotls were amputated and transcripts were sampled after 0, 5, and 14 days of regeneration. Results: Considerable similarity was observed between NR and DL transcriptional programs at 5 and 14 days post amputation (dpa). Genes with extracellular functions that are critical to wound healing were upregulated while muscle-specific genes were downregulated. Thus, many processes that are regulated during early limb regeneration do not depend upon nerve-derived factors. The majority of the transcriptional differences between NR and DL limbs were correlated with blastema formation; cell numbers increased in NR limbs after 5 dpa and this yielded distinct transcriptional signatures of cell proliferation in NR limbs at 14 dpa. These transcriptional signatures were not observed in DL limbs. Instead, gene expression changes within DL limbs suggest more diverse and protracted wound-healing responses. 454 cDNA sequencing complemented the microarray analysis by providing deeper sampling of transcriptional programs and associated biological processes. Assembly of new 454 cDNA sequences with existing expressed sequence tag (EST) contigs from the Ambystoma EST database more than doubled (3935 to 9411) the number of non-redundant human-A. mexicanum orthologous sequences. Conclusion: Many new candidate gene sequences were discovered for the first time and these will greatly enable future studies of wound healing, epigenetics, genome stability, and nerve-dependent blastema formation and outgrowth using the axolotl model.

Background: Following amputation, urodele salamander limbs reprogram somatic cells to form a blastema that self-organizes into the missing limb parts to restore the structure and function of the limb. To help understand the molecular basis of blastema formation, we used quantitative label-free liquid chromatography-mass spectrometry/mass spectrometry (LC-MS/MS)-based methods to analyze changes in the proteome that occurred 1, 4 and 7 days post amputation (dpa) through the mid-tibia/fibula of axolotl hind limbs. Results: We identified 309 unique proteins with significant fold change relative to controls (0 dpa), representing 10 biological process categories: (1) signaling, (2) Ca²⁺ binding and translocation, (3) transcription, (4) translation, (5) cytoskeleton, (6) extracellular matrix (ECM), (7) metabolism, (8) cell protection, (9) degradation, and (10) cell cycle. In all, 43 proteins exhibited exceptionally high fold changes. Of these, the ecotropic viral integrative factor 5 (EVI5), a cell cycle-related oncoprotein that prevents cells from entering the mitotic phase of the cell cycle prematurely, was of special interest because its fold change was exceptionally high throughout blastema formation. Conclusion: Our data were consistent with previous studies indicating the importance of inositol triphosphate and Ca²⁺ signaling in initiating the ECM and cytoskeletal remodeling characteristic of histolysis and cell dedifferentiation. In addition, the data suggested that blastema formation requires several mechanisms to avoid apoptosis, including reduced metabolism, differential regulation of proapoptotic and antiapoptotic proteins, and initiation of an unfolded protein response (UPR). Since there is virtually no mitosis during blastema formation, we propose that high levels of EVI5 function to arrest dedifferentiated cells somewhere in the G₁/S/G₂ phases of the cell cycle until they have accumulated under the wound epidermis and enter mitosis in response to neural and epidermal factors. Our findings indicate the general value of quantitative proteomic analysis in understanding the regeneration of complex structures.

The axolotl (Ambystoma mexicanum) provides critical models for studying regeneration, evolution, and development. However, its large genome (∼32 Gb) presents a formidable barrier to genetic analyses. Recent efforts have yielded genome assemblies consisting of thousands of unordered scaffolds that resolve gene structures, but do not yet permit large-scale analyses of genome structure and function. We adapted an established mapping approach to leverage dense SNP typing information and for the first time assemble the axolotl genome into 14 chromosomes. Moreover, we used fluorescence in situ hybridization to verify the structure of these 14 scaffolds and assign each to its corresponding physical chromosome. This new assembly covers 27.3 Gb and encompasses 94% of annotated gene models on chromosomal scaffolds. We show the assembly's utility by resolving genome-wide orthologies between the axolotl and other vertebrates, identifying the footprints of historical introgression events that occurred during the development of axolotl genetic stocks, and precisely mapping several phenotypes including a large deletion underlying the cardiac mutant. This chromosome-scale assembly will greatly facilitate studies of the axolotl in biological research.

Vertebrates harbor recognizably orthologous gene complements but vary 100-fold in genome size. How chromosomal organization scales with genome expansion is unclear, and how acute changes in gene regulation, as during axolotl limb regeneration, occur in the context of a vast genome has remained a riddle. Here, we describe the chromosome-scale assembly of the giant, 32 Gb axolotl genome. Hi-C contact data revealed the scaling properties of interphase and mitotic chromosome organization. Analysis of the assembly yielded understanding of the evolution of large, syntenic multigene clusters, including the Major Histocompatibility Complex (MHC) and the functional regulatory landscape of the Fibroblast Growth Factor 8 (Axfgf8) region. The axolotl serves as a primary model for studying successful regeneration. Significance The axolotl is an important model organism because it is a tetrapod with a similar body plan to humans. Unlike humans, the axolotl regenerates limbs and other complex tissues. Therefore, the axolotl contributes to understanding evolution, development, and regeneration. With sophisticated tools for gene modification and tissue labeling, a fully assembled genome sequence was a sorely missing resource. Assembly was difficult because the genome size is 10× that of humans. Here, we use a cross-linking strategy called Hi-C to link together fragmented genome sequences to chromosome scale. We show that gene regulation occurs over very large genomic distances and that mitotic chromosomes are packaged efficiently.

In contrast to mammals, salamanders have a remarkable ability to regenerate their spinal cord and recover full movement and function after tail amputation. To identify genes that may be associated with this greater regenerative ability, we designed an oligonucleotide microarray and profiled early gene expression during natural spinal cord regeneration in Ambystoma mexicanum. We sampled tissue at five early time points after tail amputation and identified genes that registered significant changes in mRNA abundance during the first 7 days of regeneration. A list of 1036 statistically significant genes was identified. Additional statistical and fold change criteria were applied to identify a smaller list of 360 genes that were used to describe predominant expression patterns and gene functions. Our results show that a diverse injury response is activated in concert with extracellular matrix remodeling mechanisms during the early acute phase of natural spinal cord regeneration. We also report gene expression similarities and differences between our study and studies that have profiled gene expression after spinal cord injury in rat. Our study illustrates the utility of a salamander model for identifying genes and gene functions that may enhance regenerative ability in mammals.

The molecular genetic toolkit of the Mexican axolotl, a classic model organism, has matured to the point where it is now possible to identify genes for mutant phenotypes. We used a positional cloning–candidate gene approach to identify molecular bases for two historic axolotl pigment phenotypes: white and albino. White (d/d) mutants have defects in pigment cell morphogenesis and differentiation, whereas albino (a/a) mutants lack melanin. We identified in white mutants a transcriptional defect in endothelin 3 (edn3), encoding a peptide factor that promotes pigment cell migration and differentiation in other vertebrates. Transgenic restoration of Edn3 expression rescued the homozygous white mutant phenotype. We mapped the albino locus to tyrosinase (tyr) and identified polymorphisms shared between the albino allele (tyr ^a ) and tyr alleles in a Minnesota population of tiger salamanders from which the albino trait was introgressed. tyr ^a has a 142 bp deletion and similar engineered alleles recapitulated the albino phenotype. Finally, we show that historical introgression of tyr ^a significantly altered genomic composition of the laboratory axolotl, yielding a distinct, hybrid strain of ambystomatid salamander. Our results demonstrate the feasibility of identifying genes for traits in the laboratory Mexican axolotl.