Sequencing and analysis of common bean ESTs. Building a foundation for functional genomics

Transcript profiling of common bean (Phaseolus vulgaris L.) using the GeneChip Soybean Genome Array: optimizing analysis by masking biased probes

BACKGROUND

Common bean (Phaseolus vulgaris L.) and soybean (Glycine max) both belong to the Phaseoleae tribe and share significant coding sequence homology. This suggests that the GeneChip(R) Soybean Genome Array (soybean GeneChip) may be used for gene expression studies using common bean.

RESULTS

To evaluate the utility of the soybean GeneChip for transcript profiling of common bean, we hybridized cRNAs purified from nodule, leaf, and root of common bean and soybean in triplicate to the soybean GeneChip. Initial data analysis showed a decreased sensitivity and accuracy of measuring differential gene expression in common bean cross-species hybridization (CSH) GeneChip data compared to that of soybean. We employed a method that masked putative probes targeting inter-species variable (ISV) regions between common bean and soybean. A masking signal intensity threshold was selected that optimized both sensitivity and accuracy of measuring differential gene expression. After masking for ISV regions, the number of differentially-expressed genes identified in common bean was increased by 2.8-fold reflecting increased sensitivity. Quantitative RT-PCR (qRT-PCR) analysis of 20 randomly selected genes and purine-ureide pathway genes demonstrated an increased accuracy of measuring differential gene expression after masking for ISV regions. We also evaluated masked probe frequency per probe set to gain insight into the sequence divergence pattern between common bean and soybean. The sequence divergence pattern analysis suggested that the genes for basic cellular functions and metabolism were highly conserved between soybean and common bean. Additionally, our results show that some classes of genes, particularly those associated with environmental adaptation, are highly divergent.

CONCLUSIONS

The soybean GeneChip is a suitable cross-species platform for transcript profiling in common bean when used in combination with the masking protocol described. In addition to transcript profiling, CSH of the GeneChip in combination with masking probes in the ISV regions can be used for comparative ecological and/or evolutionary genomics studies.

Cloning and genetic diversity analysis of a new P5CS gene from common bean (Phaseolus vulgaris L.)

Delta(1)-pyrroline-5-carboxylate synthetase (P5CS) is the rate-limiting enzyme involved in the biosynthesis of proline in plants. By the 3' rapid amplification of cDNA ends (3'-RACE) approach, a 2,246-bp cDNA sequence was obtained from common bean (Phaseolus vulgaris L.), denominated PvP5CS2 differing from another P5CS gene that we cloned previously from common bean (PvP5CS). The predicted amino acid sequence of PvP5CS2 has an overall 93.2% identity GmP5CS (Glycine max L. P5CS). However, PvP5CS2 shows only 83.7% identity in amino acid sequence to PvP5CS, suggesting PvP5CS2 represents a homolog of the soybean P5CS gene. Abundant indel (insertion and deletion events) and SNP (single nucleotide polymorphisms) were found in the cloned PvP5CS2 genome sequence when comparing 24 cultivated and 3 wild common bean accessions and these in turn reflected aspects of common bean evolution. Sequence alignment showed that genotypes from the same gene pool had similar nucleotide variation, while genotypes from different gene pools had distinctly different nucleotide variation for PvP5CS2. Furthermore, diversity along the gene sequence was not evenly distributed, being low in the glutamic-g-semialdehyde dehydrogenase catalyzing region, moderate in the Glu-5-kinase catalyzing region and high in the intervening region. Neutrality tests showed that PvP5CS2 was a conserved gene undergoing negative selection. A new marker (Pv97) was developed for genetic mapping of PvP5CS2 based on an indel between DOR364 and G19833 sequences and the gene was located between markers Bng126 and BMd045 on chromosome b01. The relationship of PvP5CS2 and a previously cloned pyrroline-5-carboxylate synthetase gene as well as the implications of this work on selecting for drought tolerance in common bean are discussed.

Synteny mapping between common bean and soybean reveals extensive blocks of shared loci

BACKGROUND

Understanding syntentic relationship between two species is critical to assessing the potential for comparative genomic analysis. Common bean (Phaseolus vulgaris L.) and soybean (Glycine max L.), the two most important members of the Phaseoleae legumes, appear to have a diploid and polyploidy recent past, respectively. Determining the syntentic relationship between these two species will allow researchers to leverage not only genomic resources but also genetic data for important agronomic traits to improve both of these species.

RESULTS

Genetically-positioned transcript loci of common bean were mapped relative to the recent soybean 1.01 pseudochromosome assembly. In nearly every case, each common bean locus mapped to two loci in soybean, a result consistent with the duplicate polyploidy history of soybean. Blocks of synteny averaging 32 cM in common bean and 4.9 Mb in soybean were observed for all 11 common bean linkage groups, and these blocks mapped to all 20 soybean pseudochromosomes. The median physical-to-genetic distance ratio in common bean (based on soybean physical distances) was approximately 120 kb/cM. approximately 15,000 common bean sequences (primarily EST contigs and EST singletons) were electronically positioned onto the common bean map using the shared syntentic blocks as references points.

CONCLUSION

The collected evidence from this mapping strongly supports the duplicate history of soybean. It further provides evidence that the soybean genome was fractionated and reassembled at some point following the duplication event. These well mapped syntentic relationships between common bean and soybean will enable researchers to target specific genomic regions to discover genes or loci that affect phenotypic expression in both species.

Identification of candidates for cyclotide biosynthesis and cyclisation by expressed sequence tag analysis of Oldenlandia affinis

Background

Cyclotides are a family of circular peptides that exhibit a range of biological activities, including anti-bacterial, cytotoxic, anti-HIV activities, and are proposed to function in plant defence. Their high stability has motivated their development as scaffolds for the stabilisation of peptide drugs. Oldenlandia affinis is a member of the Rubiaceae (coffee) family from which 18 cyclotides have been sequenced to date, but the details of their processing from precursor proteins have only begun to be elucidated. To increase the speed at which genes involved in cyclotide biosynthesis and processing are being discovered, an expressed sequence tag (EST) project was initiated to survey the transcript profile of O. affinis and to propose some future directions of research on in vivo protein cyclisation.

Results

Using flow cytometry the holoploid genome size (1C-value) of O. affinis was estimated to be 4,210 - 4,284 Mbp, one of the largest genomes of the Rubiaceae family. High-quality ESTs were identified, 1,117 in total, from leaf cDNAs and assembled into 502 contigs, comprising 202 consensus sequences and 300 singletons. ESTs encoding the cyclotide precursors for kalata B1 (Oak1) and kalata B2 (Oak4) were among the 20 most abundant ESTs. In total, 31 ESTs encoded cyclotide precursors, representing a distinct commitment of 2.8% of the O. affinis transcriptome to cyclotide biosynthesis. The high expression levels of cyclotide precursor transcripts are consistent with the abundance of mature cyclic peptides in O. affinis. A new cyclotide precursor named Oak5 was isolated and represents the first cDNA for the bracelet class of cyclotides in O. affinis. Clones encoding enzymes potentially involved in processing cyclotides were also identified and include enzymes involved in oxidative folding and proteolytic processing.

Conclusion

The EST library generated in this study provides a valuable resource for the study of the cyclisation of plant peptides. Further analysis of the candidates for cyclotide processing discovered in this work will increase our understanding and aid in reconstructing cyclotide production using transgenic systems and will benefit their development in pharmaceutical applications and insect-resistant crop plants.

Biochemical and molecular characterization of PvPAP3, a novel purple acid phosphatase isolated from common bean enhancing extracellular ATP utilization

Purple acid phosphatases (PAPs) play diverse physiological roles in plants. In this study, we purified a novel PAP, PvPAP3, from the roots of common bean (Phaseolus vulgaris) grown under phosphate (Pi) starvation. PvPAP3 was identified as a 34-kD monomer acting on the specific substrate, ATP, with a broad pH range and a high heat stability. The activity of PvPAP3 was insensitive to tartrate, indicating that PvPAP3 is a PAP-like protein. Amino acid sequence alignment and phylogenetic analysis suggest that PvPAP3 belongs to the group of plant PAPs with low molecular mass. Transient expression of 35S:PvPAP3-green fluorescent protein in onion (Allium cepa) epidermal cells verified that it might anchor on plasma membrane and be secreted into apoplast. Pi starvation led to induction of PvPAP3 expression in both leaves and roots of common bean, and expression of PvPAP3 was strictly dependent on phosphorus (P) availability and duration of Pi starvation. Furthermore, induction of PvPAP3 expression was more rapid and higher in a P-efficient genotype, G19833, than in a P-inefficient genotype, DOR364, suggesting possible roles of PvPAP3 in P efficiency in bean. In vivo analysis using a transgenic hairy root system of common bean showed that both growth and P uptake of bean hairy roots from the PvPAP3 overexpression transgenic lines were significantly enhanced when ATP was supplied as the sole external P source. Taken together, our results suggest that PvPAP3 is a novel PAP that might function in the adaptation of common bean to P deficiency, possibly through enhancing utilization of extracellular ATP as a P source.

Single strand conformation polymorphism based SNP and Indel markers for genetic mapping and synteny analysis of common bean (Phaseolus vulgaris L.)

BACKGROUND

Expressed sequence tags (ESTs) are an important source of gene-based markers such as those based on insertion-deletions (Indels) or single-nucleotide polymorphisms (SNPs). Several gel based methods have been reported for the detection of sequence variants, however they have not been widely exploited in common bean, an important legume crop of the developing world. The objectives of this project were to develop and map EST based markers using analysis of single strand conformation polymorphisms (SSCPs), to create a transcript map for common bean and to compare synteny of the common bean map with sequenced chromosomes of other legumes.

RESULTS

A set of 418 EST based amplicons were evaluated for parental polymorphisms using the SSCP technique and 26% of these presented a clear conformational or size polymorphism between Andean and Mesoamerican genotypes. The amplicon based markers were then used for genetic mapping with segregation analysis performed in the DOR364 x G19833 recombinant inbred line (RIL) population. A total of 118 new marker loci were placed into an integrated molecular map for common bean consisting of 288 markers. Of these, 218 were used for synteny analysis and 186 presented homology with segments of the soybean genome with an e-value lower than 7 x 10-12. The synteny analysis with soybean showed a mosaic pattern of syntenic blocks with most segments of any one common bean linkage group associated with two soybean chromosomes. The analysis with Medicago truncatula and Lotus japonicus presented fewer syntenic regions consistent with the more distant phylogenetic relationship between the galegoid and phaseoloid legumes.

CONCLUSION

The SSCP technique is a useful and inexpensive alternative to other SNP or Indel detection techniques for saturating the common bean genetic map with functional markers that may be useful in marker assisted selection. In addition, the genetic markers based on ESTs allowed the construction of a transcript map and given their high conservation between species allowed synteny comparisons to be made to sequenced genomes. This synteny analysis may support positional cloning of target genes in common bean through the use of genomic information from these other legumes.

A nomadic subtelomeric disease resistance gene cluster in common bean

The B4 resistance (R) gene cluster is one of the largest clusters known in common bean (Phaseolus vulgaris [Pv]). It is located in a peculiar genomic environment in the subtelomeric region of the short arm of chromosome 4, adjacent to two heterochromatic blocks (knobs). We sequenced 650 kb spanning this locus and annotated 97 genes, 26 of which correspond to Coiled-Coil-Nucleotide-Binding-Site-Leucine-Rich-Repeat (CNL). Conserved microsynteny was observed between the Pv B4 locus and corresponding regions of Medicago truncatula and Lotus japonicus in chromosomes Mt6 and Lj2, respectively. The notable exception was the CNL sequences, which were completely absent in these regions. The origin of the Pv B4-CNL sequences was investigated through phylogenetic analysis, which reveals that, in the Pv genome, paralogous CNL genes are shared among nonhomologous chromosomes (4 and 11). Together, our results suggest that Pv B4-CNL was derived from CNL sequences from another cluster, the Co-2 cluster, through an ectopic recombination event. Integration of the soybean (Glycine max) genome data enables us to date more precisely this event and also to infer that a single CNL moved from the Co-2 to the B4 cluster. Moreover, we identified a new 528-bp satellite repeat, referred to as khipu, specific to the Phaseolus genus, present both between B4-CNL sequences and in the two knobs identified at the B4 R gene cluster. The khipu repeat is present on most chromosomal termini, indicating the existence of frequent ectopic recombination events in Pv subtelomeric regions. Our results highlight the importance of ectopic recombination in R gene evolution.

Global changes in the transcript and metabolic profiles during symbiotic nitrogen fixation in phosphorus-stressed common bean plants

Phosphorus (P) deficiency is widespread in regions where the common bean (Phaseolus vulgaris), the most important legume for human consumption, is produced, and it is perhaps the factor that most limits nitrogen fixation. Global gene expression and metabolome approaches were used to investigate the responses of nodules from common bean plants inoculated with Rhizobium tropici CIAT899 grown under P-deficient and P-sufficient conditions. P-deficient inoculated plants showed drastic reduction in nodulation and nitrogenase activity as determined by acetylene reduction assay. Nodule transcript profiling was performed through hybridization of nylon filter arrays spotted with cDNAs, approximately 4,000 unigene set, from the nodule and P-deficient root library. A total of 459 genes, representing different biological processes according to updated annotation using the UniProt Knowledgebase database, showed significant differential expression in response to P: 59% of these were induced in P-deficient nodules. The expression platform for transcription factor genes based in quantitative reverse transcriptase-polymerase chain reaction revealed that 37 transcription factor genes were differentially expressed in P-deficient nodules and only one gene was repressed. Data from nontargeted metabolic profiles indicated that amino acids and other nitrogen metabolites were decreased, while organic and polyhydroxy acids were accumulated, in P-deficient nodules. Bioinformatics analyses using MapMan and PathExpress software tools, customized to common bean, were utilized for the analysis of global changes in gene expression that affected overall metabolism. Glycolysis and glycerolipid metabolism, and starch and Suc metabolism, were identified among the pathways significantly induced or repressed in P-deficient nodules, respectively.

Comparative analyses of genotype dependent expressed sequence tags and stress-responsive transcriptome of chickpea wilt illustrate predicted and unexpected genes and novel regulators of plant immunity

BACKGROUND

The ultimate phenome of any organism is modulated by regulated transcription of many genes. Characterization of genetic makeup is thus crucial for understanding the molecular basis of phenotypic diversity, evolution and response to intra- and extra-cellular stimuli. Chickpea is the world's third most important food legume grown in over 40 countries representing all the continents. Despite its importance in plant evolution, role in human nutrition and stress adaptation, very little ESTs and differential transcriptome data is available, let alone genotype-specific gene signatures. Present study focuses on Fusarium wilt responsive gene expression in chickpea.

RESULTS

We report 6272 gene sequences of immune-response pathway that would provide genotype-dependent spatial information on the presence and relative abundance of each gene. The sequence assembly led to the identification of a CaUnigene set of 2013 transcripts comprising of 973 contigs and 1040 singletons, two-third of which represent new chickpea genes hitherto undiscovered. We identified 209 gene families and 262 genotype-specific SNPs. Further, several novel transcription regulators were identified indicating their possible role in immune response. The transcriptomic analysis revealed 649 non-cannonical genes besides many unexpected candidates with known biochemical functions, which have never been associated with pathostress-responsive transcriptome.

CONCLUSION

Our study establishes a comprehensive catalogue of the immune-responsive root transcriptome with insight into their identity and function. The development, detailed analysis of CaEST datasets and global gene expression by microarray provide new insight into the commonality and diversity of organ-specific immune-responsive transcript signatures and their regulated expression shaping the species specificity at genotype level. This is the first report on differential transcriptome of an unsequenced genome during vascular wilt.

Development and diversity of Andean-derived, gene-based microsatellites for common bean (Phaseolus vulgaris L.)

BACKGROUND

Gene-based (genic) microsatellites are a useful tool for plant genetics and simple sequence repeat loci can often be found in coding regions of the genome. While EST sequencing can be used to discover genic microsatellites, direct screening of cDNA libraries for repeat motifs can save on overall sequencing costs. The objective of this research was to screen a large cDNA library from and Andean common bean genotype for six di-nucleotide and tri-nucleotide repeat motifs through a filter hybridization approach and to develop microsatellite markers from positive clones.

RESULTS

Robotics were used for high-throughput colony picking and to create a high-density filter of 18,432 double spotted cDNA clones which was followed by hybridization with repeat motif containing probes based on GA, CA, AAT, CAG, CAA and ACG repeats. A total of 1203 positive clones were identified by their addresses and sequenced from 5' ends and if required from 3' ends to confirm repeat motif and length. Out of 886 high quality sequences, 497 had complete microsatellite loci that were not truncated by the sequencing reaction and of these tri-nucleotide repeats were more common than di-nucleotide repeats. Different motifs were found in different frequencies in the 5' and 3' ends of the cDNAs. In a microsatellite development program, primers were designed for 248 SSR loci which were tested on a panel of 18 common bean genotypes to determine their potential as genetic markers finding higher average polymorphism information content for di-nucleotide repeat markers (0.3544) than for tri-nucleotide repeat markers (0.1536).

CONCLUSION

The present study provides a set of validated gene-based markers for common bean that are derived from G19833, an Andean landrace that is an important source of disease and abiotic stress tolerance which has been used for physical map development and as a mapping parent. Gene-based markers appear to be very efficient at separating divergent wild and cultivated accessions as well as Andean and Mesoamerican gene pools and therefore will be useful for diversity analyses and for comparative and transcript mapping in common bean.

Generation of Phaseolus vulgaris ESTs and investigation of their regulation upon Uromyces appendiculatus infection

BACKGROUND

Phaseolus vulgaris (common bean) is the second most important legume crop in the world after soybean. Consequently, yield losses due to fungal infection, like Uromyces appendiculatus (bean rust), have strong consequences. Several resistant genes were identified that confer resistance to bean rust infection. However, the downstream genes and mechanisms involved in bean resistance to infection are poorly characterized.

RESULTS

A subtractive bean cDNA library composed of 10,581 unisequences was constructed and enriched in sequences regulated by either bean rust race 41, a virulent strain, or race 49, an avirulent strain on cultivar Early Gallatin carrying the resistance gene Ur-4. The construction of this library allowed the identification of 6,202 new bean ESTs, significantly adding to the available sequences for this plant. Regulation of selected bean genes in response to bean rust infection was confirmed by qRT-PCR. Plant gene expression was similar for both race 41 and 49 during the first 48 hours of the infection process but varied significantly at the later time points (72-96 hours after inoculation) mainly due to the presence of the Avr4 gene in the race 49 leading to a hypersensitive response in the bean plants. A biphasic pattern of gene expression was observed for several genes regulated in response to fungal infection.

CONCLUSION

The enrichment of the public database with over 6,000 bean ESTs significantly adds to the genomic resources available for this important crop plant. The analysis of these genes in response to bean rust infection provides a foundation for further studies of the mechanism of fungal disease resistance. The expression pattern of 90 bean genes upon rust infection shares several features with other legumes infected by biotrophic fungi. This finding suggests that the P. vulgaris-U. appendiculatus pathosystem could serve as a model to explore legume-rust interaction.

Development of chickpea EST-SSR markers and analysis of allelic variation across related species

Despite chickpea being the third important grain legume, there is a limited availability of genomic resources, especially of the expressed sequence tag (EST)-based markers. In this study, we generated 822 chickpea ESTs from immature seeds as well as exploited 1,309 ESTs from the chickpea database, thus utilizing a total of 2,131 EST sequences for development of functional EST-SSR markers. Two hundred and forty-six simple sequence repeat (SSR) motifs were identified from which 183 primer pairs were designed and 60 validated as functional markers. Genetic diversity analysis across 30 chickpea accessions revealed ten markers to be polymorphic producing a total of 29 alleles and an observed heterozygosity average of 0.16 thereby exhibiting low levels of intra-specific polymorphism. However, the markers exhibited high cross-species transferability ranging from 68.3 to 96.6% across the six annual Cicer species and from 29.4 to 61.7% across the seven legume genera. Sequence analysis of size variant amplicons from various species revealed that size polymorphism was due to multiple events such as copy number variation, point mutations and insertions/deletions in the microsatellite repeat as well as in the flanking regions. Interestingly, a wide prevalence of crossability-group-specific sequence variations were observed among Cicer species that were phylogenetically informative. The neighbor joining dendrogram clearly separated the chickpea cultivars from the wild Cicer and validated the proximity of C. judaicum with C. pinnatifidum. Hence, this study for the first time provides an insight into the distribution of SSRs in the chickpea transcribed regions and also demonstrates the development and utilization of genic-SSRs. In addition to proving their suitability for genetic diversity analysis, their high rates of transferability also proved their potential for comparative genomic studies and for following gene introgressions and evolution in wild species, which constitute the valuable secondary genepool in chickpea.

Transcript and proteomic analysis of developing white lupin (Lupinus albus L.) roots

BACKGROUND

White lupin (Lupinus albus L.) roots efficiently take up and accumulate (heavy) metals, adapt to phosphate deficiency by forming cluster roots, and secrete antimicrobial prenylated isoflavones during development. Genomic and proteomic approaches were applied to identify candidate genes and proteins involved in antimicrobial defense and (heavy) metal uptake and translocation.

RESULTS

A cDNA library was constructed from roots of white lupin seedlings. Eight thousand clones were randomly sequenced and assembled into 2,455 unigenes, which were annotated based on homologous matches in the NCBInr protein database. A reference map of developing white lupin root proteins was established through 2-D gel electrophoresis and peptide mass fingerprinting. High quality peptide mass spectra were obtained for 170 proteins. Microsomal membrane proteins were separated by 1-D gel electrophoresis and identified by LC-MS/MS. A total of 74 proteins were putatively identified by the peptide mass fingerprinting and the LC-MS/MS methods. Genomic and proteomic analyses identified candidate genes and proteins encoding metal binding and/or transport proteins, transcription factors, ABC transporters and phenylpropanoid biosynthetic enzymes.

CONCLUSION

The combined EST and protein datasets will facilitate the understanding of white lupin's response to biotic and abiotic stresses and its utility for phytoremediation. The root ESTs provided 82 perfect simple sequence repeat (SSR) markers with potential utility in breeding white lupin for enhanced agronomic traits.

Essential role of MYB transcription factor: PvPHR1 and microRNA: PvmiR399 in phosphorus-deficiency signalling in common bean roots

Phosphorus (P), an essential element for plants, is one of the most limiting nutrients for plant growth. A few transcription factor (TF) genes involved in P-starvation signalling have been characterized for Arabidopsis thaliana and rice. Crop production of common bean (Phaseolus vulgaris L.), the most important legume for human consumption, is often limited by low P in the soil. Despite its agronomic importance, nothing is known about transcriptional regulation in P-deficient bean plants. We functionally characterized the P-deficiency-induced MYB TF TC3604 (Dana Farber Cancer Institute, Common Bean Gene Index v.2.0), ortholog to AtPHR1 (PvPHR1). For its study, we applied RNAi technology in bean composite plants. PvPHR1 is a positive regulator of genes implicated in P transport, remobilization and homeostasis. Although there are no reports on the regulatory roles of microRNAs (miRNA) in bean, we demonstrated that PvmiR399 is an essential component of the PvPHR1 signalling pathway. The analysis of DICER-like1 (PvDCL1) silenced bean composite plants suppressed for accumulation of PvmiR399 and other miRNAs suggested that miR399 is a negative regulator of the ubiquitin E2 conjugase: PvPHO2 expression. Our results set the basis for understanding the signalling for P-starvation responses in common bean and may contribute to crop improvement.

Transcriptional regulation and signaling in phosphorus starvation: what about legumes?

The availability of soil phosphorus (P), an essential element, is one of the most important requirements for plant growth and crop production. The morphological and physiological adaptations evolved by plants to cope with P starvation have been well characterized. Several P deficiency plant responses are regulated at the transcriptional level. Microarray analysis has generated valuable information on global gene expression in Arabidopsis thaliana grown under P-stress. Despite the identification of P responsive genes, little is known about the regulation of gene expression changes. Four transcription factors, PHR1, WRKY75, ZAT6 and BHLH32, involved in P starvation signaling have been characterized in Arabidopsis, and signaling pathways are deciphered. This review analyzes the current knowledge of transcriptional regulation of P starvation responses in Arabidopsis vis-à-vis legumes such as lupine, common bean and Medicago truncatula. The knowledge on regulatory and signaling mechanisms involved in P acquisition and use in legumes will be useful for improvement of these crops, which account for a large proportion of the world's crop production, providing good nutritional quality feed and food.

Expressed sequence tags of Trichinella spiralis muscle stage larvae

In order to obtain greater insight into the relevant genomic expression patterns of Trichinella spiralis, 992 expressed sequence tags (ESTs) were collected from a cDNA library of T. spiralis muscle stage larvae and assembled into 60 clusters and 385 singletons. Of them, 445 (44.7%) ESTs were annotated to their homologous genes, and small fractions were matched to known genes of nematodes. The annotated ESTs were classified into 25 eukaryotic orthologous groups (KOG). Cytochrome C oxidase (34 clones) was found to be most frequent species.

Improvement of drought tolerance and grain yield in common bean by overexpressing trehalose-6-phosphate synthase in rhizobia

Improving stress tolerance and yield in crops are major goals for agriculture. Here, we show a new strategy to increase drought tolerance and yield in legumes by overexpressing trehalose-6-phosphate synthase in the symbiotic bacterium Rhizobium etli. Phaseolus vulgaris (common beans) plants inoculated with R. etli overexpressing trehalose-6-phosphate synthase gene had more nodules with increased nitrogenase activity and higher biomass compared with plants inoculated with wild-type R. etli. In contrast, plants inoculated with an R. etli mutant in trehalose-6-phosphate synthase gene had fewer nodules and less nitrogenase activity and biomass. Three-week-old plants subjected to drought stress fully recovered whereas plants inoculated with a wild-type or mutant strain wilted and died. The yield of bean plants inoculated with R. etli overexpressing trehalose-6-phosphate synthase gene and grown with constant irrigation increased more than 50%. Macroarray analysis of 7,200 expressed sequence tags from nodules of plants inoculated with the strain overexpressing trehalose-6-phosphate synthase gene revealed upregulation of genes involved in stress tolerance and carbon and nitrogen metabolism, suggesting a signaling mechanism for trehalose. Thus, trehalose metabolism in rhizobia is key for signaling plant growth, yield, and adaptation to abiotic stress, and its manipulation has a major agronomical impact on leguminous plants.

Host genes involved in nodulation preference in common bean (Phaseolus vulgaris)-rhizobium etli symbiosis revealed by suppressive subtractive hybridization

Common bean cultivars are nodulated preferentially by Rhizobium etli lineages from the same center of host diversification. Nodulation was found to be earlier and numerous in bean plants inoculated with the cognate strain. We predicted that analysis of transcripts at early stages of the interaction between host and rhizobium would identify plant genes that are most likely to be involved in this preferential nodulation. Therefore, we applied a suppressive subtractive hybridization approach in which cDNA from a Mesoamerican cultivar inoculated with either the more- or less-efficient strain of R. etli was used as the driver and the tester, respectively. Forty-one independent tentative consensus sequences (TCs) were obtained and classified into different functional categories. Of 11 selected TCs, 9 were confirmed by quantitative reverse-transcriptase polymerase chain reaction. Two genes show high homology to previously characterized plant receptors. Two other upregulated genes encode for Rab11, a member of the small GTP-binding protein family, and HAP5, a subunit of the heterotrimeric CCAAT-transcription factor. Interestingly, one of the TCs encodes for an isoflavone reductase, which may lead to earlier Nod factor production by specific strains of rhizobia. The transcript abundance of selected cDNAs also was found to be higher in mature nodules of the more efficient interaction. Small or no differences were observed when an Andean bean cultivar was inoculated with a cognate strain, suggesting involvement of these genes in the strain-specific response. The potential role of these genes in the early preferential symbiotic interaction is discussed.

Ascorbate and homoglutathione metabolism in common bean nodules under stress conditions and during natural senescence

Ascorbate and glutathione are major antioxidants and redox buffers in plant cells but also play key functions in growth, development, and stress responses. We have studied the regulation of ascorbate and homoglutathione biosynthesis in common bean (Phaseolus vulgaris) nodules under stress conditions and during aging. The expression of five genes of the major ascorbate biosynthetic pathway was analyzed in nodules, and evidence was found that L-galactono-1,4-lactone dehydrogenase, the last committed step of the pathway, is posttranscriptionally regulated. Also, in nodules under stress conditions, gamma-glutamylcysteine synthetase was translationally regulated, but homoglutathione synthetase (mRNA and activity) and homoglutathione (content and redox state) were not affected. Most interestingly, in nodules exposed to jasmonic acid, dehydroascorbate reductase activity was posttranslationally suppressed, ascorbate oxidase showed strong transcriptional up-regulation, and dehydroascorbate content increased moderately. These changes were not due to a direct effect of jasmonic acid on the enzyme activities but might be part of the signaling pathway in the response of nodules to stress. We determined ascorbate, homoglutathione, and ascorbate-glutathione pathway enzyme activities in two senescing stages of nodules undergoing oxidative stress. When all parameters were expressed on a nodule fresh weight basis, we found that in the first stage ascorbate decreased by 60% and homoglutathione and antioxidant activities remained fairly constant, whereas in the second stage ascorbate and homoglutathione, their redox states, and their associated enzyme activities significantly decreased. The coexistence in the same plants of nodules at different senescence stages, with different ascorbate concentrations and redox states, indicates that the life span of nodules is in part controlled by endogenous factors and points to ascorbate as one of the key players.

Early developmental and stress responsive ESTs from mungbean, Vigna radiata (L.) Wilczek, seedlings

Although mungbean (Vigna radiata (L.) Wilczek) is commonly used as human food; the genomic resources of this species available in databases are limited. This study aims to develop expressed sequence tag (EST) resources for mungbean genes informative to early seedling development and chilling response. Two mungbean varieties that differ in disease resistance were found to also differ in their susceptibility to chilling temperatures. A total of 1,198 ESTs were obtained from one cDNA library and four PCR-select cDNA subtraction libraries; among these 523 were clustered into 136 contigs and 675 were singletons. The 811 non-redundant uniESTs were compared to GenBank using the Basic Local Alignment Search Tool (BLAST) and WU-BLAST algorithms, of these only 489 uniESTs had significant sequence homology, which may be involved in resuming the metabolic activity of seedlings, switching on photomorphogenesis, fuelling photosynthesis and/or initiating the unique developmental programs. Their encoded proteins may associate with regulatory proteins to trigger a direct stress response or participate in acclimation to environmental stressors. The uniEST platform reported will enrich the genomic resources of mungbean for functional genomic research on seedling development and chilling response of tropical crops and provide targets for improving the chilling tolerance of the tropical crops.

Wounding and pathogen infection induce a chloroplast-targeted lipoxygenase in the common bean (Phaseolus vulgaris L.)

Chloroplastic LOXs are implicated in the biosynthesis of oxylipins like jasmonic acid and C6 volatiles among others. In this study, we isolated the cDNA of a novel chloroplast-targeted Phaseolus vulgaris LOX, (PvLOX6). This gene is highly induced after wounding, non-host pathogen infection, and by signaling molecules as H2O2, SA, ethylene and MeJA. The phylogenetic analysis of PvLOX6 showed that it is closely related to chloroplast-targeted LOX from potato (H1) and tomato (TomLOXC); both of them are implicated in the biosynthesis of C6 volatiles. Induction of PvLOX6 mRNA by wounding ethylene and jasmonic acid on the one side, and non-host pathogen, salicylic acid on the other indicates that common bean uses the same LOX to synthesize oxylipins in response to different stresses.

Molecular cloning and characterization of phosphorus starvation responsive genes in common bean (Phaseolus vulgaris L.)

Common bean (Phaseolus vulgaris L.) is one of the most important food legumes in the world and its production is limited by low phosphate (Pi) availability in many arable soils. To gain better insight into the molecular mechanisms by which common bean adapts to low Pi availability, we generated a suppression subtractive cDNA library to identify genes involved in P starvation responses. Over 240 putative Pi starvation-responsive genes were identified. The identified clones were sequenced and BLASTx/BLASTn analysis revealed an array of 82 genes showing a high degree of sequence homology to known and unknown proteins in the database. Transcript abundance of seven genes representing different functional categories was examined by Northern blot analysis. Six genes were strongly induced/enhanced under Pi deficiency confirming the results of SSH. Full length cDNAs for three genes, representing PvIDS4-like, PvPS2, and PvPT1 were cloned and characterized. The open reading frame (ORF) of PvIDS4-like encodes a 281-amino acid protein, containing a SPX domain. The ORF of PvPS2 gene encodes a 271-amino acid protein coding for a putative phosphatase. The PvPT1 encodes a 531-amino acid protein exhibiting high homology with high affinity Pi transporters. Expression patterns of these three genes in relation to Pi availability were evaluated with two contrasting genotypes (P-inefficient Dor364 and P-efficient G19833). Both Northern and RT-PCR results showed enhanced accumulation of phosphate transporters and phosphatases in P-efficient genotype, implying that in addition to modified root morphology and architecture, increased P transport and phosphatases activity might contribute to efficient Pi acquisition and translocation in G19833 common bean genotype under limited Pi conditions.

Fast, efficient and reproducible genetic transformation of Phaseolus spp. by Agrobacterium rhizogenes

This transformation procedure generates, with high efficiency (70-90%), hairy roots in cultivars, landraces and accessions of Phaseolus vulgaris (common bean) and other Phaseolus spp. Hairy roots rapidly develop after wounding young plantlets with Agrobacterium rhizogenes, at the cotyledon node, and keeping the plants in high-humidity conditions. Callogenesis always precedes hairy-root formation, and after 15 days, when roots develop at wounded sites, the stem with the normal root is cleaved below the hairy root zone. Transgenic roots and nodules co-transformed with a binary vector can be easily identified using a reporter gene. This procedure, in addition to inducing robust transgenic hairy roots that are susceptible to being nodulated by rhizobia and to fixing nitrogen efficiently, sets the foundation for a high-throughput functional genomics approach on the study of root biology and root-microbe interactions. This protocol can be completed within 30 days.

Characterization of an 18,166 EST dataset for cassava (Manihot esculenta Crantz) enriched for drought-responsive genes

Cassava (Manihot esculenta Crantz) is a staple food for over 600 million people in the tropics and subtropics and is increasingly used as an industrial crop for starch production. Cassava has a high growth rate under optimal conditions but also performs well in drought-prone areas and on marginal soils. To increase the tools for understanding and manipulating drought tolerance in cassava, we generated expressed sequence tags (ESTs) from normalized cDNA libraries prepared from dehydration-stressed and control well-watered tissues. Analysis of a total of 18,166 ESTs resulted in the identification of 8,577 unique gene clusters (5,383 singletons and 3,194 clusters). Functional categories could be assigned to 63% of the unigenes, while another approximately 11% were homologous to hypothetical genes with unclear functions. The remaining approximately 26% were not significantly homologous to sequences in public databases suggesting that some may be novel and putatively specific to cassava. The dehydration-stressed library uncovered numerous ESTs with recognized roles in drought-responses, including those that encode late-embryogenesis-abundant proteins thought to confer osmoprotective functions during water stress, transcription factors, heat-shock proteins as well as proteins involved in signal transduction and oxidative stress. The unigene clusters were screened for short tandem repeats for further development as microsatellite markers. A total of 592 clusters contained 646 repeats, representing 3.3% of the ESTs queried. The ESTs presented here are the first dehydration stress transcriptome of cassava and can be utilized for the development of microarrays and gene-derived molecular markers to further dissect the molecular basis of drought tolerance in cassava.

Rapid evolutionary change of common bean (Phaseolus vulgaris L) plastome, and the genomic diversification of legume chloroplasts

Background

Fabaceae (legumes) is one of the largest families of flowering plants, and some members are important crops. In contrast to what we know about their great diversity or economic importance, our knowledge at the genomic level of chloroplast genomes (cpDNAs or plastomes) for these crops is limited.

Results

We sequenced the complete genome of the common bean (Phaseolus vulgaris cv. Negro Jamapa) chloroplast. The plastome of P. vulgaris is a 150,285 bp circular molecule. It has gene content similar to that of other legume plastomes, but contains two pseudogenes, rpl33 and rps16. A distinct inversion occurred at the junction points of trnH-GUG/rpl14 and rps19/rps8, as in adzuki bean [1]. These two pseudogenes and the inversion were confirmed in 10 varieties representing the two domestication centers of the bean. Genomic comparative analysis indicated that inversions generally occur in legume plastomes and the magnitude and localization of insertions/deletions (indels) also vary. The analysis of repeat sequences demonstrated that patterns and sequences of tandem repeats had an important impact on sequence diversification between legume plastomes and tandem repeats did not belong to dispersed repeats. Interestingly, P. vulgaris plastome had higher evolutionary rates of change on both genomic and gene levels than G. max, which could be the consequence of pressure from both mutation and natural selection.

Conclusion

Legume chloroplast genomes are widely diversified in gene content, gene order, indel structure, abundance and localization of repetitive sequences, intracellular sequence exchange and evolutionary rates. The P. vulgaris plastome is a rapidly evolving genome.

Identification of regulatory pathways involved in the reacquisition of root growth after salt stress in Medicago truncatula

Root growth and function are determined by the action of environmental stresses through specific genes that adapt root development to these restrictive conditions. We have defined in vitro conditions affecting the growth and recovery of Medicago truncatula roots after a salt stress. A dedicated macroarray containing 384 genes, based on a large-scale subtractive hybridization approach, was constructed and used to analyze gene expression during salt stress and recovery of root growth from this stress. Several potential regulatory genes were identified as being linked to this recovery process: a novel RNA-binding protein, a small G-protein homologous to ROP9, a receptor-like kinase, two TF IIIA-like and an AP2-like transcription factors (TF), MtZpt2-1, MtZpt2-2 and MtAp2, and a histidine kinase associated with cytokinin transduction pathways. The two ZPT2-type TFs were also rapidly induced by cold stress in roots. By analyzing transgenic M. truncatula plants showing reduced expression levels of both TFs and affected in their capacity to recover root growth after a salt stress, we identified potential target genes that were either activated or repressed in these plants. Overexpression of MtZpt2-1 in roots conferred salt tolerance and affected the expression of three putative targets in the predicted manner: a cold-regulated A (CORA) homolog, a flower-promoting factor (FPF1) homolog and an auxin-induced proline-rich protein (PRP) gene. Hence, regulatory networks depending on TFIIIA-like transcription factors are involved in the control of root adaptation to salt stress.

Phosphorus stress in common bean: root transcript and metabolic responses

Phosphorus (P) is an essential element for plant growth. Crop production of common bean (Phaseolus vulgaris), the most important legume for human consumption, is often limited by low P in the soil. Functional genomics were used to investigate global gene expression and metabolic responses of bean plants grown under P-deficient and P-sufficient conditions. P-deficient plants showed enhanced root to shoot ratio accompanied by reduced leaf area and net photosynthesis rates. Transcript profiling was performed through hybridization of nylon filter arrays spotted with cDNAs of 2,212 unigenes from a P deficiency root cDNA library. A total of 126 genes, representing different functional categories, showed significant differential expression in response to P: 62% of these were induced in P-deficient roots. A set of 372 bean transcription factor (TF) genes, coding for proteins with Inter-Pro domains characteristic or diagnostic for TF, were identified from The Institute of Genomic Research/Dana Farber Cancer Institute Common Bean Gene Index. Using real-time reverse transcription-polymerase chain reaction analysis, 17 TF genes were differentially expressed in P-deficient roots; four TF genes, including MYB TFs, were induced. Nonbiased metabolite profiling was used to assess the degree to which changes in gene expression in P-deficient roots affect overall metabolism. Stress-related metabolites such as polyols accumulated in P-deficient roots as well as sugars, which are known to be essential for P stress gene induction. Candidate genes have been identified that may contribute to root adaptation to P deficiency and be useful for improvement of common bean.

UniRef: comprehensive and non-redundant UniProt reference clusters

MOTIVATION

Redundant protein sequences in biological databases hinder sequence similarity searches and make interpretation of search results difficult. Clustering of protein sequence space based on sequence similarity helps organize all sequences into manageable datasets and reduces sampling bias and overrepresentation of sequences.

RESULTS

The UniRef (UniProt Reference Clusters) provide clustered sets of sequences from the UniProt Knowledgebase (UniProtKB) and selected UniProt Archive records to obtain complete coverage of sequence space at several resolutions while hiding redundant sequences. Currently covering >4 million source sequences, the UniRef100 database combines identical sequences and subfragments from any source organism into a single UniRef entry. UniRef90 and UniRef50 are built by clustering UniRef100 sequences at the 90 or 50% sequence identity levels. UniRef100, UniRef90 and UniRef50 yield a database size reduction of approximately 10, 40 and 70%, respectively, from the source sequence set. The reduced redundancy increases the speed of similarity searches and improves detection of distant relationships. UniRef entries contain summary cluster and membership information, including the sequence of a representative protein, member count and common taxonomy of the cluster, the accession numbers of all the merged entries and links to rich functional annotation in UniProtKB to facilitate biological discovery. UniRef has already been applied to broad research areas ranging from genome annotation to proteomics data analysis.

AVAILABILITY

UniRef is updated biweekly and is available for online search and retrieval at http://www.uniprot.org, as well as for download at ftp://ftp.uniprot.org/pub/databases/uniprot/uniref.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Agrobacterium rhizogenes transformation of the Phaseolus spp.: a tool for functional genomics

A fast, reproducible, and efficient transformation procedure employing Agrobacterium rhizogenes was developed for Phaseolus vulgaris L. wild accessions, landraces, and cultivars and for three other species belonging to the genus Phaseolus: P. coccineus, P. lunatus, and P. acutifolius. Induced hairy roots are robust and grow quickly. The transformation frequency is between 75 and 90% based on the 35-S promoter-driven green fluorescent protein and beta-glucuronidase expression reporter constructs. When inoculated with Rhizobium tropici, transgenic roots induce normal determinate nodules that fix nitrogen as efficiently as inoculated standard roots. The A. rhizogenes-induced hairy root transformation in the genus Phaseolus sets the foundation for functional genomics programs focused on root physiology, root metabolism, and root-microbe interactions.

JUICE: a data management system that facilitates the analysis of large volumes of information in an EST project workflow

Background

Expressed sequence tag (EST) analyses provide a rapid and economical means to identify candidate genes that may be involved in a particular biological process. These ESTs are useful in many Functional Genomics studies. However, the large quantity and complexity of the data generated during an EST sequencing project can make the analysis of this information a daunting task.

Results

In an attempt to make this task friendlier, we have developed JUICE, an open source data management system (Apache + PHP + MySQL on Linux), which enables the user to easily upload, organize, visualize and search the different types of data generated in an EST project pipeline. In contrast to other systems, the JUICE data management system allows a branched pipeline to be established, modified and expanded, during the course of an EST project.

The web interfaces and tools in JUICE enable the users to visualize the information in a graphical, user-friendly manner. The user may browse or search for sequences and/or sequence information within all the branches of the pipeline. The user can search using terms associated with the sequence name, annotation or other characteristics stored in JUICE and associated with sequences or sequence groups. Groups of sequences can be created by the user, stored in a clipboard and/or downloaded for further analyses.

Different user profiles restrict the access of each user depending upon their role in the project. The user may have access exclusively to visualize sequence information, access to annotate sequences and sequence information, or administrative access.

Conclusion

JUICE is an open source data management system that has been developed to aid users in organizing and analyzing the large amount of data generated in an EST Project workflow. JUICE has been used in one of the first functional genomics projects in Chile, entitled "Functional Genomics in nectarines: Platform to potentiate the competitiveness of Chile in fruit exportation". However, due to its ability to organize and visualize data from external pipelines, JUICE is a flexible data management system that should be useful for other EST/Genome projects. The JUICE data management system is released under the Open Source GNU Lesser General Public License (LGPL). JUICE may be downloaded from or .

Identification of novel drought-related mRNAs in common bean roots by differential display RT-PCR

Drought is a major constraint for the production of common bean (Phaseolus vulgaris L.). To identify molecular responses to water deficit, we performed a differential display RT-PCR (DDRT) analysis using roots of bean plants grown aeroponically and submitted to dehydration. This allowed us to visualise 1200 DDRT bands, 8.7% of which showed a clear regulation by dehydration, and to clone 42 cDNAs, called PvD1 to PvD42. Among them, 20 early-dehydration-responsive cDNAs were selected by reverse northern that were induced or repressed before detectable water status changes and induction of ABA-regulated genes. Northern analysis for 16 PvD clones confirmed these early regulations and allowed us to identify four late dehydration-responsive genes. Their putative involvement in signalling, protein turn-over and translocation, chaperones as well as root growth modulations in response to water stress is discussed.

Identification of candidate phosphorus stress induced genes in Phaseolus vulgaris through clustering analysis across several plant species

Common bean (Phaseolus vulgaris L.) is the world's most important grain legume for direct human consumption. However, the soils in which common bean predominate are frequently limited by the availability of phosphorus (P). Improving bean yield and quality requires an understanding of the genes controlling P acquisition and use, ultimately utilising these genes for crop improvement. Here we report an in silico approach for the identification of genes involved in adaptation of P. vulgaris and other legumes to P-deficiency. Some 22 groups of genes from four legume species and Arabidopsis thaliana, encoding diverse functions, were identified as statistically over-represented in EST contigs from P-stressed tissues. By combining bioinformatics analysis with available micro / macroarray technologies and clustering results across five species, we identified 52 P. vulgaris candidate genes belonging to 19 categories as induced by P-stress response. Transport-related, stress (defence and regulation) signal transduction genes are abundantly represented. Manipulating these genes through traditional breeding methodologies and / or biotechnology approaches may allow us to improve crop P-nutrition.

A novel fix- symbiotic mutant of Lotus japonicus, Ljsym105, shows impaired development and premature deterioration of nodule infected cells and symbiosomes

Nitrogen-fixing symbiosis between legume plants and rhizobia is established through complex interactions between two symbiotic partners. To identify the host legume genes that play crucial roles in such interactions, we isolated a novel Fix- mutant, Ljsym105, from a model legume Lotus japonicus MG-20. The Ljsym105 plants displayed nitrogen-deficiency symptoms after inoculation with Mesorhizobium loti under nitrogen-free conditions, but their growth recovered when supplied with nitrogen-rich nutrients. Ljsym105 was recessive and monogenic and mapped on the upper portion of chromosome 4. The mutant Ljsym105 formed an increased number of small and pale-pink nodules. Nitrogenase (acetylene reduction) activity per nodule fresh weight was low but retained more than 50% of that of the wild-type nodules. Light and electron microscopic observations revealed that the Ljsym105 nodule infected cells were significantly smaller than those of wild-type plants, contained enlarged symbiosomes with multiple bacteroids, and underwent deterioration of the symbiosomes prematurely as well as disintegration of the whole infected cell cytoplasm. These results indicate that the ineffectiveness of the Ljsym105 nodules is primarily due to impaired growth of infected cells accompanied with the premature senescence induced at relatively early stages of nodule development. These symbiotic phenotypes are discussed in respect to possible functions of the LjSym105 locus in the symbiotic interactions required for establishment of the nitrogen-fixing symbiosis.

Sorghum expressed sequence tags identify signature genes for drought, pathogenesis, and skotomorphogenesis from a milestone set of 16,801 unique transcripts

Improved knowledge of the sorghum transcriptome will enhance basic understanding of how plants respond to stresses and serve as a source of genes of value to agriculture. Toward this goal, Sorghum bicolor L. Moench cDNA libraries were prepared from light- and dark-grown seedlings, drought-stressed plants, Colletotrichum-infected seedlings and plants, ovaries, embryos, and immature panicles. Other libraries were prepared with meristems from Sorghum propinquum (Kunth) Hitchc. that had been photoperiodically induced to flower, and with rhizomes from S. propinquum and johnsongrass (Sorghum halepense L. Pers.). A total of 117,682 expressed sequence tags (ESTs) were obtained representing both 3' and 5' sequences from about half that number of cDNA clones. A total of 16,801 unique transcripts, representing tentative UniScripts (TUs), were identified from 55,783 3' ESTs. Of these TUs, 9,032 are represented by two or more ESTs. Collectively, these libraries were predicted to contain a total of approximately 31,000 TUs. Individual libraries, however, were predicted to contain no more than about 6,000 to 9,000, with the exception of light-grown seedlings, which yielded an estimate of close to 13,000. In addition, each library exhibits about the same level of complexity with respect to both the number of TUs preferentially expressed in that library and the frequency with which two or more ESTs is found in only that library. These results indicate that the sorghum genome is expressed in highly selective fashion in the individual organs and in response to the environmental conditions surveyed here. Close to 2,000 differentially expressed TUs were identified among the cDNA libraries examined, of which 775 were differentially expressed at a confidence level of 98%. From these 775 TUs, signature genes were identified defining drought, Colletotrichum infection, skotomorphogenesis (etiolation), ovary, immature panicle, and embryo.

A 6374 unigene set corresponding to low abundance transcripts expressed following fertilization in Solanum chacoense Bitt, and characterization of 30 receptor-like kinases

In order to characterize regulatory genes that are expressed in ovule tissues after fertilization we have undertaken an EST sequencing project in Solanum chacoense, a self-incompatible wild potato species. Two cDNA libraries made from ovule tissues covering embryo development from zygote to late torpedo-stage were constructed and plated at high density on nylon membranes. To decrease EST redundancy and enrich for transcripts corresponding to weakly expressed genes a self-probe subtraction method was used to select the colonies harboring the genes to be sequenced. 7741 good sequences were obtained and, from these, 6374 unigenes were isolated. Thus, the self-probe subtraction resulted in a strong enrichment in singletons, a decrease in the number of clones per contigs, and concomitantly, an enrichment in the total number of unigenes obtained (82%). To gain insights into signal transduction events occurring during embryo development all the receptor-like kinases (or protein receptor kinases) were analyzed by quantitative real-time RT-PCR. Interestingly, 28 out of the 30 RLK isolated were predominantly expressed in ovary tissues or young developing fruits, and 23 were transcriptionaly induced following fertilization. Thus, the self-probe subtraction did not select for genes weakly expressed in the target tissue while being highly expressed elsewhere in the plant. Of the receptor-like kinases (RLK) genes isolated, the leucine-rich repeat (LRR) family of RLK was by far the most represented with 25 members covering 11 LRR classes.