Skip to main content

ORIGINAL RESEARCH article

Front. Oncol., 29 March 2023
Sec. Cancer Metabolism
This article is part of the Research Topic Linking Cellular Metabolism to Hematological Malignancies, volume II View all 11 articles

Single-cell RNA sequencing depicts metabolic changes in children with aplastic anemia

Qin Zhou&#x;Qin ZhouLifen Huang&#x;Lifen HuangYong LiuYong LiuJunbin HuangJunbin HuangLuping WenLuping WenJing YangJing YangJintang Liang*Jintang Liang*Yun Chen*Yun Chen*Chun Chen*Chun Chen*
  • Pediatric Hematology Laboratory, Division of Hematology/Oncology, Department of Pediatrics, The Seventh Affiliated Hospital of Sun Yat-Sen University, Shenzhen, Guangdong, China

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 (TNFSF12) 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
www.frontiersin.org

Figure 1 Identify cells and classify subsets with single-cell analysis. Schematic representation of the experimental design and single-cell mRNA sequencing (scRNA-seq) sequencing procedure (A). The surface marker of immune cells was used to identify cells and classify subsets. Dot sizes represent percent expressed and dot colors represent the average-expressed scale (B). Fraction of cell subsets’ differences between groups. Colors indicate cell clusters with numbered labels (C).

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
www.frontiersin.org

Figure 2 Single cell type identification and enrichment analysis. Uniform manifold approximation and project visualization of the control group (A) and aplastic anemia (AA) (B) based on single-cell transcriptomes (left). Gene Ontology (GO) enrichment classification histogram depicting the number of up- or downregulated genes in the biological process, cell component, and molecular function ontologies (C). Dot plot representing top 20 enriched GO terms (ranked by Q values) based on bulk RNA-seq (D).

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
www.frontiersin.org

Figure 3 scRNA-seq revealed metabolic differences in AA. t-SNE plot of all cells representing the cell clusters analyzed by scRNA-seq. Each dot represents a single cell; each color corresponds to one cluster (A). t-SNE plot showing the metabolism differences of lysine degradation between AA patients and healthy donors. Colors indicate logarithmic-transformed P values (B). t-SNE plot showing the significantly expressed gene NENF in AA patients compared with those in healthy donors (C).

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
www.frontiersin.org

Figure 4 scRNA-seq revealed metabolic differences in AA. Bar chart showing top 20 of the KEGG pathway enriched in cell metabolism compared to those in the control group (A). Gene set enrichment analysis (GSEA) plots of valine, leucine, and isoleucine degradation differences in AA patients compared with those in healthy donors (B). Heat map of differential genetic expression handled with the z-score to normalize gene expression (C). Bar chart showing the number of up- or downregulated genes in bulk and 19 subsets. Red indicates upregulation, and green indicates downregulation (D).

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
www.frontiersin.org

Figure 5 T-lymphocyte metabolism analysis. Dot plot depicting the top 20 KEGG pathway in T lymphocytes. Dot sizes represent the gene number, and dot colors represent Q values (A). t-SNE plot showing the differently expressed gene INPP4B in T lymphocytes. All but clusters 5 and 18 are lymphocytes. Each dot represents a single cell. Colors indicate logarithmic-transformed P values (B). Dot plot showing top 20 of disease ontology enrichment in T lymphocytes (C).

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 TNFSF12 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
www.frontiersin.org

Figure 6 Natural killer (NK) cell metabolism analysis. Bar chart showing top 20 of the KEGG pathway in NK cells. Column length represents the percentage of differential genes, and column colors represent Q values (A). t-SNE plot showing differently expressed genes AKR1C3 (B) and CHST2 (C) in NK cells. Cluster 5 is an NK cell. Each dot represents a single cell; Colors indicate logarithmic transformed P values. Bar chart depicting top 20 of disease ontology enrichment in NK cells (D). Intercellular interaction network between 18 subsets of control (E) and AA (F). Molecular interaction states of 36 ligand–receptors between clusters 4 and 5 in the control (G) and AA (H) groups. Molecules in red indicate that the Pearson correlation coefficient target gene prediction ability is higher.

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.

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. Young NS, Calado RT, Scheinberg P. Current concepts in the pathophysiology and treatment of aplastic anemia. Blood (2006) 108:2509–19. doi: 10.1182/blood-2006-03-010777

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Killick SB, Bown N, Cavenagh J, Dokal I, Foukaneli T, Hill A, et al. Guidelines for the diagnosis and management of adult aplastic anaemia. Br J Haematol (2016) 172:187–207. doi: 10.1111/bjh.13853

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Medinger M, Drexler B, Lengerke C, Passweg J. Pathogenesis of acquired aplastic anemia and the role of the bone marrow microenvironment. Front Oncol (2018) 8:587. doi: 10.3389/fonc.2018.00587

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Imi T, Katagiri T, Hosomichi K, Zaimoku Y, Nguyen VH, Nakagawa N, et 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

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Nissen C, Stern M. Acquired immune mediated aplastic anemia: Is it antineoplastic? Autoimmun Rev (2009) 9:11–6. doi: 10.1016/j.autrev.2009.02.032

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Kurre P, Johnson FL, Deeg HJ. Diagnosis and treatment of children with aplastic anemia. Pediatr Blood Cancer (2005) 45:770–80. doi: 10.1002/pbc.20322

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Montane E, Ibáñez L, Vidal X, Ballarín E, Puig R, García N, et al. Epidemiology of aplastic anemia: A prospective multicenter study. Haematologica (2008) 93:518–23. doi: 10.3324/haematol.12020

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Young NS, Maciejewski J. The pathophysiology of acquired aplastic anemia. N Engl J Med (1997) 336:1365–72. doi: 10.1056/NEJM199705083361906

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Young NS. Hematopoietic cell destruction by immune mechanisms in acquired aplastic anemia. Semin Hematol (2000) 37:3–14. doi: 10.1016/s0037-1963(00)90026-x

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Bacigalupo A. How I treat acquired aplastic anemia. Blood (2017) 129:1428–36. doi: 10.1182/blood-2016-08-693481

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Desmond R, Townsley DM, Dumitriu B, Olnes MJ, Scheinberg P, Bevans M, et 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

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Li X, Sun X, Carmeliet P. Hallmarks of endothelial cell metabolism in health and disease. Cell Metab (2019) 30:414–33. doi: 10.1016/j.cmet.2019.08.011

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Torrino S, Grasset EM, Audebert S, Belhadj I, Lacoux C, Haynes M, et 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

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Vander Heiden MG, Cantley LC, Thompson CB. Understanding the warburg effect: The metabolic requirements of cell proliferation. Science (2009) 324:1029–33. doi: 10.1126/science.1160809

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Andrejeva G, Rathmell JC. 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

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Nakaya M, Xiao Y, Zhou X, Chang J-H, Chang M, Cheng X, et 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

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Chang CH, Curtis JD, Maggi LB Jr, Faubert B, Villarino AV, O'Sullivan D, et al. Posttranscriptional control of T cell effector function by aerobic glycolysis. Cell (2013) 153:1239–51. doi: 10.1016/j.cell.2013.05.016

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Shi LZ, Wang R, Huang G, Vogel P, Neale G, Green DR, et 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

PubMed Abstract | CrossRef Full Text | Google Scholar

19. O'Sullivan D, Pearce EL. Targeting T cell metabolism for therapy. Trends Immunol (2015) 36:71–80. doi: 10.1016/j.it.2014.12.004

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Calvaresi EC, Hergenrother PJ. Glucose conjugation for the specific targeting and treatment of cancer. Chem Sci (2013) 4:2319–33. doi: 10.1039/C3SC22205E

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Fastres A, Pirottin D, Fievez L, Marichal T, Desmet CJ, Bureau F, et 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

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Gao S, Wu Z, Arnold B, Diamond C, Batchu S, Giudice V, et 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

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et 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

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Chen X, Yang TT, Qiu XC, Ji ZG, Li CX, Long H, et 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

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Zhu J, Yamane H, Paul WE. Differentiation of effector CD4 T cell populations (*). Annu Rev Immunol (2010) 28:445–89. doi: 10.1146/annurev-immunol-030409-101212

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Siska PJ, Rathmell JC. T Cell metabolic fitness in antitumor immunity. Trends Immunol (2015) 36:257–64. doi: 10.1016/j.it.2015.02.007

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Kidani Y, Elsaesser H, Hock MB, Vergnes L, Williams KJ, Argus JP, et 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

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Koppenol WH, Bounds PL, Dang CV. Otto Warburg's contributions to current concepts of cancer metabolism. Nat Rev Cancer (2011) 11:325–37. doi: 10.1038/nrc3038

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Donnelly RP, Loftus RM, Keating SE, Liou KT, Biron CA, Gardiner CM, et al. mTORC1-dependent metabolic reprogramming is a prerequisite for NK cell effector function. J Immunol (2014) 193:4477–84. doi: 10.4049/jimmunol.1401558

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Michalek RD, Gerriets VA, Jacobs SR, Macintyre AN, MacIver NJ, Mason EF, et 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

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Doughty CA, Bleiman BF, Wagner DJ, Dufort FJ, Mataraza JM, Roberts MF, et 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

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Gerriets VA, Kishton RJ, Nichols AG, Macintyre AN, Inoue M, Ilkayeva O, et al. Metabolic programming and PDHK1 control CD4+ T cell subsets and inflammation. J Clin Invest (2015) 125:194–207. doi: 10.1172/JCI76012

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Leandro J, Houten SM. The lysine degradation pathway: Subcellular compartmentalization and enzyme deficiencies. Mol Genet Metab (2020) 131:14–22. doi: 10.1016/j.ymgme.2020.07.010

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Abbasov ME, Kavanagh ME, Ichu TA, Lazear MR, Tao Y, Crowley VM, et al. A proteome-wide atlas of lysine-reactive chemistry. Nat Chem (2021) 13:1081–92. doi: 10.1038/s41557-021-00765-4

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Sinclair LV, Rolf J, Emslie E, Shi Y-B, Taylor PM, Cantrell DA, et 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

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Carr EL, Kelman A, Wu GS, Gopaul R, Senkevitch E, Aghvanyan A, et 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

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Myers KC, Sauter S, Zhang X, Bleesing JJ, Davies SM, Wells SI, et 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

CrossRef Full Text | Google Scholar

38. Keppel MP, Saucier N, Mah AY, Vogel TP, Cooper MA. Activation-specific metabolic requirements for NK cell IFN-gamma production. J Immunol (2015) 194:1954–62. doi: 10.4049/jimmunol.1402099

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Buck MD, O'Sullivan D, Klein Geltink RI, Curtis JD, Chang CH, Sanin DE, et al. Mitochondrial dynamics controls T cell fate through metabolic programming. Cell (2016) 166:63–76. doi: 10.1016/j.cell.2016.05.035

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Mah-Som AY, Keppel MP, Tobin JM, Kolicheski A, Saucier N, Sexl V, et 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

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Ohta H, Kimura I, Konishi M, Itoh N. 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

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Su X, Wang Y, Li A, Zan L, Wang H. Neudesin neurotrophic factor promotes bovine preadipocyte differentiation and inhibits myoblast myogenesis. Anim (Basel) (2019) 9(12):1109. doi: 10.3390/ani9121109

CrossRef Full Text | Google Scholar

43. Kimura I, Konishi M, Asaki T, Furukawa N, Ukai K, Mori M, et 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

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Norris FA, Atkins RC, Majerus PW. 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

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Dzneladze I, He R, Woolley JF, Son MH, Sharobim MH, Greenberg SA, et 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

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Penning TM. The aldo-keto reductases (AKRs): Overview. Chem Biol Interact (2015) 234:236–46. doi: 10.1016/j.cbi.2014.09.024

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Hashimoto H, Ishino Y, Jiang W, Yoshimura T, Takeda-Uchimura Y, Uchimura K, et 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

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Brodsky RA, Jones RJ. 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

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Lim S, Phillips JB, Madeira da Silva L, Zhou M, Fodstad O, Owen LB, et al. Interplay between immune checkpoint proteins and cellular metabolism. Cancer Res (2017) 77:1245–9. doi: 10.1158/0008-5472.CAN-16-1647

PubMed Abstract | CrossRef Full Text | Google Scholar

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.

Edited by:

Jian Yu, Beihang University, China

Reviewed by:

Jian Ge, First Affiliated Hospital of Anhui Medical University, China
Hubing Shi, Sichuan University, China

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, Y2hlbmNodW5AbWFpbC5zeXN1LmVkdS5jbg==; Yun Chen, Y2hlbnk2NTNAbWFpbC5zeXN1LmVkdS5jbg==; Jintang Liang, a3RsZXVuZ0BjdWhrLmVkdS5oaw==

These authors have contributed equally to this work

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.