Abstract:Long-read sequencing (LRS) enables variant calling of high-quality structural variants (SVs). Genotypers of SVs utilize these precise call sets to increase the recall and precision of genotyping in short-read sequencing (SRS) samples. With the extensive growth in availabilty of SRS datasets in recent years, we should be able to calculate accurate population allele frequencies of SV. However, reprocessing hundreds of terabytes of raw SRS data to genotype new variants is impractical for population-scale studies, a computational challenge known as the N+1 problem. Solving this computational bottleneck is necessary to analyze new SVs from the growing number of pangenomes in many species, public genomic databases, and pathogenic variant discovery studies. To address the N+1 problem, we propose The Great Genotyper, a population genotyping workflow. Applied to a human dataset, the workflow begins by preprocessing 4.2K short-read samples of a total of 183TB raw data to create an 867GB Counting Colored De Bruijn Graph (CCDG). The Great Genotyper uses this CCDG to genotype a list of phased or unphased variants, leveraging the CCDG population information to increase both precision and recall. The Great Genotyper offers the same accuracy as the state-of-the-art genotypers with the addition of unprecedented performance. It took 100 hours to genotype 4.5M variants in the 4.2K samples using one server with 32 cores and 145GB of memory. A similar task would take months or even years using single-sample genotypers. The Great Genotyper opens the door to new ways to study SVs. We demonstrate its application in finding pathogenic variants by calculating accurate allele frequency for novel SVs. Also, a premade index is used to create a 4K reference panel by genotyping variants from the Human Pangenome Reference Consortium (HPRC). The new reference panel allows for SV imputation from genotyping microarrays. Moreover, we genotype the GWAS catalog and merge its variants with the 4K reference panel. We show 6.2K events of high linkage between the HPRC's SVs and nearby GWAS SNPs, which can help in interpreting the effect of these SVs on gene functions. This analysis uncovers the detailed haplotype structure of the human fibrinogen locus and revives the pathogenic association of a 28 bp insertion in the FGA gene with thromboembolic disorders.

The power and promise of population genomics: from genotyping to genome typing

An overview of current population genomics methods for the analysis of whole‐genome resequencing data in eukaryotes

Large-scale Genotyping of Complex DNA

The Great Genotyper: A Graph-Based Method for Population Genotyping of Small and Structural Variants

Genotyping-by-Sequencing for Populus Population Genomics: An Assessment of Genome Sampling Patterns and Filtering Approaches

Pangenome-based genome inference allows efficient and accurate genotyping across a wide spectrum of variant classes

Population genomics for wildlife conservation and management

Genotype Calling from Population-Genomic Sequencing Data

Combining population genomics and quantitative genetics: finding the genes underlying ecologically important traits

Genome-wide genetic marker discovery and genotyping using next-generation sequencing

The future is now: Amplicon sequencing and sequence capture usher in the conservation genomics era

Population Diversity at the Single-Cell Level

Genotyping Polyploids from Messy Sequencing Data

Typological thinking in human genomics research contributes to the production and prominence of scientific racism

A Fast and Flexible Statistical Model for Large-Scale Population Genotype Data: Applications to Inferring Missing Genotypes and Haplotypic Phase

Benchmarking Imputed Low Coverage Genomes in a Human Population Genetics Context

Population Genomics Along With Quantitative Genetics Provides a More Efficient Valorization of Crop Plant Genetic Diversity in Breeding and Pre-breeding Programs

Evaluation of whole exome sequencing as an alternative to BeadChip and whole genome sequencing in human population genetic analysis

Optimization Methods for Genotype Data Analysis in Epidemiological Studies

Short Tandem Repeats in the era of next-generation sequencing: from historical loci to population databases

Impact of genotyping errors on statistical power of association tests in genomic analyses: A case study.