Genotype imputation was done with the population-specific SISu v4.2 reference panel.
The reference panel variant call set was produced with the GATK HaplotypeCaller algorithm by following GATK best practices for variant calling.
Genotype-, sample- and variant-wise QC was carried out iteratively by using the Hail framework v0.2 and the resulting high-quality WGS data for 8,554 individuals were phased with Eagle 2.3.5 as described in the previous section.
Genotype imputation was carried out by using the population-specific SISu v4.2 imputation reference panel with Beagle 4.1 (version 27Jan18.7e1) as described in the following protocol: dx.doi.org/10.17504/protocols.io.xbgfijw.
Post-imputation quality control involved checking the expected conformity of the imputation INFO-value distribution, MAF differences between the target dataset and the imputation reference panel and checking chromosomal continuity of the imputed genotype calls.
SISu v4.2 consists of 8,554 WGS of Finnish individuals from 5 research cohorts from:
METSIM (PIs Markku Laakso and Mike Boehnke)
FINRISK (PI Pekka Jousilahti)
Corogene (PI Juha Sinisalo)
Biobank of Eastern Finland (PI Arto Mannermaa)
Finnish EUFAM Dyslipidemia Study (PIs Marja-Riitta Taskinen and Samuli Ripatti)
High-coverage (25x) WGS data used to develop the SISu v4.2 reference panel were generated at the McDonnell Genome Institute at Washington University (PIs Ira Hall and Nathan Stitziel).
Hail v0.2
Cromwell-42
Wdltool-0.14
Plink 1.9 and 2.0
BCFtools 1.7 and 1.9
Eagle 2.3.5
Beagle 4.1 (version 27Jan18.7e1)
R 3.4.1 (packages: data.table 1.10.4, sm 2.2-5.4)
The BCOR files were created using LDstore from the Finnish SISu panel v4.2.
The panel has been divided per chromosome. For example, to use the LD information in the first chromosome, FG_LD_chr1.bcor
would be the file to use.
number of samples: 3775
window size: 1500 kb
accuracy: low
number of threads: 96
LD threshold to include correlations: 0.05
LDstore v1.1 can be downloaded via:
And an example to extract variant range 20 Mb - 50 Mb from chromosome 7 is as follows:
It is not preferred to use these LD estimate files for e.g. fine-mapping, since many of the fine-mapping methods (e.g. SuSiE) require in-sample LD information for good results!
Chip genotype data processing and QC Samples were genotyped with Illumina (Illumina Inc., San Diego, CA, USA) and Affymetrix arrays (Thermo Fisher Scientific, Santa Clara, CA, USA).
Genotype calls were made with GenCall and zCall algorithms for Illumina and AxiomGT1 algorithm for Affymetrix data.
Chip genotyping data produced with previous chip platforms and reference genome builds were lifted over to build version 38 (GRCh38/hg38) following the protocol described here: dx.doi.org/10.17504/protocols.io.xbhfij6.
In sample-wise quality control steps, individuals with ambiguous gender, high genotype missingness (>5%), excess heterozygosity (+-4SD) and non-Finnish ancestry were excluded. In variant-wise quality control steps, variants with high missingness (>2%), low HWE P-value (<1e-6) and low minor allele count (MAC<3) were excluded.
Before imputation, chip-genotyped samples were pre-phased with Eagle 2.3.5 using the default parameters, except the number of conditioning haplotypes, which was set to 20,000.
FinnGen individuals were genotyped with Illumina and Affymetrix chip arrays (Illumina Inc., San Diego, and Thermo Fisher Scientific, Santa Clara, CA, USA).
Chip genotype data were imputed using the population-specific SISu v4.2 imputation reference panel of 8,554 whole genomes.
Merged imputed genotype data is composed of 116 data sets that include samples from multiple cohorts.
Total number of individuals: 473,681
Total number of variants (merged set): 21,311,942
Reference assembly: GRCh38/hg38