Wufuer Aini1,2,3†Limin Xie1,2,3†Wanyu Hu1,2,3Yuan Tang1,2,3Hongling Peng4,5,6Guangsen Zhang5
Tuo Deng1,2,3*- 1National Clinical Research Center for Metabolic Diseases, and Department of Metabolism and Endocrinology, The Second Xiangya Hospital of Central South University, Changsha, Hunan, China
- 2Key Laboratory of Diabetes Immunology, Ministry of Education, and Metabolic Syndrome Research Center, The Second Xiangya Hospital of Central South University, Changsha, Hunan, China
- 3Clinical Immunology Center, The Second Xiangya Hospital of Central South University, Changsha, Hunan, China
- 4Department of Hematology, The Second Xiangya Hospital of Central South University, Changsha, Hunan, China
- 5Institute of Molecular Hematopathy, The Second Xiangya Hospital of Central South University, Changsha, Hunan, China
- 6Hunan Key Laboratory of Tumor Models and Individualized Medicine, Changsha, Hunan, China
Background: Polycythemia Vera (PV) is a type of typical Myeloproliferative Neoplasms (MPNs) characterized with excessive erythropoiesis and thrombosis. Anoikis is a special programmed cell death mode induced by the adhesion disorder between cells and extracellular matrix (ECM) or adjacent cells facilitating cancer metastasis. However, few studies have focused on the role of anoikis in PV, especially on the development of PV.
Methods: The microarray and RNA-seq results were screened from the Gene Expression Omnibus (GEO) database and the anoikis-related genes (ARGs) were downloaded from Genecards. The functional enrichment analysis of intersecting differentially expressed genes (DEGs) and protein-protein interaction (PPI) network analysis were performed to discover hub genes. The hub genes expression was tested in the training (GSE136335) and validation cohort (GSE145802), and RT-qPCR was performed to verify the gene expression in PV mice.
Results: In the training GSE136335, a total of 1,195 DEGs was obtained from Myeloproliferative Neoplasm (MPN) patients compared with controls, among which 58 were anoikis-related DEGs. The significant enrichment of the apoptosis and cell adhesion pathways (i.e., cadherin binding) were shown in functional enrichment analysis. The PPI network was conducted to identify top five hub genes (CASP3, CYCS, HIF1A, IL1B, MCL1). The expression of CASP3 and IL1B were significantly upregulated both in validation cohort and PV mice and downregulated after treatment, suggesting that CASP3 and IL1B could be important indicators for disease surveillance.
Conclusion: Our research revealed a relationship between anoikis and PV for the first time by combined analysis of gene level, protein interaction and functional enrichment, allowing novel insights into mechanisms of PV. Moreover, CASP3 and IL1B may become promising indicators of PV development and treatment.
Introduction
Myeloproliferative Neoplasms (MPNs) are the abnormal clonal myeloproliferation arising from the hematopoietic stem cells with acquired somatic mutation, and constitutional systematic phenotypes (Cervantes et al., 2008; Mead and Mullally, 2017). The latest World Health Organization (WHO) classified MPNs into chronic myeloid leukemia (CML) BCR::ABL1+, polycythemia vera (PV), essential thrombocythemia (ET), and primary myelofibrosis (PMF) (Arber et al., 2016). The latter three are generally called BCR::ABL1- MPNs (hereinafter referred to as MPNs).
PV, the commonest MPN, is characterized with erythrocytosis in consequence of JAK2 exon 14 or 12 somatic mutations (Baxter et al., 2005; James et al., 2005; Kralovics et al., 2005; Levine et al., 2005; Scott et al., 2007). Unlike ET or PMF, PV is characterized by progressive increase in erythropoiesis, splenomegaly due to extramedullary hematopoiesis, and arterial and venous thrombosis (Spivak, 2017). In addition, some patients with PV are at risk of transformation to bone marrow failure, myelofibrosis, and acute leukemia (Spivak, 2017). Ruxolitinib, a typical JAK1/2 inhibitor, was used to provide systemic symptom relief in post-PV myelofibrosis (PPMF) (Verstovsek et al., 2012) and chronic-phase PV (Vannucchi et al., 2015; Passamonti et al., 2017), while did not promote apoptosis of hematopoietic stem cells (HSCs) and failed to achieve molecular remissions (Wang et al., 2014).
Anoikis, a type of programmed apoptosis, is induced by incorrect attachment between cells and extracellular matrix (ECM) or adjacent cells, playing important role in eradicating misplaced or dislodged cells to maintain tissue homeostasis (Han et al., 2021). Some studies suggested anoikis occurred through intrinsic/extrinsic pathways and could be triggered by several intracellular signals such as DNA damage and endoplasmic reticulum stress (Li et al., 2019; Zhi et al., 2022). Anoikis is essential for cancer cell survival through the detachment from ECM, and anoikis resistance is a key factor in cancer invasion and metastasis (Berezovskaya et al., 2005; Kim et al., 2012; Kakavandi et al., 2018; Xiao et al., 2019; Yu et al., 2022; Zhou et al., 2022).
Clonal cell proliferation leads to increased risk of thrombosis and threatens the life of PV and ET patients (Barbui et al., 2013). Whereas, complications of thrombosis in myelofibrosis (MF) are as frequent as in ET but less than in PV (Kc et al., 2017). Recently, the abnormal expression of adhesion molecules like Lu/BCAM and CD147 was reported in PV patients, which was reduced under interferon-α treatment but enhanced by hydroxyurea (Brusson et al., 2018). These results indicated the important role of cell adhesion in PV and PV treatment. However, few studies had systematically evaluated the effect of anoikis in hematological cancers especially MPNs.
Therefore, in this study, we focused on exploring the expression levels of anoikis-related genes in the occurrence, progression and treatment of MPNs, as well as exploring the key anoikis-related genes by means of bioinformatics analysis.
Materials and methods
Data collection
Key word “myeloproliferative neoplasms” was used to search gene expression profiles of MPNs in the Gene Expression Omnibus (GEO, http://www.ncbi.nlm.nih.gov/geo/) data portal. The following criteria were used to filter the obtained database: 1. The gene expression profiling must cover all types of chronic and progressed MPNs (PV, ET, PMF) and controls. 2. The cell used for sequencing should be CD34+ MNCs (mononuclear cell) from peripheral blood (PB) or bone marrow (BM). 3. The processed data or raw data from datasets must be provided thus could be available to reanalyze. 4. The GEO dataset numbered GSE136335, GSE145802 were selected. Further, anoikis-related genes (ARGs) were downloaded from the GeneCard database (https://www.genecards.org/).
Design of this study
The study design was shown in Figure 1. Firstly, differentially expressed genes (DEGs) were screened in MPNs and MPN sub types including PV, ET, PMF and post-PV myelofibrosis (pPVMF), as well as the intersection of the DEGs and ARGs was analyzed. And then functional enrichment GO and KEGG analyses were performed. Lastly, PPI network analysis was conducted and the top 5 hub genes were found in the network. In addition, the corresponding hub genes in different validation cohorts were validated and the hub genes of PV in our mice model were detected.
FIGURE 1. Study flowchart. GSE, Gene Expression Omnibus Series; GO, gene ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes; PV, Polycythemia Vera; BM, Bone Marrow; MEPs, megakaryocyte and erythroid progenitors.
Data processing and screening of anoikis-related DEGs
Differentially expressed genes (DEGs) between patients with MPNs and healthy controls were identified in GSE136335 dataset via the “limma” R package (Gentleman, 2005), the cutoff for adjustment: p-value <0.05 and FC (fold change) > 1.5. The “ggplot2” R package (Ito and Murphy, 2013) was used to visualize the DEGs into a volcano map, while the “pheatmap” R package (Szekely and Rizzo, 2005) was used to cluster the anoikis-related DEGs. The intersection between DEGs and ARGs was visualized by the Veen plot.
Functional enrichment analysis and protein-protein interaction network analysis of anoikis-related DEGs.
Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of anoikis-related DEGs were performed by using “clusterProfiler” R package (Wu et al., 2021). GO analysis is important to explore the biological functions, consisting of three subtypes, including biological process (BP), cellular component (CC), and molecular function (MF) (Ashburner et al., 2000). KEGG analysis is used to explore potential pathways (Kanehisa and Goto, 2000) The STRING online website (https://string-db.org/) was to conduct the protein-protein interaction (PPI) network analysis of anoikis-related DEGs. The results of STRING were imported into Cytoscape (version 3.7.2), and the hub genes was extracted by using the cytoHubba plugin. Finally, top five genes were selected with the highest scores (Degree algorithm) in the network as hub genes.
Mouse model
Animal experiments were approved by the Central South University, Changsha, Institutional Animal Care and Use committee. Animals were maintained a 12:12 h light/dark cycle with a temperature-controlled specific pathogen free (SPF) environment. 8-wk-old male mice were used in this study. For analyses in the Jak2V617F bone marrow transplantation (BMT) model, experiments were performed in strict adherence to China laws for animal welfare. This study was approved by the Ethics Committee of The Second Xiangya Hospital of Central South University (20220247).
For gene expression analysis in the conditional Jak2V617F knock-in model of PV, bone marrow (BM) from primary Vavcre Jak2V617F (Vavcre+ Jak2VF/+) mice (Tiedt et al., 2008; Mullally et al., 2010) (Jackson Lab (United States) was transplanted into lethally irradiated C57BL/6 mice (Slac Laboratory Animal Inc. Shanghai, China) (PV group), and the recipients progressed to PV before being used in the MEPs (megakaryocyte and erythroid progenitors) flow sorting experiment. As competitor transplantation, bone marrow from C57BL/6 recipient mice was transplanted into lethally irradiated C57BL/6 recipients (WT group).
6–8 weeks after BMT, BM cells of PV and WT recipients were harvested from bilateral femurs and tibias by crushing bones with mortar and pestle using PBS (Servicebio). Cells were filtered through 70 μm nylon mesh to obtain a single-cell suspension. Red blood cells were depleted by using erythrocyte lysis buffer (ACK buffer, Invitrogen).
For FACS analysis and cell sorting: EasySep™ Mouse Hematopoietic Progenitor Cell Isolation Kit (Catalog #19856, STEMCELL) was used to enrich the Lin− cells with biotinylated antibodies recognizing specific cell surface markers. The cell surface was stained using antibodies as shown in Supplementary Table S2: FITC-anti-CD4 (GK1.5, BD), FITC-anti-CD5 (53–7.3, Biolegend), FITC-anti-CD8 (53–6.72, Biolegend), FITC-anti-B220 (RA3-6B2, Biolegend), FITC-anti-TER-119 (Biolegend), FITC-Gr-1 (RB6-8C5, Biolegend), FITC-anti-Mac-1 (M1/70, Biolegend), PE-cy7-anti-Sca-1 (D7, BD), anti-PE-Kit (2B8, Biolegend), BV421-anti-CD34 (RAM34, BD), APC-Fcγ R III/II (93, Biolegend). Zombie aqua (Biolegend) was used to distinguish dead and alive cells. Live, singlet cells were selected for gating and cell sorting (Megakaryocyte Erythroid Progenitor cells, MEP: Lineage (CD4, CD5, CD8a, Mac1, B220,Ter119, Gr1)− Sca1+in− Sca-1− c-Kit+ CD34− Fcγ R III/II− Zombie aqua−) (Mead et al., 2013; Rao et al., 2019). 80,000—120,000 MEPs were isolated from PV and WT recipient mice by FACS sorting and used to further gene expression analysis. Cells were analyzed on a Fortessa Flow Cytometer and sorted on a FACSAria-II cell sorter (BD biosciences). Data was analyzed either using FACS Diva software or FlowJo (version 10.8.1) software (Treestar).
RNA extraction and quantitative real-time PCR
After extracting total RNA of flow-sorted MEP cells using TRIzol (Invitrogen, United States) according to the manufacturer’s protocol, reverse transcription high-capacity cDNA reverse transcription kit (AG, AG11728, Changsha, China) was performed. Lastly qRT-PCR using SYB Green Master quantitative PCR system (YEASEN, Shanghai, China) was conducted. Relative expression levels was calculated by using the 2-△△CT method. The sequences of the primers were listed in Supplementary Table S3.
Statistical analysis
R software (4.1.3) and GraphPad Prism 6.0 were used among all statistical analyses. Statistical significance was set at p-value <0.05 for all analyses.
Results
GEO information
The study flow chart was shown in Figure 1. The two GEO datasets numbered GSE136335 and GSE145802 were selected to study according to the criteria previously set. Information and characteristics of the three datasets were summarized in Table1. The GSE136335 was confirmed as a discovery cohort with GSE145802 paired as a validated cohort.
TABLE 1. Summary of GEO datasets involving MPN patients.
Identification of anoikis-related DEGs between MPNs and controls
1195 DEGs were identified by screening GSE136335, with 557 downregulated genes (Figure 2A, blue dot) and 638 upregulated genes (Figure 2A, red dot). 795 anoikis-related genes were downloaded from the “GeneCard” database (Supplementary Table S1). Veen diagram (Figure 2B); the left region (1,137 genes) refered to unique DEGs in GSE136335; the middle intersection region (58 genes) refers to anoikis-related DEGs; and the right region (736 genes) refered to unique anoikis genes. The heat map in Figure 2C showed the standardized expression of anoikis-related DEGs (36 upregulated and 22 downregulated). The PPI network (Figure 2D) containing 54 nodes and 144 edges was obtained. Figure 3 showed the relative expression profiles of 58 anoikis-related DEGs in the MPNs (Figure 3A) and MPN-subtypes (Figure 3B).
FIGURE 2. Screening of differentially expressed anoikis-related genes (ARGs) in MPN patients and healthy controls. (A) Volcano plot of genes differentially expressed between MPN patients and healthy controls in the GSE136335 dataset. Blue nodes represent downregulation in MPNs; red nodes represent upregulation; and grey nodes represent no significant difference from controls. (B) Intersection of differentially expressed genes (DEGs) in the GSE136335 and anoikis genes. The count on the left (1,137 genes) refers to DEGs unique to GSE136335; the count in the middle (58 genes) refers to anoikis-related DEGs; and the count on the right (736 genes) refers to unique to anoikis genes. (C) Heat map of 58 anoikis-related DEGs. (D) Protein–protein interaction (PPI) network of differentially expressed anoikis DEGs.
FIGURE 3. Expression of 58 anoikis-related DEGs in MPN (A) and subtype clusters (B).
In GO-BP analysis (Figure 4A), the major pathways were also apoptosis or oxidative stress pathways. In GO-CC analysis (Figure 4B), collagen-containing extracellular matrix was significantly enriched. The results of enrichment analysis in GO-MF (Figure 4C) were mainly the activation of cell adhesion-related pathway like cadherin binding. In the KEGG analysis (Figure 4D), the top 10 enriched pathways were mainly apoptosis, FoxO signaling pathway, TNF signaling pathway, etc.
FIGURE 4. The enrichment analysis of anoikis-related DEGs in MPNs. The dot plots display the Top10 terms enriched by anoikis-related DEGs in BPs (A), CCs (B), MFs (C) and KEGG pathways (D).
Identification of anoikis-related DEGs between PV and controls
925 DEGs were identified in PV patients compared with controls, with 513 downregulated genes and 412 upregulated genes (Figure 5A). Intersection analysis showed 36 anoikis-related DEGs (Figure 5B). The standardized expression of anoikis-related DEGs (17 upregulated and 19 downregulated) was shown in the heat map (Figure 5C). The results of enrichment analysis in GO and KEGG (Figure 5D; Supplementary Figure S1A) were also mainly focused on the apoptosis and activation of cell adhesion-related pathway like integrin binding. PPI network (Figure 6A) was constructed to select top 5 hub genes (CASP3, IL1B, HIF1A, CYCS, MCL1) (Figure 6B) according to the Degree algorithm. The expression levels of hub genes in training set were shown (Figure 6C).
FIGURE 5. Landscape of Anoikis-related Genes in Polycythemia Vera. (A)Volcano plot of genes differentially expressed between Polycythemia Vera (PV) patients and healthy controls in the GSE136335 dataset. Blue nodes represent downregulation in PV; red nodes represent upregulation; and grey nodes represent no significant difference from controls. (B) Intersection of differentially expressed genes (DEGs) in the GSE136335 and anoikis genes. The count on the left (889 genes) refers to DEGs unique to GSE136335; the count in the middle (36 genes) refers to anoikis-related DEGs; and the count on the right (758 genes) refers to unique to anoikis genes. (C) Heat map of 36 anoikis-related DEGs. (D) Top 10 GO biological process pathway.
FIGURE 6. The protein-protein interaction network of anoikis-related DEGs in Polycythemia Vera. (A) PPI network for 36 anoikis-related DEGs; (B) Subnetwork of hub genes from the PPI network. Node color reflects the degree of connectivity (red color represents a higher degree, and yellow color represents a lower degree). (C) The expression difference of five hub genes in the training set (GSE136335).
Validation analysis was conducted by using another dataset GSE145802. The expression levels of anoikis related DEGs in PV, CASP3, IL1B, HIF1A were significantly upregulated and CYCS was obviously downregulated, while the expression of MCL1 was downward trend in PV patients compared with normal controls (Figure 7A). However, the expression levels of anoikis related DEGs in PV-treated, only HIF1A was upregulate (Figure 7B).
FIGURE 7. Expression of CASP3, CYCS, HIF1A, IL1B and MCL1. The expression difference of CASP3, CYCS, HIF1A, IL1B and MCL1 in GSE145802 validation set in untreated-PV patients (PV) and healthy controls (HC) (A), and treated-PV patients (PV-treated) and HC (B).
Expression difference of hub genes in PV mice model
The expression of five PV hub genes was verified with our experimental data. PV mice model with Jak2V617F mutation (PV) and the corresponding WT mice (WT) were constructed by bone marrow transplantation (BMT) assay (Figure 8A). The spleen weight (Figure 8B) and hematocrit (HCT) (Figure 8C) of the recipient mice were measured 6–8 weeks after BMT, and the results indicated PV models successfully were constructed. PV and WT MEP cells were obtained by Mouse Hematopoietic Progenitor Cell Isolation Kit combined with FACS sorting (List of flow cytometry antibodies can be found in Supplementary Table S2) following the strict experimental procedures (Figure 8D). The relative mRNA expression (The sequences of the primers are listed in Supplementary Table S3) of Casp3 and Il1b were upregulated in MEPs of PV mice compared with WT mice (Figure 8E). However, other hub genes had no difference between these two group mice.
FIGURE 8. The expression difference of PV hub genes in PV mice model. (A) The workflow of WT and PV mice model constructed by bone marrow transplantation (BMT). (B) Spleen weight 6–8 weeks after BMT. Data are represented as mean ± SEM (***p < 0.001 versus WT group). (C) Hematocrit (HCT) 6–8 weeks after BMT. Data are represented as mean ± SEM (**p < 0.01 versus WT group). (D) The experimental procedure of obtaining WT and PV mice bone marrow megakaryocyte and erythroid progenitors (MEPs) by Lin-cell enrichment and flow sorting. (E) The validation of expression difference of five PV hub genes in WT and PV mice MEPs.
Discussion
Myeloproliferative neoplasms (MPNs) are myeloid malignancies characterized by clonal expansion of one or more differentiated myeloid lineages (Dameshek, 1951). Anoikis, a particular form of programmed cell death, was reported to perform important biological function in preventing attachment or grow of dysplastic cell for tissue homeostasis and development (Bakir et al., 2020; Li et al., 2021). Anoikis plays a crucial role in tumor invasion, migration, distant metastasis and drug resistance (Jin et al., 2017; Jiang et al., 2020; Kim et al., 2021).
In our study, a new view of differential expression anoikis-related gene sets in MPN patients and controls was provided. MPNs are the malignancies caused by dysfunctional HSC, thus CD34+ PB/BM MNCs was chosen for bioinformatics analysis. The intersection of MPN-DEGs and ARGs revealed 58 anoikis-related genes were specifically associated with MPNs. Gene set enrichment analysis showed cell adhesion and apoptosis-related pathways were fully enriched to activate the anoikis pathway. However, in MPNs, circulating cells with disease driver mutations (such as JAK2V617F), including red cells, leukocytes, platelets, as well as some vascular endothelial cells, were abnormal (Hasselbalch et al., 2021). When activated in MPNs, these abnormal cells interacted with each other and created proadhesive and prothrombotic environment in the circulation to predispose patients to thrombosis and occlusive complications (Hasselbalch et al., 2021).
PV is characterized by mutations in JAK2 exon 14 or 12 in hematopoietic stem cells, leading to erythrocytosis and systematic symptoms (Greenfield et al., 2021). PV mainly reduces the quality of life and shortens the life span of patients due to the risk of vascular events and transformation to myelofibrosis, myelodysplastic syndrome, and acute myeloid leukemia (Tefferi et al., 2013; Tiribelli et al., 2015; Frederiksen et al., 2019). Our results found that 36 DEGs in PV were associated with anoikis. Combined with the results of gene expression and gene enrichment, CD34+ mononuclear cells in PV deflected stronger adhesion effect than apoptosis, in accordance with the characteristics of PV prone to thromboembolism. Five hub genes (CASP3, IL1B, HIF1A, CYCS, MCL1) was selected with highest score through PPI network combined further analysis. The pro-apoptotic gene CASP3 was expressed higher in PV-CD34+ cells in the training and validation cohorts and in MEP cells in our PV mouse model, contrary to previous report (Steidl et al., 2005). The large variation may root in PV patients of different sources. The level of proinflammatory factor IL1B was upregulated in PV-CD34+ cells in training and validation cohorts and in MEP cells in our PV mouse model, consistent with previous report (Vaidya et al., 2012). Baumeister, Julian et al. reported higher protein level of HIF1α in 32D Jak2V617F-expressing cells than Jak2WT, and identified HIF-1 target genes involved in the Warburg effect as a possible underlying mechanism in JAK2V617F-positive MPN (Baumeister et al., 2020). The mRNA level of HIF1A was higher in PV-CD34+ cells in validation cohort. MCL1, a typical anti-apoptotic gene, plays a crucial role in the survival of HSCs and early progenitors (Opferman et al., 2005). Importantly, depletion of Mcl-1 was able to disrupt the viability of JAK2V617F mutant cells and sensitized mutant cells to JAK2 inhibitor therapy (Rubert et al., 2011). CYCS, as an apoptotic protein, is located between the inner and outer mitochondrial membranes (Ow et al., 2008; Shakeri et al., 2017) and is released into the cytoplasm in response to cell death signals and activated by pro-apoptotic Bcl2 family proteins (Wolter et al., 1997; Kennedy et al., 1999). The expressions of apoptotic-related genes MCL1 and CYCS decreased in the PV-CD34 + cells of the training and verification cohorts, while no difference was found in the PV mice MEP cells. CASP3 and IL1B were significantly increased in PV and decreased after treatment in validation cohort and expected to become indicators in PV development and treatment.
The inconformity between our results and previously reported may arise from following aspects: Firstly, the analyzed data was from peripheral blood or bone marrow mononuclear cells of MPN patients, excluding erythrocytes and granulocytes. Secondly, some vascular endothelial cells with driver mutations were also excluded. Thirdly, MPN in our study was a group of diseases rather than a single disease, and the degree of cell adhesion and apoptosis varies among subtypes. Finally, our results suggested that CD34+ mononuclear cells in bone marrow or peripheral blood contributed significantly less to thromboembolism than other types of cells.
Although MPNs are a myeloid proliferative disease, our research firstly revealed the critical role of anoikis-related genes in MPNs. CASP3, CYCS, HIF1A, IL1B, MCL1 are the anoikis-related hub genes in PV and CASP3 and IL1B may become important indicators of PV development and treatment. And our research may provide new insight into mechanisms of PV.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.
Ethics statement
The animal study was reviewed and approved by the Ethics Committee of The Second Xiangya Hospital of Central South University.
Author contributions
WA, LX, and TD designed the study. WA and LX collected the data. WA and WH completed the experimental validation, including bone marrow transplantation, flow sorting and RNA extraction and quantitative realtime PCR. LX analyzed the data and WA drafted the manuscript. LX, TD, WH, and YT made critical revisions to the manuscript for important intellectual content. TD, HP, and GZ provided guidance on statistical analysis and presentation of results. All authors approved the final manuscript and agreed to submit for publication.
Funding
The National Natural Science Foundation of China, Key Program (82130024), the Science and Technology Innovation Program of Hunan Province (2020RC4009) and the Project of Innovation-Driven Plan of Central South University (2020CX015).
Acknowledgments
We acknowledged Gene Expression Omnibus database for providing their platforms and contributors for uploading their meaningful datasets. We acknowledged contributors for uploading their valuable databases. We also acknowledged the funding of the National Natural Science Foundation of China, Key Program (82130024), the Science and Technology Innovation Program of Hunan Province (2020RC4009) and the Project of Innovation-Driven Plan of Central South University (2020CX015).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2023.1139351/full#supplementary-material
References
Arber, D. A., Orazi, A., Hasserjian, R., Thiele, J., Borowitz, M. J., Le Beau, M. M., et al. (2016). The 2016 revision to the World Health Organization classification of myeloid neoplasms and acute leukemia. Blood 127 (20), 2391–2405. doi:10.1182/blood-2016-03-643544
Ashburner, M., Ball, C. A., Blake, J. A., Botstein, D., Butler, H., Cherry, J. M., et al. (2000). Gene ontology: Tool for the unification of biology. The gene ontology consortium. Nat. Genet. 25 (1), 25–29. doi:10.1038/75556
Bakir, B., Chiarella, A. M., Pitarresi, J. R., and Rustgi, A. K. (2020). EMT, MET, plasticity, and tumor metastasis. Trends Cell Biol. 30 (10), 764–776. doi:10.1016/j.tcb.2020.07.003
Barbui, T., Finazzi, G., and Falanga, A. (2013). Myeloproliferative neoplasms and thrombosis. Blood 122 (13), 2176–2184. doi:10.1182/blood-2013-03-460154
Baumeister, J., Chatain, N., Hubrich, A., Maie, T., Costa, I. G., Denecke, B., et al. (2020). Hypoxia-inducible factor 1 (HIF-1) is a new therapeutic target in JAK2V617F-positive myeloproliferative neoplasms. Leukemia 34 (4), 1062–1074. doi:10.1038/s41375-019-0629-z
Baxter, E. J., Scott, L. M., Campbell, P. J., East, C., Fourouclas, N., Swanton, S., et al. (2005). Acquired mutation of the tyrosine kinase JAK2 in human myeloproliferative disorders. Lancet 365 (9464), 1054–1061. doi:10.1016/S0140-6736(05)71142-9
Berezovskaya, O., Schimmer, A. D., Glinskii, A. B., Pinilla, C., Hoffman, R. M., Reed, J. C., et al. (2005). Increased expression of apoptosis inhibitor protein XIAP contributes to anoikis resistance of circulating human prostate cancer metastasis precursor cells. Cancer Res. 65 (6), 2378–2386. doi:10.1158/0008-5472.CAN-04-2649
Brusson, M., De Grandis, M., Cochet, S., Bigot, S., Marin, M., Leduc, M., et al. (2018). Impact of hydroxycarbamide and interferon-alpha on red cell adhesion and membrane protein expression in polycythemia vera. Haematologica 103 (6), 972–981. doi:10.3324/haematol.2017.182303
Cervantes, F., Passamonti, F., and Barosi, G. (2008). Life expectancy and prognostic factors in the classic BCR/ABL-negative myeloproliferative disorders. Leukemia 22 (5), 905–914. doi:10.1038/leu.2008.72
Dameshek, W. (1951). Some speculations on the myeloproliferative syndromes [editorial]. BloodBlood 6127 (46), 372–375.
Frederiksen, H., Szepligeti, S., Bak, M., Ghanima, W., Hasselbalch, H. C., and Christiansen, C. F. (2019). Vascular diseases in patients with chronic myeloproliferative neoplasms - impact of comorbidity. Clin. Epidemiol. 11, 955–967. doi:10.2147/CLEP.S216787
Gentleman, R. Bioinformatics and computational biology solutions using R and Bioconductor. New York: Springer Science+Business Media; 2005. xix, 473 p. p.
Greenfield, G., McMullin, M. F., and Mills, K. (2021). Molecular pathogenesis of the myeloproliferative neoplasms. J. Hematol. Oncol. 14 (1), 103. doi:10.1186/s13045-021-01116-z
Han, H. J., Sung, J. Y., Kim, S. H., Yun, U. J., Kim, H., Jang, E. J., et al. (2021). Fibronectin regulates anoikis resistance via cell aggregate formation. Cancer Lett. 508, 59–72. doi:10.1016/j.canlet.2021.03.011
Hasselbalch, H. C., Elvers, M., and Schafer, A. I. (2021). The pathobiology of thrombosis, microvascular disease, and hemorrhage in the myeloproliferative neoplasms. Blood 137 (16), 2152–2160. doi:10.1182/blood.2020008109
Ito, K., and Murphy, D. (2013). Application of ggplot2 to pharmacometric graphics. CPT Pharmacometrics Syst. Pharmacol. 2 (10), e79. doi:10.1038/psp.2013.56
James, C., Ugo, V., Le Couedic, J. P., Staerk, J., Delhommeau, F., Lacout, C., et al. (2005). A unique clonal JAK2 mutation leading to constitutive signalling causes polycythaemia vera. Nature 434 (7037), 1144–1148. doi:10.1038/nature03546
Jiang, K., Yao, G., Hu, L., Yan, Y., Liu, J., Shi, J., et al. (2020). MOB2 suppresses GBM cell migration and invasion via regulation of FAK/Akt and cAMP/PKA signaling. Cell Death Dis. 11 (4), 230. doi:10.1038/s41419-020-2381-8
Jin, L., Chun, J., Pan, C., Alesi, G. N., Li, D., Magliocca, K. R., et al. (2017). Phosphorylation-mediated activation of LDHA promotes cancer cell invasion and tumour metastasis. Oncogene 36 (27), 3797–3806. doi:10.1038/onc.2017.6
Kakavandi, E., Shahbahrami, R., Goudarzi, H., Eslami, G., and Faghihloo, E. (2018). Anoikis resistance and oncoviruses. J. Cell Biochem. 119 (3), 2484–2491. doi:10.1002/jcb.26363
Kanehisa, M., and Goto, S. (2000). Kegg: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 28 (1), 27–30. doi:10.1093/nar/28.1.27
Kc, D., Falchi, L., and Verstovsek, S. (2017). The underappreciated risk of thrombosis and bleeding in patients with myelofibrosis: A review. Ann. Hematol. 96 (10), 1595–1604. doi:10.1007/s00277-017-3099-2
Kennedy, S. G., Kandel, E. S., Cross, T. K., and Hay, N. (1999). Akt/Protein kinase B inhibits cell death by preventing the release of cytochrome c from mitochondria. Mol. Cell Biol. 19 (8), 5800–5810. doi:10.1128/MCB.19.8.5800
Kim, H., Choi, P., Kim, T., Kim, Y., Song, B. G., Park, Y. T., et al. (2021). Ginsenosides Rk1 and Rg5 inhibit transforming growth factor-β1-induced epithelial-mesenchymal transition and suppress migration, invasion, anoikis resistance, and development of stem-like features in lung cancer. J. Ginseng Res. 45 (1), 134–148. doi:10.1016/j.jgr.2020.02.005
Kim, Y. N., Koo, K. H., Sung, J. Y., Yun, U. J., and Kim, H. (2012). Anoikis resistance: An essential prerequisite for tumor metastasis. Int. J. Cell Biol. 2012, 306879. doi:10.1155/2012/306879
Kralovics, R., Passamonti, F., Buser, A. S., Teo, S. S., Tiedt, R., Passweg, J. R., et al. (2005). A gain-of-function mutation of JAK2 in myeloproliferative disorders. N. Engl. J. Med. 352 (17), 1779–1790. doi:10.1056/NEJMoa051113
Levine, R. L., Wadleigh, M., Cools, J., Ebert, B. L., Wernig, G., Huntly, B. J., et al. (2005). Activating mutation in the tyrosine kinase JAK2 in polycythemia vera, essential thrombocythemia, and myeloid metaplasia with myelofibrosis. Cancer Cell 7 (4), 387–397. doi:10.1016/j.ccr.2005.03.023
Li, S., Chen, Y., Zhang, Y., Jiang, X., Jiang, Y., Qin, X., et al. (2019). Shear stress promotes anoikis resistance of cancer cells via caveolin-1-dependent extrinsic and intrinsic apoptotic pathways. J. Cell Physiol. 234 (4), 3730–3743. doi:10.1002/jcp.27149
Li, S., Zhang, J., Qian, S., Wu, X., Sun, L., Ling, T., et al. (2021). S100A8 promotes epithelial-mesenchymal transition and metastasis under TGF-β/USF2 axis in colorectal cancer. Cancer Commun. (Lond) 41 (2), 154–170. doi:10.1002/cac2.12130
Mead, A. J., Kharazi, S., Atkinson, D., Macaulay, I., Pecquet, C., Loughran, S., et al. (2013). FLT3-ITDs instruct a myeloid differentiation and transformation bias in lymphomyeloid multipotent progenitors. Cell Rep. 3 (6), 1766–1776. doi:10.1016/j.celrep.2013.04.031
Mead, A. J., and Mullally, A. (2017). Myeloproliferative neoplasm stem cells. Blood 129 (12), 1607–1616. doi:10.1182/blood-2016-10-696005
Mullally, A., Lane, S. W., Ball, B., Megerdichian, C., Okabe, R., Al-Shahrour, F., et al. (2010). Physiological Jak2V617F expression causes a lethal myeloproliferative neoplasm with differential effects on hematopoietic stem and progenitor cells. Cancer Cell 17 (6), 584–596. doi:10.1016/j.ccr.2010.05.015
Opferman, J. T., Iwasaki, H., Ong, C. C., Suh, H., Mizuno, S., Akashi, K., et al. (2005). Obligate role of anti-apoptotic MCL-1 in the survival of hematopoietic stem cells. Science 307 (5712), 1101–1104. doi:10.1126/science.1106114
Ow, Y. P., Green, D. R., Hao, Z., and Mak, T. W. (2008). Cytochrome c: Functions beyond respiration. Nat. Rev. Mol. Cell Biol. 9 (7), 532–542. doi:10.1038/nrm2434
Passamonti, F., Griesshammer, M., Palandri, F., Egyed, M., Benevolo, G., Devos, T., et al. (2017). Ruxolitinib for the treatment of inadequately controlled polycythaemia vera without splenomegaly (RESPONSE-2): A randomised, open-label, phase 3b study. Lancet Oncol. 18 (1), 88–99. doi:10.1016/S1470-2045(16)30558-7
Rao, T. N., Hansen, N., Hilfiker, J., Rai, S., Majewska, J. M., Lekovic, D., et al. (2019). JAK2-mutant hematopoietic cells display metabolic alterations that can be targeted to treat myeloproliferative neoplasms. Blood 134 (21), 1832–1846. doi:10.1182/blood.2019000162
Rubert, J., Qian, Z., Andraos, R., Guthy, D. A., and Radimerski, T. (2011). Bim and Mcl-1 exert key roles in regulating JAK2V617F cell survival. BMC Cancer 11, 24. doi:10.1186/1471-2407-11-24
Scott, L. M., Tong, W., Levine, R. L., Scott, M. A., Beer, P. A., Stratton, M. R., et al. (2007). JAK2 exon 12 mutations in polycythemia vera and idiopathic erythrocytosis. N. Engl. J. Med. 356 (5), 459–468. doi:10.1056/NEJMoa065202
Shakeri, R., Kheirollahi, A., and Davoodi, J. (2017). Apaf-1: Regulation and function in cell death. Biochimie 135, 111–125. doi:10.1016/j.biochi.2017.02.001
Spivak, J. L. (2017). Myeloproliferative neoplasms. N. Engl. J. Med. 376 (22), 2168–2181. doi:10.1056/NEJMra1406186
Steidl, U., Schroeder, T., Steidl, C., Kobbe, G., Graef, T., Bork, S., et al. (2005). Distinct gene expression pattern of malignant hematopoietic stem and progenitor cells in polycythemia vera. Ann. N. Y. Acad. Sci. 1044, 94–108. doi:10.1196/annals.1349.013
Szekely, G. J., and Rizzo, M. L. (2005). Hierarchical clustering via joint between-within distances: Extending ward's minimum variance method. J. Classif. 22 (2), 151–183. doi:10.1007/s00357-005-0012-9
Tefferi, A., Rumi, E., Finazzi, G., Gisslinger, H., Vannucchi, A. M., Rodeghiero, F., et al. (2013). Survival and prognosis among 1545 patients with contemporary polycythemia vera: An international study. Leukemia 27 (9), 1874–1881. doi:10.1038/leu.2013.163
Tiedt, R., Hao-Shen, H., Sobas, M. A., Looser, R., Dirnhofer, S., Schwaller, J., et al. (2008). Ratio of mutant JAK2-V617F to wild-type Jak2 determines the MPD phenotypes in transgenic mice. Blood 111 (8), 3931–3940. doi:10.1182/blood-2007-08-107748
Tiribelli, M., Barraco, D., De Marchi, F., Marin, L., Medeot, M., Damiani, D., et al. (2015). Clinical factors predictive of myelofibrotic evolution in patients with polycythemia vera. Ann. Hematol. 94 (5), 873–874. doi:10.1007/s00277-014-2257-z
Vaidya, R., Gangat, N., Jimma, T., Finke, C. M., Lasho, T. L., Pardanani, A., et al. (2012). Plasma cytokines in polycythemia vera: Phenotypic correlates, prognostic relevance, and comparison with myelofibrosis. Am. J. Hematol. 87 (11), 1003–1005. doi:10.1002/ajh.23295
Vannucchi, A. M., Kiladjian, J. J., Griesshammer, M., Masszi, T., Durrant, S., Passamonti, F., et al. (2015). Ruxolitinib versus standard therapy for the treatment of polycythemia vera. N. Engl. J. Med. 372 (5), 1670–1671. doi:10.1056/NEJMc1502524
Verstovsek, S., Mesa, R. A., Gotlib, J., Levy, R. S., Gupta, V., DiPersio, J. F., et al. (2012). A double-blind, placebo-controlled trial of ruxolitinib for myelofibrosis. N. Engl. J. Med. 366 (9), 799–807. doi:10.1056/NEJMoa1110557
Wang, X., Ye, F., Tripodi, J., Hu, C. S., Qiu, J., Najfeld, V., et al. (2014). JAK2 inhibitors do not affect stem cells present in the spleens of patients with myelofibrosis. Blood 124 (19), 2987–2995. doi:10.1182/blood-2014-02-558015
Wolter, K. G., Hsu, Y. T., Smith, C. L., Nechushtan, A., Xi, X. G., and Youle, R. J. (1997). Movement of Bax from the cytosol to mitochondria during apoptosis. J. Cell Biol. 139 (5), 1281–1292. doi:10.1083/jcb.139.5.1281
Wu, T., Hu, E., Xu, S., Chen, M., Guo, P., Dai, Z., et al. (2021). clusterProfiler 4.0: A universal enrichment tool for interpreting omics data. Innov. (Camb) 2 (3), 100141. doi:10.1016/j.xinn.2021.100141
Xiao, T., Xu, Z., Zhou, Y., Zhang, H., Geng, J., Liang, Y., et al. (2019). Loss of TP53I11 enhances the extracellular matrix-independent survival by promoting activation of AMPK. IUBMB Life 71 (2), 183–191. doi:10.1002/iub.1949
Yu, Y., Song, Y., Cheng, L., Chen, L., Liu, B., Lu, D., et al. (2022). CircCEMIP promotes anoikis-resistance by enhancing protective autophagy in prostate cancer cells. J. Exp. Clin. Cancer Res. 41 (1), 188. doi:10.1186/s13046-022-02381-7
Zhi, Z., Ouyang, Z., Ren, Y., Cheng, Y., Liu, P., Wen, Y., et al. (2022). Non-canonical phosphorylation of Bmf by p38 MAPK promotes its apoptotic activity in anoikis. Cell Death Differ. 29 (2), 323–336. doi:10.1038/s41418-021-00855-3
Keywords: myeloproliferative neoplasms, anoikis-related genes, differentially expressed genes, protein-protein interaction, hub genes
Citation: Aini W, Xie L, Hu W, Tang Y, Peng H, Zhang G and Deng T (2023) Exploration and identification of anoikis-related genes in polycythemia vera. Front. Genet. 14:1139351. doi: 10.3389/fgene.2023.1139351
Received: 06 January 2023; Accepted: 09 February 2023;
Published: 17 February 2023.
Edited by:
Hang Song, Anhui University of Chinese Medicine, ChinaReviewed by:
Jian-ping Hao, First Affiliated Hospital of Xinjiang Medical University, ChinaYanqiang Li, Boston Children’s Hospital, Harvard Medical School, United States
Jie Lv, St. Jude Children’s Research Hospital, United States
Copyright © 2023 Aini, Xie, Hu, Tang, Peng, Zhang and Deng. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Tuo Deng, dengtuo@csu.edu.cn
†These authors have contributed equally to this work