Focal adhesion signaling and therapy resistance in cancer

Identification of an m6A-Related Long Noncoding RNA Risk Model for Predicting Prognosis and Directing Treatments in Patients With Colon Adenocarcinoma

N6-methyladenosine (m6A) and lncRNAs have been implicated in the development of colon cancer, including tumorigenesis, migration, and invasion. However, the specific effect of m6A regulators on lncRNAs is not clear, and m6A-related lncRNAs may be new prognostic biomarkers and may help direct treatment and medication. We identified 29 prognostic m6A-related lncRNAs and constructed a risk model using 12 lncRNAs. The model was an independent prognostic factor and could accurately predict the prognosis. A stable and robust nomogram that combined the model and pathologic stage was constructed. A total of 2,424 differentially expressed genes (DEGs) were identified based on the model. Functional analysis of the DEGs showed that they were associated with tumor progression, helping investigate the underlying biological functions and signaling pathways of the risk model. In addition, the low-risk group based on the risk model had more sensitivity to afatinib, metformin, and GW.441756, and patients with low risk would more likely respond to immunotherapy. Moreover, patients with higher risk were more sensitive to olaparib, bexarotene, and doxorubicin.

Integrin β1 in Pancreatic Cancer: Expressions, Functions, and Clinical Implications

Simple Summary

Pancreatic cancer (PC) is a highly aggressive malignant tumor with an extremely poor prognosis. Early diagnosis and treatment are key to improving the survival rate of PC patients. Emerging studies show that integrins might contribute to the pathogenesis of PC. This review presents the various signaling pathways that are mediated by integrins in PC and emphasizes the multiple functions of integrin β1 in malignant behaviors of PC. It also discusses the clinical significance of integrin β1 as well as integrin β1-based therapy in PC patients.

Abstract

Pancreatic cancer (PC) is characterized by rapid progression and a high mortality rate. The current treatment is still based on surgical treatment, supplemented by radiotherapy and chemotherapy, and new methods of combining immune and molecular biological treatments are being explored. Despite this, the survival rate of PC patients is still very disappointing. Therefore, clarifying the molecular mechanism of PC pathogenesis and developing precisely targeted drugs are key to improving PC prognosis. As the most common β subunit of the integrin family, integrin β1 has been proved to be closely related to the vascular invasion, distant metastasis, and survival of PC patients, and treatment targeting integrin β1 in PC has gained initial success in animal models. In this review, we summarize the various signaling pathways by which integrins are involved in PC, focusing on the roles of integrin β1 in the malignant behaviors of PC. Additionally, recent studies regarding the feasibility of integrin β1 as a diagnostic and prognostic biomarker in PC are also discussed. Finally, we present the progress of several integrin β1-based clinical trials to highlight the potential of integrin β1 as a target for personalized therapy in PC.

A nanomechanical strategy involving focal adhesion kinase for overcoming drug resistance in breast cancer

Despite implementation of nanomechanical studies in cancer research, studies on the nanomechanical aspects of drug resistance in cancer are lacking. Here, we established the mechanical signatures of drug-resistant breast cancer cells using atomic force microscopy-based indentation techniques and functionalized nanopatterned substrates (NPS). Additionally, we examined the expression of proteins pertinent to focal adhesions in order to elucidate the molecular signatures responsible for the acquisition of drug resistance in breast cancer cells. Drug-resistant breast cancer cells exhibited mechanical reinforcement, increased actin stress fibers, dysfunctional mechano-reciprocal interaction with the NPS, vinculin overexpression, and improved focal adhesion kinase (FAK) activity. Owing to differences in FAK activation upon co-treatment with a FAK inhibitor, the drug-resistant breast cancer cells were eradicated more efficiently than invasive breast cancer cells having pro-survival activity. These findings demonstrated the potential of a novel co-treatment regimen using FAK inhibitors for overcoming drug resistance in breast cancer cells.

Identification of Feature Genes of a Novel Neural Network Model for Bladder Cancer

Background: The combination of deep learning methods and oncogenomics can provide an effective diagnostic method for malignant tumors; thus, we attempted to construct a reliable artificial neural network model as a novel diagnostic tool for Bladder cancer (BLCA). Methods: Three expression profiling datasets (GSE61615, GSE65635, and GSE100926) were downloaded from the Gene Expression Omnibus (GEO) database. GSE61615 and GSE65635 were taken as the train group, while GSE100926 was set as the test group. Differentially expressed genes (DEGs) were filtered out based on the logFC and FDR values. We also performed Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses to explore the biological functions of the DEGs. Consequently, we utilized a random forest algorithm to identify feature genes and further constructed a neural network model. The test group was given the same procedures to validate the reliability of the model. We also explored immune cells' infiltration degree and correlation coefficients through the CiberSort algorithm and corrplot R package. The qRT-PCR assay was implemented to examine the expression level of the feature genes in vitro. Results: A total of 265 DEGs were filtered out and significantly enriched in muscle system processes, collagen-containing and focal adhesion signaling pathways. Based on the random forest algorithm, we selected 14 feature genes to construct the neural network model. The area under the curve (AUC) of the training group was 0.950 (95% CI: 0.850-1.000), and the AUC of the test group was 0.667 (95% CI: 0.333-1.000). Besides, we observed significant differences in the content of immune infiltrating cells and the expression levels of the feature genes. Conclusion: After repeated verification, our neural network model had clinical feasibility to identify bladder cancer patients and provided a potential target to improve the management of BLCA.

ANXA1: An Important Independent Prognostic Factor and Molecular Target in Glioma

Objective: The expression, prognosis, and related mechanisms of ANXA1 are investigated in glioma, with the objective to find potential therapeutic molecular targets for glioma. Methods: We analyzed the gene expression of ANXA1 using glioma-related databases, including the Chinese Glioma Genome Atlas (CGGA) database, The Cancer Genome Atlas (TCGA) database, and the Gene Expression Omnibus (GEO) database. Moreover, we collected the sample tissues and corresponding paracancerous tissues of 23 glioma patients and then conducted a Western blot experiment to verify the expression and correlate survival of ANXA1. Moreover, we generated survival ROC curves, performing univariate and multivariate Cox analyses and the construction of the nomogram. Differential expression analysis was conducted by high and low grouping based on the median of the ANXA1 gene expression values. We conducted Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis and Gene Set Enrichment Analysis (GSEA) to explore possible mechanisms, and gene co-expression analysis was also performed. Results: The results showed that the ANXA1 expression level was higher in gliomas than in normal tissues, and a high expression level of ANXA1 in gliomas was associated with poorer prognosis. The independent prognosis analysis showed that the ANXA1 gene was an independent prognostic factor of glioma. In the analysis of KEGG and Gene Set Enrichment Analysis (GSEA), it is shown that ANXA1 may play an important role in glioma patients by affecting extracellular matrix (ECM)-receptor interaction and the focal adhesion signal pathway. The core genes, including COL1A1, COL1A2, FN1, ITGA1, and ITGB1, were screened for gene correlation and prognosis analysis. The expression level of the five genes was verified by qPCR in glioma. We concluded that these five core genes and ANXA1 could play a synergistic role in gliomas. Conclusion: The results indicated that a high expression level of ANXA1 leads to worse prognosis and ANXA1 is an independent prognostic factor and a potentially important target for the treatment of gliomas.

A methylation-driven genes prognostic signature and the immune microenvironment in epithelial ovarian cancer

Aberrant gene methylation has been implicated in the development and progression of tumors. In this study, we aimed to identity methylation driven genes involved in epithelial ovarian cancer (EOC) to establish a prognostic signature for patients with EOC. We identified and verified 6 MDGs that are closely related to the prognosis of ovarian cancer. A prognostic risk score model and nomogram for predicting the prognosis of ovarian cancer were constructed based on the six MDGs. It can also effectively reflect the immune environment and immunotherapy response of ovarian cancer. These MDGs have great significance to the implementation of individualized treatment and disease monitoring of ovarian cancer patients.

Focal Adhesion-Related Signatures Predict the Treatment Efficacy of Chemotherapy and Prognosis in Patients with Gastric Cancer

Background

The current tumor-node-metastasis (TNM) staging system is insufficient for predicting the efficacy of chemotherapy in patients with gastric cancer (GC). This study aimed to analyze the association between the focal adhesion pathway and therapeutic efficacy of chemotherapy in patients with GC.

Methods

RNA sequencing was performed on 33 clinical samples from patients who responded or did not respond to treatment prior to neoadjuvant chemotherapy. The validation sets containing 696 GC patients with RNA data from three cohorts (PKUCH, TCGA, and GSE14210) were analyzed. A series of machine learning and bioinformatics approaches was combined to build a focal adhesion-related signature model to predict the treatment efficacy and prognosis of patients with GC.

Results

Among the various signaling pathways associated with cancer, focal adhesion was identified as a risk factor related to the treatment efficacy of chemotherapy and prognosis in patients with GC. The focal adhesion-related gene model (FAscore) discriminated patients with a high FAscore who are insensitive to neoadjuvant chemotherapy in our training cohort, and the predicted value was further verified in the GSE14210 cohort. Survival analysis also demonstrated that patients with high FAscores had a relatively shorter survival compared to those with low FAscores. In addition, we found that the levels of tumor mutation burden (TMB) and microsatellite instability (MSI) increased with an increase in FAscore, and the tumor microenvironment (TME) also shifted to a pro-tumor immune microenvironment.

Conclusion

The FAscore model can be used to predict the treatment efficacy of chemotherapy and select appropriate treatment strategies for patients with GC.

Targeted Inhibition of O-Linked β-N-Acetylglucosamine Transferase as a Promising Therapeutic Strategy to Restore Chemosensitivity and Attenuate Aggressive Tumor Traits in Chemoresistant Urothelial Carcinoma of the Bladder

Acquisition of acquired chemoresistance during treatment cycles in urothelial carcinoma of the bladder (UCB) is the major cause of death through enhancing the risk of cancer progression and metastasis. Elevated glucose flux through the abnormal upregulation of O-linked β-N-acetylglucosamine (O-GlcNAc) transferase (OGT) controls key signaling and metabolic pathways regulating diverse cancer cell phenotypes. This study showed that OGT expression levels in two human UCB cell models with acquired resistance to gemcitabine and paclitaxel were significantly upregulated compared with those in parental cells. Reducing hyper-O-GlcNAcylation by OGT knockdown (KD) markedly facilitated chemosensitivity to the corresponding chemotherapeutics in both cells, and combination treatment with OGT-KD showed more severe growth defects in chemoresistant sublines. We subsequently verified the suppressive effects of OGT-KD monotherapy on cell migration/invasion in vitro and xenograft tumor growth in vivo in chemoresistant UCB cells. Transcriptome analysis of these cells revealed 97 upregulated genes, which were enriched in multiple oncogenic pathways. Our final choice of suspected OGT glycosylation substrate was VCAN, S1PR3, PDGFRB, and PRKCG, the knockdown of which induced cell growth defects. These findings demonstrate the vital role of dysregulated OGT activity and hyper-O-GlcNAcylation in modulating treatment failure and tumor aggression in chemoresistant UCB.

TGIF1 overexpression promotes glioma progression and worsens patient prognosis

Transforming growth factor β-induced factor homeobox 1 (TGIF1) reportedly promotes the pathological processes of various malignant tumors. However, few studies have investigated the role of TGIF1 in gliomas. We aimed to explore the relationship between TGIF1 expression and the clinical characteristics of patients with glioma, including their overall survival. A total of thousands transcriptome datapoints were downloaded from public databases to determine the correlations between TGIF1 and various clinicopathological features using the Wilcoxon or Kruskal-Wallis tests. The Kaplan-Meier and Cox statistical methods were used to explore the prognostic significance of TGIF1. Gene set enrichment analysis (GSEA) was used to indirectly identify the pathological mechanisms modulated by TGIF1, and compounds that inhibit its expression were determined using a connectivity map (CMap). TGIF1 was significantly overexpressed in gliomas and was correlated with unfavorable prognostic factors and shorter overall survival. Cox analysis confirmed that TGIF1 expression was a significant predictor of poor prognosis in patients with glioma. GSEA revealed that the signaling pathways associated with TGIF1 expression in glioma included extracellular matrix receptor- and cell cycle-modulating proteins. CMap analysis showed that the small molecules scriptaid, torasemide, dexpropranolol, ipratropium bromide, and harmine were potential negative regulators of TGIF1. Finally, in vitro experiments demonstrated that knockdown of TGIF1 significantly inhibited the proliferation and invasion of glioma cell. Taken together, our study, which is the first to comprehensively analyze TGIF1 in gliomas, revealed it to be a novel oncogene in terms of its association with this disease. As such, TGIF1 may be a potential therapeutic target for individualized treatment of patients with glioma.

A Predictive Model for Prognosis and Therapeutic Response in Hepatocellular Carcinoma Based on a Panel of Three MED8-Related Immunomodulators

The current tumor-node-metastasis (TNM) system is limited in predicting the survival and guiding the treatment of hepatocellular carcinoma (HCC) patients since the TNM system only focuses on the anatomical factors, regardless of the intratumoral molecule heterogeneity. Besides, the landscape of intratumoral immune genes has emerged as a prognostic indicator. The mediator complex subunit 8 (MED8) is a major polymerase regulator and has been described as an oncogene in renal cell carcinoma, but its pathophysiological significance of HCC and its contribution to the prognosis of HCC remain unclear. Here, we aimed to discuss the expression profile and clinical correlation of MED8 in HCC and construct a predictive model based on MED8-related immunomodulators as a supplement to the TNM system. According to our analyses, MED8 was overexpressed in HCC tissues and increased expression of MED8 was an indicator of poor outcome in HCC. The knockdown of MED8 weakened the proliferation, colony forming, and migration of HepG2 and Huh7 cells. Subsequently, a predictive model was identified based on a panel of three MED8-related immunomodulators using The Cancer Genome Atlas (TCGA) database and further validated in International Cancer Genome Consortium (ICGC) database. The combination of the predictive model and the TNM system could improve the performance in predicting the survival of HCC patients. High-risk patients had poor overall survival in TCGA and ICGC databases, as well as in subgroup analysis with early clinicopathology classification. It was also found that high-risk patients had a higher probability of recurrence in TCGA cohort. Furthermore, low-risk score indicated a better response to immunotherapy and drug therapy. This predictive model can be served as a supplement to the TNM system and may have implications in prognosis stratification and therapeutic guidance for HCC.

The regulatory landscape of retinoblastoma: a pathway analysis perspective

Retinoblastoma (Rb) is a rare intraocular tumour in early childhood, with an approximate incidence of 1 in 18 000 live births. Experimental studies for Rb are complex due to the challenges associated with obtaining a normal retina to contrast with diseased tissue. In this work, we reanalyse a dataset that contains normal retina samples. We identified the individual genes whose expression is different in Rb in contrast with normal tissue, determined the pathways whose global expression pattern is more distant from the global expression observed in normal tissue, and finally, we identified which transcription factors regulate the highest number of differentially expressed genes (DEGs) and proposed as transcriptional master regulators (TMRs). The enrichment of DEGs in the phototransduction and retrograde endocannabinoid signalling pathways could be associated with abnormal behaviour of the processes leading to cellular differentiation and cellular proliferation. On the other hand, the TMRs nuclear receptor subfamily 5 group A member 2 and hepatocyte nuclear factor 4 gamma are involved in hepatocyte differentiation. Therefore, the enrichment of aberrant expression in these transcription factors could suggest an abnormal retina development that could be involved in Rb origin and progression.

Marine Migrastatics: A Comprehensive 2022 Update

Metastasis is responsible for the bad prognosis in cancer patients. Advances in research on metastasis prevention focus attention on the molecular mechanisms underlying cancer cell motility and invasion to improve therapies for long-term survival in cancer patients. The so-called “migrastatics” could help block cancer cell invasion and lead to the rapid development of antimetastatic therapies, improving conventional cancer therapies. In the relentless search for migrastatics, the marine environment represents an important source of natural compounds due to its enormous biodiversity. Thus, this review is a selection of scientific research that has pointed out in a broad spectrum of in vitro and in vivo models the anti-cancer power of marine-derived products against cancer cell migration and invasion over the past five years. Overall, this review might provide a useful up-to-date guide about marine-derived compounds with potential interest for pharmaceutical and scientific research on antimetastatic drug endpoints.

Recurrence of Papillary Thyroid Cancer: A Systematic Appraisal of Risk Factors

CONTEXT

Thyroid cancer recurrence is associated with increased mortality and adverse outcomes. Recurrence risk is currently predicted using clinical tools, often restaging patients after treatment. Detailed understanding of recurrence risk at disease onset could lead to personalized and improved patient care.

OBJECTIVE

We aimed to perform a comprehensive bioinformatic and experimental analysis of 3 levels of genetic change (mRNA, microRNA, and somatic mutation) apparent in recurrent tumors and construct a new combinatorial prognostic risk model.

METHODS

We analyzed The Cancer Genome Atlas data (TCGA) to identify differentially expressed genes (mRNA/microRNA) in 46 recurrent vs 455 nonrecurrent thyroid tumors. Two exonic mutational pipelines were used to identify somatic mutations. Functional gene analysis was performed in cell-based assays in multiple thyroid cell lines. The prognostic value of genes was evaluated with TCGA datasets.

RESULTS

We identified 128 new potential biomarkers associated with recurrence, including 40 mRNAs, 39 miRNAs, and 59 genetic variants. Among differentially expressed genes, modulation of FN1, ITGα3, and MET had a significant impact on thyroid cancer cell migration. Similarly, ablation of miR-486 and miR-1179 significantly increased migration of TPC-1 and SW1736 cells. We further utilized genes with a validated functional role and identified a 5-gene risk score classifier as an independent predictor of thyroid cancer recurrence.

CONCLUSION

Our newly proposed risk model based on combinatorial mRNA and microRNA expression has potential clinical utility as a prognostic indicator of recurrence. These findings should facilitate earlier prediction of recurrence with implications for improving patient outcome by tailoring treatment to disease risk and increasing posttreatment surveillance.

Folic Acid Treatment Directly Influences the Genetic and Epigenetic Regulation along with the Associated Cellular Maintenance Processes of HT-29 and SW480 Colorectal Cancer Cell Lines

Simple Summary

Folic acid (FA) participates in DNA synthesis and in DNA methylation; hence, it has a dual role in established neoplasms. We aimed to observe this phenomenon on FA-treated colorectal cancer cell lines (HT-29, SW480). Our results demonstrated that the maintenance processes, namely cell proliferation, cell viability, and DNA repair, were altered in HT-29 cells for short-term FA supplementation, while genetic and epigenetic regulations of SW480 cells were also affected. Despite the fact that FA is a precursor molecule in methyl donor formation, DNA methylation alterations were observed in both directions, primarily influencing the pathways of carcinogenesis. Moreover, behind the great number of differentially expressed genes, other FA-related effects than promoter methylation were suspected. All of our results point beyond the attributes related to FA so far. The different response of the two cell lines is worth considering in clinical practice to facilitate the effectiveness of therapy in the case of tumor heterogeneity.

Abstract

Folic acid (FA) is a synthetic form of vitamin B9, generally used as a nutritional supplement and an adjunctive medication in cancer therapy. FA is involved in genetic and epigenetic regulation; therefore, it has a dual modulatory role in established neoplasms. We aimed to investigate the effect of short-term (72 h) FA supplementation on colorectal cancer; hence, HT-29 and SW480 cells were exposed to different FA concentrations (0, 100, 10,000 ng/mL). HT-29 cell proliferation and viability levels elevated after 100 ng/mL but decreased for 10,000 ng/mL FA. Additionally, a significant (p ≤ 0.05) improvement of genomic stability was detected in HT-29 cells with micronucleus scoring and comet assay. Conversely, the FA treatment did not alter these parameters in SW480 samples. RRBS results highlighted that DNA methylation changes were bidirectional in both cells, mainly affecting carcinogenesis-related pathways. Based on the microarray analysis, promoter methylation status was in accordance with FA-induced expression alterations of 27 genes. Our study demonstrates that the FA effect was highly dependent on the cell type, which can be attributed to the distinct molecular background and the different expression of proliferation- and DNA-repair-associated genes (YWHAZ, HES1, STAT3, CCL2). Moreover, new aspects of FA-regulated DNA methylation and consecutive gene expression were revealed.

Plasma bioscience for medicine, agriculture and hygiene applications

Nonthermal biocompatible plasma (NBP) sources operating in atmospheric pressure environments and their characteristics can be used for plasma bioscience, medicine, and hygiene applications, especially for COVID-19 and citizen. This review surveyed the various NBP sources, including a plasma jet, micro-DBD (dielectric barrier discharge) and nanosecond discharged plasma. The electron temperatures and the plasma densities, which are produced using dielectric barrier discharged electrode systems, can be characterized as 0.7 ~ 1.8 eV and (3-5) × 1014-15 cm-3, respectively. Herein, we introduce a general schematic view of the plasma ultraviolet photolysis of water molecules for reactive oxygen and nitrogen species (RONS) generation inside biological cells or living tissues, which would be synergistically important with RONS diffusive propagation into cells or tissues. Of the RONS, the hydroxyl radical [OH] and hydrogen peroxide H2O2 species would mainly result in apoptotic cell death with other RONS in plasma bioscience and medicines. The diseased biological protein, cancer, and mutated cells could be treated by using a NBP or plasma activated water (PAW) resulting in their apoptosis for a new paradigm of plasma medicine.

A Transcriptome Sequencing Study on Genome-Wide Gene Expression Differences of Lung Cancer Cells Modulated by Fucoidan

Fucoidan has received increasing attention in anti-(lung) tumors. However, the effect of fucoidan on the gene changes of lung cancer cells (LCCs) has not been examined systematically. Herein, we investigate the effect of fucoidan on the phenotypes of LCCs and their gene expression by transcriptome sequencing analysis. The phenotypes of LCCs are significantly inhibited by fucoidan. Importantly, compared to LCCs, 1 mg/ml fucoidan has no effect on the phenotypes of normal cells. Further, 6,930 differentially expressed genes (DEGs) in the transcriptome of LCCs (3,501 up-regulated and 3,429 down-regulated genes) are detected via RNA-sequencing between the fucoidan and control groups. Gene Ontology analysis confirms that DEGs are reflected in DNA replication, cell-substrate junction, regulation of cell cycle phase transition, apoptosis, focal adhesion, cadherin binding, and cell adhesion molecule binding. Thus, our findings on the transcriptomic level highlight the therapeutic potential of fucoidan for lung cancer treatment.

BGN May be a Potential Prognostic Biomarker and Associated With Immune Cell Enrichment of Gastric Cancer

Background: Biglycan (BGN) plays a role in the occurrence and progression of several malignant tumors, though its role in gastric cancer (GC) remains unclear. The objective of this study was to investigate BGN expression, its role in GC prognosis, and immune infiltration. Material and Methods: Gene expression data and corresponding clinical information were downloaded from TCGA and GTEx, respectively. We compared the expression of BGN in GC and normal tissues and verified the differential expression via Real-Time PCR and immunohistochemistry. BGN-related differentially expressed genes (DEGs) were identified. Additionally, the relationships between BGN gene expression and clinicopathological variables and survival in patients with GC were also investigated through univariate and multivariate Cox regression analyses. Finally, we established a predictive model that could well predict the probability of 1-, 3-, and 5-years survival in GC. Results: We found a significantly higher expression of BGN in GC than that in normal tissues (p < 0.001), which was verified by Real-Time PCR (p < 0.01) and immunohistochemistry (p < 0.001). The 492 identified DEGs were primarily enriched in pathways related to tumor genesis and metastasis, including extracellular matrix (ECM)-receptor interaction, focal adhesion pathway, Wnt signaling, and signaling by VEGF. BGN expression was positively correlated with the enrichment of the NK cells (r = 0.620, p < 0.001) and macrophages (r = 0.550, p < 0.001), but negatively correlated with the enrichment of Th17 cells (r = 0.250, p < 0.001). BGN expression was also significantly correlated with histologic grade (GI&G2 vs. G3, p < 0.001), histologic type (Diffuse type vs. Tubular type, p < 0.001), histologic stage (stage I vs. stage II and stage I vs. stage III, p < 0.001), T stage (T1 vs. T2, T1 vs. T3, and T1 vs. T4, p < 0.001) and Helicobacter pylori (HP) infection (yes vs. no, p < 0.05) in GC. High BGN expression showed significant association with poor overall survival (OS) in GC patients (HR = 1.53 (1.09-2.14), p = 0.013). The constructed nomogram can well predict the 1-, 3-, and 5-years overall survival probability of GC patients (C-index = 0.728). Conclusion: BGN plays an important role in the occurrence and progression of GC and is a potential biomarker for the diagnosis and treatment of GC.

Plasticity within aldehyde dehydrogenase-positive cells determines prostate cancer radiosensitivity

Tumor heterogeneity and cellular plasticity are key determinants of tumor progression, metastatic spread and therapy response driven by the cancer stem cell (CSC) population. Within the present study, we analyzed irradiation-induced plasticity within the aldehyde dehydrogenase (ALDH)-positive population in prostate cancer (PCa). The radiosensitivity of xenograft tumors derived from ALDH+ and ALDH- cells was determined with local tumor control analyses and demonstrated different dose-response profiles, time to relapse and focal adhesion signaling. The transcriptional heterogeneity was analyzed in pools of ten DU145 and PC3 cells with multiplex gene expression analyses and illustrated a higher degree of heterogeneity within the ALDH+ population that even increases upon irradiation in comparison to ALDH- cells. Phenotypic conversion and clonal competition were analyzed with fluorescence protein-labeled cells to distinguish cellular origins in competitive 3D cultures and xenograft tumors. We found that the ALDH+ population outcompetes ALDH- cells and drives tumor growth, in particular upon irradiation. The observed dynamics of the cellular state compositions between ALDH+ and ALDH- cells in vivo before and after tumor irradiation was reproduced by a probabilistic Markov compartment model that incorporates cellular plasticity, clonal competition and phenotype-specific radiosensitivities. Transcriptional analyses indicate that the cellular conversion from ALDH- into ALDH+ cells within xenograft tumors under therapeutic pressure was partially mediated through induction of the transcriptional repressor SNAI2. In summary, irradiation-induced cellular conversion events are present in xenograft tumors derived from PCa cells and may be responsible for radiotherapy failure.

Growth factor receptor and β1 integrin signaling differentially regulate basal clonogenicity and radiation survival of fibroblasts via a modulation of cell cycling

Cell adhesion to extracellular matrix proteins mediates resistance to radio- and chemotherapy by activating integrin signaling. In addition, mutual and cooperative interactions between integrin and growth factor receptor signaling contribute to the cellular radiation response. Here, we investigate to which extend the crosstalk between β1 integrins and growth factor receptor signaling determines the cellular radiation response of fibroblasts by assessing clonogenic survival and cell cycling. By utilizing growth factor signaling competent and either β1 integrin wildtype GD25β1A fibroblasts or β1 integrin mutant, signaling incompetent GD25β1B fibroblasts, we show basal clonogenic survival to depend on growth factor receptor but not integrin signaling. Our data further suggest the cooperation between β1 integrins and growth factor receptors to be critical for enhancing the radiation-induced G2/M cell cycle block leading to improved clonogenic radiation survival. By pharmacological inhibition of EGFR and PI3K, we additionally show that the essential contribution of EGFR signaling to radiogenic G2/M cell cycle arrest depends on the co-activation of the β1 integrin signaling axis, but occurs independent of PI3K. Taken together, elucidation of the signaling circuitry underlying the EGFR/β1 integrin crosstalk may support the development of advanced molecular targeted therapies for radiation oncology.

Screening and Comprehensive Analysis of Cancer-Associated tRNA-Derived Fragments

tRNA-derived fragments (tRFs) constitute a novel class of small non-coding RNA cleaved from tRNAs. In recent years, researches have shown the regulatory roles of a few tRFs in cancers, illuminating a new direction for tRF-centric cancer researches. Nonetheless, more specific screening of tRFs related to oncogenesis pathways, cancer progression stages and cancer prognosis is continuously demanded to reveal the landscape of the cancer-associated tRFs. In this work, by combining the clinical information recorded in The Cancer Genome Atlas (TCGA) and the tRF expression profiles curated by MINTbase v2.0, we systematically screened 1,516 cancer-associated tRFs (ca-tRFs) across seven cancer types. The ca-tRF set collectively combined the differentially expressed tRFs between cancer samples and control samples, the tRFs significantly correlated with tumor stage and the tRFs significantly correlated with patient survival. By incorporating our previous tRF-target dataset, we found the ca-tRFs tend to target cancer-associated genes and onco-pathways like ATF6-mediated unfolded protein response, angiogenesis, cell cycle process regulation, focal adhesion, PI3K-Akt signaling pathway, cellular senescence and FoxO signaling pathway across multiple cancer types. And cell composition analysis implies that the expressions of ca-tRFs are more likely to be correlated with T-cell infiltration. We also found the ca-tRF expression pattern is informative to prognosis, suggesting plausible tRF-based cancer subtypes. Together, our systematic analysis demonstrates the potentially extensive involvements of tRFs in cancers, and provides a reasonable list of cancer-associated tRFs for further investigations.

Exosomes of A549 Cells Induced Migration, Invasion, and EMT of BEAS-2B Cells Related to let-7c-5p and miR-181b-5p

As carriers containing abundant biological information, exosomes could deliver the property of donor cells to recipient cells. Emerging studies have shown that tumor cells could secrete a mass of exosomes into the microenvironment to regulate bystander cells. However, the underlying mechanisms of such a phenomenon remain largely unexplored. In this research, we purified and identified the exosomes of A549 cells and found that A549-cell-derived exosomes promoted BEAS-2B cells migration, invasion, and epithelial-mesenchymal transition (EMT). Importantly, we observed that let-7c-5p and miR-181b-5p were attenuated in A549-cell-derived exosomes compared to BEAS-2B-cell-derived exosomes. The analysis of miRNA expression level in BEAS-2B cells indicated that incubation with A549-cell-derived exosomes reduced the expression levels of let-7c-5p and miR-181b-5p. In transient transfections assay, we found that downregulation of let-7c-5p and miR-181b-5p simultaneously showed stronger promotion of BEAS-2B cells migration and invasion than individually. Moreover, exosomes secreted from A549 cells with upregulated expression of let-7c-5p and miR-181b-5p significantly reduce their regulatory effect on BEAS-2B cells. Bioinformatics analyses revealed that let-7c-5p and miR-181b-5p inhibit the EMT process mainly by regulating focal adhesion and mitogen-activated protein kinase (MAPK) signaling pathway. Thus, our data demonstrated that A549-cell-derived exosomal let-7c-5p and miR-181b-5p could induce migration, invasion, and EMT in BEAS-2B cells, which might be regulated through focal adhesion and MAPK signaling pathway. The expression level of let-7c-5p and miR-181b-5p may show great significance for the early diagnosis of lung cancer.

Proteomic Screening and Verification of Biomarkers in Different Stages of Mycosis Fungoides: A pilot Study

Mycosis fungoides (MF) is the most common cutaneous T-cell lymphoma; in advanced stages, it can involve multiple organs and has a poor prognosis. Early detection of the disease is still urgent, but there is no optimal therapy for advanced MF. In the present study, quantitative proteomic analyses (label-free quantitation, LFQ) were applied to tissue samples of different stages of MF and tissue samples from controls (eczema patients and healthy donors) to conduct preliminary molecular analysis to clarify the pathogenesis of the disease. Differential protein expression analysis demonstrated that 113 and 305 proteins were associated with the early and advanced stages of MF, respectively. Gene ontology (GO) enrichment analysis was conducted to determine the potential functions of the proteins, which could be classified into three categories: biological process, cellular component, and molecular function. The results revealed that a series of biological processes, including “initiation of DNA replication” and “nucleosome assembly,” were involved in the disease. Moreover, cellular components, including the “desmosome” and “integrin complex,” may affect the invasion and metastasis of MF via molecular functions, including “integrin binding” and “cadherin binding”. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis demonstrated that “focal adhesion DNA replication,” “Toll-like receptor signalling pathway” and other pathways were also involved. A parallel reaction monitoring (PRM) assay was applied to validate the identified differentially expressed proteins. In conclusion, the above proteomic findings may have great diagnostic and prognostic value in diverse malignancies, especially MF. Nevertheless, further studies are still needed to explore the precise mechanisms of MF.

Precision N-glycoproteomics reveals elevated LacdiNAc as a novel signature of intrahepatic cholangiocarcinoma

Primary liver cancer, mainly comprising hepatocellular carcinoma (HCC) and intrahepatic cholangiocarcinoma (ICC), remains a major global health problem. Although ICC is clinically different from HCC, their molecular differences are still largely unclear. In this study, precision N‐glycoproteomic analysis was performed on both ICC and HCC tumors as well as paracancer tissues to investigate their aberrant site‐specific N‐glycosylation. By using our newly developed glycoproteomic methods and novel algorithm, termed ‘StrucGP’, a total of 486 N‐glycan structures attached on 1235 glycosites were identified from 894 glycoproteins in ICC and HCC tumors. Notably, glycans with uncommon LacdiNAc (GalNAcβ1‐4GlcNAc) structures were distinguished from their isomeric glycans. In addition to several bi‐antennary and/or bisecting glycans that were commonly elevated in ICC and HCC, a number of LacdiNAc‐containing, tri‐antennary, and core‐fucosylated glycans were uniquely increased in ICC. More interestingly, almost all LacdiNAc‐containing N‐glycopeptides were enhanced in ICC tumor but not in HCC tumor, and this phenomenon was further confirmed by lectin histochemistry and the high expression of β1‐4 GalNAc transferases in ICC at both mRNA and protein expression levels. The novel N‐glycan alterations uniquely detected in ICC provide a valuable resource for future studies regarding to the discovery of ICC diagnostic biomarkers, therapeutic targets, and mechanism investigations.

Prognostic Value of SLC16A3(MCT4) in Lung Adenocarcinoma and Its Clinical Significance

Background

Currently, lung adenocarcinoma is the most common form of lung cancer. Although the pathogenesis of lung adenocarcinoma is progressing rapidly, the mortality rate of lung adenocarcinoma is still high. Therefore, it is necessary to search for a new biomarker to guide the prognosis of lung adenocarcinoma.

Methods

The significance of SLC16A3 in lung adenocarcinoma was investigated by multi-database analysis. GEPIA, UALCAN, TIMER, Cbioportal, and R software were used for research.

Results

Our study found that SLC16A3 was highly expressed in lung adenocarcinoma and was associated with poor prognosis. Further studies have shown that SLC16A3 is involved in some metabolic pathways. Not only that, SLC16A3 is associated with immune cell infiltration and tumor mutation burden (TMB).

Conclusion

SLC16A3 has good prognostic significance in lung adenocarcinoma, based on which to explore treatment options may improve the prognosis of patients.

Identification of potential core genes in gastric cancer using bioinformatics analysis

Background

Gastric cancer is the third leading cause of cancer-related mortality in China. Most patients with gastric cancer have no obvious early symptoms; thus, many of them are in the middle and late stages of gastric cancer at first diagnosis and miss the best treatment opportunity. Molecular targeted therapy is particularly important in changing this status quo.

Methods

Three microarray datasets (GSE29272, GSE33651, and GSE54129) were selected from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) were screened using GEO2R. The Database for Annotation, Visualization and Integrated Discovery (DAVID) was used to analyze the functional features of these DEGs and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment. The protein-protein interaction (PPI) of these DEGs was visualized by Cytoscape software. The expressions of hub genes were evaluated based on Gene Expression Profiling Interactive Analysis (GEPIA). Moreover, we used the online Kaplan-Meier plotter survival analysis tool to evaluate the prognostic values of hub genes. The Target Scan database was used to predict microRNAs that could regulate the target gene, collagen type IV alpha 1 chain (COL4A1). The OncomiR database was used to analyze the expression levels of three microRNAs, as well as the relationships with tumor stage, grade, and prognosis.

Results

We identified 78 DEGs, including 53 upregulated genes and 25 downregulated genes. The DEGs were mainly enriched in extracellular matrix organization, extracellular structure organization, and response to wounding. Moreover, three KEGG pathways were markedly enriched, including focal adhesion, complement and coagulation cascades, and extracellular matrix (ECM)-receptor interaction. Among these 78 genes, we selected 10 hub genes. The overexpression levels of these hub genes were closely related to poor prognosis and the development of gastric cancer (except for COL3A1, LOX, and CXCL8). Moreover, we found that microRNA-29a-3p, miR-29b-3p, and miR-29c-3p were the potential microRNAs that could regulate the target gene, COL4A1.

Conclusions

Our results showed that FN1, COL1A1, TIMP1, COL1A2, SPARC, COL4A1, and SERPINE1 could contribute to the development of novel molecular targets and biomarker-driven treatments for gastric cancer.

Acetylenic Synthetic Betulin Derivatives Inhibit Akt and Erk Kinases Activity, Trigger Apoptosis and Suppress Proliferation of Neuroblastoma and Rhabdomyosarcoma Cell Lines

Neuroblastoma (NB) and rhabdomyosarcoma (RMS), the most common pediatric extracranial solid tumors, still represent an important clinical challenge since no effective treatment is available for metastatic and recurrent disease. Hence, there is an urgent need for the development of new chemotherapeutics to improve the outcome of patients. Betulin (Bet), a triterpenoid from the bark of birches, demonstrated interesting anti-cancer potential. The modification of natural phytochemicals with evidenced anti-tumor activity, including Bet, is one of the methods of receiving new compounds for potential implementation in oncological treatment. Here, we showed that two acetylenic synthetic Bet derivatives (ASBDs), EB5 and EB25/1, reduced the viability and proliferation of SK-N-AS and TE671 cells, as measured by MTT and BrdU tests, respectively. Moreover, ASBDs were also more cytotoxic than temozolomide (TMZ) and cisplatin (cis-diaminedichloroplatinum [II], CDDP) in vitro, and the combination of EB5 with CDDP enhanced anti-cancer effects. We also showed the slowdown of cell cycle progression at S/G2 phases mediated by EB5 using FACS flow cytometry. The decreased viability and proliferation of pediatric cancers cells after treatment with ASBDs was linked to the reduced activity of kinases Akt, Erk1/2 and p38 and the induction of apoptosis, as investigated using Western blotting and FACS. In addition, in silico analyses of the ADMET profile found EB5 to be a promising anti-cancer drug candidate that would benefit from further investigation.

Cross-Talk of Focal Adhesion-Related Gene Defines Prognosis and the Immune Microenvironment in Gastric Cancer

Background: Focal adhesion, as the intermediary between tumor cells and extracellular matrix communication, plays a variety of roles in tumor invasion, migration, and drug resistance. However, the potential role of focal adhesion-related genes in the microenvironment, immune cell infiltration, and drug sensitivity of gastric cancer (GC) has not yet been revealed.

Methods: The genetic and transcriptional perspectives of focal adhesion-related genes were systematically analyzed. From a genetic perspective, the focal adhesion index (FAI) was constructed based on 18 prognosis-related focus adhesion-related genes to evaluate the immune microenvironment and drug sensitivity. Then three prognosis-related genes were used for consistent clustering to identify GC subtypes. Finally, use FLT1, EGF, COL5A2, and M2 macrophages to develop risk signatures, and establish a nomogram together with clinicopathological characteristics.

Results: Mutations in the focal adhesion-related gene affect the survival time and clinical characteristics of GC patients. FAI has been associated with a shorter survival time, immune signaling pathways, M2 macrophage infiltration, epithelial-mesenchymal transition (EMT) signaling, and diffuse type of GC. FAI recognizes ALK, cell cycle, and BMX signaling pathways inhibitors as sensitive agents for the treatment of GC. FLT1, EGF, and COL5A2 may distinguish GC subtypes. The established risk signature is of great significance to the prognostic evaluation of GC based on FLT1, EGF, and COL5A2 and M2 macrophage expression.

Conclusion: The focal adhesion-related gene is a potential biomarker for the evaluation of the immune microenvironment and prognosis. This work emphasizes the potential impact of the focal adhesion pathway in GC therapy and highlights its guiding role in prognostic evaluation.

An Update on Molecular Pathways Regulating Vasculogenic Mimicry in Human Osteosarcoma and Their Role in Canine Oncology

Canine tumors are valuable comparative models for human counterparts, especially to explore novel biomarkers and to understand pathways and processes involved in metastasis. Vasculogenic mimicry (VM) is a unique property of malignant cancer cells which promote metastasis. Thus, it represents an opportunity to investigate both the molecular mechanisms and the therapeutic targets of a crucial phenotypic malignant switch. Although this biological process has been largely investigated in different human cancer types, including osteosarcoma, it is still largely unknown in veterinary pathology, where it has been mainly explored in canine mammary tumors. The presence of VM in human osteosarcoma is associated with poor clinical outcome, reduced patient survival, and increased risk of metastasis and it shares the main pathways involved in other type of human tumors. This review illustrates the main findings concerning the VM process in human osteosarcoma, search for the related current knowledge in canine pathology and oncology, and potential involvement of multiple pathways in VM formation, in order to provide a basis for future investigations on VM in canine tumors.

A Focal Adhesion-Related Gene Signature Predicts Prognosis in Glioma and Correlates With Radiation Response and Immune Microenvironment

Background

Glioma is the most frequent brain malignancy presenting very poor prognosis and high recurrence rate. Focal adhesion complexes play pivotal roles in cell migration and act as hubs of several signaling pathways.

Methods

We used bioinformatic databases (CGGA, TCGA, and GEO) and identified a focal adhesion-related differential gene expression (FADG) signature by uniCox and LASSO regression analysis. We calculated the risk score of every patient using the regression coefficient value and expression of each gene. Survival analysis, receiver operating characteristic curve (ROC), principal component analysis (PCA), and stratified analysis were used to validate the FADG signature. Then, we conducted GSEA to identify the signaling pathways related to the FADG signature. Correlation analysis of risk scores between the immune checkpoint was performed. In addition, the correlation of risk scores and genes related with DNA repair was performed. CIBERSORT and ssGSEA were used to explore the tumor microenvironment (TME). A nomogram that involved our FADG signature was also constructed.

Results

In total, 1,726 (528 patients diagnosed with WHO II, 591 WHO III, and 603 WHO IV) cases and 23 normal samples were included in our study. We identified 29 prognosis-related genes in the LASSO analysis and constructed an eight FADG signature. The results from the survival analysis, stratified analysis, ROC curve, and univariate and multivariate regression analysis revealed that the prognosis of the high-risk group was significantly worse than the low-risk group. Correlation analysis between risk score and genes that related with DNA repair showed that the risk score was positively related with BRCA1, BRCA2, RAD51, TGFB1, and TP53. Besides, we found that the signature could predict the prognosis of patients who received radiation therapy. SsGSEA indicated that the high-risk score was positively correlated with the ESTIMATE, immune, and stromal scores but negatively correlated with tumor purity. Notably, patients in the high-risk group had a high infiltration of immunocytes. The correlation analysis revealed that the risk score was positively correlated with B7-H3, CTLA4, LAG3, PD-L1, and TIM3 but inversely correlated with PD-1.

Conclusion

The FADG signature we constructed could provide a sensitive prognostic model for patients with glioma and contribute to improve immunotherapy management guidelines.

A comprehensive analysis of different gene classes in pancreatic cancer: SIGLEC15 may be a promising immunotherapeutic target

BACKGROUND

Pancreatic cancer (PC) is one of the most lethal cancer types with an extremely poor diagnosis and prognosis. This study aimed to comprehensively analyze the relationships between PC and different gene classes.

METHODS

Numerous genes from different categories were selected from the UALCAN database. Expression and survival analysis of these genes were performed via GEPIA, starBase and Kaplan-Meier Plotter tools. The correlations between PC-related genes and frequently mutated genes in PC as well as myeloid-derived suppressor cells (MDSCs) infiltration levels were explored by TIMER tool. The associations between PC-related genes, immune checkpoints and 182 core cancer-intrinsic CTLs-evasion genes were analyzed by R software. Besides, KEGG analysis were performed for the PC-related genes.

RESULTS

14 genes were identified to be highly expressed in pancreatic cancer and significantly associated with poor prognosis. Besides, high expression of these genes were observed in patients with KRAS or TP53 mutations. Most genes were significantly positively associated with immune checkpoint SIGLEC15, however, showed negative relations to PDCD1, CTLA4, LAG3, TIGIT, PDCD1LG2. In addition, all 14 genes exhibited close relationships with MDSC infiltration levels and various core cancer-intrinsic CTLs-evasion genes, especially DNTTIP1, FADD, ARF6, BCL2L1, CEP55, GALE, PDCD6IP, and RCE1. We also explored the most related pathways with these genes to further reveal the pathogenesis and metastatic mechanisms of PC.

CONCLUSION

Our study analyzed the relationships between 14 PC-related genes and pancreatic cancer from different angles, which may contribute to a better understanding of unsolved mystery in PC.

Identification and Validation of a PPP1R12A-Related Five-Gene Signature Associated With Metabolism to Predict the Prognosis of Patients With Prostate Cancer

Background

Prostate cancer (PCa) is the most common malignant male neoplasm in the American male population. Our prior studies have demonstrated that protein phosphatase 1 regulatory subunit 12A (PPP1R12A) could be an efficient prognostic factor in patients with PCa, promoting further investigation. The present study attempted to construct a gene signature based on PPP1R12A and metabolism-related genes to predict the prognosis of PCa patients.

Methods

The mRNA expression profiles of 499 tumor and 52 normal tissues were extracted from The Cancer Genome Atlas (TCGA) database. We selected differentially expressed PPP1R12A-related genes among these mRNAs. Tandem affinity purification-mass spectrometry was used to identify the proteins that directly interact with PPP1R12A. Gene set enrichment analysis (GSEA) was used to extract metabolism-related genes. Univariate Cox regression analysis and a random survival forest algorithm were used to confirm optimal genes to build a prognostic risk model.

Results

We identified a five-gene signature (PPP1R12A, PTGS2, GGCT, AOX1, and NT5E) that was associated with PPP1R12A and metabolism in PCa, which effectively predicted disease-free survival (DFS) and biochemical relapse-free survival (BRFS). Moreover, the signature was validated by two internal datasets from TCGA and one external dataset from the Gene Expression Omnibus (GEO).

Conclusion

The five-gene signature is an effective potential factor to predict the prognosis of PCa, classifying PCa patients into high- and low-risk groups, which might provide potential novel treatment strategies for these patients.

Silencing of ER-resident oxidoreductase PDIA3 inhibits malignant biological behaviors of multidrug-resistant gastric cancer

Glycosylation is a common posttranslational modification of proteins, which plays a role in the malignant transformation, growth, progression, chemoresistance, and immune response of tumors. Disulfide isomerase family A3 (PDIA3) specifically acts on newly synthesized glycoproteins to promote the correct folding of sugar chains. Studies have shown that PDIA3 participates in multidrug-resistant gastric cancer (MDR-GC). In this study, we performed western blot analysis and immunohistochemistry to identify PDIA3 expression. Cell proliferation was assessed by CCK-8 assay. Transwell assays were used to detect the migration and invasion abilities of cells. Immunoprecipitation coupled to mass spectrometry (IP-MS) analysis was employed to identify PDIA3-interacting proteins and the associated pathways in MDR-GC cells. Glycoprotein interactions and translocation were detected by immunofluorescence assay. The results showed that PDIA3 knockdown significantly inhibited the proliferation, invasion, and migration abilities of MDR-GC cells. Kyoto Encyclopedia of Genes and Genomes analysis of the IP-MS results showed that PDIA3 was closely associated with focal adhesion pathways in MDR-GC cells. Additionally, important components of focal adhesion pathways, including fibronectin-1 (FN1) and integrin α5 (ITGA5), were identified as pivotal PDIA3-binding glycoproteins. Knockdown of PDIA3 altered the cellular locations of FN1 and ITGA5, leading to abnormal accumulation. In conclusion, our results suggest that knockdown of PDIA3 inhibited the malignant behaviors of MDR-GC cells and influenced the translocation of FN1 and ITGA5.

MicroRNA profiles in five pairs of early gastric cancer tissues and adjacent non-cancerous tissues

Approximately half of the world's gastric cancer cases and deaths occur in China. In addition, the incidence and mortality rates of gastric cancer in Gansu province in China are much higher than the average nationwide levels. The present study investigated microRNA (miRNA/miR) profiles in early gastric cancer (EGC) without specific symptoms. miRNA expression levels in five pairs of EGC tissues and adjacent non-cancerous mucosa tissues of patients from Gansu province in China were analyzed using a miRNA microarray. A total of 47 differentially expressed miRNAs (DEMs) were identified. Subsequently, mRNA expression profiles of three pairs of cancer tissues and adjacent non-cancerous tissues from 3 Asian patients with stage I or stage II gastric cancer (stage I/II; American Joint Committee on Cancer classification, Eighth Edition) were obtained from The Cancer Genome Atlas database, and differentially expressed genes (DEGs) were identified. The target genes of DEMs were filtered from the DEGs using the miRDB database and a miRNA-gene network was constructed. The functions of DEMs were evaluated using the tool for annotations of human miRNAs database, and via Gene Ontology analysis, Kyoto Encyclopedia of Genes and Genomes analysis and Gene Set Enrichment Analysis of the target genes. Finally, survival analyses of DEMs, which were in the miRNA-gene network, was performed. The results suggested that a number of miRNAs, including hsa-let-7a-5p, hsa-miR-27a-3p, hsa-miR-126-5p and hsa-miR-424-5p, may serve critical roles in EGC. The present study could provide a basis for the identification of EGC screening biomarkers. Furthermore, the present study may provide a basis for the exploration of the cause of the high incidence of gastric cancer in Gansu province from the perspective of miRNAs.

Proteomic characteristics of PM2.5 -induced differentially expressed proteins in human renal tubular epithelial cells

Human renal epithelial (HK-2) cells were treated with PM_2.5 (50 μg/mL) from Shenzhen and Taiyuan, proteomics and bioinformatics were used to screen the differentially expressed proteins (DEPs). A total of 577 DEPs were screened after HK-2 cells exposed to Shenzhen PM_2.5, of which 426 were up-regulated and 151 were down-regulated. A total of 1250 DEPs were screened in HK-2 cells after exposure to Taiyuan PM_2.5, of which 488 were up-regulated and 185 were down-regulated. The top 10 proteins with the highest number of nodes were screened using the interaction network map of DEPs. HK-2 cells exposed to Shenzhen PM_2.5 contained CYR61, CTGF, and THBS1 proteins, while HK-2 cells exposed to Taiyuan PM_2.5 contained ALB, FN1, and CYR61 proteins. Additionally, PM_2.5 components were detected, PM_2.5 samples from Shenzhen and Taiyuan induced obvious changes in DEPs expression, the difference in DEPs between the two cities was probably associated with the different PM_2.5 components.

Expression and prognostic potential of PLEK2 in head and neck squamous cell carcinoma based on bioinformatics analysis

Background

PLEK2 (pleckstrin) could bind to membrane‐bound phosphatidylinositols and further promote cell spread. Recently, several studies have noted the importance of PLEK2 in tumor metastasis. However, the role of PLEK2 in head and neck squamous cell carcinoma (HNSCC) remains to be elucidated.

Methods

The PLEK2 expression in HNSCC was identified using Oncomine, Gene Expression Omnibus (GEO), UALCAN databases, and western blot analysis. Prognosis analysis was performed using Kaplan–Meier plotter, DriverDBv3, UALCAN, UCSC Xena, and GEO databases. Single‐cell functional analysis was further performed using the cancerSEA database. The PLEK2‐related co‑expressed genes were identified, and gene set enrichment analysis was performed using LinkedOmics. Furthermore, the top 10 hub genes were identified using the cytoHubba plug‐in of Cytoscape. Then, gene enrichment analysis, pathway activity, and drug sensitivity analyses of the hub genes were performed using the R package “clusterProfiler” and GSCAlite. Finally, the UCSC Xena browser was utilized to explore the hub gene most likely to play a synergic role with PLEK2 in HNSCC.

Results

Elevated expression of PLEK2 was observed in HNSCC and even in HNSCC subgroups based on diverse clinicopathological features, portending a poor prognosis in HNSCC. PLEK2 was correlated with metastasis and hypoxia in HNSCC, and the PLEK2‐related co‐expressed genes were mainly involved in the focal adhesion pathway. The top 10 hub genes were primarily enriched in focal adhesion, HPV infection, ECM‐receptor interaction, and PI3K‐AKT signaling pathway, and epithelial–mesenchymal transition pathway was activated. Furthermore, the expression levels of the hub genes were associated with sensitivity and resistance to various small molecules and anti‐cancer drugs. Further study suggested that ITGA3 and PLEK2 might be viewed as inextricably linked in facilitating HNSCC metastasis.

Conclusions

In general, PLEK2 might serve as a potential biomarker for the diagnosis of HNSCC and guide the development of targeted therapies for HNSCC.

A Novel Six Metastasis-Related Prognostic Gene Signature for Patients With Osteosarcoma

Osteosarcoma is the most common malignant bone tumor, and although there has been significant progress in its management, metastases often herald incurable disease. Here we defined genes differentially expressed between primary and metastatic osteosarcoma as metastasis-related genes (MRGs) and used them to construct a novel six-MRG prognostic signature for overall survival of patients with osteosarcoma. Validation in internal and external datasets confirmed satisfactory accuracy and generalizability of the prognostic model, and a nomogram based on the signature and clinical variables was constructed to aid clinical decision-making. Of the six MRGs, FHIT is a well-documented tumor suppressor gene that is poorly defined in osteosarcoma. Consistent with tumor suppressor function, FHIT was downregulated in osteosarcoma cells and human osteosarcoma samples. FHIT overexpression inhibited osteosarcoma proliferation, migration, and invasion both in vitro and in vivo. Mechanistically, FHIT overexpression upregulate the epithelial marker E-cadherin while repressing the mesenchymal markers N-cadherin and vimentin. Our six-MRG signature represents a novel and clinically useful prognostic biomarker for patients with osteosarcoma, and FHIT might represent a therapeutic target by reversing epithelial to mesenchymal transition.

Data mining-based study of collagen type III alpha 1 (COL3A1) prognostic value and immune exploration in pan-cancer

The extracellular matrix (ECM) shows an essential effect during the occurrence and procession of human cancers. Type III collagen is a crucial component of ECM. Collagen Type III Alpha 1(COL3A1) is aberrantly expressed in a variety of cancers. Nevertheless, the role of COL3A1 in pan-cancer stays unidentified. In this study, we explored public databases, including Cancer Genome Atlas (TCGA), GTEx, GEPIA, cBioPortal, Oncommine, TIMER and GENEMANIA databases to identify the differential expression of COL3A1 in human cancer tissues and normal samples, followed by its prognostic value for patient survival. In addition, we explore the association between COL3A1 expression and immune infiltration. Further, we used the GeneMANIA database and Gene Set Enrichment Analysis (GSEA) to investigate Protein–Protein Interaction (PPI) and gene functional enrichment. Results show that COL3A1 expressed higher in tumor samples than in normal samples. Upregulation of COL3A1 is associated with a worse prognosis and a more advanced cancer stage. COL3A1 expression shows significant positive correlations with tumor-infiltrating immune cells (TIICs), including neutrophils, macrophages, CD8 + T cells, CD4 + T cells, dendritic cells, and B cells. Markers of TIICs demonstrated distinct patterns of COL3A1-related immune infiltration. COL3A1 expression was associated with ECM receptor interaction, regulation of actin cytoskeleton and focal adhesion pathways via GSEA analysis. In conclusion, COL3A1 may be a molecular biomarker for prognosis and immune infiltration in pan-cancer. It might act as a potential target for a new insight of human cancers management.

CAM-DR: Mechanisms, Roles and Clinical Application in Tumors

Despite the continuous improvement of various therapeutic techniques, the overall prognosis of tumors has been significantly improved, but malignant tumors in the middle and advanced stages still cannot be completely cured. It is now evident that cell adhesion-mediated resistance (CAM-DR) limits the success of cancer therapies and is a great obstacle to overcome in the clinic. The interactions between tumor cells and extracellular matrix (ECM) molecules or adjacent cells may play a significant role in initiating the intracellular signaling pathways that are associated with cell proliferation, survival upon binding to their ligands. Recent studies illustrate that these adhesion-related factors may contribute to the survival of cancer cells after chemotherapeutic therapy, advantageous to resistant cells to proliferate and develop multiple mechanisms of drug resistance. In this review, we focus on the molecular basis of these interactions and the main signal transduction pathways that are involved in the enhancement of the cancer cells’ survival. Furthermore, therapies targeting interactions between cancer cells and their environment to enhance drug response or prevent the emergence of drug resistance will also be discussed.

Toward comprehensive functional analysis of gene lists weighted by gene essentiality scores

MOTIVATION

Gene functional enrichment analysis represents one of the most popular bioinformatics methods for annotating the pathways and function categories of a given gene list. Current algorithms for enrichment computation such as Fisher's exact test and hypergeometric test totally depend on the category count numbers of the gene list and one gene set. In this case, whatever the genes are, they were treated equally. However, actually genes show different scores in their essentiality in a gene list and in a gene set. It is thus hypothesized that the essentiality scores could be important and should be considered in gene functional analysis.

RESULTS

For this purpose, here we proposed WEAT (https://www.cuilab.cn/weat/), a weighted gene set enrichment algorithm and online tool by weighting genes using essentiality scores. We confirmed the usefulness of WEAT using three case studies, the functional analysis of one aging-related gene list, one gene list involved in Lung Squamous Cell Carcinoma (LUSC), and one cardiomyopathy gene list from Drosophila model. Finally, we believe that the WEAT method and tool could provide more possibilities for further exploring the functions of given gene lists.

AVAILABILITY

The datasets generated and analyzed during the current study are available on our website at https://www.cuilab.cn/weat/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genome-wide discovery of hidden genes mediating known drug-disease association using KDDANet

Many of genes mediating Known Drug-Disease Association (KDDA) are escaped from experimental detection. Identifying of these genes (hidden genes) is of great significance for understanding disease pathogenesis and guiding drug repurposing. Here, we presented a novel computational tool, called KDDANet, for systematic and accurate uncovering the hidden genes mediating KDDA from the perspective of genome-wide functional gene interaction network. KDDANet demonstrated the competitive performances in both sensitivity and specificity of identifying genes in mediating KDDA in comparison to the existing state-of-the-art methods. Case studies on Alzheimer's disease (AD) and obesity uncovered the mechanistic relevance of KDDANet predictions. Furthermore, when applied with multiple types of cancer-omics datasets, KDDANet not only recapitulated known genes mediating KDDAs related to cancer, but also revealed novel candidates that offer new biological insights. Importantly, KDDANet can be used to discover the shared genes mediating multiple KDDAs. KDDANet can be accessed at http://www.kddanet.cn and the code can be freely downloaded at https://github.com/huayu1111/KDDANet .

Research trends for papillary thyroid carcinoma from 2010 to 2019: A systematic review and bibliometrics analysis

BACKGROUND

Thyroid carcinoma comprises the fastest rising incidence of carcinomas over the past decade. Papillary thyroid carcinoma (PTC) is the most predominant type of thyroid carcinoma. This study aimed to assess the research trends in the field of PTC.

METHODS

Publications from January 2010 to December 2019 were retrieved from the Web of Science Core Collection database using Thompson Reuters. Searching strategies were determined according to Medical Subject Heading terms. Different kinds of bibliometrics software, such as HistCite and VOSviewer, and online bibliometrics analysis platforms were utilized to evaluate and visualize the results.

RESULTS

A total of 8102 publications across 93 countries were identified, with the annual number of publications showing an increasing trend. The United States, China, and South Korea showed their dominant position in PTC publication outputs, H-index, total citations, and international collaborations. Thyroid was the most productive journal. Akira Miyauchi published the most articles, and the most productive institution was Yonsei University. The hotspots keywords proliferation, invasion and metastasis, diagnoses and prognoses, therapeutic resistance, recurrence, and microcarcinomas appeared earlier and were sustained over the last 3 years.

CONCLUSIONS

This bibliometric study provides a comprehensive analysis delineating the scientific productivity, collaboration, and research hotspots within the PTC field, which will be very helpful when focusing on the direction of research over the next few years.

TEAD4 is an Immune Regulating-Related Prognostic Biomarker for Bladder Cancer and Possesses Generalization Value in Pan-Cancer

Recent studies have revealed the significant role of TEA domain family member 4 (TEAD4) in the development and progression of cancer. However, the potential role of TEAD4 in the progression of bladder cancer (BC) remains to be explored. The aim of this study was to determine whether TEAD4 could serve as a pan-cancer predictor of the prognosis for BC. Based on data mined from public databases, expression levels and clinical value of TEAD4 were identified in BC and human pan-cancers. Reverse transcription-quantitative polymerase chain reaction (RT-qPCR) analysis was performed to detect the TEAD4 expression levels in BC cell lines. Gene Set Enrichment Analysis (GSEA) was carried out for functional analysis in BC, and the relationship between infiltrating immune cells and TEAD4 expression was evaluated by the CIBERSORT algorithm in BC and pan-cancer data. TEAD4 was overexpressed and associated with poor prognosis in BC and several types of cancers. GSEA and CIBERSORT algorithm suggested that various pathways including immune-related pathways were enriched in TEAD4 high expression group and several immunocytes infiltrated were correlated with the expression of TEAD4. This study revealed TEAD4 is an immune regulating-related predictor of prognosis for BC and has generalization value in pan-cancer.

Probing the Anti-Cancer Potency of Sulfated Galactans on Cholangiocarcinoma Cells Using Synchrotron FTIR Microspectroscopy, Molecular Docking, and In Vitro Studies

Sulfated galactans (SG) isolated from red alga Gracilaria fisheri have been reported to inhibit the growth of cholangiocarcinoma (CCA) cells, which was similar to the epidermal growth factor receptor (EGFR)-targeted drug, cetuximab. Herein, we studied the anti-cancer potency of SG compared to cetuximab. Biological studies demonstrated SG and cetuximab had similar inhibition mechanisms in CCA cells by down-regulating EGFR/ERK pathway, and the combined treatment induced a greater inhibition effect. The molecular docking study revealed that SG binds to the dimerization domain of EGFR, and this was confirmed by dimerization assay, which showed that SG inhibited ligand-induced EGFR dimer formation. Synchrotron FTIR microspectroscopy was employed to examine alterations in cellular macromolecules after drug treatment. The SR-FTIR-MS elicited similar spectral signatures of SG and cetuximab, pointing towards the bands of RNA/DNA, lipids, and amide I vibrations, which were inconsistent with the changes of signaling proteins in CCA cells after drug treatment. Thus, this study demonstrates the underlined anti-cancer mechanism of SG by interfering with EGFR dimerization. In addition, we reveal that FTIR signature spectra offer a useful tool for screening anti-cancer drugs’ effect.

Extracellular Matrix: Emerging Roles and Potential Therapeutic Targets for Breast Cancer

Increasing evidence shows that the extracellular matrix (ECM) is an important regulator of breast cancer (BC). The ECM comprises of highly variable and dynamic components. Compared with normal breast tissue under homeostasis, the ECM undergoes many changes in composition and organization during BC progression. Induced ECM proteins, including fibrinogen, fibronectin, hyaluronic acid, and matricellular proteins, have been identified as important components of BC metastatic cells in recent years. These proteins play major roles in BC progression, invasion, and metastasis. Importantly, several specific ECM molecules, receptors, and remodeling enzymes are involved in promoting resistance to therapeutic intervention. Additional analysis of these ECM proteins and their downstream signaling pathways may reveal promising therapeutic targets against BC. These potential drug targets may be combined with new nanoparticle technologies. This review summarizes recent advances in functional nanoparticles that target the ECM to treat BC. Accurate nanomaterials may offer a new approach to BC treatment.

Enantioselective toxicity and oxidative stress effects of acetochlor on earthworms (Eisenia fetida) by mediating the signaling pathway

Acetochlor (ACT) as a widely used chiral chloroacetamide herbicide is appropriate to evaluate the potential toxicity in soil ecosystems at enantiomeric level. The acute and subchronic toxicities of R-acetochlor (R-ACT) and S-acetochlor (S-ACT) on earthworms (Eisenia fetida) were investigated in the present study. Residual analyses showed that S-ACT degraded faster than R-ACT in artificial soil with half-lives of 16.5 and 21.7 d, respectively. Additionally, significant enantioselective acute toxicity in earthworms from between S-ACT and R-ACT (p < 0.05) was observed, and the acute toxicity of R-ACT were 1.9 and 1.5 times higher than those of S-ACT in the filter paper test and artificial soil test. The hydroxyl radical (OH^-) content, superoxide dismutase (SOD) and antioxidant enzyme catalase (CAT) activities, and cytochrome P450 content in earthworms significantly increased under the influence of ACT enantiomers; however, the acetylcholinesterase (AchE) activity was significantly inhibited after exposure to the two enantiomers. Moreover, lipid peroxidation and DNA damage were induced by ACT enantiomers. The results of transcriptome sequencing indicated that R-ACT induced a stronger oxidative stress effect than S-ACT in earthworms by mediating signaling pathways, which may be the primary reason for the enantioselective toxicity between S-ACT and R-ACT. Overall, the results demonstrated that R-ACT has a higher risk than S-ACT in the soil environment, which is important for understanding the enantioselective behavior of chloroacetamide pesticides.

Identification of Hub Genes Associated with Lung Adenocarcinoma Based on Bioinformatics Analysis

Lung adenocarcinoma (LUAD) is one of the malignant lung tumors. However, its pathology has not been fully understood. The purpose of this study is to identify the hub genes associated with LUAD by bioinformatics methods. Three gene expression datasets including GSE116959, GSE74706, and GSE85841 downloaded from the Gene Expression Omnibus (GEO) database were used in this study. The differentially expressed genes (DEGs) related to LUAD were screened by using the limma package. Gene Ontology (GO) and KEGG analysis of DEGs were carried out through the DAVID website. The protein-protein interaction (PPI) of differentially expressed genes was drawn by the STRING website, and the results were imported into Cytoscape for visualization. Then, the PPI network was analyzed by using MCODE, and the modules with a score greater than 5 were found by using cytoHubba. Finally, the GEPIA database and UALCAN database were used to verify and analyze the survival of hub genes. We identified 67 upregulated genes and 277 downregulated genes from three LUAD datasets. The results of GO analysis showed that the downregulated genes were significantly enriched in matrix adhesion and angiogenesis and upregulated differential genes were significantly enriched in cell adhesion and vascular development. KEGG pathway analysis showed that the differential genes of LUAD were significantly enriched in viral carcinogenesis and adhesion spots. The PPI network of differentially expressed genes consists of 269 nodes and 625 interactions. In addition, three modules with scores greater than 5 and seven hub genes, namely, MCM4, BIRC5, CDC20, CDC25C, FOXM1, GTSE1, and RFC4, playing an important role in the PPI network were screened out. In this study, we obtained the hub genes and pathways related to LUAD, revealing the molecular mechanism and pathogenesis of LUAD, which is helpful for the early detection of LUAD and provides a new idea for the treatment of LUAD.

The varied influences of cell adhesion and spreading on gene transfection of mesenchymal stem cells on a micropatterned substrate

Transmembrane transport of exogenous genes is widely investigated because of high demand for gene therapy. Both gene carriers and cellular conditions can affect gene transfection efficiency. Although cell morphology has been reported to affect cell functions, the influence of cell adhesion area and cell spreading area on the transfection of exogenous genes remains unclear because it is difficult to separate the individual influence of these areas during normal cell culture. In this study, micropatterns were prepared to separately control the adhesion and spreading areas of human bone marrow-derived mesenchymal stem cells (hMSCs). Transfection efficiency of the green fluorescent protein gene to hMSCs cultured on the micropatterns was compared. Cells with a larger adhesion area showed higher transfection efficiency, while cell spreading area hardly affected gene transfection efficiency. Cell adhesion area had dominant influence on gene transfection. Microparticle uptake and BrdU staining showed that the cellular uptake capacity and DNA synthesis activity increased with the increase in cell adhesion area, but were not affected by cell spreading area. The different influence of cell adhesion area and cell spreading area on gene transfection was correlated with their influence on cellular uptake capacity, DNA synthesis activity, focal adhesion formation, cytoskeletal mechanics, and mechanotransduction signal activation. The results suggest that cell adhesion area and cell spreading area had different influence on gene transfection; this finding should provide useful information for the manipulation of cell functions in gene therapy, protein modification, and cell reprogramming. STATEMENT OF SIGNIFICANCE: Cell adhesion and spreading are important morphological factors during the interaction of cells with biomaterial surfaces or interfaces. However, the predominant morphological factor that affects cellular functions such as gene transfection remains unclear. In the present study, special micropatterns were used to precisely control cell adhesion and spreading areas independently. Mesenchymal stem cells cultured on the micropatterns were transfected with the green fluorescent protein gene to compare the different influence of cell adhesion and spreading areas on gene transfection efficiency. Cell adhesion area showed dominant influence on gene transfection, while cell spreading area did not affect gene transfection. The dominant influence of cell adhesion area could be explained by cellular uptake capacity and DNA synthesis activity through the formation of FAs, cytoskeletal mechanics, and YAP/TAZ nuclear localization. The results provide new insights of correlation between cell morphology and cellular functions for designing functional biomaterials.

From Oncogenic Signaling Pathways to Single-Cell Sequencing of Immune Cells: Changing the Landscape of Cancer Immunotherapy

Over the past decade, there have been remarkable advances in understanding the signaling pathways involved in cancer development. It is well-established that cancer is caused by the dysregulation of cellular pathways involved in proliferation, cell cycle, apoptosis, cell metabolism, migration, cell polarity, and differentiation. Besides, growing evidence indicates that extracellular matrix signaling, cell surface proteoglycans, and angiogenesis can contribute to cancer development. Given the genetic instability and vast intra-tumoral heterogeneity revealed by the single-cell sequencing of tumoral cells, the current approaches cannot eliminate the mutating cancer cells. Besides, the polyclonal expansion of tumor-infiltrated lymphocytes in response to tumoral neoantigens cannot elicit anti-tumoral immune responses due to the immunosuppressive tumor microenvironment. Nevertheless, the data from the single-cell sequencing of immune cells can provide valuable insights regarding the expression of inhibitory immune checkpoints/related signaling factors in immune cells, which can be used to select immune checkpoint inhibitors and adjust their dosage. Indeed, the integration of the data obtained from the single-cell sequencing of immune cells with immune checkpoint inhibitors can increase the response rate of immune checkpoint inhibitors, decrease the immune-related adverse events, and facilitate tumoral cell elimination. This study aims to review key pathways involved in tumor development and shed light on single-cell sequencing. It also intends to address the shortcomings of immune checkpoint inhibitors, i.e., their varied response rates among cancer patients and increased risk of autoimmunity development, via applying the data from the single-cell sequencing of immune cells.

A novel focal adhesion related gene signature for prognostic prediction in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is a highly heterogeneous disease. Reduced expression of focal adhesion is considered as an important prerequisite for tumor cell invasion and metastasis. However, the prognostic value of focal adhesion related genes in HCC remains to be further determined. In this study, RNA expression profiles were downloaded from public databases. A five focal adhesion related gene signature model was established by the least absolute shrinkage and selection operator Cox regression analysis, which categorized patients into high- and low-risk groups. Multivariate Cox regression analysis showed that the risk score was an independent predictor for overall survival. Single-sample gene set enrichment analysis revealed that immune status was different between the two risk groups, and tumor-related pathways were enriched in high-risk group. The risk score was significantly associated with tumor grade, tumor stage, immune scores, and immune infiltrate types. Pearson correlation showed that the expression level of prognostic genes was associated with anti-tumor drug sensitivity. Besides, the mRNA and protein expression of prognostic genes was significantly different between HCC tissues and adjacent non-tumorous tissues in our separate cohort. Taken together, a novel focal adhesion related gene signature can be used for prognostic prediction in HCC, which may be a therapeutic alternative.

Prognostic Significance and Tumor Immune Microenvironment Heterogenicity of m5C RNA Methylation Regulators in Triple-Negative Breast Cancer

Purpose

The m5C RNA methylation regulators are closely related to tumor proliferation, occurrence, and metastasis. This study aimed to investigate the gene expression, clinicopathological characteristics, and prognostic value of m5C regulators in triple-negative breast cancer (TNBC) and their correlation with the tumor immune microenvironment (TIM).

Methods

The TNBC data, Luminal BC data and HER2 positive BC data set were obtained from The Cancer Genome Atlas and Gene Expression Omnibus, and 11 m5C RNA methylation regulators were analyzed. Univariate Cox regression and the least absolute shrinkage and selection operator regression models were used to develop a prognostic risk signature. The UALCAN and cBioportal databases were used to analyze the gene characteristics and gene alteration frequency of prognosis-related m5C RNA methylation regulators. Gene set enrichment analysis was used to analyze cellular pathways enriched by prognostic factors. The Tumor Immune Single Cell Hub (TISCH) and Timer online databases were used to explore the relationship between prognosis-related genes and the TIM.

Results

Most of the 11 m5C RNA methylation regulators were differentially expressed in TNBC and normal samples. The prognostic risk signature showed good reliability and an independent prognostic value. Prognosis-related gene mutations were mainly amplified. Concurrently, the NOP2/Sun domain family member 2 (NSUN2) upregulation was closely related to spliceosome, RNA degradation, cell cycle signaling pathways, and RNA polymerase. Meanwhile, NSUN6 downregulation was related to extracellular matrix receptor interaction, metabolism, and cell adhesion. Analysis of the TISCH and Timer databases showed that prognosis-related genes affected the TIM, and the subtypes of immune-infiltrating cells differed between NSUN2 and NSUN6.

Conclusion

Regulatory factors of m5C RNA methylation can predict the clinical prognostic risk of TNBC patients and affect tumor development and the TIM. Thus, they have the potential to be a novel prognostic marker of TNBC, providing clues for understanding the RNA epigenetic modification of TNBC.