Next Article in Journal
Risk of Cancer in Patients with Inflammatory Bowel Diseases and Keys for Patient Management
Next Article in Special Issue
Incidental Node Metastasis as an Independent Factor of Worse Disease-Free Survival in Patients with Papillary Thyroid Carcinoma
Previous Article in Journal
Pancreatic Cancer Incidence Trends by Race, Ethnicity, Age and Sex in the United States: A Population-Based Study, 2000–2018
Previous Article in Special Issue
Recurrence Risk Evaluation in Patients with Papillary Thyroid Carcinoma: Multicenter Machine Learning Evaluation of Lymph Node Variables
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Profile of MicroRNAs Associated with Death Due to Disease Progression in Metastatic Papillary Thyroid Carcinoma Patients

by
Ana Kober Leite
1,2,*,
Kelly Cristina Saito
3,
Thérèse Rachell Theodoro
4,
Fátima Solange Pasini
5,
Luana Perrone Camilo
2,
Carlos Augusto Rossetti
2,
Beatriz Godoi Cavalheiro
1,
Venâncio Avancini Ferreira Alves
6,
Luiz Paulo Kowalski
1,
Maria Aparecida Silva Pinhal
4,
Edna Teruko Kimura
3 and
Leandro Luongo Matos
1,2
1
Head and Neck Surgery Department, Instituto do Câncer do Estado de São Paulo, Hospital das Clínicas da Faculdade de Medicina da Universidade de São Paulo, São Paulo 05403-010, SP, Brazil
2
Faculdade Israelita de Ciências da Saúde Albert Einstein, Hospital Albert Einstein, São Paulo 05653-120, SP, Brazil
3
Department of Cell and Developmental Biology, Institute of Biomedical Sciences, University of São Paulo, São Paulo 05508-000, SP, Brazil
4
Discipline of Biochemistry, Department of Morfology and Physiology, Faculdade de Medicina Do ABC, São Paulo 09060-870, SP, Brazil
5
Centro de Investigação Translacional em Oncologia, Instituto do Câncer do Estado de São Paulo Paulo (ICESP), Hospital das Clínicas (HCFMUSP), Faculdade de Medicina da Universidade de São Paulo, Sao Paulo 01246-000, SP, Brazil
6
Pathology Department, Instituto do Câncer do Estado de São Paulo, Laboratório de Investigação Médica 14 (LIM14), Faculdade de Medicina da Universidade de São Paulo, São Paulo 01246-903, SP, Brazil
*
Author to whom correspondence should be addressed.
Cancers 2023, 15(3), 869; https://doi.org/10.3390/cancers15030869
Submission received: 29 December 2022 / Revised: 25 January 2023 / Accepted: 30 January 2023 / Published: 31 January 2023
(This article belongs to the Special Issue Advances in Understanding the Immune Network of Thyroid Cancers)

Abstract

:

Simple Summary

Papillary thyroid cancer is not an aggressive cancer, even when metastasis are present; therefore, treatment has been downgraded in recent years. However rare, mortality exists and finding the factors associated with mortality is essential and has not yet been fully accomplished; the answer probably lies in molecular factors. This study aimed to find MicroRNAs related to death in metastatic patients and found that patients who died due to progression of PTC had higher expression levels of miR-101-3p, miR-17-5p, and miR-191-5p compared with those of patients with stable metastatic disease. These findings are new in the literature and may open new doors in differentiating the patients who need more aggressive treatment from those who do not.

Abstract

Papillary thyroid carcinoma (PTC) is the most common neoplasm of the endocrine system and has an excellent long-term prognosis, with low rates of distant metastatic disease. Although infrequent, there are cases of deaths directly related to PTC, especially in patients with metastatic disease, and the factors that could be associated with this unfavorable outcome remain a major challenge in clinical practice. Recently, research into genetic factors associated with PTC has gained ground, especially mutations in the TERT promoter and BRAF gene. However, the role of microRNAs remains poorly studied, especially in those patients who have an unfavorable outcome at follow-up. This paper aims to evaluate molecular markers related to the different pathological processes of PTC, as well as the histological characteristics of the neoplasm, and to compare this profile with prognosis and death from the disease using an analysis of patients treated for metastatic disease in a single tertiary cancer center. Evaluation of microRNA expression in paraffin-embedded tumor specimens was carried out by quantitative PCR using the TaqMan® Low Density Array (TLDA) system. Metastatic patients who died from progression of PTC had higher expressions of miR-101-3p, miR-17-5p, and miR-191-5p when compared to patients with stable metastatic disease. These findings are of great importance but should be considered as preliminary because of the small sample.

1. Introduction

Papillary thyroid carcinoma (PTC) is the most common malignant neoplasm of the thyroid gland; however, it represents only 1% of all diagnosed cases of cancer and approximately 0.2% of deaths from this disease [1,2]. In most patients, it is curable and has an indolent course and excellent prognosis. The literature shows a 10-year disease-specific mortality of less than 5%. However, despite the low mortality, there are cases of deaths directly related to PTC, and factors that may be associated with unfavorable outcomes in these cases remain a major challenge in clinical practice [3].
In recent years, a series of genetic and molecular alterations related to the genesis of PTC have been studied and partially defined [4]. It is known that the initiation and progression of PTC involve the accumulation of several genetic and epigenetic modifications [5,6].
MicroRNAs (miRNAs) are small endogenous non-coding RNAs that mediate the regulation of gene expression at the post-transcriptional level. The expression of miRNAs plays an important role in various physiological processes, including cell cycle control, differentiation, proliferation, apoptosis, cell homeostasis, and organogenesis. The dysregulation of miRNA expression is considered as an important step in tumor development and progression in several human tissues. The overexpression of specific miRNAs can lead to the suppression of tumor suppressor genes, and the underexpression of other miRNAs can result in the increased expression of oncogenes; both situations affect cell proliferation, differentiation, and apoptosis, leading to tumor growth and progression. miRNAs have also been associated with several stages of the invasion–metastasis cascade in cancer, such as modification of the tumor microenvironment, local invasion, survival in the vasculature, and proliferation at distant sites [7]. However, the exact biological role of miRNAs in thyroid carcinogenesis remains uncertain. In an independent study, researchers analyzed the expression of miRNAs in different thyroid tumors and showed the dysregulation of miRNAs in tumor tissues compared with normal ones [8].
In recent decades, there has been a change in the causes of death related to PTC. Early series reported mortality as mainly related to local disease progression, whereas recent series showed distant metastasis (DM) as the main cause of death, which may be explained by the increase in aggressive resections and widespread application of radioactive iodine (RAI) therapy [1]. However, why histologically similar tumors have an indolent course in some cases but lead to tumor progression and death in other cases remains unknown. This is also observed among metastatic PTC patients; some patients have stable and asymptomatic disease for many years, whereas others have aggressive behavior, leading to death. Therefore, the present study evaluated the expression of miRNAs related to death in patients with metastatic PTC.

2. Materials and Methods

2.1. Data Survey

Among all patients with PTC treated by the Head and Neck Surgery Department at the Instituto do Câncer do Estado de São Paulo, Hospital das Clínicas of the Faculdade de Medicina, Universidade de São Paulo (Icesp, HCFMUSP), 108 patients were identified as having DM (at the time of initial diagnosis of PTC or during follow-up). These patients were fully managed at the Institution or referred from another location after the initial surgery. Initial diagnosis and surgical treatment were performed from 1986 to 2015 and, during this period, 3555 patients were surgically treated for PTC at the institution.
The pathological stage system (pTNM) was revised according to the 8th edition of the AJCC/UICC Cancer Staging Manual [9] and initial stage (at the time of diagnosis) was used.
The diagnosis of DM was based on clinical and/or positive findings on iodine-131 whole-body scanning (WBS), radiography (chest X-ray), computed tomography (CT), MRI, or tissue biopsy.
Dedifferentiation was defined based on the pathological examination of the tumor or metastatic disease and considered when there was an abrupt transformation from a well-differentiated tumor to a tumor with a high-grade morphology without the original distinctive histological features [10].
Iodine avidity was determined by visual uptake at the known site of metastatic disease on WBS after the first adjuvant treatment with RAI under hypothyroid conditions. The tumor was considered refractory to RAI therapy based on the 2015 American Thyroid Association Guidelines [11].
Patients were considered to have died from PTC if extensive, symptomatic, and progressive metastatic disease was present at the last follow-up and it was declared as the primary cause of death or as a significant contributing cause of death on the death certificate.

2.2. Selection of Samples for Molecular Study

The material used in the molecular analyses was a small part of the tumor remnant in formalin-fixed paraffin-embedded (FFPE) blocks obtained from anatomopathological examination.
For molecular analyses, all patients who died from the disease and had representative tumor material from FFPE samples at the institution were included, and among these patients, paired samples were selected at a 2:1 ratio from a second group of metastatic patients with a long follow-up period. Therefore, 16 patients who died due to PTC progression and 8 “controls” (i.e., metastatic patients with stable disease during follow-up) were included, giving a total of 24 patients (samples). The final sample consisted of 17 deceased patients and 7 controls (one patient died of DM during follow-up and at the end of the analyses).

2.3. Sequencing for BRAF and TERT Mutations

For TERT promoter and BRAF mutation analyses genomic DNA extraction was performed using a nucleic acid isolation protocol with FFPE samples as described by Kizys et al. [12]. The amplification of the TERT promoter and BRAF sequences was performed by nested PCR. In brief, an extended sequence of the gene was amplified in the first PCR (PCR I), and the product was used as a template in the second round of PCR (PCR II). The specific primer pairs for each sequence were designed by the Primer-BLAST program [13] and are presented in Table S1 (Supplementary Material). The PCR reaction mixture in a final volume of 25 µL was prepared with 400 ng of genomic DNA, 400 nM primers, 5 µL of MyTaq DNA Polymerase Buffer 5×, 0.2 mM dNTP, and 1 µL of MyTaq DNA Polymerase (Bioline®; Meridian Bioscience®, Cincinnati, OH, USA). The amplification was performed using the 2720 Thermal Cycler (Applied Biosystems, Thermo Fischer Scientific®, Waltham, MA, USA) under the following conditions: pre-denaturation at 95 °C for 5 min, followed by 40 cycles of denaturation at 94 °C for 30 s, annealing at 51–60 °C for 35–40 s, and extension at 72 °C for 35 s, and a final extension phase at 72 °C for 7 min. Therefore, 2 µL of the amplified product from PCR I was used as the template for PCR II. The final product was recovered from 1% agarose gel electrophoresis and purified using the QIAquick Gel Extraction Kit (Qiagen®, Hilden, Germany). For capillary sequencing by the Sanger method, the sample was mixed with a fluorescent nucleotide solution according to the manufacturer’s recommendations (BigDye® Terminator v3.1 Cycle Sequencing Kit; Applied Biosystems, Thermo Fischer Scientific®, Waltham, MA, USA) and amplified using the ABI 3730 DNA Analyzer (Applied Biosystems, Thermo Fischer Scientific®, Waltham, MA, USA). The generated sequences were analyzed by Sequence Scanner Software v2.0 (Applied Biosystems, Thermo Fischer Scientific®, Waltham, MA, USA), and the identification of mutational points was obtained using an alignment program in the BLAT database [14].

2.4. MicroRNA Detection Technique

Sixty-four microRNAs were selected for this study. The selection was carried out based on an extensive literature review, including miRNAs previously associated with PTC, tumor aggressiveness, metastatic disease, and invasion processes such as epithelial mesenchymal transition in thyroid cancer and other tumors.
The level of miRNA expression was assessed by reverse transcription-quantitative polymerase chain reaction (RT-qPCR) according to a previously stablished procedure [15]. RNA was extracted from 10 sections (5 μm) of paraffin-embedded tissue using the MagMAX FFPE RNA Ultra Kit (Applied Biosystems, Foster City, CA, USA) according to the manufacturer’s protocol. RT-qPCR was performed using the TaqMan Low Density Array system (Thermo Fisher Scientific®, Waltham, MA, USA). The RNA samples were subjected to reverse transcriptase assay with polyadenylation to obtain the cDNA product. All amplifications were performed in triplicate, and the cycle threshold (CT) values were determined with a cloud platform (Thermo Fisher Scientific®, Waltham, MA, USA). For methodological rigor, miRNAs with CT values > 36 were regarded as indeterminate and not included. In addition, only miRNAs detected in at least 70% of the samples were considered for statistical analysis. Therefore, the following 18 miRNAs were not included in the analyses: miR-129-5p, miR-130b-3p, miR-137-3p, miR-138-2-3p, miR-146b-3p, miR-155-3p, miR-17-3p, miR-187-3p, miR-302c-3p, miR-30e-5p, miR-34b-3p, miR-34c-5p, miR-455-3p, miR-4788, miR-506-3p, miR-654-3p, miR-9-5p, and miR-98-5p. miRNA expression was normalized between samples using the quantile method with the Expander software [16] selecting the most stable expressions among the samples (endogenous normalizers). The mean of the miRNAs let-7g-5p and miR-181a-5p was used to normalize the expression values of the remaining miRNAs. The expression levels of these miRNAs were quantified using the 2-ΔCt method to identify upregulated or downregulated miRNAs.
Commercially available primers sold by Thermo Fisher Scientific® were used and are listed in Table S2 (Supplementary Material). They are constructed based on data available at miRbase.org.

2.5. Statistical Analyses

For the statistical analyses, the values of continuous variables are presented as the mean and standard deviation (SD) or standard error (SE). For categorical variables, absolute and relative frequencies are presented. The cutoff values of quantitative variables were determined by receiver operating characteristic (ROC) analysis, and the cutoff values for qualitative variables were determined based on the association of categories. For the comparison of two sample populations, the Mann–Whitney test was used, and for three or more populations, the Kruskal–Wallis test with Dunn’s auxiliary test was employed. Comparisons of the frequency of a phenomenon between groups for categorical variables were performed with the chi-square test or Fisher’s exact test. Variables with p < 0.10 in univariate analysis were used for multivariate analysis, which was performed using linear regression models for independent quantitative variables. The inclusion of variables with p < 0.10 in the multivariate analysis excludes possible dependent and confounding variables. In addition to sensitivity and specificity, accuracy estimates (area under the ROC curve—AUC) were determined by ROC curve analysis (positive Z-scores were used in cases of upregulated miRNAs, and negative Z-scores were used in the cases of downregulated miRNAs). In all analyses, SPSS version 28.0 (IBM®, Endicott, NY, USA) was used, and the Multiple Experiment Viewer (MEV version 4.5.0; available at: http://mev.tm4.org, accessed on 22 February 2022) was used to obtain the heatmaps.

3. Results

The descriptive data of the patients included in the study are shown in Table 1. The group of patients who died from disease progression had more metastatic disease at multiple sites and demonstrated disease dedifferentiation, and there was a higher number of RAI-refractory patients. For other variables, the distribution was relatively homogeneous between the two groups, including age, stage, histopathological characteristics, and treatment received. Notably, we selected patients with a long follow-up period among metastatic patients with stable disease as we aimed to exclude follow-up time as a possible bias. Therefore, the follow-up time was approximately three times longer in the group of patients alive at the end of the follow-up compared with that of patients who died due to disease progression.
Four patients (23.5%) who died from the disease had a mutation in the BRAF gene, and eight patients (47%) had a mutation in the TERT promoter gene. Three individuals had both mutations. In the group of metastatic patients with stable disease, the rates of mutations in the BRAF and TERT genes were 28.6% (two cases) and 100%, respectively. In addition, BRAF and TERT mutations were present in 25% and 62.5% of the sample of metastatic individuals, respectively, regardless of the outcome (death or alive with stable disease).
The let-7c-5p, let-7e-5p, 101-3p, 16-5p, 17-5p, 181b-5p, 191-5p, 19a-3p, 19b-3p, 200c-3p, 20a- 5p, 29a-3p, 30b-5p, 30c-5p, 30d-5p, 31-5p, and 34a-5p miRNAs were differentially expressed between the two groups in univariate analysis, as shown in Table 2.
miRNAs with p-values lower than 0.1 were selected for multivariate linear regression analysis to identify those expressed independently. The data for the miRNAs that were significant in the multivariate analyses are shown in Table 3 and the variables that were included but not significant are shown in Table S3 (Supplementary Material). In comparison with patients with stable metastatic disease, patients who died due to disease progression had higher expression levels of miR-101-3p, miR-17-5p, and miR-191-5p. Three models based on miR-17-5p alone; miR-17-5p in association with mir-101-3p; and miR-101-3p, miR-17-5p, and miR-191-5p together were proposed, as shown in Table 3.
The results of heatmap analysis using the three proposed models with the miRNAs differentially expressed in multivariate analysis (overexpressed miR-17-5p, miR-101-3p, and miR-191-5p) are presented in Figure 1. A clear segregation between patients who died from PTC progression compared with metastatic patients without disease progression was achieved, especially when considering the association of the three selected miRNAs (model 3).
Values greater than 0.861, 2.155, and 0.455 were established as the ROC curve cutoff values of miR-101-3p, miR-17-5p, and miR-191-5p, respectively, for the evaluation of death of patients with metastatic PTC. The three models showed high accuracy in predicting the mortality of patients with metastatic PTC, as shown in Table 4. None of the patients with these microRNAs expression bellow the cutoff value died; therefore, the specificity is 100%. However, the low number of patients in both groups did not allow the analysis of the risk of death by logistic regression of the cutoff values of the three miRNAs.

4. Discussion

The present study demonstrated the overexpression of the miRNAs 17-5p, 101-3p, and 191-5p in patients who died due to metastatic PTC progression. Furthermore, the expression of the three miRNAs was lower in metastatic patients with an indolent course (defined as alive and without metastatic disease progression in a long follow-up period) than in patients with worse PTC progression. This is the first study to identify a different profile of miRNAs associated with death in PTC patients, the knowledge of which could facilitate the understanding of the biology of these tumors and contribute to accurate follow-up and specific therapeutics in the future.
Various miRNAs in PTC samples were previously evaluated in different studies in the last decade. However, most studies compare PTC with normal tissue while the present study included only metastatic PTC patients. Liu et al. [17] evaluated a large profile of miRNAs in 126 samples from PTC patients by comparing their expression levels with those in matching normal thyroid tissues in order to understand the role of miRNAs in PTC carcinogenesis. The authors identified the overexpression and underexpression of 180 and 68 miRNAs in cancer samples, respectively, including the overexpression of miR-222 (target genes KIT and AXIN2), miR-15a (target genes AXIN2 and FOXO1), and miR-221 (target gene KIT) and the underexpression of miR-206 (target gene MET), miR-299-3p (target gene ITGAV), miR-103 (target gene ITGA2), and miR-101 (target gene ITGA3), which was also downregulated in the group of patients with indolent metastatic disease in our study. In another study, Rosignolo et al. [18] evaluated the expression of 754 different miRNAs in PTC patients. The authors identified eight overexpressed miRNAs (miR-221-3p, miR-222-3p, miR-146a-5p, miR-24-3p, miR-146b-5p, miR-103a-3p, miR-28-3p, and miR-191-5p) in serum of PTC patients compared with patients with benign nodules. Similarly, miR-191-5p was overexpressed in the group of patients who died due to PTC progression in our study. Serum levels of miR-146a-5p and miR-221-3p were also associated with worse response to therapy during follow-up. Bioinformatic analysis for microarray and RNA-seq has also been used for identification of potential target miRNAs. Using these methods, Wang et al. [19] identified miR-146b-5p, miR-15a-5p, miR-21-5p, miR221-3p, and miR-222-3p as differently expressed in PTC compared to normal tissue.
There are limited studies on the three differentially expressed miRNAs identified in our study. Moreover, the evaluation of death as an outcome for metastatic PTC patients is unique in the literature, thus limiting the comparison of our findings with those of previous studies.
miR-17-5p is a member of the miR-17-92 cluster, located on chromosome 13q31-32. It was previously described as related to tumor proliferation, invasion, migration, and apoptosis inhibition in thyroid cancer, targeting PTEN and other pathways [20,21,22]. Specifically, miR-17 has been identified as upregulated in anaplastic thyroid carcinoma (ATC) compared with adjacent normal thyroid tissues [23,24]. In another study, the inhibition of the expression of this miRNA was associated with growth arrest and the induction of anaplastic cell apoptosis via caspase activation [25]. However, another group found that the inhibition of miR-17 increased MYCN and c-MYC expression, which increased pri-miR-17-92 transcripts and the expression of oncogene miRNAs (miR-18a and miR-19a), demonstrating that miR-17 inhibition did not regulate thyroid tumor growth in two different cell lines of anaplastic thyroid carcinoma [26]. The conflicting results suggest that different miRNAs have different expression levels and functions depending on the tumor microenvironment and behavior. This notion is consistent with the findings of Liu et al., who observed that miR-17-5p exhibited a dual function in tumor progression [27]. Additionally, the increased serum level of miR-17-5p was found in aggressive cases of medullary thyroid cancer compared with PTC [28].
miR-101 is an important tumor suppressor and its underexpression is associated with thyroid malignancy. It is also considered as an “anti-metastatic” miRNA in PTC, which can inhibit migration by targeting the Rac1 gene [29,30,31] and induce apoptosis by enhancing the TRAIL pathway [32]. In previous studies, miR-101-3p overexpression inhibited the invasion and migration of PTC cells in vitro through different pathways [33,34]. Zhao et al. [35] demonstrated that miR-101 reduced PTC cell proliferation, apoptosis resistance, and invasion by targeting USP22. In addition, its underexpression was associated with lymph node metastasis and poor prognosis (lower survival) among PTC patients, which is consistent with our results showing the underexpression of this miRNA in metastatic PTC patients. Nevertheless, we also observed the higher expression of miR-101 in metastatic cases that progressed to death, demonstrating that the precise mechanism of miR-101-3p in aggressive PTC remains unclear.
The miR-101 gene has two copies, hsa-miR-101.1 located in the intron of the LOC107983962 gene (hg38: chr1:65,058,434–6,505,850 [−]) and hsa-miR101.2 located in the intron of the RCL1 gene (hg38: chr9:4,850,297–4,850,375 [+]). Furthermore, the miRNA is not a single sequence but a series of multiple isomiRs with expression and sequence diversity. Guo et al. [36] estimated five isomiRs with two seeds from the miR-101-3p.2 locus, showing a significant difference in gaining or losing target mRNAs. Deregulated isomiRs may also contribute to the abnormal expression patterns of their target mRNAs, and lncRNAs may disrupt the isomiR:mRNA network, further controlling mRNA expression.
In some tumors, miR-101-3p was confirmed to target the Moloney murine leukemia virus 1 (Pim-1) oncogene [37]. The oncogenic Pim kinase proteins (Pim-1/2/3) regulate tumorigenesis through the phosphorylation of essential proteins that control the cell cycle and proliferation [38,39,40]. Therefore, miR-101-3p could be considered as a potential antineoplastic molecule. However, unlike most kinases, Pim does not have regulatory domains; once translated, it is constitutively active [41]. Taken together, the findings suggest that despite miR-101-3p promoting a decrease in the synthesis of Pim kinases, the level of these proteins may not be affected as their degradation could be inhibited. Therefore, even with the overexpression of miR-101-3p, the level of Pim kinases might be high, which would explain the worse prognosis of patients with PTC [42,43]. In addition, researchers have recently demonstrated the role of Pim kinases in immunomodulation; their increase may cause tumor immune escape and help tumor cells survive [44].
MDM4 directly binds to p53, resulting in degradation. Therefore, it is considered as a negative regulator of the p53 suppression pathway and has been identified as a target site for miR-191 [45]. The binding of miR-191 to MDM4-C allele mRNA, which could result in the increased expression of MDM4, may explain tumor progression in different types of cancers such as ovarian cancer and retinoblastoma [46,47]. The overexpression of miR-191-5p was reported only in the serum of PTC patients compared with normal thyroid tissues [18], indicating its role in thyroid carcinogenesis. Furthermore, Colamaio et al. [48] reported the high expression of miR-191 in PTC and ATC and its downregulation in follicular thyroid carcinoma and the follicular variant of PTC, suggesting the role of this miRNA in thyroid carcinogenesis. This is the first report of the association of miR-191-5p with the aggressive behavior of PTC.
The strength of the present study is that it is the first to evaluate the role of different miRNAs in death due to tumor progression in metastatic PTC patients, a rare but fatal event. The rarity of this event makes it difficult to study in PTC, and large series are not available. Nevertheless, the small sample size is a limitation of our study, and our results should be considered as preliminary. However, since commercially available primers were used, the results can be reproduced.
Some discrepancies between the overexpression of the three miRNAs identified in the present study and their underexpression in cases of poor treatment response or outcome reported in the literature, mostly in in vitro studies, may be attributed to a possible type I error; however, this is a statistically remote possibility. A more plausible interpretation of the findings would be that the role of the differentially expressed miRNAs has not been well established for this particular outcome. Another limitation is the fact that this study only analyzed a selected group of microRNAs; therefore, there could be other miRNAs significant to this outcome that where not evaluated.

5. Conclusions

The present study demonstrated that metastatic patients who died due to the progression of PTC had higher expression levels of miR-101-3p, miR-17-5p, and miR-191-5p compared with those of patients who had stable metastatic disease. Further studies on the miRNA profile of patients with unfavorable outcomes are required to better understand the role of miRNAs not only in thyroid carcinogenesis but also in the development of metastatic disease, tumor progression, and death.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cancers15030869/s1, Table S1: Specific primer pairs for each sequence used in TERT promoter and BRAF mutation detection; Table S2: Primers sequences for the studied microRNAs; Table S3: Variables not significant in the multivariate linear regression analyses.

Author Contributions

Conceptualization, A.K.L. and L.L.M.; methodology, A.K.L. and L.L.M.; formal analysis, T.R.T., F.S.P. and K.C.S.; investigation, L.P.C. and C.A.R.; data curation, A.K.L. and B.G.C.; writing—original draft preparation, A.K.L. and L.L.M.; writing—review and editing, L.P.K.; E.T.K., V.A.F.A. and M.A.S.P.; supervision, L.L.M.; project administration, A.K.L.; funding acquisition, L.L.M. and E.T.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Conselho Nacional de Desenvolvimento Científico e Tecnológico—CNPq, Public Brazilian Agency, (grants number: 407520/2018-0, 311210/2021-0, 409443/2021-2) and Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP), Public Brazilian Agency (grants numer: 2020/12459-1 and 2019/25116-8).

Institutional Review Board Statement

The project was approved by the Institutional Review Board (IRB; CAAE number 44997215.1.0000.0065).

Informed Consent Statement

All patients who were being followed up at the institution’s outpatient clinic completed the Informed Consent Form (ICF) at the time of their usual outpatient appointment. The IRB waived the requirement of the ICF for patients who died.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Leite, A.K.N.; Cavalheiro, B.G.; Kulcsar, M.A.; Hoff, A.O.; Brandao, L.G.; Cernea, C.R.; Matos, L.L. Deaths related to differentiated thyroid cancer: A rare but real event. Arch. Endocrinol. Metab. 2017, 61, 222–227. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Leite, A.K.N.; de Araujo Filho, V.J.F.; de Matos, L.L. Are All Small Papillary Thyroid Cancers Harmless? A Word of Caution. Endocr. Pract. 2019, 25, 1233–1234. [Google Scholar] [CrossRef]
  3. Nunes, K.S.; Matos, L.L.; Cavalheiro, B.G.; Magnabosco, F.F.; Tavares, M.R.; Kulcsar, M.A.; Hoff, A.O.; Kowalski, L.P.; Leite, A.K. Risk factors associated with disease-specific mortality in papillary thyroid cancer patients with distant metastases. Endocrine 2022, 75, 814–822. [Google Scholar] [CrossRef]
  4. Pacini, F.; Castagna, M.G. Approach to and treatment of differentiated thyroid carcinoma. Med. Clin. N. Am. 2012, 96, 369–383. [Google Scholar] [CrossRef] [PubMed]
  5. Bose, S.; Walts, A.E. Thyroid fine needle aspirate: A post-Bethesda update. Adv. Anat. Pathol. 2012, 19, 160–169. [Google Scholar] [CrossRef] [PubMed]
  6. Maciel, R.M.; Kimura, E.T.; Cerutti, J.M. Pathogenesis of differentiated thyroid cancer (papillary and follicular). Arq. Bras. Endocrinol. Metabol. 2005, 49, 691–700. [Google Scholar] [CrossRef] [Green Version]
  7. Jansson, M.D.; Lund, A.H. MicroRNA and cancer. Mol. Oncol. 2012, 6, 590–610. [Google Scholar] [CrossRef] [Green Version]
  8. Yuan, Z.M.; Yang, Z.L.; Zheng, Q. Deregulation of microRNA expression in thyroid tumors. J. Zhejiang Univ. Sci. B 2014, 15, 212–224. [Google Scholar] [CrossRef] [Green Version]
  9. Amin, M.B.; Edge, S.B.; Greene, F.L.; Byrd, D.R.; Brookland, R.K.; Washington, M.K.; Gershenwald, J.E.; Compton, C.C.; Hess, K.R.; Sullivan, D.C.; et al. AJCC Cancer Staging Manual, 8th ed.; Springer: Chicago, IL, USA, 2017. [Google Scholar]
  10. Leite, A.K.; Kulcsar, M.A.; de Godoi Cavalheiro, B.; de Mello, E.S.; Alves, V.A.; Cernea, C.R.; Matos, L.L. Death Related to Pulmonary Metastasis in Patients with Differentiated Thyroid Cancer. Endocr. Pract. 2017, 23, 72–78. [Google Scholar] [CrossRef]
  11. Haugen, B.R.; Alexander, E.K.; Bible, K.C.; Doherty, G.M.; Mandel, S.J.; Nikiforov, Y.E.; Pacini, F.; Randolph, G.W.; Sawka, A.M.; Schlumberger, M.; et al. 2015 American Thyroid Association Management Guidelines for Adult Patients with Thyroid Nodules and Differentiated Thyroid Cancer: The American Thyroid Association Guidelines Task Force on Thyroid Nodules and Differentiated Thyroid Cancer. Thyroid 2016, 26, 1–133. [Google Scholar] [CrossRef]
  12. Kizys, M.M.; Cardoso, M.G.; Lindsey, S.C.; Harada, M.Y.; Soares, F.A.; Melo, M.C.; Montoya, M.Z.; Kasamatsu, T.S.; Kunii, I.S.; Giannocco, G.; et al. Optimizing nucleic acid extraction from thyroid fine-needle aspiration cells in stained slides, formalin-fixed/paraffin-embedded tissues, and long-term stored blood samples. Arq. Bras. Endocrinol. Metab. 2012, 56, 618–626. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Ye, J.; Coulouris, G.; Zaretskaya, I.; Cutcutache, I.; Rozen, S.; Madden, T.L. Primer-BLAST: A tool to design target-specific primers for polymerase chain reaction. BMC Bioinform. 2012, 13, 134. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Kent, W.J.; Sugnet, C.W.; Furey, T.S.; Roskin, K.M.; Pringle, T.H.; Zahler, A.M.; Haussler, D. The human genome browser at UCSC. Genome Res. 2002, 12, 996–1006. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Matos, L.L.; Menderico Junior, G.M.; Theodoro, T.R.; Pasini, F.S.; Ishikawa, M.M.; Ribeiro, A.A.B.; de Mello, E.S.; Pinhal, M.; Moyses, R.A.; Kulcsar, M.A.V.; et al. Cancer-associated fibroblast regulation by microRNAs promotes invasion of oral squamous cell carcinoma. Oral. Oncol. 2020, 110, 104909. [Google Scholar] [CrossRef] [PubMed]
  16. Ulitsky, I.; Maron-Katz, A.; Shavit, S.; Sagir, D.; Linhart, C.; Elkon, R.; Tanay, A.; Sharan, R.; Shiloh, Y.; Shamir, R. Expander: From expression microarrays to networks and functions. Nat. Protoc. 2010, 5, 303–322. [Google Scholar] [CrossRef] [PubMed]
  17. Liu, X.; He, M.; Hou, Y.; Liang, B.; Zhao, L.; Ma, S.; Yu, Y.; Liu, X. Expression profiles of microRNAs and their target genes in papillary thyroid carcinoma. Oncol. Rep. 2013, 29, 1415–1420. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Rosignolo, F.; Sponziello, M.; Giacomelli, L.; Russo, D.; Pecce, V.; Biffoni, M.; Bellantone, R.; Lombardi, C.P.; Lamartina, L.; Grani, G.; et al. Identification of Thyroid-Associated Serum microRNA Profiles and Their Potential Use in Thyroid Cancer Follow-Up. J. Endocr. Soc. 2017, 1, 3–13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Wang, J.; Wu, L.; Jin, Y.; Li, S.; Liu, X. Identification of key miRNAs in papillary thyroid carcinoma based on data mining and bioinformatics methods. Biomed. Rep. 2020, 12, 11–16. [Google Scholar] [CrossRef] [Green Version]
  20. Liu, X.; Liu, G.; Lu, Y.; Shi, Y. Long non-coding RNA HOTAIR promotes cell viability, migration and invasion in thyroid cancer cells by sponging miR-17-5p. Neoplasma 2020, 67, 229–237. [Google Scholar] [CrossRef]
  21. Geng, X.; Sun, Y.; Fu, J.; Cao, L.; Li, Y. MicroRNA-17-5p inhibits thyroid cancer progression by suppressing Early growth response 2 (EGR2). Bioengineered 2021, 12, 2713–2722. [Google Scholar] [CrossRef]
  22. Shi, Y.P.; Liu, G.L.; Li, S.; Liu, X.L. miR-17-5p knockdown inhibits proliferation, autophagy and promotes apoptosis in thyroid cancer via targeting PTEN. Neoplasma 2020, 67, 249–258. [Google Scholar] [CrossRef] [PubMed]
  23. Das, P.K.; Asha, S.Y.; Abe, I.; Islam, F.; Lam, A.K. Roles of Non-Coding RNAs on Anaplastic Thyroid Carcinomas. Cancers 2020, 12, 3159. [Google Scholar] [CrossRef] [PubMed]
  24. Fuziwara, C.S.; Saito, K.C.; Kimura, E.T. Thyroid Follicular Cell Loss of Differentiation Induced by MicroRNA miR-17-92 Cluster Is Attenuated by CRISPR/Cas9n Gene Silencing in Anaplastic Thyroid Cancer. Thyroid 2020, 30, 81–94. [Google Scholar] [CrossRef] [Green Version]
  25. Takakura, S.; Mitsutake, N.; Nakashima, M.; Namba, H.; Saenko, V.A.; Rogounovitch, T.I.; Nakazawa, Y.; Hayashi, T.; Ohtsuru, A.; Yamashita, S. Oncogenic role of miR-17-92 cluster in anaplastic thyroid cancer cells. Cancer Sci. 2008, 99, 1147–1154. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Sweat, Y.; Ries, R.J.; Sweat, M.; Su, D.; Shao, F.; Eliason, S.; Amendt, B.A. miR-17 acts as a tumor suppressor by negatively regulating the miR-17-92 cluster. Mol. Ther. Nucleic Acids 2021, 26, 1148–1158. [Google Scholar] [CrossRef]
  27. Liu, L.; Yang, J.; Zhu, X.; Li, D.; Lv, Z.; Zhang, X. Long noncoding RNA H19 competitively binds miR-17-5p to regulate YES1 expression in thyroid cancer. FEBS J. 2016, 283, 2326–2339. [Google Scholar] [CrossRef] [Green Version]
  28. Zhang, A.; Wang, C.; Lu, H.; Chen, X.; Ba, Y.; Zhang, C.; Zhang, C.Y. Altered Serum MicroRNA Profile May Serve as an Auxiliary Tool for Discriminating Aggressive Thyroid Carcinoma from Nonaggressive Thyroid Cancer and Benign Thyroid Nodules. Dis. Mrk. 2019, 2019, 3717683. [Google Scholar] [CrossRef] [Green Version]
  29. Wang, C.; Lu, S.; Jiang, J.; Jia, X.; Dong, X.; Bu, P. Hsa-microRNA-101 suppresses migration and invasion by targeting Rac1 in thyroid cancer cells. Oncol. Lett. 2014, 8, 1815–1821. [Google Scholar] [CrossRef] [Green Version]
  30. Lin, X.; Guan, H.; Li, H.; Liu, L.; Liu, J.; Wei, G.; Huang, Z.; Liao, Z.; Li, Y. miR-101 inhibits cell proliferation by targeting Rac1 in papillary thyroid carcinoma. Biomed. Rep. 2014, 2, 122–126. [Google Scholar] [CrossRef] [Green Version]
  31. Lima, C.R.; Gomes, C.C.; Santos, M.F. Role of microRNAs in endocrine cancer metastasis. Mol. Cell Endocrinol. 2017, 456, 62–75. [Google Scholar] [CrossRef]
  32. Zhu, J.; Li, Z. Overexpression of miR-101 promotes TRAIL-induced mitochondrial apoptosis in papillary thyroid carcinoma by targeting c-met and MCL-1. Oncotarget 2017, 8, 108665–108675. [Google Scholar] [CrossRef] [Green Version]
  33. Du, Y.L.; Liang, Y.; Cao, Y.; Liu, L.; Li, J.; Shi, G.Q. LncRNA XIST Promotes Migration and Invasion of Papillary Thyroid Cancer Cell by Modulating MiR-101-3p/CLDN1 Axis. Biochem. Genet. 2021, 59, 437–452. [Google Scholar] [CrossRef] [PubMed]
  34. Chen, F.; Yang, D.; Ru, Y.; Cao, S.; Gao, A. MicroRNA-101 Targets CXCL12-Mediated Akt and Snail Signaling Pathways to Inhibit Cellular Proliferation and Invasion in Papillary Thyroid Carcinoma. Oncol. Res. 2019, 27, 691–701. [Google Scholar] [CrossRef] [PubMed]
  35. Zhao, H.; Tang, H.; Huang, Q.; Qiu, B.; Liu, X.; Fan, D.; Gong, L.; Guo, H.; Chen, C.; Lei, S.; et al. MiR-101 targets USP22 to inhibit the tumorigenesis of papillary thyroid carcinoma. Am. J. Cancer Res. 2016, 6, 2575–2586. [Google Scholar] [PubMed]
  36. Guo, L.; Dou, Y.; Yang, Y.; Zhang, S.; Kang, Y.; Shen, L.; Tang, L.; Zhang, Y.; Li, C.; Wang, J.; et al. Protein profiling reveals potential isomiR-associated cross-talks among RNAs in cholangiocarcinoma. Comput. Struct. Biotechnol. J. 2021, 19, 5722–5734. [Google Scholar] [CrossRef]
  37. Liu, X.Y.; Liu, Z.J.; He, H.; Zhang, C.; Wang, Y.L. MicroRNA-101-3p suppresses cell proliferation, invasion and enhances chemotherapeutic sensitivity in salivary gland adenoid cystic carcinoma by targeting Pim-1. Am. J. Cancer Res. 2015, 5, 3015–3029. [Google Scholar]
  38. Mondello, P.; Cuzzocrea, S.; Mian, M. Pim kinases in hematological malignancies: Where are we now and where are we going? J. Hematol. Oncol. 2014, 7, 95. [Google Scholar] [CrossRef] [Green Version]
  39. Eerola, S.K.; Kohvakka, A.; Tammela, T.L.J.; Koskinen, P.J.; Latonen, L.; Visakorpi, T. Expression and ERG regulation of PIM kinases in prostate cancer. Cancer Med. 2021, 10, 3427–3436. [Google Scholar] [CrossRef]
  40. Wen, Q.L.; Yi, H.Q.; Yang, K.; Yin, C.T.; Yin, W.J.; Xiang, F.Y.; Bao, M.; Shuai, J.; Song, Y.W.; Ge, M.H.; et al. Role of oncogene PIM-1 in the development and progression of papillary thyroid carcinoma: Involvement of oxidative stress. Mol. Cell Endocrinol. 2021, 523, 111144. [Google Scholar] [CrossRef]
  41. Warfel, N.A.; Kraft, A.S. PIM kinase (and Akt) biology and signaling in tumors. Pharmacol. Ther. 2015, 151, 41–49. [Google Scholar] [CrossRef] [Green Version]
  42. Losman, J.A.; Chen, X.P.; Vuong, B.Q.; Fay, S.; Rothman, P.B. Protein phosphatase 2A regulates the stability of Pim protein kinases. J. Biol. Chem. 2003, 278, 4800–4805. [Google Scholar] [CrossRef] [PubMed]
  43. Park, Y.K.; Obiang-Obounou, B.W.; Lee, K.B.; Choi, J.S.; Jang, B.C. AZD1208, a pan-Pim kinase inhibitor, inhibits adipogenesis and induces lipolysis in 3T3-L1 adipocytes. J. Cell Mol. Med. 2018, 22, 2488–2497. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Szydłowski, M.; Prochorec-Sobieszek, M.; Szumera-Ciećkiewicz, A.; Derezińska, E.; Hoser, G.; Wasilewska, D.; Szymańska-Giemza, O.; Jabłońska, E.; Białopiotrowicz, E.; Sewastianik, T.; et al. Expression of PIM kinases in Reed-Sternberg cells fosters immune privilege and tumor cell survival in Hodgkin lymphoma. Blood 2017, 130, 1418–1429. [Google Scholar] [CrossRef] [Green Version]
  45. Mohammad Khanlou, Z.; Pouladi, N.; Hosseinpour Feizi, M.; Pedram, N. Lack of Associations of the MDM4 rs4245739 Polymorphism with Risk of Thyroid Cancer among Iranian-Azeri Patients: A Case-Control Study. Asian Pac. J. Cancer Prev. 2017, 18, 1133–1138. [Google Scholar] [CrossRef] [PubMed]
  46. McEvoy, J.; Ulyanov, A.; Brennan, R.; Wu, G.; Pounds, S.; Zhang, J.; Dyer, M.A. Analysis of MDM2 and MDM4 single nucleotide polymorphisms, mRNA splicing and protein expression in retinoblastoma. PLoS ONE 2012, 7, e42739. [Google Scholar] [CrossRef]
  47. Wynendaele, J.; Bohnke, A.; Leucci, E.; Nielsen, S.J.; Lambertz, I.; Hammer, S.; Sbrzesny, N.; Kubitza, D.; Wolf, A.; Gradhand, E.; et al. An illegitimate microRNA target site within the 3’ UTR of MDM4 affects ovarian cancer progression and chemosensitivity. Cancer Res. 2010, 70, 9641–9649. [Google Scholar] [CrossRef] [Green Version]
  48. Colamaio, M.; Borbone, E.; Russo, L.; Bianco, M.; Federico, A.; Califano, D.; Chiappetta, G.; Pallante, P.; Troncone, G.; Battista, S.; et al. miR-191 down-regulation plays a role in thyroid follicular tumors through CDK6 targeting. J. Clin. Endocrinol. Metab. 2011, 96, E1915–E1924. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Heatmaps demonstrating the differently expressed microRNAs and hierarchical cluster analyses of all samples at the differentiation between dead patients due to papillary thyroid carcinoma progression from those alive with stable metastatic disease (without tumor progression at long follow-up), considering only the microRNAs differently expressed at the multivariate analysis (overexpression of miR-17-5p, miR-101-3p, and miR-191-5p). Yellow spots represent samples with high expression and those in blue low represent microRNAs expression.
Figure 1. Heatmaps demonstrating the differently expressed microRNAs and hierarchical cluster analyses of all samples at the differentiation between dead patients due to papillary thyroid carcinoma progression from those alive with stable metastatic disease (without tumor progression at long follow-up), considering only the microRNAs differently expressed at the multivariate analysis (overexpression of miR-17-5p, miR-101-3p, and miR-191-5p). Yellow spots represent samples with high expression and those in blue low represent microRNAs expression.
Cancers 15 00869 g001
Table 1. Descriptive data of both groups of patients with metastatic papillary thyroid carcinoma included in the study.
Table 1. Descriptive data of both groups of patients with metastatic papillary thyroid carcinoma included in the study.
VariableDeath Due to Cancer Progression n = 17 (%)Alive Metastatic Patients n = 7 (%)
Sex
Female12 (70.6)4 (57.1)
Male5 (29.4)3 (42.9)
Age
Mean ± SD (years-old)54.8 ± 13.256 ± 13.4
Histopathological data
Extrathyroidal extension9 (52.9)3 (50.0)
Vascular invasion11 (73.3)4 (66.7)
Multifocality9 (52.9)4 (57.1)
Size (mean ± SD; cm)4.7 ± 3.64.1 ± 2.6
Initial pT classification
pT14 (23.5)1 (16.7)
pT22 (11.8)1 (16.7)
pT35 (29.4)2 (33.3)
pT4a2 (11.8)2 (33.3)
pT4b4 (23.5)0
Initial pN classification
pN08 (50)2 (28.6)
pN1a3 (18.8)1 (14.3)
pN1b5 (31.3)4 (57.1)
Treatment
Total thyroidectomy16 (94.1)7 (100)
Central compartment neck dissection8 (50)5 (71.4)
Level II-V neck dissection5 (31.3)4 (57.1)
R0 surgery13 (86.7)6 (85.7)
Radioiodine therapy (RIT)12 (70.6)7 (100)
Number of RIT
One6 (50)4 (57.1)
>26 (50)3 (42.9)
Distant metastases17 (100)7 (100)
Lung16 (94.1)7 (100)
Liver4 (23.5)0
Bones12 (70.6)1 (14.3)
Multiple sites14 (82.4)1 (14.3)
Diagnosis with the primary tumor12 (70.6)4 (57.1)
Dedifferentiation8 (47.1)0
Follow-up
Radioiodine refractory disease9 (75)3 (42.9)
Regional failure7 (50)1 (14.3)
Follow-up time (mean ± SD; months)51.9 ± 37.9158.6 ± 32.9
Table 2. MicroRNAs differently expressed between death patients due to progression of metastatic papillary thyroid carcinoma and those alive with stable metastatic disease (without tumor progression at long follow-up).
Table 2. MicroRNAs differently expressed between death patients due to progression of metastatic papillary thyroid carcinoma and those alive with stable metastatic disease (without tumor progression at long follow-up).
miRDeath Due to Tumor ProgressionAlive with Metastatic Diseasep-Value
(Mann–Whitney)
MEANSEMEANSE
let-7b-5p0.0820.0360.0760.0380.413
let-7c-5p5.1853.0025.8411.1050.020
let-7d-5p0.4450.0520.4700.0640.619
let-7e-5p2.9811.2435.4901.6700.020
let-7f-5p3.8550.7484.9501.1640.494
let-7i-5p20.4633.41411.6091.2680.089
miR-1-3p17.1297.26314.5966.5620.891
miR-101-3p3.8630.5170.7760.1540.005
miR-10b-5p0.4710.1240.7740.3030.590
miR-125a-5p7.2301.5507.4221.6330.664
miR-138-5p3.9190.9061.8110.6180.081
miR-141-3p17.3033.3369.4122.4340.172
miR-146b-5p2.6860.9134.9571.4190.179
miR-16-5p5.6101.4020.9260.1440.019
miR-17-5p1.8110.4170.3510.1170.005
miR-181b-5p0.2620.0480.4740.0530.014
miR-18a-5p0.0580.0250.0960.0940.364
miR-191-5p1.0570.1810.2360.044<0.001
miR-199a-3p1.3280.5120.7520.2300.757
miR-19a-3p0.5660.2120.0530.0190.003
miR-19b-3p2.2201.1340.1250.035<0.001
miR-200a-3p3.8420.6422.1830.3860.209
miR-200b-3p38.35111.60051.3679.2580.065
miR-200c-3p16.7834.20433.1084.6730.011
miR-203a-3p0.1010.0440.0560.0320.773
miR-205-5p2.7911.1541.2030.6910.602
miR-20a-5p0.5260.1110.1240.033<0.001
miR-21-5p55.03019.91619.7515.1680.383
miR-214-3p0.9650.5030.0940.0620.100
miR-221-3p23.8345.93918.5185.1181.000
miR-222-3p1.6670.4061.7210.3120.452
miR-29a-3p6.5871.1123.2600.3990.006
miR-30a-5p5.9691.7622.2990.4330.072
miR-30b-5p15.4183.4904.9580.9640.018
miR-30c-5p26.5615.9023.6620.6940.005
miR-30d-5p1.9300.3830.6630.1710.021
miR-30e-3p0.6910.1540.5500.2010.559
miR-31-5p2.2521.9221.4030.3710.017
miR-34a-5p3.7511.0375.6680.9500.032
miR-423-5p0.4540.0930.8960.2670.091
miR-4290.3660.0970.1470.0530.178
miR-483-3p1.2980.7520.2890.1030.831
miR-92a-3p15.5747.6104.6120.8870.671
Legend: SE = standard-error. Highlighted in bold are the variables with p < 0.05.
Table 3. Multivariate analyses (linear regression) identifying microRNAs of risk of death due to tumor progression in patients with metastatic papillary thyroid carcinoma.
Table 3. Multivariate analyses (linear regression) identifying microRNAs of risk of death due to tumor progression in patients with metastatic papillary thyroid carcinoma.
ModelUnstandardized CoefficientsStandardized Coefficientsp-Value95% CI (for Beta)
BIFBetaLowerUpper
1Constant−0.0210.113 0.860−0.3100.268
miR-17-5p0.5550.0730.9590.0010.3670.743
2Constant−0.1220.079 0.199−0.3430.098
miR-17-5p0.4510.0590.7790.0020.2880.614
miR-101-3p0.0660.0230.2940.0440.0030.129
3Constant−0.1830.036 0.015−0.298−0.068
miR-17-5p0.3410.0350.5890.0020.2300.452
miR-101-3p0.0710.0100.3180.0050.0410.102
miR-191-5p0.2650.0600.2310.0210.0750.456
Legend: SE = standard-error. Highlighted in bold are the variables with p < 0.05.
Table 4. Evaluation of microRNAs 17-5p, 101-3p, and 191-5p for death due to tumor progression prediction in patients with metastatic papillary thyroid carcinoma.
Table 4. Evaluation of microRNAs 17-5p, 101-3p, and 191-5p for death due to tumor progression prediction in patients with metastatic papillary thyroid carcinoma.
Micro-RNASensitivitySpecificityAUC95% CI (AUC)
LowerUpper
Model 1
miR-17-5p ≥ 0.861
0.9231.0000.9620.8721.051
Model 2
miR-17-5p ≥ 0.861
AND miR-101-3p ≥ 2.155
0.7861.0000.8930.7501.036
Model 3
miR-17-5p ≥ 0.861
AND miR-101-3p ≥ 2.155
AND miR-191-5p ≥ 0.455
0.7141.0000.8570.6941.021
Legend: AUC = area under the ROC curve.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Leite, A.K.; Saito, K.C.; Theodoro, T.R.; Pasini, F.S.; Camilo, L.P.; Rossetti, C.A.; Cavalheiro, B.G.; Alves, V.A.F.; Kowalski, L.P.; Pinhal, M.A.S.; et al. Profile of MicroRNAs Associated with Death Due to Disease Progression in Metastatic Papillary Thyroid Carcinoma Patients. Cancers 2023, 15, 869. https://doi.org/10.3390/cancers15030869

AMA Style

Leite AK, Saito KC, Theodoro TR, Pasini FS, Camilo LP, Rossetti CA, Cavalheiro BG, Alves VAF, Kowalski LP, Pinhal MAS, et al. Profile of MicroRNAs Associated with Death Due to Disease Progression in Metastatic Papillary Thyroid Carcinoma Patients. Cancers. 2023; 15(3):869. https://doi.org/10.3390/cancers15030869

Chicago/Turabian Style

Leite, Ana Kober, Kelly Cristina Saito, Thérèse Rachell Theodoro, Fátima Solange Pasini, Luana Perrone Camilo, Carlos Augusto Rossetti, Beatriz Godoi Cavalheiro, Venâncio Avancini Ferreira Alves, Luiz Paulo Kowalski, Maria Aparecida Silva Pinhal, and et al. 2023. "Profile of MicroRNAs Associated with Death Due to Disease Progression in Metastatic Papillary Thyroid Carcinoma Patients" Cancers 15, no. 3: 869. https://doi.org/10.3390/cancers15030869

APA Style

Leite, A. K., Saito, K. C., Theodoro, T. R., Pasini, F. S., Camilo, L. P., Rossetti, C. A., Cavalheiro, B. G., Alves, V. A. F., Kowalski, L. P., Pinhal, M. A. S., Kimura, E. T., & Matos, L. L. (2023). Profile of MicroRNAs Associated with Death Due to Disease Progression in Metastatic Papillary Thyroid Carcinoma Patients. Cancers, 15(3), 869. https://doi.org/10.3390/cancers15030869

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop