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.

Computationally Inferred Genealogical Networks Uncover Long-Term Trends in Assortative Mating

The persistent homology of genealogical networks

Quantitative Analysis of Genealogy Using Digitised Family Trees

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

Statistical Analysis of Genealogical Trees for Polygamic Species

A rapid, accurate approach to inferring pedigrees in endogamous populations

Differentiating mechanism from outcome for ancestry-assortative mating in admixed human populations

Simultaneous Inference of Past Demography and Selection from the Ancestral Recombination Graph under the Beta Coalescent

Data Mining of Online Genealogy Datasets for Revealing Lifespan Patterns in Human Population

Estimating dispersal rates and locating genetic ancestors with genome-wide genealogies

Genetic similarity between relatives provides evidence on the presence and history of assortative mating

Efficient genome ancestry inference in complex pedigrees with inbreeding

Using Genetic Data to Build Intuition about Population History

Iterative Level-0: A new and fast algorithm to traverse mating networks calculating the inbreeding and relationship coefficients

An empirical approach to demographic inference with genomic data

An explanation for the sister repulsion phenomenon in Patterson’s f-statistics

Automatic Discovery of Families of Network Generative Processes

Computing the joint distribution of the total tree length across loci in populations with variable size

CONE: Community Oriented Network Estimation Is a Versatile Framework for Inferring Population Structure in Large-Scale Sequencing Data

Detecting Epistatic Selection with Partially Observed Genotype Data Using Copula Graphical Models