Revealing the biological features of the axolotl pancreas as a new research model

Introduction: The pancreas plays a crucial role in digestion and blood glucose regulation. Current animal models, primarily mice and zebrafish, have limited the exploration of pancreatic biology from an evolutionary-developmental perspective. Tetrapod vertebrate axolotl (Ambystoma mexicanum) serves as a valuable model in developmental, regenerative, and evolutionary biology. However, the fundamental biology of the axolotl pancreas remains underexplored. This study aims to characterize the unique developmental, functional, and evolutionary features of the axolotl pancreas to expand the understanding of pancreatic biology.

Methods: We conducted morphological, histological, and transcriptomic analyses to investigate the axolotl pancreas. Pancreatic development was observed using in situ hybridization and immunostaining for key pancreatic markers. RNA sequencing was performed to profile global gene expression during larva and adult stages. And differential gene expression analysis was used to characterize the conserved and unique gene patterns in the axolotl pancreas. Functional assays, including glucose tolerance tests and insulin tolerance tests, were optimized for individual axolotls. To assess pancreatic gene function, Pdx1 mutants were generated using CRISPR/Cas9-mediated gene editing, and their effects on pancreatic morphology, endocrine cell populations, and glucose homeostasis were analyzed.

Results: The axolotl pancreas contains all known pancreatic cell types and develops from dorsal and ventral buds. Both of buds contribute to exocrine and endocrine glands. The dorsal bud produces the major endocrine cell types, while the ventral bud generates α and δ cells, but not β cells. Differential gene expression analysis indicated a transition in global gene expression from pancreatic cell fate commitment and the cell cycle to glucose response, hormone synthesis, and secretion, following the development progression. Notably, the adult axolotl pancreas exhibits slower metabolic activity compared to mammals, as evidenced by the results of GTT and ITT. The mutation of Pdx1 resulted in hyperglycemia and a significant reduction in pancreatic cell mass, including a complete loss of endocrine cells, although it did not lead to a lethal phenotype.

Discussion: This study examines the axolotl pancreas, highlighting the conservation of pancreatic development. Our study highlights the unique features of the axolotl pancreas and broadens the scope of animal models available for pancreatic evolution and disease research.

1 Introduction

The pancreas is an important organ in the mammalian digestive system, composed of exocrine and endocrine cells that release enzymes for digestion and hormones to regulate blood glucose metabolism. Clustering of the endocrine cells forms the islets of Langerhans in mammals, which consist of five different endocrine cell types, Insulin (INS)-secreting β cells, Glucagon (GCG)-producing α cells, Somatostatin (SST)-secreting δ cells, Pancreatic polypeptides (PPY)-producing PP cells and ε cells that secrete Ghrelin (GHRL) (Lewis and Mao, 2023; Steiner et al., 2010; Bonner-Weir et al., 2015). These endocrine cells regulate each other’s secretion and work synergistically to maintain blood glucose level (Röder et al., 2016). Diabetes mellitus, characterized by hyperglycaemia due to absolute or relative insulin insufficiency arising from β cell impairments in human (Cole and Florez, 2020; Ghila et al., 2022), is one of the most significant challenges in public health (Hossain et al., 2024). Investigating the mechanisms of pancreatic development and regeneration using model organisms can enhance our understanding of how endocrine cells are generated and regenerated, providing insights into the fundamental principles of organ development and potential regenerative therapies for pancreatic diseases.

Currently, animal models such as mice and zebrafish are the most commonly used organisms for studying the development and regeneration of pancreatic β cells (Aguayo-Mazzucato and Bonner-Weir, 2018; Yang et al., 2020). The pancreas of rodents and zebrafish exhibits significant similarity in cell type composition and gene expression programs during development and homeostasis (Slack, 1995). However, they also exhibit certain variations that can confer advantage or disadvantage for pancreas research (Mi et al., 2024; Weng et al., 2019). While rodents are evolutionarily closer to humans, resulting in similarities in pancreatic structure and the state of exocrine and endocrine cells, they present limitations in the context of pancreatic regeneration. Notably, adult mice have a relatively low β cell replication rate and are unable to restore β cell mass under type I and type II diabetic conditions (Rankin and Kushner, 2009; Sasaki et al., 2021; Teta et al., 2005). In contrast, zebrafish exhibit a higher capacity for β cell regeneration, providing advantages over mammalian models for studies on pancreatic regeneration (Mi et al., 2024; Moss et al., 2009). Moreover, the transparency and small size of zebrafish make them ideal for live imaging of the pancreas and drug screening, including the identification of compounds that promote β cell regeneration using genetically modified strains (Li et al., 2021; Teame et al., 2019). However, the small body size of zebrafish presents challenges for examining physiological parameters of glucose metabolism, such as glucose tolerance test (GTT) and insulin tolerance test (ITT) (Kimmel and Meyer, 2016; Yang et al., 2020). Consequently, exploring new animal models positioned at crucial evolutionary nodes and exhibiting unique biological features will enrich the resources available for pancreatic research and provide additional insights into pancreatic development, regeneration, and evolution.

The axolotl, a tetrapod vertebrate, possesses remarkable regenerative capacities (Vieira et al., 2020). Previous studies have demonstrated that many organs in axolotls, including limbs, tails, brains, spinal cords, hearts and lenses, can regenerate following injury (Bryant et al., 2017; Cano-Martínez et al., 2010; Cura Costa et al., 2021; Echeverri and Tanaka, 2002; Lust et al., 2022; Suetsugu-Maki et al., 2012; Wei et al., 2022). From an evolutionary perspective, axolotls belong to amphibians and are classified between fish and rodents (Lu, 2023). Despite extensive research on canonical organs in axolotls, the pancreas has been poorly characterized, limiting its application in pancreatic biology. Recently, Sørensen and colleagues identified the presence of islets in the axolotl pancreas using Hematoxylin-Eosin (H&E) staining and immunohistochemistry with an insulin antibody on tissue sections. They also demonstrated streptozotocin-induced hyperglycemia, detected via intraperitoneal GTT (Sørensen et al., 2022). However, little is known about the complete cellular composition, developmental features, gene expression, and genetic regulation of pancreatic function in axolotls as a potential model for diabetes mellitus.

In this study, we reported that the axolotl pancreas comprises ducts, acinar cells and islets, which contain α, β, δ, and ε cells. Notably, we observed that the PP cells, which secrete PPY, co-expressed the typical α cell marker GCG. Developmentally, we found that a significant portion of the pancreas originates from the dorsal bud, while the uncinate process, derived from the ventral bud, contains α and δ cells but no β cells. At the gene expression level, we noted strong conservation of gene regulation in axolotl pancreas during development and homeostasis, although it exhibited slower metabolic activity compared to mammals. Moreover, following the mutation of Pancreatic and duodenal homeobox 1 (Pdx1), the mutant axolotls survived but exhibited a hyperglycemic phenotype, deficiencies in the exocrine pancreas, and varying degrees of loss of endocrine cells. Additionally, intraperitoneal GTT and ITT revealed a delayed response to glucose and insulin in the axolotl pancreas, consistent with our gene expression analysis and suggesting a slower rate of glucose metabolism. Overall, our study characterizes the fundamental biological features of the axolotl pancreas and lays the groundwork for using axolotl as a model for pancreatic research.

2 Materials and methods

2.1 Animal husbandry

The axolotl (Ambystoma mexicanum) d/d strain used in this study was maintained in fresh water, with individuals housed separately and fed daily. Prior to surgical procedures, axolotls were anaesthetized using a 0.01% solution of ethyl-p-aminobenzoate (Benzocaine; Sigma). All animal experiments were carried out according to the guidelines of the ethics committee of Guangdong Provincial People’s Hospital (approval number KY-Q-2022-395-02).

2.2 Generation of Pdx1 mutant axolotls

Pdx1 mutant axolotls were generated using CRISPR/Cas9 technology, following established protocols. The Pdx1 guide RNA (gRNA) sequence (5′-AAG​GAG​GAG​GAC​AAG​AAG​CG-3′) was synthesized by GenScript. A mixture containing 0.4 μg/μL Pdx1 gRNA, 0.5 μg/μL Cas9 protein, and 1× Cas9 buffer was injected into single-cell stage embryos from the d/d strain, as previously described (Fei et al., 2018). The injected embryos developed into juveniles, which were subsequently genotyped to identify the F₀ transgenic animals. Genomic DNA was extracted from axolotl tail snips using a solution containing 50 mM NaOH and 1 M Tris-HCl. The target locus of the Pdx1 gRNA was PCR amplified using primers Pdx1-KOfw1 (5′-act​tga​cga​gct​gtt​cgc​cca​ttt​c-3′) and Pdx1-KOre2 (5′-tga​gca​tgt​cct​tgc​cag​ggt​ccc-3′). The PCR products were then subjected to Sanger sequencing (Tsingke) for alignment with the original Pdx1 sequence. F₀ transgenic axolotls were further bred to generate F₁ offspring with various genotypes (Supplementary Figure S1).

2.3 Immunofluorescence

Samples were fixed in 3.7% formaldehyde freshly prepared in MEM buffer (MEMFA) overnight at 4°C. After dehydration in 30% sucrose, the samples were embedded in OCT (Sakura, 4,583). Cryosections of 10 µm were then collected for immunohistochemical analysis as previously described (Fei et al., 2017). The sections were blocked with a blocking buffer (0.3% Triton X-100, 5% fetal calf serum in PBS) for 1 h at room temperature, then incubated with primary antibodies overnight at 4°C. The primary antibodies used in this study were: Rabbit anti-Somatostatin antibody (Abclonal, A9274, 1:200), Guinea pig anti-Glucagon antibody (TAKARA, M182, 1:1,000), Mouse anti-Insulin antibody (Sino Biological, 11038-MM14-4.5, 1:1,000). Following several washes in PBST (0.3% Triton X-100 in PBS), the corresponding secondary antibodies Alexa Fluor 488-Goat anti-guinea pig (Invitrogen, A11073, 1:500), Alexa Fluor 555-Donkey anti-Rabbit (Thermofisher, A-31572, 1:500), Alexa Fluor 647-Donkey anti-Mouse (Jackson, 715605150, 1:500) and DAPI (Sigma, 28718903,1:000), diluted in blocking buffer were applied to the sections and incubated for 2 h at room temperature. The slides were then washed several times in PBST and mounted with Mowiol medium. Immunofluorescence images were acquired using a confocal microscope (ZEISS, LSM980).

2.4 H&E staining

Pancreatic tissue from adult axolotls were collected and fixed in MEMFA overnight at 4°C. After dehydration using a graded series of ethanol (50%, 70%, 95% and 100% diluted in PBS) and clearing in xylene, the samples were embedded in paraffin, and 5 μm sections were collected. For H&E staining, the sections were dried at 65°C for 1 h, treated with xylene to remove the paraffin, rehydrated in decreasing concentrations of ethanol, and finally in distilled water. The sections were incubated in hematoxylin solution for 5 min. After several washes with distilled water, the sections were sequentially incubated in 1% hydrochloric acid in ethanol for 10 s, washed gently with running water for 15 min, and then incubated in eosin solution for 1 min. The stained sections were dehydrated with ethanol and transferred into xylene, then mounted with neutral resin. Once the resin dried, the sections were imaged using an inverted microscope (Olympus, IX83).

2.5 In situ hybridization

To synthesize DIG-labelled or fluoresce-labelled antisense RNA probes, corresponding DNA templates containing the T7 promoter were amplified with two rounds of PCR, using axolotl pancreatic cDNA synthesized from total RNA of the pancreas as the template. The antisense RNA probes were transcribed in vitro using a T7 RNA polymerase transcription kit (NEB, M0251L). In situ hybridizations and fluorescence in situ hybridization were performed on 10 μm cryosections from axolotl pancreas or embryos, as previously described (Fei et al., 2017; He et al., 2020). Briefly, for in situ hybridizations, 500 ng/mL of DIG-labeled probes were hybridized at 60°C on 10 μm cryosections overnight. After hybridization, the cryosections were washed with a stringent buffer and incubated overnight at 4°C with anti-DIG-AP antibody (1:5,000, Roche, 11093274910). The signal was detected using BM Purple solution (Roche, 11442074001). For fluorescence in situ hybridization, a mixture of 300 ng/mL DIG-labeled and fluorescein-labeled RNA probes was hybridized at 60°C overnight on 10 μm cryosections. After washing with stringent buffer, the cryosections were incubated with anti-fluorescein-POD (1:500, Roche, 11426346910), and then the signal was detected using TSA fluorescein (1:50, PerkinElmer, NEL701A001KT) at room temperature. Following signal development, 3% hydrogen peroxide was used to deactivate any residual POD, and the cryosections were incubated with anti-DIG-POD (1:1,000, Roche, 11633716001) after washing. TSA Cy5 (1:50, PerkinElmer, NEL705A001KT) was used to generate a detectable signal. List of primers for RNA probe synthesis: Amy-forward: 5′-ATG​AGG​CTC​CTT​CTG​CTG​CTT​GC-3′, Amy-reverse: 5′-CTT​GTA​TTG​TTT​ATA​ACT​TG-3′, T7-Amy-reverse: 5′-TTG​AAA​TTA​ATA​CGA​CTC​ACT​ATA​GGG​CTT​GTA​TTG​TTT​ATA​ACT​TG-3′, Ins-forward: 5′-ATG​GCT​CTT​TGG​GTG​CGT​G-3′, Ins-reverse: 5′ACT​CCT​TCT​ACC​CCC​TCT​GC-3′, T7-Ins-reverse: 5′-TTG​AAA​TTA​ATA​CGA​CTC​ACT​ATA​GGG​ACT​CCT​TCT​ACC​CCC​TCT​GC-3′, Krt19-forward: 5′-ATG​TCG​TCG​AGT​ACA​GCT​TCA​GG-3′, Krt19-reverse: 5′-TTT​TCC​CCC​GCT​TTT​GGA​C-3′, T7-Krt19-reverse: 5′-ATT​ATG​CTG​AGT​GAT​ATC​CCT​CTT​TTT​CCC​CCG​CTT​TTG​GAC-3′, Ptf1a-forward: 5′-TGC​AGT​CCA​TCA​ACG​ACG​CCT​TC-3′, Ptf1a-reverse: 5′-GCG​TCA​CTG​TAA​TGA​AAG​C-3′, T7-Ptf1a-reverse: 5′-TTG​AAA​TTA​ATA​CGA​CTC​ACT​ATA​GGG​GCG​TCA​CTG​TAA​TGA​AAG​C-3′, Mln-forward: 5′-CCC​ACA​TCC​ACA​CTC​TAC​AGT​AG-3′, Mln-reverse: 5′-TCA​GTC​CAC​ATT​TTG​GTG​TTC​TG-3′, T7-Mln-reverse: 5′-TTG​AAA​TTA​ATA​CGA​CTC​ACT​ATA​GGG​TCA​GTC​CAC​ATT​TTG​GTG​TTC​TG-3′.

2.6 Whole-mount in situ hybridization

Axolotl embryos were fixed overnight in MEMFA at 4°C. Following fixation, the embryos were washed with PTW (DEPC-treated 0.1× PBS with 0.2% Tween) through a gradient of decreasing methanol concentrations (100%, 90%, 75%, 50%, and 25%). The samples were then treated with proteinase K (5 μg/mL) at room temperature for 5 min, followed by incubation at 65°C in hybridization buffer with 500 ng/mL antisense RNA probes. To remove unbound RNA probes, the samples were washed multiple times with post-hybridization buffer and MAB buffer. Next, samples were incubated with anti-DIG-AP antibody (1:5,000, Roche, 11093274910) in 1× blocking buffer (Roche, 11096176001). The signal was detected using BM Purple substrate (Roche, 11442074001). Finally, the samples were stored in 50% glycerol diluted in PBS at 4°C and imaged using an Olympus SZX10 microscope.

2.7 Transmission electron microscopy

TEM was conducted following standard procedures. Briefly, the isolated axolotl pancreas was cut into approximately 1 mm³ tissue blocks, which were then transferred to an EP tube containing fresh fixative for TEM (Servicebio; G1102) and fixed overnight at 4°C. After fixation, the samples were rinsed three times with 0.1M phosphate buffer (pH7.4) and post-fixed in 1.0% osmium tetroxide in phosphate buffer at room temperature for 2 h. The tissue was then dehydrated through a graded series of ethanol, embedded in Epon resin, and polymerized at 60°C overnight. Ultrathin sections were cut, mounted on copper grids, and stained with uranyl acetate and lead citrate. All grids were observed and imaged using a transmission electron microscope (Hitachi; HT7800/HT7700).

2.8 RNA sequencing

Total RNA was extracted from the pancreas of larvae (about 3–4 cm, 2 months old) and adult axolotls (about 25 cm, 2 years old) using TRIzol reagent (Invitrogen) and purified via LiCl precipitation. RNA sequencing was conducted on six biological samples, comprising three replicates from both the larval and adult groups. Specifically, three adults and twenty-four larvae were used for sequencing sample preparation. Three micrograms of total RNA were utilized to construct cDNA libraries according to the standard Illumina RNA-seq protocol. The generated cDNA library was sequenced using the Illumina NovaSeq6000 via a 2 × 150 bp paired-end protocol. Quality filtering and removal of residual adaptor sequences was conducted on read pairs using Trimmomatic v0.33. Salmon v1.10.0 was employed to build the axolotl genome index and guide transcript quantification. Axolotl genome assembly v6.0-DD from the axolotl-omics database (https://www.axolotl-omics.org/assemblies) (Schloissnig et al., 2021) served as the reference genome along with annotation files.

2.9 RNA sequencing analysis

Transcript expression levels were imported into R v4.3 and summarized at the gene level using the R/tximport v1.30.0. Gene count data were analyzed for differential expression using the Bioconductor software package DESeq2 v1.42.1. Gene expression was quantified using the transcript per million (TPM) values. Metascape was utilized for gene functional annotation. Mfuzz was applied to analyze the clustering of time series gene expression in larval and adult axolotl pancreas. For conservation analysis, RNA sequence data from untreated human samples (GSE205853) (Benaglio et al., 2022), Wt_2w and Wt_12w group of mice (GSE171251) (de Jesus et al., 2021) were compared to our adult axolotl RNA sequence data. For the identification of homologous genes, we obtained gene annotations from two sources: the human and mouse genome annotations from Ensembl, and the axolotl genome annotation (version 6.0-DD) from the axolotl-omics database. Genes that share annotations with mouse and human were considered homologous genes. Gene length and sequencing depth were standardized by using TPM values to compare homologous gene expression across the three species. |log2Fold Change| ≥ 1 with a False Discovery Rate (FDR) < 0.05 was utilized to capture meaningful changes in gene expression during development for the axolotl and mouse. To identify higher and lower expressed genes in axolotls compared to mice and humans, |log2Fold Change| ≥ 2 and FDR <0.01 were applied. Principal component analysis (PCA) was performed using the prcomp function on TPM-transformed raw counts of the most highly expressed genes across all pancreas samples. Pearson correlation coefficients were calculated to assess the correlations between the three species.

2.10 Glucose tolerance test and insulin tolerance test

Animals measuring about 15 cm in length were used for GTT and ITT. For the GTT experiment, axolotls were fasted for 12–16 h before being intraperitoneally injected with D-Glucose (2 mg/g) after anesthesia. Blood glucose levels were monitored at specific intervals using a glucose meter (Roche, Accu-Chek), sampling from the punctured dorsal body surface vessel. 48 h after the GTT experiment, the animals were returned to clear water for normal feeding. For the ITT experiment, axolotls were fasted for 6 h before being intraperitoneally injected with recombinant human insulin (insulin glargine; 5 μL/g) after anesthesia. Blood glucose levels were monitored at specific intervals as outlined in the GTT procedure. Similarly, following the ITT experiment, the animals were injected with 2 mg/g glucose solution and returned to clear water for normal feeding; the same concentration of glucose solution was injected for the next 2 days.

2.11 Statistical analysis

All axolotls were randomly assigned to different groups, and all experiments were repeated at least three times independently. For experimental data analysis, calculations were performed using GraphPad Prism 8.0.2. Unpaired, two-tailed Student’s t-tests were performed to analyze the P values for single comparisons between two groups, while one-way ANOVA was conducted to analyze P values for multiple comparisons. All data were presented as the mean ± s. e.m. Each experiment was performed blind, and no predetermination was made for sample size. P value < 0.05 was considered significant. *P < 0.05; **P < 0.01; ***P < 0.001; ****P < 0.0001; ns, P > 0.05.

3 Results

3.1 Axolotl pancreas exhibits an evolutionarily conserved structure and cell type composition

To systematically characterize the axolotl pancreas, we isolated pancreas from adult axolotls and examined its morphological features, composition and distribution of pancreatic cell types. Consistent with other vertebrates, the axolotl pancreas is situated beneath the stomach. The pancreas head and tail connect to the duodenum and spleen, respectively, while the body of the pancreas is attached to the jejunum via the mesentery (Figures 1A, B). The uncinate process of the pancreas wraps around the duodenum and connects to the liver and gall bladder (Figure 1C). H&E staining of paraffin sections of the pancreas showed that islets are distributed throughout the pancreatic tissue, surrounded by exocrine tissues, including acinar cells, ducts and blood vessels (Supplementary Figure S2B; Figure 1D). In situ hybridization using classical pancreatic cell markers (Muraro et al., 2016) demonstrated the presence of Ins, a marker for endocrine β cells in islets, and Amylase (Amy) in acinar cells along with Keratin 19 (Krt19) in ductal cells (Figures 1E, F) in the axolotl pancreas.

Figure 1

Figure 1. The structure of adult axolotl pancreas. (A) The pancreas model of an adult axolotl shows that the basic structure of the pancreas includes the head, body, tail, and uncinate process. (B) Anatomical structure of the ventral view of the pancreas. The pancreas head connects to the duodenum, the tail connects to the spleen, and the body is suspended from the jejunum through the mesentery. L, liver; D, duodenum; J, jejunum; P, pancreas; S, spleen; G, gallbladder. (C) Anatomy of the dorsal view of the pancreas. After turning the jejunum to the right, the uncinate process connects to the liver and gallbladder (blue line). (D) Hematoxylin-Eosin staining of the axolotl pancreas. A, acinar; I, islet; D, ductal; E, erythrocyte; CA, centroacinar. (E–J) In situ hybridization and immunofluorescence staining show the molecular features of the axolotl pancreas. (E) Fluorescence in situ hybridization for Amylase (Amy, green), and Insulin (Ins, red), combined with DAPI (blue) on pancreas paraffin sections. (F) In situ hybridization for Keratin 19 (Krt19, purple) on pancreas cryosections. (G) Immunofluorescence for INS (red), Glucagon (GCG, green), and Somatostatin (SST, white), combined with DAPI (blue) on pancreas cryosections, showing the islets are scattered throughout the pancreas, and the uncinate process has no INS⁺ β cells. (H) Immunofluorescence of a magnified islet shows that the islet consists of a cluster of INS⁺ β cells (red), surrounded by GCG⁺ α cells (green) and SST⁺ δ cells (white). (I) Fluorescence in situ hybridization for Motilin (Mln, white) and Ins (red) on pancreas cross-cryosections shows a very small number of Mln⁺ ε cells surrounding the β-cell cluster. (J) Immunofluorescence for Pancreatic polypeptide (PPY, white), GCG (green), and INS (red) combined with DAPI (blue) on pancreas cryosections shows PPY⁺ PP cells surround β cells and co-express with GCG⁺ α cells. Scale bar: (B, C) 200 μm; (D–F, H–J) 100 μm; (G) 500 μm.

We next focused on characterizing the endocrine cells. It is well-known that pancreas contains five types of endocrine cells in mammals (Baron et al., 2016). To this end, we used specific endocrine cell markers INS, GCG, and SST to identify β, α, and δ cells respectively. Immunofluorescence on pancreatic sections showed that the main body of the pancreas (head, body, and tail) consists predominantly of INS⁺ β cells-clustered islets, surrounded by sparsely distributed GCG⁺ α cells and closely associated SST⁺ δ cells (Figures 1G, H). In contrast, the uncinate process of the pancreas exclusively contained GCG⁺ α and SST⁺ δ cells, with no INS⁺ β cells (Supplementary Figure S2A; Figure 1G). Additionally, double fluorescence in situ hybridization revealed that pancreatic islets composed of Ins⁺ β cells were tightly surrounded by a limited number of Motilin (Mln)⁺ ε cells (Figure 1I). Interestingly, immunofluorescence analysis using an antibody against PPY revealed that PPY is expressed exclusively in GCG⁺ α cells (Figure 1J), suggesting the identities of α and PP cells are perhaps intermingled in axolotls. Electron microscopy showed acinar cells containing large electron-dense zymogen granules likely enriched with digestive proteins, while β cells were filled with numerous insulin granules that lacked halos, in contrast to their mammalian counterparts (Roep, 2020; Susan, 2004) (Supplementary Figures S2C, G). In addition, we also found the presence of α and δ cells in the stomach and small intestine of embryonic and juvenile axolotls (Supplementary Figure S3). In summary, axolotl pancreas, although exhibiting some unique feature, demonstrates significant structural, compositional, and subcellular similarities to other species.

3.2 Axolotl pancreas develops from ventral and dorsal buds

We next investigated the developmental origin of the adult pancreas. The differences in endocrine cell types between the uncinate process and the main portion of the axolotl pancreas suggest distinct embryonic developmental mechanism. During mammalian pancreas development, it has been reported that the ventral and dorsal pancreatic anlage give rise to the uncinate process and the main pancreas, respectively (Longnecker, 2014). To determine whether the axolotl pancreas develops from ventral and dorsal buds during the embryonic stage, we performed whole-mount and cryosection in situ hybridization on axolotl embryos. We detected the pancreatic ventral and dorsal buds as early as stage 36, as indicated by the expression of Pancreas associated transcription factor 1a (Ptf1a) and Amy (Figures 2A, D, E). By stage 43, the dorsal bud expanded and elongated in conjunction with the growth of the intestinal tube (Figures 2B, F, G). Eventually, the dorsal and ventral buds fused to form early pancreas a few days following stage 43 (Figures 2C, H, I). Consistent with the cell types identified in adult pancreas, immunohistochemistry analysis revealed that the β cells, δ cells and α cells in the dorsal bud started to develop at stages 40, 42 and 43, respectively (Figure 2J), while only α and δ cells were observed in the ventral bud (Figure 2K). These results indicate that the major body of the pancreas originates from the dorsal bud, whereas the uncinate process arises from the ventral bud, aligning with mammalian pancreas development.

Figure 2

Figure 2. The development of the pancreas in axolotl embryos. (A–C) A model of pancreatic development. At stage 36, the dorsal bud and ventral bud merge. At stage 43, the dorsal bud elongates along intestinal tube. After 5 days, the dorsal bud and ventral bud fuse at stage 44. The arrows indicate the slice direction. The light blue rectangles represent sections shown in panels (D–I). The red section indicates the dorsal bud, and the pink section indicates the ventral bud. (D) In situ hybridization for Pancreas associated transcription factor 1a (Ptf1a, purple) on paraffin sections of embryos shows that Ptf1a is expressed in both the dorsal bud (red circle) and the ventral bud (pink circle) at stage 36. (E) Whole-mount in situ hybridization shows that Amy is expressed in both the dorsal (blue) and ventral (purple) buds at stage 39. (F, G) In situ hybridization for Amy (purple) shows that the dorsal bud extends as the intestinal tube grows. (H, I) In situ hybridization for Amy (purple) shows that the dorsal and ventral buds have fused by stage 44. (J) Immunofluorescence for INS (red), GCG (green), SST (white), combined with DAPI (blue) on embryo cryosections shows that endocrine cells begin to emerge from dorsal bud. INS⁺ β cells (green) emerges at stage 40, and SST⁺ δ cells (white) and GCG⁺ α cells (red) emerging at stages 42 and 43, respectively. (K) Immunofluorescence shows that INS⁺ β cells are concentrated in the dorsal bud of the pancreas, while only GCG⁺ α cells (red) and SST⁺ δ cells (white) are found in the ventral bud. Eb, external branchia; St, stomach; Int, intestine; L, live; Lu, lung; H, heart. Scale bar: (D–I) 200 μm; (J) 50 μm; (K) 100 μm.

3.3 The transcriptome of the axolotl pancreas exhibits evolutionary conservation and distinction

To explore the molecular characteristics of both the developing and adult pancreas, we sequenced and analyzed the gene expression profiles of total RNA extracted from 2-month-old larvae and 2-year-old adult axolotls. Our analysis identified a total of 23,931 genes expressed in the developing pancreas and 19,592 genes in the adult pancreas. Of these, 18,492 genes were commonly expressed at both stages (Supplementary Figure S4A). Gene Ontology analysis of the genes uniquely expressed in the developing pancreas (5,439 genes, with 2,634 annotated genes) revealed their primarily involvement in tissue morphogenesis and digestive system development. In contrast, the genes expressed in the adult pancreas (1,100 genes, with 565 annotated genes) were mainly associated with the positive regulation of the phosphorus metabolic process (Supplementary Figure S4B).

Further analysis on the 18,492 (16,614 annotated genes) commonly expressed genes identified 1,888 upregulated genes (1,238 annotated genes) and 2,016 downregulated genes (1,210 annotated genes) in the adult pancreas compared to the developing pancreas (Supplementary Figure S4C). Gene Ontology analysis indicated that the upregulated genes were primarily linked to hormone responses, secretion regulation, and hormone level control, while the downregulated genes were typically associated with the mitotic cell cycle, DNA metabolic processes, and cell cycle regulation (Figure 3A). Moreover, clustering of time series gene expression shows that genes involved in ribonucleoprotein complex biogenesis and RNA splicing are also highly active, in addition to those related to mitotic cell cycle phase transitions in larvae. Genes associated with signal release, peptide transport, and insulin secretion are even more active in adults (Supplementary Figure S5). These findings suggest a transition of the axolotl pancreas from an actively developing state in early stages to a homeostatic state in adulthood, indicating functional maturation during development, similar to observations in mice (de Jesus et al., 2021) (Supplementary Figures S4D, E).

Figure 3

Figure 3. Differential gene expression analysis in the larval and adult pancreas of axolotls. (A) Gene Ontology (GO) analysis of the top five downregulated and upregulated genes during development shows that GO terms related to cell division, such as the regulation of mitosis, DNA metabolism, and the cell cycle, are downregulated, while GO terms related to digestion, such as response to hormone, regulation of secretion, and regulation of hormone levels, are upregulated. (B) Differential gene expression (DGE) analysis of the adult and larva pancreas after normalization. (C) Gene homology assessment of human, adult axolotls, and 12-week-old mice reveals 6,457 commonly expressed genes, which are annotated in all three species. (D) GO analysis of the 6,457 commonly expressed genes shows that, compared to mice and humans, axolotls have 131 genes with higher expression and 126 genes with lower expression. (E) Comparison of the expression levels of key pancreatic genes in adult axolotls, mice, and humans after normalization. (F) GO analysis of genes with low expression in axolotls reveals their involvement in processes such as insulin secretion, glucose homeostasis, endocrine development, and organ regeneration.

The development and functional maintenance of the pancreas require coordinated gene expression regulation. We focused on well-known genes that play crucial roles during the development of pancreatic lineages (Aigha and Abdelalim, 2020; Burlison et al., 2008; Yu et al., 2021). The expression levels of markers for most pancreatic multipotent progenitor cells (Pdx1, Ptf1a, Nkx6.1, Nkx2.2 and Foxa2), exocrine progenitor cells (Ptf1a and Sox9), and endocrine progenitor cells (Neurog3, Isl1, Neurod1 and Pax6) were significantly higher in larvae compared to adult axolotls (Figure 3B). Notably, α cell fate determiners Arx, MafB, and Peg10 exhibited comparable expression levels between larvae and adults. In contrast, β cell fate-related genes Gng12 and MafA were expressed at higher levels in larvae. Functional hormone-coding genes in the pancreas, including Ins and Amy, were expressed significantly higher in adults than in larvae, while Gcg, Sst, Ppy and Mln showed lower expression levels. This regulation pattern correlates with increased demands for food digestion and metabolic homeostasis during growth. Overall, the expression patterns and dynamics of key pancreatic genes during development in juvenile and adult axolotls are conserved with those observed in mammals (Yu and Xu, 2020).

To explore the evolutionarily distinctive features of the axolotl pancreas, we compared the transcriptome of adult axolotls with those of humans and mice (Benaglio et al., 2022; de Jesus et al., 2021). Among the 19,592 genes in the adult pancreas, 11,860 were found to have homologous counterparts in both humans and mice. Among these, 6,457 genes were commonly expressed across all three species (Supplementary Figure S4F; Figure 3C). Key pancreatic genes were present in all three species, except for Neurog3, which is absent in both humans and axolotls (Figure 3E). We identified and analyzed differentially expressed genes exhibiting at least four-fold difference in expression levels between axolotls and mice or humans. Our results revealed that 131 genes were expressed at higher levels and 126 genes at lower levels in adult axolotls compared to the other two species (Supplementary Figure S4G). The highly expressed genes were primarily involved in cytoplasmic translation, while the lowly expressed genes were related to hormone transport and negative regulation of cell population proliferation (Figure 3D). Among the 126 downregulated genes, such as Pdx1, Neurod1, Ptprn, Slc30a8, Rbp4, and Sox4, many were associated with insulin synthesis, insulin secretion, and glucose response, indicating that axolotl pancreas may display a delayed insulin secretion response to glucose compared to mammals (Figure 3F). Therefore, we examined the expression of metabolic genes in the axolotl and found that glucose sensor genes Slc2a1 and Slc2a2, insulin secretion genes Stx1a and Vamp8, and glycolytic-related kinases Pfkl, Pfkm, Pkm, Gck, Pgk1, and Ldhb are lowly expressed in the axolotl (Supplementary Figure S4H). Additionally, the genes Slc7a5, Pdx1, Cdkn1a, and Cebpb are known to be downregulated in association with positive effects on organ regeneration (Fitzner et al., 2013; Holland et al., 2005; Juran et al., 2021; Wang et al., 2022). We observed that these genes were expressed at lower levels in axolotls (Figure 3F). In conclusion, the transcriptomic characteristics of the axolotl pancreas indicate an evolutionary conservation of gene expression, suggesting that glucose regulation within pancreatic endocrine cells is less sensitive compared to mammals. Moreover, this may point to a significant regenerative potential within the axolotl pancreas.

3.4 Mutations of Pdx1 cause various levels of pancreas defects but not lethality

Evaluating gene function is crucial for understanding biological processes. Therefore, we targeted Pdx1 to assess the functional conservation and diversity of pancreatic-specific genes. Pdx1 is a transcription factor expressed in multipotent progenitor cells; it is essential for pancreatic development and the maintenance of β cell function (Zhu et al., 2021a; Hui and Perfetti, 2002; Stoffers et al., 1997). In mammals, Pdx1 knockout results in severe, life-threatening pancreatic developmental defects (Fujitani et al., 2006; Jonsson et al., 1994; Schwitzgebel et al., 2003). To investigate this further, we used CRISPR/Cas9 technology to generate Pdx1 mutants. The gRNA was designed to target the region immediately following the nuclear localization signal within the homeobox domain (Figure 4A). F₀ founders were bred with d/d axolotls to produce F₁ heterozygous embryos, in which we identified various insertions or deletions (indels) at the targeted site in germline-transmitted F₁ animals (Figure 4A). Additionally, we performed intercrosses of individual F₀ founders to obtain “homozygous” mutants, considering both alleles with frameshift mutations as “homozygous”.

Figure 4

Figure 4. The mutations of Pancreatic and duodenal homeobox 1 (Pdx1) result in varying degrees of absence in both the endocrine and exocrine cells of the axolotl pancreas, attributable to different allele mutations. (A) Schematic diagrams of the Pdx1 gene structure and PDX1 protein domains, along with a schematic illustration of the actual effects of Pdx1 mutations on the proteins. CRISPR gRNAs was designed to target the Pdx1 sequence and injected into one-cell-stage eggs, and resulted in mild and severe phenotypes. Axolotls with a complete homeobox in both Pdx1 alleles (Pdx1^s/s) exhibit a severe phenotype. In contrast, a mismatch or frameshift mutation in the homeobox region of both Pdx1 alleles (Pdx1^m/m) results in a milder pancreatic defect. (B) Fasting blood glucose levels in control, heterozygous, and homozygous axolotls show that blood glucose significantly increase in Pdx1^−/− mutants, while no significant change is observed in Pdx1^+/− mutants. (C–E) Anatomical analysis of control (C) and Pdx1 mutant axolotls (blue circle). A complete pancreatic structure, but with a small portion of the pancreas body missing (D). While a severe phenotype is characterized by the development of only a remnant of pancreatic head (E). (F) In situ hybridization for Amy (blue) on cryosections shows that Pdx1 mutants causes a deficiency in the exocrine pancreas. (G) Immunofluorescence for INS (red), GCG (green), combined with DAPI (blue) on cryosections shows INS⁺/GCG⁺ double hormone cells (arrows) appear in Pdx1^m/m mutant axolotls, while Pdx1^s/s mutant axolotls have no endocrine cells. (H) Statistical analysis of endocrine cell counts shows a decrease in the number of β, α, and δ cells following Pdx1 mutations. (I) Statistical analysis of the proportions of β, α, and δ cells in the control, Pdx1^m/m, and Pdx1^s/s groups. (J) Statistical analysis of double hormone cells in the control and Pdx1^m/m groups. (K) Statistical analysis shows that animal growth is severely inhibited in Pdx1^s/s mutant axolotls. Data are the mean ± s. e.m. Scale bar: (C–E) 200 μm; (F) 200 μm; (G) 50 μm.

The “homozygous” mutants survived but predominantly exhibited stunted growth and dwarfism (Figure 4K). Following a 24-h fasting period, we assessed the blood glucose levels of axolotls measuring 7 cm in length, subsequently conducting genotyping. The results indicated that the Pdx1 double allelic mutants, as opposed to heterozygous individuals, exhibited hyperglycaemia compared to the d/d controls (Figure 4B). We then isolated the pancreas from 8-month-old “homozygous” mutants and sibling controls to investigate pancreatic developmental phenotypes. Based on pancreatic morphology, we classified the bi-allelic “homozygous” mutants into two categories: mild (Pdx1^m/m) and severe (Pdx1^s/s) phenotypes. The pancreas of Pdx1^m/m mutants showed relatively normal in size and morphology, while that of Pdx1^s/s mutants was significantly reduced, consisting only of a dense remnant of the pancreatic head (Figures 4C–E).

mRNA in situ hybridization using an antisense Amy probe revealed a notable reduction in exocrine tissue in “homozygous” mutants, particularly pronounced in Pdx1^s/s mutants, alongside an increased proportion of ductal structures (Figure 4F). Further immunofluorescence analysis using antibodies against INS, GCG and SST demonstrated a decrease in endocrine cells in Pdx1^m/m mutants and a near-total depletion in Pdx1^s/s mutants (Figures 4G–I). Intriguingly, the presence of INS⁺/GCG⁺ dihormonal cells and a small number of GCG⁺/SST⁺ dihormonal cells, which were nearly absent in 8-month-old d/d individuals, was noted in Pdx1^m/m mutants (Figure 4J; Supplementary Figure S6). These findings underscore the critical role of Pdx1 in pancreas development and blood glucose regulation in axolotls, consistent with observations in mammals (Offield et al., 1996; Oliver-Krasinski et al., 2009). In summary, Pdx1 mutant axolotls exhibit pancreatic developmental defects and a reduction in both endocrine and exocrine cell populations, characteristics similar to those observed in mammals. However, unlike mammals (Jonsson et al., 1994), these mutant axolotls can survive for long period, despite exhibiting severe pancreatic defects.

3.5 GTT and ITT reveal a slow blood glucose metabolism and a delayed insulin response in axolotl pancreas

GTT and ITT approach are standard methods for assessing glucose metabolism and insulin response. We employed intraperitoneal GTT and ITT to characterize these processes in axolotls. Assessing these processes in smaller models like zebrafish is challenging, so it is particularly relevant to determine their feasibility in axolotls. We initially followed established procedures used in mice (Cassano et al., 2020; Qu et al., 2019; Zhu et al., 2021b), but noted unexpected fluctuations in glucose and insulin levels during the measurement periods. Considering the evidence of the slower glucose metabolism in axolotls suggested by sequencing data, we optimized our experimental protocols for intraperitoneal GTT and ITT. Blood glucose levels were measured at various time points: 0, 1, 2, 6, 12, 24, 36 and 48 h following a 2 mg/g glucose injection during GTT (Figure 5A). In comparison to blood glucose measurement sampling in mice (Small et al., 2022), the intervals for measurements in axolotls are longer. Our findings revealed that the average fasting blood glucose level in adult axolotls was 3.317 ± 0.419 mmol/L 12 h after an intraperitoneal injection of glucose (2 mg/g), blood glucose levels rose to 7.55 ± 0.369 mmol/L. By 36 h post-injection, glucose levels returned to basal level (2.867 ± 0.486 mmol/L). ITT was performed on 6 h fasted axolotls. 5 μL/g insulin glargine was injected into the axolotl, and blood glucose levels were measured at 0, 0.5, 1, 2, 4, 6 and 12 h after insulin injection (Figure 5B). The results showed that blood glucose levels decreased from 3.962 ± 0.134 mmol/L at baseline to 1.525 ± 0.207 mmol/L at the 12-h mark, where levels stabilized. In summary, axolotls exhibit lower blood glucose homeostasis compared to mammals, aligning with the notion that blood glucose levels correlate positively with metabolic rates (Simon et al., 2019; DeFronzo et al., 1979). This aligns with gene expression analysis and further indicates slower glucose metabolism in axolotls.

Figure 5

Figure 5. Functional characterization of axolotl islets. (A) The glucose tolerance test shows that blood glucose levels gradually rise between 0 and 12 h post-injection of a 2 mg/g glucose, followed by a subsequent decrease. (B) The insulin tolerance test shows that blood glucose levels gradually decrease after injection of a 5 μL/g insulin. Data are the mean ± s. e.m.

4 Discussion

In this study, we investigated the fundamental cellular and molecular features of the axolotl pancreas during development and homeostasis. Our findings reveal that axolotl pancreas exhibits typical evolutionary transitional characteristics, positioning it between the rodents and teleost. Notably, the axolotl pancreas also possesses unique traits, such as a slow glucose metabolic rate and distinct endocrine cell state. Furthermore, the ability to conduct GTT and ITT in individual axolotls, combined with their regenerative potential, makes this species a promising model for pancreas research.

During pancreatic evolution, Protochordates such as amphioxus lack a pancreas but possess enteroendocrine cells, including Ins⁺ and Sst⁺ cells (Dai et al., 2024). There are GCG⁺ and SST⁺ cells present in axolotl’s stomach and intestines, but no INS⁺ cells (Supplementary Figure S3). Basally-branching vertebrates such as lampreys exhibit a single-hormone (Ins⁺) islet organ during larval stages, which later develops into a three-hormone (Ins⁺, Sst⁺ and Ppy⁺) islet organ in adults (Youson and Al-Mahrouki, 1999). In Osteichthyes (e.g., lungfish, zebrafish), the pancreas is located near the anterior and dorsal part of the intestine, with the islets encapsulated by exocrine cells (Icardo et al., 2010; Joss, 2010). The zebrafish is a commonly used model organism widely employed in genetic and embryonic development research (Veldman and Lin, 2008). As a tetrapod amphibian, the evolution of the axolotl pancreas occupies an intermediate position between zebrafish and mammals (Boisvert et al., 2013). Unlike zebrafish, where the ventral bud solely contributes to the exocrine pancreas, both buds in axolotls can differentiate into endocrine and exocrine tissues, similar to mammals (Pan and Wright, 2011). The expression of Ptf1a, a pivotal transcription factor for exocrine formation, is localized in the ventral bud of zebrafish (Matsuda, 2018; Lin et al., 2004). However, axolotls retained the ability to express Ptf1a in both buds. It appears that the dorsal bud of the pancreas gradually acquired the capability to develop into both endocrine and exocrine glands during evolution. Although both buds contribute to pancreatic endocrine lineages, the ventral bud in the axolotl does not contribute to β cell production, indicating its transitional position between the fish and mice. Additionally, axolotls have large islets similar to those in mammals, unlike the β cells of Tilapia, which are found in Brockmann bodies or isolated islets along the mesentery, lacking surrounding exocrine tissue (Moss et al., 2009; Xu et al., 2004).

Certain pancreatic cells in axolotls exhibit distinct characteristics that differentiate them from their mammalian counterparts. For instance, PPY and GCG are co-expressed in axolotl pancreas. This phenomenon has been reported in various simple vertebrates, such as Dipnoiformes and Lepisosteus osseus (Groff and Youson, 1998). During pancreatic development in mice, PP cells and α cells share a common progenitor, indicating that PP cells and α cells differentiate into two distinct cell types during evolution. In previous studies, PPY⁺/GCG⁺ bihormonal cells have been reported (Huang et al., 2009; Aragón et al., 2015). In adult wild-type mouse islets, approximately 18.2% of PPY⁺ cells are PPY⁺/GCG⁺ bihormonal cells. In contrast, in human non-diabetic islets, only 0.1%–0.5% of PPY⁺ cells are PPY⁺/GCG⁺ bihormonal cells, a proportion lower than in mice. This suggests that as species evolve, the differentiation of the endocrine lineage becomes more refined to achieve precise regulation. The co-expression of PPY and GCG in the axolotl pancreas suggests the presence of a unique regulatory mechanism that may influence endocrine cell function and glucose homeostasis. It is well known that PPY regulates pancreatic secretion through both endocrine and exocrine tissues (Kojima et al., 2007), while GCG mainly regulates glucose metabolism (Unger, 1971). Their co-expression in the same cell may indicate a coordinated regulation of metabolic and digestive functions, particularly in response to metabolic demands or tissue repair during development or regeneration. Perez-Frances and colleagues reported that bihormonal cells are upregulated after injury and can contribute to endocrine cell populations (Perez-Frances et al., 2021). Therefore, the presence of PPY⁺/GCG⁺ bihormonal cells may serve as a reservoir of pancreatic islet stem cells, and when endocrine function is compromised, these bihormonal cells may act as progenitor cells during development, promoting the regeneration of other islet lineages.

Additionally, ε cells in axolotls do not express the marker GHRL but instead contain a GHRL-like peptide, MLN, believed to have evolved from the same ancestral gene as GHRL (Asakawa et al., 2001; Hara et al., 2018). In zebrafish, Ghrl and Mln are co-expressed in ε cells, whereas Mln is absent in rats and mice (Kawamura et al., 2019; Lavergne et al., 2020). Given the roles of both peptides in regulating gastrointestinal motility and appetite--critical for nutrient digestion and absorption (Kitazawa and Kaiya, 2019; Peeters, 2005), we designated the MLN⁺ pancreatic cells as ε cells in axolotls. Unlike the direct contact observed between GCG⁺ α cells and SST⁺ δ cells with β cells in humans (Abdulreda et al., 2013), most surrounding endocrine cells in axolotls maintain a noticeable distance from β cell clusters. Although the precise implications of this spatial arrangement are unknown, it suggests that endocrine cell functions in axolotls may not be as tightly regulated as in mammals.

Pdx1 is expressed in multipotent progenitor cells of the pancreas and is crucial for maintaining the identity and function of β cells by repressing the reprogramming of α cells in adulthood (Gao et al., 2014). In mice, a homozygous deletion of the second exon of Pdx1 results in the failure of pancreatic development and death within a few days after birth due to hyperglycemia (Offield et al., 1996; Jonsson et al., 1994). In axolotls, following Pdx1 mutation, the pancreatic area in F₀ chimeras was reduced, but the number of endocrine cells remained unaffected, and only a few individuals exhibited an increase in blood glucose. In contrast, F₁ “homozygous” mutants managed to survive despite a nearly absent pancreatic tissue structure, although they were reduced in size compared to controls, reminiscent of pancreas agenesis, neonatal diabetes, and delayed fetal development resulting from Pdx1 homozygous loss in humans (Schwitzgebel et al., 2003). The Pdx family consists of two members, Pdx1 and Pdx2, with reports suggesting that Pdx2 has been lost in tetrapod (Mulley and Holland, 2010). Our sequencing data and the previously assembled axolotl genome support the absence of Pdx2 in axolotls (Schloissnig et al., 2021; Smith et al., 2019). Thus, the loss of Pdx1 is unlikely to be compensated by a homologous gene, implying that other compensatory mechanisms may exist. Given that the knockout target site in Pdx1 is at the end of the DNA-binding domain, following the nuclear localization signal, and considering Pdx1’s many binding sites that regulate pancreatic development, we hypothesize that a complete DNA-binding domain with a 3′UTR frameshift may permit PDX1 to enter the nucleus and function in a “dominant-negative” manner, resulting in a more severe phenotype. In contrast, deletion at the end of the DNA-binding domain could inhibit PDX1 from entering the nucleus, allowing other compensatory transcription factors to function, resulting in a milder phenotype. However, these transcription factors may not regulate endocrine differentiation as precisely as PDX1, potentially leading to an increase in dihormonal cells.

Metabolism in axolotls is characterized by sluggishness, as evidenced by the attenuated expression levels of glucose sensor genes, insulin secretion genes, and glycolytic-related kinases. In addition, the expression of K⁺/Ca²⁺ ion channel protein genes Cbarp, Abcc8 and Kcnj11, and genes involved in glucose metabolism signaling pathways, such as Akt1/2, Prkaa1, and Stk11, are also downregulated in axolotls (Supplementary Figure S4H). This diminished expression renders axolotls less sensitive to glucose and insulin, leading to sluggish blood glucose regulation, as shown by our ITT and GTT results. Some studies suggest that food intake or glucose infusion can promote β cell replication and pluripotent stem cells differentiation (Alonso et al., 2007; Bonner-Weir et al., 1989; Wang et al., 2024), while insulin-resistant states lead to a compensatory increase in β cell mass (Kulkarni et al., 2004). We hypothesize that the unique characteristics of glucose metabolism observed in axolotls may enhance β cell regeneration, warranting further investigation into the mechanisms underlying β cell regeneration.

In summary, we discuss both common and unique pancreatic features between axolotl and other species. The axolotl serves as a valuable model for enriching our understanding of the relationship between ontogeny and phylogeny. Furthermore, the axolotl is well-known for its exceptional regenerative abilities, particularly in limbs, brains and spinal cord. Therefore, studying the regeneration of the axolotl pancreas is highly meaningful. Given the axolotl’s regenerative capacity, future research could investigate its ability to regenerate the pancreas following injury. Exploring stem cell or progenitor sources, the signaling pathways, and transcriptional networks involved in pancreatic regeneration in axolotls could provide valuable insights into tissue repair mechanisms and offer new avenues for regenerative therapies, especially for conditions like diabetes. This represents an exciting direction in the field of regenerative biology, building on the foundation laid by this work.

Data availability statement

The datasets presented in this study can be found in online repositories. The name of the repository is the China National Center for Bioinformation’s Genome Sequence Archive (GSA), and the accession number is CRA019924: https://ngdc.cncb.ac.cn/gsa/s/519FVvW5.

Ethics statement

Ethical approval was not required for the studies involving humans because Human RNA sequences database resource from NCBI. The studies were conducted in accordance with the local legislation and institutional requirements. The human samples used in this study were acquired from RNA sequences from untreated human samples (GSE205853) (Benaglio et al., 2022) were used in this study. Written informed consent to participate in this study was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and the institutional requirements. The animal studies were approved by Guangdong Provincial People’s Hospital. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent was obtained from the owners for the participation of their animals in this study.

Author contributions

HM: Data curation, Investigation, Writing–original draft, Writing–review and editing. GP: Data curation, Methodology, Writing–review and editing. YH: Data curation, Investigation, Methodology, Writing–review and editing. BL: Data curation, Writing–review and editing, Investigation. YZ: Data curation, Investigation, Writing–review and editing. YW: Data curation, Software, Writing–review and editing. WF: Software, Writing–review and editing. YS: Data curation, Writing–review and editing. XP: Software, Writing–review and editing. LS: Data curation, Software, Writing–review and editing. IS: Supervision, Writing–review and editing. YL: Funding acquisition, Project administration, Writing–original draft, Writing–review and editing. JF: Funding acquisition, Supervision, Writing–review and editing.

Funding

The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was supported by BGI grant (BGIRSZ20210002), the National Key R&D Program of China (2021YFA0805000; 2019YFE0106700), the National Natural Science Foundation of China (32070819; 31970782; 92268114), the High-level Hospital Construction Project of Guangdong Provincial People’s Hospital (DFJHBF202103; KJ012021012).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2025.1531903/full#supplementary-material

References

Abdulreda, M., Caicedo, A., and Berggren, P.-O. (2013). A natural body window to study human pancreatic islet cell function and survival. replacement, Regen. & reprogramming 1, 111–122.

PubMed Abstract | Google Scholar

Aguayo-Mazzucato, C., and Bonner-Weir, S. (2018). Pancreatic β cell regeneration as a possible therapy for diabetes. Cell Metab. 27, 57–67. doi:10.1016/j.cmet.2017.08.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Aigha, I. I., and Abdelalim, E. M. (2020). NKX6.1 transcription factor: a crucial regulator of pancreatic β cell development, identity, and proliferation. Stem Cell Res. & Ther. 11, 459. doi:10.1186/s13287-020-01977-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Alonso, L. C., Yokoe, T., Zhang, P., Scott, D. K., Kim, S. K., O'Donnell, C. P., et al. (2007). Glucose infusion in mice: a new model to induce beta-cell replication. Diabetes 56, 1792–1801. doi:10.2337/db06-1513

PubMed Abstract | CrossRef Full Text | Google Scholar

Aragón, F., Karaca, M., Novials, A., Maldonado, R., Maechler, P., and Rubí, B. (2015). Pancreatic polypeptide regulates glucagon release through PPYR1 receptors expressed in mouse and human alpha-cells. Biochimica Biophysica Acta (BBA) - General Subj. 1850, 343–351. doi:10.1016/j.bbagen.2014.11.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Asakawa, A., Inui, A., Kaga, O., Yuzuriha, H., Nagata, T., Ueno, N., et al. (2001). Ghrelin is an appetite-stimulatory signal from stomach with structural resemblance to motilin. Gastroenterology 120, 337–345. doi:10.1053/gast.2001.22158

PubMed Abstract | CrossRef Full Text | Google Scholar

Baron, M., Veres, A., Wolock, S. L., Faust, A. L., Gaujoux, R., Vetere, A., et al. (2016). A single-cell transcriptomic map of the human and mouse pancreas reveals inter- and intra-cell population structure. Cell Syst. 3, 346–360.e4. doi:10.1016/j.cels.2016.08.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Benaglio, P., Zhu, H., Okino, M.-L., Yan, J., Elgamal, R., Nariai, N., et al. (2022). Type 1 diabetes risk genes mediate pancreatic beta cell survival in response to proinflammatory cytokines. Cell Genom. 2, 100214. doi:10.1016/j.xgen.2022.100214

PubMed Abstract | CrossRef Full Text | Google Scholar

Boisvert, C. A., Joss, J. M. P., and Ahlberg, P. E. (2013). Comparative pelvic development of the axolotl (Ambystoma mexicanum) and the Australian lungfish (Neoceratodus forsteri): conservation and innovation across the fish-tetrapod transition. EvoDevo 4, 3. doi:10.1186/2041-9139-4-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Bonner-Weir, S., Deery, D., Leahy, J. L., and Weir, G. C. (1989). Compensatory growth of pancreatic beta-cells in adult rats after short-term glucose infusion. Diabetes. 38, 49–53. doi:10.2337/diab.38.1.49

PubMed Abstract | CrossRef Full Text | Google Scholar

Bonner-Weir, S., Sullivan, B. A., and Weir, G. C. (2015). Human islet morphology revisited: human and rodent islets are not so different after all. J. Histochem Cytochem. 63, 604–612. doi:10.1369/0022155415570969

PubMed Abstract | CrossRef Full Text | Google Scholar

Bryant, D. M., Johnson, K., DiTommaso, T., Tickle, T., Couger, M. B., Payzin-Dogru, D., et al. (2017). A tissue-mapped axolotl de novo transcriptome enables identification of limb regeneration factors. Cell Rep. 18, 762–776. doi:10.1016/j.celrep.2016.12.063

PubMed Abstract | CrossRef Full Text | Google Scholar

Burlison, J. S., Long, Q., Fujitani, Y., Wright, C. V., and Magnuson, M. A. (2008). Pdx-1 and Ptf1a concurrently determine fate specification of pancreatic multipotent progenitor cells. Dev. Biol. 316, 74–86. doi:10.1016/j.ydbio.2008.01.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Cano-Martínez, A., Vargas-González, A., Guarner-Lans, V., Prado-Zayago, E., León-Olea, M., and Nieto-Lima, B. (2010). Functional and structural regeneration in the axolotl heart (Ambystoma mexicanum) after partial ventricular amputation. Arch. Cardiol. Mex. 80, 79–86.

PubMed Abstract | Google Scholar

Cassano, V., Leo, A., Tallarico, M., Nesci, V., Cimellaro, A., Fiorentino, T. V., et al. (2020). Metabolic and cognitive effects of ranolazine in type 2 diabetes mellitus: data from an in vivo model. vivo model 12, 382. doi:10.3390/nu12020382

PubMed Abstract | CrossRef Full Text | Google Scholar

Cole, J. B., and Florez, J. C. (2020). Genetics of diabetes mellitus and diabetes complications. Nat. Rev. Nephrol. 16, 377–390. doi:10.1038/s41581-020-0278-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Cura Costa, E., Otsuki, L., Rodrigo Albors, A., Tanaka, E. M., and Chara, O. (2021). Spatiotemporal control of cell cycle acceleration during axolotl spinal cord regeneration. Elife. 10, e55665. doi:10.7554/eLife.55665

PubMed Abstract | CrossRef Full Text | Google Scholar

Dai, Y., Pan, R., Pan, Q., Wu, X., Cai, Z., Fu, Y., et al. (2024). Single-cell profiling of the amphioxus digestive tract reveals conservation of endocrine cells in chordates. Sci. Adv. 10, eadq0702. doi:10.1126/sciadv.adq0702

PubMed Abstract | CrossRef Full Text | Google Scholar

DeFronzo, R. A., Tobin, J. D., and Andres, R. (1979). Glucose clamp technique: a method for quantifying insulin secretion and resistance. Am. J. Physiol. 237, E214–E223. doi:10.1152/ajpendo.1979.237.3.E214

PubMed Abstract | CrossRef Full Text | Google Scholar

de Jesus, D. S., Mak, T. C., Wang, Y.-F., von Ohlen, Y., Bai, Y., Kane, E., et al. (2021). Dysregulation of the Pdx1/Ovol2/Zeb2 axis in dedifferentiated β-cells triggers the induction of genes associated with epithelial-mesenchymal transition in diabetes. Mol. Metab. 53, 101248. doi:10.1016/j.molmet.2021.101248

PubMed Abstract | CrossRef Full Text | Google Scholar

Echeverri, K., and Tanaka, E. M. (2002). Ectoderm to mesoderm lineage switching during axolotl tail regeneration. Science 298, 1993–1996. doi:10.1126/science.1077804

PubMed Abstract | CrossRef Full Text | Google Scholar

Fei, J.-F., Lou, W. P.-K., Knapp, D., Murawala, P., Gerber, T., Taniguchi, Y., et al. (2018). Application and optimization of CRISPR–Cas9-mediated genome engineering in axolotl (Ambystoma mexicanum). Nat. Protoc. 13, 2908–2943. doi:10.1038/s41596-018-0071-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Fei, J.-F., Schuez, M., Knapp, D., Taniguchi, Y., Drechsel, D. N., and Tanaka, E. M. (2017). Efficient gene knockin in axolotl and its use to test the role of satellite cells in limb regeneration. Proc. Natl. Acad. Sci. 114, 12501–12506. doi:10.1073/pnas.1706855114

PubMed Abstract | CrossRef Full Text | Google Scholar

Fitzner, B., Lange, A., Müller, S., and Jaster, R. (2013). Cdkn1a is a key mediator of rat pancreatic stellate cell senescence. Pancreatology 13, 254–262. doi:10.1016/j.pan.2013.03.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Fujitani, Y., Fujitani, S., Boyer, D. F., Gannon, M., Kawaguchi, Y., Ray, M., et al. (2006). Targeted deletion of a cis-regulatory region reveals differential gene dosage requirements for Pdx1 in foregut organ differentiation and pancreas formation. Genes Dev. 20, 253–266. doi:10.1101/gad.1360106

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, T., McKenna, B., Li, C., Reichert, M., Nguyen, J., Singh, T., et al. (2014). Pdx1 maintains β cell identity and function by repressing an α cell program. Cell Metab. 19, 259–271. doi:10.1016/j.cmet.2013.12.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Ghila, L., Furuyama, K., Grey, S. T., Scholz, H., and Chera, S. J. F. i. G. (2022). Beta-cell fate: from gene circuits to disease mechanisms. Lausanne, Switzerland: Frontiers Media SA.

Google Scholar

Groff, K. E., and Youson, J. H. (1998). Fine structure and immunocytochemistry of cells within the endocrine pancreas of the gar (Lepisosteus osseus). Can. J. Zoology 76, 6–18. doi:10.1139/z97-161

CrossRef Full Text | Google Scholar

Hara, Y., Yamaguchi, K., Onimaru, K., Kadota, M., Koyanagi, M., Keeley, S. D., et al. (2018). Shark genomes provide insights into elasmobranch evolution and the origin of vertebrates. Nat. Ecol. & Evol. 2, 1761–1771. doi:10.1038/s41559-018-0673-5

PubMed Abstract | CrossRef Full Text | Google Scholar

He, J., Mo, D., Chen, J., and Luo, L. (2020). Combined whole-mount fluorescence in situ hybridization and antibody staining in zebrafish embryos and larvae. Nat. Protoc. 15, 3361–3379. doi:10.1038/s41596-020-0376-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Holland, A. M., Gonez, L. J., Naselli, G., MacDonald, R. J., and Harrison, L. C. (2005). Conditional expression demonstrates the role of the homeodomain transcription factor Pdx1 in maintenance and regeneration of beta-cells in the adult pancreas. Diabetes 54, 2586–2595. doi:10.2337/diabetes.54.9.2586

PubMed Abstract | CrossRef Full Text | Google Scholar

Hossain, M. J., Al-Mamun, M., and Islam, M. R. (2024). Diabetes mellitus, the fastest growing global public health concern: early detection should be focused. Health Sci. Rep. 7, e2004. doi:10.1002/hsr2.2004

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, Y., Sun, M., Jiang, M., and Fu, B. Y. (2009). Immunohistochemical localization of glucagon and pancreatic polypeptide on rat endocrine pancreas: coexistence in rat islet cells. Eur. J. Histochem. 53, e10. doi:10.4081/ejh.2009.e10

PubMed Abstract | CrossRef Full Text | Google Scholar

Hui, H., and Perfetti, R. (2002). Pancreas duodenum homeobox-1 regulates pancreas development during embryogenesis and islet cell function in adulthood. Eur. J. Endocrinol. 146, 129–141. doi:10.1530/eje.0.1460129

PubMed Abstract | CrossRef Full Text | Google Scholar

Icardo, J. M., Wong, W. P., Colvee, E., Loong, A. M., and Ip, Y. (2010). The anatomy of the gastrointestinal tract of the African lungfish, Protopterus annectens. Protopterus annectens 293, 1146–1154. doi:10.1002/ar.21154

PubMed Abstract | CrossRef Full Text | Google Scholar

Jonsson, J., Carlsson, L., Edlund, T., and Edlund, H. (1994). Insulin-promoter-factor 1 is required for pancreas development in mice. Nature. 371, 606–609. doi:10.1038/371606a0

PubMed Abstract | CrossRef Full Text | Google Scholar

Joss, J. (2010). The lungfish endocrine system in The Biology of lungfishes (Enfield, New Hampshire: CRC Press Taylor & Francis Group, Science Publishers).

Google Scholar

Juran, C. M., Zvirblyte, J., Cheng-Campbell, M., Blaber, E. A., and Almeida, E. A. (2021). Cdkn1a deletion or suppression by cyclic stretch enhance the osteogenic potential of bone marrow mesenchymal stem cell-derived cultures. Stem Cell Res. 56, 102513. doi:10.1016/j.scr.2021.102513

PubMed Abstract | CrossRef Full Text | Google Scholar

Kawamura, T., Matsuura, B., Miyake, T., Abe, M., Ikeda, Y., and Hiasa, Y. (2019). Effects of motilin receptor agonists and ghrelin in human motilin receptor transgenic mice. Int. J. Mol. Sci. 20, 1521. doi:10.3390/ijms20071521

PubMed Abstract | CrossRef Full Text | Google Scholar

Kimmel, R., and Meyer, D. (2016). Zebrafish pancreas as a model for development and disease. Methods cell Biol. 134, 431–461. doi:10.1016/bs.mcb.2016.02.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Kitazawa, T., and Kaiya, H. (2019). Regulation of gastrointestinal motility by motilin and ghrelin in vertebrates. Front. Endocrinol. (Lausanne). 10, 278. doi:10.3389/fendo.2019.00278

PubMed Abstract | CrossRef Full Text | Google Scholar

Kojima, S., Ueno, N., Asakawa, A., Sagiyama, K., Naruo, T., Mizuno, S., et al. (2007). A role for pancreatic polypeptide in feeding and body weight regulation. Peptides 28, 459–463. doi:10.1016/j.peptides.2006.09.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Kulkarni, R. N., Jhala, U. S., Winnay, J. N., Krajewski, S., Montminy, M., and Kahn, C. R. (2004). PDX-1 haploinsufficiency limits the compensatory islet hyperplasia that occurs in response to insulin resistance. J. Clin. Invest. 114, 828–836. doi:10.1172/JCI21845

PubMed Abstract | CrossRef Full Text | Google Scholar

Lavergne, A., Tarifeño-Saldivia, E., Pirson, J., Reuter, A.-S., Flasse, L., Manfroid, I., et al. (2020). Pancreatic and intestinal endocrine cells in zebrafish share common transcriptomic signatures and regulatory programmes. BMC Biol. 18, 109–119. doi:10.1186/s12915-020-00840-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Lewis, B., and Mao, J. (2023). “Development of the pancreas and related structures,” in The pancreas: an integrated textbook of basic science, medicine, and surgery, 1–8.

Google Scholar

Li, Y., Jia, Z., Zhang, S., and He, X. (2021). Progress in gene-editing technology of zebrafish. Biomolecules 11, 1300. doi:10.3390/biom11091300

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, J. W., Biankin, A. V., Horb, M. E., Ghosh, B., Prasad, N. B., Yee, N. S., et al. (2004). Differential requirement for ptf1a in endocrine and exocrine lineages of developing zebrafish pancreas. Dev. Biol. 270, 474–486. doi:10.1016/j.ydbio.2004.02.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Longnecker, D. (2014). Anatomy and histology of the pancreas. Pancreapedia: Exocrine Pancreas Knowledge Base.

Google Scholar

Lu, B. (2023). Evolutionary insights into the relationship of frogs, salamanders, and caecilians and their adaptive traits, with an emphasis on salamander regeneration and longevity. Animals 13, 3449. doi:10.3390/ani13223449

PubMed Abstract | CrossRef Full Text | Google Scholar

Lust, K., Maynard, A., Gomes, T., Fleck, J. S., Camp, J. G., Tanaka, E. M., et al. (2022). Single-cell analyses of axolotl telencephalon organization, neurogenesis, and regeneration. Science 377, eabp9262. doi:10.1126/science.abp9262

PubMed Abstract | CrossRef Full Text | Google Scholar

Matsuda, H. (2018). Zebrafish as a model for studying functional pancreatic β cells development and regeneration. Dev. Growth Differ. 60, 393–399. doi:10.1111/dgd.12565

PubMed Abstract | CrossRef Full Text | Google Scholar

Mi, J., Ren, L., and Andersson, O. (2024). Leveraging zebrafish to investigate pancreatic development, regeneration, and diabetes. Trends Mol. Med. 30, 932–949. doi:10.1016/j.molmed.2024.05.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Moss, J. B., Koustubhan, P., Greenman, M., Parsons, M. J., Walter, I., and Moss, L. G. (2009). Regeneration of the pancreas in adult zebrafish. Diabetes 58, 1844–1851. doi:10.2337/db08-0628

PubMed Abstract | CrossRef Full Text | Google Scholar

Mulley, J. F., and Holland, P. W. (2010). Parallel retention of Pdx2 genes in cartilaginous fish and coelacanths. Mol. Biol. Evol. 27, 2386–2391. doi:10.1093/molbev/msq121

PubMed Abstract | CrossRef Full Text | Google Scholar

Muraro, M. J., Dharmadhikari, G., Grün, D., Groen, N., Dielen, T., Jansen, E., et al. (2016). A single-cell transcriptome atlas of the human pancreas. Cell Syst. 3, 385–394. e3. doi:10.1016/j.cels.2016.09.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Offield, M. F., Jetton, T. L., Labosky, P. A., Ray, M., Stein, R. W., Magnuson, M. A., et al. (1996). PDX-1 is required for pancreatic outgrowth and differentiation of the rostral duodenum. Development 122, 983–995. doi:10.1242/dev.122.3.983

PubMed Abstract | CrossRef Full Text | Google Scholar

Oliver-Krasinski, J. M., Kasner, M. T., Yang, J., Crutchlow, M. F., Rustgi, A. K., Kaestner, K. H., et al. (2009). The diabetes gene Pdx1 regulates the transcriptional network of pancreatic endocrine progenitor cells in mice. J. Clin. Invest. 119, 1888–1898. doi:10.1172/JCI37028

PubMed Abstract | CrossRef Full Text | Google Scholar

Pan, F. C., and Wright, C. (2011). Pancreas organogenesis: from bud to plexus to gland. Dev. Dyn. 240, 530–565. doi:10.1002/dvdy.22584

PubMed Abstract | CrossRef Full Text | Google Scholar

Peeters, T. (2005). Ghrelin: a new player in the control of gastrointestinal functions. Gut. 54, 1638–1649. doi:10.1136/gut.2004.062604

PubMed Abstract | CrossRef Full Text | Google Scholar

Perez-Frances, M., van Gurp, L., Abate, M. V., Cigliola, V., Furuyama, K., Bru-Tari, E., et al. (2021). Pancreatic Ppy-expressing γ-cells display mixed phenotypic traits and the adaptive plasticity to engage insulin production. Nat. Commun. 12, 4458. doi:10.1038/s41467-021-24788-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Qu, Y.-L., Deng, C.-H., Luo, Q., Shang, X.-Y., Wu, J.-X., Shi, Y., et al. (2019). Arid1a regulates insulin sensitivity and lipid metabolism. EBioMedicine 42, 481–493. doi:10.1016/j.ebiom.2019.03.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Rankin, M. M., and Kushner, J. A. J. D. (2009). Adaptive beta-cell proliferation is severely restricted with advanced age. Diabetes 58, 1365–1372. doi:10.2337/db08-1198

PubMed Abstract | CrossRef Full Text | Google Scholar

Röder, P. V., Wu, B., Liu, Y., and Han, W. (2016). Pancreatic regulation of glucose homeostasis. Exp. & Mol. Med. 48, e219. doi:10.1038/emm.2016.6

PubMed Abstract | CrossRef Full Text | Google Scholar

Roep, B. O. (2020). There is something about insulin granules. Diabetes 69, 2575–2577. doi:10.2337/dbi20-0039

PubMed Abstract | CrossRef Full Text | Google Scholar

Sasaki, H., Saisho, Y., Inaishi, J., and Itoh, H. (2021). Revisiting regulators of human β-cell mass to achieve β-cell-centric approach toward type 2 diabetes. J. Endocr. Soc. 5, bvab128. doi:10.1210/jendso/bvab128

PubMed Abstract | CrossRef Full Text | Google Scholar

Schloissnig, S., Kawaguchi, A., Nowoshilow, S., Falcon, F., Otsuki, L., Tardivo, P., et al. (2021). The giant axolotl genome uncovers the evolution, scaling, and transcriptional control of complex gene loci. Proc. Natl. Acad. Sci. 118, e2017176118. doi:10.1073/pnas.2017176118

PubMed Abstract | CrossRef Full Text | Google Scholar

Schwitzgebel, V. M., Mamin, A., Brun, T., Ritz-Laser, B., Zaiko, M., Maret, A., et al. (2003). Agenesis of human pancreas due to decreased half-life of insulin promoter factor 1. J. Clin. Endocrinol. Metab. 88, 4398–4406. doi:10.1210/jc.2003-030046

PubMed Abstract | CrossRef Full Text | Google Scholar

Simon, K., and Wittmann, I. (2019). Can blood glucose value really be referred to as a metabolic parameter? Rev. Endocr. Metabolic Disord. 20, 151–160. doi:10.1007/s11154-019-09504-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Slack, J. J. D. (1995). Developmental biology of the pancreas. Development 121, 1569–1580. doi:10.1242/dev.121.6.1569

PubMed Abstract | CrossRef Full Text | Google Scholar

Small, L., Ehrlich, A., Iversen, J., Ashcroft, S. P., Trošt, K., Moritz, T., et al. (2022). Comparative analysis of oral and intraperitoneal glucose tolerance tests in mice. Mol. Metab. 57, 101440. doi:10.1016/j.molmet.2022.101440

PubMed Abstract | CrossRef Full Text | Google Scholar

Smith, J. J., Timoshevskaya, N., Timoshevskiy, V. A., Keinath, M. C., Hardy, D., and Voss, S. R. (2019). A chromosome-scale assembly of the axolotl genome. Genome Res. 29, 317–324. doi:10.1101/gr.241901.118

PubMed Abstract | CrossRef Full Text | Google Scholar

Sørensen, P. L., Dittrich, A. L., and Lauridsen, H. (2022). A novel animal model in diabetes: regeneration of β cells in the axolotl salamander? 36.

Google Scholar

Steiner, D. J., Kim, A., Miller, K., and Hara, M. J. I. (2010). Pancreatic islet plasticity: interspecies comparison of islet architecture and composition. Islets. 2, 135–145. doi:10.4161/isl.2.3.11815

PubMed Abstract | CrossRef Full Text | Google Scholar

Stoffers, D. A., Zinkin, N. T., Stanojevic, V., Clarke, W. L., and Habener, J. F. (1997). Pancreatic agenesis attributable to a single nucleotide deletion in the human IPF1 gene coding sequence. Nat. Genet. 15, 106–110. doi:10.1038/ng0197-106

PubMed Abstract | CrossRef Full Text | Google Scholar

Suetsugu-Maki, R., Maki, N., Nakamura, K., Sumanas, S., Zhu, J., Del Rio-Tsonis, K., et al. (2012). Lens regeneration in axolotl: new evidence of developmental plasticity. BMC Biol. 10, 103–108. doi:10.1186/1741-7007-10-103

PubMed Abstract | CrossRef Full Text | Google Scholar

Susan, B.-W. (2004). “Endocrine pancreas,” in Encyclopedia of gastroenterology. Editor L. R. JOHNSON (New York: Elsevier).

Google Scholar

Teame, T., Zhang, Z., Ran, C., Zhang, H., Yang, Y., Ding, Q., et al. (2019). The use of zebrafish (Danio rerio) as biomedical models. Anim. Front. 9, 68–77. doi:10.1093/af/vfz020

PubMed Abstract | CrossRef Full Text | Google Scholar

Teta, M., Long, S. Y., Wartschow, L. M., Rankin, M. M., and Kushner, J. A. J. D. (2005). Very slow turnover of beta-cells in aged adult mice. Diabetes 54, 2557–2567. doi:10.2337/diabetes.54.9.2557

PubMed Abstract | CrossRef Full Text | Google Scholar

Unger, R. H. (1971). Glucagon physiology and pathophysiology. N. Engl. J. Med. 285, 443–449. doi:10.1056/nejm197108192850806

PubMed Abstract | CrossRef Full Text | Google Scholar

Veldman, M. B., and Lin, S. (2008). Zebrafish as a developmental model organism for pediatric research. Pediatr. Res. 64, 470–476. doi:10.1203/PDR.0b013e318186e609

PubMed Abstract | CrossRef Full Text | Google Scholar

Vieira, W. A., Wells, K. M., and McCusker, C. D.(2020). Advancements to the axolotl model for regeneration and aging. Gerontology 66, 212–222. doi:10.1159/000504294

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, C., Abadpour, S., Olsen, P. A., Wang, D., Stokowiec, J., Chera, S., et al. (2024). Glucose concentration in regulating induced pluripotent stem cells differentiation toward insulin-producing cells. Cells 37, 11900. doi:10.3389/ti.2024.11900

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J., Yang, C., Kong, F., Zhang, Z., Ji, S., Sun, G., et al. (2022). Cebpb regulates skeletal stem cell osteogenic differentiation and fracture healing via the WNT/β-Catenin pathway. Stem Cells Int. 2022, 2091615–2091616. doi:10.1155/2022/2091615

PubMed Abstract | CrossRef Full Text | Google Scholar

Wei, X., Fu, S., Li, H., Liu, Y., Wang, S., Feng, W., et al. (2022). Single-cell Stereo-seq reveals induced progenitor cells involved in axolotl brain regeneration. Science 377, eabp9444. doi:10.1126/science.abp9444

PubMed Abstract | CrossRef Full Text | Google Scholar

Weng, C.-C., Lin, Y.-C., and Cheng, K.-H. (2019). The use of genetically engineered mouse models for studying the function of mutated driver genes in pancreatic cancer. J. Clin. Med. 8, 1369. doi:10.3390/jcm8091369

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, B.-Y., Morrison, C. M., Yang, H., Wright Jr, J. R., and endocrinology, c. (2004). Tilapia islet grafts are highly alloxan-resistant. Gen. Comp. Endocrinol. 137, 132–140. doi:10.1016/j.ygcen.2004.02.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, B., Covington, B. A., and Chen, W. (2020). In vivo generation and regeneration of β cells in zebrafish. Cell Regen. 9, 9. doi:10.1186/s13619-020-00052-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Youson, J. H., and Al-Mahrouki, A. A. (1999). Ontogenetic and phylogenetic development of the endocrine pancreas (islet organ) in fishes. Gen. Comp. Endocrinol. 116 (3), 303–335. doi:10.1016/j.xpro.2021.100728

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, X.-X., Qiu, W.-L., Yang, L., Wang, Y.-C., He, M.-Y., Wang, D., et al. (2021). Sequential progenitor states mark the generation of pancreatic endocrine lineages in mice and humans. Cell Res. 31, 886–903. doi:10.1038/s41422-021-00486-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, X.-X., and Xu, C.-R. (2020). Understanding generation and regeneration of pancreatic β cells from a single-cell perspective. Development. 147, dev179051. doi:10.1242/dev.179051

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, X., Oguh, A., Gingerich, M. A., Soleimanpour, S. A., Stoffers, D. A., and Gannon, M. (2021a). Cell cycle regulation of the Pdx1 transcription factor in developing pancreas and insulin-producing β-cells. Diabetes. 70, 903–916. doi:10.2337/db20-0599

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, Y.-X., Zhou, Y.-C., Zhang, Y., Sun, P., Chang, X.-A., and Han, X. (2021b). Protocol for in vivo and ex vivo assessments of glucose-stimulated insulin secretion in mouse islet β cells. Star. Protoc. 2, 100728. doi:10.1016/j.xpro.2021.100728

PubMed Abstract | CrossRef Full Text | Google Scholar

Glossary

Nkx6.1 NK6 Homeobox 1

Nkx2.2 NK2 Homeobox 2

Foxa2 Forkhead Box A2

Sox9 SRY-Box Transcription Factor 9

Neurog3 Neurogenin 3

Isl1 ISL LIM Homeobox 1

Neurod1 Neuronal Differentiation 1

Pax6 Paired Box 6

Arx Aristaless Related Homeobox

MafB MAF BZIP Transcription Factor B

Peg10 Paternally Expressed 10

Gng12 G Protein Subunit Gamma 12

MafA MAF BZIP Transcription Factor A

Ptprn Protein Tyrosine Phosphatase Receptor Type N

Slc30a8 Solute Carrier Family 30 Member 8

Rbp4 Retinol Binding Protein 4

Sox4 SRY-Box Transcription Factor 4

Slc2a1 Solute Carrier Family 2 Member 1

Slc2a2 Solute Carrier Family 2 Member 2

Stx1a Syntaxin 1A

Vamp8 Vesicle Associated Membrane Protein 8

Pfkl Phosphofructokinase, Liver Type

Pfkm Phosphofructokinase, Muscle

Pkm Pyruvate Kinase M1/2

Gck Glucokinase

Pgk1 Phosphoglycerate Kinase 1

Ldhb Lactate Dehydrogenase B

Slc7a5 Solute Carrier Family 7 Member 5

Cdkn1a Cyclin Dependent Kinase Inhibitor 1A

Cebpb CCAAT Enhancer Binding Protein Beta

Cbarp CACN Subunit Beta Associated Regulatory Protein

Abcc8 ATP Binding Cassette Subfamily C Member 8

Kcnj11 Potassium Inwardly Rectifying Channel Subfamily J Member 11

Akt1/2 AKT Serine/Threonine Kinase 1/2

Prkaa1 Protein Kinase AMP-Activated Catalytic Subunit Alpha 1

Stk11 Serine/Threonine Kinase 11