Comparative RNA-seq analysis in the unsequenced axolotl: the oncogene burst highlights early gene expression in the blastema

hoxc12/c13 as key regulators for rebooting the developmental program in Xenopus limb regeneration

During organ regeneration, after the initial responses to injury, gene expression patterns similar to those in normal development are reestablished during subsequent morphogenesis phases. This supports the idea that regeneration recapitulates development and predicts the existence of genes that reboot the developmental program after the initial responses. However, such rebooting mechanisms are largely unknown. Here, we explore core rebooting factors that operate during Xenopus limb regeneration. Transcriptomic analysis of larval limb blastema reveals that hoxc12/c13 show the highest regeneration specificity in expression. Knocking out each of them through genome editing inhibits cell proliferation and expression of a group of genes that are essential for development, resulting in autopod regeneration failure, while limb development and initial blastema formation are not affected. Furthermore, the induction of hoxc12/c13 expression partially restores froglet regenerative capacity which is normally very limited compared to larval regeneration. Thus, we demonstrate the existence of genes that have a profound impact alone on rebooting of the developmental program in a regeneration-specific manner.

Manipulating Myc for reparative regeneration

The Myc family of proto-oncogenes is a key node for the signal transduction of external pro-proliferative signals to the cellular processes required for development, tissue homoeostasis maintenance, and regeneration across evolution. The tight regulation of Myc synthesis and activity is essential for restricting its oncogenic potential. In this review, we highlight the central role that Myc plays in regeneration across the animal kingdom (from Cnidaria to echinoderms to Chordata) and how Myc could be employed to unlock the regenerative potential of non-regenerative tissues in humans for therapeutic purposes. Mastering the fine balance of harnessing the ability of Myc to promote transcription without triggering oncogenesis may open the door to many exciting opportunities for therapeutic development across a wide array of diseases.

The Emerging Concept of Transportome: State of the Art

The array of ion channels and transporters expressed in cell membranes, collectively referred to as the transportome, is a complex and multifunctional molecular machinery; in particular, at the plasma membrane level it finely tunes the exchange of biomolecules and ions, acting as a functionally adaptive interface that accounts for dynamic plasticity in the response to environmental fluctuations and stressors. The transportome is responsible for the definition of membrane potential and its variations, participates in the transduction of extracellular signals, and acts as a filter for most of the substances entering and leaving the cell, thus enabling the homeostasis of many cellular parameters. For all these reasons, physiologists have long been interested in the expression and functionality of ion channels and transporters, in both physiological and pathological settings and across the different domains of life. Today, thanks to the high-throughput technologies of the postgenomic era, the omics approach to the study of the transportome is becoming increasingly popular in different areas of biomedical research, allowing for a more comprehensive, integrated, and functional perspective of this complex cellular apparatus. This article represents a first effort for a systematic review of the scientific literature on this topic. Here we provide a brief overview of all those studies, both primary and meta-analyses, that looked at the transportome as a whole, regardless of the biological problem or the models they used. A subsequent section is devoted to the methodological aspect by reviewing the most important public databases annotating ion channels and transporters, along with the tools they provide to retrieve such information. Before conclusions, limitations and future perspectives are also discussed.

The salamander blastema within the broader context of metazoan regeneration

Throughout the animal kingdom regenerative ability varies greatly from species to species, and even tissue to tissue within the same organism. The sheer diversity of structures and mechanisms renders a thorough comparison of molecular processes truly daunting. Are "blastemas" found in organisms as distantly related as planarians and axolotls derived from the same ancestral process, or did they arise convergently and independently? Is a mouse digit tip blastema orthologous to a salamander limb blastema? In other fields, the thorough characterization of a reference model has greatly facilitated these comparisons. For example, the amphibian Spemann-Mangold organizer has served as an amazingly useful comparative template within the field of developmental biology, allowing researchers to draw analogies between distantly related species, and developmental processes which are superficially quite different. The salamander limb blastema may serve as the best starting point for a comparative analysis of regeneration, as it has been characterized by over 200 years of research and is supported by a growing arsenal of molecular tools. The anatomical and evolutionary closeness of the salamander and human limb also add value from a translational and therapeutic standpoint. Tracing the evolutionary origins of the salamander blastema, and its relatedness to other regenerative processes throughout the animal kingdom, will both enhance our basic biological understanding of regeneration and inform our selection of regenerative model systems.

Limb blastema formation: How much do we know at a genetic and epigenetic level?

Regeneration of missing body parts is an incredible ability which is present in a wide number of species. However, this regenerative capability varies among different organisms. Urodeles (salamanders) are able to completely regenerate limbs after amputation through the essential process of blastema formation. The blastema is a collection of relatively undifferentiated progenitor cells that proliferate and repattern to form the internal tissues of a regenerated limb. Understanding blastema formation in salamanders may enable comparative studies with other animals, including mammals, with more limited regenerative abilities and may inspire future therapeutic approaches in humans. This review focuses on the current state of knowledge about how limb blastemas form in salamanders, highlighting both the possible roles of epigenetic controls in this process as well as limitations to scientific understanding that present opportunities for research.

A Practical Guide for CRISPR-Cas9-Induced Mutations in Axolotls

Clustered regularly interspaced short palindromic repeats (CRISPR) is a powerful tool that enables editing of the axolotl genome. In this chapter, we will cover how to retrieve gene sequences, confirm annotation, design CRISPR targets, analyze indels, and screen for Crispant axolotls. This is a comprehensive guide on how to use CRISPR on your favorite gene and gain insights into its function.

Now that We Got There, What Next?

As seen in the protocols in this book, the opportunities to pursue work at the cellular and molecular work in salamanders have considerably broadened over the last years. The availability of genomic information and genome editing, and the possibility to image tissues live and other methods enhance the spectrum of biological questions accessible to all researchers. Here I provide a personal perspective on what I consider exciting future questions open for investigation.

The specialist in regeneration-the Axolotl-a suitable model to study bone healing?

While the axolotl's ability to completely regenerate amputated limbs is well known and studied, the mechanism of axolotl bone fracture healing remains poorly understood. One reason might be the lack of a standardized fracture fixation in axolotl. We present a surgical technique to stabilize the osteotomized axolotl femur with a fixator plate and compare it to a non-stabilized osteotomy and to limb amputation. The healing outcome was evaluated 3 weeks, 3, 6 and 9 months post-surgery by microcomputer tomography, histology and immunohistochemistry. Plate-fixated femurs regained bone integrity more efficiently in comparison to the non-fixated osteotomized bone, where larger callus formed, possibly to compensate for the bone fragment misalignment. The healing of a non-critical osteotomy in axolotl was incomplete after 9 months, while amputated limbs efficiently restored bone length and structure. In axolotl amputated limbs, plate-fixated and non-fixated fractures, we observed accumulation of PCNA+ proliferating cells at 3 weeks post-injury similar to mouse. Additionally, as in mouse, SOX9-expressing cells appeared in the early phase of fracture healing and amputated limb regeneration in axolotl, preceding cartilage formation. This implicates endochondral ossification to be the probable mechanism of bone healing in axolotls. Altogether, the surgery with a standardized fixation technique demonstrated here allows for controlled axolotl bone healing experiments, facilitating their comparison to mammals (mice).

Normal embryonic development and neonatal digit regeneration in mice overexpressing a stem cell factor, Sall4

Sall4 encodes a transcription factor and is known to participate in the pluripotency network of embryonic stem cells. Sall4 expression is known to be high in early stage post-implantation mouse embryos. During early post-gastrulation stages, Sall4 is highly expressed in the tail bud and distal limb buds, where progenitor cells are maintained in an undifferentiated status. The expression of Sall4 is rapidly downregulated during embryonic development. We previously demonstrated that Sall4 is required for limb and posterior axial skeleton development by conditional deletion of Sall4 in the T (Brachyury) lineage. To gain insight into Sall4 functions in embryonic development and postnatal digit regeneration, we genetically overexpressed Sall4 in the mesodermal lineage by the TCre transgene and a novel knockin allele of Rosa26-loxP-stop-loxP-Sall4. In significant contrast to severe defects by Sall4 loss of function reported in previous studies, overexpression of Sall4 resulted in normal morphology and pattern in embryos and neonates. The length of limb long bones showed subtle reduction in Sall4-overexpression mice. It is known that the digit tip of neonatal mice has level-specific regenerative ability after experimental amputation. We observed Sall4 expression in the digit tip by using a sensitive Sall4-LacZ knock-in reporter expression. Sall4 overexpression did not alter the regenerative ability of the terminal phalange that normally regenerates after amputation. Moreover, Sall4 overexpression did not confer regenerative ability to the second phalange that normally does not regenerate after amputation. These genetic experiments show that overexpression of Sall4 does not alter the development of the appendicular and axial skeleton, or neonatal digit regeneration. The results suggest that Sall4 acts as a permissive factor rather than playing an instructive role.

Wnt Signaling Coordinates the Expression of Limb Patterning Genes During Axolotl Forelimb Development and Regeneration

After amputation, axolotl salamanders can regenerate their limbs, but the degree to which limb regeneration recapitulates limb development remains unclear. One limitation in answering this question is our lack of knowledge about salamander limb development. Here, we address this question by studying expression patterns of genes important for limb patterning during axolotl salamander limb development and regeneration. We focus on the Wnt signaling pathway because it regulates multiple functions during tetrapod limb development, including limb bud initiation, outgrowth, patterning, and skeletal differentiation. We use fluorescence in situ hybridization to show the expression of Wnt ligands, Wnt receptors, and limb patterning genes in developing and regenerating limbs. Inhibition of Wnt ligand secretion permanently blocks limb bud outgrowth when treated early in limb development. Inhibiting Wnt signaling during limb outgrowth decreases the expression of critical signaling genes, including Fgf10, Fgf8, and Shh, leading to the reduced outgrowth of the limb. Patterns of gene expression are similar between developing and regenerating limbs. Inhibition of Wnt signaling during regeneration impacted patterning gene expression similarly. Overall, our findings suggest that limb development and regeneration utilize Wnt signaling similarly. It also provides new insights into the interaction of Wnt signaling with other signaling pathways during salamander limb development and regeneration.

Re-building limbs, one cell at a time

New techniques for visualizing and interrogating single cells hold the key to unlocking the underlying mechanisms of salamander limb regeneration.

Assembly-free rapid differential gene expression analysis in non-model organisms using DNA-protein alignment

Background

RNA-seq is being increasingly adopted for gene expression studies in a panoply of non-model organisms, with applications spanning the fields of agriculture, aquaculture, ecology, and environment. For organisms that lack a well-annotated reference genome or transcriptome, a conventional RNA-seq data analysis workflow requires constructing a de-novo transcriptome assembly and annotating it against a high-confidence protein database. The assembly serves as a reference for read mapping, and the annotation is necessary for functional analysis of genes found to be differentially expressed. However, assembly is computationally expensive. It is also prone to errors that impact expression analysis, especially since sequencing depth is typically much lower for expression studies than for transcript discovery.

Results

We propose a shortcut, in which we obtain counts for differential expression analysis by directly aligning RNA-seq reads to the high-confidence proteome that would have been otherwise used for annotation. By avoiding assembly, we drastically cut down computational costs – the running time on a typical dataset improves from the order of tens of hours to under half an hour, and the memory requirement is reduced from the order of tens of Gbytes to tens of Mbytes. We show through experiments on simulated and real data that our pipeline not only reduces computational costs, but has higher sensitivity and precision than a typical assembly-based pipeline. A Snakemake implementation of our workflow is available at: https://bitbucket.org/project_samar/samar.

Conclusions

The flip side of RNA-seq becoming accessible to even modestly resourced labs has been that the time, labor, and infrastructure cost of bioinformatics analysis has become a bottleneck. Assembly is one such resource-hungry process, and we show here that it can be avoided for quick and easy, yet more sensitive and precise, differential gene expression analysis in non-model organisms.

Supplementary Information

The online version contains supplementary material available at (10.1186/s12864-021-08278-7).

Evidence of requirement for homologous-mediated DNA repair during Ambystoma mexicanum limb regeneration

BACKGROUND

Limb regeneration in the axolotl is achieved by epimorphosis, thus depending on the blastema formation, a mass of progenitor cells capable of proliferating and differentiating to recover all lost structures functionally. During regeneration, the blastema cells accelerate the cell cycle and duplicate its genome, which is inherently difficult to replicate because of its length and composition, thus being prone to suffer double-strand breaks.

RESULTS

We identified and characterized two remarkable components of the homologous recombination repair pathway (Amex.RAD51 and Amex.MRE11), which were heterologously expressed, biochemically characterized, and inhibited by specific chemicals. These same inhibitors were applied at different time points after amputation to study their effects during limb regeneration. We observed an increase in cellular senescent accompanied by a slight delay in regeneration at 28 days post-amputation regenerated tissues; moreover, inhibitors caused a rise in the double-strand break signaling as a response to the inhibition of the repair mechanisms.

CONCLUSIONS

We confirmed the participation and importance of homologous recombination during limb regeneration. Where the chemical inhibition induces double-strand breaks that lead to DNA damage associated senescence, or in an alternatively way, this damage could be possibly repaired by a different DNA repair pathway, permitting proper regeneration and avoiding senescence. This article is protected by copyright. All rights reserved.

Alternative Animal Models of Aging Research

Most research on mechanisms of aging is being conducted in a very limited number of classical model species, i.e., laboratory mouse (Mus musculus), rat (Rattus norvegicus domestica), the common fruit fly (Drosophila melanogaster) and roundworm (Caenorhabditis elegans). The obvious advantages of using these models are access to resources such as strains with known genetic properties, high-quality genomic and transcriptomic sequencing data, versatile experimental manipulation capabilities including well-established genome editing tools, as well as extensive experience in husbandry. However, this approach may introduce interpretation biases due to the specific characteristics of the investigated species, which may lead to inappropriate, or even false, generalization. For example, it is still unclear to what extent knowledge of aging mechanisms gained in short-lived model organisms is transferable to long-lived species such as humans. In addition, other specific adaptations favoring a long and healthy life from the immense evolutionary toolbox may be entirely missed. In this review, we summarize the specific characteristics of emerging animal models that have attracted the attention of gerontologists, we provide an overview of the available data and resources related to these models, and we summarize important insights gained from them in recent years. The models presented include short-lived ones such as killifish (Nothobranchius furzeri), long-lived ones such as primates (Callithrix jacchus, Cebus imitator, Macaca mulatta), bathyergid mole-rats (Heterocephalus glaber, Fukomys spp.), bats (Myotis spp.), birds, olms (Proteus anguinus), turtles, greenland sharks, bivalves (Arctica islandica), and potentially non-aging ones such as Hydra and Planaria.

Timing Does Matter: Nerve-Mediated HDAC1 Paces the Temporal Expression of Morphogenic Genes During Axolotl Limb Regeneration

Sophisticated axolotl limb regeneration is a highly orchestrated process that requires highly regulated gene expression and epigenetic modification patterns at precise positions and timings. We previously demonstrated two waves of post-amputation expression of a nerve-mediated repressive epigenetic modulator, histone deacetylase 1 (HDAC1), at the wound healing (3 days post-amputation; 3 dpa) and blastema formation (8 dpa onward) stages in juvenile axolotls. Limb regeneration was profoundly inhibited by local injection of an HDAC inhibitor, MS-275, at the amputation sites. To explore the transcriptional response of post-amputation axolotl limb regeneration in a tissue-specific and time course-dependent manner after MS-275 treatment, we performed transcriptome sequencing of the epidermis and soft tissue (ST) at 0, 3, and 8 dpa with and without MS-275 treatment. Gene Ontology (GO) enrichment analysis of each coregulated gene cluster revealed a complex array of functional pathways in both the epidermis and ST. In particular, HDAC activities were required to inhibit the premature elevation of genes related to tissue development, differentiation, and morphogenesis. Further validation by Q-PCR in independent animals demonstrated that the expression of 5 out of 6 development- and regeneration-relevant genes that should only be elevated at the blastema stage was indeed prematurely upregulated at the wound healing stage when HDAC1 activity was inhibited. WNT pathway-associated genes were also prematurely activated under HDAC1 inhibition. Applying a WNT inhibitor to MS-275-treated amputated limbs partially rescued HDAC1 inhibition, resulting in blastema formation defects. We propose that post-amputation HDAC1 expression is at least partially responsible for pacing the expression timing of morphogenic genes to facilitate proper limb regeneration.

From Autonomy to Integration, From Integration to Dynamically Balanced Integrated Co-existence: Non-aging as the Third Stage of Development

Reversible senescence at the cellular level emerged together with tissue specialization in Metazoans. However, this reversibility (ability to permanently rejuvenate) through recapitulation of early stages of development, was originally a part of ontogenesis, since the pressure of integrativeness was not dominant. The complication of specialization in phylogenesis narrowed this “freedom of maneuver”, gradually “truncating” remorphogenesis to local epimorphosis and further up to the complete disappearance of remorphogenesis from the ontogenesis repertoire. This evolutionary trend transformed cellular senescence into organismal aging and any recapitulation of autonomy into carcinogenesis. The crown of specialization, Homo sapiens, completed this post-unicellular stage of development, while in the genome all the potential for the next stage of development, which can be called the stage of balanced coexistence of autonomous and integrative dominants within a single whole. Here, completing the substantiation of the new section of developmental biology, we propose to call it Developmental Biogerontology.

Single-cell RNA-seq reveals novel mitochondria-related musculoskeletal cell populations during adult axolotl limb regeneration process

While the capacity to regenerate tissues or limbs is limited in mammals, including humans, axolotls are able to regrow entire limbs and major organs after incurring a wound. The wound blastema has been extensively studied in limb regeneration. However, due to the inadequate characterization of ECM and cell subpopulations involved in the regeneration process, the discovery of the key drivers for human limb regeneration remains unknown. In this study, we applied large-scale single-cell RNA sequencing to classify cells throughout the adult axolotl limb regeneration process, uncovering a novel regeneration-specific mitochondria-related cluster supporting regeneration through energy providing and the ECM secretion (COL2+) cluster contributing to regeneration through cell-cell interactions signals. We also discovered the dedifferentiation and re-differentiation of the COL1+/COL2+ cellular subpopulation and exposed a COL2-mitochondria subcluster supporting the musculoskeletal system regeneration. On the basis of these findings, we reconstructed the dynamic single-cell transcriptome of adult axolotl limb regenerative process, and identified the novel regenerative mitochondria-related musculoskeletal populations, which yielded deeper insights into the crucial interactions between cell clusters within the regenerative microenvironment.

Towards comparative analyses of salamander limb regeneration

Among tetrapods, only salamanders can regenerate their limbs and tails throughout life. This amazing regenerative ability has attracted the attention of scientists for hundreds of years. Now that large, salamander genomes are beginning to be sequenced for the first time, omics tools and approaches can be used to integrate new perspectives into the study of tissue regeneration. Here we argue the need to move beyond the primary salamander models to investigate regeneration in other species. Salamanders at first glance come across as a phylogenetically conservative group that has not diverged greatly from their ancestors. While salamanders do present ancestral characteristics of basal tetrapods, including the ability to regenerate limbs, data from fossils and data from studies that have tested for species differences suggest there may be considerable variation in how salamanders develop and regenerate their limbs. We review the case for expanded studies of salamander tissue regeneration and identify questions and approaches that are most likely to reveal commonalities and differences in regeneration among species. We also address challenges that confront such an initiative, some of which are regulatory and not scientific. The time is right to gain evolutionary perspective about mechanisms of tissue regeneration from comparative studies of salamander species.

Functional Characterization of the Lin28/let-7 Circuit During Forelimb Regeneration in Ambystoma mexicanum and Its Influence on Metabolic Reprogramming

The axolotl (Ambystoma mexicanum) is a caudate amphibian, which has an extraordinary ability to restore a wide variety of damaged structures by a process denominated epimorphosis. While the origin and potentiality of progenitor cells that take part during epimorphic regeneration are known to some extent, the metabolic changes experienced and their associated implications, remain unexplored. However, a circuit with a potential role as a modulator of cellular metabolism along regeneration is that formed by Lin28/let-7. In this study, we report two Lin28 paralogs and eight mature let-7 microRNAs encoded in the axolotl genome. Particularly, in the proliferative blastema stage amxLin28B is more abundant in the nuclei of blastemal cells, while the microRNAs amx-let-7c and amx-let-7a are most downregulated. Functional inhibition of Lin28 factors increase the levels of most mature let-7 microRNAs, consistent with an increment of intermediary metabolites of the Krebs cycle, and phenotypic alterations in the outgrowth of the blastema. In summary, we describe the primary components of the Lin28/let-7 circuit and their function during axolotl regeneration, acting upstream of metabolic reprogramming events.

Conservation analysis of core cell cycle regulators and their transcriptional behavior during limb regeneration in Ambystoma mexicanum

Ambystoma mexicanum (axolotl) has been one of the major experimental models for the study of regeneration during the past 100 years. Axolotl limb regeneration takes place through a multi-stage and complex developmental process called epimorphosis that involves diverse events of cell reprogramming. Such events start with dedifferentiation of somatic cells and the proliferation of quiescent stem cells to generate a population of proliferative cells called blastema. Once the blastema reaches a mature stage, cells undergo progressive differentiation into the diverse cell lineages that will form the new limb. Such pivotal cell reprogramming phenomena depend on the fine-tuned regulation of the cell cycle in each regeneration stage, where cell populations display specific proliferative capacities and differentiation status. The axolotl genome has been fully sequenced and released recently, and diverse RNA-seq approaches have also been generated, enabling the identification and conservatory analysis of core cell cycle regulators in this species. We report here our results from such analyses and present the transcriptional behavior of key regulatory factors during axolotl limb regeneration. We also found conserved protein interactions between axolotl Cyclin Dependent Kinases 2, 4 and 6 and Cyclins type D and E. Canonical CYC-CDK interactions that play major roles in modulating cell cycle progression in eukaryotes.

TimeMeter assesses temporal gene expression similarity and identifies differentially progressing genes

Comparative time series transcriptome analysis is a powerful tool to study development, evolution, aging, disease progression and cancer prognosis. We develop TimeMeter, a statistical method and tool to assess temporal gene expression similarity, and identify differentially progressing genes where one pattern is more temporally advanced than the other. We apply TimeMeter to several datasets, and show that TimeMeter is capable of characterizing complicated temporal gene expression associations. Interestingly, we find: (i) the measurement of differential progression provides a novel feature in addition to pattern similarity that can characterize early developmental divergence between two species; (ii) genes exhibiting similar temporal patterns between human and mouse during neural differentiation are under strong negative (purifying) selection during evolution; (iii) analysis of genes with similar temporal patterns in mouse digit regeneration and axolotl blastema differentiation reveals common gene groups for appendage regeneration with potential implications in regenerative medicine.

BMP signaling is essential for sustaining proximo-distal progression in regenerating axolotl limbs

Amputation of a salamander limb triggers a regeneration process that is perfect. A limited number of genes have been studied in this context and even fewer have been analyzed functionally. In this work, we use the BMP signaling inhibitor LDN193189 on Ambystoma mexicanum to explore the role of BMPs in regeneration. We find that BMP signaling is required for proper expression of various patterning genes and that its inhibition causes major defects in the regenerated limbs. Fgf8 is downregulated when BMP signaling is blocked, but ectopic injection of either human or axolotl protein did not rescue the defects. By administering LDN193189 treatments at different time points during regeneration, we show clearly that limb regeneration progresses in a proximal to distal fashion. This demonstrates that BMPs play a major role in patterning of regenerated limbs and that regeneration is a progressive process like development.

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.

von Willebrand factor D and EGF domains is an evolutionarily conserved and required feature of blastemas capable of multitissue appendage regeneration

Regenerative ability varies tremendously across species. A common feature of regeneration of appendages such as limbs, fins, antlers, and tails is the formation of a blastema—a transient structure that houses a pool of progenitor cells that can regenerate the missing tissue. We have identified the expression of von Willebrand factor D and EGF domains (vwde) as a common feature of blastemas capable of regenerating limbs and fins in a variety of highly regenerative species, including axolotl (Ambystoma mexicanum), lungfish (Lepidosiren paradoxa), and Polpyterus (Polypterus senegalus). Further, vwde expression is tightly linked to the ability to regenerate appendages in Xenopus laevis. Functional experiments demonstrate a requirement for vwde in regeneration and indicate that Vwde is a potent growth factor in the blastema. These data identify a key role for vwde in regenerating blastemas and underscore the power of an evolutionarily informed approach for identifying conserved genetic components of regeneration.

vwde expression is a common feature of blastemas capable of fin and limb regeneration.

vwde expression is tightly tied to regeneration‐competency.

vwde is required for axolotl limb regeneration, with transient knockdown resulting in severe endpoint phenotypes.

Eya2 promotes cell cycle progression by regulating DNA damage response during vertebrate limb regeneration

How salamanders accomplish progenitor cell proliferation while faithfully maintaining genomic integrity and regenerative potential remains elusive. Here we found an innate DNA damage response mechanism that is evident during blastema proliferation (early- to late-bud) and studied its role during tissue regeneration by ablating the function of one of its components, Eyes absent 2. In eya2 mutant axolotls, we found that DNA damage signaling through the H2AX histone variant was deregulated, especially within the proliferating progenitors during limb regeneration. Ultimately, cell cycle progression was impaired at the G1/S and G2/M transitions and regeneration rate was reduced. Similar data were acquired using acute pharmacological inhibition of the Eya2 phosphatase activity and the DNA damage checkpoint kinases Chk1 and Chk2 in wild-type axolotls. Together, our data indicate that highly-regenerative animals employ a robust DNA damage response pathway which involves regulation of H2AX phosphorylation via Eya2 to facilitate proper cell cycle progression upon injury.

Remembering where we are: Positional information in salamander limb regeneration

Fifty years ago, Lewis Wolpert defined an important question in developmental biology: how are cell fates determined by the positions of cells within a system? He proposed that cells retain positional values as if they lie within a coordinate system and that the interpretation of these values produces patterns in development. He referred to this concept as positional information. Though initially controversial, this concept of positional information has proven to be profoundly influential in developmental biology. One area in which the influence of Wolpert's theoretical work can be clearly demonstrated is the study of limb regeneration in salamanders. Here, we review the work in limb regeneration leading up to Wolpert defining the concept of positional information and how his theory has guided regeneration research over the subsequent 50 years.

Multiplex CRISPR/Cas screen in regenerating haploid limbs of chimeric Axolotls

Axolotls and other salamanders can regenerate entire limbs after amputation as adults, and much recent effort has sought to identify the molecular programs controlling this process. While targeted mutagenesis approaches like CRISPR/Cas9 now permit gene-level investigation of these mechanisms, genetic screening in the axolotl requires an extensive commitment of time and space. Previously, we quantified CRISPR/Cas9-generated mutations in the limbs of mosaic mutant axolotls before and after regeneration and found that the regenerated limb is a highfidelity replicate of the original limb (Flowers et al. 2017). Here, we circumvent aforementioned genetic screening limitations and present methods for a multiplex CRISPR/Cas9 haploid screen in chimeric axolotls (MuCHaChA), which is a novel platform for haploid genetic screening in animals to identify genes essential for limb regeneration.

Lag penalized weighted correlation for time series clustering

BACKGROUND

The similarity or distance measure used for clustering can generate intuitive and interpretable clusters when it is tailored to the unique characteristics of the data. In time series datasets generated with high-throughput biological assays, measurements such as gene expression levels or protein phosphorylation intensities are collected sequentially over time, and the similarity score should capture this special temporal structure.

RESULTS

We propose a clustering similarity measure called Lag Penalized Weighted Correlation (LPWC) to group pairs of time series that exhibit closely-related behaviors over time, even if the timing is not perfectly synchronized. LPWC aligns time series profiles to identify common temporal patterns. It down-weights aligned profiles based on the length of the temporal lags that are introduced. We demonstrate the advantages of LPWC versus existing time series and general clustering algorithms. In a simulated dataset based on the biologically-motivated impulse model, LPWC is the only method to recover the true clusters for almost all simulated genes. LPWC also identifies clusters with distinct temporal patterns in our yeast osmotic stress response and axolotl limb regeneration case studies.

CONCLUSIONS

LPWC achieves both of its time series clustering goals. It groups time series with correlated changes over time, even if those patterns occur earlier or later in some of the time series. In addition, it refrains from introducing large shifts in time when searching for temporal patterns by applying a lag penalty. The LPWC R package is available at https://github.com/gitter-lab/LPWC and CRAN under a MIT license.

Integrative Analysis of Axolotl Gene Expression Data from Regenerative and Wound Healing Limb Tissues

Axolotl (Ambystoma mexicanum) is a urodele amphibian endowed with remarkable regenerative capacities manifested in scarless wound healing and restoration of amputated limbs, which makes it a powerful experimental model for regenerative biology and medicine. Previous studies have utilized microarrays and RNA-Seq technologies for detecting differentially expressed (DE) genes in different phases of the axolotl limb regeneration. However, sufficient consistency may be lacking due to statistical limitations arising from intra-laboratory analyses. This study aims to bridge such gaps by performing an integrative analysis of publicly available microarray and RNA-Seq data from axolotl limb samples having comparable study designs using the “merging” method. A total of 351 genes were found DE in regenerative samples compared to the control in data of both technologies, showing an adjusted p-value < 0.01 and log fold change magnitudes >1. Downstream analyses illustrated consistent correlations of the directionality of DE genes within and between data of both technologies, as well as concordance with the literature on regeneration related biological processes. qRT-PCR analysis validated the observed expression level differences of five of the top DE genes. Future studies may benefit from the utilized concept and approach for enhanced statistical power and robust discovery of biomarkers of regeneration.

Advancements to the Axolotl Model for Regeneration and Aging

Loss of regenerative capacity is a normal part of aging. However, some organisms, such as the Mexican axolotl, retain striking regenerative capacity throughout their lives. Moreover, the development of age-related diseases is rare in this organism. In this review, we will explore how axolotls are used as a model system to study regenerative processes, the exciting new technological advancements now available for this model, and how we can apply the lessons we learn from studying regeneration in the axolotl to understand, and potentially treat, age-related decline in humans.

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.

Longitudinal 16S rRNA data derived from limb regenerative tissue samples of axolotl Ambystoma mexicanum

The Mexican axolotl (Ambystoma mexicanum) is a critically endangered species and a fruitful amphibian model for regenerative biology. Despite growing body of research on the cellular and molecular biology of axolotl limb regeneration, microbiological aspects of this process remain poorly understood. Here, we describe bacterial 16S rRNA amplicon dataset derived from axolotl limb tissue samples in the course of limb regeneration. The raw data was obtained by sequencing V3–V4 region of 16S rRNA gene and comprised 14,569,756 paired-end raw reads generated from 21 samples. Initial data analysis using DADA2 pipeline resulted in amplicon sequence variant (ASV) table containing a total of ca. 5.9 million chimera-removed, high-quality reads and a median of 296,971 reads per sample. The data constitute a useful resource for the research on the microbiological aspects of axolotl limb regeneration and will also broadly facilitate comparative studies in the developmental and conservation biology of this critically endangered species.

Design Type(s)	time series design • organism development design
Measurement Type(s)	rRNA 16S
Technology Type(s)	DNA sequencing
Factor Type(s)	biological replicate • temporal_interval
Sample Characteristic(s)	Ambystoma mexicanum • limb • laboratory environment

Machine-accessible metadata file describing the reported data (ISA-Tab format)

More Than Just a Bandage: Closing the Gap Between Injury and Appendage Regeneration

The remarkable regenerative capabilities of amphibians have captured the attention of biologists for centuries. The frogs Xenopus laevis and Xenopus tropicalis undergo temporally restricted regenerative healing of appendage amputations and spinal cord truncations, injuries that are both devastating and relatively common in human patients. Rapidly expanding technological innovations have led to a resurgence of interest in defining the factors that enable regenerative healing, and in coupling these factors to human therapeutic interventions. It is well-established that early embryonic signaling pathways are critical for growth and patterning of new tissue during regeneration. A growing body of research now indicates that early physiological injury responses are also required to initiate a regenerative program, and that these differ in regenerative and non-regenerative contexts. Here we review recent insights into the biophysical, biochemical, and epigenetic processes that underlie regenerative healing in amphibians, focusing particularly on tail and limb regeneration in Xenopus. We also discuss the more elusive potential mechanisms that link wounding to tissue growth and patterning.

Blastemal progenitors modulate immune signaling during early limb regeneration

Blastema formation, a hallmark of limb regeneration, requires proliferation and migration of progenitors to the amputation plane. Although blastema formation has been well described, the transcriptional programs that drive blastemal progenitors remain unknown. We transcriptionally profiled dividing and non-dividing cells in regenerating stump tissues, as well as the wound epidermis, during early axolotl limb regeneration. Our analysis revealed unique transcriptional signatures of early dividing cells and, unexpectedly, repression of several core developmental signaling pathways in early regenerating stump tissues. We further identify an immunomodulatory role for blastemal progenitors through interleukin 8 (IL-8), a highly expressed cytokine in subpopulations of early blastemal progenitors. Ectopic il-8 expression in non-regenerating limbs induced myeloid cell recruitment, while IL-8 knockdown resulted in defective myeloid cell retention during late wound healing, delaying regeneration. Furthermore, the il-8 receptor cxcr-1/2 was expressed in myeloid cells, and inhibition of CXCR-1/2 signaling during early stages of limb regeneration prevented regeneration. Altogether, our findings suggest that blastemal progenitors are active early mediators of immune support, and identify CXCR-1/2 signaling as an important immunomodulatory pathway during the initiation of regeneration.

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.

Using transcriptomics to enable a plethodontid salamander (Bolitoglossa ramosi) for limb regeneration research

Background

Tissue regeneration is widely distributed across the tree of life. Among vertebrates, salamanders possess an exceptional ability to regenerate amputated limbs and other complex structures. Thus far, molecular insights about limb regeneration have come from a relatively limited number of species from two closely related salamander families. To gain a broader perspective on the molecular basis of limb regeneration and enhance the molecular toolkit of an emerging plethodontid salamander (Bolitoglossa ramosi), we used RNA-Seq to generate a de novo reference transcriptome and identify differentially expressed genes during limb regeneration.

Results

Using paired-end Illumina sequencing technology and Trinity assembly, a total of 433,809 transcripts were recovered and we obtained functional annotation for 142,926 non-redundant transcripts of the B. ramosi de novo reference transcriptome. Among the annotated transcripts, 602 genes were identified as differentially expressed during limb regeneration. This list was further processed to identify a core set of genes that exhibit conserved expression changes between B. ramosi and the Mexican axolotl (Ambystoma mexicanum), and presumably their common ancestor from approximately 180 million years ago.

Conclusions

We identified genes from B. ramosi that are differentially expressed during limb regeneration, including multiple conserved protein-coding genes and possible putative species-specific genes. Comparative analyses reveal a subset of genes that show similar patterns of expression with ambystomatid species, which highlights the importance of developing comparative gene expression data for studies of limb regeneration among salamanders.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5076-0) contains supplementary material, which is available to authorized users.

Comparative transcriptomics of limb regeneration: Identification of conserved expression changes among three species of Ambystoma

Transcriptome studies are revealing the complex gene expression basis of limb regeneration in the primary salamander model - Ambystoma mexicanum (axolotl). To better understand this complexity, there is need to extend analyses to additional salamander species. Using microarray and RNA-Seq, we performed a comparative transcriptomic study using A. mexicanum and two other ambystomatid salamanders: A. andersoni, and A. maculatum. Salamanders were administered forelimb amputations and RNA was isolated and analyzed to identify 405 non-redundant genes that were commonly, differentially expressed 24 h post amputation. Many of the upregulated genes are predicted to function in wound healing and developmental processes, while many of the downregulated genes are typically expressed in muscle. The conserved transcriptional changes identified in this study provide a high-confidence dataset for identifying factors that simultaneous orchestrate wound healing and regeneration processes in response to injury, and more generally for identifying genes that are essential for salamander limb regeneration.

The journey of cells through regeneration

The process of building an organ, appendage, or organism requires the precise coordination of cells in space and time. Regeneration of those same tissues adds an additional element of complexity, emerging from the chaos of disease or injury to build a mass of progenitors from mature tissue. Translating insights from natural examples of tissue regeneration into engineered regenerative therapies requires a deep understanding of the journey of a cell directly following injury to its contribution to functional, scaled replacement tissue. Here we step through the chronological phases of regeneration and highlight emerging work that brings us closer to elucidating the unique intrinsic and extrinsic properties of cells during epimorphic regeneration.

Changes in the appendicular skeleton during metamorphosis in the axolotl salamander (Ambystoma mexicanum)

Axolotl salamanders (Ambystoma mexicanum) remain aquatic in their natural state, during which biomechanical forces on their diarthrodial limb joints are likely reduced relative to salamanders living on land. However, even as sexually mature adults, these amphibians can be induced to metamorphose into a weight-bearing terrestrial stage by environmental stress or the exogenous administration of thyroxine hormone. In some respects, this aquatic to terrestrial transition of axolotl salamanders through metamorphosis may model developmental and changing biomechanical skeletal forces in mammals during the prenatal to postnatal transition at birth and in the early postnatal period. To assess differences in the appendicular skeleton as a function of metamorphosis, anatomical and gene expression parameters were compared in skeletal tissues between aquatic and terrestrial axolotls that were the same age and genetically full siblings. The length of long bones and area of cuboidal bones in the appendicular skeleton, as well as the cellularity of cartilaginous and interzone tissues of femorotibial joints were generally higher in aquatic axolotls compared with their metamorphosed terrestrial siblings. A comparison of steady-state mRNA transcripts encoding aggrecan core protein (ACAN), type II collagen (COL2A1), and growth and differentiation factor 5 (GDF5) in femorotibial cartilaginous and interzone tissues did not reveal any significant differences between aquatic and terrestrial axolotls.

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.

Virtual Genome Walking across the 32 Gb Ambystoma mexicanum genome; assembling gene models and intronic sequence

Large repeat rich genomes present challenges for assembly using short read technologies. The 32 Gb axolotl genome is estimated to contain ~19 Gb of repetitive DNA making an assembly from short reads alone effectively impossible. Indeed, this model species has been sequenced to 20× coverage but the reads could not be conventionally assembled. Using an alternative strategy, we have assembled subsets of these reads into scaffolds describing over 19,000 gene models. We call this method Virtual Genome Walking as it locally assembles whole genome reads based on a reference transcriptome, identifying exons and iteratively extending them into surrounding genomic sequence. These assemblies are then linked and refined to generate gene models including upstream and downstream genomic, and intronic, sequence. Our assemblies are validated by comparison with previously published axolotl bacterial artificial chromosome (BAC) sequences. Our analyses of axolotl intron length, intron-exon structure, repeat content and synteny provide novel insights into the genic structure of this model species. This resource will enable new experimental approaches in axolotl, such as ChIP-Seq and CRISPR and aid in future whole genome sequencing efforts. The assembled sequences and annotations presented here are freely available for download from https://tinyurl.com/y8gydc6n . The software pipeline is available from https://github.com/LooseLab/iterassemble .

Regeneration in distantly related species: common strategies and pathways

While almost all animals are able to at least partially replace some lost parts, regeneration abilities vary considerably across species. Here we study gene expression patterns in distantly related species to investigate conserved regeneration strategies. To this end, we collect from the literature transcriptomic data obtained during the regeneration of three species (Hydra magnipapillata, Schmidtea mediterranea, and Apostichopus japonicus), and compare them with gene expression during regeneration in vertebrates and mammals. This allows us to identify a common set of differentially expressed genes and relevant shared pathways that are conserved across species during the early stage of the regeneration process. We also find a set of differentially expressed genes that in mammals are associated to the presence of macrophages and to the epithelial–mesenchymal transition. This suggests that features of the sophisticated wound healing strategy of mammals are already observable in earlier emerging metazoans.

All animals capable of regenerating a lost body part, from an organ or a limb to the whole organism, use a common set of genes. This is the striking discovery of a team of researchers from the Center for Complexity and Biosystems of the University of Milan, led by Caterina La Porta. They analyzed the genetic activity in regenerating tissues from widely different species—from hydra to mice. They found that some of the genes active at the beginning of the regeneration process are the same in very different species, including mammals which have lost this function during evolution, except for the restoration of the liver. The discovery of this shared genetic signature is of great importance to understand how regeneration evolved and could be useful for future regeneration therapies.

A Tissue-Mapped Axolotl De Novo Transcriptome Enables Identification of Limb Regeneration Factors

Mammals have extremely limited regenerative capabilities; however, axolotls are profoundly regenerative and can replace entire limbs. The mechanisms underlying limb regeneration remain poorly understood, partly because the enormous and incompletely sequenced genomes of axolotls have hindered the study of genes facilitating regeneration. We assembled and annotated a de novo transcriptome using RNA-sequencing profiles for a broad spectrum of tissues that is estimated to have near-complete sequence information for 88% of axolotl genes. We devised expression analyses that identified the axolotl orthologs of cirbp and kazald1 as highly expressed and enriched in blastemas. Using morpholino anti-sense oligonucleotides, we find evidence that cirbp plays a cytoprotective role during limb regeneration whereas manipulation of kazald1 expression disrupts regeneration. Our transcriptome and annotation resources greatly complement previous transcriptomic studies and will be a valuable resource for future research in regenerative biology.

Transcriptome dynamics along axolotl regenerative development are consistent with an extensive reduction in gene expression heterogeneity in dedifferentiated cells

Although in recent years the study of gene expression variation in the absence of genetic or environmental cues or gene expression heterogeneity has intensified considerably, many basic and applied biological fields still remain unaware of how useful the study of gene expression heterogeneity patterns might be for the characterization of biological systems and/or processes. Largely based on the modulator effect chromatin compaction has for gene expression heterogeneity and the extensive changes in chromatin compaction known to occur for specialized cells that are naturally or artificially induced to revert to less specialized states or dedifferentiate, I recently hypothesized that processes that concur with cell dedifferentiation would show an extensive reduction in gene expression heterogeneity. The confirmation of the existence of such trend could be of wide interest because of the biomedical and biotechnological relevance of cell dedifferentiation-based processes, i.e., regenerative development, cancer, human induced pluripotent stem cells, or plant somatic embryogenesis. Here, I report the first empirical evidence consistent with the existence of an extensive reduction in gene expression heterogeneity for processes that concur with cell dedifferentiation by analyzing transcriptome dynamics along forearm regenerative development in Ambystoma mexicanum or axolotl. Also, I briefly discuss on the utility of the study of gene expression heterogeneity dynamics might have for the characterization of cell dedifferentiation-based processes, and the engineering of tools that afforded better monitoring and modulating such processes. Finally, I reflect on how a transitional reduction in gene expression heterogeneity for dedifferentiated cells can promote a long-term increase in phenotypic heterogeneity following cell dedifferentiation with potential adverse effects for biomedical and biotechnological applications.