Skip to main content

ORIGINAL RESEARCH article

Front. Endocrinol., 31 May 2023
Sec. Thyroid Endocrinology
This article is part of the Research Topic Molecular Characterization of Thyroid Lesions in the Era of “Next Generation” Techniques: Volume II View all 8 articles

Identifying and analyzing the key genes shared by papillary thyroid carcinoma and Hashimoto’s thyroiditis using bioinformatics methods

Ting-ting LiuTing-ting Liu1De-tao Yin,,*De-tao Yin2,3,4*Nan WangNan Wang5Na LiNa Li1Gang DongGang Dong1Meng-fan PengMeng-fan Peng1
  • 1Department of Ultrasound, The First Affiliated Hospital of Zhengzhou University, Zhengzhou, China
  • 2Department of Thyroid Surgery, The First Affiliated Hospital of Zhengzhou University, Zhengzhou, China
  • 3Engineering Research Center of Multidisciplinary Diagnosis and Treatment of Thyroid Cancer of Henan Province, Zhengzhou, China
  • 4Key Medicine Laboratory of Thyroid Cancer of Henan Province, Zhengzhou, China
  • 5Department of Breast Surgery, The First Affiliated Hospital of Zhengzhou University, Zhengzhou, China

Background: Hashimoto’s thyroiditis (HT) is a chronic autoimmune disease that poses a risk factor for papillary thyroid carcinoma (PTC). The present study aimed to identify the key genes shared by HT and PTC for advancing the current understanding of their shared pathogenesis and molecular mechanisms.

Methods: HT- and PTC-related datasets (GSE138198 and GSE33630, respectively) were retrieved from the Gene Expression Omnibus (GEO) database. Genes significantly related to the PTC phenotype were identified using weighted gene co-expression network analysis (WGCNA). Differentially expressed genes (DEGs) were identified between PTC and healthy samples from GSE33630, and between HT and normal samples from GSE138198. Subsequently, functional enrichment analysis was performed using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). Transcription factors and miRNAs regulating the common genes in PTC and HT were forecasted using the Harmonizome and miRWalk databases, respectively, and drugs targeting these genes were investigated using the Drug-Gene Interaction Database (DGIdb). The key genes in both GSE138198 and GSE33630 were further identified via Receiver Operating Characteristic (ROC) analysis. The expression of key genes was verified in external validation set and clinical samples using quantitative real-time polymerase chain reaction (qRT-PCR) and immunohistochemistry (IHC).

Results: In total, 690 and 1945 DEGs were associated with PTC and HT, respectively; of these, 56 were shared and exhibited excellent predictive accuracy in the GSE138198 and GSE33630 cohorts. Notably, four genes, Alcohol Dehydrogenase 1B (ADH1B), Active BCR-related (ABR), alpha-1 antitrypsin (SERPINA1), and lysophosphatidic acid receptor 5 (LPAR5) were recognized as key genes shared by HT and PTC. Subsequently, EGR1 was identified as a common transcription factor regulating ABR, SERPINA1, and LPAR5 expression. These findings were confirmed using qRT-PCR and immunohistochemical analysis.

Conclusion: Four (ADH1B, ABR, SERPINA1, and LPAR5) out of 56 common genes exhibited diagnostic potential in HT and PTC. Notably, this study, for the first time, defined the close relationship between ABR and HT/PTC progression. Overall, this study provides a basis for understanding the shared pathogenesis and underlying molecular mechanisms of HT and PTC, which might help improve patient diagnosis and prognosis.

1 Introduction

Hashimoto’s thyroiditis (HT), also known as chronic lymphocytic thyroiditis, is an autoimmune disorder characterized by attack and destruction of the thyroid gland by the immune system (1). Although the etiology of HT is known to be complex, including genetic, environmental, and epigenetic factors, the precise associated mechanisms remain unclear (2). Moreover, HT is associated with an increased risk of various malignant tumors, particularly papillary thyroid carcinoma (PTC) (3). Indeed, multiple studies have proposed a relationship between HT and possible malignant transformation involving immunological and endocrinological pathogenic links. Thus, chronic inflammation might function as an important inducer of cellular transformation and tumor progression in the thyroid (4).

Thyroid carcinoma is the most prevalent endocrine malignant neoplasm worldwide, which includes papillary thyroid carcinoma (PTC), medullary thyroid carcinoma (MTC), anaplastic thyroid carcinoma (ATC), and follicular thyroid carcinoma (FTC) (5). The incidence of PTC has increased in recent years, making it the most common type of thyroid cancer, accounting for approximately 70% of all cases (6). Several well-known risk factors contribute to the development and progression of PTC, including exposure to ionizing radiation (7). However, other risk factors, including sex, obesity, diabetes, smoking, alcohol consumption, and genetic factors, have also been described (8). Meanwhile, compared to patients with PTC but without HT, patients with both PTC and HT exhibited slower tumor growth, less lymphatic spread, and less distant metastasis, possibly because of lymphatic infiltration into the tumor site (9). Therefore, researchers have become increasingly interested in characterizing the factors driving, or the mechanisms underlying, the transformation of chronic inflammation into thyroid malignancy (10, 11).

Since 1955, when the association between HT and PTC was first described (12), numerous epidemiological studies have confirmed the high degree of coexistence between HT and PTC, ranging from 20–85% (1316). Indeed, the incidence of PTC in patients with HT is several times higher than that in patients without HT (1316). Further, despite the clear relationship defined between HT and PTC, no such association has been reported between HT and other thyroid cancers, including MTC, ATC, and FTC (17). Paradoxically, in addition to its role in PTC development, HT has also been shown to protect against PTC progression. More specifically, the prognosis and clinicopathological features of patients with both PTC and HT are better than those of patients without HT. Approximately 18.9–23.2% of patients with PTC are reported to have HT, and show a better prognosis than that of patients without HT (14, 18). However, the association between HT and PTC remains controversial, with some arguing that patients with HT undergo more frequent screening, which confounds this association (19). Therefore, identifying the key genes shared by HT and PTC can provide important insights regarding their shared pathogenesis and the molecular mechanisms underlying the development of HT and PTC.

Bioinformatics has been widely employed to analyze diagnostic and therapeutic targets in various diseases (19, 20). In this study, we aimed to screen the key genes shared by HT and PTC, based on a series of bioinformatics analyses using the GSE33630 and GSE138198 datasets (21, 22), along with analyses of the regulatory mechanisms and potential drugs targeting these genes. Subsequently, quantitative real-time polymerase chain reaction (qRT-PCR) and immunohistochemical analysis of clinical PTC samples in HT was performed to confirm the expression trends of the key genes. We believe that this study provides theoretical support for continued investigation into the shared pathogenesis and molecular mechanisms of PTC and HT, and thereby help improves the diagnosis and prognosis of patients with these diseases.

2 Materials and methods

2.1 Dataset collection

In this study, the Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/geo/) was queried using the keywords ‘Hashimoto’s thyroiditis’ and ‘papillary thyroid carcinoma’ to retrieve two datasets for HT and PTC, respectively. In our analysis, we used the transcriptome sequencing data of 45 normal thyroid tissues and 49 PTC tissue samples from the GSE33630 dataset. Further, the transcriptome sequence data of three normal thyroid tissues, 13 HT samples, and eight PTC samples with HT disease background from the GSE138198 dataset were included in our study. RNA sequencing data in thyroid tumors related to The Cancer Genome Atlas (TCGA) cohorts was utilized as an external validation set for the key gene expression (23), including 513 tumor samples and 59 control samples.

2.2 Analysis for the differentially expressed genes (DEGs) in PTC and HT

Based on the data of 45 normal thyroid tissue and 49 PTC tissue samples from GSE33630 dataset, the R package ‘WGCNA’ (24) was used to generate a signed co-expression network for filtering PTC-related genes that were significantly relevant to the PTC phenotype using weighted gene co-expression network analysis (WGCNA). After assessing the presence of outlier samples using unsupervised clustering, we selected an optimal soft-thresholding parameter of 11 to ensure a scale-free network, where scale free R2 = 0.85 and the mean connectivity was close to 0. We then transformed the adjacency matrix into the Topological Overlap Measure (TOM), and conducted hierarchical clustering. Genes with similar expression profiles were classified into the same gene modules using the DynamicTreeCut algorithm with a minimum size of 200. Subsequently, the correlation of modules with traits was calculated and displayed as a heatmap, and the module with the strongest correlation was selected as the key module for further analysis.

Based on |log2FoldChange (FC)| > 1 and False discovery rate (FDR) < 0.05 (where p-values was adjusted with a Benjamini-Hochberg (BH) adjustment for multiple testing through the method of ‘fdr’ in the ‘p.adjust’ function) (25), we authenticated the differentially expressed PTC-related genes in GSE33630 dataset using the R package ‘limma’ package (26). Differentially expressed genes (DEGs) between 13 HT and 3 normal samples in the GSE138198 dataset were also mined using Wilcoxon test with |log2FC| > 1 and p-value < 0.05 as the filtering criteria.

2.3 Functional annotation analysis

The R package ‘clusterProfiler’ (27) was first used to explore the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment terms relevant to the DEGs; GO terms were categorized into cellular components (CC), molecular functions (MF), and biological processes (BP). The significance criterion was adjusted to a p-value < 0.05.

2.4 Recognition of common key genes in HT and PTC

The differentially expressed PTC-related genes and DEGs in HT were intersected and presented in a Venn diagram, where the overlapped genes were considered as common genes and imported into Cytoscape to map the co-expression network with the biological terms using the ClueGO plug-in. The common genes were imported into Cytoscape, and the ClueGO plug-in was used to map the co-expression network between the biological terms and these genes. Furthermore, according to the gene expression patterns, the common genes were included in subsequent Receiver Operating Characteristic (ROC) analyses for identifying their predictive accuracy in GSE33630 (PTC vs Normal) and GSE138198 (HT vs Normal) datasets using the ‘pROC’ package, where the best threshold value was selected according to the Youden index. In addition, genes with area under the curve (AUC) values ≥ 0.95 in the GSE33630 and GSE138198 datasets were identified as common key genes shared by HT and PTC.

2.5 Establishment of the miRNA/TF-HT and PTC related genes-drug network

Key transcription factor (TF)-regulating genes were predicted using the Harmonizome database (http://amp.pharm.mssm.edu/Harmonizome) (CHEA Transcription Factor Targets dataset). The miRNAs targeting the common genes as well as key genes were predicted using the miRWalk database (http://mirwalk.umm.uni-heidelberg.de/) (screening condition: score > 0.9). Drugs targeting the common genes and key genes were predicted using the DGIdb (http://www.dgidb.org) database. The final miRNA/TF-key gene-drug network was mapped using Cytoscape software (28).

2.6 RNA acquisition and real-time quantitative PCR (qRT-PCR)

This study was approved by the Ethics Committee of the First Affiliated Hospital of Zhengzhou University and was conducted according to the principles of the Declaration of Helsinki. Written informed consent was obtained from each patient before sample collection. The tumor tissues and non-tumor adjacent tissues (NAT) collected from nine patients with PTC in HT to assess the expression of key genes. All patients underwent surgery at the First Affiliated Hospital of Zhengzhou University, Henan, China, between July 2021 and March 2022, without receiving any anticancer treatment before surgery. All patients were between 21 and 60 years old and consented to the postoperative follow-up plan (Supplementary Table 1). Tissue specimens were collected within 30 minutes of surgery and immediately frozen in liquid nitrogen. Postoperative monitoring and treatment were continued according to the relevant consensus guidelines. The degree of tumor differentiation was graded according to WHO grading system(29). Total RNA from nine NAT and nine PTC samples was isolated using TRIzol reagent, following the manufacturer’s instructions (Ambion, Austin, Texas, USA). Total RNA was reverse-transcribed into cDNA using a SweScript First Strand cDNA Synthesis Kit (Service bio, Wuhan, China), according to the manufacturer’s protocol. QPCR was subsequently carried out using the 2×Universal Blue SYBR Green qPCR Master Mix (Service bio, Wuhan, China), according to the manufacturer’s instructions. The primer sequences used for PCR are listed in Table 1. The relative expression level was compared to that of the internal reference gene (GAPDH) and calculated using the 2−ΔΔCq method (30).

TABLE 1
www.frontiersin.org

Table 1 The primer sequences used for RT-qPCR.

2.7 Immunohistochemistry (IHC)

We performed immunohistochemical analysis on the tumor tissues and NAT collected from eight patients with PTC in HT. The tissues were collected and fixed in 4% paraformaldehyde, embedded in paraffin, and sectioned into 5 µm thick slices. The sections were placed in EDTA (pH 9.0) for antigen repair, washed with phosphate-buffered saline (PBS, pH 7.4), treated with 3% H2O2, blocked with bovine serum albumin (BSA), and incubated overnight at 4 °C with antibodies specific to ABR (1:200; ab224129, Affinity, Changzhou, Jiangsu, China), ADH1B (1:250; bs-10591R, Bioss, Beijing, China), LPAR5 (1:200; bs-15366R, Bioss, Beijing, China), and SERPINA1 (1:200; bs-0096R, Bioss, Beijing, China). Subsequently, the samples were treated with a secondary goat anti-rabbit IgG antibody (1:200; DAKO) for 50 minutes at 37 °C, and the positive sites were labeled with DAB (diaminobenzidine) developing solution (G1211, Service bio, Wuhan, China). Finally, hematoxylin staining (G1004, Service bio, Wuhan, China) was performed to visualize the nuclei. Images (200× magnification) were captured with a microscope (OLYMPUS), and three different visual fields were analyzed. The calculation of the immune response score (IRS) is typically done manually by visual assessment by the operator or technician. IRS is calculated by multiplying the proportion score (ranging from 0 to 4) by the staining intensity score (ranging from 0 to 3) for each cell type evaluated according to the method proposed by Rem-mele and Stegner(31). The immunoreactivity score (IRS) was determined as IRS = SI × PP (staining intensity× percentage of positive cells) and was divided into four grades as follows: 0, non-staining means no positive staining (0–5%); 1, Light yellow means weak positive; 2, Brownish yellow means medium positive; 3, tan means strong positive. The cell ratio was then graded on a 0–4 scale (positive cell ratio = positive cell count/total cell count) as follows: 0–5%, grade 0; 6–25%, grade 1, 26–50%, grade 2; 51–75%, grade 3; and 75–100%, grade 4. The resulting scores for each cell type are then summed up to obtain the total IRS score, which can range from 0 to 12. When IRS score>3 is an immune response (+).

2.8 Statistical analysis

All analyses were conducted using R programming language. Data from different groups were compared using the Wilcoxon test or Student’s t-test. Unless otherwise specified, a p-value < 0.05 was considered statistically significant.

3 Results

3.1 Differentially expressed PTC-related genes

To identify the differentially expressed PTC-related genes, WGCNA was performed using the GSE33630 dataset. The sample cluster analysis did not exclude outlier samples (Figure 1A). To ensure that the interactions between the genes maximally conformed to a scale-free distribution, the optimal soft threshold (R2 = 0.85) was selected as 11 (Figure 1B). Next, eight modules were developed based on a gene clustering tree and a dynamic tree-cutting algorithm (with a minimum of 200 genes per gene module) (Figures 1C, D). Following correlation analysis between the modules and sample traits (normal or PTC), the blue module with the highest correlation was selected as the key module (Figure 1D). Ultimately, 4235 genes in the key module were considered PTC-related genes (Supplementary Table 2). Based on these, differential analysis between the PTC and normal groups was performed. In total, 690 divergent genes, including 367 with increased expression and 323 with decreased expression in PTC samples, were characterized (Figure 1E; Supplementary Table 3).

FIGURE 1
www.frontiersin.org

Figure 1 (A) Sample clustering and detection of outliers. (B) The scale-free fit index (scale-free R2, y-axis) as a function of the soft-threshold power (x-axis) and the mean connectivity of the network. (C) Gene expression dendrogram. The color annotations provide a simple visual comparison of module assignments (branch cuttings) based on the dynamic tree-cutting method. (D) Trait-module correlation presented in a Module-Trait Relationships chart; the blue module has the highest correlation and represents the key module. (E) A total of 690 genes were differentially expressed; 367 = (red) were upregulated, and 323 (blue) were downregulated.

Functional enrichment analysis was performed to further investigate the functions of the 690 differentially expressed PTC-related genes. As shown in Supplementary Table 4, 286 GO items (229 BP items, 43 CC items, and 14 MF items) and nine KEGG pathways were identified. The top five items in each category were displayed as chordal graphs (Figure 2). These genes have been implicated in synapse-related biological processes, cell-substrate adhesion, regulation of cell adhesion, hormone metabolic processes, thyroid hormone generation, thyroid hormone metabolic processes, and immune-related biological processes. Moreover, these genes were found to be associated with transcriptional dysregulation in cancer, as well as the TGF-β signaling pathway, ECM-receptor interaction, cytokine-cytokine receptor interaction, complement and coagulation cascades, and p53 signaling pathway.

FIGURE 2
www.frontiersin.org

Figure 2 Functional enrichment analysis of differentially expressed genes (DEGs) in papillary thyroid carcinoma (PTC). (A) Chord plot of the top five Gene Ontology (GO) biological processes (BP) enriched pathways. (B) Chord plot of the top five GO cellular components (CC) enriched pathways. (C) Chord plot of the top five GO molecular functions (MF) enriched pathways. (D) Chord plot of the top five Kyoto Encyclopedia of Genes and Genomes (KEGG) enriched pathways.

3.2 DEGs in HT

Differential expression analysis conducted using the GSE138198 dataset to determine DEGs between HT and normal samples, ultimately revealed 1945 DEGs, including 1264 that were upregulated and 681 that were downregulated in HT samples (Figures 3A, B, Supplementary Table 5). Functional enrichment analysis further revealed 834 GO terms (748 BP, 59 CC, and 27 MF), and 52 KEGG pathways (Supplementary Table 6). The top ten items in each classification are displayed in chordal graphs in Figures 3C, D, these genes were primarily involved in biological processes related to the immune system, oxidative stress, and cell adhesion. Additionally, they were implicated in immune-related pathways, cell adhesion molecules, autoimmune thyroid disease, phagosome, Fc gamma R-mediated phagocytosis, the PPAR signaling pathway, and the FoxO signaling pathway.

FIGURE 3
www.frontiersin.org

Figure 3 (A) Volcano plot of differentially expressed genes (DEGs) between Hashimoto’s thyroiditis (HT) samples and normal samples in the GSE138198 dataset through Wilcoxon test; red points represent upregulated DEGs, and green points represent downregulated DEGs. (B) Heatmap of DEGs between HT samples and normal samples in the GSE138198 dataset; the red squares represent the upregulated DEGs and the blue squares represent the downregulated DEGs. (C) Bar plot of the enriched GO terms of HT DEGs. (D) Bar plot of the enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways of HT DEGs.

3.3 Analysis for common genes in PTC and HT

To explore the key genes shared between PTC and HT, we first compared the DEGs related to PTC-and HT, yielding 32 upregulated and 24 downregulated genes shared among the PTC and HT samples (Figures 4A, C; Supplementary Table 7). The functions and pathways of these 56 differential genes were further investigated by constructing a co-expression network using the ClueGO plug-in in Cytoscape software, which further confirmed that the upregulated common genes were associated with negative regulation of T cell mediated immune response to tumor cell (Figure 4B; Supplementary Table 8). Simultaneously, the downregulated common genes were related to negative regulation of phosphatidylcholine biosynthesis (Figure 4D; Supplementary Table 9).

FIGURE 4
www.frontiersin.org

Figure 4 (A) Venn diagram showing 32 common upregulated differentially expressed genes (DEGs) shared by papillary thyroid carcinoma (PTC) and Hashimoto’s thyroiditis (HT). (B) Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment terms of the 32 common upregulated DEGs shared by PTC and HT visualized using the ClueGO plugin in Cytoscape software. (C) Venn diagram showing the 24 common downregulated DGEs shared by PTC and HT. (D) Function enrichment results of 24 common downregulated DGEs shared by PTC and HT.

The prognostic value of these genes was also evaluated in subsequent ROC analysis using the GSE33630 and GSE138198 datasets. It was found that the AUC value of the 56 common genes was greater than 0.8, revealing excellent predictive accuracy in PTC and HT (Supplementary Table 10).

Considering the important regulatory roles of various miRNAs and TFs in PTC progression, we performed miRNA and TF prediction for these common genes to explore their upstream regulators. Using the miRWalk and Harmonizome databases, the miRNAs/TFs-genes networks targeting the upregulated and downregulated common genes were identified (Supplementary Tables 11, 12), respectively. As shown in Figure 5A, the expression of many upregulated genes was regulated by TF SOX2 and hsa-miR-7110-5p. Meanwhile, HNF4A was the common TF binding to various downregulated genes, such as DEPTOR and TGFBR3 (Figure 5B), and ADH1B was connected with many other genes, such as hsa-miR-6879-5p and hsa-miR-330-5p.

FIGURE 5
www.frontiersin.org

Figure 5 (A) The miRNA/TF-key gene-drug network constructed based on 32 common upregulated differentially expressed genes (DEGs) shared by papillary thyroid carcinoma (PTC) and Hashimoto’s thyroiditis (HT). (B) The miRNA/TF-key gene-drug network constructed based on 24 common downregulated DGEs shared by PTC and HT. The red arrows represent the common DEGs, the green diamonds represent miRNAs, the purple circles represent the transcription factors (TFs), and the red rectangles represent the compounds.

Next, potential therapeutic agents or molecular compounds targeting the common genes were identified for PTC and HT treatment using the DGIdB. The results revealed plenty of compounds that could target upregulated AURKB, including the HER inhibitor LAPATINIB, which can target both AURKB and CCND1 (Figure 5A). Additionally, the downregulated AKR1C3 and ADH1B could be targeted by several compounds including DOXORUBICIN and ALCOHOL (Figure 5B).

3.4 Identification of key genes in PTC and HT

Based on the screening criterion of AUC values ≥ 0.95, four genes (ADH1B, ABR, SERPINA1, and LPAR5) in both the GSE33630 and GSE138198 datasets were identified as common key genes in HT and PTC (Supplementary Table 10; Figures 6A–C), and the Youden’s index for each key gene was calculated in Figure S1.

FIGURE 6
www.frontiersin.org

Figure 6 (A) Venn diagram showing the key genes of Hashimoto’s thyroiditis (HT) and papillary thyroid carcinoma (PTC) with area under the curve (AUC) values > 0.95 in the GSE33630 and GSE138198 datasets. (B) Receiver operating characteristic (ROC) curves of the four key genes (ADH1B, ABR, SERPINA1, LPAR5) in the HT and PTC in the GSE33630 dataset. (C) ROC curves of the four key genes (ADH1B, ABR, SERPINA1, and LPAR5) in HT and PTC in the GSE138198 dataset.

A miRNA/TF-key gene-drug network was generated for the four key genes (Figure S2). Notably, most miRNAs targeted ABR, including miR-145-5p, miR-506-3p, miR-34a, and miR-449a. Among the potential TFs, EGR1 was the common target of SERPINA1, ABR, and LPAR5. YY1 and hsa-miR-297 were potential ADH1B targets, while SERPINA1 and ABR were regulated by MYC, AR, and ESR1. Additionally, SERPINA1 and LPAR5 expression may be affected by RUNX1 and SMAD3. With respect to the potential molecular compounds targeting the key genes, the results revealed four compounds that could target ADH1B (acetaldehyde, fomepizole, alcohol, and abacavir). The alpha-1 antitrypsin inhibitor, IGMESINE (C23H29N) may target SERPINA1. Furthermore, two compounds, arachidonoyl glycine and CHEMBL1630084, could target LPAR5 (Figure S2). However, information for ABR was lacking in DGIdb.

3.5 Verification of key gene expression in clinical samples

Next, the expression patterns of the four key genes were assessed in external validated cohorts (including TCGA-PTC cohorts, GSE138198 dataset and clinical samples with PTC in HT). In the TCGA-PTC dataset three key genes, ABR, LPAR5, and SERPINA1, were found upregulated in PTC samples compared with those in normal samples. In contrast, ADH1B expression was downregulated in PTC samples compared to normal samples (Figure 7A). This trend was shown by the expression data from GSE138198 cohort (Figure 7B). Notably, it can be seen that changes in gene expression in PTC samples with HT background were more pronounced than that in HT samples, indicating the potential synergies between HT and PTC progression.

FIGURE 7
www.frontiersin.org

Figure 7 (A) Expression of key genes in PTC samples compared to normal samples in the TGCA-PTC dataset. (B) Expression of key genes in three normal thyroid tissue samples, 13 Hashimoto’s thyroiditis (HT) samples and eight papillary thyroid carcinoma (PTC) samples with HT disease background in the GSE138198 dataset. The asterisks represent the significance of the difference, the more asterisks, the more significant the difference.

To further validate these expression patterns, we performed RT-qPCR on tumor tissues and NAT collected from nine patients with PTC in HT. Consistent with the bioinformatics data analysis results, ABR, LPAR5, and SERPINA1 were significantly upregulated in clinical PTC samples compared to NAT, whereas ADH1B was significantly downregulated (Figure 8). Additionally, we assessed the expression of the key genes in tumor tissues and NAT collected from eight patients with PTC in HT using IHC. In agreement with the bioinformatics data analysis results, the abundance of ABR, LPAR5, and SERPINA1 was significantly upregulated in clinical PTC samples compared to NAT, whereas ADH1B was significantly downregulated (Figure 9).

FIGURE 8
www.frontiersin.org

Figure 8 Results of the public database data analysis confirmed using qRT-PCR. In nine papillary thyroid carcinoma (PTC) samples with Hashimoto’s thyroiditis (HT) disease. (A) ABR was significantly upregulated in cancer tissues (PTC) compared to non-tumor adjacent tissue (NAT). (B) ADH1B was significantly downregulated in PTC compared to NAT. (C) LPAR5 was significantly upregulated in PTC compared to NAT. (D) SERPINA1 was significantly upregulated in PTC compared to NAT. The asterisks represent the significance of the difference, the more asterisks, the more significant the difference.

FIGURE 9
www.frontiersin.org

Figure 9 ABR, ADH1B, LPAR5, and SERPINA1 abundance in papillary thyroid carcinoma (PTC). (A) Immunohistochemical (IHC) analysis of ABR, ADH1B, LPAR5, and SERPINA1 abundance in cancer tissues (PTC) and non-tumor adjacent tissue (NAT) of eight PTC samples with Hashimoto’s thyroiditis (HT) disease. (B) IHC scores of ABR, ADH1B, LPAR5, and SERPINA1 in PTC and NAT. The asterisks represent the significance of the difference, the more asterisks, the more significant the difference.

4 Discussion

We assessed differences in gene expression between PTC, HT, and normal thyroid samples. First, we identified 690 divergent genes between PTC and normal thyroid samples. Functional enrichment analysis revealed that PTC occurrence may be associated with synapse-related biological and hormone metabolic processes, cell-substrate adhesion, cell adhesion regulation, p53 and TGF-β signaling, and immune-related biological processes. Inflammatory mediators contribute to tumor progression via myriad mechanisms (32). Subsequently, 1945 DEGs were detected between HT and normal samples in the GSE138198 dataset. These genes were implicated in immune-related and oxidative stress-related pathways, phagosomes, Fc gamma R-mediated phagocytosis, the PPAR and FoxO signaling pathways, and cell adhesion-related biological processes. Thus, the biological pathways associated with the DEGs of HT and PTC are closely related. Although HT etiology has not yet been fully elucidated, it is associated with the release of cytokines that can damage thyroid follicular cells following helper T lymphocyte activation (33, 34). Inflammatory cytokines, such as tumor necrosis factor-α, IL-6, and IL-1, are associated with COX-2 expression, which is related to various malignancies (34, 35). Moreover, cross-reactivity with lactoperoxidase, that induces chronic inflammation in HT, can promote PTC (36). Chronic inflammation in turn leads to matrix remodeling, resulting in changes in cell adhesion, DNA damage, genetic alterations (such as p53 mutations), and protective antitumor immunosuppression (3739). Therefore, we suggest that chronic inflammation represents a possible mechanism of PTC owing to its effect on the immune microenvironment and tumor suppressor gene mutation induction (32).

There were 56 common genes in PTC and HT, of which the 32 up-regulated genes were related to the negative regulation of T cell-mediated tumor cell immune response, while the 24 down-regulated genes were related to the negative regulation of phosphatidylcholine biosynthesis. HT is an autoimmune disease mainly mediated by T cells (40). Phosphatidylcholine synthesis plays an important role in HT (41), but these pathways have not been reported in PTC. Due to the strong correlation between PTC and HT, these pathways can provide new research directions for future PTC studies.

In addition to ROCAUC analysis of the 56 common genes, we further investigated the transcriptional regulation and drugs targeting these genes to gain a deeper understanding of the progression and potential therapeutic mechanisms of PTC. Increased expression of SOX2 promotes the progression of PTC (42, 43). MiR-330-5p can promote the progression of PTC (44) by upregulating ETV4. These genes were affected by competing endogenous RNA (ceRNA) regulatory mechanisms, which may provide new therapeutic targets for the treatment of PTC. Furthermore, miR-6879-5p, as a key exocrine marker, is associated with the progression of PTC (44). Here, the target gene of miR-330-5p and miR-6879-5p was ADH1B, both of which have been reported in PTC, but their interaction mechanism has not been studied, and further exploration is warranted. Previous bioinformatics analysis (45) showed that HNF4A has a significant effect on the progression of PTC, but the mechanism of its interaction with the target gene SERPINA1 has not been investigated. Our results indicate that HNF4A is a common TF-binding site for many downregulated genes. However, functional analysis of these genes awaits further study. Further investigation of these regulatory mechanisms could elucidate the underlying molecular mechanisms of PTC progression. In our drug-prediction analysis, we focused on two genes, CCND1 and AKR1C3, which have been implicated in PTC pathogenesis. LAPATINIB has been studied in the treatment of PTC drug resistance as a HER inhibitor (46). In addition, LAPATINIB has been shown to reduce the expression of CCND1 in breast cancer (47). Here, the expression of CCND1 increased, suggesting that LAPATINIB could be a potential therapeutic agent for PTC by targeting CCND1. Moreover, we analyzed the interaction between AKR1C3 and DOXORUBICIN, a commonly used chemotherapy drug for PTC, and found that targeting AKR1C3 might inhibit DOXORUBICIN drug resistance in PTC cells and enhance the efficacy of DOXORUBICIN treatment in PTC (48, 49).

ROC analysis of the 56 genes from the GSE33630 and GSE138198 datasets revealed four key genes shared by HT and PTC with AUC values of > 0.95 (i.e., ADH1B, ABR, SERPINA1, and LPAR5). ADH1B is related to many phenotypic traits, including alcohol metabolism, liver function, and cancer (50, 51). Moreover, single nucleotide polymorphisms (SNPs) in ADH1B have been correlated with esophageal square cell carcinoma (ESCC), colorectal cancer (CRC), and overall cancer. However, ADH1B SNPs are not associated with PTC or HT (5155). Furthermore, ADH1B is primarily a metabolism-related gene that is involved in the physiological processes of ethanol, uric acid, fat, and glucose metabolism as well as in the regulation of peripheral blood lymphocyte proliferation and nerve axons protection(5658). ADH1B serves as an autotarget antigen for Graves’ disease, an autoimmune inflammatory disease (59). Therefore, we posit that ADH1B may represent the key gene responsible for HT and PTC interaction.

ABR, located on chromosome 17p, is closely associated with chronic and acute myeloid and acute lymphoblastic leukemias (60, 61). ABR protein expression levels are the highest in the central nervous system and may interact simultaneously or sequentially with members of the Rho family to regulate and coordinate cellular signaling (62, 63). To the best of our knowledge, this study is the first, to report a relationship between ABR and HT or PTC. We hypothesize that ABR may influence HT or PTC development via myelination regulation.

SERPINA1, a protein associated with inflammatory disease, organismal injury, and abnormal thyroid function, is significantly lower in hyperthyroidism samples than in normal thyroid samples (64, 65). As a hub gene for PTC-DEGs, SERPINA1 is potential diagnostic gene for PTC and is related to the extracellular matrix pathway (6668). Our findings further suggest that SERPINA1 serves as an important cross-linking gene involved in the pathogenesis of HT and PTC.

Lysophosphatidic acid receptor 5 (LPAR5) overexpression is involved in mediating thyroid cancer progression; however, the mechanism underlying this process requires further elucidation (69, 70). Additionally, LPAR5 is reportedly associated with lysophosphatidic acid-induced pro-inflammatory signaling cascades in microglia (7173). Therefore, LPAR5 may be involved in a common pathway between HT and PTC, which represents a promising area of future investigation.

Furthermore, interaction networks containing 88 nodes and 100 edges were constructed to identify potential regulatory mechanisms targeting the four key genes, from which 21 miRNAs and 55 TFs targeting key genes were predicted. In the miRNA network, numerous miRNAs targeting the ABR were involved in PCT onset and development regulation, for instance; miR-34-5p (as a plasma exosomal miRNA) may assist identify benign malignant thyroid nodules and is promising for early liquid biopsy of PTC (74). Another study found that miR-34a-5p induction promotes apoptosis and inhibits human PTC cell line proliferation and migration (K1 and TPC-1) (75). The n384546/miR-145-5p/AKT3 pathway (where miR-145-5p is located), may inhibit PTC cell proliferation, invasion and migration in vitro and in vivo (76). MiR-506-3p may inhibit PTC proliferation by suppressing YAP1 expression and regulating the YAP1-CDK2/Cy clinical E1 cell cycle pathway (77). MiR-449a may partially inhibit PTC progression by downregulating metamucin (MTDH) (78). Numerous miRNAs targeting ABR hold diagnostic and therapeutic significance for PTC, and the ABR gene, as one of the prognostic genes of PTC, has not been targeted in relevant studies and deserves further exploration.

In the transcription factor network, all three prognostic genes were regulated by EGR1 (TF). Guo and Zhang (79) reported that the expression of EGR1/2 affects the proliferation of PTC cells and is related to poor prognosis(79). The transcription factor, Yin Yang 1 (YY1), that targets ADH1B, was used to identify 88% patients with PTC and is expected to play a role in the differential diagnosis of PTC (80). Co-targeting MYC that regulates SERPINA1 and ABR can promote PTC proliferation by regulating ANXA1 and repressing lncRNA PAX8-AS1:28 expression (81, 82). Androgen receptor (AR) expression changes affect PTC progression. Furthermore, AR axis suppression in PTC patients may contribute to the aggressive behavior of PTC (83). Higher ESR1 expression in PTC is associated with poorer prognosis and lower overall survival in women (84). SMAD3-activation of the TGF-ß/Smad3 pathway, which jointly targets ABR and LPAR5, could alter tumor cell function in PTC by suppressing NIS expression and altering tumor cell function (85). Another related gene, RUNX1, was upregulated in PTC tissues and expression levels correlated with PTC staging. RUNX1 knockdown could inhibit PTC proliferation, metastasis, and invasion (86).

Finally, we explored compounds that effectively target the four key genes. Currently, few studies within the DGIdb database report potential targeting compounds of these four key genes. Among compounds targeting ADH1B, alcohol intake was negatively associated with PTC cancer risk (87). Fomepizole inhibits the initial metabolism of ethylene glycol and methanol via ADH1B inhibition (88). There is no related research on acetaldehyde and abacavir that warrants further discussion. Meanwhile, human alpha-1-proteinase inhibitor (89), a potential targeting compound of SERPINA1, elicits a therapeutic effect on progressive ultimate fatal emphysema by inhibiting neutrophil elastase in the lung. Additionally, CHEMBL1630084, which blocked melanoma lung metastasis via LPA5 activation in a murine model (90). LPAR5 was highly expressed in PTC, and the lungs are important metastatic target organs related to PTC. Thus, CHEMBL1630084 and LPAR5 may provide protective effects against PTC pulmonary metastasis, and warrants further investigation. Furthermore, arachidonoylglycine has diagnostic potential for PTC (91), but the specific reason has not been studied. In this study, arachidonoylglycine was a compound targeting LPAR5, which provides insights for the study of related diagnostic mechanisms.

Although we performed a rigorous bioinformatics analysis, the number of PTC samples in the HT dataset was limited. To increase the credibility of the results, we verified the expression of key genes in clinical samples from patients with PTC in HT using qRT-PCR and IHC. However, there are few limitations in this study, including the small sample size. In the future, we will conduct a comparative analysis of PTC with HT vs HT by incorporating more HT samples, and expand our sample size to facilitate a more in-depth analysis of the genetic differences and similarities in patients with PTC and HT.

5 Conclusions

Using comprehensive bioinformatics methods, qRT-PCR, and IHC, we successfully screened, analyzed, and verified four key genes shared between HT and PTC. Specifically, ABR, LPAR5, and SERPINA1 were significantly up-regulated, while ADH1B was significantly down-regulated in clinical PTC samples compared to normal samples. Furthermore, we investigated the underlying regulatory mechanism of these four important genes as well as prospective drugs targeting them. We evaluated the common pathogenesis and molecular mechanism underlying the occurrence of HT and PTC and provided a foundation for further characterizing the relationship in these conditions. In conclusion, our findings advance the current understanding of the molecular etiopathogenesis of HT and PTC, and support the potential clinical applications of candidate genes in their treatment.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Ethics statement

The studies involving human participants were reviewed and approved by the ethics committee of scientific research and clinical trial of The First Affiliated Hospital of Zhengzhou University. The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.

Author contributions

TTL performed the data analysis and wrote the manuscript. NW, NL, and GD contributed to the data analysis and manuscript revision. MFP contributed to the literature search and data extraction. NW and DTY proofread the manuscript. TTL and DTY conceived and designed the study. All authors have contributed to the article and have approved the submitted version.

Funding

General Project of Natural Science Foundation of Henan Province; Award Number: 222300420568; Recipient, DY. Key Medical Science and Technology Project of Henan Province; Award Number: SBGJ202101014; Recipient, DY. Major Scientific Research Projects of Traditional Chinese Medicine in Henan Province; Award Number: 20-21ZYZD14; Recipient, DY. Cultivation of Young and Middle-aged Health Science and Technology Innovation Leading Talents in Henan Province; Award Number: YXKC2020015; Recipient, DY. Key University Science Research Project of Henan Province, Department of Education of Henan Province, Award Number: 23A320016; Recipient, NW.

Acknowledgments

We offer special thanks to DTY for his support and advice, GD for his review, and NW for her help with revising this manuscript. Thanks to my daughter, her hugs were what kept me going through this project.

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.2023.1140094/full#supplementary-material

Supplementary Figure 1 | (A) Receiver operating characteristic (ROC) curves of the four key genes (ADH1B, ABR, SERPINA1, and LPAR5) between papillary thyroid carcinoma (PTC) and normal samples in the GSE33630 dataset. (B) ROC curves of the four key genes between Hashimoto’s thyroiditis (HT) and normal samples in the GSE138198 dataset.

Supplementary Figure 2 | The miRNA/TF-key gene-drug network created based on common key genes in papillary thyroid carcinoma (PTC) and Hashimoto’s thyroiditis (HT). The orange arrows represent the key genes in PTC and HT, the green diamonds represent miRNAs, the purple circles represent the transcription factors (TFs), and the red rectangles represent the compounds.

References

1. Zhang QY, Ye XP, Zhou Z, Zhu CF, Li R, Fang Y, et al. Lymphocyte infiltration and thyrocyte destruction are driven by stromal and immune cell components in hashimoto's thyroiditis. Nat Commun (2022) 13(1):775. doi: 10.1038/s41467-022-28120-2

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Ralli M, Angeletti D, Fiore M, D'aguanno V, Lambiase A, Artico M, et al. Hashimoto's thyroiditis: an update on pathogenic mechanisms, diagnostic protocols, therapeutic strategies, and potential malignant transformation. Autoimmun Rev (2020) 19(10):102649. doi: 10.1016/j.autrev.2020.102649

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Chen YK, Lin CL, Cheng FT, Sung FC, Kao CH. Cancer risk in patients with hashimoto's thyroiditis: a nationwide cohort study. Br J Cancer (2013) 109(9):2496–501. doi: 10.1038/bjc.2013.597

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Hu X, Wang X, Liang Y, Chen X, Zhou S, Fei W, et al. Cancer risk in hashimoto's thyroiditis: a systematic review and meta-analysis. Front Endocrinol (Lausanne) (2022) 13:937871. doi: 10.3389/fendo.2022.937871

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Niu H, Yang J, Yang K, Huang Y. The relationship between RASSF1A promoter methylation and thyroid carcinoma: a meta-analysis of 14 articles and a bioinformatics of 2 databases (PRISMA). Med (Baltimore) (2017) 96(46):e8630. doi: 10.1097/MD.0000000000008630

CrossRef Full Text | Google Scholar

6. Li H, Weng J, Shi Y, Gu W, Mao Y, Wang Y, et al. An improved deep learning approach for detection of thyroid papillary cancer in ultrasound images. Sci Rep (2018) 8(1):6600. doi: 10.1038/s41598-018-25005-7

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Sulaieva O, Selezniov O, Shapochka D, Belemets N, Nechay O, Chereshneva Y, et al. Hashimoto's thyroiditis attenuates progression of papillary thyroid carcinoma: deciphering immunological links. Heliyon (2020) 6(1):e03077. doi: 10.1016/j.heliyon.2019.e03077

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Jauculan MC, Buenaluz-Sedurante M, Jimeno CA. Risk factors associated with disease recurrence among patients with low-risk papillary thyroid cancer treated at the university of the Philippines-Philippine general hospital. Endocrinol Metab (Seoul) (2016) 31(1):113–9. doi: 10.3803/EnM.2016.31.1.113

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Kim EY, Kim WG, Kim WB, Kim TY, Kim JM, Ryu JS, et al. Coexistence of chronic lymphocytic thyroiditis is associated with lower recurrence rates in patients with papillary thyroid carcinoma. Clin Endocrinol (Oxf) (2009) 71(4):581–6. doi: 10.1111/j.1365-2265.2009.03537.x

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Chen DS, Mellman I. Oncology meets immunology: the cancer-immunity cycle. Immunity (2013) 39(1):1–10. doi: 10.1016/j.immuni.2013.07.012

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Oh CM, Park S, Lee JY, Won YJ, Shin A, Kong HJ, et al. Increased prevalence of chronic lymphocytic thyroiditis in Korean patients with papillary thyroid cancer. PLoS One (2014) 9(6):e99054. doi: 10.1371/journal.pone.0099054

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Dailey ME, Lindsay S, Skahen R. Relation of thyroid neoplasms to hashimoto disease of the thyroid gland. AMA Arch Surg (1955) 70(2):291–7. doi: 10.1001/archsurg.1955.01270080137023

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Konturek A, Barczyński M, Wierzchowski W, Stopa M, Nowak W. Coexistence of papillary thyroid cancer with hashimoto thyroiditis. Langenbecks Arch Surg (2013) 398(3):389–94. doi: 10.1007/s00423-012-1021-x

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Lee JH, Kim Y, Choi JW, Kim YS. The association between papillary thyroid carcinoma and histologically proven hashimoto's thyroiditis: a meta-analysis. Eur J Endocrinol (2013) 168(3):343–9. doi: 10.1530/EJE-12-0903

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Lai X, Xia Y, Zhang B, Li J, Jiang Y. A meta-analysis of hashimoto's thyroiditis and papillary thyroid carcinoma risk. Oncotarget (2017) 8(37):62414–24. doi: 10.18632/oncotarget.18620

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Silva De Morais N, Stuart J, Guan H, Wang Z, Cibas ES, Frates MC, et al. The impact of hashimoto thyroiditis on thyroid nodule cytology and risk of thyroid cancer. J Endocr Soc (2019) 3(4):791–800. doi: 10.1210/js.2018-00427

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Resende De Paiva C, Grønhøj C, Feldt-Rasmussen U, Von Buchwald C. Association between hashimoto's thyroiditis and thyroid cancer in 64,628 patients. Front Oncol (2017) 7:53. doi: 10.3389/fonc.2017.00053

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Ragusa F, Fallahi P, Elia G, Gonnella D, Paparo SR, Giusti C, et al. Hashimotos' thyroiditis: epidemiology, pathogenesis, clinic and therapy. Best Pract Res Clin Endocrinol Metab (2019) 33(6):101367. doi: 10.1016/j.beem.2019.101367

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Dorris ER, Smyth P, O'leary JJ, Sheils O. MIR141 expression differentiates hashimoto thyroiditis from PTC and benign thyrocytes in Irish archival thyroid tissues. Front Endocrinol (Lausanne) (2012) 3:102. doi: 10.3389/fendo.2012.00102

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Yan X, He B, Hu L, Gao J, Chen S, Jiang G. Insight into the endocrine disrupting effect and cell response to butyltin compounds in H295R cell: evaluated with proteomics and bioinformatics analysis. Sci Total Environ (2018) 628-629:1489–96. doi: 10.1016/j.scitotenv.2018.02.165

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Dom G, Tarabichi M, Unger K, Thomas G, Oczko-Wojciechowska M, Bogdanova T, et al. A gene expression signature distinguishes normal tissues of sporadic and radiation-induced papillary thyroid carcinomas. Br J Cancer (2012) 107(6):994–1000. doi: 10.1038/bjc.2012.302

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Subhi O, Schulten HJ, Bagatian N, Al-Dayini R, Karim S, Bakhashab S, et al. Genetic relationship between hashimoto`s thyroiditis and papillary thyroid carcinoma with coexisting hashimoto`s thyroiditis. PLoS One (2020) 15(6):e0234566. doi: 10.1371/journal.pone.0234566

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Integrated genomic characterization of papillary thyroid carcinoma. Cell (2014) 159(3):676–90. doi: 10.1186/1471-2105-9-559

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Langfelder P, Horvath S. WGCNA: an r package for weighted correlation network analysis. BMC Bioinf (2008) 9:559. doi: 10.1111/j.2517-6161.1995.tb02031.x

CrossRef Full Text | Google Scholar

25. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Royal Stat Soc: Series B (Methodological) (1995) 57(1):289–300. doi: 10.1093/nar/gkv007

CrossRef Full Text | Google Scholar

26. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res (2015) 43(7):e47. doi: 10.1093/nar/gkv007

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Wu T, Hu E, Xu S, Chen M, Guo P, Dai Z, et al. clusterProfiler 4.0: a universal enrichment tool for interpreting omics data. Innovation (Camb) (2021) 2(3):100141. doi: 10.1101/gr.1239303

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res (2003) 13(11):2498–504. doi: 10.1007/s12022-022-09707-3

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Baloch ZW, Asa SL, Barletta JA, Ghossein RA, Juhlin CC, Jung CK, et al. Overview of the 2022 WHO classification of thyroid neoplasms. Endocr Pathol (2022) 33(1):27–63. doi: 10.1006/meth.2001.1262

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-delta delta C(T)) method. Methods (2001) 25(4):402–8. doi: 10.1002/cam4.4920

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Remmele W, Stegner HE. [Recommendation for uniform definition of an immunoreactive score (IRS) for immunohistochemical estrogen receptor detection (ER-ICA) in breast cancer tissue]. Pathologe (1987) 8(3):138–40.

PubMed Abstract | Google Scholar

32. Liotti F, Visciano C, Melillo RM. Inflammation in thyroid oncogenesis. Am J Cancer Res (2012) 2(3):286–97. doi: 10.14639/0392-100X-N1081

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Hanege FM, Tuysuz O, Celik S, Sakallıoglu O, Arslan Solmaz O. Hashimoto's thyroiditis in papillary thyroid carcinoma: a 22-year study. Acta Otorhinolaryngol Ital (2021) 41(2):142–5. doi: 10.1001/jamanetworkopen.2021.18526

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Xu S, Huang H, Qian J, Liu Y, Huang Y, Wang X, et al. Prevalence of hashimoto thyroiditis in adults with papillary thyroid cancer and its association with cancer recurrence and outcomes. JAMA Netw Open (2021) 4(7):e2118526. doi: 10.1038/s41419-019-1600-7

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Fu X, Zhang H, Chen Z, Yang Z, Shi D, Liu T, et al. TFAP2B overexpression contributes to tumor growth and progression of thyroid cancer through the COX-2 signaling pathway. Cell Death Dis (2019) 10(6):397. doi: 10.3389/fimmu.2017.00521

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Fröhlich E, Wahl R. Thyroid autoimmunity: role of anti-thyroid antibodies in thyroid and extra-thyroidal diseases. Front Immunol (2017) 8:521. doi: 10.1210/jc.2007-1199

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Latrofa F, Ricci D, Grasso L, Vitti P, Masserini L, Basolo F, et al. Characterization of thyroglobulin epitopes in patients with autoimmune and non-autoimmune thyroid diseases using recombinant human monoclonal thyroglobulin autoantibodies. J Clin Endocrinol Metab (2008) 93(2):591–6. doi: 10.1016/j.tem.2014.09.001

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Ehlers M, Schott M. Hashimoto's thyroiditis and papillary thyroid cancer: are they immunologically linked? Trends Endocrinol Metab (2014) 25(12):656–64. doi: 10.1038/s41598-020-67615-0

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Mohamed SY, Ibrahim TR, Elbasateeny SS, Abdelaziz LA, Farouk S, Yassin MA, et al. Clinicopathological characterization and prognostic implication of FOXP3 and CK19 expression in papillary thyroid carcinoma and concomitant hashimoto's thyroiditis. Sci Rep (2020) 10(1):10651. doi: 10.3892/etm.2016.3854

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Xia N, Chen G, Liu M, Ye X, Pan Y, Ge J, et al. Anti-inflammatory effects of luteolin on experimental autoimmune thyroiditis in mice. Exp Ther Med (2016) 12(6):4049–54. doi: 10.4158/EP-2019-0162

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Liu J, Fu J, Jia Y, Yang N, Li J, Wang G. SERUM METABOLOMIC PATTERNS IN PATIENTS WITH AUTOIMMUNE THYROID DISEASE. Endocr Pract (2020) 26(1):82–96. doi: 10.1016/j.bbrc.2019.09.138

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Li Q, Wang XJ, Jin JH. SOX2-induced upregulation of lncRNA LINC01510 promotes papillary thyroid carcinoma progression by modulating miR-335/SHH and activating hedgehog pathway. Biochem Biophys Res Commun (2019) 520(2):277–83. doi: 10.1530/ERC-19-0381

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Zhang W, Zhang H, Zhao X. circ_0005273 promotes thyroid carcinoma progression by SOX2 expression. Endocr Relat Cancer (2020) 27(1):11–21. doi: 10.3892/or.2021.8001

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Qi Y, He J, Zhang Y, Wang L, Yu Y, Yao B, et al. Circular RNA hsa_circ_0001666 sponges miR−330−5p, miR−193a−5p and miR−326, and promotes papillary thyroid carcinoma progression via upregulation of ETV4. Oncol Rep (2021) 45(4). doi: 10.1186/s12957-017-1190-8

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Chen W, Liu Q, Lv Y, Xu D, Chen W, Yu J. Special role of JUN in papillary thyroid carcinoma based on bioinformatics analysis. World J Surg Oncol (2017) 15(1):119. doi: 10.18632/oncotarget.15773

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Cheng L, Jin Y, Liu M, Ruan M, Chen L. HER inhibitor promotes BRAF/MEK inhibitor-induced redifferentiation in papillary thyroid cancer harboring BRAFV600E. Oncotarget (2017) 8(12):19843–54. doi: 10.3892/or.2016.5073

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Ohnishi Y, Yasui H, Kakudo K, Nozaki M. Lapatinib-resistant cancer cells possessing epithelial cancer stem cell properties develop sensitivity during sphere formation by activation of the ErbB/AKT/cyclin D2 pathway. Oncol Rep (2016) 36(5):3058–64. doi: 10.1111/bcp.12021

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Voon PJ, Yap HL, Ma CY, Lu F, Wong AL, Sapari NS, et al. Correlation of aldo-ketoreductase (AKR) 1C3 genetic variant with doxorubicin pharmacodynamics in Asian breast cancer patients. Br J Clin Pharmacol (2013) 75(6):1497–505. doi: 10.1016/j.biopha.2014.12.022

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Zhong T, Xu F, Xu J, Liu L, Chen Y. Aldo-keto reductase 1C3 (AKR1C3) is associated with the doxorubicin resistance in human breast cancer via PTEN loss. BioMed Pharmacother (2015) 69:317–25. doi: 10.1002/ajmg.b.32523

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Polimanti R, Gelernter J. ADH1B: from alcoholism, natural selection, and cancer to the human phenome. Am J Med Genet B Neuropsychiatr Genet (2018) 177(2):113–25. doi: 10.1002/ajmg.b.32523

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Seol JE, Kim J, Lee BH, Hwang DY, Jeong J, Lee HJ, et al. Folate, alcohol, ADH1B and ALDH2 and colorectal cancer risk. Public Health Nutr (2020), 30: 1–8. doi: 10.1053/j.gastro.2009.07.070

CrossRef Full Text | Google Scholar

52. Cui R, Kamatani Y, Takahashi A, Usami M, Hosono N, Kawaguchi T, et al. Functional variants in ADH1B and ALDH2 coupled with alcohol and smoking synergistically enhance esophageal cancer risk. Gastroenterology (2009) 137(5):1768–75. doi: 10.1093/carcin/bgy011

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Offermans NSM, Ketcham SM, Van Den Brandt PA, Weijenberg MP, Simons C. Alcohol intake, ADH1B and ADH1C genotypes, and the risk of colorectal cancer by sex and subsite in the Netherlands cohort study. Carcinogenesis (2018) 39(3):375–88. doi: 10.1042/BSR20181915

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Tan B, Ning N. Association of ADH1B Arg47His polymorphism with the risk of cancer: a meta-analysis. Biosci Rep (2019) 39(4):754–762. doi: 10.4143/crt.2020.478

CrossRef Full Text | Google Scholar

55. Choi CK, Shin MH, Cho SH, Kim HY, Zheng W, Long J, et al. Association between ALDH2 and ADH1B polymorphisms and the risk for colorectal cancer in koreans. Cancer Res Treat (2021) 53(3):754–62. doi: 10.1539/joh.45.408

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Tian Y, Ishikawa H, Piao FY, Yamamoto H, Yamauchi T, Duan ZW, et al. Micronucleus assay of human lymphocytes: a comparison of cytokinesis-block and human capillary blood lymphocytes methods. J Occup Health (2003) 45(6):408–9. doi: 10.1016/j.mrfmmm.2006.11.026

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Ishikawa H, Ishikawa T, Yamamoto H, Fukao A, Yokoyama K. Genotoxic effects of alcohol in human peripheral lymphocytes modulated by ADH1B and ALDH2 gene polymorphisms. Mutat Res (2007) 615(1-2):134–42. doi: 10.4238/gmr.15048740

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Ma L, Lu ZN. Role of ADH1B rs1229984 and ALDH2 rs671 gene polymorphisms in the development of alzheimer's disease. Genet Mol Res (2016) 15(4):747–759. doi: 10.1136/bjophthalmol-2015-306634

CrossRef Full Text | Google Scholar

59. Cheng KC, Hung CT, Cheng KY, Chen KJ, Wu WC, Suen JL, et al. Proteomic surveillance of putative new autoantigens in thyroid orbitopathy. Br J Ophthalmol (2015) 99(11):1571–6. doi: 10.1093/nar/17.21.8821

PubMed Abstract | CrossRef Full Text | Google Scholar

60. Heisterkamp N, Morris C, Groffen J. ABR, an active BCR-related gene. Nucleic Acids Res (1989) 17(21):8821–31. doi: 10.18632/oncotarget.22093

PubMed Abstract | CrossRef Full Text | Google Scholar

61. Namasu CY, Katzerke C, Bräuer-Hartmann D, Wurm AA, Gerloff D, Hartmann JU, et al. ABR, a novel inducer of transcription factor C/EBPα, contributes to myeloid differentiation and is a favorable prognostic factor in acute myeloid leukemia. Oncotarget (2017) 8(61):103626–39. doi: 10.1016/S0021-9258(19)85281-X

PubMed Abstract | CrossRef Full Text | Google Scholar

62. Heisterkamp N, Kaartinen V, Van Soest S, Bokoch GM, Groffen J. Human ABR encodes a protein with GAPrac activity and homology to the DBL nucleotide exchange factor domain. J Biol Chem (1993) 268(23):16903–6. doi: 10.1073/pnas.92.22.10282

PubMed Abstract | CrossRef Full Text | Google Scholar

63. Chuang TH, Xu X, Kaartinen V, Heisterkamp N, Groffen J, Bokoch GM. Abr and bcr are multifunctional regulators of the rho GTP-binding protein family. Proc Natl Acad Sci U S A (1995) 92(22):10282–6. doi: 10.3390/molecules25122831

PubMed Abstract | CrossRef Full Text | Google Scholar

64. Masood A, Benabdelkamel H, Ekhzaimy AA, Alfadda AA. Plasma-based proteomics profiling of patients with hyperthyroidism after antithyroid treatment. Molecules (2020) 25(12):2831. doi: 10.1097/CAD.0000000000001108

PubMed Abstract | CrossRef Full Text | Google Scholar

65. Luo L, Zheng W, Chen C, Sun S. Searching for essential genes and drug discovery in breast cancer and periodontitis via text mining and bioinformatics analysis. Anticancer Drugs (2021) 32(10):1038–45. doi: 10.1186/1755-8794-4-30

PubMed Abstract | CrossRef Full Text | Google Scholar

66. Vierlinger K, Mansfeld MH, Koperek O, Nöhammer C, Kaserer K, Leisch F. Identification of SERPINA1 as single marker for papillary thyroid carcinoma through microarray meta analysis and quantification of its discriminatory power in independent validation. BMC Med Genomics (2011) 4:30. doi: 10.7150/thno.26862

PubMed Abstract | CrossRef Full Text | Google Scholar

67. Teng H, Mao F, Liang J, Xue M, Wei W, Li X, et al. Transcriptomic signature associated with carcinogenesis and aggressiveness of papillary thyroid carcinoma. Theranostics (2018) 8(16):4345–58. doi: 10.7150/ijms.63402

PubMed Abstract | CrossRef Full Text | Google Scholar

68. Liu C, Pan Y, Li Q, Zhang Y. Bioinformatics analysis identified shared differentially expressed genes as potential biomarkers for hashimoto's thyroiditis-related papillary thyroid cancer. Int J Med Sci (2021) 18(15):3478–87. doi: 10.1042/BSR20190083

PubMed Abstract | CrossRef Full Text | Google Scholar

69. Zhang S, Wang Q, Han Q, Han H, Lu P. Identification and analysis of genes associated with papillary thyroid carcinoma by bioinformatics methods. Biosci Rep (2019) 39(4). doi: 10.1111/cas.14837

CrossRef Full Text | Google Scholar

70. Zhao WJ, Zhu LL, Yang WQ, Xu SJ, Chen J, Ding XF, et al. LPAR5 promotes thyroid carcinoma cell proliferation and migration by activating class IA PI3K catalytic subunit p110β. Cancer Sci (2021) 112(4):1624–32. doi: 10.3389/fncel.2019.00531

PubMed Abstract | CrossRef Full Text | Google Scholar

71. Plastira I, Joshi L, Bernhart E, Schoene J, Specker E, Nazare M, et al. Small-molecule lysophosphatidic acid receptor 5 (LPAR5) antagonists: versatile pharmacological tools to regulate inflammatory signaling in BV-2 microglia cells. Front Cell Neurosci (2019) 13:531. doi: 10.1080/21655979.2021.1983975

PubMed Abstract | CrossRef Full Text | Google Scholar

72. Villalpando-Rodriguez GE, Blankstein AR, Konzelman C, Gibson SB. Lysosomal destabilizing drug siramesine and the dual tyrosine kinase inhibitor lapatinib induce a synergistic ferroptosis through reduced heme oxygenase-1 (HO-1) levels. Oxid Med Cell Longevity (2019) 2019:9561281. doi: 10.3389/fgene.2020.00449

CrossRef Full Text | Google Scholar

73. Xu L, Xu C, Lin X, Lu H, Cai Y. Interference with lysophosphatidic acid receptor 5 ameliorates oxidized low-density lipoprotein-induced human umbilical vein endothelial cell injury by inactivating NOD-like receptor family, pyrin domain containing 3 inflammasome signaling. Bioengineered (2021) 12(1):8089–99. doi: 10.1186/s12576-021-00793-2

PubMed Abstract | CrossRef Full Text | Google Scholar

74. Liang M, Yu S, Tang S, Bai L, Cheng J, Gu Y, et al. A panel of plasma exosomal miRNAs as potential biomarkers for differential diagnosis of thyroid nodules. Front Genet (2020) 11:449. doi: 10.1038/s41419-019-1637-7

PubMed Abstract | CrossRef Full Text | Google Scholar

75. Sun W, Yin D. Long noncoding RNA CASC7 inhibits the proliferation and migration of papillary thyroid cancer cells by inhibiting miR-34a-5p. J Physiol Sci (2021) 71(1):9. doi: 10.1016/j.prp.2020.153231

PubMed Abstract | CrossRef Full Text | Google Scholar

76. Feng J, Zhou Q, Yi H, Ma S, Li D, Xu Y, et al. A novel lncRNA n384546 promotes thyroid papillary cancer progression and metastasis by acting as a competing endogenous RNA of miR-145-5p to regulate AKT3. Cell Death Dis (2019) 10(6):433. doi: 10.1097/CAD.0000000000001340

PubMed Abstract | CrossRef Full Text | Google Scholar

77. Chen L, Wang X, Ji C, Hu J, Fang L. MiR-506-3p suppresses papillary thyroid cancer cells tumorigenesis by targeting YAP1. Pathol Res Pract (2020) 216(12):153231. doi: 10.1007/s10528-021-10075-6

PubMed Abstract | CrossRef Full Text | Google Scholar

78. Feng Y, Yang X, Wang Y, Chi N, Yu J, Fu X. circRNA mannosidase alpha class 1A member 2 contributes to the proliferation and motility of papillary thyroid cancer cells through upregulating metadherin via absorbing microRNA-449a. Anticancer Drugs (2023) 34(1):44–56. doi: 10.1007/s12022-015-9359-6

PubMed Abstract | CrossRef Full Text | Google Scholar

79. Guo H, Zhang L. EGR1/2 inhibits papillary thyroid carcinoma cell growth by suppressing the expression of PTEN and BAX. Biochem Genet (2021) 59(6):1544–57. doi: 10.1007/s10528-021-10075-6

PubMed Abstract | CrossRef Full Text | Google Scholar

80. Arribas J, Castellví J, Marcos R, Zafón C, Velázquez A. Expression of YY1 in differentiated thyroid cancer. Endocr Pathol (2015) 26(2):111–8. doi: 10.1080/21655979.2021.1996511

PubMed Abstract | CrossRef Full Text | Google Scholar

81. Zhang Y, Li F, Chen J. MYC promotes the development of papillary thyroid carcinoma by inhibiting the expression of lncRNA PAX8−AS1:28. Oncol Rep (2019) 41(4):2511–7. doi: 10.3390/cancers12051109

PubMed Abstract | CrossRef Full Text | Google Scholar

82. Ying X, Chen L, Xie J, Hu Y, Wu Q, Cao L, et al. ANXA1 (Annexin A1) regulated by MYC (MYC proto-oncogene) promotes the growth of papillary thyroid carcinoma. Bioengineered (2021) 12(2):9251–65. doi: 10.1245/s10434-017-5780-z

PubMed Abstract | CrossRef Full Text | Google Scholar

83. Chou CK, Chi SY, Chou FF, Huang SC, Wang JH, Chen CC, et al. Aberrant expression of androgen receptor associated with high cancer risk and extrathyroidal extension in papillary thyroid carcinoma. Cancers (Basel) (2020) 12(5):1109. doi: 10.1155/2021/6840642

PubMed Abstract | CrossRef Full Text | Google Scholar

84. Yi JW, Kim SJ, Kim JK, Seong CY, Yu HW, Chai YJ, et al. Upregulation of the ESR1 gene and ESR ratio (ESR1/ESR2) is associated with a worse prognosis in papillary thyroid carcinoma: the impact of the estrogen receptor α/β expression on clinical outcomes in papillary thyroid carcinoma patients. Ann Surg Oncol (2017) 24(12):3754–62. doi: 10.21037/gs-22-219

PubMed Abstract | CrossRef Full Text | Google Scholar

85. Yan Z, Yangyanqiu W, Shuwen H, Jing M, Haihong L, Gong C, et al. Downregulation of Rap1GAP expression activates the TGF-β/Smad3 pathway to inhibit the expression of Sodium/Iodine transporter in papillary thyroid carcinoma cells. BioMed Res Int (2021) 2021:6840642. doi: 10.1007/978-3-319-98788-0_1

PubMed Abstract | CrossRef Full Text | Google Scholar

86. Xu J, Zheng G, Guo H, Meng K, Zhang W, He R, et al. Bioinformatics analysis of downstream circRNAs and miRNAs regulated by runt-related transcription factor 1 in papillary thyroid carcinoma. Gland Surg (2022) 11(5):868–81. doi: 10.21037/gs-22-219

PubMed Abstract | CrossRef Full Text | Google Scholar

87. Huang H, Zhao N, Chen Y, Deziel N, Dai M, Li N, et al. Alcohol consumption and risk of thyroid cancer: a population based case-control study in Connecticut. Adv Exp Med Biol (2018) 1032:1–14. doi: 10.1007/978-3-319-98788-0_1

PubMed Abstract | CrossRef Full Text | Google Scholar

88. Eriksson CJ, Fukunaga T, Sarkola T, Chen WJ, Chen CC, Ju JM, et al. Functional relevance of human adh polymorphism. Alcohol Clin Exp Res (2001) 25(5 Suppl ISBRA):157s–63s. doi: 10.1111/j.1530-0277.2001.tb02391.x

PubMed Abstract | CrossRef Full Text | Google Scholar

89. Karnaukhova E, Ophir Y, Golding B. Recombinant human alpha-1 proteinase inhibitor: towards therapeutic use. Amino Acids (2006) 30(4):317–32. doi: 10.1007/s00726-005-0324-4

PubMed Abstract | CrossRef Full Text | Google Scholar

90. Gupte R, Siddam A, Lu Y, Li W, Fujiwara Y, Panupinthu N, et al. Synthesis and pharmacological evaluation of the stereoisomers of 3-carba cyclic-phosphatidic acid. Bioorg Med Chem Lett (2010) 20(24):7525–8. doi: 10.1016/j.bmcl.2010.09.115

PubMed Abstract | CrossRef Full Text | Google Scholar

91. Sun Q, Zhao H, Liu Z, Wang F, He Q, Xiu C, et al. Identifying potential metabolic tissue biomarkers for papillary thyroid cancer in different iodine nutrient regions. Endocrine (2021) 74(3):582–91. doi: 10.1007/s12020-021-02773-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: Hashimoto’s thyroiditis, papillary thyroid carcinoma, key genes, SERPINA1, LPAR5, ABR

Citation: Liu TT, Yin DT, Wang N, Li N, Dong G and Peng MF (2023) Identifying and analyzing the key genes shared by papillary thyroid carcinoma and Hashimoto’s thyroiditis using bioinformatics methods. Front. Endocrinol. 14:1140094. doi: 10.3389/fendo.2023.1140094

Received: 18 January 2023; Accepted: 09 May 2023;
Published: 31 May 2023.

Edited by:

Yury O. Nunez Lopez, AdventHealth, United States

Reviewed by:

Xiaofeng Wu, Pfizer, United States
Weihang Li, Fourth Military Medical University, China
Malgorzata Trofimiuk-Muldner, Jagiellonian University Medical College, Poland

Copyright © 2023 Liu, Yin, Wang, Li, Dong and Peng. 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: De-tao Yin, detaoyin@zzu.edu.cn

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.