Peptidyl-prolyl cis/trans isomerase (Pin1), encoded by PIN1 gene with locus in chromosome 19p13, is an enzyme that catalytically induces conformational changes in proteins after phosphorylation on serine or threonine residues preceding proline (pSer/Thr-Pro motifs); in this way, it has an influence on protein interactions and intracellular localizations of proteins. The aim of the study were: 1) an assessment of PIN1 gene expression level in benign and malignant thyroid lesions; 2) the evaluation of possible correlations between gene expression and histopathological variants of papillary thyroid carcinoma (PTC) or tumour size, classified according to TNM classification of primary tumours (in case of PTC only); 3) the estimation of possible relationships between expression of the gene in question and patients' sex or age.
Seventy (70) tissue samples were analyzed: 32 cases of PTC, 7 cases of medullary thyroid carcinoma (MTC), 7 cases of follicular adenoma (FA), and 24 cases of nodular goitre (NG). In real-time polymerase chain reaction (real-time PCR), two-step RT-PCR (reverse transcriptase-polymerase chain reaction) in an ABI PRISM 7500 Sequence Detection System was employed. The PIN1 gene expression level was assessed, calculating the mean relative quantification rate (RQ rate) increase for each sample.
The level of PIN1 gene expression (compared to that in macroscopically unchanged thyroid tissue) was higher in PTC group than those in FA, MTC and/or NG groups, but the statistical significance was noted for difference between PTC and NG groups only. On the other hand, the differences of RQ rate value between different PTC variants were statistically insignificant. No correlations were found between RQ values and tumour size, as well as between RQ values and patients' sex or age in PTC group.
The PIN1 gene expression may have - in future - an important meaning in the diagnostics of PTC and in understanding its pathogenesis. However, our results - mostly due to the small number of cases - do not yet allow considering PIN1 gene as a prognostic molecular PTC marker.
The PIN1 gene encodes an essential peptidyl-prolyl cis/trans isomerase (PPIase; EC 220.127.116.11) which has a profound impact on key proteins involved in the regulation of cell cycle . The human PIN1 gene is mapped to chromosome 19p13 . Catalyzing conformational changes in substrates after phosphorylation, Pin1 recognizes phospho-serine or phospho-threonine bonds followed by proline (pSer/Thr-Pro motifs) and plays a vital role in protein folding or refolding . Through isomerization Pin1 can modulate the conformation of protein substrates, changing their enzymatic activity which leads to alteration in different cell function .
Pin1 enzyme activity is important for oncogenic and cell-signalling pathways via conformational changes following phosphorylation of many molecules (e.g. Bcl-2, p53, c-Jun, beta-catenin, cyclin D1 and RAF kinase) [4-7]. The Pin1-regulated alteration in protein interactions and stability have been shown to be associated with cell transformation and progression of certain types of human tumours (e.g., prostate, colon, lung, ovary and breast cancers, and lymphoma) [8-10]. It has also been reported that inhibition of Pin1 function in human tumour cells, using Pin1 antisense RNA, induces mitotic arrest and apoptosis . Further studies indicate that Pin1 plays a critical role in cell restriction points - particularly in G1/S phase; Pin1 modulates some important regulators of that cell cycle point, e.g. cyclin D1 and its upstream transcriptional factors. The increased Pin1 expression correlates significantly with overexpression of cyclin D1 mRNA and protein. Interestingly, Pin1 binds phosphorylated c-Jun, influencing the increase of its transcriptional activity towards cyclin D1 - by activating c-Jun/AP-1 and beta-catenin/T-cell factor (TCF) and, in this way it cooperates with Ras signalling pathway .
Pin1 is involved not only in the pathogenesis of human cancer, but - presumably - also in the pathogenesis of asthma and Alzheimer's disease .
Recently, in some reports - including the study on thyroid tumours - strong correlation between Pin1 and cyclin D1 immunoexpression and/or cyclin D1 mRNA and PIN1 expression via interaction with Wnt signalling pathway has been observed . Moreover, it has been suggested that Pin1 may promote cyclin D1 overexpression directly or through accumulation of beta-catenin in thyroid cancer cells .
There are rather scarce reports which are focused on evaluation of PIN1 mRNA expression in thyroid tumours. Therefore, in order to recognize the Pin1 role in the pathogenesis of thyroid gland, we have compared the levels of PIN1 gene expression in benign and malignant thyroid lesions. We have assessed PIN1 gene expression in real-time polymerase chain reaction (real-time PCR), two-step RT-PCR (reverse transcriptase-polymerase chain reaction) and we have evaluated the possible correlations between the gene expression and the histopathological variant of examined PTCs (classic, tall-cell, follicular), tumour size grouped according to TNM classification of primary tumours (pT1, pT2, pT3, pT4). We have also decided to estimate the possible relationships between expression of the gene in question and patients' sex or age.
Materials and methods
The Ethical Committee of the Medical University of Lodz approved the studies reported in the present article.
Thyroid tissue samples
Tumour tissue samples (100-150 mg) were obtained from 70 patients who had undergone surgery treatment at the Centre of Oncology, the Maria Skłodowska-Curie Institute in Gliwice, at the Holy Family Municipal Hospital in Lodz and at the Department of General and Colorectal Surgery, Medical University of Lodz, Poland during the years 2002-2007. Total tumour tissue samples, immediately after resection, were collected into lysis buffer (Buffer RLT, Qiagen Sciences, USA).
Among the studied samples 32 cases of primary papillary thyroid carcinoma - PTC (26 females, 6 males; mean age 41.8 ± 15.7; mean ± SD) were used. Histopathological variants of PTC were as follows: PTC classic variant (20 cases: 17 females, 3 males; mean age 36.5 ± 14.4), PTC follicular variant (9 cases: 6 females, 3 males; mean age 49.9 ± 16.0), PTC tall-cell variant (3 cases: 3 females; mean age 53.3 ± 5.1) (Table 1). Seven (7) cases of medullary thyroid carcinoma - MTC (6 females, 1 male; mean age 56.0 ± 17.6) and 7 cases of follicular adenoma - FA (6 females, 1 male; mean age 49.4 ± 14.8) were also included in the study. The control group consisted of 24 cases of nodular goitre - NG (22 females, 2 males; mean age 48.9 ± 14.7) (Table 2).
Table 1. Sex, age, histopathological diagnosis of malignant thyroid tumours in the studied patients
Table 2. Sex, age, histopathological diagnosis of benign thyroid lesions in the studied patients.
The whole study cohort comprised 70 persons: 60 females and 10 males, the mean age was 46.4 ± 15.8 years (age ranged from 16 to 75 years). Histopathological diagnoses for malignant thyroid lesions, according to WHO Classification of Tumours , have been obtained from pathomorphological reports and are presented in Table 1, together with TNM classification and AJCC stage groupings .
In the study design, macroscopically unchanged thyroid tissue served for the reference standard (calibrator).
Isolation of total RNA and reverse transcription (RT)
Total RNA was extracted from fresh tisssues, using RNeasy Protect Midi Kit, (Qiagen, Hilden, Germany), according to the manufacturer's recommendations. RNA concentration was spectrophotometrically assessed by measuring absorbance at 260 and 280 nm (Ultraspec 2000 UV/Visible Spectrophotometer, Pharmacia Biotech, Sweden).
Afterwards the reverse transcription (RT) was performed in a Personal Thermocycler (Eppendorf, Germany), using 1 μg of total RNA in the presence of 50 μM of oligo d(T)16 and Reverse Transcriptase MultiScribe™ (50 U/μl) in a total volume of 30 μl, including also: 10× TaqMan RT Buffer (containing 100 mM Tris-HCL pH 8.3, 500 mM KCL), MgCl2 solution (25 mM), dNTPs mixture (10 mM), RNAse Inhibitor (20 U/μl), and nuclease-free water (TaqMan Reverse Transcriptase Reagents, Applied Biosystems, Foster City, CA, USA). The reactions were conducted for 10 minutes at 25°C, followed by 30 minutes at 48°C, afterwards the samples were heated for 5 min to 95°C, and next placed at 4°C. The cDNA concentration and purity were spectrophotometrically assessed by measuring absorbance at 260 and 280 nm (Ultraspec 2000 UV/Visible Spectrophotometer, Pharmacia Biotech, Sweden).
Quantitative analysis of the relative amount of PIN1 mRNA by real-time PCR
An established Relative Quantification Polymerase Chain Reaction (RQ-PCR) assay for PIN1 mRNA expression was conducted in ABI PRISM 7500 Sequence Detection System (Applied Biosystems), according to the manufacturer's protocol. The PCR reactions for PIN1 gene were run with 50 ng of cDNA in a total volume of 50 μl, using TaqMan® Universal PCR Master Mix (Applied Biosystems, Foster City, CA, USA) and the predesigned and labelled primer/probe set (Assays-on-Demand™ Gene Expression assay mix, Hs01598309_ml, Applied Biosystems). After initial incubation at 50°C for 2 min to allow uracil-N-glycosylase (UNG) digestion, and at 95°C for 10 min to activate the AmpliTaq Gold®DNA polymerase, both of them were provided by the Universal PCR Master Mix, the samples were amplified through 40 biphasic cycles of 95°C for 15 sec and 60°C for 1 min (Table 3).
Table 3. Real-time PCR conditions for PIN1 gene amplification.
Macroscopically unchanged thyroid tissue was used as a calibrator (standard sample). Amplification reactions were done in triplicate for each sample (cDNA from the same PCR reaction but in separate wells). Controls with no template cDNA were used with each assay (negative control).
The expression levels of β-actin gene (ACTB) were measured, as endogenous control (reference gene), using the appropriate Assays-on-Demand™ Gene Expression product (Hs99999903_ml, Applied Biosystems Foster City, CA, USA).
Both gene expressions were measured for each thyroid lesion sample in the same PCR reaction but in separate wells.
Assays-on-Demand™ Gene Expression product consists of 20 × mix of unlabelled PCR primers (18 μM for each) and TaqMan® MGB probes (5 μM) with FAM™ (6-carboxy-fluorescein) at the 5' end as the reporter dye, and a non-fluorescent quencher (TAMRA, 6-carboxy-tetramethylrhodamine) at the 3' end. The two-minute, 50°C step was required for optimal AmpErase® UNG activity when TaqMan® Universal PCR Master Mix (P/N 4304437) was used.
The fluorescence signal was measured in real-time PCR in the extension phase of the PCR reaction, and the measurement, proportional to the quantity of sample cDNA in the reaction, was plotted as an amplification curve against cycle number. A threshold value of fluorescence in the exponential part of the amplification curve was selected, and - for each sample - the number of cycles was measured, which was needed by the signal to reach the threshold (threshold cycle, CT). The larger were the quantities of the starting material, the lower were CT values.
Data analysis of relative gene expression was performed with the Taq-Man SDS analysis software on an ABI PRISM 7500 Sequence Detection System (Applied Biosystems), and the results were exported to Excel sheets for further processing.
Fluorescence emission data were determined as CT values for each reaction and, for each sample, triplicate CT values were averaged. The average CT value for β-actin was subtracted from the average PIN1 CT value to yield the ΔCT value (ΔCT = CT target - CT reference). Normalization to the reference gene (β-actin) has been necessary to account for the variabilities in sample quantity and quality, and for the variabilities in PCR efficiency among samples. The assay, described in this report, involves the determination of a ΔΔCT value. This is calculated by: ΔΔCT = ΔCT test sample - ΔCT calibrator sample. The higher is the ΔΔCT value, the lower is the expression of PIN1 in the specimen. Fold-differences, representing relative expression results, are calculated using the following equation: relative fold increase = 2-ΔΔCT. In order to evaluate the analytical sensitivity, specificity, and accuracy of the assay, serial dilutions were performed.
The data were statistically analyzed, using Newman-Keuls' multiple comparison test followed by Kruskal-Wallis' test, in order to compare the level of expression values (RQ) among the four studied independent groups (PTC, MTC, FA, NG). Basic measures of location (i.e. mean, median), measures of dispersion (SD, SEM), and minimum and maximum values were calculated to provide detailed descriptions of gene expressions in selected groups. Additionally, box-and-whiskers plots were introduced to display the difference in gene expressions levels in different samples.
The Kruskal-Wallis' test was used to compare the level of expression values (RQ) among three variants of PTC (PTC v. classic, PTC v. tall-cell, PTC v. follicular). The Spearman's rank correlation coefficient and U Mann-Whitney's test were performed in order to correlate the level of expression of PIN1 gene (RQ values) with examined parameters (age, sex, tumour size, classified according to TNM definition of primary tumours - pT1, pT2, pT3, pT4).
The accepted level of statistical significance was at p < 0.05. The results are presented as means ± SEM and means ± SD values.
Statistica for Windows 7.0 program was applied for calculations.
The specimens were amplified in the ABI PRISM 7500 Sequence Detection System in reaction, containing primers and probes for PIN1 gene and a control gene, β-actin. The Sequence Detection System software, provided with the instrument, analyses the fluorescence data, generated during the reaction, and calculates the cycle number at which fluorescence crosses the threshold value (CT). Triplicate amplifications of the sample produced nearly identical, overlapping curves, from which CT values were calculated.
In our experiment low CT values for PIN1 gene in all malignant thyroid tissues (PTC - 32 and MTC - 7 specimens) were observed.
The assay, described in this paper, was based on the determination of a ΔΔCT value (ΔΔCT method) for each sample (the greater was the ΔΔCT value, the lower was the expression of PIN1 gene) and on calculating the difference in PIN1 expression level between tumour sample and calibrator (RQ value). The results for calculated RQ in the studied groups (PTC, MTC, FA, NG) are presented in Table 4. The RQ value for PTC, classic variant, ranged from 0.441 to 52.279 (the mean RQ was 8.43); for PTC, follicular variant, ranged from 0.804 to 9.124 (the mean RQ was 3.03) and for PTC, tall-cell variant, ranged from 0.363 to 12.616 (the mean RQ was 6.42).
Table 4. The expression level (RQ) of PIN1 gene, calculated by the ΔΔCT method in the studied groups.
Kruskal-Wallis' test, applied for RQ comparison between the four studied groups (PTC, MTC, FA, NG), showed statistically significant differences (p = 0.0233) among them (Figure 1). The results of Newman-Keuls' multiple comparison test revealed statistically significant differences in RQ values only between the PTC group (higher RQ values) and the control group - NG (p = 0.031977). In cases of other groups, the differences in RQ values were not significant (p > 0.05) (Table 5). Similarly, the results of Kruskal-Wallis' test did not show any significant differences among the particular variants of PTC (PTC v. classic, PTC v. tall-cell, PTC v. follicular) (p = 0.9501) (Figure 2). Additionally, in the whole studied group of PTCs, no correlations were found between PIN1 expression level and patients' sex, age or tumour size (p > 0.05). The results of Mann-Whitney's U test confirmed no significant differences in RQ values between the PTC group and MTC group (p = 0.6873) (Figure 3).
Table 5. Statistical analysis of results - comparison of RQ value among the four studied groups (PTC, MTC, FA, NG). Newman-Keuls' test, p - level of significance.
Figure 1. Box-and-whisker plots, representing the expression of PIN1 gene in the studied groups (PTC, MTC, FA, NG). Results are calculated as RQ values. Whiskers represent means ± SD (standard deviation) for particular groups. Boxes represent means ± SEM (standard error of mean). The data were statistically analyzed, using Newman-Keuls' multiple comparison test, followed by Kruskal-Wallis' test, in order to compare the level of expression values (RQ) among the four studied independent groups, p < 0.05; * - p < 0.031977 (Neuman-Keuls' test).
Figure 2. Box-and-whisker plots, representing the difference in the expression of PIN1 gene in the variants of PTC (PTC v. classic, PTC v. tall-cell, PTC v. follicular). Results are calculated as RQ values; no significant differences among studied groups are shown. Whiskers represent means ± SD (standard deviation) for particular groups. Boxes represent means ± SEM (standard error of mean). The results were statistically analyzed using Kruskal-Wallis' test, p = 0.9501.
Figure 3. Box-and-whisker plots, representing the difference in the expression of PIN1 gene between the PTC group and the MTC group. Results are calculated as RQ values; no significant differences between studied groups are demonstrated. Whiskers represent means ± SD (standard deviation) for particular groups. Boxes represent means ± SEM (standard error of mean). The results were statistically analyzed using Mann-Whitney's U test, p = 0.6873.
The molecular pathogenesis of thyroid cancer is very complex and still unclear. Recent molecular studies on thyroid carcinogenesis have proved that many specific genetic alterations in protooncogenes are associated with development of thyroid cancer [16,17]. The mutational activation of MAPK/ERK signalling pathway is crucial to understand the carcinogenesis in the thyroid. The mutations of genes coding effectors of MAPK/ERK pathway (RET, RAS, NTRK1, BRAF) which rarely overlap, may be identified in more than 70% of PTC cases [18-20]. Additionally, the chromosomal rearrangements (RET/PTC or Trk) appear to be of crucial importance in activation of MAPK/ERK pathway [21,22].
Despite the fact that constitutive activation of the above mentioned pathway is a very important for the development of thyroid cancer, many other molecular regulatory mechanisms, involved in thyroid cancer transformation, are recognized.
In our study, we have examined whether the PIN1 gene participates in thyroid carcinogenesis.
The phosphorylation of proteins on serine or threonine residues preceding proline (pSer/Thr-Pro motifs), catalized by the peptidyl-prolyl cis/trans isomerase Pin1, has been confirmed as a major cell transformation mechanism activating different oncogenic pathways in many types of human tumours (including thyroid cancer) [9,13,23,24]. In immunohistochemical studies of thyroid tumours, as well as of other types of cancers (e.g., oral squamous carcinoma), high Pin1 expression promotes cyclin D1 overexpression, directly or through accumulation of beta-catenin [13,25]. It has been reported that PIN1 expression is regulated by the transcription factor E2F and Neu/Ras signalling. The enhanced expression and stability of cyclin D1 by transcriptional activation and/or posttranslational stabilization by PIN1 gene or Pin1 protein may lead to cell proliferation and transformation. Therefore, PIN1 might be a crucial factor in modulating upregulation of cyclin D1 level by Neu/Ras oncogenic signalling [7,26].
Bearing in mind all these facts, in our present study, we have assessed PIN1 gene expression level in benign and malignant thyroid lesions. In our earlier study, we demonstrated that cyclin D1 gene expression levels were higher in malignant thyroid tumours (PTC, MTC), when compared to those in NG or FA groups . Our present results have revealed statistically significant differences in expression of PIN1 mRNA between PTC group and benign thyroid lesions (FA, NG) (higher RQ value in PTC). It should be stressed that our data are in accordance with observations of other authors . Nakashima et al.  have observed overexpression of cyclin D1 mRNA in PTC (thyroid tumours from a radiocontaminated area), taking course with high expression of PIN1 gene. In the quoted paper, correlation between high Pin1 protein level and cyclin D1 immunoexpression level in PTC has also been confirmed . It is important that Pin1 can promote cyclin D1 overexpression directly or through beta-catenin stabilization [6,13]. Our present results, which have focused on PIN1 mRNA expression, as well as our earlier data, concerning cyclin D1 mRNA overexpression in thyroid malignancy, may suggest that one of the major regulatory mechanisms in cell transformation in thyroid carcinogenesis is PIN1 gene activity . It is tempting to speculate that PIN1 gene positively regulates cyclin D1 gene function not only through posttranslational stabilization but also by an effect on transcriptional level.
In our present study, we have applied real-time relative quantification PCR assay for PIN1 mRNA expression and, employing this extremely sensitive and precise method, we have confirmed overexpression of PIN1 gene in PTC. These results may suggest a diagnostic utility of PIN1 mRNA expression - as a malignant lesion marker - for thyroid PTC. However, the results should be confirmed in larger number of PTC cases. The utility of PIN1 gene may be - in future - very crucial for diagnostics and treatment of thyroid cancer (especially of PTC), based on the facts that overexpression of Pin1 has been found to be an excellent prognostic marker in some other human cancers (e.g., prostate, breast, colorectal cancers) [8,10,23]. In these immunohistochemical studies, strong relationship between Pin1 protein level and clinical outcome has been observed [8,10,23]. Unfortunately, it is difficult to discuss the relationship between PIN1 mRNA and prognosis in PTC, because - according to our knowledge - there were no reported studies, focusing on the evaluation of PIN1 gene expression in association with pathological features and clinical outcome of PTC, so far. Our present data have not documented any significant differences in the expression level of PIN1 mRNA among the particular variants of PTC (classic vs. tall-cell, vs. follicular variants). Moreover, no correlations were found between PIN1 expression level and patients' sex, age or tumour size in PTC group.
In summary, PIN1 gene expression may have - in future - certain significance in PTC diagnostics and in understanding its pathogenesis. However, our results - due to the small number of cases - do not yet allow considering PIN1 gene as a prognostic molecular PTC marker.
The authors declare that they have no competing interests.
AL designed and coordinated the study and wrote the manuscript; EB and KC carried out the molecular genetic studies; JL participated in performing molecular studies; WK coordinated the collection of tissue samples; AC-M participated in performing molecular studies and drafted the manuscript.
This study was supported by the funds from the Medical University of Lodz, Poland (project no. 502-11-570).
Prog Cell Cycle Res 2003, 5:477-487. PubMed Abstract
Wulf GM, Ryo A, Wulf GG, Lee SW, Niu T, Petkova V, Lu KP: Pin1 is overexpressed in breast cancer and cooperates with Ras signaling in increasing the transcriptional activity of c-Jun towards cyclin D1.
Rippmann JF, Hobbie S, Daiber C, Guillard B, Bauer M, Birk J, Nar H, Garin-Chesa P, Rettig WJ, Schnapp A: Phosphorylation-dependent proline isomerization catalyzed by Pin1 is essential for tumor cell survival and entry into mitosis.
Nakashima M, Meirmanov S, Naruke Y, Kondo H, Saenko V, Rogounovitch T, Shimizu-Yoshida Y, Takamura N, Namba H, Ito M, Abrosimov A, Lushnikov E, Roumiantsev P, Tsyb A, Yamashita S, Sekine I: Cyclin D1 overexpression in thyroid tumours from a radio-contaminated area and its correlation with Pin1 and aberrant beta-catenin expression.
Endocr Regul 2000, 34:99-113. PubMed Abstract
Neuroendocrinol Lett 2007, 28:77-105. PubMed Abstract
Kimura ET, Nikiforova MN, Zhu Z, Knauf JA, Nikiforov YE, Fagin JA: High prevalence of BRAF mutations in thyroid cancer: genetic evidence for constitutive activation of the RET/PTC-RAS-BRAF signaling pathway in papillary thyroid carcinoma.
Soares P, Trovisco V, Rocha AS, Lima J, Castro P, Preto A, Maximo V, Botelho T, Seruca R, Sobrinho-Simoes M: BRAF mutations and RET/PTC rearrangements are alternative events in the pathogenesis of PTC.
Frattini M, Ferrario C, Bressan P, Balestra D, De Cecco L, Mondellini P, Bongarzone I, Collini P, Gariboldi M, Pilotti S, Pierotti MA, Greco A: Alternative mutations of BRAF, RET, and NTRK1 are associated with similar but distinct gene expression patterns in papillary thyroid cancer.
Brzeziańska E, Karbownik M, Migdalska-Sęk M, Pastuszak-Lewandoska D, Włoch J, Lewiński A: Molecular analysis of the RET and NTRK1 gene rearrangements in papillary thyroid carcinoma in the Polish population.
Neuroendocrinol Lett 2007, 28:341-350. PubMed Abstract