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

Genomics Associated Interventions for Heat Stress Tolerance in Cool Season Adapted Grain Legumes

Cool season grain legumes occupy an important place among the agricultural crops and essentially provide multiple benefits including food supply, nutrition security, soil fertility improvement and revenue for farmers all over the world. However, owing to climate change, the average temperature is steadily rising, which negatively affects crop performance and limits their yield. Terminal heat stress that mainly occurred during grain development phases severely harms grain quality and weight in legumes adapted to the cool season, such as lentils, faba beans, chickpeas, field peas, etc. Although, traditional breeding approaches with advanced screening procedures have been employed to identify heat tolerant legume cultivars. Unfortunately, traditional breeding pipelines alone are no longer enough to meet global demands. Genomics-assisted interventions including new-generation sequencing technologies and genotyping platforms have facilitated the development of high-resolution molecular maps, QTL/gene discovery and marker-assisted introgression, thereby improving the efficiency in legumes breeding to develop stress-resilient varieties. Based on the current scenario, we attempted to review the intervention of genomics to decipher different components of tolerance to heat stress and future possibilities of using newly developed genomics-based interventions in cool season adapted grain legumes.

A common bean truncated CRINKLY4 kinase controls gene-for-gene resistance to the fungus Colletotrichum lindemuthianum

Identifying the molecular basis of resistance to pathogens is critical to promote a chemical-free cropping system. In plants, nucleotide-binding leucine-rich repeat constitute the largest family of disease resistance (R) genes, but this resistance can be rapidly overcome by the pathogen, prompting research into alternative sources of resistance. Anthracnose, caused by the fungus Colletotrichum lindemuthianum, is one of the most important diseases of common bean. This study aimed to identify the molecular basis of Co-x, an anthracnose R gene conferring total resistance to the extremely virulent C. lindemuthianum strain 100. To that end, we sequenced the Co-x 58 kb target region in the resistant JaloEEP558 (Co-x) common bean and identified KTR2/3, an additional gene encoding a truncated and chimeric CRINKLY4 kinase, located within a CRINKLY4 kinase cluster. The presence of KTR2/3 is strictly correlated with resistance to strain 100 in a diversity panel of common beans. Furthermore, KTR2/3 expression is up-regulated 24 hours post-inoculation and its transient expression in a susceptible genotype increases resistance to strain 100. Our results provide evidence that Co-x encodes a truncated and chimeric CRINKLY4 kinase probably resulting from an unequal recombination event that occurred recently in the Andean domesticated gene pool. This atypical R gene may act as a decoy involved in indirect recognition of a fungal effector.

Genome-wide transcriptional changes triggered by water deficit on a drought-tolerant common bean cultivar

BACKGROUND

Common bean (Phaseolus vulgaris L.) is a relevant crop cultivated over the world, largely in water insufficiency vulnerable areas. Since drought is the main environmental factor restraining worldwide crop production, efforts have been invested to amend drought tolerance in commercial common bean varieties. However, scarce molecular data are available for those cultivars of P. vulgaris with drought tolerance attributes.

RESULTS

As a first approach, Pinto Saltillo (PS), Azufrado Higuera (AH), and Negro Jamapa Plus (NP) were assessed phenotypically and physiologically to determine the outcome in response to drought on these common bean cultivars. Based on this, a Next-generation sequencing approach was applied to PS, which was the most drought-tolerant cultivar to determine the molecular changes at the transcriptional level. The RNA-Seq analysis revealed that numerous PS genes are dynamically modulated by drought. In brief, 1005 differentially expressed genes (DEGs) were identified, from which 645 genes were up-regulated by drought stress, whereas 360 genes were down-regulated. Further analysis showed that the enriched categories of the up-regulated genes in response to drought fit to processes related to carbohydrate metabolism (polysaccharide metabolic processes), particularly genes encoding proteins located within the cell periphery (cell wall dynamics). In the case of down-regulated genes, heat shock-responsive genes, mainly associated with protein folding, chloroplast, and oxidation-reduction processes were identified.

CONCLUSIONS

Our findings suggest that secondary cell wall (SCW) properties contribute to P. vulgaris L. drought tolerance through alleviation or mitigation of drought-induced osmotic disturbances, making cultivars more adaptable to such stress. Altogether, the knowledge derived from this study is significant for a forthcoming understanding of the molecular mechanisms involved in drought tolerance on common bean, especially for drought-tolerant cultivars such as PS.

Identification and analysis of alternative splicing events in Phaseolus vulgaris and Glycine max

BACKGROUND

The vast diversification of proteins in eukaryotic cells has been related with multiple transcript isoforms from a single gene that result in alternative splicing (AS) of primary transcripts. Analysis of RNA sequencing data from expressed sequence tags and next generation RNA sequencing has been crucial for AS identification and genome-wide AS studies. For the identification of AS events from the related legume species Phaseolus vulgaris and Glycine max, 157 and 88 publicly available RNA-seq libraries, respectively, were analyzed.

RESULTS

We identified 85,570 AS events from P. vulgaris in 72% of expressed genes and 134,316 AS events in 70% of expressed genes from G. max. These were categorized in seven AS event types with intron retention being the most abundant followed by alternative acceptor and alternative donor, representing ~75% of all AS events in both plants. Conservation of AS events in homologous genes between the two species was analyzed where an overrepresentation of AS affecting 5'UTR regions was observed for certain types of AS events. The conservation of AS events was experimentally validated for 8 selected genes, through RT-PCR analysis. The different types of AS events also varied by relative position in the genes. The results were consistent in both species.

CONCLUSIONS

The identification and analysis of AS events are first steps to understand their biological relevance. The results presented here from two related legume species reveal high conservation, over ~15-20 MY of divergence, and may point to the biological relevance of AS.

Differentially expressed genes in mycorrhized and nodulated roots of common bean are associated with defense, cell wall architecture, N metabolism, and P metabolism

Legumes participate in two important endosymbiotic associations, with phosphorus-acquiring arbuscular mycorrhiza (AM, soil fungi) and with nitrogen-fixing bacterial rhizobia. These divergent symbionts share a common symbiotic signal transduction pathway that facilitates the establishment of mycorrhization and nodulation in legumes. However, the unique and shared downstream genes essential for AM and nodule development have not been identified in crop legumes. Here, we used ion torrent next-generation sequencing to perform comparative transcriptomics of common bean (Phaseolus vulgaris) roots colonized by AM or rhizobia. We analyzed global gene expression profiles to identify unique and shared differentially expressed genes (DEGs) that regulate these two symbiotic interactions, and quantitatively compared DEG profiles. We identified 3,219 (1,959 upregulated and 1,260 downregulated) and 2,645 (1,247 upregulated and 1,398 downregulated) unigenes that were differentially expressed in response to mycorrhizal or rhizobial colonization, respectively, compared with uninoculated roots. We obtained quantitative expression profiles of unique and shared genes involved in processes related to defense, cell wall structure, N metabolism, and P metabolism in mycorrhized and nodulated roots. KEGG pathway analysis indicated that most genes involved in jasmonic acid and salicylic acid signaling, N metabolism, and inositol phosphate metabolism are variably expressed during symbiotic interactions. These combined data provide valuable information on symbiotic gene signaling networks that respond to mycorrhizal and rhizobial colonization, and serve as a guide for future genetic strategies to enhance P uptake and N-fixing capacity to increase the net yield of this valuable grain legume.

Natural rhizobial diversity helps to reveal genes and QTLs associated with biological nitrogen fixation in common bean

Common bean is one of the most important crops for human feed, and the most important legume for direct consumption by millions of people, especially in developing countries. It is a promiscuous host legume in terms of nodulation, able to associate with a broad and diverse range of rhizobia, although the competitiveness for nodulation and the nitrogen fixation capacity of most of these strains is generally low. As a result, common bean is very inefficient for symbiotic nitrogen fixation, and nitrogen has to be supplied with chemical fertilizers. In the last years, symbiotic nitrogen fixation has received increasing attention as a sustainable alternative to nitrogen fertilizers, and also as a more economic and available one in poor countries. Therefore, optimization of nitrogen fixation of bean-rhizobia symbioses and selection of efficient rhizobial strains should be a priority, which begins with the study of the natural diversity of the symbioses and the rhizobial populations associated. Natural rhizobia biodiversity that nodulates common bean may be a source of adaptive alleles acting through phenotypic plasticity. Crosses between accessions differing for nitrogen fixation may combine alleles that never meet in nature. Another way to discover adaptive genes is to use association genetics to identify loci that common bean plants use for enhanced biological nitrogen fixation and, in consequence, for marker assisted selection for genetic improvement of symbiotic nitrogen fixation. In this review, rhizobial biodiversity resources will be discussed, together with what is known about the loci that underlie such genetic variation, and the potential candidate genes that may influence the symbiosis' fitness benefits, thus achieving an optimal nitrogen fixation capacity in order to help reduce reliance on nitrogen fertilizers in common bean.

The Class II Trehalose 6-phosphate Synthase Gene PvTPS9 Modulates Trehalose Metabolism in Phaseolus vulgaris Nodules

Legumes form symbioses with rhizobia, producing nitrogen-fixing nodules on the roots of the plant host. The network of plant signaling pathways affecting carbon metabolism may determine the final number of nodules. The trehalose biosynthetic pathway regulates carbon metabolism and plays a fundamental role in plant growth and development, as well as in plant-microbe interactions. The expression of genes for trehalose synthesis during nodule development suggests that this metabolite may play a role in legume-rhizobia symbiosis. In this work, PvTPS9, which encodes a Class II trehalose-6-phosphate synthase (TPS) of common bean (Phaseolus vulgaris), was silenced by RNA interference in transgenic nodules. The silencing of PvTPS9 in root nodules resulted in a reduction of 85% (± 1%) of its transcript, which correlated with a 30% decrease in trehalose contents of transgenic nodules and in untransformed leaves. Composite transgenic plants with PvTPS9 silenced in the roots showed no changes in nodule number and nitrogen fixation, but a severe reduction in plant biomass and altered transcript profiles of all Class II TPS genes. Our data suggest that PvTPS9 plays a key role in modulating trehalose metabolism in the symbiotic nodule and, therefore, in the whole plant.

Genetic Divergence between Camellia sinensis and Its Wild Relatives Revealed via Genome-Wide SNPs from RAD Sequencing

Tea is one of the most popular beverages across the world and is made exclusively from cultivars of Camellia sinensis. Many wild relatives of the genus Camellia that are closely related to C. sinensis are native to Southwest China. In this study, we first identified the distinct genetic divergence between C. sinensis and its wild relatives and provided a glimpse into the artificial selection of tea plants at a genome-wide level by analyzing 15,444 genomic SNPs that were identified from 18 cultivated and wild tea accessions using a high-throughput genome-wide restriction site-associated DNA sequencing (RAD-Seq) approach. Six distinct clusters were detected by phylogeny inferrence and principal component and genetic structural analyses, and these clusters corresponded to six Camellia species/varieties. Genetic divergence apparently indicated that C. taliensis var. bangwei is a semi-wild or transient landrace occupying a phylogenetic position between those wild and cultivated tea plants. Cultivated accessions exhibited greater heterozygosity than wild accessions, with the exception of C. taliensis var. bangwei. Thirteen genes with non-synonymous SNPs exhibited strong selective signals that were suggestive of putative artificial selective footprints for tea plants during domestication. The genome-wide SNPs provide a fundamental data resource for assessing genetic relationships, characterizing complex traits, comparing heterozygosity and analyzing putatitve artificial selection in tea plants.

SNP discovery in common bean by restriction-associated DNA (RAD) sequencing for genetic diversity and population structure analysis

Researchers have made great advances into the development and application of genomic approaches for common beans, creating opportunities to driving more real and applicable strategies for sustainable management of the genetic resource towards plant breeding. This work provides useful polymorphic single-nucleotide polymorphisms (SNPs) for high-throughput common bean genotyping developed by RAD (restriction site-associated DNA) sequencing. The RAD tags were generated from DNA pooled from 12 common bean genotypes, including breeding lines of different gene pools and market classes. The aligned sequences identified 23,748 putative RAD-SNPs, of which 3357 were adequate for genotyping; 1032 RAD-SNPs with the highest ADT (assay design tool) score are presented in this article. The RAD-SNPs were structurally annotated in different coding (47.00 %) and non-coding (53.00 %) sequence components of genes. A subset of 384 RAD-SNPs with broad genome distribution was used to genotype a diverse panel of 95 common bean germplasms and revealed a successful amplification rate of 96.6 %, showing 73 % of polymorphic SNPs within the Andean group and 83 % in the Mesoamerican group. A slightly increased He (0.161, n = 21) value was estimated for the Andean gene pool, compared to the Mesoamerican group (0.156, n = 74). For the linkage disequilibrium (LD) analysis, from a group of 580 SNPs (289 RAD-SNPs and 291 BARC-SNPs) genotyped for the same set of genotypes, 70.2 % were in LD, decreasing to 0.10 %in the Andean group and 0.77 % in the Mesoamerican group. Haplotype patterns spanning 310 Mb of the genome (60 %) were characterized in samples from different origins. However, the haplotype frameworks were under-represented for the Andean (7.85 %) and Mesoamerican (5.55 %) gene pools separately. In conclusion, RAD sequencing allowed the discovery of hundreds of useful SNPs for broad genetic analysis of common bean germplasm. From now, this approach provides an excellent panel of molecular tools for whole genome analysis, allowing integrating and better exploring the common bean breeding practices.

Genome and transcriptome analysis of the Mesoamerican common bean and the role of gene duplications in establishing tissue and temporal specialization of genes

BACKGROUND

Legumes are the third largest family of angiosperms and the second most important crop class. Legume genomes have been shaped by extensive large-scale gene duplications, including an approximately 58 million year old whole genome duplication shared by most crop legumes.

RESULTS

We report the genome and the transcription atlas of coding and non-coding genes of a Mesoamerican genotype of common bean (Phaseolus vulgaris L., BAT93). Using a comprehensive phylogenomics analysis, we assessed the past and recent evolution of common bean, and traced the diversification of patterns of gene expression following duplication. We find that successive rounds of gene duplications in legumes have shaped tissue and developmental expression, leading to increased levels of specialization in larger gene families. We also find that many long non-coding RNAs are preferentially expressed in germ-line-related tissues (pods and seeds), suggesting that they play a significant role in fruit development. Our results also suggest that most bean-specific gene family expansions, including resistance gene clusters, predate the split of the Mesoamerican and Andean gene pools.

CONCLUSIONS

The genome and transcriptome data herein generated for a Mesoamerican genotype represent a counterpart to the genomic resources already available for the Andean gene pool. Altogether, this information will allow the genetic dissection of the characters involved in the domestication and adaptation of the crop, and their further implementation in breeding strategies for this important crop.

Whole genome sequencing and analysis of Swarna, a widely cultivated indica rice variety with low glycemic index

Swarna is a popular cultivated indica rice variety with low glycemic index (GI) but its genetic basis is not known. The whole genome of Swarna was sequenced using Illumina’s paired-end technology, and the reads were mapped to the Nipponbare reference genome. Overall, 65,984 non-synonymous SNPs were identified in 20,350 genes, and in silico analysis predicted that 4,847 of them in 2,214 genes may have deleterious effect on protein functions. Polymorphisms were found in all the starch biosynthesis genes, except the gene for branching enzyme IIa. It was found that T/G SNP at position 246, ‘A’ at position 2,386, and ‘C’ at position 3,378 in the granule bound starch synthase I gene, and C/T SNP at position 1,188 in the glucose-6-phosphate translocator gene may contribute to the low GI phenotype in Swarna. All these variants were also found in the genome of another low GI indica rice variety from Columbia, Fedearroz 50. The whole genome analysis of Swarna helped to understand the genetic basis of GI in rice, which is a complex trait involving multiple factors.

The micro-RNA72c-APETALA2-1 node as a key regulator of the common bean-Rhizobium etli nitrogen fixation symbiosis

Micro-RNAs are recognized as important posttranscriptional regulators in plants. The relevance of micro-RNAs as regulators of the legume-rhizobia nitrogen-fixing symbiosis is emerging. The objective of this work was to functionally characterize the role of micro-RNA172 (miR172) and its conserved target APETALA2 (AP2) transcription factor in the common bean (Phaseolus vulgaris)-Rhizobium etli symbiosis. Our expression analysis revealed that mature miR172c increased upon rhizobial infection and continued increasing during nodule development, reaching its maximum in mature nodules and decaying in senescent nodules. The expression of AP2-1 target showed a negative correlation with miR172c expression. A drastic decrease in miR172c and high AP2-1 mRNA levels were observed in ineffective nodules. Phenotypic analysis of composite bean plants with transgenic roots overexpressing miR172c or a mutated AP2-1 insensitive to miR172c cleavage demonstrated the pivotal regulatory role of the miR172 node in the common bean-rhizobia symbiosis. Increased miR172 resulted in improved root growth, increased rhizobial infection, increased expression of early nodulation and autoregulation of nodulation genes, and improved nodulation and nitrogen fixation. In addition, these plants showed decreased sensitivity to nitrate inhibition of nodulation. Through transcriptome analysis, we identified 114 common bean genes that coexpressed with AP2-1 and proposed these as being targets for transcriptional activation by AP2-1. Several of these genes are related to nodule senescence, and we propose that they have to be silenced, through miR172c-induced AP2-1 cleavage, in active mature nodules. Our work sets the basis for exploring the miR172-mediated improvement of symbiotic nitrogen fixation in common bean, the most important grain legume for human consumption.

Fine-mapping of a major QTL controlling angular leaf spot resistance in common bean (Phaseolus vulgaris L.)

KEY MESSAGE

A major QTL for angular leaf spot resistance in the common bean accession G5686 was fine-mapped to a region containing 36 candidate genes. Markers have been developed for marker-assisted selection. Common bean (Phaseolus vulgaris L.) is an important grain legume and an essential protein source for human nutrition in developing countries. Angular leaf spot (ALS) caused by the pathogen Pseudocercospora griseola (Sacc.) Crous and U. Braun is responsible for severe yield losses of up to 80%. Breeding for resistant cultivars is the most ecological and economical means to control ALS and is particularly important for yield stability in low-input agriculture. Here, we report on a fine-mapping approach of a major quantitative trait locus (QTL) ALS4.1(GS, UC) for ALS resistance in a mapping population derived from the resistant genotype G5686 and the susceptible cultivar Sprite. 180 F3 individuals of the mapping population were evaluated for ALS resistance and genotyped with 22 markers distributed over 11 genome regions colocating with previously reported QTL for ALS resistance. Multiple QTL analysis identified three QTL regions, including one major QTL on chromosome Pv04 at 43.7 Mbp explaining over 75% of the observed variation for ALS resistance. Additional evaluation of 153 F4, 89 BC1F2 and 139 F4/F5/BC1F3 descendants with markers in the region of the major QTL delimited the region to 418 kbp harboring 36 candidate genes. Among these, 11 serine/threonine protein kinases arranged in a repetitive array constitute promising candidate genes for controlling ALS resistance. Single nucleotide polymorphism markers cosegregating with the major QTL for ALS resistance have been developed and constitute the basis for marker-assisted introgression of ALS resistance into advanced breeding germplasm of common bean.

Genome-wide identification of SSR and SNP markers from the non-heading Chinese cabbage for comparative genomic analyses

BACKGROUND

Non-heading Chinese cabbage (NHCC), belonging to Brassica, is an important leaf vegetable in Asia. Although genetic analyses have been performed through conventional selection and breeding efforts, the domestication history of NHCC and the genetics underlying its morphological diversity remain unclear. Thus, the reliable molecular markers representative of the whole genome are required for molecular-assisted selection in NHCC.

RESULTS

A total of 20,836 simple sequence repeats (SSRs) were detected in NHCC, containing repeat types from mononucleotide to nonanucleotide. The average density was 62.93 SSRs/Mb. In gene regions, 5,435 SSRs were identified in 4,569 genes. A total of 5,008 primer pairs were designed, and 74 were randomly selected for validation. Among these, 60 (81.08%) were polymorphic in 18 Cruciferae. The number of polymorphic bands ranged from two to five, with an average of 2.70 for each primer. The average values of the polymorphism information content, observed heterozygosity, Hardy-Weinberg equilibrium, and Shannon's information index were 0.2970, 0.4136, 0.5706, and 0.5885, respectively. Four clusters were classified according to the unweighted pair-group method with arithmetic average cluster analysis of 18 genotypes. In addition, a total of 1,228,979 single nucleotide polymorphisms (SNPs) were identified in the NHCC through a comparison with the genome of Chinese cabbage, and the average SNP density in the whole genome was 4.33/Kb. The number of SNPs ranged from 341,939 to 591,586 in the 10 accessions, and the average heterozygous SNPs ratio was ~42.53%. All analyses showed these markers were high quality and reliable. Therefore, they could be used in the construction of a linkage map and for genetic diversity studies for NHCC in future.

CONCLUSIONS

This is the first systematic and comprehensive analysis and identification of SSRs in NHCC and 17 species. The development of a large number of SNP and SSR markers was successfully achieved for NHCC. These novel markers are valuable for constructing genetic linkage maps, comparative genome analysis, quantitative trait locus (QTL) mapping, genome-wide association studies, and marker-assisted selection in NHCC breeding system research.

Analysis of BAC-end sequences in common bean (Phaseolus vulgaris L.) towards the development and characterization of long motifs SSRs

The increasing volume of genomic data on the Phaseolus vulgaris species have contributed to its importance as a model genetic species and positively affected the investigation of other legumes of scientific and economic value. To expand and gain a more in-depth knowledge of the common bean genome, the ends of a number of bacterial artificial chromosome (BAC) were sequenced, annotated and the presence of repetitive sequences was determined. In total, 52,270 BESs (BAC-end sequences), equivalent to 32 Mbp (~6 %) of the genome, were processed. In total, 3,789 BES-SSRs were identified, with a distribution of one SSR (simple sequence repeat) per 8.36 kbp and 2,000 were suitable for the development of SSRs, of which 194 were evaluated in low-resolution screening. From 40 BES-SSRs based on long motifs SSRs (≥ trinucleotides) analyzed in high-resolution genotyping, 34 showed an equally good amplification for the Andean and for the Mesoamerican genepools, exhibiting an average gene diversity (H E) of 0.490 and 5.59 alleles/locus, of which six classified as Class I showed a H E ≥ 0.7. The PCoA and structure analysis allowed to discriminate the gene pools (K = 2, FST = 0.733). From the 52,270 BESs, 2 % corresponded to transcription factors and 3 % to transposable elements. Putative functions for 24,321 BESs were identified and for 19,363 were assigned functional categories (gene ontology). This study identified highly polymorphic BES-SSRs containing tri- to hexanucleotides motifs and bringing together relevant genetic characteristics useful for breeding programs. Additionally, the BESs were incorporated into the international genome-sequencing project for the common bean.

An RNA-Seq based gene expression atlas of the common bean

BACKGROUND

Common bean (Phaseolus vulgaris) is grown throughout the world and comprises roughly 50% of the grain legumes consumed worldwide. Despite this, genetic resources for common beans have been lacking. Next generation sequencing, has facilitated our investigation of the gene expression profiles associated with biologically important traits in common bean. An increased understanding of gene expression in common bean will improve our understanding of gene expression patterns in other legume species.

RESULTS

Combining recently developed genomic resources for Phaseolus vulgaris, including predicted gene calls, with RNA-Seq technology, we measured the gene expression patterns from 24 samples collected from seven tissues at developmentally important stages and from three nitrogen treatments. Gene expression patterns throughout the plant were analyzed to better understand changes due to nodulation, seed development, and nitrogen utilization. We have identified 11,010 genes differentially expressed with a fold change ≥ 2 and a P-value < 0.05 between different tissues at the same time point, 15,752 genes differentially expressed within a tissue due to changes in development, and 2,315 genes expressed only in a single tissue. These analyses identified 2,970 genes with expression patterns that appear to be directly dependent on the source of available nitrogen. Finally, we have assembled this data in a publicly available database, The Phaseolus vulgaris Gene Expression Atlas (Pv GEA), http://plantgrn.noble.org/PvGEA/ . Using the website, researchers can query gene expression profiles of their gene of interest, search for genes expressed in different tissues, or download the dataset in a tabular form.

CONCLUSIONS

These data provide the basis for a gene expression atlas, which will facilitate functional genomic studies in common bean. Analysis of this dataset has identified genes important in regulating seed composition and has increased our understanding of nodulation and impact of the nitrogen source on assimilation and distribution throughout the plant.

Fine mapping of Co-x, an anthracnose resistance gene to a highly virulent strain of Colletotrichum lindemuthianum in common bean

The Co - x anthracnose R gene of common bean was fine-mapped into a 58 kb region at one end of chromosome 1, where no canonical NB-LRR-encoding genes are present in G19833 genome sequence. Anthracnose, caused by the phytopathogenic fungus Colletotrichum lindemuthianum, is one of the most damaging diseases of common bean, Phaseolus vulgaris. Various resistance (R) genes, named Co-, conferring race-specific resistance to different strains of C. lindemuthianum have been identified. The Andean cultivar JaloEEP558 was reported to carry Co-x on chromosome 1, conferring resistance to the highly virulent strain 100. To fine map Co-x, 181 recombinant inbred lines derived from the cross between JaloEEP558 and BAT93 were genotyped with polymerase chain reaction (PCR)-based markers developed using the genome sequence of the Andean genotype G19833. Analysis of RILs carrying key recombination events positioned Co-x at one end of chromosome 1 to a 58 kb region of the G19833 genome sequence. Annotation of this target region revealed eight genes: three phosphoinositide-specific phospholipases C (PI-PLC), one zinc finger protein and four kinases, suggesting that Co-x is not a classical nucleotide-binding leucine-rich encoding gene. In addition, we identified and characterized the seven members of common bean PI-PLC gene family distributed into two clusters located at the ends of chromosomes 1 and 8. Co-x is not a member of Co-1 allelic series since these two genes are separated by at least 190 kb. Comparative analysis between soybean and common bean revealed that the Co-x syntenic region, located at one end of Glycine max chromosome 18, carries Rhg1, a major QTL contributing to soybean cyst nematode resistance. The PCR-based markers generated in this study should be useful in marker-assisted selection for pyramiding Co-x with other R genes.

Transcriptome analysis of salt tolerant common bean (Phaseolus vulgaris L.) under saline conditions

Salinity is one of the important abiotic stress factors that limit crop production. Common bean, Phaseolus vulgaris L., a major protein source in developing countries, is highly affected by soil salinity and the information on genes that play a role in salt tolerance is scarce. We aimed to identify differentially expressed genes (DEGs) and related pathways by comprehensive analysis of transcriptomes of both root and leaf tissues of the tolerant genotype grown under saline and control conditions in hydroponic system. We have generated a total of 158 million high-quality reads which were assembled into 83,774 all-unigenes with a mean length of 813 bp and N50 of 1,449 bp. Among the all-unigenes, 58,171 were assigned with Nr annotations after homology analyses. It was revealed that 6,422 and 4,555 all-unigenes were differentially expressed upon salt stress in leaf and root tissues respectively. Validation of the RNA-seq quantifications (RPKM values) was performed by qRT-PCR (Quantitative Reverse Transcription PCR) analyses. Enrichment analyses of DEGs based on GO and KEGG databases have shown that both leaf and root tissues regulate energy metabolism, transmembrane transport activity, and secondary metabolites to cope with salinity. A total of 2,678 putative common bean transcription factors were identified and classified under 59 transcription factor families; among them 441 were salt responsive. The data generated in this study will help in understanding the fundamentals of salt tolerance in common bean and will provide resources for functional genomic studies.

Nucleotide diversity based on phaseolin and iron reductase genes in common bean accessions of different geographical origins

Discriminating genotypes within plant collections is imperative, and DNA sequence approaches for detecting single nucleotide polymorphisms (SNPs) have proved essential in any modern analysis of germplasm. By sequencing the α-Phs and PvFRO1 genes that, respectively, encode phaseolin and an iron reductase, we prospected for SNPs in exonic and intronic regions of both genes in a sample of 31 accessions of Phaseolus vulgaris from Mesoamerican and Andean gene pools, and one accession of Phaseolus lunatus, chosen as an outgroup. Sequence alignment showed 95 SNPs in α-Phs and 83 in PvFRO1, but diversity along the nucleotide sequences was not evenly distributed in both genes. Accessions from the same gene pool showed greater similarity than those from different gene pools, and the cluster patterns obtained in this study were consistent with the hierarchical organization into two P. vulgaris gene pools. The polymorphisms detected in the α-Phs gene allowed better discrimination among the accessions within each cluster than the PvFRO1 polymorphisms. Furthermore, some variations within exons changes amino acids in both predicted protein sequences. In an unprecedented result, the phaseolin-predicted amino acid variation allowed most of the accessions to be typified.

Identification and characterization of nucleotide variations in the genome of Ziziphus jujuba (Rhamnaceae) by next generation sequencing

In this study, single-nucleotide polymorphisms (SNPs) and insertions/deletions (InDels) in the genome of Ziziphus jujuba were identified using sequences generated by the Roche 454 GS-FLX sequencer. A total of, 573,141 reads were produced with an average read length of 360 bp. After quality control, 258,754 of the filtered reads were assembled into 23,864 contigs, and 293,458 remained as singletons. Using the contig assemblies as a reference, 17,160 SNPs and 478 InDels were identified. Among the SNPs, transitions occurred three times more frequently than transversions. In transitions, the number of C/T and G/A transitions was similar. Among the transversions, A/T was the most abundant, and C/G was much rarer than any of the other types of transversions, accounting for only about half the numbers of A/C, A/T and G/T transversions. For the InDels, mononucleotide changes amounted to 64.4% of the total number of InDels. In general, the frequency of detected InDels decreased as the length of the InDels increased. This study provides valuable marker resources for future genetic studies of Ziziphus spp.

Transcript profiling of common bean nodules subjected to oxidative stress

Several environmental stresses generate high amounts of reactive oxygen species (ROS) in plant cells, resulting in oxidative stress. Symbiotic nitrogen fixation (SNF) in the legume-rhizobia symbiosis is sensitive to damage from oxidative stress. Active nodules of the common bean (Phaseolus vulgaris) exposed to the herbicide paraquat (1,1'-dimethyl-4,4'-bipyridinium dichloride hydrate), which stimulates ROS accumulation, exhibited reduced nitrogenase activity and ureide content. We analyzed the global gene response of nodules subjected to oxidative stress using the Bean Custom Array 90K, which includes probes from 30,000 expressed sequence tags (ESTs). A total of 4280 ESTs were differentially expressed in stressed bean nodules; of these, 2218 were repressed. Based on Gene Ontology analysis, these genes were grouped into 42 different biological process categories. Analysis with the PathExpress bioinformatic tool, adapted for bean, identified five significantly repressed metabolic pathways related to carbon/nitrogen metabolism, which is crucial for nodule function. Quantitative reverse transcription (qRT)-PCR analysis of transcription factor (TF) gene expression showed that 67 TF genes were differentially expressed in nodules exposed to oxidative stress. Putative cis-elements recognized by highly responsive TF were detected in promoter regions of oxidative stress regulated genes. The expression of oxidative stress responsive genes and of genes important for SNF in bacteroids analyzed in stressed nodules revealed that these conditions elicited a transcriptional response.

A comparison of the molecular organization of genomic regions associated with resistance to common bacterial blight in two Phaseolus vulgaris genotypes

Resistance to common bacterial blight, caused by Xanthomonas axonopodis pv. phaseoli, in Phaseolus vulgaris is conditioned by several loci on different chromosomes. Previous studies with OAC-Rex, a CBB-resistant, white bean variety of Mesoamerican origin, identified two resistance loci associated with the molecular markers Pv-CTT001 and SU91, on chromosome 4 and 8, respectively. Resistance to CBB is assumed to be derived from an interspecific cross with Phaseolus acutifolius in the pedigree of OAC-Rex. Our current whole genome sequencing effort with OAC-Rex provided the opportunity to compare its genome in the regions associated with CBB resistance with the v1.0 release of the P. vulgaris line G19833, which is a large seeded bean of Andean origin, and (assumed to be) CBB susceptible. In addition, the genomic regions containing SAP6, a marker associated with P. vulgaris-derived CBB-resistance on chromosome 10, were compared. These analyses indicated that gene content was highly conserved between G19833 and OAC-Rex across the regions examined (>80%). However, fifty-nine genes unique to OAC Rex were identified, with resistance gene homologues making up the largest category (10 genes identified). Two unique genes in OAC-Rex located within the SU91 resistance QTL have homology to P. acutifolius ESTs and may be potential sources of CBB resistance. As the genomic sequence assembly of OAC-Rex is completed, we expect that further comparisons between it and the G19833 genome will lead to a greater understanding of CBB resistance in bean.

Phenotyping common beans for adaptation to drought

Common beans (Phaseolus vulgaris L.) originated in the New World and are the grain legume of greatest production for direct human consumption. Common bean production is subject to frequent droughts in highland Mexico, in the Pacific coast of Central America, in northeast Brazil, and in eastern and southern Africa from Ethiopia to South Africa. This article reviews efforts to improve common bean for drought tolerance, referring to genetic diversity for drought response, the physiology of drought tolerance mechanisms, and breeding strategies. Different races of common bean respond differently to drought, with race Durango of highland Mexico being a major source of genes. Sister species of P. vulgaris likewise have unique traits, especially P. acutifolius which is well adapted to dryland conditions. Diverse sources of tolerance may have different mechanisms of plant response, implying the need for different methods of phenotyping to recognize the relevant traits. Practical considerations of field management are discussed including: trial planning; water management; and field preparation.

Medicago truncatula DNF2 is a PI-PLC-XD-containing protein required for bacteroid persistence and prevention of nodule early senescence and defense-like reactions

Medicago truncatula and Sinorhizobium meliloti form a symbiotic association resulting in the formation of nitrogen-fixing nodules. Nodule cells contain large numbers of bacteroids which are differentiated, nitrogen-fixing forms of the symbiotic bacteria. In the nodules, symbiotic plant cells home and maintain hundreds of viable bacteria. In order to better understand the molecular mechanism sustaining the phenomenon, we searched for new plant genes required for effective symbiosis. We used a combination of forward and reverse genetics approaches to identify a gene required for nitrogen fixation, and we used cell and molecular biology to characterize the mutant phenotype and to gain an insight into gene function. The symbiotic gene DNF2 encodes a putative phosphatidylinositol phospholipase C-like protein. Nodules formed by the mutant contain a zone of infected cells reduced to a few cell layers. In this zone, bacteria do not differentiate properly into bacteroids. Furthermore, mutant nodules senesce rapidly and exhibit defense-like reactions. This atypical phenotype amongst Fix(-) mutants unravels dnf2 as a new actor of bacteroid persistence inside symbiotic plant cells.

Isolation and characterization of nucleotide-binding site resistance gene homologues in common bean (Phaseolus vulgaris)

Common bean production is constrained by many fungal, viral, and bacterial pathogens. Thus, the identification of resistance (R) genes is an important focal point of common bean research. The main goal of our study was to identify resistance gene homologues (RGH) in the crop, using degenerate primers designed from conserved sequences in the nucleotide-binding site (NBS) domains of R-genes from the model legume Medicago truncatula. Total DNA of the Andean common bean genotype G19833 was used for amplification of over 500 primer combinations. Sequencing of amplicons showed that 403 cloned fragments had uninterrupted open reading frames and were considered representative of functional RGH genes. The sequences were grouped at two levels of nucleotide identity (90 and 80%) and representative sequences of each group were used for phylogenetic analyses. The RGH sequence diversity of common bean was divided into TIR and non-TIR families, each with different clusters. The TIR sequences grouped into 14 clades while non-TIR sequences grouped into seven clades. Pairwise comparisons showed purifying selection, although some sequences may have been the result of diversifying selection. Knowledge about RGH genes in common bean can allow the design of molecular markers for pyramiding of resistance genes against various pathogens.

Down-regulation of PvTRE1 enhances nodule biomass and bacteroid number in the common bean

Legume-rhizobium interactions have been widely studied and characterized, and the disaccharide trehalose has been commonly detected during this symbiotic interaction. It has been proposed that trehalose content in nodules during this symbiotic interaction might be regulated by trehalase. In the present study, we assessed the role of trehalose accumulation by down-regulating trehalase in the nodules of common bean plants. We performed gene expression analysis for trehalase (PvTRE1) during nodule development. PvTRE1 was knocked down by RNA interference (RNAi) in transgenic nodules of the common bean. PvTRE1 expression in nodulated roots is mainly restricted to nodules. Down-regulation of PvTRE1 led to increased trehalose content (78%) and bacteroid number (almost one order of magnitude). In addition, nodule biomass, nitrogenase activity, and GOGAT transcript accumulation were significantly enhanced too. The trehalose accumulation, triggered by PvTRE1 down-regulation, led to a positive impact on the legume-rhizobium symbiotic interaction. This could contribute to the agronomical enhancement of symbiotic nitrogen fixation.

Dissecting Phaseolus vulgaris innate immune system against Colletotrichum lindemuthianum infection

Background

The genus Colletotrichum is one of the most economically important plant pathogens, causing anthracnose on a wide range of crops including common beans (Phaseolus vulgaris L.). Crop yield can be dramatically decreased depending on the plant cultivar used and the environmental conditions. This study aimed to identify potential genetic components of the bean immune system to provide environmentally friendly control measures against this fungus.

Methodology and Principal Findings

As the common bean is not amenable to reverse genetics to explore functionality and its genome is not fully curated, we used putative Arabidopsis orthologs of bean expressed sequence tag (EST) to perform bioinformatic analysis and experimental validation of gene expression to identify common bean genes regulated during the incompatible interaction with C. lindemuthianum. Similar to model pathosystems, Gene Ontology (GO) analysis indicated that hormone biosynthesis and signaling in common beans seem to be modulated by fungus infection. For instance, cytokinin and ethylene responses were up-regulated and jasmonic acid, gibberellin, and abscisic acid responses were down-regulated, indicating that these hormones may play a central role in this pathosystem. Importantly, we have identified putative bean gene orthologs of Arabidopsis genes involved in the plant immune system. Based on experimental validation of gene expression, we propose that hypersensitive reaction as part of effector-triggered immunity may operate, at least in part, by down-regulating genes, such as FLS2-like and MKK5-like, putative orthologs of the Arabidopsis genes involved in pathogen perception and downstream signaling.

Conclusions/Significance

We have identified specific bean genes and uncovered metabolic processes and pathways that may be involved in the immune response against pathogens. Our transcriptome database is a rich resource for mining novel defense-related genes, which enabled us to develop a model of the molecular components of the bean innate immune system regulated upon pathogen attack.

Genome-wide DNA polymorphisms in elite indica rice inbreds discovered by whole-genome sequencing

Advances in next-generation sequencing technologies have aided discovery of millions of genome-wide DNA polymorphisms, single nucleotide polymorphisms (SNPs) and insertions-deletions (InDels), which are an invaluable resource for marker-assisted breeding. Whole-genome resequencing of six elite indica rice inbreds (three cytoplasmic male sterile and three restorer lines) resulted in the generation of 338 million 75-bp paired-end reads, which provided 85.4% coverage of the Nipponbare genome. A total of 2 819 086 nonredundant DNA polymorphisms including 2 495 052 SNPs, 160 478 insertions and 163 556 deletions were discovered between the inbreds and Nipponbare, providing an average of 6.8 SNPs/kb across the genome. Distribution of SNPs and InDels in the chromosome was nonrandom with SNP-rich and SNP-poor regions being evident across the genome. A contiguous 4.3-Mb region on chromosome 5 with extremely low SNP density was identified. Overall, 83 262 nonsynonymous SNPs spanning 16 379 genes and 3620 nonsynonymous InDels in 2625 genes have been discovered which provide valuable insights into the basis underlying performance of the inbreds and the hybrids between these inbred combinations. SNPs and InDels discovered from this diverse set of indica rice inbreds not only enrich SNP resources for molecular breeding but also enable the study of genome-wide variations on hybrid performance.

Gene-based single nucleotide polymorphism markers for genetic and association mapping in common bean

Background

In common bean, expressed sequence tags (ESTs) are an underestimated source of gene-based markers such as insertion-deletions (Indels) or single-nucleotide polymorphisms (SNPs). However, due to the nature of these conserved sequences, detection of markers is difficult and portrays low levels of polymorphism. Therefore, development of intron-spanning EST-SNP markers can be a valuable resource for genetic experiments such as genetic mapping and association studies.

Results

In this study, a total of 313 new gene-based markers were developed at target genes. Intronic variation was deeply explored in order to capture more polymorphism. Introns were putatively identified after comparing the common bean ESTs with the soybean genome, and the primers were designed over intron-flanking regions. The intronic regions were evaluated for parental polymorphisms using the single strand conformational polymorphism (SSCP) technique and Sequenom MassARRAY system. A total of 53 new marker loci were placed on an integrated molecular map in the DOR364 × G19833 recombinant inbred line (RIL) population. The new linkage map was used to build a consensus map, merging the linkage maps of the BAT93 × JALO EEP558 and DOR364 × BAT477 populations. A total of 1,060 markers were mapped, with a total map length of 2,041 cM across 11 linkage groups. As a second application of the generated resource, a diversity panel with 93 genotypes was evaluated with 173 SNP markers using the MassARRAY-platform and KASPar technology. These results were coupled with previous SSR evaluations and drought tolerance assays carried out on the same individuals. This agglomerative dataset was examined, in order to discover marker-trait associations, using general linear model (GLM) and mixed linear model (MLM). Some significant associations with yield components were identified, and were consistent with previous findings.

Conclusions

In short, this study illustrates the power of intron-based markers for linkage and association mapping in common bean. The utility of these markers is discussed in relation with the usefulness of microsatellites, the molecular markers by excellence in this crop.

Characterization of two putative protein phosphatase genes and their involvement in phosphorus efficiency in Phaseolus vulgaris

Protein dephosphorylation mediated by protein phosphatases plays a major role in signal transduction of plant responses to environmental stresses. In this study, two putative protein phosphatases, PvPS2:1 and PvPS2:2 were identified and characterized in bean (Phaseolus vulgaris). The two PvPS2 members were found to be localized to the plasma membrane and the nucleus by transient expression of PvPS2:GFP in onion epidermal cells. Transcripts of the two PvPS2 genes were significantly increased by phosphate (P(i) ) starvation in the two bean genotypes, G19833 (a P-efficient genotype) and DOR364 (a P-inefficient genotype). However, G19833 exhibited higher PvPS2:1 expression levels than DOR364 in both leaves and roots during P(i) starvation. Increased transcription of PvPS2:1 in response to P(i) starvation was further verified through histochemical analysis of PvPS2:1 promoter fusion ß-glucuronidase (GUS) in transgenic Arabidopsis plants. Analysis of PvPS2:1 overexpression lines in bean hairy roots and Arabidopsis showed that PvS2:1 was involved in root growth and P accumulation. Furthermore, expression levels of two P(i) starvation responsive genes were upregulated and the APase activities were enhanced in the overexpressing PvPS2:1 Arabidopsis lines. Taken together, our results strongly suggested that PvPS2:1 positively regulated plant responses to P(i) starvation, and could be further targeted as a candidate gene to improve crop P efficiency.

Comparative analysis of PvPAP gene family and their functions in response to phosphorus deficiency in common bean

BACKGROUND

Purple acid phosphatases (PAPs) play a vital role in adaptive strategies of plants to phosphorus (P) deficiency. However, their functions in relation to P efficiency are fragmentary in common bean.

PRINCIPAL FINDINGS

Five PvPAPs were isolated and sequenced in common bean. Phylogenetic analysis showed that PvPAPs could be classified into two groups, including a small group with low molecular mass, and a large group with high molecular mass. Among them, PvPAP3, PvPAP4 and PvPAP5 belong to the small group, while the other two belong to the large group. Transient expression of 35S:PvPAPs-GFP on onion epidermal cells verified the variations of subcellular localization among PvPAPs, suggesting functional diversities of PvPAPs in common bean. Quantitative PCR results showed that most PvPAPs were up-regulated by phosphate (Pi) starvation. Among them, the expression of the small group PvPAPs responded more to Pi starvation, especially in the roots of G19833, the P-efficient genotype. However, only overexpressing PvPAP1 and PvPAP3 could result in significantly increased utilization of extracellular dNTPs in the transgenic bean hairy roots. Furthermore, overexpressing PvPAP3 in Arabidopsis enhanced both plant growth and total P content when dNTPs were supplied as the sole external P source.

CONCLUSIONS

The results suggest that PvPAPs in bean varied in protein structure, response to P deficiency and subcellular localization. Among them, both PvPAP1 and PvPAP3 might function as utilization of extracellular dNTPs.

Four allantoinase genes are expressed in nitrogen-fixing soybean

Soybean (Glycine max L. [Merr]) plants export nitrogen from the nodules as ureides during symbiotic biological nitrogen fixation. Ureides also play a role as nitrogen storage compounds in the seeds and are broken down in germinating seedlings. In this work we identified four soybean genes encoding allantoinase (E.C. 3.5.2.5), an enzyme involved in both ureide production in nodules and ureide catabolism in leaves and other sink tissues. We examined ureide content, allantoinase enzyme activity and expression of these genes, which we term GmALN1 through GmALN4, in germinating seedlings and in vegetative tissues from 45 day old soybean plants. GmALN1 and GmALN2 transcripts were measured in all tissues, but similar levels of expression of GmALN3 and GmALN4 was only observed in nodules. The soybean allantoinase genes seem to have arisen through tandem gene duplication followed by a whole genome duplication. We looked for evidence of the tandem duplication in common bean (Phaseolus vulgaris L.) and present evidence that it occured sometime in the bean lineage before these two species diverged, but before soybean became a tetraploid.

Construction and EST sequencing of full-length, drought stress cDNA libraries for common beans (Phaseolus vulgaris L.)

BACKGROUND

Common bean is an important legume crop with only a moderate number of short expressed sequence tags (ESTs) made with traditional methods. The goal of this research was to use full-length cDNA technology to develop ESTs that would overlap with the beginning of open reading frames and therefore be useful for gene annotation of genomic sequences. The library was also constructed to represent genes expressed under drought, low soil phosphorus and high soil aluminum toxicity. We also undertook comparisons of the full-length cDNA library to two previous non-full clone EST sets for common bean.

RESULTS

Two full-length cDNA libraries were constructed: one for the drought tolerant Mesoamerican genotype BAT477 and the other one for the acid-soil tolerant Andean genotype G19833 which has been selected for genome sequencing. Plants were grown in three soil types using deep rooting cylinders subjected to drought and non-drought stress and tissues were collected from both roots and above ground parts. A total of 20,000 clones were selected robotically, half from each library. Then, nearly 10,000 clones from the G19833 library were sequenced with an average read length of 850 nucleotides. A total of 4,219 unigenes were identified consisting of 2,981 contigs and 1,238 singletons. These were functionally annotated with gene ontology terms and placed into KEGG pathways. Compared to other EST sequencing efforts in common bean, about half of the sequences were novel or represented the 5' ends of known genes.

CONCLUSIONS

The present full-length cDNA libraries add to the technological toolbox available for common bean and our sequencing of these clones substantially increases the number of unique EST sequences available for the common bean genome. All of this should be useful for both functional gene annotation, analysis of splice site variants and intron/exon boundary determination by comparison to soybean genes or with common bean whole-genome sequences. In addition the library has a large number of transcription factors and will be interesting for discovery and validation of drought or abiotic stress related genes in common bean.

Identification and analysis of common bean (Phaseolus vulgaris L.) transcriptomes by massively parallel pyrosequencing

BACKGROUND

Common bean (Phaseolus vulgaris) is the most important food legume in the world. Although this crop is very important to both the developed and developing world as a means of dietary protein supply, resources available in common bean are limited. Global transcriptome analysis is important to better understand gene expression, genetic variation, and gene structure annotation in addition to other important features. However, the number and description of common bean sequences are very limited, which greatly inhibits genome and transcriptome research. Here we used 454 pyrosequencing to obtain a substantial transcriptome dataset for common bean.

RESULTS

We obtained 1,692,972 reads with an average read length of 207 nucleotides (nt). These reads were assembled into 59,295 unigenes including 39,572 contigs and 19,723 singletons, in addition to 35,328 singletons less than 100 bp. Comparing the unigenes to common bean ESTs deposited in GenBank, we found that 53.40% or 31,664 of these unigenes had no matches to this dataset and can be considered as new common bean transcripts. Functional annotation of the unigenes carried out by Gene Ontology assignments from hits to Arabidopsis and soybean indicated coverage of a broad range of GO categories. The common bean unigenes were also compared to the bean bacterial artificial chromosome (BAC) end sequences, and a total of 21% of the unigenes (12,724) including 9,199 contigs and 3,256 singletons match to the 8,823 BAC-end sequences. In addition, a large number of simple sequence repeats (SSRs) and transcription factors were also identified in this study.

CONCLUSIONS

This work provides the first large scale identification of the common bean transcriptome derived by 454 pyrosequencing. This research has resulted in a 150% increase in the number of Phaseolus vulgaris ESTs. The dataset obtained through this analysis will provide a platform for functional genomics in common bean and related legumes and will aid in the development of molecular markers that can be used for tagging genes of interest. Additionally, these sequences will provide a means for better annotation of the on-going common bean whole genome sequencing.

Nodulin 41, a novel late nodulin of common bean with peptidase activity

BACKGROUND

The legume-rhizobium symbiosis requires the formation of root nodules, specialized organs where the nitrogen fixation process takes place. Nodule development is accompanied by the induction of specific plant genes, referred to as nodulin genes. Important roles in processes such as morphogenesis and metabolism have been assigned to nodulins during the legume-rhizobium symbiosis.

RESULTS

Here we report the purification and biochemical characterization of a novel nodulin from common bean (Phaseolus vulgaris L.) root nodules. This protein, called nodulin 41 (PvNod41) was purified through affinity chromatography and was partially sequenced. A genomic clone was then isolated via PCR amplification. PvNod41 is an atypical aspartyl peptidase of the A1B subfamily with an optimal hydrolytic activity at pH 4.5. We demonstrate that PvNod41 has limited peptidase activity against casein and is partially inhibited by pepstatin A. A PvNod41-specific antiserum was used to assess the expression pattern of this protein in different plant organs and throughout root nodule development, revealing that PvNod41 is found only in bean root nodules and is confined to uninfected cells.

CONCLUSIONS

To date, only a small number of atypical aspartyl peptidases have been characterized in plants. Their particular spatial and temporal expression patterns along with their unique enzymatic properties imply a high degree of functional specialization. Indeed, PvNod41 is closely related to CDR1, an Arabidopsis thaliana extracellular aspartyl protease involved in defense against bacterial pathogens. PvNod41's biochemical properties and specific cell-type localization, in uninfected cells of the common bean root nodule, strongly suggest that this aspartyl peptidase has a key role in plant defense during the symbiotic interaction.

SNP marker diversity in common bean (Phaseolus vulgaris L.)

Single nucleotide polymorphism (SNP) markers have become a genetic technology of choice because of their automation and high precision of allele calls. In this study, our goal was to develop 94 SNPs and test them across well-chosen common bean (Phaseolus vulgaris L.) germplasm. We validated and accessed SNP diversity at 84 gene-based and 10 non-genic loci using KASPar technology in a panel of 70 genotypes that have been used as parents of mapping populations and have been previously evaluated for SSRs. SNPs exhibited high levels of genetic diversity, an excess of middle frequency polymorphism, and a within-genepool mismatch distribution as expected for populations affected by sudden demographic expansions after domestication bottlenecks. This set of markers was useful for distinguishing Andean and Mesoamerican genotypes but less useful for distinguishing within each gene pool. In summary, slightly greater polymorphism and race structure was found within the Andean gene pool than within the Mesoamerican gene pool but polymorphism rate between genotypes was consistent with genepool and race identity. Our survey results represent a baseline for the choice of SNP markers for future applications because gene-associated SNPs could themselves be causative SNPs for traits. Finally, we discuss that the ideal genetic marker combination with which to carry out diversity, mapping and association studies in common bean should consider a mix of both SNP and SSR markers.

Chromosome identification in the Andean common bean accession G19833 (Phaseolus vulgaris L., Fabaceae)

Characterization of all chromosomes of the Andean G19833 bean genotype was carried out by fluorescent in situ hybridization. Eleven single-copy genomic sequences, one for each chromosome, two BACs containing subtelomeric and pericentromeric repeats and the 5S and 45S ribosomal DNA (rDNA) were used as probes. Comparison to the Mesoamerican accession BAT93 showed little divergence, except for additional 45S rDNA sites in four chromosome pairs. Altogether, the results indicated a relative karyotypic stability during the evolution of the Andean and Mesoamerican gene pools of P. vulgaris.

Identification of SNP and SSR markers in eggplant using RAD tag sequencing

Background

The eggplant (Solanum melongena L.) genome is relatively unexplored, especially compared to those of the other major Solanaceae crops tomato and potato. In particular, no SNP markers are publicly available; on the other hand, over 1,000 SSR markers were developed and publicly available. We have combined the recently developed Restriction-site Associated DNA (RAD) approach with Illumina DNA sequencing for rapid and mass discovery of both SNP and SSR markers for eggplant.

Results

RAD tags were generated from the genomic DNA of a pair of eggplant mapping parents, and sequenced to produce ~17.5 Mb of sequences arrangeable into ~78,000 contigs. The resulting non-redundant genomic sequence dataset consisted of ~45,000 sequences, of which ~29% were putative coding sequences and ~70% were in common between the mapping parents. The shared sequences allowed the discovery of ~10,000 SNPs and nearly 1,000 indels, equivalent to a SNP frequency of 0.8 per Kb and an indel frequency of 0.07 per Kb. Over 2,000 of the SNPs are likely to be mappable via the Illumina GoldenGate assay. A subset of 384 SNPs was used to successfully fingerprint a panel of eggplant germplasm, producing a set of informative diversity data. The RAD sequences also included nearly 2,000 putative SSRs, and primer pairs were designed to amplify 1,155 loci.

Conclusion

The high throughput sequencing of the RAD tags allowed the discovery of a large number of DNA markers, which will prove useful for extending our current knowledge of the genome organization of eggplant, for assisting in marker-aided selection and for carrying out comparative genomic analyses within the Solanaceae family.

Analysis of common bean expressed sequence tags identifies sulfur metabolic pathways active in seed and sulfur-rich proteins highly expressed in the absence of phaseolin and major lectins

BACKGROUND

A deficiency in phaseolin and phytohemagglutinin is associated with a near doubling of sulfur amino acid content in genetically related lines of common bean (Phaseolus vulgaris), particularly cysteine, elevated by 70%, and methionine, elevated by 10%. This mostly takes place at the expense of an abundant non-protein amino acid, S-methyl-cysteine. The deficiency in phaseolin and phytohemagglutinin is mainly compensated by increased levels of the 11S globulin legumin and residual lectins. Legumin, albumin-2, defensin and albumin-1 were previously identified as contributing to the increased sulfur amino acid content in the mutant line, on the basis of similarity to proteins from other legumes.

RESULTS

Profiling of free amino acid in developing seeds of the BAT93 reference genotype revealed a biphasic accumulation of gamma-glutamyl-S-methyl-cysteine, the main soluble form of S-methyl-cysteine, with a lag phase occurring during storage protein accumulation. A collection of 30,147 expressed sequence tags (ESTs) was generated from four developmental stages, corresponding to distinct phases of gamma-glutamyl-S-methyl-cysteine accumulation, and covering the transitions to reserve accumulation and dessication. Analysis of gene ontology categories indicated the occurrence of multiple sulfur metabolic pathways, including all enzymatic activities responsible for sulfate assimilation, de novo cysteine and methionine biosynthesis. Integration of genomic and proteomic data enabled the identification and isolation of cDNAs coding for legumin, albumin-2, defensin D1 and albumin-1A and -B induced in the absence of phaseolin and phytohemagglutinin. Their deduced amino acid sequences have a higher content of cysteine than methionine, providing an explanation for the preferential increase of cysteine in the mutant line.

CONCLUSION

The EST collection provides a foundation to further investigate sulfur metabolism and the differential accumulation of sulfur amino acids in seed of common bean. Identification of sulfur-rich proteins whose levels are elevated in seed lacking phaseolin and phytohemagglutinin and sulfur metabolic genes may assist the improvement of protein quality.

The characterization of a new set of EST-derived simple sequence repeat (SSR) markers as a resource for the genetic analysis of Phaseolus vulgaris

Background

Over recent years, a growing effort has been made to develop microsatellite markers for the genomic analysis of the common bean (Phaseolus vulgaris) to broaden the knowledge of the molecular genetic basis of this species. The availability of large sets of expressed sequence tags (ESTs) in public databases has given rise to an expedient approach for the identification of SSRs (Simple Sequence Repeats), specifically EST-derived SSRs. In the present work, a battery of new microsatellite markers was obtained from a search of the Phaseolus vulgaris EST database. The diversity, degree of transferability and polymorphism of these markers were tested.

Results

From 9,583 valid ESTs, 4,764 had microsatellite motifs, from which 377 were used to design primers, and 302 (80.11%) showed good amplification quality. To analyze transferability, a group of 167 SSRs were tested, and the results showed that they were 82% transferable across at least one species. The highest amplification rates were observed between the species from the Phaseolus (63.7%), Vigna (25.9%), Glycine (19.8%), Medicago (10.2%), Dipterix (6%) and Arachis (1.8%) genera. The average PIC (Polymorphism Information Content) varied from 0.53 for genomic SSRs to 0.47 for EST-SSRs, and the average number of alleles per locus was 4 and 3, respectively. Among the 315 newly tested SSRs in the BJ (BAT93 X Jalo EEP558) population, 24% (76) were polymorphic. The integration of these segregant loci into a framework map composed of 123 previously obtained SSR markers yielded a total of 199 segregant loci, of which 182 (91.5%) were mapped to 14 linkage groups, resulting in a map length of 1,157 cM.

Conclusions

A total of 302 newly developed EST-SSR markers, showing good amplification quality, are available for the genetic analysis of Phaseolus vulgaris. These markers showed satisfactory rates of transferability, especially between species that have great economic and genomic values. Their diversity was comparable to genomic SSRs, and they were incorporated in the common bean reference genetic map, which constitutes an important contribution to and advance in Phaseolus vulgaris genomic research.

Gene-based SSR markers for common bean (Phaseolus vulgaris L.) derived from root and leaf tissue ESTs: an integration of the BMc series

BACKGROUND

Sequencing of cDNA libraries for the development of expressed sequence tags (ESTs) as well as for the discovery of simple sequence repeats (SSRs) has been a common method of developing microsatellites or SSR-based markers. In this research, our objective was to further sequence and develop common bean microsatellites from leaf and root cDNA libraries derived from the Andean gene pool accession G19833 and the Mesoamerican gene pool accession DOR364, mapping parents of a commonly used reference map. The root libraries were made from high and low phosphorus treated plants.

RESULTS

A total of 3,123 EST sequences from leaf and root cDNA libraries were screened and used for direct simple sequence repeat discovery. From these EST sequences we found 184 microsatellites; the majority containing tri-nucleotide motifs, many of which were GC rich (ACC, AGC and AGG in particular). Di-nucleotide motif microsatellites were about half as common as the tri-nucleotide motif microsatellites but most of these were AGn microsatellites with a moderate number of ATn microsatellites in root ESTs followed by few ACn and no GCn microsatellites. Out of the 184 new SSR loci, 120 new microsatellite markers were developed in the BMc (Bean Microsatellites from cDNAs) series and these were evaluated for their capacity to distinguish bean diversity in a germplasm panel of 18 genotypes. We developed a database with images of the microsatellites and their polymorphism information content (PIC), which averaged 0.310 for polymorphic markers.

CONCLUSIONS

The present study produced information about microsatellite frequency in root and leaf tissues of two important genotypes for common bean genomics: namely G19833, the Andean genotype selected for whole genome shotgun sequencing from race Peru, and DOR364 a race Mesoamerica subgroup 2 genotype that is a small-red seeded, released variety in Central America. Both race Peru and Mesoamerica subgroup 2 (small red beans) have been understudied in comparison to race Nueva Granada and Mesoamerica subgroup 1 (black beans) both with regards to gene expression and as sources of markers. However, we found few differences between SSR type and frequency between the G19833 leaf and DOR364 root tissue-derived ESTs. Overall, our work adds to the analysis of microsatellite frequency evaluation for common bean and provides a new set of 120 BMc markers which combined with the 248 previously developed BMc markers brings the total in this series to 368 markers. Once we include BMd markers, which are derived from GenBank sequences, the current total of gene-based markers from our laboratory surpasses 500 markers. These markers are basic for studies of the transcriptome of common bean and can form anchor points for genetic mapping studies in the future.

An EST-based analysis identifies new genes and reveals distinctive gene expression features of Coffea arabica and Coffea canephora

BACKGROUND

Coffee is one of the world's most important crops; it is consumed worldwide and plays a significant role in the economy of producing countries. Coffea arabica and C. canephora are responsible for 70 and 30% of commercial production, respectively. C. arabica is an allotetraploid from a recent hybridization of the diploid species, C. canephora and C. eugenioides. C. arabica has lower genetic diversity and results in a higher quality beverage than C. canephora. Research initiatives have been launched to produce genomic and transcriptomic data about Coffea spp. as a strategy to improve breeding efficiency.

RESULTS

Assembling the expressed sequence tags (ESTs) of C. arabica and C. canephora produced by the Brazilian Coffee Genome Project and the Nestlé-Cornell Consortium revealed 32,007 clusters of C. arabica and 16,665 clusters of C. canephora. We detected different GC3 profiles between these species that are related to their genome structure and mating system. BLAST analysis revealed similarities between coffee and grape (Vitis vinifera) genes. Using KA/KS analysis, we identified coffee genes under purifying and positive selection. Protein domain and gene ontology analyses suggested differences between Coffea spp. data, mainly in relation to complex sugar synthases and nucleotide binding proteins. OrthoMCL was used to identify specific and prevalent coffee protein families when compared to five other plant species. Among the interesting families annotated are new cystatins, glycine-rich proteins and RALF-like peptides. Hierarchical clustering was used to independently group C. arabica and C. canephora expression clusters according to expression data extracted from EST libraries, resulting in the identification of differentially expressed genes. Based on these results, we emphasize gene annotation and discuss plant defenses, abiotic stress and cup quality-related functional categories.

CONCLUSION

We present the first comprehensive genome-wide transcript profile study of C. arabica and C. canephora, which can be freely assessed by the scientific community at http://www.lge.ibi.unicamp.br/coffea. Our data reveal the presence of species-specific/prevalent genes in coffee that may help to explain particular characteristics of these two crops. The identification of differentially expressed transcripts offers a starting point for the correlation between gene expression profiles and Coffea spp. developmental traits, providing valuable insights for coffee breeding and biotechnology, especially concerning sugar metabolism and stress tolerance.

Syntenic relationships among legumes revealed using a gene-based genetic linkage map of common bean (Phaseolus vulgaris L.)

Molecular linkage maps are an important tool for gene discovery and cloning, crop improvement, further genetic studies, studies on diversity and evolutionary history, and cross-species comparisons. Linkage maps differ in both the type of marker and type of population used. In this study, gene-based markers were used for mapping in a recombinant inbred (RI) population of Phaseolus vulgaris L. P. vulgaris, common dry bean, is an important food source, economic product, and model organism for the legumes. Gene-based markers were developed that corresponded to genes controlling mutant phenotypes in Arabidopsis thaliana, genes undergoing selection during domestication in maize, and genes that function in a biochemical pathway in A. thaliana. Sequence information, including introns and 3' UTR, was generated for over 550 genes in the two genotypes of P. vulgaris. Over 1,800 single nucleotide polymorphisms and indels were found, 300 of which were screened in the RI population. The resulting LOD 2.0 map is 1,545 cM in length and consists of 275 gene-based and previously mapped core markers. An additional 153 markers that mapped at LOD <1.0 were placed in genetic bins. By screening the parents of other mapping populations, it was determined that the markers were useful for other common Mesoamerican × Andean mapping populations. The location of the mapped genes relative to their homologs in Arabidopsis thaliana (At), Medicago truncatula (Mt), and Lotus japonicus (Lj) were determine by using a tblastx analysis with the current psedouchromosome builds for each of the species. While only short blocks of synteny were observed with At, large-scale macrosyntenic blocks were observed with Mt and Lj. By using Mt and Lj as bridging species, the syntenic relationship between the common bean and peanut was inferred.

High-throughput SNP discovery and assay development in common bean

BACKGROUND

Next generation sequencing has significantly increased the speed at which single nucleotide polymorphisms (SNPs) can be discovered and subsequently used as molecular markers for research. Unfortunately, for species such as common bean (Phaseolus vulgaris L.) which do not have a whole genome sequence available, the use of next generation sequencing for SNP discovery is much more difficult and costly. To this end we developed a method which couples sequences obtained from the Roche 454-FLX system (454) with the Illumina Genome Analyzer (GA) for high-throughput SNP discovery.

RESULTS

Using a multi-tier reduced representation library we discovered a total of 3,487 SNPs of which 2,795 contained sufficient flanking genomic sequence for SNP assay development. Using Sanger sequencing to determine the validation rate of these SNPs, we found that 86% are likely to be true SNPs. Furthermore, we designed a GoldenGate assay which contained 1,050 of the 3,487 predicted SNPs. A total of 827 of the 1,050 SNPs produced a working GoldenGate assay (79%).

CONCLUSIONS

Through combining two next generation sequencing techniques we have developed a method that allows high-throughput SNP discovery in any diploid organism without the need of a whole genome sequence or the creation of normalized cDNA libraries. The need to only perform one 454 run and one GA sequencer run allows high-throughput SNP discovery with sufficient sequence for assay development to be performed in organisms, such as common bean, which have limited genomic resources.

Three highly similar formate dehydrogenase genes located in the vicinity of the B4 resistance gene cluster are differentially expressed under biotic and abiotic stresses in Phaseolus vulgaris

In higher plants, formate dehydrogenase (FDH, EC1.2.1.2.) catalyzes the NAD-linked oxidation of formate to CO(2), and FDH transcript accumulation has been reported after various abiotic stresses. By sequencing a Phaseolus vulgaris BAC clone encompassing a CC-NBS-LRR gene rich region of the B4 resistance gene cluster, we identified three FDH-encoding genes. FDH is present as a single copy gene in the Arabidopsis thaliana genome, and public database searches confirm that FDH is a low copy gene in plant genomes, since only 33 FDH homologs were identified from 27 plant species. Three independent prediction programs (Predotar, TargetP and Mitoprot) used on this large subset of 33 plant FDHs, revealed that mitochondrial localization of FDH might be the rule in higher plants. A phylogenetic analysis suggests a scenario of local FDH gene duplication in an ancestor of the Phaseoleae followed by another more recent duplication event after bean/soybean divergence. The expression levels of two common bean FDH genes under different treatments were investigated by quantitative RT-PCR analysis. FDH genes are differentially up-regulated after biotic and abiotic stresses (infection with the fungus Colletotrichum lindemuthianum, and dark treatment, respectively). The present study provides the first report of FDH transcript accumulation after biotic stress, suggesting the involvement of FDH in the pathogen resistance process.