Large-scale analysis of full-length cDNAs from the tomato (Solanum lycopersicum) cultivar Micro-Tom, a reference system for the Solanaceae genomics
- Koh Aoki†1Email author,
- Kentaro Yano†2,
- Ayako Suzuki2,
- Shingo Kawamura2,
- Nozomu Sakurai1,
- Kunihiro Suda1,
- Atsushi Kurabayashi1,
- Tatsuya Suzuki3,
- Taneaki Tsugane3,
- Manabu Watanabe3,
- Kazuhide Ooga1,
- Maiko Torii1,
- Takanori Narita4,
- Tadasu Shin-i4,
- Yuji Kohara4,
- Naoki Yamamoto2,
- Hideki Takahashi5,
- Yuichiro Watanabe6,
- Mayumi Egusa7,
- Motoichiro Kodama7,
- Yuki Ichinose8,
- Mari Kikuchi9,
- Sumire Fukushima9,
- Akiko Okabe9,
- Tsutomu Arie9,
- Yuko Sato10,
- Katsumi Yazawa10,
- Shinobu Satoh10,
- Toshikazu Omura11,
- Hiroshi Ezura11 and
- Daisuke Shibata1
© Aoki et al; licensee BioMed Central Ltd. 2010
Received: 23 September 2009
Accepted: 30 March 2010
Published: 30 March 2010
The Solanaceae family includes several economically important vegetable crops. The tomato (Solanum lycopersicum) is regarded as a model plant of the Solanaceae family. Recently, a number of tomato resources have been developed in parallel with the ongoing tomato genome sequencing project. In particular, a miniature cultivar, Micro-Tom, is regarded as a model system in tomato genomics, and a number of genomics resources in the Micro-Tom-background, such as ESTs and mutagenized lines, have been established by an international alliance.
To accelerate the progress in tomato genomics, we developed a collection of fully-sequenced 13,227 Micro-Tom full-length cDNAs. By checking redundant sequences, coding sequences, and chimeric sequences, a set of 11,502 non-redundant full-length cDNAs (nrFLcDNAs) was generated. Analysis of untranslated regions demonstrated that tomato has longer 5'- and 3'-untranslated regions than most other plants but rice. Classification of functions of proteins predicted from the coding sequences demonstrated that nrFLcDNAs covered a broad range of functions. A comparison of nrFLcDNAs with genes of sixteen plants facilitated the identification of tomato genes that are not found in other plants, most of which did not have known protein domains. Mapping of the nrFLcDNAs onto currently available tomato genome sequences facilitated prediction of exon-intron structure. Introns of tomato genes were longer than those of Arabidopsis and rice. According to a comparison of exon sequences between the nrFLcDNAs and the tomato genome sequences, the frequency of nucleotide mismatch in exons between Micro-Tom and the genome-sequencing cultivar (Heinz 1706) was estimated to be 0.061%.
The collection of Micro-Tom nrFLcDNAs generated in this study will serve as a valuable genomic tool for plant biologists to bridge the gap between basic and applied studies. The nrFLcDNA sequences will help annotation of the tomato whole-genome sequence and aid in tomato functional genomics and molecular breeding. Full-length cDNA sequences and their annotations are provided in the database KaFTom http://www.pgb.kazusa.or.jp/kaftom/ via the website of the National Bioresource Project Tomato http://tomato.nbrp.jp.
The Solanaceae family comprises 1000-2000 species that show wide morphological variability and ecological adaptability . This taxon includes a number of vegetable crops including fruit-bearing vegetables, tuber-bearing vegetables, and ornamental plants, many of which have economic importance. Tomato (Solanum lycopersicum) is a member of the Solanaceae, and it has served as a model system for fruit development  and plant defense [3, 4]. A number of studies have accumulated substantial information regarding the genetics and physiology of tomato.
To understand the physiological processes of tomato at the molecular level, it is necessary to link traits to DNA sequence. Whole-genome sequencing is one of the approaches to achieve this goal. A project for sequencing the tomato genome was launched by The International Solanaceae Project (SOL) in the year 2003, and the project is currently in the final stage of sequence polishing (March, 2010) . Accumulation of expressed sequence tags (ESTs) is an alternative approach to collecting DNA sequence information. A large number of ESTs have been produced to provide a resource for gene finding and expression studies [5–9]. For tomato ESTs, several repositories are available worldwide, such as the Dana-Farber Cancer Institute (DFCI) Tomato Gene Index, The SOL Genomics Network (SGN) , MiBASE , PlantGDB , and TomatEST , and the number of tomato ESTs accumulated in the databases is 296,957 (January 8, 2010). To elucidate the functions of individual genes, attempts have been made to organize ESTs into consensus sequences, which are referred to as unigenes or tentative consensus (TC) sequences. However, a more direct approach to collecting DNA information of individual genes is to analyze full-length cDNA sequences. Major advantages of sequencing full-length cDNA over the EST assembly include that the data is derived from a single clone rather than the assembly of multiple ESTs which can generate contigs containing sequences from multiple transcripts. Another major advantage of this approach is that full-length cDNA clones have complete sequences of transcripts including coding regions (CDSs) and untranslated regions (UTRs). This facilitates the subsequent annotation and prediction of genomic structures. The complete CDSs in full-length cDNAs allow the prediction of protein sequences, conserved domains, and conserved motifs. Furthermore, full-length cDNAs are easy to use with gene transfer systems, permitting the functional analysis of individual genes through reverse-genetics approaches. Thus, full-length cDNAs are a powerful genomic research tool.
In tomato genomics, a miniature tomato cultivar of S. lycopersicum, Micro-Tom , is regarded as an excellent model system. Micro-Tom has characteristics that make it suitable for experimental study, such as small size, short generation time, and ease of transformation [14, 15]. Recently, resources such as various mutagenized lines have been developed in the Micro-Tom-background . Thus, a collection of Micro-Tom full-length cDNAs would boost the molecular bioengineering of tomato by synergistically combining with other resources. However, no full-length cDNA collection of Micro-Tom, or of tomato, has yet been established.
In this study, we developed a large-scale collection of Micro-Tom full-length cDNAs. We constructed full-length-enriched cDNA libraries from Micro-Tom leaf, fruit, and root tissues, and we obtained full-length sequences of 13,227 cDNAs. We then checked for redundancy between sequences, for the presence of CDSs, and for the contamination of chimeric sequences, thus generating a set of 11,502 non-redundant full-length cDNA sequences (nrFLcDNAs) for sequence analysis. We report the results of UTR analyses, including length distribution, base composition; and the results of CDS analyses, including classification of functions of deduced proteins, metabolic pathway annotations, and comparisons with genes of other plants. We also report the results of mapping nrFLcDNAs onto the draft of tomato whole-genome sequence. Mapping onto genome sequence suggested that introns of tomato genes are relatively long compared to Arabidopsis and rice. We also found that the frequency of nucleotide mismatch between exon regions of Micro-Tom and those of the genome-sequencing cultivar, Heinz 1706, is 0.061%. Finally, the value of Micro-Tom full-length cDNA information will be discussed in terms of integration with the tomato genome sequence and other tomato resources.
Results and discussion
Construction of full-length-enriched cDNA libraries
Micro-Tom tissues used for RNA preparation
CMV(strain TN) and satellite RNA
Inoculated leaves, systemic leaves
16 dpi (21-day old)
ToMV (strain Lta1, L11A, L2a, Ltb1, LJB)
Inoculated leaves, systemic leaves
7, 14 dpi (28, 35-day old)
Non-pathogenic Alternaria alternata
6, 24 hpi (5-week old)
Alternaria alternate f. sp. lycopersici (As-27)
6, 24 hpi (5-week old)
Corynespora cassiicola (isolated from tomato and cucumber)
6, 24 hpi (3-week old)
Cladosporium fulvum (strain 210, 211, 217)
5 dpi (40-day old)
Fusarium oxysporum f. sp. lycopersici race 1 and 2
Inoculated to root, 7 dpi (4-week old)
Non-pathogenic Fusarium oxysporum
Inoculated to root, 7 dpi (4-week old)
Pseudomonas syringae pv. tomato (DC3000 wild type), pv tabaci isolate 6605 (wild type, ΔfliC, ΔfliD)
24 hpi (6-week old)
Flagellin from P. syringae pv. tabaci isolate 6605, fgl22
3 hps (6-week old)
Probenazole (100 μg/ml)
4 dps (4-week old)
Validamycin A (100 μg/ml)
4 dps (4-week old)
Acibenzolar-S-methyl (100 μg/ml)
4 dps (4-week old)
Validoxylamine A (100 μg/ml)
4 dps (4-week old)
Salicilic acid (0.5 mM)
2 dps (4-week old)
Methyl jasmonic acid (50 μM)
2 dps (4-week old)
40 daa, year 2003, 2004
42 daa, year 2003, 2004
48 daa, year 2003, 2004
50 daa, year 2003, 2004
Roots from plants without flower
6, 10-week old
Roots from flowered plants
10, 12- week old
Fusarium oxysporum race 2
1, 7, 14, 21 dpi, 4- to 11-week old
Full-length sequencing of cDNAs
Full-length sequencing was performed using the representative clones by the clone-by-clone primer walking method. Out of 22,900 clones, full-length inserts of 13,227 cDNA clones were sequenced (see Additional file 2: Clone number list). Mean of phred quality value of each base was 65 (i.e., one error in 3 × 106 bases), and more than 70% of bases exceeded quality value 68. This collection was registered to the high-throughput cDNA sequence (HTC) division of the DDBJ.
Generation of a non-redundant set of the full-length cDNAs
We checked for redundancy within the set of 13,227 HTCs, and then a set of 12,106 HTCs representing non-redundant transcripts was sent for further processing (see Additional file 2: Clone number list). To check whether each of the 12,106 HTCs has a full-length CDS or not, we first searched for cDNAs derived from non-coding RNAs using the NONCODE version 2.0 dataset . Three HTCs (FC25DB10, LEFL3003O13, and LEFL3054J22) were identified as cDNAs derived from non-coding RNA (see Additional file 3: cDNA derived from ncRNA) and excluded from further analysis. Next, we confirmed that all of the HTCs did not match DNA sequences from pathogens used to attack leaf and root tissues, suggesting that the probability of contamination with pathogen-derived transcripts was negligible (see Additional file 2: Clone number list, for the accession numbers of pathogen-derived sequences). We then checked for chimeric clones putatively generated during cDNA cloning process. According to the criteria described in Methods, 36 HTCs were regarded as chimeric clones (see Additional file 2: Clone number list). We then identified sets of HTCs generated by the alternative splicing event, and excluded 470 HTCs containing retained introns (see Additional file 2: Clone number list). From the remaining 11,597 HTCs, CDSs were predicted by using FrameDP  and GeneMark.hmm-E  programs, and by selecting open reading frame encoding the longest amino acid sequence. Essentially, a CDS encoding amino acid sequence that had the highest similarity to proteins registered in either the nr or the tomato SBM protein databases was selected as representative CDS of the HTC (Additional file 4: Predicted CDS of 11,597 non-redundant HTCs). HTCs with CDSs shorter than nine amino acid (aa)long were excluded from further analysis, since the shortest length tomato protein registered in UniProt was 10 aa (UniProt accession number, Q6TS30; description, ENOD40). Out of 11,597 non-redundant HTCs, CDS encoding full-length protein was predicted for 11,502 HTCs. Full-length CDS were not predicted for 73 HTCs and no coding regions were predicted for 22 HTCs. Finally, a set of 11,502 HTCs (see Additional file 2: Clone number list) was generated, hereafter referred to as nrFLcDNA. This nrFLcDNA set was subjected to subsequent sequence analyses.
Length distributions of UTRs and CDSs
Length of cDNA inserts, CDSs, and UTRs, and comparison with full-length cDNAs from other plants
Figures 3C and 3D show the distributions of the 5'- and 3'-UTR lengths, respectively. In contrast to CDS length, tomato nrFLcDNAs had longer UTRs compared to those of other plants. The average 5'- and 3'-UTR lengths were 175 and 257 bp, respectively. The median 5'- and 3'-UTR lengths were 106 and 220 bp, respectively. On both average and median bases, the UTRs of tomato are longer than those of Arabidopsis, soybean, poplar, and maize, and slightly shorter than those of rice (Table 2). The relatively long UTRs of tomato likely contribute to the regulation of mRNA transcription, translation, and stability .
Nucleotide composition for different parts of tomato full-length cDNAs
Alternative splicing and retained intron
Out of these 1206 non-redundant HTCs, 388 HTCs were aligned with SGN tomato BAC sequences (bacs.v374.seq.20081128091837 downloaded from the SGN ftp site http://sgn.cornell.edu/bulk/input.pl?mode=ftp), which allowed us to classify the types of 434 alternative splicing events (see Additional file 2: Clone number list). The retained intron-type was the most frequent (55% of all events). The alternative acceptor site- and alternative donor site-types also occurred frequently, comprising 20% and 16% of all events, respectively. The exon skip-type and other types (e.g., alternative terminal exon) were rather rare. The occurrence of alternative acceptor site- and alternative donor site-types was slightly higher than that reported previously for tomato . nrFLcDNAs with retained introns encode proteins with alternative amino acid sequences. We found that 509 nrFLcDNAs may contain retained full- or partial-length introns (see Additional file 2: Clone number list). Larger scale identification of alternative splicing events using the draft of tomato whole-genome sequence is currently underway.
Functional annotation of tomato full-length cDNA
Transcription factors found in the nrFLcDNAs
Description of IPR entry
No. of tomato nrFLcDNAs
Myb-like, SHAQKYF class
Myb transcription factor
Zinc finger, C2H2-type
CCAAT-HAP5, HAP3, and HAP2
Transcription factor CBF/NF-Y/archaeal histone
Zinc finger, PHD-type
Pathogenesis-related transcriptional factor and ERF, DNA-binding
bZIP transcription factor, bZIP-1
Heat shock factor (HSF)-type, DNA-binding
Zinc finger, CCCH-type
Transcription factor, MADS-box
Myb-related protein P, C-terminal
High mobility group box
No apical meristem (NAM) protein
Transcription factor, TCP
GRAS transcription factor
Zinc finger, GATA-type
Transcriptional factor B3
Ethylene insensitive 3
Zinc finger, Dof-type
Zinc finger, B-box, "IPR002926 Zinc finger, CONSTANS-type" was deleted from IPR entry.
ZF-HD homeobox protein Cys/His-rich dimerisation region
Protein of unknown function DUF573
Response regulator receiver
Plant regulator RWP-RK
BZR1, transcriptional repressor
Transcription factor jumonji
Zinc finger, Zim17-type
ssDNA-binding transcriptional regulator
Comparative analysis with genes of other plants
To obtain insight into functions of nrFLcDNAs that were not found in other plants, we extracted 694 nrFLcDNAs that matched none of the other plant genes (see Additional file 2: Clone number list). More than 90% of the 694 nrFLcDNAs matched the prerelease of tomato genome shotgun sequence with E-value smaller than 1e-120 and only 27 nrFLcDNAs did not match genome sequence with E-value threshold of 1e-10, suggesting that most of the 694 nrFLcDNAs are tomato-derived transcripts. By searching the InterPro database, out of the 694 nrFLcDNAs, 54 nrFLcDNAs were assigned with 75 protein domains. We also identified known protein domains by searching the similarity with the tomato SBM protein datasets. Out of 48 nrFLcDNAs matched (E-value < 1e-10) SBM proteins, 19 nrFLcDNA had known protein domains. Domains such as Cyclin-like F-box (IPR001810), Aldehyde dehydrogenase (IPR015590), and Polynucleotidyl transferase, Ribonuclease H fold (IPR012337) were identified in both searches. However, known functional protein domains were not found in most (640 out of 694) of the nrFLcDNAs that were not found in other 16 plants.
Mapping full-length cDNAs onto tomato genomic sequence
To investigate the occurrence of introns in UTRs, the positions of translation start sites and translation termination sites were mapped to the exons. The frequency of the first methionine in the internal exons was 19.9%. On the other hand, the frequency of stop codons in the internal exons was 5.7%. This result indicates that intron mapping to a UTR was more likely to occur in 5'-UTRs than in 3'-UTRs.
This mapping also provides insights into the single nucleotide polymorphism frequency between Micro-Tom and the genome-sequencing cultivar, Heinz 1706. To estimate frequency of nucleotide mismatch, we used nrFlcDNA-scaffold pairs with identity in the BLAST output equal to or larger than 99.5%. Resulting set of nrFlcDNA-scaffold pairs contained 90% of the total pairs. Frequency of nucleotide mismatches in exon regions was estimated to be 0.061% (i.e., one difference in 1640 nt). This is nearly comparable to the previous result based on the comparison of ESTs between Micro-Tom and other S. lycopersicum cultivars . This result suggests that the Micro-Tom cDNA sequence serves as a good reference for the tomato genome sequence. We note that the results presented here will be revised when annotation of the tomato genome sequence is finished in near future.
A tomato full-length cDNA database: KaFTom
Information regarding the13,227 full-length cDNA sequences including annotations and the results of similarity searches are distributed from the KaFTom database http://www.pgb.kazusa.or.jp/kaftom/. KaFTom also provides the results of mapping full-length sequences onto SGN tomato BAC sequences and predictions of exons and introns. It directly links to the clone request site, and all of the full-length cDNA clones (89,872 clones) are available from the National Bioresource Project Tomato http://tomato.nbrp.jp.
Significance of Micro-Tom full-length cDNA in tomato genomics
Full-length cDNAs serves as a valuable tool that will accelerate tomato genomics in several respects. First, the sequence information of the full-length cDNAs will be integrated into the whole-genome sequence information of tomato, and it will help genome annotation and the identification of regulatory elements in UTRs. Second, sequences of the full-length cDNAs are valuable information for generating molecular markers. Third, full-length cDNAs help our understanding of the functions of tomato genes through gain-of-function and loss-of-function analyses. Fourth, full-length cDNAs promote the use of "targeting induced local lesions in genome" (TILLING) screening. By combining the use of full-length cDNA sequence information with EMS-mutagenized lines of Micro-Tom , the TILLING approach will provide an efficient way to screen non-transgenic mutant lines as parental germplasm favorable for breeding programs.
We developed a set of 13,227 full-length cDNAs from the model tomato cultivar Micro-Tom, and then we generated a set of 11,502 nrFLcDNAs, each of which represents a non-redundant transcript coding a full-length CDS. Analysis of the 5'-UTRs, CDSs, and 3'-UTRs demonstrated that tomato transcripts have longer 5'- and 3'-UTRs than other plant species. Classification of functions of deduced proteins according to the GO annotation revealed that the nrFLcDNA set covers a broad range of proteins. Comparison of the nrFLcDNAs with genes of other plants facilitated the identification of tomato cDNAs that have very low similarity to 16 other plants tested, most of which did not have known protein domains. Mapping of the nrFLcDNAs onto tomato genome sequences revealed that tomato genes have longer introns than Arabidopsis and rice. The frequency of single nucleotide polymorphisms between exons of Micro-Tom and the genome-sequencing cultivar, Heinz 1706, was estimated to be 0.061%. In combination with other resources, the Micro-Tom full-length cDNAs will bridge the gap between basic and applied studies, providing a valuable tool not only for tomato whole-genome annotation but also for genomics studies and molecular breeding.
The miniature tomato cultivar, Micro-Tom (Solanum lycopersicum cv Micro-Tom), was used to construct full-length cDNA libraries. Total RNA samples were prepared from leaves, fruits, and roots of Micro-Tom. The tissues used for RNA preparation are summarized in Table 1. For a more detailed description of tissue samples, see Additional file 1: Micro-Tom tissues used for RNA preparation.
Leaves were treated with various pathogens. For fungal and bacterial inoculation, Micro-Tom plants were grown in sterilized soil at 20-28°C with a photoperiod of 14-16 h light (7000 lux, or 85 μmol/m2 s)/8-10 h dark. For virus inoculation, Micro-Tom plants were grown in a mixture of vermiculite and perlite (1:1 (v/v)) at 23-25°C with a photoperiod of 16 h light (8000 lux, or 97 μmol/m2 s)/8 h dark, and irrigated with a nutrient medium every three days. Pathogen treatments were performed as described previously .
Fruits were harvested at four ripening stages, mature green, breaker, turning, and red ripe stages. We used two harvests of Micro-Tom fruit for fruit RNA preparation. In the first harvest, Micro-Tom plants were grown in the year 2003 and fruit pericarp was harvested as described previously . In the second harvest, Micro-Tom plants were grown in the year 2004 and fruit pericarp was harvested as described previously .
Roots were harvested from plants before and after flowering, and from plants treated with Fusarium oxyporum race 2. For root samples from plants before and after flowering, Micro-Tom seeds were sterilized with 70% (v/v) ethanol for 1 min, and germinated on sterile MS-agar plates (3% (w/v) sucrose) in a photoperiod of 12 h light (7000 lux, or 85 μmol/m2 s)/12 h dark at 26°C. On the 10th day after germination, seedlings were transferred to a sterile mixture of vermiculite and Powersoil (1:2) (Kureha Chemical Industries, http://www.kureha.co.jp/ and Kanto Hiryo Industries, http://www.okumurashoji.co.jp/) and were grown in a photoperiod of 12 h light (7000 lux, or 85 μmol/m2 s)/12 h dark at 26°C. Roots were harvested from approximately three-month-old plants before or after flowering. For F. oxyporum race 2 treatment, Micro-Tom seeds were germinated on sterile soil, grown in sterile conditions, and treated with F. oxyporum race 2 as described previously .
Construction of the full-length-enriched cDNA library
Equal amounts of total RNA (approximately 60 μg each) from 42 leaf tissue types were mixed and then used to construct the LEFL1 cDNA library. Equal amounts of total RNA (approximately 2400 μg each) from four fruit tissues harvested in the year 2004 were mixed and then used to construct the LEFL2 cDNA library. Likewise, equal amounts of total RNA (approximately 230 μg each) from 12 root tissue types were mixed and then used to construct the LEFL3 cDNA library. Construction of the above-mentioned libraries was performed by the biotinylated CAP-trapper method as described previously . For the RNA mixture from fruits harvested in the year 2003, a full-length-enriched cDNA library, namely the FC library, was constructed by the vector-capping method described previously . Finally, cDNA inserts were cloned into the pFLCIII vector or the pGCAPzf3 vector.
5'-end sequencing and clustering
cDNA clones randomly collected from the libraries were single-pass sequenced from the 5'-end as described previously . From the 5'-end sequences, vector-derived sequences and low-quality regions (Phred quality value < 30) were trimmed using a combination of the cross_match program http://www.phrap.org and the phred program . Sequence data of lengths shorter than 50 bases were not included in further data processing. If the repetition of a single nucleotide was longer than 10% of a sequence, such sequences were not included in further data processing. The 5'-end sequences were combined with previously sequenced Micro-Tom ESTs  and ESTs registered in the SGN database (file name: tomato_species_2008_10_21.seq, obtained from the SGN ftp site ftp://ftp.sgn.cornell.edu/est_sequences/species/Tomato/), then clustered using the phrap program http://www.phrap.org/phredphrapconsed.html. Clusters containing 5'-end sequences of clones derived from the FC and LEFL libraries were chosen. Based on the sequence alignment within each cluster, an FC or LEFL clone carrying the insert that has the longest extension in the 5'-direction was selected as a representative of the cluster. The representative clones were sent for full-length sequencing. All of the 5'-end sequences were registered to the DDBJ http://www.ddbj.nig.ac.jp/index-e.html with accession numbers BW684914-BW692959, DB678259-DB727670, and FS179211-206426.
The clones selected as the representatives of the clusters were re-arrayed in 96-well plates. A small aliquot of each clone was inoculated to 2 ml of LB medium containing 50 μg/ml ampicillin, and grown overnight at 37°C in 96-deep-well plates. Plasmid DNA was prepared from the overnight culture using CosMCPrep (Agencourt Bioscience, Beverly, MS, USA) according to the manufacturer's protocol. By using the plasmid DNA as templates, full-length sequencing was performed by a clone-by-clone primer walking method using a Model 3730 × l sequence analyzer (Applied BioSystems, Forster City, CA, USA). Sequencing reaction was performed using the BigDye Terminator version 3.1 cycle sequencing kit (Applied BioSystems). Sequences were assembled using the phrap program http://www.phrap.org/phredphrapconsed.html. Based on the assembled sequences, primers corresponding to the terminal sequences were designed using the Primer3 program . The primer walking procedure was repeated until a poly(A) tail or a vector-derived sequence appeared. The full-length sequences were finished by trimming the vector-derived sequences from both ends using the cross_match program. Full-length cDNA sequences (13,227 sequences) obtained by this procedure were registered to the high-throughput cDNA sequence (HTC) division of the DDBJ with accession numbers AB211519-211522, AB211526, AK224591-AK224910, AK246135-AK248077, and AK319176-330134.
Selection of a non-redundant set of full-length cDNAs
To identify redundant HTCs, pair-wise similarity searches between members of the 13,227 HTC sequences were performed using BLASTN. Redundant HTCs were identified using the following criteria; E-value < 1e-180, nucleotide identity ≥ 97%, and the percentage of alignment length exceeded 80% of the length of each HTC. The resulting non-redundant HTCs (12,105 HTCs) were then checked for whether they contain clones derived from non-coding RNA, clones derived from pathogen transcripts, whether they contain chimeric clones, and whether they contain clones with retained introns. HTCs derived from non-coding RNAs were identified by searching for similarity of the HTCs against the "ncrna_NONCODE[v2.0].fasta" dataset obtained from the NONCODE download site http://www.noncode.org/download.htmusing BLASTN with a threshold E-value of 1e-180. To search HTCs derived from pathogen-derived transcripts, nucleotide sequences from Alternalia alternate, Cladosporium fulvum, Corynespora cassiicola, Cucumber mosaic virus, Fusarium oxysporum, Pseudomonas syringae, and Tomato mosaic virus were retrieved from NCBI nucleotide database (see Additional file 2: Clone number list for accession numbers of pathogen-derived sequences). Similarity of the HTCs to pathogen nucleotide sequences were searched using BLASTN with threshold E-value of 1e-30. When the E-value was smaller than 1e-30, the nrFLcDNA was further checked for whether or not they have similar SGN tomato unigenes or DFCI tomato TCs with smaller E-values. If so, the nrFLcDNA was regarded as cDNA derived from Micro-Tom, not from pathogens. To identify chimeric clones, the full-length cDNA sequences were searched against SGN tomato unigenes 'draft_Solanum_lycopersicum_transcript_assembly_2009_08_05.fasta' with cut-off E-value of 1e-50. HTCs that matched more than two unigenes having different functional descriptions were identified as candidate chimeric clones. Out of these candidates, if the candidate clone matched to two or more scaffolds of the prerelease of tomato genome shotgun sequence (S_lycopersicum_scaffolds_20091201.fa) downloaded from the SGN ftp site ftp://ftp.solgenomics.net/tomato_genome/wgs/assembly, the candidate was regarded as a chimeric clone. Out of the remaining candidates, if a candidate HTC had either a BamHI or XhoI site (which was used for ligation of cDNA to vectors) in between regions matched to unigenes, the candidate was regarded as a chimeric clone. To identify HTCs containing retained intron, HTC sequences were subjected to similarity searches against the HTC itself using BLASTN with threshold E-value of 1e-50. To detect introns, pair-wise alignments of each target nrFLcDNA with members of that multi-sequence group were carried out using est2genome . HTCs derived from non-coding RNAs, pathogen-derived transcripts, chimeric clones, and intron-containing clones were excluded and non-redundant 11,597 HTCs were subjected to the prediction of coding sequences.
Amino acid sequences were predicted in sense orientation using three methods. First, the CDS coding the longest amino acid sequence was predicted. Second, CDSs were predicted using FrameDP http://iant.toulouse.inra.fr/FrameDP/. Parameters were set as follows: framed_minimum_peptide length, 10; reference protein database, TAIR9_pep_20090619; e-value cut-off for considering ncbi-blastx hits, 1e-3; the method used for the first classification, GC3; the maximum number of models, 3; the maximum number of iteration, 3; cut-offs for considering ncbi-blastx HSPs, e-value 1e-4 and length 100. Third, CDSs were predicted using GeneMark.hmm eukaryotic version 3.3 . The model used was a_thaliana.mod which was supported by the program. Amino acid sequences of all possible frames were then subjected to BLASTP search (threshold E-value was 1e-10) against NCBI nr and the protein dataset (protein_sequence) of Tomato SBM Database http://www.kazusa.or.jp/tomato/. CDSs with the smallest E-value against proteins registered in either nr or SBM proteins were selected. If CDS with the smallest E-value did not cover full-length protein, CDS with the second smallest E-value were selected. If E-values of all CDSs were the same, the longest CDSs were selected. Whether or not HTCs have full-length CDS was assessed by checking if the amino acid sequences contain both start and stop codons. Resulting set of 11,502 non-redundant HTCs encoding full-length proteins was referred to as 'nrFLcDNA.'
UTRs were identified according to the selected CDS. 5'-UTRs were defined as the nucleotide sequence upstream of the start codon. 3'-UTRs were defined as the nucleotide sequence downstream of the stop codon to the poly(A) tail. Nucleotide composition analysis was performed essentially by counting A, T, C, G, N, and X using an in-house Perl program.
Identification of retained introns and alternative splicing
After removing poly(A) tails, a set of 12,106 non-redundant HTC sequences were subjected to similarity searches against non-redundant HTCs, the SGN tomato unigenes (Tomato_200607_build_1, ftp://ftp.sgn.cornell.edu/unigene_builds), and the DFCI Tomato Gene Index (release 12.0, http://compbio.dfci.harvard.edu/tgi/) using BLASTN with threshold E-value of 1e-50. Sequences with similarity to a given target nrFLcDNA were grouped into a multi-sequence group. To detect introns, pair-wise alignments of each target non-redundant HTCs with members of that multi-sequence group were carried out using est2genome. Intron detection was performed bi-directionally, that is, nrFLcDNA/TC/unigene and a target nrFLcDNA were set as "genome" and "est," and vice versa. A BAC sequence file named "bacs.v374.seq.20081128091837" was downloaded from the SGN ftp site http://sgn.cornell.edu/bulk/input.pl?mode=ftp. The target HTCs were mapped onto the BAC sequence with the threshold E-value of 1e-50.
Annotation of full-length DNA sequences
Functional annotations for the nrFLcDNAs were provided according to sequence similarity with public sequence datasets. To estimate the similarity in a protein sequence, nrFLcDNA sequences were queried against protein datasets from NCBI nr, UniProt http://www.uniprot.org/, Arabidopsis protein sequences (TAIR9, http://www.arabidopsis.org/), and rice protein sequences (RAP-DB, http://rapdb.dna.affrc.go.jp/) using the BLASTX and BLASTP. To estimate the similarity to tomato unigenes, nrFLcDNA sequences were queried against the unigene/TC datasets from SGN unigene sequences (draft_Solanum_lycopersicum_transcript_assembly_2009_08_05.fasta, ftp://ftp.sgn.cornell.edu/unigene_builds) and the DFCI Tomato Gene Index (release 12.0, http://compbio.dfci.harvard.edu/tgi/) using the BLASTN. Searches for the protein domains in the amino acid sequences predicted from six open reading frames were performed using InterProScan . Gene ontology (GO) annotations  for Arabidopsis genes with the highest similarity to nrFLcDNAs were retrieved from a TAIR GO annotation search http://www.arabidopsis.org/tools/bulk/go/index.jsp. GO annotation documents for InterProScan domains including "interpro2go," "entry.list," and "ParentChildTreeFile.txt" were obtained at the InterPro public ftp site ftp://ftp.ebi.ac.uk/pub/databases/interpro/ in May 2009. Data files of Arabidopsis transcription factors (file named AtTFDB.zip) were downloaded from AGRIS AtTFDB http://arabidopsis.med.ohio-state.edu/AtTFDB/. The difference in the frequency of GO terms was tested by Fisher's exact probability test (P < 0.01).
The SGN tomato unigene (Tomato_200607_build_1, downloaded from SGN ftp site ftp://ftp.sgn.cornell.edu/unigene_builds/) with the highest similarity to each nrFLcDNAs was searched for using BLASTN, and then SGN tomato unigenes that matched nrFLcDNAs at the threshold E-value of < 1e-180 were assigned to tomato metabolic pathways according to LycoCyc version 126.96.36.199 ftp://ftp.sgn.cornell.edu/pathways/lycocyc.dump.txt. Pathway hierarchy was retrieved on May 1, 2009 from "Hierarchical pathway ranking for Solanum lycopersicum" http://solcyc.sgn.cornell.edu/LYCO/hierarchical.html. The version of LycoCyc used in this study provides pathway annotations to 1532 SGN unigenes. The difference in the frequency of pathway annotation terms was tested by Fisher's exact probability test (P < 0.05).
Similarity to genes of other plants
The similarity of nrFLcDNAs to genes of other plants was estimated using tBLASTN with E-value cut-off of 1e-10 against DFCI Gene Indeices of Hordeum vulgare (barley), Triticum aestivum (wheat), Zea mays (maize), Pinus (pine), Picea (spruce), Populus (poplar), Lotus japonicus, Medicago truncatula, Glycine max (soybean), Citrus sinensis (orange), Malus × domestica (apple), Vitis vinifera (grape), Nicotiana tabacum (tobacco), and Solanum tuberosum (potato) (see Additional file 5: Datasets used for comparison with other plants). GO annotations for DFCI tomato TCs ("LGIGO.071608") were downloaded from the DFCI public ftp site. The difference in the frequency of GO terms was tested by Fisher's exact probability test (P < 0.05).
Mapping full-length cDNA sequences to tomato genome sequence
Poly(A) tails of the nrFLcDNA sequences were excised, and this set of nrFLcDNAs without poly(A) tails was used for similarity searches against the Prerelease of Tomato Genome Shotgun Sequence (S_lycopersicum_scaffolds_20091201.fa, ftp://ftp.solgenomics.net/tomato_genome/wgs/assembly). Tomato genome scaffold sequences with the highest similarity to a given nrFLcDNA sequence were identified using BLASTN (threshold E-value, 1e-50). Each pair of scaffold and nrFLcDNA sequences was then submitted to exon and intron prediction using the est2genome program  by setting scaffold and nrFLcDNA sequences as "genome" and "est," respectively. To analyze numbers and lengths of exons and introns, cut-off value for the alignment length was set as 90% of nrFLcDNA length. Cut-off value for nucleotide identity was set as 90%. To estimate the single nucleotide polymorphisms (SNPs) between nrFLcDNA and scaffold sequences, the number of nucleotide mismatches was counted in the nrFLcDNA-scaffold pairs having 'identity' in BLAST output equal to or larger than 99.5%.
The data presented in this study, including full-length cDNA sequences of the 13,227 HTCs, the results of similarity searches against nr, UniProt, TAIR, RAP-DB, SGN tomato unigenes, DFCI tomato TCs, and InterProScan, and mapping of full-length sequences against SGN tomato BAC sequences are available at the tomato full-length cDNA database KaFTom http://www.pgb.kazusa.or.jp/kaftom/. KaFTom is accessible directly or via the National Bioresource Project Tomato website http://tomato.nbrp.jp.
high throughput cDNA sequence
non-redundant full-length cDNA
Dana-Farber Cancer Institute
SOL genomics network
We are grateful for and acknowledge use of the draft tomato genome sequence, which was generated by the International Tomato Genome Sequencing Consortium http://solgenomics.net/tomato/. We are also grateful for and acknowledge use of the tomato SBM dataset, which was generated by Kazusa DNA Research Institute. We thank Hideki Hirakawa and Shinobu Nakayama (Kazusa DNA Res. Inst.) for assistance in GeneMark.hmm analysis. We thank Shusei Sato (Kazusa DNA Res. Inst.) for critical reading of the manuscript, and Kenta Shirasawa (Kazusa DNA Res. Inst.) for helpful discussion. We also thank Tsugumi Isozaki, Miyuki Inde, and Tsurue Aoyama (Kazusa DNA Res. Inst.) for plant care and lab assistance. We thank Hiroshi Otani (Tottori Univ.), Wataru Hasama (Oita Pref. Agric. Res. Center) and Hiroshi Shiomi (Takii and Co. Ltd.) for providing the Corynespora cassiicola isolate. This work was supported by National Bioresource Project, Genome program, "Enhancing tomato resources by sequencing Micro-Tom full-length cDNA" (2008, MEXT, Japan) to KA, by the Japan Solanaceae Consortium (JSOL), a grant from Meiji Univ. to KY, and a grant from the Kazusa DNA Res. Inst. to DS and KA.
- Knapp S: Tobacco to tomatoes: a phylogenetic perspective on fruit diversity in the Solanaceae. J Exp Bot. 2002, 53: 2001-2022. 10.1093/jxb/erf068.PubMedView ArticleGoogle Scholar
- Giovannoni JJ: Genetic regulation of fruit development and ripening. Plant Cell. 2004, 16 (Suppl): S170-180. 10.1105/tpc.019158.PubMed CentralPubMedView ArticleGoogle Scholar
- Pedley KF, Martin GB: Molecular basis of Pto-mediated resistance to bacterial speck disease in tomato. Annu Rev Phytopathol. 2003, 41: 215-243. 10.1146/annurev.phyto.41.121602.143032.PubMedView ArticleGoogle Scholar
- Mueller LA, Tanksley SD, Giovannoni JJ, van Eck J, Stack S, Choi D, Kim BD, Chen M, Cheng Z, Li C: The Tomato Sequencing Project, the First Cornerstone of the International Solanaceae Project (SOL). Comp Funct Genomics. 2005, 6: 153-158. 10.1002/cfg.468.PubMed CentralPubMedView ArticleGoogle Scholar
- Mueller LA, Solow TH, Taylor N, Skwarecki B, Buels R, Binns J, Lin C, Wright MH, Ahrens R, Wang Y: The SOL Genomics Network: a comparative resource for Solanaceae biology and beyond. Plant Physiol. 2005, 138: 1310-1317. 10.1104/pp.105.060707.PubMed CentralPubMedView ArticleGoogle Scholar
- Hoeven Van der R, Ronning C, Giovannoni J, Martin G, Tanksley S: Deductions about the number, organization, and evolution of genes in the tomato genome based on analysis of a large expressed sequence tag collection and selective genomic sequencing. Plant Cell. 2002, 14: 1441-1456. 10.1105/tpc.010478.PubMed CentralPubMedView ArticleGoogle Scholar
- Alba R, Fei Z, Payton P, Liu Y, Moore SL, Debbie P, Cohn J, D'Ascenzo M, Gordon JS, Rose JK: ESTs, cDNA microarrays, and gene expression profiling: tools for dissecting plant physiology and development. Plant J. 2004, 39: 697-714. 10.1111/j.1365-313X.2004.02178.x.PubMedView ArticleGoogle Scholar
- Fei Z, Tang X, Alba RM, White JA, Ronning CM, Martin GB, Tanksley SD, Giovannoni JJ: Comprehensive EST analysis of tomato and comparative genomics of fruit ripening. Plant J. 2004, 40: 47-59. 10.1111/j.1365-313X.2004.02188.x.PubMedView ArticleGoogle Scholar
- Yamamoto N, Tsugane T, Watanabe M, Yano K, Maeda F, Kuwata C, Torki M, Ban Y, Nishimura S, Shibata D: Expressed sequence tags from the laboratory-grown miniature tomato (Lycopersicon esculentum) cultivar Micro-Tom and mining for single nucleotide polymorphisms and insertions/deletions in tomato cultivars. Gene. 2005, 356: 127-134. 10.1016/j.gene.2005.04.026.PubMedView ArticleGoogle Scholar
- Yano K, Watanabe M, Yamamoto N, Tsugane T, Aoki K, Sakurai N, Shibata D: MiBASE: a database of a miniature tomato cultivar Micro-Tom. Plant Biotechnol. 2006, 23: 195-198.View ArticleGoogle Scholar
- Dong Q, Lawrence CJ, Schlueter SD, Wilkerson MD, Kurtz S, Lushbough C, Brendel V: Comparative plant genomics resources at PlantGDB. Plant Physiol. 2005, 139: 610-618. 10.1104/pp.104.059212.PubMed CentralPubMedView ArticleGoogle Scholar
- D'Agostino N, Aversano M, Frusciante L, Chiusano ML: TomatEST database: in silico exploitation of EST data to explore expression patterns in tomato species. Nucleic Acids Res. 2007, 35: D901-905. 10.1093/nar/gkl921.PubMed CentralPubMedView ArticleGoogle Scholar
- Scott JW, Harbaugh BK: Micro-Tom-a miniature dwarf tomato. Florida Agr Expt Sta Circ. 1989, 370: 1-6.Google Scholar
- Meissner R, Jacobson Y, Melamed S, Levyatuv S, Shalev G, Ashri A, Elkind Y, Levy AA: A new model system for tomato genetics. Plant J. 1997, 12: 1465-1472. 10.1046/j.1365-313x.1997.12061465.x.View ArticleGoogle Scholar
- Shibata D: Genome sequencing and functional genomics approaches in tomato. J Gen Plant Pathol. 2005, 71: 1-7. 10.1007/s10327-004-0150-7.View ArticleGoogle Scholar
- Matsukura C, Aoki K, Fukuda N, Mizoguchi T, Asamizu E, Saito T, Shibata D, Ezura H: Comprehensive resources for tomato functional genomics based on the miniature model tomato micro-tom. Curr Genomics. 2008, 9: 436-443. 10.2174/138920208786241225.PubMed CentralPubMedView ArticleGoogle Scholar
- Kato S, Ohtoko K, Ohtake H, Kimura T: Vector-capping: a simple method for preparing a high-quality full-length cDNA library. DNA Res. 2005, 12: 53-62. 10.1093/dnares/12.1.53.PubMedView ArticleGoogle Scholar
- Carninci P, Kvam C, Kitamura A, Ohsumi T, Okazaki Y, Itoh M, Kamiya M, Shibata K, Sasaki N, Izawa M: High-efficiency full-length cDNA cloning by biotinylated CAP trapper. Genomics. 1996, 37: 327-336. 10.1006/geno.1996.0567.PubMedView ArticleGoogle Scholar
- He S, Liu C, Skogerbo G, Zhao H, Wang J, Liu T, Bai B, Zhao Y, Chen R: NONCODE v2.0: decoding the non-coding. Nucleic Acids Res. 2008, 36: D170-172. 10.1093/nar/gkm1011.PubMed CentralPubMedView ArticleGoogle Scholar
- Gouzy J, Carrere S, Schiex T: FrameDP: sensitive peptide detection on noisy matured sequences. Bioinformatics. 2009, 25: 670-671. 10.1093/bioinformatics/btp024.PubMed CentralPubMedView ArticleGoogle Scholar
- Lomsadze A, Ter-Hovhannisyan V, Chernoff YO, Borodovsky M: Gene identification in novel eukaryotic genomes by self-training algorithm. Nucleic Acids Res. 2005, 33: 6494-6506. 10.1093/nar/gki937.PubMed CentralPubMedView ArticleGoogle Scholar
- Yamada K, Lim J, Dale JM, Chen H, Shinn P, Palm CJ, Southwick AM, Wu HC, Kim C, Nguyen M: Empirical analysis of transcriptional activity in the Arabidopsis genome. Science. 2003, 302: 842-846. 10.1126/science.1088305.PubMedView ArticleGoogle Scholar
- Umezawa T, Sakurai T, Totoki Y, Toyoda A, Seki M, Ishiwata A, Akiyama K, Kurotani A, Yoshida T, Mochida K: Sequencing and analysis of approximately 40,000 soybean cDNA clones from a full-length-enriched cDNA library. DNA Res. 2008, 15: 333-346. 10.1093/dnares/dsn024.PubMed CentralPubMedView ArticleGoogle Scholar
- Alexandrov NN, Troukhan ME, Brover VV, Tatarinova T, Flavell RB, Feldmann KA: Features of Arabidopsis genes and genome discovered using full-length cDNAs. Plant Mol Biol. 2006, 60: 69-85. 10.1007/s11103-005-2564-9.PubMedView ArticleGoogle Scholar
- Ralph SG, Chun HJ, Cooper D, Kirkpatrick R, Kolosova N, Gunter L, Tuskan GA, Douglas CJ, Holt RA, Jones SJ: Analysis of 4,664 high-quality sequence-finished poplar full-length cDNA clones and their utility for the discovery of genes responding to insect feeding. BMC Genomics. 2008, 9: 57-10.1186/1471-2164-9-57.PubMed CentralPubMedView ArticleGoogle Scholar
- Jia J, Fu J, Zheng J, Zhou X, Huai J, Wang J, Wang M, Zhang Y, Chen X, Zhang J: Annotation and expression profile analysis of 2073 full-length cDNAs from stress-induced maize (Zea mays L.) seedlings. Plant J. 2006, 48: 710-727. 10.1111/j.1365-313X.2006.02905.x.PubMedView ArticleGoogle Scholar
- Alexandrov NN, Brover VV, Freidin S, Troukhan ME, Tatarinova TV, Zhang H, Swaller TJ, Lu YP, Bouck J, Flavell RB: Insights into corn genes derived from large-scale cDNA sequencing. Plant Mol Biol. 2009, 69: 179-194. 10.1007/s11103-008-9415-4.PubMed CentralPubMedView ArticleGoogle Scholar
- Mignone F, Gissi C, Liuni S, Pesole G: Untranslated regions of mRNAs. Genome Biol. 2002, 3: REVIEWS0004-10.1186/gb-2002-3-3-reviews0004.PubMed CentralPubMedView ArticleGoogle Scholar
- Iida K, Seki M, Sakurai T, Satou M, Akiyama K, Toyoda T, Konagaya A, Shinozaki K: Genome-wide analysis of alternative pre-mRNA splicing in Arabidopsis thaliana based on full-length cDNA sequences. Nucleic Acids Res. 2004, 32: 5096-5103. 10.1093/nar/gkh845.PubMed CentralPubMedView ArticleGoogle Scholar
- Kikuchi S, Satoh K, Nagata T, Kawagashira N, Doi K, Kishimoto N, Yazaki J, Ishikawa M, Yamada H, Ooka H: Collection, mapping, and annotation of over 28,000 cDNA clones from japonica rice. Science. 2003, 301: 376-379. 10.1126/science.1081288.PubMedView ArticleGoogle Scholar
- Ner-Gaon H, Leviatan N, Rubin E, Fluhr R: Comparative cross-species alternative splicing in plants. Plant Physiol. 2007, 144: 1632-1641. 10.1104/pp.107.098640.PubMed CentralPubMedView ArticleGoogle Scholar
- Davuluri RV, Sun H, Palaniswamy SK, Matthews N, Molina C, Kurtz M, Grotewold E: AGRIS: Arabidopsis gene regulatory information server, an information resource of Arabidopsis cis-regulatory elements and transcription factors. BMC Bioinformatics. 2003, 4: 25-10.1186/1471-2105-4-25.PubMed CentralPubMedView ArticleGoogle Scholar
- Riechmann JL, Heard J, Martin G, Reuber L, Jiang C, Keddie J, Adam L, Pineda O, Ratcliffe OJ, Samaha RR: Arabidopsis transcription factors: genome-wide comparative analysis among eukaryotes. Science. 2000, 290: 2105-2110. 10.1126/science.290.5499.2105.PubMedView ArticleGoogle Scholar
- Takahashi H, Shimizu A, Arie T, Rosmalawati S, Fukushima S, Kikuchi M, Hikichi Y, Kanda A, Takahashi A, Kiba A: Catalog of Micro-Tom tomato responses to common fungal, bacterial, and viral pathogens. J Gen Plant Pathol. 2005, 71: 8-22. 10.1007/s10327-004-0168-x.View ArticleGoogle Scholar
- Tsugane T, Watanabe M, Yano K, Sakurai N, Suzuki H, Shibata D: Expressed sequence tags of full-length cDNA clones from the miniature tomato (Lycopersicon esculentum) cultivar Micro-Tom. Plant Biotechnol. 2005, 22: 161-165.View ArticleGoogle Scholar
- Iijima Y, Nakamura Y, Ogata Y, Tanaka K, Sakurai N, Suda K, Suzuki T, Suzuki H, Okazaki K, Kitayama M: Metabolite annotations based on the integration of mass spectral information. Plant J. 2008, 54: 949-962. 10.1111/j.1365-313X.2008.03434.x.PubMed CentralPubMedView ArticleGoogle Scholar
- Chomczynski P, Sacchi N: Single-step method of RNA isolation by acid guanidinium thiocyanate-phenol-chloroform extraction. Anal Biochem. 1987, 162: 156-159. 10.1016/0003-2697(87)90021-2.PubMedView ArticleGoogle Scholar
- Ewing B, Hillier L, Wendl MC, Green P: Base-calling of automated sequencer traces using phred. I. Accuracy assessment. Genome Res. 1998, 8: 175-185.PubMedView ArticleGoogle Scholar
- Green P: Phrap documentation. 1996Google Scholar
- Rozen S, Skaletsky H: Primer3 on the WWW for general users and for biologist programmers. Methods Mol Biol. 2000, 132: 365-386.PubMedGoogle Scholar
- Mott R: EST_GENOME: a program to align spliced DNA sequences to unspliced genomic DNA. Comput Appl Biosci. 1997, 13: 477-478.PubMedGoogle Scholar
- Quevillon E, Silventoinen V, Pillai S, Harte N, Mulder N, Apweiler R, Lopez R: InterProScan: protein domains identifier. Nucleic Acids Res. 2005, 33: W116-120. 10.1093/nar/gki442.PubMed CentralPubMedView ArticleGoogle Scholar
- Barrell D, Dimmer E, Huntley RP, Binns D, O'Donovan C, Apweiler R: The GOA database in 2009--an integrated Gene Ontology Annotation resource. Nucleic Acids Res. 2009, 37: D396-403. 10.1093/nar/gkn803.PubMed CentralPubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.