Exome sequencing identifies breast cancer susceptibility genes and defines the contribution of coding variants to breast cancer risk.

Linkage and candidate gene studies have identified several breast cancer susceptibility genes, but the overall contribution of coding variation to breast cancer is unclear. To evaluate the role of rare coding variants more comprehensively, we performed a meta-analysis across three large whole-exome sequencing datasets, containing 26,368 female cases and 217,673 female controls. Burden tests were performed for protein-truncating and rare missense variants in 15,616 and 18,601 genes, respectively. Associations between protein-truncating variants and breast cancer were identified for the following six genes at exome-wide significance (P < 2.5 × 10-6): the five known susceptibility genes ATM, BRCA1, BRCA2, CHEK2 and PALB2, together with MAP3K1. Associations were also observed for LZTR1, ATR and BARD1 with P < 1 × 10-4. Associations between predicted deleterious rare missense or protein-truncating variants and breast cancer were additionally identified for CDKN2A at exome-wide significance. The overall contribution of coding variants in genes beyond the previously known genes is estimated to be small.

A Comprehensive Analysis of 21 Actionable Pharmacogenes in the Spanish Population: From Genetic Characterisation to Clinical Impact.

The implementation of pharmacogenetics (PGx) is a main milestones of precision medicine nowadays in order to achieve safer and more effective therapies. Nevertheless, the implementation of PGx diagnostics is extremely slow and unequal worldwide, in part due to a lack of ethnic PGx information. We analysed genetic data from 3006 Spanish individuals obtained by different high-throughput (HT) techniques. Allele frequencies were determined in our population for the main 21 actionable PGx genes associated with therapeutical changes. We found that 98% of the Spanish population harbours at least one allele associated with a therapeutical change and, thus, there would be a need for a therapeutical change in a mean of 3.31 of the 64 associated drugs. We also identified 326 putative deleterious variants that were not previously related with PGx in 18 out of the 21 main PGx genes evaluated and a total of 7122 putative deleterious variants for the 1045 PGx genes described. Additionally, we performed a comparison of the main HT diagnostic techniques, revealing that after whole genome sequencing, genotyping with the PGx HT array is the most suitable solution for PGx diagnostics. Finally, all this information was integrated in the Collaborative Spanish Variant Server to be available to and updated by the scientific community.

POLRMT as a Novel Susceptibility Gene for Cardiotoxicity in Epirubicin Treatment of Breast Cancer Patients.

Anthracyclines are among the most used chemotherapeutic agents in breast cancer (BC). However their use is hampered by anthracycline-induced cardiotoxicity (AIC). The currently known clinical and genetic risk factors do not fully explain the observed inter-individual variability and only have a limited ability to predict which patients are more likely to develop this severe toxicity. To identify novel predictive genes, we conducted a two-stage genome-wide association study in epirubicin-treated BC patients. In the discovery phase, we genotyped over 700,000 single nucleotide variants in a cohort of 227 patients. The most interesting finding was rs62134260, located 4kb upstream of POLRMT (OR = 5.76, P = 2.23 × 10-5). We replicated this association in a validation cohort of 123 patients (P = 0.021). This variant regulates the expression of POLRMT, a gene that encodes a mitochondrial DNA-directed RNA polymerase, responsible for mitochondrial gene expression. Individuals harbouring the risk allele had a decreased expression of POLRMT in heart tissue that may cause an impaired capacity to maintain a healthy mitochondrial population in cardiomyocytes under stressful conditions, as is treatment with epirubicin. This finding suggests a novel molecular mechanism involved in the development of AIC and may improve our ability to predict patients who are at risk.

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

Association Between ABCB1 Genetic Variants and Persistent Chemotherapy-Induced Alopecia in Women With Breast Cancer.

Importance: Persistent chemotherapy-induced alopecia (pCIA) has been recently described in patients with breast cancer and in its most severe form occurs in up to 10% of these patients. Genetic risk factors associated with pCIA have not been adequately explored. Objective: To identify genetic variants associated with pCIA. Design, Setting, and Participants: In this genetic association study, 215 women with breast cancer treated with docetaxel-based chemotherapy with a follow-up of 1.5 to 10 years after the end of the treatment were recruited retrospectively through 3 hospital oncology units across Spain between 2005 and 2018. Severe pCIA was defined as lack of scalp hair recovery (Common Terminology Criteria for Adverse Events, version 3.0, grade 2) 18 months or more after the end of treatment. Patients with grade 2 pCIA were selected as cases, and those with no sign of residual alopecia 12 months after the end of docetaxel treatment were selected as controls. A genome-wide association study in a discovery phase was conducted, and logistic regression was used to identify variants associated with the risk to develop this adverse effect. The validity of the association was addressed through a replication phase. Exposures: Docetaxel-based chemotherapy. Main Outcomes and Measures: Genotypes of single-nucleotide variants associated with pCIA. Results: In total, 215 women with breast cancer (median age, 51.6 years; interquartile range, 44-60 years) were recruited (173 patients for the discovery phase and 42 patients for the replication phase). In the discovery phase, ABCB1 genetic variants were associated with risk to develop pCIA. In particular, single-nucleotide variation rs1202179, a regulatory variant located in an enhancer element that interacts with the ABCB1 promoter, was associated with the occurrence of pCIA. This finding was validated in the replication cohort (combined odds ratio, 4.05; 95% CI, 2.46-6.67; P = 3.946 × 10-8). This variant is associated with ABCB1 mRNA expression, and the risk allele was associated with decreased ABCB1 expression levels (P = 1.64 × 10-20). Conclusions and Relevance: This is the first study, to our knowledge, that identifies an association between a regulatory variant in the ABCB1 gene and the occurrence of pCIA in patients with breast cancer who were treated with docetaxel-based therapies. This finding suggests an important insight into the biological mechanisms underlying pCIA and opens the opportunity to explore personalized treatment of these patients.

Recurrent Germline DLST Mutations in Individuals with Multiple Pheochromocytomas and Paragangliomas.

Pheochromocytomas and paragangliomas (PPGLs) provide some of the clearest genetic evidence for the critical role of metabolism in the tumorigenesis process. Approximately 40% of PPGLs are caused by driver germline mutations in 16 known susceptibility genes, and approximately half of these genes encode members of the tricarboxylic acid (TCA) cycle. Taking as a starting point the involvement of the TCA cycle in PPGL development, we aimed to identify unreported mutations that occurred in genes involved in this key metabolic pathway and that could explain the phenotypes of additional individuals who lack mutations in known susceptibility genes. To accomplish this, we applied a targeted sequencing of 37 TCA-cycle-related genes to DNA from 104 PPGL-affected individuals with no mutations in the major known predisposing genes. We also performed omics-based analyses, TCA-related metabolite determination, and 13C5-glutamate labeling assays. We identified five germline variants affecting DLST in eight unrelated individuals (∼7%); all except one were diagnosed with multiple PPGLs. A recurrent variant, c.1121G>A (p.Gly374Glu), found in four of the eight individuals triggered accumulation of 2-hydroxyglutarate, both in tumors and in a heterologous cell-based assay designed to functionally evaluate DLST variants. p.Gly374Glu-DLST tumors exhibited loss of heterozygosity, and their methylation and expression profiles are similar to those of EPAS1-mutated PPGLs; this similarity suggests a link between DLST disruption and pseudohypoxia. Moreover, we found positive DLST immunostaining exclusively in tumors carrying TCA-cycle or EPAS1 mutations. In summary, this study reveals DLST as a PPGL-susceptibility gene and further strengthens the relevance of the TCA cycle in PPGL development.

Exome array analysis identifies ETFB as a novel susceptibility gene for anthracycline-induced cardiotoxicity in cancer patients.

PURPOSE: Anthracyclines are widely used chemotherapeutic drugs that can cause progressive and irreversible cardiac damage and fatal heart failure. Several genetic variants associated with anthracycline-induced cardiotoxicity (AIC) have been identified, but they explain only a small proportion of the interindividual differences in AIC susceptibility. METHODS: In this study, we evaluated the association of low-frequency variants with risk of chronic AIC using the Illumina HumanExome BeadChip array in a discovery cohort of 61 anthracycline-treated breast cancer patients with replication in a second independent cohort of 83 anthracycline-treated pediatric cancer patients, using gene-based tests (SKAT-O). RESULTS: The most significant associated gene in the discovery cohort was ETFB (electron transfer flavoprotein beta subunit) involved in mitochondrial ß-oxidation and ATP production (P = 4.16 × 10-4) and this association was replicated in an independent set of anthracycline-treated cancer patients (P = 2.81 × 10-3). Within ETFB, we found that the missense variant rs79338777 (p.Pro52Leu; c.155C > T) made the greatest contribution to the observed gene association and it was associated with increased risk of chronic AIC in the two cohorts separately and when combined (OR 9.00, P = 1.95 × 10-4, 95% CI 2.83-28.6). CONCLUSIONS: We identified and replicated a novel gene, ETFB, strongly associated with chronic AIC independently of age at tumor onset and related to anthracycline-mediated mitochondrial dysfunction. Although experimental verification and further studies in larger patient cohorts are required to confirm our finding, we demonstrated that exome array data analysis represents a valuable strategy to identify novel genes contributing to the susceptibility to chronic AIC.

Association analysis between breast cancer genetic variants and mammographic density in a large population-based study (Determinants of Density in Mammographies in Spain) identifies susceptibility loci in TOX3 gene.

BACKGROUND: Mammographic density (MD) is regarded as an intermediate phenotype in breast cancer development. This association study investigated the influence of 14 breast cancer susceptibility loci identified through previous genome-wide association studies on MD among the participants in the "Determinants of Density in Mammographies in Spain" (DDM-Spain) study. METHODS: Our study covered a total of 3348 Caucasian women aged 45-68years, recruited from seven Spanish breast cancer screening centres having DNA available. Mammographic density was blindly assessed by a single reader using a semiquantitative scale. Ordinal logistic models, adjusted for age, body mass index and menopausal status, were used to estimate the association between each genotype and MD. RESULTS: Evidence of association with MD was found for variant rs3803662 (TOX3) (Odds Ratio (OR)=1.13, 95% Confidence Interval (CI)=1.03-1.25), and marginal evidence of association for susceptibility loci rs3817198 (LSP1) (OR=1.09, 95% CI=1.00-1.20) and rs2981582 (FGFR2) (OR=0.92, 95% CI=0.84-1.01). Two other loci were associated with MD solely among pre-menopausal women, namely, rs4973768 (SLC4A7) (OR=0.83, 95% CI=0.70-1.00) and rs4415084 (MEPS30) (OR=1.22, 95% CI=1.00-1.49). CONCLUSIONS: Our findings lend some support to the hypothesis which links these susceptibility loci to MD.

Polymorphisms in the hypoxia inducible factor 1-alpha and the impact on the prognosis of early stages of oral cancer.

BACKGROUND: Hypoxia-inducible factor-1 (HIF-1) is the key regulator of cellular responses to hypoxia and presumably plays a central role in the control of tumor growth. Polymorphisms or mutations increasing its activity and stability in vitro under normoxia have recently been identified. In this study, we aimed to investigate the effect of C1772T and G1790A single nucleotide polymorphisms (SNPs), located within the exon 12 of HIF-1alpha on the prognosis of early stages of oral squamous cell carcinoma (OSCC). METHODS: The frequency of C1772T and G1790A polymorphisms was determined by PCR-RFLP in 139 DNA samples from healthy volunteers and 74 patients with surgically treated T1/2 N0 OSCC. The impact of HIF-1alpha SNPs on tumor size, invasive depth, pathological features, and histological grade was studied. Correlations between genotype and relapse and/or disease-specific survival were evaluated by Kaplan-Meier analysis and log-rank test. RESULTS: Concerning G1790A SNP, the frequencies of GA heterozygous and AA variant homozygous genotypes were significantly higher in patients than in healthy volunteers (32.8% vs. 6.5% and 4.7% vs. none, respectively) (P < .0001). Also, the presence of the variant allele A was associated to disease-relapse (P = .02) and shorter disease-free survival (P = .04). The genotype distribution of C1772T did not diverge between patients and healthy subjects, and no differences were observed with respect to disease-free or overall survival. CONCLUSIONS: Our results suggest that G1790A polymorphism in the HIF-1alpha gene might confer susceptibility to OSCC and could be a marker of disfavorable prognosis at early stages.

Combination of polymorphisms within 5' and 3' untranslated regions of thymidylate synthase gene modulates survival in 5 fluorouracil-treated colorectal cancer patients.

In the present study we explored the effect of three polymorphisms of the TS gene on overall and progression- free survival of colorectal cancer (CRC) patients subjected to 5FU chemotherapy. A 28 bp variable number of tandem repeats (VNTR), a G/C single nucleotide polymorphism (SNP), and a deletion of 6 bp at position 1494 were studied. The possible combined effect of these DNA polymorphisms on the clinical outcome of patients was also evaluated. A retrospective study was carried out on paraffin-embedded sections from 113 patients diagnosed of advanced CRC. TS genotyping methods were polymerase chain reaction (PCR) for VNTR and PCR, followed by restriction length fragment polymorphism (PCR-RFLP) for SNP and ins/del 6 bp. To study the combined effect of TS polymorphisms, four categories were defined accordingly to the level of expression attributed to SNP and ins/del 6 bp genotypes: C&allele 6-, C&6+/6+, G&allele6- and G&6+/6+. VNTR and ins/del 6 bp genotypes varied with tumour anatomical site: 2R/2R genotype was rare in left-sided tumours (7.0% vs. 26.3% of right-sided and 24.1% of rectal cancers; P<0.01), where the variant allele 6- was very frequent (69.0%). Instead, most patients with right-sided tumours were wild-type homozygous 6+/6+ (63.9%) (P<0.01). Heterozygous 6+/6- genotype was more frequent among tumours classified as C (50.0%) and D (76.5%) Dukes stages (P=0.05). None of the studied polymorphisms alone affected overall or progression-free survival (PFS). C&6+/6+ and G&6+/6+ combined genotypes were respectively associated to the best and worst PFS (P=0.03 when compared with each other), while combinations carrying the allele 6- determined an intermediate evolution that might be indicative of a variable response to chemotherapy. The rate of Dukes B stage tumours was unexpectedly high (59.1%) among patients with the unfavourable G&6+/6+ combination. In our study the combination of high TS expression genotypes G&6+/6+ identifies a group of high risk within CRC patients treated with 5FU.