Whole Genome Sequence Analysis Using JSpecies Tool Establishes Clonal Relationships between Listeria monocytogenes Strains from Epidemiologically Unrelated Listeriosis Outbreaks

Functional Analyses of the Bacillus velezensis HMB26553 Genome Provide Evidence That Its Genes Are Potentially Related to the Promotion of Plant Growth and Prevention of Cotton Rhizoctonia Damping-Off

Bacillus spp. is one kind of the important representative biocontrol agents against plant diseases and promoting plant growth. In this study, the whole genomic sequence of bacterial strain HMB26553 was obtained. A phylogenetic tree based on the genome and ANI (average nucleotide identity), as well as dDDH (digital DNA-DNA hybridization), was constructed, and strain HMB26553 was identified as Bacillus velezensis. Fourteen biosynthetic gene clusters responsible for secondary metabolite were predicted via anti-SMASH, and six secondary metabolites were identified by UHPLC-QTOF-MS/MS (ultra-high-performance liquid chromatography coupled to quadrupole-time-of-flight tandem mass spectrometry). When the phytopathogen Rhizoctonia solani was treated with B. velezensis HMB26553, the mycelial structure changed, ROS (reactive oxygen species) accumulated, and the mitochondrial membrane potential decreased. Characteristics of strain HMB26553 were predicted and confirmed by genomic information and experiments, such as producing IAA, siderophore, extracellular enzymes and biofilm, as well as moving and promoting cotton growth. All these results suggested the mechanisms by which B. velezensis HMB26553 inhibits pathogen growth and promotes cotton growth, which likely provided the potential biocontrol agent to control cotton Rhizoctonia damping-off.

The identification of the new species Nitratireductor thuwali sp. nov. reveals the untapped diversity of hydrocarbon-degrading culturable bacteria from the arid mangrove sediments of the Red Sea

Introduction

The geological isolation, lack of freshwater inputs and specific internal water circulations make the Red Sea one of the most extreme-and unique-oceans on the planet. Its high temperature, salinity and oligotrophy, along with the consistent input of hydrocarbons due to its geology (e.g., deep-sea vents) and high oil tankers traffic, create the conditions that can drive and influence the assembly of unique marine (micro)biomes that evolved to cope with these multiple stressors. We hypothesize that mangrove sediments, as a model-specific marine environment of the Red Sea, act as microbial hotspots/reservoirs of such diversity not yet explored and described.

Methods

To test our hypothesis, we combined oligotrophic media to mimic the Red Sea conditions and hydrocarbons as C-source (i.e., crude oil) with long incubation time to allow the cultivation of slow-growing environmentally (rare or uncommon) relevant bacteria.

Results and discussion

This approach reveals the vast diversity of taxonomically novel microbial hydrocarbon degraders within a collection of a few hundred isolates. Among these isolates, we characterized a novel species, Nitratireductor thuwali sp. nov., namely, Nit1536^T. It is an aerobic, heterotrophic, Gram-stain-negative bacterium with optimum growth at 37°C, 8 pH and 4% NaCl, whose genome and physiological analysis confirmed the adaptation to extreme and oligotrophic conditions of the Red Sea mangrove sediments. For instance, Nit1536^T metabolizes different carbon substrates, including straight-chain alkanes and organic acids, and synthesizes compatible solutes to survive in salty mangrove sediments. Our results showed that the Red Sea represent a source of yet unknown novel hydrocarbon degraders adapted to extreme marine conditions, and their discovery and characterization deserve further effort to unlock their biotechnological potential.

Listeria monocytogenes in foods-From culture identification to whole-genome characteristics

Listeria monocytogenes is an important foodborne pathogen, which is able to persist in the food production environments. The presence of these bacteria in different niches makes them a potential threat for public health. In the present review, the current information on the classical and alternative methods used for isolation and identification of L. monocytogenes in food have been described. Although these techniques are usually simple, standardized, inexpensive, and are routinely used in many food testing laboratories, several alternative molecular-based approaches for the bacteria detection in food and food production environments have been developed. They are characterized by the high sample throughput, a short time of analysis, and cost-effectiveness. However, these methods are important for the routine testing toward the presence and number of L. monocytogenes, but are not suitable for characteristics and typing of the bacterial isolates, which are crucial in the study of listeriosis infections. For these purposes, novel approaches, with a high discriminatory power to genetically distinguish the strains during epidemiological studies, have been developed, e.g., whole-genome sequence-based techniques such as NGS which provide an opportunity to perform comparison between strains of the same species. In the present review, we have shown a short description of the principles of microbiological, alternative, and modern methods of detection of L. monocytogenes in foods and characterization of the isolates for epidemiological purposes. According to our knowledge, similar comprehensive papers on such subject have not been recently published, and we hope that the current review may be interesting for research communities.

Persistence of Listeria monocytogenes ST5 in Ready-to-Eat Food Processing Environment

Most human listeriosis is foodborne, and ready-to-eat (RET) foods contaminated by Listeria monocytogenes during processing are found to be common vehicles. In this study, a total of four L. monocytogens STs (ST5, ST121, ST120, and ST2) have been identified in two RTE food plants from 2019 to 2020 in Shanghai, China. The L. monocytogenes ST5 was predominant in one RTE food processing plant, and it persists in the RTE meat processing plant with continued clone transmission. The genetic features of the four STs isolates were different. ST5 and ST121 had the three genes clpL, mdrL, and lde; however, ST120 and ST2 had two genes except for clpL. SSI-1was present in ST5, ST121, and ST120. Additionally, SSI-2 was present only in the ST121 isolates. ST120 had all six biofilm-forming associated genes (actA, prfA, lmo0673, recO, lmo2504 and luxS). The ST2 isolate had only three biofilm-forming associated genes, which were prfA, lmo0673, and recO. The four ST isolates had different biofilm formation abilities at different stages. The biofilm formation ability of ST120 was significantly higher when grown for one day. However, the biofilm formation ability of ST120 reduced significantly after growing for four days. In contrast, the biofilm formation ability of ST5 and ST121 increased significantly. These results suggested that ST5 and ST121 had stronger ability to adapt to stressful environments. Biofilms formed by all four STs grown over four days can be sanitized entirely by a disinfectant concentration of 500 mg/L. Additionally, only ST5 and ST121 biofilm cells survived in sub-lethal concentrations of chlorine-containing disinfectant. These results suggested that ST5 and ST121 were more resistant to chlorine-containing disinfectants. These results indicated that the biofilm formation ability of L. monocytogenes isolates changed at different stages. Additionally, the persistence in food processing environments might be verified by the biofilm formation, stress resistance, etc. Alternatively, these results underlined that disinfectants should be used at lethal concentrations. More attention should be paid to ST5 and ST121, and stronger surveillance should be taken to prevent and control the clonal spread of L. monocytogenes isolates in food processing plants in Shanghai.

Genetic Diversity and Relationships of Listeria monocytogenes Serogroup IIa Isolated in Poland

In the present study, 100 L. monocytogenes isolates of serogroup IIa from food and food production environments in Poland were characterized towards the presence of virulence, resistance, and stress response genes using whole-genome sequencing (WGS). The strains were also molecularly typed and compared with multi-locus sequence typing (MLST) and core genome MLST analyses. The present isolates were grouped into 6 sublineages (SLs), with the most prevalent SL155 (33 isolates), SL121 (32 isolates), and SL8 (28 isolates) and classified into six clonal complexes, with the most prevalent CC155 (33 strains), CC121 (32 isolates), and CC8 (28 strains). Furthermore, the strains were grouped to eight sequence types, with the most prevalent ST155 (33 strains), ST121 (30 isolates), and ST8 (28; strains) followed by 60 cgMLST types (CTs). WGS data showed the presence of several virulence genes or putative molecular markers playing a role in pathogenesis of listeriosis and involved in survival of L. monocytogenes in adverse environmental conditions. Some of the present strains were molecularly closely related to L. monocytogenes previously isolated in Poland. The results of the study showed that food and food production environments may be a source of L. monocytogenes of serogroup IIa with pathogenic potential.

Insights Into the Cultivable Bacterial Fraction of Sediments From the Red Sea Mangroves and Physiological, Chemotaxonomic, and Genomic Characterization of Mangrovibacillus cuniculi gen. nov., sp. nov., a Novel Member of the Bacillaceae Family

Mangrove forests are dynamic and productive ecosystems rich in microbial diversity; it has been estimated that microbial cells in the mangrove sediments constitute up to 91% of the total living biomass of these ecosystems. Despite in this ecosystem many of the ecological functions and services are supported and/or carried out by microorganisms (e.g., nutrient cycling and eukaryotic-host adaptation), their diversity and function are overlooked and poorly explored, especially for the oligotrophic mangrove of the Red Sea coast. Here, we investigated the cultivable fraction of bacteria associated with the sediments of Saudi Arabian Red Sea mangrove forest by applying the diffusion-chamber-based approach in combination with oligotrophic medium and long incubation time to allow the growth of bacteria in their natural environment. Cultivation resulted in the isolation of numerous representatives of Isoptericola (n = 51) and Marinobacter (n = 38), along with several less abundant and poorly study taxa (n = 25) distributed across ten genera. Within the latest group, we isolated R1DC41^T, a novel member of the Bacillaceae family in the Firmicutes phylum. It showed 16S rRNA gene similarity of 94.59–97.36% with closest relatives of Rossellomorea (which was formerly in the Bacillus genus), Domibacillus, Bacillus, and Jeotgalibacillus genera. Based on the multilocus sequence analysis (MLSA), R1DC41^T strain formed a separated branch from the listed genera, representing a novel species of a new genus for which the name Mangrovibacillus cuniculi gen. nov., sp. nov. is proposed. Genomic, morphological, and physiological characterizations revealed that R1DC41^T is an aerobic, Gram-stain-variable, rod-shaped, non-motile, endospore-forming bacterium. A reduced genome and the presence of numerous transporters used to import the components necessary for its growth and resistance to the stresses imposed by the oligotrophic and salty mangrove sediments make R1DC41^T extremely adapted to its environment of origin and to the competitive conditions present within.

Listeria monocytogenes Contamination Characteristics in Two Ready-to-Eat Meat Plants From 2019 to 2020 in Shanghai

Listeria monocytogenes is a ubiquitous foodborne pathogen that causes listeriosis and is mostly linked to consumption of ready-to-eat (RTE) foods. Lack of hygiene in food processing environments may be a primary reason for contamination by L. monocytogenes isolates. In this study, L. monocytogenes strains isolated from two RTE meat processing plants in the Shanghai municipality, China, were characterized during 2019–2020 using pulsed-field gel electrophoresis and whole-genome sequencing. Results showed that 29 samples (12.2%) out of 239 were positive for L. monocytogenes, with 21 (18.9%) and 8 (6.25%) isolates from plants A and B, respectively. The packaging room at plant A had the most contamination (14, 48.3%; p < 0.05), with a peak occurrence of 76.5% in processing environments. Nineteen L. monocytogenes isolates belonging to the pulsotype (PT) 7 group were indistinguishable (≥ 95.7%). Furthermore, core-genome multiple loci sequencing typing identified up to nine allelic differences, and the closet pairwise differences among these ST5 isolates included 0–16 small nucleotide polymorphisms. Therefore, L. monocytogenes likely persisted at plant A during 2019–2020 with ongoing clone transmission. In contrast, no L. monocytogenes isolates were identified from processing environments at plant B. Most L. monocytogenes isolates were sampled from raw materials (62.5%). Several isolates (ST378, ST8, and ST120) were detected only once in 2020 and were considered as transient isolates. However, three ST121 isolates with the same PT (PT2) were detected in 2020 and should be noted for their stronger survival ability in harsh environments. These results suggest that continuous monitoring, stringent surveillance, and source tracking are crucial to guaranteeing food safety in RTE food plants.

Comparative genomics analysis of strains from diverse sources reveals the evolutionary relationship of Aeromonas veronii

Aeromonas veronii (A. veronii, AV) strains are emerging zoonotic and aquatic pathogens, yet we know very little about their genomics. This study aims to utilize comparative genomics to investigate the intraspecific genetic diversity, differences in virulence factors and evolutionary mechanisms of A. veronii strains from diverse sources and to fundamentally demonstrate their pathogenic mechanisms. We conducted comparative genomics analysis of 39 A. veronii strains from different sources and found that 1993 core genes are shared by these strains and that these shared core genes may be necessary to maintain the basic characteristics of A. veronii. Additionally, phylogenetic relationship analysis based on these shared genes revealed that a distant relationship between the AMC34 strain and the other 38 strains but that, the genetic relationship among the 38 strains is relatively close, indicating that AMC34 may not belong to A. veronii. Furthermore, analysis of shared core genes and average nucleotide identity (ANI) values showed no obvious correlation with the location of A. veronii isolation and genetic relationship. Our research indicates the evolutionary mechanism of A. veronii from different sources and provides new insights for a deeper understanding of its pathogenic mechanism.

Mangrovivirga cuniculi gen. nov., sp. nov., a moderately halophilic bacterium isolated from bioturbated Red Sea mangrove sediment, and proposal of the novel family Mangrovivirgaceae fam. nov

A strictly aerobic, Gram-stain-negative, non-motile, rod-shaped bacterium, designated strain R1DC9^T, was isolated from sediments of a mangrove stand on the Red Sea coast of Saudi Arabia via diffusion chamber cultivation. Strain R1DC9^T grew at 20-40 °C (optimum, 37 °C), pH 6-10 (optimum, pH 8) and 3-11 % NaCl (optimum, 7-9 %) in the cultivation medium. The genome of R1DC9^T was 4 661 901 bp long and featured a G+C content of 63.1 mol%. Phylogenetic analyses based on the 16S rRNA gene sequence and whole-genome multilocus sequence analysis using 120 concatenated single-copy genes revealed that R1DC9^T represents a distinct lineage in the order Cytophagales and the phylum Bacteroidetes separated from the Roseivirgaceae and Marivirgaceae families. R1DC9^T displayed 90 and 89 % 16S rRNA gene sequence identities with Marivirga sericea DSM 4125^T and Roseivirga ehrenbergii KMM 6017^T, respectively. The predominant quinone was MK7. The polar lipids were phosphatidylethanolamine, two unknown phospholipids and two unknown lipids. The predominant cellular fatty acids were the saturated branch chain fatty acids iso-C_15 : 0, iso-C_17 : 0 3-OH and iso-C_17 : 0, along with a low percentage of the monounsaturated fatty acid C_16 : 1  ω5c. Based on differences in phenotypic, physiological and biochemical characteristics from known relatives, and the results of phylogenetic analyses, R1DC9^T (=KCTC 72349^T=JCM 33609^T=NCCB 100698^T) is proposed to represent a novel species in a new genus, and the name Mangrovivirga cuniculi gen. nov., sp. nov. is proposed. The distinct phylogenetic lineage among the families in the order Cytophagales indicates that R1DC9^T represents a new family for which the name Mangrovivirgaceae fam. nov. is proposed.

Genomic Analysis of Prophages Recovered from Listeria monocytogenes Lysogens Found in Seafood and Seafood-Related Environment

A prophage is a phage-related sequence that is integrated into a bacterial chromosome. Prophages play an important role in bacterial evolution, survival, and persistence. To understand the impact of Listeria prophages on their host genome organizations, this work sequenced two L. monocytogenes strains (134LM and 036LM), previously identified as lysogens by mitomycin C induction. Draft genomes were generated with assembly sizes of 2,953,877 bp and 3,000,399 bp. One intact prophage (39,532 bp) was inserted into the comK gene of the 134LM genome. Two intact prophages (48,684 bp and 39,488 bp) were inserted in tRNA-Lys and elongation-factor genes of the 036LM genome. The findings confirmed the presence of three corresponding induced phages previously obtained by mitomycin C induction. Comparative genomic analysis of three prophages obtained in the newly sequenced lysogens with 61 prophages found in L. monocytogenes genomes, available in public databases, identified six major clusters using whole genome-based phylogenetic analysis. The results of the comparative genomic analysis of the prophage sequences provides knowledge about the diversity of Listeria prophages and their distribution among Listeria genomes in diverse environments, including different sources or geographical regions. In addition, the prophage sequences and their insertion sites contribute to the genomic diversity of L. monocytogenes genomes. These data of prophage sequences, prophage insertion sites, and prophage sequence comparisons, together with ANIb confirmation, could be useful for L. monocytogenes classification by prophages. One potential development could be refinement of prophage typing tools for monitoring or surveillance of L. monocytogenes contamination and transmission.

Molecular Typing of Listeria monocytogenes IVb Serogroup Isolated from Food and Food Production Environments in Poland

Listeria monocytogenes is one of the most important foodborne pathogens that may be present in food and in food processing environments. In the present study, 91 L. monocytogenes isolates of serogroup IVb from raw meat, ready-to-eat food and food production environments in Poland were characterized by whole genome sequencing (WGS). The strains were also compared, using core genome multi-locus sequence typing (cgMLST) analysis, with 186 genomes of L. monocytogenes recovered worldwide from food, environments, and from humans with listeriosis. The L. monocytogenes examined belonged to three MLST clonal complexes: CC1 (10; 11.0% isolates), CC2 (70; 76.9%), and CC6 (11; 12.1%). CC1 comprised of two STs (ST1 and ST515) which could be divided into five cgMLST, CC2 covered two STs (ST2 and ST145) with a total of 20 cgMLST types, whereas CC6 consisted of only one ST (ST6) classified as one cgMLST. WGS sequences of the tested strains revealed that they had several pathogenic markers making them potentially hazardous for public health. Molecular comparison of L. monocytogenes strains tested in the present study with those isolated from food and human listeriosis showed a relationship between the isolates from Poland, but not from other countries.

Development and Validation of a Quantitative PCR Method for Species Verification and Serogroup Determination of Listeria monocytogenes Isolates

ABSTRACT

Listeria monocytogenes (Lm) is one of the leading causes of death because of foodborne illness, affecting the elderly, pregnant women, neonates, and people who are immunocompromised. Serologically, Lm can be classified into 13 serotypes, although only 4 are typically linked with food contamination and illness. Since 2000, a shift in serotypes involved in listeriosis outbreaks has been observed, suggesting that tracking of serotypes could help identify emerging trends. A PCR method developed in 2004 allowed detection of the four major serotypes as molecular serogroups, corresponding to broad phylogenetic groups. In this study, a novel quantitative PCR (qPCR) method was developed that uses two multiplex qPCRs, one to confirm the Listeria genus and Lm species and the second for Lm molecular serogrouping. This method was compared with the U.S. Food and Drug Administration Bacteriological Analytical Manual (BAM) method for Lm and the seroagglutination method, using a 208-strain panel. Comparison of the genus and species qPCR assay with the BAM methods found an equal or slightly higher accuracy for the qPCR method (>98%), compared with the BAM protocol (>96%), when evaluated against independent characterization data. Molecular serogrouping using the qPCR method (96.6%) was more accurate than the seroagglutination assay (75.6%). The qPCR method identified Lm 4bV strains, which could not be resolved using seroagglutination. The qPCR could not identify lineage III and IV serotype 4b strains but did correctly identify 16 of 18 lineage III and IV strains. The qPCR method performed genus identification for the Listeria species Lm, L. innocua, L. welshimeri, L. ivanovii, and L. seeligeri. In addition, the method performed species identification for Lm and classified Lm into six molecular serogroups: 2A, 2B, 2C, 4B, NT, and 4bV. This method provided a rapid and accurate confirmation of Lm and serogroup determinations; furthermore, it could help identify otherwise unlinked strains by enabling whole genome sequencing analysis based on broad phylogeny, independent of other information.

HIGHLIGHTS

10-Year Molecular Surveillance of Listeria monocytogenes Using Whole-Genome Sequencing in Shanghai, China, 2009-2019

Listeria monocytogenes is an etiologic agent of listeriosis, and has emerged as an important foodborne pathogen worldwide. In this study, the molecular characteristics of 155 L. monocytogenes isolates from seven food groups in Shanghai, the biggest city in China, were identified using whole-genome sequencing (WGS). Most L. monocytogenes isolates (79.3%) were obtained between May and October from 2009 to 2019. The serogroups and clonal complexes (CCs) of L. monocytogenes were found useful for identifying potential health risks linked to foods. Differences in distributions of serogroups and CCs among different food groups were analyzed using t-test. The results showed that the IIa and IVb serogroups were identified among most of food groups. However, the prevalence of serogroup IIb was significantly higher in ready-to-eat (RTE) food and raw seafood than in other food groups, similar to group IIc in raw meat and raw poultry than others. Meanwhile, the prevalence of CC9 in raw meat and raw poultry, CC8 in raw poultry, and CC87 in raw seafood significantly exceeded that of in other food groups. Specially, CC87 was the predominant CC in foodborne and clinical isolates in China, indicating that raw seafood may induce a high-risk to food safety. Also, hypervirulence pathogenicity islands LIPI-3 and LIPI-4 were found in CC3, CC1, and CC87, respectively. The clonal group CC619 carried LIPI-3 and LIPI-4, as previously reported in China. Core genome multilocus sequence typing (cgMLST) analysis suggested that CC87 isolates from the same food groups in different years had no allelic differences, indicating that L. monocytogenes could persist over years. These 10-year results in Shanghai underscore the significance of molecular epidemiological surveillance of L. monocytogenes in foodborne products in assessing the potential risk of this pathogen, and further address food safety issues in China.

Microevolution and Gain or Loss of Mobile Genetic Elements of Outbreak-Related Listeria monocytogenes in Food Processing Environments Identified by Whole Genome Sequencing Analysis

Whole genome sequencing (WGS) analyses have been instrumental in traceback investigations of Listeria monocytogenes (Lm). To demonstrate how long-read sequencing analysis can capture and describe relationships among isolates from clinical, food, and environmental sources, we analyzed 366 long-read- and shotgun-sequenced isolates from 16 Lm outbreak strains associated with cantaloupe, leafy green, stone fruit, caramel apple, mung bean sprout, multiple cheese products, multiple ice cream products, and their production environments. The analyses demonstrated that outbreak strains could be distributed in different areas and zones of food production environments through persistent or repeated contamination. Multi-strain and multi-clone contamination were common. Further, WGS could differentiate among isolates collected at different time points or from different production lines in the same facility, revealing microevolution events in processing environments. Our comparison between complete and shotgun genomes showed that isolates of the same outbreak strain diversified mostly by gain/loss of plasmids and chromosome-borne prophages that constitute 2 to 5% of the chromosome. In contrast, other genes missing in the shotgun genomes were randomly scattered, constituting ~0.5% of the chromosome. Among different outbreak strains of the same CC, most gene-scale differences were due to gain/loss of mobile genetic elements, such as plasmids, chromosome-borne prophages, a Tn916 like transposon, and Listeria Genomic Island 2. The nucleotide variations in the same prophage and the same plasmid shared among isolates of the same outbreak strain were limited, which enabled different WGS tools to unambiguously cluster isolates of the same outbreak strain. In some outbreak strains, correlation between prophage gain/loss and single nucleotide polymorphism (SNP) accumulations in the genome backbone were observed.

Virulence assessment of Listeria monocytogenes grown in different foods using a Galleria mellonella model

Various produce including cantaloupe, caramel-coated apples, and packaged salads, have been recognized in recent years as vehicles for listeriosis, a human foodborne disease caused by intracellular pathogen Listeria monocytogenes. Our knowledge regarding the role of these foods in L. monocytogenes virulence, however, is limited. Understanding their role in modulating L. monocytogenes virulence can be useful in risk assessments and for developing control measures. In this study, we employed the Galleria mellonella larvae model to evaluate virulence potential of fifteen clinical, environmental and food isolates of L. monocytogenes, related to three major outbreaks, after growth on different foods. The non-human pathogen Listeria innocua was also included in the panel. Strains were inoculated in parallel in 5ml of brain heart infusion (BHI) broth, and on the surfaces of cantaloupe and apple fragments (5g each) at about 10⁵ colony forming units (CFU)/ml/fragment. One set of inoculated broth and food fragments was incubated at 10°C for 5 days while the second set was kept at 25°C for 3 days. L. monocytogenes cells were recovered from the fruits and BHI, washed twice, re-suspended in saline, and used to inoculate G. mellonella larvae at final concentrations of 10⁶ and 10⁵ CFU/larva. The larvae were incubated at 37°C and monitored for mortality (LT₅₀—time taken to kill 50% of the larvae) and phenotypic changes over seven days. L. monocytogenes grown on cantaloupe and apple flesh surfaces resulted in higher virulence than when grown in BHI. L. monocytogenes infection at 10⁶ CFU/larvae resulted in an average LT₅₀ of ≤ 30, 36 and 47 hours on cantaloupe, apples and BHI, respectively. These results represent a 2.5–4-fold increased mortality compared with an LT₅₀ ≥120 hours in larvae infected with the same doses of L. innocua grown in corresponding matrices. Similar trends were also recorded with doses of about 10⁵ CFU /larvae.

Genomic and phenomic analysis of a marine bacterium, Photobacterium marinum J15

Photobacterium species are widely distributed in the marine environment. The overall metabolism of this genus remains largely unknown. In order to improve our knowledge on this bacterium, the relationship between the genome and phenome of the Photobacterium isolate was analyzed. The cream colored, Gram-negative, rod-shaped and motile bacterial strain, J15, was isolated from marine water of Tanjung Pelepas, Johor, Malaysia. The 5,684,538 bp genome of strain J15 comprised 3 contigs (2 chromosomes and 1 plasmid) with G + C content of 46.39 % and contained 4924 protein-coding genes including 180 tRNAs and 40 rRNAs. The phenotypic microarray (PM) as analyzed using BIOLOG showed the utilization of; i) 93 of the 190 carbon sources tested, where 61 compounds were used efficiently; ii) 41 of the 95 nitrogen sources tested, where 22 compounds were used efficiently; and iii) 3 of the 94 phosphorous and sulphur sources tested. Furthermore, high tolerance to osmotic stress, basic pH and toxic compounds as well as resistance to antibiotics of strain J15 were determined by BIOLOG PM. The ANI and kSNP analyses revealed that strain J15 to be the same species with Photobacterium marinum AK15 with ANI value of 96.93 % and bootstrapping value of 100 in kSNP. Based on the ANI and kSNP analyses, strain J15 was identified as P. marinum J15.

Comparative Genomics of Streptococcus thermophilus Support Important Traits Concerning the Evolution, Biology and Technological Properties of the Species

Streptococcus thermophilus is a major starter for the dairy industry with great economic importance. In this study we analyzed 23 fully sequenced genomes of S. thermophilus to highlight novel aspects of the evolution, biology and technological properties of this species. Pan/core genome analysis revealed that the species has an important number of conserved genes and that the pan genome is probably going to be closed soon. According to whole genome phylogeny and average nucleotide identity (ANI) analysis, most S. thermophilus strains were grouped in two major clusters (i.e., clusters A and B). More specifically, cluster A includes strains with chromosomes above 1.83 Mbp, while cluster B includes chromosomes below this threshold. This observation suggests that strains belonging to the two clusters may be differentiated by gene gain or gene loss events. Furthermore, certain strains of cluster A could be further subdivided in subgroups, i.e., subgroup I (ASCC 1275, DGCC 7710, KLDS SM, MN-BM-A02, and ND07), II (MN-BM-A01 and MN-ZLW-002), III (LMD-9 and SMQ-301), and IV (APC151 and ND03). In cluster B certain strains formed one distinct subgroup, i.e., subgroup I (CNRZ1066, CS8, EPS, and S9). Clusters and subgroups observed for S. thermophilus indicate the existence of lineages within the species, an observation which was further supported to a variable degree by the distribution and/or the architecture of several genomic traits. These would include exopolysaccharide (EPS) gene clusters, Clustered Regularly Interspaced Short Palindromic Repeats (CRISPRs)-CRISPR associated (Cas) systems, as well as restriction-modification (R-M) systems and genomic islands (GIs). Of note, the histidine biosynthetic cluster was found present in all cluster A strains (plus strain NCTC12958^T) but was absent from all strains in cluster B. Other loci related to lactose/galactose catabolism and urea metabolism, aminopeptidases, the majority of amino acid and peptide transporters, as well as amino acid biosynthetic pathways were found to be conserved in all strains suggesting their central role for the species. Our study highlights the necessity of sequencing and analyzing more S. thermophilus complete genomes to further elucidate important aspects of strain diversity within this starter culture that may be related to its application in the dairy industry.

Genomic Diversity of Common Sequence Types of Listeria monocytogenes Isolated from Ready-to-Eat Products of Animal Origin in South Africa

Listeria monocytogenes is a highly fatal foodborne causative agent that has been implicated in numerous outbreaks and related deaths of listeriosis in the world. In this study, six L. monocytogenes isolated from ready-to-eat (RTE) meat products were analysed using Whole Genome Sequencing (WGS) to identify virulence and resistance genes, prophage sequences, PCR-serogroups, and sequence types (STs). The WGS identified four different STs (ST1, ST121, ST204, and ST876) that belonged to serogroup 4b (lineage I) and 1/2a (lineage II). Core genome, and average nucleotide identity (ANI) phylogenetic analyses showed that the majority of strains from serogroup 4b (lineage I) clustered together. However, two isolates that belong to serogroup 1/2a (lineage II) grouped far from each other and the other strains. Examination of reference-guided scaffolds for the presence of prophages using the PHAge Search Tool Enhanced Release (PHASTER) software identified 24 diverse prophages, which were either intact or incomplete/questionable. The National Center for Biotechnology Information- Nucleotide Basic Local Alignment Search Tool (NCBI-BLASTn) revealed that Listeria monocytogenes strains in this study shared some known major virulence genes that are encoded in Listeria pathogenicity islands 1 and 3. In general, the resistance profiles for all the isolates were similar and encoded for multidrug, heavy metal, antibiotic, and sanitizer resistance genes. All the isolates in this study possessed genes that code for resistance to common food processing antiseptics such as Benzalkonium chloride.

Societal Implications of the Internet of Pathogens

The growth of pathogen genomics shows no signs of abating. Whole-genome sequencing of clinical viral and bacterial isolates continues to grow in nearly exponential bounds. Reductions in cost driven by new technology have created a seamless environment for generating, sharing, and analyzing pathogen genomes. The high-resolution view of infectious disease transmission dynamics offered by analyzing whole genomes from pathogens, coupled with the genomicist ethic of widespread data sharing, has created a veritable Internet of pathogens, which inadvertently produces new threats to patient privacy and protected heath information. The health care system, and society more generally, have yet to explore the far-reaching privacy concerns raised by readily accessible pathogen genomic data. The recent use of human genomic databases, the existence of freely available alternative data and metadata sources, and lax regulation of collecting publicly available genomes to identify individuals in a criminal context raise concerning parallels about what is possible with pathogen genomics. The growing ability to ascertain culpability for infectious disease transmission at a nearly individual level could change our perspective on disease outbreaks from one based on public health to one based on individual liability. These technological breakthroughs in the absence of an understanding of potential privacy and liability issues lead to questions about the dominant paradigm of better living through pathogen genomics.

Whole genome sequencing uses for foodborne contamination and compliance: Discovery of an emerging contamination event in an ice cream facility using whole genome sequencing

We review how FDA surveillance identifies several ways that whole genome sequencing (WGS) improves actionable outcomes for public health and compliance in a case involving Listeria monocytogenes contamination in an ice cream facility. In late August 2017 FDA conducted environmental sampling inside an ice cream facility. These isolates were sequenced and deposited into the GenomeTrakr databases. In September 2018 the Centers for Disease Control and Prevention contacted the Florida Department of Health after finding that the pathogen analyses of three clinical cases of listeriosis (two in 2013, one in 2018) were highly related to the aforementioned L. monocytogenes isolates collected from the ice cream facility. in 2017. FDA returned to the ice cream facility in late September 2018 and conducted further environmental sampling and again recovered L. monocytogenes from environmental subsamples that were genetically related to the clinical cases. A voluntary recall was issued to include all ice cream manufactured from August 2017 to October 2018. Subsequently, FDA suspended this food facility's registration. WGS results for L. monocytogenes found in the facility and from clinical samples clustered together by 0-31 single nucleotide polymorphisms (SNPs). The FDA worked together with the Centers for Disease Control and Prevention, as well as the Florida Department of Health, and the Florida Department of Agriculture and Consumer Services to recall all ice cream products produced by this facility. Our data suggests that when available isolates from food facility inspections are subject to whole genome sequencing and the subsequent sequence data point to linkages between these strains and recent clinical isolates (i.e., <20 nucleotide differences), compliance officials should take regulatory actions early to prevent further potential illness. The utility of WGS for applications related to enforcement of FDA compliance programs in the context of foodborne pathogens is reviewed.

Identification of Salmonella Taxon-Specific Peptide Markers to the Serovar Level by Mass Spectrometry

We present an LC-MS/MS pipeline to identify taxon-specific tryptic peptide markers for the identification of Salmonella at the genus, species, subspecies, and serovar levels of specificity. Salmonella enterica subsp. enterica serovars Typhimurium and its four closest relatives, Saintpaul, Heidelberg, Paratyphi B, and Muenchen, were evaluated. A decision-tree approach was used to identify peptides common to the five Salmonella proteomes for evaluation as genus-, species-, and subspecies-specific markers. Peptides identified for two or fewer Salmonella strains were evaluated as potential serovar markers. Currently, there are approximately 140 000 assembled bacterial genomes publicly available, more than 8500 of which are for Salmonella. Consequently, the specificity of each candidate peptide marker was confirmed across all publicly available protein sequences in the NCBI nonredundant (nr) database. The performance of a subset of candidate taxon-specific peptide markers was evaluated in a targeted mass-spectrometry method. The presented workflow offers a marked improvement in specificity over existing MALDI-TOF-based bacterial identification platforms for the identification of closely related Salmonella serovars.

New Aspects on Listeria monocytogenes ST5-ECVI Predominance in a Heavily Contaminated Cheese Processing Environment

The eradication of Listeria monocytogenes from food chains is still a great challenge for the food industry and control authorities since some clonal complexes (CCs) are either better adapted to food processing environments (FPEs) or are globally widespread. In this work, we focus on the in-house evolution of L. monocytogenes genotypes collected from a heavily contaminated FPE whose contamination pattern underwent a massive and yet unexplained change. At the beginning of the sampling in 2010, a high variety of most likely transient L. monocytogenes genotypes was detected belonging to sequence type (ST) 1, ST7, ST21, ST37. After several efforts to intensify the hygiene measures, the variability was reduced to L. monocytogenes ST5 that was dominant in the following years 2011 and 2012. We aimed to elucidate possible genetic mechanisms responsible for the high abundance and persistence of ST5 strains in this FPE. Therefore, we compared the genomes of six L. monocytogenes ST5 strains to the less frequently occurring transient L. monocytogenes ST37 and ST204 from the same FPE as well as the highly abundant ST1 and ST21 isolated in 2010. Whole genome analysis indicated a high degree of conservation among ST5 strains [average nucleotide identity (ANI) 99.93-99.99%; tetranucleotide correlation 0.99998-0.99999]. Slight differences in pulsed field gel electrophoresis (PFGE) patterns of two ST5 isolates could be explained by genetic changes in the tRNA-Arg-TCT prophages. ST5 and ST204 strains harbored virtually identical 91 kbp plasmids related to plasmid group 2 (pLM80 and pLMUCDL175). Interestingly, highly abundant genotypes present in the FPE in 2010 did not harbor any plasmids. The ST5 plasmids harbored an efflux pump system (bcrABC cassette) and heavy metal resistance genes possibly providing a higher tolerance to disinfectants. The pLM80 prototype plasmids most likely provide important genetic determinants for a better survival of L. monocytogenes in the FPE. We reveal short-term evolution of L. monocytogenes strains within the same FPE over a 3 year period and our results suggest that plasmids are important for the persistence of ST5 strains in this FPE.

Geographically widespread honeybee-gut symbiont subgroups show locally distinct antibiotic-resistant patterns

How long-term antibiotic treatment affects host bacterial associations is still largely unknown. The honeybee-gut microbiota has a simple composition, so we used this gut community to investigate how long-term antibiotic treatment affects host-associated microbiota. We investigated the phylogenetic relatedness, genomic content (GC percentage, genome size, number of genes and CRISPR) and antibiotic-resistant genes (ARG) for strains from two abundant members of the honeybee core gut microbiota (Gilliamella apicola and Snodgrassella alvi). Domesticated honeybees are subjected to geographically different management policies, so we used two research apiaries, representing different antibiotic treatment regimens in their apiculture: low antibiotic usage (Norway) and high antibiotic usage (Arizona, USA). We applied whole-genome shotgun sequencing on 48 G. apicola and 22 S. alvi. We identified three predominating subgroups of G. apicola in honeybees from both Norway and Arizona. For G. apicola, genetic content substantially varied between subgroups and distance similarity calculations showed similarity discrepancy between subgroups. Functional differences between subgroups, such as pectin-degrading enzymes (G. apicola), were also identified. In addition, we identified horizontal gene transfer (HGT) of transposon (Tn10)-associated tetracycline resistance (Tet B) across the G. apicola subgroups in the Arizonan honeybees, using interspace polymorphisms in the Tet B determinant. Our results support that honeybee-gut symbiont subgroups can resist long-term antibiotic treatment and maintain functionality through acquisition of geographically distinct antibiotic-resistant genes by HGT.

A Comprehensive Evaluation of the Genetic Relatedness of Listeria monocytogenes Serotype 4b Variant Strains

Recently, we have identified a link between four listeriosis incidents/outbreaks to a variant of Listeria monocytogenes (Lm) serotype 4b strains, 4bV. Although 4bV strains have been reported from clinical specimens as well as from foods, listeriosis outbreaks occurring in 2014–2016 were the first reported outbreaks involving 4bV in the USA. Since traditional typing methods do not detect members of this group, we undertook a systematic and retrospective analysis of all Lm in the NCBI WGS Sequence Read Archive database to investigate the burden of 4bV strains among all listeriosis cases. This analysis identified the presence of isolates causing sporadic cases as well as those associated with the aforementioned outbreaks, as determined by WGS and traditional epidemiology. In total, approximately 350 Lm 4bV strains were identified from multiple parts of the USA as well as from Australia and Chile, dating back to 2001. The genomic relatedness of these strains was compared using the CFSAN SNP Pipeline and multi-virulence-locus sequence typing (MVLST). Using the CFSAN Pipeline tool, the 4bV strains were found to group into seven clusters that were separate from 4b strains. All seven clades appeared to contain isolates from both clinical and non-clinical sources. Conversely, the MVLST analysis revealed that practically all of the strains belonged to a single clade, suggesting that 4bV strains from disparate geographic regions and sources are under varied selective pressure, restricting diversity across these six virulence loci while allowing more variability across the genome as a whole. Further evaluation of these 4bV strains identified genes potentially acquired from a lineage II source external to the lmo0733–lmo0739 region, as well as highly conserved SNPs unique to the 4bV strains when compared to those from other lineages. Taken together, these data suggest that 4bV strains have undergone adaptive responses to selective pressures that may enhance survival in the environment while maintaining the pathogenic potential of serotype 4b strains.

Whole-Genome Sequencing and Concordance Between Antimicrobial Susceptibility Genotypes and Phenotypes of Bacterial Isolates Associated with Bovine Respiratory Disease

Extended laboratory culture and antimicrobial susceptibility testing timelines hinder rapid species identification and susceptibility profiling of bacterial pathogens associated with bovine respiratory disease, the most prevalent cause of cattle mortality in the United States. Whole-genome sequencing offers a culture-independent alternative to current bacterial identification methods, but requires a library of bacterial reference genomes for comparison. To contribute new bacterial genome assemblies and evaluate genetic diversity and variation in antimicrobial resistance genotypes, whole-genome sequencing was performed on bovine respiratory disease–associated bacterial isolates (Histophilus somni, Mycoplasma bovis, Mannheimia haemolytica, and Pasteurella multocida) from dairy and beef cattle. One hundred genomically distinct assemblies were added to the NCBI database, doubling the available genomic sequences for these four species. Computer-based methods identified 11 predicted antimicrobial resistance genes in three species, with none being detected in M. bovis. While computer-based analysis can identify antibiotic resistance genes within whole-genome sequences (genotype), it may not predict the actual antimicrobial resistance observed in a living organism (phenotype). Antimicrobial susceptibility testing on 64 H. somni, M. haemolytica, and P. multocida isolates had an overall concordance rate between genotype and phenotypic resistance to the associated class of antimicrobials of 72.7% (P < 0.001), showing substantial discordance. Concordance rates varied greatly among different antimicrobial, antibiotic resistance gene, and bacterial species combinations. This suggests that antimicrobial susceptibility phenotypes are needed to complement genomically predicted antibiotic resistance gene genotypes to better understand how the presence of antibiotic resistance genes within a given bacterial species could potentially impact optimal bovine respiratory disease treatment and morbidity/mortality outcomes.

Serotype to genotype: The changing landscape of listeriosis outbreak investigations

The classical definition of a disease outbreak is the occurrence of cases of disease in excess of what would normally be expected in a community, geographical area or time period. The establishment of an outbreak then starts with the identification of an incidence of cases above the normally expected threshold during a given time period. Subsequently, the cases are examined using a variety of subtyping methods to identify potential linkages. As listeriosis disease has a long incubation period, relating a single source or multiple sources of contaminated food to clinical disease is challenging and time consuming. The vast majority of human listeriosis cases are caused by three serotypes, 1/2a, 1/2b, and 4b. Thus serotyping of isolates from suspected foods and clinical samples, although useful for eliminating some food sources, has a very limited discriminatory power. The advent of faster and more affordable sequencing technology, coupled with increased computational power, has permitted comparisons of whole Listeria genome sequences from isolates recovered from clinical, food, and environmental sources. These analyses made it possible to identify outbreaks and the source much more accurately and faster, thus leading to a reduction in number of illnesses as well as a reduction in economic losses. Initial DNA sequence information also facilitated the development of a simple molecular serotype protocol which allowed for the identification of major disease causing serotypes of L. monocytogenes, including a clade of 4b variant (4bV) strains of L. monocytogenes involved in at least 3 more recent listeriosis outbreaks in the US. Furthermore, data generated using whole genome sequence (WGS) analyses was successfully utilized to develop a pan-genomic DNA microarray as well as a single nucleotide polymorphism (SNP) based analysis. Herein, we present and compare, the two recently developed sub-typing technologies and discuss how these methods are not only important in outbreak investigations, but could also shed light on possible adaptations to different foods and environments.

A clade of Listeria monocytogenes serotype 4b variant strains linked to recent listeriosis outbreaks associated with produce from a defined geographic region in the US

Four listeriosis incidences/outbreaks, spanning 19 months, have been linked to Listeria monocytogenes serotype 4b variant (4bV) strains. Three of these incidents can be linked to a defined geographical region, while the fourth is likely to be linked. In this study, whole genome sequencing (WGS) of strains from these incidents was used for genomic comparisons using two approached. The first was JSpecies tetramer, which analyzed tetranucleotide frequency to assess relatedness. The second, the CFSAN SNP Pipeline, was used to perform WGS SNP analyses against three different reference genomes to evaluate relatedness by SNP distances. In each case, unrelated strains were included as controls. The analyses showed that strains from these incidents form a highly related clade with SNP differences of ≤101 within the clade and >9000 against other strains. Multi-Virulence-Locus Sequence Typing, a third standardized approach for evaluation relatedness, was used to assess the genetic drift in six conserved, known virulence loci and showed a different clustering pattern indicating possible differences in selection pressure experienced by these genes. These data suggest a high degree of relatedness among these 4bV strains linked to a defined geographic region and also highlight the possibility of alterations related to adaptation and virulence.

Genomic characterization of eight Ensifer strains isolated from pristine caves and a whole genome phylogeny of Ensifer (Sinorhizobium)

A total of eight Ensifer sp. strains were isolated from two pristine cave environments. One strain was isolated from a cave water pool located in the Wind Cave National Park, South Dakota, USA and the remaining seven strains were isolated from Lechuguilla Cave of Carlsbad Caverns National Park, New Mexico, USA. Whole genome sequencing and comparative genomic analyses of the eight isolates compared to various type strains from the genera Ensifer and Sinorhizobium demonstrates that although members in these genera can be phylogenetically separated into two distinct clades, the percentage of conserved proteins (POCP) between various type strains from Ensifer and Sinorhizobium are consistently higher than 50%, providing strong genomic evidence to support the classification of the genera Ensifer and Sinorhizobium into a single genus.

Core Genome Multilocus Sequence Typing for Identification of Globally Distributed Clonal Groups and Differentiation of Outbreak Strains of Listeria monocytogenes

Many listeriosis outbreaks are caused by a few globally distributed clonal groups, designated clonal complexes or epidemic clones, of Listeria monocytogenes, several of which have been defined by classic multilocus sequence typing (MLST) schemes targeting 6 to 8 housekeeping or virulence genes. We have developed and evaluated core genome MLST (cgMLST) schemes and applied them to isolates from multiple clonal groups, including those associated with 39 listeriosis outbreaks. The cgMLST clusters were congruent with MLST-defined clonal groups, which had various degrees of diversity at the whole-genome level. Notably, cgMLST could distinguish among outbreak strains and epidemiologically unrelated strains of the same clonal group, which could not be achieved using classic MLST schemes. The precise selection of cgMLST gene targets may not be critical for the general identification of clonal groups and outbreak strains. cgMLST analyses further identified outbreak strains, including those associated with recent outbreaks linked to contaminated French-style cheese, Hispanic-style cheese, stone fruit, caramel apple, ice cream, and packaged leafy green salad, as belonging to major clonal groups. We further developed lineage-specific cgMLST schemes, which can include accessory genes when core genomes do not possess sufficient diversity, and this provided additional resolution over species-specific cgMLST. Analyses of isolates from different common-source listeriosis outbreaks revealed various degrees of diversity, indicating that the numbers of allelic differences should always be combined with cgMLST clustering and epidemiological evidence to define a listeriosis outbreak.

IMPORTANCE Classic multilocus sequence typing (MLST) schemes targeting internal fragments of 6 to 8 genes that define clonal complexes or epidemic clones have been widely employed to study L. monocytogenes biodiversity and its relation to pathogenicity potential and epidemiology. We demonstrated that core genome MLST schemes can be used for the simultaneous identification of clonal groups and the differentiation of individual outbreak strains and epidemiologically unrelated strains of the same clonal group. We further developed lineage-specific cgMLST schemes that targeted more genomic regions than the species-specific cgMLST schemes. Our data revealed the genome-level diversity of clonal groups defined by classic MLST schemes. Our identification of U.S. and international outbreaks caused by major clonal groups can contribute to further understanding of the global epidemiology of L. monocytogenes.