Abstract
Introduction:
Aplastic anemia (AA) is a bone marrow hematopoietic failure syndrome mediated by immune cells. The mechanism of this immune disorder is not well understood and therapeutic strategies still need to be improved.
Methods:
Studies have found that abnormalities in metabolisms promote the survival of AA cells. In recent years, an increasing number of studies have reported the immunosuppressive therapy for the treatment of AA. In this study, we analyzed the transcriptome of AA from peripheral blood compared with healthy donors by single-cell sequencing and identified the affected metabolic pathways including lysine degradation. We demonstrated that the metabolic abnormalities of T lymphocytes mainly focus on glycolysis/gluconeogenesis. In addition, the metabolic abnormalities of natural killer cells concentrated in oxidative phosphorylation.
Results:
The key genes involved in abnormal metabolic processes were Neustein neurotrophic factor (NENF), inositol polyphosphate-4-phosphatase type II B (INPP4B), aldo-keto reductase family 1, member C3 (AKR1C3), and carbohydrate (N-acetylglucosamine-6-O) sulfotransferase 2 (CHST2) by differential gene expression analysis.
Discussion:
Molecule interaction analysis showed that tumor necrosis factor superfamily, member 12 (TNFSM12) in tumor necrosis factor (TNF) signaling was broadly activated in AA. In conclusion, we suppose that the treatment of the immune cells’ abnormal metabolic pathway may contribute to the development of novel strategies to treat AA.
Introduction
Aplastic anemia (AA) is a bone marrow hematopoietic failure syndrome caused by various etiologies. It is characterized by a decreased proliferation of bone marrow hematopoietic cells and peripheral blood pancytopenia (1). The main clinical manifestations are anemia, hemorrhage, and infection (2). The pathogenesis of AA is complicated, including the abnormality of the hematopoietic microenvironment (3), deficiency of hematopoietic stem cells/progenitor cells (4), and disorders of the immune system (5). AA occurs at any age. However, compared with adults, a large proportion of children with AA have a relatively high incidence rate of the bone marrow failure syndrome (6). The incidence of AA varies geographically, which is two-to-three times higher in Asia than in the western world (7).
Most acquired AA is considered to be the destruction of bone marrow hematopoietic cells mediated by T cells (8). Early studies showed that removing lymphocytes from the bone marrow with AA could increase the number of cell colonies in tissue culture, while adding the same lymphocytes to a normal bone marrow would inhibit in vitro hematopoiesis (9). In clinic, human leukocyte antigen (HLA)-matched sibling bone marrow transplantation is the first-line treatment for AA patients under 40 years old (10). The combined immunosuppressive therapy of eltrombopag, thymoglobulin, and cyclosporine A is the initial treatment of refractory AA, and patients’ survival rate is approximately 90% (11).
Cell metabolism promotes the absorption of nutrients and various components required for cell synthesis (12), enabling organisms to grow and reproduce, maintain their structure, and react to the external environment (13). In terms of obtaining energy and biosynthesis, aerobic glycolysis is very important for cell proliferation (14). The transformation of intracellular metabolic pathways in immune cells alters the function of immune cells (15). In immune cells, six metabolic pathways were intensively discussed: glycolysis, the tricarboxylic acid cycle, the pentose phosphate pathway, fatty acid oxidation, fatty acid synthesis, and amino acid metabolism (15). Amino acid metabolism plays an important role in regulating the function of innate and adaptive immune systems. Previous investigations indicate that the deletion of the transporter Alanine/Serine/Cysteine Transporter 2 (ASCT2) (responsible for glutamine and leucine uptake) gene in T cells will damage the function of helper T cell 1 (TH1) and helper T cell 17 (TH17) cells (16). The transition of the metabolic pathway from oxidative phosphorylation to aerobic glycolysis is a sign of T-cell activation and a key step to satisfy the metabolic requirements in the process of cell proliferation (17). Studies have shown that the glucose analog 2-deoxyglucose, an inhibitor of the glycolysis pathway, inhibits T cells from developing into TH17 cells (18). Therefore, metabolic reprogramming is critical for T-cell activation and functional execution (19). In order to promote drug uptake and enhance the delivery ability to target T-cell populations, we coupled them with glucose transporters (20). Therefore, regulating the metabolism of T cells may be a therapeutic means to treat AA.
Single-cell messenger RNA (mRNA) sequencing (scRNA-seq) is a technology for an unbiased, high-throughput, and high-resolution transcriptome analysis of cell heterogeneity in populations (21). It aggregates cells and identifies new subsets, as well as gene expression in various tissues (22). We created a thorough transcriptional map of immune cells from healthy controls and patients with AA by single-cell RNA sequencing. Then, we explored the changes of the cellular transcriptome in patients and identified the key metabolic pathways that may affect the occurrence of AA.
Materials and methods
Clinical samples
The transcriptomic profile was obtained from peripheral blood from five children with AA (AA: LJX-AA, LZL-AA, SLT-AA, WJL-AA, and XF-AA) and three healthy donors (Ctrl: CBC-Ctrl, CRC-Ctrl, and LJJ-Ctrl) (Figure S1). All samples were recruited via The Seventh Affiliated Hospital of Sun Yat-sen University. All participants provided written informed consent before inclusion in the study.
Single-cell mRNA library preparation and sequencing
The complementary DNA (cDNA)/DNA/small RNA libraries were sequenced on the Illumina sequencing platform by Genedenovo Biotechnology Co., Ltd (Guangzhou, China). Cellranger was used to remove the reads with low sequencing quality and then compare them with the reference genome to annotate as a specific gene. After unique molecular identifiers (UMI) correction and statistics, the unfiltered feature barcode matrix was obtained. The cells in the data were identified and distinguished according to the unfiltered feature barcode array. To filtrate multicellular samples, Doublet-Finder was applied to calculate the gel beads in emulsion (GEM) probability [pattern analysis and neural networks (pANN) value]. In addition, we used the following indicators to perform cell filtration: the number of genes identified in a single cell (200.0–3,600.0), the total number of UMI in a single cell (<17,000.0), and the proportion of mitochondrial gene expression in a single cell (<25.0%).
Data analysis and visualization
We conducted Harmony for data consolidation and batch effect correction. Dimension reduction, cell clustering, and differential gene expression were performed using the Seurat package. Based on the subset information of cells, we set | log2FC | ≥ 0.36 and the proportion of cells expressing target genes in each group ≥ 0.1 as the threshold. Then, we used the modeling and simulation team (MAST) obstacle model to test the significance of differences. We screened the pathway with the enrichment degree of top 20 in kyoto encylopaedia of genes and genomes (KEGG) A class Metabolism from the KEGG enrichment of clusters (Ctrl vs. AA). We calculated the number of differential genes in each term in the Gene Ontology (GO) database (http://www.geneontology.org/) and applied a hypergeometric test to find the GO entries that are significantly enriched. The Kyoto Encyclopedia of Genes and Genomes pathway was performed by Omicshare tools (http://www.omicshare.com/tools/). Gene Set Enrichment Analysis (GSEA) was performed by using the software GSEA and MSigDB (23). Disease Ontology Analysis was performed by the disease ontology (DO) database (http://disease-ontology.org/). The CellPhoneDB package was used to estimate cell–cell communication.
Statistical analysis
Data visualizing and statistical analysis were performed using GraphPad Prism 8.0 (GraphPad Software Inc, CA, USA). Differences between experimental groups were analyzed using unpaired Student’s t-test. P value < 0.05 was considered significant.
Results
Single-cell analysis and cell type identification
To interrogate the metabolism of immune cells in patients with AA, clinical peripheral blood specimens were analyzed by 10x Genomics based on scRNA-seq [single-cell tagged reverse transcription sequencing (STRT-seq)]. The cells were labeled and differential genes were analyzed by Seurat to complete the statistics and distribution mapping. Then, we used the GO database and KEGG database to analyze the enrichment of divergence genes (Figure 1A). We distinguished cell subsets by immune cells’ specific surface markers. Among the 19 subsets in the bubble plot, except for clusters 5 and 18, which are NK cells, the rest are T cells, including CD4 T cells, CD8 T cells, and Treg cells (Figure 1B). The percentage of different clusters showed that in AA, cluster 1 and cluster 2 are significantly distinct with the control group (Figure 1C).
Figure 1
The initial dimension reduction and unsupervised clustering of single-cell transcriptomes classified cells into 19 groups (Figures 2A, B). Compared with the control group, the difference of cluster 3 and cluster 5 is rather obvious in AA. At the same time, we conducted GO enrichment analysis from three aspects (biological process, cell component, and molecular function). The research found that, in the biological process part, genes upregulated in the AA group were significantly more than those downregulated (Figure 2C). Further study on GO enrichment analysis showed that, compared with the control group, most of the first 20 GO terms in the AA group were related to the regulation of cell metabolism (Figure 2D). Therefore, the abnormal regulation of cell metabolism and abnormal expression of related genes may be associated with the occurrence of AA.
Figure 2
Single-cell mRNA sequencing revealed metabolic differences in aplastic anemia
In order to further understand which metabolic processes are unusually regulated and metabolism-related genes are abnormally expressed in patients with AA, we conducted the KEGG pathway and GSEA. So as to display the distribution characteristics of metabolic pathways in different cell subsets, we manufactured the t-distributedstochastic neighbor embedding (t-SNE) map through the R package. The soft k-means clustering algorithm is used to cluster the dimension-reduced data, and cells were clustered into 19 cell types (Figure 3A). In the t-SNE difference distribution map of metabolism pathways, it was intuitively seen that lysine degradation was more common in AA (Figure 3B). We manufactured the homogenized gene expression into a t-SNE map. The results showed that the expression of the Neustein neurotrophic factor (NENF) gene in AA was significantly higher than that in the control group (Figure 3C).
Figure 3
Compared with the healthy donors, analysis found that the dysregulation of lysine degradation in patients with AA had marked statistical significance (AA vs. Ctrl, q = 0.0051) (Figure 4A). Interestingly, in the GSEA, we found that valine, leucine, and isoleucine degradation–related genes were more abundant in healthy donors (Figure 4B), which may indicate that valine, leucine, and isoleucine accumulated multiple times in the AA patients. To sum up, we concluded that the pathobolism of amino acids may play a key important role in the occurrence of AA. Next, we further analyzed from the aspect of metabolizing gene expression. In the heat map, there are significantly more abnormal regulated genes in cluster 3 and cluster 5 in AA compared with the control group (Figure 4C). Seurat software was performed to analyze the divisions between the subset of cells. The results were similar to the said GO analysis. The genes upregulated in AA were significantly more than those downregulated. At the same time, the metabolized gene expression in the third and fifth subsets was more obvious than that in other subsets (Figure 4D). These findings provided clues for our follow-up research on the mechanism of AA.
Figure 4
T-lymphocyte metabolism analysis
As previously mentioned, compared with the control group, the number of abnormally expressed genes in the fifth subset of the AA group is the highest among all subsets. Then, we analyzed the changes of metabolic pathways in the fifth subgroup through the KEGG pathway, listed the metabolic pathways of top 20, and drew dot plots. It was found that the glycolysis/gluconeogenesis metabolic pathways were the most significant (Figure 5A). Because most cases of AA are an immune system disorder disease mediated by T cells, this article mainly studies the divergences of immune cell metabolism between the AA group and the control group. Then, we drew a t-SNE map according to the expressed genes through the R package and found that the INPP4B gene was specifically expressed in T cells but almost not in NK cells (Figure 5B). To check out which diseases may be caused by abnormal metabolic T cells in the process of AA progression, we conducted Disease Ontology analysis, and the results showed that AA was most likely to be converted to leukemia (Figure 5C). The said results provide an important basis for our follow-up research on the treatment and prevention of AA.
Figure 5
Natural killer cell metabolism analysis
It was reported that the dysfunction of natural killer (NK) cells may also be related to the occurrence of AA. Then, we analyzed the changes of metabolic pathways in NK cells through the KEGG pathway and found that the oxidative phosphorylation in NK cells in AA was abnormal, which had statistical significance (Figure 6A). We continued to explore which metabolism-related genes are abnormally expressed in NK cells and drew the t-SNE map. Interestingly, the expression amount of the AKR1C3 gene and CHST2 gene in NK cells is significantly higher than that in other cell subsets (Figures 6B, C). CHST2 (carbohydrate (N-acetylglucosamine-6-O) sulfotransferase 2) gene encodes a sulfotransferase protein, which catalyzes the sulfation of non-reducing n-acetylglucosamine residues and participates in the metabolism of lymphocytes at the inflammatory sites (24). Similar to what was mentioned earlier, we also did the Disease Ontology analysis of NK cells, listed the metabolic pathways of top 20, and drew bar charts. It was found that the metabolic abnormalities of NK cells were the most likely to cause AA to revert to cancer and the disease of cell promotion (Figure 6D). The intercellular interaction network between 18 subsets showed that cluster 5 had the strongest correlation among all cell subsets in the AA group (Figure 6F), while cluster 4 had the strongest correlation among all cell subsets in the control group (Figure 6E), indicating that the communication relationship of subgroup 5 in the AA group is enhanced. Among molecule interactions, the expression of the TNF signaling pathway, mechanistic target of rapamycin (mTOR) signaling pathway, and PI3K-Akt signaling pathway in AA (Figure 6H) and control (Figure 6G) had the most obvious difference. We further observed that TNFSM12 in TNF signaling was broadly activated in AA, which might contribute to the high reduction of normal blood cells. The metabolic correlation analysis of NK cells in AA provides a new perspective for us to study the mechanism of AA.
Figure 6
Discussion
T cells were divided into cytotoxic T cells, helper T cells, regulatory T cells, and memory T cells according to their functions and surface markers. After being activated by antigens in the peripheral circulation, naive CD4 T cells proliferate and differentiate into various subsets of T helper cells, including Th1, Th2, and Th17 cells (25). The metabolic reprogramming of T cells enables them to shift from oxidative metabolism to biosynthetic metabolism to support rapid cell growth (26). Activated CD4 T cells require efficient glucose uptake, glycolysis, glutamate decomposition, and lipid synthesis to maintain cell proliferation (27). The enhancement of aerobic glycolysis in cells makes glucose and other nutrients not oxidized in mitochondria to produce ATP but used for the biosynthesis of nucleic acids, lipids, and amino acids (28). The boost of the glycolytic metabolic pathway occurs mostly in activated NK cells (29), T lymphocytes (30), and B lymphocytes (31). Previous studies have shown that pyruvate dehydrogenase is a key enzyme in T-cell glycolysis and oxidative metabolism (32). In our study, we found that the lysine degradation pathway of immune cells in AA is significantly higher than that in normal samples through the single-cell sequencing analysis of peripheral blood clinical samples. Lysine degradation is caused by ϵ-deamination or α-deamination reaction and produces two acetyl coenzyme A and several reductants (33). It has been found that the selective modification of lysine sites in proteins by aminophiles disrupts the interaction between proteins and RNA in the immune response (34). Our results showed that the disturbance of the immune system in AA was also related to the degradation of valine, leucine, and isoleucine. Previous studies indicated that L-amino acid transporters that are composed of Slc7a5 and CD98 induce leucine uptake, activate the mTOR pathway, and affect T-cell metabolism (35). During the activation of T cells and B cells, the transcription of intracellular glutamine transporters SNAT1 and SNAT2 is enhanced (36). It is covered that the NK cell count is decreased and the activity is impaired in patients with Fanconi anemia (37). NK cells mainly rely on oxidative phosphorylation to generate energy and activate downstream to produce interferon-γ (IFN-γ) (38). Oxidative phosphorylation and glycolysis are two major metabolic pathways for energy production and cell function maintenance. In immune cells, oxidative phosphorylation can regulate the formation of memory cells and related inflammatory reactions (39). Several recent studies in the US have shown that Cox10 (a gene encoding the composition of mitochondrial electron transfer chain complex IV) plays an important role in NK-cell antigen-specific amplification and murine cytomegalovirus (MCMV) infection (40).
At the transcriptional level, our research found that the NENF gene was upregulated in the immune cells of AA. Neudesin was initially identified as a secreted protein with neurotrophic activity. It has a conservative cytochrome 5-like heme/steroid binding domain and can activate intracellular signal pathways by binding to G protein–coupled receptors (41). NENF is essential in a variety of biological processes, including neural function, fat metabolism, and tumorigenesis (42). It has been found that neudesin inhibits adipogenesis in mouse embryonic fibroblasts cells 3T3-L1 (3T3-L1) cells through mitogen-activated protein kinase (MAPK) cascade reaction (43). We also detected that INPP48B gene was upregulated in the T cells of patients with AA, and AKR1C3 and CHST2 were specifically upregulated in NK cells. INPP4B was initially identified as an enzyme that preferentially hydrolyzes the 4-phosphate of phosphatidylinositol-3,4-bisphosphate (PI(3,4)P2), to generate phosphatidylinositol-3-phosphate (PI(3)P) (44). In recent studies, the overexpression of INPP4B in AML cells enhancescolony-forming potential and induces chemotherapy resistance in acute myelocytic leukemia (AML) patients (45). As a soluble enzyme of the aldehyde ketone reductase family, AKR1C3 plays an important role in regulating prostaglandin, the steroid hormone, and retinoic acid metabolism (46). Chst2 encodes a carbohydrate sulfotransferase that catalyzes the sulfation of the C6 position of GlcNAc during keratan sulfate biosynthesis (47).
Through disease ontology enrichment analysis, the abnormal metabolism of T cells is likely to cause AA to develop into leukemia, and the metabolic changes of NK cells are likely to lead to abnormal cell proliferation diseases and tumors. Relevant results have also confirmed that the secondary myelodysplastic syndrome and acute leukemia usually develop from severe AA after immunosuppressive therapy (48). These findings support the broad potential of targeting functional lysine in the human proteome. It has been discovered that drugs acting on the surface receptors of the CTLA-4 and PD-1 can limit the uptake of glucose and amino acids, so as to negatively regulate the activation of T cells (49). Detailed investigations on specific mechanisms of metabolic abnormalities in the immune cells of AA are needed in the future.
Statements
Data availability statement
The data presented in the study are deposited in the CNCB repository, submission of HRA: HRA003544, project: PRJCA012914.
Ethics statement
The studies involving human participants were reviewed and approved by Shenzhen Science and Technology Innovation Commission. Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin. Written informed consent was obtained from the individual(s), and minor(s)’ legal guardian/next of kin, for the publication of any potentially identifiable images or data included in this article.
Author contributions
LH, JH, and YL conceived the project. LH and LW performed the experiments. QZ, JH, and JY analyzed data. QZ wrote the manuscript. All authors contributed to the article and approved the submitted version.
Funding
The work was supported by Pediatric Hematological Oncology, Seventh Affiliated Hospital, Sun Yat-sen University. We also thank Shenzhen Basic Research Project (Natural Science Fund) (202103243003611) to support the design of the study, Shenzhen Science and Technology Innovation Commission (JCYJ20190814164601648 and JCYJ20210324123003009) to support sample collection and analysis, Sanming Project of Medicine in Shenzhen (No. SZSM202011004) and Free Exploration Project of Shenzhen Science and Technology Innovation Commission (JCYJ20180307150419435) to support the manuscript preparation and publication.
Acknowledgments
We are grateful to/thank Guangzhou Genedenovo Biotechnology Co., Ltd for assisting in sequencing and/or bioinformatics analysis.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2023.1075408/full#supplementary-material
References
1
YoungNSCaladoRTScheinbergP. Current concepts in the pathophysiology and treatment of aplastic anemia. Blood (2006) 108:2509–19. doi: 10.1182/blood-2006-03-010777
2
KillickSBBownNCavenaghJDokalIFoukaneliTHillAet al. Guidelines for the diagnosis and management of adult aplastic anaemia. Br J Haematol (2016) 172:187–207. doi: 10.1111/bjh.13853
3
MedingerMDrexlerBLengerkeCPasswegJ. Pathogenesis of acquired aplastic anemia and the role of the bone marrow microenvironment. Front Oncol (2018) 8:587. doi: 10.3389/fonc.2018.00587
4
ImiTKatagiriTHosomichiKZaimokuYNguyenVHNakagawaNet al. Sustained clonal hematopoiesis by HLA-lacking hematopoietic stem cells without driver mutations in aplastic anemia. Blood Adv (2018) 2:1000–12. doi: 10.1182/bloodadvances.2017013953
5
NissenCSternM. Acquired immune mediated aplastic anemia: Is it antineoplastic? Autoimmun Rev (2009) 9:11–6. doi: 10.1016/j.autrev.2009.02.032
6
KurrePJohnsonFLDeegHJ. Diagnosis and treatment of children with aplastic anemia. Pediatr Blood Cancer (2005) 45:770–80. doi: 10.1002/pbc.20322
7
MontaneEIbáñezLVidalXBallarínEPuigRGarcíaNet al. Epidemiology of aplastic anemia: A prospective multicenter study. Haematologica (2008) 93:518–23. doi: 10.3324/haematol.12020
8
YoungNSMaciejewskiJ. The pathophysiology of acquired aplastic anemia. N Engl J Med (1997) 336:1365–72. doi: 10.1056/NEJM199705083361906
9
YoungNS. Hematopoietic cell destruction by immune mechanisms in acquired aplastic anemia. Semin Hematol (2000) 37:3–14. doi: 10.1016/s0037-1963(00)90026-x
10
BacigalupoA. How I treat acquired aplastic anemia. Blood (2017) 129:1428–36. doi: 10.1182/blood-2016-08-693481
11
DesmondRTownsleyDMDumitriuBOlnesMJScheinbergPBevansMet al. Eltrombopag restores trilineage hematopoiesis in refractory severe aplastic anemia that can be sustained on discontinuation of drug. Blood (2014) 123:1818–25. doi: 10.1182/blood-2013-10-534743
12
LiXSunXCarmelietP. Hallmarks of endothelial cell metabolism in health and disease. Cell Metab (2019) 30:414–33. doi: 10.1016/j.cmet.2019.08.011
13
TorrinoSGrassetEMAudebertSBelhadjILacouxCHaynesMet al. Mechano-induced cell metabolism promotes microtubule glutamylation to force metastasis. Cell Metab (2021) 33:1342–1357.e1310. doi: 10.1016/j.cmet.2021.05.009
14
Vander HeidenMGCantleyLCThompsonCB. Understanding the warburg effect: The metabolic requirements of cell proliferation. Science (2009) 324:1029–33. doi: 10.1126/science.1160809
15
AndrejevaGRathmellJC. Similarities and distinctions of cancer and immune metabolism in inflammation and tumors. Cell Metab (2017) 26:49–70. doi: 10.1016/j.cmet.2017.06.004
16
NakayaMXiaoYZhouXChangJ-HChangMChengXet al. Inflammatory T cell responses rely on amino acid transporter ASCT2 facilitation of glutamine uptake and mTORC1 kinase activation. Immunity (2014) 40:692–705. doi: 10.1016/j.immuni.2014.04.007
17
ChangCHCurtisJDMaggiLBJrFaubertBVillarinoAVO'SullivanDet al. Posttranscriptional control of T cell effector function by aerobic glycolysis. Cell (2013) 153:1239–51. doi: 10.1016/j.cell.2013.05.016
18
ShiLZWangRHuangGVogelPNealeGGreenDRet al. HIF1alpha-dependent glycolytic pathway orchestrates a metabolic checkpoint for the differentiation of TH17 and treg cells. J Exp Med (2011) 208:1367–76. doi: 10.1084/jem.20110278
19
O'SullivanDPearceEL. Targeting T cell metabolism for therapy. Trends Immunol (2015) 36:71–80. doi: 10.1016/j.it.2014.12.004
20
CalvaresiECHergenrotherPJ. Glucose conjugation for the specific targeting and treatment of cancer. Chem Sci (2013) 4:2319–33. doi: 10.1039/C3SC22205E
21
FastresAPirottinDFievezLMarichalTDesmetCJBureauFet al. Characterization of the bronchoalveolar lavage fluid by single cell gene expression analysis in healthy dogs: A promising technique. Front Immunol (2020) 11:1707. doi: 10.3389/fimmu.2020.01707
22
GaoSWuZArnoldBDiamondCBatchuSGiudiceVet al. Single-cell RNA sequencing coupled to TCR profiling of large granular lymphocyte leukemia T cells. Nat Commun (2022) 13:1982. doi: 10.1038/s41467-022-29175-x
23
SubramanianATamayoPMoothaVKMukherjeeSEbertBLGilletteMAet al. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U.S.A. (2005) 102:15545–50. doi: 10.1073/pnas.0506580102
24
ChenXYangTTQiuXCJiZGLiCXLongHet al. Gene expression profiles of human osteosarcoma cell sublines with different pulmonary metastatic potentials. Cancer Biol Ther (2011) 11:287–92. doi: 10.4161/cbt.11.2.13966
25
ZhuJYamaneHPaulWE. Differentiation of effector CD4 T cell populations (*). Annu Rev Immunol (2010) 28:445–89. doi: 10.1146/annurev-immunol-030409-101212
26
SiskaPJRathmellJC. T Cell metabolic fitness in antitumor immunity. Trends Immunol (2015) 36:257–64. doi: 10.1016/j.it.2015.02.007
27
KidaniYElsaesserHHockMBVergnesLWilliamsKJArgusJPet al. Sterol regulatory element-binding proteins are essential for the metabolic programming of effector T cells and adaptive immunity. Nat Immunol (2013) 14:489–99. doi: 10.1038/ni.2570
28
KoppenolWHBoundsPLDangCV. Otto Warburg's contributions to current concepts of cancer metabolism. Nat Rev Cancer (2011) 11:325–37. doi: 10.1038/nrc3038
29
DonnellyRPLoftusRMKeatingSELiouKTBironCAGardinerCMet al. mTORC1-dependent metabolic reprogramming is a prerequisite for NK cell effector function. J Immunol (2014) 193:4477–84. doi: 10.4049/jimmunol.1401558
30
MichalekRDGerrietsVAJacobsSRMacintyreANMacIverNJMasonEFet al. Cutting edge: distinct glycolytic and lipid oxidative metabolic programs are essential for effector and regulatory CD4+ T cell subsets. J Immunol (2011) 186:3299–303. doi: 10.4049/jimmunol.1003613
31
DoughtyCABleimanBFWagnerDJDufortFJMatarazaJMRobertsMFet al. Antigen receptor-mediated changes in glucose metabolism in b lymphocytes: Role of phosphatidylinositol 3-kinase signaling in the glycolytic control of growth. Blood (2006) 107:4458–65. doi: 10.1182/blood-2005-12-4788
32
GerrietsVAKishtonRJNicholsAGMacintyreANInoueMIlkayevaOet al. Metabolic programming and PDHK1 control CD4+ T cell subsets and inflammation. J Clin Invest (2015) 125:194–207. doi: 10.1172/JCI76012
33
LeandroJHoutenSM. The lysine degradation pathway: Subcellular compartmentalization and enzyme deficiencies. Mol Genet Metab (2020) 131:14–22. doi: 10.1016/j.ymgme.2020.07.010
34
AbbasovMEKavanaghMEIchuTALazearMRTaoYCrowleyVMet al. A proteome-wide atlas of lysine-reactive chemistry. Nat Chem (2021) 13:1081–92. doi: 10.1038/s41557-021-00765-4
35
SinclairLVRolfJEmslieEShiY-BTaylorPMCantrellDAet al. Control of amino-acid transport by antigen receptors coordinates the metabolic reprogramming essential for T cell differentiation. Nat Immunol (2013) 14:500–8. doi: 10.1038/ni.2556
36
CarrELKelmanAWuGSGopaulRSenkevitchEAghvanyanAet al. Glutamine uptake and metabolism are coordinately regulated by ERK/MAPK during T lymphocyte activation. J Immunol (2010) 185:1037–44. doi: 10.4049/jimmunol.0903586
37
MyersKCSauterSZhangXBleesingJJDaviesSMWellsSIet al. Impaired immune function in children and adults with fanconi anemia. Pediatr Blood Cancer (2017) 64(11):10.1002/pbc.26599. doi: 10.1002/pbc.26599
38
KeppelMPSaucierNMahAYVogelTPCooperMA. Activation-specific metabolic requirements for NK cell IFN-gamma production. J Immunol (2015) 194:1954–62. doi: 10.4049/jimmunol.1402099
39
BuckMDO'SullivanDKlein GeltinkRICurtisJDChangCHSaninDEet al. Mitochondrial dynamics controls T cell fate through metabolic programming. Cell (2016) 166:63–76. doi: 10.1016/j.cell.2016.05.035
40
Mah-SomAYKeppelMPTobinJMKolicheskiASaucierNSexlVet al. Reliance on Cox10 and oxidative metabolism for antigen-specific NK cell expansion. Cell Rep (2021) 35:109209. doi: 10.1016/j.celrep.2021.109209
41
OhtaHKimuraIKonishiMItohN. Neudesin as a unique secreted protein with multi-functional roles in neural functions, energy metabolism, and tumorigenesis. Front Mol Biosci (2015) 2:24. doi: 10.3389/fmolb.2015.00024
42
SuXWangYLiAZanLWangH. Neudesin neurotrophic factor promotes bovine preadipocyte differentiation and inhibits myoblast myogenesis. Anim (Basel) (2019) 9(12):1109. doi: 10.3390/ani9121109
43
KimuraIKonishiMAsakiTFurukawaNUkaiKMoriMet al. Neudesin, an extracellular heme-binding protein, suppresses adipogenesis in 3T3-L1 cells via the MAPK cascade. Biochem Biophys Res Commun (2009) 381:75–80. doi: 10.1016/j.bbrc.2009.02.011
44
NorrisFAAtkinsRCMajerusPW. The cDNA cloning and characterization of inositol polyphosphate 4-phosphatase type II. evidence for conserved alternative splicing in the 4-phosphatase family. J Biol Chem (1997) 272:23859–64. doi: 10.1074/jbc.272.38.23859
45
DzneladzeIHeRWoolleyJFSonMHSharobimMHGreenbergSAet al. INPP4B overexpression is associated with poor clinical outcome and therapy resistance in acute myeloid leukemia. Leukemia (2015) 29:1485–95. doi: 10.1038/leu.2015.51
46
PenningTM. The aldo-keto reductases (AKRs): Overview. Chem Biol Interact (2015) 234:236–46. doi: 10.1016/j.cbi.2014.09.024
47
HashimotoHIshinoYJiangWYoshimuraTTakeda-UchimuraYUchimuraKet al. Keratan sulfate regulates the switch from motor neuron to oligodendrocyte generation during development of the mouse spinal cord. Neurochem Res (2016) 41:450–62. doi: 10.1007/s11064-016-1861-9
48
BrodskyRAJonesRJ. Riddle: What do aplastic anemia, acute promyelocytic leukemia, and chronic myeloid leukemia have in common? Leukemia (2004) 18:1740–2. doi: 10.1038/sj.leu.2403487
49
LimSPhillipsJBMadeira da SilvaLZhouMFodstadOOwenLBet al. Interplay between immune checkpoint proteins and cellular metabolism. Cancer Res (2017) 77:1245–9. doi: 10.1158/0008-5472.CAN-16-1647
Summary
Keywords
aplastic anemia (AA), cell metabolism, scRNA seq, T lymphocytes, NK cells, NENF
Citation
Zhou Q, Huang L, Liu Y, Huang J, Wen L, Yang J, Liang J, Chen Y and Chen C (2023) Single-cell RNA sequencing depicts metabolic changes in children with aplastic anemia. Front. Oncol. 13:1075408. doi: 10.3389/fonc.2023.1075408
Received
20 October 2022
Accepted
27 February 2023
Published
29 March 2023
Volume
13 - 2023
Edited by
Jian Yu, Beihang University, China
Reviewed by
Jian Ge, First Affiliated Hospital of Anhui Medical University, China; Hubing Shi, Sichuan University, China
Updates
Copyright
© 2023 Zhou, Huang, Liu, Huang, Wen, Yang, Liang, Chen and Chen.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Chun Chen, chenchun@mail.sysu.edu.cn; Yun Chen, cheny653@mail.sysu.edu.cn; Jintang Liang, ktleung@cuhk.edu.hk
†These authors have contributed equally to this work
This article was submitted to Cancer Metabolism, a section of the journal Frontiers in Oncology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.