CYP2D6 phenotypes are associated with adverse outcomes related to opioid medications

BACKGROUND

Variation in the CYP2D6 gene may affect response to opioids in both poor and ultrarapid metabolizers, but data demonstrating such associations have been mixed, and the impact of variants on toxicity-related symptoms (e.g., nausea) is unclear. Therefore, we examined the association between CYP2D6 phenotype and poor pain control or other adverse symptoms related to the use of opioids in a sample of primary care patients.

MATERIALS AND METHODS

We identified all patients in the Mayo Clinic RIGHT Protocol who were prescribed an opioid medication between July 01, 2013 and June 30, 2015, and categorized patients into three phenotypes: poor, intermediate to extensive, or ultrarapid CYP2D6 metabolizers. We reviewed the electronic health record of these patients for indications of poor pain control or adverse symptoms related to medication use. Associations between phenotype and outcomes were assessed using Chi-square tests and logistic regression.

RESULTS

Overall, 257 (25% of RIGHT Protocol participants) patients received at least one opioid prescription; of these, 40 (15%) were poor metabolizers, 146 (57%) were intermediate to extensive metabolizers, and 71 (28%) were ultrarapid metabolizers. We removed patients that were prescribed a CYP2D6 inhibitor medication (n=38). After adjusting for age and sex, patients with a poor or ultrarapid phenotype were 2.7 times more likely to experience either poor pain control or an adverse symptom related to the prescription compared to patients with an intermediate to extensive phenotype (odds ratio: 2.68; 95% CI: 1.39, 5.17; p=0.003).

CONCLUSION

Our results suggest that >30% of patients with a poor or ultrarapid CYP2D6 phenotype may experience an adverse outcome after being prescribed codeine, tramadol, oxycodone, or hydrocodone. These medications are frequently prescribed for pain relief, and ~39% of the US population is expected to carry one of these phenotypes, suggesting that the population-level impact of these gene-drug interactions could be substantial.

CYP2D6 phenotypes are associated with adverse outcomes related to opioid medications

Background

Variation in the CYP2D6 gene may affect response to opioids in both poor and ultrarapid metabolizers, but data demonstrating such associations have been mixed, and the impact of variants on toxicity-related symptoms (e.g., nausea) is unclear. Therefore, we examined the association between CYP2D6 phenotype and poor pain control or other adverse symptoms related to the use of opioids in a sample of primary care patients.

Materials and methods

We identified all patients in the Mayo Clinic RIGHT Protocol who were prescribed an opioid medication between July 01, 2013 and June 30, 2015, and categorized patients into three phenotypes: poor, intermediate to extensive, or ultrarapid CYP2D6 metabolizers. We reviewed the electronic health record of these patients for indications of poor pain control or adverse symptoms related to medication use. Associations between phenotype and outcomes were assessed using Chi-square tests and logistic regression.

Results

Overall, 257 (25% of RIGHT Protocol participants) patients received at least one opioid prescription; of these, 40 (15%) were poor metabolizers, 146 (57%) were intermediate to extensive metabolizers, and 71 (28%) were ultrarapid metabolizers. We removed patients that were prescribed a CYP2D6 inhibitor medication (n=38). After adjusting for age and sex, patients with a poor or ultrarapid phenotype were 2.7 times more likely to experience either poor pain control or an adverse symptom related to the prescription compared to patients with an intermediate to extensive phenotype (odds ratio: 2.68; 95% CI: 1.39, 5.17; p=0.003).

Conclusion

Our results suggest that >30% of patients with a poor or ultrarapid CYP2D6 phenotype may experience an adverse outcome after being prescribed codeine, tramadol, oxycodone, or hydrocodone. These medications are frequently prescribed for pain relief, and ~39% of the US population is expected to carry one of these phenotypes, suggesting that the population-level impact of these gene–drug interactions could be substantial.

Individual housing-based socioeconomic status predicts risk of accidental falls among adults

PURPOSE

Accidental falls are a major public health concern among people of all ages. Little is known about whether an individual-level housing-based socioeconomic status measure is associated with the risk of accidental falls.

METHODS

Among 12,286 Mayo Clinic Biobank participants residing in Olmsted County, Minnesota, subjects who experienced accidental falls between the biobank enrollment and September 2014 were identified using ICD-9 codes evaluated at emergency departments. HOUSES (HOUsing-based Index of SocioEconomic Status), a socioeconomic status measure based on individual housing features, was also calculated. Cox regression models were utilized to assess the association of the HOUSES (in quartiles) with accidental fall risk.

RESULTS

Seven hundred eleven (5.8%) participants had at least one emergency room visit due to an accidental fall during the study period. Subjects with higher HOUSES were less likely to experience falls in a dose-response manner (hazard ratio: 0.58; 95% confidence interval: 0.44-0.76 for comparing the highest to the lowest quartile). In addition, the HOUSES was positively associated with better health behaviors, social support, and functional status.

CONCLUSIONS

The HOUSES is inversely associated with accidental fall risk requiring emergency care in a dose-response manner. The HOUSES may capture falls-related risk factors through housing features and socioeconomic status-related psychosocial factors.

BRCA2 Hypomorphic Missense Variants Confer Moderate Risks of Breast Cancer

Breast cancer risks conferred by many germline missense variants in the BRCA1 and BRCA2 genes, often referred to as variants of uncertain significance (VUS), have not been established. In this study, associations between 19 BRCA1 and 33 BRCA2 missense substitution variants and breast cancer risk were investigated through a breast cancer case-control study using genotyping data from 38 studies of predominantly European ancestry (41,890 cases and 41,607 controls) and nine studies of Asian ancestry (6,269 cases and 6,624 controls). The BRCA2 c.9104A>C, p.Tyr3035Ser (OR = 2.52; P = 0.04), and BRCA1 c.5096G>A, p.Arg1699Gln (OR = 4.29; P = 0.009) variant were associated with moderately increased risks of breast cancer among Europeans, whereas BRCA2 c.7522G>A, p.Gly2508Ser (OR = 2.68; P = 0.004), and c.8187G>T, p.Lys2729Asn (OR = 1.4; P = 0.004) were associated with moderate and low risks of breast cancer among Asians. Functional characterization of the BRCA2 variants using four quantitative assays showed reduced BRCA2 activity for p.Tyr3035Ser compared with wild-type. Overall, our results show how BRCA2 missense variants that influence protein function can confer clinically relevant, moderately increased risks of breast cancer, with potential implications for risk management guidelines in women with these specific variants. Cancer Res; 77(11); 2789-99. ©2017 AACR.

Genetic modifiers of CHEK2*1100delC-associated breast cancer risk

PURPOSE

CHEK2*1100delC is a founder variant in European populations that confers a two- to threefold increased risk of breast cancer (BC). Epidemiologic and family studies have suggested that the risk associated with CHEK2*1100delC is modified by other genetic factors in a multiplicative fashion. We have investigated this empirically using data from the Breast Cancer Association Consortium (BCAC).

METHODS

Using genotype data from 39,139 (624 1100delC carriers) BC patients and 40,063 (224) healthy controls from 32 BCAC studies, we analyzed the combined risk effects of CHEK2*1100delC and 77 common variants in terms of a polygenic risk score (PRS) and pairwise interaction.

RESULTS

The PRS conferred odds ratios (OR) of 1.59 (95% CI: 1.21-2.09) per standard deviation for BC for CHEK2*1100delC carriers and 1.58 (1.55-1.62) for noncarriers. No evidence of deviation from the multiplicative model was found. The OR for the highest quintile of the PRS was 2.03 (0.86-4.78) for CHEK2*1100delC carriers, placing them in the high risk category according to UK NICE guidelines. The OR for the lowest quintile was 0.52 (0.16-1.74), indicating a lifetime risk close to the population average.

CONCLUSION

Our results confirm the multiplicative nature of risk effects conferred by CHEK2*1100delC and the common susceptibility variants. Furthermore, the PRS could identify carriers at a high lifetime risk for clinical actions.Genet Med advance online publication 06 October 2016.

Increased risk of hospitalization for ultrarapid metabolizers of cytochrome P450 2D6

Background

Cytochrome P450 2D6 (CYP2D6) is responsible for the metabolism of clinically used drugs and other environmental exposures, but it is unclear whether the CYP2D6 phenotype is associated with adverse health outcomes. The aim was to determine the association of CYP2D6 phenotype with the risk of hospitalization or an emergency department (ED) visit among a group of primary care patients.

Methods

In this study, 929 adult patients underwent CYP2D6 testing. The primary outcome was risk of hospitalization or an ED visit from January 2005 through September 2014. CYP2D6 genotypes were interpreted as 1 of 7 clinical phenotypes, from ultrarapid to poor metabolizer, and patients with the extensive metabolizer phenotype were used as the reference group. The hazard ratios (HRs) and 95% confidence intervals (CIs) were estimated for finding the association of CYP2D6 phenotypes with the risk of hospitalization or an ED visit by using Cox proportional hazard models and adjusting for age and sex.

Results

The median age was 49 years (interquartile range, 46–52 years); 74% of patients had 3 or fewer chronic conditions, 285 had at least 1 hospitalization, and 496 had at least 1 ED visit. The risk of hospitalization was higher among patients who were ultrarapid metabolizers compared to extensive metabolizers (47% vs 30%; HR, 1.69; 95% CI, 1.11–2.57), as was the risk of an ED visit (62% vs 49%; HR, 1.50; 95% CI, 1.05–2.14). For poor metabolizers compared to extensive metabolizers, there was no difference in the risk of hospitalization (HR, 0.95; 95% CI, 0.58–1.56), but there was an increase in the risk of an ED visit (HR, 1.38; 95% CI, 0.96–1.98) (the difference was not statistically significant).

Conclusion

We found an increased risk of hospitalization or an ED visit among ultrarapid compared to extensive CYP2D6 metabolizers. Further research identifying the mechanisms of the association and ultimate clinical utility is warranted.

Variability in assigning pathogenicity to incidental findings: insights from LDLR sequence linked to the electronic health record in 1013 individuals

Knowledge of variant pathogenicity is key to implementing genomic medicine. We describe variability between expert reviewers in assigning pathogenicity to sequence variants in LDLR, the causal gene in the majority of cases of familial hypercholesterolemia. LDLR was sequenced on the Illumina HiSeq platform (average read depth >200 × ) in 1013 Mayo Biobank participants recruited from 2012 to 2013. Variants with a minor allele frequency (MAF) <5% predicted to be functional or referenced in HGMD (Human Gene Mutation Database) or NCBI-ClinVar databases were reviewed. To assign pathogenicity, variant frequency in population data sets, computational predictions, reported observations and patient-level data including electronic health record-based post hoc phenotyping were leveraged. Of 178 LDLR variants passing quality control, 25 were selected for independent review using either an in-house protocol or a disease/gene-specific semi-quantitative framework based on the American College of Medical Genetics and Genomics-recommended lines of evidence. NCBI-ClinVar included interpretations for all queried variants with 74% (14/19) of variants with >1 submitter showing inconsistency in classification and 26% (5/19) appearing with conflicting clinical actionability. The discordance rate (one-step level of disagreement out of five classes in variant interpretation) between the reviewers was 40% (10/25). Two LDLR variants were independently deemed clinically actionable and returnable. Interpretation of LDLR variants was often discordant among ClinVar submitters and between expert reviewers. A quantitative approach based on strength of each predefined criterion in the context of specific genes and phenotypes may yield greater consistency between different reviewers.

PALB2, CHEK2 and ATM rare variants and cancer risk: data from COGS

BACKGROUND

The rarity of mutations in PALB2, CHEK2 and ATM make it difficult to estimate precisely associated cancer risks. Population-based family studies have provided evidence that at least some of these mutations are associated with breast cancer risk as high as those associated with rare BRCA2 mutations. We aimed to estimate the relative risks associated with specific rare variants in PALB2, CHEK2 and ATM via a multicentre case-control study.

METHODS

We genotyped 10 rare mutations using the custom iCOGS array: PALB2 c.1592delT, c.2816T>G and c.3113G>A, CHEK2 c.349A>G, c.538C>T, c.715G>A, c.1036C>T, c.1312G>T, and c.1343T>G and ATM c.7271T>G. We assessed associations with breast cancer risk (42 671 cases and 42 164 controls), as well as prostate (22 301 cases and 22 320 controls) and ovarian (14 542 cases and 23 491 controls) cancer risk, for each variant.

RESULTS

For European women, strong evidence of association with breast cancer risk was observed for PALB2 c.1592delT OR 3.44 (95% CI 1.39 to 8.52, p=7.1×10-5), PALB2 c.3113G>A OR 4.21 (95% CI 1.84 to 9.60, p=6.9×10-8) and ATM c.7271T>G OR 11.0 (95% CI 1.42 to 85.7, p=0.0012). We also found evidence of association with breast cancer risk for three variants in CHEK2, c.349A>G OR 2.26 (95% CI 1.29 to 3.95), c.1036C>T OR 5.06 (95% CI 1.09 to 23.5) and c.538C>T OR 1.33 (95% CI 1.05 to 1.67) (p≤0.017). Evidence for prostate cancer risk was observed for CHEK2 c.1343T>G OR 3.03 (95% CI 1.53 to 6.03, p=0.0006) for African men and CHEK2 c.1312G>T OR 2.21 (95% CI 1.06 to 4.63, p=0.030) for European men. No evidence of association with ovarian cancer was found for any of these variants.

CONCLUSIONS

This report adds to accumulating evidence that at least some variants in these genes are associated with an increased risk of breast cancer that is clinically important.

Integrating Pharmacogenomics into Clinical Practice: Promise vs Reality

BACKGROUND

Limited information is available regarding primary care clinicians' response to pharmacogenomic clinical decision support (PGx-CDS) alerts integrated in the electronic health record.

METHODS

In February 2015, 159 clinicians in the Mayo Clinic primary care practice were sent e-mail surveys to understand their perspectives on the implementation and use of pharmacogenomic testing in their clinical practice. Surveys assessed how the clinicians felt about pharmacogenomics and whether they thought electronic PGx-CDS alerts were useful. Information was abstracted on the number of CDS alerts the clinicians received between October 2013 and the date their survey was returned. CDS alerts were grouped into 2 categories: the alert recommended caution using the prescription, or the alert recommended an alternate prescription. Finally, data were abstracted regarding whether the clinician changed their prescription in response to the alert recommendation.

RESULTS

The survey response rate was 57% (n = 90). Overall, 52% of the clinicians did not expect to use or did not know whether they would use pharmacogenomic information in their future prescribing practices. Additionally, 53% of the clinicians felt that the alerts were confusing, irritating, frustrating, or that it was difficult to find additional information. Finally, only 30% of the clinicians that received a CDS alert changed their prescription to an alternative medication.

CONCLUSIONS

Our results suggest a lack of clinician comfort with integration of pharmacogenomic data into primary care. Further efforts to refine PGx-CDS alerts to make them as useful and user-friendly as possible are needed to improve clinician satisfaction with these new tools.

Multidisciplinary model to implement pharmacogenomics at the point of care

Purpose

Despite potential clinical benefits, implementation of pharmacogenomics (PGx) faces many technical and clinical challenges. These challenges can be overcome by a comprehensive and systematic implementation model.

Methods

The development and implementation of PGx was organized into eight interdependent components addressing resources, governance, clinical practice, education, testing, knowledge translation, clinical decision support (CDS) and maintenance. Several aspects of the implementation were assessed including adherence to the model, production of PGx-CDS interventions and access to educational resources.

Results

Between 8/2012 and 6/2015, 21 specific drug-gene interactions were reviewed and 18 of them were implemented in the electronic medical record as PGx-CDS interventions. There was complete adherence to the model with variable production time (98 to 392 days) and delay time (0 to 148 days). The implementation impacted approximately 1247 unique providers and 3788 unique patients. A total of 11 educational resources complementary to the drug-gene interactions and 5 modules specific for pharmacists were developed and implemented.

Conclusion

A comprehensive operational model can support PGx implementation into routine prescribing. Institutions can use this model as a roadmap to support similar efforts. However, we also identified challenges that will require major multidisciplinary and multi-institutional efforts to make PGx a universal reality.

Fine-scale mapping of 8q24 locus identifies multiple independent risk variants for breast cancer

Previous genome-wide association studies among women of European ancestry identified two independent breast cancer susceptibility loci represented by single nucleotide polymorphisms (SNPs) rs13281615 and rs11780156 at 8q24. A fine-mapping study across 2.06 Mb (chr8:127,561,724-129,624,067, hg19) in 55,540 breast cancer cases and 51,168 controls within the Breast Cancer Association Consortium was conducted. Three additional independent association signals in women of European ancestry, represented by rs35961416 (OR = 0.95, 95% CI = 0.93-0.97, conditional p = 5.8 × 10(-6) ), rs7815245 (OR = 0.94, 95% CI = 0.91-0.96, conditional p = 1.1 × 10(-6) ) and rs2033101 (OR = 1.05, 95% CI = 1.02-1.07, conditional p = 1.1 × 10(-4) ) were found. Integrative analysis using functional genomic data from the Roadmap Epigenomics, the Encyclopedia of DNA Elements project, the Cancer Genome Atlas and other public resources implied that SNPs rs7815245 in Signal 3, and rs1121948 in Signal 5 (in linkage disequilibrium with rs11780156, r(2)  = 0.77), were putatively functional variants for two of the five independent association signals. The results highlighted multiple 8q24 variants associated with breast cancer susceptibility in women of European ancestry.

Functional mechanisms underlying pleiotropic risk alleles at the 19p13.1 breast-ovarian cancer susceptibility locus

A locus at 19p13 is associated with breast cancer (BC) and ovarian cancer (OC) risk. Here we analyse 438 SNPs in this region in 46,451 BC and 15,438 OC cases, 15,252 BRCA1 mutation carriers and 73,444 controls and identify 13 candidate causal SNPs associated with serous OC (P=9.2 × 10(-20)), ER-negative BC (P=1.1 × 10(-13)), BRCA1-associated BC (P=7.7 × 10(-16)) and triple negative BC (P-diff=2 × 10(-5)). Genotype-gene expression associations are identified for candidate target genes ANKLE1 (P=2 × 10(-3)) and ABHD8 (P<2 × 10(-3)). Chromosome conformation capture identifies interactions between four candidate SNPs and ABHD8, and luciferase assays indicate six risk alleles increased transactivation of the ADHD8 promoter. Targeted deletion of a region containing risk SNP rs56069439 in a putative enhancer induces ANKLE1 downregulation; and mRNA stability assays indicate functional effects for an ANKLE1 3'-UTR SNP. Altogether, these data suggest that multiple SNPs at 19p13 regulate ABHD8 and perhaps ANKLE1 expression, and indicate common mechanisms underlying breast and ovarian cancer risk.

Aspirin use and the risk of cholangiocarcinoma

Whether aspirin use is protective against cholangiocarcinoma (CCA) remains unclear. We determined the association between aspirin use and other risk factors for each CCA subtype individually. In a hospital-based case-control study, 2395 CCA cases (1169 intrahepatic, 995 perihilar, and 231 distal) seen at the Mayo Clinic, Rochester, MN, from 2000 through 2014 were enrolled. Controls selected from the Mayo Clinic Biobank were matched two to one with cases by age, sex, race, and residence (n = 4769). Associations between aspirin use, other risk factors, and CCA risk were determined. Aspirin was used by 591 (24.7%) CCA cases and 2129 (44.6%) controls. There was a significant inverse association of aspirin use with all CCA subtypes, with adjusted odds ratios (AORs) of 0.35 (95% confidence interval [CI], 0.29-0.42), 0.34 (95% CI 0.27-0.42), and 0.29 (95% CI 0.19-0.44) for intrahepatic, perihilar, and distal CCA, respectively (P < 0.001 for all). Primary sclerosing cholangitis was more strongly associated with perihilar (AOR = 453, 95% CI 104-999) than intrahepatic (AOR = 93.4, 95% CI 27.1-322) or distal (AOR = 34.0, 95% CI 3.6-323) CCA, whereas diabetes was more associated with distal (AOR = 4.2, 95% CI 2.5-7.0) than perihilar (AOR = 2.9, 95% CI 2.2-3.8) or intrahepatic (AOR = 2.5, 95% CI 2.0-3.2) CCA. Cirrhosis not related to primary sclerosing cholangitis was associated with both intrahepatic and perihilar CCA, with similar AORs of 14. Isolated inflammatory bowel disease without primary sclerosing cholangitis was not associated with any CCA subtype.

CONCLUSIONS

Aspirin use was significantly associated with a 2.7-fold to 3.6-fold decreased risk for the three CCA subtypes; our study demonstrates that individual risk factors confer risk of different CCA subtypes to different extents. (Hepatology 2016;64:785-796).

An intergenic risk locus containing an enhancer deletion in 2q35 modulates breast cancer risk by deregulating IGFBP5 expression

Breast cancer is the most diagnosed malignancy and the second leading cause of cancer mortality in females. Previous association studies have identified variants on 2q35 associated with the risk of breast cancer. To identify functional susceptibility loci for breast cancer, we interrogated the 2q35 gene desert for chromatin architecture and functional variation correlated with gene expression. We report a novel intergenic breast cancer risk locus containing an enhancer copy number variation (enCNV; deletion) located approximately 400Kb upstream to IGFBP5, which overlaps an intergenic ERα-bound enhancer that loops to the IGFBP5 promoter. The enCNV is correlated with modified ERα binding and monoallelic-repression of IGFBP5 following oestrogen treatment. We investigated the association of enCNV genotype with breast cancer in 1,182 cases and 1,362 controls, and replicate our findings in an independent set of 62,533 cases and 60,966 controls from 41 case control studies and 11 GWAS. We report a dose-dependent inverse association of 2q35 enCNV genotype (percopy OR = 0.68 95%CI 0.55-0.83, P = 0.0002; replication OR = 0.77 95% CI 0.73-0.82, P = 2.1 × 10-19) and identify 13 additional linked variants (r2 > 0.8) in the 20Kb linkage block containing the enCNV (P = 3.2 × 10-15 - 5.6 × 10-17). These associations were independent of previously reported 2q35 variants, rs13387042/rs4442975 and rs16857609, and were stronger for ER-positive than ER-negative disease. Together, these results suggest that 2q35 breast cancer risk loci may be mediating their effect through IGFBP5.

Genetically Predicted Body Mass Index and Breast Cancer Risk: Mendelian Randomization Analyses of Data from 145,000 Women of European Descent

BACKGROUND

Observational epidemiological studies have shown that high body mass index (BMI) is associated with a reduced risk of breast cancer in premenopausal women but an increased risk in postmenopausal women. It is unclear whether this association is mediated through shared genetic or environmental factors.

METHODS

We applied Mendelian randomization to evaluate the association between BMI and risk of breast cancer occurrence using data from two large breast cancer consortia. We created a weighted BMI genetic score comprising 84 BMI-associated genetic variants to predicted BMI. We evaluated genetically predicted BMI in association with breast cancer risk using individual-level data from the Breast Cancer Association Consortium (BCAC) (cases  =  46,325, controls  =  42,482). We further evaluated the association between genetically predicted BMI and breast cancer risk using summary statistics from 16,003 cases and 41,335 controls from the Discovery, Biology, and Risk of Inherited Variants in Breast Cancer (DRIVE) Project. Because most studies measured BMI after cancer diagnosis, we could not conduct a parallel analysis to adequately evaluate the association of measured BMI with breast cancer risk prospectively.

RESULTS

In the BCAC data, genetically predicted BMI was found to be inversely associated with breast cancer risk (odds ratio [OR]  =  0.65 per 5 kg/m2 increase, 95% confidence interval [CI]: 0.56-0.75, p = 3.32 × 10-10). The associations were similar for both premenopausal (OR   =   0.44, 95% CI:0.31-0.62, p  =  9.91 × 10-8) and postmenopausal breast cancer (OR  =  0.57, 95% CI: 0.46-0.71, p  =  1.88 × 10-8). This association was replicated in the data from the DRIVE consortium (OR  =  0.72, 95% CI: 0.60-0.84, p   =   1.64 × 10-7). Single marker analyses identified 17 of the 84 BMI-associated single nucleotide polymorphisms (SNPs) in association with breast cancer risk at p < 0.05; for 16 of them, the allele associated with elevated BMI was associated with reduced breast cancer risk.

CONCLUSIONS

BMI predicted by genome-wide association studies (GWAS)-identified variants is inversely associated with the risk of both pre- and postmenopausal breast cancer. The reduced risk of postmenopausal breast cancer associated with genetically predicted BMI observed in this study differs from the positive association reported from studies using measured adult BMI. Understanding the reasons for this discrepancy may reveal insights into the complex relationship of genetic determinants of body weight in the etiology of breast cancer.

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

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

RAD51B in Familial Breast Cancer

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

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

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

OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Identification of four novel susceptibility loci for oestrogen receptor negative breast cancer

Common variants in 94 loci have been associated with breast cancer including 15 loci with genome-wide significant associations (P<5 × 10(-8)) with oestrogen receptor (ER)-negative breast cancer and BRCA1-associated breast cancer risk. In this study, to identify new ER-negative susceptibility loci, we performed a meta-analysis of 11 genome-wide association studies (GWAS) consisting of 4,939 ER-negative cases and 14,352 controls, combined with 7,333 ER-negative cases and 42,468 controls and 15,252 BRCA1 mutation carriers genotyped on the iCOGS array. We identify four previously unidentified loci including two loci at 13q22 near KLF5, a 2p23.2 locus near WDR43 and a 2q33 locus near PPIL3 that display genome-wide significant associations with ER-negative breast cancer. In addition, 19 known breast cancer risk loci have genome-wide significant associations and 40 had moderate associations (P<0.05) with ER-negative disease. Using functional and eQTL studies we implicate TRMT61B and WDR43 at 2p23.2 and PPIL3 at 2q33 in ER-negative breast cancer aetiology. All ER-negative loci combined account for ∼11% of familial relative risk for ER-negative disease and may contribute to improved ER-negative and BRCA1 breast cancer risk prediction.

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

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

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

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

Incorporating a Genetic Risk Score Into Coronary Heart Disease Risk Estimates: Effect on Low-Density Lipoprotein Cholesterol Levels (the MI-GENES Clinical Trial)

BACKGROUND

Whether knowledge of genetic risk for coronary heart disease (CHD) affects health-related outcomes is unknown. We investigated whether incorporating a genetic risk score (GRS) in CHD risk estimates lowers low-density lipoprotein cholesterol (LDL-C) levels.

METHODS AND RESULTS

Participants (n=203, 45-65 years of age, at intermediate risk for CHD, and not on statins) were randomly assigned to receive their 10-year probability of CHD based either on a conventional risk score (CRS) or CRS + GRS ((+)GRS). Participants in the (+)GRS group were stratified as having high or average/low GRS. Risk was disclosed by a genetic counselor followed by shared decision making regarding statin therapy with a physician. We compared the primary end point of LDL-C levels at 6 months and assessed whether any differences were attributable to changes in dietary fat intake, physical activity levels, or statin use. Participants (mean age, 59.4±5 years; 48% men; mean 10-year CHD risk, 8.5±4.1%) were allocated to receive either CRS (n=100) or (+)GRS (n=103). At the end of the study period, the (+)GRS group had a lower LDL-C than the CRS group (96.5±32.7 versus 105.9±33.3 mg/dL; P=0.04). Participants with high GRS had lower LDL-C levels (92.3±32.9 mg/dL) than CRS participants (P=0.02) but not participants with low GRS (100.9±32.2 mg/dL; P=0.18). Statins were initiated more often in the (+)GRS group than in the CRS group (39% versus 22%, P<0.01). No significant differences in dietary fat intake and physical activity levels were noted.

CONCLUSIONS

Disclosure of CHD risk estimates that incorporated genetic risk information led to lower LDL-C levels than disclosure of CHD risk based on conventional risk factors alone.

CLINICAL TRIAL REGISTRATION

URL: http://www.clinicaltrials.gov. Unique identifier: NCT01936675.

Improvement in Cardiovascular Risk Prediction with Electronic Health Records

The aim of this study was to compare the QRISKII, an electronic health data-based risk score, to the Framingham Risk Score (FRS) and atherosclerotic cardiovascular disease (ASCVD) score. Risk estimates were calculated for a cohort of 8783 patients, and the patients were followed up from November 29, 2012, through June 1, 2015, for a cardiovascular disease (CVD) event. During follow-up, 246 men and 247 women had a CVD event. Cohen's kappa statistic for the comparison of the QRISKII and FRS was 0.22 for men and 0.23 for women, with the QRISKII classifying more patients in the higher-risk groups. The QRISKII and ASCVD were more similar with kappa statistics of 0.49 for men and 0.51 for women. The QRISKII shows increased discrimination with area under the curve (AUC) statistics of 0.65 and 0.71, respectively, compared to the FRS (0.59 and 0.66) and ASCVD (0.63 and 0.69). These results demonstrate that incorporating additional data from the electronic health record (EHR) may improve CVD risk stratification.

Genetic predisposition to ductal carcinoma in situ of the breast

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSION

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