Skip to main content

ORIGINAL RESEARCH article

Front. Endocrinol., 04 January 2023
Sec. Thyroid Endocrinology
This article is part of the Research Topic Recent Advances in Papillary Thyroid Carcinoma: Diagnosis and Predictive Factors View all 22 articles

Whole-exome sequencing and bioinformatic analyses revealed differences in gene mutation profiles in papillary thyroid cancer patients with and without benign thyroid goitre background

Zing Hong EngZing Hong Eng1Mardiaty Iryani Abdullah,Mardiaty Iryani Abdullah1,2Khoon Leong NgKhoon Leong Ng3Azlina Abdul AzizAzlina Abdul Aziz1Nurul Hannis Arba&#x;ieNurul Hannis Arba’ie3Nurullainy Mat RashidNurullainy Mat Rashid1Sarni Mat Junit*Sarni Mat Junit1*
  • 1Department of Molecular Medicine, Faculty of Medicine, Universiti Malaya, Kuala Lumpur, Malaysia
  • 2Department of Biomedical Science, Kulliyyah of Allied Health Sciences, International Islamic University Malaysia, Kuantan, Pahang, Malaysia
  • 3Department of Surgery, Faculty of Medicine, Universiti Malaya, Kuala Lumpur, Malaysia

Background: Papillary thyroid cancer (PTC) is the most common thyroid malignancy. Concurrent presence of cytomorphological benign thyroid goitre (BTG) and PTC lesion is often detected. Aberrant protein profiles were previously reported in patients with and without BTG cytomorphological background. This study aimed to evaluate gene mutation profiles to further understand the molecular mechanism underlying BTG, PTC without BTG background and PTC with BTG background.

Methods: Patients were grouped according to the histopathological examination results: (i) BTG patients (n = 9), (ii) PTC patients without BTG background (PTCa, n = 8), and (iii) PTC patients with BTG background (PTCb, n = 5). Whole-exome sequencing (WES) was performed on genomic DNA extracted from thyroid tissue specimens. Nonsynonymous and splice-site variants with MAF of ≤ 1% in the 1000 Genomes Project were subjected to principal component analysis (PCA). PTC-specific SNVs were filtered against OncoKB and COSMIC while novel SNVs were screened through dbSNP and COSMIC databases. Functional impacts of the SNVs were predicted using PolyPhen-2 and SIFT. Protein-protein interaction (PPI) enrichment of the tumour-related genes was analysed using Metascape and MCODE algorithm.

Results: PCA plots showed distinctive SNV profiles among the three groups. OncoKB and COSMIC database screening identified 36 tumour-related genes including BRCA2 and FANCD2 in all groups. BRAF and 19 additional genes were found only in PTCa and PTCb. “Pathways in cancer”, “DNA repair” and “Fanconi anaemia pathway” were among the top networks shared by all groups. However, signalling pathways related to tyrosine kinases were the most significantly enriched in PTCa while “Jak-STAT signalling pathway” and “Notch signalling pathway” were the only significantly enriched in PTCb. Ten SNVs were PTC-specific of which two were novel; DCTN1 c.2786C>G (p.Ala929Gly) and TRRAP c.8735G>C (p.Ser2912Thr). Four out of the ten SNVs were unique to PTCa.

Conclusion: Distinctive gene mutation patterns detected in this study corroborated the previous protein profile findings. We hypothesised that the PTCa and PTCb subtypes differed in the underlying molecular mechanisms involving tyrosine kinase, Jak-STAT and Notch signalling pathways. The potential applications of the SNVs in differentiating the benign from the PTC subtypes requires further validation in a larger sample size.

1 Introduction

Abnormal focal growth of thyroid cells resulting in thyroid nodules are very common in the general population. Although the majority of them are diagnosed as benign thyroid goitre (BTG), malignancies occur in 5% to 10% of nodules (1). Thyroid cancer, the most prevalent endocrine malignancy (2), is in the top ten of the most common cancer types in the world’s female population, including Malaysia (3). Papillary thyroid cancer (PTC) is thyroid follicular cell-derived, constituting approximately 80% of all thyroid cancer cases (4). Mutations that trigger oncogenic activation of mitogen-activated protein kinase (MAPK) signalling pathway such as the BRAFV600E mutation are frequently linked to many malignancies including melanoma, colorectal cancer, and PTC (5).

Fine-needle aspiration cytology (FNAC) is a standard pre-operative, minimally invasive procedure to determine thyroid nodule malignancy status. However, up to one-fourth of the cases usually falls into the indeterminate Bethesda categories III and IV (6). Malignancy status is usually confirmed by histopathological examination (HPE) of thyroid tissue samples following partial or total thyroidectomy. HPE is the gold standard for thyroid nodule diagnosis, and yet only 40% of the indeterminate nodules were confirmed to be malignant after HPE (6, 7). Particularly, some of the confirmed PTC through HPE biopsies were found to have a BTG cytomorphological structure. Whether these PTCs are unique PTC subtypes, or the intermediate state of BTG transformation to PTC, remains unknown. Differences in tissue and serum protein profiles were reported in PTC patients with and without BTG background (8) indicative of differences in the underlying mechanisms between the two PTC subtypes and/or differences in their disease progression. This study aimed to evaluate gene mutation profiles to further understand the molecular mechanism underlying BTG, PTC without BTG background and PTC with BTG background. The presence of distinct tumour-related genetic profiles may be able to differentiate the two PTC subtypes in this cohort of patients.

2 Materials and methods

2.1 Subjects

This study was approved by the University of Malaya Medical Centre (UMMC)’s Medical Research Ethics Committee (MREC ID NO: 2019619-7540) in accordance with the ICH GCP guidelines and the Declaration of Helsinki. Written informed consent was obtained from all patients before the study was carried out.

The malignancy status of the patients was assessed through FNA cytology (FNAC) and was further confirmed by histopathological examination (HPE) of tissue specimens following partial or total thyroidectomy. All thyroid tissue specimens were placed in Allprotect tissue reagent (Qiagen, Hilden, Germany) at the time of retrieval and then stored at –80°C until further analysis. The patients were categorised into three groups based on the HPE reports: i) BTG (n = 9), ii) PTC without a BTG cytomorphological background (PTCa) (n = 8) and iii) PTC with a background of BTG (PTCb) (n = 5).

2.2 Genomic DNA extraction from thyroid tissue

Freshly excised thyroid tissue sample from each individual in the respective groups was submerged overnight in Allprotect tissue reagent (Qiagen, Hilden, Germany) at 4°C before storage at –80°C. Genomic DNA (gDNA) was extracted from the tissue samples using Qiagen AllPrep DNA/RNA/Protein Mini Kit (Qiagen, Hilden, Germany) according to the manufacturer’s protocol. The concentration and purity of the extracted gDNA were determined by Invitrogen Qubit dsDNA BR Assay kit (Thermo Fisher Scientific, Massachusetts, USA) on Qubit 2.0 Fluorometer, and by Thermo Scientific NanoDrop™ 2000c Spectrophotometer (Thermo Fisher Scientific, Massachusetts, USA), respectively. The gDNA integrity test was then performed using 1% agarose gel electrophoresis.

2.3 Whole-exome sequencing analysis

Genomic DNA (gDNA) samples of the patients; BTG (n = 9), PTCa (n = 8), and PTCb (n = 5) were sent to BGI Biotechnology Company (Shenzhen, China) for whole-exome sequencing (WES) analysis. The qualified gDNA was randomly fragmented into fragments with a range of 150-200 bp using the Adaptive Focused Acoustics® (AFA®) technology of Covaris Ltd. The adapter-ligated templates were purified using the AgencourtAMPure Solid Phase Reversible Immobilisation (SPRI) beads and fragments with an insert size of about 200 bp were excised. Extracted adapter-ligated templates were then amplified by ligation-mediated polymerase chain reaction (LM-PCR), purified, and hybridised to the SureSelect Biotinylated RNA baits for enrichment and measured using the Agilent 2100 Bioanalyzer. The captured products were then circularised before rolling circle amplification (RCA) was used to produce DNA Nanoballs (DNBs). The captured libraries were then sequenced using the BGISEQ-500 sequencing platform and processed by BGISEQ base-calling Software with default parameters. The filtered WES data in FASTQ format were aligned to Reference Genome of human genome built 37 (http://hgdownload.cse.ucsc.edu/goldenPath/hg19/bigZips/) using Burrows-Wheeler Aligner (BWA) software. Genomic variations were detected by HaplotypeCaller of Genome Analysis Toolkit (GATK) (v3.6). The hard-filtering method was applied to obtain high confident variant calls. All the SNVs were then annotated using the SnpEff program (http://snpeff.sourceforge.net/). The WES results were validated by BRAFV600E mutation screening using PCR-direct DNA sequencing (9). Primers targeting the BRAF mutation were designed using Primer3 (https://primer3.ut.ee/). The primer sequence is as follows; Forward: 5’-CTCTTCATAATGCTTGCTCTGATAG-3’; Reverse: 5’-CCTCAATTCTTACCATCCAC-3’.

2.4 Variant filtering and principal component analysis

The variants with the following criteria were retained for further analyses: (i) minor allele frequency (MAF) ≤ 0.01 in the 1000 Genomes Project control database; and (ii) nonsynonymous and splice-site variants. Principal component analysis (PCA) was carried out on the retained SNVs using Python 3.6 software from a single computer workstation equipped with 3.1 GHz Dual-Core Intel Core i5 and 8 GB RAM. All the software and Python libraries used in the study work were from open sources. The following four variables of categorical data were used to represent an SNV molecular component chosen for the PCA: sample ID (e.g. patients 1, 2 etc), mutant genes (e.g. AATF), specific mutation (e.g. NM_012138:c.G872A:p.G291D) and genotype status (heterozygous/homozygous) as well as three classes of the disease type (BTG, PTCa and PTCb). One hot encoding was then run to convert the categorical variables into the numerical interpretation. A new dimension data of 22,427 rows and 22,850 columns were produced. All data were standardised in such a manner that it has mean as 0 and standard deviation as 1. During model training, the PCA algorithm was carried out to reduce high dimensional data into a minimum number of dimensions with 90% of the retained variance.

2.5 Database search of single nucleotide variants in tumour-related genes

The variants were filtered against two cancer-related gene panels from Oncology Knowledge Base (OncoKB) and Catalogue Of Somatic Mutation In Cancer (COSMIC) Cancer Gene Census databases (10, 11). The filtered variants were then screened through the Single Nucleotide Polymorphisms Database (dbSNP, http://www.ncbi.nlm.nih.gov/snp) and COSMIC database for identification of those that were novel. To determine PTC-specific SNVs, those that were identified in more than one PTC patients were further subjected to the Allele Frequency Aggregator (ALFA) project and the genome Aggregation Database (gnomAD) through Ensembl database (https://asia.ensembl.org/index.html) of which SNVs with MAF < 0.01 were retained. The PTC-specific SNVs were then analysed for their potential functional impact using Polymorphism Phenotype 2 (PolyPhen-2, http://genetics.bwh.harvard.edu/pph2/) and Sorting Intolerant from Tolerant (SIFT, https://sift.bii.a-star.edu.sg/). The variants predicted to be deleterious were further validated using Sanger sequencing method.

2.6 Protein-protein interaction network enrichment analysis

PPI enrichment analysis was carried out through Metascape tool (http://metascape.org) with the following databases: STRING, BioGrid, OmniPath, InWeb_IM9 using default parameters. Tumour-related genes that were found to be mutated in each disease group were used as the input gene sets. STRING (physical score > 0.132) and BioGrid were used to identify the protein-protein physical interactions. The Molecular Complex Detection (MCODE) algorithm was then applied to identify densely connected network components in the PPI analysis results. Pathway and process enrichment analyses were then carried out on the MCODE modules based on the GO Biological Processes and KEGG Pathway biological pathway databases.

3 Results

3.1 Nonsynonymous and splice-site variants identification and PCA analysis

A total of 15703 nonsynonymous and splice-site variants were identified in the three cohorts of patients (Figure 1A). PTCa patients had the highest average variant count (≈ 1156) followed by PTCb (≈ 1029) and BTG (893) patients. Missense mutation was found to be the most common type of variant type among all subjects. PCA plots of the filtered WES data for the three groups are shown in Figure 1B. In general, distinctive patterns of the four variables were observed for the BTG (Figure 1B(i)), PTCa (Figure 1B(ii)) and PTCb (Figure 1B(iii)) when the data were analysed separately. A similar distinctive pattern was observed when the datasets from the three disease groups were co-analysed. The variable profiles for PTCa showed a higher degree of dispersal compared to PTCb in relation to BTG. In addition, BTG, PTCa and PTCb showed some extent of similarities in their patterns represented by the shared area at the centre of the PCA plot [Figure 1B(iv)].

FIGURE 1
www.frontiersin.org

Figure 1 Analysis of nonsynonymous and splice-site single nucleotide variants (SNVs) identified through whole exome sequencing in patients with benign thyroid goitre (BTG), papillary thyroid cancer without BTG background (PTCa) and with BTG background (PTCb). (A) Nonsynonymous SNVs identified in each patient were filtered against 1000 Genome Projects (MAF ≤ 0.01). (B) Principal Component Analysis (PCA) of SNV profiles (gene, mutation, and genotype status) of the filtered nonsynonymous SNVs in (i) BTG, (ii) PTCa in relation to BTG, (iii) PTCb in relation to BTG and (iv) all disease groups.

3.2 Analysis of mutated tumour-related genes in BTG, PTCa and PTCb

The 15703 nonsynonymous and splice-site variants were found in a total of 7875 genes in the three groups of patients. When the 7875 genes were filtered against OncoKB and COSMIC cancer gene census databases, a total of 259 tumour-related genes were identified (Figure 2A). OncoKB was able to identify more tumour-related genes in our cohort of patients (410 genes) compared to the COSMIC database (318 genes). Among the 259 mutated tumour-related genes, 149 were found in the BTG group, followed by 144 and 111 in the PTCa and PTCb groups, respectively (Figure 2B). Thirty-six tumour-related genes including ALK, AR, BRCA2, FANCD2 and CBL were found to be mutated in all the three disease groups. An additional 35 genes including ATR and IRS4 were found only in the BTG and PTCa. Another 18 tumour-related genes such as NCOR2 and PALB2 were found to be mutated in BTG and PTCb only. In addition, genes including BRAF, DCTN1 and RAD51B were among the 20 tumour-related genes that were found to be mutated only in the PTCa and PTCb patients. The list of mutated tumour-related genes in each disease group is presented in Supplementary Table 1. A representative of the Sanger sequencing result of the BRAFV600E is shown in Supplementary Figure 1.

FIGURE 2
www.frontiersin.org

Figure 2 Venn diagrams illustrating findings from the database search of mutated tumour-related genes in this cohort of patients. (A) A total of 259 mutated tumour-related genes were listed in both the OncoKB and COSMIC databases. (B) Among the 259 mutated tumour-related genes, there were 36 genes found to be mutated in all three disease groups.

3.3 Protein-protein interaction network enrichment analysis of the mutated tumour-related genes found in BTG, PTCa and PTCb

The top enriched KEGG ontology terms in the PPI networks of the three groups are shown in Figure 3A. “Pathways in cancer” (hsa05200), “Fanconi anaemia pathway” (hsa03460) and “Homologous recombination” (hsa03440) were significantly enriched in the three groups with “Pathways in cancer” and “Fanconi anaemia pathway” were the most enriched in PTCa while “Homologous recombination” was the most enriched in PTCb. “EGFR tyrosine kinase inhibitor resistance” (hsa01521) was exclusively enriched in PTCa while “Jak-STAT signalling pathway” (hsa04630), “Notch signalling pathway” (hsa04330) and “Endocrine resistance” (hsa01522) were only enriched in PTCb.

FIGURE 3
www.frontiersin.org

Figure 3 Protein-protein interaction (PPI) network analysis showing the top enriched pathways for the three disease groups based on the (A) KEGG database and (B) GO ontology database. Pathways with p-value < 0.01 (LogP < –2) was considered as significant. BTG, Benign thyroid goitre; PTCa, Papillary thyroid cancer without BTG background; PTCb, Papillary thyroid cancer with BTG background.

“DNA repair” (GO:0006281), “cellular response to DNA damage stimulus” (GO:0006974) and “double-strand break repair” (GO:0006302) were the enriched Gene Ontology (GO) terms shared among the three groups with the highest significance observed in PTCb followed by PTCa and BTG (Figure 3B). Three GO networks were only significantly enriched in the malignant groups namely “peptidyl-tyrosine phosphorylation” (GO:0018108), “peptidyl-tyrosine modification” (GO:0018212) and “regulation of kinase activity” (GO:0043549), with higher significance shown in PTCa compared to PTCb. “Mismatch repair” (GO:0006298), “gland development” (GO:0048732) and “morphogenesis of an epithelium” (GO:0002009) were highly enriched in the BTG group while “homeostasis of number of cells” (GO:0048872) and “response to steroid hormone” (GO:0048545) were only enriched in PTCb.

The enriched terms were categorised into three main modules: (a) Cancer-related pathways, (b) DNA damage and repair-related pathways, and (c) Signalling Pathways and the tumour-related genes that were involved in each module are shown in Figure 4. Although there was no shared gene among the three groups in the Cancer-related pathways, three genes were shared between the two malignant groups namely AR, HSP90AB1 and GNAS (Figure 4A). NOTCH2 was the only gene shared between BTG and PTCb. In the DNA damage and repair-related pathways, BRCA2 (FANCD1) and FANCD2 were the only two genes shared by the three groups (Figure 4B). In addition to the two genes, ERCC4, NBN, FANCA, FANCE and FANCG were shared between BTG and PTCa; WRN, POLD1 and PALB2 were shared between BTG and PTCb groups while BRIP1 and RAD51B were identified in PTCa and PTCb only. As shown in Figure 4C, none of the genes in the Signalling pathways were shared between the two malignant groups. BTG had none of the enriched term for this module.

FIGURE 4
www.frontiersin.org

Figure 4 Distribution of genes of the enriched pathways in BTG, PTCa and PTCb. The enriched ontology terms are categorised according to the MCODE algorithm modules: (A) Cancer-related pathways and (B) DNA damage and repair-related pathways and (C) Signalling pathways. BTG, Benign thyroid goitre; PTCa, Papillary thyroid cancer without BTG background; PTCb, Papillary thyroid cancer with BTG background.

3.4 Identification of potential molecular markers for PTC

A total of 108 SNVs were identified in more than one subjects. The complete list of the SNVs is presented in Supplementary Table 2. Comparative analysis against the ALFA project and gnomAD dataset through Ensembl database (https://asia.ensembl.org/index.html) identified 10 PTC-specific SNVs (Table 1). Four of the SNVs, ATR c.7817G>A (p.Arg2606Gln), IRS4 c.605A>G (p.Lys202Arg), PCM1 c.3520A>G (p.Thr1174Ala) and TRRAP c.8735G>C (p.Ser2912Thr) were only detected in PTCa. Two SNVs, DCTN1 c.2786C>G (p.Ala929Gly) and TRRAP c.8735G>C (p.Ser2912Thr) were presumed novel as they could not be found in the dbSNP and COSMIC databases. In silico functional analysis through Polyphen2 predicted that the two novel mutations were “possibly damaging” while SIFT predicted that DCTN1 c.2786C>G (p.Ala929Gly) and TRRAP c.8735G>C (p.Ser2912Thr) was “tolerated” and “deleterious” respectively. Another six previously reported SNVs namely ARID1B c.1181C>G (p.Ala394Gly), BRAF c.1799T>A (p.Val600Glu), PDE4DIP c.4073T>C (p.Ile1358Thr), ATR c.7817G>A (p.Arg2606Gln), IRS4 c.605A>G (p.Lys202Arg), and PCM1 c.3520A>G (p.Thr1174Ala) were also predicted to be deleterious. Although predicted to be of benign variants, TET2 c.2440C>T (p.Arg814Cys) and USP6 c.287G>T (p.Arg96Leu) showed PTC-specificity.

TABLE 1
www.frontiersin.org

Table 1 Ten PTC-specific SNVs and their in silico functional impact predictions.

Allelic frequency comparison of six SNPs: BRAF c.1799T>A (p.Val600Glu), TET2 c.2440C>T (p.Arg814Cys), USP6 c.287G>T (p.Arg96Leu), ATR c.7817G>A (p.Arg2606Gln), IRS4 c.605A>G (p.Lys202Arg) and PCM1 c.3520A>G (p.Thr1174Ala) with various normal populations is presented in Figure 5. The minor allele for the six PTC-specific SNVs were either totally absent (BRAF c.1799T>A (p.Val600Glu) and IRS4 c.605A>G (p.Lys202Arg) or present in less than 1% in the world (ALL), African (AFR), American (AMR), East Asian (EAS), European (EUR) and South Asian (SAS) populations. Population genetic data for the other four SNVs; ARID1B c.1181C>G (p.Ala394Gly), DCTN1 c.2786C>G (p.Ala929Gly), PDE4DIP c.4073T>C (p.Ile1358Thr) and TRRAP c.8735G>C (p.Ser2912Thr) were not available in the database for allelic frequency comparison to be made.

FIGURE 5
www.frontiersin.org

Figure 5 Comparison of allelic frequencies of ten PTC-specific SNPs in BTG, PTCa, PTCb and various populations. The data of different populations were obtained from Ensembl population genetic database (https://asia.ensembl.org/index.html). ALL, all populations; AFR, African; AMR, American; EAS, East Asian; EUR, European; SAS, South Asian. Major allele refers to the wildtype allele; minor allele refers to the mutant allele. BTG, Benign thyroid goitre; PTCa, Papillary thyroid cancer without BTG background; PTCb, Papillary thyroid cancer with BTG background.

4 Discussion

To understand the underlying genetic alterations in PTCa and PTCb, WES data of the two groups were compared to those in BTG. WES is one of the NGS platforms that can detect variations in the protein-coding area of genes (exons), which account for approximately 3.09% of the whole genome (12). Its function and application to detect disease-causing mutations have been previously demonstrated in papillary thyroid cancer cases (1316) and other clinical studies (17, 18).

In this study, nonsynonymous variants with MAF ≤ 1% were prioritised due to their direct impact on the functionality and the available tools to interrogate their pathogenic effects (19, 20). The three groups showed similarities in the type of variants of which missense variants being the most common, followed by splice-site and nonsense variants. However, among the three groups, PTCa had the highest total number of variants, followed by PTCb and BTG. DNA damage is common, and it is usually repaired by various DNA repair machineries. The equilibrium between the occurrence of DNA damage and repair, if skewed, leads to DNA damage accumulation, and might start the cancer onset (21). DNA damage accumulation might have contributed towards the higher number of variants in both the malignant PTCa and PTCb groups as compared to the BTG. The different average number of variants between PTCa and PTCb may indicate differences in the underlying disease progression mechanisms or differences in the phases/states during the transformation from the benign to malignant condition. Since the transformed WES data had superhigh dimensions, PCA was done to reduce the dimensionality to find distinct patterns within the data without losing its original information. Ninety percent of variance were set during the model training to get the best visualisation of the distinct patterns while retaining as much information in the datasets as possible (22). The 10% loss in information is to be expected during the dimensional reduction and other factors including duplication, less informative or the complex nature of the information (23). Distinctive profiles were also observed when the filtered WES data for PTCa and PTCb in relation to BTG, were further visualised through PCA plots. This may also point to differences in the underlying molecular mechanisms. Despite the distinctive profiles, some extent of similar patterns in PCA plots were also observed among the BTG, PTCa and PTCb. The pattern similarities could be interpreted that both PTCa and PTCb originate from BTG or PTCb could be an intermediate state in the BTG to PTCa transformation.

To further understand the similarities and differences in the underlying molecular mechanisms behind PTCa and PTCb in relation to BTG, the gene variants identified in all patients were subjected to comparative analyses against OncoKB and COSMIC cancer gene census databases followed by PPI pathway enrichment analyses. MCODE algorithm was applied to identify the densely connected network components in the PPI analysis. The three groups shared significant enrichment of the Cancer-related pathways and DNA damage repair-related pathways. Although no specific tumour-related gene was shared among the three groups, NOTCH2 was shared between BTG and PTCb while AR, HSP90AB1 and GNAS were shared between BTG and PTCa suggesting their possible roles in the BTG-to-PTCb and BTG-to-PTCa transformation, respectively. These observations also suggest that the mechanism of transformation from BTG to PTCa differs from that of BTG to PTCb. In general, the DNA damage and repair-related pathways were more significantly enriched in both PTCa and PTCb compared to BTG. “Fanconi anaemia pathway” was most enriched in PTCa while “Homologous recombination” and “mismatch repair” were most enriched in PTCb and BTG respectively. Notably, BRCA2 and FANCD2 were shared by three groups indicating that the two genes possibly playing key roles in triggering the transformation of BTG to both PTCa and PTCb. DNA lesions can be caused by either endogenous or exogenous agents, such as reactive oxygen species (ROS) produced during cell metabolism or ionising radiation (13). As such, cells are equipped with different DNA repair systems, including base excision repair (BER), nucleotide excision repair (NER), mismatch repair (MMR), homologous recombination (HR) and non-homologous end-joining (NHEJ) against different types of DNA damages and active throughout different cell cycle stages, maintaining our genome stability. Aberrant DNA repair might drive accumulation of the DNA damages, followed by cancer onset. DNA damage in thyroid disorders is common, due to the reliance on ROS in thyroid hormone biosynthesis (14, 15).

Fanconi anaemia (FA) pathway has been found to be important in DNA repair mechanisms (24, 25), although none of our patients had haematologic disorders. Based on the literatures, upon recognition and binding of FA core complex that comprises seven complementary groups (A, B, C, E, F, L and M) to the DNA lesions, the ubiquitination activated FANCD2-I will then localise to the damaged DNA loci. The FA pathway then regulates the subsequent repair processes including DNA replication, cell-cycle arrest, and DNA damage repair. BRCA2 is a recombination mediator that co-localises with ubiquitinated FANCD2 and facilitates the formation of Rad51 nucleofilaments, linking the FA pathways to homologous recombination (26). BRCA2 mutation has been linked to many cancer types, including breast cancer, ovarian cancer, pancreatic cancer and glioblastoma (2730). Generally, cancer patients with the BRCA2 mutations were reported to have a more aggressive phenotype compared to those with FANCD2 mutation (24). In this study, in addition to BRCA2 and FANCD2, another three genes in the FA pathways, FANCA, FANCE and FANCG, were exclusively shared between BTG and PTCa indicating their importance in the transformation of BTG to PTCa with an aggressive phenotype.

In addition to BRCA2 and FANCD2, PALB2 is another gene node which was found in both BTG and PTCb. PALB2 has been reported to be indispensable to BRCA1 and BRCA2 in DNA repair mechanisms (31, 32). In contrast to FANCD2 monoubiquitination which has a mild impact on HR, BRCA2 activity is essential for DNA double-strand break (DSB) repair (24). PALB2 and BRCAs genes have been well-reported due to their high prevalence in breast cancer where patients with PALB2 gene mutation were found to have shorter survival years (3335). It was also reported that approximately 49.1% of thyroid cancer patients were found to have breast tumours (36). PALB2 codes for a protein that serves as a bridge protein between BRCA1 and BRCA2 proteins to form a complex that initiates HR of the DSB (37, 38). Mutations in PALB2 might lead to its BRCAs binding function loss, which leads to the impaired HR pathway, and the accumulation of DNA DSB lesions. The ATM gene found in PTCb’s PPI network further strengthens our speculation. ATM gene codes for the checkpoint kinase ataxia telangiectasia mutated, which activates the intra-S checkpoint in response to DSBs to arrest DNA replication, the possible loss of this kinase function might lead to further DNA lesions accumulation in the PTCb group. Further assessment on DSBs in PTCb with concomitant presence of breast cancer is thus warranted.

While there were degrees of similarities and differences in DNA damage repair-related pathways gene enrichment, in PTCa and PTCb, the “Signalling pathways” enriched network were distinctive in the two groups. Genomic instability has been suggested to be an early event in cancer development, where the consecutive genetic alterations that affect the normal cell cycle machinery might promote progression from a relatively benign proliferative cells lump to malignant tumour (3941). Oncogenic mutations that lead to abnormal synthesis of receptors or ligands involved in signalling pathways, such as growth factor receptor tyrosine kinases (RTKs) and serine/threonine kinase, can cause hyperactive oncogenic signalling pathways and dysregulate cell cycle machinery (40). Signalling pathways related to tyrosine kinases were the most significantly enriched in PTCa while “Jak-STAT signalling pathway” and “Notch signalling pathway” were the only significantly enriched in PTCb. BTG was found to have none of the enriched term for this module. Oncogenic activation of the MAPK signalling cascade was found to be indispensable to the PTC development (31). The oncogenic activation of MAPK cascade upregulates the expression of the receptor tyrosine kinase and its binding to the ligands, such as growth factors, might lead to the oncogenic cell proliferation. However, the absence of this MCODE network in the PTCb group further revealed the different underlying mechanisms in the PTCa and PTCb disease progression. Jak-STAT and Notch signalling pathways were shown to be affected in the PTCb group. Jak-STAT signalling pathway has been reported to have tumour effects in PTC development (42). An additional study had also reported the potential role of affected STAT3 pathway and MAPK signalling cascade due to prolonged H2O2 insult in the thyroid gland and leading to the PTC progression (43, 44). The role of Notch signalling in cancer progression remains controversial, as it was found to be oncogenic and tumour suppressive in different studies. Increased Notch1 and Notch2 expression in PTC cases were found to be associated with more aggressive phenotypes (45). The upregulated Notch1 signalling pathway was linked to metastasis, increased tumour size, and led to a poor prognosis in PTC patients (45, 46). The affected Notch signalling pathway in the PTCb patients in the present study might result in the poor prognosis of PTCb where their presence of BTG cytomorphological structure is in the background. However, discordant findings were also reported (47, 48).

Taken together, findings from the in-silico PPI network enrichment analysis suggest that erroneous DNA repair mechanism might be the primary cause of thyroid tumour development of BTG. Further oncogenic mutations that affect the signalling pathways then drive the transformation of PTCa and PTCb from BTG. The transformation of BTG to PTCa follows a different mechanism which involves a more aggressive tyrosine kinase-related pathways while the transformation of BTG to PTCb is linked to presumably a less aggressive JAK-STAT/Notch signalling pathways. Further confirmation on this is required. The difference in the signalling pathways will allow a more personalised therapeutic target for PTCa and PTCb. The proposed mechanism of transformation from BTG to PTCa and PTCb is summarised in Figure 6.

FIGURE 6
www.frontiersin.org

Figure 6 The proposed mechanism of BTG transformation to PTCa and PTCb in this cohort of patients.

Through the WES sequencing analysis, there were ten SNVs showing PTC-specific characteristic where also predicted to be functionally deleterious in in-silico functional analysis BRAF c.1799T>A (p.Val600Glu) is one of ten SNVs that were found to be PTC-specific. BRAFV600E mutation, which was detected in approximately 60% of PTC cases, is associated with disease aggressiveness, recurrence and increased mortality (17, 18, 4952). Its high prevalence in PTC patients has led to the suggestion of its potential use as a diagnostic and prognostic genetic biomarker (5355). In the current study, BRAFV600E was confirmed to be PTC-specific. Its relatively low occurrence of 23.1% suggests that BRAFV600E is not a common underlying mutation in our cohort of patients. The prevalence of BRAFV600E mutation in PTC patients in Malaysia and in the rest of the Southeast Asia region remains largely unknown (56). While this mutation had been reported to be associated with intra-axial brain tumour patients in Malaysia (57), our study is the first to report the presence of this mutation in PTC patients in Malaysia. Further screening of the BRAFV600E mutation in a larger sample size will better reflect its prevalence in the Malaysian PTC population. While BRAFV600E was able to differentiate the BTG from the malignant PTC cases, it was unable to differentiate between PTCa and PTCb in this study.

Six SNPs, BRAF c.1799T>A (p.Val600Glu), TET2 c.2440C>T (p.Arg814Cys), USP6 c.287G>T (p.Arg96Leu), ATR c.7817G>A (p.Arg2606Gln), IRS4 c.605A>G (p.Lys202Arg) and PCM1 c.3520A>G (p.Thr1174Ala) were present in less than 1% of the world population as the African (AFR), American (AMR), East Asian (EAS), European (EUR) and South Asian (SAS) populations further suggesting that these were pathological instead of neutral mutations. The allelic frequency of another four SNVs {ARID1B c.1181C>G (p.Ala394Gly), DCTN1 c.2786C>G (p.Ala929Gly), PDE4DIP c.4073T>C (p.Ile1358Thr) and TRRAP c.8735G>C (p.Ser2912Thr)} in the various populations were unknown due to data unavailability. While the ten SNVs were PTC-specific, four of the mutations; ATR, IRS4, PCM1 and TRRAP were present only in PTCa. ATR encodes for serine/threonine-protein kinase ATR, is important for DNA damage sensors, activating DNA damage checkpoint against replication stress. The potential loss of function of ATR due to rs199948706 (p.Arg2606Gln) was suggested to increase the DNA damage, initiating the PTCa progression from BTG. Insulin receptor substrate 4 (IRS4) can induce hyperactivation of PI3K/AKT pathway and promote tumourigenesis, in the absence of insulin or other growth factors (56). IRS4 mRNA is expressed in various human tissues including pituitary, thyroid, ovary, prostate, fibroblasts, however, with a low expression level (58). Upregulated IRS4 was detected in various cancer cell lines, and high phosphatidylinositol triphosphate level was also identified (59, 60). PCM1 codes for pericentriolar material 1, chromosomal aberrations of this gene have been linked to various malignancies. RET/PCM1 translocation had been identified in PTC patient, however its role in tumourigenesis is not fully elucidated. TRRAP mutation which was identified in BRAF wildtype PTC patients (61), plays an important role in the recruitment of histone acetyltransferase (HAT) complexes to the chromatin, regulating transcription and DNA repair. Its FATC domain can bind to Myc, regulating Myc oncogenic activities (62, 63). In this present study, the TRRAP p.S2912T mutation is located within the FAT domain, where the function and protein-protein interaction of this domain remains largely unknown. The ten SNVs may be useful as a diagnostic tool to determine malignancy status of thyroid nodules to complement the existing diagnostic methods pending further validation in a clinical setting and in a larger sample size.

In conclusion, distinctive gene mutation patterns were detected in BTG, PTCa and PTCb which corroborated the previous findings on protein profiles in similar cohort of patients. Based on the current study, it is hypothesised that FANCD1 (BRCA2) and FANCD2 along with other genes that encode for various components of the DNA damage and repair-related pathways, play a key role in the progression of BTG to PTC. Gene mutation patterns did not indicate that PTCb is the intermediate state of transformation of BTG to PTCa. Instead, PTCa and PTCb are subtypes that differ in the underlying molecular mechanisms involving tyrosine kinase-related signalling for the former and Jak-STAT and Notch signalling pathways for the latter. However, in vitro functional analysis would need to be carried out to further validate our speculations. The potential applications of the SNVs, especially the ten PTC-specific, in differentiating the benign from the PTC subtypes requires further validation in a larger sample size.

Data availability statement

The data presented in the study are deposited in the Figshare repository, doi: 10.6084/m9.figshare.21714812.

Ethics statement

The studies involving human participants were reviewed and approved by The University of Malaya Medical Centre (UMMC)’s Medical Research Ethics Committee (MREC ID NO: 2019619-7540) in accordance with the ICH GCP guidelines and the Declaration of Helsinki. The patients/participants provided their written informed consent to participate in this study.

Author contributions

ZE participated in the study design, performed a major part of the experiments, analysed and interpreted the data, drafted the manuscript, prepared figures, and tables. MA performed part of the experiments, and reviewed manuscript drafts. KN performed surgeries, involved in specimen collection process, and contributed the clinical data of the patients and reviewed manuscript drafts. AA analysed the data, and reviewed manuscript drafts. NA collected samples, analysed the data, and reviewed drafts of the paper. NR analysed the data, and reviewed drafts of the paper. SJ conceived and designed the experiments, analysed, and interpreted the data, reviewed drafts, and critically revised the final manuscript. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the FRGS/1/2018/SKK08/UM/02/16 grant from the Ministry of Higher Education Malaysia, and the funder had no role in the study design, data collection and analysis, decision to publish, or manuscript preparation.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2022.1039494/full#supplementary-material

References

1. Patel KN, Yip ÃYL, Lubitz CC, Grubbs EG, Miller BS, Shen W, et al. The American association of endocrine surgeons guidelines for the definitive surgical management of thyroid disease in adults. Ann Surg (2020) 271(3):e21–93. doi: 10.1097/SLA.0000000000003580

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Jin M, Li Z, Sun Y, Zhang M, Chen X, Zhao H, et al. Association analysis between the interaction of RAS family genes mutations and papillary thyroid carcinoma in the han Chinese population. Int J Med Sci (2021) 18(2):441–7. doi: 10.7150/ijms.50026

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Manan AA, Basri H, Kaur N, Rahman SZA, Amir PN, Ali N, et al. Malaysia National cancer registry report (MNCRR) 2012-2016. Malaysia: National Cancer Registry (2019) p. 1–116. Available at: https://www.moh.gov.my/moh/resources/Penerbitan/Laporan/Umum/2012-2016%20(MNCRR)/MNCR_2012-2016_FINAL_(PUBLISHED_2019).pdf.

Google Scholar

4. Kim J-K, Seong CY, Bae IE, Yi JW, Yu HW, Kim S-J, et al. Comparison of immunohistochemistry and direct sequencing methods for identification of the BRAF(V600E) mutation in papillary thyroid carcinoma. Ann Surg Oncol (2018) 25(6):1775–81. doi: 10.1245/s10434-018-6460-3

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Cohn AL, Day B-M, Abhyankar S, McKenna E, Riehl T, Puzanov I. BRAF(V600) mutations in solid tumors, other than metastatic melanoma and papillary thyroid cancer, or multiple myeloma: A screening study. Onco Targets Ther (2017) 10:965–71. doi: 10.2147/OTT.S120440

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Wang H, Mehrad M, Ely KA, Liang J, Solórzano CC, Neblett WW III, et al. Incidence and malignancy rates of indeterminate pediatric thyroid nodules. Cancer Cytopathol (2019) 127(4):231–9. doi: 10.1002/cncy.22104

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Paulson VA, Rudzinski ER, Hawkins DS. Thyroid cancer in the pediatric population. Genes (Basel) (2019) 10(9):723. doi: 10.3390/genes10090723

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Abdullah MI, Lee CC, Junit SM, Ng KL, Hashim OH. Tissue and serum samples of patients with papillary thyroid cancer with and without benign background demonstrate different altered expression of proteins. PeerJ (2016) 4:e2450. doi: 10.7717/peerj.2450

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Lee CC, Harun F, Jalaludin MY, Heh CH, Othman R, Novel MJSA. Homozygous c.1502T>G (p.Val501Gly) mutation in the thyroid peroxidase gene in Malaysian sisters with congenital hypothyroidism and multinodular goiter. Int J Endocrinol (2013) 2013:987186. doi: 10.1155/2013/987186

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Chakravarty D, Gao J, Phillips S, Kundra R, Zhang H, Wang J, et al. OncoKB: A precision oncology knowledge base. JCO Precis Oncol (2017) 1:PO.17.00011. doi: 10.1200/PO.17.00011

CrossRef Full Text | Google Scholar

11. Tate JG, Bamford S, Jubb HC, Sondka Z, Beare DM, Bindal N, et al. COSMIC: The catalogue of somatic mutations in cancer. Nucleic Acids Res (2019) 47(D1):D941–7. doi: 10.1093/nar/gky1015

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Suwinski P, Ong C, Ling MHT, Poh YM, Khan AM, Ong HS. Advancing personalized medicine through the application of whole exome sequencing and big data analytics. Front Genet (2019) 10:49. doi: 10.3389/fgene.2019.00049

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Chakarov S, Petkova R, Russev GC, Zhelev N. DNA Damage and mutation. types of DNA damage. Biodiscovery (2014) 11(1):e8957. doi: 10.7750/BioDiscovery.2014.11.1

CrossRef Full Text | Google Scholar

14. Ameziane El Hassani R, Buffet C, Leboulleux S, Dupuy C. Oxidative stress in thyroid carcinomas: Biological and clinical significance. Endocr Relat Cancer (2019) 26(3):R131–43. doi: 10.1530/ERC-18-0476

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Donmez-Altuntas H, Bayram F, Bitgen N, Ata S, Hamurcu Z, Baskol G. Increased chromosomal and oxidative DNA damage in patients with multinodular goiter and their association with cancer. Int J Endocrinol (2017) 2017:2907281. doi: 10.1155/2017/2907281

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Fang Y, Ma X, Zeng J, Jin Y, Hu Y, Wang J, et al. The profile of genetic mutations in papillary thyroid cancer detected by whole exome sequencing. Cell Physiol Biochem (2018) 50(1):169–78. doi: 10.1159/000493966

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Chong JX, Buckingham KJ, Jhangiani SN, Boehm C, Sobreira N, Smith JD, et al. The genetic basis of mendelian phenotypes: Discoveries, challenges, and opportunities. Am J Hum Genet (2015) 97(2):199–215. doi: 10.1016/j.ajhg.2015.06.009

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Jeste SS, Geschwind DH. Disentangling the heterogeneity of autism spectrum disorder through genetic findings. Nat Rev Neurol (2014) 10(2):74–81. doi: 10.1038/nrneurol.2013.278

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Cooper GM, Shendure J. Needles in stacks of needles: Finding disease-causal variants in a wealth of genomic data. Nat Rev Genet (2011) 12(9):628–40. doi: 10.1038/nrg3046

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Gelfman S, Wang Q, McSweeney KM, Ren Z, La Carpia F, Halvorsen M, et al. Annotating pathogenic non-coding variants in genic regions. Nat Commun (2017) 8(1):236. doi: 10.1038/s41467-017-00141-2

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Loeb KR, Loeb LA. Significance of multiple mutations in cancer. Carcinogenesis (2000) 21(3):379–85. doi: 10.1093/carcin/21.3.379

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Cangelosi R, Goriely A. Component retention in principal component analysis with application to cDNA microarray data. Biol Direct (2007) 2(1):2. doi: 10.1186/1745-6150-2-2

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Björklund M. Be careful with your principal components. Evol (N Y). (2019) 73(10):2151–8. doi: 10.1111/evo.13835

CrossRef Full Text | Google Scholar

24. Moldovan G-L, D’Andrea AD. How the fanconi anemia pathway guards the genome. Annu Rev Genet (2009) 43:223–49. doi: 10.1146/annurev-genet-102108-134222

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Walden H, Deans AJ. The fanconi anemia DNA repair pathway: structural and functional insights into a complex disorder. Annu Rev Biophys (2014) 43:257–78. doi: 10.1146/annurev-biophys-051013-022737

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Nakanishi K, Yang Y-G, Pierce AJ, Taniguchi T, Digweed M, D’Andrea AD, et al. Human fanconi anemia monoubiquitination pathway promotes homologous DNA repair. Proc Natl Acad Sci USA (2005) 102(4):1110–5. doi: 10.1073/pnas.0407796102

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Ramus SJ, Harrington PA, Pye C, DiCioccio RA, Cox MJ, Garlinghouse-Jones K, et al. Contribution of BRCA1 and BRCA2 mutations to inherited ovarian cancer. Hum Mutat (2007) 28(12):1207–15. doi: 10.1002/humu.20599

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Tai YC, Domchek S, Parmigiani G, Chen S. Breast cancer risk among male BRCA1 and BRCA2 mutation carriers. J Natl Cancer Inst (2007) 99(23):1811–4. doi: 10.1093/jnci/djm203

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Pilarski R. The role of BRCA testing in hereditary pancreatic and prostate cancer families. Am Soc Clin Oncol Educ B (2019) 39):79–86. doi: 10.1200/EDBK_238977

CrossRef Full Text | Google Scholar

30. Raufi A, Alsharedi M, Khelfa Y, Tirona M. Bilateral triple-negative invasive breast cancer with a BRCA2 mutation, and glioblastoma: A case report and literature review. J Breast cancer (2017) 20:108–11. doi: 10.4048/jbc.2017.20.1.108

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Zaballos MA, Santisteban P. Key signaling pathways in thyroid cancer. J Endocrinol (2017) 235(2):R43–61. doi: 10.1530/JOE-17-0266

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Wu S, Zhou J, Zhang K, Chen H, Luo M, Lu Y, et al. Molecular mechanisms of PALB2 function and its role in breast cancer management. Front Oncol (2020) 10:301. doi: 10.3389/fonc.2020.00301

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Cybulski C, Kluźniak W, Huzarski T, Wokołorczyk D, Kashyap A, Jakubowska A, et al. Clinical outcomes in women with breast cancer and a PALB2 mutation: A prospective cohort analysis. Lancet Oncol (2015) 16(6):638–44. doi: 10.1016/S1470-2045(15)70142-7

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Li A, Geyer FC, Blecua P, Lee JY, Selenica P, Brown DN, et al. Homologous recombination DNA repair defects in PALB2-associated breast cancers. NPJ Breast Cancer (2019) 5(1):23. doi: 10.1038/s41523-019-0115-9

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Szczerba E, Kamińska K, Mierzwa T, Misiek M, Kowalewski J, Lewandowska MA. BRCA1/2 mutation detection in the tumor tissue from selected polish patients with breast cancer using next generation sequencing. Genes (Basel) (2021) 12(4):519. doi: 10.3390/genes12040519

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Kamihara J, LaDuca H, Dalton E, Speare V, Garber JE, Black MH. Germline mutations in cancer predisposition genes among patients with thyroid cancer. J Clin Oncol (2017) 35(15_suppl):1581. doi: 10.1200/JCO.2017.35.15_suppl.1581

CrossRef Full Text | Google Scholar

37. Zhang F, Ma J, Wu J, Ye L, Cai H, Xia B, et al. PALB2 links BRCA1 and BRCA2 in the DNA-damage response. Curr Biol (2009) 19(6):524–9. doi: 10.1016/j.cub.2009.02.018

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Xia B, Sheng Q, Nakanishi K, Ohashi A, Wu J, Christ N, et al. Control of BRCA2 cellular and clinical functions by a nuclear partner, PALB2. Mol Cell (2006) 22(6):719–29. doi: 10.1016/j.molcel.2006.05.022

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Sever R, Brugge JS. Signal transduction in cancer. Cold Spring Harb Perspect Med (2015) 5(4):a006098. doi: 10.1101/cshperspect.a006098

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Yip HYK, Papa A. Signaling pathways in cancer: Therapeutic targets, combinatorial treatments, and new developments. Cells (2021) 10(3):659. doi: 10.3390/cells10030659

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Campbell PJ, Getz G, Korbel JO, Stuart JM, Jennings JL, Stein LD, et al. Pan-cancer analysis of whole genomes. Nature (2020) 578(7793):82–93. doi: 10.1038/s41586-020-1969-6

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Couto JP, Daly L, Almeida A, Knauf JA, Fagin JA, Sobrinho-Simões M, et al. STAT3 negatively regulates thyroid tumorigenesis. Proc Natl Acad Sci USA (2012) 109(35):E2361–70. doi: 10.1073/pnas.1201232109

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Kim WW, Ha TK, Bae SK. Clinical implications of the BRAF mutation in papillary thyroid carcinoma and chronic lymphocytic thyroiditis. J Otolaryngol - Head Neck Surg (2018) 47(1):1–6. doi: 10.1186/s40463-017-0247-6

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Lee CC, Abdullah MI, Junit SM, Ng KL, Wong SY, Fatimah Ramli NS, et al. Malignant transformation of benign thyroid nodule is caused by prolonged H2O2 insult that interfered with the STAT3 pathway? Int J Clin Exp Med (2016) 9(9):18601–17.

Google Scholar

45. Gallo C, Fragliasso V, Donati B, Torricelli F, Tameni A, Piana S, et al. The bHLH transcription factor DEC1 promotes thyroid cancer aggressiveness by the interplay with NOTCH1. Cell Death Dis (2018) 9(9):871. doi: 10.1038/s41419-018-0933-y

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Piana S, Zanetti E, Bisagni A, Ciarrocchi A, Giordano D, Torricelli F, et al. Expression of NOTCH1 in thyroid cancer is mostly restricted to papillary carcinoma. Endocr Connect (2019) 8(8):1089–96. doi: 10.1530/EC-19-0303

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Xiao X, Ning L, Chen H. Notch1 mediates growth suppression of papillary and follicular thyroid cancer cells by histone deacetylase inhibitors. Mol Cancer Ther (2009) 8(2):350–6. doi: 10.1158/1535-7163.MCT-08-0585

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Yu X-M, Jaskula-Sztul R, Georgen MR, Aburjania Z, Somnay YR, Leverson G, et al. Notch1 signaling regulates the aggressiveness of differentiated thyroid cancer and inhibits SERPINE1 expression. Clin Cancer Res (2016) 22(14):3582–92. doi: 10.1158/1078-0432.CCR-15-1749

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Chang C, Chang Y, Huang H, Yeh K, Liu T, Chang J. Determination of the mutational landscape in Taiwanese patients with papillary thyroid cancer by whole-exome sequencing. Hum Pathol (2018) 78:151–8. doi: 10.1016/j.humpath.2018.04.023

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Lee WK, Lee SG, Yim SH, Kim D, Kim H, Jeong S, et al. Whole exome sequencing identifies a novel hedgehog-interacting protein G516R mutation in locally advanced papillary thyroid cancer. Int J Mol Sci (2018) 91(10):2867. doi: 10.3390/ijms19102867

CrossRef Full Text | Google Scholar

51. Xing M. BRAF V600E mutation and papillary thyroid cancer. JAMA (2013) 310(5):535. doi: 10.1001/jama.2013.8592

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Xing M, Liu R, Liu X, Murugan AK, Zhu G, Zeiger MA, et al. BRAF V600E and TERT promoter mutations cooperatively identify the most aggressive papillary thyroid cancer with highest recurrence. J Clin Oncol (2014) 32(25):2718–26. doi: 10.1200/JCO.2014.55.5094

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Nam JK, Jung CK, Song BJ, Lim DJ, Chae BJ, Lee NS, et al. Is the BRAF(V600E) mutation useful as a predictor of preoperative risk in papillary thyroid cancer? Am J Surg (2012) 203(4):436–41. doi: 10.1016/j.amjsurg.2011.02.013

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Xing M, Clark D, Guan H, Ji M, Dackiw A, Carson KA, et al. BRAF mutation testing of thyroid fine-needle aspiration biopsy specimens for preoperative risk stratification in papillary thyroid cancer. J Clin Oncol (2009) 27(18):2977–82. doi: 10.1200/JCO.2008.20.1426

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Trimboli P, Treglia G, Condorelli E, Romanelli F, Crescenzi A, Bongiovanni M, et al. BRAF-mutated carcinomas among thyroid nodules with prior indeterminate FNA report: A systematic review and meta-analysis. Clin Endocrinol (Oxf) (2016) 84(3):315–20. doi: 10.1111/cen.12806

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Rashid FA, Munkhdelger J, Fukuoka J, Bychkov A. Prevalence of BRAF(V600E) mutation in Asian series of papillary thyroid carcinoma-a contemporary systematic review. Gland Surg (2020) 9(5):1878–900. doi: 10.21037/gs-20-430

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Mohamed Yusoff AA, Abd Radzak SM, Mohd Khair SZN, Abdullah JM. Significance of BRAF(V600E) mutation in intra-axial brain tumor in Malaysian patients: Case series and literature review. Exp Oncol (2021) 43(2):159–67. doi: 10.32471/exp-oncology.2312-8852.vol-43-no-2.16076

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Sesti G, Federici M, Hribal ML, Lauro D, Sbraccia P, Lauro R. Defects of the insulin receptor substrate (IRS) system in human metabolic disorders. FASEB J (2001) 15(12):2099–111. doi: 10.1096/fj.01-0009rev

PubMed Abstract | CrossRef Full Text | Google Scholar

59. Hoxhaj G, Dissanayake K, MacKintosh C. Effect of IRS4 levels on PI 3-kinase signalling. PloS One (2013) 8(9):e73327. doi: 10.1371/journal.pone.0073327

PubMed Abstract | CrossRef Full Text | Google Scholar

60. Li X, Zhong L, Wang Z, Chen H, Liao D, Zhang R, et al. Phosphorylation of IRS4 by CK1γ2 promotes its degradation by CHIP through the ubiquitin/lysosome pathway. Theranostics (2018) 8(13):3643–53. doi: 10.7150/thno.26021

PubMed Abstract | CrossRef Full Text | Google Scholar

61. Demeure MJ, Aziz M, Rosenberg R, Gurley SD, Bussey KJ, Carpten JD. Whole-genome sequencing of an aggressive BRAF wild-type papillary thyroid cancer identified EML4-ALK translocation as a therapeutic target. World J Surg (2014) 38(6):1296–305. doi: 10.1007/s00268-014-2485-3

PubMed Abstract | CrossRef Full Text | Google Scholar

62. Park J, Kunjibettu S, McMahon SB, Cole MD. The ATM-related domain of TRRAP is required for histone acetyltransferase recruitment and myc-dependent oncogenesis. Genes Dev (2001) 15(13):1619–24. doi: 10.1101/gad.900101

PubMed Abstract | CrossRef Full Text | Google Scholar

63. Liu X, Tesfai J, Evrard YA, Dent SYR, Martinez E. C-myc transformation domain recruits the human STAGA complex and requires TRRAP and GCN5 acetylase activity for transcription activation. J Biol Chem (2003) 278(22):20405–12. doi: 10.1074/jbc.M211795200

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: benign thyroid goitre, whole-exome sequencing, thyroid tumourigenesis, single nucleotide variants, papillary thyroid cancer

Citation: Eng ZH, Abdullah MI, Ng KL, Abdul Aziz A, Arba’ie NH, Mat Rashid N and Mat Junit S (2023) Whole-exome sequencing and bioinformatic analyses revealed differences in gene mutation profiles in papillary thyroid cancer patients with and without benign thyroid goitre background. Front. Endocrinol. 13:1039494. doi: 10.3389/fendo.2022.1039494

Received: 08 September 2022; Accepted: 07 December 2022;
Published: 04 January 2023.

Edited by:

Jose Federico Carrillo, National Institute of Cancerology (INCAN), Mexico

Reviewed by:

Mudiana Muhamad, MARA University of Technology, Malaysia
Weibo Xu, Fudan University, China

Copyright © 2023 Eng, Abdullah, Ng, Abdul Aziz, Arba’ie, Mat Rashid and Mat Junit. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Sarni Mat Junit, c2FybmlAdW0uZWR1Lm15

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.