High efficiency recombineering in lactic acid bacteria

Bacterial strains, plasmids and media

Bacterial strains and plasmids used in this study are listed in Supplementary Table S1. Lactobacilli and their derivatives were cultured under anaerobic conditions at 37°C in deMan Rogosa Sharpe (MRS) medium (Difco, BD BioSciences). Lactococcus and its derivatives were cultured static at 30°C in M17-broth (Difco, BD BioSciences) that was supplemented with glucose to a final concentration of 0.5% (w/v). Electrocompetent cells of lactobacilli and Lactococcus were prepared as described before (28–30). If required, antibiotics were added as follows: 5 µg/ml erythromycin or chloramphenicol for lactobacilli and Lactococcus, 25 and 50 µg/ml rifampicin for lactobacilli and lactococci, respectively.

Reagents and enzymes

Restriction enzymes and T4 DNA ligases were purchased from New England Biolabs (NEB). PCR amplicons for cloning purposes were generated by KOD DNA hot start polymerase (Novagen), and PCR reactions for screening purposes were performed with Taq DNA polymerase (Denville Scientific). Pellet Paint Precipitant (Novagen) was used to precipitate DNA prior to restriction digestion, ligation and transformation. All oligonucleotides used in this study were purchased from Integrated DNA Technologies.

Construction of RecT expression vectors: pJP005 and pJP042

Oligonucleotides for cloning are listed in Supplementary Table S2. The RecT expression vectors (pJP005 and pJP042) were constructed as follows: pNZ8048 was digested with NcoI and KpnI, followed by Pellet Paint precipitation (Novagen). Genomic DNA of L. reuteri ATCC PTA 6475 (Wizard genomic DNA purification kit, Promega) was used as template for PCR amplification of the recT1 gene with oligonucleotide pair oJP022–oJP023, and the amplicon was digested and precipitated as described for pNZ8048. Purified DNA was quantified using Nanodrop-1000 (Thermo Scientific), and ligations were performed using a 1:1 (vector:insert) molar ratio using conditions specified by the manufacturer (NEB). The ligation mixture was transformed in L. lactis NZ9000, and Cm-resistant colonies were subjected to colony PCR with oligonucleotides that flank the multiple cloning site (oJP024–oJP025). Insertion was verified by restriction digest analysis and the integrity of the sequence was confirmed by sequence analysis. The resultant construct was named pJP005.

pSIP411 (a kind gift from Lars Axelsson, Nofima, Norway) is known to yield inducible and titratable gene expression in L. reuteri ATCC PTA 6475 (Lars Axelsson, personal communication) and was therefore used for this work. The backbone of pSIP411 was amplified with oligonucleotide pair oJP367–oJP368, followed by Pellet Paint precipitation, and digestion with NcoI and XbaI. Using oligonucleotide pair oJP369–oJP370 recT1 was amplified as described above, and the amplicon was purified followed by digestion with NcoI and XbaI. Ligation was performed as described for the construction of pJP005, and the ligation mix was transformed in L. lactis MG1363. Em-resistant colonies were screened by colony PCR for insertion of recT using oligonucleotide pair oJP415–oJP416. Integrity of the resultant clone, hereafter referred to as pJP042, was done in an identical manner as described for pJP005.

Preparation of cells for recombineering

All recombineering experiments with lactobacilli harboring pJP042 were performed with the following conditions: induction at OD₆₀₀ > 0.55 and <0.65 for 20 min with 10 ng/ml induction peptide (MAGNSSNFIHKIKQIFTHR; Peptide2.0 Inc), followed by 5 min incubation on ice and preparation of competent cells as described above. All electroporations described were performed with a BioRad Genepulser and 2-mm electroporation cuvettes (BioRad).

Lactococcus lactis NZ9000 harboring pJP005 was induced for 30 min with 0.1% (v/v) filter-sterilized nisin at OD₆₀₀ > 0.2 and < 0.3, followed by preparation of competent cells. Nisin was derived from an overnight culture of L. lactis NZ9700, and aliquots were stored at −20°C for up to 6 months.

Optimization of recombineering in L. reuteri

All oligonucleotides for recombineering are listed in Supplementary Table S2. Oligonucleotides for recombineering purposes were ordered at 100 nmol scale, desalted, and without any purifications. All recombineering oligonucleotides are complementary to the leading strand, and the non-complementary bases for all recombineering oligonucleotides are located centrally. With exception of the experiment to assess the effect of the olignucleotide concentration on the recombineering efficiency, we used 100 µg of oligonucleotide. Electrocompetent cells harboring pJP042 were prepared as described above, and for all experiments a volume of 5 µl (oligonucleotide or water) was electroporated. One microlitre of MRS was added, and bacteria were recovered for 2.5 h followed by serial dilution plating on MRS agar plates and MRS agar plates containing 25 µg/ml rifampicin. The number of rifampicin-resistant colonies is expressed per 10⁹ bacteria.

To assess the efficiency of recombineering in relation to a single or multiple mismatches we targeted the rpoB gene. Based on work in E. coli (31) we predicted the amino acid change H488R would yield a rifampicin-resistant phenotype. Oligonucleotides had a single mismatch (oJP133), or multiple mismatches (oJP311 and oJP577). No multiple mismatches could be designed which, when incorporated into a new DNA strand, would be silent mutations. As it is unknown whether an amino acid change adjacent to H488 affects the rifampicin-resistant phenotype, we identified experimentally which amino acid change would still yield a rifampicin-resistant phenotype, yielding oJP311 and oJP577.

Other optimization studies, including the effect of oligonucleotide concentration or multiple mutations in the wobble base in relation to the recombineering efficiency, were performed in an identical manner as described above.

Non-selected recombineering in L. reuteri

Competent cells of L. reuteri harboring pJP042 were transformed with 100 µg oligonucleotide oJP577 and after recovery serial dilutions were plated on MRS plates. One hundred colonies were replica plated on MRS agar plates and MRS agar plates containing 25 µg/ml rifampicin. Rifampicin-resistant colonies were subjected to PCR (using oligonucleotides oJP097–oJP098) followed by sequence analysis to confirm incorporation of the recombineering oligonucleotide.

For genes other than rpoB, a selection of genes was chosen that targeted different regions of the chromosome of L. reuteri. For each gene a 90-mer oligonucleotide was designed with multiple consecutive mismatches. With exception of oJP475 and oJP810, recombineering oligonucleotides were designed such that an in-frame stop codon and a restriction endonuclease recognition site (either EcoRI, BamHI or HindIII) were created upon incorporation. For all genes 100 µg recombineering oligonucleotide was transformed and after recovery serial dilutions were spread on MRS plates without antibiotic selection. From the pool of transformants approximately 300 colonies were subjected to screening. We applied two screening strategies. The first strategy involves identification of mutations by restriction digest when the mutation being incorporated generated a new restriction site. PCR amplification of a 1-kb fragment in which the target site is located centrally was performed. Amplicons were subsequently subjected to restriction digest analysis with the appropriate enzyme. A partial or completely digested amplicon was indicative of incorporation of the recombineering oligonucleotide. The rationale for the partial digest is that one strand is targeted by oligonucleotide recombineering, i.e. the PCR product can be derived from mixed genotypes yielding a mix of amplicons. Alternatively, a fully digested amplicon can be expected if the mixed genotypes separated during recovery. The second method is a mismatch amplification mutation analysis–PCR (MAMA–PCR) (32,33), and relies on the use of an oligonucleotide of which the 3′-end is complementary to the recombinant sequence. Upon identification of a recombinant genotype by either method, cells were streaked on plate to separate the genotypes followed by screening to identify a pure genotype. The oligonucleotide sequences related to the targets listed in Figure 3b are available upon request.

Recombineering in other LAB

We isolated natural rifampicin-resistant mutants for Lactobacillus plantarum BAA-793, Lactobacillus gasseri ATCC 33323, and L. lactis NZ9000, and determined the sequence of rpoB by sequence analysis. For lactobacilli, 10⁹ cells were plated on MRS containing 25 µg/ml rifampicin whereas for Lactococcus 10⁹ cells were plated on GM17 harboring 50 µg/ml rifampicin. Five rifampicin-resistant colonies of each strain were subjected to PCR using oligonucleotide pair oJP451–oJP458, oJP497–oJP505 and oJP528–oJP536 for L. plantarum, L. gasseri and L. lactis, respectively, and the sequence of the amplicons was determined and analyzed. Mutations that were identified in the rpoB genes were a G-T transversion at position 1438 (H480Y) for L. plantarum, a C–T transition at position 1468 (H490Y) for L. gasseri, and a C-A transversion at position 1456 (H486N) for L. lactis. Next we designed recombineering oligonucleotides targeting the rpoB gene for L. plantarum (oJP492), L. gasseri (oJP579) and L. lactis (oJP547 and oJP563). For L. gasseri we designed an oligonucleotide that, when incorporated, yields the amino acid changes T489S and H490R, and we determined experimentally that these changes yield a rifampicin-resistant phenotype. Competent cells for recombineering purposes were prepared as described above, followed by transformation of 100 µg of the corresponding oligonucleotides. The number of rifampicin-resistant colonies were determined and expressed per 10⁹ viable cells.

For construction of a double mutant by recombineering in L. lactis NZ9000, competent cells expressing RecT were transformed with 50 µg oJP563 (to yield a rifampicin-resistant phenotype) and 50 µg oJP584 which targets malG, and transformants were recovered on GM17 agar plates harboring 50 µg/ml rifampicin. A total of approximately 300 rifampicin-resistant colonies were screened by MAMA-PCR using oligonucleotides oJP575, oJP576 and oJP598. Two percent of the screened colonies yielded two amplicons, indicative for co-incorporation of oJP584. A streak plate was made to separate the genotypes, and subsequently 10 colonies were screened with the above oligonucleotides to identify a pure genotype yielding strain RPRB4004.

Mutating the active site of Ddl

Active site residues of the Ddl protein of E. coli had previously been identified and in vitro analysis with recombinant protein showed that a single amino acid change (Y216F) could convert a dipeptide ligase to a depsipeptide ligase (34). To identify the corresponding active site of Ddl from L. reuteri ATCC PTA 6475 we first obtained the predicted 3D structure using the SWISS-MODEL server (35). Next we used Chimera software (36) to perform a structural overlay of Ddl proteins derived from L. reuteri and E. coli (37) to determine the predicted active site residue of Ddl L. reuteri that corresponds to the Y216 of E. coli, and this was predicted to be F258 for Ddl-L. reuteri. E. coli crystal structures (PDB ID: 1IOV for wild-type, and 1IOW for Ddl-Y126F) were obtained from the RCSB protein data bank.

An 80-mer oligonucleotide was designed (oJP810) which, upon incorporation, yields seven base changes which resulted in a single amino acid change (F258Y). Electrocompetent cells expressing RecT1 were transformed with 100 µg oJP810. We anticipated a slower growth phenotype for Ddl (F258Y) compared to the wild-type strain. Overgrowth by wild-type cells in a mixed population makes it difficult to identify of a mutant genotype by colony screening. We therefore recovered the cells for 4 h instead of 2.5 h to allow separation of the genotypes, and plated the cells on MRS agar plates. After three days incubation, seven small colonies were subjected to PCR with oligo pair oJP697–oJP698 to yield a 500-bp product (Supplementary Figure S1), followed by ApoI restriction digest analysis according to the manufacturer specified conditions. As incorporation of oJP810 alters a single ApoI recognition site in the amplicon it is expected that an undigested PCR product is indicative for a mutant genotype (Supplementary Figure S1). Digested amplicons which yielded fragments corresponding to 300 and 200 bp are indicative for a wild-type genotype. The resultant mutant strain was named RPRB3003.

Vancomycin susceptibility testing

The susceptibility to vancomycin was assessed with Etest strips (BioMérieux), and the minimum inhibitory concentration (MIC) was expressed as micrograms per millilitre. The range of vancomycin concentrations tested on the strip was 0.016–256 µg/ml. The assay was performed as previously described (38).

Whole genome sequence analysis

We found that introduction of pSIP411 or its derivatives in L. reuteri ATCC PTA 6475 introduces a null-mutation in locus HMPREF0536_1706 in a sub-population of the transformants (unpublished results, Laura Ortiz, JPVP and RB). Prior to whole genome sequence analysis strains were selected in which no mutation was present in this locus. After identification of the mutant genotype, strains were cured from pJP042 by growing the strains for one passage without antibiotic selection followed by plating on non-selective MRS agar plates. By patch plating we identified erythromycin-sensitive colonies that are considered cured from pJP042. All strains with exception of RPRB1696 were subjected by Illumina genome sequencing by the Research Technology Support Facility (RTSF) at Michigan State University. Genomic DNA from L. reuteri and L. lactis was isolated with the Wizard Genomic DNA kit (Promega). For whole genome sequencing, genomic DNA samples were further processed with the Illumina Genomic Sample Prep kit according to the manufacturer's recommendations. Paired-end 75 base runs were performed on the Illumina GAIIx sequencer(39). Strain RPRB1696 was subjected to whole genome sequencing using a Roche 454 GS Junior sequencer. Samples were prepared according to the manufacturer's recommendations. Reads were mapped to the closed L. reuteri F275 genome sequence yielding for the Illumina data on average 139 ± 14-fold depth and 99.986 ± 0.01% genome coverage. The 454 sequencing resulted in a 36-fold depth and 99.446% genome coverage. Fold depth and coverage percentages were determined using the utility genomeCoverageBed within the BEDtools suite(40).

SNP analysis was performed with the software packages ‘From Reads to Results (R2R)’(41), MAQ (Sourceforge) and BWA(42) using standard settings. Identified SNPs in each of the recombineered strains were confirmed by PCR amplification of the mutated region followed by sequence analysis. The presence of larger INDELs and inversions was assessed by MAQ, BWA and Pindel(43).

Establishing recombineering in L. reuteri ATCC PTA 6475

Although there are several genetic tools available for LAB (44–47) there is no system to generate targeted point mutations in the chromosomes of members of the phylum Firmicutes without the need for antibiotic selection. Therefore we wanted to establish ssDNA recombineering in L. reuteri that may serve as a model for other lactobacilli and lactococci. Using E. coli recT as a query we searched the draft genome sequence of L. reuteri ATCC PTA 6475 and identified two recT homologs, recT1 (locus tag HMPREF0536_0455) and recT2 (locus tag HMPREF0536_0582).

To assess the ability of the L. reuteri RecT proteins to support recombineering we placed recT1 under the control of an inducible promoter using the sakacin-based expression vector pSIP411(48) to yield plasmid pJP042. A simple method to quantitatively assess the efficacy of recombineering is to mutate a site in the RNA polymerase gene (rpoB) that results in an amino acid change yielding a rifampicin-resistant phenotype. We designed an 80-mer oligonucleotide (oJP133) targeting the lagging strand of DNA synthesis which has a single centrally located non-complementary base that results in a C•A mismatch (Figure 2a). Incorporation of oJP133 generates the amino acid change H488R in RpoB that confers rifampicin resistance in L. reuteri. After optimization of RecT1 expression conditions (data not shown), transformation of 1 µg oJP133 yielded approximately a 10-fold increase in rifampicin-resistant colonies compared to background levels (Figure 2b). Expression of RecT2 yielded similar recombineering efficiencies (data not shown), and we used RecT1 for all future recombineering experiments in L. reuteri.

Figure 2.

Establishing and optimizing recombineering in L. reuteri. (a) The dsDNA sequence of the targeted rpoB region of L. reuteri is shown aligned with amino acids specified by each codon. On the left the leading and lagging strand are indicated and below the oligonucleotides are listed which result in different mismatches, all resulting in a rifampicin-resistant phenotype. The corresponding amino acid changes are listed on the right. (b) Titrations of the amount of oligonucleotides oJP133 (Ο) and oJP577 (Δ) were performed with 1, 5, 25 and 100 µg of oligonucleotide. Rifampicin resistant colonies derived from the control transformation are represented as a bar graph. Data are expressed as rifampicin-resistant cfu per 10⁹ cells. Data shown are the averages of three independent experiments and error bars represent standard deviation. (c) Assessment of recombineering efficiencies with 100 µg oJP133, oJP311, oJP813 and oJP577 in the wild-type strain (black bars), and in the mismatch repair deficient strains MutS1⁻ (light grey bars) and MutL⁻ (dark grey bars). Data shown are the averages of three independent experiments and error bars represent standard deviation.

High concentration of ssDNA required for efficient recombineering

In E. coli recombineering levels are saturated when using ∼115 ng of a 70 base oligonucleotide(14,49). To assess whether recombineering in L. reuteri is saturated with 1 µg, we transformed oJP133 over a two-log range with a maximum of 100 µg oligonucleotide (Figure 2b). There is a linear correlation between the number of recombinants that were obtained and the amount of oligonucleotide transformed up to 100 µg (R² = 0.98; Figure 2b). At 100 µg of oligonucleotide we observed a 1000-fold increase in recombineering efficiency over background levels. No increase in recombinants could be obtained by using >100 µg of oligonucleotide as this significantly reduced cell viability and recombineering efficiency. Co-transformation of carrier DNA or the use of oligonucleotides with phosphorothiorate linkages to prevent host nuclease activity did not increase the recombineering efficiency in L. reuteri (data not shown).

Increased recombineering efficiency using mutations that avoid mismatch repair

To assess if mismatch repair can be avoided in L. reuteri in a similar way as described for E. coli (see introduction), we designed an oligonucleotide (oJP311) with two non-complementary bases, one of which yields a C•C mismatch that is, at least in E. coli, not efficiently recognized by the MMR (Figure 2a). oJP311 also showed a linear correlation (R² = 0.99) for the increase in rifampicin-resistant colonies compared to the oligonucleotide concentration that was transformed (Figure 2b). The presence of an adjacent non-complementary base that yields a C•C mismatch did increase the recombineering efficiency, yielding on average 3 ± 1-fold more recombinants when compared to oJP133 (Figure 2b).

To assess whether the resultant C•C mismatch completely avoids mismatch repair in L. reuteri we inactivated two key players in MMR, mutS and mutL (see ‘Materials and Methods’ section for details), followed by testing the recombineering efficiencies with 100 µg of the above oligonucleotides. The MMR is only completely avoided when an oligonucleotide yields a similar number of recombinants in MutS⁻ and MutL⁻ strains compared to the wild-type. The level of recombinants obtained with oJP133 and oJP311 is on average 75-fold and 20-fold, respectively, lower in the wild-type background compared to MutS⁻ or MutL⁻ strains, demonstrating that in L. reuteri mismatch repair is not completely avoided with a single C•C mismatch (Figure 2c). We also find that, unlike in E. coli, an oligonucleotide yielding multiple mismatches at wobble positions (oJP813, Figure 2a) flanking the mutation yielding H488R in RpoB only partially avoided mismatch repair (Figure 2c).

To investigate whether the MMR can be completely avoided with an oligonucleotide with multiple adjacent mismatches we designed an oligonucleotide (oJP577) with five consecutive mismatches (Figure 2a). oJP577 yielded over a 100-fold increase in recombineering levels compared to recombineering levels obtained with an oligonucleotide with a single non-complementary base yielding a C•A mismatch. Next we assessed the recombineering levels of oJP577 in the MMR-deficient strains, showing that five adjacent mismatches yields similar levels of recombinants, showing that it completely avoids the MMR in L. reuteri (Figure 2c). Notably this level of recombineering yielded a 100 000-fold increase in recombineering compared to background levels, which indicated we had sufficient recombineering frequency to generate mutations without the need for selection.

Isolation of mutants without phenotypic selection in L. reuteri

We use the term ‘non-selected recombineering’ for the identification of mutants within a pool of cells without antibiotic or other phenotypic selection. Non-selected targeted recombineering is a powerful time-efficient approach that allows, when avoiding the MMR, construction of strain variants that are comparable to strains obtained through directed evolution. To identify if we could perform non-selected recombineering in L. reuteri we transformed oJP577, the most efficient oligonucleotide in this study, and plated the transformations on medium without rifampicin selection. One hundred colonies were picked from antibiotic-free plates and replica-plated on plates containing rifampicin. We found 5% of the transformants had incorporated the mutation into their genome (Figure 3a). These results show that targeted non-selected mutations can be made in L. reuteri.

Figure 3.

Non-selected recombineering in L. reuteri. (a) One hundred colonies derived from transformation with 100 µg oJP577, which upon incorporation yields a rifampicin-resistant phenotype (RpoB H488R), were patched onto MRS agar without antibiotic selection (left—MRS), and onto MRS agar containing rifampicin (right—MRS + RIF). (b) Overview of a selection of genes mutated by non-selected recombineering in the L. reuteri ATCC PTA 6475 chromosome. The locations of the mutated genes on the L. reuteri ATCC PTA 6475 chromosome are indicated based on the closed genome of the closely related strain L. reuteri F275 (NCBI accession number NC_010609). The numbers displayed represent the last four digits of the locus tag (HMPREF0535_XXXX) for each of the genes targeted and the recombineering frequency shown in parentheses. The rpoB gene (DNA directed RNA polymerase; locus tag HMPREF0536_0828) and the ddl gene (d-Ala-d-Ala ligase, locus tag HMPREF0536_1572) are listed using their gene names based on experimental evidence. ND: not determined; ori: origin of replication.

To show that the non-selected recombineering is not dependent on the location of rpoB, we targeted ten additional genes in L. reuteri that are found at different locations in the chromosome (Figure 3b). For each of these genes basepair changes were made without selection with recombineering rates that ranged from 0.4% to 19%. This demonstrates that non-selected recombineering can be performed throughout the genome. Currently we are investigating the basis for variation in recombineering efficiency at different loci.

Extending recombineering to other LAB

Current technologies do not allow the construction of targeted point-mutations in a high-throughput manner in L. lactis, a species with industrial and medical relevance for its application in the dairy industry in fermentation processes and its application as a biotherapeutic delivery vehicle. As pJP042 is not functional in L. lactis MG1363 (data not shown), we switched to the well-established nisin-inducible expression system for the expression of RecT1 [for a review, see (50)]. We used L. lactis NZ9000 as the expression host for pJP005, a derivative of pNZ8048 in which the L. reuteri recT1 gene was cloned under the control of the nisin-inducible promoter P_nisA. We designed oligonucleotides with one and four non-complementary bases that will yield a rifampicin resistant phenotype (Figure 4a). In L. lactis recombineering is more efficient with an oligonucleotide (oJP563) with multiple adjacent non-complementary bases as 50-fold more recombinants were obtained compared to a single non-complementary base that yields a T•C mismatch (Figure 4b). As observed in L. reuteri we find in L. lactis that recombineering efficiencies obtained with oJP563 are high enough for non-selected recombineering, yielding a 1% efficiency (data not shown).

Figure 4.

Establishing recombineering in other LAB. (a) The dsDNA sequence of the targeted rpoB region of L. lactis is shown aligned with the encoded amino acids for each codon. On the left are the leading and lagging strand indicated and below oligonucleotides are listed yielding different mismatches, all resulting in a rifampicin-resistant phenotype. The corresponding amino acid changes are listed on the right. (b) Recombineering efficiencies obtained in L. lactis transformed with 100 µg oJP547 or oJP563 and the level of natural rifampicin resistant L. lactis colonies per 10⁹ cells (control). (c) Comparison of recombineering levels between L. reuteri (oJP577), L. lactis (oJP563), L. gasseri (oJP579) and L. plantarum (oJP492). All data shown are averages from three independent experiments and error bars represent standard deviation.

We also extended recombineering to two human isolates of Lactobacillus whose genomes are publicly available: L. plantarum BAA 793 (51) and L. gasseri ATCC 33323 (52). We used pJP042 in combination with 100 µg oligonucleotide that targets rpoB to yield the amino acid changes D480Y and T489S + H490Y in L. plantarum and L. gasseri, respectively. We show, without any further optimization in these strains, that recombineering in L. gasseri yields ∼10⁴ rifampicin-resistant colonies with an oligonucleotide containing five adjacent mismatches (Figure 4c). We were only able to design an oligonucleoitde with a single mismatch in L. plantarum (oJP492) that yields a rifampicin-resistant phenotype, yielding approximately 10³ rifampicin-resistant colonies. Although recombineering levels in both strains are lower compared to L. reuteri and L. lactis, we anticipate that optimization of electroporation conditions, testing RecT from different hosts, optimizing expression levels of RecT, combined with oligonucleotide design will further increase the level of recombineering in these strains.

ssDNA recombineering did not yield untargeted genome changes in L. reuteri and L. lactis

A major concern of recombineering, which so far has not been adequately addressed in any organism, is the specificity of ssDNA recombineering and the possibility of obtaining additional untargeted mutations and genome rearrangements. In E. coli prolonged expression of the λ-Red genes yielded a hypermutagenic phenotype based on the increase of natural rifampicin-resistant colonies (53). To assess whether our approach is mutagenic we sequenced the complete genome of three L. reuteri strains: the wild-type control strain, a strain in which RecT was expressed from pJP042 but not exposed to ssDNA (RPRB0000), and a strain which has undergone a single recombineering event using 100 µg of oligonucleotide o1696 that targets a putative 6-aminohexanoate-cyclic-dimer hydrolase (locus tag HMPREF0536_1696; RPRB1696). In RPRB0000, a single base insertion was identified in the gene encoding dextransucrase (locus tag HMPREF0536_1850) resulting in a frame-shift with a consequent stop codon downstream. In RPRB1696 we confirmed the correct mutation in locus tag HMPREF0536_1696 specified by o1696. In addition we identified a single non-targeted point-mutation in a gene encoding a putative transcription regulation protein (locus tag HMPREF0536_1614), a G–T transversion in the third base of the start codon.

We also sequenced the genome of a L. lactis strain that has undergone a double recombineering event in which rpoB (locus tag LLNZ_10235) and malG (locus tag LLNZ_03840) were mutated (RPRB4004) simultaneously. Compared to the control strain (RPRB4000) no point mutations or other genome alterations were identified other than the mutations introduced by the oligonucleotides. Taken together these data indicate ssDNA recombineering is specific and not hypermutagenic in the sequenced L. reuteri and L. lactis isolates that have undergone a recombineering event.

Using recombineering to create a vancomycin-sensitive L. reuteri strain

To emphasize the utility of recombineering we converted L. reuteri from vancomycin resistant to vancomycin sensitive by mutating the active site of Ddl directly in the chromosome. Vancomycin inhibits the peptidoglycan biosynthesis by binding to muramylpentapeptides that terminate in d-alanyl-d-alanine (d-Ala-d-Ala) and blocks the addition of these precursors into the peptidoglycan backbone (54,55). Many lactobacilli incorporate d-alanyl-d-lactate (d-Ala-d-Lac) into their muramylpentapeptides, and this is the basis for their vancomycin resistance as peptides that terminate in d-Ala-d-Lac are bound 1000-fold less efficiently by vancomycin compared to peptides terminating in d-Ala-d-Ala (56,57). The enzyme responsible for the addition of d-Ala or d-lac is d-Ala-d-Ala ligase (Ddl). In E. coli Ddl is a dipeptide ligase as it catalyzes the formation of d-Ala-d-Ala whereas in L. mesenteorides, and probably in many lactobacilli, Ddl is a depsipeptide ligase as it catalyzes the formation of d-Ala-d-Lac. Superimposing the structures of Ddl-E. coli and Ddl-L. reuteri show that both proteins are structurally related (Figure 5a). In vitro analyses in E. coli and L. mesenteroides have previously identified active site residues in Ddl that upon mutation could switch from a dipeptide ligase to a depsipeptide ligase and vice versa(34,58). In Ddl-E. coli the residues Y216–E15–S150 form a hydrogen-bonding triad (Figure 5b) that is linked with the specificity of the enzyme although the exact mechanism for the shift in specificity is unknown (34). By replacing the tyrosine at position 216 with a phenylalanine (Y216F) in Ddl-E. coli there is no hydrogen bond formed between F216–E15, and as a consequence Ddl-Y216F has gained depsipeptide activity (Figure 5b) (34). To translate this to L. reuteri, we identified the predicted active site residues in Ddl- L. reuteri by superimposing Ddl-L. reuteri with Ddl-E. coli (Fig 5c). The L. reuteri triad is predicted to be F258–E17–S186, and mutating phenylalanine at position 258 to a tyrosine (F258Y) would, based on the E. coli model, yield a hydrogen-bond between Y258-E17 and change subsequently the enzyme activity from a depsipeptide ligase to a dipeptide ligase. As a result the cells render sensitivity to vancomycin.

Figure 5.

Converting vancomycin-resistant L. reuteri to vancomycin-sensitive by a single amino acid change. (a) The determined 3D structure of the d-ala-d-ala ligase (Ddl) protein of E.coli (PDB ID 1IOV; blue) and the modeled protein structure of Ddl of L. reuteri (green) were superimposed to locate the putative active site residues in Ddl of L. reuteri which, in E. coli, form a hydrogen-bonded network (Y216–E15–S150; boxed region). (b) Superimposing the determined structures of the Ddl proteins of the E. coli wild-type (blue) and a mutant (PDB ID 1IOW; grey) show that replacing a tyrosine at position 216 with a phenylalanine (Y216F) does not form a hydrogen-bond between F216–E15, and changes the enzymatic activity from a dipeptide ligase to a depsipeptide ligase. (c) We predicted the active site residues of Ddl-L. reuteri by superimposing the wild-type E. coli Ddl protein (blue) with the predicted protein structure of L. reuteri (green), and residues F258–E17–S186 in Ddl-L. reuteri correspond to the E. coli triad Y216–E15–S15, respectively. According to the E. coli model changing the phenylalanine at position 258 to a tyrosine (F258Y) in Ddl-L. reuteri would establish a hydrogen-bond between F258–E17 and may therefore yield dipeptide ligase activity which subsequently would result in a vancomycin-sensitive phenotype. (d) The dsDNA sequence of the ddl region of L. reuteri is shown aligned with the encoded amino acids for each triplet. Italic sequence represents the ApoI restriction site that is mutated by incorporation of oJP810. On the left are the leading and lagging strand indicated, and below is oJP810 showing the different mismatches with the resulting amino acid change (F258Y) on the right. (e) Susceptibility of L. reuteri wild-type (left) and the mutant RPRB3003 (right) to vancomycin using an Etest assay. The MIC is determined by identifying where bacterial growth intersects the Etest strip.

We designed an oligonucleotide that mutates two bases in the F258 codon of ddl along with bases in wobble positions in flanking codons to generate the single amino acid change F258Y (Figure 5d). The base changes disrupt an ApoI restriction endonuclease site and therefore we screened colonies by PCR amplification of the target region followed by ApoI restriction digest analysis. An undigested amplicon was obtained (Supplementary Figure S1) and we confirmed by sequence analysis that the mutation was successfully incorporated into the ddl gene. The resulting strain RPRB3003 showed increased sensitivity to vancomycin (MIC = 1.5 µg/ml) compared to the wild-type (MIC >256 µg/ml; Figure 5e). This demonstrates that non-selected ssDNA recombineering may be applied to strains used as biotherapeutic delivery vehicles to remove undesirable characteristics, such as unwanted antibiotic resistance.