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 genealogical interpretation of fixed and random effects models of complex traits

Association Test Between Haplotypes and Longitudinal Traits in Complex Pedigrees.

Mixed Linear Model Approaches of Association Mapping for Complex Traits Based on Omics Variants

Mixed Linear Model Approaches for Analyzing Genetic Models of Complex Quantitative Traits

Regression-based Approach for Testing the Association Between Multi-Region Haplotype Configuration and Complex Trait

Mapping the genetic architecture of complex traits in experimental populations.

Disentangling Linkage and Population Structure in Association Mapping

Best Prediction of the Additive Genomic Variance in Random-Effects Models

Gene-based principal component logistic regression model and its application on genome-wide association study]

Probabilistic inference of the genetic architecture underlying functional enrichment of complex traits

Linear Models for Analysis of Multiple Single Nucleotide Polymorphisms with Quantitative Traits in Unrelated Individuals

A logistic mixture model for a family-based association study

A Sparse Graph-Structured Lasso Mixed Model for Genetic Association with Confounding Correction

Causal inference and GWAS: Rubin, Pearl, and Mendelian randomization

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

A Regression-based Approach to Robust Estimation and Inference for Genetic Covariance

A scalable approach for genome-wide inference of ancestral recombination graphs

Integrative functional linear model for genome-wide association studies with multiple traits

Optimal Estimation Of Genetic Relatedness In High-Dimensional Linear Models

A general and efficient representation of ancestral recombination graphs

Simple models of non-random mating and environmental transmission bias standard human genetics statistical methods