Data for multiple common susceptibility alleles for breast cancer may be combined to identify women at different levels of breast cancer risk. Such stratification could guide preventive and screening strategies. However, empirical evidence for genetic risk stratification is lacking.

Genome-wide association studies (GWAS) have identified common variants that predispose individuals to a higher body mass index (BMI), an independent risk factor for endometrial cancer. Composite genotype risk scores (GRS) based on the joint effect of published BMI risk loci were used to explore whether endometrial cancer shares a genetic background with obesity. Genotype and risk factor data were available on 3,376 endometrial cancer case and 3,867 control participants of European ancestry from the Epidemiology of Endometrial Cancer Consortium GWAS. A BMI GRS was calculated by summing the number of BMI risk alleles at 97 independent loci. For exploratory analyses, additional GRSs were based on subsets of risk loci within putative etiologic BMI pathways. The BMI GRS was statistically significantly associated with endometrial cancer risk (P = 0.002). For every 10 BMI risk alleles a woman had a 13% increased endometrial cancer risk (95% CI: 4%, 22%). However, after adjusting for BMI, the BMI GRS was no longer associated with risk (per 10 BMI risk alleles OR = 0.99, 95% CI: 0.91, 1.07; P = 0.78). Heterogeneity by BMI did not reach statistical significance (P = 0.06), and no effect modification was noted by age, GWAS Stage, study design or between studies (P≥0.58). In exploratory analyses, the GRS defined by variants at loci containing monogenic obesity syndrome genes was associated with reduced endometrial cancer risk independent of BMI (per BMI risk allele OR = 0.92, 95% CI: 0.88, 0.96; P = 2.1 x 10-5). Possessing a large number of BMI risk alleles does not increase endometrial cancer risk above that conferred by excess body weight among women of European descent. Thus, the GRS based on all current established BMI loci does not provide added value independent of BMI. Future studies are required to validate the unexpected observed relation between monogenic obesity syndrome genetic variants and endometrial cancer risk.

Despite the high prevalence of prostate cancer in older men, the predictive value of a polygenic risk score (PRS) remains uncertain in men aged ≥70 years. We used a 6.6 million-variant PRS to predict the risk of incident prostate cancer in a prospective study of 5701 men of European descent aged ≥70 years (mean age 75 years) enrolled in the ASPirin in Reducing Events in the Elderly (ASPREE) clinical trial. The study endpoint was prostate cancer, including metastatic or non-metastatic disease, confirmed by an expert panel. After excluding participants with a history of prostate cancer at enrolment, we used a multivariable Cox proportional hazards model to assess the association between the PRS and incident prostate cancer risk, adjusting for covariates. Additionally, we examined the distribution of Gleason grade groups by PRS group to determine if a higher PRS was associated with higher grade disease. We tested for interaction between the PRS and aspirin treatment. Logistic regression was used to independently assess the association of the PRS with prevalent (pre-trial) prostate cancer, reported in medical histories. During a median follow-up time of 4.6 years, 218 of the 5701 participants (3.8%) were diagnosed with prostate cancer. The PRS predicted incident risk with a hazard ratio (HR) of 1.52 per standard deviation (SD) (95% confidence interval (CI) 1.33-1.74, p < 0.001). Men in the top quintile of the PRS distribution had an almost three times higher risk of prostate cancer than men in the lowest quintile (HR = 2.99 (95% CI 1.90-4.27), p < 0.001). However, a higher PRS was not associated with a higher Gleason grade groups. We found no interaction between aspirin treatment and the PRS for prostate cancer risk. The PRS was also associated with prevalent prostate cancer (odds ratio = 1.80 per SD (95% CI 1.65-1.96), p < 0.001).While a PRS for prostate cancer is strongly associated with incident risk in men aged ≥70 years, the clinical utility of the PRS as a biomarker is currently limited by its inability to select for clinically significant disease.

Genomic risk prediction models for breast cancer (BC) have been predominantly developed with data from women aged 40-69 years. Prospective studies of older women aged ≥70 years have been limited. We assessed the effect of a 313-variant polygenic risk score (PRS) for BC in 6339 older women aged ≥70 years (mean age 75 years) enrolled into the ASPREE trial, a randomized double-blind placebo-controlled clinical trial investigating the effect of daily 100 mg aspirin on disability-free survival. We evaluated incident BC diagnoses over a median follow-up time of 4.7 years. A multivariable Cox regression model including conventional BC risk factors was applied to prospective data, and re-evaluated after adding the PRS. We also assessed the association of rare pathogenic variants (PVs) in BC susceptibility genes (BRCA1/BRCA2/PALB2/CHEK2/ATM). The PRS, as a continuous variable, was an independent predictor of incident BC (hazard ratio (HR) per standard deviation (SD) = 1.4, 95% confidence interval (CI) 1.3-1.6) and hormone receptor (ER/PR)-positive disease (HR = 1.5 (CI 1.2-1.9)). Women in the top quintile of the PRS distribution had over two-fold higher risk of BC than women in the lowest quintile (HR = 2.2 (CI 1.2-3.9)). The concordance index of the model without the PRS was 0.62 (95% CI 0.56-0.68), which improved after addition of the PRS to 0.65 (95% CI 0.59-0.71). Among 41 (0.6%) carriers of PVs in BC susceptibility genes, we observed no incident BC diagnoses. Our study demonstrates that a PRS predicts incident BC risk in women aged 70 years and older, suggesting potential clinical utility extends to this older age group.

Genome-wide association studies (GWASs) of follicular lymphoma (FL) have previously identified human leukocyte antigen (HLA) gene variants. To identify additional FL susceptibility loci, we conducted a large-scale two-stage GWAS in 4,523 case subjects and 13,344 control subjects of European ancestry. Five non-HLA loci were associated with FL risk: 11q23.3 (rs4938573, p = 5.79 × 10(-20)) near CXCR5; 11q24.3 (rs4937362, p = 6.76 × 10(-11)) near ETS1; 3q28 (rs6444305, p = 1.10 × 10(-10)) in LPP; 18q21.33 (rs17749561, p = 8.28 × 10(-10)) near BCL2; and 8q24.21 (rs13254990, p = 1.06 × 10(-8)) near PVT1. In an analysis of the HLA region, we identified four linked HLA-DRβ1 multiallelic amino acids at positions 11, 13, 28, and 30 that were associated with FL risk (pomnibus = 4.20 × 10(-67) to 2.67 × 10(-70)). Additional independent signals included rs17203612 in HLA class II (odds ratio [OR(per-allele)] = 1.44; p = 4.59 × 10(-16)) and rs3130437 in HLA class I (OR(per-allele) = 1.23; p = 8.23 × 10(-9)). Our findings further expand the number of loci associated with FL and provide evidence that multiple common variants outside the HLA region make a significant contribution to FL risk.

Previous genome-wide association studies (GWAS) have identified six risk loci for renal cell carcinoma (RCC). We conducted a meta-analysis of two new scans of 5,198 cases and 7,331 controls together with four existing scans, totalling 10,784 cases and 20,406 controls of European ancestry. Twenty-four loci were tested in an additional 3,182 cases and 6,301 controls. We confirm the six known RCC risk loci and identify seven new loci at 1p32.3 (rs4381241, P=3.1 × 10-10), 3p22.1 (rs67311347, P=2.5 × 10-8), 3q26.2 (rs10936602, P=8.8 × 10-9), 8p21.3 (rs2241261, P=5.8 × 10-9), 10q24.33-q25.1 (rs11813268, P=3.9 × 10-8), 11q22.3 (rs74911261, P=2.1 × 10-10) and 14q24.2 (rs4903064, P=2.2 × 10-24). Expression quantitative trait analyses suggest plausible candidate genes at these regions that may contribute to RCC susceptibility.

Transcriptome-wide association studies (TWAS) test the association between traits and genetically predicted gene expression levels. The power of a TWAS depends in part on the strength of the correlation between a genetic predictor of gene expression and the causally relevant gene expression values. Consequently, TWAS power can be low when expression quantitative trait locus (eQTL) data used to train the genetic predictors have small sample sizes, or when data from causally relevant tissues are not available. Here, we propose to address these issues by integrating multiple tissues in the TWAS using sparse canonical correlation analysis (sCCA). We show that sCCA-TWAS combined with single-tissue TWAS using an aggregate Cauchy association test (ACAT) outperforms traditional single-tissue TWAS. In empirically motivated simulations, the sCCA+ACAT approach yielded the highest power to detect a gene associated with phenotype, even when expression in the causal tissue was not directly measured, while controlling the Type I error when there is no association between gene expression and phenotype. For example, when gene expression explains 2% of the variability in outcome, and the GWAS sample size is 20,000, the average power difference between the ACAT combined test of sCCA features and single-tissue, versus single-tissue combined with Generalized Berk-Jones (GBJ) method, single-tissue combined with S-MultiXcan, UTMOST, or summarizing cross-tissue expression patterns using Principal Component Analysis (PCA) approaches was 5%, 8%, 5% and 38%, respectively. The gain in power is likely due to sCCA cross-tissue features being more likely to be detectably heritable. When applied to publicly available summary statistics from 10 complex traits, the sCCA+ACAT test was able to increase the number of testable genes and identify on average an additional 400 additional gene-trait associations that single-trait TWAS missed. Our results suggest that aggregating eQTL data across multiple tissues using sCCA can improve the sensitivity of TWAS while controlling for the false positive rate.

The Nurses' Health Study (NHS), Nurses' Health Study II (NHSII), Health Professionals Follow Up Study (HPFS) and the Physicians Health Study (PHS) have collected detailed longitudinal data on multiple exposures and traits for approximately 310,000 study participants over the last 35 years. Over 160,000 study participants across the cohorts have donated a DNA sample and to date, 20,691 subjects have been genotyped as part of genome-wide association studies (GWAS) of twelve primary outcomes. However, these studies utilized six different GWAS arrays making it difficult to conduct analyses of secondary phenotypes or share controls across studies. To allow for secondary analyses of these data, we have created three new datasets merged by platform family and performed imputation using a common reference panel, the 1,000 Genomes Phase I release. Here, we describe the methodology behind the data merging and imputation and present imputation quality statistics and association results from two GWAS of secondary phenotypes (body mass index (BMI) and venous thromboembolism (VTE)). We observed the strongest BMI association for the FTO SNP rs55872725 (β = 0.45, p = 3.48x10-22), and using a significance level of p = 0.05, we replicated 19 out of 32 known BMI SNPs. For VTE, we observed the strongest association for the rs2040445 SNP (OR = 2.17, 95% CI: 1.79-2.63, p = 2.70x10-15), located downstream of F5 and also observed significant associations for the known ABO and F11 regions. This pooled resource can be used to maximize power in GWAS of phenotypes collected across the cohorts and for studying gene-environment interactions as well as rare phenotypes and genotypes.

Breast cancer in situ (BCIS) diagnoses, a precursor lesion for invasive breast cancer, comprise about 20 % of all breast cancers (BC) in countries with screening programs. Family history of BC is considered one of the strongest risk factors for BCIS.

Genome-wide association studies (GWAS) have so far reported 12 loci associated with serous epithelial ovarian cancer (EOC) risk. We hypothesized that some of these loci function through nearby transcription factor (TF) genes and that putative target genes of these TFs as identified by coexpression may also be enriched for additional EOC risk associations.

Prostate cancer is a highly heritable disease with large disparities in incidence rates across ancestry populations. We conducted a multiancestry meta-analysis of prostate cancer genome-wide association studies (107,247 cases and 127,006 controls) and identified 86 new genetic risk variants independently associated with prostate cancer risk, bringing the total to 269 known risk variants. The top genetic risk score (GRS) decile was associated with odds ratios that ranged from 5.06 (95% confidence interval (CI), 4.84-5.29) for men of European ancestry to 3.74 (95% CI, 3.36-4.17) for men of African ancestry. Men of African ancestry were estimated to have a mean GRS that was 2.18-times higher (95% CI, 2.14-2.22), and men of East Asian ancestry 0.73-times lower (95% CI, 0.71-0.76), than men of European ancestry. These findings support the role of germline variation contributing to population differences in prostate cancer risk, with the GRS offering an approach for personalized risk prediction.

Genome-wide association studies have identified several common susceptibility alleles for epithelial ovarian cancer (EOC). To further understand EOC susceptibility, we examined previously ungenotyped candidate variants, including uncommon variants and those residing within known susceptibility loci.

Genome-wide association studies (GWAS) have identified numerous common prostate cancer (PrCa) susceptibility loci. We have fine-mapped 64 GWAS regions known at the conclusion of the iCOGS study using large-scale genotyping and imputation in 25 723 PrCa cases and 26 274 controls of European ancestry. We detected evidence for multiple independent signals at 16 regions, 12 of which contained additional newly identified significant associations. A single signal comprising a spectrum of correlated variation was observed at 39 regions; 35 of which are now described by a novel more significantly associated lead SNP, while the originally reported variant remained as the lead SNP only in 4 regions. We also confirmed two association signals in Europeans that had been previously reported only in East-Asian GWAS. Based on statistical evidence and linkage disequilibrium (LD) structure, we have curated and narrowed down the list of the most likely candidate causal variants for each region. Functional annotation using data from ENCODE filtered for PrCa cell lines and eQTL analysis demonstrated significant enrichment for overlap with bio-features within this set. By incorporating the novel risk variants identified here alongside the refined data for existing association signals, we estimate that these loci now explain ∼38.9% of the familial relative risk of PrCa, an 8.9% improvement over the previously reported GWAS tag SNPs. This suggests that a significant fraction of the heritability of PrCa may have been hidden during the discovery phase of GWAS, in particular due to the presence of multiple independent signals within the same region.

Disruption in circadian gene expression, whether due to genetic variation or environmental factors (e.g., light at night, shiftwork), is associated with increased incidence of breast, prostate, gastrointestinal and hematologic cancers and gliomas. Circadian genes are highly expressed in the ovaries where they regulate ovulation; circadian disruption is associated with several ovarian cancer risk factors (e.g., endometriosis). However, no studies have examined variation in germline circadian genes as predictors of ovarian cancer risk and invasiveness. The goal of the current study was to examine single nucleotide polymorphisms (SNPs) in circadian genes BMAL1, CRY2, CSNK1E, NPAS2, PER3, REV1 and TIMELESS and downstream transcription factors KLF10 and SENP3 as predictors of risk of epithelial ovarian cancer (EOC) and histopathologic subtypes. The study included a test set of 3,761 EOC cases and 2,722 controls and a validation set of 44,308 samples including 18,174 (10,316 serous) cases and 26,134 controls from 43 studies participating in the Ovarian Cancer Association Consortium (OCAC). Analysis of genotype data from 36 genotyped SNPs and 4600 imputed SNPs indicated that the most significant association was rs117104877 in BMAL1 (OR = 0.79, 95% CI = 0.68-0.90, p = 5.59 × 10-4]. Functional analysis revealed a significant down regulation of BMAL1 expression following cMYC overexpression and increasing transformation in ovarian surface epithelial (OSE) cells as well as alternative splicing of BMAL1 exons in ovarian and granulosa cells. These results suggest that variation in circadian genes, and specifically BMAL1, may be associated with risk of ovarian cancer, likely through disruption of hormonal pathways.

Defective cellular transport processes can lead to aberrant accumulation of trace elements, iron, small molecules and hormones in the cell, which in turn may promote the formation of reactive oxygen species, promoting DNA damage and aberrant expression of key regulatory cancer genes. As DNA damage and uncontrolled proliferation are hallmarks of cancer, including epithelial ovarian cancer (EOC), we hypothesized that inherited variation in the cellular transport genes contributes to EOC risk.

Candidate variant association studies have been largely unsuccessful in identifying common breast cancer susceptibility variants, although most studies have been underpowered to detect associations of a realistic magnitude. We assessed 41 common non-synonymous single-nucleotide polymorphisms (nsSNPs) for which evidence of association with breast cancer risk had been previously reported. Case-control data were combined from 38 studies of white European women (46 450 cases and 42 600 controls) and analyzed using unconditional logistic regression. Strong evidence of association was observed for three nsSNPs: ATXN7-K264R at 3p21 [rs1053338, per allele OR = 1.07, 95% confidence interval (CI) = 1.04-1.10, P = 2.9 × 10(-6)], AKAP9-M463I at 7q21 (rs6964587, OR = 1.05, 95% CI = 1.03-1.07, P = 1.7 × 10(-6)) and NEK10-L513S at 3p24 (rs10510592, OR = 1.10, 95% CI = 1.07-1.12, P = 5.1 × 10(-17)). The first two associations reached genome-wide statistical significance in a combined analysis of available data, including independent data from nine genome-wide association studies (GWASs): for ATXN7-K264R, OR = 1.07 (95% CI = 1.05-1.10, P = 1.0 × 10(-8)); for AKAP9-M463I, OR = 1.05 (95% CI = 1.04-1.07, P = 2.0 × 10(-10)). Further analysis of other common variants in these two regions suggested that intronic SNPs nearby are more strongly associated with disease risk. We have thus identified a novel susceptibility locus at 3p21, and confirmed previous suggestive evidence that rs6964587 at 7q21 is associated with risk. The third locus, rs10510592, is located in an established breast cancer susceptibility region; the association was substantially attenuated after adjustment for the known GWAS hit. Thus, each of the associated nsSNPs is likely to be a marker for another, non-coding, variant causally related to breast cancer risk. Further fine-mapping and functional studies are required to identify the underlying risk-modifying variants and the genes through which they act.

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.

To further dissect the genetic architecture of colorectal cancer (CRC), we performed whole-genome sequencing of 1,439 cases and 720 controls, imputed discovered sequence variants and Haplotype Reference Consortium panel variants into genome-wide association study data, and tested for association in 34,869 cases and 29,051 controls. Findings were followed up in an additional 23,262 cases and 38,296 controls. We discovered a strongly protective 0.3% frequency variant signal at CHD1. In a combined meta-analysis of 125,478 individuals, we identified 40 new independent signals at P < 5 × 10-8, bringing the number of known independent signals for CRC to ~100. New signals implicate lower-frequency variants, Krüppel-like factors, Hedgehog signaling, Hippo-YAP signaling, long noncoding RNAs and somatic drivers, and support a role for immune function. Heritability analyses suggest that CRC risk is highly polygenic, and larger, more comprehensive studies enabling rare variant analysis will improve understanding of biology underlying this risk and influence personalized screening strategies and drug development.

Mammographic density is an established breast cancer risk factor with a strong genetic component and can be increased in women using menopausal hormone therapy (MHT). Here, we aimed to identify genetic variants that may modify the association between MHT use and mammographic density.

The genetic contributions to breast cancer development among Latinas are not well understood. Here we carry out a genome-wide association study of breast cancer in Latinas and identify a genome-wide significant risk variant, located 5' of the Estrogen Receptor 1 gene (ESR1; 6q25 region). The minor allele for this variant is strongly protective (rs140068132: odds ratio (OR) 0.60, 95% confidence interval (CI) 0.53-0.67, P=9 × 10(-18)), originates from Indigenous Americans and is uncorrelated with previously reported risk variants at 6q25. The association is stronger for oestrogen receptor-negative disease (OR 0.34, 95% CI 0.21-0.54) than oestrogen receptor-positive disease (OR 0.63, 95% CI 0.49-0.80; P heterogeneity=0.01) and is also associated with mammographic breast density, a strong risk factor for breast cancer (P=0.001). rs140068132 is located within several transcription factor-binding sites and electrophoretic mobility shift assays with MCF-7 nuclear protein demonstrate differential binding of the G/A alleles at this locus. These results highlight the importance of conducting research in diverse populations.