Analysis of grape ESTs: global gene expression patterns in leaf and berry

Swift metabolite changes and leaf shedding are milestones in the acclimation process of grapevine under prolonged water stress

BACKGROUND

Grape leaves provide the biochemical substrates for berry development. Thus, understanding the regulation of grapevine leaf metabolism can aid in discerning processes fundamental to fruit development and berry quality. Here, the temporal alterations in leaf metabolism in Merlot grapevine grown under sufficient irrigation and water deficit were monitored from veraison until harvest.

RESULTS

The vines mediated water stress gradually and involving multiple strategies: osmotic adjustment, transcript-metabolite alteration and leaf shedding. Initially stomatal conductance and leaf water potential showed a steep decrease together with the induction of stress related metabolism, e.g. up-regulation of proline and GABA metabolism and stress related sugars, and the down-regulation of developmental processes. Later, progressive soil drying was associated with an incremental contribution of Ca²⁺ and sucrose to the osmotic adjustment concomitant with the initiation of leaf shedding. Last, towards harvest under progressive stress conditions following leaf shedding, incremental changes in leaf water potential were measured, while the magnitude of perturbation in leaf metabolism lessened.

CONCLUSIONS

The data present evidence that over time grapevine acclimation to water stress diversifies in temporal responses encompassing the alteration of central metabolism and gene expression, osmotic adjustments and reduction in leaf area. Together these processes mitigate leaf water stress and aid in maintaining the berry-ripening program.

Proteomic and metabolic traits of grape exocarp to explain different anthocyanin concentrations of the cultivars

The role of grape berry skin as a protective barrier against damage by physical injuries and pathogen attacks requires a metabolism able to sustain biosynthetic activities such as those relating to secondary compounds (i.e., flavonoids). In order to draw the attention on these biochemical processes, a proteomic and metabolomic comparative analysis was performed among Riesling Italico, Pinot Gris, Pinot Noir, and Croatina cultivars, which are known to accumulate anthocyanins to a different extent. The application of multivariate statistics on the dataset pointed out that the cultivars were distinguishable from each other and the order in which they were grouped mainly reflected their relative anthocyanin contents. Sorting the spots according to their significance 100 proteins were characterized by LC-ESI-MS/MS. Through GC-MS, performed in Selected Ion Monitoring (SIM) mode, 57 primary metabolites were analyzed and the differences in abundance of 16 of them resulted statistically significant to ANOVA test. Considering the functional distribution, the identified proteins were involved in many physiological processes such as stress, defense, carbon metabolism, energy conversion and secondary metabolism. The trends of some metabolites were related to those of the protein data. Taken together, the results permitted to highlight the relationships between the secondary compound pathways and the main metabolism (e.g., glycolysis and TCA cycle). Moreover, the trend of accumulation of many proteins involved in stress responses, reinforced the idea that they could play a role in the cultivar specific developmental plan.

Construction of cDNA library and preliminary analysis of expressed sequence tags from tea plant [Camellia sinensis (L) O. Kuntze]

Tea is the most popular non-alcoholic and healthy beverage across the world. The understanding of the genetic organization and molecular biology of tea plant, which is very poorly understood at present, is required for quantum increase in productivity and efficient use of germplasm for either cultivation or breeding program. Single-pass sequencing of randomly selected cDNA clones is the most widely accepted technique for gene identification and cloning. In the present study, a good quality cDNA library was constructed and preliminary analysis of ESTs was carried out. The titers of unamplified and amplified libraries were 1.4 × 10(6)pfu/ml and 5.27 × 10(8)pfu/ml respectively. A total of 210 cDNA clones from the constructed cDNA library were sequenced and analyzed. A total of 84 high quality Expressed Sequence Tags (ESTs) were generated, among which 71 ESTs had significant homology with sequences in NCBI non-redundant protein database by BLAST X analysis. About 80% ESTs had poly (A) tail at 3' end indicating that the cDNAs were full length. The database-matched ESTs were classified into putative cellular roles, viz. energy-related category (corresponding to 20% of total BLAST X matched ESTs), Transcription (14.2%), protein synthesis (14.2%) cell growth and division (8.6%), cell structure (5.7%), signal transduction (5.7%), transporters (2.9%), disease and defenses (2.9%), secondary metabolism (2.9%) and gene regulation (2.9%). This study provides an overview of the mRNA expression profile and first hand information of gene sequence expressed in tender leaves and apical buds of tea plant.

The first set of EST resource for gene discovery and marker development in pigeonpea (Cajanus cajan L.)

BACKGROUND

Pigeonpea (Cajanus cajan (L.) Millsp) is one of the major grain legume crops of the tropics and subtropics, but biotic stresses [Fusarium wilt (FW), sterility mosaic disease (SMD), etc.] are serious challenges for sustainable crop production. Modern genomic tools such as molecular markers and candidate genes associated with resistance to these stresses offer the possibility of facilitating pigeonpea breeding for improving biotic stress resistance. Availability of limited genomic resources, however, is a serious bottleneck to undertake molecular breeding in pigeonpea to develop superior genotypes with enhanced resistance to above mentioned biotic stresses. With an objective of enhancing genomic resources in pigeonpea, this study reports generation and analysis of comprehensive resource of FW- and SMD- responsive expressed sequence tags (ESTs).

RESULTS

A total of 16 cDNA libraries were constructed from four pigeonpea genotypes that are resistant and susceptible to FW ('ICPL 20102' and 'ICP 2376') and SMD ('ICP 7035' and 'TTB 7') and a total of 9,888 (9,468 high quality) ESTs were generated and deposited in dbEST of GenBank under accession numbers GR463974 to GR473857 and GR958228 to GR958231. Clustering and assembly analyses of these ESTs resulted into 4,557 unique sequences (unigenes) including 697 contigs and 3,860 singletons. BLASTN analysis of 4,557 unigenes showed a significant identity with ESTs of different legumes (23.2-60.3%), rice (28.3%), Arabidopsis (33.7%) and poplar (35.4%). As expected, pigeonpea ESTs are more closely related to soybean (60.3%) and cowpea ESTs (43.6%) than other plant ESTs. Similarly, BLASTX similarity results showed that only 1,603 (35.1%) out of 4,557 total unigenes correspond to known proteins in the UniProt database (or= 5 sequences detected 102 single nucleotide polymorphisms (SNPs) in 37 contigs. As an example, a set of 10 contigs were used for confirming in silico predicted SNPs in a set of four genotypes using wet lab experiments. Occurrence of SNPs were confirmed for all the 6 contigs for which scorable and sequenceable amplicons were generated. PCR amplicons were not obtained in case of 4 contigs. Recognition sites for restriction enzymes were identified for 102 SNPs in 37 contigs that indicates possibility of assaying SNPs in 37 genes using cleaved amplified polymorphic sequences (CAPS) assay.

CONCLUSION

The pigeonpea EST dataset generated here provides a transcriptomic resource for gene discovery and development of functional markers associated with biotic stress resistance. Sequence analyses of this dataset have showed conservation of a considerable number of pigeonpea transcripts across legume and model plant species analysed as well as some putative pigeonpea specific genes. Validation of identified biotic stress responsive genes should provide candidate genes for allele mining as well as candidate markers for molecular breeding.

A framework map from grapevine V3125 (Vitis vinifera 'Schiava grossa' x 'Riesling') x rootstock cultivar 'Börner' (Vitis riparia x Vitis cinerea) to localize genetic determinants of phylloxera root resistance

Grapevine rootstock cultivar 'Börner' is a hybrid of Vitis riparia and Vitis cinerea Arnold that shows high resistance to phylloxera (Daktulosphaira vitifoliae Fitch). To localize the determinants of phylloxera root resistance, the susceptible grapevine V3125 (Vitis vinifera 'Schiava grossa' x 'Riesling') was crossed to 'Börner'. Genetic framework maps were built from the progeny. 235 microsatellite markers were placed on the integrated parental map. They cover 1,155.98 cM on 19 linkage groups with an average marker distance of 4.8 cM. Phylloxera resistance was scored by counting nodosities after inoculation of the root system. Progeny plants were triplicated and experimentally infected in 2 years. A scan of the genetic maps indicated a quantitative trait locus on linkage group 13. This region was targeted by six microsatellite-type markers newly developed from the V. vinifera model genome sequence. Two of these appear closely linked to the trait, and can be useful for marker-assisted breeding.

Crotalus durissus collilineatus venom gland transcriptome: analysis of gene expression profile

Crotalus durissus rattlesnakes are responsible for the most lethal cases of snakebites in Brazil. Crotalus durissus collilineatus subspecies is related to a great number of accidents in Southeast and Central West regions, but few studies on its venom composition have been carried out to date. In an attempt to describe the transcriptional profile of the C. durissus collilineatus venom gland, we generated a cDNA library and the sequences obtained could be identified by similarity searches on existing databases. Out of 673 expressed sequence tags (ESTs) 489 produced readable sequences comprising 201 singletons and 47 clusters of two or more ESTs. One hundred and fifty reads (60.5%) produced significant hits to known sequences. The results showed a predominance of toxin-coding ESTs instead of transcripts coding for proteins involved in all cellular functions. The most frequent toxin was crotoxin, comprising 88% of toxin-coding sequences. Crotoxin B, a basic phospholipase A(2) (PLA(2)) subunit of crotoxin, was represented in more variable forms comparing to the non-enzymatic subunit (crotoxin A), and most sequences coding this molecule were identified as CB1 isoform from Crotalus durissus terrificus venom. Four percent of toxin-related sequences in this study were identified as growth factors, comprising five sequences for vascular endothelial growth factor (VEGF) and one for nerve growth factor (NGF) that showed 100% of identity with C. durissus terrificus NGF. We also identified two clusters for metalloprotease from PII class comprising 3% of the toxins, and two for serine proteases, including gyroxin (2.5%). The remaining 2.5% of toxin-coding ESTs represent singletons identified as homologue sequences to cardiotoxin, convulxin, angiotensin-converting enzyme inhibitor and C-type natriuretic peptide, Ohanin, crotamin and PLA(2) inhibitor. These results allowed the identification of the most common classes of toxins in C. durissus collilineatus snake venom, also showing some unknown classes for this subspecies and even for C. durissus species, such as cardiotoxins and VEGF.

A blackberry (Rubus L.) expressed sequence tag library for the development of simple sequence repeat markers

BACKGROUND

The recent development of novel repeat-fruiting types of blackberry (Rubus L.) cultivars, combined with a long history of morphological marker-assisted selection for thornlessness by blackberry breeders, has given rise to increased interest in using molecular markers to facilitate blackberry breeding. Yet no genetic maps, molecular markers, or even sequences exist specifically for cultivated blackberry. The purpose of this study is to begin development of these tools by generating and annotating the first blackberry expressed sequence tag (EST) library, designing primers from the ESTs to amplify regions containing simple sequence repeats (SSR), and testing the usefulness of a subset of the EST-SSRs with two blackberry cultivars.

RESULTS

A cDNA library of 18,432 clones was generated from expanding leaf tissue of the cultivar Merton Thornless, a progenitor of many thornless commercial cultivars. Among the most abundantly expressed of the 3,000 genes annotated were those involved with energy, cell structure, and defense. From individual sequences containing SSRs, 673 primer pairs were designed. Of a randomly chosen set of 33 primer pairs tested with two blackberry cultivars, 10 detected an average of 1.9 polymorphic PCR products.

CONCLUSION

This rate predicts that this library may yield as many as 940 SSR primer pairs detecting 1,786 polymorphisms. This may be sufficient to generate a genetic map that can be used to associate molecular markers with phenotypic traits, making possible molecular marker-assisted breeding to compliment existing morphological marker-assisted breeding in blackberry.

Expressed sequence tag analysis of the response of apple (Malus x domestica'Royal Gala') to low temperature and water deficit

Leaf, bark, xylem and root tissues were used to make nine cDNA libraries from non-stressed (control) 'Royal Gala' apple trees, and from 'Royal Gala' trees exposed to either low temperature (5 degrees C for 24 h) or water deficit (45% of saturated pot mass for 2 weeks). Over 22 600 clones from the nine libraries were subjected to 5' single-pass sequencing, clustered and annotated using blastx. The number of clusters in the libraries ranged from 170 to 1430. Regarding annotation of the sequences, blastx analysis indicated that within the libraries 65-72% of the clones had a high similarity to known function genes, 6-15% had no functional assignment and 15-26% were completely novel. The expressed sequence tags were combined into three classes (control, low-temperature and water deficit) and the annotated genes in each class were placed into 1 of 10 different functional categories. The percentage of genes falling into each category was then calculated. This analysis indicated a distinct downregulation of genes involved in general metabolism and photosynthesis, while a significant increase in defense/stress-related genes, protein metabolism and energy was observed. In particular, there was a three-fold increase in the number of stress genes observed in the water deficit libraries indicating a major shift in gene expression in response to a chronic stress. The number of stress genes in response to low temperature, although elevated, was much less than the water deficit libraries perhaps reflecting the shorter (24 h) exposure to stress. Genes with greater than five clones in any specific library were identified and, based on the number of clones obtained, the fold increase or decrease in expression in the libraries was calculated and verified by semiquantitative polymerase chain reaction. Genes, of particular note, that code for the following proteins were overexpressed in the low-temperature libraries: dehydrin and metallothionein-like proteins, ubiquitin proteins, a dormancy-associated protein, a plasma membrane intrinsic protein and an RNA-binding protein. Genes that were upregulated in the water deficit libraries fell mainly into the functional categories of stress (heat shock proteins, dehydrins) and photosynthesis. With few exceptions, the overall differences in downregulated genes were nominal compared with differences in upregulated genes. The results of this apple study are similar to other global studies of plant response to stress but offer a more detailed analysis of specific tissue response (bark vs xylem vs leaf vs root) and a comparison between an acute stress (24-h exposure to low temperature) and a chronic stress (2 weeks of water deficit).

Genome-wide transcriptional analysis of grapevine berry ripening reveals a set of genes similarly modulated during three seasons and the occurrence of an oxidative burst at vèraison

Background

Grapevine (Vitis species) is among the most important fruit crops in terms of cultivated area and economic impact. Despite this relevance, little is known about the transcriptional changes and the regulatory circuits underlying the biochemical and physical changes occurring during berry development.

Results

Fruit ripening in the non-climacteric crop species Vitis vinifera L. has been investigated at the transcriptional level by the use of the Affymetrix Vitis GeneChip^® which contains approximately 14,500 unigenes. Gene expression data obtained from berries sampled before and after véraison in three growing years, were analyzed to identify genes specifically involved in fruit ripening and to investigate seasonal influences on the process. From these analyses a core set of 1477 genes was found which was similarly modulated in all seasons. We were able to separate ripening specific isoforms within gene families and to identify ripening related genes which appeared strongly regulated also by the seasonal weather conditions. Transcripts annotation by Gene Ontology vocabulary revealed five overrepresented functional categories of which cell wall organization and biogenesis, carbohydrate and secondary metabolisms and stress response were specifically induced during the ripening phase, while photosynthesis was strongly repressed. About 19% of the core gene set was characterized by genes involved in regulatory processes, such as transcription factors and transcripts related to hormonal metabolism and signal transduction. Auxin, ethylene and light emerged as the main stimuli influencing berry development. In addition, an oxidative burst, previously not detected in grapevine, characterized by rapid accumulation of H₂O₂ starting from véraison and by the modulation of many ROS scavenging enzymes, was observed.

Conclusion

The time-course gene expression analysis of grapevine berry development has identified the occurrence of two well distinct phases along the process. The pre-véraison phase represents a reprogramming stage of the cellular metabolism, characterized by the expression of numerous genes involved in hormonal signalling and transcriptional regulation. The post-véraison phase is characterized by the onset of a ripening-specialized metabolism responsible for the phenotypic traits of the ripe berry. Between the two phases, at véraison, an oxidative burst and the concurrent modulation of the anti-oxidative enzymatic network was observed. The large number of regulatory genes we have identified represents a powerful new resource for dissecting the mechanisms of fruit ripening control in non-climacteric plants.

Analysis of protein changes during grape berry ripening by 2-DE and MALDI-TOF

Grape berry, a nonclimacteric fruit, during ripening turns from green, hard and acidic to coloured, soft and sweet. Many studies have focused on dynamic changes of mRNA levels, metabolites, sugars or individual proteins, but this is the first report of a proteomic approach applied to the screening of the most prominent variations that take place during berry ripening. Vitis vinifera cv. 'Nebbiolo Lampia' berries were collected at 10-day intervals, starting 1 month after flowering to complete ripe stage; total protein extracts from deseeded berries were separated by 2-DE. A total of 730 spots were detected in the 2-DE gels. 118 protein spots, differentially expressed during berry development, were subjected to MALDI-TOF analysis. Ninety-three of them were identified, corresponding to 101 proteins. The majority of proteins were linked to metabolism, energy and protein synthesis and fate. In comparison to published surveys of major berry proteins, fewer proteins related to stress response and more proteins related to cell structure were differentially expressed. Our data confirm a general decrease of glycolysis during ripening, and an increase of PR proteins in the range of 20-35 kDa. They furthermore suggest that oxidative stress decreases during ripening while extensive cytoskeleton rearrangement takes place in this period.

Generation of ESTs in Vitis vinifera wine grape (Cabernet Sauvignon) and table grape (Muscat Hamburg) and discovery of new candidate genes with potential roles in berry development

We report the generation and analysis of a total of 77,583 expressed sequence tags (ESTs) from two grapevine (Vitis vinifera L.) cultivars, Cabernet Sauvignon (wine grape) and Muscat Hamburg (table grape) with a focus on EST sequence quality and assembly optimization. The majority of the ESTs were derived from normalized cDNA libraries representing berry pericarp and seed developmental series, pooled non-berry tissues including root, flower, and leaf in Cabernet Sauvignon, and pooled tissues of berry, seed, and flower in Muscat Hamburg. EST and unigene sequence quality were determined by computational filtering coupled with small-scale contig reassembly, manual review, and BLAST analyses. EST assembly was optimized to better discriminate among closely related paralogs using two independent grape sequence sets, a previously published set of Vitis spp. gene families and our EST dataset derived from pooled leaf, flower, and root tissues of Cabernet Sauvignon. Sequence assembly within individual libraries indicated that those prepared from pooled tissues contributed the most to gene discovery. Annotations based upon searches against multiple databases including tomato and strawberry sequences helped to identify putative functions of ESTs and unigenes, particularly with respect to fleshy fruit development. Sequence comparison among the three wine grape libraries identified a number of genes preferentially expressed in the pericarp tissue, including transcription factors, receptor-like protein kinases, and hexose transporters. Gene ontology (GO) classification in the biological process aspect showed that GO categories corresponding to 'transport' and 'cell organization and biogenesis', which are associated with metabolite movement and cell wall structural changes during berry ripening, were higher in pericarp than in other tissues in the wine grape studied. The sequence data were used to characterize potential roles of new genes in berry development and composition.

EST analysis of the heading leaf of Chinese cabbage (Brassica rapa L. ssp. pekinensis) in the early phase of the heading stage

A cDNA library was constructed from the heading leaf in the early phase of the heading stage of Chinese cabbage (Brassica rapa L. ssp. pekinensis). By sequencing the randomly selected clones, 1363 sequences longer than 200 bp were found, with better trace data. After removing the poly(A) and contamination sequences, 1162 ESTs longer than 150 bp were obtained, of which 1102 shared significant similarity with known sequences in protein and nucleotide databases of the National Center for Biotechnology Information (NCBI) as revealed by searches using the BLASTX and BLASTN engines. Functional assignment of the ESTs was based on the method used in the Arabidopsis thaliana genome-sequencing project. About 77% of the putative protein sequences with known biological functions best matched with those of A. thaliana deposited in the non-redundant database of NCBI. These data suggest that Chinese cabbage is closely related to A. thaliana. This result is different from that reported in other Brassica species. At nucleotide level, however, 51% of the ESTs were homologous to those deposited for A. thaliana when all ESTs were searched against the est-others database. In addition, 60 ESTs had no homology with any of the plant gene sequences deposited in GenBank. These ESTs are very important for understanding the unique developmental process of Chinesecabbage and elaborating its genetic mapping. Among the genes with assigned functions, the most abundant representatives were those involved in protein synthesis and energy metabolism. With the 1162 ESTs, 895 non-redundant contigs were generated after being aligned using the Seqman II module of DNAStar software at the threshold of more than 80% homology over a minimum of 40 base pairs. Of these, 723 were singletons containing only one EST sequence, indicating that many kinds of such genes are expressed in the heading leaf of Chinese cabbage. An expression profile of Chinese cabbage heading leaf with the 1162 ESTs was therefore acquired in this work. This could be very useful for uncovering the mechanism of the heading process, which is the most obvious characteristic of Chinese cabbage and perhaps other related species, such as Brassica oleracea. This work could accelerate the finding and characterization of genes specifically expressed in the heading stage of Chinese cabbage.

Major differences observed in transcript profiles of blueberry during cold acclimation under field and cold room conditions

Our laboratory has been working toward increasing our understanding of the genetic control of cold hardiness in blueberry (Vaccinium section Cyanococcus) to ultimately use this information to develop more cold hardy cultivars for the industry. Here, we report using cDNA microarrays to monitor changes in gene expression at multiple times during cold acclimation under field and cold room conditions. Microarrays contained over 2,500 cDNA inserts, approximately half of which had been picked and single-pass sequenced from each of two cDNA libraries that were constructed from cold acclimated floral buds and non-acclimated floral buds of the fairly cold hardy cv. Bluecrop (Vaccinium corymbosum L.). Two biological samples were examined at each time point. Microarray data were analyzed statistically using t tests, ANOVA, clustering algorithms, and online analytical processing (OLAP). Interestingly, more transcripts were found to be upregulated under cold room conditions than under field conditions. Many of the genes induced only under cold room conditions could be divided into three major types: (1) genes associated with stress tolerance; (2) those that encode glycolytic and TCA cycle enzymes, and (3) those associated with protein synthesis machinery. A few of the genes induced only under field conditions appear to be related to light stress. Possible explanations for these differences are discussed in physiological context. Although many similarities exist in how plants respond during cold acclimation in the cold room and in the field environment, there are major differences suggesting caution should be taken in interpreting results based only on artificial, cold room conditions.

Comparative analysis of expressed sequence tags (ESTs) of ginseng leaf

The expressed sequence tags (ESTs) referenced in this report are the first transcriptomes in a leaf from a half-shade ginseng plant. A cDNA library was constructed from samples of the leaves of 4-year-old Panax ginseng plants, which were cultured in a field. The 2,896 P. ginseng cDNA clones represent 1,576 unique sequences, consisting of 1,167 singletons and 409 contig sequences. BLAST comparisons of the cDNAs in GenBank's non-redundant databases revealed that 2,579 of the 2,896 cDNAs (89.1%) exhibited a high degree of sequence homology to genes from other organisms. The majority of the identified transcripts were found to be genes related with energy, metabolism, subcellular localization, and protein synthesis and transport. The chlorophyll a/b-binding protein ESTs in the ginseng leaf samples manifested a substantially higher level of expression than was observed in other plant leaves. The ESTs involved in ginsenoside biosynthesis were also identified and discussed.

Abundant transcripts of malting barley identified by serial analysis of gene expression (SAGE)

Serial analysis of gene expression (SAGE) was applied to the major cereal crop barley (Hordeum vulgare) to characterize the transcriptional profile of grain during the malting process. Seven SAGE libraries were generated from seed at different time points during malting, in addition to one library from dry mature seed. A total of 155,206 LongSAGE tags, representing 41,909 unique sequences, was generated. This study reports an in-depth analysis of the most abundant transcripts from each of eight specific time points in a malting barley time course. The 100 most abundant tags from each library were analysed to identify the putative functional role of highly abundant transcripts. The largest functional groups included transcripts coding for stress response and cell defence, ribosomal proteins and storage proteins. The most abundant tag represented B22EL8, a barley metallothionein, which showed significant up-regulation across the malting time course. Considerable changes in the abundance profiles of some of the highly abundant tags occurred at 24 h post-steeping, indicating that it may be an important time point for gene expression changes associated with barley seed germination.

Characterizing the grape transcriptome. Analysis of expressed sequence tags from multiple Vitis species and development of a compendium of gene expression during berry development

We report the analysis and annotation of 146,075 expressed sequence tags from Vitis species. The majority of these sequences were derived from different cultivars of Vitis vinifera, comprising an estimated 25,746 unique contig and singleton sequences that survey transcription in various tissues and developmental stages and during biotic and abiotic stress. Putatively homologous proteins were identified for over 17,752 of the transcripts, with 1,962 transcripts further subdivided into one or more Gene Ontology categories. A simple structured vocabulary, with modules for plant genotype, plant development, and stress, was developed to describe the relationship between individual expressed sequence tags and cDNA libraries; the resulting vocabulary provides query terms to facilitate data mining within the context of a relational database. As a measure of the extent to which characterized metabolic pathways were encompassed by the data set, we searched for homologs of the enzymes leading from glycolysis, through the oxidative/nonoxidative pentose phosphate pathway, and into the general phenylpropanoid pathway. Homologs were identified for 65 of these 77 enzymes, with 86% of enzymatic steps represented by paralogous genes. Differentially expressed transcripts were identified by means of a stringent believability index cutoff of > or =98.4%. Correlation analysis and two-dimensional hierarchical clustering grouped these transcripts according to similarity of expression. In the broadest analysis, 665 differentially expressed transcripts were identified across 29 cDNA libraries, representing a range of developmental and stress conditions. The groupings revealed expected associations between plant developmental stages and tissue types, with the notable exception of abiotic stress treatments. A more focused analysis of flower and berry development identified 87 differentially expressed transcripts and provides the basis for a compendium that relates gene expression and annotation to previously characterized aspects of berry development and physiology. Comparison with published results for select genes, as well as correlation analysis between independent data sets, suggests that the inferred in silico patterns of expression are likely to be an accurate representation of transcript abundance for the conditions surveyed. Thus, the combined data set reveals the in silico expression patterns for hundreds of genes in V. vinifera, the majority of which have not been previously studied within this species.

Comparative analysis of expressed sequence tags from cold-acclimated and non-acclimated leaves of Rhododendron catawbiense Michx

An expressed sequence tag (EST) analysis approach was undertaken to identify major genes involved in cold acclimation of Rhododendron, a broad-leaf, woody evergreen species. Two cDNA libraries were constructed, one from winter-collected (cold-acclimated, CA; leaf freezing tolerance -53 degrees C) leaves, and the other from summer-collected (non-acclimated, NA; leaf freezing tolerance -7 degrees C) leaves of field-grown Rhododendron catawbiense plants. A total of 862 5'-end high-quality ESTs were generated by sequencing cDNA clones from the two libraries (423 from CA and 439 from NA library). Only about 6.3% of assembled unique transcripts were shared between the libraries, suggesting remarkable differences in gene expression between CA and NA leaves. Analysis of the relative frequency at which specific cDNAs were picked from each library indicated that four genes or gene families were highly abundant in the CA library including early light-induced proteins (ELIP), dehydrins/late embryogenesis abundant proteins (LEA), cytochrome P450, and beta-amylase. Similarly, seven genes or gene families were highly abundant in the NA library and included chlorophyll a/b-binding protein, NADH dehydrogenase subunit I, plastidic aldolase, and serine:glyoxylate aminotransferase, among others. Northern blot analyses for seven selected abundant genes confirmed their preferential expression in either CA or NA leaf tissues. Our results suggest that osmotic regulation, desiccation tolerance, photoinhibition tolerance, and photosynthesis adjustment are some of the key components of cold adaptation in Rhododendron.

Comparative analysis of expressed sequence tags from different organs of Vitis vinifera L

Expressed sequence tags (ESTs) are providing a valuable approach to sampling organism-expressed genomes, especially when studying large genomes such as those of many plants. We report on the comparison of 8,647 ESTs generated from six different grape (Vitis vinifera L.) organs: berry, root, leaf, bud, shoot and inflorescence. Clustering and assembly of these ESTs resulted in 4,203 unique sequences and revealed that at this level of EST sampling, each organ shares a low percentage of transcripts with the others. To define organ relationships based on EST counts, we calculated a distance matrix of pairwise correlation coefficients between the libraries which indicated bud, inflorescence and shoot as a group distinct from the other organs considered in this study. A putative function was identified for about 85% of the unique sequences. By assigning them to specific functional classes, we were able to highlight strong differences between organs in the metabolism, protein biosynthesis and photosynthesis categories. This grape EST collection has also proven to be a valuable source for the development of 'functional' simple sequence repeats (SSRs) markers: a total of 405 SSRs have been identified. EST sequences and annotation results have been organised in the IASMA-grape database, freely available at the address http://genomics.iasma.it.

cDNA microarray analysis of developing grape (Vitis vinifera cv. Shiraz) berry skin

Microarray analysis of Vitis vinifera cv. Shiraz developing berries has revealed the expression patterns of several categories of genes. Microarray slides were constructed from 4,608 PCR-amplified cDNA clones derived from a ripening grape berry cDNA library. The mRNA expression levels of the genes represented by these cDNAs were measured in flowers, week 2 post-flowering whole berries, week 5, week 8, week 10 (véraison, green berries), week 12 and week 13 berry skin. In addition, a comparison of RNA expression in pigmented and unpigmented berry skin at véraison (week 10) was undertaken. Image and statistical analysis revealed four sets of genes with distinctive and similar expression profiles over the course of berry development. The first set was composed of genes which had maximum RNA expression in flowers, followed by a steady decrease in expression. The most prominent group within this set were genes which have a role in photosynthesis. The second set of cDNAs was dominated by genes involved in flavonoid biosynthesis and had a peak of expression week 2 post-flowering. The data indicate co-ordinate regulation of flavonoid biosynthetic genes which code for the enzymes 4-coumarate-CoA ligase, chalcone synthase, chalcone isomerase, flavonone hydroxylase, anthocyanidin reductase and cytochrome b5. The third set of cDNAs exhibited maximum expression week 5 post-flowering, midway between flowering and véraison, a period of rapid berry growth. This set of cDNAs is dominated by genes which code for structural cell wall proteins. The fourth set of genes was dramatically up-regulated at véraison and remained up-regulated until 13 weeks post-flowering. This set of genes was composed of a diverse range of genes, a reflection of the complexity of ripening, most with no known function.

Grape berry biochemistry revisited upon proteomic analysis of the mesocarp

Major soluble proteins of grapevine ripe berries were extracted from six different cultivars including non vinifera, with trichloroacetic acid acetone and resolved in two-dimensional electrophoresis (2-DE) gels. About three hundred spots were detected on the 2-DE map after colloidal blue staining. From 2-DE map of cv. Gamay mesocarp, 67 proteins were identified (p > 0.95) using matrix-assisted laser desorption/ionization-mass spectrometry analysis. About 34%, 19%, and 13% of identified proteins play, respectively, a role in energy metabolism, defense, and stress response and primary metabolism. 2-DE analysis revealed considerable accumulation of dehydrin, invertase, and a putative transcription factor in the ripe fruit, in addition to pathogenesis-related proteins such as chitinase and thaumatin-like proteins previously documented as prevalent proteins in ripe berries. Actual translation of redundant transcripts of unclear function such as Grip31, Grip32, and Grip61 recently cloned in ripe grape berries was confirmed. The relative abundance of UDP-glucose pyrophosphorylase and vacuolar invertase strongly supported a key role of the apoplastic pathway of sugar loading during ripening. Comparative analysis shows that differences between cultivars were low, but different isoforms of alcohol dehydrogenase and of a transcription factor of hexose transporter were obvious in the six cultivars. Peptide mass fingerprinting suggests that the Adh isoforms would be Adh2/Adh6 or Adh2/Adh7 dimers and unambiguously shows that considerable deletion/insertion inside Adh7 are not cloning artifacts.

Quantitative trait locus analysis of fungal disease resistance factors on a molecular map of grapevine

A full-sibling F1 population comprising 153 individuals from the cross of 'Regent' x 'Lemberger' was employed to construct a genetic map based on 429 molecular markers. The newly-bred red grapevine variety 'Regent' has multiple field-resistance to fungal diseases inherited as polygenic traits, while 'Lemberger' is a traditional fungus-susceptible cultivar. The progeny segregate quantitatively for resistances to Plasmopara viticola and Uncinula necator, fungal pathogens that threaten viticulture in temperate areas. A double pseudo-testcross strategy was employed to construct the two parental maps under high statistical stringency for linkage to obtain a robust marker frame for subsequent quantitative trait locus (QTL) analysis. In total, 185 amplified fragment length polymorphism, 137 random amplified polymorphic DNA, 85 single sequence repeat and 22 sequence characterized amplified region or cleaved amplified polymorphic sequence markers were mapped. The maps were aligned by co-dominant or doubly heterozygous dominant anchor markers. Twelve pairs of homologous linkage groups could be integrated into consensus linkage groups. Resistance phenotypes and segregating characteristics were scored as quantitative traits in three or four growing seasons. Interval mapping reproducibly localized genetic factors that correlated with fungal disease resistances to specific regions on three linkage groups of the maternal 'Regent' map. A QTL for resistance to Uncinula necator was identified on linkage group 16, and QTLs for endurance to Plasmopara viticola on linkage groups 9 and 10 of 'Regent'. Additional QTLs for the onset of berry ripening ("veraison"), berry size and axillary shoot growth were identified. Berry color segregated as a simple trait in this cross of two red varieties and was mapped as a morphological marker. Six markers derived from functional genes could be localized. This dissection of polygenic fungus disease resistance in grapevine allows the development of marker-assisted selection for breeding, the characterization of genetic resources and the isolation of the corresponding genes.

Sequence Tag Catalogs of Dust Mite-Expressed Genomes

Dust mites are a major source of indoor allergens. They contain a large number of components that react with immunoglobulin (Ig) E in individuals with allergies and are capable of inducing sensitization, and allergic respiratory and cutaneous diseases. With a significant proportion of the population affected in some way by mite allergies, it is essential that we improve our understanding of these organisms so that control strategies could be defined and its allergens better understood. Thus, we have initiated a project using the expressed sequence tagging (EST) strategy to study the major species of dust mites associated with allergic diseases, in particular, the American house dust mite, Dermatophagoides farinae, as well as Blomia tropicalis, the most prevalent mite in domestic tropical dwellings. The work has recently been expanded to include 'storage' mites such as Tyrophagus putrescentiae, Acarus siro, Lepidoglyphus destructor, Glycyphagus domesticus, Suidasia medanensis, and Aleuroglyphus ovatus. More than 50% of the initial 3000 ESTs from the D. farinae and B. tropicalis dust mites showed significant matches to known genes and were categorized into eight functional groups (such as proteins involved in metabolism, gene expression, protein synthesis, cell signaling, etc.). Of specific interest, however, were the homologs to known mite allergens, in addition to a number of sequences bearing significant homology to allergens from non-mite sources previously not known to exist in mites. The availability of these allergen sequences has facilitated their expression and subsequent characterization in our laboratory in terms of their IgE-binding reactivity. The wealth of sequence information, generated via the EST project, has also facilitated the identification of polymorphic forms of allergens, the investigation of differential gene expression under various environmental conditions via DNA microarrays, as well as the analysis of protein level expression profiling via the proteomics approach. Additionally, ESTs have also ameliorated the understanding of the phylogenetic relationships between mites, and enabled the isolation of gene products crucial for life processes so that mite control strategies can be more effectively devised. Taken together, the utilization of the EST strategy has opened up numerous new avenues by which the allergist can engage more effectively in the study of dust mites with the ultimate aim of developing appropriate treatment regimens for mite-induced allergy.

Genes associated with the end of dormancy in grapes

A grape bud EST library was constructed and 4270 ESTs sequenced. The library clones were arrayed for the purpose of investigating the level of gene expression over time, particularly leading up to the buds' release from dormancy. The arrays were hybridized with P(33)-labeled probes produced from samples of buds collected at weekly intervals. These probes covered the time from 9 weeks prior to bud burst until just after the emergence of the shoots. Expression patterns from these genes have been examined. It was found that 74% of the genes in the data set were homologous to known proteins. Genes were then assigned to functional categories according to their primary BLAST match. Of these 13% were involved with photosynthesis, 13% with disease resistance and defense, 5% energy, 12% metabolism, 20% protein production and processing, 25% cell structure and plant growth and the remaining 12% were unclassified The expression pattern of a selection of "candidate" genes retrieved from literature previously reporting an association with dormancy changes was assessed. On closer examination most of these genes relate to the oxidative processes and stress responses within the cell. The results of this study show that even in the dormant state, gene expression in the buds is high.

Molecular linkage maps of Vitis vinifera L. and Vitis riparia Mchx

Two linkage maps for grape (Vitis spp.) have been developed based on 81 F(1) plants derived from an interspecific cross between the wine cultivar Moscato bianco (Vitis vinifera L.) and a Vitis riparia Mchx. accession, a donor of pathogen resistance traits. The double pseudotest-cross mapping strategy was applied using three types of molecular markers. The efficiency of SSRs to anchor homologous linkage groups from different Vitis maps and the usefulness of AFLPs in saturating molecular linkage maps were evaluated. Moreover, the SSCP technique was developed based on sequence information in public databases concerning genes involved in flavonoid and stilbene biosynthesis. For the maternal genetic map a total of 338 markers were assembled in 20 linkage groups covering 1,639 cM, whereas 429 loci defined the 19 linkage groups of the paternal map which covers 1,518 cM. The identification of 14 linkage groups common to both maps was possible based on 21 SSR and 19 AFLP loci. The position of SSR loci in the maps presented here was consistent with other published mapping experiments in Vitis.

Identification of genes possibly related to storage root induction in sweet potato

To identify genes related to initiation of storage root development in sweet potato, a cDNA library was constructed with early stage storage roots (0.3-1 cm in diameter). Single-pass sequences of the 5' ends of 2859 sweet potato cDNA clones were assembled into 483 clusters and 442 singletons. Comparison of sweet potato expressed sequence tags (ESTs) to nodulation/tumorigenesis-related sequence databases (nodule-, tumor-, potato tuber- and development-related sequences) revealed that homologs of 39 sweet potato EST sequences potentially involved in gene regulation, signal transduction and development were present in at least one of the nodulation/tumorigenesis-related sequence databases. Northern blot analyses of these 39 sequences identified 22 differentially expressed genes in early stage storage root and fibrous root. These differentially expressed genes will be potential candidates for research to elucidate the molecular processes related to sweet potato storage root induction.

A survey of gene expression and diversity in the venom glands of the pitviper snake Bothrops insularis through the generation of expressed sequence tags (ESTs)

In order to produce a global panorama of the transcriptional activity of snake venom glands and to correlate with its venom composition, we constructed a DNA complementary to RNA library from the venom glands of the Viperidae snake Bothrops insularis for the generation of expressed sequence tags (ESTs). Sequences from 610 independent clones were grouped in 297 clusters, revealing the putative identification of 210 distinct gene products. Toxin sequences correspond to 56% of all transcripts (85 clusters), being the metalloproteinases (23%) and the bradykinin-potentiating peptides (11%) the major components. This approach revealed a new highly expressed toxin similar to vascular endothelial growth factor, which was recently reported (J. Biol. Chem. 276 (2001) 39836). Among the 125 clusters matching cellular proteins, the major part represents molecules involved in gene and protein expression, notably in disulfide bond assembly, reflecting a high specialization of this tissue for toxin synthesis. An unusual representation of retrotransposon-like sequences was also found and could be related to the occurrence and diversity of many paralogous forms of toxins in the venom gland. Our B. insularis dbEST allowed the identification of the most common classes of toxins present in Viperidae venoms, which parallels the complex hemorrhagic effects evoked by the venom on the prey. In addition, it provides the first comprehensive set of reptilian gene sequences described so far.

Genome-wide identification of nodule-specific transcripts in the model legume Medicago truncatula

The Medicago truncatula expressed sequence tag (EST) database (Gene Index) contains over 140,000 sequences from 30 cDNA libraries. This resource offers the possibility of identifying previously uncharacterized genes and assessing the frequency and tissue specificity of their expression in silico. Because M. truncatula forms symbiotic root nodules, unlike Arabidopsis, this is a particularly important approach in investigating genes specific to nodule development and function in legumes. Our analyses have revealed 340 putative gene products, or tentative consensus sequences (TCs), expressed solely in root nodules. These TCs were represented by two to 379 ESTs. Of these TCs, 3% appear to encode novel proteins, 57% encode proteins with a weak similarity to the GenBank accessions, and 40% encode proteins with strong similarity to the known proteins. Nodule-specific TCs were grouped into nine categories based on the predicted function of their protein products. Besides previously characterized nodulins, other examples of highly abundant nodule-specific transcripts include plantacyanin, agglutinin, embryo-specific protein, and purine permease. Six nodule-specific TCs encode calmodulin-like proteins that possess a unique cleavable transit sequence potentially targeting the protein into the peribacteroid space. Surprisingly, 114 nodule-specific TCs encode small Cys cluster proteins with a cleavable transit peptide. To determine the validity of the in silico analysis, expression of 91 putative nodule-specific TCs was analyzed by macroarray and RNA-blot hybridizations. Nodule-enhanced expression was confirmed experimentally for the TCs composed of five or more ESTs, whereas the results for those TCs containing fewer ESTs were variable.