# Dual epigenetic regulation of ERα36 expression in breast cancer cells

Charlène Thiébaut, Amand Chesnel, Jean-Louis Merlin, Maelle Chesnel, Agnès Leroux, Alexandre Harlé, Hélène Dumond

► To cite this version:

Charlène Thiébaut, Amand Chesnel, Jean-Louis Merlin, Maelle Chesnel, Agnès Leroux, et al. Dual epigenetic regulation of ERα36 expression in breast cancer cells. International Journal of Molecular Sciences, 2019, 20 (11), pp.2637. 10.3390/ijms20112637. [https://doi.org/10.3390/ijms20112637](https://doi.org/10.3390/ijms20112637)

HAL Id: hal-02145413

https://hal.science/hal-02145413

Submitted on 25 Oct 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Dual Epigenetic Regulation of ERα36 Expression in Breast Cancer Cells

Charlène Thiebaut<fnref n="1" />, Amand Chesnel<fnref n="1" />, Jean-Louis Merlin<fnref n="2" />, Maelle Chesnel<fnref n="1" />, Agnès Leroux<fnref n="3" />, Alexandre Harlé<fnref n="2" />* and Hélène Dumond <fnref n="1" />,*

Received: 21 March 2019; Accepted: 27 May 2019; Published: 29 May 2019

<figure><img src="image_3.png" /></figure>

check for
updates

**Abstract:** Breast cancer remains the major cause of cancer-induced morbidity and mortality in women. Among the different molecular subtypes, luminal tumors yet considered of good prognosis often develop acquired resistance to endocrine therapy. Recently, misregulation of ERα36 was reported to play a crucial role in this process. High expression of this ERα isoform was associated to preneoplastic phenotype in mammary epithelial cells, disease progression, and enhanced resistance to therapeutic agents in breast tumors. In this study, we identified two mechanisms that could together contribute to ERα36 expression regulation. We first focused on hsa-miR-136-5p, an ERα36 3'UTR-targeting microRNA, the expression of which inversely correlated to the ERα36 one in breast cancer cells. Transfection of hsa-miR136-5p mimic in MCF-7 cells resulted in downregulation of ERα36. Moreover, the demethylating agent decitabine was able to stimulate hsa-miR-136-5p endogenous expression, thus indirectly decreasing ERα36 expression and counteracting tamoxifen-dependent stimulation. The methylation status of ERα36 promoter also directly modulated its expression level, as demonstrated after decitabine treatment of breast cancer cell and confirmed in a set of tumor samples. Taken together, these results open the way to a direct and an indirect ERα36 epigenetic modulation by decitabine as a promising clinical strategy to counteract acquired resistance to treatment and prevent relapse.

**Keywords:** breast cancer; ERα36; microRNA; methylation; endocrine therapy resistance

## 1. Introduction

Despite considerable improvement of therapeutic strategies, breast cancer remains a leading cause of cancer mortality in women, mainly due to acquired resistance to treatment. Indeed, 30–50% of luminal tumors, yet considered as highly curable, develop endocrine therapy resistance through unclear and diverse mechanisms, including alterations of estrogen receptor (ER) gene sequence or epigenetic marks, alternative signaling pathway usage or variation in the pattern of ER isoform [1]. Recently, Wang and colleagues argued that one of the key actors of such acquired resistance could be the ER variant ERα36, the expression and/or activity of which become stimulated under estrogen deprivation or low dose tamoxifen exposure [2–5].

ERα36 is a 36kDa receptor encoded by the ESR1 locus and transcribed from an alternative promoter located into the first intron [6]. Compared to the canonical ERα66 protein, alternative splicing results in a shortened receptor lacking both AF-1 and AF-2 transactivation domains as well as the C-terminal part of the ligand-binding domain. ERα36 retains dimerization and DNA-binding domains and contains a unique 27 amino acid C-terminal sequence involved in cofactors interactions [7,8]. 17β-estradiol (E2)

<footnote n="1">¹ Université de Lorraine, CNRS, CRAN, F-54000 Nancy, France; charlene.thiebaut@univ-lorraine.fr (C.T.); amand.chesnel@univ-lorraine.fr (A.C.); maelle.chesnel22@gmail.com (M.C.)</footnote>

<footnote n="2">² Université de Lorraine, CNRS, CRAN, Institut de Cancérologie de Lorraine, F-54000 Nancy, France; jl.merlin@nancy.unicancer.fr (J.-L.M.); alexandre.harle@crick.ac.uk (A.H.)</footnote>

<footnote n="3">³ Institut de Cancérologie de Lorraine, F-54000 Nancy, France; a.leroux@nancy.unicancer.fr</footnote>

<footnote>* Correspondence: helene.dumond@univ-lorraine.fr</footnote>
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Strategy for the screening of microRNAs that target ERα36 but not ERα66.</figcaption></figure>

## 2.2. *hsa-miR-136-5p Under/over Expression Modulates ERα36 mRNA Expression Level*

To confirm that hsa-miR-136-5p modulate ERα36 expression, MCF-7 cells were transfected with either hsa-miR-136-5p mimic or inhibitor. ERα36 mRNA and hsa-miR-136-5p expression levels were measured by RT-qPCR analysis. The ratio of ERα36 expression in hsa-miR-136-5p mimic- or inhibitor-transfected versus control cells indicated that ERα36 mRNA expression level significantly decreased by 100% in the case of hsa-miR-136-5p mimic transfection compared to the hsa-miR-136-5p inhibitor one (Figure 2A). Transfection with hsa-miR-136-5p inhibitor led to an enhanced ERα36 expression (Figure 2A), but this tendency remained not statistically significant, probably due to experimental variation in inhibitor transfection efficacy (data not shown).

The specific interaction between hsa-miR-136-5p and the 3'UTR region of ERα36 mRNA was then assessed in both MCF-7 and MDA-MB-231 cells, which both express ERα36 but belong either to luminal A or to triple-negative molecular subtype, respectively (Figure S1). The plasmid pMIR-Report-Luc-ERα36wt (wt) contained the wild-type hsa-miR-136-5p target sequence of ERα36 3'UTR and pMIR-Report-Luc-ERα36mut (mut), a corresponding mutated sequence. Both the wt and mut sequence were cloned downstream the luciferase gene, thus driving luciferase mRNA stability in the presence of hsa-miR-136-5p. The cells were co-transfected with either the wt or the mut plasmid and hsa-miR-136-5p mimic or inhibitor. Luciferase activity was measured 24 h after transfection and normalized for transfection efficacy using the alkaline phosphatase activity from the co-transfected pEZX-PG04-Secreted alkaline phosphatase plasmid. Transfection of hsa-miR-136-5p mimic or inhibitor in MCF-7 or MDA-MB-231 cells resulted in a 120% to 150% decrease and a 180% to 240% increase of luciferase activity compared to control cells, respectively (Figure 2B). Both cell lines transfected by the mut plasmid displayed a higher basal level of luciferase activity than the wt one. There was no statistical variation between the hsa-miR-136-5p mimic or inhibitor transfected cells (Figure 2C). Therefore, we concluded that hsa-miR-136-5p directly and specifically targeted ERα36 mRNA, thus modulating its expression level in breast cancer cells.
<figure><img src="image_5.png" /><figcaption>**Figure 2.** hsa-miR-136-5p mimic/inhibitor directly modulates ERα36 expression. (a) MCF-7 cells were transfected with either hsa-miR-136-5p mimic or inhibitor. ERα36 expression was measured by RT-qPCR analysis 24 h after transfection, and the percentage of ERα36 expression variation in hsa-miR-136-5p mimic or inhibitor transfected cells versus control ones was calculated. MCF-7 or MDA-MB-231 cells were co-transfected with hsa-miR-136-5p mimic or inhibitor and (b) the pMIR-Report-Luc-ERα36wt plasmid or (c) the pMIR-Report-Luc-ERα36mut one. Luciferase activity was measured 24 h after transfection using the Bright-Glo™ luciferase assay system. To evaluate the transfection efficacy and normalize the results, cells were also transfected at the same time with the pEZX-PG04-Secreted alkaline phosphatase plasmid. The alkaline phosphatase activity was measured 24 h after transfection in the culture medium with the Secrete-Pair™ Gaussia luciferase dual and single luminescence assay Kit. Each bar represents mean ± S.E.M. (a) N = 5 (b) N = 3 and (c) N = 5. \* P < 0.05.</figcaption></figure>

## 2.3. *The Demethylating Agent DAC Stimulates hsa-miR-136-5p and Represses ERα36 Expression*

Many microRNAs are epigenetically regulated trough the methylation status of their promoter [30]. hsa-miR-136-5p is located in the DLK1-DIO3 imprinted locus of chromosome 14q32, which contains a cluster of noncoding RNAs and the antisense gene RTL1, all differentially regulated by the methylation status of the IG-DMR. Methylation changes in this locus, which may affect gene or microRNA expression associated with Temple or Kagami-Ogata syndrome and oligometastastic phenotype [31,32]. Therefore, we hypothesized that expression of hsa-miR-136-5p could be upregulated after cell exposure to the demethylating agent DAC and thus downregulate ERα36. The MCF-7 cells were treated for 48 h with DAC or DMSO as control. DAC exposure did not alter cell viability (data not shown). MicroRNAs and long mRNAs were harvested from the same cells. hsa-miR-136-5p expression level (Figure 3A) and ERα36 expression (Figure 3B) were measured using RT-qPCR. DAC treatment induced a significant 48% increase of hsa-miR-136-5p expression as well as a 59% decrease of ERα36 expression. The upregulation was independent of the hsa-miR-136-5p sponge CRNDE (colorectal neoplasia differentially expressed gene) (Figure S2).
<figure><img src="image_6.png" /><figcaption>**Figure 3.** The demethylating agent decitabine (DAC) stimulates hsa-miR-136-5p expression and represses ERα36 one. MCF-7 cells were treated for 48 h with DMSO 0.01% (as control) or the demethylating agent DAC 5 µM. (a) hsa-miR-136-5p expression level was assessed using TaqMan™ microRNA assay and (b) ERα36 expression was measured by RT-qPCR analysis. Each bar represents mean ± S.E.M. (a) N = 6 (b) N ≥ 3. * P < 0.05; ** P < 0.01.</figcaption></figure>

## 2.4. DAC Counteracts OHT-dependent Stimulation of ERα36 Expression

OHT was previously reported to either stimulate or repress ERα36 expression in MCF-7 cells exposed to 1 µM or 5 µM, respectively [5,33]. Since DAC (5 µM) and OHT (1 µM) displayed inverse effects on ERα36 expression, MCF-7 cells were exposed either to DAC or to OHT (1 µM), or co-treated with both compounds for 48 h. Exposure to 1 µM OHT-induced ERα36 expression by 8-fold although hsa-miR-136-5p expression was not affected by the treatment (data not shown). Moreover, co-treatment with DAC and OHT was able to counteract OHT-dependent stimulation of ERα36 expression (Figure 4A). Transfection of DAC/OHT co-treated cells with hsa-miR-136-5p inhibitor seemed to cancel the DAC-dependent impairment of OHT effect (Figure 4B), indicating that hsa-miR-136-5p was indeed a key element in this process. The variability of ERα36 expression observed in the [DAC/OHT + hsa-miR-136-5p inhibitor] condition was most likely due to uncontrolled experimental variations of the hsa-miR-136-5p inhibitor amount present into the cells. Taken together, these data suggested that DAC could be used to chemically repress ERα36 expression and prevent its induction by OHT through the stimulation of hsa-miR-136-5p.

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Decitabine (DAC) treatment counteracts 4-hydroxy-tamoxifen (OHT) induced ERα36 expression. (a) MCF-7 cells were treated for 48 h with DMSO (as control), OHT 1 µM or by a combination OHT 1 µM/DAC 5 µM. ERα36 expression was measured by RT-qPCR analysis. (b) MCF-7 cells were transfected or not with hsa-miR-136-5p inhibitor. Both control and hsa-miR-136-5p inhibitor transfected cells were treated for 48 h by a combination of OHT 1 µM/DAC 5 µM. ERα36 expression was measured by RT-qPCR analysis. OHT treatment was used as reference. Each bar represents mean ± S.E.M. N ≥ 3. * P < 0.05.</figcaption></figure>
## 2.5. *ERα36 Expression May be Modulated by the Methylation Status of Its Promoter, In Vivo and In Vitro*

Previous studies from our laboratory [16,34] suggested that ERα36 expression may be modulated by the methylation status of its promoter sequence, as defined by Wang's laboratory [14]. We searched for CpG islands located in the first intron of ESR1 using the EMBOSS CpG plot website (http://www.ebi.ac.uk/Tools/seqstats/emboss_cpgplot/) and identified four clusters (Figure S3). The methylation status of corresponding sequences was determined after bisulfite conversion in a subset of 17 breast primary tumors collected from the biobank of the Cancerology Institute of Lorraine (Vandoeuvre-lès-Nancy, France). Total RNAs were harvested from the same samples in order to measure ERα36 expression level by RT-qPCR. The correlation between each CpG methylation level and ERα36 expression level was addressed via a Pearson rank statistical test (Figure 5). The results clearly suggested that the methylation of 4 CpGs located on the island 4 was correlated to a high ERα36 expression.

<figure><img src="image_8.png" /><figcaption>**Figure 5.** ERα36 expression is related to the methylation status of its promoter in breast tumor samples. Seventeen breast tumor samples were collected from the tumor biobank of the Cancerology Institute of Lorraine. ERα36 promoter genomic DNA was submitted to bisulfite sequencing and ERα36 expression level was performed by real-time PCR analyses in triplicates for each tumor sample. The correlation between each CpG methylation and ERα36 expression was established by a MATLAB program (Pearson rank correlation test). \* CpG which methylation level is significantly (*P* < 0.05) correlated to ERα36 expression level.</figcaption></figure>

## 2.6. *DAC Directly Targets ERα36 Genomic Sequence*

Since the methylation status of ERα36 promoter seemed to be linked to its expression in tumor biopsies, we addressed a direct effect of the demethylating agent DAC on ERα36 expression through demethylation of its promoter. MCF-7 cells were exposed to DAC for 48 h, and the methylation status of the four CpGs significantly related to ERα36 expression was assessed after bisulfite conversion by real-time PCR. ERα66 is known to repress ERα36 expression and was reported to regulate its target genes through epigenetic modulator recruitment [35]. Thus, we compared the effect of the exposure to DAC or to the ERα66 inhibitor fulvestrant (Figure 6A). DAC triggered a 53% decrease in the methylation level of the ERα36 promoter sequence where ICI treatment had no effect. The same experiment was
performed in MCF-7 cells transfected by an ERα66 siRNA to knockdown ERα66 and exposed to DAC for 48 h (Figure 6B). The ERα66 knockdown was confirmed by Western blot (Figure 6C). DAC significantly and similarly demethylated ERα36 promoter in both sicontrol and siERα66 transfected cells. This confirmed that DAC directly acted on ERα36 promoter sequence independently of ERα66 and further suggested that the DAC-dependent repression of ERα36 expression could account for (i) an indirect effect through the hsa-miR-136-5p targeting of ERα36 3'UTR and (ii) a direct effect through demethylation of its promoter.

<figure><img src="image_9.png" /><figcaption>**Figure 6.** Decitabine directly targets ERα36 genomic sequence. (a) MCF-7 cells were treated for 48 h with DMSO (as control), DAC or ICI 182,780 (N = 4); (b) MCF-7 cells transfected with either scrambled siRNA (si control, N = 4) or siRNA targeting the exon 1 of ESR1 (siERα66, N = 3) were treated for 48 h with DMSO (as control) or DAC. The methylation level of the ERα36 promoter fourth CpG island was assessed after bisulfite conversion. A percentage of methylation level variation was calculated for DAC or ICI treated cells versus DMSO ones. (c) Representative Western blot analysis of ERα66 expression in MCF-7 cells transfected with either si control or si ERα66. α-Tubulin was used as a loading control. Each bar represents mean ± S.E.M. **P** < 0.01.</figcaption></figure>

**3. Discussion**

This study depicts, for the first time, two epigenetic mechanisms that regulate ERα36 expression in breast cancer cells. In line with the study of Li and colleagues [5], we showed that ERα36 expression increased significantly after a low dose tamoxifen exposure in MCF-7 cell line. This upregulation was previously associated with a decreased sensitivity to tamoxifen and enhanced proliferative, migratory, and invasive abilities of breast cancer cells [5,15]. Mechanistically, computational modeling and docking analysis revealed that tamoxifen is able to bind directly ERα36 through its ligand-binding domain [2]. This interaction activated various estrogen nongenomic pathways such as MAPK/Erk
and PI3K/Akt as well as the ERα36/EGFR/HER2 regulatory loop responsible for tamoxifen therapy escape [9,15,36]. These experimental results provided support and explanations to clinical studies wherein a high ERα36 expression in tumors treated with tamoxifen was significantly associated with a poor prognosis and an increased rate of metastases [2,37]. In this context, the understanding of mechanisms that regulate ERα36 expression was important to try preventing the deleterious effect of ERα36-triggered response to tamoxifen.

First, we showed in vitro and in breast tumor samples that a low methylation level of 4 CpGs localized in the promoter of ERα36 correlated with a low ERα36 mRNA expression. A demethylation of this CpG rich region was namely obtained in vitro after a 48 h DAC exposure which led to ERα36 downregulation. This supported the causal relationship between the methylation status of the promoter, and the expression level of ERα36. ERα66 has been described as the main direct transcriptional regulator of ERα36 expression and is able to regulate his target genes through the methylation level of their promoter [14,35,38,39]. When we silenced ERα66 expression or activity by either siRNAs or ICI 280 182, we showed that ERα66 did not drive the methylation-dependent repression of ERα36.

DAC is an FDA-approved hypomethylating agent, currently used to treat blood cancer. Recently, a study conducted by researchers from Florida on animal models demonstrated that a low dose of DAC prevented breast cancer cells from spreading and could also be used to treat aggressive chemoresistant triple-negative breast cancer [40,41]. As previously described in the case of ERβ silencing [42], treatment with DAC may promote ERα36 promoter demethylation and thus override the molecular mechanisms underlying ERα36-dependent tamoxifen resistance.

Second, our data indicated that DAC also acted indirectly on ERα36 expression since it upregulated hsa-miR-136-5p expression which, in turn, targeted ERα36 mRNA. hsa-miR-136-5p was described as a tumor suppressor microRNA, associated with impaired tumorigenesis and metastasis in prostate cancer by targeting MAP2K4 [43], renal cell carcinoma [23], lung squamous cell carcinoma by controlling metabolic processes [22], lung adenocarcinoma by suppressing cell adhesion gene expression [24], and glioma and cutaneous squamous cell carcinoma by targeting WNT signaling [44,45]. The downregulation of hsa-miR-136-5p has been described in a chemo-radioresistant group of patients with advanced stage cervical squamous cell carcinoma. Conversely, a high ERα36 expression was associated to estrogen nongenomic signaling driving metastatic progression, acquired resistance to hormone treatment, and poor outcome [11,13,46].

Our in vitro data demonstrated that hsa-miR-136-5p directly targeted ERα36 3'UTR and led to mRNA expression downregulation. These data suggested that the upregulation of ERα36 expression in response to OHT could be impaired either by the direct mimicking or by the DAC-dependent indirect stimulation of hsa-miR-136-5p that may consequently override ERα36-dependent resistance to OHT. Further in vivo experiments are needed to explore whether treatment with DAC, prior or in combination with OHT, would be a relevant therapeutic strategy.

In the mammary gland, hsa-miR-136-5p expression was higher in normal breast tissue than in breast carcinoma and in luminal vs. triple-negative breast cancer cell lines [18,19,47]. In line, Pascek and coworkers described a significant higher expression in luminal tumors than in TNBCs [48]. However, other studies from tumor biopsies did not report hsa-miR-136-5p expression as a reliable marker for predicting either breast cancer molecular type [49], pathological stage [50] or patient outcome [51]. Little is known about hsa-miR-136-5p regulation by DNA methylation status. The promoter sequence of hsa-miR-136-5p has not been clearly defined. A 135 bp sequence is located upstream hsa-miR-136-5p and separates hsa-miR-437-5p from hsa-miR-136-5p, but the promoter status of this sequence or the relationship between its methylation level and hsa-miR-136-5p expression are not described. Moreover, hsa-miR-136-5p is located in the 14q32.2 Dlk1-Dio3 chromosomal region, which has been described to be submitted to global and complex regulation by parental imprinting through two DMR [52]. Therefore, it would be interesting to determine precisely the regulatory elements of hsa-miR-136-5p and to set up a prospective study in order to investigate both its expression and the methylation status of its putative promoter in breast tumors.
Expression of hsa-miR-136-5p was also tightly regulated by a combination of sponge lncRNAs, such as CRNDE, HOX transcript antisense RNA (HOTAIR), lncRNA-small nucleolar RNA host gene 14 (SNHG14) or the circular RNA hsa_circ_0008309 [19,53–55]. Like in many solid tumors, the main factor to which expression could be relevant to predict patient outcome is CRNDE, the hsa-miR-136-5p sponge [19]. CRNDE was of peculiar interest in breast cancer since a recent meta-analysis confirmed this lncRNA could act as an oncogene which upregulated expression correlated with poor prognosis (late Tumor Node Metastasis grade) and tumor progression (high tumor size and lymph node metastasis) [56,57]. In vitro, CRNDE controls hsa-miR-136-5p level which, in turn, is able to target RASAL2 and Wnt/beta-catenin signaling in triple-negative breast cancer cells [18,19]. In the present study, we showed that hsa-miR-136-5p could modulate ERα36 expression independently of CRNDE. Thus, we added ERα36 as another piece to the puzzle of the complex epigenetic regulation of breast cancer progression.

# 4. Materials and Methods

## 4.1. Reagents

4-hydroxytamoxifen (OHT), fulvestrant (ICI 182, 780), and 5-Aza-2'-deoxycytidine (decitabine, DAC) were purchased from Sigma-Aldrich, France. Stock solutions (10 mM) of each compound were prepared in dimethyl sulfoxide (DMSO) and further diluted in culture medium for in vitro cell treatment. All working solutions were freshly made just before treatment, and control cells were treated with DMSO, diluted with the same factor.

## 4.2. Plasmids

Reporter plasmids containing the miR-136-5p target sequence from ERα36 3'UTR downstream the luciferase open reading frame were constructed as follows: pMIR-Report-Luc-KLF4-FL vector was purchased from the nonprofit Addgene repository (Addgene plasmid #34597) as a gift from Michael Ruppert [58]. The KLF4 sequence was removed and replaced by either the wild-type ERα36 3'UTR sequence (5'-CCACCATCAAATCAAATTGA-3') or the mutated one (5'-CCACCACTACAATAC AATTGA-3'). The corresponding plasmids were renamed pMIR-Report- Luc-ERα36wt and pMIR- Report-Luc-ERα36mut, respectively. pEZX-PG04-Secreted Alkaline Phosphatase plasmid was purchased from GeneCoppoeia™ (Tebu-bio, Le Perray-en-Yvelines, France).

## 4.3. Cell Culture

MCF-7 and MDA-MB-231 cells were purchased in 2015 from ATCC® (HTB-22™ and HTB-26™, respectively; Molsheim, France) and maintained in 10% fetal bovine serum (FBS) 1% Glutamine supplemented Dulbecco/Vogt modified Eagle's minimal essential medium (DMEM) (GIBCO) or Roswell Park Memorial Institute medium (RPMI) (GIBCO) medium, respectively.

## 4.4. Transfection

Transient transfections of MCF-7 or MDA-MB-231 cells were performed using JetPRIME transfection reagent (Ozyme, Saint-Cyr-l'École, France) in accordance with the manufacturer's instructions. Briefly, cells were plated at a density of 2.5 × 10⁵ cells per well in 6-well plates. After 24 h, transient transfection was performed with either 1 µg plasmid and/or 30 pmol (15 nM final) miRVana hsa-miR-136-5p mimic (5'-CAUCAUCAAACAAAUGGAGUTT-3') or inhibitor (5'-UCCAUCAUCAAACAAAUGGAGU-3') (Thermo Fischer, Illkirch, France) in serum-containing medium. The growth medium was renewed 4 hours after the transfection, and cells were treated with selected agents for 24 h.
### 4.5. Luciferase/Alkaline Phosphatase Assays

The luciferase activity associated with the plasmid pMIR-Report-Luc-ERα36wt and the plasmid pMIR-Report-Luc-ERα36mut was assessed using the Bright-Glo™ luciferase assay system (Promega Corporation, Madison, WI, USA) in accordance with the manufacturer's instructions. To normalize the results, pEZX-PG04-Secreted alkaline phosphatase plasmid was used, and the secreted alkaline phosphatase activity was detected on cell culture supernatant using the Secrete-Pair™ Gaussia luciferase dual and single luminescence assay kit (GeneCoppoeia™ Tebu-bio, Le Perray-en-Yvelines, France) in accordance with the manufacturer's instructions. Luminescence was measured by a microplate reader (Victor x3, Perkin-Elmer, Villebon-sur-Yvette, France).

### 4.6. MicroRNA Extraction and Reverse Transcription

Micro and long RNAs were isolated using the mirVana™ miRNA Isolation Kit (Ambion, Austin, TX, USA) in accordance with the manufacturer's instructions. The purity and quantity of RNA were assessed using the NanoDrop ND-2000C spectrophotometer (Thermo Scientific, Wilmington, DE, USA). The samples were used immediately or stored at -80 °C for future use. Reverse transcription (RT) and real-time PCR (qPCR) were carried out using TaqMan™ microRNA assay (Applied Biosystems, Foster City, CA, USA). RT reactions were performed in a volume of 15 µL, and each reaction contained 60 ng of microRNA. RT reactions were performed with a Mastercycler® Gradient (Eppendorf, Montesson, France) with the following conditions: 16 °C for 30 min, 42 °C for 30 min, 85 °C for 5 min, and 4 °C on hold. Reactions without addition of reverse transcriptase (RT (-) controls) were performed alongside with cDNA synthesis of each sample and used in subsequent procedures to control the potential genomic DNA contamination. A total of 1.4 µL of RT reaction product was added in the 20 µL qPCR reaction mix containing TaqMan™ Universal PCR Master Mix II without UNG (Applied Biosystems). All TaqMan assays were run in triplicate on a CFX96 Touch™ real-time PCR detection system (Biorad, Marnes-la-Coquette, France). Real-time PCR cycling conditions consisted of 95 °C for 10 min, followed by 50 cycles of 95 °C for 15 s and 60 °C for 1 min. TaqMan Probe (hsa-miR-136-5p: 5'-UGAGCCCUCGGAGGACUCCAUUUGUUUUGAUGAUGGAUUCUUAUGCUCCAUCAUCG UCUCAAAUGAGUCUUCUGUGGGUUCU-3'; hsa-miR-423-5p: 5'-UGAGGGGCAGAGAGCGAGAC UUUUCUAUUUUCCAAAAGCUCGGUCUGAGGCCCCUCAGU-3'). Assays were performed at least as independent experimental triplicates, and the mean values were used to calculate expression levels, using the ΔΔC(t) method referring to hsa-miR-423 housekeeping miR expression (as recommended by the manufacturer, Applied Biosystems).

### 4.7. Long RNAs RT and Real-time PCR Analysis

RT and real-time PCR analyses were performed as previously described [16]. The following primers were used for qRT-PCR: GAPDH forward (Fw) 5'-TGC-ACC-ACC-AAC-TGC-TTA-GC-3', GAPDH reverse (Rev) 5'-GGC-ATG-GAC-TGT-GGT-CAT-GAG-3', ERα36 forward (Fw) 5'-ATG-AAT-CTG-CAG-GGA-GAG-GA-3', ERα36 reverse (Rev) 5'-GGC-TTT-AGA-CAC-GAG-GAA-ACC-3', CRNDE forward (Fw) 5'-ATA-TTC-AGC-CGT-TGG-TCT-TTG-A-3', and CRNDE reverse (Rev) 5'-ATA-TTC-AGC-CGT-TGG-TCT-TTG-A-3'. Assays were performed at least in triplicate, and the mean values were used to calculate expression levels, using the ΔΔC(t) method referring to GAPDH housekeeping gene expression.

### 4.8. Western Immunoblotting

Western blots were performed as described previously [10]. We used the primary antibody anti-ERα66 (F10, Santa Cruz Biotechnology, Heidelberg, Germany). The anti-α-Tubulin (GTX102079, Euromedex, Souffelweyersheim, France) was used as a loading control. Protein expression profiles were revealed with Clarity Western ECL Substrate (Biorad), and banding quantification was performed using the Image Lab 6.0.0 software (Biorad).
## 4.9. Patients and Samples

Cryopreserved samples (Biobank registration number AC-2008-174) from 17 patients treated at Cancerology Institute of Lorraine (Nancy, France) between 2014 and 2015 with histology proven T1 to T4 invasive breast cancer were used for this study. Among the 17 patients, 14 (82%) had an invasive ductal carcinoma and 3 (18%) had a nonspecific invasive carcinoma. Breast tumors were examined by hematoxylin and eosin staining and histopathology defined in accordance with the sixth edition of TNM [59]. Hormonal receptors, HER2 and Ki-67 expressions were assessed using immunohistochemistry (Dako, Les Ulis, France). Patients' characteristics are detailed in Table S1. The investigations were carried out following the rules of the Declaration of Helsinki of 1975 (https://www.wma.net/what-we-do/medical-ethics/declaration-of-helsinki/), revised in 2013. All patients gave their informed consent for the use of their samples, and all data were anonymized prior to analysis to protect patient confidentiality. This study was performed within the framework of the tumor banking program of the Cancerology Institute of Lorraine including ethical considerations (approval AC-2013-1916 French Ministry of Research).

## 4.10. Bisulfite Sequencing

Bisulfite sequencing data of tumor samples were obtained from Active Motif's Targeted Next-Gen Bisulfite Sequencing Service (Active Motif Europe, La Hulpe, Belgium).

## 4.11. DNA Methylation Analysis or ERα36 Promoter Sequence

Genomic DNA was harvested from MCF-10A, MCF-7 or MDA-MB-231, and bisulfite conversion of DNA was performed with the bisulfite conversion kit (Active Motif Europe, La Hulpe, Belgium) in accordance with the manufacturer's instructions. Real-Time PCR analysis of methylated/unmethylated DNA was performed after conversion using specific primers (forward 5'-GAG-TTT-AAA-ATA-AGT-TAT-ATG-GAA-GTA-TAA-GTG-3'; Meth reverse 5'-CGA-CCG-CGC-TCC-TTC-CAC-AAT-AAC-TAC-GA-3'; Unmeth reverse 5'-CAA-CCA-CAC-TCC-TTC-CAC-AAT-AAC-TAC-AA-3').

## 4.12. Statistical Analysis

Statistical analyses were performed with MATLAB R2018a software (MathWorks) using Student's t-test for unpaired samples and ANOVA followed by Bonferroni post-hoc test for multigroup comparison, with a significant *p*-value threshold below 5%. \* *p* < 0.05; \*\* *p* < 0.01; \*\*\* *p* < 0.001. Standard deviations or standard errors were indicated on figures as advocated by Altman and Bland [60].

**Supplementary Materials:** Supplementary materials can be found at http://www.mdpi.com/1422-0067/20/11/2637/s1.

**Author Contributions:** Conceptualization, H.D.; Methodology, H.D., A.C. and A.L.; Validation, H.D., A.H., A.C. and J.L.M.; Investigation, C.T. and M.C.; Resources, A.L.; Writing—Original Draft Preparation, C.T. and A.H.; Writing—Review, ανδ Editing, H.D. and A.H.; Supervision, A.H. and H.D.; Project Administration, H.D.; Funding Acquisition, H.D. and J.L.M.

**Funding:** This research was funded by the Agence nationale de sécurité sanitaire de l'alimentation, de l'environnement et du travail (ANSES; PNR EST 2012-2-014), INSERM Cancer-Environnement (ENV201304) and La Ligue contre le Cancer (2016).

**Acknowledgments:** We thank Audrey Taverniers for technical assistance.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Szostakowska, M.; Trębińska-Stryjewska, A.; Grzybowska, E.A.; Fabisiewicz, A. Resistance to endocrine therapy in breast cancer: Molecular mechanisms and future goals. *Breast Cancer Res. Treat.* **2019**, *173*, 489–497. [CrossRef]
2. Wang, Q.; Jiang, J.; Ying, G.; Xie, X.-Q.; Zhang, X.; Xu, W.; Zhang, X.; Song, E.; Bu, H.; Ping, Y.-F.; et al. Tamoxifen enhances stemness and promotes metastasis of ERα36+ breast cancer by upregulating ALDH1A1 in cancer cells. *Cell Res.* **2018**, *28*, 336–358. [CrossRef]

3. Wang, Z.; Zhang, X.; Shen, P.; Loggie, B.W.; Chang, Y.; Deuel, T.F. A variant of estrogen receptor-α, hER-α36: Transduction of estrogen- and antiestrogen-dependent membrane-initiated mitogenic signaling. *Proc. Natl. Acad. Sci. USA* **2006**, *103*, 9063–9068. [CrossRef]

4. Zhang, X.; Wang, Z.-Y. Estrogen receptor-α variant, ER-α36, is involved in tamoxifen resistance and estrogen hypersensitivity. *Endocrinology* **2013**, *154*, 1990–1998. [CrossRef]

5. Li, G.; Zhang, J.; Jin, K.; He, K.; Zheng, Y.; Xu, X.; Wang, H.; Wang, H.; Li, Z.; Yu, X.; et al. Estrogen receptor-α36 is involved in development of acquired tamoxifen resistance via regulating the growth status switch in breast cancer cells. *Mol. Oncol.* **2013**, *7*, 611–624. [CrossRef] [PubMed]

6. Wang, Z.; Zhang, X.; Shen, P.; Loggie, B.W.; Chang, Y.; Deuel, T.F. Identification, cloning, and expression of human estrogen receptor-alpha36, a novel variant of human estrogen receptor-alpha66. *Biochem. Biophys. Res. Commun.* **2005**, *336*, 1023–1027. [CrossRef]

7. Chaudhri, R.A.; Hadadi, A.; Lobachev, K.S.; Schwartz, Z.; Boyan, B.D. Estrogen receptor-alpha 36 mediates the anti-apoptotic effect of estradiol in triple negative breast cancer cells via a membrane-associated mechanism. *Biochim. Biophys. Acta* **2014**, *1843*, 2796–2806. [CrossRef]

8. Omarjee, S.; Jacquemetton, J.; Poulard, C.; Rochel, N.; Dejaegere, A.; Chebaro, Y.; Treilleux, I.; Marangoni, E.; Corbo, L.; Romancer, M.L. The molecular mechanisms underlying the ERα-36-mediated signaling in breast cancer. *Oncogene* **2017**, *36*, 2503–2514. [CrossRef]

9. Chaudhri, R.A.; Schwartz, N.; Elbaradie, K.; Schwartz, Z.; Boyan, B.D. Role of ERα36 in membrane-associated signaling by estrogen. *Steroids* **2014**, *81*, 74–80. [CrossRef] [PubMed]

10. Thiebaut, C.; Chamard-Jovenin, C.; Chesnel, A.; Morel, C.; Djermoune, E.-H.; Boukhobza, T.; Dumond, H. Mammary epithelial cell phenotype disruption in vitro and in vivo through ERalpha36 overexpression. *PLoS ONE* **2017**, *12*, e0173931. [CrossRef] [PubMed]

11. Wang, Z.-Y.; Yin, L. Estrogen receptor alpha-36 (ER-α36): A new player in human breast cancer. *Mol. Cell. Endocrinol.* **2015**, *418*, 193–206. [CrossRef]

12. Soltysik, K.; Czekaj, P. ERα36–Another piece of the estrogen puzzle. *Eur. J. Cell Biol.* **2015**, *94*, 611–625. [CrossRef]

13. Chamard-Jovenin, C.; Jung, A.C.; Chesnel, A.; Abecassis, J.; Flament, S.; Ledrappier, S.; Macabre, C.; Boukhobza, T.; Dumond, H. From ERα66 to ERα36: A generic method for validating a prognosis marker of breast tumor progression. *BMC Syst. Biol.* **2015**, *9*, 28. [CrossRef]

14. Zou, Y.; Ding, L.; Coleman, M.; Wang, Z. Estrogen receptor-alpha (ER-α) suppresses expression of its variant ER-α36. *FEBS Lett.* **2009**, *583*, 1368–1374. [CrossRef]

15. Yin, L.; Wang, Z.-Y. Roles of the ER-α36-EGFR/HER2 positive regulatory loops in tamoxifen resistance. *Steroids* **2016**, *111*, 95–99. [CrossRef]

16. Chamard-Jovenin, C.; Thiebaut, C.; Chesnel, A.; Bresso, E.; Morel, C.; Smail-Tabbone, M.; Devignes, M.-D.; Boukhobza, T.; Dumond, H. Low-Dose Alkylphenol Exposure Promotes Mammary Epithelium Alterations and Transgenerational Developmental Defects, But Does Not Enhance Tumorigenic Behavior of Breast Cancer Cells. *Front. Endocrinol.* **2017**, *8*, 272. [CrossRef]

17. Piletič, K.; Kunej, T. MicroRNA epigenetic signatures in human disease. *Arch. Toxicol.* **2016**, *90*, 2405–2419. [CrossRef]

18. Yan, M.; Li, X.; Tong, D.; Han, C.; Zhao, R.; He, Y.; Jin, X. miR-136 suppresses tumor invasion and metastasis by targeting RASAL2 in triple-negative breast cancer. *Oncol. Rep.* **2016**, *36*, 65–71. [CrossRef]

19. Huan, J.; Xing, L.; Lin, Q.; Xui, H.; Qin, X. Long noncoding RNA CRNDE activates Wnt/β-catenin signaling pathway through acting as a molecular sponge of microRNA-136 in human breast cancer. *Am. J. Transl. Res.* **2017**, *9*, 1977–1989.

20. Guo, T.; Pan, G. MicroRNA-136 functions as a tumor suppressor in osteosarcoma via regulating metadherin. *Cancer Biomark. Sect. Dis. Markers* **2018**, *22*, 79–87. [CrossRef]

21. Deng, G.; Gao, Y.; Cen, Z.; He, J.; Cao, B.; Zeng, G.; Zong, S. miR-136-5p Regulates the Inflammatory Response by Targeting the IKKβ/NF-κB/A20 Pathway After Spinal Cord Injury. *Cell. Physiol. Biochem. Int. J. Exp. Cell. Physiol. Biochem. Pharmacol.* **2018**, *50*, 512–524. [CrossRef]
22. Xie, Z.-C.; Li, T.-T.; Gan, B.-L.; Gao, X.; Gao, L.; Chen, G.; Hu, X.-H. Investigation of miR-136-5p key target genes and pathways in lung squamous cell cancer based on TCGA database and bioinformatics analysis. *Pathol. Res. Pract.* **2018**, *214*, 644–654. [CrossRef]

23. Chen, P.; Zhao, L.; Pan, X.; Jin, L.; Lin, C.; Xu, W.; Xu, J.; Guan, X.; Wu, X.; Wang, Y.; et al. Tumor suppressor microRNA-136-5p regulates the cellular function of renal cell carcinoma. *Oncol. Lett.* **2018**, *15*, 5995–6002. [CrossRef]

24. Li, T.-T.; Gao, X.; Gao, L.; Gan, B.-L.; Xie, Z.-C.; Zeng, J.-J.; Chen, G. Role of upregulated miR-136-5p in lung adenocarcinoma: A study of 1242 samples utilizing bioinformatics analysis. *Pathol. Res. Pract.* **2018**, *214*, 750–766. [CrossRef]

25. Kozomara, A.; Griffiths-Jones, S. miRBase: Annotating high confidence microRNAs using deep sequencing data. *Nucleic Acids Res.* **2014**, *42*, D68–D73. [CrossRef]

26. Kozomara, A.; Griffiths-Jones, S. miRBase: Integrating microRNA annotation and deep-sequencing data. *Nucleic Acids Res.* **2011**, *39*, D152–D157. [CrossRef]

27. Griffiths-Jones, S.; Grocock, R.J.; van Dongen, S.; Bateman, A.; Enright, A.J. miRBase: microRNA sequences, targets and gene nomenclature. *Nucleic Acids Res.* **2006**, *34*, D140–D144. [CrossRef]

28. Griffiths-Jones, S.; Saini, H.K.; van Dongen, S.; Enright, A.J. miRBase: Tools for microRNA genomics. *Nucleic Acids Res.* **2008**, *36*, D154–D158. [CrossRef]

29. Panwar, B.; Omenn, G.S.; Guan, Y. miRmine: A database of human miRNA expression profiles. *Bioinform. Oxf. Engl.* **2017**, *33*, 1554–1560. [CrossRef]

30. Gulyaeva, L.F.; Kushlinskiy, N.E. Regulatory mechanisms of microRNA expression. *J. Transl. Med.* **2016**, *14*, 143. [CrossRef]

31. Eggermann, T.; Perez de Nanclares, G.; Maher, E.R.; Temple, I.K.; Tümer, Z.; Monk, D.; Mackay, D.J.G.; Grønskov, K.; Riccio, A.; Linglart, A.; et al. Imprinting disorders: A group of congenital disorders with overlapping patterns of molecular changes affecting imprinted loci. *Clin. Epigenetics* **2015**, *7*, 123. [CrossRef]

32. Uppal, A.; Wightman, S.C.; Mallon, S.; Oshima, G.; Pitroda, S.P.; Zhang, Q.; Huang, X.; Darga, T.E.; Huang, L.; Andrade, J.; et al. 14q32-encoded microRNAs mediate an oligometastatic phenotype. *Oncotarget* **2015**, *6*, 3540–3552. [CrossRef]

33. Xu, Z.; Huang, B.; Liu, J.; Wu, X.; Luo, N.; Wang, X.; Zheng, X.; Pan, X. Combinatorial anti-proliferative effects of tamoxifen and naringenin: The role of four estrogen receptor subtypes. *Toxicology* **2018**, *410*, 231–246. [CrossRef]

34. Ajj, H.; Chesnel, A.; Pinel, S.; Plénat, F.; Flament, S.; Dumond, H. An alkylphenol mix promotes seminoma derived cell proliferation through an ERalpha36-mediated mechanism. *PLoS ONE* **2013**, *8*, e61758. [CrossRef]

35. Ariazi, E.A.; Taylor, J.C.; Black, M.A.; Nicolas, E.; Slifker, M.J.; Azzam, D.J.; Boyd, J. A New Role for ERα: Silencing via DNA Methylation of Basal, Stem Cell, and EMT Genes. *Mol. Cancer Res.* **2017**, *15*, 152–164. [CrossRef]

36. Chaudhri, R.A.; Olivares-Navarrete, R.; Cuenca, N.; Hadadi, A.; Boyan, B.D.; Schwartz, Z. Membrane estrogen signaling enhances tumorigenesis and metastatic potential of breast cancer cells via estrogen receptor-α36 (ERα36). *J. Biol. Chem.* **2012**, *287*, 7169–7181. [CrossRef]

37. Shi, L.; Dong, B.; Li, Z.; Lu, Y.; Ouyang, T.; Li, J.; Wang, T.; Fan, Z.; Fan, T.; Lin, B.; et al. Expression of ER-{alpha}36, a novel variant of estrogen receptor {alpha}, and resistance to tamoxifen treatment in breast cancer. *J. Clin. Oncol.* **2009**, *27*, 3423–3429. [CrossRef]

38. Fleischer, T.; Tekpli, X.; Mathelier, A.; Wang, S.; Nebdal, D.; Dhakal, H.P.; Sahlberg, K.K.; Schlichting, E.; Børresen-Dale, A.-L.; Borgen, E.; et al. DNA methylation at enhancers identifies distinct breast cancer lineages. *Nat. Commun.* **2017**, *8*, 1379. [CrossRef]

39. Xue, X.; Yang, Y.A.; Zhang, A.; Fong, K.; Kim, J.; Song, B.; Li, S.; Zhao, J.C.; Yu, J. LncRNA HOTAIR enhances ER signaling and confers tamoxifen resistance in breast cancer. *Oncogene* **2016**, *35*, 2746–2755. [CrossRef]

40. Borges, S.; Döppler, H.; Perez, E.A.; Andorfer, C.A.; Sun, Z.; Anastasiadis, P.Z.; Thompson, E.A.; Geiger, X.J.; Storz, P. Pharmacologic reversion of epigenetic silencing of the PRKD1 promoter blocks breast tumor cell invasion and metastasis. *Breast Cancer Res.* **2013**, *15*, R66. [CrossRef]

41. Yu, J.; Qin, B.; Moyer, A.M.; Nowsheen, S.; Liu, T.; Qin, S.; Zhuang, Y.; Liu, D.; Lu, S.W.; Kalari, K.R.; et al. DNA methyltransferase expression in triple-negative breast cancer predicts sensitivity to decitabine. *J. Clin. Invest.* **2018**, *128*, 2376–2388. [CrossRef]
42. Pitta, C.A.; Papageorgis, P.; Charalambous, C.; Constantinou, A.I. Reversal of ER-β silencing by chromatin modifying agents overrides acquired tamoxifen resistance. *Cancer Lett.* **2013**, *337*, 167–176. [CrossRef]

43. Zhu, Y.; Shao, S.; Pan, H.; Cheng, Z.; Rui, X. MicroRNA-136 inhibits prostate cancer cell proliferation and invasion by directly targeting mitogen-activated protein kinase kinase 4. *Mol. Med. Rep.* **2018**, *17*, 4803–4810. [CrossRef]

44. Li, D.-X.; Fei, X.-R.; Dong, Y.-F.; Cheng, C.-D.; Yang, Y.; Deng, X.-F.; Huang, H.-L.; Niu, W.-X.; Zhou, C.-X.; Xia, C.-Y.; et al. The long non-coding RNA CRNDE acts as a ceRNA and promotes glioma malignancy by preventing miR-136-5p-mediated downregulation of Bcl-2 and Wnt2. *Oncotarget* **2017**, *8*, 88163–88178. [CrossRef]

45. Singh, A.; Willems, E.; Singh, A.; Ong, I.M.; Verma, A.K. Ultraviolet radiation-induced differential microRNA expression in the skin of hairless SKH1 mice, a widely used mouse model for dermatology research. *Oncotarget* **2016**, *7*, 84924–84937. [CrossRef]

46. Wallacides, A.; Chesnel, A.; Ajj, H.; Chillet, M.; Flament, S.; Dumond, H. Estrogens promote proliferation of the seminoma-like TCam-2 cell line through a GPER-dependent ERα36 induction. *Mol. Cell. Endocrinol.* **2012**, *350*, 61–71. [CrossRef]

47. Iorio, M.V.; Ferracin, M.; Liu, C.-G.; Veronese, A.; Spizzo, R.; Sabbioni, S.; Magri, E.; Pedriali, M.; Fabbri, M.; Campiglio, M.; et al. MicroRNA gene expression deregulation in human breast cancer. *Cancer Res.* **2005**, *65*, 7065–7070. [CrossRef]

48. Paszek, S.; Gabło, N.; Barnaś, E.; Szybka, M.; Morawiec, J.; Kołacińska, A.; Zawlik, I. Dysregulation of microRNAs in triple-negative breast cancer. *Ginekol. Pol.* **2017**, *88*, 530–536. [CrossRef]

49. Braicu, C.; Raduly, L.; Morar-Bolba, G.; Cojocneanu, R.; Jurj, A.; Pop, L.-A.; Pileczki, V.; Ciocan, C.; Moldovan, A.; Irimie, A.; et al. Aberrant miRNAs expressed in HER-2 negative breast cancers patient. *J. Exp. Clin. Cancer Res.* **2018**, *37*, 257. [CrossRef]

50. Yerukala Sathipati, S.; Ho, S.-Y. Identifying a miRNA signature for predicting the stage of breast cancer. *Sci. Rep.* **2018**, *8*, 16138. [CrossRef]

51. Cheng, D.; He, H.; Liang, B. A three-microRNA signature predicts clinical outcome in breast cancer patients. *Eur. Rev. Med. Pharmacol. Sci.* **2018**, *22*, 6386–6395.

52. Kagami, M.; O'Sullivan, M.J.; Green, A.J.; Watabe, Y.; Arisaka, O.; Masawa, N.; Matsuoka, K.; Fukami, M.; Matsubara, K.; Kato, F.; Ferguson-Smith, A.C.; Ogata, T. The IG-DMR and the MEG3-DMR at human chromosome 14q32.2: Hierarchical interaction and distinct functional properties as imprinting control centers. *PLoS Genet.* **2010**, *6*, e1000992. [CrossRef]

53. Zhu, L.; Liu, Y.; Chen, Q.; Yu, G.; Chen, J.; Chen, K.; Yang, N.; Zeng, T.; Yan, S.; Huang, A.; et al. Long-Noncoding RNA Colorectal Neoplasia Differentially Expressed Gene as a Potential Target to Upregulate the Expression of IRX5 by miR-136-5P to Promote Oncogenic Properties in Hepatocellular Carcinoma. *Cell. Physiol. Biochem.* **2018**, *50*, 2229–2248. [CrossRef]

54. Li, B.; Wang, F.; Li, X.; Sun, S.; Shen, Y.; Yang, H. Hsa_circ_0008309 May Be a Potential Biomarker for Oral Squamous Cell Carcinoma. *Dis. Markers* **2018**, *2018*, 7496890. [CrossRef]

55. Duan, C.; Liu, Y.; Li, Y.; Chen, H.; Liu, X.; Chen, X.; Yue, J.; Zhou, X.; Yang, J. Sulfasalazine alters microglia phenotype by competing endogenous RNA effect of miR-136-5p and long non-coding RNA HOTAIR in cuprizone-induced demyelination. *Biochem. Pharmacol.* **2018**, *155*, 110–123. [CrossRef]

56. Dai, M.; Li, S.; Qin, X. Colorectal neoplasia differentially expressed: A long noncoding RNA with an imperative role in cancer. *OncoTargets Ther.* **2018**, *11*, 3755–3763. [CrossRef]

57. Wang, W.; Yuan, F.; Xu, J. The prognostic role of long noncoding RNA CRNDE in cancer patients: A systematic review and meta-analysis. *Neoplasma* **2019**, *66*, 73–82. [CrossRef]

58. Lin, C.-C.; Liu, L.-Z.; Addison, J.B.; Wonderlin, W.F.; Ivanov, A.V.; Ruppert, J.M. A KLF4-miRNA-206 autoregulatory feedback loop can promote or inhibit protein translation depending upon cell context. *Mol. Cell. Biol.* **2011**, *31*, 2513–2527. [CrossRef]

59. Sobin, L.H.; Wittekind, C.H. *TNM classification of malignant tumors, 6th ed*; John Wiley & Sons: Hoboken, NJ, USA, 2002.

60. Altman, D.G.; Bland, J.M. Standard deviations and standard errors. *BMJ* **2005**, *331*, 903. [CrossRef]
<figure><img src="image_10.png" /></figure>

© 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).