Rapid identification of enteric bacteria from whole genome sequences using average nucleotide identity metrics.

Identification of enteric bacteria species by whole genome sequence (WGS) analysis requires a rapid and an easily standardized approach. We leveraged the principles of average nucleotide identity using MUMmer (ANIm) software, which calculates the percent bases aligned between two bacterial genomes and their corresponding ANI values, to set threshold values for determining species consistent with the conventional identification methods of known species. The performance of species identification was evaluated using two datasets: the Reference Genome Dataset v2 (RGDv2), consisting of 43 enteric genome assemblies representing 32 species, and the Test Genome Dataset (TGDv1), comprising 454 genome assemblies which is designed to represent all species needed to query for identification, as well as rare and closely related species. The RGDv2 contains six Campylobacter spp., three Escherichia/Shigella spp., one Grimontia hollisae, six Listeria spp., one Photobacterium damselae, two Salmonella spp., and thirteen Vibrio spp., while the TGDv1 contains 454 enteric bacterial genomes representing 42 different species. The analysis showed that, when a standard minimum of 70% genome bases alignment existed, the ANI threshold values determined for these species were ≥95 for Escherichia/Shigella and Vibrio species, ≥93% for Salmonella species, and ≥92% for Campylobacter and Listeria species. Using these metrics, the RGDv2 accurately classified all validation strains in TGDv1 at the species level, which is consistent with the classification based on previous gold standard methods.

Identification and Characterization of ten Escherichia coli Strains Encoding Novel Shiga Toxin 2 Subtypes, Stx2n as Well as Stx2j, Stx2m, and Stx2o, in the United States.

The sharing of genome sequences in online data repositories allows for large scale analyses of specific genes or gene families. This can result in the detection of novel gene subtypes as well as the development of improved detection methods. Here, we used publicly available WGS data to detect a novel Stx subtype, Stx2n in two clinical E. coli strains isolated in the USA. During this process, additional Stx2 subtypes were detected; six Stx2j, one Stx2m strain, and one Stx2o, were all analyzed for variability from the originally described subtypes. Complete genome sequences were assembled from short- or long-read sequencing and analyzed for serotype, and ST types. The WGS data from Stx2n- and Stx2o-producing STEC strains were further analyzed for virulence genes pro-phage analysis and phage insertion sites. Nucleotide and amino acid maximum parsimony trees showed expected clustering of the previously described subtypes and a clear separation of the novel Stx2n subtype. WGS data were used to design OMNI PCR primers for the detection of all known stx1 (283 bp amplicon), stx2 (400 bp amplicon), intimin encoded by eae (221 bp amplicon), and stx2f (438 bp amplicon) subtypes. These primers were tested in three different laboratories, using standard reference strains. An analysis of the complete genome sequence showed variability in serogroup, virulence genes, and ST type, and Stx2 pro-phages showed variability in size, gene composition, and phage insertion sites. The strains with Stx2j, Stx2m, Stx2n, and Stx2o showed toxicity to Vero cells. Stx2j carrying strain, 2012C-4221, was induced when grown with sub-inhibitory concentrations of ciprofloxacin, and toxicity was detected. Taken together, these data highlight the need to reinforce genomic surveillance to identify the emergence of potential new Stx2 or Stx1 variants. The importance of this surveillance has a paramount impact on public health. Per our description in this study, we suggest that 2017C-4317 be designated as the Stx2n type-strain.

Reoccurring Escherichia coli O157:H7 Strain Linked to Leafy Greens-Associated Outbreaks, 2016-2019.

Genomic characterization of an Escherichia coli O157:H7 strain linked to leafy greens-associated outbreaks dates its emergence to late 2015. One clade has notable accessory genomic content and a previously described mutation putatively associated with increased arsenic tolerance. This strain is a reoccurring, emerging, or persistent strain causing illness over an extended period.

O-Antigen Diversification Masks Identification of Highly Pathogenic Shiga Toxin-Producing Escherichia coli O104:H4-Like Strains.

Shiga toxin-producing Escherichia coli (STEC) can give rise to a range of clinical outcomes from diarrhea to the life-threatening systemic condition hemolytic-uremic syndrome (HUS). Although STEC O157:H7 is the serotype most frequently associated with HUS, a major outbreak of HUS occurred in 2011 in Germany and was caused by a rare serotype, STEC O104:H4. Prior to 2011 and since the outbreak, STEC O104:H4 strains have only rarely been associated with human infections. From 2012 to 2020, intensified STEC surveillance was performed in Germany where the subtyping of ~8,000 clinical isolates by molecular methods, including whole-genome sequencing, was carried out. A rare STEC serotype, O181:H4, associated with HUS was identified, and like the STEC O104:H4 outbreak strain, this strain belongs to sequence type 678 (ST678). Genomic and virulence comparisons revealed that the two strains are phylogenetically related and differ principally in the gene cluster encoding their respective lipopolysaccharide O-antigens but exhibit similar virulence phenotypes. In addition, five other serotypes belonging to ST678 from human clinical infection, such as OX13:H4, O127:H4, OgN-RKI9:H4, O131:H4, and O69:H4, were identified from diverse locations worldwide. IMPORTANCE Our data suggest that the high-virulence ensemble of the STEC O104:H4 outbreak strain remains a global threat because genomically similar strains cause disease worldwide but that the horizontal acquisition of O-antigen gene clusters has diversified the O-antigens of strains belonging to ST678. Thus, the identification of these highly pathogenic strains is masked by diverse and rare O-antigens, thereby confounding the interpretation of their potential risk.

Cronobacter sakazakii Infections in Two Infants Linked to Powdered Infant Formula and Breast Pump Equipment - United States, 2021 and 2022.

Cronobacter sakazakii, a species of gram-negative bacteria belonging to the Enterobacteriaceae family, is known to cause severe and often fatal meningitis and sepsis in young infants. C. sakazakii is ubiquitous in the environment, and most reported infant cases have been attributed to contaminated powdered infant formula (powdered formula) or breast milk that was expressed using contaminated breast pump equipment (1-3). Previous investigations of cases and outbreaks have identified C. sakazakii in opened powdered formula, breast pump parts, environmental surfaces in the home, and, rarely, in unopened powdered formula and formula manufacturing facilities (2,4-6). This report describes two infants with C. sakazakii meningitis reported to CDC in September 2021 and February 2022. CDC used whole genome sequencing (WGS) analysis to link one case to contaminated opened powdered formula from the patient's home and the other to contaminated breast pump equipment. These cases highlight the importance of expanding awareness about C. sakazakii infections in infants, safe preparation and storage of powdered formula, proper cleaning and sanitizing of breast pump equipment, and using WGS as a tool for C. sakazakii investigations.

Functional analysis of Niemann-Pick disease type C family protein, NPC1a, in Drosophila melanogaster.

During embryonic gonad coalescence, primordial germ cells (PGCs) follow a carefully choreographed migratory route circumscribed by guidance signals towards somatic gonadal precursor cells (SGPs). In Drosophila melanogaster, SGP-derived Hedgehog (Hh), which serves as a guidance cue for the PGCs, is potentiated by mesodermally restricted HMGCoA-reductase (Hmgcr) and the ABC transporter Multi-drug-resistant-49 (Mdr49). Given the importance of cholesterol modification in the processing and long-distance transmission of the Hh ligand, we have analyzed the involvement of the Niemann-Pick disease type C-1a (NPC1a) protein, a cholesterol transporter, in germ cell migration and Hedgehog signaling. We show that mesoderm-specific inactivation of Npc1a results in germ cell migration defects. Similar to Mdr49, PGC migration defects in the Npc1a embryos are ameliorated by a cholesterol-rich diet. Consistently, reduction in Npc1a weakens the ability of ectopic HMG Coenzyme A reductase (Hmgcr) to induce germ cell migration defects. Moreover, compromising Npc1a levels influences Hh signaling adversely during wing development, a process that relies upon long-range Hh signaling. Last, doubly heterozygous embryos (Mdr49/Npc1a) display enhanced germ cell migration defects when compared with single mutants (Npc1a/+ or Mdr49/+), supporting cooperative interaction between the two.

PacBio Genome Sequences of Eight Escherichia albertii Strains Isolated from Humans in the United States.

Escherichia albertii is an emerging pathogen that is closely related to Escherichia coli and can carry some of the same virulence genes as E. coli. Here, we report the release of Illumina-corrected PacBio sequences for eight E. albertii genomes. Two of these strains carry Shiga toxin 2f.

PacBio Genome Sequences of Escherichia coli Serotype O157:H7, Diffusely Adherent E. coli, and Salmonella enterica Strains, All Carrying Plasmids with an mcr-1 Resistance Gene.

We report here Illumina-corrected PacBio whole-genome sequences of an Escherichia coli serotype O157:H7 strain (2017C-4109), an E. coli serotype O[undetermined]:H2 strain (2017C-4173W12), and a Salmonella enterica subsp. enterica serovar Enteritidis strain (2017K-0021), all of which carried the mcr-1 resistance gene on an IncI2 or IncX4 plasmid. We also determined that pMCR-1-CTSe is identical to a previously published plasmid, pMCR-1-CT.

Interlaboratory Evaluation of the U.S. Food and Drug Administration Escherichia coli Identification Microarray for Profiling Shiga Toxin-Producing Escherichia coli.

The U.S. Food and Drug Administration Escherichia coli Identification (FDA-ECID) microarray provides rapid molecular characterization of E. coli. The effectiveness of the FDA-ECID for characterizing Shiga toxin-producing E. coli (STEC) was evaluated by three federal laboratories and one reference laboratory with a panel of 54 reference E. coli strains from the External Quality Assurance program. Strains were tested by FDA-ECID for molecular serotyping (O and H antigens), Shiga toxin subtyping, and the presence of the ehxA and eae genes for enterohemolysin and intimin, respectively. The FDA-ECID O typing was 96% reproducible among the four laboratories and 94% accurate compared with the reference External Quality Assurance data. Discrepancies were due to the absence of O41 target loci on the array and to two pairs of O types with identical target sequences. H typing was 96% reproducible and 100% accurate, with discrepancies due to two strains from one laboratory that were identified as mixed by FDA-ECID. Shiga toxin (Stx) type 1 subtyping was 100% reproducible and accurate, and Stx2 subtyping was 100% reproducible but only 64% accurate. FDA-ECID identified most Stx2 subtypes but had difficulty distinguishing among stx2a, stx2c, and stx2d genes because of close similarities of these sequences. FDA-ECID was 100% effective for detecting ehxA and eae and accurately subtyped the eae alleles. This interlaboratory study revealed that FDA-ECID for STEC characterization was highly reproducible for molecular serotyping, stx and eae subtyping, and ehxA detection. However, the array was less useful for distinguishing among the highly homologous O antigen genes and the stx2a, stx2c, and stx2d subtypes.

High-Quality Whole-Genome Sequences for 59 Historical Shigella Strains Generated with PacBio Sequencing.

Shigella spp. are enteric pathogens that cause shigellosis. We report here the high-quality whole-genome sequences of 59 historical Shigella strains that represent the four species and a variety of serotypes.

High-Quality Whole-Genome Sequences for 21 Enterotoxigenic Escherichia coli Strains Generated with PacBio Sequencing.

Enterotoxigenic Escherichia coli (ETEC) is an important diarrheagenic pathogen. We report here the high-quality whole-genome sequences of 21 ETEC strains isolated from patients in the United States, international diarrheal surveillance studies, and cruise ship outbreaks.

High-Quality Complete and Draft Genome Sequences for Three Escherichia spp. and Three Shigella spp. Generated with Pacific Biosciences and Illumina Sequencing and Optical Mapping.

Escherichia spp., including E. albertii and E. coli, Shigella dysenteriae, and S. flexneri are causative agents of foodborne disease. We report here reference-level whole-genome sequences of E. albertii (2014C-4356), E. coli (2011C-4315 and 2012C-4431), S. dysenteriae (BU53M1), and S. flexneri (94-3007 and 71-2783).

Microscale water distribution and its effects on organic carbon decomposition in unsaturated soils.

Microscale water distribution in the subsurface is key to many geochemical and biogeochemical reactions. This study investigated microscale water distribution and movement in unsaturated soils using micro-continuum hydrodynamic models, and examined the effect of microscale water distribution on organic carbon (C) decomposition using a micro-continuum biogeochemical reaction model. The micro-continuum hydrodynamic model that relates capillary pressure to porosity captured the measured water imbibition curve at the core scale, and exhibited reasonable water distribution and movement at the microscale. The simulations of organic C decomposition illustrate that microscale water distribution strongly affected the distribution of C decomposition rates by regulating the availability of dissolved organic C and oxygen. Particularly, changes in water distribution altered the location and intensity of reactive hotspots and thereby CO2 flux from soils. The microscale interactions between water content and organic C decomposition rate provide underlying mechanisms for explaining macroscale phenomenon observed in laboratory and fields. Overall, this study presents a useful tool for explicating hydro-biogeochemical behaviors in the subsurface by integrating micro-continuum hydrodynamic and biogeochemical reaction modeling.

Engineering the Autotroph Methanococcus maripaludis for Geraniol Production.

The rapid autotrophic growth of the methanogenic archaeon Methanococcus maripaludis on H2 and CO2 makes it an attractive microbial chassis to inexpensively produce biochemicals. To explore this potential, a synthetic gene encoding geraniol synthase (GES) derived from Ocimum basilicum was cloned into a M. maripaludis expression vector under selection for puromycin resistance. Recombinant expression of GES in M. maripaludis during autotrophic growth on H2/CO2 or formate yielded geraniol at 2.8 and 4.0 mg g(-1) of dry weight, respectively. The yield of geraniol decreased 2-3-fold when organic carbon sources were added to stimulate heterotrophic growth. In the absence of puromycin, geraniol production during autotrophic growth on formate increased to 4.6 mg g(-1) of dry weight. A conceptual model centered on the autotrophic acetyl coenzyme A biosynthetic pathway identified strategies to divert more autotrophic carbon flux to geraniol production.

A scalable method for molecular network reconstruction identifies properties of targets and mutations in acute myeloid leukemia.

A key aim of systems biology is the reconstruction of molecular networks. We do not yet, however, have networks that integrate information from all datasets available for a particular clinical condition. This is in part due to the limited scalability, in terms of required computational time and power, of existing algorithms. Network reconstruction methods should also be scalable in the sense of allowing scientists from different backgrounds to efficiently integrate additional data. We present a network model of acute myeloid leukemia (AML). In the current version (AML 2.1), we have used gene expression data (both microarray and RNA-seq) from 5 different studies comprising a total of 771 AML samples and a protein-protein interactions dataset. Our scalable network reconstruction method is in part based on the well-known property of gene expression correlation among interacting molecules. The difficulty of distinguishing between direct and indirect interactions is addressed by optimizing the coefficient of variation of gene expression, using a validated gold-standard dataset of direct interactions. Computational time is much reduced compared to other network reconstruction methods. A key feature is the study of the reproducibility of interactions found in independent clinical datasets. An analysis of the most significant clusters, and of the network properties (intraset efficiency, degree, betweenness centrality, and PageRank) of common AML mutations demonstrated the biological significance of the network. A statistical analysis of the response of blast cells from 11 AML patients to a library of kinase inhibitors provided an experimental validation of the network. A combination of network and experimental data identified CDK1, CDK2, CDK4, and CDK6 and other kinases as potential therapeutic targets in AML.