Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Open Access

Peer-reviewed

Research Article

Plasmid Complement of Lactococcus lactis NCDO712 Reveals a Novel Pilus Gene Cluster

  • Mariya Tarazanova,

    Affiliations NIZO food research B.V., Ede, The Netherlands, TI Food and Nutrition, Wageningen, The Netherlands, Molecular Genetics, University of Groningen, Groningen, The Netherlands

    ⨯
  • Marke Beerthuyzen,

    Affiliations NIZO food research B.V., Ede, The Netherlands, TI Food and Nutrition, Wageningen, The Netherlands

    ⨯
  • Roland Siezen,

    Affiliations TI Food and Nutrition, Wageningen, The Netherlands, Centre for Molecular and Biomolecular Informatics, Radboud UMC, Nijmegen, The Netherlands, Microbial Bioinformatics, Ede, The Netherlands

    ⨯
  • Marcela M. Fernandez-Gutierrez,

    Affiliations TI Food and Nutrition, Wageningen, The Netherlands, Host-Microbe Interactomics Group, Animal Sciences, Wageningen University, Wageningen, The Netherlands

    ⨯
  • Anne de Jong,

    Affiliations TI Food and Nutrition, Wageningen, The Netherlands, Molecular Genetics, University of Groningen, Groningen, The Netherlands

    ⨯
  • Sjoerd van der Meulen,

    Affiliations TI Food and Nutrition, Wageningen, The Netherlands, Molecular Genetics, University of Groningen, Groningen, The Netherlands

    ⨯
  • Jan Kok,

    Affiliations TI Food and Nutrition, Wageningen, The Netherlands, Molecular Genetics, University of Groningen, Groningen, The Netherlands

    ⨯
  • Herwig Bachmann

    Herwig.Bachmann@nizo.com

    Affiliations NIZO food research B.V., Ede, The Netherlands, TI Food and Nutrition, Wageningen, The Netherlands

    ⨯

Abstract

Lactococcus lactis MG1363 is an important gram-positive model organism. It is a plasmid-free and phage-cured derivative of strain NCDO712. Plasmid-cured strains facilitate studies on molecular biological aspects, but many properties which make L. lactis an important organism in the dairy industry are plasmid encoded. We sequenced the total DNA of strain NCDO712 and, contrary to earlier reports, revealed that the strain carries 6 rather than 5 plasmids. A new 50-kb plasmid, designated pNZ712, encodes functional nisin immunity (nisCIP) and copper resistance (lcoRSABC). The copper resistance could be used as a marker for the conjugation of pNZ712 to L. lactis MG1614. A genome comparison with the plasmid cured daughter strain MG1363 showed that the number of single nucleotide polymorphisms that accumulated in the laboratory since the strains diverted more than 30 years ago is limited to 11 of which only 5 lead to amino acid changes. The 16-kb plasmid pSH74 was found to contain a novel 8-kb pilus gene cluster spaCB-spaA-srtC1-srtC2, which is predicted to encode a pilin tip protein SpaC, a pilus basal subunit SpaB, and a pilus backbone protein SpaA. The sortases SrtC1/SrtC2 are most likely involved in pilus polymerization while the chromosomally encoded SrtA could act to anchor the pilus to peptidoglycan in the cell wall. Overexpression of the pilus gene cluster from a multi-copy plasmid in L. lactis MG1363 resulted in cell chaining, aggregation, rapid sedimentation and increased conjugation efficiency of the cells. Electron microscopy showed that the over-expression of the pilus gene cluster leads to appendices on the cell surfaces. A deletion of the gene encoding the putative basal protein spaB, by truncating spaCB, led to more pilus-like structures on the cell surface, but cell aggregation and cell chaining were no longer observed. This is consistent with the prediction that spaB is involved in the anchoring of the pili to the cell.

Citation: Tarazanova M, Beerthuyzen M, Siezen R, Fernandez-Gutierrez MM, de Jong A, van der Meulen S, et al. (2016) Plasmid Complement of Lactococcus lactis NCDO712 Reveals a Novel Pilus Gene Cluster. PLoS ONE 11(12): e0167970. https://doi.org/10.1371/journal.pone.0167970

Editor: Akira Ishihama, Hosei Daigaku, JAPAN

Received: August 17, 2016; Accepted: November 24, 2016; Published: December 12, 2016

Copyright: © 2016 Tarazanova et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: All data are contained in the paper and the supporting information as well as in public databases (accession numbers given in paper).

Funding: This study was funded by the Top Institute Food & Nutrition (TIFN, Program FF001, Wageningen, The Netherlands). MB and HB are also employed by NIZO Food Research. TIFN provided support in the form of salaries for authors HB, MB, RS, AdJ, SvM, JK and HB. Neither NIZO food research nor TIFN had an additional role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: A patent pertaining to the use of strains expressing lactococcal pili was filed by TIFN. The affiliations to TIFN and NIZO food research do not alter our adherence to PLOS ONE policies on sharing data and materials.

Introduction

Lactococcus lactis is a gram-positive, non-pathogenic, non-spore forming lactic acid bacterium (LAB) that is often isolated from plant material or a dairy environment [1,2]. It is widely used in the dairy industry as a starter culture for the production of cheese, buttermilk and quark. Strains of L. lactis typically contain one to eight different plasmids [3,4] ranging from 1 kb to more than 100 kb in size [5,6]. Examples are L. lactis SK11, a phage-resistant dairy strain used in cheese making that carries 5 plasmids [7], and L. lactis IL594, the 7-plasmid-containing parent of the plasmid-free laboratory strain IL1403 [8]. The plasmids often specify traits of industrial importance such as the ability to grow on lactose, to utilize milk protein or stress resistance [9,10]. Other important plasmid-encoded functions include bacteriocin production [11,12] and resistance [13,14], metal ion resistance [15], antibiotic resistance [16,17] and bacteriophage resistance [18,19]. A recent publication described a CRISPR-Cas system that was encoded by a lactococcal plasmid, although it was concluded not to be functional [20]. In addition, several genes related to lactococcal surface properties are carried on plasmids [21–23], such as aggL, a gene responsible for cell auto-aggregation, or genes responsible for adhesion to mucus [24].

One of the most intensively studied L. lactis strains is MG1363, a plasmid-cured derivative of strain NCDO712 [25,26]. NCDO712 was originally isolated from a dairy starter culture and was described to harbor 5 plasmids with molecular sizes of approximately 33, 9, 5.2, 2.5, and 1.8 MDa [26]. During plasmid curing of strain NCDO712, derivatives carrying individual plasmids were obtained, allowing a targeted analysis of plasmid-encoded functions [26]. Important biotechnological properties of the strain, namely lactose utilization and proteolysis, were linked to the 33 MDa (55 kb) plasmid pLP712 [26,27].

This lactose/protease plasmid pLP712 [27] can be transferred to other lactococcal strains by conjugation [28]. Conjugation can occur through co-integrate formation between pLP712 and a genome-encoded sex factor (SF). The co-integrate is roughly double the size of pLP712 since the SF is 59.5 kb [29]. After conjugation approximately half of the transconjugants displayed an aggregating phenotype and transferred the lactose/protease plasmid with high-frequency [28,30–33]. The aggregating phenotype is only seen in transconjugants carrying the pLP712-SF co-integrate [29,34] and was linked to the cluA gene on the SF [34,35]. This gene encodes a surface protein involved in cell-to-cell contact and cell aggregation [34,35].

Some L. lactis strains also express proteinaceous surface appendages called pili [21,22,24]. Pili are known to have different functions in bacteria, including adhesion to surfaces (type I pili) or motility (type IV pili) [36–39]. Pilus biosynthesis genes can be encoded by the chromosomal DNA [22] or by plasmids [21] and they are described to be involved in cell aggregation [40], bacterial adherence to host cells [41–43] and attachment to environmental substrates/surfaces [44].

Two plasmids of L. lactis NCDO712 have been sequenced. Plasmid pLP712 (55395 bp) harbors the genes for lactose import and catabolism, the extracellular protease PrtP, and genes encoding extracellular proteins, transposases, and hypothetical proteins [27]. Plasmid pSH71, the smallest one (2062 bp) [45], is highly similar to pWV01 [45,46] and both of them are the basis of an array of broad host-range cloning and gene expression vectors [45,47–51].

Here we sequenced all plasmids of L. lactis NCDO712 and found that contrary to earlier reports, it contains not 5, but 6 plasmids [26]. The additional plasmid encodes functional nisin immunity and copper resistance genes. Additionally, we identified a novel pilus gene cluster on plasmid pSH74, which we showed to be functional by over-expression and phenotypic analyses.

Materials and Methods

Bacterial strains, growth conditions and medium

L. lactis subsp. cremoris NCDO712 [26] and its derivatives (Table 1) were grown at 30°C in M17 (Oxoid, Thermo Scientific, Hampshire, UK). The lactose-positive L. lactis strains NDO712 and its derivatives SH4109 and MG1299 were grown in M17 containing 1% lactose (LM17), all other strains were grown in M17 supplemented with 1% glucose (GM17). When required, erythromycin (Ery; 10 μg/ml), chloramphenicol (Cm; 5 μg/ml), rifampicin (Rif; 50 μg/ml), or streptomycin (Str; 100 μg/ml) was added.

thumbnail
Download:
Table 1. Strains and plasmids used in this study.

https://doi.org/10.1371/journal.pone.0167970.t001

DNA sequencing and sequence assembly

Total DNA of L. lactis NCDO712 was isolated using phenol-chloroform extraction as previously described [53] with the following modifications. Exponentially growing cells were harvested by centrifugation (10 min at 6240 g) after which the cell pellet was re-suspended in THMS buffer (30 mM Tris-HCL (pH 8), 3 mM magnesium chloride, 25% sucrose) containing lysozyme (2 mg/ml) and 50 μg/ml RNase and incubated for 1 h at 37°C. Subsequently, the cells were treated with SDS (at 1% final concentration) for 20 min at 65°C. Next, proteinase K (0.3 mg/ml) was added and incubation was continued for 10 min at 37°C. Total DNA was extracted from the lysate using several extractions with phenol/chloroform after which the DNA was precipitated with isopropanol, was dissolved in sterile water and stored at 4°C.

The purified total DNA was sheared to fragments of approximately 500 bp using the Covaris ultrasone device (KBioscience, LGC, Köln, Germany). The paired-end NEB NExtGen library preparation kit (New England Biolabs, Inc., MA, US) was used according to the manufacturer’s instructions to prepare NGS libraries. The libraries were sequenced on an Illumina HiSeq2000 (Illumina, Inc., San Diego, CA, USA), providing paired-end sequences of 101 bases. Velvet [54,55] was used in combination with VelvetOptimiser (http://bioinformatics.net.au/software.velvetoptimiser.shtml) to perform de novo paired-end assembly of the genome. All contigs that did not map to the published nucleotide sequence of the L. lactis MG1363 genome (Accession: NC_009004.1) were assumed to be plasmid fragments; these were first scaffolded by mapping onto known L. lactis plasmids in the NCBI database. Further scaffolding was supported by PacBio sequencing (BaseClear, Leiden, The Netherlands) on a 5-kb library of NCDO712 total DNA. Remaining gaps in the plasmid sequences were closed with dedicated PCR reactions followed by amplicon sequencing (BaseClear, Leiden, The Netherlands).

Initial automatic annotation of the plasmids was performed using the RAST annotation server [56]. Manual curation of plasmid-encoded features was performed with Artemis [57,58], followed by family, domain, motif and context analyses of encoded proteins using BlastP (NCBI) and Interpro (http://www.ebi.ac.uk/interpro/). IS elements and transposase genes were identified using IS Finder (https://www-is.biotoul.fr//). The DNA sequences of the assembled plasmids were used for a BlastN search (http://blast.ncbi.nlm.nih.gov/Blast.cgi) in the NCBI plasmid database containing complete plasmids. Determination of single nucleotide polymorphisms (SNPs) in the nucleotide sequence of the chromosome of strain NCDO712 relative to that of its derivative strain MG1363 was performed using the Breseq software package [59] and the GenBank file: NC_009004 in combination with corresponding next-generation sequencing data: SRA064225 as templates.

For non-synonymous SNPs the software SIFT [60] and the UniProt-TrEMBL database (http://www.uniprot.org/) were used to predict whether an amino acid substitution would affect protein function.

The nucleotide sequences of the pSH74 and pNZ712 plasmids were deposited in the NCBI GenBank database with accession numbers KX138410 and KX138409, respectively.

Determination of nisin and copper resistance

Overnight cultures of L. lactis were diluted in fresh GM17 medium to a final optical density at 600 nm (OD600) of 0.03. To measure nisin resistance, nisin from L. lactis (N6764-5G, Sigma-Aldrich, Steinheim, Germany) was added to the medium to different end concentrations (0–20 ng/ml). L. lactis NZ9700 (Table 1) was used as a control. The strains were grown in 10 ml sterile tubes for 7 h at 30°C. The OD600 was measured after 4 h and after 7 h using a UV/Visible Ultrospec 2000 spectrophotometer (Pharmacia Biotech, Cambridge, England).

To measure copper resistance, CuSO4 (0–4.8 mM end concentrations) was added to the growth medium. A 96-well microplate with the samples was incubated for 21 h at 30°C. The OD600 was measured every 15 min with a SpectraMax spectrophotometer (Molecular Devices, Wokingham, UK). The minimum inihibitory concentration (MIC) of copper was determined as the concentration of CuSO4 (g/L) that did not result in visible growth after 15 hours of incubation.

Pilin overexpression in L. lactis

The spaCBA-srtC1-srtC2 locus (designated as pil locus) including its 300-bp upstream region was amplified with KOD Hot Start Polymerase (Merck Millipore, Madison, Wisconsin, USA) using the pilinPstI forward and the pilinXhoI reverse primers (S1 Table). The purified PCR product was digested with PstI and XhoI and ligated to similarly digested pIL253 using T4 DNA Ligase (Invitrogen, Breda, The Netherlands). The ligation mixture was used to transform [61] electro-competent [62] MG1363 cells. Transformants harbouring the anticipated pIL253pil plasmid (Table 1) were identified performing colony PCR with primers pSH74_FW and pSH74_RV (S1 Table).

An internal deletion of 1451 bp in the spaCB gene was constructed by digestion of pIL253pil with AatII followed by re-ligation and introduction of the plasmid in L. lactis MG1363. The resulting plasmid was designated pIL253pilΔ1. Transformation of plasmids into other strains was performed with the protocol described above.

Cell aggregation

L. lactis cells from a 10-ml overnight culture were washed twice with 10 ml sterile 10 mM phosphate buffer (pH 6.8) and re-suspended in the same buffer. Cell sedimentation was observed visually and cell chaining was examined using bright field microscopy.

Scanning Electron Microscopy (SEM)

Bacterial cells were cultured for 1 day on GM17 agar plates. From plates with 50–100 colonies, small pieces of agar gel carrying less than 5 colonies were cut out and placed in a microscope sample holder. All further steps of cell fixation, washing, dehydration, staining, freeze-drying, electron microscopy, and image analysis were performed according to [63]. A FEI Magellan 400 FESEM electron microscope (Wageningen Electron Microscopy Centre, The Netherlands) was used for imaging.

Conjugation

Conjugation experiment were performed as described previously [64,65] with L. lactis MG1614 as the recipient strain. Transconjugants were selected on milk agar plates containing 0,004% bromocresol purple (Merck, Darmstadt, Germany), streptomycin (100 μg/ml) and rifampicin (50 μg/ml) with the donor strains NCDO712(pIL253pil) or MG1299(pIL253pil). LM17 plates supplemented with streptomycin (100 μg/ml), rifampicin (50 μg/ml) and 1.2 mM CuSO4 were used for the conjugative transfer of pNZ712, carrying the copper resistance genes, from L. lactis NCDO712 to L. lactis MG1614.

Results and Discussion

Chromosomal differences between L. lactis strains NCDO712 and MG1363

Sequencing of the total DNA of NCDO712 allowed detecting 11 Single Nucleotide Polymorphisms (SNPs) between the chromosomes of L. lactis NCDO712 and its plasmid-cured derivative MG1363 [66]. The latter was isolated in 1983 following multiple rounds of chemical- and protoplast-induced plasmid curing [26]. Among the 11 SNPs in the chromosome of L. lactis NCDO712, three are synonymous and three are in intergenic regions, while the other five lead to amino acid changes in proteins (S2 Table). The nucleotide sequencing data also suggests the occurrence of genome re-arrangements mainly due to mobile genetic elements, but their verification was outside the scope of this study.

Only one of the three SNPs in the intergenic regions is predicted to be in a promoter region, that of the mtlA gene encoding a putative mannitol-specific PTS system EIIBC component. Differential RNA sequencing has pinpointed the transcription start site (TSS) of mtlA at position 26465 [67]. The mutation at the position 26455 suggests that the -10 box is altered from an optimal TATAAT into TACAAT. Furthermore, three of the protein sequence-affecting SNPs in the genes encoding a hypothetical protein and two transposases were predicted not to affect protein functions. These predictions were made using SIFT, an algorithm that analyzes the effect of mutations based on the degree of conservation of amino acid residues [60,68]. Mutation in the gapB and tsf genes encoding the glyceraldehyde 3-phosphate dehydrogenase and elongation factor TS, respectively, were predicted to affect protein function. Whether these mutations are caused by genetic drift or whether they might confer a fitness advantage in a laboratory environment is unclear. However, the data indicate that the number of SNPs that accumulated since the isolation of MG1363 more than 30 years ago is limited.

MG1363 is also described to have lost a ~40 kb DNA fragment during its generation from NCDO712, which rendered MG1363 insensitive to UV induced lysis [69]. Based on this information we expected to find this additional ~40 kb fragment in the genome of NCDO712 when comparing the sequence to MG1363. Because we were not able to identify such a fragment we tried to induce lysis in NCDO712 by treatment with either mitomycin C or exposure to UV radiation. With neither approach we were able to induce lysis in NCDO712, suggesting that the variant we were working with did spontaneously loose this fragment. A NCDO712 variant without this 40 kb fragment is phenotypically the same as strain SH4109. However due to potential genome rearrangements in our NCDO712 isolate and the known high genome plasticity in this strain and its derivatives, it cannot be concluded that our isolate is identical to SH4109.

L. lactis NCDO712 harbors six plasmids

Assembly of all nucleotide sequence reads that did not map onto the chromosome of L. lactis MG1363 revealed that L. lactis NCDO712 carries a total of 6 rather than the previously described 5 plasmids [26]. Using L. lactis NCDO712 derivatives harboring single plasmid species [26], we linked the plasmids identified here to the earlier described plasmids of this strain (Table 2). The published plasmid sizes did not fully correspond to the respective sizes determined here. This could be caused by the limitations of size estimation based on agarose gel electrophoresis used previously [26]. However, we cannot exclude that plasmid rearrangements have occurred during strain propagation over the years that could account for (part of) the differences, although we did not observe such genetic changes in the 55.4-kb pLP712 plasmid (see below) [27]. The additional plasmid identified in this study, designated pNZ712, has a size similar to that of pLP712, which may explain why it escaped detection in 1983 [26].

thumbnail
Download:
Table 2. Comparison of NCDO712 plasmids with earlier studies.

https://doi.org/10.1371/journal.pone.0167970.t002

The copy number of the individual plasmids varied between 2 and 4, based on mean coverage number of chromosomal DNA and plasmid coverage (Table 2). The plasmid replication mode was determined using previously described criteria (2, 53). Rolling circle replication (RCR) was identified on the basis of the presence of Rep-family protein encoding genes and a double-stranded origin (dso) of replication, while a replication initiator protein-encoding repB gene and an origin of replication (ori) are indicative of theta-type plasmid replication [9]. These analyses indicate that pSH71 replicates through a rolling circle mechanism [9,45], while the other 5 plasmids replicate using a theta-type mechanism [9].

Plasmids pLP712 and pSH71 have been sequenced and described earlier [27,45]. The nucleotide sequence of the pLP712 is identical to the one determined here, except for a single nucleotide difference, whereas the pSH71 sequence differs by 6 nucleotides (Table 2).

To investigate the relationships of plasmids pSH72, pSH73, pNZ712, and pSH74 with other known plasmids we compared their nucleotide sequences with 1955 plasmid sequences in the NCBI database (database consulted on February 1, 2015). pSH72 (3597 bp) has the highest copy number, approximately 4 copies per cell, and only appears to encode the replication genes repB, repX, and repC. Except for 3 nucleotide differences pSH72 is identical to plasmid pND324 isolated from L. lactis subsp. lactis LL57-1 [70] (Table 2). The biological function of this plasmid is unclear. pSH73 is identical to pAG6, a plasmid isolated from Lactococcus lactis ssp. cremoris 712, which is most likely the same strain as NCDO712 or a derivative of it [29]. The only SNP (A → G, pAG6 → pSH73) is at nucleotide position 1143 of the hsdS gene encoding a type-I restriction/modification system specificity subunit. This SNP is silent and does not change the amino acid sequence. Plasmid pSH73 harbors, next to the replication genes repX and repB, also cadCA genes that are predicted to encode a cadmium resistance regulatory protein and a cadmium efflux ATPase (Fig 1 and S3 Table).

thumbnail
Download:
Fig 1. Plasmid maps of pSH73, pSH74, and pNZ712.

Plasmid pNZ712 includes genes encoding functional nisin immunity (nisCIP) and copper resistance (lcoRSABC). The 16-kb plasmid pSH74 contains a novel 8-kb pilus gene cluster spaCB-spaA-srtC1-srtC2. Plasmid pSH73 harbors repX, repB and cadCA genes. The latter genes were annotated as a cadmium resistance regulatory protein and a cadmium efflux ATPase.

https://doi.org/10.1371/journal.pone.0167970.g001

The two other plasmids pSH74 and pNZ712 have partial similarity (93–99% nucleotide identity) to known L. lactis plasmids. These high identities were observed in genes encoding several functions such as replication, transposase, resolvase, and in copper resistance-associated genes and in the nisCIP genes (Fig 1, S1 Fig and S4 Table). Detailed sequence analysis of pSH74 identified putative pilus biosynthesis genes, which we annotated as spaCB-spaA-srtC1-srtC2.

Plasmid pNZ712 specifies nisin immunity and copper resistance

Nisin is a lanthionine-containing antimicrobial peptide that binds to lipid II, disrupts the cytoplasmic membrane and causes death of susceptible bacterial cells [71]. The nisin operon nisABTCIPRKFEG carries, next to the nisin structural gene nisA, genes responsible for nisin modification, transport and precursor cleavage, genes involved in nisin operon regulation and genes specifying resistance to the bacteriocin [47,72]. Plasmid pNZ712 carries nisCIP, albeit nisC is only partially present. The nisC gene product NisC, in concerted action with NisB, is involved in posttranslational modification of the nisin precursor. The nisI gene specifies nisin immunity, while nisP encodes the serine protease involved in maturation of the nisin precursor.

To determine whether pNZ712 nisI is functional, L. lactis NCDO712 was grown in LM17 medium supplemented with 0 or 20 ng/ml nisin. L. lactis NZ9700 was used as a nisin resistant control. All strains reached a maximal OD600 of 2.96±0.63 after 7 h of growth without nisin. When they were grown in the presence of 20 ng/ml nisin, the positive control strain NZ9700 reached an OD600 of 1.8±0.06. L. lactis NCDO712 and SH4109, two strains carrying all 6 plasmids including pNZ712 (nisCIP), reached an OD600 of 0.63±0.01 and 0.67±0.09, respectively. The latter strains reached lower final optical densities than strain NZ9700, which has full immunity function as it carries nisFEG genes encoding an ABC transporter [73] that contributes to nisin immunity nisI in NCDO712. The plasmid-free L. lactis MG1363 and other derivatives carrying single plasmids from L. lactis NCDO712, but not pNZ712, reached an OD600 of only 0.06±0.046. From these results we conclude that the pNZ712-encoded nisI is functional and provides nisin resistance. The production of nisin has been suggested to give an advantage to L. lactis in a plant environment where nis genes often co-occur with sucrose utilization genes [74]. The nisin operon seems less prevalent in dairy isolates, indicating a limited benefit in this environment. The constitution of the starter culture from which L. lactis NCDO712 was isolated is unknown, but if this culture would have contained a nisin producer the maintenance of the nisCIP genes could have given an advantage.

Plasmid pNZ712 also carries a set of genes, lcoRSABC, which may be involved in copper resistance. Plasmid-encoded copper resistance has been described previously in Streptococcus (now Lactococcus) lactis [75,76], but has not been described for L. lactis NCDO712. The genes lcoRS are involved in transcription regulation of lcoABC, the products of which confer copper resistance by lowering the accumulation of copper inside the lactococcal cell [7]. Furthermore, pNZ712_23 encodes a putative copper-(or silver)-translocating P-type ATPase.

To examine whether pNZ712 specifies copper resistance, L. lactis NCDO712 and several of its derivatives were grown with or without CuSO4. L. lactis MG1363 and other NCDO712 derivatives did not grow when 0.8 mM or more CuSO4 was present in the medium. L. lactis NCDO712 was the only of the tested strains that was able to grow in the presence of up to 4 mM (1 g/L) CuSO4, indicating that the lcoRSABC gene cluster and/or pNZ712_23 on pNZ712 are functional. The minimum inhibitory concentration of CuSO4 was1 g/l (S2 Fig). We also showed that pNZ712 harboring the copper- and nisin-resistance genes can be transferred from NCDO712 to MG1614 via conjugation using copper as a selective marker. The conjugation frequencies were 8.1e-8 ± 2.6e-8 transconjugants per donor cell and using PCR we confirmed that only pNZ712 was transferred to the recipient strain.

Interestingly, a recently diverged derivative of L. lactis NCDO712, strain C2 [69,77], was reported to carry 5 plasmids. The resistance of this strain to copper and other metal ions was suggested to be encoded by the lactose plasmid pLM3001 (approximately 30 MDa), of which the nucleotide sequence has not been determined [75]. Our results, showing that copper resistance and lactose utilization genes reside on different plasmids in strain NCDO712, indicate that plasmid re-arrangements may have occurred in either or both of these strains [77].

A novel pilus gene cluster is present on pSH74

Plasmid pSH74 harbors an 8-kb gene cluster, which we annotated as spaCB-spaA-srtC1-srtC2 (Fig 2). The order of genes in pilus biosynthesis clusters differs among gram-positive bacteria [39,78]. The synteny on pSH74 resembles that of the spaC-spaB-spaA-srtC1 cluster of Lactobacillus rhamnosus GG most [42,43,79]. However, spaC and spaB seem to be fused in L. lactis NCDO712 (see below) and the operon-encoded proteins display only 30 to 45% amino acid sequence identity. An additional difference to L. rhamnosus GG is that the GG pilus cluster contains only one sortase gene while the pSH74 pilin gene cluster contains two adjacent srtC genes encoding SrtC enzymes of 413 and 392 amino acid residues, respectively. These show only 38% mutual amino acid identity. Two or even three consecutive srtC genes have previously been identified in the pili biosynthesis gene clusters in Streptococcus agalactiae, S. pneumoniae, and Clostridium diphteriae [80,81].

thumbnail
Download:
Fig 2. Schematic representation of the pilus gene cluster on pSH74 and the suggested functions of its constituting genes (not to scale).

The internal deletion of 1451 bases in spaCB is indicated in red.

https://doi.org/10.1371/journal.pone.0167970.g002

The GC content in the NCDO712 pil operon (35.3%) closely matches that of the NCDO712 chromosome (35.7%) suggesting that its acquisition is not a recent event or that it has been transferred from a microbial species with a similar GC content. The latter is consistent with the presence of a highly homologous pilus gene cluster in the genome of Leuconostoc citreum (> 90% protein sequence identity) with a GC content of 38.9%.

Previously two other pilus gene clusters were identified in L. lactis. Oxaran and co-workers [22] identified a chromosomally encoded pilus gene cluster in L. lactis IL1403, yhgD, yhgE, yhhA, yhhB. This chromosomal cluster, which does not show homology at the protein or synteny level to the pSH74 pilus gene cluster, appears to be present in all lactococcal genomes known to date, including those of L. lactis strains SK11, KF147, and MG1363 (llmg_1800–1803) [22] as well as NCDO712. Similarly, a plasmid-located 6.9-kb pilus gene cluster yhgE2-srtC2-yhhB2-ORF4 detected in L. lactis TIL448 [21] does not show significant homology or synteny to the pilus gene cluster identified here. A sequence comparison between the chromosomal pilus gene cluster in NCDO712 and the plasmid pSH74-located gene cluster identified here showed no synteny. Protein sequence analysis showed that SrtC1 and SrtC2 have a 35% (over 66% of the query) and 32% (over 73% of the query) identity with the chromosomally encoded SrtC (llmg_1801) respectively. Typically less than 30% amino acid identity is observed for the other proteins if the query sequence is at least 20% of the length of the protein. The presence of two pilus gene clusters in one strain suggests that they might have different functions and it raises the question whether L. lactis NCDO712 is able to produce the two types of pili simultaneously. Previous work has shown that L. pneumophila AA100 expressed both, long and short pili, but not at the same time [82] while concurrent expression of two types of pili has been reported for E. faecium [83]. L. lactis MG1363, which contains the chromosomal pil operon, only produces pili on its cell surface when it overexpresses the pil operon from pSH74, as is clear from the electron microscopy data and phenotypic characterization presented here. Whether conditions exist in which both types of pili are expressed simultaneously in L. lactis NCDO712 remains to be determined.

We hypothesize that the functions of the pilin gene products of L. lactis NCDO712 pSH74 are similar to those in Lb. rhamnosus GG [80] and other gram-positive bacteria [36,38,84–87]. Based on this hypothesis the SpaA protein of L. lactis NCDO712, SpaALl, is presumably the major pilus backbone subunit (Fig 2). It contains the typical C-terminal LPSTGGAG motif involved in sortase C-catalyzed transpeptidation and has the characteristic YPKN “pilin motif”. The YPKN motif is involved in intermolecular peptide bond formation between the carbonyl-group carbon of the threonine residue of the pilin subunit and the side-chain of the lysine in the pilin motif of the neighboring pilin subunit [41]. These bonds lead to the formation of membrane-associated covalently-linked dimers with a pilin motif that can interact with other pilin subunits, forming an elongated pilus fibre [41]. SpaALl also carries the YVLNETKAP “E box”, which has a structural role in pilus assembly (52). This was illustrated by the involvement of the “E box” of SpaA from Corynebacterium diphtheria in the attachment of SpaB to polymerized SpaA fibres [38].

Lb. rhamnosus GG SpaB (SpaBGG) is the basal pilin subunit and SpaCGG is the pilin tip protein. While spaC and spaB are separate genes in Lb. rhamnosus GG, they are fused in L. lactis NCDO712. The first 840 amino acid residues of L. lactis NCDO712 SpaCB (SpaCBLl) correspond to SpaCLl, while the remaining 260 residues resemble SpaBLl. The SpaC segment in SpaCBLl contains an LPSTGGSG motif which is a potential cleavage site for SrtC between the threonine (T) and glycine (G) residue [41], thereby splitting SpaCBLl into the respective proteins, SpaCLl and SpaBLl. The SpaBLl resembling part of SpaCBLl contains a C-terminal LPDTGE motif (Fig 2B) that is predicted to be targeted by SrtA and serve as a peptidoglycan anchoring sequence [79,80]. Hence, SpaBLl is most likely at the base of the lactococcal pilus.

The SpaCLl segment of SpaCBLl has an “E box” (YALTETKTP), which has a structural role in pilus assembly (14, 61, 78). Possibly, after cleavage of SpaCBLl by SrtC, the “E box” is used to link SpaCLl to SpaALl. The SpaCLl segment contains a von Willebrand type-A domain (vWFA). It has been speculated that a vWFA domain, which is also present in SpaC of L. rhamnosus GG [43], may have lectin-like binding properties and could bind to sugars with high carbohydrate specificity. SpaCLl also carries a collagen-binding domain and two collagen-binding surface-protein-Cna B-type domains, which might be involved in bacterial adhesion to surfaces (IPR008970; SSF49478) (http://supfam.org). Taken together, these observations suggest that, analogous to SpaCGG in the Lb. rhamnosus GG pilus, SpaCLl might fulfill the tip protein function.

Based on protein homology searches we identified similar pilus gene clusters in L. lactis subsp. cremoris CNCM I-1631 contig_071 (accession number: AGHX01000000; LLCRE1631_01806, LLCRE1631_01807, LLCRE1631_01808, LLCRE1631_01809) (97% identity), in L. lactis subsp. lactis 1AA59 contig_056 (accession number: AZQT01000047.1 and AZQT01000000; KHE76387.1, KHE76388.1, KHE76389.1, KHE76390.1) (100% identity), as well as in the Leuconostoc citreum genome (accession numbers: WP_048699698, WP_048699696, WP_048699695, WP_048699693) (>90% identity). In all these cases, the corresponding protein sequences are annotated as encoding hypothetical proteins. The spaCB genes are fused in these strains as well as in L. lactis NCDO712, suggesting that this might be of relevance for e.g., controlling the ratio of the two proteins in the pilus structure.

Pilus formation and attachment to peptidoglycan in Lb. rhamnosus GG is governed by the pilin-specific sortase SrtC1 and the house-keeping sortase SrtA. SrtC1 specifically targets a triple glycine motif LPxTGGxG at the N-terminal end of pilin proteins, and catalyzes their assembly into pili [84]. The chromosomally-encoded enzyme SrtA targets N-terminal LPxTGE or LPxTGD motifs [80] and covalently anchors extracellular proteins, including pilin proteins, to peptidoglycan in the cell wall. The SrtA recognition sequence does not necessarily prevent its recognition by a sortase other than SrtA [80]. Thus, SrtA and SrtC1 may each be able to recognize both LPxTGD/E and LPxTGGxG motifs used for pilin polymerization and anchoring [80]. However, the protein-structural features involved in motif recognition by sortases are currently unknown. When spaA and srtC1 of Lb. rhamnosus GG were co-expressed in L. lactis NZ9000, it was observed that SrtC1 recognized and polymerized the SpaA protein. Lb. rhamnosus GG SpaC has the same LPxTG motif as SpaA [80], therefore, SrtC1 possibly recognizes the pilus tip protein SpaC similar to SpaA, while SrtA only anchored SpaB to the peptidoglycan. Based on the similarity of the pSH74-located pilus cluster to that of Lb. rhamnosus GG, we speculate the gene functions to be similar, namely that one or both of the plasmid encoded C sortases in L. lactis NCDO712 are responsible for the assembly of pilin proteins in this strain, while the chromosomally encoded SrtA (llmg_1449) covalently anchors the pili to the peptidoglycan of the cell wall [79,80]. Data obtained from RNAseq performed on total RNA isolated from L. lactis NCDO712 growing under different stress conditions revealed that srtC1 and srtC2 are usually co-expressed, while there is no correlation in expression levels between the other genes of this gene cluster [67].

The pilus gene cluster on L. lactis NCDO712 pSH74 is functional

To examine whether the identified putative pilin genes are functional, the entire pilus gene cluster with its native promoter was cloned in the medium-copy number (45–85) plasmid pIL253 [52], resulting in pIL253pil. Initial attempts to clone pil downstream of the nisin-inducible nisA promoter or the constitutive purC promoter failed. The few clones that were obtained all carried the same internal deletion in the spaCB genes. Deliberate deletion of a similar internal 1.5 kb fragment resulted in pIL253pilΔ1 (Table 1, Fig 2). The deletion lead to a frame-shift mutation shortening the SpaCBLl protein by 744 amino acid residues and resulting in a 394-residue truncated SpaC protein. This protein still carries the vWFA domain but the YPKN pilin motif, the LPDTGE motif (Fig 2), the “E box” and LPS amino acids of LPSTGGSG motif were deleted. RNAseq data revealed that the native promoter includes a leader sequence of ~200 nt upstream of spaCB [67]. Although the presence of this leader seems important for the successful over-expression, its precise role is currently not known. Both plasmids pIL253pil and pIL253pilΔ1 were introduced in L. lactis strains MG1363 and IL1403. The observed pilus over-expression from the pIL253-based vector resulted in cell aggregation and rapid sedimentation of the cultures of both strains (Fig 3). The cells also grew in chains containing more than 10 cells, while the wild type strain MG1363 and control strain MG1363pIL253 did not form such chains (Fig 4). Cells expressing the truncated version of SpaCBLl displayed neither cell aggregation (Fig 3) nor cell chaining (Fig 4). As the anchoring of pili to the cell surface is expected to be essential for their functionality the effect of the spaCB truncation is consistent with the loss of cell aggregation and chaining.

thumbnail
Download:
Fig 3. Cell aggregation after pilus gene over-expression.

A—L. lactis NCDO712, B—empty vector control MG1363(pIL253), C—MG1363(pIL253pilΔ1), D—MG1363(pIL253pil). The images were taken 3 min after re-suspension of cells from an overnight culture.

https://doi.org/10.1371/journal.pone.0167970.g003

thumbnail
Download:
Fig 4. Pilin gene overexpression in L. lactis MG1363 leads to increased chain length.

A—MG1363(pIL253), B—MG1363(pIL253pil), C—MG1363(pIL253pilΔ1).

https://doi.org/10.1371/journal.pone.0167970.g004

Scanning electron microscopy (SEM) revealed that L. lactis NCDO712 (Fig 5A) cells have a relatively rough surface while those of L. lactis IL1403 were very smooth (compare A and B in Fig 5). The introduction of pIL253pil resulted in an increased roughness of the cell surface for both for IL1403 and MG1363. Furthermore, even though L. lactis NCDO712 harbors the pil operon on pSH74, pili were not visible on the surface of these bacteria. This indicates that expression of pili from the native plasmid, estimated to be present at 3 copies per cell (Table 2), is not high enough to allow their detection. However, expression of the same operon (with the native promoter) from plasmid pIL253, which is present in higher (45–85) copy numbers [52], is sufficient to detect pilus-like structures on cells of L. lactis IL1403(pIL253pil) (Figs 5C and 6E), MG1363(pIL253pil) (Figs 5E and 6C) and NCDO712(pIL253pil) (Fig 6A). Similar to what we found in NCDO712, detection of pili from the L. lactis IL1403 chromosomal pilus gene cluster yhgD, yhgE, yhhA, yhhB could only be achieved by their overexpression [22]. However, in L. lactis isolates from vegetables, such as KF282 and NCDO2118, and in the clinical isolates 2885–86 and 810–85, pili were detectable under standard growth conditions using negative staining and TEM analysis [22]. In another study wild-type L. lactis TIL448 was shown to carry a plasmid-encoded 6.9-kb pilus gene cluster, yhgE2-srtC2-yhhB2-ORF4, that lead to the formation of short (260–440 nm) and low density pili at the cell surface detectable with AFM [21]. Overexpression of the L. lactis TIL448 pilus gene cluster yhgE2-srtC2-yhhB2-ORF4 in L lactis MG1363 (TIL1293) led to thicker (2.3 nm), longer (1.2–2.5 μm) and higher density pili at the cell surface [21].

thumbnail
Download:
Fig 5. Scanning electron microscopy at 50.000 x magnification of L. lactis strains over-expressing spaCB-spaA-srtC1-srtC2.

A—NCDO712, B—IL1403, C—IL1403(pIL253pil), D—MG1363(pIL253), E—MG1363(pIL253pil), F—MG1363(pIL253pilΔ1). White arrows indicate pili in panel C. Black bars: 2 μm in all panels.

https://doi.org/10.1371/journal.pone.0167970.g005

thumbnail
Download:
Fig 6. Scanning electron microscopy of spaCB-spaA-srtC1-srtC2 over-expressing L. lactis strains.

A—NCDO712(pIL253pil), B—NCDO712, C—MG1363(pIL253pil), D—MG1363(pIL253), E—IL1403(pIL253pil), F- MG1363(pIL253pilΔ1). Black bars: 1 μm in all panels.

https://doi.org/10.1371/journal.pone.0167970.g006

Interestingly, MG1363(pIL253pilΔ1) cells expressing the truncated version of SpaCBLl (Figs 2 and 5F) were also decorated with pilus-like structures, but these appeared to be more disoriented (Fig 5, compare E and F), which would be consistent with SpaA being the shaft pilin. This disorientation may imply that the pili produced by MG1363(pIL253pilΔ1) are not appropriately attached to the peptidoglycan, which is consistent with the truncated form of SpaCBLl lacking the predicted SrtA motif for peptidoglycan anchoring. All in all, the results reveal that the spaCB-spaA-srtC1-srtC2 pilus gene cluster from pSH74 of L. lactis NCDO712 is functional and leads to the formation of pilus-like appendices on the surface of lactococcal cells.

The Lactobacillus rhamnosus GG pilus operon spaC-spaB-spaA-srtC1 is involved in attachment to human epithelial colorectal adenocarcinoma cells [43], and an intestinal mucus-binding capacity was attributed to the tip pilin subunit SpaC [42] via mucus-binding domain. Lb. rhamnosus GG also carries a pilus gene cluster, spaFED, on its chromosome, where SpaD is the pilus backbone protein, SpaE is at the base of the pilus structure, and SpaF is the minor pilin subunit that locates to the tip of the pilus. SpaF was also shown to be responsible for adhesion of pili-carrying cells to the intestinal mucus [42,88]. Furthermore, cloning of spaCGG, spaBGG, spaDGG, spaEGG, spaFGG in Escherichia coli and assessment of the adherence of these proteins to human intestinal mucus reveled that also the SpaB pilin subunit plays a role in binding to intestinal mucus, through electrostatic contacts [42]. In the L. lactis TIL448 pilus gene cluster yhgE2-srtC2-yhhB2-ORF4, the latter is the putative pilus tip protein. It contains a lectin-like domain (PF00139) predicted to have carbohydrate-binding properties, which could be involved in binding to plant cell walls [21]. Furthermore, over-expression of TIL448 yhgE2-srtC2-yhhB2-ORF4 in L. lactis MG1363 (TIL1293) was shown to increase attachment of the bacteria to Caco-2 cells [21]. No significant differences in attachment to Caco-2 and HT29 (a human colonic carcinoma cell line) cells were observed when stationary-phase cells of L. lactis with or without pIL253pil were tested (S5 Table). Similar to our results earlier studies also described that strain NCDO712 (synonym C2 or TIL256) did not adhere to Caco-2 cells [21].

In gram-negative bacteria the length of pili is approximately 1–2 μm, and the diameter is between 1 nm and 10 nm [36,41]. In gram-positive bacteria pili are more varied with reported lengths of 0.3–3 μm, and a diameter between 3–10 nm [41]. The Lb. rhamnosus GG spaC-spaB-spaA-srtC1 pili are up to 1 μm in length, with a diameter of 1 to10 nm [43,79]. The pili encoded by yhgD, yhgE, yhhA, yhhB in L. lactis IL1403 are up to 3 μm long and have a diameter of approximately 5 nm [22]. Pili on the surface of wild type L. lactis TIL448 specified by yhgE2-srtC2-yhhB2-ORF4 were short (up to 450 nm) and thin, but the pili of the overexpressing yhgE2-srtC2-yhhB2-ORF4 cluster strain L. lactis MG1363 (TIL1293) were longer with an average length of 2 μm and diameter of approximately 2.3 nm [21]. Interestingly, the pili overexpressed by MG1363(pIL253pil) are shorter and thicker than earlier reported lactococcal pili [21,22,24], having a length of 200–240 nm and a diameter between 18 and 20 nm. The appendages produced by MG1363(pIL253pilΔ1) were thinner and had an average diameter of approximately 14 nm and the length of 120–560 nm.

Conjugation efficiency

One of the roles of certain pili in gram-negative bacteria is to enable conjugation [89]. The overexpression of the pili identified here frequently displayed cell-cell contact (Fig 6A and 6C), which led us to examine whether the L. lactis NCDO712 pil operon might play a role in DNA transfer by conjugation. L. lactis NCDO712 and MG1363 carry a 50-kb sex factor in its chromosome that is involved in co-integrate formation with and subsequent conjugal transfer of the lactose/protease plasmid pLP712 [29,35]. This process can be readily quantified using the plasmid-free, lactose-deficient L. lactis MG1614 (Lac-, PrtP-, StrepR, RifR) (Table 1) as a recipient and selecting for lactose-fermenting colonies that are also resistant to rifampicin and streptomycin. The efficiency of transfer of pLP712 from strains NCDO712 and MG1299 was up to 22 and 16 fold increased respectively, when the strains carried pIL253pil and, thus overexpressed pili (S6 Table). We did not observe co-transformation of pIL253pil during the conjugation of pLP712. Although the observed increase in the frequency of conjugation is limited, the observations indicate that the pili contribute to the efficacy of exchange of DNA between lactococcal cells. Whether this effect is caused by pili-mediated cell clumping and, with that, increased cell-cell contact (Fig 6), or whether the pili are involved in actual DNA transfer requires further investigation.

Conclusion

Total DNA sequencing of L. lactis NCDO712 revealed functional, plasmid-encoded nisin immunity and copper resistance genes as well as a novel pilus gene cluster (spaCB-spaA-srtC1-srtC2). Based on bioinformatic analyses and experimental results we predict the pilin tip protein to be SpaC, the basal subunit to be formed by SpaB and SpaA to represent the backbone protein. While surface decoration of bacterial cells with pilus-like structures is widespread and functions such as DNA transfer [37], surface attachment and biofilm formation [22] or interaction with host cells[21] have been described [36], it is still unclear what the function(s) of pili in lactococci is. Pilus expression occurs under standard laboratory conditions in lactococcal strains isolated from plant material or from humans [22], but this is not the case for any of the investigated dairy isolates [1,7,25]. Pili of plant-associated Pseudomonas syringae are involved in the attachment of the bacteria to the plant surface [44]. It remains to be determined whether pili in dairy lactococci are evolutionary remnants from their proposed plant ancestry [90] or whether they fulfill a function in the dairy environment.

Supporting Information

S1 Fig. Plasmid identity to pNZ712.

Five out of 100 plasmids with high partial sequence identity to pNZ712 are shown. The similar regions include genes encoding copper resistance genes lcoRSAB, mobilization genes mobD, mobC and a gene involved in replication (repA) with 99% (yellow) and 93% (grey) identity (see also S4 Table). L. lactis pND306 of L. lactis subsp. lactis 1252D [91] encoding the copper resistance associated lcoRSABC genes were almost identical to a similar locus (6.6 kb) encoded by pNZ712 (only 3 SNPs in the 6.6 kb lcoRSABC locus), of which also a part (lcoRSA) is present on pSK11P from L. lactis subsp. cremoris SK11 [7].

https://doi.org/10.1371/journal.pone.0167970.s001

(TIF)

S2 Fig. Growth of L. lactis NCDO712 and its derivatives carrying one or two plasmids (see Table 1) in medium with or without 0.8 mM (0.2 g/L)–4.8 mM (1.2 g/L) CuSO4.

L. lactis NCDO712 is the only strain carrying pNZ712 with lcoRSABC coding for copper resistance genes. Each curve represents the average of 3 biological replications. Error bars show standard deviation.

https://doi.org/10.1371/journal.pone.0167970.s002

(TIF)

S1 Table. Oligonucleotides used in this study.

https://doi.org/10.1371/journal.pone.0167970.s003

(PDF)

S2 Table. Single nucleotide polymorphisms between the chromosomes of L. lactis NCDO712 and its plasmids-free derivative MG1363 [66].

https://doi.org/10.1371/journal.pone.0167970.s004

(PDF)

S3 Table. Main plasmid genes and features.

https://doi.org/10.1371/journal.pone.0167970.s005

(PDF)

S4 Table. The main alignment results of pNZ712 and pSH74 to known lactococcal plasmids.

https://doi.org/10.1371/journal.pone.0167970.s006

(PDF)

S5 Table. Microbial adhesion to HT29 and Caco cells given in percent cells recovered after incubation with the cell lines.

Standard deviation is given in parenthesis (n = 3).

https://doi.org/10.1371/journal.pone.0167970.s007

(PDF)

S6 Table. Effect of the pilus gene cluster on conjugation of the lactose-fermenting ability to L. lactis MG1614.

https://doi.org/10.1371/journal.pone.0167970.s008

(PDF)

Acknowledgments

We thank Marcel Giesbers (Electron Microscopy Center at Wageningen University) for his excellent technical assistance in pili visualization using SEM and Michiel Kleerebezem for critically reading the manuscript.

Author Contributions

  1. Conceptualization: MT JK HB.
  2. Formal analysis: MT AdJ RS HB.
  3. Investigation: MT MB MFG SvM.
  4. Supervision: JK HB.
  5. Writing – original draft: MT RS JK HB.
  6. Writing – review & editing: MT JK HB.

References

  1. 1. Bolotin A, Wincker P, Mauger S, Jaillon O, Malarme K, Weissenbach J, et al. The Complete Genome Sequence of the Lactic Acid Bacterium Lactococcus lactis ssp. lactis IL1403. Genome Res. 2001;11: 731–753. pmid:11337471
  2. 2. van Hylckama Vlieg JET, Rademaker JLW, Bachmann H, Molenaar D, Kelly WJ, Siezen RJ. Natural diversity and adaptive responses of Lactococcus lactis. Curr Opin Biotechnol. 2006;17: 183–90. pmid:16517150
  3. 3. Ainsworth S, Mahony J, van Sinderen D. The plasmid complement of Lactococcus lactis UC509.9 encodes multiple bacteriophage resistance systems. Appl Environ Microbiol. 2014;80: 4341–9. pmid:24814781
  4. 4. Fallico V, Ross RP, Fitzgerald GF, McAuliffe O. Novel conjugative plasmids from the natural isolate Lactococcus lactis subspecies cremoris DPC3758: a repository of genes for the potential improvement of dairy starters. J Dairy Sci. Elsevier; 2012;95: 3593–608.
  5. 5. Kojic M, Strahinic I, Topisirovic L. Proteinase PI and lactococcin A genes are located on the largest plasmid in Lactococcus lactis subsp. lactis bv. diacetylactis S50. Can J Microbiol. 2005;51: 305–14. pmid:15980892
  6. 6. Sing WD, Klaenhammer TR. Characterization of Restriction-Modification Plasmids from Lactococcus lactis ssp. cremoris and Their EffectsWhen Combined with pTR2030. J Dairy Sci. Elsevier; 1991;74: 1133–1144.
  7. 7. Siezen R, Renckens B. Complete sequences of four plasmids of Lactococcus lactis subsp. cremoris SK11 reveal extensive adaptation to the dairy environment. Appl Environ Microbiol. 2005;71: 8370–8383.
  8. 8. Górecki RK, Koryszewska-Bagińska A, Gołębiewski M, Żylińska J, Grynberg M, Bardowski JK. Adaptative Potential of the Lactococcus Lactis IL594 Strain Encoded in Its 7 Plasmids. Biswas I, editor. PLoS One. 2011;6: e22238. pmid:21789242
  9. 9. Ainsworth S, Stockdale S, Bottacini F, Mahony J, van Sinderen D. The Lactococcus lactis plasmidome: much learnt, yet still lots to discover. FEMS Microbiol Rev. 2014;38: 1066–88. pmid:24861818
  10. 10. Cui Y, Hu T, Qu X, Zhang L, Ding Z, Dong A. Plasmids from Food Lactic Acid Bacteria: Diversity, Similarity, and New Developments. Int J Mol Sci. 2015;16: 13172–13202. pmid:26068451
  11. 11. Coakley M, Fitzgerald G, Ros RP. Application and evaluation of the phage resistance- and bacteriocin-encoding plasmid pMRC01 for the improvement of dairy starter cultures. Appl Environ Microbiol. 1997;63: 1434–40. Available: http://www.ncbi.nlm.nih.gov/pubmed/9097441 pmid:9097441
  12. 12. van Belkum MJ, Hayema BJ, Jeeninga RE, Kok J, Venema G. Organization and nucleotide sequences of two lactococcal bacteriocin operons. Appl Environ Microbiol. 1991;57: 492–8. Available: http://www.ncbi.nlm.nih.gov/pubmed/1901707 pmid:1901707
  13. 13. Grohmann E, Muth G, Espinosa M. Conjugative Plasmid Transfer in Gram-Positive Bacteria. Microbiol Mol Biol Rev. 2003;67: 277–301. pmid:12794193
  14. 14. Campelo AB, Roces C, Mohedano ML, Lopez P, Rodriguez A, Martinez B. A bacteriocin gene cluster able to enhance plasmid maintenance in Lactococcus lactis. Microb Cell Fact. 2014;13: 77. pmid:24886591
  15. 15. Mills S, McAuliffe OE, Coffey A, Fitzgerald GF, Ross RP. Plasmids of lactococci-genetic accessories or genetic necessities? [Internet]. FEMS Microbiology Reviews. 2006. pp. 243–273.
  16. 16. Mantovani HC, Russell JB. Nisin resistance of Streptococcus bovis. Appl Environ Microbiol. 2001;67: 808–13. pmid:11157247
  17. 17. Gevers D, Huys G, Swings J. In vitro conjugal transfer of tetracycline resistance from Lactobacillusisolates to other Gram-positive bacteria. FEMS Microbiol Lett. 2003;225: 125–130. pmid:12900030
  18. 18. O’Driscoll J, Glynn F, Fitzgerald GF, Van Sinderen D. Sequence analysis of the lactococcal plasmid pNP40: A mobile replicon for coping with environmental hazards. J Bacteriol. 2006;188: 6629–6639. pmid:16952955
  19. 19. Sterkenburg A, Van Leeuwen P, Wouters JTM. Loss of phage resistance encoded by plasmid pSK112 in chemostat cultures of Lactococcus lactis ssp. cremoris SK110. Biochimie. 1988;70: 451–456. pmid:3139065
  20. 20. Millen AM, Horvath P, Boyaval P, Romero D. Mobile CRISPR/Cas-mediated bacteriophage resistance in Lactococcus lactis. Janssen PJ, editor. PLoS One. 2012;7: e51663. pmid:23240053
  21. 21. Meyrand M, Guillot A, Goin M, Furlan S, Armalyte J, Kulakauskas S, et al. Surface proteome analysis of a natural isolate of Lactococcus lactis reveals the presence of pili able to bind human intestinal epithelial cells. Mol Cell Proteomics. 2013;12: 3935–47. pmid:24002364
  22. 22. Oxaran V, Ledue-Clier F, Dieye Y, Herry J-M, Péchoux C, Meylheuc T, et al. Pilus biogenesis in Lactococcus lactis: molecular characterization and role in aggregation and biofilm formation. Biswas I, editor. PLoS One. 2012;7: 1–18.
  23. 23. Kojic M, Jovcic B, Strahinic I, Begovic J, Lozo J, Veljovic K, et al. Cloning and expression of a novel lactococcal aggregation factor from Lactococcus lactis subsp. lactis BGKP1. BMC Microbiol. 2011;11: 265. pmid:22182285
  24. 24. Le DTL, Tran TL, Duviau MP, Meyrand M, Guérardel Y, Castelain M, et al. Unraveling the role of surface mucus-binding protein and pili in muco-adhesion of Lactococcus lactis. PLoS One. 2013;8: e79850. pmid:24260308
  25. 25. Linares DM, Kok J, Poolman B. Genome sequences of Lactococcus lactis MG1363 (revised) and NZ9000 and comparative physiological studies. J Bacteriol. 2010;192: 5806–12. pmid:20639323
  26. 26. Gasson MJ. Plasmid complements of Streptococcus lactis NCDO 712 and other lactic streptococci after protoplast-induced curing. J Bacteriol. 1983;154: 1–9. Available: http://www.ncbi.nlm.nih.gov/pubmed/6403500 pmid:6403500
  27. 27. Wegmann U, Overweg K, Jeanson S, Gasson M, Shearman C. Molecular characterization and structural instability of the industrially important composite metabolic plasmid pLP712. Microbiology. 2012;158: 2936–45. pmid:23023974
  28. 28. Gasson MJ, Davies FL. High-frequency conjugation associated with Streptococcus lactis donor cell aggregation. J Bacteriol. 1980;143: 1260–1264. pmid:6773926
  29. 29. Gasson MJ, Swindell S, Maeda S, Dodd HM. Molecular rearrangement of lactose plasmid DNA associated with high-frequency transfer and cell aggregation in Lactococcus lactis 712. Mol Microbiol. 1992;6: 3213–23. pmid:1453959
  30. 30. Klaenhammer TR, Sanozky RB. Conjugal transfer from Streptococcus lactis ME2 of plasmids encoding phage resistance, nisin resistance and lactose-fermenting ability: evidence for a high-frequency conjugative plasmid responsible for abortive infection of virulent bacteriophage. J Gen Microbiol. 1985;131: 1531–41. pmid:3930657
  31. 31. Anderson DG, McKay LL. Genetic and physical characterization of recombinant plasmids associated with cell aggregation and high-frequency conjugal transfer in Streptococcus lactis ML3. J Bacteriol. 1984;158: 954–62. Available: http://www.ncbi.nlm.nih.gov/pubmed/6327653 pmid:6327653
  32. 32. Walsh PM, McKay LL. Recombinant plasmid associated with cell aggregation and high-frequency conjugation of Streptococcus lactis ML3. J Bacteriol. 1981;146: 937–944. pmid:6787018
  33. 33. Van Der Lelie D, Chavarri F, Venema G, Gasson MJ. Identification of a new genetic determinant for cell aggregation associated with lactose plasmid transfer in Lactococcus lactis. Appl Environ Microbiol. 1991;57: 201–206. pmid:1903626
  34. 34. Stentz R, Jury K, Eaton T, Parker M, Narbad A, Gasson M, et al. Controlled expression of CluA in Lactococcus lactis and its role in conjugation. Microbiology. 2004;150: 2503–12. pmid:15289547
  35. 35. Godon J-J, Jury K, Shearman CA, Gasson MJ. The Lactococcus lactis sex-factor aggregation gene cluA. Mol Microbiol. Wiley Online Library; 1994;12: 655–663.
  36. 36. Proft T, Baker EN. Pili in Gram-negative and Gram-positive bacteria—Structure, assembly and their role in disease. Cellular and Molecular Life Sciences. 2009. pp. 613–635. pmid:18953686
  37. 37. Filloux A. A variety of bacterial pili involved in horizontal gene transfer. Journal of Bacteriology. 2010. pp. 3243–3245. pmid:20418394
  38. 38. Scott JR, Zähner D. Pili with strong attachments: Gram-positive bacteria do it differently. Molecular Microbiology. 2006. pp. 320–330. pmid:16978260
  39. 39. Mandlik A, Swierczynski A, Das A, Ton-That H. Pili in Gram-positive bacteria: assembly, involvement in colonization and biofilm development. Trends in Microbiology. 2008. pp. 33–40. pmid:18083568
  40. 40. Tripathi P, Beaussart A, Andre G, Rolain T, Lebeer S, Vanderleyden J, et al. Towards a nanoscale view of lactic acid bacteria. Micron. Elsevier Ltd; 2012;43: 1323–1330.
  41. 41. Telford JL, Barocchi M a, Margarit I, Rappuoli R, Grandi G. Pili in Gram-positive pathogens. Nat Rev Microbiol. 2006;4: 509–519. pmid:16778837
  42. 42. Von Ossowski I, Reunanen J, Satokari R, Vesterlund S, Kankainen M, Huhtinen H, et al. Mucosal adhesion properties of the probiotic Lactobacillus rhamnosus GG SpaCBA and SpaFED pilin subunits. Appl Environ Microbiol. 2010;76: 2049–2057. pmid:20118368
  43. 43. Kankainen M, Paulin L, Tynkkynen S, von Ossowski I, Reunanen J, Partanen P, et al. Comparative genomic analysis of Lactobacillus rhamnosus GG reveals pili containing a human- mucus binding protein. Proc Natl Acad Sci U S A. 2009;106: 17193–17198. pmid:19805152
  44. 44. Taguchi F, Ichinose Y. Role of type IV pili in virulence of Pseudomonas syringae pv. tabaci 6605: correlation of motility, multidrug resistance, and HR-inducing activity on a nonhost plant. Mol Plant Microbe Interact. 2011;24: 1001–1011. pmid:21615203
  45. 45. De Vos W. Gene cloning and expression in lactic streptococci. FEMS Microbiol Lett. 1987;46: 281–295.
  46. 46. Leenhouts KJ, Kok J, Venema G. Lactococcal plasmid pWVO1 as an integration vector for lactococci. Appl Environ Microbiol. 1991;57: 2562–2567. pmid:1768128
  47. 47. Mierau I, Kleerebezem M. 10 Years of the nisin-controlled gene expression system (NICE) in Lactococcus lactis. Applied Microbiology and Biotechnology. 2005. pp. 705–717. pmid:16088349
  48. 48. De Ruyter PGGA, Kuipers OP, De Vos WM. Controlled gene expression systems for Lactococcus lactis with the food- grade inducer nisin. Appl Environ Microbiol. 1996;62: 3662–3667. pmid:8837421
  49. 49. Mierau I, Leij P, van Swam I, Blommestein B, Floris E, Mond J, et al. Industrial-scale production and purification of a heterologous protein in Lactococcus lactis using the nisin-controlled gene expression system NICE: the case of lysostaphin. Microb Cell Fact. 2005;4: 15. pmid:15921518
  50. 50. Kok J, Van Der Vossen JMBM, Venema G. Construction of plasmid cloning vectors for lactic streptococci which also replicate in Bacillus subtilis and Escherichia coli. Appl Environ Microbiol. 1984;48: 726–731. pmid:6095756
  51. 51. Kuipers OP, De Ruyter PGG a, Kleerebezem M, De Vos WM. Controlled overproduction of proteins by lactic acid bacteria. Trends in Biotechnology. 1997. pp. 135–140. pmid:9131833
  52. 52. Simon D, Chopin A. Construction of a vector plasmid family and its use for molecular cloning in Streptococcus lactis. Biochimie. 1988;70: 559–566. pmid:2844302
  53. 53. Sambrook J, Fritsch EF, Maniatis T. Molecular cloning : a laboratory manual. Cold Spring Habor Laboratory press; 1989.
    • 54. Zerbino DR, Birney E. Velvet: Algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 2008;18: 821–829. pmid:18349386
    • 55. Zerbino DR. Using the Velvet de novo Assembler for Short-Read Sequencing Technologies. Current Protocols in Bioinformatics. Hoboken, NJ, USA: John Wiley & Sons, Inc.; 2010. pp. 1–13.
      • 56. Aziz RK, Bartels D, Best A a, DeJongh M, Disz T, Edwards R a, et al. The RAST Server: rapid annotations using subsystems technology. BMC Genomics. 2008;9: 75. pmid:18261238
      • 57. Rutherford K, Parkhill J, Crook J, Horsnell T, Rice P, Rajandream M a, et al. Artemis: sequence visualization and annotation. Bioinformatics. 2000;16: 944–945. pmid:11120685
      • 58. Carver T, Harris SR, Berriman M, Parkhill J, McQuillan J a. Artemis: An integrated platform for visualization and analysis of high-throughput sequence-based experimental data. Bioinformatics. 2012;28: 464–469. pmid:22199388
      • 59. Barrick JE, Colburn G, Deatherage DE, Traverse CC, Strand MD, Borges JJ, et al. Identifying structural variation in haploid microbial genomes from short-read resequencing data using breseq. BMC Genomics. 2014;15: 1039. pmid:25432719
      • 60. Kumar P, Henikoff S, Ng PC. Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm. Nat Protoc. 2009;4: 1073–1081. pmid:19561590
      • 61. Wells JM, Wilson PW, Le Page RW. Improved cloning vectors and transformation procedure for Lactococcus lactis. The Journal of applied bacteriology. 1993. pp. 629–36. pmid:8349525
      • 62. Holo H, Nes IF. High-frequency transformation, by electroporation, of Lactococcus lactis subsp. cremoris grown with glycine in osmotically stabilized media. Appl Environ Microbiol. 1989;55: 3119–3123. pmid:16348073
      • 63. Ingham CJ, Van Den Ende M, Wever PC, Schneeberger PM. Rapid antibiotic sensitivity testing and trimethoprim-mediated filamentation of clinical isolates of the Enterobacteriaceae assayed on a novel porous culture support. J Med Microbiol. 2006;55: 1511–1519. pmid:17030910
      • 64. McKay LL, Baldwin K a., Walsh PM. Conjugal transfer of genetic information in group N streptococci. Appl Environ Microbiol. 1980;40: 84–91. pmid:6773476
      • 65. Steele JL, McKay LL. Partial characterization of the genetic basis for sucrose metabolism and nisin production in Streptococcus lactis. Appl Environ Microbiol. 1986;51: 57–64. pmid:3082279
      • 66. Wegmann U, O’Connell-Motherway M, Zomer A, Buist G, Shearman C, Canchaya C, et al. Complete genome sequence of the prototype lactic acid bacterium Lactococcus lactis subsp. cremoris MG1363. J Bacteriol. 2007;189: 3256–3270. pmid:17307855
      • 67. van der Meulen SB, de Jong A, Kok J. Transcriptome landscape of Lactococcus lactis reveals many novel RNAs including a small regulatory RNA involved in carbon uptake and metabolism. RNA Biol. 2016;13: 353–66. pmid:26950529
      • 68. Sim N-L, Kumar P, Hu J, Henikoff S, Schneider G, Ng PC. SIFT web server: predicting effects of amino acid substitutions on proteins. Nucleic Acids Res. 2012;40: W452–7. pmid:22689647
      • 69. Le Bourgeois P, Daveran-Mingot ML, Ritzenthaler P. Genome plasticity among related Lactococcus strains: Identification of genetic events associated with macrorestriction polymorphisms. J Bacteriol. 2000;182: 2481–2491. pmid:10762249
      • 70. Duan K, Liu CQ, Liu YJ, Ren J, Dunn NW. Nucleotide sequence and thermostability of pND324, a 3.6-kb plasmid from Lactococcus lactis. Appl Microbiol Biotechnol. 1999;53: 36–42. pmid:10645623
      • 71. Montville TJ, Chen Y. Mechanistic action of pediocin and nisin: Recent progress and unresolved questions. Appl Microbiol Biotechnol. 1998;50: 511–519. pmid:9917136
      • 72. Van der Meer JR, Polman J, Beerthuyzen MM, Siezen RJ, Kuipers OP, De Vos WM. Characterization of the Lactococcus lactis nisin A operon genes nisP, encoding a subtilisin-like serine protease involved in precursor processing, and nisR, encoding a regulatory protein involved in nisin biosynthesis. J Bacteriol. 1993;175: 2578–2588. pmid:8478324
      • 73. de Ruyter PG, Kuipers OP, Beerthuyzen MM, van Alen-Boerrigter I, de Vos WM. Functional analysis of promoters in the nisin gene cluster of Lactococcus lactis. J Bacteriol. 1996;178: 3434–9. Available: http://www.ncbi.nlm.nih.gov/pubmed/8655538 pmid:8655538
      • 74. Kelly WJ, Davey GP, Ward LJ. Characterization of lactococci isolated from minimally processed fresh fruit and vegetables. Int J Food Microbiol. 1998;45: 85–92. Available: http://www.ncbi.nlm.nih.gov/pubmed/9924939 pmid:9924939
      • 75. Efstathiou JD, McKay LL. Inorganic salts resistance associated with a lactose-fermenting plasmid in Streptococcus lactis. J Bacteriol. 1977;130: 257–65. pmid:404284
      • 76. Liu CQ, Charoechai P, Khunajakr N, Deng YM, Widodo , Dunn NW. Genetic and transcriptional analysis of a novel plasmid-encoded copper resistance operon from Lactococcus lactis. Gene. 2002;297: 241–247. pmid:12384305
      • 77. Davies FL, Underwood HM, Gasson MJ. The Value of Plasmid Profiles for Strain Identification in Lactic Streptococci and the Relationship between Streptoccocus lactis 712, ML3 and C2. J Appl Bacteriol. 1981;51: 325–337.
      • 78. Ton-That H, Schneewind O. Assembly of pili in Gram-positive bacteria. Trends in Microbiology. 2004. pp. 228–234. pmid:15120142
      • 79. Reunanen J, von Ossowski I, Hendrickx AP a, Palva A, de Vosa WM. Characterization of the SpaCBA pilus fibers in the probiotic Lactobacillus rhamnosus GG. Appl Environ Microbiol. 2012;78: 2337–2344. pmid:22247175
      • 80. Douillard FP, Rasinkangas P, von Ossowski I, Reunanen J, Palva A, de Vos WM. Functional Identification of Conserved Residues Involved in Lactobacillus rhamnosus Strain GG Sortase Specificity and Pilus Biogenesis. J Biol Chem. 2014;289: 15764–15775. pmid:24753244
      • 81. Ton-That H, Marraffini L a, Schneewind O. Protein sorting to the cell wall envelope of Gram-positive bacteria. Biochim Biophys Acta. 2004;1694: 269–78. pmid:15546671
      • 82. Stone BJ, Abu Kwaik Y. Expression of multiple pili by Legionella pneumophila: identification and characterization of a type IV pilin gene and its role in adherence to mammalian and protozoan cells. Infect Immun. 1998;66: 1768–1775. pmid:9529112
      • 83. Hendrickx AP a, Bonten MJM, van Luit-Asbroek M, Schapendonk CME, Kragten AHM, Willems RJL. Expression of two distinct types of pili by a hospital-acquired Enterococcus faecium isolate. Microbiology. 2008;154: 3212–23. pmid:18832326
      • 84. Kline KA, Dodson KW, Caparon MG, Hultgren SJ. A tale of two pili: assembly and function of pili in bacteria. Trends Microbiol. 2010;18: 224–232. pmid:20378353
      • 85. Marraffini LA, DeDent AC, Schneewind O. Sortases and the Art of Anchoring Proteins to the Envelopes of Gram-Positive Bacteria. Microbiol Mol Biol Rev. 2006;70: 192–221. pmid:16524923
      • 86. Cozzi R, Nuccitelli A, D’Onofrio M, Necchi F, Rosini R, Zerbini F, et al. New insights into the role of the glutamic acid of the E-box motif in group B Streptococcus pilus 2a assembly. FASEB J. 2012;26: 2008–18. pmid:22253480
      • 87. Ton-That H, Marraffini L a, Schneewind O. Sortases and pilin elements involved in pilus assembly of Corynebacterium diphtheriae. Mol Microbiol. 2004;53: 251–261. pmid:15225319
      • 88. Rintahaka J, Yu X, Kant R, Palva A, Von Ossowski I. Phenotypical analysis of the Lactobacillus rhamnosus GG fimbrial spafed operon: Surface expression and functional characterization of recombinant SpaFED pili in Lactococcus lactis. PLoS One. 2014;9: e113922. pmid:25415357
      • 89. Clarke M, Maddera L, Harris RL, Silverman PM. F-pili dynamics by live-cell imaging. Proc Natl Acad Sci U S A. 2008;105: 17978–17981. pmid:19004777
      • 90. Bachmann H, Starrenburg MJC, Molenaar D, Kleerebezem M, Vlieg JET van H, van Hylckama Vlieg JET. Microbial domestication signatures of Lactococcus lactis can be reproduced by experimental evolution. Genome Res. 2012;22: 115–24. pmid:22080491
      • 91. Khunajakr N, Liu CQ, Charoenchai P, Dunn NW. A plasmid-encoded two-component regulatory system involved in copper-inducible transcription in Lactococcus lactis. Gene. 1999;229: 229–35. pmid:10095123