Development of a genome-wide anchored microsatellite map for common bean (Phaseolus vulgaris L.)

Mapping of adult plant recessive resistance to anthracnose in Indian common bean landrace Baspa/KRC 8

BACKGROUND

The common bean (Phaseolus vulgaris) has become the food of choice owing to its wealthy nutritional profile, leading to a considerable increase in its cultivation worldwide. However, anthracnose has been a major impediment to production and productivity, as elite bean cultivars are vulnerable to this disease. To overcome barriers in crop production, scientists worldwide are working towards enhancing the genetic diversity of crops. One way to achieve this is by introducing novel genes from related crops, including landraces like KRC 8. This particular landrace, found in the North Western Himalayan region, has shown adult plant resistance against anthracnose and also possesses a recessive resistance gene.

METHODS AND RESULTS

In this study, a population of 179 F2:9 RIL individuals (Jawala × KRC 8) was evaluated at both phenotypic and genotypic levels using over 830 diverse molecular markers to map the resistance gene present in KRC 8. We have successfully mapped a resistance gene to chromosome Pv01 using four SSR markers, namely IAC 238, IAC 235, IAC 259, and BM 146. The marker IAC 238 is closely linked to the gene with a distance of 0.29 cM, while the other markers flank the recessive resistance gene at 10.87 cM (IAC 259), 17.80 cM (BM 146), and 25.22 cM (IAC 235). Previously, a single recessive anthracnose resistance gene (co-8) has been reported in the common bean accession AB 136. However, when we performed PCR amplification with our tightly linked marker IAC 238, we got different amplicons in AB 136 and KRC 8. Interestingly, the susceptible cultivar Jawala produced the same amplicon as AB 136. This observation indicated that the recessive gene present in KRC 8 is different from co-8. As the gene is located far away from the Co-1 locus, we suggest naming the recessive gene co-Indb/co-19. Fine mapping of co-Indb in KRC 8 may provide new insights into the cloning and characterization of this recessive gene so that it can be incorporated into future bean improvement programs. Further, the tightly linked marker IAC 238 can be utilized in marker assisted introgression in future bean breeding programs.

CONCLUSION

The novel co-Indb gene present in Himalayan landrace KRC 8, showing adult plant resistance against common bean anthracnose, is independent from all the resistance genes previously located on chromosome Pv01.

Fine mapping of a new common bean anthracnose resistance gene (Co-18) to the proximal end of Pv10 in Indian landrace KRC-5

KEY MESSAGE

Mapping and fine mapping of bean anthracnose resistance genes is a continuous process. We report fine mapping of anthracnose resistance gene Co-18 which is the first anthracnose gene mapped to Pv10. The discovery of resistance gene is a major gain in the bean anthracnose pathosystem research. Among the Indian common bean landraces, KRC-5 exhibit high levels of resistance to the bean anthracnose pathogen Colletotrichum lindemuthianum. To precisely map the anthracnose resistance gene, we used a Recombinant Inbred Line (F2:9 RIL) population (KRC-5 × Jawala). The inheritance test revealed that KRC-5 carries a dominant resistance gene temporarily designated as Co-18. We discovered two RAPD markers linked to Co-18 among 287 RAPD markers. These RAPD markers were eventually developed into SCARs (Sc-OPR15 and Sc-OPF6) and flank Co-18 on chromosome Pv10 at a distance of 5.3 and 4.2 cM, respectively. At 4.0-4.1 Mb on Pv10, we detected a SNP (single-nucleotide polymorphism) signal. We synthesized 58 SSRs and 83 InDels from a pool of 135 SSRs and 1134 InDels, respectively. Five SSRs, four InDels, and two SCARs were used to generate the high-density linkage map, which led to the identification of two SSRs (SSR24 and SSR36) that are tightly linked to Co-18. These two SSRs flank the Co-18 to 178 kb genomic region with 13 candidate genes including five NLR (nucleotide-binding and leucine-rich repeat) genes. The closely linked markers SSR24 and SSR36 will be used in cloning and pyramiding of the Co-18 gene with other R genes to develop durable resistant bean varieties.

Unique Salt-Tolerance-Related QTLs, Evolved in Vigna riukiuensis (Na+ Includer) and V. nakashimae (Na+ Excluder), Shed Light on the Development of Super-Salt-Tolerant Azuki Bean (V. angularis) Cultivars

Wild relatives of crops have the potential to improve food crops, especially in terms of improving abiotic stress tolerance. Two closely related wild species of the traditional East Asian legume crops, Azuki bean (Vigna angularis), V. riukiuensis "Tojinbaka" and V. nakashimae "Ukushima" were shown to have much higher levels of salt tolerance than azuki beans. To identify the genomic regions responsible for salt tolerance in "Tojinbaka" and "Ukushima", three interspecific hybrids were developed: (A) azuki bean cultivar "Kyoto Dainagon" × "Tojinbaka", (B) "Kyoto Dainagon" × "Ukushima" and (C) "Ukushima" × "Tojinbaka". Linkage maps were developed using SSR or restriction-site-associated DNA markers. There were three QTLs for "percentage of wilt leaves" in populations A, B and C, while populations A and B had three QTLs and population C had two QTLs for "days to wilt". In population C, four QTLs were detected for Na⁺ concentration in the primary leaf. Among the F₂ individuals in population C, 24% showed higher salt tolerance than both wild parents, suggesting that the salt tolerance of azuki beans can be further improved by combining the QTL alleles of the two wild relatives. The marker information would facilitate the transfer of salt tolerance alleles from "Tojinbaka" and "Ukushima" to azuki beans.

QTL mapping and identification of genes associated with the resistance to Acanthoscelides obtectus in cultivated common bean using a high-density genetic linkage map

BACKGROUND

Common bean (Phaseolus vulgaris L.) is an important agricultural product with large nutritional value, and the insect pest Acanthoscelides obtectus (Say) seriously affects its product quality and commodity quality during storage. Few researches on genes of bruchid resistance have investigated in common bean cultivars.

RESULTS

In this study, a bruchid-resistant cultivar black kidney bean and a highly susceptible accession Longyundou3 from different gene banks were crossed to construct a recombinant inbred line population. The genetic analysis indicated a quantitative inheritance of the bruchid resistance trait controlled by polygenes. A high-density genetic map of a total map distance of 1283.68 cM with an average interval of 0.61 cM between each marker was constructed using an F₆ population of 157 recombinant inbred lines. The map has 3106 bin markers, containing 2,234,769 SNPs. Using the high-density genetic map, a new quantitative trait locus for the resistance to Acanthoscelides obtectus was identified on chromosome 6. New molecular markers based on the candidate region were developed, and this locus was further delimited to an interval of 122.3 kb between SSR markers I6-4 and I6-16 using an F₂ population. This region comprised five genes. Phvul.006G003700, which encodes a bifunctional inhibitor, may be a potential candidate gene for bruchid resistance. Sequencing analysis of candidate gene identified a 5 bp insertion-deletion in promoter of gene Phvul.006G003700 between two parents. Expression analysis of candidate gene revealed that the expression level of Phvul.006G003700 in bruchid-resistant parent was markedly higher than that in bruchid-susceptible parent both in dry seeds and leaves.

CONCLUSIONS

A high-density genetic linkage map was constructed utilizing whole-genome resequencing and one new QTL for bruchid resistance was identified on chromosome 6 in common bean cultivar. Phvul.006G003700 (encoding a bifunctional inhibitor) may be a potential candidate gene. These results may form the basis for further research to reveal the bruchid resistance molecular mechanism of common bean.

Insights into the genetic diversity of an underutilized Indian legume, Vigna stipulacea (Lam.) Kuntz., using morphological traits and microsatellite markers

Vigna stipulacea (Lam.) Kuntz., commonly known as Minni payaru is an underutilized legume species and has a great potential to be utilized as food crop. To evaluate and select the best germplasm to be harnessed in the breeding programme, we assessed the genetic diversity of V. stipulacea (94 accessions) conserved in the Indian National Genebank, based on morphological traits and microsatellite markers. Significant variation was recorded for the morphological traits studied. Euclidean distance using UPGMA method grouped all accessions into two major clusters. Accessions were identified for key agronomic traits such as, early flowering (IC331436, IC251436, IC331437); long peduncle length (IC553518, IC550531, IC553557, IC553540, IC550532, IC553564); and more number of seeds per pod (IC553529, IC622865, IC622867, IC553528). To analyse the genetic diversity among the germplasm 33 SSR primers were used anda total of 116 alleles were detected. The number of alleles varied from two to seven, with an average of 3.52 per loci. The polymorphic information content values varied from 0.20 to 0.74, with a mean of 0.40. The high number of alleles per locus and the allelic diversity in the studied germplasm indicated a relatively wider genetic base of V. stipulacea. Phylogenetic analysis clustered accessions into seven clades. Population structure analysis grouped them into five genetic groups, which were partly supported by PCoA and phylogenetic tree. Besides, PCoA and AMOVA also decoded high genetic diversity among the V. stipulacea accessions. Thus, morphological and microsatellite markers distinguished V. stipulacea accessions and assessed their genetic diversity efficiently. The identified promising accessions can be utilized in Vigna improvement programme through introgression breeding and/or can be used for domestication and enhanced utilization of V. stipulacea.

Genetic diversity and population genetic structure analysis of an extensive collection of wild and cultivated Vigna accessions

Vigna is a large, pan-tropic and highly variable group of the legumes family which is known for its > 10 cultivated species having significant commercial value for their nutritious grains and multifarious uses. The wild vignas are considered a reservoir of numerous useful traits which can be deployed for introgression of resistance to biotic and abiotic stresses, seed quality and enhanced survival capability in extreme environments. Nonetheless, for their effective utilization through introgression breeding information on their genetic diversity, population structure and crossability is imperative. Keeping this in view, the present experiment was undertaken with 119 accessions including 99 wild Vigna accessions belonging to 19 species and 18 cultivated genotypes of Vigna and 2 of Phaseolus. Total 102 polymorphic SSRs were deployed to characterize the material at molecular level which produced 1758 alleles. The genotypes were grouped into four major clusters which were further sub-divided in nine sub-clusters. Interestingly, all cultivated species shared a single cluster while no such similarities were observed for the wild accessions as these were distributed in different groups of sub-clusters. The co-dominant allelic data of 114 accessions were then utilized for obtaining status of the accessions and their hybrid forms. The model-based population structure analysis categorized 114 accessions of Vigna into 6 genetically distinct sub-populations (K = 6) following admixture-model based simulation with varying levels of admixture. 91 (79.82%) accessions resembled their hierarchy and 23 (20.18%) accessions were observed as the admixture forms. Maximum number of accessions (25) were grouped in sub-population (SP) 6 and the least accessions were grouped in SP3 and SP5 (11 each). The population genetic structure, therefore, supported genetic diversity analysis and provided an insight into the genetic lineage of these species which will help in effective use of germplasm for development of cultivars following selective prebreeding activities.

Angular Leaf Spot Resistance Loci Associated With Different Plant Growth Stages in Common Bean

Angular leaf spot (ALS) is a disease that causes major yield losses in the common bean crop. Studies based on different isolates and populations have already been carried out to elucidate the genetic mechanisms of resistance to ALS. However, understanding of the interaction of this resistance with the reproductive stages of common bean is lacking. The aim of the present study was to identify ALS resistance loci at different plant growth stages (PGS) by association and linkage mapping approaches. An BC₂F₃ inter-gene pool cross population (AND 277 × IAC-Milênio - AM population) profiled with 1,091 SNPs from genotyping by sequencing (GBS) was used for linkage mapping, and a carioca diversity panel (CDP) genotyped by 5,398 SNPs from BeadChip assay technology was used for association mapping. Both populations were evaluated for ALS resistance at the V2 and V3 PGSs (controlled conditions) and R8 PGS (field conditions). Different QTL (quantitative trait loci) were detected for the three PGSs and both populations, showing a different quantitative profile of the disease at different plant growth stages. For the three PGS, multiple interval mapping (MIM) identified seven significant QTL, and the Genome-wide association study (GWAS) identified fourteen associate SNPs. Several loci validated regions of previous studies, and Phg-1, Phg-2, Phg-4, and Phg-5, among the 5 loci of greatest effects reported in the literature, were detected in the CDP. The AND 277 cultivar contained both the Phg-1 and the Phg-5 QTL, which is reported for the first time in the descendant cultivar CAL143 as ALS10.1^UC. The novel QTL named ALS11.1^AM was located at the beginning of chromosome Pv11. Gene annotation revealed several putative resistance genes involved in the ALS response at the three PGSs, and with the markers and loci identified, new specific molecular markers can be developed, representing a powerful tool for common bean crop improvement and for gain in ALS resistance.

Assessment of the Origin and Diversity of Croatian Common Bean Germplasm Using Phaseolin Type, SSR and SNP Markers and Morphological Traits

Landraces represent valuable genetic resources for breeding programmes to produce high-yielding varieties adapted to stressful environmental conditions. Although the common bean (Phaseolus vulgaris L.) is an economically important food legume for direct human consumption worldwide, common bean production in Croatia is based almost exclusively on landraces and there is no common bean breeding program. Information on phaseolin type and results of population structure and genetic diversity obtained by analysis of SSR and SNP markers, in combination with the morphological characterization of 174 accessions of 10 common bean landraces (morphotypes), enabled thorough classification of accessions. The accessions were classified into phaseolin type H1 (“S”) of Mesoamerican origin and phaseolin types H2 (“H” or “C”) and H3 (“T”) of Andean origin. By applying distance- and model-based clustering methods to SSR markers, the accessions were classified into two clusters at K = 2 separating the accessions according to the centres of origin, while at K = 3, the accessions of Andean origin were further classified into two clusters of accessions that differed in phaseolin type (H2 and H3). Using SNP markers, model-based analysis of population structure was performed, the results of which were consistent with those of SSR markers. In addition, 122 accessions were assigned to 14 newly formed true-type morphogenetic groups derived from three different domestication events: (1) Mesoamerican (H1A) (“Biser”, “Kukuruzar”, “Tetovac”, “Trešnjevac”), (2) Andean—indeterminate type (H2B1) (“Dan noć”, “Sivi”, “Puter”, ”Sivi prošarani”, “Trešnjevac”) and (3) Andean—determinate type (H3B2) (“Bijeli”, “Dan noć”, “Puter”, “Trešnjevac”, “Zelenčec”). The rest of the accessions could represent putative hybrids between morphogenetic groups. The differences between the true-type groups of accessions were further analysed based on nine quantitative traits, and the subsets of traits that best distinguish among centres of origin (A: Mesoamerican, B: Andean) and genetic groups (H1A, H2B1, H3B2) were proposed.

Marker association study of yield attributing traits in common bean (Phaseolus vulgaris L.)

Common bean is gaining acceptance as one of the most valuable major food consumed worldwide owing to innumerable nutritional and therapeutic benefits. Comparatively less productivity in underdeveloped countries encouraged us to proceed for QTL mining of yield traits in common bean. Heretofore, multiple yield associated markers have been detected all over the world; even so, the present work is looked on as the first report on identification of novel/new potent markers by exploiting the germplasm of Northern India. A panel of one hundred and thirty five genotypes was used for morphological studies and based on preliminary molecular evaluation; a set of ninety six diverse common bean genotypes (core set) was selected for association analysis. Molecular data generated by a total of ninety eight microsatellite markers (53 genomic and 45 genic SSRs) revealed high estimation of polymorphism among the genotypes that were observed to be divided into two major sub-populations and varying levels of admixtures based on population structure analyses. By employing both MLM and GLM analysis approaches, we identified 46 and 16 significant marker-trait associations (p ≤ 0.005) respectively, few of which have already been reported and hence validate our results. PVBR213 marker was found to be strongly associated with days to bud initiation trait when analyzed with both the approaches. Phenotypic variation of identified significant markers ranged from 3.1% to 32.7% where PVBR87, PVBR213, X96999 and X57022 explain more than 30% of phenotypic variation for 100 seed weight, days to bud initiation, pods per plant and pod length traits respectively. These findings introduce highly informative markers to aid marker-assisted selection program in common bean for high yield performance along with good agronomic merit.

Current Perspectives on Introgression Breeding in Food Legumes

Food legumes are important for defeating malnutrition and sustaining agri-food systems globally. Breeding efforts in legume crops have been largely confined to the exploitation of genetic variation available within the primary genepool, resulting in narrow genetic base. Introgression as a breeding scheme has been remarkably successful for an array of inheritance and molecular studies in food legumes. Crop wild relatives (CWRs), landraces, and exotic germplasm offer great potential for introgression of novel variation not only to widen the genetic base of the elite genepool for continuous incremental gains over breeding cycles but also to discover the cryptic genetic variation hitherto unexpressed. CWRs also harbor positive quantitative trait loci (QTLs) for improving agronomic traits. However, for transferring polygenic traits, "specialized population concept" has been advocated for transferring QTLs from CWR into elite backgrounds. Recently, introgression breeding has been successful in developing improved cultivars in chickpea (Cicer arietinum), pigeonpea (Cajanus cajan), peanut (Arachis hypogaea), lentil (Lens culinaris), mungbean (Vigna radiata), urdbean (Vigna mungo), and common bean (Phaseolus vulgaris). Successful examples indicated that the usable genetic variation could be exploited by unleashing new gene recombination and hidden variability even in late filial generations. In mungbean alone, distant hybridization has been deployed to develop seven improved commercial cultivars, whereas in urdbean, three such cultivars have been reported. Similarly, in chickpea, three superior cultivars have been developed from crosses between C. arietinum and Cicer reticulatum. Pigeonpea has benefited the most where different cytoplasmic male sterility genes have been transferred from CWRs, whereas a number of disease-resistant germplasm have also been developed in Phaseolus. As vertical gene transfer has resulted in most of the useful gene introgressions of practical importance in food legumes, the horizontal gene transfer through transgenic technology, somatic hybridization, and, more recently, intragenesis also offer promise. The gains through introgression breeding are significant and underline the need of bringing it in the purview of mainstream breeding while deploying tools and techniques to increase the recombination rate in wide crosses and reduce the linkage drag. The resurgence of interest in introgression breeding needs to be capitalized for development of commercial food legume cultivars.

Genetic Architecture and Genomic Prediction of Cooking Time in Common Bean (Phaseolus vulgaris L.)

Cooking time of the common bean is an important trait for consumer preference, with implications for nutrition, health, and environment. For efficient germplasm improvement, breeders need more information on the genetics to identify fast cooking sources with good agronomic properties and molecular breeding tools. In this study, we investigated a broad genetic variation among tropical germplasm from both Andean and Mesoamerican genepools. Four populations were evaluated for cooking time (CKT), water absorption capacity (WAC), and seed weight (SdW): a bi-parental RIL population (DxG), an eight-parental Mesoamerican MAGIC population, an Andean (VEF), and a Mesoamerican (MIP) breeding line panel. A total of 922 lines were evaluated in this study. Significant genetic variation was found in all populations with high heritabilities, ranging from 0.64 to 0.89 for CKT. CKT was related to the color of the seed coat, with the white colored seeds being the ones that cooked the fastest. Marker trait associations were investigated by QTL analysis and GWAS, resulting in the identification of 10 QTL. In populations with Andean germplasm, an inverse correlation of CKT and WAC, and also a QTL on Pv03 that inversely controls CKT and WAC (CKT3.2/WAC3.1) were observed. WAC7.1 was found in both Mesoamerican populations. QTL only explained a small part of the variance, and phenotypic distributions support a more quantitative mode of inheritance. For this reason, we evaluated how genomic prediction (GP) models can capture the genetic variation. GP accuracies for CKT varied, ranging from good results for the MAGIC population (0.55) to lower accuracies in the MIP panel (0.22). The phenotypic characterization of parental material will allow for the cooking time trait to be implemented in the active germplasm improvement programs. Molecular breeding tools can be developed to employ marker-assisted selection or genomic selection, which looks to be a promising tool in some populations to increase the efficiency of breeding activities.

Diversification and genetic structure of the western-to-eastern progression of European Phaseolus vulgaris L. germplasm

BACKGROUND

Common bean (Phaseolus vulgaris L.) is the most important food legume for direct human consumption around the world, as it represents a valuable source of components with nutritional and health benefits.

RESULTS

We conducted a study to define and explain the genetic relatedness and diversification level of common bean (Phaseolus vulgaris L.) germplasm from Portugal to Ukraine, along a western-to-eastern line of southern European countries, including Poland. This was based on the P. vulgaris genetic structure, and was designed to better describe its distribution and domestication pathways in Europe. Using the multi-crop passport descriptors that include geographic origin and different phaseolin types (corresponding to the Mesoamerican and Andean gene pools), 782 accessions were obtained from nine gene banks and 12 geographic origins. We selected 33 genome/ gene-related/ gene-pool-related nuclear simple sequence repeat markers that covered the genetic diversity across the P. vulgaris genome. The overall polymorphic information content was 0.800. Without specifying geographic origin, global structure cluster analysis generated 10 genetic clusters. Among the PvSHP1 markers, the most informative for gene pool assignment of the European P. vulgaris germplasm was PvSHP1-B. Results of AMOVA show that 89% of the molecular variability is shared within the 782 accessions, with 4% molecular variability among the different geographic origins along this western-to-eastern line of southern Europe (including Poland).

CONCLUSIONS

This study shows that the diversification line of the European P. vulgaris germplasm followed from the western areas of southern Europe (Portugal, Spain, Italy, Slovenia) to the more eastern areas of southern Europe. This progression defines three geographically separated subgroups, as the northern (Poland, Ukraine, Romania), southern (Albania, Bulgaria), and central (Bosnia and Herzegovina, Serbia, Hungary) areas of eastern Europe.

Rice bean: a lesser known pulse with well-recognized potential

Required genetic resources for the improvement of agronomic, nutritional and economic value of rice bean are available in the world collection. International cooperative effort is required to utilize and conserve them. Rice bean [Vigna umbellata (Thunb.) Ohwi and Ohashi], a lesser known pulse among the Asiatic Vigna, has long been considered as a food security crop of small and marginal farmers of Southeast Asia. Considered as a nutritionally rich food and fodder, it is also a source of genes for biotic and abiotic stress tolerance including drought, soil acidity and storage pest. Although it spread from its centre of domestication in the Indo-China region to other parts around the world, it never became an important crop anywhere probably because of agronomic disadvantages. Crop improvement for determinate nature, good yield, less variable seed colour, pleasant organoleptic properties and lower antinutrients is required. Scanning of scientific literature indicates that genetic resources with desirable agronomic and nutritional traits exist within the current collection but are spread across countries. Genomic studies in the species indicate that except for insect resistance and aluminium toxicity tolerance, not much attention has been paid to decipher and utilize other stress tolerance and nutritional quality traits. Collaborative efforts towards improving farming, food, trade value and off-farm conservation of rice bean would not only help marginal farmers but will also help to preserve the yet to be explored genomic resources available in this sturdy pulse.

Construction of genetic linkage map and genome dissection of domestication-related traits of moth bean (Vigna aconitifolia), a legume crop of arid areas

The moth bean (Vigna aconitifolia), possibly the most primitive crop of the genus Vigna, is a highly drought- and heat-resistant legume grown in arid areas. Moth bean domestication involved phenotypic changes, including reduction of seed dormancy and pod shattering, increased organ size, and earlier flowering and maturity. However, the genetics of the domestication process in moth bean is not known. In this study, we constructed a genetic linkage map for moth bean and used the map to identify quantitative trait loci (QTL) for domestication-related traits of an F₂ population of 188 individuals produced from a cross of wild moth bean (TN67) and cultivated moth bean (ICPMO056). The genetic linkage map comprised 11 linkage groups (LG) of 172 simple sequence repeat markers and spanned a total length of 1016.8 centiMorgan (cM), with an average marker distance of 7.34 cM. A comparative genome analysis showed high genome synteny between moth bean and mungbean (Vigna radiata), adzuki bean (Vigna angularis), rice bean (Vigna umbellata), and yardlong bean (Vigna unguiculata). In total, 50 QTLs and 3 genes associated with 20 domestication-related traits were identified. Most of the QTLs belonged to five LGs (1, 2, 4, 7, and 10). Key traits related to domestication such as seed dormancy and pod shattering were controlled by large-effect QTLs (PVE > 20%) with one or two minor QTLs, whereas all other traits were controlled by one-seven minor QTLs, apart from seed weight, which was controlled by one major and seven minor QTLs. These results suggest that a small number of mutations with large phenotypic effects have contributed to the domestication of the moth bean. Comparative analysis of QTLs with related Vigna crops revealed that there are several domestication-related large-effect QTLs that had not been used in moth bean domestication. This study provides a basic genetic map and identified genome regions associated with domestication-related traits, which will be useful for the genetic improvement of the moth bean and related Vigna species.

Evolution of SSR diversity from wild types to U.S. advanced cultivars in the Andean and Mesoamerican domestications of common bean (Phaseolus vulgaris)

Progress in common bean breeding requires the exploitation of genetic variation among market classes, races and gene pools. The present study was conducted to determine the amount of genetic variation and the degree of relatedness among 192 selected common bean advanced cultivars using 58 simple-sequence-repeat markers (SSR) evenly distributed along the 11 linkage groups of the Phaseolus reference map. All the lines belonged to commercial seed type classes that are widely grown in the USA and include both dry bean and snap beans for the fresh and processing markets. Through population structure, principal components analyses, cluster analysis, and discriminant analysis of principal components (DAPC), Andean and Mesoamerican genotypes as well as most American commercial type classes could be distinguished. The genetic relationship among the commercial cultivars revealed by the SSR markers was generally in agreement with known pedigree data. The Mesoamerican cultivars were separated into three major groups-black, small white, and navy accessions clustered together in a distinct group, while great northern and pinto clustered in another group, showing mixed origin. The Andean cultivars were distributed in two different groups. The kidney market classes formed a single group, while the green bean accessions were distributed between the Andean and Mesoamerican groups, showing inter-gene pool genetic admixture. For a subset of 24 SSR markers, we compared and contrasted the genetic diversity of the commercial cultivars with those of wild and domesticated landrace accessions of common bean. An overall reduction in genetic diversity was observed in both gene pools, Andean and Mesoamerican, from wild to landraces to advanced cultivars. The limited diversity in the commercial cultivars suggests that an important goal of bean breeding programs should be to broaden the cultivated gene pool, particularly the genetic diversity of specific commercial classes, using the genetic variability present in common bean landraces.

Evidence for Strong Kinship Influence on the Extent of Linkage Disequilibrium in Cultivated Common Beans

Phaseolus vulgaris is an important grain legume for human consumption. Recently, association mapping studies have been performed for the species aiming to identify loci underlying quantitative variation of traits. It is now imperative to know whether the linkage disequilibrium (LD) reflects the true association between a marker and causative loci. The aim of this study was to estimate and analyze LD on a diversity panel of common beans using ordinary r² and r2 extensions which correct bias due to population structure (rS²), kinship (rV²), and both (rVS²). A total of 10,362 single nucleotide polymorphisms (SNPs) were identified by genotyping by sequencing (GBS), and polymorphisms were found to be widely distributed along the 11 chromosomes. In terms of r2, high values of LD (over 0.8) were identified between SNPs located at opposite chromosomal ends. Estimates for rV² were lower than those for rS². Results for rV² and rVS² were similar, suggesting that kinship may also include information on population structure. Over genetic distance, LD decayed to 0.1 at a distance of 1 Mb for rVS². Inter-chromosomal LD was also evidenced. This study showed that LD estimates decay dramatically according to the population structure, and especially the degree of kinship. Importantly, the LD estimates reported herein may influence our ability to perform association mapping studies on P. vulgaris.

QTL analysis of domestication syndrome in zombi pea (Vigna vexillata), an underutilized legume crop

Zombi pea (Vigna vexillata (L.) A. Rich) is an underutilized crop belonging to the genus Vigna. Two domesticated forms of zombi pea are cultivated as crop plants; seed and tuber forms. The cultivated seed form is present in Africa, while the cultivated tuber form is present in a very limited part of Asia. Genetics of domestication have been investigated in most of cultivated Vigna crops by means of quantitative trait locus (QTL) mapping. In this study, we investigated genetics of domestication in zombi pea by QTL analysis using an F2 population of 139 plants derived from a cross between cultivated tuber form of V. vexillata (JP235863) and wild V. vexillata (AusTRCF66514). A linkage map with 11 linkage groups (LGs) was constructed from this F2 population using 145 SSR, 117 RAD-seq and 2 morphological markers. Many highly segregation distorted markers were found on LGs 5, 6, 7, 8, 10 and 11. Most of the distorted markers were clustered together and all the markers on LG8 were highly distorted markers. Comparing this V. vexillata linkage map with linkage maps of other four Vigna species demonstrated several genome rearrangements in V. vexillata. QTL analysis for 22 domestication-related traits was investigated by inclusive composite interval mapping in which 37 QTLs were identified for 18 traits; no QTL was detected for 4 traits. Number of QTLs detected in each trait ranged from 1 to 5 with an average of only 2.3. Five QTLs for tuber width and three QTLs for tuber weight. Interestingly, 2 QTLs each for tuber width and tuber weight detected on LG2 and LG4 were located at similar position and wild allele increased tuber width and weight. This indicated wild germplasm having small tuber have potential to increase yield of large tuber cultivated type. Large-effect QTLs (PVE > 20%) were on LG4 (pod length), LG5 (leaf size and seed thickness), and LG7 (for seed-related traits). Comparison of domestication-related QTLs of the zombi pea with those of cowpea (Vigna unguiculata), azuki bean (Vigna angularis), mungbean (Vigna radiata) and rice bean (Vigna umbellata) revealed that there was conservation of some QTLs for seed size, pod size and leaf size between zombi pea and cowpea and that QTLs associated with seed size (weight, length, width and thickness) in each species were clustered on same linkage.

QTL analyses for tolerance to abiotic stresses in a common bean (Phaseolus vulgaris L.) population

Common bean productivity is reduced by several abiotic stress factors like drought and low soil fertility, leading to yield losses particularly in low input smallholder farming systems in the tropics. To understand the genetics of stress tolerance, and to improve adaptation of common bean to adverse environments, the BAT 881 x G21212 population of 95 recombinant inbred lines (RILs) was evaluated under different abiotic stress conditions in 15 trials across four locations in Colombia, representing two higher altitude (Darién, Popayán) and two lower altitude (Palmira, Quilichao) locations. Stress vs non-stress treatments showed that yields were reduced in drought trials in Palmira by 13 and 31%, respectively, and observed yield reductions in low phosphorus stress were 39% in Quilichao, 16% in Popayán, and 71% in Darién, respectively. Yield components and biomass traits were also reduced. Traits linked to dry matter redistribution from stems, leaves and pods to seed, such as pod harvest index and total non-structural carbohydrates, were found to be important factors contributing to yield in all conditions. In contrast, early maturity was correlated with improved yield only in lower altitude locations, whereas in higher altitudes delayed maturity promoted yield. Superior RILs that combine stress tolerance and high cross-location productivity were identified. Lines that showed good yield under strong stress conditions also performed well under non-stress conditions, indicating that breeder's selection can be applied for both conditions at the same time. Quantitative trait loci (QTL) analyses revealed a stable yield QTL on chromosome Pv04, detected individually in all locations, several stress treatments and in best linear unbiased predictions (BLUPs) across all trials. Furthermore, two QTL hotspots for maturity traits were identified on Pv01 and Pv08, which are the most stable QTL. The constitutive yield QTL could serve as a good candidate for marker development and could be used in marker assisted selection. Increased understanding of the physiology of abiotic stress tolerance, combined with the availability of superior germplasm and molecular tools, will aid breeding efforts for further improvement of these plant traits.

Uneven recombination rate and linkage disequilibrium across a reference SNP map for common bean (Phaseolus vulgaris L.)

Recombination (R) rate and linkage disequilibrium (LD) analyses are the basis for plant breeding. These vary by breeding system, by generation of inbreeding or outcrossing and by region in the chromosome. Common bean (Phaseolus vulgaris L.) is a favored food legume with a small sequenced genome (514 Mb) and n = 11 chromosomes. The goal of this study was to describe R and LD in the common bean genome using a 768-marker array of single nucleotide polymorphisms (SNP) based on Trans-legume Orthologous Group (TOG) genes along with an advanced-generation Recombinant Inbred Line reference mapping population (BAT93 x Jalo EEP558) and an internationally available diversity panel. A whole genome genetic map was created that covered all eleven linkage groups (LG). The LGs were linked to the physical map by sequence data of the TOGs compared to each chromosome sequence of common bean. The genetic map length in total was smaller than for previous maps reflecting the precision of allele calling and mapping with SNP technology as well as the use of gene-based markers. A total of 91.4% of TOG markers had singleton hits with annotated Pv genes and all mapped outside of regions of resistance gene clusters. LD levels were found to be stronger within the Mesoamerican genepool and decay more rapidly within the Andean genepool. The recombination rate across the genome was 2.13 cM / Mb but R was found to be highly repressed around centromeres and frequent outside peri-centromeric regions. These results have important implications for association and genetic mapping or crop improvement in common bean.

Gene/QTL discovery for Anthracnose in common bean (Phaseolus vulgaris L.) from North-western Himalayas

Common bean (Phaseolus vulgaris L.) is one of the most important grain legume crops in the world. The beans grown in north-western Himalayas possess huge diversity for seed color, shape and size but are mostly susceptible to Anthracnose disease caused by seed born fungus Colletotrichum lindemuthianum. Dozens of QTLs/genes have been already identified for this disease in common bean world-wide. However, this is the first report of gene/QTL discovery for Anthracnose using bean germplasm from north-western Himalayas of state Jammu & Kashmir, India. A core set of 96 bean lines comprising 54 indigenous local landraces from 11 hot-spots and 42 exotic lines from 10 different countries were phenotyped at two locations (SKUAST-Jammu and Bhaderwah, Jammu) for Anthracnose resistance. The core set was also genotyped with genome-wide (91) random and trait linked SSR markers. The study of marker-trait associations (MTAs) led to the identification of 10 QTLs/genes for Anthracnose resistance. Among the 10 QTLs/genes identified, two MTAs are stable (BM45 & BM211), two MTAs (PVctt1 & BM211) are major explaining more than 20% phenotypic variation for Anthracnose and one MTA (BM211) is both stable and major. Six (06) genomic regions are reported for the first time, while as four (04) genomic regions validated the already known QTL/gene regions/clusters for Anthracnose. The major, stable and validated markers reported during the present study associated with Anthracnose resistance will prove useful in common bean molecular breeding programs aimed at enhancing Anthracnose resistance of local bean landraces grown in north-western Himalayas of state Jammu and Kashmir.

Distinctness of Brazilian common bean cultivars with carioca and black grain by means of morphoagronomic and molecular descriptors

Common bean (Phaseolus vulgaris L.) is one of the most important legumes for human consumption and is a staple food in the diet of the population of some countries of Latin America, Africa and Asia. The distinction between cultivars is based predominantly on morphological descriptors, which proved inefficient for the differentiation of some cultivars. This study had the objective of describing, distinguishing and evaluating the agronomic potential of 39 common bean cultivars of the carioca and black grain groups registered for cultivation in Brazil, based on 49 morphoagronomic descriptors and microsatellite (simple sequence repeat -SSR) markers. The morphoagronomic traits of each commercial group were characterized in four environments. Thirty-seven SSR markers were used for the molecular description. The morphological data, analyzed by the Shanonon-Weaver index, detected low variability among cultivars for qualitative data. On the other hand, the estimates of variance analysis, relative importance of the traits and hierarchical grouping analysis applied to the quantitative variables showed that the descriptors related to plant morphology were the most important for the carioca group, and those related to seed morphology were determining for the black group. The genetic parameters estimated for SSR markers by hierarchical and Bayesian cluster analysis identified 116 alleles, with 33 and 30 polymorphic loci and 24 and 22 private alleles for the carioca and black groups, respectively. The combined use of morphoagronomic and molecular descriptors improves the distinguishability of cultivars, contributing in a more efficient way to breeding and to the protection of cultivars.

Dissection of Resistance Genes to Pseudomonas syringae pv. phaseolicola in UI3 Common Bean Cultivar

Few quantitative trait loci have been mapped for resistance to Pseudomonas syringae pv. phaseolicola in common bean. Two F₂ populations were developed from the host differential UI3 cultivar. The objective of this study was to further characterize the resistance to races 1, 5, 7 and 9 of Psp included in UI3. Using a QTL mapping approach, 16 and 11 main-effect QTLs for pod and primary leaf resistance were located on LG10, explaining up to 90% and 26% of the phenotypic variation, respectively. The homologous genomic region corresponding to primary leaf resistance QTLs detected tested positive for the presence of resistance-associated gene cluster encoding nucleotide-binding and leucine-rich repeat (NL), Natural Resistance Associated Macrophage (NRAMP) and Pentatricopeptide Repeat family (PPR) proteins. It is worth noting that the main effect QTLs for resistance in pod were located inside a 3.5 Mb genomic region that included the Phvul.010G021200 gene, which encodes a protein that has the highest sequence similarity to the RIN4 gene of Arabidopsis, and can be considered an important candidate gene for the organ-specific QTLs identified here. These results support that resistance to Psp from UI3 might result from the immune response activated by combinations of R proteins, and suggest the guard model as an important mechanism in pod resistance to halo blight. The candidate genes identified here warrant functional studies that will help in characterizing the actual defense gene(s) in UI3 genotype.

Establishing the Bases for Introducing the Unexplored Portuguese Common Bean Germplasm into the Breeding World

Common bean (Phaseolus vulgaris L.) is among the most important grain legumes for human consumption worldwide. Portugal has a potentially promising common bean germplasm, resulting from more than five centuries of natural adaptation and farmers' selection. Nevertheless, limited characterization of this resource hampers its exploitation by breeding programs. To support a more efficient conservation of the national bean germplasm and promote its use in crop improvement, we performed, for the first time, a simultaneous molecular marker (21 microsatellites and a DNA marker for phaseolin-type diversity analysis) and seed and plant morphological characterization (14 traits) of 175 accessions from Portuguese mainland and islands traditional bean-growing regions. A total of 188 different alleles were identified and an average pairwise Cavalli-Sforza and Edwards' chord genetic distance of 0.193 was estimated among accessions. To relate the Portuguese germplasm with the global common bean diversity, 17 wild relatives and representative accessions from the Andean and Mesoamerican gene pools were evaluated at the molecular level. No correlation was detected between the variability found and the geographic origin of accessions. Structure analysis divided the collection into three main clusters. Most of the Portuguese accessions grouped with the race representatives and wild relatives from the Andean region. One third of the national germplasm had admixed genetic origin and might represent putative hybrids among gene pools from the two original centers of domestication in the Andes and Mesoamerica. The molecular marker-based classification was largely congruent with the three most frequent phaseolin haplotype patterns observed in the accessions analyzed. Seed and plant morphological characterization of 150 Portuguese common bean accessions revealed a clear separation among genetic structure and phaseolin haplotype groups of accessions, with seed size and shape and the number of locules per pod the most discriminant traits. Additionally, we used molecular and morphological data to develop a series of smaller core collections that, by maximizing the genetic and morphological diversity of the original collection, represents the Portuguese common bean germplasm with minimum repetitiveness. A core collection with 37 accessions contained 100% of the genetic variation found in the entire collection. This core collection is appropriate for a more detailed characterization and should be explored, as a priority, in national and international common bean breeding efforts. Furthermore, the identified intermediate accessions (with admixed genetic origin) may have novel genetic combinations useful in future bean breeding.

Short-Term Local Adaptation of Historical Common Bean (Phaseolus vulgaris L.) Varieties and Implications for In Situ Management of Bean Diversity

Recognizing both the stakes of traditional European common bean diversity and the role farmers’ and gardeners’ networks play in maintaining this diversity, the present study examines the role that local adaptation plays for the management of common bean diversity in situ. To the purpose, four historical bean varieties and one modern control were multiplied on two organic farms for three growing seasons. The fifteen resulting populations, the initial ones and two populations of each variety obtained after the three years of multiplication, were then grown in a common garden. Twenty-two Simple Sequence Repeat (SSR) markers and 13 phenotypic traits were assessed. In total, 68.2% of tested markers were polymorphic and a total of 66 different alleles were identified. F_ST analysis showed that the genetic composition of two varieties multiplied in different environments changed. At the phenotypic level, differences were observed in flowering date and leaf length. Results indicate that three years of multiplication suffice for local adaptation to occur. The spatial dynamics of genetic and phenotypic bean diversity imply that the maintenance of diversity should be considered at the scale of the network, rather than individual farms and gardens. The microevolution of bean populations within networks of gardens and farms emerges as a research perspective.

De novo transcriptome analysis and microsatellite marker development for population genetic study of a serious insect pest, Rhopalosiphum padi (L.) (Hemiptera: Aphididae)

The bird cherry-oat aphid, Rhopalosiphum padi (L.), is one of the most abundant aphid pests of cereals and has a global distribution. Next-generation sequencing (NGS) is a rapid and efficient method for developing molecular markers. However, transcriptomic and genomic resources of R. padi have not been investigated. In this study, we used transcriptome information obtained by RNA-Seq to develop polymorphic microsatellites for investigating population genetics in this species. The transcriptome of R. padi was sequenced on an Illumina HiSeq 2000 platform. A total of 114.4 million raw reads with a GC content of 40.03% was generated. The raw reads were cleaned and assembled into 29,467 unigenes with an N50 length of 1,580 bp. Using several public databases, 82.47% of these unigenes were annotated. Of the annotated unigenes, 8,022 were assigned to COG pathways, 9,895 were assigned to GO pathways, and 14,586 were mapped to 257 KEGG pathways. A total of 7,936 potential microsatellites were identified in 5,564 unigenes, 60 of which were selected randomly and amplified using specific primer pairs. Fourteen loci were found to be polymorphic in the four R. padi populations. The transcriptomic data presented herein will facilitate gene discovery, gene analyses, and development of molecular markers for future studies of R. padi and other closely related aphid species.

Identification of novel major and minor QTLs associated with Xanthomonas oryzae pv. oryzae (African strains) resistance in rice (Oryza sativa L.)

BACKGROUND

Xanthomonas oryzae pv. oryzae (Xoo) is the causal agent of Bacterial Leaf Blight (BB), an emerging disease in rice in West-Africa which can induce up to 50 % of yield losses. So far, no specific resistance gene or QTL to African Xoo were mapped. The objectives of this study were to identify and map novels and specific resistance QTLs to African Xoo strains.

RESULTS

The reference recombinant inbred lines (RIL) mapping population derived from the cross between IR64 and Azucena was used to investigate Xoo resistance. Resistance to African and Philippine Xoo strains representing different races was assessed on the RIL population under greenhouse conditions. Five major quantitative trait loci (QTL) for resistance against African Xoo were located on different chromosomes. Loci on chromosomes 1, 7, 9, 10 and 11 explained as much as 13 %, 37 %, 13 %, 11 % and 15 % of resistance variation, respectively. A major novel QTL located on chromosome 7 explained 37 % of the phenotypic variance to the African Xoo corresponding to race A3 whereas that on chromosome 11 is effective to all African races tested. Together with genes and QTLs for resistance to bacterial blight previously described, the QTLs described here were mapped onto the reference O. sativa subs japonica (var. Nipponbare) physical map.

CONCLUSION

We characterized new resistance QTLs. While some co-localize with known resistance genes/QTLs to Asian strains, others are specific to African strains. We result with new information on genes and QTLs for resistance to bacterial blight that will be useful for controlling the disease.

Genetic Characterization of Green Bean (Phaseolus vulgaris L.) Accessions from Turkey with SCAR and SSR Markers

Characterization, conservation, and utilization of genetic resources is essential for the sustainability in agriculture. Plant genetic resources are important for breeding efforts designed for the generation of new cultivars or for the improvement of existing ones. Green bean has been cultivated extensively in Turkey giving rise to local accessions through selection over time and adaptation to various environmental conditions. The objective of the present study was to determine the genetic relationships of green bean accessions collected from Kırşehir Province of Turkey, located at the central Anatolia. Within a population of 275 green bean accessions, 50 accessions were selected on the basis of morphological observations for further evaluation with SSR and STS/SCAR markers together with 4 reference cultivars of Andean and Mesoamerican origin. SSR markers selected on the basis of high polymorphism information content revealed the genetic relatedness of selected green bean accessions. STS/SCAR markers associated with bean anthracnose, common bacterial blight, white mold, halo blight, and phaseolin protein demonstrated the inheritance of resistance traits of local accessions at the selected loci. These findings may help better utilize genetic resources and furthermore are expected to facilitate forthcoming breeding studies for the generation of novel cultivars well adapted to the region.

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

Background

Common 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.

Results

Transcript-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.

Conclusions

The 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 material

The online version of this article (doi:10.1186/s12864-016-2499-3) contains supplementary material, which is available to authorized users.

Major Contribution of Flowering Time and Vegetative Growth to Plant Production in Common Bean As Deduced from a Comparative Genetic Mapping

Determinacy growth habit and accelerated flowering traits were selected during or after domestication in common bean. Both processes affect several presumed adaptive traits such as the rate of plant production. There is a close association between flowering initiation and vegetative growth; however, interactions among these two crucial developmental processes and their genetic bases remain unexplored. In this study, with the aim to establish the genetic relationships between these complex processes, a multi-environment quantitative trait locus (QTL) mapping approach was performed in two recombinant inbred line populations derived from inter-gene pool crosses between determinate and indeterminate genotypes. Additive and epistatic QTLs were found to regulate flowering time, vegetative growth, and rate of plant production. Moreover, the pleiotropic patterns of the identified QTLs evidenced that regions controlling time to flowering traits, directly or indirectly, are also involved in the regulation of plant production traits. Further QTL analysis highlighted one QTL, on the lower arm of the linkage group Pv01, harboring the Phvul.001G189200 gene, homologous to the Arabidopsis thaliana TERMINAL FLOWER1 (TFL1) gene, which explained up to 32% of phenotypic variation for time to flowering, 66% for vegetative growth, and 19% for rate of plant production. This finding was consistent with previous results, which have also suggested Phvul.001G189200 (PvTFL1y) as a candidate gene for determinacy locus. The information here reported can also be applied in breeding programs seeking to optimize key agronomic traits, such as time to flowering, plant height and an improved reproductive biomass, pods, and seed size, as well as yield.

Genetic Dissection of ICP-Detected Nutrient Accumulation in the Whole Seed of Common Bean (Phaseolus vulgaris L.)

Nutrient transport to grain legume seeds is not well studied and can benefit from modern methods of elemental analysis including spectroscopic techniques. Some cations such as potassium (K) and magnesium (Mg) are needed for plant physiological purposes. Meanwhile, some minerals such as copper (Cu), iron (Fe), molybdenum (Mo), and zinc (Zn) are important micronutrients. Phosphorus (P) is rich in legumes, while sulfur (S) concentration is related to essential amino acids. In this research, the goal was to analyze a genetic mapping population of common bean (Phaseolus vulgaris L.) with inductively coupled plasma (ICP) spectrophotometry to determine concentrations of and to discover quantitative trait loci (QTL) for 15 elements in ground flour of whole seeds. The population was grown in randomized complete block design experiments that had been used before to analyze Fe and Zn. A total of 21 QTL were identified for nine additional elements, of which four QTL were found for Cu followed by three each for Mg, Mn, and P. Fewer QTL were found for K, Na and S. Boron (B) and calcium (Ca) had only one QTL each. The utility of the QTL for breeding adaptation to element deficient soils and association with previously discovered nutritional loci are discussed.

De novo assembly and characterization of foot transcriptome and microsatellite marker development for Paphia textile

Paphia textile is an important, aquaculture bivalve clam species distributed mainly in China, Philippines, and Malaysia. Recent studies of P. textile have focused mainly on artificial breeding and nutrition analysis, and the transcriptome and genome of P. textile have rarely been reported. In this work, the transcriptome of P. textile foot tissue was sequenced on an Illumina HiSeq™ 2000 platform. A total of 20,219,795 reads were generated, resulting in 4.08 Gb of raw data. The raw reads were cleaned and assembled into 54,852 unigenes with an N50 length of 829 bp. Of these unigenes, 38.92% were successfully annotated based on their matches to sequences in seven public databases. Among the annotated unigenes, 14,571 were assigned Gene Ontology terms, 5448 were classified to Clusters of Orthologous Groups categories, and 6738 were mapped to 228 pathways in the Kyoto Encyclopedia of Genes and Genomes database. For functional marker development, 5605 candidate simple sequence repeats were identified in the transcriptome and 80 primer pairs were selected randomly and amplified in a wild population of P. textile. A total of 36 loci that exhibited obvious repeat length polymorphisms were detected. The transcriptomic data and microsatellite markers will provide valuable resources for future functional gene analyses, genetic map construction, and quantitative trait loci mapping in P. textile.

Construction of an SSR and RAD-Marker Based Molecular Linkage Map of Vigna vexillata (L.) A. Rich

Vigna vexillata (L.) A. Rich. (tuber cowpea) is an underutilized crop for consuming its tuber and mature seeds. Wild form of V. vexillata is a pan-tropical perennial herbaceous plant which has been used by local people as a food. Wild V. vexillata has also been considered as useful gene(s) source for V. unguiculata (cowpea), since it was reported to have various resistance gene(s) for insects and diseases of cowpea. To exploit the potential of V. vexillata, an SSR-based linkage map of V. vexillata was developed. A total of 874 SSR markers successfully amplified single DNA fragment in V. vexillata among 1,336 SSR markers developed from Vigna angularis (azuki bean), V. unguiculata and Phaseolus vulgaris (common bean). An F2 population of 300 plants derived from a cross between salt resistant (V1) and susceptible (V5) accessions was used for mapping. A genetic linkage map was constructed using 82 polymorphic SSR markers loci, which could be assigned to 11 linkage groups spanning 511.5 cM in length with a mean distance of 7.2 cM between adjacent markers. To develop higher density molecular linkage map and to confirm SSR markers position in a linkage map, RAD markers were developed and a combined SSR and RAD markers linkage map of V. vexillata was constructed. A total of 559 (84 SSR and 475 RAD) markers loci could be assigned to 11 linkage groups spanning 973.9 cM in length with a mean distance of 1.8 cM between adjacent markers. Linkage and genetic position of all SSR markers in an SSR linkage map were confirmed. When an SSR genetic linkage map of V. vexillata was compared with those of V. radiata and V. unguiculata, it was suggested that the structure of V. vexillata chromosome was considerably differentiated. This map is the first SSR and RAD marker-based V. vexillata linkage map which can be used for the mapping of useful traits.

Mapping-by-sequencing of Ligon-lintless-1 (Li 1 ) reveals a cluster of neighboring genes with correlated expression in developing fibers of Upland cotton (Gossypium hirsutum L.)

Mapping-by-sequencing and SNP marker analysis were used to fine map the Ligon-lintless-1 ( Li 1 ) short fiber mutation in tetraploid cotton to a 255-kb region that contains 16 annotated proteins. The Ligon-lintless-1 (Li 1 ) mutant of cotton (Gossypium hirsutum L.) has been studied as a model for cotton fiber development since its identification in 1929; however, the causative mutation has not been identified yet. Here we report the fine genetic mapping of the mutation to a 255-kb region that contains only 16 annotated genes in the reference Gossypium raimondii genome. We took advantage of the incompletely dominant dwarf vegetative phenotype to identify 100 mutants (Li 1 /Li 1 ) and 100 wild-type (li 1 /li 1 ) homozygotes from a mapping population of 2567 F2 plants, which we bulked and deep sequenced. Since only homozygotes were sequenced, we were able to use a high stringency in SNP calling to rapidly narrow down the region harboring the Li 1 locus, and designed subgenome-specific SNP markers to test the population. We characterized the expression of all sixteen genes in the region by RNA sequencing of elongating fibers and by RT-qPCR at seven time points spanning fiber development. One of the most highly expressed genes found in this interval in wild-type fiber cells is 40-fold under-expressed at the day of anthesis (DOA) in the mutant fiber cells.  This gene is a major facilitator superfamily protein, part of the large family of proteins that includes auxin and sugar transporters. Interestingly, nearly all genes in this region were most highly expressed at DOA and showed a high degree of co-expression. Further characterization is required to determine if transport of hormones or carbohydrates is involved in both the dwarf and lintless phenotypes of Li 1 plants.

SNP Assay Development for Linkage Map Construction, Anchoring Whole-Genome Sequence, and Other Genetic and Genomic Applications in Common Bean

A total of 992,682 single-nucleotide polymorphisms (SNPs) was identified as ideal for Illumina Infinium II BeadChip design after sequencing a diverse set of 17 common bean (Phaseolus vulgaris L) varieties with the aid of next-generation sequencing technology. From these, two BeadChips each with >5000 SNPs were designed. The BARCBean6K_1 BeadChip was selected for the purpose of optimizing polymorphism among market classes and, when possible, SNPs were targeted to sequence scaffolds in the Phaseolus vulgaris 14× genome assembly with sequence lengths >10 kb. The BARCBean6K_2 BeadChip was designed with the objective of anchoring additional scaffolds and to facilitate orientation of large scaffolds. Analysis of 267 F2 plants from a cross of varieties Stampede × Red Hawk with the two BeadChips resulted in linkage maps with a total of 7040 markers including 7015 SNPs. With the linkage map, a total of 432.3 Mb of sequence from 2766 scaffolds was anchored to create the Phaseolus vulgaris v1.0 assembly, which accounted for approximately 89% of the 487 Mb of available sequence scaffolds of the Phaseolus vulgaris v0.9 assembly. A core set of 6000 SNPs (BARCBean6K_3 BeadChip) with high genotyping quality and polymorphism was selected based on the genotyping of 365 dry bean and 134 snap bean accessions with the BARCBean6K_1 and BARCBean6K_2 BeadChips. The BARCBean6K_3 BeadChip is a useful tool for genetics and genomics research and it is widely used by breeders and geneticists in the United States and abroad.

Disease resistance in rice and the role of molecular breeding in protecting rice crops against diseases

Rice diseases (bacterial, fungal, or viral) threaten food productivity. Host resistance is the most efficient, environmentally friendly method to cope with such diverse pathogens. Quantitative resistance conferred by quantitative trait loci (QTLs) is a valuable resource for rice disease resistance improvement. Although QTLs confer partial but durable resistance to many pathogen species in different crop plants, the molecular mechanisms of quantitative disease resistance remain mostly unknown. Quantitative resistance and non-host resistance are types of broad-spectrum resistance, which are mediated by resistance (R) genes. Because R genes activate different resistance pathways, investigating the genetic spectrum of resistance may lead to minimal losses from harmful diseases. Genome studies can reveal interactions between different genes and their pathways and provide insight into gene functions. Protein–protein interaction (proteomics) studies using molecular and bioinformatics tools may further enlighten our understanding of resistance phenomena.

Genetic diversity and population structure of an Italian landrace of runner bean (Phaseolus coccineus L.): inferences for its safeguard and on-farm conservation

The landraces are considered important sources of valuable germplasm for breeding activities to face climatic changes as well as to satisfy the requirement of new varieties for marginal areas. Runner bean (Phaseolus coccineus L.) is one of the most cultivated Phaseolus species worldwide, but few studies have been addressed to assess the genetic diversity and structure within and among landrace populations. In the present study, 20 different populations of a runner bean landrace from Central Italy named "Fagiolone," together with 41 accessions from Italy and Mesoamerica, were evaluated by using 14 nuclear SSRs to establish its genetic structure and distinctiveness. Results indicated that "Fagiolone" landrace can be considered as a dynamic evolving open-pollinated population that shows a significant level of genetic variation, mostly detected within populations, and the presence of two main genetic groups, of which one distinguished from other Italian runner bean landraces. Results highlighted also a relevant importance of farmers' management practices able to influence the genetic structure of this landrace, in particular the seed exchanges and selection, and the past introduction in cultivation of landraces/cultivars similar to seed morphology, but genetically rather far from "Fagiolone." The most suitable on-farm strategies for seed collection, conservation and multiplication will be defined based on our results, as a model for threatened populations of other allogamous crop species. STRUCTURE and phylogenetic analyses indicated that Mesoamerican accessions and Italian landraces belong to two distinct gene pools confirming the hypothesis that Europe could be considered a secondary diversification center for P. coccineus.

Punctuated distribution of recombination hotspots and demarcation of pericentromeric regions in Phaseolus vulgaris L

High density genetic maps are a reliable tool for genetic dissection of complex plant traits. Mapping resolution is often hampered by the variable crossover and non-crossover events occurring across the genome, with pericentromeric regions (pCENR) showing highly suppressed recombination rates. The efficiency of linkage mapping can further be improved by characterizing and understanding the distribution of recombinational activity along individual chromosomes. In order to evaluate the genome wide recombination rate in common beans (Phaseolus vulgaris L.) we developed a SNP-based linkage map using the genotype-by-sequencing approach with a 188 recombinant inbred line family generated from an inter gene pool cross (Andean x Mesoamerican). We identified 1,112 SNPs that were subsequently used to construct a robust linkage map with 11 groups, comprising 513 recombinationally unique marker loci spanning 943 cM (LOD 3.0). Comparative analysis showed that the linkage map spanned >95% of the physical map, indicating that the map is almost saturated. Evaluation of genome-wide recombination rate indicated that at least 45% of the genome is highly recombinationally suppressed, and allowed us to estimate locations of pCENRs. We observed an average recombination rate of 0.25 cM/Mb in pCENRs as compared to the rest of genome that showed 3.72 cM/Mb. However, several hot spots of recombination were also detected with recombination rates reaching as high as 34 cM/Mb. Hotspots were mostly found towards the end of chromosomes, which also happened to be gene-rich regions. Analyzing relationships between linkage and physical map indicated a punctuated distribution of recombinational hot spots across the genome.

Association mapping for five agronomic traits in the common bean (Phaseolus vulgaris L.)

BACKGROUND

The common bean is the most important grain legume and a major source of protein in many developing countries. We analysed the following traits: pod fibre (PF), seeds per pod (SPP), plant type (PT), growth habit (GH), and days to flowering (DF) for a set of diverse common bean accessions and determined whether such traits were associated with amplified fragment length polymorphism (AFLP), simple sequence repeat (SSR) and single nucleotide polymorphism (SNP) markers.

RESULTS

In this study, 66 common bean genotypes were used and genotyped with 233 AFLP, 105 SNP and 80 SSR markers. The association analysis between markers and five traits was performed using a General Linear Model (GLM) in Trait Analysis by aSSociation, Evolution and Linkage (TASSEL). The population structure was determined using the STRUCTURE software, and seven groups (K = 7) were identified among genotypes. The associations for such traits were identified and quantified; 62 markers were associated with the five traits.

CONCLUSION

This study demonstrated that association mapping using a reasonable number of markers, distributed across the genome and with the appropriate number of individuals harboured to detect DNA markers linked to the traits of PF, SPP, PT, GH and DF in common bean.

Unraveling the efficiency of RAPD and SSR markers in diversity analysis and population structure estimation in common bean

Increase in food production viz-a-viz quality of food is important to feed the growing human population to attain food as well as nutritional security. The availability of diverse germplasm of any crop is an important genetic resource to mine the genes that may assist in attaining food as well as nutritional security. Here we used 15 RAPD and 23 SSR markers to elucidate diversity among 51 common bean genotypes mostly landraces collected from the Himalayan region of Jammu and Kashmir, India. We observed that both the markers are highly polymorphic. The discriminatory power of these markers was determined using various parameters like; percent polymorphism, PIC, resolving power and marker index. 15 RAPDs produced 171 polymorphic bands, while 23 SSRs produced 268 polymorphic bands. SSRs showed a higher PIC value (0.300) compared to RAPDs (0.243). Further the resolving power of SSRs was 5.241 compared to 3.86 for RAPDs. However, RAPDs showed a higher marker index (2.69) compared to SSRs (1.279) that may be attributed to their higher multiplex ratio. The dendrograms generated with hierarchical UPGMA cluster analysis grouped genotypes into two main clusters with various degrees of sub clustering within the cluster. Here we observed that both the marker systems showed comparable accuracy in grouping genotypes of common bean according to their area of cultivation. The model based STRUCTURE analysis using 15 RAPD and 23 SSR markers identified a population with 3 sub-populations which corresponds to distance based groupings. High level of genetic diversity was observed within the population. These findings have further implications in common bean breeding as well as conservation programs.

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.

BAC-end microsatellites from intra and inter-genic regions of the common bean genome and their correlation with cytogenetic features

Highly polymorphic markers such as simple sequence repeats (SSRs) or microsatellites are very useful for genetic mapping. In this study novel SSRs were identified in BAC-end sequences (BES) from non-contigged, non-overlapping bacterial artificial clones (BACs) in common bean (Phaseolus vulgaris L.). These so called "singleton" BACs were from the G19833 Andean gene pool physical map and the new BES-SSR markers were used for the saturation of the inter-gene pool, DOR364×G19833 genetic map. A total of 899 SSR loci were found among the singleton BES, but only 346 loci corresponded to the single di- or tri-nucleotide motifs that were likely to be polymorphic (ATT or AG motifs, principally) and useful for primer design and individual marker mapping. When these novel SSR markers were evaluated in the DOR364×G19833 population parents, 136 markers revealed polymorphism and 106 were mapped. Genetic mapping resulted in a map length of 2291 cM with an average distance between markers of 5.2 cM. The new genetic map was compared to the most recent cytogenetic analysis of common bean chromosomes. We found that the new singleton BES-SSR were helpful in filling peri-centromeric spaces on the cytogenetic map. Short genetic distances between some new singleton-derived BES-SSR markers was common showing suppressed recombination in these regions compared to other parts of the genome. The correlation of singleton-derived SSR marker distribution with other cytogenetic features of the bean genome is discussed.

Fine mapping of a dominant thermo-sensitive genic male sterility gene (BntsMs) in rapeseed (Brassica napus) with AFLP- and Brassica rapa-derived PCR markers

A new thermo-sensitive dominant genic male sterility (TSDGMS) line of Brassica napus was found and mapped in this paper. Our result will greatly accelerate the map-based cloning of the BntsMs gene. TE5A is a thermo-sensitive dominant genic male sterility line originating from spontaneous mutation of the inbred line TE5 in Brassica napus and provides a promising system for the development of hybrid cultivars. Genetic analysis has revealed that the BntsMs mutant is controlled by a single, dominant gene. Here, we describe the fine mapping of BntsMs using amplified fragment length polymorphism (AFLP) and intron polymorphism (IP) methodologies. We screened 1,024 primer combinations and then identified five AFLP markers linked to the BntsMs gene, two of which were successfully converted into sequence-characterised amplified region (SCAR) markers. The linkage of the markers was identified by analysing a large BC2 population of 700 recessive-fertility individuals. Two SCAR markers were found in the flanking region of the BntsMs gene at distance of 3.5 and 4.8 cm. Based on sequence information from the previously screened AFLP markers and on genome organisation comparisons of the A genome of Brassica rapa and Arabidopsis, seven IP markers linked to the BntsMs gene were developed. By analysing the 700 recessive-fertility individuals, two IP markers, IP004 and IP470, were localised to the flanking region of the BntsMs gene at a distance of 0.3 and 0.2 cm, respectively. A comparison of the B. rapa and Arabidopsis genomes revealed 27 genes of B. rapa in the flanking region of these two IP markers. It is likely that the molecular markers developed from these investigations will greatly accelerate the positional cloning of the BntsMs gene.

Genetic assessment of common bean (Phaseolus vulgaris L.) accessions by peroxidase gene-based markers

BACKGROUND

Peroxidase, a plant-specific oxidoreductase, is a heme-containing glycoprotein encoded by a large multigenic family in plants. Plant peroxidases (POXs, EC 1.11.1.7) play important roles in many self-defense interactions in plants. Here, 67 common bean (Phaseolus vulgaris L.) genotypes were studied using a POX gene-based marker method. Comparison of POX genes could resolve evolutionary relationships in common bean.

RESULTS

Eighty fragments were obtained with 20 primer pairs that amplified one (POX8c) to eight (ATP29) bands, with a mean of four bands per primer pair. The average (polymorphic information content) PIC value for the POX products was 0.40. The maximum variation (93%) was found between Turkey (#33) and India (#52) and between Antalya (#33) and India (#53). The minimum variation (0%) was found among four pairs: Bozdag (#2) and Karadeniz (#38), Kirklareli (#11) and Turkey (#15, 16, 43), Bandirma (#13) and Turkey (#15, 16, 43), and Kirklareli (#10) and Bandirma (#22). UPGMA was used to discriminate the common bean genotypes into five clusters, while STRUCTURE software was used to investigate the genetic population structure.

CONCLUSION

The results showed that POX gene family markers can be used to study genotypic diversity and provide new information for breeding programs and common bean improvement practices.

Developing market class specific InDel markers from next generation sequence data in Phaseolus vulgaris L

Next generation sequence data provides valuable information and tools for genetic and genomic research and offers new insights useful for marker development. This data is useful for the design of accurate and user-friendly molecular tools. Common bean (Phaseolus vulgaris L.) is a diverse crop in which separate domestication events happened in each gene pool followed by race and market class diversification that has resulted in different morphological characteristics in each commercial market class. This has led to essentially independent breeding programs within each market class which in turn has resulted in limited within market class sequence variation. Sequence data from selected genotypes of five bean market classes (pinto, black, navy, and light and dark red kidney) were used to develop InDel-based markers specific to each market class. Design of the InDel markers was conducted through a combination of assembly, alignment and primer design software using 1.6× to 5.1× coverage of Illumina GAII sequence data for each of the selected genotypes. The procedure we developed for primer design is fast, accurate, less error prone, and higher throughput than when they are designed manually. All InDel markers are easy to run and score with no need for PCR optimization. A total of 2687 InDel markers distributed across the genome were developed. To highlight their usefulness, they were employed to construct a phylogenetic tree and a genetic map, showing that InDel markers are reliable, simple, and accurate.

High-resolution genetic mapping of rice bacterial blight resistance gene Xa23

Bacterial blight (BB) caused by Xanthomonas oryzae pv. oryzae (Xoo) is the most devastating bacterial disease of rice (Oryza sativa L.), a staple food crop that feeds half of the world's population. In management of this disease, the most economical and effective approach is cultivating resistant varieties. Due to rapid change of pathogenicity in the pathogen, it is necessary to identify and characterize more host resistance genes for breeding new resistant varieties. We have previously identified the BB resistance (R) gene Xa23 that confers the broadest resistance to Xoo strains isolated from different rice-growing regions and preliminarily mapped the gene within a 1.7 cm region on the long arm of rice chromosome 11. Here, we report fine genetic mapping and in silico analysis of putative candidate genes of Xa23. Based on F2 mapping populations derived from crosses between Xa23-containing rice line CBB23 and susceptible varieties JG30 or IR24, six new STS markers Lj36, Lj46, Lj138, Lj74, A83B4, and Lj13 were developed. Linkage analysis revealed that the new markers were co-segregated with or closely linked to the Xa23 locus. Consequently, the Xa23 gene was mapped within a 0.4 cm region between markers Lj138 and A83B4, in which the co-segregating marker Lj74 was identified. The corresponding physical distance between Lj138 and A83B4 on Nipponbare genome is 49.8 kb. Six Xa23 candidate genes have been annotated, including four candidate genes encoding hypothetical proteins and the other two encoding a putative ADP-ribosylation factor protein and a putative PPR protein. These results will facilitate marker-assisted selection of Xa23 in rice breeding and molecular cloning of this valuable R gene.