1471-2164-5-96 1471-2164 Research article <p>FunnyBase: a systems level functional annotation of <it>Fundulus </it>ESTs for the analysis of gene expression</p> Paschall E Justin jep67c@umkc.edu Oleksiak F Marjorie mfoleksi@ncsu.edu VanWye D Jeffrey jvanwye@rsmas.miami.edu Roach L Jennifer jlroach@rsmas.miami.edu Whitehead J Andrew awhitehead@rsmas.miami.edu Wyckoff J Gerald wyckoffg@umkc.edu Kolell J Kevin kolellk@umkc.edu Crawford L Douglas dcrawford@rsmas.miami.edu

Division of Molecular Biology and Biochemistry, 5100 Rockhill Rd., University of Missouri-Kansas City 64110, USA

Department of Environmental & Molecular Toxicology, North Carolina State University; Raleigh, NC 27695-7633 USA

Division of Marine Biology and Fisheries, NIEHS Marine and Freshwater Biomedical Sciences Center, Rosenstiel School of Marine & Atmospheric Science, University of Miami, Miami, FL 33149, USA

BMC Genomics 1471-2164 2004 5 1 96 http://www.biomedcentral.com/1471-2164/5/96 15610557 10.1186/1471-2164-5-96
17 8 2004 20 12 2004 20 12 2004 2004 Paschall et al; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

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.

Results

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.

Conclusion

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 jpaschall@bioinfobase.umkc.edu.

Background

Investigating patterns of gene expression using mouse and human microarrays has produced insights into cancer 12, cardiac diseases 3456, and metabolic disorders 789101112. 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 1314151617181920 and including simple diploblastic organisms 21. These studies within a diversity of organisms provide insights not provided by 'model' species (species that are genetically well defined or with annotated genomes 22). For example, 'non-model' organisms have provided insight into the natural variation in gene expression 23, social castes among bees 2425, hypoxia 26, and physiological responses to variation in the thermal environment 2728. To investigate adaptive variation in gene expression we use the teleost Fundulus heteroclitus (killifish) 2329.

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 30 but not frequent enough to extinguish local adaptation 3132. 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 233334. 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 35 and Unigene 36 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 37. 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) 38 assignments and membership in metabolic pathways as defined by the Kyoto Encyclopedia of Genes and Genomes (KEGG) 39. 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

FunnyBase (Fig. 1) is divided into 3 modules: Sequence Pipeline, Hierarchical Annotation, and Microarray Production and Analysis. The Sequence Pipeline takes sequences and quality output files from the sequencer, applies vector screening, quality trimming, clone tracking, and clustering (described below) to produce a set of unique sequences that are deposited in the 'Sequence Data' and 'Cluster Data' tables. Notice "unique sequences" are a combination of singletons (single unique ESTs) and clusters of overlapping sequences.

<p>Figure 1</p>

FunnyBase annotation scheme

FunnyBase annotation scheme. The integration of the three FunnyBase modules: sequence pipeline, hierarchical annotations and microarray production and analysis. Database tables are shown in cylinders, arrows are data flow, and dashed lines indicate the integration of data from multiple sources.

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 40, NCBI UniGene 36, and NCBI non-redundant NR 41) 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.

Controls

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 42 against previous sequencing results so that the identity of printed microarray spots are verified.

EST clustering

EST projects generate a number of redundant sequences due to the random selection of cDNAs from tissue libraries (Table 1). Clustering redundant sequences is a critical first step of analysis in order to identify genes to target for subtraction. The program CAP3 by Xiaong Huang 43 was used to cluster EST sequences with a 30 bp overlap and 75 percent similarity.

<p>Table 1</p>

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.

ID

Number of ESTs

Evalue

Description


1616

977

0

Vitellogenin I precursor

2262

734

2e-88

Cytochrome c oxidase polypeptide II

1348

507

1e-102

Alpha-1-antitrypsin homolog precursor

1026

481

5e-61

Zona pellucida sperm-binding protein 3 precursor

1727

401

0

Serotransferrin precursor

555

397

0

Cytochrome c oxidase polypeptide I

640

369

8e-61

Zona pellucida sperm-binding protein B precursor

1549

351

2e-77

Apolipoprotein A-I precursor

570

331

1e-111

Cytochrome c oxidase polypeptide III

1178

304

2e-45

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.

The distribution of the number of clusters with two or more ESTs is depicted in Figure 2. Although most clusters (2,581 or 68%) in FunnyBase have two or three-to-four sequences (Fig. 2), a small number of highly expressed genes form clusters with a large number of ESTs. For example, there are three clusters that contain 512 to 1,024 ESTs and ten clusters that contain 256 to 512 ESTs (bottom two bars for killifish in Fig. 2). This distribution is similar to other teleost fish EST collections (Fig 2). Notice, as more ESTs are added, clusters tend to get larger (more ESTs per cluster) rather than new small clusters growing in frequency. Of the 3,779 killifish clusters, 14% have more than eight ESTs, yet of the 14,714 Medaka clusters, 48% have eight or more ESTs. These distributions suggest that adding more EST sequences has diminishing returns.

<p>Figure 2</p>

Frequency of cluster size class in teleosts

Frequency of cluster size class in teleosts. The frequency of the number of ESTs in each cluster is shown for Fundulus heteroclitus (Killifish: total number of clusters 3,779), Oryzias latipes (Medaka: total number of clusters 7,401), Oncorhynchus mykiss (Rainbow Trout: total number of clusters 11,405) and Danio rerio (Zebrafish: total number of clusters 14,714). For example, the first bar indicates that approximately 1,000 clusters contain 2 ESTs in all four teleost fish. Clusters for other species (not F. heteroclitus) are based on NCBI UniGene.

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 4445. 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.

Gene annotation

Of the 12,776 unique ESTs (3,776 clusters and 8,991 singletons), 3,877 (30%) were annotated. The distribution of e-values for these annotations is shown in Figure 3. Most (84%) of these ESTs have e-values less than 10-5. Among the clusters, 2,265 of 3, 779 (60%) were annotated as compared to 333 of 1,131 (30%) of confirmed high quality singletons. The lower percent of annotated singletons suggests that these are either rare fish-specific genes, or represent otherwise divergent, likely non-coding (5' or 3' UTR), regions.

<p>Figure 3</p>

Distribution of e-values for annotations

Distribution of e-values for annotations. Gene annotation is based on sequence similarity. The e-values, which describe the probability of random sequence similarity, are shown as the negative log value (e.g., 10-5 = 5).

Annotations for ESTs are based on similarity using BLASTX or BLASTN 42 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.

The use of multiple databases increases the total number of annotated ESTs (Fig. 4) as compared with any one source and provides opportunity to compare annotation between all three sources for 1,841 (47%) sequences. GenBank NR provided the most number of annotations, but these tend to be less informative (see systematic functional annotations below). Human Unigene provided an additional 311 (8% of total) annotations. SwissProt provided an additional 32 annotations (1% of total) with 743 fewer annotated sequences than the NR. However, SwissProt is uniformly well annotated as compared to NR where informative functional annotation can easily be buried by numerous uninformative hits at similar e-values. Besides increasing the number of annotations, comparing the annotations from multiple databases ensures that mistakes in the curation can be detected and information such as alternative gene names can be compiled from multiple sources.

<p>Figure 4</p>

Venn diagram of annotations from three different databases

Venn diagram of annotations from three different databases. Number of unique ESTs annotated by three different databases: GenBank non-redundant (NR, total annotated = 3,534), Swiss-Prot (total annotated = 2,829), and human Unigene (total annotated = 2,390).

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 46 that are not available in the single FASTA description line returned by BLAST search.

KEGG 39 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.

Of the 3,877 annotated ESTs in FunnyBase, 588 (14%) participate in one or more pathways defined by KEGG. These 588 ESTs represent 105 different pathways. Table 2 provides a breakdown of the number of ESTs in FunnyBase for the 10 pathways associated with the largest number of distinct sequences (contigs or singletons). The extent that a given pathway is represented in FunnyBase can be used to identify metabolic differences among tissues 47 or in different species.

<p>Table 2</p>

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


Glycolysis/Gluconeogenesis

89

Oxidative phosphorylation

86

Fatty acid metabolism

70

Pyruvate metabolism

69

Tryptophan metabolism

66

Butanoate metabolism

48

Glycerolipid metabolism

47

Valine, leucine and isoleucine degradation

46

Glycine, serine and threonine metabolism

44

Propanoate metabolism

44

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 38. 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 48 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

As a case study in how functional annotation in FunnyBase can be integrated with microarray data in a rigorous manner, we used a data set based on a microarray of metabolic genes printed from ESTs annotated in FunnyBase 47. Statistical analysis of this set of 363 metabolic genes identified a set of 62 genes that showed statistically significant regression between gene expression levels and temperature along the Atlantic coast. That is, among individuals collected from different locations along the thermocline and then acclimated to common physiological conditions for at least nine months before analysis, 17% of the metabolic genes had a linear relationship between the amount of mRNA and the environmental temperature these animals evolved in. Our hypothesis was that this set of 59 genes represents a functionally different set than those genes that do not show regression with temperature. To test this hypothesis we examined the frequency of genes annotated with a given GO term in the statistically significant gene set versus the non-significant genes. Figure 5 shows the relevant proportions in each set for GO terms that are represented by 5 or more ESTs in the significant set. For example, the GO term Amine Metabolism (GO:0009308) is assigned to 14% of the 62 statistically significant genes but only 3% of the non-significant genes (those that do not show significant regression with temperature). A Fisher-exact test indicates these frequencies (14% vs. 3%) represent different underlying distributions (p < 0.001). Specifically, genes involved in amine metabolism are overrepresented in the set of genes that show regression with temperature as compared with the remaining sequences. This significant increase is found for two other non-mutually exclusive GO terms: amino acid and derivative metabolism, and amino acid metabolism (p < 0.05). Other GO terms show a reverse trend although none were statistically significant. For example, ion transport (p = 0.08) and cell growth (p = 0.16) had few genes with a clinal variation in expression. These data suggest that the functions of genes influence whether they are affected by ecologically interesting patterns of expression (Fig 5).

<p>Figure 5</p>

Distribution of significant and non-significant genes relative to GO terms

Distribution of significant and non-significant genes relative to GO terms. The relationship between gene expression from "common gardened" fish and the environment they evolved in was statistically analyzed and grouped by GO terms. Black bars represent genes whose levels of expression has a significant regression with the environmental thermal cline among populations (p < 0.05). Hatched bars represent the set of genes with no significant relationship to the thermal environment.

Web interface

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.

Conclusions

Customized species specific EST databases are available for many species 1718214950515253545556. 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 57. 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.

Methods

Organism

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 32. 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.

cDNA library

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 58 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 4445. 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).

Validation

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.

Subtraction

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

Computational work was done on an Apple G5 dual 2 GHz processor system with 4 GB of RAM. Data are stored in a MySQL database, perl scripts were used extensively for parsing and loading data, and PHP was used on an APACHE web server to construct the user interface. Additional programs available from their authors are mentioned within context. Software and databases are described in Table 3.

<p>Table 3</p>

Software and Databases. Publicly available software and databases used for FunnyBase. The version and/or download date are listed.

Resource

Version and/or Download Date


Software

Stand-alone BLAST

2.2.8

Cross match

0.990329

CAP3

January, 2004


Public Sequence Similarity Databases

Swiss-Prot

NR

44

Human RefSeq

June, 2004

Human Unigene

June, 2004

Zebrafish Unigene

June, 2004

Medaka Unigene

June, 2004

Rainbow Trout Unigene

June, 2004

June, 2004


Functional Annotation

Gene Ontology

2004-06-04

KEGG

June, 2004

Microarrays

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 59 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).

Authors' contributions

JP designed, scripted and implemented FunnyBase and provided statistical analyses of database. MFO initiated, designed protocols and provided sequences for F. heteroclitus' EST project. JDV optimized robotic interfaces for sequencing and sequenced ESTs. JLR and KJK sequenced ESTs. GJW collaborated on bioinformatics and database development. JAW provided microarray data and analyses. DLC initiated F. heteroclitus' EST project and developed the database and annotation schemes for FunnyBase. All authors read and approved the final manuscript.

Acknowledgements

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.

<p>Gene expression patterns of breast carcinomas distinguish tumor subclasses with clinical implications</p> 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 Proceedings of the National Academy of Sciences of the United States of America 2001 98 10869 10874 11553815 10.1073/pnas.191367098 58566 <p>Systematic variation in gene expression patterns in human cancer cell lines</p> 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 Nature Genetics 2000 24 227 235 10700174 10.1038/73432 <p>Expression profiling reveals distinct sets of genes altered during induction and regression of cardiac hypertrophy</p> Friddle CJ Koga T Rubin EM Bristow J PNAS 2000 97 6745 6750 10829065 10.1073/pnas.100127897 18725 <p>Identification of new genes differentially expressed in coronary artery disease by expression profiling</p> Archacki SR Angheloiu G Tian XL Tan FL DiPaola N Shen GQ Moravec C Ellis S Topol EJ Wang Q Physiol Genomics 2003 15 65 74 12902549 <p>Cardiac hypertrophy by hypertension and exercise training exhibits different gene expression of enzymes in energy metabolism</p> Iemitsu M Miyauchi T Maeda S Sakai S Fujii N Miyazaki H Kakinuma Y Matsuda M Yamaguchi I Hypertens Res 2003 26 829 837 14621187 10.1291/hypres.26.829 <p>Transcriptomal analysis of failing and nonfailing human hearts</p> Steenman M Chen YW Le Cunff M Lamirault G Varro A Hoffman E Leger JJ Physiol Genomics 2003 12 97 112 12429867 <p>Cardiac transcriptional response to acute and chronic angiotensin II treatments</p> Larkin JE Frank BC Gaspard RM Duka I Gavras H Quackenbush J Physiol Genomics 2004 18 152 66 10.1152/physiolgenomics.00057.2004 15126644 <p>Comprehensive analysis of the mouse metabolome based on the transcriptome</p> Bono H Nikaido I Kasukawa T Hayashizaki Y Okazaki Y Genome Res 2003 13 1345 1349 12819132 10.1101/gr.974603 403659 <p>Alterations in mitochondrial function in a mouse model of hypertrophic cardiomyopathy</p> Lucas DT Aryal P Szweda LI Koch WJ Leinwand LA Am J Physiol Heart Circ Physiol 2003 284 H575 83 12414446 <p>Global transcription analysis of Krebs tricarboxylic acid cycle mutants reveals an alternating pattern of gene expression and effects on hypoxic and oxidative genes</p> McCammon MT Epstein CB Przybyla-Zawislak B McAlister-Henn L Butow RA Mol Biol Cell 2003 14 958 972 12631716 10.1091/mbc.E02-07-0422 151572 <p>Genetics and functional genomics of type 2 diabetes mellitus</p> Toye A Gauguier D Genome Biology 2003 4 241 14659011 10.1186/gb-2003-4-12-241 329413 <p>Coordinated patterns of gene expression for substrate and energy metabolism in skeletal muscle of diabetic mice</p> Yechoor VK Patti ME Saccone R Kahn CR Proc Natl Acad Sci U S A 2002 99 10587 10592 12149437 10.1073/pnas.142301999 124982 <p>Transposon-mediated enhancer trapping in medaka</p> Grabher C Henrich T Sasado T Arenz A Wittbrodt J Furutani-Seiki M Gene 2003 322 57 66 14644497 10.1016/j.gene.2003.09.009 <p>A survey of expressed genes in Japanese flounder (Paralichthys olivaceus) liver and spleen</p> Inoue S Nam BH Hirono I Aoki T Molecular Marine Biology & Biotechnology 1997 6 376 380 <p>15,000 Unique zebrafish EST clusters and their future use in microarray for profiling gene expression patterns during embryogenesis</p> Lo J Lee S Xu M Liu F Ruan H Eun A He Y Ma W Wang W Wen Z Peng J Genome Research 2003 13 455 466 12618376 10.1101/gr.885403 430290 <p>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</p> Savan R Sakai M Comparative Biochemistry & Physiology Part B, Biochemistry & Molecular Biology 131B 2002 1 71 82 10.1016/S1096-4959(01)00488-2 <p>The construction of an EST database for Bombyx mori and its application</p> 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 Proceedings of the National Academy of Sciences of the United States of America 2003 100 14121 14126 14614147 10.1073/pnas.2234984100 283556 <p>Development and application of a salmonid EST database and cDNA microarray: data mining and interspecific hybridization characteristics</p> 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 Genome Res 2004 14 478 490 14962987 10.1101/gr.1687304 353236 <p>Annotated Expressed Sequence Tags and cDNA Microarrays for Studies of Brain and Behavior in the Honey Bee</p> Whitfield CW Band MR Bonaldo MF Kumar CG Liu L Pardinas JR Robertson HM Soares MB Robinson GE Genome Res 2002 12 555 566 11932240 10.1101/gr.5302 187514 <p>The utility of natural populations for microarray analyses: isolation of genes necessary for functional genomic studies.</p> Oleksiak MF Kolell K Crawford DL Marine Biotechnology 2001 3 S203 S211 14961317 10.1007/s10126-001-0043-0 <p>CnidBase: The Cnidarian Evolutionary Genomics Database</p> Ryan JF Finnerty JR Nucleic Acids Research 2003 31 159 163 12519972 10.1093/nar/gkg116 165563 <p>Functional genomics does not have to be limited to a few select organisms.</p> Crawford DL Genome Biology 2001 2 INTERACTIONS1001 http://www.genomebiology.com/2001/2/1/interactions/1001/ 10.1186/gb-2001-2-1-interactions1001 <p>Variation in gene expression within and among natural populations</p> Oleksiak MF Churchill GA Crawford DL Nature Genetics 2002 32 261 266 12219088 10.1038/ng983 <p>Gene Expression Profiles in the Brain Predict Behavior in Individual Honey Bees</p> Whitfield CW Cziko AM Robinson GE Science 2003 302 296 299 14551438 10.1126/science.1086807 <p>Expression profiles during honeybee caste determination</p> Evans JD Wheeler DE Genome Biology 2001 2 RESEARCH0001 17597 11178278 <p>Hypoxia-induced gene expression profiling in the euryoxic fish Gillichtys mirabilis</p> Gracey AY Troll JV Somero GN Proceeding of National Academy of Science, USA 2001 98 1993 1998 10.1073/pnas.98.4.1993 <p>Changes in gene expression associated with acclimation to constant temperatures and fluctuating daily temperatures in an annual killifish Austrofundulus limnaeus</p> Podrabsky JE Somero GN J Exp Biol 2004 207 2237 2254 15159429 10.1242/jeb.01016 <p>Differential gene expression in the brain of channel catfish (Ictalurus punctatus) in response to cold acclimation</p> Ju Z Dunham RA Liu Z Molecular Genetics and Genomics 2002 268 87 95 12242503 10.1007/s00438-002-0727-9 <p>Phylogenetic analysis of glycolytic enzyme expression</p> Pierce VA Crawford DL Science 1997 275 256 259 10.1126/science.276.5310.256 <p>Gene flow and mitochondrial DNA variation in the killifish Fundulus heteroclitus</p> Brown BL Chapman RW Evolution 1991 45 1147 1161 <p>Molecular basis of evolutionary adaptation at the lactate dehydrogenase-B locus in the fish Fundulus heteroclitus</p> Crawford DL Powers DA Proceedings of the National Academy of Sciences of the United States of America 1989 86 9365 9369 298496 2594773 <p>A multidisciplinary approach to the selectionist/neutralist controversy using the model teleost, Fundulus heteroclitus</p> Powers DA Smith M Gonzalez-Villasenor I DiMichelle L Crawford DL Bernardi G Lauerman T Oxford Surveys in Evolutionary Biology New York, NY, Oxford University Press Futuyma D and Antonovics J 1993 9 43 108 <p>Phylogenetic analysis of thermal acclimation of the glycolytic enzymes in the genus Fundulus</p> Pierce VA Crawford DL Physiological Zoology 1997 70 597 609 9361133 <p>Phylogenetic analysis of glycolytic enzyme expression.</p> Pierce VA Crawford DL Science 1997 276 256 259 9092475 10.1126/science.276.5310.256 <p>http://www.tigr.org/tdb/tgi/</p> TIGR <p>http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=unigene</p> UniGene <p/> dbESTs NCBI (http://wwwncbinlmnihgov/dbEST/dbEST_summaryhtml) 1 July 2004 <p>Gene ontology: tool for the unification of biology. The Gene Ontology Consortium</p> 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 Nature Genetics 2000 25 25 29 10802651 10.1038/75556 <p>The Kyoto encyclopedia of genes and genomes--KEGG</p> Wixon J Kell D Yeast 2000 17 48 55 10928937 10.1002/(SICI)1097-0061(200004)17:1<48::AID-YEA2>3.0.CO;2-H <p>UniProt: the Universal Protein knowledgebase</p> 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 Nucleic Acids Research 2004 32 Database issue D115 9 14681372 10.1093/nar/gkh131 308865 <p>NCBI NR (non-redudant) Database</p> NCBI <p>Basic local alignment search tool</p> Altschul SF Gish W Miller W Myers EW Lipman DJ Journal of Molecular Biology 1990 215 403 410 2231712 10.1006/jmbi.1990.9999 <p>CAP3: A DNA Sequence Assembly Program</p> Huang X Madan A Genome Research 1999 9 868 877 10508846 10.1101/gr.9.9.868 <p>Generation and analysis of 280,000 human expressed sequence tags</p> 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 Genome Research 1996 6 807 828 8889549 <p>Normalization and subtraction: two approaches to facilitate gene discovery</p> Bonaldo MF Lennon G Soares MB Genome Research 1996 6 791 806 8889548 <p>Online Mendelian Inheritance in Man, OMIM (TM)</p> McKusick-Nathans Institute for Genetic Medicine JHUBMDNCBINLM http://www.ncbi.nlm.nih.gov/omim/ <p>Variation in tissue-specific gene expression</p> Whitehead JA Crawford DL Genome Biology 2004 In Press <p>RefSeq and LocusLink: NCBI gene-centered resources</p> Pruitt KD Maglott DR Nucleic Acids Research 2001 29 137 140 11125071 10.1093/nar/29.1.137 29787 <p>A large database of chicken bursal ESTs as a resource for the analysis of vertebrate gene function</p> Abdrakhmanov I Lodygin D Geroth P Arakawa H Law A Plachy J Korn B Buerstedde JM Genome Res 2000 10 2062 2069 11116100 10.1101/gr.10.12.2062 <p>SGMD: the Soybean Genomics and Microarray Database</p> Alkharouf NW Matthews BF Nucleic Acids Res 2004 32 Database issue D398 400 14681442 10.1093/nar/gkh126 308860 <p>ZmDB, an integrated database for maize genome research</p> Dong Q Roy L Freeling M Walbot V Brendel V Nucleic Acids Res 2003 31 244 247 12519992 10.1093/nar/gkg082 165529 <p>Winter Flounder Expressed Sequence Tags: Establishment of an EST Database and Identification of Novel Fish Genes</p> Douglas SE Gallant JW Bullerwell CE Wolff C Munholland J Reith ME 1999 1 458 464 10525680 <p>Gene discovery using the maize genome database ZmDB</p> Gai X Lal S Xing L Brendel V Walbot V Nucleic Acids Res 2000 28 94 96 10592191 10.1093/nar/28.1.94 102455 <p>GDR (Genome Database for Rosaceae): integrated web resources for Rosaceae genomics and genetics research</p> Jung S Jesudurai C Staton M Du Z Ficklin S Cho I Abbott A Tomkins J Main D BMC Bioinformatics 2004 5 130 15357877 10.1186/1471-2105-5-130 517928 <p>Development of a porcine brain cDNA library, EST database, and microarray resource</p> Nobis W Ren X Suchyta SP Suchyta TR Zanella AJ Coussens PM Physiol Genomics 2003 16 153 159 14559975 10.1152/physiolgenomics.00099.2003 <p>New genomic resources for the honey bee(Apis mellifera L.): development of a deep-coverage BAC library and a preliminary STC database</p> 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 Genet Mol Res 2002 1 306 316 14963821 <p>Natural Variation in Cardiac Metabolism and Gene Expression in Fundulus heteroclitus</p> Oleksiak MF Roach JL Crawford DL Nature Genetics 2004 Advance_Online 28 November 2004; doi:10.1038/ng1483 <p>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</p> Karchner SI Powell WH Hahn ME Journal of Biological Chemistry 1999 274 33814 33824 10559277 10.1074/jbc.274.47.33814 <p>Amplified Rna Synthesized from Limited Quantities of Heterogeneous Complementary DNA</p> Van Gelder RN Von Zastrow ME Yool A Dement WC Barchas JD Eberwine JH Proceedings of the National Academy of Sciences of the United States of America 1990 87 1663 1667 53542 1689846