- 1Department of Cell Biology and Physiology, The University of North Carolina at Chapel Hill School of Medicine, Chapel Hill, NC, United States
- 2Lineberger Comprehensive Cancer Center, The University of North Carolina at Chapel Hill School of Medicine, Chapel Hill, NC, United States
- 3Department of Microbiology and Immunology, The University of North Carolina at Chapel Hill School of Medicine, Chapel Hill, NC, United States
- 4Department of Microbiology and Immunology, The University of Iowa, Iowa City, IA, United States
- 5Holden Comprehensive Cancer Center, The University of Iowa, Iowa City, IA, United States
- 6Department of Genetics, The University of North Carolina at Chapel Hill School of Medicine, Chapel Hill, NC, United States
B cells can act as potent suppressors of anti-tumor T cell immunity, presenting a mechanism of resistance to immunotherapy. In pancreatic ductal adenocarcinoma, B cells can display a T cell-suppressive or regulatory phenotype centered on the expression of the cytokine Interleukin 35 (IL-35). While B cell-mediated immunosuppression presents a barrier to anti-tumorigenic T cell function, it is not clear how regulatory B cell function could be targeted, and the signals that promote this suppressive phenotype in B cells are not well understood. Here we use a novel IL-35 reporter model to understand which signaling pathways are important for immunosuppressive properties in B cells. In vitro analysis of IL-35 reporter B cells revealed a synergy between the BCR and TLR4 signaling pathways is sufficient to induce IL-35 expression. However, in vivo, B cell receptor activation, as opposed to MyD88 signaling in B cells, is central to B cell-mediated suppression and promotion of pancreatic cancer growth. Further analysis identified protein kinase D2 (PKD2) as being a key downstream regulator of IL-35 expression in B cells. Regulatory B cells with an inactivating mutation in PKD2 failed to produce IL-35 or fully suppress effector T cell function in vitro. Furthermore, inhibition of PKD in B cells decreased tumor growth and promoted effector T cell function upon adoptive transfer into B cell-deficient mice. Collectively, these data provide insight into how regulatory B cell function is promoted in pancreatic cancer and identify potential therapeutic targets to restrain this function.
Introduction
Pancreatic ductal adenocarcinoma (PDAC) is an aggressive malignancy with few viable treatment options for late-stage tumors resulting in a 10% 5-year survival rate (1). Pancreatic tumorigenesis is fueled by a robust infiltration of immunosuppressive cells such as regulatory T cells, myeloid-derived suppressor cells, and M2 macrophages (2–4). We previously identified B cells as key suppressors of T cell infiltration and effector function within pancreatic tumors, predominantly through production of the cytokine IL-35 (5, 6). However, the key regulators of this suppressive activity by B cells in the tumor microenvironment are still unclear.
Immunosuppressive or ‘regulatory’ B cells were originally identified in many disease contexts outside of cancer. The contribution of these B cell populations on disease regulation has been primarily linked to expression of anti-inflammatory IL-10 by to suppress host T cell responses (7–10). We and others have shown that B cells can also suppress T-cell activity in autoimmunity, infection, and cancer via the cytokine IL-35 (5, 6, 11–13). PDAC progression specifically, is dependent on IL-35 and not IL-10 production by B cells (5, 13). IL-35, a heterodimeric cytokine from the IL-12 family of cytokines, consists of p35 (Il12a) and EBi3 (Ebi3) subunits that potently suppress T cell immunity (14). Expression of the IL-35 subunits and IL-10 by B cells has been linked to stimulation of the B cell receptor (BCR) (8, 15), Toll-like receptors (TLRs) (10, 12, 16), and/or CD40 (12, 17). However, it is difficult to dissect precise expression pattern of the IL-35 subunits in the same B cell, due to the cytokine’s heterodimeric nature. To enable studies on regulation of IL-35 expression in B cells, we generated an Il12a-GFP reporter mouse (18) to accompany a previously generated Ebi3-TdTomato reporter mouse (19).
The exact antigen receptor signals acting on B cells to elicit immunosuppression have not been identified, but plausible BCR and/or TLR specific antigens are present within the pancreatic tumor microenvironment (20–24). These signals come in the forms of tumor-associated antigens, damage-associated molecular patterns, and a tumor-associated pancreatic microbiota (25–27). The roles of the BCR and TLRs on regulatory B cell activity in PDAC and ultimately promotion of pancreatic tumor growth have not been thoroughly examined. This study addresses the respective roles of BCR and TLR4 signaling by B cells in pancreatic cancer and how they regulate expression of IL-35.
Materials and Methods
Mouse Models
All mouse protocols were reviewed and approved by the Institutional Animal Care and Use Committee of the University of North Carolina at Chapel Hill (Chapel Hill, NC, USA). Animals were maintained in a specific pathogen-free facility. Six- to 8-week-old wild-type (WT) C57Bl/6J mice were purchased from The Charles River Laboratories (strain #027). 6- to 8-week-old MD4 (strain #002595), µMT (strain #002288), Prkd2S707A/S711A (strain #017285), CD19Cre/+ (strain #006785), and Myd88Fl/Fl (strain #008888) mice on C57Bl/6J background were purchased from Jackson Laboratories. Ebi3Tom/Tom mice were obtained from D. A. Vignali (University of Pittsburgh, Pittsburgh, PA) (19). Il12aGFP/GFP mice were generated at the University of North Carolina at Chapel Hill (Chapel Hill, NC, USA) (18). Ebi3Tom/Tom; Il12aGFP/GFP mice were generated by breeding Ebi3Tom/Tom mice with Il12aGFP/GFP mice until progeny were homozygous for both reporter alleles. CD19Cre/+; Myd88Fl/Fl mice were generated by breeding CD19Cre/+ mice with Myd88Fl/Fl mice until progeny were homozygous for the floxed Myd88 allele. CD19+/+; Myd88Fl/Fl littermates lacking Cre expression were used as controls.
Cell Lines
The murine PDAC cell line KPC4662 (p48Cre/+; LSL-KrasG12D/+; LSL-Tp53R173/+)was derived from a primary pancreatic tumor of C57Bl/6J KPC mice by Dr. Robert Vonderheide’s laboratory (University of Pennsylvania Perelman School of Medicine, Philadelphia, PA.) (28). Cells were maintained at 37°C and 5% CO2 in DMEM (#11995-065, Gibco) containing 10% Fetal calf serum (Corning) and 1x penicillin– streptomycin (#15140-122, Gibco). Cells were confirmed to be Mycoplasma and endotoxin free. KPC4662 Cells that were injected orthotopically were at <16 passages from original derivation.
Orthotopic Tumor Modelling
For intrapancreatic injection of KPC4662 cancer cells, mice were anesthetized using a ketamine (100 mg/kg; Med-Vet International) and xylazine (10 mg/kg; Med-Vet International) cocktail. The depth of anesthesia was confirmed by verifying an absence of response to toe pinch. An incision in the left flank was made, and 1 x 105 KPC4662 cells in ice-cold PBS mixed at 1:1 dilution with Matrigel (#354234, Corning) in a volume of 50 µL were injected using a 28-gauge needle into the tail of the pancreas. The wound was closed in two layers, with running 5-0 Vicryl RAPIDE sutures (Ethicon) for the body wall, and 5-0 PROLENE sutures (Ethicon) for the skin. All mice were given buprenorphine (0.1 mg/kg; Med-Vet International) subcutaneously after orthotopic surgery for pain relief. Tumors were grown in mice for 21 days unless otherwise noted. Tumor volume was measured using electronic calipers after animal sacrifice and calculated using the formula: Volume = Length * (Width2)/0.52.
Depletion of CD4+ T Cells and CD8+ T Cells In Vivo
For CD4+ and CD8+ T-cell depletion studies, 200 µg of αCD4 (#BP0003-1, clone GK1.4, Bio X Cell) and 200 µg of αCD8 (#BP0004-1, clone 53-6.7, Bio X Cell) or an IgG isotype control rat IgG2b (#BE0086, clone MPC-11, Bio X Cell) and rat IgG2a (#BE0089, clone 2A3, Bio X Cell), respectively, were administered intraperitoneally daily starting 3 days prior to tumor cell injection and twice a week after tumor cell injection. Mice were sacrificed 21 days after tumor implantation. Depletion of T cells was confirmed by flow cytometry at time of animal sacrifice.
Ultrasound Imaging of Orthotopic Tumors
Monitoring of orthotopic pancreatic tumors was performed weekly on the Vevo 2100 Imaging System (VisualSonics Inc.) using the MS550D transducer. Mice were anaesthetized using isoflurane (2%) throughout the procedure. Hair was removed from left flank of each mouse by electric razor and hair removal cream (Nair, Church & Dwight). Images were taken at 11 mm image depth and measurements were calculated using the Vevo LAB software (VisualSonics Inc.).
Lymphocyte Isolation
Single-cell suspensions were prepared from tumors, spleens, and draining lymph nodes isolated from orthotopic models. Spleens and draining lymph nodes were mechanically disrupted using the plunger end of a 5 mL syringe and resuspended in 2% FCS/PBS. Pelleted samples were then depleted of red blood cells using RBC lysis buffer (eBioscience; 00-4333-57) and resuspended in 2% FCS/PBS. For isolation of tumor-infiltrating lymphocytes, tumor tissue was thoroughly minced using sterile razor blades and digested with collagenase IV (1.25 mg/mL; #LS004188, Worthington), 0.1% soybean trypsin inhibitor (#T9128, Sigma), hyaluronidase (1 mg/mL; #LS002592, Worthington), and DNase I (100 mg/mL; #LS002007, Worthington) in complete DMEM for 30 minutes at 37°C. Cell suspensions were passed through a 70µm cell strainer (#352350, Falcon) and resuspended in RPMI 1640 media (Gibco). Leukocytes were isolated from processed tumor tissues by Ficoll-Paque PLUS (#17-1440-03, GE Healthcare) density gradient centrifugation and resuspended in 2% FCS/PBS. For FACS sorted cell populations, cells were stained with fluorophore-labeled antibodies for 30 minutes on ice in FACS buffer (2% FCS/PBS). After staining, cells were washed twice with FACS buffer and resuspended in cell sorting buffer (1% FCS/PBS). Lymphocyte populations were sorted using a BD FACS ARIA III sorter to isolate CD19+CD21HiCD5+CD1dHi Breg cells, CD19+CD21LoCD5-CD1dLo Bcon cells, CD4+ T cells, and CD8+ T cells. Cells were collected in complete RPMI media containing 10% FCS with 1x penicillin–streptomycin (#15140-122, Gibco) antibiotics and 1x 2-β-Mercaptoethanol. For assays using total B cells, B cells were isolated from tissues using the EasySep Mouse B cell Isolation Kit (#19854, StemCell Technologies) according to manufacturer’s instructions. For adoptive transfer experiments, 6- to 8-week-old µMT mice were intravenously reconstituted with 1 x 107 purified DMSO or CRT0066101 (PKDi) treated WT splenic B cells 16 hours prior to orthotopic injection. B cells were purified by negative magnetic selection (#19854, StemCell Technologies) and purity was >95% by flow cytometry. Both male and female mice were used in all studies.
B Cell and T Cell In Vitro Cell Culture
Unless otherwise noted, B cells were activated with αCD40 (1 µg/mL, clone HM40-3, BioLegend), αIgM, µ chain (10µg/mL, #115-006-075, Jackson Immuno Research) and/or LPS (2 µg/mL, #L4391-1MG, Sigma) in complete RPMI for 48 hours as previously described (5, 6). For CD8+ and CD4+ T-cell culture, 1 x 105 T cells were cultured in a 96-well plate with plate bound αCD3 (1 µg/mL, #BE0001-1, Clone 145-2C11, Bio X Cell) and soluble αCD28 (2 µg/mL, #BE0015-1, Clone 37.51, Bio X Cell) in complete RPMI for 48 hours. For detection of cytokines in B cells, 1 x 105 cells were cultured using conditions indicated in figure legends for 48 hours in complete RPMI medium followed by the addition of PMA (50 ng/mL; #P8139, Sigma), Ionomycin (200 ng/mL; #9995S, Cell Signaling Technologies), and Brefeldin A (1x, #420601, BioLegend) for the final 5 hours at 37°C. For detection of cytokines in T cells, 1x 105 cells were cultured with αCD3 and αCD28 for 48 hours with 1x Brefeldin A being added for the final 5 hours of culture. For transwell co-culture assays, a total of 1 x 105 CD19+CD21HiCD5+CD1dHi B cells and 1 x 105 CD8+ T cells or 1 x 105 CD4+ T cells were co-cultured in 96-well Transwell plates (#3381; Corning) with B cells occupying the top chamber and CD4+ or CD8+ T cells in the bottom chamber for 48 hours. For experiments utilizing kinase inhibitors, splenic B cells were cultured for 45 minutes with the indicated inhibitor prior to 48-hour stimulation in vitro. The following inhibitors were used: Selleckchem - ERKi (SCH772984, #S7101), PI3Ki (Idelalisib, #S2226), BTKi (Ibrutinib, #S2680); Cell Signaling Technologies – NFATi (FK-506, #9974S); Sigma – LYNi (Dasatinib, CDS023389-25MG); NF-kBi (Compound A) was a gift from Dr. Albert Baldwin (29). Inhibitor concentration was determined by the lowest amount of inhibitor required to decrease target activation of at least 75% from DMSO treated condition. Inhibition was quantified using ImageStudio software (LI-COR).
Immunophenotyping by Flow Cytometry
Isolated immune cells were washed and blocked with αCD16/CD32 (100ug/100,000 cells, #101320, BioLegend) for 5 minutes on ice and then stained with labeled antibodies against surface markers or cytokines on ice for 30 minutes in FACS buffer (2% FCS/PBS). Intracellular staining for Foxp3 was performed using a Foxp3/Transcription Factor Staining Kit (#00-5523, eBioscience) following surface marker staining and was used according to manufacturer’s instructions. The following monoclonal antibodies directed against mouse antigens were used for flow cytometry: BioLegend - CD45 (clone 30-F11), CD3 (17A2), CD4 (GK1.5), Foxp3 (FJK-16S), CD8 (53-6.7), NK1.1 (PK136), CD11c (N418), CD19 (6D5), CD5 (53-7.3), CD21/35 (7E9), IgD (11–26), IgM (RMM-1), CD23 (B3B4), CD138 (281-2); BD Biosciences - CD1d (clone 1B1), CD93 (AA4.1), CD11b (M1/70). Viability of cells was determined by staining with either Zombie NIR (#423105, BioLegend) or Zombie Red (#423113, BioLegend) viability dye. All samples were acquired on LSRII Fortessa (BD Bioscience) at UNC Flow Cytometry Core Facility and analyzed by FlowJo version 10.2 (TreeStar, Inc.).
Detection of Intracellular Cytokines by Flow Cytometry
Cells were washed and blocked with αCD16/CD32 (100ug/100,000 cells, #101320, BioLegend) for 5 minutes on ice. Cells were then washed and stained with labeled antibodies against surface markers on ice for 30 minutes in FACS buffer (2% FCS/PBS). After surface staining, cells were fixed and permeabilized using the Cytofix/Cytoperm Kit (#554714, BD Biosciences) for 20 minutes at 4°C. Intracellular staining was performed using fluorophore-conjugated cytokine antibodies for 1 hour at 4°C. Cells were then washed twice and resuspended in cell staining buffer (#420201, BioLegend). The following antibodies directed against mouse antigens were used for flow cytometry: BioLegend - IL-10 (clone JES5-16E3), IFNy (clone XMG1.2), TNFα (clone MPX-XT22); eBioscience - p35 (clone 4D10p35); R&D Systems - EBi3 (clone 355022). All samples were acquired on LSRII Fortessa (BD Bioscience) at the UNC Flow Cytometry Core Facility and analyzed by FlowJo version 10.2 (TreeStar, Inc.).
Western Immunoblotting
Cells were lysed using RIPA Buffer (#89901, Thermo Fisher) containing HALT phosphatase/protease inhibitors (#78440, Thermo Fisher) on ice and then centrifuged for 10 minutes at 13000rpm for 10 minutes to extract supernatant. Lysates were quantified using Pierce BCA Protein Assay (#23225, Thermo Fisher) according to manufacturer’s instructions. 40µg of total protein was boiled for 5 minutes at 95°C, mixed with 6x laemmli loading buffer, and then loaded and into polyacrylamide gels. Protein was transferred onto LF-PVDF membranes (#162-0264, Bio-Rad) followed by overnight primary antibody incubation. Primary antibody was probed with fluorescent secondary antibodies and detected using the LI-COR Odyssey imager (LI-COR). Primary and secondary antibodies were diluted in 5% BSA/TBS for blocking. The following primary antibodies were used for immunoblotting: Cell Signaling Technologies – pErk1/2 Thr202/Tyr204 (D13.14.4E), Erk1/2 (L34F12), Nfat1 (D43B1), pP65 Ser536 (93H1), P65 (L846), pAkt Ser473 (D9E), Akt (40D4), pSyk Tyr525/526 (C87C1), pSrc Tyr416 (D49G4), Src (L4A1), pBtk Tyr223 (D9T6H), Btk (D6T2C), pPKD/PKCμ Ser744/748 (#2054S), β-Actin (13E5). The following secondary antibodies were used for primary antibody detection: IRDye 800CW Goat anti-Rabbit IgG (#925-32211, LI-COR), IRDye 680RD Goat anti-Mouse IgG (#925-68070, LI-COR).
Quantitative PCR Analysis of Gene Expression
RNA was prepared from splenic, tumoral, and draining lymph node B cells of CD19Cre/+; Myd88Fl/Fl and CD19+/+; Myd88Fl/Fl tumor-bearing mice using the RNeasy Micro Kit (#74004, Qiagen). cDNA was generated from 1µg of total RNA using Maxima First-Strand cDNA synthesis RT Kit (#K1672, Thermo Fisher). qPCR analysis was performed using Perfecta SYBR Green FastMix (95072-250, Quantabio) on the QuantStudio 6 instrument (Applied Biosystems). Results were normalized to the expression of Tbp as an internal control, and each sample was run in triplicate. Gene expression was determined by the ΔΔCt method (2–ΔΔCt). The following primers were used for PCR amplification: Tbp FWD 5’-AGAACAATCCAGACTAGCAGCA–3’, Tbp REV 5’-GGGAACTTCACATCACAGCTC-3’, Il10 FWD 5’-GCTCTTACTGATGGCATGAG-3’, Il10 REV 5’-CGCAGCTCTAGGAGCATGTG-3’, Ebi3 FWD 5’- CTTACAGGCTCGGTGTGGC-3’, Ebi3 REV 5’-GTGACATTTAGCATGTAGGGCA-3’, Il12a FWD 5’-CATCGATGAGCTGATGCAGT-3’ Il12a REV 5’-CAGATAGCCCATCACCCTGT-3’, Prkd1 FWD 5’-GGGGGCATCTCGTTCCATC-3’, Prkd1 REV5’- GTGCCGAAAAAGCAGGATCTT, Prkd2 FWD 5’- GGGGTCTCCTTCCATATCCAG-3’, Prkd2 REV 5’-ACGATAGAACAGGCTAGTTGC-3’, Prkd3 FWD 5’- GTCTGTCAAATGTATCTCTGCCA-3’, Prkd3 REV 5’- GGTGAGTATGTGACTCTTCACTG-3’, Il10 FWD 5’-CCCATTCCTCGTCACGATCTC-3’, Il10 REV 5’-TCAGACTGGTTTGGGATAGGTTT-3’, Tgfb1 FWD 5’-CTCCCGTGGCTTCTAGTGC-3’, Tgfb1 REV 5’-GCCTTAGTTTGGACAGGATCTG-3’, Pdcd1 FWD 5’-ACCCTGGTCATTCACTTGGG-3’, Pdcd1 REV 5’-CATTTGCTCCCTCTGACACTG-3’, Pdcd1l1 FWD 5’-GCTCCAAAGGACTTGTACGTG-3’, Pdcd1l1 REV 5’-TGATCTGAAGGGCAGCATTTC-3’, Il12b FWD 5’-TGGTTTGCCATCGTTTTGCTG-3’, Il12b REV 5’-ACAGGTGAGGTTCACTGTTTCT-3’, Il27a FWD 5’-CTGTTGCTGCTACCCTTGCTT-3’, Il27a REV 5’-CACTCCTGGCAATCGAGATTC-3’.
Statistical Analysis
Statistical analysis was performed using GraphPad Prism software and the statistical test used is indicated in the figure legends. P value of < 0.05 was considered statistically significant. The following denotations were used to further clarify statistical significance: ns=p>0.05, *=p<0.05, **=P<0.01, ***=p<0.005, ****=p<0.001. Unpaired student’s t-test was used to compare experiments with exactly two groups. One-way ANOVA was used to compare experiments with more than two groups. Two-way ANOVA was used to compare experiments with more than two groups tested under more than one condition. Sizes of animal or treatment groups are indicated in figure legends.
Results
BCR and TLR4 or CD40 Activation Is Sufficient to Drive Breg IL-35 Expression Ex Vivo
To understand how the expression of IL-35 subunits is regulated at single cell resolution, we generated a bona fide IL-35 reporter mouse by crossing our Il12a-GFP reporter mouse (18) to a previously generated and validated Ebi3-TdTomato reporter mouse (19). We tested the responsiveness of primary IL-35-reporter B cells to signals previously linked to either IL-10 and/or IL-35 expression in vitro (12, 17, 18). We previously identified the CD19+CD1dHiCD21HiCD5+ subset of B cells (regulatory B cells, Breg) as being the dominant B cell source of IL-35 expression in PDAC (5, 6). We utilized primary splenic CD19+CD1dHiCD21HiCD5+ B cells from IL-35 reporter mice for ex vivo stimulation assays featuring combinations of BCR, TLR4, and CD40 stimulation (Figure 1A). Using this model, we determined that BCR stimulation (αIgM) alone is not sufficient to induce IL-35 expression as the reporter signal intensity was not significantly increased after stimulation, due to fluctuation in the frequencies of individual TdTomato+ and GFP+ populations (Figure 1B, C). However, LPS stimulation or CD40 stimulation alone in the absence of a BCR signal was sufficient to expand a population of GFP+TdTomato+ Breg cells, with most significant changes in GFP MFI signal (Figures 1B–E). Combining LPS and CD40 stimulation led to the most significant increase in GFP+TdTomato+ cell frequency, number, and reporter signal intensity (Figures 1B–E). The use of αIgM instead of LPS as the main antigenic signal accompanying CD40 stimulation did generate a significant GFP+TdTomato+ cell population, but not to the same magnitude of LPS + αCD40 stimulation. Changes in the frequency of GFP+TdTomato+ reporter cells in the presence of LPS and CD40 activation were primarily driven by increase in TdTomato signal. We also examined the effect of combining BCR and TLR stimulation in the presence or absence of CD40 co-stimulation. Interestingly, combined BCR and TLR stimulation resulted in similar reporter signal intensity to LPS stimulation, but the combination resulted in significantly more double reporter-positive cells (p=0.02) and higher frequency (p=0.05) (Figures 1B–E). This suggests that BCR stimulation helps further drive the expansion of Bregs but requires TLR co-stimulation to do so.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/745873/fimmu-12-745873-HTML/image_m/fimmu-12-745873-g001.jpg)
Figure 1 Combined BCR/TLR signaling induces IL-35 expression in B cells. (A) Representative flow cytometry plots of splenic Il12aGFP/GFP; Ebi3Tom/Tom CD19+CD21HiCD1dHiCD5+ B cells after stimulation. (B) Frequency of GFP+TdTomato+ B cells from (A). N=3. (C) Number of GFP+TdTomato+ B cells from (A). N=3. (D) TdTomato mean fluorescent intensity (MFI) in B cells from (A). N=3. (E) GFP mean fluorescent intensity (MFI) in B cells from (A). N=3. (F) Representative flow cytometry plots of WT splenic CD19+CD21HiCD1dHiCD5+ (top) and CD19+CD21LoCD1dLoCD5- (bottom) B cells intracellularly stained for p35 and EBi3 after stimulation. (G) Frequency of p35+EBi3+ B cells from (F). N=3. (H) Frequency (top) and number (bottom) of GFP+TdTomato+ B cells from Il12aGFP/GFP; Ebi3Tom/Tom B cells derived from tumor-naïve and KPC4662 tumor-bearing mice after ex vivo stimulation. N=3. (I) TdTomato mean fluorescent intensity (MFI) in B cells from (H) (top) and GFP MFI in B cells from (H) (bottom) after indicated ex vivo stimulation. N=3. Error bars indicate SEM. P values for 1B-E, G were calculated using one-way ANOVA. P values for 1H-I were calculated using two-way ANOVA. NS, non-significant, *p<0.05, **p<0.005, ***p<0.001, ****p<0.0001. Experiments were performed using 7–8-week-old mice of indicated genotypes. Experiments were repeated 3 times.
Intracellular cytokine staining mirrored reporter findings to show that BCR stimulation alone was not sufficient to induce IL-35 (p35+EBi3+) expression in the regulatory B cell subset (Supplementary Figures 1A, B and Figures 1F, G). We observed a significantly lower percentage of stained p35+Ebi3+ cells compared to GFP+TdTomato+ B cells in each condition, but the trends in IL-35+ cell frequency was conserved between the two approaches. The addition of αCD40 stimulation to αIgM was able to induce a population of IL-35+ Breg cells (Figures 1F, G). Furthermore, BCR + LPS stimulation with or without CD40 co-stimulation was able to induce IL-35 expression in regulatory B cells, but not CD19+CD1dLoCD21LoCD5- Bcon cells (Figures 1F, G). This suggests that ex vivo, BCR stimulation promotes the expansion of double positive regulatory cells only in the presence of either TLR and/or CD40 stimulation, and that inhibition of either the BCR or TLR signaling arm could potentially impact regulatory B cell function.
We also asked whether B cells derived from healthy mice would have the same responses to different stimuli than B cells from a tumor-bearing animal. When exposed to various stimuli ex vivo, splenic B cells derived from tumor-bearing mice generated a somewhat greater frequency of GFP+TdTomato+ B cells in response to LPS or αIgM + αCD40 compared to tumor-naïve B cells (Figure 1H and Supplementary Figure 1C). Examination of GFP and TdTomato signal intensity revealed that the increased frequency is due to GFP signal and not TdTomato (Figure 1I). The only instance we observed a significant increase in TdTomato signal intensity due to the influence of PDAC was with LPS + αCD40 stimulation, which did not have an increase in GFP intensity (Figure 1I). This is very interesting because Ebi3 expression is more highly regulated of the two subunits, but PDAC is primarily modulating Il12a expression in the spleen. Overall, this does suggest that pancreatic tumors can pre-condition B cells to have augmented IL-35 responses to stimuli.
BCR Engagement Promotes PDAC Growth and Represses Effector T Cell Responses
To understand the role of BCR activation in pancreatic tumor growth in vivo, we utilized the MD4 mouse model (30). B cells derived from MD4 mice have a fixed transgenic BCR that is reactive to hen egg lysozyme, therefore incapable of engaging host or pancreatic tumor-generated antigens. Autoimmune studies in MD4 mice linked disease remittance and inflammatory T cell activity to decreased suppressive abilities of B cells, thus we hypothesized that the inability to sense antigen may affect the suppressive abilities of B cells in PDAC (8). Tumor-naïve MD4 mice did not show any significant alterations in the development of immune cell populations or B cell subsets compared to WT mice (Supplementary Figures 2A, B). Syngeneic PDAC cells orthotopically injected into the pancreas of MD4 mice formed significantly smaller tumors than WT controls, indicating that BCR signaling can promote tumor growth (Figure 2A). Examination of the B cell infiltrate revealed a significant decrease in the frequency of IL-35+ and IL-10+ B cells in MD4 mice as compared to WT, indicating that BCR signaling in response to PDAC induces an immunosuppressive B cell phenotype in vivo (Figure 2B and Supplementary Figure 3A). The frequency of the Breg cell subset was not altered in MD4 tumors suggesting that the BCR signaling plays a role in regulatory B cell function and not expansion of the regulatory B cell subset in response to PDAC (Figure 2C). We recently showed that decrease in IL-35 production by B cells in PDAC results in significantly increased effector T cell responses coupled with decreased tumor growth (6). Consistent with this, we detected an increased frequency of IFNγ+CD4+ and IFNγ+CD8+ T cells in MD4 tumors coupled with a decreased frequency of Foxp3+ T cells compared to WT tumoral immune infiltrates (Figure 2D and Supplementary Figure 3B). Furthermore, we observed increased effector CD4+ and CD8+ T cell activity in both the spleen and draining lymph nodes of MD4 mice indicating the presence of a systemic response in secondary lymphoid organs as well (Figures 2E, F and Supplementary Figures 3C, D). To understand if decreased tumor growth in MD4 mice was dependent on T cell responses, we depleted CD4+ or CD8+ T cells (Figure 2G). Depletion of CD4+ T cells from MD4 mice resulted in significantly increased pancreatic tumor growth compared to control MD4 tumors, whereas depletion of CD8+ T cells showed a non-significant trend towards increase in tumor growth compared to isotype controls (Figure 2H). Collectively these data suggest that BCR activation promotes pancreatic cancer growth largely through suppression of CD4+ T cell responses in vivo.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/745873/fimmu-12-745873-HTML/image_m/fimmu-12-745873-g002.jpg)
Figure 2 B cell receptor activation promotes pancreatic tumor growth, B cell IL-35 expression, and effector T cell activity. (A) Measured volumes and weights of orthotopic KPC4662 tumors from WT and MD4 mice 21 days post-tumor injection. N=18. (B) Frequencies of IL-35+ (p35+) and IL-10+ B cells isolated from WT and MD4 KPC4662 orthotopic tumors at day 21 post-tumor injection. N=3. (C) Frequency (left) and number (right) of CD19+CD21HiCD1dHiCD5+ B cells isolated from WT and MD4 KPC4662 orthotopic tumors at day 21 post-tumor injection. N=6. (D) Frequencies of and Foxp3+ and IFNγ+ CD4+ or CD8+ T cells isolated from KPC4662 orthotopic tumors of WT and MD4 tumor-bearing mice at day 21 post-tumor injection. N=3. (E) Frequencies of IFNγ + and TNFα+ CD4+ or CD8+ T cells isolated from spleens of WT and MD4 KPC4662 orthotopic tumor-bearing mice at day 21 post-injection. N=3. (F) Frequencies of IFN γ + and TNFα+ CD4+ or CD8+ T cells isolated from tumor-draining lymph nodes (TdLN) of WT and MD4 KPC4662 orthotopic tumor-bearing mice at day 21 post-injection. N=3. (G) Schematic of T cell depletion in MD4 mice with KPC4662 orthotopic tumors and representative flow cytometry plots of CD45+ splenocytes isolated at day 21 post-tumor injection. (H) Measured volumes and weights of orthotopic KPC4662 tumors from WT and MD4 mice with or without CD4+ or CD8+ T cell depletion 21 days post-tumor injection. N=5. Error bars indicate SEM; p values in (A–F) were calculated using unpaired t-test; p values in (H) were calculated using one-way ANOVA. NS, non-significant, *p<0.05, **p<0.005, ***p<0.001. Data are representative of three independent experiments. Experiments were performed using 7–8-week-old mice of indicated genotypes.
MyD88 Signaling in B Cells Does Not Impact IL-35+ Breg Frequency In Vivo
TLR-mediated B cell suppression of T cell responses in models of bacterial infection and autoimmune disease has been linked mainly to expression of IL-10 (10, 31), and is highly dependent on signaling through the TLR adapter protein MyD88. Recently, TLRs have also been discovered to induce IL-35 expression in B cells in the same disease models (12), leading us to question the roles of TLRs on B cells in PDAC progression and T cell suppression. To understand the role of B cell specific TLR activation in pancreatic tumor growth in vivo, we generated CD19Cre/+; Myd88Fl/Fl mice. Healthy mice with Myd88 recombination specifically in B cells did not have altered immune cell development or B cell development, congruent with a similar Cd79aCre; Myd88Fl/Fl model (Supplementary Figures 1A, B) (32). Orthotopic pancreatic tumors in CD19Cre/+; Myd88Fl/Fl mice displayed trends towards decreased growth as compared to CD19+/+; Myd88Fl/Fl controls, but the change in growth was less robust than the tumor growth change displayed in MD4 mice (Figures 3A, 2B). We observed that IL-10 expression in B cells was significantly downregulated in CD19Cre/+; Myd88Fl/Fl mice (Figure 3B). However, expression of neither IL-35 subunits was affected, indicating that B cell expression of IL-10 is at least partially dependent on MyD88 in PDAC, but IL-35 could primarily be regulated by another antigenic pathway, most likely the BCR (Figure 3B). The loss of MyD88 signaling in B cells resulted in increased CD8+ T cell effector responses in the spleen but not in the tumor nor the tumor-draining lymph node upon tumor challenge, suggesting insufficient T cell priming (Figures 3C–E). Inversely, CD4+ T cell effector responses were significantly increased in the tumor and tumor-draining lymph node of CD19Cre/+; Myd88Fl/Fl mice and unaltered in the spleen (Figures 3C–E). Collectively these data demonstrate that MyD88 signaling in B cells does not play a prominent role in promoting PDAC growth.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/745873/fimmu-12-745873-HTML/image_m/fimmu-12-745873-g003.jpg)
Figure 3 PDAC influences B cell IL-35 expression in vivo in a MyD88-independent manner. (A) Measured volumes and weights of orthotopic KPC4662 tumors from Cd19+/+; Myd88Fl/Fl (WT) and Cd19Cre/+; Myd88Fl/Fl (Myd88Fl/Fl) mice at day 21 post-tumor injection. N=17. Cumulative of 3 individual experiments. (B) RT-PCR analysis of Il12a, Ebi3, and Il10 expression from CD19+ B cells derived from tumors and spleens of Cd19+/+; Myd88Fl/Fl and Cd19Cre/+; Myd88Fl/Fl mice at day 21 post-tumor injection. N=5. (C) Intracellular flow cytometry of IFNγ and TNFα in CD8+ and CD4+ T cells derived from tumors of Cd19+/+; Myd88Fl/Fl (WT) and Cd19Cre/+; Myd88Fl/Fl (MyD88Fl/Fl) KPC4662 orthotopic tumor-bearing mice at day 21 post-tumor injection. N=4. (D) Intracellular flow cytometry of IFNγ and TNFα in CD4+ and CD8+ T cells derived from spleens of Cd19+/+; Myd88Fl/Fl and Cd19Cre/+; Myd88Fl/Fl KPC4662 orthotopic tumor-bearing mice at day 21 post-tumor injection. N=4. (E) Intracellular flow cytometry of IFNγ and TNFα in CD4+ and CD8+ T cells derived from spleens of Cd19+/+; Myd88Fl/Fl and Cd19Cre/+; Myd88Fl/Fl KPC4662 orthotopic tumor-bearing mice at day 21 post-tumor injection. N=4. Error bars indicate SEM; p values in (A–D) were calculated using t-test; p values in (B) were calculated using two-way ANOVA. NS, non-significant, *p<0.05, **p<0.005, ***p<0.001 Experiments were performed using 7–8-week-old mice of indicated genotypes.
BCR Pathway Kinase Inhibition Represses IL-35 Expression
Since BCR signaling was dominant in controlling pancreatic tumor growth in vivo, we examined how key signaling molecules in the BCR signaling pathway could contribute to IL-35 expression. To understand how the BCR could affect IL-35 expression we performed a screen with small molecule inhibitors of the BCR signaling pathway using primary IL-35-reporter B cells isolated from Il12aGFP/GFP; Ebi3Tom/Tom spleens (Figure 4A). Inhibitor efficacy was assessed by Western blotting (Supplementary Figure 4A) and analyzed for reporter fluorescence via flow cytometry. We reasoned that the effector pathways that may most significantly affect IL-35 production may include pathways that promote cross talk between BCR and CD40 signaling, mimicking B cell and T cell interactions in vivo. Protein Kinase D (PKD) is a serine/threonine kinase that in response to BCR stimulation promotes synergy between the BCR and CD40 in a TRAF2-dependent manner (33, 34). PKD is required in B cells for complete proliferation, IL-6, TNFα, and IgM secretion responses from antigen and CD40 co-stimulation (33, 34). We tested the role of PKD in IL-35 expression as well as upstream BCR kinases Bruton’s tyrosine kinase (Btk), Phosphoinositide 3-kinase delta (PI3Kδ), and Lyn. Additionally, we tested key downstream signaling effectors: extracellular signal-regulated kinase 1/2 (ERK1/2), Nuclear factor of activated T-cells (NFAT), and nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB) in immunosuppressive cytokine production. We chose to investigate the roles of these effectors in the context of not only BCR stimulation, but also TLR (LPS) stimulation as well, since recent work demonstrated that effective TLR signaling activates and requires components of the BCR signaling pathway to fully signal (35).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/745873/fimmu-12-745873-HTML/image_m/fimmu-12-745873-g004.jpg)
Figure 4 Identification of BCR pathway signaling molecules that govern IL-35 expression. (A) Heatmap displaying frequency of GFP+TdTomato+ B cells from all stimulation conditions and inhibitors 48 hours after stimulation/inhibition. N=3 per condition. Experiment repeated 2 times. (B) Heatmap displaying frequency of viable B cells from all stimulation conditions and inhibitors 48 hours after stimulation/inhibition. N=3 per condition. Experiment repeated 2 times. (C) Heatmap displaying TdTomato mean fluorescent intensity (MFI) from all stimulation conditions and inhibitors 48 hours after stimulation/inhibition. N=3 per condition. Experiment repeated 2 times. (D) Heatmap displaying GFP mean fluorescent intensity (MFI) from all stimulation conditions and inhibitors 48 hours after stimulation/inhibition. N=3 per condition. Experiment repeated 2 times. (E) Western immunoblots of WT splenic B cells treated with DMSO, BTKi (5uM), or PKDi (10uM) for 30 minutes followed by stimulation with aIgM (10µg/mL). 30ug of total protein loaded per lane. (F) Quantification of immunoblots in (E).
We observed significant decreases in the frequency of GFP+TdTomato+ B cells when several kinases were individually inhibited (Figure 4B). Additionally, we observed significant changes in the MFI of GFP and TdTomato reporter signal with different inhibitors (Figures 4C, D). Inhibition of the BCR signaling phosphatase Lyn did not result in significant alteration of GFP or TdTomato intensity in response to all stimuli and significantly increased activity in most instances where BCR (αIgM) stimulation was used (Figures 4B–D and Supplementary Figures 6D, 7D). This is likely due to Lyn’s role in the negative regulation of the BCR signaling pathway and is consistent with previous data showing expansion of immunosuppressive IL-10 producing B cells in Lyn-/- mice (36). Inhibition of BCR signaling kinases PKD, PI3Kδ and Btk, on the other hand, resulted in the most significant downregulation of GFP+TdTomato+ cell frequency and TdTomato signal intensity in almost all stimulation conditions (Figures 4B–D and Supplementary Figures 5A-C, 6A-C). In conditions containing αIgM stimulation with these kinase inhibitors, PKDi was the most efficient at reducing GFP+TdTomato+ cell frequency and TdTomato signal intensity (Figures 4B–D and Supplementary Figures 5A-C, 6A-C). At the same time, Btk or PI3Kδ inhibition was more toxic than PKD inhibition in αIgM stimulation conditions (Figure 4E and Supplementary Figures 4B, C). This is likely related to the wider range of signaling pathways supporting survival that are activated downstream of Btk or PI3Kδ than those downstream of PKD. Overall, these observations suggest that PKD inhibition may be the most promising therapeutic target for IL-35 inhibition.
With regards to the signaling effector pathways downstream of BCR and CD40 signaling, the most prominent effect on reporter signal intensity was seen with inhibition of NF-κB signaling. Inhibition of IKKα/β (NF-κBi) resulted in near complete abolition of both GFP and Tomato signals regardless of stimulation conditions (Supplementary Figures 5E, 6E, 7E). This agrees with prior reports stating the requirement of NF-kB signaling in the expression of Ebi3 in B cells (37) and Il12a in antigen presenting cells (38). However, NF-kB inhibition in vitro over 48 hours resulted in high cellular toxicity regardless of the stimulation condition, thus making it an unlikely candidate for therapeutic use. Inhibition of ERK1/2 did not result in significant alterations in GFP or TdTomato reporter signal intensity (Supplementary Figures 6F, 7F). ERK1/2 inhibition did however decrease the frequency of GFP+TdTomato+ B cells leading us to conclude that ERK signaling does not play a role in regulating IL-35 expression, but rather aids in the proliferation of IL-35 expressing cells (Supplementary Figure 5F). This is consistent with the effect of MEK inhibition on IL-10+ regulatory B cell expansion (39). Inhibition of NFAT signaling, which regulates BCR-dependent IL-10 expression (40), resulted in significant decreases in GFP+TdTomato+ cell frequency primarily through inhibition of TdTomato reporter intensity (Supplementary Figure 5G, 6G). The regulation of Ebi3 transcription by NFAT is entirely plausible as there is an NFAT binding region in the Ebi3 promoter region but was determined to be non-essential for activation in BMDCs (37). We continued to evaluate the role of PKD in B cell signaling by examining the effects on ERK, NFAT, and NF-kB activation in comparison to Btk inhibition, which achieves similar IL-35 reporter inhibition (Figure 4F). Inhibition of Btk with Ibrutinib and PKD with CRT0066101 (41) resulted in significant inhibition of ERK, NFAT, and NF-kB signaling to relatively similar extents (Figure 4G). Collectively, this data points to IL-35 expression, with inhibition of PKD pathway being the most efficient in shutting down expression of IL35 and maintaining cellular viability.
Protein Kinase D2 Function in B Cells Promotes IL-35 Expression and PDAC Growth
We focused our efforts on further evaluating the role of PKD in IL-35 expression by B cells. CRT0066101 (PKDi) inhibits all 3 PKD isoforms, but the relative contributions of each isoform on IL-35 expression is unknown. In B lymphocytes, the PKD2 isoform is predominantly expressed (42) whereas PKD3 is expressed at lower levels and PKD1 expression was not detected by RT-PCR analysis (Figure 5A). Next, we aimed to understand the relative contributions of PKD2 and PKD3 in B cell signaling. To do this, we used B cells isolated from a mouse model harboring a knock-in of S707A and S711A mutations into the WT Prkd2 locus, which result in a lack of PKD2 catalytic activity (42). Therefore, B cells derived from this model only have active PKD3 signaling. PMA stimulation of homozygous Prkd2-S707A/S711A splenic B cells resulted in an almost complete absence of pan-PKD serine phosphorylation in the catalytic domain leading us to conclude that PKD2 and not PKD3 is the dominant active PKD isoform in primary B cells (Figure 5B). PKD signaling is not only active in B cells but is also active in Kras mutant pancreatic cancer cells (43). Prior studies revealed PKD function significantly impacts pancreatic cancer cell growth and metastatic potential (41, 44). RT-PCR analysis revealed that KPC cells express all PKD isoforms, predominantly PKD1, which is consistent with other PDAC cell lines (41, 45, 46) (Figure 5C). To understand the role of PKD signaling in KPC4662, we used CRT0066101 and performed MTT proliferation assays in vitro (Figure 5D). Consistent with other PDAC lines previously tested (41), KPC cell proliferation is significantly inhibited with loss of PKD function (Figure 5D).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/745873/fimmu-12-745873-HTML/image_m/fimmu-12-745873-g005.jpg)
Figure 5 Protein Kinase D regulates PDAC growth in immunocompetent mice. (A) RT-PCR analysis of PKD isoform expression in WT splenic B cells. ND=not detectable. N=3. (B) Phospho-PKD (S744/748) immunoblot of PMA-stimulated splenic B cells that are WT (+/+), heterozygous (M/+), or homozygous (M/M) for S707A/S711A mutation in Prkd2. (C) RT-PCR analysis of PKD isoform expression in KPC4662 cells. N=3. (D) MTT proliferation assay of DMSO or CRT0066101 treated KPC4662 cells in vitro. N=5. (E) Treatment schema and study design (left). Briefly, 10 days post-KPC4662 injection, WT mice were stratified into groups and treated with CRT0066101 or vehicle QOD by oral gavage. Mice were treated for 2 weeks before sacrifice at day 24 post-tumor injection. Orthotopic KPC4662 tumor volumes of CRT0066101 and vehicle treated mice as measured by ultrasound imaging (right). (F) Frequencies of indicated CD45+ immune populations isolated from tumors (top) and cell number per mm3 of tumor (bottom) from vehicle and CRT0066101 treated mice. (G) Frequencies of indicated CD45+ immune populations isolated from tumors (top) and cell number (bottom) from vehicle and CRT0066101 treated mice. Error bars indicate SEM; p values in (A, E) were calculated using unpaired t-test; p values in (C) were calculated using one-way ANOVA; p values in (D, F, G) were calculated using two-way ANOVA. NS, non-significant, *p<0.05, **p<0.005, ***p<0.001 ****p<0.0001. Experiments were performed using 7–8-week-old mice of indicated genotypes.
In vivo, we observed a significant delay in orthotopic tumor growth of immunocompetent mice treated with CRT0066101 inhibitor (Figure 5E). However, this experiment could not discern the potential contribution of PKD inhibition in cancer cells versus immune cells on tumor growth. To parse apart the effects of immune PKD function from cancer-cell autonomous PKD function in PDAC, we orthotopically injected KPC4662 cells into Prkd2-S707A/S711A mice. We observed no significant changes in tumor growth in the absence of host PKD2 function (Supplementary Figure 8A) suggesting that cancer cell PKD function is the key driver of PDAC tumor growth. However, we hypothesized that inhibition of PKD function in T cells, may have hindered anti-tumor immunity, resulting in increased tumor growth. PKD contributes to the expression of the effector cytokine IFNγ in response to TCR stimulation (42, 47). We also observed that activated Prkd2-S707A/S711A CD4+ and CD8+ T cells have significantly reduced IFNγ production ex vivo suggesting that T cell PKD function also plays a significant role in tumor growth (Supplementary Figures 9A, B). Regarding B cells, examination of the tumor immune infiltrate in both PKDi and Prkd2-S707A/S711A tumors by flow cytometry revealed that the frequency of infiltrating B cells was significantly decreased (Figure 5F) (Supplementary Figure 8B). Inversely, we observed a substantial increase in the frequency of CD11b+ myeloid cells upon PKD pharmalogical or genetic inhibition (Figure 5F) (Supplementary Figure 8B). Suppression of PKD in myeloid cells prevents pro-inflammatory cytokine expression in response to TLR9 agonists (48) and bacterial infection (49), but the effects in response to cancer are not known. CD4+ and CD8+ T cell frequencies were not decreased with PKDi, but CD19+ B cell frequencies and overall cell number were significantly reduced (Figure 5G). We further examined if the decreased tumoral B cell infiltrate was attributable to a specific B cell subset and found that PKD inhibition primarily decreased the frequency of follicular (CD19+CD93-CD21LoIgDHi) B cells within the tumor (Supplementary Figure 10A). Collectively, this data suggests that the decrease in tumor growth is due to inhibition of cancer-cell autonomous PKD, and the role of B cell and/or T cell PKD signaling in PDAC tumor immunity needs to be further understood in vivo.
To elucidate the effects of PKD2 function on IL-35 specifically in the regulatory B cell subset, we stimulated WT and Prkd2S707A/S711A CD19+CD1dHiCD21HiCD5+ Bregs in vitro and analyzed p35 and EBi3 expression by flow cytometry (Figure 6A). Deficiency in PKD2 function led to a significant decrease in the frequency of p35+EBi3+ (IL-35+) Breg cells, confirming that PKD2 is a critical regulator of IL-35 expression in regulatory B cells (Figure 6B). Furthermore, we tested how deficiency in PKD2 function affected the expression of other IL-12 family cytokines, IL-12 and IL-27. Il27a and Il12b are not widely expressed by regulatory B cells (5). PKD2 functional deficiency does not enhance expression of Il12b under varying stimulation conditions, but it does significantly enhance Il27a under BCR and CD40 stimulation suggesting that PKD2 represses IL-27 production by regulatory B cells in certain conditions (Supplementary Figure 11A). We also tested the expression of other immunosuppressive cytokines, Il10 and Tgfb1 in the context of PKD2 functional deficiency. While Tgfb1 expression is not altered with PKD2 function, Il10 expression was significantly enhanced in certain stimulation conditions indicating a potential role for PKD2 in regulating IL-10 as well (Supplementary Figure 11B). To test the functional ability of Prkd2S707A/S711A Bregs to suppress T cell function we established in vitro transwell co-culture assays, containing activated WT and Prkd2-S707A/S711A CD19+CD1dHiCD21HiCD5+ Bregs and activated WT CD4+ and CD8+ T cells (Figure 6C). Prkd2S707A/S711A regulatory B cells were significantly less effective at suppressing effector T cell cytokine expression in both CD4+ and CD8+ T cells (Figure 6C).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/745873/fimmu-12-745873-HTML/image_m/fimmu-12-745873-g006.jpg)
Figure 6 Inhibition of B cell PKD function suppresses IL-35 expression and tumor growth. (A) Representative flow cytometry plots of intracellular EBi3 and p35 in WT and homozygous Prkd2S707A/S711A CD19+CD21HiCD1dHiCD5+ and CD19+CD21LoCD1dLoCD5- splenic B cells stimulated for 48 hours with LPS + αIgM + αCD40 ex vivo. N=4. (B) Frequency of p35+EBi3+ B cells from (A). N=4. (C) Intracellular flow cytometry of IFNγ and TNFα from CD4+ and CD8+ T cells co-cultured in contact-independent transwells with either WT or homozygous Prkd2S707A/S711A CD19+CD21HiCD1dHiCD5+ regulatory B cells. N=3. (D) Orthotopic KPC4662 tumor volumes of µMT mice reconstituted with DMSO or CRT0066101 (PKDi) treated WT splenic B cells. N=9. Cumulative of 2 independent experiments. (E) Quantification of intratumoral CD19+ B cells per mm3 of tumor by flow cytometry from mice in (D). N=5. (F) Quantification of intratumoral CD4+ T cells (left) and CD8+ T cells (right) per mm3 of tumor by flow cytometry from mice in (D). N=5. Error bars indicate SEM; p values in (B, C) were calculated using unpaired t-test; p values in (D–F) were calculated using one-way ANOVA; p values in (G) were calculated using two-way ANOVA. NS, non-significant, *p<0.05, **p<0.005, ***p<0.001, ****p<0.0001. Experiments were performed using 7–8-week-old mice of indicated genotypes.
To understand the role of B cell specific PKD2 function on tumor growth we reconstituted mice lacking mature B cells (µMT) with splenic B cells treated ex vivo with either CRT0066101 (PKDi) or DMSO followed by orthotopic injection of KPC cells. Mice reconstituted with PKDi-treated B cells had significantly smaller tumors than mice reconstituted with DMSO-treated B cells (Figure 6D). Examination of the tumor infiltrating lymphocytes revealed that the number of B cells in tumors was not altered by PKD inhibition (Figure 6E). However, PKD inhibition in B cells led to a significant increase in intratumoral CD4+ and CD8+ T cells as compared to DMSO controls (Figure 6F). Collectively this data suggests that PKD2 plays an essential role in not only IL-35 expression by B cells, but also effector T cell function and ultimately regulates balance of tolerogenic and anti-tumor immunity in pancreatic tumor growth.
Discussion
We have identified BCR signaling as an important driver in establishment of IL35+ regulatory B cell function in pancreatic cancer. Our data suggests that while both BCR and MyD88 signaling in B cells promote IL-35 expression in vitro, BCR signaling plays a dominant role in IL-35 expression in vivo in response to pancreatic cancer. We previously observed the immunological effects of targeting IL-35 in PDAC and concluded that IL-35 inhibition potently increases effector T cell responses and reduces tumor growth (5, 6). In comparison, inhibited BCR signaling in vivo (Figure 2A) results in similar potent increases in effector T cell responses and decreased tumor growth (Figures 2C–F). Inhibition of MyD88 signaling however, did not lead to potent T cell suppression both in vivo and in vitro (Figures 3C–F). This enforces the fact that B cell suppression of T cells in PDAC is primarily BCR-mediated.
One possibility is that MyD88 is not the key regulator of IL-35 expression in response to TLR stimuli. While all TLRs except TLR3 signal through MyD88, TLR4 also signals through TRIF to mediate downstream gene expression that is dependent on IRF3. Trif-/- murine dendritic cells express significantly less Ebi3 after LPS stimulation than WT controls (50), so B cells may also require TRIF for Ebi3 expression. TLR stimulation also requires components of the BCR signaling pathway to fully transduce gene expression as well. For example, Il10 and Il6 expression was significantly decreased in Syk-deficient B cells after stimulation with LPS demonstrating that TLR4 signaling is not dependent on MyD88 alone (35). Activation of the BCR pathway kinases through TLR stimulation may explain the role of TLR signaling in response to pancreatic cancer. Additionally, the BCRs ability to capture antigens could be acting as a priming mechanism for TLR signaling. Intracellular TLR responses are highly dependent on BCR activity if nothing more than to shuttle TLR antigens such as DAMPs to TLR-containing endosomes. BCR recognition of apoptotic cells is required to induce-TLR9 mediated IL-10 expression in B cells and enhances TLR4-mediated IL-10 expression (15). We primarily focused on TLR4, but the role of intracellular TLR stimulation in B cell-mediated immunosuppression will be important to investigate in the future.
Recently a few studies have looked at the role of BCR signaling proteins to understand their relevance in regulatory B cell biology. Inhibition of the kinase MEK in BCR stimulated B cells inhibited the development of regulatory B cells, but the effect on suppressive cytokine production is not known (39). Similarly, inhibition of Bruton’s tyrosine kinase (Btk), a central kinase in the BCR signaling pathway, inhibits regulatory B cell development and IL-35 mRNA expression in vitro (51). However, this finding was performed in the context of αCD40, IL-6, and IL-1β activating Btk so the role of the BCR in this process is still unknown. A key point of emphasis in future studies will be to understand the composition of the BCR in regulatory B cells. It is not yet known whether these BCRs are directed in a tumor-antigen specific manner or are more host-reactive in nature like the BCRs of marginal zone and B-1 B cells. A key regulator of IL-35 expression in vivo we did not examine is CD40. TLR stimulation alone produces significantly less IL-35 in vitro than TLR and CD40 co-stimulation (12, 18). Furthermore, B cell expression of CD40 is also required for protection in models of collagen-induced arthritis, ulcerative colitis, EAE, and salmonella infection (7–9, 12, 52). CD40 stimulation to regulatory B cells is derived from cognate CD4+ T cells in EAE (52), but it is not known in pancreatic cancer where antigen-specific T cell responses are not as prevalent. Deletion of CD40 in B cells harboring PDAC tumors, could theoretically result in two opposing results. Loss of CD40 in regulatory B cells would hypothetically inhibit their expansion, but loss of CD40 in effector B cell populations would inhibit their maturation as well.
We demonstrate that Protein Kinase D2 function in B cells is essential for expression of IL-35 and promotes pancreatic tumor growth. Pancreatic cancer cell-specific PKD function is important for cell proliferation, so the double-edged effect of inhibiting B cell and PDAC cell PKD makes it a promising target heading forward. Previous PDAC model studies in vivo utilized immunocompromised mice for evaluation of PKD inhibition (41). We have expanded that knowledge by inhibiting PKD function in immunocompetent PDAC models and discovered that tumor growth remains delayed. However, PKD inhibition seems to affect more than B cells in the immune system as T cell effector function is highly dependent on PKD function. Targeting PKD inhibitors to specific cell types will be key for future therapeutic approaches.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics Statement
The animal study was reviewed and approved by Division of Comparative Medicine, The University of North Carolina at Chapel Hill.
Author Contributions
DM and YP-G contributed to the conception and design of the study. DM, BM, and YP-G developed the methodology of the study. DM, BM, and CS contributed to the acquisition of all data, data quality assessment, and data analysis in the study. DM and YP-G contributed to the interpretation of the data. DM and YP-G were responsible for writing the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported in part by R37 CA230786 (YP-G), University Cancer Research Fund at the University of North Carolina at Chapel Hill, United States (YP-G); Concern Foundation Conquer Cancer Now Award (YP-G), F31 CA239494 (DM), and T32 CA71341-20 (DM). The UNC Flow Cytometry Core Facility and the UNC Lineberger Animal Studies Core are supported in part by P30 CA016086 Cancer Center Core Support Grant to the UNC LCCC.
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.
Acknowledgments
We thank Nancy Kren for assistance with sample harvesting. We thank Naim Rashid for discussions concerning statistical analysis. We thank the UNC animal studies core for assistance with CRT0066101 dosing. We thank the UNC Flow Cytometry Core Facility for assistance with cell sorting and technical support of data acquisition. Schematics in Figures 1A, D, and 2G, were made using BioRender.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2021.745873/full#supplementary-material
References
1. Siegel RL, Miller KD, Fuchs HE, Jemal A. Cancer Statistics, 2021. CA Cancer J Clin (2021) 71(1):7–33. doi: 10.3322/caac.21654
2. Johnson BA 3rd, Yarchoan M, Lee V, Laheru DA, Jaffee EM. Strategies for Increasing Pancreatic Tumor Immunogenicity. Clin Cancer Res (2017) 23(7):1656–69. doi: 10.1158/1078-0432.CCR-16-2318
3. Ho WJ, Jaffee EM, Zheng L. The Tumour Microenvironment in Pancreatic Cancer - Clinical Challenges and Opportunities. Nat Rev Clin Oncol (2020) 17(9):527–40. doi: 10.1038/s41571-020-0363-5
4. Leinwand J, Miller G. Regulation and Modulation of Antitumor Immunity in Pancreatic Cancer. Nat Immunol (2020) 21(10):1152–9. doi: 10.1038/s41590-020-0761-y
5. Mirlekar B, Michaud D, Searcy R, Greene K, Pylayeva-Gupta Y. IL-35 Hinders Endogenous Anti-Tumor T Cell Immunity and Responsiveness to Immunotherapy in Pancreatic Cancer. Cancer Immunol Res (2018) 6(9):1014–24. doi: 10.1158/2326-6066.CIR-17-0710
6. Mirlekar B, Michaud D, Lee SJ, Kren NP, Harris C, Greene K, et al. B Cell-Derived IL35 Drives STAT3-Dependent CD8(+) T-Cell Exclusion in Pancreatic Cancer. Cancer Immunol Res (2020) 8(3):292–308. doi: 10.1158/2326-6066.CIR-19-0349
7. Mauri C, Gray D, Mushtaq N, Londei M. Prevention of Arthritis by Interleukin 10-Producing B Cells. J Exp Med (2003) 197(4):489–501. doi: 10.1084/jem.20021293
8. Fillatreau S, Sweenie CH, McGeachy MJ, Gray D, Anderton SM. B Cells Regulate Autoimmunity by Provision of IL-10. Nat Immunol (2002) 3(10):944–50. doi: 10.1038/ni833
9. Mizoguchi A, Mizoguchi E, Takedatsu H, Blumberg RS, Bhan AK. Chronic Intestinal Inflammatory Condition Generates IL-10-Producing Regulatory B Cell Subset Characterized by CD1d Upregulation. Immunity (2002) 16(2):219–30. doi: 10.1016/S1074-7613(02)00274-1
10. Neves P, Lampropoulou V, Calderon-Gomez E, Roch T, Stervbo U, Shen P, et al. Signaling via the MyD88 Adaptor Protein in B Cells Suppresses Protective Immunity During Salmonella Typhimurium Infection. Immunity (2010) 33(5):777–90. doi: 10.1016/j.immuni.2010.10.016
11. Wang RX, Yu CR, Dambuza IM, Mahdi RM, Dolinska MB, Sergeev YV, et al. Interleukin-35 Induces Regulatory B Cells That Suppress Autoimmune Disease. Nat Med (2014) 20(6):633–41. doi: 10.1038/nm.3554
12. Shen P, Roch T, Lampropoulou V, O'Connor RA, Stervbo U, Hilgenberg E, et al. IL-35-Producing B Cells are Critical Regulators of Immunity During Autoimmune and Infectious Diseases. Nature (2014) 507(7492):366–70. doi: 10.1038/nature12979
13. Pylayeva-Gupta Y, Das S, Handler JS, Hajdu CH, Coffre M, Koralov SB, et al. IL35-Producing B Cells Promote the Development of Pancreatic Neoplasia. Cancer Discovery (2016) 6(3):247–55. doi: 10.1158/2159-8290.CD-15-0843
14. Collison LW, Workman CJ, Kuo TT, Boyd K, Wang Y, Vignali KM, et al. The Inhibitory Cytokine IL-35 Contributes to Regulatory T-Cell Function. Nature (2007) 450(7169):566–9. doi: 10.1038/nature06306
15. Miles K, Heaney J, Sibinska Z, Salter D, Savill J, Gray D, et al. A Tolerogenic Role for Toll-Like Receptor 9 is Revealed by B-Cell Interaction With DNA Complexes Expressed on Apoptotic Cells. Proc Natl Acad Sci United States America (2012) 109(3):887–92. doi: 10.1073/pnas.1109173109
16. Barr TA, Brown S, Mastroeni P, Gray D. TLR and B Cell Receptor Signals to B Cells Differentially Program Primary and Memory Th1 Responses to Salmonella Enterica. J Immunol (2010) 185(5):2783–9. doi: 10.4049/jimmunol.1001431
17. Yanaba K, Bouaziz JD, Matsushita T, Tsubata T, Tedder TF. The Development and Function of Regulatory B Cells Expressing IL-10 (B10 Cells) Requires Antigen Receptor Diversity and TLR Signals. J Immunol (2009) 182(12):7459–72. doi: 10.4049/jimmunol.0900270
18. Michaud D, Mirlekar B, Bischoff S, Cowley DO, Vignali DAA, Pylayeva-Gupta Y. Pancreatic Cancer-Associated Inflammation Drives Dynamic Regulation of P35 and Ebi3. Cytokine (2020) 125:154817. doi: 10.1016/j.cyto.2019.154817
19. Turnis ME, Sawant DV, Szymczak-Workman AL, Andrews LP, Delgoffe GM, Yano H, et al. Interleukin-35 Limits Anti-Tumor Immunity. Immunity (2016) 44(2):316–29. doi: 10.1016/j.immuni.2016.01.013
20. Kang R, Lotze MT, Zeh HJ, Billiar TR, Tang D. Cell Death and DAMPs in Acute Pancreatitis. Mol Med (2014) 20:466–77. doi: 10.2119/molmed.2014.00117
21. Hernandez C, Huebener P, Schwabe RF. Damage-Associated Molecular Patterns in Cancer: A Double-Edged Sword. Oncogene (2016) 35(46):5931–41. doi: 10.1038/onc.2016.104
22. Suthers AN, Sarantopoulos S. TLR7/TLR9- and B Cell Receptor-Signaling Crosstalk: Promotion of Potentially Dangerous B Cells. Front Immunol (2017) 8:775. doi: 10.3389/fimmu.2017.00775
23. Rubin SJS, Bloom MS, Robinson WH. B Cell Checkpoints in Autoimmune Rheumatic Diseases. Nat Rev Rheumatol (2019) 15(5):303–15. doi: 10.1038/s41584-019-0211-0
24. Jia X, Liu H, Xu C, Han S, Shen Y, Miao X, et al. MiR-15a/16-1 Deficiency Induces IL-10-Producing CD19(+) TIM-1(+) Cells in Tumor Microenvironment. J Cell Mol Med (2019) 23(2):1343–53. doi: 10.1111/jcmm.14037
25. Balachandran VP, Luksza M, Zhao JN, Makarov V, Moral JA, Remark R, et al. Identification of Unique Neoantigen Qualities in Long-Term Survivors of Pancreatic Cancer. Nature (2017) 551(7681):512–6. doi: 10.1038/nature24462
26. Ochi A, Nguyen AH, Bedrosian AS, Mushlin HM, Zarbakhsh S, Barilla R, et al. MyD88 Inhibition Amplifies Dendritic Cell Capacity to Promote Pancreatic Carcinogenesis via Th2 Cells. J Exp Med (2012) 209(9):1671–87. doi: 10.1084/jem.20111706
27. Pushalkar S, Hundeyin M, Daley D, Zambirinis CP, Kurz E, Mishra A, et al. The Pancreatic Cancer Microbiome Promotes Oncogenesis by Induction of Innate and Adaptive Immune Suppression. Cancer Discovery (2018) 8(4):403–16. doi: 10.1158/2159-8290.CD-17-1134
28. Bayne LJ, Beatty GL, Jhala N, Clark CE, Rhim AD, Stanger BZ, et al. Tumor-Derived Granulocyte-Macrophage Colony-Stimulating Factor Regulates Myeloid Inflammation and T Cell Immunity in Pancreatic Cancer. Cancer Cell (2012) 21(6):822–35. doi: 10.1016/j.ccr.2012.04.025
29. Carneiro-Lobo TC, Scalabrini LC, Magalhaes LDS, Cardeal LB, Rodrigues FS, Dos Santos EO, et al. IKKbeta Targeting Reduces KRAS-Induced Lung Cancer Angiogenesis In Vitro and In Vivo: A Potential Anti-Angiogenic Therapeutic Target. Lung Cancer (2019) 130:169–78. doi: 10.1016/j.lungcan.2019.02.027
30. Goodnow CC, Crosbie J, Adelstein S, Lavoie TB, Smith-Gill SJ, Brink RA, et al. Altered Immunoglobulin Expression and Functional Silencing of Self-Reactive B Lymphocytes in Transgenic Mice. Nature (1988) 334(6184):676–82. doi: 10.1038/334676a0
31. Lampropoulou V, Hoehlig K, Roch T, Neves P, Calderon Gomez E, Sweenie CH, et al. TLR-Activated B Cells Suppress T Cell-Mediated Autoimmunity. J Immunol (2008) 180(7):4763–73. doi: 10.4049/jimmunol.180.7.4763
32. Hou B, Saudan P, Ott G, Wheeler ML, Ji M, Kuzmich L, et al. Selective Utilization of Toll-Like Receptor and MyD88 Signaling in B Cells for Enhancement of the Antiviral Germinal Center Response. Immunity (2011) 34(3):375–84. doi: 10.1016/j.immuni.2011.01.011
33. Haxhinasto SA, Bishop GA. A Novel Interaction Between Protein Kinase D and TNF Receptor-Associated Factor Molecules Regulates B Cell Receptor-CD40 Synergy. J Immunol (2003) 171(9):4655–62. doi: 10.4049/jimmunol.171.9.4655
34. Haxhinasto SA, Bishop GA. Synergistic B Cell Activation by CD40 and the B Cell Antigen Receptor: Role of B Lymphocyte Antigen Receptor-Mediated Kinase Activation and Tumor Necrosis Factor Receptor-Associated Factor Regulation. J Biol Chem (2004) 279(4):2575–82. doi: 10.1074/jbc.M310628200
35. Schweighoffer E, Nys J, Vanes L, Smithers N, Tybulewicz VLJ. TLR4 Signals in B Lymphocytes are Transduced via the B Cell Antigen Receptor and SYK. J Exp Med (2017) 214(5):1269–80. doi: 10.1084/jem.20161117
36. Scapini P, Lamagna C, Hu Y, Lee K, Tang Q, DeFranco AL, et al. B Cell-Derived IL-10 Suppresses Inflammatory Disease in Lyn-Deficient Mice. Proc Natl Acad Sci United States America (2011) 108(41):E823–32. doi: 10.1073/pnas.1107913108
37. Wirtz S, Becker C, Fantini MC, Nieuwenhuis EE, Tubbe I, Galle PR, et al. EBV-Induced Gene 3 Transcription is Induced by TLR Signaling in Primary Dendritic Cells via NF-Kappa B Activation. J Immunol (2005) 174(5):2814–24. doi: 10.4049/jimmunol.174.5.2814
38. Grumont R, Hochrein H, O'Keeffe M, Gugasyan R, White C, Caminschi I, et al. C-Rel Regulates Interleukin 12 P70 Expression in CD8(+) Dendritic Cells by Specifically Inducing P35 Gene Transcription. J Exp Med (2001) 194(8):1021–32. doi: 10.1084/jem.194.8.1021
39. Yarchoan M, Mohan AA, Dennison L, Vithayathil T, Ruggieri A, Lesinski GB, et al. MEK Inhibition Suppresses B Regulatory Cells and Augments Anti-Tumor Immunity. PloS One (2019) 14(10):e0224600. doi: 10.1371/journal.pone.0224600
40. Matsumoto M, Fujii Y, Baba A, Hikida M, Kurosaki T, Baba Y. The Calcium Sensors STIM1 and STIM2 Control B Cell Regulatory Function Through Interleukin-10 Production. Immunity (2011) 34(5):703–14. doi: 10.1016/j.immuni.2011.03.016
41. Harikumar KB, Kunnumakkara AB, Ochi N, Tong Z, Deorukhkar A, Sung B, et al. A Novel Small-Molecule Inhibitor of Protein Kinase D Blocks Pancreatic Cancer Growth In Vitro and In Vivo. Mol Cancer Ther (2010) 9(5):1136–46. doi: 10.1158/1535-7163.MCT-09-1145
42. Matthews SA, Navarro MN, Sinclair LV, Emslie E, Feijoo-Carnero C, Cantrell DA. Unique Functions for Protein Kinase D1 and Protein Kinase D2 in Mammalian Cells. Biochem J (2010) 432(1):153–63. doi: 10.1042/BJ20101188
43. Liou GY, Doppler H, Braun UB, Panayiotou R, Scotti Buzhardt M, Radisky DC, et al. Protein Kinase D1 Drives Pancreatic Acinar Cell Reprogramming and Progression to Intraepithelial Neoplasia. Nat Commun (2015) 6:6200. doi: 10.1038/ncomms7200
44. Youssef I, Ricort JM. Deciphering the Role of Protein Kinase D1 (PKD1) in Cellular Proliferation. Mol Cancer Res (2019) 17(10):1961–74. doi: 10.1158/1541-7786.MCR-19-0125
45. Trauzold A, Schmiedel S, Sipos B, Wermann H, Westphal S, Roder C, et al. PKCmu Prevents CD95-Mediated Apoptosis and Enhances Proliferation in Pancreatic Tumour Cells. Oncogene (2003) 22(55):8939–47. doi: 10.1038/sj.onc.1207001
46. Ochi N, Tanasanvimon S, Matsuo Y, Tong Z, Sung B, Aggarwal BB, et al. Protein Kinase D1 Promotes Anchorage-Independent Growth, Invasion, and Angiogenesis by Human Pancreatic Cancer Cells. J Cell Physiol (2011) 226(4):1074–81. doi: 10.1002/jcp.22421
47. Navarro MN, Sinclair LV, Feijoo-Carnero C, Clarke R, Matthews SA, Cantrell DA. Protein Kinase D2 has a Restricted But Critical Role in T-Cell Antigen Receptor Signalling in Mature T-Cells. Biochem J (2012) 442(3):649–59. doi: 10.1042/BJ20111700
48. Park JE, Kim YI, Yi AK. Protein Kinase D1: A New Component in TLR9 Signaling. J Immunol (2008) 181(3):2044–55. doi: 10.4049/jimmunol.181.3.2044
49. Sutton JA, Rogers LM, Dixon B, Kirk L, Doster R, Algood HM, et al. Protein Kinase D Mediates Inflammatory Responses of Human Placental Macrophages to Group B Streptococcus. Am J Reprod Immunol (2019) 81(3):e13075. doi: 10.1111/aji.13075
50. Molle C, Nguyen M, Flamand V, Renneson J, Trottein F, De Wit D, et al. IL-27 Synthesis Induced by TLR Ligation Critically Depends on IFN Regulatory Factor 3. J Immunol (2007) 178(12):7607–15. doi: 10.4049/jimmunol.178.12.7607
51. Das S, Bar-Sagi D. BTK Signaling Drives CD1d(hi)CD5(+) Regulatory B-Cell Differentiation to Promote Pancreatic Carcinogenesis. Oncogene (2019) 38(17):3316–24. doi: 10.1038/s41388-018-0668-3
Keywords: regulatory B cell, IL-35, pancreatic cancer, PKD, BCR
Citation: Michaud D, Mirlekar B, Steward C, Bishop G and Pylayeva-Gupta Y (2022) B Cell Receptor Signaling and Protein Kinase D2 Support Regulatory B Cell Function in Pancreatic Cancer. Front. Immunol. 12:745873. doi: 10.3389/fimmu.2021.745873
Received: 22 July 2021; Accepted: 07 December 2021;
Published: 03 January 2022.
Edited by:
Peter Brossart, University of Bonn, GermanyReviewed by:
Ingunn Stromnes, University of Minnesota Twin Cities, United StatesAlessandra Fiore, Max Planck Institute of Biochemistry, Germany
Francesco De Sanctis, University of Verona, Italy
Copyright © 2022 Michaud, Mirlekar, Steward, Bishop and Pylayeva-Gupta. 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: Yuliya Pylayeva-Gupta, yuliyap1@email.unc.edu
†These authors have contributed equally to this work