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.

Statistical Analysis of Genealogical Trees for Polygamic Species

Genealogical trees from genetic distances

Probabilistic reconstruction of genealogies for polyploid plant species

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

Polygene: Population Genetics Analyses for Autopolyploids Based on Allelic Phenotypes

Statistical Phylogenetic Tree Analysis Using Differences of Means

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

Species, Clusters and the 'Tree of Life': A graph-theoretic perspective

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

Quantitative Analysis of Genealogy Using Digitised Family Trees

Reconciliation of Gene and Species Trees With Polytomies

Computing the probability of gene trees concordant with the species tree in the multispecies coalescent

Markovian simulation for ancestors trees

Branching trees I: Concatenation and infinite divisibility

DISCO: Species Tree Inference using Multicopy Gene Family Tree Decomposition

[The Application of Statistical Approaches in Studying Population Genetics with STR Data].

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

Statistical Inference of Allopolyploid Species Networks in the Presence of Incomplete Lineage Sorting

Distance to the stochastic part of phylogenetic varieties

Genomic data analysis in tree spaces

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