The Cohort Consortium Vitamin D Pooling Project of Rarer Cancers (VDPP), a consortium of 10 prospective cohort studies from the United States, Finland, and China, was formed to examine the associations between circulating 25-hydroxyvitamin D (25(OH)D) concentrations and the risk of rarer cancers. Cases (total n = 5,491) included incident primary endometrial (n = 830), kidney (n = 775), ovarian (n = 516), pancreatic (n = 952), and upper gastrointestinal tract (n = 1,065) cancers and non-Hodgkin lymphoma (n = 1,353) diagnosed in the participating cohorts. At least 1 control was matched to each case on age, date of blood collection (1974-2006), sex, and race/ethnicity (n = 6,714). Covariate data were obtained from each cohort in a standardized manner. The majority of the serum or plasma samples were assayed in a central laboratory using a direct, competitive chemiluminescence immunoassay on the DiaSorin LIAISON platform (DiaSorin, Inc., Stillwater, Minnesota). Masked quality control samples included serum standards from the US National Institute of Standards and Technology. Conditional logistic regression analyses were conducted using clinically defined cutpoints, with 50-<75 nmol/L as the reference category. Meta-analyses were also conducted using inverse-variance weights in random-effects models. This consortium approach permits estimation of the association between 25(OH)D and several rarer cancers with high accuracy and precision across a wide range of 25(OH)D concentrations.

Most men diagnosed with prostate cancer will experience indolent disease; hence, discovering genetic variants that distinguish aggressive from nonaggressive prostate cancer is of critical clinical importance for disease prevention and treatment. In a multistage, case-only genome-wide association study of 12,518 prostate cancer cases, we identify two loci associated with Gleason score, a pathological measure of disease aggressiveness: rs35148638 at 5q14.3 (RASA1, P=6.49 × 10(-9)) and rs78943174 at 3q26.31 (NAALADL2, P=4.18 × 10(-8)). In a stratified case-control analysis, the SNP at 5q14.3 appears specific for aggressive prostate cancer (P=8.85 × 10(-5)) with no association for nonaggressive prostate cancer compared with controls (P=0.57). The proximity of these loci to genes involved in vascular disease suggests potential biological mechanisms worthy of further investigation.

Genome-wide association studies have identified hundreds of genetic variants associated with specific cancers. A few of these risk regions have been associated with more than one cancer site; however, a systematic evaluation of the associations between risk variants for other cancers and lung cancer risk has yet to be performed.

Associations between single nucleotide polymorphisms (SNPs) at 5p15 and multiple cancer types have been reported. We have previously shown evidence for a strong association between prostate cancer (PrCa) risk and rs2242652 at 5p15, intronic in the telomerase reverse transcriptase (TERT) gene that encodes TERT. To comprehensively evaluate the association between genetic variation across this region and PrCa, we performed a fine-mapping analysis by genotyping 134 SNPs using a custom Illumina iSelect array or Sequenom MassArray iPlex, followed by imputation of 1094 SNPs in 22 301 PrCa cases and 22 320 controls in The PRACTICAL consortium. Multiple stepwise logistic regression analysis identified four signals in the promoter or intronic regions of TERT that independently associated with PrCa risk. Gene expression analysis of normal prostate tissue showed evidence that SNPs within one of these regions also associated with TERT expression, providing a potential mechanism for predisposition to disease.

Some observational studies suggest that a higher selenium status is associated with a lower risk of prostate cancer but have been generally too small to provide precise estimates of associations, particularly by disease stage and grade.

Evidence from a small number of studies suggests that longer telomere length measured in peripheral leukocytes is associated with an increased risk of non-Hodgkin lymphoma (NHL). However, these studies may be biased by reverse causation, confounded by unmeasured environmental exposures and might miss time points for which prospective telomere measurement would best reveal a relationship between telomere length and NHL risk. We performed an analysis of genetically inferred telomere length and NHL risk in a study of 10 102 NHL cases of the four most common B-cell histologic types and 9562 controls using a genetic risk score (GRS) comprising nine telomere length-associated single-nucleotide polymorphisms. This approach uses existing genotype data and estimates telomere length by weighing the number of telomere length-associated variant alleles an individual carries with the published change in kb of telomere length. The analysis of the telomere length GRS resulted in an association between longer telomere length and increased NHL risk [four B-cell histologic types combined; odds ratio (OR) = 1.49, 95% CI 1.22-1.82,P-value = 8.5 × 10(-5)]. Subtype-specific analyses indicated that chronic lymphocytic leukemia or small lymphocytic lymphoma (CLL/SLL) was the principal NHL subtype contributing to this association (OR = 2.60, 95% CI 1.93-3.51,P-value = 4.0 × 10(-10)). Significant interactions were observed across strata of sex for CLL/SLL and marginal zone lymphoma subtypes as well as age for the follicular lymphoma subtype. Our results indicate that a genetic background that favors longer telomere length may increase NHL risk, particularly risk of CLL/SLL, and are consistent with earlier studies relating longer telomere length with increased NHL risk.

Recent genome-wide association studies (GWASs) have identified common genetic variants at 5p15.33, 6p21-6p22 and 15q25.1 associated with lung cancer risk. Several other genetic regions including variants of CHEK2 (22q12), TP53BP1 (15q15) and RAD52 (12p13) have been demonstrated to influence lung cancer risk in candidate- or pathway-based analyses. To identify novel risk variants for lung cancer, we performed a meta-analysis of 16 GWASs, totaling 14 900 cases and 29 485 controls of European descent. Our data provided increased support for previously identified risk loci at 5p15 (P = 7.2 × 10(-16)), 6p21 (P = 2.3 × 10(-14)) and 15q25 (P = 2.2 × 10(-63)). Furthermore, we demonstrated histology-specific effects for 5p15, 6p21 and 12p13 loci but not for the 15q25 region. Subgroup analysis also identified a novel disease locus for squamous cell carcinoma at 9p21 (CDKN2A/p16(INK4A)/p14(ARF)/CDKN2B/p15(INK4B)/ANRIL; rs1333040, P = 3.0 × 10(-7)) which was replicated in a series of 5415 Han Chinese (P = 0.03; combined analysis, P = 2.3 × 10(-8)). This large analysis provides additional evidence for the role of inherited genetic susceptibility to lung cancer and insight into biological differences in the development of the different histological types of lung cancer.

Pathway analysis of genome-wide association studies (GWAS) offer a unique opportunity to collectively evaluate genetic variants with effects that are too small to be detected individually. We applied a pathway analysis to a bladder cancer GWAS containing data from 3,532 cases and 5,120 controls of European background (n = 5 studies). Thirteen hundred and ninety-nine pathways were drawn from five publicly available resources (Biocarta, Kegg, NCI-PID, HumanCyc, and Reactome), and we constructed 22 additional candidate pathways previously hypothesized to be related to bladder cancer. In total, 1421 pathways, 5647 genes and ∼90,000 SNPs were included in our study. Logistic regression model adjusting for age, sex, study, DNA source, and smoking status was used to assess the marginal trend effect of SNPs on bladder cancer risk. Two complementary pathway-based methods (gene-set enrichment analysis [GSEA], and adapted rank-truncated product [ARTP]) were used to assess the enrichment of association signals within each pathway. Eighteen pathways were detected by either GSEA or ARTP at P≤0.01. To minimize false positives, we used the I(2) statistic to identify SNPs displaying heterogeneous effects across the five studies. After removing these SNPs, seven pathways ('Aromatic amine metabolism' [P(GSEA) = 0.0100, P(ARTP) = 0.0020], 'NAD biosynthesis' [P(GSEA) = 0.0018, P(ARTP) = 0.0086], 'NAD salvage' [P(ARTP) = 0.0068], 'Clathrin derived vesicle budding' [P(ARTP) = 0.0018], 'Lysosome vesicle biogenesis' [P(GSEA) = 0.0023, P(ARTP)<0.00012], 'Retrograde neurotrophin signaling' [P(GSEA) = 0.00840], and 'Mitotic metaphase/anaphase transition' [P(GSEA) = 0.0040]) remained. These pathways seem to belong to three fundamental cellular processes (metabolic detoxification, mitosis, and clathrin-mediated vesicles). Identification of the aromatic amine metabolism pathway provides support for the ability of this approach to identify pathways with established relevance to bladder carcinogenesis.

Genome-wide association studies have reported 56 independently associated colorectal cancer (CRC) risk variants, most of which are non-coding and believed to exert their effects by modulating gene expression. The computational method PrediXcan uses cis-regulatory variant predictors to impute expression and perform gene-level association tests in GWAS without directly measured transcriptomes. In this study, we used reference datasets from colon (n = 169) and whole blood (n = 922) transcriptomes to test CRC association with genetically determined expression levels in a genome-wide analysis of 12,186 cases and 14,718 controls. Three novel associations were discovered from colon transverse models at FDR ≤ 0.2 and further evaluated in an independent replication including 32,825 cases and 39,933 controls. After adjusting for multiple comparisons, we found statistically significant associations using colon transcriptome models with TRIM4 (discovery P = 2.2 × 10- 4, replication P = 0.01), and PYGL (discovery P = 2.3 × 10- 4, replication P = 6.7 × 10- 4). Interestingly, both genes encode proteins that influence redox homeostasis and are related to cellular metabolic reprogramming in tumors, implicating a novel CRC pathway linked to cell growth and proliferation. Defining CRC risk regions as one megabase up- and downstream of one of the 56 independent risk variants, we defined 44 non-overlapping CRC-risk regions. Among these risk regions, we identified genes associated with CRC (P < 0.05) in 34/44 CRC-risk regions. Importantly, CRC association was found for two genes in the previously reported 2q25 locus, CXCR1 and CXCR2, which are potential cancer therapeutic targets. These findings provide strong candidate genes to prioritize for subsequent laboratory follow-up of GWAS loci. This study is the first to implement PrediXcan in a large colorectal cancer study and findings highlight the utility of integrating transcriptome data in GWAS for discovery of, and biological insight into, risk loci.

Chromosome 8q24 is a susceptibility locus for multiple cancers, including prostate cancer. Here we combine genetic data across the 8q24 susceptibility region from 71,535 prostate cancer cases and 52,935 controls of European ancestry to define the overall contribution of germline variation at 8q24 to prostate cancer risk. We identify 12 independent risk signals for prostate cancer (p < 4.28 × 10-15), including three risk variants that have yet to be reported. From a polygenic risk score (PRS) model, derived to assess the cumulative effect of risk variants at 8q24, men in the top 1% of the PRS have a 4-fold (95%CI = 3.62-4.40) greater risk compared to the population average. These 12 variants account for ~25% of what can be currently explained of the familial risk of prostate cancer by known genetic risk factors. These findings highlight the overwhelming contribution of germline variation at 8q24 on prostate cancer risk which has implications for population risk stratification.

Accurate colorectal cancer (CRC) risk prediction models are critical for identifying individuals at low and high risk of developing CRC, as they can then be offered targeted screening and interventions to address their risks of developing disease (if they are in a high-risk group) and avoid unnecessary screening and interventions (if they are in a low-risk group). As it is likely that thousands of genetic variants contribute to CRC risk, it is clinically important to investigate whether these genetic variants can be used jointly for CRC risk prediction. In this paper, we derived and compared different approaches to generating predictive polygenic risk scores (PRS) from genome-wide association studies (GWASs) including 55,105 CRC-affected case subjects and 65,079 control subjects of European ancestry. We built the PRS in three ways, using (1) 140 previously identified and validated CRC loci; (2) SNP selection based on linkage disequilibrium (LD) clumping followed by machine-learning approaches; and (3) LDpred, a Bayesian approach for genome-wide risk prediction. We tested the PRS in an independent cohort of 101,987 individuals with 1,699 CRC-affected case subjects. The discriminatory accuracy, calculated by the age- and sex-adjusted area under the receiver operating characteristics curve (AUC), was highest for the LDpred-derived PRS (AUC = 0.654) including nearly 1.2 M genetic variants (the proportion of causal genetic variants for CRC assumed to be 0.003), whereas the PRS of the 140 known variants identified from GWASs had the lowest AUC (AUC = 0.629). Based on the LDpred-derived PRS, we are able to identify 30% of individuals without a family history as having risk for CRC similar to those with a family history of CRC, whereas the PRS based on known GWAS variants identified only top 10% as having a similar relative risk. About 90% of these individuals have no family history and would have been considered average risk under current screening guidelines, but might benefit from earlier screening. The developed PRS offers a way for risk-stratified CRC screening and other targeted interventions.

Objective To determine if circulating concentrations of vitamin D are causally associated with risk of cancer.Design Mendelian randomisation study.Setting Large genetic epidemiology networks (the Genetic Associations and Mechanisms in Oncology (GAME-ON), the Genetic and Epidemiology of Colorectal Cancer Consortium (GECCO), and the Prostate Cancer Association Group to Investigate Cancer Associated Alterations in the Genome (PRACTICAL) consortiums, and the MR-Base platform).Participants 70 563 cases of cancer (22 898 prostate cancer, 15 748 breast cancer, 12 537 lung cancer, 11 488 colorectal cancer, 4369 ovarian cancer, 1896 pancreatic cancer, and 1627 neuroblastoma) and 84 418 controls.Exposures Four single nucleotide polymorphisms (rs2282679, rs10741657, rs12785878 and rs6013897) associated with vitamin D were used to define a multi-polymorphism score for circulating 25-hydroxyvitamin D (25(OH)D) concentrations.Main outcomes measures The primary outcomes were the risk of incident colorectal, breast, prostate, ovarian, lung, and pancreatic cancer and neuroblastoma, which was evaluated with an inverse variance weighted average of the associations with specific polymorphisms and a likelihood based approach. Secondary outcomes based on cancer subtypes by sex, anatomic location, stage, and histology were also examined.Results There was little evidence that the multi-polymorphism score of 25(OH)D was associated with risk of any of the seven cancers or their subtypes. Specifically, the odds ratios per 25 nmol/L increase in genetically determined 25(OH)D concentrations were 0.92 (95% confidence interval 0.76 to 1.10) for colorectal cancer, 1.05 (0.89 to 1.24) for breast cancer, 0.89 (0.77 to 1.02) for prostate cancer, and 1.03 (0.87 to 1.23) for lung cancer. The results were consistent with the two different analytical approaches, and the study was powered to detect relative effect sizes of moderate magnitude (for example, 1.20-1.50 per 25 nmol/L decrease in 25(OH)D for most primary cancer outcomes. The Mendelian randomisation assumptions did not seem to be violated.Conclusions There is little evidence for a linear causal association between circulating vitamin D concentration and risk of various types of cancer, though the existence of causal clinically relevant effects of low magnitude cannot be ruled out. These results, in combination with previous literature, provide evidence that population-wide screening for vitamin D deficiency and subsequent widespread vitamin D supplementation should not currently be recommended as a strategy for primary cancer prevention.

Pancreatic cancer (PC) is a complex disease in which both non-genetic and genetic factors interplay. To date, 40 GWAS hits have been associated with PC risk in individuals of European descent, explaining 4.1% of the phenotypic variance.

Evidence for aspirin's chemopreventative properties on colorectal cancer (CRC) is substantial, but its mechanism of action is not well-understood. We combined a proteomic approach with Mendelian randomization (MR) to identify possible new aspirin targets that decrease CRC risk.

The identification of recurrent founder variants in cancer predisposing genes may have important implications for implementing cost-effective targeted genetic screening strategies. In this study, we evaluated the prevalence and relative risk of the CHEK2 recurrent variant c.349A>G in a series of 462 Portuguese patients with early-onset and/or familial/hereditary prostate cancer (PrCa), as well as in the large multicentre PRACTICAL case-control study comprising 55,162 prostate cancer cases and 36,147 controls. Additionally, we investigated the potential shared ancestry of the carriers by performing identity-by-descent, haplotype and age estimation analyses using high-density SNP data from 70 variant carriers belonging to 11 different populations included in the PRACTICAL consortium. The CHEK2 missense variant c.349A>G was found significantly associated with an increased risk for PrCa (OR 1.9; 95% CI: 1.1-3.2). A shared haplotype flanking the variant in all carriers was identified, strongly suggesting a common founder of European origin. Additionally, using two independent statistical algorithms, implemented by DMLE+2.3 and ESTIAGE, we were able to estimate the age of the variant between 2300 and 3125 years. By extending the haplotype analysis to 14 additional carrier families, a shared core haplotype was revealed among all carriers matching the conserved region previously identified in the high-density SNP analysis. These findings are consistent with CHEK2 c.349A>G being a founder variant associated with increased PrCa risk, suggesting its potential usefulness for cost-effective targeted genetic screening in PrCa families.

Impaired lung function is often caused by cigarette smoking, making it challenging to disentangle its role in lung cancer susceptibility. Investigation of the shared genetic basis of these phenotypes in the UK Biobank and International Lung Cancer Consortium (29,266 cases, 56,450 controls) shows that lung cancer is genetically correlated with reduced forced expiratory volume in one second (FEV1: rg = 0.098, p = 2.3 × 10-8) and the ratio of FEV1 to forced vital capacity (FEV1/FVC: rg = 0.137, p = 2.0 × 10-12). Mendelian randomization analyses demonstrate that reduced FEV1 increases squamous cell carcinoma risk (odds ratio (OR) = 1.51, 95% confidence intervals: 1.21-1.88), while reduced FEV1/FVC increases the risk of adenocarcinoma (OR = 1.17, 1.01-1.35) and lung cancer in never smokers (OR = 1.56, 1.05-2.30). These findings support a causal role of pulmonary impairment in lung cancer etiology. Integrative analyses reveal that pulmonary function instruments, including 73 novel variants, influence lung tissue gene expression and implicate immune-related pathways in mediating the observed effects on lung carcinogenesis.

Germline variation and smoking are independently associated with pancreatic ductal adenocarcinoma (PDAC). We conducted genome-wide smoking interaction analysis of PDAC using genotype data from four previous genome-wide association studies in individuals of European ancestry (7,937 cases and 11,774 controls). Examination of expression quantitative trait loci data from the Genotype-Tissue Expression Project followed by colocalization analysis was conducted to determine whether there was support for common SNP(s) underlying the observed associations. Statistical tests were two sided and P < 5 × 10-8 was considered statistically significant. Genome-wide significant evidence of qualitative interaction was identified on chr2q21.3 in intron 5 of the transmembrane protein 163 (TMEM163) and upstream of the cyclin T2 (CCNT2). The most significant SNP using the Empirical Bayes method, in this region that included 45 significantly associated SNPs, was rs1818613 [per allele OR in never smokers 0.87, 95% confidence interval (CI), 0.82-0.93; former smokers 1.00, 95% CI, 0.91-1.07; current smokers 1.25, 95% CI 1.12-1.40, P interaction = 3.08 × 10-9). Examination of the Genotype-Tissue Expression Project data demonstrated an expression quantitative trait locus in this region for TMEM163 and CCNT2 in several tissue types. Colocalization analysis supported a shared SNP, rs842357, in high linkage disequilibrium with rs1818613 (r 2 = 0. 94) driving both the observed interaction and the expression quantitative trait loci signals. Future studies are needed to confirm and understand the differential biologic mechanisms by smoking status that contribute to our PDAC findings. SIGNIFICANCE: This large genome-wide interaction study identifies a susceptibility locus on 2q21.3 that significantly modified PDAC risk by smoking status, providing insight into smoking-associated PDAC, with implications for prevention.

Recessive genetic variation is thought to play a role in non-Hodgkin lymphoma (NHL) etiology. Runs of homozygosity (ROH), defined based on long, continuous segments of homozygous SNPs, can be used to estimate both measured and unmeasured recessive genetic variation. We sought to examine genome-wide homozygosity and NHL risk.

Previous studies had limited power to assess the associations of testosterone with aggressive disease as a primary endpoint. Further, the association of genetically predicted testosterone with aggressive disease is not known. We investigated the associations of calculated free and measured total testosterone and sex hormone-binding globulin (SHBG) with aggressive, overall and early-onset prostate cancer. In blood-based analyses, odds ratios (OR) and 95% confidence intervals (CI) for prostate cancer were estimated using conditional logistic regression from prospective analysis of biomarker concentrations in the Endogenous Hormones, Nutritional Biomarkers and Prostate Cancer Collaborative Group (up to 25 studies, 14 944 cases and 36 752 controls, including 1870 aggressive prostate cancers). In Mendelian randomisation (MR) analyses, using instruments identified using UK Biobank (up to 194 453 men) and outcome data from PRACTICAL (up to 79 148 cases and 61 106 controls, including 15 167 aggressive cancers), ORs were estimated using the inverse-variance weighted method. Free testosterone was associated with aggressive disease in MR analyses (OR per 1 SD = 1.23, 95% CI = 1.08-1.40). In blood-based analyses there was no association with aggressive disease overall, but there was heterogeneity by age at blood collection (OR for men aged <60 years 1.14, CI = 1.02-1.28; Phet  = .0003: inverse association for older ages). Associations for free testosterone were positive for overall prostate cancer (MR: 1.20, 1.08-1.34; blood-based: 1.03, 1.01-1.05) and early-onset prostate cancer (MR: 1.37, 1.09-1.73; blood-based: 1.08, 0.98-1.19). SHBG and total testosterone were inversely associated with overall prostate cancer in blood-based analyses, with null associations in MR analysis. Our results support free testosterone, rather than total testosterone, in the development of prostate cancer, including aggressive subgroups.

Bile acids (BAs), major regulators of the gut microbiota, may play an important role in hepatobiliary cancer etiology. However, few epidemiologic studies have comprehensively examined associations between BAs and liver or biliary tract cancer. In the Alpha-Tocopherol, Beta-Carotene Cancer Prevention (ATBC) study, we designed 1:1 matched, nested, case-control studies of primary liver cancer (n = 201 cases), fatal liver disease (n = 261 cases), and primary biliary tract cancer (n = 138 cases). Using baseline serum collected ≤30 years before diagnosis or death, we measured concentrations of 15 BAs with liquid chromatography-tandem mass spectrometry. We estimated odds ratios (ORs) and 95% confidence intervals (CIs) using multivariable conditional logistic regression models, adjusted for age, education, diabetes status, smoking, alcohol intake, and body mass index. We accounted for multiple comparisons using a false discovery rate (FDR) correction. Comparing the highest to the lowest quartile, seven BAs were positively associated with liver cancer risk, including taurocholic acid (TCA) (OR, 5.62; 95% CI, 2.74-11.52; Q trend < 0.0001), taurochenodeoxycholic acid (TCDCA) (OR, 4.77; 95% CI, 2.26-10.08; Q trend < 0.0001), and glycocholic acid (GCA) OR, 5.30; 95% CI, 2.41-11.66; Q trend < 0.0001), and 11 were positively associated with fatal liver disease risk, including TCDCA (OR, 9.65; 95% CI, 4.41-21.14; Q trend < 0.0001), TCA (OR, 7.45; 95% CI, 3.70-14.97; Q trend < 0.0001), and GCA (OR, 6.98; 95% CI, 3.32-14.68; Q trend < 0.0001). For biliary tract cancer, associations were generally >1 but not significant after FDR correction. Conjugated BAs were strongly associated with increased risk of liver cancer and fatal liver disease, suggesting mechanistic links between BA metabolism and liver cancer or death from liver disease.