Front. Immunol. Frontiers in Immunology Front. Immunol. 1664-3224 Frontiers Media S.A. 10.3389/fimmu.2023.1241632 Immunology Original Research Refined control of CRISPR-Cas9 gene editing in Clostridium sporogenes: the creation of recombinant strains for therapeutic applications Kubiak Aleksandra M. 1 2 † Claessen Luuk 1 2 Zhang Yanchao 2 Khazaie Khashayarsha 3 * Bailey Tom S. 2 * † 1 Exomnis Biotech BV, Maastricht, Netherlands 2 The M-Lab, Department of Precision Medicine, GROW - School for Oncology and Reproduction, Maastricht University, Maastricht, Netherlands 3 Department of Immunology and Cancer Biology, Mayo Clinic, Phoenix, AZ, United States

Edited by: Venkatakrishna Rao Jala, University of Louisville, United States

Reviewed by: Shigeru Kamiya, Kyorin University, Japan; Prabhat Kumar Talukdar, Washington State University, United States

*Correspondence: Khashayarsha Khazaie, Khazaie@mayo.edu; Tom S. Bailey, t.bailey@maastrichtuniversity.nl

†These authors share senior authorship

 05 10 2023 2023 14 1241632 16 06 2023 01 09 2023 Copyright © 2023 Kubiak, Claessen, Zhang, Khazaie and Bailey 2023 Kubiak, Claessen, Zhang, Khazaie and Bailey

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.

Despite considerable clinical success, the potential of cancer immunotherapy is restricted by a lack of tumour-targeting strategies. Treatment requires systemic delivery of cytokines or antibodies at high levels to achieve clinically effective doses at malignant sites. This is exacerbated by poor penetration of tumour tissue by therapeutic antibodies. High-grade immune-related adverse events (irAEs) occur in a significant number of patients (5-15%, cancer- and therapeutic-dependent) that can lead to lifelong issues and can exclude from treatment patients with pre-existing autoimmune diseases. Tumour-homing bacteria, genetically engineered to produce therapeutics, is one of the approaches that seeks to mitigate these drawbacks. The ability of Clostridium sporogenes to form spores that are unable to germinate in the presence of oxygen (typical of healthy tissue) offers a unique advantage over other vectors. However, the limited utility of existing gene editing tools hinders the development of therapeutic strains. To overcome the limitations of previous systems, expression of the Cas9 protein and the gRNA was controlled using tetracycline inducible promoters. Furthermore, the components of the system were divided across two plasmids, improving the efficiency of cloning and conjugation. Genome integrated therapeutic genes were assayed biochemically and in cell-based functional assays. The potency of these strains was further improved through rationally-conceived gene knock-outs. The new system was validated by demonstrating the efficient addition and deletion of large sequences from the genome. This included the creation of recombinant strains expressing two pro-inflammatory cytokines, interleukin-2 (IL-2) and granulocyte macrophage-colony stimulating factor (GM-CSF), and a pro-drug converting enzyme (PCE). A comparative, temporal in vitro analysis of the integrant strains and their plasmid-based equivalents revealed a substantial reduction of cytokine activity in chromosome-based constructs. To compensate for this loss, a 7.6 kb operon of proteolytic genes was deleted from the genome. The resultant knock-out strains showed an 8- to 10-fold increase in cytokine activity compared to parental strains.

 CRISPR-Cas9 Clostridium sporogenes secretion cancer cytokines immunotherapy KWF Kankerbestrijding10.13039/501100004622 section-in-acceptance Cancer Immunity and Immunotherapy Introduction

Immunotherapy has revolutionised modern clinical oncology. The concept of exploiting natural anti-microbial processes to steer the immune response against malignant tissue was exploited over a century ago by William Coley and led to attempts to enhance natural defences against cancer with immunogenised tumour cell vaccines (1). Earlier attempts with cytokine expression (2) and viral modification of tumour cells (3) produced promising results in animal models and led to clinical trials that are still ongoing (4, 5). A deeper understanding of novel immune system approaches for removing breaks through an immune checkpoint blockade rather than by stimulating immunity has introduced a new era of hope and expectation in cancer therapy (6–9). Nevertheless, major obstacles, including immune adverse effects and accelerated disease in response to therapy (10), remain to be understood and overcome. The search for complimentary approaches has revitalized attempts for targeting the tumour microenvironment through pharmacologic or immunologic strategies (3, 4, 11, 12). Tumour lytic viruses (13, 14) [reviewed in detail by Syyam et al. (15)] and bacteria (16) [reviewed in detail by Hu et al. (12)] have been at the forefront of the latter strategy in clinical attempts to hinder advanced cancers.

The observation that bacteria can colonise and persist in cancerous tumours is a long-standing one. Reports as early as 1813 observed regression of pre-existing tumours following gas gangrene infections, but more recent reports document the existence of complex communities of bacteria in tumours. This brings attention to several opportunities and raises further questions. In most cases, tumour-resident bacteria are opportunistic inhabitants rather than the initial cause of the cancer (17). The source of this “tumour microbiome” is not clear, but the gastrointestinal (GI) tract is an obvious candidate. There is evidence, both experimental and natural, of translocation from the GI tract (18, 19). Using modern RNA and DNA sequencing techniques, the profile of this community can be determined (20). This profile could be exploited to predict cancer prognosis and treatment response (21) or to improve treatment strategies (22). While a range of bacterial species, both obligately and facultatively anaerobic, have been reported to occur in tumours, the presence of Clostridium species is a recurring theme, which neatly links back to similar reports 200 years ago (23). A major advance in the use of bacteria for targeting of tumours was the creation of an attenuated strain of Clostridium novyi, an obligately anaerobic clostridia (24). Clostridium novyi-NT entered clinical trials (16) after extensive animal modelling in mice, rats, and dogs, showing promising results particularly in large animal clinical trials (25). Other examples include Salmonella vaccines (26, 27), which have been attenuated through genetic modifications (28) or self-destructive derivatives to overcome potential safety concerns (29). More recent attempts involve the use of viruses or bacteria for the intratumoural delivery of cargo, including tumour antigens (30).

Whether viral or bacterial, the challenges to using these vectors for intratumoural delivery of therapeutics are common. Appropriate selection of therapeutic agents aside, clinical efficacy will depend on the potency and dose of the agent at the site of disease. Since its inception, numerous innovations have yielded vast gains in the specificity and potency of immunotherapeutic agents, at least preclinically. Amino acid substitutions and fusion to receptor or antibody domains have enabled improvements in both the target affinity and half-life of the agents (31). Whether or not to increase the dose of the therapeutic is a decision that lies in the hands of the synthetic biologist and is something that draws on vast collective multidisciplinary knowledge. In the case of bacterial vectors, a rational approach is to optimise the steps between the DNA of the recombinant gene and the malignant target: transcription, translation, and secretion. Screening of gene promoters may be the simplest approach to increase transcription, and this has yielded impressive results in Bacteroides species (32). Increasing mRNA stability could also enhance expression, but it remains poorly characterised (33). Secretion, at least via the Sec pathway, has been studied extensively in several bacteria, including the model organism Bacillus subtilis (34–39). Equivalent systematic studies are yet to be conducted for Clostridium species.

The ability to precisely edit the genome brings significant benefits to therapeutic strain development. Through deletion of specific genes or modification of their expression, the phenotype of the host can be adjusted to suit any given application. Genome editing also enables containment of genetic elements, while plasmid-based expression cassettes introduce the risk of horizontal gene transfer, both in the environment (40) and within mammalian hosts (41, 42). Integration of therapeutic genes averts these risks and reduces the rate of mutation (43–45). In the last 10 years, numerous methods for gene editing in Clostridium species have been reported, the majority of which are based on CRISPR-Cas systems (46–50). While these illustrate that considerable progress has been made in this area, existing systems are hindered by low plasmid transfer efficiency and a lack of potency.

This paper reports the development of C. sporogenes, an obligate anaerobic, spore-forming bacteria, as a therapeutic delivery vector. Building on previous work (51–53), the aim of this study was to overcome the obstacles that have prevented this approach from becoming a viable clinical product. Specifically, our aim was to develop a synthetic strain that can produce functional product at therapeutic levels from a chromosome-based heterologous expression cassette. In pursuit of these objectives, a highly efficient and quick two-plasmid CRISPR-Cas9 system that enables deletion and insertion of large sequences (≥7 kb) was developed. Using the improved method, recombinant strains producing either cytokines or prodrug-converting enzyme (PCE) were created and validated. In addition, modification of the host’s phenotype through the targeted deletion of proteolytic genes resulted in a 7-fold increase of heterologous product.

 Materials and methods Bacterial strains, growth conditions and cell lines

Details of all bacterial strains used in this study are listed in Table S1 . C. sporogenes-NT, a non-haemolytic strain of C. sporogenes NCIMB 10696 (henceforth referred to as WT), was created previously by deleting the putative Streptolysin S (SLS) operon (51). Two Escherichia coli strains (10-beta and S17-1) were used as the expression and conjugative donor strains, respectively.

Growth of C. sporogenes strains was carried out under anaerobic conditions in an anaerobic cabinet (model MG1000 Mark II, Don Whitley Scientific Ltd.; 80% N2, 10% CO2, 10% H2) at 37°C. C. sporogenes was grown in a bovine-free version of TY media, labelled “Peptone Yeast Thioglycolate” (PYT), supplemented with D-cycloserine (Cs, 250 μg/ml), thiamphenicol (Tm, 15 μg/ml), erythromycin (Erm, 30 μg/ml), and anhydrotetracycline (aTc, 96 ng/ml) where appropriate.

E. coli strains were grown aerobically at 37°C and 200 rpm in liquid LB medium or solid LB with 1.5% agar supplemented with Erm (300 μg/ml) or chloramphenicol (Cm, 25 μg/ml for solid media and 12.5 μg/ml for liquid media) or aTc (32 or 96 ng/mL).

Murine bone marrow FDCP-1 (“Factor-Dependent Continuous cell line, Paterson Laboratories”, ACC 368, DSMZ) and cytotoxic T-cells CTLL-2 cell lines (93042610, ECACC) were cultivated according to the manufacturer’s instructions.

 Generation of plasmid-based cargo strains in <italic>C. sporogenes</italic>-NT

Details of all vectors and specific oligonucleotides used in this study are listed in Tables S1 and S2 of the Supplementary Material. For amplification of genomic or plasmid DNA, a high-fidelity polymerase was used according to the manufacturer’s protocol (Phusion Master Mix with HF buffer, M0531, New England Biolabs, NEB).

To construct plasmids harbouring three different cargos, FMN reductase - NfrA, murine interleukin 2 (mIL-2) and murine granulocyte-macrophage colony-stimulating factor (mGM-CSF), the DNA sequences were codon optimised according to the known codon usage preference of C. sporogenes and based on the reported sequence (UniProtKB: P39605, P04351, P01587). DNA fragments of desired cargo (NfrA, mIL-2 or mGM-CSF) and promoter/promoter-signal sequences were commercially synthesised (Integrated DNA Technologies, IDT Inc.) with sites recognized by type IIS restriction enzymes ( Table S3 ).

The pGG212L vector was used as described before (53) to generate three constructs via Golden Gate Assembly (GGA), using BsaI sites. Ligated plasmids were transformed by heat shock into E. coli 10-beta and plated on LB plates supplemented with Cm. Sequence-confirmed plasmids were either conjugated to C. sporogenes-NT or used as a donor DNA in the CRISPR-Cas9 integration stages. Resultant transconjugants have been confirmed by colony PCR using primers specified in Table S2 and by Sanger sequencing (Azenta Life Sciences, formerly Genewiz, Inc.).

 Creation of two-plasmid CRISPR-Cas9 editing plasmids Construction and conjugation of the tetracycline inducible Cas9 vector

The aTc inducible Cas9 vector pTetR-P IPL12 -Cas9 is schematically illustrated in Figure 1A . Using pMTL83151 as template (54), a backbone fragment was created by polymerase chain reaction (PCR) using forward and reverse primers complementary to sequences upstream of pCB102 and downstream of traJ. The codon optimised tetR gene and divergent promoter sequences, consisting of PIPL12 and PminiP4_tU (53), were designed in silico as a fusion sequence with flanking BsaI restriction sites and commercially synthesised (IDT Inc.). Using two primer pairs, the cas9 gene was amplified from Streptococcus pyogenes M1 genomic DNA as two overlapping fragments, each containing 5’ BsaI restriction sites. In order to remove a BsmBI restriction site present in the native cas9 coding sequence, a synonymous codon change (CGT to AGA) was introduced into the reverse primer of the 5’ cas9 fragment. Two new terminators—derived from the Lactococcus lactis pepN gene and the B. subtilis tyrS gene—were incorporated downstream of cas9 and tetR, respectively. The terminators were split between adjacent plasmid fragments and incorporated as 5’ extensions to primers. The resulting four fragments were ligated using the Golden Gate Assembly protocol (NEB) to form the final plasmid. The inducible pTetR-P IPL12 -Cas9 vector was transferred by conjugation to each variant of C. sporogenes.

Schematic representation of the two-plasmid gene editing system. (A) pTetR-P IPL12 -Cas9 carries the native SpCas9, modified to remove a BsmBI restriction site. The tet repressor gene (tetR) is expressed constitutively and binds to the tet operator (tetO) sequence of the synthetic P IPL12 promoter in the presence of (anhydro)tetracycline. This is the basis for control of cas9 expression. Both vectors utilised backbones from pMTL8000 vectors (54); gRNA, guide RNA consisting of spacer and scaffold sequences; LHA/RHA: left/right homology arm; lacZα, β-galactosidase counter-selection marker (for cargo cloning). (B) A representative workflow for gene editing using the CRISPR-Cas system. Conjugative transfer of pTetR-P IPL12 -Cas9 into the host is conducted prior to and separately from that of p8222F-gX-HC. In parallel, a gRNA spacer sequence and the DNA (that is to be integrated) is cloned into p8222F-gX-HC. Following conjugative transfer of p8222F-gX-HC, viable colonies are restreaked to induction plates, which contain Erm, Tm, and aTc, and incubated for 24h. Cultures able to grow on induction plates are screened by PCR to detect genomic recombination. Plasmid loss is achieved by subculture in the absence of antibiotic selection. If Erm selection is maintained, pTetR-P IPL12 -Cas9 is retained by the host and additional rounds of gene editing can be conducted. (B) created with BioRender.com.

 Construction of knock-out vectors

To construct the non-targeting vector, p8222F, an inducible ferredoxin promoter (PfetO) was designed to include XbaI at the 5’-end and SalI-AatII-AscI cloning site at the 3’-end, and it was synthesized by IDT ( Figure 1A ; Table S3 ). For final construction, the inducible promoter was inserted into the pMTL82221 vector chassis via XbaI and AscI restriction sites (54).

gRNA protospacer sequences were selected using the CRISPOR web tool (55). Cloneable gRNA fragments were created using a conserved reverse primer (corresponding to the 3’ of the gRNA handle and terminator) and a forward primer containing the spacer sequence at the 5’ and the beginning of the gRNA handle at the 3’. Complementary sequences between these primers enabled the creation of a product in a template-free “primer dimer” reaction. A SalI and AatII restriction site were incorporated up- and downstream of the gRNA product through 5’ primer sequences. To generate editing templates (repair cassette), sequences of approximately 750 bp long upstream and downstream of the chosen location were amplified from the C. sporogenes NCIMB 10696 genomic DNA, incorporating type IIS restriction sites. Fragments designated as left and right homology arms (LHA/RHA) were assembled using Golden Gate method. To demonstrate and exemplify the conjugation efficiency and tight regulation of tetracycline induction, a p8222F-sg7 vector was constructed, containing only SLS-targeting sgRNA module, without repair cassette.

 Construction of Golden Gate compatible knock-in vectors

Two Golden Gate compatible knock-in vectors (p8222F-g7-SLS-GGL and p8222F-g2-PyrE-GGL) were constructed, targeting the SLS operon and pyrE gene, respectively ( Figure 1A ; Supplementary Figure S1 ). The Golden Gate-adapted lacZa module (56) was amplified and ligated between the repair cassette flanking SLS operon or pyrE gene, respectively.

To generate CRISPR-Cas9 editing plasmids with the desired cargo, p8222F-g7-SLS-GGL or p8222F-g2-PyrE-GGL were used in a GGA reaction alongside previously constructed pGG212 harbouring NfrA, mIL-2, or mGM-CSF, respectively, with BsmBI and T4 DNA ligase. Ligated plasmids were transformed by heat shock into E. coli 10-beta and plated on LB plates supplemented with Erm. Resultant colonies were PCR-screened for correct assembly and subsequent plasmid samples were confirmed by Sanger sequencing.

 Generation of <italic>C. sporogenes</italic> knock-out and integrant strains

Constructed knock-out and knock-in vectors were heat shock transformed into E. coli S17-1 and conjugated into C. sporogenes strain harbouring pTetR-P IPL12 -Cas9 plasmid. Conjugation was carried out as previously described (54) with some additional considerations. A critical ratio of 1: 5 (donor to recipient) was followed. Prior to conjugation, the recipient strain was heat shocked for 2 minutes at 55°C (57). Donor-recipient mixture was plated on PYT media supplemented with both Erm and Tm. Transconjugant colonies were subsequently re-streaked on PYT plates with the addition of aTc inducer. The procedure is schematically illustrated in Figure 1B .

Resultant knock-out or integrant strains were screened by colony PCR using primers that flank the genome locus ( Table S2 ). The editing vector and pTetR-P IPL12 -Cas9 vector were removed from knock-out mutants by subculturing in non-selective liquid cultures (24-48h) followed by patch plating to selective and non-selective plates to observe loss of antibiotic resistance. Plasmid loss was further confirmed by PCR specific for the plasmid-based traJ gene (data not shown). Resultant positive clones were sent for sequencing. In the event where resultant recombinant strains were subjected to another round of integration with a different editing template, the pTetR-P IPL12 -Cas9 vector was maintained.

Measurement of conjugation efficiency

To calculate conjugation efficiency (CE), donor-recipient spot mixtures were incubated for 7 hours on non-selective plates before being resuspended in 0.5 ml PBS, serially diluted, and spread on selection plates containing Cs and Erm. When using C. sporogenes bearing the pTetR-P IPL12 -Cas9 as the recipient, Tm was added to plates to maintain the Cas9 plasmid in transconjugants. Additionally, the conjugation mixtures were also plated on relevant selection plates with the addition of aTc as the inducer (96 ng/ml). Approximate conjugation efficiencies were then calculated as transconjugant colony forming unit (CFU) per total C. sporogenes recipient CFU.

 Phenotype analysis of <italic>C. sporogenes</italic> knock-out and integrant strains Growth and endospore formation for integrant strains

The ability of spores of C. sporogenes integrants to return to vegetative growth was evaluated by measuring the change in OD600 over a period of 24 hours. Overnight cultures were inoculated into a fresh PYT media (1:100) and incubated at 37°C under anaerobic conditions. Additionally, uracil auxotrophy for C. sporogenes strains lacking the pyrE gene was evaluated by plate assay using agar plates with a chemically defined composition of media and supplemented with uracil (20 μg/ml) where appropriate. Endospore formation titres for C. sporogenes integrant strains were evaluated at 120h time point. Briefly, 1 ml samples were treated for 20 min at 80 ˚C to heat-kill vegetative cells and plated in serial dilutions (from 10-1 to 10-7 in sterile and anaerobically reduced PBS) on PYT plates supplemented with Cs. After 24-48h anaerobic incubation, colonies from agar plates were enumerated and spore titre calculated as CFU/ml. A spo0A mutant (CspWT-Δspo0A), in which the master regulator of sporulation has been deleted, was used as a negative control for colony formation after heat treatment.

 Detection of recombinant mIL-2 and murine GM-CSF

To quantify the levels of cytokines secreted from either plasmid-based or integrant strains of C. sporogenes, supernatant samples were obtained from three different growth time points (5h, 7h and 10h) during bacterial growth. Subcultures of mIL-2 or mGM-CSF-expressing C. sporogenes were centrifuged (10,000 x g, 10 min), filtered (0.22 μm syringe filter), and the sterile supernatant retained.

The levels of mIL-2 and mGM-CSF were determined by cytokine-specific enzyme-linked immunosorbent assay (ELISA) in accordance with the manufacturer’s instructions (88-7334; 88-7334, Invitrogen). The results were recorded in a microplate reader (iD3, SpectraMax) and calculated based on recombinant mIL-2 and mGM-CSF standards and normalised to OD600 = 1 to account for differences in growth rates.

The biological activities of the secreted cytokines were determined in a lymphocyte proliferation assay using the CTLL-2 or FDCP-1 T cell line. Supernatant-stimulated proliferation was measured using either an MTT assay (for CTLL-2) according to published methods (58) or AlamarBlue (for FDCP-1). The standard curve was prepared using purified recombinant mIL-2 or mGM-CSF (212-12 and 315-03, Peprotech). For the MTT assay, absorbance was recorded at 470 nm. AlamarBlue assay results were recorded using the fluorometric read-out at 560/590 nm excitation/emission, in a multimode microplate reader (iD3, SpectraMax).

 Menadione reductase assay

To compare the activity of NfrA nitroreductase expressed in plasmid-based and integrant strains, menadione reductase activity was determined by the method of Knox et al. (59). To do this, 7-hour samples of 2 ml cultures were harvested by centrifugation. Cell pellets were subjected to BugBuster cell lysate preparation according to manufacturer’s recommendations. The assay was performed in 1 ml plastic cuvettes and a spectrophotometer set for 37°C. The increase of absorbance at 550 nm for 1 min was recorded and the rate divided by the volume of lysate (in ml). The extinction coefficient (EmM) of cytochrome c was used at 14.79. To normalise the units of activity per weight of protein, a bicinchoninic acid assay (BCA) was performed in accordance with manufacturer’s guidelines (23225, Thermo Scientific). The final menadione NfrA activity was expressed in units per mg (U/mg).

 Gelatin zymography and quantitative analysis of proteolytic activity

To detect and compare the levels of proteases in wild-type or recombinant C. sporogenes strains, gelatin zymography was performed using 10% Novex® Zymogram Gel according to manufacturer’s instructions (MAN0005885, Novex/Invitrogen). To do this, samples from cultures of C. sporogenes-NT-ΔpyrE and C. sporogenes-NT-ΔpyrEΔPR1 were taken at four different time points (5h, 7h,10h, and 12h), centrifuged (10,000 x g, 10 min) and filter sterilised (0.22 μm syringe filter). Sterile supernatants, normalised according to cell density (OD600 = 1.2), were mixed with Tris-Glycine SDS Sample Buffer (2x), of which 20 μl was loaded under non-reducing conditions into precast gelatin zymography gels. Clostridium butyricum DSM 10702 supernatant sample, obtained from 10-hour growth, was used as a control. The zymograph was run at 125 V for 90 minutes.

To quantify the total protease activity in culture supernatants, a colorimetric protease assay was performed in 96-well plates, according to manufacturer’s recommendations (23263, Thermo Scientific™ Pierce™ Colorimetric Protease Assay Kit). The spectrophotometric measurement was recorded using a microplate reader (iD3, SpectraMax). Comparative protease activities were calculated and expressed in ng/ml using trypsin standard curve as a reference.

 Biological replicates and statistical analyses

All data presented in this manuscript represent the results of at least two independent experiments. Statistical evaluations were performed with GraphPad Prism 9 software (San Diego, CA, United States). For the ELISA, cell proliferation and protease assays, and for the conjugation efficiencies, data were analysed using unpaired t-test to compare strain variants analysed at different time points or conditions. Values of P<0.05, P<0.01, P<0.001 were considered significant (∗), highly significant (∗∗), or extremely significant (∗∗∗) respectively. Data represent means ± standard deviation (SD).

 Results Two-plasmid, inducible CRISPR-Cas9 system improves editing efficiency in <italic>C. sporogenes</italic>

In the present study, we have constructed a genome editing system based on two plasmids: a cas9-expression plasmid and a second plasmid bearing the gRNA and repair cassette. Expression of the cas9 gene is controlled using a tetracycline (tet) inducible promoter (60). A schematic representation of the plasmid is shown in Figure 1A . Prior to the creation of this plasmid, gene expression control using the tet repressor and operator was confirmed in E. coli using the gusA (β-glucuronidase, GUS) reporter ( Supplementary Figure S2 ).

In order to test the new system, a series of plasmids, designed to knock-out sequences of the C. sporogenes genome, were created. Using the labelling scheme p8222F-gX-HC, these plasmids contain the gRNA (“X”) under the control of the tet inducible Pfdx promoter (PfetO) and the target-specific homology cassette (“HC”). Four genome loci were targeted: pyrE (609 bp), spo0A (1.2 kb), spoIIE (2.4 kb), and the SLS operon (8.64 kb).

Prior to mutagenesis, the pTetR-P IPL12 -Cas9 plasmid was conjugated to the C. sporogenes strain of interest (C. sporogenes WT or C. sporogenes-NT). These cas9-bearing cultures were then used as conjugation recipients of the second plasmid. Following induction, deletion of the target sequence was detected in all clones screened by flank PCR, indicating 100% efficiency ( Supplementary Figures S3A, B ).

Poor conjugation efficiency has the potential to distort reported genome editing efficiency. To demonstrate that perceived high genome editing efficiency is not an artefact of low transconjugant yield, we sought to demonstrate the conjugation performance of this system in the presence and absence of inducer. This was achieved by normalising transconjugant CFUs based on recipient titre, referred to henceforth as conjugation efficiency (CE). In addition, we evaluated the effect of cargo size on p8222F CE. Plasmids bearing the repair cassette and/or sgRNA were targeted to the SLS locus. Following conjugation, CE was determined both for empty C. sporogenes and for C. sporogenes bearing pTetR-P IPL12 -Cas9. In order to assess the effect of Cas9 and gRNA expression on conjugation in the absence of a repair cassette, a plasmid bearing the SLS gRNA, but without a repair cassette, was constructed (p8222F-sg7).

Conjugation of p8222F, in which gRNA and repair cassette are absent, showed similar efficiency in both recipients, indicating the presence of pTetR-P IPL12 -Cas9 does not affect conjugation ( Figure 2 ). The addition of inducer did not influence p8222F CE for either recipient. The addition of gRNA (p8222F-g7) caused a non-significant decrease of CE in C. sporogenes pTetR-P IPL12 -Cas9, which was reduced further in the presence of inducer (P=0.002). This was not observed when empty C. sporogenes was used as the recipient. The addition of the repair cassette (p8222F-g7-SLS) appears to “rescue” CE in C. sporogenes pTetR-P IPL12 -Cas9, with comparable efficiencies in the presence and absence of inducer. The inclusion of 3 kb of “cargo” (p8222F-g7-SLS-3kb) reduced CE, with or without inducer, in C. sporogenes pTetR-P IPL12 -Cas9 but the difference was not statistically significant.

The effect of the size and the composition of four different vectors on the conjugation efficiencies of two recipients: empty C. sporogenes WT and strain bearing the inducible Cas9-plasmid. p8222F (4,998 bp): plasmid with only P fetO promoter; p8222F-g7 (5,128 bp): plasmid with addition of sgRNA targeting SLS operon; p8222F-g7-SLS (6,567 bp): plasmid targeting SLS operon with repair cassette; p8222F-g7-SLS-3kb (9,698 bp): plasmid targeting SLS operon with repair cassette and 3 kb λDNA cargo. The transconjugants were counted on selection plates: Erm/Cs for C. sporogenes WT (in turquoise) and Erm/Cs/Tm for C. sporogenes-pTetR-P IPL12 -Cas9 (in orange) with and without inducer (aTc, 96 ng/ml). Approximate conjugation efficiencies were calculated as obtained transconjugant CFU/total C. sporogenes recipient CFU. Each bar represents the mean and standard deviation of data collected from two experiments performed using biological duplicates. Statistical analysis was performed using unpaired t-test, ns – not significant, ** P<0.01.

 Golden Gate-compatible vectors streamline cargo cloning

In order to streamline cargo cloning, plasmids targeted to two genome loci were engineered to contain golden gate cloning sites between the homology arms (56), referred to as p8222F-g2-PyrE-GGL and p8222F-g7-SLS-GGL. These plasmids target the pyrE locus and SLS locus, respectively. Using GGA, NfrA, mIL-2 and mGM-CSF were cloned into p8222F-g7-SLS-GGL. Based on the results of previous work, the native promoter of the ferredoxin (Pfdx) gene was selected to ensure strong expression of the recombinant proteins (53), and the native secretion peptide of nprM3 to promote efficient secretion (51). Using C. sporogenes WT and C. sporogenes-ΔpyrE as hosts, knock-in mutants of NfrA, mIL-2 and mGM-CSF were created, resulting in C. sporogenes-NT and C. sporogenes-NT-ΔpyrE versions of each recombinant strain ( Supplementary Figure S4A ). Deletion of the pyrE gene has been shown to cause uracil auxotrophy (52), and could serve as a method for biocontainment of genetically engineered C. sporogenes when used in clinical applications. Although we have only utilised the SLS-targeting plasmid in this study, the creation of the p8222F-g2-PyrE-GGL plasmid enables integration at the pyrE locus in a second round of genome editing, using C. sporogenes-NT as the host. Other loci can be targeted using the same rational.

In order to test the ability of the CRISPR-Cas9 system to integrate large DNA fragments, two arbitrary fragments of lambda DNA were created and cloned between the homology arms. The largest cloned fragment was 5 kb, resulting in a total plasmid size of 11.7 kb. Larger fragments were not tested as the efficiency of both E. coli transformation and C. sporogenes conjugation is low for plasmids above 12 kb. Based on total colonies screened, integration efficiencies were 91.66-100% for cargo sizes ranging from 835 bp to 5,256 bp ( Supplementary Figure S4B ).

 Chromosome-based intracellular and secreted products are biologically active

Clones of six recombinant strains were taken forwards for biochemical and functional validation. The recombinant strains were generated in two different backgrounds: ΔSLS and ΔSLSΔpyrE. Growth on defined media agar plates with or without uracil was tested ( Supplementary Figure S5 ). Recombinant strains generated in the ΔSLSΔpyrE background were not able to grow in the absence of uracil, indicating they are auxotrophic. In contrast, strains generated in the ΔSLS background did not require uracil supplementation. To further characterise the new strains, growth over 24 hours in rich media and ability to form spores over 120 hours was tested. Growth rate and sporulation of both the recombinant strains and the knock-out strains were somewhat similar to those of C. sporogenes WT ( Figures 3A–C ). The results of these assays indicate that addition of heterologous expression cassettes does not have a significant effect in these processes under the conditions tested.

Characterisation of growth, sporulation, and enzymatic activity assay in C. sporogenes strains. (A) Growth of Clostridium strains was measured as a direct increase in absorbance at 600 nm throughout the course of 24-hour bacterial incubation in bovine free medium (PYT). The symbols represent the average of three independent experiments, and error bars indicate the standard errors of the means. (B) The summary of OD600 measurements recorded at the times of pre-determined experimental time point sample collections (5h, 7h and 10h). (C) Spore titres after 120h incubation in PYT broth. Samples of heat-treated cultures (80°C, 20 min) were plated in serial dilution on agar plates and enumerated following 24-48h incubation. Bars represent the number of CFU (colony forming unit) per ml of culture. The data represent the average of three independent experiments and error bars indicate the standard error of the mean. The sporulation-deficient Clostridium sporogenes-NT-Δspo0A mutant was used as a negative control to rule out experimental error. The detection limit for colony counts was 50 CFU/ml. (D) Menadione reductase activity assays on C. sporogenes cell lysates containing either plasmid-based or chromosome-integrated nfrA gene. Results show specific menadione reductase activities of C. sporogenes NfrA-expressing strains soluble fractions obtained from 7h time points during the growth of bacterial cultures. The C. sporogenes-NT-ΔpyrE was used as background control. Activities were normalised to total protein concentrations, determined using BCA assay. Activities are expressed in units per mg. The data represent the average of three independent experiments and error bars indicate the standard errors of the means. Statistical analysis was performed using unpaired t-test, *** P<0.001.

The nitroreductase enzyme NfrA is expressed intracellularly and can be quantified using the menadione assay, as reported previously (52). Recombinant strains expressing NfrA from a plasmid or the chromosome were generated in the ΔSLS and ΔSLSΔpyrE backgrounds. Host genotype did not affect nitroreductase activity, but plasmid-based strains produced significantly more activity than chromosome-based equivalents (P<0.0004) ( Figure 3D ).

Levels of secreted mIL-2 and mGM-CSF in culture supernatants after 5, 7, and 10 hours of growth were quantified by ELISA. Supernatants from cultures of equivalent plasmid-based constructs were taken at the same time points for comparison. Cell density was measured for each sample at each time point ( Figure 3B ), and calculated cytokine levels were normalised to OD600 = 1. In strains of both cytokines, maximum product levels occurred after 5 hours of growth, with sequential reductions at 7 and 10 hours ( Figure 4A ). Product levels were substantially (more than 10-fold) higher from plasmid-based expression cassettes at all time points compared to the same cassette integrated into the chromosome.

Validation of C. sporogenes-NT-ΔpyrE strains harbouring either plasmid-based or integrated murine cytokines (mIL-2 and mGM-CSF) at the SLS locus in the quantitative and functional assays. (A) Results of commercial ELISA tests indicating the quantities of mIL-2 and mGM-CSF cytokines present in the supernatants of cytokine bearing-C. sporogenes strains at 5-,7-, and 10-hour growth in PYT media. (B) Results of MTT and AlamarBlue functional assays following the incubation of cytokine specific T-cells in the presence of culture supernatants of C. sporogenes strains. Recombinant murine IL2 and GM-CSF were used to prepare mIL-2 and mGM-CSF standard curves., (C) denotes C. sporogenes-NT-ΔpyrE control, (BDL) - below detectable levels. Statistical analysis was performed using unpaired t-test, * P<0.05, ** P<0.01, *** P<0.001.

Cytokine function was evaluated in cell proliferation assays using the CTLL-2 and FDCP-1 cell lines for mIL-2 and GM-CSF, respectively ( Figure 4B ). In line with the ELISA results, peak levels of mIL-2 (3.9 µg/ml) and mGM-CSF (2.5 µg/ml) occurred in cultures of plasmid-based strains after 5 hours of growth. At the same time point, chromosome-based constructs produced substantially less product (0.5 µg/ml and 0.04 µg/ml, respectively). For the plasmid-based strain, levels of mIL-2 were reduced at 7 hours (-35%) and almost undetectable at 10 hours (-98%). In contrast, levels of mGM-CSF from the plasmid-based strain remained relatively high at all time points, reducing sequentially by 32% and 50% at 7 and 10 hours, respectively.

 Deletion of the 7 kb nprM operon increases the yield of secreted cytokines

Levels of secreted cytokine were significantly reduced in later phases of bacterial growth cultures. This led us to hypothesise that proteases secreted by the host bacteria were degrading the recombinant product. To test our hypothesis, we sought to investigate the impact that deleting proteolytic genes would have on the levels of secreted recombinant product. A bioinformatic screen of the C. sporogenes NCIMB 10696 genome revealed a 7.4 kb operon, annotated as six metallopeptidase genes (1,607,357 bp to 1,615,014 bp), henceforth referred to as the nprM operon. Using the same method as for pyrE and the SLS operon, the six genes were deleted using gRNA and repair cassettes targeted to the nprM operon ( Figure 5A ). This was conducted in C. sporogenes-NT-ΔpyrE and in the cognate mIL-2 and GM-CSF expressing strains. Deletion of the target operon (“ΔPR1”) was confirmed by PCR and Sanger sequencing ( Figure 5B ). The loss of these genes did not affect the host’s ability to grow or sporulate in rich media ( Figures 3A–C ).

Schematic representation of nprM proteolytic operon knock-out in C. sporogenes mediated by inducible two-plasmid CRISPR-Cas9 system. (A) Deletion of six proteolytic genes has been exemplified with C. sporogenes-NT-ΔpyrE strain and subsequently repeated using two cytokine-bearing integrant strains. Triangles represent indicative alignment of screening primers as presented in Table S2 . (B) Gel electrophoresis showing PCR screens of all three proteolytic knock-out strains conventionally labelled as ΔPR1, where S1: CspNTΔPΔPR1 (C. sporogenes-NT-ΔpyrEΔPR1), S2: CspNTΔPΔPR1::mIL-2 (C. sporogenes-NT-ΔpyrEΔPR1::mIL-2) and S3: CspNTΔPΔPR1::mGM-CSF (C. sporogenes-NT-ΔpyrEΔPR1::mGM-CSF). Left image: Deletion of the operon detected in colonies subjected to a PCR reaction using primers flanking the deletion locus. Right image: Subsequent confirmation of proteolytic operon deletion in selected clones using internal (junction) and flanking primers. (L) denotes DNA marker, (C1-C3): equivalent C. sporogenes-NT-ΔpyrE controls.

In order to visualise the effect of the nprM operon deletion on the levels of secreted proteins, a gelatin zymograph experiment was conducted ( Figure 6A ). For each time point (5, 7, 10 and 12-hour), sterile supernatants of CspNT-ΔpyrE and CspNT-ΔpyrEΔPR1 were loaded, pairwise, onto the zymograph. Following electrophoresis, we observed a distinct difference at all time points between the two strains. In relative terms, the CspNT-ΔpyrE samples showed reduced quantities of high molecular weight proteases, but distinct lower molecular weight proteases that are absent in the CspNT-ΔpyrEΔPR1 samples.

Analysis of proteolytic activity in C. sporogenes strains. (A) Qualitative analysis of proteolytic activity by gelatin zymography of the C. sporogenes-NT-ΔpyrE and C. sporogenes-NT-ΔpyrEΔPR1 supernatants at four different time points (5h, 7h, 10h and 12h). In total, 20 μl of samples were loaded on a SDS polyacrylamide gel containing 0.05% of gelatin. After electrophoresis, proteins were renaturated in the gel and the proteolytic activity was revealed after an overnight incubation required for gelatin degradation. CB10 - C. butyricum control supernatant collected at 10h growth. Samples were normalised according to OD600 of 1.2 prior to loading. M: Pre-stained Protein Ladder PageRuler™ (Thermo Scientific™); P+: C. sporogenes-NT-ΔpyrE samples (PR1 positive); P-: C. sporogenes-NT-ΔpyrEΔPR1 (PR1 negative); CB10 - C. butyricum control supernatant collected at 10h growth. (B) Quantitative analysis of proteolytic activities detected in C. sporogenes strains using Pierce™ colorimetric protease assay (Thermo Scientific™) and sample supernatants collected at three different time points (5h, 7h, and 10h). CspNTΔP (P+) and CspNTΔPΔPR1 (P-) – PR1-positive and PR1-negative backgrounds, respectively. Standard curve for the assay has been constructed using serially diluted TPCK trypsine stock solution. The absorbance measurement has been carried out at 450 nm. The data represent the average of two independent experiments, and error bars indicate the standard errors of the means. Statistical analysis was performed using unpaired t-test, *** P<0.001.

A colorimetric protease assay was conducted to quantify the impact of this deletion on extracellular proteolysis ( Figure 6B ). The Pierce Protease Assay detects proteolysis through the reaction of peptide fragments, created by hydrolysis of succinylated casein, with trinitrobenzene sulfonic acid (TNBSA). The resultant TNB-peptide adducts produce a yellow colour, which can be quantified by measuring absorbance at a 450 nm wavelength and comparison with a protease calibration curve. As for previous experiments, supernatants were collected at 5, 7, and 10h for three strains, mIL-2, GM-CSF, and parental, in two different backgrounds: ΔpyrE and ΔpyrEΔPR1. This experimental design allowed us to compare proteolysis in different contexts and across time and to do multiple pairwise analyses of proteolytic activity. A statistically significant reduction in proteolysis was observed in ΔpyrEΔPR1 strains compared to ΔpyrE strains at all time points (P<0.05). This held true for strains with and without recombinant products.

Having established that the nprM knock-out strains showed reduced proteolysis, we sought to test our hypothesis that this phenotype would translate to increased recombinant product in culture supernatants. The ELISA and cell proliferation assays, detailed in Figure 4 , were repeated using the CspNT-ΔpyrEΔPR1 expressing mIL-2 and GM-CSF. Equivalent strains in the CspNT-ΔpyrE background (used in the previous experiment) were tested again in these assays. At all time points, recombinant product was increased in the CspNT-ΔpyrEΔPR1 supernatants ( Figures 7A, B ). The difference was less pronounced at the earliest time point (5h), with increased ELISA/cell assay measurements of 244%/148% and 248%/142% for mIL-2 and GM-CSF, respectively. The difference in product levels was more pronounced at 7h but increased further to peak at 10h with increased ELISA/cell assay measurements of 1653%/1546% and 807%/363% for mIL-2 and GM-CSF, respectively. In all cases, the differences were statistically significant (P<0.05).

Validation and comparison of C. sporogenes-NT-ΔpyrE and C. sporogenes-NT-ΔpyrEΔPR1 strains harbouring integrated murine cytokines (mIL-2 and mGM-CSF) in the quantitative and functional assays. (A) Results of commercial ELISA tests indicating the quantities of mIL-2 and mGM-CSF cytokines present in the supernatants of cytokine bearing-C. sporogenes strains at 5-,7-, and 10-hour growth in PYT media. (B) Results of MTT and AlamarBlue functional assays following the incubation of cytokine specific T-cells in the presence of culture supernatants of C. sporogenes strains. Recombinant murine IL-2 and GM-CSF were used to prepare mIL-2 and mGM-CSF standard curves. (C) denotes C. sporogenes-NT-ΔpyrEΔPR1 control, (BDL) – below detectable levels. The data represent the average of three independent experiments and error bars indicate the standard errors of the means. Statistical analysis was performed using unpaired t-test, ** P<0.01, *** P<0.001.

 Discussion

The objectives of this study were two-fold: 1) to combine several improvements to CRISPR-Cas9 technology that result in the development of a cloning-friendly and efficient genetic tool for modification of Clostridium genomes and 2) to demonstrate the application of this tool in the development of oncologically relevant recombinant strains. To achieve this, a number of impediments had to be overcome. Although less of an issue in C. sporogenes, poor efficiency of transformation, by conjugation or electroporation, has been a major hindrance to strain development in related Clostridium species (61). Irrelevant of the transformation method, the inverse relationship between efficiency and plasmid size is difficult to overcome. To mitigate this, we chose to split the components of the CRISPR-Cas9 system across two plasmids. We constructed a 9.2 kb Cas9 expression plasmid (pTetR-P IPL12 -Cas9) and a second, 6.6 kb editing plasmid (p8222F-gX-HC), containing the gRNA and repair cassette. The relatively small size of p8222F-gX-HC leaves significant headroom for the addition of sequences to be integrated before the overall size of the plasmid hinders cloning and conjugation. Transfer of DNA is considerably less efficient in other species of clostridia and is likely to necessitate further measures to facilitate transformation. Improvements have been achieved by selective DNA methylation to bypass host restriction-modification systems (62) and through heat-shock of the recipient prior to conjugation (57). The equivalent single plasmid system would require a plasmid of 11 kb or larger. An additional advantage of this approach is that multiple rounds of genome editing can be conducted following conjugation of the pTetR-P IPL12 -Cas9 plasmid without the need to remove the Cas9 expression vector. To reduce the time and resources required to clone new sequences for integration, we incorporated the previously reported golden gate cloning site (56). This enables one-pot cloning reactions and blue-white colony screening. Using previously validated integration sites, such as those demonstrated in this paper, this system can be used to create libraries of integrant strains at a cost of time and resources that most academic labs could bear.

Obtaining transconjugants that bear both plasmids is significantly easier when Cas9 expression is repressed. Previous versions of this system relied on constitutive expression of gRNA or the cas9 gene or both (63). By placing both Cas9 and gRNA expression under the control of inducible promoters, the inherent toxicity of Cas9, with or without gRNA, is alleviated. However, the reduced conjugation efficiency of all p8222F-based plasmids into a pTetR-P IPL12 -Cas9-bearing host compared to plasmid-less host indicates that a degree of host toxicity remains, possibly due to incomplete repression of cas9 transcription.

In this context, utilising inducible promoters for Cas9 and/or gRNA expression risks insufficient expression of these components, and a concomitant loss of potency. Ultimately this would result in a low frequency of mutant isolation, due to inadequate Cas9-mediated DNA cleavage and subsequent homology-directed repair. Assessment of a panel of tet inducible promoters in a GUS reporter assay indicated that expression from synthetic promoter IPL12-miniP4_tU was high in the presence of inducer, but indistinguishable from background in its absence ( Supplementary Figure S2 ). Following conjugation and induction, our results indicate that all surviving colonies had recombined with the provided repair cassette to produce the desired mutant strains. This held true in tests with cargo sizes up to 5 kb. Combined with the data for CE and gene deletion, these results demonstrate that deletion of sequences up to 8.6 kb and addition of at least 5 kb can be achieved by a single user in a short amount of time without the need for high volume colony screens.

Chromosomal integration of foreign DNA offers several advantages compared to plasmid-based approaches. Transfer of extrachromosomal DNA between bacteria in multispecies communities is well documented, and maintenance of a plasmid typically requires the presence of antibiotic. The risk of introducing synthetic plasmids to other bacteria, within a patient or the environment, is likely to be a major hurdle for regulatory approval. Furthermore, plasmids are structurally and segregationally unstable, which can lead to unpredictable gene expression or inactivation of the gene. However, plasmid-based genes can achieve higher levels of expression than when they are chromosome-based.

We selected two oncologically relevant cytokines to assess the capability of the new gene editing system and to evaluate the impact of integration on biological activity. Interleukin-2, a potent T cell activator, is approved for the treatment of renal cell carcinoma and metastatic melanoma (64, 65), where high-dose regimens produce significant clinical responses (66) but are hindered by severe toxicity (67). Attempts to prevent toxicities using lower dosage regimens can cause a significant drop in therapeutic effect due to the dominant effect of immunosuppressive regulatory T cells (68). “Immune-related adverse effects” (irAEs) during immunotherapy are the result of the agent interfering with normal immune function in healthy tissues, and typically manifests as autoimmune-like symptoms. Agents that are hindered or, in some cases, unusable due to irAEs are ideal candidates for bacteria-mediated intratumoural delivery, which could, in theory, maximise the dose at the tumour site without requiring systemically high doses. As for IL-2, the influence of GM-CSF in anti-cancer immunity is dose-dependent. GM-CSF stimulates the production, maturation, and activation of innate immune cells, such as neutrophils, macrophages and dendritic cells. Thus, it is postulated that GM-CSF treatment could promote activation of the adaptive immune response indirectly by promoting tumour-reactive innate cells (69). In addition, strains expressing a pro-drug converting nitroreductase enzyme (NfrA), reported previously (52), were developed. NfrA is currently being investigated for its potential to convert PCE CP-506, entering phase 1/2 clinical trial (70). Due to the absence of a secretion peptide in the coding sequence of NfrA, this product remains inside the host cell. Characterisation of the secreted cytokines and the intracellular NfrA enabled a comparison of plasmid and chromosomal expression in two different contexts. Although the menadione reductase activity levels detected from chromosomal strains are lower than in other already investigated NTR-expressing strains (such as NmeNTR), the advantage of clinically relevant CP-506 prodrug for CDEPT approach might be of a bigger value (52).

We sought to quantify the production and secretion of the cytokines biochemically and to determine functional activity in cell culture assays. Furthermore, the effect of bacterial culture duration (as a function of growth phase) on cytokine production was assessed. For both cytokines, ELISAs indicated that plasmid-borne expression cassettes produced substantially higher levels of product in culture supernatants. By sampling cultures at multiple time points, we were able to track secretion of the cytokines across growth phases ( Figure 3B ). In previous studies, supernatants were analysed at a single time point (7-hour), which represents late-exponential phase growth in C. sporogenes (51). However, in the present study we sampled at three timepoints: 5, 7, and 10 hours, representing mid-exponential, late-exponential, and stationary growth phases. Regardless of the product or genetic construct, the highest product titres occurred at the 5-hour time point and were lower at 7 hours, decreasing further by 10 hours. An ELISA relies on target-specific antibodies and can detect the presence of correctly folded protein in its native conformation. However, it is not a direct measurement of protein function. We reasoned that testing the functionality of the recombinant cytokines in cell culture assays would yield more meaningful results. Proliferation of the cytotoxic T cell line CTLL-2, and the myeloid progenitor cell line FDCP-1, were used to quantify mIL-2 and mGM-CSF, respectively, in culture supernatants. These assays revealed a similar trend: plasmid-based constructs produced higher levels of proliferation and peak titres occurred at the earliest timepoints, diminishing over culture time.

Increased expression of plasmid-encoded products compared to those that are chromosomally integrated is widely reported in the literature (71). This is attributed to the presence of multiple copies of the plasmid in each cell, which enables a higher rate of transcription. However, the cause of the substantial variation in product titre between time points, for plasmid- and chromosome-based products of both cytokines, was unclear. We sought to investigate the cause of the decline in cytokine levels at later time points. We rationalised that there were two major factors that could influence yield: production and degradation. Some bacterial gene promoters have been shown to vary according to growth phase, such as that of the ptb gene of Clostridium acetobutylicum (72). For the promoter in question (Pfdx) we did not find equivalent reports of growth-phase dependent expression.

We turned our attention to potential causes of degradation. The genome sequence of C. sporogenes reveals several coding sequences with homology to characterised B. subtilis protease genes, including several matrix metalloproteases and bacillolysin, all of which are secreted (73). Protease secretion from C. sporogenes occurs in the late exponential phase and is repressed by glucose, ammonia, phosphate, ATP/ADP and certain amino acids, suggesting protease secretion occurs when energy sources are depleted. Previous studies have speculated that it is the proteolytic capability of C. sporogenes that enables it to effectively colonise solid tumours (74, 75). Several studies have demonstrated that genetic manipulation of secreted proteases to increase or decrease proteolysis has been demonstrated in several studies, including in B. subtilis (76, 77) and L. lactis (78). We hypothesised that co-secretion of native proteases was adversely affecting the stability of our recombinant product, which would explain the reduction of cytokine activity at later time points. To test this hypothesis, we identified and deleted an operon of six genes in the C. sporogenes NCIMB 10696 genome related to proteolysis. Following deletion of this operon in the cytokine-expressing strains, the quantitative and functional assays were repeated. Increased levels of cytokine were observed in all samples but was most pronounced in samples from later time points, corresponding to late exponential and stationary phases, growth stages in which protease expression and secretion increases due to nutrient depletion.

Through the creation of an effective CRISPR-Cas9 system, this study illustrates the potential of recombinant C. sporogenes strains as drug delivery vehicles. As our collective knowledge of the biology and genetics of clostridia increases, alongside the creation of effective genetic tools that are genuinely accessible, our ability to fully exploit these species in pursuit of medical applications will be realised. While the results of this study are promising, they only demonstrate efficacy in the simple, contrived setup of in vitro cell assays. The logical next step is to determine efficacy in vivo, that is, the therapeutic efficacy of Clostridium-mediated immunotherapy in comparison to current standard practice. Selection of potent, “next generation” agents and high-throughput experiments to further optimise expression and secretion will increases the chance of success. However, trying to gauge efficacy by comparison to systemically administered agents may be futile. Perhaps it is time to leap into the unknown.

 Data availability statement

The original contributions presented in the study are included in the article/ Supplementary Material , further inquiries can be directed to the corresponding authors. Plasmids pTetR-P IPL12 -Cas9 and p8222F-g7-SLS-GGL will be made available through Addgene following publication.

Ethics statement

Ethical approval was not required for the studies on animals in accordance with the local legislation and institutional requirements because only commercially available established cell lines were used.

 Author contributions

AK, LC, YZ, and TB co-designed the study. AK, LC, YZ, and TB performed the experiments and collected and analysed the data. AK, TB, and KK wrote the manuscript. KK assisted in the conceptualisation and evolution of ideas. AK and TB initiated the project and share senior authorship. TB wrote the funding grant. All authors contributed to data interpretation, discussions, and revision of the manuscript. All authors contributed to the article and approved the submitted version.

 Funding

This project received funding from KWF grant number 13056.

 Acknowledgments

The authors thank Natasja Lieuwes, Rianne Biemans, Jolanda Piepers and Kim Savelkous for technical support in the laboratory. The authors thank Professor Nigel Minton (SBRC, University of Nottingham, UK) for training in Clostridium research techniques, scientific guidance and support, past and present.

Conflict of interest

AK was an employee of Exomnis Biotech B.V. at the time of the study design and execution. AK is an inventor on two patents pending WO2021123391 and WO2022171904. TB is an inventor on one patent pending WO2022171904. LC was an employee of Exomnis Biotech B.V. at the time of the study design and execution 2021/2022.

The remaining 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.

 Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2023.1241632/full#supplementary-material

References Coley WB . The treatment of Malignant tumors by repeated inoculations of erysipelas: with a report of ten original cases. Am J Med Sci (1893) 105(5):487–516. doi: 10.1097/00000441-189305000-00001 Dranoff G . Cytokines in cancer pathogenesis and cancer therapy. Nat Rev Cancer (2004) 4(1):11–22. doi: 10.1038/nrc1252 Schirrmacher V van Gool S Stuecker W . Breaking therapy resistance: An update on oncolytic newcastle disease virus for improvements of cancer therapy. Biomedicines (2019) 7(3):66. doi: 10.3390/biomedicines7030066 Schirrmacher V van Gool S Stuecker W . Counteracting immunosuppression in the tumor microenvironment by oncolytic newcastle disease virus and cellular immunotherapy. Int J Mol Sci (2022) 23(21):13050. doi: 10.3390/ijms232113050 Li K Zhao Y Hu X Jiao J Wang W Yao H . Advances in the clinical development of oncolytic viruses. Am J Transl Res (2022) 14(6):4192–206. Leach DR Krummel MF Allison JP . Enhancement of antitumor immunity by CTLA-4 blockade. Science (1996) 271(5256):1734–6. doi: 10.1126/science.271.5256.1734 Juneja M Kobelt D Walther W Voss C Smith J Specker E . Statin and rottlerin small-molecule inhibitors restrict colon cancer progression and metastasis via MACC1. PLoS Biol (2017) 15(6):e2000784. doi: 10.1371/journal.pbio.2000784 Martin-Orozco N Wang YH Yagita H Dong C . Cutting Edge: Programmed death (PD) ligand-1/PD-1 interaction is required for CD8+ T cell tolerance to tissue antigens. J Immunol (2006) 177(12):8291–5. doi: 10.4049/jimmunol.177.12.8291 Dong F Wang L Davis JJ Hu W Zhang L Guo W . Eliminating established tumor in nu/nu nude mice by a tumor necrosis factor-alpha-related apoptosis-inducing ligand-armed oncolytic adenovirus. Clin Cancer Res (2006) 12(17):5224–30. doi: 10.1158/1078-0432.CCR-06-0244 Kamada T Togashi Y Tay C Ha D Sasaki A Nakamura Y . PD-1(+) regulatory T cells amplified by PD-1 blockade promote hyperprogression of cancer. Proc Natl Acad Sci U.S.A. (2019) 116(20):9999–10008. doi: 10.1073/pnas.1822001116 Bejarano L Jordao MJC Joyce JA . Therapeutic targeting of the tumor microenvironment. Cancer Discov (2021) 11(4):933–59. doi: 10.1158/2159-8290.CD-20-1808 Hu CW Chang YC Liu CH Yu YA Mou KY . Development of a TNF-alpha-mediated Trojan Horse for bacteria-based cancer therapy. Mol Ther (2022) 30(7):2522–36. doi: 10.1016/j.ymthe.2022.04.008 Ribas A Dummer R Puzanov I VanderWalde A Andtbacka RHI Michielin O . Oncolytic virotherapy promotes intratumoral T cell infiltration and improves anti-PD-1 immunotherapy. Cell (2017) 170(6):1109–19 e10. doi: 10.1016/j.cell.2017.08.027 Chesney JA Ribas A Long GV Kirkwood JM Dummer R Puzanov I . Randomized, double-blind, placebo-controlled, global phase III trial of talimogene laherparepvec combined with pembrolizumab for advanced melanoma. J Clin Oncol (2023) 41(3):528–40. doi: 10.1200/JCO.22.00343 Syyam A Nawaz A Ijaz A Sajjad U Fazil A Irfan S . Adenovirus vector system: construction, history and therapeutic applications. Biotechniques (2022) 73(6):297–305. doi: 10.2144/btn-2022-0051 Janku F Zhang HH Pezeshki A Goel S Murthy R Wang-Gillam A . Intratumoral injection of Clostridium novyi-NT spores in patients with treatment-refractory advanced solid tumors. Clin Cancer Res (2021) 27(1):96–106. doi: 10.1158/1078-0432.CCR-20-2065 Boesch M Horvath L Baty F Pircher A Wolf D Spahn S . Compartmentalization of the host microbiome: how tumor microbiota shapes checkpoint immunotherapy outcome and offers therapeutic prospects. J Immunother Cancer (2022) 10(11):e005401. doi: 10.1136/jitc-2022-005401 de Andres J Jimenez E Chico-Calero I Fresno M Fernandez L Rodriguez JM . Physiological translocation of lactic acid bacteria during pregnancy contributes to the composition of the milk microbiota in mice. Nutrients (2017) 10(1):14. doi: 10.3390/nu10010014 Cronin M Morrissey D Rajendran S El Mashad SM van Sinderen D O'Sullivan GC . Orally administered bifidobacteria as vehicles for delivery of agents to systemic tumors. Mol Ther (2010) 18(7):1397–407. doi: 10.1038/mt.2010.59 Walker SP Tangney M Claesson MJ . Sequence-based characterization of intratumoral bacteria-a guide to best practice. Front Oncol (2020) 10:179. doi: 10.3389/fonc.2020.00179 Hermida LC Gertz EM Ruppin E . Predicting cancer prognosis and drug response from the tumor microbiome. Nat Commun (2022) 13(1):2896. doi: 10.1038/s41467-022-30512-3 Zhou X Kandalai S Hossain F Zheng Q . Tumor microbiome metabolism: A game changer in cancer development and therapy. Front Oncol (2022) 12:933407. doi: 10.3389/fonc.2022.933407 Vautier AH . Vues générales sur la maladie cancéreuse. (Cand. Arsène-Hippolyte Vautier.). Méd.–Paris (1813). Dang LH Bettegowda C Huso DL Kinzler KW Vogelstein B . Combination bacteriolytic therapy for the treatment of experimental tumors. Proc Natl Acad Sci U.S.A. (2001) 98(26):15155–60. doi: 10.1073/pnas.251543698 Roberts NJ Zhang L Janku F Collins A Bai R-Y Staedtke V . Intratumoral injection of Clostridium novyi-NT spores induces antitumor responses. Sci Trans Med (2014) 6:249ra111–249ra111. doi: 10.1126/scitranslmed.3008982 Pangilinan CR Wu LH Lee CH . Salmonella impacts tumor-induced macrophage polarization, and inhibits SNAI1-mediated metastasis in melanoma. Cancers (Basel) (2021) 13(12):2894. doi: 10.3390/cancers13122894 Aganja RP Sivasankar C Hewawaduge C Lee JH . Safety assessment of compliant, highly invasive, lipid A-altered, O-antigen-defected Salmonella strains as prospective vaccine delivery systems. Vet Res (2022) 53(1):76. doi: 10.1186/s13567-022-01096-z Curtiss R 3rd Xin W Li Y Kong W Wanda SY Gunn B . New technologies in using recombinant attenuated Salmonella vaccine vectors. Crit Rev Immunol (2010) 30(3):255–70. doi: 10.1615/CritRevImmunol.v30.i3.30 Kong W . Development of antiviral vaccine utilizing self-destructing Salmonella for antigen and DNA vaccine delivery. Methods Mol Biol (2021) 2225:39–61. doi: 10.1007/978-1-0716-1012-1_3 Aliahmad P Miyake-Stoner SJ Geall AJ Wang NS . Next generation self-replicating RNA vectors for vaccines and immunotherapies. Cancer Gene Ther (2022) 30(6):785–93. doi: 10.1038/s41417-022-00435-8 Zheng X Wu Y Bi J Huang Y Cheng Y Li Y . The use of supercytokines, immunocytokines, engager cytokines, and other synthetic cytokines in immunotherapy. Cell Mol Immunol (2022) 19(2):192–209. doi: 10.1038/s41423-021-00786-6 Whitaker WR Shepherd ES Sonnenburg JL . Tunable expression tools enable single-cell strain distinction in the gut microbiome. Cell (2017) 169(3):538–46. doi: 10.1016/j.cell.2017.03.041 Cetnar DP Salis HM . Systematic quantification of sequence and structural determinants controlling mRNA stability in bacterial operons. ACS Synth Biol (2021) 10(2):318–32. doi: 10.1021/acssynbio.0c00471 Brockmeier U Caspers M Freudl R Jockwer A Noll T Eggert T . Systematic screening of all signal peptides from Bacillus subtilis: A powerful strategy in optimizing heterologous protein secretion in Gram-positive bacteria. J Mol Biol (2006) 362(3):393–402. doi: 10.1016/j.jmb.2006.07.034 Grasso S Dabene V Hendriks M Zwartjens P Pellaux R Held M . Signal peptide efficiency: from high-throughput data to prediction and explanation. ACS Synth Biol (2023) 12(2):390–404. doi: 10.1021/acssynbio.2c00328 Kruger A Welsch N Durwald A Brundiek H Wardenga R Piascheck H . A host-vector toolbox for improved secretory protein overproduction in Bacillus subtilis. Appl Microbiol Biotechnol (2022) 106(13-16):5137–51. doi: 10.1007/s00253-022-12062-2 Volkenborn K Kuschmierz L Benz N Lenz P Knapp A Jaeger KE . The length of ribosomal binding site spacer sequence controls the production yield for intracellular and secreted proteins by Bacillus subtilis. Microb Cell Fact (2020) 19(1):154. doi: 10.1186/s12934-020-01404-2 Yang H Qu J Zou W Shen W Chen X . An overview and future prospects of recombinant protein production in Bacillus subtilis. Appl Microbiol Biotechnol (2021) 105(18):6607–26. doi: 10.1007/s00253-021-11533-2 Yao D Su L Li N Wu J . Enhanced extracellular expression of Bacillus stearothermophilus alpha-amylase in Bacillus subtilis through signal peptide optimization, chaperone overexpression and alpha-amylase mutant selection. Microb Cell Fact (2019) 18(1):69. Abe K Nomura N Suzuki S . Biofilms: hot spots of horizontal gene transfer (HGT) in aquatic environments, with a focus on a new HGT mechanism. FEMS Microbiol Ecol (2020) 96(5):fiaa031. doi: 10.1093/femsec/fiaa031 Liu L Chen X Skogerbo G Zhang P Chen R He S . The human microbiome: a hot spot of microbial horizontal gene transfer. Genomics (2012) 100(5):265–70. doi: 10.1016/j.ygeno.2012.07.012 Thomas CM Nielsen KM . Mechanisms of, and barriers to, horizontal gene transfer between bacteria. Nat Rev Microbiol (2005) 3(9):711–21. doi: 10.1038/nrmicro1234 Deatherage DE Leon D Rodriguez AE Omar SK Barrick JE . Directed evolution of Escherichia coli with lower-than-natural plasmid mutation rates. Nucleic Acids Res (2018) 46(17):9236–50. doi: 10.1093/nar/gky751 Kanie S Horibata K Kawano M Isogawa A Sakai A Matsuo N . Roles of RecA protein in spontaneous mutagenesis in Escherichia coli. Genes Genet Syst (2007) 82(2):99–108. doi: 10.1266/ggs.82.99 Sakofsky CJ Grogan DW . Endogenous mutagenesis in recombinant sulfolobus plasmids. J Bacteriol (2013) 195(12):2776–85. doi: 10.1128/JB.00223-13 Canadas IC Groothuis D Zygouropoulou M Rodrigues R Minton NP . RiboCas: A universal CRISPR-based editing tool for Clostridium. ACS Synth Biol (2019) 8(6):1379–90. doi: 10.1021/acssynbio.9b00075 Wilding-Steele T Ramette Q Jacottin P Soucaille P . Improved CRISPR/Cas9 tools for the rapid metabolic engineering of Clostridium acetobutylicum . Int J Mol Sci (2021) 22(7):3704. doi: 10.3390/ijms22073704 Seys FM Rowe P Bolt EL Humphreys CM Minton NP . A gold standard, CRISPR/Cas9-based complementation strategy reliant on 24 nucleotide bookmark sequences. Genes (Basel) (2020) 11(4):458. doi: 10.3390/genes11040458 Zhang J Hong W Zong W Wang P Wang Y . Markerless genome editing in Clostridium beijerinckii using the CRISPR-Cpf1 system. J Biotechnol (2018) 284:27–30. doi: 10.1016/j.jbiotec.2018.07.040 Wasels F Jean-Marie J Collas F López-Contreras AM Lopes Ferreira N . A two-plasmid inducible CRISPR/Cas9 genome editing tool for Clostridium acetobutylicum . J Microbiological Methods (2017) 140:5–11. doi: 10.1016/j.mimet.2017.06.010 Kubiak AM Bailey TS Dubois LJ Theys J Lambin P . Efficient secretion of murine IL-2 from an attenuated strain of Clostridium sporogenes, a novel delivery vehicle for cancer immunotherapy. Front Microbiol (2021) 12:669488. doi: 10.3389/fmicb.2021.669488 Heap JT Theys J Ehsaan M Kubiak AM Dubois L Paesmans K . Spores of Clostridium engineered for clinical efficacy and safety cause regression and cure of tumors in vivo . Oncotarget (2014) 5:1761–9. doi: 10.18632/oncotarget.1761 Zhang Y Bailey TS Kubiak AM Lambin P Theys J . Heterologous gene regulation in clostridia: rationally designed gene regulation for industrial and medical applications. ACS Synth Biol (2022) 11(11):3817–28. doi: 10.1021/acssynbio.2c00401 Heap JT Pennington OJ Cartman ST Minton NP . A modular system for Clostridium shuttle plasmids. J Microbiological Methods (2009) 78:79–85. doi: 10.1016/j.mimet.2009.05.004 Concordet JP Haeussler M . CRISPOR: intuitive guide selection for CRISPR/Cas9 genome editing experiments and screens. Nucleic Acids Res (2018) 46(W1):W242–W5. doi: 10.1093/nar/gky354 Bailey TS Hittmeyer P Zhang Y Kubiak AM . Streamlined assembly of cloning and genome editing vectors for genus Clostridium. iScience (2023) 26(8):107484. doi: 10.1016/j.isci.2023.107484 Alva A Daniels GA Wong MKK Kaufman HL Morse MA McDermott DF . Contemporary experience with high-dose interleukin-2 therapy and impact on survival in patients with metastatic melanoma and metastatic renal cell carcinoma. Cancer Immunology Immunother (2016) 65(12):1533–44. doi: 10.1007/s00262-016-1910-x Soman G Yang X Jiang H Giardina S Vyas V Mitra G . MTS dye based colorimetric CTLL-2 cell proliferation assay for product release and stability monitoring of Interleukin-15: Assay qualification, standardization and statistical analysis. J Immunol Methods (2009) 348(1-2):83–94. doi: 10.1016/j.jim.2009.07.010 Knox RJ Boland MP Friedlos F Coles B Southan C Roberts JJ . The nitroreductase enzyme in Walker cells that activates 5-(aziridin-1-yl)-2,4-dinitrobenzamide (CB 1954) to 5-(aziridin-1-yl)-4-hydroxylamino-2-nitrobenzamide is a form of NAD(P)H dehydrogenase (quinone) (EC 1. 6.99.2). Biochem Pharmacol (1988) 37(24):4671–7. doi: 10.1016/0006-2952(88)90336-X Nagaraju S Davies NK Walker DJ Kopke M Simpson SD . Genome editing of Clostridium autoethanogenum using CRISPR/Cas9. Biotechnol Biofuels (2016) 9:219. doi: 10.1186/s13068-016-0638-3 Minton NP Ehsaan M Humphreys CM Little GT Baker J Henstra AM . A roadmap for gene system development in. Clostridium (2016) 41:104–12. Anaerobe. doi: 10.1016/j.anaerobe.2016.05.011 Woods C Humphreys CM Rodrigues RM Ingle P Rowe P Henstra AM . A novel conjugal donor strain for improved DNA transfer into Clostridium spp. Anaerobe (2019) 59:184–91. doi: 10.1016/j.anaerobe.2019.06.020 McAllister KN Sorg JA . CRISPR genome editing systems in the genus Clostridium: a timely advancement. J Bacteriol (2019) 201(16):e00219–19. doi: 10.1128/JB.00219-19 Payne R Glenn L Hoen H Richards B Smith JW Lufkin R . Durable responses and reversible toxicity of high-dose interleukin-2 treatment of melanoma and renal cancer in a Community Hospital Biotherapy Program. J ImmunoTherapy Cancer (2014) 14;2:13. doi: 10.1186/2051-1426-2-13 Amin A White RL . Interleukin-2 in renal cell carcinoma: A has-been or a still-viable option? J Kidney Cancer VHL (2014) 1:74. doi: 10.15586/jkcvhl.2014.18 McDermott DF Regan MM Clark JI Flaherty LE Weiss GR Logan TF . Randomized phase III trial of high-dose interleukin-2 versus subcutaneous interleukin-2 and interferon in patients with metastatic renal cell carcinoma. J Clin Oncol (2005) 23(1):133–41. doi: 10.1200/jco.2005.03.206 Buchbinder EI Dutcher JP Daniels GA Curti BD Patel SP Holtan SG . Therapy with high-dose Interleukin-2 (HD IL-2) in metastatic melanoma and renal cell carcinoma following PD1 or PDL1 inhibition. J ImmunoTherapy Cancer (2019) 7(1):49. doi: 10.1186/s40425-019-0522-3 Ahmadzadeh M Rosenberg SA . IL-2 administration increases CD4+CD25hi Foxp3 + regulatory T cells in cancer patients. Blood (2006) 107(6):2409–14. doi: 10.1182/blood-2005-06-2399 Kumar A Taghi Khani A Sanchez Ortiz A Swaminathan S . GM-CSF: a double-ddged sword in cancer immunotherapy. Front Immunol (2022) 13:901277. doi: 10.3389/fimmu.2022.901277 Almeida J . Phase 1/2 clinical trial of CP-506 (HAP) in monotherapy or with carboplatin or ICI. (2023). Available at: https://clinicaltrials.gov/study/NCT04954599 Saleski TE Chung MT Carruthers DN Khasbaatar A Kurabayashi K Lin XN . Optimized gene expression from bacterial chromosome by high-throughput integration and screening. Sci Adv (2021) 7(7):eabe1767. doi: 10.1126/sciadv.abe1767 Girbal L Mortier-Barriere I Raynaud F Rouanet C Croux C Soucaille P . Development of a sensitive gene expression reporter system and an inducible promoter-repressor system for Clostridium acetobutylicum. Appl Environ Microbiol (2003) 69(8):4985–8. doi: 10.1128/AEM.69.8.4985-4988.2003 Zhao L Ye B Zhang Q Cheng D Zhou C Cheng S . Construction of second generation protease-deficient hosts of Bacillus subtilis for secretion of foreign proteins. Biotechnol Bioeng (2019) 116(8):2052–60. doi: 10.1002/bit.26992 Kubiak AM Minton NP . The potential of clostridial spores as therapeutic delivery vehicles in tumour therapy. Res Microbiol (2015) 166:244–54. doi: 10.1016/j.resmic.2014.12.006 Barbe S Van Mellaert L Anne J . The use of clostridial spores for cancer treatment. J Appl Microbiol (2006) 101(3):571–8. doi: 10.1111/j.1365-2672.2006.02886.x Wu XC Lee W Tran L Wong SL . Engineering a Bacillus subtilis expression-secretion system with a strain deficient in six extracellular proteases. J Bacteriol (1991) 173(16):4952–8. doi: 10.1128/jb.173.16.4952-4958.1991 Aguilar Suarez R Stulke J van Dijl JM . Less is more: toward a genome-reduced Bacillus cell factory for "difficult proteins". ACS Synth Biol (2019) 8(1):99–108. doi: 10.1021/acssynbio.8b00342 van de Guchte M Kodde J van der Vossen JM Kok J Venema G . Heterologous gene expression in Lactococcus lactis subsp. lactis: synthesis, secretion, and processing of the Bacillus subtilis neutral protease. Appl Environ Microbiol (1990) 56(9):2606–11. doi: 10.1128/aem.56.9.2606-2611.1990