WHI Harmonized and Imputed GWAS Data

This substudy phs000746 WHI Imputation is the effort of a joint imputation and harmonization effort for many of the Genome Wide Association Studies (GWAS) within the Women's Health Initiative (WHI) Clinical Trials and Observational Studies. Summary level phenotypes for the WHI Cohort study participants can be viewed at the top-level study page [phs000200](./study.cgi?study_id=phs000200) WHI Cohort. Individual level phenotype data and molecular data for all WHI top-level study and substudies are available by requesting Authorized Access to the WHI Cohort study [phs000200](./study.cgi?study_id=phs000200). These studies jointly involved over 30,000 samples, alignment ("flipping") to the same reference panel, imputation to the 1000 genomes, identification of genetically related individuals, and computations of principal components and comparison with self-reported ethnicity. The harmonization/imputation effort involves 6 GWAS, as described in the table below. (WHI Study #) Hip Fracture GWAS (BA3) [SHARe](./study.cgi?study_id=phs000386) (M5) [GARNET](./study.cgi?study_id=phs000315) (M13) [WHIMS+](./study.cgi?study_id=phs000675) (W63) GECCO (AS224) MOPMAP (AS264) Hip Fracture GWAS (BA3) [GARNET](./study.cgi?study_id=phs000315) (M13) GECCO (AS224) Directly genotyped GWAS data on dbGaP No Yes Yes Data uploaded to dbGaP June 2013 Yes No No Yes Yes Study that funded GWAS WHI-BAA3 National Heart, Lung, and Blood Institute (NHLBI) NHLBI WHI-GARNET - National Human Genome Research Institute (NHGRI) NHLBI WHI-AS224 - National Cancer Institute (NCI) WHI-AS264 - National Institute of Environmental Health Sciences (NIEHS) and Univ of North Carolina WHI-BAA3 National Heart, Lung, and Blood Institute (NHLBI) WHI-GARNET - National Human Genome Research Institute (NHGRI) WHI-AS224 - National Cancer Institute (NCI) GWAS platform Illumina 550K and 610K Affymetrix 6.0 Illumina HumanOmni1-Quad v1-0 B HumanOmniExpress Exome-8v1_B Illumina 610 and Cytochip 370K Affymetrix Gene Titan, Axiom Genome-Wide Human CEU I Array Plate Illumina 550K and 610K Illumina HumanOmni1-Quad v1-0 B Illumina 610 and Cytochip 370K Design Case-control Cohort Case-control (4 case groups) Cohort Case-control Case-Control Case-control Case-control (4 case groups) Case-control Phenotype for cases Hip Fracture NA Type 2 Diabetes, Myocardial Infarction, Stroke, Venous Thrombosis NA Colorectal cancer Ventricular Ectopy (ever) Hip Fracture Type 2 Diabetes, Myocardial Infarction, Stroke, Venous Thrombosis Colorectal cancer Other sample details NA Minorities Hormone Therapy Clinical Trial Hormone Therapy Clinical Trial NA Controls selected within centers, years, seasons and visit years in which cases originated NA Hormone Therapy Clinical Trial NA Ethnicity Mostly white African American and Hispanic Mostly white White White White Mostly white Mostly white White Sample size* 3690 11992 4883 5687 2493 3069 3690 4883 2493 *The sample sizes are the number of samples after QC that are available on dbGaP. Note that there are some subjects that are in multiple studies, as detailed below. #For some of the data files these two platforms are considered different studies. **Initial QC** Initial QC had already been carried out on each of the GWAS studies, using the GENEVA protocol or protocols that were vey similar. Some of the pertinent QC parameters used for each of the studies are shown in the table below. Hip Fracture GWAS SHARe GARNET WHIMS+ GECCO MOPMAP SHARe WHIMS+ MOPMAP Minimal sample call rate 98% 95% 98% 97% 97% 95% 98% 98% 97% Minimal SNP call rate 98% 90% 98% 98% 98% 90% 98% 98% 98% Hardy Weinberg P-value cut-off below which SNPs are excluded 1e-4 1e-6 1e-4 1e-4 1e-4 1e-6 1e-4 1e-4 1e-4 Samples used for Hardy Weinberg calculations Controls of European-ancestry All samples, separate for Hispanics and African Americans Unrelated controls of European-ancestry All Controls All Controls of European-ancestry Unrelated controls of European-ancestry Controls Minimum allele frequency cut-off 1% 1% None 1% 5% 0.5% 1% None 5% **Imputation** The imputation was done using the following procedures. - Match the strand of the GWAS data with the 1KGP data by comparing the letters of the alleles (ambiguous A/T or C/G SNPs were excluded). - We used the 1kGP reference panel (1092 samples; v2.20101123 for GECCO; v3.20101123 for GARNET, HIPFX, MOPMAP, WHIMS+). - The GWAS data were first split into chunks. Each chunk has 10000 SNPs and neighboring chunks have 1000 overlapping SNPs. Then all chunks were phased using Beagle and then combined using mergebeaglechunks.jar (available from the BEAGLE website). - An autoclip file was created for minimac to specify what the range of the chunks (start and stop) and the SNPs to be imputed within the chunk (core_start and core_end) so that no SNP needs to be be imputed twice. All chunks were imputed into 1kGP using minimac. - SNPs that could not be imputed with high enough confidence (cut-off R2>0.1) were omitted for that particular study (but still appear as columns of missing data in files if they were kept in the other studies, to facilitate alignment). - We did not impute the X chromosome. The SHARe study was independently imputed to the same reference panel. The procedures used were similar to those listed above, except that MACH was used to carry out the imputation. **Harmonization** - A panel of 5665 SNPs was used for checking the pairwise concordance among all samples in GARNET, GECCO, HIPFX, SHARe, WHIMS+ and MOPMAP. - The same panel of SNPs was used for principal component check together with HapMap samples to identify ethnicity outlier. - The same panel of SNPs were used for checking IBD in plink to identify relatedness among samples. - Another PC analysis was done for combined samples (after removing of ineligible duplicates) in all studies then the resulting PCs were mapped back to samples within each study. - A netcdf file of imputed results was created for each chromosome in each study. Different studies have the same set of SNPs - SNPs that were not successfully imputed in a particular study but are in other studies are listed as missing values. - A SNP info file was also created along with each netcdf file describing the SNP name, chromosome, position, count allele, alternative allele, count allele frequency, and imputation quality for each SNP. **Duplicates** As subjects for each of these GWAS were selected independently we checked for duplicates between the studies. We removed a small number of samples that - were supposed to be duplicates but had a concordance rate smaller than 90%; and - appeared duplicates but were from unrelated individuals, who appeared not to be monozygotic twins. We kept samples that - were monozygotic twins (see relatedness below), and - duplicates between studies in our data sets. There currently are 29846 unique subjects in the data. **Relatedness** We carried out an IBD analysis using a subset of 5665 SNPS and the Plink package. We used the results to identify 42 parent off-spring pairs a