Comprehensive Expression and Prognosis Analyses of ITGB1 in AML | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research article Comprehensive Expression and Prognosis Analyses of ITGB1 in AML Xinyi Zhou, Nan Jin, Baoan Chen This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-88926/v1 This work is licensed under a CC BY 4.0 License Status: Posted Version 1 posted You are reading this latest preprint version Abstract Background Acute myeloid leukemia (AML) is a dangerous type of leukemia. The emergence of multidrug resistance (MDR) and recurrence limits the prognosis and survival of patients. In recent studies, we have known that the bone marrow microenvironment was closely related to the poor prognosis of AML. However, the underlying mechanisms are still far from fully understood. By utilizing the bioinformatics analysis, we screened out integrin β1 (ITGB1) as the hub gene, which is associated with the bone marrow microenvironment mediated changes of AML cells, with expression profile GSE73157 downloaded from National Center for Biotechnology Information-Gene Expression Omnibus (NCBI-GEO) database. Methods R studio software was used to screen out candidate hub genes and further visualize the differential expression. R package “limma” was to find out differentially expressed genes (DEGs). Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway were conducted by R package “cluster Profiler”. Furthermore, protein-protein interaction (PPI) network was also performed by online tool STRING and software Cytoscape. Last but not least, online tool PrognoScan and GEPIA was utilized for the evaluation of clinical significance of the selected hub gene. P and Cox p value <0.05 was considered to be statistical significance. Results ITGB1 was filtrated as the only hub gene in this profile. We found that patients with high expression of ITGB1 had significantly longer overall survival (OS) than those with low expression (COX p value= 0.016730). Besides, the expression of the ITGB1 gene in AML patients is lower than that in normal people significantly (p value<0.01). Conclusion We identified ITGB1 as a key gene in the bone marrow microenvironment mediated poor prognosis in AML. The down-regulated expression of ITGB1 was related to AML patients’ poor outcome. ITGB1 may be a potential marker for predicting and guiding AML treatment. Cancer Biology Oncology Acute myeloid leukemia the bone marrow microenvironment Bioinformatics analysis Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Background Leukemia is a malignant clonal disease of the hematopoietic system, which accounts for 4% of all cancers[1]. Among them, acute myeloid leukemia (AML) is a common and dangerous type of leukemia. Currently, chemotherapy is still the most important option for the treatment of AML. However, the emergence of multidrug resistance (MDR) and recurrence greatly limits the effects of chemotherapy drugs. Therefore, the research on MDR of AML cells has always been a research hotspot in hematology. In recent years, studies have found that the bone marrow microenvironment not only plays an important role in the hematopoietic process, but is closely related to the occurrence and development of tumors. The bone marrow microenvironment contains multiple cell types, such as osteoblasts, mesenchymal stem cells, endothelial cells, neurons and so on. In normal humans, the bone marrow microenvironment affects hematopoietic stem cells mostly through the regulation of blood oxygen partial pressure and the secretion of cytokines[2]. However, in cancer patients, the occurrence of MDR and recurrence is inseparable from the role of bone marrow microenvironment. Leukemia and leukemia stem cells (LSCs) interact with the bone marrow microenvironment in a different way from hematopoietic stem cells to affect the development of leukemia[3]. LSCs can hide in the bone marrow niche to escape the killing effect of chemotherapy. As a result, the bone marrow niche has also become a key part of the recurrence and MDR of leukemia[4-6]. According to recent researches, we have known that MDR mediated by the bone marrow microenvironment might stem from the following several mechanisms: 1. The influence of cytokines[7, 8]; 2. Changes in adhesion of leukemia cells [9, 10]; 3. The regulation of expression of MDR-related genes[11]; 4. The adjustment of cell cycle and metabolism[12]. In addition, many factors can also be contributed to the recurrence of leukemia, such as the secretion of cytokines or the changes with certain pathways. As the emergence and rapid development of high-throughput platforms and microarray, numerous molecular heterogeneity on AML has been acknowledged. In order to find potential key biomarker of AML related to the bone marrow microenvironment-mediated poor prognosis, through the use of biological information technology and the analysis of high-throughput sequencing chip data, we screen out the differentially expressed genes (DEGs) of AML cells cultured alone and co-cultured with stromal cells, which were all treated with different concentrations of Arsenic trioxide (ATO). Furthermore, we pick up ITGB1 as a hub gene that can be considered as the predictor of prognosis. We systematically revel the occurrence and development of AML on the level of molecule and provide a potential guidance for targeted therapy. Materials And Methods 2.1 Analysis and screening of differential expressed genes (DEGs) The GEO database ( https://www.ncbi.nlm.nih.gov/geo/ ) of the National Center for Biotechnology Information was used to select and download the gene express profile GSE73157[13]. This profile was generated on the platform of GPL17077 (Agilent-039494 SurePrint G3 Human GE v2 8*60K Microarray 039381). In order to reduce false positive results, the online analysis tool GEO2R in the GEO database was used to process the data profile. The screen criteria were adjusted p value < 0.05 and the use of Benjamini-Hochberg. Then, R package “limma” was used to clarity the DEGs between NB4 cultured alone and co-cultured with stromal cells after the treatment of ATO. The screen criteria were adjusted p value 1[14]. All genes were visualized by volcanic maps and top 50 dramatically differentially expressed genes were selected to draw a heatmap by R package “ggplot2”[15]. 2.2 Functional enrichment analysis and pathway analysis of DEGs R package “clusterProfiler” was used to conduct Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis in order to preform functional enrichment and pathway analysis of DEGs[16]. Biological process (BP), molecular function (MF) and cellular component (CC) are included in GO enrichment. Analysis results were illustrated as figures by R package “GOplot”[17]. A p value less than 0.05 was considered statistically significant. Also, the ClueGo plug-ins of Cytoscape software 3.8.1 was utilized to display the relationship between pathways[18]. A p value less than 0.05 was considered statistically significant. 2.3 Screening candidate genes through protein-protein interaction (PPI) network The PPI prediction website STRING ( https://string-db.org ) constructed protein interaction network[19]. And Cytoscape software 3.8.0 was used to visualize the interactions of DEGs[18]. Further, we used “cytoHubba” to degree the interaction of hub-gene clustering by the method of Maxima Clique Centrality (MCC). Overlapped genes were separated by online toll Bioinformatics & Evolutionary Genomics (http://bioinformatics.psb.ugent.be/webtools/Venn/) and were displayed by Wayne diagram. 2.4 Selecting candidate genes and conducting survival analysis In order to evaluate the prognostic value of candidate genes in AML patients, we utilized the PrognoScan database ( http://dna00.bio.kyutech.ac.jp/PrognoScan/ ) and online tool Gene Expression Profiling Interaction Analysis ( http://gepia.cancer-pku.cn/index.html ) to conduct the survival analysis[20, 21]. The results are displayed with hazard ratio (HR) and Cox p value from a Log-rank test. Cox p value < 0.05 was considered statistically significant. Besides, according to online tool Gene Expression Profiling Interaction Analysis (GEPIA), we compared the expression of candidate genes in AML patients and normal people. The screen criteria were adjusted p value 1 Results 3.1 Analysis and screening of differential expressed genes (DEGs) The goal of this study is to screen out the hub gene of AML related to the bone marrow microenvironment-mediated MDR (Figure 1). Perform statistical analysis on GEO database, we found that the data set GSE73157contains 32080 DEGs in AML cultured alone groups and co-cultured groups (All the following are referred to as alone groups and co-culture groups, respectively). By using R package “limma”, DEGs were classified as up-regulated genes (marked as red plots) and down-regulated genes (marked as green plots) (Figure 2a). Besides, a heatmap illustrated top 50 DEGs (Figure 2b). 3.2 Functional enrichment analysis and pathway analysis of DEGs Go and KEGG enrichment analysis was performed on all DEGs with R package “clusterProfiler”. The results of GO include three subontologies referred to as biological process (BP), cellular components (CC), molecular function (MF). The GO results show that angiogenesis pathway, cell adhesion molecule binding pathway and collagen-containing extracellular matrix pathway were selected as the most significant pathway in each subontologies, respectively (Figure 3a-c). And we chose angiogenesis pathway as our further analysis pathway as a result of comparation of p values. 177 genes were found enriched in this GO term. In addition, PI3K-Akt signaling pathway was the top enriched pathway of the DEGs in KEGG enrichment analysis (Figure 3d). After we used the ClueGo to further analyze the interrelation of the enriched pathways and the DEGs, PI3K-Akt signaling pathway was still the most significant pathway, and involved 120 DEGs (Figure 3e, f). Finally, we selected two pathways for our further research, which involved 177 and 120 DEGs, respectively (Table1). Then, after taking the intersection of the DEGs included in the two pathways, we could find that there are 26 genes included both in the two pathways (Figure 3g). They were PDGFRB, FLT1, NGFR, EFNA1, FGFR1, ITGB1, VEGFC, BRCA1, PDGFRA, PGF, EPHA2, EFNA3, PIK3CG, SYK, EREG, PTK2, ITGB8, FGF2, THBS1, FGF18, THBS2, AKT3, IL6, ANGPT1, FN1, ANGPT2. DEGs Genes Angiogenesis pathway ABL1, ACVR1, ACVRL1, ADM, AKT3, AMOTL2, ANG, ANGPTL4, ANGPTL6, ANXA2, AQP1, ARHGAP22, C3, CARD10, CAV1, CCBE1, CCL2, CD40, CDH13, CELA1, CHI3L1, CHRNA7, COL18A1, COL8A2, CYP1B1, DCN, DDAH1, DLL4, E2F7, ECM1, EFNA1, EFNA3, EGFL7, ELK3, EMP2, ENPP2, EPAS1, EPHA2, EPHB2, EREG, ESM1, FAP, FGF2, FLT1, FN1, FZD5, GATA6, GBX2, GPNMB, GPR4, GREM1, HDAC9, HEY1, HLA-G, HMOX1, HOXB13, HPSE, HTATIP2, ID1, IL18, IL1A, IL1B, IL6, ITGAX, ITGB1, ITGB2, ITGB8, KLF5, LEF1, LOXL2, LRG1, MAPK14, MDK, MFGE8, MMP14, MMP2, NFATC4, NGFR, NOTCH3, NRP1, NRP2, PARVA, PDGFRA, PDGFRB, PGF, PGK1, PLCD1, PLK2, PPARG, PRKD2, PTGS2, PTK2, PTN, PTPRM, RAMP3, RBM15, RHOB, RNF213, RORA, RRAS, SASH1, SAT, SERPINE1, SERPINF1, SFRP1, SFRP2, SMAD1, SOX17, SPARC, SRPX2, SYNJ2BP, TGFBI, THBS1, THBS2, THY1, TIE1, TNFAIP3, TNFSF12, TWIST1, VASH1, VEGFC, WNT5A, ADM2, AGGF1, ANGPT1, ANGPT2, BRCA1, CCR2, CX3CL1, CXCR2, CYBB, DAB2IP, DAG1, ECSCR, EFNB2, EMILIN1, EPHA1, F3, FBXW7, FGF18, FGFBP1, FGFR1, FOXJ2, GATA2, GPLD1, HHEX, HOXA7, IHH, JAG1, JUN, KRIT1, LEMD3, MAPK7, MED1, NCL, NOTCH4, NPPB, NRAR, PDE3B, PIK3CG, PKNOX1, POFUT1, PRKX, ROCK2, SETD2, SIRT1, SOX18, SPRED1, SPRY2, SRF, STAB1, SYK, TBX1, TERT, TGFBR1, THSD7A, UBP1 PI3K-Akt signaling pathway TCL1B, GNGT2, OSMR, GNG11, GNG10, IGF2, IGF1R, LAMB2, LAMB3, BCL2, GNG7, F2R,, LAMA1, COL9A3, MYB, LAMA4, EIF4EBP1, LAMB1, ITGA11, THBS1, PDGFRA, THBS2, TNC, GNG12, pik3ap1, INSR, FGF18, LPAR1, RXRA, MET, HRAS, LPAR5, HSP90B1, AREG, KIT, CCND2, PDGFC, FLT3, RPTOR, FLT1, PDGFRB, LPAR4, PIK3CG, PIK3R1, CCNE2, ITGA1, IKBKB, FN1, AKT3, SOS2, IRS1, CREB3, CDC37, KITLG, NRAS, IL3RA, Ifna4, CDK6, IL2RB, NGFR, CDK4, BRCA1, IL2RA, JAK2, THEM4, JAK3, EREG, ITGA2B, ITGA3, TLR4, ITGB4, MAPK3, ITGB5, CREB3L2, ITGB7, PIK3R5, IL6R, ANGPT1, TLR2, ANGPT2, EFNA1, GYS2, IL6, CDKN1A, TSC2ITGB1, PGF, YWHAB, PTK2, EFNA3, NFKB1, SPP1, SYK, VEGFC, CSF3, GNG2, FGF2, LAMC3, IFNB1 Table 1 : DEGs identified from selected pathways of GO and KEGG 3.3 Screening candidate genes through protein-protein interaction (PPI) network In order to further understand the DEGs, we displayed the protein-protein interaction (PPI) network including all the DEGs in the two selected pathways mentioned above in STRING (Figure 4a, b). Then, cytoHubba plug-ins in Cytoscape was also used to screen out top 30 candidate hub genes of each pathway according nodes rank (Figure 4c, d), and all the top hub genes has been listed in Table 2. In this part, we identified 4 common genes in the two sets of top 30 hub genes, including FN1, THBS1, ITGB1 and THBS2 as candidate hub genes. Angiogenesis pathway PI3K-Akt signaling pathway Rank Name Score Rank Name Score 1 FN1 1.68E+08 1 PTK2 7.42E+14 2 CCL2 1.67E+08 2 ITGB1 7.42E+14 3 MMP2 1.65E+08 3 ITGA3 7.42E+14 4 FGF2 1.58E+08 4 ITGA1 7.42E+14 5 SERPINE1 1.38E+08 5 ITGA2B 7.42E+14 6 THBS1 1.33E+08 6 ITGB5 7.42E+14 7 IL6 1.24E+08 6 ITGA11 7.42E+14 8 VEGFC 1.21E+08 8 ITGB7 7.42E+14 9 FLT1 1.07E+08 9 ITGB4 7.40E+14 10 PGF 1.07E+08 10 ITGB8 7.40E+14 11 ANGPT2 1.04E+08 10 ITGA10 7.40E+14 12 ANGPT1 1.00E+08 12 COL1A1 7.18E+14 13 PTGS2 6.38E+07 13 COL1A2 7.17E+14 14 COL18A1 4.98E+07 14 COL6A1 7.17E+14 15 IL1B 4.08E+07 15 COL6A3 7.15E+14 16 JUN 3.63E+07 16 COL6A2 7.14E+14 17 HMOX1 3.31E+07 17 LAMA4 3.83E+14 18 MMP14 2.35E+07 18 LAMB1 3.82E+14 19 PPARG 2.18E+07 19 THBS2 3.58E+14 20 CAV1 1.86E+07 20 COL9A3 3.56E+14 21 IL18 1.73E+07 21 LAMB2 2.37E+13 22 MAPK14 1.58E+07 22 LAMA1 2.36E+13 23 ITGB1 1.34E+07 23 LAMC3 2.10E+13 24 PDGFRB 1.30E+07 24 LAMB3 2.09E+13 25 SIRT1 1.10E+07 25 THBS1 3.95E+12 26 F3 7622520 26 SPP1 1.58E+12 27 ANG 7258344 27 FN1 1.43E+12 28 THBS2 6250378 28 TNC 9.39E+10 29 IL1A 5449080 29 PIK3CG 6.24E+09 30 CYBB 3774960 30 PIK3R1 6.24E+09 Table 2: The top 30 genes with the highest score of each pathway through the Cytoscape “cytoHubba” module analysis 3.4 Selecting candidate genes and conducting survival analysis Subsequently, the correlation between hub genes and prognosis of AML patients was analyzed through PrognoScan and GEPIA online survival website. The results showed that except for ITGB1 gene, as the result of no survival significance, the other genes were removed. We found that patients with high expression of ITGB1 had significantly longer overall survival (OS) than those with low expression (COX p value= 0.016730) (figure 5a, b). Besides, through the online website GEPIA, we performed a visual analysis of the expression of hub gene ITGB1 in the TCGA database in AML patients and normal people, showing that the expression of the ITGB1 gene in AML patients is lower than that in normal people significantly (p value<0.01) (Figure 5c). Discussion With the elucidation of molecular mechanisms and the emergence of targeted therapies, the treatment of AML has made great progress, and the overall survival rate of patients has been significantly improved. However, the resistance of leukemia cells and the recurrence after treatment are still two major problems that plague us. We now believe that the protective effect of leukemia stem cells (LSCs) mediated by the bone marrow microenvironment is also an important factor for the refractory and recurrence of leukemia. Recent studies have shown that there is a complex interaction between leukemia cells and the bone marrow microenvironment. The two induce and support each other, which together lead to the occurrence and development of the disease. In this study, we analyzed 32080 DEGs between cultured alone and co-cultured with stomal cells AML cells based on GSE73157 dataset. Then, we used GO enrichment analysis and KEGG pathway analysis to detect the dataset to explore the interaction between DEGs. GO analysis showed that these genes mainly reacted with the process of angiogenesis, collagen-containing extracellular matrix production and cytokine binding, suggesting that these DEGs may be involved in improving the ability of tumor cell to migrate and invade, promoting metastasis, recurrence and multiple drug resistance of cancer[22, 23]. KEGG analysis illustrated that 120 DEGs mainly affect the AML cells through PI3K-Akt signaling pathway. As we all known, phosphatidylinositol 3-kinase/protein kinase B/mammalian target of rapamycin (PI3K/Akt/mTOR) signaling pathway is over activated in hematological malignancies, and the development of inhibitors of PI3K/Akt/mTOR signaling pathway has become a common concern for hematology researchers in recent years[24]. Excessive activation of this pathway may promote the tumor proliferation and weaken the immune monitoring of tumor cells[25]. At the same time, we utilized the online tool STRING to display the functional interaction between proteins to discover the underlying mechanism of occurrence of development of AML. And then using the software Cytoscape 3.8.1 to analyze the protein interaction network diagram, we screened out FN1, THBS1, ITGB1, THBS2 as 4 candidate hub genes. The above 4 candidate genes not only interact closely with other genes, but may also determine the functions of other genes. After further survival analysis of the 4 selected genes, we found that only ITGB1 gene is statistically significant. We found that the expression of ITGB1 is negatively correlated with the prognosis of AML. Furthermore, compared with normal people, the expression of ITGB1 is significantly lower in AML patients. That is, ITGB1 may be a key gene for evaluating the prognosis and conducting targeted therapy of AML patients. A huge amount of evidences show that integrin family molecules are involved in the occurrence and development of tumors, and play an important role in regulating tumor cell proliferation, invasion and metastasis and other biological behaviors[26]. Integrin β1 (ITGB1) is a relatively important part of molecules in the cell adhesion molecule family, and is the main receptor that mediates the interaction between cells and extracellular matrix. This saying is consistent with our above analysis results: the bone marrow microenvironment may affect AML cells by regulating cell adhesion molecules and extracellular matrix. At present, our understanding of ITGB1 is limited to the primary stage, and its biological functions are not yet very clear. In the direction of solid tumors, there are certain studies on the direction of ITGB1, but there are few reports on the research of ITGB1 of hematological tumors in the world. In summary, ITGB1 may be a potential marker for predicting and guiding the treatment of AML, which requires further attention and research. Conclusion In summary, ITGB1 was a hub gene related to the bone microenvironment mediated poor prognosis in AML. The expression of ITGB1 was negatively correlated with prognosis of AML. This study used bioinformatics methods to analyze the relevant data of AML researches to obtain potential markers for predicting and guiding AML treatment, which provides a theoretical basis for a deeper understanding of the occurrence and development of AML, and also give the follow-up animal experiments and clinical trials research directions. Abbreviations AML: Acute myeloid leukemia MDR: Multidrug resistance ITGB1: Integrin β1 NCBI-GEO: National Center for Biotechnology Information-Gene Expression Omnibus database DEGs: Differentially expressed genes GO: Gene Ontology KEGG: Kyoto Encyclopedia of Genes and Genomes PPI: Protein-protein interaction OS: Overall survival LSCs: Leukemia stem cells ATO: Arsenic trioxide BP: Biological process MF: Molecular function CC: Cellular component MCC: Maxima Clique Centrality HR: Hazard ratio GEPIA: Gene Expression Profiling Interaction Analysis Declarations Ethics approval and consent to participate Not applicable. Consent for publication Not applicable. Availability of data and materials The dataset analysed during the current study are available in the NCBI-GEO repository, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE73157. Competing interests The authors report no conflicts of interest in this work. Funding This work was supported by The Natural Science Foundation of Jiangsu Province (BK20150639); Key Medical of Jiangsu Province (ZDXKB2016020); Six talent peaks project in Jiangsu Province (WSW-033);Innovative and entrepreneurial doctors of Jiangsu province. Authors' contributions XYZ, NJ and BAC were responsible for confirming the topic, collecting the data and analyzing article. XYZ designed and wrote the paper, edited the figures, legends and tables. All authors read and approved the final manuscript. Acknowledgements Not applicable. References Siegel, R.L., K.D. Miller, and A. Jemal, Cancer statistics, 2020. CA Cancer J Clin, 2020. 70 (1): p. 7-30. Kfoury, Y. and D.T. Scadden, Mesenchymal cell contributions to the stem cell niche. Cell Stem Cell, 2015. 16 (3): p. 239-53. Krause, D.S., et al., Differential regulation of myeloid leukemias by the bone marrow microenvironment. Nat Med, 2013. 19 (11): p. 1513-7. Frisch, B.J., et al., Functional inhibition of osteoblastic cells in an in vivo mouse model of myeloid leukemia. Blood, 2012. 119 (2): p. 540-50. Fulzele, K., et al., Myelopoiesis is regulated by osteocytes through Gsα-dependent signaling. Blood, 2013. 121 (6): p. 930-9. Welner, R.S., et al., Treatment of chronic myelogenous leukemia by blocking cytokine alterations found in normal stem and progenitor cells. Cancer Cell, 2015. 27 (5): p. 671-81. Flavell, R.A., et al., The polarization of immune cells in the tumour environment by TGFbeta. Nat Rev Immunol, 2010. 10 (8): p. 554-67. Massagué, J., TGFbeta in Cancer. Cell, 2008. 134 (2): p. 215-30. Sugiyama, T., et al., Maintenance of the hematopoietic stem cell pool by CXCL12-CXCR4 chemokine signaling in bone marrow stromal cell niches. Immunity, 2006. 25 (6): p. 977-88. Dillmann, F., et al., Plerixafor inhibits chemotaxis toward SDF-1 and CXCR4-mediated stroma contact in a dose-dependent manner resulting in increased susceptibility of BCR-ABL+ cell to Imatinib and Nilotinib. Leuk Lymphoma, 2009. 50 (10): p. 1676-86. Hazlehurst, L.A., T.H. Landowski, and W.S. Dalton, Role of the tumor microenvironment in mediating de novo resistance to drugs and physiological mediators of cell death. Oncogene, 2003. 22 (47): p. 7396-402. Colmone, A., et al., Leukemic cells create bone marrow niches that disrupt the behavior of normal hematopoietic progenitor cells. Science, 2008. 322 (5909): p. 1861-5. Barrett, T., et al., NCBI GEO: archive for functional genomics data sets--update. Nucleic Acids Res, 2013. 41 (Database issue): p. D991-5. Ritchie, M.E., et al., limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res, 2015. 43 (7): p. e47. Ito, K. and D. Murphy, Application of ggplot2 to Pharmacometric Graphics. CPT Pharmacometrics Syst Pharmacol, 2013. 2 (10): p. e79. Yu, G., et al., clusterProfiler: an R package for comparing biological themes among gene clusters. Omics, 2012. 16 (5): p. 284-7. Walter, W., F. Sánchez-Cabo, and M. Ricote, GOplot: an R package for visually combining expression data with functional analysis. Bioinformatics, 2015. 31 (17): p. 2912-4. Shannon, P., et al., Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res, 2003. 13 (11): p. 2498-504. Szklarczyk, D., et al., STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res, 2019. 47 (D1): p. D607-d613. Mizuno, H., et al., PrognoScan: a new database for meta-analysis of the prognostic value of genes. BMC Med Genomics, 2009. 2 : p. 18. Tang, Z., et al., GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res, 2017. 45 (W1): p. W98-w102. Li, T., et al., Tumor angiogenesis and anti-angiogenic gene therapy for cancer. Oncol Lett, 2018. 16 (1): p. 687-702. Terranova, V.P., D. Maslow, and G. Markus, Directed migration of murine and human tumor cells to collagenases and other proteases. Cancer Res, 1989. 49 (17): p. 4835-41. Gao, Y., C.Y. Yuan, and W. Yuan, Will targeting PI3K/Akt/mTOR signaling work in hematopoietic malignancies? Stem Cell Investig, 2016. 3 : p. 31. O'Donnell, J.S., et al., PI3K-AKT-mTOR inhibition in cancer immunotherapy, redux. Semin Cancer Biol, 2018. 48 : p. 91-103. Kapp, T.G., et al., Integrin modulators: a patent review. Expert Opin Ther Pat, 2013. 23 (10): p. 1273-95. Cite Share Download PDF Status: Posted Version 1 posted You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-88926","acceptedTermsAndConditions":true,"allowDirectSubmit":true,"archivedVersions":[],"articleType":"Research article","associatedPublications":[],"authors":[{"id":3219679,"identity":"80d89792-f7ab-48d0-91d3-f4e5789ef840","order_by":0,"name":"Xinyi Zhou","email":"","orcid":"","institution":"Southeast University Zhongda Hospital Department of Hematology","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Xinyi","middleName":"","lastName":"Zhou","suffix":""},{"id":3219680,"identity":"4e2b479d-b98a-4c5c-9828-030b459ef194","order_by":1,"name":"Nan Jin","email":"","orcid":"","institution":"Southeast University Zhongda Hospital Department of Hematology","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Nan","middleName":"","lastName":"Jin","suffix":""},{"id":3219681,"identity":"1469c18a-24d0-40b1-b57b-a30720c47075","order_by":2,"name":"Baoan Chen","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAAzElEQVRIiWNgGAWjYDACCQbGBx8qbOQMwDwDC6K0MBvOOJNmbMDADNIiQZQWNmnetsOJG8BaGIjQwj+7+THIlvTt7P1HN/wokGDgb+9OwG/JnWOGIL/k7uw5zHazB+gwiTNnN+DVYiCRYAyyJXfDjWS2GzxALQYSuYS0pH8D+SXdAKjl5h/itOSYgbQkgLTcJsoWiTtnikEOM9xw5rDZbRkDCR6CfuGf3b4R5H15g+ONz26++WMjx9/ei18LBuAhTfkoGAWjYBSMAqwAABXHS1JLAe9wAAAAAElFTkSuQmCC","orcid":"","institution":"Southeast University Zhongda Hospital Department of Hematology","correspondingAuthor":true,"submittingAuthor":false,"prefix":"","firstName":"Baoan","middleName":"","lastName":"Chen","suffix":""}],"badges":[],"createdAt":"2020-10-06 23:43:21","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-88926/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-88926/v1","draftVersion":[],"editorialEvents":[],"editorialNote":"","failedWorkflow":false,"files":[{"id":2894144,"identity":"409bbe61-fbb4-463e-9092-c952c35be979","added_by":"auto","created_at":"2020-10-09 18:46:45","extension":"jpg","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":95865,"visible":true,"origin":"","legend":"A schematic view of the procedure of the study with GSE73157","description":"","filename":"Figure1.JPG","url":"https://assets-eu.researchsquare.com/files/rs-88926/v1/d645ac1de30ef913f07ea94b.JPG"},{"id":2894145,"identity":"214fb8e2-0297-47da-9d68-841e8ade1eeb","added_by":"auto","created_at":"2020-10-09 18:46:45","extension":"jpg","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":153266,"visible":true,"origin":"","legend":"Identification of differentially expressed genes in GSE73157 dataset. Figure 2a Volcano plot of GSE73157 dataset. Red plots represent up-regulated genes and the green plots represent down-regulated ones with adjusted p value \u003c 0.05 and [log2FoldChange (log2FC)] \u003e 1. Other plots represent the remaining genes with no significant difference. Figure 2b Heatmap of the top 50 DEGs (50 up- and 50 down-regulated genes).","description":"","filename":"Figure2.JPG","url":"https://assets-eu.researchsquare.com/files/rs-88926/v1/dc3311a93c14620045f9d2f9.JPG"},{"id":2894146,"identity":"096a03dc-6087-4286-8b3c-946a07131916","added_by":"auto","created_at":"2020-10-09 18:46:46","extension":"jpeg","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":266561,"visible":true,"origin":"","legend":"GO and KEGG enrichment analysis a. The GO biological process (BP) enrichment analysis. b. The GO molecular function (MF) analysis. c. The GO cellular component (CC) enrichment analysis. d. The KEGG enrichment analysis. e. The interrelation between pathways of KEGG. f. Venn diagram showed the common gene of candidate genes. g. Numbers of genes enriched in the certain pathway. GO, Gene Ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes.","description":"","filename":"floatimage3.jpeg","url":"https://assets-eu.researchsquare.com/files/rs-88926/v1/eb6158fbaac20294f7703ea8.jpeg"},{"id":2894147,"identity":"131a79e7-4089-4db7-ade2-5e4ffbe5d186","added_by":"auto","created_at":"2020-10-09 18:46:46","extension":"jpeg","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":671557,"visible":true,"origin":"","legend":"Protein-protein interaction (PPI) analysis a. Genes identified from angiogenesis pathway. b. Genes identified from PI3K-Akt signaling pathway. c. Analysis of interaction of top 30 genes from angiogenesis pathway. d. Analysis of interaction of top 30 genes from PI3K-Akt signaling pathway.","description":"","filename":"floatimage4.jpeg","url":"https://assets-eu.researchsquare.com/files/rs-88926/v1/8a29b1c8e1b3fff15f76d1f2.jpeg"},{"id":2894148,"identity":"e809320d-8a13-4d76-889a-6de0378d081e","added_by":"auto","created_at":"2020-10-09 18:46:46","extension":"jpeg","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":96369,"visible":true,"origin":"","legend":"Analysis of selected hub gene ITGB1. a. Kaplan-Meier survival curves comparing high and low expression of ITGB1 in AML patients with PrognoScan (COX p value= 0.016730). b. Survival analysis curves comparing high and low expression of ITGB1 in AML patients with GEPIA (p value \u003c0.05). c. Comparation the expression of ITGB1 of AML patients and normal people (p value \u003c0.01).","description":"","filename":"floatimage5.jpeg","url":"https://assets-eu.researchsquare.com/files/rs-88926/v1/2ae66839ab1ae940c405bcd7.jpeg"},{"id":13601631,"identity":"0acaf8dd-4dea-4136-8659-3b86221768fa","added_by":"auto","created_at":"2021-09-17 05:48:21","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":932154,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-88926/v1/d4e82ba8-0393-427e-988e-a88d42c77af3.pdf"}],"financialInterests":"","formattedTitle":"\u003cp\u003eComprehensive Expression and Prognosis Analyses of ITGB1 in AML\u003c/p\u003e","fulltext":[{"header":"Background","content":"\u003cp\u003eLeukemia is a malignant clonal disease of the hematopoietic system, which accounts for 4% of all cancers[1]. Among them, acute myeloid leukemia (AML) is a common and dangerous type of leukemia. Currently, chemotherapy is still the most important option for the treatment of AML. However, the emergence of multidrug resistance (MDR) and recurrence greatly limits the effects of chemotherapy drugs. Therefore, the research on MDR of AML cells has always been a research hotspot in hematology.\u003c/p\u003e\n\u003cp\u003eIn recent years, studies have found that the bone marrow microenvironment not only plays an important role in the hematopoietic process, but is closely related to the occurrence and development of tumors. The bone marrow microenvironment contains multiple cell types, such as osteoblasts, mesenchymal stem cells, endothelial cells, neurons and so on. In normal humans, the bone marrow microenvironment affects hematopoietic stem cells mostly through the regulation of blood oxygen partial pressure and the secretion of cytokines[2]. However, in cancer patients, the occurrence of MDR and recurrence is inseparable from the role of bone marrow microenvironment. Leukemia and leukemia stem cells (LSCs) interact with the bone marrow microenvironment in a different way from hematopoietic stem cells to affect the development of leukemia[3]. LSCs can hide in the bone marrow niche to escape the killing effect of chemotherapy. As a result, the bone marrow niche has also become a key part of the recurrence and MDR of leukemia[4-6]. According to recent researches, we have known that MDR mediated by the bone marrow microenvironment might stem from the following several mechanisms: 1. The influence of cytokines[7, 8]; 2. Changes in adhesion of leukemia cells [9, 10]; 3. The regulation of expression of MDR-related genes[11]; 4. The adjustment of cell cycle and metabolism[12]. In addition, many factors can also be contributed to the recurrence of leukemia, such as the secretion of cytokines or the changes with certain pathways.\u003c/p\u003e\n\u003cp\u003eAs the emergence and rapid development of high-throughput platforms and microarray, numerous molecular heterogeneity on AML has been acknowledged. In order to find potential key biomarker of AML related to the bone marrow microenvironment-mediated poor prognosis, through the use of biological information technology and the analysis of high-throughput sequencing chip data, we screen out the differentially expressed genes (DEGs) of AML cells cultured alone and co-cultured with stromal cells, which were all treated with different concentrations of Arsenic trioxide (ATO). Furthermore, we pick up ITGB1 as a hub gene that can be considered as the predictor of prognosis. We systematically revel the occurrence and development of AML on the level of molecule and provide a potential guidance for targeted therapy.\u003c/p\u003e"},{"header":"Materials And Methods","content":"\u003cp\u003e\u003cstrong\u003e2.1 Analysis and screening of differential expressed genes (DEGs)\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe GEO database (\u003ca href=\"https://www.ncbi.nlm.nih.gov/geo/\"\u003ehttps://www.ncbi.nlm.nih.gov/geo/\u003c/a\u003e) of the National Center for Biotechnology Information was used to select and download the gene express profile GSE73157[13]. This profile was generated on the platform of GPL17077 (Agilent-039494 SurePrint G3 Human GE v2 8*60K Microarray 039381). In order to reduce false positive results, the online analysis tool GEO2R in the GEO database was used to process the data profile. The screen criteria were adjusted p value \u0026lt; 0.05 and the use of Benjamini-Hochberg. Then, R package \u0026ldquo;limma\u0026rdquo; was used to clarity the DEGs between NB4 cultured alone and co-cultured with stromal cells after the treatment of ATO. The screen criteria were adjusted p value \u0026lt; 0.05 and [log\u003csub\u003e2\u003c/sub\u003eFoldChange (log\u003csub\u003e2\u003c/sub\u003eFC)] \u0026gt; 1[14]. All genes were visualized by volcanic maps and top 50 dramatically differentially expressed genes were selected to draw a heatmap by R package \u0026ldquo;ggplot2\u0026rdquo;[15].\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003e2.2 Functional enrichment analysis and pathway analysis of DEGs\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eR package \u0026ldquo;clusterProfiler\u0026rdquo; was used to conduct Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis in order to preform functional enrichment and pathway analysis of DEGs[16]. Biological process (BP), molecular function (MF) and cellular component (CC) are included in GO enrichment. Analysis results were illustrated as figures by R package \u0026ldquo;GOplot\u0026rdquo;[17]. A p value less than 0.05 was considered statistically significant. Also, the ClueGo plug-ins of Cytoscape software 3.8.1 was utilized to display the relationship between pathways[18]. A p value less than 0.05 was considered statistically significant.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003e2.3 Screening candidate genes through protein-protein interaction (PPI) network\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe PPI prediction website STRING (\u003ca href=\"https://string-db.org\"\u003ehttps://string-db.org\u003c/a\u003e) constructed protein interaction network[19]. And Cytoscape software 3.8.0 was used to visualize the interactions of DEGs[18]. Further, we used \u0026ldquo;cytoHubba\u0026rdquo; to degree the interaction of hub-gene clustering by the method of Maxima Clique Centrality (MCC). Overlapped genes were separated by online toll Bioinformatics \u0026amp; Evolutionary Genomics (http://bioinformatics.psb.ugent.be/webtools/Venn/) and were displayed by Wayne diagram.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003e2.4 Selecting candidate genes and conducting survival analysis\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eIn order to evaluate the prognostic value of candidate genes in AML patients, we utilized the PrognoScan database (\u003ca href=\"http://dna00.bio.kyutech.ac.jp/PrognoScan/\"\u003ehttp://dna00.bio.kyutech.ac.jp/PrognoScan/\u003c/a\u003e) and online tool Gene Expression Profiling Interaction Analysis (\u003ca href=\"http://gepia.cancer-pku.cn/index.html\"\u003ehttp://gepia.cancer-pku.cn/index.html\u003c/a\u003e) to conduct the survival analysis[20, 21]. The results are displayed with hazard ratio (HR) and Cox p value from a Log-rank test. Cox p value \u0026lt; 0.05 was considered statistically significant. Besides, according to online tool Gene Expression Profiling Interaction Analysis (GEPIA), we compared the expression of candidate genes in AML patients and normal people. The screen criteria were adjusted p value \u0026lt; 0.01 and [log\u003csub\u003e2\u003c/sub\u003eFoldChange (log\u003csub\u003e2\u003c/sub\u003eFC)] \u0026gt; 1\u003c/p\u003e"},{"header":"Results","content":"\u003cp\u003e\u003cstrong\u003e3.1 Analysis and screening of differential expressed genes (DEGs)\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe goal of this study is to screen out the hub gene of AML related to the bone marrow microenvironment-mediated MDR (Figure 1). Perform statistical analysis on GEO database, we found that the data set GSE73157contains 32080 DEGs in AML cultured alone groups and co-cultured groups (All the following are referred to as alone groups and co-culture groups, respectively). By using R package \u0026ldquo;limma\u0026rdquo;, DEGs were classified as up-regulated genes (marked as red plots) and down-regulated genes (marked as green plots) (Figure 2a). Besides, a heatmap illustrated top 50 DEGs (Figure 2b).\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003e3.2 Functional enrichment analysis and pathway analysis of DEGs\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eGo and KEGG enrichment analysis was performed on all DEGs with R package \u0026ldquo;clusterProfiler\u0026rdquo;. The results of GO include three subontologies referred to as biological process (BP), cellular components (CC), molecular function (MF). The GO results show that angiogenesis pathway, cell adhesion molecule binding pathway and collagen-containing extracellular matrix pathway were selected as the most significant pathway in each subontologies, respectively (Figure 3a-c). And we chose angiogenesis pathway as our further analysis pathway as a result of comparation of p values. 177 genes were found enriched in this GO term. In addition, PI3K-Akt signaling pathway was the top enriched pathway of the DEGs in KEGG enrichment analysis (Figure 3d). After we used the ClueGo to further analyze the interrelation of the enriched pathways and the DEGs, PI3K-Akt signaling pathway was still the most significant pathway, and involved 120 DEGs (Figure 3e, f). Finally, we selected two pathways for our further research, which involved 177 and 120 DEGs, respectively (Table1). Then, after taking the intersection of the DEGs included in the two pathways, we could find that there are 26 genes included both in the two pathways (Figure 3g). They were PDGFRB, FLT1, NGFR, EFNA1, FGFR1, ITGB1, VEGFC, BRCA1, PDGFRA, PGF, EPHA2, EFNA3, PIK3CG, SYK, EREG, PTK2, ITGB8, FGF2, THBS1, FGF18, THBS2, AKT3, IL6, ANGPT1, FN1, ANGPT2.\u003c/p\u003e\n\u003ctable border=\"1\" width=\"557\"\u003e\n\u003ctbody\u003e\n\u003ctr\u003e\n\u003ctd width=\"113\"\u003e\n\u003cp\u003eDEGs\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"444\"\u003e\n\u003cp\u003eGenes\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"113\"\u003e\n\u003cp\u003eAngiogenesis pathway\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"444\"\u003e\n\u003cp\u003eABL1, ACVR1, ACVRL1, ADM, AKT3, AMOTL2, ANG, ANGPTL4, ANGPTL6, ANXA2, AQP1, ARHGAP22, C3, CARD10, CAV1, CCBE1, CCL2, CD40, CDH13, CELA1, CHI3L1, CHRNA7, COL18A1, COL8A2, CYP1B1, DCN, DDAH1, DLL4, E2F7, ECM1, EFNA1, EFNA3, EGFL7, ELK3, EMP2, ENPP2, EPAS1, EPHA2, EPHB2, EREG, ESM1, FAP, FGF2, FLT1, FN1, FZD5, GATA6, GBX2, GPNMB, GPR4, GREM1, HDAC9, HEY1, HLA-G, HMOX1, HOXB13, HPSE, HTATIP2, ID1, IL18, IL1A, IL1B, IL6, ITGAX, ITGB1, ITGB2, ITGB8, KLF5, LEF1, LOXL2, LRG1, MAPK14, MDK, MFGE8, MMP14, MMP2, NFATC4, NGFR, NOTCH3, NRP1, NRP2, PARVA, PDGFRA, PDGFRB, PGF, PGK1, PLCD1, PLK2, PPARG, PRKD2, PTGS2, PTK2, PTN, PTPRM, RAMP3, RBM15, RHOB, RNF213, RORA, RRAS, SASH1, SAT, SERPINE1, SERPINF1, SFRP1, SFRP2, SMAD1, SOX17, SPARC, SRPX2, SYNJ2BP, TGFBI, THBS1, THBS2, THY1, TIE1, TNFAIP3, TNFSF12, TWIST1, VASH1, VEGFC, WNT5A, ADM2, AGGF1, ANGPT1, ANGPT2, BRCA1, CCR2, CX3CL1, CXCR2, CYBB, DAB2IP, DAG1, ECSCR, EFNB2, EMILIN1, EPHA1, F3, FBXW7, FGF18, FGFBP1, FGFR1, FOXJ2, GATA2, GPLD1, HHEX, HOXA7, IHH, JAG1, JUN, KRIT1, LEMD3, MAPK7, MED1, NCL, NOTCH4, NPPB, NRAR, PDE3B, PIK3CG, PKNOX1, POFUT1, PRKX, ROCK2, SETD2, SIRT1, SOX18, SPRED1, SPRY2, SRF, STAB1, SYK, TBX1, TERT, TGFBR1, THSD7A, UBP1\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"113\"\u003e\n\u003cp\u003ePI3K-Akt signaling pathway\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"444\"\u003e\n\u003cp\u003eTCL1B, GNGT2, OSMR, GNG11, GNG10, IGF2, IGF1R, LAMB2, LAMB3, BCL2, GNG7, F2R,, LAMA1, COL9A3, MYB, LAMA4, EIF4EBP1, LAMB1, ITGA11, THBS1, PDGFRA, THBS2, TNC, GNG12, pik3ap1, INSR, FGF18, LPAR1, RXRA, MET, HRAS, LPAR5, HSP90B1, AREG, KIT, CCND2, PDGFC, FLT3, RPTOR, FLT1, PDGFRB, LPAR4, PIK3CG, PIK3R1, CCNE2, ITGA1, IKBKB, FN1, AKT3, SOS2, IRS1, CREB3, CDC37, KITLG, NRAS, IL3RA, Ifna4, CDK6, IL2RB, NGFR, CDK4, BRCA1, IL2RA, JAK2, THEM4, JAK3, EREG, ITGA2B, ITGA3, TLR4, ITGB4, MAPK3, ITGB5, CREB3L2, ITGB7, PIK3R5, IL6R, ANGPT1, TLR2, ANGPT2, EFNA1, GYS2, IL6, CDKN1A, TSC2ITGB1, PGF, YWHAB, PTK2, EFNA3, NFKB1, SPP1, SYK, VEGFC, CSF3, GNG2, FGF2, LAMC3, IFNB1\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003c/tbody\u003e\n\u003c/table\u003e\n\u003cp\u003e\u003cstrong\u003eTable 1\u003c/strong\u003e: DEGs identified from selected pathways of GO and KEGG\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003e3.3 Screening candidate genes through protein-protein interaction (PPI) network\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eIn order to further understand the DEGs, we displayed the protein-protein interaction (PPI) network including all the DEGs in the two selected pathways mentioned above in STRING (Figure 4a, b). Then, cytoHubba plug-ins in Cytoscape was also used to screen out top 30 candidate hub genes of each pathway according nodes rank (Figure 4c, d), and all the top hub genes has been listed in Table 2. In this part, we identified 4 common genes in the two sets of top 30 hub genes, including FN1, THBS1, ITGB1 and THBS2 as candidate hub genes.\u003c/p\u003e\n\u003ctable border=\"1\" width=\"553\"\u003e\n\u003ctbody\u003e\n\u003ctr\u003e\n\u003ctd colspan=\"3\" width=\"287\"\u003e\n\u003cp\u003e\u003cstrong\u003eAngiogenesis pathway\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd colspan=\"3\" width=\"266\"\u003e\n\u003cp\u003e\u003cstrong\u003ePI3K-Akt signaling pathway\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003eRank\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"121\"\u003e\n\u003cp\u003eName\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003eScore\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003eRank\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eName\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003eScore\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e1\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"121\"\u003e\n\u003cp\u003eFN1\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e1.68E+08\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e1\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003ePTK2\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e7.42E+14\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e2\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"121\"\u003e\n\u003cp\u003eCCL2\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e1.67E+08\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e2\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eITGB1\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e7.42E+14\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e3\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"121\"\u003e\n\u003cp\u003eMMP2\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e1.65E+08\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e3\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eITGA3\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e7.42E+14\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e4\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"121\"\u003e\n\u003cp\u003eFGF2\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e1.58E+08\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e4\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eITGA1\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e7.42E+14\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e5\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"121\"\u003e\n\u003cp\u003eSERPINE1\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e1.38E+08\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e5\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eITGA2B\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e7.42E+14\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e6\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"121\"\u003e\n\u003cp\u003eTHBS1\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e1.33E+08\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e6\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eITGB5\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e7.42E+14\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e7\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"121\"\u003e\n\u003cp\u003eIL6\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e1.24E+08\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e6\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eITGA11\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e7.42E+14\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e8\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"121\"\u003e\n\u003cp\u003eVEGFC\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e1.21E+08\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e8\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eITGB7\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e7.42E+14\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e9\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"121\"\u003e\n\u003cp\u003eFLT1\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e1.07E+08\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e9\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eITGB4\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e7.40E+14\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e10\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"121\"\u003e\n\u003cp\u003ePGF\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e1.07E+08\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e10\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eITGB8\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e7.40E+14\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e11\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"121\"\u003e\n\u003cp\u003eANGPT2\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e1.04E+08\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e10\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eITGA10\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e7.40E+14\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e12\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"121\"\u003e\n\u003cp\u003eANGPT1\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e1.00E+08\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e12\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eCOL1A1\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e7.18E+14\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e13\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"121\"\u003e\n\u003cp\u003ePTGS2\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e6.38E+07\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e13\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eCOL1A2\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e7.17E+14\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e14\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"121\"\u003e\n\u003cp\u003eCOL18A1\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e4.98E+07\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e14\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eCOL6A1\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e7.17E+14\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e15\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"121\"\u003e\n\u003cp\u003eIL1B\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e4.08E+07\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e15\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eCOL6A3\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e7.15E+14\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e16\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"121\"\u003e\n\u003cp\u003eJUN\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e3.63E+07\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e16\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eCOL6A2\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e7.14E+14\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e17\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"121\"\u003e\n\u003cp\u003eHMOX1\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e3.31E+07\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e17\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eLAMA4\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e3.83E+14\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e18\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"121\"\u003e\n\u003cp\u003eMMP14\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e2.35E+07\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e18\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eLAMB1\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e3.82E+14\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e19\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"121\"\u003e\n\u003cp\u003ePPARG\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e2.18E+07\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e19\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eTHBS2\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e3.58E+14\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e20\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"121\"\u003e\n\u003cp\u003eCAV1\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e1.86E+07\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e20\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eCOL9A3\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e3.56E+14\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e21\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"121\"\u003e\n\u003cp\u003eIL18\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e1.73E+07\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e21\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eLAMB2\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e2.37E+13\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e22\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"121\"\u003e\n\u003cp\u003eMAPK14\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e1.58E+07\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e22\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eLAMA1\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e2.36E+13\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e23\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"121\"\u003e\n\u003cp\u003eITGB1\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e1.34E+07\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e23\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eLAMC3\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e2.10E+13\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e24\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"121\"\u003e\n\u003cp\u003ePDGFRB\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e1.30E+07\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e24\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eLAMB3\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e2.09E+13\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e25\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"121\"\u003e\n\u003cp\u003eSIRT1\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e1.10E+07\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e25\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eTHBS1\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e3.95E+12\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e26\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"121\"\u003e\n\u003cp\u003eF3\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e7622520\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e26\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eSPP1\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e1.58E+12\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e27\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"121\"\u003e\n\u003cp\u003eANG\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e7258344\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e27\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eFN1\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e1.43E+12\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e28\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"121\"\u003e\n\u003cp\u003eTHBS2\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e6250378\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e28\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003eTNC\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e9.39E+10\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e29\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"121\"\u003e\n\u003cp\u003eIL1A\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e5449080\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e29\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003ePIK3CG\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e6.24E+09\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e30\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"121\"\u003e\n\u003cp\u003eCYBB\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e3774960\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"61\"\u003e\n\u003cp\u003e30\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"100\"\u003e\n\u003cp\u003ePIK3R1\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd width=\"105\"\u003e\n\u003cp\u003e6.24E+09\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003c/tbody\u003e\n\u003c/table\u003e\n\u003cp\u003e\u003cstrong\u003eTable 2: \u003c/strong\u003eThe top 30 genes with the highest score of each pathway through the Cytoscape \u0026ldquo;cytoHubba\u0026rdquo; module analysis\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003e3.4 Selecting candidate genes and conducting survival analysis\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eSubsequently, the correlation between hub genes and prognosis of AML patients was analyzed through PrognoScan and GEPIA online survival website. The results showed that except for ITGB1 gene, as the result of no survival significance, the other genes were removed. We found that patients with high expression of ITGB1 had significantly longer overall survival (OS) than those with low expression (COX p value= 0.016730) (figure 5a, b). Besides, through the online website GEPIA, we performed a visual analysis of the expression of hub gene ITGB1 in the TCGA database in AML patients and normal people, showing that the expression of the ITGB1 gene in AML patients is lower than that in normal people significantly (p value\u0026lt;0.01) (Figure 5c).\u003c/p\u003e"},{"header":"Discussion","content":"\u003cp\u003eWith the elucidation of molecular mechanisms and the emergence of targeted therapies, the treatment of AML has made great progress, and the overall survival rate of patients has been significantly improved. However, the resistance of leukemia cells and the recurrence after treatment are still two major problems that plague us. We now believe that the protective effect of leukemia stem cells (LSCs) mediated by the bone marrow microenvironment is also an important factor for the refractory and recurrence of leukemia. Recent studies have shown that there is a complex interaction between leukemia cells and the bone marrow microenvironment. The two induce and support each other, which together lead to the occurrence and development of the disease.\u003c/p\u003e\n\u003cp\u003eIn this study, we analyzed 32080 DEGs between cultured alone and co-cultured with stomal cells AML cells based on GSE73157 dataset. Then, we used GO enrichment analysis and KEGG pathway analysis to detect the dataset to explore the interaction between DEGs. GO analysis showed that these genes mainly reacted with the process of angiogenesis, collagen-containing extracellular matrix production and cytokine binding, suggesting that these DEGs may be involved in improving the ability of tumor cell to migrate and invade, promoting metastasis, recurrence and multiple drug resistance of cancer[22, 23]. KEGG analysis illustrated that 120 DEGs mainly affect the AML cells through PI3K-Akt signaling pathway. As we all known, \u0026nbsp;phosphatidylinositol 3-kinase/protein kinase B/mammalian target of rapamycin (PI3K/Akt/mTOR) signaling pathway is over activated in hematological malignancies, and the development of inhibitors of PI3K/Akt/mTOR signaling pathway has become a common concern for hematology researchers in recent years[24]. Excessive activation of this pathway may promote the tumor proliferation and weaken the immune monitoring of tumor cells[25]. At the same time, we utilized the online tool STRING to display the functional interaction between proteins to discover the underlying mechanism of occurrence of development of AML. And then using the software Cytoscape 3.8.1 to analyze the protein interaction network diagram, we screened out FN1, THBS1, ITGB1, THBS2 as 4 candidate hub genes. The above 4 candidate genes not only interact closely with other genes, but may also determine the functions of other genes. After further survival analysis of the 4 selected genes, we found that only ITGB1 gene is statistically significant. We found that the expression of ITGB1 is negatively correlated with the prognosis of AML. Furthermore, compared with normal people, the expression of ITGB1 is significantly lower in AML patients. That is, ITGB1 may be a key gene for evaluating the prognosis and conducting targeted therapy of AML patients.\u003c/p\u003e\n\u003cp\u003eA huge amount of evidences show that integrin family molecules are involved in the occurrence and development of tumors, and play an important role in regulating tumor cell proliferation, invasion and metastasis and other biological behaviors[26]. Integrin \u0026beta;1 (ITGB1) is a relatively important part of molecules in the cell adhesion molecule family, and is the main receptor that mediates the interaction between cells and extracellular matrix. This saying is consistent with our above analysis results: the bone marrow microenvironment may affect AML cells by regulating cell adhesion molecules and extracellular matrix. At present, our understanding of ITGB1 is limited to the primary stage, and its biological functions are not yet very clear. In the direction of solid tumors, there are certain studies on the direction of ITGB1, but there are few reports on the research of ITGB1 of hematological tumors in the world. In summary, ITGB1 may be a potential marker for predicting and guiding the treatment of AML, which requires further attention and research.\u003c/p\u003e"},{"header":"Conclusion","content":"\u003cp\u003eIn summary, ITGB1 was a hub gene related to the bone microenvironment mediated poor prognosis in AML. The expression of ITGB1 was negatively correlated with prognosis of AML. This study used bioinformatics methods to analyze the relevant data of AML researches to obtain potential markers for predicting and guiding AML treatment, which provides a theoretical basis for a deeper understanding of the occurrence and development of AML, and also give the follow-up animal experiments and clinical trials research directions.\u003c/p\u003e"},{"header":"Abbreviations","content":"\u003cp\u003eAML: Acute myeloid leukemia\u003c/p\u003e\n\u003cp\u003eMDR: Multidrug resistance\u003c/p\u003e\n\u003cp\u003eITGB1: Integrin \u0026beta;1\u003c/p\u003e\n\u003cp\u003eNCBI-GEO: National Center for Biotechnology Information-Gene Expression Omnibus database\u003c/p\u003e\n\u003cp\u003eDEGs: Differentially expressed genes\u003c/p\u003e\n\u003cp\u003eGO: Gene Ontology\u003c/p\u003e\n\u003cp\u003eKEGG: Kyoto Encyclopedia of Genes and Genomes\u003c/p\u003e\n\u003cp\u003ePPI: Protein-protein interaction\u003c/p\u003e\n\u003cp\u003eOS: Overall survival\u003c/p\u003e\n\u003cp\u003eLSCs: Leukemia stem cells\u003c/p\u003e\n\u003cp\u003eATO: Arsenic trioxide\u003c/p\u003e\n\u003cp\u003eBP: Biological process\u003c/p\u003e\n\u003cp\u003eMF: Molecular function\u003c/p\u003e\n\u003cp\u003eCC: Cellular component\u003c/p\u003e\n\u003cp\u003eMCC: Maxima Clique Centrality\u003c/p\u003e\n\u003cp\u003eHR: Hazard ratio\u003c/p\u003e\n\u003cp\u003eGEPIA: Gene Expression Profiling Interaction Analysis\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eEthics approval and consent to participate\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eNot applicable.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eConsent for publication\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eNot applicable.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAvailability of data and materials\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe dataset analysed during the current study are available in the NCBI-GEO repository, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE73157.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eCompeting interests \u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe authors report no conflicts of interest in this work.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eFunding \u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThis work was supported by The Natural Science Foundation of Jiangsu Province (BK20150639); Key Medical of Jiangsu Province (ZDXKB2016020); Six talent peaks project in Jiangsu Province (WSW-033);Innovative and entrepreneurial doctors of Jiangsu province.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAuthors' contributions\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eXYZ, NJ and BAC were responsible for confirming the topic, collecting the data and analyzing article. XYZ designed and wrote the paper, edited the figures, legends and tables. All authors read and approved the final manuscript.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAcknowledgements\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eNot applicable.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\n\u003cli\u003eSiegel, R.L., K.D. Miller, and A. Jemal, \u003cem\u003eCancer statistics, 2020.\u003c/em\u003e CA Cancer J Clin, 2020. \u003cstrong\u003e70\u003c/strong\u003e(1): p. 7-30.\u003c/li\u003e\n\u003cli\u003eKfoury, Y. and D.T. Scadden, \u003cem\u003eMesenchymal cell contributions to the stem cell niche.\u003c/em\u003e Cell Stem Cell, 2015. \u003cstrong\u003e16\u003c/strong\u003e(3): p. 239-53.\u003c/li\u003e\n\u003cli\u003eKrause, D.S., et al., \u003cem\u003eDifferential regulation of myeloid leukemias by the bone marrow microenvironment.\u003c/em\u003e Nat Med, 2013. \u003cstrong\u003e19\u003c/strong\u003e(11): p. 1513-7.\u003c/li\u003e\n\u003cli\u003eFrisch, B.J., et al., \u003cem\u003eFunctional inhibition of osteoblastic cells in an in vivo mouse model of myeloid leukemia.\u003c/em\u003e Blood, 2012. \u003cstrong\u003e119\u003c/strong\u003e(2): p. 540-50.\u003c/li\u003e\n\u003cli\u003eFulzele, K., et al., \u003cem\u003eMyelopoiesis is regulated by osteocytes through Gs\u0026alpha;-dependent signaling.\u003c/em\u003e Blood, 2013. \u003cstrong\u003e121\u003c/strong\u003e(6): p. 930-9.\u003c/li\u003e\n\u003cli\u003eWelner, R.S., et al., \u003cem\u003eTreatment of chronic myelogenous leukemia by blocking cytokine alterations found in normal stem and progenitor cells.\u003c/em\u003e Cancer Cell, 2015. \u003cstrong\u003e27\u003c/strong\u003e(5): p. 671-81.\u003c/li\u003e\n\u003cli\u003eFlavell, R.A., et al., \u003cem\u003eThe polarization of immune cells in the tumour environment by TGFbeta.\u003c/em\u003e Nat Rev Immunol, 2010. \u003cstrong\u003e10\u003c/strong\u003e(8): p. 554-67.\u003c/li\u003e\n\u003cli\u003eMassagu\u0026eacute;, J., \u003cem\u003eTGFbeta in Cancer.\u003c/em\u003e Cell, 2008. \u003cstrong\u003e134\u003c/strong\u003e(2): p. 215-30.\u003c/li\u003e\n\u003cli\u003eSugiyama, T., et al., \u003cem\u003eMaintenance of the hematopoietic stem cell pool by CXCL12-CXCR4 chemokine signaling in bone marrow stromal cell niches.\u003c/em\u003e Immunity, 2006. \u003cstrong\u003e25\u003c/strong\u003e(6): p. 977-88.\u003c/li\u003e\n\u003cli\u003eDillmann, F., et al., \u003cem\u003ePlerixafor inhibits chemotaxis toward SDF-1 and CXCR4-mediated stroma contact in a dose-dependent manner resulting in increased susceptibility of BCR-ABL+ cell to Imatinib and Nilotinib.\u003c/em\u003e Leuk Lymphoma, 2009. \u003cstrong\u003e50\u003c/strong\u003e(10): p. 1676-86.\u003c/li\u003e\n\u003cli\u003eHazlehurst, L.A., T.H. Landowski, and W.S. Dalton, \u003cem\u003eRole of the tumor microenvironment in mediating de novo resistance to drugs and physiological mediators of cell death.\u003c/em\u003e Oncogene, 2003. \u003cstrong\u003e22\u003c/strong\u003e(47): p. 7396-402.\u003c/li\u003e\n\u003cli\u003eColmone, A., et al., \u003cem\u003eLeukemic cells create bone marrow niches that disrupt the behavior of normal hematopoietic progenitor cells.\u003c/em\u003e Science, 2008. \u003cstrong\u003e322\u003c/strong\u003e(5909): p. 1861-5.\u003c/li\u003e\n\u003cli\u003eBarrett, T., et al., \u003cem\u003eNCBI GEO: archive for functional genomics data sets--update.\u003c/em\u003e Nucleic Acids Res, 2013. \u003cstrong\u003e41\u003c/strong\u003e(Database issue): p. D991-5.\u003c/li\u003e\n\u003cli\u003eRitchie, M.E., et al., \u003cem\u003elimma powers differential expression analyses for RNA-sequencing and microarray studies.\u003c/em\u003e Nucleic Acids Res, 2015. \u003cstrong\u003e43\u003c/strong\u003e(7): p. e47.\u003c/li\u003e\n\u003cli\u003eIto, K. and D. Murphy, \u003cem\u003eApplication of ggplot2 to Pharmacometric Graphics.\u003c/em\u003e CPT Pharmacometrics Syst Pharmacol, 2013. \u003cstrong\u003e2\u003c/strong\u003e(10): p. e79.\u003c/li\u003e\n\u003cli\u003eYu, G., et al., \u003cem\u003eclusterProfiler: an R package for comparing biological themes among gene clusters.\u003c/em\u003e Omics, 2012. \u003cstrong\u003e16\u003c/strong\u003e(5): p. 284-7.\u003c/li\u003e\n\u003cli\u003eWalter, W., F. S\u0026aacute;nchez-Cabo, and M. Ricote, \u003cem\u003eGOplot: an R package for visually combining expression data with functional analysis.\u003c/em\u003e Bioinformatics, 2015. \u003cstrong\u003e31\u003c/strong\u003e(17): p. 2912-4.\u003c/li\u003e\n\u003cli\u003eShannon, P., et al., \u003cem\u003eCytoscape: a software environment for integrated models of biomolecular interaction networks.\u003c/em\u003e Genome Res, 2003. \u003cstrong\u003e13\u003c/strong\u003e(11): p. 2498-504.\u003c/li\u003e\n\u003cli\u003eSzklarczyk, D., et al., \u003cem\u003eSTRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets.\u003c/em\u003e Nucleic Acids Res, 2019. \u003cstrong\u003e47\u003c/strong\u003e(D1): p. D607-d613.\u003c/li\u003e\n\u003cli\u003eMizuno, H., et al., \u003cem\u003ePrognoScan: a new database for meta-analysis of the prognostic value of genes.\u003c/em\u003e BMC Med Genomics, 2009. \u003cstrong\u003e2\u003c/strong\u003e: p. 18.\u003c/li\u003e\n\u003cli\u003eTang, Z., et al., \u003cem\u003eGEPIA: a web server for cancer and normal gene expression profiling and interactive analyses.\u003c/em\u003e Nucleic Acids Res, 2017. \u003cstrong\u003e45\u003c/strong\u003e(W1): p. W98-w102.\u003c/li\u003e\n\u003cli\u003eLi, T., et al., \u003cem\u003eTumor angiogenesis and anti-angiogenic gene therapy for cancer.\u003c/em\u003e Oncol Lett, 2018. \u003cstrong\u003e16\u003c/strong\u003e(1): p. 687-702.\u003c/li\u003e\n\u003cli\u003eTerranova, V.P., D. Maslow, and G. Markus, \u003cem\u003eDirected migration of murine and human tumor cells to collagenases and other proteases.\u003c/em\u003e Cancer Res, 1989. \u003cstrong\u003e49\u003c/strong\u003e(17): p. 4835-41.\u003c/li\u003e\n\u003cli\u003eGao, Y., C.Y. Yuan, and W. Yuan, \u003cem\u003eWill targeting PI3K/Akt/mTOR signaling work in hematopoietic malignancies?\u003c/em\u003e Stem Cell Investig, 2016. \u003cstrong\u003e3\u003c/strong\u003e: p. 31.\u003c/li\u003e\n\u003cli\u003eO'Donnell, J.S., et al., \u003cem\u003ePI3K-AKT-mTOR inhibition in cancer immunotherapy, redux.\u003c/em\u003e Semin Cancer Biol, 2018. \u003cstrong\u003e48\u003c/strong\u003e: p. 91-103.\u003c/li\u003e\n\u003cli\u003eKapp, T.G., et al., \u003cem\u003eIntegrin modulators: a patent review.\u003c/em\u003e Expert Opin Ther Pat, 2013. \u003cstrong\u003e23\u003c/strong\u003e(10): p. 1273-95.\u003c/li\u003e\n\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":true,"highlight":"","institution":"","isAcceptedByJournal":false,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"
[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true},"keywords":"Acute myeloid leukemia, the bone marrow microenvironment, Bioinformatics analysis","lastPublishedDoi":"10.21203/rs.3.rs-88926/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-88926/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eBackground\u003c/p\u003e\u003cp\u003eAcute myeloid leukemia (AML) is a dangerous type of leukemia. The emergence of multidrug resistance (MDR) and recurrence limits the prognosis and survival of patients. In recent studies, we have known that the bone marrow microenvironment was closely related to the poor prognosis of AML. However, the underlying mechanisms are still far from fully understood. By utilizing the bioinformatics analysis, we screened out integrin β1 (ITGB1) as the hub gene, which is associated with the bone marrow microenvironment mediated changes of AML cells, with expression profile GSE73157 downloaded from National Center for Biotechnology Information-Gene Expression Omnibus (NCBI-GEO) database.\u003c/p\u003e\u003cp\u003eMethods\u003c/p\u003e\u003cp\u003eR studio software was used to screen out candidate hub genes and further visualize the differential expression. R package “limma” was to find out differentially expressed genes (DEGs). Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway were conducted by R package “cluster Profiler”. Furthermore, protein-protein interaction (PPI) network was also performed by online tool STRING and software Cytoscape. Last but not least, online tool PrognoScan and GEPIA was utilized for the evaluation of clinical significance of the selected hub gene. P and Cox p value \u0026lt;0.05 was considered to be statistical significance.\u003c/p\u003e\u003cp\u003eResults\u003c/p\u003e\u003cp\u003eITGB1 was filtrated as the only hub gene in this profile. We found that patients with high expression of ITGB1 had significantly longer overall survival (OS) than those with low expression (COX p value= 0.016730). Besides, the expression of the ITGB1 gene in AML patients is lower than that in normal people significantly (p value\u0026lt;0.01).\u003c/p\u003e\u003cp\u003eConclusion\u003c/p\u003e\u003cp\u003eWe identified ITGB1 as a key gene in the bone marrow microenvironment mediated poor prognosis in AML. The down-regulated expression of ITGB1 was related to AML patients’ poor outcome. ITGB1 may be a potential marker for predicting and guiding AML treatment.\u003c/p\u003e","manuscriptTitle":"Comprehensive Expression and Prognosis Analyses of ITGB1 in AML","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2020-10-09 18:46:43","doi":"10.21203/rs.3.rs-88926/v1","editorialEvents":[{"type":"communityComments","content":0}],"status":"published","journal":{"display":true,"email":"
[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"e73c5019-ad9a-45cc-b69a-a4aed3f661bd","owner":[],"postedDate":"October 9th, 2020","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"posted","subjectAreas":[{"id":734299,"name":"Cancer Biology"},{"id":734300,"name":"Oncology"}],"tags":[],"updatedAt":"2020-10-30T21:11:54+00:00","versionOfRecord":[],"versionCreatedAt":"2020-10-09 18:46:43","video":"","vorDoi":"","vorDoiUrl":"","workflowStages":[]},"version":"v1","identity":"rs-88926","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-88926","identity":"rs-88926","version":["v1"]},"buildId":"WrCJVZZCHTDjtuVLN7oU0","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}
Text is read by the "Ask this paper" AI Q&A widget below.
Extraction quality varies by source — PMC NXML preserves structure
cleanly, OA-HTML may include some navigation residue, and OA-PDF can
have broken hyphenation. The publisher copy
(via DOI)
is the canonical version.