Whether Detection of Gene Mutations Could Identify Low- or High-Risk Papillary Thyroid Microcarcinoma? Data from 393 Cases Using the Next-Generation Sequencing

Jin, Lei; Zhou, Liang; Wang, Jian-Biao; Tao, Li; Lu, Xiao-Xiao; Yan, Na; Chen, Qian-Ming; Cao, Li-Ping; Xie, Lei

doi:https://doi.org/10.1155/2024/2470721

International Journal of Endocrinology

On this page

Abstract Introduction Methods Results Discussion Data Availability Disclosure Conflicts of Interest Authors’ Contributions Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2024 | Article ID 2470721 | https://doi.org/10.1155/2024/2470721

Whether Detection of Gene Mutations Could Identify Low- or High-Risk Papillary Thyroid Microcarcinoma? Data from 393 Cases Using the Next-Generation Sequencing

Lei Jin,¹Liang Zhou,¹Jian-Biao Wang,¹Li Tao,¹Xiao-Xiao Lu,¹Na Yan,²Qian-Ming Chen,³Li-Ping Cao,⁴and Lei Xie¹

Academic Editor: Claudio Casella

Received19 Jun 2023

Revised20 Nov 2023

Accepted30 Dec 2023

Published16 Jan 2024

Abstract

Objective. The objective of this study is to explore the utilization of next-generation sequencing (NGS) technology in evaluating the likelihood of identifying individuals with papillary thyroid microcarcinoma (PTMC ≤10 mm) who are at high or low risk. Design. NGS was used to analyze 393 formalin-fixed, paraffin-embedded tissues of PTC tumors, all of which were smaller than 15 mm. Results. The study found that bilateralism, multifocality, intrathyroidal spread, and extrathyroidal extension were present in 84 (21.4%), 153 (38.9%), 16 (4.1%), and 54 (13.7%) cases, respectively. Metastasis of cervical lymph nodes was identified in 226 (57.5%) cases and 96 (24.4%) cases with CLNM >5. Out of the total number of cases studied, 8 cases (2.3%) showed signs of tumor recurrence, all of which were localized and regional. Genetic alterations were detected in 342 cases (87.0%), with 336 cases revealing single mutations and 6 cases manifesting compound mutations. 332 cases (84.5%) had BRAF^V600E mutation, 2 cases had KRAS^Q61K mutation, 2 cases had NRAS^Q61R mutation, 8 cases had RET/PTC1 rearrangement, 3 cases had RET/PTC3 rearrangement, and 1 case had TERT promoter mutation. Additionally, six individuals harbored concurrent mutations in two genes. These mutations were of various types and combinations: BRAF^V600E and NRAS^Q61R (n = 2), BRAF^V600E and RET/PTC3 (n = 2), BRAF^V600E and RET/PTC1 (n = 1), and BRAF^V600E and TERT promoter (n = 1). The subsequent analysis did not uncover a significant distinction in the incidence of gene mutation or fusion between the cN0 and cN1 patient cohorts. The presence of BRAF^V600E mutation and CLNM incidence rates were found to be positively correlated with larger tumor size in PTMC. Our data showed that gene mutations did not appear to have much to do with high-risk papillary thyroid microcarcinoma (PTMC). However, when we looked at tumor size, we found that if the tumor was at least 5 millimeters in size, there was a higher chance of it being at high risk for PTM (, odds ratio (OR) = 2.55, 95% confidence interval (CI): 1.57–4.14). Identification of BRAF^V600E mutation was not demonstrated to be significantly correlated with advanced clinicopathological characteristics, although it was strongly associated with a bigger tumor diameter (OR = 4.92, 95% CI: 2.40–10.07, ). Conclusion. In clinical practice, BRAFV600E mutation does not consistently serve as an effective biomarker to distinguish high-risk PTMC or predict tumor progression. The size of the tumor has a significant correlation with its aggressive characteristics. PTMC with a diameter of ≤5 mm should be distinguished and targeted as a unique subset for specialized treatment.

1. Introduction

Papillary thyroid carcinoma (PTC) is a pathological subtype of thyroid cancer that is the most common of all endocrine-related cancers. It is known to have low-grade, inert biological characteristics. Papillary thyroid microcarcinoma (PTMC) is defined as PTC measured ≤10 mm in its greatest diameter. In recent years, there has been a notable upward trajectory in the prevalence of PTC, with approximately half of cases diagnosed as PTMC [1–3]. However, the majority of PTMC cases are actually occult and detected due to overscreening in healthy individuals. Based on the satisfactory results from Kuma Hospital [4, 5], active surveillance (AS) has been proposed for patients with low-risk PTMC [6]. However, how to reliably and safely identify low-risk PTMC patients has noticeably attracted thyroid surgeons’ attention.

Over the last two decades, substantial molecular genetic alteration studies have provided better insights into the understanding of the progression of PTC. Specific point mutations in BRAF or RAS, or the presence of RET chromosomal rearrangement, have the ability to activate the receptor tyrosine kinase (RTK)/mitogen-activated protein kinase (MAPK) pathway, a critical regulator of various cellular processes including proliferation, differentiation, adhesion, migration, and apoptosis [7–9]. PTC can progress to poorly or undifferentiated thyroid cancer by additional hits on TERT, tumor suppressors, or phosphoinositide 3-kinase (PI3K) pathway genes [8, 10–12]. In this study, we retrospectively analyzed the gene mutation status in 393 PTC patients using patients’ clinical and pathological characteristics by next-generation sequencing (NGS), in order to investigate the possibility of the identification of high- or low-risk PTMC patients.

2. Methods

2.1. Patients and Samples

A total of 393 formalin-fixed, paraffin-embedded (FFPE) tumor tissues from surgically removed thyroid samples were retrospectively analyzed by the NGS between August 2018 and July 2020. All of the tumors with a maximum diameter of less than 15 mm were confirmed to be PTC by pathological examination. The postoperative follow-up periods for a total of 393 patients spanned a range from 6 to 103 months, with a median duration of 20 months. Clinical records of these individuals who underwent surgical procedures between the timeframe of February 2013 and July 2020 were meticulously reviewed utilizing a comprehensive database retrospectively compiled by the Department of Head and Neck Surgery, the Affiliated Sir Run Run Shaw Hospital. The indication of thyroidectomy was strictly in accordance with the 2009 or 2015 American Thyroid Association (ATA) guideline [13, 14]. The indications of central neck dissection (CND) and modified lateral neck dissection (LND) were previously described [15]. This research project was endorsed by the Ethical Committee of the Sir Run Run Shaw Hospital Affiliated to Zhejiang University School of Medicine.

There were 300 cases with cN0, including 103 cases with a diameter of ≤5 mm, 151 cases with a diameter of 5∼10 mm, and 46 cases with a diameter of 10∼15 mm, as well as 93 cases with cN1 among 393 PTC patients. There were 320 patients with PTMC, including 218 high-risk and 102 low-risk PTMC patients (Figure 1).

2.2. Sample Preparation and Nucleic Acid Extraction

Following vendor manual instructions, genomic DNA and mRNA from FFPE samples were extracted using an AllPrep DNA/RNA FFPE tissue kit. This extraction kit is manufactured by Qiagen in Germany. Concentrations of DNA and cDNA were estimated using Qubit 3.0 from Thermo Fisher Scientific in the USA. In addition, DNA fragment distribution was assessed by a Qsep100 system from Bioptic, Taiwan, China.

2.3. NGS Library Preparation

The sequencing library should be prepared by the specific guidelines outlined in the U.S. Paragon Genomic DNA related to thyroid carcinoma was enriched by multiplex PCR, which was then connected with barcode and sequencing adapters. The amplified DNA libraries were purified using AMPure XP beads from Beckman Company of the United States. The Thermo Qubit 1 × dsDNA Assay Kit was a specific assay for quantifying next-generation sequencing libraries. The Qsep100 system developed by Biotic was used to analyze fragment size distribution. The well-known oncogenic drivers, including BRAF, RAS, TERT, and TP53, as well as gene rearrangements, such as RET and PAX-PPARG, were detected and are described in detail in Table 1.

2.4. Sequencing and Bioinformatics Analysis

The sequencing process was carried out on the NextSeq 500 sequencing platform from Illumina, using the PE150 sequencing strategy. Tumor samples demonstrated a mean coverage depth of approximately 1000X. The threshold for the initiation of allele mutations was set at 5%. Sequencing data were first converted to FASTQ format, and then, the human genome was mapped to the human reference sequence (hg19) using the Burrows-Wheeler aligner. Strelka was used to detect somatic insertions and deletions (InDel), while muTect was used to identify somatic single nucleotide variants (SNV). Genomic variant annotation was performed via the ANNOVAR 21 software.

2.5. Statistical Analysis

SPSS 20.0 software (IBM, Armonk, NY, USA) was employed for the execution of the statistical analysis. All continuous variables were reported using the mean value (standard deviation, SD) or median value (range). Categorical variables were represented as counts and corresponding percentages (%) to indicate the frequency of each category. In the event of data consistent with normal distribution and variance homogeneity, statistical analysis of differences between groups was conducted through the application of ANOVA. Mann–Whitney–Wilcoxon tests were used to evaluate differences among various groups, which were composed of data with skewed distributions. Categorical variables were assessed using either Pearson’s chi-square test for large samples or Fisher’s exact test for small sample sizes. To determine statistical significance, we adopted a two-sided significance threshold of across the distribution.

3. Results

3.1. The Clinical and Pathological and Molecular Features of Tumor Size Less than 15 mm in 393 PTC Patients

This study included 393 PTC cases with tumor diameter ≥15 mm. Detailed clinicopathological features are outlined in Table 2. In the present study, 274 women and 119 men were involved, with a mean age of 43.2 ± 12.0 years. Chronic lymphocytic thyroiditis, bilateralism, multifocality, intrathyroidal spreading, and extrathyroidal extension (ETE) were found in 71 (18.1%), 84 (21.4%), 153 (38.9%), 16 (4.1%), and 54 (13.7%) patients, respectively. Cervical lymph node metastasis (CLNM) was found in 226 (57.5%) patients, including central neck compartment in 157 (40.0%) cases, lateral neck compartment in 8 (2%) cases, and both compartments in 61 (15.5%) cases; however, CLNM more than 5 (>5 CLNMs) was found in 96 (24.4%) cases. According to the 8th edition of the American Joint Committee on Cancer (AJCC) TNM staging system, 358 (91.1%) patients were classified as stage I and 35 (8.9%) patients as stage II. During the follow-up, tumor recurrence was detected in 9 (2.3%) cases, and all of them were locoregional (2 cases in contralateral lobe and 7 cases in lateral neck compartmental lymph nodes). None of them had any history of undergoing radiation therapy for head and neck cancer before surgery and distant metastasis. In addition, 3 patients were lost to follow-up.

As presented in Table 2, 332 (84.5%) cases showed BRAF^V600E mutation, and RAS mutation was detected in 4 (1%) cases, including KRAS^Q61K and NRAS^Q61R in 2 cases equally. All of the point mutations were missense mutations. RET/PTC rearrangement was found in 11 (2.8%) cases, including RET/PTC1 in 8 cases and RET/PTC3 in 3 cases. TERT promoter gene mutation was detected in 1 (0.3%) case. No mutations in TP53 and HRAS genes were found, and PAX8/PPARG translocation was not detected.

Besides, there were 6 cases simultaneously carrying 2 gene mutations (fusions), including BRAF^V600E and NRAS^Q61R (n = 2), BRAF^V600E and RET/PTC3 (n = 2), BRAF^V600E and RET/PTC1 (n = 1), and BRAF^V600E and TERT promoter mutation (n = 1). No one had three or more gene mutations (fusions). However, it is noteworthy that 5 patients had older ages (>60 years old), and all of their PTCs had low-risk manifestations. The clinicopathological data of these 6 patients are listed in Table 3.

3.2. Comparison of Gene Mutation Status and Clinical Pathological Characteristics between cN0 and cN1 Patients with Tumor Diameter Less than 15 mm

Among 393 PTC patients, 93 (23.7%) patients showed clinical CLNM (cN1) before surgery. The 393 PTC patients whose tumor size was less than 15 mm were dispensed into two groups: cN0 group (n = 300) and cN1 (n = 93) group. Group cN1 was regarded as simulated clinical progression of observable PTC with cN0 (≤15 mm in size). Table 2 compares the clinicopathological features and genetic alterations between the cN0 and cN1 groups.

Compared with cN0 PTC, cN1 PTC showed more aggressive features, such as maximum tumor size, bilateralism, multifocality, intrathyroidal spreading, ETE, and LNM, except for TNM staging and tumor recurrence. Nevertheless, genetic mutations or fusions observed between cN0 and cN1 groups were not significantly different for BRAF^V600E, NRAS^Q61R, KRAS^Q61K, TERT promoter, and RET/PTC.

3.3. Correlation Analysis of Tumor Growth with Clinicopathological Features and Gene Mutation Status in cN0 PTC Patients with Tumor Size Less than 15 mm

In order to gain insight into the genetic mutation status during the progression of PTC, a total of 300 cN0 PTC cases were categorized into three parts based on primary tumor diameter: the first group (group A) consisted of 103 patients with tumor diameters of 5 mm or less, while the second group (group B) had a total of 151 patients with tumor diameters ranging from 5 mm to 10 mm. The third group (group C) contained 46 patients with tumor diameters ranging from 10 mm to 15 mm. Clinical progression of observable PTMCs was simulated by evaluating group C. The data are shown in Table 4.

As to the clinicopathological features, only CLNM was significantly different among the three groups (group A vs. group B, , odds ratio (OR) = 0.32; group B vs. group C, , OR = 0.49), including >5 CLNMs (group A vs. group B, , OR = 0.06; group B vs. group C, , OR = 0.65). Regarding genetic mutation, there was a significant difference only in BRAF^V600E mutation among the three groups (group A vs. group B, , OR = 0.23; group B vs. group C, , OR = 0.45). Therefore, the incidence rates of CLNM and BRAF^V600E mutation could increase during the PTMC growth.

3.4. Perform Univariate and Multivariate Logistic Regression Analyses on PTMC High-Risk Factors

All PTMC patients were divided into two groups: low-risk group (n = 102) and high-risk group (n = 218). It was found that BRAF^V600E, NRAS^Q61R, KRAS^Q61K, RET/PTC1, and RET/PTC3 showed no correlation with high-risk PTMC (Table 5). The multivariate logistic regression analysis has been further confirmed that gene mutations were not risk factors for high-risk PTMC in the present study; however, a significant correlation was identified between bigger tumor diameter and high-risk PTMC (, OR = 2.55, 95% confidence interval (CI): 1.57–4.14).

As per the ATA guidelines published in 2015, patients with more than 5 CLNMs are classified as the non-low-risk category. Therefore, the multivariate logistic regression analysis was carried out to appraise the integrated effect of histopathological and molecular characteristics on >5 CLNMs. It was revealed that gene mutations, including single and multiple mutations, were not significantly associated with >5 CLNMs, whereas male gender (, OR = 3.1, 95% CI: 1.60–5.86), older age at the time of diagnosis (40–60 years old, , OR = 3.73, 95% CI: 1.91–7.28; >60 years old, , OR = 17.89, 95% CI: 2.27–140.10), larger tumor size (, OR = 5.25, 95% CI: 2.29–12.04), intrathyroidal spreading (, OR = 10.62, 95% CI: 1.88–59.80), and bilateralism (, OR = 2.47, 95% CI: 1.18–5.18) were significantly associated with >5 CLNMs (Table 6).

3.5. The Association of the Presence of BRAF^V600E Mutation with Clinicopathological Features of PTMC Patients

As the prevalence of BRAF^V600E mutation was significantly higher than the other mutations in PTC patients, it was further attempted to analyze the association between BRAF^V600E mutation and clinicopathological features of PTMC patients. It was found that there was not a big connection between the presence of BRAF^V600E mutation and the more serious stages of the disease, such as multifocality, intrathyroidal spreading, ETE, LNM, and tumor recurrence, while it was significantly associated with the bigger tumor diameter (OR = 4.92, 95% CI: 2.40–10.07, ) (Table 7).

4. Discussion

In a 30-year longitudinal study of 140 cases, Woolner et al. established the concept of occult papillary carcinoma in 1960, referring to PTCs less than or equal to 15 mm in size [16]. Back in 1989, the World Health Organization created a new terminology to replace occult papillary carcinoma. They named it papillary microcarcinoma and decided it should be defined as papillary thyroid cancer that is less than or equal to 10 mm in diameter [17]. PTMC has been recognized as an essential variation of PTC due to exhibiting minimal malignant potential and rarely undergoing distant metastasis. These characteristics are typically found in autopsy studies at a frequency ranging from 6% to 35%. However, modern diagnostic methods have increased the detection rate of this type of cancer in the living. The best diameter of low-risk PTMCs for observation is usually between 10 mm and 15 mm. However, the Chinese Association of Thyroid Oncology (CATO) recommends that the dimensions of low-risk papillary thyroid microcarcinomas identified should not exceed 5 mm [18]. The perspective of CATO is derived from a comprehensive examination of two discriminating markers for AS in a sample of 1,001 low-risk patients diagnosed with PTMC. Compared to the findings in the Kuma low-risk PTMC group (lesions ≤10 mm), Qian et al. observed a significantly lower incidence of multifocal lesions, extrathyroidal extensions (ETEs), central LN metastasis, and progression rates in the CATO low-risk PTMC group (lesions ≤5 mm). Moreover, there was an observed improvement in the duration of disease-free survival. So we selected PTC patients with primary tumors that were less than or equal to 15 mm. We divided these patients into three groups based on their tumor size. The tumors that measured from 10 mm to 15 mm were presumed to represent a clinical progression of the visually noted PTMC.

From the study of the PTC in The Cancer Genome Atlas (TCGA) cohort that excluded poorly differentiated and anaplastic thyroid cancers, alteration of the MAPK pathway was detected in 83% of all the tested PTC samples, with the mutually exclusive activating mutations in BRAF (62%) and RAS (13%), as well as RET/PTC rearrangement (6%), and less frequently NTRK3 fusion (1.5%), NTRK1 fusion (1.3%), and ALK fusion (0.8%) [19]. Regarding PTMC, Rodrigues et al. reviewed 49 studies published from 1998 to 2015 and summarized the overall prevalence of BRAF^V600E, RAS, TERT promoter mutation, and RET/PTC rearrangement in 57.4% (5741/10004), 3.8% (4/106), 4.6% (19/409), and 44% (70/159) of cases, respectively [20]. In the present study, the prevalence of BRAF^V600E, RAS, TERT promoter mutation, and RET/PTC rearrangement in PTC samples with tumor size ≤15 mm was 84.5% (332/393), 1.0% (4/393), 0.3% (1/393), and 2.8% (11/393), respectively; however, in PTMC samples, no TERT promoter mutation was detected, and the prevalence of BRAF^V600E, RAS mutation, and RET/PTC rearrangement was 82.8% (265/320), 0.9% (3/320), and 2.5% (8/320), respectively. The difference in the spectrum of genetic alterations mentioned above was supposed to be associated with race, patient selection, the gene sequencing technique, the source of specimens, and the histologic subtype of PTC.

BRAF, RAS, and RET/PTC mutations are commonly recognized as mutually exclusive in PTC. This suggests that more than one mutation in these genes would likely not provide any additional biological benefits beyond those already conferred by a single mutation [21–23]. Zou et al. reported concomitant BRAF^V600E and KRAS^Q61R mutations in 2 classical PTCs, concomitant BRAF^V600E mutation and RET/PTC1 rearrangement in 6 classical PTCs, and 1 tall-cell variant of PTC as well as concomitant NRAS^Q61R mutation and RET/PTC1 rearrangement in 2 follicular variants of PTCs, and they demonstrated that concomitant gene mutation might be the cause of intratumor heterogeneity. Furthermore, they found that the concomitant mutations were associated with the advanced stage of the disease [24]. In the present study, there were 6 PTC cases simultaneously carrying 2 gene mutations, including BRAF^V600E and NRAS^Q61R (n = 2), BRAF^V600E and RET/PTC3 (n = 2), BRAF^V600E and RET/PTC1 (n = 1), and BRAF^V600E and TERT promoter mutation (n = 1). However, such cases showed no noticeable aggressive tumor behaviors. Among these 6 patients, it was found that 5 patients had older age (> 60 years old); thus, we hypothesized that with the increase of age, the effects of carcinogenic factors on the body were enhanced, which might lead to multiple mutations or fusion of somatic genes.

It is recommended at Kuma Hospital that low-risk patients with PTMC undergo active surveillance (AS) as the primary treatment [6]. In the event that a tumor exhibits an increase in size of 3 mm or more, a tumor measuring 12 mm in size, or the presence of lymph node metastases during the follow-up period, prompt consideration should be given to the implementation of “rescue surgery.” What role does gene mutation play in the course of disease progression? Whether gene mutation could predict the progression of PTMC? In light of these observations, we employed our dataset to model the progression of clinical symptoms and designed two distinct segments. In the first part, it was attempted to simulate LNM using two groups of cN0 and cN1, according to whether LNM was confirmed preoperatively. Expectedly, more aggressive features of PTMC were found in the cN1 group, while gene mutation showed no significant difference between the two groups. In the second part, cN0 patients were divided into three groups: ≤5 mm, 5∼10 mm, and 10∼15 mm according to the tumor size, in order to simulate tumor growth. It was found that the incidence of BRAF^V600E mutation significantly increased during the growth of PTMC, suggesting that BRAF^V600E mutation might promote tumor growth. The mechanism may include three aspects: (1) the BRAF^V600E mutation leads to sustained activation of the RAS-RAF-MEK-ERK/MAPK pathway, enhancing cell mitotic capacity, ultimately leading to abnormal cell proliferation, and inducing tumor occurrence [10]. (2) The BRAF^V600E mutation promotes methylation of the thyrotropin receptor (TSHR) gene promoter, leading to TSHR silencing or a significant decrease in expression, increasing the feedback of thyrotropin (TSH), and promoting tumor cell growth [25]. (3) The BRAF^V600E mutation upregulates TERT expression in a non-TERT promoter-dependent manner through the oncogene MYC, which also promotes the occurrence of PTC [26]. However, in clinical practice, BRAFV600E mutation could be used as a biomarker for identifying thyroid nodule malignancy, rather than predicting tumor growth.

Multiple studies have documented a favorable prognosis among patients with PTMC [27–30], with results indicating a clinical course that is generally favorable. This favorable prognosis has been corroborated by prospective clinical studies, in which patients with PTMC who chose not to undergo surgical intervention received AS primarily [31, 32]. Zaid et al. reported a series of 30,180 adult patients with PTMC. Notably, 5,628 patients (18.7%) showed advanced features, including central lymph node metastasis (8.0%), lateral lymph node metastasis (4.4%), microscopic extrathyroidal extension (6.7%), macroscopic extrathyroidal extension (0.3%), lymphovascular invasion (4.4%), and distant metastasis (0.4%). They described PTMC as a wolf in sheep’s clothing and suggested that lobectomy seemed to be a more appropriate treatment for low-risk PTMC patients rather than AS [33]. Whether gene detection can identify low- or high-risk PTMC patients preoperatively? PTMC was dispensed into two groups in our study: the high-risk group, in which PTMC had some advanced features, including ETE, LNM, multifocality, and intrathyroidal spreading, and the low-risk group (without any aggressive features). As a result, it was revealed that the larger the tumor size, the more aggressive the characteristics, and only tumor size could affect the advanced features of PTMC, rather than gene mutation, especially BRAF^V600E.

Although the primary objective of this study was to explore the correlation between gene mutations and the aggressiveness of the PTMC, it was clearly found that the tumor size, rather than gene mutation, was significantly associated with the progression of high-risk PTMC. For instance, when the tumor grew up from 5 to 10 mm, the incidence of LNM was doubled, and the nonmutation rate of BRAF^V600E was tripled. The findings of the present study were also similar to the previously reported results. After analysis of 977 PTMC cases, Zheng et al. found that tumor size (>5 mm) was an independent risk factor associated with LNM [34]. In a study looking at two screening criteria for AS, Qian et al. found that low-risk patients with less than 5 mm of PTMC had a lower chance of having multifocal lesions, ETE, central LNM, or disease progression. This means they were also less likely to relapse and had longer disease-free survival times [18]. Therefore, a tumor with a size of less than 5 mm should be regarded as a special subgroup of PTMC that showed to be less aggressive than the other PTMCs and might be treated specially with local ablation, such as radiofrequency, in order to block the tumor progression during the follow-up of PTMC.

In conclusion, our findings indicate that BRAFV600E mutation may not be a reliable biomarker to clinically distinguish high-risk PTMC patients or predict tumor progression. Tumor size was found to be highly correlated with the disease’s aggressiveness. PTMC tumors measuring ≤5 mm should be identified as a distinct subgroup and treated with specific treatment strategies.

4.1. Limitation

The present study had several limitations. First, the 393 patients in this cohort were analyzed inconsecutive from 2013 to 2020 because we were unable to obtain consent from all patients for NGS of their samples. The other reason was that not all patient samples were of suitable quality for NGS. Therefore, a selective bias might exist in the study. Second, the relatively short follow-up time may have an impact on the long-term recurrence rate. In our study, tumor recurrence was detected in 9 (2.3%) cases, and it may be better than other studies, but we believe the data are true. From January 2013 to December 2015, our medical center treated a total of 562 patients with papillary thyroid cancer (PTC) with initial surgical therapy during our prior research endeavor, which consisted of either unilateral lobectomy or total thyroidectomy, and central neck dissection. Some patients also underwent lateral neck dissection. Five recurrence cases were identified between a minimum follow-up period of 40 months and a maximum follow-up period of 75 months, resulting in a median follow-up period of 57 months [35]. We speculate that the recurrence rate is lower than others because all the patients underwent prophylactic or therapeutic cervical lymph node dissection, and the surgical team is very experienced and well-practiced. Third, samples for the NGS were from formalin-fixed, paraffin-embedded tumor tissues and histologic subtypes of PTC were not further classified; therefore, the map of gene mutation might be slightly different from other relative studies. Lastly, this NGS panel was small, but it included 9 hot spot gene mutations and rearrangements related to thyroid cancer and is a cost-effective auxiliary diagnostic test for thyroid nodules. The other PTC-associated drivers, such as EIF1AX, AKT1, and PIK3CA, will be considered and measured in future research.

Data Availability

The raw data reported in this study can be accessed and downloaded from the CNGB Nucleotide Sequence Archive (CNSA) (https://db.cngb.org/cnsa/) with accession no. CNP0004004.

Disclosure

The abstract of this article has been presented at the 2022 Annual Meeting of the American Thyroid Association, October 19–23, 2022, in Montréal, Quebec, Canada (Meeting Program and Abstracts. Thyroid: official journal of the American Thyroid Association 2022; 32(S1):P-1-A-135).

Conflicts of Interest

Na Yan is an esteemed member of the Key Laboratory of Digital Technology in Medical Diagnostics of Zhejiang Province, Dian Diagnostics Group Co., Ltd. The remaining authors solemnly state that their research was undertaken devoid of any commercial or financial influence.

Authors’ Contributions

Lei Jin, Liang Zhou, Jianbiao Wang, and Lei Xie conceptualized the article. Na Yan and Xiaoxiao Lu utilized NGS. Lei Jin, Liang Zhou, Li Tao, and Xiaoxiao Lu collected data and performed analysis. Jin Lei, Jianbiao Wang, Xiaoxiao Lu, and Lei Xie validated the article. Jin Lei and Lei Xie wrote the original draft and reviewed and edited the article.

Acknowledgments

This research was financially supported by a grant from the Bureau of Zhejiang Provincial Clinical Research Center (grant number 2022-KFKT-05).

References

L. Davies and H. G. Welch, “Increasing incidence of thyroid cancer in the United States, 1973-2002,” JAMA, vol. 295, no. 18, pp. 2164–2167, 2006.
View at: Publisher Site | Google Scholar
H. Lim, S. S. Devesa, J. A. Sosa, D. Check, and C. M. Kitahara, “Trends in thyroid cancer incidence and mortality in the United States, 1974-2013,” JAMA, vol. 317, no. 13, pp. 1338–1348, 2017.
View at: Publisher Site | Google Scholar
L. Furuya-Kanamori, K. J. L. Bell, J. Clark, P. Glasziou, and S. A. R. Doi, “Prevalence of differentiated thyroid cancer in autopsy studies over six decades: a meta-analysis,” Journal of Clinical Oncology, vol. 34, no. 30, pp. 3672–3679, 2016.
View at: Publisher Site | Google Scholar
Y. Ito, A. Miyauchi, M. Kihara, T. Higashiyama, K. Kobayashi, and A. Miya, “Patient age is significantly related to the progression of papillary microcarcinoma of the thyroid under observation,” Thyroid, vol. 24, no. 1, pp. 27–34, 2014.
View at: Publisher Site | Google Scholar
A. Miyauchi, Y. Ito, and H. Oda, “Insights into the management of papillary microcarcinoma of the thyroid,” Thyroid, vol. 28, no. 1, pp. 23–31, 2018.
View at: Publisher Site | Google Scholar
Y. Ito and A. Miyauchi, “Active surveillance of low-risk papillary thyroid microcarcinomas,” Gland Surgery, vol. 9, no. 5, pp. 1663–1673, 2020.
View at: Publisher Site | Google Scholar
M. R. Haroon Al Rasheed and B. Xu, “Molecular alterations in thyroid carcinoma,” Surgical Pathology Clinics, vol. 12, no. 4, pp. 921–930, 2019.
View at: Publisher Site | Google Scholar
I. Landa, T. Ibrahimpasic, L. Boucai et al., “Genomic and transcriptomic hallmarks of poorly differentiated and anaplastic thyroid cancers,” Journal of Clinical Investigation, vol. 126, no. 3, pp. 1052–1066, 2016.
View at: Publisher Site | Google Scholar
H. Duan, Y. Li, P. Hu et al., “Mutational profiling of poorly differentiated and anaplastic thyroid carcinoma by the use of targeted next-generation sequencing,” Histopathology, vol. 75, no. 6, pp. 890–899, 2019.
View at: Publisher Site | Google Scholar
J. A. Fagin and S. A. Wells Jr., “Biologic and clinical perspectives on thyroid cancer,” New England Journal of Medicine, vol. 375, no. 11, pp. 1054–1067, 2016.
View at: Publisher Site | Google Scholar
R. Liu and M. Xing, “TERT promoter mutations in thyroid cancer,” Endocrine-Related Cancer, vol. 23, no. 3, pp. R143–R155, 2016.
View at: Publisher Site | Google Scholar
J. M. Oh and B. C. Ahn, “Molecular mechanisms of radioactive iodine refractoriness in differentiated thyroid cancer: impaired sodium iodide symporter (NIS) expression owing to altered signaling pathway activity and intracellular localization of NIS,” Theranostics, vol. 11, no. 13, pp. 6251–6277, 2021.
View at: Publisher Site | Google Scholar
D. S. Cooper, G. M. Doherty, B. R. Haugen et al., “Revised American Thyroid Association management guidelines for patients with thyroid nodules and differentiated thyroid cancer,” Thyroid, vol. 19, no. 11, pp. 1167–1214, 2009.
View at: Publisher Site | Google Scholar
B. R. Haugen, E. K. Alexander, K. C. Bible 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, vol. 26, no. 1, pp. 1–133, 2016.
View at: Publisher Site | Google Scholar
J. B. Wang, K. Wu, L. H. Shi, Y. Y. Sun, F. B. Li, and L. Xie, “In situ preservation of the inferior parathyroid gland during central neck dissection for papillary thyroid carcinoma,” British Journal of Surgery, vol. 104, no. 11, pp. 1514–1522, 2017.
View at: Publisher Site | Google Scholar
L. B. Woolner, M. L. Lemmon, O. H. Beahrs, B. M. Black, and F. R. Keating Jr., “Occult papillary carcinoma of the thyroid gland: a study of 140 cases observed in a 30-year period,” Journal of Clinical Endocrinology and Metabolism, vol. 20, no. 1, pp. 89–105, 1960.
View at: Publisher Site | Google Scholar
C. Hedinger, E. Dillwyn Williams, and L. H. Sobin, “The WHO histological classification of thyroid tumors: a commentary on the second edition,” Cancer, vol. 63, no. 5, pp. 908–911, 1989.
View at: Publisher Site | Google Scholar
K. Qian, K. Guo, X. Zheng et al., “Contrastive study of two screening criteria for active surveillance in patients with low-risk papillary thyroid microcarcinoma: a retrospective analysis of 1001 patients,” Oncotarget, vol. 8, no. 39, pp. 65836–65846, 2017.
View at: Publisher Site | Google Scholar
N. Agrawal, R. Akbani, B. Aksoy et al., “Integrated genomic characterization of papillary thyroid carcinoma,” Cell, vol. 159, no. 3, pp. 676–690, 2014.
View at: Publisher Site | Google Scholar
A. C. Rodrigues, G. Penna, E. Rodrigues, P. Castro, M. Sobrinho-Simões, and P. Soares, “The genetics of papillary microcarcinomas of the thyroid: diagnostic and prognostic implications,” Current Genomics, vol. 18, no. 3, pp. 244–254, 2017.
View at: Publisher Site | Google Scholar
E. T. Kimura, M. N. Nikiforova, Z. Zhu, J. A. Knauf, Y. E. Nikiforov, and J. A. Fagin, “High prevalence of BRAF mutations in thyroid cancer: genetic evidence for constitutive activation of the RET/PTC-RAS-BRAF signaling pathway in papillary thyroid carcinoma,” Cancer Research, vol. 63, no. 7, pp. 1454–1457, 2003.
View at: Google Scholar
M. Xing, “BRAF mutation in papillary thyroid cancer: pathogenic role, molecular bases, and clinical implications,” Endocrine Reviews, vol. 28, no. 7, pp. 742–762, 2007.
View at: Publisher Site | Google Scholar
R. M. Melillo, M. D. Castellone, V. Guarino et al., “The RET/PTC-RAS-BRAF linear signaling cascade mediates the motile and mitogenic phenotype of thyroid cancer cells,” Journal of Clinical Investigation, vol. 115, no. 4, pp. 1068–1081, 2005.
View at: Publisher Site | Google Scholar
M. Zou, E. Y. Baitei, A. S. Alzahrani et al., “Concomitant RAS, RET/PTC, or BRAF mutations in advanced stage of papillary thyroid carcinoma,” Thyroid, vol. 24, no. 8, pp. 1256–1266, 2014.
View at: Publisher Site | Google Scholar
F. Orim, A. Bychkov, M. Shimamura et al., “Thyrotropin signaling confers more aggressive features with higher genomic instability on BRAF(V600E)-induced thyroid tumors in a mouse model,” Thyroid, vol. 24, no. 3, pp. 502–510, 2014.
View at: Publisher Site | Google Scholar
FI Moulana, AAH Priyani, MVC de Silva, and RS Dassanayake, “BRAF-oncogene-induced senescence and the role of thyroid-stimulating hormone signaling in the progression of papillary thyroid carcinoma,” Horm Cancer, vol. 9, no. 1, pp. 1–11, 2018.
View at: Publisher Site | Google Scholar
Y. W. Jeon, H. G. Gwak, S. T. Lim, J. Schneider, and Y. J. Suh, “Long-term prognosis of unilateral and multifocal papillary thyroid microcarcinoma after unilateral lobectomy versus total thyroidectomy,” Annals of Surgical Oncology, vol. 26, no. 9, pp. 2952–2958, 2019.
View at: Publisher Site | Google Scholar
J. Lee, J. H. Park, C. R. Lee, W. Y. Chung, and C. S. Park, “Long-term outcomes of total thyroidectomy versus thyroid lobectomy for papillary thyroid microcarcinoma: comparative analysis after propensity score matching,” Thyroid, vol. 23, no. 11, pp. 1408–1415, 2013.
View at: Publisher Site | Google Scholar
W. Zheng, J. Li, P. Lv, Z. Chen, and P. Fan, “Treatment efficacy between total thyroidectomy and lobectomy for patients with papillary thyroid microcarcinoma: a systemic review and meta-analysis,” European Journal of Surgical Oncology, vol. 44, no. 11, pp. 1679–1684, 2018.
View at: Publisher Site | Google Scholar
F. I. Macedo and V. K. Mittal, “Total thyroidectomy versus lobectomy as initial operation for small unilateral papillary thyroid carcinoma: a meta-analysis,” Surgical Oncology, vol. 24, no. 2, pp. 117–122, 2015.
View at: Publisher Site | Google Scholar
Y. Ito, A. Miyauchi, and H. Oda, “Low-risk papillary microcarcinoma of the thyroid: a review of active surveillance trials,” European Journal of Surgical Oncology, vol. 44, no. 3, pp. 307–315, 2018.
View at: Publisher Site | Google Scholar
I. Sugitani, Y. Ito, D. Takeuchi et al., “Indications and strategy for active surveillance of adult low-risk papillary thyroid microcarcinoma: consensus statements from the Japan association of endocrine surgery task force on management for papillary thyroid microcarcinoma,” Thyroid, vol. 31, no. 2, pp. 183–192, 2021.
View at: Publisher Site | Google Scholar
Z. Al-Qurayshi, N. Nilubol, R. P. Tufano, and E. Kandil, “Wolf in sheep's clothing: papillary thyroid microcarcinoma in the US,” Journal of the American College of Surgeons, vol. 230, no. 4, pp. 484–491, 2020.
View at: Publisher Site | Google Scholar
X. Zheng, S. Wei, Y. Han et al., “Papillary microcarcinoma of the thyroid: clinical characteristics and BRAF(V600E) mutational status of 977 cases,” Annals of Surgical Oncology, vol. 20, no. 7, pp. 2266–2273, 2013.
View at: Publisher Site | Google Scholar
L. Shi, L. Zhou, J. Wang et al., “The effect of the area proportion of the metastatic lesion within the central metastatic lymph node on response to therapy in papillary thyroid carcinoma,” Oncology Letters, vol. 21, no. 4, p. 284, 2021.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2024 Lei Jin et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

183

Downloads

218

Citations

International Journal of Endocrinology

Whether Detection of Gene Mutations Could Identify Low- or High-Risk Papillary Thyroid Microcarcinoma? Data from 393 Cases Using the Next-Generation Sequencing

Abstract

1. Introduction

2. Methods

2.1. Patients and Samples

2.2. Sample Preparation and Nucleic Acid Extraction

2.3. NGS Library Preparation

2.4. Sequencing and Bioinformatics Analysis

2.5. Statistical Analysis

3. Results

3.1. The Clinical and Pathological and Molecular Features of Tumor Size Less than 15 mm in 393 PTC Patients

3.2. Comparison of Gene Mutation Status and Clinical Pathological Characteristics between cN0 and cN1 Patients with Tumor Diameter Less than 15 mm

3.3. Correlation Analysis of Tumor Growth with Clinicopathological Features and Gene Mutation Status in cN0 PTC Patients with Tumor Size Less than 15 mm

3.4. Perform Univariate and Multivariate Logistic Regression Analyses on PTMC High-Risk Factors

3.5. The Association of the Presence of BRAFV600E Mutation with Clinicopathological Features of PTMC Patients

4. Discussion

4.1. Limitation

Data Availability

Disclosure

Conflicts of Interest

Authors’ Contributions

Acknowledgments

References

Copyright

3.5. The Association of the Presence of BRAF^V600E Mutation with Clinicopathological Features of PTMC Patients