Following the flank tumor protocol and experimental timeline (Figure 1), the tumor growth and immune response to a targeted tyrosine kinase inhibitor (TKI) therapy and nivolumab combination treatment was studied in two distinct human colorectal cancer (CRC) PDXs. The TKI drugs have been studied in immunodeficient hosts to evaluate tumor growth only29. This model enabled the study of changes in the immune response of the TKI alone, and more importantly, in combination with anti-PD-1. This study was focused on the combination treated cohorts in two distinct experiments that represent a successful evaluation using HIS-BRGS mice and a technically flawed experiment. For PDX CRC307P, the combination treatment slowed the growth of the tumor, as determined by tumor growth volumes over time, tumor weights at harvest, and the SGR (Figure 2A-C). On the other hand, the growth of a PDX developed from a metastatic tumor from the same patient, PDX CRC307M, tested in a different cohort of HIS-BRGS mice, was less affected by the same combination treatment in HIS-BRGS mice (Figure 2D).
Despite the allocation of HIS-BRGS mice into equivalent experimental groups based on overall human (hCD45+) and human T cell (hCD3+) chimerism in the blood prior to tumor implantation (Figure 3A-C), both parameters increased in the peripheral immune system (spleens) and the tumor-infiltrating leukocytes (TIL) in combination treated CRC307P-bearing mice, but not the CRC307M model (Figure 3D-F). Notably, although both HIS-BRGS cohorts had comparable human and T cell chimerism in the blood prior to tumor implantation, the CRC307M cohort had very little lymph node chimerism and failed to develop appreciable levels of T cells in the spleen in the treated cohort (Figure 3D-F). The CRC307M experiment represents a technically flawed example due to overall insufficient splenic T cell chimerism and lymph node chimerism. We suggest the exclusion of HIS mice that have <20% hCD3+ chimerism in the spleens and/or <1.5 x 106 hCD45+ cells in the lymph nodes at the end of the study. For the CRC307M model, the majority of mice were excluded, mostly due to very small lymph nodes, leaving only two mice per cohort.
Further investigation of the human T cells (Figure 4A) revealed more activated (HLA-DR+) T cells in the CR307P tumors, but not the lymph, from combination-treated mice (Figure 4B). In the CRC307M "negative" experiment, there was no significant increase of HLA-DR+ T cells in the TILs of treated HIS-BRGS mice when analyzing all the mice, suggesting this non-optimal HIS-BRGS cohort influenced the data significance due to HIS mice with very few T cells and no activation (Figure 4B). Indeed, the increase in HLA-DR+ T cells in the TILs in the CRC307M model did reach significance when excluding the low T cell chimerism HIS-BRGS mice (Figure 4B). In addition, there were more effector memory CD8+ T cells and fewer TIM-3+ (terminally exhausted) T cells in the combination-treated CRC307P tumors, whereas this difference was not noted in the CRC307M model (Figure 4C,D). In this experiment, no changes in the frequencies of cytotoxic T cells (Granzyme B+ or IFNγ+TNFα+) populations were observed among the combination-treated mice, although higher cytotoxic T cells were observed in untreated (Figure 4E) or treated (data not shown) tumors relative to lymph nodes. Importantly, although there were no differences in frequencies among the T cell populations, higher numbers of cytotoxic T cells were observed in the tumors of treated HIS-CRC307P-BRGS mice, with higher frequencies (Figure 3B) and numbers of human T cells in the tumors. On the other hand, this combination treatment showed no effect on the frequencies of Tregs in either CRC307P lymph organs or tumors, although the CRC307M data showed a trend of reduced Tregs that would need to be validated in another experiment (Figure 4F).
In addition to interrogating the human immune system, immune-related changes on tumor cells were also evaluated using flow cytometry (Figure 5A). By gating on EpCAM+ cells, as described in the protocol, increased expression of both MHC Class I (HLA-ABC) and Class II (HLA-DR) was found on the CRC307P tumor cells excised from combination-treated HIS-BRGS mice (Figure 5B). In the CRC307M model, this same drug treatment induced HLA Class II expression on the tumor cells, although to a lesser degree than in the CRC307P model (Figure 5C). Thus, the combination treatment appears to induce upregulation of MHC class II, independent of T cell infiltration (Figure 5B,C). Notably, MHC expression on the tumor cells was lower than that on human immune cells, a finding consistent with human tumor reports (Figure 5B). Similarly, the combination-treatment resulted in increased PD-L1 expression on the EpCAM+ CRC307P tumor cells (Figure 5D).
Finally, correlations of immune responses with tumor growth were investigated by plotting the SGR of the tumor versus immune parameters. Although the frequency of CD4+ T cells in the tumors of untreated mice showed no correlation with tumor growth, increased CD4+ T cells showed a significant (Figure 6A) correlation with smaller tumor growth, and more specifically the HLA-DR+ activated T cells (Figure 6B), in the combination treated HIS mice.

Figure 1: Schematic illustrating the generation of HIS-BRGS mice and implantation of human CDX or PDX tumors for cancer immunotherapy studies. The timeline is important to ensure the presence of T cells that take months to engraft in the CB-derived HIS mice. Created with Biorender.com. Please click here to view a larger version of this figure.

Figure 2: Analysis of tumor growth. Tumor growth measurements for CRC307P PDX in control (black) or combination treated (red) HIS-BRGS mice quantified by (A) volume over time, (B) weights at the end of the study, or (C) specific growth rate (SGR). (D) SGRs for CRC307M PDX. Data include tumors implanted into both flanks of six vehicle HIS-BRGS mice, seven combination-treated BRGS mice for CRC307P, and only two vehicle and combination-treated mice for CRC307M due to high exclusion rates. Statistical analyses between two independent groups were performed using unpaired, parametric two-group Welch's t-test. **p < 0.01, ****p < 0.0001. Please click here to view a larger version of this figure.

Figure 3: Human immune and T cell chimerism in lymph organs and tumors of CRC PDX tumor-bearing HIS-BRGS mice. (A) Representative flow cytometry analysis of human (hCD45) and mouse (mCD45) hematopoietic and human T (CD3) cells in peripheral blood (PBMC) prior to tumor implantation, and in LNs and tumors (TIL) at the end of the study. (B,C) Equivalent human and T cell chimerism in blood at 14 weeks in HIS-BRGS mice subsequently injected with (B) CRC307P or (C) CRC307M PDXs and untreated or treated with combination immunotherapy (Tx). (D-F) Increased human T cells in lymph organs and tumors (TIL) of combination-treated HIS-BRGS-CRC307P mice but not HIS-BRGS-CRC307M mice. Data show human and T cell frequencies on the Y-axis as a percentage of the parent population in (D) lymph nodes (LNs), (E) spleens (SP), and (F) tumors of individual mice at the end of the study. Statistical analyses between two independent groups were performed using unpaired, parametric two-group Welch's t-test. *p < 0.05, **p < 0.01, ***p < 0.001. Please click here to view a larger version of this figure.

Figure 4: Analysis of T cell activation by flow cytometry in peripheral lymph organs and tumors of CRC PDX bearing HIS-BRGS mice. (A) Representative flow cytometry analysis of T cells in LNs, spleens (SPs), and tumors (TILs) excised from HIS-BRGS mice measuring the following populations: activated T cells (HLA-DR+), naïve T cells (CD45RA+ CCR7+), Tem (CD45RA-, CCR7-), Tregs (CD25+, FoxP3+), and cytotoxic T cells (Granzyme B+, TNFα+, and/or IFNγ+). (B-D) Frequencies of indicated T cell populations in lymph organs or TILs of HIS-CRC307P-BRGS and HIS-CRC307M-BRGS mice: (B) HLA-DR+ activated T-cells, (C) CCR7-CD45RA-Tem CD8+ T cells, (D) inhibitory receptor TIM-3 CD8+ T cells, (E) Granzyme B CD8+ T cells, and (F) CD25+FoxP3+ CD4+ Tregs. For the HIS-CRC307M-BRGS mice in (B), the first data set shows all mice analyzed at harvest, whereas the second data set includes only those with sufficient LN and splenic T cell chimerism. In all graphs, the filled symbols for 307M represent mice with sufficient chimerism, whereas open symbols are excluded HIS mice. Statistical analyses between two independent groups were performed using unpaired, parametric two-group Welch's t-test. *p < 0.05, **p < 0.01, ***p < 0.001. Please click here to view a larger version of this figure.

Figure 5: Flow cytometry analysis of MHC class I, II, and PD-L1 on human tumors in HIS-BRGS mice. (A) Representative flow cytometry dot plots, illustrating gating strategy for measuring expression levels of MHC Class I (HLA-ABC), Class II (HLA-DR), and PD-L1 on epithelial (EpCAM+) human tumors in HIS-BRGS mice. (B-D) Mean fluorescence intensity (MFI) of HLA-ABC (B), HLA-DR (B,C), and PD-L1 (D) on human (hCD45+) and tumor (EpCAM+) cells from dissociated CRC307P (B,D) or CRC307M (C) PDXs excised from HIS-BRGS mice. Statistical analyses between two independent groups were performed using unpaired, parametric two-group Welch's t-test. *p <0.05, **p<0.01, ****p<0.0001. Abbreviation: Tx = Treatment Please click here to view a larger version of this figure.

Figure 6: Analysis of immunocorrelates of response. The SGRs of CRC307P tumor growth in the flanks of untreated (Veh) or combination immunotherapy-treated (Tx) HIS-BRGS mice were plotted versus the frequency of (A) CD4+ or (B) HLA-DR+ T cells as an indicator of immune parameters that correlate with reduced tumor growth (i.e., smaller SGR). A simple linear regression analysis was performed to indicate significant correlations. R2 values indicate degree of correlation. *p < 0.05. Please click here to view a larger version of this figure.
Table 1: Surface staining panel worksheet. Please click here to download this Table.
Table 2: Cell staining panels and flow cytometry gating worksheet. Please click here to download this Table.
Supplemental File 1: Recipes for media and solutions used in the study. Please click here to download this File.
Supplemental File 2: Representative flow cytometry gating for each stain. Please click here to download this File.