FunnyBase: a systems level functional annotation of Fundulus ESTs for the analysis of gene expression
© Paschall et al; licensee BioMed Central Ltd. 2004
Received: 17 August 2004
Accepted: 20 December 2004
Published: 20 December 2004
While studies of non-model organisms are critical for many research areas, such as evolution, development, and environmental biology, they present particular challenges for both experimental and computational genomic level research. Resources such as mass-produced microarrays and the computational tools linking these data to functional annotation at the system and pathway level are rarely available for non-model species. This type of "systems-level" analysis is critical to the understanding of patterns of gene expression that underlie biological processes.
We describe a bioinformatics pipeline known as FunnyBase that has been used to store, annotate, and analyze 40,363 expressed sequence tags (ESTs) from the heart and liver of the fish, Fundulus heteroclitus. Primary annotations based on sequence similarity are linked to networks of systematic annotation in Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) and can be queried and computationally utilized in downstream analyses. Steps are taken to ensure that the annotation is self-consistent and that the structure of GO is used to identify higher level functions that may not be annotated directly. An integrated framework for cDNA library production, sequencing, quality control, expression data generation, and systems-level analysis is presented and utilized. In a case study, a set of genes, that had statistically significant regression between gene expression levels and environmental temperature along the Atlantic Coast, shows a statistically significant (P < 0.001) enrichment in genes associated with amine metabolism.
The methods described have application for functional genomics studies, particularly among non-model organisms. The web interface for FunnyBase can be accessed at http://genomics.rsmas.miami.edu/funnybase/super_craw4/. Data and source code are available by request at firstname.lastname@example.org.
Investigating patterns of gene expression using mouse and human microarrays has produced insights into cancer [1, 2], cardiac diseases [3–6], and metabolic disorders [7–12]. These and many other functional genomics studies rely on full genomic sequence to establish well-annotated databases. Yet, microarrays based on EST collections are increasingly being used for diverse species, from honey bees to fish [13–20] and including simple diploblastic organisms . These studies within a diversity of organisms provide insights not provided by 'model' species (species that are genetically well defined or with annotated genomes ). For example, 'non-model' organisms have provided insight into the natural variation in gene expression , social castes among bees [24, 25], hypoxia , and physiological responses to variation in the thermal environment [27, 28]. To investigate adaptive variation in gene expression we use the teleost Fundulus heteroclitus (killifish) [23, 29].
The killifish Fundulus heteroclitus are distributed along the eastern coast of North America which has one of the steepest thermal clines in the world: northern populations have environmental temperatures more than 12°C below southern populations across 12 degrees of latitude. Migration among populations is sufficient to minimize random genetic drift  but not frequent enough to extinguish local adaptation [31, 32]. Populations are large (>10,000) and affected by historical, demographic and selective constraints, providing a framework for the partitioning of variation in gene expression within and among populations. Additionally, the well-established phylogenetic relationship among Fundulus species can be used to discern adaptive changes [23, 33, 34]. These characteristics make F. heteroclitus an ideal species to investigate adaptive variation in gene expression.
Microarrays from diverse EST collections offer opportunities to address many biological problems, but to effectively use this information often requires a locally generated bioinformatics approach. Tools like the TIGR Gene index  and Unigene  provide significant information on many species, yet these databases do not meet the needs of functional genomics projects for many non-model species. Currently, TIGR and NCBI provide gene indices for 28 and 23 animal species, respectively. Yet, there are 63 animal species with more than 10,000 ESTs . The number of species with ESTs >10,000 has continued to grow, and there was approximately a 20% increase in the preceding three months. While annotation from these resources can be accessed through web-based homology searches, for many laboratory collections of ESTs it is difficult to use existing tools to achieve a systems-level view of gene functions and relationships. Rather than simply browsing functional information over the web for a different group's project, laboratories that produce novel EST collections and microarrays require customized databases providing access to integrated functional annotation as expression data are being analyzed.
We have developed FunnyBase to meet these functional genomics needs. FunnyBase provides functional information for >40,000 ESTs from the teleost fish Fundulus heteroclitus, provides the means to quickly process, evaluate, and store annotation based on similarity searches of public resources, and integrates these data with species-specific clustering and microarray analysis. Perhaps ironically, the greatest challenge for functional annotation based on similarity searches is an overabundance of data. There are a number of databases to chose from, and often the single best hit from a given database search is not the most informative. FunnyBase implements a strategy to make maximum use of systems-level functional information from Gene Ontology (GO)  assignments and membership in metabolic pathways as defined by the Kyoto Encyclopedia of Genes and Genomes (KEGG) . Specifically, several sequence databases are queried and results integrated to maximize the number of annotated sequences. Alignments and scores for all homology based associations are tracked, allowing further evaluation and statistical studies.
Microarray data using genes annotated in FunnyBase can be systematically analyzed in the context of biological functions. We present a case study to illustrate how assessment of systems-level annotation can identify statistically significant functional differences among sets of genes.
Results and discussion
The Hierarchical Annotation module uses the consensus sequences from the clusters or singletons, and integrates primary annotation such as gene name and description with associated pathways and systems-level functional annotation. This may include gene function (e.g., enzyme catalyst), metabolic or signal pathway (e.g., oxidative phosphorylation), or biological function (e.g., protein translation). Sequence data from the first module and functional annotation from the second are matched using database similarity searches (BLASTX and BLASTN). E-values, bit scores and local alignments are stored in the 'Similarity Data' table for all significant matches. One of the strengths of FunnyBase is the use of different sequence databases (SwissProt , NCBI UniGene , and NCBI non-redundant NR ) to provide separate annotations. Although these databases are not completely independent, the three separate annotations provide verification of gene names.
The third module, Microarray Production & Analysis provides a list of unique genes to be printed and integrates expression data from microarray experiments with the Hierarchical Annotation module. This provides functional annotation for expression data. FunnyBase annotation is accessible through the web or through local SQL queries and data-mining scripts.
EST isolation and sequencing
The overall strategy used to isolate and sequence thousands of Fundulus cDNAs was (1) generate a high quality unidirectional cDNA library, (2) normalize the library, (3) randomly pick colonies and amplify by PCR the cDNA within the vector, (4) sequence and identify PCR products, and (5) after approximately every 1,000 clones, subtract these from the normalized library and repeat steps 3–5. Details for all protocols are provided at http://crawford.rsmas.miami.edu/ and were used in the Comparative Functional Genomic course at Mount Desert Island, ME 2000.
We sequenced 46,433 ESTs and 40,043 of these are available in the dbEST database at NCBI (dbEST identification numbers: 23,480,307 to 23,515,306; 23,520,047 to 23,525,409 and 24,320,128 to 24,320,184) as of June 26, 2004. Sequences in FunnyBase are identified by a number series: unique sequence number, array number, plate number and well identifier (example: 23434_125_001_H04). The remaining 6,966 un-submitted sequences failed to meet one of the sequence quality parameters. Two criteria are used for defining "good" sequences: 1) >100 bp of sequence with Phred score >20 or 2) form an overlapping cluster with other sequences. Of the 19,937 sequences processed with the current version of FunnyBase, 17,893 (90%) passed one of these quality parameters. In earlier iterations of FunnyBase, visual inspection and later a sequencer-specific quality measure equivalent to a Phred score of 15 were used as filters resulting in 3,603 of the first 5,760 sequences (63%) and 13,922 of the next 15,168 sequences (92%) meeting quality standards, respectively. Re-sequences account for 5,668 sequences, and 4,625 of these were submitted.
One of the most important steps for producing microarrays from cDNA libraries is being able to associate the bacteria containing the cDNAs of interest with the EST annotation. High-throughput procedures are highly prone to tracking errors including: loading plates into an automatic sequencer in the wrong order, orienting symmetric plates in the wrong direction, or mislabeling of plates. The ability to identify these types of mistakes requires controls for identifying plates and plate orientation. The FunnyBase system has a number of integrated quality control steps. First, a Ctenophore cDNA (NCBI: accession number: CN992733) that is unlike anything else in GenBank is used as a control. Controls are placed in 96-well plates in wells corresponding to the plate number and two orientation wells (A5 and F9). Sequences from 96-well plates are automatically scanned for these controls so that the identity and orientation are confirmed and a report is generated for manual review. Secondly, all clones used for microarray production are re-sequenced. This is necessary because individual cDNAs are cherry picked, re-grown and re-amplified, and each of these steps has the potential to introduce or magnify an error. For example, for a 6,000 gene array, a 5% error rate would result in 300 incorrect clones. Re-sequenced array plates are compared using pair-wise BLAST  against previous sequencing results so that the identity of printed microarray spots are verified.
The Ten Most Frequently Annotated ESTs. Clusters with the greatest number of annotated ESTs, the sequence id, number of ESTs that cluster together, and the e-value (probability of similarity) are listed.
Number of ESTs
Vitellogenin I precursor
Cytochrome c oxidase polypeptide II
Alpha-1-antitrypsin homolog precursor
Zona pellucida sperm-binding protein 3 precursor
Cytochrome c oxidase polypeptide I
Zona pellucida sperm-binding protein B precursor
Apolipoprotein A-I precursor
Cytochrome c oxidase polypeptide III
ATP synthase a chain
FunnyBase contains a total of 40,043 EST sequences from F. heteroclitus heart and liver. Clustering with CAP3 yields 3,776 clusters that contain 30,688 ESTs (77%). The remaining 8,991 ESTs (23%) are singletons. By storing the results of clustering with annotation, FunnyBase easily identifies these genes and aids in the selection of genes to be used for library subtraction with the goal of picking less common transcripts.
The 10 annotated clusters with the most sequences are listed in Table 1. In microarray experiments these genes tend to be highly expressed and the fluorescent signal tends to saturate the photomultiplier tube. These genes also serve to verify microarray printing because the predicted spots for these genes have the strongest signal.
One of the objectives of EST projects is to isolate most, if not all, genes expressed in a tissue or organism. The increasing size of larger clusters with more sequencing efforts indicates that strategies to increase the probability of isolating new genes need to be employed. We used two strategies. First, we normalized the library to reduce the differences among expressed genes to less than 10-fold among rare and abundant mRNAs [44, 45]. Using this technique we were able to reduce redundancy in annotated genes from 33% in the non-normalized library to 11% after normalization. Second, we targeted specific sequences for subtraction: annotated cDNAs with high frequencies were targeted in order to focus effort on picking new, rare sequences. Through subtraction we were able to increase the rate of discovery of new annotatable sequences from 24% to 36%. However, analysis of these results indicate that a set of highly expressed sequences, some of which were not subtracted because they were not in the set of annotated genes, still make up much of the EST library and should be the focus of future subtractions.
Annotations for ESTs are based on similarity using BLASTX or BLASTN  to sequences in one of six-public databases: Swiss-Prot, Human UniGene, Danio rerio UniGene, Oncorhynchus mykiss Unigene, Oryzias latipes Unigene, or GenBank NR. FunnyBase includes locally parsed copies of these databases in a relational format. Thus, all searches are done locally and annotation features beyond the FASTA description can be queried. Consensus sequences from the Fundulus EST clusters as well as high quality single unique sequences (singletons) are used as query sequences for BLAST searches against these public databases. The use of consensus sequences, when available, allows sequences that do not contain regions of significant similarity with known protein (e.g., 5' or 3' noncoding regions) to be annotated if they are members of an annotated cluster. All hits with e-value less than 0.001 and their associated alignments are stored in the database and tracked with any associated functional annotation. Users can specify a custom level of significance when assessing the validity of homology based annotation. This record, which goes beyond storing a certain number of 'best hits', is critical because in many cases additional results may have a negligibly lower alignment score, but provide much more useful functional data.
Systematic functional annotation: KEGG and Gene Ontology
In conjunction with performing similarity searches by BLAST, FunnyBase includes locally parsed representations of public databases such as SWISS-PROT in a relational database format. These databases provide additional information that cross-references other public resources such as GO, KEGG or OMIM  that are not available in the single FASTA description line returned by BLAST search.
KEGG  is a unique tool that represents metabolic and signal-transduction pathways both visually and computationally. FunnyBase links annotated genes to enzymes in KEGG pathways based on enzyme commission (EC) numbers. These pathway associations are stored and queries can readily identify genes from a given pathway that show specific patterns of expression. For visual inspection of the pathway, the web interface http://genomics.rsmas.miami.edu/funnybase/super_craw4/ links directly to the graphical KEGG pathways in which a gene occurs.
Number of distinct sequences in the Top 10 most common KEGG pathways. The KEGG pathway name and number of distinct sequences (clusters or singletons) from FunnyBase are presented. There are more distinct sequences than enzymes in a pathway because many enzymes have several protein subunits and many proteins have several different loci encoding the same subunit (e.g., NADH dehydrogenase, a protein complex of oxidative phosphorylation, has 26 protein subunits and 42 loci for these subunits).
Sequence count for TOP 10 pathways
Fatty acid metabolism
Valine, leucine and isoleucine degradation
Glycine, serine and threonine metabolism
The Gene Ontology project (GO) has produced a structured vocabulary in the form of an acyclic directed graph that biologists can use to annotate genes in a systematic manner . FunnyBase includes two non-trivial steps to make the best possible use of GO terms. First, many GO annotations are lost if only the single 'best hit' from a homology search is considered because GO annotation is applied most often to a few model species such as human that may not appear as the single 'best hit' in a list of BLAST results. FunnyBase identifies the gene name associated with the 'best hit' BLAST result and then uses all GO annotation associated with hits from the complete BLAST results that have the same gene name as the 'best hit' and an e-value of e < 10-12. The goal of this approach is to identify annotation associated with a single 'best hit' gene based on results that may come from multiple species (orthologous genes) and therefore may have varying degrees of sequence similarity due to phylogenetic distance, but to avoid the problem of selecting an inconsistent set of GO terms arising from gene families that share regions of sequence similarity but may have different functions.
Secondly, GO annotation in public databases tends to annotate sequences with only the most specific GO term available, for example RNA polymerase II transcription factor activity, enhancer binding (GO:0003705) rather than the more general parent term transcription regulator activity (GO:0030528). However, in functional genomic analysis, significant patterns of expression may exist at the more general level of functional description. FunnyBase takes advantage of the connected parent-child relationship of GO terms provided by using the relational database version of GO available for download at http://www.geneontology.org to identify such relationships. These data are used to extract the tree of more general GO terms related to those provided by public databases. A FunnyBase script then re-annotates genes with this more complete set of GO terms.
Of the 3,877 annotated genes, 1,912 (54%) are assigned one or more GO terms with a total 6,728 GO assignments being made directly based on information in public databases such as SwissProt. Using parent-child GO term relationship backtracking, an additional 34,112 GO term assignments were made, resulting in a final count of GO assignments of 36,024 excluding the most general terms that divide GO into three categories. Thus, on average, 19 GO terms are assigned to each of 1,912 annotated genes.
Gene scaffolding: clustering of clusters
Humans have approximately 30,000 expressed genes, yet there are over 1,000,000 human UniGenes (NCBI). Clearly, these clusters of cDNAs greatly overestimate the number of unique genes. Similarly, FunnyBase has multiple clusters for the same gene: 15 apolipoprotein I, 10 cytochrome oxidase I, and 53 vitellogenin clusters. To provide a more precise estimate of the number of unique genes, consensus sequences were queried against the 27,695 sequences from the Human RefSeq  database, then grouped by identical gene symbol. Of the 2,376 Fundulus clusters that were similar to a sequence in Human RefSeq (e-value < 10-10), 1,818 (76%) had distinct gene annotations. This method of clustering clusters by similarity to well-annotated reference sequences provides a method to more accurately define the number of unique genes represented by an EST set.
Case study: using functional annotation for microarray analysis
The web interface http://genomics.rsmas.miami.edu/funnybase/super_craw4/ provides public access to the FunnyBase system and dataset. Searches can query by keyword in annotation, gene name, GO term, metabolic pathway, clone or plate id, and BLAST homology search. All data including raw sequences, cluster memberships, cluster alignments, and alignments with homologous sequences are provided for the user to examine the source of annotations. Links associated with each annotation are made to external resources such as GO's AMIGO browser, KEGG pathways, SwissProt, and NCBI records.
Other Fundulus sequences
FunnyBase was constructed to annotate sequences for the analysis of gene expression. It provides identification and annotation for genes in the Crawford laboratory with a primary goal of identifying clones useful for the construction of microarrays. As such, other Fundulus sequences in Genbank are not included. However, FunnyBase forms the basis of the TIGR Killifish gene index http://www.tigr.org/tigr-scripts/tgi/T_index.cgi?species=killifish that includes publicly available F. heteroclitus sequences.
Customized species specific EST databases are available for many species [17, 18, 21, 49–56]. FunnyBase provides an integrated method to annotate ESTs with the most biologically relevant set of associations and provides several innovations for the production of ESTs for microarrays. Control sequences are identified in each 96-well plate so that mislabelled or inverted plates are automatically detected. Annotations are based upon several different public databases. The multiple annotations provide greater assurance about gene description and greater frequency of annotation than any one database. The most functionally informative innovation of FunnyBase is the process of culling through numerous primary similarity search results in order to identify links to systematic functional databases in GO and KEGG. These provide a discrete set of terms that can be analyzed statistically and that are organized into networks that represent biological knowledge of higher-level functional and pathway associations. The range of databases queried by similarity search and the tracking of homology beyond a single 'best' hit maximizes the opportunity to obtain this annotation. A richer set of GO terms is achieved by using all hits with e-values less than 10-11 that represent the same gene as the 'best hit'. Additional GO terms that represent more general functions than those found in public annotation are derived through the parent-child relationship of the Gene Ontology. EC numbers provide links, via KEGG, to metabolic pathways and these stored terms can be used to investigate the relationship between gene expression in specific metabolic pathways including cardiac metabolism . To provide a more accurate accounting of the number of unique genes, consensus sequence from clusters of ESTs were queried against the Human RefSeq database and those sequences sharing the same gene symbol are grouped based on this scaffolding information. These approaches use publicly available bioinformatics tools (BLAST, CAP3, Phred, Cross-Match, Perl, and the MySQL database management system). The application of theses tools in an appropriate framework as outlined in FunnyBase can be used to create a systems level functional genomics annotation system useful for EST databases to study biological processes among a rich diversity of organisms.
The animal protocols used in the present study have been approved by the University of Miami Institute Animal Care and Use Committee. The teleost fish Fundulus heteroclitus used for ESTs were collected from two sites: Scorton Creek in Sandwich, MA, and Stone Harbor, NJ. These populations are in the central portion of the thermal cline and have relatively high levels of heterozygosity . These fish were subjected to the following environmental regime before tissues were harvested for mRNA extraction: kept in controlled temperature and aeration conditions, and acclimated to common conditions (20°C, 15 ppt salinity) in re-circulating aquaria for at least nine months before experiments. Following this common acclimation a subset of fish were subjected to one of several stresses: 4°C, 34°C, hypoxia, or a complex mix of hydrocarbons.
To effectively isolate and sequence thousands of cDNAs for the production of microarrays, a unidirectional cDNA library with few non-recombinants was required. We created four cDNA libraries: heart libraries from non-stressed and stressed fish and liver libraries from non-stressed and stressed fish. The non-stressed F. heteroclitus cardiac and liver libraries were provided by Drs. S. Karchner and M. Hahn, WHOI  and were constructed using the UniZap λ cDNA Gigpack Gold cloning kit (Stratagene, La Jolla, CA, USA). The cardiac library was produced from 27 fish hearts (both sexes) sampled from Scorton Creek in Sandwich, MA. The cDNAs in these libraries are oriented such that the 5'end of each cDNA is ligated to EcoR1 and 3' poly A is ligated to XhoI. These libraries had less than 1% non-recombinants, i.e. 2 of 300 random clones from a non-normalized library had no inserts. The stressed libraries included 4 fish subjected to the four stressors (above) and 4 non-stressed individuals. Unidirectional heart and liver libraries were constructed such that the 5'end of each cDNA is ligated to EcoR1 and 3' poly A is ligated to XhoI of the plasmid vector pSmart (Lucigen, Middleton, WI, USA). The pSmart-cDNA vector was designed for EST work. The vector expresses kanamycin-resistance and has a terminator on both sides of the cDNA insertion site preventing expression of cDNA. These two attributes (non-expression and Kan-resistance) increase the stability of different genes in the library versus cDNA libraries in Amp libraries with Lac promoters (Crawford, unpublished). These libraries had less than 1% non-recombinants.
Normalization of cDNA libraries reduces the differences among expressed genes to less than 10-fold among rare and abundant mRNAs [44, 45]. Normalized libraries were produced by isolating cDNAs from approximately 1012 plasmids. The cDNAs were isolated using PCR amplification with vector specific primers immediately 5' and 3' to the insertion site (EcoRI and XhoI sites). These PCR products (PCR-cDNAs) were denatured and hybridized to single stranded plasmids from the cardiac cDNA library. Taking advantage of Cot values, the most abundant cDNAs were annealed to the more abundant PCR products and were removed selectively by hydroxyapatite-column chromatography. The single-stranded plasmids in the flow-through were converted to double strands using the Sequenase DNA polymerase (Amersham, Piscataway, NJ, USA). DH10s E. coli (BRL) were transformed with these double-stranded plasmids by electroporation. The number of recovered plasmids and the resulting complexity of the normalized library depended on the duration of hybridization or Cot values. Two normalized libraries were made using either a 12 or 24 hour hybridization. The library from the 12-hour hybridization yielded 250,000 plasmids. The library from the 24 hour hybridization yielded 3,000 plasmids and had a greater representation of rare mRNAs and greater frequency of non-recombinants.
Isolation and sequencing of cDNAs
Characterization of cDNAs (growth of individual bacterial colonies containing plasmids, PCRs, purification of PCR products, sequencing reactions) used 96 well plates and octopipettes. To characterize cDNAs, 96 individual bacterial colonies from the normalized library were randomly chosen, and each was grown in 1.25 ml of Superbroth in 2 ml-96 well plates. After 18 hours of growth, two 250 ul bacterial glycerol stocks were made and stored in 96 well plates at -80°C. One microliter of these bacterial growths was used for PCR reactions using forward and reverse plasmid specific primers: (PucF = CGCCAGGGTTTTCCCAGTCACG, PucR = GAGCGGATAACAATTTCACACAGGAAA). PCR reactions had 0.2 mM dNTPs, 10 pmoles of each primer, 1 unit of Promega Taq (0.2 ul), and reaction buffer with detergents and DMSO (final concentrations: 50 mM Tris HCl, pH 9.2 (25°C), 16 mM (NH4)2SO4, 2.25 mM MgCl2, 2% (v/v) DMSO, 0.1% (v/v) Tween 20). Two-step thermal cycle conditions were used (94°C for 10 seconds; then 32 cycles of 94°C for 30 seconds followed by 70°C for 5 minutes; then 72°C for 15 minutes). PCR products were purified manually in 96 well format using Sephadex G-50 in a deep well plate with a 0.2 microfilter (Millipore, Billerica, USA) or robotically using AmPure (Agencourt, Beverly, MA, USA) and EvolP3 96 pipetting liquid handling system (PerkinElmer Life Sciences Inc., Boston, MA, USA).
PCR products were sequenced from the 5' end (relative to the mRNA) on an ABI 373 or ABI 3730 sequencer using ABI "Big Dye" reaction mix. We typically used 1/16 the amount of reaction mix, yielding 300 to 400 unambiguous bases. Sequences were purified using biotin primers and streptavidin coated magnetic beads (for the ABI 373) or Agencort CleanSeq (for the ABI 3730).
We used three procedures to verify that the correct sequence was associated with each cDNA. 1) Each 96-well plate had three wells with a "marker cDNA" (Ctenophore cDNA #5, a random cDNA with no similarity to any sequence in GenBank). Two wells (#40 and #67) always contained the marker cDNA, and thus any misloading or mislabeling of sequencing lanes was identifiable. The third marker cDNA was placed in a well that corresponds to the plate number (e.g., plate 2 had the marker in well 2). 2) After the production of 12 plates, one row (8 wells) from each plate was re-sequenced. Thus, 8/96 or ~8% of all sequences and their locations were confirmed. 3) cDNAs used for microarrays were re-sequenced. These measures are important to ensure that the correct and known cDNAs are printed.
The complexity of the normalized library was reduced by subtracting the characterized cDNAs previously isolated from the normalized library. Subtraction greatly reduced the probability of isolating the same cDNA and thus improved the efficiency of screening the library for unique clones. Subtraction used a 100-fold molar excess of biotin-labeled antisense cDNAs produced by PCR using all the characterized cDNAs as substrates and vector-specific primers in which the 3' primer was labeled with biotin. These PCR products were hybridized to the cDNA libraries in the presence of oligo-dA and vector-specific oligos (that prevented non-specific hybridization to oligo-dT or vector sequences). After a 24 hour hybridization, genes in the library that bound to these biotin-labeled PCR products were removed with the use of magnetized, streptavidin coated beads. DH10s E. coli were transformed with the subtracted library by electroporation.
Hardware and software
Software and Databases. Publicly available software and databases used for FunnyBase. The version and/or download date are listed.
Version and/or Download Date
Public Sequence Similarity Databases
Rainbow Trout Unigene
Microarrays were printed using a select 384 cDNAs from F. heteroclitus cardiac library encoding essential proteins for cellular metabolism isolated from over 40,000 expressed sequences http://genomics.rsmas.miami.edu/FunnyBase/super_craw4/. These 384 cDNAs were amplified with amine-linked primers and printed on 3-D Link Activated slides (Surmodics Inc., Eden Prairie, MN, USA) using GeneMachine OminGrider, and blocked following slide manufacturer protocols. The suite of 384 amplified cDNAs was printed as a group in four spatially separated replicates. Four hybridization zones of these four replicate arrays were printed per slide, with each zone set separated by a hydrophobic barrier. Samples were hybridized twice; once with Cy3 and once with Cy5 resulting in overall technical replication of 8-fold per sample.
Sample preparation and hybridization
RNA was extracted from tissue homogenate in a chaotropic buffer using phenol/cholorform/isoamyl alcohol and RNA quality was assessed using the Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA). RNA for hybridization was prepared by amplification using a modified Eberwine protocol  using the Ambion Amino Allyl MessageAmp aRNA Kit. Cy3 and Cy5 were hybridized to slides and incubated 12–18 hours at 42°C. Following hybridization, slides were scanned using the Packard Bioscience ScanArray Express microarray scanner (PerkinElmer Life Sciences Inc., Boston, MA, USA) and images processed using ImaGene (Biodiscovery Inc., Marina del Rey, CA, USA).
Initial support for the F. heteroclitus EST project was provided by NSF BioInformatics post-doctoral fellowship 0074520 to MFO and NSF/IBN grant 9986602 to DLC. We would like to thank Kristin Horgan of M.J. research for arranging the loan of Tetrad-thermal cycler for the Compartive Functional Genomic Course. Additionally, we would like to thank AP-Biotech and specifically Dr. Robert Feldman for use of MegaBace used to re-sequence cDNAs.
Current EST isolation, sequencing and bioinformatics is supported by NSF/OCE grant 0221879 and NIH/NHLBI R01 HL65470 to DLC and NIH/NIEHS ES011588 to MFO.
- Sorlie T, Perou CM, Tibshirani R, Aas T, Geisler S, Johnsen H, Hastie T, Eisen MB, van RM, Jeffrey SS, Thorsen T, Quist H, Matese JC, Brown PO, Botstein D, Lonning PE, Borresen-Dale AL: Gene expression patterns of breast carcinomas distinguish tumor subclasses with clinical implications. Proceedings of the National Academy of Sciences of the United States of America. 2001, 98: 10869-10874. 10.1073/pnas.191367098.PubMed CentralView ArticlePubMedGoogle Scholar
- Ross D, Scherf U, Eisen M, Perou C, Rees C, Spellman P, Iyer V, Jeffrey S, Van RM, Waltham M, Pergamenschikov A, Lee J, Lashkari D, Shalon D, Myers T, Weinstein J, Botstein D, Brown P: Systematic variation in gene expression patterns in human cancer cell lines. Nature Genetics. 2000, 24: 227-235. 10.1038/73432.View ArticlePubMedGoogle Scholar
- Friddle CJ, Koga T, Rubin EM, Bristow J: Expression profiling reveals distinct sets of genes altered during induction and regression of cardiac hypertrophy. PNAS. 2000, 97: 6745-6750. 10.1073/pnas.100127897.PubMed CentralView ArticlePubMedGoogle Scholar
- Archacki SR, Angheloiu G, Tian XL, Tan FL, DiPaola N, Shen GQ, Moravec C, Ellis S, Topol EJ, Wang Q: Identification of new genes differentially expressed in coronary artery disease by expression profiling. Physiol Genomics. 2003, 15: 65-74.View ArticlePubMedGoogle Scholar
- Iemitsu M, Miyauchi T, Maeda S, Sakai S, Fujii N, Miyazaki H, Kakinuma Y, Matsuda M, Yamaguchi I: Cardiac hypertrophy by hypertension and exercise training exhibits different gene expression of enzymes in energy metabolism. Hypertens Res. 2003, 26: 829-837. 10.1291/hypres.26.829.View ArticlePubMedGoogle Scholar
- Steenman M, Chen YW, Le Cunff M, Lamirault G, Varro A, Hoffman E, Leger JJ: Transcriptomal analysis of failing and nonfailing human hearts. Physiol Genomics. 2003, 12: 97-112.View ArticlePubMedGoogle Scholar
- Larkin JE, Frank BC, Gaspard RM, Duka I, Gavras H, Quackenbush J: Cardiac transcriptional response to acute and chronic angiotensin II treatments. Physiol Genomics. 2004, 152-66. 10.1152/physiolgenomics.00057.2004. 18
- Bono H, Nikaido I, Kasukawa T, Hayashizaki Y, Okazaki Y: Comprehensive analysis of the mouse metabolome based on the transcriptome. Genome Res. 2003, 13: 1345-1349. 10.1101/gr.974603.PubMed CentralView ArticlePubMedGoogle Scholar
- Lucas DT, Aryal P, Szweda LI, Koch WJ, Leinwand LA: Alterations in mitochondrial function in a mouse model of hypertrophic cardiomyopathy. Am J Physiol Heart Circ Physiol. 2003, 284: H575-83.View ArticlePubMedGoogle Scholar
- McCammon MT, Epstein CB, Przybyla-Zawislak B, McAlister-Henn L, Butow RA: Global transcription analysis of Krebs tricarboxylic acid cycle mutants reveals an alternating pattern of gene expression and effects on hypoxic and oxidative genes. Mol Biol Cell. 2003, 14: 958-972. 10.1091/mbc.E02-07-0422.PubMed CentralView ArticlePubMedGoogle Scholar
- Toye A, Gauguier D: Genetics and functional genomics of type 2 diabetes mellitus. Genome Biology. 2003, 4: 241-10.1186/gb-2003-4-12-241.PubMed CentralView ArticlePubMedGoogle Scholar
- Yechoor VK, Patti ME, Saccone R, Kahn CR: Coordinated patterns of gene expression for substrate and energy metabolism in skeletal muscle of diabetic mice. Proc Natl Acad Sci U S A. 2002, 99: 10587-10592. 10.1073/pnas.142301999.PubMed CentralView ArticlePubMedGoogle Scholar
- Grabher C, Henrich T, Sasado T, Arenz A, Wittbrodt J, Furutani-Seiki M: Transposon-mediated enhancer trapping in medaka. Gene. 2003, 322: 57-66. 10.1016/j.gene.2003.09.009.View ArticlePubMedGoogle Scholar
- Inoue S, Nam BH, Hirono I, Aoki T: A survey of expressed genes in Japanese flounder (Paralichthys olivaceus) liver and spleen. Molecular Marine Biology & Biotechnology. 1997, 6: 376-380.Google Scholar
- Lo J, Lee S, Xu M, Liu F, Ruan H, Eun A, He Y, Ma W, Wang W, Wen Z, Peng J: 15,000 Unique zebrafish EST clusters and their future use in microarray for profiling gene expression patterns during embryogenesis. Genome Research. 2003, 13: 455-466. 10.1101/gr.885403.PubMed CentralView ArticlePubMedGoogle Scholar
- Savan R, Sakai M: Analysis of expressed sequence tags (EST) obtained from common carp, Cyprinus carpio L., head kidney cells after stimulation by two mitogens, lipopolysaccharide and concanavalin-A. Comparative Biochemistry & Physiology Part B, Biochemistry & Molecular Biology 131B. 2002, 1: 71-82. 10.1016/S1096-4959(01)00488-2.View ArticleGoogle Scholar
- Mita K, Morimyo M, Okano K, Koike Y, Nohata J, Kawasaki H, Kadono-Okuda K, Yamamoto K, Suzuki MG, Shimada T, Goldsmith MR, Maeda S: The construction of an EST database for Bombyx mori and its application. Proceedings of the National Academy of Sciences of the United States of America. 2003, 100: 14121-14126. 10.1073/pnas.2234984100.PubMed CentralView ArticlePubMedGoogle Scholar
- Rise ML, von Schalburg KR, Brown GD, Mawer MA, Devlin RH, Kuipers N, Busby M, Beetz-Sargent M, Alberto R, Gibbs AR, Hunt P, Shukin R, Zeznik JA, Nelson C, Jones SR, Smailus DE, Jones SJ, Schein JE, Marra MA, Butterfield YS, Stott JM, Ng SH, Davidson WS, Koop BF: Development and application of a salmonid EST database and cDNA microarray: data mining and interspecific hybridization characteristics. Genome Res. 2004, 14: 478-490. 10.1101/gr.1687304.PubMed CentralView ArticlePubMedGoogle Scholar
- Whitfield CW, Band MR, Bonaldo MF, Kumar CG, Liu L, Pardinas JR, Robertson HM, Soares MB, Robinson GE: Annotated Expressed Sequence Tags and cDNA Microarrays for Studies of Brain and Behavior in the Honey Bee. Genome Res. 2002, 12: 555-566. 10.1101/gr.5302.PubMed CentralView ArticlePubMedGoogle Scholar
- Oleksiak MF, Kolell K, Crawford DL: The utility of natural populations for microarray analyses: isolation of genes necessary for functional genomic studies. Marine Biotechnology. 2001, 3: S203-S211. 10.1007/s10126-001-0043-0.View ArticlePubMedGoogle Scholar
- Ryan JF, Finnerty JR: CnidBase: The Cnidarian Evolutionary Genomics Database. Nucleic Acids Research. 2003, 31: 159-163. 10.1093/nar/gkg116.PubMed CentralView ArticlePubMedGoogle Scholar
- Crawford DL: Functional genomics does not have to be limited to a few select organisms. Genome Biology. 2001, 2: INTERACTIONS1001 http://www.genomebiology.com/2001/2/1/interactions/1001/-10.1186/gb-2001-2-1-interactions1001.View ArticleGoogle Scholar
- Oleksiak MF, Churchill GA, Crawford DL: Variation in gene expression within and among natural populations. Nature Genetics. 2002, 32: 261-266. 10.1038/ng983.View ArticlePubMedGoogle Scholar
- Whitfield CW, Cziko AM, Robinson GE: Gene Expression Profiles in the Brain Predict Behavior in Individual Honey Bees. Science. 2003, 302: 296-299. 10.1126/science.1086807.View ArticlePubMedGoogle Scholar
- Evans JD, Wheeler DE: Expression profiles during honeybee caste determination. Genome Biology. 2001, 2: RESEARCH0001-PubMed CentralPubMedGoogle Scholar
- Gracey AY, Troll JV, Somero GN: Hypoxia-induced gene expression profiling in the euryoxic fish Gillichtys mirabilis. Proceeding of National Academy of Science, USA. 2001, 98: 1993-1998. 10.1073/pnas.98.4.1993.View ArticleGoogle Scholar
- Podrabsky JE, Somero GN: Changes in gene expression associated with acclimation to constant temperatures and fluctuating daily temperatures in an annual killifish Austrofundulus limnaeus. J Exp Biol. 2004, 207: 2237-2254. 10.1242/jeb.01016.View ArticlePubMedGoogle Scholar
- Ju Z, Dunham RA, Liu Z: Differential gene expression in the brain of channel catfish (Ictalurus punctatus) in response to cold acclimation. Molecular Genetics and Genomics. 2002, 268: 87-95. 10.1007/s00438-002-0727-9.View ArticlePubMedGoogle Scholar
- Pierce VA, Crawford DL: Phylogenetic analysis of glycolytic enzyme expression. Science. 1997, 275: 256-259. 10.1126/science.276.5310.256.View ArticleGoogle Scholar
- Brown BL, Chapman RW: Gene flow and mitochondrial DNA variation in the killifish Fundulus heteroclitus. Evolution. 1991, 45: 1147-1161.View ArticleGoogle Scholar
- Crawford DL, Powers DA: Molecular basis of evolutionary adaptation at the lactate dehydrogenase-B locus in the fish Fundulus heteroclitus. Proceedings of the National Academy of Sciences of the United States of America. 1989, 86: 9365-9369.PubMed CentralView ArticlePubMedGoogle Scholar
- Powers DA, Smith M, Gonzalez-Villasenor I, DiMichelle L, Crawford DL, Bernardi G, Lauerman T: A multidisciplinary approach to the selectionist/neutralist controversy using the model teleost, Fundulus heteroclitus. Oxford Surveys in Evolutionary Biology. Edited by: Futuyma D and Antonovics J. 1993, New York, NY, Oxford University Press, 9: 43-108.Google Scholar
- Pierce VA, Crawford DL: Phylogenetic analysis of thermal acclimation of the glycolytic enzymes in the genus Fundulus. Physiological Zoology. 1997, 70: 597-609.View ArticlePubMedGoogle Scholar
- Pierce VA, Crawford DL: Phylogenetic analysis of glycolytic enzyme expression. Science. 1997, 276: 256-259. 10.1126/science.276.5310.256.View ArticlePubMedGoogle Scholar
- TIGR: http://www.tigr.org/tdb/tgi/.
- UniGene: http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=unigene.
- dbESTs NCBI: . (http://wwwncbinlmnihgov/dbEST/dbEST_summaryhtml). 1 July 2004:
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nature Genetics. 2000, 25: 25-29. 10.1038/75556.PubMed CentralView ArticlePubMedGoogle Scholar
- Wixon J, Kell D: The Kyoto encyclopedia of genes and genomes--KEGG. Yeast. 2000, 17: 48-55. 10.1002/(SICI)1097-0061(200004)17:1<48::AID-YEA2>3.0.CO;2-H.View ArticlePubMedGoogle Scholar
- Apweiler R, Bairoch A, Wu CH, Barker WC, Boeckmann B, Ferro S, Gasteiger E, Huang H, Lopez R, Magrane M, Martin MJ, Natale DA, O'Donovan C, Redaschi N, Yeh LS: UniProt: the Universal Protein knowledgebase. Nucleic Acids Research. 2004, 32 Database issue: D115-9. 10.1093/nar/gkh131.View ArticleGoogle Scholar
- NCBI: NCBI NR (non-redudant) Database.
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. Journal of Molecular Biology. 1990, 215: 403-410. 10.1006/jmbi.1990.9999.View ArticlePubMedGoogle Scholar
- Huang X, Madan A: CAP3: A DNA Sequence Assembly Program. Genome Research. 1999, 9: 868-877. 10.1101/gr.9.9.868.PubMed CentralView ArticlePubMedGoogle Scholar
- Hillier LD, Lennon G, Becker M, Bonaldo MF, Chiapelli B, Chissoe S, Dietrich N, DuBuque T, Favello A, Gish W, Hawkins M, Hultman M, Kucaba T, Lacy M, Le M, Le N, Mardis E, Moore B, Morris M, Parsons J, Prange C, Rifkin L, Rohlfing T, Schellenberg K, Marra M: Generation and analysis of 280,000 human expressed sequence tags. Genome Research. 1996, 6: 807-828.View ArticlePubMedGoogle Scholar
- Bonaldo MF, Lennon G, Soares MB: Normalization and subtraction: two approaches to facilitate gene discovery. Genome Research. 1996, 6: 791-806.View ArticlePubMedGoogle Scholar
- McKusick-Nathans Institute for Genetic Medicine JHUBMDNCBINLM: Online Mendelian Inheritance in Man, OMIM (TM). [http://www.ncbi.nlm.nih.gov/omim/]
- Whitehead JA, Crawford DL: Variation in tissue-specific gene expression. Genome Biology. 2004, In Press:Google Scholar
- Pruitt KD, Maglott DR: RefSeq and LocusLink: NCBI gene-centered resources. Nucleic Acids Research. 2001, 29: 137-140. 10.1093/nar/29.1.137.PubMed CentralView ArticlePubMedGoogle Scholar
- Abdrakhmanov I, Lodygin D, Geroth P, Arakawa H, Law A, Plachy J, Korn B, Buerstedde JM: A large database of chicken bursal ESTs as a resource for the analysis of vertebrate gene function. Genome Res. 2000, 10: 2062-2069. 10.1101/gr.10.12.2062.PubMed CentralView ArticlePubMedGoogle Scholar
- Alkharouf NW, Matthews BF: SGMD: the Soybean Genomics and Microarray Database. Nucleic Acids Res. 2004, 32 Database issue: D398-400. 10.1093/nar/gkh126.View ArticleGoogle Scholar
- Dong Q, Roy L, Freeling M, Walbot V, Brendel V: ZmDB, an integrated database for maize genome research. Nucleic Acids Res. 2003, 31: 244-247. 10.1093/nar/gkg082.PubMed CentralView ArticlePubMedGoogle Scholar
- Douglas SE, Gallant JW, Bullerwell CE, Wolff C, Munholland J, Reith ME: Winter Flounder Expressed Sequence Tags: Establishment of an EST Database and Identification of Novel Fish Genes. 1999, 1: 458-464.Google Scholar
- Gai X, Lal S, Xing L, Brendel V, Walbot V: Gene discovery using the maize genome database ZmDB. Nucleic Acids Res. 2000, 28: 94-96. 10.1093/nar/28.1.94.PubMed CentralView ArticlePubMedGoogle Scholar
- Jung S, Jesudurai C, Staton M, Du Z, Ficklin S, Cho I, Abbott A, Tomkins J, Main D: GDR (Genome Database for Rosaceae): integrated web resources for Rosaceae genomics and genetics research. BMC Bioinformatics. 2004, 5: 130-10.1186/1471-2105-5-130.PubMed CentralView ArticlePubMedGoogle Scholar
- Nobis W, Ren X, Suchyta SP, Suchyta TR, Zanella AJ, Coussens PM: Development of a porcine brain cDNA library, EST database, and microarray resource. Physiol Genomics. 2003, 16: 153-159. 10.1152/physiolgenomics.00099.2003.View ArticlePubMedGoogle Scholar
- Tomkins JP, Luo M, Fang GC, Main D, Goicoechea JL, Atkins M, Frisch DA, Page RE, Guzman-Novoa E, Yu Y, Hunt G, Wing RA: New genomic resources for the honey bee(Apis mellifera L.): development of a deep-coverage BAC library and a preliminary STC database. Genet Mol Res. 2002, 1: 306-316.PubMedGoogle Scholar
- Oleksiak MF, Roach JL, Crawford DL: Natural Variation in Cardiac Metabolism and Gene Expression in Fundulus heteroclitus. Nature Genetics. 2004, Advance_Online 28 November 2004; doi:10.1038/ng1483:Google Scholar
- Karchner SI, Powell WH, Hahn ME: Identification and functional characterization of two highly divergent aryl hydrocarbon receptors (AHR1 and AHR2) in the teleost Fundulus heteroclitus. Evidence for a novel subfamily of ligand-binding basic helix loop helix-Per-ARNT-Sim (bHLH-PAS) factors. Journal of Biological Chemistry. 1999, 274: 33814-33824. 10.1074/jbc.274.47.33814.View ArticlePubMedGoogle Scholar
- Van Gelder RN, Von Zastrow ME, Yool A, Dement WC, Barchas JD, Eberwine JH: Amplified Rna Synthesized from Limited Quantities of Heterogeneous Complementary DNA. Proceedings of the National Academy of Sciences of the United States of America. 1990, 87: 1663-1667.PubMed CentralView ArticlePubMedGoogle Scholar