Posaconazole-ibrutinib interaction cannot be avoided by staggered dosing: How to optimize ibrutinib dose during posaconazole treatment

AIMS

Ibrutinib is used in the treatment of certain B-cell malignancies. Due to its CYP3A4-mediated metabolism and highly variable pharmacokinetics, it is prone to potentially harmful drug-drug interactions.

METHODS

In a randomized, placebo-controlled, three-phase crossover study, we examined the effect of the CYP3A4-inhibiting antifungal posaconazole on ibrutinib pharmacokinetics. Eleven healthy participants ingested repeated doses of 300 mg of posaconazole either in the morning or in the evening, or placebo. A single dose of ibrutinib (30, 70 or 140 mg, respectively) was administered at 9 AM, 1 or 12 h after the preceding posaconazole/placebo dose.

RESULTS

On average, morning posaconazole increased the dose-adjusted geometric mean area under the plasma concentration-time curve from zero to infinity (AUC0-∞ ) and peak plasma concentration (Cmax ) of ibrutinib 9.5-fold (90% confidence interval [CI] 6.3-14.3, P < 0.001) and 8.5-fold (90% CI 5.7-12.8, P < 0.001), respectively, while evening posaconazole increased those 10.3-fold (90% CI 6.7-16.0, P < 0.001) and 8.2-fold (90% CI 5.2-13.2, P < 0.001), respectively. Posaconazole had no significant effect on the half-life of ibrutinib, but substantially reduced the metabolite PCI-45227 to ibrutinib AUC0-∞ ratio. There were no significant differences in ibrutinib pharmacokinetics between morning and evening posaconazole phases.

CONCLUSIONS

Posaconazole increases ibrutinib exposure substantially, by about 10-fold. This interaction cannot be avoided by dosing the drugs 12 h apart. In general, a 70-mg daily dose of ibrutinib should not be exceeded during posaconazole treatment to avoid potentially toxic systemic ibrutinib concentrations.

Solanidine is a sensitive and specific dietary biomarker for CYP2D6 activity

BACKGROUND

Individual assessment of CYP enzyme activities can be challenging. Recently, the potato alkaloid solanidine was suggested as a biomarker for CYP2D6 activity. Here, we aimed to characterize the sensitivity and specificity of solanidine as a CYP2D6 biomarker among Finnish volunteers with known CYP2D6 genotypes.

RESULTS

Using non-targeted metabolomics analysis, we identified 9152 metabolite features in the fasting plasma samples of 356 healthy volunteers. Machine learning models suggested strong association between CYP2D6 genotype-based phenotype classes with a metabolite feature identified as solanidine. Plasma solanidine concentration was 1887% higher in genetically poor CYP2D6 metabolizers (gPM) (n = 9; 95% confidence interval 755%, 4515%; P = 1.88 × 10-11), 74% higher in intermediate CYP2D6 metabolizers (gIM) (n = 89; 27%, 138%; P = 6.40 × 10-4), and 35% lower in ultrarapid CYP2D6 metabolizers (gUM) (n = 20; 64%, - 17%; P = 0.151) than in genetically normal CYP2D6 metabolizers (gNM; n = 196). The solanidine metabolites m/z 444 and 430 to solanidine concentration ratios showed even stronger associations with CYP2D6 phenotypes. Furthermore, the areas under the receiver operating characteristic and precision-recall curves for these metabolic ratios showed equal or better performances for identifying the gPM, gIM, and gUM phenotype groups than the other metabolites, their ratios to solanidine, or solanidine alone. In vitro studies with human recombinant CYP enzymes showed that solanidine was metabolized mainly by CYP2D6, with a minor contribution from CYP3A4/5. In human liver microsomes, the CYP2D6 inhibitor paroxetine nearly completely (95%) inhibited the metabolism of solanidine. In a genome-wide association study, several variants near the CYP2D6 gene associated with plasma solanidine metabolite ratios.

CONCLUSIONS

These results are in line with earlier studies and further indicate that solanidine and its metabolites are sensitive and specific biomarkers for measuring CYP2D6 activity. Since potato consumption is common worldwide, this biomarker could be useful for evaluating CYP2D6-mediated drug-drug interactions and to improve prediction of CYP2D6 activity in addition to genotyping.

Ticagrelor Increases Exposure to the Breast Cancer Resistance Protein Substrate Rosuvastatin

Ticagrelor and rosuvastatin are often used concomitantly after atherothrombotic events. Several cases of rhabdomyolysis during concomitant ticagrelor and rosuvastatin have been reported, suggesting a drug-drug interaction. We showed recently that ticagrelor inhibits breast cancer resistance protein (BCRP) and organic anion transporting polypeptide (OATP) 1B1, 1B3, and 2B1-mediated rosuvastatin transport in vitro. The aim of this study was to investigate the effects of ticagrelor on rosuvastatin pharmacokinetics in humans. In a randomized, crossover study, 9 healthy volunteers ingested a single dose of 90 mg ticagrelor or placebo, followed by a single 10 mg dose of rosuvastatin 1 hour later. Ticagrelor 90 mg or placebo were additionally administered 12, 24, and 36 hours after their first dose. Ticagrelor increased rosuvastatin area under the plasma concentration-time curve (AUC) and peak plasma concentration 2.6-fold (90% confidence intervals: 1.8-3.8 and 1.7-4.0, P = 0.001 and P = 0.003), and prolonged its half-life from 3.1 to 6.6 hours (P = 0.009). Ticagrelor also decreased the renal clearance of rosuvastatin by 11% (3%-19%, P = 0.032). The N-desmethylrosuvastatin:rosuvastatin AUC0-10h ratio remained unaffected by ticagrelor. Ticagrelor had no effect on the plasma concentrations of the endogenous OATP1B substrates glycodeoxycholate 3-O-glucuronide, glycochenodeoxycholate 3-O-glucuronide, glycodeoxycholate 3-O-sulfate, and glycochenodeoxycholate 3-O-sulfate, or the sodium-taurocholate cotransporting polypeptide substrate taurocholic acid. These data indicate that ticagrelor increases rosuvastatin concentrations more than twofold in humans, probably mainly by inhibiting intestinal BCRP. Because the risk for rosuvastatin-induced myotoxicity increases along with rosuvastatin plasma concentrations, using ticagrelor concomitantly with high doses of rosuvastatin should be avoided.

NTHL1 is a recessive cancer susceptibility gene

In search of novel breast cancer (BC) risk variants, we performed a whole-exome sequencing and variant analysis of 69 Finnish BC patients as well as analysed loss-of-function variants identified in DNA repair genes in the Finns from the Genome Aggregation Database. Additionally, we carried out a validation study of SERPINA3 c.918-1G>C, recently suggested for BC predisposition. We estimated the frequencies of 41 rare candidate variants in 38 genes by genotyping them in 2482-4101 BC patients and in 1273-3985 controls. We further evaluated all coding variants in the candidate genes in a dataset of 18,786 BC patients and 182,927 controls from FinnGen. None of the variants associated significantly with cancer risk in the primary BC series; however, in the FinnGen data, NTHL1 c.244C>T p.(Gln82Ter) associated with BC with a high risk for homozygous (OR = 44.7 [95% CI 6.90-290], P = 6.7 × 10-5) and a low risk for heterozygous women (OR = 1.39 [1.18-1.64], P = 7.8 × 10-5). Furthermore, the results suggested a high risk of colorectal, urinary tract, and basal-cell skin cancer for homozygous individuals, supporting NTHL1 as a recessive multi-tumour susceptibility gene. No significant association with BC risk was detected for SERPINA3 or any other evaluated gene.

CYP2C19 loss-of-function alleles and use of omeprazole or esomeprazole increase the risk of cardiovascular outcomes in patients using clopidogrel

Our aim was to investigate in a real-life prospective patient cohort how CYP2C19 loss-of-function (LOF) variants and CYP2C19 inhibitor omeprazole or esomeprazole influence the incidence of cardiovascular events in patients using clopidogrel. Data based simultaneously on these factors are conflicting and sparse. A cohort of prospective patients (n = 1972) with acute coronary syndrome (n = 1302) or symptomatic chronic coronary disease (n = 656) was followed for 365 days after hospitalization with information on purchased prescription drugs, hospital discharge, death, and genotype for CYP2C19*2, CYP2C19*3, and CYP2C19*8 LOF variants. The primary study outcome measurement was cardiovascular death or recurring myocardial infarction or stroke. Altogether, 608 patients (30.8%) carried CYP2C19 LOF alleles. During the 365-day follow-up 252 patients (12.8%) had an ischemic vascular event. Cardiovascular events were significantly more frequent in carriers of CYP2C19 LOF alleles (14.8%, 95% confidence interval [CI], 11.7-17.8) than in non-carriers (10.8%, 95% CI, 9.0-12.6, p = 0.0159). Omeprazole or esomeprazole use was similar among LOF allele carriers (n = 131, 21.5%) and non-carriers (n = 250, 18.3%, p = 0.185). Cardiovascular events were significantly more common in a composite group consisting of all CYP2C19 LOF carriers regardless of proton pump inhibitor use status and non-carriers using omeprazole or esomeprazole than in non-carriers not using omeprazole or esomeprazole (14.8%, 95% CI, 12.2-17.3 vs. 9.9%, 95% CI, 8.0-11.9, p = 0.00173). We observed significantly more cardiovascular events in carriers of CYP2C19 LOF variants and in non-carriers using omeprazole or esomeprazole. For optimal patient care, both genetics and concomitant medication should be considered.

Real-world pharmacogenetics of statin intolerance: effects of SLCO1B1, ABCG2 , and CYP2C9 variants

OBJECTIVE

The association of SLCO1B1 c.521T>C with simvastatin-induced muscle toxicity is well characterized. However, different statins are subject to metabolism and transport also by other proteins exhibiting clinically meaningful genetic variation. Our aim was to investigate associations of SLCO1B1 c.521T>C with intolerance to atorvastatin, fluvastatin, pravastatin, rosuvastatin, or simvastatin, those of ABCG2 c.421C>A with intolerance to atorvastatin, fluvastatin, or rosuvastatin, and that of CYP2C9*2 and *3 alleles with intolerance to fluvastatin.

METHODS

We studied the associations of these variants with statin intolerance in 2042 patients initiating statin therapy by combining genetic data from samples from the Helsinki Biobank to clinical chemistry and statin purchase data.

RESULTS

We confirmed the association of SLCO1B1 c.521C/C genotype with simvastatin intolerance both by using phenotype of switching initial statin to another as a marker of statin intolerance [hazard ratio (HR) 1.88, 95% confidence interval (CI) 1.08-3.25, P  = 0.025] and statin switching along with creatine kinase measurement (HR 5.44, 95% CI 1.49-19.9, P  = 0.011). No significant association was observed with atorvastatin and rosuvastatin. The sample sizes for fluvastatin and pravastatin were relatively small, but SLCO1B1 c.521T>C carriers had an increased risk of pravastatin intolerance defined by statin switching when compared to homozygous reference T/T genotype (HR 2.11, 95% CI 1.01-4.39, P  = 0.047).

CONCLUSION

The current results can inform pharmacogenetic statin prescribing guidelines and show feasibility for the methodology to be used in larger future studies.

CYP2C19 loss-of-function alleles and the use of omeprazole or esomeprazole increase the risk of cardiovascular outcomes in patients using clopidogrel – a real-world study

Background

Patients carrying loss-of-function (LOF) variants of CYP2C19 or using the CYP2C19-inhibitors omeprazole or esomeprazole cannot fully bioactivate clopidogrel, and thus are lacking optimal thrombosis prevention. However, clopidogrel is a relevant alternative to ticagrelor due to its lower incidence of hemorrhagic events and lower cost.

Purpose

The purpose of this study was to investigate in a real-life prospective patient cohort how loss-of-function (LOF) CYP2C19 variants and omeprazole or esomeprazole influence the incidence of cardiovascular events (cardiovascular deaths, reinfarction and strokes) in patients using clopidogrel.

Methods

A cohort of prospective patients (N=1972) with acute coronary syndrome (ACS) (N=1302) or symptomatic chronic coronary disease (N=670) were followed for 365 days. Information of purchased prescription drugs, hospital discharge, and death registry was available for all patients. The patients were retrospectively genotyped for the CYP2C19*2, CYP2C19*3, and CYP2C19*8 LOF variants.

Results

A total of 608 (30.8%) patients carried CYP2C19 LOF alleles. During the 365 days follow-up, 252 patients had a recurring ischemic vascular event, either myocardial infarction (N=189), stroke (N=24) or cardiovascular death (N=39). Vascular events were significantly more frequent among carriers of CYP2C19 LOF alleles [14.8% (95% CI 11.7–17.8] than in non-carriers [10.8% (95% CI 9.0–12.6)], p=0.016. Omeprazole or esomeprazole use was similar among the LOF allele carriers and non-carriers, but inclusion of their use in the analysis strengthened the association. Vascular events were significantly more common in a composite group of carriers of CYP2C19 LOF alleles and non-carriers who were using omeprazole or esomeprazole than in non-carriers who were not using omeprazole or esomeprazole 14.8% (95% CI 12.2–17.3) vs 9.9% (95% CI 8.0–11.9) p=0.002. There was a tendency towards lower risk of hemorrhagic complications in carriers of LOF alleles than in non-carriers (p=0.07).

Conclusions

Patients with impaired clopidogrel bioactivation due to loss-of-function CYP2C19 variants or concomitant omeprazole or esomeprazole use are at increased risk of cardiovascular events. For relevant patient care both genetics and concomitant medication should be considered both in everyday patient care and further randomised studies.

Funding Acknowledgement

Type of funding sources: None.

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.

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.

Association of Genomic Domains in BRCA1 and BRCA2 with Prostate Cancer Risk and Aggressiveness

Pathogenic sequence variants (PSV) in BRCA1 or BRCA2 (BRCA1/2) are associated with increased risk and severity of prostate cancer. We evaluated whether PSVs in BRCA1/2 were associated with risk of overall prostate cancer or high grade (Gleason 8+) prostate cancer using an international sample of 65 BRCA1 and 171 BRCA2 male PSV carriers with prostate cancer, and 3,388 BRCA1 and 2,880 BRCA2 male PSV carriers without prostate cancer. PSVs in the 3' region of BRCA2 (c.7914+) were significantly associated with elevated risk of prostate cancer compared with reference bin c.1001-c.7913 [HR = 1.78; 95% confidence interval (CI), 1.25-2.52; P = 0.001], as well as elevated risk of Gleason 8+ prostate cancer (HR = 3.11; 95% CI, 1.63-5.95; P = 0.001). c.756-c.1000 was also associated with elevated prostate cancer risk (HR = 2.83; 95% CI, 1.71-4.68; P = 0.00004) and elevated risk of Gleason 8+ prostate cancer (HR = 4.95; 95% CI, 2.12-11.54; P = 0.0002). No genotype-phenotype associations were detected for PSVs in BRCA1. These results demonstrate that specific BRCA2 PSVs may be associated with elevated risk of developing aggressive prostate cancer. SIGNIFICANCE: Aggressive prostate cancer risk in BRCA2 mutation carriers may vary according to the specific BRCA2 mutation inherited by the at-risk individual.

The Spectrum of FANCM Protein Truncating Variants in European Breast Cancer Cases

Germline protein truncating variants (PTVs) in the FANCM gene have been associated with a 2–4-fold increased breast cancer risk in case-control studies conducted in different European populations. However, the distribution and the frequency of FANCM PTVs in Europe have never been investigated. In the present study, we collected the data of 114 European female breast cancer cases with FANCM PTVs ascertained in 20 centers from 13 European countries. We identified 27 different FANCM PTVs. The p.Gln1701* PTV is the most common PTV in Northern Europe with a maximum frequency in Finland and a lower relative frequency in Southern Europe. On the contrary, p.Arg1931* seems to be the most common PTV in Southern Europe. We also showed that p.Arg658*, the third most common PTV, is more frequent in Central Europe, and p.Gln498Thrfs*7 is probably a founder variant from Lithuania. Of the 23 rare or unique FANCM PTVs, 15 have not been previously reported. We provide here the initial spectrum of FANCM PTVs in European breast cancer cases.

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.

Recurrent moderate-risk mutations in Finnish breast and ovarian cancer patients

Mutations in BRCA1 and BRCA2 genes predispose to breast and ovarian cancer (BC/OC) with a high lifetime risk, whereas mutations in PALB2, CHEK2, ATM, FANCM, RAD51C and RAD51D genes cause a moderately elevated risk. In the Finnish population, recurrent mutations have been identified in all of these genes, the latest being CHEK2 c.319+2T>A and c.444+1G>A. By genotyping 3,156 cases and 2,089 controls, we estimated the frequencies of CHEK2 c.319+2T>A and c.444+1G>A in Finnish BC patients. CHEK2 c.319+2T>A was detected in 0.7% of the patients, and it was associated with a high risk of BC in the unselected patient group (OR = 5.40 [95% CI 1.58-18.45], p = 0.007) and similarly in the familial patient group. CHEK2 c.444+1G>A was identified in 0.1% of all patients. Additionally, we evaluated the combined prevalence of recurrent moderate-risk gene mutations in 2,487 BC patients, 556 OC patients and 261 BRCA1/2 carriers from 109 families. The overall frequency of the mutations was 13.3% in 1,141 BRCA1/2-negative familial BC patients, 7.5% in 1,727 unselected BC patients and 7.2% in 556 unselected OC patients. At least one moderate-risk gene mutation was found in 12.5% of BRCA1 families and 7.1% of BRCA1 index patients, as well as in 17.0% of BRCA2 families and 11.3% of BRCA2 index patients, and the mutations were associated with an additional risk in the BRCA1/2 index patients (OR = 2.63 [1.15-5.48], p = 0.011). These results support gene panel testing of even multiple members of BC families where several mutations may segregate in different individuals.

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.

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 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.

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.

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.

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.

Prediction of Breast and Prostate Cancer Risks in Male BRCA1 and BRCA2 Mutation Carriers Using Polygenic Risk Scores

Purpose BRCA1/2 mutations increase the risk of breast and prostate cancer in men. Common genetic variants modify cancer risks for female carriers of BRCA1/2 mutations. We investigated-for the first time to our knowledge-associations of common genetic variants with breast and prostate cancer risks for male carriers of BRCA1/ 2 mutations and implications for cancer risk prediction. Materials and Methods We genotyped 1,802 male carriers of BRCA1/2 mutations from the Consortium of Investigators of Modifiers of BRCA1/2 by using the custom Illumina OncoArray. We investigated the combined effects of established breast and prostate cancer susceptibility variants on cancer risks for male carriers of BRCA1/2 mutations by constructing weighted polygenic risk scores (PRSs) using published effect estimates as weights. Results In male carriers of BRCA1/2 mutations, PRS that was based on 88 female breast cancer susceptibility variants was associated with breast cancer risk (odds ratio per standard deviation of PRS, 1.36; 95% CI, 1.19 to 1.56; P = 8.6 × 10-6). Similarly, PRS that was based on 103 prostate cancer susceptibility variants was associated with prostate cancer risk (odds ratio per SD of PRS, 1.56; 95% CI, 1.35 to 1.81; P = 3.2 × 10-9). Large differences in absolute cancer risks were observed at the extremes of the PRS distribution. For example, prostate cancer risk by age 80 years at the 5th and 95th percentiles of the PRS varies from 7% to 26% for carriers of BRCA1 mutations and from 19% to 61% for carriers of BRCA2 mutations, respectively. Conclusion PRSs may provide informative cancer risk stratification for male carriers of BRCA1/2 mutations that might enable these men and their physicians to make informed decisions on the type and timing of breast and prostate cancer risk management.

Identification of 12 new susceptibility loci for different histotypes of epithelial ovarian cancer

To identify common alleles associated with different histotypes of epithelial ovarian cancer (EOC), we pooled data from multiple genome-wide genotyping projects totaling 25,509 EOC cases and 40,941 controls. We identified nine new susceptibility loci for different EOC histotypes: six for serous EOC histotypes (3q28, 4q32.3, 8q21.11, 10q24.33, 18q11.2 and 22q12.1), two for mucinous EOC (3q22.3 and 9q31.1) and one for endometrioid EOC (5q12.3). We then performed meta-analysis on the results for high-grade serous ovarian cancer with the results from analysis of 31,448 BRCA1 and BRCA2 mutation carriers, including 3,887 mutation carriers with EOC. This identified three additional susceptibility loci at 2q13, 8q24.1 and 12q24.31. Integrated analyses of genes and regulatory biofeatures at each locus predicted candidate susceptibility genes, including OBFC1, a new candidate susceptibility gene for low-grade and borderline serous EOC.

FANCM c.5101C>T mutation associates with breast cancer survival and treatment outcome

Breast cancer (BC) is a heterogeneous disease, and different tumor characteristics and genetic variation may affect the clinical outcome. The FANCM c.5101C > T nonsense mutation in the Finnish population associates with increased risk of breast cancer, especially for triple‐negative breast cancer patients. To investigate the association of the mutation with disease prognosis, we studied tumor phenotype, treatment outcome, and patient survival in 3,933 invasive breast cancer patients, including 101 FANCM c.5101C > T mutation carriers and 3,832 non‐carriers. We also examined association of the mutation with nuclear immunohistochemical staining of DNA repair markers in 1,240 breast tumors. The FANCM c.5101C > T mutation associated with poor 10‐year breast cancer‐specific survival (hazard ratio (HR)=1.66, 95% confidence interval (CI) 1.09–2.52, p = 0.018), with a more pronounced survival effect among familial cases (HR = 2.93, 95% CI 1.5–5.76, p = 1.80 × 10⁻³). Poor disease outcome of the carriers was also found among the estrogen receptor (ER) positive subgroup of patients (HR = 1.8, 95% CI 1.09–2.98, p = 0.021). Reduced survival was seen especially among patients who had not received radiotherapy (HR = 3.43, 95% CI 1.6–7.34, p = 1.50 × 10⁻³) but not among radiotherapy treated patients (HR = 1.35, 95% CI 0.82–2.23, p = 0.237). Significant interaction was found between the mutation and radiotherapy (p = 0.040). Immunohistochemical analyses show that c.5101C > T carriers have reduced PAR‐activity. Our results suggest that FANCM c.5101C > T nonsense mutation carriers have a reduced breast cancer survival but postoperative radiotherapy may diminish this survival disadvantage.

What's new?

Variations in DNA repair genes can predispose individuals to breast cancer, with one example being FANCM c.5101C > T, a nonsense mutation in the Fanconi Anemia DNA repair pathway. In previous work, FANCM c.5101C > T was associated with increased breast cancer risk in the Finnish population. Here, the mutation is further shown to be associated with adverse breast cancer outcome. Mutation‐positive Finnish patients exhibited reduced long‐term survival and increased risk of disease recurrence. Survival was worse particularly for patients who were not treated with radiotherapy, indicating that FANCM c.5101C>T may interact with radiotherapy to improve disease outcome in mutation carriers.

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.

RAD51, XRCC3, and XRCC2 mutation screening in Finnish breast cancer families

Majority of the known breast cancer susceptibility genes have a role in DNA repair and the most important high-risk genes BRCA1 and BRCA2 are specifically involved in the homologous recombination repair (HRR) of DNA double-strand breaks. A central player in HRR is RAD51 that binds DNA at the damage site. The RAD51 paralogs RAD51B, RAD51C, RAD51D, XRCC2, and XRCC3 facilitate the binding of RAD51 to DNA. While germline mutations in RAD51C and RAD51D are associated with high ovarian cancer risk and RAD51B polymorphisms with breast cancer, the contribution of RAD51, XRCC3, and XRCC2 is more unclear. To investigate the role of RAD51, XRCC3, and XRCC2 in breast cancer predisposition and to identify putative recurrent founder mutations in the Finnish population where such mutations have been observed in most of the currently known susceptibility genes, we screened 182 familial Finnish breast or ovarian cancer patients for germline variation in the RAD51and XRCC3 genes and 342 patients for variation in XRCC2, with a subset of the patients selected on the basis of decreased RAD51 protein expression on tumors. We also performed haplotype analyses for 1516 breast cancer cases and 1234 controls to assess the common variation in these genes. No pathogenic mutations were detected in any of the genes and the distribution of haplotypes was similar between cases and controls. Our results suggest that RAD51, XRCC3, and XRCC2 do not substantially contribute to breast cancer predisposition in the Finnish population.