Dataset Open Access

Relate-estimated coalescence rates, allele ages, and selection p-values for the 1000 Genomes Project

Speidel, Leo; Forest, Marie; Shi, Sinan; Myers, Simon R.

Overview

Coalescence rates, allele ages, and p-values for evidence of positive selection calculated for 2478 samples of the 1000 Genomes Project using Relate.

We estimated the joint genealogy of all 1000 GP populations and then extracted the embedded genealogy for each population.
For the genealogy of each population, we jointly estimated the population size history and branch lengths. 
Variants segregating in more than one population therefore have correlated but different allele ages in each population.

Please refer to Speidel et al. Nature Genetics (2019) for more details or email leo.speidel@outlook.com for any queries.

Coalescence rates

The zipped directory coalescence_rates.zip contains coalescence rates for 26 populations in the 1000 Genomes Project data set.

  • The .coal files show the haploid coalescence rates, please refer to the Relate documentation for the file format.
  • The popsize.RData file is an R data frame storing the diploid population sizes (0.5/coalescence rate) calculated using the .coal files. The columns of this data frame, named "pop_size", are
    • gens_ago: Time in generations at which epoch starts. (To get years from generations, we multiply by 28.)
    • population_size: Diploid population size in this epoch.
    • population: Name of population 
    • region: Name of region (AFR, AMR, EAS, EUR, SAS)

Allele ages and selection p-values

The zipped directories allele_ages_*.zip contain R data frames for each 1000GP population storing allele ages and selection p-values.
Please note that only mutations that segregate in the population and map to a unique branch in the Relate-estimated marginal trees are included. Selection p-values are only provided for mutations of DAF > 2 that pass quality filters (see Speidel et al., 2019). 

To get an age estimate for a neutral mutation, use 0.5*(lower_age + upper_age). To get years from generations, we multiply by 28.

The columns of these data frames, named "allele_ages", are

  • CHR: chromosome index
  • BP: base-pair position (GRCh37)
  • ID: id of SNP
  • lower_age: Age in generations of coalescence event at the lower end of the branch onto which the mutation maps
  • upper_age: Age in generations of coalescence event at the upper end of the branch onto which the mutation maps
  • ancestral/derived: Ancestral/derived allele
  • upstream: Upstream (5') allele
  • downstream: Downstream (3') allele
  • DAF: Derived-allele frequency
  • pvalue: log10 p-value for selection evidence

For R object files, use load() to load data frames into R.
Files (7.4 GB)
Name Size
allele_ages_AFR.zip
md5:c3d94e2084205dcb7101bd51d3660409
2.8 GB Download
allele_ages_AMR.zip
md5:f7da0238962e45a2446fbe9d88b6fef8
1.1 GB Download
allele_ages_EAS.zip
md5:f1e5241291fe674fb5c3d4390c9d6663
1.0 GB Download
allele_ages_EUR.zip
md5:c637c78c2248ca46144eb40d2c4f6c4d
1.2 GB Download
allele_ages_SAS.zip
md5:9ee3fada4ae59271992d8288ee8a82d1
1.2 GB Download
coalescence_rates.zip
md5:b7b247836e1d078de755824fcecfc75b
19.7 kB Download
  • Speidel et al., Nature Genetics 2019, A method for genome-wide genealogy estimation for thousands of samples. https://doi.org/10.1038/s41588-019-0484-x

256
115
views
downloads
All versions This version
Views 256256
Downloads 115115
Data volume 155.4 GB155.4 GB
Unique views 237237
Unique downloads 4545

Share

Cite as