Asymptotic tests for Hardy–Weinberg equilibrium in hexaploids

Abstract

Hexaploids, a group of organisms containing three complete sets of chromosomes in a single nucleus, are of utmost importance to evolutionary studies and breeding programs. Many studies have focused on hexaploid linkage analysis and QTL mapping in controlled crosses, but little methodology has been developed to reveal how hexaploids diversify and evolve in natural populations. We formulate a general framework for studying the pattern of genetic variation in autohexaploid populations through testing deviation from Hardy– Weinberg equilibrium (HWE) at individual molecular markers. We confirm that hexaploids cannot reach exact HWE but can approach asymptotic HWE at 8–9 generations of random mating. We derive a statistical algorithm for testing HWE and the occurrence of double reduction for autopolyploids, a phenomenon that affects population variation during long evolutionary processes. We perform computer simulation to validate the statistical behavior of our test procedure and demonstrate its usefulness by analyzing a real data set for autohexaploid chrysanthemum. When extended to allohexaploids, our test procedure will provide a generic tool for illustrating the genome structure of hexaploids in the quest to infer their evolutionary status and design association studies of complex traits.

References

[1] Bever JD, Felber F . The theoretical population genetics of autopolyploidy. Oxf Surv Evol Biol. 1992; 8: 185-217.

[2] Otto SP, Whitton J . Polyploid incidence and evolution. Annu Rev Genet. 2000; 34: 401-37.

[3] Soltis DE, Soltis PS, Tate JA . Advances in the study of polyploidy since plant speciation. New Phytol. 2004; 161: 173-91.

[4] Van de Peer Y, Ashman TL, Soltis PS et al. Polyploidy: an evolutionary and ecological force in stressful times. Plant Cell. 2021; 33: 11-26.

[5] Periyasamy S, Schwarz E, Popejoy AB et al. Genome-wide association studies in ancestrally diverse populations: opportunities, methods, pitfalls, and recommendations. Cell. 2019; 179: 589-603.

[6] Hardy GH . Mendelian proportions in a mixed population. Science. 1908; 28: 49-50.

[7] Boyer SH . Classic papers on genetics. (book reviews: papers on human genetics). Science. 1963; 142: 1646.

[8] Stern C . The Hardy-Weinberg law. Science. 1943; 97: 137-8.

[9] Crow JF . Eighty years ago: the beginnings of population genetics. Genetics. 1988; 119: 473-6.

[10] Mayo O . A century of Hardy-Weinberg equilibrium. Twin Res Hum Genet. 2008; 11: 249-56.

[11] Engels WR . Exact tests for Hardy-Weinberg proportions. Genetics. 2009; 183: 1431-41.

[12] Wang J, Shete S . Testing departure from Hardy-Weinberg proportions. Methods Mol Biol. 2012; 850: 77-102.

[13] Leal SM . Detection of genotyping errors and pseudo-SNPs via deviations from Hardy-Weinberg equilibrium. Genet Epidemiol. 2005; 29: 204-14.

[14] Hosking L, Lumsden S, Lewis K et al. Detection of genotyping errors by Hardy-Weinberg equilibrium testing. Eur J Hum Genet. 2004; 12: 395-9.

[15] Dufresne F, Stift M, Vergilino R et al. Recent progress and challenges in population genetics of polyploid organisms: an overview of current state-of-the-art molecular and statistical tools. Mol Ecol. 2014; 23: 40-69.

[16] Meirmans PG, Liu S, Tienderen PV . The analysis of polyploid genetic data. J Hered. 2018; 109: 283-96.

[17] Sun L, Gan J, Jiang L et al. Recursive test of Hardy-Weinberg equilibrium in tetraploids. Trends Genet. 2021; 37: 504-13.

[18] Geiringer H. Chromatid segregation of tetraploids and hexaploids. Genetics. 1949; 34: 665-84.

[19] Darlington CD . Chromosome behaviour and structural hybridity in the tradescantiae. Genetics. 1929; 21: 207-86.

[20] Mather K. Segregation and linkage in autotetraploids. J Genet. 1936; 32: 287-314.

[21] Templeton AR . Population Genetics and Microevolutionary Theory . Oxford, UK: John Wiley & Sons, 2006.

[22] Sumitomo K, Shirasawa K, Isobe S et al. Genome-wide association study overcomes the genome complexity in autohexaploid chrysanthemum and tags SNP markers onto the flower color genes. Sci Rep. 2019; 9: 1-9.

[23] Wu RL, Gallo-Meagher M, Littell RC et al. A general polyploid model for analyzing gene segregation in outcrossing tetraploid species. Genetics. 2001; 159: 869-82.

[24] Ma CX, Casella G, Shen ZJ et al. A unified framework for mapping quantitative trait loci in bivalent tetraploids using single-dose restriction fragments: a case study from alfalfa. Genome Res. 2002; 12: 1974-81.

[25] Wu RL, Ma CX . A general framework for statistical linkage analysis in multivalent tetraploids. Genetics. 2005; 170: 899-907.

[26] Huang K, Dunn DW, Li Z et al. Inference of individual ploidy level using codominant markers. Mol Ecol Resour. 2019; 19: 1218-29.

[27] Huang K, Dunn DW, Ritland K et al. POLYGENE: population genetics analyses for autopolyploids based on allelic phenotypes. Methods Ecol Evol. 2020; 11: 448-56.

[28] Yang D, Li F, Wang J et al. A framework to model a web of linkage disequilibria for natural allotetraploid populations. Methods Ecol Evol. 2021; 00: 1-9.

[29] Butruille DV, Boiteux LS . Selection-mutation balance in polysomic tetraploids: impact of double reduction and gametophytic selection on the frequency and subchromosomal localization of deleterious mutations. Proc Natl Acad Sci U S A. 2000; 97: 6608-13.

[30] Lu Y, Yang X, Tong C et al. A multivalent three-point linkage analysis model of autotetraploids. Brief Bioinform. 2013; 14: 460-8.

[31] Jiang LB, Ren XY, Wu RL . Computational characterization of double reduction in autotetraploid natural populations. Plant J. 2021; 105: 1703-9.

[32] Okada Y, Monden Y, Nokihara K et al. Genome-wide association studies (GWAS) for yield and weevil resistance in sweet potato (Ipomoea batatas (L.) lam) . Plant Cell Rep. 2019; 38: 1383-92.

[33] Bararyenya A, Olukolu BA, Tukamuhabwa P et al. Genome-wide association study identified candidate genes controlling continuous storage root formation and bulking in hexaploid sweetpotato. BMC Plant Biol. 2020; 20: 3.

[34] Kumar D, Sharma S, Sharma R et al. Genome-wide association study in hexaploid wheat identifies novel genomic regions associated with resistance to root lesion nematode (Pratylenchus thornei) . Sci Rep. 2021; 11: 3572.