Artificial Intelligence-driven Prediction Revealed CFTR Associated With Therapy Outcome Of Breast Cancer: A Feasibility Study.

INTRODUCTION: In silico tools capable of predicting the functional consequences of genomic differences between individuals, many of which are AI-driven, have been the most effective over the past two decades for non-synonymous single nucleotide variants (nsSNVs). When appropriately selected for the purpose of the study, a high predictive performance can be expected. In this feasibility study, we investigate the distribution of nsSNVs with an allele frequency below 5%. To classify the putative functional consequence, a tier-based filtration led by AI-driven predictors and scoring system were implemented to the overall-decision making process, resulting in a list of prioritised genes. METHODS: The study has been conducted on breast cancer patients of homogeneous ethnicity. Germline rare variants have been sequenced in genes that influence pharmacokinetic parameters of anticancer drugs or molecular signalling pathways in cancer. After AI-driven functional pathogenicity classification and data mining in pharmacogenomic (PGx) databases, variants were collapsed to the gene level and ranked according to their putative deleterious role. RESULTS: In breast cancer patients, seven of the twelve genes prioritised based on the predictions were found to be associated with response to oncotherapy, histological grade, and tumour subtype. Most importantly, we showed that the group of patients with at least one rare nsSNVs in Cystic Fibrosis Transmembrane Conductance Regulator (CFTR) had significantly reduced disease-free (Log Rank, p=0.002) and overall survival (Log Rank, p=0.006). CONCLUSION: AI-driven in silico analysis with PGx data mining provided an effective approach navigating for functional consequences across germline genetic background, which can be easily integrated into the overall decision-making process for future studies. The study revealed a statistically significant association with numerous clinicopathological parameters, including treatment response. Our study indicates that CFTR may be involved in the processes influencing the effectiveness of oncotherapy or in the malignant progression of the disease itself.

Integrative analysis of mRNA and miRNA expression profiles and somatic variants in oxysterol signaling in early-stage luminal breast cancer.

Oxysterols, oxidized derivatives of cholesterol, act in breast cancer (BC) as selective estrogen receptor modulators and affect cholesterol homeostasis, drug transport, nuclear and cell receptors, and other signaling proteins. Using data from three highly overlapping sets of patients (N = 162 in total) with early-stage estrogen-receptor-positive luminal BC-high-coverage targeted DNA sequencing (113 genes), mRNA sequencing, and full micro-RNA (miRNA) transcriptome microarrays-we describe complex oxysterol-related interaction (correlation) networks, with validation in public datasets (n = 538) and 11 databases. The ESR1-CH25H-INSIG1-ABCA9 axis was the most prominent, interconnected through miR-125b-5p, miR-99a-5p, miR-100-5p, miR-143-3p, miR-199b-5p, miR-376a-3p, and miR-376c-3p. Mutations in SC5D, CYP46A1, and its functionally linked gene set were associated with multiple differentially expressed oxysterol-related genes. STARD5 was upregulated in patients with positive lymph node status. High expression of hsa-miR-19b-3p was weakly associated with poor survival. This is the first study of oxysterol-related genes in BC that combines DNA, mRNA, and miRNA multiomics with detailed clinical data. Future studies should provide links between intratumoral oxysterol signaling depicted here, circulating oxysterol levels, and therapy outcomes, enabling eventual clinical exploitation of present findings.

Single Nucleotide Variants in KIF14 Gene May Have Prognostic Value in Breast Cancer.

INTRODUCTION: Human kinesin 14 (KIF14) is one of the 70 prognostic marker genes (so-called Amsterdam profile) previously identified by the microarray of breast carcinomas, and its high transcript expression in tumor specimens indicates a poor prognosis for patients. We performed a pilot study to explore the prognostic and predictive meaning of KIF14 germline genetic variability in breast cancer patients. METHODS: KIF14 coding sequence, including 5' and 3' untranslated regions and overlaps to introns for identification of splicing sites, was analyzed using next-generation sequencing in the testing set of blood DNA samples from 105 breast cancer patients with clinical follow-up. After rigorous evaluation of major allele frequency, haplotype blocks, in silico predicted functional aspects, expression quantitative trait loci, and clinical associations, eight single nucleotide variants were subsequently validated in the evaluation set of 808 patients. RESULTS: Carriers of minor alleles G (rs17448931) or T (rs3806362) had significantly shorter overall survival than wild type homozygotes (p = 0.010 and p = 0.023, respectively) thus successfully replicating the results of the testing set. Both associations remained significant in the multivariate Cox regression analysis, including molecular subtype and stage as covariates (hazard ratio, HR = 1.7, 95% confidence interval (CI) = 1.1-2.8 for rs17448931 and HR = 1.9, CI 1.2-3.0 for rs3806362). DISCUSSION: In conclusion, our preliminary data suggest that minor alleles in rs17448931 and rs3806362 of KIF14 represent candidate biomarkers of poor prognosis of breast cancer patients. After pending validation in independent populations and eventual functional characterization, these candidates might become useful biomarkers in the clinics.

Germline and somatic genetic variability of oxysterol-related genes in breast cancer patients with early disease of the luminal subtype.

Oxysterols, oxidized derivatives of cholesterol, have been implicated in multiple pathologies, including cancer. In breast cancer, the link is especially strong due to interactions between oxysterols and estrogen receptor activity. Here, we provide the first dedicated study of 113 oxysterol-related genes in breast cancer patients of the luminal subtype, in terms of both their somatic and germline variability, using targeted high-throughput DNA sequencing of 100 normal-tumor pairs with very high coverage. In the full cohort, or subsets of patients stratified by therapy, we found 12 germline variants in ABCA1, ABCA8, ABCC1, GPR183, LDLR, MBTPS1, NR1I2, OSBPL2, OSBPL3, and OSBPL5 to associate with poor survival of patients and variants in ABCA8, ABCG2, and HSD3B7 (three in total) associated with better survival. However, no associations remained significant after correction for multiple tests. Analysis of somatic variants revealed significantly (after FDR correction) poorer survival in patients mutated in CYP46A1 and 9 interacting (according to STRING analysis) genes, as well as in OSBPL3 and a set of 20 genes that collectively associated with the progesterone receptor status of patients. We propose further exploration of these genes in an integrative manner together with gene expression and epigenomic data.

Genetic variations in 3'UTRs of SMUG1 and NEIL2 genes modulate breast cancer risk, survival and therapy response.

Breast cancer (BC) is the most frequent malignancy in women accounting for approximately 2 million new cases worldwide annually. Several genetic, epigenetic and environmental factors are known to be involved in BC development and progression, including alterations in post-transcriptional gene regulation mediated by microRNAs (miRNAs). Single nucleotide polymorphisms (SNPs) located in miRNA binding sites (miRSNPs) in 3'-untranslated regions of target genes may affect miRNA-binding affinity and consequently modulate gene expression. We have previously reported a significant association of miRSNPs in the SMUG1 and NEIL2 genes with overall survival in colorectal cancer patients. SMUG1 and NEIL2 are DNA glycosylases involved in base excision DNA repair. Assuming that certain genetic traits are common for solid tumours, we have investigated wherever variations in SMUG1 and NEIL2 genes display an association with BC risk, prognosis, and therapy response in a group of 673 BC patients and 675 healthy female controls. Patients with TC genotype of NEIL2 rs6997097 and receiving only hormonal therapy displayed markedly shorter overall survival (HR = 4.15, 95% CI = 1.7-10.16, P = 0.002) and disease-free survival (HR = 2.56, 95% CI = 1.5-5.7, P = 0.02). Our results suggest that regulation of base excision repair glycosylases operated by miRNAs may modulate the prognosis of hormonally treated BC.

Role of Genetic Variation in Cytochromes P450 in Breast Cancer Prognosis and Therapy Response.

Breast cancer is the most frequent cancer in the female population worldwide. The role of germline genetic variability in cytochromes P450 (CYP) in breast cancer prognosis and individualized therapy awaits detailed elucidation. In the present study, we used the next-generation sequencing to assess associations of germline variants in the coding and regulatory sequences of all human CYP genes with response of the patients to the neoadjuvant cytotoxic chemotherapy and disease-free survival (n = 105). A total of 22 prioritized variants associating with a response or survival in the above evaluation phase were then analyzed by allelic discrimination in the large confirmation set (n = 802). Associations of variants in CYP1B1, CYP4F12, CYP4X1, and TBXAS1 with the response to the neoadjuvant cytotoxic chemotherapy were replicated by the confirmation phase. However, just association of variant rs17102977 in CYP4X1 passed the correction for multiple testing and can be considered clinically and statistically validated. Replicated associations for variants in CYP4X1, CYP24A1, and CYP26B1 with disease-free survival of all patients or patients stratified to subgroups according to therapy type have not passed a false discovery rate test. Although statistically not confirmed by the present study, the role of CYP genes in breast cancer prognosis should not be ruled out. In conclusion, the present study brings replicated association of variant rs17102977 in CYP4X1 with the response of patients to the neoadjuvant cytotoxic chemotherapy and warrants further research of genetic variation CYPs in breast cancer.

SLC46A1 Haplotype with Predicted Functional Impact has Prognostic Value in Breast Carcinoma.

BACKGROUND AND OBJECTIVE: Membrane solute carrier transporters play an important role in the transport of a wide spectrum of substrates including anticancer drugs and cancer-related physiological substrates. This study aimed to assess the prognostic relevance of gene expression and genetic variability of selected solute carrier transporters in breast cancer. METHODS: Gene expression was determined by quantitative real-time polymerase chain reaction. All SLC46A1 and SLCO1A2 exons and surrounding non-coding sequences in DNA extracted from the blood of patients with breast cancer (exploratory phase) were analyzed by next-generation sequencing technology. Common variants (minor allele frequency ≥ 5%) with in silico-predicted functional relevance were further analyzed in a large cohort of patients with breast cancer (n = 815) and their prognostic and predictive potential was estimated (validation phase). RESULTS: A gene expression and bioinformatics analysis suggested SLC46A1 and SLCO1A2 to play a putative role in the prognosis of patients with breast cancer. In total, 135 genetic variants (20 novel) were identified in both genes in the exploratory phase. Of these variants, 130 were non-coding, three missense, and two synonymous. One common variant in SLCO1A2 and four variants in SLC46A1 were predicted to be pathogenic by in silico programs and subsequently validated. A SLC46A1 haplotype block composed of rs2239911-rs2239910-rs8079943 was significantly associated with ERBB2/HER2 status and disease-free survival of hormonally treated patients. CONCLUSIONS: This study revealed the prognostic value of a SLC46A1 haplotype block for breast cancer that should be further studied.

Role of Genetic Variation in ABC Transporters in Breast Cancer Prognosis and Therapy Response.

Breast cancer is the most common cancer in women in the world. The role of germline genetic variability in ATP-binding cassette (ABC) transporters in cancer chemoresistance and prognosis still needs to be elucidated. We used next-generation sequencing to assess associations of germline variants in coding and regulatory sequences of all human ABC genes with response of the patients to the neoadjuvant cytotoxic chemotherapy and disease-free survival (n = 105). A total of 43 prioritized variants associating with response or survival in the above testing phase were then analyzed by allelic discrimination in the large validation set (n = 802). Variants in ABCA4, ABCA9, ABCA12, ABCB5, ABCC5, ABCC8, ABCC11, and ABCD4 associated with response and variants in ABCA7, ABCA13, ABCC4, and ABCG8 with survival of the patients. No association passed a false discovery rate test, however, the rs17822931 (Gly180Arg) in ABCC11, associating with response, and the synonymous rs17548783 in ABCA13 (survival) have a strong support in the literature and are, thus, interesting for further research. Although replicated associations have not reached robust statistical significance, the role of ABC transporters in breast cancer should not be ruled out. Future research and careful validation of findings will be essential for assessment of genetic variation which was not in the focus of this study, e.g., non-coding sequences, copy numbers, and structural variations together with somatic mutations.

Chromosomal damage and telomere length in peripheral blood lymphocytes of cancer patients.

Accumulation of non­specific structural chromosomal aberrations (CAs) and telomere shortening contribute to genome instability, which constitutes as one of the hallmarks of cancer. CAs arise due to direct DNA damage or telomere shortening. CAs in peripheral blood lymphocytes (PBL), which are considered to be markers of exposure, have been previously reported to serve a role in the pathophysiology and progression of cancer through mechanisms that are poorly understood. In addition, the prognostic relevance of telomere length (TL) in patients with cancer remains to be elucidated. In the present study, CAs and TL in PBL isolated from patients with newly diagnosed cancer (151 breast, 96 colorectal, 90 lung) and 335 cancer­free control individuals were investigated. These results were then correlated with clinicopathological factors and follow­up data. The accumulation of CAs in PBL was observed with increased susceptibility to breast and lung cancer (P<0.0001), while individuals with longer TL were found to be at a higher risk of breast cancer (P<0.0001). Increased chromatid­type aberrations were also revealed to be associated with lower overall survival of patients with breast and colorectal cancers using a multivariate model. Compared with control individuals, no association was observed between TL and CAs or age in patients with cancer. In conclusion, the present study demonstrates the association between CAs/TL in PBL and the susceptibility, prognosis and survival of patients with breast, colorectal and lung cancer.

Estrogen Receptor Status Oppositely Modifies Breast Cancer Prognosis in BRCA1/BRCA2 Mutation Carriers Versus Non-Carriers.

Breast cancer (BC) prognosis in BRCA1 and BRCA2 mutation carriers has been reported contradictorily, and the significance of variables influencing prognosis in sporadic BC is not established in BC patients with hereditary BRCA1/BRCA2 mutations. In this retrospective cohort study, we analyzed the effect of clinicopathological characteristics on BC prognosis (disease-free survival [DFS] and disease-specific survival [DSS]) in hereditary BRCA1/BRCA2 mutation carriers. We enrolled 234 BRCA1/BRCA2 mutation carriers and 899 non-carriers, of whom 191 carriers and 680 non-carriers, with complete data, were available for survival analyses. We found that patients with ER-positive tumors developed disease recurrence 2.3-times more likely when they carried a BRCA1/BRCA2 mutation (23/60; 38.3% ER-positive carriers vs. 74/445; 16.6% ER-positive non-carriers; p < 0.001). ER-positive mutation carriers also had a 3.4-times higher risk of death due to BC compared with ER-positive non-carriers (13/60; 21.7% vs. 28/445; 6.3%; p < 0.001). Moreover, prognosis in ER-negative BRCA1/BRCA2 mutation carriers was comparable with that in ER-positive non-carriers. Our study demonstrates that ER-positivity worsens BC prognosis in BRCA1/BRCA2 mutation carriers, while prognosis for carriers with ER-negative tumors (including early-onset) is significantly better and comparable with that in ER-positive, older BC non-carriers. These observations indicate that BRCA1/BRCA2 mutation carriers with ER-positive BC represent high-risk patients.

Use of Germline Genetic Variability for Prediction of Chemoresistance and Prognosis of Breast Cancer Patients.

The aim of our study was to set up a panel for targeted sequencing of chemoresistance genes and the main transcription factors driving their expression and to evaluate their predictive and prognostic value in breast cancer patients. Coding and regulatory regions of 509 genes, selected from PharmGKB and Phenopedia, were sequenced using massive parallel sequencing in blood DNA from 105 breast cancer patients in the testing phase. In total, 18,245 variants were identified of which 2565 were novel variants (without rs number in dbSNP build 150) in the testing phase. Variants with major allele frequency over 0.05 were further prioritized for validation phase based on a newly developed decision tree. Using emerging in silico tools and pharmacogenomic databases for functional predictions and associations with response to cytotoxic therapy or disease-free survival of patients, 55 putative variants were identified and used for validation in 805 patients with clinical follow up using KASPTM technology. In conclusion, associations of rs2227291, rs2293194, and rs4376673 (located in ATP7A, KCNAB1, and DFFB genes, respectively) with response to neoadjuvant cytotoxic therapy and rs1801160 in DPYD with disease-free survival of patients treated with cytotoxic drugs were validated and should be further functionally characterized.

Genotype and Haplotype Analyses of TP53 Gene in Breast Cancer Patients: Association with Risk and Clinical Outcomes.

Variations in the TP53 gene have been suggested to play a role in many cancers, including breast. We previously observed an association between TP53 haplotypes based on four polymorphisms (rs17878362, rs1042522, rs12947788, and rs17884306) and the risk of colorectal and pancreatic cancer. Based on these results, in the present study, we have investigated the same polymorphisms and their haplotypes in 705 breast cancer cases and 611 healthy controls in relation to the disease risk, histopathological features of the tumor and clinical outcomes. In comparison to the most common haplotype A1-G-C-G, all the other identified haplotypes were globally associated with a significantly decreased breast cancer risk (P = 0.006). In particular, the A2-G-C-G haplotype was associated with a marked decreased risk of breast cancer when compared with the common haplotype (P = 0.0001). Moreover, rs1042522 in patients carrying the GC genotype and receiving only the anthracycline-based chemotherapy was associated with both overall and disease-free survival (recessive model for overall survival HR = 0.30 95% CI 0.11-0.80, P = 0.02 and for disease-free survival HR = 0.42 95% CI 0.21-0.84, P = 0.01). Present results suggest common genetic features in the susceptibility to breast and gastrointestinal cancers in respect to TP53 variations. In fact, similar haplotype distributions were observed for breast, colorectal, and pancreatic patients in associations with cancer risk. Rs1042522 polymorphism (even after applying the Dunn-Bonferroni correction for multiple testing) appears to be an independent prognostic marker in breast cancer patients.

The expression profile of ATP-binding cassette transporter genes in breast carcinoma.

AIM: ATP-binding cassette (ABC) transporters contribute to development of resistance to anticancer drugs via ATP-dependent drug efflux. A major goal of our study was to investigate associations between the expression of ABC transporters and outcome of breast carcinoma patients. PATIENTS & METHODS: Transcript levels of all 49 human ABC transporters were determined in post-treatment tumor and non-neoplastic tissue samples from 68 breast carcinoma patients treated by neoadjuvant chemotherapy. Six ABC transporters were then evaluated in independent series of 100 pretreatment patients. RESULTS: ABCA5/6/8/9/10, ABCB1/5/11, ABCC6/9, ABCD2/4, ABCG5 and ABCG8 were significantly downregulated and ABCA2/3/7/12, ABCB2/3/8/9/10, ABCC1/4/5/10/11/12, ABCD1/3, ABCE1, ABCF1/2/3 and ABCG1 were upregulated in post-treatment tumors compared with non-neoplastic tissues. Significant associations of intratumoral levels of ABCC1 and ABCC8 with grade and expression of hormonal receptors were found in both sets of patients. ABCA12, ABCA13 and ABCD2 levels were significantly associated with the response to neoadjuvant chemotherapy in post-treatment patients. Protein expression of ABCA12, ABCC8 and ABCD2 in tumor tissues of patients with breast carcinoma was observed by immunoblotting for the first time. CONCLUSION: ABCA12, ABCA13, ABCC1, ABCC8 and ABCD2 present potential modifiers of progression and response to the chemotherapy of breast carcinoma.

Detection of frequent ABCB1 polymorphisms by high-resolution melting curve analysis and their effect on breast carcinoma prognosis.

BACKGROUND: The ABCB1 gene encodes P-glycoprotein implicated in the development of cellular drug resistance. The aim of this study was to develop high-resolution melting (HRM) analysis for determination of ABCB1 polymorphisms and evaluate their associations with clinical data of breast carcinoma patients. METHODS: HRM analysis was designed to assess five single nucleotide polymorphisms (SNPs) in ABCB1 (rs2214102, rs1128503, rs2032582, rs2032583 and rs1045642) in genomic DNA from 103 breast carcinoma patients. Results were confirmed by direct DNA sequencing. RESULTS: HRM analysis revealed distinct patterns of melting curves for the respective genotypes of all followed SNPs. Sensitivity of HRM analysis compared with direct DNA sequencing was superior (97.1% vs. 93.9%). The overall accuracy of HRM was 97.6%. The coefficients of variation in replicate experiments encompassed the range 0.002%-0.038%. On the basis of the examined SNPs, one strong haplotype block containing rs2032582 and rs1128503 SNPs was identified. Significant associations of rs2032582 SNP with tumor size, negative HER-2/neu status, and family history of breast carcinoma were found. Patients carrying the ancestral homozygous genotype (GG) in rs2214102 had significantly worse progression-free survival in comparison with carriers of the non-ancestral allele (A) in the adjuvant set (p=0.005). CONCLUSIONS: A rapid, accurate, low-cost and time-effective method for screening ABCB1 SNPs was developed. Significant associations of ABCB1 rs2032582 and rs2214102 SNPs with prognostic factors and survival of patients were found.