Epidermal growth factor receptor activation confers resistance to lenvatinib in thyroid cancer cells

Exploration and validation of key genes associated with early lymph node metastasis in thyroid carcinoma using weighted gene co-expression network analysis and machine learning

Background

Thyroid carcinoma (THCA), the most common endocrine neoplasm, typically exhibits an indolent behavior. However, in some instances, lymph node metastasis (LNM) may occur in the early stages, with the underlying mechanisms not yet fully understood.

Materials and methods

LNM potential was defined as the tumor's capability to metastasize to lymph nodes at an early stage, even when the tumor volume is small. We performed differential expression analysis using the 'Limma' R package and conducted enrichment analyses using the Metascape tool. Co-expression networks were established using the 'WGCNA' R package, with the soft threshold power determined by the 'pickSoftThreshold' algorithm. For unsupervised clustering, we utilized the 'ConsensusCluster Plus' R package. To determine the topological features and degree centralities of each node (protein) within the Protein-Protein Interaction (PPI) network, we used the CytoNCA plugin integrated with the Cytoscape tool. Immune cell infiltration was assessed using the Immune Cell Abundance Identifier (ImmuCellAI) database. We applied the Least Absolute Shrinkage and Selection Operator (LASSO), Support Vector Machine (SVM), and Random Forest (RF) algorithms individually, with the 'glmnet,' 'e1071,' and 'randomForest' R packages, respectively. Ridge regression was performed using the 'oncoPredict' algorithm, and all the predictions were based on data from the Genomics of Drug Sensitivity in Cancer (GDSC) database. To ascertain the protein expression levels and subcellular localization of genes, we consulted the Human Protein Atlas (HPA) database. Molecular docking was carried out using the mcule 1-click Docking server online. Experimental validation of gene and protein expression levels was conducted through Real-Time Quantitative PCR (RT-qPCR) and immunohistochemistry (IHC) assays.

Results

Through WGCNA and PPI network analysis, we identified twelve hub genes as the most relevant to LNM potential from these two modules. These 12 hub genes displayed differential expression in THCA and exhibited significant correlations with the downregulation of neutrophil infiltration, as well as the upregulation of dendritic cell and macrophage infiltration, along with activation of the EMT pathway in THCA. We propose a novel molecular classification approach and provide an online web-based nomogram for evaluating the LNM potential of THCA (http://www.empowerstats.net/pmodel/?m=17617_LNM). Machine learning algorithms have identified ERBB3 as the most critical gene associated with LNM potential in THCA. ERBB3 exhibits high expression in patients with THCA who have experienced LNM or have advanced-stage disease. The differential methylation levels partially explain this differential expression of ERBB3. ROC analysis has identified ERBB3 as a diagnostic marker for THCA (AUC=0.89), THCA with high LNM potential (AUC=0.75), and lymph nodes with tumor metastasis (AUC=0.86). We have presented a comprehensive review of endocrine disruptor chemical (EDC) exposures, environmental toxins, and pharmacological agents that may potentially impact LNM potential. Molecular docking revealed a docking score of -10.1 kcal/mol for Lapatinib and ERBB3, indicating a strong binding affinity.

Conclusion

In conclusion, our study, utilizing bioinformatics analysis techniques, identified gene modules and hub genes influencing LNM potential in THCA patients. ERBB3 was identified as a key gene with therapeutic implications. We have also developed a novel molecular classification approach and a user-friendly web-based nomogram tool for assessing LNM potential. These findings pave the way for investigations into the mechanisms underlying differences in LNM potential and provide guidance for personalized clinical treatment plans.

Cuproptosis-related lncRNAs predict prognosis and immune response of thyroid carcinoma

Objective: To estimate the survival and prognosis of patients with thyroid carcinoma (THCA) based on the Long non-coding RNA (lncRNA) traits linked to cuproptosis and to investigate the connection between the immunological spectrum of THCA and medication sensitivity. Methods: RNA-Seq data and clinical information for THCA were obtained from the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. We built a risk prognosis model by identifying and excluding lncRNAs associated with cuproptosis using Cox regression and LASSO methods. Both possible biological and immune infiltration functions were investigated using Principal Component Analysis (PCA), Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and immunoassays. The sensitivity of the immune response to possible THCA medicines was assessed using ratings for tumor immune dysfunction and exclusion (TIDE) and tumor mutational burden (TMB). Results: Seven cuproptosis-related lncRNAs were used to construct our prognostic prediction model: AC108704.1, DIO3OS, AL157388.1, AL138767.3, STARD13-AS, AC008532.1, and PLBD1-AS1. Using data from TCGA's training, testing, and all groups, Kaplan-Meier and ROC curves demonstrated this feature's adequate predictive validity. Different clinical characteristics have varying effects on cuproptosis-related lncRNA risk models. Further analysis of immune cell infiltration and single sample Gene Set Enrichment Analysis (ssGSEA) supported the possibility that cuproptosis-associated lncRNAs and THCA tumor immunity were closely connected. Significantly, individuals with THCA showed a considerable decline in survival owing to the superposition effect of patients in the high-risk category and high TMB. Additionally, the low-risk group had a higher TIDE score compared with the high-risk group, indicating that these patients had suboptimal immune checkpoint blocking responses. To ensure the accuracy and reliability of our results, we further verified them using several GEO databases. Conclusion: The clinical and risk aspects of cuproptosis-related lncRNAs may aid in determining the prognosis of patients with THCA and improving therapeutic choices.

Anti-EGFR/BRAF-Tyrosine Kinase Inhibitors in Thyroid Carcinoma

Alterations in significant genes located on chromosome 7 - including epidermal growth factor receptor (EGFR) and also v-Raf murine sarcoma viral oncogene homolog B (BRAF) as a mitogen-activated protein kinase (MAPK)  - combined or not with numerical imbalances of the whole chromosome (aneuploidy-polysomy) are crucial genetic events involved in the development and progression of malignancies. Identification of EGFR/BRAF-dependent specific somatic mutations and other mechanisms of deregulation (i.e., amplification) is critical for applying targeted therapeutic approaches [tyrosine kinase inhibitors (TKIs] or monoclonal antibodies (mAbs). Thyroid carcinoma is a specific pathological entity characterized by a variety of histological sub-types. Follicular thyroid carcinoma (FTC), papillary thyroid carcinoma (PTC), medullary thyroid carcinoma (MTC), and anaplastic thyroid carcinoma (ATC) represent its main sub-types. In the current review, we explore the role of EGFR/BRAF alterations in thyroid carcinoma in conjunction with the corresponding anti-EGFR/BRAF TKI-based novel therapeutic strategies for patients with specific genetic signatures.