High miR-30 Expression Associates with Improved Breast Cancer Patient Survival and Treatment Outcome.

Deregulated miRNA expression has been suggested in several stages of breast cancer pathogenesis. We have studied the miR-30 family, in particular miR-30d, in relation to breast cancer patient survival and treatment outcomes. With tumor specimens from 1238 breast cancer patients, we analyzed the association of miR-30d expression with tumor characteristics with the 5-year occurrence of breast cancer-specific death or distant metastasis (BDDM), and with 10-year breast cancer survival (BCS). We conducted a two-stage drug-screen to investigate the impact of miR-30 family members (miR-30a-30e) on sensitivity to doxorubicin and lapatinib in six breast cancer cell lines HCC1937, HCC1954, MDA-MB-361, MCF7, MDA-MB-436 and CAL-120, using drug sensitivity scores (DSS) to compare the miR-30 family mimics to their specific inhibitors. The study was complemented with Ingenuity Pathway Analysis (IPA) with the METABRIC data. We found that while high miR-30d expression is typical for aggressive tumors, it predicts better metastasis-free (pBDDM = 0.035, HR = 0.63, 95% CI = 0.4-0.9) and breast cancer-specific survival (pBCS = 0.018, HR = 0.61, 95% CI = 0.4-0.9), especially in HER2-positive (pBDDM = 0.0009), ER-negative (pBDDM = 0.003), p53-positive (pBDDM = 0.011), and highly proliferating (pBDDM = 0.0004) subgroups, and after adjuvant chemotherapy (pBDDM = 0.035). MiR-30d predicted survival independently of standard prognostic markers (pBDDM = 0.0004). In the drug-screening test, the miR-30 family sensitized the HER2-positive HCC1954 cell line to lapatinib (p < 10-2) and HCC1937, MDA-MB-361, MDA-MB-436 and CAL120 to doxorubicin (p < 10-4) with an opposite impact on MCF7. According to the pathway analysis, the miR-30 family has a suppressive effect on cell motility and metastasis in breast cancer. Our results suggest prognostic and predictive potential for the miR-30 family, which warrants further investigation.

TP53-based interaction analysis identifies cis-eQTL variants for TP53BP2, FBXO28, and FAM53A that associate with survival and treatment outcome in breast cancer.

TP53 overexpression is indicative of somatic TP53 mutations and associates with aggressive tumors and poor prognosis in breast cancer. We utilized a two-stage SNP association study to detect variants associated with breast cancer survival in a TP53-dependent manner. Initially, a genome-wide study (n = 575 cases) was conducted to discover candidate SNPs for genotyping and validation in the Breast Cancer Association Consortium (BCAC). The SNPs were then tested for interaction with tumor TP53 status (n = 4,610) and anthracycline treatment (n = 17,828). For SNPs interacting with anthracycline treatment, siRNA knockdown experiments were carried out to validate candidate genes.In the test for interaction between SNP genotype and TP53 status, we identified one locus, represented by rs10916264 (p(interaction) = 3.44 × 10-5; FDR-adjusted p = 0.0011) in estrogen receptor (ER) positive cases. The rs10916264 AA genotype associated with worse survival among cases with ER-positive, TP53-positive tumors (hazard ratio [HR] 2.36, 95% confidence interval [C.I] 1.45 - 3.82). This is a cis-eQTL locus for FBXO28 and TP53BP2; expression levels of these genes were associated with patient survival specifically in ER-positive, TP53-mutated tumors. Additionally, the SNP rs798755 was associated with survival in interaction with anthracycline treatment (p(interaction) = 9.57 × 10-5, FDR-adjusted p = 0.0130). RNAi-based depletion of a predicted regulatory target gene, FAM53A, indicated that this gene can modulate doxorubicin sensitivity in breast cancer cell lines.If confirmed in independent data sets, these results may be of clinical relevance in the development of prognostic and predictive marker panels for breast cancer.

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

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

SNP-SNP interaction analysis of NF-κB signaling pathway on breast cancer survival.

In breast cancer, constitutive activation of NF-κB has been reported, however, the impact of genetic variation of the pathway on patient prognosis has been little studied. Furthermore, a combination of genetic variants, rather than single polymorphisms, may affect disease prognosis. Here, in an extensive dataset (n = 30,431) from the Breast Cancer Association Consortium, we investigated the association of 917 SNPs in 75 genes in the NF-κB pathway with breast cancer prognosis. We explored SNP-SNP interactions on survival using the likelihood-ratio test comparing multivariate Cox' regression models of SNP pairs without and with an interaction term. We found two interacting pairs associating with prognosis: patients simultaneously homozygous for the rare alleles of rs5996080 and rs7973914 had worse survival (HRinteraction 6.98, 95% CI=3.3-14.4, P=1.42E-07), and patients carrying at least one rare allele for rs17243893 and rs57890595 had better survival (HRinteraction 0.51, 95% CI=0.3-0.6, P = 2.19E-05). Based on in silico functional analyses and literature, we speculate that the rs5996080 and rs7973914 loci may affect the BAFFR and TNFR1/TNFR3 receptors and breast cancer survival, possibly by disturbing both the canonical and non-canonical NF-κB pathways or their dynamics, whereas, rs17243893-rs57890595 interaction on survival may be mediated through TRAF2-TRAIL-R4 interplay. These results warrant further validation and functional analyses.

Germline variation in TP53 regulatory network genes associates with breast cancer survival and treatment outcome.

Germline variation in the TP53 network genes PRKAG2, PPP2R2B, CCNG1, PIAS1 and YWHAQ was previously suggested to have an impact on drug response in vitro. Here, we investigated the effect on breast cancer survival of germline variation in these genes in 925 Finnish breast cancer patients and further analyzed five single nucleotide polymorphisms (SNPs) in PRKAG2 (rs1029946, rs4726050, rs6464153, rs7789699) and PPP2R2B (rs10477313) for 10-year survival in breast cancer patients, interaction with TP53 R72P and MDM2-SNP309, outcome after specific adjuvant therapy and correlation to tumor characteristics in 4,701 invasive cases from four data sets. We found evidence for carriers of PRKAG2-rs1029946 and PRKAG2-rs4726050 having improved survival in the pooled data (HR 0.53, 95% CI 0.3-0.9; p = 0.023 for homozygous carriers of the rare G-allele and HR 0.85, 95% CI 0.7-0.9; p = 0.049 for carriers of the rare G allele, respectively). PRKAG2-rs4726050 showed a significant interaction with MDM2-SNP309, with PRKAG2-rs4726050 rare G-allele having a dose-dependent effect for better breast cancer survival confined only to MDM2 SNP309 rare G-allele carriers (HR 0.45, 95% CI 0.2-0.7; p = 0.001). This interaction also emerged as an independent predictor of better survival (p = 0.047). PPP2R2B-rs10477313 rare A-allele was found to predict better survival (HR 0.82, 95% CI 0.6-0.9; p = 0.018), especially after hormonal therapy (HR 0.66, 95% CI 0.5-0.9; p = 0.048). These findings warrant further studies and suggest that genetic markers in TP53 network genes such as PRKAG2 and PPP2R2B might affect prognosis and treatment outcome in breast cancer patients.

NQO1 expression correlates inversely with NFκB activation in human breast cancer.

NQO1 participates in cellular defense against oxidative stress and regulates apoptosis via p53- and NFκB-mediated pathways. We have previously found that homozygous missense variant NQO1*2 (rs1800566) predicts poor survival among breast cancer patients, particularly after anthracycline-based adjuvant chemotherapy. Here, we investigated NQO1 and NFκB protein expression and global gene expression profiles in breast tumors with correlation to tumor characteristics and survival after adjuvant chemotherapy. We used immunohistochemical analysis of tissue microarrays to study NQO1 and NFκB expression in two series of tumors: 1000 breast tumors unselected for treatment and 113 from a clinical trial comparing chemotherapy regimens after anthracycline treatment in advanced breast cancer. We used gene expression arrays to define genes co-expressed with NQO1 and NFκB. NQO1 and nuclear NFκB were expressed in 83% and 11% of breast tumors, and correlated inversely (P = 0.012). NQO1 protein expression was associated with estrogen receptor (ER) expression (P = 0.011), whereas 34.5% of NFκB-nuclear/activated tumors were ER negative (P = 0.001). NQO1 protein expression and NFκB activation showed only trends, but no statistical significance for patient survival or outcome after anthracycline treatment. Gene expression analysis highlighted 193 genes that significantly correlated with both NQO1 and NFκB in opposite directions, consistent with the expression patterns of the two proteins. Inverse correlation was found with genes related to oxidation/reduction, lipid biosynthesis and steroid metabolism, immune response, lymphocyte activation, Jak-STAT signaling and apoptosis. The inverse relationship between NQO1 protein expression and NFκB activation, underlined also by inverse patterns of association with ER and gene expression profiles of tumors, suggests that NQO1-NFκB interaction in breast cancer is different from several other tissue types, possibly due to estrogen receptor signaling in breast cancer. Neither NQO1 nor NFκB protein expression appear as significant prognostic or predictive markers in breast cancer.

MiR-34a expression has an effect for lower risk of metastasis and associates with expression patterns predicting clinical outcome in breast cancer.

MiR-34a acts as a candidate tumour suppressor gene, and its expression is reduced in several cancer types. We aimed to study miR-34a expression in breast cancer and its correlation with tumour characteristics and clinical outcome, and regulatory links with other genes. We analysed miR-34a expression in 1,172 breast tumours on TMAs. 25% of the tumours showed high, 43% medium and 32% low expression of miR-34a. High miR-34a expression associated with poor prognostic factors for breast cancer: positive nodal status (p = 0.006), high tumour grade (p<0.0001), ER-negativity (p = 0.0002), HER2-positivity (p = 0.0002), high proliferation rate (p<0.0001), p53-positivity (p<0.0001), high cyclin E (p<0.0001) and γH2AX (p<0.0001). However, multivariate analysis adjusting for conventional prognostic factors indicated that high miR-34a expression in fact associated with a lower risk of recurrence or death from breast cancer (HR = 0.63, 95% CI = 0.41-0.96, p = 0.031). Gene expression analysis by differential miR-34a expression revealed an expression signature with an effect on both the 5-year and 10-year survival of the patients (p<0.001). Functional genomic analysis highlighted a novel regulatory role of the transcription factor MAZ, apart from the known control by p53, on the expression of miR-34a and a number of miR-34a targets. Our findings suggest that while miR-34a expression activation is a marker of aggressive breast tumour phenotype it exerts an independent effect for a lower risk of recurrence or death from breast cancer. We also present an expression signature of 190 genes associated with miR-34a expression. Our analysis for regulatory loops suggest that MAZ and p53 transcription factors co-operate in modulating miR-34a, as well as miR-34a targets involved in several cellular pathways. Taken together, these results suggest that the network of genes co-regulated with and targeted by miR-34a form a group of down-stream effectors that maybe of use in predicting clinical outcome, and that highlight novel regulatory mechanisms in breast cancer.