Downregulation of 15-hydroxyprostaglandin dehydrogenase during acquired tamoxifen resistance and association with poor prognosis in ERα-positive breast cancer

Aim: Tamoxifen (TAM) resistance remains a clinical issue in breast cancer. The authors previously reported that 15-hydroxyprostaglandin dehydrogenase (HPGD) was significantly downregulated in tamoxifen-resistant (TAMr) breast cancer cell lines. Here, the authors investigated the relationship between HPGD expression, TAM resistance and prediction of outcome in breast cancer. Methods: HPGD overexpression and silencing studies were performed in isogenic TAMr and parental human breast cancer cell lines to establish the impact of HPGD expression on TAM resistance. HPGD expression and clinical outcome relationships were explored using immunohistochemistry and in silico analysis. Results: Restoration of HPGD expression and activity sensitised TAMr MCF-7 cells to TAM and 17β-oestradiol, whilst HPGD silencing in parental MCF-7 cells reduced TAM sensitivity. TAMr cells released more prostaglandin E2 (PGE2) than controls, which was reduced in TAMr cells stably transfected with HPGD. Exogenous PGE2 signalled through the EP4 receptor to reduce breast cancer cell sensitivity to TAM. Decreased HPGD expression was associated with decreased overall survival in ERα-positive breast cancer patients. Conclusions: HPGD downregulation in breast cancer is associated with reduced response to TAM therapy via PGE2-EP4 signalling and decreases patient survival. The data offer a potential target to develop combination therapies that may overcome acquired tamoxifen resistance.


Introduction
Due to its ability to bind to and modulate oestrogen receptor alpha (ERα) activity, Tamoxifen (TAM) was the first targeted therapy for breast cancer [1], with its widespread use in the clinic now extending to four decades.While aromatase inhibitors are now the preferred first line endocrine treatment in the post-menopausal setting [2], TAM is favoured in pre-menopausal breast cancer patients [3].However, approximately 70% of patients initially respond to TAM, but most eventually acquire resistance [2].Acquisition of TAM-resistance (TAMr) continues to be a major limitation for long-term management of breast cancer patients.However, the mechanism(s) responsible for TAM insensitivity is not yet fully understood.
In the last three decades, various groups have developed cell line models in order to help unravel the mechanistic basis of TAMr.Most of these have been achieved through long term culture of the hormonesensitive MCF-7 human breast cancer cells in sub-lethal doses of TAM, during which resistant sub-clones develop after 3-6 months [4][5][6][7][8].Using two independently-derived ER positive MCF-7 human breast cancer cell models of acquired TAMr developed in our laboratory [9], an Affymetrix ® microarray study showed that 131 genes were upregulated and 156 were downregulated by at least 3-fold, in TAMr MCF-7 cells relative to isogenic control cells [10].One of the genes that were downregulated was 15-hydroxyprostaglandin dehydrogenase (HPGD).HPGD (EC 1.1.1.141)is a member of the short-chain alcohol (OH) dehydrogenase family and is the key NAD + -dependent enzyme responsible for the biological inactivation of prostaglandins, including prostaglandin E 2 (PGE 2 ), which is synthesized via the cyclooxygenase (COX) pathway.Downregulation of HPGD has been shown in several malignancies including lung, colon, bladder, endometrial, and gastric cancer and has been shown to have a tumour suppressor roles in some settings [11][12][13][14][15][16].Studies in breast cancer are contradictory; higher HPGD expression has been reported in ERαpositive MCF-7 cell with reduced expression in ERα-negative MDA-MB-231 cells, where up-regulation was observed following treatment with demethylating agents [17].In clinically more aggressive primary breast cancers, particularly those with a triple negative phenotype, and in breast cancer metastases, HPGD is overexpressed, with high HPGD expression associated with poor prognosis and reduced survival [18].In the same study, HPGD silencing induced a mesenchymal-epithelial transition resulting in a less migratory phenotype.More recently, data mining from publicly available breast cancer gene expression datasets from the Cancer Genome Atlas (TCGA, https://www.cancer.gov/about-nci/organization/ccg/research/structuralgenomics/tcga)and Oncomine (https://www.oncomine.org/),revealed decreased expression of HPGD mRNA in breast cancer samples compared with normal, healthy breast tissue [19].Furthermore, HPGD expression was lower still in more aggressive basal/triple negative and HER2-positive breast cancers.Conversely, in rare apocrine breast carcinomas, HPGD was highly expressed [20].Application of a triclustering algorithm, δ-TRIMAX to microarray datasets characterized HPGD as one of ten so-called hub-genes associated with breast cancer [21], lending support for a role in breast carcinogenesis.
Building on our previous data [10], the aim of this study was to investigate a causal link between HPGD expression and function and TAMr, as well as explore the relationship between HPGD expression in breast cancer samples and clinical outcomes.The effects of HPGD expression on outcome were examined through in silico analysis and by immunohistochemistry in a breast cancer cohort with long term follow up, including breast cancer outcomes during adjuvant TAM.

Cell lines and culture
The two stable TAMr MCF-7 human breast cancer cell lines have been described previously [9,10].All MCF-7 isogenic lines were cultured in phenol-red-free RPMI 1640 containing L-glutamine (Invitrogen, UK) supplemented with 5% charcoal-stripped steroid-depleted foetal calf serum (FCS; Harlan SeraLab, UK), 100 U/mL penicillin and 100 U/mL streptomycin and 100 nM 4-hydroxytamoxifen (TAM; Sigma-Aldrich, UK) for 12-24 months.Parental cells [termed wild type (WT) MCF-7] were cultured in the same medium, but with 0.01% (v/v) ethanol vehicle.Experiments were conducted in phenol-red-free RPMI 1640 supplemented with 5% charcoal-stripped steroid-depleted FCS + TAM (0.1-100 nM) or 17β-oestradiol (E2; 0.001-1 nM).Bi-monthly mycoplasma checks were consistently negative and annual short tandem repeat (STR) profiling confirmed cell provenance, both carried out as a service in Leeds.To minimize genetic drift, cell stocks were frozen at low passage and experimental cultures replaced from these stocks every 3-6 months.We also studied LCC1, -2 and -9 breast carcinoma cells, which were derived originally in the Clarke laboratory at Georgetown, Washington DC and cultured as described [7,8,22].LCC1 cells are estrogen-independent and -responsive TAM-sensitive breast cancer cells, which are derived from an estrogen-independent variant of MCF-7 cells (MCF7/MIII) through in vivo selection in oophorectomized nude mice (with circulating estrogen levels similar to a postmenopausal woman), and subsequently re-cultured in vitro to become a stable cell line.LCC2 cells are stable, ER-positive, estrogen-independent, TAM-resistant, and respond to fulvestrant [7].They were derived from the selection of LCC1 for TAM resistance in vitro.LCC9 cells exhibits cross resistance to TAM and are unresponsive to fulvestrant [7].They were derived from the selection of LCC1 for fulvestrant resistance in vitro.MCF7:5C and MCF-7:WS8 are oestrogen-independent and dependent, respectively, ERpositive, progesterone receptor (PR)-negative, and TAM-resistant breast cancer cells [23].They were derived from WT MCF-7 cells following long-term estrogen deprivation (LTED).The main characteristics of the cell lines used are detailed in Table 1.

Microarray analysis of LCC1, -2 and -9
Cells were seeded at 300, 000 cells per well in 6-well plates (Corning, UK).They were treated the next day with 1 nM E2 for the indicated durations.All experiments were performed in triplicate.Subsequently, cells were harvested in TRIzol (Life Technologies, UK) and total RNA was isolated using the Qiagen RNeasy RNA purification kit (Qiagen, UK).Isolated RNA (1 µg) was then processed on the Affymetrix Clariom S microarray platform (Affymetrix, UK) and used to determine the significance of differential expression of HPGD at basal levels.

LCC2
Stepwise selection against TAM

PGE 2 ELISA and HPGD enzyme activity assay
Both assays were performed using previously protocols published [26].For PGE 2 , serum-free conditioned medium was collected from cells after overnight incubation, and PGE 2 levels were measured using a competitive immunoassay (Amersham Biosciences, UK), and normalised to total cellular protein content as measured by Bradford assay.HPGD activity was measured in lysates of WT MCF-7 and TAMr cells by assessing the transfer of tritium from 3 H-PGE 2 to glutamate, catalysed by HPGD and glutamate dehydrogenase [27].

EP4 and HPGD knockdown
Small-interfering (si)-RNAs targeting EP4 or HPGD (ThermoFisher Scientific, UK) were reverse transfected into cells using Lipofectamine 2000 according to the manufacturer's instructions.Briefly, MCF-7 cells were mixed with media containing small-interfering RNA (siRNA)/lipid (final siRNA concentration 10 nM) and seeded into 96 well plates at 0.4 x 10 4 cells/well and allowed to establish overnight before incubation with drug or vehicle control.Response to treatment was assessed by MTT assay.EP4 or HPGD knockdown was assessed by real-time quantitative polymerase chain reaction (qRT-PCR) as described previously [28].

In silico transcriptomic analysis
The relationship between HPGD expression and cancer outcomes was analyzed using Kaplan-Meier Plotter (KMplot, https://kmplot.com/analysis/)[30].Using the 2017 release of the database, a cohort of ERαpositive breast tumours, previously treated with TAM only and no adjuvant chemotherapy, was analyzed and the association between HPGD expression and overall survival (OS).Patients were dichotomized as high or low HPGD expression using lower tertile as a cut-off.Multivariate analysis was done using in-built software with Ki67, ER and HER2 as covariates.

In silico analyses for the identification of miRNAs targeting HPGD
We applied miRWalk version 3.0 [31] and TargetScan version 7.2 [32] to a previously described Affymetrix miRNA dataset of endocrine sensitive MCF-7:WS8 and resistant MCF-7:5C cells [33] in order to identify miRNA-target interactions sites for major HPGD transcripts using default settings.OncoLnc (http://www.oncolnc.org/) was used to retrieve HPGD and miR-3200-3p expression data from matching breast cancer cases in the TCGA database [34].Kaplan-Meier Plotter [30] was used to investigate the significance of miR-3200-3p on survival in ERα+ breast cancer treated with any endocrine therapy (not limited to TAM).Cases were dichotomised as high or low miR-3200-3p expression using the median expression value as a cut-off.

Statistical analysis
For in vitro analyses, one-way ANOVA was performed.The log rank test was used to compare patient survival in the primary breast cancer cohort.Analyses were performed using GraphPad Prism version 7.03 (GraphPad Software, La Jolla California, USA).For Affymetrix analyses, ANOVA P-values were calculated using the Affymetrix Transcriptome Analysis Console (Affymetrix).The time course plots were plotted as fold-change with the baseline value set at 1.

HPGD is downregulated in TAMr MCF-7 MMU2 cells
We demonstrated that HPGD was downregulated in two independent TAMr MCF-7 cell lines, which had been cultured continuously in 100 nM TAM for 12 (TAMr 1) and 21 (TAMr 2) months [10], respectively by Western blot (Figure 1a) and real time PCR (Figure 1b) compared with parental MCF-7 cells.Microarray analysis of the independent LCC cell line series and MCF7:5C cells showed that only the oestrogen-independent but oestrogen-responsive TAMs cell line LCC1 showed increased expression of HPGD in response to E2 over time (Figure 1c).

Effect of HPGD overexpression in MCF-7 TAMr 2 cells
We then explored the effects of introduction of HPGD into MCF-7 TAMr 2 cells.Successful expression of HPGD protein following stable transfection was confirmed by Western blotting (Figure 2a).This restored activity of HPGD in each of the transfected cell lines to levels like WT MCF-7 (Figure 2a and b).This also sensitized these cells to the inhibitory effects of TAM to approximately 60% of the parental MCF-7 cell response (Figure 2c), whilst restoring sensitivity to E2 almost completely (Figure 2d).

Modulating the PGE 2 axis influences cellular response to TAM
Since HPGD is responsible for the biological inactivation of prostaglandins including PGE 2 , we tested whether increased PGE 2 signaling might explain TAMr.We first examined PGE 2 levels in serum-free conditioned medium collected from MCF-7 and MCF-7 TAMr 2 cells.TAMr cells released higher amounts of PGE 2 compared with parental MCF-7 cells and this could be reduced to levels approaching that of MCF-7 by HPGD overexpression in these cells (Figure 3a).Next, we tested the effect of adding exogenous PGE 2 (1 µM) on TAM response in two different ERα-positive breast cancer cell lines, MCF-7 and T47D.PGE 2 reduced sensitivity of both cell lines to TAM (Figure 3b).As PGE 2 -signalling is mediated by G-protein coupled receptors, we examined if the stimulatory EP2 or EP4 receptors might be responsible for mediating the effects of PGE 2 on TAM sensitivity.We demonstrated that PGE 1 -OH (an EP4 agonist) was able to mimic the effects PGE 2 in decreasing TAM sensitivity (Figure 3c; P < 0.015), albeit to a lesser extent than PGE 2 .However, butaprost (EP2 agonist) had no effect at equivalent concentrations (Figure 3d).This suggested that the EP4 receptor might be responsible for mediating the effect of PGE 2 on TAM sensitivity.Partial silencing of EP4 in WT MCF-7 cells was achieved using one siRNA (Figure 3e).EP4 silencing using siRNA 1 inhibited PGE 2 induced TAM resistance in these cells (Figure 3f).

Loss of HPGD expression predicts worse outcomes in breast cancer
We analysed the relationship between HPGD expression and clinical outcome for breast cancer patients using the on-line resource Kmplot [30] and from our own cohort of breast cancer patient's data in separate survival analyses.In the Kmplot cohort, we found that ER+ patients who received adjuvant TAM only and had 'low' expression of HPGD, had reduced OS by univariate analysis [hazard ratio (HR) = 0.28 (0.13-0.59),P = 0.0007] (Figure 4a).This remained significant on multivariate analysis (P = 0.001, Figure 4b).We then conducted a retrospective immunohistochemical study of HPGD expression in our own cohort of breast patients, who all received adjuvant TAM therapy and related HPGD expression to survival.As these TMAs have been used extensively in several other studies [24,29,35], the number of viable cores that could be reliably evaluated was only 144 out of 350 cases.The level of HPGD staining was generally weak to moderate (Figure 4c) in positive cases, however cases which completely lacked HPGD were associated with worse OS [HR = 0.3 (0.15-0.69),P = 0.047] (Figure 4d), in line with the in silico analysis (Figure 4a).The small number of TMA cases available precluded multivariate analysis.

In silico identification of miRNA-HPGD target interactions
With growing evidence that HPGD is inactivated epigenetically, we sought to identify miRNAs that may post-transcriptionally modulate HPGD expression.We interrogated a global miRNA data set established in an independent TAM resistance model MCF-7:5C [33].These TAMr cells were derived from WT MCF-7 cells (Table 1).In silico analysis revealed numerous miRNAs that were differentially expressed in TAMr MCF-7:5C cells compared with TAMs MCF-7:WS8 cells and that are predicted to interact with the three prime untranslated region (3'-UTR) of the HPGD transcripts (Table 2).We then tested the association between each of the listed miRNAs and HPGD mRNA levels in 987 breast cancer samples from the TCGA database using OncoLnc [34] where both miRNA and mRNA expression data were available.Only one miRNA, namely miR-3200-3p, had an inverse relationship between its expression levels and HPGD mRNA levels in this dataset (Spearman, ρ = -0.25,P < 0.0001).We used the Kmplot online resource [30] to analyse the association between miR-3200-3p expression and OS in ER+ breast cancer patients treated with any endocrine therapy, irrespective of grade, molecular subtype, and lymph node status.High miR-3200-3p expression was associated with significantly reduced OS (Figure 5a) and a weak negative correlation (R 2 = 0.25) was observed between miR-3200-3p and HPGD (Figure 5b) using data from TCGA analysed using the OncoLnc platform [34].

Discussion
We have demonstrated a downregulation of HPGD expression in TAM-resistant breast cancer both in vitro in multiple series of TAMs and TAMr isogenic cell line pairs and in clinical samples using retrospective immunohistochemistry studies and in silico approaches.We have also shown that TAM resistance is associated with PGE 2 signalling via the EP4 receptor in vitro.A key role of HPGD is to inactivate PGE 2 .Having confirmed its downregulation in TAMr cells using various experimental approaches and in different cell lineages, we showed that stable transfection of human HPGD into MCF-7 TAMr 2MCF-7 cells restored HPGD enzyme activity, at the same time significantly sensitising these cells to TAM, and suggesting a functional role of HPGD downregulation in acquired TAMr.Our hypothesis that HPGD mediated its effects via this axis was strengthened by the observation that treatment of 2 different ER+ TAM-sensitive cell lines with exogenous PGE 2 increased the resistance of both cell lines to the growthinhibitory effects of TAM.PGE 2 exerts its effects through binding to cell surface G-protein coupled EP receptors.Four subtypes have been identified, EP1, -2, -3 and -4, with each coupled to different intracellular signalling pathways.All are stimulatory except for EP3, which exerts inhibitory effects.We used a pharmacological approach to target EP2 and EP4, treating parental MCF-7 and T47D cells with PGE 1 -OH (EP4 receptor agonist) or with butaprost (EP2 receptor agonist).PGE 1 -OH was able to mimic the effects of PGE 2 (albeit to a lesser extent), whereas butaprost had no effect at equivalent concentrations.We suggest that the EP4 receptor is responsible for mediating the TAMr effect of PGE 2 .While our knockdown experiment implied a role for EP4, we were only able to achieve partial knockdown with a single siRNA, which limits the conclusions from our data.Consistent with our data, a study has shown that EP4 gene expression was upregulated in TAMr MCF-7 cells, also generated from long term oestrogen deprivation, and in patients resistant to aromatase inhibitor therapy.Furthermore, the exposure to two separate EP4-specific antagonists, GW627368X and ONO-AE3-208, decreased estrogen-independent cell growth in LTED MCF-7 cells [36].These authors proposed a molecular mechanism by which EP4 signalling may drive ligand-independent ERα activation by driving the binding of the coactivator-associated arginine methyltransferase 1 (CARM1), an ERα co-factor, to the receptor, leading to its activated status [36].
There is evidence that HPGD is post-transcriptionally regulated via miRNAs, with miR-21 identified as targeting HPGD in oral tongue squamous cell carcinoma [37].Silencing of miR-21 expression in human breast cancer cell lines has been shown to increase TAM sensitivity [38].miR-21 was found to be overexpressed in TAMr cell line compared to TAMs cell lines [39], which may contribute to HPGD downregulation in TAM resistant breast cancer.However, miR-21 was not identified as differentially expressed in TAMr MCF-7:5C cells compared with TAMs MCF-7:WS8 cells.We identified several differentially expressed miRNAs predicted to interact with the 3'UTR of HPGD transcripts.Of these, miR-3200-3p expression was found to be inversely correlated with HPGD mRNA level, suggesting a possible role of miR-3200-3p in the regulation of HPGD expression.Using an in silico approach we found that increased miR-3200-3p levels were significantly associated with decreased OS in ER positive breast cancer patients treated with endocrine therapy.There are no data available thus far establishing a direct interaction between miR-3200-3p and HPGD mRNA, however correlation between HPGD and miR-3200-3p in breast cancer cases in the TCGA database using the OncoLnc platform [34] revealed a weak negative correlation.A key next step would be to validate that the 3'-UTR of HPGD is a target of miR-3200-3p in a luciferase reporter assay in vitro, with in vivo validation also warranted.Interestingly, while this work was under review, miR-3200-3p was identified using data from TCGA as a so-called miRNA master regulator (MMR) in breast cancer and was one of 61 MMRs which displayed high oncogenic activity [40].
The immunohistochemistry on human breast cancer samples and the in silico analysis demonstrated that reduced expression of HPGD was associated with worse patient outcomes, consistent with a study demonstrating that low levels of HPGD predicted reduced recurrence free survival and OS in a cohort of breast cancer patients, split into normal weight, overweight and obese groups, in which few received endocrine therapy (15%) [41].In our study, all cases received adjuvant TAM, with HPGD downregulation associated with reduced response to endocrine treatment and worse clinical outcomes.Some studies of HPGD expression in breast tissues have demonstrated higher expression in normal breast tissue than breast cancer [19], supporting a tumour suppressive role [17,42].However, this is an inconsistent finding.One study showed that high, rather than low, HPGD expression was significantly associated with poor outcome in triple negative breast cancer [18].However, tumours lacking expression of either ERα, PR and HER2 are recognised to have worse prognosis.Our cohort was predominantly ER-positive (77%) and the in silico analysis was restricted to ER-positive cases only.Such context-dependency has been supported through a data mining approach where HPGD expression was found to be differentially expressed across the various molecular subtypes of breast cancer [19].However, in this study [19], HPGD expression was significantly decreased in triple negative breast cancer compared to other molecular subtypes.
Data extrapolated from clinical trials is beginning to indicate that the risk of breast cancer relapse is reduced by extending time on endocrine therapy from 5 to 10 years [43].While this is good for patient outcome, there will likely be a lag period before these benefits are seen in the clinic.Understanding the mechanisms which underline endocrine therapy resistance remains an important biological question.This study demonstrated that HPGD downregulation in ER-positive breast cancer is associated with reduced response to adjuvant TAM therapy via PGE 2 -EP4 signalling and decreased patient survival.These data offer a potential target to develop combination therapies that may overcome TAM acquired resistance.

Table 1 .
Characteristics of TAM-resistant breast cancer cell lines derived from MCF-7

Figure 1 .Figure 2 .
Figure 1.HPGD is downregulated in TAMr derivatives of MCF-7 cells.Western blot (a) showing HPGD protein expression in WT MCF-7 and its loss in two, independently derived, TAMr MCF-7 cell lines, TAMr 1 and TAMr 2; qRT-PCR (b) shows relative expression of HPGD mRNA in these cell lines.Expression was determined using RPLP0 as the reference gene and further normalized to baseline at time 0 h (± SD); in (c), expression of HPGD mRNA in LCC1 cells (TAMs) versus LCC2 cells (TAMr), LCC9 (TAMr and fulvestrant resistant), and MCF7:5C (TAMr) shows an upregulation over time in LCC1 cells only, with the resistant variants remaining largely unchanged during treatment with 1 nM E2

Figure 3 .
Figure 3. Modulating the PGE 2 axis influences cellular response to TAM.(a) Overexpression of HPGD reduced PGE 2 production in MCF-7 TAMr 2 cells to similar levels to basal production in WT MCF-7 cells; (b) exogenous PGE 2 (1 μM, broken lines) reduced sensitivity of WT MCF-7 and T47D human breast cancer cells to TAM; (c) PGE 1 -OH mimicked the effects PGE 2 in decreasing TAM sensitivity in WT MCF-7 and T47D cells, albeit to a lesser extent; (d) butaprost was ineffective at equivalent concentrations; (e) EP4 siRNA1 was selected for effective transient silencing of EP4 expression in WT MCF-7; (f) transient silencing of EP4 reduced sensitivity to TAM, only in the presence of PGE 2 .* denotes P < 0.05; ** denotes P < 0.002; *** denotes P < 0.005 normalised to respective vehicle controls

Figure 4 .
Figure 4. Kaplan-Meier survival analysis of the relationship between HPGD expression and breast cancer outcomes using in silico data mining and immunohistochemistry. Using KMplot, HPGD expression was categorised as 'high' (red) or 'low' (black), with the lower tertile expression used as cut-off.(a) Low HPGD expression is associated significantly reduced OS [HR 0.28 (0.13-0.59),P = 0.00041] in ERα breast tumours treated with TAM; (b) the association between HPGD expression and OS remained significant by multivariate analysis; (c) examples of HPGD staining (brown, DAB staining) in breast cancer tissue sections counterstained with haematoxylin (blue), the scale bar indicating 60 µm with some cases lacking any appreciable HPGD (top panels), and others showing weak to moderate staining (bottom panel); (d) a retrospective immunohistochemical analysis of HPGD expression in 130 breast cancers represented on TMAs showed HPGD-negative cases (black) had significantly worse OS [HR = 0.3 (0.15-0.69),P = 0.047] compared with HPGD-positive cases (red)

Figure 5 .
Figure 5. Relationship between HPGD and miR-3200-3p expression in breast cancer.Using KMplot, miR-3200-3p expression was categorized as high or low, with median expression used as cut-off.Kaplan-Meier survival analysis demonstrated a statistically significant effect of miR-3200-3p on the OS of ERα-positive breast cancer patients treated with any endocrine therapy (a); correlation analysis of matched gene expression levels of HPGD and miR-3200-3p in 987 breast cancer cases obtained from the TCGA revealed a weak negative correlation [R 2 = 0.25; (b)]

Table 2 .
[32]ilico prediction of miRNAs differentially expressed in TAMr MCF-7:5C vs. TAMs MCF-7:WS8 that are predicted to interact with 3'UTR of major HPGD transcripts (continued) Listed are miRNAs with fold change > 1.5; miRNAs underlined refer to the DLK-DIO3 locus on chromosome 14; miRWalk 3.0[31]and TargetScan 7.2[32]tools have been used to identify miRNA-target interactions; a major targeted transcripts are listed with respective Ensembl id; b binding probability and longest consecutive pairings are listed according to miRWalk 3.0 predictions; c additional miRNAs predicted only with TargetScan 7.2 tool (on most prevalent ENST00000296522.6 transcript)