论著 · 基础研究

成人与儿童急性髓系白血病患者肿瘤免疫相关的差异表达基因分析

  • 李玲玲 ,
  • 李倩 ,
  • 李明玉 ,
  • 刘峥 ,
  • 沈倩诚
展开
  • 1.上海市儿童医院,上海交通大学附属儿童医院中心实验室,上海 200062
    2.上海交通大学基础医学院医药生物信息学中心,上海 200025
    3.上海工程技术大学管理学院,上海 201620
李玲玲(1988—),女,实验技术员,硕士生;电子信箱:lingling117@sjtu.edu.cn

网络出版日期: 2021-05-27

基金资助

上海市卫生和计划生育委员会青年项目(20184Y0220)

Analysis of tumor immune-related differentially expressed genes in adults and children with acute myeloid leukemia

  • Ling-ling LI ,
  • Qian LI ,
  • Ming-yu LI ,
  • Zheng LIU ,
  • Qian-cheng SHEN
Expand
  • 1.Department of Central Laboratory, Shanghai Children's Hospital, Shanghai Jiao Tong University, Shanghai 200062, China
    2.Medicinal Bioinformatics Center, Shanghai Jiao Tong University College of Basic Medical Sciences, Shanghai 200025, China
    3.School of Management, Shanghai University of Engineering Science, Shanghai 201620, China

Online published: 2021-05-27

Supported by

Youth Project of Shanghai Municipal Commission of Health and Family Planning(20184Y0220)

摘要

目的·分析成人与儿童急性髓系白血病(acute myeloid leukemia,AML)患者骨髓中肿瘤免疫相关的差异表达基因(differentially expressed genes,DEGs)。方法·从GEO(Gene Expression Omnibus)数据库下载GSE134589数据集,选取初次确诊/复发状态的患者,按年龄分为儿童组(0~16岁,34例)、中青年组(17~59岁,62例)和老年组(60~80岁,62例),用R语言程序包筛选不同组患者骨髓样本中肿瘤免疫相关的DEGs。选取中青年组与儿童组,老年组与儿童组共同的DEGs,与完全缓解状态的成人与儿童患者的DEGs进行比对,并进行功能富集分析。利用Kaplan-Meier法筛选与预后显著相关的基因,通过构建蛋白质相互作用(protein-protein interaction,PPI)网络筛选核心调控基因,将以上2种方法筛选出的基因视为关键基因。用GEPIA服务器对比关键基因在AML、弥漫大B细胞淋巴瘤(diffuse large B cell lymphoma,DLBCL)和胸腺癌的肿瘤样本与正常人样本的表达量,在GEXC网站分析关键基因在AML患者肿瘤干细胞和健康人的原始造血细胞中mRNA的表达情况。结果·GSE134589数据集中,中青年组与儿童组比,上调的DEGs有51个,下调的有21个;老年组与儿童组比,上调的DEGs有47个,下调的有20个;而中青年组与老年组比,没有发现DEG。筛选了中青年组和老年组共同的肿瘤免疫相关DEGs 57个,其中上调有39个,下调有18个,仅有3个基因的表达水平差异在疾病完全缓解时仍有统计学意义。57个共同DEGs主要富集在白细胞迁移和细胞因子介导的信号通路,其中白细胞介素2受体α亚基(interleukin-2 receptor subunit alpha,IL2RA)、FMS-样酪氨酸激酶3(FMS-like tyrosine kinase 3,FLT3)高表达的患者与低表达的患者相比总体生存期显著缩短(均P<0.05),补体成分3a受体1(complement component 3a receptor 1,C3AR1)是PPI网络的核心调控基因。这3个基因作为关键基因,均在AML肿瘤样本特异性高表达(均P<0.05),IL2RA还在DLBCL患者样本中显著高表达(P<0.05)。IL2RA在AML肿瘤干细胞和健康人的原始造血细胞中均低表达,FLT3均高表达,而C3AR1的表达在AML肿瘤干细胞特异性升高。结论·成人与儿童AML患者预后的差异可能与骨髓中肿瘤免疫相关基因表达的差异有关,其中IL2RA、FLT3C3AR1可能是发挥重要作用的关键基因。

本文引用格式

李玲玲 , 李倩 , 李明玉 , 刘峥 , 沈倩诚 . 成人与儿童急性髓系白血病患者肿瘤免疫相关的差异表达基因分析[J]. 上海交通大学学报(医学版), 2021 , 41(5) : 579 -587 . DOI: 10.3969/j.issn.1674-8115.2021.05.004

Abstract

Objective

·To analyze tumor immune-related differentially expressed genes (DEGs) in the bone marrow of adults and children with acute myeloid leukemia (AML).

Methods

·The GSE134589 data set was downloaded from the GEO (Gene Expression Omnibus) database, and the patients in initial diagnosis/relapse were selected and divided into children group (0?16 years old, 34 cases), young and middle-aged group (17?59 years old, 62 cases) and elders group (60?80 years old, 62 cases). The tumor immune-related DEGs in the bone marrow samples of different groups of patients were screened by R language packages. The common DEGs of young and middle-aged group vs children group and elders group vs children group were compared with the DEGs of adults group vs children group in complete remission, and the functional enrichment analysis was conducted with the common DEGs. The Kaplan-Meier method was used to screen the DEGs that were significantly related to prognosis, and the core regulatory DEGs were screened by constructing a protein-protein interaction (PPI) network. The genes screened by the above two methods were regarded as key genes. The expression levels of the key genes in AML, diffuse large B cell lymphoma (DLBCL) and thymoma tumor samples and normal human samples were compared by GEPIA server, and the mRNA expression levels of key genes in the human AML tumor stem cells and the original hematopoietic cells of healthy people were analyzed by GEXC website.

Results

·In the GSE134589 data set, there were 51 DEGs up-regulated and 21 down-regulated between the young and middle-aged group and the children group; 47 DEGs were up-regulated and 20 down-regulated between the elders group and the children group; no DEG was found between the young and middle-aged group and the elders group. Fifty-seven tumor immune-related DEGs were screened in both the young and middle-aged group and the elders group, of which 39 were up-regulated and 18 were down-regulated; in these DEGs only 3 genes showed statistically significant difference in expression level when the disease was incomplete remission. The 57 common DEGs were mainly enriched in leukocyte migration and cytokine-related signal pathways. The patients with high expression of interleukin-2 receptor subunit alpha (IL2RA) and FMS-like tyrosine kinase 3 (FLT3) had significantly shorter overall survival than those with low expression (both P<0.05), and complement component 3a receptor 1 (C3AR1) was the core regulatory gene of the PPI network. The three DEGs were selected as key genes. They were all specifically highly expressed in the AML tumor samples (all P<0.05), and IL2RA was also significantly highly expressed in the DLBCL samples (P<0.05). The expression level of IL2RA was low both in the AML tumor stem cells and the original group of hematopoietic cells, but FLT3 was highly expressed in these cells. The expression level of C3AR1 was specifically high in the AML tumor stem cells.

Conclusion

·The difference in the prognosis between adults and children with AML may be related to the differences in the expression of tumor immune-related genes in bone marrow, in which IL2RA, FLT3 and C3AR1 may be the key genes.

参考文献

1 Vago L, Gojo I. Immune escape and immunotherapy of acute myeloid leukemia[J]. J Clin Invest, 2020, 130(4): 1552-1564.
2 D?hner H, Weisdorf DJ, Bloomfield CD. Acute myeloid leukemia[J]. N Engl J Med, 2015, 373(12): 1136-1152.
3 Swartz MA, Iida N, Roberts EW, et al. Tumor microenvironment complexity: emerging roles in cancer therapy[J]. Cancer Res, 2012, 72(10): 2473-2480.
4 Edgar R, Domrachev M, Lash AE. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository[J]. Nucleic Acids Res, 2002, 30(1): 207-210.
5 Vadakekolathu J, Minden MD, Hood T, et al. Immune landscapes predict chemotherapy resistance and immunotherapy response in acute myeloid leukemia[J]. Sci Transl Med, 2020, 12(546): eaaz0463.
6 Raney BJ, Dreszer TR, Barber GP, et al. Track data hubs enable visualization of user-defined genome-wide annotations on the UCSC Genome Browser[J]. Bioinformatics, 2014, 30(7): 1003-1005.
7 Ritchie ME, Phipson B, Wu D, et al. Limma powers differential expression analyses for RNA-sequencing and microarray studies[J]. Nucleic Acids Res, 2015, 43(7): e47.
8 Zhou Y, Zhou B, Pache L, et al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets[J]. Nat Commun, 2019, 10(1): 1523.
9 Tang Z, Li C, Kang B, et al. GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses[J]. Nucleic Acids Res, 2017, 45(W1): W98-W102.
10 Szklarczyk D, Franceschini A, Wyder S, et al. STRING v10: protein-protein interaction networks, integrated over the tree of life[J]. Nucleic Acids Res, 2015, 43(Database issue): D447-D452.
11 Shannon P, Markiel A, Ozier O, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks[J]. Genome Res, 2003, 13(11): 2498-2504.
12 Chin CH, Chen SH, Wu HH, et al. cytoHubba: identifying hub objects and sub-networks from complex interactome[J]. BMC Syst Biol, 2014, 8(): S11.
13 Bandettini WP, Kellman P, Mancini C, et al. MultiContrast Delayed Enhancement (MCODE) improves detection of subendocardial myocardial infarction by late gadolinium enhancement cardiovascular magnetic resonance: a clinical validation study[J]. J Cardiovasc Magn Reson, 2012, 14: 83.
14 Gentles AJ, Plevritis SK, Majeti R, et al. Association of a leukemic stem cell gene expression signature with clinical outcomes in acute myeloid leukemia[J]. JAMA, 2010, 304(24): 2706-2715.
15 Klemm F, Joyce JA. Microenvironmental regulation of therapeutic response in cancer[J]. Trends Cell Biol, 2015, 25(4): 198-213.
16 Junttila MR, de Sauvage FJ. Influence of tumour micro-environment heterogeneity on therapeutic response[J]. Nature, 2013, 501(7467): 346-354.
17 Bennett JM, Catovsky D, Daniel MT, et al. Proposals for the classification of the acute leukaemias[J]. Br J Haematol, 1976, 33(4): 451-458.
18 Olkhanud PB, Damdinsuren B, Bodogai M, et al. Tumor-evoked regulatory B cells promote breast cancer metastasis by converting resting CD4? T cells to T-regulatory cells[J]. Cancer Res, 2011, 71(10): 3505-3515.
19 Wejksza K, Lee-Chang C, Bodogai M, et al. Cancer-produced metabolites of 5-lipoxygenase induce tumor-evoked regulatory B cells via peroxisome proliferator-activated receptor α[J]. J Immunol, 2013, 190(6): 2575-2584.
20 Yamamoto Y, Kiyoi H, Nakano Y, et al. Activating mutation of D835 within the activation loop of FLT3 in human hematologic malignancies[J]. Blood, 2001, 97(8): 2434-2439.
21 Reca R, Mastellos D, Majka M, et al. Functional receptor for C3a anaphylatoxin is expressed by normal hematopoietic stem/progenitor cells, and C3a enhances their homing-related responses to SDF-1[J]. Blood, 2003, 101(10): 3784-3793.
22 Mo?hle R, Bautz F, Rafii S, et al. The chemokine receptor CXCR-4 is expressed on CD34+ hematopoietic progenitors and leukemic cells and mediates transendothelial migration induced by stromal cell-derived factor-1[J]. Blood, 1998, 91(12): 4523-4530.
23 Lin Y, Ma Q, Li L, et al. The CXCL12-CXCR4 axis promotes migration, invasiveness, and EMT in human papillary thyroid carcinoma B-CPAP cells via NF-κB signaling[J]. Biochem Cell Biol, 2018, 96(5): 619-626.
24 Martin GH, Roy N, Chakraborty S, et al. CD97 is a critical regulator of acute myeloid leukemia stem cell function[J]. J Exp Med, 2019, 216(10): 2362-2377.
25 Zhao K, Yin LL, Zhao DM, et al. IL1RAP as a surface marker for leukemia stem cells is related to clinical phase of chronic myeloid leukemia patients[J]. Int J Clin Exp Med, 2014, 7(12): 4787-4798.
26 Warren CFA, Wong-Brown MW, Bowden NA. BCL-2 family isoforms in apoptosis and cancer[J]. Cell Death Dis, 2019, 10(3): 177.
27 Ferlin M, Noraz N, Hertogh C, et al. Insulin-like growth factor induces the survival and proliferation of myeloma cells through an interleukin-6-independent transduction pathway[J]. Br J Haematol, 2000, 111(2): 626-634.
28 Di Pietro A, Conseil G, Pérez-Victoria JM, et al. Modulation by flavonoids of cell multidrug resistance mediated by P-glycoprotein and related ABC transporters[J]. Cell Mol Life Sci, 2002, 59(2): 307-322.
29 Chiu HY, Sun KH, Chen SY, et al. Autocrine CCL2 promotes cell migration and invasion via PKC activation and tyrosine phosphorylation of paxillin in bladder cancer cells[J]. Cytokine, 2012, 59(2): 423-432.
文章导航

/