Research

Open Access

# ESTuber db: an online database for *Tuber borchii* EST sequences

Barbara Lazzari*†¹,³, Andrea Caprera†²,³, Cristian Cosentino¹,⁵,
Alessandra Stella³, Luciano Milanesi⁴ and Angelo Viotti¹

Address: ¹Istituto di Biologia e Biotecnologia Agraria, via Bassini 15, 20133 Milan, Italy, ²CISI, Via Fratelli Cervi 93, 20090 Segrate (MI), Italy, ³Parco Tecnologico Padano, Via Einstein – Località Cascina Codazza, 26900 Lodi, Italy, ⁴Istituto Tecnologie Biomediche, Via Fratelli Cervi 93, 20090 Segrate (MI), Italy and ⁵Darmstadt University of Technology, Institute of Botany, Schnittspahnstrasse 3-5, 64287 Darmstadt, Germany

Email: Barbara Lazzari* - lazzari@ibba.cnr.it; Andrea Caprera - andrea.caprera@itb.cnr.it; Cristian Cosentino - cosentino@bio.tu-darmstadt.de; Alessandra Stella - alessandra.stella@tecnoparco.org; Luciano Milanesi - luciano.milanesi@itb.cnr.it; Angelo Viotti - viotti@ibba.cnr.it

* Corresponding author †Equal contributors

from Italian Society of Bioinformatics (BITS): Annual Meeting 2006
Bologna, Italy. 28–29 April, 2006

Published: 8 March 2007

BMC Bioinformatics 2007, 8(Suppl 1):S13 doi:10.1186/1471-2105-8-S1-S13

This article is available from: http://www.biomedcentral.com/1471-2105/8/S1/S13

© 2007 Lazzari 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:** The ESTuber database (http://www.itb.cnr.it/estuber) includes 3,271 *Tuber borchii* expressed sequence tags (EST). The dataset consists of 2,389 sequences from an in-house prepared cDNA library from truffle vegetative hyphae, and 882 sequences downloaded from GenBank and representing four libraries from white truffle mycelia and ascocarps at different developmental stages. An automated pipeline was prepared to process EST sequences using public software integrated by in-house developed Perl scripts. Data were collected in a MySQL database, which can be queried via a php-based web interface.

**Results:** Sequences included in the ESTuber db were clustered and annotated against three databases: the GenBank nr database, the UniProtKB database and a third in-house prepared database of fungi genomic sequences. An algorithm was implemented to infer statistical classification among Gene Ontology categories from the ontology occurrences deduced from the annotation procedure against the UniProtKB database. Ontologies were also deduced from the annotation of more than 130,000 EST sequences from five filamentous fungi, for intra-species comparison purposes.

Further analyses were performed on the ESTuber db dataset, including tandem repeats search and comparison of the putative protein dataset inferred from the EST sequences to the PROSITE database for protein patterns identification. All the analyses were performed both on the complete sequence dataset and on the contig consensus sequences generated by the EST assembly procedure.

**Conclusion:** The resulting web site is a resource of data and links related to truffle expressed genes. The Sequence Report and Contig Report pages are the web interface core structures which, together with the Text search utility and the Blast utility, allow easy access to the data stored in the database.
**Background**

Despite being dependent on other living organisms for energy supply, fungi have colonized every biotic or abiotic environment. Deciphering their genetic information is important *per se* but also provides integrative knowledge on the genetic interactions among fungi and other organisms. In this work the hypogeous symbiotic white truffle *Tuber borchii* (*Tb*) was considered. This species has recently attracted attention for its *in vitro* manipulation properties [1], as a possible resource of antiviral agents [2] and for other interesting features [3,4].

A cDNA library was prepared from *in vitro* cultured mycelia and extensive EST sequencing was performed. A total of 2,389 successful sequences was obtained and submitted to the EBI databank for GenBank accession number assignment. Sequences were given the accession numbers from AM165128 to AM167516. Additional 882 sequences from *Tb* vegetative hyphae and ascocarps were downloaded from GenBank and added to the dataset, thus obtaining a collection of 3,271 *Tb* ESTs extracted from five libraries. The aim was to produce an extensive and easily accessible database for truffle, providing sequence analysis and annotation, as well as links to related resources. To date, the ESTuber db [5] is the most complete information repository for truffle EST sequences. Together with the truffle.org project databases [6], ESTuber db can offer information of wide interest to truffle investigators.

**Construction and content**

**Library preparation and clone management**

*Tb* vegetative hyphae were grown on PDA agar medium for 20 days, collected and lyophylized. Total RNA was extracted by lithium chloride differential precipitation and poly(A)⁺ RNA was purified on a cellulose oligo-dT column. cDNA was synthesized using the Stratagene's (La Jolla, CA) cDNA synthesis and cloning kit, directionally ligated into the Uni-zap XR (λ zap II) vector (site 1, EcoRI; site 2, XhoI) and packaged using Gigapack Gold packaging extracts. A library (named 'T' in the ESTuber db web interface) was obtained representing more than 7 × 10⁷ plaque forming units. Plasmids (pBluescript SK) containing cDNA inserts were mass-excised from phage stocks using ExAssit Interference Resistant Helper Phage and propagated in SOLR cells, according to the manufacturer's instructions.

More than 12,000 plasmid clones were selected and replicated twice in 384 well microplates by the colony picking QPix workstation (Genetix). 3,068 clones were subsequently amplified in 96 well microplates by the TemplPhi Amplification kit (Amersham), then lyophylized and submitted to Macrogen Inc. (Seoul, Korea) for direct 3' sequencing. Sequencing was performed with the

BigDye Terminator method (Applied Biosystems) on capillary sequencers, using the M13 reverse primer.

**Sequence introduction in the database**

From sequencing of the T library, a total of 3,068 sequences was obtained. Electropherograms were read with the Phred program [7] and multifasta sequence and quality files were created and used to feed the program Lucy [8]. Both low quality regions and vector sequences were removed, according to Lucy standard parameters, thus obtaining 2,389 successful sequences longer than 100 bp. As 3' sequencing was performed, vector-free, high quality sequences were reverse-complemented with the program REVSEQ [9].

882 *Tb* sequences from two non-normalised ascocarp libraries (named 'MF' and 'IF' in the web interface, obtained from mature and immature ascocarps, respectively) [10] and from two non-normalised mycelium libraries (named 'SM' and 'LM', obtained from vegetative hyphae grown on solid or liquid medium, respectively) [11], were downloaded from GenBank and added to the dataset. Full details on the five libraries are included in the Library details page of the ESTuber web site. All the sequence files, the corresponding GenBank accession numbers and the available quality files and electropherograms were stored in the MySQL database.

**Sequence processing and annotation pipeline**

An automated pipeline (Fig. 1) was prepared to process EST sequences using public software integrated by a number of in-house developed Perl scripts. The complete dataset in multifasta format was used as input for the CAP3 program [12]: parameters were set to -p 98 -o 100 -f 30, for appropriate EST clustering. The CAP3 output was parsed and relevant data were stored in the database. A unigene dataset was defined, including all the singlets and each contig's longest sequence.

The 3,271 EST sequences and all the contig consensus sequences were annotated locally with a triple procedure.

The main annotation was performed with blastx [13] versus the GenBank nr database (referred to as 'NCBI blast' in the web interface). The GenBank nr database encompasses all non-redundant GenBank CDS translations + RefSeq Proteins + PDB + SwissProt + PIR + PRF for a total of more than 2,450,000 sequences, at present. All possible frames were considered for EST translation. Low complexity filters were disabled, while the remaining blastx parameters were the defaults and no threshold was set. The BLAST output was parsed with an in-house prepared parser. The most important BLAST output values, as well as the complete BLAST output pages, were stored in the database.
<figure><img src="image_1.png" /><figcaption>**Figure 1**
The ESTuber db pipeline. Data flow among the main programs included in the ESTuber db pipeline. Accessory programs were excluded from the chart.</figcaption></figure>

The second annotation (referred to as 'genomic blast' in the web interface) was performed against an in-house prepared database encompassing 42,029 fungal genomic sequences from four filamentous fungi (Aspergillus ssp, Fusarium ssp, Magnaporthe ssp, Neurospora ssp) and a dimorphic fungus (Saccharomyces ssp), downloaded from the GenBank CoreNucleotide db (January 2006). Complete chromosome sequences as well as whole genome shotgun sequences were excluded from the database. The blastn algorithm was used, low complexity sequence filtering was performed and the E-value threshold of 1e-10 was set.

The third annotation was performed with blastx versus the UniProtKB database [14] – encompassing at present more than 2,700,000 entries – that is the reference database for the Gene Ontology (GO) project [15] (referred to as 'GO blast' in the web interface). The E-value threshold of 1e-10 was adopted.

Based on data contained in the available database of associations among proteins and GO elements [16], a Perl

script was developed to relate ESTuber GO blast best blast hits to Gene Ontology categories. All the GO elements related to each best blast hit protein were considered, even those having indirect relations. In these cases, the GO 'is_a' relations were taken into account. Matching ontologies were stored into the database, to allow successive dynamic creation of statistics for the ontology occurrences. 130,459 EST sequences from five fungi (the same organisms considered for the preparation of the genomic blast database) were downloaded from GenBank and processed for GO annotation to provide intra-species comparisons of the ontology distribution in GO categories.

The annotation procedures were carried out either on a Linux 9-biprocessor AMDx64 cluster with 2 GB RAM on each node with the mpiBLAST program [17], that allows parallelization of the BLAST procedure, or through the BioinfoGRID [18] computing resource, that uses a series of calculation sites to perform very time-consuming elaborations in a distributed way, parallelizing their executions and exploiting the use of remote resources efficiently.

The entire EST dataset, as well as the contig consensus sequences, were processed for tandem repeat search with the Tandem Repeats Finder program (parameters: 2, 7, 7, 50, 2000) [19]. Prior to repeat search PolyA and PolyT tail containing sequences were trimmed with the TrimEST program (EMBOSS) (parameters: -minlength 4 -mis-matches 3). Links to the Tandem Repeats Finder output files were introduced in the web interface, when present.

The ESTate package programs [20] were used to infer putative protein sequences from the whole ESTtuber dataset and from the contig consensus sequences. Fasta2usage (parameters: -w 6 -j 3) and wordprob were used to prepare specific word probability files and FrameFinder was used for sequence translation and selection of the most probable putative protein sequences. During the protein prediction procedure, comparisons were made among protein sets obtained with the same parameters but with different word probability files. As a very limited number of Tuber coding sequences was present in public databases, there was no possibility to generate a truffle specific word probability file. Two different in-house generated files were tested: the former obtained from 111,936 fungi mRNA sequences downloaded from the NCBI RefSeq database, and the latter specific for filamentous fungi, obtained from 41,190 mRNA sequences from Magnaporthe ssp, Aspergillus ssp and Neurospora ssp, downloaded from RefSeq. No Tuber ssp and Fusarium ssp mRNA sequences were present in RefSeq to date. Comparisons on the derivative protein datasets were based on sequence length and on the number of sequences finding their best translation on
a lower frame, with respect to the expected composition of
forward/reverse oriented sequences of the dataset. When
considering these parameters, the best results were
obtained with the first of the above-mentioned word
probability files. Furthermore, the number of PROSITE
matches found for the protein set created according to the
fungi specific word probability file was 10% higher than
the one retrieved when translations generated according
to the filamentous fungi specific file were used. A supple-
mentary control on the automatic protein prediction pro-
cedure was carried out by BLAST. Contig consensus
sequences were annotated by blastx against the GenBank
nr database and six frames translations were considered.
Translations of the contig consensus sequences obtained
by FrameFinder were also annotated by blastp against the
same database and results were compared. Differences
were observed for 15 sequences out of 356 (4.2%), that
were annotated on a different translation frame with
respect to the one selected by FrameFinder. Based on these
observations, the FrameFinder algorithm was considered
suitable for this work purposes, and protein sequences
predicted according to the fungi word probability file
were inserted in the database and used as queries against
PROSITE [21].

Comparisons to the PROSITE database – version as of 24/01/2006 – were performed locally with the ScanPROSITE program [22] and matching patterns were retrieved and linked to the corresponding sequences in the web interface. Frequently matching patterns and profiles were excluded from the analysis.

Further details on software usage and parameters setting
are provided at the Processing, Assembly and Annotation
Protocol page of the cited web site.

Statistics on the sequence analysis and on the database
current status were inferred from data stored in the data-
base and are reported in Table 1.

The database

The ESTuber MySQL database is a dynamic structure
where new sequences will be added and new features are
being implemented, so the database structure is subject to
modifications.

The 'Sequences' table and the 'Contigs' table are the core
structures of the database, together with the
'Cont_Seq_relations' table, where relations among
sequences and contigs are defined (Fig. 2).

The web interface

The ESTuber db web interface is based on the php lan-
guage and manages all the incoming queries as well as all
the graphical outputs dynamic creation.

Contig graphical display and GO statistics graphical bars
are prepared on-the-fly in response to the users' requests
and are not stored in any part of the database (Fig. 3). The
dynamic management of these data provides up-to-date
displays without the need for refreshing corresponding
database tables.

Utility and discussion

Data flow in the ESTuber db pipeline

Sequence introduction in the ESTuber db was kept apart
from the main automated pipeline as chromatograms or
sequence quality files were not always available. EST
sequences downloaded from public repositories lack all
the upstream information and were introduced as such.
Nonetheless, when possible, quality data were stored in
the database for possible further analysis with programs
that require these data as mandatory.

The main pipeline is fed with a single multifasta file con-
taining all the available EST sequences; a number of in-
house developed Perl scripts allows automatic data flow
among the main programs ensuring full compatibility
among all the programs' input/output files. Program out-
puts are parsed and all relevant data are stored in the data-
base, so that at the end of the pipeline the database tables
are completely filled in. The pipeline structure is modular:
each main program can be added or removed from the
analysis easily and quickly. This allowed the authors to
use this same structure for different purposes, with differ-
ent datasets ([23] and A. Caprera, B. Lazzari, A. Stella, P.
Mariani, unpublished results).

Data availability and retrieval

The ESTuber db was structured to be a repository of data
obtained from EST sequence analysis and to provide links
to related external resources. Three main interfaces were
set up to access the database contents: the Sequence/Con-
tig report pages, the Text search utility and the local Blast
utility.

The Sequence report and Contig report pages are the start-
ing points to access sequence and contig related data. Data
are presented in summary tables and are linked to internal
and external sources. Furthermore, each sequence and
each contig were assigned a detailed page where the nucle-
otide sequence is given. In each contig page, the contig
graphical display is provided, as well as the contig align-
ment. Original blast outputs are retrievable for each
sequence and each contig consensus sequence.

The text search utility is structured to allow performing
combined keyword searches on all the fields of the
Sequence/Contig report tables. As NCBI blast results are
presented for each sequence, a user-defined significance
threshold based on the E-value can be set in the search to
<table><caption>Table 1: Statistics on the ESTuber db status. Sequences are considered annotated for Best Blast Hit E-values ≤ 1e-10.</caption><tr><td>Total number of sequences</td><td>3,271</td></tr><tr><td>Average base count</td><td>530.43</td></tr><tr><td>Number of singletons</td><td>2,071</td></tr><tr><td>Number of contigs</td><td>356</td></tr><tr><td>Number of putative unigenes</td><td>2,427</td></tr><tr><td>Sequences annotated versus GenBank nr db (%)</td><td>60.19</td></tr><tr><td>Sequences annotated versus GO db (%)</td><td>48.88</td></tr><tr><td>Sequences annotated versus genomic sequences db (%)</td><td>8.6</td></tr><tr><td>Repeats-containing sequences (%)</td><td>7.64</td></tr><tr><td>Sequences containing PROSITE matches (%)</td><td>16.84</td></tr></table>

retrieve only significantly annotated sequences. Searches can be restricted to sequence subsets (unigenes/not unigenes, singlets/contig related sequences or repeats containing sequences). Query outputs have the same format of the sequence/contig report tables, presenting only the entries matching the searched query, and can be downloaded both as multifasta sequence files and as text files with tab separated values, containing all the Sequence/Contig report page fields. This format can be easily imported in the most common spreadsheets.

The local blast interface allows users to blast their own sequences against either the nucleic ESTuber db dataset or the derivative putative protein database.

Additional information on the sequence dataset is provided by the two statistics web pages: the first reporting statistics on sequences and on the database status, and the second on the ontologies distribution in the database. GO blast annotations were used to infer the distribution of the annotated sequences among Gene Ontology functional categories, and statistics on the ontology occurrences were created for all the sequences of the dataset, for library-specific subsets and for the mycelium and ascocarps subsets. Additional statistics are provided concerning a dimorphic and four filamentous fungi EST datasets. Ontologies can be hierarchically browsed and are extensively linked to the AmiGO web pages [24].

Download of the ESTuber sequences and of the contig consensus sequences is allowed in multiple formats, and CAP3 alignments and .ace files are also available. Sequence user-defined subsets can be easily downloaded as query outputs.

The web interface includes a detailed help page to assist users in browsing the ESTuber db and in interpreting program outputs.

Sequence annotations are periodically updated and the consistency of external links is constantly verified. We report here on the second version of the ESTuber db: the previous release, including only the T library sequences, is maintained and is accessible via the link given in the ESTuber home page.

### Dataset analysis

The ESTuber db unigene dataset consists of 2,427 sequences. Based on sequence annotation versus the Gen-Bank nr db, 1,400 *Tuber borchii* unigenes were related to known sequences (BLAST E-value = 1e-10). Therefore, the remaining 1,427 unigenes are novel sequences. Among

<figure><img src="image_2.png" /><figcaption>**Figure 2**
The ESTuber database. Main tables of the ESTuber db.</figcaption></figure>
<figure><img src="image_3.png" /><figcaption>**Figure 3**
The ESTuber db contig page. A detailed page for each contig is included in the database, where a graphical display, the contig consensus sequence and the contig alignment are given, together with any other supplementary information.</figcaption></figure>

these, the 73.16 % belong to the in-house prepared mycelium cDNA library.

The assembly procedure generated 356 contigs. Analysis of their composition revealed that only 25 contigs contain both ascocarp and mycelium sequences, while the others are tissue-specific, 269 being composed only by mycelium sequences and 62 being composed only by ascocarp sequences. The most abundant transcripts in the database are assembled in four contigs, two composed by ascocarp sequences (Contig 86, 35 sequences, and Contig 41, 26 sequences) and two by mycelium sequences (Contig 55, 28 sequences, and Contig 210, 25 sequences). None of these contigs is significantly annotated, thus confirming the specificity of these sequences for truffle.

### Tandem repeats analysis

Tandem Repeats Finder identified 250 repeats-containing sequences in the complete EST dataset as well as 32 repeats-containing sequences in contig consensus sequences. Due to the fact that a number of EST sequences contained more than one repeat, 331 repeats were actually identified in the complete dataset. The same analysis was performed on a dataset including all the singlet sequences and the contig consensus sequences, and 33 repeats longer than 24 basepairs were identified. Blastn analysis of these repeats against the nr nucleic GenBank database (parameters: -w 7 -e 1000) did not identify any significant homology with regions associated to known features. The repeats consensus sequences were translated with Transeq (EMBOSS) and the outcome aminoacid sequences were used to scan the PROSITE database. Three known patterns were identified: the ubiquitin domain signature and pro-

file (PROSITE entries PS00299 and PS50053, found in Contig22), the bipartite nuclear targeting sequence (PROSITE entry PS00015, found in sequence T15O19) and the ankyrin repeat region circular profile (PROSITE entry PS50297, found in sequences T12B4 and SMVA56).

The Tandem Repeats Finder algorithm was conceived to identify tandem repeats, but due to the fact that short sequence repetitions are considered as contiguous repeats with coinciding period size and repeat size, the algorithm identifies simple sequence repeats (SSR), too. In the ESTuber db sequence collection, the repeats consensus sequence length varied from 1 basepair (single base stretches), to 230 basepairs (in the polyubiquitin sequence T10J12), and 51 repeats (15.4 % of the total) were SSR, with period sizes shorter than 5 basepairs.

### Protein prediction and protein domains search

When PROSITE was queried with the ESTuber putative protein dataset, 551 positive hits were recovered, while 111 hits were reported when PROSITE was queried with contig consensus sequences translations. All the matching patterns were introduced in the ESTuber db and links to the correspondent PROSITE entries were added.

### Ontologies distribution analysis

Statistics on the Tuber ontology occurrences are provided for the whole sequence dataset, for library-specific subsets and for the ascocarps and mycelia subsets. Data from two libraries (T and SM), both obtained from hyphae grown on solid medium, are also merged in a further subset.

<figure><img src="image_4.png" /><figcaption>**Figure 4**
Discovery rates for the ascocarp and mycelium datasets. Organ unique sequences plotted against the number of ESTs. The discovery rate is represented by the angular coefficient of the line connecting two consecutive points. The straight line ($y=x$) corresponds to a discovery rate of 1.</figcaption></figure>
<table><caption>Table 2: Ontologies distribution in GO categories.</caption><thead><tr><th></th><th>Tb ascocarps</th><th>Tb mycelia</th><th>Aspergillus</th><th>Fusarium</th><th>Magnaporthe</th><th>Neurospora</th><th>Saccharomyces</th></tr></thead><tbody><tr><td>Total number of sequences</td><td>683</td><td>2,588</td><td>35,924</td><td>10,302</td><td>53,102</td><td>28,089</td><td>3,041</td></tr><tr><td>Unigenes percentage</td><td>62.22</td><td>84.89</td><td>91.6</td><td>91.67</td><td>94.78</td><td>56.28</td><td>99.44</td></tr><tr><td>GO annotation percentage</td><td>55.49</td><td>46.14</td><td>61.09</td><td>46.89</td><td>44.54</td><td>44.15</td><td>78.72</td></tr><tr><td><strong>MOLECULAR FUNCTION</strong></td><td>39.97</td><td>42.19</td><td>54.2</td><td>43.74</td><td>40.77</td><td>40.14</td><td>67.91</td></tr><tr><td>binding</td><td>20.94</td><td>21.14</td><td>25.41</td><td>19.78</td><td>18.84</td><td>19.5</td><td>43.47</td></tr><tr><td>catalytic activity</td><td>17.72</td><td>20.83</td><td>32.1</td><td>30.35</td><td>20.39</td><td>21.98</td><td>30.71</td></tr><tr><td>structural molecule activity</td><td>6.59</td><td>6.45</td><td>5.11</td><td>2.3</td><td>5.57</td><td>7.77</td><td>16.38</td></tr><tr><td>transporter activity</td><td>3.51</td><td>3.71</td><td>6.43</td><td>6.57</td><td>5.76</td><td>3.9</td><td>6.91</td></tr><tr><td>molecular function unknown</td><td>3.37</td><td>0.97</td><td>0.88</td><td>0.62</td><td>0.72</td><td>0.54</td><td>1.84</td></tr><tr><td>signal transducer activity</td><td>0.15</td><td>0.39</td><td>1.05</td><td>0.49</td><td>0.6</td><td>0.61</td><td>0.72</td></tr><tr><td>translation regulator activity</td><td>1.02</td><td>1.43</td><td>1.2</td><td>1.48</td><td>1.82</td><td>0.95</td><td>1.51</td></tr><tr><td>transcription regulator activity</td><td>0</td><td>0.46</td><td>1.79</td><td>0.79</td><td>2.39</td><td>1.08</td><td>2.7</td></tr><tr><td>enzyme regulator activity</td><td>0.15</td><td>0.7</td><td>0.39</td><td>0.19</td><td>0.34</td><td>0.2</td><td>1.71</td></tr><tr><td>antioxidant activity</td><td>0.73</td><td>0.35</td><td>0.62</td><td>0.14</td><td>0.55</td><td>0.48</td><td>1.25</td></tr><tr><td><strong>BIOLOGICAL PROCESS</strong></td><td>42.02</td><td>35.55</td><td>43.59</td><td>34.49</td><td>32.45</td><td>34.18</td><td>68.33</td></tr><tr><td>physiological process</td><td>41.43</td><td>32.92</td><td>40.93</td><td>32.93</td><td>28.91</td><td>32.4</td><td>67.41</td></tr><tr><td>cellular process</td><td>37.92</td><td>30.87</td><td>35.44</td><td>30.36</td><td>25.61</td><td>28.93</td><td>66</td></tr><tr><td>regulation of biological process</td><td>0.73</td><td>2.09</td><td>3.7</td><td>1.96</td><td>4.04</td><td>2.21</td><td>6.68</td></tr><tr><td><strong>CELLULAR COMPONENT</strong></td><td>15.96</td><td>15.92</td><td>15.56</td><td>8.83</td><td>11.89</td><td>12.84</td><td>55.87</td></tr><tr><td>organelle</td><td>13.76</td><td>11.67</td><td>11.81</td><td>6.88</td><td>8.97</td><td>10.5</td><td>50.51</td></tr><tr><td>protein complex</td><td>13.91</td><td>11.01</td><td>8.79</td><td>4.46</td><td>7.11</td><td>7.73</td><td>28.74</td></tr></tbody></table>

Ontologies distribution in the three main GO categories and in their subclasses of the Tb ascocarps and mycelia subsets and of the five fungi EST datasets considered for comparison purposes. Data on EST collection redundancies and on the GO annotation percentages are given at the top of the table. Matches with GO entries are expressed as percentage of matching sequences with respect to the total number of sequences of the considered dataset. Only GO categories where at least one of the considered organisms had a percentage of matches > 1 were included in the table.

Although single library contribution to the whole dataset is variable, the T library being far more represented than the others, the GO annotation percentage varies from

57.68 in the MF library to 46.84 in the T library. The NCBI annotation results reflect the same trend, the ascocarp libraries being more annotated than the mycelium ones and spanning from 70.33% of annotation in the MF library to 57.47% of annotation in the T library (data available on the ESTuber web site). Furthermore, ascocarp libraries are more redundant than mycelium libraries, as evident from the unigene percentage which is much lower for ascocarp sequences (62,22%) than for mycelium sequences (84,89%). This is confirmed by the analysis of the discovery rates obtained for the pools of ascocarp and mycelium sequences (Fig. 4).

Changes in the distribution of the ontology occurrences between sequences from ascocarps and from cultured hyphae are evident mainly in GO categories related to biological processes (Table 2). These differences are also present in comparisons between sequences from ascocarps at different maturation stages, the 'biological process' GO category being over-represented in mature ascocarps with respect to immature ascocarps (data available in the 'Statistics on GO annotation' page of the database web site).

<figure><img src="image_5.png" /><figcaption>Figure 5</figcaption></figure>

Pie representation of the ontology occurrences. Distribution of the matching ontologies in the second-level GO categories for the different considered organisms. Only GO categories where at least one of the organisms had a percentage of matches > 1 were included in pies.
Statistical analysis of the ontologies distribution is pro-
vided also for other five organisms (Tab. 2). Four of these
(Aspergillus, Fusarium, Magnaporthe and Neurospora) are
classified as filamentous fungi, while the fifth, Saccharo-
myces, has a dimorphic behaviour, exhibiting both yeast
forms and pseudohyphal phases. Dimorphism in S. cere-
visiae is associated with the expression of specific proteins
[25]. Saccharomyces EST sequences were considered in the
ontologies analysis, to investigate for a possible differen-
tial distribution of ontologies occurrences with respect to
the other monomorphic filamentous fungi.

With the aim to investigate the redundancy level of the
different EST collections, species-specific assembly proce-
dures were performed with CAP3 and the number of uni-
genes was determined for each dataset (Tab. 2). CAP3
parameters were the same as in the ESTuber pipeline.

Comparisons among different fungi EST collections reveal
a very high GO annotation percentage of *Saccharomyces*
sequences. Furthermore, the *Saccharomyces* dataset is
mostly composed of unigenes. Actually, *Aspergillus* and
*Saccharomyces* are more represented than the other fungi
in the UniProtKB database, while *Magnaporthe* is under-
represented. As a consequence, the high annotation per-
centage of *Saccharomyces* sequences could be almost in
part attributed to the comparison database composition.
Pie graphical representation of the ontology occurrences
in the different organisms provides a presentation of the
data that is independent from the total annotation per-
centage (Fig. 5). Accordingly, differences in representation
of specific GO categories are evident among the presented
organisms.

Conclusion

The ESTuber db version as of January 2006 encompasses
3,271 EST sequences from ascocarps and from *in-vitro* cul-
tured hyphae. At the moment none of the represented
libraries is specifically oriented to the study of the proc-
esses that regulate the establishment of symbiosis
between the host plant and the fungus. This is actually the
most interesting aspect of the truffle project itself. As truf-
fle is a typical niche product the interest in its biology is
concentrated in the few producing countries, but a deep
insight in the mechanisms involved in the symbiosis
establishment could be of wider interest. The efforts of the
truffle research community will, hopefully, lead to the
production of new sequences and other related data. Even
if at the moment the limited number of sequences doesn't
allow any other significant analysis, the ESTuber db is
structured to be easily integrated with new functionalities.
New tools will be added in response to the future needs of
the project.

**Availability and requirements**

The ESTuber db is available at http://www.itb.cnr.it/estu-
ber.

Authors' contributions

BL prepared the T library, defined the pipeline structure
and parameters and drafted the manuscript. AC structured
the database and the web interface and wrote all the pipe-
line accessory programs. CC cared about clone and
sequence management. AS participated in the design of
the study and critically revised the manuscript. LM coordi-
nated the integration of bioinformatical resources. AV
guided and coordinated the execution of the project. All
authors read and approved the final manuscript.

Acknowledgements

We acknowledge Prof. Paola Bonfante and Prof. Francis Martin for allowing us to introduce their sequences in the ESTuber db; Dr. Ivan Merelli and Dr.Salvatore Di Nardo for helping us in establishing the multi-processor mpiBLAST procedure; Prof. Francesco Salamini and Dr. Riccardo Percudani for critical reading of the manuscript, the Max Planck Institute for collaboration in library plates preparation and Dr. Carlo Pozzi for cooperation in structuring the database.

This work was supported by: CNR Strategic Project "TUBER: biotecnologia della micorrizzazione"; Parco Tecnologico Padano; Laboratory of Interdisciplinary Technologies in Bioinformatics (LITBIO) and CISI; MIUR "Functional Genomics", "Bioinformatics for Genome and Proteome" projects; EU BioinfoGRID project and FIRB projects.

This article has been published as part of *BMC Bioinformatics* Volume 8, Supplement 1, 2007: Italian Society of Bioinformatics (BITS): Annual Meeting 2006. The full contents of the supplement are available online at [http://www.biomedcentral.com/1471-2105/8?issue=S1](http://www.biomedcentral.com/1471-2105/8?issue=S1).

**References**

1. Grimaldi B, de Raaf MA, Filetici P, Ottonello S, Ballario P: **Agrobacterium-mediated gene transfer and enhanced green fluorescent protein visualization in the mycorrhizal ascomycete Tuber borchii: a first step towards truffle genetics. Curr Genet** 2005, 48:69-74.

2. Percudani R, Montanini B, Ottonello S: **The anti-HIV cyanovirin-N domain is evolutionarily conserved and occurs as a protein module in eukaryotes. Proteins** 2005, 60:670-678.

3. Montanini B, Viscomi AR, Bolchi A, Martin Y, Siverio JM, Balestrini R, Bonfante P, Ottonello S: **Functional properties and differential mode of regulation of the nitrate transporter from a plant symbiotic ascomycete. Biochem J** 2006, 15:125-134.

4. Barbieri E, Bertini L, Rossi I, Ceccaroli P, Saltarelli R, Guidi C, Zam-bonelli A, Stocchi V: **New evidence for bacterial diversity in the ascoma of the ectomycorrhizal fungus Tuber borchii Vittad. FEMS Microbiol Lett** 2005, 247:23-35.

5. ESTuber database [http://www.itb.cnr.it/estuber]

6. Truffle.org [http://www.truffle.org/](http://www.truffle.org/)

7. Ewing B, Hiller L, Wendl M, Green P: **Basecalling of automated sequence traces using phred. I. Accuracy assessment. Genome Research** 1998, 8:175-185.

8. Chou H-H, Holmes MH: **DNA sequence quality trimming and vector removal. Bioinformatics** 2001, 17:1093-1104.

9. **EMBOSS** [http://emboss.sourceforge.net/apps/](http://emboss.sourceforge.net/apps/)

10. Gabella S, Abba S, Duplessis S, Montanini B, Martin F, Bonfante P: **Transcript profiling reveals novel marker genes involved in fruiting body formation in Tuber borchii. Eukaryot Cell** 2005, 4:1599-1602.

11. Lacourt I, Duplessis S, Abba S, Bonfante P, Martin F: **Isolation and characterization of differentially expressed genes in the**
mycelium and fruit body of Tuber borchii. *Appl Environ Microbiol* 2002, 68:4574-4582.

12. Huan X, Madan A: **CAP3:** A DNA sequence assembly program. *Genome Research* 1999, 9:868-877.

13. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: **Basic local alignment search tool**. *J Mol Biol* 1990, **215**:403-410.

14. UniProt [http://www.ebi.ac.uk/uniprot/]

15. The Gene Ontology Consortium: **Gene Ontology: tool for the unification of biology.** *Nature Genet* 2000, **25**:25-29.

16. GOA project [http://www.ebi.ac.uk/GOA]

17. mpiBLAST [http://mpiblast.lanl.gov/]

18. BioinfoGRID [http://www.bioinfogrid.eu]

19. Benson G: **Tandem repeats finder: a program to analyze DNA sequences**. *Nucleic Acid Res* 1999, **27**:573-580.

20. **ESTate - Expressed Sequence Tag Analysis Tools** [http://bioweb.pasteur.fr/docs/man/man/ESTate.I.html]

21. Falquet L, Pagni M, Bucher P, Hulo N, Sigrist CJ, Hofmann K, Bairoch A: **The PROSITE database, its status in 2002**. *Nucleic Acids Res* 2002, **30**:235-238.

22. Gattiker A, Gasteiger E, Bairoch A: **ScanPROSITE: a reference implementation of a PROSITE scanning tool**. *Appl Bioinformatics* 2002, **1**:107-108.

23. Lazzari B, Caprera A, Vecchietti A, Stella A, Milanesi L, Pozzi C: **ESTree db: a tool for peach functional genomics**. *BMC Bioinformatics* 2005, **6**(Suppl 4):S16.

24. AmiGO [http://www.godatabase.org/]

25. Viard B, Kuriyama H: **Phase-Specific Protein Expression in the Dimorphic Yeast Saccharomyces cerevisiae**. *Biochem Biophys Res Commun* 1997, **233**:480-486.

Publish with **BioMed Central** and every
scientist can read your work free of charge

"BioMed Central will be the most significant development for disseminating the results of biomedical research in our lifetime."

Sir Paul Nurse, Cancer Research UK

Your research papers will be:

* available free of charge to the entire biomedical community

* peer reviewed and published immediately upon acceptance

* cited in PubMed and archived on PubMed Central

* yours — you keep the copyright

Submit your manuscript here:
http://www.biomedcentral.com/info/publishing_adv.asp

<figure><img src="image_6.png" /></figure>

BioMedcentral