SummaryLegumes played central roles in the development of agriculture and civilization, and today account for approximately one-third of the world's primary crop production. Unfortunately, most cultivated legumes are poor model systems for genomic research. Therefore, Medicago truncatula, which has a relatively small diploid genome, has been adopted as a model species for legume genomics. To enhance its value as a model, we have generated a gene expression atlas that provides a global view of gene expression in all major organ systems of this species, with special emphasis on nodule and seed development. The atlas reveals massive differences in gene expression between organs that are accompanied by changes in the expression of key regulatory genes, such as transcription factor genes, which presumably orchestrate genetic reprogramming during development and differentiation. Interestingly, many legume-specific genes are preferentially expressed in nitrogen-fixing nodules, indicating that evolution endowed them with special roles in this unique and important organ. Comparative transcriptome analysis of Medicago versus Arabidopsis revealed significant divergence in developmental expression profiles of orthologous genes, which indicates that phylogenetic analysis alone is insufficient to predict the function of orthologs in different species. The data presented here represent an unparalleled resource for legume functional genomics, which will accelerate discoveries in legume biology.

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.