Improved Imputation of Common and Uncommon Single Nucleotide Polymorphisms (SNPs) with a New Reference Set

Statistical imputation of genotype data is an important technique for analysis of genome-wide association studies (GWAS). We have built a reference dataset to improve imputation accuracy for studies of individuals of primarily European descent using genotype data from the Hap1, Omni1, and Omni2.5 human SNP arrays (Illumina). Our dataset contains 2.5-3.1 million variants for 930 European, 157 Asian, and 162 African/African-American individuals. Imputation accuracy of European data from Hap660 or OmniExpress array content, measured by the proportion of variants imputed with R^2^>0.8, improved by 34%, 23% and 12% for variants with MAF of 3%, 5% and 10%, respectively, compared to imputation using publicly available data from 1,000 Genomes and International HapMap projects. The improved accuracy with the use of the new dataset could increase the power for GWAS by as much as 8% relative to genotyping all variants. This reference dataset is available to the scientific community through the NCBI dbGaP portal. Future versions will include additional genotype data as well as non-European populations

Common single nucleotide polymorphisms in genes related to immune function and risk of papillary thyroid cancer

Accumulating evidence suggests that alterations in immune function may be important in the etiology of papillary thyroid cancer (PTC). To identify genetic markers in immune-related pathways, we evaluated 3,985 tag single nucleotide polymorphisms (SNPs) in 230 candidate gene regions (adhesion-extravasation-migration, arachidonic acid metabolism/eicosanoid signaling, complement and coagulation cascade, cytokine signaling, innate pathogen detection and antimicrobials, leukocyte signaling, TNF/NF-kB pathway or other) in a case-control study of 344 PTC cases and 452 controls. We used logistic regression models to estimate odds ratios (OR) and calculate one degree of freedom P values of linear trend (P(SNP-trend)) for the association between genotype (common homozygous, heterozygous, variant homozygous) and risk of PTC. To correct for multiple comparisons, we applied the false discovery rate method (FDR). Gene region- and pathway-level associations (P(Region) and P(Pathway)) were assessed by combining individual P(SNP-trend) values using the adaptive rank truncated product method. Two SNPs (rs6115, rs6112) in the SERPINA5 gene were significantly associated with risk of PTC (P(SNP-FDR)/P(SNP-trend) = 0.02/6×10(−6) and P(SNP-FDR)/P(SNP-trend) = 0.04/2×10(−5), respectively). These associations were independent of a history of autoimmune thyroiditis (OR = 6.4; 95% confidence interval: 3.0–13.4). At the gene region level, SERPINA5 was suggestively associated with risk of PTC (P(Region-FDR)/P(Region = )0.07/0.0003). Overall, the complement and coagulation cascade pathway was the most significant pathway (P(Pathway) = 0.02) associated with PTC risk largely due to the strong effect of SERPINA5. Our results require replication but suggest that the SERPINA5 gene, which codes for the protein C inhibitor involved in many biological processes including inflammation, may be a new susceptibility locus for PTC

Comprehensive resequence analysis of a 97 kb region of chromosome 10q11.2 containing the MSMB gene associated with prostate cancer

Genome-wide association studies of prostate cancer have identified single nucleotide polymorphism (SNP) markers in a region of chromosome 10q11.2, harboring the microseminoprotein-β (MSMB) gene. Both the gene product of MSMB, the prostate secretory protein 94 (PSP94) and its binding protein (PSPBP), have been previously investigated as serum biomarkers for prostate cancer progression. Recent functional work has shown that different alleles of the significantly associated SNP in the promoter of MSMB found to be associated with prostate cancer risk, rs10993994, can influence its expression in tumors and in vitro studies. Since it is plausible that additional variants in this region contribute to the risk of prostate cancer, we have used next-generation sequencing technology to resequence a ~97-kb region that includes the area surrounding MSMB (chr10: 51,168,025–51,265,101) in 36 prostate cancer cases, 26 controls of European origin, and 8 unrelated CEPH individuals in order to identify additional variants to investigate in functional studies. We identified 241 novel polymorphisms within this region, including 142 in the 51-kb block of linkage disequilibrium (LD) that contains rs10993994 and the proximal promoter of MSMB. No sites were observed to be polymorphic within the exons of MSMB

A multi-stage genome-wide association study of bladder cancer identifies multiple susceptibility loci.

We conducted a multi-stage, genome-wide association study of bladder cancer with a primary scan of 591,637 SNPs in 3,532 affected individuals (cases) and 5,120 controls of European descent from five studies followed by a replication strategy, which included 8,382 cases and 48,275 controls from 16 studies. In a combined analysis, we identified three new regions associated with bladder cancer on chromosomes 22q13.1, 19q12 and 2q37.1: rs1014971, (P = 8 × 10⁻¹²) maps to a non-genic region of chromosome 22q13.1, rs8102137 (P = 2 × 10⁻¹¹) on 19q12 maps to CCNE1 and rs11892031 (P = 1 × 10⁻⁷) maps to the UGT1A cluster on 2q37.1. We confirmed four previously identified genome-wide associations on chromosomes 3q28, 4p16.3, 8q24.21 and 8q24.3, validated previous candidate associations for the GSTM1 deletion (P = 4 × 10⁻¹¹) and a tag SNP for NAT2 acetylation status (P = 4 × 10⁻¹¹), and found interactions with smoking in both regions. Our findings on common variants associated with bladder cancer risk should provide new insights into the mechanisms of carcinogenesis

A multi-stage genome-wide association study of bladder cancer identifies multiple susceptibility loci.

We conducted a multi-stage, genome-wide association study of bladder cancer with a primary scan of 591,637 SNPs in 3,532 affected individuals (cases) and 5,120 controls of European descent from five studies followed by a replication strategy, which included 8,382 cases and 48,275 controls from 16 studies. In a combined analysis, we identified three new regions associated with bladder cancer on chromosomes 22q13.1, 19q12 and 2q37.1: rs1014971, (P = 8 × 10⁻¹²) maps to a non-genic region of chromosome 22q13.1, rs8102137 (P = 2 × 10⁻¹¹) on 19q12 maps to CCNE1 and rs11892031 (P = 1 × 10⁻⁷) maps to the UGT1A cluster on 2q37.1. We confirmed four previously identified genome-wide associations on chromosomes 3q28, 4p16.3, 8q24.21 and 8q24.3, validated previous candidate associations for the GSTM1 deletion (P = 4 × 10⁻¹¹) and a tag SNP for NAT2 acetylation status (P = 4 × 10⁻¹¹), and found interactions with smoking in both regions. Our findings on common variants associated with bladder cancer risk should provide new insights into the mechanisms of carcinogenesis

Fine mapping the KLK3 locus on chromosome 19q13.33 associated with prostate cancer susceptibility and PSA levels

Measurements of serum prostate-specific antigen (PSA) protein levels form the basis for a widely used test to screen men for prostate cancer. Germline variants in the gene that encodes the PSA protein (KLK3) have been shown to be associated with both serum PSA levels and prostate cancer. Based on a resequencing analysis of a 56 kb region on chromosome 19q13.33, centered on the KLK3 gene, we fine mapped this locus by genotyping tag SNPs in 3,522 prostate cancer cases and 3,338 controls from five case–control studies. We did not observe a strong association with the KLK3 variant, reported in previous studies to confer risk for prostate cancer (rs2735839; P = 0.20) but did observe three highly correlated SNPs (rs17632542, rs62113212 and rs62113214) associated with prostate cancer [P = 3.41 × 10−4, per-allele trend odds ratio (OR) = 0.77, 95% CI = 0.67–0.89]. The signal was apparent only for nonaggressive prostate cancer cases with Gleason score <7 and disease stage <III (P = 4.72 × 10−5, per-allele trend OR = 0.68, 95% CI = 0.57–0.82) and not for advanced cases with Gleason score >8 or stage ≥III (P = 0.31, per-allele trend OR = 1.12, 95% CI = 0.90–1.40). One of the three highly correlated SNPs, rs17632542, introduces a non-synonymous amino acid change in the KLK3 protein with a predicted benign or neutral functional impact. Baseline PSA levels were 43.7% higher in control subjects with no minor alleles (1.61 ng/ml, 95% CI = 1.49–1.72) than in those with one or more minor alleles at any one of the three SNPs (1.12 ng/ml, 95% CI = 0.96–1.28) (P = 9.70 × 10−5). Together our results suggest that germline KLK3 variants could influence the diagnosis of nonaggressive prostate cancer by influencing the likelihood of biopsy