Genomic selection with fixed-effect markers improves the prediction accuracy for Capsaicinoid contents in Capsicum annuum

Abstract

Capsaicinoids provide chili peppers (Capsicum spp.) with their characteristic pungency. Several structural and transcription factor genes are known to control capsaicinoid contents in pepper. However, many other genes also regulating capsaicinoid contents remain unknown, making it difficult to develop pepper cultivars with different levels of capsaicinoids. Genomic selection (GS) uses genome-wide random markers (including many in undiscovered genes) for a trait to improve selection efficiency. In this study, we predicted the capsaicinoid contents of pepper breeding lines using several GS models trained with genotypic and phenotypic data from a training population. We used a core collection of 351 Capsicum accessions and 96 breeding lines as training and testing populations, respectively. To obtain the optimal number of single nucleotide polymorphism (SNP) markers for GS, we tested various numbers of genome-wide SNP markers based on linkage disequilibrium. We obtained the highest mean prediction accuracy (0.550) for different models using 3294 SNP markers. Using this marker set, we conducted GWAS and selected 25 markers that were associated with capsaicinoid biosynthesis genes and quantitative trait loci for capsaicinoid contents. Finally, to develop more accurate prediction models, we obtained SNP markers from GWAS as fixed-effect markers for GS, where 3294 genome-wide SNPs were employed. When four to five fixed-effect markers from GWAS were used as fixed effects, the RKHS and RR-BLUP models showed accuracies of 0.696 and 0.689, respectively. Our results lay the foundation for developing pepper cultivars with various capsaicinoid levels using GS for capsaicinoid contents.

References

[1] Fraenkel L, Bogardus ST, Concato J et al. Treatment options in knee osteoarthritis. Arch Intern Med. 2004; 164: 1299.

[2] Yang ZH, Wang XH, Wang HP et al. Capsaicin mediates cell death in bladder cancer T24 cells through reactive oxygen species production and mitochondrial depolarization. Urology. 2010; 75: 735-41.

[3] Urbina SL, Roberts MD, Kephart WC et al. Effects of twelve weeks of capsaicinoid supplementation on body composition, appetite and self-reported caloric intake in overweight individuals. Appetite. 2017; 113: 264-73.

[4] Tanaka Y, Hosokawa M, Miwa T et al. Newly mutated putative-aminotransferase in nonpungent pepper (Capsicum annuum) results in biosynthesis of capsinoids, capsaicinoid analogues . J Agric Food Chem. 2010; 58: 1761-7.

[5] Stewart C, Kang BC, Liu K et al. The Pun1 gene for pungency in pepper encodes a putative acyltransferase. Plant J. 2005; 42: 675-88.

[6] Mazourek M, Pujar A, Borovsky Y et al. A dynamic Interface for Capsaicinoid systems biology. Plant Physiol. 2009; 150: 1806-21.

[7] Han K, Jang S, Lee JH et al. A MYB transcription factor is a candidate to control pungency in Capsicum annuum . Theor Appl Genet. 2019; 132: 1235-46.

[8] Sun B, Zhou X, Chen C et al. Coexpression network analysis reveals an MYB transcriptional activator involved in capsaicinoid biosynthesis in hot peppers. Hortic Res. 2020; 7: 162.

[9] Zhu Z, Sun B, Cai W et al. Natural variations in the MYB transcription factor MYB31 determine the evolution of extremely pungent peppers. New Phytol. 2019; 223: 922-38.

[10] Medina-Lara F, Echevarría-Machado I, Pacheco-Arjona R et al. Influence of nitrogen and potassium fertilization on fruiting and capsaicin content in habanero pepper (Capsicum chinense Jacq.) . HortScience. 2008; 43: 1549-54.

[11] Kim JS, Park M, Lee DJ et al. Characterization of putative capsaicin synthase promoter activity. Molecules and Cells. 2009; 28: 331-9.

[12] Arce-Rodríguez ML, Ochoa-Alejo N . An R2R3-MYB transcription factor regulates Capsaicinoid biosynthesis. Plant Physiol. 2017; 174: 1359-70.

[13] Han K, Lee HHY, Ro NNY et al. QTL mapping and GWAS reveal candidate genes controlling capsaicinoid content in Capsicum. Plant Biotechnol J. 2018; 16: 1546-58.

[14] Goddard M . Genomic selection: prediction of accuracy and maximisation of long term response. Genetica. 2009; 136: 245-57.

[15] Meuwissen THE, Hayes BJ, Goddard ME . Prediction of Total genetic value using genome-wide dense marker maps. Genetics. 2001; 157: 1819-29.

[16] Schaeffer LRR . Strategy for applying genome-wide selection in dairy cattle. J Anim Breed Genet. 2006; 123: 218-23.

[17] Bernardo R, Yu J . Prospects for Genomewide selection for quantitative traits in maize. Crop Sci. 2007; 47: 1082-90.

[18] Zhao Y, Mette MF, Reif JC . Genomic selection in hybrid breeding. Plant Breed. 2015; 134: 1-10.

[19] Bernardo R . Genomewide selection when major genes are known. Crop Sci. 2014; 54: 68-75.

[20] Sarinelli JM, Murphy JP, Tyagi P et al. Training population selection and use of fixed effects to optimize genomic predictions in a historical USA winter wheat panel. Theor Appl Genet. 2019; 132: 1247-61.

[21] Spindel JE, Begum H, Akdemir D et al. Genome-wide prediction models that incorporate de novo GWAS are a powerful new tool for tropical rice improvement. Heredity (Edinb). 2016; 116: 395-408.

[22] Hong JP, Ro N, Lee HY et al. Genomic selection for prediction of fruit-related traits in pepper (Capsicum spp.). Frontiers in Plant Sci. 2020; 11: 1-11, 570871.

[23] Ben-Chaim A, Borovsky Y, Falise M et al. QTL analysis for capsaicinoid content in Capsicum. Theor Appl Genet. 2006; 113: 1481-90.

[24] Yarnes SC, Ashrafi H, Reyes-Chin-Wo S et al. Identification of QTLs for capsaicinoids, fruit quality, and plant architecture-related traits in an interspecific Capsicum RIL population. Genome. 2013; 56: 61-74.

[25] Lee J, Park SJ, Hong SC et al. QTL mapping for capsaicin and dihydrocapsaicin content in a population of Capsicum annuum ‘NB1’ × Capsicum chinense ‘Bhut Jolokia.’ . Plant Breed. 2016; 135: 376-83.

[26] Jarret RL, Dang P . Revisiting the waxy locus and the CAPSICUM ANNUUM l. COMPLEX. Georgia Journal of Science. 2004; 62: 118-32.

[27] Jarret RL . DNA barcoding in a crop Genebank: the Capsicum annuum species complex. The Open Biology Journal. 2008; 1: 35-42.

[28] Wang J, Zhang Z . GAPIT version 3: boosting power and accuracy for genomic association and prediction. Genomics, Proteomics & Bioinformatics. 2021; 19: 629-40.

[29] Yu J, Pressoir G, Briggs WH et al. A unified mixed-model method for association mapping that accounts for multiple levels of relatedness. Nat Genet. 2006; 38: 203-8.

[30] Li M, Liu X, Bradbury P et al. Enrichment of statistical power for genome-wide association studies. BMC Biol. 2014; 12: 73.

[31] Wang H, Ye M, Fu Y et al. Modeling genome-wide by environment interactions through omnigenic interactome networks. Cell Rep. 2021; 35: 1-13, 109114.

[32] Endelman JB . Ridge regression and other kernels for genomic selection with R package rrBLUP. The Plant Genome. 2011; 4: 250-5.

[33] Montesinos López OA, Montesinos López A, Crossa J . Random Forest for Genomic Prediction. In: Multivariate Statistical Machine Learning Methods for Genomic Prediction.Cham: Springer, 2022, 633-81.

[34] Zou H, Hastie T . Regularization and variable selection via the elastic net. Journal of the Royal Statistical Society: Series B (Statistical Methodology). 2005; 67: 301-20.

[35] Friedman J, Hastie T, Tibshirani R . Regularization paths for generalized linear models via coordinate descent. J Stat Softw. 2010; 33: 1-22.

[36] Ogutu JO, Schulz-Streeck T, Piepho HP . Genomic selection using regularized linear regression models: ridge regression, lasso, elastic net and their extensions. BMC Proc. 2012; 6: S10.

[37] Habier D, Fernando RL, Kizilkaya K et al. Extension of the bayesian alphabet for genomic selection. BMC Bioinformatics. 2011; 12: 186.

[38] Park T, Casella G . The Bayesian lasso. J Am Stat Assoc. 2008; 103: 681-6.

[39] Mutshinda CM, Sillanpää MJ . Extended Bayesian LASSO for multiple quantitative trait loci mapping and unobserved phenotype prediction. Genetics. 2010; 186: 1067-75.

[40] Stewart-Brown BB, Song Q, Vaughn JN et al. Genomic selection for yield and seed composition traits within an applied soybean breeding program. G3 Genes|Genomes|Genetics. 2019; 9: 2253-65.

[41] An Y, Chen L, Li YX et al. Genome-wide association studies and whole-genome prediction reveal the genetic architecture of KRN in maize. BMC Plant Biol. 2020; 20: 490.

[42] Ravelombola WS, Qin J, Shi A et al. Genome-wide association study and genomic selection for soybean chlorophyll content associated with soybean cyst nematode tolerance. BMC Genomics. 2019; 20: 904.

[43] Garcia AAF, Wang S, Melchinger AE et al. Quantitative trait loci mapping and the genetic basis of Heterosis in maize and Rice. Genetics. 2008; 180: 1707-24.

[44] Chan DI, Vogel HJ . Current understanding of fatty acid biosynthesis and the acyl carrier protein. Biochem J. 2010; 430: 1-19.

[45] Lee HY, Ro NY, Jeong HJ et al. Genetic diversity and population structure analysis to construct a core collection from a large Capsicum germplasm . BMC Genet. 2016; 17: 142.

[46] Wimmer V, Albrecht T, Auinger HJ et al. Synbreed: a framework for the analysis of genomic prediction data using R. Bioinformatics. 2012; 28: 2086-7.

[47] Han K, Jeong HJ, Sung J et al. Biosynthesis of capsinoid is controlled by the Pun1 locus in pepper. Mol Breed. 2013; 31: 537-48.

[48] Bates D, Mächler M, Bolker B et al. Fitting linear mixed-effects models using lme4. J Stat Softw. 2015; 67: 1-48.

[49] Breiman L . Random forests. Mach Learn. 2001; 45: 5-32.

[50] Gianola D, van Kaam JBCHM . Reproducing kernel Hilbert spaces regression methods for genomic assisted prediction of quantitative traits. Genetics. 2008; 178: 2289-303.