# Integrating single-cell RNA-seq and bulk RNA-seq to construct prognostic signatures to explore the role of glutamine metabolism in breast cancer

**Authors:** Shengbin Pei, Pengpeng Zhang, Huilin Chen, Shuhan Zhao, Yuhan Dai, Lili Yang, Yakun Kang, Mingjie Zheng, Yiqin Xia, Hui Xie

PMC · DOI: 10.3389/fendo.2023.1135297 · 2023-02-10

## TL;DR

This study combines single-cell and bulk RNA sequencing to develop a new prognostic model for breast cancer based on glutamine metabolism genes, improving survival prediction and treatment response.

## Contribution

A novel prognostic risk signature (PRS) based on glutamine metabolism-related genes is developed for breast cancer prognosis and immunotherapy response.

## Key findings

- The PRS model outperforms traditional clinicopathological features in predicting breast cancer survival.
- SNX3 gene knockdown significantly reduces cancer cell proliferation and migration in vitro.
- PRS subgroups reveal distinct immune status and mutational patterns in breast cancer.

## Abstract

Although breast cancer (BC) treatment has entered the era of precision therapy, the prognosis is good in the case of comprehensive multimodal treatment such as neoadjuvant, endocrine, and targeted therapy. However, due to its high heterogeneity, some patients still cannot benefit from conventional treatment and have poor survival prognoses. Amino acids and their metabolites affect tumor development, alter the tumor microenvironment, play an increasingly obvious role in immune response and regulation of immune cell function, and are involved in acquired and innate immune regulation; therefore, amino acid metabolism is receiving increasing attention.

Based on public datasets, we carried out a comprehensive transcriptome and single-cell sequencing investigation. Then we used 2.5 Weighted Co-Expression Network Analysis (WGCNA) and Cox to evaluate glutamine metabolism-related genes (GRGs) in BC and constructed a prognostic model for BC patients. Finally, the expression and function of the signature key gene SNX3 were examined by in vitro experiments.

In this study, we constituted a risk signature to predict overall survival (OS) in BC patients by glutamine-related genes. According to our risk signature, BC patients can obtain a Prognostic Risk Signature (PRS), and the response to immunotherapy can be further stratified according to PRS. Compared with traditional clinicopathological features, PRS demonstrated robust prognostic power and accurate survival prediction. In addition, altered pathways and mutational patterns were analyzed in PRS subgroups. Our study sheds some light on the immune status of BC. In in vitro experiments, the knockdown of SNX3, an essential gene in the signature, resulted in a dramatic reduction in proliferation, invasion, and migration of MDA-MB-231 and MCF-7 cell lines.

We established a brand-new PRS consisting of genes associated with glutamine metabolism. It expands unique ideas for the diagnosis, treatment, and prognosis of BC.

## Linked entities

- **Genes:** SNX3 (sorting nexin 3) [NCBI Gene 8724]
- **Diseases:** breast cancer (MONDO:0004989)

## Full-text entities

- **Genes:** MMADHC (metabolism of cobalamin associated D) [NCBI Gene 27249] {aka C2orf25, CL25022, HMAD, MACD, MAHCD, cblD}, PIK3CA (phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit alpha) [NCBI Gene 5290] {aka CCM4, CLAPO, CLOVE, CWS5, HMH, MCAP}, LIMCH1 (LIM and calponin homology domains 1) [NCBI Gene 22998] {aka LIMCH1A, LMO7B}, UQCRFS1 (ubiquinol-cytochrome c reductase, Rieske iron-sulfur polypeptide 1) [NCBI Gene 7386] {aka ISP, MC3DN10, RIP1, RIS1, RISP, UQCR5}, EREG (epiregulin) [NCBI Gene 2069] {aka EPR, ER, Ep}, TMEM14B (transmembrane protein 14B) [NCBI Gene 81853], FABP7 (fatty acid binding protein 7) [NCBI Gene 2173] {aka B-FABP, BLBP, FABPB, MRG}, SLC3A2 (solute carrier family 3 member 2) [NCBI Gene 6520] {aka 4F2, 4F2HC, 4T2HC, CD98, CD98HC, MDU1}, USH2A (usherin) [NCBI Gene 7399] {aka RP39, US2, USH2, dJ1111A8.1}, CD4 (CD4 molecule) [NCBI Gene 920] {aka CD4mut, IMD79, Leu-3, OKT4D, T4}, MYC (MYC proto-oncogene, bHLH transcription factor) [NCBI Gene 4609] {aka MRTL, MYCC, bHLHe39, c-Myc}, CYSTM1 (cysteine rich transmembrane module containing 1) [NCBI Gene 84418] {aka C5orf32, ORF1-FL49}, DMD (dystrophin) [NCBI Gene 1756] {aka BMD, CMD3B, DXS142, DXS164, DXS206, DXS230}, BTG1 (BTG anti-proliferation factor 1) [NCBI Gene 694] {aka APRO2}, NFKB1 (nuclear factor kappa B subunit 1) [NCBI Gene 4790] {aka CVID12, EBP-1, KBF1, NF-kB, NF-kB1, NF-kappa-B1}, TLE5 (TLE family member 5, transcriptional modulator) [NCBI Gene 166] {aka AES, AES-1, AES-2, ESP1, GRG, GRG5}, CCNA2 (cyclin A2) [NCBI Gene 890] {aka CCN1, CCNA}, SLC38A2 (solute carrier family 38 member 2) [NCBI Gene 54407] {aka ATA2, PRO1068, SAT2, SNAT2}, CDH1 (cadherin 1) [NCBI Gene 999] {aka Arc-1, BCDS1, CD324, CDHE, ECAD, LCAM}, JAK1 (Janus kinase 1) [NCBI Gene 3716] {aka AIIDE, JAK1A, JAK1B, JTK3}, TP53 (tumor protein p53) [NCBI Gene 7157] {aka BCC7, BMFS5, LFS1, P53, TRP53}, SNX3 (sorting nexin 3) [NCBI Gene 8724] {aka Grd19, MCOPS8, SDP3}, GATA3 (GATA binding protein 3) [NCBI Gene 2625] {aka HDR, HDRS}, PTEN (phosphatase and tensin homolog) [NCBI Gene 5728] {aka 10q23del, BZS, CWS1, DEC, GLM2, MHAM}, PKIB (cAMP-dependent protein kinase inhibitor beta) [NCBI Gene 5570] {aka PRKACN2}, RBP1 (retinol binding protein 1) [NCBI Gene 5947] {aka CRABP-I, CRBP, CRBP1, CRBPI, RBPC, hCRBP1}, SLC7A11 (solute carrier family 7 member 11) [NCBI Gene 23657] {aka CCBR1, xCT}, KDELR2 (KDEL endoplasmic reticulum protein retention receptor 2) [NCBI Gene 11014] {aka ELP-1, ERD2.2, OI21}, IFNG (interferon gamma) [NCBI Gene 3458] {aka IFG, IFI, IMD69}, KMT2C (lysine methyltransferase 2C) [NCBI Gene 58508] {aka HALR, KLEFS2, MLL3}, GLUL (glutamate-ammonia ligase) [NCBI Gene 2752] {aka DEE116, GLNS, GS, PIG43, PIG59}, EI24 (EI24 autophagy associated transmembrane protein) [NCBI Gene 9538] {aka EPG4, PIG8, TP53I8}, LINC01614 (long intergenic non-protein coding RNA 1614) [NCBI Gene 105373869] {aka LCAL4}, GAPDH (glyceraldehyde-3-phosphate dehydrogenase) [NCBI Gene 2597] {aka G3PD, GAPD, HEL-S-162eP}, STK17A (serine/threonine kinase 17a) [NCBI Gene 9263] {aka DRAK1}, HMCN1 (hemicentin 1) [NCBI Gene 83872] {aka ARMD1, FBLN6, FIBL-6, FIBL6}, E2F1 (E2F transcription factor 1) [NCBI Gene 1869] {aka E2F-1, RBAP1, RBBP3, RBP3}, ERRFI1 (ERBB receptor feedback inhibitor 1) [NCBI Gene 54206] {aka GENE-33, MIG-6, MIG6, RALT}, RPL17 (ribosomal protein L17) [NCBI Gene 6139] {aka DBA22, L17, PD-1, RPL23, uL22}, STAT3 (signal transducer and activator of transcription 3) [NCBI Gene 6774] {aka ADMIO, ADMIO1, APRF, HIES}, MAP3K1 (mitogen-activated protein kinase kinase kinase 1) [NCBI Gene 4214] {aka MAPKKK1, MEKK, MEKK 1, MEKK1, SRXY6}, GLS (glutaminase) [NCBI Gene 2744] {aka AAD20, CASGID, DEE71, EIEE71, GAC, GAM}, SLC38A1 (solute carrier family 38 member 1) [NCBI Gene 81539] {aka ATA1, NAT2, SAT1, SNAT1}, ICOSLG (inducible T cell costimulator ligand) [NCBI Gene 23308] {aka B7-H2, B7H2, B7RP-1, B7RP1, B7h, CD275}, TNF (tumor necrosis factor) [NCBI Gene 7124] {aka DIF, IMD127, TNF-alpha, TNFA, TNFSF2, TNLG1F}, MMP7 (matrix metallopeptidase 7) [NCBI Gene 4316] {aka MMP-7, MPSL1, PUMP-1}, GPX4 (glutathione peroxidase 4) [NCBI Gene 2879] {aka GPx-4, GSHPx-4, MCSP, PHGPx, SMDS, snGPx}, NDUFB9 (NADH:ubiquinone oxidoreductase subunit B9) [NCBI Gene 4715] {aka B22, CI-B22, LYRM3, MC1DN24, UQOR22}, IGKC (immunoglobulin kappa constant) [NCBI Gene 3514] {aka HCAK1, IGKCD, Km}, MTOR (mechanistic target of rapamycin kinase) [NCBI Gene 2475] {aka FRAP, FRAP1, FRAP2, RAFT1, RAPT1, SKS}, CD274 (CD274 molecule) [NCBI Gene 29126] {aka ADMIO5, B7-H, B7H1, PD-L1, PDCD1L1, PDCD1LG1}
- **Diseases:** precancerous (MESH:D011230), neurodegenerative diseases (MESH:D019636), colorectal cancer (MESH:D015179), liver cancer (MESH:D006528), osteosarcoma (MESH:D012516), lung cancer (MESH:D008175), TNBC (MESH:D064726), autoimmune diseases (MESH:D001327), metastasis (MESH:D009362), inflammatory (MESH:D007249), gastric cancer (MESH:D013274), glutamine (MESH:C536832), TIDE (MESH:D007154), respiratory diseases (MESH:D012140), EdU (OMIM:616087), BC (MESH:D001943), Cancer (MESH:D009369)
- **Chemicals:** platinum (MESH:D010984), branched-chain amino acids (MESH:D000597), penicillin (MESH:D010406), streptomycin (MESH:D013307), lipid (MESH:D008055), CO2 (MESH:D002245), SYBR Green (MESH:C098022), glutamate (MESH:D018698), valine (MESH:D014633), Triton-X-100 (MESH:D017830), cysteine (MESH:D003545), nitrogen (MESH:D009584), Lipofectamine 3000 (-), TRIzol (MESH:C411644), PBS (MESH:D007854), Apollo (MESH:C062106), glutathione (MESH:D005978), Amino acids (MESH:D000596), carbon (MESH:D002244), palivizumab (MESH:D000069455), atezolizumab (MESH:C000594389), paraformaldehyde (MESH:C003043), glucose (MESH:D005947), iron (MESH:D007501), Gln (MESH:D005973), leucine (MESH:D007930), NADPH (MESH:D009249), 5-Ethynyl-2'-deoxyuridine (MESH:C031086), Crystal violet (MESH:D005840)
- **Species:** Homo sapiens (human, species) [taxon 9606]
- **Mutations:** glutamine to glutamate
- **Cell lines:** MCF-7 — Homo sapiens (Human), Invasive breast carcinoma of no special type, Cancer cell line (CVCL_0031), MDA-MB-231 — Homo sapiens (Human), Breast adenocarcinoma, Cancer cell line (CVCL_0062)

## Figures

10 figures with captions in the complete paper: https://tomesphere.com/paper/PMC9950399/full.md

---
Source: https://tomesphere.com/paper/PMC9950399