- Original research article
- Open Access
Serum microRNA profiling to distinguish papillary thyroid cancer from benign thyroid masses
Journal of Otolaryngology - Head & Neck Surgery volume 44, Article number: 33 (2015)
Papillary thyroid cancer (PTC) is increasing in incidence. Fine needle aspiration is the gold standard for diagnosis, but results can be indeterminate. Identifying tissue and serum biomarkers, like microRNA, is therefore desirable. We sought to identify miRNA that is differentially expressed in the serum of patients with PTC.
Serum miRNA was quantified in 31 female thyroidectomy patients: 13 with benign disease and 18 with PTC. qPCR results were compared for significant fold-changes in 175 miRNAs, against a pooled control.
128 miRNA qualified for analysis. There were identifiable fold-changes in miRNA levels between benign and control, and between PTC and control. There were statistically significant fold changes in the level of four miRNAs between benign and PTC: hsa-miR-146a-5p and hsa-miR-199b-3p were down-regulated, while hsa-let7b-5p and hsa-miR-10a-5p were up-regulated.
MicroRNA is differentially expressed in the serum of patients with PTC. Serum miRNA has the potential to aid in thyroid cancer diagnosis.
Papillary thyroid carcinoma (PTC) accounts for up to 90 % of thyroid carcinomas, yet distinguishing carcinoma from benign thyroid conditions such as follicular adenomas, cysts and goiter can be difficult . Thirty percent of fine needle aspirate specimens from thyroid nodules are indeterminate, and up to 75 % of patients undergo a hemithyroidectomy for benign disease unnecessarily [2–5]. As the incidence of PTC increases, along with the detection of incidental thyroid nodules on CT and ultrasound, research into the development of more sophisticated means of diagnosis is needed.
Recent research on the development of such diagnostic tools has focused on determining thyroid carcinoma gene expression profiles to identify markers of benign or malignant pathology . Such markers, traditionally tested in tumor tissue, have been used to stratify indeterminate fine needle aspiration results to those more likely to be PTC versus benign nodules. Examples of potential markers that appear to be consistently up-regulated in PTC tissue samples are chitinase 3 like-1 (YKL-40), galectin-3 (Gal-3), cytokeratin 19 (CK19), tissue inhibitor of metalloproteinases-1 (TIMP-1) and angiopoetin-1 (Ang −1) [7–10].
An alternative to tumor tissue sampling is serum-based diagnostic testing. Serum testing would potentially allow for minimally invasive, safe and repeatable diagnosis, requiring less technical skill for collection. Although the feasibility of potential blood-based diagnostic testing has been shown , no biomarkers have been identified that reliably reflect elevated tumor tissue levels.
MicroRNAs (miRNAs) are non-coding RNA molecules involved in the regulation of gene expression . miRNAs bind to messenger RNA (mRNA) and promote their degradation or prevent their translation into proteins. Over 2000 miRNA have been identified in the human genome, and more are being identified at a rapid rate . The first miRNA was incidentally identified in 1993; however, its function as a positive and negative bioregulator was not recognized until the early 2000s [14, 15]. Deregulation of miRNA has been associated with a number of diseases, including malignancies such as chronic lymphocytic leukemia, colon cancer, glioblastoma and astrocytoma [16–19]. Traditionally, miRNAs have been detected in tissue specimens; however, recent studies have also shown that miRNAs remains relatively stable and detectable in blood samples [20, 21].
Numerous independent studies have identified miRNA expression profiles in differentiated thyroid carcinoma tissue, suggesting they may play a role in the development and progression of this carcinoma. miRNA profiles identified in thyroid cancer have included miRNA-221, −222 and -181b [22, 23]. It is possible that these unique miRNA may allow for improved diagnosis, prognosis and monitoring response to treatment in thyroid carcinoma.
The objective of this study was to identify potential miRNA serum markers in patients with PTC. We hypothesized that, given the stability of miRNA in serum compared to traditional biomarkers, patients with PTC would show a measureable change in the select miRNA levels relative to patients with benign lesions.
Materials and methods
This study was approved by the Capital Health Research Ethics Board and informed consent was obtained from all participants.
Blood samples (approximately 15 mL) were collected in EDTA tubes from patients at the time of intravenous catheter placement in the operating room, prior to the administration of any anesthetic. A dedicated research nurse prepared serum samples from whole blood by centrifugation to separate serum from other blood components. Serum was aliquoted into 1.5 mL Eppendorf tubes, which were de-identified, barcoded and stored at −80 °C for future use. This process generated a biobank of 205 samples. From these, 45 samples (26 PTC and 19 benign) were randomly selected from female patients with thyroid masses greater than 1 cm. Separate RNA extractions were performed on three aliquots of a commercially available serum samples from a pool of twenty healthy individuals (Precision Biologics, Dartmouth, NS, Canada) for use as controls. These samples were screened for the presence of hemolyzed red cells using the NanoDrop™ 1000 spectrophotometer (NanoDrop Products, Wilmington, DE, USA) to detect optical density at 415 nm. Samples with evidence of hemolysis were excluded prior to further analysis, leaving 31 samples (13 benign, 18 PTC) (Table 1).
Total RNA was extracted from 200 μL serum using the miRNeasy® Mini Kit (Qiagen, Mississauga, ON, Canada) as described by the manufacturer. A mixture containing 1.25 μg/mL of MS2 bacteriophage RNA, and RNA spike-ins was added to the serum. RNA was purified on an RNeasy® Mini spin column (Qiagen), eluted in 50 μL of RNase-free water, and stored at −80 °C.
Sample quality control using spike-ins
To ensure that the quality of the input RNA was sufficiently high for effective reverse transcription (RT), cDNA synthesis and amplification, two types of RNA spike-ins were used. RNA isolation controls (UniSp2 and UniSp4) were added to the sample prior to purification and used to detect differences in extraction efficiency and the presence of inhibitors. The cDNA synthesis control (UniSp6) was added in the RT reaction to evaluate cDNA synthesis.
Reverse Transcription (RT) and RT-qPCR
For quality assessment, qPCR assays for miRNA-103, −191, −23a, −30c, and −451 were performed. These miRNAs are expressed at a predictable level in the majority of sample types and are used to evaluate baseline miRNA content in the samples. A second hemolysis test was performed at this time using two of these microRNA: miR-451, which is increased in hemolyzed red blood cells compared to miR-23a, which is relatively stable in serum and plasma. A difference in Ct values (miR-23a minus miR-451) <8.0 qualified samples for analysis. Negative controls excluding template from the RT reaction were performed and profiled with the samples. Total RNA (2 μL) was reverse transcribed in 10 μL reactions using the miRCURY LNA™ Universal RT microRNA PCR, Polyadenylation and cDNA synthesis kit (Exiqon, Woburn, MA, USA). Each RT was performed in duplicate, including an artificial RNA spike-in (UniSp6). cDNA was diluted 50-fold and assayed in 10 μL qPCR reactions according to the protocol for miRCURY LNA™ Universal RT microRNA PCR (Exiqon).
miRNAs were profiled using the microRNA Ready-to-Use PCR human serum/plasma panel containing 175 miRNAs plus controls (Exiqon). 5 μL RNA was reverse transcribed in 25 μL reactions and assayed as described above. Each microRNA was assayed once by qPCR and negative controls excluding template from the RT reaction were included. PCR was performed in a LightCycler® 480 Real-Time PCR System (Roche Applied Science, Laval, QC, Canada) in 384-well plates. The amplification curves were analyzed using the Roche LC software, both for determination of expression values (2nd derivative Ct method) and for melting curve analysis.
Data analysis and statistics
The raw data was extracted from the LightCycler 480 software. For quality control samples, an average Ct was calculated for the duplicate RTs, and evaluation of expression levels was performed based on raw Ct-values. The amplification efficiency was calculated using algorithms similar to the LinReg software . All assays were inspected for distinct melting curves to ensure that the Tm was within the known specifications for the assay.
Statistics (limma expression analysis, including t-tests and false-discovery rates) were done on the R platform using the HTqPCR package . Features where values were missing in more than 50 % of samples from any group (PTC, benign, control) and features where the negative controls were within 5 Ct values of the experimental Ct values were removed from the analysis. Data for microRNAs that were present in all 31 samples were normalized using norm rank-invariant normalization because it yielded the lowest standard deviations in the control replicates. For differential expression analysis, features with standard deviation >0.5 in the control group were removed. Random forest analysis was performed in the R platform using the random Forest package  and the ROCR package .
Gene targets for differentially expressed microRNAs were identified using a microRNA target prediction and annotation database (mirDB) using a prediction score of ≥80 (www.mirDB.org).
The consistent levels of all assays using RNA spike-ins (Uni-Sp2, Uni-Sp4 and Uni-Sp6) show that the RNA extraction as well as RT efficiency was similar for all samples, and that both RT and qPCR were successful. The baseline level, which is derived from the no-template control RT reaction, also shows that none of the samples contain inhibitors. Quality assessment of all positive and negative controls using qPCR assays for the commonly expressed miR-103, −191, −23a, −30c, and −451 showed Ct values between 22–24 for the positive control and ~40 for the negative control, indicating that the RT reaction worked well. The expression levels for all of the tested samples for these miRNAs were comparable to each other and to what is seen for other samples of similar types. Furthermore, the difference in Ct between miR-23a and miR-451 are 8.0 or lower in all samples, indicating minimal signs of hemolysis. Based on these quality control data we proceeded with the miR-panel profiling of all passed samples.
On average, 136 miRNAs were detected per sample and 128 miRNAs qualified for analysis, of which 64 were detected in all 31 samples profiled. After rank invariant normalization followed by Random Forest, t-Test and limma statistical analyses, two miRNAs (hsa-miR-150-5p and hsa-miR-342-3p) were down-regulated and two (hsa-let7b-5p and hsa-miR-191-5p) were up-regulated in benign vs. control samples (Fig. 1a; p < 0.001, FDR <5 %). For PTC vs. control samples, one additional microRNA (hsa-miR-146a-5p) was down-regulated and one additional microRNA was up-regulated (hsa-miR-93-5p) (Fig. 1b; p < 0.001, FDR <1 %). Importantly, between PTC and benign samples, there was statistically significant down-regulation of hsa-miR-146a-5p and hsa-miR-199b-3p, and up-regulation of hsa-let7b-5p and hsa-miR-10a-5p (Fig. 1c, p < 0.05), although the false discovery rate was high (50 %) in this comparison. Random Forest analysis produced a classifier consisting of 11 microRNAs (Fig. 2a), of which miRNAs hsa-miR-10a-5p, hsa-miR-146a-5p and hsa-miR-199b-3p were the most informative. Using this classifier, an Area Under the Curve (AUC) of 0.851 was obtained (Fig. 2b). Results of the limma and t-Test analysis are shown in Table 2.
We performed additional statistical testing by adjusting the p-values for multiple comparisons. The miRNA differences between control and benign samples were statistically significant (adjusted p < 0.05 for all miRNA), as were the comparison between PTC and control. When this additional statistical analysis was performed, the difference in miRNA between benign and malignant miRNA was no longer significant (p > 0.05).
By searching miRDB with the 11 microRNAs in the classifier, 186 unique gene targets were identified. These targets were mostly associated with metabolic process (GO:0008152, 100 targets), cellular process (GO:0009987, 70 targets) and cell communication (GO:0007154, 54 targets) (Fig. 3). Six microRNAs were down-regulated (hsa-miR-146a-5p, −199-3p, −376a-3p, 339-5p, 28-3p, and let-7d-3p), which would result in enhanced expression of the target genes. Five microRNAs were up-regulated (hsa-miR-375, −10b-5p, −10a-5p, −505-3p, and -let7b -5p), which would result in repression of the target genes.
This study identified that the serum of patients with benign thyroid tumors showed statistically significant down-regulation of miR-150-5p and miR-342-3p, and up-regulation of miR-let7b-5p and miR-191-5p compared to controls. The serum of patients with PTC also showed deregulation of these miRNA, with additional down-regulation of miR-146a-5p and up-regulation of miR-93-5p. Most importantly, there was statistically significant down-regulation of miR-146a-5p and miR-199b-3p and up-regulation of let7b-5p and miR-10a-5p when comparing PTC serum to that of patients with benign tumors using limma and t-tests. These are the standard tests used in miRNA literature. Unfortunately, these differences between miRNA between benign and PTC patients did not reach significance with the addition of adjusted p-value for multiple comparators.
The differences in expression levels between the control group and both the benign and PTC groups were the greatest, with average fold changes of 3.5 and 3.0 respectively. However, the average fold change between miRNAs that were differentially regulated between the patients with PTC and those with benign thyroid tumors was only 1.5. The practical impact of this modest difference is a relatively high false discovery rate of 50 %. In other words, two of the four miRNA are likely false positives. Therefore, this panel of miRNAs would require further validation in an independent cohort of patients.
Random Forest plots are used for classification and regression to create predictive models for the variable importance of large numbers of variables. A perfect predictor would have an area under the curve (AUC) of 1, whereas a predictor with an AUC of 0.5 has a predictive value no better than chance. The AUC in the Random Forest plot in PTC versus benign using a set of the 11 most informative miRNA was 0.851. These miRNA target genes were confirmed to be involved in metabolic and cellular processes and cell communication (Fig. 3).
Numerous studies have identified tissue miRNAs involved in PTC tumorigenesis. The majority of these compare PTC tissue to normal thyroid tissue – either from comparative controls or from adjacent “normal” thyroid tissue. He et al. were among the first to identify up-regulation of miRNA 146b, 221 and 222 in 20 patients with PTC . miRNA-146b is closely related to miRNA-146a, which was identified in our study. MiR-146a is a known tumor suppressor and its down-regulation appears to enhance tumorigenesis. It plays a role in the classical NF-κB signal transduction pathway in PTC by modulating expression of protein kinase C epsilon (PKCε) and Ras/Raf-1 signaling. It also targets TRAF6 (mirDB score 100) and IRAK1 (mirDB score 87), which form part of a negative feedback loop in NF-κB signaling . Pallate et al. also identified the down-regulation of miRNA-181b in PTC tissue, corroborating the results of He’s group . Tetzlaff et al. compared PTC tissue to multinodular goiter using formalin fixed paraffin embedded tissue for analysis rather than fresh frozen tissue, and identified 13 up-regulated and 26 down-regulated miRNAs, including those stated above .
The let7 family of miRNA was the second group to be identified during the early study of miRNA . Let7b-5p was up-regulated in our study, both in benign and malignant tumors, and exhibited a statistically significant increase when they were compared. It is often significantly up-regulated in various cancers and exerts its effects by enhancing expression of oncogenes such as Ras, c-Myc, cyclin D and HMG [33–35]. Elevated let-7b expression in PTC may have a similar effect on these genes [36–38].
The most significant limitation of our study was its small sample size (13 benign, 18 PTC and 3 controls), limiting our statistical power. This pilot study does have the potential to guide future research into these and other miRNAs to aid in the diagnosis of PTC.
Our study only included the PTC subtype of thyroid carcinoma. Although the diagnostic accuracy for anaplastic and medullary carcinoma is greater on cytopathology alone, follicular carcinoma still remains in the differential diagnosis in non-diagnostic FNA. A limited number of studies have published data on tissue miRNAs in FTC patients. The incidence of FTC is significantly less than that of PTC, therefore obtaining sufficient numbers for meaningful data analysis is more difficult [39, 40].
Although the application of tissue tumor markers such as miRNA and gene profiling may potentially improve diagnostic accuracy, these must be obtained by FNA biopsy and therefore remain invasive, unpleasant and inconvenient for patients. This study suggests miRNA in serum is stable enough for detection, supporting previous evidence to this effect. This would potentially provide clinicians with a novel, non-invasive method to improve diagnostic accuracy while avoiding discomfort to the patients. Serum miRNA for diagnosis would be particularly useful in patients with multinodular goiter. In these patients, there are often numerous thyroid nodules >15 mm and it can be difficult to select which nodules should be biopsied . Serum-based testing would allow thyroid malignancy to be identified much more simply in these patients. A previous study suggests Canadian wait times for thyroid surgery are inappropriately long [42, 43]. Preoperative risk stratification for malignancy will help manage the limited operative time available in our Canadian health care system to expedite care for those more likely to have malignancy.
Serum-based testing in cancer is gaining in popularity as its potential feasibility is recognized. By March 2015, more than 10 papers on this topic had been published in a diversity of tumor types. These include non-small cell lung carcinoma [44, 45], hepatocellular carcinoma , ovarian cancer , and others. Studies in thyroid cancer have been limited to date. A study by Yu et al. in 2012 was the first to examine serum miRNAs in PTC and identified three that were significantly up-regulated in Chinese patients with PTC relative to benign and control patients . This study was limited in that it only aimed to detect 5-fold-changes or more, and only looked for those miRNA that were up-regulated in sera from patients with PTC relative to those with benign disease. A second study in Caucasian patients identified a single miRNA as being significantly up-regulated and one down-regulated in PTC compared to benign sera . These targets, however, were different than the miRNAs identified in the study focusing on Chinese patients.
The previous studies have failed to tightly stratify patients based on clinical characteristics, which can have a significant effect on miRNA expression profiles, although this has not universally been found to be the case . Each study has identified different miRNAs that are differentially expressed in sera of patients with PTC. The populations have all been quite disparate from one another. We have chosen to include exclusively female patients with tumors >1 cm in order to more tightly control differences in baseline miRNA profiles. Nonetheless, the differences in miRNA expression identified between studies suggests a need for a large-scale, tightly controlled study which would serve as a baseline for studies in other populations with different clinical characteristics.
Clinical application of miRNAs may extend to disease monitoring and prognostication. The prognosis of patients with PTC is generally good, with a 10-year survival of 85 % . Some thyroid tumors, such as the tall cell variant, can present aggressively. miRNA detection may identify important molecular information that may determine tumor clinicopathological characteristics and predict aggressiveness of disease and patients at risk of poorer outcomes, and future studies are planned to stratify based on histological subtypes. A previous study has already shown promise in detecting serum miRNAs that may be increased in recurrent PTC . Future miRNA research may even aid in management of PTC patients, given the role of miRNA in cell proliferation, differentiation and apoptosis.
This study identifies four microRNAs that exhibit statistically significant changes when comparing the serum of patients with PTC to those with benign nodules. These results demonstrate the potential for serum-based miRNA assays to improve diagnostic accuracy of PTC.
Chappell H, Mery L, Pritwish D, Dryer D, Ellison L, MacIntyre M, Marrett L, Wier HK: Canadian Cancer Society’s steering committee on cancer statistics. In Book Canadian Cancer Society’s steering committee on cancer statistics (Editor ed.^eds.). City: Canadian Cancer Society; 2012
Hartl DM, Travagli JP. The updated american thyroid association guidelines for management of thyroid nodules and differentiated thyroid cancer: a surgical perspective. Thyroid. 2009;19(11):1149–51.
Wiseman SM, Baliski C, Irvine R, Anderson D, Wilkins G, Filipenko D, et al. Hemithyroidectomy: the optimal initial surgical approach for individuals undergoing surgery for a cytological diagnosis of follicular neoplasm. Ann Surg Oncol. 2006;13:425–32.
Isaac A, Jeffery CC, Seikaly H, Al-Marzouki H, Harris JR, O’Connell DA. Predictors of non-diagnostic cytology in surgeon-performed ultrasound guided fine needle aspiration of thyroid nodules. Isaac et al. J Otolaryngol Head Neck Surgery. 2014;43:48. http://www.journalotohns.com/content/43/1/48.
Williams BA, Bullock MJ, Trites JR, Taylor SM, Hart RD. Rates of thyroid malignancy by FNA diagnostic criteria. J Otolaryngol Head Neck Surgery. 2013;42:61. http://www.journalotohns.com/content/42/1/61.
Yip L. Molecular diagnostic testing and the indeterminate thyroid nodule. Curr Opin Oncol. 2014;26:8–13.
Park YJ, Kwak SH, Kim DC, Kim H, Choe G, Park Do J, et al. Diagnostic value of galectin-3, HBME-1, cytokeratin 19, high molecular weight cytokeratin, cyclin D1 and p27(kip1) in the differential diagnosis of thyroid nodules. J Korean Med Sci. 2007;22:621–8.
Jarzab B, Wiench M, Fujarewicz K, Simek K, Jarzab M, Oczko-Wojciechowska M, et al. Gene expression profile of papillary thyroid cancer: sources of variability and diagnostic implications. Cancer Res. 2005;65:1587–97.
Shi Y, Parhar RS, Zou M, Hammami MM, Akhtar M, Lum ZP, et al. Tissue inhibitor of metalloproteinases-1 (TIMP-1) mRNA is elevated in advanced stages of thyroid carcinoma. Br J Cancer. 1999;79:1234–9.
Hsueh C, Lin JD, Wu IC, Chao TC, Yu JS, Liou MJ, et al. Vascular endothelial growth factors and angiopoietins in presentations and prognosis of papillary thyroid carcinoma. J Surg Oncol. 2011;103:395–9.
Makki FM, Taylor SM, Shahnavaz A, Leslie A, Gallant J, Douglas S, et al. Serum biomarkers of papillary thyroid cancer. Otolaryngol Head Neck Surg. 2013;42:16.
He L, Hannon GJ. MicroRNAs: Small RNAs with big roles in gene regulation. Nat Rev Genet. 2004;5:522–31.
Friedlander MR, Lizano E, Houben AJS, Bezdan D, Báñez Coronel M, Kudla G, et al. Evidence for the biogenesis of more than 1000 novel human microRNAs. Genome Biol. 2014;15:R57. doi:10.1186/gb-2014-15-4-r57.
Lee RC, Feinbaum RL, Ambros V. The C. elegans heterochronic gene lin-4 encodes small RNAs with antisense complementarity to lin-14. Cell. 1993;75(5):843–54.
Cuellar TL, McManus MT. MicroRNAs and endocrine biology. J Endocrinol. 2005;187(3):327–32.
Mraz M, Pospisilova S. MicroRNAs in chronic lymphocytic leukemia: From causality to associations and back. Expert Rev Hematol. 2012;5(6):579–81.
He L, Thomson JM, Hemann MT, Hernando-Monge E, Mu D, Goodson S, et al. A microRNA polycistron as a potential human oncogene. Nature. 2005;435(7043):828–33.
Mraz M, Pospisilova S, Malinova K, Slapak I, Mayer J. MicroRNAs in chronic lymphocytic leukemia pathogenesis and disease subtypes. Leuk Lymphoma. 2009;50(3):506–9. doi:10.1080/10428190902763517.
Moller H, Rasmussen A, Anderson H, Johensen K, Henriksen M, Duroux M. A systemic review of microRNA in glioblastoma multiforme: micromodulators in the mesenchymal mode of migration and invasion. Mol Neurobio. 2013;47:131–44.
Blondal T, Jensby Nielsen S, Baker A, Andreasen D, Mouritzen P, Wrang Teilum M, et al. Assessing sample and miRNA profile quality in serum and plasma or other biofluids. Methods. 2013;59:81–6.
Cortez M, Calin G. MicroRNA identification in plasma and serum: a new tool to diagnose and monitor diseases. Expert Opin Biol Ther. 2009;9:703–11.
Li X, Abdel-Mageed A, Mondal D, Kandil E. MicroRNA expression profiles in differentiated thyroid cancer, a review. Int J Clin Exp Med. 2013;6:74–80.
Chapelle A, Jazdzewski K. MicroRNAs in thyroid cancer. J Clin Endocrinol Metabl. 2011;96:3326–36.
Ruijter JM, Ramakers C, Hoogaars WMH, Karlen Y, Bakker O, van den Hoff MJB, et al. Amplification efficiency: linking baseline and bias in the analysis of quantitative PCR data. Nucleic Acids Res. 2009;37(6):e45.
Dvinge H, Bertone P. HTqPCR: High-throughput analysis and visualization of quantitative real-time PCR data in R. Bioinformatics. 2009;25(24):3325.
Liaw A, Wiener M. Classification and regression by randomForest. R News. 2002;2(3):18–22.
Sing T, Sander O, Beerenwinkel N, Lengauer T. ROCR: visualizing classifier performance in R. Bioinformatics. 2005;21(20):7881.
He H, Jazdzewski K, Li W, Liyanarachchi S, Nagy R, Volinia S, et al. The role of microRNA genes in papillary thyroid carcinoma. Proc Natl Acad Sci U S A. 2005;102:19075–80.
Jazdzewski K, Liyanarachchi S, Swierniak M, Pachucki J, Ringel MD, Jarzab B, et al. Polymorphic mature microRNAs from passenger strand of pre-miR-146a contribute to thyroid cancer. Proc Natl Acad Sci U S A. 2009;106:1502–5.
Pallante P, Visone R, Ferracin M, Ferraro A, Berlingieri MT, Troncone G, et al. MicroRNA deregulation in human thyroid papillary carcinomas. Endocr Relat Cancer. 2006;13:497–508.
Tetzlaff MT, Liu A, Xu X, Master SR, Baldwin DA, Tobias JW, et al. Differential expression of miRNAs in papillary thyroid carcinoma compared to multinodular goiter using formalin fixed paraffin embedded tissues. Endocr Pathol. 2007;18:163–73.
Ambros V. MicroRNAs: tiny regulators with great potential. Cell. 2001;107(7):823–6.
Ricarte-Filho JC, Fuziwara CS, Yamashita AS, Rezende E, Da Silva MJ, Kimura ET. Effects of let-7 microRNA on cell growth and differentiation of papillary thyroid cancer. Transl Oncol. 2009;2(4):236–41.
Esquela-Kerscher A, Slack FJ. Oncomirs - microRNAs with a role in cancer. Nat Rev Cancer. 2006;6:259–69.
Ruzzo A, Canestrari E, Galluccio N, Santini D, Vincenzi B, Tonini G, et al. Role of KRAS let-7 LCS6 SNP in metastatic colorectal cancer patients. Ann Oncol. 2010;22(1):234–5.
Osada H, Takahashi T. Review Article: let-7 and miR-17-92: Small-sized major players in lung cancer development. Cancer Sci. 2010;102(1):9–17.
Nie K, Zhang T, Allawi H, Gomez M, Liu Y, Chadburn A, et al. Epigenetic down-regulation of the tumor suppressor gene PRDM1/blimp-1 in diffuse large B cell lymphomas : a potential role of the MicroRNA Let-7. Am J Pathol. 2010;177(3):1470–9.
Lee ST, Chu K, Oh HJ, Im WS, Lim JY, Kim SK, et al. Let-7 microRNA inhibits the proliferation of human glioblastoma cells. J Neurooncol. 2010;102(1):19–24.
Rossing M, Borup R, Henao R, Winther O, Vikesaa J, Niazi O, et al. Down-regulation of microRNAs controlling tumourigenic factors in follicular thyroid carcinoma. J Mol Endocrinol. 2012;48:11–23.
Weber F, Teresi RE, Broelsch CE, Frilling A, Eng C. A limited set of human MicroRNA is deregulated in follicular thyroid carcinoma. J Clin Endocrinol Metab. 2006;91:3584–91.
Luo J, Mcmanus C, Chen H, Sippel RS. Are there predictors of malignancy in patients with multinodular goiter? J Surg Res. 2012;174:209–10.
Brake MK, Moore P, Taylor SM, Trites J, Murray S, Hart R. Expectantly waiting: a survey of thyroid surgery wait times among Canadian Otolaryngologists. J Otolaryngol Head Neck Surgery. 2013;42:47. http://www.journalotohns.com/content/42/1/47.
Merdad M, Eskander A, De Almeida J, Freeman J, Rotstein L, Ezzat S, et al. Current management of papillary thyroid microcarcinoma in Canada. J Otolaryngol Head Neck Surgery. 2014;43:32. http://www.journalotohns.com/content/43/1/32.
Singh RK, Bethune DC, Xu Z, Douglas SE. Role of microRNAs in progression and recurrence of early-stage lung adenocarcinoma. Pulm Res Respir Med. 2015;2(1):52–62.
Wang RJ, Zheng YH, Wang P, Zhang JZ. Serum miR-125a-5p, miR-145 and miR-146a as diagnostic biomarkers in non-small cell lung cancer. Int J Clin Exp Pathol. 2015;8(1):765–71.
He S, Zhang DC, Wei C. MicroRNAs as biomarkers for hepatocellular carcinoma diagnosis and prognosis. Clin Res Hepatol Gastroenterol 2015; http://dx.doi.org/10.1016/j.clinre.2015.01.006.
Liang H, Jiang Z, Xie G, Lu Y. Serum microRNA-145 as a novel biomarker in human ovarian cancer. Tumor Biol. doi:10.1007/s13277-015-3191-y. Published online February 27, 2015.
Yu S, Liu Y, Wang J, Guo Z, Zhang Q, Yu F, et al. Circulating microRNA profiles as potential biomarkers for diagnosis of papillary thyroid carcinoma. J Clin Endocrinol Metab. 2012;97(6):2084–92.
Cantara S, Pilli T, Sebastiani G, Cevenini G, Busonero G, Cardinale S, et al. Circulating miRNA 95 and miRNA190 are sensitive markers for the differential diagnosis or thyroid nodules in a Caucasian population. J Clin Endocrinol Metab. 2014;99(11):4190–8.
Lee JC, Zhao JT, Clifton-Bligh RJ, Gill A, Gundara JS, Ip JC, et al. MicroRNA-222 and micro-RNA-146b are tissue and circulating biomarkers of recurrent papillary thyroid cancer. Cancer. 2013; doi:10.1002/cncr.28254.
Busaidy NL, Cabanillas ME. Differentiated thyroid cancer: management of patients with radioiodine nonresponsive disease. J Thyroid Res. 2012;2012:618985.
The authors would like to acknowledge Susanne Penny and Andrew Leslie, project bioinformaticians, for their work in statistical analysis. This project was funded by the Dalhousie Department of Surgery Seed Grant and the National Research Council of Canada (NRC).
The authors have no financial disclosure or conflicts of interest. This material has not previously been published and is not currently under evaluation in any other peer-reviewed publication.
MERG participated in sample preparation, specimen processing, assisted in bench-top analysis, data compilation and analysis, and drafted the manuscript. RDH was involved in study design, subject recruitment, specimen collection and processing and manuscript preparation. FM submitted the ethics proposal and assisted in specimen collection and processing. RS supervised miRNA analysis and assisted in study design, sample stratification and GQ, data analysis and preparation of the manuscript. DP was involved in study design, statistical analysis and manuscript preparation. AB assisted in manuscript preparation. MB was responsible for specimen preparation and was involved in study design. MHR assisted in statistical analysis and manuscript preparation. JRBT and SMT assisted in subject recruitment, specimen collection, and study design. SD was involved in study design and manuscript preparation. All authors read and approved the final manuscript.
Robert D. Hart and Rama Singh are co-senior authors on this manuscript