M.D Akihito Harusato

Emilie Viennois

Lucie Etienne-Mesmin

Shingo Matsuyama

Hirohito Abo

Satoru Osuka

Nicholas W Lukacs

Yuji Naito

Yoshito Itoh

Jian-Dong Li

Didier Merlin

Andrew T Gewirtz

Ph.D Timothy L Denning
email: tdenning@gsu.edu

Early-life microbiota exposure restricts myeloid-derived suppressor cell-driven colonic tumorigenesis

Keywords: microbiota restricts colonic tumorigenesis Early-life microbiota restricts colonic tumorigenesis microbiota, colitis-associated cancer, chemokines, myeloid-derived suppressor cell

come

Introduction

The gut microbiota consists of trillions of bacteria that co-develop with the host and foster physiological and immunological homeostasis [START_REF] Ley | Worlds within worlds: evolution of the vertebrate gut microbiota[END_REF][START_REF] Chu | Innate immune recognition of the microbiota promotes hostmicrobial symbiosis[END_REF][START_REF] Honda | The microbiota in adaptive immune homeostasis and disease[END_REF]. Microbiota-induced immune responses are particularly important early in life [START_REF] Hasegawa | Transitions in oral and intestinal microflora composition and innate immune receptordependent stimulation during mouse development[END_REF][START_REF] Palmer | Development of the human infant intestinal microbiota[END_REF] when the maternal microbiota can influence prenatal and early post-natal host development [START_REF] Gomez De Aguero | The maternal microbiota drives early postnatal innate immune development[END_REF]. Exposure to microbiota and associated metabolites during this time may have long-lasting effects on immune and non-immune homeostasis of the gut and other organs, and ultimately disease outcomes. Because colitis-associated cancer (CAC) is influenced by the microbiota and inflammatory signaling pathways [START_REF] Gagliani | The fire within: microbes inflame tumors[END_REF], we explored the contribution of early-life microbiota exposure to the development of CAC later in life.

Here, we show that mice raised under germ-free conditions until weaning and then conventionalized with complex gut microbiota into adulthood (exGF mice) were hypersusceptible to CAC when compared to specific pathogen-free (SPF) mice. The increased number and size of colon tumors and intestinal epithelial cell (IEC) proliferation in exGF mice were associated with augmented intratumoral chemokine expression and granulocytic myeloid-derived suppressor cell (G-MDSC) accumulation. Accordingly, CXCR2 neutralization in exGF mice reversed the hyper-susceptibility to CAC. Collectively, these findings highlight a crucial role for early-life microbiota in establishing gut homeostasis that restrains colon cancer in adulthood.

Methods

Mice

C57BL/6 mice were obtained from the Jackson Laboratories and housed in SPF conditions (SPF mice). GF C57BL/6 mice were kept under GF conditions in Park Bioservices isolators.

Conventionalization of GF mice was performed after weaning at the age of 3 weeks [START_REF] Zhan | Gut microbiota protects against gastrointestinal tumorigenesis caused by epithelial injury[END_REF] transferring the GF mice out of GF isolators and co-housing them with SPF mice for 3 weeks (exGF mice). Animal protocols were approved by the IACUC of Georgia State University.

Colitis-associated cancer model

CAC was induced as previously described [START_REF] Viennois | Dietary Emulsifier-Induced Low-Grade Inflammation Promotes Colon Carcinogenesis[END_REF]. Briefly, as shown in Fig. 1A, SPF and exGF mice were intraperitoneally injected with AOM (10 mg/kg body weight) (Sigma) diluted in PBS. Mice were then subjected to two cycles of dextran sodium sulfate (DSS) treatment (MP Biomedicals), in which each cycle consisted of 2.5% DSS for 7 days followed by 14 days with regular water.

After the colitis-associated cancer protocol was complete (day 122 of life; day 80 following initiation of the AOM/DSS regimen), mice were euthanized, and colon length, colon weight, and spleen weights were measured. At day 98 and 122 of life, mice were fasted for 15 hours, and the colonoscopy procedure was performed under general anesthesia with isoflurane (Karl Storz Endoskope, El Segundo, CA).

CXCR2 neutralization in vivo

Mice were injected with anti-CXCR2 (polyclonal goat anti-mouse; a gift from Dr. Lukacs) or goat IgG (isotype control; Bio X Cell) via the intraperitoneal route (500 g twice a week) beginning 4 days after the AOM injection. Antibodies against CXCR2 were generated by immunization of a goat with murine CXCR2 peptide in multiple intradermal sites and purified over a protein A column as previously described [START_REF] Mehrad | CXC chemokine receptor-2 ligands are necessary components of neutrophil-mediated host defense in invasive pulmonary aspergillosis[END_REF]. Isolation of colonic LP cells was performed as previously described [START_REF] Harusato | Phenotypic and functional profiling of mouse intestinal antigen presenting cells[END_REF]. Briefly, colonic tissue was cut into pieces 0.5 cm in length and placed in an orbital shaker for 20 minutes at 37°C in Hank's Balanced Salt Solution (HBSS) with 10% fetal bovine serum (FBS) and 2 mM EDTA in order to remove epithelial cells. This shaking step was repeated and then the remaining tissue was minced and placed in RPMI with 10% FBS, collagenase IV (1 mg/mL; Sigma), and DNase I (40 U/mL; Roche) and shaken for 15 minutes at 37°C. The contents are then filtered through a 100 μm cell strainer directly into a 50-mL conical tube. Each 50-mL conical tube was topped off with HBSS with FBS and centrifuged at 4 °C twice. The cell pellet was resuspended into a cell suspension, which was stained with specific antibodies and analyzed by flow cytometry.

Isolation of colonic lamina propria (LP) cells and flow cytometry

Fluorescence-labeled antibodies were from Thermo Fisher Scientific: CD11b (M1/70), CD11c (N418), CD127 (A7R43), F4/80 (BM8), CD90.2 (53-2.1), Ly6C (HK1.4), Ly6G (1A8), Ki67 (16A8); from Becton Dickinson (BD): MHC-II (M5/114.15.2), CD3 (145-2C11), CD4 (L3T4), CD8 (5H10), CD45 (30F11), E-cadherin (DECMA-1); Biolegend: NK1.1 (PK136), TCR (H57-597). Fc block (2.4G2) was purchased from BD. LIVE/DEAD fixable aqua dead cell stain (Thermo Scientific) was used as a viability dye.

16S rRNA and mRNA sequencing

Ilumina 16S ribosomal (r)RNA gene sequencing was performed on feces isolated from SPF and exGF mice at days 42, 98, and 122, as previously described [START_REF] Chassaing | Dietary emulsifiers impact the mouse gut microbiota promoting colitis and metabolic syndrome[END_REF]. Briefly, quality filtering was performed using the Quantitative Insights Into Microbial Ecology (QIIME; http://qiime.org) software package [START_REF] Caporaso | QIIME allows analysis of high-throughput community sequencing data[END_REF], and sequences were aligned to operational taxonomic units (OTUs) using the UCLUST algorithm [START_REF] Edgar | Search and clustering orders of magnitude faster than BLAST[END_REF] and classified taxonomically using the Greengenes reference database 13_8 [START_REF] Mcdonald | An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea[END_REF]. A single representative sequence for each OTU was aligned, and a phylogenetic tree was built using FastTree [START_REF] Price | FastTree: computing large minimum evolution trees with profiles instead of a distance matrix[END_REF]. The phylogenetic tree was used for computing the unweighted UniFrac distances between samples, rarefaction was performed, and used to compare abundances of OTUs across samples. Rarefaction was performed at 12962 sequences per sample, for which a total of 5650 OTUs were observed (818 +/-203 per sample). Principal coordinates analysis (PCoA) plots were generated through QIIME and used to assess the variation between experimental groups (beta diversity). The PCoA plot representations are 3D and the dot sizes are dependent on the relative localization. Alpha diversity was analyzed using the Shannon index. LEfSe (LDA Effect Size) was used to investigate bacterial members that drive differences between groups, with a significance level of 2 used as a threshold [START_REF] Segata | Metagenomic biomarker discovery and explanation[END_REF]

TUNEL assay

To quantitate the number of apoptotic cells in colonic epithelial cells, mouse proximal colon devoid of tumor were fixed in 10%-buffered formalin for 24 hours at room temperature, embedded in paraffin, sectioned at 5-μm thickness, deparaffinized, and stained for apoptotic nuclei according to the manufacturer's instructions using the In-Situ Cell Death Detection Kit (Sigma Aldrich, St. Louis, MO). TUNEL-positive cells overlapping with DAPI nuclear staining were counted per crypt by using a BZ-X710 microscope (Keyence).

RNA isolation and qPCR

Total RNA was isolated from murine colonic tissues using the Qiagen RNeasy Mini Kit and QIAcube with on-column DNase digestion. The RNA concentration was measured using NanoDrop 8000 (Thermo Scientific). cDNA was generated from 1g of RNA using the Superscript First-Strand Synthesis System for RT-PCR and random hexamer primers (Invitrogen). qPCR was performed with SYBR Green on a StepOnePlus real-time PCR system (Applied Biosystems). Gene expression was normalized to Gapdh and analyzed using 2 1 of [START_REF] Chen | Chitinase 3like-1 expression in colonic epithelial cells as a potentially novel marker for colitisassociated neoplasia[END_REF]. Microarray analysis was performed as previously described [START_REF] Harusato | IL-36gamma signaling controls the induced regulatory T cell-Th9 cell balance via NFkappaB activation and STAT transcription factors[END_REF].

Microarray assay were performed on 20 RNA samples isolated from colon mucosa of 20 patients, including 5 normal controls, 4 quiescent UC, and 11 UC with neoplasia. Inclusion criteria is as follows: a previous clinical diagnosis of UC confirmed by an expert GI pathologist, a disease duration > 7 years, and an extent of disease >20 cm proximal to the anal verge. For microarray analysis, expression and raw expression data were summarized and normalized using the Robust Multi-array Average algorithm from the Bioconductor library for the R statistical programming system. All samples were clustered using Ward's method in Spotfire (TIBCO Software, Palo Alto, CA).

Statistics

Results

Absence of microbiota prior to weaning exacerbates colitis-associated cancer in adult mice

To investigate the contribution of early-life microbiota exposure to CAC, we compared SPF mice to exGF mice (born and raised in GF conditions until weaning at day 21 when they were co-housed with age-and sex-matched SPF mice for an additional 21 days). Both groups of mice were subsequently administered AOM at day 42 of life followed by two 7-day cycles of DSS to induce CAC (Fig. 1A). We observed that exGF mice lost significantly more body weight compared to SPF mice, beginning during the first cycle of DSS (Fig. 1B). ExGF mice also developed significantly larger colonic tumors, as assessed by endoscopy at day 98 and 122 (Fig. 1C) and macroscopic examination (Fig. 1D). Similar differences between exGF and SPF mice were observed for total tumor surface area (Fig. 1E), tumor surface area per individual tumor at day 122 (Fig. 1F), and in the total number of colonic tumors, specifically in the 1-2mm and 2-3mm size range (Fig. 1G). We did not observe significant differences in spleen weights and colon lengths between the groups (Fig. 1H). These data demonstrated that the lack of microbiota before weaning favors CAC development in adulthood.

Enhanced colonic tumor size and epithelial turnover in AOM/DSS-treated exGF mice

Histological examination confirmed the presence of larger adenomas particularly in the middle to distal colonic region and rectum in AOM/DSS-treated exGF mice compared to SPF controls (Supplementary Fig. S1A-B). Given this increased tumor burden observed in the exGF group, we investigated whether colonic epithelial cell death and proliferation were altered in these mice.

We observed a significantly increased number of TUNEL + cells in tissue sections (Supplementary Fig. S1C-D (IECs) (Supplementary Fig. S1E-F) among exGF mice compared to SPF controls. These results suggest that increased IEC turnover in exGF mice may enhance CAC.

ExGF mice exhibit altered microbiota composition and colonic gene expression

We next analyzed the microbiota and host gene expression that may be associated with and/or contribute to the enhanced colonic tumor burden in exGF mice. To analyze microbiota composition, 16S rRNA gene sequencing was performed on feces isolated from SPF and exGF mice at days 42, 98, and 122. Principal coordinate analysis (PCoA) of the unweighted Unifrac distance revealed separate clustering of SPF (blue dots) and exGF samples (red dots) at day 42 preceding initiation of AOM/DSS (Fig. 2A). We observed distinct clustering between SPF and exGF mice (blue and red dots, respectively; Supplementary Fig. S2A) at day 42 (red dots; Supplementary Fig. S2B), with this microbiota clustering becoming more distinct at day 98 and day 122 (blue and orange dots, respectively; Supplementary Fig. S2B). We observed two subgroups for both SPF and exGF animals at day 42 (Fig. 2A), which could be explained by the samples originating from two independent experiments. Shannon index analysis demonstrated that the alpha diversity of exGF mice was not significantly altered compared to SPF mice (Supplementary Fig. S2C). We next analyzed the unweighted Unifrac distance separating these groups in order to quantitate the differences, which confirmed that alterations in the microbiota composition between SPF and exGF mice was present at day 42 and significantly exacerbated at day 98 and 122 (Fig. 2B). Microbiota alteration in exGF mice at day 42 was characterized by increases in members of the families Verrucomicrobiaceae, Rikenellaceae, and Erysipelotrichaceae and reductions in Dehalobacteriaceae, Peptococcaceae, and Ruminococcaceae (Fig. 2C, Supplementary Fig. S2D-E).

on February 21, 2019. © 2019 American Association for Cancer Research.

We next performed RNA-seq analysis of total colon tissue at days 7, 14, 21, and 42 of life, all prior to the initiation of AOM-DSS. Volcano plots indicated that the number of genes significantly upregulated or downregulated in exGF mice by at least 2-fold at day 14 (96 genes up in exGF; 253 genes down in exGF) and day 21 (375 genes up in exGF; 216 genes down in exGF) were higher than at day 7 (96 genes up in exGF; 82 genes down in exGF) and day 42 (97 genes up in exGF; 100 genes down in exGF) (Fig. 2D, Supplementary Fig. S3A). These data are consistent with the expansion of microbiota in SPF mice, but not GF mice, around the time of transitioning away from breast milk to eating solid food (~day 14) and the time of weaning (day 21). At day 42, just prior to the initiation of AOM/DSS, analysis of specific genes that were significantly upregulated in exGF mice revealed higher expression of the monocyte chemokine Ccl8, the antimicrobial peptide Reg3g and several pro-proliferative genes (Mmp3, Wnt2, HOXD3, Kif26b). Alternatively, lower expression of anti-inflammatory genes (Il2, Unc5c) and tumor suppressor genes (Nod2, Lgi3, Klf12) were observed in exGF mice (Fig. 2D).

We further performed heatmap analyses of SPF and exGF colon tissue samples from day 42 and found that 551 genes were differentially expressed (Supplementary Fig. S3B). Pie chart representation of pathways was used to gain insight into specific signaling pathways that may have been altered by early-life microbiota by analyzing genes with a more than 2-fold increase expression in exGF mice compared to SPF mice. This approach identified the 'inflammationmediated by chemokine and cytokine signaling pathway' as one of the most enriched pathways (Supplementary Fig. S3C). Collectively, these data indicated a steady-state shift towards a proproliferative, pro-inflammatory milieu in the colons of exGF mice, which is associated with distinct microbiota composition and enhanced CAC.

Enhanced CXCR2 ligands and granulocytic MDSCs in AOM/DSS-treated exGF mice

In light of enhanced pro-inflammatory gene expression and the alteration in the expression of cytokine/chemokine signaling pathways in exGF mice at day 42, we next explored how specific chemokine signaling may contribute to increased CAC. We focused on the CXCL2/CXCR2 axis because CXCR2-dependent recruitment of myeloid-derived suppressor cells (MDSCs) has been implicated in CAC [START_REF] Katoh | CXCR2-expressing myeloid-derived suppressor cells are essential to promote colitis-associated tumorigenesis[END_REF]. Using qRT-PCR we observed higher CXCL1 and CXCL2 mRNA expression in tumor tissues (T), but not non-tumor tissue (N), obtained from exGF mice compared to SPF mice, and CXCL5 mRNA followed a similar trend (Fig. 3A). To examine whether increased CXCR2 ligand expression was associated with enhanced MDSCs, we next analyzed granulocytic MDSCs (G-MDSCs; CD11b + Ly6G high ) and monocytic MDSCs (M-MDSCs; CD11b + Ly6C + ) [START_REF] Veglia | Myeloid-derived suppressor cells coming of age[END_REF] in colons from AOM/DSS-treated exGF and SPF mice. Following treatment with AOM/DSS, FACS analyses of colonic tissues revealed that exGF mice displayed significantly increased accumulation of G-MDSCs, but not M-MDSCs, compared to SPF mice (Fig. 3B-E).

We next investigated whether CXCR2 ligands were also enhanced during CAC in humans using datasets for colonic tissues obtained from healthy controls (HC) and ulcerative colitis (UC) patients with or without remote colonic neoplasia. In UC patients, the expression of the CXCR2, as well as the CXCR2 ligands CXCL1, CXCL2, CXCL3, CXCL5, and CXCL8, were all elevated in colonic tissue isolated from patients harboring remote colonic neoplasia compared those without remote neoplasia or HC (Supplementary Fig. S4). Taken together, these data suggest that a lack of microbiota in early-life may enhance expression of CXCR2 ligands accumulation of G-MDSCs in colonic tissue.

CXCR2 blockade inhibits colonic tumorigenesis in AOM/DSS-treated exGF mice

We further explored whether in vivo neutralization of CXCR2 could limit the enhanced recruitment of MDSCs, as well as CAC observed in these mice. Indeed, anti-CXCR2 treatment inhibited the development of AOM/DSS-induced macroscopic colonic tumors in exGF mice to a similar extent as that observed in anti-CXCR2-treated SPF mice, when compared to either group treated with control IgG (Fig. 4A). Similar effects of anti-CXCR2 treatment were observed with regards to inhibiting total tumor surface area (Fig. 4B), tumor surface area per individual tumor at day 122 (Fig. 4C), and the total number of colonic tumors (Fig. 4D), whereas no significant differences in spleen weights and colon lengths between the two groups were seen (Fig. 4E). We confirmed that injection of anti-CXCR2 into AOM/DSS-treated SPF and exGF mice prevented G-MDSC accumulation in colonic tissue (Fig. 4F).

Discussion

A previous report demonstrated that AOM/DSS-treated GF mice develop more CAC compared to SPF mice and linked this phenotype to impaired inflammatory responses and delayed repair processes [START_REF] Zhan | Gut microbiota protects against gastrointestinal tumorigenesis caused by epithelial injury[END_REF]. Our complementary findings show that exGF mice had an altered microbiota, enhanced chemokine expression and G-MDSCs in the colon, and were hyper-susceptible to CAC. These data indicate that early-life exposure to microbiota induced enduring microbial and host changes that predisposed mice to CAC.

The microbiota differences in exGF mice were characterized by increases in Verrucomicrobiaceae, Rikenellaceae, and Erysipelotrichaceae. These results are consistent with reports that Akkermansia, Bacteroidales, and Erysipelotrichaceae correlate with increased CAC [START_REF] Baxter | Structure of the gut microbiome following colonization with human feces determines colonic tumor burden[END_REF][START_REF] Chen | NLRP12 attenuates colon inflammation by maintaining colonic microbial diversity and promoting protective commensal bacterial growth[END_REF]. Both Bacteroides and Akkermansia are mucin degraders, and an overabundance of these bacteria may impair mucosal barrier defenses leading to increased inflammation and tumorigenesis in exGF mice. Alternatively, Dehalobacteriaceae, Peptococcaceae, and Ruminococcaceae, which have been linked to protection from CAC [START_REF] Richard | Mucosaassociated microbiota dysbiosis in colitis associated cancer[END_REF][START_REF] Fernandez | A Galacto-Oligosaccharides Preparation Derived From Lactulose Protects Against Colorectal Cancer Development in an Animal Model[END_REF][START_REF] Alteber | The anti-inflammatory IFITM genes ameliorate colitis and partially protect from tumorigenesis by changing immunity and microbiota[END_REF], were significantly reduced in exGF mice. These findings provide insight into how the altered microbiota in exGF mice may lead to enhanced CAC.

Collectively, our results suggest that CAC susceptibility in adulthood may be imprinted in the prenatal and/or early postnatal period by maternal-derived microbiota [START_REF] Macpherson | How nutrition and the maternal microbiota shape the neonatal immune system[END_REF]. Additional insights into how early-life exposure to microbiota regulates intestinal homeostasis, inflammation, and tumorigenesis may lead to opportunities for targeted interventions [START_REF] Torow | The Neonatal Window of Opportunity: Setting the Stage for Life-Long Host-Microbial Interaction and Immune Homeostasis[END_REF][START_REF] Tamburini | 8 males and 7 females) injected with control antibody (IgG) or a CXCR2 neutralizing antibody (500 g twice a week; SPF/cont IgG, n=8; SPF/aCXCR2, n=8; exGF/cont IgG[END_REF].

Disclosure of Potential Conflicts of Interest

No potential conflicts of interests were disclosed by any of the authors. values corresponding to increased expression in exGF mice compared to SPF mice and negative values corresponding to decreased expression in exGF mice compared to SPF mice. Significance between the two groups is indicated by -log10 P-value on the y-axis. Red dots correspond to genes with a P-value <0.05 between exGF and SPF mice. Blue dots correspond to genes with at least a 2-fold decreased or increased expression in exGF and SPF mice. Green dots correspond to genes with at least a 2-fold decreased or increased abundance in exGF and SPF mice and a P-

  ), and increased Ki67 staining of isolated intestinal epithelial cells on February 21, 2019. © 2019 American Association for Cancer Research.

Figure 1 .

 1 Figure 1. Absence of microbiota prior to weaning exacerbates CAC in adult mice. A,

Figure 2 .

 2 Figure 2. ExGF mice exhibit altered microbiota composition and intestinal gene expression.

  value <0.05. One-way ANOVA and Bonferroni post-hoc test was used to determine significance. Data are from (C,D) one or (A,B) two independent experiments.

Figure 3 .

 3 Figure 3. CXCR2 ligands and G-MDSCs are enhanced in colons of AOM/DSS-treated

Figure 4 .

 4 Figure 4. Blockade of CXCR2 ameliorates carcinogenesis in exGF mice. A, Gross view of

  days) by on February 21, 2019. © 2019 American Association for Cancer Research. cancerimmunolres.aacrjournals.org Downloaded from Author manuscripts have been peer reviewed and accepted for publication but have not yet been edited. Author Manuscript Published OnlineFirst on February 19, 2019; DOI: 10.1158/2326-6066.CIR-18-0444

  Author manuscripts have been peer reviewed and accepted for publication but have not yet been edited. Author Manuscript Published OnlineFirst on February 19, 2019; DOI: 10.1158/2326-6066.CIR-18-0444

									. The
	mRNA-seq experiments were performed by Novogene (Beijing, China). Volcano plots were
	generated	using	R,	and	heatmaps	were	generated	using	Morpheus

(https://software.broadinstitute.org/morpheus/). Pie chart representation of pathways was performed using Panther Classification System (http://pantherdb.org) to gain insight into specific signaling pathways. The accession numbers for the unprocessed transcriptomic sequencing data and the unprocessed 16S sequencing data reported in this paper are in the European Nucleotide Archive: PRJEB30607 and PRJEB30606, respectively.

Histology

Mouse distal colons with tumors were fixed in 10%-buffered formalin for 24 hours at room temperature and subsequently embedded in paraffin. Sectioning, hematoxylin/eosin staining, slide scanning, and assessment of histology were subsequently performed at the Institute for Biomedical Sciences histology core at Georgia State University by using a BX61 microscope (OLYMPUS). on February 21, 2019. © 2019 American Association for Cancer Research. cancerimmunolres.aacrjournals.org Downloaded from

  Author manuscripts have been peer reviewed and accepted for publication but have not yet been edited. Author Manuscript Published OnlineFirst on February 19, 2019; DOI: 10.1158/2326-6066.CIR-18-0444 Microarray datasets (GSE37283; (17)) were downloaded from the Gene Expression Omnibus database (ttps://www.ncbi.nlm.nih.gov/geo/). Detailed patient demographic information is available in Table

	Human microarray analysis
	-ΔΔCT
	method with duplicates.
	Primers used were:
	cxcl1 (F: GCTTGAAGGTGTTGCCCTCAG; R: AAGCCTCGCGACCATTCTTG)
	cxcl2 (F: CGCTGTCAATGCCTGAAG; R: GGCGTCACACTCAAGCTCT)
	cxcl5 (F: GGTCCACAGTGCCCTACG; R: GCGAGTGCATTCCGCTTA)
	gapdh (F: TGGCAAAGTGGAGATTGTTGCC; R: AAGATGGTGATGGGCTTCCCG)

on February 21, 2019. © 2019 American Association for Cancer Research. cancerimmunolres.aacrjournals.org Downloaded from

  Statistical analyses were performed with GraphPad Prism software, version 6.0b (Graphpad Software). Student's t test or one-way ANOVA and Tukey's Multiple Comparison Test was used to determine significance. *P < 0.05, **P < 0.01, ***P < 0.001; ns, not significant. Author manuscripts have been peer reviewed and accepted for publication but have not yet been edited. Author Manuscript Published OnlineFirst on February 19, 2019; DOI: 10.1158/2326-6066.CIR-18-0444

on February 21, 2019. © 2019 American Association for Cancer Research. cancerimmunolres.aacrjournals.org Downloaded from

  on February 21, 2019. © 2019 American Association for Cancer Research. cancerimmunolres.aacrjournals.org Downloaded from Author manuscripts have been peer reviewed and accepted for publication but have not yet been edited. Author Manuscript Published OnlineFirst on February 19, 2019; DOI: 10.1158/2326-6066.CIR-18-0444

  on February 21, 2019. © 2019 American Association for Cancer Research. cancerimmunolres.aacrjournals.org Downloaded from

  on February 21, 2019. © 2019 American Association for Cancer Research. cancerimmunolres.aacrjournals.org Downloaded from Author manuscripts have been peer reviewed and accepted for publication but have not yet been edited. Author Manuscript Published OnlineFirst on February 19, 2019; DOI: 10.1158/2326-6066.CIR-18-0444

on February 21, 2019. © 2019 American Association for Cancer Research. cancerimmunolres.aacrjournals.org Downloaded from Author manuscripts have been peer reviewed and accepted for publication but have not yet been edited. Author Manuscript Published OnlineFirst on February 19, 2019; DOI: 10.1158/2326-6066.CIR-18-0444

Acknowledgements

This work was supported by: a National Institutes of Health grant DK097256 (to T.L.D.); a Crohn's and Colitis Foundation Research Fellowship Award (to A.H.) and a Career Development Award (to E.V.); Japanese Society for the Promotion of Science KAKENHI Grants 17H07013 and 18K15128 (to A.H.); and a Senior Research Career Scientist Award BX00447601 from the Veterans Health Administration (to D.M.). on February 21, 2019. © 2019 American Association for Cancer Research. cancerimmunolres.aacrjournals.org Downloaded from Author manuscripts have been peer reviewed and accepted for publication but have not yet been edited. Author Manuscript Published OnlineFirst on February 19, 2019; DOI: 10.1158/2326-6066.CIR-18-0444