Evaluation of diversity among common beans (Phaseolus vulgaris L.) from two centers of domestication using 'omics' technologies
© Mensack et al; licensee BioMed Central Ltd. 2010
Received: 10 June 2010
Accepted: 2 December 2010
Published: 2 December 2010
Genetic diversity among wild accessions and cultivars of common bean (Phaseolus vulgaris L.) has been characterized using plant morphology, seed protein allozymes, random amplified polymorphic DNA, restriction fragment length polymorphisms, DNA sequence analysis, chloroplast DNA, and microsatellite markers. Yet, little is known about whether these traits, which distinguish among genetically distinct types of common bean, can be evaluated using omics technologies.
Three 'omics' approaches: transcriptomics, proteomics, and metabolomics were used to qualitatively evaluate the diversity of common bean from two Centers of Domestication (COD). All three approaches were able to classify common bean according to their COD using unsupervised analyses; these findings are consistent with the hypothesis that differences exist in gene transcription, protein expression, and synthesis and metabolism of small molecules among common bean cultivars representative of different COD. Metabolomic analyses of multiple cultivars within two common bean gene pools revealed cultivar differences in small molecules that were of sufficient magnitude to allow identification of unique cultivar fingerprints.
Given the high-throughput and low cost of each of these 'omics' platforms, significant opportunities exist for their use in the rapid identification of traits of agronomic and nutritional importance as well as to characterize genetic diversity.
Common bean (Phaseolus vulgaris L.) is one of the oldest cultivated crops in the Americas and is the most important grain legume for human consumption with production more than double that of the second most important grain legume, chickpea (Cicer arietinum L.) . Common bean was domesticated in the Americas by indigenous people during pre-Colombian times. Archeological data suggest that bean was independently domesticated in different regions of the Americas including, the Andean region of South America , Argentina , and Mexico [4, 5]. The oldest domesticated beans were found at archeological sites in each of these regions between 4300 and 8000 B.P. [6, 7] Changes in bean plant phenotype as a result of domestication include but are not limited to growth habit, seed size, seed retention, and maturity[3, 8]. However, the molecular events that underlie these differences in agronomic traits have not been elucidated.
Current tools used to characterize genetic diversity within common bean include plant morphology , seed protein allozymes [11–14], random amplified polymorphic DNA, restriction fragment length polymorphisms , DNA sequence analysis , chloroplast DNA , and microsatellite markers [19, 20]. While effective for evaluating genetic diversity of parent plants, these techniques are of limited use in breeding programs focused on crop improvement. Methods that can assess potential differences in biological function among cultivars on a genome wide basis have become available in high-throughput and low-cost formats. These 'omics' platforms provide an unprecedented opportunity to ultimately identify the underlying molecular mechanisms that account for traits of agronomic and nutritional importance.
The research presented herein utilized three 'omics' platforms (transcriptomics, proteomics, and metabolomics) to determine whether beans from Middle American and Andean COD differ in their patterns of gene transcription, protein expression, and/or small molecule synthesis and metabolism. To our knowledge, this is the first report that shows differences in gene transcription, protein expression, and metabolite profiles between beans from two COD using 'omics' techniques. Moreover the metabolomic approach was used to evaluate whether common bean cultivars within the same gene pool, referred to as a market class, differ in their profile of metabolites. These proof-in-principle experiments herald significant opportunities for utilizing 'omic' technologies in the rapid identification of traits with agronomic and nutritional importance that could serve to guide common bean breeding programs.
Three 'omics' platforms to distinguish common bean COD
Experiments were conducted to determine whether or not the three 'omics' technology platforms (transcriptomics, proteomics, and metabolomics) would be capable of distinguishing dry bean cultivars from each of the two COD. The approach used was sequential; transcriptomics experiments were performed first, followed by proteomics, and then metabolomics, with increasing complexity of the experimental design as work progressed to each new platform. The white kidney bean and navy bean market classes representing the Andean and Middle American gene pools were selected as representative of the respective COD. Furthermore, both market classes are commercially important, and had a white pigmented seed coat thus limiting the likelihood of identifying qualitative differences between COD that were due solely to synthesis of different pigments.
Phaseolus vulgaris genes differentially expressed between navy and white kidney bean seeds two weeks after flowering.
Annotation (BLAST hit and organism)
calmodulin 2 [Medicago truncatula]
annexin [Medicago truncatula]
peroxidase [Glycine max]
RUBISCO small chain 1 precursor [Glycine max]
nuM1 [Medicago sativa]
ubiquitin [Glycine max]
ubiquitin [Lycopersicon esculentum]
pectin methylesterase-like protein [Arabidopsis thaliana]
ubiquitin [Glycine max]
polyubiquitin containing 7 ubiquitin monomers
putative transcriptional coactivator [Brassica rapa]
heat shock factor protein 3 (HSF3)/heat shock transcription factor 3 (HSTF3) [Arabidopsis thaliana]
glutamyl-tRNA synthetase [Arabidopsis thaliana]
translation initiation factor IF-1 [Glycine max]
putative small nuclear ribonucleoprotein polypeptide E [Oryza sativa (japonica cultivar-group)]
elongation factor 1-alpha, putative [Arabidopsis thaliana]
ribosomal protein [Petunia × hybrida]
40 S ribosomal protein S21 (RPS21C) [Arabidopsis thaliana]
putative ribosomal protein [Capsicum annuum]
putative protein [Arabidopsis thaliana]
60 S ribosomal protein L24 [Prunus avium]
40 s ribosomal protein S23 [Euphorbia esula]
putative ribosomal protein S29 [Oryza sativa (japonica cultivar-group)]
ubiquitin fusion protein/40 S ribosomal protein S27a [Zea mays]
zinc finger (C3HC4-type RING finger) protein-related [Arabidopsis thaliana]
unknown protein [Arabidopsis thaliana]
expressed protein [Arabidopsis thaliana]
farnesylated protein [Cicer arietinum]
putative leucine-rich repeat protein [Arabidopsis thaliana]
LIM domain-containing protein [Arabidopsis thaliana]
ripening-induced protein [Fragaria vesca]
Phaseolus vulgaris genes differentially expressed between navy and white kidney bean seeds three weeks after flowering.
Ratio WK:NBa, b
Annotation (BLAST hit and organism)
nucleoside diphosphate kinase [Glycine max]
nitrite transporter [Cucumis sativus]
oxygen evolving complex 33 kDa photosystem II protein [Nicotiana tabacum]
phototropic-responsive NPH3 family protein [Arabidopsis thaliana]
signal transducer activity
unknown [Arabidopsis thaliana]
unnamed protein product [Arabidopsis thaliana]
Based on the transcript profiling data reported in the previous section, it was predicted that the profile of expressed protein would differ between cultivars representative of COD. Protein expression was assessed using a two-dimensional gel system instead of an LC-MS platform, so that a visual picture of similarities and differences in types and amounts of proteins expressed between COD could be obtained. Moreover, rather than limiting the analysis to one cultivar from each market class-COD category, three cultivars from each genepool were evaluated. The list of cultivars studied can be found in the Materials and Methods section. Master gels of the two dry bean germplasm (white kidney and navy) were created for purposes of exploratory analysis. Close visual inspection showed distinct differences in both the absolute number of proteins as well as differences in concentration (spot density) between COD. With the false discovery rate set at zero, 2186 spots were selected from the master gel using the PDQuest software.
Cultivated varieties distinguished by metabolomic fingerprinting
The elucidation of molecular differences underlying traits that permit the division of common bean germplasm into distinct centers of domestication has been largely unexplored by high-throughput 'omic' technologies. A triad of 'omic' technologies: transcriptomics, proteomics, and metabolomics, were used herein to determine whether differences in gene transcription, protein expression, or synthesis/metabolism of small molecules were of sufficient magnitude to distinguish among bean cultivars within and between two COD (Figure 1). As shown in Figure 2b, 3c, and 5a, distinct differences were observed among cultivars and between the two COD in all three 'omics' signatures.
The transcript profiling experiments revealed that chronological age of the seed from anthesis has a significant impact on the number of differentially expressed genes and the functions (gene ontology) associated with those genes. However, examination of the gene lists in Table 1 failed to provide a clear picture of the physiological distinctions between COD. One limitation to consider is that distinct P. vulgaris genes may exist without a soybean homologue to utilize during hybridization and thus a cross-species microarray may exclude unique dry bean transcripts. Reciprocally, the chance of false positives was increased in the cross-species use of the cDNA spotted arrays because of reduced hybridization stringency. The drawbacks associated with using a soybean array for common bean transcriptomics stem from reduced synteny between P. vulgaris and G. max. The most recent common ancestor for P. vulgaris and G. max was 19 Ma. At some time point following divergence, diploidization of the G. max genome occurred making its chromosome number 2n = 80 compared to 2n = 40 of P. vulgaris. Genomic rearrangement including translocation or gene loss during diploidization of the soybean leads to reduced synteny between the two species and must be considered when attempting to determine chromosome overlap between the soybean and the common bean for genetic mapping. The use of the Affymetrix array platform as reported in Yang et al.  may serve to overcome some of the noted limitations, and it is worthwhile to note that the Whole-Genome Sequence of Common Bean project and the Ibero-American whole genome sequencing project include both COD . With the entire common bean genome sequenced, a complete dry bean microarray can be synthesized and annotated to advance the transcriptomics platform Challenges in gene list interpretation also have been encountered using other organisms and have led to the development of alternative approaches to data analysis methods such as Gene Set Enrichment Analysis (GSEA) . GSEA permits the investigator to use existing knowledge to test specific hypotheses about differences in pathways and networks of gene expression between two or more treatment conditions. Efforts to incorporate GSEA for market class stratification were limited by cross-species hybridization and the lack of dry bean and soybean gene lists. Consequently, no useful data was gained. This situation underscores the importance of current efforts to sequence the Phaseolus genome. When that work is completed, it will be possible to generate hypotheses for additional testing based on observations about plant morphology  or seed protein allozymes [11–13] that have been used to characterize genetic diversity within common bean. In addition, as work on metabolomic analyses continues and new compound libraries are constructed, metabolomic data such as that shown in Figure 4 and 5 will identify specific small molecules that distinguish cultivars within each COD and within a market class-COD grouping. That information will permit future hypothesis generation about the biosynthetic pathways that are induced or repressed to account for differences in metabolite profiles. Since differences in small molecule biosynthesis are likely to result not only from differences in gene transcription but also arise from differences in the translation of transcripts and post translational modifications of proteins, the potential contributions of proteomic analyses such as those shown in Figure 3 are clear. For future analyses, we recommend that an HPLC or UPLC separation approach be employed rather than 2DGE to interrogate the Phaseolus proteome because of the speed, throughput, and relative cost of the LC-based platform.
As one considers quantitative differences among the three 'omics' signatures that distinguished between the Andean and Middle American COD, it is readily apparent that marked differences exist in the number of features that distinguish between COD when tallying results across platforms. However we caution the reader that the purpose of this experimental approach was primarily the qualitative evaluation of distinguishing differences in transcript, protein, and/or metabolite expression between COD. While the next generation of experiments would be well served to use seeds matched for physiological age across all omics platforms, these initial experiments were performed using available mature seed with the exception of the initial work on transcript expression which required embryos to be harvested in the field and immediately frozen in liquid nitrogen.
Many factors other than the biology of the system could account, at least in part, for the numeric differences observed. For example, the use of the soy-array to identify differentially expressed transcripts or the use of only one of an almost endless number of assay conditions that have been developed for proteomic or metabolomic analyses are likely to make the numeric differences noted unreliable. Nonetheless, the magnitude of the observed differences indicates the potential merit of these approaches as a fertile source of information for plant breeders.
Overall, the 'omics' chosen for a particular set of experiments should depend on the question being asked. While considerable information can be gleaned from each 'platform, the biological activity assessed is interactive, and accordingly, the greatest insights are likely to emerge through the integration of complimentary data sets. Transcriptomics, proteomics, and metabolomics are just three of these .
We have shown that two genetically diverse dry bean germplasm can be easily differentiated using a suite of three 'tools (transcriptomics, proteomics, and metabolomics). Using this panel of techniques we have provided a glimpse into a signature which can be used to determine COD. Furthermore, known genetic distances among cultivars and between COD were validated using proteomics and metabolomics. 'Omic' signatures unique to P. vulgaris germplasm may be useful to assess complex traits or reactions to biotic and abiotic stress and to incorporate genetic diversity in breeding efforts. This work compliments other techniques currently employed by breeders to assess dry bean genetic diversity. As signatures are further developed in future work, they have the potential to guide those engaged in crop improvement for selecting traits of agronomic and nutritional importance.
Navy bean (cv. Seahawk, cv. Norstar, and cv. Vista) and white kidney bean (cv. Beluga, cv. Silver Cloud, and cv. Lassen) cultivars representative of the two market classes were selected for this study to represent the Middle American and Andean COD, respectively. Seeds from each of the market classes used in this study were obtained from field grown plants at the Colorado State University Agricultural Research, Development and Education Center (ARDEC), Fort Collins, CO. The bean seed used for analysis within each omics approach was of the same physiological age but differed across omics platforms. Seeds for transcriptomics analyses were collected from plants at two and three weeks after anthesis (plant flowering), immediately frozen in liquid nitrogen and stored at -80°C. Seeds for proteomics and metabolomics were harvested at full maturity and air dried and stored at room temperature (RT, 22 ± 2°C) until use.
RNA was extracted from navy and white kidney bean seeds from field grown plants in 2007. Seeds were collected in the field and immediately frozen as described above. Two cultivars, Silver Cloud (white kidney) and Norstar (navy) were selected for spotted array gene expression analysis. After storage, frozen seeds were separated from pods and RNA was isolated according to the TRIzol reagent protocol (Life Technologies, Gaithersburg, MD). The RNA sample was purified with the RNeasy Mini Kit (Qiagen, Valencia, CA) according to Affymetrix (Santa Clara, CA) instructions. RNA concentrations and purity were determined using an ND-1000 NanoDrop spectrophotometer (Thermo Scientific NanoDrop, Wilmington, DE). RNA integrity was evaluated by the Experion Bioanalyzer Automated Electrophoresis System (Bio-Rad Laboratories, Hercules, CA).
Since a commercial Phaseolus vulgaris microarray was not available, a soybean (Glycine max L. Merr.) array was selected because it was the closest phylogenetic relative among available arrays . Although Glycine max and Phaseolus vulgaris differ in chromosome number and genome size (the soybean genome is twice as large as common dry bean), linkage mapping of DNA markers found an average conserved block length of 13.9 cM between the two genomes indicating high conversation and preservation [29, 30]. Twenty-five μg of RNA was used to prepare complimentary DNA for spotted 2-color microarray analysis. Complimentary DNA (cDNA) was synthesized using the Superscript III kit (Invitrogen, Carlsbad, CA) and hybridized to a microarray developed for soybean  using the Genisphere Array 50 kit (Genisphere, Hatfield, PA). Briefly, RNA (29.5 μL) was mixed with dye-appropriate RT primer (1.5 μL) and heated to 80°C for 10 min, and transferred to ice for 2 min. Superase-in RNase inhibitor (1 μL) was added before 18 μL of reaction mix composed of 5× Superscript III first strand buffer (10 μL), 10 mM dNTP mix (2 μL), 0.1 M dithiotreitol (4 μL) and Superscript III enzyme (4 μL). After incubation at RT for 5 min, tubes were placed in an iQ iCycler (Bio-Rad) for 2 h at 50°C. The reaction was stopped with 0.5 M NaOH/50 mM EDTA (7 μL) followed by incubation at 65°C for 10 min and neutralized with 1 M Tris-HCl (10 μL). cDNA from the navy and white kidney bean samples were combined and purified using Geneclean Turbo kit (Qbiogene, Carlsbad, CA). The purified cDNA was stored at -20°C until use. Arrays were cross-linked by exposure to 65 mJ of UV irradiation and blocked by incubation in 250 mL prehybridization solution [5× SSC (18% Sodium Chloride, 9% Sodium Citrate) buffer, 0.1% SDS and 0.01% BSA solution] at 42°C for 1-2 h. The arrays were dried by centrifugation at RT for 5 min and placed into hybridization chambers and incubated at 42°C for 15 min. Formamide-based buffer (2×, 50 μL) and LNA dT blocker (2 μL) were added to the cDNA sample and heated for 10 min at 80°C. A cover slip was placed on the array and the cDNA mix was introduced between the cover slip and the array surface. The solution was distributed evenly on the array by capillary action. Water (15 μL) was added to the wells of the array chambers followed by incubation for 16 h at 42°C. The array was vigorously washed 3 times for 2 min, first with 2× SSC and 0.2% SDS (250 mL) followed by 2× SSC (250 mL) and then with 0.2× SSC (250 mL). The arrays were centrifuged for 3 min (1000 × g) to dry and warmed at 42°C for 15 min. A mix of 2× formamide-based buffer (42 μL), water (35 μL) and 3.5 μL each of Cy3 and Cy5 was prepared in the dark and incubated at 80°C for 10 min. The hybridization mix was then injected onto the soybean cDNA microarray. Arrays were incubated at 42°C for 3 h. After the second hybridization step, the arrays were washed and dried as described above and immediately scanned with a Genepix 4000B scanner (Molecular Devices, Sunnyvale, CA).
Scanned images were analyzed using GenePix Pro 6.0 software (Molecular Devices) where spot features were aligned for annotation. Poorly hybridized spots were discarded from analysis using the software parameters and the intensities of each dye were quantified. Fluorescence intensities were normalized and exported for statistical analyses.
Protein isolation and two-dimensional gel electrophoresis (2DGE)
2DGE was performed with bean extracts according to previously published methods . Approximately 15 mg of bean seed samples from each of the six cultivars listed above were suspended in sample lysis buffer, containing 7 M urea, 2 M thiourea, 40 mM DTT, 2% CHAPS, 1% Pharmalyte pH 3-10 (GE Healthcare, Piscataway, NJ), and trace amount of bromophenol blue (BPB). Insoluble debris was removed after centrifugation for 30 min (1000 × g, 15°C). The supernatant was recovered and used for analysis. Protein concentration was quantified using the Bradford assay.
Bean seed protein (200 μg) was applied to reswelled, immobilized pH gradient (IPG) strips (pH 4-7L, 24 cm; GE Healthcare). The isoelectrofocusing conditions were as follows: 10V to 300 V for 3 h and at 5000 V for a total of 95 kVh at 20°C. After isoelectrofocusing, the strips were equilibrated in a buffer containing 30% glycerol, 6 M urea, 2% SDS, 10 mg/mL dithiothreitol (DTT) and 0.05 M Tris-HCl (pH 6.8) for 15 min and then for an additional 15 min in equilibration buffer in which 42.5 mg/mL iodoacetamide replaced the DTT. The strips were positioned at the top of 13-16% gradient polyacrylamide gels with 0.5% agarose containing Laemmli sample buffer . SDS-PAGE was performed in Laemmli electrophoresis buffer at 150 V at RT. Proteins were stained with alkaline ammoniacal silver staining  scanned with a GS-800 Calibrated Densitometer (Bio-Rad) and analyzed using PDQuest v7.1.1 software (Bio-Rad).
Analysis of gels
Samples were analyzed using 2DGE to obtain quantitative protein profiles within the molecular weight range of 25 to 150 kDa. Each bean cultivar was run in duplicate. Protein spots were automatically detected using PDQuest v7.1.1 software. All spots were also manually confirmed. Images of duplicate gels were superimposed and a master gel generated for each of the six cultivars. Proteins were quantified using spot densitometry.
Comparisons of the 2DGE protein patterns generated an inclusion list for only those proteins that differed significantly between the two COD (greater than 2-fold change). Protein patterns were also compared between cultivars within each COD. Although out of the scope of this work, these spots can be used in the future for protein identification using MALDI TOF MS.
Metabolite extraction and analysis
It is imperative that the majority of the proteins within the dry bean samples are removed prior to metabolomic fingerprinting in order to avoid confounded results from possible protein fragments. The dried bean seeds were boiled to denature proteins and freeze dried for storage. The proteins were then precipitated during the metabolite extraction process using cold ethanol. The boiling procedure is as follows: 0.5 kg of bean was soaked in distilled water for 3 h at RT. After 3 h the beans were drained, rinsed thoroughly with deionized water, and blanched for 5 min at 93°C. Beans were boiled for 60 min in a 1.5% KCl solution using a pressure cooker. Finally, beans were freeze dried (Genesis SQ25LL, Virtis Company, Gardiner, NY) to powder form and stored at -80°C until use.
Metabolites were extracted using cold ethanol (65%, -20°C). Approximately 2.5 g of freeze dried bean powder was added to ethanol (65%, 50 mL) and vortexed to ensure complete mixing. The mixture was sonicated at RT for 2 h and centrifuged (1000 × g, 10 min) to separate the insoluble material from the ethanol extract. The extract was decanted into a clean conical tube and stored at -20°C until analysis up to 1 mo.
Ultra Performance Liquid chromatography-Mass Spectrometry (UPLC-MS)
Sample separation was performed using an Acquity UPLC under the control of MassLynx software (Waters, Millford, MA, USA). The sample set was randomized and held in an 8°C sample manager during the analysis. For each chromatographic run, a 1 μL sample injection was loaded to a 1.0 × 100 mm Waters Acquity UPLC BEH C18 column with 1.7 μm particle size held at 40°C. Separation was performed by reverse phase chromatography at a flow rate of 0.15 ml/min. The eluent consisted of water and methanol (Fisher, Optima® LC/MS grade) supplemented with formic acid (Fluka, LC/MS grade) in the following proportions: Solvent A = 95:5 water:methanol + 0.1% formic acid; Solvent B = 5:95 water:methanol + 0.1% formic acid. The separation method is described as follows: 0.1 min hold at 100% A, 14.9 min linear gradient to 100% B, 5 min hold at 100% B, 1 min linear gradient to 100% A, and 1 min hold at 100% A. A blank injection of water and 15 min chromatographic run was preformed between samples to eliminate possible carryover of analytes and to re-equilibrate the column. This cleaning method is described as follows: 0.1 min hold at 100% A, 2.9 min linear gradient to 100% B, 1 min hold at 100% B, 3 min linear gradient to 100% A, and 8 min hold at 100% A for equilibration. The flow rate for all steps was held at 0.15 ml/min.
Eluate was directed to a Q-TOF Micro quadrupole orthogonal acceleration time-of-flight mass spectrometer controlled with MassLynx software (Waters/MicroMass, Millford, MA, USA) using electrospray ionization in the positive mode (ESI+). Mass data were collected between 50 and 1000 m/z at a rate of two scans per second with a 0.1 second interscan delay. The voltage and temperature parameters were tuned for general profiling as follows: capillary = 3000 V; sample cone = 30 V; extraction cone = 2.0 V; desolvation temperature = 300°C; and source temperature = 130°C. Mass spectral peaks were centered during acquisition producing centroid data. Leucine Enkephalin was infused via an orthogonal ESI probe and baffle system (LockMass) which allowed reference ions to be detected for a single half-second scan every 10 s in an independent data collection channel. The standard mass was averaged across 10 scans providing a continuous reference for mass correction of analyte data.
Chromatographic and spectral LC-MS peaks were detected, extracted, and aligned using MarkerLynx software (Waters, Millford, MA, USA). Chromatographic peaks were detected between 0 and 18 min with a retention time error window of 0.1 min. Apex track peak detection parameters were used, automatically detecting peak width and baseline noise. No smoothing was applied. To reduce the detection and inclusion of noise as data, an intensity threshold value of 40 and a noise elimination value of 6 were used. Mass spectral peaks were detected between 50 and 1000 m/z with a mass error window of 0.07 m/z. The de-isotoping function was enabled to eliminate the inclusion of isotopic peaks. A matrix of features as defined by retention time and mass was generated, and the relative intensity of all features, as determined by area, was calculated for all individual samples. Potential effects of technical variability were minimized by normalizing the intensity values to the total ion current (TIC) such that the summation of all feature intensities in each individual sample were equal.
For microarray analysis, data were imported into Partek Discovery Suite software (Partek, St. Louis, MO), PC and 1-way ANOVA (random effects) (p < 0.01) analyses were performed. For 2DGE, PC, hierarchal clustering based on Euclidian distances, and 1-way ANOVA (random effects) (p < 0.05) analyses were carried out using Partek Discovery Suite software. Finally, the LC-MS feature matrix was mean centered and imported into SIMCA-P+ software (Umetrics, Inc., Umeå, Sweeden). PC analysis was performed using Pareto scaling. Hierarchal clustering based on Euclidian distance and 1-way ANOVA (random effects) (p < 0.05) analyses were performed using Partek Discovery Suite software. Fold change for both proteomics and metabolomics was calculated using Partek Discovery Suite.
We thank Anne Hess for guidance regarding transcriptomics data transformation and analysis, Sophie Hermann for development of the microarray experiment protocols, SungGu Lee for his technical assistance in proteomic analyses, and Erica Daniell for harvesting the beans from ARDEC fields and RNA extraction for transcriptomics analysis. The work presented here is funded in part by the Colorado Agricultural Experiment Station, Fort Collins, Colorado and by AICR (grant #08A032).
- Gepts P, Aragao FL, deBarros E, Blair MW, Brondani R, Broughton W, Galasso I, Hernandez G, Kami J, Lariguet P, et al: Genomics of Phaseolus Beans, a Major Source of Dietary Protein and Micronutrients in the Tropics. Genomics of Tropical Crop Plants. Edited by: Moore PH, Ming R. 2008, Philadelphia: Springer, 113-143. full_text.View Article
- Kaplan L, Kaplan LN: Phaseolus in archaeology. Genetic Resources of Phaseolus Beans. Edited by: Gepts P Dordrecht. 1988, Netherlands: Kluwer Academic Publishers, 125-142.View Article
- Gepts P, Debouck DG: Origin, domestication, and evolution of the common bean, Phaseolus vulgaris. Common Beans: Research for Crop Improvement. Edited by: Voysest O, Van Schoonhoven A Oxon. 1991, UK: CAB International, 7-53.
- Kaplan L: Archeology and domestication in American Phaseolus (beans). Economic Botany. 1965, 19: 358-368.View Article
- Kaplan L, McNeish RS: Prehistoric bean remains from caves in the Ocampo region of Tamaulipas, Mexico. 1960, 19: 33-35.
- Kaplan L, Lynch TF: Phaseolus (Fabaceae) in archaeology: AMS radiocarbon dates and their significance for pre-Colombian agriculture. Economic Botany. 1999, 53 (3): 261-272.View Article
- Piperno DR, Dillehay TD: Starch grains on human teeth reveal early broad crop diet in northern Peru. Proceedings of the National Academy of Sciences of the United States of America. 2008, 105 (50): 19622-19627. 10.1073/pnas.0808752105.PubMed CentralPubMedView Article
- Koinange EMK, Singh SP, Gepts P: Genetic control of the domestication syndrome in common bean. Crop Science. 1996, 36 (4): 1037-1045. 10.2135/cropsci1996.0011183X003600040037x.View Article
- Singh SP, Gepts P, Debouck DG: Races of Common Bean (Phaseolus-Vulgaris, Fabaceae). Economic Botany. 1991, 45 (3): 379-396.View Article
- Kwak M, Kami JA, Gepts P: The Putative Mesoamerican Domestication Center of Phaseolus vulgaris Is Located in the Lerma-Santiago Basin of Mexico. Crop Science. 2009, 49 (2): 554-563. 10.2135/cropsci2008.07.0421.View Article
- Singh SP, Nodari R, Gepts P: Genetic Diversity in Cultivated Common Bean .1. Allozymes. Crop Science. 1991, 31 (1): 19-23. 10.2135/cropsci1991.0011183X003100010004x.View Article
- Gepts P, Bliss FA: Phaseolin Variability Among Wild and Cultivated Common Beans (Phaseolus-Vulgaris) from Colombia. Economic Botany. 1986, 40 (4): 469-478.View Article
- Gepts P, Kmiecik K, Pereira P, Bliss FA: Dissemination Pathways of Common Bean (Phaseolus-Vulgaris, Fabaceae) Deduced from Phaseolin Electrophoretic Variability .1. the Americas. Economic Botany. 1988, 42 (1): 73-85.View Article
- Koenig R, Gepts P: Allozyme Diversity in Wild Phaseolus-Vulgaris - Further Evidence for 2 Major Centers of Genetic Diversity. Theoretical and Applied Genetics. 1989, 78 (6): 809-817. 10.1007/BF00266663.PubMedView Article
- Freyre R, Rios R, Guzman L, Debouck DG, Gepts P: Ecogeographic distribution of Phaseolus spp (Fabaceae) in Bolivia. Economic Botany. 1996, 50 (2): 195-215.View Article
- Velasquez VLB, Gepts P: Rflp Diversity of Common Bean (Phaseolus-Vulgaris) in Its Centers of Origin. Genome. 1994, 37 (2): 256-263. 10.1139/g94-036.PubMedView Article
- Kami J, Velasquez VB, Debouck DG, Gepts P: Identification of Presumed Ancestral Dna-Sequences of Phaseolin in Phaseolus-Vulgaris. Proceedings of the National Academy of Sciences of the United States of America. 1995, 92 (4): 1101-1104. 10.1073/pnas.92.4.1101.PubMed CentralPubMedView Article
- Chacon SMI, Pickersgill B, Debouck DG, Arias JS: Phylogeographic analysis of the chloroplast DNA variation in wild common bean (Phaseolus vulgaris L.) in the Americas. Plant Systematics and Evolution. 2007, 266 (3): 175-195. 10.1007/s00606-007-0536-z.View Article
- Blair MW, Giraldo MC, Buendia HF, Tovar E, Duque MC, Beebe SE: Microsatellite marker diversity in common bean (Phaseolus vulgaris L.). Theoretical and Applied Genetics. 2006, 113 (1): 100-109. 10.1007/s00122-006-0276-4.PubMedView Article
- Blair MW, Pedraza F, Buendia HF, Gaitan-Solis E, Beebe SE, Gepts P, Tohme J: Development of a genome-wide anchored microsatellite map for common bean (Phaseolus vulgaris L.). Theoretical and Applied Genetics. 2003, 107 (8): 1362-1374. 10.1007/s00122-003-1398-6.PubMedView Article
- Yang SS, Valdes-Lopez O, Xu WW, Bucciarelli B, Gronwald JW, Hernandez G, Vance CP: Transcript profiling of common bean (Phaseolus vulgaris L.) using the GeneChip Soybean Genome Array: optimizing analysis by masking biased probes. BMC Plant Biol. 2010, 10: 85-10.1186/1471-2229-10-85.PubMed CentralPubMedView Article
- Bar-Or C, Czosnek H, Koltai H: Cross-species microarray hybridizations: a developing tool for studying species diversity. Trends Genet. 2007, 23 (4): 200-207. 10.1016/j.tig.2007.02.003.PubMedView Article
- Vodkin LO, Khanna A, Shealy R, Clough SJ, Gonzalez DO, Philip R, Zabala G, Thibaud-Nissen F, Sidarous M, Stromvik MV, et al: Microarrays for global expression constructed with a low redundancy set of 27,500 sequenced cDNAs representing an array of developmental stages and physiological conditions of the soybean plant. BMCGenomics. 2004, 5 (1): 73-
- Wink M: Plant secondary metabolism: Diversity, function and its evolution. Nat Prod Commun. 2008, 3 (8): 1205-1216.
- McClean P, Lavin M, Gepts P, Jackson S: Phaseolus vulgaris: A Diploid Model for Soybean. Genetics and Genomics of Soybean. Edited by: Stacey G. 2008, New York: Springer, 55-76. full_text.View Article
- Common Bean Coordinated Agricultural Project. [http://www.reeis.usda.gov/web/crisprojectpages/219849.html]
- Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, et al: Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proceedings of the National Academy of Sciences of the United States of America. 2005, 102 (43): 15545-15550. 10.1073/pnas.0506580102.PubMed CentralPubMedView Article
- Joyce AR, Palsson BO: The model organism as a system: integrating 'omics' data sets. NatRevMolCell Biol. 2006, 7 (3): 198-210.
- Choi HK, Mun JH, Kim DJ, Zhu H, Baek JM, Mudge J, Roe B, Ellis N, Doyle J, Kiss GB, et al: Estimating genome conservation between crop and model legume species. ProcNatlAcadSciUSA. 2004, 101 (43): 15289-15294.View Article
- Boutin SR, Young ND, Olson TC, Yu ZH, Vallejos CE, Shoemaker RC: Genome conservation among three legume genera detected with DNA markers. Genome. 1995, 38 (5): 928-937. 10.1139/g95-122.PubMedView Article
- Lee SG, Mykles DL: Proteomics and signal transduction in the crustacean molting gland. Integrative and Comparative Biology. 2006, 46 (6): 965-977. 10.1093/icb/icl047.PubMedView Article
- Oakley BR, Kirsch DR, Morris NR: A Simplified Ultrasensitive Silver Stain for Detecting Proteins in Polyacrylamide Gels. Analytical Biochemistry. 1980, 105 (2): 361-363. 10.1016/0003-2697(80)90470-4.PubMedView Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<url>http://creativecommons.org/licenses/by/2.0</url>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.