Cell Cycle Signaling Research Articles

Structure of nascent 5S RNPs at the crossroad between ribosome assembly and MDM2-p53 pathways.

The 5S ribonucleoprotein (RNP) is assembled from its three components (5S rRNA, Rpl5/uL18 and Rpl11/uL5) before being incorporated into the pre-60S subunit. However, when ribosome synthesis is disturbed, a free 5S RNP can enter the MDM2-p53 pathway to regulate cell cycle and apoptotic signaling. Here we reconstitute and determine the cryo-electron microscopy structure of the conserved hexameric 5S RNP with fungal or human factors. This reveals how the nascent 5S rRNA associates with the initial nuclear import complex Syo1-uL18-uL5 and, upon further recruitment of the nucleolar factors Rpf2 and Rrs1, develops into the 5S RNP precursor that can assemble into the pre-ribosome. In addition, we elucidate the structure of another 5S RNP intermediate, carrying the human ubiquitin ligase Mdm2, which unravels how this enzyme can be sequestered from its target substrate p53. Our data provide molecular insight into how the 5S RNP can mediate between ribosome biogenesis and cell proliferation.

Metastatic Infiltration of Nervous Tissue and Periosteal Nerve Sprouting in Multiple Myeloma-Induced Bone Pain in Mice and Human.

Multiple myeloma (MM) is a neoplasia of B plasma cells that often induces bone pain. However, the mechanisms underlying myeloma-induced bone pain (MIBP) are mostly unknown. Using a syngeneic MM mouse model, we show that periosteal nerve sprouting of calcitonin gene-related peptide (CGRP+) and growth associated protein 43 (GAP43+) fibers occurs concurrent to the onset of nociception and its blockade provides transient pain relief. MM patient samples also showed increased periosteal innervation. Mechanistically, we investigated MM induced gene expression changes in the dorsal root ganglia (DRG) innervating the MM-bearing bone of male mice and found alterations in pathways associated with cell cycle, immune response and neuronal signaling. The MM transcriptional signature was consistent with metastatic MM infiltration to the DRG, a never-before described feature of the disease that we further demonstrated histologically. In the DRG, MM cells caused loss of vascularization and neuronal injury, which may contribute to late-stage MIBP. Interestingly, the transcriptional signature of a MM patient was consistent with MM cell infiltration to the DRG. Overall, our results suggest that MM induces a plethora of peripheral nervous system alterations that may contribute to the failure of current analgesics and suggest neuroprotective drugs as appropriate strategies to treat early onset MIBP.SIGNIFICANCE STATEMENT Multiple myeloma (MM) is a painful bone marrow cancer that significantly impairs the quality of life of the patients. Analgesic therapies for myeloma-induced bone pain (MIBP) are limited and often ineffective, and the mechanisms of MIBP remain unknown. In this manuscript, we describe cancer-induced periosteal nerve sprouting in a mouse model of MIBP, where we also encounter metastasis to the dorsal root ganglia (DRG), a never-before described feature of the disease. Concomitant to myeloma infiltration, the lumbar DRGs presented blood vessel damage and transcriptional alterations, which may mediate MIBP. Explorative studies on human tissue support our preclinical findings. Understanding the mechanisms of MIBP is crucial to develop targeted analgesic with better efficacy and fewer side effects for this patient population.

Glycosyltransferase GLT8D1 and GLT8D2 serve as potential prognostic biomarkers correlated with Tumor Immunity in Gastric Cancer

BackgroundGlycosylation involved in various biological function, aberrant glycosylation plays an important role in cancer development and progression. Glycosyltransferase 8 domain containing 1 (GLT8D1) and GLT8D2, as members of the glycosyltransferase family proteins, are associated with transferase activity. However, the association between GLT8D1/2 and gastric cancer (GC) remains unclear. We aimed to investigate the potential prognostic value and oncogenic role of GLT8D1/2 in GC.MethodsThe relationship between GLT8D1/2 and GC was evaluated through comprehensive bioinformatics approaches. A series of factors like gene expression patterns, Kaplan-Meier survival analyses, Cox regression analyses, prognostic nomogram, calibration curves, ROC curves, function enrichment analyses, tumor immunity association, genetic alterations, and DNA methylation were included. Data and statistical analyses were performed using R software (v3.6.3).ResultsBoth GLT8D1 and GLT8D2 expression were significantly upregulated in GC tissues(n = 414) compared with normal tissues(n = 210), and high expression of GLT8D1/2 was remarkably correlated with poor prognosis for GC patients. Cox regression analyses implied that GLT8D1/2 could act as independent prognostic factors in GC. Furthermore, gene function analyses indicated that multiple signaling pathways involving tumor oncogenesis and development enriched, such as mTOR, cell cycle, MAPK, Notch, Hedgehog, FGF, and PI3K-Akt signaling pathways. Moreover, GLT8D1/2 was significantly associated with immune cell infiltration, immune checkpoint genes, and immune regulators TMB/MSI.ConclusionGLT8D1/2 may serve as potential prognostic markers of poor prognosis in GC correlated with tumor immunity. The study provided an insight into identifying potential biomarkers and targets for prognosis, immunotherapy response, and therapy in GC.

Effect of Borrelia burgdorferi on the Expression of miRNAs in Breast Cancer and Normal Mammary Epithelial Cells.

Breast cancer is one of the leading causes of death in women worldwide. Recent studies have demonstrated that inflammation due to infections with microorganisms could play a role in breast cancer development. One of the known human pathogens, Borrelia burgdorferi, the causative agent of Lyme disease, has been shown to be present in various types of breast cancer and is associated with poor prognosis. We reported that B. burgdorferi can invade breast cancer cells and affect their tumorigenic phenotype. To better understand the genome-wide genetic changes caused by B. burgdorferi, we evaluated the microRNA (miRNA or miR) expression profiles of two triple-negative breast cancer cell lines and one non-tumorigenic mammary cell line before and after B. burgdorferi infection. Using a cancer-specific miRNA panel, four miRNAs (miR-206, 214-3p, 16-5p, and 20b-5p) were identified as potential markers for Borrelia-induced changes, and the results were confirmed by quantitative real-time reverse transcription (qRT-PCR). Among those miRNAs, miR-206 and 214 were the most significantly upregulated miRNAs. The cellular impact of miR-206 and 214 was evaluated using DIANA software to identify related molecular pathways and genes. Analyses showed that the cell cycle, checkpoints, DNA damage-repair, proto-oncogenes, and cancer-related signaling pathways are mostly affected by B. burgdorferi infection. Based on this information, we have identified potential miRNAs which could be further evaluated as biomarkers for tumorigenesis caused by pathogens in breast cancer cells.

Effect of hypoxic status on survival and therapeutic responses in gastric cancer.

e16003 Background: As a significant factor of tumor microenvironment (TME), hypoxia is closely related to cancer metastasis, immune escape and drug resistance. The potential prognosis of hypoxia, as well as its influence on the TME, and therapeutic response have not yet been comprehensively studied in gastric cancer (GC). Here, we developed a prognosis risk assessment model (hypoxia_pred model) based on genes related to hypoxia, and further explored the exact role an impact of hypoxia in GC. Methods: Clinical information and gene expression profile were obtained from Gene Expression Omnibus (GEO) and the Cancer Genome Atlas (TCGA). 375 patients from TCGA-STAD cohort were used as training set. A 21-hypoxic genes expression dataset from TCGA-STAD was performed an unsupervised clustering analysis for different molecular subtypes of GC. Genes associated prognosis were screened by univariate Cox analysis with a cutoff of pvalue＜0.05. The prognostic risk model based on 9-genes transcriptome data of both hypoxia and prognosis related genes was constructed by Least Absolute Shrinkage and Selection Operator regression analysis. Another two gene expression datasets including 409 patients from GEO (GSE62254 and GSE26901) were employed for verifying the forecast performance of this model. In addition, we used ESTATE and single sample Gene Set Enrichment Analysis (ssGSEA) to explore tumor immune infiltration and potential molecular mechanism of hypoxia in GC. We also evaluated the influence of hypoxia on the prognosis, targeted therapy and immunotherapy of GC patients by survival analysis, targeted drug sensitivity analysis and the assessments of the Cancer Immunome Atlas database. Results: We found hypoxia_pred model is a valuable predictor for the survival of GC in validation cohort (AUC > 0.7). AUC for 1- and 3-year overall survival of GC patients in the TCGA cohort were 0.71 and 0.73. In GC, it was supposed that the aggravation of hypoxia was accompanied by the increase of immune rejection phenotype and the infiltration of immunosuppressive cells. In the validation set, patients predicted high risk by the model have worse immunotherapeutic response and prognosis, and may profit from drugs that inhibit cell cycle signaling pathways. In addition, the outcomes by hypoxia_pred for changes in TME, efficacy of immunotherapy and drug sensitivity were quite hypoxia_pred. targeted drug sensitivity, immune microenvironment, biological function, immunotherapeutic response, and immune cell infiltration were found to be significantly different in those molecular subtypes and risk groups. Conclusions: The model based on hypoxic environment provides theoretical foundation for predicting tumor prognosis and clinical treatment, and may improve the prognosis of GC patients by recommending personalized immunotherapy.

Pathway enrichment analysis in tumors with high mutation burdens and interferon-γ signature expression: A pancancer analysis.

3138 Background: Tumor mutation burden (TMB) correlates with immunotherapy response, but outcomes still vary in high TMB cancers, and many high TMB tumors lack T cell infiltration. Here, we assessed gene expression signatures of high TMB, including both inflamed and non-inflamed tumors (defined by interferon-γ [IFN-γ] gene signatures). Methods: Gene set enrichment analysis (GSEA) was assessed from RNA-sequencing data in tumors in TCGA. TMB and IFN-γ signature score were assessed by previously published methods. We performed GSEA on high vs low TMB tumors, and further stratified by high vs low IFN-γ signature expression within high TMB tumors. We validated results in Tempus RNA-sequencing data from 264 patients at Vanderbilt. Results: 4020 TCGA cases in non-small cell lung, melanoma, bladder, head and neck, kidney, and colorectal cancer were included. Major signatures upregulated in the high TMB group included cell cycle signaling, nuclear protein/mRNA export, cholesterol biosynthetics, and ubiquitin signaling. High TMB/high IFN-γ signature tumors were enriched in immune-related pathways (neutrophil and monocyte chemotaxis, T cell activation and proliferation, and interleukin-1 secretion), and non-immune pathways including cAMP metabolic process and phosphatidylinositol signaling. High TMB/low IFN-γ signature tumors had enriched pathways in nonsense mediated decay, electron transport chain/oxidative phosphorylation, and mitochondrial translational disassembly. 264 tumors from Vanderbilt were used for validation; pathways including cholesterol biosynthetics also correlated with immunotherapy response. Conclusions: High TMB tumors were enriched in pathways such as ubiquitin signaling, which could regulate immune cell function in antitumor immunity. Within high TMB cancers, tumors with higher IFN-γ signatures have enriched immune-related, and metabolic pathways, which could suggest novel therapeutic targets to make high TMB/ low IFN-γ tumors more immunogenic.

Identification of small-molecule inhibitors against the interaction of RNA-binding protein PSF and its target RNA for cancer treatment.

Diverse cellular activities are modulated through a variety of RNAs, including long noncoding RNAs (lncRNAs), by binding to certain proteins. The inhibition of oncogenic proteins or RNAs is expected to suppress cancer cell proliferation. We have previously demonstrated that PSF interaction with its target RNAs, such as androgen-induced lncRNA CTBP1-AS, is critical for hormone therapy resistance in prostate and breast cancers. However, the action of protein-RNA interactions remains almost undruggable to date. High-throughput screening (HTS) has facilitated the discovery of drugs for protein-protein interactions. In the present study, we developed an in vitro alpha assay using Flag peptide-conjugated lncRNA, CTBP1-AS, and PSF. We then constructed an effective HTS screening system to explore small compounds that inhibit PSF-RNA interactions. Thirty-six compounds were identified and dose-dependently inhibited PSF-RNA interaction in vitro. Moreover, chemical optimization of these lead compounds and evaluation of cancer cell proliferation revealed two promising compounds, N-3 and C-65. These compounds induced apoptosis and inhibited cell growth in prostate and breast cancer cells. By inhibiting PSF-RNA interaction, N-3 and C-65 up-regulated signals that are repressed by PSF, such as the cell cycle signals by p53 and p27. Furthermore, using a mouse xenograft model for hormone therapy-resistant prostate cancer, we revealed that N-3 and C-65 can significantly suppress tumor growth and downstream target gene expression, such as the androgen receptor (AR). Thus, our findings highlight a therapeutic strategy through the development of inhibitors for RNA-binding events in advanced cancers.

Key genes involved with prognosis were identified in lung adenocarcinoma by integrated bioinformatics analysis

ObjectiveBy screening the core genes in lung adenocarcinoma (LUAD) with bioinformatics, our study evaluated its prognosis value and role in infiltration process of immune cells. MethodsUsing GEO database, we screened 5 gene chips, including GSE11072, GSE32863, GSE43458, GSE115002, and GSE116959. Then, we obtained the corresponding differentially expressed genes by analyzed 5 gene chips online by GEO2R (P < 0.05, |logFC| > 1). Then, through DAVID online platform, Cytoscape 3.6.1 software and PPI network analysis, the network was visualized and obtain the final core genes. Next, we plan to use the GEPIA, UALCAN, Kaplan–Meier plotter and Time 2.0 database for corresponding analysis. The GEPIA database was used to verify the expression of core genes in LUAD and normal lung tissues, and survival analysis was used to evaluate the value of core genes in the prognosis of LUAD patients. UALCAN was used to verify the expression of the LUAD core gene and promoter methylation status, and the predictive value of core genes was evaluated in LUAD patients by the Kaplan–Meier plotter online tool. Then, we used the Time 2.0 database to identify the relationship to immune infiltration in LUAD. Finally, we used the human protein atlas (HPA) database for online immunohistochemical analysis of the expressed proteins. ResultsThe expression of CCNB2 and CDC20 in LUAD were higher than those in normal lung tissues, their increased expression was negatively correlated with the overall survival rate of LUAD, and they were involved in cell cycle signal transduction, oocyte meiosis signal transduction as well as the infiltration process of immune cells in LUAD. The expression proteins of CCNB2 and CDC20 were also different in lung cancer tissue and normal lung tissue. Therefore, CCNB2 and CDC20 were identified as the vital core genes. ConclusionCCNB2 and CDC20 are essential genes that may constitute prognostic biomarkers in LUAD, they also participate the immune infiltration process and protein expression process of LUAD, and might provides basis for clinical anti-tumor drug research.

Identification of Key Genes Related to Immune Cells in Patients with COVID-19 Via Integrated Bioinformatics-Based Analysis.

COVID-19 has spread all over the world which poses a serious threat to social economic development and public health. Despite enormous progress has been made in the prevention and treatment of COVID-19, the specific mechanism and biomarker related to disease severity or prognosis have not been clarified yet. Our study intended to further explore the diagnostic markers of COVID-19 and their relationship with serum immunology by bioinformatics analysis. The datasets about COVID-19 were downloaded from the Gene Expression Omnibus (GEO) dataset. The differentially expressed genes (DEGs) were selected via the limma package. Then, weighted gene co-expression network analysis (WGCNA) was conducted to identify the critical module associated with the clinic status. The intersection DEGs were processed for further enrichment analysis. The final diagnostic genes for COVID-19 were selected and verified through special bioinformatics algorithms. There were significant DEGs between the normal and COVID-19 patients. These genes were mainly enriched in cell cycle, complement and coagulation cascade, extracellular matrix (ECM) receptor interaction, and the P53 signaling pathway. As much as 357 common intersected DEGs were selected in the end. These DEGs were enriched in organelle fission, mitotic cell cycle phase transition, DNA helicase activity, cell cycle, cellular senescence, and P53 signaling pathway. Our study also identified CDC25A, PDCD6, and YWAHE were potential diagnostic markers of COVID-19 with the AUC (area under curve), 0.958 (95% CI 0.920-0.988), 0.941(95% CI 0.892-0.980), and 0.929 (95% CI 0.880-0.971). Moreover, CDC25A, PDCD6, and YWAHE were correlated with plasma cells, macrophages M0, T cells CD4 memory resting, T cells CD8, dendritic cells, and NK cells. Our study discovered that CDC25A, PDCD6, and YWAHE can be used as diagnostic markers for COVID-19. Moreover, these biomarkers were also closely associated with immune cell infiltration, which plays a pivotal role in the diagnosis and progression of COVID-19.

Bioinformatics analysis of aging-related genes in thoracic aortic aneurysm and dissection.

Thoracic aortic aneurysm and dissection (TAAD) is a cardiovascular disease with a high mortality rate. Aging is an important risk factor for TAAD. This study explored the relationship between aging and TAAD and investigated the underlying mechanisms, which may contribute to the diagnosis and treatment of TAAD. Human aging genes were obtained from the Aging Atlas official website. Various datasets were downloaded from the GEO database:the human TAAD dataset GSE52093 were used for screening differentially expressed genes (DEGs); GSE137869, GSE102397 and GSE153434 were used as validation sets, and GSE9106 was used for diagnostic prediction of receiver operating characteristic (ROC) curves. Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), Gene Set Enrichment Analysis (GSEA), and protein-protein interaction (PPI) network analysis were used to screen differentially co-expressed genes from human aging genes and TAAD. Using five methods of the cytoHubba plugin in Cytoscape software (Degree, Closeness, EPC, MNC, Radiality), hub genes were identified from the differentially co-expressed genes. Single-cell RNA sequencing was used to verify the expression levels of hubgenes in different cell types of aortic tissue. ROC curves were used to further screen for diagnostic genes. A total of 70 differentially co-expressed genes were screened from human aging genes and DEGs in human TAAD dataset GSE52093. GO enrichment analysis revealed that the DEGs played a major role in regulating DNA metabolism and damaged DNA binding. KEGG enrichment analysis revealed enrichment in the longevity regulating pathway, cellular senescence, and HIF-1 signaling pathway. GSEA indicated that the DEGs were concentrated in the cell cycle and aging-related p53 signaling pathway. The five identified hubgenes were MYC, IL6, HIF1A, ESR1, and PTGS2. Single-cell sequencing of the aging rat aorta showed that hubgenes were expressed differently in different types of cells in aortic tissue. Among these five hubgenes, HIF1A and PTGS2 were validated in the aging dataset GSE102397; MYC, HIF1A and ESR1 were validated in the TAAD dataset GSE153434. The combined area under the diagnostic ROC curve (AUC) values for the five hub genes were >0.7 in the testing and training sets of the dataset GSE9106. The combined AUC values of MYC and ESR1 were equal to the combin ed AUC values of the five hub genes. The HIF-1 signaling pathway may play an important role in TAAD and aging. MYC and ESR1 may have diagnostic value for aging-related TAAD.

The Prognostic Value and Potential Mechanism of Tumor-Nutrition-Inflammation Index and Genes in Patients with Advanced Lung Cancer.

Lung cancer (LC) has the highest mortality rate all over the world. It is necessary to search for novel potential biomarkers that are easily accessible and inexpensive in identifying patients with LC at early stage. A total of 195 patients with advanced LC who have received first-line chemotherapy were involved in this study. The optimized cut-off values of AGR and SIRI (AGR = albumin/globulin; SIRI = neutrophil ∗ monocyte/lymphocyte) were determined by survival function analysis based on R software. COX regression analysis was performed to obtain the independent factors for establishing the nomogram model. A nomogram model comprising these independent prognostic parameters was built for the TNI (tumor-nutrition-inflammation index) score calculation. The predictive accuracy was demonstrated through ROC curve and calibration curves after index concordance. The optimized cut-off values of AGR and SIRI were 1.22 and 1.60, respectively. It was revealed that liver metastasis, SCC, AGR, and SIRI were independent prognostic factors in advanced lung cancer by Cox analysis. Afterwards, the nomogram model comprised of these independent prognostic parameters was built for TNI scores calculation. Based on the TNI quartile values, patients were divided into four groups. And it was indicated that higher TNI had worse OS (P < 0.05) via Kaplan-Meier analysis and log-rank test. Moreover, the C-index and 1-year AUC area were 0.756 (0.723-0.788) and 75.62, respectively. There was high consistency shown in the calibration curves between predicted and actual survival proportions in the TNI model. In addition, tumor-nutrition-inflammation index and genes play an important role in LC development that might affect some pathways related to tumor development including cell cycle, homologous recombination, and P53 signaling pathway from a molecular level. TNI might be an analytical tool which was practical and precise for survival prediction of patients with advanced LC. Tumor-nutrition-inflammation index and genes play an important role in LC development. A preprint has previously been published [1].

Bioinformatics-Based Identification of CircRNA-MicroRNA-mRNA Network for Calcific Aortic Valve Disease.

Calcific aortic valve disease (CAVD) is the most common native valve disease. Valvular interstitial cell (VIC) osteogenic differentiation and valvular endothelial cell (VEC) dysfunction are key steps in CAVD progression. Circular RNA (circRNAs) is involved in regulating osteogenic differentiation with mesenchymal cells and is associated with multiple disease progression, but the function of circRNAs in CAVD remains unknown. Here, we aimed to investigate the effect and potential significance of circRNA-miRNA-mRNA networks in CAVD. Two mRNA datasets, one miRNA dataset, and one circRNA dataset of CAVD downloaded from GEO were used to identify DE-circRNAs, DE-miRNAs, and DE-mRNAs. Based on the online website prediction function, the common mRNAs (FmRNAs) for constructing circRNA-miRNA-mRNA networks were identified. GO and KEGG enrichment analyses were performed on FmRNAs. In addition, hub genes were identified by PPI networks. Based on the expression of each data set, the circRNA-miRNA-hub gene network was constructed by Cytoscape (version 3.6.1). 32 DE-circRNAs, 206 DE-miRNAs, and 2170 DE-mRNAs were identified. Fifty-nine FmRNAs were obtained by intersection. The KEGG pathway analysis of FmRNAs was enriched in pathways in cancer, JAK-STAT signaling pathway, cell cycle, and MAPK signaling pathway. Meanwhile, transcription, nucleolus, and protein homodimerization activity were significantly enriched in GO analysis. Eight hub genes were identified based on the PPI network. Three possible regulatory networks in CAVD disease were obtained based on the biological functions of circRNAs including: hsa_circ_0026817-hsa-miR-211-5p-CACNA1C, hsa_circ_0007215-hsa-miR-1252-5p-MECP2, and hsa_circ_0007215-hsa-miR-1343-3p- RBL1. The present bionformatics analysis suggests the functional effect for the circRNA-miRNA-mRNA network in CAVD pathogenesis and provides new targets for therapeutics.

The SMC5/6 complex subunit MMS21 regulates stem cell proliferation in rice.

SMC5/6 complex subunit OsMMS21 is involved in cell cycle and hormone signaling and required for stem cell proliferation during shoot and root development in rice. The structural maintenance of chromosome (SMC)5/6 complex is required for nucleolar integrity and DNA metabolism. Moreover, METHYL METHANESULFONATE SENSITIVITY GENE 21 (MMS21), a SUMO E3 ligase that is part of the SMC5/6 complex, is essential for the root stem cell niche and cell cycle transition in Arabidopsis. However, its specific role in rice remains unclear. Here, OsSMC5 and OsSMC6 single heterozygous mutants were generated using CRISPR/Cas9 technology to elucidate the function of SMC5/6 subunits, including OsSMC5, OsSMC6, and OsMMS21, in cell proliferation in rice. ossmc5/ + and ossmc6/ + heterozygous single mutants did not yield homozygous mutants in their progeny, indicating that OsSMC5 and OsSMC6 both play necessary roles during embryo formation. Loss of OsMMS21 caused severe defects in both the shoot and roots in rice. Transcriptome analysis showed a significant decrease in the expression of genes involved in auxin signaling in the roots of osmms21 mutants. Moreover, the expression levels of the cycB2-1 and MCM genes, which are involved the cell cycle, were significantly lower in the shoots of the mutants, indicating that OsMMS21 was involved in both hormone signaling pathways and the cell cycle. Overall, these findings indicate that the SUMO E3 ligase OsMMS21 is required for both shoot and root stem cell niches, improving the understanding of the function of the SMC5/6 complex in rice.

Enteric nervous system damage caused by abnormal intestinal butyrate metabolism may lead to functional constipation.

Functional constipation (FC) is a high morbidity gastrointestinal disease for which dysfunction in the enteric nervous system is a major pathogenesis mechanism. To enhance our understanding of the involvement of intestinal microbiota and its metabolites in the pathogenesis of FC, we conducted a shotgun metagenomic sequencing analysis of gut microbiota and serum short-chain fatty acids (SCFAs) analysis in 460 Chinese women with different defecation frequencies. We observed that the abundance ofFusobacterium_varium, a butyric acid-producing bacterium, was positively correlated (P = 0.0096) with the frequency of defecation; however, the concentrations of serum butyric acid was negatively correlated (P = 3.51E-05) with defecation frequency. These results were verified in an independent cohort (6 patients with FC and 6 controls). To further study the effects of butyric acid on intestinal nerve cells, we treated mouse intestinal neurons in vitro with various concentrations of butyrate (0.1, 0.5, 1, and 2.5 mM). We found that intestinal neurons treated with 0.5 mM butyrate proliferated better than those in the other treatment groups, with significant differences in cell cycle and oxidative phosphorylation signal pathways. We suggest that the decreased butyrate production resulting from the reduced abundance of Fusobacterium in gut microbiota affects the proliferation of intestinal neurons and the energy supply of intestinal cells. However, with FC disease advancing, the consumption and excretion of butyric acid reduce, leading to its accumulation in the intestine. Moreover, the accumulation of an excessively high amount of butyric acid inhibits the proliferation of nerve cells and subsequently exacerbates the disease.

A prognostic model based on the Augmin family genes for LGG patients

Gliomas are the most prevalent primary tumors in the central nervous system. Despite some breakthroughs in the treatment of glioma in recent years, survival rates remain low. Although genes of the Augmin family play a key role in microtubule nucleation, the role they play in gliomas is unclear. Transcriptome data were extracted from UCSC XENA and GTEx for low-grade glioma (LGG) and normal tissues, respectively. The protein interaction network associated with Augmin family genes was established using STRING and GeneMANIA databases. Enrichment analysis of gene-related functions and pathways was used to explore potential biological pathways and TIMER to assess immune cell infiltration. Regression analysis and Kaplan–Meier analysis were used to look at the clinical characteristics of the Augmin family genes and the association with the prognosis of patients with glioma. The results showed that the mRNA expression of Augmin family genes was significantly elevated in LGG tissues, except for HAUS7. Immunoregulation, cell cycle, apoptosis and other signaling pathways may be involved in the development and progression of LGG. Except for HAUS4 and HAUS7, the expression of all genes was positively correlated with immune cell infiltration. High expression of HAUS1, HAUS3, HAUS5, HAUS7, HAUS8 and low expression of HAUS4, HAUS6 in LGG was associated with poor prognosis. The risk models constructed based on the pivotal genes HAUS2, HAUS4 and HAUS8 were validated by nomogram and confirmed to be clinically useful for predicting the prognosis of LGG.

Zinc finger and SCAN domain-containing protein 18 is a potential DNA methylation-modified tumor suppressor and biomarker in breast cancer.

Zinc finger and SCAN domain-containing protein 18 (ZSCAN18) has been investigated as a putative biomarker of multiple human cancers. However, the expression profile, epigenetic modification, prognostic value, transcription regulation, and molecular mechanism of ZSCAN18 in breast cancer (BC) remain unknown. In the study, we present an integrated analysis of ZSCAN18 in BC based on public omics datasets with the use of multiple bioinformatics tools. Genes potentially regulated through restoration of ZSCAN18 expression in MDA-MB-231 cells were investigated to identify pathways associated with BC. We observed that ZSCAN18 was downregulated in BC and mRNA expression was significantly correlated with clinicopathological parameters. Low expression of ZSCAN18 was found in the HER2-positive and TNBC subtypes. High expression of ZSCAN18 was associated with good prognosis. As compared to normal tissues, the extent of ZSCAN18 DNA methylation was greater with fewer genetic alterations in BC tissues. ZSCAN18 was identified as a transcription factor that might be involved in intracellular molecular and metabolic processes. Low ZSCAN18 expression was associated with the cell cycle and glycolysis signaling pathway. Overexpression of ZSCAN18 inhibited mRNA expression of genes associated with the Wnt/β-catenin and glycolysis signaling pathways, including CTNNB1, BCL9, TSC1, and PFKP. ZSCAN18 expression was negatively correlated with infiltrating B cells and dendritic cells (DCs), as determined by the TIMER web server and reference to the TISIDB. ZSCAN18 DNA methylation was positively correlated with activated B cells, activated CD8+ and CD4+ T cells, macrophages, neutrophils, and activated DCs. Moreover, five ZSCAN18-related hub genes (KDM6B, KAT6A, KMT2D, KDM1A, and HSPBP1) were identified. ZSCAN18, ZNF396, and PGBD1 were identified as components of a physical complex. ZSCAN18 is a potential tumor suppressor in BC, as expression is modified by DNA methylation and associated with patient survival. In addition, ZSCAN18 plays important roles in transcription regulation, the glycolysis signaling pathway, and the tumor immune microenvironment.

Clinical significance and biology of circulating tumor DNA in high-risk early-stage HER2-negative breast cancer receiving neoadjuvant chemotherapy

Circulating tumor DNA (ctDNA) analysis may improve early-stage breast cancer treatment via non-invasive tumor burden assessment. To investigate subtype-specific differences in the clinical significance and biology of ctDNA shedding, we perform serial personalized ctDNA analysis in hormone receptor (HR)-positive/HER2-negative breast cancer and triple-negative breast cancer (TNBC) patients receiving neoadjuvant chemotherapy (NAC) in the I-SPY2 trial. ctDNA positivity rates before, during, and after NAC are higher in TNBC than in HR-positive/HER2-negative breast cancer patients. Early clearance of ctDNA 3weeks after treatment initiation predicts a favorable response to NAC in TNBC only. Whereas ctDNA positivity associates with reduced distant recurrence-free survival in both subtypes. Conversely, ctDNA negativity after NAC correlates with improved outcomes, even in patients with extensive residual cancer. Pretreatment tumor mRNA profiling reveals associations between ctDNA shedding and cell cycle and immune-associated signaling. On the basis of these findings, the I-SPY2 trial will prospectively test ctDNA for utility in redirecting therapy to improve response and prognosis.

Genomic and transcriptomic features between primary and paired metastatic fumarate hydratase–deficient renal cell carcinoma

BackgroundFumarate hydratase–deficient renal cell carcinoma (FH-RCC) is a rare highly aggressive subtype of kidney cancer for which the distinct genomic, transcriptomic, and evolutionary relationships between metastatic and primary lesions are still unclear.MethodsIn this study, whole-exome, RNA-seq, and DNA methylation sequencing were performed on primary-metastatic paired specimens from 19 FH-RCC cases, including 23 primary and 35 matched metastatic lesions. Phylogenetic and clonal evolutionary analyses were used to investigate the evolutionary characteristics of FH-RCC. Transcriptomic analyses, immunohistochemistry, and multiple immunofluorescence experiments were performed to identify the tumor microenvironmental features of metastatic lesions.ResultsPaired primary and metastatic lesions generally showed similar characteristics of tumor mutation burden, tumor neoantigen burden, microsatellite instability score, CNV burden, and genome instability index. Notably, we identified an FH-mutated founding MRCA (the most recent common ancestor) clone that dominated the early evolutionary trajectories in FH-RCC. Although both primary and metastatic lesions manifested high immunogenicity, metastatic lesions exhibited higher enrichment of T effector cells and immune-related chemokines, together with upregulation of PD-L1, TIGIT, and BTLA. In addition, we found that concurrent NF2 mutation may be associated with bone metastasis and upregulation of cell cycle signature in metastatic lesions. Furthermore, although in FH-RCC metastatic lesions in general shared similar CpG island methylator phenotype with primary lesions, we found metastatic lesions displaying hypomethylated chemokine and immune checkpoints related genomic loci.ConclusionsOverall, our study demonstrated the genomic, epigenomic, and transcriptomic features of metastatic lesions in FH-RCC and revealed their early evolutionary trajectory. These results provided multi-omics evidence portraying the progression of FH-RCC.

Abstract B012: Effects of adagrasib on oncogenic signaling, immune cell regulation and biomarkers of response in preliminary clinical analyses

Abstract The recent approval of KRAS G12C inhibitors has provided novel treatment options for lung, colorectal and other cancer patients harboring this mutation and has changed the landscape of research for these cancers. Current priorities in the field include characterizing the mechanism of action of KRAS G12C inhibitors and the mechanisms underlying resistance to KRAS inhibition. Utilizing access to pre and post-treatment patient samples from KRYSTAL-1, preliminary gene expression analyses from baseline and post-adagrasib (Cycle 1, Day 8) treated patient samples demonstrated multiple oncogenic pathway gene signatures, including MYC, mTOR, cell cycle, EMT and inflammation-related signatures, were significantly regulated by adagrasib. KRAS/MAPK pathway and cell cycle-related genes including ETV4, CCND1, DUSP6, TOP1, and CENPA, as well as the Singh KRAS dependency signature were significantly downregulated in Cycle1, Day 8 compared to baseline tumor biopsies. Genes implicated in inflammation and the immune response were significantly upregulated in Cycle 1, Day 8 compared to baseline tumor biopsies. In agreement with the gene expression data as well as findings from syngeneic mouse models, preliminary immunohistochemical analysis on matched baseline, post-adagrasib treated (Cycle 1, Day 8), and end of treatment FFPE patient samples revealed marked changes in tumor cell mechanistic biomarkers and immune cell types following adagrasib treatment. These observations included decreased tumor Ki67 and increased tumor PD-L1, as well as marked alterations in tumor immune cell composition, including increased CD8+ T cells and decreased MDSCs following adagrasib treatment in several patients. Finally, flow cytometry and TCRb sequencing data from adagrasib and pembrolizumab-treated patient blood samples demonstrated an increase in several activated CD8+ T cell populations and the emergence of new T cell clones in a subset of patients after combination treatment. Together, these data suggest robust evaluation of patient biopsies pre and post adagrasib treatment may better predict outcomes that can be achieved in patients harboring KRASG12C mutant cancers and can guide early clinical development strategies, including support for combining adagrasib with immune checkpoint inhibitors. Citation Format: Jill Hallin, Laura Hover, Julio Fernandez-Benet, Adam Pavlicek, Kenna Anderes, Pasi A. Jänne, Gregory J. Riely, Alexander I. Spira, Jun Zhang, Peter Olson, James G. Christensen. Effects of adagrasib on oncogenic signaling, immune cell regulation and biomarkers of response in preliminary clinical analyses [abstract]. In: Proceedings of the AACR Special Conference: Targeting RAS; 2023 Mar 5-8; Philadelphia, PA. Philadelphia (PA): AACR; Mol Cancer Res 2023;21(5_Suppl):Abstract nr B012.

The Prognostic and Drug-targeting Value of Lymphoid Enhancer-binding Factor-1 in Hepatocellular Carcinoma.

Lymphoid Enhancer-Binding Factor-1 (LEF1) was previously reported to contribute to a variety of malignancies, including Hepatocellular Carcinoma (HCC). However, its role in HCC is poorly understood. To explore the role of LEF1 in HCC, including its prognostic and drug-targeting value. The LEF1 expression and patient characteristics were investigated. The associations between clinical characteristics and LEF1 were analyzed using both univariate and multivariate logistic regression. Cox regression and Kaplan-Meier curves were used to explore the clinicopathological factors related to overall survival in patients with HCC. A nomogram to predict the survival rate was constructed and validated. The Kyoto Encyclopedia of Genes and Genomes database (KEGG) was used to explore the function of LEF1. Gene Set Enrichment Analysis (GSEA) was also performed using The Cancer Genome Atlas dataset. Furthermore, compounds that may have the potential to be targeted drugs in the treatment of LEF1-overexpressing HCC were identified using the Comparative Toxicogenomics Database (CTD), patents about these drugs in HCC were also reviewed through Worldwide Espacenet® and Patentscope®. Increased expression of LEF1 was significantly associated with high histological grade of HCC (odds ratio (OR) = 2.521 for grade (G) 2 vs. G1, OR = 2.550 for G3 vs. G1, OR = 7.081 for G4 vs. G1, all P < 0.05). A Kaplan-Meier survival curve showed that HCC patients with LEF1 overexpression had a poor prognosis compared with those with normal LEF1 expression (P = 0.025). Multivariate Cox regression analysis revealed that LEF1 is an independent prognostic factor for the overall survival of patients with HCC (Hazard Ratio (HR) = 1.095; P = 0.04). The constructed nomogram to predict the survival rate produced a statistically significant prediction (area under the curve (AUC) = 86.68). In addition, Gene Ontology (GO) and KEGG analysis of genes co-expressed with the protein showed that LEF1 was associated with transcriptional regulation. GSEA suggested that the cell cycle, the WNT signaling pathway, and the NOTCH signaling pathway may be the key pathways regulated by LEF1 in HCC. Furthermore, the Comparative Toxicogenomics Database (CTD) identified nine compounds that may have the potential to be targeted drugs in the treatment of LEF1-overexpressing HCC. Patent reviews suggested that these drugs may show some efficacy in HCC, but whether these drugs interact with LEF1 and improve the prognosis for patients with HCC remains to be explored. LEF1 is a latent prognostic molecular biomarker of HCC. The cell cycle, and WNT and NOTCH signaling pathways are regulated by LEF1 in HCC. LEF1 could be a potential drug target for HCC.