Abstract
Background:
Thyroid cancer diagnosis in the United States has increased by 2.3-folds in the last three decades. Up to 30% of thyroid fine-needle aspiration biopsy (FNAB) results are inconclusive. Several differentially expressed microRNAs (miRNAs) have been identified as candidate diagnostic markers for thyroid nodules. We hypothesized that these differentially expressed miRNAs may improve the accuracy of FNAB in difficult to diagnose thyroid nodules.
Methods:
Expression levels of four miRNAs (miR-7, -126, -374a, and let-7g) were analyzed using quantitative real-time reverse transcription-polymerase chain reaction in 95 FNAB samples as the training set. A predictor model was formulated based on the most differentially expressed miRNA (miR-7) ΔCt value and the model was applied on a separate cohort of 59 FNAB samples as the validation set.
Results:
miR-7 was the best predictor to distinguish benign from malignant thyroid FNAB samples. The other three miRNAs were co-expressed and did not significantly contribute to the predictor model. miR-7 had a sensitivity of 100%, specificity of 29%, positive predictive value (PPV) of 36%, negative predictive value (NPV) of 100%, and overall accuracy of 76% when applied to the validation set. In subgroup analysis of preoperative nondiagnostic, indeterminate, or suspicious FNAB samples, the predictor model had an overall accuracy of 37% with sensitivity of 100%, specificity of 20%, PPV of 25%, and NPV of 100%.
Conclusions:
miR-7 may be a helpful adjunct marker to thyroid FNAB in tumor types which are inconclusive. Given the high NPV of miR-7, a patient with a benign result based on the predictor model may be followed as opposed to performing an immediate diagnostic thyroidectomy. Future prospective clinical trials evaluating its accuracy in a larger cohort are warranted to determine its clinical utility.
Introduction
There have been several factors which could influence the need for diagnostic thyroidectomy, such as repeat FNAB or secondary review of FNAB samples (8 –11). There has been a significant effort focused on combining the FNAB result with clinical information or combining the information from other diagnostic modalities such as ultrasonography, immunocytochemical markers, and molecular markers to improve the accuracy of thyroid FNAB (12 –19). Within the area of molecular markers there have been several studies evaluating the potential of small noncoding RNA called microRNAs (miRNAs) as diagnostic markers for thyroid cancer.
miRNAs are ∼22 nucleotides long and are responsible for gene regulation at the mRNA level (20). miRNA regulate gene expression by binding to the 3’untranslated region of their target mRNAs and can either lead to repression of translation or mRNA degradation depending on the level of the complementarity to their target mRNA sequence (21). There have been several studies which have investigated the miRNA signatures in thyroid cancer. Some of them have identified dysregulated miRNAs in papillary thyroid cancer compared with normal thyroid tissue or multinodular goiter (22 –26) and some studies have demonstrated the feasibility of testing these dysregulated miRNAs in FNAB samples (22,24,26,27). A study by Weber et al. addressed in particular the dysregulated miRNAs in follicular thyroid carcinoma (FTC) compared with follicular adenoma (FA) and found that two miRNAs (miR-197 and -346) were significantly overexpressed in FTC compared with their benign counterparts (28). Another study by Sheu et al. studied the diagnostic utility of miR-146b, -181b, -21, -221, and -222 in follicular variant of papillary thyroid cancer (FPTC) and FTC and they found that these five miRNAs can distinguish common variants of papillary thyroid cancer from FA/multinodular goiter (29). Several miRNA expression levels have also been implicated in thyroid tumorigenesis (24,28,30) and miR-146b in particular has been found to be associated with thyroid cancer with high risk features (31).
Recently our group found downregulation of miR-126 and miR-7 in thyroid cancer (32). We compared the two broad categories of benign and malignant tumor samples, which included the challenging histologic subtypes which are difficult to distinguish by FNAB (i.e., FA/carcinoma, Hürthle cell adenoma/carcinoma, and FPTC). Discovery of these miRNAs in difficult to diagnose histologic subtypes led to our current effort in validating this information in FNAB samples. Given the evidence suggesting suboptimal applicability of tissue expression data to FNAB samples, the goal of this study was to formulate a model that would predict benign from malignant samples based on a group of FNAB samples as the training set and validate this model on a separate cohort of FNAB samples as a validation set.
Materials and Methods
FNAB samples and patient information
FNAB samples, patient demographics, and histopathological information were obtained at the Johns Hopkins University School of Medicine, Baltimore, MD under an institutional review board approved protocol after written informed consent. All the FNA samples were obtained under ultrasound guidance, intraoperatively, before an incision was made. The samples were placed immediately into RNAlater® (Applied biosystems/Ambion, Foster City, CA) and stored at −80°C until they were ready for RNA extraction. The samples were classified based on The Bethesda System for Reporting Thyroid Cytopathology. A total of 95 FNAB samples were used as the training set to derive a predictor model and the model was applied to 59 FNAB samples as the validation set. The 59 FNAB validation samples were initially tested blinded without the cytological or histological diagnosis. Patient demographic and clinical information are summarized in Table 1.
Group 1 consisted of 93 patients with 95 nodules.
BSRTC, Bethesda System for Reporting Thyroid Cytology; TNM, tumor–node–metastasis; SD, standard deviation; NA, not available.
RNA isolation
RNA extraction of FNAB samples was performed using miRCURY™ RNA isolation kit (Exiqon, Woburn, MA) according to the manufacturer's protocol. RNA quantity and quality were determined using NanoDrop (Thermo Scientific, Wilmington, DE) and RNA was stored at −80°C until they were used for quantitative real-time reverse transcription-polymerase chain reaction (qRT-PCR).
miRNA expression analysis using qRT-PCR
qRT-PCR was performed on the 95 FNAB training cohort for miRNA-7 (assay #000268) and miRNA-126 (#002228) based on the model previously proposed by our group (32). Two additional miRNAs (miR-374a [assay #000563] and let-7g [assay #002282]), which were the third and fourth most significantly expressed miRNAs in our prior study, were also included in this analysis. qRT-PCR was performed on the 59 FNAB validation cohort for miRNA-7 (assay #000268) and miRNA-126 (#002228). TaqMan® microRNA assay probes were purchased from Applied Biosystems (Foster City, CA). cDNA was created using 5ng of total RNA and miRNA specific hairpin primers with TaqMan microRNA Reverse Transcription Kit (PN 4366597; Applied Biosystems). Five nanograms of total RNA were used in a 15 μL reaction to create the cDNA template. We used 2.5 μL of cDNA in a 10 μL reaction and all qRT-PCR reactions were performed in triplicates using miRNA specific primers (TaqMan microRNA assay) and TaqMan Universal PCR Master Mix on an ABI 7900HT Fast Real-Time PCR system (Applied Biosystems). miR-625* (assay #002432) was used as the endogenous control based on our previous data, which showed the least variability across thyroid tumor types (32). miRNA expression levels were calculated as the delta cycle threshold (ΔCt), which was calculated from cycle threshold (Ct) of the miRNA of interest subtracted by the Ct value of the endogenous control.
Statistical analysis
An initial screening procedure was performed on the ΔCt values of 95 FNAB training set samples in which the distribution of each of the miRNAs was compared between the two groups of interest (benign vs. malignant) using an exact Wilcoxon rank sum test. Any miRNAs that were found to have possible differences in distribution, as exhibited by having univariate p-values from the screening procedure of <0.05, were further individually evaluated and then jointly in a logistic regression model. All p-values were analyzed as two-tailed tests. The Mann–Whitney U test was used to compare miRNA expression levels between benign and malignant samples.
Results
Patient demographics and clinicopathological information
The clinical and demographical data for FNAB samples are summarized in Table 1. All the clinicopathological features were relatively similar between the training and the validation cohort, except there were slightly more patients with a history of radiation in the validation cohort. Final histology of FNAB samples was primarily papillary thyroid cancer and FPTC. The majority of malignant samples were comprised of American Joint Committee on Cancer tumor–node–metastasis classification system stage I tumors.
miRNA expression levels and predictor model derivation based on 95 FNAB training set
Expression levels of miR-7, -126, -374a, and let-7g were tested using qRT-PCR (Fig. 1). Only miR-7 was significantly differentially expressed between benign and malignant samples. miR-7 was significantly downregulated in malignant samples compared with benign with a p-value of <0.0001.

Scatter-plot of miRNAs [miR-7
Based on the exact Wilcoxon rank sum test and logistic regression model, miR-7 was the only miRNA associated with the ability to classify tumors as benign or malignant. The cutoff value for malignant classification was miR-7 ΔCt ≥−3.40456; classify benign if miR-7 ΔCt <−3.40456. This model had a sensitivity of 82%, specificity of 73%, positive predictive value (PPV) of 56%, negative predictive value (NPV) of 91%, and overall accuracy of 76%.
miRNA expression levels and result of the predictor model in the FNAB validation set
Expression levels of miR-7 and -126 were tested using qRT-PCR. Both miR-7 and miR-126 were significantly downregulated in malignant samples compared with benign ones with p-values of 0.0015 and 0.0076, respectively (Fig. 2).

Scatter-plot of the 2 most differentially expressed miRNAs [miR-7
When the predictor model derived from the FNAB training set was applied to the 59 samples, it had a sensitivity of 100%, specificity of 29%, PPV of 36%, NPV of 100%, and overall accuracy of 76%. When the model was specifically applied to the subgroup preoperatively designated as nondiagnostic, indeterminate, or suspicious it had a sensitivity of 100%, specificity of 20%, PPV of 25%, and NPV of 100%.
Discussion
In this study, we formulated a predictor model that would distinguish benign from malignant samples based on miRNA expression levels. We tested the expression levels of four reported miRNAs (miR-7, -126, -374a, and let-7g) that were most significantly differentially expressed based on tissue microarray and qRT-PCR (32) analysis in 95 FNAB training set samples. We then applied the model on a separate cohort of 59 FNAB samples as the validation set. The application of the model based on the 95 FNAB training set on the 59 FNAB validation set demonstrated a sensitivity of 100% with NPV of 100% but with suboptimal specificity and PPV. When the model was applied to the subgroup of patients with inconclusive FNAB findings, it retained high sensitivity and NPV of 100% but with lower specificity and PPV.
The introduction of routine use of FNAB has significantly reduced the number of unnecessary thyroidectomies for benign disease while increasing the yield of malignancy in thyroidectomy specimen (33 –35). However, indeterminate lesions which may comprise up to 30% of FNAB result still pose significant challenge to clinicians because many of these patients require a thyroidectomy for definitive histologic diagnosis. There have been numerous studies trying to improve the accuracy of FNAB to obviate the need for unnecessary thyroid surgeries and to allow for more complete initial thyroid operations. Among these efforts, several studies have focused on miRNA expression profiling of thyroid cancer. Within these studies, miR-222, -221, and -146b are consistently demonstrated to be upregulated in papillary thyroid cancer (22 –24). However, in our initial microarray analysis of difficult-to-diagnose thyroid nodules on FNA cytology, we did not find these miRs to be differentially expressed on quantitative RT-PCR validation of the microarray data (32). Numerous studies having identified differentially expressed miRNAs in thyroid cancer, specific miRNA(s) which would distinguish benign from malignant thyroid tumor in indeterminate FNAB samples have not yet been validated in clinical FNAB samples.
We found miR-7 to be the most accurate miRNA to distinguish benign from malignant thyroid FNAB samples. There currently is no available data on miR-7 in thyroid cancer or thyroid tumorigenesis. miR-7 has been characterized as a tumor suppressor in some malignancies [e.g., tongue squamous cell carcinoma (36) and glioblastoma (37)] but also has been characterized as an oncogenic miRNA “oncomiR” in lung cancer (38). It has also been associated with aggressiveness of lymph-node negative, estrogen-receptor positive human breast cancers (39). Further, miR-7 has been found to target several different proto-oncogenes such as epidermal growth factor receptor (EGFR) (40,41), insulin-like growth factor 1 receptor (IGF1R) (36), and p2-activated kinase 1 (Pak1) (42). There are some studies that suggest that miR-7 acts upstream of EGFR as its regulator (40,41). However, other investigators suggest that miR-7 is downstream of EGFR, suggesting that EGFR regulates miR-7 expression (38). Nonetheless, there appears to be an intricate relationship between miR-7 expression levels with EGFR and with other proto-oncogenes that will require further functional studies to understand its role in thyroid cancer initiation and/or progression.
Our previous study is the first to demonstrate the downregulation of miR-7 in thyroid cancer and this was demonstrated in both thyroid tissue and FNAB samples (32). To our knowledge, our study is also the first to evaluate the miRNA expression profile of the specific difficult to diagnose histologic subgroups on FNAB and validate these candidate diagnostic miRNAs in FNAB samples. Many patients undergo a diagnostic thyroidectomy for inconclusive thyroid FNAB results. Although the accuracy of our predictor model was suboptimal, it had a sensitivity and NPV of 100% in both the validation cohort and the subgroup of indeterminate FNAB samples. This implies that when the predictor model predicts negative for malignancy based on miR-7 measurement such patients may be closely followed rather than undergoing immediate diagnostic thyroidectomy. However, due to its low specificity and PPV, a positive result would not be accurate enough to make a decision on the appropriate extent of the initial thyroidectomy. We recognized additional studies with even larger sample size, showing similarly high sensitivity and NPV of 100%, will be necessary for clinical implementation of miR-7 expression level analysis as an adjunct to FNAB.
In summary, miR-7 was identified to be the most significantly differentially expressed miRNAs between benign and malignant thyroid FNAB samples and was the best predictor for distinguishing benign from malignant FNAB samples. Our predictor model for miR-7 has excellent NPV. Future prospective clinical trials should be considered and conducted in a larger cohort of FNAB samples to validate its clinical utility.
Footnotes
Acknowledgment
This research was supported by the Intramural Research Program, Center for Cancer Research, National Cancer Institute, National Institutes of Health.
Author Disclosure Statement
No competing financial interests exist.
