Comprehensive MicroRNA Expression Profiling Identifies Novel Markers in Follicular Variant of Papillary Thyroid Carcinoma

Abstract

Background:

Follicular variant of papillary thyroid carcinoma (FVPTC) shares features of papillary (PTC) and follicular (FTC) thyroid carcinomas on a clinical, morphological, and genetic level. MicroRNA (miRNA) deregulation was extensively studied in PTCs and FTCs. However, very limited information is available for FVPTC. The aim of this study was to assess miRNA expression in FVPTC with the most comprehensive miRNA array panel and to correlate it with the clinicopathological data.

Methods:

Forty-four papillary thyroid carcinomas (17 FVPTC, 27 classic PTC) and eight normal thyroid tissue samples were analyzed for expression of 748 miRNAs using Human Microarray Assays on the ABI 7900 platform (Life Technologies, Carlsbad, CA). In addition, an independent set of 61 tumor and normal samples was studied for expression of novel miRNA markers detected in this study.

Results:

Overall, the miRNA expression profile demonstrated similar trends between FVPTC and classic PTC. Fourteen miRNAs were deregulated in FVPTC with a fold change of more than five (up/down), including miRNAs known to be upregulated in PTC (miR-146b-3p, -146-5p, -221, -222 and miR-222-5p) and novel miRNAs (miR-375, -551b, 181-2-3p, 99b-3p). However, the levels of miRNA expression were different between these tumor types and some miRNAs were uniquely dysregulated in FVPTC allowing separation of these tumors on the unsupervised hierarchical clustering analysis. Upregulation of novel miR-375 was confirmed in a large independent set of follicular cell derived neoplasms and benign nodules and demonstrated specific upregulation for PTC. Two miRNAs (miR-181a-2-3p, miR-99b-3p) were associated with an adverse outcome in FVPTC patients by a Kaplan–Meier (p<0.05) and multivariate Cox regression analysis (p<0.05).

Conclusions:

Despite high similarity in miRNA expression between FVPTC and classic PTC, several miRNAs were uniquely expressed in each tumor type, supporting their histopathologic differences. Highly upregulated miRNA identified in this study (miR-375) can serve as a novel marker of papillary thyroid carcinoma, and miR-181a-2-3p and miR-99b-3p can predict relapse-free survival in patients with FVPTC thus potentially providing important diagnostic and predictive value.

Introduction

P apillary thyroid carcinoma (PTC) is the most common type of thyroid malignancy (1). It comprises several variants that are classified on the basis of their morphologic features. Follicular variant of papillary carcinoma (FVPTC) is the most common variant, and is diagnosed in about 15–30% of PTC (1,2). Historically, all thyroid tumors with follicular growth were diagnosed as follicular thyroid carcinomas (FTC). The first histologic description of FVPTC occurred in 1960, but nevertheless, it was not until 1977 when this entity was further characterized in a landmark paper by Chem and Rosai and found its way into daily clinical practice (2,3). Since then, the scientific debate about this remarkable entity is ongoing. While these tumors have characteristic nuclear features of PTC, they demonstrate a follicular growth pattern and complete lack of well-formed papillae (2). Moreover, the metastatic pattern of these tumors often resembles classic PTC with cervical lymph-node metastasis, but some cases may demonstrate distant hematologic metastases, a finding more common in FTC (4).

At the molecular level, FVPTC also appears as an entity located between PTC and FTC. It frequently harbors RAS mutations and the PAX8/PPARγ rearrangement can also be detected; both events are common in follicular tumors and are rare in classic PTC (5). On the other hand, it may have BRAF gene mutations—but, in contrast to the usual BRAF V600E found in classic PTC, FVPTC has a distinct BRAF K601E mutation. BRAF K601E is very rare in FTC and was only described in one case in previous publications (6,7). In addition, the overall gene expression profile in FVPTC is distinct from classic PTC, underlining these differences (8,9). In summary, FVPTC shares features of classic PTC and FTC on the morphological, genetic, and prognostic level.

MicroRNAs (miRNAs) are small endogenous, noncoding RNAs that regulate gene expression in practically every major cellular function, and as a consequence, deregulated miRNAs have been linked to a variety of different cancers throughout the human body, including thyroid carcinomas (10,11). Thyroid tumors encompass different miRNA profiles depending on the histological subtype (12,13). Some miRNAs like miR-221 or -222 have been reported in both classic PTC and FTC, while others like miR-146b and miR-31 have been reported in PTC only (14).

The goal of this study was to examine miRNA expression in a large series of FVPTC using the most comprehensive miRNA profile available and correlate it with the clinicopathological data.

Materials and Methods

Tissue samples and patient characteristics

For miRNA expression profiling, 52 neoplastic and non-neoplastic formalin-fixed paraffin embedded (FFPE) thyroid tissue samples were analyzed including 17 FVPTC, 27 classic PTC, and 8 normal thyroid tissues. The tissues were received from the University Hospital Zürich, Switzerland, and surrounding pathology institutes.

Two groups of patients were included in this study: (a) patients with an adverse clinical outcome (ACO) defined as recurrent tumor relapse or patient death, and (b) a control group of patients that had a normal (favorable) outcome without tumor relapse or death. For the ACO group, the records were searched at the Department of Nuclear Medicine, Canton Zürich, for thyroid carcinoma cases that had been operated between 1990 and 2006 and had one of the following characteristics: tumor relapse after first radioiodine therapy, distant metastases, or tumor-associated death. All patients from this group were matched by age, stage, and sex to the group of patients with a favorable outcome. Patient characteristics are summarized in Table 1. The study was approved by the Cantonal Research Ethics Board (STV 28-2006). All tumors were classified according to widely accepted diagnostic histologic criteria by two board certified pathologists, blinded to the clinicopathological outcome (M.D. and A.P.) (1). Validation of novel miRNAs was performed on 61 snap-frozen thyroid tissues (22 FVPTC, 24 PTC, 8 FTC, 2 hyperplastic nodules, and 5 normal tissues) from surgically removed thyroid samples collected at the Department of Pathology, University of Pittsburgh Medical Center following Institutional Review Board (IRB) approval.

Table 1.

Patient Characteristics

	Classic PTC	FVPTC	Total
Sex (n)
Female	18	10	28
Male	9	7	16
Age (years)
Years±SE Mean	47.8±3	50.9±4.3	49±3
pT
pT1	4	1	5
pT2	3	0	3
pT3	20	15	35
pT4	0	1	1
Status OS (n)
Censored	24	13	37
Death	3	4	7
t follow-up OS (months)	74.6±47	80.8±47.4	77±46.7
Status TSS (n)
Censored	25	16	41
Death	2	1	3
t follow-up TSS (months)	74.6±47	80.8±47.4	77±46.7
Status relapse (n)
Censored	8	8	16
Relapse	19	9	28
t follow-up RFS (months)	25.7±29.6	55.9±54	37.4±42.8
Total (n)	27	17	44

PTC, papillary thyroid carcinoma; FVPTC, follicular variant of papillary thyroid carcinoma; OS, overall survival; TSS, tumor-specific survival; RFS, relapse-free survival.

RNA isolation

RNA was isolated from FFPE tissue samples using the RecoverAll kit (Ambion, Life Technologies) according to the manufacturer's instructions. Each FFPE tissue specimen was stained with hematoxylin and eosin (H&E) to ensure that characteristic features of the thyroid tumor were present. Areas with high density and purity (>80%) of tumor cells were marked for microdissection of adjacent sections to minimize contamination from surrounding healthy thyroid tissue or infiltrating cells. Overlapping areas in up to six adjacent slides were manually microdissected from 15 μm unstained histological sections under the guidance of a H&E-stained slide using an Olympus SZ61 stereomicroscope (Olympus, Hamburg, Germany). RNA quality and quantity was assessed with a spectrophotometer (NanoDrop 1000; Thermo Scientific, Wilmington, DE). Total RNA was extracted from snap-frozen surgical specimens using Trizol reagent (Invitrogen, Life Technologies) as previously described (15).

miRNA expression analysis

Quantitation of mature miRNA expression levels in thyroid tumors and normal thyroid tissue was performed by RT-PCR using TaqMan^® Human Microarray Arrays v3.0 (Applied Biosystems, Life Technologies), which was designed to detect 754 human miRNAs. The array was performed on an ABI 7900 platform (Applied Biosystems, Life Technologies). First, 150 ng of total RNA was reverse transcribed using a high-capacity cDNA archive kit (Applied Biosystems, Life Technologies) and then followed by preamplification on an ABI 7500 Real-Time PCR System (Applied Biosystems, Life Technologies). The following PCR was run on an ABI 7900 Real-Time PCR System (Applied Biosystems, Life Technologies). Out of three different endogenous controls (RNU44, RNU48, U6 snRNA), the two most stable ones were used for normalization of RNA input: RNU44 and U6 snRNA. Nonhuman miRNA ath-miR159a was used as a negative control.

MiRNA expression levels were calculated by relative quantitation using DataAssist v3.0 software (Applied Biosystems, Life Technologies) and the fold-expression changes were determined by the 2^-ΔΔCt method (16). Outliers among replicates were excluded and p-values were adjusted using Benjamini–Hochberg false discovery rate. The data are presented as the fold change of miRNA expression in tumors relative to normal thyroid tissues after normalization to endogenous controls. Maximum allowed Ct values were 38. Expression of individual miRNAs was analyzed using the TaqMan^® individual miRNA assays (Applied Biosystems, Life Technologies) according to the manufacturer's instructions. A RNA input of 10 ng for high-capacity cDNA archive kit (Applied Biosystems, Life Technologies) was reverse transcribed followed by amplification on an ABI 7500 Real-Time PCR System (Applied Biosystems, Life Technologies).

Patient follow-up

Complete follow-up data were collected using chart reports and the cancer registry of the Canton Zürich and recorded as overall survival (OS), tumor-specific survival (TSS), and relapse-free survival (RFS). RFS was defined as the time between a patient's thyroidectomy and tumor relapse after the first radioiodine therapy.

Statistical analysis

Descriptive statistics, Mann–Whitney U-test, Kaplan–Meier analysis, and Cox regression were calculated with SPSS V21 (IBM, Armonk, NY). Fold changes were calculated with Dataassyst V3.1 (Applied Biosystems, Life Technologies), as was the unsupervised hierarchical clustering.

Results

Seventeen FVPTC, 27 classic PTC, and 8 normal thyroid tissues were analyzed for expression of 754 miRNAs using the most updated miRNA array panel, which covers the complete Sanger Database v.14. Overall, some miRNAs were similarly expressed in FVPTC and classic PTC. However, the levels of miRNA expression were different between these tumor types, and some miRNAs were uniquely dysregulated in FVPTC and classic PTC. The unsupervised hierarchical clustering analysis showed clear separation of the two entities from each other, demonstrating individual clusters for FVPTC, PTC, and normal thyroid tissue (Fig. 1), supporting their molecular and histopathologic differences.

FIG. 1.

Unsupervised hierarchical clustering (distance measure: Euklidian Distance; clustering method: complete linkage) of all reported miRNAs. Follicular variant of papillary thyroid carcinoma (FVPTC), classic papillary thyroid carcinoma (PTC), and normal thyroid tissue form three distinct clusters. Color images available online at www.liebertpub.com/thy

Comparison of miRNA expression between FVPTC and classic PTC

We identified 21 miRNAs significantly (p<0.05) deregulated (more than twofold up/down) in FVPTC as compared to normal tissue, with 14 of these deregulated more than fivefold up or down (Table 2). Strongly upregulated miRNAs in FVPTC were miR-146b-3p, -146-5p, -221, -222, and miR-222-5p, -181a, which were previously reported in PTC (14,17,18). Interestingly, the level of expression for miR-146b-3p and miR-146b-5p was similar between FVPTC and classic PTC (30–33-fold for miR-146b-3p and 73–77-fold for miR-146b-5p respectively). However, expression of miR-221 and miR-222-5p that are known to be upregulated in both PTC and FTC was twice as high in FVPTC compared to classic PTC.

Table 2.

Fold Changes and p-Values of Significant Deregulated miRNAs

miRNA	FVPTC vs. NT fold change >2 (p<0.05)	FVPTC vs. NT fold change >5 (p<0.05)	FVPTC vs. PTC fold change >5 (p<0.05)	FVPTC (RQ)	FVPTC (p-value)	Classic PTC (RQ)	Classic PTC(p-value)
hsa-let-7i-3p	×			3.2283	0.0138	1.9262	0.1875
hsa-miR-1179	×	×		0.0928	0.0269	0.0291	0.0007
hsa-miR-125a-3p			×	12.5594	0.075	1.1127	0.9453
hsa-miR-1271			×	2.8553	0.3094	0.4106	0.4034
hsa-miR-138	×	×		0.1564	0.0091	0.2647	0.0013
hsa-miR-144-5p	×			0.2472	0.0143	0.1392	0.0004
hsa-miR-146b-3p	×	×		29.4838	0.0143	32.8455	0.0138
hsa-miR-146b-5p	×	×		77.1832	0	73.2796	0.0002
hsa-miR-153			×	2.7257		0.1184
hsa-miR-15a-3p	×			3.8019	0.0138	3.3222	0.0148
hsa-miR-181a-2-3p	×	×		6.2228	0.0039	3.3402	0.0253
hsa-miR-204	×	×		0.1209	0.0053	0.058	0
hsa-miR-21	×			4.1654	0.0085	4.9069	0.0009
hsa-miR-219-5p	×	×		0.0372	0.0121	0.0185	0.0305
hsa-miR-221	×	×		16.1265	0	7.2353	0.0004
hsa-miR-222	×	×		16.2601	0.0003	13.5943	0.0008
hsa-miR-222-5p	×	×		9.7092	0.0138	3.396	0.1571
hsa-miR-34a	×			3.285	0.0137	3.05	0.0022
hsa-miR-34a-3p	×			2.5873	0.043	1.8738	0.0588
hsa-miR-375	×	×		37.094	0.0003	20.1274	0.0009
hsa-miR-451	×			0.3051	0.0252	0.1747	0.0004
hsa-miR-542-5p	×	×		7.2879	0.0185	2.1064	0.4096
hsa-miR-551b	×	×		22.5078	0.0138	15.7151	0.0166
hsa-miR-623			×	0.4794	0.7667	2.3888	0.6962
hsa-miR-99b-3p	×	×		7.0874	0.0509	1.3553	0.8164

NT, normal thyroid; ×, miRNAs with significant difference in expression between studied tissues (PTC, FVPTC, and NT).

We identified a number of downregulated miRNAs in both FVPTC and classic PTC, including miR-1179, -138, -144-5p, -199b-5p, -204, -219-5p, and miR-451 (Table 2). The levels of downregulation ranged between 3- and 54-fold. Four miRNAs were expressed significantly differently in FVPTC compared to classic PTC with a fold change more than five (miR-125a-3p, -1271, -153, and -623; Table 2).

In addition, two novel miRNAs miR-375 and miR-551b were found to be highly upregulated in FVPTC. MiR-375 showed more than a 37-fold upregulation as compared to normal thyroid tissue and miR-551b was upregulated more than 22-fold (Table 2). Analysis of their expression in classic PTC also revealed their upregulation but at a lower level (20-fold and 16-fold respectively).

Validation of miR-375 expression

MiR-375 was highly upregulated in FVPTC but at a lower level in classic PTC using the array approach. To confirm array data and to determine whether miR-375 upregulation is specific for papillary carcinoma, we studied its expression in an independent set of 61 follicular cell-derived neoplasms and normal thyroid tissue, including 22 FVPTC, 24 classic PTC, 8 FTC, 2 hyperplastic nodules and 5 normal tissues. In addition, we reanalyzed all arrayed FVPTC and classic PTC samples for miR-375 expression using individual RT-PCR reactions. MiR-375 was strongly upregulated in FVPTC and classic PTC but not in FTC, HN, or normal thyroid. A high correlation between array data and individual assays was observed (r=0.6).

Correlation of miRNA expression with outcome in FVPTC patients

The up/downregulated miRNAs identified in all FVPTC samples were correlated with the survival data. Patients with an adverse clinical outcome, that is, tumor relapse after first radioiodine therapy, distant metastases, or tumor-related death were compared to the control group. Because FVPTC are generally indolent neoplasms, overall survival (OS) is not a good clinical endpoint. Tumor-specific death (TSS) is a very rare event due to the low aggressiveness of these neoplasms, and our study did not reach statistical power for this outcome. Recurrent tumor relapse imposes a significant reduction in patient quality of life and therefore was used for outcome analysis.

Upregulation of miR-181a-2-3p and of miR-99b-3p and downregulation of miR-222 were associated with a favorable outcome in terms of tumor relapse in FVPTC (p<0.05). In the first two, Kaplan–Meier analysis reached statistical significance (p<0.05; Fig. 2). Both miRNAs were confirmed in two separate multivariate Cox regression analyses, including age, sex, and tumor stage [miR-181a-2-3p: p<0.05; Exp(B) 0.136; miR-99b-3p: Exp(B) 0.170], the only parameters that reached statistical significance.

FIG. 2.

Kaplan–Meier survival plots: relapse-free survival (RFS) is reduced for all three shown miRNAs, but only miR-181a-2-3p and miR-99b reach statistical significance (log rank p<0.002 and p<0.018).

Discussion

We assessed the miRNA expression in a large series of follicular variant and classic PTCs with the most comprehensive currently available miRNA coverage. We correlated miRNA expression with clinicopathological data and patient survival, revealing potential important prognostic markers for FVPTC. Finally, we reported a novel miRNA, miR-375, which is extremely consistent and highly upregulated in FVPTC and PTC but not in follicular carcinomas, hyperplastic nodules, and normal thyroid tissue.

MiRNAs have been demonstrated to play important roles in a variety of fundamental cellular processes, including cell proliferation, differentiation, and cell death (10,14). As a result, miRNAs are directly involved in the development and progression of cancer, representing an important regulatory mechanism. In this study, we assessed the miRNAome with the most comprehensive currently available coverage of Sanger miRBase v.14 in a large series of FVPTC samples.

There are several well-characterized miRNAs known to be upregulated (miR-146b, -21, -221, -222) (14,18,19) and downregulated (miR-1, miR-138) (18,20) in PTC. In the present work, we have confirmed these findings, which underline the robustness of this analysis and also strengthen the previously published work by ourselves and others. Interestingly, the level of expression for miR-146b-3p and miR-146b-5p was similar between FVPTC and classic PTC. However, expression of miR-221 and miR-222 that are known to be upregulated in both PTC and FTC was twice as high in FVPTC compared to classic PTC. Potentially, this may be due to the fact that FVPTCs share features of both classic PTC and FTC, and miR-221 and miR-222 are important in regulation of both tumor types. In addition, we previously demonstrated a correlation between miRNA expression and the mutational status of PTC (14). We showed that both miR-221 and miR-222 were upregulated more strongly in PTC with RAS mutations, which are commonly found in FVPTC. It would be potentially interesting to determine the prevalence of RAS mutations in this cohort and correlate it with miRNA expression.

We identified several new, never before reported miRNAs being significantly upregulated in FVPTC and classic PTC. MiR-375 was upregulated more than 35-fold in FVPTC and 20-fold in classic PTC, and not found to be upregulated in normal thyroid tissue, hyperplastic nodules, and follicular carcinomas. High-level and PTC-specific upregulation may indicate a potential diagnostic utility for this miRNA. MiR-375 is known to play a role in other human malignancies, including lung and breast cancer (12,21,22). Another miRNA found to be upregulated in both tumor types was miR-551b. However, this miRNA has never been shown to play a role in any human malignancy and therefore awaits confirmation in other studies and organ systems.

PTCs are known to be a heterogenous group of tumors, sharing in general a favorable outcome. Specific mutation profiles are found to be associated with distinct histological subtypes of PTC, such as the BRAF V600E mutations with a classic PTC phenotype, whereas RAS gene mutations and the BRAF K601E substitution are more frequently found in FVPTC (23,24). The present study demonstrates that this is also true for miRNA expression. Overall, the miRNA profile demonstrated similar trends in expression between FVPTC and classic PTC. However, some miRNAs were uniquely dysregulated in FVPTC and classic PTC. For example, miR-125a-3p, -1271, and -153 showed upregulation in FVPTC but were expressed at normal level or even downregulated in classic PTC, whereas the opposite pattern of expression was found for miR-623. These miRNAs have never been reported in thyroid cancer but are known to be dysregulated in other human malignancies, including esophageal, gastric, colon, hepatocellular, pancreas, breast cancer, neuroblastomas, and leukemias (21,25 –36). The unsupervised hierarchical clustering analysis showed clear separation of the two types of PTC from each other, demonstrating individual clusters for FVPTC and classic PTC, and supporting their molecular, histopathologic, and clinical differences.

FVPTC are known to have a favorable outcome. Identifying the small subset of patients who will suffer from recurrent disease is of major clinical interest. One of the important players in this scenario is the BRAF V600E mutation, which can be detected in about 40–45% of PTC (37) but at much lower frequency in FVPTC. Many studies demonstrated an association of this mutation with aggressive tumor behavior, although not all studies could confirm these findings (38). It is evident that a subset of tumors harboring the BRAF V600E mutation still has a good prognosis, underscoring the need for new molecular markers. Recently, a correlation between deregulation of miRNA expression in miR-146b, -222, and -130b and aggressiveness of PTC was reported (20,39). In this study, we were able to confirm upregulation of miR-222 and also identified novel markers of aggressiveness in FVPTC. Several miRNAs reached statistical significance when looking into RFS, including miR-181a-2-3p, -222, and -99b-3p. The Kaplan–Meier survival analysis was able to confirm a significant result for miR-181a-2-3p and miR-99b-3p. Patients had an adverse outcome if miR-222 was expressed at a higher level and miR-181a-2-3p and miR-99b-3p were expressed at a lower level and as compared to normal thyroid tissue. A multivariate analysis including age, sex, and tumor stage showed that the two latter miRNAs were the only significant parameters, predicting patient relapse. MiR-181a-2-3p has not been reported in human cancer before, but miR-99b-3p has been linked to prostate cancer (40,41). Although further studies are needed to to confirm this finding in a larger cohort of patients, these miRNAs bear the potential to serve as good molecular markers helping to stratify patients better in the future and facilitate guiding appropriate clinical management.

In conclusion, we report a comprehensive miRNA profile of FVPTC, the most common histologic type of PTC. We discovered a new miRNA (miR-375) that is upregulated in FVPTC and classic PTC, and identified several miRNAs that can predict patient relapse free survival. Further studies will have to prove their potential diagnostic and predictive clinical value.

Footnotes

Acknowledgments

M.D. was supported by the Fondation pour la recherche Nuovo-Soldati, the Gertrud-Hagmann-Stiftung für Malignomforschung, and the Research Support Foundation, and Y.E.N. was supported by the National Institute of Health grant R01 CA88041. We thank Prof. Steinert, Dr. Valenta, Dr. Haldemann, and Dr. Meili for providing us with clinicopathological data from their nuclear medicine departments. We also thank Dr. Flury, Cantonal Hospital Winterthur, Dr. Moll, Cantonal Hospital Münsterlingen, Prof. Diebold, Cantonal Hospital Luzern, Prof. Hochstetter, Institut Enge, Dr. Scheidegger, Viollier Basel, Dr. Baltisser, Institute Regenbogen, and Dr. Riehle, Institute Arnaboldi, for their collaboration. We also thank Alexiy Nikiforov for help with editing this manuscript and valuable comments on the submitted work.

Author Disclosure Statement

No competing financial interests exist.

References

DeLellis

, Lloyd

, Heitz

, Eng

. 2004. Pathology, Genetics of Tumours of Endocrine Organs. DeLellis

, Lloyd

, Heitz

, Eng

. IARC Press, Lyon [World Health Organization classification of tumours].

Nikiforov

YEB

, Paul

, Thompson

. 2009. Diagnostic Pathology and Molecular Genetics of the Thyroid. First edition Lippincott Williams & Wilkins: Philadelphia.

Chem

, Rosai

. 1977. Follicular variant of thyroid papillary carcinoma: a clinicopathologic study of six cases. Am J Surg Pathol, 1:123–130.

Chetty

. 2011. Follicular patterned lesions of the thyroid gland: a practical algorithmic approach. J Clin Pathol, 64:737–741.

Adeniran

, Zhu

, Gandhi

, Steward

, Fidler

, Giordano

, Biddinger

, Nikiforov

. 2006. Correlation between genetic alterations and microscopic features, clinical manifestations, and prognostic characteristics of thyroid papillary carcinomas. Am J Surg Pathol, 30:216–222.

Nikiforov

. 2011. Molecular analysis of thyroid tumors. Mod Pathol, 24:S34–43.

Pennelli

, Vianello

, Barollo

, Pezzani

, Merante Boschin

, Pelizzo

, Mantero

, Rugge

, Mian

. 2011. BRAF(K601E) mutation in a patient with a follicular thyroid carcinoma. Thyroid, 21:1393–1396.

Giordano

, Kuick

, Thomas

, Misek

, Vinco

, Sanders

, Zhu

, Ciampi

, Roh

, Shedden

, Gauger

, Doherty

, Thompson

, Hanash

, Koenig

, Nikiforov

. 2005. Molecular classification of papillary thyroid carcinoma: distinct BRAF, RAS, and RET/PTC mutation-specific gene expression profiles discovered by DNA microarray analysis. Oncogene, 24:6646–6656.

Igci

, Arslan

, Akarsu

, Erkilic

, Igci

, Oztuzcu

, Cengiz

, Gogebakan

, Cakmak

, Demiryurek

. 2011. Differential expression of a set of genes in follicular and classic variants of papillary thyroid carcinoma. Endocr Pathol, 22:86–96.

10.

, Getz

, Miska

, Alvarez-Saavedra

, Lamb

, Peck

, Sweet-Cordero

, Ebert

, Mak

, Ferrando

, Downing

, Jacks

, Horvitz

, Golub

. 2005. MicroRNA expression profiles classify human cancers. Nature, 435:834–838.

11.

Bartel

. 2004. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell, 116:281–297.

12.

Nikiforova

, Tseng

, Steward

, Diorio

, Nikiforov

. 2008. MicroRNA expression profiling of thyroid tumors: biological significance and diagnostic utility. J Clin Endocrinol Metab, 93:1600–1608.

13.

Dettmer

, Vogetseder

, Durso

, Moch

, Komminoth

, Perren

, Nikiforov

, Nikiforova

. 2013. MicroRNA expression array identifies novel diagnostic markers for conventional and oncocytic follicular thyroid carcinomas. J Clin Endocrinol Metab, 98:E1–E7.

14.

Nikiforova

, Chiosea

, Nikiforov

. 2009. MicroRNA expression profiles in thyroid tumors. Endocr Pathol, 20:85–91.

15.

Nikiforova

, Caudill

, Biddinger

, Nikiforov

. 2002. Prevalence of RET/PTC rearrangements in Hashimoto's thyroiditis and papillary thyroid carcinomas. Int J Surg Pathol, 10:15–22.

16.

Livak

, Schmittgen

. 2001. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods, 25:402–408.

17.

Pallante

, Visone

, Ferracin

, Ferraro

, Berlingieri

, Troncone

, Chiappetta

, Liu

, Santoro

, Negrini

, Croce

, Fusco

. 2006. MicroRNA deregulation in human thyroid papillary carcinomas. Endocrine-related cancer, 13:497–508.

18.

, Jazdzewski

, Li

, Liyanarachchi

, Nagy

, Volinia

, Calin

, Liu

, Franssila

, Suster

, Kloos

, Croce

, de la Chapelle

. 2005. The role of microRNA genes in papillary thyroid carcinoma. Proc Natl Acad Sci USA, 102:19075–19080.

19.

Tetzlaff

, Liu

, Xu

, Master

, Baldwin

, Tobias

, Livolsi

, Baloch

. 2007. Differential expression of miRNAs in papillary thyroid carcinoma compared to multinodular goiter using formalin fixed paraffin embedded tissues. Endocr Pathol, 18:163–173.

20.

Yip

, Kelly

, Shuai

, Armstrong

, Nikiforov

, Carty

, Nikiforova

. 2011. MicroRNA signature distinguishes the degree of aggressiveness of papillary thyroid carcinoma. Ann Surg Oncol, 18:2035–2041.

21.

Giricz

, Reynolds

, Ramnauth

, Liu

, Wang

, Stead

, Childs

, Rohan

, Shapiro

, Fineberg

, Kenny

, Loudig

. 2012. Hsa-miR-375 is differentially expressed during breast lobular neoplasia and promotes loss of mammary acinar polarity. J Pathol, 226:108–119.

22.

Nishikawa

, Osada

, Okazaki

, Arima

, Tomida

, Tatematsu

, Taguchi

, Shimada

, Yanagisawa

, Yatabe

, Toyokuni

, Sekido

, Takahashi

. 2011. miR-375 is activated by ASH1 and inhibits YAP1 in a lineage-dependent manner in lung cancer. Cancer Res, 71:6165–6173.

23.

Trovisco

, Soares

, Preto

, de Castro

, Lima

, Castro

, Maximo

, Botelho

, Moreira

, Meireles

, Magalhaes

, Abrosimov

, Cameselle-Teijeiro

, Sobrinho-Simoes

. 2005. Type and prevalence of BRAF mutations are closely associated with papillary thyroid carcinoma histotype and patients' age but not with tumour aggressiveness. Virchows Arch, 446:589–595.

24.

Adeniran

, Zhu

, Gandhi

, Steward

, Fidler

, Giordano

, Biddinger

, Nikiforov

. 2006. Correlation between genetic alterations and microscopic features, clinical manifestations, and prognostic characteristics of thyroid papillary carcinomas. Am J Surg Pathol, 30:216–222.

25.

Volinia

, Galasso

, Sana

, Wise

, Palatini

, Huebner

, Croce

. 2012. Breast cancer signatures for invasiveness and prognosis defined by deep sequencing of microRNA. Proc Natl Acad Sci USA, 109:3024–3029.

26.

Lam

, Wang

, Shin

, Zhang

, Morrison

, Sun

, Ng

, Yu

, Jin

. 2011. A microRNA contribution to aberrant Ras activation in gastric cancer. Am J Trans Res, 3:209–218.

27.

Yao

, Suo

, Li

, Liu

, Tian

, Ni

, Zhang

, Nan

, Song

, Huang

. 2009. MicroRNA profiling of human gastric cancer. Mol Med Report, 2:963–970.

28.

Favreau

, Sathyanarayana

. 2012. miR-590-5p, miR-219-5p, miR-15b and miR-628-5p are commonly regulated by IL-3, GM-CSF and G-CSF in acute myeloid leukemia. Leuk Res, 36:334–341.

29.

, Chang

, Meng

, Wang

, Xie

, Tang

, Li

, Song

, Lin

. 2012. MicroRNA-375 targets AEG-1 in hepatocellular carcinoma and suppresses liver cancer cell growth in vitro and in vivo. Oncogene, 31:3357–3369.

30.

Lin

, Chen

, Huang

, Gao

, Ye

, Song

, Shen

. 2011. MicroRNA expression profiles in human colorectal cancers with liver metastases. Oncol Rep, 25:739–747.

31.

Bray

, Tivnan

, Bryan

, Foley

, Watters

, Tracey

, Davidoff

, Stallings

. 2011. MicroRNA-542-5p as a novel tumor suppressor in neuroblastoma. Cancer Lett, 303:56–64.

32.

Tanic

, Zajac

, Gomez-Lopez

, Benitez

, Martinez-Delgado

. 2012. Integration of BRCA1-mediated miRNA and mRNA profiles reveals microRNA regulation of TRAF2 and NFkappaB pathway. Breast Cancer Res Treat, 134:41–51.

33.

Feber

, Xi

, Pennathur

, Gooding

, Bandla

, Wu

, Luketich

, Godfrey

, Litle

. 2011. MicroRNA prognostic signature for nodal metastases and survival in esophageal adenocarcinoma. Ann Thorac Surg, 91:1523–1530.

34.

Kalimutho

, Blanco

, Di Cecilia

, Sileri

, Cretella

, Pallone

, Federici

, Bernardini

. 2011. Differential expression of miR-144* as a novel fecal-based diagnostic marker for colorectal cancer. Journal of Gastroenterology, 46:1391–1402.

35.

Zhang

, Chen

, Lin

, Lwin

, Wright

, Moscinski

, Dalton

, Seto

, Wright

, Sotomayor

, Tao

. 2012. Myc represses miR-15a/miR-16-1 expression through recruitment of HDAC3 in mantle cell and other non-Hodgkin B-cell lymphomas. Oncogene, 31:3002–3008.

36.

Kovalchuk

, Filkowski

, Meservy

, Ilnytskyy

, Tryndyak

, Chekhun

, Pogribny

. 2008. Involvement of microRNA-451 in resistance of the MCF-7 breast cancer cells to chemotherapeutic drug doxorubicin. Mol Cancer Ther, 7:2152–2159.

37.

Xing

. 2005. BRAF mutation in thyroid cancer. Endocr Relat Cancer, 12:245–262.

38.

Xing

. 2007. BRAF mutation in papillary thyroid cancer: pathogenic role, molecular bases, and clinical implications. Endocr Rev, 28:742–762.

39.

Chou

, Chen

, Chou

, Chang

, Chen

, Lee

, Yang

, Cheng

, Huang

, Liu

. 2010. miR-146b is highly expressed in adult papillary thyroid carcinomas with high risk features including extrathyroidal invasion and the BRAF(V600E) mutation. Thyroid, 20:489–494.

40.

Sun

, Lee

, Malhotra

, Kim

, Matecic

, Evans

, Jensen

, Moskaluk

, Dutta

. 2011. miR-99 family of MicroRNAs suppresses the expression of prostate-specific antigen and prostate cancer cell proliferation. Cancer Res, 71:1313–1324.

41.

Maurel

, Jalvy

, Ladeiro

, Combe

, Vachet

, Sagliocco

, Bioulac-Sage

, Pitard

, Jacquemin-Sablon

, Zucman-Rossi

, Laloo

, Grosset

. 2013. A functional screening identifies five micrornas controlling glypican-3: role of mir-1271 down-regulation in hepatocellular carcinoma. Hepatology, 57:195–204.