Next Article in Journal
Polyoxypregnane Ester Derivatives and Lignans from Euphorbia gossypina var. coccinea Pax.
Previous Article in Journal
Cytotoxic and Genotoxic Evaluation of Biosynthesized Silver Nanoparticles Using Moringa oleifera on MCF-7 and HUVEC Cell Lines
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Plants
Volume 11
Issue 10
10.3390/plants11101298
plants-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Development of a Model for Genomic Prediction of Multiple Traits in Common Bean Germplasm, Based on Population Structure

by
Jing Shao
1,2,3,
Yangfan Hao
2,
Lanfen Wang
2,
Yuxin Xie
2,
Hongwei Zhang
2,
Jiangping Bai
1,3,*,
Jing Wu
1,2,* and
Junjie Fu
2,*
1
Department of Crop Genetics and Breeding, College of Agronomy, Gansu Agricultural University, Lanzhou 730070, China
2
Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing 100081, China
3
Gansu Provincial Key Laboratory of Aridland Crop Science, Lanzhou 730070, China
*
Authors to whom correspondence should be addressed.
Plants 2022, 11(10), 1298; https://doi.org/10.3390/plants11101298
Submission received: 25 March 2022 / Revised: 8 May 2022 / Accepted: 9 May 2022 / Published: 12 May 2022
(This article belongs to the Section Plant Genetics, Genomics and Biotechnology)
Browse Figures
Review Reports Versions Notes

Abstract

:
Due to insufficient identification and in-depth investigation of existing common bean germplasm resources, it is difficult for breeders to utilize these valuable genetic resources. This situation limits the breeding and industrial development of the common bean (Phaseolus vulgaris L.) in China. Genomic prediction (GP) is a breeding method that uses whole-genome molecular markers to calculate the genomic estimated breeding value (GEBV) of candidate materials and select breeding materials. This study aimed to use genomic prediction to evaluate 15 traits in a collection of 628 common bean lines (including 484 landraces and 144 breeding lines) to determine a common bean GP model. The GP model constructed by landraces showed a moderate to high predictive ability (ranging from 0.59–0.88). Using all landraces as a training set, the predictive ability of the GP model for most traits was higher than that using the landraces from each of two subgene pools, respectively. Randomly selecting breeding lines as additional training sets together with landrace training sets to predict the remaining breeding lines resulted in a higher predictive ability based on principal components analysis. This study constructed a widely applicable GP model of the common bean based on the population structure, and encouraged the development of GP models to quickly aggregate excellent traits and accelerate utilization of germplasm resources.

1. Introduction

The common bean (Phaseolus vulgaris L.) is an edible legume with the widest and largest cultivation area in the world [1]. It is an important source of plant protein for humans and is high in protein, has a medium starch content, is low in fat, and is rich in mineral nutrients [2]. The generally low yield of the common bean limits its advantages; therefore, increasing the production will have an important impact on improving the nutrition and health of consumers, especially in developing countries [2]. Hence, genetic research and breeding improvement of bean germplasms are essential. The genetic resources of the common bean mainly consist of two highly differentiated gene pools: the Andean (An) and the Mesoamerican (M) [3,4]. These two gene pools contain abundant germplasm resources, including landraces and breeding lines.
During the domestication of the common bean, landraces on which humans depend for cultivation and diet are often grown in low-input production systems. They have a wide planting range and strong adaptability [5]. Alleles for endemic disease resistance and tolerance to major climatic stresses have also been reported [6]. A broad genetic basis for landraces has been identified in germplasm evaluation studies, and the diversity of breeding lines has narrowed through breeding and domestication [7,8]. Therefore, common bean landraces are a resource for increasing the genetic diversity of breeding lines. This suggests the possibility of predicting and identifying important agronomic traits and effective gene loci [9] using the genetic variation present in the existing common bean germplasm. A population with rich genetic diversity could serve as an effective training set in genetic assessment models based on molecular genetic markers.
Genomic prediction (GP) [10] has been widely used for food crops, such as barley [11], wheat [12], maize [13], and rice [14]. The construction of prediction models provides a good method for selecting promising varieties. It has three main advantages: it increases the genetic gain, shortens the breeding period, and reduces capital costs. In recent years, an increasing number of studies have focused on the GP of common bean. More advanced research mainly focuses on the breeding and development of the common bean, and the construction model of related single nucleotide polymorphisms (SNPs) is determined by combining genome-wide association analysis (GWAS) and GP [15,16,17,18]. Other studies have mainly focused on the genomic prediction of a few traits in common bean cultivars using different models [19,20,21]. Several studies have been conducted to improve model predictive ability by considering genotype-by-environment interactions (G×E) [22,23,24]. In particular, the common bean is a species with a strong population structure. In existing studies, population genetic structure has been considered in the GP of several types of crops with well-defined genetic population distribution and less domestication, such as wheat [25,26], rice [25], maize [27], and oats [28], to improve the predictive power. However, studies about modeling the distribution of common bean germplasm resources are limited.
In this study, we optimized the model’s precision to predict breeding lines by exploring the composition of the training set in the GP of the common bean (Figure 1). Fifteen traits, i.e., days to flowering (DF), days to maturity (DM), plant height (PH), stem node number (SNN), branch number (BN), pod per plant (PP), pod length (PL), pod width (PW), pod height (PDH), seeds per pod (SP), grain yield per plant (GY), 100-grain weight (GW), seed length (SL), seed width (SW), and seed height (SH), were used. The description of each trait was included in Supplementary Table S1 [29]. First, prediction abilities were measured in two gene pools of common bean landraces using a ten-fold cross-validation repeated 100 times. Second, we analyzed the ability to combine or assign the landraces of the two gene pools to predict the breeding lines of both gene pools. Third, we considered the ability to add parts of breeding lines as a training set to predict the value of the remaining breeding lines. This study provides new perspectives for common bean breeding and offers a reference strategy for other crops with a strong population structure.

2. Results

2.1. Evaluation of Population Size and Marker Number on the Model’s Precision

Fifteen traits were used to evaluate the model’s predictive ability among 484 common bean landraces from the germplasm collection at the Chinese Academy of Agriculture Sciences (CAAS). We performed cross-validation using the ridge regression best linear unbiased prediction (RR-BLUP) model for landraces from two common bean gene pools, respectively. Multiple gradients of population size and number of SNPs were designed to examine the predictive ability of the combined landraces for multiple traits. This study found that the predictive ability of the landraces from the two gene pools differed greatly in different traits (0.35–0.82) (Supplementary Figure S1, Supplementary Table S2), and that of the combined landraces was higher (0.59–0.88) (Figure 2, Supplementary Figure S2). As shown in these figures, the predictive ability could be improved with an increase in the population size and number of SNPs, but it no longer showed significant improvement when it increased to a certain number. The optimal population size (opt_pop) and optimal number of SNPs (opt_snp) differed for each trait. In particular, the opt_pop of DF is 40% and opt_snp is 2000 (0.67), the opt_pop of PH is 40% and opt_snp is 1000 (0.84), the opt_pop of SP is 40% and opt_snp is 500 (0.77), and the opt_pop of SH is 40% and opt_snp is 500 (0.81) (Supplementary Table S3). However, the predictive ability of each trait ranged from 0.56 for GY to 0.85 for GW at the opt_pop and opt_snp. These results illustrated that landraces as a training set could obtain better predictive ability for all investigated traits, which lays the foundation for population selection to optimize the prediction of breeding lines.

2.2. Ability of Landrace Subgene Pools in Predicting Breeding Lines

We designed a pattern of landrace predicting breeding lines based on two gene pools (An and M) of the common bean. The ability of landraces to predict breeding lines ranged from 0.54 to 0.82, except for two traits that were more affected by photoperiod, DF, and DM (0.32 and 0.21, respectively) (Table 1). A few differences were found when comparing the ability of landraces (An + M) and landraces (An) to predict breeding lines (An). This ability was also comparable when landraces (An + M) and landraces (M) were used to predict breeding lines (M; Supplementary Figure S3). Subsequently, we counted the accessions corresponding with the predicted and observed values from the top 30 and top 30% for each trait. The number of the same accessions was recorded as the ratio (predicted/observed) of the top 30 and the top 30% of the number of materials. The results showed that the ratios differed for each trait. However, there was no significant difference in the ratios of predicted breeding lines (An) using landrace (An + M), or predicted breeding lines (An) using landrace (An) for the same trait. Similar results were found for the Mesoamerican gene pools (Supplementary Table S4). Since combining and assigning gene pools had no significant effect on the landraces’ ability to predict breeding lines, we do not recommend assigning landrace pools as the training set in the subsequent prediction studies.

2.3. Optimization of Landrace Training Sets by Adding Breeding Lines

The genetic data of the common bean were transformed using principal components analysis (PCA) into components explaining most of the genetic variations. The first three axes of PCA explained 63.23% of the variability in the entire collection. The two gene pools of the common bean, Andean and Mesoamerican, were undoubtedly distinguished as two clusters (Figure 3). Based on the common bean germplasm, we designed a two-step scenario to optimize the composition of the training set by adding 50% of the breeding lines (Table 2). Firstly, scenario 1 investigated the impact of adding breeding lines as the training set on the predictive ability of the two gene pools. In the prediction using 50% Andean gene pool breeding lines (A2) as the testing set, the predictive ability of the model was improved in all traits to varying degrees compared to when landraces and 50% breeding lines (A1 or A1 + M1) were used as the training set. Adding 50% breeding lines (M1 or A1 + M1) to the prediction using 50% Mesoamerican gene pool breeding lines (M2) as the testing set also improved the model’s precision. However, the addition of breeding lines was more stable for improving the Andean model’s precision (Supplementary Figure S4). Then, scenario 2 made a prediction comparison on how to select 50% of the individuals in the breeding lines. The results showed that compared with the 50% breeding lines randomly selected as the additional training set, the 50% breeding lines (A1 + M1) extracted from the subgene pool as the additional training set had a more stable and improved predictive ability of the model, especially for the traits PL (+0.12) and PW (+0.25) (Table 3).

3. Discussion

Landraces are accessions with rich genetic diversity and are widely used in the genomic prediction of crop germplasms. Moderate to high prediction abilities were obtained for accessions such as wheat [30,31], maize [32], and white lupin [33]. Our results are in concordance with those of previous GP studies on landraces. The model’s precision ranged from 0.59 for GY to 0.88 for PH. We evaluated the model’s precision when the landraces were used as training sets and breeding lines were used as testing sets. The GP model built by combining the gene pool landraces (An + M) did not significantly improve the ability to predict breeding lines (An) or breeding lines (M). We speculated that the main reason could be that the common bean is domesticated from wild relatives, which are characterized by their inhabitation in a relatively narrow ecological niche [34]. Andean and Mesoamerican gene pools are probably two independent origin centers of the common bean before domestication, which resulted in the inability to communicate genetic information between the two gene pools [35]. Therefore, the model’s precision was not substantially improved when the training set consisted of landraces (An + M) that combined the two gene pools. Subsequently, the kinship between the training set (landraces) and the testing set (breeding lines) was examined using PCA (Figure 3). When a part of the breeding lines was added to the training set as an additional training set, the model’s precision was improved across multiple traits. It was observed that the breeding lines selected based on the gene pool resulted in more stable and more accurate improvements of the model than the randomly selected breeding lines.
The genetic relationship between the training sets and testing sets is crucial when using GP to select promising breeding materials from a germplasm collection. Understanding the population structure of germplasm resources and ensuring the distribution of genetically related individuals in the training sets and testing sets could solve the problem of genetic relationships in structured populations and maximize predictive ability [26,27,36]. In this study, scenario 1 indicated that we artificially increased the relevance between the training sets and the testing sets by adding a part of the breeding lines as the additional training sets. The GP model predicted the remaining breeding lines more accurately. In the Andean gene pool, compared to G1, the prediction abilities of G2 and G3 were improved, especially in GW (+0.27, +0.27) and SW (+0.23, +0.24). In the Mesoamerican gene pool, the prediction abilities of G5 and G6 were improved for traits such as DF (+0.18, +0.17) and SNN (+0.14, +0.14) compared to that of G4 (Table 3). The PCA indicated that accessions within the Andean gene pool were more highly related to each other than accessions within the Mesoamerican gene pool (Figure 3). The prediction abilities for different traits in the two gene pools were related to the genetic correlation between the germplasm in each subgene pool. Similarly, in a genomic prediction study of four traits (heading, height, biomass, and yield) in oats, it has been found that the multi-group training sets will obtain higher predictive ability (0.32–0.87) than across-group scenarios (−0.55–0.27) [28].
The population structure plays a critical role in optimizing the composition of the training set. In a rice population with a strong structure, stratified sampling has higher accuracy than other random sampling methods for four traits: florets per panicle, flowering time, plant height, and protein content [25]. The stratified sampling here was based on several clusters of the rice population structure in PCA. This was the same as the sub-gene pool sampling of the common bean population in scenario 2 in this study. For all traits, adding extracted breeding lines from the subgene pool as an additional training set is beneficial for improving the ability to predict the remaining breeding lines. The predictive ability of G7 and G8 showed a remarkable increase for two traits, PL and PW, compared to that of other traits. Our findings are consistent with those of previous studies of GP in similarly structured populations. For example, in spring wheat germplasm, Muleta et al. [26] purposefully divided accessions into two major subpopulations through population structure analysis. On this basis, the training sets were selected to predict the testing sets, and promising results (≥0.7) were obtained for disease resistance traits such as infection type and disease severity. In the oat population, the prediction abilities of the testing sets with different training set organizations showed significant differences. For oat yield [28], selection based on population structure could increase the predictive ability of the training set to the testing set when compared with random selection. Yield prediction was improved by 0.35 and 0.03 in two predicted subpopulations of oats.
However, there is evidence that considering the population structure to optimize the training set can improve the model’s precision [37,38,39]. For the first time, the optimized training set was used in two major gene pools to predict the breeding lines in common bean germplasm resources. This study provides new ideas for utilizing common bean germplasm resources in China and a new breeding strategy for similar small crops in the primary domestication stage. The effective use of the collected genetic diversity data provided help for further breeding and research on the common bean.

4. Materials and Methods

4.1. Plant Materials

The common bean dataset used in this study, i.e., the Chinese Academy of Agriculture Sciences (CAAS) collection determined and evaluated by Wu et al. [29], comprised 628 germplasm resources. This included 484 landraces (including 223 from the Andean gene pool and 261 from the Mesoamerican gene pool) and 144 breeding lines (including 60 from the Andean gene pool and 84 from the Mesoamerican gene pool).

4.2. Genotyping

Genotyping data for 4.8 million SNPs distributed over all 11 chromosomes published in an article by Wu et al. [29] were used to investigate the model’s precision for 628 common bean landraces and breeding lines. In this study, we have processed these data for genomic prediction analysis. After filtering on the minor allele frequency (>5%) and the percentage of missing values (<50%) using VCFtools version 0.1.13 (https://github.com/vcftools/vcftools, accessed 2 March 2021), the sample function in R version 4.0.2 (https://www.r-project.org/) was used to sample 11 chromosomes according to the proportion of SNPs, and 9781 SNP markers were obtained.

4.3. Phenotyping

Phenotypic data for 15 traits in the common bean dataset from the published article by Wu et al. [29] were analyzed: days to flowering (DF), days to maturity (DM), plant height (PH), stem node number (SNN), branch number (BN), pod per plant (PP), pod length (PL), pod width (PW), pod height (PDH), seeds per pod (SP), grain yield per plant (GY), 100-grain weight (GW), seed length (SL), seed width (SW), and seed height (SH). The experiments were carried out from 2014 to 2016 at three locations (i.e., Bijie, Harbin, and Sanya in China) [29]. Phenotypic data for each trait were calculated from the average of nine location–year combinations.

4.4. Genomic Prediction Model and Analysis

The genomic estimated breeding values (GEBVs) were calculated using the ridge regression best linear unbiased prediction (RR-BLUP) of the GP model. The model assumes that the effects of all SNPs have a common variance, and a stable predictive ability could be obtained in most traits (especially traits influenced by a large number of minor genes), which is implemented in the R package “rrBLUP”. The rrBLUP package [40] can solve any mixed model of the form:
y = X β + Z u + ε
where y is an N × 1 vector of observed phenotypes (N is the size of the training set), X is the identity matrix for the fixed effects β , Z   is an N × M genotypic matrix for the random effects u with u ~ N 0 , I σ u 2 (M is the number of SNPs), and ε is the vector of residuals with ε ~ N 0 , I σ u 2 , where I is the identity matrix. In the Z matrix, {0,1,2} for biallelic SNP markers of genotypes (i.e., 0 and 2 for homozygotes and 1 for heterozygotes). The common bean is a highly self-pollinated crop [2] and its genotype data are composed of {0,1}.
In this study, we conducted three GP analyses based on this model. (1) Seven different population sizes (10, 20, 40, 50, 60, 80, and 100%) and eleven different numbers of SNPs (100, 500, 1000, 1500, 2000, 3000, 4000, 5000, 6000, 7000, and 9781) were used in a ten-fold cross-validation with 100 replications for 15 agronomic traits in 484 common bean landraces. (2) Combining gene pool prediction: landraces (An + M) were used as the training set to predict breeding lines (An) and breeding lines (M). Assigning gene pool prediction: landraces (An) were used as the training set to predict breeding lines (An), and landraces (M) were used as the training set to predict breeding lines (M). (3) Two scenarios were designed: adding 50% breeding lines in different gene pools as additional training sets to predict the remaining 50% of breeding lines in each gene pool, selecting by gene pool; or randomly selecting the added breeding lines to predict the remaining individuals. The predictive ability was calculated using the mean of Pearson’s correlation coefficient (r) between the estimated and observed values in the testing set [41]. The closer the r-value is to 1, the higher the model’s predictive ability. The figures and tables in the article were generated using the R package ggplot2 (https://cran.r-project.org/package=ggplot2, accessed on 2 March 2021) and Microsoft Excel 2019.

5. Conclusions

This study used existing common bean germplasm resources to perform genomic prediction in multiple traits. The research showed the following. (1) The number of markers and the size of the training population required to build a GP model for different traits of the common bean differ according to their heritability levels. For traits with high heritability, the number of markers required by the GP model and the size of the training population was smaller than those with low heritability. (2) Using the common bean as an example, the construction of the GP model for germplasm resource prediction using landraces and the optimization of the training population were analyzed. The main result was that adding a part of the breeding lines as an additional training set could improve the predictive ability of the GP model built by landraces on the remaining breeding lines. By analyzing the population structure of all accessions, the training population was optimized and the traits were predicted more accurately. This investigation provides an argument for mining and analyzing germplasm resources through GP for crops with a population structure similar to that of the common bean and provides a reference for germplasm selection for breeders.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/plants11101298/s1, Table S1. Abbreviations and description for fifteen traits. Figure S1. Ten-fold cross-validation with 100 replications was performed for 15 agronomic traits in landraces from two common bean gene pools, respectively. Table S2. Mean values of prediction abilities for 15 agronomic traits in the landraces from two common bean gene pools, respectively. Table S3. The model’s precision of 484 landraces using different numbers of SNPs and different population size for 15 traits. Figure S2. Different number of SNPs (100, 500, 1000, 1500, 2000, 3000, 4000, 5000, 6000, 7000, and 9781) and different population size (10, 20, 40, 50, 60, 80, and 100%) were used in ten-fold cross-validation with 100 replications for 11 agronomic traits in 484 common bean landraces. Figure S3. Histograms of predictive ability of landraces as a training set to predict breeding lines among 15 traits. Table S4. The table records the number of identical materials in the material corresponding to the top 30 (as well as the top 30%) predicted and observed values in the testing set. Figure S4. Histograms of the model’s precision for 8 groups of two scenarios among 15 traits.

Author Contributions

Conceptualization, Y.H. and J.F.; Data curation, L.W. and J.W.; Methodology, Y.H. and H.Z.; Software, J.S.; Validation, J.S. and Y.H.; Formal Analysis, J.F.; Investigation, Y.H. and J.S.; Resources, J.W. and L.W.; Writing—Original Draft Preparation, J.S.; Writing—Review and Editing, J.S., Y.H., Y.X., J.W. and J.F.; Visualization, J.S.; Supervision, J.F. and J.B.; Project Administration, J.F.; Funding Acquisition, J.F. and J.W. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Key Research and Development Program of China (2020YFE0202300); China Agriculture Research System of MOF and MARA—Food Legumes (CARS-08); and Agricultural Science and Technology Innovation Program of CAAS.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The drafter is most grateful to J.B. from Gansu Agricultural University and J.F., J.W. and L.W. from the Chinese Academy of Agricultural Sciences for the opportunity to study. Additionally, special thanks to the author Y.H. for her wonderful collaboration and patient support. Finally, the sincerest thanks go to the author Y.X. for her fluent use of grammar and careful revision.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kelly, J.D. Advances in common bean improvement: Some case histories with broader applications. Acta Hortic. 2004, 637, 99–122. [Google Scholar] [CrossRef]
  2. Broughton, W.J.; Hernandez, G.; Blair, M.; Beebe, S.; Gepts, P.; Vanderleyden, J. Beans (Phaseolus spp.)—Model food legumes. Plant Soil 2003, 252, 55–128. [Google Scholar] [CrossRef] [Green Version]
  3. Koenig, R.; Gepts, P. Segregation and linkage of genes for seed proteins, isozymes, and morphological traits in common bean (Phaseolus vulgaris). J. Hered. 1989, 80, 455–459. [Google Scholar] [CrossRef]
  4. Bitocchi, E.; Rau, D.; Bellucci, E.; Rodriguez, M.; Murgia, M.L.; Gioia, T.; Santo, D.; Nanni, L.; Alttene, G.; Papa, R. Beans (Phaseolus ssp) as a Model for Understanding Crop Evolution. Front. Plant Sci. 2017, 8, 722. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Tiranti, B.; Negri, V. Selective microenvironmental effects play a role in shaping genetic diversity and structure in a Phaseolus vulgaris L landrace: Implications for on-farm conservation. Mol. Ecol. 2007, 16, 4942–4955. [Google Scholar] [CrossRef]
  6. Carović-Stanko, K.; Liber, Z.; Vidak, M.; Barešić, A.; Grdiša, M.; Lazarević, B.; Šatović, Z. Genetic Diversity of Croatian Common Bean Landraces. Front. Plant Sci. 2017, 8, 604. [Google Scholar] [CrossRef] [Green Version]
  7. Gioia, T.; Logozzo, G.; Marzario, S.; Zeuli, P.S.; Gepts, P. Evolution of SSR diversity from wild types to US. Advanced cultivars in the Andean and Mesoamerican domestications of common bean (Phaseolus vulgaris). PLoS ONE 2019, 14, e0211342. [Google Scholar] [CrossRef] [Green Version]
  8. Zhou, Z.; Jiang, Y.; Wang, Z.; Gou, Z.; Lyu, J.; Li, W.; Yu, Y.; Shu, L.; Zhao, Y.; Ma, Y.; et al. Resequencing 302 wild and cultivated accessions identifies genes related to domestication and improvement in soybean. Nat. Biotechnol. 2015, 33, 408–414. [Google Scholar] [CrossRef] [Green Version]
  9. Wu, J. Research Progress on Common Bean Genomics Study. J. Sichuan Agric. Univ. 2021, 39, 4–10. [Google Scholar]
  10. Meuwissen, T.H.E.; Hayes, B.J.; Goddard, M.E. Prediction of total genetic value using genome-wide dense marker maps. Genetics 2001, 157, 1819–1829. [Google Scholar] [CrossRef]
  11. Zhong, S.; Dekkers, J.C.M.; Fernando, R.L.; Jannink, J.-L. Factors Affecting Accuracy From Genomic Selection in Populations Derived From Multiple Inbred Lines: A Barley Case Study. Genetics 2009, 182, 355–364. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Poland, J.A.; Endelman, J.; Dawson, J.; Rutkoski, J.; Wu, S.; Manes, Y.; Dreisigacker, S.; Crossa, J.; Sánchez-Villeda, H.; Sorrells, M.; et al. Genomic Selection in Wheat Breeding using Genotyping-by-Sequencing. Plant Genome 2012, 5, 103–113. [Google Scholar] [CrossRef] [Green Version]
  13. Riedelsheimer, C.; Czedik-Eysenberg, A.; Grieder, C.; Lisec, J.; Technow, F.; Sulpice, R.; Altmann, T.; Stitt, M.; Willmitzer, L.; Melchinger, A.E. Genomic and metabolic prediction of complex heterotic traits in hybrid maize. Nat. Genet. 2012, 44, 217–220. [Google Scholar] [CrossRef] [PubMed]
  14. Spindel, J.; Begum, H.; Akdemir, D.; Virk, P.; Collard, B.; Redona, E.; Atlin, G.; Jannink, J.L.; McCouch, S.R. Genomic Selection and Association Mapping in Rice (Oryza sativa): Effect of Trait Genetic Architecture, Training Population Composition, Marker Number and Statistical Model on Accuracy of Rice Genomic Selection in Elite, Tropical Rice Breeding Lines. PLoS Genet. 2015, 11, e1004982. [Google Scholar]
  15. Wen, L.; Chang, H.-X.; Brown, P.J.; Domier, L.L.; Hartman, G.L. Genome-wide association and genomic prediction identifies soybean cyst nematode resistance in common bean including a syntenic region to soybean Rhg1 locus. Hortic. Res. 2019, 6, 9. [Google Scholar] [CrossRef] [Green Version]
  16. Diaz, L.M.; Arredondo, V.; Ariza-Suarez, D.; Aparicio, J.; Buendia, H.F.; Cajiao, C.; Mosquera, G.; Beebe, S.E.; Mukankusi, C.M.; Raatz, B. Genetic Analyses and Genomic Predictions of Root Rot Resistance in Common Bean Across Trials and Populations. Front. Plant Sci. 2021, 12, 629221. [Google Scholar] [CrossRef]
  17. Diaz, S.; Ariza-Suarez, D.; Ramdeen, R.; Aparicio, J.; Arunachalam, N.; Hernandez, C.; Diaz, H.; Ruiz, H.; Piepho, H.-P.; Raatz, B. Genetic Architecture and Genomic Prediction of Cooking Time in Common Bean (Phaseolus vulgaris L.). Front. Plant Sci. 2021, 11, 2257. [Google Scholar] [CrossRef]
  18. Shi, A.; Gepts, P.; Song, Q.; Xiong, H.; Michaels, T.E.; Chen, S. Genome-Wide Association Study and Genomic Prediction for Soybean Cyst Nematode Resistance in USDA Common Bean (Phaseolus vulgaris) Core Collection. Front. Plant Sci. 2021, 12, 624156. [Google Scholar] [CrossRef]
  19. Barili, L.D.; Vale, N.M.D.; Silva, F.F.E.; Carneiro, J.E.D.S.; De Oliveira, H.R.; Vianello, R.P.; Valdisser, P.A.M.R.; Nascimento, M. Genome prediction accuracy of common bean via Bayesian models. Ciênc. Rural 2018, 48. [Google Scholar] [CrossRef]
  20. Nascimento, A.C.; Nascimento, M.; Azevedo, C.; Silva, F.; Barili, L.; Vale, N.; Carneiro, J.E.; Cruz, C.; Carneiro, P.C.; Serão, N. Quantile Regression Applied to Genome-Enabled Prediction of Traits Related to Flowering Time in the Common Bean. Agronomy 2019, 9, 796. [Google Scholar] [CrossRef] [Green Version]
  21. Rosado, R.; Cruz, C.; Barili, L.; Carneiro, J.D.S.; Carneiro, P.; Carneiro, V.; Da Silva, J.; Nascimento, M. Artificial Neural Networks in the Prediction of Genetic Merit to Flowering Traits in Bean Cultivars. Agriculture 2020, 10, 638. [Google Scholar] [CrossRef]
  22. Keller, B.; Ariza-Suarez, D.; De la Hoz, J.; Aparicio, J.S.; Portilla-Benavides, A.E.; Buendia, H.F.; Mayor, V.M.; Studer, B.; Raatz, B. Genomic Prediction of Agronomic Traits in Common Bean (Phaseolus vulgaris L.) Under Environmental Stress. Front. Plant Sci. 2020, 11, 1001. [Google Scholar] [CrossRef] [PubMed]
  23. Melo, V.L.D.; Marçal, T.D.S.; Rocha, J.R.A.S.D.C.; dos Anjos, R.S.R.; Carneiro, P.C.S.; Carneiro, J.E.D.S. Modeling (co)variance structures for genetic and non-genetic effects in the selection of common bean progenies. Euphytica 2020, 216, 1–13. [Google Scholar] [CrossRef]
  24. De Souza, M.H.; Pereira Júnior, J.D.; Steckling, S.D.M.; Mencalha, J.; Dias, F.D.S.; Rocha, J.R.D.A.S.D.C.; Carneiro, P.C.S.; Carneiro, J.E.D.S. Adaptability and stability analyses of plants using random regression models. PLoS ONE 2020, 15, e0233200. [Google Scholar]
  25. Isidro, J.; Jannink, J.-L.; Akdemir, D.; Poland, J.; Heslot, N.; Sorrells, M.E. Training set optimization under population structure in genomic selection. Theor. Appl. Genet. 2015, 128, 145–158. [Google Scholar] [CrossRef] [Green Version]
  26. Muleta, K.T.; Bulli, P.; Zhang, Z.; Chen, X.; Pumphrey, M. Unlocking Diversity in Germplasm Collections via Genomic Selection: A Case Study Based on Quantitative Adult Plant Resistance to Stripe Rust in Spring Wheat. Plant Genome 2017, 10, 1–15. [Google Scholar] [CrossRef] [Green Version]
  27. Cao, S.; Song, J.; Yuan, Y.; Zhang, A.; Ren, J.; Liu, Y.; Qu, J.; Hu, G.; Zhang, J.; Wang, C.; et al. Genomic Prediction of Resistance to Tar Spot Complex of Maize in Multiple Populations Using Genotyping-by-Sequencing SNPs. Front. Plant Sci. 2021, 12, 1438. [Google Scholar] [CrossRef]
  28. Rio, S.; Gallego-Sánchez, L.; Montilla-Bascón, G.; Canales, F.J.; Sánchez, J.I.Y.; Prats, E. Genomic prediction and training set optimization in a structured Mediterranean oat population. Theor. Appl. Genet. 2021, 134, 1–15. [Google Scholar] [CrossRef]
  29. Wu, J.; Wang, L.; Fu, J.; Chen, J.; Wei, S.; Zhang, S.; Zhang, J.; Tang, Y.; Chen, M.; Zhu, J.; et al. Resequencing of 683 common bean genotypes identifies yield component trait associations across a north-south cline. Nat. Genet. 2020, 52, 118–125. [Google Scholar] [CrossRef]
  30. Crossa, J.; Jarquín, D.; Franco, J.; Pérez-Rodríguez, P.; Burgueño, J.; Saint-Pierre, C.; Vikram, P.; Sansaloni, C.; Petroli, C.; Akdemir, D.; et al. Genomic Prediction of Gene Bank Wheat Landraces. G3-Genes Genomes Genet. 2016, 6, 1819–1834. [Google Scholar] [CrossRef] [Green Version]
  31. Kehel, Z.; Sanchez-Garcia, M.; El Baouchi, A.; Aberkane, H.; Tsivelikas, A.; Charles, C.; Amri, A. Predictive Characterization for Seed Morphometric Traits for Genebank Accessions Using Genomic Selection. Front. Ecol. Evol. 2020, 8, 32. [Google Scholar] [CrossRef] [Green Version]
  32. Gorjanc, G.; Jenko, J.; Hearne, S.J.; Hickey, J.M. Initiating maize pre-breeding programs using genomic selection to harness polygenic variation from landrace populations. BMC Genom. 2016, 17, 1–15. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Annicchiarico, P.; Nazzicari, N.; Ferrari, B.; Harzic, N.; Carroni, A.M.; Romani, M.; Pecetti, L. Genomic prediction of grain yield in contrasting environments for white lupin genetic resources. Mol. Breed. 2019, 39, 142. [Google Scholar] [CrossRef]
  34. Assefa, T.; Mahama, A.A.; Brown, A.V.; Cannon, E.K.S.; Rubyogo, J.C.; Rao, I.M.; Blair, M.W.; Cannon, S.B. A review of breeding objectives, genomic resources, and marker-assisted methods in common bean (Phaseolus vulgaris L.). Mol. Breed. 2019, 39, 20. [Google Scholar] [CrossRef] [Green Version]
  35. Kwak, M.; Gepts, P. Structure of genetic diversity in the two major gene pools of common bean (Phaseolus vulgaris L., Fabaceae). Theor. Appl. Genet. 2009, 118, 979–992. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Roos, A.D.; Hayes, B.J.; Goddard, M.E. Reliability of Genomic Predictions Across Multiple Populations. Genetics 2009, 183, 1545–1553. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Guo, Z.; Tucker, D.M.; Basten, C.J.; Gandhi, H.; Ersoz, E.; Guo, B.; Xu, Z.; Wang, D.; Gay, G. The impact of population structure on genomic prediction in stratified populations. Theor. Appl. Genet. 2014, 127, 749–762. [Google Scholar] [CrossRef]
  38. Norman, A.; Taylor, J.; Edwards, J.; Kuchel, H. Optimising Genomic Selection in Wheat: Effect of Marker Density, Population Size and Population Structure on Prediction Accuracy. G3-Genes Genomes Genet. 2018, 8, 2889–2899. [Google Scholar] [CrossRef] [Green Version]
  39. Sapkota, S.; Boyles, R.; Cooper, E.; Brenton, Z.; Myers, M.; Kresovich, S. Impact of sorghum racial structure and diversity on genomic prediction of grain yield components. Crop Sci. 2020, 60, 132–148. [Google Scholar] [CrossRef] [Green Version]
  40. Endelman, J.B. Ridge Regression and other Kernels for Genomic Selection with R Package rrBLUP. Plant Genome 2011, 4, 250–255. [Google Scholar] [CrossRef] [Green Version]
  41. Shikha, M.; Kanika, A.; Rao, A.R.; Mallikarjuna, M.G.; Gupta, H.S.; Nepolean, T. Genomic Selection for Drought Tolerance Using Genome-Wide SNPs in Maize. Front. Plant Sci. 2017, 8, 550. [Google Scholar] [CrossRef] [PubMed]
Plants 11 01298 g001 550
Figure 1. A workflow for genomic prediction of common bean germplasm. The above three ways are represented by the flow of the three numbers (1, 2, 3) in the figure. (1) The landraces from each of two subgene pools were used as training sets to predict breeding lines, respectively. (2) All landraces were used as training sets to predict breeding lines. (3) A part of the breeding lines was selected and included in the training sets.
Figure 1. A workflow for genomic prediction of common bean germplasm. The above three ways are represented by the flow of the three numbers (1, 2, 3) in the figure. (1) The landraces from each of two subgene pools were used as training sets to predict breeding lines, respectively. (2) All landraces were used as training sets to predict breeding lines. (3) A part of the breeding lines was selected and included in the training sets.
Plants 11 01298 g001
Plants 11 01298 g002 550
Figure 2. Different number of SNPs and different population sizes were used in ten-fold cross-validation using 100 replications for 15 agronomic traits in 484 common bean landraces (The figure shows four traits: DF, days to flowering; PH, plant height; SP, seeds per pod; and SH, seed height. Other traits are shown in Supplementary Figure S2).
Figure 2. Different number of SNPs and different population sizes were used in ten-fold cross-validation using 100 replications for 15 agronomic traits in 484 common bean landraces (The figure shows four traits: DF, days to flowering; PH, plant height; SP, seeds per pod; and SH, seed height. Other traits are shown in Supplementary Figure S2).
Plants 11 01298 g002
Plants 11 01298 g003 550
Figure 3. The three principal components of a PCA performed on 9781 SNPs markers among 628 common bean accessions. The four colors (shapes) represent landraces and breeding lines in the two gene pools.
Figure 3. The three principal components of a PCA performed on 9781 SNPs markers among 628 common bean accessions. The four colors (shapes) represent landraces and breeding lines in the two gene pools.
Plants 11 01298 g003
Table
Table 1. The predictive ability of landraces as the training set to predict breeding lines using 15 traits.
Table 1. The predictive ability of landraces as the training set to predict breeding lines using 15 traits.
TraitsLandraces (An + M)Landraces (An)Landraces (M)
Breeding Line (An + M)Breeding Line (An)Breeding Line (M)Breeding Line (An)Breeding Line (M)
DF0.32040.16990.20630.18320.2043
DM0.21150.3917−0.00450.4370−0.1070
PH0.67610.42010.65360.40350.6257
SNN0.64920.34920.38650.27470.2640
BN0.53650.45730.59070.44140.5954
PP0.77970.19740.68560.27310.6967
PL0.76460.56880.84260.57860.8441
PW0.67020.59270.68480.55740.6374
PDH0.72670.48960.52550.50880.5085
SP0.78180.24220.49950.26450.3562
GY0.59310.20960.47990.18100.4886
GW0.76160.33870.61060.33250.6046
SL0.77050.44580.70700.47440.7222
SW0.71950.46720.62540.40070.6146
SH0.81550.74010.69240.73300.6917
DF, days to flowering; DM, days to maturity; PH, plant height; SNN, stem node number; BN, branch number; PP, pod per plant; PL, pod length; PW, pod width; PDH, pod height; SP, seeds per pod; GY, grain yield per plant; GW, 100-grain weight; SL, seed length; SW, seed width; SH, seed height.
Table
Table 2. Two-step scenario for optimizing the training set by adding 50% of the breeding lines.
Table 2. Two-step scenario for optimizing the training set by adding 50% of the breeding lines.
ScenariosScenario 1Scenario 2
GroupsG1G2G3G4G5G6G7G8
Training setsLandraces An + M(484)Landraces An + M(484)Landraces An + M(484)Landraces An + M(484)Landraces An + M(484)Landraces An + M(484)Landraces An + M(484)Landraces An + M(484)
Additional
training sets		 Breeding lines A1 (30)Breeding lines A1 + M1 (72) Breeding lines M1 (42)Breeding lines A1 + M1 (72)Breeding lines A1 + M1 (72)Breeding lines random (72)
Testing setsBreeding lines A2 (30)Breeding lines A2 (30)Breeding lines A2 (30)Breeding lines M2 (42)Breeding lines M2 (42)Breeding lines M2 (42)Breeding lines A2 + M2 (72)Breeding lines remaining (72)
The 50% breeding lines (An) were randomly selected as breeding lines A1 (30). The remaining breeding lines (An) are used as breeding lines A2 (30). The 50% breeding lines (M) were randomly selected as breeding lines M1 (42). The remaining breeding lines (M) are used as breeding lines M2 (42). The 50% breeding lines were randomly selected as breeding lines random (72). The remaining breeding lines are used as breeding lines remaining (72).
Table
Table 3. The model’s precision for 8 groups across two scenarios among 15 traits.
Table 3. The model’s precision for 8 groups across two scenarios among 15 traits.
ScenariosScenario 1Scenario 2
GroupsG1G2G3G4G5G6G7G8
DF0.25660.35900.36250.21090.39340.38000.46640.3686
DM0.51080.64420.6237−0.08460.05470.03230.38680.4308
PH0.49580.62940.64460.64630.78180.77950.79040.8260
SNN0.39200.53990.55490.44970.58690.58650.72550.6850
BN0.55380.69500.68200.60530.69210.69960.67790.6922
PP0.24650.27710.26750.73800.73560.73980.79860.7385
PL0.67400.81030.81110.89340.93900.94560.89290.7724
PW0.65590.78400.78540.76620.89750.89660.88740.6407
PDH0.36260.37720.36310.33120.33550.33160.68370.5964
SP0.16110.31680.35040.53890.58550.58770.78640.7056
GY0.24680.24960.25150.57950.61910.62100.64400.5802
GW0.26130.52680.52740.62380.61650.61550.78040.6944
SL0.38760.46740.47750.73640.72850.72260.76470.6956
SW0.43140.66140.66800.70610.75960.76240.81370.7460
SH0.71140.79880.79880.66500.73890.72900.85350.8186
DF, days to flowering; DM, days to maturity; PH, plant height; SNN, stem node number; BN, branch number; PP, pod per plant; PL, pod length; PW, pod width; PDH, pod height; SP, seeds per pod; GY, grain yield per plant; GW, 100-grain weight; SL, seed length; SW, seed width; SH, seed height.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Shao, J.; Hao, Y.; Wang, L.; Xie, Y.; Zhang, H.; Bai, J.; Wu, J.; Fu, J. Development of a Model for Genomic Prediction of Multiple Traits in Common Bean Germplasm, Based on Population Structure. Plants 2022, 11, 1298. https://doi.org/10.3390/plants11101298

AMA Style

Shao J, Hao Y, Wang L, Xie Y, Zhang H, Bai J, Wu J, Fu J. Development of a Model for Genomic Prediction of Multiple Traits in Common Bean Germplasm, Based on Population Structure. Plants. 2022; 11(10):1298. https://doi.org/10.3390/plants11101298

Chicago/Turabian Style

Shao, Jing, Yangfan Hao, Lanfen Wang, Yuxin Xie, Hongwei Zhang, Jiangping Bai, Jing Wu, and Junjie Fu. 2022. "Development of a Model for Genomic Prediction of Multiple Traits in Common Bean Germplasm, Based on Population Structure" Plants 11, no. 10: 1298. https://doi.org/10.3390/plants11101298

APA Style

Shao, J., Hao, Y., Wang, L., Xie, Y., Zhang, H., Bai, J., Wu, J., & Fu, J. (2022). Development of a Model for Genomic Prediction of Multiple Traits in Common Bean Germplasm, Based on Population Structure. Plants, 11(10), 1298. https://doi.org/10.3390/plants11101298

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop