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.

A spectral graph approach to discovering genetic ancestry

Novel probabilistic models of spatial genetic ancestry with applications to stratification correction in genome-wide association studies

European American Stratification in Ovarian Cancer Case Control Data: the Utility of Genome-Wide Data for Inferring Ancestry

Principal Component Analyses in Anthropological Genetics

A geographic history of human genetic ancestry

Ancestry inference using principal component analysis and spatial analysis: a distance-based analysis to account for population substructure

GemTools: A fast and efficient approach to estimating genetic ancestry

GRAF-pop: A Fast Distance-Based Method To Infer Subject Ancestry from Multiple Genotype Datasets Without Principal Components Analysis

Power Analysis of Principal Components Regression in Genetic Association Studies.

Analysis of East Asia genetic substructure using genome-wide SNP arrays

SNVstory: inferring genetic ancestry from genome sequencing data

Detecting Genomic Signatures of Natural Selection with Principal Component Analysis: Application to the 1000 Genomes Data

Manifold Learning for Human Population Structure Studies

Addressing Ancestry Disparities in Genomic Medicine: A Geographic-aware Algorithm

Systematic analyses of AISNPs screening and classification algorithms based on genome-wide data for forensic biogeographic ancestry inference

Computation of ancestry scores with mixed families and unrelated individuals

Fast and reliable ancestral reconstruction on ancient genotype data with non-negative Least square and Principal Component Analysis

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

That’s not a Hybrid: How to Distinguish Patterns of Admixture and Isolation-by-Distance

TeraPCA: a fast and scalable software package to study genetic variation in tera-scale genotypes

Genetic ancestry plays a central role in population pharmacogenomics