A Practical Framework for Constructing Population-Specific and Alternate-Contig-Aware Genome References: A case study of Vietnam

Current studies in human genomics typically rely on the standard genome reference GRCh38 which is known to be biased toward populations of European ancestry and therefore has limitations when applied to other populations. Although various graph-based pangenome references were constructed for several populations to deal with this bias, their usage in practice is currently still limited compared to linear genome references. Here we present a framework for constructing a population-specific genome reference using GRCh38 as backbone with alternate-contig awareness to enhance genomic data analysis in the target population. We demonstrated the advantages of our framework using both public and in-house Vietnamese whole-genome sequencing (WGS) datasets. Genomic variants derived from high-coverage WGS data of the 1000 Vietnamese Genomes Project (VN1K) were imported into our framework to build a Vietnamese-specific Genome Reference (VGR). VGR was then compared to GRCh38 in read alignment and variant calling using high-coverage WGS data of 99 Vietnamese individuals (KHV) from the 1000 Genomes Project (1kGP). Using Omni array genotyping data from 99 KHV samples as an independent benchmark, we found that VGR improved variant-calling precision and reduced false-positive calls compared to GRCh38. Our framework could be easily used for other populations as long as they have a variant database similar to VN1K. Our code is publicly available at github.com/VinGenome/VGR