Combined genetic and splicing analysis of BRCA1 c.[594-2A>C; 641A>G] highlights the relevance of naturally occurring in-frame transcripts for developing disease gene variant classification algorithms

A recent analysis using family history weighting and co-observation classification modeling indicated that BRCA1 c.594-2A > C (IVS9-2A > C), previously described to cause exon 10 skipping (a truncating alteration), displays characteristics inconsistent with those of a high risk pathogenic BRCA1 variant. We used large-scale genetic and clinical resources from the ENIGMA, CIMBA and BCAC consortia to assess pathogenicity of c.594-2A > C. The combined odds for causality considering case-control, segregation and breast tumor pathology information was 3.23 × 10-8 Our data indicate that c.594-2A > C is always in cis with c.641A > G. The spliceogenic effect of c.[594-2A > C;641A > G] was characterized using RNA analysis of human samples and splicing minigenes. As expected, c.[594-2A > C; 641A > G] caused exon 10 skipping, albeit not due to c.594-2A > C impairing the acceptor site but rather by c.641A > G modifying exon 10 splicing regulatory element(s). Multiple blood-based RNA assays indicated that the variant allele did not produce detectable levels of full-length transcripts, with a per allele BRCA1 expression profile composed of ≈70-80% truncating transcripts, and ≈20-30% of in-frame Δ9,10 transcripts predicted to encode a BRCA1 protein with tumor suppression function.We confirm that BRCA1c.[594-2A > C;641A > G] should not be considered a high-risk pathogenic variant. Importantly, results from our detailed mRNA analysis suggest that BRCA-associated cancer risk is likely not markedly increased for individuals who carry a truncating variant in BRCA1 exons 9 or 10, or any other BRCA1 allele that permits 20-30% of tumor suppressor function. More generally, our findings highlight the importance of assessing naturally occurring alternative splicing for clinical evaluation of variants in disease-causing genes.

RAD51B in Familial Breast Cancer

Common variation on 14q24.1, close to RAD51B, has been associated with breast cancer: rs999737 and rs2588809 with the risk of female breast cancer and rs1314913 with the risk of male breast cancer. The aim of this study was to investigate the role of RAD51B variants in breast cancer predisposition, particularly in the context of familial breast cancer in Finland. We sequenced the coding region of RAD51B in 168 Finnish breast cancer patients from the Helsinki region for identification of possible recurrent founder mutations. In addition, we studied the known rs999737, rs2588809, and rs1314913 SNPs and RAD51B haplotypes in 44,791 breast cancer cases and 43,583 controls from 40 studies participating in the Breast Cancer Association Consortium (BCAC) that were genotyped on a custom chip (iCOGS). We identified one putatively pathogenic missense mutation c.541C>T among the Finnish cancer patients and subsequently genotyped the mutation in additional breast cancer cases (n = 5259) and population controls (n = 3586) from Finland and Belarus. No significant association with breast cancer risk was seen in the meta-analysis of the Finnish datasets or in the large BCAC dataset. The association with previously identified risk variants rs999737, rs2588809, and rs1314913 was replicated among all breast cancer cases and also among familial cases in the BCAC dataset. The most significant association was observed for the haplotype carrying the risk-alleles of all the three SNPs both among all cases (odds ratio (OR): 1.15, 95% confidence interval (CI): 1.11-1.19, P = 8.88 x 10-16) and among familial cases (OR: 1.24, 95% CI: 1.16-1.32, P = 6.19 x 10-11), compared to the haplotype with the respective protective alleles. Our results suggest that loss-of-function mutations in RAD51B are rare, but common variation at the RAD51B region is significantly associated with familial breast cancer risk.

Association of genetic susceptibility variants for type 2 diabetes with breast cancer risk in women of European ancestry

PURPOSE

Type 2 diabetes (T2D) has been reported to be associated with an elevated risk of breast cancer. It is unclear, however, whether this association is due to shared genetic factors.

METHODS

We constructed a genetic risk score (GRS) using risk variants from 33 known independent T2D susceptibility loci and evaluated its relation to breast cancer risk using the data from two consortia, including 62,328 breast cancer patients and 83,817 controls of European ancestry. Unconditional logistic regression models were used to derive adjusted odds ratios (ORs) and 95 % confidence intervals (CIs) to measure the association of breast cancer risk with T2D GRS or T2D-associated genetic risk variants. Meta-analyses were conducted to obtain summary ORs across all studies.

RESULTS

The T2D GRS was not found to be associated with breast cancer risk, overall, by menopausal status, or for estrogen receptor positive or negative breast cancer. Three T2D associated risk variants were individually associated with breast cancer risk after adjustment for multiple comparisons using the Bonferroni method (at p < 0.001), rs9939609 (FTO) (OR 0.94, 95 % CI = 0.92-0.95, p = 4.13E-13), rs7903146 (TCF7L2) (OR 1.04, 95 % CI = 1.02-1.06, p = 1.26E-05), and rs8042680 (PRC1) (OR 0.97, 95 % CI = 0.95-0.99, p = 8.05E-04).

CONCLUSIONS

We have shown that several genetic risk variants were associated with the risk of both T2D and breast cancer. However, overall genetic susceptibility to T2D may not be related to breast cancer risk.

No evidence that protein truncating variants in BRIP1 are associated with breast cancer risk: implications for gene panel testing

BACKGROUND

BRCA1 interacting protein C-terminal helicase 1 (BRIP1) is one of the Fanconi Anaemia Complementation (FANC) group family of DNA repair proteins. Biallelic mutations in BRIP1 are responsible for FANC group J, and previous studies have also suggested that rare protein truncating variants in BRIP1 are associated with an increased risk of breast cancer. These studies have led to inclusion of BRIP1 on targeted sequencing panels for breast cancer risk prediction.

METHODS

We evaluated a truncating variant, p.Arg798Ter (rs137852986), and 10 missense variants of BRIP1, in 48 144 cases and 43 607 controls of European origin, drawn from 41 studies participating in the Breast Cancer Association Consortium (BCAC). Additionally, we sequenced the coding regions of BRIP1 in 13 213 cases and 5242 controls from the UK, 1313 cases and 1123 controls from three population-based studies as part of the Breast Cancer Family Registry, and 1853 familial cases and 2001 controls from Australia.

RESULTS

The rare truncating allele of rs137852986 was observed in 23 cases and 18 controls in Europeans in BCAC (OR 1.09, 95% CI 0.58 to 2.03, p=0.79). Truncating variants were found in the sequencing studies in 34 cases (0.21%) and 19 controls (0.23%) (combined OR 0.90, 95% CI 0.48 to 1.70, p=0.75).

CONCLUSIONS

These results suggest that truncating variants in BRIP1, and in particular p.Arg798Ter, are not associated with a substantial increase in breast cancer risk. Such observations have important implications for the reporting of results from breast cancer screening panels.

No clinical utility of KRAS variant rs61764370 for ovarian or breast cancer

OBJECTIVE

Clinical genetic testing is commercially available for rs61764370, an inherited variant residing in a KRAS 3' UTR microRNA binding site, based on suggested associations with increased ovarian and breast cancer risk as well as with survival time. However, prior studies, emphasizing particular subgroups, were relatively small. Therefore, we comprehensively evaluated ovarian and breast cancer risks as well as clinical outcome associated with rs61764370.

METHODS

Centralized genotyping and analysis were performed for 140,012 women enrolled in the Ovarian Cancer Association Consortium (15,357 ovarian cancer patients; 30,816 controls), the Breast Cancer Association Consortium (33,530 breast cancer patients; 37,640 controls), and the Consortium of Modifiers of BRCA1 and BRCA2 (14,765 BRCA1 and 7904 BRCA2 mutation carriers).

RESULTS

We found no association with risk of ovarian cancer (OR=0.99, 95% CI 0.94-1.04, p=0.74) or breast cancer (OR=0.98, 95% CI 0.94-1.01, p=0.19) and results were consistent among mutation carriers (BRCA1, ovarian cancer HR=1.09, 95% CI 0.97-1.23, p=0.14, breast cancer HR=1.04, 95% CI 0.97-1.12, p=0.27; BRCA2, ovarian cancer HR=0.89, 95% CI 0.71-1.13, p=0.34, breast cancer HR=1.06, 95% CI 0.94-1.19, p=0.35). Null results were also obtained for associations with overall survival following ovarian cancer (HR=0.94, 95% CI 0.83-1.07, p=0.38), breast cancer (HR=0.96, 95% CI 0.87-1.06, p=0.38), and all other previously-reported associations.

CONCLUSIONS

rs61764370 is not associated with risk of ovarian or breast cancer nor with clinical outcome for patients with these cancers. Therefore, genotyping this variant has no clinical utility related to the prediction or management of these cancers.

PALB2: research reaching to clinical outcomes for women with breast cancer

PALB2 has taken its place with bona fide breast cancer susceptibility genes. It is now well established that women who carry loss-of-function mutations in the PALB2 gene are at similarly elevated breast cancer risks to those who carry mutations in BRCA2. Information about PALB2 is now being used in breast cancer clinical genetics practice and is routinely included in breast cancer predisposition gene panel tests. Tens of thousands of women worldwide have now had genetic tests for PALB2 mutations in the context of breast cancer susceptibility. However, prospective data related to the clinical outcomes of PALB2 mutation carriers is lacking and very little information (beyond mutation penetrance) is available to guide current clinical management for carriers (affected and unaffected by cancer). In addition, clinical classification of the vast array of non-loss-of-function genetic variants identified in PALB2 is in its infancy. These are key areas of current research efforts and are important foundations on which to move information about PALB2 into the precision public health arena.

Genetically predicted longer telomere length is associated with increased risk of B-cell lymphoma subtypes

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.

Genes associated with histopathologic features of triple negative breast tumors predict molecular subtypes

Distinct subtypes of triple negative (TN) breast cancer have been identified by tumor expression profiling. However, little is known about the relationship between histopathologic features of TN tumors, which reflect aspects of both tumor behavior and tumor microenvironment, and molecular TN subtypes. The histopathologic features of TN tumors were assessed by central review and 593 TN tumors were subjected to whole genome expression profiling using the Illumina Whole Genome DASL array. TN molecular subtypes were defined based on gene expression data associated with histopathologic features of TN tumors. Gene expression analysis yielded signatures for four TN subtypes (basal-like, androgen receptor positive, immune, and stromal) consistent with previous studies. Expression analysis also identified genes significantly associated with the 12 histological features of TN tumors. Development of signatures using these markers of histopathological features resulted in six distinct TN subtype signatures, including an additional basal-like and stromal signature. The additional basal-like subtype was distinguished by elevated expression of cell motility and glucose metabolism genes and reduced expression of immune signaling genes, whereas the additional stromal subtype was distinguished by elevated expression of immunomodulatory pathway genes. Histopathologic features that reflect heterogeneity in tumor architecture, cell structure, and tumor microenvironment are related to TN subtype. Accounting for histopathologic features in the development of gene expression signatures, six major subtypes of TN breast cancer were identified.

Atlas of prostate cancer heritability in European and African-American men pinpoints tissue-specific regulation

Although genome-wide association studies have identified over 100 risk loci that explain ∼33% of familial risk for prostate cancer (PrCa), their functional effects on risk remain largely unknown. Here we use genotype data from 59,089 men of European and African American ancestries combined with cell-type-specific epigenetic data to build a genomic atlas of single-nucleotide polymorphism (SNP) heritability in PrCa. We find significant differences in heritability between variants in prostate-relevant epigenetic marks defined in normal versus tumour tissue as well as between tissue and cell lines. The majority of SNP heritability lies in regions marked by H3k27 acetylation in prostate adenoc7arcinoma cell line (LNCaP) or by DNaseI hypersensitive sites in cancer cell lines. We find a high degree of similarity between European and African American ancestries suggesting a similar genetic architecture from common variation underlying PrCa risk. Our findings showcase the power of integrating functional annotation with genetic data to understand the genetic basis of PrCa.

Breast cancer risk variants at 6q25 display different phenotype associations and regulate ESR1, RMND1 and CCDC170

We analyzed 3,872 common genetic variants across the ESR1 locus (encoding estrogen receptor α) in 118,816 subjects from three international consortia. We found evidence for at least five independent causal variants, each associated with different phenotype sets, including estrogen receptor (ER(+) or ER(-)) and human ERBB2 (HER2(+) or HER2(-)) tumor subtypes, mammographic density and tumor grade. The best candidate causal variants for ER(-) tumors lie in four separate enhancer elements, and their risk alleles reduce expression of ESR1, RMND1 and CCDC170, whereas the risk alleles of the strongest candidates for the remaining independent causal variant disrupt a silencer element and putatively increase ESR1 and RMND1 expression.

Meta-analysis of genome-wide association studies discovers multiple loci for chronic lymphocytic leukemia

Chronic lymphocytic leukemia (CLL) is a common lymphoid malignancy with strong heritability. To further understand the genetic susceptibility for CLL and identify common loci associated with risk, we conducted a meta-analysis of four genome-wide association studies (GWAS) composed of 3,100 cases and 7,667 controls with follow-up replication in 1,958 cases and 5,530 controls. Here we report three new loci at 3p24.1 (rs9880772, EOMES, P=2.55 × 10(-11)), 6p25.2 (rs73718779, SERPINB6, P=1.97 × 10(-8)) and 3q28 (rs9815073, LPP, P=3.62 × 10(-8)), as well as a new independent SNP at the known 2q13 locus (rs9308731, BCL2L11, P=1.00 × 10(-11)) in the combined analysis. We find suggestive evidence (P<5 × 10(-7)) for two additional new loci at 4q24 (rs10028805, BANK1, P=7.19 × 10(-8)) and 3p22.2 (rs1274963, CSRNP1, P=2.12 × 10(-7)). Pathway analyses of new and known CLL loci consistently show a strong role for apoptosis, providing further evidence for the importance of this biological pathway in CLL susceptibility.

Germline mutations in PMS2 and MLH1 in individuals with solitary loss of PMS2 expression in colorectal carcinomas from the Colon Cancer Family Registry Cohort

OBJECTIVES

Immunohistochemistry for DNA mismatch repair proteins is used to screen for Lynch syndrome in individuals with colorectal carcinoma (CRC). Although solitary loss of PMS2 expression is indicative of carrying a germline mutation in PMS2, previous studies reported MLH1 mutation in some cases. We determined the prevalence of MLH1 germline mutations in a large cohort of individuals with a CRC demonstrating solitary loss of PMS2 expression.

DESIGN

This cohort study included 88 individuals affected with a PMS2-deficient CRC from the Colon Cancer Family Registry Cohort. Germline PMS2 mutation analysis (long-range PCR and multiplex ligation-dependent probe amplification) was followed by MLH1 mutation testing (Sanger sequencing and multiplex ligation-dependent probe amplification).

RESULTS

Of the 66 individuals with complete mutation screening, we identified a pathogenic PMS2 mutation in 49 (74%), a pathogenic MLH1 mutation in 8 (12%) and a MLH1 variant of uncertain clinical significance predicted to be damaging by in silico analysis in 3 (4%); 6 (9%) carried variants likely to have no clinical significance. Missense point mutations accounted for most alterations (83%; 9/11) in MLH1. The MLH1 c.113A> G p.Asn38Ser mutation was found in 2 related individuals. One individual who carried the MLH1 intronic mutation c.677+3A>G p.Gln197Argfs*8 leading to the skipping of exon 8, developed 2 tumours, both of which retained MLH1 expression.

CONCLUSIONS

A substantial proportion of CRCs with solitary loss of PMS2 expression are associated with a deleterious MLH1 germline mutation supporting the screening for MLH1 in individuals with tumours of this immunophenotype, when no PMS2 mutation has been identified.

Genetic predisposition to ductal carcinoma in situ of the breast

BACKGROUND

Ductal carcinoma in situ (DCIS) is a non-invasive form of breast cancer. It is often associated with invasive ductal carcinoma (IDC), and is considered to be a non-obligate precursor of IDC. It is not clear to what extent these two forms of cancer share low-risk susceptibility loci, or whether there are differences in the strength of association for shared loci.

METHODS

To identify genetic polymorphisms that predispose to DCIS, we pooled data from 38 studies comprising 5,067 cases of DCIS, 24,584 cases of IDC and 37,467 controls, all genotyped using the iCOGS chip.

RESULTS

Most (67 %) of the 76 known breast cancer predisposition loci showed an association with DCIS in the same direction as previously reported for invasive breast cancer. Case-only analysis showed no evidence for differences between associations for IDC and DCIS after considering multiple testing. Analysis by estrogen receptor (ER) status confirmed that loci associated with ER positive IDC were also associated with ER positive DCIS. Analysis of DCIS by grade suggested that two independent SNPs at 11q13.3 near CCND1 were specific to low/intermediate grade DCIS (rs75915166, rs554219). These associations with grade remained after adjusting for ER status and were also found in IDC. We found no novel DCIS-specific loci at a genome wide significance level of P < 5.0x10(-8).

CONCLUSION

In conclusion, this study provides the strongest evidence to date of a shared genetic susceptibility for IDC and DCIS. Studies with larger numbers of DCIS are needed to determine if IDC or DCIS specific loci exist.

Heterogeneity of luminal breast cancer characterised by immunohistochemical expression of basal markers

BACKGROUND

Luminal A breast cancer defined as hormone receptor positive and human epidermal growth factor receptor 2 (HER2) negative is known to be heterogeneous. Previous study showed that luminal A tumours with the expression of basal markers ((cytokeratin (CK) 5 or CK5/6) or epidermal growth factor receptor (EGFR)) were associated with poorer prognosis compared with those that stained negative for basal markers. Prompted by this study, we assessed whether tumour characteristics and risk factors differed by basal marker status within luminal A tumours.

METHODS

We pooled 5040 luminal A cases defined by immunohistochemistry (4490 basal-negative ((CK5 (or CK5/6))- and EGFR-) and 550 basal-positive ((CK5 (or CK5/6+)) or EGFR+)) from eight studies participating in the Breast Cancer Association Consortium. Case-case comparison was performed using unconditional logistic regression.

RESULTS

Tumour characteristics and risk factors did not vary significantly by the expression of basal markers, although results suggested that basal-positive luminal tumours tended to be smaller and node negative, and were more common in women with a positive family history and lower body mass index.

CONCLUSIONS

Most established breast cancer risk factors were similar in basal-positive and basal-negative luminal A tumours. The non-significant but suggestive differences in tumour features and family history warrant further investigations.

Breast Cancer Risk Prediction Using Clinical Models and 77 Independent Risk-Associated SNPs for Women Aged Under 50 Years: Australian Breast Cancer Family Registry

BACKGROUND

The extent to which clinical breast cancer risk prediction models can be improved by including information on known susceptibility SNPs is not known.

METHODS

Using 750 cases and 405 controls from the population-based Australian Breast Cancer Family Registry who were younger than 50 years at diagnosis and recruitment, respectively, Caucasian and not BRCA1 or BRCA2 mutation carriers, we derived absolute 5-year risks of breast cancer using the BOADICEA, BRCAPRO, BCRAT, and IBIS risk prediction models and combined these with a risk score based on 77 independent risk-associated SNPs. We used logistic regression to estimate the OR per adjusted SD for log-transformed age-adjusted 5-year risks. Discrimination was assessed by the area under the receiver operating characteristic curve (AUC). Calibration was assessed using the Hosmer-Lemeshow goodness-of-fit test. We also constructed reclassification tables and calculated the net reclassification improvement.

RESULTS

The ORs for BOADICEA, BRCAPRO, BCRAT, and IBIS were 1.80, 1.75, 1.67, and 1.30, respectively. When combined with the SNP-based score, the corresponding ORs were 1.96, 1.89, 1.80, and 1.52. The corresponding AUCs were 0.66, 0.65, 0.64, and 0.57 for the risk prediction models, and 0.70, 0.69, 0.66, and 0.63 when combined with the SNP-based score.

CONCLUSIONS

By combining a 77 SNP-based score with clinical models, the AUC for predicting breast cancer before age 50 years improved by >20%.

IMPACT

Our estimates of the increased performance of clinical risk prediction models from including genetic information could be used to inform targeted screening and prevention.

BRCA2 Polymorphic Stop Codon K3326X and the Risk of Breast, Prostate, and Ovarian Cancers

BACKGROUND

The K3326X variant in BRCA2 (BRCA2*c.9976A>T; p.Lys3326*; rs11571833) has been found to be associated with small increased risks of breast cancer. However, it is not clear to what extent linkage disequilibrium with fully pathogenic mutations might account for this association. There is scant information about the effect of K3326X in other hormone-related cancers.

METHODS

Using weighted logistic regression, we analyzed data from the large iCOGS study including 76 637 cancer case patients and 83 796 control patients to estimate odds ratios (ORw) and 95% confidence intervals (CIs) for K3326X variant carriers in relation to breast, ovarian, and prostate cancer risks, with weights defined as probability of not having a pathogenic BRCA2 variant. Using Cox proportional hazards modeling, we also examined the associations of K3326X with breast and ovarian cancer risks among 7183 BRCA1 variant carriers. All statistical tests were two-sided.

RESULTS

The K3326X variant was associated with breast (ORw = 1.28, 95% CI = 1.17 to 1.40, P = 5.9x10(-) (6)) and invasive ovarian cancer (ORw = 1.26, 95% CI = 1.10 to 1.43, P = 3.8x10(-3)). These associations were stronger for serous ovarian cancer and for estrogen receptor-negative breast cancer (ORw = 1.46, 95% CI = 1.2 to 1.70, P = 3.4x10(-5) and ORw = 1.50, 95% CI = 1.28 to 1.76, P = 4.1x10(-5), respectively). For BRCA1 mutation carriers, there was a statistically significant inverse association of the K3326X variant with risk of ovarian cancer (HR = 0.43, 95% CI = 0.22 to 0.84, P = .013) but no association with breast cancer. No association with prostate cancer was observed.

CONCLUSIONS

Our study provides evidence that the K3326X variant is associated with risk of developing breast and ovarian cancers independent of other pathogenic variants in BRCA2. Further studies are needed to determine the biological mechanism of action responsible for these associations.

Strategies for Integrated Analysis of Genetic, Epigenetic, and Gene Expression Variation in Cancer: Addressing the Challenges

The development and progression of cancer, a collection of diseases with complex genetic architectures, is facilitated by the interplay of multiple etiological factors. This complexity challenges the traditional single-platform study design and calls for an integrated approach to data analysis. However, integration of heterogeneous measurements of biological variation is a non-trivial exercise due to the diversity of the human genome and the variety of output data formats and genome coverage obtained from the commonly used molecular platforms. This review article will provide an introduction to integration strategies used for analyzing genetic risk factors for cancer. We critically examine the ability of these strategies to handle the complexity of the human genome and also accommodate information about the biological and functional interactions between the elements that have been measured-making the assessment of disease risk against a composite genomic factor possible. The focus of this review is to provide an overview and introduction to the main strategies and to discuss where there is a need for further development.

Genetic variants within the hTERT gene and the risk of colorectal cancer in Lynch syndrome

Lynch syndrome is an inherited cancer-predisposing disorder caused by germline mutations in the DNA mismatch repair (MMR) genes but there is a high degree of variability in cancer risk observed among carriers, suggesting the existence of modifying factors. Our aim was to investigate variants within the hTERT gene as a potential colorectal cancer (CRC) risk modifier for MMR gene mutation carriers. We identified 1098 MMR gene mutation carriers (420 MLH1, 481 MSH2, 126 MSH6, 53 PMS2 and 18 EPCAM) from 330 families recruited from either family cancer clinics or population cancer registries of the Australasian Colorectal Cancer Family Registry between 1997 and 2012. Using weighted Cox regression after adjusting for ascertainment bias, we estimated associations between 23 SNPs within the hTERT gene and CRC risk. During 46,836 person-years observation, 392 (36%) carriers were diagnosed with CRC at a mean age of 42.2 (standard deviation 11.4) years. There was no evidence of association between any of the hTERT SNPs and CRC risk, overall and stratified by sex and MMR gene mutated, after adjustment for multiple testing. Our findings suggest no evidence for clinical utility of the SNPs within the hTERT gene in Lynch syndrome.

Genetic variation in the immunosuppression pathway genes and breast cancer susceptibility: a pooled analysis of 42,510 cases and 40,577 controls from the Breast Cancer Association Consortium

Immunosuppression plays a pivotal role in assisting tumors to evade immune destruction and promoting tumor development. We hypothesized that genetic variation in the immunosuppression pathway genes may be implicated in breast cancer tumorigenesis. We included 42,510 female breast cancer cases and 40,577 controls of European ancestry from 37 studies in the Breast Cancer Association Consortium (2015) with available genotype data for 3595 single nucleotide polymorphisms (SNPs) in 133 candidate genes. Associations between genotyped SNPs and overall breast cancer risk, and secondarily according to estrogen receptor (ER) status, were assessed using multiple logistic regression models. Gene-level associations were assessed based on principal component analysis. Gene expression analyses were conducted using RNA sequencing level 3 data from The Cancer Genome Atlas for 989 breast tumor samples and 113 matched normal tissue samples. SNP rs1905339 (A>G) in the STAT3 region was associated with an increased breast cancer risk (per allele odds ratio 1.05, 95 % confidence interval 1.03-1.08; p value = 1.4 × 10(-6)). The association did not differ significantly by ER status. On the gene level, in addition to TGFBR2 and CCND1, IL5 and GM-CSF showed the strongest associations with overall breast cancer risk (p value = 1.0 × 10(-3) and 7.0 × 10(-3), respectively). Furthermore, STAT3 and IL5 but not GM-CSF were differentially expressed between breast tumor tissue and normal tissue (p value = 2.5 × 10(-3), 4.5 × 10(-4) and 0.63, respectively). Our data provide evidence that the immunosuppression pathway genes STAT3, IL5, and GM-CSF may be novel susceptibility loci for breast cancer in women of European ancestry.

Global measures of peripheral blood-derived DNA methylation as a risk factor in the development of mature B-cell neoplasms

AIM

To examine whether peripheral blood methylation is associated with risk of developing mature B-cell neoplasms (MBCNs).

MATERIALS & METHODS

We conducted a case-control study nested within a large prospective cohort. Peripheral blood was collected from healthy participants. Cases of MBCN were identified by linkage to cancer registries. Methylation was measured using the Infinium(®) HumanMethylation450.

RESULTS

During a median of 10.6-year follow-up, 438 MBCN cases were evaluated. Global hypomethylation was associated with increased risk of MBCN (odds ratio: 2.27, [95% CI: 1.59-3.25]). Within high CpG promoter regions, hypermethylation was associated with increased risk (odds ratio: 1.76 [95% CI: 1.25-2.48]). Promoter hypermethylation was observed in HOXA9 and CDH1 genes.

CONCLUSION

Aberrant global DNA methylation is detectable in peripheral blood collected years before diagnosis and is associated with increased risk of MBCN, suggesting changes to DNA methylation are an early event in MBCN development.

Hypomethylation of smoking-related genes is associated with future lung cancer in four prospective cohorts

DNA hypomethylation in certain genes is associated with tobacco exposure but it is unknown whether these methylation changes translate into increased lung cancer risk. In an epigenome-wide study of DNA from pre-diagnostic blood samples from 132 case-control pairs in the NOWAC cohort, we observe that the most significant associations with lung cancer risk are for cg05575921 in AHRR (OR for 1 s.d.=0.37, 95% CI: 0.31-0.54, P-value=3.3 × 10(-11)) and cg03636183 in F2RL3 (OR for 1 s.d.=0.40, 95% CI: 0.31-0.56, P-value=3.9 × 10(-10)), previously shown to be strongly hypomethylated in smokers. These associations remain significant after adjustment for smoking and are confirmed in additional 664 case-control pairs tightly matched for smoking from the MCCS, NSHDS and EPIC HD cohorts. The replication and mediation analyses suggest that residual confounding is unlikely to explain the observed associations and that hypomethylation of these CpG sites may mediate the effect of tobacco on lung cancer risk.

Analysis of Heritability and Shared Heritability Based on Genome-Wide Association Studies for Thirteen Cancer Types

BACKGROUND

Studies of related individuals have consistently demonstrated notable familial aggregation of cancer. We aim to estimate the heritability and genetic correlation attributable to the additive effects of common single-nucleotide polymorphisms (SNPs) for cancer at 13 anatomical sites.

METHODS

Between 2007 and 2014, the US National Cancer Institute has generated data from genome-wide association studies (GWAS) for 49 492 cancer case patients and 34 131 control patients. We apply novel mixed model methodology (GCTA) to this GWAS data to estimate the heritability of individual cancers, as well as the proportion of heritability attributable to cigarette smoking in smoking-related cancers, and the genetic correlation between pairs of cancers.

RESULTS

GWAS heritability was statistically significant at nearly all sites, with the estimates of array-based heritability, hl (2), on the liability threshold (LT) scale ranging from 0.05 to 0.38. Estimating the combined heritability of multiple smoking characteristics, we calculate that at least 24% (95% confidence interval [CI] = 14% to 37%) and 7% (95% CI = 4% to 11%) of the heritability for lung and bladder cancer, respectively, can be attributed to genetic determinants of smoking. Most pairs of cancers studied did not show evidence of strong genetic correlation. We found only four pairs of cancers with marginally statistically significant correlations, specifically kidney and testes (ρ = 0.73, SE = 0.28), diffuse large B-cell lymphoma (DLBCL) and pediatric osteosarcoma (ρ = 0.53, SE = 0.21), DLBCL and chronic lymphocytic leukemia (CLL) (ρ = 0.51, SE =0.18), and bladder and lung (ρ = 0.35, SE = 0.14). Correlation analysis also indicates that the genetic architecture of lung cancer differs between a smoking population of European ancestry and a nonsmoking Asian population, allowing for the possibility that the genetic etiology for the same disease can vary by population and environmental exposures.

CONCLUSION

Our results provide important insights into the genetic architecture of cancers and suggest new avenues for investigation.

Epithelial-Mesenchymal Transition (EMT) Gene Variants and Epithelial Ovarian Cancer (EOC) Risk

Epithelial-mesenchymal transition (EMT) is a process whereby epithelial cells assume mesenchymal characteristics to facilitate cancer metastasis. However, EMT also contributes to the initiation and development of primary tumors. Prior studies that explored the hypothesis that EMT gene variants contribute to epithelial ovarian carcinoma (EOC) risk have been based on small sample sizes and none have sought replication in an independent population. We screened 15,816 single-nucleotide polymorphisms (SNPs) in 296 genes in a discovery phase using data from a genome-wide association study of EOC among women of European ancestry (1,947 cases and 2,009 controls) and identified 793 variants in 278 EMT-related genes that were nominally (P < 0.05) associated with invasive EOC. These SNPs were then genotyped in a larger study of 14,525 invasive-cancer patients and 23,447 controls. A P-value <0.05 and a false discovery rate (FDR) <0.2 were considered statistically significant. In the larger dataset, GPC6/GPC5 rs17702471 was associated with the endometrioid subtype among Caucasians (odds ratio (OR) = 1.16, 95% CI = 1.07-1.25, P = 0.0003, FDR = 0.19), whereas F8 rs7053448 (OR = 1.69, 95% CI = 1.27-2.24, P = 0.0003, FDR = 0.12), F8 rs7058826 (OR = 1.69, 95% CI = 1.27-2.24, P = 0.0003, FDR = 0.12), and CAPN13 rs1983383 (OR = 0.79, 95% CI = 0.69-0.90, P = 0.0005, FDR = 0.12) were associated with combined invasive EOC among Asians. In silico functional analyses revealed that GPC6/GPC5 rs17702471 coincided with DNA regulatory elements. These results suggest that EMT gene variants do not appear to play a significant role in the susceptibility to EOC.

SNP-SNP interaction analysis of NF-κB signaling pathway on breast cancer survival

In breast cancer, constitutive activation of NF-κB has been reported, however, the impact of genetic variation of the pathway on patient prognosis has been little studied. Furthermore, a combination of genetic variants, rather than single polymorphisms, may affect disease prognosis. Here, in an extensive dataset (n = 30,431) from the Breast Cancer Association Consortium, we investigated the association of 917 SNPs in 75 genes in the NF-κB pathway with breast cancer prognosis. We explored SNP-SNP interactions on survival using the likelihood-ratio test comparing multivariate Cox' regression models of SNP pairs without and with an interaction term. We found two interacting pairs associating with prognosis: patients simultaneously homozygous for the rare alleles of rs5996080 and rs7973914 had worse survival (HRinteraction 6.98, 95% CI=3.3-14.4, P=1.42E-07), and patients carrying at least one rare allele for rs17243893 and rs57890595 had better survival (HRinteraction 0.51, 95% CI=0.3-0.6, P = 2.19E-05). Based on in silico functional analyses and literature, we speculate that the rs5996080 and rs7973914 loci may affect the BAFFR and TNFR1/TNFR3 receptors and breast cancer survival, possibly by disturbing both the canonical and non-canonical NF-κB pathways or their dynamics, whereas, rs17243893-rs57890595 interaction on survival may be mediated through TRAF2-TRAIL-R4 interplay. These results warrant further validation and functional analyses.

Height and Breast Cancer Risk: Evidence From Prospective Studies and Mendelian Randomization

BACKGROUND

Epidemiological studies have linked adult height with breast cancer risk in women. However, the magnitude of the association, particularly by subtypes of breast cancer, has not been established. Furthermore, the mechanisms of the association remain unclear.

METHODS

We performed a meta-analysis to investigate associations between height and breast cancer risk using data from 159 prospective cohorts totaling 5216302 women, including 113178 events. In a consortium with individual-level data from 46325 case patients and 42482 control patients, we conducted a Mendelian randomization analysis using a genetic score that comprised 168 height-associated variants as an instrument. This association was further evaluated in a second consortium using summary statistics data from 16003 case patients and 41335 control patients.

RESULTS

The pooled relative risk of breast cancer was 1.17 (95% confidence interval [CI] = 1.15 to 1.19) per 10cm increase in height in the meta-analysis of prospective studies. In Mendelian randomization analysis, the odds ratio of breast cancer per 10cm increase in genetically predicted height was 1.22 (95% CI = 1.13 to 1.32) in the first consortium and 1.21 (95% CI = 1.05 to 1.39) in the second consortium. The association was found in both premenopausal and postmenopausal women but restricted to hormone receptor-positive breast cancer. Analyses of height-associated variants identified eight new loci associated with breast cancer risk after adjusting for multiple comparisons, including three loci at 1q21.2, DNAJC27, and CCDC91 at genome-wide significance level P < 5×10(-8).

CONCLUSIONS

Our study provides strong evidence that adult height is a risk factor for breast cancer in women and certain genetic factors and biological pathways affecting adult height have an important role in the etiology of breast cancer.