Robust Chi-square Test to Distinguish Population Stratification from Other Reasons of Hardy-Weinberg Disequilibrium
摘要
The Hardy-Weinberg Equilibrium (HWE) test is a standard component of quality control pipelines for genomic studies and databases. However, applying the HWE test to large, multi-ethnic datasets raises substantial challenges: natural departures from HWE due to population stratification and admixture can mask technical errors. We extended the classical chi-square goodness-of-fit test to account for the presence of population stratification. Assuming genotype datasets can be partitioned into relatively homogeneous strata, we changed the null hypothesis of the classical HWE test to assume that all strata are in HWE, and created a test statistic similar to the classical chi-square but applicable to stratified populations. This test statistic, together with a few others, can distinguish HWE deviations caused by population stratification from those resulting from other reasons. We have implemented this method in the software package PopHWE. Test results using simulated data and 1000 Genomes data showed that PopHWE outperforms the classical HWE meta-analysis using Stouffer's method, while being much faster, easier to implement, and more robust to small sample sizes, large numbers of subpopulations, and the presence of ancestry grouping errors. Using PopHWE, HWE deviations caused by population stratification and those caused by genotyping errors can be displayed in different clusters on scatterplots. Furthermore, we demonstrated that the combined z-value of meta-analysis is biased toward excess heterozygosity for genotype datasets with up to moderately large sample sizes, and showed that replacing the combined z-value with its squared value could significantly improve the performance.
参考文献
Abramovs, N., Brass, A., & Tassabehji, M. (2020). Hardy-Weinberg Equilibrium in the Large Scale Genomic Sequencing Era. Front Genet, 11, 210. https://doi.org/10.3389/fgene.2020.00210
Adams, M. W. D., Ewart, K. M., Carlile, N., Rose, H. A., Walker, J. A., Hutton, I., Ho, S. Y. W., & Lo, N. (2026). Pleistocene Sea-Level Fluctuation Shapes Archipelago-Wide Population Structure in the Critically Endangered Lord Howe Island Cockroach Panesthia lata. Ecol Evol, 16(3), e72760. https://doi.org/10.1002/ece3.72760
Atsawawaranunt, K., Stuart, K. C., Whibley, A., Ewart, K. M., Rollins, L. A., Johnson, R. N., & Santure, A. W. (2026). Missing or Mis-Telling the Story? Trade-Offs for Restriction-Site Associated Compared to Whole Genome Sequencing. Mol Ecol, 35(2), e70211. https://doi.org/10.1111/mec.70211
Attia J, Thakkinstian A, McElduff P, Milne E, Dawson S, Scott RJ, Klerk Nd, Armstrong B, Thompson J. Detecting genotyping error using measures of degree of Hardy-Weinberg disequilibrium. Stat Appl Genet Mol Biol. 2010;9:Article 5. doi: 10.2202/1544-6115.1463. Epub 2010 Jan 6. PMID: 20196755.
Bycroft, C., Freeman, C., Petkova, D., Band, G., Elliott, L. T., Sharp, K., Motyer, A., Vukcevic, D., Delaneau, O., O'Connell, J., Cortes, A., Welsh, S., Young, A., Effingham, M., McVean, G., Leslie, S., Allen, N., Donnelly, P., & Marchini, J. (2018). The UK Biobank resource with deep phenotyping and genomic data. Nature, 562(7726), 203-209. https://doi.org/10.1038/s41586-018-0579-z
Chang, C. C., Chow, C. C., Tellier, L. C., Vattikuti, S., Purcell, S. M., & Lee, J. J. (2015). Second-generation PLINK: rising to the challenge of larger and richer datasets. Gigascience, 4, 7. https://doi.org/10.1186/s13742-015-0047-8
Danecek, P., Auton, A., Abecasis, G., Albers, C. A., Banks, E., DePristo, M. A., Handsaker, R. E., Lunter, G., Marth, G. T., Sherry, S. T., McVean, G., Durbin, R., & Genomes Project Analysis, G. (2011). The variant call format and VCFtools. Bioinformatics, 27(15), 2156-2158. https://doi.org/10.1093/bioinformatics/btr330
Freeberg MA, Fromont LA, D'Altri T, et al. The European Genome-phenome Archive in 2021. Nucleic Acids Research. 2022;50(D1):D980-D987. doi:10.1093/nar/gkab1059
Faure, N., Vilcot, M., Pichot, F., Riutort, J. J., Barroil, A., Holon, F., Tomasi, N., Mouillot, D., Deter, J., & Manel, S. (2025). A Fragile Stronghold: Genomics Reveal Angelshark Population Vulnerability in Corsica, a Key Mediterranean Refuge. Ecol Evol, 15(10), e72275. https://doi.org/10.1002/ece3.72275
Fritsche LG, Igl W, Bailey JN, Grassmann F, Sengupta S, et al. (2016). A large genome-wide association study of age-related macular degeneration highlights contributions of rare and common variants. Nat Genet. 48:134–143.
Garg E, Romain J, Sun L, Paterson AD (2026) Assessing Hardy-Weinberg equilibrium in T2T-aligned 1000 genomes project. PLoS Genet 22(9): e1012280. https://doi.org/10.1371/journal.pgen.1012280
Gaziano, J. M., Concato, J., Brophy, M., Fiore, L., Pyarajan, S., Breeling, J., Whitbourne, S., Deen, J., Shannon, C., Humphries, D., Guarino, P., Aslan, M., Anderson, D., LaFleur, R., Hammond, T., Schaa, K., Moser, J., Huang, G., Muralidhar, S., . . . O'Leary, T. J. (2016). Million Veteran Program: A mega-biobank to study genetic influences on health and disease. J Clin Epidemiol, 70, 214-223. https://doi.org/10.1016/j.jclinepi.2015.09.016
Graffelman, J., & Weir, B. S. (2016). Testing for Hardy-Weinberg equilibrium at biallelic genetic markers on the X chromosome. Heredity (Edinb), 116(6), 558-568. https://doi.org/10.1038/hdy.2016.20
Grover, V.K., Cole, D.E.C. and Hamilton, D.C. (2010), Attributing Hardy-Weinberg Disequilibrium to Population Stratification and Genetic Association in Case-Control Studies. Annals of Human Genetics, 74: 77-87. https://doi.org/10.1111/j.1469-1809.2009.00552.x
Hao, W., Song, M., & Storey, J. D. (2016). Probabilistic models of genetic variation in structured populations applied to global human studies. Bioinformatics, 32(5), 713-721. https://doi.org/10.1093/bioinformatics/btv641
Hao, W., & Storey, J. D. (2019). Extending Tests of Hardy-Weinberg Equilibrium to Structured Populations. Genetics, 213(3), 759-770. https://doi.org/10.1534/genetics.119.302370
Hemstrom, W., Grummer, J. A., Luikart, G., & Christie, M. R. (2024). Next-generation data filtering in the genomics era. Nat Rev Genet, 25(11), 750-767. https://doi.org/10.1038/s41576-024-00738-6
Hernández J. L., Weir B. S. A disequilibrium coefficient approach to Hardy-Weinberg testing. Biometrics. 1989 Mar;45(1):53–70.
Hsu, Y. H. (2026). Hardy-Weinberg Equilibrium Filtering in Population Genomics: Empirical Review and Decision Framework for Improved Practice. Ecol Evol, 16(1), e72688. https://doi.org/10.1002/ece3.72688
Jin, Y., Wang, H., Naj, A. C., Wang, L. S., & Lee, W. P. (2026). GrafAnc: Reliable and reproducible inference of continental and regional population structure. HGG Adv, 7(1), 100530. https://doi.org/10.1016/j.xhgg.2025.100530
Kwong, A. M., Blackwell, T. W., LeFaive, J., de Andrade, M., Barnard, J., Barnes, K. C., Blangero, J., Boerwinkle, E., Burchard, E. G., Cade, B. E., Chasman, D. I., Chen, H., Conomos, M. P., Cupples, L. A., Ellinor, P. T., Eng, C., Gao, Y., Guo, X., Irvin, M. R., . . . Kang, H. M. (2021). Robust, flexible, and scalable tests for Hardy-Weinberg equilibrium across diverse ancestries. Genetics, 218(1). https://doi.org/10.1093/genetics/iyab044
Lindley D. (1998). Statistical inference concerning Hardy-Weinberg equilibrium. In: Bernardo JM BJ, Dawid AP, and Smith AFM, editor. Bayesian Statistics. 3 ed. Oxford: Oxford University Press. p 307-26.
Marees A.T., de Kluiver H., Stringer S., Vorspan F., Curis E., Marie-Claire C., Derks E.M. A tutorial on conducting genome-wide association studies: quality control and statistical analysis. Int. J. Methods Psychiatr. Res. 2018; 27:e1608.
Meisner, J., & Albrechtsen, A. (2019). Testing for Hardy-Weinberg equilibrium in structured populations using genotype or low-depth next generation sequencing data. Mol Ecol Resour, 19(5), 1144-1152. https://doi.org/10.1111/1755-0998.13019
Meisner, J., Albrechtsen, A., & Hanghoj, K. (2021). Detecting selection in low-coverage high-throughput sequencing data using principal component analysis. BMC Bioinformatics, 22(1), 470. https://doi.org/10.1186/s12859-021-04375-2
Mon, S. L. Y., Lwin, M., Maw, A. A., Htun, L. L., Bawm, S., Nagano, Y., Nagano, A. J., Kawabe, K., Wada, Y., Okamoto, S., & Shimogiri, T. (2025). Exploring genetic diversity and population structure of Myanmar indigenous chickens using double digest restriction site-associated DNA sequencing. Anim Genet, 56(4), e70038. https://doi.org/10.1111/age.70038
Neamatzadeh, H., Dastgheib, S. A., Mazaheri, M., Masoudi, A., Shiri, A., Omidi, A., Rahmani, A., Golshan-Tafti, A., Aghasipour, M., Yeganegi, M., Bahrami, M., Aghili, K., Khajehnoori, S., & Mosavi Jarrahi, A. (2024). Hardy-Weinberg Equilibrium in Meta-Analysis Studies and Large-Scale Genomic Sequencing Era. Asian Pac J Cancer Prev, 25(7), 2229-2235. https://doi.org/10.31557/APJCP.2024.25.7.2229
Novembre, J., & Peter, B. M. (2016). Recent advances in the study of fine-scale population structure in humans. Curr Opin Genet Dev, 41, 98-105. https://doi.org/10.1016/j.gde.2016.08.007
Pearman, W. S., Urban, L., & Alexander, A. (2022). Commonly used Hardy-Weinberg equilibrium filtering schemes impact population structure inferences using RADseq data. Mol Ecol Resour, 22(7), 2599-2613. https://doi.org/10.1111/1755-0998.13646
Purcell, S., Neale, B., Todd-Brown, K., Thomas, L., Ferreira, M. A., Bender, D., Maller, J., Sklar, P., de Bakker, P. I., Daly, M. J., & Sham, P. C. (2007). PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet, 81(3), 559-575. https://doi.org/10.1086/519795
Ramirez, A. H., Sulieman, L., Schlueter, D. J., Halvorson, A., Qian, J., Ratsimbazafy, F., Loperena, R., Mayo, K., Basford, M., Deflaux, N., Muthuraman, K. N., Natarajan, K., Kho, A., Xu, H., Wilkins, C., Anton-Culver, H., Boerwinkle, E., Cicek, M., Clark, C. R., . . . All of Us Research, P. (2022). The All of Us Research Program: Data quality, utility, and diversity. Patterns (N Y), 3(8), 100570. https://doi.org/10.1016/j.patter.2022.100570
Rothman K. (2002). Random error and the role of statistics. Epidemiology: An introduction. Oxford: Oxford University press.
Royo, J. L. (2020). Hardy Weinberg Equilibrium Disturbances in Case-Control Studies Lead to Non-Conclusive Results. Cell J, 22(4), 572-574. https://doi.org/10.22074/cellj.2021.7195
Schaid, D. J., Batzler, A. J., Jenkins, G. D., & Hildebrandt, M. A. (2006). Exact tests of Hardy-Weinberg equilibrium and homogeneity of disequilibrium across strata. Am J Hum Genet, 79(6), 1071-1080. https://doi.org/10.1086/510257
Sha, Q., & Zhang, S. (2011). A test of Hardy-Weinberg equilibrium in structured populations. Genet Epidemiol, 35(7), 671-678. https://doi.org/10.1002/gepi.20617
Shoemaker J, Painter I, Weir BS. A Bayesian characterization of Hardy-Weinberg disequilibrium. Genetics. 1998 Aug;149(4):2079-88. doi: 10.1093/genetics/149.4.2079
Shyr, D., Dey, R., Li, X., Zhou, H., Boerwinkle, E., Buyske, S., Daly, M., Gibbs, R. A., Hall, I., Matise, T., Reeves, C., Stitziel, N. O., Zody, M., Neale, B. M., & Lin, X. (2024). Semi-supervised machine learning method for predicting homogeneous ancestry groups to assess Hardy-Weinberg equilibrium in diverse whole-genome sequencing studies. Am J Hum Genet, 111(10), 2129-2138. https://doi.org/10.1016/j.ajhg.2024.08.018
Stouffer SA. 1949. The American Soldier. Princeton: Princeton University Press.
Troendle JF, Yu KF. A note on testing the Hardy-Weinberg law across strata. Ann Hum Genet. 1994 Oct;58(4):397-402. doi: 10.1111/j.1469-1809.1994.tb00735.x. PMID: 7864594
Wahlund, S. (1928). Composition of population and correlation phenomenon from the viewpoint of heredity. Hereditas, 11, 65-106.
Weir BS . (1996) Genetic Data Analysis II. Sinauer Associates: Massachusetts.
Wigginton, J. E., Cutler, D. J., & Abecasis, G. R. (2005). A note on exact tests of Hardy-Weinberg equilibrium. Am J Hum Genet, 76(5), 887-893. https://doi.org/10.1086/429864
Zeng, W., Ye, Z., Liu, X., Lin, H., & Wu, J. (2026). Genotyping-by-Sequencing Reveals Low Genetic Diversity and Pronounced Geographic Structuring in the Endangered Medicinal Plant Coptis chinensis var. brevisepala. Plants (Basel), 15(3). https://doi.org/10.3390/plants15030371
Zhou, W., Qin, H., Wang, H., Zhi, H., Huang, R., Hou, S., Wang, J., Chen, L., Tian, X., Qiao, Z., Mu, Z., & Diao, X. (2026). Genome-wide association analysis and candidate gene identification for plant height in Shanxi local foxtail millet varieties. Theor Appl Genet, 139(7). https://doi.org/10.1007/s00122-026-05303-9
指标
DOI：
Submission ID：
下载次数
已发布
如何引用
下载引用
利益冲突声明
作者声明无任何需要披露的利益冲突。
Copyright
本预印本的版权持有者为作者/资助方。
This work is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License.