The ethnogeographic variability of genetic factors underlying G6PD deficiency

Glucose-6-phosphate dehydrogenase (G6PD) deficiency caused by genetic variants in the G6PD gene, constitutes the most common enzymopathy worldwide, affecting approximately 5% of the global population. While carriers are mostly asymptomatic, they are at substantial risk of acute hemolytic anemia upon certain infections or exposure to various medications. As such, information about G6PD activity status in a given patient can constitute an important parameter to guide clinical decision-making. Here, we leveraged whole genome sequencing data from 142,069 unrelated individuals across seven human populations to provide a global comprehensive map of G6PD variability. By integrating established functional classifications with stringent computational predictions using 13 partly orthogonal algorithms for uncharacterized and novel variants, we reveal the large extent of ethnogeographic variability in G6PD deficiency and highlight its population-specific genetic composition. Overall, estimated disease prevalence in males ranged between 12.2% in Africans, 2.7-3.5% across Asia and 2.1% in Middle Easterners to 202A/376 G (minor allele frequency 11.6%) and A-968 C/376 G (0.5%). In contrast, G6PD deficiency in Middle Easterners was primarily due to the Mediterranean allele (1.3%) and the population-specific Cairo variant (0.4%). In South Asia, the most prevalent deficient alleles were Mediterranean (1.7%), Kerala (1.1%), Gond (0.9%) and Orissa (0.2%), whereas in East Asian populations the Canton (1.1%), Kaiping (0.7%) and Viangchan (0.3%) variants were predominant. Combined, our analyses provide a large dataset of G6PD variability across major ethnogeographic groups and can instruct population-specific genotyping strategies to optimize genetically guided therapeutic interventions.