Integrated bioinformatics analysis and screening of hub genes in papillary thyroid carcinoma

Identification of potential targetable genes in papillary, follicular, and anaplastic thyroid carcinoma using bioinformatics analysis

PURPOSE

To perform an extensive exploratory analysis to build a deeper insight into clinically relevant molecular biomarkers in Papillary, Follicular, and Anaplastic thyroid carcinomas (PTC, FTC, ATC).

METHODS

Thirteen Thyroid Cancer (THCA) datasets incorporating PTC, FTC, and ATC were derived from the Gene Expression Omnibus. Genes differentially expressed (DEGs) between THCA and normal were identified and subjected to GO and KEGG analyses. Multiple topological properties were harnessed and protein-protein interaction (PPI) networks were constructed to identify the hub genes followed by survival analysis and validation.

RESULTS

There were 70, 87, and 377 DEGs, and 23, 27, and 53 hub genes for PTC, FTC, and ATC samples, respectively. Survival analysis detected 39 overall and 49 relapse-free survival-relevant hub genes. Six hub genes, BCL2, FN1, ITPR1, LYVE1, NTRK2, TBC1D4, were found common to more than one THCA type. The most significant hub genes found in the study were: BCL2, CD44, DCN, FN1, IRS1, ITPR1, MFAP4, MKI67, NTRK2, PCLO, TGFA. The most enriched and significant GO terms were Melanocyte differentiation for PTC, Extracellular region for FTC, and Extracellular exosome for ATC. Prostate cancer for PTC was the most significantly enriched KEGG pathway. The results were validated using TCGA data.

CONCLUSIONS

The findings unravel potential biomarkers and therapeutic targets of thyroid carcinomas.

Integrated transcriptome sequencing and weighted gene co-expression network analysis reveals key genes of papillary thyroid carcinomas

Objective

Papillary thyroid carcinoma (PTC) accounts for the majority of thyroid cancers and has a high recurrence rate. We aimed to screen key genes involved in PTC to provide novel insights into the mechanisms of PTC.

Methods

Seven microarray datasets of PTC were downloaded from gene expression omnibus database. Differentially expressed genes (DEGs) between PTC and normal samples were screened in the merged dataset. Then, protein-protein interaction (PPIs) functional modules analysis and weighted gene co-expression network analysis (WGCNA) were utilized to identify PTC-associated key genes. The identified key genes were then characterized from various aspects, including gene set enrichment analysis (GSEA) and the associations with immune infiltration, methylation levels and prognosis.

Results

A large numbers of DEGs were identified, and these DEGs are involved in several cancer pathways. Nine key genes (including down-regulated genes GNA14, AVPR1A, and WFS1, and up-regulated genes LAMB3, PLAU, MET, MFGE8, PRSS23, and SERPINA1) were identified. Patients in the AVPR1A and GNA14 high expression groups had better disease-free survival (DFS) than those in the low expression group. Key genes were mainly involved in P53 pathway, estrogen response, apoptosis, glycolysis, NOTCH signaling, epithelial mesenchymal transition, WNT_beta catenin signaling, and inflammatory response. The expression of key genes was associated with immune cell infiltration and corresponding methylation levels. The verification results of key gene proteins and mRNA expression levels using external validation datasets were consistent with our expectations, implying the involvements of key genes in PTC.

Conclusion

The key genes may serve as potential therapeutic targets for PTC. This study provides novel insights into the mechanisms underlying PTC development.

Immunological characteristics of immunogenic cell death genes and malignant progression driving roles of TLR4 in anaplastic thyroid carcinoma

Anaplastic thyroid carcinoma (ATC) was a rare malignancy featured with the weak immunotherapeutic response. So far, disorders of immunogenic cell death genes (ICDGs) were identified as the driving factors in cancer progression, while their roles in ATC remained poorly clear. Datasets analysis identified that most ICDGs were high expressed in ATC, while DE-ICDGs were located in module c1_112, which was mainly enriched in Toll-like receptor signalings. Subsequently, the ICD score was established to classify ATC samples into the high and low ICD score groups, and function analysis indicated that high ICD score was associated with the immune characteristics. The high ICD score group had higher proportions of specific immune and stromal cells, as well as increased expression of immune checkpoints. Additionally, TLR4, ENTPD1, LY96, CASP1 and PDIA3 were identified as the dynamic signature in the malignant progression of ATC. Notably, TLR4 was significantly upregulated in ATC tissues, associated with poor prognosis. Silence of TLR4 inhibited the proliferation, metastasis and clone formation of ATC cells. Eventually, silence of TLR4 synergistically enhanced paclitaxel-induced proliferation inhibition, apoptosis, CALR exposure and release of ATP. Our findings highlighted that the aberrant expression of TLR4 drove the malignant progression of ATC, which contributed to our understanding of the roles of ICDGs in ATC.

Presentation of potential genes and deleterious variants associated with non-syndromic hearing loss: a computational approach

Non-syndromic hearing loss (NSHL) is a common hereditary disorder. Both clinical and genetic heterogeneity has created many obstacles to understanding the causes of NSHL. The present study has attempted to ravel the genetic aetiology in NSHL progression and to screen out potential target genes using computational approaches. The reported NSHL target genes (2009-2020) have been studied by analyzing different biochemical and signaling pathways, interpretation of their functional association network, and discovery of important regulatory interactions with three previously established miRNAs in the human inner ear as well as in NSHL such as miR-183, miR-182, and miR-96. This study has identified SMAD4 and SNAI2 as the most putative target genes of NSHL. But pathogenic and deleterious non-synonymous single nucleotide polymorphisms discovered within SMAD4 is anticipated to have an impact on NSHL progression. Additionally, the identified deleterious variants in the functional domains of SMAD4 added a supportive clue for further study. Thus, the identified deleterious variant i.e., rs377767367 (G491V) in SMAD4 needs further clinical validation. The present outcomes would provide insights into the genetics of NSHL progression.