Abstract:Recent breakthroughs have enabled the inference of genealogies from large sequencing data-sets, accurately reconstructing local trees that describe genetic ancestry at each locus. These genealogies should also capture the correlation structure of local trees along the genome, reflecting historical recombination events and factors like demography and natural selection. However, whether reconstructed genealogies do accurately capture this correlation structure has not been rigorously explored. This is important to address, since uncovering regions that depart from expectations can drive the discovery of new biological phenomena. Addressing this is crucial, as uncovering regions that deviate from expectations can reveal new biological phenomena, such as the suppression of recombination allowing linked selection over broad regions, evidenced in humans and in adaptive introgression events in various species. We use a theoretical framework to characterise properties of genealogies, such as the distribution of genomic spans of clades and edges, and demonstrate that our theoretical results match observations in various simulated scenarios. Testing genealogies reconstructed using leading approaches, we find departures from theoretical expectations for all methods. However, for the method Relate, a set of simple corrections results in almost complete recovery of the target distributions. Applying these corrections to genealogies reconstructed using Relate for 2504 human genomes, we observe an excess of clades with unexpectedly long genomic spans (125 with p < 1 \cdot 10^{-12} clustering into 50 regions), indicating localised suppression of historical recombination. The strongest signal corresponds to a known inversion on chromosome 17, while the second strongest represents a previously unknown inversion on chromosome 10, which is most common (21%) in S.~Asians and correlates with GWAS hits for a range of phenotypes including immunological traits. Other signals suggest additional large inversions (4), copy number changes (2), and complex rearrangements or other variants (12), as well as 28 regions with strong support but no clear classification. Our approach can be readily applied to other species, and show that genealogies offer previously untapped potential to study structural variation and its impacts at a population level, revealing new phenomena impacting evolution.

Inverse population genetic problems with noise: inferring extent and structure of haplotype blocks from point allele frequencies

Statistical Properties of Nucleotide Clusters in DNA Sequences

Distribution of Recombination Crossovers and the Origin of Haplotype Blocks: the Interplay of Population History, Recombination, and Mutation

The Effect of Haplotype-Block Definitions on Inference of Haplotype-Block Structure and Htsnps Selection

Haplotype Block Partition with Limited Resources and Applications to Human Chromosome 21 Haplotype Data.

A Dynamic Programming Algorithm for Haplotype Block Partitioning and Its Application in Association Studies.

Dynamic programming algorithms for haplotype block partitioning: applications to human chromosome 21 haplotype data

Estimation of Population Parameters in Genetic Polymorphism Equilibrium by Overdominance Model

Inference of missing SNPs and information quantity measurements for haplotype blocks.

Randomly Distributed Crossovers May Generate Block-Like Patterns of Linkage Disequilibrium: an Act of Genetic Drift.

Dynamic Programming Algorithms for Haplotype Block Partitioning and Tag SNP Selection Using Haplotype Data or Genotype Data

Forward-backward algorithms with a biallelic mutation-drift model: Orthogonal polynomials, and a coalescent/urn-model based approach

Defining Loci in Restriction-Based Reduced Representation Genomic Data from Nonmodel Species: Sources of Bias and Diagnostics for Optimal Clustering

Genome-wide Compatible SNP Intervals and Their Properties

Fast and accurate haplotype inference with hidden markov model

The length of haplotype blocks and signals of structural variation in reconstructed genealogies

An Overview of the Haplotype Problems and Algorithms.

Integer programming framework for pangenome-based genome inference

Genotype Sequence Segmentation: Handling Constraints and Noise

Haplotype parsing: methods for extracting information from human genetic variations

High fraction of silent recombination in a finite population two-locus neutral birth-death-mutation model