Attachment of the RNA degradosome to the bacterial inner cytoplasmic membrane prevents wasteful degradation of rRNA in ribosome assembly intermediates.

RNA processing and degradation shape the transcriptome by generating stable molecules that are necessary for translation (rRNA and tRNA) and by facilitating the turnover of mRNA, which is necessary for the posttranscriptional control of gene expression. In bacteria and the plant chloroplast, RNA degradosomes are multienzyme complexes that process and degrade RNA. In many bacterial species, the endoribonuclease RNase E is the central component of the RNA degradosome. RNase E-based RNA degradosomes are inner membrane proteins in a large family of gram-negative bacteria (ß- and γ-Proteobacteria). Until now, the reason for membrane localization was not understood. Here, we show that a mutant strain of Escherichia coli, in which the RNA degradosome is localized to the interior of the cell, has high levels of 20S and 40S particles that are defective intermediates in ribosome assembly. These particles have aberrant protein composition and contain rRNA precursors that have been cleaved by RNase E. After RNase E cleavage, rRNA fragments are degraded to nucleotides by exoribonucleases. In vitro, rRNA in intact ribosomes is resistant to RNase E cleavage, whereas protein-free rRNA is readily degraded. We conclude that RNA degradosomes in the nucleoid of the mutant strain interfere with cotranscriptional ribosome assembly. We propose that membrane-attached RNA degradosomes in wild-type cells control the quality of ribosome assembly after intermediates are released from the nucleoid. That is, the compact structure of mature ribosomes protects rRNA against cleavage by RNase E. Turnover of a proportion of intermediates in ribosome assembly explains slow growth of the mutant strain. Competition between mRNA and rRNA degradation could be the cause of slower mRNA degradation in the mutant strain. We conclude that attachment of the RNA degradosome to the bacterial inner cytoplasmic membrane prevents wasteful degradation of rRNA precursors, thus explaining the reason for conservation of membrane-attached RNA degradosomes throughout the ß- and γ-Proteobacteria.

Redefining the bacteriophage mv4 site-specific recombination system and the sequence specificity of its attB and core-attP sites.

Through their involvement in the integration and excision of a large number of mobile genetic elements, such as phages and integrative and conjugative elements (ICEs), site-specific recombination systems based on heterobivalent tyrosine recombinases play a major role in genome dynamics and evolution. However, despite hundreds of these systems having been identified in genome databases, very few have been described in detail, with none from phages that infect Bacillota (formerly Firmicutes). In this study, we reanalyzed the recombination module of Lactobacillus delbrueckii subsp. bulgaricus phage mv4, previously considered atypical compared with classical systems. Our results reveal that mv4 integrase is a 369 aa protein with all the structural hallmarks of recombinases from the Tn916 family and that it cooperatively interacts with its recombination sites. Using randomized DNA libraries, NGS sequencing, and other molecular approaches, we show that the 21-bp core-attP and attB sites have structural similarities to classical systems only if considering the nucleotide degeneracy, with two 7-bp inverted regions corresponding to mv4Int core-binding sites surrounding a 7-bp strand-exchange region. We also examined the different compositional constraints in the core-binding regions, which define the sequence space of permissible recombination sites.

When translation elongation is impaired, the mRNA is uniformly destabilized by the RNA degradosome, while the concentration of mRNA is altered along the molecule.

mRNA sits at the crossroads of transcription, translation and mRNA degradation. Many questions remain about the coupling of these three processes in Escherichia coli and, in particular, how translation may have an effect on mRNA degradation and transcription. To characterize the interplay between mRNA degradation and translation while accounting for transcription, we altered the translation initiation or elongation and measured the effects on mRNA stability and concentration. Using a mapping method, we analysed mRNA concentration and stability at the local scale all along the transcript. We showed that a decrease in translation initiation efficiency destabilizes the mRNA and leads to a uniform decrease in mRNA concentration throughout the molecule. Prematurely terminating translation elongation by inserting a stop codon is associated with a drop in local mRNA concentration downstream of the stop codon, due to the uncoupling of transcription and translation. In contrast, this translation alteration uniformly destabilizes the coding and ribosome-free regions, in a process triggered by RNase E activity, and its ability to form the RNA degradome. These results demonstrate how ribosomes protect mRNA molecules and highlight how translation, mRNA degradation and transcription are deeply interconnected in the quality control process that avoids unproductive gene expression in cells.

Lactococcus lactis CNCM I-5388 versus NCDO2118 by its GABA hyperproduction ability, counteracts faster stress-induced intestinal hypersensitivity in rats.

Irritable bowel syndrome (IBS) is a functional gastrointestinal disorder characterized by its main symptom, visceral hypersensitivity (VH), which is aggravated by stress. Gut-brain interactions and gut bacteria may alleviate IBS symptoms, including VH. γ-amino butyric acid (GABA), produced notably by lactic acid bacteria (LAB), shows promising result in IBS symptoms treatment. In bacteria, GABA is generated through glutamate decarboxylase (GAD) metabolism of L-glutamic acid, maintaining intracellular pH. In mammals, GABA acts as an inhibitory neurotransmitter, modulating pain, stress, and anxiety. Therefore, utilizing GABA-producing LAB as a therapeutic approach might be beneficial. Our previous work showed that a GABA-producing Lactococcus lactis strain, NCDO2118, reduced VH induced by acute stress in rats after a 10-day oral treatment. Here, we identified the strain CNCM I-5388, with a four-fold higher GABA production rate under the same conditions as NCDO2118. Both strains shared 99.1% identical GAD amino acid sequences and in vitro analyses revealed the same optimal pH for GAD activity; however, CNCM I-5388 exhibited 17 times higher intracellular GAD activity and increased resistance to acidic pH. Additionally, in vivo experiments have demonstrated that CNCM I-5388 has faster anti-VH properties in rats compared with NCDO2118, starting from the fifth day of treatment. Finally, CNCM I-5388 anti-VH effects partially persisted after 5-day treatment interruption and after a single oral treatment. These findings highlight CNCM I-5388 as a potential therapeutic agent for managing VH in IBS patients.

Natural diversity of lactococci in γ-aminobutyric acid (GABA) production and genetic and phenotypic determinants.

BACKGROUND: γ-aminobutyric acid (GABA) is a bioactive compound produced by lactic acid bacteria (LAB). The diversity of GABA production in the Lactococcus genus is poorly understood. Genotypic and phenotypic approaches were therefore combined in this study to shed light on this diversity. A comparative genomic study was performed on the GAD-system genes (gadR, gadC and gadB) involved in GABA production in 36 lactococci including L. lactis and L. cremoris species. In addition, 132 Lactococcus strains were screened for GABA production in culture medium supplemented with 34 mM L-glutamic acid with or without NaCl (0.3 M). RESULTS: Comparative analysis of the nucleotide sequence alignments revealed the same genetic organization of the GAD system in all strains except one, which has an insertion sequence element (IS981) into the PgadCB promoter. This analysis also highlighted several deletions including a 3-bp deletion specific to the cremoris species located in the PgadR promoter, and a second 39-bp deletion specific to L. cremoris strains with a cremoris phenotype. Phenotypic analysis revealed that GABA production varied widely, but it was higher in L. lactis species than in L. cremoris, with an exceptional GABA production of up to 14 and 24 mM in two L. lactis strains. Moreover, adding chloride increased GABA production in some L. cremoris and L. lactis strains by a factor of up to 16 and GAD activity correlated well with GABA production. CONCLUSIONS: This genomic analysis unambiguously characterized the cremoris phenotype of L. cremoris species and modified GadB and GadR proteins explain why the corresponding strains do not produce GABA. Finally, we found that glutamate decarboxylase activity revealing GadB protein amount, varied widely between the strains and correlated well with GABA production both with and without chloride. As this protein level is associated to gene expression, the regulation of GAD gene expression was identified as a major contributor to this diversity.

Detachment of the RNA degradosome from the inner membrane of Escherichia coli results in a global slowdown of mRNA degradation, proteolysis of RNase E and increased turnover of ribosome-free transcripts.

The reason for RNase E attachment to the inner membrane is largely unknown. To understand the cell biology of RNA degradation, we have characterized a strain expressing RNase E lacking the membrane attachment site (cytoplasmic RNase E). Genome-wide data show a global slowdown in mRNA degradation. There is no correlation between mRNA stabilization and the function or cellular location of encoded proteins. The activity of cRNase E is comparable to the wild-type enzyme in vitro, but the mutant protein is unstable in vivo. Autoregulation of cRNase E synthesis compensates for protein instability. cRNase E associates with other proteins to assemble a cytoplasmic RNA degradosome. CsrB/C sRNAs, whose stability is regulated by membrane-associated CsrD, are stabilized. Membrane attachment of RNase E is thus necessary for CsrB/C turnover. In contrast to mRNA stability, ribosome-free transcripts are sensitive to inactivation by cRNase E. Our results show that effects on RNA degradation are not due to the differences in the activity or level of cRNase E, or failure to assemble the RNA degradosome. We propose that membrane attachment is necessary for RNase E stability, functional interactions with membrane-associated regulatory factors and protection of ribosome-free transcripts from premature interactions with RNase E in the nucleoid.

Competitive effects in bacterial mRNA decay.

In living organisms, the same enzyme catalyses the degradation of thousands of different mRNAs, but the possible influence of competing substrates has been largely ignored so far. We develop a simple mechanistic model of the coupled degradation of all cell mRNAs using the total quasi-steady-state approximation of the Michaelis-Menten framework. Numerical simulations of the model using carefully chosen parameters and analyses of rate sensitivity coefficients show how substrate competition alters mRNA decay. The model predictions reproduce and explain a number of experimental observations on mRNA decay following transcription arrest, such as delays before the onset of degradation, the occurrence of variable degradation profiles with increased non linearities and the negative correlation between mRNA half-life and concentration. The competition acts at different levels, through the initial concentration of cell mRNAs and by modifying the enzyme affinity for its targets. The consequence is a global slow down of mRNA decay due to enzyme titration and the amplification of its apparent affinity. Competition happens to stabilize weakly affine mRNAs and to destabilize the most affine ones. We believe that this mechanistic model is an interesting alternative to the exponential models commonly used for the determination of mRNA half-lives. It allows analysing regulatory mechanisms of mRNA degradation and its predictions are directly comparable to experimental data.

The stability of an mRNA is influenced by its concentration: a potential physical mechanism to regulate gene expression.

Changing mRNA stability is a major post-transcriptional way of controlling gene expression, particularly in newly encountered conditions. As the concentration of mRNA is the result of an equilibrium between transcription and degradation, it is generally assumed that at constant transcription, any change in mRNA concentration is the consequence of mRNA stabilization or destabilization. However, the literature reports many cases of opposite variations in mRNA concentration and stability in bacteria. Here, we analyzed the causal link between the concentration and stability of mRNA in two phylogenetically distant bacteria Escherichia coli and Lactococcus lactis. Using reporter mRNAs, we showed that modifying the stability of an mRNA had unpredictable effects, either higher or lower, on its concentration, whereas increasing its concentration systematically reduced stability. This inverse relationship between the concentration and stability of mRNA was generalized to native genes at the genome scale in both bacteria. Higher mRNA turnover in the case of higher concentrations appears to be a simple physical mechanism to regulate gene expression in the bacterial kingdom. The consequences for bacterial adaptation of this control of the stability of an mRNA by its concentration are discussed.

PNPase is involved in the coordination of mRNA degradation and expression in stationary phase cells of Escherichia coli.

BACKGROUND: Exoribonucleases are crucial for RNA degradation in Escherichia coli but the roles of RNase R and PNPase and their potential overlap in stationary phase are not well characterized. Here, we used a genome-wide approach to determine how RNase R and PNPase affect the mRNA half-lives in the stationary phase. The genome-wide mRNA half-lives were determined by a dynamic analysis of transcriptomes after transcription arrest. We have combined the analysis of mRNA half-lives with the steady-state concentrations (transcriptome) to provide an integrated overview of the in vivo activity of these exoribonucleases at the genome-scale. RESULTS: The values of mRNA half-lives demonstrated that the mRNAs are very stable in the stationary phase and that the deletion of RNase R or PNPase caused only a limited mRNA stabilization. Intriguingly the absence of PNPase provoked also the destabilization of many mRNAs. These changes in mRNA half-lives in the PNPase deletion strain were associated with a massive reorganization of mRNA levels and also variation in several ncRNA concentrations. Finally, the in vivo activity of the degradation machinery was found frequently saturated by mRNAs in the PNPase mutant unlike in the RNase R mutant, suggesting that the degradation activity is limited by the deletion of PNPase but not by the deletion of RNase R. CONCLUSIONS: This work had identified PNPase as a central player associated with mRNA degradation in stationary phase.

Estimation of time-varying growth, uptake and excretion rates from dynamic metabolomics data.

MOTIVATION: Technological advances in metabolomics have made it possible to monitor the concentration of extracellular metabolites over time. From these data, it is possible to compute the rates of uptake and excretion of the metabolites by a growing cell population, providing precious information on the functioning of intracellular metabolism. The computation of the rate of these exchange reactions, however, is difficult to achieve in practice for a number of reasons, notably noisy measurements, correlations between the concentration profiles of the different extracellular metabolites, and discontinuties in the profiles due to sudden changes in metabolic regime. RESULTS: We present a method for precisely estimating time-varying uptake and excretion rates from time-series measurements of extracellular metabolite concentrations, specifically addressing all of the above issues. The estimation problem is formulated in a regularized Bayesian framework and solved by a combination of extended Kalman filtering and smoothing. The method is shown to improve upon methods based on spline smoothing of the data. Moreover, when applied to two actual datasets, the method recovers known features of overflow metabolism in Escherichia coli and Lactococcus lactis , and provides evidence for acetate uptake by L. lactis after glucose exhaustion. The results raise interesting perspectives for further work on rate estimation from measurements of intracellular metabolites. AVAILABILITY AND IMPLEMENTATION: The Matlab code for the estimation method is available for download at https://team.inria.fr/ibis/rate-estimation-software/ , together with the datasets. CONTACT: eugenio.cinquemani@inria.fr. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The post-transcriptional regulatory system CSR controls the balance of metabolic pools in upper glycolysis of Escherichia coli.

Metabolic control in Escherichia coli is a complex process involving multilevel regulatory systems but the involvement of post-transcriptional regulation is uncertain. The post-transcriptional factor CsrA is stated as being the only regulator essential for the use of glycolytic substrates. A dozen enzymes in the central carbon metabolism (CCM) have been reported as potentially controlled by CsrA, but its impact on the CCM functioning has not been demonstrated. Here, a multiscale analysis was performed in a wild-type strain and its isogenic mutant attenuated for CsrA (including growth parameters, gene expression levels, metabolite pools, abundance of enzymes and fluxes). Data integration and regulation analysis showed a coordinated control of the expression of glycolytic enzymes. This also revealed the imbalance of metabolite pools in the csrA mutant upper glycolysis, before the phosphofructokinase PfkA step. This imbalance is associated with a glucose-phosphate stress. Restoring PfkA activity in the csrA mutant strain suppressed this stress and increased the mutant growth rate on glucose. Thus, the carbon storage regulator system is essential for the effective functioning of the upper glycolysis mainly through its control of PfkA. This work demonstrates the pivotal role of post-transcriptional regulation to shape the carbon metabolism.

Adaptation of Propionibacterium freudenreichii to long-term survival under gradual nutritional shortage.

BACKGROUND: Propionibacterium freudenreichii is an Actinobacterium widely used in the dairy industry as a ripening culture for Swiss-type cheeses, for vitamin B12 production and some strains display probiotic properties. It is reportedly a hardy bacterium, able to survive the cheese-making process and digestive stresses. RESULTS: During this study, P. freudenreichii CIRM-BIA 138 (alias ITG P9), which has a generation time of five hours in Yeast Extract Lactate medium at 30 °C under microaerophilic conditions, was incubated for 11 days (9 days after entry into stationary phase) in a culture medium, without any adjunct during the incubation. The carbon and free amino acids sources available in the medium, and the organic acids produced by the strain, were monitored throughout growth and survival. Although lactate (the preferred carbon source for P. freudenreichii) was exhausted three days after inoculation, the strain sustained a high population level of 9.3 log10 CFU/mL. Its physiological adaptation was investigated by RNA-seq analysis and revealed a complete disruption of metabolism at the entry into stationary phase as compared to exponential phase. CONCLUSIONS: P. freudenreichii adapts its metabolism during entry into stationary phase by down-regulating oxidative phosphorylation, glycolysis, and the Wood-Werkman cycle by exploiting new nitrogen (glutamate, glycine, alanine) sources, by down-regulating the transcription, translation and secretion of protein. Utilization of polyphosphates was suggested.

Oligomerized backbone pilin helps piliated Lactococcus lactis to withstand shear flow.

The present work focuses on the role of pili present at the cell surface of Lactococcus lactis in bacterial adhesion to abiotic (hydrophobic polystyrene) and biotic (mucin-coated polystyrene) surfaces. Native pili-displaying strains and isogenic derivatives in which pilins or sortase C structural genes had been modified were used. Surface physico-chemistry, morphology and shear-flow-induced detachment of lactococcal cells were evaluated. The involvement of pili in L. lactis adhesion was clearly demonstrated, irrespective of the surface characteristics (hydrophobic/hydrophilic, presence or not of specific binding sites). The accessory pilin, PilC, and the backbone pilin, PilB, were revealed to play a major role in adhesion, provided that the PilB was present in its polymerized form. Within the population fraction that remained attached to the surface under increasing shear flow, different association behaviors were observed, showing that pili could serve as anchoring sites thus hampering the effect of shear flow on cell orientation and detachment.

Dual role of transcription and transcript stability in the regulation of gene expression in Escherichia coli cells cultured on glucose at different growth rates.

Microorganisms extensively reorganize gene expression to adjust growth rate to changes in growth conditions. At the genomic scale, we measured the contribution of both transcription and transcript stability to regulating messenger RNA (mRNA) concentration in Escherichia coli. Transcriptional control was the dominant regulatory process. Between growth rates of 0.10 and 0.63 h(-1), there was a generic increase in the bulk mRNA transcription. However, many transcripts became less stable and the median mRNA half-life decreased from 4.2 to 2.8 min. This is the first evidence that mRNA turnover is slower at extremely low-growth rates. The destabilization of many, but not all, transcripts at high-growth rate correlated with transcriptional upregulation of genes encoding the mRNA degradation machinery. We identified five classes of growth-rate regulation ranging from mainly transcriptional to mainly degradational. In general, differential stability within polycistronic messages encoded by operons does not appear to be affected by growth rate. We show here that the substantial reorganization of gene expression involving downregulation of tricarboxylic acid cycle genes and acetyl-CoA synthetase at high-growth rates is controlled mainly by transcript stability. Overall, our results demonstrate that the control of transcript stability has an important role in fine-tuning mRNA concentration during changes in growth rate.

Acetate Exposure Determines the Diauxic Behavior of Escherichia coli during the Glucose-Acetate Transition.

UNLABELLED: Growth of Escherichia coli on glucose in batch culture is accompanied by the excretion of acetate, which is consumed by the cells when glucose is exhausted. This glucose-acetate transition is classically described as a diauxie (two successive growth stages). Here, we investigated the physiological and metabolic properties of cells after glucose exhaustion through the analysis of growth parameters and gene expression. We found that E. coli cells grown on glucose in batch culture produce acetate and consume it after glucose exhaustion but do not grow on acetate. Acetate is catabolized, but key anabolic genes--such as the genes encoding enzymes of the glyoxylate shunt--are not upregulated, hence preventing growth. Both the induction of the latter anabolic genes and growth were observed only after prolonged exposure to low concentrations of acetate and could be accelerated by high acetate concentrations. We postulate that such decoupling between acetate catabolism and acetate anabolism might be an advantage for the survival of E. coli in the ever-changing environment of the intestine. IMPORTANCE: The glucose-acetate transition is a valuable experimental model for comprehensive investigations of metabolic adaptation and a current paradigm for developing modeling approaches in systems microbiology. Yet, the work reported in our paper demonstrates that the metabolic behavior of Escherichia coli during the glucose-acetate transition is much more complex than what has been reported so far. A decoupling between acetate catabolism and acetate anabolism was observed after glucose exhaustion, which has not been reported previously. This phenomenon could represent a strategy for optimal utilization of carbon resources during colonization and persistence of E. coli in the gut and is also of significant interest for biotechnological applications.

Genome-wide investigation of mRNA lifetime determinants in Escherichia coli cells cultured at different growth rates.

BACKGROUND: Changes to mRNA lifetime adjust mRNA concentration, facilitating the adaptation of growth rate to changes in growth conditions. However, the mechanisms regulating mRNA lifetime are poorly understood at the genome-wide scale and have not been investigated in bacteria growing at different rates. RESULTS: We used linear covariance models and the best model selected according to the Akaike information criterion to identify and rank intrinsic and extrinsic general transcript parameters correlated with mRNA lifetime, using mRNA half-life datasets for E. coli, obtained at four growth rates. The principal parameter correlated with mRNA stability was mRNA concentration, the mRNAs most concentrated in the cells being the least stable. However, sequence-related features (codon adaptation index (CAI), ORF length, GC content, polycistronic mRNA), gene function and essentiality also affected mRNA lifetime at all growth rates. We also identified sequence motifs within the 5'UTRs potentially related to mRNA stability. Growth rate-dependent effects were confined to particular functional categories (e.g. carbohydrate and nucleotide metabolism). Finally, mRNA stability was less strongly correlated with the amount of protein produced than mRNA concentration and CAI. CONCLUSIONS: This study provides the most complete genome-wide analysis to date of the general factors correlated with mRNA lifetime in E. coli. We have generalized for the entire population of transcripts or excluded determinants previously defined as regulators of stability for some particular mRNAs and identified new, unexpected general indicators. These results will pave the way for discussions of the underlying mechanisms and their interaction with the growth physiology of bacteria.

A genome-scale integration and analysis of Lactococcus lactis translation data.

Protein synthesis is a template polymerization process composed by three main steps: initiation, elongation, and termination. During translation, ribosomes are engaged into polysomes whose size is used for the quantitative characterization of translatome. However, simultaneous transcription and translation in the bacterial cytosol complicates the analysis of translatome data. We established a procedure for robust estimation of the ribosomal density in hundreds of genes from Lactococcus lactis polysome size measurements. We used a mechanistic model of translation to integrate the information about the ribosomal density and for the first time we estimated the protein synthesis rate for each gene and identified the rate limiting steps. Contrary to conventional considerations, we find significant number of genes to be elongation limited. This number increases during stress conditions compared to optimal growth and proteins synthesized at maximum rate are predominantly elongation limited. Consistent with bacterial physiology, we found proteins with similar rate and control characteristics belonging to the same functional categories. Under stress conditions, we found that synthesis rate of regulatory proteins is becoming comparable to proteins favored under optimal growth. These findings suggest that the coupling of metabolic states and protein synthesis is more important than previously thought.

The significance of translation regulation in the stress response.

BACKGROUND: The stress response in bacteria involves the multistage control of gene expression but is not entirely understood. To identify the translational response of bacteria in stress conditions and assess its contribution to the regulation of gene expression, the translational states of all mRNAs were compared under optimal growth condition and during nutrient (isoleucine) starvation. RESULTS: A genome-scale study of the translational response to nutritional limitation was performed in the model bacterium Lactococcus lactis. Two measures were used to assess the translational status of each individual mRNA: the fraction engaged in translation (ribosome occupancy) and ribosome density (number of ribosomes per 100 nucleotides). Under isoleucine starvation, half of the mRNAs considered were translationally down-regulated mainly due to decreased ribosome density. This pattern concerned genes involved in growth-related functions such as translation, transcription, and the metabolism of fatty acids, phospholipids and bases, contributing to the slowdown of growth. Only 4% of the mRNAs were translationally up-regulated, mostly related to prophagic expression in response to stress. The remaining genes exhibited antagonistic regulations of the two markers of translation. Ribosome occupancy increased significantly for all the genes involved in the biosynthesis of isoleucine, although their ribosome density had decreased. The results revealed complex translational regulation of this pathway, essential to cope with isoleucine starvation.To elucidate the regulation of global gene expression more generally, translational regulation was compared to transcriptional regulation under isoleucine starvation and to other post-transcriptional regulations related to mRNA degradation and mRNA dilution by growth. Translational regulation appeared to accentuate the effects of transcriptional changes for down-regulated growth-related functions under isoleucine starvation although mRNA stabilization and lower dilution by growth counterbalanced this effect. CONCLUSIONS: We show that the contribution of translational regulation to the control of gene expression is significant in the stress response. Post-transcriptional regulation is complex and not systematically co-directional with transcription regulation. Post-transcriptional regulation is important to the understanding of gene expression control.

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.

The carbohydrate metabolism signature of lactococcus lactis strain A12 reveals its sourdough ecosystem origin.

Lactococcus lactis subsp. lactis strain A12 was isolated from sourdough. Combined genomic, transcriptomic, and phenotypic analyses were performed to understand its survival capacity in the complex sourdough ecosystem and its role in the microbial community. The genome sequence comparison of strain A12 with strain IL1403 (a derivative of an industrial dairy strain) revealed 78 strain-specific regions representing 23% of the total genome size. Most of the strain-specific genes were involved in carbohydrate metabolism and are potentially required for its persistence in sourdough. Phenotype microarray, growth tests, and analysis of glycoside hydrolase content showed that strain A12 fermented plant-derived carbohydrates, such as arabinose and α-galactosides. Strain A12 exhibited specific growth rates on raffinose that were as high as they were on glucose and was able to release sucrose and galactose outside the cell, providing soluble carbohydrates for sourdough microflora. Transcriptomic analysis identified genes specifically induced during growth on raffinose and arabinose and reveals an alternative pathway for raffinose assimilation to that used by other lactococci.