Augmented annotation and orthologue analysis for Oryctolagus cuniculus: Better Bunny
© Craig et al.; licensee BioMed Central Ltd. 2012
Received: 13 December 2011
Accepted: 25 April 2012
Published: 8 May 2012
The rabbit is an important model organism used in a wide range of biomedical research. However, the rabbit genome is still sparsely annotated, thus prohibiting extensive functional analysis of gene sets derived from whole-genome experiments. We developed a web-based application that provides augmented annotation and orthologue analysis for rabbit genes. Importantly, the application allows comprehensive functional analysis through the use of orthologous relationships.
Using data extracted from several public bioinformatics repositories we created Better Bunny, a database and query tool that extensively augments the available functional annotation for rabbit genes. Using the complete set of target genes from a commercial rabbit gene expression microarray as our benchmark, we are able to obtain functional information for 88 % of the genes on the microarray. Previously, functional information was available for fewer than 10 % of the rabbit genes.
We have developed a freely available, web-accessible bioinformatics tool that enables investigators to quickly and easily perform extensive functional analysis of rabbit genes (http://cptweb.cpt.wayne.edu). The software application fills a critical void for a wide range of biomedical research that relies on the rabbit model and requires characterization of biological function for large sets of genes.
KeywordsRabbit Genome Annotation Orthologue Microarray Functional analysis
The rabbit (Oryctolagus cuniculus) is an important model organism used extensively in biomedical research. Offering advantages over other organisms for a variety of physiological systems, the rabbit has contributed to studies in a wide range of disciplines, including embryology, ophthalmology, toxicology, pulmonary and cardiovascular research, and neurology [1–4]. For example, the rabbit has proven to be invaluable in neurodevelopmental research, largely due to a temporal pattern of oligodendrocyte maturation and myelination that closely parallels that of humans. In rabbit brain development structural and functional changes occur most rapidly during the perinatal period, generally starting several days before birth and continuing in the postnatal period . This pattern is very similar to what occurs in humans, where immature oligodendrocytes increase rapidly in number in the third trimester followed by myelination that initiates around the time of birth, with maximum increases during the first year of life . Thus, structural development of the rabbit brain is very similar to that of humans, but within a more compressed time-frame. In contrast, rodents are postnatal brain developers, while sheep, pigs and monkeys are prenatal brain developers. Consequently, the rabbit offers a preferred means to investigate the effects of perinatal events and environmental exposures on brain development. This model has been used extensively in one of our own laboratories to study the relationship between maternal infection and neonatal brain injury [7, 8], and by other investigators to study the effects of various prenatal and perinatal insults to the brain .
Despite the demonstrated importance of rabbits in biomedical research, the rabbit genome is still poorly annotated. This presents a significant challenge where whole-genome technologies are used to characterize molecular and cellular events. Using high-throughput genomic methods such as microarrays, investigators often obtain lengthy lists of differentially expressed genes, and a commonly employed analysis paradigm is to perform functional annotation analysis to identify enriched ontologies, pathways, motifs, and other molecular details associated with a set of genes. A variety of software is available for functional analysis of common model organisms; however, this approach is currently impeded when using rabbit data due to sparse annotation that is scattered among various bioinformatics resources. Only 1076 rabbit genes are annotated with biological process information per the current database release of the Gene Ontology Consortium (GO) . This compares to more than 27,000 genes with an annotated biological process in human, and approximately 16,000 and 24,000 genes in rat and mouse respectively. Additionally, the rabbit is not well represented in other informatics resources. For example, rabbit is not included in the Kyoto Encyclopedia of Genes and Genomes (KEGG) which is widely used for biological pathway analysis . These limitations are further exacerbated by the absence of rabbit in Homologene , an NCBI resource that provides orthologous relationships among model organisms. Such relationships can be effectively used to infer functional details for genes having little annotation, where an orthologue is better characterized in another organism. In the case of rabbit, one could utilize the well-annotated and highly homologous genes from human, mouse, and rat for inferred functional annotation. However, at this time investigators who depend on the rabbit model for biomedical studies are confronted by the obstacle of a time-consuming informatics effort to assemble pertinent annotation and orthologue data, consequently they risk overlooking potentially important findings in their experiments.
The need for a freely available tool to perform functional analysis of rabbit genes is highlighted by a recent report that used a rabbit model to study the global gene expression response to Mycobacterium tuberculosis infection and investigate the influence of immune system modulation on treatment efficacy . The authors used Agilent rabbit gene expression microarrays to measure genome-wide expression changes in lung tissues. After obtaining lists of differentially expressed genes through statistical analysis of the microarray data, the authors sought to perform functional classification of the identified genes. However, they noted that functional annotation of most rabbit genes is currently unavailable. So they performed extensive computational work to extract orthologous relationships from a commercial bioinformatics database, which were then used to infer the function of rabbit genes based on annotation of their human, mouse, and rat counterparts. While the authors of the study were successful using this approach, many research groups do not have access to such proprietary data and/or the computational expertise to derive a functional annotation database. To enable investigators to perform functional analysis of rabbit genes, we have created a web-accessible functional annotation tool using data derived from public repositories. The application is easy to use and freely available.
To address the paucity and dispersed state of available rabbit annotations, with the ultimate goal of functional analysis, we have created a web-based application for augmented annotation and analysis of rabbit genes: “Better Bunny.” Access to the web server is freely available at http://cptweb.cpt.wayne.edu. The server provides annotation and orthologous relationships assembled from public bioinformatics resources, including NCBI  and Ensembl  databases. Importantly, comprehensive pathway, ontology, and functional analysis can be instantly performed on rabbit gene lists using orthologous relationships provided through Better Bunny.
One of the most comprehensive and widely used applications for functional annotation analysis is the Database for Annotation, Visualization and Integrated Discovery (DAVID) [17, 18]. This well-maintained resource provides extensive functional analysis of input gene sets, including enrichment of pathways and ontologies, annotation clustering, and identification of protein interactions and domains. We have seamlessly integrated DAVID analysis into the Better Bunny application. The output from Better Bunny includes a “David Annotations” link at the top of the Ensembl Gene ID column for each species selected. This link automatically submits the gene list from the selected species to DAVID and opens a new browser window with the results. Adjacent to the “David Annotations” link is a “?” that provides information on using DAVID for functional annotation. Readers are encouraged to explore the ample help available at the DAVID web site for additional information on this valuable resource. While very little functional information is typically obtained when rabbit gene identifiers are submitted to DAVID, using orthologues via Better Bunny quickly provides a rich source of inferred functional information for most rabbit genes in a list.
Results and discussion
To assess the value of Better Bunny in the context of a whole-genome experiment, we examined the annotations associated with the entire set of probes on the Agilent rabbit gene expression oligonucleotide microarray. There are 43,603 probes on the microarray, with many genes targeted by redundant probes. The following discussion reflects annotation of the unique set of target genes. Of the 12,118 unique genes represented on the microarray (based on Ensembl gene ID), only 916 genes have a gene name per the vendor’s annotation, 2324 have a Unigene identifier, and no gene ontology information is provided. The remaining genes include many uncharacterized clones. Using the augmented annotation available through Better Bunny, more than 10,000 genes are provided a name, largely through Ensembl resources, and 9470 genes have a molecular function assigned from the Gene Ontology Consortium, most established by Better Bunny through orthologous relationships.
Of great utility, Better Bunny identifies 11,026 putative human orthologues with at least 50 % identity to the corresponding rabbit sequence, and a similar number are available for mouse and rat. This enables instant and extensive functional annotation for most genes on the microarray via the integrated DAVID analysis available in Better Bunny. To demonstrate the functional information that is gained through Better Bunny, we submitted the probe list for all genes on the Agilent microarray to Better Bunny and specified the output to include human orthologues having a minimum of 50 % identity. The set of 11,026 orthologous human genes was then submitted to DAVID using the integrated link available in Better Bunny. Using default settings in DAVID, functional annotation was obtained for 10,731 genes, representing a wide range of functional information. For example, SwissProt (UniProt) comments were obtained for 10,221 genes; Gene Ontology Biological Process annotation was obtained for 8285 genes; 3184 genes were associated with KEGG pathways; and Interpro protein classification information was obtained for 9697 genes. This compares to only one rabbit gene with functional annotation when rabbit gene identifiers are submitted to DAVID. The dramatic difference in results obtained for the two species reflects the scarcity of rabbit gene annotation available in the underlying databases, as exemplified by the GO biological process annotation for rabbit which currently covers only approximately 1,000 genes, and the results demonstrate the benefit of using orthologous relationships for functional annotation.
In summary, we have developed a freely available and web-accessible bioinformatics tool that enables investigators to quickly and easily perform extensive functional characterization of lengthy lists of rabbit genes. The software application fills a critical void for investigators who employ the rabbit model in their research and who wish to characterize the biological function and cellular role associated with sets of genes. Using the complete set of target genes from a commercial rabbit microarray, we were able to obtain functional information for 88 % of the genes, whereas functional information was previously available for fewer than 10 % of the rabbit genes on the microarray. The augmented gene annotation, orthologue identification, and integrated functional analysis available through Better Bunny are expected to greatly enhance the knowledge gained from a wide variety of biomedical research projects using the rabbit model.
Availability and requirements
Project name: Better Bunny
Project home page: http://cptweb.cpt.wayne.edu
Operating system(s): Platform independent
Programming language: PHP / Python / MySQL
Other requirements: none
Any restrictions to use by non-academics: none
We would like to thank Dr. David Svinarich for helpful suggestions and Brad Sherman for assistance with DAVID integration. This work was supported in part by the National Institutes of Health NICHD [K08 HD50562-01A1 to S.K.].
- Puschel B, Daniel N, Bitzer E, Blum M, Renard JP, Viebahn C: The rabbit (Oryctolagus cuniculus): a model for mammalian reproduction and early embryology. Cold Spring Harb Protoc 2010. 2010:pdb emo139Google Scholar
- Popp MP, Liu L, Timmers A, Esson DW, Shiroma L, Meyers C, Berceli S, Tao M, Wistow G, Schultz GS, Sherwood MB: Development of a microarray chip for gene expression in rabbit ocular research. Mol Vis 2007, 13: 164–173.PubMed CentralPubMedGoogle Scholar
- Coico R, Woodruff-Pak DS: Immunotherapy for Alzheimer's disease: harnessing our knowledge of T cell biology using a cholesterol-fed rabbit model. J Alzheimers Dis 2008, 15: 657–671.PubMedGoogle Scholar
- Subbian S, Tsenova L, O'Brien P, Yang G, Koo MS, Peixoto B, Fallows D, Dartois V, Muller G, Kaplan G: Phosphodiesterase-4 inhibition alters gene expression and improves isoniazid-mediated clearance of Mycobacterium tuberculosis in rabbit lungs. PLoS Pathog 2011, 7: e1002262. 10.1371/journal.ppat.1002262PubMed CentralView ArticlePubMedGoogle Scholar
- Drobyshevsky A, Song SK, Gamkrelidze G, Wyrwicz AM, Derrick M, Meng F, Li L, Ji X, Trommer B, Beardsley DJ, et al.: Developmental changes in diffusion anisotropy coincide with immature oligodendrocyte progression and maturation of compound action potential. J Neurosci 2005, 25: 5988–5997. 10.1523/JNEUROSCI.4983-04.2005View ArticlePubMedGoogle Scholar
- Kinney HC, Brody BA, Kloman AS, Gilles FH: Sequence of central nervous system myelination in human infancy. II. Patterns of myelination in autopsied infants. J Neuropathol Exp Neurol 1988, 47: 217–234. 10.1097/00005072-198805000-00003View ArticlePubMedGoogle Scholar
- Kannan S, Saadani-Makki F, Balakrishnan B, Dai H, Chakraborty PK, Janisse J, Muzik O, Romero R, Chugani DC: Decreased cortical serotonin in neonatal rabbits exposed to endotoxin in utero. J Cereb Blood Flow Metab 2011, 31: 738–749. 10.1038/jcbfm.2010.156PubMed CentralView ArticlePubMedGoogle Scholar
- Kannan S, Saadani-Makki F, Muzik O, Chakraborty P, Mangner TJ, Janisse J, Romero R, Chugani DC: Microglial activation in perinatal rabbit brain induced by intrauterine inflammation: detection with 11 C-(R)-PK11195 and small-animal PET. J Nucl Med 2007, 48: 946–954. 10.2967/jnumed.106.038539View ArticlePubMedGoogle Scholar
- Derrick M, Luo NL, Bregman JC, Jilling T, Ji X, Fisher K, Gladson CL, Beardsley DJ, Murdoch G, Back SA, Tan S: Preterm fetal hypoxia-ischemia causes hypertonia and motor deficits in the neonatal rabbit: a model for human cerebral palsy? J Neurosci 2004, 24: 24–34. 10.1523/JNEUROSCI.2816-03.2004View ArticlePubMedGoogle Scholar
- The Gene Ontology Consortium[http://www.geneontology.org/]
- Kyoto Encyclopedia of Genes and Genomes[http://www.genome.jp/kegg/]
- National Center for Biotechnology Information[http://www.ncbi.nlm.nih.gov/]
- Vilella AJ, Severin J, Ureta-Vidal A, Heng L, Durbin R, Birney E: EnsemblCompara GeneTrees: Complete, duplication-aware phylogenetic trees in vertebrates. Genome Research 2009, 19: 327–335.PubMed CentralView ArticlePubMedGoogle Scholar
- Sangar V, Blankenberg DJ, Altman N, Lesk AM: Quantitative sequence-function relationships in proteins based on gene ontology. BMC Bioinformatics 2007, 8: 294. 10.1186/1471-2105-8-294PubMed CentralView ArticlePubMedGoogle Scholar
- Huang DW, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protocols 2008, 4: 44–57. 10.1038/nprot.2008.211View ArticleGoogle Scholar
- The Database for Annotation, Visualization and Integrated Discovery (DAVID)[http://david.abcc.ncifcrf.gov/]
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.