Next Article in Journal
The Role of ROS as a Double-Edged Sword in (In)Fertility: The Impact of Cancer Treatment
Next Article in Special Issue
Medullary Thyroid Carcinoma Mutational Spectrum Update and Signaling-Type Inference by Transcriptional Profiles: Literature Meta-Analysis and Study of Tumor Samples
Previous Article in Journal
Immunophenotypic Analysis of Acute Megakaryoblastic Leukemia: A EuroFlow Study
Previous Article in Special Issue
Gene Editing with CRISPR/Cas Methodology and Thyroid Cancer: Where Are We?
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cancers
Volume 14
Issue 6
10.3390/cancers14061584
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

APOBEC SBS13 Mutational Signature—A Novel Predictor of Radioactive Iodine Refractory Papillary Thyroid Carcinoma

by
Sarah Siraj
1,
Tariq Masoodi
1,
Abdul K. Siraj
1,
Saud Azam
1,
Zeeshan Qadri
1,
Sandeep K. Parvathareddy
1,
Rong Bu
1,
Khawar S. Siddiqui
2,
Saif S. Al-Sobhi
3,
Mohammed AlDawish
4 and
Khawla S. Al-Kuraya
1,*
1
Human Cancer Genomic Research, King Faisal Specialist Hospital and Research Centre, P.O. Box 3354, Riyadh 11211, Saudi Arabia
2
Department of Pediatric Hematology-Oncology, King Faisal Specialist Hospital and Research Center, P.O. Box 3354, Riyadh 11211, Saudi Arabia
3
Department of Surgery, King Faisal Specialist Hospital and Research Center, P.O. Box 3354, Riyadh 11211, Saudi Arabia
4
Department of Endocrinology and Diabetes, Prince Sultan Military Medical City, P.O. Box 261370, Riyadh 11342, Saudi Arabia
*
Author to whom correspondence should be addressed.
Cancers 2022, 14(6), 1584; https://doi.org/10.3390/cancers14061584
Submission received: 7 February 2022 / Revised: 3 March 2022 / Accepted: 15 March 2022 / Published: 21 March 2022
(This article belongs to the Special Issue Advances in Thyroid Carcinoma)
Browse Figures
Review Reports Versions Notes

Abstract

:

Simple Summary

Around 15% of papillary thyroid carcinoma (PTC) patients are not cured using standard surgery followed by radioactive iodine (RAI) therapy, and instead develop refractory disease. The aim of this study was to help identify RAI-refractory PTC patients early and guide precision medicine by performing a clinical and genomic characterization of RAI-refractory and avid PTCs. RAI-refractory PTCs had a more aggressive clinical presentation, a higher number of mutations, harbored more TERT promoter (TERTp) mutations, and were enriched with APOBEC-related mutations. Notably, the APOBEC single-base substitution (SBS) mutational signature, SBS13, and TERTp mutations were revealed to be independent predictors of RAI refractoriness in PTC. Although SBS13 and TERTp mutations alone highly predicted RAI refractoriness, when combined, they formed a stronger predictor of RAI refractoriness in PTC. This study highlights the APOBEC SBS13 mutational signature as a novel independent predictor of RAI refractoriness in a distinct subgroup of PTC.

Abstract

Standard surgery followed by radioactive iodine (131I, RAI) therapy are not curative for 5–20% of papillary thyroid carcinoma (PTC) patients with RAI refractory disease. Early predictors indicating therapeutic response to RAI therapy in PTC are yet to be elucidated. Whole-exome sequencing was performed (at median depth 198x) on 66 RAI-refractory and 92 RAI-avid PTCs with patient-matched germline. RAI-refractory tumors were significantly associated with distinct aggressive clinicopathological features, including positive surgical margins (p = 0.016) and the presence of lymph node metastases at primary diagnosis (p = 0.012); higher nonsilent tumor mutation burden (p = 0.011); TERT promoter (TERTp) mutation (p < 0.0001); and the enrichment of the APOBEC-related single-base substitution (SBS) COSMIC mutational signatures 2 (p = 0.030) and 13 (p < 0.001). Notably, SBS13 (odds ratio [OR] 30.4, 95% confidence intervals [CI] 1.43–647.22) and TERTp mutation (OR 41.3, 95% CI 4.35–391.60) were revealed to be independent predictors of RAI refractoriness in PTC (p = 0.029 and 0.001, respectively). Although SBS13 and TERTp mutations alone highly predicted RAI refractoriness, when combined, they significantly increased the likelihood of predicting RAI refractoriness in PTC. This study highlights the APOBEC SBS13 mutational signature as a novel independent predictor of RAI refractoriness in a distinct subgroup of PTC.

1. Introduction

Papillary thyroid carcinoma (PTC) is the most prevalent endocrine malignancy both worldwide [1,2] and in the Middle East [3,4]. Under standard treatment (surgical intervention followed by adjuvant radioactive iodine [131I; RAI] to destroy residual disease, prevent tumor growth and/or recurrence), PTC has excellent prognosis [5,6]. Nevertheless, 5–20% of patients become refractory to RAI therapy, seemingly no longer responding or showing resistance to RAI therapy [7,8,9]. RAI refractoriness in PTC has been associated with poorer clinical outcome, increased mortality [10,11] and limited treatment options that are typically accompanied by considerable toxicity [12,13]. The factors contributing to RAI refractoriness are yet to be fully elucidated. Hence, the clinical and genomic profiling of RAI-refractory PTCs are essential to guide and develop more effective diagnostic/prognostic biomarkers and alternative therapies.
Mutational signatures can reflect a predisposition towards specific mutations due to active endogenous and/or exogenous biological processes operating on the human cancer genome [14,15]. Accumulated mutations may be resultant of multiple mutagenic processes over the lifespan of a tumor [16]. Mathematical approaches allow individual mutational signatures to be deconstructed from somatic mutation catalogues, whilst approximating the number of mutations with the probability that they were caused by each mutational signature [17,18,19,20,21]. Consequently, mutational processes such as aging and the pathological activation of apolipoprotein B mRNA-editing enzymes (APOBEC) have recently been implicated in the heterogeneous nature of PTC [16,22,23,24]. Importantly, mutational signatures can serve as biomarkers of drug response and disease prognosis, directly impacting patient management and treatment options [25,26,27]. Thus, it is necessary to assess the clinical potential of mutational signatures in RAI-refractory PTC.
This study aims to identify the characteristic clinicopathological and genomic features in RAI-refractory PTCs that can help predict RAI refractory disease early (following primary surgical excision). Whole-exome sequencing (WES) was performed on 158 (66 RAI refractory and 92 avid) adult (>18 years) PTCs and patient-matched germline. RAI refractoriness was defined based on the recently published joint consensus from the American Thyroid Association, the European Association of Nuclear Medicine, the European Thyroid Association, the Society of Nuclear Medicine and Molecular Imaging on Current Diagnostic and Theranostic Approaches, and current literature [6,28,29,30]. All PTC samples were sequenced prior to the classification of RAI refractory disease.
Overall, this study highlights a distinct aggressive clinical nature and select genomic features characteristic of RAI-refractory PTC, whilst proposing the APOBEC SBS13 mutational signature as a novel independent predictor of RAI refractoriness in PTC. The detection of RAI refractory disease following standard surgical removal of PTCs and more suitable therapeutic interventions may potentially improve patient outcomes.

2. Materials and Methods

2.1. Clinical Cohort

One-hundred and fifty-eight (66 RAI refractory and 92 avid) PTC samples, normal whole blood and clinicopathological data were collected from adult (>18 years) patients, from King Faisal Specialist Hospital and Research Center and Prince Sultan Military Medical City. The tumors were not subjected to any therapeutic regimen prior to surgical excision. Upon optimized surgery, all patients received ≥1 dose of adjuvant RAI therapy. Patients were followed up at 3–6-month intervals for neck ultrasound, thyroid function tests, serum thyroglobulin levels and thyroglobulin antibody analyses. Progression-free survival was defined as the time, in months, from the date of the initial surgery to the detection of any tumor recurrence (local, regional, or distant). In case of no recurrence, the date of last follow-up was considered the study endpoint. Clinical and treatment data were collected from medical records and are summarized in Supplementary Table S1 and Figure S1.
This study was conducted in accordance with the Declaration of Helsinki and approved by the Institutional Review Board of the King Faisal Specialist Hospital & Research Center, under the Project Research Advisory Committee (RAC) #2110 031 & #2211 168 on PTC archival clinical samples. Written consent was obtained from all patients included in the study.

2.2. Imaging Analysis

Whole body 131I scans were performed via SPECT (single-photon emission computerized tomography) post-surgery to assess uptake status, and post-131I therapy to assess response to therapy. Local recurrences and distant metastases were confirmed by 131I scan, CT (computed tomography), MRI (magnetic resonance imaging), 18F-FDG PET/CT (fluorodeoxyglucose-positron emission tomography/computed tomography), or bone scan, where relevant.

2.3. Definition of Radioactive Iodine Refractoriness

Based on the recently published joint consensus from the American Thyroid Association, the European Association of Nuclear Medicine, the European Thyroid Association, the Society of Nuclear Medicine and Molecular Imaging on Current Diagnostic and Theranostic Approaches, and current literature [6,28,29,30], PTCs were classified as RAI refractory if any of the following were fulfilled:
  • 131I uptake absent on diagnostic 131I scan of locoregional recurrence or distant metastasis.
  • 131I uptake absent on a 131I scan performed several days after 131I treatment.
  • 131I uptake present in some, but not all tumor foci.
  • Disease progression despite a cumulative 131I activity of ≥600 mCi.
  • Metastatic disease progression despite 131I uptake.
  • Rising serum thyroglobulin levels ≥6 months after 131I treatment.
  • Structural disease progression after 131I treatment.

2.4. Sample Processing

DNA was isolated from formalin-fixed, paraffin-embedded (FFPE) tumor and blood samples using a Gentra DNA isolation kit (Gentra, Minneapolis, MN, USA), following the manufacturer’s recommendations as described previously [31]. DNA was quantified by Qubit (Invitrogen, Waltham, MA, USA).

2.5. Whole-Exome Sequencing

Preceding the classification of RAI refractoriness, whole-exome sequencing (WES) was performed using SureSelectXT Target Enrichment (Agilent, Santa Clara, CA, USA) on Illumina NovaSeq 6000, at median depth 198× (range 84–543) for tumors, and 204× (range 94–665) for germline. The sequencing reads were aligned to the human reference genome hg19 using Burrows–Wheeler Aligner (BWA) v0.7.15 algorithm [32], followed by local realignment and PCR duplicate marking via Picard tools (v1.119). Base-quality recalibration was performed with GATK v3.8.0 [33]. All quality metrics were obtained using GATK and FastQC.

2.6. Mutation Calling

Somatic single nucleotide variants (SNVs) were identified via MuTect v1.1.7 [34], while somatic small insertions and deletions (indels) were identified via VarScan v2.3.9. Variants were annotated using ANNOVAR [35]. SNVs were filtered using the ”KEEP” parameter, and indels were filtered with a somatic p-value filter of ≤0.001. Germline-supporting reads had to be ≥8 with a VAF <5%. Common SNPs with a minor allele frequency of >0.01 in dbSNP, the National Heart, Lung, and Blood Institute exome sequencing project, 1000 Genomes, Exome Aggregation Consortium (ExAC) or in-house exome sequencing data of ~800 normal samples, were excluded. Variants present at segmental duplication regions that had a variant allele frequency (VAF) <5% or a regional sequencing depth of the variant <8 reads were also removed.

2.7. Gene Copy Number Profiling

FACETS v0.5.13 [36] was used to determine copy number variants (CNVs), purity and regions of loss of heterozygosity (LOH), while the outputs generated by MuTect and VarScan were utilized to calculate mean allelic frequency (MAF). The CNVs were further defined as gains, losses, amplifications and deletions in relation to average ploidy [37]. CNVs with <10 total markers, and <5 heterozygous markers were excluded from the analysis.

2.8. Telomerase Reverse Transcriptase Promoter Mutations

Although the telomerase reverse transcriptase promoter (TERTp) region is not covered in WES, TERTp hotspot mutations were explored using polymerase chain reaction (PCR) and Sanger sequencing, as previously described [38]. All tumors were included, except three cases (Supplementary Figure S1) for which analysis could not be performed due to insufficient DNA availability. The sequencing traces were analyzed using Mutation Surveyor v5.0.1 (Soft Genetics, LLC, State College, PA, USA). NCBI GenBank reference sequences were used for the mutation analysis.

2.9. Driver Events

All nonsilent mutations and CNVs encompassing genes were assessed for pathogenicity using the Cancer Genomics Interpreter via the Thyroid papillary (“THP”) cancer type [39]. Mutations identified as “TIER 1” due to well-documented evidence of promoting oncogenic transformation in cancer or “TIER 2” with strong but still emerging evidence of having roles in cancer, and CNVs categorized as known driver amplifications or deletions in cancer-related genes, were all classified as driver events. All putative driver mutations were manually validated using Integrated Genomics Viewer (IGV) v2.4.10. Variants with >2 mapping quality zero reads were considered false positives and removed from the analysis. Hotspot TERTp mutations were also considered driver mutations due to their known role in PTC tumor progression [40,41,42].

2.10. Mutational Signature Analysis

Mutational signatures were predicted via the deconstructSigs package in R [21], using the published Catalogue of Somatic Mutations in Cancer (COSMIC) 30 single base substitution (SBS) mutational signatures [43] as a reference. The tumors were initially assessed for all 30 SBS mutational signatures and for associations with RAI avidity. Thereafter, only 20 COSMIC mutational signatures were selected for subsequent analyses due to known biological relevance [16,22,23,24] or high levels and/or prevalence in the analyzed PTC cohort (n = 113) (Supplementary Table S2).

2.11. Statistical Analysis

All statistical analyses were executed on IBM SPSS Statistics (v.21). The Mann–Whitney U Test, the Chi-Squared Test or Fisher’s Exact Test were utilized to compare the continuous and categorical variables, where relevant. Spearman’s rank correlation was used to determine associations. Binary logistic regression modelling was used for the RAI refractory predictive analysis via the “Enter” method in SPSS, where the nonsilent tumor mutation burden (TMB; high/low), mutational signatures (present/absent), BRAF (present/absent) and TERTp mutation (present/absent) were taken as categorical variables. High TMB was defined as a TMB > the cohort median of 8, whereas a TMB ≤ cohort median was considered low. For all the statistical tests performed, the analyses were two-tailed, with p < 0.05 considered statistically significant. For the multivariable analyses (unless specified otherwise), the p-values were adjusted for confounding variables, including age, tumor size and tumor purity, as continuous variables. To decrease the false discovery rate (FDR), the p-values were corrected using the Benjamini–Hochberg method adjusted at FDR <10%, where relevant.

3. Results

3.1. Clinicopathological Associations

Potential clinicopathological features associated with RAI refractoriness were investigated using multivariable analysis in the 66 RAI-refractory versus 92 avid PTCs. Although RAI-refractory patients appeared to be slightly older in age (median 44, range 20–83) than avid patients (median 39, range 20–75), this was not statistically significant (p = 0.327; Mann–Whitney U Test) (Table 1). In contrast, RAI-refractory tumors were significantly associated with positive surgical margins and lymph node metastases present at the time of diagnosis (p = 0.016 and 0.012, respectively; Chi-Squared Test). No other associations were observed.

3.2. Genomic Landscape

A total of 1654 somatic nonsilent mutations (median 8, range 0–161) and 162 copy number variants (CNVs) (median 0, range 0–27) were identified across the PTC cohort (Supplementary Tables S3 and S4). In RAI-refractory PTCs, a median of 11 mutations (range 1–161) and a median of 0 CNVs (range 0–10) were identified, whereas RAI-avid PTCs had a median of 7 mutations (range 1–19) and 0 CNVs (range 0–27) (Supplementary Figure S1). Although the CNV burden did not differ (p = 0.797; Mann–Whitney U test), a significantly higher tumor mutational burden (TMB) was observed in RAI-refractory compared to avid tumors (p < 0.0001; Mann–Whitney U Test). Despite mutations significantly correlating with age (rho = 0.494, p < 0.0001; Spearman’s rank correlation), the significantly higher TMB in RAI-refractory tumors was retained upon multivariable analysis (p = 0.011) (Figure 1A).
While CNVs were sparse in both RAI-refractory and avid PTCs, amplifications on chromosome 1q were found exclusively in RAI-refractory tumors (4.5%, 3/66; p = 0.076; Fisher’s Exact Test), and deletions on chromosome 2q32.1 were only found in RAI-avid tumors (4.3%, 4/92; p = 0.146; Fisher’s Exact Test) (Supplementary Table S4).
Putative driver events (mutations and CNVs) were determined in the PTC cohort. A total of 248 driver events were identified in 118 patients (median 1, range 0–60) (Supplementary Figure S2). On average, significantly more driver events were found in RAI-refractory (1.56 ± 1.67) compared to avid tumors (0.93 ± 0.74) (p = 0.010; Mann–Whitney U Test). Driver mutations in BRAF (48.5%, 32/66 RAI refractory versus 53.3%, 49/92 avid; p = 0.553; Chi-Squared Test) and RAS (7.6%, 5/66 RAI refractory versus 10.9%, 10/92 avid; p = 0.487; Chi-Squared Test) were amongst the most recurrent and mutually exclusive in both RAI-refractory and avid tumors. Whereas TERTp mutations were significantly enriched in RAI-refractory tumors (31.8%, 21/66 versus 2.2%, 2/92 avid; p < 0.0001; Fisher’s Exact Test). To test previously reported BRAFV600E and TERTp mutation (in combination and alone) [44] as predictors of RAI refractoriness, the tumors were divided according to genotype (Supplementary Table S5). By BRAFV600E alone, RAI refractoriness was only detected in 29.2% (odds ratio [OR] 0.7, 95% confidence intervals [CI] 0.3–1.4) of tumors, by TERTp mutation alone in 100% (OR 31.7, 95% CI 1.8–567.5), and by concurrent BRAFV600E and TERTp mutations in 85.7% (OR 11.8, 95% CI 2.5–56.1), versus 37.3% (25/67) of tumors negative for both markers (p = 0.326, 0.019 and 0.002, respectively).
In contrast, AKT3 and MDM4 amplifications were the only driver CNV events found in RAI-refractory PTCs (both in 4.5%, 3/66 of cases; p = 0.076; Fisher’s Exact Test), and only a single RAI-avid tumor (PTC-501) exhibited multiple driver CNV events (p = 1.000; Fisher’s Exact Test).

3.3. Frequency of Mutational Signatures

The mutational processes active in RAI-refractory and avid PTCs were investigated using the 30 published single-base substitution (SBS) mutational signatures from the Catalogue of Somatic Mutations in Cancer (COSMIC) database [43]. Due to the relatively low mutational burden of PTC compared to other cancer types [24], mutational signature analysis could only be performed on 81.8% (54/66) of RAI-refractory and 64.1% (59/92) of avid tumors with ≥10 single nucleotide variants (SNVs) (Supplementary Figure S3). Overall, the median number of SNVs in the 113 tumors was 20 (range 10–345), with 24 (range 10–345) in RAI-refractory and 17 (range 10–36) in RAI-avid tumors (Figure 2A).
The number of SNVs contributing to each mutational signature was investigated in the RAI-refractory and avid PTCs. No mutational signatures were found to be enriched in RAI-avid tumors. However, SNVs in RAI-refractory tumors were significantly enriched for age-related mutations (SBS1) and APOBEC activity, via both SBS2 and SBS13 (Benjamini–Hochberg-adjusted p < 0.001, 0.030 and <0.001, respectively; Mann–Whitney U Test) (Figure 1B).
In contrast with a previous study [24], APOBEC mutational signatures were not associated with BRAFV600E mutation. A total of 60.9% (14/23) of tumors were found harboring concurrent APOBEC mutational signature (SBS2 and/or SBS13) and BRAFV600E mutations, compared to 52.5% (47/90) of the APOBEC-negative tumors found to be harboring a BRAFV600E mutation (p = 0.473; Chi-Squared Test). Similarly, APOBEC mutational signatures were not associated with RAS mutation (APOBEC and RAS-positive tumors: 8.7%, 2/23 versus APOBEC-negative but RAS-positive tumors: 20%, 18/90; p = 0.364; Fisher’s Exact Test). In contrast, APOBEC mutational signatures were significantly associated with TERTp mutation, with concurrent APOBEC and TERTp mutations being harbored in 34.8% (8/23), compared to 14.4% (13/90) of APOBEC-negative tumors harboring a TERTp mutation (p = 0.025; Chi-Squared Test). To determine whether APOBEC-related mutagenesis was the biological process responsible for inducing driver mutations in APOBEC-positive tumors, a probabilistic approach [25] was used. From the 45 driver mutations found in the APOBEC-positive tumors, only seven (15.6%) in three samples comprised of C > T or C > G substitutions in a TCN trinucleotide context. Upon analysis, only four of the seven driver mutations showed a high probability of being induced by APOBEC activity. Notably, the driver mutations MED12S745L c.C2234T (97.1% probability; in case THY-380), PPP2R1AS296F c.C230T (90.6% probability; in case THY-380) and STAG2Q773X c.C2317T (88.5% probability; in case PTC-212) were highly likely to have been induced by APOBEC SBS2 activity. Driver mutation FXR1S400X c.C1199G (in case THY-380) was highly likely, with 100% probability, to have been induced by APOBEC SBS13. Therefore, only two cases were likely to have incurred driver mutations as a direct result of aberrant APOBEC activity. Most driver events, including TERTp (C228T and C250T) and BRAFV600E mutations, were most likely consequent of other biological processes or selection, not APOBEC activity.

3.4. Patterns of Mutational Signatures

Using unsupervised clustering, the overall patterns of mutational signatures between RAI-refractory and avid PTCs were investigated. Ten clusters were identified: four main clusters (C1–C4) and six relatively smaller clusters (C5–C10) (Figure 2B). Cluster C1 predominantly exhibited the clock-like SBS5 mutational signature of unknown etiology. Cluster C2 was characterized by the dominant homologous recombination-based DNA repair (HR; SBS3) mutational signature. Cluster 3 mainly contained age-related mutations (SBS1). Cluster 4 presented with relatively high contributions of SBS16 mutations (of unknown etiology). Cluster 5 involved a high proportion of mutations related to aflatoxin exposure (SBS24). Cluster C6 contained tumors with a mix of defective DNA mismatch-repair (MMR) signatures, including SBS6, SBS15, SBS20, SBS21 and/or SBS26. Cluster C7 predominantly presented the defective DNA base excision repair (BER; SBS30) mutational signature. Two very small clusters showed high contributions of SBS8 (C8) and SBS19 (C9) mutations, both of unknown etiology. Cluster 10 was distinctly characterized by the APOBEC mutational signatures SBS2 and SBS13. Furthermore, while clusters C1–C9 were found in both RAI-refractory and avid tumors, C10 was the principal cluster found in RAI-refractory tumors, which was also coupled with concurrent TERTp mutations (in 55.6%) and higher TMBs (in 77.8%).

3.5. Prediction of RAI Refractoriness

Multivariable logistic regression was used to assess if any of the mutational signatures, TMB, BRAFV600E or TERTp mutations could be used to predict RAI refractoriness independently in PTC. Only the APOBEC SBS13 mutational signature (OR 30.4, 95% CI 1.4–647.2) and TERTp (OR 41.3, 95% CI 4.4–391.6) mutations significantly predicted RAI refractoriness in PTC (p = 0.029 and 0.001, respectively) (Figure 3). In contrast, BRAFV600E (OR 0.3, 95% CI 0.1–1.0) and the mutational signature SBS23 of unknown etiology (OR 0.2, 95% CI 0.03–0.9) mutations significantly predicted a lower risk for developing RAI refractory disease in our PTC cohort (p = 0.043 and 0.034, respectively).
To further investigate the clinical utility of the APOBEC SBS13 mutational signature and TERTp mutation as independent biomarkers to predict RAI refractoriness in PTC, the cohort was divided according to genotype (Table 2). RAI refractoriness was predicted in 87.5% (OR 14.5, 95% CI 1.7–124.0) of tumors with SBS13 alone, 85.7% (OR 12.4, 95% CI 2.6–59.6) of tumors with TERTp mutation alone, and 100% (OR 30.8, 95% CI 1.7–559.4) of tumors harboring concurrent SBS13 and TERTp mutations, versus 32.5% (27/83) of tumors negative for both markers (p = 0.015, 0.002 and 0.021, respectively).

4. Discussion

In a relatively large cohort, the clinical and genomic landscape of RAI-refractory and avid PTCs were systematically profiled, revealing important differences in clinicopathological presentation, TMB, TERTp mutation and mutational signatures.
Whereby, the clinicopathological presentation of Saudi RAI-refractory PTC patients was noticeably more aggressive than avid patients, with multivariable analysis showing significant involvement of surgical margins and presence of lymph node metastases at the time of diagnosis in RAI-refractory PTCs. This was in accordance with RAI-refractory PTC patients of African ancestry [28]. In contrast, Shobab et al. [45] did not find any differences between RAI-refractory and avid patients, possibly be due to their smaller sample size and/or differences in population ethnicity, with >50% of their cohort being of “white” ethnicity.
Despite PTC having a widely accepted lower TMB than most cancers, the average TMB of the PTC cohort was still lower (0.25 nonsynonymous mutations/Mb) compared to TCGA-Thyroid (The Cancer Genome Atlas; 0.41 nonsynonymous mutations/Mb) and previous reporting [46,47]. However, as mutations correlate with age in PTC [24], this can partially be explained by a younger median age at diagnosis of PTC in the Saudi population [3]. Nevertheless, a significantly higher TMB was observed in RAI-refractory compared to avid tumors, though TMB was not found to be an independent predictor of RAI refractoriness.
Notably, the APOBEC SBS13 mutational signature and previously identified TERTp mutations [48,49] were found to be highly predictive as independent biomarkers of RAI refractoriness in a distinct aggressive subgroup of PTC in our cohort, whereas BRAFV600E and SBS23 mutations were more predictive of RAI avidity rather than refractoriness. Even when only considering BRAFV600E and TERTp mutations (alone or in combination), only TERTp mutations alone were predictive of RAI refractoriness in PTC. In contrast to previous reporting [44], BRAFV600E instead suppressed the predictive power when co-existing with TERTp mutation. However, the status of SBS13 was not considered by the authors. Despite reports of mutated BRAFV600E in RAI-refractory patients, the effects of this mutation are inconsistent. Some studies report repression of the plasma glycoprotein involved with iodide uptake, sodium iodine symporter (NIS), due to BRAF mutations in PTCs [50,51], suggesting BRAF mutation is associated with RAI-refractory PTC. Other studies revealed no significant difference in BRAF mutations between RAI-refractory and RAI-avid PTC patients [52,53]. It is known that the genetic alterations involved in PTC tumorigenesis vary across different ethnicities, indicating that genetic and environmental backgrounds can influence the prevalence and effect of genetic alterations in a population. In our study, BRAF mutation was more predictive of RAI avidity than refractoriness, which might partly be explained by differences in patient ethnicity and environmental backgrounds, but this needs to be validated in a larger study.
Furthermore, in contrast to previous reporting [24], APOBEC mutational signatures were not associated with BRAFV600E mutation in our PTC cohort, but were significantly associated with TERTp mutation. APOBEC mutagenesis, specifically through kataegis (when mutations cluster at rearrangement breakpoints) during telomere crisis can promote tumor growth, but not indefinitely [54]. However, when combined with TERTp, TERT re-activating mutations can sustain cellular lifespan by healing short telomeres [55], to help restore their capacity for DNA proliferation [56]. Whether APOBEC mutations were occurring via kataegis clustering or throughout the genome in cases with concurrent TERTp could not be determined due to the limited scope of this study. Future whole-genome mutational signature analysis and functional studies would be required.
Since only SBS13 and TERTp were found to be positive predictors of RAI refractoriness, we focused on testing their ability alone or in combination to predict RAI refractoriness. Both SBS13 and TERTp mutations alone highly predicted RAI refractoriness in PTC. However, when these two markers co-existed, the likelihood of predicting RAI-refractory PTC was significantly increased.
The SBS13 APOBEC mutational signature produces distinctive C > G and C > A mutations (in a TCN trinucleotide context) generated by the uracil excision of error-prone polymerases or abasic sites by uracilDNA glycosylase, in contrast to the characteristic C > T mutations of SBS2, which are triggered by uracil replication following APOBEC cytidine deamination [43,57].
As previously reported [22,24,57], SBS2 and SBS13 in our PTC cohort were also observed in TMB-high thyroid cancers. RAI-refractory tumors formed a unique cluster enriched with APOBEC mutational signatures, high TMB and TERTp mutations.
As standard RAI therapy may not be appropriate for all PTC patients, the development of further treatment options is imperative. As such, ATR (ataxia telangiectasia and Rad3-related kinase) and immune-checkpoint inhibitors could potentially be alternatives for the majority of our RAI-refractory tumors, due to the collective high TMB and/or involvement of aberrant APOBEC activity, via synthetic lethality [25,26,58]. It would be relevant for such patients to be included in clinical trials investigating the benefit of these therapies in RAI-refractory PTC.
Due to the uniqueness of our study and the controversy surrounding the definitions of RAI refractoriness in PTC, the available scientific literature for direct comparison was limited. Still, despite the need of further validation, our findings may help identify patients more likely to be refractory to standard RAI therapy based on their mutational signatures obtained from WES data of routinely excised PTC tumors.

5. Conclusions

This study highlights distinct clinical and genomic features in RAI-refractory PTC, proposing the APOBEC SBS13 mutational signature as a novel independent predictor of RAI refractoriness in a more aggressive subgroup of PTC. Where the detection of RAI refractory disease following the standard surgical removal of PTCs, and more suitable therapeutic interventions may potentially improve patient outcomes.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cancers14061584/s1, Figure S1: Genomic landscape of RAI (A) refractory (n = 66) and (B) avid (n = 92) PTC cohort (total n = 158); Figure S2: Driver events in RAI (A) refractory (n = 50) and (B) avid (n = 66) PTCs; Figure S3: Twenty selected mutational signatures in individual RAI (A) refractory (n = 54) and (B) avid tumors (n = 59); Table S1: Clinical and treatment follow-up data for radioactive iodine (RAI) refractoriness classification of PTC cohort; Table S2: Justification for the 20 selected COSMIC mutational signatures for all subsequent analyses. Overview of distribution and prevalence of 30 COSMIC mutational signatures in the PTC cohort (n = 113), with outcomes for statistical analysis and justifications; Table S3: Somatic mutations in the PTC cohort; Table S4: Copy number variations predicted by Facets algorithm; Table S5: Relationship between BRAFV600E alone, TERT promoter (TERTp) mutation alone, their co-occurrence, and radioactive iodine (RAI) refractoriness in PTC.

Author Contributions

Conceptualization, A.K.S. and K.S.A.-K.; Data curation, S.S., T.M., S.A. and Z.Q.; Formal analysis, S.S., T.M., S.A., Z.Q., S.K.P. and K.S.S.; Funding acquisition, K.S.A.-K.; Investigation, S.S., T.M., S.A., Z.Q., S.K.P., R.B., S.S.A.-S. and M.A.; Methodology, S.S., T.M. and A.K.S.; Project administration, A.K.S. and K.S.A.-K.; Resources, S.S.A.-S. and M.A.; Software, S.S., T.M., S.A., Z.Q., R.B. and K.S.S.; Supervision, A.K.S. and K.S.A.-K.; Validation, S.S., T.M., S.A., Z.Q., R.B. and K.S.S.; Visualization, S.S. and S.A.; Writing—original draft, S.S.; Writing—review and editing, S.S., S.K.P., R.B. and K.S.A.-K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Saudi Government for the Thyroid Cancer Genome Project as part of the International Cancer Genome Consortium—Accelerate Research in Genomic Oncology (ICGC-ARGO) under the Project RAC #2110 031.

Institutional Review Board Statement

This study was conducted in accordance with the Declaration of Helsinki, and approved by the Institutional Review Board of the King Faisal Specialist Hospital & Research Center, under Project RAC #2110 031 & #2211 168 on human PTC archival clinical samples.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study. All clinical data were de-identified.

Data Availability Statement

The whole-exome sequencing data of this study are part of the International Cancer Genome Consortium (ICGC), which have been deposited to the European Genome-phenome Archive (EGA) under the accession number: EGAS00001001788.

Acknowledgments

The authors would like to acknowledge and thank Allianah Benito, Maria Angelita Sabido and Nabil Siraj for their technical assistance.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Wiltshire, J.J.; Drake, T.M.; Uttley, L.; Balasubramanian, S.P. Systematic review of trends in the incidence rates of thyroid cancer. Thyroid 2016, 26, 1541–1552. [Google Scholar] [CrossRef] [PubMed]
  2. Siegel, R.L.; Miller, K.D.; Jemal, A. Cancer statistics, 2020. CA Cancer J. Clin. 2020, 70, 7–30. [Google Scholar] [CrossRef] [PubMed]
  3. Saudi Cancer Registry. Cancer Incidence Report. Saudi Arabia. 2015. Available online: https://nhic.gov.sa/eServices/Documents/ESCRfinal6NOV.pdf (accessed on 4 September 2020).
  4. Al-Zaher, N.; Al-Salam, S.; El Teraifi, H. Thyroid carcinoma in the United Arab Emirates: Perspectives and experience of a tertiary care hospital. Hematol. Oncol. Stem Cell Ther. 2008, 1, 14–21. [Google Scholar] [CrossRef] [Green Version]
  5. Grant, C.S. Recurrence of papillary thyroid cancer after optimized surgery. Gland. Surg. 2015, 4, 52. [Google Scholar]
  6. Haugen, B.R.; Alexander, E.K.; Bible, K.C.; Doherty, G.M.; Mandel, S.J.; Nikiforov, Y.E.; Pacini, F.; Randolph, G.W.; Sawka, A.M.; Schlumberger, M.; 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 2016, 26, 1–133. [Google Scholar] [CrossRef] [Green Version]
  7. Worden, F. Treatment strategies for radioactive iodine-refractory differentiated thyroid cancer. Ther. Adv. Med. Oncol. 2014, 6, 267–279. [Google Scholar] [CrossRef] [Green Version]
  8. Xing, M.; Haugen, B.R.; Schlumberger, M. Progress in molecular-based management of differentiated thyroid cancer. Lancet 2013, 381, 1058–1069. [Google Scholar] [CrossRef] [Green Version]
  9. Sciuto, R.; Romano, L.; Rea, S.; Marandino, F.; Sperduti, I.; Maini, C.L. Natural history and clinical outcome of differentiated thyroid carcinoma: A retrospective analysis of 1503 patients treated at a single institution. Ann. Oncol. 2009, 20, 1728–1735. [Google Scholar] [CrossRef]
  10. Nixon, I.J.; Whitcher, M.M.; Palmer, F.L.; Tuttle, R.M.; Shaha, A.R.; Shah, J.P.; Patel, S.G.; Ganly, I. The impact of distant metastases at presentation on prognosis in patients with differentiated carcinoma of the thyroid gland. Thyroid 2012, 22, 884–889. [Google Scholar] [CrossRef] [Green Version]
  11. Durante, C.; Haddy, N.; Baudin, E.; Leboulleux, S.; Hartl, D.; Travagli, J.P.; Caillou, B.; Ricard, M.; Lumbroso, J.D.; De Vathaire, F.; et al. Long-term outcome of 444 patients with distant metastases from papillary and follicular thyroid carcinoma: Benefits and limits of radioiodine therapy. J. Clin. Endocrinol. Metab. 2006, 91, 2892–2899. [Google Scholar] [CrossRef]
  12. Faugeras, L.; Pirson, A.S.; Donckier, J.; Michel, L.; Lemaire, J.; Vandervorst, S.; D’Hondt, L. Refractory thyroid carcinoma: Which systemic treatment to use? Ther. Adv. Med. Oncol. 2018, 10, 1758834017752853. [Google Scholar] [CrossRef]
  13. Aashiq, M.; Silverman, D.A.; Na’ara, S.; Takahashi, H.; Amit, M. Radioiodine-Refractory Thyroid Cancer: Molecular Basis of Redifferentiation Therapies, Management, and Novel Therapies. Cancers 2019, 11, 1382. [Google Scholar] [CrossRef] [Green Version]
  14. Temko, D.; Tomlinson, I.P.M.; Severini, S.; Schuster-Bockler, B.; Graham, T.A. The effects of mutational processes and selection on driver mutations across cancer types. Nat. Commun. 2018, 9, 1857. [Google Scholar] [CrossRef] [Green Version]
  15. Alexandrov, L.B.; Stratton, M.R. Mutational signatures: The patterns of somatic mutations hidden in cancer genomes. Curr. Opin. Genet. Dev. 2014, 24, 52–60. [Google Scholar] [CrossRef] [Green Version]
  16. Alexandrov, L.B.; Kim, J.; Haradhvala, N.J.; Huang, M.N.; Tian Ng, A.W.; Wu, Y.; Boot, A.; Covington, K.R.; Gordenin, D.A.; Bergstrom, E.N.; et al. The repertoire of mutational signatures in human cancer. Nature 2020, 578, 94–101. [Google Scholar] [CrossRef] [Green Version]
  17. Alexandrov, L.B.; Nik-Zainal, S.; Wedge, D.C.; Campbell, P.J.; Stratton, M.R. Deciphering signatures of mutational processes operative in human cancer. Cell Rep. 2013, 3, 246–259. [Google Scholar] [CrossRef] [Green Version]
  18. Fischer, A.; Illingworth, C.J.; Campbell, P.J.; Mustonen, V. EMu: Probabilistic inference of mutational processes and their localization in the cancer genome. Genome Biol. 2013, 14, R39. [Google Scholar] [CrossRef] [Green Version]
  19. Gehring, J.S.; Fischer, B.; Lawrence, M.; Huber, W. SomaticSignatures: Inferring mutational signatures from single-nucleotide variants. Bioinformatics 2015, 31, 3673–3675. [Google Scholar] [CrossRef] [Green Version]
  20. Shiraishi, Y.; Tremmel, G.; Miyano, S.; Stephens, M. A Simple Model-Based Approach to Inferring and Visualizing Cancer Mutation Signatures. PLoS Genet. 2015, 11, e1005657. [Google Scholar] [CrossRef] [Green Version]
  21. Rosenthal, R.; McGranahan, N.; Herrero, J.; Taylor, B.S.; Swanton, C. DeconstructSigs: Delineating mutational processes in single tumors distinguishes DNA repair deficiencies and patterns of carcinoma evolution. Genome Biol. 2016, 17, 31. [Google Scholar] [CrossRef] [Green Version]
  22. Masoodi, T.; Siraj, A.K.; Siraj, S.; Azam, S.; Qadri, Z.; Parvathareddy, S.K.; Al-Sobhi, S.S.; AlDawish, M.; Alkuraya, F.S.; Al-Kuraya, K.S. Evolution and Impact of Subclonal Mutations in Papillary Thyroid Cancer. Am. J. Hum. Genet. 2019, 105, 959–973. [Google Scholar] [CrossRef]
  23. Masoodi, T.; Siraj, A.K.; Siraj, S.; Azam, S.; Qadri, Z.; Albalawy, W.N.; Parvathareddy, S.K.; Al-Sobhi, S.S.; Al-Dayel, F.; Alkuraya, F.S.; et al. Whole-Exome Sequencing of Matched Primary and Metastatic Papillary Thyroid Cancer. Thyroid 2020, 30, 42–56. [Google Scholar] [CrossRef] [PubMed]
  24. Pozdeyev, N.; Gay, L.M.; Sokol, E.S.; Hartmaier, R.; Deaver, K.E.; Davis, S.; French, J.D.; Borre, P.V.; LaBarbera, D.V.; Tan, A.C.; et al. Genetic Analysis of 779 Advanced Differentiated and Anaplastic Thyroid Cancers. Clin. Cancer Res. 2018, 24, 3059–3068. [Google Scholar] [CrossRef] [Green Version]
  25. Brady, S.W.; Gout, A.M.; Zhang, J. Therapeutic and prognostic insights from the analysis of cancer mutational signatures. Trends Genet. 2022, 38, 194–208. [Google Scholar] [CrossRef]
  26. Ma, J.; Setton, J.; Lee, N.Y.; Riaz, N.; Powell, S.N. The therapeutic significance of mutational signatures from DNA repair deficiency in cancer. Nat. Commun. 2018, 9, 3292. [Google Scholar] [CrossRef]
  27. Maura, F.; Degasperi, A.; Nadeu, F.; Leongamornlert, D.; Davies, H.; Moore, L.; Royo, R.; Ziccheddu, B.; Puente, X.S.; Avet-Loiseau, H.; et al. A practical guide for mutational signature analysis in hematological malignancies. Nat. Commun. 2019, 10, 2969. [Google Scholar] [CrossRef] [PubMed]
  28. Hurst, Z.; Liyanarachchi, S.; He, H.; Brock, P.; Sipos, J.; Nabhan, F.; Kebebew, E.; Green, P.; Cote, G.J.; Sherman, S.; et al. Risk Haplotypes Uniquely Associated with Radioiodine-Refractory Thyroid Cancer Patients of High African Ancestry. Thyroid 2019, 29, 530–539. [Google Scholar] [CrossRef]
  29. Tuttle, R.M.; Ahuja, S.; Avram, A.M.; Bernet, V.J.; Bourguet, P.; Daniels, G.H.; Dillehay, G.; Draganescu, C.; Flux, G.; Fuhrer, D.; et al. Controversies, Consensus, and Collaboration in the Use of (131)I Therapy in Differentiated Thyroid Cancer: A Joint Statement from the American Thyroid Association, the European Association of Nuclear Medicine, the Society of Nuclear Medicine and Molecular Imaging, and the European Thyroid Association. Thyroid 2019, 29, 461–470. [Google Scholar] [CrossRef] [Green Version]
  30. Gulec, S.A.; Ahuja, S.; Avram, A.M.; Bernet, V.J.; Bourguet, P.; Draganescu, C.; Elisei, R.; Giovanella, L.; Grant, F.; Greenspan, B.; et al. A Joint Statement from the American Thyroid Association, the European Association of Nuclear Medicine, the European Thyroid Association, the Society of Nuclear Medicine and Molecular Imaging on Current Diagnostic and Theranostic Approaches in the Management of Thyroid Cancer. Thyroid 2021, 31, 1009–1019. [Google Scholar] [CrossRef]
  31. Siraj, A.K.; Masoodi, T.; Bu, R.; Beg, S.; Al-Sobhi, S.S.; Al-Dayel, F.; Al-Dawish, M.; Alkuraya, F.S.; Al-Kuraya, K.S. Genomic Profiling of Thyroid Cancer Reveals a Role for Thyroglobulin in Metastasis. Am. J. Hum. Genet. 2016, 98, 1170–1180. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Li, H.; Durbin, R. Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics 2010, 26, 589–595. [Google Scholar] [CrossRef] [Green Version]
  33. McKenna, A.; Hanna, M.; Banks, E.; Sivachenko, A.; Cibulskis, K.; Kernytsky, A.; Garimella, K.; Altshuler, D.; Gabriel, S.; Daly, M.; et al. The Genome Analysis Toolkit: A MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010, 20, 1297–1303. [Google Scholar] [CrossRef] [Green Version]
  34. Cibulskis, K.; Lawrence, M.S.; Carter, S.L.; Sivachenko, A.; Jaffe, D.; Sougnez, C.; Gabriel, S.; Meyerson, M.; Lander, E.S.; Getz, G. Sensitive detection of somatic point mutations in impure and heterogeneous cancer samples. Nat. Biotechnol. 2013, 31, 213–219. [Google Scholar] [CrossRef]
  35. Wang, K.; Li, M.; Hakonarson, H. ANNOVAR: Functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010, 38, e164. [Google Scholar] [CrossRef]
  36. Shen, R.; Seshan, V.E. FACETS: Allele-specific copy number and clonal heterogeneity analysis tool for high-throughput DNA sequencing. Nucleic Acids Res. 2016, 44, e131. [Google Scholar] [CrossRef] [Green Version]
  37. Jamal-Hanjani, M.; Wilson, G.A.; McGranahan, N.; Birkbak, N.J.; Watkins, T.B.K.; Veeriah, S.; Shafi, S.; Johnson, D.H.; Mitter, R.; Rosenthal, R.; et al. Tracking the Evolution of Non-Small-Cell Lung Cancer. N. Engl. J. Med. 2017, 376, 2109–2121. [Google Scholar] [CrossRef] [Green Version]
  38. Bu, R.; Siraj, A.K.; Divya, S.P.; Kong, Y.; Parvathareddy, S.K.; Al-Rasheed, M.; Al-Obaisi, K.A.S.; Victoria, I.G.; Al-Sobhi, S.S.; Al-Dawish, M.; et al. Telomerase reverse transcriptase mutations are independent predictor of disease-free survival in Middle Eastern papillary thyroid cancer. Int. J. Cancer 2018, 142, 2028–2039. [Google Scholar] [CrossRef] [Green Version]
  39. Tamborero, D.; Rubio-Perez, C.; Deu-Pons, J.; Schroeder, M.P.; Vivancos, A.; Rovira, A.; Tusquets, I.; Albanell, J.; Rodon, J.; Tabernero, J.; et al. Cancer Genome Interpreter annotates the biological and clinical relevance of tumor alterations. Genome Med. 2018, 10, 25. [Google Scholar] [CrossRef]
  40. Liu, C.; Liu, Z.; Chen, T.; Zeng, W.; Guo, Y.; Huang, T. TERT promoter Mutation and Its Association with Clinicopathological Features and Prognosis of Papillary Thyroid Cancer: A Meta-analysis. Sci. Rep. 2016, 6, 36990. [Google Scholar] [CrossRef]
  41. McKelvey, B.A.; Umbricht, C.B.; Zeiger, M.A. Telomerase Reverse Transcriptase (TERT) Regulation in Thyroid Cancer: A Review. Front. Endocrinol. 2020, 11, 485. [Google Scholar] [CrossRef]
  42. Donati, B.; Ciarrocchi, A. Telomerase and Telomeres Biology in Thyroid Cancer. Int. J. Mol. Sci. 2019, 20, 2887. [Google Scholar] [CrossRef] [Green Version]
  43. Helleday, T.; Eshtad, S.; Nik-Zainal, S. Mechanisms underlying mutational signatures in human cancers. Nat. Rev. Genet. 2014, 15, 585–598. [Google Scholar] [CrossRef] [PubMed]
  44. Liu, J.; Liu, R.; Shen, X.; Zhu, G.; Li, B.; Xing, M. The Genetic Duet of BRAF V600E and TERT Promoter Mutations Robustly Predicts Loss of Radioiodine Avidity in Recurrent Papillary Thyroid Cancer. J. Nucl. Med. 2020, 61, 177–182. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Shobab, L.; Gomes-Lima, C.; Zeymo, A.; Feldman, R.; Jonklaas, J.; Wartofsky, L.; Burman, K.D. Clinical, Pathological, and Molecular Profiling of Radioactive Iodine Refractory Differentiated Thyroid Cancer. Thyroid 2019, 29, 1262–1268. [Google Scholar] [CrossRef] [PubMed]
  46. Yoo, S.K.; Song, Y.S.; Lee, E.K.; Hwang, J.; Kim, H.H.; Jung, G.; Kim, Y.A.; Kim, S.J.; Cho, S.W.; Won, J.K.; et al. Integrative analysis of genomic and transcriptomic characteristics associated with progression of aggressive thyroid cancer. Nat. Commun. 2019, 10, 2764. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Cancer Genome Atlas Research Network. Integrated genomic characterization of papillary thyroid carcinoma. Cell 2014, 159, 676–690. [Google Scholar] [CrossRef] [Green Version]
  48. Meng, Z.; Matsuse, M.; Saenko, V.; Yamashita, S.; Ren, P.; Zheng, X.; Jia, Q.; Tan, J.; Li, N.; Zheng, W.; et al. TERT promoter mutation in primary papillary thyroid carcinoma lesions predicts absent or lower (131) i uptake in metastases. IUBMB Life 2019, 71, 1030–1040. [Google Scholar] [CrossRef]
  49. Yang, X.; Li, J.; Li, X.; Liang, Z.; Gao, W.; Liang, J.; Cheng, S.; Lin, Y. TERT Promoter Mutation Predicts Radioiodine-Refractory Character in Distant Metastatic Differentiated Thyroid Cancer. J. Nucl. Med. 2017, 58, 258–265. [Google Scholar] [CrossRef] [Green Version]
  50. Riesco-Eizaguirre, G.; Rodríguez, I.; De la Vieja, A.; Costamagna, E.; Carrasco, N.; Nistal, M.; Santisteban, P. The BRAFV600E oncogene induces transforming growth factor β secretion leading to sodium iodide symporter repression and increased malignancy in thyroid cancer. Cancer Res. 2009, 69, 8317–8325. [Google Scholar] [CrossRef] [Green Version]
  51. Dong, H.; Shen, W.-z.; Yan, Y.-j.; Yi, J.-l.; Zhang, L. Effects of BRAFV600E mutation on Na+/I-symporter expression in papillary thyroid carcinoma. J. Huazhong Univ. Sci. Technol. Med. Sci. 2016, 36, 77–81. [Google Scholar] [CrossRef]
  52. Anekpuritanang, T.; Uataya, M.; Claimon, A.; Laokulrath, N.; Pongsapich, W.; Pithuksurachai, P. The Association Between Radioiodine Refractory in Papillary Thyroid Carcinoma, Sodium/Iodide Symporter Expression, and BRAFV600E Mutation. OncoTargets Ther. 2021, 14, 3959. [Google Scholar] [CrossRef]
  53. Zoghlami, A.; Roussel, F.; Sabourin, J.-C.; Kuhn, J.-M.; Marie, J.-P.; Dehesdin, D.; Choussy, O. BRAF mutation in papillary thyroid carcinoma: Predictive value for long-term prognosis and radioiodine sensitivity. Eur. Ann. Otorhinolaryngol. Head Neck Dis. 2014, 131, 7–13. [Google Scholar] [CrossRef] [Green Version]
  54. Maciejowski, J.; Li, Y.; Bosco, N.; Campbell, P.J.; de Lange, T. Chromothripsis and Kataegis Induced by Telomere Crisis. Cell 2015, 163, 1641–1654. [Google Scholar] [CrossRef] [Green Version]
  55. Chiba, K.; Lorbeer, F.K.; Shain, A.H.; McSwiggen, D.T.; Schruf, E.; Oh, A.; Ryu, J.; Darzacq, X.; Bastian, B.C.; Hockemeyer, D. Mutations in the promoter of the telomerase gene TERT contribute to tumorigenesis by a two-step mechanism. Science 2017, 357, 1416–1420. [Google Scholar] [CrossRef] [Green Version]
  56. von Morgen, P.; Maciejowski, J. The ins and outs of telomere crisis in cancer. Genome Med. 2018, 10, 89. [Google Scholar] [CrossRef]
  57. Roberts, S.A.; Gordenin, D.A. Hypermutation in human cancer genomes: Footprints and mechanisms. Nat. Rev. Cancer 2014, 14, 786–800. [Google Scholar] [CrossRef] [Green Version]
  58. Ishaque, N.; Abba, M.L.; Hauser, C.; Patil, N.; Paramasivam, N.; Huebschmann, D.; Leupold, J.H.; Balasubramanian, G.P.; Kleinheinz, K.; Toprak, U.H.; et al. Whole genome sequencing puts forward hypotheses on metastasis evolution and therapy in colorectal cancer. Nat. Commun. 2018, 9, 4782. [Google Scholar] [CrossRef] [Green Version]
Cancers 14 01584 g001 550
Figure 1. Box plots showing (A) the absolute number of alterations, including mutations (tumor mutation burden; TMB) and copy number variations (CNV) in RAI-refractory and avid PTC tumors, and (B) the number of mutations contributing to selected single-base substitution (S; SBS) mutational signatures, with Benjamini–Hochberg-adjusted p-values (Mann–Whitney Test; alpha level 0.05). For each SBS mutational signature, only tumors with ≥1 mutation contributing to the SBS mutational signature were included for box plot representation. S1 represents age-related mutations; S2 and S13 represent APOBEC-related mutagenesis; S6, S15, S20, S21 and S26 represent defective DNA mismatch-repair (MMR) mutations; S3 represents defective homologous recombination-related DNA repair (HR); S30 represents defective base-nucleotide excision repair (BER); S4 represents tobacco smoking; S7 represents ultraviolet (UV) exposure; S24 represents aflatoxin exposure; S29 represents chewing tobacco; and S5, S8, S12, S16, S19 and S23 are of unknown etiology.
Figure 1. Box plots showing (A) the absolute number of alterations, including mutations (tumor mutation burden; TMB) and copy number variations (CNV) in RAI-refractory and avid PTC tumors, and (B) the number of mutations contributing to selected single-base substitution (S; SBS) mutational signatures, with Benjamini–Hochberg-adjusted p-values (Mann–Whitney Test; alpha level 0.05). For each SBS mutational signature, only tumors with ≥1 mutation contributing to the SBS mutational signature were included for box plot representation. S1 represents age-related mutations; S2 and S13 represent APOBEC-related mutagenesis; S6, S15, S20, S21 and S26 represent defective DNA mismatch-repair (MMR) mutations; S3 represents defective homologous recombination-related DNA repair (HR); S30 represents defective base-nucleotide excision repair (BER); S4 represents tobacco smoking; S7 represents ultraviolet (UV) exposure; S24 represents aflatoxin exposure; S29 represents chewing tobacco; and S5, S8, S12, S16, S19 and S23 are of unknown etiology.
Cancers 14 01584 g001
Cancers 14 01584 g002 550
Figure 2. (A) Overview of absolute number of mutations contributing to each mutational signature per sample. (B) Distribution of each mutational signature per sample. Patients are ordered by hierarchical clustering. An alpha level of 0.05 was used for statistical significance. S1 represents age-related mutations; S2 and S13 represent APOBEC-related mutagenesis; S6, S15, S20, S21 and S26 represent defective DNA mismatch-repair (MMR) mutations; S3 represents defective homologous recombination-related DNA repair (HR); S30 represents defective base-nucleotide excision repair (BER); S4 represents tobacco smoking; S7 represents ultraviolet (UV) exposure; S24 represents aflatoxin exposure; S29 represents chewing tobacco; and S5, S8, S12, S16, S19 and S23 are of unknown etiology.
Figure 2. (A) Overview of absolute number of mutations contributing to each mutational signature per sample. (B) Distribution of each mutational signature per sample. Patients are ordered by hierarchical clustering. An alpha level of 0.05 was used for statistical significance. S1 represents age-related mutations; S2 and S13 represent APOBEC-related mutagenesis; S6, S15, S20, S21 and S26 represent defective DNA mismatch-repair (MMR) mutations; S3 represents defective homologous recombination-related DNA repair (HR); S30 represents defective base-nucleotide excision repair (BER); S4 represents tobacco smoking; S7 represents ultraviolet (UV) exposure; S24 represents aflatoxin exposure; S29 represents chewing tobacco; and S5, S8, S12, S16, S19 and S23 are of unknown etiology.
Cancers 14 01584 g002
Cancers 14 01584 g003 550
Figure 3. Multivariable logistic regression model (when adjusting for age, tumor size and tumor purity as continuous variables), using mutational signatures, BRAFV600E mutation, TERT promoter mutation and TMB to predict RAI refractoriness in PTC. An alpha level of 0.05 was used for statistical significance. S1 represents age-related mutations; S2 and S13 represent APOBEC-related mutagenesis; S6, S15, S20, S21 and S26 represent defective DNA mismatch-repair (MMR) mutations; S3 represents defective homologous recombination-related DNA repair (HR); S30 represents defective base-nucleotide excision repair (BER); S4 represents tobacco smoking; S7 represents ultraviolet (UV) exposure; S24 represents aflatoxin exposure; S29 represents chewing tobacco; and S5, S8, S12, S16, S19 and S23 are of unknown etiology.
Figure 3. Multivariable logistic regression model (when adjusting for age, tumor size and tumor purity as continuous variables), using mutational signatures, BRAFV600E mutation, TERT promoter mutation and TMB to predict RAI refractoriness in PTC. An alpha level of 0.05 was used for statistical significance. S1 represents age-related mutations; S2 and S13 represent APOBEC-related mutagenesis; S6, S15, S20, S21 and S26 represent defective DNA mismatch-repair (MMR) mutations; S3 represents defective homologous recombination-related DNA repair (HR); S30 represents defective base-nucleotide excision repair (BER); S4 represents tobacco smoking; S7 represents ultraviolet (UV) exposure; S24 represents aflatoxin exposure; S29 represents chewing tobacco; and S5, S8, S12, S16, S19 and S23 are of unknown etiology.
Cancers 14 01584 g003
Table
Table 1. Clinicopathological characteristics of radioactive iodine (RAI) refractory (n = 66) and avid (n = 92) PTC cases.
Table 1. Clinicopathological characteristics of radioactive iodine (RAI) refractory (n = 66) and avid (n = 92) PTC cases.
Clinicopathological VariablesRAI Refractory
n (%)		RAI Avid
n (%)		Odds Ratio
(95% CI 1)		Multivariable p-Value
Age 0.327
Median44391.0 (0.9–1.0)
Range20–8320–75
Sex 0.483
Male18 (27.3)21 (22.8)0.7 (0.2–2.1)
Female48 (72.7)71 (77.2)
Histopathological subtype
Classical variant33 (55.9)45 (57.7)1
Follicular variant12 (20.3)24 (30.8)0.9 (0.2–4.6)0.918
Tall cell variant14 (23.7)9 (11.5)0.9 (0.2–3.9)0.932
Extrathyroidal extension 0.155
Absent26 (39.4)63 (68.5)1
Present40 (60.6)29 (31.5)4.5 (0.6–35.7)
Lymphovascular invasion 0.239
Absent25 (50.0)47 (55.3)1
Present25 (50.0)38 (44.7)0.5 (0.2–1.6)
Surgical margin involvement 0.016
Absent20 (34.5)66 (73.3)1
Present38 (65.5)24 (26.7)4.8 (1.3–17.0)
Tumour focality 0.577
Unifocal39 (59.1)41 (44.6)
Multifocal27 (40.9)51 (55.4)1.4 (0.5–4.1)
Stage
I34 (51.5)65 (70.7)1
II7 (10.6)9 (9.8)0.3 (0.01–8.3)0.484
III10 (15.2)13 (14.1)2.2 (0.3–15.0)0.413
IV15 (22.7)5 (5.4)3.3 (0.1–72.2)0.453
pT
T18 (12.1)30 (32.6)1
T212 (18.2)23 (25.0)1.2 (0.2–6.3)0.864
T335 (53.0)36 (39.1)0.7 (0.1–6.2)0.771
T411 (16.7)3 (3.3)0.5 (0.01–24.3)0.745
pN 0.012
N024 (37.5)58 (69.9)1
N140 (62.5)25 (30.1)4.4 (1.4–14.3)
pM
M054 (81.8)91 (98.9)
M112 (18.2)1 (1.1)
Types of surgery 0.059
Hemi-thyroidectomy6 (9.1)7 (7.6)1
Total thyroidectomy60 (90.9)85 (92.4)0.1 (0.01–1.1)
1 CI: confidence interval.
Table
Table 2. Relationship between the APOBEC SBS13 mutational signature alone, TERT promoter (TERTp) mutation alone, their co-occurrence, and radioactive iodine (RAI)-refractoriness in PTC.
Table 2. Relationship between the APOBEC SBS13 mutational signature alone, TERT promoter (TERTp) mutation alone, their co-occurrence, and radioactive iodine (RAI)-refractoriness in PTC.
GenotypeTotal Cases
n		RAI Refractory
n (%)		RAI Avid
n (%)		p-Value
(Alpha = 0.05)		Odds Ratio
(95% Confidence Intervals)
SBS13 (−) and TERTp (−)8327 (32.5)56 (67.5)
SBS13 (+) and TERTp (−)87 (87.5)1 (12.5)0.01514.5 (1.7–124.0)
SBS13 (−) and TERTp (+)1412 (85.7)2 (14.3)0.00212.4 (2.6–59.6)
SBS13 (+) and TERTp (+)77 (100.0)0 (0.0)0.02130.8 (1.7–559.4)
SBS13: single base substitution (SBS) mutational signature 13 mutations; TERTp: TERT promoter mutations; (−): negative for/absent; (+): positive for/present.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Siraj, S.; Masoodi, T.; Siraj, A.K.; Azam, S.; Qadri, Z.; Parvathareddy, S.K.; Bu, R.; Siddiqui, K.S.; Al-Sobhi, S.S.; AlDawish, M.; et al. APOBEC SBS13 Mutational Signature—A Novel Predictor of Radioactive Iodine Refractory Papillary Thyroid Carcinoma. Cancers 2022, 14, 1584. https://doi.org/10.3390/cancers14061584

AMA Style

Siraj S, Masoodi T, Siraj AK, Azam S, Qadri Z, Parvathareddy SK, Bu R, Siddiqui KS, Al-Sobhi SS, AlDawish M, et al. APOBEC SBS13 Mutational Signature—A Novel Predictor of Radioactive Iodine Refractory Papillary Thyroid Carcinoma. Cancers. 2022; 14(6):1584. https://doi.org/10.3390/cancers14061584

Chicago/Turabian Style

Siraj, Sarah, Tariq Masoodi, Abdul K. Siraj, Saud Azam, Zeeshan Qadri, Sandeep K. Parvathareddy, Rong Bu, Khawar S. Siddiqui, Saif S. Al-Sobhi, Mohammed AlDawish, and et al. 2022. "APOBEC SBS13 Mutational Signature—A Novel Predictor of Radioactive Iodine Refractory Papillary Thyroid Carcinoma" Cancers 14, no. 6: 1584. https://doi.org/10.3390/cancers14061584

APA Style

Siraj, S., Masoodi, T., Siraj, A. K., Azam, S., Qadri, Z., Parvathareddy, S. K., Bu, R., Siddiqui, K. S., Al-Sobhi, S. S., AlDawish, M., & Al-Kuraya, K. S. (2022). APOBEC SBS13 Mutational Signature—A Novel Predictor of Radioactive Iodine Refractory Papillary Thyroid Carcinoma. Cancers, 14(6), 1584. https://doi.org/10.3390/cancers14061584

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