Prediction of key genes and pathways involved in trastuzumab-resistant gastric cancer

KIF26B and CREB3L1 Derived from Immunoscore Could Inhibit the Progression of Ovarian Cancer

Background

Ovarian cancer (OV) is characteristic of high incidence rate and fatality rate in the malignant tumors of female reproductive system. Researches on pathogenesis and therapeutic targets for OV need to be continued. This study mainly analyzed the immune-related pathogenesis and discovered the key immunotherapy targets for OV.

Methods

WGCNA was used for excavating hub gene modules and hub genes related to the immunity of OV. Enrichment analysis was aimed to analyze the related pathways of hub gene modules. Biological experiments were used for exploring the effect of hub genes on SKOV3 cells.

Results

We identified two hub gene modules related to the immunoscore of OV and found that these genes in the modules were related to the extracellular matrix and viral infections. At the same time, we also discovered six hub genes related to the immunity of OV. Among them, KIF26B and CREB3L1 can affect the proliferation, migration, and invasion of SKOV3 cells by the Wnt/β-catenin pathway.

Conclusions

The local infection or inflammation of ovarian may affect the immunity of OV. KIF26B and CREB3L1 are expected to be potential targets for the immunotherapy of OV.

OASL knockdown inhibits the progression of stomach adenocarcinoma by regulating the mTORC1 signaling pathway

The present study investigated the effects of 2'-5' oligoadenylate synthetase-like (OASL) on the biological functions of stomach adenocarcinoma (STAD) cells and tumor formation in nude mice. The differential expression levels of OASL in the different cancer types from TCGA dataset were analyzed using gene expression profiling interactive analysis. Overall survival and the receiver operating characteristic were analyzed using the KM plotter and R, respectively. Furthermore, OASL expression and its effects on the biological functions of STAD cells were detected. The possible upstream transcription factors of OASL were predicted using JASPAR. The downstream signaling pathways of OASL were analyzed using GSEA. Tumor formation experiments were performed to evaluate the effect of OASL on tumor formation in nude mice. The results showed that OASL was highly expressed in STAD tissues and cell lines. OASL knockdown markedly inhibited cell viability, proliferation, migration, and invasion and accelerated STAD cell apoptosis. Conversely, OASL overexpression had the opposite effect on STAD cells. JASPAR analysis revealed that STAT1 is an upstream transcription factor of OASL. Furthermore, GSEA showed that OASL activated the mTORC1 signaling pathway in STAD. The protein expression levels of p-mTOR and p-RPS6KB1 were suppressed by OASL knockdown and promoted by OASL overexpression. The mTOR inhibitor, rapamycin, markedly reversed the effect of OASL overexpression on STAD cells. Additionally, OASL promoted tumor formation and increased tumor weight and volume in vivo. In conclusion, OASL knockdown suppressed the proliferation, migration, invasion, and tumor formation of STAD cells by inhibiting the mTOR signaling pathway.

Prediction of the Endocrine disruption profile of fluorinated biphenyls and analogues: An in silico study

Fluorinated biphenyls and their analogues (FBAs) are considered new persistent organic pollutants, but their endocrine-disrupting effects are still unknown. To fill this gap, the binding probability of 44 FBAs to different nuclear hormone receptors (NHRs) was predicted using Endocrine Disruptome. And molecular similarity and network toxicology analysis were used to strengthen the docking screening. The docking results showed that FBAs could have high binding potential for various NHRs, such as estrogen receptors β antagonism (ERβ an), liver X receptors α (LXRα), estrogen receptors α (ERα), and liver X receptors β (LXRβ). The similarity analysis found that the degree of overlap of the NHR repertoire was related to the Tanimoto coefficient of FBAs. Network toxicology verified a part of docking screening results and identified endocrine-disrupting pathways worthy of attention. This study found out potential endocrine-disrupting FBAs and their vulnerable, and developed a workflow that would leverage in silico approaches including molecular docking, similarity, and network toxicology for risk prioritization of potential endocrine-disrupting compounds.

Systematic Analysis of Molecular Subtypes Based on the Expression Profile of Immune-Related Genes in Pancreatic Cancer

Immunotherapy has a good therapeutic effect and provides a new approach for cancer treatment. However, only limited studies have focused on the use of molecular typing to construct an immune characteristic index for gene expression in pancreatic adenocarcinoma (PAAD) and to assess the effectiveness of immunotherapy in patients with PAAD. Clinical follow-up data and gene expression profile of PAAD patients were retrieved from The Cancer Genome Atlas (TCGA) database. Based on 184 immune features, molecular subtypes of pancreatic cancer were found by the "ConsensusClusterPlus" package, and the association between clinical features and immune cell subtype distribution was analysed. In addition, the relationship between the proportion of immune subtypes and the expression of immune checkpoints was analysed. The CIBERSORT algorithm was introduced to evaluate the immune scores of different molecular subtypes. We used the tumor immune dysfunction and exclusion (TIDE) algorithm to assess the potential clinical effect of immunotherapy interventions on single-molecule subtypes. In addition, the oxidative stress index was constructed by linear discriminant analysis DNA (LDA), and weighted correlation network analysis was performed to identify the core module of the index and its characteristic genes. Expression of hub genes was verified by immunohistochemical analysis in an independent PAAD cohort. Pancreatic cancer is divided into three molecular subtypes (IS1, IS2, and IS3), with significant differences in prognosis between multiple cohorts. Expression of immune checkpoint-associated genes was significantly reduced in IS3 and higher in IS1 and IS2, suggesting that the three subgroups have different responsiveness to immunotherapy interventions. The results of the CIBERSORT analysis showed that IS1 exhibited the highest levels of immune infiltration, whereas the results of the TIDE analysis showed that the T-cell dysfunction score of IS1 was higher than that of IS2 and IS3. Furthermore, IS3 was found to be more sensitive to 5-FU and to have a higher immune signature index than IS1 and IS2. Based on WGCNA analysis, 10 potential gene markers were identified, and their expression at the protein level was verified by immunohistochemical analysis. Specific molecular expression patterns in pancreatic cancer can predict the efficacy of immunotherapy and influence the prognosis of patients.

Overexpression of OAS1 Is Correlated With Poor Prognosis in Pancreatic Cancer

Background

OAS1 expression in pancreatic cancer has been confirmed by many studies. However, the prognostic value and mechanism of OAS1 in pancreatic cancer have not been analyzed.

Methods

The RNA-seq in pancreatic cancer were obtained by UCSC XENA and GEO database. In addition, immunohistochemical validation and analysis were performed using samples from the 900th hospital. The prognosis of OAS1 was evaluated by timeROC package, Cox regression analysis, and Kaplan-Meier survival curves. Then, the main functional and biological signaling pathways enrichment and its relationship with the abundance of immune cells were analyzed by bioinformatics.

Results

OAS1 was highly expressed in pancreatic cancer compared with normal pancreatic tissue. High OAS1 expression was associated with poor overall survival (p<0.05). The OAS1 was significantly correlated to TNM staging (p=0.014). The timeROC analysis showed that the AUC of OAS1 was 0.734 for 3-year OS. In addition, the expression of OAS1 was significantly correlated with the abundance of a variety of immune markers. GSEA showed that enhanced signaling pathways associated with OAS1 include Apoptosis, Notch signaling pathway, and P53 signaling pathway.

Conclusions

OAS1 is a valuable prognostic factor in pancreatic cancer. Moreover, it may be a potential immunotherapeutic target.

Ontologies and Knowledge Graphs in Oncology Research

The complexity of cancer research stems from leaning on several biomedical disciplines for relevant sources of data, many of which are complex in their own right. A holistic view of cancer—which is critical for precision medicine approaches—hinges on integrating a variety of heterogeneous data sources under a cohesive knowledge model, a role which biomedical ontologies can fill. This study reviews the application of ontologies and knowledge graphs in cancer research. In total, our review encompasses 141 published works, which we categorized under 14 hierarchical categories according to their usage of ontologies and knowledge graphs. We also review the most commonly used ontologies and newly developed ones. Our review highlights the growing traction of ontologies in biomedical research in general, and cancer research in particular. Ontologies enable data accessibility, interoperability and integration, support data analysis, facilitate data interpretation and data mining, and more recently, with the emergence of the knowledge graph paradigm, support the application of Artificial Intelligence methods to unlock new knowledge from a holistic view of the available large volumes of heterogeneous data.

Role of helical structure and dynamics in oligoadenylate synthetase 1 (OAS1) mismatch tolerance and activation by short dsRNAs

The 2'-5'-oligoadenylate synthetases (OAS) are innate immune sensors of cytosolic double-stranded RNA (dsRNA) that play a critical role in limiting viral infection. How these proteins are able to avoid aberrant activation by cellular RNAs is not fully understood, but adenosine-to-inosine (A-to-I) editing has been proposed to limit accumulation of endogenous RNAs that might otherwise cause stimulation of the OAS/RNase L pathway. Here, we aim to uncover whether and how such sequence modifications can restrict the ability of short, defined dsRNAs to activate the single-domain form of OAS, OAS1. Unexpectedly, we find that all tested inosine-containing dsRNAs have an increased capacity to activate OAS1, whether in a destabilizing (I•U) or standard Watson-Crick-like base pairing (I-C) context. Additional variants with strongly destabilizing A•C mismatches or stabilizing G-C pairs also exhibit increased capacity to activate OAS1, eliminating helical stability as a factor in the relative ability of the dsRNAs to activate OAS1. Using thermal difference spectra and molecular dynamics simulations, we identify both increased helical dynamics and specific local changes in helical structure as important factors in the capacity of short dsRNAs to activate OAS1. These helical features may facilitate more ready adoption of the distorted OAS1-bound conformation or stabilize important structures to predispose the dsRNA for optimal binding and activation of OAS1. These studies thus reveal the molecular basis for the greater capacity of some short dsRNAs to activate OAS1 in a sequence-independent manner.

A comprehensive gene expression profile analysis of prostate cancer cells resistant to paclitaxel and the potent target to reverse resistance

Background: Paclitaxel resistance is the major clinical obstacle in the chemotherapy of prostate cancer (PCa), but the resistant mechanism is less investigated.Purpose: To establish two paclitaxel-resistant PCa cells, provide a comprehensive gene expression profile analysis of resistant cells and the potential target to reverse resistance.Methods: Two Paclitaxel-resistant PCa cells (PC3/PR, LNcap/PR) were established by gradually increasing drug concentration. MTT and transwell assays were performed to detect drug sensitivity, cell proliferation and migration abilities. RNA-Sequencing (RNA-seq) and bioinformatic analyses were performed to identify abnormally expressed genes (AEGs) in resistant cells, and annotate the biological functions of AEGs. The role of the candidate AEG, TLR-4, on the resistant phenotypes was further investigated.Results: The resistance index of resistant cells was 2-3, and they showed a slower proliferation and increased migration ability. 4741 AEGs were screened out (Log2fold change absolute: log2FC(abs)  > 1) in the resistant cells, and they were enriched in 2'-5'-oligoadenylate synthetase activity and chemical carcinogenesis. A number of AEGs, CCND2, IGFBP3, FOS, SHH, ZEB2, and members of FGF, FGFR and WNT families were also identified to be involved in cancer- and resistant phenotype-related processes. Finally, TLR-4 was validated significantly increased in resistant cells, and knockdown of TLR-4 increased drug-sensitivity, inhibited the proliferation and migration abilities.Conclusions: The study provided a comprehensive gene expression profile of paclitaxel-resistant PCa cells, and TLR-4 could be a potential target to reverse paclitaxel resistance.

Engineering an anti-HER2 biparatopic antibody with a multimodal mechanism of action

The receptor tyrosine kinase HER2 acts as oncogenic driver in numerous cancers. Usually, the gene is amplified, resulting in receptor overexpression, massively increased signaling and unchecked proliferation. However, tumors become frequently addicted to oncogenes and hence are druggable by targeted interventions. Here, we design an anti-HER2 biparatopic and tetravalent IgG fusion with a multimodal mechanism of action. The molecule first induces HER2 clustering into inactive complexes, evidenced by reduced mobility of surface HER2. However, in contrast to our earlier binders based on DARPins, clusters of HER2 are thereafter robustly internalized and quantitatively degraded. This multimodal mechanism of action is found only in few of the tetravalent constructs investigated, which must target specific epitopes on HER2 in a defined geometric arrangement. The inhibitory effect of our antibody as single agent surpasses the combination of trastuzumab and pertuzumab as well as its parental mAbs in vitro and it is effective in a xenograft model.

Detailed Molecular Mechanism and Potential Drugs for COL1A1 in Carboplatin-Resistant Ovarian Cancer

Carboplatin resistance in ovarian cancer (OV) is a major medical problem. Thus, there is an urgent need to find novel therapeutic targets to improve the prognosis of patients with carboplatin-resistant OV. Accumulating evidence indicates that the gene COL1A1 (collagen type I alpha 1 chain) has an important role in chemoresistance and could be a therapeutic target. However, there have been no reports about the role of COL1A1 in carboplatin-resistant OV. This study aimed to establish the detailed molecular mechanism of COL1A1 and predict potential drugs for its treatment. We found that COL1A1 had a pivotal role in carboplatin resistance in OV by weighted gene correlation network analysis and survival analysis. Moreover, we constructed a competing endogenous RNA network (LINC00052/SMCR5-miR-98-COL1A1) based on multi-omics data and experiments to explore the upstream regulatory mechanisms of COL1A1. Two key pathways involving COL1A1 in carboplatin resistance were identified by co-expression analysis and pathway enrichment: the "ECM-receptor interaction" and "focal adhesion" Kyoto Encyclopedia of Genes and Genomes pathways. Furthermore, combining these results with those of cell viability assays, we proposed that ZINC000085537017 and quercetin were potential drugs for COL1A1 based on virtual screening and the TCMSP database, respectively. These results might help to improve the outcome of OV in the future.

Prognostic characterization of OAS1/OAS2/OAS3/OASL in breast cancer

BACKGROUND

Prognostic biomarkers remain a focus in breast cancer during last decades. More reliable predictors to adequately characterize the prognosis of breast cancer are essential. The 2'-5'-oligoadenylate synthetases (OAS), composing of OAS1, OAS2, OAS3, and OAS-like (OASL), are interferon (IFN)-induced antiviral enzymes, with their prognostic roles remain to be characterized.

METHODS

Prognostic values of OAS family members were assessed by multiple public available resources.

RESULTS

High mRNA expression of OAS1 and OAS3 were correlated with worse prognosis for all breast cancer patients, whereas OAS2 was associated with favorable prognosis. The prognostic values of OAS family in different clinicopathologic subtypes were also characterized. In DNA methylation level, cg12560128 in OAS2, cg06800840 and cg26328872 in OASL showed significant prognostic values. The mRNA expression of OAS members signature in high/low risk overall survival groups was opposite to the high/low risk recurrence free survival groups. Neutrophil cell exhibited highest correlation with all OAS members in tumor immune infiltrating estimation.

CONCLUSIONS

This study provided new insight into the prognostic roles of OAS in breast cancer with potential mechanistic values.

Identification and Validation of an Individualized EMT-Related Prognostic Risk Score Formula in Gastric Adenocarcinoma Patients

Background

The epithelial-mesenchymal transition (EMT) is a pivotal process for fibrotic disease, embryonic development, and wound healing. Moreover, some evidence has proven that the disorder of EMT also plays an important role in carcinogenesis, especially invasion and metastasis of various tumors (Ritchie et al., 2015). Additionally, gastric adenocarcinoma (GAC) is a common gastrointestinal malignancy which is the fourth most commonly diagnosed tumor. Our study is aimed at identifying the prognostic value of EMT-related genes in gastric adenocarcinoma.

Methods

Firstly, high-throughput and clinical data were downloaded from The Cancer Genome Atlas (TCGA) database. 99 differentially expressed EMT-related genes (ERGs) were obtained in these gastric adenocarcinoma data. Secondly, GO and KEGG enrichment analyses show that EMT may promote gastric carcinogenesis. Next, 10 ERGs associated with prognosis of gastric adenocarcinoma patients are screened out by univariate Cox regression, and 6 pivotal prognostic ERGs (MMP8, MMP11, TFDP3, MYB, F2, and CNTN1) are identified through multivariate Cox regression. These 6 genes are confirmed with significant prognostic value in gastric adenocarcinoma through overall survival (OS) analysis. Finally, a risk score formula is constructed and tested in another gastric adenocarcinoma cohort from GEO.

Results

99 differentially expressed EMT-related genes (ERGs) and their enriched pathways are identified. 10 ERGs are strongly related to the prognosis of GAC patients. A risk score formula of 6 prognosis-related ERGs used to predict the prognosis of gastric adenocarcinoma patients is identified and tested (risk score = 0.448115∗expression value of MMP8 + 0.378892∗expression value of MMP11 − 0.3226∗expression value of MYB + 1.322812∗expression value of TFDP3 + 0.325063∗expression value of F2 + 0.334197∗expression value of CNTN1).

Conclusion

This study provides a potential prognostic signature for predicting prognosis of gastric adenocarcinoma patients and molecular insights of EMT in gastric adenocarcinoma, and the formula focusing on the prognosis of gastric adenocarcinoma can be effective.

Early Modulation of Circulating MicroRNAs Levels in HER2-Positive Breast Cancer Patients Treated with Trastuzumab-Based Neoadjuvant Therapy

Circulating microRNA (ct-miRNAs) are able to identify patients with differential response to HER2-targeted therapy. However, their dynamics are largely unknown. We assessed 752 miRNAs from 52 NeoALTTO patients with plasma pairs prior and two weeks after trastuzumab. Increased levels of ct-miR-148a-3p and ct-miR-374a-5p were significantly associated with pathological complete response (pCR) (p = 0.008 and 0.048, respectively). At a threshold ≥ the upper limit of the 95%CI of the mean difference, pCR resulted 45% (95%CI 24%–68%), and 44% (95%CI 22%–69%) for ct-miR-148a-3p and ct-miR-374a-5p, respectively. Notably, ct-miR-148a-3p retained its predictive value (OR 3.42, 95%CI 1.23–9.46, p = 0.018) in bivariate analysis along with estrogen receptor status. Combined information from ct-miR-148a-3p and ct-miR140-5p, which we previously reported to identify trastuzumab-responsive patients, resulted in greater predictive capability over each other, with pCR of 54% (95%CI 25%–81%) and 0% (95%CI 0%–31%) in ct-miR-148a/ct-miR-140-5p high/present and low/absent, respectively. GO and KEGG analyses showed common enriched terms between the targets of these ct-miRNAs, including cell metabolism regulation, AMPK and MAPK signaling, and HCC progression. In conclusion, early modulated ct-miR-148-3p may inform on the functional processes underlying treatment response, integrate the information from already available predictive biomarkers, and identify patients likely to respond to single agent trastuzumab-based neoadjuvant therapy.

Identification of potential microRNAs and their targets in promoting gefitinib resistance by integrative network analysis

Background

Non-small cell lung cancer (NSCLC) accounts for about 80-85% of lung cancers. Epidermal growth factor receptor tyrosine kinase inhibitors (EGFR-TKIs) such as gefitinib are considered the best choice for first-line treatment for the patients with NSCLC harboring EGFR-activating alterations. Nonetheless, 10-30% of patients may not obtain an objective response and may also experience rapid progression. The aim of our research, based on the integrative bioinformatics review, was to identify the possible miRNAs involved in gefitinib resistance.

Method

A gefitinib-resistant network composed of 15 miRNAs and 34 targets were constructed by using the bioinformatics analyses of three microarray datasets. Of these miRNAs, effects of miR-342-3p on gefitinib resistance were investigated on a gefitinib-resistant cell model (A549/GR and PC/GR cells).

Results

We reported that over-expression of miR-342-3p could significantly increase the resistance to gefitinib of A549/GR and PC9/GR cells and vice versa. Then, we recognized CPA4 as a target of hsa-miR-342-3p by a luciferase reporter assay. The increase in hsa-miR-342-3p levels led to a significant reduction in CPA4 protein expression. However, the opposite results were observed upon miR-342-3p knockdown. Finally, we found that enforced CPA4 expression partially reversed miR-342-3p effects in A549/GR cells.

Conclusions

Collectively, these findings suggest that the upregulation of miR-342-3p contributes to gefitinib resistance by targeting CPA4, which may serve as a potential treatment option to overcome gefitinib resistance in patients with NSCLC.