BACKGROUND
saffron  is a triploid, sterile plant, probably derived from the wild species crocus cartwrightianus.
it has been propagated and used as a spice and medicinal plant in the mediterranean area for thousands of years  <cit> .
the domestication of saffron probably occurred in the greek-minoan civilization between  <dig>  and  <dig>  b.c.
a fresco depicting saffron gatherers, dating back to  <dig>  b.c.
has been unearthed on the island of santorini, greece.
saffron is commonly considered the most expensive spice on earth.
nowadays, the main producing countries are iran, greece, spain, italy, and india .
apart from the commercial and historical aspects, several other characteristics make saffron an interesting biological system: the spice is derived from the stigmas of the flower , which are harvested manually and subjected to desiccation.
the main colors of saffron, crocetin and crocetin glycosides, and the main flavors, picrocrocin and safranal, are derived from the oxidative cleavage of the carotenoid, zeaxanthin  <cit>  .
saffron belongs to the iridaceae  with poorly characterized genomes of relatively large size.
the characterization of the transcriptome of saffron stigmas is likely to shed light on several important biological phenomena: the molecular basis of flavor and color biogenesis in spices, the biology of the gynoecium, and the genomic organization of iridaceae.
for these reasons, we have undertaken the sequencing and bioinformatics characterization of expressed sequence tags  from saffron stigmas.
RESULTS
sequencing and assembly
an oriented cdna library from mature saffron stigmas in lambda uni-zap  <cit>  was kindly provided by prof. bilal camara, university of strasbourg.
the library was subjected to automated excision, and the cdna inserts were subjected to pcr amplification and sequenced from the 5' end.
<dig>  electropherograms were analyzed with the phred program  <cit> .
low quality sequences were removed from the 5' and 3' ends, and the sequences were further processed to remove vector contaminations and to mask low complexity and/or repeat sub-sequences.
this process reduced the original dataset to  <dig>  high-quality sequences longer than  <dig> nucleotides.
only  <dig>  est fragments whose length is greater than or equal to  <dig> nucleotides were considered for the submission to the ncbi dbest division.
they are accessible under the accession numbers from ex <dig> to ex <dig>
the est dataset was subjected to a clustering/assembling procedure  <cit> , in order to group ests putatively derived from the same gene and to generate a tentative consensus sequence  per putative transcript.
the total number of clusters generated are  <dig> .
each cluster should correspond to a unique gene, i.e.
it represents a gene index.
<dig>  clusters are made up of a single est and are therefore classified as singletons.
the remaining  <dig> clusters are made up of  <dig>  ests, assembled into  <dig> tcs .
in  <dig> clusters, ests are assembled so that multiple tcs are defined .
multiple tcs in a cluster have common regions of high similarity that may be due to possible alternative transcripts, to paralogy or to domain sharing.
the gc content distribution in the dataset is reported in figure  <dig>  the average gc content is around 44%.
the database and the web interface
the dataset was used to construct the saffron genes database  <cit> .
the database architecture consists of a main mysql relational database where all the data generated are deposited, and two satellite databases mygo and mykegg.
a user-friendly web interface is created using html and php scripts.
a pre-defined query system supports data retrieval; html-tree graphical display is implemented to browse enzyme classes and metabolic pathways.
transcripts, which correspond to criteria defined by the user, can be mapped on-the-fly onto the kegg metabolic maps, which are accessible as gif images  <cit> .
the electropherograms of the single ests can be downloaded to re-check sequence quality.
automated functional annotation
in order to assign a preliminary function to each transcript, the tcs and singletons were compared using blastx to the uniprotkb/swiss-prot database.
of  <dig>  transcripts,  <dig>   have no hits, while the remaining  <dig>  have at least one significant match in the protein database.
within this latter set,  <dig>  are described as hypothetical, unknown or expressed proteins thus not confirming an effective functional role of the transcript product.
gene ontology terms were assigned automatically to those  <dig> transcripts matching a protein in the uniprotkb/swiss-prot database whose accession numbers are present into the satellite database mygo .
in many cases, multiple gene ontology terms could be assigned to the same transcript, resulting in  <dig> assignments to the molecular function,  <dig> to the biological process and finally  <dig>  to the cellular component class.
to give a broad overview of the ontology content, the entire set of the ontologies was mapped onto the plant go slims terms.
in the molecular function ontology class, the most represented terms describe catalytic  and hydrolase activity  .
the remaining categories are less represented.
considering the biological process class, the vast majority of the go assignments corresponds to the more general transport category  .
finally, for the cellular component class the assignments were mainly given to the plastid , mitochondrion , and cytoplasmic membrane-bound vesicle  components .
<dig> transcripts are associated to  <dig> distinct enzymes as they are classified and described into the enzyme repository  <cit> .
<dig> out of the  <dig> enzymes had mappings to  <dig> kegg biochemical pathways  <cit> .
as we know, some enzymes can occur in more than one pathway; on the other hand there are  <dig> enzymes which only act in a single pathway, that were classified as pathway-specific .
genes expressed in crocus stigmas
est abundance in a contig can be indicative of the mrna relative abundance in the stigma tissue.
we identified the tcs that are composed of ≥  <dig> ests .
the most highly expressed tc, cl000057: <dig> , bears homology to short chain dehydrogenases .
this protein family comprises members involved in hormone biosynthesis, like the aba <dig> gene of arabidopsis which catalyzes the conversion of xanthoxin into aba aldehyde  <cit> , or in sexual organ identity, like the tasselseed <dig>  gene of maize .
ts <dig> is expressed in pistil primordia cells of maize, where it activates a cell death process eliminating these cells from male reproductive organs  <cit> .
biochemical studies suggest that the ts <dig> protein is a hydroxysteroid dehydrogenase  <cit> .
it will be interesting to determine the function and substrate specificity of the saffron cl000057: <dig> product.
a large number of cytochrome p <dig> sequences are expressed in saffron stigmas, some of which at very high levels .
also, lipid metabolism seems to be very active, judging from the tcs encoding proteins involved in this process .
several tcs encode putative carotenoid metabolism enzymes : cl000944: <dig> encodes non-heme -β-carotene-hydroxylase, which is highly expressed in saffron stigmas  <cit> .
cl000627: <dig> encodes a putative glucosyltransferase, very similar to ugtcs <dig>  which is able to glycosylate crocetin in vitro  <cit>  .
cl001532: <dig> and cl001032: <dig> also, encode putative isoprenoid gtases, one of which could represent the still missing enzyme responsible for the glycosylation of picrocrocin .
cl001432: <dig> encodes a protein similar to plastid terminal oxidase, involved in phytoene desaturation  <cit> , while est cr36_b <dig> encodes a protein similar to fibrillin, which is a carotenoid-binding protein in pepper chromoplasts  <cit> .
cl <dig> encodes a carboxyl methyltransferase very similar to the one catalyzing the synthesis of bixin  <cit>  .
this tc seems to encode a "short" form of the annatto and crocus methyltransferases from genbank, possibly derived from alternative splicing .
although a methyltransferase reaction has not been described in saffron stigmas, the biosynthesis of bixin and that of crocin share some features in common, since both pigments are derived from the oxidative cleavage of a carotenoid  <cit> .
finally, cl000045: <dig> encodes a protein highly similar to the cauliflower or gene product, a plastid-associated protein with a cysteine-rich dnaj domain.
a dominant or mutation induces β-carotene accumulation in cauliflower inflorescences, suggesting that or is somehow involved in the control of chromoplast differentiation  <cit> .
several tcs encode putative transcription factors .
the most abundantly expressed, cl000348: <dig>  encodes a myb-like protein with high similarity to lhmyb  myb <dig>  and myb <dig> .
all three factors are highly expressed in flowers.
also highly expressed is cl001329: <dig>  encoding a putative mads box transcription factor.
this protein shows high similarity to aodef, a b-functional transcription factor from asparagus expressed in stamens and inner tepals  <cit>  and to lmads <dig>  a lily protein whose ectopic expression in dominant negative form causes an ap3-like phenotype in arabidopsis  <cit> .
finally, several tcs – cl000209: <dig>  cl000582: <dig>  cl001827: <dig>  and cl <dig> – show similarity to potyviral sequences, indicating that the sequenced library likely derives from virus-infected tissue.
potyviruses like iris mild mosaic virus are known to infect crocus  <cit> .
the sequences of these tcs will prove useful for diagnostic and phytosanitary purposes.
CONCLUSIONS
the saffron genes database  <cit>  has been designed to manage and to explore the est collection from saffron stigmas, providing a reference for the expression pattern analysis in this tissue as well as a primary view of the genomic properties of this species, representative of iridaceae.
the complete set of raw est sequences, as well as of their electopherograms, are maintained in the database allowing users investigate on library qualities and on single est structural features .
annotation is provided for single ests as well as for their assemblies , to evaluate the consistency of the automated functional assignments.
the putative transcripts determined to be associated to enzymes are organized into classes and can be viewed also in terms of enzyme assignments to metabolic pathways.
this represents a straightforward way to investigate the properties of the stigma transriptome.
as discussed above, this transcriptome contains a series of interesting sequences, whose function can now be tested using in vivo or in vitro approaches.

 saffron is a triploid, sterile plant, probably derived from the wild species crocus cartwrightianus. it has been propagated and used as a spice and medicinal plant in the mediterranean area for thousands of years. the domestication of saffron probably occurred in the greek-minoan civilization between and b.c. a fresco depicting saffron gatherers, dating back to b.c. has been unearthed on the island of santorini, greece. saffron is commonly considered the most expensive spice on earth. nowadays, the main producing countries are iran, greece, spain, italy, and india. apart from the commercial and historical aspects, several other characteristics make saffron an interesting biological system: the spice is derived from the stigmas of the flower, which are harvested manually and subjected to desiccation. the main colors of saffron, crocetin and crocetin glycosides, and the main flavors, picrocrocin and safranal, are derived from the oxidative cleavage of the carotenoid, zeaxanthin. saffron belongs to the iridaceae with poorly characterized genomes of relatively large size. the characterization of the transcriptome of saffron stigmas is likely to shed light on several important biological phenomena: the molecular basis of flavor and color biogenesis in spices, the biology of the gynoecium, and the genomic organization of iridaceae. for these reasons, we have undertaken the sequencing and bioinformatics characterization of expressed sequence tags from saffron stigmas.  sequencing and assembly an oriented cdna library from mature saffron stigmas in lambda uni-zap was kindly provided by prof. bilal camara, university of strasbourg. the library was subjected to automated excision, and the cdna inserts were subjected to pcr amplification and sequenced from the 5' end. electropherograms were analyzed with the phred program. low quality sequences were removed from the 5' and 3' ends, and the sequences were further processed to remove vector contaminations and to mask low complexity and/or repeat sub-sequences. this process reduced the original dataset to high-quality sequences longer than nucleotides. only est fragments whose length is greater than or equal to nucleotides were considered for the submission to the ncbi dbest division. they are accessible under the accession numbers from ex to ex the est dataset was subjected to a clustering/assembling procedure, in order to group ests putatively derived from the same gene and to generate a tentative consensus sequence per putative transcript. the total number of clusters generated are. each cluster should correspond to a unique gene, i.e. it represents a gene index. clusters are made up of a single est and are therefore classified as singletons. the remaining clusters are made up of ests, assembled into tcs. in clusters, ests are assembled so that multiple tcs are defined. multiple tcs in a cluster have common regions of high similarity that may be due to possible alternative transcripts, to paralogy or to domain sharing. the gc content distribution in the dataset is reported in figure the average gc content is around 44%. the database and the web interface the dataset was used to construct the saffron genes database. the database architecture consists of a main mysql relational database where all the data generated are deposited, and two satellite databases mygo and mykegg. a user-friendly web interface is created using html and php scripts. a pre-defined query system supports data retrieval; html-tree graphical display is implemented to browse enzyme classes and metabolic pathways. transcripts, which correspond to criteria defined by the user, can be mapped on-the-fly onto the kegg metabolic maps, which are accessible as gif images. the electropherograms of the single ests can be downloaded to re-check sequence quality. automated functional annotation in order to assign a preliminary function to each transcript, the tcs and singletons were compared using blastx to the uniprotkb/swiss-prot database. of transcripts, have no hits, while the remaining have at least one significant match in the protein database. within this latter set, are described as hypothetical, unknown or expressed proteins thus not confirming an effective functional role of the transcript product. gene ontology terms were assigned automatically to those transcripts matching a protein in the uniprotkb/swiss-prot database whose accession numbers are present into the satellite database mygo. in many cases, multiple gene ontology terms could be assigned to the same transcript, resulting in assignments to the molecular function, to the biological process and finally to the cellular component class. to give a broad overview of the ontology content, the entire set of the ontologies was mapped onto the plant go slims terms. in the molecular function ontology class, the most represented terms describe catalytic and hydrolase activity. the remaining categories are less represented. considering the biological process class, the vast majority of the go assignments corresponds to the more general transport category. finally, for the cellular component class the assignments were mainly given to the plastid, mitochondrion, and cytoplasmic membrane-bound vesicle components. transcripts are associated to distinct enzymes as they are classified and described into the enzyme repository. out of the enzymes had mappings to kegg biochemical pathways. as we know, some enzymes can occur in more than one pathway; on the other hand there are enzymes which only act in a single pathway, that were classified as pathway-specific. genes expressed in crocus stigmas est abundance in a contig can be indicative of the mrna relative abundance in the stigma tissue. we identified the tcs that are composed of ≥ ests. the most highly expressed tc, cl000057:, bears homology to short chain dehydrogenases. this protein family comprises members involved in hormone biosynthesis, like the aba gene of arabidopsis which catalyzes the conversion of xanthoxin into aba aldehyde, or in sexual organ identity, like the tasselseed gene of maize. ts is expressed in pistil primordia cells of maize, where it activates a cell death process eliminating these cells from male reproductive organs. biochemical studies suggest that the ts protein is a hydroxysteroid dehydrogenase. it will be interesting to determine the function and substrate specificity of the saffron cl000057: product. a large number of cytochrome p sequences are expressed in saffron stigmas, some of which at very high levels. also, lipid metabolism seems to be very active, judging from the tcs encoding proteins involved in this process. several tcs encode putative carotenoid metabolism enzymes: cl000944: encodes non-heme -β-carotene-hydroxylase, which is highly expressed in saffron stigmas. cl000627: encodes a putative glucosyltransferase, very similar to ugtcs which is able to glycosylate crocetin in vitro. cl001532: and cl001032: also, encode putative isoprenoid gtases, one of which could represent the still missing enzyme responsible for the glycosylation of picrocrocin. cl001432: encodes a protein similar to plastid terminal oxidase, involved in phytoene desaturation, while est cr36_b encodes a protein similar to fibrillin, which is a carotenoid-binding protein in pepper chromoplasts. cl encodes a carboxyl methyltransferase very similar to the one catalyzing the synthesis of bixin. this tc seems to encode a "short" form of the annatto and crocus methyltransferases from genbank, possibly derived from alternative splicing. although a methyltransferase reaction has not been described in saffron stigmas, the biosynthesis of bixin and that of crocin share some features in common, since both pigments are derived from the oxidative cleavage of a carotenoid. finally, cl000045: encodes a protein highly similar to the cauliflower or gene product, a plastid-associated protein with a cysteine-rich dnaj domain. a dominant or mutation induces β-carotene accumulation in cauliflower inflorescences, suggesting that or is somehow involved in the control of chromoplast differentiation. several tcs encode putative transcription factors. the most abundantly expressed, cl000348: encodes a myb-like protein with high similarity to lhmyb myb and myb. all three factors are highly expressed in flowers. also highly expressed is cl001329: encoding a putative mads box transcription factor. this protein shows high similarity to aodef, a b-functional transcription factor from asparagus expressed in stamens and inner tepals and to lmads a lily protein whose ectopic expression in dominant negative form causes an ap3-like phenotype in arabidopsis. finally, several tcs – cl000209: cl000582: cl001827: and cl – show similarity to potyviral sequences, indicating that the sequenced library likely derives from virus-infected tissue. potyviruses like iris mild mosaic virus are known to infect crocus. the sequences of these tcs will prove useful for diagnostic and phytosanitary purposes.  the saffron genes database has been designed to manage and to explore the est collection from saffron stigmas, providing a reference for the expression pattern analysis in this tissue as well as a primary view of the genomic properties of this species, representative of iridaceae. the complete set of raw est sequences, as well as of their electopherograms, are maintained in the database allowing users investigate on library qualities and on single est structural features. annotation is provided for single ests as well as for their assemblies, to evaluate the consistency of the automated functional assignments. the putative transcripts determined to be associated to enzymes are organized into classes and can be viewed also in terms of enzyme assignments to metabolic pathways. this represents a straightforward way to investigate the properties of the stigma transriptome. as discussed above, this transcriptome contains a series of interesting sequences, whose function can now be tested using in vivo or in vitro approaches. 
the complete set of raw est sequences, as well as of their electopherograms, are maintained in the database, allowing users to investigate sequence qualities and est structural features. saffron flowers have been used as a spice and medicinal plant ever since the greek-minoan civilization. the saffron genes database represents the first reference collection for the genomics of iridaceae, for the molecular biology of stigma biogenesis, as well as for the metabolic pathways underlying saffron secondary metabolism. we produced high quality expressed sequence tags from a saffron stigma cdna library. the saffron stigma transcriptome contains a series of interesting sequences.
 saffron flowers have been used as a spice and medicinal plant ever since the greek-minoan civilization. the edible part – the stigmas – are commonly considered the most expensive spice in the world and are the site of a peculiar secondary metabolism, responsible for the characteristic color and flavor of saffron.  we produced high quality expressed sequence tags from a saffron stigma cdna library. this collection is accessible and searchable through the saffron genes database http://www.saffrongenes.org. the ests have been grouped into clusters, each corresponding to a different expressed gene, and annotated. the complete set of raw est sequences, as well as of their electopherograms, are maintained in the database, allowing users to investigate sequence qualities and est structural features. the saffron stigma transcriptome contains a series of interesting sequences.  the saffron genes database represents the first reference collection for the genomics of iridaceae, for the molecular biology of stigma biogenesis, as well as for the metabolic pathways underlying saffron secondary metabolism. 