Loci discovery, network-guided approach, and genomic prediction for drought tolerance index in a multi-parent advanced generation intercross (MAGIC) cowpea population

Abstract

Cowpea is a nutrient-dense legume that significantly contributes to the population’s diet in sub-Saharan Africa and other regions of the world. Improving cowpea cultivars to be more resilient to abiotic stress such as drought would be of great importance. The use of a multi-parent advanced generation intercross (MAGIC) population has been shown to be efficient in increasing the frequency of rare alleles that could be associated with important agricultural traits. In addition, drought tolerance index has been reported to be a reliable parameter for assessing crop tolerance to water-deficit conditions. Therefore, the objectives of this study were to evaluate the drought tolerance index for plant growth habit, plant maturity, flowering time, 100-seed weight, and grain yield in a MAGIC cowpea population, to conduct genome-wide association study (GWAS) and identify single nucleotide polymorphism (SNP) markers associated with the drought tolerance indices, to investigate the potential relationship existing between the significant loci associated with the drought tolerance indices, and to conduct genomic selection (GS). These analyses were performed using the existing phenotypic and genotypic data published for the MAGIC population which consisted of 305 F₈ recombinant inbred lines (RILs) developed at University of California, Riverside. The results indicated that: (1) large variation in drought tolerance indices existed among the cowpea genotypes, (2) a total of 14, 18, 5, 5, and 35 SNPs were associated with plant growth habit change due to drought stress, and drought tolerance indices for maturity, flowering time, 100-seed weight, and grain yield, respectively, (3) the network-guided approach revealed clear interactions between the loci associated with the drought tolerance traits, and (4) the GS accuracy varied from low to moderate. These results could be applied to improve drought tolerance in cowpea through marker-assisted selection (MAS) and genomic selection (GS). To the best of our knowledge, this is the first report on marker loci associated with drought tolerance indices in cowpea.

References

[1] Weng, Y. et al. A rapid method for measuring seed protein content in cowpea (Vigna unguiculata (L.) Walp). Am. J. Plant Sci. 8, 2387-2396 (2017).

[2] Perrino, P., Laghetti, G., Zeuli, P. L. S. & Monti, L. M. Diversification of cowpea in the mediterranean and other centres of cultivation. Genet. Resour. Crop. Evol. 40, 121-132 (1993).

[3] Frota, K. M. G., Soares, R. A. M. & Arêas, J. A. G. Chemical composition of cowpea (Vigna unguiculata L. Walp), BRS-Milênio cultivar. Food Sci. Technol. 28, 470-476 (2008).

[4] Moreira-Araújo, R. S. R., Sampaio, G. R., Soares, R. A. M., Silva, C. P. & Arêas, J. A. G. Identification and quantification of antioxidant compounds in cowpea. Rev. Cienc. Agron. 48, 799-805 (2017).

[5] Qin, J. et al. Population structure analysis and association mapping of seed antioxidant content in USDA cowpea (Vigna unguiculata L. Walp.) core collection using SNPs. Canad. J. Plant Sci. 96, 1026-1036 (2016).

[6] Wison, C., Liu, X., Lesch, S. M. & Suarez, D. L. Growth response of major U.S. cowpea cultivars. I. Biomass accumulation and salt tolerance. HortScience 41, 225-230 (2006).

[7] Singh, B. B., Ajeigbe, H. A., Tarawali, S. A., Fernandez-Rivera, S. & Abubakar, M. Improving the production and utilization of cowpea as food and fodder. Field Crops Res. 84, 169-177 (2006).

[8] Agbicodo, E. M., Fatokun, C. A., Muranaka, S., Visser, R. G. F. & Linden van der, C. G. Breeding drought tolerant cowpea: constraints, accomplishments, and future prospects. Euphytica 167, 353-370 (2009).

[9] Fatokun, C. A., Boukar, O. & Muranaka, S. Evaluation of cowpea (Vigna unguiculata (L.) Walp.) germplasm lines for tolerance to drought. Plant Genet. Resour. 10, 171-176 (2012).

[10] Burridge, J. D. et al. Genome-wide association mapping and agronomic impact of cowpea root architecture. Theor. Appl. Genet. 130, 419-431 (2017).

[11] Carvalho, M., Lino-Neto, T., Rosa, E. & Carnide, V. Cowpea: a legume crop for a challenging environment. J. Sci. Food Agric. 97, 4273-4284 (2017).

[12] Verbree, D. A., Singh, B. B. & Payne, W. A. Genetics and heritability of shoot drought tolerance in cowpea seedlings. Crop Sci. 55, 146-153 (2015).

[13] Timko, M. P. et al. Sequencing and analysis of the gene-rich space of cowpea. BMC Genomics 9, 103 (2008).

[14] Huang, B. E. et al. MAGIC populations in crops: current status and future prospects. Theor. Appl. Genet. 128, 999-1017 (2015).

[15] Ram, R., Mehta, M., Balmer, L., Gatti, D. M. & Morahan, G . Rapid identification of major-effect genes using the collaborative cross. Genetics 198, 75-86 (2014).

[16] Kover, P. X. et al. A multiparent advanced generation inter-cross to fine-map quantitative traits in Arabidopsis thaliana. PLoS Genet. 5, e1000551 (2009).

[17] Huang, B. E. et al. A multiparent advanced generation inter-cross population for genetic analysis in wheat. Plant Biotechnol. J. 10, 826-839 (2012).

[18] Bandillo, N. et al. Multi-parent advanced generation inter-cross (MAGIC) populations in rice: progress and potential for genetics research and breeding. Rice 6, 11 (2013).

[19] Gaur, P. M., Jukanti, A. K. & Varshney, R. K. Impact of genomic technologies on chickpea breeding strategies. Agronomy 2, 199-221 (2012).

[20] Huynh, B. L. et al. A multi-parent advanced generation inter-cross MAGIC population for genetic analysis and improvement of cowpea (Vigna Unguiculata L. Walp.). Plant J. 93, 1129-1142 (2018).

[21] Muñoz-Amatriaín, M. et al. Genome resources for climate-resilient cowpea, an essential crop for food security. Plant J. 89, 1042-1054 (2017).

[22] Olatoye, M. O., Hu, Z. & Aikpokpodion, P. O. Epistasis detection and modeling for genomic selection in cowpea (Vigna unguiculata L. Walp.). Front. Genet. 10, 677 (2019).

[23] Price, A. H. Believe it or not, QTLs are accurate! Trends Plant Sci. 11, 213-216 (2006).

[24] Hamblin, M. T., Buckler, E. S. & Jannink, J. L. Population genetics of genomics-based crop improvement methods. Trends Genet. 27, 98-106 (2011).

[25] Ravelombola, W. S. et al. Investigation on various aboveground traits to identify drought tolerance in cowpea seedlings. HortScience 53, 1757-1765 (2018).

[26] Saad, F. F., El-Mohsen, A. A. A., Abd, M. A. & Al-Soudan, I. H. Effective selection criteria for evaluating some barley crosses for water stress tolerance. Adv. Agric. Biol. 1, 112-123 (2014).

[27] Hayes, B. J., Bowman, P. J., Chamberlain, A. C., Verbyla, K. & Goddard, M. E. Accuracy of genomic breeding values in multi-breed dairy cattle populations. Genet. Sel. Evol. 41, 51 (2009).

[28] Heffner, E. L., Sorrells, M. E. & Jannink, J. L. Genomic selection for crop improvement. Crop Sci. 49, 1-12 (2009).

[29] Muchero, W., Ehlers, J. D., Close, T. J., Roberts, P. A. & Mapping, Q. T. L. for drought stress-induced premature senescence and maturity in cowpea [Vigna unguiculata (L.) Walp.]. Theor. Appl. Genet. 118, 849-863 (2009).

[30] Ehlers, J. D., Hall, A. E., Patel, P. N., Roberts, P. A. & Matthews, W. C. Registration of `California Blackeye 27’ cowpea. Crop Sci. 40, 854 (2000).

[31] Huynh, B. L. et al. A major QTL corresponding to the Rk locus for resistance to root-knot nematodes in cowpea (Vigna unguiculata L. Walp.). Theor. Appl. Genet. 129, 87-95 (2016).

[32] Lucas, M. R., Ehlers, J. D., Roberts, P. A. & Timothy, J. C. Markers for quantitative inheritance of resistance to foliar thrips in cowpea. Crop Sci. 52, 2075-2081 (2012).

[33] Muchero, W., Ehlers, J. D., Close, T. J., Roberts, P. A. & Genic, S. N. P. markers and legume synteny reveal candidate genes underlying QTL for Macrophomina phaseolina resistance and maturity in cowpea [Vigna unguiculata (L) Walp.]. BMC Genomics 12, 8 (2011).

[34] Ouédraogo, J. T., Tignegre, J. B., Timko, M. P. & Belzile, F. J. AFLP markers linked to resistance against Striga gesnerioides race 1 in Cowpea (Vigna unguiculata). Genome 45, 787-793 (2002).

[35] Pottorff, M. O., Li, G., Ehlers, J. D., Close, T. J. & Roberts, P. A. Genetic mapping, synteny, and physical location of two loci for Fusarium oxysporum f. sp. tracheiphilum race 4 resistance in cowpea [Vigna unguiculata (L.) Walp]. Mol. Breed. 33, 779-791 (2014).

[36] Team R Development Core . R: a language and environment for statistical computing (R foundation for statistical computing, Vienna, 2011).

[37] Huang, M., Liu, X., Zhou, Y., Summers, R. M. & Zhang, Z. BLINK: a package for the next level of genome-wide association studies with both individuals and markers in the millions. GigaScience 8, giy154 (2019).

[38] Kaler, A. S. et al. Genome-wide association mapping of carbon isotope and oxygen isotope ratios in diverse soybean genotypes. Crop Sci. 57, 3085-3100 (2017).

[39] Shin, J. H. et al. Package ‘Ldheatmap’. https://sfustatgen.github.io/LDheatmap/ (2019).

[40] Fang, C. et al. Genome-wide association studies dissect the genetic networks underlying agronomical traits in soybean. Genome Biol. 18, 161 (2017).

[41] Purcell, S. et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 81, 559-575 (2007).

[42] Otasek, D., Morris, J. H., Bouças, J., Pico, A. R. & Demchak, B. Cytoscape automation: empowering workflow-based network analysis. Genome Biol. 20, 1-15 (2019).

[43] Meuwissen, T. H. E., Hayes, B. J. & Goddard, M. E. Prediction of total genetic value using genome-wide dense marker maps. Genetics 257, 1819-1829 (2001).

[44] Endelman, J. B. Ridge regression and other kernels for genomic selection with R package rrBLUP. Plant Genome J. 4, 250-255 (2011).

[45] Shikha, M. et al. genomic selection for drought tolerance using genome-wide SNPs in maize. Front. Plant Sci. 8, 550 (2017).

[46] Singh, B. B., Mai-Kodomi, Y. & Terao, T. A simple screening method for drought tolerance in cowpea. Indian J. Genet. Plant Breed. 59, 211-220 (1999).

[47] Assefa, T. et al. Genome-Wide associations and epistatic interactions for internode number, plant height, seed weight and seed yield in soybean. BMC Genomics 20, 527 (2019).

[48] Diers, B. W. et al. Genetic Architecture of soybean yield and agronomic traits. G3 8, 3367-3375 (2018).

[49] Jarzyniak, K. M. & Jasiński, M. Membrane transporters and drought resistance-a complex issue. Front. Plant Sci. 5, 687 (2014).

[50] Butt, H. I. et al. GaMYB85, an R2R3 MYB Gene, in transgenic Arabidopsis plays an important role in drought tolerance. BMC Plant Biol. 17, 142 (2017).

[51] Tang, Y. et al. Overexpression of a MYB family gene, OsMYB6, Increases drought and salinity stress tolerance in transgenic rice. Front. Plant Sci. 10, 168 (2019).

[52] Stracke, R., Werber, M. & Weisshaar, B. The R2R3-MYB gene family in Arabidopsis thaliana. Curr. Opin. Plant Biol. 4, 447-456 (2001).

[53] Ravelombola, W. S. et al. Genome-wide association study and genomic selection for soybean chlorophyll content associated with soybean cyst nematode tolerance. BMC Genomics 20, 904 (2019).

[54] Lozada, D. N., Mason, R. E., Sarinelli, J. M. & Brown-Guedira, G. Accuracy of genomic selection for grain yield and agronomic traits in soft red winter wheat. BMC Genet. 20, 82 (2019).