Linkage disequilibrium and population structure in wild and domesticated populations of Phaseolus vulgaris L

Evidence for introduction bottleneck and extensive inter-gene pool (Mesoamerica x Andes) hybridization in the European common bean (Phaseolus vulgaris L.) germplasm

Common bean diversity within and between Mesoamerican and Andean gene pools was compared in 89 landraces from America and 256 landraces from Europe, to elucidate the effects of bottleneck of introduction and selection for adaptation during the expansion of common bean (Phaseolus vulgaris L.) in Europe. Thirteen highly polymorphic nuclear microsatellite markers (nuSSRs) were used to complement chloroplast microsatellite (cpSSRs) and nuclear markers (phaseolin and Pv-shatterproof1) data from previous studies. To verify the extent of the introduction bottleneck, inter-gene pool hybrids were distinguished from "pure" accessions. Hybrids were identified on the basis of recombination of gene pool specific cpSSR, phaseolin and Pv-shatterproof1 markers with a Bayesian assignments based on nuSSRs, and with STRUCTURE admixture analysis. More hybrids were detected than previously, and their frequency was almost four times larger in Europe (40.2%) than in America (12.3%). The genetic bottleneck following the introduction into Europe was not evidenced in the analysis including all the accessions, but it was significant when estimated only with "pure" accessions, and five times larger for Mesoamerican than for Andean germplasm. The extensive inter-gene pool hybridization generated a large amount of genotypic diversity that mitigated the effects of the bottleneck that occurred when common bean was introduced in Europe. The implication for evolution and the advantages for common bean breeding are discussed.

Genetic Diversity and Population Structure of Tetraploid Wheats (Triticum turgidum L.) Estimated by SSR, DArT and Pedigree Data

Levels of genetic diversity and population genetic structure of a collection of 230 accessions of seven tetraploid Triticum turgidum L. subspecies were investigated using six morphological, nine seed storage protein loci, 26 SSRs and 970 DArT markers. The genetic diversity of the morphological traits and seed storage proteins was always lower in the durum wheat compared to the wild and domesticated emmer. Using Bayesian clustering (K = 2), both of the sets of molecular markers distinguished the durum wheat cultivars from the other tetraploid subspecies, and two distinct subgroups were detected within the durum wheat subspecies, which is in agreement with their origin and year of release. The genetic diversity of morphological traits and seed storage proteins was always lower in the improved durum cultivars registered after 1990, than in the intermediate and older ones. This marked effect on diversity was not observed for molecular markers, where there was only a weak reduction. At K >2, the SSR markers showed a greater degree of resolution than for DArT, with their identification of a greater number of groups within each subspecies. Analysis of DArT marker differentiation between the wheat subspecies indicated outlier loci that are potentially linked to genes controlling some important agronomic traits. Among the 211 loci identified under selection, 109 markers were recently mapped, and some of these markers were clustered into specific regions on chromosome arms 2BL, 3BS and 4AL, where several genes/quantitative trait loci (QTLs) are involved in the domestication of tetraploid wheats, such as the tenacious glumes (Tg) and brittle rachis (Br) characteristics. On the basis of these results, it can be assumed that the population structure of the tetraploid wheat collection partially reflects the evolutionary history of Triticum turgidum L. subspecies and the genetic potential of landraces and wild accessions for the detection of unexplored alleles.

Inconsistent genetic structure among members of a multitrophic system: did bruchid parasitoids (Horismenus spp.) escape the effects of bean domestication?

Anthropogenic range expansion and cultural practices have modified the distribution, abundance and genetic diversity of domesticated organisms, thereby altering multitrophic assemblages through space and time. The putative Mesoamerican domestication centre of the common bean, Phaseolus vulgaris L., in Mexico allows investigating the effects of plant domestication on the genetic structure of members of a multitrophic system. The aim of this study was to compare the evolutionary history of Horismenus parasitoids (Hymenoptera: Eulophidae) to those of their bruchid beetle hosts (Coleoptera: Bruchidae) and their domesticated host plant (P. vulgaris), in the context of traditional agriculture in Mexico. We analyzed the population genetic structure of four Horismenus species in Mexico using mitochondrial COI haplotype data. The two most abundant parasitoid species were Horismenus depressus and Horismenus missouriensis. Horismenus missouriensis were infected by Wolbachia endosymbionts and had little to no population differentiation (F(ST) = 0.06). We suspect the mitochondrial history of H. missouriensis to be blurred by Wolbachia, because differentiation among infected vs. non-infected individuals exists (F(ST) = 0.11). Populations of H. depressus were found to be highly differentiated (F(ST) = 0.34), but the genetic structuring could not be explained by tested spatial components. We then compared the genetic structure observed in this parasitoid species to previously published studies on bruchid beetles and their host plants. Despite extensive human-mediated migration and likely population homogenization of its two Acanthoscelides bruchid beetle hosts, H. depressus populations are structured like its host plant, by a recent dispersal from a diverse ancestral gene pool. Distinct evolutionary dynamics may explain inconsistent patterns among trophic levels. Parasitoids likely migrate from wild bean populations and are poorly adapted to bean storage conditions similar to their bruchid beetle hosts. Integrating several trophic levels to the study of evolutionary history has proven to be fruitful in detecting different ecological responses to human-mediated disturbances and host parasite interactions.

Demographic factors shaped diversity in the two gene pools of wild common bean Phaseolus vulgaris L

Wild common bean (Phaseolus vulgaris L.) is distributed throughout the Americas from Mexico to northern Argentina. Within this range, the species is divided into two gene pools (Andean and Middle American) along a latitudinal gradient. The diversity of 24 wild common bean genotypes from throughout the geographic range of the species was described by using sequence data from 13 loci. An isolation-migration model was evaluated using a coalescent analysis to estimate multiple demographic parameters. Using a Bayesian approach, Andean and Middle American subpopulations with high percentage of parentages were observed. Over all loci, the Middle American gene pool was more diverse than the Andean gene pool (π(sil)=0.0089 vs 0.0068). The two subpopulations were strongly genetically differentiated over all loci (F(st)=0.29). It is estimated that the two current wild gene pools diverged from a common ancestor ∼111 000 years ago. Subsequently, each gene pool underwent a bottleneck immediately after divergence and lasted ∼40 000 years. The Middle American bottleneck population size was ∼46% of the ancestral population size, whereas the Andean was 26%. Continuous asymmetric gene flow was detected between the two gene pools with a larger number of migrants entering Middle American gene pool from the Andean gene pool. These results suggest that because of the complex population structure associated with the ancestral divergence, subsequent bottlenecks in each gene pool, gene pool-specific domestication and intense selection within each gene pool by breeders; association mapping would best be practised within each common bean gene pool.

European Phaseolus coccineus L. landraces: population structure and adaptation, as revealed by cpSSRs and phenotypic analyses

Relatively few studies have extensively analysed the genetic diversity of the runner bean through molecular markers. Here, we used six chloroplast microsatellites (cpSSRs) to investigate the cytoplasmic diversity of 331 European domesticated accessions of the scarlet runner bean (Phaseolus coccineus L.), including the botanical varieties albiflorus, bicolor and coccineus, and a sample of 49 domesticated and wild accessions from Mesoamerica. We further explored the pattern of diversity of the European landraces using 12 phenotypic traits on 262 individuals. For 158 European accessions, we studied the relationships between cpSSR polymorphisms and phenotypic traits. Additionally, to gain insights into the role of gene flow and migration, for a subset of 115 accessions, we compared and contrasted the results obtained by cpSSRs and phenotypic traits with those obtained in a previous study with 12 nuclear microsatellites (nuSSRs). Our results suggest that both demographic and selective factors have roles in the shaping of the population genetic structure of the European runner bean. In particular, we infer the existence of a moderate-to-strong cytoplasmic bottleneck that followed the expansion of the crop into Europe, and we deduce multiple domestication events for this species. We also observe an adaptive population differentiation in the phenology across a latitudinal gradient, which suggests that selection led to the diversification of the runner bean in Europe. The botanical varieties albiflorus, bicolor and coccineus, which are based solely on flower colour, cannot be distinguished based on these cpSSRs and nuSSRs, nor according to the 12 quantitative traits.

Genetic Diversity and Population Structure of Tetraploid Wheats (Triticum turgidum L.) Estimated by SSR, DArT and Pedigree Data

Levels of genetic diversity and population genetic structure of a collection of 230 accessions of seven tetraploid Triticum turgidum L. subspecies were investigated using six morphological, nine seed storage protein loci, 26 SSRs and 970 DArT markers. The genetic diversity of the morphological traits and seed storage proteins was always lower in the durum wheat compared to the wild and domesticated emmer. Using Bayesian clustering (K = 2), both of the sets of molecular markers distinguished the durum wheat cultivars from the other tetraploid subspecies, and two distinct subgroups were detected within the durum wheat subspecies, which is in agreement with their origin and year of release. The genetic diversity of morphological traits and seed storage proteins was always lower in the improved durum cultivars registered after 1990, than in the intermediate and older ones. This marked effect on diversity was not observed for molecular markers, where there was only a weak reduction. At K >2, the SSR markers showed a greater degree of resolution than for DArT, with their identification of a greater number of groups within each subspecies. Analysis of DArT marker differentiation between the wheat subspecies indicated outlier loci that are potentially linked to genes controlling some important agronomic traits. Among the 211 loci identified under selection, 109 markers were recently mapped, and some of these markers were clustered into specific regions on chromosome arms 2BL, 3BS and 4AL, where several genes/quantitative trait loci (QTLs) are involved in the domestication of tetraploid wheats, such as the tenacious glumes (Tg) and brittle rachis (Br) characteristics. On the basis of these results, it can be assumed that the population structure of the tetraploid wheat collection partially reflects the evolutionary history of Triticum turgidum L. subspecies and the genetic potential of landraces and wild accessions for the detection of unexplored alleles.

Molecular analysis of the parallel domestication of the common bean (Phaseolus vulgaris) in Mesoamerica and the Andes

We have studied the nucleotide diversity of common bean, Phaseolus vulgaris, which is characterized by two independent domestications in two geographically distinct areas: Mesoamerica and the Andes. This provides an important model, as domestication can be studied as a replicate experiment. We used nucleotide data from five gene fragments characterized by large introns to analyse 214 accessions (102 wild and 112 domesticated). The wild accessions represent a cross-section of the entire geographical distribution of P. vulgaris. A reduction in genetic diversity in both of these gene pools was found, which was three-fold greater in Mesoamerica compared with the Andes. This appears to be a result of a bottleneck that occurred before domestication in the Andes, which strongly impoverished this wild germplasm, leading to the minor effect of the subsequent domestication bottleneck (i.e. sequential bottleneck). These findings show the importance of considering the evolutionary history of crop species as a major factor that influences their current level and structure of genetic diversity. Furthermore, these data highlight a single domestication event within each gene pool. Although the findings should be interpreted with caution, this evidence indicates the Oaxaca valley in Mesoamerica, and southern Bolivia and northern Argentina in South America, as the origins of common bean domestication.

Diversification and population structure in common beans (Phaseolus vulgaris L.)

Wild accessions of crops and landraces are valuable genetic resources for plant breeding and for conserving alleles and gene combinations in planta. The primary genepool of cultivated common beans includes wild accessions of Phaseolus vulgaris. These are of the same species as the domesticates and therefore are easily crossable with cultivated accessions. Molecular marker assessment of wild beans and landraces is important for the proper utilization and conservation of these important genetic resources. The goal of this research was to evaluate a collection of wild beans with fluorescent microsatellite or simple sequence repeat markers and to determine the population structure in combination with cultivated beans of all known races. Marker diversity in terms of average number of alleles per marker was high (13) for the combination of 36 markers and 104 wild genotypes that was similar to the average of 14 alleles per marker found for the 606 cultivated genotypes. Diversity in wild beans appears to be somewhat higher than in cultivated beans on a per genotype basis. Five populations or genepools were identified in structure analysis of the wild beans corresponding to segments of the geographical range, including Mesoamerican (Mexican), Guatemalan, Colombian, Ecuadorian-northern Peruvian and Andean (Argentina, Bolivia and Southern Peru). The combined analysis of wild and cultivated accessions showed that the first and last of these genepools were related to the cultivated genepools of the same names and the penultimate was found to be distinct but not ancestral to the others. The Guatemalan genepool was very novel and perhaps related to cultivars of race Guatemala, while the Colombian population was also distinct. Results suggest geographic isolation, founder effects or natural selection could have created the different semi-discrete populations of wild beans and that multiple domestications and introgression were involved in creating the diversity of cultivated beans.

Nucleotide diversity patterns at the drought-related DREB2 encoding genes in wild and cultivated common bean (Phaseolus vulgaris L.)

Common beans are an important food legume faced with a series of abiotic stresses the most severe of which is drought. The crop is interesting as a model for the analysis of gene phylogenies due to its domestication process, race structure, and origins in a group of wild common beans found along the South American Andes and the region of Mesoamerica. Meanwhile, the DREB2 transcription factors have been implicated in controlling non-ABA dependent responses to drought stress. With this in mind our objective was to study in depth the genetic diversity for two DREB2 genes as possible candidates for association with drought tolerance through a gene phylogenetic analysis. In this genetic diversity assessment, we analyzed nucleotide diversity at the two candidate genes Dreb2A and Dreb2B, in partial core collections of 104 wild and 297 cultivated common beans with a total of 401 common bean genotypes from world-wide germplasm analyzed. Our wild population sample covered a range of semi-mesic to very dry habitats, while our cultivated samples presented a wide spectrum of low to high drought tolerance. Both genes showed very different patterns of nucleotide variation. Dreb2B exhibited very low nucleotide diversity relative to neutral reference loci previously surveyed in these populations. This suggests that strong purifying selection has been acting on this gene. In contrast, Dreb2A exhibited higher levels of nucleotide diversity, which is indicative of adaptive selection and population expansion. These patterns were more distinct in wild compared to cultivated common beans. These approximations suggested the importance of Dreb2 genes in the context of drought tolerance, and constitute the first steps towards an association study between genetic polymorphism of this gene family and variation in drought tolerance traits. We discuss the utility of allele mining in the DREB gene family for the discovery of new drought tolerance traits from wild common bean.

Population structure and linkage disequilibrium in Lupinus albus L. germplasm and its implication for association mapping

White lupin (Lupinus albus L.) has been around since 300 B.C. and is recognized for its ability to grow on poor soils and application as green manure in addition to seed harvest. The seed has very high levels of protein (33-47 %) and oil (6-13 %). It also has many secondary metabolites that are potentially of nutraceutical value to animals and humans. Despite such a great potential, lupins role in modern agriculture began only in the twentieth century. Although a large collection of Lupinus germplasm accessions is available worldwide, rarely have they been genetically characterized. Additionally, scarce genomic resources in terms of recombinant populations and genome information have been generated for L. albus. With the advancement in association mapping methods, the natural populations have the potential to replace the recombinant populations in gene mapping and marker-trait associations. Therefore, we studied the genetic similarity, population structure and marker-trait association in a USDA germplasm collection for their current and future application in this crop improvement. A total of 122 PI (Plant Inventory) lines were screened with 18 AFLP primer pairs that generated 2,277 fragments. A subset of 892 polymorphic markers with MAF >0.05 (minor allele frequency) were used for association mapping. The cluster analysis failed to group accessions on the basis of their passport information, and a weak structure and low linkage disequilibrium (LD) were observed indicating the usefulness of the collection for association mapping. Moreover, we were also able to identify two markers (a p value of 1.53 × 10(-4) and 2.3 × 10(-4)) that explained 22.69 and 20.5 % of seed weight variation determined using R (LR) (2) . The implications of lack of geographic clustering, population structure, low LD and the ability of AFLP to map seed weight trait using association mapping and the usefulness of the PI collections in breeding programs are discussed.

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

Background

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

Results

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

Conclusions

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

Genetic structure and linkage disequilibrium in landrace populations of barley in Sardinia

Multilocus digenic linkage disequilibria (LD) and their population structure were investigated in eleven landrace populations of barley (Hordeum vulgare ssp. vulgare L.) in Sardinia, using 134 dominant simple-sequence amplified polymorphism markers. The analysis of molecular variance for these markers indicated that the populations were partially differentiated (F(ST) = 0.18), and clustered into three geographic areas. Consistent with this population pattern, STRUCTURE analysis allocated individuals from a bulk of all populations into four genetic groups, and these groups also showed geographic patterns. In agreement with other molecular studies in barley, the general level of LD was low (13% of locus pairs, with P < 0.01) in the bulk of 337 lines, and decayed steeply with map distance between markers. The partitioning of multilocus associations into various components indicated that genetic drift and founder effects played a major role in determining the overall genetic makeup of the diversity in these landrace populations, but that epistatic homogenising or diversifying selection was also present. Notably, the variance of the disequilibrium component was relatively high, which implies caution in the pooling of barley lines for association studies. Finally, we compared the analyses of multilocus structure in barley landrace populations with parallel analyses in both composite crosses of barley on the one hand and in natural populations of wild barley on the other. Neither of these serves as suitable mimics of landraces in barley, which require their own study. Overall, the results suggest that these populations can be exploited for LD mapping if population structure is controlled.

Mesoamerican origin of the common bean (Phaseolus vulgaris L.) is revealed by sequence data

Knowledge about the origins and evolution of crop species represents an important prerequisite for efficient conservation and use of existing plant materials. This study was designed to solve the ongoing debate on the origins of the common bean by investigating the nucleotide diversity at five gene loci of a large sample that represents the entire geographical distribution of the wild forms of this species. Our data clearly indicate a Mesoamerican origin of the common bean. They also strongly support the occurrence of a bottleneck during the formation of the Andean gene pool that predates the domestication, which was suggested by recent studies based on multilocus molecular markers. Furthermore, a remarkable result was the genetic structure that was seen for the Mesoamerican accessions, with the identification of four different genetic groups that have different relationships with the sets of wild accessions from the Andes and northern Peru-Ecuador. This finding implies that both of the gene pools from South America originated through different migration events from the Mesoamerican populations that were characteristic of central Mexico.

Chloroplast Microsatellite Diversity in Phaseolus vulgaris

Evolutionary studies that are aimed at defining the processes behind the present level and organization of crop genetic diversity represent the fundamental bases for biodiversity conservation and use. A Mesoamerican origin of the common bean Phaseolus vulgaris was recently suggested through analysis of nucleotide polymorphism at the nuclear level. Here, we have used chloroplast microsatellites to investigate the origin of the common bean, on the basis of the specific characteristics of these markers (no recombination, haploid genome, uniparental inheritance), to validate these recent findings. Indeed, comparisons of the results obtained through analysis of nuclear and cytoplasmic DNA should allow the resolution of some of the contrasting information available on the evolutionary processes. The main outcomes of the present study are: (i) confirmation at the chloroplast level of the results obtained through nuclear data, further supporting the Mesoamerican origin of P. vulgaris, with central Mexico representing the cradle of its diversity; (ii) identification of a putative ancestral plastidial genome, which is characteristic of a group of accessions distributed from central Mexico to Peru, but which have not been highlighted beforehand through analyses at the nuclear level. Finally, the present study suggests that when a single species is analyzed, there is the need to take into account the complexity of the relationships between P. vulgaris and its closely related and partially intercrossable species P. coccineus and P. dumosus. Thus, the present study stresses the importance for the investigation of the speciation processes of these taxa through comparisons of both plastidial and nuclear variability. This knowledge will be fundamental not only from an evolutionary point of view, but also to put P. coccineus and P. dumosus germplasm to better use as a source of useful diversity for P. vulgaris breeding.

Two methods to easily obtain nucleotide sequences from AFLP loci of interest

Genome scans based on anonymous Amplified Fragment Length Polymorphism (AFLP) markers scattered throughout the genome are becoming an increasingly popular approach to study the genetic basis of adaptation and speciation in natural populations. A shortcoming of this approach is that despite its efficiency to detect signatures of selection, it can hardly help pinpoint the specific genomic region(s), gene(s), or mutation(s) targeted by selection. Here, we present two methods to be undertaken after performing an AFLP-based genome scan to easily obtain the sequences of AFLP loci detected as outliers by population genomics approaches. The first one is based on the gel excision of the target AFLP fragment, after simplification of the AFLP fingerprint and separation of the fragments by migration. The second one is a combination of classical AFLP protocol and 454 pyrosequencing.

Investigation of the domestication of common bean (Phaseolus vulgaris) using multilocus sequence data

Multilocus sequence data collected from domesticated and related wild relatives provides a rich source of information on the effect of human selection on the diversity and adaptability of a species to complex environments. To evaluate the domestication history of common bean (Phaseolus vulgaris L.), multilocus sequence data from landraces representing the various races within the Middle American (MA) and Andean gene pools was evaluated. Across 13 loci, nucleotide diversity was similar between landraces and wild germplasm in both gene pools. The diversity data were evaluated using the approximate Bayesian computation approach to test multiple domestication models and estimate population demographic parameters. A model with a single domestication event coupled with bidirectional migration between wild and domesticated genotypes fitted the data better than models consisting of two or three domestication events in each genepool. The effective bottleneck population size was ~50% of the base population in each genepool. The bottleneck began ~8200 and ~8500 years before present and ended at ~6300 and ~7000 years before present in MA and Andean gene pools respectively. Linkage disequilibrium decayed to a greater extent in the MA genepool. Given the (1) geographical adaptation bottleneck in each wild gene pool, (2) a subsequent domestication bottleneck within each gene pool, (3) differentiation into gene-pool specific races and (4) variable extents of linkage disequilibrium, association mapping experiments for common bean would more appropriately be performed within each genepool.

Nucleotide diversity of a genomic sequence similar to SHATTERPROOF (PvSHP1) in domesticated and wild common bean (Phaseolus vulgaris L.)

Evolutionary studies in plant and animal breeding are aimed at understanding the structure and organization of genetic variations of species. We have identified and characterized a genomic sequence in Phaseolus vulgaris of 1,200 bp (PvSHP1) that is homologous to SHATTERPROOF-1 (SHP1), a gene involved in control of fruit shattering in Arabidopsis thaliana. The PvSHP1 fragment was mapped to chromosome Pv06 in P. vulgaris and is linked to the flower and seed color gene V. Amplification of the PvSHP1 sequence from the most agronomically important legume species showed a high degree of interspecies diversity in the introns within the Phaseoleae, while the coding region was conserved across distant taxa. Sequencing of the PvSHP1 sequence in a sample of 91 wild and domesticated genotypes that span the geographic distribution of this species in the centers of origin showed that PvSHP1 is highly polymorphic and, therefore, particularly useful to further investigate the origin and domestication history of P. vulgaris. Our data confirm the gene pool structure seen in P. vulgaris along with independent domestication processes in the Andes and Mesoamerica; they provide additional evidence for a single domestication event in Mesoamerica. Moreover, our results support the Mesoamerican origin of this species. Finally, we have developed three indel-spanning markers that will be very useful for bean germplasm characterization, and particularly to trace the distribution of the domesticated Andean and Mesoamerican gene pools.

Genetic diversity analysis of common beans based on molecular markers

A core collection of the common bean (Phaseolus vulgaris L.), representing genetic diversity in the entire Mexican holding, is kept at the INIFAP (Instituto Nacional de Investigaciones Forestales, Agricolas y Pecuarias, Mexico) Germplasm Bank. After evaluation, the genetic structure of this collection (200 accessions) was compared with that of landraces from the states of Oaxaca, Chiapas and Veracruz (10 genotypes from each), as well as a further 10 cultivars, by means of four amplified fragment length polymorphisms (AFLP) +3/+3 primer combinations and seven simple sequence repeats (SSR) loci, in order to define genetic diversity, variability and mutual relationships. Data underwent cluster (UPGMA) and molecular variance (AMOVA) analyses. AFLP analysis produced 530 bands (88.5% polymorphic) while SSR primers amplified 174 alleles, all polymorphic (8.2 alleles per locus). AFLP indicated that the highest genetic diversity was to be found in ten commercial-seed classes from two major groups of accessions from Central Mexico and Chiapas, which seems to be an important center of diversity in the south. A third group included genotypes from Nueva Granada, Mesoamerica, Jalisco and Durango races. Here, SSR analysis indicated a reduced number of shared haplotypes among accessions, whereas the highest genetic components of AMOVA variation were found within accessions. Genetic diversity observed in the common-bean core collection represents an important sample of the total Phaseolus genetic variability at the main Germplasm Bank of INIFAP. Molecular marker strategies could contribute to a better understanding of the genetic structure of the core collection as well as to its improvement and validation.

The heterogeneous levels of linkage disequilibrium in white spruce genes and comparative analysis with other conifers

In plants, knowledge about linkage disequilibrium (LD) is relevant for the design of efficient single-nucleotide polymorphism arrays in relation to their use in population and association genomics studies. Previous studies of conifer genes have shown LD to decay rapidly within gene limits, but exceptions have been reported. To evaluate the extent of heterogeneity of LD among conifer genes and its potential causes, we examined LD in 105 genes of white spruce (Picea glauca) by sequencing a panel of 48 haploid megagametophytes from natural populations and further compared it with LD in other conifer species. The average pairwise r(2) value was 0.19 (s.d.=0.19), and LD dropped quickly with a half-decay being reached at a distance of 65 nucleotides between sites. However, LD was significantly heterogeneous among genes. A first group of 29 genes had stronger LD (mean r(2)=0.28), and a second group of 38 genes had weaker LD (mean r(2)=0.12). While a strong relationship was found with the recombination rate, there was no obvious relationship between LD and functional classification. The level of nucleotide diversity, which was highly heterogeneous across genes, was also not significantly correlated with LD. A search for selection signatures highlighted significant deviations from the standard neutral model, which could be mostly attributed to recent demographic changes. Little evidence was seen for hitchhiking and clear relationships with LD. When compared among conifer species, on average, levels of LD were similar in genes from white spruce, Norway spruce and Scots pine, whereas loblolly pine and Douglas fir genes exhibited a significantly higher LD.

Biodiversity studies in Phaseolus species by DNA barcoding

The potential of DNA barcoding was tested as a system for studying genetic diversity and genetic traceability in bean germplasm. This technique was applied to several pure lines of Phaseolus vulgaris L. belonging to wild, domesticated, and cultivated common beans, along with some accessions of Phaseolus coccineus L., Phaseolus lunatus L., and Vigna unguiculata (L.) Walp. A multilocus approach was exploited using three chloroplast genic regions (rbcL, trnL, and matK), four intergenic spacers (rpoB-trnC, atpBrbcL, trnT-trnL, and psbA-trnH), and nuclear ITS1 and ITS2 rDNA sequences. Our main goals were to identify the markers and SNPs that show the best discriminant power at the variety level in common bean germplasm, to examine two methods (tree based versus character based) for biodiversity analysis and traceability assays, and to evaluate the overall utility of chloroplast DNA barcodes for reconstructing the origins of modern Italian varieties. Our results indicate that the neighbor-joining method is a powerful approach for comparing genetic diversity within plant species, but it is relatively uninformative for the genetic traceability of plant varieties. In contrast, the character-based method was able to identify several distinct haplotypes over all target regions corresponding to Mesoamerican or Andean accessions; Italian accessions originated from both gene pools. On the whole, our findings raise some concerns about the use of DNA barcoding for intraspecific genetic diversity studies in common beans and highlights its limitations for resolving genetic relationships between landraces and varieties.

Association mapping of common bacterial blight resistance QTL in Ontario bean breeding populations

BACKGROUND

Common bacterial blight (CBB), incited by Xanthomonas axonopodis pv. phaseoli (Xap), is a major yield-limiting factor of common bean (Phaseolus vulgaris L.) production around the world. Host resistance is practically the most effective and environmentally-sound approach to control CBB. Unlike conventional QTL discovery strategies, in which bi-parental populations (F2, RIL, or DH) need to be developed, association mapping-based strategies can use plant breeding populations to synchronize QTL discovery and cultivar development.

RESULTS

A population of 469 dry bean lines of different market classes representing plant materials routinely developed in a bean breeding program were used. Of them, 395 lines were evaluated for CBB resistance at 14 and 21 DAI (Days After Inoculation) in the summer of 2009 in an artificially inoculated CBB nursery in south-western Ontario. All lines were genotyped using 132 SNPs (Single Nucleotide Polymorphisms) evenly distributed across the genome. Of the 132 SNPs, 26 SNPs had more than 20% missing data, 12 SNPs were monomorphic, and 17 SNPs had a MAF (Minor Allelic Frequency) of less than 0.20, therefore only 75 SNPs were used for association study, based on one SNP per locus. The best possible population structure was to assign 36% and 64% of the lines into Andean and Mesoamerican subgroups, respectively. Kinship analysis also revealed complex familial relationships among all lines, which corresponds with the known pedigree history. MLM (Mixed Linear Model) analysis, including population structure and kinship, was used to discover marker-trait associations. Eighteen and 22 markers were significantly associated with CBB rating at 14 and 21 DAI, respectively. Fourteen markers were significant for both dates and the markers UBC420, SU91, g321, g471, and g796 were highly significant (p ≤ 0.001). Furthermore, 12 significant SNP markers were co-localized with or close to the CBB-QTLs identified previously in bi-parental QTL mapping studies.

CONCLUSIONS

This study demonstrated that association mapping using a reasonable number of markers, distributed across the genome and with application of plant materials that are routinely developed in a plant breeding program can detect significant QTLs for traits of interest.

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

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

Microsatellite diversity and genetic structure among common bean (Phaseolus vulgaris L.) landraces in Brazil, a secondary center of diversity

Brazil is the largest producer and consumer of common bean (Phaseolus vulgaris L.), which is the most important source of human dietary protein in that country. This study assessed the genetic diversity and the structure of a sample of 279 geo-referenced common bean landraces from Brazil, using molecular markers. Sixty-seven microsatellite markers spread over the 11 linkage groups of the common bean genome, as well as Phaseolin, PvTFL1y, APA and four SCAR markers were used. As expected, the sample showed lower genetic diversity compared to the diversity in the primary center of diversification. Andean and Mesoamerican gene pools were both present but the latter gene pool was four times more frequent than the former. The two gene pools could be clearly distinguished; limited admixture was observed between these groups. The Mesoamerican group consisted of two sub-populations, with a high level of admixture between them leading to a large proportion of stabilized hybrids not observed in the centers of domestication. Thus, Brazil can be considered a secondary center of diversification of common bean. A high degree of genome-wide multilocus associations even among unlinked loci was observed, confirming the high level of structure in the sample and suggesting that association mapping should be conducted in separate Andean and Mesoamerican Brazilian samples.

Beans in Europe: origin and structure of the European landraces of Phaseolus vulgaris L

This study focuses on the expansion of Phaseolus vulgaris in Europe. The pathways of distribution of beans into and across Europe were very complex, with several introductions from the New World that were combined with direct exchanges between European and other Mediterranean countries. We have analyzed here six chloroplast microsatellite (cpSSR) loci and two unlinked nuclear loci (for phaseolin types and Pv-shatterproof1). We have assessed the genetic structure and level of diversity of a large collection of European landraces of P. vulgaris (307) in comparison to 94 genotypes from the Americas that are representative of the Andean and Mesoamerican gene pools. First, we show that most of the European common bean landraces (67%) are of Andean origin, and that there are no strong differences across European regions for the proportions of the Andean and Mesoamerican gene pools. Moreover, cytoplasmic diversity is evenly distributed across European regions. Secondly, the cytoplasmic bottleneck that was due to the introduction of P. vulgaris into the Old World was very weak or nearly absent. This is in contrast to evidence from nuclear analyses that have suggested a bottleneck of greater intensity. Finally, we estimate that a relatively high proportion of the European bean germplasm (about 44%) was derived from hybridization between the Andean and Mesoamerican gene pools. Moreover, although hybrids are present everywhere in Europe, they show an uneven distribution, with high frequencies in central Europe, and low frequencies in Spain and Italy. On the basis of these data, we suggest that the entire European continent and not only some of the countries therein can be regarded as a secondary diversification center for P. vulgaris. Finally, we outline the relevance of these inter-gene pool hybrids for plant breeding.

Genetic differentiation among Maconellicoccus hirsutus (Hemiptera: Pseudococcidae) populations living on different host plants

The pink hibiscus mealybug Maconellicoccus hirsutus (Green) is a dangerous pest that damages a wide variety of agricultural, horticultural, and forestry crops. Amplified fragment length polymorphism (AFLP) fingerprints were used to characterize the genetic variation of 11 M. hirsutus populations infesting three plant species in Nayarit, Mexico. Analysis was carried out using four primers combinations, producing 590 polymorphic bands. Cluster analysis, as well as bootstrap dendrogram and nonmetric multidimensional scaling analysis, grouped M. hirsutus populations according to their host plant. The estimated F(ST) values indicated a high differentiation in M. hirsutus populations among the three host plant species. These results were also supported by a Bayesian analysis, which indicated a population clustering robustness according to their host plant. Genetic variation among populations is not caused by geographic distances, as shown by a Mantel test.

Integration of genome and phenotypic scanning gives evidence of genetic structure in Mesoamerican common bean (Phaseolus vulgaris L.) landraces from the southwest of Europe

Southwestern Europe has been considered as a secondary centre of genetic diversity for the common bean. The dispersal of domesticated materials from their centres of origin provides an experimental system that reveals how human selection during cultivation and adaptation to novel environments affects the genetic composition. In this paper, our goal was to elucidate how distinct events could modify the structure and level of genetic diversity in the common bean. The genome-wide genetic composition was analysed at 42 microsatellite loci in individuals of 22 landraces of domesticated common bean from the Mesoamerican gene pool. The accessions were also characterised for phaseolin seed protein and for nine allozyme polymorphisms and phenotypic traits. One of this study's important findings was the complementary information obtained from all the polymorphisms examined. Most of the markers found to be potentially under the influence of selection were located in the proximity of previously mapped genes and quantitative trait loci (QTLs) related to important agronomic traits, which indicates that population genomics approaches are very efficient in detecting QTLs. As it was revealed by outlier simple sequence repeats, loci analysis with STRUCTURE software and multivariate analysis of phenotypic data, the landraces were grouped into three clusters according to seed size and shape, vegetative growth habit and genetic resistance. A total of 151 alleles were detected with an average of 4 alleles per locus and an average polymorphism information content of 0.31. Using a model-based approach, on the basis of neutral markers implemented in the software STRUCTURE, three clusters were inferred, which were in good agreement with multivariate analysis. Geographic and genetic distances were congruent with the exception of a few putative hybrids identified in this study, suggesting a predominant effect of isolation by distance. Genomic scans using both markers linked to genes affected by selection (outlier) and neutral markers showed advantages relative to other approaches, since they help to create a more complete picture of how adaptation to environmental conditions has sculpted the common bean genomes in southern Europe. The use of outlier loci also gives a clue about what selective forces gave rise to the actual phenotypes of the analysed landraces.

Linkage disequilibrium at the APA insecticidal seed protein locus of common bean (Phaseolus vulgaris L.)

BACKGROUND

An interesting seed protein family with a role in preventing insect herbivory is the multi-gene, APA family encoding the alpha-amylase inhibitor, phytohemagglutinin and arcelin proteins of common bean (Phaseolus vulgaris). Variability for this gene family exists and has been exploited to breed for insect resistance. For example, the arcelin locus has been successfully transferred from wild to cultivated common bean genotypes to provide resistance against the bruchid species Zabrotes subfasciatus although the process has been hampered by a lack of genetic tools for and understanding about the locus. In this study, we analyzed linkage disequilibrium (LD) between microsatellite markers at the APA locus and bruchid resistance in a germplasm survey of 105 resistant and susceptible genotypes and compared this with LD in other parts of the genome.

RESULTS

Microsatellite allele diversity was found to vary with each of the eight APA-linked markers analyzed, and two markers within the APA locus were found to be diagnostic for bruchid resistance or susceptibility and for the different arcelin alleles inherited from the wild accessions. Arc1 was found to provide higher levels of resistance than Arc5 and the markers in the APA locus were highly associated with resistance showing that introgression of this gene-family from wild beans provides resistance in cultivated beans. LD around the APA locus was found to be intermediate compared to other regions of the genome and the highest LD was found within the APA locus itself for example between the markers PV-atct001 and PV-ag004.

CONCLUSIONS

We found the APA locus to be an important genetic determinant of bruchid resistance and also found that LD existed mostly within the APA locus but not beyond it. Moderate LD was also found for some other regions of the genome perhaps related to domestication genes. The LD pattern may reflect the introgression of arcelin from the wild into the cultivated background through breeding. LD and association studies for the arcelin gene, linked genes and other members of the APA family are essential for breaking linkage drag while maintaining high levels of bruchid resistance in common bean.