Physical and genetic map of the Lactococcus lactis subsp. cremoris MG1363 chromosome: comparison with that of Lactococcus lactis subsp. lactis IL 1403 reveals a large genome inversion

Harnessing diversity of Lactococcus lactis from raw goat milk: Design of an indigenous starter for the production of Rocamadour, a French PDO cheese

Twenty-four strains of Lactococcus lactis isolated from raw goat milk collected in the Rocamadour PDO area were analysed by MLST typing and phenotypic characterisation. The strains were combined to design an indigenous starter for the production of Rocamadour PDO cheese. The strains were divided into three classes based on their technological properties: acidifying and proteolytic strains in class I (12/24 strains), slightly acidifying and non-proteolytic strains in class II (2/24 strains), and non-acidifying and non-proteolytic strains in class III (10/24 strains). Interestingly, all but three strains (21/24) produced diacetyl/acetoin despite not having citrate metabolism genes, as would classically be expected for the production of these aroma compounds. Three strains (EIP07A, EIP13D, and EIP20B) were selected for the indigenous starter based on the following inclusion/exclusion criteria: (i) no negative interactions between included strains, (ii) ability to metabolize lactose and at least one strain with the prtP gene and/or capable of producing diacetyl/acetoin, and (iii) selected strains derived from different farms to maximise genetic and phenotypic diversity. Despite consisting exclusively of L. lactis strains, the designed indigenous starter allowed reproducible cheese production with performances similar to those obtained with an industrial starter and with the sensory qualities expected of Rocamadour PDO cheese.

Large-scale chromosome flip-flop reversible inversion mediates phenotypic switching of expression of antibiotic resistance in lactococci

Bacteria can gain resistance to antimicrobials by acquiring and expressing genetic elements that encode resistance determinants such as efflux pumps and drug-modifying enzymes, thus hampering treatment of infection. Previously we showed that acquisition of spectinomycin resistance in a lactococcal strain was correlated with a reversible genomic inversion, but the precise location and the genes affected were unknown. Here we use long-read whole-genome sequencing to precisely define the genomic inversion and we use quantitative PCR to identify associated changes in gene expression levels. The boundaries of the inversion fall within two identical copies of a prophage-like sequence, located on the left and right replichores; this suggests possible mechanisms for inversion through homologous recombination or prophage activity. The inversion is asymmetrical in respect of the axis between the origin and terminus of the replication and modulates the expression of a SAM-dependent methyltransferase, whose heterologous expression confers resistance to spectinomycin in lactococci and that is up-regulated on exposure to spectinomycin. This study provides one of the first examples of phase variation via large-scale chromosomal inversions that confers a switch in antimicrobial resistance in bacteria and the first outside of Staphylococcus aureus.

Genetics of Lactococci

Lactococcus lactis is the best characterized species among the lactococci, and among the most consumed food-fermenting bacteria worldwide. Thanks to their importance in industrialized food production, lactococci are among the lead bacteria understood for fundamental metabolic pathways that dictate growth and survival properties. Interestingly, lactococci belong to the Streptococcaceae family, which includes food, commensal and virulent species. As basic metabolic pathways (e.g., respiration, metal homeostasis, nucleotide metabolism) are now understood to underlie virulence, processes elucidated in lactococci could be important for understanding pathogen fitness and synergy between bacteria. This chapter highlights major findings in lactococci and related bacteria, and covers five themes: distinguishing features of lactococci, metabolic capacities including the less known respiration metabolism in Streptococcaceae, factors and pathways modulating stress response and fitness, interbacterial dialogue via metabolites, and novel applications in health and biotechnology.

Genome-wide detection of conservative site-specific recombination in bacteria

The ability of clonal bacterial populations to generate genomic and phenotypic heterogeneity is thought to be of great importance for many commensal and pathogenic bacteria. One common mechanism contributing to diversity formation relies on the inversion of small genomic DNA segments in a process commonly referred to as conservative site-specific recombination. This phenomenon is known to occur in several bacterial lineages, however it remains notoriously difficult to identify due to the lack of conserved features. Here, we report an easy-to-implement method based on high-throughput paired-end sequencing for genome-wide detection of conservative site-specific recombination on a single-nucleotide level. We demonstrate the effectiveness of the method by successfully detecting several novel inversion sites in an epidemic isolate of the enteric pathogen Clostridium difficile. Using an experimental approach, we validate the inversion potential of all detected sites in C. difficile and quantify their prevalence during exponential and stationary growth in vitro. In addition, we demonstrate that the master recombinase RecV is responsible for the inversion of some but not all invertible sites. Using a fluorescent gene-reporter system, we show that at least one gene from a two-component system located next to an invertible site is expressed in an on-off mode reminiscent of phase variation. We further demonstrate the applicability of our method by mining 209 publicly available sequencing datasets and show that conservative site-specific recombination is common in the bacterial realm but appears to be absent in some lineages. Finally, we show that the gene content associated with the inversion sites is diverse and goes beyond traditionally described surface components. Overall, our method provides a robust platform for detection of conservative site-specific recombination in bacteria and opens a new avenue for global exploration of this important phenomenon.

Comparative Genomics Reveals the Diversity of Restriction-Modification Systems and DNA Methylation Sites in Listeria monocytogenes

Listeria monocytogenes is a bacterial pathogen that is found in a wide variety of anthropogenic and natural environments. Genome sequencing technologies are rapidly becoming a powerful tool in facilitating our understanding of how genotype, classification phenotypes, and virulence phenotypes interact to predict the health risks of individual bacterial isolates. Currently, 57 closed L. monocytogenes genomes are publicly available, representing three of the four phylogenetic lineages, and they suggest that L. monocytogenes has high genomic synteny. This study contributes an additional 15 closed L. monocytogenes genomes that were used to determine the associations between the genome and methylome with host invasion magnitude. In contrast to previous findings, large chromosomal inversions and rearrangements were detected in five isolates at the chromosome terminus and within rRNA genes, including a previously undescribed inversion within rRNA-encoding regions. Each isolate's epigenome contained highly diverse methyltransferase recognition sites, even within the same serotype and methylation pattern. Eleven strains contained a single chromosomally encoded methyltransferase, one strain contained two methylation systems (one system on a plasmid), and three strains exhibited no methylation, despite the occurrence of methyltransferase genes. In three isolates a new, unknown DNA modification was observed in addition to diverse methylation patterns, accompanied by a novel methylation system. Neither chromosome rearrangement nor strain-specific patterns of epigenome modification observed within virulence genes were correlated with serotype designation, clonal complex, or in vitro infectivity. These data suggest that genome diversity is larger than previously considered in L. monocytogenes and that as more genomes are sequenced, additional structure and methylation novelty will be observed in this organism.

IMPORTANCE

Listeria monocytogenes is the causative agent of listeriosis, a disease which manifests as gastroenteritis, meningoencephalitis, and abortion. Among Salmonella, Escherichia coli, Campylobacter, and Listeria-causing the most prevalent foodborne illnesses-infection by L. monocytogenes carries the highest mortality rate. The ability of L. monocytogenes to regulate its response to various harsh environments enables its persistence and transmission. Small-scale comparisons of L. monocytogenes focusing solely on genome contents reveal a highly syntenic genome yet fail to address the observed diversity in phenotypic regulation. This study provides a large-scale comparison of 302 L. monocytogenes isolates, revealing the importance of the epigenome and restriction-modification systems as major determinants of L. monocytogenes phylogenetic grouping and subsequent phenotypic expression. Further examination of virulence genes of select outbreak strains reveals an unprecedented diversity in methylation statuses despite high degrees of genome conservation.

PFGE protocols to distinguish subspecies of Lactococcus lactis

Pulsed-field gel electrophoresis (PFGE), developed in the mid-1980s, rapidly became a "gold standard" method for analyzing bacterial chromosomes. Today, although outcompeted in resolution by alternative methods, such as optical mapping, and not applicable for high-throughput analyses, PFGE remains a valuable method for bacterial strain typing. Here, we describe optimized protocols for macrorestriction fingerprinting, characterization of plasmid complements, and gene localization by DNA-DNA hybridization of Lactococcus lactis genomes.

Novel exopolysaccharides produced by Lactococcus lactis subsp. lactis, and the diversity of epsE genes in the exopolysaccharide biosynthesis gene clusters

To characterize novel variations of exopolysaccharides (EPSs) produced by dairy strains of Lactococcus lactis subsp. lactis and subsp. cremoris, the EPSs of five dairy strains of L. lactis were purified. Sugar composition analysis showed two novel EPSs produced by strains of L. lactis subsp. lactis. One strain produced EPS lacking galactose, and the other produced EPS containing fucose. Among the eps gene clusters of these strains, the highly conserved epsD and its neighboring epsE were sequenced. Sequence and PCR analysis revealed that epsE genes were strain-specific. By Southern blot analysis using epsD, the eps gene cluster in each strain was found to locate to the chromosome or a very large plasmid. This is the first report on the identification of two novel EPSs in L. lactis subsp. lactis. The strains can be detected among other strains by using epsE genes specific to them.

Natural genetic transformation generates a population of merodiploids in Streptococcus pneumoniae

Partial duplication of genetic material is prevalent in eukaryotes and provides potential for evolution of new traits. Prokaryotes, which are generally haploid in nature, can evolve new genes by partial chromosome duplication, known as merodiploidy. Little is known about merodiploid formation during genetic exchange processes, although merodiploids have been serendipitously observed in early studies of bacterial transformation. Natural bacterial transformation involves internalization of exogenous donor DNA and its subsequent integration into the recipient genome by homology. It contributes to the remarkable plasticity of the human pathogen Streptococcus pneumoniae through intra and interspecies genetic exchange. We report that lethal cassette transformation produced merodiploids possessing both intact and cassette-inactivated copies of the essential target gene, bordered by repeats (R) corresponding to incomplete copies of IS861. We show that merodiploidy is transiently stimulated by transformation, and only requires uptake of a ∼3-kb DNA fragment partly repeated in the chromosome. We propose and validate a model for merodiploid formation, providing evidence that tandem-duplication (TD) formation involves unequal crossing-over resulting from alternative pairing and interchromatid integration of R. This unequal crossing-over produces a chromosome dimer, resolution of which generates a chromosome with the TD and an abortive chromosome lacking the duplicated region. We document occurrence of TDs ranging from ∼100 to ∼900 kb in size at various chromosomal locations, including by self-transformation (transformation with recipient chromosomal DNA). We show that self-transformation produces a population containing many different merodiploid cells. Merodiploidy provides opportunities for evolution of new genetic traits via alteration of duplicated genes, unrestricted by functional selective pressure. Transient stimulation of a varied population of merodiploids by transformation, which can be triggered by stresses such as antibiotic treatment in S. pneumoniae, reinforces the plasticity potential of this bacterium and transformable species generally.

Merodiploids are defined as cells possessing a partial duplication of their genetic material, which potentially allows evolution of new genes. Historically, some have been observed in studies of natural genetic transformation. Transformation allows the bacteria to take up foreign DNA and incorporate it into their genome by homology. It is key to the high diversity observed in the human pathogen Streptococcus pneumoniae (the pneumococcus). Here we show that transformation with self DNA generates a population of merodiploids with varied chromosomal duplications, up to almost half a genome in size. We show that formation of merodiploids by transformation only requires uptake of a small donor DNA fragment partially repeated in the chromosome. The donor repeat recombines with an alternative repeat on one arm of a replicating chromosome, whilst the non-repeated part recombines with its complement on the other arm, bridging the two. Subsequent recombination events generate a merodiploid chromosome with the region between the two repeats duplicated. Our results demonstrate that transformation, which is induced by stresses such as antibiotic treatments, transiently increases the ability of a population to form merodiploids. We suggest that creating a variety of merodiploids at a time of stress maximizes the adaptive potential of this pathogen.

Genotypic and phenotypic analysis of dairy Lactococcus lactis biodiversity in milk: volatile organic compounds as discriminating markers

The diversity of nine dairy strains of Lactococcus lactis subsp. lactis in fermented milk was investigated by both genotypic and phenotypic analyses. Pulsed-field gel electrophoresis and multilocus sequence typing were used to establish an integrated genotypic classification. This classification was coherent with discrimination of the L. lactis subsp. lactis bv. diacetylactis lineage and reflected clonal complex phylogeny and the uniqueness of the genomes of these strains. To assess phenotypic diversity, 82 variables were selected as important dairy features; they included physiological descriptors and the production of metabolites and volatile organic compounds (VOCs). Principal-component analysis (PCA) demonstrated the phenotypic uniqueness of each of these genetically closely related strains, allowing strain discrimination. A method of variable selection was developed to reduce the time-consuming experimentation. We therefore identified 20 variables, all associated with VOCs, as phenotypic markers allowing discrimination between strain groups. These markers are representative of the three metabolic pathways involved in flavor: lipolysis, proteolysis, and glycolysis. Despite great phenotypic diversity, the strains could be divided into four robust phenotypic clusters based on their metabolic orientations. Inclusion of genotypic diversity in addition to phenotypic characters in the classification led to five clusters rather than four being defined. However, genotypic characters make a smaller contribution than phenotypic variables (no genetic distances selected among the most contributory variables). This work proposes an original method for the phenotypic differentiation of closely related strains in milk and may be the first step toward a predictive classification for the manufacture of starters.

Genome analysis of food grade lactic Acid-producing bacteria: from basics to applications

Whole-genome sequencing has revolutionized and accelerated scientific research that aims to study the genetics, biochemistry and molecular biology of bacteria. Lactic acid-producing bacteria, which include lactic acid bacteria (LAB) and bifidobacteria, are typically Gram-positive, catalase-negative organisms, which occupy a wide range of natural plant- and animal-associated environments. LAB species are frequently involved in the transformation of perishable raw materials into more stable, pleasant, palatable and safe fermented food products. LAB and bifidobacteria are also found among the resident microbiota of the gastrointestinal and/or genitourinary tracts of vertebrates, where they are believed to exert health-promoting effects. At present, the genomes of more than 20 LAB and bifidobacterial species have been completely sequenced. Their genome content reflects its specific metabolism, physiology, biosynthetic capabilities, and adaptability to varying conditions and environments. The typical LAB/bifidobacterial genome is relatively small (from 1.7 to 3.3 Mb) and thus harbors a limited assortment of genes (from around 1,600 to over 3,000). These small genomes code for a broad array of transporters for efficient carbon and nitrogen assimilation from the nutritionally-rich niches they usually inhabit, and specify a rather limited range of biosynthetic and degrading capabilities. The variation in the number of genes suggests that the genome evolution of each of these bacterial groups involved the processes of extensive gene loss from their particular ancestor, diversification of certain common biological activities through gene duplication, and acquisition of key functions via horizontal gene transfer. The availability of genome sequences is expected to revolutionize the exploitation of the metabolic potential of LAB and bifidobacteria, improving their use in bioprocessing and their utilization in biotechnological and health-related applications.

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.

Regulation of flpA, flpB and rcfA promoters in Lactococcus lactis

E. coli fumarate nitrate reductase (FNR) binds to conserved FNR sites to regulate transcription under anaerobic condition. L. lactis subsp. cremoris MG1363 strain contains two FNR-like proteins (FlpA and B) encoded by flpA and flpB genes and the rcfA gene-encoded RcfA in L. lactis subsp. lactis IL1403 strain. Potential FNR-binding sites were located upstream of these genes. The flpA promoter is expressed in MG1363 anaerobically and aerobically. The flpB and rcfA promoters have typical class II FNR-dependent promoters and are activated anaerobically in MG1363 and IL1403, respectively. Despite their strong homology, the Flp and RcfA proteins cannot substitute for each other and control these promoters in the heterologous strains. The flpA and flpB promoters require FlpA and FlpB for activation in the MG1363 background. This was confirmed by expressing FlpB under nisin control in flp mutants and monitoring flpA promoter expression. In flpB- backgrounds, both FlpA and FlpB were required for flpA promoter expression. FlpB could not complement for the lack of FlpA protein in flpA- backgrounds.

Large chromosomal inversion correlated with spectinomycin resistance in Lactococcus lactis subsp. lactis bv. diacetylactis S50

A large chromosomal inversion that confers resistance to high concentrations of the antibiotic spectinomycin in Lactococcus lactis subsp. lactis bv. diacetylactis S50 was identified by pulsed field gel electrophoresis. The same type of inversion was identified in 4 independent experiments and in 4 different derivatives of strain S50, indicating the same position and the same mechanism of recombination as a response to antibiotic selective pressure in all derivatives. An analysis of ribosomal operons in strain S50 and mutants revealed that ribosomal operons are not endpoints of the recombination. Spectinomycin-resistant mutants appeared in a population of S50 derivatives at a high frequency of 2 × 10−7. These spectinomycin-resistant mutants were not able to compete successfully with the wild-type strain during 25 generations (48 h) of co-culture in vitro, indicating that inversion had a significant fitness cost. Results demonstrate that as a mechanism of genome plasticity, inversion can be directly involved in one-step development of the adaptation to a high concentration of spectinomycin.

Single independent operator sites are involved in the genetic switch of the Lactobacillus delbrueckii bacteriophage mv4

The lysogeny region of the Lactobacillus delbrueckii bacteriophage mv4 contains two divergently oriented ORFs coding for the Rep (221 aa) and Tec (64 aa) proteins. The transcription of these two genes was analysed by primer extension and Northern blot experiments on lysogenic strains. The location of the transcription initiation sites of rep and tec in the intergenic region allowed the identification of the divergently oriented non overlapping promoters P(rep) and P(tec). Transcriptional fusions analysis showed that Rep negatively regulates the P(tec) promoter and activates its own transcription, and that Tec is a negative regulator of the two promoters. As demonstrated by gel mobility shift assays, the repressor Rep binds to a single specific 17 bp site located between the P(tec) -10 and -35 regions whereas Tec binds to a single specific 40 bp long complex operator site located between the two promoters. The presence of a single specific operator site for each repressor in the intergenic region is an unusual feature.

The unconventional Xer recombination machinery of Streptococci/Lactococci

Homologous recombination between circular sister chromosomes during DNA replication in bacteria can generate chromosome dimers that must be resolved into monomers prior to cell division. In Escherichia coli, dimer resolution is achieved by site-specific recombination, Xer recombination, involving two paralogous tyrosine recombinases, XerC and XerD, and a 28-bp recombination site (dif) located at the junction of the two replication arms. Xer recombination is tightly controlled by the septal protein FtsK. XerCD recombinases and FtsK are found on most sequenced eubacterial genomes, suggesting that the Xer recombination system as described in E. coli is highly conserved among prokaryotes. We show here that Streptococci and Lactococci carry an alternative Xer recombination machinery, organized in a single recombination module. This corresponds to an atypical 31-bp recombination site (dif(SL)) associated with a dedicated tyrosine recombinase (XerS). In contrast to the E. coli Xer system, only a single recombinase is required to recombine dif(SL), suggesting a different mechanism in the recombination process. Despite this important difference, XerS can only perform efficient recombination when dif(SL) sites are located on chromosome dimers. Moreover, the XerS/dif(SL) recombination requires the streptococcal protein FtsK(SL), probably without the need for direct protein-protein interaction, which we demonstrated to be located at the division septum of Lactococcus lactis. Acquisition of the XerS recombination module can be considered as a landmark of the separation of Streptococci/Lactococci from other firmicutes and support the view that Xer recombination is a conserved cellular function in bacteria, but that can be achieved by functional analogs.

Complete genome sequence of the prototype lactic acid bacterium Lactococcus lactis subsp. cremoris MG1363

Lactococcus lactis is of great importance for the nutrition of hundreds of millions of people worldwide. This paper describes the genome sequence of Lactococcus lactis subsp. cremoris MG1363, the lactococcal strain most intensively studied throughout the world. The 2,529,478-bp genome contains 81 pseudogenes and encodes 2,436 proteins. Of the 530 unique proteins, 47 belong to the COG (clusters of orthologous groups) functional category "carbohydrate metabolism and transport," by far the largest category of novel proteins in comparison with L. lactis subsp. lactis IL1403. Nearly one-fifth of the 71 insertion elements are concentrated in a specific 56-kb region. This integration hot-spot region carries genes that are typically associated with lactococcal plasmids and a repeat sequence specifically found on plasmids and in the "lateral gene transfer hot spot" in the genome of Streptococcus thermophilus. Although the parent of L. lactis MG1363 was used to demonstrate lysogeny in Lactococcus, L. lactis MG1363 carries four remnant/satellite phages and two apparently complete prophages. The availability of the L. lactis MG1363 genome sequence will reinforce its status as the prototype among lactic acid bacteria through facilitation of further applied and fundamental research.

Localization, mobility and fidelity of retrotransposed Group II introns in rRNA genes

We previously showed that the group II Lactococcus lactis Ll.LtrB intron could retrotranspose into ectopic locations on the genome of its native host. Two integration events, which had been mapped to unique sequences, were localized in the present study to separate copies of the six L.lactis 23S rRNA genes, within operon B or D. Although further movement within the bacterial chromosome was undetectable, the retrotransposed introns were able to re-integrate into their original homing site provided on a plasmid. This finding indicates not only that retrotransposed group II introns retain mobility properties, but also that movement occurs back into sequence that is heterologous to the sequence of the chromosomal location. Sequence analysis of the retrotransposed introns and the secondary mobility events back to the homing site showed that the introns retain sequence integrity. These results are illuminating, since the reverse transcriptase (RT) of the intron-encoded protein, LtrA, has no known proofreading function, yet the mobility events have a low error rate. Enzymatic digests were used to monitor sequence changes from the wild-type intron. The results indicate that retromobility events have approximately 10(-5) misincorporations per nucleotide inserted. In contrast to the high RT error rates for retroviruses that must escape host defenses, the infrequent mutations of group II introns would ensure intron spread through retention of sequences essential for mobility.

Projector 2: contig mapping for efficient gap-closure of prokaryotic genome sequence assemblies

With genome sequencing efforts increasing exponentially, valuable information accumulates on genomic content of the various organisms sequenced. Projector 2 uses (un)finished genomic sequences of an organism as a template to infer linkage information for a genome sequence assembly of a related organism being sequenced. The remaining gaps between contigs for which no linkage information is present can subsequently be closed with direct PCR strategies. Compared with other implementations, Projector 2 has several distinctive features: a user-friendly web interface, automatic removal of repetitive elements (repeat-masking) and automated primer design for gap-closure purposes. Moreover, when using multiple fragments of a template genome, primers for multiplex PCR strategies can also be designed. Primer design takes into account that, in many cases, contig ends contain unreliable DNA sequences and repetitive sequences. Closing the remaining gaps in prokaryotic genome sequence assemblies is thereby made very efficient and virtually effortless. We demonstrate that the use of single or multiple fragments of a template genome (i.e. unfinished genome sequences) in combination with repeat-masking results in mapping success rates close to 100%. The web interface is freely accessible at http://molgen.biol.rug.nl/websoftware/projector2.

Chromosomal constraints in Gram-positive bacteria revealed by artificial inversions

We used artificial chromosome inversions to investigate the chromosomal constraints that preserve genome organization in the Gram-positive bacterium Lactococcus lactis. Large inversions, 80-1260 kb in length, disturbing the symmetry of the origin and terminus of the replication axis to various extents, were constructed using the site-specific Cre-loxP recombination system. These inversions were all mechanistically feasible and fell into various classes according to stability and effect on cell fitness. The L. lactis chromosome supports only to some extent unbalance in length of its replication arms. The location of detrimental inversions allowed identification of two constrained chromosomal regions: a large domain covering one fifth of the genome that encompasses the origin of replication (Ori domain), and a smaller domain located at the opposite of the chromosome (Ter domain).

A bacterial group II intron favors retrotransposition into plasmid targets

Group II introns, widely believed to be the ancestors of nuclear pre-mRNA introns, are catalytic RNAs found in bacteria, archaea, and eukaryotes. They are mobile genetic elements that move via an RNA intermediate. They retrohome to intronless alleles and retrotranspose to ectopic sites, aided by an intron-encoded protein with reverse transcriptase, maturase, and endonuclease activities. Many group II introns identified in bacteria reside on plasmid genomes rather than bacterial chromosomes, implying that plasmids are havens for these retroelements. This study demonstrates that almost one-fourth of retrotransposition events of the Ll.LtrB intron in Lactococcus lactis are into the plasmid donor. This level is more than twice that predicted based on target size and plasmid copy number relative to the chromosome. In particular, the fraction of such plasmid targeting events was elevated to more than one-third of retrotransposition events by mutation of the intron-encoded endonuclease, a situation that may resemble most bacterial group II introns, which lack the endonuclease. Target-site sequences on the plasmid are more relaxed than those on the chromosome, likely accounting for preferred integration into plasmid replicons. Furthermore, the direction of integration relative to promoters and origins of replication is consistent with group II intron retrotransposition into single-stranded DNA at replication forks. This work provides mechanistic rationales for the prevalence of group II introns in natural plasmid populations and underscores that targeting to plasmids, which are themselves mobile elements, could promote intron spread.

Increased production of folate by metabolic engineering of Lactococcus lactis

The dairy starter bacterium Lactococcus lactis is able to synthesize folate and accumulates large amounts of folate, predominantly in the polyglutamyl form. Only small amounts of the produced folate are released in the extracellular medium. Five genes involved in folate biosynthesis were identified in a folate gene cluster in L. lactis MG1363: folA, folB, folKE, folP, and folC. The gene folKE encodes the biprotein 2-amino-4-hydroxy-6-hydroxymethyldihydropteridine pyrophosphokinase and GTP cyclohydrolase I. The overexpression of folKE in L. lactis was found to increase the extracellular folate production almost 10-fold, while the total folate production increased almost 3-fold. The controlled combined overexpression of folKE and folC, encoding polyglutamyl folate synthetase, increased the retention of folate in the cell. The cloning and overexpression of folA, encoding dihydrofolate reductase, decreased the folate production twofold, suggesting a feedback inhibition of reduced folates on folate biosynthesis.

Retrotransposition of the Ll.LtrB group II intron proceeds predominantly via reverse splicing into DNA targets

Catalytic group II introns are mobile retroelements that invade cognate intronless genes via retrohoming, where the introns reverse splice into double-stranded DNA (dsDNA) targets. They can also retrotranspose to ectopic sites at low frequencies. Whereas our previous studies with a bacterial intron, Ll.LtrB, supported frequent use of RNA targets during retrotransposition, recent experiments with a retrotransposition indicator gene indicate that DNA, rather than RNA, is a prominent target, with both dsDNA and single-stranded DNA (ssDNA) as possibilities. Thus retrotransposition occurs in both transcriptional sense and antisense orientations of target genes, and is largely independent of homologous DNA recombination and of the endonuclease function of the intron-encoded protein, LtrA. Models based on both dsDNA and ssDNA targeting are presented. Interestingly, retrotransposition is biased toward the template for lagging-strand DNA synthesis, which suggests the possibility of the replication folk as a source of ssDNA. Consistent with some use of ssDNA targets, many retrotransposition sites lack nucleotides critical for the unwinding of target duplex DNA. Moreover, in vitro the intron reverse spliced into ssDNA more efficiently than dsDNA substrates for some of the retrotransposition sites. Furthermore, many bacterial group II introns reside on the lagging-strand template, hinting at a role for DNA replication in intron dispersal in nature.

Lactococcus lactis DPC5598, a plasmid-free derivative of a commercial starter, provides a valuable alternative host for culture improvement studies

AIMS

To generate a plasmid-free derivative of an extensively used industrial starter strain Lactococcus lactis DPC4268, which could be used as a backbone strain for starter improvement programmes.

METHODS AND RESULTS

DPC4268 containing four large plasmids was subjected to high temperature plasmid curing resulting in derivatives, each with a different plasmid complement of one, two or three different plasmids in addition to a plasmid-free derivative. Industrially relevant phenotypes were assigned to each plasmid on the basis of detailed phenotypic and genetic analyses and these were (a) proteinase activity (Prt, 60 kb) (b) lactose fermentation (Lac, 55 kb) (c) bacteriophage adsorption inhibition (Ads, 44 kb) and (d) type I restriction/modification (R/M, 40 kb). The plasmid-free variant of DPC4268 was shown to be transformable at frequencies comparable to the common laboratory strain L. lactis MG1614. Furthermore its genome was demonstrated to be significantly different from the laboratory strains L. lactis MG1614 and the recently sequenced L. lactis IL1403 genomes by pulsed-field gel electrophoresis.

CONCLUSIONS

This study produced an easily transformable plasmid-free derivative which was genomically different from both MG1614 and IL1403. In addition, important plasmid-borne industrial traits, including two phage-resistance mechanisms, were identified in DPC4268.

SIGNIFICANCE AND IMPACT OF THE STUDY

L. DPC4268 is a vitally important commercial strain used in the manufacture of Cheddar cheese. The generation of a plasmid-free derivative may provide an important backbone strain as a basis for future strain improvement purposes.

Discovering lactic acid bacteria by genomics

This review summarizes a collection of lactic acid bacteria that are now undergoing genomic sequencing and analysis. Summaries are presented on twenty different species, with each overview discussing the organisms fundamental and practical significance, environmental habitat, and its role in fermentation, bioprocessing, or probiotics. For those projects where genome sequence data were available by March 2002, summaries include a listing of key statistics and interesting genomic features. These efforts will revolutionize our molecular view of Gram-positive bacteria, as up to 15 genomes from the low GC content lactic acid bacteria are expected to be available in the public domain by the end of 2003. Our collective view of the lactic acid bacteria will be fundamentally changed as we rediscover the relationships and capabilities of these organisms through genomics.

Cre-loxP recombination system for large genome rearrangements in Lactococcus lactis

We have used a new genetic strategy based on the Cre-loxP recombination system to generate large chromosomal rearrangements in Lactococcus lactis. Two loxP sites were sequentially integrated in inverse order into the chromosome either at random locations by transposition or at fixed points by homologous recombination. The recombination between the two chromosomal loxP sites was highly efficient (approximately 1 x 10(-1)/cell) when the Cre recombinase was provided in trans, and parental- or inverted-type chromosomal structures were isolated after removal of the Cre recombinase. The usefulness of this approach was demonstrated by creating three large inversions of 500, 1,115, and 1,160 kb in size that modified the lactococcal genome organization to different extents. The Cre-loxP recombination system described can potentially be used for other gram-positive bacteria without further modification.

The complete genome sequence of the lactic acid bacterium Lactococcus lactis ssp. lactis IL1403

Lactococcus lactis is a nonpathogenic AT-rich gram-positive bacterium closely related to the genus Streptococcus and is the most commonly used cheese starter. It is also the best-characterized lactic acid bacterium. We sequenced the genome of the laboratory strain IL1403, using a novel two-step strategy that comprises diagnostic sequencing of the entire genome and a shotgun polishing step. The genome contains 2,365,589 base pairs and encodes 2310 proteins, including 293 protein-coding genes belonging to six prophages and 43 insertion sequence (IS) elements. Nonrandom distribution of IS elements indicates that the chromosome of the sequenced strain may be a product of recent recombination between two closely related genomes. A complete set of late competence genes is present, indicating the ability of L. lactis to undergo DNA transformation. Genomic sequence revealed new possibilities for fermentation pathways and for aerobic respiration. It also indicated a horizontal transfer of genetic information from Lactococcus to gram-negative enteric bacteria of Salmonella-Escherichia group.

The Complete Genome Sequence of the Lactic Acid Bacterium Lactococcus lactis ssp. lactis IL1403

Lactococcus lactis is a nonpathogenic AT-rich gram-positive bacterium closely related to the genus Streptococcus and is the most commonly used cheese starter. It is also the best-characterized lactic acid bacterium. We sequenced the genome of the laboratory strain IL1403, using a novel two-step strategy that comprises diagnostic sequencing of the entire genome and a shotgun polishing step. The genome contains 2,365,589 base pairs and encodes 2310 proteins, including 293 protein-coding genes belonging to six prophages and 43 insertion sequence (IS) elements. Nonrandom distribution of IS elements indicates that the chromosome of the sequenced strain may be a product of recent recombination between two closely related genomes. A complete set of late competence genes is present, indicating the ability of L. lactis to undergo DNA transformation. Genomic sequence revealed new possibilities for fermentation pathways and for aerobic respiration. It also indicated a horizontal transfer of genetic information fromLactococcus to gram-negative enteric bacteria ofSalmonella-Escherichia group.

[The sequence data described in this paper has been submitted to the GenBank data library under accession no. AE005176.]

Distinctive features of homologous recombination in an 'old' microorganism, Lactococcus lactis

Homologous recombination is needed to assure faithful inheritance of DNA material, especially under stress conditions. The same enzymes that repair broken chromosomes via recombination also generate biodiversity. Their activities may result in intrachromosomal rearrangements, assimilation of foreign DNA, or a combination of these events. It is generally supposed that homologous recombination systems are conserved, and function the same way everywhere as they do in Escherichia coli, the accepted paradigm. Studies in an 'older' microorganism, the gram-positive bacterium of the low GC branch Lactococcus lactis, confirm that many enzymes are conserved across species lines. However, the main components of the double strand break (DSB) repair system, an exonuclease/helicase (Exo/hel) and a short DNA modulator sequence Chi, differ markedly between bacteria, especially when compared to the gram-negative analogues. Based on our studies, a model is proposed for the functioning of the two-subunit Exo/hel of L. lactis and other gram-positive bacteria, which differs from that of the three-subunit E. coli enzyme. The differences between bacterial DSB repair systems may underlie a selection for diversity when dealing with DSB. These and other features of homologous recombination in L. lactis are discussed.

Novel sucrose transposons from plant strains of Lactococcus lactis

Lactococcus lactis strains isolated from vegetable products transferred the ability to ferment sucrose in conjugation experiments with the recipient strain L. lactis MG1614. Nisin production and sucrose fermentation were transferred together from two strains, but transfer also occurred from several other strains which did not produce nisin. Pulsed-field gel electrophoresis analysis showed that all transconjugants had acquired large chromosomal insertions at two main sites. Nisin sucrose transconjugants had gained inserts of 70 kb, while those that fermented sucrose without nisin production contained inserts of between 50 and 110 kb. Transconjugants from one donor had acquired a separate insertion of 55 kb which correlated with enhanced bacteriophage resistance, but contained neither nisin nor sucrose fermentation genes.

Short communication: salt extends the upper temperature limit for growth of Lactococcus lactis ssp. cremoris on solid M17 medium

We have determined conditions for plating of the Lactococcus lactis ssp. cremoris laboratory strain MG1363 on solid M17 broth at 38 degrees C, which is required for the optimal use of the pGhost plasmids. The addition of 1% NaCl (or KCl, potassium acetate, or sucrose at 170 mM) to M17 agar plates results in extension of the upper temperature limit for growth from 37 to 40 degrees C; no decrease in plating efficiency was detected from 30 to 39 degrees C.

Genome plasticity among related ++Lactococcus strains: identification of genetic events associated with macrorestriction polymorphisms

The genomic diversity of nine strains of the Lactococcus lactis subsp. cremoris (NCDO712, NCDO505, NCDO2031, NCDO763, MMS36, C2, LM0230, LM2301, and MG1363) was studied by macrorestriction enzyme analysis using pulsed-field gel electrophoresis. These strains were considered adequate for the investigation of genomic plasticity because they have been described as belonging to the same genetic lineage. Comparison of ApaI and SmaI genome fingerprints of each strain revealed the presence of several macrorestriction fragment length polymorphisms (RFLPs), despite a high degree of similarity of the generated restriction patterns. The physical map of the MG1363 chromosome was used to establish a genome map of the other strains and allocate the RFLPs to five regions. Southern hybridization analysis correlated the polymorphic regions with genetic events such as chromosomal inversion, integration of prophage DNA, and location of the transposon-like structures carrying conjugative factor or oligopeptide transport system.

Physical and genetic map of Streptococcus mutans GS-5 and localization of five rRNA operons

The physical map of the 2.1 megabase chromosome of Streptococcus mutans GS-5 has been refined by including all ApaI and SmaI fragments of 5 kbp or greater, and by positioning the fragments generated by the endonuclease I-CeuI. Sixty-three new genetic loci have been added to the map, so that it now contains 90 loci. The new loci include those for 35 cloned streptococcal genes of established function and for 23 S. mutans genes of putative function. In addition, five rrn operons were identified and placed on the map of the chromosome. The presence of a SmaI site in each of the rrn operons allowed the direction of transcription of each operon to be deduced. The orientation of the rrn loci indicates that their transcription is directed away from a small region of the chromosome, identifying a possible region for the initiation of chromosome replication.

The diverse and dynamic structure of bacterial genomes

Bacterial genome sizes, which range from 500 to 10,000 kbp, are within the current scope of operation of large-scale nucleotide sequence determination facilities. To date, 8 complete bacterial genomes have been sequenced, and at least 40 more will be completed in the near future. Such projects give wonderfully detailed information concerning the structure of the organism's genes and the overall organization of the sequenced genomes. It will be very important to put this incredible wealth of detail into a larger biological picture: How does this information apply to the genomes of related genera, related species, or even other individuals from the same species? Recent advances in pulsed-field gel electrophoretic technology have facilitated the construction of complete and accurate physical maps of bacterial chromosomes, and the many maps constructed in the past decade have revealed unexpected and substantial differences in genome size and organization even among closely related bacteria. This review focuses on this recently appreciated plasticity in structure of bacterial genomes, and diversity in genome size, replicon geometry, and chromosome number are discussed at inter- and intraspecies levels.

Characterization of multiple regions involved in replication and mobilization of plasmid pNZ4000 coding for exopolysaccharide production in Lactococcus lactis

We characterized the regions involved in replication and mobilization of the 40-kb plasmid pNZ4000, encoding exopolysaccharide (EPS) production in Lactococcus lactis NIZO B40. The plasmid contains four highly conserved replication regions with homologous rep genes (repB1, repB2, repB3, and repB4) that belong to the lactococcal theta replicon family. Subcloning of each replicon individually showed that all are functional and compatible in L. lactis. Plasmid pNZ4000 and genetically labeled derivatives could be transferred to different L. lactis strains by conjugation, and pNZ4000 was shown to be a mobilization plasmid. Two regions involved in mobilization were identified near two of the replicons; both included an oriT sequence rich in inverted repeats. Conjugative mobilization of the nonmobilizable plasmid pNZ124 was promoted by either one of these oriT sequences, demonstrating their functionality. One oriT sequence was followed by a mobA gene, coding for a trans-acting protein, which increased the frequency of conjugative transfer 100-fold. The predicted MobA protein and the oriT sequences show protein and nucleotide similarity, respectively, with the relaxase and with the inverted repeat and nic site of the oriT from the Escherichia coli plasmid R64. The presence on pNZ4000 of four functional replicons, two oriT sequences, and several insertion sequence-like elements strongly suggests that this EPS plasmid is a naturally occurring cointegrate.

A natural large chromosomal inversion in Lactococcus lactis is mediated by homologous recombination between two insertion sequences

Comparative analysis of chromosomal macrorestriction polymorphism of the two closely related Lactococcus lactis subsp. cremoris strains MG1363 and NCDO763 revealed the presence of a large inversion covering half of the genome. To determine what kind of genetic element could be implicated in this rearrangement, the two inversion junctions of MG1363 and NCDO763 chromosomes were cloned and characterized. Nucleotide sequence analysis showed the presence of one copy of the lactococcal IS905 element in each junction. Each copy of this element contained the same nucleotide mutation that inactivates the putative transposase. Comparison of the sequences surrounding the insertion sequence demonstrated that the large inversion arose from a single-step homologous recombination event between the two defective copies of the IS905 element. The large inversion presumably conferred no selective disadvantage on strain NCDO763 because this rearrangement did not alter the oriC-terC symmetry of the chromosome and the local genetic environment.

Combinational variation of restriction modification specificities in Lactococcus lactis

Three genes coding for a type I R-M system related to the class C enzymes have been identified on the chromosome of Lactococcus lactis strain IL1403. In addition, plasmids were found that encode only the HsdS subunit that directs R-M specificity. The presence of these plasmids in IL1403 conferred a new R-M phenotype on the host, indicating that the plasmid-encoded HsdS is able to interact with the chromosomally encoded HsdR and HsdM subunits. Such combinational variation of type I R-M systems may facilitate the evolution of their specificity and thus reinforce bacterial resistance against invasive foreign unmethylated DNA.

The citrate transport system of Lactococcus lactis subsp. lactis biovar diacetylactis is induced by acid stress

Citrate transport in Lactococcus lactis subsp. lactis biovar diacetylactis is catalyzed by citrate permease P (CitP), which is encoded by the plasmidic citP gene. We have shown previously that citP is included in the citQRP operon, which is mainly transcribed from the P1 promoter in L. lactis subsp. lactis biovar diacetylactis. furthermore, transcription of citQRP and citrate transport are not induced by the presence of citrate in the growth medium. In this work, we analyzed the influence of the extracellular pH on the expression of citP. The citrate transport system is induced by natural acidification of the medium during cell growth and by a shift to media buffered at acidic pHs. This inducible response to acid stress takes place at the transcriptional level and seems to be due to increased utilization of the P1 promoter. Increased transcription correlates with increased synthesis of CitP and results in higher citrate transport activity catalyzed by the cells. Finally, this acid stress response seems to provide L. lactis subsp. lactis biovar diacetylactis with a selective advantage resulting from cometabolism of glucose and citrate at low pHs.

A chloride-inducible acid resistance mechanism in Lactococcus lactis and its regulation

Previously, a promoter was identified in Lactococcus lactis that is specifically induced by chloride. Here, we describe the nucleotide sequence and functional analysis of two genes transcribed from this promoter, gadC and gadB. GadC is homologous to putative glutamate-gamma-aminobutyrate antiporters of Escherichia coli and Shigella flexneri and contains 12 putative membrane-spanning domains. GadB shows similarity to glutamate decarboxylases. A L. lactis gadB mutant and a strain that is unable to express both gadB and gadC was more sensitive to low pH than the wild type when NaCl and glutamate were present. Expression of gadCB in L. lactis in the presence of chloride was increased when the culture pH was allowed to decrease to low levels by omitting buffer from the medium, while glutamate also stimulated gadCB expression. Apparently, these genes encode a glutamate-dependent acid resistance mechanism of L. lactis that is optimally active under conditions in which it is needed to maintain viability. Immediately upstream of the chloride-dependent gadCB promoter Pgad, a third gene encodes a protein (GadR) that is homologous to the activator Rgg from Streptococcus gordonii. gadR expression is chloride and glutamate independent. A gadR mutant did not produce the 3kb gadCB mRNA that is found in wild-type cells in the presence of NaCl, indicating that GadR is an activator of the gadCB operon.

Large genome rearrangements discovered by the detailed analysis of 21 Pseudomonas aeruginosa clone C isolates found in environment and disease habitats

In order to determine primary genetic events which occur during the diversification of a Pseudomonas aeruginosa clone in natural habitats, comparative genome analysis of 21 isolates of a predominant clone, called clone C, derived mainly from patients with cystic fibrosis (CF) and the aquatic environment, was carried out. Physical chromosome maps were constructed for the restriction enzymes SpeI, PacI, SwaI and I-CeuI by one and two-dimensional pulsed-field gel electrophoresis and by comparison with the existing strain C map. The positioning of 26 genes generated the genetic maps. Chromosome size varied between 6345 and 6606 kilobase-pairs (kb). A plasmid of 95 kb was detected in the strains of non-CF origin and, in addition, was found to be integrated into the chromosome of all strains but one CF isolate. Four subgroups of clone C strains were discriminated by the acquisition and loss of large blocks of DNA that could cover more than 10% of the chromosome size. The exchange of DNA blocks which ranged in size from 1 kb to 214 kb occurred preferentially around the terminus of replication region which is poor in biosynthetic genes. Genetic material which was additionally introduced into strain C in comparison with strain PAO seems to be a target of mutational processes in clone C strains. Within and among subgroups CF isolates frequently exhibited large inversions affecting the whole chromosomal structure. We concluded that the exchange of DNA blocks by mechanisms of horizontal transfer and large chromosomal inversions are major factors leading to the divergence of a clone in the species P. aeruginosa.

Genome structure and phylogeny in the genus Brucella

PacI and SpeI restriction maps were obtained for the two chromosomes of each of the six species of the genus Brucella: B. melitensis, B. abortus, B. suis, B. canis, B. ovis, and B. neotomae. Three complementary techniques were used: hybridization with the two replicons as probes, cross-hybridization of restriction fragments, and a new mapping method. For each type strain, a unique I-SceI site was introduced in each of the two replicons, and the location of SpeI sites was determined by linearization at the unique site, partial digestion, and end labeling of the fragments. The restriction and genetic maps of the six species were highly conserved. However, numerous small insertions or deletions, ranging from 1 to 34 kb, were observed by comparison with the map of the reference strain of the genus, B. melitensis 16M. A 21-kb Spel fragment specific to B. ovis was found in the small chromosome of this species. A 640-kb inversion was demonstrated in the B. abortus small chromosome. All of these data allowed the construction of a phylogenetic tree, which reflects the traditional phenetic classification of the genus.

Primary and secondary structures of rRNA spacer regions in enterococci

The 16S-23S and 23S-5S rRNA spacer DNA regions (spacer regions 1 and 2, respectively) from Enterococcus faecalis, Enterococcus faecium, Enterococcus hirae, Enterococcus durans and Enterococcus mundtii were amplified by PCR. Their nucleotide sequences were established and a secondary structure model showing the interaction between the two spacer regions was built. Whereas lactococci and Streptococcus sensu stricto are characterized by a single type of spacer region 1, the enterococci show a high degree of variability in this region; thus the spacer regions 1 with and without tRNA(Ala) were characterized. However, as shown for lactococci and Streptococcus sensu stricto, the tRNA(Ala) gene does not encode the 3'-terminal CCA trinucleotide. A putative antitermination signal is found downstream from the tRNA(Ala) gene. Based on comparison with Lactococcus lactis and Streptococcus thermophilus, a double-stranded processing stem is proposed. In E, hirae, one of the three different types of spacer region 1 contains no tRNA(Ala), but displays a 107 nt insertion that forms a long stem-loop structure. A similar insertion (115 nt in length) was found in E. faecium and base compensatory mutations preserve the ability to form the long stem-loop structure. Such insertions may correspond to mobile intervening sequences, as found in the 23S rRNA coding sequences of some Gram-negative bacteria. The spacer regions 1 and 2 from the three subgroups of streptococci were compared, and except for the tRNA(Ala) gene and the double-stranded processing sites, little similarity was found, which opens large possibilities for future development of DNA-based typing methods.

Plasmid integration in a wide range of bacteria mediated by the integrase of Lactobacillus delbrueckii bacteriophage mv4

Bacteriophage mv4 is a temperate phage infecting Lactobacillus delbrueckii subsp. bulgaricus. During lysogenization, the phage integrates its genome into the host chromosome at the 3' end of a tRNA(Ser) gene through a site-specific recombination process (L. Dupont et al., J. Bacteriol., 177:586-595, 1995). A nonreplicative vector (pMC1) based on the mv4 integrative elements (attP site and integrase-coding int gene) is able to integrate into the chromosome of a wide range of bacterial hosts, including Lactobacillus plantarum, Lactobacillus casei (two strains), Lactococcus lactis subsp. cremoris, Enterococcus faecalis, and Streptococcus pneumoniae. Integrative recombination of pMC1 into the chromosomes of all of these species is dependent on the int gene product and occurs specifically at the pMC1 attP site. The isolation and sequencing of pMC1 integration sites from these bacteria showed that in lactobacilli, pMC1 integrated into the conserved tRNA(Ser) gene. In the other bacterial species where this tRNA gene is less or not conserved; secondary integration sites either in potential protein-coding regions or in intergenic DNA were used. A consensus sequence was deduced from the analysis of the different integration sites. The comparison of these sequences demonstrated the flexibility of the integrase for the bacterial integration site and suggested the importance of the trinucleotide CCT at the 5' end of the core in the strand exchange reaction.

The ldh phylogeny for environmental isolates of Lactococcus lactis is consistent with rRNA genotypes but not with phenotypes

Lactate dehydrogenase (ldh) gene sequences, levels of 16S rRNA group-specific probe binding, and phenotypic characteristics were compared for 45 environmental isolates and four commercial starter strains of Lactococcus lactis to identify evolutionary groups best suited to cheddar cheese manufacture, ldh sequences from the environmental isolates showed high similarity to those from two groups of L. lactis used for industrial fermentations, L. lactis subsp. cremoris and subsp. lactis. Within each phylogenetically defined subspecies, ldh sequence similarities were greater than 99.1%. Strains with phenotypic traits formerly diagnostic for both subspecies were found in each ldh similarity group, but only strains belonging to L. lactis subsp. cremoris by both the newer, genetic and the older, superseded phenotypic criteria were judged potentially suitable for the commercial production of cheddar cheese. Identical evolutionary relationships were inferred from ldh sequences and from binding of subspecies-specific, 16S rRNA-directed oligonucleotide probes. However, groups defined according to these chromosomal traits bore no relationship to patterns of arginine deamination, carbon substrate utilization, or bacteriophage sensitivity, which may be encoded by cryptic genes or sexually transmissible genetic elements. Fourteen new L. lactis subsp. cremoris isolates were identified as suitable candidates for cheddar cheese manufacture, and 10 of these were completely resistant to three different batteries of commercial bacteriophages known to reduce starter activity.