Amplicon-Based High-Throughput Sequencing Method Capable of Species-Level Identification of Coagulase-Negative Staphylococci in Diverse Communities

Long-read metagenomics gives a more accurate insight into the microbiota of long-ripened gouda cheeses

Metagenomic studies of the Gouda cheese microbiota and starter cultures are scarce. During the present study, short-read metagenomic sequencing (Illumina) was applied on 89 Gouda cheese and processed milk samples, which have been investigated before concerning their metabolite and taxonomic composition, the latter applying amplicon-based, high-throughput sequencing (HTS) of the full-length 16S rRNA gene. Selected samples were additionally investigated using long-read metagenomic sequencing (Oxford Nanopore Technologies, ONT). Whereas the species identified by amplicon-based HTS and metagenomic sequencing were identical, the relative abundances of the major species differed significantly. Lactococcus cremoris was more abundant in the metagenomics-based taxonomic analysis compared to the amplicon-based one, whereas the opposite was true for the non-starter lactic acid bacteria (NSLAB). This discrepancy was related to a higher fragmentation of the lactococcal DNA compared with the DNA of other species when applying ONT. Possibly, a higher fragmentation was linked with a higher percentage of dead or metabolically inactive cells, suggesting that full-length 16S rRNA gene amplicon-based HTS might give a more accurate view on active cells. Further, fungi were not abundantly present in the Gouda cheeses examined, whereas about 2% of the metagenomic sequence reads was related to phages, with higher relative abundances in the cheese rinds and long-ripened cheeses. Intraspecies differences found by short-read metagenomic sequencing were in agreement with the amplicon sequence variants obtained previously, confirming the ability of full-length 16S rRNA gene amplicon-based HTS to reach a taxonomic assignment below species level. Metagenome-assembled genomes (MAGs) were retrieved for 15 species, among which the starter cultures Lc. cremoris and Lactococcus lactis and the NSLAB Lacticaseibacillus paracasei, Loigolactobacillus rennini, and Tetragenococcus halophilus, although obtaining MAGs from Lc. cremoris and Lc. lactis was more challenging because of a high intraspecies diversity and high similarity between these species. Long-read metagenomic sequencing could not improve the retrieval of lactococcal MAGs, but, overall, MAGs obtained by long-read metagenomic sequencing solely were superior compared with those obtained by short-read metagenomic sequencing solely, reaching a high-quality draft status of the genomes.

A tool to automatically design multiplex PCR primer pairs for specific targets using diverse templates

Multiplex PCR is an increasingly popular method for identifying species, investigating environmental diversity, and conducting phylogenetic analysis. The complexity and increasing availability of diverse templates necessitate a highly automated approach to design degenerate primer pairs for specific targets with multiple sequences. Existing tools for degenerate primer design suffer from poor maintenance, semi-automation, low adaptability, and low tolerance for gaps. We developed PMPrimer, a Python-based tool for automated design and evaluation of multiplex PCR primer pairs for specific targets using diverse templates. PMPrimer automatically designs optimal multiplex PCR primer pairs using a statistical-based template filter; performs multiple sequence alignment, conserved region identification, and primer design; and evaluates the primers based on template coverage, taxon specificity, and target specificity. PMPrimer identifies conserved regions using Shannon's entropy method, tolerates gaps using a haplotype-based method, and evaluates multiplex PCR primer pairs based on template coverage and taxon specificity. We tested PMPrimer using datasets with diverse levels of conservation, sizes, and applications, including tuf genes of Staphylococci, hsp65 genes of Mycobacteriaceae, and 16S ribosomal RNA genes of Archaea. PMPrimer showed outstanding performance compared with existing tools and experimental validated primers. PMPrimer is available as a Python package at https://github.com/AGIScuipeng/PMPrimer .

Challenge tests reveal limited outgrowth of proteolytic Clostridium botulinum during the production of nitrate- and nitrite-free fermented sausages

Nitrate and nitrite salts perform a versatile role in fermented meats, including the inhibition of food pathogens (in particular proteolytic group I Clostridium botulinum). Despite the increasing interest in clean-label products, little is known about the behaviour of this pathogen in response to the removal of chemical preservatives from fermented meat formulations. Therefore, challenge tests with a cocktail of nontoxigenic group I C. botulinum strains were performed to produce nitrate/nitrite-free fermented sausages under different acidification conditions and starter culture formulations, including the use of an anticlostridial Mammaliicoccus sciuri strain. Results showed limited outgrowth of C. botulinum, even in the absence of acidification. The anticlostridial starter culture did not lead to an additional inhibitory effect. The selective plating procedure adopted within this study proofed robust to follow germination and growth of C. botulinum, inhibiting common fermentative meat microbiota. The challenge tests constitute a suitable tool to assess the behaviour of this food pathogen within fermented meats upon nitrate- and nitrite omission.

Bacterial species diversity of traditionally ripened sheep legs from the Faroe Islands (skerpikjøt)

Skerpikjøt is a traditionally ripened sheep leg product from the Faroe Islands, constituting a relatively underexplored microbial ecosystem. The objective of this study is to achieve a deeper understanding of the microbial composition of this artisanal product. Nine ripened hind legs, obtained from three different producers, were assessed regarding their bacterial communities and contents of biogenic amines, including both surface and core samples. Biogenic amine concentrations were generally low, although one sample had a somewhat elevated concentration of cadaverine. Bacterial diversity was investigated by culture-dependent and culture-independent techniques. Gram-positive catalase-positive cocci (GCC) constituted the most abundant group. Within this group, Staphylococcus equorum was the most prevailing species, followed by Kocuria sp., Mammaliicoccus vitulinus, and Staphylococcus saprophyticus. Lactic acid bacteria prevailed in only one sample and were mainly represented by Latilactobacillus curvatus. Enterobacterial communities were characterised by the prevalence of Serratia proteamaculans. Despite the majority of GCC, Clostridium putrefaciens was the most abundant bacterial species in some core samples. Taken together, the culture-dependent and culture-independent identification methods gave complementary results.

Genome-based characterization of a plasmid-associated micrococcin P1 biosynthetic gene cluster and virulence factors in Mammaliicoccus sciuri IMDO-S72

Analysis of the de novo assembled genome of Mammaliicoccus sciuri IMDO-S72 revealed the genetically encoded machinery behind its earlier reported antibacterial phenotype and gave further insight into the repertoire of putative virulence factors of this recently reclassified species. A plasmid-encoded biosynthetic gene cluster was held responsible for the antimicrobial activity of M. sciuri IMDO-S72, comprising genes involved in thiopeptide production. The compound encoded by this gene cluster was structurally identified as micrococcin P1. Further examination of its genome highlighted the ubiquitous presence of innate virulence factors mainly involved in surface colonization. Determinants contributing to aggressive virulence were generally absent, with exception of a plasmid-associated ica cluster. The native antibiotic resistance genes sal(A) and mecA were detected within the genome, amongst others, but were not consistently linked with a resistant phenotype. While mobile genetic elements were identified within the genome, such as an untypeable SCC element, they proved to be generally free of virulence- and antibiotic-related genes. These results further suggest a commensal lifestyle of M. sciuri and indicate the association of antibiotic resistance determinants with mobile genetic elements, as an important factor in conferring antibiotic resistance, in addition to their unilateral annotation. Importance Mammaliicoccus sciuri has been put forward as an important carrier of virulence and antibiotic resistance genes, which can be transmitted to clinically important staphylococcal species such as Staphylococcus aureus. As a common inhabitant of mammal skin, this species is believed to have a predominant commensal lifestyle although it has been reported as an opportunistic pathogen in some cases. This study provides an extensive genome-wide description of its putative virulence potential taking into consideration the genomic context in which these genes appear, an aspect that is often overlooked during virulence analysis. Additional genome and biochemical analysis linked M. sciuri with the production of micrococcin P1, gaining further insight to which extent these biosynthetic gene cluster are distributed amongst different related species. The frequent plasmid-associated character hints that these traits can be horizontally transferred and might confer a competitive advantage to its recipient within its ecological niche.

Molecular Screening of Microorganisms Associated with Discolored Wood in Dead European Beech Trees Suffered from Extreme Drought Event Using Next Generation Sequencing

Drought events weaken trees and make them vulnerable to attacks by diverse plant pathogens. Here, we propose a molecular method for fast screening of microorganisms associated with European beech decline after an extreme drought period (2018) in a forest of Thuringia, Germany. We used Illumina sequencing with a recent bioinformatics approach based on DADA2 to identify archaeal, bacterial, and fungal ASVs (amplicon sequence variants) based on bacterial and archaeal 16S and fungal ITS genes. We show that symptomatic beech trees are associated with both bacterial and fungal plant pathogens. Although the plant pathogen sequences were detected in both discolored and non-discolored wood areas, they were highly enriched in the discolored wood areas. We show that almost each individual tree was associated with a different combination of pathogens. Cytospora spp. and Neonectria coccinea were among the most frequently detected fungal pathogens, whereas Erwinia spp. and Pseudomonas spp. were the dominant bacterial plant pathogens. We demonstrate that bacterial plant pathogens may be of major importance in beech decline.

Metabolic potential of microbial community and distribution mechanism of Staphylococcus species during broad bean paste fermentation

Although the microbial diversity and structure in bean-based fermented foods have been widely studied, systematic studies on functional microbiota and mechanism of community forms in multi-microbial fermentation systems were still lacking. In this work, the metabolic pathway and functional potential of microbial community in broad bean paste (BBP) were investigated by metagenomics approach, and Staphylococcus, Bacillus, Weissella, Aspergillus and Zygosaccharomyces were found to be the potential predominant populations responsible for substrate alteration and flavor biosynthesis. Among them, Staphylococcus was the most abundant and widespread functional microbe, and closely related Staphylococcus species were diverse and ubiquitously distributed, with the opportunistic pathogen S. gallinarum being the most abundant Staphylococcus specie isolated from BBP. To explain the dominance status of S. gallinarum and species distributions of Staphylococcus genus, we tested the effects of abiotic and biotic factors on three Staphylococcus species using a tractable BBP model, demonstrating that adaptation to environmental conditions (environmental parameters and other functional microbes) led to the dominant position and species coexistence of Staphylococcus, and congeneric competition among Staphylococcus species further shaped ecological distributions of closely related Staphylococcus species. In general, this work revealed the metabolic potential of microbial community and distribution mechanism of Staphylococcus species during BBP fermentation, which could help traditional factories to more precisely control the safety and quality of bean-based fermented foods.

Comparison of three amplicon sequencing approaches to determine staphylococcal populations on human skin

Background

Staphylococci are important members of the human skin microbiome. Many staphylococcal species and strains are commensals of the healthy skin microbiota, while few play essential roles in skin diseases such as atopic dermatitis. To study the involvement of staphylococci in health and disease, it is essential to determine staphylococcal populations in skin samples beyond the genus and species level. Culture-independent approaches such as amplicon next-generation sequencing (NGS) are time- and cost-effective options. However, their suitability depends on the power of resolution.

Results

Here we compare three amplicon NGS schemes that rely on different targets within the genes tuf and rpsK, designated tuf1, tuf2 and rpsK schemes. The schemes were tested on mock communities and on human skin samples. To obtain skin samples and build mock communities, skin swab samples of healthy volunteers were taken. In total, 254 staphylococcal strains were isolated and identified to the species level by MALDI-TOF mass spectrometry. A subset of ten strains belonging to different staphylococcal species were genome-sequenced. Two mock communities with nine and eighteen strains, respectively, as well as eight randomly selected skin samples were analysed with the three amplicon NGS methods. Our results imply that all three methods are suitable for species-level determination of staphylococcal populations. However, the novel tuf2-NGS scheme was superior in resolution power. It unambiguously allowed identification of Staphylococcus saccharolyticus and distinguish phylogenetically distinct clusters of Staphylococcus epidermidis.

Conclusions

Powerful amplicon NGS approaches for the detection and relative quantification of staphylococci in human samples exist that can resolve populations to the species and, to some extent, to the subspecies level. Our study highlights strengths, weaknesses and pitfalls of three currently available amplicon NGS approaches to determine staphylococcal populations. Applied to the analysis of healthy and diseased skin, these approaches can be useful to attribute host-beneficial and -detrimental roles to skin-resident staphylococcal species and subspecies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12866-021-02284-1.

High-throughput amplicon sequencing to assess the impact of processing factors on the development of microbial communities during spontaneous meat fermentation

During spontaneous meat fermentation, diverse microbial communities develop over time. These communities consist mainly of lactic acid bacteria (LAB) and coagulase-negative staphylococci (CNS), of which the species composition is influenced by the fermentation temperature and the level of acidification. Recent development and application of amplicon-based high-throughput sequencing (HTS) methods have allowed to gain deeper insights into the microbial communities of fermented meats. The aim of the present study was to investigate the effect of different fermentation temperatures and acidification profiles on the CNS communities during spontaneous fermentation, using a previously developed amplicon-based HTS method targeting both the 16S rRNA and tuf genes. Spontaneous fermentations were performed with five different lots of meat to assess inter-lot variability. The process influence was investigated by fermenting the meat batters for seven days at different fermentation temperatures (23 °C, 30 °C, and 37 °C) and in the absence or presence of added glucose to simulate different acidification levels. Additionally, the results were compared with a starter culture-initiated fermentation process. The data revealed that the fermentation temperature was the most influential processing condition in shaping the microbial communities during spontaneous meat fermentation processes, whereas differences in pH were only responsible for minor shifts in the microbial profiles. Furthermore, the CNS communities showed a great level of variability, which depended on the initial microbial communities present and their competitiveness.

Application of a High-Throughput Amplicon Sequencing Method to Chart the Bacterial Communities that Are Associated with European Fermented Meats from Different Origins

Insight into the microbial species diversity of fermented meats is not only paramount to gain control over quality development, but also to better understand the link with processing technology and geographical origin. To study the composition of the microbial communities, the use of culture-independent methods is increasingly popular but often still suffers from drawbacks, such as a limited taxonomic resolution. This study aimed to apply a previously developed high-throughput amplicon sequencing (HTS) method targeting the 16S rRNA and tuf genes to characterize the bacterial communities in European fermented meats in greater detail. The data obtained broadened the view on the microbial communities that were associated with the various products examined, revealing the presence of previously underreported subdominant species. Moreover, the composition of these communities could be linked to the specificities of individual products, in particular pH, salt content, and geographical origin. In contrast, no clear links were found between the volatile organic compound profiles of the different products and the country of origin, distinct processing conditions, or microbial communities. Future application of the HTS method offers the potential to further unravel complex microbial communities in fermented meats, as well as to assess the impact of different processing conditions on microbial consortia.