In Swiss-type cheeses, characteristic nut-like and sweet flavor develops during the cheese ripening due to the metabolic activities of cheese microbiota. Temperature changes during warm and cold room ripening, and duration of ripening can significantly change the gene expression of the cheese microbiota, which can affect the flavor formation. In this study, a metagenomic and metatranscriptomic analysis of Swiss-type Maasdam cheese was performed on samples obtained during ripening in the warm and cold rooms. We reconstructed four different bacterial genomes (Lactococcus lactis, Lactobacillus rhamnosus, Lactobacillus helveticus, and Propionibacterium freudenreichii subsp. shermanii strain JS) from the Maasdam cheese to near completeness. Based on the DNA and RNA mean coverage, Lc. lactis strongly dominated (~80-90%) within the cheese microbial community. Genome annotation showed the potential for the presence of several flavor forming pathways in these species, such as production of methanethiol, free fatty acids, acetoin, diacetyl, acetate, ethanol, and propionate. Using the metatranscriptomic data, we showed that, with the exception of Lc. lactis, the central metabolism of the microbiota was downregulated during cold room ripening suggesting that fewer flavor compounds such as acetoin and propionate were produced. In contrast, Lc. lactis genes related to the central metabolism, including the vitamin biosynthesis and homolactic fermentation, were upregulated during cold room ripening.

Background: High-pressure processing (HPP) is a commonly used technique in the food industry to inactivate pathogens, including L. monocytogenes. It has been shown that L. monocytogenes is able to recover from HPP injuries and can start to grow again during long-term cold storage. To date, the gene expression profiling of L. monocytogenes during HPP damage recovery at cooling temperature has not been studied. In order identify key genes that play a role in recovery of the damage caused by HPP treatment, we performed RNA-sequencing (RNA-seq) for two L. monocytogenes strains (barotolerant RO15 and barosensitive ScottA) at nine selected time points (up to 48 h) after treatment with two pressure levels (200 and 400 MPa). Results: The results showed that a general stress response was activated by SigB after HPP treatment. In addition, the phosphotransferase system (PTS; mostly fructose-, mannose-, galactitol-, cellobiose-, and ascorbate-specific PTS systems), protein folding, and cobalamin biosynthesis were the most upregulated genes during HPP damage recovery. We observed that cell-division-related genes (divIC, dicIVA, ftsE, and ftsX) were downregulated. By contrast, peptidoglycan-synthesis genes (murG, murC, and pbp2A) were upregulated. This indicates that cell-wall repair occurs as a part of HPP damage recovery. We also observed that prophage genes, including anti-CRISPR genes, were induced by HPP. Interestingly, a large amount of RNA-seq data (up to 85%) was mapped to Rli47, which is a non-coding RNA that is upregulated after HPP. Thus, we predicted that Rli47 plays a role in HPP damage recovery in L. monocytogenes. Moreover, gene-deletion experiments showed that amongst peptidoglycan biosynthesis genes, pbp2A mutants are more sensitive to HPP. Conclusions: We identified several genes and mechanisms that may play a role in recovery from HPP damage of L. monocytogenes. Our study contributes to new information on pathogen inactivation by HPP.

Background: Psychrotrophic lactic acid bacteria (LAB) species are the dominant species in the microbiota of cold-stored modified-atmosphere-packaged food products and are the main cause of food spoilage. Despite the importance of psychrotrophic LAB, their response to cold or heat has not been studied. Here, we studied the transcriptome-level cold- and heat-shock response of spoilage lactic acid bacteria with time-series RNA-seq for Le. gelidum, Lc. piscium, and P. oligofermentans at 0 °C, 4 °C, 14 °C, 25 °C, and 28 °C. Results: We observed that the cold-shock protein A (cspA) gene was the main cold-shock protein gene in all three species. Our results indicated that DEAD-box RNA helicase genes (cshA, cshB) also play a critical role in cold-shock response in psychrotrophic LAB. In addition, several RNase genes were involved in cold-shock response in Lc. piscium and P. oligofermentans. Moreover, gene network inference analysis provided candidate genes involved in cold-shock response. Ribosomal proteins, tRNA modification, rRNA modification, and ABC and efflux MFS transporter genes clustered with cold-shock response genes in all three species, indicating that these genes could be part of the cold-shock response machinery. Heat-shock treatment caused upregulation of Clp protease and chaperone genes in all three species. We identified transcription binding site motifs for heat-shock response genes in Le. gelidum and Lc. piscium. Finally, we showed that food spoilage-related genes were upregulated at cold temperatures. Conclusions: The results of this study provide new insights on the cold- and heat-shock response of psychrotrophic LAB. In addition, candidate genes involved in cold- and heat-shock response predicted using gene network inference analysis could be used as targets for future studies.

Background: High pressure processing (HPP; i.e. 100–600 MPa pressure depending on product) is a non-thermal preservation technique adopted by the food industry to decrease significantly foodborne pathogens, includingListeria monocytogenes, from food.However, susceptibility towards pressure differs among diverse strains ofL. monocytogenesand it is unclear if this is due to their intrinsic characteristics related to genomic content. Here, we tested the barotolerance of 10 differentL. monocytogenesstrains, from food and food processing environments and widely used reference strains including clinical isolate, to pressure treatments with 400 and 600 MPa. Genome sequencing and genome comparison of the testedL. monocytogenesstrains were performed to investigate the relation between genomic profile and pressure tolerance.Results: None of the tested strains were tolerant to 600 MPa. A reduction of more than 5 log₁₀was observed for all strains after 1 min 600 MPa pressure treatment.L. monocytogenesstrain RO15 showed no significant reduction in viable cell counts after 400 MPa for 1 min and was therefore defined as barotolerant. Genome analysis of so far unsequencedL. monocytogenesstrain RO15, 2HF33, MB5, AB199, AB120, C7, and RO4 allowed us to compare the gene content of all strains tested. This revealed that the three most pressure tolerant strains had more than one CRISPR system with self-targeting spacers. Furthermore, several anti-CRISPR genes were detected in these strains. Pan-genome analysis showed that 10 prophage genes were significantly associated with the three most barotolerant strains.Conclusions: L. monocytogenesstrain RO15 was the most pressure tolerant among the selected strains. Genome comparison suggests that there might be a relationship between prophages and pressure tolerance inL. monocytogenes.

The human gut microbiome composition has been linked to Parkinson’s disease (PD). However, knowledge of the gut microbiota on the genome level is still limited. Here we performed deep metagenomic sequencing and binning to build metagenome-assembled genomes (MAGs) from 136 human fecal microbiomes (68 PD samples and 68 control samples). We constructed 952 non-redundant high-quality MAGs and compared them between PD and control groups. Among these MAGs, there were 22 different genomes of Collinsella and Prevotella, indicating high variability of those genera in the human gut environment. Microdiversity analysis indicated that Ruminococcus bromii was statistically significantly (p < 0.002) more diverse on the strain level in the control samples compared to the PD samples. In addition, by clustering all genes and performing presence-absence analysis between groups, we identified several control-specific (p < 0.05) related genes, such as speF and Fe-S oxidoreductase. We also report detailed annotation of MAGs, including Clusters of Orthologous Genes (COG), Cas operon type, antiviral gene, prophage, and secondary metabolites biosynthetic gene clusters, which can be useful for providing a reference for future studies.

The Glanville fritillary (Melitaea cinxia) butterfly is a model system for metapopulation dynamics research in fragmented landscapes. Here, we provide a chromosome-level assembly of the butterfly's genome produced from Pacific Biosciences sequencing of a pool of males, combined with a linkage map from population crosses. The final assembly size of 484 Mb is an increase of 94 Mb on the previously published genome. Estimation of the completeness of the genome with BUSCO indicates that the genome contains 92–94% of the BUSCO genes in complete and single copies. We predicted 14,810 genes using the MAKER pipeline and manually curated 1,232 of these gene models. The genome and its annotated gene models are a valuable resource for future comparative genomics, molecular biology, transcriptome, and genetics studies on this species.