The PATRIC Bioinformatics Resource Center: expanding data and analysis capabilities

Isolation of Brucella inopinata from a White's tree frog (Litoria caerulea): pose exotic frogs a potential risk to human health?

Introduction

Cold-blooded hosts, particularly exotic frogs, have become a newly recognized reservoir for atypical Brucella species and strains worldwide, but their pathogenicity to humans remains largely unknown. Here we report the isolation and molecular characterization of a B. inopinata strain (FO700662) cultured from clinical samples taken from a captive diseased White's Tree Frog (Litoria caerulea) in Switzerland. The isolation of B. inopinata from a frog along with other reports of human infection by atypical Brucella raises the question of whether atypical Brucella could pose a risk to human health and deserves further attention.

Methods

The investigations included histopathological analysis of the frog, bacterial culture and in-depth molecular characterization of strain FO700662 based on genome sequencing data.

Results and Discussion

Originally identified as Ochrobactrum based on its rapid growth and biochemical profile, strain FO700622 was positive for the Brucella- specific markers bcsp31 and IS711. It showed the specific banding pattern of B. inopinata in conventional Bruce-ladder multiplex PCR and also had identical 16S rRNA and recA gene sequences as B. inopinata. Subsequent genome sequencing followed by core genome-based MLST (cgMLST) analysis using 2704 targets (74% of the total chromosome) revealed only 173 allelic differences compared to the type strain of B. inopinata BO1^T, while previously considered the closest related strain BO2 differed in 2046 alleles. The overall average nucleotide identity (ANI) between the type strain BO1^T and FO700622 was 99,89%, confirming that both strains were almost identical. In silico MLST-21 and MLVA-16 also identified strain FO700662 as B. inopinata. The nucleotide and amino acid-based phylogenetic reconstruction and comparative genome analysis again placed the isolate together with B. inopinata with 100% support. In conclusion, our data unequivocally classified strain FO700622, isolated from an exotic frog, as belonging to B. inopinata.

Bacteriophage therapy against pathological Klebsiella pneumoniae ameliorates the course of primary sclerosing cholangitis

Primary sclerosing cholangitis (PSC) is characterized by progressive biliary inflammation and fibrosis. Although gut commensals are associated with PSC, their causative roles and therapeutic strategies remain elusive. Here we detect abundant Klebsiella pneumoniae (Kp) and Enterococcus gallinarum in fecal samples from 45 PSC patients, regardless of intestinal complications. Carriers of both pathogens exhibit high disease activity and poor clinical outcomes. Colonization of PSC-derived Kp in specific pathogen-free (SPF) hepatobiliary injury-prone mice enhances hepatic Th17 cell responses and exacerbates liver injury through bacterial translocation to mesenteric lymph nodes. We developed a lytic phage cocktail that targets PSC-derived Kp with a sustained suppressive effect in vitro. Oral administration of the phage cocktail lowers Kp levels in Kp-colonized germ-free mice and SPF mice, without off-target dysbiosis. Furthermore, we demonstrate that oral and intravenous phage administration successfully suppresses Kp levels and attenuates liver inflammation and disease severity in hepatobiliary injury-prone SPF mice. These results collectively suggest that using a lytic phage cocktail shows promise for targeting Kp in PSC.

The basis of antigenic operon fragmentation in Bacteroidota and commensalism

The causes for variability of pro-inflammatory surface antigens that affect gut commensal/opportunistic dualism within the phylum Bacteroidota remain unclear (1, 2). Using the classical lipopolysaccharide/O-antigen 'rfb operon' in Enterobacteriaceae as a surface antigen model (5-gene-cluster rfbABCDX), and a recent rfbA-typing strategy for strain classification (3), we characterized the architecture/conservancy of the entire rfb operon in Bacteroidota. Analyzing complete genomes, we discovered that most Bacteroidota have the rfb operon fragmented into non-random gene-singlets and/or doublets/triplets, termed 'minioperons'. To reflect global operon integrity, duplication, and fragmentation principles, we propose a five-category (infra/supernumerary) cataloguing system and a Global Operon Profiling System for bacteria. Mechanistically, genomic sequence analyses revealed that operon fragmentation is driven by intra-operon insertions of predominantly Bacteroides-DNA (thetaiotaomicron/fragilis) and likely natural selection in specific micro-niches. Bacteroides-insertions, also detected in other antigenic operons (fimbriae), but not in operons deemed essential (ribosomal), could explain why Bacteroidota have fewer KEGG-pathways despite large genomes (4). DNA insertions overrepresenting DNA-exchange-avid species, impact functional metagenomics by inflating gene-based pathway inference and overestimating 'extra-species' abundance. Using bacteria from inflammatory gut-wall cavernous micro-tracts (CavFT) in Crohn's Disease (5), we illustrate that bacteria with supernumerary-fragmented operons cannot produce O-antigen, and that commensal/CavFT Bacteroidota stimulate macrophages with lower potency than Enterobacteriaceae, and do not induce peritonitis in mice. The impact of 'foreign-DNA' insertions on pro-inflammatory operons, metagenomics, and commensalism offers potential for novel diagnostics and therapeutics.

Pseudomonas quebecensis sp. nov., a bacterium isolated from root-zone soil of a native legume, Amphicarpaea bracteata (L.) Fernald, in Quebec, Canada

Four bacterial strains (S1Bt3, S1Bt7, S1Bt30 and S1Bt42^T) isolated from soil collected from the rhizosphere of a native legume, Amphicarpaea bracteata, were investigated using a polyphasic approach. Colonies were fluorescent, white-yellowish, circular and convex with regular margins on King's B medium. Cells were Gram-reaction-negative, aerobic, non-spore-forming rods. Oxidase- and catalase-positive. The optimal growth temperature of the strains was 37 °C. Phylogenetic analysis of the 16S rRNA gene sequences placed the strains within the genus Pseudomonas. Analysis of the 16S rRNA-rpoD-gyrB concatenated sequences clustered the strains and well separated from Pseudomonas rhodesiae CIP 104664^T and Pseudomonas grimontii CFM 97-514^T with the type strains of the closest species. Phylogenomic analysis of 92 up-to-date bacterial core gene and matrix-assisted laser desorption/ionization-time-of-flight MS biotyper data confirmed the distinct clustering pattern of these four strains. Digital DNA-DNA hybridization (41.7 %-31.2 %) and average nucleotide identity (91.1 %-87.0 %) values relative to closest validly published Pseudomonas species were below the species delineation thresholds of 70 and 96 %, respectively. Fatty acid composition results validated the taxonomic position of the novel strains in the genus Pseudomonas. Phenotypic characteristics from carbon utilization tests differentiated the novel strains from closely related Pseudomonas species. In silico prediction of secondary metabolite biosynthesis gene clusters in the whole-genome sequences of the four strains revealed the presence of 11 clusters involved in the production of siderophore, redox-cofactor, betalactone, terpene, arylpolyene and nonribosomal peptides. Based on phenotypic and genotypic data, strains S1Bt3, S1Bt7, S1Bt30 and S1Bt42^T represent a novel species for which the name Pseudomonas quebecensis sp. nov. is proposed. The type strain is S1Bt42^T (=DOAB 746^T=LMG 32141^T=CECT 30251^T). The genomic DNA G+C content is 60.95 mol%.

Reconstructor: a COBRApy compatible tool for automated genome-scale metabolic network reconstruction with parsimonious flux-based gap-filling

MOTIVATION

Genome-scale metabolic network reconstructions (GENREs) are valuable for understanding cellular metabolism in silico. Several tools exist for automatic GENRE generation. However, these tools frequently (i) do not readily integrate with some of the widely-used suites of packaged methods available for network analysis, (ii) lack effective network curation tools, (iii) are not sufficiently user-friendly, and (iv) often produce low-quality draft reconstructions.

RESULTS

Here, we present Reconstructor, a user-friendly, COBRApy-compatible tool that produces high-quality draft reconstructions with reaction and metabolite naming conventions that are consistent with the ModelSEED biochemistry database and includes a gap-filling technique based on the principles of parsimony. Reconstructor can generate SBML GENREs from three input types: annotated protein .fasta sequences (Type 1 input), a BLASTp output (Type 2), or an existing SBML GENRE that can be further gap-filled (Type 3). While Reconstructor can be used to create GENREs of any species, we demonstrate the utility of Reconstructor with bacterial reconstructions. We demonstrate how Reconstructor readily generates high-quality GENRES that capture strain, species, and higher taxonomic differences in functional metabolism of bacteria and are useful for further biological discovery.

AVAILABILITY AND IMPLEMENTATION

The Reconstructor Python package is freely available for download. Complete installation and usage instructions and benchmarking data are available at http://github.com/emmamglass/reconstructor.

Complete genome sequencing of Bacillus cabrialesii TE3T: A plant growth-promoting and biological control agent isolated from wheat (Triticum turgidum subsp. durum) in the Yaqui Valley

Bacillus cabrialesii TE3^T is a strictly aerobic and Gram-stain-positive plant growth-promoting bacterium, motile and catalase-positive. In addition, strain TE3^T was also recently described as a biological control agent. Here, we present the complete circularized genome of this type strain, as well as a whole genome analysis identifying genes of agricultural interest. Thus, a hybrid assembly method was performed using short-read sequencing through the Illumina MiSeq platform, and long-read sequencing through the MinION sequencing technology by Oxford Nanopore Technology (ONT). This assembly method showed a closed circular chromosome of 4,125,766 bp and 44.2% G + C content. The strain TE3^T genome annotation, based on the RAST platform, presented 4,282 Coding DNA sequences (CDS) distributed in 335 subsystems, from which 4 CDS are related to the promotion of plant growth and 28 CDS to biological control. Also, Prokka (Rapid Prokaryotic Genome Annotation) predicted a total of 119 RNAs composed of 87 tRNAs, 31 rRNA, and 1 tmRNA; and the PGAP (Prokaryotic Genome Annotation Pipeline) predicted a total of 4,212 genes (3,991 CDS). Additionally, seven putative biosynthetic gene clusters were identified by antiSMASH, such as Fengycin, Bacilysin, Subtilosin A, Bacillibactin, Bacillaene, Surfactin, and Rizocticin A, which are related to antimicrobial and antifungal properties, whose gene presence was further supported by the Prokaryotic Genome Annotation Pipeline (PGAP) annotation. Thus, the complete genome of Bacillus cabrialesii TE3^T showed promising bioactivities for the use of this type strain to bioformulate bacterial inoculants for sustainable agriculture.

Chromosome folding and prophage activation reveal specific genomic architecture for intestinal bacteria

BACKGROUND

Bacteria and their viruses, bacteriophages, are the most abundant entities of the gut microbiota, a complex community of microorganisms associated with human health and disease. In this ecosystem, the interactions between these two key components are still largely unknown. In particular, the impact of the gut environment on bacteria and their associated prophages is yet to be deciphered.

RESULTS

To gain insight into the activity of lysogenic bacteriophages within the context of their host genomes, we performed proximity ligation-based sequencing (Hi-C) in both in vitro and in vivo conditions on the 12 bacterial strains of the OMM¹² synthetic bacterial community stably associated within mice gut (gnotobiotic mouse line OMM¹²). High-resolution contact maps of the chromosome 3D organization of the bacterial genomes revealed a wide diversity of architectures, differences between environments, and an overall stability over time in the gut of mice. The DNA contacts pointed at 3D signatures of prophages leading to 16 of them being predicted as functional. We also identified circularization signals and observed different 3D patterns between in vitro and in vivo conditions. Concurrent virome analysis showed that 11 of these prophages produced viral particles and that OMM¹² mice do not carry other intestinal viruses.

CONCLUSIONS

The precise identification by Hi-C of functional and active prophages within bacterial communities will unlock the study of interactions between bacteriophages and bacteria across conditions (healthy vs disease). Video Abstract.

Shine: A novel strategy to extract specific, sensitive and well-conserved biomarkers from massive microbial genomic datasets

BACKGROUND

Concentrations of the pathogenic microorganisms' DNA in biological samples are typically low. Therefore, DNA diagnostics of common infections are costly, rarely accurate, and challenging. Limited by failing to cover updated epidemic testing samples, computational services are difficult to implement in clinical applications without complex customized settings. Furthermore, the combined biomarkers used to maintain high conservation may not be cost effective and could cause several experimental errors in many clinical settings. Given the limitations of recent developed technology, 16S rRNA is too conserved to distinguish closely related species, and mosaic plasmids are not effective as well because of their uneven distribution across prokaryotic taxa.

RESULTS

Here, we provide a computational strategy, Shine, that allows extraction of specific, sensitive and well-conserved biomarkers from massive microbial genomic datasets. Distinguished with simple concatenations with blast-based filtering, our method involves a de novo genome alignment-based pipeline to explore the original and specific repetitive biomarkers in the defined population. It can cover all members to detect newly discovered multicopy conserved species-specific or even subspecies-specific target probes and primer sets. The method has been successfully applied to a number of clinical projects and has the overwhelming advantages of automated detection of all pathogenic microorganisms without the limitations of genome annotation and incompletely assembled motifs. Using on our pipeline, users may select different configuration parameters depending on the purpose of the project for routine clinical detection practices on the website https://bioinfo.liferiver.com.cn with easy registration.

CONCLUSIONS

The proposed strategy is suitable for identifying shared phylogenetic markers while featuring low rates of false positive or false negative. This technology is suitable for the automatic design of minimal and efficient PCR primers and other types of detection probes.

Comparative genomics of Edwardsiella anguillarum and Edwardsiella piscicida isolated in Taiwan enables the identification of distinctive features and potential virulence factors using Oxford-Nanopore MinION® sequencing

Edwardsiella tarda (ET) and Edwardsiella anguillarum (EA) are the most harmful bacterial fish pathogens in Taiwan. However, there is confusion regarding the genotypic identification of E. tarda and E. piscicida (EP). Therefore, we used a novel Nanopore MinION MK1C platform to sequence and compare the complete genomes of E. piscicida and E. anguillarum. The number of coding genes, rRNA, and tRNA recorded for E. anguillarum and E. piscicida were 8322, 25, and 98, and 5458, 25, and 98, respectively. Ribosomal multilocus sequence typing (rMLST) for E. piscicida indicated 35 rps. The shared clusters between E. anguillarum and E. piscicida indicated several unique clusters for the individual genomes. The phylogenetic tree analysis for all complete genomes indicated that E. anguillarum and E. piscicida were placed into two species-specific genotypes. Distribution of subsystems for annotated genomes found that genes related to virulence, defence, and disease for E. anguillarum were 103 and those for E. piscicida were 60 and pathogenic islands (PI) were 498 and 225, respectively. Vaccine candidates were identified in silico from the core genes using high antigenic, solubility, and secretion probabilities. Altogether, the genome data revealed distinctive features between E. anguillarum and E. piscicida, which suggest different pathogenicity and thus the need for separate preventive strategies.

Mini-review: Current and Future Perspectives on Microbially Focused Restoration Strategies in Tallgrass Prairies

Ecosystem restoration is a critical conservation strategy, especially for increasing resilience and resistance to climate change. Current restoration efforts that convert reclaimed agricultural land to native tallgrass prairies typically focus on aboveground communities, but it can take decades to restore soil microbial biodiversity and function using these strategies, if they recover at all. This incomplete restoration can have detrimental impacts on longer-term restoration goals, such as supporting late-successional plant species and facilitating soil carbon sequestration. Soil microorganisms are key components in determining the fate of organic material that enters the soil. They mediate decomposition rates and contribute to plant-microbe-soil interactions, produce microbial biomass, necromass, and metabolic products, and physically protect soil carbon through aggregation. Interactions with plants and controls over soil carbon vary widely depending on the specific microbial taxa present, their physiology, their functional capabilities, and their responses to environmental stressors. Thus, the ability for new restorations, prairie conservation corridors, and prairies planted in marginal lands to act as carbon sinks and help balance greenhouse gas emissions can depend on the success of microbial restoration. Next-generation sequencing approaches can support novel methods for evaluating existing restoration practices and developing microbially focused management strategies. This review summarizes the growing body of literature describing microbially focused tallgrass prairie restoration and considers when and how integrating next-generation sequencing approaches into management efforts can be beneficial. We provide a roadmap for future restoration efforts where microbial ecologists, restoration ecologists, and land managers can work together to meet their goals to promote climate-ready restored ecosystems.

Composition and function of the Galapagos penguin gut microbiome vary with age, location, and a putative bacterial pathogen

Microbial colonization plays a direct role in host health. Understanding the ecology of the resident microbial community for a given host species is thus an important step for detecting population vulnerabilities like disease. However, the idea of integrating microbiome research into conservation is still relatively new, and wild birds have received less attention in this field than mammals or domesticated animals. Here we examine the composition and function of the gut microbiome of the endangered Galapagos penguin (Spheniscus mendiculus) with the goals of characterizing the normal microbial community and resistome, identifying likely pathogens, and testing hypotheses of structuring forces for this community based on demographics, location, and infection status. We collected fecal samples from wild penguins in 2018 and performed 16S rRNA gene sequencing and whole genome sequencing (WGS) on extracted DNA. 16S sequencing revealed that the bacterial phyla Fusobacteria, Epsilonbacteraeota, Firmicutes, and Proteobacteria dominate the community. Functional pathways were computed from WGS data, showing genetic functional potential primarily focused on metabolism-amino acid metabolism, carbohydrate metabolism, and energy metabolism are the most well-represented functional groups. WGS samples were each screened for antimicrobial resistance, characterizing a resistome made up of nine antibiotic resistance genes. Samples were screened for potential enteric pathogens using virulence factors as indicators; Clostridium perfringens was revealed as a likely pathogen. Overall, three factors appear to be shaping the alpha and beta diversity of the microbial community: penguin developmental stage, sampling location, and C. perfringens. We found that juvenile penguins have significantly lower alpha diversity than adults based on three metrics, as well as significantly different beta diversity. Location effects are minimal, but one site has significantly lower Shannon diversity than the other primary sites. Finally, when samples were grouped by C. perfringens virulence factors, we found dramatic changes in beta diversity based on operational taxonomic units, protein families, and functional pathways. This study provides a baseline microbiome for an endangered species, implicates both penguin age and the presence of a potential bacterial pathogen as primary factors associated with microbial community variance, and reveals widespread antibiotic resistance genes across the population.

Beyond the spore, the exosporium sugar anthrose impacts vegetative Bacillus anthracis gene regulation in cis and trans

The Bacillus anthracis exosporium nap is the outermost portion of spore that interacts with the environment and host systems. Changes to this layer have the potential to impact wide-ranging physiological and immunological processes. The unique sugar, anthrose, normally coats the exosporium nap at its most distal points. We previously identified additional mechanisms rendering B. anthracis anthrose negative. In this work, several new ant ^- B. anthracis strains are identified and the impact of anthrose negativity on spore physiology is investigated. We demonstrate that live-attenuated Sterne vaccines as well as culture filtrate anthrax vaccines generate antibodies targeting non-protein components of the spore. The role of anthrose as a vegetative B. anthracis Sterne signaling molecule is implicated by luminescent expression strain assays, RNA-seq experiments, and toxin secretion analysis by western blot. Pure anthrose and the sporulation-inducing nucleoside analogue decoyinine had similar effects on toxin expression. Co-culture experiments demonstrated gene expression changes in B. anthracis depend on intracellular anthrose status (cis) in addition to anthrose status of extracellular interactions (trans). These findings provide a mechanism for how a unique spore-specific sugar residue affects physiology, expression and genetics of vegetative B. anthracis with impacts on the ecology, pathogenesis, and vaccinology of anthrax.

The genome of Bacillus tequilensis EA-CB0015 sheds light into its epiphytic lifestyle and potential as a biocontrol agent

Different Bacillus species have successfully been used as biopesticides against a broad range of plant pathogens. Among these, Bacillus tequilensis EA-CB0015 has shown to efficiently control Black sigatoka disease in banana plants, presumably by mechanisms of adaptation that involve modifying the phyllosphere environment. Here, we report the complete genome of strain EA-CB0015, its precise taxonomic identity, and determined key genetic features that may contribute to its effective biocontrol of plant pathogens. We found that B. tequilensis EA-CB0015 harbors a singular 4 Mb circular chromosome, with 3,951 protein-coding sequences. Multi-locus sequence analysis (MLSA) and average nucleotide identity (ANI) analysis classified strain EA-CB0015 as B. tequilensis. Encoded within its genome are biosynthetic gene clusters (BGCs) for surfactin, iturin, plipastatin, bacillibactin, bacilysin, subtilosin A, sporulation killing factor, and other natural products that may facilitate inter-microbial warfare. Genes for indole-acetic acid (IAA) synthesis, the use of diverse carbon sources, and a multicellular lifestyle involving motility, biofilm formation, quorum sensing, competence, and sporulation suggest EA-CB0015 is adept at colonizing plant surfaces. Defensive mechanisms to survive invading viral infections and preserve genome integrity include putative type I and type II restriction modification (RM) and toxin/antitoxin (TA) systems. The presence of bacteriophage sequences, genomic islands, transposable elements, virulence factors, and antibiotic resistance genes indicate prior occurrences of genetic exchange. Altogether, the genome of EA-CB0015 supports its function as a biocontrol agent against phytopathogens and suggest it has adapted to thrive within phyllosphere environments.

Draft Genome Sequence of Afifella sp. Strain JA880, Isolated from a Salt Pond

We report the draft genome sequence of Afifella sp. strain JA880, which was isolated from a saltwater pond near Pata, Gujarat, India. The genome assembly contains 3,794,364 bp, with a GC content of 63.5%. The genome sequence provides insights into the metabolic potential of Afifella sp. strain JA880.

Blastopirellula sediminis sp. nov. a new member of Pirellulaceae isolated from the Andaman and Nicobar Islands

Two cream-coloured strains (JC732^T, JC733) of Gram-stain negative, mesophilic, catalase and oxidase positive, aerobic bacteria which divide by budding, form crateriform structures, and cell aggregates were isolated from marine habitats of Andaman and Nicobar Islands, India. Both strains had genome size of 7.1 Mb and G + C content of 58.9%. Both strains showed highest 16S rRNA gene-based similarity with Blastopirellula retiformator Enr8^T (98.7%). Strains JC732^T and JC733 shared 100% identity of 16S rRNA gene and genome sequences. The coherence of both strains with the genus Blastopirellula was supported by the 16S rRNA gene based and the phylogenomic trees. Further, the chemo-taxonomic characters and the genome relatedness indices [ANI (82.4%), AAI (80.4%) and dDDH (25.2%)] also support the delineation at the species level. Both strains have the capability to degrade chitin and genome analysis shows the ability to fix N₂. Based on the phylogenetic, phylogenomic, comparative genomic, morphological, physiological, and biochemical characteristics, strain JC732^T is described as a new species of the genus Blastopirellula for which the name Blastopirellula sediminis sp. nov. is proposed, with strain JC733 as an additional strain.

Revisiting the taxonomy of the genus Rhodopirellula with the proposal for reclassification of the genus to Rhodopirellula sensu stricto, Aporhodopirellula gen. nov., Allorhodopirellula gen. nov. and Neorhodopirellula gen. nov

The current genus Rhodopirellula consists of marine bacteria which belong to the family Pirellulaceae of the phylum Planctomycetota. Members of the genus Rhodopirellula are aerobic, mesophiles and chemoheterotrophs. The here conducted analysis built on 16S rRNA gene sequence and multi-locus sequence analysis based phylogenomic trees suggested that the genus is subdivided into four clades. Existing Rhodopirellula species were studied extensively based on phenotypic, genomic and chemotaxonomic parameters. The heterogeneity was further confirmed by overall genome-related indices (OGRI) including digital DNA-DNA hybridization (dDDH), average nucleotide identity (ANI), average amino acid identity (AAI), and percentage of conserved proteins (POCP). AAI and POCP values between the clades of the genus Rhodopirellula were 62.2-69.6% and 49.5-62.5%, respectively. Comparative genomic approaches like pan-genome analysis and conserved signature indels (CSIs) also support the division of the clades. The genomic incoherence of the members of the genus is further supported by variations in phenotypic characteristics. Thus, with the here applied integrated comparative genomic and polyphasic approaches, we propose the reclassification of the genus Rhodopirellula to three new genera: Aporhodopirellula gen. nov., Allorhodopirellula gen. nov., and Neorhodopirellula gen. nov.

Contrasting Genetic Diversity of Listeria Pathogenicity Islands 3 and 4 Harbored by Nonpathogenic Listeria spp

Listeria monocytogenes causes the severe foodborne disease listeriosis. Several clonal groups of L. monocytogenes possess the pathogenicity islands Listeria pathogenicity island 3 (LIPI-3) and LIPI-4. Here, we investigated the prevalence and genetic diversity of LIPI-3 and LIPI-4 among 63 strains of seven nonpathogenic Listeria spp. from the natural environment, i.e., wildlife (black bears [Ursus americanus]) and surface water. Analysis of the whole-genome sequence data suggested that both islands were horizontally acquired but differed considerably in their incidence and genetic diversity. LIPI-3 was identified among half of the L. innocua strains in the same genomic location as in L. monocytogenes (guaA hot spot) in a truncated form, with only three strains harboring full-length LIPI-3, and a highly divergent partial LIPI-3 was observed in three Listeria seeligeri strains, outside the guaA hot spot. Premature stop codons (PMSCs) and frameshifts were frequently noted in the LIPI-3 gene encoding listeriolysin S. On the other hand, full-length LIPI-4 without any PMSCs was found in all Listeria innocua strains, in the same genomic location as L. monocytogenes and with ~85% similarity to the L. monocytogenes counterpart. Our study provides intriguing examples of genetic changes that pathogenicity islands may undergo in nonpathogenic bacterial species, potentially in response to environmental pressures that promote either maintenance or degeneration of the islands. Investigations of the roles that LIPI-3 and LIPI-4 play in nonpathogenic Listeria spp. are warranted to further understand the differential evolution of genetic elements in pathogenic versus nonpathogenic hosts of the same genus. IMPORTANCE Listeria monocytogenes is a serious foodborne pathogen that can harbor the pathogenicity islands Listeria pathogenicity island 3 (LIPI-3) and LIPI-4. Intriguingly, these have also been reported in nonpathogenic L. innocua from food and farm environments, though limited information is available for strains from the natural environment. Here, we analyzed whole-genome sequence data of nonpathogenic Listeria spp. from wildlife and surface water to further elucidate the genetic diversity and evolution of LIPI-3 and LIPI-4 in Listeria. While the full-length islands were found only in L. innocua, LIPI-3 was uncommon and exhibited frequent truncation and genetic diversification, while LIPI-4 was remarkable in being ubiquitous, albeit diversified from L. monocytogenes. These contrasting features demonstrate that pathogenicity islands in nonpathogenic hosts can evolve along different trajectories, leading to either degeneration or maintenance, and highlight the need to examine their physiological roles in nonpathogenic hosts.

Yersiniomics, a Multi-Omics Interactive Database for Yersinia Species

The genus Yersinia includes a large variety of nonpathogenic and life-threatening pathogenic bacteria, which cause a broad spectrum of diseases in humans and animals, such as plague, enteritis, Far East scarlet-like fever (FESLF), and enteric redmouth disease. Like most clinically relevant microorganisms, Yersinia spp. are currently subjected to intense multi-omics investigations whose numbers have increased extensively in recent years, generating massive amounts of data useful for diagnostic and therapeutic developments. The lack of a simple and centralized way to exploit these data led us to design Yersiniomics, a web-based platform allowing straightforward analysis of Yersinia omics data. Yersiniomics contains a curated multi-omics database at its core, gathering 200 genomic, 317 transcriptomic, and 62 proteomic data sets for Yersinia species. It integrates genomic, transcriptomic, and proteomic browsers, a genome viewer, and a heatmap viewer to navigate within genomes and experimental conditions. For streamlined access to structural and functional properties, it directly links each gene to GenBank, the Kyoto Encyclopedia of Genes and Genomes (KEGG), UniProt, InterPro, IntAct, and the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) and each experiment to Gene Expression Omnibus (GEO), the European Nucleotide Archive (ENA), or the Proteomics Identifications Database (PRIDE). Yersiniomics provides a powerful tool for microbiologists to assist with investigations ranging from specific gene studies to systems biology studies. IMPORTANCE The expanding genus Yersinia is composed of multiple nonpathogenic species and a few pathogenic species, including the deadly etiologic agent of plague, Yersinia pestis. In 2 decades, the number of genomic, transcriptomic, and proteomic studies on Yersinia grew massively, delivering a wealth of data. We developed Yersiniomics, an interactive web-based platform, to centralize and analyze omics data sets on Yersinia species. The platform allows user-friendly navigation between genomic data, expression data, and experimental conditions. Yersiniomics will be a valuable tool to microbiologists.

Exploring the Potential Molecular Mechanisms of Interactions between a Probiotic Consortium and Its Coral Host

Beneficial microorganisms for corals (BMCs) have been demonstrated to be effective probiotics to alleviate bleaching and mitigate coral mortality in vivo. The selection of putative BMCs is traditionally performed manually, using an array of biochemical and molecular tests for putative BMC traits. We present a comprehensive genetic survey of BMC traits using a genome-based framework for the identification of alternative mechanisms that can be used for future in silico selection of BMC strains. We identify exclusive BMC traits associated with specific strains and propose new BMC mechanisms, such as the synthesis of glycine betaine and ectoines. Our roadmap facilitates the selection of BMC strains while increasing the array of genetic targets that can be included in the selection of putative BMC strains to be tested as coral probiotics. IMPORTANCE Probiotics are currently the main hope as a potential medicine for corals, organisms that are considered the marine "canaries of the coal mine" and that are threatened with extinction. Our experiments have proved the concept that probiotics mitigate coral bleaching and can also prevent coral mortality. Here, we present a comprehensive genetic survey of probiotic traits using a genome-based framework. The main outcomes are a roadmap that facilitates the selection of coral probiotic strains while increasing the array of mechanisms that can be included in the selection of coral probiotics.

First insights into the gut microbiome of Diatraea saccharalis: From a sugarcane pest to a reservoir of new bacteria with biotechnological potential

A country’s biodiversity is a key resource for the development of a sustainable bioeconomy. However, often the most biodiverse countries on the planet hardly profit from their biological diversity. On the contrary, occasionally components of that biodiversity become a threat to society and its food sustainability. That is the case of the sugarcane borer Diatraea saccharalis. Here, the analysis of the bacteria associated with the digestive tract of D. saccharalis reveals a rich and diverse microbiota. Two types of diets were analyzed under laboratory conditions. The metataxonomic analysis revealed a number of taxa common to most of the larval pools analyzed with relative abundances exceeding 5%, and five families of bacteria which have also been reported in the gut of another Lepidoptera. A large fraction of microorganisms detected by amplicon sequencing were considered to be rare and difficult to cultivate. However, among the cultivable microorganisms, 12 strains with relevant biotechnological features were identified. The strain that showed the highest cellulolytic activity (GCEP-101) was genome sequenced. The analysis of the GCEP-101 complete genome revealed that the values of 16S rRNA identity, the Average Nucleotide Identity, and the digital DNA–DNA hybridization place the strain as a candidate for a new species within the genus Pseudomonas. Moreover, the genome annotation of the putative new species evidenced the presence of genes associated with cellulose degradation, revealing the hidden potential of the pest as a reservoir of biotechnologically relevant microorganisms.

Recombination Drives Evolution of Carbapenem-Resistant Klebsiella pneumoniae Sequence Type 11 KL47 to KL64 in China

Carbapenem-resistant Enterobacteriaceae, especially carbapenemase-producing Klebsiella pneumoniae, is an urgent problem in health care facilities worldwide. K. pneumoniae isolates classified as sequence type 11 (ST11) are largely responsible for the spread of carbapenem-resistant K. pneumoniae (CRKP) in China. Our previous phylogenetic reconstruction suggested that CRKP ST11 capsular locus 64 (KL64) was derived from an ST11-KL47 ancestor through recombination. However, the molecular origin of KL64 remains largely unknown, and our understanding of the recombination is incomplete. Here, we screened a global sample of 22,600 K. pneumoniae genomes and searched for KL64-harboring STs, which were found to be ST1764, ST3685, ST1764-1LV, ST30, ST505, ST147, and ST11, wherein ST1764, ST3685, ST1764-1LV, and ST30 belonged to a clonal complex, CC1764. We compared the genetic structures of representative strains from ST11-KL47, ST11-KL64, CC1764-KL64, ST505-KL64, and ST147-KL64 and further performed phylogenetic analysis and single-nucleotide polymorphism analysis among 248 isolates from all these STs. The results suggested a recombination event has occurred in a homologous ~154-kb region covering KL and the lipopolysaccharide biosynthesis locus (OL) between a recipient ST11-KL47-OL101 and a donor CC1764 (except ST30), giving rise to ST11-KL64-O2v1 strains (recombination I). Furthermore, we also found an infrequent ST11-KL64-O2v1 subclone which was not produced by recombination I but was hybridized from ST11-KL47-OL101 and ST147-KL64-O2v1 strains through recombination of a homologous ~485-kb region covering KL and OL (recombination II). Our findings provide important insights into the role of recombination in the evolution of clinical strains and the diversity of capsule and lipopolysaccharide of widely distributed KPC-associated ST11 K. pneumoniae in China. IMPORTANCE Chromosomal recombination events are considered to contribute to the genetic diversification and ultimate success of many bacterial pathogens. A previous study unravelled the molecular evolution history of ST258 strains, which have been largely responsible for the spread of KPC in the United States. Here, we used comparative genomic analyses to discover two recombination events in ST11 CRKP strains, which is a predominant KPC-associated CRKP clone in China. Two new ST11-CRKP subclones with altered capsule and lipopolysaccharide, which are two primary determinants of antigenicity and antigenic diversity among K. pneumoniae, were produced through these two recombination events, respectively. Horizontal transfer of the KL and OL appears to be a crucial element driving the molecular evolution of K. pneumoniae strains. These findings not only extend our understanding of the molecular evolutionary history of ST11 but also are an important step toward the development of preventive, diagnostic, and therapeutic strategies for CRKP.

Soil Collected from a Single Great Smoky Mountains Trail Contains a Diversity of Listeria monocytogenes and Listeria spp

Listeria monocytogenes, a foodborne pathogen, and other Listeria spp. are present in natural environments. Isolating and characterizing strains from natural reservoirs can provide insight into the prevalence and diversity of Listeria spp. in these environments, elucidate their contribution to contamination of agricultural and food processing environments and food products, and lead to the discovery of novel species. In this study, we evaluated the diversity of Listeria spp. isolated from soil in a small region of the Great Smoky Mountains National Park, the most biodiverse national park in the U.S. National Park system. Of the 17 Listeria isolates recovered, whole-genome sequencing revealed that 14 were distinct strains. The strains represented a diversity of Listeria species (L. monocytogenes [n = 9], L. cossartiae subsp. cossartiae [n = 1], L. marthii [n = 1], L. booriae [n = 1], and a potentially novel Listeria sp. [n = 2]), as well as a diversity of sequence types based on multilocus sequence typing (MLST) and core genome MLST, including many novel designations. The isolates were not closely related (≥99.99% average nucleotide identity) to any isolates in public databases (NCBI, PATRIC), which also indicated novelty. The Listeria samples isolated in this study were collected from high-elevation sites near a creek that ultimately leads to the Mississippi River; thus, Listeria present in this natural environment could potentially travel downstream to a large region that includes portions of nine southeastern and midwestern U.S. states. This study provides insight into the diversity of Listeria spp. in the Great Smoky Mountains and indicates that this environment is a reservoir of novel Listeria spp. IMPORTANCE Listeria monocytogenes is a foodborne pathogen that can cause serious systemic illness that, although rare, usually results in hospitalization and has a relatively high mortality rate compared to other foodborne pathogens. Identification of novel and diverse Listeria spp. can provide insights into the genomic evolution, ecology, and evolution and variance of pathogenicity of this genus, especially in natural environments. Comparing L. monocytogenes and Listeria spp. isolates from natural environments, such as those recovered in this study, to contamination and/or outbreak strains may provide more information about the original natural sources of these strains and the pathways and mechanisms that lead to contamination of food products and agricultural or food processing environments.

Horizontal Gene Transfer and Loss of Serotype-Specific Genes in Listeria monocytogenes Can Lead to Incorrect Serotype Designations with a Commonly-Employed Molecular Serotyping Scheme

Listeria monocytogenes is a Gram-positive, facultative intracellular foodborne pathogen capable of causing severe, invasive illness (listeriosis). Three serotypes, 1/2a, 1/2b, and 4b, are leading contributors to human listeriosis, with 4b including the major hypervirulent clones. The multiplex PCR scheme developed by Doumith and collaborators employs primers targeting specific lineages (e.g., lineage II-specific lmo0737, lineage I-specific LMOf2365_2059) or serotypes (e.g., serotype 4b-specific LMOf2365_1900). The Doumith scheme (DS) is extensively employed for molecular serotyping of L. monocytogenes due to its high accuracy, relative ease, and affordability. However, for certain strains, the DS serotype designations are in conflict with those relying on antibody-based schemes or whole-genome sequence (WGS) analysis. In the current study, all 27 tested serotype 4b strains with sequence type 782 (ST782) within the hypervirulent clonal complex 2 (CC2) were designated 1/2b/3b using the DS. These strains lacked the serotype 4b-specific gene LMOf2365_1900, while retaining LMOf2365_2059, which, together with prs, yields the DS 1/2b/3b profile. Furthermore, 15 serotype 1/2a strains of four STs, mostly from water, were designated 1/2b/3b using the DS. These strains lacked the lmo0737 cassette but harbored genomic islands with LMOf2365_2059, thus yielding the DS 1/2b/3b profile. Lastly, we investigated a novel, dual 1/2a-1/2b profile obtained using the DS with 21 serotype 1/2a strains of four STs harboring both the lmo0737 cassette and genomic islands with LMOf2365_2059. The findings suggest that for certain strains and clones of L. monocytogenes the DS designations should be viewed with caution and complemented with alternative tools, e.g., traditional serotyping or WGS analysis. IMPORTANCE Listeria monocytogenes is a foodborne pathogen responsible for severe illness (listeriosis), especially in pregnant women and their fetuses, immunocompromised individuals, and the elderly. Three serotypes, 1/2a, 1/2b, and 4b, account for most human listeriosis, with certain serotype 4b clonal complexes (CCs) overrepresented in human disease. Serotyping remains extensively employed in Listeria epidemiologic investigations, and a multiplex PCR-based serotyping scheme is widely used. However, the PCR gene targets can be lost or gained via horizontal gene transfer, leading to novel PCR profiles without known serotype designations or to incorrect serotype assignments. Thus, an entire serotype 4b clone of the hypervirulent CC2 would be misidentified as serotype 1/2b, and several strains of serotype 1/2a would be identified as serotype 1/2b. Such challenges are especially common in novel clones from underexplored habitats, e.g., wildlife and surface water. The findings suggest caution in application of molecular serotyping, while highlighting Listeria's diversity and potential for horizontal gene transfer.

Genomic Epidemiological Analysis of Antimicrobial-Resistant Bacteria with Nanopore Sequencing

Antimicrobial-resistant (AMR) bacterial infections caused by clinically important bacteria, including ESKAPE pathogens (Enterococcus faecium, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa, and Enterobacter species) and mycobacteria (Mycobacterium tuberculosis and nontuberculous mycobacteria), have become a global public health threat. Their epidemic and pandemic clones often accumulate useful accessory genes in their genomes, such as AMR genes (ARGs) and virulence factor genes (VFGs). This process is facilitated by horizontal gene transfer among microbial communities via mobile genetic elements (MGEs), such as plasmids and phages. Nanopore long-read sequencing allows easy and inexpensive analysis of complex bacterial genome structures, although some aspects of sequencing data calculation and genome analysis methods are not systematically understood. Here we describe the latest and most recommended experimental and bioinformatics methods available for the construction of complete bacterial genomes from nanopore sequencing data and the detection and classification of genotypes of bacterial chromosomes, ARGs, VFGs, plasmids, and other MGEs based on their genomic sequences for genomic epidemiological analysis of AMR bacteria.

Genomic analysis of the international high-risk clonal lineage Klebsiella pneumoniae sequence type 395

BACKGROUND

Klebsiella pneumoniae, which is frequently associated with hospital- and community-acquired infections, contains multidrug-resistant (MDR), hypervirulent (hv), non-MDR/non-hv as well as convergent representatives. It is known that mostly international high-risk clonal lineages including sequence types (ST) 11, 147, 258, and 307 drive their global spread. ST395, which was first reported in the context of a carbapenemase-associated outbreak in France in 2010, is a less well-characterized, yet emerging clonal lineage.

METHODS

We computationally analyzed a large collection of K. pneumoniae ST395 genomes (n = 297) both sequenced in this study and reported previously. By applying multiple bioinformatics tools, we investigated the core-genome phylogeny and evolution of ST395 as well as distribution of accessory genome elements associated with antibiotic resistance and virulence features.

RESULTS

Clustering of the core-SNP alignment revealed four major clades with eight smaller subclades. The subclades likely evolved through large chromosomal recombination, which involved different K. pneumoniae donors and affected, inter alia, capsule and lipopolysaccharide antigen biosynthesis regions. Most genomes contained acquired resistance genes to extended-spectrum cephalosporins, carbapenems, and other antibiotic classes carried by multiple plasmid types, and many were positive for hypervirulence markers, including the siderophore aerobactin. The detection of "hybrid" resistance and virulence plasmids suggests the occurrence of the convergent ST395 pathotype.

CONCLUSIONS

To the best of our knowledge, this is the first study that investigated a large international collection of K. pneumoniae ST395 genomes and elucidated phylogenetics and detailed genomic characteristics of this emerging high-risk clonal lineage.

Solicola gregarius gen. nov., sp. nov., a soil actinobacterium isolated after enhanced cultivation with Micrococcus luteus culture supernatant

An actinobacterial strain, designated A5X3R13^T, was isolated from a compost soil suspension supplemented with extracellular material from a Micrococcus luteus-culture supernatant. The strain was cultured on tenfold-diluted reasoner's 2A agar. The cells were ovoid-to-rod shaped, non-motile, Gram-stain-positive, oxidase-negative, catalase-positive and had a width of 0.5 µm and a length of 0.8-1.2 µm. The results of both 16S rRNA-based phylogenetic and whole-genome analyses indicate that A5X3R13^T forms a distinct lineage within the family Nocardioidaceae (order Propionibacteriales). On the basis of the 16S rRNA gene sequence, A5X3R13^T was closely related to Aeromicrobium terrae CC-CFT486^T (96.2 %), Nocardioides iriomotensis IR27-S3^T (96.2 %), Nocardioides guangzhouensis 130^T (95.6 %), Marmoricola caldifontis YIM 730233^T (95.5 %), Aeromicrobium alkaliterrae KSL-107^T (95.4 %), Aeromicrobium choanae 9H-4^T (95.4 %), Aeromicrobium panaciterrae Gsoil 161^T (95.3 %), and Nocardioides jensenii NBRC 14755^T (95.2 %). The genome had a length of 4 915 757 bp, and its DNA G+C content was 68.5 mol %. The main fatty acids were 10-methyl C_17 : 0, C_16 : 0, C_15 : 0, C_18 : 0, C_17 : 0 and iso-C_16 : 0. The main polar lipids were phosphatidylglycerol, diphosphatidylglycerol, phosphatidylinositol and two unidentified phospholipids. MK-9(H₄) was the predominant respiratory quinone. The peptidoglycan type was A3γ (A41.1) and contained alanine, glycine, glutamic acid and ll-diaminopimelic acid in a molar ratio of 1.2 : 0.9 : 1.0 : 0.8. On the basis of the results of the phylogenetic and phenotypic analyses and comparisons with other members of the family Nocardioidaceae, strain A5X3R13^T is proposed to represent a novel species within a novel genus, for which the name Solicola gregarius gen. nov., sp. nov. is proposed. The type strain is A5X3R13^T (=DSM 112953^T=NCCB 100840^T).

Streptomyces guryensis sp. nov. exhibiting antimicrobial activity, isolated from riverside soil

A novel Gram-stain-positive, aerobic actinobacterial strain designated NR30^T was isolated from riverside soil. A polyphasic approach was employed for the taxonomic characterization of NR30^T. The strain developed extensively branched light brown to light pink substrate mycelia and light grey aerial mycelia, and produced spiny spores in loose spiral spore chains on ISP 3 and 4 agars. NR30^T grew at 10-40°C (optimum, 30°C), at pH 6.0-9.0 (optimum, pH 8.0) and in the presence of 0-3 % NaCl (optimum, 0 %). Analysis of 16S rRNA gene sequences indicated that NR30^T represents a member of the genus Streptomyces. NR30^T shared the highest 16S rRNA gene sequence similarity with Streptomyces cyaneus NRRL B-2296^T (98.6 %). On the basis of orthologous average nucleotide identity, NR30^T was most closely related to Streptomyces panaciradicis NBRC 109811^T with 86.3 % identity. The results of the digital DNA-DNA hybridization analysis also indicated low levels of relatedness with other species, as the highest value was observed with Streptomyces panaciradicis NBRC 109811^T (31.1 %). The major fatty acids of the strain were anteiso-C_15 : 0, C_16 : 0, iso-C_16 : 0 and anteiso-C_17 : 0. The major respiratory quinones were MK-9(H₈) and MK-9(H₆). The diagnostic polar lipids were diphosphatidylglycerol, phosphatidylethanolamine and phosphatidylinositol mannoside. The major cell wall diamino acid was ll-diaminopimelic acid, and the characteristic whole-cell sugars were rhamnose, ribose and glucose. The DNA G+C content was 70.3 mol %. NR30^T exhibited antimicrobial activity against several Gram-negative bacteria and yeasts. On the basis of the results of both phenotypic and phylogenetic analyses, strain NR30^T evidently represents a novel species of the genus Streptomyces, and the name Streptomyces guryensis sp. nov. (type strain=NR30^T =KCTC 49653^T=LMG 32476^T) is proposed accordingly.

Phylogenomic analysis of the genus Alcanivorax: proposal for division of this genus into the emended genus Alcanivorax and two novel genera Alloalcanivorax gen. nov. and Isoalcanivorax gen. nov

The members of the genus Alcanivorax are key players in the removal of petroleum hydrocarbons from polluted marine environments. More than half of the species were described in the last decade using 16S rRNA gene phylogeny and genomic-based metrics. However, the 16S rRNA gene identity (<94 %) between some members of the genus Alcanivorax suggested their imprecise taxonomic status. In this study, we examined the taxonomic positions of Alcanivorax species using 16S rRNA phylogeny and further validated them using phylogenomic-related indexes such as digital DNA-DNA hybridization (dDDH), average nucleotide identity (ANI), average amino acid identity (AAI), percentage of conserved proteins (POCP) and comparative genomic studies. ANI and dDDH values confirmed that all the Alcanivorax species were well described at the species level. The phylotaxogenomic analysis showed that Alcanivorax species formed three clades. The inter-clade values of AAI and POCP were less than 70 %. The pan-genome evaluation depicted that the members shared 1223 core genes and its number increased drastically when analysed clade-wise. Therefore, these results necessitate the transfer of clade II and clade III members into Isoalcanivorax gen. nov. and Alloalcanivorax gen. nov., respectively, along with the emended description of the genus Alcanivorax sensu stricto.

Pyridoxal 5'-phosphate synthesis and salvage in Bacteria and Archaea: predicting pathway variant distributions and holes

Pyridoxal 5’-phosphate or PLP is a cofactor derived from B₆ vitamers and essential for growth in all known organisms. PLP synthesis and salvage pathways are well characterized in a few model species even though key components, such as the vitamin B₆ transporters, are still to be identified in many organisms including the model bacteria Escherichia coli or Bacillus subtilis. Using a comparative genomic approach, PLP synthesis and salvage pathways were predicted in 5840 bacterial and archaeal species with complete genomes. The distribution of the two known de novo biosynthesis pathways and previously identified cases of non-orthologous displacements were surveyed in the process. This analysis revealed that several PLP de novo pathway genes remain to be identified in many organisms, either because sequence similarity alone cannot be used to discriminate among several homologous candidates or due to non-orthologous displacements. Candidates for some of these pathway holes were identified using published TnSeq data, but many remain. We find that ~10 % of the analysed organisms rely on salvage but further analyses will be required to identify potential transporters. This work is a starting point to model the exchanges of B₆ vitamers in communities, predict the sensitivity of a given organism to drugs targeting PLP synthesis enzymes, and identify numerous gaps in knowledge that will need to be tackled in the years to come.

PorinPredict: In Silico Identification of OprD Loss from WGS Data for Improved Genotype-Phenotype Predictions of P. aeruginosa Carbapenem Resistance

The increasing integration of genomics into routine clinical diagnostics requires reliable computational tools to identify determinants of antimicrobial resistance (AMR) from whole-genome sequencing data. Here, we developed PorinPredict, a bioinformatic tool that predicts defects of the Pseudomonas aeruginosa outer membrane porin OprD, which are strongly associated with reduced carbapenem susceptibility. PorinPredict relies on a database of intact OprD variants and reports inactivating mutations in the coding or promoter region. PorinPredict was validated against 987 carbapenemase-negative P. aeruginosa genomes, of which OprD loss was predicted for 454 out of 522 (87.0%) meropenem-nonsusceptible and 46 out of 465 (9.9%) meropenem-susceptible isolates. OprD loss was also found to be common among carbapenemase-producing isolates, resulting in even further increased MICs. Chromosomal mutations in quinolone resistance-determining regions and OprD loss commonly co-occurred, likely reflecting the restricted use of carbapenems for multidrug-resistant infections as recommended in antimicrobial stewardship programs. In combination with available AMR gene detection tools, PorinPredict provides a robust and standardized approach to link P. aeruginosa phenotypes to genotypes. IMPORTANCE Pseudomonas aeruginosa is a major cause of multidrug-resistant nosocomial infections. The emergence and spread of clones exhibiting resistance to carbapenems, a class of critical last-line antibiotics, is therefore closely monitored. Carbapenem resistance is frequently mediated by chromosomal mutations that lead to a defective outer membrane porin OprD. Here, we determined the genetic diversity of OprD variants across the P. aeruginosa population and developed PorinPredict, a bioinformatic tool that enables the prediction of OprD loss from whole-genome sequencing data. We show a high correlation between predicted OprD loss and meropenem nonsusceptibility irrespective of the presence of carbapenemases, which are a second widespread determinant of carbapenem resistance. Isolates with resistance determinants to other antibiotics were disproportionally affected by OprD loss, possibly due to an increased exposure to carbapenems. Integration of PorinPredict into genomic surveillance platforms will facilitate a better understanding of the clinical impact of OprD modifications and transmission dynamics of resistant clones.

Galaxy ASIST: A web-based platform for mapping and assessment of global standards of antimicrobial susceptibility: A case study in Acinetobacter baumannii genomes

Introduction

Antimicrobial susceptibility testing (AST) is used to determine the susceptibility of an organism to antibiotics. The determination of susceptibility is based on MIC breakpoints and is provided by EUCAST and CLSI. Likewise, phenotypic classification criteria developed by CDC/ECDC are used for the classification of pathogens into susceptible, multidrug-resistant, extremely drug-resistant, or totally drug-resistant categories. Whole-genome sequencing (WGS)-based diagnosis is now supplementing existing gold-standard microbiology methods for rapid and more precise AST, and therefore, EUCAST recommended quality criteria to assess whole-genome sequence for reporting the same. In this study, these three global standards, MIC breakpoints, phenotypic classification, and genome quality, are applied to the largest publicly available data for Acinetobacter baumannii (AB), the most critical priority pathogen identified by WHO.

Materials and Methods

The drug sensitivity profile and genomes for isolates of AB were obtained from PATRIC and evaluated with respect to AST standards (CLSI and EUCAST). Whole genome quality assessment and antimicrobial resistance mapping is performed with QUAST and ABRicate, respectively. Four in-house methods are developed for mapping standards and are integrated into a Galaxy workflow based system, Galaxy-ASIST. Analysis of the extent of agreement between CLSI 2022 and EUCAST 2022 for antibiotics was carried out using Cohen's kappa statistics.

Results and Discussion

An automated pipeline, Galaxy-ASIST, is designed and developed for the characterization of clinical isolates based on these standards. Evaluation of over 6,500 AB strains using Galaxy-ASIST indicated that only 10% of the publicly available datasets have metadata to implement these standards. Furthermore, given that CLSI and EUCAST have different MIC breakpoints, discrepancies are observed in the classification of resistant and susceptible isolates following these standards. It is, therefore, imperative that platforms are developed that allow the evaluation of ever increasing phenotypic and genome sequence datasets for AST. Galaxy-ASIST offers a centralized repository and a structured metadata architecture to provide a single globally acceptable framework for AST profiling of clinical isolates based on global standards. The platform also offers subsequent fine mapping of antimicrobial-resistant determinants. Galaxy-ASIST is freely available at https://ab-openlab.csir.res.in/asist.

Pre-epidemic evolution of the MRSA USA300 clade and a molecular key for classification

Introduction

USA300 has remained the dominant community and healthcare associated methicillin-resistant Staphylococcus aureus (MRSA) clone in the United States and in northern South America for at least the past 20 years. In this time, it has experienced epidemic spread in both of these locations. However, its pre-epidemic evolutionary history and origins are incompletely understood. Large sequencing databases, such as NCBI, PATRIC, and Staphopia, contain clues to the early evolution of USA300 in the form of sequenced genomes of USA300 isolates that are representative of lineages that diverged prior to the establishment of the South American epidemic (SAE) clade and North American epidemic (NAE) clade. In addition, historical isolates collected prior to the emergence of epidemics can help reconstruct early events in the history of this lineage.

Methods

Here, we take advantage of the accrued, publicly available data, as well as two newly sequenced pre-epidemic historical isolates from 1996, and a very early diverging ACME-negative NAE genome, to understand the pre-epidemic evolution of USA300. We use database mining techniques to emphasize genomes similar to pre-epidemic isolates, with the goal of reconstructing the early molecular evolution of the USA300 lineage.

Results

Phylogenetic analysis with these genomes confirms that the NAE and SAE USA300 lineages diverged from a most recent common ancestor around 1970 with high confidence, and it also pinpoints the independent acquisition events of the of the ACME and COMER loci with greater precision than in previous studies. We provide evidence for a North American origin of the USA300 lineage and identify multiple introductions of USA300 into South and North America. Notably, we describe a third major USA300 clade (the pre-epidemic branching clade; PEB1) consisting of both MSSA and MRSA isolates circulating around the world that diverged from the USA300 lineage prior to the establishment of the South and North American epidemics. We present a detailed analysis of specific sequence characteristics of each of the major clades, and present diagnostic positions that can be used to classify new genomes.

Draft Genome Sequences of Probiotic Candidate Schleiferilactobacillus harbinensis Isolated from Fermented Milk "Laban"

Lactic acid bacteria are a major component of dairy products, especially species belonging to the genus Lactobacillus. This study reports the whole-genome sequence of Schleiferilactobacillus harbinensis isolated from laban, indigenous fermented milk of Saudi Arabia. The genome sequence is 3,023,618 bp long, has 179 contigs, and has a G+C content of 53.3%.

Complete Genomic Sequence of Bacillus coagulans Strain JBI-YZ6.3: a Natural Spore-Forming Isolate from Food-Grade Tapioca Starch

Bacillus coagulans strain JBI-YZ6.3 is a safe probiotic bacterium isolated from food-grade tapioca starch. The complete genome of B. coagulans JBI-YZ6.3 comprises one circular chromosome of 3.5 Mb and contains no toxigenic and antibiotic resistance genes, providing molecular information to support the strain's safety and usage as a probiotic.

Characterization and comprehensive genome analysis of novel bacteriophage, vB_Kpn_ZCKp20p, with lytic and anti-biofilm potential against clinical multidrug-resistant Klebsiella pneumoniae

Introduction

The rise of infections by antibiotic-resistant bacterial pathogens is alarming. Among these, Klebsiella pneumoniae is a leading cause of death by hospital-acquired infections, and its multidrug-resistant strains are flagged as a global threat to human health, which necessitates finding novel antibiotics or alternative therapies. One promising therapeutic alternative is the use of virulent bacteriophages, which specifically target bacteria and coevolve with them to overcome potential resistance. Here, we aimed to discover specific bacteriophages with therapeutic potential against multiresistant K. pneumoniae clinical isolates.

Methods and Results

Out of six bacteriophages that we isolated from urban and medical sewage, phage vB_Kpn_ZCKp20p had the broadest host range and was thus characterized in detail. Transmission electron microscopy suggests vB_Kpn_ZCKp20p to be a tailed phage of the siphoviral morphotype. In vitro evaluation indicated a high lytic efficiency (30 min latent period and burst size of ∼100 PFU/cell), and extended stability at temperatures up to 70°C and a wide range of (2-12) pH. Additionally, phage vB_Kpn_ZCKp20p possesses antibiofilm activity that was evaluated by the crystal violet assay and was not cytotoxic to human skin fibroblasts. The whole genome was sequenced and annotated, uncovering one tRNA gene and 33 genes encoding proteins with assigned functions out of 85 predicted genes. Furthermore, comparative genomics and phylogenetic analysis suggest that vB_Kpn_ZCKp20p most likely represents a new species, but belongs to the same genus as Klebsiella phages ZCKP8 and 6691. Comprehensive genomic and bioinformatics analyses substantiate the safety of the phage and its strictly lytic lifestyle.

Conclusion

Phage vB_Kpn_ZCKp20p is a novel phage with potential to be used against biofilm-forming K. pneumoniae and could be a promising source for antibacterial and antibiofilm products, which will be individually studied experimentally in future studies.

Molecular insights into novel environmental strains of Klebsiella quasipneumoniae harboring different antimicrobial-resistance genes

Introduction

The emergence of bacterial pathogens in environmental hosts represents a major risk to public health. This study aimed at characterizing seven novel environmental strains of K. quasipneumoniae using a genomic approach which was misidentified by phenotypic methods in a previous batch of 27 species thought to be K. pneumoniae.

Methods

Whole-genome sequencing was performed using the Illumina platform, and the generated raw reads were de novo assembled. Comparative genomic, resistome, virulome, mobilome, and phylogeny were then investigated using dierent bioinformatics tools.

Results

Six strains were identified as K. quasipneumoniae subsp similipneumoniae and one as K. quasipneumoniae subsp. quasipneumoniae. All isolates were resistant to ampicillin, cephalexin, and amoxicillin-clavulanic acid and harbored the fosA, bla _OKP types, oqxB, and oqxA genes. One isolate additionally harbored a gene cassettes consisting of bla _SHV-1, bla _OXA-1, aac(6')-Ib-cr, catB genes. The aminoglycoside-modifying enzyme gene aph(3")-Ia was bracketed by two insertion elements. Plasmid analyses showed that IncFIB_K was the most prevalent plasmid, circulating in six isolates, while one isolate exhibited seven different plasmids. The isolates have virulence genes responsible for capsule formation, lipopolysaccharide, iron uptake aerobactin (iutA), salmochelins (iroE, iroN), enterobactin siderophore, adherence, and biofilm formation (mrkA, mrkB, mrkC, mrkD, mrkF, and mrkH).

Conclusion

Our study highlights the ecology and transmission of K. quasipneumoniae (which have the ability to disseminate to other environmental sources including animals) outside the clinical setting and the contribution of water, vegetables, and table surfaces as potential reservoirs of farm-to-fork transmission of disease via local markets in Khartoum, Sudan.

Exploring antibiotic resistance genes, mobile gene elements, and virulence gene factors in an urban freshwater samples using metagenomic analysis

Antibiotic resistance genes (ARGs) and antimicrobial resistance elements (AMR) are novel environmental contaminants that pose a significant risk to human health globally. Freshwater contains a variety of microorganisms that might affect human health; its quality must be assessed before use. However, the dynamics of mobile genetic elements (MGEs) and ARG propagation in freshwater have rarely been studied in Singapore. Therefore, this study used metagenomics to compare diversity, virulence factor composition, and ARG and MGE co-occurrence with bacterial communities in paired (n = 8) environmental freshwater samples. KneadData, FMAP, and Kraken2 were used for bioinformatics analysis and R (v4.1.1) for statistical analysis. Sequence reads with a total of 9043 species were taxonomically classified into 66 phyla, 130 classes, 261 orders, 584 families, and 2477 genera. Proteobacteria, Bacteroidetes, Actinobacteria, and Firmicutes were found the Phyla in all samples. Analysis of QIIME output by PICRUSt and ß-diversity showed unique clusters and functional microbial community structures. A total of 2961 ARGs were found that conferred resistance to multidrug, aminoglycosides, tetracyclines, elfamycins, and more. The classified ARG mechanism revealed significant distribution of virulence factors in bacterial cells. Transposes and transposon were highly correlated to ARG gene transfer. Co-occurrence network analysis showed several MGEs appear to use the same ARGs (intI and rho) and were dominant in all samples. Furthermore, ARGs are also highly correlated with bacteria like Campylobacter and Escherichia. This study enhances the understanding of antibiotic risk assessment and provides a new perspective on bacterial assembly contamination and the functional prevalence of ARGs and MGEs with antibiotic resistance bacteria. Moreover, it raises public awareness because these contaminants put people's lives at risk of acquiring bacterial infections. In addition, it can also help propose hybrid water treatment approaches.

Genomic potential for exopolysaccharide production and differential polysaccharide degradation in closely related Alteromonas sp. PRIM-21 and Alteromonas fortis 1T

Members of the genus Alteromonas are widely distributed in diverse marine environments and are often associated with marine organisms. Their ability to produce exopolysaccharides (EPS) and depolymerize sulfated algal polysaccharides has provided industrial importance to some species. Here, we describe the draft genome of an algae-associated strain namely, Alteromonas sp. PRIM-21 isolated from the southwest coast of India to understand the EPS biosynthetic pathways as well as polysaccharide depolymerization system in comparison to the closely related strain Alteromonas fortis 1^T that shares 99.8% 16S rRNA gene sequence similarity. Whole-genome shotgun sequencing of Alteromonas sp. PRIM-21 yielded 50 contigs with a total length of 4,638,422 bp having 43.86% GC content. The resultant genome shared 95.9% OrthoANI value with A. fortis 1 ^T, and contained 4125 predicted protein-coding genes, 71 tRNA and 10 rRNA genes. Genes involved in Wzx/Wzy-, ABC transporter- and synthase-dependent pathways for EPS production and secretion were common in both Alteromonas sp. PRIM-21 and A. fortis 1^T. However, the distribution of carbohydrate-active enzymes (CAZymes) was heterogeneous. The strain PRIM-21 harbored polysaccharide lyases for the degradation of alginate, ulvan, arabinogalactan and chondroitin. This was further validated from the culture-based assays using seven different polysaccharides. The depolymerizing ability of the bacteria may be useful in deriving nutrients from the biopolymers produced in the algal host while the EPS biosynthesis may provide additional advantages for life in the stressful marine environment. The results also highlight the genetic heterogeneity in terms of polysaccharide utilization among the closely related Alteromonas strains.

Robinsoniella peoriensis: an emerging pathogen with few virulence factors

Robinsoniella peoriensis is a Gram-positive bacterium which is anaerobic, spore-forming, and non-motile. It was initially isolated and characterized from feces and swine manure. Strains of this species have since been identified from different mammalian and non-mammalian gastrointestinal tracts. Strains have also been isolated from a variety of human infections, such as bacteremia, bone infections, and skin structures. R. peoriensis has recently been reported as causative for pyometra, which could result in death in the absence of sufficient antimicrobial treatment. However, to the author's knowledge, there has not been a single virulence factor identified. A major challenge of modern medicine is the failure of conventional procedures to characterize the capability of an emerging pathogen to cause disease. The goal of this study is to initially characterize the pathogenicity of this bacterium using a pathogenomics approach.

Isolation and characterization of a novel Tenacibaculum species and a corresponding bacteriophage from a Mediterranean fish hatchery: Description of Tenacibaculum larymnensis sp. nov. and Tenacibaculum phage Larrie

Tenacibaculum larymnensis sp. nov., a novel species of the Tenacibaculum genus was isolated from a commercial fish hatchery in Greece. The novel species is phylogenetically close to T. discolor and was biochemically and genetically characterized. The genome of T. larymnensis has 3.66 Mbps length, 31.83% GC content and the genomic analysis demonstrated that it harbors a wide enzymatic repertoire suggestive of increased degrading capacity but also several virulence factors including hemolysins, secretion systems, transporters, siderophores, pili and extracellular proteins. Using the novel strain, a virulent bacteriophage designated as Tenacibaculum phage Larrie was isolated and characterized. Larrie is a novel Siphovirus with relatively large genome, 77.5 kbps with 111 ORFs, a GC content of 33.7% and an exclusively lytic lifestyle. The new phage-host system can serve as an efficient model to study microbial interactions in the aquatic environment which contribute to the nutrient cycling.

The K-mer antibiotic resistance gene variant analyzer (KARGVA)

Characterization of antibiotic resistance genes (ARGs) from high-throughput sequencing data of metagenomics and cultured bacterial samples is a challenging task, with the need to account for both computational (e.g., string algorithms) and biological (e.g., gene transfers, rearrangements) aspects. Curated ARG databases exist together with assorted ARG classification approaches (e.g., database alignment, machine learning). Besides ARGs that naturally occur in bacterial strains or are acquired through mobile elements, there are chromosomal genes that can render a bacterium resistant to antibiotics through point mutations, i.e., ARG variants (ARGVs). While ARG repositories also collect ARGVs, there are only a few tools that are able to identify ARGVs from metagenomics and high throughput sequencing data, with a number of limitations (e.g., pre-assembly, a posteriori verification of mutations, or specification of species). In this work we present the k-mer, i.e., strings of fixed length k, ARGV analyzer - KARGVA - an open-source, multi-platform tool that provides: (i) an ad hoc, large ARGV database derived from multiple sources; (ii) input capability for various types of high-throughput sequencing data; (iii) a three-way, hash-based, k-mer search setup to process data efficiently, linking k-mers to ARGVs, k-mers to point mutations, and ARGVs to k-mers, respectively; (iv) a statistical filter on sequence classification to reduce type I and II errors. On semi-synthetic data, KARGVA provides very high accuracy even in presence of high sequencing errors or mutations (99.2 and 86.6% accuracy within 1 and 5% base change rates, respectively), and genome rearrangements (98.2% accuracy), with robust performance on ad hoc false positive sets. On data from the worldwide MetaSUB consortium, comprising 3,700+ metagenomics experiments, KARGVA identifies more ARGVs than Resistance Gene Identifier (4.8x) and PointFinder (6.8x), yet all predictions are below the expected false positive estimates. The prevalence of ARGVs is correlated to ARGs but ecological characteristics do not explain well ARGV variance. KARGVA is publicly available at https://github.com/DataIntellSystLab/KARGVA under MIT license.

Stingray epidermal microbiomes are species-specific with local adaptations

Marine host-associated microbiomes are affected by a combination of species-specific (e.g., host ancestry, genotype) and habitat-specific features (e.g., environmental physiochemistry and microbial biogeography). The stingray epidermis provides a gradient of characteristics from high dermal denticles coverage with low mucus to reduce dermal denticles and high levels of mucus. Here we investigate the effects of host phylogeny and habitat by comparing the epidermal microbiomes of Myliobatis californica (bat rays) with a mucus rich epidermis, and Urobatis halleri (round rays) with a mucus reduced epidermis from two locations, Los Angeles and San Diego, California (a 150 km distance). We found that host microbiomes are species-specific and distinct from the water column, however composition of M. californica microbiomes showed more variability between individuals compared to U. halleri. The variability in the microbiome of M. californica caused the microbial taxa to be similar across locations, while U. halleri microbiomes were distinct across locations. Despite taxonomic differences, Shannon diversity is the same across the two locations in U. halleri microbiomes suggesting the taxonomic composition are locally adapted, but diversity is maintained by the host. Myliobatis californica and U. halleri microbiomes maintain functional similarity across Los Angeles and San Diego and each ray showed several unique functional genes. Myliobatis californica has a greater relative abundance of RNA Polymerase III-like genes in the microbiome than U. halleri, suggesting specific adaptations to a heavy mucus environment. Construction of Metagenome Assembled Genomes (MAGs) identified novel microbial species within Rhodobacteraceae, Moraxellaceae, Caulobacteraceae, Alcanivoracaceae and Gammaproteobacteria. All MAGs had a high abundance of active RNA processing genes, heavy metal, and antibiotic resistant genes, suggesting the stingray mucus supports high microbial growth rates, which may drive high levels of competition within the microbiomes increasing the antimicrobial properties of the microbes.

Shotgun Metagenomics of Gastric Biopsies Reveals Compositional and Functional Microbiome Shifts in High- and Low-Gastric-Cancer-Risk Populations from Colombia, South America

Along with Helicobacter pylori infection, the gastric microbiota is hypothesized to modulate stomach cancer risk in susceptible individuals. Whole metagenomic shotgun sequencing (WMS) is a sequencing approach to characterize the microbiome with advantages over traditional culture and 16S rRNA sequencing including identification of bacterial and non-bacterial taxa, species/strain resolution, and functional characterization of the microbiota. In this study, we used WMS to survey the microbiome in extracted DNA from antral gastric biopsy samples from Colombian patients residing in the high-risk gastric cancer town Túquerres (n = 10, H. pylori-positive = 7) and low-risk town of Tumaco (n = 10, H. pylori-positive = 6). Kraken2/Bracken was used for taxonomic classification and abundance. Functional gene profiles were inferred by InterProScan and KEGG analysis of assembled contigs and gene annotation. The most abundant taxa represented bacteria, non-human eukaryota, and viral genera found in skin, oral, food, and plant/soil environments including Staphylococus, Streptococcus, Bacillus, Aspergillus, and Siphoviridae. H. pylori was the predominant taxa present in H. pylori-positive samples. Beta diversity was significantly different based on H. pylori-status, risk group, and sex. WMS detected more bacterial taxa than 16S rRNA sequencing and aerobic, anaerobic, and microaerobic culture performed on the same gastric biopsy samples. WMS identified significant differences in functional profiles found between H. pylori-status, but not risk or sex groups. H. pylori-positive samples were significantly enriched for H. pylori-specific genes including virulence factors such as vacA, cagA, and urease, while carbohydrate and amino acid metabolism genes were enriched in H. pylori-negative samples. This study shows WMS has the potential to characterize the taxonomy and function of the gastric microbiome as risk factors for H. pylori-associated gastric disease. Future studies will be needed to compare and validate WMS versus traditional culture and 16S rRNA sequencing approaches for characterization of the gastric microbiome.

Mucilaginibacter aquariorum sp. nov., Isolated from Fresh Water

A Gram-stain-negative, rod-shaped bacterial strain, JC4^T, was isolated from a freshwater sample and determined the taxonomic position. Initial identification based on 16S rRNA gene sequences revealed that strain JC4^T is affiliated to the genus Mucilaginibacter with a sequence similarity of 97.97% to Mucilaginibacter rigui WPCB133^T. The average nucleotide identity and digital DNA-DNA hybridization values between strain JC4^T and Mucilaginibacter species were estimated below 80.92% and 23.9%, respectively. Strain JC4^T contained summed feature 3 (C_16:1 ω7c and/or C_16:1 ω6c) and iso-C_15:0 as predominant cellular fatty acids. The dominant polar lipids were identified as phosphatidylethanolamine, one unidentified aminophospholipid, one unidentified phospholipid, and two unidentified lipids. The respiratory quinone was MK-7. The genomic DNA G+C content of strain JC4^T was determined to be 42.44%. The above polyphasic evidences support that strain JC4^T represents a novel species of the genus Mucilaginibacter, for which the name Mucilaginibacter aquariorum sp. nov. is proposed. The type strain is JC4^T (= KCTC 92230^T = LMG 32715^T).

SourceFinder: a Machine-Learning-Based Tool for Identification of Chromosomal, Plasmid, and Bacteriophage Sequences from Assemblies

High-throughput genome sequencing technologies enable the investigation of complex genetic interactions, including the horizontal gene transfer of plasmids and bacteriophages. However, identifying these elements from assembled reads remains challenging due to genome sequence plasticity and the difficulty in assembling complete sequences. In this study, we developed a classifier, using random forest, to identify whether sequences originated from bacterial chromosomes, plasmids, or bacteriophages. The classifier was trained on a diverse collection of 23,211 chromosomal, plasmid, and bacteriophage sequences from hundreds of bacterial species. In order to adapt the classifier to incomplete sequences, each complete sequence was subsampled into 5,000 nucleotide fragments and further subdivided into k-mers. This three-class classifier succeeded in identifying chromosomes, plasmids, and bacteriophages using k-mer distributions of complete and partial genome sequences, including simulated metagenomic scaffolds with minimum performance of 0.939 area under the receiver operating characteristic curve (AUC). This classifier, implemented as SourceFinder, has been made available as an online web service to help the community with predicting the chromosomal, plasmid, and bacteriophage sources of assembled bacterial sequence data (https://cge.food.dtu.dk/services/SourceFinder/). IMPORTANCE Extra-chromosomal genes encoding antimicrobial resistance, metal resistance, and virulence provide selective advantages for bacterial survival under stress conditions and pose serious threats to human and animal health. These accessory genes can impact the composition of microbiomes by providing selective advantages to their hosts. Accurately identifying extra-chromosomal elements in genome sequence data are critical for understanding gene dissemination trajectories and taking preventative measures. Therefore, in this study, we developed a random forest classifier for identifying the source of bacterial chromosomal, plasmid, and bacteriophage sequences.

Type IV Pili Are a Critical Virulence Factor in Clinical Isolates of Paenibacillus thiaminolyticus

Hydrocephalus, the leading indication for childhood neurosurgery worldwide, is particularly prevalent in low- and middle-income countries. Hydrocephalus preceded by an infection, or postinfectious hydrocephalus, accounts for up to 60% of hydrocephalus in these areas. Since many children with hydrocephalus suffer poor long-term outcomes despite surgical intervention, prevention of hydrocephalus remains paramount. Our previous studies implicated a novel bacterial pathogen, Paenibacillus thiaminolyticus, as a causal agent of neonatal sepsis and postinfectious hydrocephalus in Uganda. Here, we report the isolation of three P. thiaminolyticus strains, Mbale, Mbale2, and Mbale3, from patients with postinfectious hydrocephalus. We constructed complete genome assemblies of the clinical isolates as well as the nonpathogenic P. thiaminolyticus reference strain and performed comparative genomic and proteomic analyses to identify potential virulence factors. All three isolates carry a unique beta-lactamase gene, and two of the three isolates exhibit resistance in culture to the beta-lactam antibiotics penicillin and ampicillin. In addition, a cluster of genes carried on a mobile genetic element that encodes a putative type IV pilus operon is present in all three clinical isolates but absent in the reference strain. CRISPR-mediated deletion of the gene cluster substantially reduced the virulence of the Mbale strain in mice. Comparative proteogenomic analysis identified various additional potential virulence factors likely acquired on mobile genetic elements in the virulent strains. These results provide insight into the emergence of virulence in P. thiaminolyticus and suggest avenues for the diagnosis and treatment of this novel bacterial pathogen. IMPORTANCE Postinfectious hydrocephalus, a devastating sequela of neonatal infection, is associated with increased childhood mortality and morbidity. A novel bacterial pathogen, Paenibacillus thiaminolyticus, is highly associated with postinfectious hydrocephalus in an African cohort. Whole-genome sequencing, RNA sequencing, and proteomics of clinical isolates and a reference strain in combination with CRISPR editing identified type IV pili as a critical virulence factor for P. thiaminolyticus infection. Acquisition of a type IV pilus-encoding mobile genetic element critically contributed to converting a nonpathogenic strain of P. thiaminolyticus into a pathogen capable of causing devastating diseases. Given the widespread presence of type IV pilus in pathogens, the presence of the type IV pilus operon could serve as a diagnostic and therapeutic target in P. thiaminolyticus and related bacteria.

Classification of bacterial plasmid and chromosome derived sequences using machine learning

Plasmids are important genetic elements that facilitate horizonal gene transfer between bacteria and contribute to the spread of virulence and antimicrobial resistance. Most bacterial genome sequences in the public archives exist in draft form with many contigs, making it difficult to determine if a contig is of chromosomal or plasmid origin. Using a training set of contigs comprising 10,584 chromosomes and 10,654 plasmids from the PATRIC database, we evaluated several machine learning models including random forest, logistic regression, XGBoost, and a neural network for their ability to classify chromosomal and plasmid sequences using nucleotide k-mers as features. Based on the methods tested, a neural network model that used nucleotide 6-mers as features that was trained on randomly selected chromosomal and plasmid subsequences 5kb in length achieved the best performance, outperforming existing out-of-the-box methods, with an average accuracy of 89.38% ± 2.16% over a 10-fold cross validation. The model accuracy can be improved to 92.08% by using a voting strategy when classifying holdout sequences. In both plasmids and chromosomes, subsequences encoding functions involved in horizontal gene transfer-including hypothetical proteins, transporters, phage, mobile elements, and CRISPR elements-were most likely to be misclassified by the model. This study provides a straightforward approach for identifying plasmid-encoding sequences in short read assemblies without the need for sequence alignment-based tools.

Characterizing the role of phosphatidylglycerol-phosphate phosphatases in Acinetobacter baumannii cell envelope biogenesis and antibiotic resistance

The dissemination of multi-drug resistant Acinetobacter baumannii threatens global healthcare systems and necessitates the development of novel therapeutic options. The Gram-negative bacterial cell envelope provides a first defensive barrier against antimicrobial assault. Essential components of this multi-layered complex are the phospholipid-rich membranes. Phosphatidylglycerol phosphate (PGP) phosphatases are responsible for a key step in the biosynthesis of a major phospholipid species, phosphatidylglycerol (PG), but these enzymes have also been implicated in the biogenesis of other cell envelope components. Our bioinformatics analyses identified two putative PGP candidates in the A. baumannii genome, PgpA and PgpB. Phospholipid analyses of isogenic pgpA mutants in two distinct A. baumannii strains revealed a shift in the desaturation levels of phosphatidylethanolamine (PE) phospholipid species, possibly due to the activation of the phospholipid desaturase DesA. We also investigated the impact of the inner membrane phosphatases on other cell envelope components, which revealed a role of PgpB in the maintenance of the A. baumannii peptidoglycan layer, and consequently carbapenem resistance. Collectively, this work provides novel insights into the roles of PGP phosphatases on the global lipidomic landscape of A. baumannii and their interconnectivity with the biogenesis of other cell envelope components. The non-essentiality of these candidates exemplifies metabolic versatility of A. baumannii, which is believed to be key to its success as global pathogen.

Genomic landscape of prominent XDR Acinetobacter clonal complexes from Dhaka, Bangladesh

BACKGROUND

Acinetobacter calcoaceticus-A. baumannii (ACB) complex pathogens are known for their prevalence in nosocomial infections and extensive antimicrobial resistance (AMR) capabilities. While genomic studies worldwide have elucidated the genetic context of antibiotic resistance in major international clones (ICs) of clinical Acinetobacter spp., not much information is available from Bangladesh. In this study, we analysed the AMR profiles of 63 ACB complex strains collected from Dhaka, Bangladesh. Following this, we generated draft genomes of 15 of these strains to understand the prevalence and genomic environments of AMR, virulence and mobilization associated genes in different Acinetobacter clones.

RESULTS

Around 84% (n = 53) of the strains were extensively drug resistant (XDR) with two showing pan-drug resistance. Draft genomes generated for 15 strains confirmed 14 to be A. baumannii while one was A. nosocomialis. Most A. baumannii genomes fell under three clonal complexes (CCs): the globally dominant CC1 and CC2, and CC10; one strain had a novel sequence type (ST). AMR phenotype-genotype agreement was observed and the genomes contained various beta-lactamase genes including bla_OXA-23 (n = 12), bla_OXA-66 (n = 6), and bla_NDM-1 (n = 3). All genomes displayed roughly similar virulomes, however some virulence genes such as the Acinetobactin bauA and the type IV pilus gene pilA displayed high genetic variability. CC2 strains carried highest levels of plasmidic gene content and possessed conjugative elements carrying AMR genes, virulence factors and insertion sequences.

CONCLUSION

This study presents the first comparative genomic analysis of XDR clinical Acinetobacter spp. from Bangladesh. It highlights the prevalence of different classes of beta-lactamases, mobilome-derived heterogeneity in genetic architecture and virulence gene variability in prominent Acinetobacter clonal complexes in the country. The findings of this study would be valuable in understanding the genomic epidemiology of A. baumannii clones and their association with closely related pathogenic species like A. nosocomialis in Bangladesh.

Actinomyces ruminis sp. nov., an obligately anaerobic bacterium isolated from the rumen of cattle

An obligately anaerobic, rod-shaped, Gram-stain-positive, non-spore-forming, non-motile bacterial strain; designated as CtC72^T was isolated from the rumen of cattle. The 16S rRNA gene sequence similarity of less than 98.65% revealed the strain as a member of the genus Actinomyces, nearest to but distinct from Actinomyces qiguomingii DSM 106201^T, Actinomyces ruminicola DSM 27982^T, Actinomyces procaprae JCM 33484^T, Actinomyces succiniciruminis TISTR 2317, Actinomyces glycerinitolerans TISTR 2318. The low values of digital DNA-DNA hybridization (< 70%) and average nucleotide identity (< 95%) further highlighted the distinctive nature of strain CtC72^T from its closest relatives. The strain CtC72^T could grow at temperatures between 30 and 50 °C (optimum 40 °C), pH between 6.0 and 9.0 (optimum 7.5-8.0), and NaCl between 0 and 1.5% (optimum 0%). The strain hydrolysed cellulose and xylan and utilised a range of mono-, di-, and oligo-saccharides as a source of carbon and energy. Glucose fermentation resulted in acetic acid and formic acid as major metabolic products, while propionic acid, lactic acid, and ethanol as minor products along with CO₂ production. The DNA G + C content of strain CtC72^T was 68.40 (mol%, T_m) and 68.05 (%, digital). Major cellular fatty acids (> 10%) were C_16:0, C_18:1 ω9c_, and C_18:1 ω9c DMA. Based on these data, we propose that strain CtC72^T be classified as a novel species, Actinomyces ruminis sp. nov., under the genus Actinomyces. The type strain is CtC72^T (= KCTC 15726^T = JCM 32641^T = MCC 3500^T).