Skip to main content

Comparative iTRAQ proteomics revealed proteins associated with lobed fin regeneration in Bichirs



Polypterus senegalus can fully regenerate its pectoral lobed fins, including a complex endoskeleton, with remarkable precision. However, despite the enormous potential of this species for use in medical research, its regeneration mechanisms remain largely unknown.


To identify the differentially expressed proteins (DEPs) during the early stages of lobed fin regeneration in P. senegalus, we performed a differential proteomic analysis using isobaric tag for relative and absolute quantitation (iTRAQ) approach based quantitative proteome from the pectoral lobed fins at 3 time points. Furthermore, we validated the changes in protein expression with multiple-reaction monitoring (MRM) analysis.


The experiment yielded a total of 3177 proteins and 15,091 unique peptides including 1006 non-redundant (nr) DEPs. Of these, 592 were upregulated while 349 were downregulated after lobed fin amputation when compared to the original tissue. Bioinformatics analyses showed that the DEPs were mainly associated with Ribosome and RNA transport, metabolic, ECM-receptor interaction, Golgi and endoplasmic reticulum, DNA replication, and Regulation of actin cytoskeleton.


To our knowledge, this is the first proteomic research to investigate alterations in protein levels and affected pathways in bichirs’ lobe-fin/limb regeneration. In addition, our study demonstrated a highly dynamic regulation during lobed fin regeneration in P. senegalus. These results not only provide a comprehensive dataset on differentially expressed proteins during the early stages of lobe-fin/limb regeneration but also advance our understanding of the molecular mechanisms underlying lobe-fin/limb regeneration.


The mystery of limb regeneration in vertebrates has attracted the attention of scientists for many years. The newts and salamanders have one iconic capacity of complete regeneration of all appendages [1]. Anurans have only restricted regeneration capacity, and usually lose it after metamorphosis [2]. Mammals have the capacity of regeneration of the digit tip, but not limb [3]. Polypterus regenerates its lobed pectoral fins with remarkable accuracy [4]. Advances in the study of the development of regenerated limbs have been reported, but why certain species are able to regenerate their limbs is still unknown. In addition, the molecular mechanisms in bichirs lobed fins (endoskeleton) regeneration process remain unknown. It is possible that if we comprehensively analyze the molecular mechanisms involved in the regeneration of lobed pectoral fins in bichirs, we can more easily understand the process of complete lobe-fin/limb regeneration ability.

Salamanders seem to be the only living tetrapods that can regenerate full limbs [5]. However, among vertebrates, the paired lobe-fins/limbs of adult lungfishes and bichirs are also able to fully regenerate [4, 6, 7]. It is notable that the lobe-fins/limbs of bichirs include a bony endoskeleton that differs from the dermal exoskeleton of the fins of teleost fishes, although bichirs were always classified as Actinopterygii [8]. In addition, all three of these species live in water-land transitional environments.

Polypterus senegalus, also known as the Senegal bichir or the gray bichir, is a prototypical fish species in the Polypterus genus. Although the oldest fossil polypteriformes records are only known from 112 to 99 Mya [9], some phylogenetic hypotheses predict a divergence age of 385 to 426.8 Mya [8, 10]. Polypterus shares many characteristics with amphibians in body structure, such as paired lungs [11,12,13,14]. In addition, it has been observed crawling on land both in the lab [15] and in the open [16] by its complex lobed fins, which have equivalent radius, ulna, radiale and ulnare [4, 14, 17].

Newt is remarkable in terms of its regenerative ability. However, it is not as convenient as other model organisms. To some degree, such a phenomenon is caused by their enormous genome size, about 25–50 pg (C value) that is approximately 10 times of that of humans (C value 3.5 pg). In addition, the large introns in newt genes may also reduce the value of this group for research compared to other tetrapods [18]. Previous studies have found that P. senegalus is able to initiate lobe-fin/limb regeneration similar to amphibian, with the formation of a blastema [4, 19]. P. senegalus has a relative small genome size (~ 3.4 Gbp), and short time of reproductive cycle (could be as short as 6 months) [11], making it a better candidate as a model organism for limb/lobe-fin regeneration in vertebrates. Polypteriformes are characterized with derived and plesiomorphic features, which make these fishes attractive subjects for evolutionary and developmental comparisons of lobe-fin/limb regeneration. The research herein was targeted at identifying candidate genes and molecular mechanisms that might have played a vital role in lobe-fin/limb regeneration.


Growth curve

The overall workflow and animal sample for the present study were displayed in Fig. 1a and b, respectively. The green curve in Fig. 1c represents the growth of the regenerating forelimbs of P. senegalus over 15 weeks, and the sampling time points of 0 dpa (0D_L and 0D_R), 4 dpa (4D_B) and 12 dpa (12D_D) are marked. The regenerating forelimbs of P. senegalus exhibited a rapid growth rate in the early stage, and after 9 weeks, it plateaued and persisted for a long time. During the plateaued stage, the growth rate was low (only 0.35 mm over 5 weeks). Photos of regenerating lobed fins at key time points were obtained (Fig. 1d). In order to provide an intuitive understanding of the lobed fin regeneration process, the histologic features of stump tissue in 0 dpa, 4 dpa, and 12 dpa were shown in Fig. 1e.

Fig. 1
figure 1

Study sampling and experimental scheme. a The overall workflow of this study. b A photo of a live specimen of P. senegalus with a homemade fish-holder. c Growth curve of P. senegalus pectoral lobed fin after amputation. The sampling time points of 0 dpa (0D_L and 0D_R), 4 dpa (4D_B) and 12 dpa (12D_D) are labeled. d P. senegalus lobed fins at different regeneration stages. e Tissues of three time points were shown by Alcian blue- and alizarin red- staining, and feulgen staining on paraffin tissue sections. The position of the section was shown on the X-ray photo. Scale bars: 0.5 mm. dpa: days post amputation

Overview of the library and quantitative proteomics analysis

An 8-label isobaric tags for relative and absolute quantification (iTRAQ) experiment was performed for eight samples of the pectoral fins of bichirs at three time points, using P. senegalus genome data as background library. The overall proteomic results, such as basic descriptive statistics, peptide length distribution, protein coverage distribution, and mass errors of all of the identified peptides, are presented in Additional file 1. A total of 3177 proteins and 15,091 unique peptides were quantified in our experiment (Additional file 2), representing 15.4% of the 20,582 predicted proteins in the P. senegalus genome. We determined that 1006 non-redundant (nr) proteins were differentially regulated during the regeneration progress. Among these proteins, 592 (4D_B/0D_L) and 87 (12D_D/4D_B) were upregulated while 349 (4D_B/0D_L) and 55 (12D_D/4D_B) were downregulated after lobed fin amputation (Additional file 3).

Comparison of replicate datasets—Using the quantitative value of proteins in the duplicate data, the degree of variation and the distribution of statistical variation were calculated. All of the mean CV values were less than 0.15, revealing a relatively high repeatability among multiple samples (Additional file 4: A), and the protein ratio distribution is shown in Additional file 4: B. On the basis of protein ratios transformed by log2, we conducted a linear regression analysis to compare the two experimental replicates, so as to evaluate the reproducibility of the iTRAQ proteomic results. R2 values showed the strong linear correlation of the two experimental replicates for both 4D_B and 12D_D (Additional file 4: C). As indicated by these results, the level of reproducibility between replicate datasets was high.

Functional classification of differentially expressed proteins

It was shown from Fig. 2 (only p < 0.001 and top 10 terms were showed) that at both 4D_B and 12D_D, five KEGG pathways were regulated differently, which was displayed in the analysis result of KEGG pathway enrichment: “Ribosome” (ko03010), “Oxidative phosphorylation” (ko00190), “Parkinson’s disease” (ko05012), “Huntington’s disease” (ko05016), “Alzheimer’s disease” (ko05010) and “Metabolic pathways” (ko01100) (Additional file 5: A).

Fig. 2
figure 2

Enrichment analysis of the differentially expressed proteins in 4D_B and 12D_D. GO terms and KEGG pathways enriched in the differentially expressed proteins at 4 and 12 days post amputation. 12D_D VS 0D_L was shown in red. 4D_B VS 0D_L was shown in green

After GO term enrichment analysis, we got several biological processes, including “cellular respiration translational termination”, “translation”, “translational initiation”, “oxidation-reduction process” and “generation of precursor metabolites and energy” were significantly enriched in both 4D_B and 12D_D, suggesting that the processes above are very active during the regeneration process (Additional file 5: B).

The GO molecular function terms, including “ribosome structural constituent,” “structural molecule activity,” “oxidoreductase activity” and “cofactor binding,” and GO cellular component terms, including “mitochondrial inner membrane,” “organelle inner membrane,” “mitochondrion,” “cytosolic ribosome,” “ribosome” and “respiratory chain,” were significantly enriched in both 4D_B and 12D_D (Additional file 5: C and D).

Protein−protein interaction network analysis

To further understand the biological process or pathway regulated in the early stage of lobe fin regeneration, a protein interaction network of DEPs in 4D_B/0D_L was performed by STRING with Xenopus Silurana as background organism. As a result, 697 nodes and 8263 edges were retrieved. After the computing by MCODE in Cytoscape, 16 clusters were obtained (Additional files 6 and 7). And the top 10 clusters were displayed here. Cluster 1 and 9 consists of proteins involved in ribosome and RNA transport. Cluster 2–6 consists of proteins involved in oxidative phosphorylation and metabolic pathways. Cluster 7–8 consists of proteins involved in ECM-receptor interaction. Cluster 10 consists of proteins involved in Golgi and endoplasmic reticulum. It is worth noting that cluster 7 contains three integrins, and the main signaling pathways involved are DNA replication, Biosynthesis of amino acids, Regulation of actin cytoskeleton, and ECM-receptor interaction (Fig. 3).

Fig. 3
figure 3

Protein-protein interaction network of DEPs in cluster 7 by STRING and MCODE. Proteins that are associated with each other are linked by an edge. The up-regulated DEPs (4D_B/0D_L) in 4D_B sample were in red, and the down-regulated DEPs were in green. Integrins were displayed in circles. The main pathway descriptions of these DEPs were from STRING

Validation of changes in protein expression using multiple-reaction monitoring analysis

Multiple-reaction monitoring (MRM) was used to provide further independent validation of the expression status of 30 of the quantified proteins. Finally, 27 results from MRM were in line with the quantitative proteomic results (Fig. 4). Details of the dynamic MRM results, including the fold changes, normalized intensities and p-values of 30 proteins are given in Additional file 8.

Fig. 4
figure 4

Validation of some proteins. Correlation between MRM and iTRAQ data of 30 proteins. 4D/0D, 4 days post amputation / 0 day post amputation


The lobe fin regeneration proteome

This study, for the first time, obtained a quantitative proteome profile of lobed fin (endoskeleton) regeneration in bichirs. The growth curve in P. senegalus (Fig. 1c) shows a difference with the growth curve in N. viridescens [20]. Iten and Bryant’s research on newt shows that amputation location information does not affect the developmental sequence of limb regeneration events, but it does affect the rate of growth of regenerates [20]. Although the lengths of the resections were different, their amputation positions were located at the proximal end of the ulna along the longitudinal axis. Therefore, it is still speculated that the growth rate of bichir appendage is higher in the early stage than the newt limb.

The P. senegalus limb/lobe-fin regeneration proteome contains proteins that are expressed in the cancer stem cells, nerve cells, muscle cells, and cartilage cells. As a consequence, we found cancer stem cell-associated proteins, such as mmp2, mmp9, mmp13, mmp14, paxillin, Ras, src and c-myc [21, 22], as well as genes expressed in the immune system, such as C5 [23], C9 [24], and HMGB1 [25]. Besides, numerous renowned and supposed growth regulation proteins were detected. Among those were c-met [26], FGF2 [27] and TGFBI [28]. We also identified a MARCKS-like protein (MLP) that was previously reported to induce the initial cell cycle response during limb regeneration in axolotls [29]. Furthermore, considerable proteins that were reported to be engaged in neurodegeneration mechanisms during limb regeneration were found, such as CASP3 [30], CASP7 [31], and its related proteins ERK1/2 [32], Clathrin [33], Cdk5 [34], SOD2 [35], and ApoE [36, 37]. In addition to stem cell pathways, we also identified wnt/β-Catenin signaling proteins which are reported to play an important role in regulating vertebrate limb regeneration, such as GSK-3β, mTOR and the β-Catenin [38, 39]. As demonstrated from these results, within the first 12 days, a wide variety of protein families that have different functions expressed in the epimorphosis process of the lobed fin are included in this proteome.

Blastema formation

In a recent study, through the iTRAQ method, Tang, J., et al. found that differentially expressed proteins were associated with wound healing, immune response, cellular process, metabolism and binding in the salamanders’ blastema formation [40]. By similar iTRAQ method, Geng et al. found that in the early stages of newt’s blastema formation, differentially expressed proteins were concentrated in the categories such as signaling, Ca2+ binding and translocation, transcription and translation, immune response, cell death, cytoskeleton, metabolism [41]. Correspondingly, in the current study, DEPs are mainly concentrated in ribosome and RNA transport, metabolic, ECM-receptor interaction, Golgi and endoplasmic reticulum, DNA replication, and Regulation of actin cytoskeleton, suggesting that these biological processes are active during the appendage regeneration process.

Integrin signaling plays an indispensable role in the regeneration process of the spinal cord [42]. In the current study, we found that the integrins-related signaling pathway was significantly enriched in bichirs lobe-fin regeneration. It was demonstrated that αvβ5 integrin plays a key role in the dedifferentiation of chondrocytes by activating ERK signaling [43]. Previous studies have shown that integrin is involved in spontaneous axonal regeneration after peripheral nerve injury, and axonal regeneration would be inhibited if integrins were inactivated by axon-repulsive molecules [43,44,45]. It is concluded that the regeneration process in vivo is always associated with increased integrin [42]. However, our study found that during the blastema formation of bichirs, two integrins were up-regulated, but one integrin was down-regulated. This means that integrins may play a complex role in the process of appendage regeneration.

Differentially regulated proteins during limb/lobe-fin regeneration

Using our iTRAQ method, hundreds of genes that were differentially regulated in various periods of limb/lobe-fin regeneration were identified. For instance, MMP13a and VTN were strongly upregulated during the first four days. MMP13 has been implicated as a contributor to skeletal muscle regeneration and is critical for myoblast migration [46], and VTN is a regulator of multimerization and collagen binding during liver regeneration [47], as well as a growth factor complex for wound repair and tissue regeneration [48]. It is possible that MMP13a and VTN also play an important role in the regulation of lobed fin regeneration in bichirs. We also found many other kinds of endopeptidases, showing that the remodeling of the extracellular matrix is essential for regeneration.

Another gene from the same cluster, Chit1, which is a biochemical marker of macrophage activation [49], was strongly upregulated during the first four days. Macrophages remove tissue debris and activate stem cell populations [50] and are required for adult salamander limb regeneration [51].

The decorin protein gene DCN was downregulated during the first four days but was later upregulated in the days that followed, and this expression pattern is consistent with this gene playing two different roles during regeneration. During the early stages, DCN may be involved in the response to dedifferentiation because a low level of decorin is conducive to the growth of various tumor cell lines and to an increased abundance of anti-inflammatory molecules [52]. During the later stages, DCN may be involved in the response to redifferentiation and growth because a high level of decorin promotes muscle cell differentiation and muscle regeneration [53], suppresses scar formation and promotes axon growth [54]. Finally, it is possible that the level of expression of decorin play a key role in lobe-fin/limb regeneration.

HSPB1 (HSP27) is a protein that is related to regeneration, which could speeds up axonal growth in vitro after peripheral nerve injury [55]. The HSPB1 protein was downregulated during the first four days but was later upregulated in the days that followed. During the early stages, the low level of HSPB1 may have occurred in response to apoptosis because HSPB1 is involved in protection against necrotic and apoptotic cell death. The high level of HSPB1 during the later stages was consistent with its late-expression property, which is different from most axonal injury-regulated and growth-associated genes [56].

Integrin beta 4 (ITGB4) is also called CD104. In the biology of invasive carcinoma, it is likely to play a crucial role. Integrins adjust cell-extracellular matrix (ECM) or cell-cell interactions, and transformed signals that manage the growth of cell and expression of gene. The forced activation of integrins can overcome inhibition and increase axon regeneration [57]. Integrin-beta 1 (ITGB1) regulates chondrocyte proliferation and apoptosis [58], and the knockdown and knockout of ITGB1 in hepatocytes impair liver regeneration [59]. The ITGB4 protein was upregulated during the first four days and later returned to normal levels in the days that followed. During the early stages, ITGB4 may be involved in the response to dedifferentiation and cell migration because the upregulation of ITGB4 can restore the regenerative performance of adult neurons [60], and promote epithelial-mesenchymal transition, cell scattering, cell motility, and vimentin expression [61].

Membrane transporters, synaptic transmission, and regulators of nerve development were found in the considerable enrichment, supporting a model in which nerve signaling might be crucial for early stages of limb regeneration. As shown from recent researches, the regeneration polarity in planarians can be regulated by the signaling of ventral nerve cords [62]. And limb regeneration in axolotl can be enhanced by nerve signaling [63, 64].

In summary, as revealed from our time-course analysis, the limb/lobe-fin regeneration is regulated by a dynamic, complicated proteomic network that is rapidly caused after amputation (Fig. 5), which might be triggered by nerve deviation and ECM interaction and transducted by integrins that launch a cascade of regenerative processes, such as cell migration, dedifferentiation and regrowth, as well as scar suppression.

Fig. 5
figure 5

Proposed model depicting the molecular mechanism based on the P. senegalus proteome. The combination of all of these mechanisms regulates the expression of hundreds of proteins and promotes cell migration, dedifferentiation and regrowth, and scar suppression

Validation of the proteome result

To validate the results of the proteome obtained from the current study, we used the Multiple reaction monitoring (MRM) approach. Multiple reaction monitoring (MRM) using mass spectrometry is a highly sensitive and selective method for the targeted quantitation of protein/peptide abundances in complex biological samples. Previous studies have shown that targeted MRM proteomics can be used as a verification tool for candidates in the context of a comprehensive proteome [65, 66]. Among the 30 selected target proteins, the expression trend of 27 proteins was verified, and the difference in numerical values may be attributed to differences in statistical models or biological variation between the samples. Three candidate proteins, such as MMP14, ITGB1, and KCTD12, show opposite expression trends in iTRAQ and MRM results, which may be due to biological variation between the samples.


In this study, we constructed a lobed fin regeneration proteome for P. senegalus and generated time-resolved profiles of the expressed genes based on the growth curve. Functional analysis indicated that the differentially expressed proteins were associated with Ribosome and RNA transport, metabolic, ECM-receptor interaction, Golgi and endoplasmic reticulum, DNA replication, and Regulation of actin cytoskeleton. These results not only provide a comprehensive dataset on differentially expressed proteins during the early stages of lobe-fin/limb regeneration but also advance our understanding of the molecular mechanisms underlying lobe-fin/limb regeneration.


Tissue collection and image acquisition

The animal experiment for this study was approved by the Institutional Animal Care and Use Committee of Institute of Hydrobiology (Approval ID: Y21304501). Adult specimens of P. senegalus were commercially obtained and maintained in glass tanks at 26 °C. The bichirs were anesthetized in 0.1% Tricaine mesylate (MS-222) before image collection and sampling. The bichirs were then washed with pure water and dried with filter paper to clean the pectoral fins, and amputations were performed with small surgical scissors and forceps along the cross section of the middle of the base of the pectoral lobe fin. For iTRAQ, tissues (or regenerates) were collected at 0, 4 and 12 days post-amputation (dpa), and these experiments were conducted in two biological replicates (distinguished by _C and _E), with each replicate consists of 10-pooled biological tissues for each stage. The 0 dpa sample was from the amputated fin itself, and divided into the 0 dpa lobe-fin (0D_L) sample and the 0 dpa ray-fin (0D_R) sample. The 4 dpa sample which consists of blastema was named 4D_B. The 12 dpa sample which consists of re-differentiation tissues was named 12D_D. The samples were flash frozen by liquid nitrogen for protein preparation. For MRM, the same method was used to collect tissues (or regenerates), with three replicates.

The regeneration progress was observed to develop the lobed fins regeneration growth curve. A vernier caliper and a camera were used to measure the length of the regenerates and to take photos, respectively, over 15 weeks. For histological analysis, stump tissues of 0 dpa, 4 dpa, and 12 dpa were collected and immersed in the 4% paraformaldehyde fixative.

iTRAQ process and MRM process

iTRAQ process—The tissues were treated using the same method as previously described [67]. Briefly, the tissues were split and then the proteins were reduced. Then, Bradford method was used to determine the protein concentration. Total protein was digested and processed with 8-plex iTRAQ reagent (Applied Biosystems) as previously described with minor modifications [67]. Briefly, the samples were labeled with the iTRAQ tags as follows: Sample 0D_R_C (117 tag), Sample 0D_R_E (113 tag), Sample 0D_L_E (114 tag), Sample 0D_L_C (118 tag), Sample 4D_B_E (115 tag), Sample 4D_B_C (119 tag), Sample 12D_D_E (116 tag), and Sample 12D_D_C (121 tag).

To separate the iTRAQ-labeled peptide mixtures, load sample and acquire data, the same method was used as previously described [67]. Briefly, peptide supernatant was loaded into an LC-20 AD nano HPLC and then eluted onto a 10 cm analytical C18 column. Data acquisition was performed with a Triple TOF 5600 System fitted with a Nano spray III source and a pulled quartz tip as the emitter. The raw data files acquired from 5600 were converted into MGF files using 5600-ms converter, and the MGF files were searched. Protein identification was performed using the Mascot search engine (Matrix Science, London, UK; version 2.3.02) against the P. senegalus genome database containing 20,582 sequences.

Proteins were identified and quantified as previously described with minor modifications [68]. Briefly, peptides with significance scores (≥20) at the 99% confidence interval by a Mascot probability analysis greater than “identity” were counted as identified. Only p values < 0.05, fold changes larger than 1.4 were considered significant differential expressions of proteins.

MRM process—For protein digestion, the same method was used as in iTRAQ process. Then, the samples were spiked with 50 fmol of β-galactosidase for data normalization, and MRM analyses were performed on a QTRAP 5500 mass spectrometer (AB SCIEX, Foster City, CA) equipped with a Waters Nano-Acquity Ultra Performance LC system. The method was as previously described with minor modifications [69]. Briefly, the peptides were separated on a BEH130 C18 column (0.075 × 200-mm column, 1.7 μm; Waters) at 300 nL/min and eluted with a gradient of 2–40% solvent B for 30 min, 40–60% solvent B for 3 min, a 2-min linear gradient to 80% solvent B and maintenance at 80% for 5 min. For the QTRAP 5500 MS, a spray voltage of 2100 V, nebulizer gas at 20 psi, and dwell time of 10 ms were used. We use MSstats with the linear mixed-effects model the P values were adjusted to control the FDR at a cutoff of 0.05. After data analysis with Skyline software, all of the proteins with a p value < 0.05 and a fold change > 1.5 were considered significant.

Bioinformatics analysis

Functional annotations of the proteins were conducted using Blast2GO [70] program against the non-redundant protein database (NR; NCBI; v2.5; animal; 8,428,593 sequences). Orthologous proteins were clustered based on the Cluster of Orthologous Groups (COGs) of proteins database [71]. KEGG (animal) [72] was used to identify the molecular interaction and reaction networks of these proteins. Hypergeometric tests were employed to perform GO enrichment and KEGG pathway enrichment, with reference dataset of all identified proteins in the whole proteome. The protein interaction network analysis was conducted using STRING ( Then the network was visualized by Cytoscape v3.2.1 and further analyzed for densely connected regions by Molecular Complex Detection (MCODE) v1.4.1 [73].

Availability of data and materials

All raw data of proteome has been deposited in the Peptide Atlas database ( with the accession PASS00842.



0 dpa lobe-fin sample


0 dpa ray-fin sample


12 dpa sample which consist of re-differentiation tissues


4 dpa sample which consist of blastema


Apolipoprotein E


Complement component 5


Complement component 9




Cyclin dependent kinase 5




Cluster of Orthologous Groups


Coefficient of variation




Differentially expressed proteins


days post amputation


Extracellular matrix


Extracellular signal–regulated kinases


Fibroblast growth factor


Gene ontology


Glycogen synthase kinase 3 beta


High mobility group box


Heat shock protein beta-1


Integrin Subunit Beta 4


isobaric tag for relative and absolute quantitation


Kyoto Encyclopedia of Genes and Genomes


MARCKS-like protein


matrix metallopeptidase 13a


matrix metallopeptidase 14


matrix metallopeptidase 2


matrix metallopeptidase 9


Multiple-reaction monitoring


Mass spectrometry


mammalian target of rapamycin


Million years ago




Superoxide dismutase 2


Transforming Growth Factor Beta Induced




  1. Brockes JP, Kumar A. Appendage regeneration in adult vertebrates and implications for regenerative medicine. Science. 2005;310:1919–23.

    Article  CAS  PubMed  Google Scholar 

  2. Galis F, Wagner GP, Jockusch EL. Why is limb regeneration possible in amphibians but not in reptiles, birds, and mammals? Evol Dev. 2003;5:208–20.

    Article  PubMed  Google Scholar 

  3. Han M, Yang X, Lee J, Allan CH, Muneoka K. Development and regeneration of the neonatal digit tip in mice. Dev Biol. 2008;315:125–35.

    Article  CAS  PubMed  Google Scholar 

  4. Cuervo R, Hernandez-Martinez R, Chimal-Monroy J, Merchant-Larios H, Covarrubias L. Full regeneration of the tribasal Polypterus fin. Proc Natl Acad Sci U S A. 2012;109:3838–43.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Nogueira AF, Costa CM, Lorena J, Moreira RN, Frota-Lima GN, Furtado C, Robinson M, Amemiya CT, Darnet S, Schneider I. Tetrapod limb and sarcopterygian fin regeneration share a core genetic programme. Nat Commun. 2016;7:13364.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Conant EB. Regeneration in the African lungfish, Protopterus. 3. Regeneration during fasting and estivation. Biol Bull. 1973;144:248–61.

    Article  CAS  PubMed  Google Scholar 

  7. Nikiforova AI, Golichenkov VA. Characteristics of the reparative regeneration of fins in the polypterid fish (Polypteridae, Actinopterygii). Ontogenez. 2012;43:136–42.

    CAS  PubMed  Google Scholar 

  8. Betancur-R R, Broughton RE, Wiley EO, Carpenter K, López JA, Li C, Holcroft NI, Arcila D, Sanciangco M, Cureton Ii JC, et al. The tree of life and a new classification of bony fishes. PLoS Currents. 2013;5 ecurrents.tol.53ba26640df26640ccaee26675bb26165c26648c26288.

  9. Mank JE, Avise JC. Comparative phylogenetic analysis of male alternative reproductive tactics in ray-finned fishes. Evolution. 2006;60:1311–6.

    Article  PubMed  Google Scholar 

  10. Broughton RE, Betancur-R R, Li C, Arratia G, Ortí G. Multi-locus phylogenetic analysis reveals the pattern and tempo of bony fish evolution. PLoS Currents. 2013;5 ecurrents.tol.2ca8041495ffafd8041490c8092756e75247483e.

  11. Bartsch P, Gemballa S, Piotrowski T. The embryonic and larval development of Polypterus senegalus Cuvier, 1829: its staging with reference to external and skeletal features, behaviour and Locomotory habits. Acta Zool. 1997;78:309–28.

    Article  Google Scholar 

  12. Clack JA. Earliest known tetrapod braincase and the evolution of the stapes and fenestra ovalis. Nature. 1994;369:392–4.

    Article  Google Scholar 

  13. Coates MI, Clack JA. Fish-like gills and breathing in the earliest known tetrapod. Nature. 1991;352:234–6.

    Article  Google Scholar 

  14. Budgett JS. On the structure of the larval Polypterus. Trans Zool Soc London. 1902;16:315–46.

    Article  Google Scholar 

  15. Standen EM, Du TY, Larsson HC. Developmental plasticity and the origin of tetrapods. Nature. 2014;513:54–8.

    Article  CAS  PubMed  Google Scholar 

  16. Klaatsch H. Die Brustflosse der Crossopterygier: Ein Beitrag zur Anwendung der Archipterygium-Theorie auf die Gliedmaassen der Landwirbelthiere. Leipzig; 1896.

  17. Wilhelm BC, Du TY, Standen EM, Larsson HC. Polypterus and the evolution of fish pectoral musculature. J Anat. 2015;226:511–22.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Smith JJ, Putta S, Zhu W, Pao GM, Verma IM, Hunter T, Bryant SV, Gardiner DM, Harkins TT, Voss SR. Genic regions of a large salamander genome contain long introns and novel genes. BMC Genomics. 2009;10:19.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Lu S, Yang L, Jiang H, Chen J, Yu G, Chen Z, Xia X, He S. Bichirs employ similar genetic pathways for limb regeneration as are used in lungfish and salamanders. Gene. 2018.

  20. Iten LE, Bryant SV. Forelimb regeneration from different levels of amputation in the newt, Notophthalmus viridescens: length, rate, and stages. Wilhelm Roux’ Archiv für Entwicklungsmechanik der Organismen. 1973;173:263–82.

    Article  PubMed  Google Scholar 

  21. Roomi MW, Monterrey JC, Kalinovsky T, Rath M, Niedzwiecki A. Patterns of MMP-2 and MMP-9 expression in human cancer cell lines. Oncol Rep. 2009;21:1323–33.

    CAS  PubMed  Google Scholar 

  22. Rha SY, Yang WI, Kim JH, Roh JK, Min JS, Lee KS, Kim BS, Chung HC. Different expression patterns of MMP-2 and MMP-9 in breast cancer. Oncol Rep. 1998;5:875–9.

    CAS  PubMed  Google Scholar 

  23. Godwin JW, Brockes JP. Regeneration, tissue injury and the immune response. J Anat. 2006;209:423–32.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Rutkowski MJ, Sughrue ME, Kane AJ, Ahn BJ, Fang SN, Parsa AT. The complement cascade as a mediator of tissue growth and regeneration. Inflamm Res. 2010;59:897–905.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Sattler S, Rosenthal N. The neonate versus adult mammalian immune system in cardiac repair and regeneration. Biochim Biophys Acta. 1863;2016:1813–21.

    Google Scholar 

  26. Huh CG, Factor VM, Sanchez A, Uchida K, Conner EA, Thorgeirsson SS. Hepatocyte growth factor/c-met signaling pathway is required for efficient liver regeneration and repair. Proc Natl Acad Sci U S A. 2004;101:4477–82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Giampaoli S, Bucci S, Ragghianti M, Mancino G, Zhang F, Ferretti P. Expression of FGF2 in the limb blastema of two Salamandridae correlates with their regenerative capability. Proc Biol Sci. 2003;270:2197–205.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Lenkowski JR, Qin Z, Sifuentes CJ, Thummel R, Soto CM, Moens CB, Raymond PA. Retinal regeneration in adult zebrafish requires regulation of TGFbeta signaling. Glia. 2013;61:1687–97.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Sugiura T, Wang H, Barsacchi R, Simon A, Tanaka EM. MARCKS-like protein is an initiating molecule in axolotl appendage regeneration. Nature. 2016;531:237–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Young C, Roth KA, Klocke BJ, West T, Holtzman DM, Labruyere J, Qin YQ, Dikranian K, Olney JW. Role of caspase-3 in ethanol-induced developmental neurodegeneration. Neurobiol Dis. 2005;20:608–14.

    Article  CAS  PubMed  Google Scholar 

  31. Gillingwater TH, Wishart TM. Mechanisms underlying synaptic vulnerability and degeneration in neurodegenerative disease. Neuropathol Appl Neurobiol. 2013;39:320–34.

    Article  CAS  PubMed  Google Scholar 

  32. Subramaniam S, Unsicker K. ERK and cell death: ERK1/2 in neuronal death. FEBS J. 2010;277:22–9.

    Article  CAS  PubMed  Google Scholar 

  33. Rudinskiy N, Grishchuk Y, Vaslin A, Puyal J, Delacourte A, Hirling H, Clarke PG, Luthi-Carter R. Calpain hydrolysis of alpha- and beta2-adaptins decreases clathrin-dependent endocytosis and may promote neurodegeneration. J Biol Chem. 2009;284:12447–58.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Cheung ZH, Ip NY. Cdk5: a multifaceted kinase in neurodegenerative diseases. Trends Cell Biol. 2012;22:169–75.

    Article  CAS  PubMed  Google Scholar 

  35. Flynn JM, Melov S. SOD2 in mitochondrial dysfunction and neurodegeneration. Free Radic Biol Med. 2013;62:4–12.

    Article  CAS  PubMed  Google Scholar 

  36. Masliah E, Samuel W, Veinbergs I, Mallory M, Mante M, Saitoh T. Neurodegeneration and cognitive impairment in apoE-deficient mice is ameliorated by infusion of recombinant apoE. Brain Res. 1997;751:307–14.

    Article  CAS  PubMed  Google Scholar 

  37. Krasemann S, Madore C, Cialic R, Baufeld C, Calcagno N, El Fatimy R, Beckers L, O'Loughlin E, Xu Y, Fanek Z, et al. The TREM2-APOE pathway drives the transcriptional phenotype of dysfunctional microglia in neurodegenerative diseases. Immunity. 2017;47:566–81 e569.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Kawakami Y, Esteban CR, Raya M, Kawakami H, Marti M, Dubova I, Belmonte JCI. Wnt/beta-catenin signaling regulates vertebrate limb regeneration. Genes Dev. 2006;20:3232–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Arioka M, Takahashi-Yanaga F, Sasaki M, Yoshihara T, Morimoto S, Takashima A, Mori Y, Sasaguri T. Acceleration of bone development and regeneration through the Wnt/beta-catenin signaling pathway in mice heterozygously deficient for GSK-3beta. Biochem Biophys Res Commun. 2013;440:677–82.

    Article  CAS  PubMed  Google Scholar 

  40. Tang J, Yu Y, Zheng H, Yin L, Sun M, Wang W, Cui J, Liu W, Xie X, Chen F. ITRAQ-based quantitative proteomic analysis of Cynops orientalis limb regeneration. BMC Genomics. 2017;18:750.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  41. Geng XF, Guo JL, Zang XY, Sun JY, Li PF, Zhang FC, Xu CS. iTRAQ-based proteomic analysis of adaptive response in the regenerating limb of the Cynops orientalis newt. Int J Dev Biol. 2015;59:487–96.

    Article  CAS  PubMed  Google Scholar 

  42. Lemons ML, Condic ML. Integrin signaling is integral to regeneration. Exp Neurol. 2008;209:343–52.

    Article  CAS  PubMed  Google Scholar 

  43. Fukui N, Ikeda Y, Tanaka N, Wake M, Yamaguchi T, Mitomi H, Ishida S, Furukawa H, Hamada Y, Miyamoto Y, et al. alphavbeta5 integrin promotes dedifferentiation of monolayer-cultured articular chondrocytes. Arthritis Rheum. 2011;63:1938–49.

    Article  CAS  PubMed  Google Scholar 

  44. Huo Y, Yin XL, Ji SX, Zou H, Lang M, Zheng Z, Cai XF, Liu W, Chen CL, Zhou YG, et al. Amino-Nogo inhibits optic nerve regeneration and functional recovery via the integrin alphav signaling pathway in rats. Cell Physiol Biochem. 2015;35:616–26.

    Article  CAS  PubMed  Google Scholar 

  45. Hu F, Strittmatter SM. The N-terminal domain of Nogo-a inhibits cell adhesion and axonal outgrowth by an integrin-specific mechanism. J Neurosci. 2008;28:1262–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Lei H, Leong D, Smith LR, Barton ER. Matrix metalloproteinase 13 is a new contributor to skeletal muscle regeneration and critical for myoblast migration. Am J Physiol Cell Physiol. 2013;305:C529–38.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Sano K, Asanuma-Date K, Arisaka F, Hattori S, Ogawa H. Changes in glycosylation of vitronectin modulate multimerization and collagen binding during liver regeneration. Glycobiology. 2007;17:784–94.

    Article  CAS  PubMed  Google Scholar 

  48. Upton Z, Cuttle L, Noble A, Kempf M, Topping G, Malda J, Xie Y, Mill J, Harkin DG, Kravchuk O, et al. Vitronectin: growth factor complexes hold potential as a wound therapy approach. J Invest Dermatol. 2008;128:1535–44.

    Article  CAS  PubMed  Google Scholar 

  49. Di Rosa M, Malaguarnera G, De Gregorio C, Drago F, Malaguarnera L. Evaluation of CHI3L-1 and CHIT-1 expression in differentiated and polarized macrophages. Inflammation. 2013;36:482–92.

    Article  CAS  PubMed  Google Scholar 

  50. Arnold L, Henry A, Poron F, Baba-Amer Y, van Rooijen N, Plonquet A, Gherardi RK, Chazaud B. Inflammatory monocytes recruited after skeletal muscle injury switch into antiinflammatory macrophages to support myogenesis. J Exp Med. 2007;204:1057–69.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Godwin JW, Pinto AR, Rosenthal NA. Macrophages are required for adult salamander limb regeneration. Proc Natl Acad Sci U S A. 2013;110:9415–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Moscatello DK, Santra M, Mann DM, McQuillan DJ, Wong AJ, Iozzo RV. Decorin suppresses tumor cell growth by activating the epidermal growth factor receptor. J Clin Invest. 1998;101:406–12.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Li Y, Li J, Zhu J, Sun B, Branca M, Tang Y, Foster W, Xiao X, Huard J. Decorin gene transfer promotes muscle cell differentiation and muscle regeneration. Mol Ther. 2007;15:1616–22.

    Article  CAS  PubMed  Google Scholar 

  54. Davies JE, Tang X, Denning JW, Archibald SJ, Davies SJ. Decorin suppresses neurocan, brevican, phosphacan and NG2 expression and promotes axon growth across adult rat spinal cord injuries. Eur J Neurosci. 2004;19:1226–42.

    Article  PubMed  Google Scholar 

  55. Ma CH, Omura T, Cobos EJ, Latremoliere A, Ghasemlou N, Brenner GJ, van Veen E, Barrett L, Sawada T, Gao F, et al. Accelerating axonal growth promotes motor recovery after peripheral nerve injury in mice. J Clin Invest. 2011;121:4332–47.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Costigan M, Mannion RJ, Kendall G, Lewis SE, Campagna JA, Coggeshall RE, Meridith-Middleton J, Tate S, Woolf CJ. Heat shock protein 27: developmental regulation and expression after peripheral nerve injury. J Neurosci. 1998;18:5891–900.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Eva R, Fawcett J. Integrin signalling and traffic during axon growth and regeneration. Curr Opin Neurobiol. 2014;27:179–85.

    Article  CAS  PubMed  Google Scholar 

  58. Zhang LQ, Zhao GZ, Xu XY, Fang J, Chen JM, Li JW, Gao XJ, Hao LJ, Chen YZ. Integrin-beta1 regulates chondrocyte proliferation and apoptosis through the upregulation of GIT1 expression. Int J Mol Med. 2015;35:1074–80.

    Article  CAS  PubMed  Google Scholar 

  59. Speicher T, Siegenthaler B, Bogorad RL, Ruppert R, Petzold T, Padrissa-Altes S, Bachofner M, Anderson DG, Koteliansky V, Fassler R, Werner S. Knockdown and knockout of beta1-integrin in hepatocytes impairs liver regeneration through inhibition of growth factor signalling. Nat Commun. 2014;5:3862.

    Article  CAS  PubMed  Google Scholar 

  60. Condic ML. Adult neuronal regeneration induced by transgenic integrin expression. J Neurosci. 2001;21:4782–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Masugi Y, Yamazaki K, Emoto K, Effendi K, Tsujikawa H, Kitago M, Itano O, Kitagawa Y, Sakamoto M. Upregulation of integrin beta4 promotes epithelial-mesenchymal transition and is a novel prognostic marker in pancreatic ductal adenocarcinoma. Lab Investig. 2015;95:308–19.

    Article  CAS  PubMed  Google Scholar 

  62. Oviedo NJ, Morokuma J, Walentek P, Kema IP, Gu MB, Ahn JM, Hwang JS, Gojobori T, Levin M. Long-range neural and gap junction protein-mediated cues control polarity during planarian regeneration. Dev Biol. 2010;339:188–99.

    Article  CAS  PubMed  Google Scholar 

  63. Endo T, Bryant SV, Gardiner DM. A stepwise model system for limb regeneration. Dev Biol. 2004;270:135–45.

    Article  CAS  PubMed  Google Scholar 

  64. Satoh A, Gardiner DM, Bryant SV, Endo T. Nerve-induced ectopic limb blastemas in the axolotl are equivalent to amputation-induced blastemas. Dev Biol. 2007;312:231–44.

    Article  CAS  PubMed  Google Scholar 

  65. Prasad B. Targeted MRM proteomics is a better protein quantification method over Western-blotting. J Anal Bio Tech. 2013.

  66. Gillette MA, Carr SA. Quantitative analysis of peptides and proteins in biomedicine by targeted mass spectrometry. Nat Methods. 2013;10:28–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  67. Xu CP, Li X, Hu YJ, Cui Z, Wang L, Liang L, Zhou YL, Yang YJ, Yu B. Quantitative proteomics reveals ELP2 as a regulator to the inhibitory effect of TNF-alpha on osteoblast differentiation. J Proteome. 2015;114:234–46.

    Article  CAS  Google Scholar 

  68. Yan G, Lestari R, Long B, Fan Q, Wang Z, Guo X, Yu J, Hu J, Yang X, Chen C, et al. Comparative proteomics analysis reveals L-arginine activates ethanol degradation pathways in HepG2 cells. Sci Rep. 2016;6:23340.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  69. Su Y, Xu L, Wang Z, Peng Q, Yang Y, Chen Y, Que Y. Comparative proteomics reveals that central metabolism changes are associated with resistance against Sporisorium scitamineum in sugarcane. BMC Genomics. 2016;17:800.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  70. Conesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21:3674–6.

    Article  CAS  PubMed  Google Scholar 

  71. Tatusov RL, Galperin MY, Natale DA, Koonin EV. The COG database: a tool for genome-scale analysis of protein functions and evolution. Nucleic Acids Res. 2000;28:33–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  72. Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  73. Bader GD, Hogue CW. An automated method for finding molecular complexes in large protein interaction networks. BMC Bioinformatics. 2003;4:2.

    Article  PubMed  PubMed Central  Google Scholar 

Download references


We thank Dr. David M. Gardiner and Dr. Susan V. Bryant for helpful discussions.


This work was supported by grants from the National Natural Science Foundation of China (91731301) and from the Strategic Priority Research Program (XDB13020100 and XDB060101) to S.H., and from the National Natural Science Foundation of China (31972866 and 31601858) and Institute of Hydrobiology, Chinese Academy of Sciences (Y55Z09 and Y85E03) to L.Y., and from the National Natural Science Foundation of China (31702016) to Y.W..

Author information

Authors and Affiliations



SL, LY, FG and SH designed the experiments. SL, YW, XG and XW performed the experiments. SL, LY, KD and QX analyzed the data. SL wrote the paper. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Shunping He.

Ethics declarations

Ethics approval

All experiments were carried out in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1.

All of the proteins identified in the present study with annotations.

Additional file 2.

Overview of the proteomic results. (A) Basic information statistics. (B) Identified peptide distribution. (C) Protein coverage. (D) Mass errors were determined for all identified peptides.

Additional file 3.

Differentially expressed proteins at different time points.

Additional file 4.

Congruence of evaluation among different data set. (A) CV distribution. a, 0D_L-VS-4D_B. b, 0D_L-VS-12D_D. (B) Protein ratio distribution. c, 0D_L-VS-4D_B. d, 0D_L-VS-12D_D. (C) Linear regression analysis based on the log2-transformed protein ratios. e, linear correlation between the two experiment replicates for 4D_B. f, linear correlation between the two experiment replicates for 12D_D.

Additional file 5.

GO and pathway enrichment analysis of differentially expressed proteins.

Additional file 6.

Protein-protein interaction networks of DEPs in 4D/0D.

Additional file 7.

Some protein-protein interaction clusters that calculated by the MCODE software.

Additional file 8.

Details of the dynamic MRM results.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Lu, S., Xiong, Q., Du, K. et al. Comparative iTRAQ proteomics revealed proteins associated with lobed fin regeneration in Bichirs. Proteome Sci 17, 6 (2019).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: