Identification of miR-10b, miR-26a, miR-146a and miR-153 as potential triple-negative breast cancer biomarkers

Background Familial triple-negative breast cancers are often linked to mutations in the BRCA1 tumor suppressor gene. In sporadic triple-negative breast cancers BRCA1 is frequently inactivated at the transcriptional level, and it has been reported that this inactivation may be brought about by promoter methylation. More recently, it was found that BRCA1 may also be regulated at the post-transcriptional level by miRNAs. Here, we explored the expression of putative BRCA1-regulating miRNAs in sporadic human triple-negative breast cancer cells. Methods Nine sporadic human breast cancer-derived cell lines and one benign breast epithelium-derived cell line were assessed for their hormone receptor, growth factor receptor and cytokeratin status by immunocytochemistry. The expression of 5 selected miRNAs predicted to target BRCA1 was assessed using qRT-PCR in the 10 cell lines. In addition, expression profiles of 84 known breast cancer-associated miRNAs were established in these 10 cell lines using PCR Array and qRT-PCR, respectively. The putative role of pre-selected candidate miRNAs in breast cancer development was assessed through exogenous expression of these miRNAs and their anti-miRNAs (‘antagomirs’) in MDA-MB-231 and MCF7 breast cancer-derived cells. Results Based on our expression profiling results, four candidate miRNAs (miR-10b, miR-26a, miR-146a and miR-153) were selected as being potentially involved in triple-negative breast cancer development. Exogenous expression assays revealed that miR-10b and miR-26a, but not miR-146a, can down-regulate the expression of BRCA1 in both triple-negative MDA-MB-231 and luminal epithelial MCF7 breast cancer-derived cells, whereas miR-153 could down-regulate BRCA1 expression only in MCF7 cells. In silico analysis of The Cancer Genome Atlas (TCGA) data confirmed that miR-146a is significantly higher expressed in triple-negative breast tumors compared to other (non triple-negative) breast tumors. Conclusion Our work provides evidence for the involvement of specific miRNAs in triple-negative breast cancer development through regulating BRCA1 expression. Electronic supplementary material The online version of this article (doi:10.1007/s13402-015-0239-3) contains supplementary material, which is available to authorized users.


Introduction
Breast cancer is the most commonly diagnosed cancer in women (23 % of all cancers), with an estimated diagnosis of 1.38 million new cases each year worldwide [1][2][3][4]. At the molecular level, breast cancers can be classified in four main subtypes, depending on the expression of hormone and growth factor receptors. The luminal subtypes A and B express both the estrogen receptor (ER) and the progesteron receptor (PR), the HER2-enriched subtype shows amplification and over-expression of the human epidermal growth factor receptor 2 (HER2), and the triple-negative subtype does Electronic supplementary material The online version of this article (doi:10.1007/s13402-015-0239-3) contains supplementary material, which is available to authorized users. not express any of these receptors. This classification correlates well with the prognosis, ranging from good for luminal A tumors to poor for the triple-negative tumors [5,6]. The latter category, representing~10 to 15 % of all breast cancers, is generally considered aggressive [7,8].
Triple-negative breast cancers are often linked to inactivation of the BRCA1 tumor suppressor gene. In familial breast cancer, 85 % of the BRCA1 mutated tumors are triple-negative [9,10]. In sporadic triplenegative breast cancers, BRCA1 is frequently inactivated at the transcriptional level, and it has been shown that BRCA1 inactivation may be due to methylation of its promoter brought about by ID4 [11,12]. More recently it has been found that BRCA1 may also be regulated at the post-transcriptional level by microRNAs [13].
MicroRNAs (miRNAs) are small non-coding 19-25 nucleotide-long RNAs that can post-transcriptionally regulate gene expression by binding to the 3'untranslated regions (3′-UTR) of target messenger RNAs (mRNAs), thereby leading to mRNA degradation or translational repression [14][15][16]. miRNAs have been shown to be involved in diverse biological processes [17][18][19]. In human breast cancer it has been shown that they can act either as tumor suppressors (i.e., miR-206, miR-17-5p, miR-125a, miR-125b, miR-200, let-7, miR-34 and miR-31) or as oncogenes (i.e., miR-21, miR-1 55 , m i R -10 b , m i R -37 3 a n d m i R -5 2 0c ) [ 2 0] . Previously, it has been shown that BRCA1 gene expression is relatively low in ER-negative and high-grade breast cancers [21] and that BRCA1 gene expression is significantly down-regulated in triple-negative breast cancers [22]. The number of miRNAs that may regulate the expression of BRCA1 or that may serve as transcriptional targets of BRCA1 is rapidly increasing. Considering the tumor-suppressive role of BRCA1, any perturbation in this regulatory function is likely to have an effect on BRCA1-mediated tumor development [23]. Moskwa et al. reported that miR-182-mediated downregulation of BRCA1 can impede DNA repair and, as such, affect breast cancer therapy [24]. It has also been shown that miR-146a and miR-146b-5p can down-regulate the expression of the BRCA1 gene in triple-negative sporadic breast cancers [18]. In contrast, a recent study reported that miR-146a expression levels increased simultaneously with BRCA1 expression levels. In addition, it was suggested that post-transcriptional regulation of epidermal growth factor receptor (EGFR) expression by BRCA1 may be mediated by miR-146a [25].
Here, we investigated miRNA expression profiles in nine human breast cancer-derived cell lines and one benign breast epithelium-derived cell line. By doing so, we identified a set of miRNAs putatively involved in BRCA1 gene expression regulation and breast cancer development.

Cell culture
Nine selected human breast cancer-derived cell lines and one benign breast epithelium-derived cell line (Table 1) were cultured in a humified incubator at 37°C containing 5 % CO2. Cell line MCF10a was obtained from the American Type Culture Collection (ATCC; Manassas, VA, USA) and was cultured in DMEM/F12 medium (Invitrogen Life Technologies, Carlsbad, CA, USA) supplemented with 10 % horse serum, 20 ng/ml EGF, 100 ng/ml cholera toxin, 500 ng/ml hydrocortisone, 2 mM L-glutamine and 20 ng/ml gentamicin.
Cell lines MCF7, T47D, MDA-MB-231, HCC1937 and MDA-MB-436 were also obtained from the ATCC and were cultured in RPMI-1640 medium (Invitrogen Life Technologies, Carlsbad, CA, USA) supplemented with 10 % fetal bovine serum (FBS), 2 mM L-glutamine and 20 ng/ml gentamicin. Cell lines SUM149PT and SUM1315MO2 were obtained from Asterand (Royston, Hertfordshire, UK) and were grown in Ham's F12 medium according to the supplier's instructions. Cell lines SUM1315-LXSN and SUM1315-BRCA1 were derived from the BRCA1 mutated (185delAG) SUM1315MO2 cell line, stably transfected with empty LXSN or LXSN-BRCA1 plasmids, respectively. The BRCA1 mutation status of the other cell lines has previously been reported by Elstrodt et al. [26].

Immunocytochemistry
Cells were fixed in Preservcyt solution (Thinprep) and cytoblocks were prepared using a Shandon Cytoblock kit (Thermo Scientific). The estrogen receptor (ER), progesterone receptor (PR), human epidermal growth factor receptor 2 (HER2), cytokeratin5/6, cytokeratin14, epidermal growth factor receptor (EGFR) and tyrosine kinase receptor cKIT protein status was determined by immunocytochemistry on 3 μmthick sections. Immunostainings were performed in a BenchmarkXT fully automatized stainer (Ventana Medical Systems) and scored semi-quantitatively by an expert pathologist under a light microscope.

MiRNA target prediction algorithms
Five available algorithms predicting miRNA-mRNA interactions were used to select miRNAs targeting the 3′-UTR of the BRCA1 gene transcript: TargescanHuman  For all algorithms we used the default parameters and we selected miRNAs that putatively target BRCA1 by at least three algorithms.

RNA extraction, reverse transcription and quantitative real-time PCR
Total cell RNA was extracted using QIAzol reagent (Qiagen miRNeasy mini kit) according to the manufacturer's instructions. The RNA quality and yields were assessed using an Agilent Bioanalyzer 2100 (Kit Agilent RNA 6000 nano) in conjunction with a spectrophotometer. One microgram of total RNA was reverse-transcribed using a miscript II RT kit (Qiagen, France). 15 nanograms of cDNA was used in triplicate for quantitative real-time PCR using a miscript SYBER Green PCR kit (Qiagen, France) on an ABI7900HT system and an Applied Biosystems ViiA™ 7 Real-Time PCR system. Selected miRNA primers were obtained from miScript Primer Assay (Qiagen, France) and are listed in supplementary  Table S1. U6 (Hs-RNU6-2-1 miscript Primer Assay) and 18S rRNA were used as internal controls. All samples were normalized to the internal controls and fold changes were calculated through relative quantification (RQ = 2 -ΔΔCT ).

PCR array-based profiling
The Human Breast Cancer miScript miRNA PCR Array (MTHS-109ZE-4) was used according to the manufacturer's protocol (Qiagen, France) for the profiling of 84 miRNAs known or predicted to alter in expression during breast cancer initiation and/or progression. The protocol enables real-time PCR profiling of mature miRNA on an ABI7900HT system in combination with a miScript SYBR Green PCR kit, which contains a miScript Universal Primer (reverse primer), and a QuantiTect SYBR Green PCR Master Mix. 250 nanogram of total RNA was reverse-transcribed using a miscript II RT kit (Qiagen).

Transfection assays
Selected miRNA mimics and miRNA inhibitors ('antagomirs') were obtained from miScript miRNA Mimic and miScript miRNA Inhibitor, respectively (Qiagen, France) and are listed in supplementary Table S2. The siRNA used for inhibiting BRCA1 gene expression (siBRCA1) was synthesized by Thermo Scientific, Dharmacon (M-003,461-02-0005, siGENOME SMART pool, Human BRCA1 (672), 5 nmol). Cells were seeded at a density of 300,000 cells per well in 6-well plates. Twenty-four hours later, miRNAs and siRNAs were transfected using Lipofectamine 2000 reagent (Invitrogen, CA, USA) following the manufacturer's instructions. As a mock control, cells were transfected with transfection reagent alone (Tmock). Three days after transfection, the cells were washed with PBS and recovered for further analysis.

2.8
In silico analysis of miRNAs using TCGA Both clinical and miRNA sequencing data of invasive breast cancers were downloaded from The Cancer Genome Atlas (TCGA) database. A total of 519 patients with information on ER, PR and HER2 status were selected to compare the expression profiles of four selected miRNAs (miR-10b, miR-26a, miR-146a and miR-153) in the respective tumors. 88 cases were found to have a negative ER, PR and HER2 (i.e., triple-negative) phenotype, whereas 431 cases were positive for at least one of these receptors.

Statistical analysis
Student's t-test was used to assess statistical differences in mean expression between groups. A p-value ≤ 0.05 was considered statistically significant.

Results
3.1 miR-10b and miR-26a are preferentially expressed in triple-negative breast cancer-derived cell lines In order to identify miRNAs implicated in triple-negative breast cancer, we used a miScript miRNA PCR Array to profile the expression of 84 miRNAs considered to be involved in breast tumorigenesis. Ten cell lines were analyzed: one benign breast epithelium-derived cell line (MCF10a) The cell lines were characterized by immunocytochemistry as presented in Table 1. In order to reveal whether particular miRNAs are specifically involved in triple-negative breast cancers, we compared their expression in the respective cell line subgroups. For this analysis we focused on six miRNAs that were found to be differently expressed in at least two subgroups by PCR Array, i.e., miR-10b, miR-15b, miR-26a, miR-155, miR-206 and miR-485-5p ( Table 2, Table S3). To validate the PCR Array results, we performed qRT-PCR analysis of the six selected miRNAs in the 10 cell lines. Only two of these miRNAs (i.e., miR-10b and miR-26a) were found to be well-expressed in these cell lines. The other miRNAs were expressed at low levels or were not detectable at all. miR-10b was highly expressed in 3 triple-negative breast cancer cell lines (MDA-MB-231, MDA-MB-436 and SUM1315MO2) (Fig. 1a). miR-26a was found to be expressed in all 10 cell lines, and to be higher expressed in the triple-negative cell lines than the luminal cell lines (MCF7, T47D) (Fig. 1b).
3.2 miR-146a and miR-153 target the 3′-UTR of BRCA1 BRCA1 gene expression is frequently down-regulated in triple-negative breast cancers. We hypothesized that the identification of miRNAs targeting the 3′-UTR of BRCA1 might be used to uncover miRNAs involved in the development of triple-negative breast cancers. To this end, we performed an in silico computational search using five available algorithms: TargescanHuman [27], Tarbase [28], Miranda [29], MicroCosmTargets [30] and PicTar [31]. For all algorithms we used default parameters, and we selected miRNAs predicted to target BRCA1 by at least three algorithms. From the selected miRNAs (listed in Table S4), we chose five miRNAs, i.e., two miRNAs (miR-146a and miR-146b-5p) that were previously found to be overexpressed in breast tumors and to down-regulate BRCA1 expression [18] and three miRNAs (miR-132, miR-212, miR-153) that were most frequently predicted to interact with the 3′-UTR of BRCA1. We validated the expression of these five selected miRNAs by qRT-PCR in the 10 cell lines and found that miR-146a was higher expressed in 3 triple-negative breast cancer cell lines: MDA-MB-231, MDA-MB-436 and SUM1315-BRCA1 (Fig. 1c), that miR-153 was higher expressed in 4 triple-negative breast cancer cell lines: MDA-MB-436, SUM1315MO2, SUM149PT and HCC1937 (Fig. 1d), and that miR-146b-5p, miR-132 and miR-212 were expressed in all 10 cell lines at different levels ( Fig. S1) without a significant difference between the triple-negative, luminal or benign cell lines. Based on these results we focused our further analyses, next to miR-10b and miR-26a, on miR-146a and miR-153.

BRCA1 expression is regulated by miR-10b, miR-26a and miR-153
The above expression profiling and in silico analyses allowed us to select 4 miRNAs, miR-10b, miR-26a, miR-146a and miR-153 for further analysis. To determine the effect of these miRNAs on endogenous BRCA1 expression, we transfected two BRCA1 wild-type cell lines (MDA-MB-231 and MCF7) with the respective miRNAs (miRNA mimics) and anti-miRNAs (inhibitors or 'antagomirs'). The efficiency of the miRNA mimics and inhibitors was verified by miRNA quantification using qRT-PCR (Fig. S2), revealing sufficient efficiencies only for miR-10b, miR-26a, miR-146a, anti-miR-146a and miR-153. The BRCA1 siRNA efficiency was also tested and found to be sufficient (Fig. S3). Since MCF7 is a luminal cell line and MDA-MB-231 a triple-negative cell line, our data can potentially reveal whether the miRNAs tested may have a different effect on BRCA1 in these two subtypes of breast cancer.
To determine whether miR-146a affects endogenous BRCA1 expression, we compared the expression level in MDA-MB-231 and MCF7 cells after transfection with miR-146a, anti-miR-146a and a mock control (transfection reagent alone). By doing so, we found that miR-146a does not appear to regulate the expression of BRCA1 (Fig. 2). Next, we tested the effect of miR-153 on the expression of BRCA1. The results obtained (Fig. 2) indicate that miR-153 can down-regulate the expression of BRCA1 in MCF7 cells. In contrast, we found that in MDA-MB-231 cells miR-153 induced BRCA1 up-regulation. After miR-10b transfection into MDA-MB-231 and MCF7, the BRCA1expression level was found to be decreased in both of them (Fig. 2). These results suggest that miR-10b can down-regulate the expression of BRCA1 in MDA-MB-231 and MCF7 cells. After miR-26a transfection into MDA-MB-231 and MCF7 cells, the expression level of BRCA1 was again found to be decreased in both of them (Fig. 2), indicating that also miR-26a can down-regulate the expression of BRCA1 in MDA-MB-231 and MCF7 cells.

EGFR expression is not regulated by miR-146a in MDA-MB-231 and MCF7 cells
Kumaraswamy et al. reported that BRCA1 can induce miR-146a expression which, in turn, can repress EGFR expression in mammary cells [25]. To determine whether miR-146a can affect endogenous EGFR expression in our cell lines, we compared its expression level in MDA-MB-231 and MCF7 cells after transfection with miR-146a, anti-miR-146a, Tmock (transfection reagent alone) and siBRCA1 (Fig. S4). After doing so, we found that in these cell lines miR-146a had no effect on the expression of EGFR.

Stimulation and inhibition of cell proliferation by miRNAs
To assess the effect of miR-10b, miR-26a, miR146a, anti-miR-146a and miR-153 on MDA-MB-231 and MCF7 cell proliferation, the cells were transfected with the respective (anti-) miRNAs, as also siBRCA1 and Tmock (transfection reagent alone). In MDA-MB-231 cells, the proliferation assay (see materials and methods) showed that miR-26a significantly stimulated and that miR-146a inhibited proliferation compared to mock-transfected cells. In MCF7 Bold entries are miRNAs that were found to be differently expressed in at least two subgroups by PCR Array miR-10b, miR-26a, miR-146a, miR-153 as potentiel TNBC biomarkers cells, miR-153, miR-10b and miR-26a significantly inhibited proliferation compared to mock-transfected cells (Fig. 3).
3.6 Differential expression of miR-10b, miR-26a, miR-146a and miR-153 in TCGA TCGA was queried for the expression of miR-10b, miR-26a, miR-146a and miR-153 in different primary breast cancer tissues. As shown in Fig. 4 the average expression of miR-153 and miR-10b was found to be significantly lower in triplenegative tumors (p = 0.00045 and p = 0.00038, respectively) than in luminal tumors, while the expression of miR-146a was found to be significantly higher in triple-negative tumors (p = 0.000005). For miR-26a, the difference in average expression was not found to be significant in these tumors (p = 0.6234).

Discussion
Triple-negative breast cancers, representing~10 to 15 % of all human breast cancers, have gained growing interest in recent years [32]. In addition, the role of microRNAs (miRNAs) in the epigenetic regulation of many cellular processes has increasingly become recognized as an important way to finetune gene expression [33]. Here, expression profiling, transfection and in silico assays were performed to identify potential triple-negative breast cancer miRNA biomarkers, i.e., miR-10b, miR-26a, miR-146a and miR-153. We found that miR-10b was higher expressed in triplenegative breast cancer-derived cell lines compared to luminal breast cancer-derived cell lines. TCGA data, however, indicated that the average expression of miR-10b was significantly lower in triple-negative tumors than in non triple-negative tumors. O'Day and Lal reported that miR-10b was over-expressed only in metastatic cancer cells, and was found to promote tumor cell migration, invasion and metastasis in vivo [20]. Whereas miR-10b has been found to be highly expressed in several human cancers [34], Biagioni et al. found that miRNA-10b was downregulated in tumor tissues compared to their normal matched counterparts due to hypermethylation of CpG islands upstream from the miR-10b/10b* locus [35]. We found that miR-10b can down-regulate the expression of BRCA1 in MDA-MB-231 and MCF7 cells. High miR-10b expression may thus be a way to silence BRCA1 expression in BRCA1 wild-type cells. Recently, it was reported that miR-10b may also play a critical role in TGF-induced epithelial-mesenchymal transition (EMT) in breast cancer and, as such, may be considered as a possible therapeutic target [36]. Moreover, miR-10b was shown to be able to promote the invasion and metastasis of tumor cells through post-transcriptional regulation of HOXD10 [37]. Overall, these results indicate that miR-10b may act as an oncogenic miRNA in breast cancer cells. Nevertheless, Biagioni et al. reported that exogenous miR-10b expression reduced both  In silico expression analysis of miRNAs using TCGA data. Clinical and miRNA expression data for breast cancer were downloaded from The Cancer Genome Atlas (TCGA) database. The expression of four miRNAs (miR-10b, miR-26a, miR-146a and miR-153) was compared in 88 breast tumors with a negative ER, PR and HER2 status (i.e., triple negative phenotype) and in 431 breast tumors that were positive for at least one of the receptors. Student's t-test was used to assess statistical differences in mean expression levels between these two groups the in vitro and in vivo proliferative capacities of MCF7 and MDA-MB-468 cells [35]. Our results are compatible with those of Biagioni et al., as we found that miR-10b significantly inhibited the proliferation of MCF7 cells (p = 0.006). Thus, restoration of miR-10b expression may hold a therapeutic promise for breast cancer treatment.
Another potential biomarker of triple-negative breast cancer that we identified is miR-26a, and we found that miR-26a can down-regulate the expression of BRCA1 in MDA-MB-231 and MCF7 cells. These results suggest an oncogenic function of this miRNA in breast cancer. In contrast, Gao et al. reported that miR-26a could inhibit the proliferation and migration of breast cancer cells through repression of MCL-1, and that miR-26a could increase the sensitivity of breast cancer cells to paclitaxel [38]. Others found that miR-26a can promote ovarian cancer cell proliferation by targeting ERα [39]. We found that miR-26a significantly stimulated the proliferation of MDA-MB-231 cells and inhibited the proliferation of MCF7 cells. We also found that the average expression of miR-26a was not significantly different between triple-negative and non triple-negative breast tumors in the TCGA database. The role of miR-26a thus appears to be complex and may dependent on the tissue and/or tumor context.
Our cell line and TCGA data analyses show that miR-146a is significantly over-expressed in triple-negative breast cancers. Similar results were reported by others [18,25]. We also found that miR-146a does not affect the expression of BRCA1 in MDA-MB-231 and MCF7 cells. In contrast, Gracia et al. reported that miR-146a can down-regulate the expression of BRCA1 in triple-negative sporadic breast cancers [18]. This discrepancy may be due to the different cell lines used for the exogenous expression assays, i.e., Gracia et al. used HeLa cells and three mammary cell lines (MDA-MB-468, MDA-MB-157 and MDA-MB-436). Another study reported that miR-146a can inhibit the invasion and migration of the highly metastatic human breast cancer cell line MDA-MB-231 [40]. These data are in line with our results showing that miR-146a can inhibit the proliferation of MDA-MB-231 cells. In a recent study it was hypothesized that BRCA1 may down-regulate the expression of the EGFR gene by increasing the miR-146a level in breast cancer-derived cells (HCC1937) [25]. Our results are compatible with this hypothesis, as we found an increase in miR-146a expression in BRCA1 wild-type SUM1315-BRCA1 cells compared to BRCA1 mutated SUM1315-LXSN cells. We also found that exogenous expression of miR-146a decreased the proliferation of MDA-MB-231 cells that highly express EGFR, but not in MCF7 cells that do not express EGFR. qRT-PCR assessment of EGFR expression in exogenous miR-146a expressing MDA-MB-231 and MCF7 cells revealed that the EGFR gene was not significantly regulated by this miRNA. Taken together, these results suggest a complex feedback regulation between miR-146a and BRCA1, which again most likely depends on the cellular context.
Using bioinformatic tools, we identified a binding site for miR-153 in the 3′-UTR of BRCA1. Subsequently, we found that miR-153 can induce BRCA1 down-regulation in MCF7 cells and up-regulation in MDA-MB-231 cells. These contradictory results reflect those reported in the literature. On one hand, Anaya-Ruiz et al. reported that silencing of miR-153 significantly inhibited the growth, reduced the proliferation and induced apoptosis in the triple-negative sporadic breast cancer-derived cell line MDA-MB-231. These results indicate that miR-153 may function as an oncogenic miRNA, whose deregulation could be involved in the initiation and/or development of human breast cancer [41]. In another report, Wu et al. suggested that miR-153 may play an important role in promoting the proliferation of human prostate cancer cells and, as such, may represent a novel mechanism of miRNAmediated PTEN silencing in prostate cancer cells [42]. These data are compatible with our finding that BRCA1 upregulation by miR-153 had no effect on the proliferation of MDA-MB-231 cells. On the other hand, Zhao et al. found that transient transfection of miR-153 into glioblastoma multiforme stem cells (GBM-SCs) can inhibit their stemness properties, repress their growth potential and induce apoptosis [43]. These latter data are consistent with our results in MCF7 cells: although BRCA1 was down-regulated by miR-153, it inhibit proliferation. This difference in response observed between MDA-MB-231 and MCF7 cells may be explained by differences in the endogenous levels of miR-153, i.e., higher in luminal MCF7 cells and lower in triple-negative MDA-MB-231 cells. TCGA data analysis also indicated that the average expression of miR-153 is significantly lower in triple-negative breast cancers. We conclude that also the function of miR-153 is complex, and identifying targets of miR-153 would be an interesting means to better understand its role in cell proliferation and carcinogenesis-related processes.
In conclusion, we identified miR-10b, miR-26a, miR-146a and miR-153 as potential triple-negative breast cancer biomarkers. These miRNAs may be instrumental for the future design of novel targeted treatment options of these breast cancers. This 'from bench-top to bed-side' translation remains, however, challenging [44] and additional research is required to better understand the function(s) of these miRNAs.
Acknowledgments This research was funded by grants from the ministry of Higher Education and Scientific Research in Tunisia, the Laboratory of Biochemistry Research unit UR 12ES08 Cell Signaling and Disease, Faculty of Medicine of Monastir, the Molecular Oncology Laboratory at Center Jean Perrin in Clermont-Ferrand and the Fuda (Founding of the University of Auvergne).
Compliance with ethical standards We wish to confirm that there are no known conflicts of interest associated with this publication and that there has been no significant financial support for this work that could have influenced its outcome.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http:// creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.