Integrated bioinformatics analysis of the crucial candidate genes and pathways associated with glucocorticoid resistance in acute lymphoblastic leukemia

Synergistic combination of isogarcinol isolated from edible fruits of Garcinia multiflora and dexamethasone to overcome leukemia glucocorticoid resistance

Isogarcinol (ISO), a cytotoxic polycyclic polyprenylated acylphloroglucinol isolated from the edible fruits of Garcinia multiflora. However, synergistic combination of ISO and dexamethasone (DEX) to overcome leukemia glucocorticoid resistance has never been investigated. Therefore, in this study, the effects of ISO in combination with DEX was conducted on leukemia in vivo and glucocorticoid resistance in vitro. As a result, the combination of the two compounds could efficiently inhibit leukemia progression in mice and reverse DEX resistance in acute lymphoblastic leukemia (ALL) Jurkat cells. Significantly, our findings indicated that c-Myc may be a potential target of ISO, as it is involved in cell cycle arrest and apoptosis by the combination of ISO and DEX in Jurkat cells. Furthermore, western blot analysis revealed that ISO and DEX inhibits the PI3K/Akt/mTOR signaling pathway and promotes the nuclear translocation of glucocorticoid receptor (GR), which activates target genes NR3C1 and TSC22D3, leading to apoptosis in Jurkat cells. Hence, our results suggest that ISO, as a safe and effective food-derived agent, can enhance the anti-leukemia effects of DEX.

Identification of Hub Genes Associated with Resistance to Prednisolone in Acute Lymphoblastic Leukemia Based on Weighted Gene Co-expression Network Analysis

Resistance against glucocorticoids which are used to reduce inflammation and treatment of a number of diseases, including leukemia, is known as the first stage of treatment failure in acute lymphoblastic leukemia. Since these drugs are the essential components of chemotherapy regimens for ALL and play an important role in stop of cell growth and induction of apoptosis, it is important to identify genes and the molecular mechanism that may affect glucocorticoid resistance. In this study, we used the GSE66705 dataset and weighted gene co-expression network analysis (WGCNA) to identify modules that correlated more strongly with prednisolone resistance in type B lymphoblastic leukemia patients. The PPI network was built using the DEGs key modules and the STRING database. Finally, we used the overlapping data to identify hub genes. out of a total of 12 identified modules by WGCNA, the blue module was find to have the most statistically significant correlation with prednisolone resistance and Nine genes including SOD1, CD82, FLT3, GART, HPRT1, ITSN1, TIAM1, MRPS6, MYC were recognized as hub genes Whose expression changes can be associated with prednisolone resistance. Enrichment analysis based on the MsigDB repository showed that the altered expressed genes of the blue module were mainly enriched in IL2_STAT5, KRAS, MTORC1, and IL6-JAK-STAT3 pathways, and their expression changes can be related to cell proliferation and survival. The analysis performed by the WGCNA method introduced new genes. The role of some of these genes was previously reported in the resistance to chemotherapy in other diseases. This can be used as clues to detect treatment-resistant (drug-resistant) cases in the early stages of diseases.

Regulation and Functions of α6-Integrin (CD49f) in Cancer Biology

Over the past decades, our knowledge of integrins has evolved from being understood as simple cell surface adhesion molecules to receptors that have a complex range of intracellular and extracellular functions, such as delivering chemical and mechanical signals to cells. Consequently, they actively control cellular proliferation, differentiation, and apoptosis. Dysregulation of integrin signaling is a major factor in the development and progression of many tumors. Many reviews have covered the broader integrin family in molecular and cellular studies and its roles in diseases. Nevertheless, further understanding of the mechanisms specific to an individual subunit of different heterodimers is more useful. Thus, we describe the current understanding of and exploratory investigations on the α6-integrin subunit (CD49f, VLA6; encoded by the gene itga6) in normal and cancer cells. The roles of ITGA6 in cell adhesion, stemness, metastasis, angiogenesis, and drug resistance, and as a diagnosis biomarker, are discussed. The role of ITGA6 differs based on several features, such as cell background, cancer type, and post-transcriptional alterations. In addition, exosomal ITGA6 also implies metastatic organotropism. The importance of ITGA6 in the progression of a number of cancers, including hematological malignancies, suggests its potential usage as a novel prognostic or diagnostic marker and useful therapeutic target for better clinical outcomes.

A Predictive Disease Risk Model for Ankylosing Spondylitis: Based on Integrated Bioinformatic Analysis and Identification of Potential Biomarkers Most Related to Immunity

Background

The pathogenesis of ankylosing spondylitis (AS) is still not clear, and immune-related genes have not been systematically explored in AS. The purpose of this paper was to identify the potential early biomarkers most related to immunity in AS and develop a predictive disease risk model with bioinformatic methods and the Gene Expression Omnibus database (GEO) to improve diagnostic and therapeutic efficiency.

Methods

To identify differentially expressed genes and create a gene coexpression network between AS and healthy samples, we downloaded the AS-related datasets GSE25101 and GSE73754 from the GEO database and employed weighted gene coexpression network analysis (WGCNA). We used the GSVA, GSEABase, limma, ggpubr, and reshape2 packages to score immune data and investigated the links between immune cells and immunological functions by using single-sample gene set enrichment analysis (ssGSEA). The value of the core gene set and constructed model for early AS diagnosis was investigated by using receiver operating characteristic (ROC) curve analysis.

Results

Biological function and immune score analyses identified central genes related to immunity, key immune cells, key related pathways, gene modules, and the coexpression network in AS. Granulysin (GNLY), Granulysin (GZMK), CX3CR1, IL2RB, dysferlin (DYSF), and S100A12 may participate in AS development through NK cells, CD8⁺ T cells, Th1 cells, and other immune cells and represent potential biomarkers for the early diagnosis of AS occurrence and progression. Furthermore, the T cell coinhibitory pathway may be involved in AS pathogenesis.

Conclusion

The AS disease risk model constructed based on immune-related genes can guide clinical diagnosis and treatment and may help in the development of personalized immunotherapy.

The regulatory role of PDE4B in the progression of inflammatory function study

Inflammation is a response of the body to external stimuli (eg. chemical irritants, bacteria, viruses, etc.), and when the stimuli are persistent, they tend to trigger chronic inflammation. The presence of chronic inflammation is an important component of the tumor microenvironment produced by a variety of inflammatory cells (eg. macrophages, neutrophils, leukocytes, etc.). The relationship between chronic inflammation and cancer development has been widely accepted, and chronic inflammation has been associated with the development of many cancers, including chronic bronchitis and lung cancer, cystitis inducing bladder cancer. Moreover, chronic colorectitis is more likely to develop into colorectal cancer. Therefore, the specific relationship and cellular mechanisms between inflammation and cancer are a hot topic of research. Recent studies have identified phosphodiesterase 4B (PDE4B), a member of the phosphodiesterase (PDEs) protein family, as a major cyclic AMP (cAMP) metabolizing enzyme in inflammatory cells, and the therapeutic role of PDE4B as chronic inflammation, cancer. In this review, we will present the tumors associated with chronic inflammation, and PDE4B potential clinical application.

Population genetics of PDE4B (Phosphodiesterase-4B) in neglected native americans: implications for cancer pharmacogenetics

PDE4B (phosphodiesterase‐4B) has an important role in cancer and in pharmacology of some disorders, such as inflammatory diseases. Remarkably in Native Americans, PDE4B variants are associated with acute lymphoblastic leukemia (ALL) relapse, as this gene modulates sensitivity of glucocorticoids used in ALL chemotherapy. PDE4B allele rs6683977.G, associated with genomic regions of Native American origin in US‐Hispanics (admixed among Native Americans, Europeans, and Africans), increases ALL relapse risk, contributing to an association between Native American ancestry and ALL relapse that disappeared with an extra‐phase of chemotherapy. This result insinuates that indigenous populations along the Americas may have high frequencies of rs6683977.G, but this has never been corroborated. We studied ancestry and PDE4B diversity in 951 healthy individuals from nine Latin American populations. In non‐admixed Native American populations rs6683977.G has frequencies greater than 90%, is in linkage disequilibrium with other ALL relapse associated and regulatory variants in PDE4B‐intron‐7, conforming haplotypes showing their highest worldwide frequencies in Native Americans (>0.82). Our findings inform the discussion on the pertinence of an extra‐phase of chemotherapy in Native American populations, and exemplifies how knowledge generated in US‐Hispanics is relevant for their even more neglected and vulnerable Native American ancestors along the American continent.

Prognostic Relevance of Expression of EMP1, CASP1, and NLRP3 Genes in Pediatric B-Lineage Acute Lymphoblastic Leukemia

Glucocorticoid (GC), such as prednisolone, is an essential component of multidrug chemotherapy regimen for pediatric acute lymphoblastic leukemia (ALL). Resistance to GC in leukemia cells is associated with disease progression and poor prognosis. Despite the extensive use of GC for many years, molecular mechanisms underlying its resistance in ALL have not been fully uncovered. Recent studies have shown a potential role of EMP1, CASP1, and NLRP3 genes in prednisolone response. In this study on 148 pediatric B-ALL patients, we studied these three genes to assess their association with prednisolone response measured by day 8 blast count after 7 days of induction therapy with prednisolone. Intriguingly, ALL samples exhibited higher expression of EMP1 along with a low expression of CASP1 and NLRP3 compared to disease free normal bone marrow collected from patients with solid tumors. Among the three analyzed genes, only EMP1 was found to be overexpressed in prednisolone poor responders (p=0.015). Further, a comparison of gene expression between cytogenetic subtypes revealed higher expression of EMP1 in BCR-ABL subtype. Expression of EMP1 in multiple gene expression datasets was used for gene set enrichment analysis, which revealed TNF-α, IL-2-STAT5 signaling, inflammatory responses and hypoxia as the major positively associated pathways and E2F targets as negatively associated pathways. Interestingly, the clinical remission rate was higher in CASP1 high patients (p=0.048). In univariate survival analysis, higher EMP1 expression was associated with poor prognostic measures while higher expression of NLRP3 and CASP1 was associated with better prognostic measures in our data. Further, multivariate analysis revealed an independent association of high CASP1 and NLRP3 with a better prognosis. This study strengthens the available evidence that mRNA expression of EMP1, CASP1, and NLRP3 may serve as potential biomarkers for risk stratification of pediatric B-ALL patients.

Identification of Key Genes Associated with Changes in the Host Response to Severe Burn Shock: A Bioinformatics Analysis with Data from the Gene Expression Omnibus (GEO) Database

Background

Patients with severe burns continue to display a high mortality rate during the initial shock period. The precise molecular mechanism underlying the change in host response during severe burn shock remains unknown. This study aimed to identify key genes leading to the change in host response during burn shock.

Methods

The GSE77791 dataset, which was utilized in a previous study that compared hydrocortisone administration to placebo (NaCl 0.9%) in the inflammatory reaction of severe burn shock, was downloaded from the Gene Expression Omnibus (GEO) database and analyzed to identify the differentially expressed genes (DEGs). Functional enrichment analyses of Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were performed. The protein–protein interaction (PPI) network of DEGs was constructed using the Search Tool for the Retrieval of Interacting Genes (STRING) database and then visualized in Cytoscape. In addition, important modules in this network were selected using the Molecular Complex Detection (MCODE) algorithm, and hub genes were identified in cytoHubba, a Cytoscape plugin.

Results

A total of 1059 DEGs (508 downregulated genes and 551 upregulated genes) were identified from the dataset. The DEGs enriched in GO terms and KEGG pathways were related to immune response. The PPI network contained 439 nodes and 2430 protein pairs. Finally, important modules and hub genes were identified using the different Cytoscape plugins. The key genes in burn shock were identified as arginase 1 (ARG1), cytoskeleton-associated protein (CKAP4), complement C3a receptor (C3AR1), neutrophil elastase (ELANE), gamma-glutamyl hydrolase (GGH), orosomucoid (ORM1), and quiescin sulfhydryl (QSOX1).

Conclusion

The DEGs, functional terms and pathways, and hub genes identified in the present study can help shed light on the molecular mechanism underlying the changes in host response during burn shock and provide potential targets for early detection and treatment of burn shock.

Gene regulation network analysis reveals core genes associated with survival in glioblastoma multiforme

Glioblastoma multiforme (GBM) is a very serious mortality of central nervous system cancer. The microarray data from GSE2223, GSE4058, GSE4290, GSE13276, GSE68848 and GSE70231 (389 GBM tumour and 67 normal tissues) and the RNA-seq data from TCGA-GBM dataset (169 GBM and five normal samples) were chosen to find differentially expressed genes (DEGs). RRA (Robust rank aggregation) method was used to integrate seven datasets and calculate 133 DEGs (82 up-regulated and 51 down-regulated genes). Subsequently, through the PPI (protein-protein interaction) network and MCODE/ cytoHubba methods, we finally filtered out ten hub genes, including FOXM1, CDK4, TOP2A, RRM2, MYBL2, MCM2, CDC20, CCNB2, MYC and EZH2, from the whole network. Functional enrichment analyses of DEGs were conducted to show that these hub genes were enriched in various cancer-related functions and pathways significantly. We also selected CCNB2, CDC20 and MYBL2 as core biomarkers, and further validated them in CGGA, HPA and CCLE database, suggesting that these three core hub genes may be involved in the origin of GBM. All these potential biomarkers for GBM might be helpful for illustrating the important role of molecular mechanisms of tumorigenesis in the diagnosis, prognosis and targeted therapy of GBM cancer.

Integrated bioinformatics analysis of the crucial candidate genes and pathways associated with glucocorticoid resistance in acute lymphoblastic leukemia

Glucocorticoids (GC) are the foundation of the chemotherapy regimen in acute lymphoblastic leukemia (ALL). However, resistance to GC is observed more frequently than resistance to other chemotherapy agents in patients with ALL relapse. Moreover, the mechanism underlying the development of GC resistance in ALL has not yet been fully uncovered. In this study, we used bioinformatic analysis methods to integrate the candidate genes and pathways participating in GC resistance in ALL and subsequently verified the bioinformatics findings with in vitro cell experiments. Ninety‐nine significant common differentially expressed genes (DEGs) associated with GC resistance were determined by integrating two gene profile datasets, including GC‐sensitive and ‐resistant samples. Using Kyoto Encyclopedia of Genes and Genomes (KEGG) and REACTOME pathways analysis, the signaling pathways in which DEGs were significantly enriched were clustered. The GC resistance‐related biologically functional interactions were visualized as DEG‐associated Protein–Protein Interaction (PPI) network complexes, with 98 nodes and 127 edges. MYC, a node which displayed the highest connectivity in all edges, was highlighted as the core gene in the PPI network. Increased C‐MYC expression was observed in adriamycin‐resistant BALL‐1/ADR cells, which we demonstrated was also resistant to dexamethasone. These results outlined a panorama in which the solitary and scattered experimental results were integrated and expanded. The potential promising target of the candidate pathways and genes involved in GC resistance of ALL was concomitantly revealed.