Analysis Of Single-cell RNA Sequencing Data Research Articles

Improving cell type identification with Gaussian noise-augmented single-cell RNA-seq contrastive learning.

Cell type identification is an important task for single-cell RNA-sequencing (scRNA-seq) data analysis. Many prediction methods have recently been proposed, but the predictive accuracy of difficult cell type identification tasks is still low. In this work, we proposed a novel Gaussian noise augmentation-based scRNA-seq contrastive learning method (GsRCL) to learn a type of discriminative feature representations for cell type identification tasks. A large-scale computational evaluation suggests that GsRCL successfully outperformed other state-of-the-art predictive methods on difficult cell type identification tasks, while the conventional random genes masking augmentation-based contrastive learning method also improved the accuracy of easy cell type identification tasks in general.

Cytoplasmatic Localization of Six1 in Male Testis and Spermatogonial Stem Cells.

Sine oculis homeobox 1 (Six1) is an important factor for embryonic development and carcinoma malignancy. However, the localization of Six1 varies due to protein size and cell types in different organs. In this study, we focus on the expression and localization of Six1 in male reproductive organ via bioinformatics analysis and immunofluorescent detection. The potential interacted proteins with Six1 were also predicted by protein-protein interactions (PPIs) and Enrichr analysis. Bioinformatic data from The Cancer Genome Atlas and Genotype-Tissue Expression project databases showed that SIX1 was highly expressed in normal human testis, but low expressed in the testicular germ cell tumor sample. Human Protein Atlas examination verified that SIX1 level was higher in normal than that in cancer samples. The sub-localization of SIX1 in different reproductive tissues varies but specifically in the cytoplasm and membrane in testicular cells. In mouse cells, single cell RNA-sequencing data analysis indicated that Six1 expression level was higher in mouse spermatogonial stem cells (mSSCs) and differentiating spermatogonial than in other somatic cells. Immunofluorescence staining showed the cytoplasmic localization of Six1 in mouse testis and mSSCs. Further PPIs and Enrichr examination showed the potential interaction of Six1 with bone morphogenetic protein 4 (Bmp4) and catenin Beta-1 (CtnnB1) and stem cell signal pathways. Cytoplasmic localization of Six1 in male testis and mSSCs was probably associated with stem cell related proteins Bmp4 and CtnnB1 for stem cell development.

Integrative analysis of single-cell and bulk RNA-sequencing data revealed T cell marker genes based molecular sub-types and a prognostic signature in lung adenocarcinoma

Immunotherapy has emerged as a promising modality for addressing advanced or conventionally drug-resistant malignancies. When it comes to lung adenocarcinoma (LUAD), T cells have demonstrated significant influence on both antitumor activity and the tumor microenvironment. However, their specific contributions remain largely unexplored. This investigation aimed to delineate molecular subtypes and prognostic indicators founded on T cell marker genes, thereby shedding light on the significance of T cells in LUAD prognosis and precision treatment. The cellular phenotypes were identified by scrutinizing the single-cell data obtained from the GEO repository. Subsequently, T cell marker genes derived from single-cell sequencing analyses were integrated with differentially expressed genes from the TCGA repository to pinpoint T cell-associated genes. Utilizing Cox analysis, molecular subtypes and prognostic signatures were established and subsequently verified using the GEO dataset. The ensuing molecular and immunological distinctions, along with therapy sensitivity between the two sub-cohorts, were examined via the ESTIMATE, CIBERSORT, and ssGSEA methodologies. Compartmentalization, somatic mutation, nomogram development, chemotherapy sensitivity prediction, and potential drug prediction analyses were also conducted according to the risk signature. Additionally, real-time qPCR and the HPA database corroborated the mRNA and protein expression patterns of signature genes in LUAD tissues. In summary, this research yielded an innovative T cell marker gene-based signature with remarkable potential to prognosis and anticipate immunotherapeutic outcomes in LUAD patients.

Predicting mitophagy-related genes and unveiling liver endothelial cell heterogeneity in hepatic ischemia-reperfusion injury.

Hepatic Ischemia-Reperfusion Injury (HIRI) is a major complication in liver transplants and surgeries, significantly affecting postoperative outcomes. The role of mitophagy, essential for removing dysfunctional mitochondria and maintaining cellular balance, remains unclear in HIRI. To unravel the role of mitophagy-related genes (MRGs) in HIRI, we assembled a comprehensive dataset comprising 44 HIRI samples alongside 44 normal control samples from the Gene Expression Omnibus (GEO) database for this analysis. Using Random Forests and Support Vector Machines - Recursive Feature Elimination (SVM-RFE), we pinpointed eight pivotal genes and developed a logistic regression model based on these findings. Further, we employed consensus cluster analysis for classifying HIRI patients according to their MRG expression profiles and conducted weighted gene co-expression network analysis (WGCNA) to identify clusters of genes that exhibit high correlation within different modules. Additionally, we conducted single-cell RNA sequencing data analysis to explore insights into the behavior of MRGs within the HIRI. We identified eight key genes (FUNDC1, VDAC1, MFN2, PINK1, CSNK2A2, ULK1, UBC, MAP1LC3B) with distinct expressions between HIRI and controls, confirmed by PCR validation. Our diagnostic model, based on these genes, accurately predicted HIRI outcomes. Analysis revealed a strong positive correlation of these genes with monocytic lineage and a negative correlation with B and T cells. HIRI patients were divided into three subclusters based on MRG profiles, with WGCNA uncovering highly correlated gene modules. Single-cell analysis identified two types of endothelial cells with different MRG scores, indicating their varied roles in HIRI. Our study highlights the critical role of MRGs in HIRI and the heterogeneity of endothelial cells. We identified the macrophage migration inhibitory factor (MIF) and cGAS-STING (GAS) pathways as regulators of mitophagy's impact on HIRI. These findings advance our understanding of mitophagy in HIRI and set the stage for future research and therapeutic developments.

A High Immune-Related Index with the Suppression of cGAS-STING Pathway is a Key Determinant to Herceptin Resistance in HER2+ Breast Cancer.

Resistance to HER2-targeted therapy is the major cause of treatment failure in patients with HER2+ breast cancer (BC). Given the key role of immune microenvironment in tumor development, there is a lack of an ideal prognostic model that fully accounts for immune infiltration. In this study, WGCNA analysis was performed to discover the relationship between immune-related signaling and prognosis of HER2+ BC. After Herceptin-resistant BC cell lines established, transcriptional profiles of resistant cell line and RNA-sequencing data from GSE76360 cohort were analyzed for candidate genes. 85 samples of HER2+ BC from TCGA database were analyzed by the Cox regression, XGBoost and Lasso algorithm to generalize a credible immune-related prognostic index (IRPI). Correlations between the IRPI signature and tumor microenvironment were further analyzed by multiple algorithms, including single-cell RNA sequencing data analysis. Patients with high IRPI had suppressive tumor immune microenvironment and worse prognosis. The suppression of type I interferon signaling indicated by the IRPI in Herceptin-resistant HER2+ BC was validated. And we elucidated that the suppression of cGAS-STING pathway is the key determinant underlying immune escape in Herceptin-resistant BC with high IRPI. A combination of STING agonist and DS-8201 could serve as a new strategy for Herceptin-resistant HER2+ BC.

Artesunate attenuates serum amyloid A-induced M1 macrophage differentiation through the promotion of PHGDH

Clinical studies indicated that Serum Amyloid A (SAA) might be a promising biomarker for forecasting the activity, severity, and adverse prognosis of systemic lupus erythematosus (SLE). Simultaneously, a positive correlation has been observed between macrophages, Th17 cells, and SLE disease activity, with both these immune cells being affected by SAA. Presently, the relationship between SAA and the aforementioned immune cell types in SLE remains to be elucidated. To discern the immune cell type most closely associated with SAA, we undertook a single-cell RNA sequencing data analysis via the GEO database. Subsequent results revealed a strong association between macrophages and SAA, a relationship further validated through flow cytometry of spleen macrophages in the MRL/lpr model. We discovered that SAA stimulate M1 macrophage differentiation along with the upregulation of pro-inflammatory cytokines such as IL-6 and IL-1β. Our findings suggest that SAA may promote M1 macrophage differentiation via the downregulation of phosphoglycerate dehydrogenase (PHGDH). Artesunate (ART), primarily utilized for malaria treatment, was shown to inhibit M1 macrophage differentiation and pro-inflammatory cytokine levels via upregulating the PHGDH expression, thereby attenuating the disease activity in SLE.

Current annotation strategies for T cell phenotyping of single-cell RNA-seq data.

Single-cell RNA sequencing (scRNA-seq) has become a popular technique for interrogating the diversity and dynamic nature of cellular gene expression and has numerous advantages in immunology. For example, scRNA-seq, in contrast to bulk RNA sequencing, can discern cellular subtypes within a population, which is important for heterogenous populations such as T cells. Moreover, recent advancements in the technology allow the parallel capturing of the highly diverse T-cell receptor (TCR) sequence with the gene expression. However, the field of single-cell RNA sequencing data analysis is still hampered by a lack of gold-standard cell phenotype annotation. This problem is particularly evident in the case of T cells due to the heterogeneity in both their gene expression and their TCR. While current cell phenotype annotation tools can differentiate major cell populations from each other, labelling T-cell subtypes remains problematic. In this review, we identify the common automated strategy for annotating T cells and their subpopulations, and also describe what crucial information is still missing from these tools.

Senescent skeletal muscle fibroadipogenic progenitors recruit and promote M2 polarization of macrophages.

Senescent cells compromise tissue structure and function in older organisms. We recently identified senescent fibroadipogenic progenitors (FAPs) with activated chemokine signaling pathways in the skeletal muscle of old mice, and hypothesized these cells may contribute to the age-associated accumulation of immune cells in skeletal muscle. In this study, through cell-cell communication analysis of skeletal muscle single-cell RNA-sequencing data, we identified unique interactions between senescent FAPs and macrophages, including those mediated by Ccl2 and Spp1. Using mouse primary FAPs invitro, we verified increased expression of Ccl2 and Spp1 and secretion of their respective proteins in the context of both irradiation- and etoposide-induced senescence. Compared to non-senescent FAPs, the medium of senescent FAPs markedly increased the recruitment of macrophages in an invitro migration assay, an effect that was mitigated by preincubation with antibodies to either CCL2 or osteopontin (encoded by Spp1). Further studies demonstrated that the secretome of senescent FAPs promotes polarization of macrophages toward an M2 subtype. These data suggest the unique secretome of senescent FAPs may compromise skeletal muscle homeostasis by recruiting and directing the behavior of macrophages.

Stiffness-induced cancer-associated fibroblasts are responsible for immunosuppression in a platelet-derived growth factor ligand-dependent manner.

Pancreatic ductal adenocarcinoma (PDAC) is associated with a vast stromal reaction that arises mainly from cancer-associated fibroblasts (CAFs) and promotes both immune escape and tumor growth. Here, we used a mouse model with deletion of the activin A receptor ALK4 in the context of the KrasG12D mutation, which strongly drives collagen deposition that leads to tissue stiffness. By ligand-receptor analysis of single-cell RNA-sequencing data, we identified that, in stiff conditions, neoplastic ductal cells instructed CAFs through sustained platelet-derived growth factor (PDGF) signaling. Tumor-associated tissue rigidity resulted in the emergence of stiffness-induced CAFs (siCAFs) in vitro and in vivo. Similar results were confirmed in human data. siCAFs were able to strongly inhibit CD8+ T-cell responses in vitro and in vivo, promoting local immunosuppression. More importantly, targeting PDGF signaling led to diminished siCAF and reduced tumor growth. Our data show for the first time that early paracrine signaling leads to profound changes in tissue mechanics, impacting immune responses and tumor progression. Our study highlights that PDGF ligand neutralization can normalize the tissue architecture independent of the genetic background, indicating that finely tuned stromal therapy may open new therapeutic avenues in pancreatic cancer.

Comprehensive analysis of single-cell and bulk RNA-sequencing data identifies B cell marker genes signature that predicts prognosis and analysis of immune checkpoints expression in head and neck squamous cell carcinoma

Recent studies have shown that B cells and the associated tertiary lymphoid structures (TLS) correlate with the response of patients to immune checkpoint inhibitors (ICIs) and predict overall survival (OS) in cancer patients. We screened 145 B cell marker genes (BCMG) by a comprehensive analysis of single-cell RNA-sequencing (scRNA-seq) data of head and neck squamous cell carcinoma (HNSC) from the Gene Expression Omnibus (GEO) database. The BCMG signature (BCMGS) was established using The Cancer Genome Atlas (TCGA) dataset of HNSC and verified in four independent datasets. The multivariate Cox regression analysis identified the signature as an independent prognostic factor. A prognostic nomogram was constructed with independent prognostic factors using the TCGA dataset. GO and KEGG analysis revealed the underlying signaling pathways related to this signature. Study of immune profiles showed that patients in the low-risk group presented discriminative immune-cell infiltrations. Furthermore, the low-risk group was featured by higher TCR and BCR diversity, which suggested that low-risk patients may be more sensitive to ICIs. Immunohistochemistry was performed, and we found that high expression of FTH1 was significantly correlated with poor OS (P = 0.025). The expression of TIM-3, LAG-3 and PD-1 was positively correlated and associated with better OS in HNSC. However, there was no statistically significant difference between PD-L1, PD-L2, CTLA-4, TIGIT and prognosis. The BCMGS was a promising prognostic biomarker in HNSC, which may help to interpret the responses to immunotherapy and provide a new perspective for future research on the treatment in HNSC.

Gut microbiome dysbiosis in men who have sex with men increases HIV infection risk through immunity homeostasis alteration.

Recent studies pointed out that gut microbiome dysbiosis in HIV infection was possibly confounded in men who have sex with men (MSM), but there is a lack of evidence. It also remained unclear how MSM-associated gut microbiome dysbiosis affected human health. This study aimed to compare the differences in gut microbiome changes between HIV and MSM and reveal the potential impacts of MSM-associated gut microbiome dysbiosis on the immune system. We searched available studies based on the PubMed database, and all gut microbiome changes associated with HIV infection and MSM were extracted from the enrolled studies. The gutMgene database was used to identify the target genes and metabolites of the gut microbiome. Bioinformatic technology and single-cell RNA sequencing data analysis were utilized to explore the impacts of these gut microbiome changes on human immunity. The results showed significant overlaps between the gut microbiome associated with HIV and that of MSM. Moreover, bioinformatic analysis revealed that gut microbiome dysbiosis in MSM had an impact on several pathways related to immunity, including the IL-17 signaling pathway and Th17 cell differentiation. Additionally, target genes of MSM-associated gut microbiome were found to be highly expressed in monocytes and lymphocytes, suggesting their potential regulatory role in immune cells. Furthermore, we found that MSM-associated gut microbiome could produce acetate and butyrate which were reported to increase the level of inflammatory factors. In conclusion, this study highlighted that MSM-associated gut microbiome dysbiosis might increase the risk of HIV acquisition by activating the immune system. Further studies are expected to elucidate the mechanism by which gut microbiome dysbiosis in MSM modulates HIV susceptibility.

Evaluation of deep learning-based feature selection for single-cell RNA sequencing data analysis

BackgroundFeature selection is an essential task in single-cell RNA-seq (scRNA-seq) data analysis and can be critical for gene dimension reduction and downstream analyses, such as gene marker identification and cell type classification. Most popular methods for feature selection from scRNA-seq data are based on the concept of differential distribution wherein a statistical model is used to detect changes in gene expression among cell types. Recent development of deep learning-based feature selection methods provides an alternative approach compared to traditional differential distribution-based methods in that the importance of a gene is determined by neural networks.ResultsIn this work, we explore the utility of various deep learning-based feature selection methods for scRNA-seq data analysis. We sample from Tabula Muris and Tabula Sapiens atlases to create scRNA-seq datasets with a range of data properties and evaluate the performance of traditional and deep learning-based feature selection methods for cell type classification, feature selection reproducibility and diversity, and computational time.ConclusionsOur study provides a reference for future development and application of deep learning-based feature selection methods for single-cell omics data analyses.

A novel missense variant in CDK5RAP2 associated with non-obstructive azoospermia

ObjectiveThe most severe type of male infertility is non-obstructive azoospermia (NOA), where there is no sperm in the ejaculate due to failure of spermatogenesis, affecting 10%–20% of infertile men with azoospermia. Genetic studies have identified dozens of NOA genes. The main aim of the present study is to identify a novel monogenic mutation that may cause NOA. Materials and methodsWe studied the pedigree of a consanguineous family with three NOA and one fertile brother by a family-based exome-sequencing, segregation analysis, insilico protein modeling and single-cell RNA sequencing data analysis. ResultsBioinformatics analysis followed by sanger sequencing revealed that three NOA brothers were homozygous for a rare missense variant in Cyclin Dependent Kinase Regulatory Subunit Associated Protein 2 (Centrosomin) CDK5RAP2 (NM_018249:exon26:c.A4003T:p.R1335W, rs761196443). Protein modeling demonstrated that CDK5RAP2, Arg1335Trp resided nearby the Microtubule Associated Protein RP/EB Family Member 1 (EB1/MAPRE1) interaction site. As a consequence of the R1335W mutation, the positively charged Arginine was replaced by to the hydrophobic tryptophan residue, possibly leading to local instability in the structure and perturbation in the CDK5RAP2-MAPRE1 interaction. ConclusionOur study reports a novel missense variant of CDK5RAP2 that segregates in homozygosity with male infertility and NOA in a consanguineous family. In silico structural predictions and gene expression data indicate a potential role of the CDK5RAP2 variant in causing defective centrosomic maturation during spermatogenesis.

Integrated analysis of single-cell and bulk RNA-sequencing data reveals the prognostic value and molecular function of THSD7A in gastric cancer.

The biological role and prognostic value of thrombospondin domain-containing 7A (THSD7A) in gastric cancer remain unclear. Our purpose was to determine the molecular mechanisms underlying the functioning of THSD7A and its prognostic value in gastric cancer. Gastric cancer-associated single cell and bulk RNA sequencing data obtained from two databases, were analyzed. We used bulk RNA sequencing to examine the differential expression of THSD7A in gastric cancer and normal gastric tissues and explored the relationship between THSD7A expression and clinicopathological characteristics. Kaplan-Meier survival and Cox analyses revealed the prognostic value of THSD7A. Gene set enrichment and immune infiltration analyses were used to determine the cancer-promoting mechanisms of THSD7A and its effect on the immune microenvironment. We explored the relationship between THSD7A expression and sensitivity of anti-tumor drugs and immune checkpoint levels. Biological functions of THSD7A were validated at single-cell and in vitro levels. THSD7A expression was significantly increased in gastric cancer samples. High THSD7A expression was associated with poor clinical phenotypes and prognoses. Cox analysis showed that THSD7A was an independent risk factor for patients with gastric cancer. Enrichment analysis suggested that epithelial-mesenchymal transition and inflammatory responses may be potential pro-cancer mechanisms of THSD7A. Upregulation of THSD7A promoted infiltration by M2 macrophages and regulatory T cells. High THSD7A expression suppressed the sensitivity of patients with gastric cancer to drugs, such as 5-fluorouracil, bleomycin, and cisplatin, and upregulated immune checkpoints, such as HAVCR2, PDCD1LG2, TIGIT, and CTLA4. At the single cell level, THSD7A was an endothelial cell-associated gene and endothelial cells overexpressing THSD7A showed unique pro-oncogenic effects. In vitro experiments confirmed that THSD7A was overexpressed in gastric cancer samples and cells, and that knocking out THSD7A significantly inhibited gastric cancer cell proliferation and invasion. THSD7A overexpression may be a unique prognostic marker and therapeutic target in gastric cancer. Therefore, our study provides a new perspective on the precise treatment of gastric cancer.

Comparative immunohistochemical analysis suggests a conserved role of EPS8L1 in epidermal and hair follicle barriers of mammals.

The mammalian skin and its appendages depend on tightly coordinated differentiation of epithelial cells. Epidermal growth factor receptor (EGFR) pathway substrate 8 (EPS8) like 1 (EPS8L1) is enriched in the epidermis among human tissues and has also been detected in the epidermis of lizards. Here, we show by the analysis of single-cell RNA-sequencing data that EPS8L1 mRNA is co-expressed with filaggrin and loricrin in terminally differentiated human epidermal keratinocytes. Comparative genomics indicated that EPS8L1 is conserved in all main clades of mammals, whereas the orthologous gene has been lost in birds. Using a polyclonal antibody against EPS8L1, we performed an immunohistochemical screening of skin from diverse mammalian species and immuno-electron microscopy of human skin. EPS8L1 was detected predominantly in the granular layer of the epidermis in monotremes, marsupial, and placental mammals. The labeling was partly associated with cell membranes, and it was evident along the perimeter of keratinocytes at the transition with the cornified layer of the epidermis, similar to involucrin distribution. Basal, spinous, and the fully mature cornified layers lacked immunolabeling of EPS8L1. In addition to the epidermis, the hair follicle inner root sheath (IRS) was immunolabeled. Both epidermal granular layer and IRS contribute to the barrier function of the skin, suggesting that EPS8L1 is involved in the regulation of these barriers.

Prediction of cell-cell communication patterns of dorsal root ganglion cells: single-cell RNA sequencing data analysis.

Abstract JOURNAL/nrgr/04.03/01300535-202406000-00042/inline-graphic1/v/2023-10-24T010719Z/r/image-tiff Dorsal root ganglion neurons transmit peripheral somatic information to the central nervous system, and dorsal root ganglion neuron excitability affects pain perception. Dorsal root ganglion stimulation is a new approach for managing pain sensation. Knowledge of the cell-cell communication among dorsal root ganglion cells may help in the development of new pain and itch management strategies. Here, we used the single-cell RNA-sequencing (scRNA-seq) database to investigate intercellular communication networks among dorsal root ganglion cells. We collected scRNA-seq data from six samples from three studies, yielding data on a total of 17,766 cells. Based on genetic profiles, we identified satellite glial cells, Schwann cells, neurons, vascular endothelial cells, immune cells, fibroblasts, and vascular smooth muscle cells. Further analysis revealed that eight types of dorsal root ganglion neurons mediated proprioceptive, itch, touch, mechanical, heat, and cold sensations. Moreover, we predicted several distinct forms of intercellular communication among dorsal root ganglion cells, including cell-cell contact, secreted signals, extracellular matrix, and neurotransmitter-mediated signals. The data mining predicted that Mrgpra3-positive neurons robustly express the genes encoding the adenosine Adora2b (A2B) receptor and glial cell line-derived neurotrophic factor family receptor alpha 1 (GFRα-1). Our immunohistochemistry results confirmed the coexpression of the A2B receptor and GFRα-1. Intrathecal injection of the A2B receptor antagonist PSB-603 effectively prevented histamine-induced scratching behaviour in a dose-dependent manner. Our results demonstrate the involvement of the A2B receptor in the modulation of itch sensation. Furthermore, our findings provide insight into dorsal root ganglion cell-cell communication patterns and mechanisms. Our results should contribute to the development of new strategies for the regulation of dorsal root ganglion excitability.

ZINBMM: a general mixture model for simultaneous clustering and gene selection using single-cell transcriptomic data

Clustering is a critical component of single-cell RNA sequencing (scRNA-seq) data analysis and can help reveal cell types and infer cell lineages. Despite considerable successes, there are few methods tailored to investigating cluster-specific genes contributing to cell heterogeneity, which can promote biological understanding of cell heterogeneity. In this study, we propose a zero-inflated negative binomial mixture model (ZINBMM) that simultaneously achieves effective scRNA-seq data clustering and gene selection. ZINBMM conducts a systemic analysis on raw counts, accommodating both batch effects and dropout events. Simulations and the analysis of five scRNA-seq datasets demonstrate the practical applicability of ZINBMM.

Restoring gluconeogenesis by TEF inhibited proliferation and promoted apoptosis and immune surveillance in kidney renal clear cell carcinoma

BackgroundKidney renal clear cell carcinoma (KIRC) is the major histological subtype of kidney tumor which covers approximately 80% of the cases. Although various therapies have been developed, the clinical outcome remains unsatisfactory. Metabolic dysregulation is a key feature of KIRC, which impacts progression and prognosis of the disease. Therefore, understanding of the metabolic changes in KIRC is of great significance in improving the treatment outcomes.MethodsThe glycolysis/gluconeogenesis genes were analyzed in the KIRC transcriptome from the Cancer Genome Atlas (TCGA) by the different expression genes (DEGs) test and survival analysis. The gluconeogenesis-related miRNAs were identified by ImmuLncRNA. The expression levels of indicated genes and miRNAs were validated in KIRC tumor and adjunct tissues by QPCR. The effects of miR-4477b and PCK1 on cell proliferation and apoptosis were examined using the cell viability assay, cell apoptosis assay, and clone information. The interaction of miR-4477b with TEF was tested by the luciferase report gene assay. The different gluconeogenesis statuses of tumor cells and related signatures were investigated by single-cell RNA sequencing (scRNA-seq) analysis.ResultsThe 11 gluconeogenesis genes were found to be suppressed in KIRC (referring as PGNGs), and the less suppression of PGNGs indicated better survival outcomes. Among the 11 PGNGs, we validated four rate-limiting enzyme genes in clinical tumor patients. Moreover, restoring gluconeogenesis by overexpressing PCK1 or TEF through miR-4477b inhibition significantly inhibited tumor cell proliferation, colony formation, and induced cell apoptosis in vitro. Independent single-cell RNA sequencing (scRNA-seq) data analysis revealed that the tumor cells had high levels of PGNG expression (PGNG + tumor cells) represented a phenotype of early stage of neoplasia and prompted immune surveillance.ConclusionsOur study suggests that the deficiency of gluconeogenesis is a key metabolic feature of KIRC, and restoring gluconeogenesis could effectively inhibit the proliferation and progression of KIRC cells.

Stationary-to-migratory transition in glioblastoma stem-like cells driven by a fatty acid-binding protein 7-RXRα neurogenic pathway.

Glioblastoma (GBM) stem-like cells (GSCs) are crucial drivers of treatment resistance and tumor recurrence. While the concept of "migrating" cancer stem cells was proposed a decade ago, the roles and underlying mechanisms of the heterogeneous populations of GSCs remain poorly defined. Cell migration using GBM cell lines and patient-derived GSCs was examined using Transwell inserts and the scratch assay. Single-cell RNA sequencing data analysis were used to map GSC drivers to specific GBM cell populations. Xenografted mice were used to model the role of brain-type fatty acid-binding protein 7 (FABP7) in GBM infiltration and expansion. The mechanism by which FABP7 and its fatty acid ligands promote GSC migration was examined by gel shift and luciferase gene reporter assays. A subpopulation of FABP7-expressing migratory GSCs was identified, with FABP7 upregulating SOX2, a key modulator for GBM stemness and plasticity, and ZEB1, a prominent factor in GBM epithelial-mesenchymal transition and invasiveness. Our data indicate that GSC migration is driven by nuclear FABP7 through activation of RXRα, a nuclear receptor activated by polyunsaturated fatty acids (PUFAs). Infiltrative progression in GBM is driven by migratory GSCs through activation of a PUFA-FABP7-RXRα neurogenic pathway.

Uptake of long-chain fatty acids from the bone marrow suppresses CD8+ T-cell metabolism and function in multiple myeloma

AbstractT cells demonstrate impaired function in multiple myeloma (MM) but suppressive mechanisms in the bone marrow microenvironment remain poorly defined. We observe that bone marrow CD8+ T-cell function is decreased in MM compared with controls, and is also consistently lower within bone marrow samples than in matched peripheral blood samples. These changes are accompanied by decreased mitochondrial mass and markedly elevated long-chain fatty acid uptake. In vitro modeling confirmed that uptake of bone marrow lipids suppresses CD8+ T function, which is impaired in autologous bone marrow plasma but rescued by lipid removal. Analysis of single-cell RNA-sequencing data identified expression of fatty acid transport protein 1 (FATP1) in bone marrow CD8+ T cells in MM, and FATP1 blockade also rescued CD8+ T-cell function, thereby identifying this as a novel target to augment T-cell activity in MM. Finally, analysis of samples from cohorts of patients who had received treatment identified that CD8+ T-cell metabolic dysfunction resolves in patients with MM who are responsive to treatment but not in patients with relapsed MM, and is associated with substantial T-cell functional restoration.