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.

Non-alignment comparison of human and high primate genomes

An investigation into inter- and intragenomic variations of graphic genomic signatures

Success of alignment-free oligonucleotide (k-mer) analysis confirms relative importance of genomes not genes in speciation and phylogeny

Human-chimpanzee alignment: Ortholog Exponentials and Paralog Power Laws

Comparative analyses of multi-species sequences from targeted genomic regions

Comparative Gene Prediction Based on Gene Structure Conservation.

Large-Scale Genome Comparison Based on Cumulative Fourier Power and Phase Spectra: Central Moment and Covariance Vector

Hierarchical Structure Analysis Describing Abnormal Base Composition of Genomes.

Peculiar k-mer Spectra Are Correlated with 3D Contact Frequencies and Breakpoint Regions in the Human Genome

Comparative Genomics Approaches to Study Organism Similarities and Differences

Towards a core genome: pairwise similarity searches on interspecific genomic data

Mapping the Space of Genomic Signatures

Comparative analysis of the nucleotide composition biases in exons and introns of human genes

An advanced approach for DNA sequencing and similarities analysis on the basis of groupings of nucleotide bases

Prokaryote Phylogeny Without Sequence Alignment: from Avoidance Signature to Composition Distance.

Alignment-Free Sequence Comparison Based on Next Generation Sequencing Reads: Extended Abstract.

DrosOCB: a high resolution map of conserved non coding sequences in Drosophila

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

Use of 3D chaos game representation to quantify DNA sequence similarity with applications for hierarchical clustering

Determination of Ancestral Alleles for Human Single-Nucleotide Polymorphisms Using High-Density Oligonucleotide Arrays.

Exploring geometry of genome space via Grassmann manifolds