Pan-genome analysis of the emerging foodborne pathogen Cronobacter spp. suggests a species-level bidirectional divergence driven by niche adaptation

Genomic Analysis of Cronobacter condimenti s37: Identification of Resistance and Virulence Genes and Comparison with Other Cronobacter and Closely Related Species

Cronobacter condimenti are environmental commensals that have not been associated with any clinical infections. To date, they are the least understood and described Cronobacter species within the genus. The objective of this study was to use a draft genome sequence (DGS) of the Cronobacter condimenti strain s37 to screen for genes encoding for antibiotic resistance, virulence, response to environmental stress, and biofilm formation. The strain was isolated in Poland from commercial small radish sprouts. This is the second genome of this species available in the GenBank database. The comparative genome analysis (cgMLST) of C. condimenti s37 with other Cronobacter spp. including the pathogenic species C. sakazakii and the plant-associated closely related genera Franconibacter and Siccibacter was also performed. The assembled and annotated genome of the C. condimenti s37 genome was 4,590,991 bp in length, with a total gene number of 4384, and a GC content of 55.7%. The s 37 genome encoded for genes associated with resistance to stressful environmental conditions (metal resistance genes: zinc, copper, osmotic regulation, and desiccation stress), 17 antimicrobial resistance genes encoding resistance to various classes of antibiotics and 50 genes encoding for the virulence factors. The latter were mainly genes associated with adhesion, chemotaxis, hemolysis, and biofilm formation. Cg-MLST analysis (3991 genes) revealed a greater similarity of C. condimenti s37 to S. turicensis, F. pulveris, and C. dublinensis than to other species of the genus Cronobacter. Studies on the diversity, pathogenicity, and virulence of Cronobacter species isolated from different sources are still insufficient and should certainly be continued. Especially the analysis of rare strains such as s37 is very important because it provides new information on the evolution of these bacteria. Comparative cgMLST analysis of s37 with other Cronobacter species, as well as closely related genera Franconibacter and Siccibacter, complements the knowledge on their adaptability to specific environments such as desiccation.

Cronobacter spp. in foods of plant origin: occurrence, contamination routes, and pathogenic potential

Cronobacter is an emerging bacterial pathogen associated with infections such as necrotizing enterocolitis, sepsis, and meningitis in neonates and infants, related to the consumption of powdered infant formula. In addition, this bacterium can also cause infections in adults by the ingestion of other foods. Thus, this review article aims to report the occurrence and prevalence of Cronobacter spp. in foods of plant origin, as well as the possible sources and routes of contamination in these products, and the presence of pathogenic strains in these foods. Cronobacter was present in a wide variety of cereal-based foods, vegetables, herbs, spices, ready-to-eat foods, and foods from other categories. This pathogen was also found in cultivation environments, such as soils, compost, animal feces, rice and vegetable crops, as well as food processing industries, and domestic environments, thus demonstrating possible contamination routes. Furthermore, sequence types (ST) involved in clinical cases and isolates resistant to antibiotics were found in Cronobacter strains isolated from food of plant origin. The identification of Cronobacter spp. in plant-based foods is of great importance to better elucidate the vehicles and routes of contamination in the primary production chain and processing facility, until the final consumption of the food, to prevent infections.

Phylogenomics of the Liquorilactobacillus Genus

The genus Liquorilactobacillus is a new genus commonly found in wine and plants. Despite its significance, previous studies on Liquorilactobacillus are primarily focused on phenotypic experiments, with limited genome-level studies. This study used comparative genomics to analyze 24 genomes from the genus Liquorilactobacillus, including two novel sequenced strains (IMAU80559 and IMAU80777). A phylogenetic tree of 24 strains was constructed based on 122 core genes and divided into two clades, A and B. Significant differences in GC content were observed between the two clades (P = 10e^-4). Additionally, change revealed to suggests that clade B has more exposure to prophage infection having an upgraded immune system. Further analysis of functional annotation and selective pressure suggests that clade A was subjected to greater selection pressure than B clade (P = 3.9e^-6) and had higher number of functional types annotated than clade B (P = 2.7e^-3), while clade B had a lower number of pseudogenes than clade A (P = 1.9e^-2). The findings suggest that differently prophages and environmental stress may have influenced the common ancestor of clades A and B during evolution, leading to the development of two distinct clades.

Comparative phylogenomics of Trueperella pyogenes reveals host-based distinction of strains

Trueperella pyogenes, an opportunistic pathogen causes various ailments in different animals. Different strains from different animals have distinct characters phenotypically and genotypically. Hence understanding the strains in a particular geographical location helps in framing the preventive measures. Comparative genomics of all the available T. pyogenes genome in the NCBI was conducted to understand the relatedness among strains. Whole genome phylogeny showed host associated clustering of strains recovered from swine lungs. Core genome phylogeny also showed host associated clustering mimicking whole genome phylogeny results. MLST analysis showed that there was higher diversity among cattle strains. Multidimensional scaling revealed five swine clusters, two cattle and buffalo clusters. Pangenome analysis also showed that T. pyogenes had an open genome with 57.09% accessory genome. Host specific genes were identified by pangenome analysis, and (R)-citramalate synthase was specific for swine strains of Asian origin. Host specifc genes identified by pangenome analysis can be exploited for developing a molecular assay to specifically identify the strains. The study shows that MLST having higher discriminatory power can be used as an epidemiological tool for strain discrimination of T. pyogenes.

Halocatena salina sp. nov., a filamentous halophilic archaeon isolated from Aiding Salt Lake

A filamentous cell-shaped halophilic archaeon (strain AD-1^T) was isolated from Aiding Salt Lake, PR China. Its colonies on HCM7 agar plates were pinkish white, 1-4 mm (diameter), elevated and round. The optimum conditions for growth were observed at 42 °C, 4.3 M NaCl, 0.01 M MgCl₂ and pH 7. Strain AD-1^T could hydrolyse Tween 60, Tween 80, starch and gelatin. Phylogenetic analysis based on 16S rRNA gene, rpoB' and the concatenated 484 single-copy orthologous proteins revealed that strain AD-1^T formed a clade with Halocatena pleomorpha SPP-AMP-1^T. The average nucleotide identity and in silico DNA-DNA hybridization values between strain AD-1^T and Halocatena pleomorpha SPP-AMP-1^T were both below the species delineation thresholds (95~96 and 70 %, respectively). The major phospholipids of strain AD-1^T were phosphatidic acid, phosphatidylglycerol and phosphatidylglycerol phosphate methyl ester, while the major glycolipids were sulphated galactosyl mannosyl glucosyl diether, galactosyl mannosyl glucosyl diether and glucosyl mannosyl glucosyl diether. The phenotypic, phylogenetic and genome-based analyses suggested that strain AD-1^T (=CGMCC 1.13724^T=JCM 32960^T) represents a novel species, for which the name Halocatena salina sp. nov. is proposed.

Genome insights from the identification of a novel Pandoraea sputorum isolate and its characteristics

In this study, we sequenced a bacteria isolate Pandoraea sp. 892iso isolated from a Phytophthora rubi strain which is an important plant pathogenic oomycete, identified through genome and combined the data with existing genomic data from other 28 the genus of Pandoraea species. Next, we conducted a comparative genomic analysis of the genome structure, evolutionary relationships, and pathogenic characteristics of Pandoraea species. Our results identified Pandoraea sp. 892iso as Pandoraea sputorum at both the genome and gene levels. At the genome level, we carried out phylogenetic analysis of single-copy, gene co-linearity, ANI (average nucleotide identity) and AAI (average amino acid identity) indices, rpoB similarity, MLSA phylogenetic analysis, and genome-to-genome distance calculator calculations to identify the relationship between Pandoraea sp. 892iso and P. sputorum. At the gene level, the quorum sensing genes ppnI and ppnR and the OXA-159 gene were assessed. It is speculated that Pandoraea sp. 892iso is the endosymbiont of the Oomycetes strain of Phytophthora rubi.

Virulence, antimicrobial susceptibility, and phylogenetic analysis of Cronobacter sakazakii isolates of food origins from Jordan

AIMS

to characterize a collection of Cronobacter sakazakii isolates collected from various origins in Jordan.

METHODS AND RESULTS

the isolates were characterized using 16S rRNA sequencing, DNA microarray, multi-locus sequence typing (MLST), O-serotyping, virulence gene identification, and antibiotic susceptibility testing. The identities and phylogenetic relatedness revealed that C. sakazakii sequence type 4 (ST4) and Csak O:1 serotype was the most prevalent STs and serovars among these C. sakazakii strains. PCR screening of putative virulence genes showed that the siderophore-interacting protein gene (sip) and iron acquisition gene clusters (eitCBAD and iucABCD/iutA) were the most detected genes with noticeable variability in the type 6 secretion system (T6SS) and filamentous hemagglutinin/adhesion (FHA) gene loci. The antibiotic resistance profiles revealed that the majority of the isolates were susceptible to all antibiotics used despite harboring a class C β-lactamase resistance gene.

CONCLUSIONS

the results described in this report provide additional insights about the considerable genotypic and phenotypic heterogeneity within C. sakazakii.

SIGNIFICANCE AND IMPACT OF THE STUDY

the information reported in this study might be of great value in understanding the origins of C. sakazakii isolates, in addition to their diversity and variability, which might be helpful in preventing future outbreaks of this pathogen.

Characterization of Cronobacter sakazakii Strains Originating from Plant-Origin Foods Using Comparative Genomic Analyses and Zebrafish Infectivity Studies

Cronobacter sakazakii continues to be isolated from ready-to-eat fresh and frozen produce, flours, dairy powders, cereals, nuts, and spices, in addition to the conventional sources of powdered infant formulae (PIF) and PIF production environments. To understand the sequence diversity, phylogenetic relationship, and virulence of C. sakazakii originating from plant-origin foods, comparative molecular and genomic analyses, and zebrafish infection (ZI) studies were applied to 88 strains. Whole genome sequences of the strains were generated for detailed bioinformatic analysis. PCR analysis showed that all strains possessed a pESA3-like virulence plasmid similar to reference C. sakazakii clinical strain BAA-894. Core genome analysis confirmed a shared genomic backbone with other C. sakazakii strains from food, clinical and environmental strains. Emerging nucleotide diversity in these plant-origin strains was highlighted using single nucleotide polymorphic alleles in 2000 core genes. DNA hybridization analyses using a pan-genomic microarray showed that these strains clustered according to sequence types (STs) identified by multi-locus sequence typing (MLST). PHASTER analysis identified 185 intact prophage gene clusters encompassing 22 different prophages, including three intact Cronobacter prophages: ENT47670, ENT39118, and phiES15. AMRFinderPlus analysis identified the CSA family class C β-lactamase gene in all strains and a plasmid-borne mcr-9.1 gene was identified in three strains. ZI studies showed that some plant-origin C. sakazakii display virulence comparable to clinical strains. Finding virulent plant-origin C. sakazakii possessing significant genomic features of clinically relevant STs suggests that these foods can serve as potential transmission vehicles and supports widening the scope of continued surveillance for this important foodborne pathogen.

Comparative Genomic Analysis Reveals Intestinal Habitat Adaptation of Ligilactobacillus equi Rich in Prophage and Degrading Cellulase

Ligilactobacillus equi is common in the horse intestine, alleviates the infection of Salmonella, and regulates intestinal flora. Despite this, there have been no genomic studies on this species. Here, we provide the genomic basis for adaptation to the intestinal habitat of this species. We sequenced the genome of L. equi IMAU81196, compared this with published genome information from three strains in NCBI, and analyzed genome characteristics, phylogenetic relationships, and functional genes. The mean genome size of L. equi strains was 2.08 ± 0.09 Mbp, and the mean GC content was 39.17% ± 0.19%. The genome size of L. equi IMAU81196 was 1.95 Mbp, and the GC content was 39.48%. The phylogenetic tree for L. equi based on 1454 core genes showed that the independent branch of strain IMAU81196 was far from the other three strains. In terms of genomic characteristics, single-nucleotide polymorphism (SNP) sites, rapid annotation using subsystem technology (RAST), carbohydrate activity enzymes (CAZy), and predictions of prophage, we showed that strain L. equi JCM 10991T and strain DSM 15833T are not equivalent strains.It is worth mentioning thatthestrain of L. equi has numerous enzymes related to cellulose degradation, and each L. equi strain investigated contained at least one protophage. We speculate that this is the reason why these strains are adapted to the intestinal environment of horses. These results provide new research directions for the future.

Invited review: Stress resistance of Cronobacter spp. affecting control of its growth during food production

Members of the Cronobacter genus include food-borne pathogens that can cause infections in infants, with a mortality rate as high as 40 to 80%. The high fatality rate of Cronobacter and its isolation from numerous types of food, especially from powdered infant formula, demonstrate the serious nature of this organism. The source tracking of Cronobacter spp. and the analysis of high-frequency species from different sources are helpful for a more targeted control. Furthermore, the persistence during food processing and storage may be attributed to strong resistance of Cronobacter spp. to environment stresses such as heat, pH, and desiccation. There are many factors that support the survival of Cronobacter spp. in harsh environments, such as some genes, regulatory systems, and biofilms. Advanced detection technology is helpful for the strict monitoring of Cronobacter spp. In addition to the traditional heat treatment, many new control techniques have been developed, and the ability to control Cronobacter spp. has been demonstrated. The control of this bacteria is required not only during manufacture, but also through the selection of packaging methods to reduce postprocessing contamination. At the same time, the effect of inactivation methods on product quality and safety must be considered. This review considers the advances in our understanding of environmental stress response in Cronobacter spp. with special emphasis on its implications in food processing.

Comparative genomics reveals environmental adaptation differences between Cronobacter species

The genus Cronobacter is an opportunistic food-borne pathogen which is able to adapt to diverse environments and shows considerable genetic diversity. Genomic analysis can be used to reveal the variation across the genus with respect to virulence, drug resistance and factors involved in horizontal gene transfer mechanisms, such as integrons, conjugative plasmids, and recombinases. In this study, whole-genome comparative analysis of 27 Cronobacter genomes (12 existing and 15 newly assembled genomes) was performed. A total of 110,010 protein-coding genes were grouped into 11,262 clusters, including 2637 core genes, 4814 strain-specific genes and 3811 dispensable genes. Clusters of Orthologous Groups (COG) analysis indicated that 97.35% of the core genes were for substrate transport and metabolism, and the antibiotic resistance genetic determinants were classified into 136 antibiotic resistance ontologies (AROs). A total of 80 genomic islands (GIs) were identified which contained several type VI secretion system gene clusters, and these were likely to have been acquired by horizontal gene transfer. This study has generated a comprehensive characterization of the genus Cronobacter, leading to a better understanding of the mechanisms and ecological functions among the genome features, speciation, and environmental adaptation strategies.

Novel phage vB_CtuP_B1 for controlling Cronobacter malonaticus and Cronobacter turicensis in ready-to-eat lettuce and powered infant formula

Cronobacter spp. are important foodborne pathogens that are a threat to people of all ages, but especially neonates and infants. Bacteriophages are biological agents that are potentially useful for the control of foodborne pathogens. However, there has been little research on the control of C. malonaticus and C. turicensis using bacteriophages. In the present study, a novel lytic phage vB_CtuP_B1 (hereafter referred to as B1)-which can simultaneously lyse C. malonaticus and C. turicensis- was isolated from river water in Guangzhou, China, and was used in the control of Cronobacter contaminated food. The phage has a short tail, and has been identified as a new species of Kayfunavirus based on genomic and phylogenetic analyses. One-step growth and stability assays revealed that phage B1 has a very short latent period (<5 min) and a large burst size (4006 pfu/cell), and is highly stable between 25 and 60 °C and between pH 5 and 11. Its genome encodes two lytic proteins, but does not contain any genes responsible for antibiotic resistance and virulence factors. In broth, the phage B1 completely inhibited the growth of C. malonaticus cro2475W and C. turicensis cro1541A1-1 for up to 6 h. On lettuce, phage B1 reduced the viable count of C. turicensis cro1541A1-1 to below the detection limit for bacteria on lettuce (<10 cfu/mL) after 6 h at 4 °C and 2 h at 25 °C, and also significantly reduced the viable count of C. malonaticus cro2475W at those temperatures. In powdered infant formula, the viable counts of both the phage-treated bacterial hosts were significantly reduced after 2 h of storage at 4 °C or 37 °C. Furthermore, phage B1 reduced the viable count of C. turicensis cro1541A1-1 to below the detection limit (<10 cfu/mL) from 4 h to 24 h at 37 °C. It significantly inhibited the growth of C. turicensis cro1541A1-1 than that of C. malonaticus cro2475W (P < 0.05). In conclusion, phage B1 with high stability and strong lytic ability is potentially useful for controlling C. malonaticus and C. turicensis.

Diversity of Cronobacter genus isolated between 1970 and 2019 on the American continent and genotyped using multi-locus sequence typing

This study aimed to evaluate the Cronobacter spp. strains isolated on the American continent and characterized using multi-locus sequence typing (MLST) available in the PubMLST database and current literature. From 465 Cronobacter spp. strains, the majority (n = 267, 57.4%) was from North America, mainly from USA (n = 234) and 198 (42.6%) were from South America, mainly from Brazil (n = 196). A total of 232 (49.9%) were isolated from foods, 102 (21.9%) from environmental, 87 (18.7%) from clinical, 27 (5.8%) from PIF, one from water (0.2%) and 16 (3.5%) from unknown sources. A total of five species were represented: Cronobacter sakazakii (374, 80.4%), Cronobacter malonaticus (41, 8.8%), Cronobacter dublinensis (29, 6.2%), Cronobacter turicensis (16, 3.5%) and Cronobacter muytjensii (5, 1.1%). The strains with complete MLST profile (n = 345) were assigned to 98 STs, a ratio of 3.5 strain by ST found and the calculated Simpson`s index was 0.93. The strains showed a high diversity and after eBURST analysis, 30 STs (n = 189) formed 12 single and/or double-locus variant clonal complexes (CC). A total of 38 STs (38.7%) were associated with clinical cases of infection, including well established C. sakazakii CC 1, 4, 8 and 83; C. malonaticus ST60, 307, 394 and 440; and C. sakazakii ST 12 and 494.

Screening of specific nucleic acid targets for Cronobacter sakazakii and visual detection by loop-mediated isothermal amplification and lateral flow dipstick method in powdered infant formula

Due to the lack of specific genes for rapid detection methods of Cronobacter sakazakii in food samples, whole genome sequence analysis was performed in this investigation using the basic local alignment search tool. Forty-two DNA fragments unique to C. sakazakii were mined, then primers were designed and screened by PCR and loop-mediated isothermal amplification (LAMP). Two primer sets, CS1 and CS31, were found as specific and stable primers, with their corresponding nucleic acid targets the CSK29544_00235 gene and CSK29544_03484 gene, respectively. Furthermore, compared with 3 genes reported previously, these 2 genes were verified as more specific to C. sakazakii among Cronobacter species, by sequence similarity alignment using Cronobacter MLST databases (http://pubmlst.org/cronobacter). The specificity of the LAMP reaction approached 100% by using 48 bacterial strains, which included 22 C. sakazakii strains. Subsequently, LAMP was combined with visual lateral flow dipstick (LFD) based on the above 2 nucleic acid targets, and was demonstrated as a rapid, efficient method with high specificity. Finally, the detection sensitivity of this assay system for pure cultures and artificially contaminated milk was measured as 4.5 × 10⁰ cfu/mL and 5.7 × 10¹ cfu/g, respectively. Total time to detection for this assay was within 2 h. Thus, the establishment of this LAMP-LFD method shows great significance and potential for rapid detection of C. sakazakii in powdered infant formula.

Isolation and characterization of new phage vB_CtuP_A24 and application to control Cronobacter spp. in infant milk formula and lettuce

Cronobacter spp. are opportunistic pathogenic bacteria that cause severe diseases in neonates and infants. Bacteriophages are novel antibacterial agents with a potential to control this pathogen. In the current study, a novel lytic Cronobacter phage, vB_CtuP_A24, was isolated from a river in Guangzhou, China. The phage was characterized by a short, non-contractile tail and a long head and identified as a new member of the family Podoviridae. Complete genome sequence analysis of this phage indicated that its genome contained 75,106 bp of DNA, an average GC content of 44.05%, and 108 predicted open reading frames (ORFs). The annotated ORFs were associated with phage structure, packaging, host lysis, DNA manipulation, transcription, and additional functions. Genes responsible for antibiotic resistance, virulence, and toxic effects were not present. Cronobacter phage vB_CtuP_A24 is a novel lytic phage that can lyse five Cronobacter spp. It was stable over a wide range of temperatures (25-60 °C) and pH values (pH 4-11) and displayed a short latent period (approximately 10 min) and a large burst size (500 plaque-forming units (PFUs)/cell). In Luria-Bertani (LB) broth, phage A24 effectively inhibited the growth of C. dublinensis cro280B, C. sakazakii 465G, and C. malonaticus cro695W up to 9, 10, and 12 h, respectively, while in infant milk formulas, it inhibited two strains, C. sakazakii 465G and C. malonaticus cro695W, up to 24 h at 37 °C, with maximum reduction levels of approximately 5.12 ± 3.95 and 7.38 ± 3.03 log10 colony-forming unit (CFU)/mL, respectively. In lettuce, the reduction of all three strains was highly significant at 25 °C. However, the growth of C. dublinensis cro280B and C. sakazakii 465G were not significantly inhibited at 4 °C. In conclusion, Cronobacter spp. phage vB_CtuP_A24, which has lytic ability against five Cronobacter species, stability under different environments, and shows potential as a promising biocontrol agent against Cronobacter spp. in food production.

Pseudonocardia broussonetiae sp. nov., an endophytic actinomycete isolated from the roots of Broussonetia papyrifera

A novel endophytic actinomycete, designated strain Gen 01^T, was isolated from the roots of Broussonetia papyrifera and characterized by using a polyphasic approach. The predominant cellular fatty acids were iso-C_16 : 0, summed feature 3, iso H-C_16 : 1, C_16 : 0 and iso-C_14 : 0. The polar lipid profile consisted of diphosphatidylglycerol, phosphatidylethanolamine, phosphatidylinositol mannosides, phospholipids of unknown structure containing glucosamine inositol, phosphatidylinositol and unidentified phospholipids. The major menaquinone was MK-8 (H₄). The DNA G+C content of the genome sequence, consisting of 7 177 725 bp, was 74.5 mol%. Phylogenetic analysis of the full-length 16S rRNA gene sequences showed that strain Gen 01^T belongs to the genus Pseudonocardia with the highest sequence similarity to Pseudonocardia petroleophila CGMCC 4.1532^T (98.9 %) and lower than 98.7 % similarity to other species of the genus Pseudonocardia with validly published names. The average nucleotide identity and digital DNA-DNAhybridization values between strain Gen 01^T and P. petroleophila CGMCC 4.1532^T were 84.6 and 30.9 %, respectively. Furthermore, the morphological, physiological and biochemical characteristics were sufficient to categorize strain Gen 01^T as being distinct from P. petroleophila CGMCC 4.1532^T. Consequently, based on phenotypic and genotypic characteristics, strain Gen 01^T represents a novel species of the genus Pseudonocardia, for which the name Pseudonocardia broussonetiae sp. nov. is proposed. The type strain is Gen 01^T (=CICC 24820^T=JCM 33840^T).

Steptomyces fagopyri sp. nov., a novel actinomycete isolated from rhizospheric soil of Fagopyrum dibotrys

A novel actinomycete, designated strain QMT-28^T, was isolated from rhizosphere soil of Fagopyrum dibotrys collected from Shuangfeng, Hunan Province, PR China. Strain QMT-28^T grew well on International Streptomyces Project series media and formed well-developed, branched substrate hyphae and aerial mycelium that differentiated into loose spiral spore chains consisting of cylindrical spores with smooth surfaces. The diagnostic diamino acid was ll-diaminopimelic acid and the whole-cell sugars were galactose and glucose. The predominant fatty acids were C_18 : 1 cis9, summed feature 6 (C_18 : 2 cis 9,12/C_18 : 0 a) and C_16 : 0. The polar lipids included diphosphatidylglycerol, hydroxy phospatidylethanolamine, phosphatidylethanolamine, phosphatidylinositol, phosphatidylinositol mannosides, phospholipids of unknown structure containing glucosamine and several unidentified phospholipids. The major menaquinones were MK-9, MK-9(H₂), MK-9(H₄), MK-9(H₆) and MK-9(H₈). The genome size of strain QMT-28^T was about 8.7 Mbp with a G+C content of 71.2 mol%. Phylogenetic analysis showed that the novel strain was closely related to Streptomyces olivochromogenes DSM 40451^T (99.5 % similarity), Streptomyces mirabilis NBRC 13450^T (98.9 %), Streptomyces kanamyceticus NBRC 13414^T (98.9 %), Streptomyces kaempferi I37^T (98.9 %) and Streptomyces arcticus ZLN234^T (98.8 %). However, the average nucleotide identity values, the digital DNA-DNA hybridization values and the multilocus sequence analysis evolutionary distances between this strain and closely related strains showed that it belonged to a distinct species. In addition, these results were also supported by differences in the phenotypic characteristics between QMT-28^T and five closely related type strains. Consequently, strain QMT-28^T should represent a novel species of the genus Streptomyces, with the suggested name Streptomyces fagopyri sp. nov. The type strain is QMT-28^T (=CICC 24808^T=JCM 33796^T).

Streptomyces phaeolivaceus sp. nov. and Streptomyces broussonetiae sp. nov., isolated from the leaves and rhizosphere soil of Broussonetia papyrifera

Two novel actinobacteria, designated strains GY16^T and T44^T, were isolated from the leaves and rhizosphere soil of Broussonetia papyrifera, respectively. A polyphasic approach was used for determining their taxonomic position. Results of 16S rRNA gene sequence analysis indicated that strain GY16^T exhibited highest similarities to Streptomyces cinereoruber subsp. fructofermentans CGMCC 4.1593^T (98.82 %), Streptomyces deccanensis KCTC 19241^T (98.76 %), Streptomyces scabiei NRRL B-16523^T (98.69 %), Streptomyces europaeiscabiei KACC 20186^T (98.69 %) and Streptomyces rishiriensis NBRC 13407^T (98.69 %), and strain T44^T showed 99.2, 99.1, 99.1 and <98.7 % sequence similarities to Streptomyces filipinensis CGMCC 4.1452^T, Streptomyces achromogenes subsp. achromogenes DSM 40028^T, Streptomyces durhamensis DSM 40539^T and other Streptomyces species, respectively. Phylogenetic analysis based on 16S rRNA gene sequences indicated that strain GY16^T formed an independent subclade, which indicated that strain GY16^T should belong to a potential novel species; and strain T44^T was closely related to S. filipinensis CGMCC 4.1452^T, S. achromogenes subsp. achromogenes DSM 40028^T, S. durhamensis DSM 40539^T and S. yokosukanensis DSM 40224^T. However, the multilocus sequence analysis evolutionary distance, average nucleotide identity and DNA-DNA hybridization values between closely related relatives were far from the species-level thresholds. In addition, phenotypic and chemotaxonomic characteristics further confirmed that strains GY16^T and T44^T belonged to two distinct species. Based on these results, it is concluded that the isolated strains represent novel species within the genus Streptomyces, for which the names Streptomyces phaeolivaceus sp. nov. (type strain GY16^T=CICC 24807^T=KCTC 49326^T) and Streptomyces broussonetiae sp. nov. (type strain T44^T=CICC 24819^T=JCM 33918^T) are proposed.

Analysis of the Molecular Diversity Among Cronobacter Species Isolated From Filth Flies Using Targeted PCR, Pan Genomic DNA Microarray, and Whole Genome Sequencing Analyses

Cronobacter species are opportunistic pathogens capable of causing life-threatening infections in humans, with serious complications arising in neonates, infants, immuno-compromised individuals, and elderly adults. The genus is comprised of seven species: Cronobacter sakazakii, Cronobacter malonaticus, Cronobacter turicensis, Cronobacter muytjensii, Cronobacter dublinensis, Cronobacter universalis, and Cronobacter condimenti. Despite a multiplicity of genomic data for the genus, little is known about likely transmission vectors. Using DNA microarray analysis, in parallel with whole genome sequencing, and targeted PCR analyses, the total gene content of two C. malonaticus, three C. turicensis, and 14 C. sakazaki isolated from various filth flies was assessed. Phylogenetic relatedness among these and other strains obtained during surveillance and outbreak investigations were comparatively assessed. Specifically, microarray analysis (MA) demonstrated its utility to cluster strains according to species-specific and sequence type (ST) phylogenetic relatedness, and that the fly strains clustered among strains obtained from clinical, food and environmental sources from United States, Europe, and Southeast Asia. This combinatorial approach was useful in data mining for virulence factor genes, and phage genes and gene clusters. In addition, results of plasmidotyping were in agreement with the species identity for each strain as determined by species-specific PCR assays, MA, and whole genome sequencing. Microarray and BLAST analyses of Cronobacter fly sequence datasets were corroborative and showed that the presence and absence of virulence factors followed species and ST evolutionary lines even though such genes were orthologous. Additionally, zebrafish infectivity studies showed that these pathotypes were as virulent to zebrafish embryos as other clinical strains. In summary, these findings support a striking phylogeny amongst fly, clinical, and surveillance strains isolated during 2010–2015, suggesting that flies are capable vectors for transmission of virulent Cronobacter spp.; they continue to circulate among United States and European populations, environments, and that this “pattern of circulation” has continued over decades.

All ANIs are not created equal: implications for prokaryotic species boundaries and integration of ANIs into polyphasic taxonomy

In prokaryotic taxonomy, a set of criteria is commonly used to delineate species. These criteria are generally based on cohesion at the phylogenetic, phenotypic and genomic levels. One such criterion shown to have promise in the genomic era is average nucleotide identity (ANI), which provides an average measure of similarity across homologous regions shared by a pair of genomes. However, despite the popularity and relative ease of using this metric, ANI has undergone numerous refinements, with variations in genome fragmentation, homologue detection parameters and search algorithms. To test the robustness of a 95-96 % species cut-off range across all the commonly used ANI approaches, seven different methods were used to calculate ANI values for intra- and interspecies datasets representing three classes in the Proteobacteria. As a reference point, these methods were all compared to the widely used blast-based ANI (i.e. ANIb as implemented in JSpecies), and regression analyses were performed to investigate the correlation of these methods to ANIb with more than 130000 individual data points. From these analyses, it was clear that ANI methods did not provide consistent results regarding the conspecificity of isolates. Most of the methods investigated did not correlate perfectly with ANIb, particularly between 90 and 100% identity, which includes the proposed species boundary. There was also a difference in the correlation of methods for the different taxon sets. Our study thus suggests that the specific approach employed needs to be considered when ANI is used to delineate prokaryotic species. We furthermore suggest that one would first need to determine an appropriate cut-off value for a specific taxon set, based on the intraspecific diversity of that group, before conclusions on conspecificity of isolates can be made, and that the resulting species hypotheses be confirmed with analyses based on evolutionary history as part of the polyphasic approach to taxonomy.

The Secretion of Toxins and Other Exoproteins of Cronobacter: Role in Virulence, Adaption, and Persistence

: Cronobacter species are considered an opportunistic group of foodborne pathogenic bacteria capable of causing both intestinal and systemic human disease. This review describes common virulence themes shared among the seven Cronobacter species and describes multiple exoproteins secreted by Cronobacter, many of which are bacterial toxins that may play a role in human disease. The review will particularly concentrate on the virulence factors secreted by C. sakazakii, C. malonaticus, and C. turicensis, which are the primary human pathogens of interest. It has been discovered that various species-specific virulence factors adversely affect a wide range of eukaryotic cell processes including protein synthesis, cell division, and ion secretion. Many of these factors are toxins which have been shown to also modulate the host immune response. These factors are encoded on a variety of mobile genetic elements such as plasmids and transposons; this genomic plasticity implies ongoing re-assortment of virulence factor genes which has complicated our efforts to categorize Cronobacter into sharply defined genomic pathotypes.

Prevalence, Distribution, and Phylogeny of Type Two Toxin-Antitoxin Genes Possessed by Cronobacter Species where C. sakazakii Homologs Follow Sequence Type Lineages

Cronobacter species are a group of foodborne pathogenic bacteria that cause both intestinal and systemic human disease in individuals of all age groups. Little is known about the mechanisms that Cronobacter employ to survive and persist in foods and other environments. Toxin-antitoxin (TA) genes are thought to play a role in bacterial stress physiology, as well as in the stabilization of horizontally-acquired re-combinatorial elements such as plasmids, phage, and transposons. TA systems have been implicated in the formation of a persistence phenotype in some bacterial species including Escherichia coli and Salmonella. This project's goal was to understand the phylogenetic relatedness among TA genes present in Cronobacter. Preliminary studies showed that two typical toxin genes, fic and hipA followed species evolutionary lines. A local database of 22 TA homologs was created for Cronobacter sakazakii and a Python version 3 shell script was generated to extract TA FASTA sequences present in 234 C. sakazakii genomes previously sequenced as part of Center for Food Safety and Applied Nutrition's (CFSAN) GenomeTrakr project. BLAST analysis showed that not every C. sakazakii strain possessed all twenty-two TA loci. Interestingly, some strains contained either a toxin or an antitoxin component, but not both. Five common toxin genes: ESA_00258 (parDE toxin-antitoxin family), ESA_00804 (relBE family), ESA_01887 (relBE family), ESA_03838 (relBE family), and ESA_04273 (YhfG-Fic family) were selected for PCR analysis and the primers were designed to detect these genes. PCR analysis showed that 55 of 63 strains possessed three of these genes Sequence analysis identified homologs of the target genes and some of the strains were PCR-negative for one or more of the genes, pointing to potential nucleotide polymorphisms in those loci or that these toxin genes were absent. Phylogenetic studies using a Cronobacter pan genomic microarray showed that for the most part TAs follow species evolutionary lines except for a few toxin genes possessed by some C. malonaticus and C. universalis strains; this demonstrates that some TA orthologues share a common phylogeny. Within the C. sakazakii strains, the prevalence and distribution of these TA homologs by C. sakazakii strain BAA-894 (a powdered infant formula isolate) followed sequence-type evolutionary lineages. Understanding the phylogeny of TAs among the Cronobacter species is essential to design future studies to realize the physiological mechanisms and roles for TAs in stress adaptation and persistence of Cronobacter within food matrices and food processing environments.

Occurrence of virulence factors in Cronobacter sakazakii and Cronobacter malonaticus originated from clinical samples

BACKGROUND

Cronobacter spp. are Gram-negative, facultative-anaerobic, non-spore forming, enteric coliform bacteria, which belongs to the Enterobacteriaceae family. Cronobacter spp. are opportunistic pathogens that have brought rare but life-threatening infections such as meningitis, necrotizing enterocolitis and bloodstream infections in neonates and infants. Information on the diversity, pathogenicity and virulence of Cronobacter species obtained from various sources is still relatively scarce and fragmentary. The aim of this study was to examine and analyse different pathogenicity and virulence factors among C. sakazakii and C. malonaticus strains isolated from clinical samples.

METHODS

The thirty-six clinical Cronobacter strains have been used in this study. This bacterial collection consists of 25 strains of C. sakazakii and 11 strains of C. malonaticus, isolated from different clinical materials. Seven genes (ompA, inv, sip, aut, hly, fliC, cpa) were amplified by PCR. Moreover, the motility and the ability of these strains to adhere and invade human colorectal adenocarcinoma (HT-29) and mouse neuroblastoma (N1E-115) cell lines were investigated.

RESULTS

Our results showed that all tested strains were able to adhere to both used cell lines, HT-29 and N1E-115 cells. The invasion assay showed that 66.7% (24/36) of isolates were able to invade N1-E115 cells while 83% (30/36) of isolates were able to invade HT-29 cells. On the average, 68% of the C. sakazakii strains exhibited seven virulence factors and only 18% in C. malonaticus. All strains amplified ompA and fliC genes. The other genes were detected as follow: sip 97% (35/36), hlyA 92% (33/36), aut 94% (34/36), cpa 67% (24/36), and inv 69% (25/36).

CONCLUSIONS

C. sakazakii and C malonaticus strains demonstrate the diversity of the virulence factors present among these pathogens. It is necessary to permanently monitor the hospital environment to appropriately treat and resolve cases associated with disease. Furthermore, in-depth knowledge is needed about the source and transmission vehicles of pathogens in hospitals to adopt pertinent prevention measures.

Draft genomes of Cronobacter sakazakii strains isolated from dried spices bring unique insights into the diversity of plant-associated strains

Cronobacter sakazakii is a Gram-negative opportunistic pathogen that causes life- threatening infantile infections, such as meningitis, septicemia, and necrotizing enterocolitis, as well as pneumonia, septicemia, and urinary tract and wound infections in adults. Here, we report 26 draft genome sequences of C. sakazakii, which were obtained from dried spices from the USA, the Middle East, China, and the Republic of Korea. The average genome size of the C. sakazakii genomes was 4393 kb, with an average of 4055 protein coding genes, and an average genome G + C content of 56.9%. The genomes contained genes related to carbohydrate transport and metabolism, amino acid transport and metabolism, and cell wall/membrane biogenesis. In addition, we identified genes encoding proteins involved in osmotic responses such as DnaJ, Aquaproin Z, ProQ, and TreF, as well as virulence-related and heat shock-related proteins. Interestingly, a metabolic island comprised of a variably-sized xylose utilization operon was found within the spice-associated C. sakazakii genomes, which supports the hypothesis that plants may serve as transmission vectors or alternative hosts for Cronobacter species. The presence of the genes identified in this study can support the remarkable phenotypic traits of C. sakazakii such as the organism's capabilities of adaptation and survival in response to adverse growth environmental conditions (e.g. osmotic and desiccative stresses). Accordingly, the genome analyses provided insights into many aspects of physiology and evolutionary history of this important foodborne pathogen.

Reconstituting the History of Cronobacter Evolution Driven by Differentiated CRISPR Activity

Cronobacter strains harboring the CRISPR-Cas system are important foodborne pathogens causing serious neonatal infections. However, the specific role of the CRISPR-Cas system in bacterial evolution remains relatively unexplored. In this study, we investigated the impact of the CRISPR-Cas system on Cronobacter evolution and obtained 137 new whole-genome Cronobacter sequences by next-generation sequencing technology. Among the strains examined (n = 240), 90.6% (193/213) of prevalent species Cronobacter sakazakii, Cronobactermalonaticus, and Cronobacterdublinensis strains had intact CRISPR-Cas systems. Two rare species, Cronobactercondimenti (n = 2) and Cronobacteruniversalis (n = 6), lacked and preserved the CRISPR-Cas system at a low frequency (1/6), respectively. These results suggest that the presence of one CRISPR-Cas system is important for a Cronobacter species to maintain genome homeostasis for survival. The Cronobacter ancestral strain is likely to have harbored both subtype I-E and I-F CRISPR-Cas systems; during the long evolutionary process, subtype I-E was retained while subtype I-F selectively degenerated in Cronobacter species and was even lost by the major Cronobacter pathovars. Moreover, significantly higher CRISPR activity was observed in the plant-associated species Cdublinensis than in the virulence-related species C. sakazakii and Cmalonaticus Similar spacers of CRISPR arrays were rarely found among species, suggesting intensive change through adaptive acquisition and loss. Differentiated CRISPR activity appears to be the product of environmental selective pressure and might contribute to the bidirectional divergence and speciation of CronobacterIMPORTANCE This study reports the evolutionary history of Cronobacter under the selective pressure of the CRISPR-Cas system. One CRISPR-Cas system in Cronobacter is important for maintaining genome homeostasis, whereas two types of systems may be redundant and not conducive to acquiring beneficial DNA for environmental adaptation and pathogenicity. Differentiated CRISPR activity has contributed to the bidirectional divergence and genetic diversity of Cronobacter This perspective makes a significant contribution to the literature by providing new insights into CRISPR-Cas systems in general, while further expanding the roles of CRISPR beyond conferring adaptive immunity and demonstrating a link to adaptation and species divergence in a genus. Moreover, our study provides new insights into the balance between genome homeostasis and the uptake of beneficial DNA related to CRISPR-based activity in the evolution of Cronobacter.

Genomic characterization of malonate positive Cronobacter sakazakii serotype O:2, sequence type 64 strains, isolated from clinical, food, and environment samples

Background

Malonate utilization, an important differential trait, well recognized as being possessed by six of the seven Cronobacter species is thought to be largely absent in Cronobacter sakazakii (Csak). The current study provides experimental evidence that confirms the presence of a malonate utilization operon in 24 strains of sequence type (ST) 64, obtained from Europe, Middle East, China, and USA; it offers explanations regarding the genomic diversity and phylogenetic relatedness among these strains, and that of other C. sakazakii strains.

Results

In this study, the presence of a malonate utilization operon in these strains was initially identified by DNA microarray analysis (MA) out of a pool of 347 strains obtained from various surveillance studies involving clinical, spices, milk powder sources and powdered infant formula production facilities in Ireland and Germany, and dried dairy powder manufacturing facilities in the USA. All ST64 C. sakazakii strains tested could utilize malonate. Zebrafish embryo infection studies showed that C. sakazakii ST64 strains are as virulent as other Cronobacter species. Parallel whole genome sequencing (WGS) and MA showed that the strains phylogenetically grouped as a separate clade among the Csak species cluster. Additionally, these strains possessed the Csak O:2 serotype. The nine-gene, ~ 7.7 kbp malonate utilization operon was located in these strains between two conserved flanking genes, gyrB and katG. Plasmidotyping results showed that these strains possessed the virulence plasmid pESA3, but in contrast to the USA ST64 Csak strains, ST64 Csak strains isolated from sources in Europe and the Middle East, did not possess the type six secretion system effector vgrG gene.

Conclusions

Until this investigation, the presence of malonate-positive Csak strains, which are associated with foods and clinical cases, was under appreciated. If this trait was used solely to identify Cronobacter strains, many strains would likely be misidentified. Parallel WGS and MA were useful in characterizing the total genome content of these Csak O:2, ST64, malonate-positive strains and further provides an understanding of their phylogenetic relatedness among other virulent C. sakazakii strains.

Electronic supplementary material

The online version of this article (10.1186/s13099-018-0238-9) contains supplementary material, which is available to authorized users.

DNA-Sequence Based Typing of the Cronobacter Genus Using MLST, CRISPR-cas Array and Capsular Profiling

The Cronobacter genus is composed of seven species, within which a number of pathovars have been described. The most notable infections by Cronobacter spp. are of infants through the consumption of contaminated infant formula. The description of the genus has greatly improved in recent years through DNA sequencing techniques, and this has led to a robust means of identification. However some species are highly clonal and this limits the ability to discriminate between unrelated strains by some methods of genotyping. This article updates the application of three genotyping methods across the Cronobacter genus. The three genotyping methods were multilocus sequence typing (MLST), capsular profiling of the K-antigen and colanic acid (CA) biosynthesis regions, and CRISPR-cas array profiling. A total of 1654 MLST profiled and 286 whole genome sequenced strains, available by open access at the PubMLST Cronobacter database, were used this analysis. The predominance of C. sakazakii and C. malonaticus in clinical infections was confirmed. The majority of clinical strains being in the C. sakazakii clonal complexes (CC) 1 and 4, sequence types (ST) 8 and 12 and C. malonaticus ST7. The capsular profile K2:CA2, previously proposed as being strongly associated with C. sakazakii and C. malonaticus isolates from severe neonatal infections, was also found in C. turicensis, C. dublinensis and C. universalis. The majority of CRISPR-cas types across the genus was the I-E (Ecoli) type. Some strains of C. dublinensis and C. muytjensii encoded the I-F (Ypseudo) type, and others lacked the cas gene loci. The significance of the expanding profiling will be of benefit to researchers as well as governmental and industrial risk assessors.

Multiplexed Lateral Flow Test for Detection and Differentiation of Cronobacter sakazakii Serotypes O1 and O2

The ubiquitous and opportunistic pathogen Cronobacter sakazakii is responsible for severe meningitis, sepsis, and necrotizing enterocolitis in neonates and infants associated with ingestion of contaminated powdered infant formula (PIF). The current ISO method for isolation and detection of Cronobacter spp. is laborious, time-consuming and expensive. In this study, a multiplexed lateral flow test strip was developed to rapidly detect and simultaneously serotype O1 and O2 C. sakazakii serotypes. The assay is based on two monoclonal antibodies (MAb) that specifically bind to the lipopolysaccharides (LPS) of these pathogens. The test strip provides results very quickly; C. sakazakii could be detected in pure culture within 15 min with a sensitivity of 107 CFU/ml. After non-selective enrichment for 18 h as low as one Cronobacter cell per g PIF could be detected. Moreover, the established lateral flow assay (LFA) offers excellent specificity showing no cross-reactivity with other C. sakazakii serotypes, Cronobacter species or Enterobacteriaceae tested. These characteristics, together with several advantages such as speed, simplicity in performance, low analysis cost, and no requirement of specialized skills or sophisticated equipment make the developed multiplexed LFA suitable for reliable detection and serotyping of C. sakazakii serotypes O1 and O2.

Metagenetic analysis of the bacterial communities of edible insects from diverse production cycles at industrial rearing companies

Despite the continuing development of new insect-derived food products, microbial research on edible insects and insect-based foods is still very limited. The goal of this study was to increase the knowledge on the microbial quality of edible insects by comparing the bacterial community composition of mealworms (Tenebrio molitor) and crickets (Acheta domesticus and Gryllodes sigillatus) from several production cycles and rearing companies. Remarkable differences in the bacterial community composition were found between different mealworm rearing companies and mealworm production cycles from the same company. In comparison with mealworms, the bacterial community composition of the investigated crickets was more similar among different companies, and was highly similar between both cricket species investigated. Mealworm communities were dominated by Spiroplasma and Erwinia species, while crickets were abundantly colonised by (Para)bacteroides species. With respect to food safety, only a few operational taxonomic units could be associated with potential human pathogens such as Cronobacter or spoilage bacteria such as Pseudomonas. In summary, our results implicate that at least for cricket rearing, production cycles of constant and good quality in terms of bacterial composition can be obtained by different rearing companies. For mealworms however, more variation in terms of microbial quality occurs between companies.

Draft Genome Sequence of Cronobacter sakazakii GP1999, Sequence Type 145, an Epiphytic Isolate Obtained from the Tomato's Rhizoplane/Rhizosphere Continuum

We present here the draft genome of Cronobacter sakazakii GP1999, a sequence type 145 strain isolated from the rhizosphere of tomato plants. Assembly and annotation of the genome resulted in a genome of 4,504,670 bp in size, with 4,148 coding sequences, and a GC content of 56.8%.

Use of a Pan-Genomic DNA Microarray in Determination of the Phylogenetic Relatedness among Cronobacter spp. and Its Use as a Data Mining Tool to Understand Cronobacter Biology

Cronobacter (previously known as Enterobacter sakazakii) is a genus of Gram-negative, facultatively anaerobic, oxidase-negative, catalase-positive, rod-shaped bacteria of the family Enterobacteriaceae. These organisms cause a variety of illnesses such as meningitis, necrotizing enterocolitis, and septicemia in neonates and infants, and urinary tract, wound, abscesses or surgical site infections, septicemia, and pneumonia in adults. The total gene content of 379 strains of Cronobacter spp. and taxonomically-related isolates was determined using a recently reported DNA microarray. The Cronobacter microarray as a genotyping tool gives the global food safety community a rapid method to identify and capture the total genomic content of outbreak isolates for food safety, environmental, and clinical surveillance purposes. It was able to differentiate the seven Cronobacter species from one another and from non-Cronobacter species. The microarray was also able to cluster strains within each species into well-defined subgroups. These results also support previous studies on the phylogenic separation of species members of the genus and clearly highlight the evolutionary sequence divergence among each species of the genus compared to phylogenetically-related species. This review extends these studies and illustrates how the microarray can also be used as an investigational tool to mine genomic data sets from strains. Three case studies describing the use of the microarray are shown and include: (1) the determination of allelic differences among Cronobacter sakazakii strains possessing the virulence plasmid pESA3; (2) mining of malonate and myo-inositol alleles among subspecies of Cronobacter dublinensis strains to determine subspecies identity; and (3) lastly using the microarray to demonstrate sequence divergence and phylogenetic relatedness trends for 13 outer-membrane protein alleles among 240 Cronobacter and phylogenetically-related strains. The goal of this review is to describe microarrays as a robust tool for genomics research of this assorted and important genus, a criterion toward the development of future preventative measures to eliminate this foodborne pathogen from the global food supply.

Comparative Genome Analyses of Vibrio anguillarum Strains Reveal a Link with Pathogenicity Traits

Comparative genome analysis of strains of a pathogenic bacterial species can be a powerful tool to discover acquisition of mobile genetic elements related to virulence. Here, we compared 28 V. anguillarum strains that differed in virulence in fish larval models. By pan-genome analyses, we found that six of nine highly virulent strains had a unique core and accessory genome. In contrast, V. anguillarum strains that were medium to nonvirulent had low genomic diversity. Integration of genomic and phenotypic features provides insights into the evolution of V. anguillarum and can also be important for survey and diagnostic purposes.

Vibrio anguillarum is a marine bacterium that can cause vibriosis in many fish and shellfish species, leading to high mortalities and economic losses in aquaculture. Although putative virulence factors have been identified, the mechanism of pathogenesis of V. anguillarum is not fully understood. Here, we analyzed whole-genome sequences of a collection of V. anguillarum strains and compared them to virulence of the strains as determined in larval challenge assays. Previously identified virulence factors were globally distributed among the strains, with some genetic diversity. However, the pan-genome revealed that six out of nine high-virulence strains possessed a unique accessory genome that was attributed to pathogenic genomic islands, prophage-like elements, virulence factors, and a new set of gene clusters involved in biosynthesis, modification, and transport of polysaccharides. In contrast, V. anguillarum strains that were medium to nonvirulent had a high degree of genomic homogeneity. Finally, we found that a phylogeny based on the core genomes clustered the strains with moderate to no virulence, while six out of nine high-virulence strains represented phylogenetically separate clusters. Hence, we suggest a link between genotype and virulence characteristics of Vibrio anguillarum, which can be used to unravel the molecular evolution of V. anguillarum and can also be important from survey and diagnostic perspectives.

IMPORTANCE Comparative genome analysis of strains of a pathogenic bacterial species can be a powerful tool to discover acquisition of mobile genetic elements related to virulence. Here, we compared 28 V. anguillarum strains that differed in virulence in fish larval models. By pan-genome analyses, we found that six of nine highly virulent strains had a unique core and accessory genome. In contrast, V. anguillarum strains that were medium to nonvirulent had low genomic diversity. Integration of genomic and phenotypic features provides insights into the evolution of V. anguillarum and can also be important for survey and diagnostic purposes.

Analysis and Characterization of Proteins Associated with Outer Membrane Vesicles Secreted by Cronobacter spp

Little is known about secretion of outer membrane vesicles (OMVs) by Cronobacter. In this study, OMVs isolated from Cronobacter sakazakii, Cronobacter turicensis, and Cronobacter malonaticus were examined by electron microscopy (EM) and their associated outer membrane proteins (OMP) and genes were analyzed by SDS-PAGE, protein sequencing, BLAST, PCR, and DNA microarray. EM of stained cells revealed that the OMVs are secreted as pleomorphic micro-vesicles which cascade from the cell's surface. SDS-PAGE analysis identified protein bands with molecular weights of 18 kDa to >100 kDa which had homologies to OMPs such as GroEL; OmpA, C, E, F, and X; MipA proteins; conjugative plasmid transfer protein; and an outer membrane auto-transporter protein (OMATP). PCR analyses showed that most of the OMP genes were present in all seven Cronobacter species while a few genes (OMATP gene, groEL, ompC, mipA, ctp, and ompX) were absent in some phylogenetically-related species. Microarray analysis demonstrated sequence divergence among the OMP genes that was not captured by PCR. These results support previous findings that OmpA and OmpX may be involved in virulence of Cronobacter, and are packaged within secreted OMVs. These results also suggest that other OMV-packaged OMPs may be involved in roles such as stress response, cell wall and plasmid maintenance, and extracellular transport.

The driving force of prophages and CRISPR-Cas system in the evolution of Cronobacter sakazakii

Cronobacter sakazakii is an important foodborne pathogens causing rare but life-threatening diseases in neonates and infants. CRISPR-Cas system is a new prokaryotic defense system that provides adaptive immunity against phages, latter play an vital role on the evolution and pathogenicity of host bacteria. In this study, we found that genome sizes of C. sakazakii strains had a significant positive correlation with total genome sizes of prophages. Prophages contributed to 16.57% of the genetic diversity (pan genome) of C. sakazakii, some of which maybe the potential virulence factors. Subtype I-E CRISPR-Cas system and five types of CRISPR arrays were found in the conserved site of C. sakazakii strains. CRISPR1 and CRISPR2 loci with high variable spacers were active and showed potential protection against phage attacks. The number of spacers from two active CRISPR loci in clinical strains was significant less than that of foodborne strains, it maybe a reason why clinical strains were found to have more prophages than foodborne strains. The frequently gain/loss of prophages and spacers in CRISPR loci is likely to drive the quick evolution of C. sakazakii. Our study provides a new insight into the co-evolution of phages and C. sakazakii.

Transcriptome analysis of Cronobacter sakazakii ATCC BAA-894 after interaction with human intestinal epithelial cell line HCT-8

Cronobacter spp. are opportunistic pathogens that are responsible for infections including severe meningitis, septicemia, and necrotizing enterocolitis in neonates and infants. To date, questions still remain regarding the mechanisms of pathogenicity and virulence determinants for each bacterial strain. In this study, we established an in vitro model for Cronobacter sakazakii ATCC BAA-894 infection of HCT-8 human colorectal epithelial cells. The transcriptome profile of C. sakazakii ATCC BAA-894 after interaction with HCT-8 cells was determined using high-throughput whole-transcriptome sequencing (RNA sequencing (RNA-seq)). Gene expression profiles indicated that 139 genes were upregulated and 72 genes were downregulated in the adherent C. sakazakii ATCC BAA-894 strain on HCT-8 cells compared to the cultured bacteria in the cell-free medium. Expressions of some flagella genes and virulence factors involved in adherence were upregulated. High osmolarity and osmotic stress-associated genes were highly upregulated, as well as genes responsible for the synthesis of lipopolysaccharides and outer membrane proteins, iron acquisition systems, and glycerol and glycerophospholipid metabolism. In sum, our study provides further insight into the mechanisms underlying C. sakazakii pathogenesis in the human gastrointestinal tract.

Linking Genomo- and Pathotype: Exploiting the Zebrafish Embryo Model to Investigate the Divergent Virulence Potential among Cronobacter spp

Bacteria belonging to the genus Cronobacter have been recognized as causative agents of life-threatening systemic infections primarily in premature, low-birth weight and immune-compromised neonates. Apparently not all Cronobacter species are linked to infantile infections and it has been proposed that virulence varies among strains. Whole genome comparisons and in silico analysis have proven to be powerful tools in elucidating potential virulence determinants, the presence/absence of which may explain the differential virulence behaviour of strains. However, validation of these factors has in the past been hampered by the availability of a suitable neonatal animal model. In the present study we have used zebrafish embryos to model Cronobacter infections in vivo using wild type and genetically engineered strains. Our experiments confirmed the role of the RepF1B-like plasmids as “virulence plasmids” in Cronobacter and underpinned the importantce of two putative virulence factors—cpa and zpx—in in vivo pathogenesis. We propose that by using this model in vivo infection studies are now possible on a large scale level which will boost the understanding on the virulence strategies employed by these pathogens.

Inside the Pan-genome - Methods and Software Overview

The number of genomes that have been deposited in databases has increased exponentially after the advent of Next-Generation Sequencing (NGS), which produces high-throughput sequence data; this circumstance has demanded the development of new bioinformatics software and the creation of new areas, such as comparative genomics. In comparative genomics, the genetic content of an organism is compared against other organisms, which helps in the prediction of gene function and coding region sequences, identification of evolutionary events and determination of phylogenetic relationships. However, expanding comparative genomics to a large number of related bacteria, we can infer their lifestyles, gene repertoires and minimal genome size. In this context, a powerful approach called Pan-genome has been initiated and developed. This approach involves the genomic comparison of different strains of the same species, or even genus. Its main goal is to establish the total number of non-redundant genes that are present in a determined dataset. Pan-genome consists of three parts: core genome; accessory or dispensable genome; and species-specific or strain-specific genes. Furthermore, pan-genome is considered to be “open” as long as new genes are added significantly to the total repertoire for each new additional genome and “closed” when the newly added genomes cannot be inferred to significantly increase the total repertoire of the genes. To perform all of the required calculations, a substantial amount of software has been developed, based on orthologous and paralogous gene identification.

Rapid detection of Cronobacter sakazakii by real-time PCR based on the cgcA gene and TaqMan probe with internal amplification control

Cronobacter sakazakii is a severe virulent strain that is frequently detected in powdered infant formula (PIF). Therefore, it is necessary to develop a fast and specific detection method. The specificity of our newly developed quantitative real-time PCR (qRT–PCR) was validated with DNA from 46 strains. Among them, 12 C. sakazakii strains were correctly amplified, whereas no positive florescent signal was observed from 34 nontarget controls. The detection limit of C. sakazakii was about 110 CFU/mL in broth and 1100 CFU/g in PIF. After enrichment in buffered peptone water for 6 h, our developed qRT–PCR assay could reliably detect C. sakazakii when the inoculation level was as low as 2 CFU/25 g (0.08 CFU/g) in PIF. The growth of C. sakazakii could be inhibited by the presence of Lactobacillus pentosus and Bacillus cereus, which used a longer enrichment period before the isolation was accomplished. However, at 5 and 50 CFU/25 g inoculation levels of C. sakazakii in the presence of 4 × 106 CFU L. pentosus/25 g or of 2 × 104 CFU B. cereus/25 g, the qRT–PCR assay could detect the presence of Cronobacter even though these artificially spiked samples were negative in culture. Therefore, our results indicated that the qRT–PCR assay could detect samples containing inhibitors and could avoid false negatives by using an internal amplification control.

Comparative Genomics of Environmental and Clinical Stenotrophomonas maltophilia Strains with Different Antibiotic Resistance Profiles

Stenotrophomonas maltophilia, a ubiquitous Gram-negative γ-proteobacterium, has emerged as an important opportunistic pathogen responsible for nosocomial infections. A major characteristic of clinical isolates is their high intrinsic or acquired antibiotic resistance level. The aim of this study was to decipher the genetic determinism of antibiotic resistance among strains from different origins (i.e., natural environment and clinical origin) showing various antibiotic resistance profiles. To this purpose, we selected three strains isolated from soil collected in France or Burkina Faso that showed contrasting antibiotic resistance profiles. After whole-genome sequencing, the phylogenetic relationships of these 3 strains and 11 strains with available genome sequences were determined. Results showed that a strain's phylogeny did not match their origin or antibiotic resistance profiles. Numerous antibiotic resistance coding genes and efflux pump operons were revealed by the genome analysis, with 57% of the identified genes not previously described. No major variation in the antibiotic resistance gene content was observed between strains irrespective of their origin and antibiotic resistance profiles. Although environmental strains generally carry as many multidrug resistant (MDR) efflux pumps as clinical strains, the absence of resistance-nodulation-division (RND) pumps (i.e., SmeABC) previously described to be specific to S. maltophilia was revealed in two environmental strains (BurA1 and PierC1). Furthermore the genome analysis of the environmental MDR strain BurA1 showed the absence of SmeABC but the presence of another putative MDR RND efflux pump, named EbyCAB on a genomic island probably acquired through horizontal gene transfer.

Analysis of the cellulose synthase operon genes, bcsA, bcsB, and bcsC in Cronobacter species: Prevalence among species and their roles in biofilm formation and cell-cell aggregation

Cronobacter species are emerging food-borne pathogens that cause severe sepsis, meningitis, and necrotizing entercolitis in neonates and infants. Bacterial pathogens such as Escherichia coli and Salmonella species produce extracellular cellulose which has been shown to be involved in rugosity, biofilm formation, and host colonization. In this study the distribution and prevalence of cellulose synthase operon genes (bcsABZC) were determined by polymerase chain reaction (PCR) analysis in 231 Cronobacter strains isolated from clinical, food, environmental, and unknown sources. Furthermore, bcsA and bcsB isogenic mutants were constructed in Cronobacter sakazakii BAA894 to determine their roles. In calcofluor binding assays bcsA and bcsB mutants did not produce cellulose, and their colonial morphotypes were different to that of the parent strain. Biofilm formation and bacterial cell-cell aggregation were significantly reduced in bcsA and bcsB mutants compared to the parental strain. bcsA or bcsAB PCR-negative strains of C. sakazakii did not bind calcofluor, and produced less biofilm and cell-cell aggregation compared to strains possessing bcsAB genes. These data indicated that Cronobacter bcsABZC were present in all clinical isolates and most of food and environmental isolates. bcsA and bcsB genes of Cronobacter were necessary to produce cellulose, and were involved in biofilm formation and cell-cell aggregation.

Comparative Genotypic and Phenotypic Analysis of Cronobacter Species Cultured from Four Powdered Infant Formula Production Facilities: Indication of Pathoadaptation along the Food Chain

Cronobacter species are opportunistic pathogens commonly found in the environment. Among the seven Cronobacter species, Cronobacter sakazakii sequence type 4 (ST-4) is predominantly associated with recorded cases of infantile meningitis. This study reports on a 26-month powdered infant formula (PIF) surveillance program in four production facilities located in distinct geographic regions. The objective was to identify the ST(s) in PIF production environments and to investigate the phenotypic features that support their survival. Of all 168 Cronobacter isolates, 133 were recovered from a PIF production environment, 31 were of clinical origin, and 4 were laboratory type strains. Sequence type 1 (n = 84 isolates; 63.9%) was the dominant type in PIF production environments. The majority of these isolates clustered with an indistinguishable pulsotype and persisted for at least an 18-month period. Moreover, DNA microarray results identified two phylogenetic lineages among ST-4 strains tested. Thereafter, the ST-1 and -4 isolates were phenotypically compared. Differences were noted based on the phenotypes expressed by these isolates. The ST-1 PIF isolates produced stronger biofilms at both 28°C and 37°C, while the ST-4 clinical isolates exhibited greater swimming activity and increased binding to Congo red dye. Given the fact that PIF is a low-moisture environment and that the clinical environment provides for an interaction between the pathogen and its host, these differences may be consistent with a form of pathoadaptation. These findings help to extend our current understanding of the epidemiology and ecology of Cronobacter species in PIF production environments.

Genome Sequence of Cronobacter sakazakii Serogroup O:4, Sequence Type 4 Strain CDC 2009-03746, Isolated from a Fatal Case of Infantile Meningitis

We report the draft genome sequence of a Cronobacter sakazakii serogroup O:4, sequence type 4 strain, CDC 2009-03746 (=NM1240=2009-06-01), isolated from a fatal case of infantile meningitis. The draft genome has a size of 4,492,904 bp and a G+C% content of 56.7.