Transcriptome landscape of Lactococcus lactis reveals many novel RNAs including a small regulatory RNA involved in carbon uptake and metabolism

Ribosome profiling reveals downregulation of UMP biosynthesis as the major early response to phage infection

Bacteria have evolved diverse defense mechanisms to counter bacteriophage attacks. Genetic programs activated upon infection characterize phage-host molecular interactions and ultimately determine the outcome of the infection. In this study, we applied ribosome profiling to monitor protein synthesis during the early stages of sk1 bacteriophage infection in Lactococcus cremoris. Our analysis revealed major changes in gene expression within 5 minutes of sk1 infection. Notably, we observed a specific and severe downregulation of several pyr operons which encode enzymes required for uridine monophosphate biosynthesis. Consistent with previous findings, this is likely an attempt of the host to starve the phage of nucleotides it requires for propagation. We also observed a gene expression response that we expect to benefit the phage. This included the upregulation of 40 ribosome proteins that likely increased the host's translational capacity, concurrent with a downregulation of genes that promote translational fidelity (lepA and raiA). In addition to the characterization of host-phage gene expression responses, the obtained ribosome profiling data enabled us to identify two putative recoding events as well as dozens of loci currently annotated as pseudogenes that are actively translated. Furthermore, our study elucidated alterations in the dynamics of the translation process, as indicated by time-dependent changes in the metagene profile, suggesting global shifts in translation rates upon infection. Additionally, we observed consistent modifications in the ribosome profiles of individual genes, which were apparent as early as 2 minutes post-infection. The study emphasizes our ability to capture rapid alterations of gene expression during phage infection through ribosome profiling.

IMPORTANCE

The ribosome profiling technology has provided invaluable insights for understanding cellular translation and eukaryotic viral infections. However, its potential for investigating host-phage interactions remains largely untapped. Here, we applied ribosome profiling to Lactococcus cremoris cultures infected with sk1, a major infectious agent in dairy fermentation processes. This revealed a profound downregulation of genes involved in pyrimidine nucleotide synthesis at an early stage of phage infection, suggesting an anti-phage program aimed at restricting nucleotide availability and, consequently, phage propagation. This is consistent with recent findings and contributes to our growing appreciation for the role of nucleotide limitation as an anti-viral strategy. In addition to capturing rapid alterations in gene expression levels, we identified translation occurring outside annotated regions, as well as signatures of non-standard translation mechanisms. The gene profiles revealed specific changes in ribosomal densities upon infection, reflecting alterations in the dynamics of the translation process.

The growth-survival trade-off is hard-wired in the Lactococcus lactis gene regulation network

Most microbes reside in oligotrophic environments for extended periods of time, requiring survival strategies that maintain proliferative capacity. We demonstrate that the non-spore-forming Lactococcus lactis KF147 progressively activates the expression of stress-associated functions in response to the declining growth rate elicited by prolonged retentostat cultivation, which coincides with up to 10⁴ -fold increased stress tolerance. Our findings provide a quantified view of the transcription and stress-tolerance adaptations underlying the growth-survival trade-off in L. lactis, and exemplify the hard-wiring of this trade-off in the lactococcal gene regulation network.

Omics Approaches to Assess Flavor Development in Cheese

Cheese is characterized by a rich and complex microbiota that plays a vital role during both production and ripening, contributing significantly to the safety, quality, and sensory characteristics of the final product. In this context, it is vital to explore the microbiota composition and understand its dynamics and evolution during cheese manufacturing and ripening. Application of high-throughput DNA sequencing technologies have facilitated the more accurate identification of the cheese microbiome, detailed study of its potential functionality, and its contribution to the development of specific organoleptic properties. These technologies include amplicon sequencing, whole-metagenome shotgun sequencing, metatranscriptomics, and, most recently, metabolomics. In recent years, however, the application of multiple meta-omics approaches along with data integration analysis, which was enabled by advanced computational and bioinformatics tools, paved the way to better comprehension of the cheese ripening process, revealing significant associations between the cheese microbiota and metabolites, as well as their impact on cheese flavor and quality.

Riboswitch RSthiT as a molecular tool in Lactococcus lactis

Previous RNA sequencing has allowed the identification of 129 long 5′ untranslated regions (UTRs) in the Lactococcus lactis MG1363 transcriptome. These sequences potentially harbor cis-acting riboswitches. One of the identified extended 5′ UTRs is a putative thiamine pyrophosphate (TPP) riboswitch. It is located immediately upstream of the thiamine transporter gene thiT (llmg_0334). To confirm this assumption, the 5′-UTR sequence was placed upstream of the gene encoding the superfolder green fluorescent protein (sfGFP), sfgfp, allowing the examination of the expression of sfGFP in the presence or absence of thiamine in the medium. The results show that this sequence indeed represents a thiamine-responsive TPP riboswitch. This RNA-based genetic control device was used to successfully restore the mutant phenotype of an L. lactis strain lacking the major autolysin gene, acmA. The L. lactisthiT TPP riboswitch (RS_thiT) is a useful molecular genetic tool enabling the gradual downregulation of the expression of genes under its control by adjusting the thiamine concentration.

IMPORTANCE The capacity of microbes with biotechnological importance to adapt to and survive under quickly changing industrial conditions depends on their ability to adequately control gene expression. Riboswitches are important RNA-based elements involved in rapid and precise gene regulation. Here, we present the identification of a natural thiamine-responsive riboswitch of Lactococcus lactis, a bacterium used worldwide in the production of dairy products. We used it to restore a genetic defect in an L. lactis mutant and show that it is a valuable addition to the ever-expanding L. lactis genetic toolbox.

Profiles of Small Regulatory RNAs at Different Growth Phases of Streptococcus thermophilus During pH-Controlled Batch Fermentation

Small regulatory RNA (sRNA) has been shown to play an important role under various stress conditions in bacteria, and it plays a vital role in regulating growth, adaptation and survival through posttranscriptional control of gene expression in bacterial cells. Streptococcus thermophilus is widely used as a starter culture in the manufacture of fermented dairy products. However, the lack of reliable information on the expression profiles and potential physiological functions of sRNAs in this species hinders our understanding of the importance of sRNAs in S. thermophilus. The present study was conducted to assess the expression profiles of sRNAs in S. thermophilus and to identify sRNAs that exhibited significant changes. A total of 530 potential sRNAs were identified, including 198 asRNAs, 135 sRNAs from intergenic regions, and 197 sRNAs from untranslated regions (UTRs). Significant changes occurred in the expression of 238, 83, 194, and 139 sRNA genes during the lag, early exponential growth, late exponential growth, and stationary phases, respectively. The expression of 14 of the identified sRNAs was verified by qRT-PCR. Predictions of the target genes of these candidate sRNAs showed that the primary metabolic pathways targeted were involved in carbon metabolism, biosynthesis of amino acids, ABC transporters, the metabolism of amino and nucleotide sugars, purine metabolism, and the phosphotransferase system. The expression of the predicted target genes was further analyzed to better understand the roles of sRNAs during different growth stages. The results suggested that these sRNAs play crucial roles by regulating biological pathways during different growth phases of S. thermophilus. According to the results, sRNAs sts141, sts392, sts318, and sts014 are involved in the regulation of osmotic stress. sRNAs sts508, sts087, sts372, sts141, sts375, and sts119 are involved in the regulation of starvation stress. sRNAs sts129, sts226, sts166, sts231, sts204, sts145, and sts236 are involved in arginine synthesis. sRNAs sts033, sts341, sts492, sts140, sts230, sts172, and sts377 are involved in the ADI pathway. The present study provided valuable information for the functional study of sRNAs in S. thermophilus and indicated a future research direction for sRNA in S. thermophilus. Overall, our results provided new insights for understanding the complex regulatory network of sRNAs in S. thermophilus.

Insights into 6S RNA in lactic acid bacteria (LAB)

BACKGROUND

6S RNA is a regulator of cellular transcription that tunes the metabolism of cells. This small non-coding RNA is found in nearly all bacteria and among the most abundant transcripts. Lactic acid bacteria (LAB) constitute a group of microorganisms with strong biotechnological relevance, often exploited as starter cultures for industrial products through fermentation. Some strains are used as probiotics while others represent potential pathogens. Occasional reports of 6S RNA within this group already indicate striking metabolic implications. A conceivable idea is that LAB with 6S RNA defects may metabolize nutrients faster, as inferred from studies of Echerichia coli. This may accelerate fermentation processes with the potential to reduce production costs. Similarly, elevated levels of secondary metabolites might be produced. Evidence for this possibility comes from preliminary findings regarding the production of surfactin in Bacillus subtilis, which has functions similar to those of bacteriocins. The prerequisite for its potential biotechnological utility is a general characterization of 6S RNA in LAB.

RESULTS

We provide a genomic annotation of 6S RNA throughout the Lactobacillales order. It laid the foundation for a bioinformatic characterization of common 6S RNA features. This covers secondary structures, synteny, phylogeny, and product RNA start sites. The canonical 6S RNA structure is formed by a central bulge flanked by helical arms and a template site for product RNA synthesis. 6S RNA exhibits strong syntenic conservation. It is usually flanked by the replication-associated recombination protein A and the universal stress protein A. A catabolite responsive element was identified in over a third of all 6S RNA genes. It is known to modulate gene expression based on the available carbon sources. The presence of antisense transcripts could not be verified as a general trait of LAB 6S RNAs.

CONCLUSIONS

Despite a large number of species and the heterogeneity of LAB, the stress regulator 6S RNA is well-conserved both from a structural as well as a syntenic perspective. This is the first approach to describe 6S RNAs and short 6S RNA-derived transcripts beyond a single species, spanning a large taxonomic group covering multiple families. It yields universal insights into this regulator and complements the findings derived from other bacterial model organisms.

Involvement of non-coding RNAs during infection of rice by Acidovorax oryzae

The expression of non-coding RNAs (ncRNAs) has been observed in a variety of bacteria. However, the function of ncRNAs and their regulatory targets are largely unknown, and few ncRNAs are found to be associated with bacterial virulence. The bacterial brown stripe pathogen Acidovorax oryzae (Ao) RS-1 shows a high level of condition-dependent differential expression of ncRNA, which we identified in a genome wide screen. We experimentally validated 66 differentially expressed ncRNAs using an integrative analysis of conservative genome sequences and transcriptomic data during in vivo interaction of the bacterial pathogen with the rice plant. To test the relevance of the differentially expressed ncRNAs, we chose four with different positions within the genome, and with different secondary structures and promoter activities. The results show that the overexpression of the four ncRNAs caused a significant change in virulence-related phenotypes, resistance to various environmental stresses, expression of secretion systems and effector proteins, while changing the expression of ncRNA putative target genes. We conclude that these ncRNAs are examples for the inherent regulatory roles for many of the observed ncRNAs in response to changing conditions such as host interaction or environmental adaption.

Deciphering the Regulation of the Mannitol Operon Paves the Way for Efficient Production of Mannitol in Lactococcus lactis

Lactococcus lactis has great potential for high-yield production of mannitol, which has not yet been fully realized. In this study, we characterize how the mannitol genes in L. lactis are organized and regulated and use this information to establish efficient mannitol production. Although the organization of the mannitol genes in L. lactis was similar to that in other Gram-positive bacteria, mtlF and mtlD, encoding the enzyme IIA component (EIIA^mtl) of the mannitol phosphotransferase system (PTS) and the mannitol-1-phosphate dehydrogenase, respectively, were separated by a transcriptional terminator, and the mannitol genes were found to be organized in two transcriptional units: an operon comprising mtlA, encoding the enzyme IIBC component (EIIBC^mtl) of the mannitol PTS, mtlR, encoding a transcriptional activator, and mtlF, as well as a separately expressed mtlD gene. The promoters driving expression of the two transcriptional units were somewhat similar, and both contained predicted catabolite responsive element (cre) genes. The presence of carbon catabolite repression was demonstrated and was shown to be relieved in stationary-phase cells. The transcriptional activator MtlR (mtlR), in some Gram-positive bacteria, is repressed by phosphorylation by EIIA^mtl, and when we knocked out mtlF, we indeed observed enhanced expression from the two promoters, which indicated that this mechanism was in place. Finally, by overexpressing the mtlD gene and using stationary-phase cells as biocatalysts, we attained 10.1 g/liter mannitol with a 55% yield, which, to the best of our knowledge, is the highest titer ever reported for L. lactis. Summing up, the results of our study should be useful for improving the mannitol-producing capacity of this important industrial organism. IMPORTANCE Lactococcus lactis is the most studied species of the lactic acid bacteria, and it is widely used in various food fermentations. To date, there have been several attempts to persuade L. lactis to produce mannitol, a sugar alcohol with important therapeutic and food applications. Until now, to achieve mannitol production in L. lactis with significant titer and yield, it has been necessary to introduce and express foreign genes, which precludes the use of such strains in foods, due to their recombinant status. In this study, we systematically characterize how the mannitol genes in L. lactis are regulated and demonstrate how this impacts mannitol production capability. We harnessed this information and managed to establish efficient mannitol production without introducing foreign genes.

A Small Regulatory RNA Generated from the malK 5' Untranslated Region Targets Gluconeogenesis in Vibrio Species

Vsr217 is a small RNA from Vibrio tasmaniensis LGP32, a pathogen associated with mortality events affecting juvenile oysters. The vsr217 gene is located within the 5' untranslated region (UTR) of malK, encoding the ATPase component of the maltose importer, and is conserved within the genus Vibrio. In the presence of maltose, vsr217 is regulated by MalT, the positive regulator of the maltose regulon. vsr217 is required in cis for the full expression of malK. In addition, Vsr217 acts in trans to downregulate the expression of fbp encoding fructose-1,6-bisphosphatase, an enzyme involved in gluconeogenesis. Thus, in the presence of maltose, the induction of Vsr217 is expected to promote glycolysis by negatively regulating the expression of a key enzyme of gluconeogenesis. IMPORTANCE Juvenile pacific oysters have been subject in recent years to summer mortality episodes with deep economic consequences. The pathogen Vibrio tasmaniensis has been associated with such mortality events. For bacterial pathogens, survival within the host requires profound metabolic adaptations according to available resources. All kinds of regulatory elements, including noncoding RNAs, orchestrate this response. Oysters are rich in glycogen, a precursor of maltose, and we previously reported that V. tasmaniensis maltose-regulated genes are strongly induced during oyster infection. Here, we report the dual mechanism by which a small regulatory RNA, generated from the 5' untranslated region of a gene belonging to the maltose regulon, acts both in cis and trans. In cis, it stimulates growth on maltose, and in trans, it downregulates the expression of a gene associated with gluconeogenesis, thus coordinating maltose utilization with central carbon metabolism.

Structure of a Core Promoter in Bifidobacterium longum NCC2705

Bacterial promoters consist of core sequence motifs termed -35 and -10 boxes. The consensus motifs are TTGACA and TATAAT, respectively, which were identified from leading investigations on Escherichia coli However, the consensus sequences are not likely to fit genetically divergent bacteria. The sigma factor of the genus Bifidobacterium has a characteristic polar domain in the N terminus, suggesting the possibility of specific promoter recognition. We reevaluated the structure of Bifidobacterium longum NCC2705 promoters and compared them to other bacteria. Transcriptional start sites (TSSs) of the B. longum NCC2705 strain were identified using transcriptome sequencing (RNA-Seq) analysis to extract promoter regions. Conserved motifs of a bifidobacterial promoter were determined using regions upstream of TSSs and a hidden Markov model. As a result, consensus motifs of the -35 and -10 boxes were TTGTGC and TACAAT, respectively. To assess each base of both motifs, we constructed 37 plasmids based on pKO403-TPCTcon, including the hup promoter connected with a chloramphenicol acetyltransferase as a reporter gene. This reporter assay showed two optimal motifs of the -35 and -10 boxes, namely, TTGNNN and TANNNT, respectively. We further analyzed spacer lengths between the -35 and -10 boxes via a bioinformatics approach. The spacer lengths predominant in bacteria have been generally reported to be approximately 17 bp. In contrast, the predominant spacer lengths in the genus Bifidobacterium and related species were 11 bp, in addition to 17 bp. A reporter assay to assess the spacer lengths indicated that the 11-bp spacer length produced unusually high activity.IMPORTANCE The structures of sigma factors vary among bacterial strains, indicating that recognition rules may also vary. Therefore, we investigated the promoter structure of Bifidobacterium longum NCC2705 using a bioinformatics approach and wet analyses. The most frequent and optimal motifs were similar to other bacterial consensus motifs. The optimal spacer length between the two boxes was reported to be 17 bp. It is widely applied to a bioinformatics approach for other bacteria. Unexpectedly, conserved spacer lengths were 11 bp as well as 17 bp in the genus Bifidobacterium Moreover, the sigma factor of the genus Bifidobacterium has a characteristic domain in the N terminus which may contribute to the additional functions. Hence, it would be valuable to reevaluate the promoter in other organisms.

The Signal and the Noise: Characteristics of Antisense RNA in Complex Microbial Communities

This study systematically evaluated the global patterns of microbial antisense expression across various environments and provides a bird’s-eye view of general patterns observed across data sets, which can provide guidelines in our understanding of antisense expression as well as interpretation of metatranscriptomic data in general. This analysis highlights that in some environments, antisense expression from microbial communities can dominate over regular gene expression. We explored some potential drivers of antisense transcription, but more importantly, this study serves as a starting point, highlighting topics for future research and providing guidelines to include antisense expression in generic bioinformatic pipelines for metatranscriptomic data.

High-throughput sequencing has allowed unprecedented insight into the composition and function of complex microbial communities. With metatranscriptomics, it is possible to interrogate the transcriptomes of multiple organisms simultaneously to get an overview of the gene expression of the entire community. Studies have successfully used metatranscriptomics to identify and describe relationships between gene expression levels and community characteristics. However, metatranscriptomic data sets contain a rich suite of additional information that is just beginning to be explored. Here, we focus on antisense expression in metatranscriptomics, discuss the different computational strategies for handling it, and highlight the strengths but also potentially detrimental effects on downstream analysis and interpretation. We also analyzed the antisense transcriptomes of multiple genomes and metagenome-assembled genomes (MAGs) from five different data sets and found high variability in the levels of antisense transcription for individual species, which were consistent across samples. Importantly, we challenged the conceptual framework that antisense transcription is primarily the product of transcriptional noise and found mixed support, suggesting that the total observed antisense RNA in complex communities arises from the combined effect of unknown biological and technical factors. Antisense transcription can be highly informative, including technical details about data quality and novel insight into the biology of complex microbial communities.

IMPORTANCE This study systematically evaluated the global patterns of microbial antisense expression across various environments and provides a bird’s-eye view of general patterns observed across data sets, which can provide guidelines in our understanding of antisense expression as well as interpretation of metatranscriptomic data in general. This analysis highlights that in some environments, antisense expression from microbial communities can dominate over regular gene expression. We explored some potential drivers of antisense transcription, but more importantly, this study serves as a starting point, highlighting topics for future research and providing guidelines to include antisense expression in generic bioinformatic pipelines for metatranscriptomic data.

A dual-chain assembly pathway generates the high structural diversity of cell-wall polysaccharides in Lactococcus lactis

In Lactococcus lactis, cell-wall polysaccharides (CWPSs) act as receptors for many bacteriophages, and their structural diversity among strains explains, at least partially, the narrow host range of these viral predators. Previous studies have reported that lactococcal CWPS consists of two distinct components, a variable chain exposed at the bacterial surface, named polysaccharide pellicle (PSP), and a more conserved rhamnan chain anchored to, and embedded inside, peptidoglycan. These two chains appear to be covalently linked to form a large heteropolysaccharide. The molecular machinery for biosynthesis of both components is encoded by a large gene cluster, named cwps In this study, using a CRISPR/Cas-based method, we performed a mutational analysis of the cwps genes. MALDI-TOF MS-based structural analysis of the mutant CWPS combined with sequence homology, transmission EM, and phage sensitivity analyses enabled us to infer a role for each protein encoded by the cwps cluster. We propose a comprehensive CWPS biosynthesis scheme in which the rhamnan and PSP chains are independently synthesized from two distinct lipid-sugar precursors and are joined at the extracellular side of the cytoplasmic membrane by a mechanism involving a membrane-embedded glycosyltransferase with a GT-C fold. The proposed scheme encompasses a system that allows extracytoplasmic modification of rhamnan by complex substituting oligo-/polysaccharides. It accounts for the extensive diversity of CWPS structures observed among lactococci and may also have relevance to the biosynthesis of complex rhamnose-containing CWPSs in other Gram-positive bacteria.

Small RNAs in cell-to-cell communications during bacterial infection

Intercellular communication is a widespread phenomenon in all domains of life. Bacteria have developed many ways of communicating with one another and with other species, either prokaryotic or eukaryotic. RNA has been a key molecule since the beginning of life on Earth, and is one of the carriers of information. Given the current antibiotic crisis, understanding the way in which pathogens communicate can lead towards improved ways to control infections when antimicrobial therapy is not possible. Different subspecies of RNA, non-coding, and of small size, designated here as ncRNAs, have been in recent years the subject of a great research effort, and results have contributed to a growing field of knowledge. This review focuses on four different aspects of ncRNA involvement in cell-to-cell communications during bacterial infections: pathogen recognition by the host, alteration of host microRNA profiles, production of domestic and secreted forms of ncRNAs and subversion of the host responses. The current review article focuses on the most recent discoveries in the field and gives an integrative idea based on the discussed studies.

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.

Deep genome annotation of the opportunistic human pathogen Streptococcus pneumoniae D39

A precise understanding of the genomic organization into transcriptional units and their regulation is essential for our comprehension of opportunistic human pathogens and how they cause disease. Using single-molecule real-time (PacBio) sequencing we unambiguously determined the genome sequence of Streptococcus pneumoniae strain D39 and revealed several inversions previously undetected by short-read sequencing. Significantly, a chromosomal inversion results in antigenic variation of PhtD, an important surface-exposed virulence factor. We generated a new genome annotation using automated tools, followed by manual curation, reflecting the current knowledge in the field. By combining sequence-driven terminator prediction, deep paired-end transcriptome sequencing and enrichment of primary transcripts by Cappable-Seq, we mapped 1015 transcriptional start sites and 748 termination sites. We show that the pneumococcal transcriptional landscape is complex and includes many secondary, antisense and internal promoters. Using this new genomic map, we identified several new small RNAs (sRNAs), RNA switches (including sixteen previously misidentified as sRNAs), and antisense RNAs. In total, we annotated 89 new protein-encoding genes, 34 sRNAs and 165 pseudogenes, bringing the S. pneumoniae D39 repertoire to 2146 genetic elements. We report operon structures and observed that 9% of operons are leaderless. The genome data are accessible in an online resource called PneumoBrowse (https://veeninglab.com/pneumobrowse) providing one of the most complete inventories of a bacterial genome to date. PneumoBrowse will accelerate pneumococcal research and the development of new prevention and treatment strategies.

The protein regulator ArgR and the sRNA derived from the 3'-UTR region of its gene, ArgX, both regulate the arginine deiminase pathway in Lactococcus lactis

Small regulatory RNAs (sRNAs) and their enormous potential and versatility have provided us with an astounding insight in the complexity of bacterial transcriptomes. sRNAs have been shown to be involved in a variety of cellular processes that range from stress to general metabolism. Here we report that the gene encoding the transcriptional regulator ArgR is immediately followed by the gene of the small regulatory RNA ArgX. The latter is transcribed from its own promoter. The production of ArgX is induced by increasing arginine concentrations and repressed by CcpA. Previously, ArgR was shown to act as a transcriptional repressor of the catabolic arginine deiminase pathway (arc operon) by binding in the promoter region of arcA. Here we demonstrate that ArgX downregulates arc mRNA levels. Furthermore, ArgX putatively blocks the translation of one of the genes in the operon, arcC1, a process that would redirect an intermediate in arginine degradation, carbamoyl phosphate, towards pyrimidine synthesis. Our findings exemplify, for the first time, the combinatorial power of a transcription factor and a small regulatory RNA derived from the 3’-UTR region. The regulators ArgR and ArgX share a common target, but act on transcription and on RNA level, respectively.

Omics-based comparative analysis of putative mobile genetic elements in Lactococcus lactis

Lactococcus lactis is globally used in food fermentation. Genomics is useful to investigate speciation and differential occurrence of (un)desired gene functions, often related to mobile DNA. This study investigates L. lactis for putative chromosomal mobile genetic elements through comparative genomics, and analyses how they contribute to chromosomal variation at strain level. Our work identified 95 loci that may range over 10% of the chromosome size when including prophages, and the loci display a marked differential occurrence in the analysed strains. Analysis of differential transcriptomics data revealed how mobile genetic elements may impact the host physiology in response to conditional changes. This insight in the genetic variation of mobile genetic elements in L. lactis holds potential to further identify important functions related to food and biotechnology applications within this important species.

The Primary Antisense Transcriptome of Halobacterium salinarum NRC-1

Antisense RNAs (asRNAs) are present in diverse organisms and play important roles in gene regulation. In this work, we mapped the primary antisense transcriptome in the halophilic archaeon Halobacterium salinarum NRC-1. By reanalyzing publicly available data, we mapped antisense transcription start sites (aTSSs) and inferred the probable 3′ ends of these transcripts. We analyzed the resulting asRNAs according to the size, location, function of genes on the opposite strand, expression levels and conservation. We show that at least 21% of the genes contain asRNAs in H. salinarum. Most of these asRNAs are expressed at low levels. They are located antisense to genes related to distinctive characteristics of H. salinarum, such as bacteriorhodopsin, gas vesicles, transposases and other important biological processes such as translation. We provide evidence to support asRNAs in type II toxin–antitoxin systems in archaea. We also analyzed public Ribosome profiling (Ribo-seq) data and found that ~10% of the asRNAs are ribosome-associated non-coding RNAs (rancRNAs), with asRNAs from transposases overrepresented. Using a comparative transcriptomics approach, we found that ~19% of the asRNAs annotated in H. salinarum belong to genes with an ortholog in Haloferax volcanii, in which an aTSS could be identified with positional equivalence. This shows that most asRNAs are not conserved between these halophilic archaea.

A computational system for identifying operons based on RNA-seq data

An operon is a set of neighboring genes in a genome that is transcribed as a single polycistronic message. Genes that are part of the same operon often have related functional roles or participate in the same metabolic pathways. The majority of all bacterial genes are co-transcribed with one or more other genes as part of a multi-gene operon. Thus, accurate identification of operons is important in understanding co-regulation of genes and their functional relationships. Here, we present a computational system that uses RNA-seq data to determine operons throughout a genome. The system takes the name of a genome and one or more files of RNA-seq data as input. Our method combines primary genomic sequence information with expression data from the RNA-seq files in a unified probabilistic model in order to identify operons. We assess our method's ability to accurately identify operons in a range of species through comparison to external databases of operons, both experimentally confirmed and computationally predicted, and through focused experiments that confirm new operons identified by our method. Our system is freely available at https://cs.wellesley.edu/~btjaden/Rockhopper/.

The Many Facets of the Small Non-coding RNA RsaE (RoxS) in Metabolic Niche Adaptation of Gram-Positive Bacteria

Small regulatory RNAs (sRNAs) are increasingly recognized as players in the complex regulatory networks governing bacterial gene expression. RsaE (synonym RoxS) is an sRNA that is highly conserved in bacteria of the Bacillales order. Recent analyses in Bacillus subtilis, Staphylococcus aureus and Staphylococcus epidermidis identified RsaE/RoxS as a potent riboregulator of central carbon metabolism and energy balance with many molecular RsaE/RoxS functions and targets being shared across species. Similarities and species-specific differences in cellular processes modulated by RsaE/RoxS suggest that this sRNA plays a prominent role in the adaptation of Gram-positive bacteria to niches with varying nutrient availabilities and environmental cues. This review summarizes recent findings on the molecular function of RsaE/RoxS and its interaction with mRNA targets. Special emphasis will be on the integration of RsaE/RoxS into metabolic regulatory circuits and, derived from this, the role of RsaE/RoxS as a putative driver to generate phenotypic heterogeneity in bacterial populations. In this respect, we will particularly discuss heterogeneous RsaE expression in S. epidermidis biofilms and its possible contribution to metabolic niche diversification, programmed bacterial lysis and biofilm matrix production.

The genome and transcriptome of Lactococcus lactis ssp. lactis F44 and G423: Insights into adaptation to the acidic environment

Nisin, as a common green (environmentally friendly), nontoxic antibacterial peptide secreted by Lactococcus lactis, is widely used to prevent the decomposition of meat and dairy products and maintains relatively high stability at low pH. However, the growth of Lc. lactis is frequently inhibited by high lactic acid concentrations produced during fermentation. This phenomenon has become a great challenge in enhancing the nisin yield for this strain. Here, the shuffled strain G423 that could survive on a solid plate at pH 3.7 was generated through protoplast fusion-mediated genome shuffling. The nisin titer of G423 peaked at 4,543 IU/mL, which was 59.9% higher than that of the same batch of the initial strain Lc. lactis F44. The whole genome comparisons between G423 and F44 indicated that 6 large fragments (86,725 bp) were inserted in G423 compared with that of Lc. lactis F44. Transcriptome data revealed that 4 novel noncoding transcripts, and the significantly upregulated genes were involved in multiple processes in G423. In particular, the expression of genes involved in cell wall and membrane biosynthesis was obviously perturbed under acidic stress. Quantitative real-time PCR analysis showed that the transcription of noncoding small RNA NC-1 increased by 2.35-fold at pH 3.0 compared with that of the control (pH 7.0). Overexpression assays indicated that small RNA NC-1 could significantly enhance the acid tolerance and nisin production of G423 and F44. Our work provided new insights into the sophisticated genetic mechanisms involved in Lc. lactis in an acidic environment, which might elucidate its potential application in food and dairy industries.

Further Elucidation of Galactose Utilization in Lactococcus lactis MG1363

Since the 1970s, galactose metabolism in Lactococcus lactis has been in debate. Different studies led to diverse outcomes making it difficult to conclude whether galactose uptake was PEP- or ATP- dependent and decide what the exact connection was between galactose and lactose uptake and metabolism. It was shown that some Lactococcus strains possess two galactose-specific systems – a permease and a PTS, even if they lack the lactose utilization plasmid, proving that a lactose-independent PTS^Gal exists. However, the PTS^Gal transporter was never identified. Here, with the help of transcriptome analyses and genetic knock-out mutants, we reveal the identities of two low-affinity galactose PTSs. A novel plant-niche-related PTS component Llmg_0963 forming a hybrid transporter Llmg_0963PtcBA and a glucose/mannose-specific PTS are shown to be involved in galactose transport in L. lactis MG1363.

Antisense transcription in Pseudomonas aeruginosa

A large number of antisense transcripts have been detected in diverse microbial genomes and considerable effort has been devoted to elucidating the functional role of antisense transcription. In this study, we reanalysed extensive RNA sequencing data from the opportunistic pathogen Pseudomonas aeruginosa and found that the majority of genes have a propensity for antisense transcription. Although antisense transcripts were found in more than 80 % of the genes of the P. aeruginosa genome, the majority of sequencing reads were mapping sense and only a minority (<2 %) were mapping antisense to genes. Similarly to the sense expression levels, the antisense expression levels varied under different environmental conditions, with the sense and antisense expression levels often being inversely regulated and modulated by the activity of alternative sigma factors. Environment-modulated antisense transcription showed a bias towards being antisense to genes within regions of genomic plasticity and to those encoding small regulatory RNAs. In the future, the validation and functional characterization of antisense transcripts, and novel transcripts that are antisense to small regulatory RNAs in particular, have the potential to contribute to our understanding of the various levels of transcriptional regulation and its dynamics in the bacterial pathogen P. aeruginosa.

The Evolution of gene regulation research in Lactococcus lactis

Lactococcus lactis is a major microbe. This lactic acid bacterium (LAB) is used worldwide in the production of safe, healthy, tasteful and nutritious milk fermentation products. Its huge industrial importance has led to an explosion of research on the organism, particularly since the early 1970s. The upsurge in the research on L. lactis coincided not accidentally with the advent of recombinant DNA technology in these years. The development of methods to take out and re-introduce DNA in L. lactis, to clone genes and to mutate the chromosome in a targeted way, to control (over)expression of proteins and, ultimately, the availability of the nucleotide sequence of its genome and the use of that information in transcriptomics and proteomics research have enabled to peek deep into the functioning of the organism. Among many other things, this has provided an unprecedented view of the major gene regulatory pathways involved in nitrogen and carbon metabolism and their overlap, and has led to the blossoming of the field of L. lactis systems biology. All of these advances have made L. lactis the paradigm of the LAB. This review will deal with the exciting path along which the research on the genetics of and gene regulation in L. lactis has trodden.

Enhancing nisin yield by engineering a small noncodding RNA anti41 and inhibiting the expression of glnR in Lactococcus lactis F44

OBJECTIVES

To engineer a small nonconding RNA anti41 to enhance nisin yield by inhibiting the expression of glnR in Lactococcus lactis F44.

RESULTS

We constructed a screening library to determine appropriate artificial sRNAs and obtained a sRNA anti41 that can produce approximately three fold of the inhibitory effect on GlnR. Moreover, the transcription levels of the direct inhibitory targets of GlnR (glnP, glnQ, amtB, and glnK) were dramatically upregulated in the anti41 overexpression strain (F44-anti41), thereby confirming the inhibitory effect of anti41 on GlnR. In addition, anti41 overexpression improved the survival rate of cells by approximately three fold under acid stress, promoted cell growth, and increased nisin yield by 29.83%.

CONCLUSIONS

We were able to provide a novel strategy for the construction of robust high-producing industrial strains.

Dynamic sporulation gene co-expression networks for Bacillus subtilis 168 and the food-borne isolate Bacillus amyloliquefaciens: a transcriptomic model

Sporulation is a survival strategy, adapted by bacterial cells in response to harsh environmental adversities. The adaptation potential differs between strains and the variations may arise from differences in gene regulation. Gene networks are a valuable way of studying such regulation processes and establishing associations between genes. We reconstructed and compared sporulation gene co-expression networks (GCNs) of the model laboratory strain Bacillus subtilis 168 and the food-borne industrial isolate Bacillus amyloliquefaciens. Transcriptome data obtained from samples of six stages during the sporulation process were used for network inference. Subsequently, a gene set enrichment analysis was performed to compare the reconstructed GCNs of B. subtilis 168 and B. amyloliquefaciens with respect to biological functions, which showed the enriched modules with coherent functional groups associated with sporulation. On basis of the GCNs and time-evolution of differentially expressed genes, we could identify novel candidate genes strongly associated with sporulation in B. subtilis 168 and B. amyloliquefaciens. The GCNs offer a framework for exploring transcription factors, their targets, and co-expressed genes during sporulation. Furthermore, the methodology described here can conveniently be applied to other species or biological processes.

Disruption of a Transcriptional Repressor by an Insertion Sequence Element Integration Leads to Activation of a Novel Silent Cellobiose Transporter in Lactococcus lactis MG1363

Lactococcus lactis subsp. cremoris strains typically carry many dairy niche-specific adaptations. During adaptation to the milk environment these former plant strains have acquired various pseudogenes and insertion sequence elements indicative of ongoing genome decay and frequent transposition events in their genomes. Here we describe the reactivation of a silenced plant sugar utilization cluster in an L. lactis MG1363 derivative lacking the two main cellobiose transporters, PtcBA-CelB and PtcBAC, upon applying selection pressure to utilize cellobiose. A disruption of the transcriptional repressor gene llmg_1239 by an insertion sequence (IS) element allows expression of the otherwise silent novel cellobiose transporter Llmg_1244 and leads to growth of mutant strains on cellobiose. Llmg_1239 was labeled CclR, for cellobiose cluster repressor.IMPORTANCE Insertion sequences (ISs) play an important role in the evolution of lactococci and other bacteria. They facilitate DNA rearrangements and are responsible for creation of new genetic variants with selective advantages under certain environmental conditions. L. lactis MG1363 possesses 71 copies in a total of 11 different types of IS elements. This study describes yet another example of an IS-mediated adaptive evolution. An integration of IS981 or IS905 into a gene coding for a transcriptional repressor led to activation of the repressed gene cluster coding for a plant sugar utilization pathway. The expression of the gene cluster allowed assembly of a novel cellobiose-specific transporter and led to cell growth on cellobiose.

Early Transcriptome Response of Lactococcus lactis to Environmental Stresses Reveals Differentially Expressed Small Regulatory RNAs and tRNAs

Bacteria can deploy various mechanisms to combat environmental stresses. Many genes have previously been identified in Lactococcus lactis that are involved in sensing the stressors and those that are involved in regulating and mounting a defense against the stressful conditions. However, the expression of small regulatory RNAs (sRNAs) during industrially relevant stress conditions has not been assessed yet in L. lactis, while sRNAs have been shown to be involved in many stress responses in other bacteria. We have previously reported the presence of hundreds of putative regulatory RNAs in L. lactis, and have used high-throughput RNA sequencing (RNA-seq) in this study to assess their expression under six different stress conditions. The uniformly designed experimental set-up enabled a highly reliable comparison between the different stress responses and revealed that many sRNAs are differentially expressed under the conditions applied. The primary stress responses of L. lactis NCDO712 was benchmarked to earlier work and, for the first time, the differential expression was assessed of transfer RNAs (tRNAs) and the genes from the six recently sequenced plasmids of NCDO712. Although, we only applied stresses for 5 min, the majority of the well-known specific stress-induced genes are already differentially expressed. We find that most tRNAs decrease after all stresses applied, except for a small number, which are increased upon cold stress. Starvation was shown to induce the highest differential response, both in terms of number and expression level of genes. Our data pinpoints many novel stress-related uncharacterized genes and sRNAs, which calls for further assessment of their molecular and cellular function. These insights furthermore could impact the way parameters are designed for bacterial culture production and milk fermentation, as we find that very short stress conditions already greatly alter gene expression.

Construction of a fur null mutant and RNA-sequencing provide deeper global understanding of the Aliivibrio salmonicida Fur regulon

BACKGROUND

The ferric uptake regulator (Fur) is a transcription factor and the main regulator of iron acquisition in prokaryotes. When bound to ferric iron, Fur recognizes its DNA binding site and generally executes its function by repressing transcription of its target genes. Due to its importance in virulence, the Fur regulon is well studied for several model bacteria. In our previous work, we used computational predictions and microarray to gain insights into Fur-regulation in Aliivibrio salmonicida, and have identified a number of genes and operons that appear to be under direct control of Fur. To provide a more accurate and deeper global understanding of the biological role of Fur we have now generated an A. salmonicida fur knock-out strain and used RNA-sequencing to compare gene expression between the wild-type and fur null mutant strains.

RESULTS

An A. salmonicida fur null mutant strain was constructed. Biological assays demonstrate that deletion of fur results in loss of fitness, with reduced growth rates, and reduced abilities to withstand low-iron conditions, and oxidative stress. When comparing expression levels in the wild-type and the fur null mutant we retrieved 296 differentially expressed genes distributed among 18 of 21 functional classes of genes. A gene cluster encoding biosynthesis of the siderophore bisucaberin represented the highest up-regulated genes in the fur null mutant. Other highly up-regulated genes all encode proteins important for iron acquisition. Potential targets for the RyhB sRNA was predicted from the list of down-regulated genes, and significant complementarities were found between RyhB and mRNAs of the fur, sodB, cysN and VSAL_I0422 genes. Other sRNAs with potential functions in iron homeostasis were identified.

CONCLUSION

The present work provides by far the most comprehensive and deepest understanding of the Fur regulon in A. salmonicida to date. Our data also contribute to a better understanding of how Fur plays a key role in iron homeostasis in bacteria in general, and help to show how Fur orchestrates iron uptake when iron levels are extremely low.

The novel sRNA s015 improves nisin yield by increasing acid tolerance of Lactococcus lactis F44

Nisin, a polycyclic antibacterial peptide produced by Lactococcus lactis, is stable at low pH. Improving the acid tolerance of L. lactis could thus enhance nisin yield. Small non-coding RNAs (sRNAs) play essential roles in acid tolerance by regulating their target mRNAs at the post-transcriptional level. In this study, a novel sRNA, s015, was identified in L. lactis F44 via the use of RNA sequencing, qRT-PCR analysis, and Northern blotting. s015 improved the acid tolerance of L. lactis and boosted nisin yield at low pH. In silico predictions enabled us to construct a library of possible s015 target mRNAs. Statistical analysis and validation suggested that s015 contains a highly conserved region (5'-GAAAAAAAC-3') that likely encompasses the regulatory core of the sRNA. atpG, busAB, cysD, ilvB, tcsR, ung, yudD, and ywdA were verified as direct targets of s015, and the interactions between s015 and its target genes were elucidated. This work provided new insight into the adaptation mechanism of L. lactis under acid stress.

Promoting acid resistance and nisin yield of Lactococcus lactis F44 by genetically increasing D-Asp amidation level inside cell wall

Nisin fermentation by Lactococcus lactis requires a low pH to maintain a relatively higher nisin activity. However, the acidic environment will result in cell arrest, and eventually decrease the relative nisin production. Hence, constructing an acid-resistant L. lactis is crucial for nisin harvest in acidic nisin fermentation. In this paper, the first discovery of the relationship between D-Asp amidation-associated gene (asnH) and acid resistance was reported. Overexpression of asnH in L. lactis F44 (F44A) resulted in a sevenfold increase in survival capacity during acid shift (pH 3) and enhanced nisin desorption capacity compared to F44 (wild type), which subsequently contributed to higher nisin production, reaching 5346 IU/mL, 57.0% more than that of F44 in the fed-batch fermentation. Furthermore, the engineered F44A showed a moderate increase in D-Asp amidation level (from 82 to 92%) compared to F44. The concomitant decrease of the negative charge inside the cell wall was detected by a newly developed method based on the nisin adsorption amount onto cell surface. Meanwhile, peptidoglycan cross-linkage increased from 36.8% (F44) to 41.9% (F44A), and intracellular pH can be better maintained by blocking extracellular H⁺ due to the maintenance of peptidoglycan integrity, which probably resulted from the action of inhibiting hydrolases activity. The inference was further supported by the acmC-overexpression strain F44C, which was characterized by uncontrolled peptidoglycan hydrolase activity. Our results provided a novel strategy for enhancing nisin yield through cell wall remodeling, which contributed to both continuous nisin synthesis and less nisin adsorption in acidic fermentation (dual enhancement).

Spore Heat Activation Requirements and Germination Responses Correlate with Sequences of Germinant Receptors and with the Presence of a Specific spoVA2mob Operon in Foodborne Strains of Bacillus subtilis

Spore heat resistance, germination, and outgrowth are problematic bacterial properties compromising food safety and quality. Large interstrain variation in these properties makes prediction and control of spore behavior challenging. High-level heat resistance and slow germination of spores of some natural Bacillus subtilis isolates, encountered in foods, have been attributed to the occurrence of the spoVA^2mob operon carried on the Tn1546 transposon. In this study, we further investigate the correlation between the presence of this operon in high-level-heat-resistant spores and their germination efficiencies before and after exposure to various sublethal heat treatments (heat activation, or HA), which are known to significantly improve spore responses to nutrient germinants. We show that high-level-heat-resistant spores harboring spoVA^2mob required higher HA temperatures for efficient germination than spores lacking spoVA^2mob The optimal spore HA requirements additionally depended on the nutrients used to trigger germination, l-alanine (l-Ala), or a mixture of l-asparagine, d-glucose, d-fructose, and K⁺ (AGFK). The distinct HA requirements of these two spore germination pathways are likely related to differences in properties of specific germinant receptors. Moreover, spores that germinated inefficiently in AGFK contained specific changes in sequences of the GerB and GerK germinant receptors, which are involved in this germination response. In contrast, no relation was found between transcription levels of main germination genes and spore germination phenotypes. The findings presented in this study have great implications for practices in the food industry, where heat treatments are commonly used to inactivate pathogenic and spoilage microbes, including bacterial spore formers.IMPORTANCE This study describes a strong variation in spore germination capacities and requirements for a heat activation treatment, i.e., an exposure to sublethal heat that increases spore responsiveness to nutrient germination triggers, among 17 strains of B. subtilis, including 9 isolates from spoiled food products. Spores of industrial foodborne isolates exhibited, on average, less efficient and slower germination responses and required more severe heat activation than spores from other sources. High heat activation requirements and inefficient, slow germination correlated with elevated resistance of spores to heat and with specific genetic features, indicating a common genetic basis of these three phenotypic traits. Clearly, interstrain variation and numerous factors that shape spore germination behavior challenge standardization of methods to recover highly heat-resistant spores from the environment and have an impact on the efficacy of preservation techniques used by the food industry to control spores.

If it transcribes, we can sequence it: mining the complexities of host-pathogen-environment interactions using RNA-seq

Host-pathogen interactions are exceedingly complex because they involve multiple host tissues, often occur in the context of normal microflora, and can span diverse microenvironments. Although decades of gene expression studies have provided detailed insights into infection processes, technical challenges have restricted experiments to single pathogenic species or host tissues. RNA-sequencing (RNA-seq) has revolutionized the study of gene expression because in addition to quantifying transcriptional output, it allows detection and characterization of all transcripts in a genome. Here, we review how refined approaches to RNA-seq are used to map the transcriptional networks that control host-pathogen interactions. These enhanced techniques include dRNA-seq and term-seq for the fine-scale mapping of transcriptional start and termination sites, and dual RNA-seq for simultaneous sequencing of host and bacterial pathogen transcriptomes. Dual RNA-seq experiments are currently limited to in vitro infection systems that do not fully reflect the complexities of the in vivo environment, thus a challenge is to develop in vivo model systems and experimental approaches that address the biological heterogeneity of host environments, followed by the integration of RNA-seq with other genome-scale datasets to identify the transcriptional networks that mediate host-pathogen interactions.

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

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

Stress Physiology of Lactic Acid Bacteria

Lactic acid bacteria (LAB) are important starter, commensal, or pathogenic microorganisms. The stress physiology of LAB has been studied in depth for over 2 decades, fueled mostly by the technological implications of LAB robustness in the food industry. Survival of probiotic LAB in the host and the potential relatedness of LAB virulence to their stress resilience have intensified interest in the field. Thus, a wealth of information concerning stress responses exists today for strains as diverse as starter (e.g., Lactococcus lactis), probiotic (e.g., several Lactobacillus spp.), and pathogenic (e.g., Enterococcus and Streptococcus spp.) LAB. Here we present the state of the art for LAB stress behavior. We describe the multitude of stresses that LAB are confronted with, and we present the experimental context used to study the stress responses of LAB, focusing on adaptation, habituation, and cross-protection as well as on self-induced multistress resistance in stationary phase, biofilms, and dormancy. We also consider stress responses at the population and single-cell levels. Subsequently, we concentrate on the stress defense mechanisms that have been reported to date, grouping them according to their direct participation in preserving cell energy, defending macromolecules, and protecting the cell envelope. Stress-induced responses of probiotic LAB and commensal/pathogenic LAB are highlighted separately due to the complexity of the peculiar multistress conditions to which these bacteria are subjected in their hosts. Induction of prophages under environmental stresses is then discussed. Finally, we present systems-based strategies to characterize the "stressome" of LAB and to engineer new food-related and probiotic LAB with improved stress tolerance.