Next Article in Journal
An Exploratory Analysis of Changes in Circulating Plasma Protein Profiles Following Image-Guided Ablation of Renal Tumours Provides Evidence for Effects on Multiple Biological Processes
Next Article in Special Issue
Advances in Biomarker-Driven Targeted Therapies in Thyroid Cancer
Previous Article in Journal
Is There a Benefit of Oxaliplatin in Combination with Neoadjuvant Chemoradiotherapy for Locally Advanced Rectal Cancer? An Updated Meta-Analysis
Previous Article in Special Issue
PPARγ Targets-Derived Diagnostic and Prognostic Index for Papillary Thyroid Cancer
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cancers
Volume 13
Issue 23
10.3390/cancers13236038
cancers-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The BRAFV600E Mutation Is Not a Risk Factor for More Aggressive Tumor Behavior in Radiogenic and Sporadic Papillary Thyroid Carcinoma at a Young Age

by
Liudmyla Zurnadzhy
1,2,†,
Tetiana Bogdanova
1,2,†,
Tatiana I. Rogounovitch
3,*,
Masahiro Ito
4,
Mykola Tronko
1,
Shunichi Yamashita
5,6,
Norisato Mitsutake
3,
Serhii Chernyshov
1,
Sergii Masiuk
7 and
Vladimir A. Saenko
2
1
State Institution “V.P. Komisarenko Institute of Endocrinology and Metabolism of the National Academy of Medical Sciences of Ukraine”, 69 Vyshgorodska Str., 04114 Kyiv, Ukraine
2
Department of Radiation Molecular Epidemiology, Atomic Bomb Disease Institute, Nagasaki University, 1-12-4 Sakamoto, Nagasaki 852-8523, Japan
3
Department of Radiation Medical Sciences, Atomic Bomb Disease Institute, Nagasaki University, Nagasaki 852-8523, Japan
4
Nagasaki Medical Center, 2-1001-1 Kubara, Omura 856-8562, Japan
5
Fukushima Medical University, Hikarigaoka 1, Fukushima 960-1295, Japan
6
National Institute of Radiological Sciences, National Institutes for Quantum Science and Technology, 4-9-1 Anagawa, Chiba 263-8555, Japan
7
State Institution “National Research Center for Radiation Medicine of the National Academy of Medical Sciences of Ukraine”, 53 Illienka Str., 04050 Kyiv, Ukraine
*
Author to whom correspondence should be addressed.
†
These two authors contribute equally to this work.
Cancers 2021, 13(23), 6038; https://doi.org/10.3390/cancers13236038
Submission received: 25 October 2021 / Revised: 26 November 2021 / Accepted: 29 November 2021 / Published: 30 November 2021
(This article belongs to the Special Issue Biomarkers of Thyroid Cancer)
Browse Figures
Versions Notes

Abstract

:

Simple Summary

Analysis of the groups of young Ukrainian patients (aged ≤28 years) with radiogenic and sporadic papillary thyroid carcinomas (PTCs) showed that the frequency of BRAFV600E was increasing with patient age, consistently remaining lower in radiogenic PTCs. In both etiopathogenic groups, the BRAFV600E-positive PTCs more frequently had a dominant papillary growth pattern, smaller tumor size, higher Ki67 labeling index, and a frequency of the major indicators of tumor invasiveness that is lower than or equal to that of the BRAFV600E-negative tumors. Comparison of the BRAFV600E-positive PTCs across the groups found a virtual absence of differences, while the BRAFV600E-negative tumors differed markedly and displayed a higher frequency of invasive tumor features in the radiogenic PTCs. Hence, there is evidence that BRAFV600E does not confer a more aggressive course of PTC in young patients regardless of tumor etiology.

Abstract

Histopathological changes in the fusion oncogene-driven papillary thyroid carcinomas (PTCs) from children and adolescents exposed to Chernobyl fallout have been extensively studied. However, characteristics of the radiogenic BRAFV600E-positive PTCs, whose proportion is growing with time, are not well described yet. We analyzed the relationship between the BRAFV600E status (determined immunohistochemically with the VE1 antibody) and the clinicopathological features of 247 radiogenic and 138 sporadic PTCs from young Ukrainian patients aged ≤28 years. The frequency of BRAFV600E was increasing with patient age, consistently remaining lower in radiogenic PTCs. In both etiopathogenic groups, the BRAFV600E-positive PTCs more frequently had a dominant papillary growth pattern, smaller tumor size, higher Ki67 labeling index, and a frequency of the major indicators of tumor invasiveness that is lower than or equal to that of the BRAFV600E-negative tumors. Comparison of the BRAFV600E-positive PTCs across the groups found a virtual absence of differences. In contrast, the BRAFV600E-negative radiogenic PTCs displayed less frequent dominant papillary and more frequent solid growth patterns, lower Ki67 labeling index, and higher invasiveness than the BRAFV600E-negative sporadic tumors. Thus, BRAFV600E is not associated with a more aggressive course of PTC in young patients regardless of etiology. The major clinicopathological differences between the radiogenic and sporadic PTCs are observed among the BRAFV600E-negative tumors.

1. Introduction

A sharp increase in the incidence of thyroid cancer among subjects aged ≤18 years in 1986 is the major health consequence of the Chernobyl accident in the exposed population. A number of studies addressing epidemiology, molecular characteristics, and clinical and pathological features of radiation-related thyroid cancer have been conducted to date.
Clinicopathological and molecular characteristics of the Chernobyl thyroid cancer, principally papillary thyroid carcinoma (PTC), have been shown to evolve with time. Tumors developing after the shorter latency (1st decade after the accident) were frequently more aggressive, especially in children and adolescents; many tumors had a solid growth pattern [1,2,3]. In contrast, PTCs diagnosed in the 2nd or 3rd decades after the accident after the longer period of latency displayed a less aggressive phenotype; the dominant papillary morphology was more common [3]. Furthermore, radiation-related PTCs in young patients were found to display pathological features of tumor aggressiveness more frequently than PTCs in non-exposed patients matched for the place of residence and age group [4].
Time-dependent changes in the clinicopathological characteristics of pediatric Chernobyl thyroid cancer were paralleled by changes in the genetic events underlying PTC. The early-onset tumors were frequently driven by fusion oncogenes such as various types of RET/PTC, more often RET/PTC3 [5,6,7] followed by a shift to RET/PTC1 after the first decade post-accident [8]. PTCs diagnosed later had a different spectrum of gene rearrangements, among which the rearranged BRAF, ETV6-NTRK3 and other gene fusions were discovered [9,10,11].
Studies of point mutations in pediatric Chernobyl PTC demonstrated the absence of genetic alterations in the RAS family genes [10,12,13]. The mutant BRAF was not found in the early-onset Chernobyl PTC, then its prevalence reached approximately 10% in the mid-1990s [14] and grew to about 15% after the 2000s [9,10,15]. One of these studies claimed the frequency of mutant BRAF was related to patient age but not to a history of radiation exposure [15]. Of interest, the frequencies of point mutations in radiation-related PTC seems to be declining with increasing 131I thyroid dose, and, in contrast, an uptrend is observed for fusion genes [9,11,16]. In sporadic pediatric PTC, BRAFV600E was not detected in a small group of non-exposed children from Ukraine [15]. However, more recent works reported BRAFV600E in about one-fourth of non-exposed children and adolescents with PTC [10,17]. Each of these studies had certain strengths and limitations; the latter mostly due to the relatively small sample sizes and frequent unavailability of the appropriate control groups.
Thus, on the one hand, radiogenic PTC displays time-dependent changes in morphology and aggressiveness, and on the other, there are changes in its molecular landscape. The genotype–phenotype relationships are better understood for the tumors harboring fusion oncogenes while less is known about clinicopathological correlations for the BRAFV600E-driven PTCs in young patients.
The aim of this study was to analyze: (i) the frequency of the BRAFV600E mutation in radiogenic and sporadic PTCs from the young Ukrainian patients aged 4–28 years at diagnosis in the whole groups and age subgroups of children (aged ≤14), adolescents (aged 15–≤18) and young adults (aged ≥19 to 28); (ii) to examine the relationship of the BRAFV600E to 131I thyroid dose and the duration of the period of latency; and (iii) to compare the structural and invasive characteristics, proliferative activity and the frequency of postoperative recurrent lymph node metastases in the BRAFV600E-positive and BRAFV600E-negative PTCs within and across radiogenic and sporadic PTC series.

2. Results

2.1. Frequency of BRAFV600E in Radiogenic and Sporadic PTC

Positive IHC reactions with the antibody to the mutant BRAFV600E protein were detected in the radiogenic series in only one of 104 PTC cases (1.0%) in children, in 6/52 of cases (11.5%) in adolescents, and 19/91 of cases (20.9%) in young adults (Table 1). In contrast, BRAFV600E was detected in five of 39 cases (12.8%) in children, 12/37 of cases (32.4%) in adolescents, and 25/62 cases (40.3%) in adults. In each age group, the frequency of BRAFV600E in radiogenic PTC was significantly lower than in sporadic PTC. Both etiopathogenic series displayed age-related uptrends for the BRAFV600E frequency (ptrend = 5.97 × 10−6 and ptrend = 0.004 for radiogenic and sporadic series, respectively). Logistic regression adjusted for sex also confirmed the increasing with age chance of BRAFV600E-positivity in radiogenic (OR = 1.179, p = 2.40 × 10−5) and sporadic (1.118, p = 0.007) series. Thus, the frequency of the BRAFV600E mutation in PTC is age-dependent regardless of a radiation history.

2.2. BRAFV600E Associations in Radiogenic PTC

The results of univariate and multivariate statistical analyses are presented in Table 2. Patients with the BRAFV600E-positive PTC were significantly older at the time of operation than those with the BRAFV600E-negative tumors (b = 6.012, p = 3.00 × 10−6), and the BRAFV600E-positive PTCs were characterized by the longer period of latency (b = 6.133, p = 2.00 × 10−6). The latter observation is unique for the radiation-related cancers since knowledge of time points of exposure to radiation and of diagnosis enables accurate determination of latency. We found that, while a conventional proportional hazard model did not perform well (Figure S1a–c), an extended Cox model in which the BRAF status was introduced as a time-dependent variable provided reasonable fit (Figure S1d). Model parameters adequately reflected the delaying development of the BRAFV600E-positive PTCs (HR = 0.021, p = 3.67 × 10−5) and an increasing rate of such tumors in time (i.e., with increasing latency, HR = 1.194, p = 8.77 × 10−5 for the BRAF × Latency interaction term).
The 131I thyroid radiation doses in patients with BRAFV600E-positive PTCs were significantly lower than in patients with BRAFV600E-negative PTCs on univariate analysis (median 200 mGy vs. 350 mGy, respectively, p = 0.019) and in a multivariate model (b = −0.295, p = 0.031).
BRAFV600E was statistically significantly associated with several histopathological characteristics, such as smaller tumor size ((b = −0.528, p = 0.013), paralleled by its higher frequency in microcarcinomas and in pT1a tumors (OR = 3.475, p = 0.008 for both), more frequent dominant papillary (OR = 5.735, p = 1.93 × 10−4; Figure 1a,b) and less frequent solid-trabecular dominant growth pattern than in the BRAFV600E-negative PTCs (OR = 0.289, p = 0.032; Figure 1d,e).
On multivariate analysis, BRAFV600E did not confer a statistically significant association with any feature of tumor aggressiveness, including multifocality, lymphatic/vascular invasion, extrathyroidal extension, regional or distant metastasis (p = 0.064 for the strongest signal). For many instances, we again observed the influence of patient age at the time of operation, a powerful confounder, which rendered the effect of BRAFV600E non-significant. If only sex was accounted for, a number of associations would be significant such as the lower frequencies of lymphatic/vascular invasion, extrathyroidal extension and N1b. Altogether, these observations suggest that BRAFV600E is unlikely to be associated with higher tumor aggressiveness in the analyzed series of young PTC patients.
Interestingly, despite the smaller tumor size, the BRAFV600E positivity was associated with the higher Ki67 LI (b = 1.674, p = 0.006), mostly due to a higher proportion of tumors with Ki67 LI ranging from 5% to 10% (OR = 3.974, p = 0.004) (Figure 1c,f)). The lower Ki67 LI (0–5%) was characterized by the significantly less frequent BRAFV600E (OR = 0.268, p = 0.004). Ki67 LI was associated with smaller tumor size in the whole radiogenic group (b = −0.672, p = 0.030), but no statistical significance could be reached in the BRAFV600E-positive (b = −0.125, p = 0.848) and the BRAFV600E-negative (b = −0.523, p = 0.135) PTCs separately.
Repeated operations for recurrent PTC in regional lymph nodes performed at least 6 months after the primary surgery occurred only in the BRAFV600E-negative group in 8/210 (3.8%, p = 0.574) cases with available follow-up information. We confirmed that reoperated recurrent lymph node metastases were BRAFV600E-negative in all eight cases as with the primary tumors (Cohen’s κ = 1.000, p = 0.001). The absence of recurrences in patients with BRAFV600E-positive PTCs, however, did not result in statistically significant differences to the BRAFV600E-negative subgroup in a multivariate model (HR = 0.523, p = 0.708). Perhaps the small total number of recurrent events did not provide sufficient statistical power to detect the difference if any existed.

2.3. BRAFV600E Associations in Sporadic PTC

Associations of the BRAFV600E-positivity with clinicopathological features in sporadic PTCs generally paralleled those in the radiogenic series (Table 3). Patients with the BRAFV600E-positive PTCs were significantly older at the time of operation (b = 2.516, p = 0.006). The BRAFV600E-positive tumors were smaller in size (b = −8.706, p = 1.52 × 10−4), and among those, the microcarcinomas and pT1a PTCs were more frequent (OR = 3.499, p = 0.003 for both). The BRAFV600E-positive PTCs more frequently displayed dominant papillary (OR = 5.462, p = 7.60 × 10−5) and, less frequently, follicular growth patterns (OR = 0.143, p = 0.001), and had more frequent oncocytic changes (OR = 2.488, p = 0.020).
Except for the lower frequency in tumors with lymphovascular invasion (OR = 0.304, p = 0.004), BRAFV600E did not associate with the major indices of tumor invasiveness including multifocality, extrathyroidal extension, and regional and distant metastases (multivariate p = 0.068, the strongest). The integrative invasiveness score or any of its categories did not differ between the BRAFV600E-positive and BRAFV600E-negative PTCs (p = 0.217, the strongest).
In concordance with observations in radiogenic PTCs, sporadic BRAFV600E-positive tumors were associated with significantly higher proliferative activity (b = 2.221, p = 0.002), due to the tumors with Ki67 LI from 5 to 10% (OR = 3.109, p = 0.006). In the tumors with low Ki67 LI, the frequency of BRAFV600E was significantly lower (OR = 0.231, p = 5.49 × 10−4) indicative of a preferential association of lower Ki67 LI with the BRAFV600E-negative status. Again, Ki67 LI was associated with smaller tumor size in the whole sporadic PTC group (b = −0.785, p = 0.005), but no statistical significance was found in the BRAFV600E-positive (b = −0.711, p = 0.056) and the BRAFV600E-negative (b = −0.517, p = 0.158) PTCs in separate analyses.
There were four recurrences in this group that were detected and reoperated six or more months after the initial surgery. Three primary tumors and reoperated metastases were BRAFV600E-positive, and the other and its recurring metastasis were BRAFV600E-negative. There was a perfect agreement in the BRAF status between the primary and recurring tumors (Cohen’s κ = 1.000, p = 0.001). The three BRAF-positive recurrent metastases were refractory to radioiodine treatment (not shown in the Table 3). Despite the relatively small number of patients with follow-up information in this group, the BRAFV600E-positivity was associated with a higher chance of recurrence on both univariate (p = 0.007 by the log-rank test) and multivariate analyzes (HR = 19.042, p = 0.031).

2.4. Comparison of the BRAFV600E-Positive or the BRAFV600E-Negative PTCs across the Radiogenic and Sporadic Series

Finally, we addressed the role of tumor etiology in PTCs with the same BRAF status. The BRAFV600E-positive PTCs from radiogenic and sporadic series were largely similar in their characteristics (Table 4). Only three parameters were statistically different: the older age at operation (b = 2.838, p = 0.018), lower frequency of oncocytic changes (OR = 0.242, p = 0.012) and longer follow-up period (b = 13.324, p = 8.82 × 10−7) in the BRAFV600E-positive radiogenic PTCs. Of note, there were no recurrences in the BRAFV600E-positive radiogenic PTCs (0/26), while 3/42 tumors recurred in the BRAFV600E-positive sporadic PTC subgroup. This difference was statistically significant on univariate analysis (p = 0.027 by the log-rank test) but not in a multivariate model (HR = 0.031, p = 0.147). No evidence of differences for other clinicopathological characteristics was found.
In contrast, the BRAFV600E-negative PTCs displayed a number of differences (Table 4). The radiogenic BRAFV600E-negative PTCs were more likely to be diagnosed in male patients (OR = 2.952, p = 0.001), more frequently displayed a solid-trabecular dominant growth pattern (OR = 2.350, p = 0.002), gross extrathyroidal extension (OR = 5.219, p = 0.008) and corresponding pT3b category (OR = 5.599, p = 0.006), higher invasiveness score (OR = 1.644, p = 0.027), and were followed-up for a longer time (b = 6.413, p = 3.81 × 10−18). Several characteristics were less frequent in the radiogenic BRAFV600E-negative PTCs: tumor encapsulation (OR = 0.315, p = 2.15 × 10−4), dominant papillary growth pattern (OR = 0.416, p = 0.002), oncocytic changes (OR = 0.298, p = 4.00 × 10−5), pT2 category (0.385, p = 6.49 × 10−4), zero invasiveness score (OR = 0.549, p = 0.046); the tumors also had an overall lower Ki67 LI (b = −1.714, p = 8.00 × 10−6), in particular 5–10% LI (OR = 0.488, p = 0.022).
The obtained results indicated that the major difference in clinicopathological characteristics between the radiogenic and sporadic PTCs were seen between the BRAFV600E-negative PTCs. Since the absence of statistically significant differences between the BRAFV600E-positive PTCs might be due in part to the relatively small sample size, thus insufficient statistical power, we performed a correspondence analysis to visualize the similarities or differences between the etiological groups of PTCs with different BRAF status. This type of statistical analysis employs mathematical and computation apparatuses different from those in the regression analysis and conceptually relates to principal component analysis but intended to explore and graphically present relationships between categorical (qualitative) variables in a low-dimension plot. The results of correspondence analysis (Figure 2) clearly demonstrated an obvious similarity between the BRAFV600E-positive radiogenic and sporadic PTCs (note an acute angle between these groups). These tumors were concordantly associated with older patient age, frequent micro-PTC, dominant papillary growth pattern, Ki67 LI from 5% to 10%, invasiveness score 0 or 1, and with the absence of regional (N0) and distant (M0) metastases, and of lymphovascular invasion.
The BRAFV600E-negative PTCs displayed marked differences in their characteristics. The radiogenic BRAFV600E-negative PTCs were the most distinct group of tumors associated with younger patient age, solid-trabecular growth pattern, absence of oncocytic changes, higher invasiveness scores, the presence of distant metastases and low Ki67 LI. The sporadic BRAFV600E-negative PTCs did not display obvious association with a particular growth pattern and appeared to be less invasive. These observations were well in line with and corroborate the results of our regression analyses presented above in Table 2, Table 3 and Table 4.
Cumulatively, our analyses suggest a clinicopathological resemblance of the BRAFV600E-positive radiogenic and sporadic PTCs and indicate that such tumors are not associated with the higher tumor aggressiveness in young patients.

3. Discussion

Here, we analyzed the clinicopathological characteristics of PTCs in the groups of young patients exposed or not to Chernobyl radiation in Ukraine with regard to the BRAF mutational status. The availability of two groups of patients with PTCs of different etiologies made the analyses within and across the groups possible.
First, our study confirmed the notion that the frequency of the BRAFV600E mutation in PTC increases with patient age (see Table 1). Literature data on age correlations of BRAFV600E are sometimes controversial, either confirming [15,16,18,19,20,21,22,23,24,25,26,27,28] (meta-analysis [29]) this or not [30,31,32,33] (meta-analyses [34,35]). The controversy may likely be explained by the rarity of PTC in early childhood and thus insufficient statistical power in some studies to detect age relatedness due to the underrepresentation of very young patients (i.e., lower age variability) in the groups available for analysis. Analyses of the BRAFV600E association with age in young patients from Fukushima, whose age distribution is relatively close to that in the present study, clearly indicated the increasing BRAF mutational frequency with patient age [24,28]. Our data also showed that the frequency of BRAFV600E in PTC rapidly increases with patients’ age transition from childhood to adolescence, and that this change is observed in both exposed and non-exposed patients. Of importance, however, the frequency of BRAFV600E remains consistently higher in sporadic PTCs as compared to radiogenic PTCs in all age subgroups. The latter observation is in line with a recent report on molecular findings in childhood PTC from Belarus [17].
To the best of our knowledge, the group of radiation-exposed children aged ≤14 years at operation for PTC presented in this study is the largest in the world (104 patients) who have been examined for the presence of BRAFV600E in the tumor. Earlier works were concordant in detecting no BRAFV600E in PTCs from exposed children of this age [10,14,15]. Our study demonstrates that BRAFV600E may exist in PTCs from such patients, albeit exceptionally rare. The BRAFV600E-positive radiogenic PTCs have a significantly longer period of latency (the time between exposure and tumor detection, see Table 2) resulting in the older age of patients at diagnosis. Naturally, because of the longer latency, the vast majority of children already migrate to older age groups by the time a BRAFV600E-positive tumor becomes detectable. In our study, the shortest period of latency of a BRAFV600E-positive PTC was 11.3 years; the tumor was diagnosed in a girl exposed to radiation at the age of 2.2 years. This example illustrates that only a combination of two rare factors, a very young age at exposure and a relatively short period of latency (as compared to the median of 22.1 years) for a BRAFV600E-positive PTC may create a condition when the patient’s age at diagnosis remains below 14 years old. Therefore, most BRAFV600E-positive PTCs occur in adolescents and adults in the radiogenic group. Our analysis of the BRAFV600E-positive and BRAFV600E-negative PTC onset in time after exposure also attests to the shift of the BRAFV600E-positive tumor development toward the longer periods of latency (see Figure S1).
In the radiogenic group, we also confirmed that the BRAFV600E mutation was associated with a lower radiation dose to the thyroid reported earlier in the Ukrainian Chernobyl PTCs [9,11,16]. Note that our study is independent of the previous works with minimal, if any, overlap in included cases.
Our analysis of clinicopathological correlations of the BRAFV600E mutation demonstrated a good concordance in the association pattern between the two etiological groups (see Table 2 and Table 3). There were a few clinicopathological characteristics’ associations which were statistically significant in one or another group, but we do not interpret those as etiology-specific since OR estimates corresponded to each other (i.e., ORS were concordantly greater or less than the value of 1). In contrast, the difference in hazard ratios for the recurrence between sporadic and radiogenic groups (HR = 19.042, p = 0.031 and HR = 0.523, p = 0.708, respectively) is suggestive that the BRAFV600E mutation may be a risk factor for recurrence at least in sporadic PTC even in young patients, which corresponds to the results of meta-analyses in adult patients [29,36,37,38,39,40].
The association of BRAFV600E with the higher Ki67 LI is of particular interest. Such a correlation was reported in adult PTC patients [18,41] but no data are available for young patients. In general, the higher Ki67 LI was associated with the greater size of PTCs in adult patients [18,42,43,44] and poorer prognosis [30,45,46,47]. The BRAFV600E mutation was also associated with the greater tumor size in some studies (meta-analyses [34,35]), although equivocally [24,27,28] (meta-analysis [39]), and with higher tumor aggressiveness (meta-analyses [29,34,35,36,38,39,48]). Our findings, therefore, provide somewhat counterintuitive yet robust evidence that, in young patients, the BRAFV600E-positive tumors, which are characterized by the higher Ki67 LI and smaller size, do not display a more aggressive phenotype. These results are in line with several previous reports of the lack of overly aggressive features in the BRAFV600E-positive PTC in pediatric patients [17,24,49,50,51,52,53].
It is tempting to speculate that, biologically, the BRAFV600E-positive PTCs with higher Ki67 LI, smaller size and diagnosed in older patients may comprise the tumors that were developing after the longer “silent” period (thus, older patient age and/or longer period of latency in exposed patients) after entering into the phase of more active growth (thus, elevated KI67 LI) when they were detected, but at this stage having not yet acquired the aggressive features observed in the BRAFV600E-positive PTCs in adult/elderly patients. In this regard, it would be interesting to investigate age-related changes of tumor aggressiveness and long-term outcomes in the BRAFV600E-positive PTCs in a special group of patients aged from children to old age. The absence of correlations between BRAFV600E and the aggressive features of PTCs in young patients does not rule out that at least some of these tumors might develop a more aggressive phenotype with time (and thus in more aged patients). Alternatively, some tumors in young patients may not be progressive, as was seen in Fukushima patients [54] in whom the frequency of the BRAFV600E mutation is about 70% [24,28]. Perhaps some BRAFV600E-positive PTCs from our series with lower Ki67 LI (about 40% of all BRAFV600E-positive tumors had KI67 LI ≤5%, see Table 2 and Table 3) might be those that would not progress.
Finally, we compared the BRAFV600E-positive and BRAFV600E-negative PTCs across the radiogenic and sporadic groups. Since a pronounced similarity was observed between the clinicopathological association patterns of the BRAFV600E-positive tumors between the radiogenic and sporadic PTC groups mentioned above, it was somewhat expected that the BRAFV600E-positive PTCs did not display many differences upon their comparisons in the two etiological groups (see Table 4).
In contrast, the BRAFV600E-negative PTCs displayed a number of differences, both in tumor morphology and the frequencies of invasive features, pointing at the higher aggressiveness of radiogenic tumors, in line with our previous report [4]. The BRAFV600E-negative PTCs are most likely driven by fusion oncogenes, whose histopathological associations have been established in earlier works [6,7,9,11,15,55]. Of importance, however, there are differences in the distribution of fusion oncogene types in the tumors of different etiological groups. The radiogenic PTCs, especially those detected after the shorter latency, had a very high frequency of RET/PTC3 rearrangements, which is associated with more aggressive tumor phenotype and a solid-trabecular growth pattern [5,6,7,8,56,57]. The frequency of RET/PTC3 is lower in sporadic PTC, even from young children, and other types of activated oncogenes do not seem to confer very high tumor aggressiveness in young patients [17,58].
Our work has several strengths and limitations. The strengths include a relatively large number of young patients in the study, the availability of two etiological groups of PTCs enabling the analyses within and across the groups, the high quality of demographic and clinicopathological information, knowledge of individual radiation doses to the thyroid, and data on the BRAF status, which was the focus of this investigation. The major limitation is the lack of molecular analyses for other oncogenic drivers, knowledge of which could facilitate more sophisticated and more detailed assessment of clinicopathological relationships in radiogenic and sporadic PTCs in young patients. Overcoming this shortcoming, however, is connected to several technical and practical issues, including the unavailability of tissues suitable for molecular analyses in some cases or the insufficient quality of nucleic acids from formalin-fixed paraffin-embedded tissues after a long-term storage.

4. Materials and Methods

4.1. Patients

The radiogenic cases included PTCs from 247 patients aged 6–28 years at diagnosis who were operated at the State Institution “V.P. Komisarenko Institute of Endocrinology and Metabolism of the National Academy of Medical Sciences of Ukraine” (IEM), Kyiv during the period from 1990 to 2014, that is, since the significant increase in thyroid cancer incidence after the Chernobyl accident [4,59,60]. Given that the highest thyroid cancer risk was observed in the youngest children who lived in the most contaminated by 131I regions of northern Ukraine [60,61], we defined radiogenic cases as those diagnosed in children aged ≤4 years in April, 1986, who lived in the Kyiv, Chernihiv and Zhytomyr regions. Additionally, the patients were divided into children (≤14), adolescents (15–18) and young adults (19–28 years old at diagnosis) for subgroup analysis.
The group of comparison, sporadic PTC cases, were from 138 patients aged 4–28 years at diagnosis born after the Chernobyl accident (from 1 January 1987 or later, i.e., not affected by 131I), who corresponded by age at operation and regions of residence (including Kyiv-city) to those of patients with radiogenic PTC. Patients were operated for PTC from 1997 to 2015. Since the oldest patient with sporadic PTC at the time of this study was aged 28 years, the upper limit for the age of patients with radiogenic PTC was also set to 28 years.
All cases in this work are a subset of our earlier histopathological study that included 301 radiogenic and 194 sporadic PTCs from patients aged up to 28 years [4]. Only those cases with sufficient tumor tissues remaining in the paraffin blocks were included in the current investigation to enable immunohistochemical staining.

4.2. Histopathology

The histopathological examination of hematoxylin/eosin-stained paraffin sections was performed by two experienced pathologists of IEM (TB and LZ). The pathological diagnosis was based on the 4th edition of the WHO histological classification [62]. Most cases were also reviewed by the international pathology panel of the CTB project [63,64]. The diagnosis of PTC was confirmed in all analyzed cases. TNM categories were determined according to the 8th edition of the pTNM classification [65]. Tumors were classified according to the dominant histological growth pattern into three categories: papillary, follicular or solid-trabecular, when the corresponding structural component exceeded 50% of a tumor section surface, and was also evaluated for the presence of oncocytic (oxyphilic/Hurtle) cell changes in the tumor epithelium.
As in our previous work [60], in addition to conventional clinicopathological features, we used an integrative variable, the “invasiveness score”, which is the arithmetic sum of every instance of multifocality, lymphatic/vascular invasion, any extrathyroidal extension (i.e., minimal or gross), N1 and M1 (commonly detected by diagnostic imaging), either isolated or in combination with other(s), for each tumor. Thus defined, the invasiveness score ranged from 0 (no invasive feature presents) to 5 (all features present).

4.3. Immunohistochemistry

Immunohistochemical (IHC) staining for BRAFV600E was performed according to the protocol used in the Department of Radiation Molecular Epidemiology of the Atomic Bomb Diseases Institute, Nagasaki University (LZ, TIR) using mouse monoclonal anti-BRAF (mutated V600E) antibody (VE1) ab228461 (Abcam, Tokyo, Japan) at a 1:100 dilution applied for one hour at 37 °C. The Novolink Polymer Detection System (250T) (Leica RE7140-K) was used to detect the IHC reaction product, which included treatment with the secondary rabbit anti-mouse antibody (IgG), attachment of the peroxidase label and visualization with DAB diluted in the Novolink DAB Substrate Buffer according to the manufacturer’s recommendations. Cell nuclei were stained with Mayer’s hematoxylin. The IHC reaction was considered positive (expression of the BRAFV600E mutant protein) in the presence of the brown color of the cytoplasm of tumor epithelial cells. Sections of a formalin-fixed paraffin-embedded tumor tissue from a patient not related to this study with PTC with the BRAFV600E mutation confirmed by Sanger sequencing were used as a positive control.
We considered a positive IHC reaction for BRAFV600E as indicating the presence of the BRAFV600E mutation at the DNA level. A high concordance between IHC methods employing the VE1 antibody and PCR- or direct sequencing-based techniques was confirmed in a recent meta-analysis [66], and also by our group using formalin-fixed paraffin-embedded material [67].
The proliferative activity of tumors was evaluated by IHC using Ki67 antibody (clone MIB-1; DAKO, Glostrup, Denmark, 1:100 dilution) in a Ventana BenchMark ULTRA instrument. The Ki67 LI was determined with the image-analyzing software (CountσCell, Ki67 antigen Semi Auto Counter, Seiko Tec LTD, Fukuoka, Japan) by counting Ki67-positive and -negative PTC nuclei in the tumor areas with the highest number of immunoreactive nuclei, and calculating the proportion of Ki67-positive nuclei. For each case, a total of approximately 1000 PTC cells were analyzed (LZ). Image analysis was performed in a blind for the BRAFV600E status manner.

4.4. Thyroid Dosimetry

131I thyroid radiation doses were calculated for each patient from the radiogenic series in the Dosimetry department of the State Institution “National Research Center for Radiation Medicine of the National Academy of Medical Sciences of Ukraine”, Kyiv, using an ecological dosimetric model, which includes the system of ecological iodine transport and biokinetic models of iodine (“TD-CTB”) [62].

4.5. Statistical Analysis

The Fisher’s exact test, the Fisher–Freeman–Halton exact test and the Cochran–Armitage test were used for univariate analysis of categorical data; the Mann–Whitney test was used to compare continuous data between any two groups. Logistic regression models were adjusted for age at operation and sex; models with very small numbers of outcomes (<5 per cell) were conducted using Firth’s approach to bias-reducing penalized maximum likelihood fit. Multivariable linear regression models were used for continuous dependent variables. The occurrence of the BRAFV600E mutation in relation to the period of latency was assessed using survival analysis methods. The Kaplan–Meier method, and the proportional hazard (Cox) and extended proportional hazard models were used. Computation and plotting of the results of the model with time-varying coefficients was performed with an SAS macro “coxtvc” [68]. Calculations were performed using IBM SPSS Statistics Version 24 software (International Business Machines Corp., Armonk, NY, USA) or the 9.4 version of SAS (SAS Institute, Cary, NC, USA). Correspondence analysis was performed in R with the “ca” package [69]. The “colgreen” option was used to calculate biplot principal coordinates for the four PTC groups and contribution coordinates (the standard coordinates multiplied by the square root of the corresponding masses) for categorical clinicopathological features. All tests were two-sided; p < 0.05 was considered statistically significant.

5. Conclusions

Our study demonstrates that the BRAFV600E mutation is more frequent in sporadic than in radiogenic PTCs in all age groups of patients, and that mutational frequency significantly increases with patient age in both radiogenic and sporadic PTCs. The thyroid radiation dose in patients with the BRAFV600E-positive PTC was significantly lower, and the period of latency was significantly longer compared to that in patients with BRAFV600E-negative PTC. The major histopathological differences between the radiogenic and sporadic PTCs were observed among the BRAFV600E-negative tumors; the radiogenic PTCs displayed morphological features of tumor aggressiveness more frequently than the sporadic ones; the latter had a somewhat milder clinical phenotype. The BRAFV600E-positive PTCs displayed similar clinicopathological association patterns in the radiogenic and sporadic series, including a high frequency of microcarcinomas, dominant papillary growth pattern, high Ki67 LI, and the presence of oncocytic changes in tumor epithelial cells. The BRAFV600E-positive PTCs were characterized by invasive properties that were lower than or comparable to those in the BRAFV600E-negative tumors in both radiogenic and sporadic series, indicating that the BRAFV600E mutation is not associated with more aggressive tumor behavior in patients of young age regardless of PTC etiology. Further studies, ideally addressing all driver oncogenes and other cancer genes at the genomic and transcriptomic levels, would be desired to determine whether the phenotype and prognosis of the BRAFV600E-positive radiogenic PTCs will be acquiring patient age-related changes similarly to those described in sporadic PTCs.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/cancers13236038/s1, Figure S1: Checking the proportional hazard assumption for the BRAF status variable in a proportional hazard and in an extended Cox model of PTC onset in time after exposure (Latency, years) for the radiogenic PTCs.

Author Contributions

Conceptualization, L.Z., T.B., T.I.R., M.T., S.Y., N.M. and V.A.S.; methodology, L.Z., T.I.R., T.B. and M.I.; clinical data L.Z., T.B., S.C.; pathological data L.Z., T.B.; formal analysis, L.Z., T.B., T.I.R. and M.I.; investigation, L.Z., T.B. and T.I.R.; dosimetry, S.M.; data curation, L.Z., T.B., T.I.R. and V.A.S.; statistical analysis T.B. and V.A.S.; writing—original draft preparation, T.B. and V.A.S.; writing—review and editing, L.Z., T.B., T.I.R., M.I., M.T., S.Y., N.M., S.C., S.M. and V.A.S.; supervision, T.B., M.T., N.M. and V.A.S.; funding acquisition, L.Z., T.I.R., M.I. and V.A.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported in part by the Program of the Network-Type Joint Usage/Research Center for Radiation Disaster Medical Science, intramurally by the Atomic Bomb Disease Institute, Nagasaki University, and the Japan Society for the Promotion of Science (JSPS), KAKENHI Grant Numbers 19K07471 and 19KK02670001.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the IEM Bioethics Committee (protocols N 22-KE of 26 April 2018 and N 31-KE of 27 February 2020), the Chernobyl Tissue Bank (CTB, project N001-2020) and the Ethics Committee of Nagasaki University (protocol 20130401–7 of 1 July 2021, the latest update).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study or their guardians (for minors).

Data Availability Statement

All data generated or analyzed during this study are included in this published article.

Acknowledgments

The authors acknowledge the commitment of the staff of the Laboratory of Morphology of Endocrine System and of the Department of Surgery of Endocrine System of IEM, who prepared all pathological material and operated the patients, respectively. The authors gratefully acknowledge the confirmation of Ukrainian diagnoses provided by the International Pathology Panel of the Chernobyl Tissue Bank, which is supported by NCI grant number U24CA082102: A. Abrosimov, T. Bogdanova, G. Fadda, J. Hunt, M. Ito, V. Livolsi, J. Rosai, E.D. Williams, N. Dvinskyh and L. Zurnadzhy.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tronko, M.D.; Bogdanova, T.I.; Komissarenko, I.V.; Epstein, O.V.; Oliynyk, V.; Kovalenko, A.; Likhtarev, I.A.; Kairo, I.; Peters, S.B.; LiVolsi, V.A. Thyroid carcinoma in children and adolescents in Ukraine after the Chernobyl nuclear accident: Statistical data and clinicomorphologic characteristics. Cancer 1999, 86, 149–156. [Google Scholar] [CrossRef]
  2. Williams, E.D.; Abrosimov, A.; Bogdanova, T.; Demidchik, E.P.; Ito, M.; LiVolsi, V.; Lushnikov, E.; Rosai, J.; Sidorov, Y.; Tronko, M.D.; et al. Thyroid carcinoma after Chernobyl latent period, morphology and aggressiveness. Br. J. Cancer 2004, 90, 2219–2224. [Google Scholar] [CrossRef] [Green Version]
  3. Bogdanova, T.; Zurnadzhy, L.; LiVolsi, V.A.; Williams, E.D.; Ito, M.; Nakashima, M.; Thomas, G.A. Thyroid cancer pathology in Ukraine after Chernobyl. In Thyroid Cancer in Ukraine after Chernobyl: Dosimetry, Epidemiology, Pathology, Molecular Biology; Tronko, M., Bogdanova, T., Saenko, V., Thomas, G.A., Likhtarov, I., Yamashita, S., Eds.; IN-TEX: Nagasaki, Japan, 2014; pp. 65–108. [Google Scholar]
  4. Bogdanova, T.I.; Saenko, V.A.; Brenner, A.V.; Zurnadzhy, L.Y.; Rogounovitch, T.I.; Likhtarov, I.A.; Masiuk, S.V.; Kovgan, L.M.; Shpak, V.M.; Thomas, G.A.; et al. Comparative Histopathologic Analysis of “Radiogenic” and “Sporadic” Papillary Thyroid Carcinoma: Patients Born Before and After the Chernobyl Accident. Thyroid 2018, 28, 880–890. [Google Scholar] [CrossRef]
  5. Nikiforov, Y.E.; Rowland, J.M.; Bove, K.E.; Monforte-Munoz, H.; Fagin, J.A. Distinct pattern of ret oncogene rearrangements in morphological variants of radiation-induced and sporadic thyroid papillary carcinomas in children. Cancer Res. 1997, 57, 1690–1694. [Google Scholar]
  6. Thomas, G.A.; Bunnell, H.; Cook, H.A.; Williams, E.D.; Nerovnya, A.; Cherstvoy, E.D.; Tronko, N.D.; Bogdanova, T.I.; Chiappetta, G.; Viglietto, G.; et al. High prevalence of RET/PTC rearrangements in Ukrainian and Belarussian post-Chernobyl thyroid papillary carcinomas: A strong correlation between RET/PTC3 and the solid-follicular variant. J. Clin. Endocrinol. Metab. 1999, 84, 4232–4238. [Google Scholar] [CrossRef] [PubMed]
  7. Santoro, M.; Thomas, G.A.; Vecchio, G.; Williams, G.H.; Fusco, A.; Chiappetta, G.; Pozcharskaya, V.; Bogdanova, T.I.; Demidchik, E.P.; Cherstvoy, E.D.; et al. Gene rearrangement and Chernobyl related thyroid cancers. Br. J. Cancer 2000, 82, 315–322. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Rabes, H.M.; Demidchik, E.P.; Sidorow, J.D.; Lengfelder, E.; Beimfohr, C.; Hoelzel, D.; Klugbauer, S. Pattern of radiation-induced RET and NTRK1 rearrangements in 191 post-chernobyl papillary thyroid carcinomas: Biological, phenotypic, and clinical implications. Clin. Cancer Res. 2000, 6, 1093–1103. [Google Scholar] [PubMed]
  9. Leeman-Neill, R.J.; Brenner, A.V.; Little, M.P.; Bogdanova, T.I.; Hatch, M.; Zurnadzy, L.Y.; Mabuchi, K.; Tronko, M.D.; Nikiforov, Y.E. RET/PTC and PAX8/PPARgamma chromosomal rearrangements in post-Chernobyl thyroid cancer and their association with iodine-131 radiation dose and other characteristics. Cancer 2013, 119, 1792–1799. [Google Scholar] [CrossRef] [Green Version]
  10. Ricarte-Filho, J.C.; Li, S.; Garcia-Rendueles, M.E.; Montero-Conde, C.; Voza, F.; Knauf, J.A.; Heguy, A.; Viale, A.; Bogdanova, T.; Thomas, G.A.; et al. Identification of kinase fusion oncogenes in post-Chernobyl radiation-induced thyroid cancers. J. Clin. Investig. 2013, 123, 4935–4944. [Google Scholar] [CrossRef] [Green Version]
  11. Efanov, A.A.; Brenner, A.V.; Bogdanova, T.I.; Kelly, L.M.; Liu, P.; Little, M.P.; Wald, A.I.; Hatch, M.; Zurnadzy, L.Y.; Nikiforova, M.N.; et al. Investigation of the Relationship Between Radiation Dose and Gene Mutations and Fusions in Post-Chernobyl Thyroid Cancer. JNCI J. Natl. Cancer. Inst. 2017, 110, 371–378. [Google Scholar] [CrossRef]
  12. Nikiforov, Y.E.; Nikiforova, M.N.; Gnepp, D.R.; Fagin, J.A. Prevalence of mutations of ras and p53 in benign and malignant thyroid tumors from children exposed to radiation after the Chernobyl nuclear accident. Oncogene 1996, 13, 687–693. [Google Scholar] [PubMed]
  13. Kumagai, A.; Namba, H.; Saenko, V.A.; Ashizawa, K.; Ohtsuru, A.; Ito, M.; Ishikawa, N.; Sugino, K.; Ito, K.; Jeremiah, S.; et al. Low frequency of BRAFT1796A mutations in childhood thyroid carcinomas. J. Clin. Endocrinol. Metab. 2004, 89, 4280–4284. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Nikiforova, M.N.; Ciampi, R.; Salvatore, G.; Santoro, M.; Gandhi, M.; Knauf, J.A.; Thomas, G.A.; Jeremiah, S.; Bogdanova, T.I.; Tronko, M.D.; et al. Low prevalence of BRAF mutations in radiation-induced thyroid tumors in contrast to sporadic papillary carcinomas. Cancer Lett. 2004, 209, 1–6. [Google Scholar] [CrossRef]
  15. Powell, N.; Jeremiah, S.; Morishita, M.; Dudley, E.; Bethel, J.; Bogdanova, T.; Tronko, M.; Thomas, G. Frequency of BRAF T1796A mutation in papillary thyroid carcinoma relates to age of patient at diagnosis and not to radiation exposure. J. Pathol. 2005, 205, 558–564. [Google Scholar] [CrossRef]
  16. Morton, L.M.; Karyadi, D.M.; Stewart, C.; Bogdanova, T.I.; Dawson, E.T.; Steinberg, M.K.; Dai, J.; Hartley, S.W.; Schonfeld, S.J.; Sampson, J.N.; et al. Radiation-related genomic profile of papillary thyroid carcinoma after the Chernobyl accident. Science 2021, 372, 725–729. [Google Scholar] [CrossRef] [PubMed]
  17. Rogounovitch, T.I.; Mankovskaya, S.V.; Fridman, M.V.; Leonova, T.A.; Kondratovitch, V.A.; Konoplya, N.E.; Yamashita, S.; Mitsutake, N.; Saenko, V.A. Major Oncogenic Drivers and Their Clinicopathological Correlations in Sporadic Childhood Papillary Thyroid Carcinoma in Belarus. Cancers 2021, 13, 3374. [Google Scholar] [CrossRef] [PubMed]
  18. Nakayama, H.; Yoshida, A.; Nakamura, Y.; Hayashi, H.; Miyagi, Y.; Wada, N.; Rino, Y.; Masuda, M.; Imada, T. Clinical significance of BRAF (V600E) mutation and Ki-67 labeling index in papillary thyroid carcinomas. Anticancer Res. 2007, 27, 3645–3649. [Google Scholar]
  19. Nikiforova, M.N.; Kimura, E.T.; Gandhi, M.; Biddinger, P.W.; Knauf, J.A.; Basolo, F.; Zhu, Z.; Giannini, R.; Salvatore, G.; Fusco, A.; et al. BRAF mutations in thyroid tumors are restricted to papillary carcinomas and anaplastic or poorly differentiated carcinomas arising from papillary carcinomas. J. Clin. Endocrinol. Metab. 2003, 88, 5399–5404. [Google Scholar] [CrossRef]
  20. Rosenbaum, E.; Hosler, G.; Zahurak, M.; Cohen, Y.; Sidransky, D.; Westra, W.H. Mutational activation of BRAF is not a major event in sporadic childhood papillary thyroid carcinoma. Mod. Pathol. 2005, 18, 898–902. [Google Scholar] [CrossRef]
  21. Trovisco, V.; Soares, P.; Preto, A.; de Castro, I.V.; Lima, J.; Castro, P.; Maximo, V.; Botelho, T.; Moreira, S.; Meireles, A.M.; et al. Type and prevalence of BRAF mutations are closely associated with papillary thyroid carcinoma histotype and patients’ age but not with tumour aggressiveness. Virchows Arch. 2005, 446, 589–595. [Google Scholar] [CrossRef]
  22. Fugazzola, L.; Puxeddu, E.; Avenia, N.; Romei, C.; Cirello, V.; Cavaliere, A.; Faviana, P.; Mannavola, D.; Moretti, S.; Rossi, S.; et al. Correlation between B-RAFV600E mutation and clinico-pathologic parameters in papillary thyroid carcinoma: Data from a multicentric Italian study and review of the literature. Endocr. Relat. Cancer 2006, 13, 455–464. [Google Scholar] [CrossRef]
  23. Kebebew, E.; Weng, J.; Bauer, J.; Ranvier, G.; Clark, O.H.; Duh, Q.Y.; Shibru, D.; Bastian, B.; Griffin, A. The prevalence and prognostic value of BRAF mutation in thyroid cancer. Ann. Surg. 2007, 246, 466–471. [Google Scholar] [CrossRef]
  24. Mitsutake, N.; Fukushima, T.; Matsuse, M.; Rogounovitch, T.; Saenko, V.; Uchino, S.; Ito, M.; Suzuki, K.; Suzuki, S.; Yamashita, S. BRAF(V600E) mutation is highly prevalent in thyroid carcinomas in the young population in Fukushima: A different oncogenic profile from Chernobyl. Sci. Rep. 2015, 5, 16976. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Nikita, M.E.; Jiang, W.; Cheng, S.M.; Hantash, F.M.; McPhaul, M.J.; Newbury, R.O.; Phillips, S.A.; Reitz, R.E.; Waldman, F.M.; Newfield, R.S. Mutational Analysis in Pediatric Thyroid Cancer and Correlations with Age, Ethnicity, and Clinical Presentation. Thyroid 2016, 26, 227–234. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Oishi, N.; Kondo, T.; Nakazawa, T.; Mochizuki, K.; Inoue, T.; Kasai, K.; Tahara, I.; Yabuta, T.; Hirokawa, M.; Miyauchi, A.; et al. Frequent BRAF (V600E) and Absence of TERT Promoter Mutations Characterize Sporadic Pediatric Papillary Thyroid Carcinomas in Japan. Endocr. Pathol. 2017, 28, 103–111. [Google Scholar] [CrossRef]
  27. Yan, C.; Huang, M.; Li, X.; Wang, T.; Ling, R. Relationship between BRAF V600E and clinical features in papillary thyroid carcinoma. Endocr. Connect. 2019, 8, 988–996. [Google Scholar] [CrossRef] [Green Version]
  28. Iwadate, M.; Mitsutake, N.; Matsuse, M.; Fukushima, T.; Suzuki, S.; Matsumoto, Y.; Ookouchi, C.; Mizunuma, H.; Nakamura, I.; Nakano, K.; et al. The Clinicopathological Results of Thyroid Cancer with BRAFV600E Mutation in the Young Population of Fukushima. J. Clin. Endocrinol. Metab. 2020, 105, e4328–e4336. [Google Scholar] [CrossRef]
  29. Wang, Z.; Chen, J.Q.; Liu, J.L.; Qin, X.G. Clinical impact of BRAF mutation on the diagnosis and prognosis of papillary thyroid carcinoma: A systematic review and meta-analysis. Eur. J. Clin. Investig. 2016, 46, 146–157. [Google Scholar] [CrossRef] [PubMed]
  30. Ito, Y.; Miyauchi, A.; Kakudo, K.; Hirokawa, M.; Kobayashi, K.; Miya, A. Prognostic significance of ki-67 labeling index in papillary thyroid carcinoma. World J. Surg. 2010, 34, 3015–3021. [Google Scholar] [CrossRef]
  31. Lee, J.H.; Lee, E.S.; Kim, Y.S. Clinicopathologic significance of BRAF V600E mutation in papillary carcinomas of the thyroid: A meta-analysis. Cancer 2007, 110, 38–46. [Google Scholar] [CrossRef]
  32. Lupi, C.; Giannini, R.; Ugolini, C.; Proietti, A.; Berti, P.; Minuto, M.; Materazzi, G.; Elisei, R.; Santoro, M.; Miccoli, P.; et al. Association of BRAF V600E mutation with poor clinicopathological outcomes in 500 consecutive cases of papillary thyroid carcinoma. J. Clin. Endocrinol. Metab. 2007, 92, 4085–4090. [Google Scholar] [CrossRef] [Green Version]
  33. Kim, S.J.; Lee, K.E.; Myong, J.P.; Park, J.H.; Jeon, Y.K.; Min, H.S.; Park, S.Y.; Jung, K.C.; Koo, D.H.; Youn, Y.K. BRAF V600E mutation is associated with tumor aggressiveness in papillary thyroid cancer. World J. Surg. 2012, 36, 310–317. [Google Scholar] [CrossRef]
  34. Li, C.; Lee, K.C.; Schneider, E.B.; Zeiger, M.A. BRAF V600E mutation and its association with clinicopathological features of papillary thyroid cancer: A meta-analysis. J. Clin. Endocrinol. Metab. 2012, 97, 4559–4570. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Ma, Y.J.; Deng, X.L.; Li, H.Q. BRAF(V600E) mutation and its association with clinicopathological features of papillary thyroid microcarcinoma: A meta-analysis. J. Huazhong Univ. Sci. Technolog. Med. Sci. 2015, 35, 591–599. [Google Scholar] [CrossRef] [PubMed]
  36. Tufano, R.P.; Teixeira, G.V.; Bishop, J.; Carson, K.A.; Xing, M. BRAF mutation in papillary thyroid cancer and its value in tailoring initial treatment: A systematic review and meta-analysis. Medicine 2012, 91, 274–286. [Google Scholar] [CrossRef] [PubMed]
  37. Chen, Y.; Sadow, P.M.; Suh, H.; Lee, K.E.; Choi, J.Y.; Suh, Y.J.; Wang, T.S.; Lubitz, C.C. BRAF(V600E) Is Correlated with Recurrence of Papillary Thyroid Microcarcinoma: A Systematic Review, Multi-Institutional Primary Data Analysis, and Meta-Analysis. Thyroid 2016, 26, 248–255. [Google Scholar] [CrossRef]
  38. Liu, C.; Chen, T.; Liu, Z. Associations between BRAF(V600E) and prognostic factors and poor outcomes in papillary thyroid carcinoma: A meta-analysis. World J. Surg. Oncol. 2016, 14, 241. [Google Scholar] [CrossRef] [Green Version]
  39. Zhang, Q.; Liu, S.Z.; Zhang, Q.; Guan, Y.X.; Chen, Q.J.; Zhu, Q.Y. Meta-Analyses of Association Between BRAF(V600E) Mutation and Clinicopathological Features of Papillary Thyroid Carcinoma. Cell Physiol. Biochem. 2016, 38, 763–776. [Google Scholar] [CrossRef] [PubMed]
  40. Vuong, H.G.; Duong, U.N.; Altibi, A.M.; Ngo, H.T.; Pham, T.Q.; Tran, H.M.; Gandolfi, G.; Hassell, L. A meta-analysis of prognostic roles of molecular markers in papillary thyroid carcinoma. Endocr. Connect. 2017, 6, R8–R17. [Google Scholar] [CrossRef] [Green Version]
  41. Siironen, P.; Nordling, S.; Louhimo, J.; Haapiainen, R.; Haglund, C. Immunohistochemical expression of Bcl-2, Ki-67, and p21 in patients with papillary thyroid cancer. Tumour Biol. 2005, 26, 50–56. [Google Scholar] [CrossRef] [PubMed]
  42. Wang, M.; Wu, W.D.; Chen, G.M.; Chou, S.L.; Dai, X.M.; Xu, J.M.; Peng, Z.H. Could Tumor Size Be a Predictor for Papillary Thyroid Microcarcinoma: A Retrospective Cohort Study. Asian Pac. J. Cancer Prev. 2015, 16, 8625–8628. [Google Scholar] [CrossRef]
  43. Zhou, Y.; Jiang, H.G.; Lu, N.; Lu, B.H.; Chen, Z.H. Expression of ki67 in papillary thyroid microcarcinoma and its clinical significance. Asian Pac. J. Cancer Prev. 2015, 16, 1605–1608. [Google Scholar] [CrossRef] [Green Version]
  44. Avdalyan, A.M.; Ivanov, A.A.; Lushnikova, E.L.; Molodykh, O.P.; Vikhlyanov, I.V. The Relationship of Immunoexpression of Ki-67 and Hsp70 with Clinical and Morphological Parameters and Prognosis of Papillary Thyroid Cancer. Bull. Exp. Biol. Med. 2020, 168, 688–693. [Google Scholar] [CrossRef] [PubMed]
  45. Miyauchi, A.; Kudo, T.; Hirokawa, M.; Ito, Y.; Kihara, M.; Higashiyama, T.; Yabuta, T.; Masuoka, H.; Shindo, H.; Kobayashi, K.; et al. Ki-67 labeling index is a predictor of postoperative persistent disease and cancer growth and a prognostic indicator in papillary thyroid carcinoma. Eur. Thyroid J. 2013, 2, 57–64. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Hirokawa, M.; Kudo, T.; Ota, H.; Suzuki, A.; Miyauchi, A. Pathological characteristics of low-risk papillary thyroid microcarcinoma with progression during active surveillance. Endocr. J. 2016, 63, 805–810. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Tanaka, A.; Matsuse, M.; Saenko, V.; Nakao, T.; Yamanouchi, K.; Sakimura, C.; Yano, H.; Nishihara, E.; Hirokawa, M.; Suzuki, K.; et al. TERT mRNA Expression as a Novel Prognostic Marker in Papillary Thyroid Carcinomas. Thyroid 2019, 29, 1105–1114. [Google Scholar] [CrossRef]
  48. Ma, B.; Wang, Y.; Yang, S.; Ji, Q. Predictive factors for central lymph node metastasis in patients with cN0 papillary thyroid carcinoma: A systematic review and meta-analysis. Int. J. Surg. 2016, 28, 153–161. [Google Scholar] [CrossRef] [PubMed]
  49. Sassolas, G.; Hafdi-Nejjari, Z.; Ferraro, A.; Decaussin-Petrucci, M.; Rousset, B.; Borson-Chazot, F.; Borbone, E.; Berger, N.; Fusco, A. Oncogenic alterations in papillary thyroid cancers of young patients. Thyroid 2012, 22, 17–26. [Google Scholar] [CrossRef]
  50. Givens, D.J.; Buchmann, L.O.; Agarwal, A.M.; Grimmer, J.F.; Hunt, J.P. BRAF V600E does not predict aggressive features of pediatric papillary thyroid carcinoma. Laryngoscope 2014, 124, E389–E393. [Google Scholar] [CrossRef]
  51. Alzahrani, A.S.; Murugan, A.K.; Qasem, E.; Alswailem, M.; Al-Hindi, H.; Shi, Y. Single Point Mutations in Pediatric Differentiated Thyroid Cancer. Thyroid 2017, 27, 189–196. [Google Scholar] [CrossRef]
  52. Geng, J.; Wang, H.; Liu, Y.; Tai, J.; Jin, Y.; Zhang, J.; He, L.; Fu, L.; Qin, H.; Song, Y.; et al. Correlation between BRAF (V600E) mutation and clinicopathological features in pediatric papillary thyroid carcinoma. Sci. China Life Sci. 2017, 60, 729–738. [Google Scholar] [CrossRef]
  53. Mostoufi-Moab, S.; Labourier, E.; Sullivan, L.; LiVolsi, V.; Li, Y.; Xiao, R.; Beaudenon-Huibregtse, S.; Kazahaya, K.; Adzick, N.S.; Baloch, Z.; et al. Molecular Testing for Oncogenic Gene Alterations in Pediatric Thyroid Lesions. Thyroid 2018, 28, 60–67. [Google Scholar] [CrossRef]
  54. Midorikawa, S.; Ohtsuru, A.; Murakami, M.; Takahashi, H.; Suzuki, S.; Matsuzuka, T.; Shimura, H.; Ohira, T.; Suzuki, S.I.; Yasumura, S.; et al. Comparative Analysis of the Growth Pattern of Thyroid Cancer in Young Patients Screened by Ultrasonography in Japan After a Nuclear Accident: The Fukushima Health Management Survey. JAMA Otolaryngol. Head Neck Surg. 2018, 144, 57–63. [Google Scholar] [CrossRef]
  55. Bogdanova, T.I.; Zurnadzhy, L.Y.; Nikiforov, Y.E.; Leeman-Neill, R.J.; Tronko, M.D.; Chanock, S.; Mabuchi, K.; Likhtarov, I.A.; Kovgan, L.M.; Drozdovitch, V.; et al. Histopathological features of papillary thyroid carcinomas detected during four screening examinations of a Ukrainian-American cohort. Br. J. Cancer 2015, 113, 1556–1564. [Google Scholar] [CrossRef] [Green Version]
  56. Smida, J.; Salassidis, K.; Hieber, L.; Zitzelsberger, H.; Kellerer, A.M.; Demidchik, E.P.; Negele, T.; Spelsberg, F.; Lengfelder, E.; Werner, M.; et al. Distinct frequency of ret rearrangements in papillary thyroid carcinomas of children and adults from Belarus. Int. J. Cancer 1999, 80, 32–38. [Google Scholar] [CrossRef]
  57. Wang, Y.L.; Zhang, R.M.; Luo, Z.W.; Wu, Y.; Du, X.; Wang, Z.Y.; Zhu, Y.X.; Li, D.S.; Ji, Q.H. High frequency of level II-V lymph node involvement in RET/PTC positive papillary thyroid carcinoma. Eur. J. Surg. Oncol. 2008, 34, 77–81. [Google Scholar] [CrossRef] [PubMed]
  58. Yamashita, S.; Saenko, V. Mechanisms of Disease: Molecular genetics of childhood thyroid cancers. Nat. Clin. Pract. Endocrinol. Metab. 2007, 3, 422–429. [Google Scholar] [CrossRef] [PubMed]
  59. Tronko, M.; Shpak, V.; Bogdanova, T.; Saenko, V.; Yamashita, S. Epidemiology of thyroid cancer in Ukraine after Chernobyl. In Thyroid Cancer in Ukraine after Chernobyl: Dosimetry, Epidemiology, Pathology, Molecular Biology; Tronko, M., Bogdanova, T., Saenko, V., Thomas, G.A., Likhtarov, I., Yamashita, S., Eds.; IN-TEX: Nagasaki, Japan, 2014; pp. 39–64. [Google Scholar]
  60. Bogdanova, T.I.; Saenko, V.A.; Hashimoto, Y.; Hirokawa, M.; Zurnadzhy, L.Y.; Hayashi, T.; Ito, M.; Iwadate, M.; Mitsutake, N.; Rogounovitch, T.I.; et al. Papillary Thyroid Carcinoma in Ukraine After Chernobyl and in Japan After Fukushima: Different Histopathological Scenarios. Thyroid 2021, 31, 1322–1334. [Google Scholar] [CrossRef]
  61. Likhtarov, I.; Thomas, G.; Kovgan, L.; Masiuk, S.; Chepurny, M.; Ivanova, O.; Gerasymenko, V.; Tronko, M.; Bogdanova, T.; Bouville, A. Reconstruction of individual thyroid doses to the Ukrainian subjects enrolled in the Chernobyl Tissue Bank. Radiat. Prot. Dosim. 2013, 156, 407–423. [Google Scholar] [CrossRef] [PubMed]
  62. Lloyd, R.V.; Osamura, R.Y.; Kloppel, G.; Rosai, J. WHO Classification of Tumours of Endocrine Organs, 4th ed.; IARC Press: Lyon, France, 2017. [Google Scholar]
  63. Thomas, G.A.; Williams, E.D.; Becker, D.V.; Bogdanova, T.I.; Demidchik, E.P.; Lushnikov, E.; Nagataki, S.; Ostapenko, V.; Pinchera, A.; Souchkevitch, G.; et al. Chernobyl tumor bank. Thyroid 2000, 10, 1126–1127. [Google Scholar] [CrossRef]
  64. Thomas, G.A. The Chernobyl Tissue Bank: Integrating research on radiation-induced thyroid cancer. J. Radiol. Prot. 2012, 32, N77–N80. [Google Scholar] [CrossRef] [PubMed]
  65. Brierley, J.D.; Gospodarowich, M.K.; Wittekind, C. TNM Classification of Malignant Tumours, 8th ed.; Wiley-Blackwell: Oxford, UK, 2017. [Google Scholar]
  66. Parker, K.G.; White, M.G.; Cipriani, N.A. Comparison of Molecular Methods and BRAF Immunohistochemistry (VE1 Clone) for the Detection of BRAF V600E Mutation in Papillary Thyroid Carcinoma: A Meta-Analysis. Head Neck Pathol. 2020, 14, 1067–1079. [Google Scholar] [CrossRef]
  67. Nakao, T.; Matsuse, M.; Saenko, V.; Rogounovitch, T.; Tanaka, A.; Suzuki, K.; Higuchi, M.; Sasai, H.; Sano, T.; Hirokawa, M.; et al. Preoperative detection of the TERT promoter mutations in papillary thyroid carcinomas. Clin. Endocrinol. 2021, 95, 790–799. [Google Scholar] [CrossRef] [PubMed]
  68. Thomas, L.; Reyes, E.M. Tutorial: Survival Estimation for Cox Regression Models with Time-Varying Coefficients Using SAS and R. J. Stat. Softw. 2014, 61, 1–23. [Google Scholar] [CrossRef] [Green Version]
  69. Nenadic, O.; Greenacre, M. Correspondence Analysis in R, with two- and three-dimensional graphics: The ca package. J. Stat. Softw. 2007, 20, 1–13. [Google Scholar]
Cancers 13 06038 g001 550
Figure 1. Immunohistochemical staining for BRAFV600E and Ki67. (a–c) the BRAFV600E-positive, and (d–f) the BRAFV600E-negative radiogenic PTCs. (a)—papillary dominant growth pattern, hematoxylin-eosin, ×200; (b)—positive IHC reaction with the VE1 anti-BRAF (mutated V600E) antibody, ×200; (c)—IHC reaction with Ki67 (Clone MIB-1) antibody (Ki67 LI 8.7%), ×200; (d)—solid dominant growth pattern, hematoxylin-eosin, ×200; (e)—negative IHC reaction with the VE1 anti-BRAF (mutated V600E) antibody,×200; (f)—IHC reaction with Ki67 (Clone MIB-1) antibody (Ki67 LI 2.4%), ×200.
Figure 1. Immunohistochemical staining for BRAFV600E and Ki67. (a–c) the BRAFV600E-positive, and (d–f) the BRAFV600E-negative radiogenic PTCs. (a)—papillary dominant growth pattern, hematoxylin-eosin, ×200; (b)—positive IHC reaction with the VE1 anti-BRAF (mutated V600E) antibody, ×200; (c)—IHC reaction with Ki67 (Clone MIB-1) antibody (Ki67 LI 8.7%), ×200; (d)—solid dominant growth pattern, hematoxylin-eosin, ×200; (e)—negative IHC reaction with the VE1 anti-BRAF (mutated V600E) antibody,×200; (f)—IHC reaction with Ki67 (Clone MIB-1) antibody (Ki67 LI 2.4%), ×200.
Cancers 13 06038 g001
Cancers 13 06038 g002 550
Figure 2. Correspondence analysis of the associations of the BRAFV600E-positive and BRAFV600E-negative PTCs of different etiology with major histopathological characteristics and tumor invasive features. The biplot displays column variables (the four PTC groups) in principal coordinates and row variables (categorical clinicopathological variables) in contribution coordinates. These coordinates present, to some extent, the association between column and row variables (e.g., the BRAFV600E-positive PTCs have a papillary growth pattern more frequently than the BRAFV600E-negative PTCs). Arrows are shown for the four PTC groups (column variables) in the graph for easier visualization of the angles between them. The smaller angle between column points (when connected to the origin) indicates the stronger correlation (e.g., the radiogenic and sporadic BRAFV600E-positive PTCs are rather similar), right or obtuse angle indicates the points are uncorrelated (e.g., the radiogenic BRAFV600E-negative PTCs are quite different from all other groups). For row points (i.e., clinicopathological characteristics), the smaller angle indicates similarity in response pattern (e.g., M0 and N0 tumors would be expected to frequently coexist, and are likely to occur in tumors with invasiveness scores 0 and 1). Dimensions 1 and 2 accounted for 96.4% of variance, and Dimension 3—for the remaining 3.6%.
Figure 2. Correspondence analysis of the associations of the BRAFV600E-positive and BRAFV600E-negative PTCs of different etiology with major histopathological characteristics and tumor invasive features. The biplot displays column variables (the four PTC groups) in principal coordinates and row variables (categorical clinicopathological variables) in contribution coordinates. These coordinates present, to some extent, the association between column and row variables (e.g., the BRAFV600E-positive PTCs have a papillary growth pattern more frequently than the BRAFV600E-negative PTCs). Arrows are shown for the four PTC groups (column variables) in the graph for easier visualization of the angles between them. The smaller angle between column points (when connected to the origin) indicates the stronger correlation (e.g., the radiogenic and sporadic BRAFV600E-positive PTCs are rather similar), right or obtuse angle indicates the points are uncorrelated (e.g., the radiogenic BRAFV600E-negative PTCs are quite different from all other groups). For row points (i.e., clinicopathological characteristics), the smaller angle indicates similarity in response pattern (e.g., M0 and N0 tumors would be expected to frequently coexist, and are likely to occur in tumors with invasiveness scores 0 and 1). Dimensions 1 and 2 accounted for 96.4% of variance, and Dimension 3—for the remaining 3.6%.
Cancers 13 06038 g002
Table
Table 1. Frequency of BRAFV600E in radiogenic and sporadic PTC in different age groups at the time of operation.
Table 1. Frequency of BRAFV600E in radiogenic and sporadic PTC in different age groups at the time of operation.
Age GroupsRadiogenic PTC (n = 247)Sporadic PTC (n = 138)p-Value 1
Number%Number%
Children aged ≤14 years1/1041.05/3912.80.006
Adolescents aged 15–18 years6/5211.512/3732.40.030
Young adults aged 19–28 years19/9120.925/6240.30.011
Total26/24710.542/13830.42.67 × 10−6
Age trend, p-value 25.97 × 10−60.004
Age association 31.179 (1.092–1.273); 2.40 × 10−51.118 (1.031–1.212); 0.007
1 The Fisher’s exact test. 2 The Cochran-Armitage test for trend, two-sided. 3 OR (95% CI) adjusted for sex; p-value. Numbers in bold indicate statistical significance.
Table
Table 2. Characteristics of the BRAFV600E-positive and BRAFV600E-negative radiogenic PTCs.
Table 2. Characteristics of the BRAFV600E-positive and BRAFV600E-negative radiogenic PTCs.
ParametersBRAFV600E(+) (n = 26)BRAFV600E(−) (n = 221)p-ValueOR, b or HR (95%CI)p-Value
Number or Value (% or IQR)Number or Value (% or IQR)UnivariateMultivariate 1
Sex (F/M, %M, F:M ratio; ref = F)20/6 (23.1%; 3.3:1)141/80 (36.2%; 1.8:1)0.2760.706 (0.260–1.914) 20.494
Age at operation, years24.1 (18.7–27.3)15.4 (12.0–21.3)2.03 × 10−66.012 (3.552–8.472)33.00 × 10−6
Age at exposure, years2.0 (1.0–3.0)2.0 (1.0–3.0)0.924−0.019 (−0.519–0.481) 30.941
Period of latency, years22.1 (17.5–24.9)13.1 (9.3–18.5)1.24 × 10−66.133 (3.679–8.588)32.00 × 10−6
Radiation dose to the thyroid, mGy200 (134–390)350 (173–825)0.019−0.295 (−0.563–0.027)30.031
Tumor size, mm, median10 (6–20)16 (12–30)4.57 × 10−4−0.528 (−0.943–0.113)0.013
≤10 mm (microcarcinoma)14 (53.8%)36 (16.3%)5.44 × 10−53.475 (1.380–8.751)0.008
Complete tumor capsule5 (19.2%)32 (14.5%)0.5600.721 (0.239–2.174)0.561
Dominant growth pattern 1.98 × 10−60.334 (0.180–0.622)0.001
papillary17 (65.4%)40 (18.1%)1.12 × 10−65.735 (2.289–14.366)1.93 × 10−4
follicular5 (19.2%)75 (33.9%)0.1830.410 (0.141–1.190)0.101
solid-trabecular 4 (15.4%)106 (48.0%)0.0010.289 (0.093–0.901)0.032
Oncocytic changes8 (30.8%)31 (14.0%)0.0421.671 (0.629–4.441)0.303
Multifocality6 (23.1%)30 (13.6%)0.2351.124 (0.395–3.194)0.827
Lymphatic/vascular invasion9 (34.6%)159 (71.9%)1.91 × 10−40.417 (0.165–1.052)0.064
Extrathyroidal extension (any)7 (26.9%)119 (53.8%)0.0120.565 (0.216–1.473)0.243
N category (N1)9 (34.6%)122 (55.2%)0.0610.649 (0.264–1.592)0.345
N1a6 (23.1%)43 (19.5%)0.6111.210 (0.434–3.376)0.715
N1b3 (11.5%)79 (35.7%)0.0140.392 (0.109–1.406)0.151
M category (M1)1 (3.8%)38 (17.2%)0.0910.340 (0.042–2.723)0.309
pT 0.0780.605 (0.310–1.176)0.139
pT121 (80.8%)126 (57.0%)0.0212.437 (0.854–6.959)0.096
pT1a14 (53.8%)36 (16.3%)5.40 × 10−53.475 (1.380–8.751)0.008
pT1b7 (26.9%)90 (40.7%)0.2060.658 (0.255–1.700)0.387
pT22 (7.7%)42 (19.0%)0.1850.338 (0.074–1.542)0.161
pT33 (11.5%)53 (24.0%)0.2160.640 (0.174–2.350)0.501
pT3a1 (3.8%)17 (7.7%)0.7030.634 (0.075–5.385)0.676
pT3b2 (7.7%)36 (16.3%)0.3890.678 (0.144–3.201)0.623
Invasiveness score1 (0–2)2 (1–3)0.0720.749 (0.531–1.056)0.100
010 (38.5%)39 (17.6%)0.0181.836 (0.737–4.575)0.192
14 (15.4%)37 (16.7%)1.0000.603 (0.187–1.941)0.396
28 (30.8%)60 (27.1%)0.6511.306 (0.511–3.339)0.577
34 (15.4%)44 (19.9%)0.7941.321 (0.394–4.423)0.652
4029 (13.1%)0.0520.148 (0.008–2.724)0.199
5012 (5.4%)0.6210.444 (0.017–11.596)0.626
Ki-67 labeling index, median4.7 (3.8–6.3)n = 210; 2.3 (1.3–3.9)7.64 × 10−61.647 (0.467–2.827)0.006
0–5%13 (50.0%)173 (78.3%)0.0030.268 (0.109–0.659)0.004
>5–10%11 (42.3%)29 (23.1%)7.00 × 10−43.974 (1.565–10.095)0.004
>10%2 (7.7%)8 (3.6%)0.2841.351 (0.252–7.253)0.726
Follow-up, years9.9 (8.3–12.2)13.5 (8.5–18.7)0.0160.469 (−1.978–2.916)0.706
LN recurrence (reoperated after 6 mo)0n = 210; 8 (3.8%)0.574 40.523 (0.004–5.410) 50.708
1 Adjusted for age at operation and sex unless otherwise specified; characteristics of sporadic PTCs were used as references; OR from logistic regression, b-coefficient from linear regression, HR from proportional hazard (Cox) regression. 2 Adjusted for age at operation. 3 Adjusted for sex. 4 The log-rank test. 5 The Firth’s penalized proportional hazard model. Numbers in bold indicate statistical significance.
Table
Table 3. Characteristics of the BRAFV600E-positive and BRAFV600 E-negative sporadic PTCs.
Table 3. Characteristics of the BRAFV600E-positive and BRAFV600 E-negative sporadic PTCs.
ParametersBRAFV600E(+) (n = 42)BRAFV600E(−) (n = 96)p-ValueOR, b or HR (95%CI)p-Value
Number or Value (% or IQR)Number or Value (% or IQR)1UnivariateMultivariate 1
Sex (F/M, %M, F:M ratio; ref = F)34/8 (19.0%; 4.3:1)81/15 (15.6%; 5.4:1)0.6261.378 (0.518–3.670) 20.521
Age at operation, years21.0 (16.4–24.3)17.1 (14.0–21.7)0.0062.516 (0.748–4.284)30.006
Tumor size, mm11 (8–15)21 (13–31)2.64 × 10−5−8.706 (−13.122–4.290)1.52 × 10−4
≤10 mm (microcarcinoma)19 (45.2%)18 (18.8%)0.0023.499 (1.531–7.966)0.003
Complete tumor capsule6 (14.3%)32 (33.3%)0.0230.231 (0.083–0.648)0.005
Dominant growth pattern 1.72 × 10−40.427 (0.248–0.736)0.002
papillary31 (73.8%)35 (36.5%)7.91 × 10−55.462 (2.355–12.666)7.60 × 10−5
follicular4 (9.5%)33 (34.4%)0.0030.143 (0.045–0.460)0.001
solid-trabecular 7 (16.7%)28 (29.2%)0.1410.565 (0.219–1.460)0.238
Oncocytic changes26 (61.9%)36 (37.5%)0.0102.488 (1.155–5.358)0.020
Multifocality10 (23.8%)11 (11.5%)0.0751.905 (0.712–5.095)0.199
Lymphatic/vascular invasion13 (31.0%)62 (64.6%)7.46 × 10−60.304 (0.134–0.688)0.004
Extrathyroidal extension (any)11 (26.2%)33 (34.4%)0.5430.897 (0.382–2.106)0.803
N category (N1)15 (35.7%)42 (43.8%)0.4540.793 (0.362–1.739)0.562
N1a11 (26.2%)16 (16.7%)0.2441.696 (0.682–4.214)0.255
N1b4 (9.5%)26 (27.1%)0.0250.343 (0.109–1.084)0.068
M category (M1)1 (2.4%)5 (5.2%)0.4070.445 (0.049–3.998)0.469
pT 0.0780.259 (0.115–0.580)0.001
pT135 (83.3%)48 (50.0%)2.72 × 10−45.098 (2.017–12.883)0.001
pT1a19 (45.2%)18 (18.8%)0.0023.499 (1.531–7.996)0.003
pT1b16 (38.1%)30 (31.3%)0.4401.427 (0.651–3.130)0.374
pT26 (14.3%)36 (37.5%)0.0080.269 (0.101–0.719)0.009
pT31 (2.4%)12 (12.5%)0.1090.175 (0.021–1.422)0.103
pT3a1 (2.4%)9 (9.4%)0.2820.227 (0.027–1.904)0.172
pT3b03 (3.1%)0.5530.336 (0.009–12.891)0.558
Invasiveness score1 (0–2)2 (0–3)0.4240.835 (0.607–1.147)0.266
015 (35.7%)27 (28.1%)0.4231.215 (0.543–2.719)0.635
113 (31.0%)20 (20.8%)0.2031.601 (0.688–3.728)0.275
27 (16.7%)21 (21.9%)0.6460.705 (0.267–1.866)0.482
35 (11.9%)21 (21.9%)0.2370.506 (0.171–1.491)0.217
42 (4.8%)7 (7.3%)0.7221.175 (0.205–6.750)0.856
500NDNDND
Ki-67 labeling indexn = 40; 5.9 (4.3–8.7)4.1 (2.5–6.9)0.0032.221 (0.842–3.600)0.002
0–5%14 (35.0%)62 (64.6%)0.0020.231 (0.101–0.531)5.49 × 10−4
>5–10%20 (50.0%)25 (26.0%)0.0093.109 (1.395–6.930)0.006
>10%6 (15.0%)9 (9.4%)0.3742.322 (0.716–7.534)0.161
Follow-up, years4.9 (2.4–8.6)5.5 (2.6–9.3)0.343−0.144 (−1.787–1.499)0.863
LN recurrence (reoperated after 6 mo)n = 39; 3 (7.7%)n = 88; 1 (1.1%)0.007419.042 (1.299–279.067)50.031
1 Adjusted for age at operation and sex unless otherwise specified; characteristics of the sporadic PTCs were used as references; OR from logistic regression, b-coefficient from linear regression, HR from proportional hazard (Cox) regression. 2 Adjusted for age at operation. 3 Adjusted for sex. 4 The log-rank test. 5 The proportional hazard model. Numbers in bold indicate statistical significance.
Table
Table 4. Statistical comparison of the BRAFV600E-positive or BRAFV600E-negative PTCs across the radiogenic and sporadic.
Table 4. Statistical comparison of the BRAFV600E-positive or BRAFV600E-negative PTCs across the radiogenic and sporadic.
CharacteristicsBRAFV600E(+) (n = 26/42) 1BRAFV600E(−) (n = 221/96) 1
UnivariateOR, b or HR (95%CI) 2MultivariateUnivariateOR, b or HR (95%CI) 2Multivariate
p-Valuep-Valuep-Valuep-Value
Sex (ref = F)0.7621.357 (0.389–4.731) 30.6321.80 × 10−42.952 (1.590–5.482)30.001
Age at operation, years0.0152.838 (0.496–5.180)40.0180.070−0.593 (−2.015–0.829) 40.412
Tumor size, mm0.5190.264 (−4.730–5.258)0.9160.224−0.865 (−4.149–2.419)0.605
≤10 mm (microcarcinoma)0.6191.262 (0.447–3.561)0.6610.6270.781 (0.404–1.508)0.461
Complete tumor capsule0.7371.516 (0.383–5.998)0.5532.22 × 10−40.315 (0.171–0.581)2.15 × 10−4
Dominant growth pattern0.5390.955 (0.303–3.004)0.9375.88 × 10−42.340 (1.472–3.721)3.25 × 10−4
papillary0.5850.938 (0.290–3.034)0.9140.0010.416 (0.239–0.724)0.002
follicular0.2871.779 (0.388–8.147)0.45810.882 (0.522–1.493)0.641
solid-trabecular10.622 (0.140–2.755)0.5310.0022.350 (1.370–4.032)0.002
Oncocytic changes0.0240.242 (0.080–0.733)0.0129.00 × 10−60.298 (0.167–0.531)4.00 × 10−5
Multifocality10.612 (0.170–2.209)0.4540.7171.281 (0.600–2.737)0.523
Lymphatic/vascular invasion0.7941.915 (0.594–6.172)0.2770.2311.364 (0.787–2.366)0.269
Extrathyroidal extension10.946 (0.294–3.043)0.9250.0012.115 (1.244–3.595)0.006
minimal 10.804 (0.242–2.667)0.7210.2531.255 (0.737–2.138)0.403
gross0.3823.224 (0.137–60.094)0.4330.0015.219 (1.532–17.781)0.008
N category (N1)10.791 (0.264–2.372)0.6750.0671.438 (0.871–2.374)0.155
N1a10.683 (0.201–2.321)0.5410.6391.187 (0.623–2.265)0.602
N1b11.198 (0.229–6.262)0.8310.1541.327 (0.765–2.302)0.315
M category (M1)12.504 (0.128–49.133)0.5450.0552.139 (0.929–4.924)0.074
pT0.3061.462 (0.389–5.503)0.5740.0010.976 (0.608–1.567)0.920
pT110.770 (0.200–2.957)0.7030.2701.361 (0.830–2.231)0.222
pT1a0.6191.262 (0.447–3.581)0.6610.6270.781 (0.404–1.508)0.461
pT1b0.4330.649 (0.214–1.974)0.4470.1311.555 (0.923–2.621)0.097
pT20.7010.563 (0.099–3.211)0.5176.39 × 10−40.385 (0.223–0.666)6.49 × 10−4
pT30.1525.074 (0.465–55.415)0.1830.0232.155 (1.072–4.335)0.031
pT3a12.504 (0.128–49.133)0.5450.6580.842 (0.354–2.006)0.698
pT3b0.1431.476 (0.036–60.523)0.8376.30 × 10−45.599 (1.651–18.985)0.006
Invasiveness score (any Ex)0.4011.063 (0.422–2.678)0.8970.0241.644 (1.059–2.552)0.027
011.105 (0.379–3.222)0.8540.0500.549 (0.305–0.990)0.046
10.2490.482 (0.133–1.754)0.2680.4270.816 (0.436–1.526)0.524
20.2311.940 (0.574–6.558)0.2860.4011.378 (0.770–2.464)0.28
30.7231.036 (0.228–4.714)0.9640.7620.713 (0.384–1.323)0.283
40.5210.495 (0.027–8.990)0.6340.1771.918 (0.792–4.644)0.149
5ND 5NDND0.02110.119 (0.618–165.760)0.105
Ki-67 labeling index0.112 6−1.114 (−3.135–0.907)0.2753.84 × 10−9 7−1.714 (−2.455–0.973)8.00 × 10−6
0–5%0.3062.057 (0.702–6.031)0.1890.0012.283 (1.304–3.995)0.004
>5–10%0.6180.560 (0.191–1.645)0.2920.0150.488 (0.264–0.903)0.022
>10%0.4640.710 (0.117–4.312)0.710.0610.473 (0.174–1.283)0.142
Follow-up, years (n = 26/39)3.39 × 10−513.324 (8.434–18.213)8.82 × 10−71.50 × 10−146.413 (5.048–7.777)3.81 × 10−18
LN recurrences (reoperated after 6 mo)0.02780.031 (0.000–0.870) 90.1470.872 80.694 (0.074–6.535) 100.750
1 The number of radiogenic/sporadic cases. 2 Adjusted for age at operation and sex unless otherwise specified; characteristics of the sporadic PTCs were used as references; OR from logistic regression, b-coefficient from linear regression, HR from proportional hazard (Cox) regression. 3 Adjusted for age at operation. 4 Adjusted for sex. 5 Not determined. 6 The number of cases is 26/40. 7 The number of cases is 210/96. 8 The log-rank test. 9 The Firth’s penalized proportional hazard model. 10 The proportional hazard model. Numbers in bold indicate statistical significance.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zurnadzhy, L.; Bogdanova, T.; Rogounovitch, T.I.; Ito, M.; Tronko, M.; Yamashita, S.; Mitsutake, N.; Chernyshov, S.; Masiuk, S.; Saenko, V.A. The BRAFV600E Mutation Is Not a Risk Factor for More Aggressive Tumor Behavior in Radiogenic and Sporadic Papillary Thyroid Carcinoma at a Young Age. Cancers 2021, 13, 6038. https://doi.org/10.3390/cancers13236038

AMA Style

Zurnadzhy L, Bogdanova T, Rogounovitch TI, Ito M, Tronko M, Yamashita S, Mitsutake N, Chernyshov S, Masiuk S, Saenko VA. The BRAFV600E Mutation Is Not a Risk Factor for More Aggressive Tumor Behavior in Radiogenic and Sporadic Papillary Thyroid Carcinoma at a Young Age. Cancers. 2021; 13(23):6038. https://doi.org/10.3390/cancers13236038

Chicago/Turabian Style

Zurnadzhy, Liudmyla, Tetiana Bogdanova, Tatiana I. Rogounovitch, Masahiro Ito, Mykola Tronko, Shunichi Yamashita, Norisato Mitsutake, Serhii Chernyshov, Sergii Masiuk, and Vladimir A. Saenko. 2021. "The BRAFV600E Mutation Is Not a Risk Factor for More Aggressive Tumor Behavior in Radiogenic and Sporadic Papillary Thyroid Carcinoma at a Young Age" Cancers 13, no. 23: 6038. https://doi.org/10.3390/cancers13236038

APA Style

Zurnadzhy, L., Bogdanova, T., Rogounovitch, T. I., Ito, M., Tronko, M., Yamashita, S., Mitsutake, N., Chernyshov, S., Masiuk, S., & Saenko, V. A. (2021). The BRAFV600E Mutation Is Not a Risk Factor for More Aggressive Tumor Behavior in Radiogenic and Sporadic Papillary Thyroid Carcinoma at a Young Age. Cancers, 13(23), 6038. https://doi.org/10.3390/cancers13236038

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

Article Metrics

Back to TopTop