BackgroundCommon bean is a legume of social and nutritional importance as a food crop, cultivated worldwide especially in developing countries, accounting for an important source of income for small farmers. The availability of the complete sequences of the two common bean genomes has dramatically accelerated and has enabled new experimental strategies to be applied for genetic research. DArTseq has been widely used as a method of SNP genotyping allowing comprehensive genome coverage with genetic applications in common bean breeding programs.ResultsUsing this technology, 6286 SNPs (1 SNP/86.5 Kbp) were genotyped in genic (43.3%) and non-genic regions (56.7%). Genetic subdivision associated to the common bean gene pools (K = 2) and related to grain types (K = 3 and K = 5) were reported. A total of 83% and 91% of all SNPs were polymorphic within the Andean and Mesoamerican gene pools, respectively, and 26% were able to differentiate the gene pools. Genetic diversity analysis revealed an average H E of 0.442 for the whole collection, 0.102 for Andean and 0.168 for Mesoamerican gene pools (F ST = 0.747 between gene pools), 0.440 for the group of cultivars and lines, and 0.448 for the group of landrace accessions (F ST = 0.002 between cultivar/line and landrace groups). The SNP effects were predicted with predominance of impact on non-coding regions (77.8%). SNPs under selection were identified within gene pools comparing landrace and cultivar/line germplasm groups (Andean: 18; Mesoamerican: 69) and between the gene pools (59 SNPs), predominantly on chromosomes 1 and 9. The LD extension estimate corrected for population structure and relatedness (r2 SV) was ~ 88 kbp, while for the Andean gene pool was ~ 395 kbp, and for the Mesoamerican was ~ 130 kbp.ConclusionsFor common bean, DArTseq provides an efficient and cost-effective strategy of generating SNPs for large-scale genome-wide studies. The DArTseq resulted in an operational panel of 560 polymorphic SNPs in linkage equilibrium, providing high genome coverage. This SNP set could be used in genotyping platforms with many applications, such as population genetics, phylogeny relation between common bean varieties and support to molecular breeding approaches.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-017-3805-4) contains supplementary material, which is available to authorized users.

New carioca bean cultivars are being introduced into the market necessitating their evaluation under trade conditions, which often require storage under ambient conditions. We therefore evaluated the darkening and hardening processes of six carioca bean genotypes each representing regular and slow darkening trait during storage under ambient conditions for five months to elucidate their relationship as a breeding strategy. Storage time adversely affected color characteristics (L*, a*, b*, C* and ΔE) depending on bean genotype, whereas hardness and resistance to cooking increased during storage independent of the lignification process. Bean darkening and hardening occurred during storage at different intensities in each genotype and were not always correlated. BRSMG-Madrepé rola, a slow darkening genotype, was unaffected (resistant to storage conditions), whereas BRS-Pontal with regular tegument darkening, was highly susceptible to storage conditions reflected in extended cooking time and darkening (low L* values). Principal component and cluster analyses on 8 constituents analyzed in this study demonstrate the difference in color characteristics, cooking time and hardness as major factors in segregating the bean genotypes. Seed coat color is an important but inappropriate single parameter for predicting the resistance to cooking or hardness induced by storage of carioca beans under ambient conditions. Development of carioca bean genotypes resistant to storage conditions is essential in reducing food losses during postharvest.

No abstract

Single-cell transcriptome analysis has been extensively applied in humans and animal models to uncover gene expression heterogeneity between the different cell types of a tissue or an organ. It demonstrated its capability to discover key regulatory elements that determine cell fate during developmental programs. Single-cell analysis requires the isolation and labeling of the messenger RNA (mRNA) derived from each cell. These challenges were primarily addressed in mammals by developing microfluidic-based approaches. For plant species whose cells contain cell walls, these approaches have generally required the generation of isolated protoplasts. Many plant tissues’ secondary cell wall hinders enzymatic digestion required for individual protoplast isolation, resulting in an unequal representation of cell types in a protoplast population. This limitation is especially critical for cell types located in the inner layers of a tissue or the inner tissues of an organ. Consequently, single-cell RNA sequencing (scRNA-seq) studies using microfluidic approaches in plants have mainly been restricted to Arabidopsis roots, for which well-established procedures of protoplast isolation are available. Here we present a simple alternative approach to generating high-quality protoplasts from plant tissue by characterizing the mRNA extracted from individual nuclei instead of whole cells. We developed the protocol using two different plant materials with varying cellular complexity levels and cell wall structure, Populus shoot apices, and more lignified stems. Using the 10× Genomics Chromium technology, we show that this procedure results in intact mRNA isolation and limited leakage, with a broad representation of individual cell transcriptomes.

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.

Drought stress, particularly during the flowering and grain-filling stages of growth, contributes to serious yield loss in common bean (Phaseolus vulgaris L.). The aim of this study was to identify genes induced in response to drought stress using transcriptome analysis of contrasting genotypes. Using leaf tissues of tolerant (BAT 477) and susceptible common bean genotypes (Pérola), collected at the flowering and grain-filling stages, four complementary deoxyribonucleic acid representational difference analysis subtractive libraries were constructed and then sequenced. A total of 7,203 (77.6 %) sequences with an average sequence size of 570 bp were considered valid, for a combined 4 Mbp sequence. According to a differential display analysis, 802 expressed sequence tags, distributed across 67 contigs, were differentially expressed by the tolerant (37 contigs) and susceptible genotypes (30 contigs) after identification under drought conditions during the two investigated plant developmental stages. Of these differential contigs, the 13 most frequent genes were exclusive to the tolerant genotype. Based on BLAST2GO, 73 % of the gene sequences were annotated and 12 % showed mapping results, with the highest similarity rate corresponding to Glycine max (41 %). According to gene ontology functional analysis, 48 % of the sequences were attributed to cell metabolic processes. Overall, 8.3 % of the transcribed sequences exhibited similarity to transcription factors, predominantly those of the AP2-EREBP family (97.8 %). Of the target sequences validated by quantitative real-time polymerase chain reaction, most genes showed an expression level that agreed with that predicted by in silico analysis. Thus, the drought transcriptome dataset is a valuable Electronic supplementary material The online version of this article

Genes related to the response to drought stress in leaf and root tissue of drought-susceptible (DS) and tolerant (DT) genotypes were characterized by RNA-Seq. In total, 54,750 transcripts, representative of 28,590 genes, were identified; of these, 1,648 were of high-fidelity (merge of 12 libraries) and described for the first time in the Andean germplasm. From the 1,239 differentially expressed genes (DEGs), 458 were identified in DT, with a predominance of genes in categories of oxidative stress, response to stimulus and kinase activity. Most genes related to oxidation-reduction terms in roots were early triggered in DT (T75) compared to DS (T150) suggestive of a mechanism of tolerance by reducing the damage from ROS. Among the KEGG enriched by DEGs up-regulated in DT leaves, two related to the formation of Sulfur-containing compounds, which are known for their involvement in tolerance to abiotic stresses, were common to all treatments. Through qPCR, 88.64% of the DEGs were validated. A total of 151,283 variants were identified and functional effects estimated for 85,780. The raw data files were submitted to the NCBI database. A transcriptome map revealed new genes and isoforms under drought. These results supports a better understanding of the drought tolerance mechanisms in beans.

No abstract