The Plasmid Complement of Lactococcus lactis UC509.9 Encodes Multiple Bacteriophage Resistance Systems

Lactococcus lactis subsp. cremoris strains are used globally for the production of fermented dairy products, particularly hard cheeses. Believed to be of plant origin, L. lactis strains that are used as starter cultures have undergone extensive adaptation to the dairy environment, partially through the acquisition of extrachromosomal DNA in the form of plasmids that specify technologically important phenotypic traits. Here, we present a detailed analysis of the eight plasmids of L. lactis UC509.9, an Irish dairy starter strain. Key industrial phenotypes were mapped, and genes that are typically associated with lactococcal plasmids were identified. Four distinct, plasmid-borne bacteriophage resistance systems were identified, including two abortive infection systems, AbiB and AbiD1, thereby supporting the observed phage resistance of L. lactis UC509.9. AbiB escape mutants were generated for phage sk1, which were found to carry mutations in orf6, which encodes the major capsid protein of this phage.     

Identification and Characterization of Lactococcal-Prophage-Carried Superinfection Exclusion Genes

Superinfection exclusion (Sie) proteins are prophage-encoded phage resistance systems. In this study, genes encoding Sie systems were identified on the genomes of Lactococcus lactis subsp. cremoris MG1363 and L. lactis subsp. lactis IL1403. These Sie systems are genetically distinct and yet were shown to act specifically against a particular subset of the 936 phage group. Each of the systems allows normal phage adsorption while affecting plasmid transduction and intracellular phage DNA replication, which points to the blocking of phage DNA injection as their common mode of action. Sie-specifying genes found on the MG1363 prophages are also present in various lactococcal strains, whereas the prophage-encoded Sie systems of IL1403 do not appear to be as widely disseminated.     

Mobile CRISPR/Cas-Mediated Bacteriophage Resistance in Lactococcus lactis

Lactococcus lactis is a biotechnological workhorse for food fermentations and potentially therapeutic products and is therefore widely consumed by humans. It is predominantly used as a starter microbe for fermented dairy products, and specialized strains have adapted from a plant environment through reductive evolution and horizontal gene transfer as evidenced by the association of adventitious traits with mobile elements. Specifically, L. lactis has armed itself with a myriad of plasmid-encoded bacteriophage defensive systems to protect against viral predation. This known arsenal had not included CRISPR/Cas (clustered regularly interspaced short palindromic repeats/CRISPR-associated proteins), which forms a remarkable microbial immunity system against invading DNA. Although CRISPR/Cas systems are common in the genomes of closely related lactic acid bacteria (LAB), none was identified within the eight published lactococcal genomes. Furthermore, a PCR-based search of the common LAB CRISPR/Cas systems (Types I and II) in 383 industrial L. lactis strains proved unsuccessful. Here we describe a novel, Type III, self-transmissible, plasmid-encoded, phage-interfering CRISPR/Cas discovered in L. lactis. The native CRISPR spacers confer resistance based on sequence identity to corresponding lactococcal phage. The interference is directed at phages problematic to the dairy industry, indicative of a responsive system. Moreover, targeting could be modified by engineering the spacer content. The 62.8-kb plasmid was shown to be conjugally transferrable to various strains. Its mobility should facilitate dissemination within microbial communities and provide a readily applicable system to naturally introduce CRISPR/Cas to industrially relevant strains for enhanced phage resistance and prevention against acquisition of undesirable genes.     

Identification of Lactococcus lactis Genes Required for Bacteriophage Adsorption

The aim of this work was to identify genes in Lactococcus lactis subsp. lactis IL1403 and Lactococcus lactis subsp. cremoris Wg2 important for adsorption of the 936-species phages bIL170 and 645, respectively. Random insertional mutagenesis of the two L. lactis strains was carried out with the vector pGh9:ISS1, and integrants that were resistant to phage infection and showed reduced phage adsorption were selected. In L. lactis IL1403 integration was obtained in the ycaG and rgpE genes, whereas in L. lactis Wg2 integration was obtained in two genes homologous to ycbC and ycbB of L. lactis IL1403. rgpE and ycbB encode putative glycosyltransferases, whereas ycaG and ycbC encode putative membrane-spanning proteins with unknown functions. Interestingly, ycaG, rgpE, ycbC, and ycbB are all part of the same operon in L. lactis IL1403. This operon is probably involved in biosynthesis and transport of cell wall polysaccharides (WPS). Binding and infection studies showed that 645 binds to and infects L. lactis Wg2, L. lactis IL1403, and L. lactis IL1403 strains with pGh9:ISS1 integration in ycaG and rgpE, whereas bIL170 binds to and infects only L. lactis IL1403 and cannot infect Wg2. These results indicate that 645 binds to a WPS structure present in both L. lactis IL1403 and L. lactis Wg2, whereas bIL170 binds to another WPS structure not present in L. lactis Wg2. Binding of bIL170 and 645 to different WPS structures was supported by alignment of the receptor-binding proteins of bIL170 and 645 that showed no homology in the C-terminal part.     

Identification, Cloning, and Characterization of a Lactococcus lactis Branched-Chain α-Keto Acid Decarboxylase Involved in Flavor Formation

The biochemical pathway for formation of branched-chain aldehydes, which are important flavor compounds derived from proteins in fermented dairy products, consists of a protease, peptidases, a transaminase, and a branched-chain α-keto acid decarboxylase (KdcA). The activity of the latter enzyme has been found only in a limited number of Lactococcus lactis strains. By using a random mutagenesis approach, the gene encoding KdcA in L. lactis B1157 was identified. The gene for this enzyme is highly homologous to the gene annotated ipd, which encodes a putative indole pyruvate decarboxylase, in L. lactis IL1403. Strain IL1403 does not produce KdcA, which could be explained by a 270-nucleotide deletion at the 3′ terminus of the ipd gene encoding a truncated nonfunctional decarboxylase. The kdcA gene was overexpressed in L. lactis for further characterization of the decarboxylase enzyme. Of all of the potential substrates tested, the highest activity was observed with branched-chain α-keto acids. Moreover, the enzyme activity was hardly affected by high salinity, and optimal activity was found at pH 6.3, indicating that the enzyme might be active under cheese ripening conditions.     

Projector: automatic contig mapping for gap closure purposes

Projector was designed for automatic positioning of contigs from an unfinished prokaryotic genome onto a template genome of a closely related strain or species. Projector mapped 84 contigs of Lactococcus lactis MG1363 (corresponding to 81% of the assembly nucleotides) against the genome of L.lactis IL1403. Ninety three percent of subsequent gap closure PCRs were successful. Moreover, a significant improvement in the N50 and N80 values (describing the assembly quality) was observed after the use of Projector. Because increasing numbers of bacterial genomes are being sequenced, Projector provides an efficient method to close a significant number of remaining gaps in the late stages of a genome sequencing project.     

Expression of prophage-encoded endolysins contributes to autolysis of Lactococcus lactis

Analysis of autolysis of derivatives of Lactococcus lactis subsp. cremoris MG1363 and subsp. lactis IL1403, both lacking the major autolysin AcmA, showed that L. lactis IL1403 still lysed during growth while L. lactis MG1363 did not. Zymographic analysis revealed that a peptidoglycan hydrolase activity of around 30 kDa is present in cell extracts of L. lactis IL1403 that could not be detected in strain MG1363. A comparison of all genes encoding putative peptidoglycan hydrolases of IL1403 and MG1363 led to the assumption that one or more of the 99 % homologous 27.9-kDa endolysins encoded by the prophages bIL285, bIL286 and bIL309 could account for the autolysis phenotype of IL1403. Induced expression of the endolysins from bIL285, bIL286 or bIL309 in L. lactis MG1363 resulted in detectable lysis or lytic activity. Prophage deletion and insertion derivatives of L. lactis IL1403 had a reduced cell lysis phenotype. RT-qPCR and zymogram analysis showed that each of these strains still expressed one or more of the three phage lysins. A homologous gene and an endolysin activity were also identified in the natural starter culture L. lactis subsp. cremoris strains E8, Wg2 and HP, and the lytic activity could be detected under growth conditions that were identical as those used for IL1403. The results presented here show that these endolysins of L. lactis are expressed during normal growth and contribute to autolysis without production of (lytic) phages. Screening for natural strains expressing homologous endolysins could help in the selection of strains with enhanced autolysis and, thus, cheese ripening properties.

     

A Virulent Phage Infecting Lactococcus garvieae,with Homology to Lactococcus lactis Phages

A new virulent phage belonging to the Siphoviridae family and able to infect Lactococcus garvieae strains was isolated from compost soil. Phage GE1 has a prolate capsid (56 by 38 nm) and a long noncontractile tail (123 nm). It had a burst size of 139 and a latent period of 31 min. Its host range was limited to only two L. garvieae strains out of 73 tested. Phage GE1 has a double-stranded DNA genome of 24,847 bp containing 48 predicted open reading frames (ORFs). Putative functions could be assigned to only 14 ORFs, and significant matches in public databases were found for only 17 ORFs, indicating that GE1 is a novel phage and its genome contains several new viral genes and encodes several new viral proteins. Of these 17 ORFs, 16 were homologous to deduced proteins of virulent phages infecting the dairy bacterium Lactococcus lactis, including previously characterized prolate-headed phages. Comparative genome analysis confirmed the relatedness of L. garvieae phage GE1 to L. lactis phages c2 (22,172 bp) and Q54 (26,537 bp), although its genome organization was closer to that of phage c2. Phage GE1 did not infect any of the 58 L. lactis strains tested. This study suggests that phages infecting different lactococcal species may have a common ancestor.     

Diversity analysis of dairy and non-dairy strains of Lactococcus lactis ssp. lactis by multilocus sequence analysis (MLSA)

The genetic diversity of 31 identified strains of Lactococcus lactis ssp. lactis isolated from different dairy and non-dairy sources were investigated at gene level using multilocus sequence analysis (MLSA) and PCR-RFLP based on the differences in four selected partial protein coding gene sequences: araT, encoding aromatic amino acid-specific aminotransferase; dtpT, encoding di/tri peptide transporter; yueF, encoding non-proteolytic protein, peptidase, M16 family; and pdhA, encoding pyruvate dehydrogenase E1 component α-subunit. A set of seven test strains from different isolation sources and one reference strain, L. lactis ssp. lactis NCDC 094, were analyzed by MLSA. The strains showed distinct diversity among themselves and exhibited a greater percent similarity with reference strains L. lactis ssp. lactis CV56 (CP002365.1), IL1403 (AE005176.1), and KF147 (CP001834.1) in comparison with L. lactis ssp. cremoris NZ9000 (CP002094.1), MG1363 (AM406671.1), and SK11 (CP00425.1). The MLSA revealed one distinct genomic lineage within strains exclusively of L. lactis ssp. lactis. This analysis also revealed no source-wise genetic relationship in the test strains analyzed. Further, PCR-RFLP of araT, dtpT, yueF and pdhA also characterized the single genomic lineage exclusively of L. lactis ssp. lactis within a total of 24 test strains.     

Transformation of Folate-Consuming Lactobacillus gasseri into a Folate Producer

Five genes essential for folate biosynthesis in Lactococcus lactis were cloned on a broad-host-range lactococcal vector and were transferred to the folate auxotroph Lactobacillus gasseri. As a result L. gasseri changed from a folate consumer to a folate producer. This principle can be used to increase folate levels in many fermented food products.     

pSEUDO, a Genetic Integration Standard for Lactococcus lactis

Plasmid pSEUDO and derivatives were used to show that llmg_pseudo_10 in Lactococcus lactis MG1363 and its homologous locus in L. lactis IL1403 are suitable for chromosomal integrations. L. lactis MG1363 and IL1403 nisin-induced controlled expression (NICE) system derivatives (JP9000 and IL9000) and two general stress reporter strains (NZ9000::PhrcA-GFP and NZ9000::PgroES-GFP) enabling in vivo noninvasive monitoring of cellular fitness were constructed.     

Cloning,Expression, and Functional Characterization of Secondary Amino Acid Transporters of Lactococcus lactis

Fourteen genes encoding putative secondary amino acid transporters were identified in the genomes of Lactococcus lactis subsp. cremoris strains MG1363 and SK11 and L. lactis subsp. lactis strains IL1403 and KF147, 12 of which were common to all four strains. Amino acid uptake in L. lactis cells overexpressing the genes revealed transporters specific for histidine, lysine, arginine, agmatine, putrescine, aromatic amino acids, acidic amino acids, serine, and branched-chain amino acids. Substrate specificities were demonstrated by inhibition profiles determined in the presence of excesses of the other amino acids. Four knockout mutants, lacking the lysine transporter LysP, the histidine transporter HisP (formerly LysQ), the acidic amino acid transporter AcaP (YlcA), or the aromatic amino acid transporter FywP (YsjA), were constructed. The LysP, HisP, and FywP deletion mutants showed drastically decreased rates of uptake of the corresponding substrates at low concentrations. The same was observed for the AcaP mutant with aspartate but not with glutamate. In rich M17 medium, the deletion of none of the transporters affected growth. In contrast, the deletion of the HisP, AcaP, and FywP transporters did affect growth in a defined medium with free amino acids as the sole amino acid source. HisP was essential at low histidine concentrations, and AcaP was essential in the absence of glutamine. FywP appeared to play a role in retaining intracellularly synthesized aromatic amino acids when these were not added to the medium. Finally, HisP, AcaP, and FywP did not play a role in the excretion of accumulated histidine, glutamate, or phenylalanine, respectively, indicating the involvement of other transporters.     

Lactococcus garvieae: Where Is It From? A First Approach to Explore the Evolutionary History of This Emerging Pathogen

The population structure and diversity of Lactococcus garvieae, an emerging pathogen of increasing clinical significance, was determined at both gene and genome level. Selected lactococcal isolates of various origins were analyzed by a multi locus sequence typing (MLST). This gene-based analysis was compared to genomic characteristics, estimated through the complete genome sequences available in database. The MLST identified two branches containing the majority of the strains and two branches bearing one strain each. One strain was particularly differentiated from the other L. garvieae strains, showing a significant genetic distance. The genomic characteristics, correlated to the MLST-based phylogeny, indicated that this “separated strain” appeared first and could be considered the evolutionary intermediate between Lactococcus lactis and L. garvieae main clusters. A preliminary genome analysis of L. garvieae indicated a pan-genome constituted of about 4100 genes, which included 1341 core genes and 2760 genes belonging to the dispensable genome. A total of 1491 Clusters of Orthologous Genes (COGs) were found to be specific to the 11 L. garvieae genomes, with the genome of the “separated strain” showing the highest presence of unique genes.     

Genes but not genomes reveal bacterial domestication of Lactococcus lactis

Background

The population structure and diversity of Lactococcus lactis subsp. lactis, a major industrial bacterium involved in milk fermentation, was determined at both gene and genome level. Seventy-six lactococcal isolates of various origins were studied by different genotyping methods and thirty-six strains displaying unique macrorestriction fingerprints were analyzed by a new multilocus sequence typing (MLST) scheme. This gene-based analysis was compared to genomic characteristics determined by pulsed-field gel electrophoresis (PFGE).

Methodology/Principal Findings

The MLST analysis revealed that L. lactis subsp. lactis is essentially clonal with infrequent intra- and intergenic recombination; also, despite its taxonomical classification as a subspecies, it displays a genetic diversity as substantial as that within several other bacterial species. Genome-based analysis revealed a genome size variability of 20%, a value typical of bacteria inhabiting different ecological niches, and that suggests a large pan-genome for this subspecies. However, the genomic characteristics (macrorestriction pattern, genome or chromosome size, plasmid content) did not correlate to the MLST-based phylogeny, with strains from the same sequence type (ST) differing by up to 230 kb in genome size.

Conclusion/Significance

The gene-based phylogeny was not fully consistent with the traditional classification into dairy and non-dairy strains but supported a new classification based on ecological separation between “environmental” strains, the main contributors to the genetic diversity within the subspecies, and “domesticated” strains, subject to recent genetic bottlenecks. Comparison between gene- and genome-based analyses revealed little relationship between core and dispensable genome phylogenies, indicating that clonal diversification and phenotypic variability of the “domesticated” strains essentially arose through substantial genomic flux within the dispensable genome.     

A one step genotypic identification of Lactococcus lactis subspecies at the species/strain levels

A Lactococcus lactis subspecies-specific primer was designed based on their repetitive genome sequences. This primer enabled L. lactis subspecies to be identified simultaneously at both the species level and also the strain level. Based on studies using 70 strains of L. lactis and 60 strains of other non-target bacteria, the identification completely matched that obtained by the sequence of the 16S rRNA gene. However, inconsistency between phenotypic and genotypic characteristics was observed in some strains isolated from milk.