Identification of Potential Biomarkers for Thyroid Cancer Using Bioinformatics Strategy: A Study Based on GEO Datasets

DEL-Thyroid: deep ensemble learning framework for detection of thyroid cancer progression through genomic mutation

Genes, expressed as sequences of nucleotides, are susceptible to mutations, some of which can lead to cancer. Machine learning and deep learning methods have emerged as vital tools in identifying mutations associated with cancer. Thyroid cancer ranks as the 5th most prevalent cancer in the USA, with thousands diagnosed annually. This paper presents an ensemble learning model leveraging deep learning techniques such as Long Short-Term Memory (LSTM), Gated Recurrent Units (GRUs), and Bi-directional LSTM (Bi-LSTM) to detect thyroid cancer mutations early. The model is trained on a dataset sourced from asia.ensembl.org and IntOGen.org, consisting of 633 samples with 969 mutations across 41 genes, collected from individuals of various demographics. Feature extraction encompasses techniques including Hahn moments, central moments, raw moments, and various matrix-based methods. Evaluation employs three testing methods: self-consistency test (SCT), independent set test (IST), and 10-fold cross-validation test (10-FCVT). The proposed ensemble learning model demonstrates promising performance, achieving 96% accuracy in the independent set test (IST). Statistical measures such as training accuracy, testing accuracy, recall, sensitivity, specificity, Mathew's Correlation Coefficient (MCC), loss, training accuracy, F1 Score, and Cohen's kappa are utilized for comprehensive evaluation.

Exosomal noncoding RNAs: decoding their role in thyroid cancer progression

Exosomes, as pivotal entities within the tumor microenvironment, orchestrate intercellular communication through the transfer of diverse molecules, among which non-coding RNAs (ncRNAs) such as miRNAs, lncRNAs, and circRNAs play a crucial role. These ncRNAs, endowed with regulatory functions, are selectively incorporated into exosomes. Emerging evidence underscores the significance of exosomal ncRNAs in modulating key oncogenic processes in thyroid cancer (TC), including proliferation, metastasis, epithelial-mesenchymal transition (EMT), angiogenesis, and immunoediting. The unique composition of exosomes shields their cargo from enzymatic and chemical degradation, ensuring their integrity and facilitating their specific expression in plasma. This positions exosomal ncRNAs as promising candidates for novel diagnostic and prognostic biomarkers in TC. Moreover, the potential of exosomes in the therapeutic landscape of TC is increasingly recognized. This review aims to elucidate the intricate relationship between exosomal ncRNAs and TC, fostering a deeper comprehension of their mechanistic involvement. By doing so, it endeavors to propel forward the exploration of exosomal ncRNAs in TC, ultimately paving the way for innovative diagnostic and therapeutic strategies predicated on exosomes and their ncRNA content.

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.

Role of prognostic gene DKK1 in oral squamous cell carcinoma

Oral squamous cell carcinoma (OSCC) is one of the most common squamous cell carcinomas of the head and neck. The morbidity and mortality rates of OSCC have increased in recent years. However, the pathogenesis of this disease remains unknown. The present study aimed to identify predictive biomarkers and therapeutic targets for OSCC. Bioinformatics screening of differentially expressed genes in OSCC was performed based on data from The Cancer Genome Atlas and Gene Expression Omnibus databases. Dickkopf Wnt signaling pathway inhibitor 1 (DKK1) was identified to be associated with survival, tumor immunity and DNA repair in OSCC. Furthermore, the effects of DKK1 were evaluated by the knockdown of DKK1 in two OSCC cell lines. The proliferation, clonogenicity, migration and invasion of the cells were assessed in vitro using Cell Counting Kit-8, colony formation, wound healing and Transwell assays, respectively, and were found to be inhibited by DKK1 knockdown. The present study suggests that DKK1 may be used in the prognosis of patients with OSCC and that targeting DKK1 is a potential strategy for OSCC therapy.

Epigenetic and Long-Term Effects of Nicotine on Biology, Behavior, and Health

Tobacco and nicotine use are associated with disease susceptibility and progression. Health challenges associated with nicotine and smoking include developmental delays, addiction, mental health and behavioral changes, lung disease, cardiovascular disease, endocrine disorders, diabetes, immune system changes, and cancer. Increasing evidence suggests that nicotine-associated epigenetic changes may mediate or moderate the development and progression of a myriad of negative health outcomes. In addition, nicotine exposure may confer increased lifelong susceptibility to disease and mental health challenges through alteration of epigenetic signaling. This review examines the relationship between nicotine exposure (and smoking), epigenetic changes, and maladaptive outcomes that include developmental disorders, addiction, mental health challenges, pulmonary disease, cardiovascular disease, endocrine disorders, diabetes, immune system changes, and cancer. Overall, findings support the contention that nicotine (or smoking) associated altered epigenetic signaling is a contributing factor to disease and health challenges.

Identification and validation of eight estrogen-related genes for predicting prognosis of papillary thyroid cancer

Papillary thyroid cancer (PTC) is one of the most common malignant tumors in female, and estrogen can affect its progression. However, the targets and mechanisms of estrogen action in PTC remain unclear. Therefore, this study focuses on the relationship between estrogen-related genes (ERGs) expression and prognosis in PTC, particularly neuropeptide U (NMU), and its important role in tumor progression. Based on The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases, differentially expressed genes (DEGs) predominantly enriched in ERGs were identified between PTC and normal tissue. Then, we identified ERGs that contributed most to PTC prognosis, including Transducer of ERBB2 1 (TOB1), trefoil factor 1 (TFF1), phospholipase A and acyltransferase 3 (PLAAT3), NMU, kinesin family member 20A (KIF20A), DNA topoisomerase II alpha (TOP2A), tetraspanin 13 (TSPAN13), and carboxypeptidase E (CPE). In addition, we confirmed that NMU was highly expressed in PTC and explored the effect of NMU on PTC cells proliferation in vitro and in vivo. The results showed that the proliferative capacity of PTC cells was significantly reduced with NMU knockdown. Moreover, the phosphorylation levels of the Kirsten rat sarcoma virus (KRAS) signaling pathway were significantly lower with NMU knockdown. These results suggest that ERGs, especially NMU, may be novel prognostic indicators in PTC.

Bioinformatics identify the role of chordin-like 1 in thyroid cancer

The abnormal expression of chordin-like 1 (CHRDL1) is identified in many cancers, while the effect of CHRDL1 in thyroid cancer (THCA) remains unclear. The University of California Santa Cruz, Gene Expression Profiling Interactive Analysis, University of Alabama at Birmingham Cancer, and Gene Expression Omnibus database (GSE33570, GSE33630, and GSE60542) were used for determining the mRNA and methylation expression of CHRDL1 in tumor and normal tissues. Human Protein Atlas was used for exploring the protein expression level of CHRDL1. The genes correlated to CHRDL1 were assessed by cBioPortal database. The prognostic value of CHRDL1 was evaluated through Kaplan-Meier method, cox regression, and nomogram analysis. Kyoto Encyclopedia of Genes and Genomes, Gene Ontology, and gene set enrichment analysis were used for predicting potential function of CHRDL1. The relationship between CHRDL1 and immune cell infiltration was determined by Pearson method. The downregulated mRNA and protein expressions of CHRDL1 were identified in THCA through the analysis of data from The Cancer Genome Atlas, Gene Expression Omnibus, and Human Protein Atlas database. The survival analysis showed that the CHRDL1 expression significantly affected disease-free interval (DFI) and progression-free interval, and CHRDL1 was an independent predictor of DFI. Besides, we found that C-C motif chemokine ligand 21 could significantly affect DFI time when it was co-expressed with CHRDL1. Additionally, the function of CHRDL1 was enriched in cell migration, apoptosis, and immune cell receptor. The downregulated expression of CHRDL1 was observed in THCA and caused poor prognosis. CHRDL1 may be involved in signal pathway related to cancer development and immune response, which suggested it could be a potential biomarker.

Identification of potential biomarkers for papillary thyroid carcinoma by comprehensive bioinformatics analysis

To perform bioinformatics analysis on the papillary thyroid carcinoma (PTC) gene chip dataset to explore new biological markers for PTC. The gene expression profiles of GSE3467 and GSE6004 chip data were collected by GEO2R, and the differentially expressed genes (DEGs) were selected for Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. Protein-protein interaction (PPI) relationship analysis was achieved using STRING, and the hub genes were obtained using the Cytoscape software. GEPIA was used to validate the expressions of the hub genes in the normal and tumor tissues and to conduct survival analyses. Pertinent genetic pathology results were fetched using the HPA database. Finally, the key genes were clinically verified by reverse transcription-polymerase chain reaction. 97 genes were jointly up-regulated and 107 genes were jointly down-regulated in GSE3467 and GSE6004. GO function enrichment analysis revealed that the DEGs were involved in the regulation of calcium ion transport into cytosol, integrin binding, and cell adhesion molecule binding. KEGG pathway enrichment analysis indicated that the DEGs were chiefly associated with thyroid cancer and non-small cell lung cancer. According to the PPI network, 30 key target genes were identified. Only the expressions of ANK2, TLE1, and TCF4 matched between the normal and tumor tissues, and were associated with disease prognosis. When compared with the normal thyroid tissues, the protein and mRNA expressions of ANK2, TLE1, and TCF4 were down-regulated in PTC. Significant differences exist in overall gene expression between the thyroid tissues of patients with PTC and those of healthy people. Furthermore, the differential genes ANK2, TLE1, and TCF4 are expected to be reliable molecular markers for the mechanism study and diagnosis of PTC.

Integrative Analysis of Potential Biomarkers Involved in the Progression of Papillary Thyroid Cancer

This study aims to explore key prognostic and diagnostic biomarkers involved in the pathogenesis of papillary thyroid cancer (PTC) which is one of the most common endocrine cancers and whose occurrence is rapidly increasing. Papillary thyroid cancer datasets containing normal and tumor samples were collected from Gene Expression Omnibus. Protein-protein interaction (PPI) network for common upregulated differentially expressed genes (DEGs) was constructed, and hub genes were studied. Gene ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis were performed to identify the vital biological behaviors and pathways involved in PTC. PPI network analysis demonstrated the interaction between 134 common upregulated DEGs, and top 15 pivotal genes with highest degree of connectivity were retrieved. Three of the hub genes (DPP4, ITGA2, FN1) were linked to the prognosis of PTC patients and considered clinically relevant core genes via survival analysis. We suggest that the identification of key genes associated with PTC development help us in understanding molecular mechanisms related to disease. These genes could also be considered the diagnostic biomarkers or as therapeutic targets in the future treatment for PTC.

Evolution of intra-tumoral heterogeneity across different pathological stages in papillary thyroid carcinoma

BACKGROUND

Intra-tumor heterogeneity (ITH) results from the continuous accumulation of mutations during disease progression, thus impacting patients' clinical outcome. How the ITH evolves across papillary thyroid carcinoma (PTC) different tumor stages is lacking.

METHODS

We used the whole-exome sequencing data from The Cancer Genome Atlas Thyroid Cancer (TCGA-THCA) cohort to track the ITH and assessed its relationship with clinical features through different stages of the PTC progression. We further assayed the expression levels of the specific genes in papillary thyroid cancer cell lines compared to an immortalized normal thyroid epithelial cell line by qRT-PCR.

RESULTS

We revealed the timing of mutational processes and the dynamics of the temporal acquisition of somatic events during the lifetime of the PTC. ITH significantly influences the PTC patient's survival rate and, as genetic heterogeneity increases, the prognosis gets worse in advanced tumor stages. ITH also affects the mutational architecture of each clinical stage which is subject to periodic fluctuations. Different mutational processes may cooperate to shape a stage-specific mutational spectrum during the progression from early to advanced tumor stages. Moreover, different evolutionary paths characterize PTC progression across pathological stages due to both mutations recurrently occurring in all stages in hotspot positions and distinct codon changes dominating in different stages. A different expression level of specific genes also exists in different thyroid cancer cell lines.

CONCLUSIONS

Our findings suggest ITH as a potential unfavorable prognostic factor in PTC and highlight the dynamic changes in different clinical stages of PTC, providing some clues for the precision medicine and suggesting different diagnostic decisions depending on the clinical stages of patients. Finally, complete clear guidelines to define risk stratification of PTC patients are lacking; thus, this work could contribute to defining patients who need more aggressive treatments and, in turn, could reduce the social burden of this cancer.

Chordin-Like 1 Regulates Epithelial-to-Mesenchymal Transition and Metastasis via the MAPK Signaling Pathway in Oral Squamous Cell Carcinoma

Background

Accumulating evidence suggests that dysregulation of Chordin-like 1 (CHRDL1) is associated with malignant biological behaviors in multiple cancers. However, the exact function and molecular mechanism of CHRDL1 in oral squamous cell carcinoma (OSCC) remain unclear.

Methods

The expression levels of CHRDL1 in OSCC tissues and CAL27 cells were determined by RT-qPCR. Immunohistochemical staining was applied to detect CHRDL1 protein expression in sample tissues from OSCC patients. Gain of function and knockdown by lentivirus were further used to examine the effects of CHRDL1 on cell proliferation, migration, invasion, and adhesion in OSCC. Tail vein injection of CAL27 cells with dysregulated CHRDL1 expression was further used to examine the effect of CHRDL1 on lung colonization. RNA sequencing was performed to explore the molecular mechanisms of CHRDL1 that underlie the progression of OSCC.

Results

CHRDL1 was significantly downregulated in OSCC tissues and CAL27 cells compared to controls. CHRDL1 knockdown enhanced migration, invasion, adhesion, and EMT, but not proliferation, in CAL27 cells. Overexpression of CHRDL1 had the opposite effects. Moreover, CHRDL1 was proven to inhibit tumor metastasis in vivo. Mechanistically, MAPK signaling pathway components, including ERK1/2, p38, and JNK, were found to regulate the malignant biological behaviors of CAL27 cells.

Conclusions

Our results suggest that CHRDL1 has an inhibitory effect on OSCC metastasis via the MAPK signaling pathway, which provides a new possible potential therapeutic target against OSCC.

TMEM139 prevents NSCLC metastasis by inhibiting lysosomal degradation of E-cadherin

Non‐small‐cell lung cancer (NSCLC) accounts for approximately 85% of all lung cancer cases and has the highest mortality rate among all solid tumors. It is characterized by early metastasis, and investigations of the molecular mechanisms underlying the progression and metastasis of NSCLC are urgently needed for the development of therapeutic targets. Here, we report that the transmembrane protein TMEM139 is significantly downregulated in NSCLC and that reduced expression of TMEM139 is correlated with a poor prognosis in NSCLC patients. Mechanistically, we found that TMEM139 directly interacts with E‐cadherin at the plasma membrane and at focal adhesion sites. Moreover, TMEM139 can prevent the lysosomal degradation of E‐cadherin, which inhibits epithelial‐mesenchymal transition, migration and invasion of NSCLC cells both in vitro and in vivo. Our study not only expands our understanding of NSCLC metastasis but also provides a foundation to develop novel therapeutic strategies.

Identification of a Four-Gene Signature for Determining the Prognosis of Papillary Thyroid Carcinoma by Integrated Bioinformatics Analysis

Purpose

Although well-differentiated papillary thyroid carcinoma (PTC) has an indolent nature and usually an excellent prognosis, some patients experience disease recurrence or death. The aim of this study was to identify prognostic markers to stratify PTC patients.

Patients and Methods

Eight gene-expression profiles (GSE3467, GSE3678, GSE5364, GSE27155, GSE33630, GSE53157, GSE60542, and GSE104005) were obtained from the Gene Expression Omnibus and used to analyze differentially expressed genes (DEGs) between PTC tissues and non-tumor tissues. Univariable Cox regression survival analysis and Lasso-penalized Cox regression analysis were performed to identify prognostic genes and establish a risk-score model based on the integrated DEGs. Kaplan–Meier (KM) and receiver operating characteristic (ROC) curves were used to validate the prognostic performance of the risk score. A nomogram was constructed based on The Cancer Genome Atlas dataset and Multivariable Cox regression analysis.

Results

A total of 165 upregulated and 207 downregulated DEGs were screened. A four-gene signature including PAPSS2, PCOLCE2, PTX3, and TGFBR3 was identified. The risk-score model showed a strong diagnosis performance for identifying patients with a poor prognosis. KM analysis showed that patients with low risk scores had a significantly more favorable overall survival (OS) than those with high risk scores (p = 0.0002). ROC curves based on the four-gene signature showed better performances in predicting 1-, 3-, and 5-year survival than did the American Joint Committee on Cancer staging system (area under the curve: 0.86 vs 0.84, 0.80 vs 0.63, and 0.79 vs 0.73, respectively). Furthermore, when combined with age and tumor status from the nomogram, the four-gene signature achieved a good performance in guiding postoperative follow-up surveillance of patients with PTC.

Conclusion

The four-gene signature was found to be a novel and reliable biomarker with great potential for clinical application in risk stratification and OS prediction in patients with PTC.

Chordin-like 1 is a novel prognostic biomarker and correlative with immune cell infiltration in lung adenocarcinoma

Chordin-like 1 (CHRDL1), an inhibitor of bone morphogenetic proteins(BMPs), has been recently reported to participate in the progression of numerous tumors, however, its role in lung adenocarcinoma (LUAD) remains unclear. Our study aimed to demonstrate relationship between CHRDL1 and LUAD based on data from The Cancer Genome Atlas (TCGA). Among them, CHRDL1 expression revealed promising power for distinguishing LUAD tissues form normal sample. Low CHRDL1 was correlated with poor clinicopathologic features, including high T stage (OR=0.45, P<0.001), high N stage (OR=0.57, P<0.003), bad treatment effect (OR=0.64, P=0.047), positive tumor status (OR=0.63, P=0.018), and TP53 mutation (OR=0.49, P<0.001). The survival curve illustrated that low CHRDL1 was significantly correlative with a poor overall survival (HR=0.60, P<0.001). At multivariate Cox regression analysis, CHRDL1 remained independently correlative with overall survival. GSEA identified that the CHRDL1 expression was related to cell cycle and immunoregulation. Immune infiltration analysis suggested that CHRDL1 was significantly correlative with 7 kinds of immune cells. Immunohistochemical validation showed that CHRDL1 was abnormally elevated and negatively correlated with Th2 cells in LUAD tissues. In conclusion, CHRDL1 might become a novel prognostic biomarker and therapy target in LUAD. Moreover, CHRDL1 may improve the effectiveness of immunotherapy by regulating immune infiltration.

Noncoding RNAs in Papillary Thyroid Cancer: Interaction with Cancer-Associated Fibroblasts (CAFs) in the Tumor Microenvironment (TME) and Regulators of Differentiation and Lymph Node Metastasis

A large majority of all thyroid cancers are papillary thyroid carcinomas (PTC), named for the specific papillary architecture observed histologically. Despite the high rate of success with modern diagnostic and therapeutic algorithms, there are significant areas where the management of PTC can be improved. Aggressive PTC subtypes that are refractory to radioactive iodine (RAI) therapy carry a more severe prognosis and account for most of PTC-related deaths. As lymph node metastasis is present in roughly 40% of all adult PTC cases, higher specificity in these tests is a clinical need, especially since lymph node metastases are associated with reduced survival and higher recurrence rates. Additionally, this cancer can progress to more dedifferentiated and aggressive variants, such as poorly differentiated papillary thyroid cancer (PDPTC) and anaplastic thyroid cancer (ATC). Therefore, development of more sensitive and specific detection methods that allow unnecessary surgeries to be avoided is of the utmost importance. The body of large-scale, unbiased gene expression analysis in PTC has focused on the coding transcriptome, specifically mRNAs and microRNAs. However, there have been implications for the potential use of long noncoding RNAs (lncRNAs) in PTC diagnosis, prognosis, and treatment via the utilization of genome-wide studies of patient samples. lncRNAs have diverse regulatory potential in gene expression, alternative splicing, posttranscriptional mRNA modification, and epigenomic alterations. Many lncRNAs have tissue-specific expression and are demonstrated to play key roles in cancer progression and prognosis. However, lncRNAs are not being exploited as biomarkers or therapeutic targets currently, despite their elucidated effects on oncogenesis. These potent biomarkers would be revolutionary in detection at early stages, as this significantly increases the chances of survival. Their aberrant expression in cancer and correlation with steps in tumorigenesis as well as their role in differentiation would allow for a promising role as a prognostic and diagnostic biomarker in thyroid cancer. This would help prevent the more aggressive ATC that derives from dedifferentiation of the less aggressive PTC and FTC. The targeting of the specific lncRNAs could also pose a valuable treatment option via preventing or reversing this dedifferentiation process and making this usually refractory form of thyroid cancer more responsive to standard treatment options.

Prognostic biomarkers for predicting papillary thyroid carcinoma patients at high risk using nine genes of apoptotic pathway

Aberrant expressions of apoptotic genes have been associated with papillary thyroid carcinoma (PTC) in the past, however, their prognostic role and utility as biomarkers remains poorly understood. In this study, we analysed 505 PTC patients by employing Cox-PH regression techniques, prognostic index models and machine learning methods to elucidate the relationship between overall survival (OS) of PTC patients and 165 apoptosis related genes. It was observed that nine genes (ANXA1, TGFBR3, CLU, PSEN1, TNFRSF12A, GPX4, TIMP3, LEF1, BNIP3L) showed significant association with OS of PTC patients. Five out of nine genes were found to be positively correlated with OS of the patients, while the remaining four genes were negatively correlated. These genes were used for developing risk prediction models, which can be utilized to classify patients with a higher risk of death from the patients which have a good prognosis. Our voting-based model achieved highest performance (HR = 41.59, p = 3.36x10^-4, C = 0.84, logrank-p = 3.8x10^-8). The performance of voting-based model improved significantly when we used the age of patients with prognostic biomarker genes and achieved HR = 57.04 with p = 10⁻⁴ (C = 0.88, logrank-p = 1.44x10^-9). We also developed classification models that can classify high risk patients (survival ≤ 6 years) and low risk patients (survival > 6 years). Our best model achieved AUROC of 0.92. Further, the expression pattern of the prognostic genes was verified at mRNA level, which showed their differential expression between normal and PTC samples. Also, the immunostaining results from HPA validated these findings. Since these genes can also be used as potential therapeutic targets in PTC, we also identified potential drug molecules which could modulate their expression profile. The study briefly revealed the key prognostic biomarker genes in the apoptotic pathway whose altered expression is associated with PTC progression and aggressiveness. In addition to this, risk assessment models proposed here can help in efficient management of PTC patients.

Development and validation of a novel 14-gene signature for predicting lymph node metastasis in papillary thyroid carcinoma

Background

There is still no reasonably accurate method of preoperatively predicting central lymph node metastasis (LNM), and it is essential to develop an effective evaluation model for predicting LNM in papillary thyroid carcinoma (PTC) patients.

Methods

PTC samples were collected from The Cancer Genome Atlas database. Candidate genes were identified as continuously upregulated or downregulated genes in the process of N0 to N1a and N1a to N1b. The least absolute shrinkage and selection operator (LASSO) regression analysis was used to construct the predictive model for LNM. Multivariate logistic regression analysis was performed to screen the potential factors related to LNM, and a nomogram was established. The risk score of the gene signature model for predicting disease-free survival (DFS) was evaluated by Kaplan-Meier analysis.

Results

A 14-gene signature was developed by LASSO regression for predicting LNM based on 69 differential expression genes (DEGs) that were continuously upregulated or downregulated in the progress of PTC. The receiver operating characteristic (ROC) curves of the 14-gene signature predicting LNM, central LNM and lateral LNM were generated. The area under the ROC (AUC) values were 0.806 [95% confidence interval (CI): 0.7608-0.8815], 0.755 (95% CI: 0.6839-0.8263) and 0.821 (95% CI: 0.7608-0.8815). The nomogram's C-index value, including the 14-gene signature and other potential risk factors, was 0.786 (95% CI: 0.7296-0.8425), and the calibration exhibited fairly good consistency with the perfect prediction. Based on the 14-gene risk score, high-risk PTC patients had a worse DFS.

Conclusions

A novel 14-gene signature was developed for predicting LNM in PTC patients. The risk score also correlated with DFS in PTC patients.

ASMdb: a comprehensive database for allele-specific DNA methylation in diverse organisms

DNA methylation is known to be the most stable epigenetic modification and has been extensively studied in relation to cell differentiation, development, X chromosome inactivation and disease. Allele-specific DNA methylation (ASM) is a well-established mechanism for genomic imprinting and regulates imprinted gene expression. Previous studies have confirmed that certain special regions with ASM are susceptible and closely related to human carcinogenesis and plant development. In addition, recent studies have proven ASM to be an effective tumour marker. However, research on the functions of ASM in diseases and development is still extremely scarce. Here, we collected 4400 BS-Seq datasets and 1598 corresponding RNA-Seq datasets from 47 species, including human and mouse, to establish a comprehensive ASM database. We obtained the data on DNA methylation level, ASM and allele-specific expressed genes (ASEGs) and further analysed the ASM/ASEG distribution patterns of these species. In-depth ASM distribution analysis and differential methylation analysis conducted in nine cancer types showed results consistent with the reported changes in ASM in key tumour genes and revealed several potential ASM tumour-related genes. Finally, integrating these results, we constructed the first well-resourced and comprehensive ASM database for 47 species (ASMdb, www.dna-asmdb.com).

Comprehensive Pan-Cancer Analysis and the Regulatory Mechanism of ASF1B, a Gene Associated With Thyroid Cancer Prognosis in the Tumor Micro-Environment

Background

The incidence of thyroid cancer, whose local recurrence and metastasis lead to death, has always been high and the pathogenesis of papillary thyroid carcinoma (PTC) has not been clearly elucidated. Therefore, the research for more accurate prognosis-related predictive biomarkers is imminent, and a key gene can often be a prognostic marker for multiple tumors.

Methods

Gene expression profiles of various cancers in the TCGA and GTEx databases were downloaded, and genes significantly associated with the prognosis of THCA were identified by combining differential analysis with survival analysis. Then, a series of bioinformatics tools and methods were used to analyze the expression of the gene in each cancer and the correlation of each expression with prognosis, tumor immune microenvironment, immune neoantigens, immune checkpoints, DNA repair genes, and methyltransferases respectively. The possible biological mechanisms were also investigated by GSEA enrichment analysis.

Results

656 differentially expressed genes were identified from two datasets and 960 DEGs that were associated with disease-free survival in THCA patients were screened via survival analysis. The former and the latter were crossed to obtain 7 key genes, and the gene with the highest risk factor, ASF1B, was selected for this study. Differential analysis of multiple databases showed that ASF1B was commonly and highly expressed in pan-cancer. Survival analysis showed that high ASF1B expression was significantly associated with poor patient prognosis in multiple cancers. In addition, ASF1B expression levels were found to be associated with tumor immune infiltration in THCA, KIRC, LGG, and LIHC, and with tumor microenvironment in BRCA, LUSC, STAD, UCEC, and KIRC. Further analysis of the relationship between ASF1B expression and immune checker gene expression suggested that ASF1B may regulate tumor immune patterns in most tumors by regulating the expression levels of specific immune checker genes. Finally, GSEA enrichment analysis showed that ASF1B high expression was mainly enriched in cell cycle, MTORC1 signaling system, E2F targets, and G2M checkpoints pathways.

Conclusions

ASF1B may be an independent prognostic marker for predicting the prognosis of THCA patients. The pan-cancer analysis suggested that ASF1B may play an important role in the tumor micro-environment and tumor immunity and it has the potential of serving as a predictive biomarker for multiple cancers.

Long Non-coding RNA RP11-395G23.3 Acts as a Competing Endogenous RNA of miR-124-3p to Regulate ROR1 in Anaplastic Thyroid Carcinoma

Anaplastic thyroid carcinoma (ATC) is one of the most aggressive human malignancies with poor prognosis. However, the underlying mechanisms of ATC remain to be elucidated. Recently, increasing studies have focused on competitive endogenous RNA (ceRNA) to discover valuable biomarkers for the diagnosis of ATC. The present study identified 705 differentially expressed mRNAs and 47 differentially expressed lncRNAs. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were also conducted. Additionally, an lncRNA/miRNA/mRNA network was constructed which included 1103 regulatory relations. The upregulation of RP11-395G23.3 in ATC cells was confirmed by quantitative reverse transcription polymerase chain reaction (qRT-PCR). In the loss of function assays, results suggested silencing of RP11-395G23.3 inhibited cell proliferation and induced cell apoptosis. Mechanically, RP11-395G23.3 could increase ROR1 via sponging miR-124-3p as a ceRNA. Moreover, ROR1 expression was decreased with the downregulation of RP11-395G23.3, but was rescued by the co-transfection of the miR-124-3p inhibitor in ATC cells. Our research suggested that the RP11-395G23.3/miR-124-3p/ROR1 axis potentially acted as a potential target for the diagnosis of ATC.

Identification of Hub Genes Associated with the Pathogenesis of Intracranial Aneurysm via Integrated Bioinformatics Analysis

BACKGROUND

At present, the pathogenesis of intracranial aneurysms (IA) remains unclear, which significantly hinders the development of novel strategies for the clinical treatment. In this study, bioinformatics methods were used to identify the potential hub genes and pathways associated with the pathogenesis of IA.

METHODS

The gene expression datasets of patients with intracranial aneurysm were downloaded from the Gene Expression Database (GEO), and the different data sets were integrated by the robust rank aggregation (RRA) method to identify the differentially expressed genes between patients with intracranial aneurysm and the controls. The functional enrichment analyses of the significant differentially expressed genes (DEGs) were performed and the protein-protein interaction (PPI) network was constructed; thereafter, the hub genes were screened by cytoHubba plug-in of Cytoscape, and finally sequencing dataset GSE122897 was used to verify the hub genes.

RESULTS

The GSE15629, GSE75436, GSE26969, and GSE6551 expression profiles have been included in this study, including 34 intracranial aneurysm samples and 26 control samples. The four datasets obtained 136 significant DEGs (45 up-regulated, 91 down-regulated). Enrichment analysis showed that the extracellular matrix structural constituent and the ECM-receptor interaction were closely related to the occurrence of IA. It was finally determined that eight hub genes associated with the development of IA, including VCAN, COL1A1, COL11A1, COL5A1, COL5A2, POSTN, THBS2, and CDH2.

CONCLUSION

The discovery of potential hub genes and pathways could enhance the understanding of the molecular mechanisms associated with the development of IA. These hub genes may be potential therapeutic targets for the management and new biomarker for the diagnosis of IA.

Prognostic Biomarker-Based Identification of Drugs for Managing the Treatment of Endometrial Cancer

INTRODUCTION

Uterine corpus endometrial carcinoma (UCEC) causes thousands of deaths per year. To improve the overall survival of patients with UCEC, there is a need to identify prognostic biomarkers and potential drugs.

OBJECTIVES

The aim of this study was twofold: the identification of prognostic gene signatures from expression profiles of pattern recognition receptor (PRR) genes and identification of the most effective existing drugs using the prognostic gene signature.

METHODS

This study was based on the expression profile of PRR genes of 541 patients with UCEC obtained from The Cancer Genome Atlas. Key prognostic signatures were identified using various approaches, including survival analysis, network, and clustering. Hub genes were identified by constructing a co-expression network. Representative genes were identified using k-means and k-medoids-based clustering. Univariate Cox proportional hazard (PH) analysis was used to identify survival-associated genes. 'cmap2' was used to identify potential drugs that can suppress/enhance the expression of prognostic genes.

RESULTS

Models were developed using hub genes and achieved a maximum hazard ratio (HR) of 1.37 (p = 0.294). Then, a clustering-based model was developed using seven genes (HR 9.14; p = 1.49 × 10^-12). Finally, a nine gene-based risk stratification model was developed (CLEC1B, CLEC3A, IRF7, CTSB, FCN1, RIPK2, NLRP10, NLRP9, and SARM1) and achieved HR 10.70; p = 1.1 × 10^-12. The performance of this model improved significantly in combination with the clinical stage and achieved HR 15.23; p = 2.21 × 10^-7. We also developed a model for predicting high-risk patients (survival ≤ 4.3 years) and achieved an area under the receiver operating characteristic curve (AUROC) of 0.86.

CONCLUSION

We identified potential immunotherapeutic agents based on prognostic gene signature: hexamethonium bromide and isoflupredone. Several novel candidate drugs were suggested, including human interferon-α-2b, paclitaxel, imiquimod, MESO-DAP1, and mifamurtide. These biomolecules and repurposed drugs may be utilised for prognosis and treatment for better survival.

Low expression of CHRDL1 and SPARCL1 predicts poor prognosis of lung adenocarcinoma based on comprehensive analysis and immunohistochemical validation

Purpose

Exploring the molecular mechanisms of lung adenocarcinoma (LUAD) is beneficial for developing new therapeutic strategies and predicting prognosis. This study was performed to select core genes related to LUAD and to analyze their prognostic value.

Methods

Microarray datasets from the GEO (GSE75037) and TCGA-LUAD datasets were analyzed to identify differentially coexpressed genes in LUAD using weighted gene coexpression network analysis (WGCNA) and differential gene expression analysis. Functional enrichment analysis was conducted, and a protein–protein interaction (PPI) network was established. Subsequently, hub genes were identified using the CytoHubba plug-in. Overall survival (OS) analyses of hub genes were performed. The Clinical Proteomic Tumor Analysis Consortium (CPTAC) and the Human Protein Atlas (THPA) databases were used to validate our findings. Gene set enrichment analysis (GSEA) of survival-related hub genes were conducted. Immunohistochemistry (IHC) was carried out to validate our findings.

Results

We identified 486 differentially coexpressed genes. Functional enrichment analysis suggested these genes were primarily enriched in the regulation of epithelial cell proliferation, collagen-containing extracellular matrix, transforming growth factor beta binding, and signaling pathways regulating the pluripotency of stem cells. Ten hub genes were detected using the maximal clique centrality (MCC) algorithm, and four genes were closely associated with OS. The CPTAC and THPA databases revealed that CHRDL1 and SPARCL1 were downregulated at the mRNA and protein expression levels in LUAD, whereas SPP1 was upregulated. GSEA demonstrated that DNA-dependent DNA replication and catalytic activity acting on RNA were correlated with CHRDL1 and SPARCL1 expression, respectively. The IHC results suggested that CHRDL1 and SPARCL1 were significantly downregulated in LUAD.

Conclusions

Our study revealed that survival-related hub genes closely correlated with the initiation and progression of LUAD. Furthermore, CHRDL1 and SPARCL1 are potential therapeutic and prognostic indicators of LUAD.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12935-021-01933-9.

Identification and analysis of key genes associated with acute myocardial infarction by integrated bioinformatics methods

BACKGROUND

Acute myocardial infarction (AMI) is a common disease leading threat to human health around the world. Here we aimed to explore new biomarkers and potential therapeutic targets in AMI through adopting integrated bioinformatics tools.

METHODS

The gene expression Omnibus (GEO) database was used to obtain genes data of AMI and no-AMI whole blood. Furthermore, differentially expressed genes (DEGs) were screened using the "Limma" package in R 3.6.1 software. Functional and pathway enrichment analyses of DEGs were performed via "Bioconductor" and "GOplot" package in R 3.6.1 software. In order to screen hub DEGs, the STRING version 11.0 database, Cytoscape and molecular complex detection (MCODE) were applied. Correlation among the hub DEGs was evaluated using Pearson's correlation analysis.

RESULTS

By performing DEGs analysis, 289 upregulated and 62 downregulated DEGs were successfully identified from GSE66360, respectively. And they were mainly enriched in the terms of neutrophil activation, immune response, cytokine, nuclear factor kappa-B (NF-κB) signaling pathway, IL-17 signaling pathway, and tumor necrosis factor (TNF) signaling pathway. Based on the data of protein-protein interaction (PPI), the top 10 hub genes were ranked, including interleukin-8 (CXCL8), TNF, N-formyl peptide receptor 2 (FPR2), growth-regulated alpha protein (CXCL1), transcription factor AP-1 (JUN), interleukin-1 beta (IL1B), platelet basic protein (PPBP), matrix metalloproteinase-9 (MMP9), toll-like receptor 2 (TLR2), and high affinity immunoglobulin epsilon receptor subunit gamma (FCER1G). What's more, the results of correlation analysis demonstrated that there was positive correlation between the 10 hub DEGs.

CONCLUSION

Ten DEGs were identified as potential candidate diagnostic biomarkers for patients with AMI in present study. However, further experiments are needed to confirm the functional pathways and hub genes associated with AMI.

A Prognostic Nomogram Model Based on mRNA Expression of DNA Methylation-Driven Genes for Gastric Cancer

PURPOSE

The exploration and interpretation of DNA methylation-driven genes might contribute to molecular classification, prognostic prediction and therapeutic choice. In this study, we built a prognostic risk model via integrating analysis of the transcriptome and methylation profile for patients with gastric cancer (GC).

METHODS

The mRNA expression profiles, DNA methylation profiles and corresponding clinicopathological information of 415 GC patients were downloaded from The Cancer Genome Atlas (TCGA). Differential expression and correlation analysis were performed to identify DNA methylation-driven genes. The candidate genes were selected by univariate Cox regression analyses followed by the least absolute shrinkage and selection operator (LASSO) regression. A prognostic risk nomogram model was then built together with clinicopathological parameters.

RESULTS

5 DNA methylation-driven genes (CXCL3, F5, GNAI1, GAMT and GHR) were identified by integrated analyses and selected to construct the prognostic risk model with clinicopathological parameters. High expression and low DNA hypermethylation of F5, GNAI1, GAMT and GHR, as well as low expression and high DNA hypomethylation of CXCL3 were significantly associated with poor prognosis rates, respectively. The high-risk group showed a significantly shorter prognosis than the low-risk group in the TCGA dataset (HR = 0.212, 95% CI = 0.139-0.322, P = 2e-15). The final nomogram model showed high predictive efficiency and consistency in the training and validation group.

CONCLUSION

We construct and validate a prognostic nomogram model for GC based on five DNA methylation-driven genes with high performance and stability. This nomogram model might be a powerful tool for prognosis evaluation in the clinic and also provided novel insights into the epigenetics in GC.