A Pangenomic Approach to Improve Population Genetics Analysis and Reference Bias in Underrepresented Middle Eastern and Horn of Africa Populations.

Abstract

Genomics plays a crucial role in addressing health disparities, yet most studies rely on the hg38 linear reference genome, limiting the potential of pangenomic approaches, particularly for underrepresented populations. In this study, we focus on characterising East African populations, particularly Somalis, by constructing a variation graph using Mozabites from the Human Genome Diversity Project (HGDP) given their ancestral affinity with Somalis. We evaluated the effectiveness of this graph-based reference in estimating effective population sizes (Ne) in Bedouins compared to the hg38 reference and examined its impact on allele frequencies and genome-wide association studies (GWAS). Applying a coalescent model to the graph-based reference produced a Ne estimate of approximately 17 for the Bedouin population, which was significantly lower than the estimate from the hg38 reference (approximately 79,000). Only the graph-based estimate fell within the 95% confidence interval in simulations, indicating improved accuracy. Moreover, graph variants exhibited significantly lower allele frequencies (p-value < 2.2 × 10^-16), suggesting potential effects on the interpretation and power of GWAS. Notably, GWAS variants specific to Bedouins derived from the graph showed lower frequencies (p = 0.023) than those obtained from the linear reference. These findings suggest that a pangenomic approach, informed by populations with ancestral affinities such as the Mozabites, provides more accurate estimates of Ne and allele frequencies. This highlights the importance of pangenomic strategies to better capture genetic diversity in underrepresented populations, a critical step towards improving population genetics studies, personalised medicine, and equitable healthcare.