Multivariable regression analysis of febrile neutropenia occurrence in early breast cancer patients receiving chemotherapy assessing patient-related, chemotherapy-related and genetic risk factors.

BACKGROUND: Febrile neutropenia (FN) is common in breast cancer patients undergoing chemotherapy. Risk factors for FN have been reported, but risk models that include genetic variability have yet to be described. This study aimed to evaluate the predictive value of patient-related, chemotherapy-related, and genetic risk factors. METHODS: Data from consecutive breast cancer patients receiving chemotherapy with 4-6 cycles of fluorouracil, epirubicin, and cyclophosphamide (FEC) or three cycles of FEC and docetaxel were retrospectively recorded. Multivariable logistic regression was carried out to assess risk of FN during FEC chemotherapy cycles. RESULTS: Overall, 166 (16.7%) out of 994 patients developed FN. Significant risk factors for FN in any cycle and the first cycle were lower platelet count (OR = 0.78 [0.65; 0.93]) and haemoglobin (OR = 0.81 [0.67; 0.98]) and homozygous carriers of the rs4148350 variant T-allele (OR = 6.7 [1.04; 43.17]) in MRP1. Other significant factors for FN in any cycle were higher alanine aminotransferase (OR = 1.02 [1.01; 1.03]), carriers of the rs246221 variant C-allele (OR = 2.0 [1.03; 3.86]) in MRP1 and the rs351855 variant C-allele (OR = 2.48 [1.13; 5.44]) in FGFR4. Lower height (OR = 0.62 [0.41; 0.92]) increased risk of FN in the first cycle. CONCLUSIONS: Both established clinical risk factors and genetic factors predicted FN in breast cancer patients. Prediction was improved by adding genetic information but overall remained limited. Internal validity was satisfactory. Further independent validation is required to confirm these findings

Intermediate Molecular Phenotypes to Identify Genetic Markers of Anthracycline-Induced Cardiotoxicity Risk.

Cardiotoxicity due to anthracyclines (CDA) affects cancer patients, but we cannot predict who may suffer from this complication. CDA is a complex trait with a polygenic component that is mainly unidentified. We propose that levels of intermediate molecular phenotypes (IMPs) in the myocardium associated with histopathological damage could explain CDA susceptibility, so variants of genes encoding these IMPs could identify patients susceptible to this complication. Thus, a genetically heterogeneous cohort of mice (n = 165) generated by backcrossing were treated with doxorubicin and docetaxel. We quantified heart fibrosis using an Ariol slide scanner and intramyocardial levels of IMPs using multiplex bead arrays and QPCR. We identified quantitative trait loci linked to IMPs (ipQTLs) and cdaQTLs via linkage analysis. In three cancer patient cohorts, CDA was quantified using echocardiography or Cardiac Magnetic Resonance. CDA behaves as a complex trait in the mouse cohort. IMP levels in the myocardium were associated with CDA. ipQTLs integrated into genetic models with cdaQTLs account for more CDA phenotypic variation than that explained by cda-QTLs alone. Allelic forms of genes encoding IMPs associated with CDA in mice, including AKT1, MAPK14, MAPK8, STAT3, CAS3, and TP53, are genetic determinants of CDA in patients. Two genetic risk scores for pediatric patients (n = 71) and women with breast cancer (n = 420) were generated using machine-learning Least Absolute Shrinkage and Selection Operator (LASSO) regression. Thus, IMPs associated with heart damage identify genetic markers of CDA risk, thereby allowing more personalized patient management.J.P.L.’s lab is sponsored by Grant PID2020-118527RB-I00 funded by MCIN/AEI/10.13039/ 501100011039; Grant PDC2021-121735-I00 funded by MCIN/AEI/10.13039/501100011039 and by the “European Union Next Generation EU/PRTR”, the Regional Government of Castile and León (CSI144P20). J.P.L. and P.L.S. are supported by the Carlos III Health Institute (PIE14/00066). AGN laboratory and human patients’ studies are supported by an ISCIII project grant (PI18/01242). The Human Genotyping unit is a member of CeGen, PRB3, and is supported by grant PT17/0019 of the PE I + D + i 2013–2016, funded by ISCIII and ERDF. SCLl is supported by MINECO/FEDER research grants (RTI2018-094130-B-100). CH was supported by the Department of Defense (DoD) BCRP, No. BC190820; and the National Cancer Institute (NCI) at the National Institutes of Health (NIH), No. R01CA184476. Lawrence Berkeley National Laboratory (LBNL) is a multi-program national laboratory operated by the University of California for the DOE under contract DE AC02-05CH11231. The Proteomics Unit belongs to ProteoRed, PRB3-ISCIII, supported by grant PT17/0019/0023 of the PE I + D +i, 2017–2020, funded by ISCIII and FEDER. RCC is funded by fellowships from the Spanish Regional Government of Castile and León. NGS is a recipient of an FPU fellowship (MINECO/FEDER). hiPSC-CM studies were funded in part by the “la Caixa” Banking Foundation under the project code HR18-00304 and a Severo Ochoa CNIC Intramural Project (Exp. 12-2016 IGP) to J.J.S

Assessment of variation in immunosuppressive pathway genes reveals TGFBR2 to be associated with prognosis of estrogen receptor-negative breast cancer after chemotherapy

Introduction: Tumor lymphocyte infiltration is associated with clinical response to chemotherapy in estrogen receptor (ER) negative breast cancer. To identify variants in immunosuppressive pathway genes associated with prognosis after adjuvant chemotherapy for ER-negative patients, we studied stage I-III invasive breast cancer patients of European ancestry, including 9,334 ER-positive (3,151 treated with chemotherapy) and 2,334 ER-negative patients (1,499 treated with chemotherapy). Methods: We pooled data from sixteen studies from the Breast Cancer Association Consortium (BCAC), and employed two independent studies for replications. Overall 3,610 single nucleotide polymorphisms (SNPs) in 133 genes were genotyped as part of the Collaborative Oncological Gene-environment Study, in which phenotype and clinical data were collected and harmonized. Multivariable Cox proportional hazard regression was used to assess genetic associations with overall survival (OS) and breast

Genetic predisposition to in situ and invasive lobular carcinoma of the breast.

Invasive lobular breast cancer (ILC) accounts for 10-15% of all invasive breast carcinomas. It is generally ER positive (ER+) and often associated with lobular carcinoma in situ (LCIS). Genome-wide association studies have identified more than 70 common polymorphisms that predispose to breast cancer, but these studies included predominantly ductal (IDC) carcinomas. To identify novel common polymorphisms that predispose to ILC and LCIS, we pooled data from 6,023 cases (5,622 ILC, 401 pure LCIS) and 34,271 controls from 36 studies genotyped using the iCOGS chip. Six novel SNPs most strongly associated with ILC/LCIS in the pooled analysis were genotyped in a further 516 lobular cases (482 ILC, 36 LCIS) and 1,467 controls. These analyses identified a lobular-specific SNP at 7q34 (rs11977670, OR (95%CI) for ILC = 1.13 (1.09-1.18), P = 6.0 × 10(-10); P-het for ILC vs IDC ER+ tumors = 1.8 × 10(-4)). Of the 75 known breast cancer polymorphisms that were genotyped, 56 were associated with ILC and 15 with LCIS at P<0.05. Two SNPs showed significantly stronger associations for ILC than LCIS (rs2981579/10q26/FGFR2, P-het = 0.04 and rs889312/5q11/MAP3K1, P-het = 0.03); and two showed stronger associations for LCIS than ILC (rs6678914/1q32/LGR6, P-het = 0.001 and rs1752911/6q14, P-het = 0.04). In addition, seven of the 75 known loci showed significant differences between ER+ tumors with IDC and ILC histology, three of these showing stronger associations for ILC (rs11249433/1p11, rs2981579/10q26/FGFR2 and rs10995190/10q21/ZNF365) and four associated only with IDC (5p12/rs10941679; rs2588809/14q24/RAD51L1, rs6472903/8q21 and rs1550623/2q31/CDCA7). In conclusion, we have identified one novel lobular breast cancer specific predisposition polymorphism at 7q34, and shown for the first time that common breast cancer polymorphisms predispose to LCIS. We have shown that many of the ER+ breast cancer predisposition loci also predispose to ILC, although there is some heterogeneity between ER+ lobular and ER+ IDC tumors. These data provide evidence for overlapping, but distinct etiological pathways within ER+ breast cancer between morphological subtypes

Combined Associations of a Polygenic Risk Score and Classical Risk Factors With Breast Cancer Risk.

We evaluated the joint associations between a new 313-variant PRS (PRS313) and questionnaire-based breast cancer risk factors for women of European ancestry, using 72 284 cases and 80 354 controls from the Breast Cancer Association Consortium. Interactions were evaluated using standard logistic regression and a newly developed case-only method for breast cancer risk overall and by estrogen receptor status. After accounting for multiple testing, we did not find evidence that per-standard deviation PRS313 odds ratio differed across strata defined by individual risk factors. Goodness-of-fit tests did not reject the assumption of a multiplicative model between PRS313 and each risk factor. Variation in projected absolute lifetime risk of breast cancer associated with classical risk factors was greater for women with higher genetic risk (PRS313 and family history) and, on average, 17.5% higher in the highest vs lowest deciles of genetic risk. These findings have implications for risk prevention for women at increased risk of breast cancer

Functional annotation of the 2q35 breast cancer risk locus implicates a structural variant in influencing activity of a long-range enhancer element

A combination of genetic and functional approaches has identified three independent breast cancer risk loci at 2q35. A recent fine-scale mapping analysis to refine these associations resulted in 1 (signal 1), 5 (signal 2), and 42 (signal 3) credible causal variants at these loci. We used publicly available in silico DNase I and ChIP-seq data with in vitro reporter gene and CRISPR assays to annotate signals 2 and 3. We identified putative regulatory elements that enhanced cell-type-specific transcription from the IGFBP5 promoter at both signals (30-to 40-fold increased expression by the putative regulatory element at signal 2, 2- to 3-fold by the putative regulatory element at signal 3). We further identified one of the five credible causal variants at signal 2, a 1.4 kb deletion (esv3594306), as the likely causal variant; the deletion allele of this variant was associated with an average additional increase in IGFBP5 expression of 1.3-fold (MCF-7) and 2.2-fold (T-47D). We propose a model in which the deletion allele of esv3594306 juxtaposes two transcription factor binding regions (annotated by estrogen receptor alpha ChIP-seq peaks) to generate a single extended regulatory element. This regulatory element increases cell-type-specific expression of the tumor suppressor gene IGFBP5 and, thereby, reduces risk of estrogen receptor-positive breast cancer (odds ratio = 0.77, 95% CI 0.74-0.81, p = 3.1 x 10(-31)).Peer reviewe

PLoS One

Age-related macular degeneration (AMD) is a common, progressive multifactorial vision-threatening disease and many genetic and environmental risk factors have been identified. The risk of AMD is influenced by lifestyle and diet, which may be reflected by an altered metabolic profile. Therefore, measurements of metabolites could identify biomarkers for AMD, and could aid in identifying high-risk individuals. Hypothesis-free technologies such as metabolomics have a great potential to uncover biomarkers or pathways that contribute to disease pathophysiology. To date, only a limited number of metabolomic studies have been performed in AMD. Here, we aim to contribute to the discovery of novel biomarkers and metabolic pathways for AMD using a targeted metabolomics approach of 188 metabolites. This study focuses on non-advanced AMD, since there is a need for biomarkers for the early stages of disease before severe visual loss has occurred. Targeted metabolomics was performed in 72 patients with early or intermediate AMD and 72 control individuals, and metabolites predictive for AMD were identified by a sparse partial least squares discriminant analysis. In our cohort, we identified four metabolite variables that were most predictive for early and intermediate stages of AMD. Increased glutamine and phosphatidylcholine diacyl C28:1 levels were detected in non-advanced AMD cases compared to controls, while the rate of glutaminolysis and the glutamine to glutamate ratio were reduced in non-advanced AMD. The association of glutamine with non-advanced AMD corroborates a recent report demonstrating an elevated glutamine level in early AMD using a different metabolomics technique. In conclusion, this study indicates that metabolomics is a suitable method for the discovery of biomarker candidates for AMD. In the future, larger metabolomics studies could add to the discovery of novel biomarkers in yet unknown AMD pathways and expand our insights in AMD pathophysiology

Ophthalmology

OBJECTIVE: In the current study we aimed to identify metabolites associated with age-related macular degeneration (AMD) by performing the largest metabolome association analysis in AMD to date. In addition, we aimed to determine the effect of AMD-associated genetic variants on metabolite levels, and aimed to investigate associations between the identified metabolites and activity of the complement system, one of the main AMD-associated disease pathways. DESIGN: Case-control assocation analysis of metabolomics data. SUBJECTS: 2,267 AMD cases and 4,266 controls from five European cohorts. METHODS: Metabolomics was performed using a high-throughput H-NMR metabolomics platform, which allows the quantification of 146 metabolite measurements and 79 derivative values. Metabolome-AMD associations were studied using univariate logistic regression analyses. The effect of 52 AMD-associated genetic variants on the identified metabolites was investigated using linear regression. In addition, associations between the identified metabolites and activity of the complement pathway (defined by the C3d/C3 ratio) were investigated using linear regression. MAIN OUTCOME MEASURES: Metabolites associated with AMD RESULTS: We identified 60 metabolites that were significantly associated with AMD, including increased levels of large and extra-large HDL subclasses and decreased levels of VLDL, amino acids and citrate. Out of 52 AMD-associated genetic variants, seven variants were significantly associated with 34 of the identified metabolites. The strongest associations were identified for genetic variants located in or near genes involved in lipid metabolism (ABCA1, CETP, APOE, LIPC) with metabolites belonging to the large and extra-large HDL subclasses. In addition, 57 out of 60 metabolites were significantly associated with complement activation levels, and these associations were independent of AMD status. Increased large and extra-large HDL levels and decreased VLDL and amino acid levels were associated with increased complement activation. CONCLUSIONS: Lipoprotein levels were associated with AMD-associated genetic variants, while decreased essential amino acids may point to nutritional deficiencies in AMD. We observed strong associations between the vast majority of the AMD-associated metabolites and systemic complement activation levels, independent of AMD status. This may indicate biological interactions between the main AMD disease pathways, and suggests that multiple pathways may need to be targeted simultaneously for successful treatment of AMD

Assessment of variation in immunosuppressive pathway genes reveals TGFBR2 to be associated with prognosis of estrogen receptor-negative breast cancer after chemotherapy

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.Abstract Introduction Tumor lymphocyte infiltration is associated with clinical response to chemotherapy in estrogen receptor (ER) negative breast cancer. To identify variants in immunosuppressive pathway genes associated with prognosis after adjuvant chemotherapy for ER-negative patients, we studied stage I-III invasive breast cancer patients of European ancestry, including 9,334 ER-positive (3,151 treated with chemotherapy) and 2,334 ER-negative patients (1,499 treated with chemotherapy). Methods We pooled data from sixteen studies from the Breast Cancer Association Consortium (BCAC), and employed two independent studies for replications. Overall 3,610 single nucleotide polymorphisms (SNPs) in 133 genes were genotyped as part of the Collaborative Oncological Gene-environment Study, in which phenotype and clinical data were collected and harmonized. Multivariable Cox proportional hazard regression was used to assess genetic associations with overall survival (OS) and breast cancer-specific survival (BCSS). Heterogeneity according to chemotherapy or ER status was evaluated with the log-likelihood ratio test. Results Three independent SNPs in TGFBR2 and IL12B were associated with OS (P  C) (per allele hazard ratio (HR) 1.54 (95% confidence interval (CI) 1.22 to 1.95), P = 3.08 × 10−4) was not found in ER-negative patients without chemotherapy or ER-positive patients with chemotherapy (P for interaction  A) with poorer OS (HR 1.50 (95% CI 1.21 to 1.86), P = 1.81 × 10−4), and rs2853694 (A > C) with improved OS (HR 0.73 (95% CI 0.61 to 0.87), P = 3.67 × 10−4). Similar associations were observed with BCSS. Association with TGFBR2 rs1367610 but not IL12B variants replicated using BCAC Asian samples and the independent Prospective Study of Outcomes in Sporadic versus Hereditary Breast Cancer Study and yielded a combined HR of 1.57 ((95% CI 1.28 to 1.94), P = 2.05 × 10−5) without study heterogeneity. Conclusions TGFBR2 variants may have prognostic and predictive value in ER-negative breast cancer patients treated with adjuvant chemotherapy. Our findings provide further insights into the development of immunotherapeutic targets for ER-negative breast cancer