Identification of a Nervous System Gene Expression Signature in Colon Cancer
Total Page:16
File Type:pdf, Size:1020Kb
bioRxiv preprint doi: https://doi.org/10.1101/2021.02.02.428317; this version posted February 3, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-NC-ND 4.0 International license. Regan et al. 2020 Neural Crest Gene Regulation of CSCs 1 Identification of a Nervous System Gene Expression Signature in Colon Cancer 2 Stem Cells Reveals a Role for Neural Crest Regulators EGR2 and SOX2 in 3 Tumorigenesis 4 5 Joseph L. Regan1,2,13*, Dirk Schumacher3,4, Stephanie Staudte1,2, Andreas Steffen1, Ralf 6 Lesche1,5, Joern Toedling1,5, Thibaud Jourdan1, Johannes Haybaeck6,7, Nicole Golob- 7 Schwarzl7,8, Dominik Mumberg1, David Henderson1, Balázs Győrffy9,10, Christian R.A. 8 Regenbrecht3,11,12, Ulrich Keilholz2, Reinhold Schäfer2,3,4, Martin Lange1,5 9 10 Affiliations 11 1Bayer AG, Drug Discovery, Pharmaceuticals, 13342 Berlin, Germany 12 2Charité Comprehensive Cancer Center, Charité - Universitätsmedizin Berlin, 10117 13 Berlin, Germany 14 3Laboratory of Molecular Tumor Pathology, Charité Universitätsmedizin Berlin, 10117 15 Berlin, Germany 16 4German Cancer Consortium (DKTK), DKFZ, 69120 Heidelberg, Germany 17 5Nuvisan ICB GmbH, 13353 Berlin, Germany 18 6Department of Pathology, Neuropathology and Molecular Pathology, Medical 19 University of Innsbruck, Austria 20 7Diagnostic & Research Center for Molecular Biomedicine, Institute of Pathology, 21 Medical University of Graz, Austria 22 8Department of Dermatology and Venereology, Medical University of Graz, Austria 23 9Department of Bioinformatics, Semmelweis University, 1094 Budapest, Hungary 1 bioRxiv preprint doi: https://doi.org/10.1101/2021.02.02.428317; this version posted February 3, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-NC-ND 4.0 International license. Regan et al. 2020 Neural Crest Gene Regulation of CSCs 24 10TTK Cancer Biomarker Research Group, Institute of Enzymology, 1117 Budapest, 25 Hungary 26 11CELLphenomics GmbH, 13125 Berlin, Germany 27 12Institute of Pathology, University Medical Center Göttingen, 37075 Göttingen, 28 Germany 29 13Lead Contact 30 *Correspondence: [email protected] (J.L.R.) 31 32 SUMMARY 33 Recent data support a hierarchical model of colon cancer driven by a population of cancer 34 stem cells (CSCs). Greater understanding of the mechanisms that regulate CSCs may 35 therefore lead to more effective treatments. Serial limiting dilution xenotransplantation 36 assays of colon cancer patient-derived tumors demonstrated ALDHPositive cells to be 37 enriched for tumorigenic self-renewing CSCs. In order to identify CSC modulators, we 38 performed RNA-sequencing analysis of ALDHPositive CSCs from a panel of colon cancer 39 patient-derived organoids (PDOs) and xenografts (PDXs). These studies demonstrated 40 CSCs to be enriched for embryonic and neural development gene sets. Functional 41 analyses of genes differentially expressed in both ALDHPositive PDO and PDX CSCs 42 demonstrated the neural crest stem cell (NCSC) regulator and wound response gene 43 EGR2 to be required for CSC tumorigenicity and to control expression of homeobox 44 superfamily embryonic master transcriptional regulator HOX genes and the embryonic 45 and neural stem cell regulator SOX2. In addition, we identify EGR2, HOXA2, HOXA4, 2 bioRxiv preprint doi: https://doi.org/10.1101/2021.02.02.428317; this version posted February 3, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-NC-ND 4.0 International license. Regan et al. 2020 Neural Crest Gene Regulation of CSCs 46 HOXA5, HOXA7, HOXB2, HOXB3 and the tumor suppressor ATOH1 as new prognostic 47 biomarkers in colorectal cancer. 48 49 INTRODUCTION 50 Colorectal cancer (CRC), the third most common cancer and fourth most common cause 51 of cancer deaths worldwide1, is a heterogeneous tumor driven by a subpopulation of 52 CSCs, that may also be the source of relapse following treatment2–5. Elucidation of the 53 mechanisms that regulate CSC survival and tumorgenicity may therefore lead to novel 54 treatments and improved patient outcomes. 55 56 CSCs are undifferentiated cancer cells that share many of the attributes of stem cells, 57 such as multipotency, self-renewal and the ability to produce daughter cells that 58 differentiate2,6,7. Stem cells are controlled by core gene networks that include the 59 embryonic master transcriptional regulator HOX genes8,9 and SOX210, whose 60 misregulation can result in aberrant stem cell function, developmental defects and 61 cancer11,12. These genes are crucial for embryonic development and their expression is 62 maintained in adult tissue stem cells, where they regulate self-renewal and 63 differentiation9,13–15. HOX genes and SOX2 are aberrantly expressed in several cancers, 64 including CRC, and emerging evidence demonstrates their involvement in the 65 transformation of tissue stem cells into CSCs11,16–23. Modulation of HOX genes and SOX2 66 could therefore provide novel therapeutic strategies to block tumorigenesis and overcome 67 therapy resistance in CRC and other CSC driven cancers. 68 3 bioRxiv preprint doi: https://doi.org/10.1101/2021.02.02.428317; this version posted February 3, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-NC-ND 4.0 International license. Regan et al. 2020 Neural Crest Gene Regulation of CSCs 69 During embryonic development of the neural crest, which gives rise to the peripheral 70 nervous system (PNS) and several non-neuronal cell types24, HOX and SOX genes are 71 regulated by retinoic acid25,26, a product of the normal tissue stem cell and CSC marker 72 aldehyde dehydrogenase (ALDH1A1, ALDH1A2, ALDH1A3)8,26–29, and by the neural 73 crest stem cell (NCSC) zinc finger transcription factor and wound response gene EGR2 74 (KROX20)30–41. 75 76 Here we carried out whole transcriptome analysis of functionally tested ALDHPositive CSCs 77 from a panel of colon PDOs and PDX models and show that colon CSCs and Lgr5Positive 78 intestinal stem cells (ISCs) are highly enriched for nervous system development and 79 neural crest genes. Furthermore, we demonstrate that the neural crest stem cell (NCSC) 80 gene EGR2 is a marker of poor prognosis in CRC and modulates expression of HOX 81 genes and SOX2 in CSCs to regulate tumorigenicity and differentiation. 82 83 RESULTS 84 Colon cancer PDOs are heterogeneous and enriched for ALDHPositive self-renewing 85 CSCs 86 Colon cancer PDO models were established from freshly isolated primary tumors and 87 metastases from colon cancer patients (Table S1) by embedding in growth-factor reduced 88 Matrigel and cultivating in serum free media, as previously described42–44. 89 Immunostaining of PDOs for the structural proteins EZRIN and EPCAM demonstrated 90 that PDOs retain the apical-basal polarity and structural adhesion of the normal intestine 91 (Figure 1A). Immunostaining of PDOs and equivalent PDX models for stem cell regulator 4 bioRxiv preprint doi: https://doi.org/10.1101/2021.02.02.428317; this version posted February 3, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-NC-ND 4.0 International license. Regan et al. 2020 Neural Crest Gene Regulation of CSCs 92 Wnt signaling protein BETA-CATENIN demonstrated differences in nuclear localization 93 of BETA-CATENIN and confirmed previous data demonstrating heterogeneous Wnt 94 signaling activity within the tumors43 (Figure 1B). Increased aldehyde dehydrogenase 95 (ALDH) activity, as measured using the Aldefluor™ assay, is a marker of CSCs in colon 96 cancer and many other cancer types29. We previously carried out limiting dilution serial 97 xenotransplantation of ALDHNegative and ALDHPositive cells and demonstrated that colon 98 CSCs are ALDHPositive and enriched for Wnt signaling activity43. However, ALDHNegative 99 cells also gave rise to tumors when transplanted at higher cell numbers. In order to 100 determine if ALDHNegative and ALDHPositive cells maintained their self-renewal and 101 tumorigenic capacity, we performed additional rounds of limiting dilution serial 102 xenotransplantation of ALDHNegative and ALDHPositive cells (Figure 1E). These data 103 confirmed that PDOs are enriched for ALDHPositive cells compared to equivalent PDX 104 models (Figure 1C and D) and that ALDHPositive CSCs self-renew to maintain their 105 tumorigenic capacity over extended rounds of xenotransplantation, but that ALDHNegative 106 cells do not (Figure 1E). 107 108 Colon CSCs are enriched for embryonic and nervous system development gene 109 expression signatures 110 In order to identify modulators of colon CSCs, ALDHNegative cells and ALDHPositive CSCs 111 were isolated from PDO and PDX models and subjected to whole transcriptome analysis 112 by RNA-sequencing. ALDH1A1 is a marker of poor prognosis in several cancer 113 types27,29,45–49 and has been reported to be responsible for the aldehyde dehydrogenase 114 activity that defines the ALDHPositive cell fraction in the AldefluorTM assay50. However, 5 bioRxiv preprint doi: https://doi.org/10.1101/2021.02.02.428317; this version posted February 3, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-NC-ND 4.0 International license. Regan et al. 2020 Neural Crest Gene Regulation of CSCs 115 nineteen different isoforms of ALDH exist and several of these, including ALDH1A2, 116 ALDH1A3 and ALDH2 have also been reported to be involved in the AldefluorTM assay51– 117 53.