<title>Background</title> <p>The genus <italic>Actinidia </italic>Lindl. is large, containing between 50 and 70 species of climbing plants originating mainly in southern China [##UREF##0##1##]. Over the past 30 years kiwifruit has developed into an important horticultural crop, firstly in New Zealand, and subsequently in other countries such as Chile, China and Italy [##UREF##1##2##,##UREF##2##3##]. Currently cultivars from three species are grown commercially; the green-fleshed kiwifruit, <italic>Actinidia deliciosa </italic>(A. Chev.) C.F.Liang et A.R.Ferguson, the closely related yellow-fleshed <italic>A. chinensis </italic>Planch. and the kiwiberry, <italic>A. arguta </italic>(Sieb. et Zucc.) Planch. ex Miq. Most of the kiwifruit cultivars grown commercially are seedling selections and there has been little systematic breeding [##UREF##1##2##]. Consequently there are still many characteristics within the genus (Fig. ##FIG##0##1A##) that could be incorporated into commercial cultivars [##UREF##0##1##] and to do this efficiently requires a better knowledge of how these characteristics are regulated. For the consumer, critical desirable attributes of kiwifruit are flavor and fragrance, appearance, healthful components and convenience.</p> <p>Flavor and fragrance are determined by the chemical composition of the fruit. A wide range of compounds has been detected within the <italic>Actinidia </italic>genus leading to distinct and different flavors in the fruit [##UREF##3##4##]. These compounds include polyphenolics, acids, alcohols and volatile compounds such as terpenes and esters. Flavor is also influenced by the sugar to acid balance, with some growers currently receiving a premium for sweeter fruit. <italic>Actinidia </italic>fruit color ranges from the original green kiwifruit, newer yellow varieties as well as red, purple and orange fruit (Fig. ##FIG##0##1B##). The extensive range of colors is caused by the presence or absence of chlorophyll, anthocyanins and carotenoids [##UREF##3##4##]. Healthful attributes of kiwifruit include its high ascorbic acid levels [##UREF##1##2##], quinic acid levels [##REF##15619581##5##], and the presence of triterpenes and folic acid. On the other hand, allergens are identified as undesirable health attributes in kiwifruit and in many other fruit [##UREF##3##4##]. Convenience includes such factors as the ability to store the fruit for a long time [##UREF##1##2##] and a long shelf life, as well as easy determination of eating ripeness and having an edible or peelable skin. Little is currently known about the genetic diversity and gene-based control of these major attributes in kiwifruit.</p> <p>EST databases as resources for uncovering genetic diversity have been established for many plant species (Additional file ##SUPPL##0##1##). Most comprehensively surveyed are <italic>Arabidopsis </italic>(<italic>Arabidopsis thaliana</italic>; 1,279,945 in GenBank) and rice (<italic>Oryza sativa</italic>; 1,211,418 in GenBank), both of which have been fully genome-sequenced. Fruit crops have been less extensively surveyed, but recently there have been reports on EST projects from fruit, including tomato (<italic>Lycopersicon esculentum </italic>[##REF##15361140##6##], 258,408 sequences in GenBank), grape (<italic>Vitis vinifera</italic>, which has also had its genome fully sequenced [##REF##17721507##7##]; [##REF##17761391##8##,##REF##15856347##9##], 363,365 ESTs in GenBank), apple (<italic>Malus </italic>× <italic>domestica </italic>[##REF##16531485##10##], 255,103 sequences in GenBank) and pineapple (<italic>Ananas comosus </italic>[##REF##15520025##11##], 5649 ESTs in GenBank). However, there are few ESTs sequenced from the Ericales (the order to which <italic>Actinidia </italic>belongs) registered in GenBank (Additional file ##SUPPL##0##1##). These significant fruit EST resources have been used to identify genes likely to be involved in the ripening process in tomato [##REF##15361140##6##] and in the generation of aroma in apple [##REF##17556515##12##]. In addition, ESTs are useful sources of simple sequence repeats (SSRs) and single-nucleotide polymorphisms (SNPs), both of which are useful markers for creating genetic maps in plants [##REF##12647067##13##, ####REF##15067386##14##, ##REF##16211378##15##, ##REF##11999831##16##, ##REF##16928278##17##, ##REF##12579429##18##, ##REF##15197579##19####15197579##19##].</p> <p>Knowledge of the <italic>Actinidia </italic>genome or transcriptome is currently restricted to the 511 sequences available in GenBank (dbEST, January 2008). To improve and develop new varieties of kiwifruit with the desired flavor, health and convenience attributes, it is useful to use modern genomics techniques in conjunction with breeding tools such as marker-assisted selection and genetic transformation. For this reason, we undertook a major EST sequencing project in <italic>Actinidia </italic>to develop a basic genetic resource covering a range of different species and tissues. In this paper we analyze more than 130,000 ESTs, derived mainly from four species in the <italic>Actinidia </italic>genus (Table ##TAB##0##1##). This information is presented with reference to the chemical composition of these species and using specific examples where we have increased our understanding of the genetic diversity and gene-based control of critical attributes.</p>
<title>Methods</title> <title>Plant material</title> <p>Tissues were collected from <italic>Actinidia </italic>species growing in New Zealand at HortResearch research orchards in Auckland, Bay of Plenty (Te Puke) and Northland (Kerikeri) from 1999 to 2003. Table ##TAB##0##1## and Additional file ##SUPPL##1##2a## provide details of tissues, species and treatments.</p> <title>Library construction and EST sequencing</title> <p>Total RNA was extracted from <italic>Actinidia </italic>tissues by established methods [##UREF##34##87##,##UREF##35##88##]. Messenger RNA was isolated from total RNA by passage through oligo(dT)-cellulose columns (GE Healthcare, USA), and cDNA cloned into either phage (ZapcDNA Synthesis Kit and Zap-cDNA Gigapack III Gold Cloning Kit; Stratagene, USA) or plasmid-based libraries (SuperScript System for cDNA Synthesis and Cloning; Invitrogen, USA). In some libraries, subtractive or other enrichment techniques were used (Additional file ##SUPPL##1##2b##).</p> <p>Plasmids from the phage cDNA libraries were mass excised, according to the manufacturer's recommendations (Stratagene). Plasmid extractions were then undertaken on individual bacterial colonies of either the phage-derived or the plasmid-derived cDNA libraries and the corresponding cDNA inserts sequenced predominantly from the 5' end. Big Dye Terminator sequencing reactions were resolved on ABI377, ABI3100, or ABI3700 sequencers, according to the manufacturer's instructions (Applied Biosystems, USA). For determination of the complete sequence of cDNA clones, M13R and M13F or T3 and T7 primers were used for 5' and 3' end sequencing. EST-specific primers were used to determine the complete sequence of cDNA clones. In situations where EST clones had long poly(A) tails (generally 40 nucleotides) and, therefore, failed to yield good quality sequence with standard sequencing primers, an anchored T24VN primer was used. Resulting sequences were edited manually and assembled using Sequencher software, version 4.0.5 (GeneCodes). Sequencing progress for each cDNA library was assessed manually for clone length, redundancy and sequence quality.</p> <title>Bioinformatics</title> <p>EST sequences were automatically trimmed of vector, adapter, and low quality sequence regions, and uploaded to a relational database. Automatic annotation was performed using the HortResearch BioView sequence annotation pipeline (BioView – an enterprise bioinformatics system for automated analysis and annotation of non-genomic DNA sequence (Crowhurst R, Davy M, Deng C, unpublished)) that utilizes a relational database (MySQL; <ext-link ext-link-type="uri" xlink:href="http://www.mysql.com"/>). The EST clustering phase was performed using The Institute for Genomic Research (TIGR) gene indices clustering tools <ext-link ext-link-type="uri" xlink:href="http://www.tigr.org/tdb/tgi/software"/>. The representation of protein families, domains, and functional sites within the <italic>Actinidia </italic>NRs was determined using InterProScan [##UREF##36##89##]. The proteome for <italic>Arabidopsis </italic>(<italic>Arabidopsis thaliana</italic>) was obtained from TAIR [##REF##12444417##90##], and comparisons to proteins from <italic>Arabidopsis </italic>using BLASTx were used to identify <italic>Actinidia </italic>NRs with similarity to <italic>Arabidopsis </italic>proteins.</p> <p>Detection of SSRs was undertaken using a PERL program within BioView that identified tandem repetition of sequence words in target sequences. SSRs were characterized by repeat type (di-, tri-, or tetra-nucleotide repeat units), repeat length, and position. Only repeats longer than 10 bp were included in the analyses. When reporting the frequency of repeat classes, different di- and tri-nucleotide sequences were combined by type; for example, AG repeats also encompassed repeats identified as GA and their complementary sequences CT or TC repeats. Prediction of SNPs and insertion/deletions and sequencing errors was performed using PERL scripts within BioView that parsed the output of contig sequences generated by the CAP3 DNA sequence assembly program [##REF##10508846##91##].</p> <p>Codon usage tables were derived from cDNA sequences encoding predicted full-length proteins. Clones were predicted to be full length only if they started with an ATG codon and terminated with a stop codon at positions equivalent to those of other plant genes. Codon usage was calculated from sequences using the CUSP program implemented within EMBOSS [##REF##10827456##92##]. For functional analysis using Mapman [##REF##14996223##27##], each <italic>Actinidia </italic>NR was assigned the nearest <italic>Arabidopsis </italic>BLASTx match (E &lt; 1.0e<sup>-10</sup>) if there was no conflicting evidence from domains or families detected by InterProScan. The NR and its EST members were then assigned a Mapman classification bin and bin name based on that assigned to the <italic>Arabidopsis </italic>match.</p> <title>Homology searching and phylogenetic analysis</title> <p>Previously published protein sequences from GenBank were used to identify genes in the <italic>Actinidia </italic>EST database using tBLASTn searches [##REF##2231712##93##], typically with a cutoff value of E &lt; 1.0e<sup>-20</sup>. Identified genes were then manually checked to ensure accuracy. Amino acid alignments of predicted proteins were constructed using Clustal X. Phylogenetic analysis was carried out using the PHYLIP suite of programs [##UREF##37##94##]. Distances were calculated using protdist, and the Fitch method was used to construct a tree. Bootstrap analysis was conducted using 1000 bootstrap replicates using seqboot [##UREF##37##94##]. Treeview (v.1.6.6) was used to display resulting trees [##REF##8902363##95##].</p> <title>Microarray construction</title> <p>For each predicted NR in the <italic>Actinidia </italic>EST database, a 45–55 mer oligonucleotide was created, using PERL scripts within BioView that selected all possible oligonucleotides within a gene based on uniqueness of sequence, lack of repetitive regions, and a constant melting temperature. For each NR the most 3' oligonucleotide that passed these criteria was selected for synthesis. As the EST libraries were sequenced from different <italic>Actinidia </italic>species, the resulting oligonucleotides were derived from genes from different <italic>Actinidia </italic>species. A total of 17,472 oligonucleotides were made, with the majority coming from <italic>A. chinensis </italic>and <italic>A. deliciosa</italic>. A summary of the oligonucleotides used to construct the <italic>Actinidia </italic>microarray slides by library and species is given in Additional file ##SUPPL##3##4##. Genomic DNA from <italic>A. deliciosa </italic>'Hayward' and <italic>A. eriantha </italic>(genotype 11-6-15e) was labeled and each microarray hybridized as described previously [##REF##17556515##12##]. Data were normalized using quantile normalization in the Bioconductor package Limma [##REF##14597310##96##].</p> <title>Chemical analysis of actinidia flowers and fruit</title> <title>Volatile flavor and aroma compounds</title> <p>Both headspace and solvent sampling was used to obtain volatiles. For headspace sampling of flowers, whole flowers (2.5–5 g) that were 50–75% open were placed in a 250 ml Quickfit™ Erlenmeyer flask to which was fitted a volatile trapping cartridge packed with 100 mg of Chromosorb 105™ absorbent and sampled and analyzed as described earlier [##REF##12737978##29##]. For fruit, ten fruit (eating-ripe) were cut longitudinally in half and one half of each fruit was used for headspace sampling and analysis [##REF##12737978##29##]. The other halves were frozen in liquid N<sub>2 </sub>for solvent extraction.</p> <p>For solvent extraction, either ~1 g of flower petals were rinsed in 2 ml of purified 50:50 pentane/Et<sub>2</sub>O or 20 g of fruit was pulped, and gently shaken in a sealed test tube with 10 mL of pentane:Et<sub>2</sub>O (8:2) for 1–2 min, four times over a 2 h period at room temperature. Samples were processed and analyzed by GC-MS as described earlier [##REF##12737978##29##].</p> <p>Quantification of compounds was carried out using an average FID response based on methyl butanoate, ethyl butanoate, hexanol and methyl benzoate. Component identification was by comparison with spectra in the Mass Spectral Database (1998 NIST and an in-house database), retention indices (in-house database) and in some cases direct GC-MS comparison with authentic standards.</p> <title>Analysis of fruit acids and sugars</title> <p>Tissue was ground in liquid N<sub>2 </sub>and acids and sugars were extracted from a known mass of tissue (1.4–2.5 g) into 10 mL of 80% EtOH (with adonitol and tartaric acid added as internal standards) at 60°C for 1 h. The extracts were processed, sugars and organic acids derivatized and then quantified by GC [##UREF##38##97##]. Identification of the acids and sugars was confirmed by GC-MS comparison with mixtures of authentic compounds.</p> <title>Sequence data</title> <p>Sequence data for the ESTs described in this paper can be found at GenBank under accession numbers <ext-link ext-link-type="gen" xlink:href="FG396013">FG396013</ext-link> – <ext-link ext-link-type="gen" xlink:href="FG528589">FG528589</ext-link>.</p>
<title>Results</title> <title>Overview</title> <p><italic>Actinidia </italic>cDNA libraries were constructed mainly from four <italic>Actinidia </italic>species: <italic>A. deliciosa, A. chinensis, A. eriantha </italic>and <italic>A. arguta </italic>(Table ##TAB##0##1##, Additional file ##SUPPL##1##2##). The libraries were constructed from petals, fruits, buds and leaves, with a small number from roots and cell culture and sequenced from the 5' end. The average edited sequence length of the 132,577 ESTs was 503 bases. Clustering these sequences using a 95% threshold resulted in 18,070 tentative consensus (TC) sequences (average length 577 bases), with 23,788 ESTs remaining as singletons. The combination of TCs and singletons are referred to as non redundant (NR) clusters. These numbers are similar to those observed with apple [##REF##16531485##10##]. Over 95% of TCs had fewer than 18 EST members, and the largest single TC had 758 members. It should be noted that <italic>A. chinensis </italic>and <italic>A. deliciosa </italic>are more closely related to each other than the other species mentioned in this paper, and that only a restricted number of genotypes within each species was used in making EST libraries compared to the total numbers available.</p> <title>Sequence analysis</title> <title>Codon usage</title> <p>Knowledge of the GC content of a genome and codon usage is useful when devising PCR-based strategies for mapping and gene isolation, as well as for hybridization studies by microarray. Based on the coding regions of 302 <italic>A. deliciosa</italic>, 319 <italic>A. chinensis </italic>and 84 <italic>A. eriantha </italic>full-length cDNAs, the GC content in the third base position was estimated to be 55%, 49% and 58%, respectively. These values are higher than the overall GC ratio from the sequences of the complete EST datasets of <italic>A. deliciosa </italic>(46%), <italic>A. chinensis </italic>(46%) and <italic>A. eriantha </italic>(48%), indicating that there is some pressure, particularly in <italic>A. deliciosa </italic>and <italic>A. eriantha</italic>, towards an increased GC ratio in coding regions compared with non-translated regions. Overall the codon usage of the three <italic>Actinidia </italic>species is similar, although not identical (Additional file ##SUPPL##2##3##). <italic>A. deliciosa </italic>and <italic>A. eriantha </italic>differ in only their preference for aspartate (GAT/GAC) and serine (TCT/TCC). In contrast, comparative codon preferences between <italic>A. deliciosa </italic>and <italic>A. chinensis </italic>show differences for alanine (GCC/GCT), glycine (GGC/GGA), isoleucine (ATC/ATT), leucine (CTC/TTG), asparagine (AAC/AAT), glutamine (CAG/CAA), threonine (ACC/ACT) and valine (GTG/GTT). In most of these codons, <italic>A. chinensis </italic>shows a greater preference for an A or T in the third base position, accounting for its lower GC third base percentage when compared with the other two species.</p> <title>Polymorphisms and genetic markers</title> <p>Both SNPs and microsatellites or SSRs are valuable tools for genetic mapping within breeding populations of many crops. A total of 32,764 biallelic SNPs were detected from the overlapping regions of 3,901 (21.6%) of the 18,070 TCs, at a rate of one SNP every 417 bp. This frequency is higher than that reported in apple using a dataset of comparable size (one SNP per 706 bp from ~150,000 ESTs [##REF##16531485##10##]). As several <italic>Actinidia </italic>species were used to construct the cDNA libraries, while the apple cDNA libraries were from different varieties of the same species, this increased frequency of SNPs could reflect the greater genetic diversity in the surveyed <italic>Actinidia </italic>EST libraries. For example the frequency in apple increased significantly to 1 in 149 bp when a greater diversity of genotypes was added (Chagne <italic>et al </italic>unpublished results). The polyploidy nature of several of the <italic>Actinidia </italic>species used may also have affected SNP frequencies. As a result of this, some of the SNPs identified may not be allelic in nature but due to homoeologous or paralogous sequences clustering in the same TC.</p> <p>EST-derived SSRs have already been shown to be valuable mapping tools in <italic>Actinidia chinensis </italic>[##REF##15067386##14##] where 150 SSRs with more than 10 dinucleotide repeat units were tested as markers in an <italic>A. chinensis </italic>mapping population. More than 90% of the SSR markers were polymorphic and segregated within the population. Subsequently, 20 of these SSRs were shown to be transportable across multiple species of <italic>Actinidia</italic>, showing the value of this resource [##REF##16211378##15##].</p> <p>Compared with apple [##REF##16531485##10##], where less than 20% of the NRs contained a microsatellite, over 30% of <italic>Actinidia </italic>NRs were found to have at least one. The other major difference between apple and <italic>Actinidia </italic>was that, while di- and tri-nucleotide repeats were equally frequent in apple at 7 to 8% frequency, di-nucleotide repeats (18%) were twice as frequent as tri-nucleotide repeats (9%) in <italic>Actinidia</italic>. Tetra-nucleotide repeats were found at similar frequencies in NRs for the two genera. The lengths of the SSRs were similar in both apple and <italic>Actinidia</italic>, with 50% of the repeats having between 12 and 14 bases. As in apple, AG repeats were the most frequent (16% of NRs, significantly more frequent than in apple) and were followed by AT (1.2%) and AC (0.8%) while GC repeats were very infrequent.</p> <p>The position of SSRs in relation to the putative initiation ATG was very similar when comparing <italic>Actinidia </italic>with other plants such as <italic>Arabidopsis </italic>and apple [##REF##16531485##10##,##REF##14764542##20##], with di-nucleotide repeats being predominantly (93%) in the 5' UTR, and tri-nucleotide repeats being more evenly distributed along the gene (37%).</p> <title>Comparative genomic DNA hybridizations by microarray</title> <p>A 17,472 feature <italic>Actinidia </italic>oligonucleotide microarray was constructed based on sequence data from the <italic>Actinidia </italic>EST database. The 44–55 mer oligonucleotides were derived from gene sequences of different <italic>Actinidia </italic>species (Additional file ##SUPPL##3##4##). To test the cross hybridization between the different genomes of the <italic>Actinidia </italic>species, genomic DNA from two <italic>Actinidia </italic>species (<italic>A. deliciosa </italic>and <italic>A. eriantha</italic>) was hybridized to two different microarray slides. A plot of average signal intensity for the two species against each other is shown in Fig. ##FIG##1##2##. Of the 13,443 informative features, the majority (98.4%) showed a similar level of hybridization between the two species (less than an average of 2-fold change), with 172 showing a higher signal intensity in <italic>A. deliciosa </italic>(plus symbols, Fig. ##FIG##1##2##), and 44 showing a higher intensity in <italic>A. eriantha </italic>(triangles, Fig. ##FIG##1##2##). There was no correlation between the species from which the oligonucleotide was derived and the increased or decreased level of hybridization (Table ##TAB##1##2##). Consequently, the increases are probably due to separate gene duplications in each species.</p> <title>Functional analysis</title> <title>Highly populated tentative consensus sequences</title> <p>The TC with the greatest number of ESTs was a cysteine proteinase (actinidin family, 758 members) followed by a metallothionein (547 members) (Additional file ##SUPPL##4##5##). However, when similar proteins were taken into account (i.e. different TCs that all matched a common <italic>Arabidopsis </italic>protein as their closest homolog with E &lt; 1.0e<sup>-100</sup>), there were 1,266 ESTs matching cysteine proteinases followed by 769 metallothionein protein-encoding ESTs.</p> <p>It was not surprising to identify the cysteine protease actinidin as the most highly abundant TC cluster in <italic>Actinidia</italic>. Actinidin can constitute up to 50% of soluble protein in mature <italic>A. deliciosa </italic>fruit at harvest [##UREF##4##21##]. Recent research has shown that this cysteine protease exists in basic and acidic forms and the amount of each form varies between species, with fruit of the important <italic>A. chinensis </italic>cultivar 'Hort16A' being almost devoid of acidic actinidin [##UREF##5##22##]. ESTs for actinidin occur throughout the fruit libraries that make up the database, and these cluster into 10 distinct genes/alleles [##UREF##5##22##]. As two of the key fruit libraries were subtracted for actinidin, the number of ESTs in the TC is likely to be an under-representation of the actual number. The function of this protease may be related to insect defense [##UREF##6##23##].</p> <p>The metallothionein TC cluster contains a gene (pKIWI504) previously identified in kiwifruit [##REF##8075403##24##] as being highly expressed in young fruit with reduced expression in the later stages of fruit development. However, the <italic>in vivo </italic>function of the metallothionein TC cluster is currently unknown. Interestingly, pineapple also had a high proportion of metallothionein ESTs [##REF##15520025##11##].</p> <p>Fourteen TCs in <italic>Actinidia </italic>with their proteins annotated by BiolView as "not assigned-unknown" proteins (i.e. matches to sequences with unknown function from <italic>Arabidopsis </italic>or other species) were highly expressed with over 100 EST members. Four of the unassigned TCs contained more than 300 EST members. These four TCs included kiwellin [##REF##16328735##25##], which is a homolog to a grape ripening-related protein [##REF##10712544##26##] with no <italic>Arabidopsis </italic>homologs (527 EST members); a homolog of At5g11420, a protein of unknown function (339 EST members); a senescence-related protein (no <italic>Arabidopsis </italic>matches, 311 EST members); and an ABA stress-related protein (no <italic>Arabidopsis </italic>matches, 302 EST members).</p> <title>Functional analysis using Mapman and InterPro</title> <p>Functional analysis of ESTs and NRs using Mapman [##REF##14996223##27##] is shown in Table ##TAB##2##3##. The analysis was conducted by identifying each NR's nearest homolog to <italic>Arabidopsis</italic>, using BLASTx and identifying domains/families that are indicative for a given function using InterProScan. Then that NR and its member ESTs were assigned to the corresponding bin of the <italic>Arabidopsis </italic>match. Of the 41,858 NRs, 28,345 had sufficient homology to an <italic>Arabidopsis </italic>sequence to assign it a functional classification in Mapman. Thus, 32% of NRs from <italic>Actinidia </italic>had no <italic>Arabidopsis </italic>homolog, but a proportion would have homologs in other crops (e.g. genes related to fruit ripening).</p> <p>A comparison of the percentage of NRs assigned to the top level Mapman classification bins and the percentage of <italic>Arabidopsis </italic>genes classified in the same bin showed a strong relationship with an r<sup>2 </sup>value of 0.96 and a slope of 0.87, near to the expected 1:1 slope. The same relationship when calculated for second level bins gave an r<sup>2 </sup>of 0.98 and a slope of 1.14 (second level bin data not shown). Taking only those top level bins that contained more than 100 NRs to remove chance as a major factor, only two bins appeared to be over-represented (i.e. &gt; 2.0 times higher than the mean) in <italic>Actinidia</italic>, namely bin 4 (glycolysis) and bin 8 (TCA cycle). This may reflect a high level of primary metabolism in fruit. Interestingly, when all NRs are included, bin 5 (fermentation) and bin 6 (gluconeogenesis) were added to the list of over-represented NRs. The same relationships for potato showed no anomalies (Table ##TAB##2##3##). This analysis demonstrates that the clustering of ESTs generated a fair representation of the distribution of genes.</p> <p>Overall, the average number of ESTs per NR in the database was 4.0. There were only two bins with a high number of ESTs per NR, bin 15 (metal handling) with over 17 ESTs per NR and bin 22 (polyamine metabolism) with over 9. However, bin 22 had a low number of NR and EST members, suggesting this result may have occurred by chance. In contrast, bin 15 included the metallothioneins and the very highly expressed EST discussed earlier (Additional file ##SUPPL##4##5##).</p> <p>TCs with high numbers of ESTs (&gt; 100) came from 18 of the 35 Mapman bins [##REF##14996223##27##] with approximately three quarters falling into the bin names "not assigned, no ontology", "protein degradation", "metal handling binding chelation and storage", "stress biotic", "cell wall degradation" and "secondary metabolism isoprenoids".</p> <p>The most common InterPro families [##REF##17202162##28##] were also analyzed (Additional file ##SUPPL##5##6##). There were 3,111 InterPro families (flagged as true by InterPro) represented in the <italic>Actinidia </italic>NRs, with the most NRs found in the protein kinases grouping. Following that were thioredoxins, while cysteine proteinases, which had the highest number of ESTs in the database, were represented as the 7<sup>th </sup>most frequently occurring NR.</p> <title>Genes related to flavor and fragrance</title> <p>Understanding the relationship between fruit flavor and the genetic diversity present in the EST database requires a detailed analysis of the flavor and fragrance compounds present in the <italic>Actinidia </italic>genus. Fruit and flower samples from the major <italic>Actinidia </italic>varieties that were used to produce the EST database were characterized using both headspace sampling and solvent extraction followed by GC-FID/MS analysis. The results of this compound analysis are presented in Additional file ##SUPPL##6##7## for esters (over 55 identified compounds), acids, alcohols and terpenoids.</p> <title>Compound analysis</title> <p>The most abundant volatile compounds in the flowers of these <italic>Actinidia </italic>species were alcohols and terpenes. <italic>A. chinensis </italic>flowers contained 71% sesquiterpenes, largely dominated by farnesols (Additional file ##SUPPL##6##7##). <italic>A. deliciosa </italic>flowers had the highest volatile ester content at around 12%, and aromatic (6–9%) and straight-chain alcohols (9%). These flowers were notable for their α-farnesene (12–26%) and germacrene D (0.3–12%) content. <italic>A. arguta </italic>flowers contained branched-chain (13%) and aromatic alcohols (14–50%), and monoterpenes (12–24%), which were mostly linalool oxidation products reported previously [##REF##12737978##29##]. <italic>A. eriantha </italic>flowers contained few volatile compounds, the notable ones being 2-phenylethanol (8.5%), 6-methyl-5-hepten-2-one (also 8.5%), and β-myrcene (5%). The flowers from two further <italic>Actinidia </italic>species showed distinct volatile profiles. <italic>A. polygama </italic>flowers contained linalool oxides (36.5%) and dehydroiridodial (37%), a monoterpene related to the <italic>Nepeta cataria </italic>(catnip) compounds [##UREF##7##30##,##UREF##8##31##]. The flowers of <italic>A. indochinensis </italic>flowers were dominated by 2-phenylethanol (25.8%) and a couple of terpenes, the major one being linalool (23.6%). The aromas from these flowers have also had descriptors attached (McNeilage et al. unpublished observations). For example <italic>A. eriantha </italic>has been described as coconut and berry, and <italic>A. arguta </italic>with descriptors including sweet, lime, rose and vanilla [##REF##12737978##29##]. <italic>A chinensis </italic>and <italic>A. deliciosa </italic>are described as having tea-rose aromas, sweet, fresh, lilac and violet, while <italic>A. polygama </italic>has gardenia and lime descriptors. These descriptors reflect the diversity of volatile compounds found in the <italic>Actinidia </italic>flowers.</p> <p>In fruit, esters accounted for up to 85% of the total volatiles. The main <italic>A. chinensis </italic>esters were the fresh, fruity ethyl butanoate which, depending upon analytical method and genotype, ranged from 3 to 54%, and butyl butanoate, both of which are common fruit esters [##UREF##9##32##]. <italic>A. chinensis </italic>'Hort16A' fruit contained significant amounts of eucalyptol (0.2 to 21%) which has been described as having a fresh, pungent cooling taste [##UREF##9##32##]. Acetaldehyde, hexanal, E-2-hexenal and ethyl butanoate are known to be important contributors to Hort16A aromas [##REF##17616207##33##]. Ethyl butanoate along with methyl butanoate and ethyl acetate were the major (~60%) compounds in the green <italic>A. deliciosa </italic>fruit and are known to be major volatile contributors to its flavor [##UREF##10##34##]. The other major compound was ethanol (16%), with only trace levels of other alcohols and terpenes. Whilst at more modest levels in <italic>A. arguta </italic>fruit (2–8%), ethyl butanoate was still substantial. <italic>A. arguta </italic>fruit were notable for their content of methyl- and ethyl benzoate, which are responsible for the distinctive aroma of feijoa (<italic>Feijoa sellowiana </italic>O. Berg) fruit [##UREF##11##35##], and also camphor (7%). The dominant ester in <italic>A. eriantha </italic>fruit was butyl acetate (12%), a flavor compound found in many fruit, which has a strong, fruity odor and a taste reminiscent of pineapple [##UREF##9##32##,##UREF##12##36##]. However, the total volatile level in <italic>A. eriantha </italic>is lower than other species and the number of esters is also lower (Additional file ##SUPPL##6##7##).</p> <title>Ester biosynthesis</title> <p>Esters are synthesized from their acid and alcohol precursors by acyltransferases (ATs). We identified 30 ATs in the <italic>Actinidia EST </italic>database, of which 25 are full length. These sequences were analyzed for their phylogenetic relationship with other plant ATs of known function (Fig. ##FIG##2##3##). Two clades were identified that contained enzymes involved in the synthesis of flavor-related esters (Clades 1 and 2 in Fig. ##FIG##2##3##). The first contained alcohol acyltransferases (AATs) from apple, melon and banana and nine AATs from <italic>Actinidia </italic>(AdAT1, AdAT2, AeAT9, AcAT15, AcAT16, AdAT17, AaAT18, AcAT20, AcAT23), while the second clade contained AATs from strawberry and three AATs from <italic>Actinidia </italic>(AdAT6, AdAT22, AcAT24). An <italic>Actinidia </italic>AT predicted to be involved in esterifying anthocyanins (AcAT25) and two ATs predicted to be anthranilate acyltransferases (AdAT8, AcAT21) were also identified. Three <italic>Actinidia </italic>ATs were closely related to the carboxylesterase-related protein 2 CER2 and therefore predicted to be elongases (AcAT30, AcAT27 and AcAT26) possibly involved in wax production.</p> <p>Multiple sequences of the same genes were recovered for some genes reflecting orthologs, alleles or sequences from different genomes within polyploid genomes. For example, the three ATs, AcAT16, AdAT17 and AdAT1, and the two ATs, AaAT18 and AdAT2, are probable variants of just two different AATs. <italic>A. deliciosa </italic>is a hexaploid (6n) possibly explaining this multiplicity of genes. We may have uncovered at least one instance where, for a single AAT, all three genomic versions have been isolated (AcAT16, AdAT17 and AdAT1).</p> <p>Carboxylesterases conduct the opposite reaction to ATs, i.e. the hydrolysis of esters into acids and alcohols [##UREF##13##37##]. As well as ester hydrolysis, plant carboxylesterases have been implicated in isoflavonoid biosynthesis [##REF##15734910##38##], plant defense [##REF##8012404##39##], and hormone regulation [##REF##17644730##40##]. Analysis of <italic>Actinidia </italic>flowers and fruits showed a wide range of acids and alcohols (Additional file ##SUPPL##6##7##).</p> <p>Previously published plant carboxylesterase (CXE) genes from GenBank were used to identify CXE genes in the <italic>Actinidia </italic>EST database using BLAST searches. The sequence alignment of CXEs revealed amino acid motifs and secondary structural features characteristic for members of the α/β hydrolase fold superfamily and more particularly of CXEs such as an active site serine surrounded by a GXSXG motif.</p> <p>In <italic>Actinidia </italic>19 CXEs were identified, 18 of which are full length (Fig. ##FIG##3##4##). These include CXE members likely to encode gibberellin receptors (AdCXE9, AdCXE17, AdCXE15), an isoflavanone dehydratase (AeCXE11) and an ortholog of a plant defense-associated CXE (AcCXE12). Other CXEs from <italic>Actinidia </italic>align well with <italic>Arabidopsis </italic>CXEs and can be associated with the classes identified elsewhere [##REF##14738307##41##]. The structure of an <italic>Actinidia </italic>CXE is now available (AeCXE1), the first plant CXE to have its structure solved [##REF##17256879##42##]. AeCXE1 can hydrolyze a range of substrates from C2 to C16 esters, with a preference for C4 moieties. There are two CXEs from <italic>Actinidia </italic>(AdCXE16, AaCXE2) that may be orthologs to AeCXE1 from <italic>A. eriantha</italic>.</p> <title>Terpene biosynthesis</title> <p><italic>Actinidia </italic>terpenoid compounds such as monoterpenes (C10), sesquiterpenes (C15) diterpenes (C20), triterpenes (C30), carotenoids, sterols, phytols and quinones are derived from two common precursor molecules, isopentenyl diphosphate and dimethylallyl diphosphate. These precursors are produced via the cytoplasmic mevalonate pathway or the chloroplastic 1-deoxy-D-xylulose-5-phosphate (mevalonate-independent) pathway [##REF##11078528##43##]. Genes for enzymes in both these pathways are present in the <italic>Actinidia </italic>EST database, with the results for analysis of ESTs in the early steps of the mevalonate pathway shown in Fig. ##FIG##4##5##. The conversion of these precursors into other intermediates of the terpenoid biosynthetic pathway, including geranyl diphosphate, farnesyl diphosphate, and geranylgeranyl diphosphate, is carried out by polyisoprene synthase genes (EC 2.5.1.x). Over 200 ESTs with sequence homology to known polyisoprene synthase genes were found in the <italic>Actinidia </italic>EST database including 116 ESTs with similarity to dimethyallyltranstransferase (EC 2.5.1.1), 66 ESTs with homology to known monoterpene and sesquiterpene synthases and 38 ESTs homologous to squalene synthase (EC 2.2.1.21). No ESTs for phytoene synthase (EC 2.5.1.32) were identified.</p> <p>The EST database collection has been used to identify a multifunctional terpene synthase gene from <italic>A. deliciosa </italic>flowers that produces the sesquiterpene germacrene D and a range of other sesquiterpene products at lower abundance [##UREF##14##44##]. This gene is represented by 10 sequences in the EST database.</p> <title>Genes related to color</title> <title>Chlorophyll catabolism</title> <p>Green fleshed commercial kiwifruit such as <italic>A. deliciosa </italic>'Hayward' are characterized by the retention of chlorophylls in the flesh of ripe fruit [##REF##11754554##45##]. In contrast, the newer yellow-fleshed fruit cultivars (e.g. <italic>A. chinensis </italic>'Hort16A' and 'Jintao') degrade the chlorophyll during ripening and the consequent loss of the green color uncovers the underlying yellow carotenoid pigmentation.</p> <p>The degradation of chlorophylls (leading to the production of colorless breakdown products) is controlled by catabolic enzymes in the chlorophyll degradation pathway [##REF##16669755##46##]. <italic>Actinidia </italic>ESTs discovered in the pathway of chlorophyll breakdown are shown in Fig. ##FIG##5##6##. Chlorophyllase is the first enzyme of the catabolic pathway and it removes the phytol chain from chlorophyll <italic>a</italic>. Homologs of this enzyme have been mainly identified in libraries of vegetative tissues (breaking buds, leaves and petals) of <italic>A. deliciosa </italic>and <italic>A. chinensis</italic>. Homologs of pheophorbide <italic>a </italic>oxygenase have also been identified, with eight ESTs representing two genes from fruit libraries of <italic>A. chinensis</italic>. This enzyme is responsible for opening the tetrapyrrolic ring of the chlorophyll molecules, the critical step for the loss of green color. Surprisingly, we did not identify any ESTs of red chlorophyll catabolite reductase, an enzyme associated with pheophorbide <italic>a </italic>oxygenase and required to complete the porphyrin ring opening and the degradation of chlorophyll.</p> <title>Carotenoid biosynthesis</title> <p>The fruit of many <italic>Actinidia </italic>species contain a range of carotenoid pigments, including β-carotene and lutein [##REF##11754554##45##]. The yellow-fleshed kiwifruit <italic>A. chinensis </italic>'Hort16A' develops as a green-fleshed fruit, but on ripening the chlorophyll is degraded exposing the carotenoids in the yellow flesh. ESTs for most steps in the biosynthesis of carotenoids [##REF##15012246##47##, ####REF##16669779##48##, ##REF##16669781##49####16669781##49##] are present in the <italic>Actinidia </italic>EST database (Fig. ##FIG##6##7##), except the first step, phytoene synthase. Ninety eight ESTs for phytoene desaturase were identified, with 94 coming from <italic>A. deliciosa </italic>libraries, and the remainder from <italic>A. chinensis </italic>libraries. Of these 98 ESTs, 90 came from a single <italic>A. deliciosa </italic>fruit library with eight from bud and meristem libraries. Of the 56 ESTs found for lycopene β-cyclase, 23 were from petal tissues and 32 from fruit libraries. These ESTs were largely found in libraries from <italic>A. deliciosa </italic>(31) and <italic>A. chinensis </italic>(24), with only one EST from <italic>A. eriantha</italic>. In contrast, only a single EST was found for lycopene ε-cyclase, which is the enzyme required with lycopene β-cyclase to convert lycopene to α-carotene [##REF##8837512##50##]. There were five ESTs of β-carotene hydroxylase found in petal (3), leaf (1) and dormant bud (1) libraries, but only a single EST was present for ε-carotene hydroxylase. Both of these enzymes are required for lutein synthesis [##REF##16492736##51##].</p> <title>Flavonoid and anthocyanin biosynthesis</title> <p>Flavonoids are polyphenolic secondary metabolites synthesized through the phenylpropanoid pathway [##UREF##15##52##]. Anthocyanins are a subgroup of the flavonoids, which appear red to blue, depending on the pH, due to additions to their phenolic rings. In red fruit of <italic>A. chinensis </italic>and <italic>A. deliciosa</italic>, the predominant anthocyanins are cyanidin-based and the preferred glycosylation is in the 3' position [##REF##16302772##53##]. While many of the kiwifruit libraries came from tissue with low levels of anthocyanins, several libraries were derived from tissue containing anthocyanins, including red fruits, brown skins and pink petals, leaves and buds (Additional file ##SUPPL##1##2a##).</p> <p>The anthocyanin biosynthetic pathway was analyzed by BLAST searching for ESTs representing different enzymes in the pathway (Fig. ##FIG##7##8##). Chalcone synthase (CHS) was a highly represented gene family in the EST database with ~300 ESTs. There appeared to be three distinct CHS NRs, possibly representing three genes, that are found in mostly young <italic>A. eriantha </italic>fruit libraries, as well as in bud and petal libraries. There was only one apparent chalcone isomerase (CHI) NR and over half the ESTs for this gene were found in shoot and bud libraries. Two distinct flavanol 3-hydroxylase (F3H) NRs were found, with ESTs predominantly found in bud and cell culture libraries. Only one flavonoid 3',5'-hydroxylase (F3'5'H) NR was found, from <italic>A. eriantha </italic>fruit, while ESTs for two distinct flavonoid 3'-hydroxylase (F3'H) NRs were predominantly found in shoot and bud libraries. Of the four potential dihydroflavonol reductase (DFR) NRs, the most highly represented was found in bud and young fruit. Three NRs, perhaps representing three gene families, of leucoanthocyanidin dioxygenase (LDOX) were found in the EST database, with most ESTs in bud libraries. The ESTs with highest sequence similarity to UDP-glucose flavonoid 3-O-glucosyltransferase (UFGT) are also predominantly found in developing buds of <italic>A. deliciosa</italic>.</p> <p>Many of the ESTs for anthocyanin biosynthesis were sequenced from young fruit libraries (particularly <italic>A. eriantha</italic>) and in <italic>A. deliciosa </italic>buds. Developing buds from <italic>A</italic>. <italic>deliciosa </italic>appear to be a darker red color and therefore presumably accumulate more anthocyanins than buds from <italic>A. chinensis </italic>and <italic>A. eriantha</italic>, so it is not surprising to find ESTs for these later biosynthetic steps in bud tissue from <italic>A. deliciosa</italic>. However, ESTs for genes encoding enzymes from the early steps in the flavonoid pathway (e.g. CHS and F3H) are found in unpigmented <italic>A. eriantha </italic>fruit suggesting these tissues may accumulate a different class of flavonoids.</p> <p>Most of the enzymes for anthocyanin biosynthesis are well represented in the EST database except for F3'5'H. This is not surprising given the only anthocyanins identified in <italic>A. deliciosa </italic>and <italic>A. chinensis </italic>are glycosylated in the 3' position. There are, however, species within the genus <italic>Actinidia</italic>, such as <italic>A. melanandra</italic>, that accumulate high levels of delphinidin, a different aglycone that would require F3'5'H activity. No ESTs from these species were sequenced.</p> <title>Genes related to healthful components</title> <title>Ascorbic acid-related genes</title> <p>Kiwifruit are well known to be high in ascorbic acid (vitamin C), with values in commercial varieties ranging from 80 to 120 mg per 100 g fresh weight of fruit [##UREF##16##54##], and with <italic>A. eriantha </italic>having up to 10 times as much again. Ascorbic acid is synthesized in plants from glucose, through three potential pathways: the L-galactose pathway, the galacturonate pathway and through <italic>myo</italic>-inositol [##UREF##17##55##]. <italic>Myo</italic>-inositol is a major sugar alcohol found in kiwifruit (Additional file ##SUPPL##6##7##; [##UREF##18##56##,##UREF##19##57##]). The only pathway that is well established is the L-galactose pathway for which all the enzymes have now been identified in plants [##UREF##17##55##,##REF##17485667##58##]. Many of the genes are known for the other two potential pathways, but several steps are still unknown. The <italic>Actinidia </italic>EST database was searched for ESTs similar to genes known to encode steps in these pathways. The numbers of ESTs homologous to various identified genes in the three pathways of ascorbic acid biosynthesis are shown in Fig. ##FIG##8##9##. Surprisingly for these high vitamin C plants, while the number of ESTs found in early steps in all the pathways is relatively large, the number of ESTs found in the two last steps is much lower, with no examples of ESTs encoding the last step of two of the pathways, that catalyzed by L-galactono lactone dehydrogenase.</p> <title>Quinic acid metabolism</title> <p>Fruit from <italic>Actinidia </italic>species have a relatively high total acid content (1–3% w/w) of which 40–60% can be quinic acid, 40–60% citric acid and 10% malic acid (Additional file ##SUPPL##6##7##). Quinic acid comprises an even higher proportion of the total acids during early fruit development [##UREF##20##59##], and over 25-fold variation in quinic acid content can be found in fruit of a single <italic>A. chinensis </italic>cross [##UREF##21##60##]. Kiwifruit are unusual in having such a high content of quinic acid as a fresh fruit, although levels of quinic acid in cranberry juice can be around 1% [##REF##12405790##61##] and peaches can have up to 0.2% quinic acid [##UREF##22##62##].</p> <p>Quinic acid is a key intermediate in lignin biosynthesis, folic acid metabolism, aromatic acid synthesis, anthranilate biosynthesis, and purine metabolism. The molecular and enzymatic control of quinic acid storage and metabolism may affect all these pathways. However, little is known about quinic acid metabolism in <italic>Actinidia </italic>species. ESTs for all the enzymes in the quinic acid biosynthetic pathway to shikimate (Fig. ##FIG##9##10##), except for a quinate dehydrogenase, have been identified in the <italic>Actinidia </italic>EST database.</p> <p>Of the 13 ESTs identified for dehydroquinate synthase, eight were found in fruit or petal libraries, the rest were sequenced from bud or leaf libraries. Of the 16 ESTs identified for the bifunctional dehydroquinate dehydratase/shikimate dehydrogenase, 11 were from fruit or petal libraries and five from bud or leaf libraries. These ESTs were found mainly in libraries made from <italic>A. deliciosa </italic>and <italic>A. chinensis </italic>tissues and to a lesser extent in <italic>A. arguta </italic>and <italic>A. eriantha </italic>libraries, reflecting the distribution of ESTs among the species. An enzyme has been reported with predominantly quinate dehydrogenase activity [##UREF##23##63##], but the gene for this enzyme has not been cloned.</p> <title>Allergens</title> <p>Food allergies are estimated to affect ~6% of young children and ~3% of adults [##REF##15842229##64##]. Kiwifruit have been recognized as a potentially allergenic fruit for over 20 years and reported allergies to kiwifruit are increasing [##REF##14675467##65##]. Most allergy symptoms to kiwifruit are quite mild, but severe reactions have been reported, particularly in young children [##REF##15248859##66##]. Allergies have been recorded in all three commercial species of <italic>Actinidia </italic>[##REF##16481086##67##].</p> <p>The <italic>Actinidia </italic>EST database contains sequences with homology to many known plant allergen proteins [##UREF##24##68##]. Potential allergens include the 2S albumin proteins, lipid transfer proteins, thaumatin-like protein (TLP), α-amylase/trypsin inhibitors, latex allergens, plant chitinases, profilins, cystatins, Bet v 1 homologous proteins and plant seed globulin allergens. However, only a small number of potential kiwifruit allergens have been directly confirmed by immunological testing, including the cysteine protease actinidin (Act d 1; [##REF##9564807##69##]), an unidentified 43 kDa protein (Act d 2; [##UREF##25##70##], a TLP [##UREF##25##70##] and kiwellin [##REF##16328735##25##].</p> <p>ESTs for actinidin occurred throughout the 37 libraries from seven species that make up the database, and these clustered into acidic and basic forms of 10 distinct genes/alleles [##UREF##5##22##]. No matches were found by BLAST searching the EST database for the 43 kDa protein Act d 2. Kiwellin ESTs were found in libraries from most <italic>Actinidia </italic>species but were particularly abundant in <italic>A. eriantha </italic>ripe fruit skin, accounting for ~40% of ESTs in that library. The EST database contained 24 NRs with homology to TLPs from other plants (e.g. Pru a 2 from cherry and thaumatin 1 from <italic>Thaumatococcus daniellii</italic>), and all but one EST contained the 16 conserved cysteine residues characteristic of this protein class [##REF##15080826##71##]. The various NRs could be aligned and clustered into 4 or 5 groups (data not shown), with two NRs accounting for the bulk of the ESTs (46 and 29%).</p> <p>Globulin proteins (11S and 7S) have been reported from a wide range of seeds, including cereals [##REF##3985947##72##,##UREF##26##73##] and legumes [##REF##7786405##74##]. Globulin-like proteins are also found in <italic>A. thaliana </italic>[##REF##12417707##75##] and a gene family of 10 globulin-like genes are present in <italic>Arabidopsis </italic>gene databases (The <italic>Arabidopsis </italic>Information Resource (TAIR), <ext-link ext-link-type="uri" xlink:href="http://arabidopsis.org/"/>). 11S globulins are composed of sets of polymorphic subunits derived from a multigene family. They exist as hexamers with molecular mass between 300 and 400 kDa, or as trimers. 11S globulin subunits consist of two polypeptide chains linked by at least one S-S bridge between cysteine residues at highly conserved positions in the acidic alpha-chain and basic beta-chain. Both chains are post-translationally generated from a common precursor protein, which represents the product of one member of the multigene family [##UREF##27##76##]. 11S globulins have been shown to be related to the 7S globulins through sequence comparisons [##UREF##28##77##].</p> <p>A BLAST search of the <italic>Actinidia </italic>EST sequence database with a range of previously identified globulin protein sequences identified two ESTs with a predicted protein sequence matching 11S globulin, both from a library made from <italic>A. chinensis </italic>whole fruit, which included seeds. Other related sequences from <italic>Actinidia </italic>fruit libraries belonged to the 7S vicilins-like family (38 ESTs), which were found in only fruit-based libraries from <italic>A. deliciosa</italic>, <italic>A. chinensis </italic>and <italic>A. arguta</italic>. These libraries would have all included seeds. A third grouping of sequences was identified that we have named other globulins (OG) (201 ESTs). The OG group appears to be more closely related to the 7S globulins than to the 11S (Fig ##FIG##10##11##). All three groupings are supported by the clustering of the <italic>Arabidopsis </italic>globulin genes as well as globulin genes from other species. While most OG ESTs were found in fruit libraries (81%), which included seeds, some members of the OG family were found in fruit skin, buds and petal libraries, suggesting a role other than as seed storage proteins. These OG ESTs were found in libraries derived from <italic>A. deliciosa </italic>(100), <italic>A. chinensis </italic>(42), <italic>A. arguta </italic>(39), and <italic>A. eriantha </italic>(16).</p> <title>Genes related to fruit softening</title> <p>Some kiwifruit cultivars exhibit outstanding storage characteristics e.g. <italic>A. deliciosa </italic>'Hayward' fruit can be stored for 4–6 months at 0°C. Other <italic>Actinidia </italic>species exhibit a range of ripening and softening behaviors [##UREF##29##78##]. <italic>A. eriantha </italic>develops a peelable skin as the fruit ripen, <italic>A. arguta </italic>ripens in less than 10 days at 20°C versus 20–25 days for <italic>A. deliciosa </italic>and <italic>A. chinensis </italic>[##UREF##30##79##], whilst some small-fruited <italic>Actinidia </italic>genotypes tend to remain firm even towards the end of the ripening process [##UREF##29##78##].</p> <p>The main chemical changes occurring in the cell wall during kiwifruit softening are pectin solubilization and degradation, reduction in the molecular weight of xyloglucan, and galactose loss from pectin side chains. Microscopically, the cell wall shows extensive swelling, until at the end of the softening process dissolution of middle lamellae can be observed [##UREF##31##80##]. These processes eventually lead to disintegration of the cell wall. Although similar chemical changes occur in other fruit [##REF##11554479##81##], these changes occur concurrently with ethylene production and the respiratory climacteric. In kiwifruit, however, most of the softening process is well separated temporally from the climacteric and ethylene production (shown schematically in Fig. ##FIG##11##12##).</p> <p>The genes for many enzymes involved in the key chemical changes in the cell wall during ripening are represented in the <italic>Actinidia </italic>EST database. Three enzymes involved in pectin degradation (pectin lyase, pectin methylesterase and polygalacturonase) are particularly abundant with 244, 611 and 414 ESTs, respectively. Genes encoding expansin, a protein implicated in pectin solubilization early in the kiwifruit softening process [##UREF##31##80##], are also highly abundant (468 ESTs). Galactose loss from fruit cell walls during the softening process has been attributed to the action of β-galactosidases [##REF##9625694##82##]. A β-galactosidase from ripe kiwifruit has been purified. However, the activity of the enzyme <italic>in vitro </italic>against synthetic or purified native substrates from kiwifruit was far too low to account for the amount of galactose loss observed during softening [##UREF##32##83##]. Ninety β-galactosidase ESTs were observed in fruit libraries in the EST database and are good candidates for genetic manipulation and to study the role of this enzyme. During softening, xyloglucan is also hydrolyzed. As xyloglucan and cellulose create the major load-bearing network of the cell wall, any xyloglucan degradation is thought to weaken the cell wall, resulting in softening of the fruit. In kiwifruit, hydrolysis of xyloglucan occurs mainly during the rapid softening phase. Endoglucanases and xyloglucan transglucosylase/hydrolase (XTH) enzymes are implicated in this process. XTHs catalyze both hydrolytic and transglucosylation reactions, and an XTH enzyme capable of carrying out both reactions has been isolated from kiwifruit [##REF##9487728##84##]. This particular XTH, however, is expressed late in kiwifruit ripening, when xyloglucan degradation has already come to an end. The 96 XTH ESTs found in fruit libraries will allow the identification of genes involved in hydrolysis of xyloglucan during the rapid softening phase. ESTs for enzymes that act on minor kiwifruit cell wall components such as galactoglucomannans or glucuronoarabinoxylans [##UREF##31##80##] were less abundant in the EST database (e.g. mannan transglycosylase/hydrolase (MTH), 10 ESTs found in fruit).</p>
<title>Discussion</title> <p>We report on a significant resource of over 130,000 ESTs derived from a range of <italic>Actinidia </italic>species (Table ##TAB##0##1##). We targeted tissues and developmental stages in order to sample genes involved in physiological and biochemical processes including fruit ripening, flavor development, control of color and the synthesis of chemicals with health-related attributes. For this reason, the two most widely cultivated species of kiwifruit, <italic>A. chinensis </italic>and <italic>A. deliciosa</italic>, are well represented with together over 100,000 ESTs (Table ##TAB##0##1##). In addition, fruit and bud libraries are also well represented, with over 38,000 and 50,000 ESTs, respectively. <italic>A. chinensis </italic>and <italic>A. deliciosa </italic>are so closely related, as is <italic>A. setosa </italic>(Li) C.F. Liang et A.R.Ferguson, that they are variously treated as being distinct species or as varieties of the one species. The other two main species studied, <italic>A. arguta </italic>and <italic>A. eriantha </italic>Benth., also have commercial potential but are more distantly related [##UREF##1##2##].</p> <p>The genus <italic>Actinidia </italic>is unusual in how much inter-taxal and intra-taxal variation in ploidy and in the wild, there is a structured reticulate pattern of diploids, tetraploids, hexaploids, and octoploids in diminishing frequency, associated, in at least some taxa, with geographic separation of ploidy races. <italic>A. deliciosa </italic>is hexaploid, <italic>A. setosa </italic>is diploid, and there are diploid and tetraploid races of <italic>A. chinensis</italic>, the tetraploids apparently coming from a restricted part of the natural distribution of the species. Most evidence suggests that diploid <italic>A. chinensis </italic>was a progenitor of tetraploid <italic>A. chinensis </italic>and hexaploid <italic>A. deliciosa </italic>but it is not clear whether genomes from other species have contributed. The basic chromosome number (n = 29) is high and it seems increasingly likely that diploid <italic>A. chinensis </italic>is itself a rediploidized palaeopolyploid [##UREF##1##2##].</p> <p>As is common in EST sequencing projects (e.g., [##REF##15361140##6##,##REF##17761391##8##,##REF##16531485##10##], there is a high degree of redundancy in the ESTs, with clustering reducing the number of unique sequences from over 132,000 to 41,858 NRs (18,070 TCs, 23,788 singletons). We would expect this number of NRs to be an overestimate of the number of genes in <italic>Actinidia</italic>, especially given that the database contains sequences from multiple species of <italic>Actinidia</italic>. Using the same correction used in the apple EST paper [##REF##16531485##10##], we expect an <italic>Actinidia </italic>genome to have around 27,000 genes.</p> <p>On average 20% (± 2% standard error) of the sequences from each library with over 1000 ESTs were singletons suggesting a high degree of novelty in these libraries. On average 28% (± 4%) of sequences did not have a homolog in the various public databases based on BLAST searches with an E value &gt; 1.0e<sup>-10</sup>. An average of 16% (± 3%) of ESTs were identified as 3' UTR candidates based on the presence of a poly(A) tail within 40 bp of the start (taking into account reverse sequences). These 3' sequences would not be expected to be identified by BLAST searches and so would affect the novelty of a library. Less than 12% of NRs did not have BLAST matches (E &gt; 10) in the <italic>Arabidopsis </italic>proteome, Uniref, NCBI ref or SwissProt databases.</p> <p>There was only a small degree of overlap in NRs between libraries. Libraries from different species and different tissues showed a 5 to 9% overlap in NRs, libraries from different species but the same tissue showed an 8 to 10% overlap and libraries from the same species but different tissue showed a 7 to 13% commonality in NRs. These comparisons were made over five large libraries with more than 9,000 EST members each and an average of 2.1 ESTs per NR. These results suggest that there were more NRs in common between libraries made from the same tissue or from the same species, but this tendency was not particularly marked.</p> <p>Detecting SNPs using an automatically assembled EST database is a cost effective way to discover new DNA polymorphisms and develop novel markers, although it can be a challenging task, especially in polyploid <italic>Actinidia </italic>species. A significant proportion of the sequence variants predicted from overlapping ESTs within an NR will correspond to "real" SNPs, which means the sequence differences found are allelic variants of a given locus and not sequencing errors or differences between paralogs, homoeologs or orthologs. Homoeolog SNPs could be particularly common in the polyploid accessions of species such as <italic>A. deliciosa </italic>and <italic>A. arguta </italic>that make up a large proportion of this database, but are also possible in diploids as a result of conserved gene pairs of paleopolyploid origin. Allelic SNPs can be used directly and converted into molecular markers for genetic mapping, population genetics and linkage disequilibrium studies or for marker-assisted selection. A SNP marker for determining the sex of kiwifruit seedlings [##UREF##33##85##] has already been successfully utilized. Since the database contains sequence data from multiple species, and ~40% of TCs are made up of more than one species, several SNPs were detected in the <italic>Actinidia </italic>EST database corresponding to sequence between orthologous loci from different <italic>Actinidia </italic>species. Hence, they cannot fully be considered as allelic SNPs, but more as species-specific variations. However, as kiwifruit breeding programs often use controlled crosses between different species, the interspecific SNPs will segregate in the progeny and be useful as markers.</p> <p>The incidence of SSRs in NRs was higher in <italic>Actinidia </italic>(30%) than in apple (20%), and the frequency of di-nucleotide and tri-nucleotide SSRs differed between these two species. This increase was evident in all of the sequence classes but greatest in AG and AC (double the incidence among apple NRs). Even though the <italic>Actinidia </italic>genome EST resource represents several species, while apple came mainly from one species, this would not explain these differences. Perhaps the longer period of domesticity in apple, based on a narrow genetic basis compared to kiwifruit, may explain the difference. Alternatively it may reflect that a greater proportion of homoeologs have grouped into TCs in the polyploid kiwifruit data than in the apple dataset.</p> <p>Overall the codon usage of the three <italic>Actinidia </italic>species shares many similarities with that of other dicotyledons represented in the codon usage database [##REF##10592250##86##]. Comparisons with <italic>Arabidopsis </italic>codon usage showed that <italic>A. deliciosa </italic>and <italic>A. eriantha </italic>differ markedly for 15 and 17 amino acids, respectively, whereas <italic>A. chinensis </italic>differs in its preference for eight particular amino acids. Further comparisons with apple, grape, pear, peach, loblolly pine, tomato, citrus, potato and tobacco showed that the codon preference of the <italic>Actinidia </italic>species is most similar to that of apple [##REF##16531485##10##]. <italic>A. deliciosa </italic>differs from apple only in its codon preference for aspartate, glycine, isoleucine and leucine. <italic>A. eriantha </italic>also differs from apple for these four amino acids and also serine. The codon preference for <italic>A. chinensis </italic>and apple also differ for only four amino acids, these being asparagine, glutamine, threonine and valine. CpG suppression is also evident in <italic>Actinidia </italic>species with an XCG/XCC ratio of between 0.68 and 0.71 for the three species evaluated. This modest level of suppression of the CpG di-nucleotides is similar to that of apple (0.64) and differs markedly from that of <italic>Arabidopsis </italic>which shows nearly no suppression (0.92) and from the high level found in grape (0.35). This may well reflect different levels of methylation in the coding sequences used by different species of plants.</p> <p>Mapman was used to assign function to the <italic>Actinidia </italic>NRs and thus to their constitutive ESTs. Only 32% of the ESTs did not have an <italic>Arabidopsis </italic>homolog at E &lt; 1.0e<sup>-10</sup>. In general, the functional distribution of NRs was very similar to the functional distribution of <italic>Arabidopsis </italic>proteins (Table ##TAB##2##3##) suggesting that the sampling of <italic>Actinidia </italic>ESTs well represented the major functional classes of plant genes. This is surprising given the biased selection of libraries with virtually no root ESTs sequenced. However, the high number of bud meristem libraries meant that genes expressed in metabolically active dividing tissue were sampled.</p> <p>Fruit of the <italic>Actinidia </italic>genus show several characteristics that distinguish them from other fruit species. These include flesh color (green is the most common, but yellow, orange and red fruit also occur in the genus), chemical composition including high vitamin C and quinic acid contents, and a novel aroma composition (Additional file ##SUPPL##6##7##), characterized by the abundant esters. In addition, kiwifruit has been identified as a fruit with a potential to cause allergenicity among consumers, although this is a problem common to many other fruit. For this reason, we analyzed the <italic>Actinidia </italic>EST database to identify genes involved in these pathways and products. These analyses showed the depth and usefulness of the database for selecting candidate genes for most steps in the selected pathways. The other useful characteristic of the <italic>Actinidia </italic>EST database is the wide range of genetic and phenotypic diversity sampled across the <italic>Actinidia </italic>genus (Fig. ##FIG##0##1##) and the value of using this diversity to discover novel traits through functional genomics and through mapping and positional cloning approaches.</p>
<title>Conclusion</title> <p>This paper describes an EST resource in the <italic>Actinidia </italic>genus and discusses many of the properties of this collection. However, there is still a tremendous challenge in understanding the molecular basis of the genetic diversity of this genus, and we expect putting this EST resource into the public domain will enhance future understanding of the genetic basis of the many divergent traits in this fruit.</p>
<p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <title>Background</title> <p>Kiwifruit (Actinidia spp.) are a relatively new, but economically important crop grown in many different parts of the world. Commercial success is driven by the development of new cultivars with novel consumer traits including flavor, appearance, healthful components and convenience. To increase our understanding of the genetic diversity and gene-based control of these key traits in Actinidia, we have produced a collection of 132,577 expressed sequence tags (ESTs).</p> <title>Results</title> <p>The ESTs were derived mainly from four <italic>Actinidia </italic>species (<italic>A. chinensis, A. deliciosa, A. arguta </italic>and <italic>A. eriantha</italic>) and fell into 41,858 non redundant clusters (18,070 tentative consensus sequences and 23,788 EST singletons). Analysis of flavor and fragrance-related gene families (acyltransferases and carboxylesterases) and pathways (terpenoid biosynthesis) is presented in comparison with a chemical analysis of the compounds present in <italic>Actinidia </italic>including esters, acids, alcohols and terpenes. ESTs are identified for most genes in color pathways controlling chlorophyll degradation and carotenoid biosynthesis. In the health area, data are presented on the ESTs involved in ascorbic acid and quinic acid biosynthesis showing not only that genes for many of the steps in these pathways are represented in the database, but that genes encoding some critical steps are absent. In the convenience area, genes related to different stages of fruit softening are identified.</p> <title>Conclusion</title> <p>This large EST resource will allow researchers to undertake the tremendous challenge of understanding the molecular basis of genetic diversity in the <italic>Actinidia </italic>genus as well as provide an EST resource for comparative fruit genomics. The various bioinformatics analyses we have undertaken demonstrates the extent of coverage of ESTs for genes encoding different biochemical pathways in <italic>Actinidia</italic>.</p>
<title>Authors' contributions</title> <p>EFW, EHAR, EJFS, GSR, JHB, KCS, MDT, MW, RMD, and RS all contributed to bioinformatics analysis. ACA, CAD, KBM, LLB, MM, MR, RDN, RGA, RPH, and SMB all undertook bioinformatics analysis and manuscript drafting. RNC undertook bioinformatics programming, HLB, and KK, did chemical analysis. AJM did chemical analysis and interpretation and manuscript writing. APG and EAM undertook design and manuscript drafting. BJJ, BN, DW, EG, KRL, MYW, RE, RW, YKY, and YYW did the EST sequencing and clone selection. ACR, ARF, LGF and MAN provided basic genetic material and analysis and manuscript drafting. WAL did manuscript drafting and bioinformatics analysis. BU contributed significantly by doing Mapman analysis and manuscript drafting. MWD and RJS did microarray analysis and interpretation. DC did SNP analysis and manuscript writing.</p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>We thank Ian Ferguson and Richard Forster for support in this project. Libraries were sequenced by Genesis Research &amp; Development Corporation Limited, Auckland, New Zealand. This analysis project was supported financially by the Foundation for Science, Research and Technology (contract C06X0403) and by HortResearch.</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p>Fruit diversity in the genus <italic>Actinidia</italic>. Fruit of species used to make EST libraries are identified by letters. A is <italic>A. arguta</italic>, C is <italic>A. chinensis</italic>, D is <italic>A. deliciosa</italic>, E is <italic>A. eriantha</italic>, I is <italic>A. indochinensis</italic>, P is <italic>A. polygama </italic>and S is <italic>A. setosa</italic>. A. <italic>hemsleyana </italic>is not in the photos.</p></caption></fig> <fig position="float" id="F2"><label>Figure 2</label><caption><p>Plot of <italic>Actinidia deliciosa </italic>versus <italic>A. eriantha </italic>microarray intensities. Genomic DNA was labeled with cy3 and cy5 in two separate reactions, combined and hybridized to a single microarray slide. Spots not found or flagged "bad" were removed and each channel was normalized using a quantile normalization. Of the two arrays, 15,024 were flagged "good" using <italic>A. deliciosa </italic>DNA and 13,785 were flagged "good" using <italic>A. eriantha </italic>DNA. Of these spots, 13,443 were common to both species. Black spots represent oligonucleotides that showed similar intensity of hybridization between the two species. Triangles represent oligonucleotides that showed an increase in intensity with <italic>A. eriantha </italic>genomic DNA and crosses represent oligonucleotides that showed an increase using <italic>A. deliciosa </italic>genomic DNA.</p></caption></fig> <fig position="float" id="F3"><label>Figure 3</label><caption><p>Phylogenetic tree of plant acyltransferases (ATs) of known function and <italic>Actinidia </italic>members of the acyltransferase family. Previously published plant acyltransferase sequences from GenBank were used to identify genes from the <italic>Actinidia </italic>EST database using BLAST searches. Abbreviations for species and AT names are as follows: AdAT, <italic>Actinidia deliciosa </italic>AT; AeAT, <italic>A. eriantha </italic>AT; AcAT, <italic>A. chinensis </italic>AT; AaAT, <italic>A arguta </italic>AT; DAT, <italic>Catharanthus roseus </italic>deacetylvindoline 4-0-acetyltransferase, GenBank Accession No. <ext-link ext-link-type="gen" xlink:href="AF053307">AF053307</ext-link>; MAT, <italic>C. roseus </italic>minovincinine 19-hydroxy-O-acetyltransferase, <ext-link ext-link-type="gen" xlink:href="AAO13736">AAO13736</ext-link>; BEAT, <italic>Clarkia breweri </italic>acetyl-CoA:benzylalcohol acetyltransferase, <ext-link ext-link-type="gen" xlink:href="AF043464">AF043464</ext-link>; SALAT, <italic>Papaver somniferum </italic>salutaridinol 7-O-acetyltransferase, <ext-link ext-link-type="gen" xlink:href="AF339913">AF339913</ext-link>; BEBT, <italic>Clarkia breweri </italic>benzoyl-CoA:benzyl alcohol benzoyl transferase, <ext-link ext-link-type="gen" xlink:href="AF500200">AF500200</ext-link>; MpAAT1, <italic>Malus pumila </italic>alcohol acyltransferase; <ext-link ext-link-type="gen" xlink:href="AY707098">AY707098</ext-link>; CM-AAT1, <italic>Cucumis melo </italic>alcohol acyltransferase, <ext-link ext-link-type="gen" xlink:href="CAA94432">CAA94432</ext-link>; SAAT, <italic>Fragaria × ananassa </italic>alcohol acyltransferase, <ext-link ext-link-type="gen" xlink:href="AAG13130">AAG13130</ext-link>; HCBT, <italic>Dianthus caryophyllus </italic>anthranilate N-hydroxycinnamoyl benzoyltransferase, <ext-link ext-link-type="gen" xlink:href="Z84383">Z84383</ext-link>; AnthocyaninAT <italic>Petunia frutescens </italic>anthocyanin acyltransferase, <ext-link ext-link-type="gen" xlink:href="BAA93453">BAA93453</ext-link>; VAAT, <italic>Fragaria vesca </italic>alcohol acyltransferase, <ext-link ext-link-type="gen" xlink:href="AX025504">AX025504</ext-link>; BanAAT, <italic>Musa acuminata </italic>alcohol acyltransferase, <ext-link ext-link-type="gen" xlink:href="AX025506">AX025506</ext-link>; CER2, <italic>A. thaliana </italic>CER2 gene, <ext-link ext-link-type="gen" xlink:href="X93080">X93080</ext-link>; HCT, <italic>Nicotiana tabacum </italic>hydroxyl-cinnamoyl transferase, <ext-link ext-link-type="gen" xlink:href="AJ507825">AJ507825</ext-link>; HQT, <italic>Lycopersicon esculentum </italic>hydroxycinnamoyl CoA quinate transferase, <ext-link ext-link-type="gen" xlink:href="AJ582652">AJ582652</ext-link>. Percentage bootstrap values (1000 bootstrap replicates) for groupings are given by each branch. Transf. on cladogram is transferase.</p></caption></fig> <fig position="float" id="F4"><label>Figure 4</label><caption><p>Phylogenetic tree of plant carboxylesterases of known function and <italic>Actinidia </italic>members of the carboxylesterase family (CXE). Previously published plant carboxylesterase sequences from GenBank were used to identify genes in the <italic>Actinidia </italic>EST database using BLAST searches. Abbreviations for species and CXE names are as follows: AtCXE, <italic>Arabidopsis thaliana</italic>, CXE; AaCXE, <italic>A. arguta </italic>CXE; AcCXE, <italic>A. chinensis </italic>CXE; AdCXE, <italic>A. deliciosa </italic>CXE; AeCXE, <italic>A. eriantha </italic>CXE; PepEST, <italic>Capsicum annuum </italic>Esterase, GenBank Accession No. <ext-link ext-link-type="gen" xlink:href="AAF77578.1">AAF77578.1</ext-link>; Gibberellin, Gibberellin receptor GID1, Q6L545; LeHSR203J, <italic>Lycopersicon esculentum </italic>HSR203J, BAA74434.1; NtHSR203J, <italic>Nicotiana tabacum </italic>HSR203J, AAF62404.1; PsHSR203j, <italic>Pisum sativum </italic>HSR203J, BAA85654.1; PsE86, <italic>P. sativum </italic>E86, <ext-link ext-link-type="gen" xlink:href="BAA85654.1">BAA85654.1</ext-link>; GeHIDM, <italic>Glycyrrhiza echinata </italic>HIDM, <ext-link ext-link-type="gen" xlink:href="BAD80839.1">BAD80839.1</ext-link>; GmHIDH, <italic>Glycine max </italic>HIDH, BAD80840.1; PrMC3, <italic>Pinus radiata </italic>MC3, <ext-link ext-link-type="gen" xlink:href="AAD04946.2">AAD04946.2</ext-link>; VvBIG8.1, <italic>Vitis vinifera </italic>BIG8.1, <ext-link ext-link-type="gen" xlink:href="&gt;AF48726-1">&gt;AF48726-1</ext-link>. Percentage bootstrap values (1000 bootstrap replicates) for groupings are given by each branch.</p></caption></fig> <fig position="float" id="F5"><label>Figure 5</label><caption><p>Distribution of <italic>Actinidia </italic>ESTs in the terpene biosynthesis pathway. Previously published sequences in GenBank belonging to the terpene biosynthetic pathway from mevalonate [##REF##11078528##43##] were used to identify genes in the <italic>Actinidia </italic>EST database using BLAST searches. ESTs refers to the number of ESTs found in the <italic>Actinidia </italic>EST database for each step of the pathway.</p></caption></fig> <fig position="float" id="F6"><label>Figure 6</label><caption><p>Distribution of <italic>Actinidia </italic>ESTs in the chlorophyll breakdown pathway. Previously published sequences in GenBank belonging to the chlorophyll breakdown pathway of higher plants [##REF##16669755##46##] were used to identify genes in the <italic>Actinidia </italic>EST database using BLAST searches. ESTs refers to the number of ESTs found in the <italic>Actinidia </italic>EST database for each step of the pathway. Where more than one number of ESTs is given, each number represents ESTs that matched a different gene for that step.</p></caption></fig> <fig position="float" id="F7"><label>Figure 7</label><caption><p>Distribution of <italic>Actinidia </italic>ESTs in the carotenoid biosynthetic pathway. Previously published sequences in GenBank belonging to the carotenoid biosynthetic pathway [##REF##15012246##47##] were used to identify genes in the <italic>Actinidia </italic>EST database using BLAST searches. ESTs refers to the number of ESTs found in the <italic>Actinidia </italic>EST database for each step of the pathway.</p></caption></fig> <fig position="float" id="F8"><label>Figure 8</label><caption><p>Distribution of <italic>Actinidia </italic>ESTs in the flavonoid biosynthetic pathway. Previously published sequences in GenBank belonging to the flavonoid biosynthetic pathway [##REF##16669781##49##] were used to identify genes in the <italic>Actinidia </italic>EST database using BLAST searches. Enzyme names are followed in brackets by the number of ESTs in different families of genes encoding these enzymes. See text for abbreviations.</p></caption></fig> <fig position="float" id="F9"><label>Figure 9</label><caption><p>Distribution of <italic>Actinidia </italic>ESTs in three ascorbic acid biosynthetic pathways. Previously published sequences in GenBank belonging to the three known ascorbic acid biosynthetic pathways in plants [##UREF##17##55##] were used to identify genes in the <italic>Actinidia </italic>EST database using BLAST searches. The number of ESTs found in the Actinidia EST database for each step of each pathway is shown. Steps where no ESTs are noted have not had genes identified.</p></caption></fig> <fig position="float" id="F10"><label>Figure 10</label><caption><p>Distribution of <italic>Actinidia </italic>ESTs in the quinate biosynthetic pathway. Previously published sequences in GenBank belonging to the quinate biosynthetic pathway [##REF##15012217##98##] were used to identify genes in the <italic>Actinidia </italic>EST database. The number of ESTs found in the <italic>Actinidia </italic>EST database for each step of each pathways is shown. Dehydroquinate dehydratase and shikimate dehydrogenase are two activities of a single bifunctional enzyme encoded by one gene [##REF##16784230##99##].</p></caption></fig> <fig position="float" id="F11"><label>Figure 11</label><caption><p>Phylogenetic tree of plant globulins of known function and <italic>Actinidia </italic>members of the globulin family. Previously published plant globulin sequences from GenBank were used to identify genes from the <italic>Actinidia </italic>EST database using BLAST searches. Abbreviations for species are as follows: Aa, <italic>Actinidia arguta</italic>; Ac, <italic>A. chinensis</italic>; Ad, <italic>A. deliciosa</italic>; At, <italic>Arabidopsis thaliana</italic>; Gm is <italic>Glycine max </italic>(soybean) followed by the Genbank accession number. Percentage bootstrap values (1000 bootstrap replicates) for groupings are given by each branch.</p></caption></fig> <fig position="float" id="F12"><label>Figure 12</label><caption><p>Distribution of <italic>Actinidia </italic>ESTs involved in postharvest kiwifruit softening. Schematic representation of postharvest kiwifruit softening in relation to the timing of key events in the softening process of ethylene-treated 'Hayward' fruit. Diagram modified from [##UREF##31##80##]). Previously published sequences in GenBank involved in postharvest softening processes were used to identify genes in the <italic>Actinidia </italic>EST database using BLAST searches. The first number in each set of brackets represents the total number of ESTs found in the <italic>Actinidia </italic>database and the second number is the number of ESTs found in fruit libraries.</p></caption></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p>Numbers of ESTs obtained from different <italic>Actinidia </italic>tissues by species. </p></caption><table frame="hsides" rules="groups"><thead><tr><td/><td align="center" colspan="6">Species</td></tr></thead><tbody><tr><td align="left">Tissue</td><td align="right"><italic>A. arguta</italic></td><td align="right"><italic>A. chinensis</italic></td><td align="right"><italic>A. deliciosa</italic></td><td align="right"><italic>A. eriantha</italic></td><td align="right">Other<sup>a</sup></td><td align="right">Total<sup>b</sup></td></tr><tr><td colspan="7"><hr/></td></tr><tr><td align="left">Bud</td><td align="right"><sup>c</sup></td><td align="right">15,689</td><td align="right">3,4519</td><td/><td/><td align="right">50,208</td></tr><tr><td align="left">Cell<sup>d</sup></td><td/><td align="right">4,851</td><td/><td/><td/><td align="right">4,851</td></tr><tr><td align="left">Fruit</td><td align="right">5,421</td><td align="right">8,453</td><td align="right">13,282</td><td align="right">11,259</td><td/><td align="right">38,415</td></tr><tr><td align="left">Leaf</td><td/><td align="right">17,325</td><td/><td/><td/><td align="right">17,325</td></tr><tr><td align="left">Petal</td><td align="right">1,836</td><td align="right">1,061</td><td align="right">9,950</td><td align="right">1,388</td><td align="right">1,422</td><td align="right">15,657</td></tr><tr><td align="left">Root</td><td/><td/><td/><td/><td align="right">5,101</td><td align="right">5,101</td></tr><tr><td align="left">Stem</td><td/><td/><td/><td/><td align="right">1,020</td><td align="right">1,020</td></tr><tr><td colspan="7"><hr/></td></tr><tr><td align="left">Total</td><td align="right">7,257</td><td align="right">47,379</td><td align="right">57,751</td><td align="right">12,647</td><td align="right">7,543</td><td align="right">132,577</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T2"><label>Table 2</label><caption><p>Comparative genomic hybridizations of different <italic>Actinidia </italic>species.</p></caption><table frame="hsides" rules="groups"><thead><tr><td/><td align="center"><italic>A. arguta</italic></td><td align="center"><italic>A. chinensis</italic></td><td align="center"><italic>A. deliciosa</italic></td><td align="center"><italic>A. eriantha</italic></td><td align="center"><italic>A. hemsleyana</italic></td><td align="center"><italic>A. indo-chinensis</italic></td><td align="center"><italic>A. poly-gama</italic></td><td align="center"><italic>A. setosa</italic></td></tr></thead><tbody><tr><td align="left">No. of oligos<sup>a</sup></td><td align="center">479<sup>b</sup></td><td align="center">8,957</td><td align="center">6,606</td><td align="center">978</td><td align="center">148</td><td align="center">12</td><td align="center">188</td><td align="center">104</td></tr><tr><td align="left">A</td><td align="center">2</td><td align="center">58</td><td align="center">95</td><td align="center">15</td><td align="center">0</td><td align="center">0</td><td align="center">1</td><td align="center">1</td></tr><tr><td colspan="9"><hr/></td></tr><tr><td align="left">B</td><td align="center">0</td><td align="center">12</td><td align="center">27</td><td align="center">5</td><td align="center">0</td><td align="center">0</td><td align="center">0</td><td align="center">0</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T3"><label>Table 3</label><caption><p>Functional classification of ESTs and NRs from <italic>Actinidia</italic>.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="center">Bin Code<break/></td><td align="left">Bin Name<break/></td><td align="center">% <italic>Actinidia </italic><break/>ESTs<sup>a</sup></td><td align="center">% <italic>Actinidia </italic>NR <break/>sequences<sup>b</sup></td><td align="center">% <italic>Arabidopsis </italic><break/>sequences</td><td align="center">% Potato <break/>NRs</td></tr></thead><tbody><tr><td align="center">1</td><td align="left">PS</td><td align="center">2.18</td><td align="center">1.22</td><td align="center">0.73</td><td align="center">0.75</td></tr><tr><td align="center">2</td><td align="left">major CHO metabolism</td><td align="center">0.86</td><td align="center">0.72</td><td align="center">0.37</td><td align="center">0.54</td></tr><tr><td align="center">3</td><td align="left">minor CHO metabolism</td><td align="center">1.02</td><td align="center">0.75</td><td align="center">0.45</td><td align="center">0.41</td></tr><tr><td align="center">4</td><td align="left">glycolysis</td><td align="center">1.03</td><td align="center">0.60</td><td align="center">0.23</td><td align="center">0.27</td></tr><tr><td align="center">5</td><td align="left">fermentation</td><td align="center">0.38</td><td align="center">0.21</td><td align="center">0.05</td><td align="center">0.07</td></tr><tr><td align="center">6</td><td align="left">Gluconeogenesis/glyoxylate cycle</td><td align="center">0.14</td><td align="center">0.11</td><td align="center">0.04</td><td align="center">0.02</td></tr><tr><td align="center">7</td><td align="left">OPP</td><td align="center">0.25</td><td align="center">0.22</td><td align="center">0.11</td><td align="center">0.12</td></tr><tr><td align="center">8</td><td align="left">TCA/org. transformation</td><td align="center">0.64</td><td align="center">0.55</td><td align="center">0.27</td><td align="center">0.27</td></tr><tr><td align="center">9</td><td align="left">mitochondrial electron transport/ATP synthesis</td><td align="center">0.82</td><td align="center">0.71</td><td align="center">0.49</td><td align="center">0.34</td></tr><tr><td align="center">10</td><td align="left">cell wall</td><td align="center">3.20</td><td align="center">2.27</td><td align="center">1.90</td><td align="center">1.11</td></tr><tr><td align="center">11</td><td align="left">lipid metabolism</td><td align="center">1.78</td><td align="center">1.91</td><td align="center">1.49</td><td align="center">1.32</td></tr><tr><td align="center">12</td><td align="left">N-metabolism</td><td align="center">0.08</td><td align="center">0.16</td><td align="center">0.09</td><td align="center">0.08</td></tr><tr><td align="center">13</td><td align="left">amino acid metabolism</td><td align="center">2.64</td><td align="center">2.00</td><td align="center">1.07</td><td align="center">1.18</td></tr><tr><td align="center">14</td><td align="left">S-assimilation</td><td align="center">0.10</td><td align="center">0.08</td><td align="center">0.05</td><td align="center">0.05</td></tr><tr><td align="center">15</td><td align="left">metal handling</td><td align="center">2.53</td><td align="center">0.59</td><td align="center">0.34</td><td align="center">0.20</td></tr><tr><td align="center">16</td><td align="left">secondary metabolism</td><td align="center">3.02</td><td align="center">2.09</td><td align="center">1.64</td><td align="center">1.42</td></tr><tr><td align="center">17</td><td align="left">hormone metabolism</td><td align="center">2.16</td><td align="center">1.99</td><td align="center">2.27</td><td align="center">1.69</td></tr><tr><td align="center">18</td><td align="left">co-factor and vitamin metabolism</td><td align="center">0.10</td><td align="center">0.14</td><td align="center">0.16</td><td align="center">0.15</td></tr><tr><td align="center">19</td><td align="left">tetrapyrrole synthesis</td><td align="center">0.28</td><td align="center">0.26</td><td align="center">0.17</td><td align="center">0.16</td></tr><tr><td align="center">20</td><td align="left">stress</td><td align="center">4.08</td><td align="center">3.33</td><td align="center">3.48</td><td align="center">2.64</td></tr><tr><td align="center">21</td><td align="left">redox. regulation</td><td align="center">1.02</td><td align="center">1.00</td><td align="center">0.70</td><td align="center">0.64</td></tr><tr><td align="center">22</td><td align="left">polyamine metabolism</td><td align="center">0.29</td><td align="center">0.12</td><td align="center">0.05</td><td align="center">0.07</td></tr><tr><td align="center">23</td><td align="left">nucleotide metabolism</td><td align="center">0.59</td><td align="center">0.68</td><td align="center">0.57</td><td align="center">0.49</td></tr><tr><td align="center">24</td><td align="left">biodegradation of xenobiotics</td><td align="center">0.04</td><td align="center">0.08</td><td align="center">0.09</td><td align="center">0.06</td></tr><tr><td align="center">25</td><td align="left">C1-metabolism</td><td align="center">0.19</td><td align="center">0.19</td><td align="center">0.14</td><td align="center">0.12</td></tr><tr><td align="center">26</td><td align="left">miscellaneous</td><td align="center">4.72</td><td align="center">4.59</td><td align="center">5.38</td><td align="center">3.69</td></tr><tr><td align="center">27</td><td align="left">RNA</td><td align="center">7.86</td><td align="center">9.86</td><td align="center">10.88</td><td align="center">6.82</td></tr><tr><td align="center">28</td><td align="left">DNA</td><td align="center">1.60</td><td align="center">2.05</td><td align="center">1.99</td><td align="center">1.25</td></tr><tr><td align="center">29</td><td align="left">protein</td><td align="center">18.61</td><td align="center">16.94</td><td align="center">13.31</td><td align="center">9.59</td></tr><tr><td align="center">30</td><td align="left">signalling</td><td align="center">3.42</td><td align="center">4.77</td><td align="center">4.46</td><td align="center">3.10</td></tr><tr><td align="center">31</td><td align="left">cell</td><td align="center">3.23</td><td align="center">3.04</td><td align="center">2.57</td><td align="center">1.68</td></tr><tr><td align="center">33</td><td align="left">development</td><td align="center">1.29</td><td align="center">1.58</td><td align="center">2.24</td><td align="center">1.33</td></tr><tr><td align="center">34</td><td align="left">transport</td><td align="center">3.62</td><td align="center">3.99</td><td align="center">3.49</td><td align="center">2.58</td></tr><tr><td align="center">35</td><td align="left">not assigned</td><td align="center">26.23</td><td align="center">31.21</td><td align="center">38.72</td><td align="center">55.79</td></tr></tbody></table></table-wrap>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional file 1</title><p>Additional Table 1. Expressed sequence tags (ESTs) sequenced in Angiosperm orders and families.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S2"><caption><title>Additional file 2</title><p>Additional Table 2a. Summary of library names, descriptions and statistics for sequenced <italic>Actinidia </italic>ESTs. Additional Table 2b. <italic>Actinidia </italic>libraries that were subtracted before sequencing to reduce redundant sequences.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S3"><caption><title>Additional file 3</title><p>Additional Table 3 Codon usage calculated using 707 full-length <italic>Actinidia </italic>cDNA sequences.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S4"><caption><title>Additional file 4</title><p>Additional Table 4. Source of oligonucleotides used to construct the <italic>Actinidia </italic>microarray slides by library and species</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S5"><caption><title>Additional file 5</title><p>Additional Table 5. Highly expressed ESTs in the <italic>Actinidia </italic>EST database.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S6"><caption><title>Additional file 6</title><p>Additional Table 6. Most frequent InterPro families found in <italic>Actinidia </italic>NRs.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S7"><caption><title>Additional file 7</title><p>Additional Table 7. Volatiles, acids, and sugars in <italic>Actinidia </italic>flowers and fruit</p></caption></supplementary-material>
<table-wrap-foot><p>The number of ESTs was calculated by amalgamating libraries from the same tissues.</p><p><sup>a </sup>Other miscellaneous species were <italic>Actinidia setosa </italic>(1,020 from stem), <italic>A. hemsleyana </italic>(5,101 from roots), <italic>A. polygama </italic>(1,348 from petal) and <italic>A. indochinensis </italic>(74 from petal).</p><p><sup>b </sup>See Additional file ##SUPPL##1##2a## for full details of libraries and tissues.</p><p><sup>c </sup>Empty cells had no ESTs sequenced.</p><p><sup>d </sup>Cell refers to cell culture.</p></table-wrap-foot> <table-wrap-foot><p>A # features where the signal is increased in <italic>A. deliciosa</italic>.</p><p>B Signal increased in <italic>A. eriantha</italic>.</p><p><sup>a </sup>The number of oligonucleotide primers on the microarray derived from each Actinidia species.</p><p><sup>b </sup>Oligonucleotide microarrays were hybridized with labeled genomic DNA from A. deliciosa or A. eriantha. The signals were normalized and selected as differentially changing between the two species using a threshold of two.</p></table-wrap-foot> <table-wrap-foot><p><sup>a </sup>Actinidia NRs were assigned to a Mapman (Thimm et al., 2004) bin according to their nearest Arabidopsis BLAST hit. ESTs belonging to the Actinidia NR were then assigned to the same bin.</p><p><sup>b </sup>Potato and Arabidopsis functional assignments were as described [##REF##14996223##27##,##REF##16649112##100##]</p></table-wrap-foot>
<graphic xlink:href="1471-2164-9-351-1"/> <graphic xlink:href="1471-2164-9-351-2"/> <graphic xlink:href="1471-2164-9-351-3"/> <graphic xlink:href="1471-2164-9-351-4"/> <graphic xlink:href="1471-2164-9-351-5"/> <graphic xlink:href="1471-2164-9-351-6"/> <graphic xlink:href="1471-2164-9-351-7"/> <graphic xlink:href="1471-2164-9-351-8"/> <graphic xlink:href="1471-2164-9-351-9"/> <graphic xlink:href="1471-2164-9-351-10"/> <graphic xlink:href="1471-2164-9-351-11"/> <graphic xlink:href="1471-2164-9-351-12"/>
<media xlink:href="1471-2164-9-351-S1.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-351-S2.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-351-S3.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-351-S4.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-351-S5.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-351-S6.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-351-S7.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media>
[{"surname": ["Ferguson"], "given-names": ["AR"], "article-title": ["The need for characterisation and evaluation of germplasm: kiwifruit as an example"], "source": ["Euphytica"], "year": ["2007"], "volume": ["154"], "fpage": ["371"], "lpage": ["382"]}, {"surname": ["Ferguson", "Huang"], "given-names": ["AR", "HW"], "article-title": ["Genetic resources of kiwifruit: domestication and breeding"], "source": ["Hort Rev"], "year": ["2007"], "volume": ["33"], "fpage": ["1"], "lpage": ["121"]}, {"surname": ["Nishiyama", "Taylor SL"], "given-names": ["I"], "article-title": ["Fruits of the "], "italic": ["Actinidia"], "source": ["Advances in Food and Nutrition Research Vol 52"], "year": ["2007"], "volume": ["52"], "publisher-name": [" Academic Press"], "fpage": ["293"], "lpage": ["324"]}, {"surname": ["Atkinson", "MacRae", "Pua EC, Davey MR"], "given-names": ["RG", "EA"], "article-title": ["Kiwifruit"], "source": ["Transgenic Crops V Biotechnology in Agriculture and Forestry"], "year": ["2007"], "volume": ["60"], "publisher-name": ["Berlin, Heidelberg, Germany  , Springer-Verlag"], "fpage": ["329"], "lpage": ["346"]}, {"surname": ["Praekelt", "McKee", "Smith"], "given-names": ["UM", "RA", "H"], "article-title": ["Molecular analysis of actinidin, the cysteine proteinase of "], "italic": ["Actinidia chinensis"], "source": ["Plant Mol Biol"], "year": ["1988"], "volume": ["10"], "fpage": ["193"], "lpage": ["202"]}, {"surname": ["Nieuwenhuizen", "Beuning", "Sutherland", "Sharma", "Cooney", "Bieleski", "Schr\u00f6der", "MacRae", "Atkinson"], "given-names": ["NJ", "LL", "PW", "NN", "JM", "LRF", "R", "EA", "RG"], "article-title": ["Identification and characterisation of acidic and novel basic forms of actinidin, the highly abundant cysteine protease from kiwifruit"], "source": ["Funct Plant Biol"], "year": ["2007"], "volume": ["34"], "fpage": ["946"], "lpage": ["961"]}, {"surname": ["Malone", "Todd", "Burgess", "Philip", "Christeller"], "given-names": ["LA", "JH", "EPJ", "BA", "JT"], "article-title": ["Effects of kiwifruit ("], "italic": ["Actinidia deliciosa", "Lepidoptera: Noctuidae"], "source": ["NZ J Crop Hort Sci"], "year": ["2005"], "volume": ["33"], "fpage": ["99\u2013105"]}, {"surname": ["Hyeon", "Isoe", "Sakan"], "given-names": ["SB", "S", "T"], "article-title": ["The structure of neomatatabiol, the potent attractant for Chrysopa from "], "italic": ["Actinidia polygama"], "source": ["Tet Lett"], "year": ["1968"], "volume": [" 51"], "fpage": ["5325"], "lpage": ["5326"]}, {"surname": ["Yoshihara", "Sakai", "Sakan"], "given-names": ["K", "T", "T"], "article-title": ["Dehydroiridodial, the pungent principle of "], "italic": ["Actinidia polygama"], "source": ["Chem Lett"], "year": ["1978"], "volume": ["7"], "fpage": ["433"], "lpage": ["434"]}, {"surname": ["Burdock"], "given-names": ["GA"], "source": ["Fenaroli's Handbook of Flavor Ingredients"], "year": ["2005"], "publisher-name": ["Boca Raton , CRC Press"]}, {"surname": ["Paterson", "MacRae", "Young"], "given-names": ["VJ", "EA", "H"], "article-title": ["Relationships between sensory properties and chemical composition of kiwifruit ("], "italic": ["Actinidia deliciosa"], "source": ["J Sci Food Agric"], "year": ["1991"], "volume": ["57"], "fpage": ["235"], "lpage": ["251"]}, {"surname": ["Shaw", "Allen", "Yates", "Franich"], "given-names": ["GJ", "JM", "MK", "FA"], "article-title": ["Volatile flavour constituents of feijoa ("], "italic": ["Feijoa sellowiana"], "source": ["J Sci Food Agric"], "year": ["1990"], "volume": ["50"], "fpage": ["357"], "lpage": ["361"]}, {"surname": ["Flath", "Black", "Guadagni", "McFadden", "Schultz"], "given-names": ["RA", "DR", "DG", "WH", "TH"], "article-title": ["Identification and organoleptic evaluation of compounds in Delicious apple essence "], "source": ["J Agric Food Chem"], "year": ["1967"], "volume": ["15"], "fpage": ["29"], "lpage": ["35"]}, {"surname": ["Vollhardt", "Schore"], "given-names": ["KPC", "NE"], "source": [" Organic Chemistry"], "year": ["1994"], "edition": ["Second"], "publisher-name": ["New York , W. H. Freeman and Company"]}, {"surname": ["Green", "Friel", "Beuning", "Matich", "MacRae"], "given-names": ["S", "EN", "L", "A", "E"], "article-title": ["A multifunctional germacrene D synthase. Patent # WO04058814A1"], "year": ["2004"]}, {"surname": ["Ikegami", "Eguchi", "Kitajima", "Inoue", "Yonemori"], "given-names": ["A", "S", "A", "K", "K"], "article-title": ["Identification of genes involved in proanthocyanidin biosynthesis of persimmon ("], "italic": ["Diospyros kaki"], "source": ["Plant Sci"], "year": ["2007"], "volume": ["172"], "fpage": ["1037"], "lpage": ["1047"]}, {"surname": ["Ferguson", "MacRae"], "given-names": ["AR", "EA"], "article-title": ["Vitamin C in "], "italic": ["Actinidia"], "source": ["Acta Hortic"], "year": ["1992"], "volume": ["297"], "publisher-name": ["Wageningen , International Society for Horticultural Science"], "fpage": ["481"], "lpage": ["487"]}, {"surname": ["Ishikawa", "Dowdle", "Smirnoff"], "given-names": ["T", "J", "N"], "article-title": ["Progress in manipulating ascorbic acid biosynthesis and accumulation in plants"], "source": ["Physiol Plant"], "year": ["2006"], "volume": ["126"], "fpage": ["343"], "lpage": ["355"]}, {"surname": ["Boldingh", "Smith", "Klages"], "given-names": ["H", "GS", "K"], "article-title": ["Seasonal concentrations of non-structural carbohydrates of five "], "italic": ["Actinidia"], "source": ["Ann Bot"], "year": ["2000"], "volume": ["85"], "fpage": ["469"], "lpage": ["476"]}, {"surname": ["Klages", "Donnison", "Boldingh", "MacRae"], "given-names": ["K", "H", "H", "E"], "article-title": ["myo-Inositol is the major sugar in "], "italic": ["Actinidia arguta"], "source": ["Aust J Plant Physiol"], "year": ["1998"], "volume": ["25"], "fpage": ["61 "], "lpage": ["668"]}, {"surname": ["Walton", "De Jong"], "given-names": ["EF", "TM"], "article-title": ["Growth and compositional changes in kiwifruit berries from three Californian locations"], "source": ["Ann Bot"], "year": ["1990"], "volume": ["66"], "fpage": ["285"], "lpage": ["298"]}, {"surname": ["Cheng", "Seal", "Boldingh", "Marsh", "MacRae", "Murphy", "Ferguson"], "given-names": ["C", "A", "H", "K", "EA", "S", "AR"], "article-title": ["Inheritance of fruit characters and fruit size in a diploid "], "italic": ["Actinidia Chinensis"], "source": ["Euphytica"], "year": ["2004"], "volume": ["138"], "fpage": ["185"], "lpage": ["195"]}, {"surname": ["Dirlewanger", "Moing", "Rothan", "Svanella", "Pronier", "Guye", "Plomion", "Monet"], "given-names": ["E", "A", "C", "L", "V", "A", "C", "R"], "article-title": ["Mapping QTLs controlling fruit quality in peach ("], "italic": ["Prunus persica"], "source": ["Theoret Appl Genet"], "year": ["1999"], "volume": ["98"], "fpage": ["18"], "lpage": ["31"]}, {"surname": ["Ossipov", "Bonner", "Ossipova", "Jensen"], "given-names": ["V", "C", "S", "R"], "article-title": ["Broad-specificity quinate (shikimate) dehydrogenase from "], "italic": ["Pinus taeda"], "source": ["Plant Physiol Biochem"], "year": ["2000"], "volume": ["38"], "fpage": ["923"], "lpage": ["928"]}, {"surname": ["Mills", "Shewry"], "given-names": ["CE", "PR"], "source": ["Plant food allergens"], "year": ["2004"], "publisher-name": ["Oxford , Blackwell Science "]}, {"surname": ["Moller", "Kayma", "Steinhart", "Paschke"], "given-names": ["M", "M", "H", "A"], "article-title": ["Isolation and characterization of a major allergen in kiwi fruit"], "source": ["Z Lebensm Unters Forsch"], "year": ["1997"], "volume": [" 205"], "fpage": ["364\u2013369"]}, {"surname": ["Singh", "Shepherd", "Langridge", "Gruen"], "given-names": ["NK", "KW", "P", "LC"], "article-title": ["Purification and biochemical characterization of triticin, a legumin-like protein in wheat endosperm"], "source": ["J Cereal Sci"], "year": ["1991"], "volume": ["13"], "fpage": ["207"], "lpage": ["219"]}, {"surname": ["Muntz", "Gueguen J, Popineau Y"], "given-names": ["K"], "article-title": ["Globulins from legume seeds: Structure and function during storage and reactivation"], "source": ["Plant proteins from European crops : food and non-food applications"], "year": ["1998"], "publisher-name": ["Berlin ; New York , Springer, c1998. xxvi, 339 p. : ill."], "fpage": ["3"], "lpage": ["12"]}, {"surname": ["Shutov", "Vaintraub"], "given-names": ["AD", "IA"], "article-title": ["Degradation of storage proteins in germinating seeds"], "source": ["Phytochem"], "year": ["1987"], "volume": ["26"], "fpage": ["1557"], "lpage": ["1566"]}, {"surname": ["White", "de Silva", "Requejo-Tapia", "Harker"], "given-names": ["A", "HN", "C", "FR"], "article-title": ["Evaluation of softening characteristics of fruit from 14 species of Actinidia"], "source": ["Postharv Biol and Tech"], "year": ["2004"], "volume": ["35"], "fpage": ["143\u2013151."]}, {"surname": ["Jackson", "Harker"], "given-names": ["PJ", "FR"], "article-title": ["Changes in the firmness of the outer pericarp, inner pericarp, and core of Actinidia species during ripening"], "source": ["NZ J Crop Hort Sci"], "year": ["1997"], "volume": ["25"], "fpage": ["185\u2013189"]}, {"surname": ["Schr\u00f6der", "Atkinson"], "given-names": ["R", "RG"], "article-title": ["Kiwifruit cell walls: towards an understanding of softening?"], "source": ["NZ J For Sci"], "year": ["2006"], "volume": ["36"], "fpage": ["112"], "lpage": ["129"]}, {"surname": ["Ross", "Redgwell", "MacRae"], "given-names": ["GS", "RJ", "EA"], "article-title": ["Kiwifruit \u03b2-galactosidase: Isolation and activity against specific fruit cell-wall polysaccharides."], "source": ["Planta"], "year": ["1993"], "volume": ["189"], "fpage": ["499\u2013506"]}, {"surname": ["Gill", "Harvey", "Gardner", "Fraser"], "given-names": ["GP", "CF", "RC", "LG"], "article-title": ["Development of sex-linked PCR markers for gender identification in Actinidia"], "source": ["Theoret Appl Genet"], "year": ["1998"], "volume": ["97"], "fpage": ["439"], "lpage": ["445"]}, {"surname": ["Chang", "Puryear", "Cairney"], "given-names": ["S", "J", "J"], "article-title": ["A simple and efficient method for isolating RNA from pine trees"], "source": ["Plant Mol Biol Rep"], "year": ["1993"], "volume": ["11"], "fpage": ["113"], "lpage": ["116"]}, {"surname": ["Lopez-Gomez", "Gomez-Lim"], "given-names": ["R", "MA"], "article-title": ["A method for extracting intact RNA from fruits rich in polysaccharides using ripe mango mesocarp"], "source": ["HortSci"], "year": ["1992"], "volume": ["27"], "fpage": ["440"], "lpage": ["442"]}, {"surname": ["Zdobnov", "Apweiler"], "given-names": ["EM", "R"], "article-title": ["InterProScan--an integration platform for the signature-recognition methods in InterPro"], "source": ["Bioinformat"], "year": ["2001"], "volume": ["17"], "fpage": ["847"], "lpage": ["848"]}, {"surname": ["Felsenstein"], "given-names": ["J"], "source": ["PHYLIP (Phylogeny Inference Package), Version 3.6a2"], "year": ["1993"], "edition": [" 3.6a2"], "publisher-name": [" Department of Genetics, University of Washington, Seattle."]}, {"surname": ["Richardson", "Marsh", "Boldingh", "Pickering", "Bulley", "Frearson", "Ferguson", "Thornber", "Bolitho", "Macrae"], "given-names": ["AC", "KB", "HL", "AH", "SM", "NJ", "AR", "SE", "KM", "EA"], "article-title": ["High growing temperatures reduce fruit carbohydrate and vitamin C in kiwifruit"], "source": ["Plant, Cell & Environment"], "year": ["2004"], "volume": ["27"], "fpage": ["423"], "lpage": ["435"]}]
BMC Genomics. 2008 Jul 27; 9:351