Open Access Highly Accessed Research article

Gene expression profiling of intestinal regeneration in the sea cucumber

Pablo A Ortiz-Pineda1, Francisco Ramírez-Gómez1, Judit Pérez-Ortiz1, Sebastián González-Díaz1, Francisco Santiago-De Jesús1, Josue Hernández-Pasos1, Cristina Del Valle-Avila1, Carmencita Rojas-Cartagena1, Edna C Suárez-Castillo2, Karen Tossas1, Ana T Méndez-Merced3, José L Roig-López1, Humberto Ortiz-Zuazaga4 and José E García-Arrarás1*

Author Affiliations

1 University of Puerto Rico, Rio Piedras, Department of Biology, San Juan, PR, USA

2 Yale University, Department of Molecular, Cellular and Developmental Biology, New Haven, CT, USA

3 Universidad del Este- SUAGM, School of Science and Technology, Carolina, PR, USA

4 University of Puerto Rico, HPCf, High Performance Computing facility, San Juan, PR, USA

For all author emails, please log on.

BMC Genomics 2009, 10:262  doi:10.1186/1471-2164-10-262

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2164/10/262


Received:23 January 2009
Accepted:8 June 2009
Published:8 June 2009

© 2009 Ortiz-Pineda et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Among deuterostomes, the regenerative potential is maximally expressed in echinoderms, animals that can quickly replace most injured organs. In particular, sea cucumbers are excellent models for studying organ regeneration since they regenerate their digestive tract after evisceration. However, echinoderms have been sidelined in modern regeneration studies partially because of the lack of genome-wide profiling approaches afforded by modern genomic tools.

For the last decade, our laboratory has been using the sea cucumber Holothuria glaberrima to dissect the cellular and molecular events that allow for such amazing regenerative processes. We have already established an EST database obtained from cDNA libraries of normal and regenerating intestine at two different regeneration stages. This database now has over 7000 sequences.

Results

In the present work we used a custom-made microchip from Agilent with 60-mer probes for these ESTs, to determine the gene expression profile during intestinal regeneration. Here we compared the expression profile of animals at three different intestinal regeneration stages (3-, 7- and 14-days post evisceration) against the profile from normal (uneviscerated) intestines. The number of differentially expressed probes ranged from 70% at p < 0.05 to 39% at p < 0.001. Clustering analyses show specific profiles of expression for early (first week) and late (second week) regeneration stages. We used semiquantitative reverse transcriptase polymerase chain reaction (RT-PCR) to validate the expression profile of fifteen microarray detected differentially expressed genes which resulted in over 86% concordance between both techniques. Most of the differentially expressed ESTs showed no clear similarity to sequences in the databases and might represent novel genes associated with regeneration. However, other ESTs were similar to genes known to be involved in regeneration-related processes, wound healing, cell proliferation, differentiation, morphological plasticity, cell survival, stress response, immune challenge, and neoplastic transformation. Among those that have been validated, cytoskeletal genes, such as actins, and developmental genes, such as Wnt and Hox genes, show interesting expression profiles during regeneration.

Conclusion

Our findings set the base for future studies into the molecular basis of intestinal regeneration. Moreover, it advances the use of echinoderms in regenerative biology, animals that because of their amazing properties and their key evolutionary position, might provide important clues to the genetic basis of regenerative processes.

Background

All living organisms exhibit, to some extent, regenerative properties that allow them to deal with environmental events, physical trauma or diseases. Regenerative capacities have been studied in terms of stem cell recruitment, cell dedifferentiation, proliferation and migration, provision of specific regulatory/trophic factors, and expression or re-expression of the developmental program in adult animals [1]. However, the molecular/genetic basis of regeneration remains obscure.

Central to regeneration studies is the choice of the model organism, since metazoan species can show large variability in their regenerative capacities. In general, the ability to replace complex body parts decreases as one moves from the basal to the more highly derived taxa. However, even within the same phylum, not all animals are able to regenerate body parts, and not all tissues within a body can be equally repaired [2,3]. One problem in elucidating the molecular basis of regeneration has been that organisms with high regenerative capacities do not lend themselves easily to traditional experimental genetics. Only recently has progress been made to make possible the use of some of these model systems to dissect the genetic basis of regeneration [4]. Among these, invertebrate systems such as Hydra and planaria, and vertebrates, such as ascidians and amphibian urodeles, have gained particular attention. Nonetheless, regeneration research remains underpopulated, and there are whole phyla of organisms, showing very interesting regenerative behaviors where little molecular research has been performed [5]

Among deuterostomes, the regenerative potential is maximally expressed in echinoderms, animals that can quickly replace most injured organs. In particular, sea cucumbers (holothurians) are excellent models for studying organ regeneration since they can regenerate many of their organs and appendages. Our laboratory has been using an echinoderm, Holothuria glaberrima as a model organism for the last decade to dissect the process of how the digestive tract regenerates once it is eliminated by auto-evisceration. We have shown that following evisceration the new intestine regenerates from the free end of the remaining mesentery [6]. Initially, the mesenterial tip thickens forming a continuous rod-like structure that extends from the esophagus to the cloaca. During the second regeneration week, luminal epithelial cells from the esophagus and from the cloaca migrate into this tube, forming the mucosal layer and giving rise to the intestinal lumen. We have performed extensive studies at the cellular level showing the involvement of cell division, dedifferentiation, and migration in the regeneration of the intestine [6-11] as well as of events associated with the remodeling of the extracellular matrix [12].

We have also used the sea cucumber to explore the role of the genes that allow for such extraordinary regenerative processes [13]. Our approach has been to focus on target genes that have been associated with regenerative processes or identified in the regenerating tissues [14,15]. Similar gene by gene approaches have been used to study regeneration processes in other echinoderms, particularly in brittle stars and crinoids [16-20]. However, the regenerative capacities of the echinoderms have yet to be explored systematically using a large number of molecular tools [21]. In fact, one of the reasons that echinoderms have been sidelined in modern regeneration studies is the lack of genome-wide profiling approaches afforded by modern genomics [21]

We have now overcome this problem by determining the profile of gene activity during intestinal regeneration in H. glaberrima using microarray technology. For this, microarray slides were made using over 7000 ESTs that we previously identified from normal and regenerating intestine cDNA libraries [13]. Here we compare the gene expression profile of animals at 3, 7 and 14 days of regeneration, following the evisceration process typical of these organisms. Results show a large number of differentially expressed genes associated with intestinal regeneration. Some of these genes are homologues to metazoan genes associated with regenerative processes while many others might be novel sequences with little or no similarities to sequences in the databanks. To explore our results, we have focused on those sequences that are differentially expressed during early intestinal regeneration and that are associated with events previously shown to occur at the cellular level, namely developmental processes, cytoskeletal transformations and extracellular matrix remodeling. Our results make available the first overview of gene expression patterns in the regenerating intestine. More importantly, our results provide the basis for molecular studies aimed at exploring the molecular basis of regeneration in a group of animals that because of their remarkable properties and their key evolutionary position will provide important clues to the genetic basis of regenerative processes.

Results

Technical Analyses

Technical analysis showed that the microarrays performed as expected. First, the intensities of the 536 internal controls of the Agilent array perfectly matched the expected ones as shown in the regression model [see additional file 1]. Second, no significant variation was observed between raw data and normalized MA plots [Additional file 2], which means that there was no significant dye bias. Third, dye swap results showed extremely high similarity in their labeling pattern [additional file 3]. Since the dye swaps were made with different animal samples (biological dye-swaps), this implies that not only did dyes produced similar results, but that results from animals at the same regeneration stage were equivalent. Fourth, analyses of genes spotted in duplicate within the same array (technical replicates) showed almost identical fold changes, evidencing the strong reproducibility of our data.

Additional file 1. Microarray internal controls performance. Numerical and graphical representation of the internal control signals.

Format: ZIP Size: 172KB Download fileOpen Data

Additional file 2. Normalization procedure. Example of two random selected arrays before and after normalization.

Format: ZIP Size: 371KB Download fileOpen Data

Additional file 3. Dye-swap hybridization. Bottom left corner of one array and the same corner in the corresponding dye-swap microarray

Format: ZIP Size: 249KB Download fileOpen Data

Gene expression profiles

A. Venn Diagrams

Microarray results showed dramatic differences in gene expression between regenerating and normal intestine and also among animals at different regeneration stages. Of the 14352 H. glaberrima probes in our microarray, 5915 showed similar expression in all groups, including normal (uneviscerated) and 3-, 7- and 14-days post-evisceration (dpe) regenerating intestines.

Most of the differential expression was observed when normal and regenerating intestines were compared. A large number of probes, 8437 (aprox 58%) showed differential expression at a significance level of p < 0.01 (Figure 1). As expected, the number of differentially expressed spots at p < 0.05 was larger (73%) but even at p < 0.001 a significant number of genes (39%) still were found to be differentially expressed (Figure 2). Differential gene expression was also observed when comparisons were made among regenerating stages. The largest differences were observed in probes that changed specifically in the 3-dpe stage, whether they were differentially expressed only at 3-dpe (2347), at both 3- and 7-dpe (3950) or at 3-, 7- and 14-dpe (1190) stages compared to normal. A smaller number of probes was differentially expressed only at 7-dpe (835).

thumbnailFigure 1. Venn Diagram of the EST distribution. 5915 probes do not show significant changes in expression when compared with normal tissues. The highest differential expression was found in both 3-dpe and 7-dpe stages (4439 differentially expressed spots); p < 0.01.

thumbnailFigure 2. Percentage of differentially expressed sequences at p < 0.001, 0.01 and 0.05. At p < 0.001 only 39% of the sequences were significantly up or down regulated. P-values are shown in logarithmic scale.

B. Volcano Plots

In order to have a graphical representation of the gene expression levels and their significance, we analyzed our results using volcano plots (Figure 3).

thumbnailFigure 3. Microarrays volcano plots. Distribution of sequence expression between normal and normal injected with LPS, 3-, 7- and 14-dpe regenerating stages. Each plot shows the logarithm of the probability of the t-test as a function of the logarithm of fold change for each EST probe. Horizontal lines in each plot represent the nominal significant level of 0.001 for the t-student under the assumption that each gene has a unique variance. The vertical lines represent the limit of significance of the change in expression (fold change >2). Differentially expressed genes are located on the right (over-expressed) and left (under-expressed) top quadrants.

Results showed not only that a large number of genes are differentially expressed, but a high level of change and significance, mainly during the first week of regeneration. For example, at 3-dpe and 7-dpe the range of odds was significantly larger than at 14-dpe. Comparison with LPS showed a small number of genes being differentially expressed.

The top differentially expressed genes for 3-, 7- and 14-dpe stages of regeneration are shown on Tables 1, 2 and 3, respectively. Given that many probes have high significant p-values, we selected the top genes by two criteria: First, those probes with the highest p-value. Second, the fold change of the ratio of the intensities. The sequences were named either with their accession numbers (singlets) or with local identifications for multisequence contigs. All nucleotide sequences have been submitted to NCBI (taxonomy ID: 31192). Genes are shown in order of fold change and probability.

Table 1. Top 10 differentially expressed genes at 3dpe

Table 2. Top 10 differentially expressed genes at 7dpe.

Table 3. Top 10 differentially expressed genes at 14dpe.

Among the genes in this top 10 list that show known identities, were cyclophilin, laminin, alpha tubulin, tenascin-R, ferritin. murinoglobulin and NM23. Interestingly, most of the differentially expressed genes showed no significant similarity (e-value < 1.00E-07) to other sequences in publicly available nucleotide or protein databases. Unknown genes, such as Contig 4863-1 and 5300-1, were up-regulated at 3-, 7- and 14-dpe, while the singlet ES728532 showed down-regulation at all regeneration stages. On the other hand, some genes were specific for one stage such as Contig 3933-1 and Contig 4702-1 at 3-dpe. However, these genes were significantly differentially expressed but with higher p-values at other time points.

An alternate strategy to focus on gene groups that are differentially expressed at particular regeneration stages is to compare gene expression among regeneration stages (3-, 7- and 14dpe) and not directly to the normal intestine. To do this, we used the expression profile at 7-dpe as reference (up-regulated genes at 7dpe are considered as those whose expression is down-regulated at normal, 3 and 14dpe- stages; contrary, down-regulated genes are those whose expression is up-regulated at normal, 3 and 14 dpe). Seventeen genes (corresponding to 26 ESTs represented by 55 probes in the microarray) were found to be up-regulated at 3-dpe when compared with 7dpe (Table 4). The number of genes with up-regulated expression at 7- and 14-dpe stages was 24 and 13, respectively (Table 5 and 6).

Table 4. Upregulated genes specific to 3dpe

Table 5. Upregulated genes specific to 7dpe

Table 6. Upregulated genes specific to 14dpe

C. Validation of Microarray Results by RT-PCR

Gene-specific relative RT-PCRs were performed to validate the expression profiles of fifteen genes identified in the microarray as differentially expressed. The sequences of the primers used are shown on Table S1. Validation included genes with putative identity (Wnt-9, Tensc-R, MMP-15, MMP-11, MMP-14, Actin 1, Actin 2, Actin 3 and Hox12). Most of these genes have been previously associated with wound-healing, development and regeneration in other animal species. The second group consisted of genes without identity (putative novel genes-Unknowns E, F and G, Contig 4874-1, Contig 5242-1 and Contig 5501-1). The validated genes include over-regulated and down-regulated sequences of the microarray.

The expression profiles of the 15 genes that were validated using RT-PCR confirmed the robustness of the microarray results (Figure 4 and 5). Eleven out of the 15 genes showed exact correlation in their expression profiles between RT-PCR and microarrays and very similar significance values in the statistical analyses. For example, unknown F showed up-regulation at 3-dpe and decreased gradually at 7- and 14-dpe. The unknowns 5501, 4874 and 5242 showed an exact match between the microarray and the PCR results. On the other hand, the correlation was not perfect for four genes (actin 1, 2 and 3, and Hox 12). These genes showed some correlations at particular stages, such as observed for actin 1 and 2 and Hox 12 at 3-dpe. At other stages some of the RT-PCR results although not significant do showed a trend that resembled the microarray results (see actin 3 at 3-dpe). Overall the correlation of RT-PCR and microarray results was 84%, with only 7 of the 45 validated time points not showing the same profile between the microarrays and the RT-PCR. A direct comparison of the RT-PCR and the microarray results is shown on Table 7.

Table 7. Significance obtained from microarray and PCR at 3, 7 and 14dpe

thumbnailFigure 4. RT-PCR profiles of Wnt-9, Hox12, Tenascin-R, MMP-11, MMP-14, MMP-15, Actin 1, Actin 2 and Actin 3 in normal and regenerating intestines. The values were normalized against NADH. This set of genes showed up-regulation in the first week of regeneration except for Actin 3 that showed down-regulation during regeneration. Interestingly, all extracellular matrix (ECM) related genes showed the peak of expression at 7-dpe. Moreover, Wnt-9, a developmental gene, was up-regulated during all stages. Values represent the mean ± S.E.M. of 7 biological samples (individual animals) per stage. Paired t-test for mean comparison was used. Significance levels were *:p < 0.05, and **:p < 0.01. Images at the base of the bar are RT-PCR pictures.

thumbnailFigure 5. RT-PCR profiles of 6 novel gene sequences in normal and regenerating intestines. The values were normalized against NADH. This set of genes was mainly differentially expressed in the first week of regeneration, except by unknown-G that keeps down-regulated at 14dpe. Three of the four up-regulated genes showed the peak of expression at 3dpe and decrease the expression until 14dpe. Only unknown 5242-1 showed the peak at 7dpe but even at 3dpe the up-regulation was highly significant when compared with normal tissues. Values represent the mean ± S.E.M. of 7 samples per stage. Paired t-test for mean comparison was used. Significance levels were: *:p < 0.05, and **: p < 0.01. Images in the base of the bar are RT-PCR pictures

D. Cluster Analysis

Cluster analysis was performed on the microarray results using hierarchical (not shown) and non-hierarchical methods in order to reveal several distinct profiles of expression (Figure 6).

thumbnailFigure 6. Cluster analysis of the microarray expression profiles for 15208 probes (All minus 536 controls). The 5 columns show the following comparisons: LPS (Red: "R") vs. Normal (Green: "G"): 3dpe(R) vs. Normal(G); 7dpe(R) vs. Normal(G); 7dpe(R) vs. 3dpe(G); 14dpe(R) vs. 7dpe(G). The number of genes in each cluster is shown above the bars. Since arrays were normalized against the LPS vs Normal results, the first column appears empty. Colors represent the level of expression of the gene, red for up-regulated and green for down-regulated; White spots represent the average of all the genes per cluster with "y" axis as a relative measure of expression levels.

Cluster # A contains probes that were highly over-expressed at 3- and 7-dpe when compared to normal intestine, and whose expression levels are slightly higher in 3- than in 7-dpe and are further decreased at 14-dpe.

Cluster # B contains probes that were moderately over-expressed at 3-dpe when compared to normal and whose expression levels decreased at 7-dpe. The expression levels of these probes at 14-dpe were similar to those at 7-dpe.

Cluster # C contains probes that were over-expressed at 3-dpe but whose expression peaked at 7-dpe, followed by a decrease in expression levels at 14-dpe.

Cluster # D contains probes that showed no significant differences in expression at any stage when compared to normal intestine.

Cluster # E contains probes that were slightly under-expressed at 3-dpe when compared to normal. The expression level of these probes was higher at 7- than at 3-dpe and somewhat similar between 7- and 14-dpe.

Cluster # F contains probes that were slightly under-expressed at 3- and 7-dpe when compared to normal and whose expression levels increase at 14-dpe.

Cluster # G contains probes that were moderately under-expressed at 3-dpe when compared to normal and whose expression levels slightly increase at 7- and 14-dpe.

Cluster # H contains probes that were highly under-expressed at 3- and 7-dpe and whose expression levels showed a moderate increase at 14-dpe.

E. Functional category analysis

We submitted all sea cucumber sequences to BlastX and parsed the results using a threshold of significance e-value < 1.00E-05. Approximately 58.1% of the sequences (3651) were "unknowns" that showed no similarity to subjects in the SwissProt database. On the other hand, 41.9% of the sequences (2636) had similarity with database sequences. Of these, 169 showed similarity with hypothetical proteins, most of them from the sea urchin S. purpuratus. GO searches http://www.geneontology.org webcite were performed with all the blast results and 1514 unique sequences were found with at least one reported ontology (Figure 7A).

thumbnailFigure 7. Distribution of gene function (GO) for ESTs analyzed in the microarrays. GO values of all ESTs printed on the microarray (A). Percentage of ESTs differentially expressed at 3-dpe (B) or 7-dpe (C). When compared with all the GOs, the 3- and 7-dpe distribution GOs showed no significant differences other than a small decrease in the number of GO results. The number of GOs results for each particular stage is shown in the bottom-left corner. The color code starts at 12 o'clock in the list order.

In order to determine if particular gene functions could be associated exclusively with specific regeneration stages, diagrams were done with only those genes that showed differential expression between 3-, 7- or 14-dpe and normal. At 3-dpe the GO search showed 1231 results while 1068 were found for 7-dpe (Figure 7B and 7C). The distribution of the GO was highly similar between all sequences and regeneration stages. For example, in all cases sequences associated with metabolic and biosynthetic processes corresponded to half of the GO results.

F. Individual gene expression

In addition to the global view of gene expression provided by microarray analyses, individual gene expression patterns can be determined. We have focused on various genes that have been associated with regenerative processes. For this we only chose those genes whose similarity to genes from other animals was sufficiently high (e-value < 1.00E-07) to be confident that they correspond to holothurian homologues. These genes are classified and presented in three main groups in relation with their function: a) Developmental genes, b) Cytoskeletal-related genes and c) Extracellular matrix-related genes.

a). Developmental genes

Many developmental gene pathways have been found to occur or be activated during regenerative processes [22,23]. Thus, we have focused on the expression profile of development-associated genes to explore their roles during intestinal regeneration. Table 8 shows the expression profile of well-known developmental genes in the microarray.

Table 8. Expression profile of developmental genes in the sea cucumber H. glaberrima during regeneration.

In general, differential expression of these genes mainly occurred at early stages of regeneration. Most of the genes are activated or inhibited at 3-dpe (8 genes, being 5 of them highly significant) and 7-dpe (6 genes, being 5 highly significant). In contrast, at 14-dpe only two developmental genes were differentially expressed, both just at p < 0.05. Hox genes stand out among the group of development-associated genes, and our group has previously identified them in the holothurians [24]. Four of these genes appeared to be over-expressed, although their expression profile differed. Hox 9, 10 and 12 were over-expressed at 3- and 7-dpe and returned to normal levels at 14-dpe. However, while Hox 10 displayed the same levels of over-expression at 3- and 7-dpe, Hox 12 and 9 showed peaks of expression at 3-dpe and 7-dpe, respectively. Hox 5 showed a different pattern of expression altogether being over-expressed only at 14-dpe. Other development-associated genes were also found to be differentially expressed in our microarrays during intestinal regeneration. Once again the expression profile differed among them. For example, Wnt 14 was found to be over-expressed at 3-, 7- and 14-dpe, while BMP-1 was only up-regulated at 7-dpe. In contrast, forkhead box K1 was down regulated at 3-dpe and, similarly, myotrophin was down-regulated at 3- and 7-dpe. Finally, not all development-associated genes showed changes in expression. For example, Kruppel-like factor 2 remained at similar levels of expression during the first two weeks of regeneration.

To validate the developmental gene results we chose to study the expression of Hox 12 and Wnt 9 by RT-PCR. The sequences for these genes were obtained from the 7-dpe library and the complete sequence of Hox 12 has been previously characterized by our laboratory [25] RT-PCR showed that genes were over expressed during early regeneration and their levels decreased during the second week of regeneration.

b). Cytoskeletal Proteins

Regeneration induces cellular changes that have been previously reported by our group. Among these are processes of dedifferentiation and migration that are associated with the formation of the new intestinal primordium during early regeneration [7-9]. Such processes are associated with changes in the expression of cytoskeletal proteins. Thus, it is expected that these changes would be reflected in the expression of mRNAs coding for cytoskeletal proteins. We have focused on the expression of seven cytoskeletal genes found within our EST database (Table 9). These include three actin family members, four tubulin family members, two myosin genes and one gelsolin gene. Although actin and tubulin isoforms are highly similar at the nucleotide level of the open reading frame (ORF), we identified ours based on differences in their 3'-unstranslated (UTR) region. The expression profile of these cytoskeletal mRNAs differed widely. Some were over-expressed during the first 2 weeks of regeneration (Actin 1, Actin 2, Tubulin alfa and alfa-1), others were under-expressed (Actin 3, Myosin-11 and gelsolin precursor) while others showed no significant changes in their expression (tubulin-beta). Of particular interest were those sequences that showed a specific change at 3-dpe and appeared to only change slightly or not at all at other regeneration stages (tubulin 3 and myosin light 3F chain), since these are more likely to play a role in specific processes occurring during intestinal regeneration.

Table 9. Expression profile of cytoskeletal genes in the sea cucumber H. glaberrima during regeneration.

Some of these sequences were validated by RT-PCR. In particular, the differential expression profiles of actins 1 and 2 that are over-expressed during early regeneration stages, while actin 3 is down-regulated, was clearly documented in our RT-PCRs.

c). ECM genes

We have previously described dramatic remodeling of the ECM and the involvement of MMPs during intestinal regeneration [12]. In this respect we searched for sequences in our database that might be associated with the ECM changes. Seven such genes were found, four of them corresponding to ECM molecules (Echinonectin, Collagen alfa-1, Laminin alpha and Tenascin-R) and three corresponding to matrix metalloproteases (MMPs) (15, 11, and 14) (Table 10). In this group, all analyzed genes showed up-regulation when compared with normal tissues. The mRNAs for the four ECM molecules showed over-expression during the first two regeneration weeks. Different to ECM molecules, MMPs were over-expressed at 3- and 7-dpe but not at 14-dpe, with MMP 15 showing low significance at 3- but high at 7-dpe. The four genes from this group that were validated (Figure 5) showed strong correlation between RT-PCR and microarray results.

Table 10. Expression profile of ECM-related genes in the sea cucumber H. glaberrima during regeneration.

d). Novel sequences

As mentioned earlier many of the sequences in our microarray slides were unknowns. That is, sequences that showed no similarities to those in databanks. In many cases this was due to the sequence being either too short or being part of the mRNA's UTR. However, in some cases the sequence was long enough to contain a significant ORF that, in principle, should have shown similarities to other proteins when submitted to BLAST searches. These are the most interesting sequences to our research group, particularly those that were differentially expressed during regeneration, since they might correspond to novel holothurian-specific sequences associated with their striking regenerative capacities.

One strategy that we explored to study these unknowns was to focus on contigs. That is, those ESTs that shared overlapping sequence identity among them, and that together formed a composed sequence. The number of ESTs from a specific library (stage of regeneration) found to form a contig can be correlated with the expression level for that sequence at the stage from where the cDNA library was made. Therefore, an indirect measurement of the level of gene expression can be inferred from the number of identified ESTs from each library that formed contigs. We identified 596 multi-sequence contigs in our EST databank that had at least two assembled ESTs. Of these, 85.4% had ORFs larger than 100 nt (~33aa), and 46.5% had ORFs larger than 195 nt (~65aa). Only 27 contigs had no ORFs made up of at least 30 nt. We focused on contigs that showed over-expression at 3- and 7-dpe that were composed by at least three different ESTs and where all the EST probes showed similar results in the microarray. Sixteen of these unknown contigs are shown on Table 11. Their putative ORFs ranged from 98 to 794 nucleotides in length and they were represented by three to twenty-nine different ESTs. As expected, down-regulated contigs at 3-dpe appeared in the H and E clusters, while the up-regulated contigs at 3- and 7-dpe were located mostly in the A and B clusters. Moreover, some of these genes were validated with RT-PCR (Figure 6) showing high consistency in their expression profiles (Table 3).

Table 11. Novel genes differentially expressed during regeneration in H. glaberrima.

e). Other sequences

Our microarray chip included the sequences of 75 ESTs from another sea cucumber A. japonicus, and 329 sequences from genes of the sea urchin S. purpuratus. Only five of the sea cucumber (Collagen Alpha-1 (CA1); Cytoskeleton actin (cAct); Fibrinogen B precursor (Fib-Bp), Heat Shock protein 90b (HSP-90) and Senescence associated protein (SAP)) showed significant differences in expression in regenerating animals. Two of A. japonicus differentially expressed ESTs were underexpressed during regeneration (Fib-Bp and SAP), while the remaining three were overexpressed. Interestingly, four were similar to the H. glaberrima ESTs and followed the same expression profile (Fib-Bp was the different one). In all cases the probes of these other organisms showed similarity with the C-terminal region of the H. glaberrima sequences. The only sequence from the sea urchin that changed during regeneration was a hypothetical protein (Glean3-07946) showing a decrease in its expression at 3- and 7-dpe. The 24 sequences from zebra fish, human, rat, mice frog and axolotl behaved as expected; none showed any significant labeling or differential expression.

Discussion

We have now shown by using microarray technology that a large number of genes, represented by EST sequences, are differentially expressed during sea cucumber intestinal regeneration. Nonetheless, microarray results are not fool-proof and must be validated by other methods. An indirect way of validating our results is to compare this microarray data to previous work from our laboratory showing individual gene expression patterns. Two examples of these are the expression of serum amyloid A (SAA) [14] and of ependymin [15], both of which have been well studied in the laboratory using a range of molecular techniques. SAA was shown by Northern blot analysis to be over-expressed during regeneration, with a peak during the second week[14]. Our microarray results also show SAA over-expression during intestinal regeneration, although the expression peak somewhat differs, occurring during the first week. This difference might be due to the fact that SAA is also associated with the immune status of the animals [26] therefore this status might differ between the animals used in the present experiments from those used for the Northern experiments performed years ago. On the other hand, the expression of ependymin, a gene associated with regeneration and plasticity, was previously shown to increase during intestinal regeneration [15]. These results were obtained with the highly accurate technique of quantitative Real Time PCR, showing a peak in expression at 7-dpe and no change at 3- or 14-dpe. These results are closely reproduced in our microarray analyses which show significant differences between 7-dpe animals and uneviscerated, but no differences at 3- or 14-dpe.

The comparison of RT-PCR results and microarray data obtained for fifteen selected ESTs provides a direct corroboration of their differential expression during intestinal regeneration. These results show that 11/15 of the sequences behave essentially the same whether measured by microarray or RT-PCR. Moreover, even the statistical significance of the result is within the same ballpark figure. The differences between RT-PCR and microarray results of the four remaining sequences (Hox 12, actins 1, 2 and 3), is mainly due to the RT-PCR results not being statistically significant, although in some cases they follow the same trend. It is interesting that these sequences are among the ones known to have closely related sequences in the genome (multiple actin genes and Hox genes), thus suggesting that the lack of correlation at some stages might be due to primers hybridizing with other sequences and increasing the background noise, making the RT-PCR results less sensitive than the microarray.

• Why do so many genes show differential expression?

One of the most striking results from our microarray experiments, is that a high percentage of the spotted ESTs were found to be differentially expressed. Even when considering a very high level of significance (p < 0.001) more than a third of the ESTs showed differential expression. This contrasts with other regeneration studies where a smaller number of genes were found. For example, in studies of zebrafish heart regeneration, only about 5% of the assayed genes were found to be differentially regulated [27] and between 8 and 14% in studies of liver regeneration in rodents [28,29].

The reasons for these differences might be due to two factors. First, most of our ESTs were obtained from cDNA libraries of regenerating intestine, thus the original pool of ESTs is probably biased toward those genes that are highly expressed in regenerating tissues. This was confirmed by the results themselves: the set of probes differentially expressed at 3dpe and 7dpe were 81.25% from the 3- and 7-dpe libraries and only 18.75% from the Normal library. In contrast, 42% of the down-regulated sequences come from the Normal library. Second, the process of intestinal regeneration that we are studying is complex and involves organogenesis, not merely wound healing. Furthermore, the intestine is an organ with multiple cell types and embryological origins. The intestine has components of the three embryonic germ layers, the enteric nervous system being an ectodermal derivative, the luminal epithelium an endodermal derivative and the muscular layer and submucosa being mesodermal derivatives. Therefore, when compared to a process such as heart or liver regeneration, where fewer cell types or more limited regeneration is involved, a larger number of genes might be necessary to be modulated to achieve the formation of the new intestine. Moreover, in the study of zebra fish heart, regeneration was limited to the ventricle and only to a portion of the organ since only 20% of the ventricle was removed. This contrasts largely from our studies in the sea cucumber where the new intestine must be completely formed from the remaining mesentery. In fact, a large number of differentially expressed genes was also found in another organism with high regenerative potential, the salamander Ambyostoma mexicanum. In this species when the tail, including the spinal cord segment is amputated, 76% of the genes in the microarray were differentially expressed [30]. It is interesting that they also found that even when increasing the stringency of the statistical analyses the number of differentially expressed genes (35%) was still significant. Thus, regeneration in deuterostomes imply the upregulation and downregulation of a large number of genes, making the search for "master regeneration genes", much more difficult than anticipated.

In view of this difficulty, we have developed certain strategies aimed at singling out those genes that might play important or crucial roles in regeneration. First, is to focus on genes that are differentially expressed at the 3-dpe stage. Previous experiments from our laboratory have shown that the cellular processes associated with intestinal regeneration occur early in intestinal regeneration. By 7-dpe, the blastema-like structure that will be the primordium for the new intestine is already formed [6]. Thus, genes associated with initiating regeneration must be regulated earlier than at 7-dpe. In this respect, some of the most likely candidates will be those genes that are differentially expressed between the 3- and 7-dpe stages. Second, focus on those genes that are most highly differentially expressed, such as those that were identified on our Top ten analyses. Third, look at sequences that form contigs. Although the presence of multiple contigs might be explained by several factors (some of them technique-related), they could be reflecting the abundance of gene expression in the original cDNA libraries, and therefore at the stage from where the library was made. Fourth, is to aim at those genes that have been previously associated with organogenesis in other species (see the following section). Finally, we expect that novel genes associated with regeneration will be characterized from our data. Thus, those sequences, and in particular those contigs that show ORFs with no similarity to sequences in databases are the strongest candidates to be pursued.

• How do microarray results correlate with organogenesis processes that occur during embryonic development?

Developmental Genes

The finding of genes associated with embryonic development being differentially expressed during regeneration is not surprising since developmental gene pathways have been found to occur or be activated during regenerative processes [22,23]. Thus, these genes and signaling pathways provide excellent targets to study the molecular events of intestinal regeneration.

One of the most studied groups of genes is the Hox genes, which encode transcription factors that specify the body axis during embryonic development. Hox genes have been implicated in various regeneration processes, including the regeneration of amphibian tails and limbs [31-33], hydra [34], rat liver [35], and zebrafish tail [36] among others. This expression can be spatially or temporally specific in regenerating tissues [37,38]. It is particularly interesting that Hox12 had previously been found in our regenerating intestine cDNA library and that there is some evidence for homeobox gene involvement in echinoderm regeneration [20]. The finding that the most posterior genes (Hox 9, 10 and 12) are the ones that showed an increase in expression during intestinal regeneration coincides with findings in the mammalian digestive tract where Hox genes are known to be expressed in a rostral-caudal gradient, and the Hox genes near the 5' end are expressed in the posterior digestive tract, i.e. intestine [39]. These are the same genes that when mutated cause defects in the development of the posterior digestive tract [40,41].

Similarly, the finding that a Wnt homologue (Wnt-9) was overexpressed in the regenerating intestine during the first two weeks of regeneration coincides with published data on this molecular family. Wnt molecules are extracellular molecules that bind to their receptors and activate a signaling pathway that includes other known proteins, such as B-catenin. In recent years, Wnt pathways have been increasingly associated with regenerative phenomena. Wnt was found to be involved in blastema formation of the regenerating limbs of anuran tadpoles [42] and in lens regeneration in newts [43]. In mammals, Wnt has been studied in bone [44], hair follicle [45] and deer antler regeneration [46] among others. Liver regeneration was retarded in the absence of beta-catenin [47]. Wnt apparently plays a key role in the control of intestinal stem cell proliferation and differentiation [48]. Thus, it will be of great interest when the holothurian Wnt-9 expression is localized to determine which intestinal cell types are expressing this protein. Wnt pathways have also been involved in other invertebrate regeneration models. In planaria, Wnt is necessary for proper brain pattern formation [49] and B-catenin for antero-posterior axis formation during regeneration [50].

Lastly, the family of bone morphogenetic proteins (BMPs) has been increasingly associated with regenerative phenomena. In newts, BMPs have been associated with the potential of the ventral iris to regenerate a new lens [51]. And in zebrafish, at least three members of the BMP family have been associated with fin regeneration [52]. Inhibition of BMP activity using the BMP antagonist, noggin, caused a reduction in blastema cell proliferation and a reduction in bone matrix deposition. Similarly, in experiments using transgenic Xenopus tadpole with increase noggin expression, both tail and limb regeneration were inhibited [53]. Finally, BMP has been shown to be essential for the correct establishment of the dorso-ventral axis during planarian regeneration [54,55]. BMPs have also been associated with normal gastrointestinal tract development in vertebrates, and alterations in their expression caused morphological abnormalities in the intestinal tissues ([56,57]). BMP have also been shown to be required for proper cloaca formation [58]. Therefore, the overexpression of BMP during intestinal regeneration observed in our system is consistent with its role in the proper formation of the digestive tract organs.

• How do microarray results correlate with cellular processes that occur during intestinal regeneration?

In order to understand how changes in gene expression are associated with changes at the tissue and organ levels it is important to have a clear understanding of the cellular events that underlie the regenerative process. For this reason, our laboratory has described the intestinal regeneration process in H. glaberrima from the cellular perspective [6-12]. Two events in particular are to be highlighted in view of their possible association with the differential gene expression profile, one is the change in cytoskeletal proteins and the second the remodeling of the extracellular matrix.

Cytoskeletal Proteins

Previous studies from our laboratory have shown dramatic changes in the cytoskeleton of cells in the regenerating intestine. Many of these changes were associated with de-differentiation of the mesenterial muscle during the initial regeneration stages and with the processes of myogenesis, innervation and luminal epithelial formation in later regenerative stages [6,8,10,11]. In fact, during the regenerative process muscle cells actively extrude their contractile apparatus as spindle-like structures (SLSs) [8]. This has also been shown in other echinoderm regenerative phenomena [1,59]. We can speculate that the under-expression of some actin and myosin isoforms is associated with this dedifferentiation process. Similarly, the over expression of some actin isoforms could be associated with the migration of cells into the regenerating tissue [7,9], or with the eventual differentiation of cells from committed precursors [11]. Nonetheless, these hypotheses will need to be verified in future experiments using in situ hybridization and immunochemical analyses.

One of the most interesting findings is the under expression of gelsolin in early regenerative stages. Results by Cowin and colleagues [60] have shown that gelsolin mRNA was upregulated and gelsolin protein was associated with actin filaments in the fetal skin of embryonic day 19 rats, but not at embryonic day 17. This is the age gap when fetal wounds properties change from being scar-free to scar-forming. In culture, skin from 17 day embryonic rats can epithelialize an excisional wound in culture, but this ability is lost by embryonic day 19. Thus, our finding of a correlation between gelsolin gene under-expression and early intestinal regeneration suggests that gelsolin might play an important inhibitory role in the process of wound healing and regeneration.

ECM

One of the events that is common to most regenerative processes is the remodeling of the extracellular matrix. These changes, which have been well studied particularly during epithelial wound healing, involve the replacement of the ECM by a transient ECM that facilitates the cellular events associated with the regenerative response [61]. In the holothurian intestinal system we have previously shown a dramatic extracellular remodeling during the early stages of regeneration [12]. During this time collagen and laminin immunoreactivities disappear from the mesentery that gives rise to the new organ. Similarly there is a rise in the activity of MMPs that coincides with the major changes in the ECM. Moreover, we showed that inhibition of MMPs at this time point causes an inhibition of the regeneration process [12]. At the microscopic level, the changes of ECM and MMP activity coincide with an increase presence of phagocytic amoebocytes in the regenerating areas [8]. Therefore, the differential expression of various molecules associated with ECM remodeling, particularly the MMPs provide specific targets to define their role in the regeneration process.

• How do our results serve to advance the regeneration field?

The cellular and molecular processes of organ regeneration have received increasing attention in recent years [5]. In fact, the question of what controls organ regeneration? has been highlighted as one of the top scientific issues that deserves to be studied with high priority [62]. However, organ regeneration remains one of the least understood biological processes, particularly at the molecular level. Progress toward characterizing the molecular basis of regenerative processes has been modest for at least three reasons [63,64]. First, the regenerative capacity of most vertebrate animals is relatively limited [3], so popular biomedical models are not very useful. For example, some mammalian tissues (e.g., intestinal epithelia) and organs (e.g., skin and liver) can readily regenerate, but others, like heart and nervous tissues, have little or no regenerative capacity. Second, although development is a process closely akin to regeneration, most animals used extensively to study the molecular aspects of development (worms, fruitflies, zebrafish, chickens and mice) lack robust regenerative capacities. Third, those animals with remarkable regenerative capacities, such as coelenterates, flatworms, tunicates, newts and echinoderms have been little studied at the molecular/genetic level [1,20,65].

This has started to change and ongoing projects promise a bright future for the identification of regeneration-associated genes. Among the studies/models that need to be highlighted in terms of EST/genomic are those in amphibians [30,66,67], zebrafish [27,68,69], ascidians [70,71] planaria [4,55,72-77] and Hydra [78-81]. However, the only echinoderm studies are those of the sea urchin genome [82], which ironically is the echinoderm group with the least regenerative capacities [83].

Analyses of regeneration-associated genes in multiple models will lead to important comparative studies among species with highly developed regenerative capacities or between closely related species where good and poor regenerators are found. The importance of comparative studies using a large range of animal models was highlighted by Sanchez Alvarado and Tsonis [4] who described the differences and similarities in the regeneration mechanisms among diverse animal species. Other investigators pointed out that the echinoderms were largely excluded from this analysis [84], a fact explained by one of the authors [21] as being due to the lack of modern molecular tools to study the regeneration phenomena. Nonetheless, all investigators agreed that the application of genome-wide profiling approaches (as well as other tools) to the echinoderms would provide significant contributions to the understanding of regeneration. Thus, the present results represent an important step toward the molecular study of the amazing regenerative capacities of echinoderms and more specifically toward their use to better understand how organs can be regenerated.

Conclusion

Despite the apparent simplicity of the echinoderms, here we show that the organ regeneration is a complex process that involves the up-regulation and down-regulation of thousands of genes. Consistent with previous findings, these genes have been reported as involved in regeneration, development, ECM rearrangement, cytoskeleton reorganization and wound healing such as Wnt-9, BMP-1, Hox12, Tubulin, Tens-R, MMPs and gelsolin respectively. Moreover, a large number of unknown genes were also found to be differentially expressed at different stages of regeneration and might represent novel genes. Finally, here we show that intestine regeneration in the sea cucumber is a novel and important model for studies to identify and characterize the molecular basis of regenerative processes.

Methods

Microchip preparation

We have used the Agilent platform to design, perform and analyze custom-made H. glaberrima arrays. For this we used the "eArray" server from Agilent to design 60-mer probes that were synthesized on the microarray. A design with 15744 spots on each array and 8 arrays per slide was selected. All sequences from the cDNA libraries of normal and regenerating intestines were represented at least twice. Thus, 14352 (91.2%) probes were from 6287 clean sequences from H. glaberrima. We also spotted other sequences. These sequences represented genes associated to development, regeneration, proliferation, wound healing and growth. This pool included 150 probes from the sea cucumber Apostichopus japonicus (reported by Zheng et al. 2006), 8 from other holothurians that could be found in the NCBI database and 658 from the sea urchin S. purputatus (whose genome was recently sequenced [82]). Finally, 40 probes were from vertebrate organisms including mice, zebrafish and axolotl. Technical controls were 536 (3.4%).

Animals and treatments

Adult sea cucumbers (10~12 cm long) were collected from the north-eastern rocky shores of Puerto Rico. Animals were eviscerated at day 0, by intra-coelomic injection of 0.35 M KCl. They were left undisturbed in the aquaria to allow regeneration for 3, 7 and 14 days. At each time-point, animals were anesthetized by immersion in ice-cold sea water for 45 min. A dorsal incision was made to completely expose the internal cavity and allow the dissection of the intestinal primordium (3-, 7- or 14-dpe). After dissection, the regenerating intestines were placed in RNAlater® (Applied biosystems/Ambion, Austin TX) solution for RNA extraction. Uneviscerated (normal) animals were used as controls. One group was injected once with 1 mg lipopolysaccharides (LPS) diluted in 0.5 mL filtered sea water as reported by [85]. Intestines of the LPS treated animals were dissected 48 hrs after injection.

RNA Extraction

RNA was extracted using a combination of the [86] method using Tri-reagent® (N.93289, Sigma, St Louis, MO) and the RNAeasy mini kit from Quiagen (Valencia, CA) The concentration and integrity of RNA were determined using the NanoDrop-1000 Spectrophotometer (NanoDrop Technologies, Rockland, DE) and the 2100 Bioanalyzer (Agilent) with an RNA 6000 Nano LabChip® Kit at the Functional Genomics Research center (FGRc -UPR, Puerto Rico).

Synthesis and labeling of cRNA

RNA samples were amplified and labeled using the Low Input Fluorescent Linear Amplification kit (Agilent). Briefly, cDNA was first synthesized combining 300 ng of RNA with T7 promoter primers, 5× First Strand Buffer, 0.1 M DTT, 10 mM dNTP mix, MMLV RT and RNaseOUT. This reaction mix was incubated at 40°C for 2 h. Then, cRNA synthesis and dye incorporation were performed as follows: cDNA of the first reaction was combined with cyanine-3-CTP (Cy3, 10 mM) or cyanine 5-CTP (Cy5,10 mM) and the Transcription Master Mix (4× transcription buffer, 0.1 M DTT, NTP mix, 50% PEG, RNAseOUT, inorganic pyrophosphase and T7 RNA polymerase) at 40°C for 2 h. The fluorescently labeled cRNA products were purified using RNeasy mini spin columns (QIAGEN), subsequently analyzed for yield and dye incorporation using the NanoDrop, and finally stored at -80°C until needed. Only samples that presented a concentration of more than 8 pmol/mg (measurement representing efficiency of dye incorporation) were used for hybridization as recommended by the manufacturer.

Hybridization & Scanning

Cy3 and Cy5-labeled cRNAs were combined, fragmented, and hybridized to an 8 × 15K 60-mer oligo microarray (Agilent) for 17.5 h at 60°C with continuing rotation at 4 rpm. After hybridization, slide arrays were subjected to two successive washes (wash solution 1: 6× SSPE, 0.005% N-Lauroylsarcosine, and wash solution 2: 0.06× SSPE, 0.005% N-Lauroylsarcosine) and dried with the Stabilization and Drying Solution (Agilent). Arrays were immediately scanned in a G250B Microarray Scanner (Agilent) to obtain fluorescence intensities and Cy5/Cy3 ratios for each gene on the array (see below). Scanning parameters were used as recommended by the manufacturer. The resolution settings for scanning were 5 μm using minimum (10%) and maximum (100%) Photo Multiplier Tube detection sensitivities.

Microarray Experimental Design

The experimental design consisted of comparisons between Normal (non-eviscerated) intestinal tissues and 3-, 7-, and 14-dpe regenerating intestines. To determine differences in gene expression between regenerating stages we also compared directly 3- vs. 7-dpe and 7- vs. 14-dpe (Figure 8). Comparisons were also made between the intestinal tissues of normal and immune activated animals (injected LPS). There were only 48 genes differentially expressed between normal and LPS-immune activated intestines and these genes will be discussed elsewhere. (F. Ramirez-Gomez et al. in prep.) These genes were subtracted from the regeneration profile to ensure that genes associated with the regeneration phenomenon were not differentially expressed due any possible immune activation (e.g., because of bacterial invasion during the evisceration process).

Technically, for each comparison, the microarray experiment was done twice, one using different biological samples, each consisting of a sample pool from 3 animals. Of the two biological replicates one was always performed as a dye-swap in order to correct any dye bias.

thumbnailFigure 8. Microarray experimental design. Normal (Nml) and regenerating intestines at 3, 7 and 14 dpe (days post-evisceration) were compared among them. Each experiment was replicated as a dye swap with different biological samples. LPS treated tissues where used to normalize the expression for immune activated genes. Arrows show direct comparisons.

Microarray data analysis

Microarray data was extracted with Agilent's Feature Extraction software (version 5.1) and analyzed using the Limma package of Bioconductor. The data discussed in this publication have been deposited in NCBI's Gene Expression Omnibus [87] and are accessible through GEO Series accession number GSE16182 http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE16182 webcite. The analysis of the microarray data consisted of the following steps: 1) within-array and between-array normalizations; 2) fitting the data to a linear model; and 3) computing differential gene expression. For normalization purposes [88] MA-plots were generated representing the (R, G) data (R = red for Cy5 and G = green for Cy3), in which the log ratio of R versus G (M value = log2 R/G) was plotted against the overall intensity of each spot (A value = log2 (R + G)/2. Within-array normalization was first applied and M-values were normalized within each array using the Global Loess Normalization method. Aquantile normalization was then applied to the A-values as a method for between-array normalization, to assure that the intensities and log-ratios had similar distributions across arrays. To estimate the average M-value for each gene and assess differential gene expression, a simple linear model was fit to the data, and M-value averages and standard deviations for each gene were obtained. To find genes with significant expression changes between groups, empirical Bayes statistics were applied to the data by moderating the standard errors of the estimated M-values. P-values were obtained from the moderated t-statistic and corrected for multiple testing with the [89] method. The null hypothesis, that there is no differential expression of genes between regeneration stages compared with normal tissues, was rejected for p-values lower than 0.05. Thus, the change in expression is given by the fold change while the believability of the change is given by the odds. Alternative to this method, the GenePix Pro 6.0.1.26 and Acuity 4.0.0.60 software (from Molecular Devices, Sunnyvale CA) was also used to analyze gene expression values (log ratios) and perform clustering analyses. Methods used to visualize the expression profiles included hierarchical [90] and non-hierarchical clustering such as Self-Organizing Maps (SOM, [91]).

Gene expression analysis

The differential expression of the genes of interest was semi-quantitatively determined by gene-specific relative RT-PCR. Primers were designed for optimal performance using the Primer3 and Net primer web-servers [see additional file 4 for the list of primer sequences]. All primers were synthesized by AlphaDNA (Montreal, Canada). PCRs were performed in the Mastercycler ep (Eppendorf, Westbury, NY) using the Promega's Taq-Polimerase kit (Madison, USA). Cycling conditions for the amplified products were as follow: 94°C × 1 min, 55°C (± 3°C) × 1 min, 72°C × 1 min; all performed for 35~40 cycles. Optical densities of the PCR products were measured using Photoshop (Adobe). Values were extracted for background and then normalized against NADH-dehydrogenase unit 5 (NADH, 241 bp). Each stage of regeneration was analyzed with at least an N = 7, each sample represented the intestine of one individual. Statistical test of mean comparison was carried in R http://cran.r-project.org webcite under the null hypothesis that there is no differential expression between different stages of regeneration and normal tissues. Significance of the paired t-test was assumed at p < 0.05.

Additional file 4. Primers. List of primers used for validation.

Format: PDF Size: 27KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

List of abbreviations

cDNA: complementary DNA; EST: Expressed sequence tag; Dpe: Days post evisceration; RT-PCR: reverse transcriptase polymerase chain reaction; LPS: lipopolysaccharides; MMP: Matrix metalloproteases; GO: Gene Ontology; ORF: Open Reading Frame; UTR: Untranslated region; mRNA: Messenger RNA; SLS: spindle like structure; ECM: Extra-Cellular Matrix.

Authors' contributions

POP designed and carried out the study, analyzed and interpreted the data, and drafted the manuscript. FRG helped in the design and participated in the molecular experiments. CRC, JRL, ESC, ATMM and KT participated in cDNA library construction, EST characterization and analysis and gene identification and sequencing. JPO, SGD, FSD, JH, and CDA performed most of the validation experiments and participated in the analysis of data. HOZ helped in bioinformatics analysis. JEGA conceived the study, participated in its design, data analyses and helped in the manuscript preparation. All authors have given final approval of the version to be published.

Acknowledgements

This work was supported by NIH (1SC1GM084770-01), NSF (IBN-0110692) and NIH- (S06GM08102). We also acknowledge partial support from NIH INBRE (P20RR16470), DEGI-PBDT (doctoral dissertation fellowship to F. Ramírez-Gómez and E. Suarez-Castillo), the Sequencing and Genotyping facility of the UPR and the University of Puerto Rico.

References

  1. Candia Carnevali MD, Bonasoro F: Introduction to the biology of regeneration in echinoderms.

    Microsc Res Tech 2001, 55(6):365-368. PubMed Abstract | Publisher Full Text OpenURL

  2. Brockes JP: Amphibian limb regeneration: rebuilding a complex structure.

    Science 1997, 276(5309):81-87. PubMed Abstract | Publisher Full Text OpenURL

  3. Tsonis PA: Regeneration in vertebrates.

    Dev Biol 2000, 221(2):273-284. PubMed Abstract | Publisher Full Text OpenURL

  4. Sanchez Alvarado A, Tsoni PA: Bridging the regeneration gap: genetic insights from diverse animal models.

    Nat Rev Genet 2006, 7(11):873-884. PubMed Abstract | Publisher Full Text OpenURL

  5. Slack JM: Regeneration research today.

    Dev Dyn 2003, 226(2):162-166. PubMed Abstract | Publisher Full Text OpenURL

  6. Garcia-Arraras JE, Estrada-Rodgers L, Santiago R, Torres II, Diaz-Miranda L, Torres-Avillan I: Cellular mechanisms of intestine regeneration in the sea cucumber, Holothuria glaberrima Selenka (Holothuroidea:Echinodermata).

    J Exp Zool 1998, 281(4):288-304. PubMed Abstract | Publisher Full Text OpenURL

  7. Cabrera-Serrano A, Garcia-Arraras JE: RGD-containing peptides inhibit intestinal regeneration in the sea cucumber Holothuria glaberrima.

    Dev Dyn 2004, 231(1):171-178. PubMed Abstract | Publisher Full Text OpenURL

  8. Candelaria AG, Murray G, File SK, Garcia-Arraras JE: Contribution of mesenterial muscle dedifferentiation to intestine regeneration in the sea cucumber Holothuria glaberrima.

    Cell Tissue Res 2006, 325(1):55-65. PubMed Abstract | Publisher Full Text OpenURL

  9. Garcia-Arraras JE, Schenk C, Rodrigues-Ramirez R, Torres II, Valentin G, Candelaria AG: Spherulocytes in the echinoderm Holothuria glaberrima and their involvement in intestinal regeneration.

    Dev Dyn 2006, 235(12):3259-3267. PubMed Abstract | Publisher Full Text OpenURL

  10. Garcia-Arraras JE, Torres-Avillan I: Developmental expression of galanin-like immunoreactivity by members of the avian sympathoadrenal cell lineage.

    Cell Tissue Res 1999, 295(1):33-41. PubMed Abstract | Publisher Full Text OpenURL

  11. Murray G, Garcia-Arraras JE: Myogenesis during holothurian intestinal regeneration.

    Cell Tissue Res 2004, 318(3):515-524. PubMed Abstract | Publisher Full Text OpenURL

  12. Quinones JL, Rosa R, Ruiz DL, Garcia-Arraras JE: Extracellular matrix remodeling and metalloproteinase involvement during intestine regeneration in the sea cucumber Holothuria glaberrima.

    Dev Biol 2002, 250(1):181-197. PubMed Abstract | Publisher Full Text OpenURL

  13. Rojas-Cartagena C, Ortiz-Pineda P, Ramirez-Gomez F, Suarez-Castillo EC, Matos-Cruz V, Rodriguez C, Ortiz-Zuazaga H, Garcia-Arraras JE: Distinct profiles of expressed sequence tags during intestinal regeneration in the sea cucumber Holothuria glaberrima.

    Physiol Genomics 2007, 31(2):203-215. PubMed Abstract | Publisher Full Text OpenURL

  14. Santiago P, Roig-Lopez JL, Santiago C, Garcia-Arraras JE: Serum amyloid A protein in an echinoderm: its primary structure and expression during intestinal regeneration in the sea cucumber Holothuria glaberrima.

    J Exp Zool 2000, 288(4):335-344. PubMed Abstract | Publisher Full Text OpenURL

  15. Suarez-Castillo EC, Medina-Ortiz WE, Roig-Lopez JL, Garcia-Arraras JE: Ependymin, a gene involved in regeneration and neuroplasticity in vertebrates, is overexpressed during regeneration in the echinoderm Holothuria glaberrima.

    Gene 2004, 334:133-143. PubMed Abstract | Publisher Full Text OpenURL

  16. Bannister R, McGonnell IM, Graham A, Thorndyke MC, Beesley PW: Afuni, a novel transforming growth factor-beta gene is involved in arm regeneration by the brittle star Amphiura filiformis.

    Dev Genes Evol 2005, 215(8):393-401. PubMed Abstract | Publisher Full Text OpenURL

  17. Bannister SR, Powell CA: Foreign body reaction to anorganic bovine bone and autogenous bone with platelet-rich plasma in guided bone regeneration.

    J Periodontol 2008, 79(6):1116-1120. PubMed Abstract | Publisher Full Text OpenURL

  18. Patruno M, McGonnell I, Graham A, Beesley P, Candia Carnevali MD, Thorndyke M: Anbmp2/4 is a new member of the transforming growth factor-beta superfamily isolated from a crinoid and involved in regeneration.

    Proc Biol Sci 2003, 270(1522):1341-1347. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  19. Patruno M, Smertenko A, Candia Carnevali MD, Bonasoro F, Beesley PW, Thorndyke MC: Expression of transforming growth factor beta-like molecules in normal and regenerating arms of the crinoid Antedon mediterranea: immunocytochemical and biochemical evidence.

    Proc Biol Sci 2002, 269(1502):1741-1747. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  20. Thorndyke MC, Chen WC, Beesley PW, Patruno M: Molecular approach to echinoderm regeneration.

    Microsc Res Tech 2001, 55(6):474-485. PubMed Abstract | Publisher Full Text OpenURL

  21. Sanchez Alvarado A: Stem cells and the Planarian Schmidtea mediterranea.

    C R Biol 2007, 330(67):498-503. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. Carlson ME, O'Connor MS, Hsu M, Conboy IM: Notch signaling pathway and tissue engineering.

    Front Biosci 2007, 12:5143-5156. PubMed Abstract | Publisher Full Text OpenURL

  23. Stocum DL: Stem cells in CNS and cardiac regeneration.

    Adv Biochem Eng Biotechnol 2005, 93:135-159. PubMed Abstract OpenURL

  24. Mendez AT, Roig-Lopez JL, Santiago P, Santiago C, Garcia-Arraras JE: Identification of Hox Gene Sequences in the Sea Cucumber Holothuria glaberrima Selenka (Holothuroidea: Echinodermata).

    Mar Biotechnol (NY) 2000, 2(3):231-240. PubMed Abstract OpenURL

  25. Mendez AT: Characterization of homeobox genes in the sea cucumber, Holothuria glaberrima Selenka (holothuroidea:Echinodermata) and expression during the regeneration of gastrointestinal tract. San Juan: University of Puerto Rico; 2001.

  26. Santiago-Cardona PG, Berrios CA, Ramirez F, Garcia-Arraras JE: Lipopolysaccharides induce intestinal serum amyloid A expression in the sea cucumber Holothuria glaberrima.

    Dev Comp Immunol 2003, 27(2):105-110. PubMed Abstract | Publisher Full Text OpenURL

  27. Lien CL, Schebesta M, Makino S, Weber GJ, Keating MT: Gene expression analysis of zebrafish heart regeneration.

    PLoS Biol 2006, 4(8):e260. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  28. Arai M, Yokosuka O, Chiba T, Imazeki F, Kato M, Hashida J, Ueda Y, Sugano S, Hashimoto K, Saisho H, et al.: Gene expression profiling reveals the mechanism and pathophysiology of mouse liver regeneration.

    J Biol Chem 2003, 278(32):29813-29818. PubMed Abstract | Publisher Full Text OpenURL

  29. Fukuhara Y, Hirasawa A, Li XK, Kawasaki M, Fujino M, Funeshima N, Katsuma S, Shiojima S, Yamada M, Okuyama T, et al.: Gene expression profile in the regenerating rat liver after partial hepatectomy.

    J Hepatol 2003, 38(6):784-792. PubMed Abstract | Publisher Full Text OpenURL

  30. 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(1):27-40. PubMed Abstract | Publisher Full Text OpenURL

  31. Beauchemin M, Noiseux N, Tremblay M, Savard P: Expression of Hox A11 in the limb and the regeneration blastema of adult newt.

    Int J Dev Biol 1994, 38(4):641-649. PubMed Abstract | Publisher Full Text OpenURL

  32. Carlson ER: Bone grafting the jaws in the 21st century: the use of platelet-rich plasma and bone morphogenetic protein.

    Alpha Omegan 2000, 93(3):26-30. PubMed Abstract OpenURL

  33. Savard P, Gates PB, Brockes JP: Position dependent expression of a homeobox gene transcript in relation to amphibian limb regeneration.

    Embo J 1988, 7(13):4275-4282. PubMed Abstract | PubMed Central Full Text OpenURL

  34. Schummer M, Scheurlen I, Schaller C, Galliot B: HOM/HOX homeobox genes are present in hydra (Chlorohydra viridissima) and are differentially expressed during regeneration.

    Embo J 1992, 11(5):1815-1823. PubMed Abstract | PubMed Central Full Text OpenURL

  35. Mizuta I, Ogasawara N, Yoshikawa H, Sakoyama Y: Identification of homeobox genes expressed during the process of rat liver regeneration after partial hepatectomy.

    Biochem Genet 1996, 34(1–2):1-15. PubMed Abstract | Publisher Full Text OpenURL

  36. Thummel R, Ju M, Sarras MP Jr, Godwin AR: Both Hoxc13 orthologs are functionally important for zebrafish tail fin regeneration.

    Dev Genes Evol 2007, 217(6):413-420. PubMed Abstract | Publisher Full Text OpenURL

  37. Brown R, Brockes JP: Identification and expression of a regeneration-specific homeobox gene in the newt limb blastema.

    Development 1991, 111(2):489-496. PubMed Abstract | Publisher Full Text OpenURL

  38. Beck F, Tata F, Chawengsaksophak K: Homeobox genes and gut development.

    Bioessays 2000, 22(5):431-441. PubMed Abstract | Publisher Full Text OpenURL

  39. Kawazoe Y, Sekimoto T, Araki M, Takagi K, Araki K, Yamamura K: Region-specific gastrointestinal Hox code during murine embryonal gut development.

    Dev Growth Differ 2002, 44(1):77-84. PubMed Abstract | Publisher Full Text OpenURL

  40. Wolgemuth DJ, Behringer RR, Mostoller MP, Brinster RL, Palmiter RD: Transgenic mice overexpressing the mouse homoeobox-containing gene Hox1.4 – exhibit abnormal gut development.

    Nature 1989, 337(6206):464-467. PubMed Abstract | Publisher Full Text OpenURL

  41. Kondo T, Dolle P, Zakany J, Duboule D: Function of posterior HoxD genes in the morphogenesis of the anal sphincter.

    Development 1996, 122(9):2651-2659. PubMed Abstract | Publisher Full Text OpenURL

  42. Yokoyama H, Ogino H, Stoick-Cooper CL, Grainger RM, Moon RT: Wnt/beta-catenin signaling has an essential role in the initiation of limb regeneration.

    Dev Biol 2007, 306(1):170-178. PubMed Abstract | Publisher Full Text OpenURL

  43. Hayashi T, Mizuno N, Takada R, Takada S, Kondoh H: Determinative role of Wnt signals in dorsal iris-derived lens regeneration in newt eye.

    Mech Dev 2006, 123(11):793-800. PubMed Abstract | Publisher Full Text OpenURL

  44. Kim JB, Leucht P, Lam K, Luppen C, Ten Berge D, Nusse R, Helms JA: Bone regeneration is regulated by wnt signaling.

    J Bone Miner Res 2007, 22(12):1913-1923. PubMed Abstract | Publisher Full Text OpenURL

  45. Ito M, Yang Z, Andl T, Cui C, Kim N, Millar SE, Cotsarelis G: Wnt-dependent de novo hair follicle regeneration in adult mouse skin after wounding.

    Nature 2007, 447(7142):316-320. PubMed Abstract | Publisher Full Text OpenURL

  46. Mount JG, Muzylak M, Allen S, Althnaian T, McGonnell IM, Price JS: Evidence that the canonical Wnt signalling pathway regulates deer antler regeneration.

    Dev Dyn 2006, 235(5):1390-1399. PubMed Abstract | Publisher Full Text OpenURL

  47. Tan X, Behari J, Cieply B, Michalopoulos GK, Monga SP: Conditional deletion of beta-catenin reveals its role in liver growth and regeneration.

    Gastroenterology 2006, 131(5):1561-1572. PubMed Abstract | Publisher Full Text OpenURL

  48. Yen TH, Wright NA: The gastrointestinal tract stem cell niche.

    Stem Cell Rev 2006, 2(3):203-212. PubMed Abstract | Publisher Full Text OpenURL

  49. Kobayashi C, Saito Y, Ogawa K, Agata K: Wnt signaling is required for antero-posterior patterning of the planarian brain.

    Dev Biol 2007, 306(2):714-724. PubMed Abstract | Publisher Full Text OpenURL

  50. Gurley KA, Rink JC, Sanchez Alvarado A: Beta-catenin defines head versus tail identity during planarian regeneration and homeostasis.

    Science 2008, 319(5861):323-327. PubMed Abstract | Publisher Full Text OpenURL

  51. Grogg MW, Call MK, Okamoto M, Vergara MN, Del Rio-Tsonis K, Tsonis PA: BMP inhibition-driven regulation of six-3 underlies induction of newt lens regeneration.

    Nature 2005, 438(7069):858-862. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  52. Smith A, Avaron F, Guay D, Padhi BK, Akimenko MA: Inhibition of BMP signaling during zebrafish fin regeneration disrupts fin growth and scleroblasts differentiation and function.

    Dev Biol 2006, 299(2):438-454. PubMed Abstract | Publisher Full Text OpenURL

  53. Beck CW, Christen B, Barker D, Slack JM: Temporal requirement for bone morphogenetic proteins in regeneration of the tail and limb of Xenopus tadpoles.

    Mech Dev 2006, 123(9):674-688. PubMed Abstract | Publisher Full Text OpenURL

  54. Molina MD, Salo E, Cebria F: The BMP pathway is essential for re-specification and maintenance of the dorsoventral axis in regenerating and intact planarians.

    Dev Biol 2007, 311(1):79-94. PubMed Abstract | Publisher Full Text OpenURL

  55. Reddien PW, Bermange AL, Kicza AM, Sanchez Alvarado A: BMP signaling regulates the dorsal planarian midline and is needed for asymmetric regeneration.

    Development 2007, 134(22):4043-4051. PubMed Abstract | Publisher Full Text OpenURL

  56. Batts LE, Polk DB, Dubois RN, Kulessa H: Bmp signaling is required for intestinal growth and morphogenesis.

    Dev Dyn 2006, 235(6):1563-1570. PubMed Abstract | Publisher Full Text OpenURL

  57. De Santa Barbara P, Williams J, Goldstein AM, Doyle AM, Nielsen C, Winfield S, Faure S, Roberts DJ: Bone morphogenetic protein signaling pathway plays multiple roles during gastrointestinal tract development.

    Dev Dyn 2005, 234(2):312-322. PubMed Abstract | Publisher Full Text OpenURL

  58. Pyati UJ, Cooper MS, Davidson AJ, Nechiporuk A, Kimelman D: Sustained Bmp signaling is essential for cloaca development in zebrafish.

    Development 2006, 133(11):2275-2284. PubMed Abstract | Publisher Full Text OpenURL

  59. Dolmatov I: Regeneration of the aquapharyngeal complex in the holothurian Eupentacta fraudatrix (Holothuroidea, Dendrochirota).

    Monogr Dev Biol 1992, 23:40-50. PubMed Abstract OpenURL

  60. Cowin AJ, Hatzirodos N, Teusner JT, Belford DA: Differential effect of wounding on actin and its associated proteins, paxillin and gelsolin, in fetal skin explants.

    J Invest Dermatol 2003, 120(6):1118-1129. PubMed Abstract | Publisher Full Text OpenURL

  61. Martin P: Wound healing – aiming for perfect skin regeneration.

    Science 1997, 276(5309):75-81. PubMed Abstract | Publisher Full Text OpenURL

  62. Davenport RJ: What controls organ regeneration?

    Science 2005, 309(5731):84. PubMed Abstract | Publisher Full Text OpenURL

  63. Sanchez Alvarado A: Regeneration in the metazoans: why does it happen?

    Bioessays 2000, 22(6):578-590. PubMed Abstract | Publisher Full Text OpenURL

  64. Sanchez Alvarado A: Regeneration and the need for simpler model organisms.

    Philos Trans R Soc Lond B Biol Sci 2004, 359(1445):759-763. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  65. Brockes JP, Kumar A: Plasticity and reprogramming of differentiated cells in amphibian regeneration.

    Nat Rev Mol Cell Biol 2002, 3(8):566-574. PubMed Abstract | Publisher Full Text OpenURL

  66. Habermann B, Bebin AG, Herklotz S, Volkmer M, Eckelt K, Pehlke K, Epperlein HH, Schackert HK, Wiebe G, Tanaka EM: An Ambystoma mexicanum EST sequencing project: analysis of 17,352 expressed sequence tags from embryonic and regenerating blastema cDNA libraries.

    Genome Biol 2004, 5(9):R67. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  67. Putta S, Smith JJ, Walker JA, Rondet M, Weisrock DW, Monaghan J, Samuels AK, Kump K, King DC, Maness NJ, et al.: From biomedicine to natural history research: EST resources for ambystomatid salamanders.

    BMC Genomics 2004, 5(1):54. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  68. Andreasen EA, Mathew LK, Lohr CV, Hasson R, Tanguay RL: Aryl hydrocarbon receptor activation impairs extracellular matrix remodeling during zebra fish fin regeneration.

    Toxicol Sci 2007, 95(1):215-226. PubMed Abstract | Publisher Full Text OpenURL

  69. Cameron DA, Gentile KL, Middleton FA, Yurco P: Gene expression profiles of intact and regenerating zebrafish retina.

    Mol Vis 2005, 11:775-791. PubMed Abstract | Publisher Full Text OpenURL

  70. Azumi K, Takahashi H, Miki Y, Fujie M, Usami T, Ishikawa H, Kitayama A, Satou Y, Ueno N, Satoh N: Construction of a cDNA microarray derived from the ascidian Ciona intestinalis.

    Zoolog Sci 2003, 20(10):1223-1229. PubMed Abstract | Publisher Full Text OpenURL

  71. Rinkevi Y, Douekch J, Haber O, Rinkevich B, Reshef R: Urochordate whole body regeneration inaugurates a diverse innate immune signaling profile.

    Dev Biol 2007, 312(1):131-146. PubMed Abstract | Publisher Full Text OpenURL

  72. Mineta K, Nakazawa M, Cebria F, Ikeo K, Agata K, Gojobori T: Origin and evolutionary process of the CNS elucidated by comparative genomics analysis of planarian ESTs.

    Proc Natl Acad Sci USA 2003, 100(13):7666-7671. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  73. Nakazawa M, Cebria F, Mineta K, Ikeo K, Agata K, Gojobori T: Search for the evolutionary origin of a brain: planarian brain characterized by microarray.

    Mol Biol Evol 2003, 20(5):784-791. PubMed Abstract | Publisher Full Text OpenURL

  74. Newmark PA, Sanchez Alvarado A: Not your father's planarian: a classic model enters the era of functional genomics.

    Nat Rev Genet 2002, 3(3):210-219. PubMed Abstract | Publisher Full Text OpenURL

  75. Rossi L, Salvetti A, Marincola FM, Lena A, Deri P, Mannini L, Batistoni R, Wang E, Gremigni V: Deciphering the molecular machinery of stem cells: a look at the neoblast gene expression profile.

    Genome Biol 2007, 8(4):R62. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  76. Zayas RM, Hernandez A, Habermann B, Wang Y, Stary JM, Newmark PA: The planarian Schmidtea mediterranea as a model for epigenetic germ cell specification: analysis of ESTs from the hermaphroditic strain.

    Proc Natl Acad Sci USA 2005, 102(51):18491-18496. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  77. Cebria F, Kudome T, Nakazawa M, Mineta K, Ikeo K, Gojobori T, Agata K: The expression of neural-specific genes reveals the structural and molecular complexity of the planarian central nervous system.

    Mech Dev 2002, 116(1–2):199-204. PubMed Abstract | Publisher Full Text OpenURL

  78. Chera S, Kaloulis K, Galliot B: The cAMP response element binding protein (CREB) as an integrative HUB selector in metazoans: clues from the hydra model system.

    Biosystems 2007, 87(2–3):191-203. PubMed Abstract | Publisher Full Text OpenURL

  79. Galliot B, Miljkovic-Licina M, de Rosa R, Chera S: Hydra, a niche for cell and developmental plasticity.

    Semin Cell Dev Biol 2006, 17(4):492-502. PubMed Abstract | Publisher Full Text OpenURL

  80. Galliot B, Miljkovic-Licina M, Ghila L, Chera S: RNAi gene silencing affects cell and developmental plasticity in hydra.

    C R Biol 2007, 330(6–7):491-497. PubMed Abstract | Publisher Full Text OpenURL

  81. Stout T, McFarland T, Appukuttan B: Suppression subtractive hybridization identifies novel transcripts in regenerating Hydra littoralis.

    J Biochem Mol Biol 2007, 40(2):286-289. PubMed Abstract | Publisher Full Text OpenURL

  82. Sea-Urchin-Genome-Sequence-Consortium: The genome of the sea urchin Strongylocentrotus purpuratus.

    Science 2006, 314(5801):941-952. PubMed Abstract | Publisher Full Text OpenURL

  83. Dubois P, Ameye L: Regeneration of spines and pedicellariae in echinoderms: a review.

    Microsc Res Tech 2001, 55(6):427-437. PubMed Abstract | Publisher Full Text OpenURL

  84. Dupont S, Thorndyke MC: Growth or differentiation? Adaptive regeneration in the brittlestar Amphiura filiformis.

    J Exp Biol 2006, 209(Pt 19):3873-3881. PubMed Abstract | Publisher Full Text OpenURL

  85. Ramirez-Gomez F, Ortiz-Pineda PA, Rojas-Cartagena C, Suarez-Castillo EC, Garcia-Arraras JE: Immune-related genes associated with intestinal tissue in the sea cucumber Holothuria glaberrima.

    Immunogenetics 2008, 60(1):57-71. PubMed Abstract | Publisher Full Text OpenURL

  86. Chomczynski P: A reagent for the single-step simultaneous isolation of RNA, DNA and proteins from cell and tissue samples.

    Biotechniques 1993, 15(3):532-534.

    536-537.

    PubMed Abstract OpenURL

  87. Edgar R, Domrachev M, Lash AE: Gene Expression Omnibus: NCBI gene expression and hybridization array data repository.

    Nucleic Acids Res 2002, 30(1):207-210. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  88. Smyth GK, Speed T: Normalization of cDNA microarray data.

    Methods 2003, 31(4):265-273. PubMed Abstract | Publisher Full Text OpenURL

  89. Benjamini Y, Hochber Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing.

    J R Stat Soc 1995, 57:289-300. OpenURL

  90. Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns.

    Proc Natl Acad Sci USA 1998, 95(25):14863-14868. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  91. Kohonen T: Cortical maps.

    Nature 1990, 346(6279):24. PubMed Abstract | Publisher Full Text OpenURL