Production of a reference transcriptome and transcriptomic database (PocilloporaBase) for the cauliflower coral, Pocillopora damicornis
© Traylor-Knowles et al; licensee BioMed Central Ltd. 2011
Received: 16 February 2011
Accepted: 29 November 2011
Published: 29 November 2011
Motivated by the precarious state of the world's coral reefs, there is currently a keen interest in coral transcriptomics. By identifying changes in coral gene expression that are triggered by particular environmental stressors, we can begin to characterize coral stress responses at the molecular level, which should lead to the development of more powerful diagnostic tools for evaluating the health of corals in the field. Furthermore, the identification of genetic variants that are more or less resilient in the face of particular stressors will help us to develop more reliable prognoses for particular coral populations. Toward this end, we performed deep mRNA sequencing of the cauliflower coral, Pocillopora damicornis, a geographically widespread Indo-Pacific species that exhibits a great diversity of colony forms and is able to thrive in habitats subject to a wide range of human impacts. Importantly, P. damicornis is particularly amenable to laboratory culture. We collected specimens from three geographically isolated Hawaiian populations subjected to qualitatively different levels of human impact. We isolated RNA from colony fragments ("nubbins") exposed to four environmental stressors (heat, desiccation, peroxide, and hypo-saline conditions) or control conditions. The RNA was pooled and sequenced using the 454 platform.
Both the raw reads (n = 1, 116, 551) and the assembled contigs (n = 70, 786; mean length = 836 nucleotides) were deposited in a new publicly available relational database called PocilloporaBase http://www.PocilloporaBase.org. Using BLASTX, 47.2% of the contigs were found to match a sequence in the NCBI database at an E-value threshold of ≤.001; 93.6% of those contigs with matches in the NCBI database appear to be of metazoan origin and 2.3% bacterial origin, while most of the remaining 4.1% match to other eukaryotes, including algae and amoebae.
P. damicornis now joins the handful of coral species for which extensive transcriptomic data are publicly available. Through PocilloporaBase http://www.PocilloporaBase.org, one can obtain assembled contigs and raw reads and query the data according to a wide assortment of attributes including taxonomic origin, PFAM motif, KEGG pathway, and GO annotation.
Over the long-term, populations can respond to stressful environmental conditions via adaptive evolution. However, over the short-term, organisms under stress must alter their physiology or behavior, and doing so commonly involves changes in gene expression. Therefore, if we are to understand how organisms respond to stress, we must understand how stress alters gene expression.
In recent years, a number of microarray studies have been undertaken to reconstruct the gene expression profiles of corals under stress [1–6]. However, the chip-based approaches used to date have not been able to evaluate the full scope of the stress response because only a fraction of potential transcripts have been represented on the chips. In the absence of a tiling array, whose production must await the sequencing and complete assembly of a coral genome, next-generation sequencing technologies are the only available approach for characterizing the full scope of the transcriptional response to stress. Currently, there is no complete coral genome available, and the most closely related genome that is publicly available is that of the starlet sea anemone Nematostella vectensis.
In the environment, corals are exposed to a variety of natural and anthropogenic stressors. Once the stressors reach a tipping point, dramatic physiological changes can occur very abruptly, including bleaching, which involves the expulsion of all symbiotic zooxanthellae [8–15]. While the loss of symbionts may be adaptive in the short-term, once bleaching occurs, if a coral colony continues to experience stressful conditions, it will rapidly succumb--scleractinian corals cannot persist indefinitely as heterotrophs. Therefore, a critical component of coral resilience is their ability to stave off bleaching through molecular stress-response mechanisms.
To date, the identification of "stress response genes" in corals and other cnidarians has been based largely on homology to functionally characterized genes in other animals (e.g., ). This approach has validity, as there is accumulating evidence that many cnidarian stress responses are likely to be largely homologous to those of triploblastic animals. However, given the long history of the Cnidaria as an independent evolutionary lineage and the many unique aspects of their biology, the cnidarian stress response repertoire is also certain to differ in key respects from that of other animals.
Published transcriptomic datasets for scleractinians
Source of RNA
Avg. Length (nt)
larvae; heat-stressed larvae; larvae treated with settlement inducer
developmental stages, incl. adults; with and without symbionts
developmental stages, incl. adults; with and without symbionts
adult colonies subject to a battery of stressors
Construction and content
PocilloporaBase Database Construction
PocilloporaBase is a relational database constructed in PostgreSQL (version 8.4.4). It houses the P. damicornis contigs generated in this study in addition to the results from a number of bioinformatics analyses performed on these contigs. The database structure and entity relationships are depicted in Additional File 1.
The database comprises thirteen tables, eight entity tables and five relations tables. The "Contigs" table houses output from the assembly including the nucleotide sequence of each contig and other key features of the contig, including its length in nucleotides, possible start sites, the total number of reads used to build the contig, and the average sequencing coverage of each nucleotide position within the contig. Through the "Hmm" table, each contig ID (CID) is linked to 0, 1 or more protein family IDs (PFID) based on a search of Pfam (an HMM search implemented in Perl was used to compare each of the P. damicornis contigs to the conserved protein domains housed at Pfam; ftp://ftp.sanger.ac.uk/pub/databases/Pfam/Tools/OldPfamScan/pfam_scan.pl). The "Pfam" table links each PFID to the name and description of the relevant protein domain. BLASTX output is summarized in the "Blast" table. Here, each CID is associated with up to five protein sequences in NCBI that produced significant BLASTX hits to the relevant P. damicornis contig. Various metrics from the BLAST hits are linked to each protein identifier (PID) including the BLAST alignment length, the bit score, the expect score, and the number of identical amino acids. Through an "Annotation" table, the PIDs are cross-referenced with gene ontology IDs (GOID). The "Ontology" table links each GOID to the name and description of the relevant gene ontology term. The "Species" table relates each PID to the species from which the corresponding protein sequence was derived and indicates whether that species is a known symbiont of P. damicornis. The "Accession" table cross-references each PID with the corresponding mRNA ID (MID) if this information is available in NCBI. The "mRNA" table links these MIDs to the corresponding mRNA name and nucleotide sequence. Through the "Tapp" table, each PID may be cross-referenced to 0, 1, or more biochemical pathways housed in the KEGG database. The "Kegg" table associates each Kegg ID (KID) with the name and definition of the corresponding biochemical pathway.
The database can be accessed via an HTML interface at http://www.PocilloporaBase.org. Users can query the database via any one of ten different search terms (described in detail in the Utility and Discussion section). Database queries are generated using the Python programming language, version 2.6.5, and the PyGreSQL DB-API module. Web pages summarizing the search results are produced as needed via Python scripts that generate HTML output. The database is hosted on a computer running Ubuntu Linux 10.04.2 using the Apache HTTP Server (version 2.2.14).
Adult coral RNA extraction, and sequencing
Assembly, Identification of contigs, and Pfam domain searches
After sequencing, short reads (< 40 nucleotides) and low quality reads that did not overlap with other sequencing reads were discarded, and the remaining 955, 910 sequencing reads were assembled at Beckman Coulter using MIRA3 . High quality singletons were included with the contigs in all subsequent analyses. Contig sequences were blasted against the adaptor sequences used for both the library preparation and the sequencing to insure that none of the contig sequences were contaminated with adaptors. Adaptor sequences were trimmed, and the assembled contigs were used to sequentially query the NCBI non-redundant protein database using BLASTX with an E-value cut-off of 0.001. The top five gene hits were assigned to each contig. All five of the top hits usually agreed on gene ontology and taxonomy, but where they disagreed, we associated multiple GO terms and multiple possible taxonomic affinities with a given contig. In order to identify conserved protein domains, all six open reading frames were blasted against the protein domain database at Pfam . Hits were retained only if they had an e-value cut off lower than 0.001.
GO term and KEGG analysis
Using the top five hits from the BLASTX search described above, each contig was assigned a list of associated protein GI numbers. GI numbers were converted to Entrez Gene IDs using the gene2accession conversion file from NCBI. Gene2go was then used to obtain relevant GO annotation for the five top BLASTX hits to each P. damicornis contig, and the GO term(s) were then associated with their respective contig. Protein GI numbers were cross-referenced to species-specific KEGG pathways . Using KEGG's lists of annotated plants and animals, these pathways were organized into corresponding lists, and generalized KEGG pathway Ids were obtained. KEGG pathway analysis was then performed, and individual contigs were then mapped to different biochemical pathways using IPath . The top 5 hits were chosen to increase the probability of finding a hit that would allow pairing each contig with its corresponding GO category and with the existent KEGG data. Using photosynthesis as an example, the top hits matched to Symbiodinium sequences, not plants. Because Symbiodinium is not represented within the KEGG database, the contigs did not appear to represent enzymes involved in photosynthesis, when, in fact, they did.
Utility and Discussion
Sequencing yielded 1, 116, 551 raw reads with an average length of 379 nucleotides (range: 29-2, 025 nt; SD = 152 nt). Reads less than 40 nucleotides in length and low quality reads that did not overlap with other reads were discarded. The remaining 955, 105 reads were assembled into 70, 786 contigs with an average length of 836 nt (range: 40-10, 512 nt; SD = 464 nt; Additional File 2; 3). These data are compared with other published scleractinian transcriptomic data sets in Table 1.
Taxonomic affinity of the sequences
Summary of Symbiodinium blast hits
Dna J-like protein 1
Dna J-like protein 2
Dna J-like protein 3
Heat shock protein 70
Heat shock protein 90
Peridinin chlorophyll-a binding protein apoprotein precursor
Ribulose bisphosphate carboxylase
Ubiquitin ligase 2
Ubiquitin-specific protease 1
The taxonomy of Pocillopora is currently regarded as tenuous, e.g., some 16 species have been defined on the basis of morphological features, but cladistic groupings defined by molecular sequence data are not always congruent with these morphologically defined taxa (e.g., ). We used reciprocal blast searches to investigate whether the sequences we generated in this study most closely resemble P. damicornis sequences in NCBI, rather than sequences from other closely related corals, including other Pocillopora species. We obtained 400 "P. damicornis" ESTs obtained from NCBI, and we used these to query the contigs housed at PocilloporaBase (blastn with an E-value cut-off of 0.001). Because a substantial fraction of the P. damicornis sequences currently housed at NCBI represent multiple copies of the same gene (generated in population genetics studies,) many of them matched to the same contig in PocilloporaBase. Overall, the 400 sequences from NCBI matched to 21 unique contigs at PocilloporaBase. We blasted these 21 contigs back against all nucleotide sequences at NCBI classified as scleractinian (Search "Scleractinia[Organism]"). The results are provided in Additional File 6. In 6 of 21 instances, the only match in the database was to a P. damicornis sequence. In 12 instances, there were matches to other corals in addition to P. damicornis, but the highest degree of sequence identity was to a sequence from P. damicornis. In one instance, a contig from PocilloporaBase exhibited equal percent identity to sequences from P. damicornis and P. meandrina. Finally, there were two contigs that exhibited a slightly higher resemblance to a sequence from a coral other than P. damicornis. On balance, these results clearly suggest that the "P. damicornis" populations sampled in this study exhibit greater sequence similarity to the P. damicornis sequences housed at NCBI than to any other coral species represented in that database. However, the taxonomic uncertainty of greatest concern pertains specifically to the genus Pocillopora. Here, the results are less decisive. Only seven of the twenty-one contig sequences we blasted against the NCBI database produced hits to sequences from P. damicornis and another Pocillopora species. In five of these seven instances, the top hit was to a P. damicornis sequence, but one contig matches better to P. molokensis and another matches equally well to P. damicornis and P. meandrina.
KEGG Pathway Analysis
The contigs were subjected to the KEGG Pathway analysis, for human and plant separately . Based on this analysis, the P. damicornis sequences were mapped to metabolic pathways on the interactive tree of life . Components of most metabolic pathways were identified, including photosynthesis (KEGG id: map00195 and 00196) and lipid metabolism. However, some pathways were largely or completely absent from the P. damicornis contigs, including aminosugars metabolism (00530), lipopolysaccharide biosynthesis (00540), peptidoglycan biosynthesis (00550), glycosphingolipid biosynthesis (00602), methane metabolism (00680), androgen and estrogen metabolism (00150), and biodegradation of most xenobiotics (Additional Files 8, 9).
PocilloporaBase: capabilities and functions
PocilloporaBase http://www.PocilloporaBase.org was modeled after StellaBase, a genomic and transcriptomic database for the starlet sea anemone, Nematostella vectensis[33, 34]. Both species-specific databases were designed to integrate with CnidBase , a phylum-wide database meant to facilitate cross-species comparisons among cnidarians. All of the raw sequencing reads generated in this study as well as the assembled contigs can be downloaded there. At present, unlike StellaBase, PocilloporaBase houses only transcriptome data and no genomic sequence data. This repository will grow as future mRNA and genomic sequencing projects generate additional data for Pocillopora damicornis.
The Gene Ontology Search allows you to query sequences using either GO id or GO description terms . Each successful query returns a table that contains the Protein ID and Protein Name, as well as GO id, GO description, and GO type. The Protein ID links to the protein's entry at NCBI. The Protein Name links to a Protein Lookup on PocilloporaBase, which returns a list of P. damicornis contigs that generated significant BLAST hits to the protein in question. The GO id links to the corresponding gene ontology page at the amiGO database. Clicking on the GO description link performs a more specific GO lookup if any gene ontology terms are children of the parent term used to conduct the original search. If too many results are returned, searches can be restricted to one of the three principal gene ontology types: biological process, cellular component, or molecular function. If the user is more interested in the number rather than the identity of the genes in the database that map to a particular gene ontology term, the "Counts only" box can be checked.
The KEGG Pathway Search allows you to query sequences using a KEGG ID number (e.g., hsa00010) or KEGG Pathway Description (e.g., Glycolysis/Gluconeogenesis). Each successful query returns a table with individual contigs identified by their Contig ID, the E-value of their match to a protein in the KEGG database, the KEGG Organism ID, homologous protein name and protein accession ID.
The presence of conserved protein motifs in one or more transcripts can be investigated by searching the data for matches to the conserved protein motifs housed at Pfam. The Pfam Protein Family Classification search allows you to query sequences by Pfam Accession number (e.g., PF00006), Motif Name (e.g.: PAX, actin, DNA_methylase) or Protein Description key words. If there is a match to a conserved protein motif at Pfam, the search returns a table of Pocillopora damicornis contigs encoding that motif sorted by E-value.
It is also possible to search for matches to a query sequence using the complete set of BLAST options. BLAST searches return contig id, the sequence for that contig, as well as the NCBI gene ID, gene name, and the gene sequence for any gene sequence found to match the original blast query. A gene search page allows for quick retrieval of gene and species information in the database.
We used the 454 sequencing platform to generate a reference transcriptome for the cauliflower coral, P. damicornis. A taxonomic analysis of the sequence data indicates that we have captured some of the diversity of the coral holobiont, as many of the sequences appear to be derived from non-metazoan taxa including bacteria, fungi, viruses, and unicellular algae of the genus Symbiodinium. The data have been organized into a publicly available relational database that will be updated and expanded as new P. damicornis sequencing data become available.
Availability and requirements
This database can be accessed using a web browser at http://www.pocilloporabase.org.
Acknowledgements and funding
The authors would like to thank Jonathan Martinez, Luc Rougee, Bob Richmond and Laetitia Hédouin for assistance with permits and animal collections, Derek Stefanik for assistance with database conception, Mark Martindale and members of the Martindale Lab of Kewalo Marine Labs for logistical support and supplies. This research was supported by a grant from the Conservation International Marine Management Area Science Program to JRF. We are grateful to three anonymous reviewers for helping us to significantly improve the clarity of the manuscript and utility of the associated database.
- DeSalvo MK, Sunagawa S, Fisher PL, Voolstra CR, Iglesias-Prieto R, Medina M: Coral host transcriptomic states are correlated with Symbiodinium genotypes. Mol Ecol. 2010, 19 (6): 1174-1186. 10.1111/j.1365-294X.2010.04534.x.View ArticlePubMedGoogle Scholar
- Edge SE, Morgan MB, Gleason DF, Snell TW: Development of a coral cDNA array to examine gene expression profiles in Montastraea faveolata exposed to environmental stress. Mar Pollut Bull. 2005, 51 (5-7): 507-523. 10.1016/j.marpolbul.2005.07.007.View ArticlePubMedGoogle Scholar
- Foret S, Kassahn KS, Grasso LC, Hayward DC, Iguchi A, Ball EE, Miller DJ: Genomic and microarray approaches to coral reef conservation biology. Coral Reefs. 2007, 26 (3): 475-486. 10.1007/s00338-007-0206-1.View ArticleGoogle Scholar
- Grasso LC, Maindonald J, Rudd S, Hayward DC, Saint R, Miller DJ, Ball EE: Microarray analysis identifies candidate genes for key roles in coral development. BMC Genomics. 2008, 9: 540-10.1186/1471-2164-9-540.View ArticlePubMedPubMed CentralGoogle Scholar
- Bay LK, Ulstrup KE, Nielsen HB, Jarmer H, Goffard N, Willis BL, Miller DJ, Van Oppen MJ: Microarray analysis reveals transcriptional plasticity in the reef building coral Acropora millepora. Mol Ecol. 2009, 18 (14): 3062-3075. 10.1111/j.1365-294X.2009.04257.x.View ArticlePubMedGoogle Scholar
- Levy O, Kaniewska P, Alon S, Eisenberg E, Karako-Lampert S, Bay LK, Reef R, Rodriguez-Lanetty M, Miller DJ, Hoegh-Guldberg O: Complex diel cycles of gene expression in coral-algal symbiosis. Science. 2011, 331 (6014): 175-10.1126/science.1196419.View ArticlePubMedGoogle Scholar
- Putnam NH, Srivastava M, Hellsten U, Dirks B, Chapman J, Salamov A, Terry A, Shapiro H, Lindquist E, Kapitonov VV, et al: Sea anemone genome reveals ancestral eumetazoan gene repertoire and genomic organization. Science. 2007, 317 (5834): 86-94. 10.1126/science.1139158.View ArticlePubMedGoogle Scholar
- Jackson JB: Colloquium paper: ecological extinction and evolution in the brave new ocean. Proc Natl Acad Sci USA. 2008, 105 (Suppl 1): 11458-11465.View ArticlePubMedPubMed CentralGoogle Scholar
- Pandolfi JM, Jackson JB, Baron N, Bradbury RH, Guzman HM, Hughes TP, Kappel CV, Micheli F, Ogden JC, Possingham HP, et al: Ecology. Are U.S. coral reefs on the slippery slope to slime?. Science. 2005, 307 (5716): 1725-1726. 10.1126/science.1104258.View ArticlePubMedGoogle Scholar
- Sandin SA, Smith JE, Demartini EE, Dinsdale EA, Donner SD, Friedlander AM, Konotchick T, Malay M, Maragos JE, Obura D, et al: Baselines and degradation of coral reefs in the Northern Line Islands. PLoS One. 2008, 3 (2): e1548-10.1371/journal.pone.0001548.View ArticlePubMedPubMed CentralGoogle Scholar
- Ostrander GK, Armstrong KM, Knobbe ET, Gerace D, Scully EP: Rapid transition in the structure of a coral reef community: the effects of coral bleaching and physical disturbance. Proc Natl Acad Sci USA. 2000, 97 (10): 5297-5302. 10.1073/pnas.090104897.View ArticlePubMedPubMed CentralGoogle Scholar
- Mydlarz LD, Holthouse SF, Peters EC, Harvell CD: Cellular responses in sea fan corals: granular amoebocytes react to pathogen and climate stressors. PLoS One. 2008, 3 (3): e1811-10.1371/journal.pone.0001811.View ArticlePubMedPubMed CentralGoogle Scholar
- Brown BE: Coral bleaching: causes and consequences. Coral Reefs. 1997, 16: S129-S138. 10.1007/s003380050249.View ArticleGoogle Scholar
- Van Oppen MJH, Gates RD: Conservation genetics and the resilience of reef-building corals. Mol Ecol. 2006, 15 (13): 3863-3883. 10.1111/j.1365-294X.2006.03026.x.View ArticlePubMedGoogle Scholar
- Green EP, Bruckner AW: The significance of coral disease epizootiology for coral reef conservation. Biol Conserv. 2000, 96 (3): 347-361. 10.1016/S0006-3207(00)00073-2.View ArticleGoogle Scholar
- Reitzel AM, Sullivan JC, Traylor-Knowles N, Finnerty JR: Genomic survey of candidate stress-response genes in the estuarine anemone Nematostella vectensis. Biol Bull. 2008, 214 (3): 233-254. 10.2307/25470666.View ArticlePubMedGoogle Scholar
- Schwarz JA, Brokstein PB, Voolstra C, Terry AY, Manohar CF, Miller DJ, Szmant AM, Coffroth MA, Medina M: Coral life history and symbiosis: functional genomic resources for two reef building Caribbean corals, Acropora palmata and Montastraea faveolata. BMC Genomics. 2008, 9: 97-10.1186/1471-2164-9-97.View ArticlePubMedPubMed CentralGoogle Scholar
- Sunagawa S, DeSalvo MK, Voolstra CR, Reyes-Bermudez A, Medina M: Identification and gene expression analysis of a taxonomically restricted cysteine-rich protein family in reef-building corals. PLoS One. 2009, 4 (3): e4865-10.1371/journal.pone.0004865.View ArticlePubMedPubMed CentralGoogle Scholar
- Voolstra CR, Schwarz JA, Schnetzer J, Sunagawa S, Desalvo MK, Szmant AM, Coffroth MA, Medina M: The host transcriptome remains unaltered during the establishment of coral-algal symbioses. Mol Ecol. 2009, 18 (9): 1823-1833. 10.1111/j.1365-294X.2009.04167.x.View ArticlePubMedGoogle Scholar
- Meyer E, Aglyamova GV, Wang S, Buchanan-Carter J, Abrego D, Colbourne JK, Willis BL, Matz MV: Sequencing and de novo analysis of a coral larval transcriptome using 454 GSFlx. BMC Genomics. 2009, 10: 219-10.1186/1471-2164-10-219.View ArticlePubMedPubMed CentralGoogle Scholar
- Portune KJ, Voolstra CR, Medina M, Szmant AM: Development and heat stress-induced transcriptomic changes during embryogenesis of the scleractinian coral Acropora palmata. Mar Genom. 2010, 3 (1): 51-62. 10.1016/j.margen.2010.03.002.View ArticleGoogle Scholar
- Kerr AM: Molecular and morphological supertree of stony corals (Anthozoa: Scleractinia) using matrix representation parsimony. Biol Rev Camb Philos Soc. 2005, 80 (4): 543-558. 10.1017/S1464793105006780.View ArticlePubMedGoogle Scholar
- Hunter CL, Evans CW: Coral-Reefs in Kaneohe Bay, Hawaii - 2 Centuries of Western Influence and 2 Decades of Data. B Mar Sci. 1995, 57 (2): 501-515.Google Scholar
- Grigg RW: Coral reefs in an urban embayment in Hawaii: A complex case history controlled by natural and anthropogenic stress. Coral Reefs. 1995, 14 (4): 253-266.View ArticleGoogle Scholar
- Rio DC, Ares M, Hannon GJ, Nilsen TW: Purification of RNA using TRIzol (TRI reagent). Cold Spring Harb Protoc. 2010, 2010 (6): pdb prot5439Google Scholar
- Rothberg JM, Leamon JH: The development and impact of 454 sequencing. Nat Biotechnol. 2008, 26 (10): 1117-1124. 10.1038/nbt1485.View ArticlePubMedGoogle Scholar
- Mardis ER: Next-generation DNA sequencing methods. Annu Rev Genomics Hum Genet. 2008, 9: 387-402. 10.1146/annurev.genom.9.081307.164359.View ArticlePubMedGoogle Scholar
- Chevreux B, Pfisterer T, Drescher B, Driesel AJ, Muller WE, Wetter T, Suhai S: Using the miraEST assembler for reliable and automated mRNA transcript assembly and SNP detection in sequenced ESTs. Genome Res. 2004, 14 (6): 1147-1159. 10.1101/gr.1917404.View ArticlePubMedPubMed CentralGoogle Scholar
- Finn RD, Mistry J, Tate J, Coggill P, Heger A, Pollington JE, Gavin OL, Gunasekaran P, Ceric G, Forslund K: The Pfam protein families database. Nucleic Acids Res. 2010, D211-222. 38 DatabaseGoogle Scholar
- Kanehisa M: The KEGG database. Novartis Found Symp. 2002, 247: 91-101. discussion 101-103, 119-128, 244-152View ArticlePubMedGoogle Scholar
- Pinzon JH, LaJeunesse TC: Species delimitation of common reef corals in the genus Pocillopora using nucleotide sequence phylogenies, population genetics and symbiosis ecology. Mol Ecol. 2011, 20 (2): 311-325. 10.1111/j.1365-294X.2010.04939.x.View ArticlePubMedGoogle Scholar
- Letunic I, Yamada T, Kanehisa M, Bork P: iPath: interactive exploration of biochemical pathways and networks. Trends Biochem Sci. 2008, 33 (3): 101-103. 10.1016/j.tibs.2008.01.001.View ArticlePubMedGoogle Scholar
- Sullivan JC, Reitzel AM, Finnerty JR: Upgrades to StellaBase facilitate medical and genetic studies on the starlet sea anemone, Nematostella vectensis. Nucleic Acids Res. 2008, D607-611. 36 DatabaseGoogle Scholar
- Sullivan JC, Ryan JF, Watson JA, Webb J, Mullikin JC, Rokhsar D, Finnerty JR: StellaBase: the Nematostella vectensis Genomics Database. Nucleic Acids Res. 2006, D495-499. 34 DatabaseGoogle Scholar
- Ryan JF, Finnerty JR: CnidBase: The Cnidarian Evolutionary Genomics Database. Nucleic Acids Res. 2003, 31 (1): 159-163. 10.1093/nar/gkg116.View ArticlePubMedPubMed CentralGoogle Scholar
- The Gene Ontology Project. [http://www.geneontology.org/]
- Wang S, Zhang L, Meyer E, Matz MV: Construction of a high-resolution genetic linkage map and comparative genome analysis for the reef-building coral Acropora millepora. Genome Biol. 2009, 10 (R216):Google Scholar
This article is published under license to 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.