A new in vivo model of intestinal colonization using Zophobas morio larvae: testing hyperepidemic ESBL- and carbapenemase-producing Escherichia coli clones.

Finding strategies for decolonizing gut carriers of multidrug-resistant Escherichia coli (MDR-Ec) is a public-health priority. In this context, novel approaches should be validated in preclinical in vivo gut colonization models before being translated to humans. However, the use of mice presents limitations. Here, we used for the first time Zophobas morio larvae to design a new model of intestinal colonization (28-days duration, T28). Three hyperepidemic MDR-Ec producing extended-spectrum ß-lactamases (ESBLs) or carbapenemases were administered via contaminated food to larvae for the first 7 days (T7): Ec-4901.28 (ST131, CTX-M-15), Ec-042 (ST410, OXA-181) and Ec-050 (ST167, NDM-5). Growth curve analyses showed that larvae became rapidly colonized with all strains (T7, ~106-7 CFU/mL), but bacterial load remained high after the removal of contaminated food only in Ec-4901.28 and Ec-042 (T28, ~103-4 CFU/mL). Moreover, larvae receiving a force-feeding treatment with INTESTI bacteriophage cocktail (on T7 and T10 via gauge needle) were decolonized by Ec-4901.28 (INTESTI-susceptible); however, Ec-042 and Ec-050 (INTESTI-resistant) did not. Initial microbiota (before administering contaminated food) was very rich of bacterial genera (e.g., Lactococcus, Enterococcus, Spiroplasma), but patterns were heterogeneous (Shannon diversity index: range 1.1-2.7) and diverse to each other (Bray-Curtis dissimilarity index ≥30%). However, when larvae were challenged with the MDR-Ec with or without administering bacteriophages the microbiota showed a non-significant reduction of the diversity during the 28-day experiments. In conclusion, the Z. morio larvae model promises to be a feasible and high-throughput approach to study novel gut decolonization strategies for MDR-Ec reducing the number of subsequent confirmatory mammalian experiments.

Detection of blaCTX-M and blaDHA genes in stool samples of healthy people: comparison of culture- and shotgun metagenomic-based approaches.

We implemented culture- and shotgun metagenomic sequencing (SMS)-based methods to assess the gut colonization with extended-spectrum cephalosporin-resistant Enterobacterales (ESC-R-Ent) in 42 volunteers. Both methods were performed using native and pre-enriched (broth supplemented with cefuroxime) stools. Native culture screening on CHROMID® ESBL plates resulted in 17 positive samples, whereas the pre-enriched culture (gold-standard) identified 23 carriers. Overall, 26 ESC-R-Ent strains (24 Escherichia coli) were identified: 25 CTX-M and 3 DHA-1 producers (2 co-producing CTX-Ms). Using the SMS on native stool ("native SMS") with thresholds ≥60% for both identity and coverage, only 7 of the 23 pre-enriched culture-positive samples resulted positive for blaCTX-M/blaDHA genes (native SMS reads mapping to blaCTX-M/blaDHAs identified in gold-standard: sensitivity, 59.0%; specificity 100%). Moreover, an average of 31.5 and 24.6 antimicrobial resistance genes (ARGs) were detected in the 23 pre-enriched culture-positive and the 19 negative samples, respectively. When the pre-enriched SMS was implemented, more blaCTX-M/blaDHA genes were detected than in the native assay, including in stools that were pre-enriched culture-negative (pre-enriched SMS reads mapping to blaCTX-M/blaDHAs identified in gold-standard: sensitivity, 78.3%; specificity 75.0%). In addition, the pre-enriched SMS identified on average 38.6 ARGs/sample, whereas for the corresponding native SMS it was 29.4 ARGs/sample. Notably, stools resulting false-negative by using the native SMS had lower concentrations of ESC-R-Ent (average: ~105 vs. ~107 CFU/g) and E. coli classified reads (average: 193,959 vs. 1.45 million) than those of native SMS positive samples. Finally, the detection of blaCTX-M/blaDHA genes was compared with two well-established bioinformatic tools. In conclusion, only the pre-enriched SMS assured detection of most carriers of ESC-R-Ent. However, its performance was not comparable to the pre-enriched culture-based approach.

Evaluation of Phenotypic Tests to Detect Extended-Spectrum ß-Lactamase (ESBL)-Producing Klebsiella oxytoca Complex Strains.

Klebsiella oxytoca complex (KoC) species may overproduce their chromosomal class A OXY ß-lactamases, conferring reduced susceptibility to piperacillin-tazobactam, expanded-spectrum cephalosporins and aztreonam. Moreover, since clavulanate maintains its ability to inhibit these enzymes, the resulting resistance phenotype may falsely resemble the production of acquired extended-spectrum ß-lactamases (ESBLs). In this work, a collection of 44 KoC strains of human and animal origin was characterized with whole-genome sequencing (WGS) and broth microdilution (BMD) susceptibility testing. Comparison of ESBL producers (n = 11; including CTX-M-15 [n = 6] and CTX-M-1 [n = 5] producers) and hyperproducers of OXYs (n = 21) showed certain phenotypic differences: piperacillin-tazobactam (MIC90s: 16 versus >64 µg/mL), cefotaxime (MIC90s: 64 versus 4 µg/mL), ceftazidime (MIC90s: 32 versus 4 µg/mL), cefepime (MIC90s: 8 versus 4 µg/mL) and associated resistance to non-ß-lactams (e.g., trimethoprim-sulfamethoxazole: 90.9% versus 14.3%, respectively). However, a clear phenotype-based distinction between the two groups was difficult. Therefore, we evaluated 10 different inhibitor-based confirmatory tests to allow such categorization. All tests showed a sensitivity of 100%. However, only combination disk tests (CDTs) with cefepime/cefepime-clavulanate and ceftazidime/ceftazidime-clavulanate or the double-disk synergy test (DDST) showed high specificity (100%, 95.5%, and 100%, respectively). All confirmatory tests in BMD or using the MIC gradient strip did not perform well (specificity, ≤87.5%). Of note, ceftazidime/ceftazidime-avibactam tests also exhibited low specificity (CDT, 87.5%; MIC gradient strip, 77.8%). Our results indicate that standard antimicrobial susceptibility profiles can raise some suspicion, but only the use of cefepime/cefepime-clavulanate CDT or DDST can guarantee distinction between ESBL-producing KoC strains and those hyperproducing OXY enzymes.

Intestinal colonization with multidrug-resistant Enterobacterales: screening, epidemiology, clinical impact, and strategies to decolonize carriers.

The clinical impact of infections due to extended-spectrum ß-lactamase (ESBL)- and/or carbapenemase-producing Enterobacterales (Ent) has reached dramatic levels worldwide. Infections due to these multidrug-resistant (MDR) pathogens-especially Escherichia coli and Klebsiella pneumoniae-may originate from a prior asymptomatic intestinal colonization that could also favor transmission to other subjects. It is therefore desirable that gut carriers are rapidly identified to try preventing both the occurrence of serious endogenous infections and potential transmission. Together with the infection prevention and control countermeasures, any strategy capable of effectively eradicating the MDR-Ent from the intestinal tract would be desirable. In this narrative review, we present a summary of the different aspects linked to the intestinal colonization due to MDR-Ent. In particular, culture- and molecular-based screening techniques to identify carriers, data on prevalence and risk factors in different populations, clinical impact, length of colonization, and contribution to transmission in various settings will be overviewed. We will also discuss the standard strategies (selective digestive decontamination, fecal microbiota transplant) and those still in development (bacteriophages, probiotics, microcins, and CRISPR-Cas-based) that might be used to decolonize MDR-Ent carriers.

Complete Genome Sequence of Entomomonas sp. Strain E2T0, Isolated from the Darkling Beetle Zophobas morio Larvae.

Here, we present the complete genome sequence of Entomomonas sp. E2T0, a strain isolated from larvae of the darkling beetle Zophobas morio. The isolate was fully resistant to aztreonam and possessed a novel class D ß-lactamase gene. The 3,325,929-bp genome consists of a chromosome and a 9,996-bp plasmid.

Genome stability during serial subculturing in hyperepidemic multidrug-resistant Klebsiella pneumoniae and Escherichia coli.

OBJECTIVES: Core-genome single nucleotide variant (cgSNV) analysis represents a powerful tool for epidemiological investigations of multidrug-resistant (MDR) bacteria. However, cgSNV thresholds to confirm whether isolates are the same clone are not formally defined. METHODS: We implemented hybrid whole-genome sequencing to study the genomic changes of four MDR isolates belonging to hyperepidemic sequence types (STs) during 20 propagation steps (T20) on MacConkey and CHROMID(R) ESBL plates. The following strains were analyzed: Klebsiella pneumoniae AE-2247421 (OXA-48/NDM-1-producing, ST101), K. pneumoniae MCL-2017-2 (CTX-M-15-producing, ST307), Escherichia coli Ec-042 (OXA-181-producing, ST410), and E. coli Ec-050 (NDM-5-producing, ST167). The genome assembly at T5 and T20 was compared to that at time point zero (T0) and to two reference genomes. RESULTS: At T20, AE-2247421 lost the IncL blaOXA-48-carrying plasmid when grown on CHROMID(R) ESBL plates, while a large fragment encompassing blaNDM-1 was lost from its IncC plasmid when grown on both plates. In contrast, no structural changes were noted for the other three strains. Regarding the cgSNVs, the following results were obtained at T5 and T20 (ranges considering the different agar plates and reference genomes): AE-2247421 (1-8 and 2-12 cgSNVs), MCL-2017-2 (both 1-2 cgSNVs), Ec-042 (both 0 cgSNVs), and Ec-050 (0-6 and 0-9 cgSNVs). CONCLUSION: We showed that structural changes and accumulation of cgSNVs can occur in few propagation steps under laboratory conditions. These changes might also arise in the clinical context in a short time, especially under antibiotics treatment. This phenomenon should be carefully considered because it might affect the final interpretation of epidemiological genomic analyses.

Carbapenemase-producing Klebsiella pneumoniae strains in Switzerland: human and non-human settings may share high-risk clones.

BACKGROUND: The spread of carbapenemase-producing Klebsiella pneumoniae (CP-Kp) strains belonging to high-risk sequence types (STs) is a concern. For Switzerland, national data about the molecular features (especially the STs) of CP-Kp of human origin is not available. In veterinary clinics, ST11 and ST307 blaOXA-48-possessing K. pneumoniae strains have been recently reported. METHODS: We analysed a collection of 285 K. pneumoniae genomes (170 were CP-Kp) isolated in Switzerland from human and non-human sources during 2006-2020. Whole-genome sequencing, core genome phylogenies and public databases were used to present a detailed overview regarding carbapenemases, STs and plasmids. RESULTS: The top five STs were (main carbapenemase gene) ST512 (blaKPC-3), ST258 (blaKPC-2) and ST101 (blaOXA-48), consisting of strains of human origin only, and ST11 (blaOXA-48) and ST307 (blaOXA-48) strains isolated from human, animal and environmental sources. However, during 2016-2020, the main STs for CP-Kp were ST11 (17.6%), ST307 and ST101 (both 14.7%), whereas ST258 (5.9%) and ST512 (4.4%) significantly declined. Most carbapenemase genes were carried on plasmids already described. Core genome analysis revealed that ST11 K. pneumoniae of animal and human origin were closely related, whereas those of ST307 were distant. CONCLUSIONS: We described, for the first time, the features of the CP-Kp circulating in Switzerland in human and non-human settings. Our genomic analysis revealed that the emerging high-risk ST11 and ST307 lineages were often isolated from non-human settings. This study provided a baseline for further whole-genome sequencing-based One-Health surveillance of CP-Kp and emphasized the need for metadata to track dissemination routes between the different settings.

Repatriation of a patient with COVID-19 contributed to the importation of an emerging carbapenemase producer.

OBJECTIVES: Patients hospitalised abroad can become colonised with multidrug-resistant bacteria and import them to their home countries. In this study, we characterised an OXA-484 carbapenemase-producing Escherichia coli strain from a Swiss patient infected by SARS-CoV-2 and repatriated from India. METHODS: At admission to Switzerland (April 2021), the patient undertook a nasopharyngeal swab to search for SARS-CoV-2 and a rectal swab to detect multidrug-resistant bacteria. Both SARS-CoV-2 and E. coli isolates were whole-genome sequenced and analysed for phylogenetic relatedness. RESULTS: The patient was infected with the SARS-CoV-2 B.1.617.2 lineage (VOC Delta), a lineage that began to be reported across Switzerland at that time. He was also colonised with a sequence type 410 (ST410) E. coli strain (L3452210II) producing OXA-484, a single amino acid variant of OXA-181. The blaOXA-484 gene was carried by a 51.5 kb IncX3 plasmid identical to those described in blaOXA-181-harbouring ST410 E. coli strains. Core genome analysis showed that L3452210II was identical (ΔSNV ≤23) to two ST410 OXA-484 producers recently reported in Qatar and Germany, but differed from other ST410 OXA-181 producers reported worldwide. CONCLUSION: The patient was infected by an emerging SARS-CoV-2 variant and also imported an E. coli producing OXA-484, an OXA-48-like carbapenemase not yet reported in Switzerland. The genetic background of L3452210II indicated that blaOXA-484 shared the same plasmid as blaOXA-181, but its bacterial host differed from most of the pandemic OXA-181-producing ST410 strains reported previously. This case description underlines that the COVID-19 crisis can contribute to the worldwide spread of emerging carbapenemase producers.

Antimicrobial-Resistant Escherichia coli Strains and Their Plasmids in People, Poultry, and Chicken Meat in Laos.

Antimicrobial resistant (AMR) Enterobacterales are widely distributed among the healthy population of the Indochinese peninsula, including Laos. However, the local reservoir of these pathogens are currently not known and possible sources such as agricultural settings and food have rarely been analyzed. In this work, we investigated the extended-spectrum cephalosporin- (ESC-) and colistin-resistant Escherichia coli strains (CST-R-Ec) isolated from the gut of local people, feces of poultry, and from chicken meat (60 samples each group) in Laos. Whole-genome sequencing (WGS) analysis based on both short- and long-read sequencing approaches were implemented. The following prevalence of ESC-R-Ec and CST-R-Ec were recorded, respectively: local people (70 and 15%), poultry (20 and 23.3%), and chicken meat (21.7 and 13.3%). Core-genome analysis, coupled with sequence type (ST)/core-genome ST (cgST) definitions, indicated that no common AMR-Ec clones were spreading among the different settings. ESC-R-Ec mostly possessed bla CTX-M-15 and bla CTX-M-55 associated to ISEcp1 or IS26. The majority of CST-R-Ec carried mcr-1 on IncX4, IncI2, IncP1, and IncHI1 plasmids similar or identical to those described worldwide; strains with chromosomal mcr-1 or possessing plasmid-mediated mcr-3 were also found. These results indicate a high prevalence of AMR-Ec in the local population, poultry, and chicken meat. While we did not observe the same clones among the three settings, most of the bla CTX-Ms and mcr-1/-3 were associated with mobile-genetic elements, indicating that horizontal gene transfer may play an important role in the dissemination of AMR-Ec in Laos. More studies should be planned to better understand the extent and dynamics of this phenomenon.

Complete Genome Sequence of a Third- and Fourth-Generation Cephalosporin-Resistant Comamonas kerstersii Isolate.

Here, we report the complete genome sequence of Comamonas kerstersii 3132976, a strain isolated from a human rectal swab sample in Switzerland. The isolate was resistant to third- and fourth-generation cephalosporins and possessed a novel class A ß-lactamase gene. The complete genome is 3,693,404 bp long with a GC content of 59.4%.

Exploring the Global Spread of Klebsiella grimontii Isolates Possessing blaVIM-1 and mcr-9.

The spread of plasmid-mediated carbapenemases within Klebsiella oxytoca is well-documented. In contrast, data concerning the closely related species Klebsiella grimontii are scarce. In fact, despite the recent report of the first blaKPC-2-producing K. grimontii, nothing is known about its clonality and antibiotic resistance patterns. In a retrospective search in our collection, we identified 2 blaVIM-positive K. oxytoca strains. Whole-genome sequencing with both Illumina and Nanopore indicated that our strains actually belonged to K. grimontii and were of sequence type 172 (ST172) and ST189. Moreover, the two strains were associated with 297-kb IncHI2/HI2A-pST1 and 90.6-kb IncFII(Yp) plasmids carrying blaVIM-1 together with mcr-9 and blaVIM-1, respectively. In the IncHI2/HI2A plasmid, blaVIM-1 was located in a class 1 integron (In110), while mcr-9 was associated with the qseC-qseB-like regulatory elements. Overall, this plasmid was shown to be very similar to those carried by other Enterobacterales isolated from food and animal sources (e.g., Salmonella and Enterobacter spp. detected in Germany and Egypt). The IncFII(Yp) plasmid was unique, and its blaVIM-1 region was associated with a rare integron (In1373). Mapping of In1373 indicated a possible origin in Austria from an Enterobacter hormaechei carrying a highly similar plasmid. Core-genome phylogenies indicated that the ST172 K. grimontii belonged to a clone of identical Swedish and Swiss strains (≤15 single nucleotide variants [SNVs] to each other), whereas the ST189 strain was sporadic. Surveillance of carbapenemase-producing K. oxytoca strains should be reinforced to detect and prevent the dissemination of new species belonging to the Klebsiella genus.

A Patient With Multiple Carbapenemase Producers Including an Unusual Citrobacter sedlakii Hosting an IncC bla NDM-1- and armA-carrying Plasmid.

BACKGROUND: Patients colonized with multiple species of carbapenemase-producing Enterobacterales (CPE) are increasingly observed. This phenomenon can be due to the high local prevalence of these pathogens, the presence of important host risk factors, and the great genetic promiscuity of some carbapenemase genes. METHODS: We analyzed 4 CPE (Escherichia coli, Klebsiella pneumoniae, Providencia stuartii, Citrobacter sedlakii), 1 extended-spectrum cephalosporin-resistant K. pneumoniae (ESC-R-Kp), and 1 carbapenemase-producing Acinetobacter baumannii simultaneously isolated from a patient transferred from Macedonia. Susceptibility tests were performed using a microdilution MIC system. The complete genome sequences were obtained by using both short-read and long-read whole-genome sequencing technologies. RESULTS: All CPE presented high-level resistance to all aminoglycosides due to the expression of the armA 16S rRNA methylase. In C. sedlakii and E. coli (ST69), both the carbapenemase bla NDM-1 and armA genes were located on an identical IncC plasmid of type 1a. The K. pneumoniae (ST268) and P. stuartii carried chromosomal bla NDM-1 and bla OXA-48, respectively, while the ESC-R-Kp (ST395) harbored a plasmid-located bla CTX-M-15. In the latter 3 isolates, armA-harboring IncC plasmids similar to plasmids found in C. sedlakii and E. coli were also detected. The A. baumannii strain possessed the bla OXA-40 carbapenemase gene. CONCLUSIONS: The characterization of the genetic organization of IncC-type plasmids harbored by 3 different species from the same patient offered insights into the evolution of these broad-host-range plasmids. Moreover, we characterized here the first complete genome sequence of a carbapenemase-producing C. sedlakii strain, providing a reference for future studies on this rarely reported species.

Travellers returning from the island of Zanzibar colonized with MDR Escherichia coli strains: assessing the impact of local people and other sources.

OBJECTIVES: Many travellers to low-income countries return home colonized at the intestinal level with extended-spectrum cephalosporin-resistant (ESC-R) and/or colistin-resistant (CST-R) Escherichia coli (Ec) strains. However, nothing is known about the local sources responsible for the transmission of these pathogens to the travellers. METHODS: We compared the ESC-R- and CST-R-Ec strains found in the pre- (n = 23) and post-trip (n = 37) rectal swabs of 37 travellers from Switzerland to Zanzibar with those (i) contemporarily isolated from local people, poultry, retailed chicken meat (n = 31), and (ii) from other sources studied in the recent past (n = 47). WGS and core-genome analyses were implemented. RESULTS: Twenty-four travellers returned colonized with ESC-R- (n = 29) and/or CST-R- (n = 8) Ec strains. Almost all ESC-R-Ec were CTX-M-15 producers and belonged to heterogeneous STs/core-genome STs (cgSTs), while mcr-positive strains were not found. Based on the strains' STs/cgSTs, only 20 subjects were colonized with ESC-R- and/or CST-R-Ec that were not present in their gut before the journey. Single nucleotide variant (SNV) analysis showed that three of these 20 travellers carried ESC-R-Ec (ST3489, ST3580, ST361) identical (0-20 SNVs) to those found in local people, chicken meat, or poultry. Three further subjects carried ESC-R-Ec (ST394, ST648, ST5173) identical or highly related (15-55 SNVs) to those previously reported in local people, fish, or water. CONCLUSIONS: This is the first known study comparing the ESC-R- and/or CST-R-Ec strains obtained from travellers and local sources using solid molecular methods. We showed that for at least one-third of the returning travellers the acquired antibiotic-resistant Ec had a corresponding strain among resident people, food, animal and/or environmental sources.

Rapid Increase of CTX-M-Producing Shigella sonnei Isolates in Switzerland Due to Spread of Common Plasmids and International Clones.

The Swiss Centre for Antibiotic Resistance (ANRESIS) has recently noted an increase of extended-spectrum cephalosporin-resistant (ESC-R) Shigella sonnei isolates nationwide (3.8% in 2016 versus 37.5% in 2019). To understand this phenomenon, we analyzed 25 representative isolates (of which 14 were ESC-R) collected in Switzerland during 2016 to 2019. Whole-genome sequencing was achieved using both the Illumina and the Nanopore platforms. Both ESC-R and extended-spectrum cephalosporin-susceptible isolates belonged to sequence type 152 (ST152). The ESC-R isolates carried blaCTX-M-3 in IncI1-pST57 (n = 5), blaCTX-M-15 in IncFII (F2:A-:B-) (n = 5), blaCTX-M-15 in IncI1-pST16, and blaCTX-M-27, blaCTX-M-55, or blaCTX-M-134 in other IncFII plasmids (n = 1 each). Plasmids having the same bla and Inc group exhibited high degrees of genetic identity to each other but also to plasmids previously reported in other Enterobacterales Core-genome analysis showed that there were 4 main clusters, each of which included strains that differed by <58 single nucleotide variants (SNVs) and that consisted of both blaCTX-M-positive and blaCTX-M-negative isolates. Moreover, most isolates belonging to the same cluster shared an identical core-genome sequence type (cgST). For instance, cluster 1 included 4 isolates of cgST113036, of which only 3 harbored the IncI1-pST57 blaCTX-M-3-positive plasmid. The 25 S. sonnei isolates were also subjected to phylogenetic comparison with deposited international strains. As a result, matching isolates (isolates that had the same cgST and that differed by <8 SNVs) have been reported in the United Kingdom, the United States, France, and the Netherlands. Overall, our results suggest that some common S. sonnei clusters can spread between continents and can be imported into other nations after international trips. Such clusters include, in part, isolates that do not possess blaESBL-harboring plasmids, indicating their tendency to acquire them from other Enterobacterales.

Long-read based de novo assembly of low-complexity metagenome samples results in finished genomes and reveals insights into strain diversity and an active phage system.

BACKGROUND: Complete and contiguous genome assemblies greatly improve the quality of subsequent systems-wide functional profiling studies and the ability to gain novel biological insights. While a de novo genome assembly of an isolated bacterial strain is in most cases straightforward, more informative data about co-existing bacteria as well as synergistic and antagonistic effects can be obtained from a direct analysis of microbial communities. However, the complexity of metagenomic samples represents a major challenge. While third generation sequencing technologies have been suggested to enable finished metagenome-assembled genomes, to our knowledge, the complete genome assembly of all dominant strains in a microbiome sample has not been demonstrated. Natural whey starter cultures (NWCs) are used in cheese production and represent low-complexity microbiomes. Previous studies of Swiss Gruyère and selected Italian hard cheeses, mostly based on amplicon metagenomics, concurred that three species generally pre-dominate: Streptococcus thermophilus, Lactobacillus helveticus and Lactobacillus delbrueckii. RESULTS: Two NWCs from Swiss Gruyère producers were subjected to whole metagenome shotgun sequencing using the Pacific Biosciences Sequel and Illumina MiSeq platforms. In addition, longer Oxford Nanopore Technologies MinION reads had to be generated for one to resolve repeat regions. Thereby, we achieved the complete assembly of all dominant bacterial genomes from these low-complexity NWCs, which was corroborated by a 16S rRNA amplicon survey. Moreover, two distinct L. helveticus strains were successfully co-assembled from the same sample. Besides bacterial chromosomes, we could also assemble several bacterial plasmids and phages and a corresponding prophage. Biologically relevant insights were uncovered by linking the plasmids and phages to their respective host genomes using DNA methylation motifs on the plasmids and by matching prokaryotic CRISPR spacers with the corresponding protospacers on the phages. These results could only be achieved by employing long-read sequencing data able to span intragenomic as well as intergenomic repeats. CONCLUSIONS: Here, we demonstrate the feasibility of complete de novo genome assembly of all dominant strains from low-complexity NWCs based on whole metagenomics shotgun sequencing data. This allowed to gain novel biological insights and is a fundamental basis for subsequent systems-wide omics analyses, functional profiling and phenotype to genotype analysis of specific microbial communities.

Population Dynamics of Lactobacillus helveticus in Swiss Gruyère-Type Cheese Manufactured With Natural Whey Cultures.

Lactobacillus helveticus, a ubiquitous bacterial species in natural whey cultures (NWCs) used for Swiss Gruyère cheese production, is considered to have crucial functions for cheese ripening such as enhancing proteolysis. We tracked the diversity and abundance of L. helveticus strains during 6 months of ripening in eight Swiss Gruyère-type cheeses using a culture-independent typing method. The study showed that the L. helveticus population present in NWCs persisted in cheese and demonstrated a stable multi-strain coexistence during cheese ripening. With regard to proteolysis, one of the eight L. helveticus populations exhibited less protein degradation during ripening.

Comparative Genomics of Completely Sequenced Lactobacillus helveticus Genomes Provides Insights into Strain-Specific Genes and Resolves Metagenomics Data Down to the Strain Level.

Although complete genome sequences hold particular value for an accurate description of core genomes, the identification of strain-specific genes, and as the optimal basis for functional genomics studies, they are still largely underrepresented in public repositories. Based on an assessment of the genome assembly complexity for all lactobacilli, we used Pacific Biosciences' long read technology to sequence and de novo assemble the genomes of three Lactobacillus helveticus starter strains, raising the number of completely sequenced strains to 12. The first comparative genomics study for L. helveticus-to our knowledge-identified a core genome of 988 genes and sets of unique, strain-specific genes ranging from about 30 to more than 200 genes. Importantly, the comparison of MiSeq- and PacBio-based assemblies uncovered that not only accessory but also core genes can be missed in incomplete genome assemblies based on short reads. Analysis of the three genomes revealed that a large number of pseudogenes were enriched for functional Gene Ontology categories such as amino acid transmembrane transport and carbohydrate metabolism, which is in line with a reductive genome evolution in the rich natural habitat of L. helveticus. Notably, the functional Clusters of Orthologous Groups of proteins categories "cell wall/membrane biogenesis" and "defense mechanisms" were found to be enriched among the strain-specific genes. A genome mining effort uncovered examples where an experimentally observed phenotype could be linked to the underlying genotype, such as for cell envelope proteinase PrtH3 of strain FAM8627. Another possible link identified for peptidoglycan hydrolases will require further experiments. Of note, strain FAM22155 did not harbor a CRISPR/Cas system; its loss was also observed in other L. helveticus strains and lactobacillus species, thus questioning the value of the CRISPR/Cas system for diagnostic purposes. Importantly, the complete genome sequences proved to be very useful for the analysis of natural whey starter cultures with metagenomics, as a larger percentage of the sequenced reads of these complex mixtures could be unambiguously assigned down to the strain level.

Amplicon Sequencing of the slpH Locus Permits Culture-Independent Strain Typing of Lactobacillus helveticus in Dairy Products.

The advent of massive parallel sequencing technologies has opened up possibilities for the study of the bacterial diversity of ecosystems without the need for enrichment or single strain isolation. By exploiting 78 genome data-sets from Lactobacillus helveticus strains, we found that the slpH locus that encodes a putative surface layer protein displays sufficient genetic heterogeneity to be a suitable target for strain typing. Based on high-throughput slpH gene sequencing and the detection of single-base DNA sequence variations, we established a culture-independent method to assess the biodiversity of the L. helveticus strains present in fermented dairy food. When we applied the method to study the L. helveticus strain composition in 15 natural whey cultures (NWCs) that were collected at different Gruyère, a protected designation of origin (PDO) production facilities, we detected a total of 10 sequence types (STs). In addition, we monitored the development of a three-strain mix in raclette cheese for 17 weeks.