NANOG is required to establish the competence for germ-layer differentiation in the basal tetrapod axolotl.

Pluripotency defines the unlimited potential of individual cells of vertebrate embryos, from which all adult somatic cells and germ cells are derived. Understanding how the programming of pluripotency evolved has been obscured in part by a lack of data from lower vertebrates; in model systems such as frogs and zebrafish, the function of the pluripotency genes NANOG and POU5F1 have diverged. Here, we investigated how the axolotl ortholog of NANOG programs pluripotency during development. Axolotl NANOG is absolutely required for gastrulation and germ-layer commitment. We show that in axolotl primitive ectoderm (animal caps; ACs) NANOG and NODAL activity, as well as the epigenetic modifying enzyme DPY30, are required for the mass deposition of H3K4me3 in pluripotent chromatin. We also demonstrate that all 3 protein activities are required for ACs to establish the competency to differentiate toward mesoderm. Our results suggest the ancient function of NANOG may be establishing the competence for lineage differentiation in early cells. These observations provide insights into embryonic development in the tetrapod ancestor from which terrestrial vertebrates evolved.

Downregulation of Yap1 during limb regeneration results in defective bone formation in axolotl.

The Hippo pathway plays an imperative role in cellular processes such as differentiation, regeneration, cell migration, organ growth, apoptosis, and cell cycle. Transcription coregulator component of Hippo pathway, YAP1, promotes transcription of genes involved in cell proliferation, migration, differentiation, and suppressing apoptosis. However, its role in epimorphic regeneration has not been fully explored. The axolotl is a well-established model organism for developmental biology and regeneration studies. By exploiting its remarkable regenerative capacity, we investigated the role of Yap1 in the early blastema stage of limb regeneration. Depleting Yap1 using gene-specific morpholinos attenuated the competence of axolotl limb regeneration evident in bone formation defects. To explore the affected downstream pathways from Yap1 down-regulation, the gene expression profile was examined by employing LC-MS/MS technology. Based on the generated data, we provided a new layer of evidence on the putative roles of increased protease inhibition and immune system activities and altered ECM composition in diminished bone formation capacity during axolotl limb regeneration upon Yap1 deficiency. We believe that new insights into the roles of the Hippo pathway in complex structure regeneration were granted in this study.

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.

The giant axolotl genome uncovers the evolution, scaling, and transcriptional control of complex gene loci.

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.

Geography is more important than life history in the recent diversification of the tiger salamander complex.

The North American tiger salamander species complex, including its best-known species, the Mexican axolotl, has long been a source of biological fascination. The complex exhibits a wide range of variation in developmental life history strategies, including populations and individuals that undergo metamorphosis; those able to forego metamorphosis and retain a larval, aquatic lifestyle (i.e., paedomorphosis); and those that do both. The evolution of a paedomorphic life history state is thought to lead to increased population genetic differentiation and ultimately reproductive isolation and speciation, but the degree to which it has shaped population- and species-level divergence is poorly understood. Using a large multilocus dataset from hundreds of samples across North America, we identified genetic clusters across the geographic range of the tiger salamander complex. These clusters often contain a mixture of paedomorphic and metamorphic taxa, indicating that geographic isolation has played a larger role in lineage divergence than paedomorphosis in this system. This conclusion is bolstered by geography-informed analyses indicating no effect of life history strategy on population genetic differentiation and by model-based population genetic analyses demonstrating gene flow between adjacent metamorphic and paedomorphic populations. This fine-scale genetic perspective on life history variation establishes a framework for understanding how plasticity, local adaptation, and gene flow contribute to lineage divergence. Many members of the tiger salamander complex are endangered, and the Mexican axolotl is an important model system in regenerative and biomedical research. Our results chart a course for more informed use of these taxa in experimental, ecological, and conservation research.

SpecHap: a diploid phasing algorithm based on spectral graph theory.

Haplotype phasing plays an important role in understanding the genetic data of diploid eukaryotic organisms. Different sequencing technologies (such as next-generation sequencing or third-generation sequencing) produce various genetic data that require haplotype assembly. Although multiple diploid haplotype phasing algorithms exist, only a few will work equally well across all sequencing technologies. In this work, we propose SpecHap, a novel haplotype assembly tool that leverages spectral graph theory. On both in silico and whole-genome sequencing datasets, SpecHap consumed less memory and required less CPU time, yet achieved comparable accuracy with state-of-art methods across all the test instances, which comprises sequencing data from next-generation sequencing, linked-reads, high-throughput chromosome conformation capture, PacBio single-molecule real-time, and Oxford Nanopore long-reads. Furthermore, SpecHap successfully phased an individual Ambystoma mexicanum, a species with gigantic diploid genomes, within 6 CPU hours and 945MB peak memory usage, while other tools failed to yield results either due to memory overflow (40GB) or time limit exceeded (5 days). Our results demonstrated that SpecHap is scalable, efficient, and accurate for diploid phasing across many sequencing platforms.

Gene and transgenics nomenclature for the laboratory axolotl-Ambystoma mexicanum.

The laboratory axolotl (Ambystoma mexicanum) is widely used in biological research. Recent advancements in genetic and molecular toolkits are greatly accelerating the work using axolotl, especially in the area of tissue regeneration. At this juncture, there is a critical need to establish gene and transgenic nomenclature to ensure uniformity in axolotl research. Here, we propose guidelines for genetic nomenclature when working with the axolotl.

A chromosome-scale assembly of the axolotl genome.

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.

Characterization of axolotl lampbrush chromosomes by fluorescence in situ hybridization and immunostaining.

The lampbrush chromosomes (LBCs) in oocytes of the Mexican axolotl (Ambystoma mexicanum) were identified some time ago by their relative lengths and predicted centromeres, but they have never been associated completely with the mitotic karyotype, linkage maps or genome assembly. We identified 9 of the axolotl LBCs using RNAseq to identify actively transcribed genes and 13 BAC (bacterial artificial clone) probes containing pieces of active genes. Using read coverage analysis to find candidate centromere sequences, we developed a centromere probe that localizes to all 14 centromeres. Measurements of relative LBC arm lengths and polymerase III localization patterns enabled us to identify all LBCs. This study presents a relatively simple and reliable way to identify each axolotl LBC cytologically and to anchor chromosome-length sequences (from the axolotl genome assembly) to the physical LBCs by immunostaining and fluorescence in situ hybridization. Our data will facilitate a more detailed transcription analysis of individual LBC loops.

Small circRNAs with self-cleaving ribozymes are highly expressed in diverse metazoan transcriptomes.

Ribozymes are catalytic RNAs present in modern genomes but regarded as remnants of a prebiotic RNA world. The paradigmatic hammerhead ribozyme (HHR) is a small self-cleaving motif widespread from bacterial to human genomes. Here, we report that most of the classical type I HHRs frequently found in the genomes of animals are contained within a novel family of non-autonomous non-LTR retrotransposons of the retrozyme class. These retroelements are expressed as abundant linear and circular RNAs of ∼170-400 nt in different animal tissues. Bioinformatic and in vitro analyses indicate an efficient self-cleavage of the HHRs harboured in most invertebrate retrozymes, whereas HHRs in retrozymes of vertebrates, such as the axolotl and other amphibians, require to act as dimeric motifs to reach higher self-cleavage rates. Ligation assays of retrozyme RNAs with a protein ligase versus HHR self-ligation indicate that, most likely, tRNA ligases and not the ribozymes are involved in the step of RNA circularization. Altogether, these results confirm the existence of a new and conserved pathway in animals and, likely, eukaryotes in general, for the efficient biosynthesis of RNA circles through small ribozymes, which opens the door for the development of new tools in the emerging field of study of circRNAs.

Deep evolutionary origin of limb and fin regeneration.

Salamanders and lungfishes are the only sarcopterygians (lobe-finned vertebrates) capable of paired appendage regeneration, regardless of the amputation level. Among actinopterygians (ray-finned fishes), regeneration after amputation at the fin endoskeleton has only been demonstrated in polypterid fishes (Cladistia). Whether this ability evolved independently in sarcopterygians and actinopterygians or has a common origin remains unknown. Here we combine fin regeneration assays and comparative RNA-sequencing (RNA-seq) analysis of Polypterus and axolotl blastemas to provide support for a common origin of paired appendage regeneration in Osteichthyes (bony vertebrates). We show that, in addition to polypterids, regeneration after fin endoskeleton amputation occurs in extant representatives of 2 other nonteleost actinopterygians: the American paddlefish (Chondrostei) and the spotted gar (Holostei). Furthermore, we assessed regeneration in 4 teleost species and show that, with the exception of the blue gourami (Anabantidae), 3 species were capable of regenerating fins after endoskeleton amputation: the white convict and the oscar (Cichlidae), and the goldfish (Cyprinidae). Our comparative RNA-seq analysis of regenerating blastemas of axolotl and Polypterus reveals the activation of common genetic pathways and expression profiles, consistent with a shared genetic program of appendage regeneration. Comparison of RNA-seq data from early Polypterus blastema to single-cell RNA-seq data from axolotl limb bud and limb regeneration stages shows that Polypterus and axolotl share a regeneration-specific genetic program. Collectively, our findings support a deep evolutionary origin of paired appendage regeneration in Osteichthyes and provide an evolutionary framework for studies on the genetic basis of appendage regeneration.

Large-scale variation in single nucleotide polymorphism density within the laboratory axolotl (Ambystoma mexicanum).

BACKGROUND: Recent efforts to assemble and analyze the Ambystoma mexicanum genome have dramatically improved the potential to develop molecular tools and pursue genome-wide analyses of genetic variation. RESULTS: To better resolve the distribution and origins of genetic variation with A mexicanum, we compared DNA sequence data for two laboratory A mexicanum and one A tigrinum to identify 702 million high confidence polymorphisms distributed across the 32 Gb genome. While the wild-caught A tigrinum was generally more polymorphic in a genome-wide sense, several multi-megabase regions were identified from A mexicanum genomes that were actually more polymorphic than A tigrinum. Analysis of polymorphism and repeat content reveals that these regions likely originated from the intentional hybridization of A mexicanum and A tigrinum that was used to introduce the albino mutation into laboratory stocks. CONCLUSIONS: Our findings show that axolotl genomes are variable with respect to introgressed DNA from a highly polymorphic species. It seems likely that other divergent regions will be discovered with additional sequencing of A mexicanum. This has practical implications for designing molecular probes and suggests a need to study A mexicanum phenotypic variation and genome evolution across the tiger salamander clade.

The rax homeobox gene is mutated in the eyeless axolotl, Ambystoma mexicanum.

BACKGROUND: Vertebrate eye formation requires coordinated inductive interactions between different embryonic tissue layers, first described in amphibians. A network of transcription factors and signaling molecules controls these steps, with mutations causing severe ocular, neuronal, and craniofacial defects. In eyeless mutant axolotls, eye morphogenesis arrests at the optic vesicle stage, before lens induction, and development of ventral forebrain structures is disrupted. RESULTS: We identified a 5-bp deletion in the rax (retina and anterior neural fold homeobox) gene, which was tightly linked to the recessive eyeless (e) axolotl locus in an F2 cross. This frameshift mutation, in exon 2, truncates RAX protein within the homeodomain (P154fs35X). Quantitative RNA analysis shows that mutant and wild-type rax transcripts are equally abundant in E/e embryos. Translation appears to initiate from dual start codons, via leaky ribosome scanning, a conserved feature among gnathostome RAX proteins. Previous data show rax is expressed in the optic vesicle and diencephalon, deeply conserved among metazoans, and required for eye formation in other species. CONCLUSION: The eyeless axolotl mutation is a null allele in the rax homeobox gene, with primary defects in neural ectoderm, including the retinal and hypothalamic primordia.

Exploring the role of microRNAs in axolotl regeneration.

The axolotl, Ambystoma mexicanum, is used extensively for research in developmental biology, particularly for its ability to regenerate and restore lost organs, including in the nervous system, to full functionality. Regeneration in mammals typically depends on the healing process and scar formation with limited replacement of lost tissue. Other organisms, such as spiny mice (Acomys cahirinus), salamanders, and zebrafish, are able to regenerate some damaged body components. Blastema is a tissue that is formed after tissue injury in such organisms and is composed of progenitor cells or dedifferentiated cells that differentiate into various cell types during regeneration. Thus, identifying the molecules responsible for initiation of blastema formation is an important aspect for understanding regeneration. Introns, a major source of noncoding RNAs (ncRNAs), have characteristic sizes in the axolotl, particularly in genes associated with development. These ncRNAs, particularly microRNAs (miRNAs), exhibit dynamic regulation during regeneration. These miRNAs play an essential role in timing and control of gene expression to order and organize processes necessary for blastema creation. Master keys or molecules that underlie the remarkable regenerative abilities of the axolotl remain to be fully explored and exploited. Further and ongoing research on regeneration promises new knowledge that may allow improved repair and renewal of human tissues.

Introducing www.axolotl-omics.org - an integrated -omics data portal for the axolotl research community.

Genomic resources are indispensable for biological investigations in model organisms. In recent years, a number of genomic resources including a full genome assembly, extensive transcriptomic data, as well as genome editing has been developed for the axolotl, a classical model organism for developmental, neurobiological and regeneration studies, making the axolotl a highly versatile system. Here we describe the Axolotl-omics website that allows rapid ortholog searches, and access to genome and transcriptomic resources.

Identification of immune and non-immune cells in regenerating axolotl limbs by single-cell sequencing.

Immune cells are known to be critical for successful limb regeneration in the axolotl (Ambystoma mexicanum), but many details regarding their identity, behavior, and function are yet to be resolved. We isolated peripheral leukocytes from the blood of adult axolotls and then created two samples for single-cell sequencing: 1) peripheral leukocytes (N = 7889) and 2) peripheral leukocytes with presumptive macrophages from the intraperitoneal cavity (N = 4998). Using k-means clustering, we identified 6 cell populations from each sample that presented gene expression patterns indicative of erythrocyte, thrombocyte, neutrophil, B-cell, T-cell, and myeloid cell populations. A seventh, presumptive macrophage cell population was identified uniquely from sample 2. We then isolated cells from amputated axolotl limbs at 1 and 6 days post-amputation (DPA) and performed single cell sequencing (N = 8272 and 9906 cells respectively) to identify immune and non-immune cell populations. Using k-means clustering, we identified 8 cell populations overall, with the majority of cells expressing erythrocyte-specific genes. Even though erythrocytes predominated, we used an unbiased approach to identify infiltrating neutrophil, macrophage, and lymphocyte populations at both time points. Additionally, populations expressing genes for epidermal cells, fibroblast-like cells, and endothelial cells were also identified. Consistent with results from previous experimental studies, neutrophils were more abundant at 1 DPA than 6 DPA, while macrophages and non-immune cells exhibited inverse abundance patterns. Of note, we identified a small population of fibroblast-like cells at 1 DPA that was represented by considerably more cells at 6 DPA. We hypothesize that these are early progenitor cells that give rise to the blastema. The enriched gene sets from our work will aid future single-cell investigations of immune cell diversity and function during axolotl limb regeneration.

Fgf- and Bmp-signaling regulate gill regeneration in Ambystoma mexicanum.

Gill regeneration has not been well studied compared to regeneration of other appendages, such as limb and tail regeneration. Here, we focused on axolotl gill regeneration and found that Fgf- and Bmp-signaling are involved in their gill regeneration mechanism. Axolotls have three pairs of gill rami, and each gill ramus has multiple gill filaments. The gills consist of mesenchyme rich in extracellular matrix and epidermis. The gill nerves are supplied from the trigeminal ganglia located in the head. Denervation resulted in no gill regeneration responses. Nerves and gills express Bmp and Fgf genes, and treating animals with Fgf- and Bmp-signaling inhibitors results in phenotypes similar to those seen in denervated gills. Inducing an accessory appendage is a standard assay in amphibian regeneration research. In our study, an accessory gill could be induced by lateral wounding, suggesting that thin axon fibers and mesenchymal Fgfs and Bmps contributed to the induction of the accessory structure. Such accessory gill induction was inhibited by the denervation. Exogenous Fgf2+Fgf8+Bmp7, which have been determined to function as a regeneration inducer in urodele amphibians, could compensate for the effects denervation has on accessory blastema formation. Our findings suggest that regeneration of appendages in axolotls is regulated by common Fgf- and Bmp-signaling cascades.

Lissamphibian limbs and the origins of tetrapod hox domains.

The expression and function of hox genes have played a key role in the debate on the evolution of limbs from fins. As an early branching tetrapod lineage, lissamphibians may provide information on the origin of the limb's hox domains and particularly how the plesiomorphic tetrapod pattern compares to the hox pattern present in fish fins. Here, we comparatively investigated the expression of hox genes in the developing limbs of axolotl and Xenopus laevis as well as in the fins of the direct developing cichlid Astatotilapia burtoni. In contrast to axolotl, which has only very low digital expression of hoxd11, Xenopus limbs recapitulate the reverse collinear hoxd expression pattern known from amniotes with clearly defined proximal and distal hoxd11 expression domains. For hoxa genes, we observe that in Xenopus limbs, as in axolotl, a clear distal domain of hoxa11 expression is present, although in the presence of a hoxa11 antisense transcript. Investigation of fins reveals the presence of hoxa11 antisense transcription in the developing fin rays in a domain similar to that of hoxa13 and overlapping with hoxa11 sense transcription. Our results indicate that full exclusion of hoxa11 from the autopod only became firmly established in amniotes. The distal antisense transcription of hoxa11, however, appears to predate the evolution of the limb, but likely originated without the concurrent implementation of the transcriptional suppression mechanism that causes mutually exclusive hoxa11 and hoxa13 domains in amniotes.

Variation under domestication in animal models: the case of the Mexican axolotl.

BACKGROUND: Species adaptation to laboratory conditions is a special case of domestication that has modified model organisms phenotypically and genetically. The characterisation of these changes is crucial to understand how this variation can affect the outcome of biological experiments. Yet despite the wide use of laboratory animals in biological research, knowledge of the genetic diversity within and between different strains and populations of some animal models is still scarce. This is particularly the case of the Mexican axolotl, which has been bred in captivity since 1864. RESULTS: Using gene expression data from nine different projects, nucleotide sequence variants were characterised, and distinctive genetic background of the experimental specimens was uncovered. This study provides a catalogue of thousands of nucleotide variants along predicted protein-coding genes, while identifying genome-wide differences between pigment phenotypes in laboratory populations. CONCLUSIONS: Awareness of the genetic variation could guide a better experimental design while helping to develop molecular tools for monitoring genetic diversity and studying gene functions in laboratory axolotls. Overall, this study highlights the cross-taxa utility that transcriptomic data might have to assess the genetic variation of the experimental specimens, which might help to shorten the journey towards reproducible research.

Advances in Decoding Axolotl Limb Regeneration.

Humans and other mammals are limited in their natural abilities to regenerate lost body parts. By contrast, many salamanders are highly regenerative and can spontaneously replace lost limbs even as adults. Because salamander limbs are anatomically similar to human limbs, knowing how they regenerate should provide important clues for regenerative medicine. Although interest in understanding the mechanics of this process has never wavered, until recently researchers have been vexed by seemingly impenetrable logistics of working with these creatures at a molecular level. Chief among the problems has been the very large size of salamander genomes, and not a single salamander genome has been fully sequenced to date. Recently the enormous gap in sequence information has been bridged by approaches that leverage mRNA as the starting point. Together with functional experimentation, these data are rapidly enabling researchers to finally uncover the molecular mechanisms underpinning the astonishing biological process of limb regeneration.