Skip to main content

TRAF3 suppression encourages B cell recruitment and prolongs survival of microbiome-intact mice with ovarian cancer



Ovarian cancer (OC) is known for exhibiting low response rates to immune checkpoint inhibitors that activate T cells. However, immunotherapies that activate B cells have not yet been extensively explored and may be a potential target, as B cells that secrete immunoglobulins have been associated with better outcomes in OC. Although the secretion of immunoglobulins is often mediated by the microbiome, it is still unclear what role they play in limiting the progression of OC.


We conducted an in-vivo CRISPR screen of immunodeficient (NSG) and immune-intact wild type (WT) C57/BL6 mice to identify tumor-derived immune-escape mechanisms in a BRAC1- and TP53-deficient murine ID8 OC cell line (designated ITB1). To confirm gene expression and signaling pathway activation in ITB1 cells, we employed western blot, qPCR, immunofluorescent staining, and flow cytometry. Flow cytometry was also used to identify immune cell populations in the peritoneum of ITB1-bearing mice. To determine the presence of IgA-coated bacteria in the peritoneum of ITB1-bearing mice and the ascites of OC patients, we employed 16S sequencing. Testing for differences was done by using Deseq2 test and two-way ANOVA test. Sequence variants (ASVs) were produced in Qiime2 and analyzed by microeco and phyloseq R packages.


We identified tumor necrosis factor receptor-associated factor 3 (TRAF3) as a tumor-derived immune suppressive mediator in ITB1 cells. Knockout of TRAF3 (TRAF3KO) activated the type-I interferon pathway and increased MHC-I expression. TRAF3KO tumors exhibited a growth delay in WT mice vs. NSG mice, which was correlated with increased B cell infiltration and activation compared to ITB1 tumors. B cells were found to be involved in the progression of TRAF3KO tumors, and B-cell surface-bound and secreted IgA levels were significantly higher in the ascites of TRAF3KO tumors compared to ITB1. The presence of commensal microbiota was necessary for B-cell activation and for delaying the progression of TRAF3KO tumors in WT mice. Lastly, we observed unique profiles of IgA-coated bacteria in the ascites of OC-bearing mice or the ascites of OC patients.


TRAF3 is a tumor-derived immune-suppressive modulator that influences B-cell infiltration and activation, making it a potential target for enhancing anti-tumor B-cell responses in OC.


For many decades, ovarian cancer (OC) has consistently been ranked among the 10 most common and lethal cancers in women worldwide [1]. The majority of OC cases are classified as high-grade serous cancers, with approximately 50% of these cases involving alterations in homologous recombination (HR) repair genes [2]. Since most OC patients are diagnosed at stages III (51%) or IV (29%) [3], currently approved treatment modalities, such as surgery, chemotherapy, and/or targeted therapies (VEGF and PARP inhibitors, and the FRα antibody–drug conjugate), have contributed to an improvement in progression-free survival, but have not made a significant impact on overall survival [4]. The response rate to immunotherapies, including agents that block programmed cell death protein 1 (PD-1), its ligand PD-L1, or cytotoxic T-lymphocyte antigen 4 (CTLA-4), has also been low (5.9–12.3%) in OC patients [5, 6]. Therefore, delineating new immune targets for the treatment of OC patients is urgently needed.

While clinical studies focusing on the role of T cells in OC have yielded disappointing results, the role of B cells in the OC tumor microenvironment (TME) has gained increasing attention [7]. The infiltration of regulatory interleukin (IL)-10-producing B cells into the ascites of OC patients is correlated with poor overall survival [8], while infiltration of immunoglobulin-producing plasma cells (IgGs [9, 10] and IgAs [11]) has been positively correlated with improved patient survival. The importance of B cells in OC is thus gaining recognition, but it is still not fully understood how various factors affect the infiltration and activation of specific B cell populations in the TME.

It is known that B cells express IgA in response to antigen stimulation and costimulatory signals and that they function primarily against viruses and bacteria [12]. As such, bacteria–B cell interactions have been widely studied in various diseases (reviewed in [13]), including a study of breast cancer, in which bacteria coated with mucosal IgA were found to be significantly associated with estrogen levels [14]. In OC, a few studies have linked certain bacterial species to tumor progression in BRCA1-mutated OC patients [15, 16], although the specific mechanism remains unclear.

In this study, using an ID8 murine OC cell line model, genetically modified with a loss of function in the TP53 and BRCA1 genes, designated here as ITB1, we identified a new tumor-derived immune suppressive mediator, tumor necrosis factor receptor-associated factor 3 (TRAF3). Through an unbiased targeted genetic screen, we showed that TRAF3 acts as an immunomodulator that promotes immune escape by suppressing the anti-tumor activity of B cells. We also showed that the anti-tumor activity of B cells is intimately associated with the gut microbiota.


Cell line

A modified version of the ID8 mouse ovarian surface epithelial cell line, ID8-gTRP53-gBRCA1 (ITB1), was obtained from Dr. Iain A. McNeish [17] under an MTA agreement. This cell line was maintained at 37 °C in a humidified atmosphere with 5% CO2 in DMEM medium supplemented with 10% fetal bovine serum (FBS), 1% l-glutamine, 200 mM, and 100 units each of penicillin and streptomycin. Mycoplasma infection was monitored regularly by PCR, and antibiotics (De-Plasma, TOKU-E, D022) were added to the culture medium as needed.

gRNA pool library production

The murine ID8-gTRP53-gBRCA1-Cas9 cell line was generated by co-transfecting a lentiviral Cas9-Blast vector (Addgene #52,962) with the packaging plasmids pCMV-dR8.91 and pCMV-VSV-G (Addgene #8454) into HEK293T cells. Following 48 h of transfection, the virus was harvested and its titer was determined. The Cas9-Blast lentivirus was then stored at -80 °C until it was used for overnight infection of ID8-gTRP53-gBRCA1 cells. Selection with 5 μg/mL of blasticidin was performed two days later. To acquire clones with high Cas9 activity, single-cell sorting into 96-well plates was performed on the ID8-gTRP53-gBRCA1-Cas9 cells. Subsequently, infection with a lentivirus driving expression of a gRNA specific for CD274 and mCherry was carried out on multiple clones. Ten days post-infection, each clone was stimulated with 10 ng/mL of interferon (IFN)γ for 24 h, and the expression of PD-L1 was determined by FACS using an anti-CD274 antibody. The efficiency of Cas9 editing was determined by measuring the percentage of PD-L1 negative cells in the transduced population.

The mouse CRISPR Brie lentiviral pooled libraries consisting of 79,637 gRNAs were co-transfected with the packaging plasmids psPAX2 #12,260 and pCMVVSV-G #8454 into HEK293T cells. Six hours after transfection, the medium was replaced with a medium supporting virus production (DMEM supplemented with 20% of FBS). After 48 h, the lentiviral medium was harvested and stored at -80℃. Before use in the in-vivo experiments, the ID8-gTRP53-gBRCA1-Cas9 cells were infected with the Brie lentivirus at a multiplicity of infection (MOI) of 0.06 and selected with puromycin for at least 10 days before injection, as described below in the section ‘In-vivo experiments.’

Generating knockout cell lines

To create knockout lines, the genome-scale CRISPR-Cas9 knockout (GeCKO) protocol was followed [18]. Oligos were designed using the same sequences as those provided in the gRNA library pool (described in ' gRNA pool library production') and were phosphorylated and annealed with T4 PNK ligase. The lentiviral CRISPR plasmid was digested and dephosphorylated using the BsmBI restriction enzyme, and the annealed oligos were ligated into the lenti CRISPRv2 vector. The plasmid was transformed into stbl2 bacteria, and positive colonies were selected using ampicillin selection (100 µg/mL). The positive colonies were expanded, and plasmid DNA was isolated using the QIAGEN® Plasmid Plus Midi Kit. The plasmid DNA was used to create lentiviruses and to infect ITB1 cell line to generate new lines expressing the desired genetic modifications. To this end, two TRAF3 guides, 5-CAGGTTCACGTGCTGTACCG-3 and 5-CGGTACAGCACGTGAACCTG-3, were used to created two stable ITB1 TRAF3 knockout (KO) clones, designated TRAF3KO1 and TRAF3KO2, respectively.

Western blot analysis

All cells were harvested using a cell scraper and washed with ice-cold PBS. A lysis buffer supplemented with phosphatase inhibitor cocktails (BioTools, B15001A/B) and a protease inhibitor (Millipore Sigma, P2714-1BTL) was used to lyse the cells, which were then placed on ice for 30 min, followed by 3 min of ultrasonic cell disruption. The lysed cells were centrifuged at 14,000 rpm for 10 min at 4 °C, and the supernatant was collected. The protein concentration was determined using the Bradford assay (Bio-Rad, Protein Assay, cat# 500–0006), and 1 μg/μL of total cellular protein from each sample was used to run a PAGE and transferred onto a PVDF membrane (Bio-Rad, 1,704,157). The membranes were soaked in a blocking solution of 5% BSA (AMRESCO, 0332-TAM) in TBS-T [Tris-buffered saline (TBS)-Tween 20 (0.1%)] for 1 h and then hybridized overnight with primary antibodies diluted in 5% BSA and 0.1% Tween. The following day, the membranes were incubated with horseradish peroxidase (HRP)-conjugated secondary antibodies (1:20,000, Jackson), diluted in the blocking solution. A chemiluminescent reaction was developed using enhanced chemiluminescence (ECL) (Westar Nova 2.0 Cyanagen XLS071.0250), and images were captured using an Azure Biosystems camera system. For nuclear and cytoplasmic separation, the Rockland Nuclear & Cytoplasmic Extract Protocol [19] was followed. In brief, the cells were lysed with a low-salt lysis buffer, which allowed the nuclear proteins to remain intact while the cytoplasmic proteins were lysed. Then, using a high-salt lysis buffer, the nuclear membrane and the nuclear proteins were lysed. The separation quality was verified by using an antibody against histone H3, which should be present only in the nuclear fraction.

RNA isolation and real-time quantitative PCR

The isolation of total RNA and its conversion to cDNA were performed using the ISOLATE II RNA Mini Kit (Bioline, BIO-52073) and the qScript cDNA synthesis kit (Quanta Bioscience, 95,047–100), respectively, as per the manufacturers' protocols. A mix of cDNA and Fast SYBR qPCR Master mix (BioGate, EZ60) along with custom primer sets from IDT were used for qPCR analysis on a Roche light cycler 480 II machine.

Flow cytometry

Isolated cells were blocked with anti-CD16/32 (anti-FcγIII/II receptor, clone 2.4G2) for 10 min at 4 °C. Zombie Aqua™ fixable viability dye was then used to identify dead cells by incubating the cells with the dye for 15 min at room temperature. Surface markers were then stained by incubating the cells with the desired antibodies for 20–30 min at 4 °C. The stained cells were analyzed using a CytoFLEX Flow Cytometer, and the results were analyzed using the CytExpert software.

For intracellular cytokine staining, cells were first cultured for 5–6 h with PMA (25 ng/mL), ionomycin (1 μM) and brefeldin A (5 µg/mL) to increase the intracellular cytokine staining signals. Surface markers were then stained, and the cells were fixed using 4% paraformaldehyde for 20 min at room temperature. The cells were permeabilized using a permeabilization buffer (1 ×) (eBioscience, 00–8333-56), and intracellular staining was performed by incubating the cells for 20 min. The cells were then analyzed on a CytoFLEX flow cytometer.

In-vivo experiments

The in-vivo experiments were conducted using 6- to 8-week-old NSG mice (NOD.Cg-Prkdcscid Il2rgtm1Wjl/SzJ, Jackson labs) and wild-type (WT) C57/BL/6 mice (Envigo, Huntingdon, UK). The mice were housed in air-filtered laminar flow cabinets and provided with food and water ad libitum. Animal experiments were conducted in compliance with protocols established by the Institutional Animal Care and Use Committee (IACUC) of Ben-Gurion University of the Negev for ensuring animal welfare and minimizing discomfort. The animal ethical clearance protocol number used for the study was IL-23–05-2020(E).

Two different experiments were performed with the WT and NSG mice, as described in the Results section. In an experiment designed to investigate the role of anti-tumor immunity in the progression of ITB1-generated disease, WT and NSG mice were each injected intraperitoneally (i.p.) with 4–5 million ITB1 cells. In a separate experiment, which was designed to identify potential immunomodulators in ITB1 cells, WT and NSG mice were injected with 4–5 million ITB1 cells that expressed a gRNA library.

Two experiments were performed only with WT mice. In the first, designed to examine the involvement of the microbiome in OC, one group of WT mice was pretreated with an antibiotic cocktail (ABX) for 2 weeks to deplete the total microbiome before TRAF3KO cell injection. This was achieved by adding ampicillin (1 g/mL), vancomycin (0.5 g/mL), neomycin sulfate (1 g/mL), and metronidazole (1 g/mL) to the drinking water and changing the water every two days [20]. The supplementation of the drinking water ABX was continued throughout the entire experiment. Two other groups of WT mice were injected with ITB1 or TRAF3KO cells and received normal drinking water (vehicle) throughout the entire experiment. For the FACS experiments, WT mice were sacrificed 2 weeks after the injection, and the peritoneal cavity was washed out with phosphate-buffered saline (PBS) to collect the contents.

For all the survival experiments, the body weight of the mice was monitored once a week, and the mice were sacrificed when a 20% increase in body weight was observed, indicating the development of ascites. The ascites fluid was collected for the isolation of bacteria or tumor cells. The digestive and reproductive systems were also collected for histology analysis.

Bacterial DNA extraction

Bacterial DNA was isolated from ascites fluid or PBS wash using a modified phenol–chloroform protocol [21]. In brief, samples were lysed with a lysis buffer containing phenol, chloroform, and isoamyl alcohol (IAA) and then homogenized by bead-beating with bacterial disruption beads (RPI, 9833) for two cycles of 1 min each on a Precellys instrument. The aqueous and phenol phases were separated by centrifugation, and the aqueous phase was re-suspended in IAA in a sterile tube. The bacterial DNA was then concentrated and isolated using isopropanol precipitation and ethanol washes. To minimize contamination, sterile conditions were used, and negative controls were taken by sampling the environment and collecting tools.

Sorting of IgA + and IgA − bacteria

The above-described PBS washes were collected from co-housed mice and then centrifuged at 8,000 × g for 5 min at 4 °C. The pellet was then resuspended in 100 μL of a blocking buffer containing 20% normal rat serum and incubated for 20 min on ice. Thereafter, the samples were stained with 100 μL of staining buffer containing phycoerythrin (PE)-conjugated anti-mouse IgA (1:12.5; eBioscience clone mA-6E1) for 30 min on ice. Following this step, the samples were washed 3 times with 1 mL of staining buffer. The anti-IgA-stained bacteria were then incubated in 1 mL of staining buffer containing 50 μL of anti-PE magnetic activated cell sorting (MACS®) beads (Miltenyi Biotec) for 15 min at 4 °C. Finally, the samples were washed twice with 1 mL of staining buffer at 8,000 × g for 5 min at 4 °C, and sorted using MACS. The IgA-positive and the IgA-negative samples were then stored at -80 °C for future use.

Sequencing and sequence analysis

The genome-scale CRISPR-Cas9 library was sequenced on an Illumina MiSeq platform using a 250-bp paired-end chip. The sequences were mapped using a modified version of the caRpools R-package [22]. Differentially expressed guides were identified by loading the normalized gRNA count table into Deseq2 and determining the top differential genes based on the mean log2 fold change and false discovery rate (FDR). To identify significant pathways that were enriched or depleted in the screen, gene set enrichment analysis, including genes with a threshold of log2 fold change < (-1), was conducted using the ReactomePA [23] R package.

To amplify the V3-V4 region of the bacterial DNA isolated from the ascites of mice, a 2-step PCR-based protocol published by Holm et al. [24] was followed. The amplicons were then sequenced using the Illumina MiSeq platform (MiSeq Reagent Kit v2, 250-bp paired-end). Quality-checking of the demultiplexed sequencing reads was carried out using FastQC (Version 0.11.8) [25] and MultiQC (Version v1.7) [26]. Adapters, primers, and low-quality reads were removed using Fastp (Version 0.2, defaults parameters) [27]. The sequences were denoised by using Qiime2 version 2021.11 [28], and an amplicon sequence variant (ASV) table was created with taxonomic assignment using the SILVA database version 138, 99% NR [29]. Contaminant sequences were identified and removed based on control samples and sequence prevalence using the decontam package (Version 1.12) [30] and squeegee tool (Version 0.1.3) [31]. Further analysis of the ASV table was done using the MicroEco package (Version 0.11) [32] and phyloseq package (Version 1.38) [33] in R (Version 4.1.1).

A similar proof-of-concept analysis was performed for ascites from women with OC (see ‘Results section’).

Statistical analysis

The in-vitro experiments were repeated a minimum of two to three times, with representative data/images being presented in the ‘Results’ section. For the in-vivo experiments, a minimum of 5 mice per group were used. Statistical analysis was conducted using GraphPad Prism 9 software and results are presented as means ± SEM. For experiments with two or more groups, a two-way ANOVA with Tukey's multiple comparison tests was used. Significance levels of p-values of 0.05, 0.01, 0.001, or 0.0001 were calculated and are indicated in the figures by *, **, ***, or ****, respectively.


CRISPR screen identifies TRAF3 as an immunomodulator of BRCA1-mutated OC

To investigate the role of anti-tumor immunity in the progression of ITB1-generated disease in-vivo, we injected ITB1 cells i.p. into both immunocompromised NSG and immunocompetent WT mice and then monitored the survival of these tumor-bearing mice (Fig. 1A). The ITB1 cells did not form a single tumor mass, but rather caused a spreading disease in the peritoneum (as can be seen by immunohistochemistry staining), which also generated a significant accumulation of ascites (Fig. S1A). Although there was a significant difference in the rates at which the WT and NSG mice developed malignant disease, the median survival differed by only 8 days (41 vs. 49 days median survival). We, therefore, concluded that the ITB1 cells may have intrinsic immune escape mechanisms that suppress anti-tumor activity in WT mice, enabling rapid tumor progression.

Fig. 1
figure 1

Identification of TRAF3 as a potential immunomodulator. A Kaplan–Meier plot showing the survival of two groups of ITB1-injected mice: pink line—NSG mice and cyan line—WT mice. B Schematic representation of library preparation, injection, and library DNA isolation. C Volcano plot showing enriched guides in ITB1-injected WT mice (pink) compared to ITB1-injected NSG mice (cyan). Guides that are significant appear above the horizontal dotted line. D Reactome pathways enriched with the enriched guides in the NSG group. The size of the dots represents the number of genes related to the pathway, and the color represents the p-value

To identify potential immunomodulators in ITB1 cells that may be responsible for such immune escape, we knocked out a panel of immune-related genes by using a targeted guide RNA library of known immunomodulators [34] and then performed an in-vivo screen in WT and NSG mice (Fig. 1B). Specifically, we sought to identify genes that would facilitate the elimination of the tumor cells in the immunocompetent WT mice but not in the NSG mice. To this end, we first transduced the ITB1 cells with a lentiviral vector that expressed Cas9 and generated single-cell clones with the aim of identifying clones with a Cas9 editing efficiency of at least 95%. To test the editing efficiency, single clones overexpressing gRNA targeting PDL1 following IFNγ stimulation were stained for PDL1 (Fig. S1B). Clone 8, identified with an efficiency of 96.5%, was expanded and transduced with the complete gRNA library. The distribution of the library was verified by sequencing the amplicons produced from this ITB1 library-expressing clone. Out of 3878 guides, consisting of 1878 guides targeting 313 immune-related genes (six guides per gene) and 2000 guides as a control, only 500 guides were sequenced less than four times. After confirming that the majority of the guides were sequenced more than 10 times (Fig. S1C), we injected the gRNA library transduced ITB1 cells i.p. into NSG and WT mice. As soon as the production of ascites had resulted in a 20% increase in the weight of the mice, we drained the ascites, collected the cells, and isolated the total DNA for gRNA sequencing and analysis (Fig. 1B). By comparing the expression of gRNAs in the tumor cells growing in NSG and WT mice, we found six significantly differentially expressed genes in NSG mice compared to WT mice. Of these, two genes, NF1 and NF2, were significantly increased in the WT mice, and four genes, ARID4A, GM46125, OPRK1, and TRAF3, were significantly increased in the NSG mice (Fig. 1C). Notably, we included only those genes that were represented by two or more guides (out of the six represented in the library).

In view of our search for new tumor-derived targets that could sensitize ITB1 cells to anti-tumor immunity in mice, we focused on gRNAs and their signaling pathways that were increased in NSG mice compared to WT mice. We therefore conducted a Reactome pathway analysis on all genes whose absolute log2 fold change value was higher than one. The analysis revealed that TRAF3 was involved in six of the eight enriched pathways shown in Fig. 1D. Since TRAF3 had the highest expression difference and has been identified as a negative regulator of major histocompatibility complex class I (MHC-I) proteins and as being associated with a poor response to immune checkpoint blockade [35], we decided to further investigate its role in mediating immune escape in the ITB1-OC model.

TRAF3 depletion activates the type-I interferon pathway and MHC-I expression and enhances the immunogenicity of ITB1 tumors

To determine the effect of tumor-derived TRAF3 on immune escape, we created two stable ITB1 TRAF3 knockout (KO) lines, designated TRAF3KO1 and TRAF3KO2 (Fig. 2A). We first compared the proliferation and migration abilities of ITB1 and TRAF3KO1/2 cells in vitro and observed a slight reduction in the proliferation rate of TRAF3KO1/2 cells, but no differences in their ability to close the scratch in a wound-healing assay (Fig. 2B and S2A). As TRAF3 is known to regulate the NF-κB pathway [36], we compared the nuclear protein levels of p65 and RelB, the total levels of p100 and p52, and the mRNA levels of IkBa, IL6, and RelB in ITB1 and TRAF3KO1/2 cells. We observed no significant activation of the NF-κB pathway in TRAF3KO1/2 cells, except for a minimal activation of the non-canonical NF-κB pathway, as indicated by Relb expression (Fig. S2B and C). Since TRAF3 is also known to mediate type I IFN signaling [36], we compared the mRNA levels of several IFN-stimulated genes in ITB1 and TRAF3KO1/2 cell lines and found a significant increase in the two TRAF3KO cell lines in the levels of three major genes in the type-I IFN pathway, namely, IFN-β, ISG15, and IRF7 (Fig. 2C). We also found the upregulation of type-I IFN to be associated with activation of the stimulator of interferon genes (STING) pathway, as indicated by the nuclear levels of TBK1 and STING upon western blot analysis (Fig. 2D and Fig. S2D) and the nuclear levels of STING upon immunofluorescence analysis (Fig. 2E and Fig. S2E). The strong activation of type-I IFN in TRAF3KO1/2 cell lines and the known association of TRAF3 with MHC-I expression and the activation of tumor immunity prompted us to examine the levels of MHC-I expression on the surface of the tumor cells. Indeed, TRAF3KO1/2 cells expressed higher levels of MHC-I compared to ITB1 cells, and the expression of MHC-I on the surface of TRAF3KO1/2 cells was almost at the same as that of IFN-stimulated ITB1 cells (Fig. 2F).

Fig. 2
figure 2

TRAF3 knockout (TRAF3KO) cells activate type-I interferon (IFN)-related pathways. A Western blot of TRAF3KO clones demonstrating the depletion of TRAF3 compared to the wild-type (WT) TRAF3 of ITB1 lysates. B Growth of ITB1 and TRAF3KO cells in culture was monitored every 8 h until confluency. The cyan line represents the growth of ITB1 cells, and the pink lines represent the growth of TRAF3KO cells. C Quantitative PCR shows the mRNA levels of genes related to the type-I IFN pathway (IFNA, IFNB, IRF7, IFIT1, ISG15). Cyan and pink bars represent the levels in ITB1 and TRAF3KO cells, respectively. D Western blot of proteins in the STING pathway in ITB1 and TRAF3KO lysates. Actin and histone H3 were used as loading controls. Numbers in parentheses represent the protein size in kDa. E Immunofluorescence images showing the expression of STING (Cy3) in the nucleus (DAPI) of ITB1 and TRAF3KO cells. F Flow cytometry results showing the surface levels of MHC-I using an antibody conjugated to the FITC fluorophore. The mean fluorescence intensities (MFI) were compared between each condition and plotted on a bar graph. G Kaplan–Meier plot showing the survival of mice in four groups: NSG mice injected with ITB1 (solid pink line), NSG mice injected with TRAF3KO cells (dotted pink line), WT mice injected with ITB1 (solid cyan line), and WT mice injected with TRAF3KO cells (dotted cyan line)

We then validated the potential role of TRAF3 as a tumor-immune suppressive mediator by injecting TRAF3KO1 (henceforth referred to as TRAF3KO) and ITB1 cells into NSG and WT mice and monitoring ascites development and the survival of the mice. In NSG mice, there was no significant difference in the growth delay of tumors bearing TRAF3KO (33 days) or ITB1 (29 days) tumors, indicating that the tumors developed at similar rates in the absence of an active immune system. In contrast, in WT mice, the median survival of mice bearing ITB1 tumors was 38 days, while mice bearing TRAF3KO tumors survived for almost three times as long—a median of 98 days (Fig. 2G). These survival differences in WT but not NSG mice further support the hypothesis that TRAF3 acts an immune-suppressive modulator in ITB1 OC cells, most probably by controlling the levels of type-I IFN.

TRAF3 suppression in OC cells induces an IgA-mediated B cell immune response

To investigate the potential mechanisms by which immune cells eliminate or reduce the growth of TRAF3KO tumors, we injected ITB1 or TRAF3KO cells i.p. into WT mice, washed out the abdominal cavity two weeks after injection, and characterized the immune populations in the washing solution by flow cytometry. Analysis of the myeloid population, determined by the expression of CD45+CD11b+, showed a decrease in mice bearing TRAF3KO tumors compared to ITB1 tumors. Classification of the myeloid subsets stained with F4/80, CD11c, and LY6C/G showed that the percentage of macrophages (CD45+CD11b+F4/80+) was significantly lower in TRAF3KO tumors, while the percentages of other populations, such as dendritic cells (CD45+CD11b+CD11c+MHC-II+), monocytic cells (CD45+CD11b+Ly6C+Ly6G), and polymorphonuclear (CD45+CD11b+Ly6CLy6G+) myeloid-derived suppressor cells (MDSCs), were not significantly different between the two groups of mice (Fig. 3A). In the lymphocyte population, there were no differences in the percentages of CD3+, CD8+ T cells, and NK cells. However, we did observe increases in CD19+ B cells and CD4+ T cells (Fig. 3B).

Fig. 3
figure 3

TRAF3 knockout (TRAF3KO) cells alter the immune response. Profile of total myeloid cells (CD45+CD11b+), macrophages (CD45+CD11b+F4/80+), dendrites (CD45+CD11b+CD11c+MHCII+), monocytic (CD45+CD11b+Ly6C+Ly6G-), and polymorphonuclear (CD45+CD11b+Ly6C-Ly6G+) MDSCs in ITB1- (cyan) and TRAF3KO (pink)-injected mice. B Profile of total T cells (CD45+CD3+), B cells (CD45+CD19+), helper T-cells (CD45+CD3+CD4+), cytotoxic T-cells (CD45+CD3+CD8+), and NK cells (CD45+CD3-NK1.1+) in ITB1-injected (cyan) and TRAF3KO-injected (pink) mice. C Expression of immunoglobulin markers (IgA, IgD, IgM) on the surface of B cells in ITB1- and TRAF3KO-injected mice. D Volcano plot of a multiple t-test analysis comparing TRAF3KO and ITB1 washes. The plot shows the difference (x-axis) and p-value (y-axis) of the mean fluorescence intensity of each immunoglobulin scanned with an Ig isotyping array, E Heatmap showing that the ovarian cancer TCGA cohort is clustered according to the B cell pathway. F Survival curve showing the difference between patients with high (cyan) and low (pink) levels of IGHA1 in the TCGA data

To better characterize the B cell population in mice bearing TRAF3KO tumors, we compared the maturation of B cells by co-staining CD19 with the surface forms of IgM, IgA, and IgD. In these three immunoglobulin subclasses, there was significantly higher surface expression of IgA on B cells in the ascites of mice bearing TRAF3KO tumors compared to those bearing ITB1 tumors (Fig. 3C). Notably, there were no significant differences in the expression of the costimulatory protein CD40 or the percentage of regulatory B cells determined by the intracellular expression of IL-10 (Fig. S3A). To confirm the increased expression of CD19+IgA+ cells in mice bearing TRAF3KO tumors, we used an immunoglobulin isotyping array to compare the levels of seven different immunoglobulins in the peritoneum of mice injected with ITB1 or TRAF3KO cells. Three of the seven immunoglobulins were more abundant in the peritoneal cavities of TRAF3KO-injected mice than in those of ITB1-injected mice, with IgA being the most significantly expressed immunoglobulin (Fig. 3D and Fig. S3B and C).

To explore whether our findings have any clinical relevance in OC patients, we examined the gene expression profiles and clinical data of 380 high grade serous OC cases from the Cancer Genome Atlas (TCGA) portal. These profiles were subjected to single-sample gene set enrichment analysis against an immune gene set of 28 different immune cell types [37]. In this analysis, we visualized the associations between the TCGA samples and the different data sets and split the samples by k-means clustering, creating two distinct clusters (Fig. 3E and Fig. S3D). We noticed that the factor that differed to the greatest extent between the two clusters was the B cell signature (Fig. S3E) and that one of the major genes in the B cell signature was immunoglobulin heavy constant alpha 1 (IGHA1), which is part of both the monomeric form of IgA and the secretory dimeric IgA complex [38]. Comparing the survival probability of patients based on their B cell signature score showed that the median survival of patients with a high score was significantly longer than that of patients with a low score (Fig. S3F), and much longer when differentiating the patients by expression levels of IGHA1 (Fig. 3F). These results further support the hypothesis that the enrichment of the TME with B cells and secreted IgAs may reflect an enhanced overall anti-tumor immune response, which prolongs the survival of OC patients.

IgA-mediated B cell immune response is affected by the commensal microbiome

IgA-secreting B cells are key players in maintaining a symbiotic relationship between the host's immune system and the commensal microbiome [39]. Disruption of this symbiotic relationship and hence of the interaction between bacteria and B cells has been implicated in several diseases, including some cancers [40, 41]. To study the effect of bacteria on the anti-tumor B cells response, we depleted the microbiome by treating WT mice with a cocktail of four different antibiotics (ABX), as described in the Methods section. We verified the microbiome depletion by alterations in the cecum size [42, 43] and by monitoring the growth of colonies isolated from the feces of treated and untreated mice on blood agar plates (Fig. S4A and B). We then injected the ABX-treated mice with TRAF3KO cells and monitored their survival in comparison with untreated mice injected with ITB1 or TRAF3KO cells. ABX treatment shortened the median survival of TRAF3KO-injected mice to 35 days, compared to 92 days for untreated TRAF3KO-injected mice (Fig. 4A). To determine the response of the immune system to ABX treatment, we collected and profiled the immune cells from the abdominal cavities of ABX-treated TRAF3KO-injected mice, compared to those from the abdominal cavities of ITB1- or TRAF3KO-injected untreated mice. For TRAF3KO-injected mice, the comparison showed overexpression of total CD19+ cells, and particularly CD19+MHC-II+ cells, in the untreated mice but not in the ABX-treated mice (Fig. 4B). We then looked for associations between the abundance of IgA-coated bacteria and the expression of TRAF3. To this end, we isolated and purified IgA-coated bacteria from the peritoneal washes of both TRAF3KO-injected and ITB1-injected mice, as previously described [44]. We then applied 16S sequencing to characterize the taxa-specific coated bacteria of each group of mice. We found higher abundance of two genera of IgA-coated bacteria (Flavobacterium and Novosphingobium) in the ITB1 group and one genus (Sphingobium) in the TRAF3KO group (Fig. 4C and Fig. S4C). Finally, using the TCGA microbiome database [45], we looked for the same associations in humans regarding IGHA1 expression. We found two bacterial genera that correlated positively and two bacterial genera that correlated negatively with IGHA1 expression (Fig. 4D).

Fig. 4
figure 4

Bacteria are necessary for the B cell anti-tumor response. A Kaplan–Meier plot showing the survival of mice in four groups: untreated WT mice injected with ITB1 (cyan line), untreated WT mice injected with TRAF3KO cells (pink line), antibiotic (ABX)-treated WT mice injected with ITB1 (cyan dashed line), and ABX-treated WT mice injected with TRAF3KO cells (pink dashed line). B Bar graph representing the percentage of CD19 + cells and MHCII + cells in ITB1 (cyan circles), TRAF3KO (pink triangles), and ABX-treated TRAF3KO- (pink open triangles) injected mice. C Bar plot showing the abundance of each genus in each group. Correlation curves showing the positive correlation of Flammeovirga and Lachnoclostridium, and the negative correlation of Aeromonas and Lewinella, with IGHA1 expression

Overall, these results emphasize the contribution of bacteria in the TME to an enhanced (or otherwise) anti-tumor B cell response.

IgA-seq of ascites can distinguish ovarian cancer patients

To evaluate the association between microbial communities and IgA in human OC, we examined, as a proof of concept, 11 ascites samples from OC patients and 4 ascites samples from patients with cirrhosis. First, we isolated the total bacterial DNA from each sample, and then amplified and sequenced the 16S rRNA. After ruling out contaminations, as described in the Methods section, we analyzed the obtained sequencing data and estimated the diversity and abundance of microbial communities in each group. Using the Chao1 model, we did not observe significant differences in the α-diversity between OC and cirrhosis patients (Fig. 5A). An ordination plot based on Bray–Curtis dissimilarities (β-diversity) partially revealed distinct microbial compositions that differed between OC and cirrhosis patients (Fig. 5B). To search for differentially abundant taxa in the ascites microbiome of OC versus cirrhosis (Fig. 5C), we used the Deseq2 method to run an enrichment analysis. As ascites samples are considered low biomass samples and some amplicon sequence variants (ASVs) were found only in one sample, we ran the enrichment analysis once with ASVs expressed in more than one sample and once without a sample threshold. The enrichment analysis of the former revealed significant differences at the genus level, with one genus, Reyranella, being enriched in cirrhosis and one genus, Brevundimonas, being enriched in OC, with a log2 fold change of > 1.2 (Fig. 5D). Analysis without threshold revealed three and five significant ASVs in cirrhosis and OC, respectively (Fig. S5A).

Fig. 5
figure 5

Ascites IgA-seq can differentiate ovarian cancer patients from cirrhosis patients. A Alpha diversity analysis using the Chao1 model shows species richness in each group. Cyan samples are from cirrhosis patients, and pink samples, from ovarian cancer (OC) patients. B Beta diversity analysis shows the differences in species richness between groups. C Bar plot showing the abundance of each phylum in the ascites of each group. D Bar plot showing the log2 fold change (x axis), discovered using Deseq2, of genera expressed at least in two or more patients. Genera whose expression was higher in cirrhosis or OC are shown in cyan or pink, respectively. Genera written in black are significantly expressed in one of the groups, while those written in gray are not significantly expressed. E Bar plot showing the abundance of each phylum in the IgA-coated bacteria of the ascites of each group. F Bar plot showing the log2 fold change (x axis), discovered using Deseq2, of IgA-coated genera expressed at least in two or more patients. Genera whose expression was higher in cirrhosis or OC are shown in cyan or pink, respectively. Genera written in black are significantly expressed in one of the groups, while those written in gray are not significantly expressed

As we found an increase of IgA and infiltration of B cells in mice injected with TRAF3KO cells, we also isolated IgA-coated bacteria from the crude ascites from OC and cirrhosis patients, as described before [44], and amplified and sequenced the 16S rRNA (scheme of the process in Fig. S5B, created with Analysis of the ASVs showed a complete shift in the relative abundance of taxa at the phylum level in OC compared to cirrhosis (Fig. 5E). The enrichment analysis revealed significant differences at the genus level, with one genus, Acidovorax, enriched in cirrhosis and two genera, Enhydrobacter and Asinibacterium, enriched in OC (Fig. 5F). Without a threshold, we found six ASVs significantly enriched in cirrhosis and five ASVs significantly enriched in OC (Fig. S5C). These results indicate that ascites contain IgA-coated bacteria and the differences in IgA-coated bacteria between cirrhosis and OC suggest potential role(s) of B cells and IgA in controlling disease progression.


In this study, we identified TRAF3 as a potential immunomodulator in OC cells. The involvement of TRAF3 in cancer pathogenesis has been extensively studied (reviewed in [46]) and, while most evidence identifies TRAF3 as a tumor suppressor, as found in an HPV-positive head and neck cancer cohort [47], TRAF3 was recently found to have an immune-suppressive role and to limit the response to immune checkpoint blockade in melanoma [35]. In contrast to previous findings in solid cancers that TRAF3 deficiency leads mainly to the activation of both canonical and noncanonical NF-κB pathways [35, 47], we showed that knocking out TRAF3 in the ITB1 cell line did not affect NF-kB pathway activation, but rather hyperactivated the type-I IFN pathway. These contrasting findings emphasize the notion that TRAF3 functions differently in different cell types, particularly when it mediates the activity of type-I IFN [48]. Furthermore, the link between type-I IFN hyperactivation and TRAF3 can be dictated by other factors, such as the inactivation of BRCA1. DNA repair defects due to BRCA1/2 mutation induce immune signaling through the cGAS/STING/type-I IFN pathway [49]. TRAF3 inactivation may enhance the activation of STING and type-I IFN pathways in BRCA1-mutated models, but may not be sufficient for it by itself in BRCA1-wild-type models.

Type I IFNs have been shown to directly affect the proliferation, apoptosis, and other functions of tumor cells and also to indirectly affect these functions by increasing MHC-I expression [50]. Our study showed that TRAF3KO cells and ITB1 cells cultured with IFNs overexpressed MHC-I, consistent with the findings of other studies [35]. Higher levels of MHC-I enhance the expression of antigens, leading to increased recognition of antigens by antigen-presenting cells (APCs) and antigen presentation to cytotoxic T cells [51]. While we did not observe changes in the levels of cytotoxic CD8+ T cells, we did observe changes in the levels of macrophages and B cells in the mice bearing TRAF3KO tumors. These mice showed decreased levels of macrophages, which have previously been linked to the promotion of tumor metastasis, to enhanced angiogenesis and chemoresistance, and to poor prognosis in OC patients [52]. Likewise, the presence of plasma and memory B cells, which are often associated with improved prognosis in OC patients [11, 53, 54], were increased in TRAF3KO tumor-bearing mice compared to mice bearing ITB1 tumors.

The anti-tumor B cell response is largely achieved through the production of tumor-specific immunoglobulins that target tumor antigens and facilitate the recruitment and priming of APCs, leading to immune cytotoxic responses [11, 55]. In our model, we observed excess production and secretion of IgA in TRAF3KO tumor-bearing mice, further emphasizing the anti-tumor role of B cells. B cell activation and differentiation can be triggered through either T-dependent or T-independent mechanisms. In both cases, the initial signal that initiates the activation cascade relies on the recognition of antigens via the B cell receptor in T-dependent activation, or via Toll-like receptors in T-independent activation [56]. These antigens are typically derived from foreign proteins, polysaccharides, or metabolites from viral, fungal, or bacterial sources [57]. Therefore, in our study, we sought to investigate the correlation between the expression and function of IgA, bacterial genera, and survival in patients with OC. We found that the presence of certain bacterial genera correlated with the expression of IGHA1 in OC patients. We demonstrated experimentally in mice that antibiotics could prevent B cell-mediated growth delay of TRAF3KO OC cells.

The presence of microbiota in the TME has been shown to play a role in many malignancies, either directly or indirectly through secreted metabolites. While this link has been studied primarily in the context of colorectal cancer using stool samples (as reviewed in [58]), some studies have linked several bacterial genera to OC in humans [15, 16, 59] and mice [60]. Although our cohort of patients was small, we found three genera in significant abundance in OC ascites – Brevundimonas, Enhydrobacter, and Asinibacterium. These genera have previously been detected in ascites fluid and even in peritonitis [61,62,63], making their presence in the ascites of OC patients plausible.

We also characterized the level of IgA coating on these bacteria, since bacteria–IgA interactions are known to have a significant impact on the anti-tumor immune response in several ways. First, microbiota-derived MHC-I-restricted peptides can stimulate an immune response, which may cross-react with cancer antigens [64]. Second, microbiome-derived short-chain fatty acids can modulate the expression of key factors that influence B cell differentiation and antibody responses [65]. These specific immunoglobulin responses against commensal bacteria may indirectly affect tumor growth by altering the serum metabolome and cytokine repertoire [66, 67] or directly engage myeloid cells through FcR-dependent antibody-dependent cell-mediated cytotoxicity (ADCC) of tumor cells [11]. Therefore, to gain a better understanding of how the secreted IgAs and changes in the commensal microbiome may affect tumor cells directly or indirectly, it may be useful to examine the levels of different FcRs and cytokines.


Overall, our research has demonstrated that TRAF3 in OC cells is an immunosuppressive modulator that down-regulates MHC class I and IFN-I signaling, limits B-cell activation, and reduces the anti-tumor immunity of B cells. We further show that the presence of an intact commensal microbiome is necessary for an effective anti-tumor B cell response. (Scheme in Fig. 6A, created with By sorting and 16S sequencing IgA coated bacteria of the ascites microbiota of cirrhosis and OC patients, we showed that the bacteria coating pattern is highly disease-associated (Scheme in Fig. 6B, created with

Fig. 6
figure 6

Schematic illustrations demonstrating the key steps of our study. A shows (i) the injection of ITB1 cells into mice, resulting in rapid tumor progression and mice dying within ~ three months. (ii) TRAF3KO cells are more immunogenic with higher expression of MHC-I (depicted in purple) and type-I IFN related genes (represented by blue circles) compared to ITB1. Injection of TRAF3KO cells into wild-type mice leads to the infiltration of IgA-secreting B cell populations, which are associated with prolonged mouse survival compared to mice injected with ITB1. (iii) The presence of microbiota is required for the in-vivo immunogenicity of TRAF3KO, as elimination of the microbiota using antibiotics enhances tumor progression. B displays 16S sequencing of the microbiota and IgA-bound bacteria in the ascites of OC and cirrhosis patients


Despite the promising results obtained in our study, there are several limitations that we should be aware of. Firstly, our study utilized a genetically modified version of the ID8 cell line, which is the only widely-used syngeneic model of OC. Hence, our findings will need to be confirmed in future studies using other tumorigenic OC cell lines with BRCA1 mutations. Secondly, while the analysis of the TME showed an increase in B cells and no changes in effector cells like CD8 T or NK cells, further analysis is required to identify the immune subset of cells that directly eliminate TRAF3KO cells. Thirdly, the analysis of a small cohort consisting of only 11 ascites samples from OC patients and 4 from cirrhosis patients provided only limited insights into the ascites bacterial abundance of these diseases. Moreover, this small cohort did not produce any clinically relevant results, and we could not associate genera with enhanced anti-tumor immunity. Finally, while 16 s sequencing provided general information about the bacterial sample abundance and allowed us to show the overall trend and to emphasize the advantages of using IgA-seq to better characterize OC patients, deep sequencing may be necessary to provide more comprehensive information regarding the function of these genera.

Availability of data and materials

Datasets generated and analyzed during the current study are available from the corresponding author upon reasonable request.


APCs :

Antigen-presenting cells


Amplicon sequence variant



ITB1 :


OC :

Ovarian cancer


Tumor microenvironment


Tumor necrosis factor receptor-associated factor 3




Wild type


  1. Global Cancer Observatory, (n.d.). (Accessed 23 Nov 2021).

  2. Konstantinopoulos PA, Ceccaldi R, Shapiro GI, D’Andrea AD. Homologous Recombination Deficiency: Exploiting the Fundamental Vulnerability of Ovarian Cancer. Cancer Discov. 2015;5:1137–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Ferlay J, Soerjomataram I, Dikshit R, Eser S, Mathers C, Rebelo M, Parkin DM, Forman D, Bray F. Cancer incidence and mortality worldwide: Sources, methods and major patterns in GLOBOCAN 2012. Int J Cancer. 2015.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Borella F, Ghisoni E, Giannone G, Cosma S, Benedetto C, Valabrega G, Katsaros D. Immune Checkpoint Inhibitors in Epithelial Ovarian Cancer: An Overview on Efficacy and Future Perspectives. Diagnostics. 2020;10:146.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Postow MA, Callahan MK, Wolchok JD. Immune Checkpoint Blockade in Cancer Therapy. J Clin Oncol. 2015;33:1974.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. González-Martín A, Sánchez-Lorenzo L. Immunotherapy with checkpoint inhibitors in patients with ovarian cancer: Still promising? Cancer. 2019;125:4616–22.

    Article  PubMed  Google Scholar 

  7. Ning F, Cole CB, Annunziata CM. Driving Immune Responses in the Ovarian Tumor Microenvironment. Front Oncol. 2021;10:3060.

    Article  Google Scholar 

  8. Wei X, Jin Y, Tian Y, Zhang H, Wu J, Lu W, Lu X. Regulatory B cells contribute to the impaired antitumor immunity in ovarian cancer patients. Tumor Biol. 2016;37:6581–8.

    Article  CAS  Google Scholar 

  9. Kroeger DR, Milne K, Nelson BH. Tumor-Infiltrating Plasma Cells Are Associated with Tertiary Lymphoid Structures, Cytolytic T-Cell Responses, and Superior Prognosis in Ovarian Cancer. Clin Cancer Res. 2016;22:3005–15.

    Article  CAS  PubMed  Google Scholar 

  10. Montfort A, Pearce O, Maniati E, Vincent BG, Bixby L, Böhm S, Dowe T, Wilkes EH, Chakravarty P, Thompson R, Topping J, Cutillas PR, Lockley M, Serody JS, Capasso M, Balkwill FR. A strong B-cell response is part of the immune landscape in human high-grade serous ovarian metastases. Clin Cancer Res. 2017.

    Article  PubMed  Google Scholar 

  11. Biswas S, Mandal G, Payne KK, Anadon CM, Gatenbee CD, Chaurio RA, Costich TL, Moran C, Harro CM, Rigolizzo KE, Mine JA, Trillo-Tinoco J, Sasamoto N, Terry KL, Marchion D, Buras A, Wenham RM, Yu X, Townsend MK, Tworoger SS, Rodriguez PC, Anderson AR, Conejo-Garcia JR. IgA transcytosis and antigen recognition govern ovarian cancer immunity. Nat. 2021;591(7850):464–70.

    Article  CAS  Google Scholar 

  12. Cerutti A. The regulation of IgA class switching. Nat Rev Immunol. 2008;8(6):421–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Slack E, Balmer ML, Macpherson AJ. B cells as a critical node in the microbiota–host immune system network. Immunol Rev. 2014;260:50–66.

    Article  CAS  PubMed  Google Scholar 

  14. Goedert JJ, Hua X, Bielecka A, Okayasu I, Milne GL, Jones GS, Fujiwara M, Sinha R, Wan Y, Xu X, Ravel J, Shi J, Palm NW, Feigelson HS. Postmenopausal breast cancer and oestrogen associations with the IgA-coated and IgA-noncoated faecal microbiota. Br J Cancer. 2018;118(4):471–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Zhou B, Sun C, Huang J, Xia M, Guo E, Li N, Lu H, Shan W, Wu Y, Li Y, Xu X, Weng D, Meng L, Hu J, Gao Q, Ma D, Chen G. The biodiversity Composition of Microbiome in Ovarian Carcinoma Patients. Sci Reports. 2019;9(1):1–11.

    Article  CAS  Google Scholar 

  16. Nené NR, Reisel D, Leimbach A, Franchi D, Jones A, Evans I, Knapp S, Ryan A, Ghazali S, Timms JF, Paprotka T, Bjørge L, Zikan M, Cibula D, Colombo N, Widschwendter M. Association between the cervicovaginal microbiome, BRCA1 mutation status, and risk of ovarian cancer: a case-control study. Lancet Oncol. 2019;20:1171–82.

    Article  PubMed  Google Scholar 

  17. Walton JB, Farquharson M, Mason S, Port J, Kruspig B, Dowson S, Stevenson D, Murphy D, Matzuk M, Kim J, Coffelt S, Blyth K, McNeish IA. CRISPR/Cas9-derived models of ovarian high grade serous carcinoma targeting Brca1, Pten and Nf1, and correlation with platinum sensitivity. Sci Rep. 2017.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Shalem O, Sanjana NE, Hartenian E, Shi X, Scott DA, Mikkelsen TS, Heckl D, Ebert BL, Root DE, Doench JG, Zhang F. Genome-Scale CRISPR-Cas9 Knockout Screening in Human Cells. Science. 2014;343:84.

    Article  CAS  PubMed  Google Scholar 

  19. Nuclear & Cytoplasmic Extract Protocol | Rockland, (n.d.). (Accessed 6 Sept 2022).

  20. Ochoa-Repáraz J, Mielcarz DW, Ditrio LE, Burroughs AR, Foureau DM, Haque-Begum S, Kasper LH. Role of Gut Commensal Microflora in the Development of Experimental Autoimmune Encephalomyelitis. J Immunol. 2009;183:6041–50.

    Article  CAS  PubMed  Google Scholar 

  21. Anahtar MN, Bowman BA, Kwon DS. Efficient nucleic acid extraction and 16s rRNA gene sequencing for bacterial community characterization. J Vis Exp. 2016.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Winter J, Breinig M, Heigwer F, Brügemann D, Leible S, Pelz O, Zhan T, Boutros M. caRpools: an R package for exploratory data analysis and documentation of pooled CRISPR/Cas9 screens. Bioinformatics. 2016;32:632–4.

    Article  CAS  PubMed  Google Scholar 

  23. Yu G, He QY. ReactomePA: An R/Bioconductor package for reactome pathway analysis and visualization. Mol Biosyst. 2016;12:477–9.

    Article  CAS  PubMed  Google Scholar 

  24. Holm JB, Humphrys MS, Robinson CK, Settles ML, Ott S, Fu L, et al. Ultrahigh-Throughput Multiplexing and Sequencing of >500-Base-Pair Amplicon Regions on the Illumina HiSeq 2500 Platform, MSystems. 2019;4.

  25. Babraham Bioinformatics - FastQC A Quality Control tool for High Throughput Sequence Data, (n.d.). (Accessed 5 Jan 2023).

  26. Ewels P, Magnusson M, Lundin S, Käller M. MultiQC: summarize analysis results for multiple tools and samples in a single report. Bioinformatics. 2016;32:3047–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Chen S, Zhou Y, Chen Y, Gu J. fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics. 2018;34:i884–90.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Bolyen E, Rideout JR, Dillon MR, Bokulich NA, Abnet CC, Al-Ghalith GA, et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat Biotechnol. 2019;37(8):852–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, Peplies J, Glöckner FO. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013;41:D590–6.

    Article  CAS  PubMed  Google Scholar 

  30. Davis NM, Proctor DM, Holmes SP, Relman DA, Callahan BJ. Simple statistical identification and removal of contaminant sequences in marker-gene and metagenomics data. Microbiome. 2018;6:1–14.

    Article  Google Scholar 

  31. Liu Y, Elworth RAL, Jochum MD, Aagaard KM, Treangen TJ. De novo identification of microbial contaminants in low microbial biomass microbiomes with Squeegee, Nat. Commun. 2022 131. 13 (2022) 1–14. doi:

  32. Liu C, Cui Y, Li X, Yao M. Microeco: An R package for data mining in microbial community ecology. FEMS Microbiol Ecol. 2021;97.

  33. McMurdie PJ, Holmes S. Phyloseq: An R Package for Reproducible Interactive Analysis and Graphics of Microbiome Census Data. PLoS ONE. 2013.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Pan D, Kobayashi A, Jiang P, De Andrade LF, Tay RE, Luoma AM, Tsoucas D, Qiu X, Lim K, Rao P, Long HW, Yuan GC, Doench J, Brown M, Liu XS, Wucherpfennig KW. A major chromatin regulator determines resistance of tumor cells to T cell-mediated killing. Science (80-). 2018;359:770–5.

    Article  CAS  Google Scholar 

  35. Gu SS, Zhang W, Wang X, Jiang P, Traugh N, Li Z, et al. Therapeutically Increasing MHC-I Expression Potentiates Immune Checkpoint Blockade. Cancer Discov. 2021.

  36. Häcker H, Tseng PH, Karin M. Expanding TRAF function: TRAF3 as a tri-faced immune regulator. Nat Rev Immunol. 2011.

    Article  PubMed  Google Scholar 

  37. Bindea G, Mlecnik B, Tosolini M, Kirilovsky A, Waldner M, Obenauf AC, Angell H, Fredriksen T, Lafontaine L, Berger A, Bruneval P, Fridman WH, Becker C, Pagès F, Speicher MR, Trajanoski Z, Galon JÔ. Spatiotemporal Dynamics of Intratumoral Immune Cells Reveal the Immune Landscape in Human Cancer. Immunity. 2013;39:782–95.

    Article  CAS  PubMed  Google Scholar 

  38. [The primary structure of a monoclonal IgA-immunoglobulin (IgA Tro.), II. The amino acid sequence of the H-chain, alpha-type, subgroup III; structure of the complete IgA-molecule (author’s transl)] - PubMed, (n.d.). (Accessed 26 Jan 2023).

  39. Zheng D, Liwinski T, Elinav E. Interaction between microbiota and immunity in health and disease. Cell Res. 2020;30(6):492–506.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Botía-Sánchez M, Alarcón-Riquelme ME, Galicia G. B Cells and Microbiota in Autoimmunity. Int J Mol Sci. 2021;22:4846.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Noble A, Pring ET, Durant L, Man R, Dilke SM, Hoyles L, James SA, Carding SR, Jenkins JT, Knight SC. Altered immunity to microbiota, B cell activation and depleted γδ/resident memory T cells in colorectal cancer. Cancer Immunol Immunother. 2022;71:2619–29.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Savage DC, Dubos R. Alterations in the mouse cecum and its flora produced by antibacterial drugs. J Exp Med. 1968;128:97–110.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Gabay O, Vicenty J, Smith D, Tiffany L, Ascher J, Curry T, Dennis J, Clouse KA. Using a Model of Germ-Free Animals to Study the Impact of Gut Microbiome in Research: A Step by Step Sterility Setting and Management. Methods Protoc. 2020;3:18.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Palm NW, De Zoete MR, Cullen TW, Barry NA, Stefanowski J, Hao L, Degnan PH, Hu J, Peter I, Zhang W, Ruggiero E, Cho JH, Goodman AL, Flavell RA. Immunoglobulin A Coating Identifies Colitogenic Bacteria in Inflammatory Bowel Disease. Cell. 2014;158:1000–10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Poore GD, Kopylova E, Zhu Q, Carpenter C, Fraraccio S, Wandro S, Kosciolek T, Janssen S, Metcalf J, Song SJ, Kanbar J, Miller-Montgomery S, Heaton R, Mckay R, Patel SP, Swafford AD, Knight R. Microbiome analyses of blood and tissues suggest cancer diagnostic approach. Nat. 2020;579(7800):567–74.

    Article  CAS  Google Scholar 

  46. Zhu S, Jin J, Gokhale S, Lu AM, Shan H, Feng J, Xie P. Genetic alterations of TRAF proteins in human cancers. Front Immunol. 2018;9:2111.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Zhang J, Chen T, Yang X, Cheng H, Spath SS, Clavijo PE, Chen J, Silvin C, Issaeva N, Su X, Yarbrough WG, Annunziata CM, Chen Z, Van Waes C. Attenuated TRAF3 fosters activation of alternative NF-kB and reduced expression of antiviral interferon, TP53, and RB to promote HPV-positive head and neck cancers. Cancer Res. 2018;78:4613–26.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Xie X, Jin J, Zhu L, Jie Z, Li Y, Zhao B, Cheng X, Li P, Sun SC. Cell type-specific function of TRAF2 and TRAF3 in regulating type i IFN induction 11 Medical and Health Sciences 1107 Immunology. Cell Biosci. 2019;9:1–10.

    Article  Google Scholar 

  49. van Vugt MATM, Parkes EE. When breaks get hot: inflammatory signaling in BRCA1/2-mutant cancers. Trends in Cancer. 2022;8:174–89.

    Article  PubMed  Google Scholar 

  50. Yu R, Zhu B, Chen D. Type I interferon-mediated tumor immunity and its role in immunotherapy. Cell Mol Life Sci. 2022;79:1–24.

    Article  Google Scholar 

  51. Parker BS, Rautela J, Hertzog PJ. Antitumour actions of interferons: implications for cancer therapy. Nat Rev Cancer. 2016;16(3):131–44.

    Article  CAS  PubMed  Google Scholar 

  52. Schweer D, McAtee A, Neupane K, Richards C, Ueland F, Kolesar J. Tumor-Associated Macrophages and Ovarian Cancer: Implications for Therapy. Cancers. 2022;14:2220.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Santoiemma PP, Reyes C, Wang LP, McLane MW, Feldman MD, Tanyi JL, Powell DJ. Systematic evaluation of multiple immune markers reveals prognostic factors in ovarian cancer. Gynecol Oncol. 2016;143:120–7.

    Article  CAS  PubMed  Google Scholar 

  54. Montfort A, Pearce O, Maniati E, Vincent BG, Bixby L, Böhm S, Dowe T, Wilkes EH, Chakravarty P, Thompson R, Topping J, Cutillas PR, Lockley M, Serody JS, Capasso M, Balkwill FR. A strong B-cell response is part of the immune landscape in human high-grade serous ovarian metastases. Clin Cancer Res. 2017;23:250–62.

    Article  CAS  PubMed  Google Scholar 

  55. Montfort A, Pearce O, Maniati E, Vincent BG, Bixby L, B€ Ohm S, Dowe T, Wilkes E.H, Chakravarty P, Thompson R, Topping J, Cutillas PR, Lockley M, Serody JS, Capasso M, Balkwill FR. Biology of Human Tumors A Strong B-cell Response Is Part of the Immune Landscape in Human High-Grade Serous Ovarian Metastases. Clin Cancer Res. 2017;23.

  56. Yu B, Wang L, Chu Y. Gut microbiota shape B cell in health and disease settings. J Leukoc Biol. 2021;110:271–81.

    Article  CAS  PubMed  Google Scholar 

  57. Heesters BA, Van Der Poel CE, Das A, Carroll MC. Antigen Presentation to B Cells. Trends Immunol. 2016;37:844–54.

    Article  CAS  PubMed  Google Scholar 

  58. Wong SH, Yu J. Gut microbiota in colorectal cancer: mechanisms of action and clinical applications. Nat Rev Gastroenterol Hepatol. 2019;16(11):690–704.

    Article  CAS  PubMed  Google Scholar 

  59. Asangba AE, Chen J, Goergen KM, Larson MC, Oberg AL, Casarin J, Multinu F, Kaufmann SH, Mariani A, Chia N, Walther-Antonio MRS. Diagnostic and prognostic potential of the microbiome in ovarian cancer treatment response. Sci Reports. 2023;13(1):1–18.

    Article  CAS  Google Scholar 

  60. Chambers LM, Esakov Rhoades EL, Bharti R, Braley C, Tewari S, Trestan L, et al. Disruption of the Gut Microbiota Confers Cisplatin Resistance in Epithelial Ovarian Cancer. Cancer Res. 2022;82:4654–69.

    Article  CAS  PubMed  Google Scholar 

  61. Enomoto H, Inoue SI, Matsuhisa A, Nishiguchi S. Diagnosis of Spontaneous Bacterial Peritonitis and an In Situ Hybridization Approach to Detect an “Unidentified” Pathogen. Int J Hepatol. 2014;(2014).

  62. Lee SR, Lee JC, Kim SH, Oh YS, Chae HD, Seo H, Kang CS, Shin TS. Altered Composition of Microbiota in Women with Ovarian Endometrioma: Microbiome Analyses of Extracellular Vesicles in the Peritoneal Fluid. Int J Mol Sci. 2021;22.

  63. Menard J, Goggs R, Mitchell P, Yang Y, Robbins S, Franklin-Guild RJ, Thachil AJ, Altier C, Anderson R, Putzel GG, McQueary H, Goodman LB. Effect of antimicrobial administration on fecal microbiota of critically ill dogs: dynamics of antimicrobial resistance over time. Anim Microbiome. 2022;4.

  64. Fluckiger A, Daillère R, Sassi M, Sixt BS, Liu P, Loos F, Richard C, Rabu C, Alou MT, Goubet AG, Lemaitre F, Ferrere G, Derosa L, Duong CPM, Messaoudene M, Gagné A, Joubert P, de Sordi L, Debarbieux L, Simon S, Scarlata CM, Ayyoub M, Palermo B, Facciolo F, Boidot R, Wheeler R, Boneca IG, Sztupinszki Z, Papp K, Csabai I, Pasolli E, Segata N, Lopez-Otin C, Szallasi Z, Andre F, Iebba V, Quiniou V, Klatzmann D, Boukhalil J, Khelaifia S, Raoult D, Albiges L, Escudier B, Eggermont A, Mami-Chouaib F, Nistico P, Ghiringhelli F, Routy B, Labarrière N, Cattoir V, Kroemer G, Zitvogel L. Cross-reactivity between tumor MHC class I–restricted antigens and an enterococcal bacteriophage. Science. 2020;369:936–42.

    Article  CAS  PubMed  Google Scholar 

  65. White CA, Pone EJ, Lam T, Tat C, Hayama KL, Li G, Zan H, Casali P. Histone Deacetylase Inhibitors Upregulate B Cell microRNAs That Silence AID and Blimp-1 Expression for Epigenetic Modulation of Antibody and Autoantibody Responses. J Immunol. 2014;193:5933–50.

    Article  CAS  PubMed  Google Scholar 

  66. Iida N, Dzutsev A, Stewart CA, Smith L, Bouladoux N, Weingarten RA, Molina DA, Salcedo R, Back T, Cramer S, Dai RM, Kiu H, Cardone M, Naik S, Patri AK, Wang E, Marincola FM, Frank KM, Belkaid Y, Trinchieri G, Goldszmid RS. Commensal bacteria control cancer response to therapy by modulating the tumor microenvironment. Science (80-). 2013;342:967–70.

    Article  CAS  PubMed Central  Google Scholar 

  67. Davar D, Dzutsev AK, McCulloch JA, Rodrigues RR, Chauvin JM, Morrison RM, Deblasio RN, Menna C, Ding Q, Pagliano O, Zidi B, Zhang S, Badger JH, Vetizou M, Cole AM, Fernandes MR, Prescott S, Costa RGF, Balaji AK, Morgun A, Vujkovic-Cvijin I, Wang H, Borhani AA, Schwartz MB, Dubner HM, Ernst SJ, Rose A, Najjar YG, Belkaid Y, Kirkwood JM, Trinchieri G, Zarour HM. Fecal microbiota transplant overcomes resistance to anti-PD-1 therapy in melanoma patients. Science (80-). 2021;371:595–602.

    Article  CAS  Google Scholar 

Download references


We would like to thank Prof. Iain A. McNeish and his members of his lab for supplying us with the ITB1 cells. We appreciate the effort of the physicians in both hospitals, Ziv and Rabin, for collecting the ascites samples of the patients.


This work was funded by the Israel Science Foundation (ISF, 302/21 and 700/16) (to ME), the Israel Cancer Research Foundation (ICRF, 17–1693-RCDA) (to ME), United States-Israel Binational Science Foundation (BSF, 2021055) (to ME), German Cancer Research Center- Ministry of Science and Technology (DKFZ-MOST #001192) and Israel Cancer Association (ICA, 20220009) (to JMG).

Author information

Authors and Affiliations



JZ, YM, RDM, YKC, ZS, JM, JMG and ME were involved in experiment design and execution. JZ, ME, YM, and JGM wrote the manuscript. All authors reviewed the final version. The author(s) read and approved the final manuscript.

Corresponding author

Correspondence to Moshe Elkabets.

Ethics declarations

Ethics approval and consent to participate

Animal experiments were conducted in compliance with protocols established by the Institutional Animal Care and Use Committee (IACUC) of Ben-Gurion University of the Negev for ensuring animal welfare and minimizing discomfort. The animal ethical clearance protocol number used for the study was IL-23–05-2020(E).

The acquisition and analysis of human ascites samples was approved by the ethics in research committee of the Rabin Medical Center and Ziv Medical Center under protocol 0450–16-RMC and 201913481, respectively. Informed consent was obtained from all participants prior to their recruitment to the study.

Consent for publication

Not applicable.

Competing interests

All authors have no potential conflicts of interest to declare.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Fig. S1.

Establishing a library-expressing cells. A WT mouse harboring ITB1-related ascites and small tumors, indicated by the white arrows or the black arrows in the IHC-stained tissues. MLN- mesenchymal lymph node. B Cas9-infected ITB monoclonal isolated colonies express an active (>95%) Cas9. In all flow cytometry experiments, we used a PD-1 guide and determined its expression before and after stimulation with IFN-γ. C Bar graph showing the number of reads per guide. Guides with the same number of reads are grouped together and indicated by the same color. Fig. S2. Knocking out TRAF3 does not affect the NF-κB pathway. A Scratch assay showing the migration of all cells over 48 h into the scratched area. Data was acquired using a live cell imaging system (Julistage). B Western blot of ITB1 and TRAF3KO lysates comparing the levels of proteins of the NF-κB pathway. The levels of p65 and RelB were monitored in the nucleus, and p100 and p52 were monitored in the whole lysates. Actin and histone H3 were used as loading controls. C Quantitative PCR shows the mRNA levels of genes related to the NF-κB pathway (RelB, IL-6, and IκBα). Cyan bars represent ITB1, solid pink bars represent TRAF3KO1, and dotted pink bars represent TRAF3KO2. D MFI of the blot on the main figure, calculated using ImageJ E Bar graph showing the differences in mean fluorescence intensity of the Cy3 channel in the nucleus between ITB1 (cyan) and TRAF3KO cells (pink). Fig. S3. B cells in TRAF3KO-injected mice express and secrete more IgAs than in ITB1-injected mice. A Expression of the costimulatory protein CD40 and interleukin-10 (IL-10) on B cells in ITB1-injected and TRAF3KO-injected mice. B, C Ig isotyping array (B) and its quantification (C). D Scatter plot showing the differentiation between the two clusters of the heatmap. E Heatmap showing the average score of each cluster for each immune pathway. F Survival curve showing the difference between patients with high B cell score (cyan) and low B cell score (dark pink). Fig. S4. Effect of ABX treatment. A Mouse treated with ABX for 10 days and its enlarged cecum. B Fecal samples were collected from both ABX-treated and naïve mice and seeded on a blood agar plate. Growth of colonies can be seen only on the left side of the plate where the naïve cells was seeded. C Difference in the abundance of three genera, compared using the Wilcoxon signed-rank test. Fig. S5. Unique bacteria in the ascites of ovarian cancer patients. A Bar plot showing the log2 fold change (x axis), discovered using Deseq2, of genera expressed at least in one or more patients. Genera expressed more in cirrhosis or OC are in cyan or pink, respectively. Genus with black name are expressed significantly in one of the groups, while the grey ones are not significant. B a schematic illustration describing the separation of IgA-coated bacteria from the ascites fluid. C Bar plot showing the log2 fold change (x axis), discovered using Deseq2, of IgA-coated genera expressed at least in one or more patients. Genera whose expression was higher in cirrhosis or OC are shown in cyan or pink, respectively. Genera written in black are significantly expressed in one of the groups, while those written in gray are not significantly expressed.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zorea, J., Motro, Y., Mazor, R.D. et al. TRAF3 suppression encourages B cell recruitment and prolongs survival of microbiome-intact mice with ovarian cancer. J Exp Clin Cancer Res 42, 107 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: