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.

Number of ancestors and length of identity-by-descent tracks over time

The ancestry of a genome

Genomic runs of homozygosity record population history and consanguinity

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

Identity-by-descent in large samples

Parental relatedness through time revealed by runs of homozygosity in ancient DNA

Inferring demographic history from a spectrum of shared haplotype lengths

Runs of homozygosity: windows into population history and trait architecture

The existence and abundance of ghost ancestors in biparental populations

The geography of recent genetic ancestry across Europe

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

On the number of genealogical ancestors tracing to the source groups of an admixed population

Bursts of coalescence within population pedigrees whenever big families occur

Population Genetics Models of Local Ancestry

Detectability of runs of homozygosity is influenced by analysis parameters and population-specific demographic history

Runs of homozygosity in sub-Saharan African populations provide insights into complex demographic histories

Counting the genetic ancestors from source populations in members of an admixed population

Genealogies of rapidly adapting populations

Estimating effective population size trajectories from time-series Identity-by-Descent (IBD) segments

Age of an allele and gene genealogies of nested subsamples for populations admitting large offspring numbers