# Analysis of 13000 unique Citrus clusters associated with fruit quality, production and salinity tolerance

Javier Terol, Ana Conesa, Jose M Colmenero, Manuel Cercos, Francisco

Tadeo, Javier Agustí, Enriqueta Alós, Fernando Andrés, Guillermo Soler,
Javier Brumos, et al.

► To cite this version:

Javier Terol, Ana Conesa, Jose M Colmenero, Manuel Cercos, Francisco Tadeo, et al. Analysis of 13000 unique Citrus clusters associated with fruit quality, production and salinity tolerance. BMC Genomics, 2007, 8, 10.1186/1471-2164-8-31. doi:10.1186/1471-2164-8-31

HAL Id: hal-01605709

https://hal.science/hal-01605709

Submitted on 31 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
Research article

Open Access

# Analysis of 13000 unique *Citrus* clusters associated with fruit quality, production and salinity tolerance

Javier Terol¹, Ana Conesa¹, Jose M Colmenero¹, Manuel Cercos¹,
Francisco Tadeo¹, Javier Agustí¹, Enriqueta Alós¹, Fernando Andres¹,
Guillermo Soler¹, Javier Brumos¹, Domingo J Iglesias¹, Stefan Götz²,
Francisco Legaz¹, Xavier Argout³, Brigitte Courtois³, Patrick Ollitrault⁴,
Carole Dossat⁵, Patrick Wincker⁵, Raphael Morillon⁴ and Manuel Talon*¹

Address: ¹Centro de Genomica, Instituto Valenciano de Investigaciones Agrarias. Carretera Moncada - Náquera, Km. 4.5 Moncada (Valencia)
E46113, Spain, ²BET-ITACA, Universidad Politécnica de Valencia, Camino de Vera, s/n 46022 Valencia, Spain, ³CIRAD AMIS, UMR PIA - Avenue Agropolis - TA 40/03 34398 Montpellier Cedex 5, France, ⁴Genoscope-Centre National de Séquençage and Centre National de la Recherche Scientifique (CNRS) Unité Mixte de Recherche (UMR)-8030, 91000 Evry, France and ⁵CIRAD FLHOR, UPR "Amélioration génétique d'espèces à multiplication végétative", Avenue Agropolis - TA 40/03 34398 Montpellier Cedex 5, France

Email: Javier Terol - jterol@ivia.es; Ana Conesa - aconesa@ivia.es; Jose M Colmenero - chroma@ivia.es; Manuel Cercos - mcercos@ivia.es; Francisco Tadeo - ftadeo@ivia.es; Javier Agustí - jagusti@ivia.es; Enriqueta Alós - ealos@ivia.es; Fernando Andres - fandres@ivia.es; Guillermo Soler - gsoler@ivia.es; Javier Brumos - jbrumos@ivia.es; Domingo J Iglesias - iglesias@ivia.es; Stefan Götz - stefang@fis.upv.es; Francisco Legaz - flegaz@ivia.es; Xavier Argout - xavier.argout@cirad.fr; Brigitte Courtois - courtois@cirad.fr; Patrick Ollitrault - ollitrault@cirad.fr; Carole Dossat - cdossat@genoscope.cns.fr; Patrick Wincker - pwincker@genoscope.cns.fr; Raphael Morillon - morillon@cirad.fr; Manuel Talon* - mtalon@ivia.es

Published: 25 January 2007

Received: 19 July 2006

BMC Genomics 2007, 8:31 doi:10.1186/1471-2164-8-31

Accepted: 25 January 2007

This article is available from: http://www.biomedcentral.com/1471-2164/8/31

© 2007 Terol 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:** Improvement of *Citrus*, the most economically important fruit crop in the world, is extremely slow and inherently costly because of the long-term nature of tree breeding and an unusual combination of reproductive characteristics. Aside from disease resistance, major commercial traits in *Citrus* are improved fruit quality, higher yield and tolerance to environmental stresses, especially salinity.

**Results:** A normalized full length and 9 standard cDNA libraries were generated, representing particular treatments and tissues from selected varieties (*Citrus clementina* and *C. sinensis*) and rootstocks (*C. reshni*, and *C. sinensis* × *Poncirus trifoliata*) differing in fruit quality, resistance to abscission, and tolerance to salinity. The goal of this work was to provide a large expressed sequence tag (EST) collection enriched with transcripts related to these well appreciated agronomical traits. Towards this end, more than 54000 ESTs derived from these libraries were analyzed and annotated. Assembly of 52626 useful sequences generated 15664 putative transcription units distributed in 7120 contigs, and 8544 singletons. BLAST annotation produced significant hits for more than 80% of the hypothetical transcription units and suggested that 647 of these might be *Citrus* specific unigenes. The unigene set, composed of ~13000 putative different transcripts, including more than 5000 novel *Citrus* genes, was assigned with putative functions based on similarity, GO annotations and protein domains

<footnote>* Corresponding author</footnote>
**Conclusion:** Comparative genomics with Arabidopsis revealed the presence of putative conserved orthologs and single copy genes in *Citrus* and also the occurrence of both gene duplication events and increased number of genes for specific pathways. In addition, phylogenetic analysis performed on the ammonium transporter family and glycosyl transferase family 20 suggested the existence of *Citrus* paralogs. Analysis of the *Citrus* gene space showed that the most important metabolic pathways known to affect fruit quality were represented in the unigene set. Overall, the similarity analyses indicated that the sequences of the genes belonging to these varieties and rootstocks were essentially identical, suggesting that the differential behaviour of these species cannot be attributed to major sequence divergences. This *Citrus* EST assembly contributes both crucial information to discover genes of agronomical interest and tools for genetic and genomic analyses, such as the development of new markers and microarrays.

**Background**

Citrus fruits are the first fruit crop in international trade in terms of economic value (FAO, 2004). Citrus fruits are typically grown in 140 countries located in tropical and subtropical areas with "Mediterranean" type climates, often facing severe abiotic stresses such as salinity, drought and iron chlorosis. *Citrus* species also suffer from different diseases and pests that considerably affect tree growth and fruit crop. The survival of the citrus industry is today critically dependent on genetically superior cultivars. However, citrus improvement through traditional techniques is unfortunately very difficult due to the unusual combination of biological characteristics of *Citrus* species, their low genetic diversity and the long-term nature of tree breeding. Thus, *Citrus* show many biological characteristics such as gametophytic self- and cross-incompatibility, apomixy, juvenility, heterozygosis, dormancy, and surprising root/shoot interactions, that strongly hamper Citrus breeding. On the other hand, genetic and allelic diversity in *Citrus* cultivars is very scarce. The global linkage disequilibrium in the cultivated citrus that probably originated from three major taxa, may be the result of an initial allopatric evolution and further limitation for predominant apomixy [1]. The fact that only mutational and/or epigenetic events are apparently involved in the diversification of secondary species, combined with human selection, have strongly reduced global genetic diversity, restricting opportunity for genetic advance.

Genomics has provided new tools for crop improvement, helping to identify and select candidate genes responsible of agronomic characters of interest, and allowing the development of fast methods to incorporate these characters into crop plants. After the completion of the Arabidopsis genome sequence [2] and the publication of sequences of indica [3] and japonica [4] rice, plant researchers have been able to scan these genomes to identify and compare genes of interest. The completion of the poplar genome sequence [5] will supply a model for tree life forms.

EST sequencing projects have facilitated appropriate strategies for gene discovery [6], molecular markers identification [7,8], and many other functional genomic developments and tools, e.g. microarray approaches [9-11]. In *Citrus*, previous EST sequencing projects have released more than 130000 ESTs to Gen Bank, mainly from *C. sinensis*, *C. unshiu* and *C. clementina* [12-16]. This information has been used to develop two different microarray platforms, based on cDNA and short oligo sequences [12,17]. In this work, main efforts have been specifically focussed on the study of pivotal traits for *Citrus* breeding, such as fruit quality, productivity and salinity tolerance. In citrus there are many aspects of fruit quality such as fruit size, shape, colour, texture, flavour and aroma compounds, and several other organoleptic properties that are acquired during ripening and earlier stages of growth [18-20]. Regarding productivity, pivotal traits to be improved are the capacity for fruit set and the resistance to abscission. Clementine mandarin (*Citrus clementina*), for example, is a self-incompatible cultivar that shows elevated ovary and fruitlet abscission while sweet orange varieties (*Citrus sinensis*), that in general exhibit standard fruit-set ratios [21,22], may lose most of the yield during ripening. In addition, the quantity and quality of water can become a limiting factor to economically viable production. In *Citrus*, it is notorious for example, that during the periods of drought, leaves and fruits remain attached to the tree until water is available and soon afterwards these organs abscise [23,24]. It is also well known that salt excess affects the size and quality of the production. The capability of *Citrus* to tolerate salinity is mostly related to the ability of the rootstock to exclude chloride, although the nature of this mechanism remains unresolved [25]. Tolerant Cleopatra mandarin rootstock (*C. reshni*), for example, accumulated lower chloride amounts than sensitive Carrizo (*C. sinensis* × *Poncirus trifoliatata*) [26,27]. The scion variety also plays a role in salinity damage, and more tolerant varieties such as Clementine are generally preferred to sweet oranges [28]. Toward this goal, cDNA libraries were prepared from the following selected genotypes: the varieties *Citrus clementina* (cv.
Clementina de Nules; elevated fruit quality, low setting) and *C. sinensis* (cv Washington Navel and Navelina; pre-harvest abscission and low salt tolerance, respectively), and the rootstocks *C. reshni* (cv *Cleopatra*; salt tolerant), and *C. sinensis × Poncirus trifoliata* (cv Carrizo; salt sensitive). The information generated with this effort, complementary to the Spanish *Citrus* Functional Genomics Project [12], has also been used for the construction of a second generation cDNA microarray of recent release. In the study, special attention has been paid on the methodological aspects, in order to obtain accurate estimates of the number of different transcripts and precise predictions of their function.

The results of this joint initiative of the IVIA, CIRAD and Genoscope designed to provide new information and useful tools for *Citrus* improvement, will speed up the discovery of genes of major agronomic interest and facilitate the development of new markers and methods to rapidly identify improved genotypes.

## Results

### Library Construction

Samples were harvested from 4 different *Citrus* species: the varieties *Citrus clementina* (cv. Clementina de Nules; elevated fruit quality, low setting) and *C. sinensis* (cv Washington Navel and Navelina; pre-harvest abscission and low salt tolerance, respectively), and the rootstocks *C. reshni* (cv *Cleopatra*; salt tolerant), and *C. sinensis × Poncirus trifoliata* (cv Carrizo; salt sensitive).

EST sequences were derived from 9 cDNA libraries constructed with standard methods and 1 additional normalized full length cDNA library (Table 1). Standard libraries, prepared from all five cultivars, were designed to represent three main agronomic traits of interest for *Citrus* improvement, i.e. fruit quality, production and salt tolerance. Since fruit quality characteristics are acquired not only during ripening but also along fruit growth [18], 2 libraries covering all pivotal developmental stages were constructed from fruit pulp and peel tissues of the high quality variety, Clementine (FruitTF, and PhII-IIIVesicles1). Abscission zones and surrounding tissues from developing ovaries, fruitlets, leaves and ripe fruits were represented in 4 libraries (AbsAOv1, Abs Cov1, AbsDev, and AbsCFruit1) from Clementine, the genotype with impaired fruit set [29] and Washington Navel, the variety with higher pre-harvest abscission. To study the response to abiotic stresses, 3 libraries were prepared with samples derived from leaves and roots subjected to salinity or dehydration (LSH, KCl-Salt1, and EHR) from the salt tolerant rootstock, Cleopatra, the sensitive hybrid, Carrizo and the sensitive scion, Navelina [27,28].

The full length library was constructed in Clementine plants, the species of main interest, grown in open field and in controlled greenhouse conditions. A broad variety of tissues and organs at different developmental stages were harvested from healthy plants and from plants subjected to many biotic (viruses, insects, fungus...setc) and abiotic (drought, salinity, ozone, alkaline-calcareous soils, flooding.etc) stresses. This strategy was planned to obtain the widest representation of the *Citrus clementina* transcriptome, including low abundant cDNAs, and to facilitate identification of genes of agronomic interest. A detailed description of the libraries is given in Materials and Methods.

### EST Sequencing and Clustering

A total number of 54136 clones were single-pass sequenced from their 5' end. After low quality and vector trimming, 52626 sequences longer than 100 bp were assembled with the CAP3 program [30]. All these sequences are available at the EST section of the GenBank. The assembly grouped 44082 ESTs into 7120 contigs, while 8544 sequences remained as singletons. The combined set of contigs and singletons resulted in 15664 unigenes representing different putative transcripts from *Citrus* species. The average length of the unigenes was 1071 bp, and 9877 of them (63%) were longer than 1000 bp. The distribution of ESTs in the contigs was the expected, with many clusters composed of 10 or less ESTs (Fig 1A). To estimate redundancy levels, the 15664 unigenes were compared with each other with the BLASTN program [31]. Sequences with at least 98% nucleotide identity over a minimum of 100 bp were assumed to be derived from the same transcript and therefore were clustered in *supercontigs* using an in-house R algorithm. Clustering of the unigenes resulted in 1135 *supercontigs* and 12759 singletons, indicating ~26% redundancy. Thus, the real number of identified expressed genes was close to 13900.

The unigene consensus sequences were used as queries in a BLASTN search against a database including 130400 ESTs from *Citrus* species retrieved from GenBank. An e value of 1e-13 was used as a cut off to ensure that only almost identical sequences were detected. The results showed that more than 5159 unigenes (33%) did not produce a significant hit, indicating that these sequences had not previously been described in *Citrus* EST collections. However, the possibility that other parts of the same parental genes were present in these collections can not be discarded. Most novel sequences (4673) were derived from the normalized library and might represent transcripts expressed at very low levels. This idea was supported by the fact that 75% of these sequences were singletons. On the other hand, no major divergences or differences at the sequence level were observed between
<table><caption>Table 1: Summary of *Citrus* EST libraries</caption><thead><tr><th>Species</th><th>Library</th><th>Library description</th><th>Type</th><th>ESTs<sup>a</sup></th><th>Exclusive unigenes<sup>b</sup></th></tr></thead><tbody><tr><td rowspan="6">C. clementina</td><td>AbsAOvI</td><td>Abscission zone A and surrounding tissues of ethylene-treated ovary explants</td><td>Standard</td><td>900</td><td>229</td></tr><tr><td>AbsCOvI</td><td>Abscission zone C and surrounding tissues of ethylene-treated ovary explants</td><td>Standard</td><td>947</td><td>240</td></tr><tr><td>AbsDev</td><td>Laminar abscission zone and surrounding tissues (petiole and blade) of developing leaves</td><td>Standard</td><td>877</td><td>258</td></tr><tr><td>FruitTF</td><td>Flavedo and juice vesicles from fruits at different developmental stages</td><td>Standard</td><td>3999</td><td>1263</td></tr><tr><td>PhII-III/VesiclesI</td><td>Juice vesicles from fruits at phases II and III</td><td>Standard</td><td>1020</td><td>70</td></tr><tr><td>NFL</td><td>Organs and tissues at different stages. Plants under normal culture practices or subjected to biotic and abiotic stresses</td><td>Normalized Full length</td><td>41288</td><td>8440</td></tr><tr><td rowspan="3">C. sinensis</td><td>AbsCFruitI</td><td>Abscission zone C and surrounding tissues of mature fruits</td><td>Standard</td><td>783</td><td>218</td></tr><tr><td>LSH</td><td>Leaves from prolonged salt-treated plants</td><td>Standard</td><td>1009</td><td>340</td></tr><tr><td>KCl-SaltI</td><td>Roots subjected salinity (Cl-) treatments</td><td>Standard</td><td>960</td><td>253</td></tr><tr><td>C. reshni</td><td>EHR</td><td>Roots subjected to water stress and re-hydration</td><td>Standard</td><td>849</td><td>121</td></tr><tr><td colspan="2">Citrus sinensis × Poncirus trifoliata</td><td></td><td>Total</td><td>52632</td><td>11432</td></tr></tbody></table>

all *Citrus* species analyzed, including the five cultivars used in this study and those with sequences submitted to GenBank (*C. sinensis*, *Citrus × paradisi*, and *C. unshiu*, mainly).

The analysis of the contigs (Fig 1B) indicated that many of them were composed of ESTs from a single library (67%), while only 195 contigs (2%) included ESTs from 4 or more libraries. The unigene Contig6498 that displayed high similarity with ATPase-like proteins, for example, was found in 7 libraries, although it was composed of only 9 ESTs. The number of exclusive unigenes of a given cDNA library was obtained adding the number of singletons and the contigs composed of ESTs from this library (Table 1). It was determined that 11432 (72.9%) clusters were exclusive of a particular library. The standard cDNA libraries provided ~21% of the assembled ESTs while the number of exclusive clusters from these libraries represented more than 26% of the unigenes.

## Sequence Annotation

Annotation of the assembled sequences was initially based on primary sequence homology searches. A BLASTX search performed against the GenBank non redundant protein database [32] with an e value cut off ≤ 1e-10, yielded 13339 unigenes with significant hits. A number of 4541 protein homologs were annotated as unknown, unnamed, hypothetical or expressed proteins.

BLASTX searches were also performed against the complete protein sets of *Arabidopsis thaliana* [33] and *Oryza sativa* [34]. The results were similar, since 12336 and

11996 significant hits were found in Arabidopsis and rice, respectively. BLAST results were parsed, determining for each *Citrus* unigene, the best hit name and description, the extent of the aligned region, the percentage of similarity and the e value. Sequences were classified based on the ORF conservation: 5595 unigenes had very high similarities (80%-100%), 5729 clusters showed high similarities (60 - 80%), while only 1883 unigenes had moderate similarities (40%-60%) and 132 sequences displayed low similarities (30%-40%). No similarities below 30% were obtained. These results indicated that a large number of unigenes (40%) had a very strong match (sequence conservation ≥ 80%) with the top-score hit of the BLASTX results.

The extent of the region of similarity between a given unigene and its best hit protein was also determined, including all High-scoring Segment Pairs (HSPs) for one hit. To this end, the following assumptions were taken: when similarity regions expanded along the complete length of the hit protein, unigenes were supposed to include a complete ORF; if HSPs matched the amino-terminal region of the hit sequence, the cDNA clones from which these ESTs were derived probably contained a complete mRNA partially sequenced; finally, when HSPs located at the carboxy-terminal region of the hit protein, cDNA clones were assumed to correspond to truncated mRNAs. The *Citrus* unigenes were classified according to these criteria, resulting in 4065 complete ORFs, 6082 complete cDNA clones, and 4132 partial or truncated cDNA clones. Taken together the complete ORFs and clones, the *Citrus* EST collection contained at least, 10147 complete cDNA clones.

<footnote><sup>a</sup> Total number of ESTs assembled</footnote>

<footnote><sup>b</sup> Unigenes (singletons and contigs) from a single library</footnote>
<figure><img src="image_2.png" /></figure>

<figure><img src="image_3.png" /><figcaption>**Figure 1**
**EST assembly results.** A – Distribution of ESTs in contigs. B – Number of contigs containing ESTs from 1 or more cDNA libraries.</figcaption></figure>

The species of the best hit sequences produced in the BLASTX search against the non redundant database were registered and classified (data not shown). Not surprisingly, 7725 hits (55.2%) corresponded to *A. thaliana* and 2610 (18.6%) to other eudicots species. About 400 unigenes produced significant best hits from species other than plants (mainly bacteria, fungi, and insects), with a high degree of conservation (= 80%). 298 of these clusters had no significant hits from Arabidopsis or rice, and a BLASTN search performed against the non redundant and EST GenBank databases, confirmed they might be contaminant sequences not originated from *Citrus* species.

A considerable number of sequences (10 to 15%) did not produce a significant hit in the BLAST searches performed. The correlation between the number of sequences that produced no hits and the length of the unigene (Fig 2), showed that 607 out of 954 sequences shorter than 500 bp (63%) did not produce a significant hit, while the percentage was only 7% for longer sequences (1067 out of 14710 unigenes). Since the BLASTX searches were performed against protein databases, the high number of no hits found for the clusters shorter than 500 bp, may indicate that these regions did not contain coding sequences, corresponding to mRNA UTRs. The absence of significant hits was certainly not due to low quality of the shorter

<figure><img src="image_4.png" /><figcaption>**Figure 2**
**BLASTX Annotation.** Relationship between the length (bp) and both the number of unigenes producing no hits (black columns) and the average similarity (grey column) with respect the best hits.</figcaption></figure>

sequences, because the level of similarity of the unigenes was always high regardless the size of the sequence (Fig 2), indicating that shorter sequences had the same quality than longer ones.

The sequences that did not produce significant hits were used as queries in a BLASTN search performed against the GenBank non redundant nucleotide database [32]. Only 40 unigenes produced significant hits showing high similarity levels (80 - 100%). A total of 24 sequences presented similarity with the complete sequence of *Poncirus trifoliata* *Citrus* tristeza virus resistance gene locus, the only genomic BAC clone (282699 bp long) from a *Citrus* species that has been sequenced [35].

The same set of unigenes produced 647 significant hits (conserved region longer than 100 bp and similarity ≥ 80%) in a further BLASTN search against the GenBank EST database [32], indicating that similar transcripts were previously isolated. Further analysis showed that most of the hits derived only from *Citrus* species, suggesting that these clusters might be putative *Citrus* exclusive genes.

**Protein translation and annotation**

A search for domains associated with a Hidden Markov Model profile was intended to improve annotation of the EST collection. To obtain better templates for annotation, the translation of the unigenes consensus sequences into polypeptide ones was carried out with the prot4EST prediction pipeline, which produces robust translations from EST sequences [36]. A BLASTP search was performed with the polypeptide sequences against the GenBank protein database [32], and 77% of them produced the same hit than the original DNA sequences, showing the accuracy of the translations. All polypeptides shorter than 30 aa, and
those shorter than 100 aa without a significant hit were discarded, and the final number of useful protein translations was 14782. The parsing of the BLASTP results confirmed that the unigenes initially classified as complete ORFs coded for complete proteins.

Once protein sequences were obtained and tested, searches against pattern or signature databases were performed. The InterPro database was chosen for these searches, as it unites secondary databases that contain overlapping information on protein families, domains and functional sites [37]. The standalone version of the interProScan tool, that combines the protein function recognition methods of the member databases of InterPro into one application, was used for the analysis. The total number of unigenes that produced hits and the number of different motifs identified varied highly depending on the database size and the analysis method. Motif search against the Pfam database produced the largest number of protein hits (7141) and also identified a high number of different motifs (1666). Table 2 shows the 20 most abundant conserved domains found with Pfam. These motifs were responsible of the 34% of the total hits. Further analysis was carried out on the subset of 2482 proteins predicted as complete with the BLASTX analyses that produced Pfam hits (Table 3). Most of the proteins, 76.15%, displayed a single type of functional domain; proteins combining 2 different motifs accounted for the

20.43%, while polypeptides with 3 or 4 different conserved domains were in a low number (3.42%).

## Gene Ontology Annotation

The Gene Ontology (GO) annotation of the *Citrus* unigene set was performed with BLAST2GO (B2G) [38]. B2G assigns GO annotations through a 3-steps procedure: BLAST against protein databases, retrieval of all GO annotations for a specified number of BLAST hits (Mapping), and GO assignment through an evaluated annotation rule (Annotation). Figure 3A shows the intensity of GO annotation. A total of 10842 unigenes were annotated with 39173 annotations, distributed among the main Gene Ontology categories: Biological Process (11868), Molecular Function (14686) and Cellular Component (12604). There were 5515 sequences annotated with all three GO categories, and 8576 that had at least two annotations. Failure in GO term assignment was due to either a negative result in the BLAST search (NoBLAST, 53%), the absence of GO annotation in any of the BLAST hits (NoMapping, 7%) or because the sequence did not fulfil quality parameters of trustable annotation (NoAnnotation, 40%) (Figure 3B). Noteworthy was that most of these last 47% of sequences, i.e., sequences with a BLAST result without GO term assignment, had best BLAST hits to unknown or hypothetical proteins, indicating the uncertainty in the functional characterization of the hit sequences.

<table><caption>Table 2: The 20 most abundant conserved motifs found with Pfam</caption><thead><tr><th>Pfam Motif Description</th><th>Pfam Motif</th><th>N<sup>a</sup></th><th>Molecular Function <sup>b</sup></th></tr></thead><tbody><tr><td>Leucine Rich Repeat</td><td>PF00560.19</td><td>457</td><td>protein-protein interaction</td></tr><tr><td>WD domain, G-beta repeat</td><td>PF00400.18</td><td>341</td><td>protein-protein interaction</td></tr><tr><td>PPR repeat</td><td>PF01535.9</td><td>218</td><td>RNA binding</td></tr><tr><td>RNA recognition motif. (a.k.a. RRM, RBD, or RNP domain)</td><td>PF00076.10</td><td>190</td><td>nucleic acid binding(GO:0003676)<sup>c</sup></td></tr><tr><td>Protein kinase domain</td><td>PF00069.13</td><td>156</td><td>protein kinase activity (GO:0004672)<sup>c</sup> ATP binding (GO:0005524)<sup>c</sup></td></tr><tr><td>Ankyrin repeat</td><td>PF00023.16</td><td>149</td><td>protein-protein interaction</td></tr><tr><td>EF hand</td><td>PF00036.18</td><td>111</td><td>calcium ion binding (GO:0005509)<sup>c</sup></td></tr><tr><td>Mitochondrial carrier protein</td><td>PF00153.13</td><td>97</td><td>binding (GO:0005488)<sup>c</sup></td></tr><tr><td>Myb-like DNA-binding domain</td><td>PF00249.17</td><td>87</td><td>DNA binding (GO:0003677)<sup>c</sup></td></tr><tr><td>F-box domain</td><td>PF00646.19</td><td>79</td><td>protein-protein interaction</td></tr><tr><td>Tetratricopeptide repeat</td><td>PF00515.14</td><td>77</td><td>protein-protein interaction</td></tr><tr><td>Armadillo/beta-catenin-like repeat</td><td>PF00514.10</td><td>77</td><td>protein-protein interaction</td></tr><tr><td>Zinc finger, C3HC4 type (RING finger)</td><td>PF00097.11</td><td>74</td><td>ubiquitin-protein ligase (GO:0004842)<sup>c</sup> zinc ion binding (GO:0008270)<sup>c</sup></td></tr><tr><td>Zinc knuckle</td><td>PF00098.10</td><td>67</td><td>nucleic acid binding(GO:0003676)<sup>c</sup></td></tr><tr><td>Zinc finger C-x8-C-x5-C-x3-H type</td><td>PF00642.13</td><td>67</td><td>nucleic acid binding(GO:0003676)</td></tr><tr><td>AP2 domain</td><td>PF00847.9</td><td>58</td><td>transcription factor (GO:0003700)<sup>c</sup></td></tr><tr><td>XYPPX repeat</td><td>PF02162.6</td><td>55</td><td>unknown</td></tr><tr><td>Kelch motif</td><td>PF01344.13</td><td>49</td><td>unknown</td></tr><tr><td>Ubiquitin family</td><td>PF00240.12</td><td>45</td><td>protein modification (GO:0006464)<sup>c</sup></td></tr><tr><td>IQ calmodulin-binding motif</td><td>PF00612.14</td><td>43</td><td>calmodulin binding</td></tr></tbody></table>

<footnote><sup>a</sup> Number of motif repeats found in the Citrus protein set.</footnote>

<footnote><sup>b</sup> Molecular function associated with the functional domain.</footnote>

<footnote><sup>c</sup> Molecular function according to the Gene Ontology annotation system.</footnote>
<table><caption>Table 3: Distribution of PFAM motifs in the hypothetical complete proteins</caption><thead><tr><th>Different motifs per protein<sup>a</sup></th><th>Motif number per protein<sup>b</sup></th><th>Protein number<sup>b</sup></th><th>%</th><th>Total Protein number (%)<sup>d</sup></th></tr></thead><tbody><tr><td rowspan="8">1</td><td>1</td><td>1766</td><td>71.1%</td><td rowspan="8">1890 (76.15%)</td></tr><tr><td>2</td><td>69</td><td>2.78%</td></tr><tr><td>3</td><td>13</td><td>0.52%</td></tr><tr><td>4</td><td>21</td><td>0.85%</td></tr><tr><td>5</td><td>7</td><td>0.28%</td></tr><tr><td>6</td><td>8</td><td>0.32%</td></tr><tr><td>7</td><td>4</td><td>0.16%</td></tr><tr><td>9</td><td>2</td><td>0.08%</td></tr><tr><td rowspan="6">2</td><td>2</td><td>452</td><td>18.2%</td><td rowspan="6">507 (20.43%)</td></tr><tr><td>3</td><td>25</td><td>1.01%</td></tr><tr><td>4</td><td>18</td><td>0.73%</td></tr><tr><td>5</td><td>2</td><td>0.08%</td></tr><tr><td>6</td><td>5</td><td>0.20%</td></tr><tr><td>8</td><td>5</td><td>0.20%</td></tr><tr><td rowspan="5">3</td><td>20</td><td>1</td><td>0.04%</td><td rowspan="5">73 (2.94%)</td></tr><tr><td>3</td><td>64</td><td>2.58%</td></tr><tr><td>4</td><td>5</td><td>0.20%</td></tr><tr><td>5</td><td>1</td><td>0.04%</td></tr><tr><td>6</td><td>2</td><td>0.08%</td></tr><tr><td>4</td><td>4</td><td>12</td><td>0.48%</td><td>12 (0.48%)</td></tr></tbody></table>

To provide a general representation of the distribution of *Citrus* gene ontology annotation, the Slim GO Classification for Plants developed at TAIR [33] was obtained, and sets of genes according to broad GO ontology categories were produced [39]. All functional categories in the Biological Processes classification were well represented in the *Citrus* unigene set (Figure 4). Similar results were obtained in the other main GO classes, Molecular function and Cellular Localization (data not shown).

## Characterizing the Citrus Gene Space

In an attempt to characterize the *Citrus* gene space, a first analysis was performed to study the biological context of the novel unigenes. Further studies were addressed to identify candidate genes for molecular markers, gene duplications and conserved gene families.

### Novelty and biological context

The results presented in this work identified more than 5159 sequences that had not been included before in any *Citrus* EST collection, and could be, therefore, novel *Citrus* unigenes. Most of these sequences were derived from the normalized full-length library (4673), that contained a mixed of reproductive and vegetative tissues particularly enriched with fruit tissues, abscission zones and salinity

samples. Thus, the unigene set of this library probably included many low abundant transcripts related to several physiological and developmental processes, including fruit quality, productivity and salinity, the targets of this study. Interestingly, 20% of these novel unigenes corresponded to unknown proteins. The number of novel unigenes included into standard libraries was 486, while 148 of them were annotated as unknown genes. To estimate the input of the standard libraries in terms of gene novelty, sets of unigenes known (or assumed) to be involved in the three processes of interest were selected and the libraries containing them identified.

Fruit quality includes many physical attributes and chemical characteristics of the fruit, such as sugar and acid content, flavour and aroma compounds (organoleptic properties). In addition, *Citrus* fruits contain an extensive array of secondary compounds with pivotal nutritional properties. These traits that are acquired along fruit growth are controlled by primary, intermediate, and secondary metabolic pathways. In order to identify genes with relevant roles in fruit quality, homologs of structural enzymes from some of these pivotal metabolic pathways were searched in *Citrus*. The pathways of known topology involved in the biosynthesis of flavonoids and their pre-

<footnote><sup>a</sup> Number of different conserved motifs found in a protein</footnote>

<footnote><sup>b</sup> Total number of motifs found in a protein</footnote>

<footnote><sup>c</sup> Number of proteins displaying the number of conserved motifs shown in the 2<sup>nd</sup> column</footnote>

<footnote><sup>d</sup> Number of proteins displaying the number of different conserved motifs shown in the 1<sup>st</sup> column</footnote>
<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /><figcaption>**Figure 3**
Gene Ontology Annotation. A – Intensity of GO annotation. The number of unigenes with GO annotations for each of the main Gene Ontology categories, biological process (BP), molecular function (MF) and cellular component (CC) is shown. The total column includes unigenes with GO terms from the three categories. B – GO annotation in contigs (black columns) and singletons (grey columns). Annotated = sequences with functional GO annotation; NoBLAST = sequences with no BLAST results; NoMapping = sequences that produced BLAST hits without ontology annotations; NoAnnotation = sequences that produced BLAST hits with not significant ontology annotations; Total = total number of analyzed sequences.</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>**Figure 4**

**Slim GO Annotation of the *Citrus* unigene set.** Grey columns indicate the number of *Citrus* unigenes included in the different slim GO annotation categories.</figcaption></figure>

cursors, the glycolysis, the Krebs cycle, the oxydative/non-oxydative pentose phosphate pathway, the aromatic aminoacid pathway and the general phenyl propanoid and flavonoid pathways were targeted. To identify the *Citrus* unigenes that putatively coded for structural enzymes of the selected pathways, the EC number [40] of the *Citrus* predicted proteins was determined. The results indicated that these pathways were fully represented (Table 4), and that most putative enzymatic activities were present with a significant degree of redundancy, both in terms of EST (expression) and Unigenes (putative paralogs) numbers. From a total of 220 unigenes (1026 ESTs) related to the enzymatic activities of the pathways analyzed in Table 4, it was found that about 15% of them corresponded to novel unigenes. For example, 7 out of the 12 unigenes assigned to the fructose 1,6 biphosphate aldolase activity in the glycolitic pathway are described in this work for the first time. Furthermore, 36% of the 58 enzymatic activities implicated in the processes displayed novel unigenes.

The lignin biosynthesis pathway was analyzed in detail, since lignin is an important component of the dietary fiber present in the fruit [41], that has important benefits

for human health [42,43]. This pathway has been recently described in *Arabidopsis thaliana* [44,45], and is also found at the AraCyc database [33]. ESTs with significant similarity to any of the Arabidopsis enzymes involved in this route were selected and reassembled with the Staden Gap4 program [46]. The number of unigenes was estimated by considering only those contigs showing non identical consensus sequences that overlapped significantly. For each enzymatic activity described in the lignin pathway a phylogenetic analysis that included the *A. thaliana* and *Citrus* proteins was carried out. These analyses redefined the number of *Citrus* unigenes that could be considered orthologs of the Arabidopsis genes. The number of *Citrus* ESTs and Unigenes related to the Arabidopsis proteins of the lignin pathway are shown in Table 5. All the enzymatic activities involved in this pathway were represented, and for 8 particular Arabidopsis proteins, *Citrus* possessed 2, 3 or even 5 orthologs according to the phylogenetic analysis. A search for enzymes of this pathway performed on the annotation database of the draft sequence of the genome of *Populus trichocarpa* [47], produced similar results, with 22 gene models assigned to the 4-coumarate coenzyme A ligase activity, and 5 gene
<table><caption>Table 4: Primary, Intermediate, and Secondary Metabolic Pathways in Citrus</caption><thead><tr><th>Pathway</th><th>Enzymatic Activity</th><th>EC reaction</th><th>Unigenes<sup>a</sup></th><th>ESTs<sup>b</sup></th></tr></thead><tbody><tr><td rowspan="8">Glycolysis</td><td>Hexokinase</td><td>2.7.1.1</td><td>4 (2)<sup>c</sup></td><td>6 (4)<sup>d</sup></td></tr><tr><td>Glucose-6-phosphoisomerase</td><td>5.3.1.9</td><td>6 (1)</td><td>16 (1)</td></tr><tr><td>6-Phosphofructokinase</td><td>2.7.1.11</td><td>10</td><td>35</td></tr><tr><td>Fructose 1,6-biphosphate aldolase</td><td>4.1.2.13</td><td>12 (7)</td><td>71 (9)</td></tr><tr><td>Triose phosphate isomerase</td><td>5.3.1.1</td><td>8 (1)</td><td>22 (1)</td></tr><tr><td>Glyceraldehyde 3-P dehydrogenase</td><td>1.2.1.12</td><td>8 (2)</td><td>52 (8)</td></tr><tr><td>Phosphoglycerate kinase</td><td>2.7.2.3</td><td>6 (1)</td><td>15 (3)</td></tr><tr><td>Phosphoglycerate mutase</td><td>5.4.2.1</td><td>1</td><td>5</td></tr><tr><td rowspan="2">Enolase</td><td>4.2.1.11</td><td>7 (4)</td><td>28 (5)</td></tr><tr><td>Pyruvate kinase</td><td>2.7.1.40</td><td>10 (1)</td><td>69 (9)</td></tr><tr><td rowspan="8">Tricarboxylic acid cycle (Krebs cycle)</td><td>Citrate synthase</td><td>4.1.3.7</td><td>2</td><td>13</td></tr><tr><td>Aconitase</td><td>4.2.1.3</td><td>6 (1)</td><td>31 (1)</td></tr><tr><td>Isocitrate dehydrogenase</td><td>1.1.1.42</td><td>3</td><td>59</td></tr><tr><td>a-Ketoglutarate dehydrogenase complex</td><td>1.2.4.2</td><td>No hits</td><td>No hits</td></tr><tr><td>Succinyl-CoA synthetase</td><td>6.2.1.5</td><td>2</td><td>14</td></tr><tr><td>Succinate dehydrogenase</td><td>1.3.5.1</td><td>2</td><td>20</td></tr><tr><td>Fumarase</td><td>4.2.1.2</td><td>3 (2)</td><td>30 (5)</td></tr><tr><td>Malate dehydrogenase</td><td>1.1.99.16</td><td>4</td><td>24</td></tr><tr><td rowspan="7">Oxidative/nonoxidative pentose phosphate pathway</td><td>Glucose 6-P-1-dehydrogenase</td><td>1.1.1.49</td><td>4</td><td>4</td></tr><tr><td>6-Phosphogluconolactonase</td><td>3.1.1.31</td><td>2</td><td>10</td></tr><tr><td>6-Phosphogluconate dehydrogenase</td><td>1.1.1.44</td><td>15</td><td>47</td></tr><tr><td>Ribose-5-P isomerase</td><td>5.3.1.6</td><td>3</td><td>12</td></tr><tr><td>Ribose-5-P 3-epimerase</td><td>5.1.3.1</td><td>No hits</td><td>No hits</td></tr><tr><td>Transketolase</td><td>2.2.1.1</td><td>3 (2)</td><td>22 (11)</td></tr><tr><td>Transaldolase</td><td>2.2.1.2</td><td>4 (2)</td><td>20 (17)</td></tr><tr><td rowspan="16">Aromatic amino acid biosynthesis</td><td>3-Deoxy-D-arabino-heptulosonate 7-P synthase</td><td>4.1.2.15</td><td>No hits</td><td>No hits</td></tr><tr><td>3-Dehydroquinate synthase</td><td>4.2.3.4</td><td>2</td><td>2</td></tr><tr><td>3-Dehydroquinate dehydratase</td><td>4.2.1.10</td><td>2</td><td>2</td></tr><tr><td>Shikimate 5 dehydrogenase</td><td>1.1.1.25</td><td>1</td><td>14</td></tr><tr><td>Shikimate kinase</td><td>2.7.1.71</td><td>7</td><td>10</td></tr><tr><td>5-Enolpyruvoylshikimate 3-P synthase</td><td>2.5.1.19</td><td>No hits</td><td>No hits</td></tr><tr><td>Chorismate synthase</td><td>4.2.3.5</td><td>2</td><td>2</td></tr><tr><td>Chorismate mutase</td><td>5.4.99.5</td><td>1</td><td>1</td></tr><tr><td>Prephenate dehydratase</td><td>4.2.1.51</td><td>3</td><td>5</td></tr><tr><td>Prephenate dehydrogenase</td><td>1.3.1.12</td><td>2 (1)</td><td>18 (2)</td></tr><tr><td>Aromatic amino acid transaminase</td><td>2.6.1.57</td><td>6</td><td>18</td></tr><tr><td>Anthranilate synthase</td><td>4.1.3.27</td><td>6 (1)</td><td>46 (29)</td></tr><tr><td>Anthranilate phosphoribosyl transferase</td><td>2.4.2.18</td><td>4</td><td>6</td></tr><tr><td>Phosphoribosylanthranilate synthase</td><td>5.3.1.24</td><td>No hits</td><td>No hits</td></tr><tr><td>Indol-3-glycerol phosphate synthase</td><td>4.1.1.48</td><td>4 (2)</td><td>6 (2)</td></tr><tr><td>Trp synthase</td><td>4.2.1.20</td><td>8 (1)</td><td>30 (1)</td></tr><tr><td>Phe ammonia-lyase</td><td>4.3.1.5</td><td>5</td><td>13</td></tr><tr><td>Cinnamate 4-hydroxilase</td><td>1.14.13.11</td><td>1</td><td>3</td></tr><tr><td>4-Coumarate coenzyme A ligase</td><td>6.2.1.12</td><td>10 (1)</td><td>40 (1)</td></tr><tr><td rowspan="6">Flavonol biosynthesis</td><td>Naringenin chalcone synthase</td><td>2.3.1.74</td><td>6</td><td>48</td></tr><tr><td>Chalcone isomerase</td><td>5.5.1.6</td><td>3</td><td>17</td></tr><tr><td>Flavanone 3-hydroxylase</td><td>1.14.11.9</td><td>2 (1)</td><td>12 (1)</td></tr><tr><td>Flavonol 3-hydroxylase</td><td>1.14.13.21</td><td>6</td><td>7</td></tr><tr><td>Flavonol synthase</td><td>No E.C.</td><td>2</td><td>19</td></tr><tr><td>Flavonol 3-O-glucosyltransferase</td><td>No E.C.</td><td>No hits</td><td>No hits</td></tr><tr><td rowspan="2">Flavonoid biosynthetic pathway: anthocyanin biosynthesis</td><td>Dihydroflavonol 4-reductase</td><td>1.1.1.219</td><td>No hits</td><td>No hits</td></tr><tr><td>Leucoanthocyanidin dioxygenase</td><td>No E.C.</td><td>2</td><td>33</td></tr></tbody></table>
<table><caption>Table 4: Primary, Intermediate, and Secondary Metabolic Pathways in Citrus (Continued)</caption><tbody><tr><td></td><td>UDP-flavonol 3-O-glucosyltransferase</td><td>2.4.1.-</td><td>No hits</td><td>No hits</td></tr><tr><td></td><td>Anthocyanin 5-O-glucosyltransferase</td><td>NoE.C.</td><td>2</td><td>9</td></tr><tr><td></td><td>Anthocyanin 5-aromatic acyltransferase</td><td>2.3.4.153</td><td>6</td><td>36</td></tr><tr><td></td><td>Anthocyanin permease</td><td>NoE.C.</td><td>No hits</td><td>No hits</td></tr><tr><td rowspan="2"><b>Proanthocyanidin biosynthesis</b></td><td>Anthocyanidin reductase</td><td>1.3.1.77</td><td>No hits</td><td>No hits</td></tr><tr><td>Leucoanthocyanidin reductase</td><td>No E.C.</td><td>2 (1)</td><td>4 (2)</td></tr></tbody></table>

models to the caffeoyl-CoA O-methyltransferase activity. These data suggest that multiple orthologs of all Arabidopsis lignin proteins could be apparently found if enough number of *Citrus* ESTs were analyzed, indicating the relevance of the lignification process in *Citrus*. The comparison of the unigenes involved in lignin biosynthesis (Table 5) against *Citrus* ESTs from GenBank, showed that 7 out of the 22 orthologs produced not hits, indicating that they were novel genes reported in this work for the first time.

Additional new unigenes implicated in fruit quality were selected based on published information [18,19,48] from relevant pathways of lipids and fatty acid metabolism and degradation [GenBank:DY258371, GenBank:DY258372, GenBank:DY258373, GenBank:DY258374, Contig0424, Contig4859, Contig5406, GenBank:DY258378, GenBank:DY258379, GenBank:DY258380, GenBank:DY258381, Contig0330], synthesis and accumulation of citric acid [GenBank:DY258383, GenBank:DY258384, Contig5931], sugar [GenBank:DY258396] and nitrate transport [Contig3203, Contig4271, GenBank:DY258401, GenBank:DY258402], and chlorophyll synthesis [GenBank:DY258395]. The analyses indicated that 64% of these novel genes were found in the normalized library, while 14% of them were isolated from FruitTF, one of the two fruit specific libraries. The other unigenes belonged to stress and abscission libraries.

The analysis of genes related to productivity was initially focused in 3 families that had been previously associated with the abscission process, the auxin responsive factors (9 novel unigenes), the receptor protein kinases (35 novel unigenes) and the EREBP (ethylene responsive element binding protein, 4 novel unigenes). The results showed that only one singleton [GenBank:EH405902] of the first family, Contig5401 of the second and Contig5227 of the third one, derived from standard abscission libraries (AbsAOv1, AbsDev, and AbsCOv1), while the remaining members (45), were isolated from the normalized library. Since many of the processes implicated in abscission are controlled by the selective removal of short-lived regul-

tory proteins, we also analyzed the occurrence of the ubiquitin/26S proteasome pathway [49] among the novel sequences. This component, deeply involved in protein degradation, has not yet been related to abscission and, therefore, no previous information is available in this regard. Interestingly, the only member of E2s Ub-conjugating enzymes [GenBank:DY258370] found in the set of novel unigenes was detected in the AbsDev library. In addition, 4 members out of 23 E3s Ub-ligases [Contigs 5267 and 5546, GenBank:DY257277 and GenBank:DY258093], were exclusively obtained from the abscission-related libraries. These putative unigenes may participate in the removal of repressor elements during the organ separation [23,50,51]. The other Ub-ligases were mostly present in the normalized library. Several cell wall structural proteins [GenBank:DY256701, GenBank:DY257041, GenBank:DY258445 and GenBank:DY258901] and two specific glycosyl hydrolases [GenBank:DY257803 and GenBank:DY258004] were exclusive of the abscission-related libraries, suggesting that abscission may also implicate active remodelling of cell walls [52]. Aside from the normalized library, the abscission libraries that were strongly enriched with specific abscission zone tissues, showed a relatively high number of novel exclusive unigenes (162).

For the analyses of novel citrus genes potentially involved in abiotic stress, the following well established biological functions were investigated: sodium Na+/H+ antiporters [GenBank:DY258370, GenBank:DY305688 and GenBank:DY300954], that are probably involved in sodium detoxification [53]; the Calcineurin B gene homolog (Contig 5589) [54], stress-induced and/or -activated protein kinases [GenBank:DY278709, Contig0907, Contig1053 and GenBank:DY291464] and the mechanical ion channel-domain containing protein [GenBank:DY262334], that are likely implicated in NaCl-associated signal transduction mechanisms; aldehyde dehydrogenases involved in detoxification [GenBank:DY301300] [55]; two genes of the inositol metabolism [GenBank:DY304982, GenBank:DY260177, GenBank:DY261021, GenBank:DY270505]; genes associated with lipid metabolism such as the phosphoinositide-

<footnote>ªNumber of citrus unigenes annotated with the emzymatic activity</footnote>

<footnote>ᵇTotal number of ESTs corresponding with the unigenes</footnote>

<footnote>The number between brackets indicates the number of unigenes (c) or ESTs (d) that are first reported in this work</footnote>
<table><caption>Table 5: Lignin Biosynthesis Pathway</caption><thead><tr><th>Enzymatic Activity</th><th>Reaction EC</th><th>Gene Name</th><th>Citrus Unigenes<sup>a</sup></th><th>Phyl. Anal.<sup>b</sup></th><th>ESTs</th></tr></thead><tbody><tr><td rowspan="7">4-coumarate<br>coenzyme A ligase</td><td rowspan="7">6.2.1.12</td><td>At1g51680</td><td>1</td><td>1</td><td>4</td></tr><tr><td>At1g62940</td><td>1</td><td>1</td><td>5</td></tr><tr><td>At1g65060</td><td>1 (1)</td><td>1</td><td>7</td></tr><tr><td>At3g21240</td><td>1 (1)</td><td>1</td><td>9</td></tr><tr><td>At4g05160</td><td>1</td><td>1</td><td>1</td></tr><tr><td>At5g63380</td><td>3 (1)</td><td>3</td><td>13</td></tr><tr><td>At1g20510</td><td>1</td><td>1</td><td>5</td></tr><tr><td colspan="2"><strong>Total</strong></td><td><strong>7</strong></td><td></td><td><strong>9</strong></td><td></td></tr><tr><td>5-hydroxy<br>coniferaldehyde o-<br>methyltransferase</td><td>2.1.1.68</td><td>At5g54160</td><td>5 (1)</td><td>5</td><td>41</td></tr><tr><td>caffeoyl-CoA O-<br>methyltransferase</td><td>2.1.1.104</td><td>At4g34050</td><td>3</td><td>2</td><td>14</td></tr><tr><td rowspan="4">cinnamoyl-CoA<br>reductase</td><td rowspan="4">1.2.1.44</td><td>At1g80820</td><td>1</td><td>1</td><td>2</td></tr><tr><td>At2g23910</td><td>1</td><td>1</td><td>10</td></tr><tr><td>At4g30470</td><td>1</td><td>1</td><td>7</td></tr><tr><td>At5g14700</td><td>1</td><td>1</td><td>10</td></tr><tr><td colspan="2"><strong>Total</strong></td><td><strong>4</strong></td><td></td><td><strong>4</strong></td><td></td></tr><tr><td rowspan="6">cinnamyl-alcohol<br>dehydrogenase</td><td rowspan="6">1.1.1.195</td><td>At1g09490</td><td>1</td><td>1</td><td>1</td></tr><tr><td>At1g51410</td><td>1</td><td>1</td><td>2</td></tr><tr><td>At4g27250</td><td>1</td><td>1</td><td>1</td></tr><tr><td>At5g19440</td><td>2 (1)</td><td>2</td><td>1</td></tr><tr><td>At3g19450</td><td>1</td><td>1</td><td>6</td></tr><tr><td>At4g34230</td><td>1</td><td>1</td><td>1</td></tr><tr><td colspan="2"><strong>Total</strong></td><td><strong>6</strong></td><td></td><td><strong>7</strong></td><td></td></tr><tr><td>coniferyl aldehyde 5-<br>hydroxylase</td><td>NIL</td><td>At4g36220</td><td>3 (2)</td><td>3</td><td>15</td></tr><tr><td>coumarate 3-<br>hydroxylase</td><td>1.14.13.36</td><td>At2g40890</td><td>2</td><td>2</td><td>25</td></tr><tr><td>feruloyl coenzyme A<br>reductase</td><td>1.2.1.44</td><td>At1g15950</td><td>2</td><td>2</td><td>32</td></tr><tr><td rowspan="3">UDP-glucose 4-<br>epimerase</td><td rowspan="3">1.2.1.44</td><td>At2g02400</td><td>2</td><td>2</td><td>17</td></tr><tr><td>At2g33590</td><td>1</td><td>1</td><td>13</td></tr><tr><td>At5g58490</td><td>1</td><td>1</td><td>62</td></tr><tr><td colspan="2"><strong>Total</strong></td><td><strong>3</strong></td><td></td><td><strong>4</strong></td><td></td></tr></tbody></table>

<footnote>a Number *Citrus* unigenes annotated assigned to the enzymatic activity</footnote>

<footnote>b Number of Unigenes that clustered with the Arabidopsis proteins associated with the enzymatic activity, in the phlyogenetic analysis.</footnote>
specific phospholipase C [GenBank: DY260755] and the lipoxygenases [GenBank: DY258546 and Contig5406]; and the membrane-associated salt-inducible protein [GenBank: DY301464]. Moreover, the following biological functions related to acclimatization to osmotic shock were searched: two different NCED4 genes involved in ABA biosynthesis (Contig0189 and Contig0309) [24]; and two genes involved in trehalose metabolism [GenBank: DY280731 and GenBank: DY294040]. Lastly, cell tolerance mechanisms universally linked to different abiotic stresses, represented by heat shock proteins and molecular chaperons [GenBank: DY261699, GenBank: DY258174, GenBank: DY303256, GenBank: DY270445, GenBank: DY258682, GenBank: DY258174, GenBank: DY271306, GenBank: DY303256, GenBank: DY270445, GenBank: DY257994, GenBank: DY260682][56]; and uncharacterized stress-responsive genes [Contig0907, Contig2771, GenBank: DY270558 and GenBank: DY260600] were also analyzed. The results of this search indicated that 42% of these abiotic stress related genes were detected in the normalized library, while 34% of them was found in fruit libraries (FruitTF and PhII-IIIVesicles1). Only 1 unigene, a putative sodium Na+/H+ antiporter [GenBank: DY305688], was found as a singleton in a salinity-related library, LSH, whereas only Contig5589 (Calcineurin B gene) contained ESTs exclusively derived from KCl-Salt1, another salinity library.

Overall, these preliminary estimates showed that most of the novel genes, presumably implicated in fruit quality, abscission, and salinity responses, were effectively recovered in the normalized library. The contribution of the fruit and abscission libraries to the set of unigenes related to fruit quality or abscission, could be roughly estimated between 5 and 15%. The lower contribution of the salinity standard libraries (less than 5%) may be due to the abundance of unspecific cross-responses among multiple abiotic and biotic stresses. For an accurate estimation of these figures, however, confirmation of gene specificity appears to be mandatory.

### Molecular markers

The use of genetic and molecular markers is crucial to facilitate the identification and cloning of genes of agro-nomic interest [57], while single copy genes are usually good candidates to be used as markers. To identify conserved A. thaliana orthologs present as single copy genes in both Arabidopsis and *Citrus*, a database with 3700 Arabidopsis single copy genes was obtained from the Compositae Genome Project Database [58], and used in a BLAST search with the *Citrus* unigene set as queries. A total of 726 *Citrus* sequences showed an unambiguous single strong BLAST hit, and reciprocal BLAST searches (Arabidopsis single copy genes versus the EST assemblies) pro-

duced the same results. The outcome of this BLAST search was compared with that obtained in the BLAST search performed against all *Citrus* ESTs from GenBank. The results showed that 129 unigenes did not generate any hit, while 445 clusters only produced hits with similarities higher than 95%, suggesting that these were ESTs probably derived from the same transcript. Although this analysis is not conclusive, the absence of hits or the occurrence of extremely high similarities, suggested that these 574 unigenes are strong candidates to be conserved orthologs of Arabidopsis single copy genes.

### Gene duplications

The BLAST search with the Arabidopsis single copy genes also produced 234 sequences with 2 or more strong *Citrus* hits. These cases were further investigated as they might be indicative of gene duplication events produced in the *Citrus* genome. In many cases, Unigenes showing the same Arabidopsis hit did not overlap, indicating that they may derive from the same transcript but were not assembled in the same contig, and therefore cannot be considered to be different genes. Finally, 18 Arabidopsis single copy genes showed strong similarity with two overlapping *Citrus* unigenes. These clusters presented the same Arabidopsis protein as their best hit, supporting the hypothesis that they are paralog genes in *Citrus* [see Additional file 1].

### Gene Family analysis

Comparative genomics was used to characterize the conserved gene families in A. thaliana and *Citrus* species. There are currently 930 gene families, comprising 6399 genes, described at the *Arabidopsis thaliana* Information Resource database [59]. The presence of these gene families in the *Citrus* unigene set was explored, allocating the *Citrus* clusters in the gene families based on the best Arabidopsis significant hit obtained. About 3000 *Citrus* unigenes were assigned to 724 families, and 52 super families, showing that 78% and 92% of the Arabidopsis families and superfamilies were represented in the *Citrus* EST collection. To exemplify the potential for *Citrus* improvement of the information included in the EST collection, two gene families with relevant agronomic interest were selected and analyzed in detail: the ammonium transporter family intimately related to plant nutritional efficiency and the glycoside hydrolase family 20, implicated in sugar synthesis in fruit.

The Arabidopsis high-affinity ammonium transporter family is composed of six members: five proteins that form the AtAMT1 subfamily [60] and a member, AtAMT2, that is distantly related to the AtAMT1 subfamily [61]. Similarity searches showed that 60 *Citrus* ESTs were significantly similar to the Arabidopsis ammonium transporters. These ESTs that corresponded to 6 unigenes (4 contigs and 2 singletons) were translated into protein with prot4EST. A
multiple alignment of the *A. thaliana* and *Citrus* sequences was performed with Clustal X, and a phylogenetic tree was constructed with the neighbor joining method (Figure 5A). The tree showed that 3 *Citrus* unigenes closely related to AtAMT2 grouped together in a cluster supported by very high bootstrap values. Thus, the AMT2 subfamily of ammonium transporters in *Citrus* comprised, at least, 3 proteins, suggesting the occurrence of a number of duplication events.

Glycosyltransferase family 20 is composed of proteins with known α, α-trehalose-phosphate synthase UDP-forming activity, and in *A. thaliana* comprises proteins At1G05590 and At3G55260. A total of 15 *Citrus* ESTs displaying significant similarity with the glycosyltransferase family 20 proteins were assembled into 4 unigenes grouped in 2 contigs and 2 singletons. As above, phylogenetic analysis showed that the *Citrus* and Arabidopsis proteins clustered together, with high bootstrap values supporting the clade (Figure 5B). This analysis suggests that the glycosyltransferase family 20 included 4 members in *Citrus* species while in *A. thaliana* it contained only two proteins.

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

<figure><img src="image_9.png" /><figcaption>**Figure 5**
**Phylogenetic analysis.** A – Phylogenetic tree of the ammonium transporter family from *A. thaliana* and *Citrus* species. B – Phylogenetic tree of the glycosyltransferase family 20 from *A. thaliana* and *Citrus* species. Glycosyltransferase family 19 cluster was collapsed in a black triangle and used as an outside group.</figcaption></figure>

## Discussion

Citrus is the main fruit tree crop in the world. However, traditional breeding for *Citrus* cultivar improvement faces many serious impediments due to the unusual combination of biological characteristics of *Citrus*, their low genetic diversity and the long-term nature of tree breeding. Genomic technology can overcome these limitations providing new tools, for example, to produce more efficient varieties and rootstocks and to identify new genes, alleles or genotypes of agronomic relevance. Improvement of knowledge of the transcriptome is one of the first tasks that have to be developed in order to understand the developmental biology of the plants and how these respond to environmental stresses. This work that pursues this goal provides a deep insight into the *Citrus* transcriptome specifically related to three major commercial traits i.e improved fruit quality, higher yield and tolerance to environmental stresses, especially salinity.

Towards this objective, 10 cDNA libraries representing particular treatments and tissues from selected varieties and rootstocks differing in fruit quality, resistance to abscission and tolerance to salinity were generated to provide a large and enriched expressed sequence tag collection. The assembly of these sequences, more than 52600 ESTs, allowed the identification of 15660 transcription units. The results of this analysis are comparable to previous reports in *Solanum tuberosum*, that detected 19892 unigenes from 61949 ESTs [62], or in *Sorghum bicolor*, with 16801 unique transcripts derived from 55783 ESTs [63]. The data showed that all sequences from the *Citrus* species analyzed, from both this study and databases, were almost identical, suggesting that the differential behaviour of these cultivars during normal fruit growth or when facing environmental adverse conditions is more likely associated with differences in gene regulation rather than with sequence divergence. This result is not unexpected since *Citrus* posses a high level of phenotypic diversity while global genetic diversity, analyzed with molecular markers, appears to be very low or practically null. A large effort was made to determine the real number of different transcripts represented in the EST collection. It is well known that the accuracy of EST clustering is affected by various error sources, such as sequencing mistakes, contaminant sequences and the presence of products of chimeric splicing. The most common error occurs when different ESTs from the same gene are falsely separated into two or more clusters [64]. To overcome this difficulty the level of redundancy was estimated comparing all unigenes with each other, and clustering them in *supercontigs*, that are more likely to correspond to real transcripts. The level of redundancy was estimated to be 26%, a value similar to that obtained in sugarcane, for example [65]. This first restriction suggests that the likely number
of unigenes in the *Citrus* collection is closer to 13900 rather than to 15660.

It was also crucial to determine the occurrence of contaminant sequences, mostly from microorganisms, since many samples were taken from open field. The presence of contaminant sequences (mainly from bacteria and fungi) is a general problem not attended in any of the EST projects we have examined. For instances, a BLASTN search performed with several contaminant sequences found in this work against the viridiplantae section of the GenBank EST database revealed a considerable number of ESTs regarded as plant sequences that really corresponded to fungi species (data not shown). Thus, the analysis reported in this work may help to prevent the presence of sequences from contaminant species in the databases. Determining the species of the unigenes best hit sequences helped to identify putative contaminants, allowing not only a more precise estimation of the real number of *Citrus* transcripts but also criteria for microarray EST selection. Since about 400 Unigenes were believed to be contaminant sequences from other species, the number of *Citrus* expressed genes was reduced to 13500.

A relevant observation of this work is that more than 38% of the 13500 unigenes (5159 sequences) are novel *Citrus* unigenes. EST sequences were obtained from two kinds of cDNA libraries, normalized full length and standard libraries. The normalized library was generated with a wide variety of reproductive and vegetative tissues, enriched with developing fruits, abscission zones and salinity samples. In the first strategy, the normalization process very effectively increased low abundant transcripts, since the bulk of the novel unigenes described (4673) derived from this library. The standard libraries that were constructed from either samples of fruits, abscission zones or salt-treated organs, were generated with the idea of providing transcripts specifically expressed at these particular tissues and organs without increasing redundancy. To estimate the contribution of the standard libraries in terms of gene novelty, a set of unigenes presumably involved in the three processes of interest was selected and the libraries containing them identified. Although for an accurate evaluation of these contributions confirmation of tissue specificity appears to be mandatory, these preliminary estimates showed that most of the novel genes were certainly recovered in the normalized library. Moreover, the contribution of the specific standard libraries to the unigene set maybe roughly estimated to be between 4 and 15%.

The primary homology searches performed against different databases allowed annotation of most unigenes, with more than 73% of them displaying a similarity degree

higher than 60%. These results agreed those obtained in previous *Citrus* [12] or sugarcane [65] EST projects. It was also shown that most of the sequences that did not produce significant hits in the BLASTX searches were shorter than 500 bp (Fig 2) and probably did not carry coding sequences. Additional efforts were performed to characterize these sequences, with supplementary BLASTN searches against the non redundant and EST nucleotide databases. These analyses gave rise to the suggestion that 647 ESTs of the *Citrus* unigene set may correspond to *Citrus* exclusive genes since the significant hits they produced were only for *Citrus* sequences, in spite of the more than 8.5 million EST sequences derived from plant species deposited at the GenBank,

Further improvement of the annotation was carried out through searches performed against secondary databases, composed of patterns or signatures. Although these prediction methods can work with DNA sequences, the error prone nature of ESTs, mainly shifts in the reading frame (missing or inserted bases) or ambiguous bases, may result in inaccuracies and loss of information. Thus, a crucial step in annotation is the robust translation of the ESTs to yield predicted polypeptides. Polypeptide sequences possesses a better template for almost all annotation tools, including InterPro and Pfam, and allow the assembly of more accurate multiple sequence alignments. High quality polypeptide predictions can be applied to functional annotation and post-genomic study in a similar way to those available for completed genomes. In the work presented here, the protein translation was performed with Prot4EST, a prediction pipeline that incorporates freely available software (ESTscan, Decoder, HSP tiling) to produce final translations that are more accurate than those derived from any single method [36]. The use of the inter-ProScan tool [37], allowed simultaneous search for motifs against 9 databases. This search produced significant results for almost 11000 predicted proteins, including 342 unigenes that did not have significant hits in the BLAST searches. From the 20 most abundant motifs found in *Citrus* with Pfam, 7 of them also were included in the top 20 list at the Pfam database, enforcing the accuracy of the analysis and the representativity of the *Citrus* EST collection. The molecular functions associated with these 20 motifs can be grouped in 4 categories: protein-protein interaction (47.26%), nucleic acid binding (27.51%), protein modification and binding (14.9%), and calcium metabolism (6.17%), which indicates the relative significance of these cellular functions in *Citrus*.

The distribution of motifs on the polypeptide sequences predicted to be complete proteins showed that the bulk of sequences displayed a single motif (76%). Proteins carrying 2 or more motifs showed unlike signatures rather than repeats of the same motif. For instance, the number
of proteins with 2 different signatures (452) was six times higher than the number of proteins that had 2 identical motifs (69). Similar relationships were found for other number of repeats.

The Gene Ontology annotation of the *Citrus* unigenes was performed with BLAST2GO (B2G), a recently developed BLAST-based GO annotation software [38]. The B2G approach uses multiple BLAST hits to search for functional annotations and assigns GO terms to the query sequence applying an annotation algorithm that considers HSP length, e-value, percentage of similarity, Evidence Code of the source annotations and the topology of the Gene Ontology. This is in contrast to most EST projects that perform annotation solely by direct assignment of the GO terms to the best hit of BLAST searches [12,66,67]. The B2G method has shown to have a high annotation recall and has been used in other EST projects in eukaryotes [68].

Metabolic pathways responsible of important agronomic traits were further surveyed to determine the extent of representation of these pathways within the *Citrus* unigene set. In addition to the finding that most enzymatic steps were represented by *Citrus* homologs, a preliminary estimate of gene duplication based on the occurrence of paralogous sequences was also provided. Defining such relationships and understanding functional diversification of paralogs is an important field of research in genomics-assisted crop improvement. Lignin biosynthetic pathway was the object of a deeper analysis since *Citrus* fruits are very rich in products with beneficial effects in preventing cancer, diabetes, ...etc such as fiber [42,69]. Dietary fiber, that consists of non digested structural and storage polysaccharides and lignin, lowers cholesterol levels and helps to normalize blood glucose and insulin levels [70]. The detailed analysis of lignin biosynthesis pathway, carried out in *Citrus*, indicated that in comparison with Arabidopsis, *Citrus* possessed at least 9 additional enzymatic activities involved in lignin synthesis (Table 5). Furthermore, the results obtained for gene models from *Populus trichocarpa* [47], also appears to support the idea that the extensive formation of secondary xylem in tree species, requiring high levels of lignin synthesis may have been the origin of the expansion of the genes involved in this pathway.

More than 570 unigenes have been suggested to be possible conserved orthologs of Arabidopsis single copy genes. Recent studies have indicated that ancient polyploidy is common across angiosperm lineages and in fact, the genomes of all angiosperms may have been influenced by at least one genome-wide duplication event [71]. Despite such events, single-copy, apparently orthologous gene sets have been identified in a broad range of angiosperms

[72,73]. Selection against duplicates may be maintaining these genes as single copy, and therefore are precious markers for comparative genetic and physical mapping, and also for phylogenetic analyses [74]. Identification of such genes in *Citrus* species is mandatory to perform this kind of analyses [75]. The study also revealed a number of genes that might be duplicated in the *Citrus* genome, while remained as single copy genes in Arabidopsis, although the possibility of finding additional copies of these genes could not be discarded, when the whole transcriptome of *Citrus* is available. If these duplications are the result of individual events or were caused by a genome-wide duplication cannot be answered with the current information.

Comparative genomics was also used to obtain an overview of conserved gene families in *Citrus*. All Arabidopsis gene families studied were well represented in the *Citrus* EST collection, although the number of their members was generally smaller, probably because the unigene set was only a partial representation. For the same reason, the finding of families that in *Citrus* clearly outnumbered their Arabidopsis counterparts is highly significant. The phylogenomic analysis performed on the gene families of ammonium transporters and glycosyltransferases supported this idea confirming the occurrence of additional members in the *Citrus* families. Ammonium is one of the prevalent nitrogen sources for growth and development of higher plants including *Citrus*. The ammonium transporter family is composed in Arabidopsis of 5 AMT1 related genes and AMT2, which is more closely related to ammonium transporters from prokaryotes than to AMT1. AtAMT2 is likely to play a significant role in moving ammonium between the apoplast and symplast of cells throughout the plant [61]. Interestingly, there are three AMT2 like genes in *Citrus* (Fig 6A). Glycosyltransferases are a ubiquitous group of enzymes that catalyse the transfer of a sugar moiety from an activated sugar donor onto saccharide or non-saccharide acceptors. Although many glycosyltransferases catalyse chemically similar reactions, they display remarkable diversity in their donor, acceptor and product specificity and thereby generate a potentially infinite number of glycoconjugates, oligo- and polysaccharides. [76]. Thus, the additional members found in this family, might be related to the complexity of sugar synthesis that takes place in the *Citrus* fruits.

## Conclusion

The assembly of more than 54000 *Citrus* ESTs from five cultivars differing in basic fruit developmental aspects, such as major traits for fruit quality and production, and in the responses to environmental conditions, provides an unprecedented insight of the *Citrus* transcriptome. This study contributes new tools for *Citrus* genetic and genomic analyses. The unigene set, composed of ~13000
putative different transcripts, including more than 5000 novel Citrus genes, was assigned with putative functions based on similarity, GO annotations and protein domains. In addition, comparative genomics was used to analyze the *Citrus* transcriptome, and evidences for numerous cases of gene duplication events were presented. The similarity analyses indicated that the sequences of the genes belonging to the varieties and rootstocks studied were essentially identical suggesting that the differential behaviour of these species cannot be attributed to main sequence divergences. This set of processed EST sequences has greatly contributed to the development of a new *Citrus* microarray.

# Methods

## 1. Plant material

The *Citrus* genotypes used to generate the cDNA libraries were the varieties *Citrus clementina*, (cv Clementina de Nules), and *C. sinensis* (cvs Navelina and Washington Navel), and the rootstocks *C. reshni* (cv *Cleopatra mandarin*) and *C. sinensis × Poncirus trifoliata* (cv Carrizo citrange). Their characteristics are as follows. Clementine is a mandarin of elevated fruit quality, high ovary and fruitlet abscission and moderate salt tolerance. Washington Navel is a late sweet orange that generally shows pre-harvest abscission. In contrast, Navelina, an early orange variety, exhibits low fruit abscission but higher salt sensitivity. *Cleopatra mandarin* is an efficient salt tolerant rootstock while the hybrid Carrizo citrange shows high salt sensitive).

## 2. Normalized Full Length Library (NFL)

### Tissue Samples and Treatments Description

All samples included in the normalized full-length library were harvested from *Citrus clementina* (cv Clementina de Nules). They were composed of the following tissues and organs: developing vegetative tips and buds, dormant buds, developing leaves, shoots, internodes and roots, abscission zones from leaves, flowers, ovaries and fruits, flowers and inflorescences, growing and senescent ovaries, developing fruitlets (stages I & II), flavedo from growing, ripening and, senescent fruits and fruit flesh (juice sacs, stages I, II & III). The library also included leaves subjected to different treatments: short- and long-term salinity, drought and rehydration, mineral deficiencies, alkaline and calcareous soils, low and high temperature, flooding, oxidative stress, wounding, insect attacks, and elicitors (harpin) treatments. All tissues were frozen in liquid nitrogen and equal amounts of homogenized tissues were mixed in a single sample for total RNA extraction.

### Library Construction

Full-length cDNA synthesis was carried out with Invitrogen proprietary RNase H reduction reverse transcriptase "cocktail" for mRNA isolation, 5' cap full-length enrichment, and the reduction of oligo(dT)-priming. Normalization was carried out by self-subtraction, with Invitrogen technology, as described by manufacturer. Normalization produced a 24 fold average reduction of the abundant clones, (from 0.16% abundance to 0.0065%). PCMVSSport6.1 was used as a cloning vector.

## 3. Standard Libraries

### Tissue Samples and Treatments Description

- **Fruit-TF:** parthenocarpic fruits of *Citrus clementina* (cv Clementina de Nules) mandarin were harvested from adult trees grown grafted onto Carrizo citrange rootstock (*Citrus sinensis × Poncirus trifoliata*) in a homogeneous orchard under normal culture practices. Flavedo (exocarp) samples were isolated from fruits collected on July 28 (69 days post anthesis, dpa), July 24 (85 dpa), August 2 (94 dpa), October 11 (164 dpa), November 18 (202 dpa), November 25 (209 dpa), December 13 (227 dpa) and January 9 (254 dpa). Samples of fruit flesh, consisting of juice vesicles (endocarp) including the segments with their membranes and vascular bundles, were taken from fruits collected on May 13 (13 dpa) and June 10 (41 dpa). Samples were frozen under liquid nitrogen and stored at -80 °C until RNA isolation. Mixtures of equal amounts of poly-A+ RNA from the samples were used.

- **PhII-IIIIVesicles1:** fruit juice vesicles from Clementine grafted onto Carrizo citrange were taken at one month intervals: July 8 (69 dpa), August 2 (94 dpa), September 12 (135 dpa), October 16 (169 dpa), November 14 (198 dpa) and December 17 (231 dpa). A mixture of equal amounts of poly-A+ RNA from the six samples was used.

- **AbsDev:** laminar abscission zone and surrounding tissues (petiole and blade) of developing leaves were harvested from Clementine on Carrizo citrange.

- **AbsCFruit1:** abscission zone C and surrounding tissues of ripe fruits were harvested from *Citrus sinensis* (cv Washington navel) scions on Carrizo.

- **AbsCOv:** abscission zone C and surrounding tissues of ethylene-treated ovary explants were harvested from Clementinescions on Carrizo.

- **AbsAOv1:** abscission zone A and surrounding tissues of ethylene-treated ovary explants at "petal fall" stage were harvested from Clementine scions on Carrizo.

- **LSH:** leaves were harvested from one-year-old *Citrus sinensis* (cv Navelina) scions grafted onto Cleopatra (*Citrus reshni*) rootstock cultured under salinity conditions. Potted plants were grown in greenhouse conditions and subjected to regular irrigations (three times per week) with 25 mM NaCl:CaCl2 solutions for 60 days.
**KCl-Salt1:** non-suberized roots, enriched in distal (actively growing) root portions were harvested from 1 year-old Cleopatra mandarin seedlings. Potted plants grown in greenhouse conditions were subjected to Cl-starvation and resupply treatments at different times with 50 and 100 mM KCl.

**EHR:** young roots from Carrizo citrange were collected 3, 6, 12, and 24 hours after water stress treatment and 1, 6, and 10 hours after re-watering. A mixture of equal amounts of poly-A+ RNA from the different samples was used.

### RNA Extraction

For AbsDev, AbsCFRuit1, AbsCOv1, AbsAOv1, and EHR libraries, total RNA was isolated from frozen tissue using the standard guanidine protocol [77]. For FruitTF library, total RNA was isolated from frozen tissue using the RNeasy Plant Mini Kit (Qiagen) and treated with RNase-free DNase (Qiagen) through column purification according to the manufacturer's instructions. For KCL-Salt1 library, total RNA was isolated from frozen tissue using acid phenol extraction and Lithium Chloride precipitation method [78]. In all cases, RNA quality was assessed by spectrophotometry and gel electrophoresis [77].

### Poly(A)+ RNA Isolation

Poly(A)+RNA was isolated from a mixture of equal amounts of total RNA from all samples using the Oligotex mRNA Mini Kit (Qiagen) following manufacturer's instructions.

### Construction protocols and cloning vectors

KCl-Salt1, AbsDev, AbsCFRuit1, and FruitTF cDNA libraries were constructed using the CloneMiner cDNA Library Construction Kit (Invitrogen) with the pDONRTM222 vector. AbsCOv1, AbsAOv1, and EHR cDNA libraries were constructed with SMART cDNA Library Construction Kit (Clontech) and pTriplEx2 as the cloning vector. SLH library was constructed with Stratagene cDNA synthesis kit and the pBluescript SK (-) vector. PhII-III-Vesicles1, and EHR libraries were constructed using the UNI-ZAP XR and Gigapack III Gold kits from Stratagene and ë-ZAP II cloning vector.

## 4. EST assembly and annotation

DNA templates were prepared using the 96-well alkaline lysis DNA method. Sequencing was performed using the ABI Big Dye Terminator Cycle Sequence Ready Reaction as described by manufacturer, with the T7 forward primer in 96 well plates in an automatic ABI 3730.

The software phred was used for base calling, and Cross-match for vector masking [79]. Reading assembly was performed with the CAP3 program [30], using read quality

and defaults parameters. Similarity searches were performed with the standalone version of BLAST [31], against the NCBI non redundant protein, nucleotide and EST databases [32], the *Arabidopsis thaliana* protein set from TAIR database [33] and the *Oryza sativa* protein set from TIGR Rice database [80]. Parsing of the BLAST results was performed with the Bio::SearchIO module [81] from the Bioperl package [82].

Protein translation was performed with prot4EST polypeptide prediction pipeline [36], which combines different methods like ESTscan [83], DECODER [84] and similarity search results (BLASTX) to produce accurate translations. Motif search was performed with the standalone version of the interProScan tool that combines the protein function recognition methods of the database members of InterPro into one single application [37]. The InterPro database unites the following secondary databases: Uniprot [85], Panther [86], PROSITE [87], PRINTS [88], Pfam [89], ProDom [90], SMART [91], TIGRFAMs [92], PIR [93] and SUPERFAMILY [94].

Gene Ontology annotation of unigenes was performed with BLAST2GO [38]. Blast2GO is a user adjustable tool that utilizes BLAST to find homologous sequences for a set of query sequences and returns an evaluated annotation from the gene ontology annotations present in the BLAST hits of each sequence. B2G parameters were: NCBI non-redundant DB for BLAST search, 20 hits maximum for BLAST result, 100 nt as minimum HSP-length to retain putative annotating hits and default Evidence Code Weights for Gene Ontology annotation that assigns high ECWs to experimental-based and curated annotations while penalized electronic and non-curated annotations. Minimum values for BLAST e-value and % similarity of the BLAST result were e-06 and 55% respectively and ultimate annotation cut-off value was set to 55. This set of parameters was shown to provide the most reliable results in the annotation of Arabidopsis sequences [38]. GOSlim annotations of the *Citrus* unigenes were also generated with the B2G software using the plant GOSlim mapping provided in TAIR.

## 5. Gene space analysis

Single copy gene set from *A. thaliana* was obtained from the Compositae Genome Project Database [58] and used as query in a BLAST search against a database generated with the *Citrus* unigenes. These results were compared with those obtained with the BLASTX search performed with the *Citrus* unigenes against the Arabidopsis complete protein set. Only Unigenes with a unique Arabidopsis significant hit that matched the results obtained with the first BLAST search were considered to be putative orthologs of the Arabidopsis single copy genes. A similar approach was used to detect possible gene duplications, selecting those
*Citrus* unigenes that had as significant hits the same *A. thaliana* protein. ESTs corresponding to the selected unigenes were reassembled with the Staden package [46] to confirm that the Unigenes were overlapping rather than identical sequences.

The *A. thaliana* gene family dataset was obtained from TAIR [33] and the Arabidopsis best hit of the *Citrus* unigenes was used to find *Citrus* representatives of relevant families. ESTs from the overrepresented *Citrus* families were selected and reassembled with the Staden package [46], and only overlapping non-identical consensus sequences were considered for further analysis. For the phylogenetic study of the ammonium transporter and glycosyltransferase 20 families, a multiple proteins sequence alignment was carried out with ClustalX [95], genetic distances were calculated with the protein correction for the poisson method [96], and phylogenetic trees were constructed with the neighbor joining tree method [97], using the Molecular Evolutionary Genetic Analysis software package MEGA3 [98].

**List of abbreviations**

aa, amino acid

bp, base pair

EST, expressed sequence tag

GO, gene ontology

nt, nucleotide

TAIR, The Arabidopsis Information Resource

TIGR, The Institute for Genomic Research

**Authors' contributions**

JT carried out bioinformatics (assembly, annotation and phylogenetic analyses) and drafted the manuscript. AC participated in the GO annotation of the sequences and helped to draft the manuscript. MC participated in sample collection and construction of cDNA libraries. JMC participated in sample collection, construction of cDNA libraries, and helped to draft the manuscript. FT participated in sample collection and construction of cDNA libraries. JA participated in sample collection and construction of cDNA libraries. EA participated in sample collection and construction of cDNA libraries. FA participated in sample collection and construction of cDNA libraries. XA contributed to the annotation of the sequences. JB participated in sample collection and construction of cDNA libraries. GS participated in sample collection and construction of cDNA libraries. BC performed EST sequencing. CD performed EST sequencing. SG participated in the GO anno-

tation of the sequences. DJI generated samples for library construction. FL generated samples for library construction. RM contributed to the coordination and design of the project. PO contributed to the coordination and design of the project. PW performed EST sequencing. MT conceived the study, coordinated the project and drafted the manuscript. All authors read and approved the final manuscript.

**Additional material**

**Additional file 1**

*Putative gene duplications in *Citrus*. The table contains the *A. thaliana* loci that have been found duplicated in *Citrus*.*

Click here for file

[http://www.biomedcentral.com/content/supplementary/1471-2164-8-31-S1.doc]

**Acknowledgements**

Most of the sequencing work was developed at Genoscope through a "Sequencage a grande echelle" 2003 grant. Work at Genoscope was funded by CNRG. Additional funding from Spanish Ministerio de Educación y Ciencia through grants GEN2001-4885-c05-03 and AGL2003-08502-C04-01, from Instituto Nacional de Investigaciones Agrarias through grants RTA04-013 and RTA05-247, from Conselleria de Agricultura, Pesca y Alimentación de la Generalitat Valenciana through IVIA grant 5309, and from Commission of the European Communities through contract 015453 is gratefully acknowledged. Samples from Clementine roots were kindly provided by Dr. L Navarro at IVIA.

**References**

1. Ollitrault P, Jacquemond C, Dubois C, Luro F: **Citrus.** In *Genetic diversity of cultivated tropical plants* Edited by: Hamon P, Seguin M, Perrier X, Glaszmann X. Montpellier, CIRAD; 2003:193-197.

2. The Arabidopsis Genome Initiative AGI: **Analysis of the genome sequence of the flowering plant Arabidopsis thaliana.** *Nature* 2000, **408**(6814):796-815.

3. Yu J, Hu S, Wang J, Wong GK, Li S, Liu B, Deng Y, Dai L, Zhou Y, Zhang X, Cao M, Liu J, Sun J, Tang J, Chen Y, Huang X, Lin W, Ye C, Tong W, Cong L, Geng J, Han Y, Li L, Li W, Hu G, Li J, Liu Z, Qi Q, Li T, Wang X, Lu H, Wu T, Zhu M, Ni P, Han H, Dong W, Ren X, Feng X, Cui P, Li X, Wang H, Xu X, Zhai W, Xu Z, Zhang J, He S, Xu J, Zhang K, Zheng X, Dong J, Zeng W, Tao L, Ye J, Tan J, Chen X, He J, Liu D, Tian W, Tian C, Xia H, Bao Q, Li G, Gao H, Cao T, Zhao W, Li P, Chen W, Zhang Y, Hu J, Liu S, Yang J, Zhang G, Xiong Y, Li Z, Mao L, Zhou C, Zhu Z, Chen R, Hao B, Zheng W, Chen S, Guo W, Tao M, Zhu L, Yuan L, Yang H: **A draft sequence of the rice genome (Oryza sativa L. ssp. indica).** *Science* 2002, **296**(5565):79-92.

4. Goff SA, Ricke D, Lan TH, Presting G, Wang R, Dunn M, Glazebrook J, Sessions A, Oeller P, Varma H, Hadley D, Hutchison D, Martin C, Katagiri F, Lange BM, Moughamer T, Xia Y, Budworth P, Zhong J, Miguel T, Paszkowski U, Zhang S, Colbert M, Sun WL, Chen L, Cooper B, Park S, Wood TC, Mao L, Quail P, Wing R, Dean R, Yu Y, Zharkikh A, Shen R, Sahasrabudhe S, Thomas A, Cannings R, Gutin A, Pruss D, Reid J, Tavtigian S, Mitchell J, Eldredge G, Scholl T, Miller RM, Bhatnagar S, Adey N, Rubano T, Tusneem N, Robinson R, Feldhaus J, Macalma T, Oliphant A, Briggs S: **A draft sequence of the rice genome (Oryza sativa L. ssp. japonica).** *Science* 2002, **296**(5565):92-100.

5. Tuskan GA, DiFazio S, Jansson S, Bohlmann J, Grigoriev I, Hellsten U, Putnam N, Ralph S, Rombauts S, Salamov A, Schein J, Sterck L, Aerts A, Bhalerao RR, Bhalerao RP, Blaudez D, Boerjan W, Brun A, Brunner A, Busov V, Campbell M, Carlson J, Chalot M, Chapman J, Chen GL,
Cooper D, Coutinho PM, Couturier J, Covert S, Cronk Q, Cunningham R, Davis J, Degroeve S, Dejardin A, dePhamphilis C, Detter J, Dirks B, Dubchak I, Duplessis S, Ehlting J, Ellis B, Gendler K, Goodstein D, Gribskov M, Grimwood J, Groover A, Gunter L, Hamberger B, Heinze B, Helariutta Y, Henrissat B, Holligan D, Holt R, Huang W, Islam-Faridi N, Jones S, Jones-Rhoades M, Jorgensen R, Joshi C, Kangasjarvi J, Karlsson J, Kelleher C, Kirkpatrick R, Kirst M, Kohler A, Kalluri U, Larimer F, Leebens-Mack J, Leple JC, Locascio P, Lou Y, Lucas S, Martin F, Montanini B, Napoli C, Nelson DR, Nelson C, Nieminen K, Nilsson O, Pereda V, Peter G, Philippe R, Pilate G, Poliakov A, Razumovskaya J, Richardson P, Rinaldi C, Ritland K, Rouze P, Ryaboy D, Schmutz J, Schrader J, Segerman B, Shin H, Siddiqui A, Sterky F, Terry A, Tsai CJ, Uberbacher E, Unnebe P, Vahala J, Wall K, Wessler S, Yang G, Yin T, Douglas C, Marra M, Sandberg G, Van de Peer Y, Rokhsar D: *The Genome of Black Cottonwood, Populus trichocarpa* (Torr. & Gray). *Science* 2006, 313(5793):1596-1604.

6. Ewing RM, Kahla AB, Poirot O, Lopez F, Audic S, Claverie JM: Large-scale statistical analyses of rice ESTs reveal correlated patterns of gene expression. *Genome Res* 1999, 9(10):950-959.

7. Dirlewanger E, Graziano E, Joobeur T, Garriga-Caldere F, Cosson P, Howad W, Arus P: *Comparative mapping and marker-assisted selection in Rosaceae fruit crops*. *Proc Natl Acad Sci U S A* 2004, 101(26):9891-9896.

8. Feingold S, Lloyd J, Norero N, Bonierbale M, Lorenzen J: *Mapping and characterization of new EST-derived microsatellites for potato (Solanum tuberosum L.)*. *Theor Appl Genet* 2005, 111(3):456-466.

9. Lu C, Hawkesford MJ, Barraclough PB, Poulton PR, Wilson ID, Barker GL, Edwards KJ: *Markedly different gene expression in wheat grown with organic or inorganic fertilizer*. *Proc Biol Sci* 2005, 272(1575):1901-1908.

10. Firnhaber C, Puhler A, Kuster H: EST sequencing and time course microarray hybridizations identify more than 700 Medicago truncatula genes with developmental expression regulation in flowers and pods. *Planta* 2005, 222(2):269-283.

11. Baxter CJ, Sabar M, Quick WP, Sweetlove LJ: *Comparison of changes in fruit gene expression in tomato introgression lines provides evidence of genome-wide transcriptional changes and reveals links to mapped QTLs and described traits*. *J Exp Bot* 2005, 56(416):1591-1604.

12. Forment J, Gadea J, Huerta L, Abizanda L, Agusti J, Alamar S, Alos E, Andres F, Arribas R, Beltran JP, Berbel A, Blazquez MA, Brumos J, Canas LA, Cercos M, Colmenero-Flores JM, Conesa A, Estables B, Gandia M, Garcia-Martinez JL, Gimeno J, Gisbert A, Gomez G, Gonzalez-Candelas L, Granell A, Guerri J, Lafuente MT, Madueno F, Marcos JF, Marques MC, Martinez F, Martinez-Godoy MC, Miralles S, Moreno P, Navarro L, Pallas V, Perez-Amador MA, Perez-Valle J, Pons C, Rodrigo I, Rodriguez PL, Royo C, Serrano R, Soler G, Tadeo F, Talon M, Terol J, Trenor M, Vaello L, Vicente O, Vidal C, Zacarias L, Conejero V: *Development of a citrus genome-wide EST collection and cDNA microarray as resources for genomic studies*. *Plant Molecular Biology* 2005, 57(3):375-391.

13. Fujii H, Shimada T, Eendo T, Shimizu T, Omura M: 29,228 Citrus ESTs- Collection and Analysis Toward The Functional Genomics Phase. In *Plant & Animal Genomes XIV Conference Town & Country Convention Center*. San Diego, CA, USA; 2006.

14. Machado MA, Souza AA, Targon ML, Takita MA, Freitas-Astua J, Filho HC, Amaral AM, Palmieri DA, Boscariol-Camargo R, Cristofani M, Carlos EF, Reis MS: Current Situation Of Citrus Genome Project In Brazil (CitEST). In *Plant & Animal Genomes XIV Conference Town & Country Convention Center*. San Diego, CA, USA; 2006.

15. Roose ML, Federici CT, Lyon MP, Fenton RD, Wanamaker S, Close TJ: *Citrus EST sequencing and prospects for a high-density microarray*. In *Plant & Animal Genomes XII Conference Town & Country Convention Center*. San Diego, CA, USA; 2004.

16. Bausher M, Shatters R, Chaparro J, Dang P, W. H, Niedz R: An expressed sequence tag (EST) set from Citrus sinensis L. Osbeck whole seedlings and the implications of further perennial source investigations. *Plant Science* 2003, 165(2):415-422.

17. Close TJ, Wanamaker S, Lyon M, Mei G, Davies C, Roose ML: A GeneChip® For Citrus. In *Plant & Animal Genomes XIV Conference Town & Country Convention Center*. San Diego, CA, USA; 2006.

18. Cercos M, Soler G, Iglesias DJ, Gadea J, Forment J, Talon M: Global Analysis of Gene Expression During Development and Ripening of Citrus Fruit Flesh. A Proposed Mechanism for Citric Acid Utilization. *Plant Mol Biol* 2006.

19. Alos E, Cercos M, Rodrigo MJ, Zacarias L, Talon M: Regulation of color break in citrus fruits. Changes in pigment profiling and gene expression induced by gibberellins and nitrate, two ripening retardants. *J Agric Food Chem* 2006, 54(13):4888-4895.

20. Iglesias DJ, Tadeo FR, Legaz F, Primo-Millo E, Talon M: In vivo sucrose stimulation of colour change in citrus fruit epicarps: Interactions between nutritional and hormonal signals. *Physiol Plant* 2001, 112(2):244-250.

21. Ben-Cheikh W, Perez-Botella J, Tadeo FR, Talon M, Primo-Millo E: Pollination Increases Gibberellin Levels in Developing Ovaries of Seeded Varieties of Citrus. *Plant Physiol* 1997, 114(2):557-564.

22. Talon M, Hedden P, Primo-Millo E: Gibberellins in Citrus sinensis: A comparison between seeded and seedless varieties. *Journal of Plant Growth Regulation* 1990, 9(1):201-206.

23. Gomez-Cadenas A, Tadeo FR, Talon M, Primo-Millo E: Leaf Abscission Induced by Ethylene in Water-Stressed Intact Seedlings of Cleopatra Mandarin Requires Previous Abscisic Acid Accumulation in Roots. *Plant Physiol* 1996, 112(1):401-408.

24. Agusti J, Zapater M, Iglesias DJ, Cercos M, Tadeo FR, Talon M: Differential expression of putative 9-cis-epoxycarotenoid dioxigenases and abscisic acid accumulation in water stressed vegetative and reproductive tissues of citrus. *Plant Science* 2006, In press:

25. Moya JL, Primo-Millo E, Talon M: Morphological factors determining salt tolerance in citrus seedlings: the shoot to root ratio modulates passive root uptake of chloride ions and their accumulation in leaves. *Plant, Cell and Environment* 1999, 22(11):1425-1433.

26. Romero-Aranda R, Moya JL, Tadeo FR, Legaz F, Primo-Millo E, Talon M: Physiological and anatomical disturbances induced by chloride salts in sensitive and tolerant citrus: beneficial and detrimental effects of cations. *Plant, Cell and Environment* 1998, 21(12):1243-1253.

27. Moya JL, Gomez-Cadenas A, Primo-Millo E, Talon M: Chloride absorption in salt-sensitive Carrizo citrange and salt-tolerant Cleopatra mandarin citrus rootstocks is linked to water use. *J Exp Bot* 2003, 54(383):825-833.

28. Iglesias DJ, Levy Y, Gomez-Cadenas A, Tadeo FR, Primo-Millo E, Talon M: Nitrate improves growth in salt-stressed citrus seedlings through effects on photosynthetic activity and chloride accumulation. *Tree Physiol* 2004, 24(9):1027-1034.

29. Talon M, Zacarias L, Primo-Millo E: Hormonal changes associated with fruit set and development in mandarins differing in their parthenocarpic ability. *Physiologia Plantarum* 1990, 79(2):400-406.

30. Huang X, Madan A: CAP3: A DNA sequence assembly program. *Genome Res* 1999, 9(9):868-877.

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

32. National Center for Biotechnology Information [http://www.ncbi.nlm.nih.gov/](http://www.ncbi.nlm.nih.gov/)

33. The Arabidopsis Information Resource [http://www.arabidopsis.org/](http://www.arabidopsis.org/)

34. TIGR Rice Genome Annotation [http://www.tigr.org/tdb/e2k1/osa1/](http://www.tigr.org/tdb/e2k1/osa1/)

35. Yang ZN, Ye XR, Molina J, Roose ML, Mirkov TE: Sequence analysis of a 282-kilobase region surrounding the citrus Tristeza virus resistance gene (Ctv) locus in Poncirus trifoliata L. Raf. *Plant Physiol* 2003, 131(2):482-492.

36. Wasmuth J, Blaxter M: prot4EST: Translating Expressed Sequence Tags from neglected genomes. BMC Bioinformatics 2004, 5(1):187.

37. Quevillon E, Silventoinen V, Pillai S, Harte N, Mulder N, Apweiler R, Lopez R: InterProScan: protein domains identifier. Nucl Acids Res 2005, 33(suppl_2):W116-120.

38. Conesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, Robles M: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. *Bioinformatics* 2005, 21(18):3674-3676.

39. Berardini TZ, Mundodi S, Reiser L, Huala E, Garcia-Hernandez M, Zhang P, Mueller LA, Yoon J, Doyle A, Lander G, Moseyko N, Yoo D, Xu I, Zoeckler B, Montoya M, Miller N, Weems D, Rhee SY: *Func-
tional Annotation of the Arabidopsis Genome Using Controlled Vocabularies. *Plant Physiol* 2004, 135(2):745-755.

40. Bairoch A: The ENZYME database in 2000. *Nucleic Acids Res* 2000, 28(1):304-305.

41. Kay RM: Dietary fiber. *J Lipid Res* 1982, 23(2):221-242.

42. Reddy BS: Prevention of colon carcinogenesis by components of dietary fiber. *Anticancer Research* 1999, 19(5A):3681-3683.

43. Aggarwal BB, Shishodia S: Molecular targets of dietary agents for prevention and therapy of cancer. *Biochem Pharmacol* 2006, 71(10):1397-1421.

44. Boerjan W, Ralph J, Baucher M: Lignin biosynthesis. *Annu Rev Plant Biol* 2003, 54:519-546.

45. Humphreys JM, Chapple C: Rewriting the lignin roadmap. *Curr Opin Plant Biol* 2002, 5(3):224-229.

46. Staden R: The Staden sequence analysis package. *Mol Biotechnol* 1996, 5(3):233-241.

47. The International Populus Genome Consortium [http://www.ornl.gov/sci/ipgc/].

48. Iglesias DJ, Lliso I, Tadeo FR, Talon M: Regulation of photosynthesis through source: sink imbalance in citrus is mediated by carbohydrate content in leaves. *Physiologia Plantarum* 2002, 116(4):563-572.

49. Smalle J, Vierstra RD: The Ubiquitin 26S proteasome proteolytic pathway. *Annual Review of Plant Biology* 2004, 55(1):555-590.

50. Iglesias DJ, Tadeo FR, Primo-Millo E, Talon M: Fruit set dependence on carbohydrate availability in citrus trees. *Tree Physiol* 2003, 23(3):199-204.

51. Gomez-Cadenas A, Tadeo FR, Primo-Millo E, Talon M: Involvement of abscisic acid and ethylene in the responses of citrus seedlings to salt shock. *Plant Physiology* 1998, 103:475-484.

52. Sexton R, Roberts JA: Cell Biology of Abscission. *Annual Review of Plant Physiology* 1982, 33(1):133-162.

53. Pardo JM, Cubero B, Leidi EO, Quintero FJ: Alkali cation exchangers: roles in cellular homeostasis and stress tolerance. *J Exp Bot* 2006, 57(5):1181-1199.

54. Mendoza I, Quintero FJ, Bressan RA, Hasegawa PM, Pardo JM: Activated calcineurin confers high tolerance to ion stress and alters the budding pattern and cell morphology of yeast cells. *J Biol Chem* 1996, 271(38):23061-23067.

55. Sunkar R, Bartels D, Kirch HH: Overexpression of a stress-inducible aldehyde dehydrogenase gene from Arabidopsis thaliana in transgenic plants improves stress tolerance. *Plant J* 2003, 35(4):452-464.

56. Vierling E: The Roles of Heat Shock Proteins in Plants. *Annual Review of Plant Physiology and Plant Molecular Biology* 1991, 42(1):579-620.

57. Varshney RK, Graner A, Sorrells ME: Genomics-assisted breeding for crop improvement. *Trends Plant Sci* 2005, 10(12):621-630.

58. Compositae Genome Project Database [http://cgpdb.ucdavis.edu].

59. Rhee SY, Beavis W, Berardini TZ, Chen G, Dixon D, Doyle A, Garcia-Hernandez M, Huala E, Lander G, Montoya M, Miller N, Mueller LA, Mundodi S, Reiser L, Tacklind J, Weems DC, Wu Y, Xu I, Yoo D, Yoon J, Zhang P: The Arabidopsis Information Resource (TAIR): a model organism database providing a centralized, curated gateway to Arabidopsis biology, research materials and community. *Nucleic Acids Res* 2003, 31(1):224-228.

60. Gazzarrini S, Lejay L, Gojon A, Ninnemann O, Frommer WB, von Wiren N: Three Functional Transporters for Constitutive, Diurnally Regulated, and Starvation-Induced Uptake of Ammonium into Arabidopsis Roots. *Plant Cell* 1999, 11(5):937-948.

61. Sohlenkamp C, Wood CC, Roeb GW, Udvardi MK: Characterization of Arabidopsis AtAMT2, a High-Affinity Ammonium Transporter of the Plasma Membrane. *Plant Physiol* 2002, 130(4):1788-1796.

62. Ronning CM, Stegalkina SS, Ascenzi RA, Bougri O, Hart AL, Utterbach TR, Vanaken SE, Riedmuller SB, White JA, Cho J, Pertea GM, Lee Y, Karamycheva S, Sultana R, Tsai J, Quackenbush J, Griffiths HM, Restrepo S, Smart CD, Fry WE, Van Der Hoeven R, Tanksley S, Zhang P, Jin H, Yamamoto ML, Baker BJ, Buell CR: Comparative analyses of potato expressed sequence tag libraries. *Plant Physiol* 2003, 131(2):419-429.

63. Pratt LH, Liang C, Shah M, Sun F, Wang H, Reid SP, Gingle AR, Paterson AH, Wing R, Dean R, Klein R, Nguyen HT, Ma H, Zhao X, Morishige DT, Mullet JE, Cordonnier-Pratt MM: Sorghum

Expressed Sequence Tags Identify Signature Genes for Drought, Pathogenesis, and Skotomorphogenesis from a Milestone Set of 16,801 Unique Transcripts. *Plant Physiol* 2005, 139(2):869-884.

64. Wang JPZ, Lindsay BG, Leebens-Mack J, Cui L, Wall K, Miller WC, dePamphilis CW: EST clustering error evaluation and correction. *Bioinformatics* 2004, 20(17):2973-2984.

65. Vettore AL, da Silva FR, Kemper EL, Souza GM, da Silva AM, Ferro MIT, Henrique-Silva F, Giglioti EA, Lemos MVF, Coutinho LL, Nobrega MP, Carrer H, Franca SC, Bacci M Jr., Goldman MHS, Gomes SL, Nunes LR, Camargo LEA, Siqueira WJ, Van Sluys MA, Thiemann OH, Kuramae EE, Santelli RV, Marino CL, Targon MLPN, Ferro JA, Silveira HCS, Marini DC, Lemos EGM, Monteiro-Vitorello CB, Tambor JHM, Carraro DM, Roberto PG, Martins VG, Goldman GH, de Oliveira RC, Truffi D, Colombo CA, Rossi M, de Araujo PG, Sculaccio SA, Angella A, Lima MMA, de Rosa VE Jr, Siviero F, Coscrato VE, Machado MA, Grivet L, Di Mauro SMZ, Nobrega FG, Menck CFM, Braga MDV, Telles GP, Cara FAA, Pedrosa G, Meidanis J, Arruda P: Analysis and Functional Annotation of an Expressed Sequence Tag Collection for Tropical Crop Sugarcane. *Genome Res* 2003, 13(12):2725-2735.

66. Udall JA, Swanson JM, Haller K, Rapp RA, Sparks ME, Hatfield J, Yu Y, Wu Y, Dowd C, Arpat AB, Sickler BA, Wilkins TA, Guo JY, Chen XY, Scheffler J, Taliercio E, Turley R, McFadden H, Payton P, Klueva N, Allen R, Zhang D, Haigler C, Wilkerson C, Suo J, Schulze SR, Pierce ML, Essenberg M, Kim HR, Llewellyn DJ, Dennis ES, Kudrna D, Wing R, Paterson AH, Soderlund C, Wendel JF: A global assembly of cotton ESTs. *Genome Res* 2006, 16(3):441-450.

67. Moser C, Segala C, Fontana P, Salakhudtinov I, Gatto P, Pindo M, Zyprian E, Toepfer R, Grando MS, Velasco R: Comparative analysis of expressed sequence tags from different organs of Vitis vinifera L. *Funct Integr Genomics* 2005, 5(4):208-217.

68. Ma J, Morrow D, Fernandes J, Walbot V: Comparative profiling of the sense and antisense transcriptome of maize lines. *Genome Biology* 2006, 7(3):R22.

69. Sugiura M, Ohshima M, Ogawa K, Yano M: Chronic administration of Satsuma mandarin fruit (Citrus unshiu Marc.) improves oxidative stress in streptozotocin-induced diabetic rat liver. *Biol Pharm Bull* 2006, 29(3):588-591.

70. Marlett JA, McBurney MI, Slavin JL: Position of the American Dietetic Association: health implications of dietary fiber. *J Am Diet Assoc* 2002, 102(7):993-1000.

71. Adams KL, Wendel JF: Polyploidy and genome evolution in plants. *Curr Opin Plant Biol* 2005, 8(2):135-141.

72. Wang CJR, Harper L, Cande WZ: High-Resolution Single-Copy Gene Fluorescence in Situ Hybridization and Its Use in the Construction of a Cytogenetic Map of Maize Chromosome 9. *Plant Cell* 2006, 18(3):529-544.

73. Fransz PF, Stam M, Montijn B, Hoopen RT, Wiegant J, Kooter JM, Oud O, Nanninga N: Detection of single-copy genes and chromosome rearrangements in Petunia hybrida by fluorescence in situ hybridization. *The Plant Journal* 1996, 9(5):767-774.

74. Soltis D, Carlson J, Farmerie W, Wall PK, Ilut D, Solow T, Mueller L, Landherr L, Hu Y, Buzgo M, Kim S, Yoo MJ, Frohlich M, Perl-Treves R, Schlarbaum S, Bliss B, Zhang X, Tanksley S, Oppenheimer D, Soltis P, Ma H, dePamphilis C, Leebens-Mack J: Floral gene resources from basal angiosperms for comparative genomics research. *BMC Plant Biology* 2005, 5(1):5.

75. Tyagi AK, Khurana JP: Plant molecular biology and biotechnology research in the post-recombinant DNA era. *Adv Biochem Eng Biotechnol* 2003, 84:91-121.

76. Coutinho PM, Deleury E, Davies GJ, Henrissat B: An evolving hierarchical family classification for glycosyltransferases. *J Mol Biol* 2003, 328(2):307-317.

77. Sambrook J, Fritsch E, Maniatis T: Molecular Cloning. In *A Laboratory Manual* 2nd edition. Cold Spring Harbor Laboratory Press, Cold Spring Harbor, N.Y.; 1989.

78. Ecker JR, Davis RW: Plant Defense Genes are Regulated by Ethylene. *PNAS* 1987, 84(15):5202-5206.

79. Ewing B, Green P: Base-Calling of Automated Sequencer Traces Using Phred. *II Error Probabilities*. 1998, 8(3):186-194.

80. The Institute for Genomic Research [http://www.tigr.org/].

81. BioPerl [http://www.bioperl.org/wiki/Main_Page].

82. Stajich JE, Block D, Boulez K, Brenner SE, Chervitz SA, Dagdigian C, Fuellen G, Gilbert JG, Korf I, Lapp H, Lehvaslaiho H, Matsalla C, Mun-
gall CJ, Osborne BI, Pocock MR, Schattner P, Senger M, Stein LD, Stupka E, Wilkinson MD, Birney E: **The Bioperl toolkit: Perl modules for the life sciences.** Genome Res 2002, 12(10):1611-1618.

83. Lottaz C, Iseli C, Jongeneel CV, Bucher P: **Modeling sequencing errors by combining Hidden Markov models.** Bioinformatics 2003, 19(90002):ii103-112.

84. Fukunishi Y, Hayashizaki Y: **Amino acid translation program for full-length cDNA sequences with frameshift errors.** Physiol Genomics 2001, 5(2):81-87.

85. Bairoch A, Apweiler R, 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 LSL: **The Universal Protein Resource (UniProt).** Nucl Acids Res 2005, 33(suppl 1):D154-159.

86. Mi H, Lazareva-Ulitsky B, Loo R, Kejariwal A, Vandergriff J, Rabkin S, Guo N, Muruganujan A, Doremieux O, Campbell MJ, Kitano H, Thomas PD: **The PANTHER database of protein families, sub-families, functions and pathways.** Nucl Acids Res 2005, 33(suppl 1):D284-288.

87. Hulo N, Sigrist CJA, Le Saux V, Langendijk-Genevaux PS, Bordoli L, Gattiker A, De Castro E, Bucher P, Bairoch A: **Recent improvements to the PROSITE database.** Nucl Acids Res 2004, 32(90001):D134-137.

88. Attwood TK, Bradley P, Flower DR, Gaulton A, Maudling N, Mitchell AL, Moulton G, Nordle A, Paine K, Taylor P, Uddin A, Zygouri C: **PRINTS and its automatic supplement, prePRINTS.** Nucl Acids Res 2003, 31(1):400-402.

89. Bateman A, Coin L, Durbin R, Finn RD, Hollich V, Griffiths-Jones S, Khanna A, Marshall M, Moxon S, Sonnhammer ELL, Studholme DJ, Yeats C, Eddy SR: **The Pfam protein families database.** Nucl Acids Res 2004, 32(90001):D138-141.

90. Bru C, Courcelle E, Carrere S, Beausse Y, Dalmar S, Kahn D: **The ProDom database of protein domain families: more emphasis on 3D.** Nucl Acids Res 2005, 33(suppl 1):D212-215.

91. Letunic I, Copley RR, Schmidt S, Ciccarelli FD, Doerks T, Schultz J, Ponting CP, Bork P: **SMART 4.0: towards genomic data integration.** Nucl Acids Res 2004, 32(90001):D142-144.

92. Haft DH, Selengut JD, White O: **The TIGRFAMs database of protein families.** Nucl Acids Res 2003, 31(1):371-373.

93. Wu CH, Yeh LSL, Huang H, Arminski L, Castro-Alvear J, Chen Y, Hu Z, Kourtesis P, Ledley RS, Suzek BE, Vinayaka CR, Zhang J, Barker WC: **The Protein Information Resource.** Nucl Acids Res 2003, 31(1):345-347.

94. Gough J, Karplus K, Hughey R, Chothia C: **Assignment of homology to genome sequences using a library of hidden Markov models that represent all proteins of known structure.** J Mol Biol 2001, 313(4):903-919.

95. Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: **The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools.** Nucleic Acids Res 1997, 25(24):4876-4882.

96. Nei M, Chakraborty R: **Empirical relationship between the number of nucleotide substitutions and interspecific identity of amino acid sequences in some proteins.** J Mol Evol 1976, 7(4):313-323.

97. Saitou N, Nei M: **The neighbor-joining method: a new method for reconstructing phylogenetic trees.** Mol Biol Evol 1987, 4(4):406-425.

98. Kumar S, Tamura K, Nei M: **MEGA3: Integrated software for Molecular Evolutionary Genetics Analysis and sequence alignment.** Brief Bioinform 2004, 5(2):150-163.

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_10.png" /></figure>

BioMedcentral