Liver X Receptor β Related to Tumor Progression and Ribosome Gene Expression in Papillary Thyroid Cancer
Article information
Abstract
Background
Intracellular lipid deposition has been reported in thyroid glands in obese animal and human. To understand the regulatory mechanism of lipid metabolism in thyroid cancer, we investigated the expression status of liver X receptor (LXR) and analyzed its clinicopathological characteristics and molecular biological features.
Methods
Expression status of LXR and its transcriptional targets in human cancers were analyzed using The Cancer Genome Atlas (TCGA). The gene-sets related to high LXRβ expression was investigated by gene set enrichment analysis (GSEA) using Kyoto Encyclopedia of Genes and Genomes (KEGG) signaling pathways and gene ontology biologic process. Quantitative reverse transcription polymerase chain reaction was performed in thyroid cancer samples using our validation cohort.
Results
In contrast to low expression of LXRα, LXRβ was highly expressed in thyroid cancer compared to the other types of human cancers. High LXRβ expression was correlated with the expression of LXRβ transcriptional targets genes, such as apolipoprotein C1 (APOC1), APOC2, apolipoprotein E (APOE), ATP binding cassette subfamily G member 8 (ABCG8), sterol regulatory elementbinding protein 1c (SREBP1c), and SPOT14. Furthermore, High LXRβ expression group indicated poor clinicopathological characteristics and aggressive molecular biological features independently from the drive mutation status. Mechanistically, high LXRβ expression was coordinately related to ribosome-related gene sets.
Conclusion
The mechanistic link between LXRβ and ribosomal activity will be addressed to develop new diagnostic and therapeutic targets in thyroid cancers.
INTRODUCTION
Thyroid cancer is the most common endocrine malignancy [1,2]. The incidence of thyroid cancer has strikingly increased in the last couple of decades, especially in Korea [3]. The cause of this increase might be the early detection of differentiated thyroid cancer, mainly papillary thyroid microcarcinoma, by the introduction of high-sensitivity neck ultrasound [4,5]. Although an increase in the incidence of thyroid cancer has been postulated, its underlying cause has not been well-recognized [6,7].
The most well-known etiology of thyroid cancer is radiation exposure [8,9]. However, radiation-related papillary thyroid carcinoma (PTC) especially occurred after the Chernobyl nuclear accident caused RET/PTC rearrangements, whereas in Korea most of the PTC cases report B-Raf proto-oncogene, serine/ threonine kinase (BRAF) V600E mutation as its driver mutation [10-12]. In addition, there is no strong evidence to support that Korean individuals are subjected to increased radiation exposure [13]. Recently, other causative factors, such as iodine intake and obesity, have also been suggested [14-16]. However, there is also no evidence of increased iodine intake in the Korean population. Obesity is likely to be associated with an increased incidence of thyroid cancer in Korea [17,18]. As seen in other developed countries, the prevalence of obesity and metabolic diseases in Korea has been increasing, thus leading to an increase in the incidence of obesity-related cancer [19]. In this context, some epidemiological studies have reported that obesity might be associated with thyroid cancer, but substantial evidence is lacking.
The effects of obesity on the normal thyroid gland in animal models and humans have been reported [14]. Abnormalities in thyroid function and morphology have been observed in genetically modified obese animals such as ob/ob and db/db mice; moreover, recent studies have reported that obesity causes lipid accumulation in the thyroid gland of humans and animals [20,21]. In addition, a high fat diet generated a more aggressive phenotype of thyroid cancer in Thrb(PV/PV)Pten(+/–) mice [22, 23]. In human thyroid cancer, the body mass index has been suggested as a factor related to aggressive tumor behavior [24,25]. Moreover, an oncocytic variant of PTC was reported to be associated with obesity [26]. These data from studies in mice and humans collectively suggest that obesity might affect the normal thyroid physiology and thyroid cancer behavior. However, although an association between obesity and thyroid cancer has been suggested, the role of lipid metabolism in thyroid cancer remains to be elucidated [18,23].
Liver X receptor (LXR), recently deorphanized by the discovery of oxysterol as its endogenous ligand, is a member of the nuclear receptor family [27-29]. Liver X receptors (LXRα and LXRβ) are important regulators of cholesterol and lipid metabolism [30,31]. In this study, we investigated the expression pattern of LXRα and LXRβ in human thyroid cancer to improve our current understanding of lipid metabolism in thyroid cancer.
METHODS
Public data acquisition
Subsets of normalized raw data for RNA expression (RNA reads per kilobase per million [RPKM] converted to transcripts per million [TPM]) and clinicopathological features were provided by The Cancer Genome Atlas (TCGA, https://tcga-data.nci.nih.gov/tcga/). Using these data sets, twelve types of cancer were analyzed for their mRNA expression status of LXRα and LXRβ: breast invasive carcinoma (BRCA; 1,092 cases), cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC; 304 cases), colon adenocarcinoma (COAD; 385 cases), liver hepatocellular carcinoma (LIHC; 361 cases), lung adenocarcinoma (LUAD; 516 cases), ovary serous cystadenocarcinoma (OV; 303 cases), pancreatic adenocarcinoma (PAAD; 147 cases), prostate adenocarcinoma (PRAD; 947 cases), skin cutaneous melanoma (SKCM; 94 cases), stomach adenocarcinoma (STAD; 415 cases), thyroid carcinoma (THCA; 505 cases), and uterine corpus endometrial carcinoma (UCEC; 545 cases). LXRα and LXRβ expression levels were compared in various types of human cancers and further analyzed in thyroid cancer by comparisons with paired or unpaired normal thyroid samples. In thyroid cancer, the major transcriptional targets of LXR were also analyzed to investigate the relation of LXRβ expression with the expression of its target genes. The protein expression level was analyzed by calculating the frequency of LXRα- and LXRβ-positive samples provided by Human Protein Atlas (HPA, www.proteinatlas.org) including 12 cancer types: BRCA 11 cases, CESC 12 cases, COAD 11 cases, LIHC 11 cases, LUAD eight cases, OV 11 cases, PAAD 12 cases, PRAD 11 cases, SKCM 12 cases, STAD nine cases, THCA four cases, and UCEC 12 cases. We also selected GSE83520, GSE77039, GSE41223, GSE96545 from the Gene Expression Omnibus (GEO) database as in vitro validation set in addition to the TCGA data.
Patients and tissue samples
PTC tissue samples along with matched normal tissues from the same patients were collected from patients who underwent thyroidectomy at Yonsei Cancer Center from April to October 2010 (Seoul, South Korea). All tissues were frozen in liquid nitrogen immediately after surgery and stored at –80°C before use. This research has been approved by the Institutional Review Board of Severance Medical Center (Seoul, South Korea) on 30 September 2013 (No. 4-2013-0546); the requirement for informed consent was waived because of its retrospective nature of the study.
RNA isolation and quantitative reverse-transcription polymerase chain reaction analysis
Total RNA was extracted according to the manufacturer’s protocol (Trizol reagent, Qiagen, CA, USA). Complementary DNA (cDNA) was synthesized from 1 μg of total RNA using the QuantiTect Reverse Transcription Kit (Qiagen). The primers used in this study are as follows: LXRβ, 5´-TGC AAC AAA CGC TCC TTC TC-3´ and 5´-CCG AGA AGA TGT TGA TGG CG-3´; RPL15, 5´-TAC GGC AAG CCT GTC CAT CAT G-3´ and 5´- GTA TGT GGA ATC TTC ACC AAC CC-3´; RPL11, 5´-AGA GTG GAG ACA GACTGA CGC G-3´ and 5´-CGG ATG CCA AAG GAT CTG ACA G-3´; RPL19, 5´-TCA CAG CCT GTA CCT GAA GGT G-3´ and 5´-CGT GCT TCC TTG GTC TTA GAC C-3´; RPS30, 5´-CGT GCT GGA AAA GTG AGA GGT C-3´ and 5´-TGC CAA AGG TGG GCA CAA CGT T-3´; glyceraldehyde-3-phosphate dehydrogenase (GAPDH), 5´-ACC ACA GTC CAT GCC ATC AC-3´ and 5´-TCC ACC ACC CTG TTG CTG TA-3´. Reverse transcription-polymerase chain reaction (RT-PCR) was performed on cDNA using the Solg 2X PCR Smart mix Kit (Solgent, Daejeon, Korea). PCR products were separated on 1.0% agarose gel in 1X TAE buffer and RedSafe (Intron, Seongnam, Korea) was used for visualization of the PCR band. The image of the gel was confirmed by MINIBIS PRO (DNR, Jerusalem, Israel). All experiments were repeated three times, and each experiment was performed in triplicate.
Gene set enrichment analysis
Gene set enrichment analysis (GSEA) is software that analyzes microarray and other data with different expressions based on biological information. Enrichment score (ES) values can be calculated using weighted Kolmogorov-Smirnov statistics. The most important gene group can be identified by identifying the location (leading edge subset) where the calculated ES value represents the maximum value. Normalized enrichment score (NES) is obtained by normalizing ES values to gene set size. The expression profile was selected as the top and bottom groups by the top 25% and bottom 25%, depending on the LXRβ expression (n=126 for each group). We selected significant gene sets based on the P value ≤0.05 and the false discovery rate (FDR) q-value <0.25.
Statistical analyses
The results of the clinicopathological analyses were expressed as mean±standard deviation (SD). Mean comparisons were performed using Student’s t test and group comparisons were performed using a 2-tailed chi-square test, Fisher’s exact test, or linear-by-linear associations. Pearson’s correlation coefficient was used to examine the association of genes. The experimental data are expressed as mean±standard error of the mean and comparisons were performed using Mann-Whitney U test. Statistical analysis was performed using Prism (GraphPad Software, SanDiego, CA, USA) and SPSS version 25.0 for Windows (IBM Corp., Armonk, NY, USA). All P values are twosided. Image analyses was conducted by ImageJ (National Institutes of Health) [32].
RESULTS
Comparisons of LXRα and LXRβ expression in human cancers
To investigate the LXRα and LXRβ expression status in human cancers, we collected and compared the mRNA levels of LXRα and LXRβ in twelve types of cancers from TCGA. LXRα showed a relatively high expression in liver cancer (59.88± 23.64), but low expression in thyroid cancer (16.36±8.09), suggesting LXRα to be a major isoform of LXRs in liver cancer (Fig. 1A). In addition, pancreatic, lung, colorectal, stomach and breast cancers also showed relatively high expression of LXRα. Interestingly, LXRβ was highly expressed in thyroid cancer (121.6±27.02) whereas LXRβ expression was relatively low in liver and breast cancer (70.73±28.28 and 66.29±24.87, respectively), also suggesting that LXRβ might be a dominant isoform of LXRs in thyroid cancer (Fig. 1B). In addition to thyroid cancer, melanoma, endometrial, prostate, cervical and lung cancer samples also showed high LXRβ mRNA expression. Compatible with our findings from mRNA expression analyses in TCGA, the protein levels analyzed using HPA (Fig. 1C, D) revealed that LXRβ expression was more frequently observed in thyroid cancer than in other types of human cancers. The mRNA and protein expression levels appeared to be somewhat different in certain types of cancers (Fig. 1C, D). For example, LXRα mRNA expression in cervical cancer seemed to be relatively low but protein expression of LXRα was frequently observed (Fig. 1A, C). However, mRNA and protein levels of LXRβ in thyroid cancer were higher than those in other types of cancers (Fig. 1B, D).
To verify that LXRβ is a dominant isoform of LXRs in thyroid cancer, we compared the mRNA expression levels of LXRα and LXRβ between normal thyroid tissues and PTC samples. As shown in Fig. 1E, LXRα expression in PTC was lower than that in unpaired and paired normal thyroid tissues (unpaired t test, P<0.0001; paired t test, P<0.01). On the contrary, LXRβ expression in PTC was higher than that in unpaired and paired normal thyroid tissues (unpaired t test, P<0.0001; paired t test, P<0.0001) (Fig. 1F). Taken together, LXRβ might be a major isoform of LXRs in thyroid cancer and its expression might be increased in a tumor-specific manner.
Clinicopathological and molecular feature analyses of patients with PTC according to LXRβ expression
To understand the role of increased LXRβ expression in patients with PTC, we compared the clinicopathological characteristics of the patients according to LXRβ expression levels (Table 1). For this analysis, we divided THCA data into four groups according to LXRβ expression and used the upper and lower quartiles as groups with high and low LXRβ expression, respectively. Interestingly, male patients were more frequently observed in the group with high LXRβ expression. In addition, TNM stage was more advanced in the group with high LXRβ expression. These data suggested that LXRβ expression might be related to aggressive clinical features of PTC. In line with clinicopathological analyses, we also compared the molecular features of patients with PTC according to the LXRβ expression levels (Table 2). A high mRNA cluster number, which has been known to be associated with aggressive molecular features, was related to high LXRβ expression. Interestingly, although the frequency of RAS proto-oncogene, GTPase (RAS) and BRAF driver mutations was not different between the two groups, the BRAF/RAS score was lower in the group with high LXRβ expression, suggesting that this group has BRAF-like features, a molecular feature of aggressive thyroid cancer. In support of this finding, the differentiation score was also lower in the group with high LXRβ expression (P<0.0001). In summary, these clinicopathological and molecular profile analyses suggest that LXRβ might have a role in tumor progression of patients with PTC.
Transcriptional target analyses of LXRβ in PTC
To gain a deeper understanding of the molecular mechanisms by which LXRβ influences tumor behavior, we performed correlation analyses between the expression status of LXRβ and its 24 well-known target genes regulating lipid metabolism. Consistent with previous reports, this analysis indicated that LXRβ expression was correlated with the expression of LXRβ target genes [33], especially apolipoprotein C1 (APOC1), APOC2, apolipoprotein E (APOE), ATP binding cassette subfamily G member 8 (ABCG8), sterol regulatory element-binding protein 1c (SREBP1c; sterol regulatory element binding transcription factor 1 [SREBF1]), and thyroid hormone responsive (THRSP, SPOT14) (Fig. 2A). Among these positively correlated genes, APOE and SREBP1c showed positive correlation co-efficiency (Pearson R=0.31, P<0.0001; and Pearson R=0.30, P<0.0001, respectively). To support the results of our correlation analyses, we compared the expression of the target genes between normal thyroid tissues and cancer samples. Four out of six representative target genes, including APOC1, APOE, SREBP1c, and SPOT14 showed increased expression in cancer samples compared with that in their unpaired (Fig. 2B) and paired normal tissues (data not shown). In line with our correlation analyses, all six representative target genes showed high expression patterns in the group with high LXRβ expression; APOC1, APOE, and SREBP1c showed a strong statistical significance (Fig. 2C).
After confirming the positive relationship of LXRβ with its target genes, we investigated the other genes correlated with LXRβ by GSEA using TCGA THCA. Using gene sets provided by Kyoto Encyclopedia of Genes and Genomes (KEGG) signaling pathways, two gene sets (base excision repair and ribosome) were coordinately enriched in the group with high LXRβ expression. Interestingly, using gene sets listed by the gene ontology (GO) term biological process, all gene sets enriched in the group with high LXRβ expression were related to the ribosome, such as large ribosomal subunit, ribosomal subunit, ribosome and ribosome assembly (Fig. 3A). Taken together, GSEA using two representative lists of gene sets consistently indicated that high LXRβ expression was correlated with high expression of ribosome-related genes. Examination of individual genes related to ribosomes in KEGG and GO gene lists indicated that ribosomal protein L15 (RPL15), RPL11, RPL19, and ribosomal protein S30 (RPS30) were the significantly up-regulated genes (Fig. 3B).
To verify our findings, we determined the Pearson correlation coefficient between these genes, which showed a positive relationship of LXRβ with RPL15, RPL11, RPL19, and RPS30 (Fig. 4A). In addition, the selected target gene was found to be increased in tumor than normal through THCA, GSE83520 data set (Fig. 4B). To verify a positive relationship of LXRβ with RPL15, RPL11, RPL19, and RPS30 in vitro, we collected in vitro datasets from GEO and analyzed GSE77039 (LXR agonist), GSE41223 (LXR agonist), and GSE96545 (LXR knockout mice). All these in vitro data analyses consistently supported our findings (Supplemental Fig. S1). Using our own validation sets from PTC patients who underwent thyroidectomy, LXRβ expression was up-regulated in cancer samples compared to that in matched contra-lateral normal tissues. In addition, RPL15, RPL11, RPL19, and RPS30 expression was also increased in cancer samples; this was confirmed by our quantitative analyses using ImageJ (Fig. 4C, D).
DISCUSSION
The incidence of obesity and obesity-related diseases, such as dyslipidemia and diabetes mellitus, has been increasing worldwide [34]. Obesity has been suggested as a cause of various cancers such as colorectal cancer, endometrial cancer, and esophageal cancer [15]. Various mechanisms have been suggested to explain how these obesity-related cancers develop. The most studied biological mechanisms are those related to insulin and insulin-like growth factors, sex steroid and sex-steroid binding globulin, and adipokines such as adiponectin and leptin [35]. In addition, obesity-related inflammatory cytokines, nuclear factor κB system, altered immune response, and oxidative stress have been suggested as alternative mechanistic explanations [36]. Hypertension, a representative disease of metabolic syndrome, increases the risk of renal cancer by exerting high mechanical pressure on the renal parenchyma [37,38]. Acid reflux induced by increased abdominal pressure could also increase the risk of esophageal adenocarcinoma [39]. Excessive intake of iodide has been suggested as a cause of thyroid cancer [16]. However, these traditional approaches for the mechanistic explanation of obesity-related cancer have not investigated the significance of intracellular lipid accumulation and metabolic changes in cancer cells. A recent study on nonalcoholic steatohepatitis (NASH) revealed that changes in lipid metabolism caused by ectopic lipid accumulation resulted in hepatocellular carcinoma (HCC) [40]. The mechanism explaining the development of HCC from NASH involves lipid-laden hepatocytes that produce lipid paradoxically; moreover, genes related to lipogenesis, such as SREBP1c, play an important role in this process [41]. As a result, lipid-regulating mechanisms, such as de novo lipogenesis, are important in NASH and cancer development.
Thyroid cancer has also been suggested as an obesity-related cancer according to recent epidemiologic, experimental, and clinical studies [15,17,19,23-25]. The mechanism by which obesity can increase the risk of thyroid cancer has remained unelucidated, although excessive intake of iodide has been suggested as one of the causes. Similar to the other internal organs such as the liver and heart, the thyroid gland can also be subject to ectopic fat accumulation, which causes pathological changes in lipid metabolism of the thyroid gland [14,20]. In line with these previous studies, lipid metabolism in normal thyroid follicular cells and thyroid cancer cells might be important research targets to understand thyroid cancer metabolism. In our study, we first investigated the expression levels of LXRα and LXRβ in human cancers. LXRs are members of the nuclear receptor family of transcription factors and important regulators of lipid and glucose metabolism; they are closely related to peroxisome proliferator activated receptors, farnesoid X receptors, and retinoid X receptors, all of which are important players in intracellular lipid metabolism [42,43]. Therefore, we postulated that gaining insights into LXR regulation might be an initial step to understand the intracellular lipid metabolism in human cancers. In liver, lung, and colorectal cancers, LXRα mRNA expression was remarkably increased compared to that in other cancers; LXRα expression is restricted to the liver, intestine, lung, and so forth, and is the highest in the liver, hence the name LXRα. LXRβ has been known to be expressed in almost every tissue, and was hence named a ubiquitous receptor [3,44]. Interestingly, LXRβ expression was higher in thyroid cancer than in other types of cancers, suggesting that intracellular lipid metabolism in thyroid cancer might play an important role in tumor behavior. Supporting our presumption, the clinicopathological analyses according to LXRβ expression status revealed that high LXRβ expression was correlated with high frequency of male gender and advanced TNM stages although their P values were marginal. In addition, the molecular features of the group with high LXRβ expression revealed interesting findings. Although the prevalence of BRAF V600E mutation was not different, high LXRβ expression was correlated with a high mRNA cluster number, a lower BRAF/RAS score and a lower differentiation score, which implied a more aggressive and de-differentiated state. Thus, the characteristics of lipid metabolism might have a prognostic and biological impact on thyroid cancer, independent of the type of driver mutation.
Although we could not verify the mechanistic link between LXRβ and ribosomes, GSEA and our validation experiments clearly showed the coordinate correlation between LXRβ and ribosome-related gene sets. Traditionally, lipid metabolism in cancer has been known to be involved in the synthesis of cell membranes essential for cell proliferation [45]. Our data indicating the close relationship of LXRβ with ribosome, suggested that increased lipid metabolism might be coupled with increased rate of protein synthesis, which is also essential process of cancer cell proliferation. Taken together, increased LXRβ expression might indicate that cancer cells were in high proliferative activity.
In fact, we analyzed overall survival and recurrence rate according to LXRβ expression status but unfortunately could not find any significant difference between high and low LXRβ expression groups (data not shown). This data indicated that high LXRβ expression could not be used as a useful prognostic factor. However, as we described above, because clinical and biological features of high LXRβ expression group consistently showed aggressive features, the prognostic impact of LXRβ expression needs further study in the future. In this study, we used American Joint Committee on Cancer (AJCC) 7th edition to analyze clinicopathological features of TCGA THCA although individual pathologic reports could be obtained as source files. This might be the limitation of this paper to understand the prognostic impact of LXRβ expression according to the currently used AJCC 8th edition.
These data suggest that lipid metabolism in thyroid cancer might be closely related to translational activity which is a known signature of aggressive tumors. In future, the crosslink between lipid metabolism and master regulators of ribosomal activity, such as mammalian target of rapamycin (mTOR)/S6K pathways, might be an interesting research topic in thyroid cancer metabolism [46-48].
SUPPLEMENTARY MATERIAL
Acknowledgements
This work was funded by the German Research Foundation (STE 1096/1-3) and the German Federal Ministry of Education and Research to the German Centre of Diabetes Research.
Notes
CONFLICTS OF INTEREST
No potential conflict of interest relevant to this article was reported.
AUTHOR CONTRIBUTIONS
Conception or design: J.L., Y.S.J. Acquisition, analysis, or interpretation of data: S.J., H.K., M.J.C. Drafting the work or revising: S.J., I.K.K., J.L., Y.S.J. Final approval of the manuscript: S.J., I.K.K., H.K., M.J.C., J.L., Y.S.J.