Evaluating Individual Genome Similarity with a Topic Model.

MOTIVATION Evaluating genome similarity among individuals is an essential step in data analysis. Advanced sequencing technology detects more and rarer variants for massive individual genomes, thus enabling individual-level genome similarity evaluation. However, the current methodologies such as the principal component analysis (PCA) lack the capability to fully leverage rare variants, and are also difficult to interpret in terms of population genetics. RESULTS Here we introduce a probabilistic topic model, latent Dirichlet allocation, to evaluate individual genome similarity. A total of 2,535 individuals from the 1,000 Genomes Project (KGP) were used to demonstrate our method. Various aspects of variant choice and model parameter selection were studied. We found that relatively rare (0.001 20,000 bp) variants are more efficient for genome similarity evaluation. At least 100,000 such variants are necessary. In our results, the populations show significantly less mixed and more cohesive visualization than the PCA results. The global similarities among the KGP genomes are consistent with known geographical, historical and cultural factors. AVAILABILITY The source code and data access are available at: https://github.com/lrjuan/LDA_genome. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.