Comparative gene expression profiling between optic nerve and spinal cord injury in Xenopus laevis reveals a core set of genes inherent in successful regeneration of vertebrate central nervous system axons
BMC Genomics volume 21, Article number: 540 (2020)
The South African claw-toed frog, Xenopus laevis, is uniquely suited for studying differences between regenerative and non-regenerative responses to CNS injury within the same organism, because some CNS neurons (e.g., retinal ganglion cells after optic nerve crush (ONC)) regenerate axons throughout life, whereas others (e.g., hindbrain neurons after spinal cord injury (SCI)) lose this capacity as tadpoles metamorphose into frogs. Tissues from these CNS regions (frog ONC eye, tadpole SCI hindbrain, frog SCI hindbrain) were used in a three-way RNA-seq study of axotomized CNS axons to identify potential core gene expression programs for successful CNS axon regeneration.
Despite tissue-specific changes in expression dominating the injury responses of each tissue, injury-induced changes in gene expression were nonetheless shared between the two axon-regenerative CNS regions that were not shared with the non-regenerative region. These included similar temporal patterns of gene expression and over 300 injury-responsive genes. Many of these genes and their associated cellular functions had previously been associated with injury responses of multiple tissues, both neural and non-neural, from different species, thereby demonstrating deep phylogenetically conserved commonalities between successful CNS axon regeneration and tissue regeneration in general. Further analyses implicated the KEGG adipocytokine signaling pathway, which links leptin with metabolic and gene regulatory pathways, and a novel gene regulatory network with genes regulating chromatin accessibility at its core, as important hubs in the larger network of injury response genes involved in successful CNS axon regeneration.
This study identifies deep, phylogenetically conserved commonalities between CNS axon regeneration and other examples of successful tissue regeneration and provides new targets for studying the molecular underpinnings of successful CNS axon regeneration, as well as a guide for distinguishing pro-regenerative injury-induced changes in gene expression from detrimental ones in mammals.
The capacity for regeneration is an ancient characteristic of metazoa. Although many invertebrates can reconstitute an entire body from residual tissue fragments [7, 60], vertebrates regenerate only tissues and organs. This capacity becomes more restricted both phylogenetically in the progression from anamniote to amniote and developmentally during the transition from embryo to adult (Fig. 1a). For example, most vertebrates are able to regenerate peripheral nervous system (PNS) axons throughout life and central nervous system (CNS) axons as embryos, but only anamniotes can also regenerate CNS axons as adults. Anuran amphibians, such as the frog Xenopus laevis, occupy a transition point in this progression . Like other anamniotes, Xenopus regenerates optic axons sufficiently to restore vision throughout life [12, 39, 135], but like amniotes, it loses its capacity to functionally regenerate spinal cord axons developmentally, because of the surge of thyroid hormone that drives metamorphosis in frogs  and late fetal development in mammals . Thus, Xenopus offers a unique opportunity for exploring within the same organism why some regions of the CNS lose their ability to regenerate axons in development, whereas others retain it. In this study, we exploited these features of Xenopus in a novel three-way comparison of RNA-seq data from two regions of the CNS that regenerate axons, one throughout life (juvenile frog eye after optic nerve crush (ONC) [39, 127, 135]) and one that does so only transiently, before and after this transition (tadpole vs. juvenile frog hindbrain after spinal cord injury (SCI) [11, 43]), to identify potential core features that distinguish CNS axon-regenerative responses from failed ones.
With the exception of an mRNA-seq study in lamprey, which demonstrated that the response of brain and spinal cord to SCI indeed differ , previous genome-wide expression studies on spinal cord regeneration have focused on either the spinal cord itself, with the lesion at the epicenter (Xenopus laevis tadpole and frog , turtle , zebrafish ), or on a regenerating tail (e.g., salamander , Xenopus tropicalis tadpole [21, 55, 80, 102]). Such regeneration involves not only axon regeneration, but also considerable wound repair and tissue restoration. Xenopus laevis hindbrain, which is the principal source of descending axons that project from brain to spinal cord [43, 143], provides an opportunity to compare transcriptional responses to SCI between axon-regenerative and non-regenerative states of CNS neurons and their support cells, separately from the lesion site and its attendant wound repair and tissue restoration. Moreover, hindbrain neurons that regenerate descending axons after SCI (i.e., reticular and raphe nuclei) in the tadpole are the very same cells that fail to do so in the frog . The Xenopus eye after ONC provides a similar opportunity to study an injury response, separately from the injury site (Fig. 1d), but after metamorphosis. To date, genome wide studies in the regenerating Xenopus visual system have been limited to SDS PAGE [126, 128] and TRAP (RNA-Seq after Translating Ribosome Affinity Purification of mRNAs) studies . Thus, both systems would benefit from baseline RNA-Seq studies of their responses to axotomy.
In performing such an RNA-Seq study, we found that the injury responses of the two regenerative tissues (frog eye after ONC and tadpole hindbrain after SCI) shared features in common with each other not seen in the non-regenerative tissue (frog hindbrain after SCI). Injury-induced genes uniquely shared between the two regenerative tissues included many previously implicated in promoting cell survival and regeneration of neural and non-neural tissues alike, and revealed cellular and physiological processes potentially at the core of a successful regenerative response to CNS injury. We offer these results and the insights gained from their analyses as a resource to stimulate new avenues of investigation into the molecular underpinnings of successful recovery from CNS traumatic injury.
Experimental design – animal procedures, sample and data collection, and Bioinformatic analyses
Studies documenting regeneration of CNS axons with return of function after ONC in frogs and the transient nature of functional recovery from SCI injury in tadpoles go back decades [36, 125]. Although initially done in other species, subsequent studies extended these findings to Xenopus laevis [11, 39]. Electrophysiological [40, 128], anatomical [33, 128, 152, 160], and ultrastructural  studies of X. laevis frog optic axon regeneration, as well as behavioral and anatomical studies of tadpole SCI [11, 42, 43], both established the time course of recovery and demonstrated true regeneration of cut axons from surviving neurons, laying the foundation for much subsequent work (e.g., reviewed in [68, 69, 127]). From these studies, three time points can be gleaned for making appropriate comparisons between the two systems (Fig. 1c): 1) an early phase of recovery at three days, when histological studies observed regenerating axons in both systems first penetrate the lesion [11, 160]; 2) the peak phase of regenerative axon outgrowth at seven days for SCI and eleven days for ONC, when the density of regenerating hindbrain axons crossing the lesion nears its peak [11, 42], and when regenerating optic axons fill the optic tract but have not yet arrived at their targets in the tectum , respectively; 3) a late recovery phase at three weeks, when regenerating optic axons cover the tectum, but have not yet fully restored a retinotopic map [128, 160], and behavioral recovery is largely but not yet complete in SCI [11, 42]. These same time points were used for the non-regenerative juvenile frog SCI to obtain data on kinetic differences.
In planning the workflow for this study (Fig. 1d), we used the entire tissue encompassing the axotomized neurons, because the cells responding to CNS injuries also included interneurons , glia , macrophages and other myeloid cells [45, 152], as well as the axotomized neurons themselves. Analyzing complex tissues requires efforts to reduce technical variation as much as possible. By using animals of the same strain and from the same supplier, preparing RNA and cDNA libraries in parallel, and then mixing and sequencing all samples together on the same flow cells, we sought to minimize technical variability, which is otherwise inherent in post hoc comparisons derived from separate studies. To balance budgetary tradeoffs between read depth and the statistical power gained from multiple biological replicates, we sequenced at a nominal depth of 30 million reads per sample, each of which was derived from either five pooled hindbrains or six pooled eyes from co-reared animals, with three such pooled biological replicates for each condition at each time point . This yielded 51 samples representing 17 conditions (Fig. 1e, F3). Because SCI tadpoles were still developing, the unoperated controls for each time point were age-matched and co-reared. For juvenile frog SCI hindbrain, a single control group sufficed for all three time points. For ONC, the left, contralateral unoperated eye served as controls for each time point to further distinguish effects directly associated with axon regeneration from indirect injury responses [126, 128]. This is important because the unoperated contralateral eye also responds to a unilateral nerve crush through indirect effects on retinal circuitry [3, 41, 127]. To assay this indirect response, we included for eye an extra control group of surgically naive, unoperated animals for future studies; this group was included in the PCA but was not used otherwise for differential gene expression analysis (see below).
Fastq files containing the sorted raw sequences were aligned against the Xenopus laevis genome (Xenbase v9.1; http://www.xenbase.org/ RRID:SCR_003280) using Bowtie2 (v2.2.9)  and TopHat (v2.1.1) , and annotated using the Mayball gene model [24, 106, 108]. The alignments yielded 34.4 ± 3.1 (S.D.) million successfully aligned reads per sample (Fig. 1f1), with only 9.6 ± 2.9%, (S.D.) reads initially flagged as potentially duplicate alignments. Such potential duplicates can occur in Xenopus laevis due to its ancestrally (~ 30 Mya) duplicated genome . The vast majority of these potential duplicates were resolved by assigning the alignment with the higher score to separate genes on different chromosomes (referred to in X. laevis as S and L homeologs). For the small remainder (< 10% of potentially duplicate alignments), reads were randomly distributed between the two homeologs. To confirm the accuracy of this procedure, we visualized a subset directly (Integrative Genomics Viewer (IGV), v2.3.88 [115, 136]). The fastq files and gene counts are both available in Gene Expression Omnibus (GEO; https://www.ncbi.nlm.nih.gov/geo/) under accession number GSE137844. Note here that gene nomenclature in Xenopus laevis is based on prospective orthology with human genes . Due to the ancient genome duplication and long evolutionary divergence time between X. laevis and human (> 300 Mya), the L and S X. laevis homeologs may have more or fewer paralogs than their human orthologs. For the initial gene counts, we enumerated these separately, because duplicate genes can have different expression patterns, potentially leading to different functions [96, 116]. Generally speaking, L homeologs are expressed at higher levels than S homeologs , and not surprisingly, they appeared more frequently than S homeologs in our gene lists, since efficiency of gene calls tends to correlate with expression level. Functional comparisons necessarily combined homeologs and paralogs under a single gene term, since most functional studies are done in other species and we cannot a priori evaluate the significance of homeologs and paralogs for these genes. In addition, for this study, we looked only at genes that have been annotated in X. laevis, since the principal aim was to evaluate these genes with regard to what is known from other models about function. We left the analysis of fastq files for the biological significance of multiple homeologs and paralogs, as well as of novel genes lacking human orthologs and other gene features such as long non-coding RNAs, to future studies. Nonetheless, for the interested reader, we have noted in Additional_File4_DESR_Data.xlsm the expression behavior of remaining homeologs and paralogs that were not differentially expressed in our current analysis.
To generate gene lists for downstream data analysis of differentially expressed genes (Fig. 1e), we used Cufflinks/CuffDiff2 [113, 114, 137, 139, 140] and the associated utilities of CummeRbund (v2.16.0) . We made this choice primarily because CuffDiff2 uses a paired means test to estimate Type I errors (p-values) for each gene pair, independently of other genes. For each pair-wise comparison, these unadjusted p-values were then ranked among the other pair-wise comparisons to assess the probability of false discovery (q) . Because samples were complex tissues, with multiple cell types potentially responding differently to the lesion, we felt that this procedure of evaluating each gene pair independently was biologically more appropriate for our study than other programs, such as DESeq2, which use a difference of variance test that incorporates the variance of all the genes in a sample into its calculations [4, 79]. In addition, unlike DESeq2, CuffDiff2 normalizes genes for transcript length, making it possible to compare expression levels (FPKM) among different genes. Although we used the CuffDiff2 outputs for downstream analyses (Additional_File1_Differential_Expression_Analysis_by_Cuffdif.xlsm), we provide DESeq2 outputs for the interested reader (Additional_File2_Differential_Expression_Analysis_by_DESeq2.xlsm). Note that the overall trends in the data with respect to temporal patterns of differentially expressed genes and the degree of overlap between different experiments at each time point matched well between the two programs (seen by comparing the CuffDiff2 data of Fig. 2 with that of DESeq2 in Additional_File3_DESeq2_Metadata.pdf), except that DESeq2 identified fewer differentially expressed genes than CuffDiff2, with > 85% of them also called by CuffDiff2.
We used an FDR ≤ 0.05 as the criterion for selecting differentially expressed genes without regard to fold-change, since cell types responding to axotomy represent only a subset of cells in each tissue, and regulatory genes typically exhibit smaller fold-changes in expression than structural genes. We did, however, filter out genes with expressions that fell below the minimal threshold in both the experimental and paired control samples (FPKM = 0.45; Fig. 1f2 ). Under these criteria and sampling conditions, differentially expressed genes generally had |log2(fold-changes)| > 0.3, and unadjusted p values of < 0.02. The percentage of differentially expressed genes in each sample ranged from 1.4% (3 day, juvenile frog ONC eye) to 11.6% (3 day juvenile frog SCI hindbrain) of the total annotated X. laevis genome (Fig. 2b; 100% = 24,382 genes).
Identifying shared features of injury-induced gene expression between axon-regenerative vs. non-regenerative CNS
The principal goal of this study was to identify cellular processes potentially shared between the two regenerative cases that were not shared with the non-regenerative case. The first indication that such processes existed came from the temporal patterns of numbers of differentially expressed genes (Fig. 2a). The numbers of such genes in the axon-regenerative tadpole SCI hindbrain and juvenile ONC eye both reached their maximum at the peak phase of regenerative axon outgrowth (7/11 days), whereas that of the non-regenerative, juvenile frog SCI hindbrain did so at 3 days.
To assess the degree to which these shared responses involved the same genes, we examined at each time point, the overlap among the three tissues (Fig. 2c,d,e). This analysis identified 324 genes that were differentially expressed in the two axon-regenerative CNS regions but not in the non-regenerative one. We termed these genes “Differentially Expressed in Successful Regeneration” (DESR genes: 253 up-regulated; 71 down-regulated, with some appearing at multiple time points). Although the differentially expressed genes in each individual tissue had a predicted FDR of < 0.05, because DESR genes were the result of combining data from two independent experiments, their expected FDR was actually much smaller (< 0.052, which for 324 genes would be ~ 1 gene). Overall, DESR genes comprised approximately 2.7% (324/12,059) of all the injury-induced, differentially expressed genes (Fig. 2c,d,e). This relatively small fraction indicated that the vast majority of differentially expressed genes in each tissue represented tissue-specific responses to axotomy. The DESR expression data, segregated by time point and including information about related paralogs and homeologs, is provided in Additional_File4_DESR_Data.xlsm. Meta analysis conducted on DESeq2 data yielded very similar results (Additional_File3_DESeq2_Metadata.pdf), indicating that these findings were not simply an artifact of the differential expression algorithm.
We next used Principal Component Analysis to obtain an overview of the relationships among gene expression profiles of the various samples (Fig. 3). As implemented in the Tuxedo Protocol (of which CuffDiff2 is a part), PCA generates eigenvectors for each set of pooled biological replicates, with the angles between eigenvectors indicating how well gene expression profiles are correlated between conditions. Diminishing acute angles indicate more highly correlated expression profiles, and orthogonal and increasing divergent angles indicate poorly and negatively correlated profiles, respectively. When all 17 conditions were analyzed (Fig. 3a), the close clustering of the eigenvectors of the eye samples and their large divergence from the closely clustered hindbrain vectors confirmed the conclusion from the overlap analysis that tissue of origin exerted the strongest influence on expression profiles. Developmental stage, injury condition, and time elapsed since injury exerted much less influence. Separate PCA of only the SCI hindbrain expression profiles (Fig. 3b) revealed that within this tissue, regenerative capacity exerted the next strongest influence, with axon-regenerative tadpole eigenvectors clustering separately from those of the non-regenerative juvenile frog. As indicated by the angles between eigenvectors, the injury responses of the axon-regenerative tadpole hindbrains differed more between the operated samples and their age-matched, unoperated controls than was the case for the non-regenerative frog hindbrains, which were more tightly clustered. Moreover, the large divergence between the SCI eigenvectors of their respective peak time points for differential expression (7 days for tadpole and 3 days for frog) indicated that these peak responses involved mostly different genes rather than differences in the timing of differential expression among the same genes. Similarly to tadpole SCI, the eigenvectors of the operated frog ONC eye also diverged extensively from those of their paired, contralateral unoperated eye controls (Fig. 3c). Also similar to tadpole SCI, the close proximity of the 3-week operated and contralateral unoperated eigenvectors with that of eyes from unoperated animals (termed surgically naïve) further indicated that by 3 weeks, the injury response had nearly, but not quite, returned to the pre-injury state. In addition, the divergence of the eigenvectors of all the contralateral unoperated controls diverged from that of surgically naive eye, consistent with the expectation arising from biochemical studies demonstrating that the unoperated eye in a unilateral ONC also responds to the injury [3, 41, 127]. This observation further supports our decision to compare operated eye with its unoperated contralateral control to better emphasize changes directly associated with regenerating an axon, as opposed to more generalized responses associated with the traumatic disruption of visual neural circuitry . Because the bilateral nature of SCI abrogated such a distinction in hindbrain, this comparison in eye acted as an additional filter in the three-way comparison, to resolve those injury-induced differences most directly related to axon regeneration from those related more generally to trauma. Scatterplot PCAs (Additional_File5_PCA_Scatterplot.pdf) generally confirmed the eigenvector PCA analyses, while revealing additional details about the variability of individual biological replicates. For example, the known variability of the ONC injury response at its beginning [3, 41, 151] was readily visible in the PCA scatterplots at 3 days. Similarly in the scatterplots of tadpole SCI hindbrain at 3 weeks, two of three replicates of the operated animals clustered closely with their age-matched unoperated controls, suggesting that animals in these two groups had recovered more fully than had those of the third. Such variability is consistent with the large amount of variability seen from behavioral studies of full recovery times from SCI in tadpoles [11, 42].
Identifying functional processes in successful CNS axon regeneration
An initial examination of DESR genes supported their biological relevance to CNS axon regeneration, since many had been previously implicated in aspects of neural injury, including (but not exclusive to) axon regeneration [e.g., socs3 (Suppressor of Cytokine Signaling 3 [32, 107, 134, 151]), leptin [9, 13], fabp7 (fatty acid binding protein 7, a.k.a. brain-lipid binding protein [6, 42, 52, 111]), mapk8 (JNK1 [51, 97]), sox11 ; sdcbp (syntenin ), and prph (peripherin [41, 84], among others]. The hypothesis that CNS axon regeneration also shares functional processes with other forms of tissue regeneration was supported by some DESR genes having been previously implicated in successful regeneration in other contexts, including axolotl limb [e.g., tmsb4x [73, 109], anxa5 , and crabp2 )] and mammalian liver [e.g., top2A (DNA topoisomerase 2a , C9 (complement C9 [112, 133]), and rrm2 (ribonucleoside reductase regulatory subunit M2 )].
Thus, we felt justified in using DESR genes to further explore potential cellular processes in successful CNS axon regeneration. For this analysis, we combined two approaches. First, we conducted an automated gene ontology analysis (GO term; Metascape ) of the DESR genes at the peak period of regenerative axon outgrowth (1 week/ 11 days for SCI hindbrain and ONC eye, respectively), for significant enrichment (−log10(P) > 4) of processes associated with these genes. This analysis primarily yielded GO terms related to cell division and mitotic control, inflammatory responses, innate immunity, wound healing, and stress responses, as well as terms associated with RNA trafficking and turnover, JAK/STAT signaling, development and differentiation, myeloid cell activity, control of gene expression and chromatin remodeling, microfilament/microtubule dynamics, and cellular metabolism. Because relatively few DESR genes were represented at the early and late time points, GO term analysis for these time points was less fruitful, since many genes were omitted. Thus, we complemented the GO term analysis with a manual survey of every DESR gene’s known functions. For this stage of the analysis, we consulted PubMed and Gene Cards for the known functions related to each gene’s orthologs across species, paying special attention to functions related to injury (both traumatic and disease-related) and regeneration (both neural and non-neural), as well as processes identified by the GO term analysis. These functions are listed for each gene individually in Additional_File4_DESR_Data.xlsm. We then grouped DESR genes into categories. The results of this combined summary analysis, with genes arranged according to functional categories and time of differential expression in axon-regenerative CNS regions, are presented in Additional_File6_DESR_Functional_Categories.pdf.
In this analysis, DESR genes fell into eleven such categories (Fig. 4; ordered from largest to smallest): Inflammatory Response and Wound Healing, Cytoskeletal, Cell Signaling, Intracellular Transport, Post-transcriptional Regulation of gene expression, Epigenetic control of gene expression, Axon Outgrowth (tropic and trophic effectors and inhibitors), DNA Replication/Repair, Lipid Metabolism, Transcription Factors, and Cellular Metabolism. All eleven categories were represented during the peak regenerative phase, reflecting the wide range of cellular processes that are active then. Five of these eleven were limited to just the peak phase, indicating these processes were relatively inactive at other time points: 1) Cell Signaling, 2) Intracellular Transport, 3) Axon Outgrowth, 4) Lipid Metabolism, and 5) Cellular Metabolism.
DNA Replication/Repair DESR genes were found at only the early and middle time points. At the early time point, one gene in this category was down-regulated in successful but up-regulated in unsuccessful CNS axon regeneration, suggesting it is important for determining regenerative success (hp1bp3, a mitotic checkpoint protein involved in maintaining heterochromatin integrity during the G1/S transition [30, 31]). At the middle time point, 5 additional mitotic checkpoint proteins were upregulated, suggesting mitotic checkpoint-related genes are especially important for promoting the peak phase of regenerative outgrowth. Because injured neurons are post-mitotic, these checkpoint genes may have other functions in regeneration than regulating mitosis. In that regard, mitotic checkpoint genes are increasingly understood as important for giving cells time to rearrange their chromosomes in preparation for major changes in gene expression . The remaining twelve up-regulated genes at 7/11 days are typically considered to be more directly involved in DNA replication and repair. They may thus reflect the proliferation of glia and myeloid cells that influence the injury response, since the axotomized neurons themselves are post-mitotic. None of these DNA Replication/Repair genes were found at 3 weeks, indicating that by this time, such activities were no longer active.
The remaining five categories spanned all three time points, providing insights into how these functions evolved throughout recovery. The largest category comprised Inflammatory Response DESR genes. These genes represented functions related to regulating inflammation, innate immunity, wound healing, and cell survival responses, as well as tissue repair and regeneration. Such genes are known principally from mammalian studies, wherein inflammation is crucial for initiating and maintaining pro-regenerative responses to injury . However, because inflammation also damages tissues, it has been difficult from mammalian studies to parse inflammation-related pro-regenerative genes from those detrimental to recovery. Their preferential up- or down-regulation in Xenopus axon-regenerative CNS tissues thus provides clues as to which genes are beneficial to regeneration. Moreover, the involvement of many of these genes in wound responses and inflammation is known from studies that either included the injury site or involved non-neural tissues, such as liver or fish tail fin. Thus, their identification in this study also shows their relevance extends more broadly, beyond tissue repair, to CNS axon regeneration. At the early time point, up-regulated Inflammatory Response DESR genes (i.e., recovery promoting) were dominated by ones typically considered pro-inflammatory (7 genes) and genes activated by JAK/STAT signaling through cytokine receptors (3 genes). The two down-regulated genes in this category (i.e., recovery inhibiting) at this early time point were a cytokeratin that is otherwise up-regulated at mammalian wound sites (krt6a) and a pro-inflammatory, calcium-binding protein secreted by macrophages (ocm2). At the middle time point marking the peak axon-regenerative response, there were two up-regulated pro-inflammatory genes held over from 3 days, and seven new ones, along with many additional up-regulated genes implicated in other inflammatory and wound-healing processes: protein ubiquitination and turnover (11 genes), chaperones (4 genes), cell survival and tissue repair (10 genes), stress response (7 genes), myeloid cell activity (5 genes), and the transition from M1 to M2 macrophages (3 genes). Conversely, down-regulated DESR genes (i.e., regeneration inhibiting) at the middle time point included genes previously associated with exacerbating inflammation, cell death, and scar formation (10 genes), genes associated with oxidative stress (2 genes), myeloid cell activities (2 genes), and with maintaining the blood-brain barrier (1 genes), plus a heat-shock protein (1 gene). At the late time point (3 weeks), all nine Inflammatory Response DESR genes were up-regulated genes previously implicated in protecting cells from detrimental hyper-inflammation and in promoting cellular repair and regeneration, mostly in contexts not previously directly associated with CNS axon regeneration.
The second largest category of DESR genes found at all three time points were those associated with the cytoskeleton. At three days, all four cytoskeletal-related DESR genes were down-regulated: two cytokeratins and two specific actin subtypes. One of the cytokeratins, krt6a (previously mentioned in wound healing), has been previously associated with epithelial lesions in mice . Knocking it out makes such wounds more fragile , suggesting that in the CNS, this gene helps remodel tissues to promote early stages of axon outgrowth. Moreover, because our sampled tissues did not include the lesion, this function is likely important even at distances removed from the injury site itself. At the peak phase of regenerative axon outgrowth, all the cytoskeletal-related DESR genes were specific tubulin and actin subtypes, along with proteins associated with regulating microtubule- and microfilament-mediated transport and dynamics. All but four were up-regulated (27 total). Collectively, they emphasized the importance of microfilament and microtubule-associated dynamics for axon outgrowth, intracellular transport, cell motility, and cellular proliferation for successful axon-regeneration, while at the same time confirming and extending the identification of specific genes associated with these processes in CNS axon regeneration. For example, tuba1a and tubb2b, which are neuronal α- and β-tubulin subtypes, had been previously associated with developing and regenerating CNS axons [88, 144], whereas dynnl2, a minus end-directed microtubule motor protein involved in the retrograde transport of proteasomes, is novel to axon regeneration. At the late time point, there were only two cytoskeletal-related DESR genes, both of which were up-regulated. They were an intermediate filament gene, prph (peripherin), which was already known to be up-regulated in both reactive glia and regenerating optic axons in Xenopus [41, 84]. Prph was also up-regulated in both regenerative and non-regenerative cases at earlier time points, but its preferential up-regulation at 3 weeks in axon-regenerative CNS over non-regenerative CNS suggests its importance in successful regeneration persists into late recovery. The remaining up-regulated cytoskeletal-related gene at 3 weeks, ebf3, has dual functions in promoting microtubule bundling and as a transcription co-factor that inhibits gliogenesis [17, 72], suggesting novel functions for promoting recovery during the late phase.
The third category of DESR genes represented at all three time points comprised genes involved in post-transcriptional control of gene expression. Post-transcriptional control of gene expression is already known to be important for successful CNS axon regeneration [51, 76], but the specific genes involved are still being discovered. This study suggests new genes. For example, at the early time point, DESR genes in this category primarily reflected changes in translational control of mRNAs. This is an understudied, potentially important function for promoting successful recovery from CNS injury, because stressed cells utilize cap-independent mRNA translation to ensure that proteins needed for survival are synthesized, while they simultaneously inhibit cap-dependent translation to facilitate cellular reprogramming in preparation for the next phase of the stress response . Our data identified eif5b, a translation initiation factor that promotes IRES-dependent mRNA translation , as one such gene potentially important for CNS axon regeneration. Another Post-transcriptional Regulation DESR gene at 3 days was a 60S ribosomal protein, rplp1, which is essential for brain development due to its effects on cyclin and p63 protein expression . During neural development, specific ribosomal proteins are often needed for selective translation of individual mRNAs , and this gene was down-regulated in the two regenerative tissues but up-regulated in non-regenerative frog SCI hindbrain, suggesting a purposeful role in determining regenerative success. The middle time point saw the up-regulation, and no down-regulation, of multiple Post-transcriptional Regulation DESR genes, spanning a range of RNA-related functions, including eight splicing factors, four additional ribosomal subunits, five mRNA translation initiation and elongation factors, and three regulators of mRNA trafficking and turnover. These most likely represent post-translational aspects of the onset of expression of the many genes involved in CNS axon regeneration. Most of these DESR genes were newly implicated here in CNS axon regeneration, although a few had been implicated previously in related contexts. For example, two up-regulated splicing factors, snrpd3 and snrpn, have been associated with neurodegeneration in Spinal Muscular Atrophy  and with developmental axon outgrowth , respectively. Prmt1 is the primary methylase targeting hnRNP K , which is an RNA-binding protein that regulates nuclear export and translation of cytoskeletal-related mRNAs essential for optic axon regeneration in Xenopus . Its selective up-regulation in the two regenerative CNS regions suggests a previously unrecognized role for methylation of hnRNP K, as well as its histone targets, in recovery from CNS trauma. The late recovery phase saw only the up-regulation of two transcripts of RNA-binding proteins that regulate trafficking, translation, and turnover of specific mRNAs: AldoA stabilizes mRNAs of neurofilaments ), which consolidate axonal growth and expand axon caliber once growing axons contact their targets [118, 120, 146]. Thus, this aldoA increase is consistent with 3 weeks being a phase of consolidation and refinement of synaptic connections . The other up-regulated mRNA, mex3a, encodes an RNA-binding protein implicated in regulating neurogenesis and degeneration , suggestive of a previously unrecognized role for mex3a target-transcripts in CNS axon regeneration.
The fourth most numerous category with representation at all 3 time points (26 genes) comprised genes implicated in epigenetic control of gene expression, which is increasingly recognized as crucial for regeneration of organs and tissues in other contexts (reviewed in ). At the 3-day time point, the two DESR genes in this category were down-regulated, specialized histone variants: hist1h4k, an H4 gene variant implicated in protecting cells from DNA damage by facilitating DNA double strand break repair , suggesting it may play a role here in protecting injured cells; and hist2h2ab, an H2 gene variant involved in nucleosome repositioning in preparation for transcription ). Nucleosome repositioning is a theme that persisted into the peak phase of regenerative axon outgrowth, with hist2hab at this next phase now switching from a down-regulated to an up-regulated DESR gene, and five additional genes known to be involved in nucleosome repositioning (four up-regulated; one down-regulated) appearing. The remaining 19 DESR genes in the Epigenetic category (13 up-regulated, 6 down-regulated) are all known to play various roles in enzymatically modifying histones and DNA through acetylation and methylation. Two of these were down-regulated components of the Polycomb Repressive Complex (PRC), newly implicating it here in successful CNS axon regeneration, although the PRC has already been implicated in tissue regeneration in other contexts [23, 27, 81, 122]. One of these two PRC genes, jarid2, persisted as the sole epigenetic-related gene in the late recovery phase, when it continued to be down-regulated.
The final DESR category represented at all three time points comprised genes associated with transcriptional control. At the earliest time point, only one such gene was up-regulated, ddit3, which is a C/EBP-related transcription factor needed to activate pro-inflammatory signals . Ddit3 has been previously reported as up-regulated in response to optic nerve injury in mouse  but had not been previously implicated in successful CNS axon regeneration. No transcriptional control DESR genes were down-regulated at this early time point. The peak phase of regenerative axon outgrowth saw eight up-regulated and three down-regulated transcription factors and co-factors. Two of the up-regulated ones were already known from previous studies in mammals to play important roles in SCI – hes5 and sox11 [56, 95, 149]; two additional ones (bcl6 and mllt11) were known from cancer studies as modulators of STAT3 [89, 99], which is a transcription factor known to play an important role in regulating inflammation-mediated responses during tissue regeneration. However, to our knowledge, neither of these two regulatory STAT3 co-factors had previously been implicated in successful CNS axon regeneration. Another down-regulated gene in this category at the middle time point was znf395, a zinc-finger transcription factor that activates pro-inflammatory cytokines , newly implicating it here in limiting aspects of inflammation that are detrimental to CNS axon regeneration. The two sole transcription factor DESR genes at 3 weeks were both up-regulated (ebf3 and irf8). Each has known roles in regulating genes involved in apoptosis and in suppressing a hyper-immune response in macrophages, respectively [53, 72]. These functions are consistent with limiting inflammation and its detrimental effects continuing to be important well into the late stages of recovery.
The KEGG Adipocytokine signaling pathway comprised multiple DESR genes
Gene regulatory networks often comprise genes with different cellular functions. Therefore, to identify gene regulatory networks potentially involved in successful CNS axon regeneration, we analyzed KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways for enrichment in DESR genes. In one such analysis we input all the DESR genes found at more than one time point, with the aim of highlighting regulatory networks that were active throughout regeneration. This analysis identified a single network, the KEGG Adipocytokine signaling pathway (map04920; Fig. 5). This network contained four up-regulated DESR genes: lep (leptin), socs3 (suppressor of cytokine signaling 3), mapk8 (JNK1) and acsbg2 (Acetyl-CoA synthetase bubblegum family member 2; labeled as FACS in Fig. 5). The up-regulation of these genes suggested they were all pro-regenerative in this context, which was surprising for socs3, because in mammalian studies, it is generally considered inhibitory to CNS axon regeneration  (but see  for a potential resolution to this paradox). In mammalian studies, this pathway interconnects the control of cellular metabolism with inflammatory and stress responses, integrating downstream pathways known from studies across different systems to be involved in various aspects of tissue regeneration. These sub-pathways include JAK/STAT signaling [15, 16, 107, 134], control of axonal transport and synthesis of axonal cytoskeletal proteins via JNK signaling [51, 97], cap-dependent mRNA translation through mTOR [1, 82, 100], and the metabolic regulation of both DNA methylation and histone acetylation/methylation [148, 150]. In addition to these four DESR genes, a further 60% (69/115) of the annotated X. laevis genes belonging to this pathway were found to be potentially differentially expressed (either up- or down-regulated) upon injury (p ≤ 0.02) during at least one time point, in at least one of the three tissues (Fig. 5). Such complex, temporally variable behavior of these additional network elements suggested that maintaining flexibility of this network throughout regeneration is an important feature for CNS axon regeneration. Collectively, these properties implicate this network as a key hub within the larger network of CNS injury response genes.
The majority (55%) of DESR genes exhibiting opposing expression between successful and unsuccessful regeneration were functionally inter-connected
We next reasoned that genes that were significantly, differentially expressed in opposing directions in the two regenerative cases (R) vs. the non-regenerative case (NR) (i.e., up-regulated in the first two and down-regulated in the last, and vice versa) would be especially instructive in providing insights into changes in gene expression that are purposeful for determining the success or failure CNS axon regeneration. Thirty-three DESR genes exhibited such opposing expression (Fig. 6a; also green and red fonts in Additional_File6_DESR_Functional_Categories.pdf). Genes up-regulated in the regenerative cases and down-regulated in the non-regenerative case, and vice versa, were classified as regeneration-promoting and -inhibiting, respectively. The identities and known functions of thirteen of these genes supported this classification. Nine regeneration-promoting genes were previously known to promote either tissue regeneration and axon outgrowth in general, or cell survival in trauma and degenerative diseases: ltf (lactoferrin, an iron-binding, neuroprotective gene ), otop3 (otopetrin 3, a member of a gene family that down-regulates inflammation ), mcm6 (minichromosome maintenance complex component 6, implicated in Müller cell-derived neurogenesis in retinal regeneration ), tubb2b and tuba3d (neuronal tubulins upregulated in axonogenesis [86, 88]), ttl (tubulin tyrosine ligase, required for retrograde transport of pro-regenerative signals after axotomy ), fabp3 (fatty acid binding protein 3, linked to partially successful recovery from spinal cord injury in opossum neonates ), fabp7 (brain lipid binding protein, a marker for radial glial endfeet, which guide growing axons ), and fads1 (fatty acid desaturase 1, linked with suppressing inflammation in liver disease ). The functions of the previously known four inhibitory genes were all associated with aggravating inflammation, promoting degeneration, and inhibiting tissue regeneration in various contexts: enpp2 (autotaxin, a stimulator of inflammation [22, 63]), slc9a3R2 (a suppressor of STAT3, pro-healing signaling in colon cancer ), plp1 (proteolipid protein 1, a major component of mature oligodendrocytes, which inhibit CNS axon regeneration [65, 155]), and znf395 (zinc finger protein 395, a transcriptional activator of pro-inflammatory cytokines ). The remaining twenty oppositely expressed DESR genes had not previously been associated with processes potentially relevant for CNS injury, inviting further study.
To further explore potential functional inter-relationships among these 33 oppositely expressed DESR genes, we performed STRING analysis for known/predicted protein-protein interactions and functional inter-relationships  (Fig. 6b). Eighteen of the 33 genes had such interactions, generating a network of five interconnected clusters (confidence = 0.300): 1) two regulators of mRNA translation; 2) three specialized histones; 3) six other genes associated with epigenetic changes in gene expression, including one such regulated gene that is part of and therefore potentially connects this network to the adipocytokine signaling pathway, abcb1 ; 4) three metabolic enzymes, one of which (idh1) regulates the TET enzymes that modulate DNA methylation states ; 5) and four microtubule-related genes, which serve to drive cell motility, axonal transport, and axon outgrowth (e.g., ). These genes thus form a nascent gene regulatory network potentially linking changes in cellular metabolism with epigenetic reprogramming and post-transcriptionally regulated changes in gene expression, converging upon cytoskeletal structural genes to drive axon outgrowth. At the regulatory heart of this network were genes involved in the epigenetic control of gene expression (dotted-lined box, Fig. 6).
In situ hybridization of retina to identify cell types expressing key DESR genes
To begin to identify cell types responding to axotomy through expression changes in key DESR genes, we performed in situ hybridization on retina at the peak phase of regenerative axon outgrowth for select genes, focusing primarily, but not exclusively, on the two networks just discussed (Fig. 7). We chose retina over hindbrain for these studies because the anatomical organization of retina greatly simplifies identifying cell types. For example, because frog retina receives no external afferents, ONC damages only the optic axons, which originate from retinal ganglion cells. These cells constitute the vast majority of cells in the ganglion cell layer, making them anatomically distinct from other local circuitry neurons, such as bipolar cells, and Müller cell astroglia. Both cell types can be indirectly affected by optic nerve injury and their cell bodies occupy other retinal layers [2, 84]. The situation is more complex in hindbrain, which harbors not only the axotomized neuronal cell bodies but also ascending axons damaged by the SCI. In addition, the axotomized hindbrain neurons are anatomically intermixed with other local neurons, making it essential to distinguish them through applying retrograde tracers . Finally, because the injured and control hindbrain are necessarily in different animals, it is harder to ensure that hybridization conditions are precisely matched than it is for retina, where both the operated and unoperated eye occupy the same section.
For key up-regulated DESR genes, we chose leptin and socs3 from the Adipocytokine Signaling pathway, prmt1, which methylates both hnRNP K and histones, and idh1 and ezh2, which are two epigenetic-related genes from the network of Fig. 6b. For down-regulated genes, we chose jarid2 and suz12, which belong to the Polycomb Repressive Complex 2 and are also part of the same network as idh1 and ezh2 (Fig. 6b). These genes were also chosen because they represented a wide range of expression levels, from abundantly expressed (FPKM > 100) to moderately and poorly expressed genes (FPKM < 50), and because they spanned the range of FDRs from <.002 to 0.05. Relative differences in expression from experimental to control (|log2(fold-change)|) also spanned the full spectrum – from 0.3 (the smallest significant change) to 2.8 (one of the largest). (It should be noted here that because probes for in situ hybridization do not distinguish between homeologs, the RNA-seq data in Fig. 7 represent values combined for the two homeologs). In all cases, the overall intensity of labeling and the relative difference between operated and unoperated eyes reflected the magnitudes seen in the RNA-seq data. Increases in the two up-regulated adipocytokine signaling pathway genes, and in idh1 and ezh2 were confined to the retinal ganglion cell layer. The changes in expression of the remaining up-regulated gene, prmt1, as well as the two down-regulated components of the Polycomb Repressive Complex 2 (jarid2, suz12), included but also extended beyond the retinal ganglion cell layer, indicating that these changes in expression occurred more globally throughout the retina.
We used a novel three-way comparison between a region of the CNS that successfully regenerates axons throughout life (retina) against another that successfully regenerates axons only transiently in development (hindbrain), before and after the transition, to serve as a biological data filter for parsing core features of a successful CNS axon-regenerative response from the thousands of injury-induced changes in gene expression that otherwise represented region-specific responses. The first indication that the two regenerative cases uniquely shared core features was the temporal pattern of the number of injury-induced genes. Whereas this number reached a maximum in the two regenerative cases during the peak phase of axonal regrowth, it did so in the non-regenerative case at the earliest time point. PCA corroborated the conclusion that gene expression profiles were largely tissue-specific, but closer inspection identified shared features between the two regenerative tissues. For example, PCA eigenvectors (Fig. 3) and scatter plots (Additional_File5_PCA_Scatterplot.pdf) differed more between injury and controls for both regenerative cases than was the case for the non-regenerative case. The lack of overlap between regenerative tadpole and non-regenerative frog SCI hindbrain, even for different time points, further indicated that the genes induced by SCI were quite different between these two tissues throughout regeneration. Individual DESR genes (324) far out numbered injury-induced genes shared by all three tissues (19 genes) at all time points, as well as the number shared between regenerative and non-regenerative SCI hindbrain at the peak phase of axonal regrowth (284 DESR genes vs. 101 shared regenerative/non-regenerative SCI hindbrain genes). Collectively, these data supported the conclusion that beginning relatively soon after injury, having the capacity to mount a successful axon-regenerative response stimulates a markedly different injury response from that of a non-regenerative tissue, which in the latter, is more focused on attempting rapid repair for survival than on axon regeneration.
Further insights into the nature of the cellular processes uniquely shared between the two regenerative cases came from the identities and subsequent analyses of the DESR genes. In support of these genes being relevant to CNS axon regeneration, DESR genes included multiple genes previously implicated in CNS injury, although mostly in other contexts, such as studies that included the injury site, peripheral axon regeneration, and resistance to neurodegenerative disease. This observation emphasized the universality of those genes in other types of CNS injury, as well as the phylogenetic conservation of the processes they represented, thereby building confidence in the data set. The data also included a large number of DESR genes shared across a broad spectrum of non-neural regenerating tissues, including zebrafish tail fin, axolotl limbs, and mammalian liver, arguing further that their importance extends beyond CNS axon regeneration to vertebrate tissue regeneration in general, and providing clues into aspects of a successful regenerative response that are universally shared among all vertebrate regenerating tissues. The many additional, novel genes implicated here for the first time in CNS axon regeneration provide new targets for investigation.
Surveying the literature for what is known about each DESR gene’s function sorted them into eleven broad functional categories. Temporal differences in the composition of these categories provided further insights into how core processes associated with successful CNS axon regeneration shifted as axon regeneration progressed. For example, the onset of axon regeneration (3 days) was dominated by up-regulated pro-inflammatory genes, several down-regulated cytoskeletal genes potentially associated with scarring, and by genes indicative of cellular reprogramming events. The latter included an up-regulated gene associated with promoting non-cap dependent mRNA translation (eif5b) and several down-regulated ribosomal subunits, two down-regulated histone subtypes associated with transcriptional activation, a down-regulated mitotic check point protein (hp1bp3), and an up-regulated C/EBP-related transcription factor (ddit3) that is regulated by stress factors associated with bacterial infections and with mammalian optic nerve injury but which had not previously been implicated in successful CNS axon regeneration [34, 101]. At the peak phase of regenerative axon outgrowth, all eleven functional categories were represented, indicating high levels of activity associated with these functional categories at the most active time point. While the presence of DESR genes previously implicated in regenerative and developmental axon outgrowth helped to provide confidence in the relevancy of the data to CNS axon regeneration, our finding additional genes associated more generally with cellular processes previously associated with wound healing and tissue repair in other contexts provides new gene targets for studying how the cellular processes with which they are connected function specifically in CNS axon regeneration. These processes included modulators of Wnt-signaling , macrophage invasion and the complement cascade [35, 74, 104], and modulators of JAK/STAT signaling [54, 134]). Still more genes offered important clues about which aspects of an inflammatory response are pro-regenerative (i.e., up-regulated DESR genes) as opposed to detrimental (i.e., down-regulated DESR genes) to CNS axon regeneration, a distinction that has been difficult to make from mammalian studies . Other DESR genes implicated cellular processes that are new to CNS recovery (e.g., regulators of mitotic checkpoints). The late time point, when axonal regrowth was largely finished but synaptic connections were still being refined, was dominated by up-regulated genes known to promote cellular repair, suppress hyper-inflammation, and stabilize the cytoskeleton. These functions are understandable during a time when the initial effects of trauma have been overcome, terminal axonal arbors are being refined and consolidated, and physiological gene functions are returning; the identification of new, specific genes involved in these processes provide new entry points for studies aimed at understanding and promoting these beneficial processes.
To begin to elucidate gene regulatory networks potentially occupying important hubs in the network of gene responses promoting successful CNS axon regeneration, we engaged in two exercises. In the first, we surveyed KEGG pathways enriched with DESR genes, especially genes that were differentially expressed at multiple time points. This analysis identified the KEGG Adipocytokine signaling pathway as enriched both for DESR genes and for injury-induced differentially expressed genes in general, independent of tissue regenerative capacity. The four DESR genes in this pathway were uniformly up-regulated, suggesting that their activation within this network promotes regeneration. However, multiple other differentially expressed genes belonging to this network exhibited both up- and down-regulation at various times in each of the tissues, further suggesting that the dynamic behavior of this network influences regenerative success. The importance of individual components of this pathway for wound repair and regeneration had already been noted across a range of systems, both neural and non-neural. For example, leptin in mammals has been implicated in the regenerative success of both liver and muscle [67, 91, 145], and in zebrafish, it has been implicated in both restoration of a damaged retina  and in recovery from heat-induced stress . In Xenopus, up-regulation of leptin has also been noted in the regenerative response of the spinal cord itself to SCI , and during development, it is involved in activating canonical Wnt-signaling in neurogenic regions of the brain . Our linking leptin to axon-regenerative regions of the CNS that are otherwise uninjured, and localizing the increased expression to the axotomized retinal ganglion cells, which do not themselves undergo active proliferation after axotomy, suggests that these cells use leptin to signal nearby support cells, such as invading macrophages [45, 152] and reactive Müller glia . Canonically, leptin and the adipocytokine signaling pathway play a crucial role in thermoregulation and dietary intake in mammals [18, 26]. Thus, our data from Xenopus adds to a growing body of evidence supporting a phylogenetically ancient linkage in vertebrates between the regulation of cellular metabolism, the stress response, and the success or failure of organ and tissue regeneration both within and outside the CNS, where the axotomized neurons themselves may be actively involved in the signaling.
Next, we examined the subset of DESR genes that exhibited opposing expression patterns between regenerative and non-regenerative tissues, reasoning that these genes could provide insights into purposeful differences between the two. For example, genes that were up-regulated in the regenerative cases but down-regulated in the non-regenerative case would likely be pro-regenerative and vice versa. Of the 33 oppositely expressed DESR genes, nine pro-regenerative (i.e., upregulated in axon-regenerative CNS and down-regulated in non-regenerative CNS) and four regeneration inhibitory (i.e., down-regulated in axon-regenerative CNS and up-regulated in non-regenerative CNS) genes had previously been implicated in either promoting or inhibiting regeneration and cell survival, respectively, in one or more contexts, thereby helping to validate the list while expanding its relevance to include CNS axon regeneration. Whereas KEGG analysis fits gene sets into known biological pathways, STRING analysis can be used to discover new gene networks from known interactions, both physical and functional. STRING analysis of these 33 oppositely expressed DESR genes placed the majority (18) into five interconnected clusters, with three clusters comprising specialized histones and genes associated with epigenetic control of gene expression at the center. These five clusters could thus conceivably comprise the core of a nascent gene regulatory network for successful CNS axon regeneration. In situ hybridization of a subset of genes in this network found that all the pro-regenerative genes were specifically up-regulated in the axotomized retinal ganglion cells, whereas the two regeneration-inhibitory components of the network, which are part of the Polycomb Repressive Complex 2, were down-regulated in multiple retinal layers. Surprisingly, another member of Polycomb Repressive Complex 2, ezh2 was up-regulated in axotomized retinal ganglion cells. Since ezh2 functions extend beyond the complex to include methylation of DNA as well as histone sites , its pro-regenerative activity is therefore likely to be separate from its functions within the Polycomb Repressive Complex. Examination of the known functions of the genes in the central three clusters, together with other epigenetic-related DESR genes (25 total), strongly implicate increasing the accessibility of DNA for transcription as being of central importance to determining the success or failure of CNS axon regeneration.
A novel three-way comparison identified candidate genes belonging to core processes shared by two CNS regions that successfully regenerate axons, that distinguished them from one that does not (i.e., juvenile frog ONC eye vs. tadpole SCI hindbrain vs. juvenile frog SCI hindbrain, respectively). Unlike earlier studies of SCI that included the lesion in the analyzed tissues, by analyzing the tissues of origin of the regenerating axons, this study emphasized genes relevant for CNS axon regeneration over wound repair and tissue restoration. Even so, many of the genes identified in our study were previously implicated in wound healing, tissue repair, and cell survival, demonstrating that these processes are shared with successful CNS axon regeneration. Because these studies were mostly conducted in other tissues and species, our study both demonstrated their phylogenetic conservation and identified many new individual genes associated with these processes as worthy targets of study in CNS axon regeneration. Finally, because many of these genes were implicated in mammals, where recovery is often incomplete, their identification in this study provides clues as to which are beneficial to recovery and which ones are detrimental. Network analyses of DESR genes that were up-regulated at multiple time points and of DESR genes that were oppositely expressed in axon-regenerative vs. non-regenerative CNS regions, respectively, identified the Adipocytokine signaling pathway and a novel gene regulatory network with epigenetic control of gene expression at its core, as important hubs in the larger network of injury-response genes associated with successful CNS axon regeneration. In situ hybridization in ONC retina of representative genes from these networks showed that the injury response of some elements of these networks were restricted to the axotomized neurons, whereas others extended to neighboring cells, demonstrating that a successful CNS axon-regenerative response to axotomy indirectly invokes responses in neighboring cells that could influence its success. This study thus provides both a resource and a starting point for future investigations into the molecular underpinnings of successful CNS axon regeneration.
Animal and surgical procedures
All animal procedures were approved by the Institutional Animal Care and Use Committees (IACUC) of the University at Albany (optic nerve crush) and Morehead State University (spinal cord transection), in accordance with the National Institutes of Health Guide for the Care and Use of Laboratory Animals. For the sake of consistency, Xenopus laevis tadpoles and juvenile frogs were all from an albino strain obtained from the same supplier (Xenopus Express, Brooksville, FL). Prior to surgery, animals were acclimated to the lab for at least a week, under a 12-h light:12-h dark photoperiod, in tanks of conditioned water at 22 °C, where they were fed every other day (Xenopus Express SFF and TP, for frogs and tadpoles, respectively). Optic nerve crush (ONC) surgeries were performed on fully anesthetized (0.1% MS222, Sigma-Aldrich), 1–3 month, post-metamorphic albino frogs (unsexed juveniles), as described . Briefly, the surgically exposed optic nerve on the right side of the animal was crushed at the orbit, approximately 0.5–1 mm from the eye, using #5 Dumont forceps (Fig. 1d), 2–4 times, until the nerve turned visibly clear within the epineural sheath. The contralateral left eye served as the unoperated control. Spinal cord transections (SCI) were performed on fully anesthetized NF  stage 53 tadpoles (0.02–0.04% MS222) and 1–3 month post-metamorphic juvenile frogs (Fig. 1d) of the same age as those used for ONC [42, 44]. The surgically exposed spinal cord was completely transected at the mid-thoracic level, and the success of the surgery was confirmed at the time of surgery by passing a Minutien pin (Fine Science Tools) through the injury site while observing the spinal cord. To further confirm that spinal cord transections were complete, animals were tested the next day for residual movement. After surgery, animals were maintained in tanks filled with aerated, sterile-filtered, conditioned rearing water containing 4 mg/l gentamicin sulfate (Sigma-Aldrich) to guard against infection.
For each biological replicate, either five hindbrains (SCI) or six eyes (ONC) were rapidly dissected from fully anesthetized (0.1% MS22) animals [3, 41, 42]. To maximize RNA integrity, hindbrains were submerged in RNAlater (ThermoFisher Scientific) at room temperature, and eyes were homogenized immediately. Tissues were homogenized (Polytron® PT 3000) at 23 °C in 500 μl of Lysis Buffer Q (Norgen Biotek Corp., Ontario, Canada), supplemented with 5 μl β-mercaptoethanol. Homogenates were then cleared by centrifugation (23 °C, 15,000 rpm, 15 min), and total RNA was isolated [RNA/DNA/Protein Purification Plus Kit, Norgen Biotek Corp (cat# 47700); RNase-free DNase I kit, Norgen Biotek (cat# 25710)]. A 2 μl aliquot of this RNA solution was quantified (Nanodrop® ND-1000 Spectrophotometer), and RNA quality was preliminarily assessed by agarose gel electrophoresis (E-Gel® EX Agarose, E-Gel® iBase™ Power System, ThermoFisher Scientific).
Library preparation and sequencing
Purified total RNA from all 51 samples was simultaneously sent for library preparation and sequencing (Center for Genomics and Molecular Medicine, University of Louisville). There, RNA was again quantified (Qubit™ Fluorometer) and assessed for quality [2100 Bioanalyzer, Agilent Technologies; RINʼs (RNA Integrity Number): 9.0 ± 0.7 (mean ± SD; n = 51), range 7.6–10.0]. Ribosomal RNA-depleted RNA (0.5 μg) was subsequently used to make barcoded cDNA libraries [TruSeq Stranded Total RNA LT Sample Prep Kit - Set A and B (Illumina, cat# RS-122-2302), with Ribo-Zero® Gold (Illumina)], which were pooled for sequencing [NextSeq 500; NextSeq 500/550 75 cycle High Output Kit v. 2 (Illumina)] in two runs to nominally yield 30 × 106, 75 base-pair (bp), single-end reads.
Read alignment and differential expression analysis
After assessing the raw fastq files for quality (FastQC v0.11.5 ), reads were aligned against the Xenopus laevis genome (Xenbase v9.1; http://www.xenbase.org/ RRID:SCR_003280) using Bowtie2 (v2.2.9)  and TopHat (v2.1.1) , and then annotated using the Mayball gene model [24, 106, 108]. The alignments yielded 34.4 ± 3.1 (S.D.) million successfully aligned reads per sample (Fig. 1f1), with only 9.6 ± 2.9%, (S.D.) reads initially flagged as potentially duplicate alignments, which can occur in Xenopus laevis due to its ancestrally (~ 30 Mya) duplicated genome . These potential duplicate alignments were resolved for the vast majority by using the alignment with the higher score to assign them to separate genes on different chromosomes (referred to in X. laevis as S and L homeologs). For the small remainder (< 10% of potentially duplicate alignments), reads were randomly distributed between the two homeologs. To confirm the accuracy of the alignments, a subset was visualized directly (Integrative Genomics Viewer (IGV) (v2.3.88) [115, 136]).
Differential expression and associated analyses [dispersion plots, gene density histograms, and eigenvector principal component analysis(PCA)] were carried out using Cufflinks/Cuffdiff2 (v.2.2.1; per-condition dispersion with a minimum count of 10) [113, 114, 137, 139, 140] and associated utilities of CummeRbund (v2.16.0) . The gene density histograms were used to set the threshold (FPKM = 0.45) for active gene expression (Fig. 1f2), and the dispersion plots were used to affirm that samples being compared statistically had similar dispersions (Fig. 1f3). In CuffDiff2, reads are normalized for transcript length, which enables comparing expression levels among different genes, and a difference of means test is performed across the biological replicates on the log2(fold-changes) [i.e., log2(FPKMinjury/FPKMcontrol)] for each gene separately to obtain a p value. These p values are then ranked to calculate an FDR (q) . Differences in individual tissues were considered significant for FDR < 0.05, without regard to the magnitude of the change.
The same read alignments and gene counts, after they were rlog-transformed, were also analyzed by DESeq2 [4, 79], and scatter-plot PCA representations were generated using the plotPCA function of DESeq2. For the latter, multiple PCA plots were generated over the range of 500 to 20,000 of the most highly-expressed genes, for which no qualitative differences were observed. DESeq2, unlike CuffDiff2, does not normalize genes for transcript length and uses a variance test, which incorporates the variance of all the genes in a sample, instead of a means test that treats the variance of each gene separately. Because the tissues being analyzed comprised a complex mix of cell types, each of which may respond differently to the lesion, we used the CuffDiff2 data sets for subsequent downstream analyses.
Gene ontologies of differentially expressed genes
For each time point, overlapping differentially expressed genes across the three conditions were identified using Awk scripts. To identify DESR genes, the script first identified differentially expressed genes in both regenerative tissues, then removed any that were also present in the non-regenerative one. The final DESR gene list was then analyzed for gene ontology (GO term) enrichment (Metascape (v3.0) ). Then, a literature search (PubMed; GeneCards) was conducted for each DESR gene for additional information related to its function that was relevant to processes potentially involved in CNS regeneration [i.e., wound healing, inflammation, innate immunity, tissue regeneration or degeneration, and axon outgrowth], and as suggested by the GO term analysis, to yield the final classification of genes into 11 functional categories.
KEGG pathway and STRING analysis
DESR genes appearing at multiple time points were entered into the KEGG pathway analysis tool (Kyoto Encyclopedia of Genes and Genomes [57,58,59]; https://www.genome.jp/kegg/). Oppositely regulated DESR genes (i.e., up-regulated in both regenerative tissues but down-regulated in the non-regenerative tissue, and vice versa), were analyzed by the Search Tool for Retrieval of Interacting Genes/Proteins to identify potential functional interactions (STRING, v10.5 [129, 131]; https://string-db.org; confidence level, 0.300; K-means clustering = 5). The Human Protein Complex Map (hu.MAP; http://proteincomplexes.org/) was used to identify proteins that physically interact within macromolecular complexes .
In situ hybridization
cDNA templates for in vitro transcription of antisense cRNA probes for in situ hybridization were prepared by reverse transcriptase (Super Script IV VILO, Invitrogen) polymerase chain reaction (Phusion High Fidelity PCR, New England Biolabs) from total RNA isolated from juvenile frog brain, with the following primers: Leptin (LEP.L; NCBI Accession number XM_018252815; length 278 nt) forward – GAT CCA AGG ACG AGC TAT AAA AAC T, reverse – GTA ACA GAC TGC GGA GGT TCT; SOCS3.S (NCBI Accession number NM_001087305; length 611) forward– ATG GTA ACG CAG AGC AAG TTC CCG, reverse– CGT TTT CTT TGT CTA CAC TGG GGA; EZH2.L (NCBI Accession number XM_018266394; length 251 nt) forward– GTA CAT GCG CTT ACG GCA AC, reverse– AGG CTA CAG CAG TGA GTG TT; IDH1.L (NCBI Accession number NM-001094553; length 250 nt) forward– AAC GCC AGG ATG TCC AAG AA, reverse– TCT CAT CAG GTG TAA TAG TGG CA; PRMT1.L (NCBI Accession number NM_001089302.1; length 275 nt) forward– GAG GCG AAG ACC TGC AAC AT, reverse– ACT TTC TTG GCA CCA GCC TT; SUZ12.S (NCBI Accession number NM_001130944; length 273 nt) forward– CTG TCA AAC CTG CAC AGA CAA, reverse– TGT CCT CTT TGG TCA CAT AGT TG; JARID2.L (NCBI Accession number BC086634.1; length 258 nt) forward – TGT GTT TTG CTT GGA GTG TGC, reverse – CAG GAT GAA GCA CTT TTG GAC A. Except for socs3, T7 RNA polymerase promoter sequences were added to the 5′ end of reverse primers for in vitro transcription directly from PCR products. PCR product sizes were confirmed on 1% agarose gels, and the products further purified (Monarch PCR and DNA Cleanup Kit, New England Biolabs). The socs3 RT-PCR product was cloned into a plasmid having a T7 promoter  for in vitro transcription. Digoxigenin-labeled antisense cRNA probes were generated by in vitro transcription (DIG RNA Labeling Kit (SP6/T7); Roche-Millipore Sigma).
For in situ hybridization, (3) fully anesthetized frogs (0.1% MS-222) were intracardially perfused first with Ringer’s solution containing MS-222 (112 mM NaCl, 2 mM KCl, 1 mM CaCl2, 1.2 mM NaHCO3, 0.1% MS-222, pH 7.4) and then with 4% paraformaldehyde in 0.1 M sodium phosphate (pH 7.4; PB) for fixation. Heads were then removed, cryoprotected in increasing concentrations of sucrose solutions (10, 20, 30%) in PB, then embedded in Tissue Freezing Medium (Triangle Biomedical Sciences) and cryo-sectioned (Leica CM1950) to obtain 18 μm thick transverse sections, which contained both eyes. These were thaw-mounted onto (+)-charged glass slides (Tissue Tak, Polysciences), which were subsequently stored at − 20 °C. For staining, slides were post-fixed in 4% paraformaldehyde in PB (20 min; 4 °C), then treated with 10 μg/ml proteinase K (Sigma-Aldrich; 6 min; 37 °C) and re-fixed in 4% paraformaldehyde (20 min, 4 °C). Slides were then incubated in 0.2 M HCl (30 min, room temperature) to inactivate endogenous alkaline phosphatase. After several washes in PB, slides were treated with 0.1 M triethanolamine (pH 8)/0.25% acetic anhydride, to further block non-specific hybridization, and subsequently washed in 2x saline-sodium citrate (SSC) buffer. Following a pre-hybridization incubation [50% formamide, 1x Denhardts, 20 mM tris (pH 7.4), 1 mM EDTA, 0.3 M NaCl, 0.1 g/ml dextran sulfate, 0.1 mg/ml herring sperm DNA, and 0.25 mg/ml yeast RNA], slides were hybridized with digoxigenin-labeled cRNA probes (100 ng in 200 μl/slide, 20 h, 50 °C). After hybridization, slides were treated with RNase A (20 μg/ml RNase A, 0.5 M NaCl, 0.25 mM EDTA, and 10 mM tris, pH 8) to further remove unhybridized probe. Slides were next washed in 2x SSC (room temperature), followed by a high stringency wash in 0.1x SSC (65 °C). Digoxigenin signals were detected by immunostaining with alkaline phosphatase-conjugated Fab fragments (Roche-Millipore Sigma), then developed with 4-nitro blue tetrazolium chloride/5-bromo-4-chloro-3-indoyl phosphate (Roche Millipore Sigma). The stained sections were then mounted under #1 coverglass in Fluoromount G (ThermoFisher Scientific). They were subsequently imaged using a Leitz Laborlux S (25X PL Fluotar, 0.6NA) microscope and a Nikon DS-Ri1 camera, using NIS-Elements BR 4.5 software (Nikon).
Availability of data and materials
The datasets generated and/or analyzed during the current study are available either in the Gene Expression Omnibus (GEO) repository (FASTQ files) [https://www.ncbi.nlm.nih.gov/gds][GSE 137844] or are included in this published article and its supplementary information files. Genomic resources for Xenopus laevis are available from Xenbase (http://www.xenbase.org, RRID:SCR_003280) and Xenopus Genome Project (http://www.marcottelab.org/index.php/Xenopus_Genome_Project).
Central nervous system
Differentially expressed in successful regeneration
False discovery rate
Fragments per Kilobase Million
Kyoto Encyclopedia of Genes and Genomes
million years ago
Optic nerve crush
Optic nerve injury
Principal component analysis
Peripheral nervous system
Spinal cord injury
- X. laevis :
Abe N, Borson SH, Gambello MJ, Wang F, Cavalli V. Mammalian target of rapamycin (mTOR) activation increases axonal growth capacity of injured peripheral nerves. J Biol Chem. 2010;285:28034–43.
Agostinone J, Di Polo A. Retinal ganglion cell dendrite pathology and synapse loss: implications for glaucoma. Prog Brain Res. 2015;220:199–216.
Ananthakrishnan L, Gervasi C, Szaro BG. Dynamic regulation of middle neurofilament RNA pools during optic nerve regeneration. Neuroscience. 2008;153:144–53.
Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11:R106.
Andrews S. FASTQC: a quality control tool for high throughput sequence data. www.bioinformatics.babraham.ac.uk/projects/fastqc. 2010.
Anthony TE, Mason HA, Gridley T, Fishell G, Heintz N. Brain lipid-binding protein is a direct target of notch signaling in radial glial cells. Genes Dev. 2005;19:1028–33.
Arimoto A, Tagawa K. Regeneration in the enteropneust hemichordate, Ptychdera flava, and its evolutionary implications. Develop Growth Differ. 2018;60:400–8.
Avci HX, Lebrun C, Wehrle R, Doulazmi M, Chatonnet F, Morei MP, Ema M, Vodjdani G, Sotelo C, Filamant F, Dusart I. Thyroid hormone triggers the developmental loss of axonal regenerative capacity via thyroid hormone receptor α1 and krüppel-like factor 9 in Purkinje cells. Proc Natl Acad Sci U S A. 2012;109:14206–11.
Avraham Y, Davidi N, Lassri V, Vorobiev L, Kabesa M, Dayan M, Chernoguz D, Berry E, Leker RR. Leptin induces neuroprotection, neurogenesis, and angiogenesis after stroke. Curr Neurovasc Res. 2011;8:313–22.
Baumgart M, Groth M, Priebe S, Savino A, Testa G, Dix A, Ripa R, Spallotta F, Gaetano C, Ori M, Terzibasi Tozzini E, Guthke R, Platzer M, Cellerino A. RNA-seq of the aging brain in the short-lived fish N. furzeri - conserved pathways and novel genes associated with neurogenesis. Aging Cell. 2014;13:965–74.
Beattie MS, Bresnahan JC, Lopate G. Metamorphosis alters the response to spinal cord transection in Xenopus laevis frogs. J Neurobiol. 1990;21:1108–22.
Beaver RS, Dunlop SA, Harman AM, Stirling RV, Easter SS Jr, Roberts JD, Beazley LD. Continued neurogenesis is not a pre-requisite for regeneration of a topographic retino-tectal projection. Vis Res. 2001;41:1765–70.
Bender MC, Sifuentes CJ, Denver RJ. Leptin induces mitosis and activates the canonical Wnt/b-catenin signaling pathway in neurogenic regions of Xenopus brain. Front Endocrinol (Lausanne). 2017;8:99.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Royal Stat Soc Ser B. 1995;57:289–300.
Benowitz LI, He Z, Goldberg JL. Reaching the brain: advances in optic nerve regeneration. Exp Neurol. 2017;287(Pt 3):365–73.
Benowitz LI, Popovich PG. Inflammation and axon regeneration. Curr Opin Neurol. 2011;24:577–83.
Bu W, Su LK. Regulation of microtubule assembly by human EB1 family proteins. Oncogene. 2001;20:3185–92.
Campfield LA, Smith FJ, Guisez Y, Devos R, Burn P. Recombinant mouse OB protein: evidence for a peripheral signal linking adiposity and central neural networks. Science. 1995;269:546–9.
Cañete-Soler R, Reddy KS, Tolan DR, Zhai J. Aldolases a and C are ribonucleolytic components of a neuronal complex that regulates the stability of the light neurofilament mRNA. J Neurosci. 2005;25:4353–64.
Chan JY, Hsieh TY, Liu ST, Chou WY, Chung MH, Huang SM. Physical and functional interactions between hnRNP K and PRMT family proteins. FEBS Lett. 2009;583:281–6.
Chang J, Baker J, Wills A. Transcriptional dynamics of tail regeneration in Xenopus tropicalis. Genesis. 2017;55:e23015.
Cholia RP, Nayyar H, Kumar R, Mantha AK. Understanding the multifaceted role of ectonucleotide pyrophosphatase/phosphodiesterase 2 (ENPP2) and its altered behaviour in human diseases. Curr Mol Med. 2015;15:932–43.
Chou J, Ferris AC, Chen T, Seok R, Yoon D, Suzuki Y. Roles of polycomb group proteins enhancer of zeste E(z) and Polycomb (pc) during metamorphosis and larval leg regeneration in the flour beetle Tribolium castaneum. Dev Biol. 2019;450:34–46.
Chung MI, Kwon T, Tu F, Brooks ER, Gupta R, Meyer M, Baker JC, Marcotte EM, Wallingford JB. Coordinated genomic control of ciliogenesis and cell movement by RFX2. eLife. 2014;3:e01439.
Conesa A, Madrigal P, Tarazona S, Gomez-Cabrero D, Cervera A, McPherson A, Szczesniak MW, Gaffney DJ, Elo L, Zhang X, Mortazavi A. A survey of best practices for RNA-seq data analysis. Genome Biol. 2016;17:13–9.
Deem JD, Muta K, Ogimoto K, Nelson JT, Velasco KR, Kaiyala KJ, Morton GJ. Leptin regulation of core body temperature involves mechanisms independent of the thyroid axis. Am J Physiol Endocrinol Metab. 2018;315:E552–64.
Dobreva G, Braun T. The yin and yang of polycomb repression in regenerating muscle. Cell Stem Cell. 2010;7:422–4.
Dovey M, Patton EE, Bowman T, North T, Goessing W, Zhou Y, Zon LI. Topoisomerase II alpha is required for embryonic development and liver regeneration in zebrafish. Mol Cell Biol. 2018;29:3746–53.
Drew K, Lee C, Huizar RL, Tu F, Borgeson B, McWhite CD, Ma Y, Wallingford JB, Marcotte EM. Integration of over 9,000 mass spectrometry experiments builds a global map of human protein complexes. Mol Systems Biol. 2017;13:932.
Dutta B, Ren Y, Hao P, Sim KH, Cheow E, Adav S, Tam JP, Sze SK. Profiling of the chromatin-associated proteome identifies HP1BP3 as a novel regulator of cell cycle progression. Mol Cell Proteomics. 2014;13:2183–97.
Dutta B, Ren Y, Lim SK, Tam JP, Sze SK. Quantitative profiling of chromatome dynamics reveals a novel role for HP1BP3 in hypoxia-induced oncogenesis. Mol Cell Proteomics. 2014;13:3236–49.
Elsaeidi F, Bemben MA, Zhao XF, Goldman D. Jak/stat signaling stimulates zebrafish optic nerve regeneration and overcomes the inhibitory actions of Socs3 and Sfpq1. J Neurosci. 2014;34:2632–44.
Fawcett JW, Gaze RM. The organization of regenerating axons in the Xenopus optic nerve. Brain Res. 1981;229:487–90.
Fernandes KA, Harder JM, Kim J, Libby RT. JUN regulates early transcriptional responses to axonal injury in retinal ganglion cells. Exp Eye Res. 2013;112:106–17.
Fernandez-Godino R, Bujakowska KM, Pierce EA. Changes in extracellular matrix cause RPE cells to make basal deposits and activate the alternative complement pathway. Hum Mol Genet. 2018;27:147–59.
Forehand CJ, Farel PB. Anatomical and behavioral recovery from the effects of spinal cord transection: dependence on metamorphosis in anuran larvae. J Neurosci. 1982;2:654–62.
Friesen WJ, Dreyfuss G. Specific sequences of the Sm and Sm-like (Lsm) proteins mediate their functions with the spinal muscular atrophy disease gene product (SMN). J Biol Chem. 2000;275:26370–5.
Galmozzi E, Aghemo A, Colombo M. Eukaryotic initiation factor 5B: a new player for the anti-hepatitis C virus effect of ribavirin? Med Hypotheses. 2012;79:471–3.
Gaze RM. Regeneration of the optic nerve in Xenopus laevis. Q J Exp Physiol Cogn Med Sci. 1959;44:290–308.
Gaze RM, Jacobson M. A study of the retino-tectal projection during regeneration of the optic nerve in the frog. Proc R Soc Lond B. 1963;157:420–48.
Gervasi C, Thyagarajan A, Szaro BG. Increased expression of multiple neurofilament mRNAs during regeneration of vertebrate central nervous system axons. J Comp Neurol. 2003;461:262–75.
Gibbs KM, Chittur SV, Szaro BG. Metamorphosis and the regenerative capacity of spinal cord axons in Xenopus laevis. Eur J Neurosci. 2011;33:9–25.
Gibbs KM, Szaro BG. Regeneration of descending projections in Xenopus laevis tadpole spinal cord demonstrated by retrograde double labeling. Brain Res. 2006;1088:68–72.
Gibbs KM, Szaro BG. Tracing Central Nervous System Axon Regeneration in Xenopus. Cold Spring Harb Protoc. 2018. https://doi.org/10.1101/pdb.prot101030.
Goodbrand IA, Gaze RM. Microglia in tadpoles of Xenopus laevis: normal distribution and the response to optic nerve injury. Anat Embryol. 1991;184:71–82.
Hall D, Poussin C, Velaqapudi VR, Empsen C, Joffraud M, Beckmann JS, Geerts AE, Ravussin Y, Ifferson M, Oresic M, Thorens B. Peroxisomal and microsomal lipid pathways associated with resistance to hepatic steatosis and reduced pro-inflammatory state. J Biol Chem. 2010;285:31011–23.
Hart T, Komori HK, LaMere S, Podshivalova K, Salomon DR. Finding the active genes in deep RNA-seq gene expression studies. BMC Genomics. 2013;14:778.
Herman PE, Papatheodorou A, Bryant SA, Waterbury CKM, Herdy JR, Arcese AA, Buxbaum JD, Smith JJ, Morgan JR, Bloom O. Highly conserved molecular pathways, including Wnt signaling, promote functional recovery from spinal cord injury in lampreys. Sci Rep. 2018;8:742.
Herwartz C, Castillo-Juarez P, Schröder L, Barron BL, Steger G. The transcription factor ZNF395 is required for the maximal hypoxic induction of proinflammatory cytokines in U87-MG cells. Mediat Inflamm. 2015;2015:804264.
Hui SP, Sengupta D, Lee SG, Sen T, Kundu S, Mathavan S, Ghosh S. Genome wide expression profiling during spinal cord regeneration identifies comprehensive cellular responses in zebrafish. PLoS One. 2014;9:e84212.
Hutchins EJ, Szaro BG. C-Jun N-terminal kinase phosphorylation of heterogeneous nuclear ribonucleoprotein K regulates vertebrate axon outgrowth via a posttranscriptional mechanism. J Neurosci. 2013;33:14666–80.
Inés M, Rehermann I, Marichal N, Russo RE, Trujillo-Cenóz O. Neural reconnection in the transected spinal cord of the freshwater turtle Trachemys dorbignyi. J Comp Neurol. 2009;515:197–214.
Ivashkiv LB, Zhao B, Park-Min KH, Takami M. Feedback inhibition of osteoclastogenesis during inflammation by IL-10, M-CSF receptor shedding, and induction of IRF8. Ann N Y Acad Sci. 2011;1237:88–94.
Jin D, Liu Y, Sun F, Wang X, Liu X, He Z. Restoration of skilled locomotion by sprouting corticospinal axons induced by co-deletion of PTEN and SOCS3. Nat Commun. 2015;6:8074.
Kakebeen AD, Chitsazan D, Williams MC, Saunders LM, Wills AE. Chromatin accessibility dynamics and single cell RNA-Seq reveal new regulators of regeneration in neural progenitors. eLife. 2020;9:e52648.
Kamei N, Kwon SM, Ishikawa M, Li M, Nakanishi K, Yamada K, Hozumi K, Kawamoto A, Ochi M, Asahara T. Endothelial progenitor cells promote astrogliosis following spinal cord injury through Jagged1-dependent notch signaling. J Neurotrauma. 2012;29:1758–69.
Kanehisa M, Furumichi M, Tanabe M, Sato Y, Morishima K. KEGG: new perspectives on genomes, pathways, diseases, and drugs. Nucl Acids Res. 2017;45:D353–61.
Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucl Acids Res. 2000;28:27–30.
Kanehisa M, Sato Y, Furumichi M, Tanabe M. New approach for understanding genome variations in KEGG. Nucl Acids Res. 2019;47:D590–5.
Kasmer SH, Rodriguez D, De Tomaso AW. Colonial ascidians as model organisms for the study of germ cells, fertility, whole body regeneration, vascular biology and aging. Curr Opin Genet Dev. 2016;39:101–6.
Kath C, Paloma GO, Muller R, Schultz C, Haucke V, Eickholt B, Schmoranzer J. PTEN suppresses axon outgrowth by down-regulating the level of detyrosinated microtubules. PLoS One. 2018;13:e0193257.
King MW, Neff AW, Mescher AL. Proteomics analysis of regenerating amphibian limbs: changes during the onset of regeneration. Int J Dev Biol. 2009;53:955–69.
Knowlden S, Georas SN. The autotaxin-LPA axis emerges as a novel regulator of lymphocyte homing and differentiation. J Immunol. 2014;192:851–7.
Lacerda R, Menezes J, Romão L. More than just scanning: the importance of cap-independent mRNA translation initiation for cellular stress response and cancer. Cell Mol Life Sci. 2017;74:1659–80.
Lang DM, Rubin BP, Schwab ME, Stuermer CAO. CNS myelin and oligodendrocytes of the Xenopus spinal cord - but not optic nerve - are nonpermissive for axon growth. J Neurosci. 1995;15:99–109.
Langmead B, Salzberg S. Fast gapped-read alignment with bowtie 2. Nat Methods. 2012;9:357–9.
Leclercq IA, Field J, Farrell GC. Leptin-specific mechanisms for impaired liver regeneration in ob/ob mice after toxic injury. Gastroenterology. 2003;124:1451–64.
Lee-Liu D, Edwards-Faret G, Tapia VS, Larrain J. Spinal cord regeneration: lessons for mammals from non-mammalian vertebrates. Genesis. 2013;51:529–44.
Lee-Liu D, Mendez-Olivos EE, Muñoz R, Larrain J. The African clawed frog Xenopus laevis: a model organism to study regeneration of the central nervous system. Neurosci Lett. 2016;652:82–93.
Lee-Liu D, Moreno M, Almonacid LI, Tapia VS, Muñoz R, von Marees J, Gaete M, Melo F, Larrain J. Genome-wide expression profile of the response to spinal cord injury in Xenopus laevis reveals extensive differences between regenerative and non-regenerative stages. Neural Dev. 2014;9:12.
Li X, Tian R, Gao R, Yang Y, BRG W, Gantier MP, McMillan NAJ, Xu D, Hu Y, Gao Y. Identification of a histone family gene signature for predicting the prognosis of cervical cancer patients. Sci Rep. 2017;7:16495.
Liao D. Emerging roles of the EBF family of transcription factors in tumor suppression. Mol Cancer Res. 2009;7:1893–901.
Lin G, Chen Y, Slack JMW. Imparting regenerative capacity to limbs by progenitor cell transplantation. Dev Cell. 2013;24:41–51.
Liu L, Aldskogius H, Svensson M. Ultrastructural localization of immunoglobulin G and complement c9 in the brain stem and spinal cord following peripheral nerve injury: an immunoelectron microscopic study. J Neurocytol. 1998;27:737–48.
Liu X, Li C, Zhang R, Xiao W, Niu X, Ye X, Li Z, Guo Y, Tan J, Li Y. The EZH2- H3K27me3- DNMT1 complex orchestrates epigenetic silencing of the wwc1 gene, a hippo/yap pathway upstream effector, in breast cancer epithelial cells. Cell Signal. 2018;51:243–56.
Liu Y, Yu H, Deaton SK, Szaro BG. Heterogeneous ribonucleoprotein K, an RNA-binding protein, is required for optic axon regeneration in Xenopus laevis. J Neurosci. 2012;32:3563–74.
Lonati C, Carlin A, Leonardi P, Valenza F, Bosari S, Catania A, Gatti S. Modulatory effects of NDP-MSH in the regenerating liver after partial hepatectomy in rats. Peptides. 2013;50:145–52.
Long Y, Li L, Li Q, He X, Cui Z. Transcriptomic characterization of temperature stress responses in larval zebrafish. PLoS One. 2013;7:e37209.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.
Love NR, Chen Y, Bonev B, Gilchrist MJ, Fairclough L, Lea R, Mohun TJ, Paredes R, Zeef LA, Amaya E. Genome-wide analysis of gene expression during Xenopus tropicalis tadpole tail regeneration. BMC Dev Biol. 2011;11:70.
Ma KH, Duong P, Moran JJ, Janaidi N, Svaren J. Polycomb repression regulates Schwann cell proliferation and axon regeneration after nerve injury. Glia. 2018;66:2487–502.
Ma XM, Blenis J. Molecular mechanisms of mTOR-mediated translational control. Nature Rev Mol Cell Biol. 2009;10:307–18.
Maki N, Kimura H. Epigenetics and regeneration. Curr Top Microbiol Immunol. 2013;367:237–52.
Martinez-De Luna RI, Ku RY, Aruck AM, Santiago F, Viczian AS, San Mauro D, Zuber ME. Müller glia reactivity follows retinal injury despite the absence of the glial fibrillary acidic protein gene in Xenopus. Dev Biol. 2017;426:219–35.
McEwan J, Lynch J, Beck CW. Expression of key retinoic acid modulating genes suggests active regulation during development and regeneration of the amphibian limb. Dev Dyn. 2011;240:1259–70.
Minoura I. Towards an understanding of the isotype-specific functions of tubulin in neurons: technical advances in tubulin expression and purification. Neurosci Res. 2017;122:1–8.
Monaghan JR, Walker JA, Page RB, Putta S, Beachy CK, Voss SR. Early gene expression during natural spinal cord regeneration in the salamander Ambystoma mexicanum. J Neurochem. 2007;101:27–40.
Moody SA, Miller V, Spanos A, Frankfurter A. Developmental expression of a neuron-specific beta-tubulin in frog (Xenopus laevis): a marker for growing axons during the embryonic period. J Comp Neurol. 1996;364:219–30.
Munoz J, Dhillon N, Janku F, Watowich SS, Hong DS. STAT3 inhibitors: finding a home in lymphoma and leukemia. Oncologist. 2014;19:536–44.
Myung K, Smith S, Kolodner RD. Mitotic checkpoint function in the formation of gross chromosomal rearrangements in Saccharomyces cerevisiae. Proc Natl Acad Sci U S A. 2004;101:15980–5.
Nguyen MH, Cheng M, Koh TJ. Impaired muscle regeneration in Ob/Ob and db/db/ mice. Sci World J. 2011;11:1525–35.
Nguyen VQ, Co C, Irie K, Li JJ. Clb/Cdc28 kinases promote nuclear export of the replication initiator proteins Mcm2-7. Curr Biol. 2000;10:195–205.
Nieuwkoop PD, Faber J. Normal table of Xenopus laevis (Daudin). New York: Garland publishing; 1994.
Noor NM, Steer DL, Wheaton BJ, Ek CJ, Truettner JS, Dietrich WD, Dziegielewska KM, Richardson SJ, Smith AI, BandeBerg JL, Saunders NR. Age-dependent changes in the proteome following complete spinal cord transection in a postnatal south American opossum (Monodelphis domestica). PLoS One. 2011;6:e27465.
Norsworthy MW, Bei F, Kawaguchi R, Sanes JR, Coppola G, He Z. Sox11 expression promotes regeneration of some retinal ganglion cell types but kills others. Neuron. 2017;94:1112–20.
Ohno S. Evolution by gene duplication New York: springer; 1970.
Oliva AA Jr, Atkins CM, Copenagle L, Banker GA. Activated c-Jun N-terminal kinase is required for axon formation. J Neurosci. 2006;26:9462–70.
Ostberg A, Norden J. Ultrastructural study of degeneration and regeneration in the amphibian tectum. Brain Res. 1979;168:441–55.
Park J, Kim S, Joh J, Remick SC, Miller DM, Yan J, Kanaan Z, Chao JH, Krem MM, Basu SK, Hagiwara S, Kenner L, Moriggl R, Bunting KD, Tse W. MLLT11/AF1q boosts oncogenic STAT3 activity through Src-PDGFR tyrosine kinase signaling. Oncotarget. 2018;2016:43960–73.
Park KK, Liu K, Hu Y, Smith PD, Wang C, Cai B, Xu B, Connolly L, Kramvis I, Sahin M, He Z. Promoting axon regeneration in the adult CNS by modulation of the PTEN/mTOR pathway. Science. 2008;322:963–6.
Park SH, Yu M, Kim J, Moon Y. C/EBP homologous protein promotes NSAID-activated gene 1-linked pro-inflammatory signals and enterocyte invasion by enteropathogenic Escherichia col. Microbes Infect. 2017;19:110–21.
Pelzer D, Phipps LS, Thuret R, Baker SM, Dorey K. Foxm1 regulates neuronal progenitor fate during spinal cord regeneration. bioRxiv. 2020. https://doi.org/10.1101/2020.02.26.962977v1.
Perucho L, Artero-Castro A, Guerrero S, Ramon y Cajal S, Lleonart ME, Wang ZQ. RPLP1, a crucial ribosomal protein for embryonic development of the nervous system. PLoS One. 2014;9:e99956.
Peterson SL, Anderson AJ. Complement and spinal cord injury: traditional and non-traditional aspects of complement cascade function in the injured spinal cord microenvironment. Exp Neurol. 2014;258:35–47.
Phipps LS, Marshall L, Dorey K, Amaya E. Model systems for regeneraiton: Xenopus. Development. 2020;147:dev180844.
Popov IK, Kwon T, Crossman DK, Crowley MR, Wallingford JB, Chang C. Identification of new regulators of embryonic patterning and morphogenesis in Xenopus gastrulae by RNA sequencing. Dev Biol. 2017;426:429–41.
Priscilla R, Szaro BG. Comparisons of SOCS mRNA and protein levels in Xenopus provide insights into optic nerve regenerative success. Brain Res. 1704;2019:150–60.
Quigley IK, Kintner C. Rfx2 stabilizes Foxj1 binding at chromatin loops to enable multiciliated cell gene expression. PLoS Genet. 2017;13:e1006538.
Rao N, Song F, Jhamb D, Wang M, Milner D, Price NM, Belecky-Adams TL, Palakai MJ, Cameron JA, Li B, Chen X, Stocum DL. Proteomic analysis of fibroblastema formation in regenerating hind limbs of Xenopus laevis froglets and comparison to axolotl. BMC Dev Biol. 2014;14:32.
Rasmussen KD, Helin K. Role of TET enzymes in DNA methylation, development, and cancer. Genes Dev. 2016;30:733–50.
Raymond PA, Barthel LK, Bernardos RL, Perkowski JJ. Molecular characterization of retinal stem cells and their niches in adult zebrafish. BMC Dev Biol. 2006;6:36.
Rensen SS, Slaats Y, Driessen A, Peutz-Koostra CJ, Nijhuis J, Steffensen R, Greve JW, Buurman WA. Activation of the complement system in human nonalcoholic fatty liver disease. Hepatology. 2009;50:1809–17.
Roberts A, Pimentel H, Trapnell C, Pachter L. Identification of novel transcripts in annotated genomes using RNA-Seq. Bionformatics. 2011;27:2325–9.
Roberts A, Trapnell C, Donaghey J, Rinn JL, Pachter L. Improving RNA-Seq expression estimates by correcting for fragment bias. Genome Biol. 2011;12:R22.
Robinson JT, Thorvaldsdottir H, Winckler W, Guttman M, Lander ES, Getz G, Mesirov JP. Integrative genomics viewer. Nature Biotech. 2011;29:24–6.
Roosa JR, Gervasi C, Szaro BG. Structure, biological activity of the upstream regulatory sequence, and conserved domains of a middle molecular mass neurofilament gene of Xenopus laevis. Mol Brain Res. 2000;82:35–51.
Rousseau E, Michel PP, Hirsch EC. The iron-binding protein lactoferrin protects vulnerable dopamine neurons from degeneration by preserving mitochondrial calcium homeostasis. Mol Pharmacol. 2013;84:888–98.
Sakaguchi T, Okada M, Kitamura T, Kawasaki K. Reduced diameter and conduction velocity of myelinated fibers in the sciatic nerve of a neurofilament-deficient mutant quail. Neurosci Lett. 1993;153:65–8.
Sandgren J, Holm S, Marino AM, Asmundsson J, Grillner P, Nister M, Diaz de Stahl T. Whole exome- and mRNA-sequencing of an AT/RT case reveals few somatic mutations and several deregulated signalling pathways in the context of SMARCB1 deficiency. Biomed Res Int. 2015;2015:862039.
Schwartz ML, Shneidman PS, Bruce J, Schlaepfer WW. Stabilization of neurofilament transcripts during postnatal development. Mol Brain Res. 1994;27:215–20.
Session AM, Uno Y, Kwon T, Chapman JA, Toyoda A, Takahashi S, Fukui A, Hisosaka A, Suzuki A, Kondo M, Van Heeringen SI, Quigley I, Heinz S, Ongino H, Ochi H, Hellsten U, Lyons JB, Simakov O, Putnam N, Stites J, Kruoki Y, Tanaka T, Michiue T, Watanabe M, Bogdanovic O, Lister R, Georgiou G, Paranjpe SS, van Kruijsbergen I, Shu S, Carlson J, Konoshita T, Ohta Y, Mawaribuchi S, Jenkins J, Grimwood J, Schmutz J, Mitros T, Mozaffari SV, Suzuki Y, Haramoto Y, Yamamoto TS, Takagi C, Heald R, Miller K, Haudenshicl C, Kitzman J, Nakayama T, Izutsu Y, Robert J, Fortriede J, Burns K, Lotay V, Karimi K, Yasuoka Y, Dichmann DS, Flajnik MS, Houston DW, Shendure J, DuPasquier L, Vize PD, Zorn AM, Ito M, Marcotte EM, Wallingford JB, Tio Y, Asashima M, Ueno N, Matsuda Y, Veenstra GJ, Fujiyama A, Harland RM, Taira M, Rokhsar DS. Genome evolution in the allotetraploid frog Xenopus laevis. Nature. 2016;538:336–43.
Shaw T, Martin P. Epigenetic reprogramming during wound healing: loss of polycomb-mediated silencing may enable upregulation of repair genes. EMBO Rep. 2009;10:881–6.
Smith PD, Sun F, Park KK, Cai B, Wang C, Kuwako K, Martinez-Carrasco I, Connolly L, He Z. SOCS3 deletion promotes optic nerve regeneration in vivo. Neuron. 2009;64:617–23.
Song W, Cho Y, Watt D, Cavalli V. Tubulin-tyrosine ligase (TTL)-mediated increase in tyrosinated α-tubulin in injured axons is required for retrograde injury signaling and axon regeneration. J Biol Chem. 2015;290:13344–53.
Sperry RW. Optic nerve regeneration with return of vision in anurans. J Neurophysiol. 1944;7:351–61.
Szaro BG, Faulkner LA, Hunt RK, Loh YP. Axonal transport of [35S] methionine labeled proteins in Xenopus optic nerve: phases of transport and the effects of nerve crush on protein patterns. Brain Res. 1984;297:337–5.
Szaro BG, Loh YP. Changes in axonal transport and glial proteins during optic nerve regeneration in Xenopus laevis. Curr Top Dev Biol. 1987;21:217–54.
Szaro BG, Loh YP, Hunt RK. Specific changes in axonally transported proteins during regeneration of the frog (Xenopus laevis) optic nerve. J Neurosci. 1985;5:192–208.
Szklarczyk D, Franceschini A, Wyder S, Forslund K, Heller D, Huerta-Cepas J, Simonovic M, Roth A, Santos A, Tsafou KP, Kuhn M, Bork P, Jensen LJ, von Mering C. STRING v10: protein-protein interaction networks, integrated over the tree of life. Nucl Acids Res. 2015;43(Database issue):D447–52.
Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, Simonovic M, Doncheva NT, Morris JH, Bork P, Jensen LJ, Mering CV. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucl Acids Res. 2019;47(D1):D607–13.
Szklarczyk D, Morris JH, Cook H, Kuhn M, Wyder S, Simonovic M, Santos A, Doncheva NT, Roth A, Bork P, Jensen LJ, von Mering C. The STRING database in 2017: quality-controlled protein-protein association networks, made broadly accessible. Nucleic Acids Res. 2017;45(D1):D362–8.
Tang YP, Wade J. Sexually dimorphic expression of the genes encoding ribosomal proteins L17 and L37 in the song control nuclei of juvenile zebra finches. Bran Res. 2006;1126:102–8.
Tao Y, Wang M, Chen E, Tang H. Liver regeneration: analysis of the main relevant signaling molecules. Mediat Inflamm. 2017;2017:4256352.
Tapia VS, Herrera-Rojas M, Larrain J. JAK/STAT pathway activation in response to spinal cord injury in regenerative and non-regenerative stages of Xenopus laevis. Regeneration (Oxf). 2017;4:21–35.
Taylor JSH, Jack JL, Easter SS Jr. Is the capacity for optic nerve regeneration related to continued retinal ganglion cell production in the frog? A test of the hypothesis that neurogenesis and axon regeneration are obligatorily linked. Eur J Neurosci. 1989;1:626–38.
Thorvaldsdottir H, Robinson JT, Mesirov JP. Integrative genomics viewer (IGV): high-performance genomics data visualization and exploration. Brief Bioinform. 2013;14:178–92.
Trapnell C, Hendrickson DG, Sauvageau M, Goff L, Rinn JL, Pachter L. Differential analysis of gene regulation at transcript resolution with RNA-seq. Nature Biotech. 2013;31:46–54.
Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009;25:1105–11.
Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, Pimentel H, Salzberg SL, Rinn JL, Pachter L. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and cufflinks. Nat Protoc. 2012;7:562–78.
Trapnell C, Williams B, Pertea G, Mortazavi A, Kwan G, van Baren J, Salzberg S, Wold B, Pachter L. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nature Biotech. 2010;28:511–5.
Tripathi S, Pohl MO, Zhou Y, Rodriguez-Frandsen A, Wang G, Stein DA, Moulton HM, DeJesus P, Che J, Mulder LC, Yanguez E, Anbdenmatten D, Pache L, et al. Meta- and orthogonal integration of influenza “OMICS” data defines a role for UBR4 in virus budding. Cell Host Microbe. 2015;18:723–35.
Valentin-Kahan A, Garcia-Tejedor GB, Robello C, Trujillo-Cenoz O, Russo RE, Alvarez-Valain F. Gene expression profiling in the injured spinal cord of Trachemys scripta elegans: an amniote with self-repair capabilities. Front Mol Neurosci. 2017;10:17.
van Mier P, ten Donkelaar HJ. Early development of descending pathways from the brain stem to the spinal cord in Xenopus laevis. Anat Embryol. 1984;170:295–306.
Veldman MB, Bemberi MA, Goldman D. Tuba1a gene expression is regulated by KLF7/7 and is necessary for CNS development and regeneration in zebrafish. Mol Cell Neurosci. 2010;43:370–83.
Wagatsuma A. Adipogenic potential can be activated during muscle regeneration. Mol Cell Biochem. 2007;304:25–33.
Walker KL, Yoo H-K, Undamatla J, Szaro BG. Loss of neurofilaments alters axonal growth dynamics. J Neurosci. 2001;21:9655–66.
Wang GX, Cho KW, Uhm M, Hu CR, Li S, Cozacov Z, Xu AE, Cheng JX, Saltiel AR, Lumenq CN, Lin JD. Otopetrin 1 protects mice from obesity-associated metabolic dysfunction through attenuating adipose tissue inflammation. Diabetes. 2014;63:1340–52.
Wang Z, Long H, Chang C, Zhao M, Lu Q. Crosstalk between metabolism and epigenetic modifications in autoimmune diseases: a comprehensive review. Cell Mol Life Sci. 2018;75(18):3353–69.
Wang Z, Reynolds A, Kirry A, Nienhaus C, Blackmore MG. Overexpression of Sox11 promotes corticospinal tract regeneration after spinal injury while interfering with functional recovery. J Neurosci. 2015;35:3139–45.
Weng YL, An R, Cassin J, Joseph J, Mi R, Wang C, Zhong C, Jin SG, Pfeifer GP, Bellacosa A, Dong X, Hoke A, He Z, Song H, Ming GL. An intrinsic epigenetic barrier for functional axon regeneration. Neuron. 2017;94:337–46.
Whitworth GB, Misaghi BC, Rosenthal DM, Mills EA, Heinen DJ, Watson AH, Ives CW, Ali SH, Bezold K, Marsh-Armstrong N, Watson FL. Translational profiling of retinal ganglion cell optic nerve regeneration in Xenopus laevis. Dev Biol. 2017;426:360–73.
Wilson MA, Gaze RM, Goodbrand LA, Taylor JSH. Regeneration in the Xenopus laevis tadpole optic nerve is preceded by a massive macrophage/microglial response. Anat Embryol. 1992;186:75–89.
Wong P, Coulombe PA. Loss of keratin 6 (K6) proteins reveals a function for intermediate filaments during wound repair. J Cell Biol. 2003;163:327–37.
Yang Y, Zhang J, Wu T, Xu X, Cao G, Li H, Chen X. Histone deacetylase 2 regulates the doxorubicin (Dox) resistance of hepatocarcinoma cells and transcription of ABCB1. Life Sci. 2019;216:200–6.
Yoshida M. Oligodendrocyte maturation in Xenopus laevis. J Neurosci Res. 1997;50:169–76.
Yoshida M, Zhao L, Grigoryan G, Shim H, He P, Yun CC. Deletion of Na+/H+ exchanger regulatory factor 2 represses colon cancer progress by suppression of Stat3 and CD24. Am J Physiol Gastrointest Liver Physiol. 2016;10:G576–98.
Yu Y, Chi B, Xia W, Gangopadhyay J, Yamazaki T, Winkelbauer-Hurt ME, Yin S, Eliasse Y, Adams E, Shaw CE, Reed R. U1 snRNP is mislocalized in ALS patient fibroblasts bearing NLS mutations in FUS and is required for motor neuron outgrowth in zebrafish. Nucleic Acids Res. 2015;43:3208–18.
Yu Y, Schachner M. Syntenin-a promotes spinal cord regeneration following injury in adult zebrafish. Eur J Neurosci. 2013;38:2280–9.
Zhao XF, Wan J, Powell C, Ramachandran R, Myers MG Jr, Goldman D. Leptin and IL-6 family cytokines synergize to stimulate Müller glia reprogramming and retina regeneration. Cell Rep. 2014;9:272–84.
Zhao Y, Szaro BG. The return of phosphorylated and nonphosphorylated epitopes of neurofilament proteins to the regenerating optic nerve of Xenopus laevis. J Comp Neurol. 1994;343:158–72.
The authors thank Dr. Christine Gervasi (University at Albany, Albany, NY) for animal rearing support, and Michael Monforte (Schenectady, NY) and Dr. Sergei Reverdatto (Univ. at Albany) for technical and bioinformatics support. Part of this work was performed with the assistance of the University of Louisville Genomics Facility, which is supported by NIH P20GM103436 (KY IDeA Networks of Biomedical Research Excellence), NIH P30GM106396 (U of L JG Brown Cancer Center Phase III CoBRE), the JG Brown Foundation, and user fees. We also thank the instructors and staff of the 2016 Xenopus Bioinformatics workshop (National Xenopus Resource Center, Marine Biological Laboratories, Woods Hole, MA; RRID:SCR_013731) for guidance in conducting RNA-seq analysis in Xenopus.
Financial support came from the New York Department of Health Spinal Cord Injury Research Board (NY DOH SCIRB C32091GG, C30837GG, and C32249GG) and the National Science Foundation (IOS 1257449) to B.G.S., and the National Institutes of Health to K.M.G. (NIH 1R15HD076643 and GM103436) and M.A.S. (NIH R15GM128049). These agencies played no role in either the design of the study, the collection, analysis, and interpretation of the data, or the writing of the manuscript.
Ethics approval and consent to participate
All animal procedures were approved by the Institutional Animal Care and Use Committees (IACUC) of the University at Albany (optic nerve crush) or Morehead State University (spinal cord transection) in accordance with the National Institutes of Health Guide for the Care and Use of Laboratory Animals.
Consent for publication
All authors report no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Cuffdiff2 differential expression data. Tabs are for each injury condition and time point vs. their relevant controls (N = 3 biological replicates, each). Tabs PW1–3) Tadpole spinal cord injury (SCI) hindbrain vs. age-matched, unoperated control hindbrain at 3 days, 1 week, and 3 weeks, respectively; PW4–6) Juvenile frog SCI hindbrain at 3 days, 1 week, and 3 weeks vs. unoperated control hindbrain; PW7–9) juvenile frog optic nerve crush (ONC) eye vs. the contralateral unoperated eye at 3 days, 11 days, and 3 weeks, respectively. Abbreviations: chr, chromosome; cntrl, control; FPKM, fragments per kilobase of transcript per million mapped reads; ONC, optic nerve crush; p, the calculated probability for a single given gene that the observed difference represented no change in expression (i.e., log2(fold change) = 0); q, false discovery rate-adjusted (FDR) p-value ; significant, “yes” if q < 0.05; Tx, spinal cord-transected; Unop, unoperated.
DESeq2 differential expression data. Tabs are for each injury condition and time point vs. their relevant controls (N = 3 biological replicates, each). Tabs PW1–3) Tadpole spinal cord injury (SCI) hindbrain vs. age-matched, unoperated control hindbrain at 3 days, 1 week, and 3 weeks, respectively; PW4–6) Juvenile frog SCI hindbrain at 3 days, 1 week, and 3 weeks vs. unoperated control hindbrain; PW7–9) juvenile frog optic nerve crush (ONC) eye vs. the contralateral unoperated eye at 3 days, 11 days, and 3 weeks, respectively. Columns: A - Gene Name; B, baseMean - mean counts of all samples, normalized for sequencing depth. It does not take into account gene length and is used for estimating the dispersion of a gene; C, log2FoldChange - log2(control/experimental); in this analysis, negative values represent injury induced increases relative to control); D, lfcSE - Standard Error of the log2 fold change; E, stat - log2FoldChange/lfcSE; F, pvalue - the probability of the null hypothesis for an individual gene (i.e., variance among replicates = variance between conditions); G, padj - false discovery rate-adjusted (FDR) p-value.
Temporal patterns of gene expression and shared injury-response genes between regenerative vs. non-regenerative tissues, as determined by DESeq2. A) Regenerative tissues [i.e., SCI tadpole hindbrain (SCI Tadpole) and ONC juvenile frog eye (ONC Juvenile)] shared similar temporal patterns of numbers of significant (FDR < 0.05) differentially expressed genes, which differed markedly from that of the non-regenerative tissue [SCI juvenile frog hindbrain (SCI Juvenile)]. Whereas the expression response of the two regenerative tissues peaked during the mid recovery phase (1 week/11 days), that of the non-regenerative tissue peaked at the early, post trauma phase (3 days). Up- and down-regulated genes are shown in green and red, respectively; S & L gene homeologs were tallied separately. B) Plot illustrating the percentage of annotated genes that were significantly (FDR < 0.05) differentially expressed with injury (100% = 24,382 genes). Additional_File2_Differential_Expression_Analysis_by_DESeq2.xlsm contains the DESeq2 output files from which A and B were derived. C - E) UpSet plots illustrating the number of genes overlapping between the samples indicated by the circles below each bar at 3 days (C), 7/11 days (D), and 3 weeks (E) after injury. Numbers of shared up- and down-regulated genes are indicated above and below each bar, respectively. The maximum number of overlapping genes between the two successfully regenerative tissues (DESR: Differentially Expressed in Successful Regeneration) occurred during the peak phase of regenerative CNS axon outgrowth.
Manually Curated DESR Gene Data (Cuffdiff2). DESR gene lists were manually curated and parsed into eleven functional categories using a literature search. Tabs: 3 Day DESR, 11 Day 1 WK DESR, 3 WK DESR: Manually curated spreadsheets for each time point, 3 days, 11 days/1 week, and 3 weeks, respectively, providing Cuffdiff2 data for each gene for tadpole SCI hindbrain (even numbered rows) and juvenile frog ONC (odd numbered rows). Column A, surgical condition; Columns B - N, relevant Cuffdiff2 data for each gene under each surgical condition (see Additional_File1_Differential_Expression_Analysis_by_Cuffdif.xlsm for abbreviations and column annotations). Gene symbols (Column B) in green and red were up-regulated and down-regulated, respectively, in successfully regenerative tissues relative to controls. Column O, primary functional category (NA, gene has not been assigned a human ortholog); Column P, secondary functional category, if relevant; Column Q, Manually curated notes, providing additional information on the expression behaviors of the other homeologs of this gene (even numbered rows), and a statement on the function of the gene that lead to its being placed in a particular functional category, followed by the literature reference supporting it.
Scatterplot representation of the Principal Component Analyses (PCA) of gene expression profiles. Ellipses group biological replicates for each indicated condition (experimentals, solid squares; controls, empty squares), indicating variation among samples. A-C, PCA of tadpole and juvenile hindbrain after spinal cord injury (SCI), and of juvenile frog after optic nerve crush (ONC), respectively. In C, expression profiles from the operated eye were compared with those of the contralateral, unoperated eyes within the same animals. D, PCA of all 17 conditions combined, supporting the tissue-specific nature of gene expression profiles. Conditions were the same as in A-C, except that data from eyes of animals receiving no injury was included (open triangles, Frog Eye, Unop). E, PCA of tadpole and juvenile frog hindbrain samples after spinal cord injury, supporting the conclusion that the differences in gene expression profiles between the time points at which numbers of differentially expressed genes reached their peaks (3 days in juvenile frog hindbrain and 1 week in tadpole hindbrain) were more than just a kinetic difference in the timing of expression of the same differentially expressed genes. Abbreviations: ONC, optic nerve crush; PC1, principal component axis 1; PC2, Principal Component axis 2; SCI, spinal cord injured; TX, spinal cord transected; Unop ONC - unoperated eye, contralateral to the operated eye; Frog Eye, Unop - eyes from unoperated animals; wk., week.
A summary list of all DESR genes, sorted by primary functional category and time point of differential expression, along with relevant references supporting the functional assignment. Note that genes were assigned their human gene ortholog symbols, without reference to whether they were L or S homeologs, or paralogs.
KEGG pathway analysis of the DESR gene lists revealed that 60% (69/115) of the Adipocytokine signaling pathway genes were differentially expressed in at least one tissue at one time point (p < 0.02). Tab. 1, Pathway: the Adipocytokine Signaling Pathway (obtained from KEGG) and the complete associated gene list; Tab. 2, Gene Expression Levels: Column B provides the list of genes and their synonyms that are in this KEGG pathway (the underlined name is the one provided in the figure in Tab. 1. Column C provides the gene symbol and tissue for each gene. ΔFPKM [log2(fold change), injury condition/control, Column D] for each differentially expressed gene in the pathway (115 genes, S & L homeologs listed separately), as well as for other isoforms for each gene. Columns D – L, the ΔFPKMs for each tissue at all three time points are listed. For significant changes, either by p or q (or both), and the raw FPKMs are listed in a cell comment for the ΔFPKMs. The meanings of the different color highlights is given in the legend at the top of the table. Tab. 3, Genes Not Found: is a list of 10 genes (including synonyms) from the KEGG adipocytokine signaling pathway that were not found in the Mayball annotations of the X. laevis transcriptome. Abbreviations are as in Additional_File1_Differential_Expression_Analysis_by_Cuffdif.xlsm.
About this article
Cite this article
Belrose, J.L., Prasad, A., Sammons, M.A. et al. Comparative gene expression profiling between optic nerve and spinal cord injury in Xenopus laevis reveals a core set of genes inherent in successful regeneration of vertebrate central nervous system axons. BMC Genomics 21, 540 (2020). https://doi.org/10.1186/s12864-020-06954-8
- Xenopus laevis
- Spinal cord injury
- Optic nerve injury
- Axon regeneration
- Central nervous system