| Literature DB >> 30718226 |
Harrison M Penrose1, Chloe Cable1, Sandra Heller1, Nathan Ungerleider1, Hani Nakhoul1, Melody Baddoo1, Alifiani B Hartono1, Sean B Lee1, Matthew E Burow2, Erik F Flemington1, Susan E Crawford3, Suzana D Savkovic4.
Abstract
BACKGROUND & AIMS: Diminished forkhead box O3 (FOXO3) function drives inflammation and cancer growth; however, mechanisms fostering these pathobiologies are unclear. Here, we aimed to identify in colon loss of FOXO3-dependent cellular and molecular changes that facilitate inflammation-mediated tumor growth.Entities:
Keywords: Colon Cancer; FOXO3; Inflammation; RNAseq
Mesh:
Substances:
Year: 2018 PMID: 30718226 PMCID: PMC6360252 DOI: 10.1016/j.jcmgh.2018.10.003
Source DB: PubMed Journal: Cell Mol Gastroenterol Hepatol ISSN: 2352-345X
Figure 1In human IBD and cancer, FOXO3 mRNA levels are reduced and associated with poor survival and partial PI3K activating mutations. (A) In human IBD tissue samples, relative expression (microarray) of FOXO3 transcript was significantly reduced compared with normal matched control tissue (GSE4183, n = 23; Affymetrix probe ID: 204131_s_at; *P < .05, IBD compared with normal colon, Student t test). (B) mRNA expression (RNAseq, RSEM) of FOXO3 across TCGA cohorts (TIMER, https://cistrome.shinyapps.io/timer/). (C) Representative Kaplan-Meier survival plot demonstrating decreased FOXO3 mRNA expression is associated with reduced colon cancer patient survival (GSE16125, n = 32; https://precog.stanford.edu/index.php). (D) PRECOG meta-z analysis indicating increased FOXO3 expression as prognostic for favorable overall cancer survival (P < .05, https://precog.stanford.edu/index.php). (E) In human colon cancer (TCGA), the 10 patients expressing the lowest FOXO3 mRNA (compared with average FOXO3 level in tumor) (RNAseq, HTseq) and their 3 (p.H1047R, p.H1048R, p.R88Q) corresponding PI3K activating mutations (MuTect2 Variant Browser) (UCSC Xena Genome Browser, http://xena.ucsc.edu).
Figure 2In mouse colon FOXO3 deficiency increases colonic tumor burden and activation of inflammatory and tumorigenic transcriptional programs. (A and B) Gross histology of WT and FOXO3 KO colonic tumors in the AOM/DSS mouse model. Colonic tumor incidence, size, and tumor burden were determined between WT and FOXO3 KO (n = 8 per experimental group, representative graph from 2 independent experiments; *P < .05 compared with WT, Student t test). (C) Venn diagrams depict the number of transcripts commonly or uniquely expressed in mouse FOXO3 KO colon and tumors (WT or KO in AOM/DSS model) relative to WT (n = 3 per experimental group, FDR < 0.05, EBseq). (D) Disease and function analysis of colonic FOXO3 KO transcriptome revealed association with inflammatory and cancer pathobiology (n = 3 per experimental group, FDR < 0.05, IPA). (E) Top canonical pathways activated in mouse FOXO3 KO colon (n = 6) compared with both mouse inflamed and dysplastic colonic tissue (GSE31106, n = 15) (FDR < 0.05, IPA).
Figure 3Loss of FOXO3 leads to altered transcript expression similarly observed in human intestinal inflammation and colon cancer. (A) Canonical pathway analysis revealed high degree of transcriptional similarity between mouse FOXO3 KO colon (n = 6) and human IBD (GSE4183, n = 23 patient samples) (FDR < 0.05, IPA). (B) Pathway analysis revealed that transcriptional signatures from mouse FOXO3 KO colon and AOM/DSS tumors shared high degree of similarity with human colon cancer, especially those with MSI (TCGA) (n = 3 per experimental group, FDR < 0.05, IPA).
Figure 4FOXO3 deficiency promotes infiltration of immune cells in mouse colon and tumor microenvironment. (A and B) CIBERSORT relative abundances of select immune cell infiltrate in FOXO3 KO colon relative to WT (fold-change) (n = 3 per experimental group; *P < .05 compared with WT, analysis of variance and Student-Newman-Keuls post-test). (C) IPA upstream regulator analysis of FOXO3 KO colon and tumor (WT and KO) transcriptomes (n = 3 per experimental group, FDR < 0.05, IPA).
Figure 5In mouse colon FOXO3 deficiency leads to activation of transcriptional programs associated with inflammation, tumorigenesis, and bacterial signaling. (A) GSEA of FOXO3 KO colon transcriptome revealed activation of inflammatory and tumorigenic transcriptional motifs (n = 6, *P < .05). (B) Heatmap profiling revealed transcripts significantly elevated in FOXO3 KO colon and tumors (WT and KO) from inflammatory and cancer signaling (n = 3 per experimental group, FDR < 0.05, IPA). (C) Total RNA was extracted from WT or FOXO3 KO colon, and TLR4 mRNA levels were quantified by using qPCR (n = 4 per group from 2 independent experiments; *P < .05 compared with WT colon, Student t test).
DE Transcripts in FOXO3 KO Colon Similarly Up- and Down-Regulated in Tumors (AOM/DSS)
| Gene | Gene name | Colon FOXO3 KO | Tumor (AOM/DSS) |
|---|---|---|---|
| SGSM3 | Small G protein signaling modulator 3 | 595.078 | 91.142 |
| FAM103A1 | Family with sequence similarity 103, member A1 | 565.122 | 152.628 |
| ZBED6 | Zinc finger, BED type containing 6 | 106.973 | 61.437 |
| MCPT4 | Mast cell protease 4 | 41.723 | 19.823 |
| CFD | Complement factor D | 33.03 | 7.816 |
| ADIPOQ | Adiponectin, C1Q and collagen domain | 30.559 | 10.799 |
| LEP | Leptin | 25.853 | 16.14 |
| CPA3 | Carboxypeptidase A3, mast cell | 25.089 | 8.962 |
| ARFIP1 | ADP ribosylation factor interacting protein 1 | 20.479 | 3.18 |
| S100A8 | S100 calcium binding protein A8 | 15.034 | 792.785 |
| AMD1 | S-adenosylmethionine decarboxylase 1 | 14.95 | 2.282 |
| RBM12 | RNA binding protein 12 | 12.773 | 11.412 |
| MOAP1 | Modulator of apoptosis 1 | 12.399 | 4.448 |
| CHSY3 | Chondroitin sulfate synthase 3 | 10.147 | 4.858 |
| KLF12 | Kruppel like factor 12 | 7.802 | 6.82 |
| WNT16 | Wnt family member 16 | 6.814 | 37.196 |
| UPRT | Uracil phosphoribosyltransferase | 6.174 | 2.941 |
| LCOR | Ligand dependent nuclear receptor Corepressor | 6.089 | 3.089 |
| SYNE3 | Spectrin repeat containing, nuclear envelope family member 3 | 5.691 | 1.515 |
| PAPLN | Papilin, proteoglycan-like sulfated glycoprotein | 5.62 | 3.152 |
| CBL | Cbl proto-oncogene | 5.582 | 2.108 |
| LMBRD2 | LMBR1 domain containing 2 | 5.479 | 2.988 |
| LNPEP | Leucyl/cystinyl aminopeptidase | 5.468 | 2.236 |
| NT5C1A | 5'-nucleotidase, cytosolic IA | 5.194 | 3.413 |
| ATP7A | ATPase, copper transporting alpha | 5.098 | 3.102 |
| CCDC85C | Coiled-coil domain containing 85C | 4.92 | 2.431 |
| PCOLCE2 | Procollagen C-endopeptidase enhancer 2 | 4.846 | 1.806 |
| SLFN5 | Schlafen 5 | 4.763 | 3.99 |
| SAMD12 | Sterile alpha motif domain containing 12 | 4.647 | 5.439 |
| ITGA2 | Integrin alpha 2 | 4.356 | 20.059 |
| EXPH5 | Exophilin 5 | 4.348 | 2.457 |
| PM20D2 | Peptidase M20 domain containing 2 | 4.105 | 6.56 |
| PTAR1 | Protein prenyltransferase alpha subunit repeat containing 1 | 4.006 | 2.109 |
| VPS13C | Vacuolar protein sorting 13C | 4 | 1.736 |
| ADAMTS12 | ADAM metallopeptidase with thrombospondin type 1 motif 12 | 3.927 | 16.032 |
| TCHH | Trichohyalin | 3.897 | 5.262 |
| XKR4 | XK related 4 | 3.895 | 3.658 |
| MIA2 | Melanoma inhibitory activity 2 | 3.847 | 2.227 |
| PDPR | Pyruvate dehydrogenase phosphatase regulatory subunit | 3.784 | 1.999 |
| PCDHGC3 | Protocadherin gamma subfamily C, 3 | 3.771 | 3.844 |
| IGFBP3 | Insulin-like growth factor binding protein 3 | 3.711 | 10.152 |
| STON2 | Stonin 2 | 3.64 | 1.98 |
| ADAMTS3 | ADAM metallopeptidase with thrombospondin type 1 motif 3 | 3.574 | 2.979 |
| ARSB | Arylsulfatase B | 3.568 | 1.626 |
| TNKS | Tankyrase | 3.502 | 2.947 |
| SERPINA3 | Serpin family A member 3 | 3.48 | 7.076 |
| SLC25A30 | Solute carrier family 25, member 30 | 3.436 | 4.609 |
| SNX29 | Sorting nexin 29 | 3.356 | 1.737 |
| WDFY2 | WD repeat and FYVE domain containing 2 | 3.308 | 3.602 |
| LPAR5 | Lysophosphatidic acid receptor 5 | 3.295 | 2.073 |
| LENEP | Lens epithelial protein | –85.425 | –69.786 |
| USP9Y | Ubiquitin specific peptidase 9, y-linked | –42.229 | –22.166 |
| DDX54 | DEAD box polypeptide 54 | –21.688 | –2.797 |
| UGT1A4 | UDP glucuronosyltransferase 1 family, polypeptide A4 | –21.403 | –14.754 |
| LIPF | Lipase F | –9.735 | –5.841 |
| RPGRIP1 | Retinitis pigmentosa GTPase regulator Interacting protein 1 | –5.034 | –3.058 |
| RPL41 | Ribosomal protein L41 | –4.189 | –2.691 |
| CES3 | Carboxylesterase 3 | –3.743 | –6.789 |
| RAB33A | RAB33A, member RAS oncogene family | –3.669 | –4.316 |
| RGS1 | Regulator of G-protein signaling 1 | –3.665 | –5.723 |
| NTAN1 | N-terminal Asn amidase | –3.566 | –2.115 |
| UGT1A1 | UDP glucuronosyltransferase 1 family, polypeptide A1 | –3.531 | –2.425 |
| UGT2B17 | UDP glucuronosyltransferase family 2, member B17 | –3.505 | –3.248 |
| MFSD9 | Major facilitator superfamily domain containing 9 | –3.397 | –2.942 |
| SPDYA | Speedy/RINGO cell cycle regulator family, member A | –3.128 | –1.815 |
| SLC17A1 | Solute carrier family 17, member 1 | –2.86 | –2.822 |
| CFAP52 | Cilia and flagella associated protein 52 | –2.763 | –1.652 |
| DPPA5 | Developmental pluripotency associated 5 | –2.698 | –3.306 |
| NR1I3 | Nuclear receptor subfamily 1 group I member 3 | –2.674 | –6.541 |
| SERPINC1 | Serpin family C member 1 | –2.57 | –1.717 |
| CDA | Cytidine deaminase | –2.547 | –1.518 |
| TAT | Tyrosine aminotransferase | –2.402 | –3.944 |
| NLRP9 | NLR family pyrin domain containing 9 | –2.339 | –5.906 |
| GDPD2 | Glycerophosphodiester phosphodiesterase domain containing 2 | –2.28 | –2.539 |
| NEK3 | NIMA related kinase 3 | –2.235 | –3.583 |
| CCDC62 | Coiled-coil domain containing 62 | –2.234 | –3.062 |
| ST8SIA5 | ST8 aspha-N-acetyl-neuraminide alpha-2,8-sialytransferase 5 | –2.232 | –6.472 |
| CHKB | Choline kinase beta | –2.219 | –2.888 |
| IL5RA | Interleukin 5 receptor subunit alpha | –2.218 | –2.015 |
| RPL10A | Ribosomal protein L10A | –2.2 | –1.514 |
| RPS13 | Ribosomal protein s13 | –2.191 | –1.329 |
| PIGH | Phosphatidylinositol glycan anchor biosynthesis, class H | –2.159 | –2.072 |
| LINGO4 | Leucine rich repeat and lg domain containing 4 | –2.159 | –5.562 |
| CCDC152 | Coiled-coil domain containing 152 | –2.131 | –3.069 |
| RSRP1 | Arginine/serine rich protein 1 | –2.13 | –1.756 |
| ALS2CL | Als2 C-terminal like | –2.081 | –1.219 |
| B3GAT2 | Beta-1,3-glucuronyltransferase 2 | –2.074 | –1.656 |
| Aph1c | Aph1 homolog, gamma secretase subunit | –2.059 | –3.541 |
| CYP2C44 | Cytochrome family 2, subfamily c, polypeptide 23 | –2.049 | –5.493 |
| GAL | Galanin | –2.045 | –1.863 |
| SUN3 | Sad1 and UNC84 domain containing 3 | –2.044 | –2.582 |
| FTH1 | Ferritin heavy polypeptide 1 | –2.033 | –1.886 |
| PLSCR4 | Phospholipid scramblase 4 | –2.028 | –3.39 |
| TCF7L2 | Transcription factor 7 like 2 | –2.015 | –2.63 |
| PNP | Purine-nucleoside phosphorylase | –2.008 | –2.019 |
| RHBDL1 | Rhomboid like 1 | –1.997 | –2.056 |
| ARGLU1 | Arginine and glutamate rich 1 | –1.991 | –1.564 |
| GSDMC | Gasdermin C | –1.987 | –5.246 |
| OOEP | Oocyte expressed protein | –1.981 | –5.338 |
| AQP8 | Aquaporin 8 | –1.967 | –9.629 |
Indicates colon cancer biomarker (FDR < 0.05, fold change >|1.5|, IPA.
Figure 6Altered expression of ITGA2, ADAMTS12, and ST8SIA5 transcripts in mouse colon and tumors. (A) Altered levels of select, novel transcripts in FOXO3 KO colon and tumors was confirmed by qPCR (n = 3 per group from 2 independent experiments; *,#P < .05 *relative to WT colon, #relative to WT tumor; analysis of variance and Student-Newman-Keuls post-test). (B) Immunohistostainings of ITGA2, ADAMTS12, ST8SIA5, and Ki67 in colon and tumors of WT and FOXO3 KO mice. In WT and FOXO3 KO colon, gray boxes represent optical zoom images from the crypt base, midsection, and apical region as displayed as panels on the right. Insets in tumor panels (WT and FOXO3 KO) depict higher magnification of tumor and stromal sections (representative staining from colon and tumors of 3 mice from 2 independent immunohistostainings, scale bar 25 μm).
Figure 7Expression of ITGA2, ADAMTS12, and ST8SIA5 transcripts across human IBD and colon cancer. (A) In human IBD tissue samples, levels of ITGA2, ADAMTS12, and ST8SIA5 were significantly altered relative to normal control (GSE4183, n = 23 patient samples, Agilent Microarray Probe ID: ITGA2 (22734_at), ADAMTS12 (226997_at), ST8SIA5 (217514_at); *P < .05 compared with healthy control, Student t test). (B) OncoPrint alteration summary of ADAMTS12, ITGA2, and ST8SIA5 in human colon cancer (TCGA) (cbioportal.org, P < .05). (C and D) mRNA expression of select, novel transcripts ITGA2, ADAMTS12, and ST8SIA5 were significantly altered in human colon cancer compared with matched control tissue (TCGA) (TIMER, ***P < .001 compared with normal tissue) and were associated with poor patient survival as measured by KM estimate (PROGgeneV2) (*P < .05).
Figure 8In human colon cancer HCT116 cells FOXO3 silencing leads to increased ITGA2 levels that promote growth. (A) Human colon cancer HCT116 clones with (dox)-inducible (2 μg/mL, 48 hours) FOXO3-specific shRNA (shFOXO3) were assessed for reduced FOXO3 mRNA levels (n = 4 from 2 independent experiments and clones; *,#P < .05, *compared with control (shCon) (-dox), #compared with shFOXO3 (-dox); analysis of variance and Student-Newman-Keuls post-test). (B) Increased ITGA2 expression in response to FOXO3 knockdown (shFOXO3) in colon cancer cells was determined by qPCR (n = 4 from 2 independent experiments and clones; *,#P < .05, *compared with shCon (-dox), #compared with shFOXO3 (-dox); analysis of variance and Student-Newman-Keuls post-test). (C) Immunoblot of elevated ITGA2 and reduced FOXO3 protein levels in shFOXO3 colon cancer cells (n = 4 from 2 independent experiments and clones; *,#P < .05, *compared with shCon (-dox), #compared with shFOXO3 (-dox); analysis of variance and Student-Newman-Keuls post-test). (D) MTS proliferation assay of colon cancer cells left untreated (control, Con) or transfected with non-specific scramble (NS) siRNA or ITGA2 specific siRNA (n = 8 per condition from 2 independent experiments; *,#P < .05, *compared with Con, #compared with NS Con; analysis of variance and Student-Newman-Keuls post-test).