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.

K-mer analysis of long-read alignment pileups for structural variant genotyping

A graph clustering algorithm for detection and genotyping of structural variants from long reads

HapKled: a Haplotype-Aware Structural Variant Calling Approach for Oxford Nanopore Sequencing Data

SVarp: pangenome-based structural variant discovery

Comprehensive Evaluation of Structural Variant Genotyping Methods Based on Long-Read Sequencing Data.

GGTyper: genotyping complex structural variants using short-read sequencing data

Benchmarking long-read aligners and SV callers for structural variation detection in Oxford nanopore sequencing data

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

Kled: an Ultra-Fast and Sensitive Structural Variant Detection Tool for Long-Read Sequencing Data.

Comprehensive and deep evaluation of structural variation detection pipelines with third-generation sequencing data

Tradeoffs in alignment and assembly-based methods for structural variant detection with long-read sequencing data

A Comparison of Structural Variant Calling from Short-Read and Nanopore-Based Whole-Genome Sequencing Using Optical Genome Mapping as a Benchmark

PanSVR: Pan-Genome Augmented Short Read Realignment for Sensitive Detection of Structural Variations

Sawfish: Improving long-read structural variant discovery and genotyping with local haplotype modeling

Long-read sequencing and structural variant characterization in 1,019 samples from the 1000 Genomes Project

FocalSV: target region-based structural variant assembly and refinement using single-molecule long read sequencing data

Personalizing pangenome graphs with k -mers

kngMap: Sensitive and Fast Mapping Algorithm for Noisy Long Reads Based on the K-Mer Neighborhood Graph

Seamless, rapid, and accurate analyses of outbreak genomic data using split k-mer analysis

STIX: Long-reads based Accurate Structural Variation Annotation at Population Scale

Haplotype-aware variant calling with PEPPER-Margin-DeepVariant enables high accuracy in nanopore long-reads