Genetic and molecular biology of breast cancer among Iranian patients

Abstract Background, Breast cancer (BC) is one of the leading causes of cancer related deaths in Iran. This high ratio of mortality had a rising trend during the recent years which is probably associated with late diagnosis. Main body Therefore it is critical to define a unique panel of genetic markers for the early detection among our population. In present review we summarized all of the reported significant genetic markers among Iranian BC patients for the first time, which are categorized based on their cellular functions. Conclusions This review paves the way of introducing a unique ethnic specific panel of diagnostic markers among Iranian BC patients. Indeed, this review can also clarify the genetic and molecular bases of BC progression among Iranians.


Background
Cancer is the second and third leading cause of deaths in developed and developing countries, respectively [1]. The China, India, and Japan with 59.7% of Asian breast cancer (BC) cases have the highest incidences in Asia [2]. BC is the most common cancer among Iranian females with mean age of between 46 and 49 years old [3]. Unfortunately, it has been shown that the BC mortality rate had a noticeable raising trend in 1990 to 2010 from 0.96 to 4.33 per 100,000 cases [4]. There is a variation about the BC incidence based on geographical districts in Iran in which the central provinces such as Tehran and Isfahan has the highest incidences [5]. Although, Iran is among the low incidence countries in the case of BC, the Iranian patients have lower ages and advanced stage tumors in comparison with another patients in western countries [6]. There are various BC risk factors which can be categorized as genetic and non-genetic. The non-genetic factors are age, life style, early menarche, late menopause, family history, weight, smoking, diet, socioeconomic condition, and air pollution. Age of diagnosis is about 40-50 years in Asian and African patients whereas it is about 60-70 years in Western countries [7,8]. Ethnic is also important in age of diagnosis in which the black women have lower ages in comparison with the white women (58 vs. 62 years) [9]. Early Menarche/Late Menopause also increases the risk of breast cancer. Delayed pregnancy and decreased period of breastfeeding are also associated with elevated BC risk. Smoking before the menopause increases the risk of BC. There is also a direct correlation between alcohol consumption and BC risk. A regular physical activity and fruits and vegetables consumptions also decreases the BC risk. In the case of genetic risk factors, the BRCA1 and BRCA 2 mutations are among the main factors which are reported and involved in BC [10,11]. Beside the BRCA family, the variety of genes are also reported during BC progression such as p53 and PTEN tumor suppressors [12,13]. Members of DNA repair system and cell cycle regulators such as ATM, CHEK2, NBN, PALB2, and RAD50 are also correlated with risk of BC [14]. Late diagnosis in advanced stages of tumor is one of the main reason of high mortality among the BC patients. Therefore, finding new diagnostic markers will improve the early detection and better management of BC among patients. Regarding the differences between populations in the case of genetic factors, it is required to introduce a population specific panel of genetic markers for every population. Despite various studies about the role of different genes in BC progression among Iranians, there is not still any unique panel of markers for the early detection among Iranian patients. Therefore, for the first time in the present study we summarized the significant genetic markers which are reported until now among the Iranian BC patients (Table 1). We categorized these reported markers in different cellular processes based on their cellular functions (Fig. 1).

Non coding RNAs and epigenetic modifications
MiRNAs as short non-coding RNAs can regulate gene expression post-transcriptionally through 3-UTR in the target mRNAs which are involved in various cellular functions including apoptosis and cell proliferation [15]. MiRNAs polymorphisms can change the function of miRNAs through modulation of their expression, processing, or binding sequences [16]. It has been observed that there was a significant correlation between rs11134527 variant in miR-218-2 and rs384262 variant in miR-301b and higher risk of BC among a sub population of Iranian patients [17]. SOX2 self renewal factor is one of the targets of miR-140 and miR-140/SOX2 interaction can be associated with cell survival in BC [18]. Another group also evaluated the serum levels of miR-140-3p expression in Iranian BC patients compared with healthy cases. They showed significant miR-140-3p over expression in BC in comparison with the normal cases. Moreover, there was a correlation between mir-140-3p expression and age, in which the premenopausal and ≤ 48 years old female with BC have shown mir-140-3p over expression. They introduced this marker as a diagnostic factor in premenopausal and in ≤ 48 years old females with BC [19]. BRCA1/BRCA2 are the targets of MiR-146a. It has been shown that the miR-146a is involved in BC. Moreover, there was a significant correlation between rs16917496 of miR-502 and higher risk of BC among a sub population of Iranian patients [20]. It has been shown that there was a correlation between hsa-miR-499 rs3746444 CC and C genotypes and high risk of BC among a subpopulation of Iranian subjects [21]. Mir-608 targets several factors such as INSR, IL-1A, GHR, and TP53 which are involved in tumor growth and apoptosis [22]. It has been observed that there was a significant correlation between rs4919510 polymorphism and BC, in which such polymorphism significantly decreased the BC susceptibility among a sub population of Iranian BC patients [23]. MiR-21 is involved in regulation of cell proliferation and apoptosis through targeting various factors such as PTEN and PDCD4 [24,25]. Another study also showed that there was a high levels of miR-21 expression in aggressive breast tumors. Moreover, a significant under expression of miR-205 was observed in ER/PR/HER2 negative breast tumors compared with normal samples. The miR-342 under and over expressions were also observed in ER/PR negative and ER/PR positive BC tissues respectively in an Iranian population [26]. Although, MiR-34a and let-7 are as tumor suppressors are observed in normal cells, their under expression is also observed in tumors. MiR-34a is one of the targets of p53. MiR-34a also targets the BCL2 and caspase 3 [27]. It has been reported that the p53, miR34a, and let-7a were significantly down regulated in BC in comparison with normal samples, however, BCL2 was up regulated in a sub population of Iranian patients [28]. MiR-328 inhibits the BCRP/ABCG2 expression which is an ABC transporter protein involved in cellular defence and multi-drug resistance. It has been observed that the tumor tissues had higher levels of miR-328 expression in comparison with normal margins among a sub population of Iranian BC patients [29]. Significant miR-132/ miR-212 downregulation and miR-22 upregulation were reported in BC samples compared with corresponding normal margins among Iranian cases. Moreover, highgrade samples had decreased expression of miR-132, miR-212, and miR-22, whereas the advanced stage (III) tumors had increased expression of these markers [30]. HOX10 transcription factor is one of the targets of mir196a2 [31]. It has been shown that the miR-196a expression in BC patients harboring C/C rs11614913 was higher than that in the cases with TT genotype. The rs11614913 polymorphism was introduced an early detection marker among a sub population of Iranian BC patients [32]. Long non-coding RNAs (lncRNA) are a class of transcripts longer than 200 nucleotides. Beside the aberration in coding RNAs, the deregulation of lncRNA is also related with tumor progression. They exerts their post-transcriptional regulatory functions through interaction with DNA, RNA, and protein molecules. CCAT2 is a lncRNA that exerts its oncogenic function through regulation of WNT/β-catenin pathway [33,34]. It has been shown that there was lower levels of CCAT2 expression in tumor tissues compared with normal margins in a group of Iranian BC patients. Moreover, significant converse correlations were also observed between CCAT2 expression, stage, and lymph node involvement in which the early stage tumors without metastatic lymph nodes had higher levels of CCAT2 expression [35]. HOX transcript antisense intergenic RNA (HOTAIR) is also a lncRNA, which binds to LSD1 and PRC2 to promote HOXD repression through histone modification and chromatin remodeling [36]. It has been shown that the HOTAIR rs920778 and rs12826786 polymorphisms were significantly correlated with high risk of   Malek-Hosseini [209] 2017 SDC1 Familial/ Sporadic BC, and rs1899663 variant also showed a significant correlation with low risk of BC among a sub population of Iranian cases [37]. Therefore, HOTAIR have probably an important role in BC progression via polycomb complexes and can be introduced as efficient prognostic and diagnostic marker for the BC. Deregulation of translation machinery and ribosome biogenesis are common features of tumor cells in post transcriptional level and epigenetic. eRF1 and eRF3 are the main terminators in eukaryotic translation process, in which the RF3 functions as an enhancer of RF1 through GTPase activity [38]. Moreover, RF3 is involved in poly (A) tail stabilization via interaction with PABP [39]. RF3 has two isoforms including RF3a and RF3b that are encoded by eRF3a/GSPT1 and eRF3b/GSPT2 genes, respectively [40,41]. The association between the GGC repeat polymorphism in exon 1 of eRF3a/GSPT1 and the potential genetic susceptibility for BC progression was evaluated among a sub population of Iranian subjects. It has been shown that the Iranian patients harboring the 12-GGC and 7-GGC alleles of eRF3a/GSPT1 have a high risk of BC. Moreover, the homozygote 7-GGC repeat can be correlated with the age of BC progression among a sub population of Iranian females [42]. The tumor progression is related with genetic and epigenetic changes. DNA methylation is one of the most important epigenetic mechanism during the neoplastic transformation. Role of DNA methylation in tumorigenesis is mainly related to the promoter hypomethylation and hypermethylation of oncogenes and tumor suppressor genes respectively. Methylenetetrahydrofolate reductase (MTHFR) regulates folate metabolism, DNA methylation, synthesis, and repair [43,44]. There is probably a correlation between MTHFR 1298AA and 677CC genotype and higher risk of BC. Moreover, the frequency of MTHFR 677 CT, TT and 1298 CT, AA variants in BC patients were lower than that in the healthy controls in an Iranian population [45]. It has been also observed that the C677T TT variant of MTHFR was significantly correlated with elevated risk of familial BC in Iranian patients [46]. Tamoxifen has been used as a treatment option for hormone responsive breast tumors. However, about onethird of patients are resistant and face with relapse [47]. The expression of DNMT1, DNMT3A, and DNMT3B DNA methyltransferases in tamoxifen resistant patients were higher than that in the tamoxifen sensitive Iranian BC patients. There were significant associations between their protein expression and tumor grade, and also DNMT3B expression and lymph node metastasis in tamoxifen sensitive group. There was significant correlation between DNMT3A and DNMT3B expression and high histologic grade in tamoxifen resistant group. Generally they showed that the DNMTs overexpression has a negative association with survival in Iranian tamoxifen-treated BC cases [48]. Another group showed that, the DNMT3B 149 C → T polymorphism was significantly associated with reduced BC risk among patients in Southern Iran [49].

Apoptosis, DNA repair, and cell cycle
Apoptosis is a critical process through the organ development and is a defense mechanism against deregulated cell proliferation and tumorigenesis. It includes several pathways such as extrinsic and intrinsic. Death receptors are the activators of extrinsic pathway which subsequently activates the caspase-8. Then, the caspase-8 initiates a cascade by cleaving caspase-3, caspase-7, or Bid [50]. A group evaluated the probable association between CASP8 6N del variant and BC susceptibility in Iranian patients. They showed that there was a significant higher frequency of del allele in the control group compared with BC cases, which highlights the anti tumori role of the CASP8 −652 6N del variant among Iranian BC cases [51]. Another group also assessed the levels of CASP8 expression in tumor in comparison with normal margins of Iranian BC patients, which showed decreased CASP8 expression in tumour tissue. Moreover, there was a significant correlation between CASP8 expression and HR status [52]. TP53 as a tumor suppressor regulates the cell cycle and apoptosis [53]. It was shown that the TP53 16-bp INS/DEL polymorphism can be associated with risk of BC among a sub population of Iranian patients [54]. There was also significant correlations between p53 expression, lymph node involvement, age, and tumor grade in Kurdish BC patients in Western Iran [55]. The mitochondrial pathway of apoptosis process is associated with BCL-2 which is also a prognostic factor in BC for its relation with defence mechanisms against the genotoxic damages [56,57]. It has been observed that there was higher levels of BCL-2 mRNA expression in tumors compared with normal margins in a sub population of Iranian BC subjects. BCL-2 over expression was positively associated with tumor grade and ER/PR, and negatively with tumor size. There was also a significant correlation between high levels of BCL-2 mRNA expression and survival rates [58]. PTEN is a tumor suppressor via the stimulation of apoptosis in tumor cells and suppression of cell proliferation. Moreover, It is a negative regulator of intracellular PIP3 and AKT/PKB signaling pathway. PTEN regulates the p53 stability through AKT-MDM2 pathway or direct interaction [59]. Moreover, PTEN-p300 complex maintains P53 acetylation following the DNA damage [60]. A significant correlation between PTEN expression and stage has been reported among a group of Iranian BC patients, in which advanced stage tumors had lower levels of PTEN expression. Moreover, the tumors without PTEN expression had lower survival rates compared with tumors with PTEN expression, which can be introduced as a prognostic marker [61]. It has been observed that the majority of BC patients had PTEN promoter methylation which was not observed among the normal cases. Moreover, there was a correlation between PTEN hypermethylation and lymph node metastasis and advanced stages tumors among a sub population of Iranian patients [62]. Another study have reported there was association between PTEN promoter methylation and BC susceptibility among a sub population of Iranian patients [63]. BIRC5 is an anti apoptotic factor and cell cycle regulator which is only expressed at the embryonic stage. It has been reported that the BIRC5 is associated with early onset BC progression, in which there was a significant correlation between BIRC5 up regulation and Iranian patients who were younger than 40 years [64]. DNA repair and cell cycle regulation are important processes related to the cancer susceptibility and chemotherapeutic responses. Cell proliferation is regulated by signaling pathways to limit the proliferation in adequate ratio. The cell cycle is regulated by several factors such as cyclins, CDK, and CDKIs. G1 to S progression is related with cyclin D and CDK 4/6, which is repressed by INK4 family [65]. P16INK4A is belonged to the INK4 family [66], and its deregulation has been reported in different tumors [67]. The p16INK4A is responsible for the G0 and G1 cell cycle arrest through suppression of cyclin D-CDK4/6, and it also has a critical role as an inhibitor of pRb phosphorylation [68,69]. It has been reported that there was a correlation between p16INK4A promoter methylation and primary BC, in which the hypermethylated cases were in primary tumors with stages of I, II and III. Therefore, they concluded that the p16INK4A promoter methylation can be introduced as an efficient marker of early detection among a sub population of Iranian sporadic BC patients [70]. Cyclin D1 is one of the main cell cycle regulators that binds with CDK4/6 to phosphorylate pRb through the cell cycle. The ataxia telangiectasia-mutated (ATM) is also involved in DNA damage recognition [71]. A study has evaluated the ATM and cyclin D1 expression in a sub population of Iranian BC patients. They showed about 21.6% of cyclin D1 under-expression which were mostly in HER-2-positive tumors. In the case of ATM, about 50% had ATM under-expression. The levels of ATM expression in normal margins were lower than normal controls which can be related to the epigenetic regulation of ATM expression. Therefore, they concluded that the ATM and cyclin D1 under and over expressions respectively can be associated with tumor progression among the Iranian BC patients [72]. It has been observed that there was a significant difference of D1853N ATM polymorphism between patients and healthy controls, introducing this polymorphism as a diagnostic marker among a sub population of Iranian cases [73]. Cyclins and CDKs activate the E2F transcription factor through inhibition of retinoblastoma protein. Cyclin E/CDK2 complex is associated with DNA replication and chromosomal stability [74]. A correlation between CCNE1 rs1406 C/A polymorphism and elevated risk of BC has been shown among an Iranian population in southeast of Iran [75]. DBC2 inhibits cell proliferation via CCND1 down-regulation [76], and is associated with ubiquitination, cell cycle, and apoptosis [77]. Methylation status of DBC2 promoter sequence was assessed in blood and BC tissue samples of an Iranian sub population, showed a correlation between aberrant DBC2 methylation and tumor progression in sporadic cases [78]. The AURKA is involved in cell cycle regulation, G2/M transition, centrosome division, and genomic instability. It has been observed that the BC cases had higher frequency of Ile/Ile variant at F31I locus of AURKA in comparison with the healthy cases. Moreover there was significant correlation between Ile/Ile F31I polymorphism and stage of BC among a sub population of Iranian patients [79]. ERCC1 as a member of nucleotide excision repair system is vital for the genome integrity. Association between ERCC1 C8092A and BC was assessed and showed a correlation between T/T of ERCC1 genotype and increased risk of BC among a sub population of Iranian patients [80]. XRCC1 is one of important factors in base excision repair to maintain the DNA against the carcinogens and oxidative species [81,82]. It has been reported that the Arg194Trp XRCC1 polymorphism was more frequent among patients compared with normal cases. Moreover, this variant was also correlated with stage. They introduced the Trp194 allele of XRCC1 as a risk factor among Iranian Kurdish BC patients [83]. Flap endonuclease 1 (FEN1) tumor suppressor is a critical factor in DNA replication, repair, and genome stability [84]. It is involved in base-excision repair and formation of Okazaki fragments during DNA replication [85]. There was a correlation between FEN1 haplotypes and BC susceptibility among a sub population of Iranian patients in Southeastern Iran [86]. BRCA family mediates double stranded DNA break repairs and chromatin remodeling [87]. EMSY is an inhibitor of the BRCA2 transactivity. It has been shown that the majority of a group of Iranian BC patients had EMSY protein expression. The sporadic cases had higher levels of EMSY expression compared with familial cases. There was a direct correlation between EMSY expression, larger tumor size, higher relapse, and increased lymph node involvement [88]. Homologues Recombination (HR) is one of the main processes during the DSB repair. RAD51 is a key factor in HR regulation. DSB activates the ATM that subsequently phosphorylate and activate the BRCA1 which is associated with RAD51 [89,90]. It has been shown that there were significant ATM, RAD51, and BRCA1 down regulations in BC tissues compared with normal margins in a sub population of Iranian subjects. ATM under expression was significantly correlated with tumor stage, vascular invasion, and malignant behavior. Down regulation of RAD51 was also significantly correlated with lymph nodes metastasis, age, and grade of tumor [91]. BRCA1 is a nuclear phosphoprotein involving in DSB repair, cell cycle regulation, and apoptosis [92,93]. There was significant correlations between BRCA1 downregulation, high grade, and lymph node metastasis. Low levels of BRCA1 expression were observed in majority of triple-negative and luminal tumors in comparison with normal tissues among a sub population of Iranian cases [94]. The altered BRCA1 expression is more frequent in primary BC, and loss or decreased BRCA1 expression in high-grade were more frequent in comparison with differentiated tumors. Moreover, there was a significant inverse association between the levels of BRCA1 and CD44 expressions in Iranian BC patients [95]. APOBEC3 is a cytidine deaminase involved in RNA editing, cell cycle regulation, and repression of the retrovirus replication. It has been observed that the APOBEC3 deletion had a significant correlation with increased risk of BC among a sub population of Iranian cases [96]. The FANCA is a tumor suppressor related to DNA repair [97]. Duplication in FANCA promoter sequence can be associated with increased risk of BC among Iranian cases [98]. Telomeres are repeated sequences at the end of linear chromosomes to maintain them against the loss of genetic information through the cell proliferation [99]. Telomere shortening leads the cells up to the senescence and apoptosis, which is an intra cellular defence mechanism against genomic instability and tumorigenesis [100]. hTERT is the main subunit of telomerase complex to maintain the telomere stability [101]. The hTERT rs2736098 variant increased the risk of BC in a subpopulation of Iranian patients in southeast Iran [102].

Immune response
The immune system is able to find tumor cells through tumor-specific antigens or other induced factors by cellular stress [103]. However the microenvironment changes by tumor cells can be resulted to tumor escape from immune response [104]. HLA system has important role in tumor antigen presentation for the immune system [105,106]. HLA-G has immunosuppressive functions during the tumor progression [107]. A mutational analysis on a sub population of BC patients in northwestern-Iran showed a slightly higher allele frequency of HLA-G 14 bp deletion among patients. It seems that the high HLA-G expression induces tumor cells toward progressive stages. They concluded that the HLA-G 14 bp InDel polymorphism is probably a risk factor of BC progression among this population [108]. HLA-l and 2 present the antigens to the CD8 + T and CD4 + T cells, respectively. HLA-DRB1/DQB1 and HLA-DQA1 are belonged to the HLA-2 beta and alpha chains paralogues, respectively. It has been shown that the HLA-DQA1*0301 variant is correlated with elevated risk of early onset BC progression. Whereas, the HLADQA1*0505 and HLA-DQB1*0602 variants were associated with decreased risk of BC among young patients. They concluded that the HLA-II polymorphisms are risk factors of BC among a sub population of Iranian patients [109]. It has been observed that there was a significant association between HLADRB1* 12 allele and risk of BC among a sub population of Iranian patients [110]. Interleukin-13 is mainly secreted by activated TH2 lymphocytes and Natural Killer T cells [111,112]. It has been shown that there was a correlation between CCA and ACA haplotypes of IL13 gene and risk of BC. Moreover, −1512C and −1055C alleles are involved in BC progression among a sub population of Iranian patients [113]. IL-27 is a heterodimeric cytokine comprising of EBI3 and P28 subunits that induces the CD4 + T cell proliferation through activation of STAT1 and 3 transcription factors [114,115]. Moreover, IL-27 in association with TGF-β has a critical function in generation of IL-10-producing Tr1 cells [116]. IL-23 also regulates the memory CD4 + T cells proliferation through STAT3 activation [117,118]. It has been reported that there was significantly higher levels of IL-27 and IL-23 expressions in BC peripheral bloods compared with normal cases, which can be introduced for the immunotherapy in Iranian BC patients [119]. CCR7 chemokine receptor is a G-protein coupled receptor involved in migration of immune cells to the lymphoid organs. It has been observed that the CCR7 expression was significantly associated with lymph node involvement, garde, and stage. Therefore, they introduced CCR7 as a prognostic marker among a sub population of Iranian BC patients [120]. CCL22 chemokine is produced by macrophages and DCs following the activation with microbial infection or CD40. The CCL22 production is regulated by Th1/2 cytokines such as IL-4, IL-5, and IFN-γ [121]. Therefore, CCL22 is associated with Th2-and Treg cells through inhibition of the immune responses against microbial infections and cancer cells [122]. It has been shown that the BC patients had significantly higher serum levels of CCL22 compared with healthy cases. There was also a direct correlation between the levels of CCL22 expression and stage of tumor. The frequencies of CC and C genotypes at rs223818 of CCL22 were significantly higher in BC patients in comparison with normal cases. Carriers of such genotypes also had higher serum levels of CCL22 compared with cases harboring GG or G variants, which highlighted the role of SNP rs223818 in BC risk via upregulation of the CCL22 [123]. CXCL12 is a chemokine that binds to its receptor (CXCR4) to activate MAPK1/MAPK3 through intracellular calcium ion which is involved in migration and adhesion. It has been shown that there was a correlation between HER2 and CXCR4 expression in a group of Iranian BC patients in which the HER2 positive tumors significantly had higher levels of CXCR4 expression in comparison with the HER2 negative tumors. They suggested the HER2 as probable inducer of CXCR4 among BC patients [124]. T-regulatory cells express cytotoxic T-lymphocyte antigen-4 (CTLA-4) and FoxP3 transcription factor [125,126]. It has been shown that there was higher levels of FOXP3 and CTLA-4 mRNA expressions in blood samples of early stage and low grade BC patients compared with healthy subjects in an Iranian population. They highlighted the importance of regulatory T cells inhibitory mechanisms through the tumor progression which prepares opportunity for the tumor cells to hide from immune system [127]. A significant correlation has been also reported between FOXP3 and OX40 expressions and tumor stage in peripheral blood samples of Iranian BC patients [128]. Obesity is one of the main socio-epidemiological health problems globally and a negative BC prognostic factor [129]. The high ratio of adipose tissue increases the risk of tumor relapse in BC [130]. Negative effect of obesity on BC survival is associated with inflammatory cytokines such as leptin and IL-6 which are secreted by adipose tissue [131]. Leptin is an adipocytokine which is important in tumor progression through the activation of AKT signaling pathway [132]. Leptin is involved in cell proliferation, migration, and survival via regulation of various factors such as AP1, ERK2, and MAPK [133]. It has been reported that the leptin G-2548A (rs7799039) polymorphism was significantly correlated with risk of BC among a sub population of Iranian patients. They concluded that the females harboring AA genotype have an earlier menarche age and higher risk of breast cancer [134]. It has been also shown that the obese BC had higher serum leptin concentration compared with non-obese Iranian BC cases. Therefore, the inflammation has a key function in distant metastasis, since the adipose is the main source of leptin the weight loss through the physical activity can be an efficient option to decrease the negative role of leptin in BC progression [135]. LEP-2548G/A frequency was evaluated and showed the Iranian patients harboring LEP-2548G variant had a noticeable increased risk of BC [136]. Chronic stress is involved in tumorigenesis through secretion of neurotransmitters like during lymphocyte migration and angiogenesis. Dopamine is a neurotransmitter which has proliferative effects and induces the T cells for secretion of anti immune cytokines [137]. Dopamine receptors are from G-protein family which play important roles to induce or suppress intracellular cAMP [138]. The stress can alter dopamine receptor expression. Dopamine receptors (DRD1-DRD5) are categorised into D1 (stimulatory receptors) and D2 (inhibitory receptors) families [138,139]. It has been reported that there was significantly higher levels of DRD2-DRD4 in a group of Iranian BC cases compared with healthy controls [140]. The transforming growth factor (TGF)-β is a pleiotropic cytokine expressed in several cells including peripheral blood mononuclear cells and platelets that plays critical roles in regulation of cell proliferation, apoptosis, and angiogenesis. The TGF-beta superfamily are involved in cell proliferation through activation of SMAD transcription factors. It has an anti cancer function in normal cells, whereas it promotes the tumor progression in advanced stage tumors. It has been observed that there was a significant difference in frequency of TGFβ1 "GTG CCG C" haplotype between BC and healthy controls among an Iranian population [141]. Moreover, it has been shown that there was significant correlations between TT and AA genotype of −509T/T and −800G/A TGF-β respectively and BC risk among a sub population of Iranian patients [142].

Kinases and growth factors
The human epidermal growth factor 2 (HER-2) is a transmembrane tyrosine kinase receptor which is involved in cell proliferation. It was shown that fourteen out of 60 (23%) had HER-2 overexpression in a subpopulation of malignant BC patients in central Iran [143]. There was also an inverse correlation between HER-2 expression and ER/PR levels in a group of Iranian BC cases [144].
Another study on BC patients in Kermanshah province of Iran have shown a significant association between age, stage, and HER2-neu expression. Therefore, they introduced the HER2-neu as an efficient prognostic and diagnostic marker in that population [145]. ErbB4 is a tyrosine kinase receptor belonging to the EGFR subfamily. Frequency of rs11895168 variant was analyzed in a subpopulation of central Iranian BC patients, and observed that none of the normal cases harbored the rs11895168 C allele. They introduced a significant correlation between C allele and higher risk of breast cancer, which can be related to the miR-1276 binding changes [146]. Another study has shown that there was a significant correlation between rs13423759 allele C of ERBB4 and increased risk of BC which can be associated with role of this allele in strengthening the miR-548as-3p/ErbB4 interaction. The rs13423759 allele A is also significantly associated with decreased risk of breast cancer. Therefore, they highlighted the role of rs13423759 as an efficient diagnostic biomarker among a group of Iranian BC patients [147]. The rs1836724 polymorphism of ErbB4 was correlated with risk of BC, in which the T allele carriers had decreased expression of ErbB4 and higher susceptibility for BC among a group of Iranian patients [148]. NGX6 has an epidermal growth factor (EGF)-like domain which functions as a negative growth factor through apoptosis induction and cell proliferation arrest [149,150]. NGX6 exerts its negative role on cell cycle via down regulation of cyclin D1, A, and E [151]. It has been shown that the levels of NGX6 mRNA expression in Iranian BC tumors were significantly lower than that in the normal samples [152]. Neuropilin-1 (NRP1) as a receptor for VEGF and PlGF in endothelial cells plays key roles in cell survival, angiogenesis, and EMT [153]. It has been shown that there was significant correlation between NRP1 expression and lymph node involvement in which the patients with metastatic lymph nodes had higher levels of NRP1 expression in comparison with the lymph node-negative cases. Moreover, the tumors with higher stages and grades had higher levels of NRP1 expressions. Therefore, they introduced NRP1 as a metastatic diagnostic marker among a sub population of Iranian BC patients [154]. Insulin-like growth factor 1 (IGF-1) is one of the activators of the AKT pathway, which is a cell proliferation inducer and cell death suppressor. There is a polymorphic cytosine-adenine (CA) repeat located in 1 kb upstream of the IGF-1 start site which is associated with cancer [155]. Length of this dinucleotide repeat was assessed in a group of Iranian BC patients. They showed that the patients harboring an allele longer than 19 and shorted than 20 have higher and lower risks of breast cancer, respectively. Therefore, there was a correlation between BC and length of the IGF-1 CA repeat which seems related to the ethnicity [156]. Phosphatidylinositol-3 kinase (PI3K) is a heterodimeric protein comprising of regulatory and catalytic subunits encoded by PIK3R1 and PIK3CA respectively, which is involved in apoptosis and cell proliferation [157]. Mutational analysis of PIK3CA catalytic subunit was done in a group of Iranian BC patients, and showed 45% of 80 cases had mutated PIK3CA that was higher than the recent reports. Majority of mutations were observed in three hotspots belonging to exon 9 and 20. There was a significant correlation between PIK3CA mutation and low grade tumors [158]. KRAS oncogene is a member of the small GTPase superfamily, which is also associated with PI3-kinase signaling pathway. It has been reported that the rs61764370 and rs712 polymorphisms in KRAS increased and decreases the BC risk in a sub population of Iranian patients [159]. SPHK1 catalyzes the phosphorylation of sphingosine to sphingosine-1-phosphate (S1P), which is a ligand of G protein-coupled receptors and regulates cell proliferation. SPHK1 exerts its inflammatory and anti apoptotic roles through several signaling pathways such as TNFalpha and NF-kappa-B. It has been observed that there was a positive correlation between SPHK1 expression in ER/PR negative tumors compared with positive BC cases. Moreover, there was an association between BMI and SPHK1 expression among the Iranian ER negative BC patients [160]. The Class I PI3K family phosphorylates PIP2 to PIP3 which subsequently activates AKT as a regulator of cell proliferation, adhesion, and migration [161]. PIK3CA is the catalytic subunit of class I PI3Ks [162] that is deregulated in various cancers [163]. A highly polymorphic intronic GT dinucleotide repeat in PIK3CA had direct correlation with risk of BC in a group of Iranian patients, in which the harbors of short allele (17 repeats) and long allele (more than 16 repeats) had significantly higher and lower risk of BC progression, respectively [164]. OCC1 is involved in differentiation of adipocytes. APPL2 is one of the small GTPase RAB5A/Rab5 effectors, which are associated with endosomal membranes and can interact with the NuRD/MeCP1 complex which is required for cell proliferation. Levels of OCC1 and APPL2 expression were assessed in BC tissues, and showed significant under and over expressions of APPL2 and OCC1 in tumors compared with corresponding normal margins, respectively. Therefore, they suggested that the OCC-1D variant suppresses the APPL2 and probably regulate the cell cycle through PI3K/AKT signaling pathway [165]. RIP1K and RIP3K are belonged to RIPK family and involved in inflammation and cell death, in response to tissue damage. It has been observed that the benign and malignant tumors had increased levels of RIP1K compared with normal tissues. Whereas the malignant tumors had decreased levels of RIP3K expression compared with normal and benign samples. Moreover the RIP1K and RIP3K had decreased and increased expressions among old Iranian BC patients with malignant tumors, respectively. The tumors with higher grade and larger size had increased levels of RIP1K expression. Whereas the RIP3K had decreased expression in malignant tumors with grade III [166].

Intracellular defense systems and drug resistance
Ubiquitin proteasome machinery is one of the main intracellular defence systems which is involved in degradation of aberrant and defective proteins. Moreover, it regulates the stability of cell cycle regulators. Ubiquitinproteasome pathway involves several members such as Ubiquitin-activating (E1) and conjugating (E2) enzymes (E1) and ubiquitin ligases (E3). Since the lifetime regulation of cellular regulatory proteins including cell cycle and transcription factors is governed by the ubiquitin process, ubiquitine proteasome deregulation can be resulted in tumor progression [167]. UBE2Q2 expression was assessed in Iranian BC cases, and showed that the majority of tumor tissues had UBE2Q2 over expression compared with their corresponding normal margins. They introduced UBE2Q2 as a novel diagnostic marker among the Iranian BC patients [168]. Glutathione S-transferases (GSTs) as the cytoprotective factors detoxifie the carcinogens and chemotherapeutic drugs through the glutathione binding. Therefore, the human polymorphisms affects the individual ability to respond toward drugs and stresses. Although, the GSTM1 is expressed in a wide spectrum of tissues such as liver, stomach, and breast, GSTT1 enzyme is mainly observed in liver and erythrocytes [169]. The association between GSTM1, GSTT1 and GSTP1 genotypes and BC were assessed in a subpopulation of Iranian cases. It was reported that the GSTM1 null and GSTP1 Ile105Val genotypes were significantly correlated with the increased risk of BC. It seems that the cases with GSTM1 null, GSTT1 null, or GSTP1 Ile105Val Genotypes have low GST detoxification capacity, which can be resulted in higher accumulation of toxic compounds in such cases [170]. GSTO1 A140D and GSTO2 N142D polymorphisms were also associated with high risk of BC in a sub population of Iranian patients [171]. Cytochrome 450 (CYP450) is involved in oxidative catalysis of anticancer drugs and hormone synthesis. Therefore, assessment of CYP450 polymorphisms is required in anticancer therapy outcomes. CYP1A1 is involved in metabolism of xenobiotics and reactive epoxide production from aromatic compounds and steroid hormones [172]. It has been reported that the heterozygote carriers of 2455G (462Val) of CYP1A1 had significantly higher risk of BC in comparison with other genotypes. Moreover, heterozygote variant was more frequent in pre-menopausal patients, highlighting the role of this variation in early onset of BC among a group of Iranian patients in Eastern Azerbaijan [173]. Non functional CYP2D6*4 and 10 alleles were also frequent among BC patients in Northern Iran [174]. Tamoxifen is treatment option of ER-positive BC in premenopausal women [175,176]. However, some patients do not response to the tamoxifen and die from tumor recurrence [177]. Tamoxifen is metabolized by cytochrome P450 2D6 (CYP2D6) in endoplasmic reticulum [178]. It has been shown that the HER2-neu positive Tamoxifen treated patients with CYP2D6 polymorphism had lower risk of tumor relapse in a sub population of Iranian BC patients [179]. Heat shock proteins (HSPs) are induced in response to anticancer drugs to maintain the cells in lethal conditions. Their cytoprotective function involve the correction of misfolded proteins and anti-apoptotic functions. Hsp27 regulates ability of cells to respond to the heat shock and oxidative stresses [180,181]. Moreover, Hsp27 has an anti-apoptotic function through suppression of caspasedependent apoptosis [182,183]. It has been observed that there was a significant higher levels of serum anti-Hsp27 antibodies in BC compared with healthy Iranian cases [184]. The ATP-binding cassette (ABC) family regulates the drugs passage across cell membranes, which alters the kinetics of drugs. The ABCG2 is related with drug resistance in tumor cells and xenobiotic detoxification [185,186]. It has been shown that there was a significant association between A allele of ABCG2 C421A variant and response to the Anthracyclines and Paclitaxol, and also increased BC risk in a group of Iranian BC patients [187].

Cell adhesion, cytoskeletal factors, and ECM
Alterations in adhesion properties facilitate the tumor cell detachment from the primary tumor bulk and formation of secondary tumors in distant organs. The CEACAM subfamily are belonged to the CEA family which are involved in cell adhesion, differentiation, apoptosis, and angiogenesis [188]. Levels of CEACAM19 expression was evaluated among a group of Iranian BC patients, showing a higher levels of CEACAM19 mRNA expression in tumors compared with normal tissues [189]. The MUC1 is a mucin glycoproteins expressing on the apical surface of epithelial breast cells and is involved in cell adhesion [190]. CA15-3 serum assay determines the MUC1 protein. CEA is one of the members of immunoglobulin family and is associated with tumors [191]. Serum levels of CEA and CA15-3 were assessed in BC compared with healthy controls among Iranian cases and showed a significant increase of CA15-3 and CEA in BC in comparison with normal cases. There were significant correlations between CEA and grade of tumor, and also between CA15-3 and tumor size [192]. MUC1 overexpression was also correlated with BC progression and metastasis among Iranian cases [193]. CDH1 tumor suppressor is a cell-cell adhesion glycoprotein which is associated with cell differentiation and polarity [194,195]. A significant correlation between CDH1 −160C/A polymorphism and BC susceptibility was reported among a group of Iranian Kurdish patients in which the A allele of CDH1 −160C/A variant was significantly associated with BC metastasis [196]. CD44 is involved normally in cell adhesion, migration, angiogenesis, and apoptosis [197]. The CD44 protein increases cell proliferation and migration in BC through an interaction with hyaluronan. It has been reported that there was a significant correlation between A > G intronic polymorphism and grade 3 of tumors in Iranian BC patients, which can be related to generation of a new SC35 binding site and modulation of splicing process. Low grade (1/2) tumors were also associated with AA genotypes [198]. Epithelial cell adhesion molecule (EpCAM) is a membrane glycoprotein associated with cell-cell adhesions and is involved in cell proliferation in different cancers [199]. EpCAM induces the Lef1 transcription factor through presenilin-2 to bind with promoter sequences of target genes such as c-myc and cyclin A/E genes [200]. C-myc protooncogene has important functions in various cellular processes such as cell cycle, angiogenesis, and apoptosis [201]. It has been observed that there was a significant correlation between EpCAM/c-myc over expression and grade III among a group of Iranian BC patients [202]. Matrix metalloproteinases are involved in ECM degradation to facilitate the cell migration and metastasis during the tumor progression. CD10 is a zinc-dependant metalloproteinase that is highly expressed in kidney and lung tissues and participates in ECM degradation [203]. Regulation of cell adhesion facilitates the tumor metastasis. Therefore, tumor cells regulate the cell adhesion through proteolytic enzymes or stimulation of fibroblasts to secrete such enzymes. It has been observed that there were inverse significant associations between stromal CD10 expression, larger tumor size, metastatic lymph nodes, and grade among Iranian BC patients [204]. Caveola are structural proteins in formation of caveola which are mainly expressed in endothelial, adipose, and smooth muscle cells [205]. CAV-1 is associated with tyrosine phosphorylation and observed in a wide spectrum of cells [206]. Role of several SNPs in CAV-1 was evaluated in BC progression among a group of Iranian patients, and showed that the cases harboring homozygous AA genotype in T29107A, G21985A, and G14713A polymorphisms have higher risk of BC progression compared with the normal cases. In contrast, the C521A homozygous variant had a protective role, and it had a significant correlation with BMI. Therefore, they introduced C521A, G14713A, G21985A, and T29107A as efficient diagnostic markers of BC progression among Iranian cases [207]. The SYCP3 is a structural component of the synaptonemal complex associated with recombination and meiotic chromosome segregation. There was a noticeable SYCP3 over expression in tumors compared with normal tissues. They introduced the SYCP3 as a CTA in BC which can be used as an early detection marker among Iranian BC patients [208]. Syndecan-1(SDC1) is an ECM receptor belonging to the heparan sulfate proteoglycan family. The extracellular domain binds to the heparan and chondroitin sulfate whereas the cytoplasmic domain is associated with cytoskeleton during the cell migration. It has been reported that there was elevated serum levels of SDC1 in BC in comparison with the healthy cases, which was also positively associated with size of BC tumors among Iranian patients [209].

Cell signaling and self renewal
AMP-activated protein kinase (AMPK) is a cellular energy sensor, which is activated with elevated AMP/ ATP ratio to induces the ATP generation. AMPK also acts as a tumor suppressor through cell cycle regulation. AMPK activates p53, which induces G1/S cell cycle arrest in response to the lack of glucose. Moreover, AMPK blocks CDKs via phosphorylation of p27 [210]. Levels of cAMP and cGMP as the second messengers are regulated via cyclic nucleotide phosphodiesterases (PDEs) [211]. PDE5 and PDE9 regulate the cGMP, which activates various downstream targets such as PKGs and ion channels [212]. Significant PDE5 and PDE9 over expressions has been observed in malignant tumors compared with benign tumors and normal margins. Moreover, there were significant correlations between phosphodiesterase expression and stage, grade and lymph node metastasis. Furthermore, there was a converse correlation between levels of expressions and age of patients. They concluded that the Iranian BC patients with higher levels of phosphodiesterases may have higher risks of malignant tumors [213]. Cyclic guanosine monophosphate (cGMP) participates in various signaling pathways through Protein Kinase G (PKG). PKGI and PKGII are serine/threonine kinases mediating the effects of cGMP [214]. It has been reported that the BC tissues had lower levels of PKG expression in comparison with the normal samples, highlighting the PKG as a tumor suppressor in a group of Iranian cases [215]. Notch is one of the main signaling pathways involved in cell proliferation, migration, and self renewal [216]. This pathway is a cell to cell contact process which is activated via the ligand binding with a family of trans membrane receptors (NOTCH1-4).
Ligand binding releases the intracellular domain of Notch (ICN) into the cytoplasm which eventually enters to the nucleus where it regulated the expression of NOTCH target genes [217]. Levels of NOTCH1 mRNA expression as one of the NOTCH receptors was assessed in BC patients and showed a significant increase in invasive ductal types compared with other histopathological types. They introduced the NOTCH1 as a prognostic marker of IDC among Iranian patients [218]. WNT pathway is also another important process to maintain a normal balance between cell differentiation and proliferation. It is initiated through the cell surface receptors resulting in activation of cytoplasmic b-catenin which activates the LEF/ TCF/PYGO2 transcriptional complex [219]. The WISP1 is belonged to the CTGF family and involved in WNT pathway. It also functions as an anti apoptotic factor through up regulation of BCL-X and inhibition of the P53 mediated apoptosis following the DNA damage via AKT activation. It has been reported that the metastatic BC tissues had significantly reduced levels of WISP1 mRNA expression in comparison with the normal subjects. Moreover, WISP1 expression was correlated with age and tumor size [220]. Cancer stem cells (CSC) are a small subpopulation of cells with self renewal and chemo-radiotherapeutic resistance abilities [221]. Aldehyde dehydrogenase 1 (ALDH1) has been suggested as a marker of mammary stem cells and prognostic marker [222,223]. Epithelial and stromal expression of ALDH1 was analysed to assess the probable correlation between ALDH1 and clinicopathological features of breast tumors in a group of Iranian patients. They observed a significant converse correlation between ALDH1 and BRCA1 expressions. Moreover there was a correlation between ALDH1 +/BRCA1 low expression phenotype and high grade tumours [224]. BC cases are routinely undergone the different treatments such as surgery, radiotherapy, and chemotherapy. Hormone therapy is also one of the treatment options for the BC cases, in which the tamoxifen targets the positive estrogen receptor tumors through apoptosis induction and ER signaling inhibition [225,226]. However, it has been observed that a noticeable ratio of tamoxifen treated cases have tumor relapse [227]. PAX2 is one of the members of PAX family of transcription factors which are involved in embryonic development of different organs [228]. PAX2 protein expression was assessed in tamoxifen treated BC patients, and showed that there was a correlation between PAX2 over expression and better survival. Moreover, the tamoxifen responsive cases had significantly higher levels of PAX2 expression compared with tamoxifen resistant patients [229]. MYC is an oncogene involving in various cellular processes, such as cell proliferation and apoptosis [230]. NDRG family is one of the targets of MYC which has four members including NDRG1-4 [231]. NDRG3 has high levels of expression in various tissues including prostate, testis, brain, heart, and kidney [232]. It has been suggested that the NDRG3 has a tumor suppressor role in BC patients, in which the triple negative and advanced stage BC tumors had the lowest levels of NDRG3 expression. Moreover, there was an inverse correlation between NDRG3 expression and survival rate, in which the higher levels of NDRG3 expression resulted in better survival rates among the Iranian BC patients [233].

Vitamin D and steroids
Steroid hormones exert their regulatory role on transcription of target genes through binding with steroid receptors including estrogen receptor (ER) and progesterone receptor (PR). ERα has an important function in BC progression through cell growth and proliferation induction. It has been shown that there was an increased expression of ERα in the levels of mRNA and protein in a sub population of Iranian BC tissues in comparison with normal margins and it was associated with tumor size [234]. A mutational analysis has also shown that the presence of the ERα gene A908G mutation among a sub population of Iranian cases with invasive BC. Moreover, they observed a significant correlation between ER-α A908G mutation and familial BC [235]. P450 family is one the most important genes which are involved in steroid hormones metabolism. It has been observed that there was a significant association between rs10,046 polymorphism of CYP19 gene and BC among a group of Iranian patients [236]. TNF-α is involved in estrogen synthesis [237,238] and down-regulation of ER [239]. Steroid hormone receptor is a prognostic factors to identify sensitive patients for hormone therapy, and also is associated with higher survival rates [240]. It has been shown that there was a significant correlation between TNFA −308 A/G polymorphism and levels of ER expression in a sub population of Iranian BC patients. Moreover, the cases harboring TNFA2 and TNFB2 alleles had lower risk of PR positive tumors in comparison with the TNFA1/A1 and TNFB1/B1 homozygous genotypes. They highlighted the correlation between immune and endocrine systems in BC progression, in which higher levels of TNF-α inhibits the PR expression [241]. TNFSF4 is also a glycoprotein belonged to the TNF family which is expressed on activated T cells and various antigen presenting cells such as dendritic cells and B cells [242][243][244]. It has been shown that the rs3850641 G allele is significantly correlated with an increased risk of BC among Iranian cases [245]. The inverse relation between 25-hydroxyvitamin D serum levels and BC can be related to the apoptotic and antiproliferative roles of vitamin D [246]. VDR is a nu-clear receptor for the steroid hormones participating in regulation target genes which are involved in cell proliferation [247,248]. VDR/RXR complex binds to the promoter sequences of target genes to regulate their expression [249]. The length of VDR poly(A) microsatellite was assessed in a sub population of Iranian BC cases, and showed a significant correlation between increased BC risk and VDR poly(A) L variant [250]. This region is associated with mRNA stability, which highlights the key role of poly(A) microsatellite in regulation of VDR mRNA expression. CYP27B1 is involved in conversion of 25-hydroxyvitamin D3 to 1,25-dihydroxyvitamin D3 which finally regulates the transcription of target genes through VDR/RXR complex [251,252]. This complex regulates the cell growth, apoptosis, and inflammation [253]. Finally the vitamin D pathway is terminated via CYP24A1 [254]. It has been reported that the levels of CYP27B1 mRNA expression in normal was higher than that in the tumor samples. In contrast, the tumors had significantly elevated levels of CYP24A1 expression compared with normal tissues. The CYP24A1 and CYP27B1 were significantly correlated with age and stage, respectively [255]. Another groups have also shown that there was a significant correlation between VDR BsmI polymorphism and elevated BC risk [256,257]. CYP17 is a steroid monooxygenase which is involved in synthesis of glucocorticoids and sex steroids [258,259]. It has been shown that the CYP17 MspA1 Polymorphism is probably age specific and is associated with increased risk of early onset BC among Iranian subjects [260].

Cancer testis antigens
Cancer-testis antigens (CTA) are commonly expressed in testis and also with a lower ratio in ovarian germ cells [261]. Tumor cells also express antigens which are restricted to the germ cells. CTAs participate in EMT, self renewal, and tumor invasion [262]. AKAP as a group of CTAs has critical roles in sperm function which exerts their role through the protein kinase A (PKA). The levels of AKAP3 mRNA expression was assessed in tumor compared with normal margins of a sub population of Iranian BC cases, and showed a significant lower expressions in tumors compared with normal tissues. However, there was an inverse correlation between the levels of AKAP3 expressions and tumor sizes and stages, highlighting that as a probable inhibitor of proliferation [263]. WW-binding protein 2 (WBP2) as a CT antigen also regulates the expression of ERα/PR target genes through activation of ERα/PR expression in BC [264]. It has been demonstrated that the WBP2NL is probably associated with cell proliferation in embryonic stem cells and tumor cell lines [265]. WBP2NL over expression was also observed among the Iranian BC patients, highlighting that as a novel prognostic factor [266]. It has been reported that there was a significant correlation between MAGE-1 expression and lymph node involvement, and also between positive nuclear MAGE-1 expression and tumor size among Iranian BC patients [267]. Expressional analysis of TSGA10 CTA was performed in Iranian BC patients and showed its expression in 70% of patients while the only 12% of patients had immune response against TSGA10 [268]. ODF is involved in sperm tail maintenance and centrosome matrix [269] which facilitates the tumor cell proliferation. The RHOXF2 is associated with cell transformation [270] and cell to cell contacts [271]. The ODF4 and RHOXF2 cancer testis antigens had overexpression in tumors compared with normal marines in a sub population of Iranian BC patients [272].

Mitochondrial DNA
Mitochondria as the powerhouse in the cells are involved in redox homeostasis, innate immunity, and apoptosis. Aberrant cellular energetics is one of the most critical risk factors during the tumor progression. Mitochondrial DNA mutations and enzyme defects are the main dysfunctions which cause the cellular energetic problems [273,274]. The mitochondrial genome encodes various factors of respiratory chain and mitochondrial non coding RNAs [275]. ND4 as one of the respiratory chain components in mitochondria complex I (initiation of electron transport chain) is involved in various cancers [276]. Through the mitochondrial respiratory chain, reactive oxygen species (ROS) are generated due to the electron transfer, therefore it is really important to maintain the production and degradation of ROS, since there is not any protective factor in mitochondrial genome [277]. Mutational analysis of ND4 in a sub population of Iranian BC patients showed that there was a correlation between ND4 alteration and BC [278]. The mtDNA alterations are associated with the increased oxidative stresses and apoptotic resistance [279,280]. The ATPase6 is belonged to the complex V genes, and contributed to the preservation of mtDNA [281]. The high frequency of ATPase6 mitochondrial mutations were observed in BC patients and highlighted the importance of mitochondrial gene variants in BC progression through modulation of metabolism, and can be introduced as molecular biomarkers among Iranian BC cases [282].

Conclusions
For the first time in the current review we have summarized all of the recent significant genetic markers among Iranian BC patients. It seems that the non-coding RNAs, epigenetic modifications, and immune responses are the most common reported cell and molecular processes which are involved in tumor progression among Iranian BC patients. Moreover, we concluded that there were fifteen out of 117 reported genes which were assessed in higher numbers of BC patients and controls and they can be considered as the main high risk genes among Iranian BC patients (Table 2). Indeed, it is believed that the current review will be useful to pave the way of introducing a population based diagnostic panel markers for the early detection of BC among Iranians. Moreover, current review clarifies the molecular and genetic bases of BC progression in this population.