| | 1 | = Impute2 pipeline = |
| | 2 | [[br]] |
| | 3 | This page describes the Imputate2 pipeline developed by the GoNL - Impute team. The pipeline performs the following steps: |
| | 4 | * Pre-processing of study data |
| | 5 | 1. Liftover PED/MAP files from build 36 to build 37 |
| | 6 | 1. Quality control of study data |
| | 7 | * Analysis |
| | 8 | 1. Phasing the study data |
| | 9 | 1. Imputing the study data |
| | 10 | A list of used software and versions can be found on the bottom of this page. |
| | 11 | |
| | 12 | All analysis jobs are generated using MOLGENIS compute, more information about MOLGENIS compute and the other analysis pipelines can be found here: https://github.com/molgenis/molgenis-pipelines |
| | 13 | [[br]] |
| | 14 | [[br]] |
| | 15 | == Reference data == |
| | 16 | * ALL_1000G_phase1integrated_v3_impute.tgz from the impute website, this dataset contains 1,092 individuals |
| | 17 | * GoNL release 4, 499 unrelated individuals, 998 haplotypes |
| | 18 | [[br]] |
| | 19 | [[br]] |
| | 20 | == Pre processing of study data == |
| | 21 | === Liftover PED/MAP files from build 36 to build 37 === |
| | 22 | * Location: https://github.com/molgenis/molgenis-pipelines/blob/master/compute4/Liftover_genome_build_PEDMAP/protocols/liftoverPEDMAP.ftl |
| | 23 | [[br]]This protocol applies the following steps: |
| | 24 | * Create a bed file based on MAP file. |
| | 25 | * Map the bed file to build 37 using UCSVC liftover tool. |
| | 26 | * Create list of unmapped SNPs. |
| | 27 | * Create plink mappings file. |
| | 28 | * Create new plink data without unmapped SNPs using Plink. |
| | 29 | [[br]] |
| | 30 | === Quality control of study data === |
| | 31 | * Location: https://github.com/molgenis/molgenis-pipelines/blob/master/compute4/Imputation_imputationtool_studyQC/protocols/studyQC.ftl |
| | 32 | [[br]]This protocol applies QC to the study data using imputationTool. This tool employs a binary format, called [http://genenetwork.nl/wordpress/trityper/ TriTyper] for rapid loading of big genotypic data. ImputationTool performs the following checks: |
| | 33 | 1. Assesses strand alignment of alleles and swap SNPs if needed. For example, if a SNP which is in LD with multiple SNPs has a negative score the alleles are swapped and LD is calculated again. If the score of the SNP is still negative the SNP is removed from the study data. |
| | 34 | 1. Regular Quality Checks done during routine processing of GWAS data. These checks include: Hardy-Weinberq equilibrium should be higher that 10^-4, minor allele frequency should be higher than 0.01 and call rate should be higher than 0.95. If any of these criteria fails the SNP is removed from the study panel. These SNPs are removed because of the high likelihood that they contain erroneous genotypes. |
| | 35 | 1. Simple sanity checks like check if the SNP is present in the reference panel or if it has null alleles. In both cases the SNP is removed from the study panel. |
| | 36 | 1. Check if there are significant differences between the allele frequencies in the two panels. Difference higher than 25% indicates that there is an important qualitative difference that guides this contrast. A possible imputation of this SNP is prone to introduce invalid information. For this reason these SNPs are removed from the study panel. |
| | 37 | 1. Assesses if the haplotype structure is comparable between reference and GWAS data. This is performed by pairwise comparison of r-squared between SNPs in both reference and GWAS. For SNPs in LD (r-squared > 0.1), the allele frequencies are compared. SNPs are removed from the GWAS data when the major allele differs more often than it is identical. |
| | 38 | [[br]] |
| | 39 | [[br]] |
| | 40 | == Analysis == |
| | 41 | === Phasing the study data === |
| | 42 | * Location: https://github.com/molgenis/molgenis-pipelines/blob/master/compute4/Imputation_shapeit_phasing/protocols/shapeitPhasing.ftl |
| | 43 | [[br]]This protocol phases the study data using shapeit, study data can be in the following formats: Ped/Map, Bed/Bim, Gen/Haps. |
| | 44 | [[br]]This protocol applies the following steps: |
| | 45 | * Gather the map files as provided in the ALL_1000G_phase1integrated_v3_impute.tgz from the impute website. |
| | 46 | * Phase the study data using the map files. |
| | 47 | [[br]] |
| | 48 | === Imputing study data === |
| | 49 | * Location: https://github.com/molgenis/molgenis-pipelines/blob/master/compute4/Imputation_impute2/protocols/impute2Imputation.ftl |
| | 50 | [[br]]This protocol imputes the study data in chromosome bins of 5million bases using Impute. Afterwards the results per chromosome bin are merged into full chromosomes. |
| | 51 | [[br]]This protocol applies the following steps: |
| | 52 | * Imputation using Impute2 the phased study data, 1000G map files and the following two parameters over 5million base bins per chromosome: |
| | 53 | * k_hap: 1500, this parameter sets the number of reference haplotypes to use (default: 500). |
| | 54 | * Ne: 20000, this parameter controls the effective population size in the population genetic-model (default: 20000). |
| | 55 | * Concatenating all chromosome bins into a full chromosome. |
| | 56 | [[br]] |
| | 57 | [[br]] |
| | 58 | == Used software == |
| | 59 | Below is a list of used software and versions. |
| | 60 | |
| | 61 | ||Software ||Version || |
| | 62 | ||[http://www.bbmriwiki.nl/svn/ebiogrid/modules/liftOverUcsc/20120905/liftOverUcsc.20120905.tgz UCSC liftOver tool] ||20120905 || |
| | 63 | ||[http://www.bbmriwiki.nl/svn/ebiogrid/modules/plink/1.07-x86_64/plink-1.07-x86_64.tgz Plink] ||v1.07 || |
| | 64 | ||[http://www.bbmriwiki.nl/svn/ebiogrid/scripts/ImputationTool-20120912.zip ImputationTool] ||20120912 || |
| | 65 | ||[http://www.bbmriwiki.nl/svn/ebiogrid/modules/shapeit/v2.r644.linux.x86_64/shapeit.v2.r644.linux.x86_64.tgz Shapeit] ||v2.r644 || |
| | 66 | ||[http://www.bbmriwiki.nl/svn/ebiogrid/modules/impute/v2.3.0_x86_64_static/impute_v2.3.0_x86_64_static.tgz Impute] ||v2.3.0 || |
| | 67 | |
| | 68 | [[br]] |
| | 69 | [[br]] |
| | 70 | For help with trac wiki formatting: http://trac.edgewall.org/wiki/WikiFormatting |
| | 71 | [[br]] |
