Two origins of blastemal progenitors define blastemal regeneration of zebrafish lower jaw

Zebrafish Models for Skeletal and Extraskeletal Osteogenesis Imperfecta Features: Unveiling Pathophysiology and Paving the Way for Drug Discovery

In the last decades, the easy genetic manipulation, the external fertilization, the high percentage of homology with human genes and the reduced husbandry costs compared to rodents, made zebrafish a valid model for studying human diseases and for developing new therapeutical strategies. Since zebrafish shares with mammals the same bone cells and ossification types, it became widely used to dissect mechanisms and possible new therapeutic approaches in the field of common and rare bone diseases, such as osteoporosis and osteogenesis imperfecta (OI), respectively. OI is a heritable skeletal disorder caused by defects in gene encoding collagen I or proteins/enzymes necessary for collagen I synthesis and secretion. Nevertheless, OI patients can be also characterized by extraskeletal manifestations such as dentinogenesis imperfecta, muscle weakness, cardiac valve and pulmonary abnormalities and skin laxity. In this review, we provide an overview of the available zebrafish models for both dominant and recessive forms of OI. An updated description of all the main similarities and differences between zebrafish and mammal skeleton, muscle, heart and skin, will be also discussed. Finally, a list of high- and low-throughput techniques available to exploit both larvae and adult OI zebrafish models as unique tools for the discovery of new therapeutic approaches will be presented.

Zebrafish Feed Intake: A Systematic Review for Standardizing Feeding Management in Laboratory Conditions

Zebrafish are one of the most used animal models in biological research and a cost-effective alternative to rodents. Despite this, nutritional requirements and standardized feeding protocols have not yet been established for this species. This is important to avoid nutritional effects on experimental outcomes, and especially when zebrafish models are used in preclinical studies, as many diseases have nutritional confounding factors. A key aspect of zebrafish nutrition is related to feed intake, the amount of feed ingested by each fish daily. With the goal of standardizing feeding protocols among the zebrafish community, this paper systematically reviews the available data from 73 studies on zebrafish feed intake, feeding regimes (levels), and diet composition. Great variability was observed regarding diet composition, especially regarding crude protein (mean 44.98 ± 9.87%) and lipid content (9.91 ± 5.40%). Interestingly, the gross energy levels of the zebrafish diets were similar across the reviewed studies (20.39 ± 2.10 kilojoules/g of feed). In most of the reviewed papers, fish received a predetermined quantity of feed (feed supplied). The authors fed the fish according to the voluntary intake and then calculated feed intake (FI) in only 17 papers. From a quantitative point of view, FI was higher than when a fixed quantity (pre-defined) of feed was supplied. Also, the literature showed that many biotic and abiotic factors may affect zebrafish FI. Finally, based on the FI data gathered from the literature, a new feeding protocol is proposed. In summary, a daily feeding rate of 9-10% of body weight is proposed for larvae, whereas these values are equal to 6-8% for juveniles and 5% for adults when a dry feed with a proper protein and energy content is used.

Hemifacial microsomia is linked to a rare homozygous variant V162I in FRK and validated in zebrafish

OBJECTIVES

Hemifacial microsomia (HFM) is a common birth defect involving the first and second branchial arch derivatives. Although several chromosomal abnormalities and causal gene variants have been identified, genetic etiologies in a majority of cases with HFM remain unknown. This study aimed to identify genetic mutations in affected individuals with HFM.

METHODS

Whole-exome sequencing and bioinformatics analysis were performed for 16 affected individuals and their family members. Sanger sequencing was applied for confirmation of selected mutations. Zebrafish embryos were used for in situ hybridization of candidate gene, microinjection with antisense morpholino, and cartilage staining.

RESULTS

A homozygous missense mutation (c.484G > A; p.V162I) in the FRK gene was identified in an 18-year-old girl with HFM and dental abnormalities. Heterozygous mutation of this mutation was identified in her parents, who are first cousins in a consanguineous family. FRK is highly expressed in the Meckel's cartilage during embryonic development in mouse and zebrafish. Knockdown of frk in zebrafish showed a lower length and width ratio of Meckel's cartilage, abnormal mandibular jaw joint, and disorganized ceratobranchial cartilage and bone.

CONCLUSIONS

We identified a recessive variant in the FRK gene as a novel candidate gene for a patient with HFM and mandibular hypoplasia and revealed its effects on craniofacial and embryonic development in zebrafish.

Cartilage regeneration in zebrafish depends on Nrg1/ErbB signaling pathway

Objective: Cartilage, as the majority of adult mammalian tissues, has limited regeneration capacity. Cartilage degradation consecutive to joint injury or aging then leads to irreversible joint damage and diseases. In contrast, several vertebrate species such as the zebrafish have the remarkable capacity to spontaneously regenerate skeletal structures after severe injuries. The objective of our study was to test the regenerative capacity of Meckel's cartilage (MC) upon mechanical injury in zebrafish and to identify the mechanisms underlying this process. Methods and Results: Cartilage regenerative capacity in zebrafish larvae was investigated after mechanical injuries of the lower jaw MC in TgBAC(col2a1a:mCherry), to visualize the loss and recovery of cartilage. Confocal analysis revealed the formation of new chondrocytes and complete regeneration of MC at 14 days post-injury (dpi) via chondrocyte cell cycle re-entry and proliferation of pre-existing MC chondrocytes near the wound. Through expression analyses, we showed an increase of nrg1 expression in the regenerating lower jaw, which also expresses Nrg1 receptors, ErbB3 and ErbB2. Pharmacological inhibition of the ErbB pathway and specific knockdown of Nrg1 affected MC regeneration indicating the pivotal role of this pathway for cartilage regeneration. Finally, addition of exogenous NRG1 in an in vitro model of osteoarthritic (OA)-like chondrocytes induced by IL1β suggests that Nrg1/ErbB pathway is functional in mammalian chondrocytes and alleviates the increased expression of catabolic markers characteristic of OA-like chondrocytes. Conclusion: Our results show that the Nrg1/ErbB pathway is required for spontaneous cartilage regeneration in zebrafish and is of interest to design new therapeutic approaches to promote cartilage regeneration in mammals.

Skeletal Biology and Disease Modeling in Zebrafish

Zebrafish are teleosts (bony fish) that share with mammals a common ancestor belonging to the phylum Osteichthyes, from which their endoskeletal systems have been inherited. Indeed, teleosts and mammals have numerous genetically conserved features in terms of skeletal elements, ossification mechanisms, and bone matrix components in common. Yet differences related to bone morphology and function need to be considered when investigating zebrafish in skeletal research. In this review, we focus on zebrafish skeletal architecture with emphasis on the morphology of the vertebral column and associated anatomical structures. We provide an overview of the different ossification types and osseous cells in zebrafish and describe bone matrix composition at the microscopic tissue level with a focus on assessing mineralization. Processes of bone formation also strongly depend on loading in zebrafish, as we elaborate here. Furthermore, we illustrate the high regenerative capacity of zebrafish bones and present some of the technological advantages of using zebrafish as a model. We highlight zebrafish axial and fin skeleton patterning mechanisms, metabolic bone disease such as after immunosuppressive glucocorticoid treatment, as well as osteogenesis imperfecta (OI) and osteopetrosis research in zebrafish. We conclude with a view of why larval zebrafish xenografts are a powerful tool to study bone metastasis. © 2021 American Society for Bone and Mineral Research (ASBMR).

Tissue regeneration during lower jaw restoration in zebrafish shows some features of epimorphic regeneration

Zebrafish have the ability to regenerate skeletal structures, including the fin, skull roof, and jaw. Although fin regeneration proceeds by epimorphic regeneration, it remains unclear whether this process is involved in other skeletal regeneration in zebrafish. Initially in epimorphic regeneration, the wound epidermis covers the wound surface. Subsequently, the blastema, an undifferentiated mesenchymal mass, forms beneath the epidermis. In the present study, we re-examined the regeneration of the zebrafish lower jaw in detail, and investigated whether epimorphic regeneration is involved in this process. We performed amputation of the lower jaw at two different positions; the proximal level (presence of Meckel's cartilage) and the distal level (absence of Meckel's cartilage). In both manipulations, a blastema-like cellular mass was initially formed. Subsequently, cartilaginous aggregates were formed in this mass. In the proximal amputation, the cartilaginous aggregates were then fused with Meckel's cartilage and remained as a skeletal component of the regenerated jaw, whereas in the distal amputation, the cartilaginous aggregates disappeared as regeneration progressed. Two molecules that were observed during epimorphic regeneration, Laminin and msxb, were expressed in the regenerating lower jaw, although the domain of msxb expression was out of the main plain of the aggregate formation. Administration of an inhibitor of Wnt/β-catenin signaling, a pathway associated with epimorphic regeneration, showed few effects on lower jaw regeneration. Our finding suggests that skeletal regeneration of the lower jaw mainly progresses through tissue regeneration that is dependent on the position in the jaw, and epimorphic regeneration plays an adjunctive role in this regeneration.

Ability of prebiotic polysaccharides to activate a HIF1α-antimicrobial peptide axis determines liver injury risk in zebrafish

Natural polysaccharides have received much attention for their ability to ameliorate hepatic steatosis induced by high-fat diet. However, the potential risks of their use have been less investigated. Here, we show that the exopolysaccharides (EPS) from Lactobacillus rhamnosus GG (LGG) and L. casei BL23 reduce hepatic steatosis in zebrafish fed a high-fat diet, while BL23 EPS, but not LGG EPS, induce liver inflammation and injury. This is due to the fact that BL23 EPS induces gut microbial dysbiosis, while LGG EPS promotes microbial homeostasis. We find that LGG EPS, but not BL23 EPS, can directly activate intestinal HIF1α, and increased HIF1α boosts local antimicrobial peptide expression to facilitate microbial homeostasis, explaining the distinct compositions of LGG EPS- and BL23 EPS-associated microbiota. Finally, we find that liver injury risk is not confined to Lactobacillus-derived EPS but extends to other types of commonly used natural polysaccharides, depending on their HIF1α activation efficiency.

Evolutionarily conserved vasa identifies embryonic and gonadal germ cells in spinyhead croaker Collichthys lucidus

In this study, a 2198 bp full-length cDNA of spinyhead croaker Collichthys lucidus vasa gene encoding 616 amino-acid residues was obtained. Multiple alignment revealed that C. lucidus vasa has eight conserved characteristic motifs of the DEAD box protein family and has the highest identity to large yellow croaker Larimichthys croceas. Reverse-transcription (RT)-PCR and Western blot analyses indicated that the vasa messenger (m)RNA and Vasa protein are specifically expressed in the gonads in both sexes. In situ hybridisation (ISH) demonstrated that vasa RNA is exclusively detected in the germ cells in C. lucidus gonads and its temporospatial expression reveals a dynamic pattern during oogenesis. Surprisingly, C. lucidus vasa 3'UTR can direct stable and specific GFP expression in the primordial germ cells (PGC) of medaka Oryzias latipes embryos. Taken together, these results suggest that because C. lucidus vasa expression delineates critical stages of oogenesis, it may be a useful molecular marker for the identification of gonadal germ cells, facilitating the isolation and utilization of germ cells in future study.

Regenerative Models for the Integration and Regeneration of Head Skeletal Tissues

Disease of, or trauma to, the human jaw account for thousands of reconstructive surgeries performed every year. One of the most popular and successful treatment options in this context involves the transplantation of bone tissue from a different anatomical region into the affected jaw. Although, this method has been largely successful, the integration of the new bone into the existing bone is often imperfect, and the integration of the host soft tissues with the transplanted bone can be inconsistent, resulting in impaired function. Unlike humans, several vertebrate species, including fish and amphibians, demonstrate remarkable regenerative capabilities in response to jaw injury. Therefore, with the objective of identifying biological targets to promote and engineer improved outcomes in the context of jaw reconstructive surgery, we explore, compare and contrast the natural mechanisms of endogenous jaw and limb repair and regeneration in regenerative model organisms. We focus on the role of different cell types as they contribute to the regenerating structure; how mature cells acquire plasticity in vivo; the role of positional information in pattern formation and tissue integration, and limitations to endogenous regenerative and repair mechanisms.

Robust and local positional information within a fin ray directs fin length during zebrafish regeneration

It has been proposed that cells are regulated to form specific morphologies and sizes according to positional information. However, the entity and nature of positional information have not been fully understood yet. The zebrafish caudal fin has a characteristic V-shape; dorsal and ventral fin rays are longer than the central ones. This fin shape regenerates irrespective of the sites or shape of fin amputation. It is thought that reformation of tissue occurs according to positional information. In this study, we developed a novel transplantation procedure for grafting a whole fin ray to an ectopic position and examined whether the information that specifies fin length exists within each fin ray. Intriguingly, when long and short fin rays were swapped, they regenerated to form longer or shorter fin rays than the adjacent host fin rays, respectively. Further, the abnormal fin ray lengths were maintained for a long time, more than 5 months, and after further re-amputation. In contrast to intra-fin grafting, when fin ray grafting was performed between fish, cells in the grafts disappeared due to immune rejection, and the grafted fin rays adapted to the host position to form a normal fin. Together, our data suggest that the information that directs fin length does exist in cells within a single fin ray and that it has a robust property-it is stable for a long time and is hard to rewrite. Our study highlighted a novel positional information mechanism for directing regenerating fin length.

Dietary Safety Assessment of Flk1-Transgenic Fish

Genetic engineering, also called genetic modification, is facing with growing demands of aquaculture and aquatic products. Although various genetically modified (GM) aquatics have been generated, it is important to evaluate biosafety of GM organisms on the human health before entering into our food chain. For this purpose, we establish a zebrafish wild adult feeding Flk1-transgenic larvae model to examine the predatory fish's histology in multiple tissues, and the global gene expression profile in the liver. 180 days of feeding trial show that there are no significantly morphological changes in intestine, liver, kidney, and sex gonads between fish fed with Flk1 transgenic fish diet (TFD) and fish fed with regular food meal (RFM). However, a characteristic skin spot and autofluorescence increase in the theca of follicle are observed in F1 generation of TFD fish. Liver RNA-sequencing analyses demonstrate that 53 out of 56712 genes or isoforms are differentially transcribed, and mostly involved in proteolysis in extracellular region. According to GO enrichment terms these deregulated genes function in catalytic activity, steroid storing, lipid metabolic process and N-Glycan biosynthesis. These results suggest that a long term of Flk1-transgenic fish diet could alter certain metabolic pathways and possibly cause related tissue deformation. Compared to the previous reports, our feasible transgenic dietary assess system could evaluate subchronic and potential health impact of transgenic fish diet by combining multi-tissue histology and liver transcriptome analyses.

Oral-Facial Tissue Reconstruction in the Regenerative Axolotl

Absence of large amounts of orofacial tissues caused by cancerous resections, congenital defects, or trauma results in sequelae such as dysphagia and noticeable scars. Oral-neck tissue regeneration was studied in the axolotl (regenerative amphibian) following a 2.5-mm punch biopsy that simultaneously removed skin, connective tissue, muscle, and cartilage in the tongue and intermandibular region. The untreated wound was studied macroscopically and histologically at 17 different time points ranging from 0 to180 days (N = 120 axolotls). At 12 hr, the wound's surface was smoothened and within 1mm, internal lingual muscular modifications occurred; at the same distance, between days 4-7 lingual muscle degradation was complete. Immunofluorescence indicates complete keratinocytes migration by 48 hr. These cells with epidermal Leydig cells, appearing yellow, lead the chin's deep tissue outgrowth until its closure on the 14th day. Regeneration speeds varied and peaked in time for each tissue, (1) deep chin 84.3 μm/hr from 24 to 96 hr, (2) superficial chin 71.1 μm/hr from 7-14 days, and (3) tongue 86.0 μm/hr between 48 hr and 7 days. Immunofluorescence to Col IV showed basement membrane reconnected between days 30-45 coinciding with the chin's dermal tissue's surface area recovery. New muscle appeared at 21 days and was always preceded by the formation of a collagen bed. Both chin tissues regain all surface area and practically all components while the lingual structure lacks some content but is generally similar to the original. The methodology and high-resolution observations described here are the first of its kind for this animal model and could serve as a basis for future studies in oral and facial regenerative research.

SIRT1 Functions as a Negative Regulator of Eukaryotic Poly(A)RNA Transport

Most eukaryotic mRNAs are polyadenylated in the nucleus, and the poly(A)-tail is required for efficient mRNA export and translation. However, mechanisms governing mRNA transport remain unclear. Here, we report that the nicotinamide adenine dinucleotide (NAD)-dependent deacetylase SIRT1 acts as an energy sensor and negatively regulates poly(A)RNA transport via deacetylating a poly(A)-binding protein, PABP1. Upon energy starvation, SIRT1 interacts with and deacetylates PABP1 and deactivates its poly(A)RNA binding, leading to nuclear accumulation of PABP1 and poly(A)RNA and thus facilitating eukaryotic cells to attenuate protein synthesis and energy consumption to adapt to energy stress. Moreover, AMPK-directed SIRT1 phosphorylation is required for energy starvation-induced PABP1-SIRT1 association, PABP1 deacetylation, and poly(A)RNA nuclear retention. In addition, the SIRT1-PABP1 association is not specific to energy starvation but represents a common stress response. These observations provide insights into dynamic modulation of eukaryotic mRNA transport and translation, suggesting that the poly(A)-tail also provides a basis for eukaryotes to effectively shut down mature mRNA transport and thereby tailor protein synthesis to maintain energy homeostasis under stress conditions.

Tissue Extract Fractions from Starfish Undergoing Regeneration Promote Wound Healing and Lower Jaw Blastema Regeneration of Zebrafish

Natural bioactive materials provide an excellent pool of molecules for regenerative therapy. In the present study, we amputate portions of the arms of Archaster typicus starfish, extract and separate the active biomaterials, and compare the effects of each fraction on in vitro wound healing and in vivo lower jaw regeneration of zebrafish. Compared with crude extract, normal hexane fractions (NHFs) have a remarkable effect on cellular proliferation and collective migration, and exhibit fibroblast-like morphology, while methanol-water fractions (MWFs) increase cell size, cell-cell adhesion, and cell death. Relative to moderate mitochondrialand lysosomal aggregation in NHFs-cultured cells, MWFs-cultured cells contain more and bigger lysosomal accumulations and clump detachment. The in vivo zebrafish lower jaw regeneration model reveals that NHFs enhance blastema formation and vasculogenesis, while MWFs inhibit fibrogenesis and induce cellular transformation. Gene expression analyses indicate that NHFs and MWFs separately activate blastema-characteristic genes as well as those genes-related to autophagy, proteasome, and apoptosis either during cell scratch healing or ganciclovir-induced apoptosis. Our results suggest that bioactive compounds from NHFs and MWFs could induce blastema formation and remodeling, respectively, and prevent tissue overgrowth.

Lessons on skeletal cell plasticity from studying jawbone regeneration in zebrafish

Three major mesenchymal cell types have important roles in determining the shapes of vertebrate animals: bone-producing osteoblasts, cartilage-producing chondrocytes, and fat-producing adipocytes. Although often considered discrete cell types, accumulating evidence is revealing mesenchymal cells of intermediate identities and interconversion of cell types. Such plasticity is particularly evident during adult skeletal repair. In this Review, we highlight recent work in zebrafish showing a role for hybrid cartilage-bone cells in large-scale regeneration of the adult jawbone, as well as their origins in the periosteum. An emerging theme is that the unique mechanical and signaling environment of the adult wound causes skeletal cell differentiation to diverge from the discrete lineages seen during development, which may aid in rapid and extensive regeneration of bone.

Ihha induces hybrid cartilage-bone cells during zebrafish jawbone regeneration

The healing of bone often involves a cartilage intermediate, yet how such cartilage is induced and utilized during repair is not fully understood. By studying a model of large-scale bone regeneration in the lower jaw of adult zebrafish, we show that chondrocytes are crucial for generating thick bone during repair. During jawbone regeneration, we find that chondrocytes co-express genes associated with osteoblast differentiation and produce extensive mineralization, which is in marked contrast to the behavior of chondrocytes during facial skeletal development. We also identify the likely source of repair chondrocytes as a population of Runx2(+)/Sp7(-) cells that emanate from the periosteum, a tissue that normally contributes only osteoblasts during homeostasis. Analysis of Indian hedgehog homolog a (ihha) mutants shows that the ability of periosteal cells to generate cartilage in response to injury depends on a repair-specific role of Ihha in the induction as opposed to the proliferation of chondrocytes. The large-scale regeneration of the zebrafish jawbone thus employs a cartilage differentiation program distinct from that seen during development, with the bone-forming potential of repair chondrocytes potentially due to their derivation from osteogenic cells in the periosteum.

Zebrafish Craniofacial Development: A Window into Early Patterning

The formation of the face and skull involves a complex series of developmental events mediated by cells derived from the neural crest, endoderm, mesoderm, and ectoderm. Although vertebrates boast an enormous diversity of adult facial morphologies, the fundamental signaling pathways and cellular events that sculpt the nascent craniofacial skeleton in the embryo have proven to be highly conserved from fish to man. The zebrafish Danio rerio, a small freshwater cyprinid fish from eastern India, has served as a popular model of craniofacial development since the 1990s. Unique strengths of the zebrafish model include a simplified skeleton during larval stages, access to rapidly developing embryos for live imaging, and amenability to transgenesis and complex genetics. In this chapter, we describe the anatomy of the zebrafish craniofacial skeleton; its applications as models for the mammalian jaw, middle ear, palate, and cranial sutures; the superior imaging technology available in fish that has provided unprecedented insights into the dynamics of facial morphogenesis; the use of the zebrafish to decipher the genetic underpinnings of craniofacial biology; and finally a glimpse into the most promising future applications of zebrafish craniofacial research.

Sox9-related signaling controls zebrafish juvenile ovary-testis transformation

In almost all vertebrates, the downstream of the sox9 signaling axis is well conserved for testis differentiation. The upstream genes of this pathway vary from species to species during evolution. Yet, little is known about how these signaling cascades are regulated and what cellular processes are dominant in ovary–testis transformation in juvenile zebrafish. In this study, we find that the transforming gonads undergo activation of sox9a-expressing stromal cells with increased deposition of extracellular matrix and formation of degenerative compartments. This leads to follicle disassembly, oocyte degeneration, follicle cell-cyp19a1a-amh conversions, and, eventually, formation of the testis cord. In vitro primary culture of juvenile ovary tissue in gonadotropins increases cytoplasmic accumulation of sox9a and p-Erk1/2, and induces mesenchymal morphology. MAPK inhibitors (MKI), a mixture of PD98059 and U0216, eliminate the cytoplasmic distribution but do not eradicate nuclear localization of sox9a and p-Erk1/2. Nuclear p53 are relatively increased in MKI-treated cells that exhibit less spreading and reduced proliferation. Despite uniform nuclear condensation, only a fraction of cells displayed the apoptotic phenotype. These results suggest that high levels of cytoplasmic sox9a and p-Erk1/2 activity activate stromal cells and enhance the production of extracellular matrix required for testis cord formation, whereas deregulation and translocation of sox9a and p-Erk1/2 induce follicle disassembly and incomplete apoptosis associated with nuclear p53. Together with the established FSH/cAMP/MAPK/AMH pathway in mammalian granulosa and Sertoli cells, we demonstrated that the sox9 axis signaling that determines testis formation in mammals also induces zebrafish ovary–testis transition, and adds to its conserved role in sex reversal.

The zebrafish as a model for complex tissue regeneration

For centuries, philosophers and scientists have been fascinated by the principles and implications of regeneration in lower vertebrate species. Two features have made zebrafish an informative model system for determining mechanisms of regenerative events. First, they are highly regenerative, able to regrow amputated fins, as well as a lesioned brain, retina, spinal cord, heart, and other tissues. Second, they are amenable to both forward and reverse genetic approaches, with a research toolset regularly updated by an expanding community of zebrafish researchers. Zebrafish studies have helped identify new mechanistic underpinnings of regeneration in multiple tissues and, in some cases, have served as a guide for contemplating regenerative strategies in mammals. Here, we review the recent history of zebrafish as a genetic model system for understanding how and why tissue regeneration occurs.

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

The salamander has the remarkable ability to regenerate its limb after amputation. Cells at the site of amputation form a blastema and then proliferate and differentiate to regrow the limb. To better understand this process, we performed deep RNA sequencing of the blastema over a time course in the axolotl, a species whose genome has not been sequenced. Using a novel comparative approach to analyzing RNA-seq data, we characterized the transcriptional dynamics of the regenerating axolotl limb with respect to the human gene set. This approach involved de novo assembly of axolotl transcripts, RNA-seq transcript quantification without a reference genome, and transformation of abundances from axolotl contigs to human genes. We found a prominent burst in oncogene expression during the first day and blastemal/limb bud genes peaking at 7 to 14 days. In addition, we found that limb patterning genes, SALL genes, and genes involved in angiogenesis, wound healing, defense/immunity, and bone development are enriched during blastema formation and development. Finally, we identified a category of genes with no prior literature support for limb regeneration that are candidates for further evaluation based on their expression pattern during the regenerative process.