Many individuals with multiple or large colorectal adenomas or early-onset colorectal cancer (CRC) have no detectable germline mutations in the known cancer predisposition genes. Using whole-genome sequencing, supplemented by linkage and association analysis, we identified specific heterozygous POLE or POLD1 germline variants in several multiple-adenoma and/or CRC cases but in no controls. The variants associated with susceptibility, POLE p.Leu424Val and POLD1 p.Ser478Asn, have high penetrance, and POLD1 mutation was also associated with endometrial cancer predisposition. The mutations map to equivalent sites in the proofreading (exonuclease) domain of DNA polymerases ɛ and δ and are predicted to cause a defect in the correction of mispaired bases inserted during DNA replication. In agreement with this prediction, the tumors from mutation carriers were microsatellite stable but tended to acquire base substitution mutations, as confirmed by yeast functional assays. Further analysis of published data showed that the recently described group of hypermutant, microsatellite-stable CRCs is likely to be caused by somatic POLE mutations affecting the exonuclease domain.

The contribution of germline copy number variants (CNVs) to risk of developing cancer in individuals with pathogenic BRCA1 or BRCA2 variants remains relatively unknown. We conducted the largest genome-wide analysis of CNVs in 15,342 BRCA1 and 10,740 BRCA2 pathogenic variant carriers. We used these results to prioritise a candidate breast cancer risk-modifier gene for laboratory analysis and biological validation. Notably, the HR for deletions in BRCA1 suggested an elevated breast cancer risk estimate (hazard ratio (HR) = 1.21), 95% confidence interval (95% CI = 1.09-1.35) compared with non-CNV pathogenic variants. In contrast, deletions overlapping SULT1A1 suggested a decreased breast cancer risk (HR = 0.73, 95% CI 0.59-0.91) in BRCA1 pathogenic variant carriers. Functional analyses of SULT1A1 showed that reduced mRNA expression in pathogenic BRCA1 variant cells was associated with reduced cellular proliferation and reduced DNA damage after treatment with DNA damaging agents. These data provide evidence that deleterious variants in BRCA1 plus SULT1A1 deletions contribute to variable breast cancer risk in BRCA1 carriers.

Melanoma, the most aggressive skin cancer type, is responsible for 75% of skin cancer related deaths worldwide. Given that New Zealand (NZ) has the world's highest melanoma incidence, we sought to determine the frequency of mutations in NZ melanomas in recurrently mutated genes. NZ melanomas were from localities distributed between North (35°S-42°S) and South Islands (41°S-47°S). A total of 529 melanomas were analyzed for BRAF exon 15 mutations by Sanger sequencing, and also by Sequenom MelaCarta MassARRAY. While, a relatively low incidence of BRAFV600E mutations (23.4%) was observed overall in NZ melanomas, the incidence of NRAS mutations in South Island melanomas was high compared to North Island melanomas (38.3% vs. 21.9%, P=0.0005), and to The Cancer Genome Atlas database (TCGA) (38.3% vs. 22%, P=0.0004). In contrast, the incidence of EPHB6G404S mutations was 0% in South Island melanomas, and was 7.8% in North Island (P=0.0002). Overall, these data suggest that melanomas from geographically different regions in NZ have markedly different mutation frequencies, in particular in the NRAS and EPHB6 genes, when compared to TCGA or other populations. These data have implications for the causation and treatment of malignant melanoma in NZ.

In genome-wide association studies (GWASs) of colorectal cancer, we have identified two genomic regions in which pairs of tagging-single nucleotide polymorphisms (tagSNPs) are associated with disease; these comprise chromosomes 1q41 (rs6691170, rs6687758) and 12q13.13 (rs7163702, rs11169552). We investigated these regions further, aiming to determine whether they contain more than one independent association signal and/or to identify the SNPs most strongly associated with disease. Genotyping of additional sample sets at the original tagSNPs showed that, for both regions, the two tagSNPs were unlikely to identify a single haplotype on which the functional variation lay. Conversely, one of the pair of SNPs did not fully capture the association signal in each region. We therefore undertook more detailed analyses, using imputation, logistic regression, genealogical analysis using the GENECLUSTER program and haplotype analysis. In the 1q41 region, the SNP rs11118883 emerged as a strong candidate based on all these analyses, sufficient to account for the signals at both rs6691170 and rs6687758. rs11118883 lies within a region with strong evidence of transcriptional regulatory activity and has been associated with expression of PDGFRB mRNA. For 12q13.13, a complex situation was found: SNP rs7972465 showed stronger association than either rs11169552 or rs7136702, and GENECLUSTER found no good evidence for a two-SNP model. However, logistic regression and haplotype analyses supported a two-SNP model, in which a signal at the SNP rs706793 was added to that at rs11169552. Post-GWAS fine-mapping studies are challenging, but the use of multiple tools can assist in identifying candidate functional variants in at least some cases.

Five single nucleotide polymorphisms (SNPs) associated with thyroid cancer (TC) risk have been reported: rs2910164 (5q24); rs6983267 (8q24); rs965513 and rs1867277 (9q22); and rs944289 (14q13). Most of these associations have not been replicated in independent populations and the combined effects of the SNPs on risk have not been examined. This study genotyped the five TC SNPs in 781 patients recruited through the TCUKIN study. Genotype data from 6122 controls were obtained from the CORGI and Wellcome Trust Case-Control Consortium studies. Significant associations were detected between TC and rs965513A (p=6.35×10(-34)), rs1867277A (p=5.90×10(-24)), rs944289T (p=6.95×10(-7)), and rs6983267G (p=0.016). rs6983267 was most strongly associated under a recessive model (P(GG vs GT + TT)=0.004), in contrast to the association of this SNP with other cancer types. However, no evidence was found of an association between rs2910164 and disease under any risk model (p>0.7). The rs1867277 association remained significant (p=0.008) after accounting for genotypes at the nearby rs965513 (p=2.3×10(-13)) and these SNPs did not tag a single high risk haplotype. The four validated TC SNPs accounted for a relatively large proportion (∼11%) of the sibling relative risk of TC, principally owing to the large effect size of rs965513 (OR 1.74).

Increasing BRCA1 and BRCA2 (collectively termed herein as BRCA) gene testing is required to improve cancer management and prevent BRCA-related cancers.

Genome-wide association studies (GWAS) have identified 14 tagging single nucleotide polymorphisms (tagSNPs) that are associated with the risk of colorectal cancer (CRC), and several of these tagSNPs are near bone morphogenetic protein (BMP) pathway loci. The penalty of multiple testing implicit in GWAS increases the attraction of complementary approaches for disease gene discovery, including candidate gene- or pathway-based analyses. The strongest candidate loci for additional predisposition SNPs are arguably those already known both to have functional relevance and to be involved in disease risk. To investigate this proposition, we searched for novel CRC susceptibility variants close to the BMP pathway genes GREM1 (15q13.3), BMP4 (14q22.2), and BMP2 (20p12.3) using sample sets totalling 24,910 CRC cases and 26,275 controls. We identified new, independent CRC predisposition SNPs close to BMP4 (rs1957636, P = 3.93×10(-10)) and BMP2 (rs4813802, P = 4.65×10(-11)). Near GREM1, we found using fine-mapping that the previously-identified association between tagSNP rs4779584 and CRC actually resulted from two independent signals represented by rs16969681 (P = 5.33×10(-8)) and rs11632715 (P = 2.30×10(-10)). As low-penetrance predisposition variants become harder to identify-owing to small effect sizes and/or low risk allele frequencies-approaches based on informed candidate gene selection may become increasingly attractive. Our data emphasise that genetic fine-mapping studies can deconvolute associations that have arisen owing to independent correlation of a tagSNP with more than one functional SNP, thus explaining some of the apparently missing heritability of common diseases.

Genome-wide association studies (GWAS) have identified ten loci harboring common variants that influence risk of developing colorectal cancer (CRC). To enhance the power to identify additional CRC risk loci, we conducted a meta-analysis of three GWAS from the UK which included a total of 3,334 affected individuals (cases) and 4,628 controls followed by multiple validation analyses including a total of 18,095 cases and 20,197 controls. We identified associations at four new CRC risk loci: 1q41 (rs6691170, odds ratio (OR) = 1.06, P = 9.55 × 10⁻¹⁰ and rs6687758, OR = 1.09, P = 2.27 × 10⁻⁹, 3q26.2 (rs10936599, OR = 0.93, P = 3.39 × 10⁻⁸), 12q13.13 (rs11169552, OR = 0.92, P = 1.89 × 10⁻¹⁰ and rs7136702, OR = 1.06, P = 4.02 × 10⁻⁸) and 20q13.33 (rs4925386, OR = 0.93, P = 1.89 × 10⁻¹⁰). In addition to identifying new CRC risk loci, this analysis provides evidence that additional CRC-associated variants of similar effect size remain to be discovered.