Common Bean Genome Research Articles

Integrating genetic maps in bambara groundnut [Vigna subterranea (L) Verdc.] and their syntenic relationships among closely related legumes

BackgroundBambara groundnut [Vigna subterranea (L) Verdc.] is an indigenous legume crop grown mainly in subsistence and small-scale agriculture in sub-Saharan Africa for its nutritious seeds and its tolerance to drought and poor soils. Given that the lack of ex ante sequence is often a bottleneck in marker-assisted crop breeding for minor and underutilised crops, we demonstrate the use of limited genetic information and resources developed within species, but linked to the well characterised common bean (Phaseolus vulgaris) genome sequence and the partially annotated closely related species; adzuki bean (Vigna angularis) and mung bean (Vigna radiata). From these comparisons we identify conserved synteny blocks corresponding to the Linkage Groups (LGs) in bambara groundnut genetic maps and evaluate the potential to identify genes in conserved syntenic locations in a sequenced genome that underlie a QTL position in the underutilised crop genome.ResultsTwo individual intraspecific linkage maps consisting of DArTseq markers were constructed in two bambara groundnut (2n = 2x = 22) segregating populations: 1) The genetic map of Population IA was derived from F2 lines (n = 263; IITA686 x Ankpa4) and covered 1,395.2 cM across 11 linkage groups; 2) The genetic map of Population TD was derived from F3 lines (n = 71; Tiga Nicuru x DipC) and covered 1,376.7 cM across 11 linkage groups. A total of 96 DArTseq markers from an initial pool of 142 pre-selected common markers were used. These were not only polymorphic in both populations but also each marker could be located using the unique sequence tag (at selected stringency) onto the common bean, adzuki bean and mung bean genomes, thus allowing the sequenced genomes to be used as an initial ‘pseudo’ physical map for bambara groundnut. A good correspondence was observed at the macro synteny level, particularly to the common bean genome. A test using the QTL location of an agronomic trait in one of the bambara groundnut maps allowed the corresponding flanking positions to be identified in common bean, mung bean and adzuki bean, demonstrating the possibility of identifying potential candidate genes underlying traits of interest through the conserved syntenic physical location of QTL in the well annotated genomes of closely related species.ConclusionsThe approach of adding pre-selected common markers in both populations before genetic map construction has provided a translational framework for potential identification of candidate genes underlying a QTL of trait of interest in bambara groundnut by linking the positions of known genetic effects within the underutilised species to the physical maps of other well-annotated legume species, without the need for an existing whole genome sequence of the study species. Identifying the conserved synteny between underutilised species without complete genome sequences and the genomes of major crops and model species with genetic and trait data is an important step in the translation of resources and information from major crop and model species into the minor crop species. Such minor crops will be required to play an important role in future agriculture under the effects of climate change.

Genome-wide in silico identification, characterization and transcriptional analysis of the family of growth-regulating factors in common bean (Phaseolus vulgaris L.) subjected to polyethylene glycol-induced drought stress

According to most recent findings, growth regulating factors (GRFs) are plant-specific transcription factors (TFs) that play important roles in many processes, including abiotic and biotic stress response mechanisms. Completion of the common bean (Phaseolus vulgaris) genome project has provided researchers with the opportunity to identify all GRF genes in this species. With this aim, a genome-wide in silico study was performed and 10 GRF proteins (called PhvGRFs) were identified in the common bean genome. Conserved and mandatory motifs (QLQ and WRC) were confirmed in all identified PhvGRFs and two segmental duplication events were determined. Most of the PhvGRFs were found to be more similar to Arabidopsis thaliana GRFs than to Zea mays GRFs in a phylogenetic tree. According to the expression analysis of 10 PhvGRFs, inversely related expression patterns were observed in the roots of Yakutiye and Zulbiye cultivars based on their capacity to adopt to drought stress. After drought treatment of the Zulbiye cultivar, a drought-sensitive common bean cultivar, PhvGRF1, PhvGRF2, PhvGRF3, PhvGRF5, PhvGRF6, PhvGRF9 and PhvGRF10 genes were upregulated 2- to 4-fold in root tissues, as compared to the untreated control. The trend of PhvGRF1, PhvGRF2, PhvGRF3, PhvGRF5, PhvGRF6, PhvGRF7, PhvGRF9 and PhvGRF10 genes showed a consistent decline of 2- to 6-fold in root tissues of the drought-tolerant Yakutiye cultivar subjected to 24 h of drought stress. We demonstrated that the expression patterns of the identified PhvGRFs correlated with the drought-stress response in a cultivar-specific manner in the common bean. We suggest that members of the GRF family can also be used for genetic engineering applications in the common bean.

In silico analysis of galactinol synthase genes in common bean (Phaseolus vulgaris L.)

Galactinol synthases (GolS, EC 2.4.1.123) belong to the glycosyltransferase 8 family. These enzymes catalyzes galactosylation of myo-inositol to form galactinol, precursor of raffinose which is solube sucrose playing an important role in stress tolerance in plants. We identified and analyzed a total of three GolS encoded genes in the whole genome of common bean. These genes possessed two or three introns. The predicted proteins contained from 322 to 340 amino acids. The common bean GolS proteins contained all conseved motifs that were characterized for other previously known GolSs. These proteins were acidic and weakly hydrophilic. The phylogenic tree constructed from GolS proteins of Arabidopsis, poplar and common bean showed that only one gene duplication event was detected in the commone bean genome. Common bean GolS genes expressed in all vegetative and reproductive tissues. However, the transcript abundances of each gene varied in different tissues. When comparing the expression levels of these three genes as defined RNAseq, we found that the expression level of PvGolS2 was the lowest, PvGolS1 expressed the strogest leaves and stem, while PvGolS3 expressed in roots and nodules. We discovered that Rhizobium incubation affected on the expression of GolS genes for the first time. In the reproductive tissues, PvGolS1 expressed in most studied tissues. Its expression level correlatively increased with the tissue development.

Development of a high-density genetic linkage map and identification of flowering time QTLs in adzuki bean (Vigna angularis).

A high-density linkage map is crucial for the identification of quantitative trait loci (QTLs), positional cloning, and physical map assembly. Here, we report the development of a high-density linkage map based on specific length amplified fragment sequencing (SLAF-seq) for adzuki bean and the identification of flowering time-related QTLs. Through SLAF library construction and Illumina sequencing of a recombinant inbred line (RIL) population, a total of 4425 SLAF markers were developed and assigned to 11 linkage groups (LGs). After binning the SLAF markers that represented the same genotype, the final linkage map of 1628.15 cM contained 2032 markers, with an average marker density of 0.80 cM. Comparative analysis showed high collinearity with two adzuki bean physical maps and a high degree of synteny with the reference genome of common bean (Phaseolus vulgaris). Using this map, one major QTL on LG03 and two minor QTLs on LG05 associated with first flowering time (FLD) were consistently identified in tests over a two-year period. These results provide a foundation that will be useful for future genomic research, such as identifying QTLs for other important traits, positional cloning, and comparative mapping in legumes.

RNAseq Analysis of Endornavirus-Infected vs. Endornavirus-Free Common Bean (Phaseolus vulgaris) Cultivar Black Turtle Soup.

Common bean (Phaseolus vulgaris L.) is the most important grain legume for direct human consumption worldwide and represents a rich source of protein, vitamins, minerals, and fiber (Broughton et al., 2003). The recent sequencing of the common bean genome, together with the availability of genomic and transcriptomic data have provided useful information to common bean breeders that will help in the development of genotypes with desirable characteristics (Schmutz et al., 2014; Vlasova et al., 2016). Endornaviruses are persistent viruses with a non-encapsidated RNA genome that ranges from 9.8 to 17.6 kb, infect plants, fungi, and oomycetes, are transmitted only via gametes, and do not cause apparent symptoms (Stielow et al., 2011; Fukuhara and Gibbs, 2012). Although endornaviruses have been reported in several economically important plant species, little is known about the effect they have on their hosts. One of the major obstacles to study their effect to the host is the lack of a transmission method. In plants, endornaviruses do not move from cell to cell and spread only during cell division. Recently, Khankhum et al. (2015) reported that most common bean genotypes of Mesoamerican origin are double-infected with Phaseolus vulgaris endornavirus 1 (PvEV1) and Phaseolus endornavirus 2 (PvEV2); in contrast, genotypes of Andean origin are often endornavirus-free. Black Turtle Soup (BTS), a cultivar of Mesoamerican origin has been reported to be double-infected by these two endornaviruses (Okada et al., 2013). A BTS endornavirus-free selection (BTS−), obtained from an endornavirus-infected BTS (BTS+) seed lot has been reported by Okada et al. (2013). To establish the bases for future research on the role that endornaviruses play in the common bean plant, and the effect these viruses have on the host gene expression, we conducted RNAseq on two BTS lines: one endornavirus-infected and the other endornavirus-free.

Comprehensive analysis and discovery of drought-related NAC transcription factors in common bean.

BackgroundCommon bean (Phaseolus vulgaris L.) is an important warm-season food legume. Drought is the most important environmental stress factor affecting large areas of common bean via plant death or reduced global production. The NAM, ATAF1/2 and CUC2 (NAC) domain protein family are classic transcription factors (TFs) involved in a variety of abiotic stresses, particularly drought stress. However, the NAC TFs in common bean have not been characterized.ResultsIn the present study, 86 putative NAC TF proteins were identified from the common bean genome database and located on 11 common bean chromosomes. The proteins were phylogenetically clustered into 8 distinct subfamilies. The gene structure and motif composition of common bean NACs were similar in each subfamily. These results suggest that NACs in the same subfamily may possess conserved functions. The expression patterns of common bean NAC genes were also characterized. The majority of NACs exhibited specific temporal and spatial expression patterns. We identified 22 drought-related NAC TFs based on transcriptome data for drought-tolerant and drought-sensitive genotypes. Quantitative real-time PCR (qRT-PCR) was performed to confirm the expression patterns of the 20 drought-related NAC genes.ConclusionsBased on the common bean genome sequence, we analyzed the structural characteristics, genome distribution, and expression profiles of NAC gene family members and analyzed drought-responsive NAC genes. Our results provide useful information for the functional characterization of common bean NAC genes and rich resources and opportunities for understanding common bean drought stress tolerance mechanisms.Electronic supplementary materialThe online version of this article (doi:10.1186/s12870-016-0882-5) contains supplementary material, which is available to authorized users.

Dynamics of a Novel Highly Repetitive CACTA Family in Common Bean (Phaseolus vulgaris).

Transposons are ubiquitous genomic components that play pivotal roles in plant gene and genome evolution. We analyzed two genome sequences of common bean (Phaseolus vulgaris) and identified a new CACTA transposon family named pvCACTA1. The family is extremely abundant, as more than 12,000 pvCACTA1 elements were found. To our knowledge, this is the most abundant CACTA family reported thus far. The computational and fluorescence in situ hybridization (FISH) analyses indicated that the pvCACTA1 elements were concentrated in terminal regions of chromosomes and frequently generated AT-rich 3 bp target site duplications (TSD, WWW, W is A or T). Comparative analysis of the common bean genomes from two domesticated genetic pools revealed that new insertions or excisions of pvCACTA1 elements occurred after the divergence of the two common beans, and some of the polymorphic elements likely resulted in variation in gene sequences. pvCACTA1 elements were detected in related species but not outside the Phaseolus genus. We calculated the molecular evolutionary rate of pvCACTA1 transposons using orthologous elements that indicated that most transposition events likely occurred before the divergence of the two gene pools. These results reveal unique features and evolution of this new transposon family in the common bean genome.

Sequence-Based Introgression Mapping Identifies Candidate White Mold Tolerance Genes in Common Bean.

White mold, caused by the necrotrophic fungus (Lib.) de Bary, is a major disease of common bean ( L.). WM7.1 and WM8.3 are two quantitative trait loci (QTL) with major effects on tolerance to the pathogen. Advanced backcross populations segregating individually for either of the two QTL, and a recombinant inbred (RI) population segregating for both QTL were used to fine map and confirm the genetic location of the QTL. The QTL intervals were physically mapped using the reference common bean genome sequence, and the physical intervals for each QTL were further confirmed by sequence-based introgression mapping. Using whole-genome sequence data from susceptible and tolerant DNA pools, introgressed regions were identified as those with significantly higher numbers of single-nucleotide polymorphisms (SNPs) relative to the whole genome. By combining the QTL and SNP data, WM7.1 was located to a 660-kb region that contained 41 gene models on the proximal end of chromosome Pv07, while the WM8.3 introgression was narrowed to a 1.36-Mb region containing 70 gene models. The most polymorphic candidate gene in the WM7.1 region encodes a BEACH-domain protein associated with apoptosis. Within the WM8.3 interval, a receptor-like protein with the potential to recognize pathogen effectors was the most polymorphic gene. The use of gene and sequence-based mapping identified two candidate genes whose putative functions are consistent with the current model of pathogenicity.

English

Expressed sequences are important sources in the development of functional heterologous microsatellite markers in phylogenetic related groups, that is, soybean and commom bean. The objective of this work was to identify and characterize expressed sequence tag (EST)-simple sequence repeats (SSRs) in silico candidates of resistance to Meloidogyne spp. in common bean. Seven DNA sequences from soybean associated with genetic resistance have been identified and obtained in the NCBI database. Its homology in common bean genome was verified using the BLAST tool. The cellular processes involved were also checked using the Blast2GO program. The identification of microsatellite markers and design of the primer pairs was performed using the SSRLocator and Primer3 programs, respectively. The transferability rate of common beans the target sequences identified was 86%, demonstrating the power of success of this method. All the cellular processes involved in the original DNA sequences were verified from EST on beans, with E-value between 0 and 2.9×10-168. Fifteen EST-SSRs candidates for common bean resistance were identified, which have proved to be suitable for their amplification by PCR. The transferability analysis of ESTs related to resistance to Meloidogyne spp., especially among soybeans and common beans is efficient. Based on this study, 15 EST-SSRs candidates are available for validation and later use. Key words: simple sequence repeats (SSRs) markers, assisted selection, ontology, pre-breeding.

Common Bean: A Legume Model on the Rise for Unraveling Responses and Adaptations to Iron, Zinc, and Phosphate Deficiencies.

Common bean (Phaseolus vulgaris) was domesticated ∼8000 years ago in the Americas and today is a staple food worldwide. Besides caloric intake, common bean is also an important source of protein and micronutrients and it is widely appreciated in developing countries for their affordability (compared to animal protein) and its long storage life. As a legume, common bean also has the economic and environmental benefit of associating with nitrogen-fixing bacteria, thus reducing the use of synthetic fertilizers, which is key for sustainable agriculture. Despite significant advances in the plant nutrition field, the mechanisms underlying the adaptation of common bean to low nutrient input remains largely unknown. The recent release of the common bean genome offers, for the first time, the possibility of applying techniques and approaches that have been exclusive to model plants to study the adaptive responses of common bean to challenging environments. In this review, we discuss the hallmarks of common bean domestication and subsequent distribution around the globe. We also discuss recent advances in phosphate, iron, and zinc homeostasis, as these nutrients often limit plant growth, development, and yield. In addition, iron and zinc are major targets of crop biofortification to improve human nutrition. Developing common bean varieties able to thrive under nutrient limiting conditions will have a major impact on human nutrition, particularly in countries where dry beans are the main source of carbohydrates, protein and minerals.

Gene-based SNP discovery in tepary bean (Phaseolus acutifolius) and common bean (P. vulgaris) for diversity analysis and comparative mapping.

BackgroundCommon bean (Phaseolus vulgaris) is an important grain legume and there has been a recent resurgence in interest in its relative, tepary bean (P. acutifolius), owing to this species’ ability to better withstand abiotic stresses. Genomic resources are scarce for this minor crop species and a better knowledge of the genome-level relationship between these two species would facilitate improvement in both. High-throughput genotyping has facilitated large-scale single nucleotide polymorphism (SNP) identification leading to the development of molecular markers with associated sequence information that can be used to place them in the context of a full genome assembly.ResultsTranscript-based SNPs were identified from six common bean and two tepary bean accessions and a subset were used to generate a 768-SNP Illumina GoldenGate assay for each species. The tepary bean assay was used to assess diversity in wild and cultivated tepary bean and to generate the first gene-based map of the tepary bean genome. Genotypic analyses of the diversity panel showed a clear separation between domesticated and cultivated tepary beans, two distinct groups within the domesticated types, and P. parvifolius was confirmed to be distinct. The genetic map of tepary bean was compared to the common bean genome assembly to demonstrate high levels of collinearity between the two species with differences limited to a few intra-chromosomal rearrangements.ConclusionsThe development of the first set of genomic resources specifically for tepary bean has allowed for greater insight into the structure of this species and its relationship to its agriculturally more prominent relative, common bean. These resources will be helpful in the development of efficient breeding strategies for both species and will facilitate the introgression of agriculturally important traits from one crop into the other.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-2499-3) contains supplementary material, which is available to authorized users.

Genome-Wide Association Studies of Anthracnose and Angular Leaf Spot Resistance in Common Bean (Phaseolus vulgaris L.)

The common bean (Phaseolus vulgaris L.) is the world’s most important legume for human consumption. Anthracnose (ANT; Colletotrichum lindemuthianum) and angular leaf spot (ALS; Pseudocercospora griseola) are complex diseases that cause major yield losses in common bean. Depending on the cultivar and environmental conditions, anthracnose and angular leaf spot infections can reduce crop yield drastically. This study aimed to estimate linkage disequilibrium levels and identify quantitative resistance loci (QRL) controlling resistance to both ANT and ALS diseases of 180 accessions of common bean using genome-wide association analysis. A randomized complete block design with four replicates was performed for the ANT and ALS experiments, with four plants per genotype in each replicate. Association mapping analyses were performed for ANT and ALS using a mixed linear model approach implemented in TASSEL. A total of 17 and 11 significant statistically associations involving SSRs were detected for ANT and ALS resistance loci, respectively. Using SNPs, 21 and 17 significant statistically associations were obtained for ANT and angular ALS, respectively, providing more associations with this marker. The SSR-IAC167 and PvM95 markers, both located on chromosome Pv03, and the SNP scaffold00021_89379, were associated with both diseases. The other markers were distributed across the entire common bean genome, with chromosomes Pv03 and Pv08 showing the greatest number of loci associated with ANT resistance. The chromosome Pv04 was the most saturated one, with six markers associated with ALS resistance. The telomeric region of this chromosome showed four markers located between approximately 2.5 Mb and 4.4 Mb. Our results demonstrate the great potential of genome-wide association studies to identify QRLs related to ANT and ALS in common bean. The results indicate a quantitative and complex inheritance pattern for both diseases in common bean. Our findings will contribute to more effective screening of elite germplasm to find resistance alleles for marker-assisted selection in breeding programs.

Genome-wide investigation and expression analysis of AP2-ERF gene family in salt tolerant common bean.

Apetala2-ethylene-responsive element binding factor (AP2-ERF) superfamily with common AP2-DNA binding domain have developmentally and physiologically important roles in plants. Since common bean genome project has been completed recently, it is possible to identify all of the AP2-ERF genes in the common bean genome. In this study, a comprehensive genome-wide in silico analysis identified 180 AP2-ERF superfamily genes in common bean (Phaseolus vulgaris). Based on the amino acid alignment and phylogenetic analyses, superfamily members were classified into four subfamilies: DREB (54), ERF (95), AP2 (27) and RAV (3), as well as one soloist. The physical and chemical characteristics of amino acids, interaction between AP2-ERF proteins, cis elements of promoter region of AP2-ERF genes and phylogenetic trees were predicted and analyzed. Additionally, expression levels of AP2-ERF genes were evaluated by in silico and qRT-PCR analyses. In silico micro-RNA target transcript analyses identified nearly all PvAP2-ERF genes as targets of by 44 different plant species' miRNAs were identified in this study. The most abundant target genes were PvAP2/ERF-20-25-62-78-113-173. miR156, miR172 and miR838 were the most important miRNAs found in targeting and BLAST analyses. Interactome analysis revealed that the transcription factor PvAP2-ERF78, an ortholog of Arabidopsis At2G28550, was potentially interacted with at least 15 proteins, indicating that it was very important in transcriptional regulation. Here we present the first study to identify and characterize the AP2-ERF transcription factors in common bean using whole-genome analysis, and the findings may serve as a references for future functional research on the transcription factors in common bean.

Development of molecular markers linked to disease resistance genes in common bean based on whole genome sequence

Common bean (Phaseolus vulgaris) is the most important grain legume for direct human consumption in the world, particularly in developing countries where it constitutes the main source of protein. Unfortunately, common bean yield stability is constrained by a number of pests and diseases. As use of resistant genotypes is the most economic and ecologically safe means for controlling plant diseases, efforts have been made to genetically characterize resistance genes (R genes) in common bean. Despite its agronomic importance, genomic resources available in common bean were limited until the recent sequencing of common bean genome (Andean genotype G19833). Besides allowing the annotation of Nucleotide Binding-Leucine Rich Repeat (NB-LRR) encoding gene family, which is the prevalent class of disease R genes in plants, access to the whole genome sequence of common bean can be of great help for intense selection to increase the overall efficiency of crop improvement programs using marker-assisted selection (MAS). This review presents the state of the art of common bean NB-LRR gene clusters, their peculiar location in subtelomeres and correlation with genetically characterized monogenic R genes, as well as how the availability of the whole genome sequence can boost the development of molecular markers for MAS.

Genome-wide identification and in silico characterisation of microRNAs, their targets and processing pathway genes in Phaseolus vulgaris L.

Common bean (Phaseolus vulgaris L., Fabaceae) is a globally important staple crop, which is an important source of calories, protein and essential micronutrients. At the genomic level little is known regarding the small non-coding RNAs within the common bean genome. One of the most important classes of such small non-coding RNAs is microRNAs (miRNAs), which control mRNA and protein expression levels in many eukaryotes. Computational methods have been applied to identify putative miRNAs in the genomes of different organisms. In this study, our objective was to comprehensively identify and characterise miRNAs from the genome and transcriptome of P. vulgaris, including both mature and precursor miRNA forms. We also sought to identify the putative proteins involved in miRNA processing and the likely target genes of common bean miRNAs. We identified 221 mature miRNAs and 136 precursor miRNAs distributed across 52 different miRNA families in the P.vulgaris genome. Amongst these, we distinguished 129 novel mature miRNAs and 123 miRNA precursors belonging to 24 different miRNA families. We also identified 31 proteins predicted to participate in the miRNA-processing pathway in P.vulgaris. Finally, we also identified 483 predicted miRNA targets, including many which corroborate results from other species, suggesting that miRNA regulatory systems are evolutionarily conserved and important for plant development. Our results expand the study of miRNAs and their target genes in common bean, and provide new opportunities to understand their roles in the biology of this important staple crop.

Genome-wide identification of the Phaseolus vulgaris sRNAome using small RNA and degradome sequencing.

BackgroundMiRNAs and phasiRNAs are negative regulators of gene expression. These small RNAs have been extensively studied in plant model species but only 10 mature microRNAs are present in miRBase version 21, the most used miRNA database, and no phasiRNAs have been identified for the model legume Phaseolus vulgaris. Thanks to the recent availability of the first version of the common bean genome, degradome data and small RNA libraries, we are able to present here a catalog of the microRNAs and phasiRNAs for this organism and, particularly, we suggest new protagonists in the symbiotic nodulation events.ResultsWe identified a set of 185 mature miRNAs, including 121 previously unpublished sequences, encoded by 307 precursors and distributed in 98 families. Degradome data allowed us to identify a total of 181 targets for these miRNAs. We reveal two regulatory networks involving conserved miRNAs: those known to play crucial roles in the establishment of nodules, and novel miRNAs present only in common bean, suggesting a specific role for these sequences. In addition, we identified 125 loci that potentially produce phased small RNAs, with 47 of them having all the characteristics of being triggered by a total of 31 miRNAs, including 14 new miRNAs identified in this study.ConclusionsWe provide here a set of new small RNAs that contribute to the broader knowledge of the sRNAome of Phaseolus vulgaris. Thanks to the identification of the miRNA targets from degradome analysis and the construction of regulatory networks between the mature microRNAs, we present here the probable functional regulation associated with the sRNAome and, particularly, in N2-fixing symbiotic nodules.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-015-1639-5) contains supplementary material, which is available to authorized users.

Investigation of genes encoding calcineurin B-like protein family in legumes and their expression analyses in chickpea (Cicer arietinum L.).

Calcium ion (Ca2+) is a ubiquitous second messenger that transmits various internal and external signals including stresses and, therefore, is important for plants’ response process. Calcineurin B-like proteins (CBLs) are one of the plant calcium sensors, which sense and convey the changes in cytosolic Ca2+-concentration for response process. A search in four leguminous plant (soybean, Medicago truncatula, common bean and chickpea) genomes identified 9 to 15 genes in each species that encode CBL proteins. Sequence analyses of CBL peptides and coding sequences (CDS) suggested that there are nine original CBL genes in these legumes and some of them were multiplied during whole genome or local gene duplication. Coding sequences of chickpea CBL genes (CaCBL) were cloned from their cDNAs and sequenced, and their annotations in the genome assemblies were corrected accordingly. Analyses of protein sequences and gene structures of CBL family in plant kingdom indicated its diverse origin but showed a remarkable conservation in overall protein structure with appearance of complex gene structure in the course of evolution. Expression of CaCBL genes in different tissues and in response to different stress and hormone treatment were studied. Most of the CaCBL genes exhibited high expression in flowers. Expression profile of CaCBL genes in response to different abiotic stresses and hormones related to development and stresses (ABA, auxin, cytokinin, SA and JA) at different time intervals suggests their diverse roles in development and plant defence in addition to abiotic stress tolerance. These data not only contribute to a better understanding of the complex regulation of chickpea CBL gene family, but also provide valuable information for further research in chickpea functional genomics.

Genome-wide identification and characterization of aquaporin gene family in common bean (Phaseolus vulgaris L.).

Plant aquaporins are a large and diverse family of water channel proteins that are essential for several physiological processes in living organisms. Numerous studies have linked plant aquaporins with a plethora of processes, such as nutrient acquisition, CO2 transport, plant growth and development, and response to abiotic stresses. However, little is known about this protein family in common bean. Here, we present a genome-wide identification of the aquaporin gene family in common bean (Phaseolus vulgaris L.), a legume crop essential for human nutrition. We identified 41 full-length coding aquaporin sequences in the common bean genome, divided by phylogenetic analysis into five sub-families (PIPs, TIPs, NIPs, SIPs and XIPs). Residues determining substrate specificity of aquaporins (i.e., NPA motifs and ar/R selectivity filter) seem conserved between common bean and other plant species, allowing inference of substrate specificity for these proteins. Thanks to the availability of RNA-sequencing datasets, expression levels in different organs and in leaves of wild and domesticated bean accessions were evaluated. Three aquaporins (PvTIP1;1, PvPIP2;4 and PvPIP1;2) have the overall highest mean expressions, with PvTIP1;1 having the highest expression among all aquaporins. We performed an EST database mining to identify drought-responsive aquaporins in common bean. This analysis showed a significant increase in expression for PvTIP1;1 in drought stress conditions compared to well-watered environments. The pivotal role suggested for PvTIP1;1 in regulating water homeostasis and drought stress response in the common bean should be verified by further field experimentation under drought stress.

Genetic divergence of common bean cultivars.

The aim of this study was to evaluate genetic divergence in the 'Carioca' (beige with brown stripes) common bean cultivar used by different institutions and in 16 other common bean cultivars used in the Rede Cooperativa de Pesquisa de Feijão (Cooperative Network of Common Bean Research), by using simple sequence repeats associated with agronomic traits that are highly distributed in the common bean genome. We evaluated 22 polymorphic loci using bulks containing DNA from 30 plants. There was genetic divergence among the Carioca cultivar provided by the institutions. Nevertheless, there was lower divergence among them than among the other cultivars. The cultivar used by Instituto Agronômico do Paraná was the most divergent in relation to the Carioca samples. The least divergence was observed among the samples used by Universidade Federal de Lavras and by Embrapa Arroz e Feijão. Of all the cultivars, 'CNFP 10104' and 'BRSMG Realce' showed the greatest dissimilarity. The cultivars were separated in two groups of greatest similarity using the Structure software. Genetic variation among cultivars was greater than the variation within or between the groups formed. This fact, together with the high estimate of heterozygosity observed and the genetic divergence of the samples of the Carioca cultivar in relation to the original provided by Instituto Agronômico de Campinas, indicates a mixture of cultivars. The high divergence among cultivars provides potential for the utilization of this genetic variability in plant breeding.

Candidate gene analysis for determinacy in pigeonpea (Cajanus spp.).

Key message We report a likely candidate gene, CcTFL1, for determinacy in pigeonpea through candidate gene sequencing analysis, mapping, QTL analysis together with comparative genomics and expression profiling. Pigeonpea (Cajanus cajan) is the sixth most important legume crop grown on ~5 million hectares globally. Determinacy is an agronomically important trait selected during pigeonpea domestication. In the present study, seven genes related to determinacy/flowering pattern in pigeonpea were isolated through a comparative genomics approach. Single nucleotide polymorphism (SNP) analysis of these candidate genes on 142 pigeonpea lines found a strong association of SNPs with the determinacy trait for three of the genes. Subsequently, QTL analysis highlighted one gene, CcTFL1, as a likely candidate for determinacy in pigeonpea since it explained 45–96 % of phenotypic variation for determinacy, 45 % for flowering time and 77 % for plant height. Comparative genomics analysis of CcTFL1 with the soybean (Glycine max) and common bean (Phaseolus vulgaris) genomes at the micro-syntenic level further enhanced our confidence in CcTFL1 as a likely candidate gene. These findings have been validated by expression analysis that showed down regulation of CcTFL1 in a determinate line in comparison to an indeterminate line. Gene-based markers developed in the present study will allow faster manipulation of the determinacy trait in future breeding programs of pigeonpea and will also help in the development of markers for these traits in other related legume species.Electronic supplementary materialThe online version of this article (doi:10.1007/s00122-014-2406-8) contains supplementary material, which is available to authorized users.