Patterns of genetic diversity in the Andean gene pool of common bean reveal a candidate domestication gene

Most studies on the genetic diversity of common bean (Phaseolus vulgaris L.) have focussed on accessions from the Mesoamerican gene pool compared to the Andean gene pool. A deeper knowledge of the genetic structure of Argentinian germplasm would enable researchers to determine how the Andean domestication event affected patterns of genetic diversity in domesticated beans and to identify candidates for genes targeted by selection during the evolution of the cultivated common bean. A collection of 116 wild and domesticated accessions representing the diversity of the Andean bean in Argentina was genotyped by means of 114 simple sequence repeat (SSR) markers. Forty-seven Mesoamerican bean accessions and 16 Andean bean accessions representing the diversity of Andean landraces and wild accessions were also included. Using the Bayesian algorithm implemented in the software STRUCTURE we identified five major groups that correspond to Mesoamerican and Argentinian wild accessions and landraces and a group that corresponds to accessions from different Andean and Mesoamerican countries. The neighbour-joining algorithm and principal coordinate clustering analysis confirmed the genetic relationships among accessions observed with the STRUCTURE analysis. Argentinian accessions showed a substantial genetic variation with a considerable number of unique haplotypes and private alleles, suggesting that they may have played an important role in the evolution of the species. The results of statistical analyses aimed at identifying genomic regions with consistent patterns of variation were significant for 35 loci (~20 % of the SSRs used in the Argentinian accessions). One of these loci mapped in or near the genomic region of the glutamate decarboxylase gene. Our data characterize the population structure of the Argentinian germplasm. This information on its diversity will be very valuable for use in introgressing Argentinian genes into commercial varieties because the majority of present-day common bean varieties are of Andean origin.     

Structure of genetic diversity in the two major gene pools of common bean (Phaseolus vulgaris L., Fabaceae)

Domesticated materials with well-known wild relatives provide an experimental system to reveal how human selection during cultivation affects genetic composition and adaptation to novel environments. In this paper, our goal was to elucidate how two geographically distinct domestication events modified the structure and level of genetic diversity in common bean. Specifically, we analyzed the genome-wide genetic composition at 26, mostly unlinked microsatellite loci in 349 accessions of wild and domesticated common bean from the Andean and Mesoamerican gene pools. Using a model-based approach, implemented in the software STRUCTURE, we identified nine wild or domesticated populations in common bean, including four of Andean and four of Mesoamerican origins. The ninth population was the putative wild ancestor of the species, which was classified as a Mesoamerican population. A neighbor-joining analysis and a principal coordinate analysis confirmed genetic relationships among accessions and populations observed with the STRUCTURE analysis. Geographic and genetic distances in wild populations were congruent with the exception of a few putative hybrids identified in this study, suggesting a predominant effect of isolation by distance. Domesticated common bean populations possessed lower genetic diversity, higher F _ST, and generally higher linkage disequilibrium (LD) than wild populations in both gene pools; their geographic distributions were less correlated with genetic distance, probably reflecting seed-based gene flow after domestication. The LD was reduced when analyzed in separate Andean and Mesoamerican germplasm samples. The Andean domesticated race Nueva Granada had the highest F _ST value and widest geographic distribution compared to other domesticated races, suggesting a very recent origin or a selection event, presumably associated with a determinate growth habit, which predominates in this race. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

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.     

Identification of loci controlling adaptation in Chinese soya bean landraces via a combination of conventional and bioclimatic GWAS

Landraces often contain genetic diversity that has been lost in modern cultivars, including alleles that confer enhanced local adaptation. To comprehensively identify loci associated with adaptive traits in soya bean landraces, for example flowering time, a population of 1938 diverse landraces and 97 accessions of the wild progenitor of cultivated soya bean, Glycine soja was genotyped using tGBS^®. Based on 99 085 high‐quality SNPs, landraces were classified into three sub‐populations which exhibit geographical genetic differentiation. Clustering was inferred from STRUCTURE, principal component analyses and neighbour‐joining tree analyses. Using phenotypic data collected at two locations separated by 10 degrees of latitude, 17 trait‐associated SNPs (TASs) for flowering time were identified, including a stable locus Chr12:5914898 and previously undetected candidate QTL/genes for flowering time in the vicinity of the previously cloned flowering genes, E1 and E2. Using passport data associated with the collection sites of the landraces, 27 SNPs associated with adaptation to three bioclimatic variables (temperature, daylength, and precipitation) were identified. A series of candidate flowering genes were detected within linkage disequilibrium (LD) blocks surrounding 12 bioclimatic TASs. Nine of these TASs exhibit significant differences in flowering time between alleles within one or more of the three individual sub‐populations. Signals of selection during domestication and/or subsequent landrace diversification and adaptation were detected at 38 of the 44 flowering and bioclimatic TASs. Hence, this study lays the groundwork to begin breeding for novel environments predicted to arise following global climate change.     

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.     

Tagging the signatures of domestication in common bean (Phaseolus vulgaris) by means of pooled DNA samples

BACKGROUND AND AIMS: The main aim of this study was to use an amplified fragment length polymorphism (AFLP)-based, large-scale screening of the whole genome of Phaseolus vulgaris to determine the effects of selection on the structure of the genetic diversity in wild and domesticated populations. METHODS: Using pooled DNA samples, seven each of wild and domesticated populations of P. vulgaris were studied using 2506 AFLP markers (on average, one every 250 kb). About 10 % of the markers were also analysed on individual genotypes and were used to infer allelic frequencies empirically from bulk data. In both data sets, tests were made to determine the departure from neutral expectation for each marker using an F(ST)-based method. KEY RESULTS: The most important outcome is that a large fraction of the genome of the common bean (16 %; P < 0.01) appears to have been subjected to effects of selection during domestication. Markers obtained in individual genotypes were also mapped and classified according to their proximities to known genes and quantitative trait loci (QTLs) of the domestication syndrome. Most of the markers that were found to be potentially under the effects of selection were located in the proximity of previously mapped genes and QTLs related to the domestication syndrome. CONCLUSIONS: Overall, the results indicate that in P. vulgaris a large portion of the genome appears to have been subjected to the effects of selection, probably because of linkage to the loci selected during domestication. As most of the markers that are under the effects of selection are linked to known loci related to the domestication syndrome, it is concluded that population genomics approaches are very efficient in detecting QTLs. A method based on bulk DNA samples is presented that is effective in pre-screening for a large number of markers to determine selection signatures.     

Population structure and genetic differentiation associated with breeding history and selection in tomato (Solanum lycopersicum L.)

Tomato (Solanum lycopersicum L.) has undergone intensive selection during and following domestication. We investigated population structure and genetic differentiation within a collection of 70 tomato lines representing contemporary (processing and fresh-market) varieties, vintage varieties and landraces. The model-based Bayesian clustering software, STRUCTURE, was used to detect subpopulations. Six independent analyses were conducted using all marker data (173 markers) and five subsets of markers based on marker type (single-nucleotide polymorphisms, simple sequence repeats and insertion/deletions) and location (exon and intron sequences) within genes. All of these analyses consistently separated four groups predefined by market niche and age into distinct subpopulations. Furthermore, we detected at least two subpopulations within the processing varieties. These subpopulations correspond to historical patterns of breeding conducted for specific production environments. We found no subpopulation within fresh-market varieties, vintage varieties and landraces when using all marker data. High levels of admixture were shown in several varieties representing a transition in the demarcation between processing and fresh-market breeding. The genetic clustering detected by using the STRUCTURE software was confirmed by two statistics, pairwise F(st) (θ) and Nei's standard genetic distance. We also identified a total of 19 loci under positive selection between processing, fresh-market and vintage germplasm by using an F(st)-outlier method based on the deviation from the expected distribution of F(st) and heterozygosity. The markers and genome locations we identified are consistent with known patterns of selection and linkage to traits that differentiate the market classes. These results demonstrate how human selection through breeding has shaped genetic variation within cultivated tomato.     

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.     

Determination of the population structure of common bean (Phaseolus vulgaris L.) accessions using lipoxygenase and resistance gene analog markers

The common bean (Phaseolus vulgaris L.) is an important food legume throughout the world. Because of the conservation across different plant species, it is possible to evaluate the degree of genetic diversity in the common bean using gene-based marker techniques. The lipoxygenase (LOX) and resistance gene analog (RGA) genes play an important role in the response to biotic and abiotic stresses. Eighty-six common bean accessions were genotyped using gene-based LOX and RGA markers. The total number of polymorphic bands ranged from 193 for LOX to 17 for RGA markers. We detected considerable diversity with a mean of 8.7 alleles per primer for the LOX analysis. For the RGA markers, the number of alleles per polymorphic locus varied from 1 to 4 with an average allele number of 2.8. The genetic similarity between the accessions based on the LOX and RGA markers ranged from 0.12 to 0.55. Using STRUCTURE, 3 groups were revealed among the accessions. The results of this study should provide valuable information for future studies on the genetic diversity of common bean accessions and for association mapping studies examining the relationships between the genotypic and phenotypic traits related to the stress response.     

Explorative genome scan to detect candidate loci for adaptation along a gradient of altitude in the common frog (Rana temporaria)

Today, with the rapid development of population genomics, the genetic basis of adaptation can be unraveled directly at the genome level, without any prerequisites about the selectively advantageous genes or traits. For nonmodel species, it is now possible to screen many markers randomly scattered across the genome and to distinguish between the neutral genetic background and outlier loci displaying an atypical behavior (e.g., a higher differentiation between populations). This study investigated the genetic frame of adaptation to a gradient of altitude in the common frog (Rana temporaria) by means of a genome scan based on 392 amplified fragment length polymorphism markers. Using two outlier detection methods never applied to dominant data so far, we sought for loci with a genetic differentiation diverging from neutral expectations when comparing populations from different altitudes. All the detected loci were sorted out according to their most probable cause for outlier behavior and classified as false positives, outliers due to local effects, or outliers associated with altitude. Altogether, eight good candidate loci were identified as potentially involved in adaptation to altitude because they were picked out in several independent interaltitude comparisons. This result illustrated the potential of genome-wide surveys to reveal selection signatures along selection gradients, where the association between environmental variables and fitness-related traits may be complex and/or cryptic. In this article, we also underlined the need for confirmation of the selection footprints for the outlier loci. Finally, we provided some preliminary insights into the genetic basis of adaptation along an altitudinal cline in the common frog.     

AFLP genome scans suggest divergent selection on colour patterning in allopatric colour morphs of a cichlid fish

Genome scan-based tests for selection are directly applicable to natural populations to study the genetic and evolutionary mechanisms behind phenotypic differentiation. We conducted AFLP genome scans in three distinct geographic colour morphs of the cichlid fish Tropheus moorii to assess whether the extant, allopatric colour pattern differentiation can be explained by drift and to identify markers mapping to genomic regions possibly involved in colour patterning. The tested morphs occupy adjacent shore sections in southern Lake Tanganyika and are separated from each other by major habitat barriers. The genome scans revealed significant genetic structure between morphs, but a very low proportion of loci fixed for alternative AFLP alleles in different morphs. This high level of polymorphism within morphs suggested that colour pattern differentiation did not result exclusively from neutral processes. Outlier detection methods identified six loci with excess differentiation in the comparison between a bluish and a yellow-blotch morph and five different outlier loci in comparisons of each of these morphs with a red morph. As population expansions and the genetic structure of Tropheus make the outlier approach prone to false-positive signals of selection, we examined the correlation between outlier locus alleles and colour phenotypes in a genetic and phenotypic cline between two morphs. Distributions of allele frequencies at one outlier locus were indeed consistent with linkage to a colour locus. Despite the challenges posed by population structure and demography, our results encourage the cautious application of genome scans to studies of divergent selection in subdivided and recently expanded populations.     

Marker-based linkage map of Andean common bean (Phaseolus vulgaris L.) and mapping of QTLs underlying popping ability traits

ABSTRACT: BACKGROUND: Nuna bean is a type of ancient common bean (Phaseolus vulgaris L.) native to the Andean region of South America, whose seeds possess the unusual property of popping. The nutritional features of popped seeds make them a healthy low fat and high protein snack. However, flowering of nuna bean only takes place under short-day photoperiod conditions, which means a difficulty to extend production to areas where such conditions do not prevail. Therefore, breeding programs of adaptation traits will facilitate the diversification of the bean crops and the development of new varieties with enhanced healthy properties. Although the popping trait has been profusely studied in maize (popcorn), little is known about the biology and genetic basis of the popping ability in common bean. To obtain insights into the genetics of popping ability related traits of nuna bean, a comprehensive quantitative trait loci (QTL) analysis was performed to detect single-locus and epistatic QTLs responsible for the phenotypic variance observed in these traits. RESULTS: A mapping population of 185 recombinant inbred lines (RILs) derived from a cross between two Andean common bean genotypes was evaluated for three popping related traits, popping dimension index (PDI), expansion coefficient (EC), and percentage of unpopped seeds (PUS), in five different environmental conditions. The genetic map constructed included 193 loci across 12 linkage groups (LGs), covering a genetic distance of 822.1 cM, with an average of 4.3 cM per marker. Individual and multi-environment QTL analyses detected a total of nineteen single-locus QTLs, highlighting among them the co-localized QTLs for the three popping ability traits placed on LGs 3, 5, 6, and 7, which together explained 24.9, 14.5, and 25.3 % of the phenotypic variance for PDI, EC, and PUS, respectively. Interestingly, epistatic interactions among QTLs have been detected, which could have a key role in the genetic control of popping. CONCLUSIONS: The QTLs here reported constitute useful tools for marker assisted selection breeding programs aimed at improving nuna bean cultivars, as well as for extending our knowledge of the genetic determinants and genotype x environment interaction involved in the popping ability traits of this bean crop.     

Interpreting genotype × environment interaction in tropical maize using linked molecular markers and environmental covariables

An understanding of the genetic and environmental basis of genotype×environment interaction (GEI) is of fundamental importance in plant breeding. In mapping quantitative trait loci (QTLs), suitable genetic populations are grown in different environments causing QTLs×environment interaction (QEI). The main objective of the present study is to show how Partial Least Squares (PLS) regression and Factorial Regression (FR) models using genetic markers and environmental covariables can be used for studying QEI related to GEI. Biomass data were analyzed from a multi-environment trial consisting of 161 lines from a F_3:4 maize segregating population originally created with the purpose of mapping QTLs loci and investigating adaptation differences between highland and lowland tropical maize. PLS and FR methods detected 30 genetic markers (out of 86) that explained a sizeable proportion of the interaction of maize lines over four contrasting environments involving two low-altitude sites, one intermediate-altitude site, and one high-altitude site for biomass production. Based on a previous study, most of the 30 markers were associated with QTLs for biomass and exhibited significant QEI. It was found that marker loci in lines with positive GEI for the highland environments contained more highland alleles, whereas marker loci in lines with positive GEI for intermediate and lowland environments contained more lowland alleles. In addition, PLS and FR models identified maximum temperature as the most-important environmental covariable for GEI. Using a stepwise variable selection procedure, a FR model was constructed for GEI and QEI that exclusively included cross products between genetic markers and environmental covariables. Higher maximum temperature in low- and intermediate-altitude sites affected the expression of some QTLs, while minimum temperature affected the expression of other QTLs. Received: 10 January 1999 / Accepted: 12 March 1999     

Mapping quantitative trait loci for a common bean (Phaseolus vulgaris L.) ideotype.

Breeding a model plant that encompasses individual traits thought to enhance yield potential, known as ideotype breeding, has traditionally focused on phenotypic selection of plants with desirable morphological traits. Broadening this breeding method to the molecular level through the use of molecular markers would avoid the environmental interactions associated with phenotypic selection. A population of 110 F5 recombinant inbred lines (RILs), derived from the cross between WO3391 and 'OAC Speedvale', was used to develop a genetic linkage map consisting of 105 random amplified polymorphic DNA (RAPD), simple sequence repeat (SSR), and sequence-tagged site (STS) markers. The map has a total length of 641 cM distributed across 8 linkage groups (LGs). Five of them were aligned on the core linkage map of bean. Twenty-one quantitative trait loci (QTLs) were identified over three environments for eight agronomic and architectural traits previously defined for a bean (Phaseolus vulgaris L.) ideotype. The QTLs were mapped to seven LGs with several regions containing QTLs for multiple traits. At least one QTL was located for each trait and a maximum of four were associated with lodging. Total explained phenotypic variance ranged from 10.6% for hypocotyl diameter to 45.4% for maturity. Some of the QTLs identified will be useful for early generation selection of tall, upright, high-yielding lines in a breeding program.     

Toward an Integrated Linkage Map of Common Bean. III. Mapping Genetic Factors Controlling Host-Bacteria Interactions

Restriction fragment length polymorphism (RFLP)-based genetic linkage maps allow us to dissect the genetic control of quantitative traits (QT) by locating individual quantitative trait loci (QTLs) on the linkage map and determining their type of gene action and the magnitude of their contribution to the phenotype of the QT. We have performed such an analysis for two traits in common bean, involving interactions between the plant host and bacteria, namely Rhizobium nodule number (NN) and resistance to common bacterial blight (CBB) caused by Xanthomonas campestris pv. phaseoli. Analyses were conducted in the progeny of a cross between BAT93 (fewer nodules; moderately resistant to CBB) and Jalo EEP558 (more nodules; susceptible to CBB). An RFLP-based linkage map for common bean based on 152 markers had previously been derived in the F(2) of this cross. Seventy F(2)-derived F(3) families were inoculated in separate greenhouse experiments with Rhizobium tropici strain UMR1899 or X. c. pv. phaseoli isolate isolate W18. Regression and interval mapping analyses were used to identify genomic regions involved in the genetic control of these traits. These two methods identified the same genomic regions for each trait, with a few exceptions. For each trait, at least four putative QTLs were identified, which accounted for approximately 50% and 75% of the phenotypic variation in NN and CBB resistance, respectively. A chromosome region on linkage group D7 carried factor(s) influencing both traits. In all other cases, the putative QTLs affecting NN and CBB were located in different linkage groups or in the same linkage group, but far apart (more than 50 cM). Both BAT93 and Jalo EEP558 contributed alleles associated with higher NN, whereas CBB resistance was always associated with BAT93 alleles. Further investigations are needed to determine whether the QTLs for NN and CBB on linkage group D7 represent linked genes or the same gene with pleiotropic effects. Identification of the QTLs raises the possibility of initiating map-based cloning and marker-assisted selection for these traits.     

Outlier analyses to test for local adaptation to breeding grounds in a migratory arctic seabird

Investigating the extent (or the existence) of local adaptation is crucial to understanding how populations adapt. When experiments or fitness measurements are difficult or impossible to perform in natural populations, genomic techniques allow us to investigate local adaptation through the comparison of allele frequencies and outlier loci along environmental clines. The thick‐billed murre (Uria lomvia) is a highly philopatric colonial arctic seabird that occupies a significant environmental gradient, shows marked phenotypic differences among colonies, and has large effective population sizes. To test whether thick‐billed murres from five colonies along the eastern Canadian Arctic coast show genomic signatures of local adaptation to their breeding grounds, we analyzed geographic variation in genome‐wide markers mapped to a newly assembled thick‐billed murre reference genome. We used outlier analyses to detect loci putatively under selection, and clustering analyses to investigate patterns of differentiation based on 2220 genomewide single nucleotide polymorphisms (SNPs) and 137 outlier SNPs. We found no evidence of population structure among colonies using all loci but found population structure based on outliers only, where birds from the two northernmost colonies (Minarets and Prince Leopold) grouped with birds from the southernmost colony (Gannet), and birds from Coats and Akpatok were distinct from all other colonies. Although results from our analyses did not support local adaptation along the latitudinal cline of breeding colonies, outlier loci grouped birds from different colonies according to their non‐breeding distributions, suggesting that outliers may be informative about adaptation and/or demographic connectivity associated with their migration patterns or nonbreeding grounds.     

Dissection of the genetic architecture of body weight in chicken reveals the impact of epistasis on domestication traits

In this contribution, we study the genetic mechanisms leading to differences in the observed growth patterns of domesticated White Leghorn chickens and their wild ancestor the red jungle fowl. An epistatic QTL analysis for several body-weight measures from hatch to adulthood confirms earlier findings that polymorphisms at >15 loci contribute to body-weight determination in an F(2) intercross between these populations and that many loci are involved in complex genetic interactions. Here, we use a new genetic model to decompose the genetic effects of this multilocus epistatic genetic network. The results show how the functional modeling of genetic effects provides new insights into how genetic interactions in a large set of loci jointly contribute to phenotypic expression. By exploring the functional effects of QTL alleles, we show that some alleles can display temporal shifts in the expression of genetic effects due to their dependencies on the genetic background. Our results demonstrate that the effects of many genes are dependent on genetic interactions with other loci and how their involvement in the domestication process relies on these interactions.     

Identification of single nucleotide polymorphisms and haplotypes associated with yield and yield components in soybean (<Emphasis Type="Italic">Glycine max</Emphasis>) landraces across multiple environments

Genome-wide association analysis is a powerful approach to identify the causal genetic polymorphisms underlying complex traits. In this study, we evaluated a population of 191 soybean landraces in five environments to detect molecular markers associated with soybean yield and its components using 1,536 single-nucleotide polymorphisms (SNPs) and 209 haplotypes. The analysis revealed that abundant phenotypic and genetic diversity existed in the studied population. This soybean population could be divided into two subpopulations and no or weak relatedness was detected between pair-wise landraces. The level of intra-chromosomal linkage disequilibrium was about 500 kb. Genome-wide association analysis based on the unified mixed model identified 19 SNPs and 5 haplotypes associated with soybean yield and yield components in three or more environments. Nine markers were found co-associated with two or more traits. Many markers were located in or close to previously reported quantitative trait loci mapped by linkage analysis. The SNPs and haplotypes identified in this study will help to further understand the genetic basis of soybean yield and its components, and may facilitate future high-yield breeding by marker-assisted selection in soybean.     

Genetic Variation and Breeding Signature in Mass Selection Lines of the Pacific Oyster (Crassostrea gigas) Assessed by SNP Markers

In breeding industries, a challenging problem is how to keep genetic diversity over generations. To investigate genetic variation and identify breeding signatures in mass selected lines of Pacific oyster (Crassostrea gigas), three sixth-generation selected lines and four wild populations were assessed using 103 single nucleotide polymorphism (SNP) markers. The genetic diversity data indicated that the selected lines exhibited a significant reduction in the observed heterozygosity and observed number of alleles per locus compared with the wild populations (P≤0.05), indicating the selected lines tended to lose genetic diversity contrasted with the wild populations. The unweighted pair-group method with arithmetic mean (UPGMA) analysis showed that the wild populations and selected lines were not separated into two groups. Using four outlier tests, a total of 17 loci were found under selection at two levels. The global outlier detection suggested that 4 common outlier loci were subject to selection using both the hierarchical island model and Bayesian likelihood approaches. At regional level, 3 SNPs were detected as outlier using at least two outlier tests and one outlier SNP (CgSNP309) was overlapped in the two wild-selected population comparisons. The candidate outlier SNPs provide valuable resources for future association studies in C. gigas.