Transcriptome Characterization of Developing Bean (Phaseolus vulgaris L.) Pods from Two Genotypes with Contrasting Seed Zinc Concentrations.

Dry bean (Phaseolus vulgaris L.) seeds are a rich source of dietary zinc, especially for people consuming plant-based diets. Within P. vulgaris there is at least two-fold variation in seed Zn concentration. Genetic studies have revealed seed Zn differences to be controlled by a single gene in two closely related navy bean genotypes, Albion and Voyager. In this study, these two genotypes were grown under controlled fertilization conditions and the Zn concentration of various plant parts was determined. The two genotypes had similar levels of Zn in their leaves and pods but Voyager had 52% more Zn in its seeds than Albion. RNA was sequenced from developing pods of both genotypes. Transcriptome analysis of these genotypes identified 27,198 genes in the developing bean pods, representing 86% of the genes in the P. vulgaris genome (v 1.0 DOE-JGI and USDA-NIFA). Expression was detected in 18,438 genes. A relatively small number of genes (381) were differentially expressed between Albion and Voyager. Differentially expressed genes included three genes potentially involved in Zn transport, including zinc-regulated transporter, iron regulated transporter like (ZIP), zinc-induced facilitator (ZIF) and heavy metal associated (HMA) family genes. In addition 12,118 SNPs were identified between the two genotypes. Of the gene families related to Zn and/or Fe transport, eleven genes were found to contain SNPs between Albion and Voyager.

The Phaseolus vulgaris ZIP gene family: identification, characterization, mapping, and gene expression.

Zinc is an essential mineral for humans and plants and is involved in many physiological and biochemical processes. In humans, Zn deficiency has been associated with retarded growth and reduction of immune response. In plants, Zn is an essential component of more than 300 enzymes including RNA polymerase, alkaline phosphatase, alcohol dehydrogenase, Cu/Zn superoxidase dismutase, and carbonic anhydrase. The accumulation of Zn in plants involves many genes and characterization of the role of these genes will be useful in biofortification. Here we report the identification and phlyogenetic and sequence characterization of the 23 members of the ZIP (ZRT, IRT like protein) family of metal transporters and three transcription factors of the bZIP family in Phaseolus vulgaris L. Expression patterns of seven of these genes were characterized in two bean genotypes (G19833 and DOR364) under two Zn treatments. Tissue analyzed included roots and leaves at vegetative and flowering stages, and pods at 20 days after flowering. Four of the genes, PvZIP12, PvZIP13, PvZIP16, and Pv bZIP1, showed differential expression based on tissue, Zn treatment, and/or genotype. PvZIP12 and PvZIP13 were both more highly expressed in G19833 than DOR364. PvZIP12 was most highly expressed in vegetative leaves under the Zn (-) treatment. PvZIP16 was highly expressed in leaf tissue, especially leaf tissue at flowering stage grown in the Zn (-) treatment. Pv bZIP1 was most highly expressed in leaf and pod tissue. The 23 PvZIP genes and three bZIP genes were mapped on the DOR364 × G19833 linkage map. PvZIP12, PvZIP13, and PvZIP18, Pv bZIP2, and Pv bZIP3 were located near QTLs for Zn accumulation in the seed. Based on the expression and mapping results, PvZIP12 is a good candidate gene for increasing seed Zn concentration and increase understanding of the role of ZIP genes in metal uptake, distribution, and accumulation of zinc in P. vulgaris.

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.

Saturation of an intra-gene pool linkage map: towards a unified consensus linkage map for fine mapping and synteny analysis in common bean.

Map-based cloning and fine mapping to find genes of interest and marker assisted selection (MAS) requires good genetic maps with reproducible markers. In this study, we saturated the linkage map of the intra-gene pool population of common bean DOR364 × BAT477 (DB) by evaluating 2,706 molecular markers including SSR, SNP, and gene-based markers. On average the polymorphism rate was 7.7% due to the narrow genetic base between the parents. The DB linkage map consisted of 291 markers with a total map length of 1,788 cM. A consensus map was built using the core mapping populations derived from inter-gene pool crosses: DOR364 × G19833 (DG) and BAT93 × JALO EEP558 (BJ). The consensus map consisted of a total of 1,010 markers mapped, with a total map length of 2,041 cM across 11 linkage groups. On average, each linkage group on the consensus map contained 91 markers of which 83% were single copy markers. Finally, a synteny analysis was carried out using our highly saturated consensus maps compared with the soybean pseudo-chromosome assembly. A total of 772 marker sequences were compared with the soybean genome. A total of 44 syntenic blocks were identified. The linkage group Pv6 presented the most diverse pattern of synteny with seven syntenic blocks, and Pv9 showed the most consistent relations with soybean with just two syntenic blocks. Additionally, a co-linear analysis using common bean transcript map information against soybean coding sequences (CDS) revealed the relationship with 787 soybean genes. The common bean consensus map has allowed us to map a larger number of markers, to obtain a more complete coverage of the common bean genome. Our results, combined with synteny relationships provide tools to increase marker density in selected genomic regions to identify closely linked polymorphic markers for indirect selection, fine mapping or for positional cloning.

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.

Variation and inheritance of iron reductase activity in the roots of common bean (Phaseolus vulgaris L.) and association with seed iron accumulation QTL.

BACKGROUND: Iron deficiency anemia is a global problem which often affects women and children of developing countries. Strategy I plants, such as common bean (Phaseolus vulgaris L.) take up iron through a process that involves an iron reduction mechanism in their roots; this reduction is required to convert ferric iron to ferrous iron. Root absorbed iron is critical for the iron nutrition of the plant, and for the delivery of iron to the shoot and ultimately the seeds. The objectives of this study were to determine the variability and inheritance for iron reductase activity in a range of genotypes and in a low × high seed iron cross (DOR364 x G19833), to identify quantitative trait loci (QTL) for this trait, and to assess possible associations with seed iron levels. RESULTS: The experiments were carried out with hydroponically grown plants provided different amounts of iron varying between 0 and 20 µM Fe(III)-EDDHA. The parents, DOR364 and G19833, plus 13 other cultivated or wild beans, were found to differ in iron reductase activity. Based on these initial experiments, two growth conditions (iron limited and iron sufficient) were selected as treatments for evaluating the DOR364 × G19833 recombinant inbred lines. A single major QTL was found for iron reductase activity under iron-limited conditions (1 µM Fe) on linkage group b02 and another major QTL was found under iron sufficient conditions (15 µM Fe) on linkage group b11. Associations between the b11 QTL were found with several QTL for seed iron. CONCLUSIONS: Genes conditioning iron reductase activity in iron sufficient bean plants appear to be associated with genes contributing to seed iron accumulation. Markers for bean iron reductase (FRO) homologues were found with in silico mapping based on common bean synteny with soybean and Medicago truncatula on b06 and b07; however, neither locus aligned with the QTL for iron reductase activity. In summary, the QTL for iron reductase activity under iron limited conditions may be useful in environments where beans are grown in alkaline soils, while the QTL for iron reductase under sufficiency conditions may be useful for selecting for enhanced seed nutritional quality.

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

BACKGROUND: Expressed sequence tags (ESTs) are an important source of gene-based markers such as those based on insertion-deletions (Indels) or single-nucleotide polymorphisms (SNPs). Several gel based methods have been reported for the detection of sequence variants, however they have not been widely exploited in common bean, an important legume crop of the developing world. The objectives of this project were to develop and map EST based markers using analysis of single strand conformation polymorphisms (SSCPs), to create a transcript map for common bean and to compare synteny of the common bean map with sequenced chromosomes of other legumes. RESULTS: A set of 418 EST based amplicons were evaluated for parental polymorphisms using the SSCP technique and 26% of these presented a clear conformational or size polymorphism between Andean and Mesoamerican genotypes. The amplicon based markers were then used for genetic mapping with segregation analysis performed in the DOR364 x G19833 recombinant inbred line (RIL) population. A total of 118 new marker loci were placed into an integrated molecular map for common bean consisting of 288 markers. Of these, 218 were used for synteny analysis and 186 presented homology with segments of the soybean genome with an e-value lower than 7 x 10-12. The synteny analysis with soybean showed a mosaic pattern of syntenic blocks with most segments of any one common bean linkage group associated with two soybean chromosomes. The analysis with Medicago truncatula and Lotus japonicus presented fewer syntenic regions consistent with the more distant phylogenetic relationship between the galegoid and phaseoloid legumes. CONCLUSION: The SSCP technique is a useful and inexpensive alternative to other SNP or Indel detection techniques for saturating the common bean genetic map with functional markers that may be useful in marker assisted selection. In addition, the genetic markers based on ESTs allowed the construction of a transcript map and given their high conservation between species allowed synteny comparisons to be made to sequenced genomes. This synteny analysis may support positional cloning of target genes in common bean through the use of genomic information from these other legumes.