Using Tumor-Infiltrating Immune Cells and a ceRNA Network Model to Construct a Prognostic Analysis Model of Thyroid Carcinoma

Elucidating immunological characteristics of the adenoma-carcinoma sequence in colorectal cancer patients in South Korea using a bioinformatics approach

Colorectal cancer (CRC) is one of the top five most common and life-threatening malignancies worldwide. Most CRC develops from advanced colorectal adenoma (ACA), a precancerous stage, through the adenoma-carcinoma sequence. However, its underlying mechanisms, including how the tumor microenvironment changes, remain elusive. Therefore, we conducted an integrative analysis comparing RNA-seq data collected from 40 ACA patients who visited Dongguk University Ilsan Hospital with normal adjacent colons and tumor samples from 18 CRC patients collected from a public database. Differential expression analysis identified 21 and 79 sequentially up- or down-regulated genes across the continuum, respectively. The functional centrality of the continuum genes was assessed through network analysis, identifying 11 up- and 13 down-regulated hub-genes. Subsequently, we validated the prognostic effects of hub-genes using the Kaplan-Meier survival analysis. To estimate the immunological transition of the adenoma-carcinoma sequence, single-cell deconvolution and immune repertoire analyses were conducted. Significant composition changes for innate immunity cells and decreased plasma B-cells with immunoglobulin diversity were observed, along with distinctive immunoglobulin recombination patterns. Taken together, we believe our findings suggest underlying transcriptional and immunological changes during the adenoma-carcinoma sequence, contributing to the further development of pre-diagnostic markers for CRC.

TIM3 activates the ERK1/2 pathway to promote invasion and migration of thyroid tumors

BACKGROUND

This study aims to study the possible action mechanism of T-cell immunoglobulin and mucin domain 3 (TIM3) on the migratory and invasive abilities of thyroid carcinoma (TC) cells.

METHODS

GSE104005 and GSE138198 datasets were downloaded from the GEO database for identifying differentially expressed genes (DEGs). Functional enrichment analysis and protein-protein interaction (PPI) analysis were performed on the common DEGs in GSE104005 and GSE138198 datasets. Subsequently, in order to understand the effect of a common DEG (TIM3) on TC cells, we performed in vitro experiments using FRO cells. The migratory and invasive abilities of FRO cells were detected by wound scratch assay and Transwell assay. Proteins expression levels of the phosphorylated (p)-extracellular signal-regulated kinase (ERK)1/2, matrix metalloproteinase-2 (MMP-2) and MMP-9 were determined via Western blotting after ERK1/2 inhibition in TIM3-NC group and TIM3-mimic group.

RESULTS

316 common DEGs were identified in GSE104005 and GSE138198 datasets. These DEGs were involved in the biological process of ERK1 and ERK2 cascade. TIM3 was significantly up-regulated in TC. In vitro cell experiments showed that TIM3 could promote migration and invasion of TC cells. Moreover, TIM3 may affect the migration, invasive abilities of TC cells by activating the ERK1/2 pathway.

CONCLUSION

The above results indicate that TIM3 may affect the migratory and invasive of TC cells by activating the ERK1/2 pathway.

Tumor-infiltrating immune cell score as an independent prognostic predictor for endometrial carcinoma: Insights from a comprehensive analysis of the immune landscape

BACKGROUND

Immune cells are crucial components in the tumor microenvironment and have a significant impact on the outcomes of patients.

AIMS

Here, we aimed to establish a prognostic score based on different types of tumor-infiltrating immune cells for Endometrial Carcinoma (EC).

METHODS AND RESULTS

We enrolled and analyzed 516 EC patients from The Cancer Genome Atlas. The relative abundance of 22 immune cells were estimated by using the CIBERSORTx algorithm. Cox regression was performed to identify potential prognostic immune cells, which were used to develop a Tumor-infiltrating Immune Cell Score (TICS). The prognostic and incremental value of TICS for overall survival were compared with traditional prognostic factors using the C-index and decision curves. Clustering analysis using all immune cells identified three immune landscape subtypes, which had weak correlation with survival. A TICS was constructed using CD8T cells, resting memory CD4 T cells, activated NK and activated DCs, and classified patients as low-, moderate- and high-risk subgroups. The low-risk subgroup had higher tumor mutation burden and activation of IL2/STAT5, IL2/STAT3 and IFN-gamma response pathways. Conversely, the high-risk subgroup was associated with DNA copy number variation, hypoxia and EMT process. The TICS subgroups significantly predicted overall survival, which was independent of patient age, tumor stage, grade and molecular classification. Moreover, we developed a nomogram incorporating TICS and clinicopathologic factors, which significantly improved the predictive accuracy compared to the clinicopathologic model alone.

CONCLUSION

The TICS is an effective and independent prognostic predictor for EC patients and may serve as a useful supplement to clinicopathological factors and molecular subtyping.

New Insights into Immune Cells and Immunotherapy for Thyroid Cancer

Thyroid cancer (TC) is the most common endocrine malignancy worldwide, and the incidence of TC has gradually increased in recent decades. Differentiated thyroid cancer (DTC) is the most common subtype and has a good prognosis. However, advanced DTC patients with recurrence, metastasis and iodine refractoriness, as well as more aggressive subtypes such as poorly differentiated thyroid cancer (PDTC) and anaplastic thyroid cancer (ATC), still pose a great challenge for clinical management. Therefore, it is necessary to continue to explore the inherent molecular heterogeneity of different TC subtypes and the global landscape of the tumor immune microenvironment (TIME) to find new potential therapeutic targets. Immunotherapy is a promising therapeutic strategy that can be used alone or in combination with drugs targeting tumor-driven genes. This article focuses on the genomic characteristics, tumor-associated immune cell infiltration and immune checkpoint expression of different subtypes of TC patients to provide guidance for immunotherapy.

Machine learning-based immune prognostic model and ceRNA network construction for lung adenocarcinoma

PURPOSE

Lung adenocarcinoma (LUAD) is a malignant tumor with a high lethality rate. Immunotherapy has become a breakthrough in cancer treatment and improves patient survival and prognosis. Therefore, it is necessary to find new immune-related markers. However, the current research on immune-related markers in LUAD is not sufficient. Therefore, there is a need to find new immune-related biomarkers to help treat LUAD patients.

METHODS

In this study, a bioinformatics approach combined with a machine learning approach screened reliable immune-related markers to construct a prognostic model to predict the overall survival (OS) of LUAD patients, thus promoting the clinical application of immunotherapy in LUAD. The experimental data were obtained from The Cancer Genome Atlas (TCGA) database, including 535 LUAD and 59 healthy control samples. Firstly, the Hub gene was screened using a bioinformatics approach combined with the Support Vector Machine Recursive Feature Elimination algorithm; then, a multifactorial Cox regression analysis by constructing an immune prognostic model for LUAD and a nomogram to predict the OS rate of LUAD patients. Finally, the regulatory mechanism of Hub genes in LUAD was analyzed by ceRNA.

RESULTS

Five genes, ADM2, CDH17, DKK1, PTX3, and AC145343.1, were screened as potential immune-related genes in LUAD. Among them, ADM2 and AC145343.1 had a good prognosis in LUAD patients (HR < 1) and were novel markers. The remaining three genes screened were associated with poor prognosis in LUAD patients (HR > 1). In addition, the experimental results showed that patients in the low-risk group had better OS rates than those in the high-risk group (P < 0.001).

CONCLUSION

In this paper, we propose an immune prognostic model to predict OS rate in LUAD patients and show the correlation between five immune genes and the level of immune-related cell infiltration. It provides new markers and additional ideas for immunotherapy in patients with LUAD.

B cells in the tumor microenvironment: Multi-faceted organizers, regulators, and effectors of anti-tumor immunity

Our understanding of tumor-infiltrating lymphocytes (TILs) is rapidly expanding beyond T cell-centric perspectives to include B cells and plasma cells, collectively referred to as TIL-Bs. In many cancers, TIL-Bs carry strong prognostic significance and are emerging as key predictors of response to immune checkpoint inhibitors. TIL-Bs can perform multiple functions, including antigen presentation and antibody production, which allow them to focus immune responses on cognate antigen to support both T cell responses and innate mechanisms involving complement, macrophages, and natural killer cells. In the stroma of the most immunologically "hot" tumors, TIL-Bs are prominent components of tertiary lymphoid structures, which resemble lymph nodes structurally and functionally. Additionally, TIL-Bs participate in a variety of other lympho-myeloid aggregates and engage in dynamic interactions with the tumor stroma. Here, we summarize our current understanding of TIL-Bs in human cancer, highlighting the compelling therapeutic opportunities offered by their unique tumor recognition and effector mechanisms.

Common Genetic Factors and Pathways in Alzheimer's Disease and Ischemic Stroke: Evidences from GWAS

Alzheimer's disease (AD) and ischemic stroke (IS) are common neurological disorders, and the comorbidity of these two brain diseases is often seen. Although AD and IS were regarded as two distinct disease entities, in terms of different etiologies and clinical presentation, recent genome-wide association studies (GWASs) revealed that there were common risk genes between AD and IS, indicating common molecular pathways and their common pathophysiology. In this review, we summarize AD and IS risk single nucleotide polymorphisms (SNPs) and their representative genes from the GWAS Catalog database, and find thirteen common risk genes, but no common risk SNPs. Furthermore, the common molecular pathways associated with these risk gene products are summarized from the GeneCards database and clustered into inflammation and immunity, G protein-coupled receptor, and signal transduction. At least seven of these thirteen genes can be regulated by 23 microRNAs identified from the TargetScan database. Taken together, the imbalance of these molecular pathways may give rise to these two common brain disorders. This review sheds light on the pathogenesis of comorbidity of AD and IS, and provides molecular targets for disease prevention, manipulation, and brain health maintenance.

Overexpression of lncRNA TUG1 enhances the efficacy of DC-CIK immunotherapy in neuroblastoma in vitro and in vivo

BACKGROUND

Long non-coding RNA (LncRNA) TUG1 plays a critical role in the development of human cancers. This study explored whether TUG1 is involved in the cytotoxicity of dendritic cells and cytokine-induced killer cells (DCs-CIK), an immunotherapy approach, in neuroblastoma.

METHODS

A TUG1 expression plasmid was transfected into DCs. Neuroblastoma SK-N-SH cells were incubated with CIK cells, DCs-CIK cells, and TUG1-overexpressing DCs-CIK cells, with or without irradiation. SK-N-SH cell viability, colony formation, migration, and apoptosis were analyzed using CCK-8, colony formation assay, transwell assay, and flow cytometry, respectively. Production of IL-12, IL-2 and IFN-γ in the supernatants was determined using ELISA. A dual luciferase activity assay was performed to confirm the molecular interactions between TUG1 and miR-204. Tumor-bearing mice were established by injection of SK-N-SH cells followed by stimulation with CIK cells, DC-CIK cells, and TUG1-overexpressing DCs-CIK cells.

RESULTS

Compared to CIK alone or DC-CIK therapy, overexpression of TUG1 significantly suppressed tumor cell proliferation, colony formation, and migration of neuroblastoma cells. Moreover, upregulation of TUG1 robustly induced apoptosis and altered key molecules associated with apoptosis and epithelial-mesenchymal transition. Contents of IL-12, IL-2 and IFN-γ were dramatically elevated in the supernatants in the coculturing system upon transfection with TUG1. In addition, TUG1 was found to be act as miR-204 sponge. Furthermore, in vivo experiments demonstrated that upregulation of TUG1 potentiated the antitumor activity of DC-CIK immunotherapy.

CONCLUSION

Overexpression of TUG1 promotes DC maturation and enhances CIK cytotoxicity, suggesting that TUG1 may be a novel target for enhancing DC-CIK based immunotherapy for neuroblastoma.

The development of a novel signature based on the m6A RNA methylation regulator-related ceRNA network to predict prognosis and therapy response in sarcomas

Background: N⁶ methyladenosine (m⁶A)-related noncoding RNAs (including lncRNAs and miRNAs) are closely related to the development of cancer. However, the gene signature and prognostic value of m⁶A regulators and m⁶A-associated RNAs in regulating sarcoma (SARC) development and progression remain largely unexplored. Therefore, further research is required.

Methods: We obtained expression data for RNA sequencing (RNA-seq) and miRNAs of SARC from The Cancer Genome Atlas (TCGA) datasets. Correlation analysis and two target gene prediction databases (miRTarBase and LncBase v.2) were used to deduce m⁶A-related miRNAs and lncRNAs, and Cytoscape software was used to construct ceRNA-regulating networks. Based on univariate Cox regression and least absolute shrinkage and selection operator (LASSO) Cox regression analyses, an m⁶A-associated RNA risk signature (m6Ascore) model was established. Prognostic differences between subgroups were explored using Kaplan–Meier (KM) analysis. Risk score-related biological phenotypes were analyzed in terms of functional enrichment, tumor immune signature, and tumor mutation signature. Finally, potential immunotherapy features and drug sensitivity predictions for this model were also discussed.

Results: A total of 16 miRNAs, 104 lncRNAs, and 11 mRNAs were incorporated into the ceRNA network. The risk score was obtained based on RP11-283I3.6, hsa-miR-455-3p, and CBLL1. Patients were divided into two risk groups using the risk score, with patients in the low-risk group having longer overall survival (OS) than those in the high-risk group. The receiver operating characteristic (ROC) curves indicated that risk characteristic performed well in predicting the prognosis of patients with SARC. In addition, lower m6Ascore was also positively correlated with the abundance of immune cells such as monocytes and mast cells activated, and several immune checkpoint genes were highly expressed in the low-m6Ascore group. According to our analysis, lower m6Ascore may lead to better immunotherapy response and OS outcomes. The risk signature was significantly associated with the chemosensitivity of SARC. Finally, a nomogram was constructed to predict the OS in patients with SARC. The concordance index (C-index) for the nomogram was 0.744 (95% CI: 0.707–0.784). The decision curve analysis (DCA), calibration plot, and ROC curve all showed that this nomogram had good predictive performance.

Conclusion: This m6Ascore risk model based on m⁶A RNA methylation regulator-related RNAs may be promising for clinical prediction of prognosis and might contain potential biomarkers for treatment response prediction for SARC patients.

SuperSonic shear imaging for the differentiation between benign and malignant thyroid nodules: a meta-analysis

PURPOSE

To assess the diagnostic value of SuperSonic shear imaging (SSI) for the differentiation between benign and malignant thyroid nodules through meta-analysis.

METHODS

Online database searches were performed on PubMed, EMBASE, the Cochrane Library, and the Web of Science until 31 July 2021. The Quality Assessment of Diagnostic Accuracy Studies-2 tool was used to assess the quality of the included studies. Three measures of diagnostic test performance were used to examine the value of SSI, including the summary area under the receiver operating characteristic curve (AUROC), the summary diagnostic odds ratio (DOR), and the summary sensitivity and specificity. Heterogeneity was explored using meta-regression and subgroup analyses.

RESULTS

Finally, 21 studies with 3376 patients were included in this study. There were a total of 4296 thyroid nodules, in which 1806 malignant nodules and 2490 benign ones were involved. Thyroid nodules exhibited a malignancy rate of 42.0% (range 5.6-79.8%), 95.1% of which were of papillary variant. SSI showed a summary sensitivity of 74% [95% confidence interval (CI) 67-79%], specificity of 82% (95% CI 77-87%) and AUROC of 0.85 (95% CI 0.82-0.88) for the differentiation between benign and malignant thyroid nodules. The summary positive likelihood ratio (LR), negative LR, and DOR were 4.2 (95% CI 3.3-5.3), 0.32 (95% CI 0.26-0.40), and 13 (95% CI 9-18), respectively.

CONCLUSIONS

SSI showed high accuracy in the diagnostic differentiation between benign and malignant thyroid nodules and can be served as a noninvasive and important adjunct for thyroid nodule evaluation.

Identification of novel rheumatoid arthritis-associated MiRNA-204-5p from plasma exosomes

Rheumatoid arthritis (RA) is an autoimmune disease characterized by infiltration of immune cells in the synovium. However, the crosstalk of immune cells and synovial fibroblasts is still largely unknown. Here, global miRNA screening in plasma exosomes was carried out with a custom microarray (RA patients vs. healthy controls = 9:9). A total of 14 exosomal miRNAs were abnormally expressed in the RA patients. Then, downregulated expression of exosomal miR-204-5p was confirmed in both the replication (RA patients vs. healthy controls = 30:30) and validation groups (RA patients vs. healthy controls = 56:60). Similar to the findings obtained in humans, a decreased abundance of exosomal miR-204-5p was observed in mice with collagen-induced arthritis (CIA). Furthermore, Spearman correlation analysis indicated that plasma exosomal miR-204-5p expression was inversely correlated with disease parameters of RA patients, such as rheumatoid factor, erythrocyte sedimentation rate, and C-reactive protein. In vitro, our data showed that human T lymphocytes released exosomes containing large amounts of miR-204-5p, which can be transferred into synovial fibroblasts, inhibiting cell proliferation. Overexpression of miR-204-5p in synovial fibroblasts suppressed synovial fibroblast activation by targeting genes related to cell proliferation and invasion. In vivo assays found that administration of lentiviruses expressing miR-204-5p markedly alleviated the disease progression of the mice with CIA. Collectively, this study identified a novel RA-associated plasma exosomal miRNA-204-5p that mediates the communication between immune cells and synovial fibroblasts and can be used as a potential biomarker for RA diagnosis and treatment.

Construction of a Tumor Immune Microenvironment-Related Prognostic Model in BRAF-Mutated Papillary Thyroid Cancer

BRAF mutation is a representative oncogenic mutation, with a frequency of 60% in papillary thyroid carcinoma (PTC), but the reasons for the poor prognosis and more aggressive course of BRAF-mutated PTC are controversial. Tumor immune microenvironment (TIME) is an essential factor permitting the development and progression of malignancy, but whether TIME participates in the prognosis of BRAF-mutated PTC has not yet been reported. The primary goal of the present study was to provide a comprehensive TIME-related prognostic model to increase the predictive accuracy of progression-free survival (PFS) in patients with BRAF-mutated PTC. In this study, we analyzed the mRNA-seq data and corresponding clinical data of PTC patients obtained from the TCGA database. By calculating the TIME scores (immune score, stromal score and ESTIMATE score), the BRAF mutation group (n=237) was dichotomized into the high- and low-score groups. By functional analysis of differentially expressed genes (DEGs) in different high/low score groups, we identified 2 key TIME-related genes, HTR3A and NIPAL4, which affected PFS in BRAF-mutated PTC. A risk scoring system was developed by multivariate Cox analysis based on the abovementioned 2 TIME-related genes. Then, the BRAF-mutated cohort was divided into the high- and low-risk groups using the median risk score as a cutoff. A high risk score correlated positively with a higher HTR3A/NIPAL4 expression level but negatively with PFS in BRAF-mutated PTC. Ultimately, a nomogram was constructed by combining risk score with clinical parameter (Tumor stage), and the areas under the ROC curve (AUCs) of the nomogram for predicting 1-, 3- and 5-year PFS were then calculated and found to be 0.694, 0.707 and 0.738, respectively, indicating the improved accuracy and clinical utility of the nomogram versus the risk score model in the BRAF-mutated PTC cohort. Moreover, we determined the associations between prognostic genes or risk score and immune cell infiltration by two-way ANOVA. In the high-risk score, high HTR3A expression, and high NIPAL4 expression groups, higher infiltration of immune cells was found. Collectively, these findings confirm that the nomogram is effective in predicting the outcome of BRAF-mutated PTC and will add a spatial dimension to the developing risk stratification system.

Integrated Analysis of Circular RNA-Associated ceRNA Network Reveals Potential circRNA Biomarkers in Human Breast Cancer

Circular RNA (circRNA) is closely related to tumorigenesis and cancer progression. Yet, the roles of cancer-specific circRNAs in the circRNA-related ceRNA network of breast cancer (BRCA) remain unclear. The aim of this study was to construct a ceRNA network associated with circRNA and to explore new therapeutic and prognostic targets and biomarkers for breast cancer. We downloaded the circRNA expression profile of BRCA from Gene Expression Omnibus (GEO) microarray datasets and downloaded the miRNA and mRNA expression profiles of BRCA from The Cancer Genome Atlas (TCGA) database. Differentially expressed mRNAs (DEmRNAs), differentially expressed miRNAs (DEmiRNAs), and differentially expressed circRNAs (DEcircRNAs) were identified, and a competitive endogenous RNA (ceRNA) regulatory network was constructed based on circRNA-miRNA pairs and miRNA-mRNA pairs. Gene ontology and pathway enrichment analyses were performed on mRNAs regulated by circRNAs in ceRNA networks. Survival analysis and correlation analysis of all mRNAs and miRNAs in the ceRNA network were performed. A total of 72 DEcircRNAs, 158 DEmiRNAs, and 2762 DE mRNAs were identified. The constructed ceRNA network contains 60 circRNA-miRNA pairs and 140 miRNA-mRNA pairs, including 40 circRNAs, 30 miRNAs, and 100 mRNAs. Functional enrichment indicated that DEmRNAs regulated by DEcircRNAs in ceRNA networks were significantly enriched in the PI3K-Akt signaling pathway, microRNAs in cancer, and proteoglycans in cancer. Survival analysis and correlation analysis of all mRNAs and miRNAs in the ceRNA network showed that 13 mRNAs and 6 miRNAs were significantly associated with overall survival, and 48 miRNA-mRNA interaction pairs had a significant negative correlation. A PPI network was established, and 21 hub genes were determined from the network. This study provides an effective bioinformatics basis for further understanding of the molecular mechanisms and predictions of breast cancer. A better understanding of the circRNA-related ceRNA network in BRCA will help identify potential biomarkers for diagnosis and prognosis.

Upregulation of ADAM12 Is Associated With a Poor Survival and Immune Cell Infiltration in Colon Adenocarcinoma

BACKGROUND

A disintegrin and metalloprotease 12 (ADAM12) is a member of the multidomain protein family, but the mechanisms by which it affects prognosis and immune cell infiltration in patients with colon adenocarcinoma (COAD) remain unclear. Here, our study aimed to analyze the prognostic value of ADAM12 and investigate the correlation between ADAM12 expression and immune cell infiltration in patients with COAD.

METHODS

Differential expression analyses were performed using the Oncomine and UALCAN databases, and prognostic analyses were conducted using PrognoScan, Gene Expression Profiling Interactive Analysis (GEPIA), and Kaplan-Meier Plotter. Then, the cBioPortal database was used to analyze alterations in the ADAM12 gene, and the STRING and Metascape websites were used to conduct Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses. Additionally, relationships between ADAM12 and the immune microenvironment were evaluated based on the TIMER, GEPIA, and TISIDB databases.

RESULTS

ADAM12 was overexpressed in COAD tissues, and higher ADAM12 expression correlated with a worse prognosis for patients with COAD. The gene regulatory network suggested that ADAM12 was mainly enriched in extracellular matrix (ECM) organization, ECM proteoglycans, skeletal system development, and ossification, among other pathways. Moreover, ADAM12 expression significantly correlated with the abundance of CD4+ T cells, B cells, CD8+ T cells, neutrophils, macrophages, dendritic cells, and their markers, as well as lymphocytes, immunomodulators, and chemokines.

CONCLUSIONS

In colorectal tumors, ADAM12 may play vital roles in regulating the ECM and the recruitment of immune cells, and we suggest that ADAM12 will become a reliable biomarker for determining response to immunotherapy and the prognosis of patients with COAD.