Figures
Abstract
Thyroid hormone receptor beta (THRB) gene is commonly deregulated in cancers and, as strengthened by animal models, postulated to play a tumor-suppressive role. Our previous studies revealed downregulation of THRB in clear cell renal cell carcinoma (ccRCC), but the culpable mechanisms have not been fully elucidated. Since epigenetic regulation is a common mechanism influencing the expression of tumor suppressors, we hypothesized that downregulation of THRB in renal cancer results from epigenetic aberrances, including CpG methylation and microRNA-dependent silencing. Our study revealed that ccRCC tumors exhibited a 56% decrease in THRB and a 37% increase in DNA methyltransferase 1 (DNMT1) expression when compared with paired non-neoplastic control samples. However, THRB CpG methylation analysis performed using BSP, SNaPshot and MSP-PCR consistently revealed no changes in methylation patterns between matched tumor and control samples. In silico analysis resulted in identification of four microRNAs (miR-155, miR-425, miR-592, and miR-599) as potentially targeting THRB transcript. Luciferase assay showed direct binding of miR-155 and miR-425 to 3′UTR of THRB, and subsequent in vivo analyses revealed that transfection of UOK171 cell line with synthetic miR-155 or miR-425 resulted in decreased expression of endogenous TRHB by 22% and 64%, respectively. Finally, real-time PCR analysis showed significant upregulation of miR-155 (354%) and miR-425 (162%) in ccRCC when compared with matched controls. Moreover, microRNA levels were negatively correlated with the amount of THRB transcript in tissue samples. We conclude that CpG methylation is not the major mechanism contributing to decreased THRB expression in ccRCC. In contrast, THRB is targeted by microRNAs miR-155 and miR-425, whose increased expression may be responsible for downregulation of THRB in ccRCC tumors.
Citation: Wojcicka A, Piekielko–Witkowska A, Kedzierska H, Rybicka B, Poplawski P, Boguslawska J, et al. (2014) Epigenetic Regulation of Thyroid Hormone Receptor Beta in Renal Cancer. PLoS ONE 9(5): e97624. https://doi.org/10.1371/journal.pone.0097624
Editor: Javier S. Castresana, University of Navarra, Spain
Received: December 18, 2013; Accepted: April 23, 2014; Published: May 21, 2014
Copyright: © 2014 Wojcicka et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: The study was supported by the National Science Centre grants NN401611940 and NN401047638 (to AN), the UNESCO/FEBS Collaborative Experimental Scholarship for Central & Eastern Europe, 2010 (to AW), and the National Science Centre Grant 2012/05/B/NZ5/01541 (to APW). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Thyroid hormone receptors (TRs): TRα and TRβ, are ligand (3,5,3′-triiodothyronine, T3)-inducible transcription factors that mediate the cellular effects of thyroid hormones, namely its active form – T3. Since T3-dependent genes include numerous important regulators of cell cycle, such as mdm2, p53 and retinoblastoma [1],[2], actions of TRs contribute to maintenance of key cellular processes including proliferation, differentiation, apoptosis and metabolism [3]. TRs include 3 functional proteins: TRα1, TRβ1 and TRβ2, encoded by THRA and THRB genes. TRα1 and TRβ1 receptors are expressed in virtually all tissues, but their expression levels and functionality depend on the cell type and developmental stage of an organism. Disturbed activity of TRs, resulting from aberrances in expression or sequences of genes coding for TRs is a common phenomenon in human cancers. It leads to disturbed expression of T3-dependent genes and, in consequence, to severe impairment in cellular homeostasis. Those observations led to hypothesis on the tumor suppressive role of TRβ, further confirmed in several elegant studies performed in cell lines and mouse models [4],[5].
One of the cancers in which aberrances in TRβ function are frequently observed, is the most common type of renal tumors - clear cell renal cell carcinoma (ccRCC). Most interestingly, THRB gene resides within 3p21-25 chromosomal region, known as a hot spot for mutations in genes involved in ccRCC pathogenesis [6]. The culpable mechanisms identified so far include mutations, aberrant expression and deregulated splicing [7],[8],[9].
It has been shown that ccRCC is accompanied by epigenetic aberrances that may directly contribute to tumorigenic process, acting at early and precancerous phase of the multistage renal tumorigenesis [10]. Numerous recent studies indicate that epigenetic deregulation is among major causes of derailed actions of tumor suppressors in cancers, and most frequently identified mechanisms include DNA methylation and microRNAs. DNA methylation consists in addition of a methyl group to a cytosine preceding guanine (CpG) in DNA sequence, and is catalyzed by DNA methyltransferases (DNMTs): DNMT1, DNMT3A and DNMT3B. DNMT3A and DNMT3B are de novo methyltransferases, expressed prevalently during early stages of development [11] that sustain parental methylation patterns. DNMT1 is referred to as the maintenance methyltransferase, as it is responsible for maintaining methylation patterns during DNA replication [12]. Aberrant methylation in cancer is caused by impaired expression and function of DNMTs. Cancer cells usually exhibit overall genomic hypomethylation, leading to microsatellite instability and activation of genes involved in metastatic changes [13],[14]. Simultaneously, severe hypermethylation of regulatory regions of tumor supressors, DNA repair and cell cycle control genes leads to progression of cancer [15]. For ccRCC it was shown that numerous cases of downregulated expression of VHL tumor suppressor, inactivated in approx. 50% of ccRCC patients, result from hypermethylation of the gene’s promoter [16].
MicroRNAs (miRNAs) are short RNAs that inhibit expression of protein-coding genes by binding to complementary sequences in 3′UTRs (Untranslated regions) of their transcripts. The sequence responsible for this recognition encompasses nucleotides 2–8 of a miRNA and has to be fully complementary to the target sequence in mRNA [17]. Numerous cancers exhibit aberrant expression of microRNAs, leading to severe changes in the transcriptome and proteome of a cancer cell. For ccRCC, a number of microRNAs were shown to be deregulated in tumor [18]–[21]. Both DNA methylation and microRNA –dependent regulation have recently brought the attention as two mechanism with greatest potential for possible therapeutic interventions [22]. It was shown that tumor suppressor genes, silenced due to DNA hypermethylation, can be reactivated with demethylating agents. Furthermore, microRNA mimics and inhibitors are currently evaluated as therapeutic molecules that show promise in personalized anticancer medicine [23].
The role of epigenetic changes in THRB disturbances in cancer has been investigated in a few studies. The potential role of microRNA-dependent regulation in THRB silencing was suggested for ccRCC [9] and proven for papillary thyroid carcinoma [24]. Several other studies suggested DNA hypermethylation as a mechanism contributing to THRB silencing in leukemia and cancers of breast, lung, and thyroid [25]–[30].
Those studies as well as frequent epigenetic deregulations observed in ccRCC [31], prompted us to analyze the possible impact of DNA methylation and microRNA-dependent regulation on THRB expression in renal cancer.
Materials and Methods
Material
Tissue samples were obtained with the permission of the Bioethics Committee of the Centre of Postgraduate Medical Education in Warsaw from patients with clear cell renal cell carcinoma (ccRCC) (Table S1). Written informed consent was obtained from all patients. The samples were divided into two groups: cancer tissues (n = 35, T) and control tissues (paired normal tissue from the opposite pole of the malignant kidney (>5 cm from tumor) with no histological evidence of tumor; n = 35, C). ccRCC was diagnosed by histology according to WHO criteria [32].
The following cell lines were used in the study: 1) HK-2 (ATCC no.: CRL-2190), a proximal tubular cell line derived from normal kidney, immortalized by transduction with human papilloma virus 16 (HPV-16) E6/E7 genes. 2) UOK171 (UOB Tumor Cell Line Repository, provided by Dr. W Marston Linehan, NIH, NCI, Bethesda, USA, patent E-033-2010/0): derived from IV stage ccRCC tumor, spontaneously immortalized. 3) HeLa (ATTC nr CCL-cervix adenocarcinoma. All cell lines were cultured according to provider’s instructions.
Real-time PCR
RNA was isolated and reverse transcribed as described previously [33]. 200 ng of RNA was used for reverse transcription. Real-time PCR was performed as described previously [33] using primers specific for transcripts of THRB (THRB-RT-F: GAACAGTCGTCGCCACATC and THRB-RT-R: GCTCGTCCTTGTCTAAGTAAC) and DNMT1 (DNMT1-RT-F: TTATCCGAGGAGGGCTAC and DNMT1-RT-R: GGCTTCACTTCTTGCTTG). Reverse transcription and real-time PCR of microRNAs was performed as described previously [34] using Taq-man probes specific for hsa-miR-155 (Life Technologies, Carlsbad, CA, USA, cat. no. 002623) and hsa-miR-425 (Life Technologies, Carlsbad, CA, USA, cat. no. 001516). Relative quantification of each expressed miRNA was calculated using the standard 2−ΔCt method. Association between the microRNA and THRB expression in tissue samples was calculated using the R environment [35] and the multiple correlation analysis according to the formula:
5-aza-2′deoxycytidine Treatment
UOK-171 cells were seeded onto 12-well plates (Corning, NY, USA) in amount of 5×10∧4 per well. The cells were cultured 24 hours in standard conditions, followed by addition of 100 mM or 10 mM 5-aza-2′deoxycytidine (5-aza-dC) (Sigma-Aldrich, Saint Louis, MO, USA). After 24 hours the cells were washed with PBS and used for RNA isolation.
Analysis of THRB Methylation
For methylation analysis, we used the sequence of previously published THRB promoter [36], GenBank Acc. no. S37458.1), updated according to the sequence of chromosome 3 genomic contig (GenBank Acc. no. NT_022517.18). Due to disconcordance between the two deposited sequences, we cloned and directly sequenced DNA region encompassing THRB promoter.
To achieve this, genomic DNA was isolated from a non-tumorous kidney sample as described previously [37]. The region containing THRB promoter was amplified using primers THRB-HindIII-PromU (CGTAAAGCTTCATATGGGTAACACTGAGGGCATAGC) and THRB-HindIII-PromL (CGTAAAGCTTACTCCTTGGGCGAAGAGAGGC) and Hot Start Perpetual OptiTaq (EURx, Gdansk, Poland), under the following conditions: 95°C-10 min., 5 cycles: [96°C-35 s, 62°C-30 s, 73°C-80 s], 40 cycles: [96°C-40 s, 59°C-30 s, 72°C-80 s], final elongation: 72°C, 10 min. The obtained PCR product was cloned into pGEM-T vector (Promega, Madison, WI, USA) and sequenced. The obtained sequence was deposited in GenBank (Acc no. KF669869).
To identify CpG islands, we employed in silico analysis using CpG Plot [38] and CpG Island Searcher [39], CpG-rich region was detected in a sequence encompassing the promoter and 5′ UTR (exon 1) of TRβ transcript.
To analyze methylation of THRB CpG islands, 800 ng of genomic DNA was bisulphite-converted using Imprint DNA Modification Kit (Sigma-Aldrich, Saint Louis, MO, USA) according to manufecturer’s instructions. To predict the nucleotide changes resulting from bisulphite conversion, THRB sequence was in silico converted using Snake charmer (http://insilico.ehu.es/restriction/two_seq/snake_charmer.html) and used for design of primers.
To perform Bisulfite-Sequencing PCR (BSP), 100 ng of bisulphite-converted DNA was used, together with Perpetual Taq HOT START (EURx, Gdansk, Poland) polymerase and primers given in Table S2, under the following conditions: initial denaturation: 95°C, 10 min, 38 cycles: [denaturation: 95°C-15 s, annealing: 56°C or 57°C-15 s, elongation: 68°C-30 s], final elongation: 75°C, 15 min, final incubation: 4°C. The annealing temperature varied from 56°C to 67°C, depending on the used primers (Table S2). The obtained PCR products were directly sequenced by a commercial service (Genomed, Warsaw, Poland).
Methylation-specific PCR (MSP-PCR) was performed as described previously [26]. The sequences of primers are given in Table S3.
SNaPshot analysis was performed using SNaPshot Multiplex Kit (Applied Biosystems, Foster City, CA, USA) and primers given in Table S3 according to manufacturer’s instructions. The products were analyzed by a commercial service (Biote21, Krakow, Poland).
The specificity of primers used for methylation analysis was analyzed using MethPrimer [40].
Analysis of microRNAs Targeting THRB Transcript
microRNAs potentially binding THRB 3′UTR were predicted using miRecords and DIANA–mirPath [41],[42] (Table S4). Next, PubMed was searched to find information on the expression of the predicted microRNAs in ccRCC. Only microRNAs with both reported increased expression in renal tumors and potentially targeting the THRB transcript were taken for further analysis.
To analyze the effect of microRNAs on indigenous THRB expression, UOK171 cells were seeded onto 12-well plates (Corning, NY, USA) using 5×10∧4 cells per well and transfected 24 hours later. Transfection mix was prepared as follows: 1) 2 µl Lipofectamine 2000 (Invitrogen/Life Technologies, Carlsbad, CA, USA) were mixed with 125 µl Opti-MEM (Gibco/Life Technologies, Carlsbad, CA, USA), and 2) 50 pmol microRNA mimics or inhibitors (Ambion/Life Technologies, Foster City, CA, USA, Table S5) were mixed with 125 µl Opti-MEM. After 10 minute incubation at room temperature the solutions were combined, incubated for additional 10 min, and added to the wells. After 6 hours transfection mixes were replaced with full medium. After 24 hours, the cells were harvested for RNA isolation.
For luciferase reporter assays, HeLa were chosen due to low endogenous expression of microRNAs (Figure S1). The cells were seeded onto 12-well plates (Corning, NY, USA) using 5×10∧4 cells per well and transfected 24 hours later. Transfection mix was prepared as follows: 1) 4 µl PEI 1µg/µl (Polyethylenimine, Polysciences, Warrington, PA, USA) were mixed with 125 µl Opti-Mem, and 2) 1 µg pGL3-luc-3′UTR-THRB plasmid (Jazdzewski et al., 2011) was mixed with 125 µl Opti-Mem. For reference plasmid, the following mixes were prepared: 1) 2 µl PEI with 62.5 µl Opti-Mem, and 2) 100 ng pRL-TK plasmid (Promega, Madison, WI, USA) with 62.5 µl Opti-Mem. After 20 minute incubation at room temperature the solutions were combined, incubated for additional 20 min and added to culture wells containing 700 µl medium without FBS. After 6 hours, the medium was replaced with transfection mixes containing microRNA mimics- prepared as described above. Information on microRNA mimics and inhibitors is given in Table S5. After 6 hours, transfection mixes were replaced with full medium; cells were lysed after 24 hours and analyzed in Dual Luciferase Reporter Assay (Promega, Madison, WI, USA) using Synergy 2 luminometer (BioTek, Winooski, VT, USA).
Results
The Expression of THRB and DNMT1 is Changed in Renal Cancer
Real-time PCR analysis performed on ccRCC and paired non-tumorous control samples revealed significant changes in THRB mRNA expression (Figure 1). The expression of THRB was decreased by 56% (p<0.0001) in tumor when compared with control samples. This decrease of expression was found in 32 out of the 35 (91.4%) analyzed paired samples.
A. THRB and DNMT1 mRNA expression in tissue samples: non-tumorous controls (C, n = 35) and paired ccRCC samples (T, n = 35). The data are shown as percent of control (C). Real-time PCR for each sample was performed in triplicates. Statistical analysis was performed using Wilcoxon matched pair test. **** p<0.0001. B. The expression of THRB and DNMT1 mRNA in cell lines: HK2: proximal tubular cell line derived from normal kidney, UOK171: ccRCC cell line derived from IV stage tumor, spontaneously immortalized. The plots show results of five independent biological experiments, measured in triplicates. Statistical analysis was performed using t-test. ** p<0.01.
Real-time PCR analysis performed on cDNA from kidney cell lines confirmed decreased expression of THRB in renal cancer. THRB expression was 63% (p = 0.0036) lower in ccRCC UOK171 cell line when compared with non-cancerous HK2 cells (Figure 1).
Previous studies analyzing the expression of DNA methyltransferase 1 (DNMT1) in renal cancer brought conflicting results [43]–[45]. Therefore, we decided to analyze DNMT1 expression in our study. DNMT1 mRNA levels were increased in 22 out of 35 analyzed ccRCC samples when compared with paired controls (Figure 1). The mean DNMT1 expression in tumor tissues was increased by 38% when compared with paired control tissue samples (p = 0.0098). In cell lines, DNMT1 expression was increased by 27% in UOK171 when compared with non-tumorous HK2 cell line (p = 0.0059).
These results suggested that increased DNMT1 levels in ccRCC samples could possibly result in elevated CpG methylation leading to downregulation of THRB expression.
Analysis of THRB CpG Methylation in Renal Cancer
To check whether the decreased THRB expression could result from CpG hypermethylation, UOK171 cells were cultured in the presence or absence of DNMT1 inhibitor, 5-aza-dC. Treatment with 5-aza-dC resulted in a statistically significant (p = 0.0124), 37% increase of THRB expression when compared with control, untreated cells (Figure 2A).
A. The effect of 5-aza-2′ deoxycytidine (5-aza-dC) on THRB mRNA expression in UOK171 cell line. The results are shown as percent of control (the cells cultured without 5-aza-dC supplementation). The plot shows results of three independent biological experiments, measured in triplicates. Statistical analysis was performed using t-test. *p<0.05. B. The sequence of THRB gene with the promoter region (GenBank Acc.no. KF669869). CpG island, encompassing the promoter and 1st exon (region −873 to +355) is shadowed blue. TATA box is bolded. CpG dinucleotides are shadowed gray. Lower case letters indicate 1st exon of the TRβ transcript. C. Representative electrophoretic analysis of MSP-PCR. U: PCR products obtained with primers specific to unmethylated sequence; M: PCR products obtained with primers specific to methylated sequence. C: control samples, T: tumor samples.
To determine whether restoration of THRB expression in renal cancer is directly regulated by promoter methylation, THRB sequences were analyzed using CpG Plot and CpG Island Searcher. Both programs predicted a CpG-rich region within the promoter and 5′UTR of THRB gene, encompassing nucleotides −873 to +355 (Figure 2B). Next, the predicted region was analyzed using bisulphite sequencing PCR (BSP). To this end, DNA samples from 15 pairs of ccRCC tissues and corresponding non-tumorous control kidney samples were bisulphite-converted and directly sequenced. Efficient bisulphite conversion was confirmed in PCR reactions with the use of primers designed to target native THRB promoter and not amplify bisulphite-converted DNA. Control amplification reactions performed on bisulphite converted DNA produced no amplification product (control primer sequences are given in Table S2). This analysis revealed lack of differences in methylation patterns between paired control and tumor samples from the same patient (Figure S2). To verify the results of BSP, the bisulphite converted DNA was additionally analyzed using SNaPshot and MSP-PCR, a highly sensitive method allowing for detection of 0.1% methylated alleles of a given CpG island locus [46] (Figure 2C). These analyses confirmed that there were no changes in CpG methylation patterns in ccRCC samples when compared with paired control samples. Interestingly, MSP-PCR results suggested that methylation patterns varied between different patients, being rather patient-specific than influenced by the disease status.
In conclusion, these results suggested that although changes in DNMT1 activity may contribute to alterations of THRB expression in cell cultures, the disturbed expression of THRB in ccRCC tissue samples is rather not a result of aberrant, tumor-specific hypermethylation of the gene’s promoter.
microRNAs miR-155 and miR-425 Directly Target THRB 3′UTR
To predict microRNAs potentially influencing THRB expression in renal cancer, two-step analysis was performed. Firstly, bioinformatic analysis using miRecords (a resource that combines 11 bioinformatic programs) and Diana-mirPath predicted nearly 600 microRNAs potentially targeting THRB 3′UTR. However, based on the assumption that silencing efficiency is the highest for microRNAs that lie within the closest proximity to either end of the 3′UTR [47] we specifically focused on microRNAs that complied with that rule. Subsequent search in PubMed allowed for selection of microRNAs that were reported as overexpressed in renal cancer [18]–[21]. Using the two-step approach, four microRNAs (miR-155, miR-425, miR-592, and miR-599) were selected for further analysis.
In order to find whether the chosen microRNAs indeed bind the 3′UTR of THRB, HeLa cells were transfected using synthetic microRNAs and a reporter construct pGL3-luc-3′UTR-THRB (Figure 3). HeLa cells were chosen for this analysis due to low endogenous expression of THRB and analyzed microRNAs. Transfection of cells with pre-miR-155 or pre-miR-425 resulted in 32% or 17% decrease of luciferase activity, respectively, when compared with control miRNA mimics (p<0.001). No changes in luciferase activity were observed in cells transfected with pre-miR-592 or pre-miR-599. Pre-miR-221, previously reported to target THRB 3′UTR [25] was used as positive control and confirmed the efficiency of miRNA-mediated silencing (Figure 3).
A. The effect of microRNA mimics on the activity of luciferase expressed from a reporter plasmid pGL3-luc-3′UTR-THRB. HeLa cells were transfected using the reporter plasmid and respective microRNA mimics: pre-miR-155, pre-miR-425, pre-miR-599, or pre-miR-221. The results are shown as percentage of control (Ctrl, cells transfected with non-targeting pre-miR). The plot shows results of three independent biological experiments, measured in triplicates. The relative activity of Firefly luciferase was normalized to Renilla luciferase activity. Statistical analysis was performed using ANOVA followed by Dunnett’s multiple comparison test. *** p<0.001. B. The effect of miR-155 and miR-425 on endogenous THRB expression in renal cancer cells. UOK171 cells were transfected using respective microRNA mimics (pre-miRs) or inhibitors (anti-miRs), and THRB mRNA expression was analyzed using real-time PCR. The results are shown as percent of control (the cells transfected with non-targeting pre-miR). The plots shows results of three independent biological experiments measured in triplicates. Statistical analysis was performed using ANOVA followed by Dunnett’s multiple comparison test. *p<0.05, *** p<0.001. C. The expression of miR-155 and miR-425 in nontumorous controls (C, n = 35) and paired ccRCC tissue samples (T, n = 35). The results are shown as percent of control. Real-time PCR for each sample was performed in triplicates. Statistical analysis was performed using Wilcoxon matched pair test. *p<0.05, **p<0.01.
miR-155 and miR-452 Influence THRB Expression in Renal Cancer Cells
To subsequently analyze the effect of miR-155 and miR-425 on endogenous THRB expression, renal cancer-derived UOK171 cell line was transfected using either microRNA mimics (pre-miR-155 or pre-miR-425) or microRNA inhibitors (anti-miR-155 or anti-miR-425) (Figure 3). Transfection resulted in a statistically significant decrease in THRB expression by 22% for pre-miR-155, p<0.05, and by 64% for pre-miR-425, p<0.001, when compared with control microRNA (Figure 3). Transfection of cell lines with microRNA inhibitors did not result in a statistically significant change of THRB expression, although a trend for increased expression was seen (Figure 3). This observation can result from the fact that levels of miR-155 and miR-425 are highly elevated in UOK171 cell line, thus inspite of numerous concentrations tested, the amount of microRNA inhibitors could be insufficient for efficient blocking of microRNA activity.
These studies suggested that miR-155 and miR-425 directly influence THRB expression in renal cancer cells.
The Expression of miR-155 and miR-425 is Elevated in Renal Cancer and Correlates with Decreased Expression of THRB
In order to check whether tumor-specific changes of miR-155 and miR-425 could affect THRB expression in ccRCC tumors, the expression of both microRNAs was analyzed in 25 ccRCC tissue samples and 25 corresponding non-tumorous kidney samples (Figure 3). The expression of both microRNAs was statistically significantly increased (by 354%, p = 0.0072 for miR-155 and by 62%, p = 0.041 for miR-425). Simultaneously, expression of both miRs was negatively correlated with expression of THRB. Although the correlation coefficient was weak (R = 0.245, p = 0.0075) in control samples, the correlation was significantly stronger in tumors (R = 0.469, p = 0.0225). Concomitantly, there was also a significant correlation between the T/N ratio of miRNA and THRB expression (R = 0.396, p = 0.002). These results indicate that THRB expression levels depend on the miR-155 and miR-425 expression and deregulation of miRs is linked with lowered THRB levels in ccRCC (Figure S3).
Discussion
In this study, using three different methods, we showed that in the analyzed samples of renal cancer aberrant THRB expression does not result from tumor-specific changes in DNA methylation of THRB promoter region. Rather, the expression of THRB in ccRCC is affected by microRNAs, miR-155 and miR-425 that directly bind to THRB 3′UTR, as suggested by the observation that elevated expression of these microRNAs in ccRCC is accompanied by downregulation of THRB.
Hypermethylation of THRB promoter was reported in several neoplasms, including breast, thyroid and lung cancers and leukemia [25]–[30]. Interestingly, the studies showed that methylation rates of THRB were highly variable among patients, since hypermethylated THRB was detected in 25–80% of analyzed tumors. Furthermore, only half of the studies on THRB methylation was performed on both tumor and paired non-tumorous control tissues. Iwasaki et al. examined 116 non-small cell lung cancers, including 6 samples for which paired adjacent control tissues were analyzed. Out of the 116 tumor samples, 54 revealed hypermethylation of THRB. When 6 paired tumor versus control samples were compared, THRB was methylated in all tumors and unmethylated in all non-tumorous control samples. In two other studies analyzing breast cancer tumors and cell lines [25],[29] hypermethylated THRB was detected in 80%-100% of tumors and c.a. 37% of paired samples. Moreover, the degree and pattern of methylation varied between samples derived from different patients. Similar interindividual variability was also observed in a study performed on thyroid tumors [30], in which paired normal tissue samples were not available and therefore not analyzed. Our study also revealed that THRB methylation significantly varied between tissue samples derived from different patients. These results suggest that lack of paired non-tumorous control samples may lead to a bias and result in classification of a patient-specific variations in THRB methylation as changes specific for cancer tissue.
Although THRB expression did not seem to be influenced by DNA methylation in our renal cancer samples, treatment of UOK171 cell line with 5-aza-2′ deoxycytidine resulted in increased expression of THRB gene. This suggests that THRB expression may be indirectly influenced by DNMT1 activity. According to previous studies, DNA methylation can indirectly affect transcription of a target gene in several mechanisms [48] including gene regulation by a reactivated transcription factor or signal transduction pathway that is affected by DNA methylation. Which of these mechanisms can affect THRB expression in renal cancer is currently unknown. Furthermore, we can not exclude the possibility that in other ccRCC patients and in the other ccRCC-derived cell lines, the THRB promoter methylation could influence THRB expression. This, hopefully, will be verified by future studies on ccRCC.
Another epigenetic mechanism that influences the expression of target genes is microRNA-dependent regulation of expression [49]. A single gene can be targeted by multiple microRNAs and single microRNA can regulate numerous target genes. In consequence, deregulation of miRNA expression that is associated with numerous cancers leads to a severe disruption of the cellular proteome and transcriptome. In this study we identified two microRNAs whose direct binding with THRB 3′UTR was confirmed in luciferase assay performed in HeLa cell line, used due to relatively low expression levels of both analyzed miRNAs. Analysis of expression in tissue samples revealed that upregulation of miR-155 and miR-425 in clear cell renal cell carcinoma is accompanied by downregulation of THRB. Insufficient specificity of the commercially offered anti-TRβ antibodies and extremely low expression levels of the TRβ protein did not allow for showing the miR-mediated effect on the TRβ protein level. It was however proven both in luciferase assay and in vivo study, in which the miRs downregulated indigenous THRB transcription.
A number of microRNAs regulating THRB expression were already described [9],[24],[50] and Jazdzewski et al. showed that microRNAs targeting THRB are strongly upregulated in papillary thyroid carcinoma, contributing to silencing of THRB in tumor. Importantly, the expression of microRNAs is highly specific for a tissue type, or characteristic for a specific disease state. Thus, microRNA signatures of ccRCC significantly differ from those of other cancers, including PTC, and this explains the fact why the obtained results were not overlapping. Since miR-155 and miR-425 analyzed in our study directly target 3′UTR of THRB and are overexpressed in ccRCC, it is possible that they contribute to decreased THRB expression in this cancer. The coefficient of multiple correlation, which is a measure of how well a given variable can be predicted using a linear function of a set of other variables, showed that THRB expression in tissue samples depends on the levels of both microRNAs. This regulation, together with previously reported altered splicing and mutations in THRB could directly lead to deregulation of thyroid receptor beta in renal cancer.
Interestingly, it was recently shown that DNA methylation and microRNA-dependent regulation are the mechanisms that complementarily regulate gene expression [51]. The study revealed a negative correlation between the level of promoter region methylation and a number of miRNA target sites in 3′UTR of a gene. Apparently, THRB belongs to the group of genes with long 3′UTRs, indeed targeted by numerous microRNAs [9],[24],[50]. This supports the notion on the important role of microRNA-dependent regulation of THRB, and on the possible role of deregulation of this mechanism in neoplastic process.
Both microRNAs analyzed in our study have well established role in tumorigenesis. For instance, aberrant expression of miR-425 was reported in hyperdiploid multiple myeloma, glioblastoma, breast, thyroid and prostate cancers [52]–[56]. Those aberrations in miR-425 expression can directly influence tumorigenesis, e.g. by affecting the expression of oncogenes [52]. Finally, in accordance with our study, upregulated expression of miR-425 was reported in renal cancer [57] and was proposed as one of markers differentiating between normal kidney and renal tumors [20].
MiR-155 is a well known oncogene which is overexpressed in a broad panel of neoplastic lesions, including cancers of breast, colon, cervix, pancreas, lung, and kidney, as well as in lymphomas and leukemias [19],[58],[59]. Most interestingly, it was recently shown that miR-155 targets VHL tumor suppressor and promotes angiogenesis in breast cancer [60]. VHL plays a pivotal role in ccRCC tumorigenesis and is mutated or silenced in more than half of sporadic clear cell renal cell carcinomas [61]. Intriguingly, VHL itself regulates the expression of miR-155 in ccRCC [62], suggesting a tumor suppressor-oncogene feedback regulation between VHL and miR-155. Moreover, the role of miR-155 in silencing of a tumor suppressor, revealed in this study, additionally supports the notion on the oncogenic role of this miRNA.
In conclusion, we identified two novel microRNAs, miR-155 and miR-425, that target THRB transcript and downregulate its expression in renal cancer. In contrast, CpG methylation is rather not the main mechanism directly contributing to deregulated THRB expression in renal cancer. The possible indirect effects of changes in DNMT1 activity on THRB expression need further evaluation.
Supporting Information
Figure S1.
Analysis of expression of microRNAs miR-155 and miR-425 in HeLa and UOK171 cells. The results show analysis of expression from cells cultured in three 25 cm2 bottles, normalized to U6 snRNA. Real-time PCR for each sample was performed in triplicates. Statistical analysis was performed using unpaired t test. **p<0.01, *p<0.001.
https://doi.org/10.1371/journal.pone.0097624.s001
(TIF)
Figure S2.
Analysis of THRB CpG methylation. Representative result of BSP. The chromatograms show results of sequencing performed on tumor (upper panel) and control (lower panel) samples taken from the same patient.
https://doi.org/10.1371/journal.pone.0097624.s002
(PDF)
Figure S3.
Correlation between the ratio of miR-155 and miR-425 and THRB expression in tumor vs control samples. Multiple correlation plots showing the association between the changes of microRNA and THRB levels in ccRCC samples. Lowered THRB expression in ccRCC when compared to paired control tissue is correlated with increased levels of miR-155 and miR-425.
https://doi.org/10.1371/journal.pone.0097624.s003
(TIF)
Table S1.
Patient characteristics. Information on 35 ccRCC patients included in the study: histopathological diagnosis, age at disease onset and gender.
https://doi.org/10.1371/journal.pone.0097624.s004
(DOCX)
Table S2.
Primers used for BSP of THRB promoter.
https://doi.org/10.1371/journal.pone.0097624.s005
(DOCX)
Table S3.
Primers used in MSP-PCR and SNaPshot.
https://doi.org/10.1371/journal.pone.0097624.s006
(DOCX)
Table S4.
microRNAs predicted to bind 3′UTR of THRB based on miRecords database.
https://doi.org/10.1371/journal.pone.0097624.s007
(XLS)
Table S5.
microRNA mimics and inhibitors used in the study.
https://doi.org/10.1371/journal.pone.0097624.s008
(DOCX)
Acknowledgments
UOK171 was provided by UOB Tumor Cell Line Repository, led by Dr. W Marston Linehan, NIH, NCI, Bethesda, USA, patent E-033-2010/0. The authors wish to thank the researchers from the Laboratory of Genomic Medicine of the Medical University of Warsaw: Dr. Krystian Jazdzewski for the kind gift of pGL3-luc-3′UTR-THRB plasmid and Anna Kubiak for the miR-THRB correlation analysis.
Author Contributions
Conceived and designed the experiments: AW APW HK BR JB PP AM AN. Performed the experiments: AW APW HK BR JB PP AM. Analyzed the data: AW APW HK JB PP AM. Contributed reagents/materials/ analysis tools: AW AN. Wrote the paper: AW APW HK BR JB PP AM AN.
References
- 1. Qi JS, Yuan Y, Desai-Yajnik V, Samuels HH (1999) Regulation of the mdm2 oncogene by thyroid hormone receptor. Mol Cell Biol 19: 864–872.
- 2. Dinda S, Sanchez A, Moudgil V (2002) Estrogen-like effects of thyroid hormone on the regulation of tumor suppressor proteins, p53 and retinoblastoma, in breast cancer cells. Oncogene 21: 761–768.
- 3. Puzianowska-Kuznicka M, Pietrzak M, Turowska O, Nauman A (2006) Thyroid hormones and their receptors in the regulation of cell proliferation. Acta Biochim Pol 53: 641–650.
- 4. Martínez-Iglesias O, Garcia-Silva S, Tenbaum SP, Regadera J, Larcher F, et al. (2009) Thyroid hormone receptor beta1 acts as a potent suppressor of tumor invasiveness and metastasis. Cancer Res 69: 501–509.
- 5. Kim WG, Cheng SY (2013) Thyroid hormone receptors and cancer. Biochim Biophys Acta 1830: 3928–3936.
- 6. Hagenkord JM, Gatalica Z, Jonasch E, Monzon FA (2011) Clinical genomics of renal epithelial tumors. Cancer Genet 204: 285–297.
- 7. Kamiya Y, Puzianowska-Kuznicka M, McPhie P, Nauman J, Cheng SY, et al. (2002) Expression of mutant thyroid hormone nuclear receptors is associated with human renal clear cell carcinoma. Carcinogenesis 23: 25–33.
- 8. Puzianowska-Kuznicka M, Nauman A, Madej A, Tanski Z, Cheng S, et al. (2000) Expression of thyroid hormone receptors is disturbed in human renal clear cell carcinoma. Cancer Lett 155: 145–152.
- 9. Master A, Wójcicka A, Piekiełko-Witkowska A, Bogusławska J, Popławski P, et al. (2010) Untranslated regions of thyroid hormone receptor beta 1 mRNA are impaired in human clear cell renal cell carcinoma. Biochim Biophys Acta 1802: 995–1005.
- 10. Arai E, Kanai Y (2011) Genetic and epigenetic alterations during renal carcinogenesis. Int J Clin Exp Pathol 4: 58–73.
- 11. Okano M, Bell DW, Haber DA, Li E (2009) DNA Methyltransferases DNMT3a and DNMT3b are Essential for De Novo Methylation and Mammalian Development. Cell 99: 247–257.
- 12. McCabe MT, Davis JN, Day ML (2005) Regulation of DNA methyltransferase 1 by the pRb/E2F1 pathway. Cancer Res 65: 3624–3632.
- 13. Ehrlich M (2002) DNA methylation in cancer: too much, but also too little. Oncogene 21: 5400–5413.
- 14. Tryndyak VP, Kovalchuk O, Pogribny IP (2006) Loss of DNA methylation and histone H4 lysine 20 trimethylation in human breast cancer cells is associated with aberrant expression of DNA methyltransferase 1, Suv4–20h2 histone methyltransferase and methyl-binding proteins. Cancer Biol Ther 5: 65–70.
- 15. Esteller M, Garcia-Foncillas J, Andion E, Goodman SN, Hidalgo OF, et al. (2000) Inactivation of the DNA-repair gene MGMT and the clinical response of gliomas to alkylating agents. N Engl J Med 343: 1350–1354.
- 16. Koul H, Huh JS, Rove KO, Crompton L, Koul S, et al. (2011) Molecular aspects of renal cell carcinoma: a review. Am J Cancer Res 1: 240–254.
- 17. Filipowicz W, Bhattacharyya SN, Sonenberg N (2008) Mechanisms of post-transcriptional regulation by microRNAs: are the answers in sight? Nat Rev Genet 9: 102–114.
- 18.
Liu H, Brannon AR, Reddy AR, Alexe G, Seiler MW, et al.. (2010) Identifying mRNA targets of microRNA dysregulated in cancer: with application to clear cell Renal Cell Carcinoma. BMC Syst Biol 27 4: 51. Available: http://www.biomedcentral.com/1752-0509/4/51. Accessed: 06 Aug 2013.
- 19. Juan D, Alexe G, Antes T, Liu H, Madabhushi A, et al. (2010 ) Identification of a microRNA panel for clear-cell kidney cancer. Urology 75: 835–841.
- 20. Youssef YM, White NM, Grigull J, Krizova A, Samy C, et al. (2011) Accurate molecular classification of kidney cancer subtypes using microRNA signature. Eur Urol. 59: 721–730.
- 21.
Osanto S, Qin Y, Buermans HP, Berkers J, Lerut E, et al.. (2012) Genome-Wide MicroRNA Expression Analysis of Clear Cell Renal Cell Carcinoma by Next Generation Deep Sequencing. PLoS One; 7(6): e38298. Available: http://www.plosone.org/article/info%3Adoi%2F10.1371%2Fjournal.pone.0038298. Accessed: 06 Aug 2013.
- 22. Hatziapostolou M, Iliopoulos D (2011) Epigenetic aberrations during oncogenesis. Cell Mol Life Sci 68: 1681–1702.
- 23. Bader AG, Brown D, Stoudemire J, Lammers P (2011) Developing therapeutic microRNAs for cancer. Gene Ther 18: 1121–1126.
- 24. Jazdzewski K, Boguslawska J, Jendrzejewski J, Liyanarachchi S, Pachucki J, et al. (2010) Thyroid hormone receptor beta (THRB) is a major target gene for microRNAs deregulated in papillary thyroid carcinoma (PTC). J Clin Endocrinol Metab 96: E546–53.
- 25. Li Z, Meng ZH, Chandrasekaran R, Kuo WL, Collins CC, et al. (2002) Biallelic inactivation of the thyroid hormone receptor beta1 gene in early stage breast cancer. Cancer Res 62: 1939–1943.
- 26. Joseph B, Ji M, Liu D, Hou P, Xing M (2007) Lack of mutations in the thyroid hormone receptor (TR) alpha and beta genes but frequent hypermethylation of the TRbeta gene in differentiated thyroid tumors. J Clin Endocrinol Metab 92: 4766–4770.
- 27. Dunwell TL, Hesson LB, Pavlova T, Zabarovska V, Kashuba V, et al. (2009) Epigenetic analysis of childhood acute lymphoblastic leukemia. Epigenetics 4: 185–193.
- 28. Iwasaki Y, Sunaga N, Tomizawa Y, Imai H, Iijima H, et al. (2010 ) Epigenetic inactivation of the thyroid hormone receptor beta1 gene at 3p24.2 in lung cancer. Ann Surg Oncol 17: 2222–2228.
- 29. Ling Y, Xu X, Hao J, Ling X, Du X, et al. (2010) Aberrant methylation of the THRB gene in tissue and plasma of breast cancer patients. Cancer Genet Cytogenet 196: 140–145.
- 30. Kim WG, Zhu X, Kim DW, Zhang L, Kebebew E, et al. (2013) Reactivation of the silenced thyroid hormone receptor β gene expression delays thyroid tumor progression. Endocrinology 154: 25–35.
- 31. Larkin J, Goh XY, Vetter M, Pickering L, Swanton C (2012 ) Epigenetic regulation in RCC: opportunities for therapeutic intervention? Nat Rev Urol 9: 147–55.
- 32. Lopez-Beltran A, Scarpelli M, Montironi R, Kirkali Z (2006) 2004 WHO classification of the renal tumors of the adults. Eur Urol 49: 798–805.
- 33.
Piekielko-Witkowska A, Wiszomirska H, Wojcicka A, Poplawski P, Boguslawska J, et al.. (2010) Disturbed expression of splicing factors in renal cancer affects alternative splicing of apoptosis regulators, oncogenes, and tumor suppressors. PLoS One 5: e13690. Available: http://www.plosone.org/article/info%3Adoi%2F10.1371%2Fjournal.pone.0013690. Accessed: 06 Aug 2013.
- 34.
Boguslawska J, Wojcicka A, Piekielko-Witkowska A, Master A, Nauman A (2011) MiR-224 targets the 3′UTR of type 1 5′-iodothyronine deiodinase possibly contributing to tissue hypothyroidism in renal cancer. PLoS One 6: e24541. Available: http://www.plosone.org/article/info%3Adoi%2F10.1371%2Fjournal.pone.0024541. Accessed: 06 Aug 2013.
- 35.
R Development Core Team (2008) R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. ISBN 3–900051–07–0.
- 36. Sakurai A, Miyamoto T, DeGroot LJ (1992) Cloning and characterization of the human thyroid hormone receptor beta 1 gene promoter. Biochem Biophys Res Commun 185: 78–84.
- 37. Laird PW, Zijderveld A, Linders K, Rudnicki MA, Jaenisch R, et al. (1991) Simplified mammalian DNA isolation procedure. Nucleic Acids Res 19: 4293.
- 38. Rice P, Longden I, Bleasby A (2000) EMBOSS: the European molecular biology open software suite. Trends Genet 16: 276–277.
- 39. Takai D, Jones PA (2003) The CpG island searcher: a new WWW resource. In Silico Biol 3: 235–240.
- 40. Li LC, Dahiya R (2002) MethPrimer: designing primers for methylation PCRs. Bioinformatics 18: 1427–1431.
- 41. Xiao F, Zuo Z, Cai G, Kang S, Gao X, et al. (2009) miRecords: an integrated resource for microRNA-target interactions. Nucleic Acids Res 37: D105–D110.
- 42. Papadopoulos GL, Alexiou P, Maragkakis M, Reczko M, Hatzigeorgiou AG (2009) DIANA-mirPath: Integrating human and mouse microRNAs in pathways. Bioinformatics 25: 1991–1993.
- 43. Robertson KD, Uzvolgyi E, Liang G, Talmadge C, Sumegi J, et al. (1999) The human DNA methyltransferases (DNMTs) 1, 3a and 3b: coordinate mRNA expression in normal tissues and overexpression in tumors. Nucleic Acids Res 27: 2291–2298.
- 44. Kimura F, Seifert HH, Florl AR, Santourlidis S, Steinhoff C, et al. (2003) Decrease of DNA methyltransferase 1 expression relative to cell proliferation in transitional cell carcinoma. Int J Cancer 104: 568–578.
- 45. Arai E, Kanai Y, Ushijima S, Fujimoto H, Mukai K, et al. (2006) Regional DNA hypermethylation and DNA methyltransferase (DNMT) 1 protein overexpression in both renal tumors and corresponding nontumorous renal tissues. Int J Cancer 119: 288–296.
- 46. Herman JG, Graff JR, Myöhänen S, Nelkin BD, Baylin SB (1996) Methylation-specific PCR: a novel PCR assay for methylation status of CpG islands. Proc Natl Acad Sci U S A 93: 9821–9826.
- 47.
Virtue A, Wang H, Yang XF (2012) MicroRNAs and toll-like receptor/interleukin-1 receptor signaling. J Hematol Oncol 5: 66. Available: http://www.jhoonline.org/content/5/1/66. Accessed: 06 Aug 2013.
- 48. Karpf AR, Jones DA (2002) Reactivating the expression of methylation silenced genes in human cancer. Oncogene 21: 5496–5503.
- 49. Sharma S, Kelly TK, Jones PA (2010) Epigenetics in cancer. Carcinogenesis 31: 27–36.
- 50. Nishi H, Ono K, Horie T, Nagao K, Kinoshita M, et al. (2011) MicroRNA-27a regulates beta cardiac myosin heavy chain gene expression by targeting thyroid hormone receptor beta1 in neonatal rat ventricular myocytes. Mol Cell Biol 31: 744–755.
- 51.
Su Z, Xia J, Zhao Z (2011) Functional complementation between transcriptional methylation regulation and post-transcriptional microRNA regulation in the human genome. BMC Genomics 12: S15. Available: http://www.biomedcentral.com/1471-2164/12/S5/S15. Accessed: 06 Aug 2013.
- 52. Rio-Machin A, Ferreira BI, Henry T, Gómez-López G, Agirre X, et al. (2013) Downregulation of specific miRNAs in hyperdiploid multiple myeloma mimics the oncogenic effect of IgH translocations occurring in the non-hyperdiploid subtype. Leukemia 27: 925–931.
- 53. Gal H, Pandi G, Kanner AA, Ram Z, Lithwick-Yanai G, et al. (2008) MIR-451 and imatinib mesylate inhibit tumor growth of glioblastoma stem cells. Biochem Biophys Res Commun 376: 86–90.
- 54. Van der Auwera I, Limame R, van Dam P, Vermeulen PB, Dirix LY, et al. (2010) Integrated miRNA and mRNA expression profiling of the inflammatory breast cancer subtype. Br J Cancer 103: 532–541.
- 55. Jacques C, Guillotin D, Fontaine JF, Franc B, Mirebeau-Prunier D, et al. (2013) DNA Microarray and miRNA Analyses Reinforce the Classification of Follicular Thyroid Tumors. J Clin Endocrinol Metab 98: E981–989.
- 56.
Watahiki A, Wang Y, Morris J, Dennis K, O’Dwyer HM, et al.. (2011) MicroRNAs associated with metastatic prostate cancer. PLoS One 6: e24950. Available http://www.plosone.org/article/info%3Adoi%2F10.1371%2Fjournal.pone.0024950. Accessed 06 August 2013.
- 57. White NM, Bao TT, Grigull J, Youssef YM, Girgis A, et al. (2011) miRNA profiling for clear cell renal cell carcinoma: biomarker discovery and identification of potential controls and consequences of miRNA dysregulation. J Urol 186: 1077–1083.
- 58. Faraoni I, Antonetti FR, Cardone J, Bonmassar E (2009) miR-155 gene: a typical multifunctional microRNA. Biochim Biophys Acta 1792: 497–505.
- 59. Jung M, Mollenkopf HJ, Grimm C, Wagner I, Albrecht M, et al. (2009) MicroRNA profiling of clear cell renal cell cancer identifies a robust signature to define renal malignancy. J Cell Mol Med 13: 3918–3928.
- 60. Kong W, He L, Richards EJ, Challa S, Xu CX, et al. (2014) Upregulation of miRNA-155 promotes tumour angiogenesis by targeting VHL and is associated with poor prognosis and triple-negative breast cancer. Oncogene 2014 33: 679–689.
- 61. Kaelin WG Jr (2004) The von Hippel-Lindau tumor suppressor gene and kidney cancer. Clin Cancer Res 10: 6290S–6295S.
- 62.
Neal CS, Michael MZ, Rawlings LH, Van der Hoek MB, Gleadle JM (2010) The VHL-dependent regulation of microRNAs in renal cancer. BMC Med 8: 64. Available: http://www.biomedcentral.com/1741-7015/8/64. Accessed: 06 Aug 2013.