CpG-biomarkers in tumor tissue and prediction models for the survival of colorectal cancer: A systematic review and external validation study

The research aimed to identify previously published CpG-methylation-based prognostic biomarkers and prediction models for colorectal cancer (CRC) prognosis and validate them in a large external cohort. A systematic search was conducted, analyzing 298 unique CpGs and 12 CpG-based prognostic models from 28 studies. After adjustment for clinical variables, 48 CpGs and five prognostic models were confirmed to be associated with survival. However, the discrimination ability of the models was insufficient, with area under the receiver operating characteristic curves ranging from 0.53 to 0.62. Calibration accuracy was mostly poor, and no significant added prognostic value beyond traditional clinical variables was observed. All prognostic models were rated at high risk of bias. While a fraction of CpGs showed potential clinical utility and generalizability, the CpG-based prognostic models performed poorly and lacked clinical relevance.

FANCM missense variants and breast cancer risk: a case-control association study of 75,156 European women

Evidence from literature, including the BRIDGES study, indicates that germline protein truncating variants (PTVs) in FANCM confer moderately increased risk of ER-negative and triple-negative breast cancer (TNBC), especially for women with a family history of the disease. Association between FANCM missense variants (MVs) and breast cancer risk has been postulated. In this study, we further used the BRIDGES study to test 689 FANCM MVs for association with breast cancer risk, overall and in ER-negative and TNBC subtypes, in 39,885 cases (7566 selected for family history) and 35,271 controls of European ancestry. Sixteen common MVs were tested individually; the remaining rare 673 MVs were tested by burden analyses considering their position and pathogenicity score. We also conducted a meta-analysis of our results and those from published studies. We did not find evidence for association for any of the 16 variants individually tested. The rare MVs were significantly associated with increased risk of ER-negative breast cancer by burden analysis comparing familial cases to controls (OR = 1.48; 95% CI 1.07-2.04; P = 0.017). Higher ORs were found for the subgroup of MVs located in functional domains or predicted to be pathogenic. The meta-analysis indicated that FANCM MVs overall are associated with breast cancer risk (OR = 1.22; 95% CI 1.08-1.38; P = 0.002). Our results support the definition from previous analyses of FANCM as a moderate-risk breast cancer gene and provide evidence that FANCM MVs could be low/moderate risk factors for ER-negative and TNBC subtypes. Further genetic and functional analyses are necessary to clarify better the increased risks due to FANCM MVs.

Incorporating progesterone receptor expression into the PREDICT breast prognostic model

BACKGROUND

Predict Breast (www.predict.nhs.uk) is an online prognostication and treatment benefit tool for early invasive breast cancer. The aim of this study was to incorporate the prognostic effect of progesterone receptor (PR) status into a new version of PREDICT and to compare its performance to the current version (2.2).

METHOD

The prognostic effect of PR status was based on the analysis of data from 45,088 European patients with breast cancer from 49 studies in the Breast Cancer Association Consortium. Cox proportional hazard models were used to estimate the hazard ratio for PR status. Data from a New Zealand study of 11,365 patients with early invasive breast cancer were used for external validation. Model calibration and discrimination were used to test the model performance.

RESULTS

Having a PR-positive tumour was associated with a 23% and 28% lower risk of dying from breast cancer for women with oestrogen receptor (ER)-negative and ER-positive breast cancer, respectively. The area under the ROC curve increased with the addition of PR status from 0.807 to 0.809 for patients with ER-negative tumours (p = 0.023) and from 0.898 to 0.902 for patients with ER-positive tumours (p = 2.3 × 10-6) in the New Zealand cohort. Model calibration was modest with 940 observed deaths compared to 1151 predicted.

CONCLUSION

The inclusion of the prognostic effect of PR status to PREDICT Breast has led to an improvement of model performance and more accurate absolute treatment benefit predictions for individual patients. Further studies should determine whether the baseline hazard function requires recalibration.

Common variants in breast cancer risk loci predispose to distinct tumor subtypes

BACKGROUND

Genome-wide association studies (GWAS) have identified multiple common breast cancer susceptibility variants. Many of these variants have differential associations by estrogen receptor (ER) status, but how these variants relate with other tumor features and intrinsic molecular subtypes is unclear.

METHODS

Among 106,571 invasive breast cancer cases and 95,762 controls of European ancestry with data on 173 breast cancer variants identified in previous GWAS, we used novel two-stage polytomous logistic regression models to evaluate variants in relation to multiple tumor features (ER, progesterone receptor (PR), human epidermal growth factor receptor 2 (HER2) and grade) adjusting for each other, and to intrinsic-like subtypes.

RESULTS

Eighty-five of 173 variants were associated with at least one tumor feature (false discovery rate < 5%), most commonly ER and grade, followed by PR and HER2. Models for intrinsic-like subtypes found nearly all of these variants (83 of 85) associated at p < 0.05 with risk for at least one luminal-like subtype, and approximately half (41 of 85) of the variants were associated with risk of at least one non-luminal subtype, including 32 variants associated with triple-negative (TN) disease. Ten variants were associated with risk of all subtypes in different magnitude. Five variants were associated with risk of luminal A-like and TN subtypes in opposite directions.

CONCLUSION

This report demonstrates a high level of complexity in the etiology heterogeneity of breast cancer susceptibility variants and can inform investigations of subtype-specific risk prediction.

Mendelian randomisation study of smoking exposure in relation to breast cancer risk

BACKGROUND

Despite a modest association between tobacco smoking and breast cancer risk reported by recent epidemiological studies, it is still equivocal whether smoking is causally related to breast cancer risk.

METHODS

We applied Mendelian randomisation (MR) to evaluate a potential causal effect of cigarette smoking on breast cancer risk. Both individual-level data as well as summary statistics for 164 single-nucleotide polymorphisms (SNPs) reported in genome-wide association studies of lifetime smoking index (LSI) or cigarette per day (CPD) were used to obtain MR effect estimates. Data from 108,420 invasive breast cancer cases and 87,681 controls were used for the LSI analysis and for the CPD analysis conducted among ever-smokers from 26,147 cancer cases and 26,072 controls. Sensitivity analyses were conducted to address pleiotropy.

RESULTS

Genetically predicted LSI was associated with increased breast cancer risk (OR 1.18 per SD, 95% CI: 1.07-1.30, P = 0.11 × 10-2), but there was no evidence of association for genetically predicted CPD (OR 1.02, 95% CI: 0.78-1.19, P = 0.85). The sensitivity analyses yielded similar results and showed no strong evidence of pleiotropic effect.

CONCLUSION

Our MR study provides supportive evidence for a potential causal association with breast cancer risk for lifetime smoking exposure but not cigarettes per day among smokers.

Author Correction: A case-only study to identify genetic modifiers of breast cancer risk for BRCA1/BRCA2 mutation carriers

A Correction to this paper has been published: https://doi.org/10.1038/s41467-021-23162-4

A case-only study to identify genetic modifiers of breast cancer risk for BRCA1/BRCA2 mutation carriers

Breast cancer (BC) risk for BRCA1 and BRCA2 mutation carriers varies by genetic and familial factors. About 50 common variants have been shown to modify BC risk for mutation carriers. All but three, were identified in general population studies. Other mutation carrier-specific susceptibility variants may exist but studies of mutation carriers have so far been underpowered. We conduct a novel case-only genome-wide association study comparing genotype frequencies between 60,212 general population BC cases and 13,007 cases with BRCA1 or BRCA2 mutations. We identify robust novel associations for 2 variants with BC for BRCA1 and 3 for BRCA2 mutation carriers, P < 10-8, at 5 loci, which are not associated with risk in the general population. They include rs60882887 at 11p11.2 where MADD, SP11 and EIF1, genes previously implicated in BC biology, are predicted as potential targets. These findings will contribute towards customising BC polygenic risk scores for BRCA1 and BRCA2 mutation carriers.

Breast Cancer Risk Genes - Association Analysis in More than 113,000 Women

BACKGROUND

Genetic testing for breast cancer susceptibility is widely used, but for many genes, evidence of an association with breast cancer is weak, underlying risk estimates are imprecise, and reliable subtype-specific risk estimates are lacking.

METHODS

We used a panel of 34 putative susceptibility genes to perform sequencing on samples from 60,466 women with breast cancer and 53,461 controls. In separate analyses for protein-truncating variants and rare missense variants in these genes, we estimated odds ratios for breast cancer overall and tumor subtypes. We evaluated missense-variant associations according to domain and classification of pathogenicity.

RESULTS

Protein-truncating variants in 5 genes (ATM, BRCA1, BRCA2, CHEK2, and PALB2) were associated with a risk of breast cancer overall with a P value of less than 0.0001. Protein-truncating variants in 4 other genes (BARD1, RAD51C, RAD51D, and TP53) were associated with a risk of breast cancer overall with a P value of less than 0.05 and a Bayesian false-discovery probability of less than 0.05. For protein-truncating variants in 19 of the remaining 25 genes, the upper limit of the 95% confidence interval of the odds ratio for breast cancer overall was less than 2.0. For protein-truncating variants in ATM and CHEK2, odds ratios were higher for estrogen receptor (ER)-positive disease than for ER-negative disease; for protein-truncating variants in BARD1, BRCA1, BRCA2, PALB2, RAD51C, and RAD51D, odds ratios were higher for ER-negative disease than for ER-positive disease. Rare missense variants (in aggregate) in ATM, CHEK2, and TP53 were associated with a risk of breast cancer overall with a P value of less than 0.001. For BRCA1, BRCA2, and TP53, missense variants (in aggregate) that would be classified as pathogenic according to standard criteria were associated with a risk of breast cancer overall, with the risk being similar to that of protein-truncating variants.

CONCLUSIONS

The results of this study define the genes that are most clinically useful for inclusion on panels for the prediction of breast cancer risk, as well as provide estimates of the risks associated with protein-truncating variants, to guide genetic counseling. (Funded by European Union Horizon 2020 programs and others.).

CYP3A7*1C allele: linking premenopausal oestrone and progesterone levels with risk of hormone receptor-positive breast cancers

BACKGROUND

Epidemiological studies provide strong evidence for a role of endogenous sex hormones in the aetiology of breast cancer. The aim of this analysis was to identify genetic variants that are associated with urinary sex-hormone levels and breast cancer risk.

METHODS

We carried out a genome-wide association study of urinary oestrone-3-glucuronide and pregnanediol-3-glucuronide levels in 560 premenopausal women, with additional analysis of progesterone levels in 298 premenopausal women. To test for the association with breast cancer risk, we carried out follow-up genotyping in 90,916 cases and 89,893 controls from the Breast Cancer Association Consortium. All women were of European ancestry.

RESULTS

For pregnanediol-3-glucuronide, there were no genome-wide significant associations; for oestrone-3-glucuronide, we identified a single peak mapping to the CYP3A locus, annotated by rs45446698. The minor rs45446698-C allele was associated with lower oestrone-3-glucuronide (-49.2%, 95% CI -56.1% to -41.1%, P = 3.1 × 10-18); in follow-up analyses, rs45446698-C was also associated with lower progesterone (-26.7%, 95% CI -39.4% to -11.6%, P = 0.001) and reduced risk of oestrogen and progesterone receptor-positive breast cancer (OR = 0.86, 95% CI 0.82-0.91, P = 6.9 × 10-8).

CONCLUSIONS

The CYP3A7*1C allele is associated with reduced risk of hormone receptor-positive breast cancer possibly mediated via an effect on the metabolism of endogenous sex hormones in premenopausal women.

Circulating cell-free DNA variables as marker of ovarian cancer patients: A pilot study

BACKGROUND

Minimal invasive blood-based molecular markers are evaluated as promising biomarkers in malignant diseases these days.

OBJECTIVE

In this pilot study, we investigated the potential of cell-free DNA (cfDNA) concentration and cell-free DNA Integrity (cfDI) as blood-based diagnostic markers for ovarian cancer patients in a retrospective study cohort.

METHODS

cfDNA concentration and cfDI were determined in the plasma of 37 ovarian cancer patients and 28 healthy controls, by measuring ALU and LINE1 repetitive DNA elements using quantitative real-time PCR.

RESULTS

A high correlation was observed between the results of ALU and LINE1. The correlated co-efficiency between the values of cfDNA concentration and cfDI was 0.86 and 0.71. As for the results between cases and controls, no or just borderline significant difference was observed in cfDI after age adjustment (P= 0.40 for ALU and P= 0.05 for LINE1) while cfDNA concentration showed a significant difference between ovarian cancer patients and healthy controls groups (P= 0.03 for ALU and P= 3.00 E-03 for LINE1). cfDNA concentration of ALU and LINE1 had an AUC of 0.81 (0.70-0.91). ALU and LINE1 cfDI reached an AUC of 0.60 (95% CI: 0.46-0.73). The combination of these markers reached the best diagnostic power with an AUC of 0.84.

CONCLUSIONS

cfDNA variables might be potentially diagnostic biomarkers in ovarian cancer, in combination with additional molecular markers. However, further studies are needed to confirm the diagnostic ability of cfDNA variables (cfDNA concentration and cfDI).

Genetic Variants in the Regulatory T cell-Related Pathway and Colorectal Cancer Prognosis

BACKGROUND

High numbers of lymphocytes in tumor tissue, including T regulatory cells (Treg), have been associated with better colorectal cancer survival. Tregs, a subset of CD4+ T lymphocytes, are mediators of immunosuppression in cancer, and therefore variants in genes related to Treg differentiation and function could be associated with colorectal cancer prognosis.

METHODS

In a prospective German cohort of 3,593 colorectal cancer patients, we assessed the association of 771 single-nucleotide polymorphisms (SNP) in 58 Treg-related genes with overall and colorectal cancer-specific survival using Cox regression models. Effect modification by microsatellite instability (MSI) status was also investigated because tumors with MSI show greater lymphocytic infiltration and have been associated with better prognosis. Replication of significant results was attempted in 2,047 colorectal cancer patients of the International Survival Analysis in Colorectal Cancer Consortium (ISACC).

RESULTS

A significant association of the TGFBR3 SNP rs7524066 with more favorable colorectal cancer-specific survival [hazard ratio (HR) per minor allele: 0.83; 95% confidence interval (CI), 0.74-0.94; P value: 0.0033] was replicated in ISACC (HR: 0.82; 95% CI, 0.68-0.98; P value: 0.03). Suggestive evidence for association was found with two IL7 SNPs, rs16906568 and rs7845577. Thirteen SNPs with differential associations with overall survival according to MSI in the discovery analysis were not confirmed.

CONCLUSIONS

Common genetic variation in the Treg pathway implicating genes such as TGFBR3 and IL7 was shown to be associated with prognosis of colorectal cancer patients.

IMPACT

The implicated genes warrant further investigation.

Breast Cancer Polygenic Risk Score and Contralateral Breast Cancer Risk

Previous research has shown that polygenic risk scores (PRSs) can be used to stratify women according to their risk of developing primary invasive breast cancer. This study aimed to evaluate the association between a recently validated PRS of 313 germline variants (PRS313) and contralateral breast cancer (CBC) risk. We included 56,068 women of European ancestry diagnosed with first invasive breast cancer from 1990 onward with follow-up from the Breast Cancer Association Consortium. Metachronous CBC risk (N = 1,027) according to the distribution of PRS313 was quantified using Cox regression analyses. We assessed PRS313 interaction with age at first diagnosis, family history, morphology, ER status, PR status, and HER2 status, and (neo)adjuvant therapy. In studies of Asian women, with limited follow-up, CBC risk associated with PRS313 was assessed using logistic regression for 340 women with CBC compared with 12,133 women with unilateral breast cancer. Higher PRS313 was associated with increased CBC risk: hazard ratio per standard deviation (SD) = 1.25 (95%CI = 1.18-1.33) for Europeans, and an OR per SD = 1.15 (95%CI = 1.02-1.29) for Asians. The absolute lifetime risks of CBC, accounting for death as competing risk, were 12.4% for European women at the 10th percentile and 20.5% at the 90th percentile of PRS313. We found no evidence of confounding by or interaction with individual characteristics, characteristics of the primary tumor, or treatment. The C-index for the PRS313 alone was 0.563 (95%CI = 0.547-0.586). In conclusion, PRS313 is an independent factor associated with CBC risk and can be incorporated into CBC risk prediction models to help improve stratification and optimize surveillance and treatment strategies.

Polymorphisms in the Angiogenesis-Related Genes EFNB2, MMP2 and JAG1 Are Associated with Survival of Colorectal Cancer Patients

An individual’s inherited genetic variation may contribute to the ‘angiogenic switch’, which is essential for blood supply and tumor growth of microscopic and macroscopic tumors. Polymorphisms in angiogenesis-related genes potentially predispose to colorectal cancer (CRC) or affect the survival of CRC patients. We investigated the association of 392 single nucleotide polymorphisms (SNPs) in 33 angiogenesis-related genes with CRC risk and survival of CRC patients in 1754 CRC cases and 1781 healthy controls within DACHS (Darmkrebs: Chancen der Verhütung durch Screening), a German population-based case-control study. Odds ratios and 95% confidence intervals (CI) were estimated from unconditional logistic regression to test for genetic associations with CRC risk. The Cox proportional hazard model was used to estimate hazard ratios (HR) and 95% CIs for survival. Multiple testing was adjusted for by a false discovery rate. No variant was associated with CRC risk. Variants in EFNB2, MMP2 and JAG1 were significantly associated with overall survival. The association of the EFNB2 tagging SNP rs9520090 (p < 0.0001) was confirmed in two validation datasets (p-values: 0.01 and 0.05). The associations of the tagging SNPs rs6040062 in JAG1 (p-value 0.0003) and rs2241145 in MMP2 (p-value 0.0005) showed the same direction of association with overall survival in the first and second validation sets, respectively, although they did not reach significance (p-values: 0.09 and 0.25, respectively). EFNB2, MMP2 and JAG1 are known for their functional role in angiogenesis and the present study points to novel evidence for the impact of angiogenesis-related genetic variants on the CRC outcome.

Transcriptome-wide association study of breast cancer risk by estrogen-receptor status

Previous transcriptome-wide association studies (TWAS) have identified breast cancer risk genes by integrating data from expression quantitative loci and genome-wide association studies (GWAS), but analyses of breast cancer subtype-specific associations have been limited. In this study, we conducted a TWAS using gene expression data from GTEx and summary statistics from the hitherto largest GWAS meta-analysis conducted for breast cancer overall, and by estrogen receptor subtypes (ER+ and ER-). We further compared associations with ER+ and ER- subtypes, using a case-only TWAS approach. We also conducted multigene conditional analyses in regions with multiple TWAS associations. Two genes, STXBP4 and HIST2H2BA, were specifically associated with ER+ but not with ER- breast cancer. We further identified 30 TWAS-significant genes associated with overall breast cancer risk, including four that were not identified in previous studies. Conditional analyses identified single independent breast-cancer gene in three of six regions harboring multiple TWAS-significant genes. Our study provides new information on breast cancer genetics and biology, particularly about genomic differences between ER+ and ER- breast cancer.

Blood-derived DNA methylation predictors of mortality discriminate tumor and healthy tissue in multiple organs

Evidence has shown that certain methylation markers derived from blood can mirror corresponding methylation signatures in internal tissues. In the current study, we aimed to investigate two strong epigenetic predictors for life span, derived from blood DNA methylation data, in tissue samples of solid cancer patients. Using data from the Cancer Genome Atlas (TCGA) and the German DACHS study, we compared a mortality risk score (MRscore) and DNAmPhenoAge in paired tumor and adjacent normal tissue samples of patients with lung (N = 69), colorectal (n = 299), breast (n = 90), head/neck (n = 50), prostate (n = 50), and liver (n = 50) cancer. To explore the concordance across tissue and blood, we additionally assessed the two markers in blood samples of colorectal cancer (CRC) cases and matched controls (n = 93) in the DACHS+ study. The MRscore was significantly elevated in tumor tissues compared to normal tissues of all cancers except prostate cancer, for which an opposite pattern was observed. DNAmPhenoAge was consistently higher in all tumor tissues. The MRscore discriminated lung, colorectal, and prostate tumor tissues from normal tissues with very high accuracy [AUCs of 0.87, 0.99 (TCGA) /0.94 (DACHS), and 0.92, respectively]. DNAmPhenoAge accurately discriminated five types of tumor tissues from normal tissues (except prostate cancer), with AUCs of 0.82–0.93. The MRscore was also significantly higher in blood samples of CRC cases than in controls, with areas under the curve (AUC) of 0.74, whereas DNAmPhenoAge did not distinguish cases from controls, with AUC of 0.54. This study provides compelling evidence that blood‐derived DNAm markers could reflect methylation changes in less accessible tissues. Further research should explore the potential use of these findings for cancer diagnosis and early detection.

Germline HOXB13 mutations p.G84E and p.R217C do not confer an increased breast cancer risk

In breast cancer, high levels of homeobox protein Hox-B13 (HOXB13) have been associated with disease progression of ER-positive breast cancer patients and resistance to tamoxifen treatment. Since HOXB13 p.G84E is a prostate cancer risk allele, we evaluated the association between HOXB13 germline mutations and breast cancer risk in a previous study consisting of 3,270 familial non-BRCA1/2 breast cancer cases and 2,327 controls from the Netherlands. Although both recurrent HOXB13 mutations p.G84E and p.R217C were not associated with breast cancer risk, the risk estimation for p.R217C was not very precise. To provide more conclusive evidence regarding the role of HOXB13 in breast cancer susceptibility, we here evaluated the association between HOXB13 mutations and increased breast cancer risk within 81 studies of the international Breast Cancer Association Consortium containing 68,521 invasive breast cancer patients and 54,865 controls. Both HOXB13 p.G84E and p.R217C did not associate with the development of breast cancer in European women, neither in the overall analysis (OR = 1.035, 95% CI = 0.859-1.246, P = 0.718 and OR = 0.798, 95% CI = 0.482-1.322, P = 0.381 respectively), nor in specific high-risk subgroups or breast cancer subtypes. Thus, although involved in breast cancer progression, HOXB13 is not a material breast cancer susceptibility gene.

Genome-wide association study identifies 32 novel breast cancer susceptibility loci from overall and subtype-specific analyses

Breast cancer susceptibility variants frequently show heterogeneity in associations by tumor subtype1-3. To identify novel loci, we performed a genome-wide association study including 133,384 breast cancer cases and 113,789 controls, plus 18,908 BRCA1 mutation carriers (9,414 with breast cancer) of European ancestry, using both standard and novel methodologies that account for underlying tumor heterogeneity by estrogen receptor, progesterone receptor and human epidermal growth factor receptor 2 status and tumor grade. We identified 32 novel susceptibility loci (P < 5.0 × 10-8), 15 of which showed evidence for associations with at least one tumor feature (false discovery rate < 0.05). Five loci showed associations (P < 0.05) in opposite directions between luminal and non-luminal subtypes. In silico analyses showed that these five loci contained cell-specific enhancers that differed between normal luminal and basal mammary cells. The genetic correlations between five intrinsic-like subtypes ranged from 0.35 to 0.80. The proportion of genome-wide chip heritability explained by all known susceptibility loci was 54.2% for luminal A-like disease and 37.6% for triple-negative disease. The odds ratios of polygenic risk scores, which included 330 variants, for the highest 1% of quantiles compared with middle quantiles were 5.63 and 3.02 for luminal A-like and triple-negative disease, respectively. These findings provide an improved understanding of genetic predisposition to breast cancer subtypes and will inform the development of subtype-specific polygenic risk scores.

A Post-Processing Algorithm for miRNA Microarray Data

One of the main disadvantages of using DNA microarrays for miRNA expression profiling is the inability of adequate comparison of expression values across different miRNAs. This leads to a large amount of miRNAs with high scores which are actually not expressed in examined samples, i.e., false positives. We propose a post-processing algorithm which performs scoring of miRNAs in the results of microarray analysis based on expression values, time of discovery of miRNA, and correlation level between the expressions of miRNA and corresponding pre-miRNA in considered samples. The algorithm was successfully validated by the comparison of the results of its application to miRNA microarray breast tumor samples with publicly available miRNA-seq breast tumor data. Additionally, we obtained possible reasons why miRNA can appear as a false positive in microarray study using paired miRNA sequencing and array data. The use of DNA microarrays for estimating miRNA expression profile is limited by several factors. One of them consists of problems with comparing expression values of different miRNAs. In this work, we show that situation can be significantly improved if some additional information is taken into consideration in a comparison.

A network analysis to identify mediators of germline-driven differences in breast cancer prognosis

Identifying the underlying genetic drivers of the heritability of breast cancer prognosis remains elusive. We adapt a network-based approach to handle underpowered complex datasets to provide new insights into the potential function of germline variants in breast cancer prognosis. This network-based analysis studies ~7.3 million variants in 84,457 breast cancer patients in relation to breast cancer survival and confirms the results on 12,381 independent patients. Aggregating the prognostic effects of genetic variants across multiple genes, we identify four gene modules associated with survival in estrogen receptor (ER)-negative and one in ER-positive disease. The modules show biological enrichment for cancer-related processes such as G-alpha signaling, circadian clock, angiogenesis, and Rho-GTPases in apoptosis.

Fine-mapping of 150 breast cancer risk regions identifies 191 likely target genes

Genome-wide association studies have identified breast cancer risk variants in over 150 genomic regions, but the mechanisms underlying risk remain largely unknown. These regions were explored by combining association analysis with in silico genomic feature annotations. We defined 205 independent risk-associated signals with the set of credible causal variants in each one. In parallel, we used a Bayesian approach (PAINTOR) that combines genetic association, linkage disequilibrium and enriched genomic features to determine variants with high posterior probabilities of being causal. Potentially causal variants were significantly over-represented in active gene regulatory regions and transcription factor binding sites. We applied our INQUSIT pipeline for prioritizing genes as targets of those potentially causal variants, using gene expression (expression quantitative trait loci), chromatin interaction and functional annotations. Known cancer drivers, transcription factors and genes in the developmental, apoptosis, immune system and DNA integrity checkpoint gene ontology pathways were over-represented among the highest-confidence target genes.

DNA methylation profiling to explore colorectal tumor differences according to menopausal hormone therapy use in women

Aim: Use of menopausal hormone therapy (MHT) has been associated with a reduced risk for colorectal cancer, but mechanisms underlying this relationship are not well understood. In the colon, MHT appears to act through estrogen receptor β (ERβ) which may influence DNA methylation by binding to DNA. Using genome-wide methylation profiling data, we aimed to identify genes that may be differentially methylated according to MHT use. Materials & methods: DNA methylation was measured using Illumina HumanMethylation450k arrays in two independent tumor sample sets of colorectal cancer patients. Differential methylation was determined using R/limma. Results: In the discovery analysis, two CpG sites showed differential DNA methylation according to MHT use, both were not replicated. In stratified analyses, 342 CpG sites were associated with current MHT use only in ERβ-positive tumors. Conclusion: The suggestive findings of differential methylation according to current MHT use in ERβ-positive tumors warrant further investigation.

The FANCM:p.Arg658* truncating variant is associated with risk of triple-negative breast cancer

Breast cancer is a common disease partially caused by genetic risk factors. Germline pathogenic variants in DNA repair genes BRCA1, BRCA2, PALB2, ATM, and CHEK2 are associated with breast cancer risk. FANCM, which encodes for a DNA translocase, has been proposed as a breast cancer predisposition gene, with greater effects for the ER-negative and triple-negative breast cancer (TNBC) subtypes. We tested the three recurrent protein-truncating variants FANCM:p.Arg658*, p.Gln1701*, and p.Arg1931* for association with breast cancer risk in 67,112 cases, 53,766 controls, and 26,662 carriers of pathogenic variants of BRCA1 or BRCA2. These three variants were also studied functionally by measuring survival and chromosome fragility in FANCM -/- patient-derived immortalized fibroblasts treated with diepoxybutane or olaparib. We observed that FANCM:p.Arg658* was associated with increased risk of ER-negative disease and TNBC (OR = 2.44, P = 0.034 and OR = 3.79; P = 0.009, respectively). In a country-restricted analysis, we confirmed the associations detected for FANCM:p.Arg658* and found that also FANCM:p.Arg1931* was associated with ER-negative breast cancer risk (OR = 1.96; P = 0.006). The functional results indicated that all three variants were deleterious affecting cell survival and chromosome stability with FANCM:p.Arg658* causing more severe phenotypes. In conclusion, we confirmed that the two rare FANCM deleterious variants p.Arg658* and p.Arg1931* are risk factors for ER-negative and TNBC subtypes. Overall our data suggest that the effect of truncating variants on breast cancer risk may depend on their position in the gene. Cell sensitivity to olaparib exposure, identifies a possible therapeutic option to treat FANCM-associated tumors.

Publisher Correction: Shared heritability and functional enrichment across six solid cancers

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Two truncating variants in FANCC and breast cancer risk

Fanconi anemia (FA) is a genetically heterogeneous disorder with 22 disease-causing genes reported to date. In some FA genes, monoallelic mutations have been found to be associated with breast cancer risk, while the risk associations of others remain unknown. The gene for FA type C, FANCC, has been proposed as a breast cancer susceptibility gene based on epidemiological and sequencing studies. We used the Oncoarray project to genotype two truncating FANCC variants (p.R185X and p.R548X) in 64,760 breast cancer cases and 49,793 controls of European descent. FANCC mutations were observed in 25 cases (14 with p.R185X, 11 with p.R548X) and 26 controls (18 with p.R185X, 8 with p.R548X). There was no evidence of an association with the risk of breast cancer, neither overall (odds ratio 0.77, 95%CI 0.44-1.33, p = 0.4) nor by histology, hormone receptor status, age or family history. We conclude that the breast cancer risk association of these two FANCC variants, if any, is much smaller than for BRCA1, BRCA2 or PALB2 mutations. If this applies to all truncating variants in FANCC it would suggest there are differences between FA genes in their roles on breast cancer risk and demonstrates the merit of large consortia for clarifying risk associations of rare variants.

A prognostic CpG score derived from epigenome-wide profiling of tumor tissue was independently associated with colorectal cancer survival

Background

Results of previous studies on the association of the CpG island methylator phenotype (CIMP) with colorectal cancer (CRC) prognosis were inconsistent and mostly based on different CIMP definitions. The current study aimed to comprehensively investigate the associations between DNA methylation on genes previously used to define CIMP status with CRC survival.

Results

Patients with CRC followed up for a median of 5.2 years were divided into a study cohort (n = 568) and a validation cohort (n = 308). DNA methylation was measured in tumor tissue using the Illumina Infinium HumanMethylation450 BeadChip and restricted to 43 genes used to define CIMP status in previous studies. Cox proportional hazard regression models were used to estimate adjusted hazard ratios (HR) and 95% confidence intervals (CI) of survival after CRC, including adjustment for tumor stage, microsatellite instability, and BRAF mutation status. In the study cohort, ten CpG sites were identified to be associated with CRC survival. Seven of these ten CpG sites were also associated with CRC survival in the validation cohort and were used to construct a prognostic score. CRC patients with a prognostic score of the lowest methylation level showed poorer disease-specific survival compared with patients with the highest methylation level in both the study cohort and the validation cohort (HR = 3.11 and 95% CI = 1.97–4.91, and HR = 3.06 and 95% CI = 1.71–5.45, respectively).

Conclusions

A CpG panel consisting of seven CpG sites was found to be strongly associated with CRC survival, independent from important clinical factors and mutations associated with CIMP. Further studies are required to confirm these findings.

Electronic supplementary material

The online version of this article (10.1186/s13148-019-0703-4) contains supplementary material, which is available to authorized users.

Genome-wide association and transcriptome studies identify target genes and risk loci for breast cancer

Genome-wide association studies (GWAS) have identified more than 170 breast cancer susceptibility loci. Here we hypothesize that some risk-associated variants might act in non-breast tissues, specifically adipose tissue and immune cells from blood and spleen. Using expression quantitative trait loci (eQTL) reported in these tissues, we identify 26 previously unreported, likely target genes of overall breast cancer risk variants, and 17 for estrogen receptor (ER)-negative breast cancer, several with a known immune function. We determine the directional effect of gene expression on disease risk measured based on single and multiple eQTL. In addition, using a gene-based test of association that considers eQTL from multiple tissues, we identify seven (and four) regions with variants associated with overall (and ER-negative) breast cancer risk, which were not reported in previous GWAS. Further investigation of the function of the implicated genes in breast and immune cells may provide insights into the etiology of breast cancer.

Genome-wide DNA methylation differences according to oestrogen receptor beta status in colorectal cancer

Involvement of sex hormones in colorectal cancer (CRC) development has been linked to oestrogen receptor β (ERβ). Expression of ERβ is found reduced in tumour tissue and inversely related to mortality. However, mechanisms are not well understood. Our study aimed to detect differentially methylated genes associated with ERβ expression, which could point to mechanisms by which ERβ could influence risk and prognosis of CRC. Epigenome-wide DNA methylation profiling was performed using Illumina HumanMethylation450k BeadChip arrays in two independent tumour sample sets of CRC patients recruited in 2003-2010 by the German DACHS study (discovery cohort n = 917, replication cohort n = 907). ERβ expression was measured using immunohistochemistry and scored as negative, moderate and high. Differentially methylated CpG sites and genomic regions were determined using limma in the R-package RnBeads. For the comparison of tumours with moderate/high ERβ versus negative expression, differentially methylated CpG sites were identified but not confirmed by replication. Comparing tumours of high with tumours of negative ERβ expression revealed 2,904 differentially methylated CpG sites of which 403 were replicated (FDR adjusted p-value<0.05). Replicated CpGs were annotated to genes such as CD36, HK1 or LRP5. A survival analysis indicates that 30 of the replicated CpGs are also associated with overall survival (FDR-adjusted p-value<0.05). The regional analysis identified 60 differentially methylated promotor regions. The epigenome-wide analysis identified both novel genes as well as genes already implicated in CRC. Follow-up mechanistic studies to better understand the regulatory role of ERβ could inform potential targets for improving treatment or prevention of CRC.

Genome-wide association study of germline variants and breast cancer-specific mortality

BACKGROUND

We examined the associations between germline variants and breast cancer mortality using a large meta-analysis of women of European ancestry.

METHODS

Meta-analyses included summary estimates based on Cox models of twelve datasets using ~10.4 million variants for 96,661 women with breast cancer and 7697 events (breast cancer-specific deaths). Oestrogen receptor (ER)-specific analyses were based on 64,171 ER-positive (4116) and 16,172 ER-negative (2125) patients. We evaluated the probability of a signal to be a true positive using the Bayesian false discovery probability (BFDP).

RESULTS

We did not find any variant associated with breast cancer-specific mortality at P < 5 × 10-8. For ER-positive disease, the most significantly associated variant was chr7:rs4717568 (BFDP = 7%, P = 1.28 × 10-7, hazard ratio [HR] = 0.88, 95% confidence interval [CI] = 0.84-0.92); the closest gene is AUTS2. For ER-negative disease, the most significant variant was chr7:rs67918676 (BFDP = 11%, P = 1.38 × 10-7, HR = 1.27, 95% CI = 1.16-1.39); located within a long intergenic non-coding RNA gene (AC004009.3), close to the HOXA gene cluster.

CONCLUSIONS

We uncovered germline variants on chromosome 7 at BFDP < 15% close to genes for which there is biological evidence related to breast cancer outcome. However, the paucity of variants associated with mortality at genome-wide significance underpins the challenge in providing genetic-based individualised prognostic information for breast cancer patients.

Shared heritability and functional enrichment across six solid cancers

Quantifying the genetic correlation between cancers can provide important insights into the mechanisms driving cancer etiology. Using genome-wide association study summary statistics across six cancer types based on a total of 296,215 cases and 301,319 controls of European ancestry, here we estimate the pair-wise genetic correlations between breast, colorectal, head/neck, lung, ovary and prostate cancer, and between cancers and 38 other diseases. We observed statistically significant genetic correlations between lung and head/neck cancer (rg = 0.57, p = 4.6 × 10-8), breast and ovarian cancer (rg = 0.24, p = 7 × 10-5), breast and lung cancer (rg = 0.18, p =1.5 × 10-6) and breast and colorectal cancer (rg = 0.15, p = 1.1 × 10-4). We also found that multiple cancers are genetically correlated with non-cancer traits including smoking, psychiatric diseases and metabolic characteristics. Functional enrichment analysis revealed a significant excess contribution of conserved and regulatory regions to cancer heritability. Our comprehensive analysis of cross-cancer heritability suggests that solid tumors arising across tissues share in part a common germline genetic basis.

The associations of DNA methylation alterations in oxidative stress-related genes with cancer incidence and mortality outcomes: a population-based cohort study

Background

Reactive oxygen species may be involved in epigenetic gene activation or silencing. We aimed to identify CpG sites, at which DNA methylation is related to urinary 8-isoprostane levels (biomarker of lipid peroxidation) and cancer or mortality outcomes. This investigation was based on a German, population-based cohort with linkage to cancer and mortality registry data (2000–2016).

Results

Blood DNA methylation in promoter regions of 519 genes, known to be involved in pathways from oxidative stress (OS) to cancer, was obtained at the cohort's baseline examination. Inverse associations of DNA methylation at cg25365794 (ALOXE3) and cg08862778 (MTOR) with 8-isoprostane levels were observed in a derivation set (n = 1000) and validated in two independent subsets of the cohort (n = 548 and n = 741). Multivariate regression models were used to evaluate the associations of DNA methylation at the two CpG sites with lung, colorectal, prostate, breast, and overall cancer incidence as well as CVD, cancer, and all-cause mortality. DNA methylation at cg25365794 (ALOXE3) was inversely associated with lung and prostate cancer incidence. DNA methylation at cg08862778 (MTOR) was associated with a 43% lower breast cancer incidence in the top vs. bottom tertile.

Conclusion

The finding for ALOXE3 may not be causal. As ALOXE3 is mainly expressed in skin tissue, the observed association might reflect the fact that both DNA methylation at the ALOXE3 gene and urinary 8-isoprostane concentrations depend on the level of OS in tissues. Contrarily, the finding for the MTOR gene and breast cancer is biologically plausible because the MTOR protein plays an important role in PI3K/Akt signaling, which is a pathway related to cancer development and cell senescence.

Electronic supplementary material

The online version of this article (10.1186/s13148-018-0604-y) contains supplementary material, which is available to authorized users.

A plasma metabolite panel as biomarkers for early primary breast cancer detection

In recent years, metabolites have attracted substantial attention as promising novel biomarkers of various diseases. However, breast cancer plasma metabolite studies are still in their infancy. Here, we investigated the potential of metabolites to serve as minimally invasive, early detection markers of primary breast cancer. We profiled metabolites extracted from the plasma of primary breast cancer patients and healthy controls using tandem mass spectrometry (UHPLC-MS/MS and FIA-MS/MS). Two metabolites were found to be upregulated, while 16 metabolites were downregulated in primary breast cancer patients compared to healthy controls in both the training and validation cohorts. A panel of seven metabolites was selected by LASSO regression analysis. This panel could differentiate primary breast cancer patients from healthy controls, with an AUC of 0.87 (95% CI: 0.81 ~ 0.92) in the training cohort and an AUC of 0.80 (95% CI: 0.71 ~ 0.87) in the validation cohort. These significantly differentiated metabolites are mainly involved in the amino acid metabolism and breast cancer cell growth pathways. In conclusion, using a metabolomics approach, we identified metabolites that have potential value for development of a multimarker blood-based test to complement and improve early breast cancer detection. The panel identified herein might be part of a prescreening tool, especially for younger women or for closely observing women with certain risks, to facilitate decision making regarding which individuals should undergo further diagnostic tests. In the future, the combination of metabolites and other blood-based molecular marker sets, such as DNA methylation, microRNA, and cell-free DNA mutation markers, will be an attractive option.

Polygenic Risk Scores for Prediction of Breast Cancer and Breast Cancer Subtypes

Stratification of women according to their risk of breast cancer based on polygenic risk scores (PRSs) could improve screening and prevention strategies. Our aim was to develop PRSs, optimized for prediction of estrogen receptor (ER)-specific disease, from the largest available genome-wide association dataset and to empirically validate the PRSs in prospective studies. The development dataset comprised 94,075 case subjects and 75,017 control subjects of European ancestry from 69 studies, divided into training and validation sets. Samples were genotyped using genome-wide arrays, and single-nucleotide polymorphisms (SNPs) were selected by stepwise regression or lasso penalized regression. The best performing PRSs were validated in an independent test set comprising 11,428 case subjects and 18,323 control subjects from 10 prospective studies and 190,040 women from UK Biobank (3,215 incident breast cancers). For the best PRSs (313 SNPs), the odds ratio for overall disease per 1 standard deviation in ten prospective studies was 1.61 (95%CI: 1.57-1.65) with area under receiver-operator curve (AUC) = 0.630 (95%CI: 0.628-0.651). The lifetime risk of overall breast cancer in the top centile of the PRSs was 32.6%. Compared with women in the middle quintile, those in the highest 1% of risk had 4.37- and 2.78-fold risks, and those in the lowest 1% of risk had 0.16- and 0.27-fold risks, of developing ER-positive and ER-negative disease, respectively. Goodness-of-fit tests indicated that this PRS was well calibrated and predicts disease risk accurately in the tails of the distribution. This PRS is a powerful and reliable predictor of breast cancer risk that may improve breast cancer prevention programs.

Genome-wide DNA methylation analysis reveals a prognostic classifier for non-metastatic colorectal cancer (ProMCol classifier)

OBJECTIVE

Pathological staging used for the prediction of patient survival in colorectal cancer (CRC) provides only limited information.

DESIGN

Here, a genome-wide study of DNA methylation was conducted for two cohorts of patients with non-metastatic CRC (screening cohort (n=572) and validation cohort (n=274)). A variable screening for prognostic CpG sites was performed in the screening cohort using marginal testing based on a Cox model and subsequent adjustment of the p-values via independent hypothesis weighting using the methylation difference between 34 pairs of tumour and normal mucosa tissue as auxiliary covariate. From the 1000 CpG sites with the smallest adjusted p-value, 20 CpG sites with the smallest Brier score for overall survival (OS) were selected. Applying principal component analysis, we derived a prognostic methylation-based classifier for patients with non-metastatic CRC (ProMCol classifier).

RESULTS

This classifier was associated with OS in the screening (HR 0.51, 95% CI 0.41 to 0.63, p=6.2E-10) and the validation cohort (HR 0.61, 95% CI 0.45 to 0.82, p=0.001). The independent validation of the ProMCol classifier revealed a reduction of the prediction error for 3-year OS from 0.127, calculated only with standard clinical variables, to 0.120 combining the clinical variables with the classifier and for 4-year OS from 0.153 to 0.140. All results were confirmed for disease-specific survival.

CONCLUSION

The ProMCol classifier could improve the prognostic accuracy for patients with non-metastatic CRC.

Identification of nine new susceptibility loci for endometrial cancer

Endometrial cancer is the most commonly diagnosed cancer of the female reproductive tract in developed countries. Through genome-wide association studies (GWAS), we have previously identified eight risk loci for endometrial cancer. Here, we present an expanded meta-analysis of 12,906 endometrial cancer cases and 108,979 controls (including new genotype data for 5624 cases) and identify nine novel genome-wide significant loci, including a locus on 12q24.12 previously identified by meta-GWAS of endometrial and colorectal cancer. At five loci, expression quantitative trait locus (eQTL) analyses identify candidate causal genes; risk alleles at two of these loci associate with decreased expression of genes, which encode negative regulators of oncogenic signal transduction proteins (SH2B3 (12q24.12) and NF1 (17q11.2)). In summary, this study has doubled the number of known endometrial cancer risk loci and revealed candidate causal genes for future study.

Pathway analysis of genetic variants in folate-mediated one-carbon metabolism-related genes and survival in a prospectively followed cohort of colorectal cancer patients

Folate-mediated one-carbon metabolism (FOCM) is a key pathway essential for nucleotide synthesis, DNA methylation, and repair. This pathway is a critical target for 5-fluorouracil (5-FU), which is predominantly used for colorectal cancer (CRC) treatment. A comprehensive assessment of polymorphisms in FOCM-related genes and their association with prognosis has not yet been performed. Within 1,739 CRC cases aged ≥30 years diagnosed from 2003 to 2007 (DACHS study), we investigated 397 single nucleotide polymorphisms (SNPs) and 50 candidates in 48 FOCM-related genes for associations with overall- (OS) and disease-free survival (DFS) using multiple Cox regression (adjusted for age, sex, stage, grade, BMI, and alcohol). We investigated effect modification by 5-FU-based chemotherapy and assessed pathway-specific effects. Correction for multiple testing was performed using false discovery rates (FDR). After a median follow-up time of 5.0 years, 585 patients were deceased. For one candidate SNP in MTHFR and two in TYMS, we observed significant inverse associations with OS (MTHFR: rs1801133, C677T: HRhet  = 0.81, 95% CI: 0.67-0.97; TYMS: rs1001761: HRhet  = 0.82, 95% CI: 0.68-0.99 and rs2847149: HRhet  = 0.82, 95% CI: 0.68-0.99). After FDR correction, one polymorphism in paraoxonase 1 (PON1; rs3917538) was significantly associated with OS (HRhet  = 1.28, 95% CI: 1.07-1.53; HRhzv  = 2.02, 95% CI:1.46-2.80; HRlogAdd  = 1.31, pFDR  = 0.01). Adjusted pathway analyses showed significant associations for pyrimidine biosynthesis (P = 0.04) and fluorouracil drug metabolism (P < 0.01) with significant gene-chemotherapy interactions, including PON1 rs3917538. This study supports the concept that FOCM-related genes could be associated with CRC survival and may modify effects of 5-FU-based chemotherapy in genes in pyrimidine and fluorouracil metabolism, which are relevant targets for therapeutic response and prognosis in CRC. These results require confirmation in additional clinical studies.

A transcriptome-wide association study of 229,000 women identifies new candidate susceptibility genes for breast cancer

The breast cancer risk variants identified in genome-wide association studies explain only a small fraction of the familial relative risk, and the genes responsible for these associations remain largely unknown. To identify novel risk loci and likely causal genes, we performed a transcriptome-wide association study evaluating associations of genetically predicted gene expression with breast cancer risk in 122,977 cases and 105,974 controls of European ancestry. We used data from the Genotype-Tissue Expression Project to establish genetic models to predict gene expression in breast tissue and evaluated model performance using data from The Cancer Genome Atlas. Of the 8,597 genes evaluated, significant associations were identified for 48 at a Bonferroni-corrected threshold of P < 5.82 × 10-6, including 14 genes at loci not yet reported for breast cancer. We silenced 13 genes and showed an effect for 11 on cell proliferation and/or colony-forming efficiency. Our study provides new insights into breast cancer genetics and biology.

Genetic overlap between endometriosis and endometrial cancer: evidence from cross-disease genetic correlation and GWAS meta-analyses

Epidemiological, biological, and molecular data suggest links between endometriosis and endometrial cancer, with recent epidemiological studies providing evidence for an association between a previous diagnosis of endometriosis and risk of endometrial cancer. We used genetic data as an alternative approach to investigate shared biological etiology of these two diseases. Genetic correlation analysis of summary level statistics from genomewide association studies (GWAS) using LD Score regression revealed moderate but significant genetic correlation (rg  = 0.23, P = 9.3 × 10-3 ), and SNP effect concordance analysis provided evidence for significant SNP pleiotropy (P = 6.0 × 10-3 ) and concordance in effect direction (P = 2.0 × 10-3 ) between the two diseases. Cross-disease GWAS meta-analysis highlighted 13 distinct loci associated at P ≤ 10-5 with both endometriosis and endometrial cancer, with one locus (SNP rs2475335) located within PTPRD associated at a genomewide significant level (P = 4.9 × 10-8 , OR = 1.11, 95% CI = 1.07-1.15). PTPRD acts in the STAT3 pathway, which has been implicated in both endometriosis and endometrial cancer. This study demonstrates the value of cross-disease genetic analysis to support epidemiological observations and to identify biological pathways of relevance to multiple diseases.

The BRCA2 c.68-7T > A variant is not pathogenic: A model for clinical calibration of spliceogenicity

Although the spliceogenic nature of the BRCA2 c.68-7T > A variant has been demonstrated, its association with cancer risk remains controversial. In this study, we accurately quantified by real-time PCR and digital PCR (dPCR), the BRCA2 isoforms retaining or missing exon 3. In addition, the combined odds ratio for causality of the variant was estimated using genetic and clinical data, and its associated cancer risk was estimated by case-control analysis in 83,636 individuals. Co-occurrence in trans with pathogenic BRCA2 variants was assessed in 5,382 families. Exon 3 exclusion rate was 4.5-fold higher in variant carriers (13%) than controls (3%), indicating an exclusion rate for the c.68-7T > A allele of approximately 20%. The posterior probability of pathogenicity was 7.44 × 10-115 . There was neither evidence for increased risk of breast cancer (OR 1.03; 95% CI 0.86-1.24) nor for a deleterious effect of the variant when co-occurring with pathogenic variants. Our data provide for the first time robust evidence of the nonpathogenicity of the BRCA2 c.68-7T > A. Genetic and quantitative transcript analyses together inform the threshold for the ratio between functional and altered BRCA2 isoforms compatible with normal cell function. These findings might be exploited to assess the relevance for cancer risk of other BRCA2 spliceogenic variants.

Publisher Correction: Evidence that breast cancer risk at the 2q35 locus is mediated through IGFBP5 regulation

This corrects the article DOI: 10.1038/ncomms5999.

Evaluation of Promoter Methylation of RASSF1A and ATM in Peripheral Blood of Breast Cancer Patients and Healthy Control Individuals

Breast cancer (BC) is the most common cancer among women and has high mortality rates. Early detection is supposed to be critical for the patient’s prognosis. In recent years, several studies have investigated global DNA methylation profiles and gene-specific DNA methylation in blood-based DNA to develop putative screening markers for cancer. However, most of the studies have not yet been validated. In our study, we analyzed the promoter methylation of RASSF1A and ATM in peripheral blood DNA of 229 sporadic patients and 151 healthy controls by the MassARRAY EpiTYPER assay. There were no significant differences in DNA methylation levels of RASSF1A and ATM between the sporadic BC cases and the healthy controls. Furthermore, we performed the Infinium HumanMethylation450 BeadChip (450K) array analysis using 48 sporadic BC cases and 48 healthy controls (cases and controls are the same from those of the MassARRAY EpiTYPER assay) and made a comparison with the published data. No significant differences were presented in DNA methylation levels of RASSF1A and ATM between the sporadic BC cases and the healthy controls. So far, the evidence for powerful blood-based methylation markers is still limited and the identified markers need to be further validated.

Circulating free DNA integrity and concentration as independent prognostic markers in metastatic breast cancer

PURPOSE

Non-invasive blood-based molecular markers have been investigated for cancer diagnosis and prognosis. Circulating free or cell-free DNA (cfDNA) variables have been shown to be putative markers in breast cancer prognosis.

METHODS

Here, we investigated the potential prognostic ability of cfDNA concentration and cfDNA integrity (cfDI) in a study cohort of 268 patients by quantitative PCR. We compared cfDNA concentration and cfDI at baseline and after one cycle of therapy in metastatic breast cancer (MBC) patients.

RESULTS

A significantly increased cfDI (P = 1.21E-7 for ALU and P = 1.87E-3 for LINE1) and decreased cfDNA concentration (P = 1.17E-3 for ALU and P = 1.60E-2 for LINE1) in both repetitive DNA elements after one cycle of therapy was observed. A multiple Cox regression model indicated that cfDI and cfDNA concentration can serve as independent prognostic markers in patients at baseline with HR (95% CI) of 0.70 (0.48-1.01) for ALU cfDI, 0.63 (0.44-0.92) for LINE1 cfDI, 2.44 (1.68-3.53) for ALU cfDNA concentration, and 2.12 (1.47-3.06) for LINE1 cfDNA concentration and after one cycle of therapy with HR (95% CI) of 0.59 (0.42-0.84) for ALU cfDI, 0.51 (0.36-0.74) for LINE1 cfDI, 1.59 (1.31-1.92) for ALU cfDNA concentration, and 1.30 (1.17-1.45) for LINE1 cfDNA concentration, respectively. By comparing integrated prediction error of different models, cfDNA variables were shown to improve the prognostic power of the CTC status.

CONCLUSIONS

We hereby show that cfDNA variables, especially in combination with other markers, can serve as attractive prognostic markers for MBC patients at baseline and during the systematic therapy.

Circulating microRNA's as a diagnostic tool for hepatocellular carcinoma in a hyper endemic HIV setting, KwaZulu-Natal, South Africa: a case control study protocol focusing on viral etiology

Background

A wide range of studies has investigated the diagnostic proficiency of extracellular microRNAs (miRNAs) in hepatocellular cancer (HCC). HCC is expected to increase in Sub-Saharan Africa (SSA), due to endemic levels of viral infection (HBV/HIV), ageing and changing lifestyles. This unique aetiological background provides an opportunity for investigating potentially novel circulating miRNAs as biomarkers for HCC in a prospective study in South Africa.

Methods

This study will recruit HCC patients from two South African cancer hospitals, situated in Durban and Pietermaritzburg in the province of KwaZulu-Natal. These cases will include both HBV mono-infected and HBV/HIV co-infected HCC cases. The control group will consist of two (2) age and sex-matched healthy population controls per HCC case randomly selected from a Durban based laboratory. The controls will exclude patients if they have any evidence of chronic liver disease. A standardised reporting approach will be adopted to detect, quantify and normalize the level of circulating miRNAs in the blood sera of HCC cases and their controls. Reverse transcription quantitative polymerase chain reaction (RT-qPCR) will be employed to quantity extracellular miRNAs. Differences in concentration of relevant miRNA by case/control status will be assessed using the Wilcoxon rank-sum (Mann-Whitney U) test. Adjustment for multiple testing (Bonferroni correction), receiver operating curves (ROC) and optimal breakpoint analyses will be employed to identify potential thresholds for the differentiation of miRNA levels of HCC cases and their controls.

Discussion

Although there is a growing base of literature regarding the role of circulating miRNAs as biomarkers, this promising field remains a ‘work in progress’. The aetiology of HBV infection in HCC is well understood, as well as it’s role in miRNA deregulation, however, the mediating role of HIV infection is unknown. HCC incidence in SSA, including South Africa, is expected to increase significantly in the next decade. A combination of factors, therefore, offers a unique opportunity to identify candidate circulating miRNAs as potential biomarkers for HBV/HIV infected HCC.

Electronic supplementary material

The online version of this article (10.1186/s12885-017-3915-z) contains supplementary material, which is available to authorized users.

Body mass index and breast cancer survival: a Mendelian randomization analysis

Background

There is increasing evidence that elevated body mass index (BMI) is associated with reduced survival for women with breast cancer. However, the underlying reasons remain unclear. We conducted a Mendelian randomization analysis to investigate a possible causal role of BMI in survival from breast cancer.

Methods

We used individual-level data from six large breast cancer case-cohorts including a total of 36 210 individuals (2475 events) of European ancestry. We created a BMI genetic risk score (GRS) based on genotypes at 94 known BMI-associated genetic variants. Association between the BMI genetic score and breast cancer survival was analysed by Cox regression for each study separately. Study-specific hazard ratios were pooled using fixed-effect meta-analysis.

Results

BMI genetic score was found to be associated with reduced breast cancer-specific survival for estrogen receptor (ER)-positive cases [hazard ratio (HR) = 1.11, per one-unit increment of GRS, 95% confidence interval (CI) 1.01-1.22, P = 0.03). We observed no association for ER-negative cases (HR = 1.00, per one-unit increment of GRS, 95% CI 0.89-1.13, P = 0.95).

Conclusions

Our findings suggest a causal effect of increased BMI on reduced breast cancer survival for ER-positive breast cancer. There is no evidence of a causal effect of higher BMI on survival for ER-negative breast cancer cases.

Identification of ten variants associated with risk of estrogen-receptor-negative breast cancer

Most common breast cancer susceptibility variants have been identified through genome-wide association studies (GWAS) of predominantly estrogen receptor (ER)-positive disease. We conducted a GWAS using 21,468 ER-negative cases and 100,594 controls combined with 18,908 BRCA1 mutation carriers (9,414 with breast cancer), all of European origin. We identified independent associations at P < 5 × 10-8 with ten variants at nine new loci. At P < 0.05, we replicated associations with 10 of 11 variants previously reported in ER-negative disease or BRCA1 mutation carrier GWAS and observed consistent associations with ER-negative disease for 105 susceptibility variants identified by other studies. These 125 variants explain approximately 16% of the familial risk of this breast cancer subtype. There was high genetic correlation (0.72) between risk of ER-negative breast cancer and breast cancer risk for BRCA1 mutation carriers. These findings may lead to improved risk prediction and inform further fine-mapping and functional work to better understand the biological basis of ER-negative breast cancer.

PHIP - a novel candidate breast cancer susceptibility locus on 6q14.1

Most non-BRCA1/2 breast cancer families have no identified genetic cause. We used linkage and haplotype analyses in familial and sporadic breast cancer cases to identify a susceptibility locus on chromosome 6q. Two independent genome-wide linkage analysis studies suggested a 3 Mb locus on chromosome 6q and two unrelated Swedish families with a LOD >2 together seemed to share a haplotype in 6q14.1. We hypothesized that this region harbored a rare high-risk founder allele contributing to breast cancer in these two families. Sequencing of DNA and RNA from the two families did not detect any pathogenic mutations. Finally, 29 SNPs in the region were analyzed in 44,214 cases and 43,532 controls from BCAC, and the original haplotypes in the two families were suggested as low-risk alleles for European and Swedish women specifically. There was also some support for one additional independent moderate-risk allele in Swedish familial samples. The results were consistent with our previous findings in familial breast cancer and supported a breast cancer susceptibility locus at 6q14.1 around the PHIP gene.

Association analysis identifies 65 new breast cancer risk loci

Breast cancer risk is influenced by rare coding variants in susceptibility genes, such as BRCA1, and many common, mostly non-coding variants. However, much of the genetic contribution to breast cancer risk remains unknown. Here we report the results of a genome-wide association study of breast cancer in 122,977 cases and 105,974 controls of European ancestry and 14,068 cases and 13,104 controls of East Asian ancestry. We identified 65 new loci that are associated with overall breast cancer risk at P < 5 × 10-8. The majority of credible risk single-nucleotide polymorphisms in these loci fall in distal regulatory elements, and by integrating in silico data to predict target genes in breast cells at each locus, we demonstrate a strong overlap between candidate target genes and somatic driver genes in breast tumours. We also find that heritability of breast cancer due to all single-nucleotide polymorphisms in regulatory features was 2-5-fold enriched relative to the genome-wide average, with strong enrichment for particular transcription factor binding sites. These results provide further insight into genetic susceptibility to breast cancer and will improve the use of genetic risk scores for individualized screening and prevention.