Identification of Protein Secretion Systems in Bacterial Genomes Using MacSyFinder

Comparative genomic analysis of Flavobacterium species causing columnaris disease of freshwater fish in Thailand: insights into virulence and resistance mechanisms

BACKGROUND

Columnaris disease, a prevalent disease among farmed and wild freshwater fish, is caused by the Flavobacterium columnare group, which includes four distinct species: F. columnare, F. oreochromis, F. covae, and F. davisii. Among these, F. oreochromis, F. covae, and F. davisii are particularly prevalent in farmed freshwater fish in Thailand. In this study, a comparative genomic analysis of 22 isolates was conducted to elucidate virulence factors, antibiotic resistance genes (ARGs), genomic islands (GIs), phages, insertion elements (ISs), and clustered regularly interspaced short palindromic repeats (CRISPRs).

RESULTS

A total of 212 putative virulence genes were predicted across three species with F. oreochromis exhibiting the highest number of unique virulence genes, followed by F. davisii, and F. covae. Moreover, 195 genes were predicted as ARGs, with F. oreochromis and F. covae showing an abundance of unique genes associated with resistance to quinolone, fluoroquinolone, and tetracycline antibiotics. Antimicrobial susceptibility testing, assessed with epidemiological cut-off values (ECVs), revealed decreased susceptibility to quinolones, fluoroquinolones and tetracycline in several isolates of F. oreochromis and F. covae. F. oreochromis and F. covae exhibited notable decreased susceptibility to quinolones, with mutations observed in the quinolone resistance-determining regions (QRDRs) of gyrA, including Ser83Phe, Ser83Val, Ser83Ala, and Asp87Tyr, the latter representing a novel mutation among isolates from Thailand. As a result, these findings suggest that gyrA is major target for quinolone in F. oreochromis, F. covae, and F. davisii, while gyrB, parC, parE might be less important to the decreased phenotypic susceptibility to this class of antimicrobials. Moreover, a tetracycline resistance gene (tetA_2) was found in only one F. covae isolate, which exhibited decreased phenotypic susceptibility to this drug, marking the first report of decreased susceptibility in this species.

CONCLUSIONS

This study provides insights into the genetic and pathogenic diversity of Flavobacterium species, aiding in the development of strategies to manage columnaris disease in farmed freshwater fish in Thailand.

Emergence of the mobile RND-type efflux pump gene cluster tmexCD1-toprJ1 in Klebsiella pneumoniae clinical isolates in Japan

BACKGROUND

Tigecycline is an antimicrobial agent with a broad spectrum of activity against both Gram-positive and Gram-negative bacteria. However, mobile tigecycline resistance gene clusters, such as tnfxB-tmexCD-toprJ, have spread globally. The prevalence of tigecycline-resistant Enterobacterales in clinical settings in Japan is unknown.

OBJECTIVES

To investigate the tnfxB-tmexCD-toprJ gene cluster in the genome sequences of Enterobacterales clinical isolates in Japan.

METHODS

We investigated the tnfxB-tmexCD-toprJ cluster from the genome sequences of 5143 Enterobacterales isolates collected from 175 hospitals around Japan between 2019 and 2020 as part of a national genomic surveillance program for antimicrobial-resistant bacteria.

RESULTS

The tnfxB1-tmexCD1-toprJ1 cluster was detected in two Klebsiella pneumoniae isolates in 2019. One isolate possessed a 299.4 kb IncFIB(K) plasmid, pJBBGAAF19431, and the other possessed a 224.9 kb IncHI1B/IncFIB(K) hybrid plasmid, pJBEAACG19501, co-carrying multiple antimicrobial resistance genes, including extended-spectrum β-lactamase genes, blaOXA-1 and blaCTX-M-27, respectively, along with tnfxB1-tmexCD1-toprJ1. The genetic context of the tnfxB1-tmexCD1-toprJ1-surrounding structure on pJBBGAAF19431 was similar to that of a K. pneumoniae plasmid pHNAH8I-1 from a chicken in China in 2017, and the cluster was embedded in an apparently intact mobile DNA element: strand-biased circularizing integrative element. The tnfxB1-tmexCD1-toprJ1 on pJBEAACG19501 was embedded in a Tn3 family transposon related to TnAs1. The plasmid pJBEAACG19501 was highly similar to that of K. pneumoniae, isolated from humans in China in 2021.

CONCLUSIONS

tmexCD-toprJ was present in Japan as of 2019. Even in Japan, where the clinical use of tigecycline is significantly rare, tmexCD-toprJ-harbouring multidrug-resistant Enterobacterales is a public health threat and requires continuous monitoring.

New, Complete Circularized Genomes of Xanthomonas citri pv. mangiferaeindicae Produced from Short- and Long-Read Co-Assembly Shed Light on Strains that Emerged a Decade Ago on Mango and Cashew in Burkina Faso

We report high-quality genomes of three strains of Xanthomonas citri pv. mangiferaeindicae, the causal agent of mango bacterial canker disease, including the pathotype strain of this pathovar and two strains from Burkina Faso that emerged a decade ago. These strains hosted two to three plasmids of sizes ranging from 19 to 86 kb. Genome mining revealed the presence of several secretion systems and effectors involved in the virulence of xanthomonads with (i) a type I secretion system of the hlyDB group; (ii) xps and xcs type II secretion systems; (iii) a type III secretion system with several type III effectors, including transcription activator-like effectors; (iv) several type IV secretion systems associated with plasmid or integrative conjugative elements mobility; (v) three type V secretion system subclasses (Va, Vb, and Vc); and (vi) a single i3* type VI secretion system. The two strains isolated in Burkina Faso from mango (Mangifera indica) and cashew (Anacardium occidentale) differed by only 14 single-nucleotide polymorphisms and shared identical secretion systems and type III effector repertoires. Several transcription activator-like effectors were identified in each strain, some of which may target plant genes previously found implicated in disease development in other xanthomonad-associated pathosystems. These results support the emergence in Burkina Faso a decade ago of very closely related strains that became epidemic on mango and cashew (i.e., two distinct host genera of a same plant family). These new genomic resources will contribute to better understanding the biology and evolution of this agriculturally major crop pathogen.

Investigation of Citrobacter freundii clinical isolates in a Chinese hospital during 2020-2022 revealed genomic characterization of an extremely drug-resistant C. freundii ST257 clinical strain GMU8049 co-carrying blaNDM-1 and a novel blaCMY variant

The emergence of multidrug-resistant Citrobacter freundii poses a significant threat to public health. C. freundii isolates were collected from clinical patients in a Chinese hospital during 2020-2022. An unusual strain, GMU8049, was not susceptible to any of the antibiotics tested, including the novel β-lactam/β-lactamase inhibitor combination ceftazidime-avibactam. Whole-genome sequencing (WGS) revealed that GMU8049 harbors a circular chromosome belonging to the rare ST257 and an IncX3 resistance plasmid. Genomic analysis revealed the coexistence of two β-lactamase genes, including plasmid-mediated blaNDM-1 and chromosomal blaCMY encoding a novel CMY variant, combined with an outer membrane porin deficiency, which may account for the extreme resistance to β-lactams. Conjugation experiment confirmed that the blaNDM-1 resistance gene located on pGMU8049 could be successfully transferred to Escherichia coli EC600. The novel CMY variant had an amino acid substitution at position 106 (N106S) compared to the closely related CMY-51. Additionally, a GMU8049-specific truncation in an OmpK37 variant that produces a premature stop codon. Moreover, a variety of chromosome-located efflux pump coding genes and virulence-related genes were also identified. Analysis of strain GMU8049 in the context of other C. freundii strains reveals an open pan-genome and the presence of mobile genetic elements that can mediate horizontal gene transfer of antimicrobial resistance and virulence genes. Our work provides comprehensive insights into the genetic mechanisms of highly resistant C. freundii, highlighting the importance of genomic surveillance of this opportunistic pathogen as a high-risk population for emerging resistance and pathogenicity.IMPORTANCEEmerging pathogens exhibiting multi-, extremely, and pan-drug resistance are a major concern for hospitalized patients and the healthcare community due to limited antimicrobial treatment options and the potential for spread. Genomic technologies have enabled clinical surveillance of emerging pathogens and modeling of the evolution and transmission of antimicrobial resistance and virulence. Here, we report the genomic characterization of an extremely drug-resistant ST257 Citrobacter freundii clinical isolate. Genomic analysis of GMU8049 with a rare ST type and unusual phenotypes can provide information on how this extremely resistant clinical isolate has evolved, including the acquisition of blaNDM-1 via the IncX3 plasmid and accumulation through chromosomal mutations leading to a novel CMY variant and deficiency of the outer membrane porin OmpK37. Our work highlights that the emergence of extremely resistant C. freundii poses a significant challenge to the treatment of clinical infections. Therefore, great efforts must be made to specifically monitor this opportunistic pathogen.

Unveiling intraspecific diversity and evolutionary dynamics of the foodborne pathogen Bacillus paranthracis through high-quality pan-genome analysis

Understanding the evolutionary dynamics of foodborne pathogens throughout host-associated habitats is of utmost importance. Bacterial pan-genomes, as dynamic entities, are strongly influenced by ecological lifestyles. As a phenotypically diverse species in the Bacillus cereus group, Bacillus paranthracis is recognized as an emerging foodborne pathogen and a probiotic simultaneously. This poorly understood species is a suitable study model for adaptive pan-genome evolution. In this study, we determined the biogeographic distribution, abundance, genetic diversity, and genotypic profiles of key genetic elements of B. paranthracis. Metagenomic read recruitment analyses demonstrated that B. paranthracis members are globally distributed and abundant in host-associated habitats. A high-quality pan-genome of B. paranthracis was subsequently constructed to analyze the evolutionary dynamics involved in ecological adaptation comprehensively. The open pan-genome indicated a flexible gene repertoire with extensive genetic diversity. Significant divergences in the phylogenetic relationships, functional enrichment, and degree of selective pressure between the different components demonstrated different evolutionary dynamics between the core and accessory genomes driven by ecological forces. Purifying selection and gene loss are the main signatures of evolutionary dynamics in B. paranthracis pan-genome. The plasticity of the accessory genome is characterized by horizontal gene transfer (HGT), massive gene losses, and weak purifying or positive selection, which might contribute to niche-specific adaptation. In contrast, although the core genome dominantly undergoes purifying selection, its association with HGT and positively selected mutations indicates its potential role in ecological diversification. Furthermore, host fitness-related dynamics are characterized by the loss of secondary metabolite biosynthesis gene clusters (BGCs) and CAZyme-encoding genes and the acquisition of antimicrobial resistance (AMR) and virulence genes via HGT. This study offers a case study of pan-genome evolution to investigate the ecological adaptations reflected by biogeographical characteristics, thereby advancing the understanding of intraspecific diversity and evolutionary dynamics of foodborne pathogens.

Beav: a bacterial genome and mobile element annotation pipeline

Comprehensive and accurate genome annotation is crucial for inferring the predicted functions of an organism. Numerous tools exist to annotate genes, gene clusters, mobile genetic elements, and other diverse features. However, these tools and pipelines can be difficult to install and run, be specialized for a particular element or feature, or lack annotations for larger elements that provide important genomic context. Integrating results across analyses is also important for understanding gene function. To address these challenges, we present the Beav annotation pipeline. Beav is a command-line tool that automates the annotation of bacterial genome sequences, mobile genetic elements, molecular systems and gene clusters, key regulatory features, and other elements. Beav uses existing tools in addition to custom models, scripts, and databases to annotate diverse elements, systems, and sequence features. Custom databases for plant-associated microbes are incorporated to improve annotation of key virulence and symbiosis genes in agriculturally important pathogens and mutualists. Beav includes an optional Agrobacterium-specific pipeline that identifies and classifies oncogenic plasmids and annotates plasmid-specific features. Following the completion of all analyses, annotations are consolidated to produce a single comprehensive output. Finally, Beav generates publication-quality genome and plasmid maps. Beav is on Bioconda and is available for download at https://github.com/weisberglab/beav.

IMPORTANCE

Annotation of genome features, such as the presence of genes and their predicted function, or larger loci encoding secretion systems or biosynthetic gene clusters, is necessary for understanding the functions encoded by an organism. Genomes can also host diverse mobile genetic elements, such as integrative and conjugative elements and/or phages, that are often not annotated by existing pipelines. These elements can horizontally mobilize genes encoding for virulence, antimicrobial resistance, or other adaptive functions and alter the phenotype of an organism. We developed a software pipeline, called Beav, that combines new and existing tools for the comprehensive annotation of these and other major features. Existing pipelines often misannotate loci important for virulence or mutualism in plant-associated bacteria. Beav includes custom databases and optional workflows for the improved annotation of plant-associated bacteria. Beav is designed to be easy to install and run, making comprehensive genome annotation broadly available to the research community.

Pan-genome insights into adaptive evolution of bacterial symbionts in mixed host-microbe symbioses represented by human gut microbiota Bacteroides cellulosilyticus

Animal hosts harbor diverse assemblages of microbial symbionts that play crucial roles in the host's lifestyle. The link between microbial symbiosis and host development remains poorly understood. In particular, little is known about the adaptive evolution of gut bacteria in host-microbe symbioses. Recently, symbiotic relationships have been categorized as open, closed, or mixed, reflecting their modes of inter-host transmission and resulting in distinct genomic features. Members of the genus Bacteroides are the most abundant human gut microbiota and possess both probiotic and pathogenic potential, providing an excellent model for studying pan-genome evolution in symbiotic systems. Here, we determined the complete genome of an novel clinical strain PL2022, which was isolated from a blood sample and performed pan-genome analyses on a representative set of Bacteroides cellulosilyticus strains to quantify the influence of the symbiotic relationship on the evolutionary dynamics. B. cellulosilyticus exhibited correlated genomic features with both open and closed symbioses, suggesting a mixed symbiosis. An open pan-genome is characterized by abundant accessory gene families, potential horizontal gene transfer (HGT), and diverse mobile genetic elements (MGEs), indicating an innovative gene pool, mainly associated with genomic islands and plasmids. However, massive parallel gene loss, weak purifying selection, and accumulation of positively selected mutations were the main drivers of genome reduction in B. cellulosilyticus. Metagenomic read recruitment analyses showed that B. cellulosilyticus members are globally distributed and active in human gut habitats, in line with predominant vertical transmission in the human gut. However, existence and/or high abundance were also detected in non-intestinal tissues, other animal hosts, and non-host environments, indicating occasional horizontal transmission to new niches, thereby creating arenas for the acquisition of novel genes. This case study of adaptive evolution under a mixed host-microbe symbiosis advances our understanding of symbiotic pan-genome evolution. Our results highlight the complexity of genetic evolution in this unusual intestinal symbiont.

Phenotypic and Genomic Characterization of Pseudomonas wuhanensis sp. nov., a Novel Species with Promising Features as a Potential Plant Growth-Promoting and Biocontrol Agent

Plant growth-promoting rhizobacterial strain FP607T was isolated from the rhizosphere of beets in Wuhan, China. Strain FP607T exhibited significant antagonism toward several phytopathogenic bacteria, indicating that FP607T may produce antimicrobial metabolites and has a stronger biocontrol efficacy against plant pathogens. Growth-promoting tests showed that FP607T produced indole-3-acetic acid (IAA), NH3, and ferritin. The genome sequence of strain FP607T was 6,590,972 bp long with 59.0% G + C content. The optimum temperature range was 25-30 °C, and the optimum pH was 7. The cells of strain FP607T were Gram-negative, short, and rod-shaped, with polar flagella. The colonies on the King's B (KB) agar plates were light yellow, smooth, and circular, with regular edges. A phylogenetic analysis of the 16S rRNA sequence and a multilocus sequence analysis (MLSA) showed that strain FP607T was most closely related to the type of strain Pseudomonas farris SWRI79T. Based on a polyphasic taxonomic approach, strain FP607T was identified as a novel species within the genus Pseudomonas, for which the name Pseudomonas wuhanensis sp. nov. was proposed. The type of strain used was FP607T (JCM 35688, CGMCC 27743, and ACCC 62446).

A bioinformatic approach to identify confirmed and probable CRISPR-Cas systems in the Acinetobacter calcoaceticus-Acinetobacter baumannii complex genomes

Introduction

The Acinetobacter calcoaceticus-Acinetobacter baumannii complex, or Acb complex, consists of six species: Acinetobacter baumannii, Acinetobacter calcoaceticus, Acinetobacter nosocomialis, Acinetobacter pittii, Acinetobacter seifertii, and Acinetobacter lactucae. A. baumannii is the most clinically significant of these species and is frequently related to healthcare-associated infections (HCAIs). Clustered regularly interspaced short palindromic repeat (CRISPR) arrays and associated genes (cas) constitute bacterial adaptive immune systems and function as variable genetic elements. This study aimed to conduct a genomic analysis of Acb complex genomes available in databases to describe and characterize CRISPR systems and cas genes.

Methods

Acb complex genomes available in the NCBI and BV-BRC databases, the identification and characterization of CRISPR-Cas systems were performed using CRISPRCasFinder, CRISPRminer, and CRISPRDetect. Sequence types (STs) were determined using the Oxford scheme and ribosomal multilocus sequence typing (rMLST). Prophages were identified using PHASTER and Prophage Hunter.

Results

A total of 293 genomes representing six Acb species exhibited CRISPR-related sequences. These genomes originate from various sources, including clinical specimens, animals, medical devices, and environmental samples. Sequence typing identified 145 ribosomal multilocus sequence types (rSTs). CRISPR-Cas systems were confirmed in 26.3% of the genomes, classified as subtypes I-Fa, I-Fb and I-Fv. Probable CRISPR arrays and cas genes associated with CRISPR-Cas subtypes III-A, I-B, and III-B were also detected. Some of the CRISPR-Cas systems are associated with genomic regions related to Cap4 proteins, and toxin-antitoxin systems. Moreover, prophage sequences were prevalent in 68.9% of the genomes. Analysis revealed a connection between these prophages and CRISPR-Cas systems, indicating an ongoing arms race between the bacteria and their bacteriophages. Furthermore, proteins associated with anti-CRISPR systems, such as AcrF11 and AcrF7, were identified in the A. baumannii and A. pittii genomes.

Discussion

This study elucidates CRISPR-Cas systems and defense mechanisms within the Acb complex, highlighting their diverse distribution and interactions with prophages and other genetic elements. This study also provides valuable insights into the evolution and adaptation of these microorganisms in various environments and clinical settings.

Capsules and their traits shape phage susceptibility and plasmid conjugation efficiency

Bacterial evolution is affected by mobile genetic elements like phages and conjugative plasmids, offering new adaptive traits while incurring fitness costs. Their infection is affected by the bacterial capsule. Yet, its importance has been difficult to quantify because of the high diversity of confounding mechanisms in bacterial genomes such as anti-viral systems and surface receptor modifications. Swapping capsule loci between Klebsiella pneumoniae strains allowed us to quantify their impact on plasmid and phage infection independently of genetic background. Capsule swaps systematically invert phage susceptibility, revealing serotypes as key determinants of phage infection. Capsule types also influence conjugation efficiency in both donor and recipient cells, a mechanism shaped by capsule volume and conjugative pilus structure. Comparative genomics confirmed that more permissive serotypes in the lab correspond to the strains acquiring more conjugative plasmids in nature. The least capsule-sensitive pili (F-like) are the most frequent in the species' plasmids, and are the only ones associated with both antibiotic resistance and virulence factors, driving the convergence between virulence and antibiotics resistance in the population. These results show how traits of cellular envelopes define slow and fast lanes of infection by mobile genetic elements, with implications for population dynamics and horizontal gene transfer.

Biocontrol of Escherichia coli O157:H7 by Enterobacter asburiae AEB30 on intact cantaloupe melons

Escherichia coli O157:H7 causes >73,000 foodborne illnesses in the United States annually, many of which have been associated with fresh ready-to-eat produce including cantaloupe melons. In this study, we created a produce-associated bacterial (PAB) library containing >7500 isolates and screened them for the ability to inhibit the growth of E. coli O157:H7 using an in vitro fluorescence-based growth assay. One isolate, identified by 16S and whole-genome sequence analysis as Enterobacter asburiae, was able to inhibit the growth of E. coli by ~30-fold in vitro and produced zones of inhibition between 13 and 21 mm against 12 E. coli outbreak strains in an agar spot assay. We demonstrated that E. asburiae AEB30 was able to grow, persist and inhibit the growth of E. coli on cantaloupe melons under simulated pre- and post-harvest conditions. Analysis of the E. asburiae AEB30 genome revealed an operon encoding a contact-dependent growth inhibition (CDI) system that when mutated resulted in the loss of E. coli growth inhibition. These data suggest that E. asburiae AEB30 is a potential biocontrol agent to prevent E. coli contamination of cantaloupe melons in both pre- and post-harvest environments and that its mode of action is via a CDI system.

Genomic analysis of Ralstonia pickettii reveals the genetic features for potential pathogenicity and adaptive evolution in drinking water

Ralstonia pickettii, the most critical clinical pathogen of the genus Ralstonia, has been identified as a causative agent of numerous harmful infections. Additionally, Ralstonia pickettii demonstrates adaptability to extreme environmental conditions, such as those found in drinking water. In this study, we conducted a comprehensive genomic analysis to investigate the genomic characteristics related to potential pathogenicity and adaptive evolution in drinking water environments of Ralstonia pickettii. Through phylogenetic analysis and population genetic analysis, we divided Ralstonia pickettii into five Groups, two of which were associated with drinking water environments. The open pan-genome with a large and flexible gene repertoire indicated a high genetic plasticity. Significant differences in functional enrichment were observed between the core- and pan-genome of different groups. Diverse mobile genetic elements (MGEs), extensive genomic rearrangements, and horizontal gene transfer (HGT) events played a crucial role in generating genetic diversity. In drinking water environments, Ralstonia pickettii exhibited strong adaptability, and the acquisition of specific adaptive genes was potentially facilitated by genomic islands (GIs) and HGT. Furthermore, environmental pressures drove the adaptive evolution of Ralstonia pickettii, leading to the accumulation of unique mutations in key genes. These mutations may have a significant impact on various physiological functions, particularly carbon metabolism and energy metabolism. The presence of virulence-related elements associated with macromolecular secretion systems, virulence factors, and antimicrobial resistance indicated the potential pathogenicity of Ralstonia pickettii, making it capable of causing multiple nosocomial infections. This study provides comprehensive insights into the potential pathogenicity and adaptive evolution of Ralstonia pickettii in drinking water environments from a genomic perspective.

Genomic insights into Spiroplasma endosymbionts that induce male-killing and protective phenotypes in the pea aphid

The endosymbiotic bacteria Spiroplasma (Mollicutes) infect diverse plants and arthropods, and some of which induce male killing, where male hosts are killed during development. Male-killing Spiroplasma strains belong to either the phylogenetically distant Citri-Poulsonii or Ixodetis groups. In Drosophila flies, Spiroplasma poulsonii induces male killing via the Spaid toxin. While Spiroplasma ixodetis infects a wide range of insects and arachnids, little is known about the genetic basis of S. ixodetis-induced male killing. Here, we analyzed the genome of S. ixodetis strains in the pea aphid Acyrthosiphon pisum (Aphididae, Hemiptera). Genome sequencing constructed a complete genome of a male-killing strain, sAp269, consisting of a 1.5 Mb circular chromosome and an 80 Kb plasmid. sAp269 encoded putative virulence factors containing either ankyrin repeat, ovarian tumor-like deubiquitinase, or ribosome inactivating protein domains, but lacked the Spaid toxin. Further comparative genomics of Spiroplasma strains in A. pisum biotypes adapted to different host plants revealed their phylogenetic associations and the diversity of putative virulence factors. Although the mechanisms of S. ixodetis-induced male killing in pea aphids remain elusive, this study underlines the dynamic genome evolution of S. ixodetis and proposes independent acquisition events of male-killing mechanisms in insects.

Xanthomonas protegens sp. nov., a novel rice seed-associated bacterium, provides in vivo protection against X. oryzae pv. oryzae, the bacterial leaf blight pathogen

Historically, Xanthomonas species are primarily known for their pathogenicity against plants, but recently, there have been more findings of non-pathogenic xanthomonads. In the present study, we report isolates from healthy rice seeds that belong to a new species, Xanthomonas protegens, a protector of the rice plants against a serious pathogenic counterpart, i.e. X. oryzae pv. oryzae upon leaf clip co-inoculation. The new member species is non-pathogenic to rice and lacks a type III secretion system. The pangenome investigation revealed a large number of unique genes, including a novel lipopolysaccharide biosynthetic gene cluster, that might be important in its adaptation. The phylo-taxonogenomic analysis revealed that X. protegens is a taxonomic outlier species of X. sontii, a core, vertically transmitted rice seed endophyte with numerous probiotic properties. Interestingly, X. sontii is also reported as a keystone species of healthy rice seed microbiome. The findings and resources will help in the development of unique gene markers and evolutionary studies of X. sontii as a successful symbiont and X. oryzae as a serious pathogen. Here, we propose X. protegens sp. nov. as a novel species of the genus Xanthomonas with PPL118 = MTCC 13396 = CFBP 9164 = ICMP 25181 as the type strain. PPL117, PPL124, PPL125, and PPL126 are other strains of the species.

Mobilizable plasmids drive the spread of antimicrobial resistance genes and virulence genes in Klebsiella pneumoniae

BACKGROUND

Klebsiella pneumoniae is a notorious clinical pathogen and frequently carries various plasmids, which are the main carriers of antimicrobial resistance and virulence genes. In comparison to self-transmissible conjugative plasmids, mobilizable plasmids have received much less attention due to their defects in conjugative elements. However, the contribution of mobilizable plasmids to the horizontal transfer of antimicrobial resistance genes and virulence genes of K. pneumoniae remains unclear. In this study, the transfer, stability, and cargo genes of the mobilizable plasmids of K. pneumoniae were examined via genetic experiments and genomic analysis.

METHODS

Carbapenem-resistant (CR) plasmid pHSKP2 and multidrug-resistant (MDR) plasmid pHSKP3 of K. pneumoniae HS11286, virulence plasmid pRJF293 of K. pneumoniae RJF293 were employed in conjugation assays to assess the transfer ability of mobilizable plasmids. Mimic mobilizable plasmids and genetically modified plasmids were constructed to confirm the cotransfer models. The plasmid morphology was evaluated through XbaI and S1 nuclease pulsed-field gel electrophoresis and/or complete genome sequencing. Mobilizable plasmid stability in transconjugants was analyzed via serial passage culture. In addition, in silico genome analysis of 3923 plasmids of 1194 completely sequenced K. pneumoniae was performed to investigate the distribution of the conjugative elements, the cargo genes, and the targets of the CRISPR-Cas system. The mobilizable MDR plasmid and virulence plasmid of K. pneumoniae were investigated, which carry oriT but lack other conjugative elements.

RESULTS

Our results showed that mobilizable MDR and virulence plasmids carrying oriT but lacking the relaxase gene were able to cotransfer with a helper conjugative CR plasmid across various Klebsiella and Escherichia coli strains. The transfer and stability of mobilizable plasmids rather than conjugative plasmids were not interfered with by the CRISPR-Cas system of recipient strains. According to the in silico analysis, the mobilizable plasmids carry about twenty percent of acquired antimicrobial resistance genes and more than seventy-five percent of virulence genes in K. pneumoniae.

CONCLUSIONS

Our work observed that a mobilizable MDR or virulence plasmid that carries oriT but lacks the relaxase genes transferred with the helper CR conjugative plasmid and mobilizable plasmids escaped from CRISPR-Cas defence and remained stable in recipients. These results highlight the threats of mobilizable plasmids as vital vehicles in the dissemination of antibiotic resistance and virulence genes in K. pneumoniae.

A systematic approach to classify and characterize genomic islands driven by conjugative mobility using protein signatures

Genomic islands (GIs) play a crucial role in the spread of antibiotic resistance, virulence factors and antiviral defense systems in a broad range of bacterial species. However, the characterization and classification of GIs are challenging due to their relatively small size and considerable genetic diversity. Predicting their intercellular mobility is of utmost importance in the context of the emerging crisis of multidrug resistance. Here, we propose a large-scale classification method to categorize GIs according to their mobility profile and, subsequently, analyze their gene cargo. We based our classification decision scheme on a collection of mobility protein motif definitions available in publicly accessible databases. Our results show that the size distribution of GI classes correlates with their respective structure and complexity. Self-transmissible GIs are usually the largest, except in Bacillota and Actinomycetota, accumulate antibiotic and phage resistance genes, and favour the use of a tyrosine recombinase to insert into a host's replicon. Non-mobilizable GIs tend to use a DDE transposase instead. Finally, although tRNA genes are more frequently targeted as insertion sites by GIs encoding a tyrosine recombinase, most GIs insert in a protein-encoding gene. This study is a stepping stone toward a better characterization of mobile GIs in bacterial genomes and their mechanism of mobility.

Pseudomonas forestsoilum sp. nov. and P. tohonis biocontrol bacterial wilt by quenching 3-hydroxypalmitic acid methyl ester

Bacterial wilt caused by Ralstonia solanacearum ranks the second top important bacterial plant disease worldwide. It is also the most important bacterial disease threatening the healthy development of Casuarina equisetifolia protection forest. 3-hydroxypalmitic acid methyl ester (3-OH PAME) functions as an important quorum sensing (QS) signal regulating the expression of virulence genes in R. solanacearum, and has been regarded as an ideal target for disease prevention and control. To screen native microorganisms capable of degrading 3-OH PAME, samples of C. equisetifolia branches and forest soil were collected and cultured in the medium containing 3-OH PAME as the sole carbon source. Bacteria with over 85% degradation rates of 3-OH PAME after 7-day incubation were further separated and purified. As a result, strain Q1-7 isolated from forest soil and strain Q4-3 isolated from C. equisetifolia branches were obtained and identified as Pseudomonas novel species Pseudomonas forestsoilum sp. nov. and P. tohonis, respectively, according to whole genome sequencing results. The degradation efficiencies of 3-OH PAME of strains Q1-7 and Q4-3 were 95.80% and 100.00% at 48 h, respectively. Both strains showed high esterase activities and inhibited R. solanacearum exopolysaccharide (EPS) and cellulase production. Application of strains Q1-7 and Q4-3 effectively protects C. equisetifolia, peanut and tomato plants from infection by R. solanacearum. Findings in this study provide potential resources for the prevention and control of bacterial wilt caused by R. solanacearum, as well as valuable materials for the identification of downstream quenching genes and the research and development of quenching enzymes for disease control.

Distinct genomic routes underlie transitions to specialised symbiotic lifestyles in deep-sea annelid worms

Bacterial symbioses allow annelids to colonise extreme ecological niches, such as hydrothermal vents and whale falls. Yet, the genetic principles sustaining these symbioses remain unclear. Here, we show that different genomic adaptations underpin the symbioses of phylogenetically related annelids with distinct nutritional strategies. Genome compaction and extensive gene losses distinguish the heterotrophic symbiosis of the bone-eating worm Osedax frankpressi from the chemoautotrophic symbiosis of deep-sea Vestimentifera. Osedax's endosymbionts complement many of the host's metabolic deficiencies, including the loss of pathways to recycle nitrogen and synthesise some amino acids. Osedax's endosymbionts possess the glyoxylate cycle, which could allow more efficient catabolism of bone-derived nutrients and the production of carbohydrates from fatty acids. Unlike in most Vestimentifera, innate immunity genes are reduced in O. frankpressi, which, however, has an expansion of matrix metalloproteases to digest collagen. Our study supports that distinct nutritional interactions influence host genome evolution differently in highly specialised symbioses.

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.

Dynamics of rice microbiomes reveal core vertically transmitted seed endophytes

BACKGROUND

Plants and their associated microbiota constitute an assemblage of species known as holobionts. The plant seed microbiome plays an important role in nutrient uptake and stress attenuation. However, the core vertically transmitted endophytes remain largely unexplored.

RESULTS

To gain valuable insights into the vertical transmission of rice seed core endophytes, we conducted a large-scale analysis of the microbiomes of two generations of six different rice varieties from five microhabitats (bulk soil, rhizosphere, root, stem, and seed) from four geographic locations. We showed that the microhabitat rather than the geographic location and rice variety was the primary driver of the rice microbiome assemblage. The diversity and network complexity of the rice-associated microbiome decreased steadily from far to near the roots, rice exterior to interior, and from belowground to aboveground niches. Remarkably, the microbiomes of the roots, stems, and seeds of the rice interior compartments were not greatly influenced by the external environment. The core bacterial endophytes of rice were primarily comprised of 14 amplicon sequence variants (ASVs), 10 of which, especially ASV_2 (Pantoea) and ASV_48 (Xanthomonas), were identified as potentially vertically transmitted taxa because they existed across generations, were rarely present in exterior rice microhabitats, and were frequently isolated from rice seeds. The genome sequences of Pantoea and Xanthomonas isolated from the parental and offspring seeds showed a high degree of average nucleotide and core protein identity, indicating vertical transmission of seed endophytes across generations. In silico prediction indicated that the seed endophytes Pantoea and Xanthomonas possessed streamlined genomes with short lengths, low-complexity metabolism, and various plant growth-promoting traits. We also found that all strains of Pantoea and Xanthomonas exhibited cellulase activity and produced indole-3-acetic acid. However, most strains exhibited insignificant antagonism to the major pathogens of rice, such as Magnaporthe oryzae and X. oryzae pv. oryzae.

CONCLUSION

Overall, our study revealed that microhabitats, rather than site-specific environmental factors or host varieties, shape the rice microbiome. We discovered the vertically transmitted profiles and keystone taxa of the rice microbiome, which led to the isolation of culturable seed endophytes and investigation of their potential roles in plant-microbiome interactions. Our results provide insights on vertically transmitted microbiota and suggest new avenues for improving plant fitness via the manipulation of seed-associated microbiomes.  Video Abstract.

Comparative Pathogenomic Analysis of Two Banana Pathogenic Dickeya Strains Isolated from China and the Philippines

Dickeya is a major and typical member of soft rot Pectobacteriaceae (SRP) with a wide range of plant hosts worldwide. Previous studies have identified D. zeae as the causal agent of banana soft rot disease in China. In 2017, we obtained banana soft rot pathogen strain FZ06 from the Philippines. Genome sequencing and analysis indicated that FZ06 can be classified as D. dadantii and represents a novel subspecies of D. dadantii, which we propose to name as subsp. paradisiaca. Compared with Chinese banana soft rot pathogenic strain D. zeae MS2, strain FZ06 has a similar host range but different virulence; FZ06 is significantly less virulent to banana and potato but more virulent to Chinese cabbage and onion. Characterization of virulence factors revealed obviously less production of pectate lyases (Pels), polygalacturonases (Pehs), proteases (Prts), and extrapolysaccharides (EPSs), as well as lower swimming and swarming motility and biofilm formation in strain FZ06. Genomic comparison of the two strains revealed five extra gene clusters in FZ06, including one Stt-type T2SS, three T4SSs, and one T4P. Expression of cell wall degrading enzyme (CWDE)-encoding genes is significantly lower in FZ06 than in MS2.

Assessing the genomic composition, putative ecological relevance and biotechnological potential of plasmids from sponge bacterial symbionts

Plasmid-mediated transfer of genes can have direct consequences in several biological processes within sponge microbial communities. However, very few studies have attempted genomic and functional characterization of plasmids from marine host-associated microbial communities in general and those of sponges in particular. In the present study, we used an endogenous plasmid isolation method to obtain plasmids from bacterial symbionts of the marine sponges Stylissa carteri and Paratetilla sp. and investigated the genomic composition, putative ecological relevance and biotechnological potential of these plasmids. In total, we isolated and characterized three complete plasmids, three plasmid prophages and one incomplete plasmid. Our results highlight the importance of plasmids to transfer relevant genetic traits putatively involved in microbial symbiont adaptation and host-microbe and microbe-microbe interactions. For example, putative genes involved in bacterial response to chemical stress, competition, metabolic versatility and mediation of bacterial colonization and pathogenicity were detected. Genes coding for enzymes and toxins of biotechnological potential were also detected. Most plasmid prophage coding sequences were, however, hypothetical proteins with unknown functions. Overall, this study highlights the ecological relevance of plasmids in the marine sponge microbiome and provides evidence that plasmids of sponge bacterial symbionts may represent an untapped resource of genes of biotechnological interest.

Diverse Genomic Traits Differentiate Sinking-Particle-Associated versus Free-Living Microbes throughout the Oligotrophic Open Ocean Water Column

Bacteria and archaea are central to the production, consumption, and remineralization of dissolved and particulate organic matter and contribute critically to carbon delivery, nutrient availability, and energy transformations in the deep ocean. To explore environmentally relevant genomic traits of sinking-particle-associated versus free-living microbes, we compared habitat-specific metagenome-assembled genomes recovered throughout the water column in the North Pacific Subtropical Gyre. The genomic traits of sinking-particle-associated versus free-living prokaryotes were compositionally, functionally, and phylogenetically distinct. Substrate-specific transporters and extracellular peptidases and carbohydrate-active enzymes were more enriched and diverse in particle-associated microbes at all depths than in free-living counterparts. These data indicate specific roles for particle-attached microbes in particle substrate hydrolysis, uptake, and remineralization. Shallow-water particle-associated microbes had elevated genomic GC content and proteome nitrogen content and reduced proteome carbon content in comparison to abyssal particle-associated microbes. An inverse trend was observed for their sympatric free-living counterparts. These different properties of attached microbes are postulated to arise in part due to elevated organic and inorganic nitrogen availability inside sinking particles. Particle-attached microbes also were enriched in genes for environmental sensing via two-component regulatory systems, and cell-cell interactions via extracellular secretion systems, reflecting their surface-adapted lifestyles. Finally, particle-attached bacteria had greater predicted maximal growth efficiencies than free-living bacterioplankton at all depths. All of these particle-associated specific genomic and proteomic features appear to be driven by microhabitat-specific elevated nutrient and energy availability as well as surface-associated competitive and synergistic ecological interactions. Although some of these characteristics have been previously postulated or observed individually, we report them together here in aggregate via direct comparisons of cooccurring free-living and sinking-particle-attached microbial genomes from the open ocean.

Biosynthetic potential of the global ocean microbiome

Natural microbial communities are phylogenetically and metabolically diverse. In addition to underexplored organismal groups¹, this diversity encompasses a rich discovery potential for ecologically and biotechnologically relevant enzymes and biochemical compounds^2,3. However, studying this diversity to identify genomic pathways for the synthesis of such compounds⁴ and assigning them to their respective hosts remains challenging. The biosynthetic potential of microorganisms in the open ocean remains largely uncharted owing to limitations in the analysis of genome-resolved data at the global scale. Here we investigated the diversity and novelty of biosynthetic gene clusters in the ocean by integrating around 10,000 microbial genomes from cultivated and single cells with more than 25,000 newly reconstructed draft genomes from more than 1,000 seawater samples. These efforts revealed approximately 40,000 putative mostly new biosynthetic gene clusters, several of which were found in previously unsuspected phylogenetic groups. Among these groups, we identified a lineage rich in biosynthetic gene clusters (‘Candidatus Eudoremicrobiaceae’) that belongs to an uncultivated bacterial phylum and includes some of the most biosynthetically diverse microorganisms in this environment. From these, we characterized the phospeptin and pythonamide pathways, revealing cases of unusual bioactive compound structure and enzymology, respectively. Together, this research demonstrates how microbiomics-driven strategies can enable the investigation of previously undescribed enzymes and natural products in underexplored microbial groups and environments.

Global ocean microbiome survey reveals the bacterial family ‘Candidatus Eudoremicrobiaceae’, which includes some of the most biosynthetically diverse microorganisms in the ocean environment.

Comprehensive Analysis Reveals the Genetic and Pathogenic Diversity of Ralstonia solanacearum Species Complex and Benefits Its Taxonomic Classification

Ralstonia solanacearum species complex (RSSC) is a diverse group of plant pathogens that attack a wide range of hosts and cause devastating losses worldwide. In this study, we conducted a comprehensive analysis of 131 RSSC strains to detect their genetic diversity, pathogenicity, and evolution dynamics. Average nucleotide identity analysis was performed to explore the genomic relatedness among these strains, and finally obtained an open pangenome with 32,961 gene families. To better understand the diverse evolution and pathogenicity, we also conducted a series of analyses of virulence factors (VFs) and horizontal gene transfer (HGT) in the pangenome and at the single genome level. The distribution of VFs and mobile genetic elements (MGEs) showed significant differences among different groups and strains, which were consistent with the new nomenclatures of the RSSC with three distinct species. Further functional analysis showed that most HGT events conferred from Burkholderiales and played a great role in shaping the genomic plasticity and genetic diversity of RSSC genomes. Our work provides insights into the genetic polymorphism, evolution dynamics, and pathogenetic variety of RSSC and provides strong supports for the new taxonomic classification, as well as abundant resources for studying host specificity and pathogen emergence.

Comparative Genomic Analysis of Vibrio cincinnatiensis Provides Insights into Genetic Diversity, Evolutionary Dynamics, and Pathogenic Traits of the Species

Vibrio cincinnatiensis is a poorly understood pathogenic Vibrio species, and the underlying mechanisms of its genetic diversity, genomic plasticity, evolutionary dynamics, and pathogenicity have not yet been comprehensively investigated. Here, a comparative genomic analysis of V. cincinnatiensis was constructed. The open pan-genome with a flexible gene repertoire exhibited genetic diversity. The genomic plasticity and stability were characterized by the determinations of diverse mobile genetic elements (MGEs) and barriers to horizontal gene transfer (HGT), respectively. Evolutionary divergences were exhibited by the difference in functional enrichment and selective pressure between the different components of the pan-genome. The evolution on the Chr I and Chr II core genomes was mainly driven by purifying selection. Predicted essential genes in V. cincinnatiensis were mainly found in the core gene families on Chr I and were subject to stronger evolutionary constraints. We identified diverse virulence-related elements, including the gene clusters involved in encoding flagella, secretion systems, several pili, and scattered virulence genes. Our results indicated the pathogenic potential of V. cincinnatiensis and highlighted that HGT events from other Vibrio species promoted pathogenicity. This pan-genome study provides comprehensive insights into this poorly understood species from the genomic perspective.

Pangenome Evolution Reconciles Robustness and Instability of Rhizobial Symbiosis

Root nodulating rhizobia are nearly ubiquitous in soils and provide the critical service of nitrogen fixation to thousands of legume species, including staple crops. However, the magnitude of fixed nitrogen provided to hosts varies markedly among rhizobia strains, despite host legumes having mechanisms to selectively reward beneficial strains and to punish ones that do not fix sufficient nitrogen. Variation in the services of microbial mutualists is considered paradoxical given host mechanisms to select beneficial genotypes. Moreover, the recurrent evolution of non-fixing symbiont genotypes is predicted to destabilize symbiosis, but breakdown has rarely been observed. Here, we deconstructed hundreds of genome sequences from genotypically and phenotypically diverse Bradyrhizobium strains and revealed mechanisms that generate variation in symbiotic nitrogen fixation. We show that this trait is conferred by a modular system consisting of many extremely large integrative conjugative elements and few conjugative plasmids. Their transmissibility and propensity to reshuffle genes generate new combinations that lead to uncooperative genotypes and make individual partnerships unstable. We also demonstrate that these same properties extend beneficial associations to diverse host species and transfer symbiotic capacity among diverse strains. Hence, symbiotic nitrogen fixation is underpinned by modularity, which engenders flexibility, a feature that reconciles evolutionary robustness and instability. These results provide new insights into mechanisms driving the evolution of mobile genetic elements. Moreover, they yield a new predictive model on the evolution of rhizobial symbioses, one that informs on the health of organisms and ecosystems that are hosts to symbionts and that helps resolve the long-standing paradox.

Long-Term Stability of Bacterial Associations in a Microcosm of Ostreococcus tauri (Chlorophyta, Mamiellophyceae)

Phytoplankton-bacteria interactions rule over carbon fixation in the sunlit ocean, yet only a handful of phytoplanktonic-bacteria interactions have been experimentally characterized. In this study, we investigated the effect of three bacterial strains isolated from a long-term microcosm experiment with one Ostreococcus strain (Chlorophyta, Mamiellophyceae). We provided evidence that two Roseovarius strains (Alphaproteobacteria) had a beneficial effect on the long-term survival of the microalgae whereas one Winogradskyella strain (Flavobacteriia) led to the collapse of the microalga culture. Co-cultivation of the beneficial and the antagonistic strains also led to the loss of the microalga cells. Metagenomic analysis of the microcosm is consistent with vitamin B12 synthesis by the Roseovarius strains and unveiled two additional species affiliated to Balneola (Balneolia) and Muricauda (Flavobacteriia), which represent less than 4% of the reads, whereas Roseovarius and Winogradskyella recruit 57 and 39% of the reads, respectively. These results suggest that the low-frequency bacterial species may antagonize the algicidal effect of Winogradskyella in the microbiome of Ostreococcus tauri and thus stabilize the microalga persistence in the microcosm. Altogether, these results open novel perspectives into long-term stability of phytoplankton cultures.

Pan-Genome Analysis of Delftia tsuruhatensis Reveals Important Traits Concerning the Genetic Diversity, Pathogenicity, and Biotechnological Properties of the Species

Delftia tsuruhatensis strains have long been known to promote plant growth and biological control. Recently, it has become an emerging opportunistic pathogen in humans. However, the genomic characteristics of the genetic diversity, pathogenicity, and biotechnological properties have not yet been comprehensively investigated. Here, a comparative pan-genome analysis was constructed. The open pan-genome with a large and flexible gene repertoire exhibited a high degree of genetic diversity. The purifying selection was the main force to drive pan-genome evolution. Significant differences were observed in the evolutionary relationship, functional enrichment, and degree of selective pressure between the different components of the pan-genome. A high degree of genetic plasticity was characterized by the determinations of diverse mobile genetic elements (MGEs), massive genomic rearrangement, and horizontal genes. Horizontal gene transfer (HGT) plays an important role in the genetic diversity of this bacterium and the formation of genomic traits. Our results revealed the occurrence of diverse virulence-related elements associated with macromolecular secretion systems, virulence factors associated with multiple nosocomial infections, and antimicrobial resistance, indicating the pathogenic potential. Lateral flagellum, T1SS, T2SS, T6SS, Tad pilus, type IV pilus, and a part of virulence-related genes exhibited general properties, whereas polar flagellum, T4SS, a part of virulence-related genes, and resistance genes presented heterogeneous properties. The pan-genome also harbors abundant genetic traits related to secondary metabolism, carbohydrate active enzymes (CAZymes), and phosphate transporter, indicating rhizosphere adaptation, plant growth promotion, and great potential uses in agriculture and biological control. This study provides comprehensive insights into this uncommon species from the genomic perspective.

IMPORTANCED. tsuruhatensis is considered a plant growth-promoting rhizobacterium (PGPR), an organic pollutant degradation strain, and an emerging opportunistic pathogen to the human. However, the genetic diversity, the evolutionary dynamics, and the genetic basis of these remarkable traits are still little known. We constructed a pan-genome analysis for D. tsuruhatensis and revealed extensive genetic diversity and genetic plasticity exhibited by open pan-genome, diverse mobile genetic elements (MGEs), genomic rearrangement, and horizontal genes. Our results highlight that horizontal gene transfer (HGT) and purifying selection are important forces in D. tsuruhatensis genetic evolution. The abundant virulence-related elements associated with macromolecular secretion systems, virulence factors, and antimicrobial resistance could contribute to the pathogenicity of this bacterium. Therefore, clinical microbiologists need to be aware of D. tsuruhatensis as an opportunistic pathogen. The genetic profiles of secondary metabolism, carbohydrate active enzymes (CAZymes), and phosphate transporter could provide insight into the genetic armory of potential applications for agriculture and biological control of D. tsuruhatensis in general.

Sulfur cycling and host-virus interactions in Aquificales-dominated biofilms from Yellowstone's hottest ecosystems

Modern linkages among magmatic, geochemical, and geobiological processes provide clues about the importance of thermophiles in the origin of biogeochemical cycles. The aim of this study was to identify the primary chemoautotrophs and host-virus interactions involved in microbial colonization and biogeochemical cycling at sublacustrine, vapor-dominated vents that represent the hottest measured ecosystems in Yellowstone National Park (~140 °C). Filamentous microbial communities exposed to extreme thermal and geochemical gradients were sampled using a remotely operated vehicle and subjected to random metagenome sequencing and microscopic analyses. Sulfurihydrogenibium (phylum Aquificae) was the predominant lineage (up to 84% relative abundance) detected at vents that discharged high levels of dissolved H₂, H₂S, and CO₂. Metabolic analyses indicated carbon fixation by Sulfurihydrogenibium spp. was powered by the oxidation of reduced sulfur and H₂, which provides organic carbon for heterotrophic community members. Highly variable Sulfurihydrogenibium genomes suggested the importance of intra-population diversity under extreme environmental and viral pressures. Numerous lytic viruses (primarily unclassified taxa) were associated with diverse archaea and bacteria in the vent community. Five circular dsDNA uncultivated virus genomes (UViGs) of ~40 kbp length were linked to the Sulfurihydrogenibium metagenome-assembled genome (MAG) by CRISPR spacer matches. Four UViGs contained consistent genome architecture and formed a monophyletic cluster with the recently proposed Pyrovirus genus within the Caudovirales. Sulfurihydrogenibium spp. also contained CRISPR arrays linked to plasmid DNA with genes for a novel type IV filament system and a highly expressed β-barrel porin. A diverse suite of transcribed secretion systems was consistent with direct microscopic analyses, which revealed an extensive extracellular matrix likely critical to community structure and function. We hypothesize these attributes are fundamental to the establishment and survival of microbial communities in highly turbulent, extreme-gradient environments.

Characterization, Pathogenicity, Phylogeny, and Comparative Genomic Analysis of Pseudomonas tolaasii Strains Isolated from Various Mushrooms in China

Since 2016, devastating bacterial blotch affecting the fruiting bodies of Agaricus bisporus, Cordyceps militaris, Flammulina filiformis, and Pleurotus ostreatus in China has caused severe economic losses. We isolated 102 bacterial strains and characterized them polyphasically. We identified the causal agent as Pseudomonas tolaasii and confirmed the pathogenicity of the strains. A host range test further confirmed the pathogen's ability to infect multiple hosts. This is the first report in China of bacterial blotch in C. militaris caused by P. tolaasii. Whole-genome sequences were generated for three strains: Pt11 (6.48 Mb), Pt51 (6.63 Mb), and Pt53 (6.80 Mb), and pangenome analysis was performed with 13 other publicly accessible P. tolaasii genomes to determine their genetic diversity, virulence, antibiotic resistance, and mobile genetic elements. The pangenome of P. tolaasii is open, and many more gene families are likely to emerge with further genome sequencing. Multilocus sequence analysis using the sequences of four common housekeeping genes (glns, gyrB, rpoB, and rpoD) showed high genetic variability among the P. tolaasii strains, with 115 strains clustered into a monophyletic group. The P. tolaasii strains possess various genes for secretion systems, virulence factors, carbohydrate-active enzymes, toxins, secondary metabolites, and antimicrobial resistance genes that are associated with pathogenesis and adapted to different environments. The myriad of insertion sequences, integrons, prophages, and genome islands encoded in the strains may contribute to genome plasticity, virulence, and antibiotic resistance. These findings advance understanding of the determinants of virulence, which can be targeted for the effective control of bacterial blotch disease.

Microbial Diversity Analysis and Genome Sequencing Identify Xanthomonas perforans as the Pathogen of Bacterial Leaf Canker of Water Spinach (Ipomoea aquatic)

Ipomoea aquatica is a leafy vegetable widely cultivated in tropical Asia, Africa, and Oceania. Bacterial leaf canker disease has been attacking the planting fields and seriously affecting the quality of I. aquatica in epidemic areas in China. This study examined the microbial composition of I. aquatica leaves with classical symptoms of spot disease. The results showed that Xanthomonas was overwhelmingly dominant in all four diseased leaf samples but rarely present in rhizospheric soil or irrigation water samples. In addition, Pantoea was also detected in two of the diseased leaf samples. Pathogen isolation, identification, and inoculation revealed that both Xanthomonas sp. TC2-1 and P. ananatis were pathogenic to the leaves of I. aquatic, causing crater-shaped ulcerative spots and yellowing with big brown rot lesions on leaves, respectively. We further sequenced the whole genome of strain TC2-1 and showed that it is a member of X. perforans. Overall, this study identified X. perforans as the causal pathogen of I. aquatica bacterial leaf canker, and P. ananatis as a companion pathogen causing yellowing and brown rot on leaves. The correct identification of the pathogens will provide important basis for future efforts to formulate targeted application strategy for bacterial disease control.

Interplay between the cell envelope and mobile genetic elements shapes gene flow in populations of the nosocomial pathogen Klebsiella pneumoniae

Mobile genetic elements (MGEs) drive genetic transfers between bacteria using mechanisms that require a physical interaction with the cellular envelope. In the high-priority multidrug-resistant nosocomial pathogens (ESKAPE), the first point of contact between the cell and virions or conjugative pili is the capsule. While the capsule can be a barrier to MGEs, it also evolves rapidly by horizontal gene transfer (HGT). Here, we aim at understanding this apparent contradiction by studying the covariation between the repertoire of capsule genes and MGEs in approximately 4,000 genomes of Klebsiella pneumoniae (Kpn). We show that capsules drive phage-mediated gene flow between closely related serotypes. Such serotype-specific phage predation also explains the frequent inactivation of capsule genes, observed in more than 3% of the genomes. Inactivation is strongly epistatic, recapitulating the capsule biosynthetic pathway. We show that conjugative plasmids are acquired at higher rates in natural isolates lacking a functional capsular locus and confirmed experimentally this result in capsule mutants. This suggests that capsule inactivation by phage pressure facilitates its subsequent reacquisition by conjugation. Accordingly, capsule reacquisition leaves long recombination tracts around the capsular locus. The loss and regain process rewires gene flow toward other lineages whenever it leads to serotype swaps. Such changes happen preferentially between chemically related serotypes, hinting that the fitness of serotype-swapped strains depends on the host genetic background. These results enlighten the bases of trade-offs between the evolution of virulence and multidrug resistance and caution that some alternatives to antibiotics by selecting for capsule inactivation may facilitate the acquisition of antibiotic resistance genes (ARGs).

A study of how the complex interaction between capsules and mobile genetic elements shapes gene flow in populations of Klebsiella pneumoniae reveals that capsule inactivation by phage pressure facilitates its subsequent re-acquisition by conjugation, and this loss and re-gain process influences the gene flow towards other lineages whenever it leads to serotype changes.

Pangenomics reveals alternative environmental lifestyles among chlamydiae

Chlamydiae are highly successful strictly intracellular bacteria associated with diverse eukaryotic hosts. Here we analyzed metagenome-assembled genomes of the "Genomes from Earth's Microbiomes" initiative from diverse environmental samples, which almost double the known phylogenetic diversity of the phylum and facilitate a highly resolved view at the chlamydial pangenome. Chlamydiae are defined by a relatively large core genome indicative of an intracellular lifestyle, and a highly dynamic accessory genome of environmental lineages. We observe chlamydial lineages that encode enzymes of the reductive tricarboxylic acid cycle and for light-driven ATP synthesis. We show a widespread potential for anaerobic energy generation through pyruvate fermentation or the arginine deiminase pathway, and we add lineages capable of molecular hydrogen production. Genome-informed analysis of environmental distribution revealed lineage-specific niches and a high abundance of chlamydiae in some habitats. Together, our data provide an extended perspective of the variability of chlamydial biology and the ecology of this phylum of intracellular microbes.

Genomic and transcriptomic insights into the habitat adaptation of the diazotrophic paddy-field cyanobacterium Nostoc sphaeroides

Nitrogen-fixing cyanobacteria are common in paddy fields, one of the most productive wetland ecosystems. Here, we present the complete genome of Nostoc sphaeroides, a paddy-field diazotroph used for food and medicine for more than 1700 years and deciphered the transcriptional regulation during the developmental transition from hormogonia to vegetative filaments with heterocysts. The genome of N. sphaeroides consists of one circular chromosome (6.48 Mb), one of the largest ever reported megaplasmids (2.34 Mb), and seven plasmids. Multiple gene families involved in the adaption to high solar radiation and water fluctuation conditions were found expanded, while genes involved in anoxic adaptation and phosphonate utilization are located on the megaplasmid, suggesting its indispensable role in environmental adaptation. Distinct gene expression patterns were observed during the light-intensity-regulated transition from hormogonia to vegetative filaments, specifically, genes encoding proteins involved in photosynthetic light reaction, carbon fixation, nitrogen metabolism and heterocyst differentiation were significantly upregulated, whereas genes related to cell motility were down-regulated. Our results provide genomic and transcriptomic insights into the adaptation of a filamentous nitrogen-fixing cyanobacterium to the highly dynamic paddy-field habitat, suggesting N. sphaeroides as an excellent system to understand the transition from aquatic to terrestrial habitats and to support sustainable rice production.

GspD, The Type II Secretion System Secretin of Leptospira, Protects Hamsters against Lethal Infection with a Virulent L. interrogans Isolate

The wide variety of pathogenic Leptospira serovars and the weak protection offered by the available vaccines encourage the search for protective immunogens against leptospirosis. We found that the secretin GspD of the type II secretion system (T2S) of Leptospira interrogans serovar Canicola was highly conserved amongst pathogenic serovars and was expressed in vivo during infection, as shown by immunohistochemistry. Convalescent sera of hamsters, dogs, and cows showed the presence of IgG antibodies, recognizing a recombinant version of this protein expressed in Escherichia coli (rGspDLC) in Western blot assays. In a pilot vaccination study, a group of eight hamsters was immunized on days zero and 14 with 50 µg of rGspDLC mixed with Freund’s incomplete adjuvant (FIA). On day 28 of the study, 1,000 LD50 (Lethal Dose 50%) of a virulent strain of Leptospira interrogans serovar Canicola (LOCaS46) were inoculated by an intraoral submucosal route (IOSM). Seventy-five percent protection against disease (p = 0.017573, Fisher’s exact test) and 50% protection against infection were observed in this group of vaccinated hamsters. In contrast, 85% of non-vaccinated hamsters died six to nine days after the challenge. These results suggest the potential usefulness of the T2S secretin GspD of Leptospira as a protective recombinant vaccine against leptospirosis.

Roles of the Hcp family proteins in the pathogenicity of Salmonella typhimurium 14028s

The type VI secretion system (T6SS) is a new secretion system that is widely distributed among Gram-negative bacteria. The core component hemolysin-coregulated protein (Hcp) can be used as both its structural protein and secretory protein or chaperone protein. Studies on Hcp are important to elucidate the overall virulence mechanism of T6SS. Salmonella typhimurium is an important foodborne pathogen. There are three copies of hcp genes identified in S. Typhimurium 14028s. This study aimed to characterize the functions of the three Hcp family proteins and to elucidate the interactions among them. The hcp gene deletion mutants were constructed by λ Red-based recombination system. Effects of hcp mutation on the pathogenicity of 14028s were studied by bacterial competition assays, Dictyostelium discoideum assays and mouse model. The three Hcp family proteins were found to play different roles. Hcp1 can affect the transcription of rpoS and type 2 flagellar gene and influence the motility of 14028s. It is also involved in the intracellular survival of 14028s in Dictyostelium discoideum; Hcp2 is involved in the early proliferative capacity of 14028s in mice and can prevent its excessive proliferation; Hcp3 did not show direct functions in these assays. Hcp1 can interact with Hcp2 and Hcp3. Deletion of one hcp gene can result in a transcription level variation in the other two hcp genes. Our findings elucidated the functions of the three Hcp family proteins in S.Typhimurium and illustrated that there are interactions between different Hcp proteins. This study will be helpful to fully understand how T6SS actions in an organism.

Comparative Genomic Analysis Reveals Genetic Mechanisms of the Variety of Pathogenicity, Antibiotic Resistance, and Environmental Adaptation of Providencia Genus

The bacterial genus Providencia is Gram-negative opportunistic pathogens, which have been isolated from a variety of environments and organisms, ranging from humans to animals. Providencia alcalifaciens, Providencia rettgeri, and Providencia stuartii are the most common clinical isolates, however, these three species differ in their pathogenicity, antibiotic resistance and environmental adaptation. Genomes of 91 isolates of the genus Providencia were investigated to clarify their genetic diversity, focusing on virulence factors, antibiotic resistance genes, and environmental adaptation genes. Our study revealed an open pan-genome for the genus Providencia containing 14,720 gene families. Species of the genus Providencia exhibited different functional constraints, with the core genes, accessory genes, and unique genes. A maximum-likelihood phylogeny reconstructed with concatenated single-copy core genes classified all Providencia isolates into 11 distant groups. Comprehensive and systematic comparative genomic analyses revealed that specific distributions of virulence genes, which were highly homologous to virulence genes of the genus Proteus, contributed to diversity in pathogenicity of Providencia alcalifaciens, Providencia rettgeri, and Providencia stuartii. Furthermore, multidrug resistance (MDR) phenotypes of isolates of Providencia rettgeri and Providencia stuartii were predominantly due to resistance genes from class 1 and 2 integrons. In addition, Providencia rettgeri and Providencia stuartii harbored more genes related to material transport and energy metabolism, which conferred a stronger ability to adapt to diverse environments. Overall, our study provided valuable insights into the genetic diversity and functional features of the genus Providencia, and revealed genetic mechanisms underlying diversity in pathogenicity, antibiotic resistance and environmental adaptation of members of this genus.

Horizontal Gene Transfer Clarifies Taxonomic Confusion and Promotes the Genetic Diversity and Pathogenicity of Plesiomonas shigelloides

The taxonomic position of P. shigelloides has been the subject of debate for a long time, and until now, the evolutionary dynamics and pathogenesis of P. shigelloides were unclear. In this study, pan-genome analysis indicated extensive genetic diversity and the presence of large and variable gene repertoires. Our results revealed that horizontal gene transfer was the focal driving force for the genetic diversity of the P. shigelloides pan-genome and might have contributed to the emergence of novel properties. Vibrionaceae and Aeromonadaceae were found to be the predominant donor taxa for horizontal genes, which might have caused the taxonomic confusion historically. Comparative genomic analysis revealed the potential of P. shigelloides to cause intestinal and invasive diseases. Our results could advance the understanding of the evolution and pathogenesis of P. shigelloides, particularly in elucidating the role of horizontal gene transfer and investigating virulence-related elements.

Plesiomonas shigelloides is an emerging pathogen that has been shown to be involved in gastrointestinal diseases and extraintestinal infections in humans. However, the taxonomic position, evolutionary dynamics, and pathogenesis of P. shigelloides remain unclear. We reported the draft genome sequences of 12 P. shigelloides strains representing different serogroups. We were able to determine a clear distinction between P. shigelloides and other members of Enterobacterales via core genome phylogeny, Neighbor-Net network, and average genome identity analysis. The pan-genome analysis of P. shigelloides revealed extensive genetic diversity and presented large flexible gene repertoires, while the core genome phylogeny exhibited a low level of clonality. The discordance between the core genome phylogeny and the pan-genome phylogeny indicated that flexible accessory genomes account for an important proportion of the evolution of P. shigelloides, which was subsequently characterized by determinations of hundreds of horizontally transferred genes (horizontal genes), massive gene expansions and contractions, and diverse mobile genetic elements (MGEs). The apparently high levels of horizontal gene transfer (HGT) in P. shigelloides were conferred from bacteria with novel properties from other taxa (mainly Vibrionaceae and Aeromonadaceae), which caused the historical taxonomic confusion and shaped the virulence gene pools. Furthermore, P. shigelloides genomes contain many macromolecular secretion system genes, virulence factor genes, and resistance genes, indicating its potential to cause intestinal and invasive infections. Collectively, our work provides insights into the phylogenetic position, evolutionary dynamic, and pathogenesis of P. shigelloides at the genomic level, which could facilitate the observation and research of this important pathogen.

IMPORTANCE The taxonomic position of P. shigelloides has been the subject of debate for a long time, and until now, the evolutionary dynamics and pathogenesis of P. shigelloides were unclear. In this study, pan-genome analysis indicated extensive genetic diversity and the presence of large and variable gene repertoires. Our results revealed that horizontal gene transfer was the focal driving force for the genetic diversity of the P. shigelloides pan-genome and might have contributed to the emergence of novel properties. Vibrionaceae and Aeromonadaceae were found to be the predominant donor taxa for horizontal genes, which might have caused the taxonomic confusion historically. Comparative genomic analysis revealed the potential of P. shigelloides to cause intestinal and invasive diseases. Our results could advance the understanding of the evolution and pathogenesis of P. shigelloides, particularly in elucidating the role of horizontal gene transfer and investigating virulence-related elements.

Chromosome Architecture and Gene Content of the Emergent Pathogen Acinetobacter haemolyticus

Acinetobacter haemolyticus is a Gammaproteobacterium that has been involved in serious diseases frequently linked to the nosocomial environment. Most of the strains causing such infections are sensitive to a wide variety of antibiotics, but recent reports indicate that this pathogen is acquiring very efficiently carbapenem-resistance determinants like the blaNDM-1 gene, all over the world. With this work we contribute with a collection set of 31 newly sequenced nosocomial A. haemolyticus isolates. Genome analysis of these sequences and others collected from RefSeq indicates that their chromosomes are organized in 12 syntenic blocks that contain most of the core genome genes. These blocks are separated by hypervariable regions that are rich in unique gene families, but also have signals of horizontal gene transfer. Genes involved in virulence or encoding different secretion systems are located inside syntenic regions and have recombination signals. The relative order of the synthetic blocks along the A. haemolyticus chromosome can change, indicating that they have been subject to several kinds of inversions. Genomes of this microorganism show large differences in gene content even if they are in the same clade. Here we also show that A. haemolyticus has an open pan-genome.

Repeated horizontal gene transfers triggered parallel evolution of magnetotaxis in two evolutionary divergent lineages of magnetotactic bacteria

Under the same selection pressures, two genetically divergent populations may evolve in parallel toward the same adaptive solutions. Here, we hypothesized that magnetotaxis (i.e., magnetically guided chemotaxis) represents a key adaptation to micro-oxic habitats in aquatic sediments and that its parallel evolution homogenized the phenotypes of two evolutionary divergent clusters of freshwater spirilla. All magnetotactic bacteria affiliated to the Magnetospirillum genus (Alphaproteobacteria class) biomineralize the same magnetic particle chains and share highly similar physiological and ultrastructural features. We looked for the processes that could have contributed at shaping such an evolutionary pattern by reconciling species and gene trees using newly sequenced genomes of Magnetospirillum related bacteria. We showed that repeated horizontal gene transfers and homologous recombination of entire operons contributed to the parallel evolution of magnetotaxis. We propose that such processes could represent a more parsimonious and rapid solution for adaptation compared with independent and repeated de novo mutations, especially in the case of traits as complex as magnetotaxis involving tens of interacting proteins. Besides strengthening the idea about the importance of such a function in micro-oxic habitats, these results reinforce previous observations in experimental evolution suggesting that gene flow could alleviate clonal interference and speed up adaptation under some circumstances.

Structure and Evolution of Acinetobacter baumannii Plasmids

Acinetobacter baumannii is an emergent bacterial pathogen that provokes many types of infections in hospitals around the world. The genome of this organism consists of a chromosome and plasmids. These plasmids vary over a wide size range and many of them have been linked to the acquisition of antibiotic-resistance genes. Our bioinformatic analyses indicate that A. baumannii plasmids belong to a small number of plasmid lineages. The general structure of these lineages seems to be very stable and consists not only of genes involved in plasmid maintenance functions but of gene sets encoding poorly characterized proteins, not obviously linked to survival in the hospital setting, and opening the possibility that they improve the parasitic properties of plasmids. An analysis of genes involved in replication, suggests that members of the same plasmid lineage are part of the same plasmid incompatibility group. The same analysis showed the necessity of classifying the Rep proteins in ten new groups, under the scheme proposed by Bertini et al. (2010). Also, we show that some plasmid lineages have the potential capacity to replicate in many bacterial genera including those embracing human pathogen species, while others seem to replicate only within the limits of the Acinetobacter genus. Moreover, some plasmid lineages are widely distributed along the A. baumannii phylogenetic tree. Despite this, a number of them lack genes involved in conjugation or mobilization functions. Interestingly, only 34.6% of the plasmids analyzed here possess antibiotic resistance genes and most of them belong to fourteen plasmid lineages of the twenty one described here. Gene flux between plasmid lineages appears primarily limited to transposable elements, which sometimes carry antibiotic resistance genes. In most plasmid lineages transposable elements and antibiotic resistance genes are secondary acquisitions. Finally, broad host-range plasmids appear to have played a crucial role.

MicroScope: an integrated platform for the annotation and exploration of microbial gene functions through genomic, pangenomic and metabolic comparative analysis

Large-scale genome sequencing and the increasingly massive use of high-throughput approaches produce a vast amount of new information that completely transforms our understanding of thousands of microbial species. However, despite the development of powerful bioinformatics approaches, full interpretation of the content of these genomes remains a difficult task. Launched in 2005, the MicroScope platform (https://www.genoscope.cns.fr/agc/microscope) has been under continuous development and provides analysis for prokaryotic genome projects together with metabolic network reconstruction and post-genomic experiments allowing users to improve the understanding of gene functions. Here we present new improvements of the MicroScope user interface for genome selection, navigation and expert gene annotation. Automatic functional annotation procedures of the platform have also been updated and we added several new tools for the functional annotation of genes and genomic regions. We finally focus on new tools and pipeline developed to perform comparative analyses on hundreds of genomes based on pangenome graphs. To date, MicroScope contains data for >11 800 microbial genomes, part of which are manually curated and maintained by microbiologists (>4500 personal accounts in September 2019). The platform enables collaborative work in a rich comparative genomic context and improves community-based curation efforts.

Identification of integrative and conjugative elements in pathogenic and commensal Neisseriaceae species via genomic distributions of DNA uptake sequence dialects

Mobile genetic elements (MGEs) are key factors responsible for dissemination of virulence determinants and antimicrobial-resistance genes amongst pathogenic bacteria. Conjugative MGEs are notable for their high gene loads donated per transfer event, broad host ranges and phylogenetic ubiquity amongst prokaryotes, with the subclass of chromosomally inserted integrative and conjugative elements (ICEs) being particularly abundant. The focus on a small number of model systems has biased the study of ICEs towards those conferring readily selectable phenotypes to host cells, whereas the identification and characterization of integrated cryptic elements remains challenging. Even though antimicrobial resistance and horizontally acquired virulence genes are major factors aggravating neisserial infection, conjugative MGEs of Neisseria gonorrhoeae and Neisseria meningitidis remain poorly characterized. Using a phenotype-independent approach based on atypical distributions of DNA uptake sequences (DUSs) in MGEs relative to the chromosomal background, we have identified two groups of chromosomally integrated conjugative elements in Neisseria: one found almost exclusively in pathogenic species possibly deriving from the genus Kingella, the other belonging to a group of Neisseria mucosa-like commensals. The former element appears to enable transfer of traditionally gonococcal-specific loci such as the virulence-associated toxin-antitoxin system fitAB to N. meningitidis chromosomes, whilst the circular form of the latter possesses a unique attachment site (attP) sequence seemingly adapted to exploit DUS motifs as chromosomal integration sites. In addition to validating the use of DUS distributions in Neisseriaceae MGE identification, the >170 identified ICE sequences provide a valuable resource for future studies of ICE evolution and host adaptation.

Identifying Conjugative Plasmids and Integrative Conjugative Elements with CONJscan

We present a computational method to identify conjugative systems in plasmids and chromosomes using the CONJscan module of MacSyFinder. The method relies on the identification of the protein components of the system using hidden Markov model profiles and then checking that the composition and genetic organization of the system is consistent with that expected from a conjugative system. The method can be assessed online using the Galaxy workflow or locally using a standalone software. The latter version allows to modify the models of the module (i.e., to change the expected components, their number, and their organization).CONJscan identifies conjugative systems, but when the mobile genetic element is integrative (ICE), one often also wants to delimit it from the chromosome. We present a method, with a script, to use the results of CONJscan and comparative genomics to delimit ICE in chromosomes. The method provides a visual representation of the ICE location. Together, these methods facilitate the identification of conjugative elements in bacterial genomes.

Comparative Genomic Analysis of Citrobacter and Key Genes Essential for the Pathogenicity of Citrobacter koseri

Citrobacter species are opportunistic bacterial pathogens that have been implicated in both nosocomial and community-acquired infections. Among the genus Citrobacter, Citrobacter koseri is often isolated from clinical material, and has been known to cause meningitis and brain abscess in neonates and immunocompromised individuals. The virulence determinants of Citrobacter, however, remain largely unknown. Based on traditional methods, the genus Citrobacter has been divided into 11 species, but this has been problematic. Here, we determined an improved, detailed, and more accurate phylogeny of the genus Citrobacter based on whole genome sequence (WGS) data from 129 Citrobacter genomes, 31 of which were sequenced in this study. A maximum likelihood (ML) phylogeny constructed with core genome single-nucleotide polymorphisms (SNPs) classified all Citrobacter isolates into 11 distinct groups, with all C. koseri strains clustering into a single group. For comprehensive and systematic comparative genomic analyses, we investigated the distribution of virulence factors, resistance genes, and macromolecular secretion systems among the Citrobacter genus. Moreover, combined with group-specific genes analysis, we identified a key gene cluster for iron transport, which is present in the C. koseri group, but absent in other the groups, suggesting that the high-pathogenicity island (HPI) cluster may be important for the pathogenicity of C. koseri. Animal experiments showed that loss of the HPI cluster significantly decreased C. koseri virulence in mice and rat. Further, we provide evidence to explain why Citrobacter freundii is less susceptible than C. koseri to several antibiotics in silico. Overall, our data reveal novel virulence clusters specific to the predominantly pathogenic C. koseri strains, which form the basis for elucidating the virulence mechanisms underlying these important pathogens.

Ecological and functional capabilities of an uncultured Kordia sp

Cultivable bacteria represent only a fraction of the diversity in microbial communities. However, the official procedures for classification and characterization of a novel prokaryotic species still rely on isolates. Nevertheless, due to single cell genomics, it is possible to retrieve genomes from environmental samples by sequencing them individually, and to assign specific genes to a specific taxon, regardless of their ability to grow in culture. In this study, a complete description was performed for uncultured Kordia sp. TARA_039_SRF, a proposed novel species within the genus Kordia, using culture-independent techniques. The type material was a high-quality draft genome (94.97% complete, 4.65% gene redundancy) co-assembled using ten nearly identical single amplified genomes (SAGs) from surface seawater in the North Indian Ocean during the Tara Oceans Expedition. The assembly process was optimized to obtain the best possible assembly metrics and a less fragmented genome. The closest relative of the species was Kordia periserrulae, which shared 97.56% similarity of the 16S rRNA gene, 75% orthologs and 89.13% average nucleotide identity. The functional potential of the proposed novel species included proteorhodopsin, the ability to incorporate nitrate, cytochrome oxidases with high affinity for oxygen, and CAZymes that were unique features within the genus. Its abundance at different depths and size fractions was also evaluated together with its functional annotation, revealing that its putative ecological niche could be particles of phytoplanktonic origin. It could putatively attach to these particles and consume them while sinking to the deeper and oxygen depleted layers of the North Indian Ocean.

A Mesophilic Aeromonas salmonicida Strain Isolated from an Unsuspected Host, the Migratory Bird Pied Avocet

Aeromonas salmonicida is a Gram-negative bacterium, known as a fish pathogen since its discovery. Although the species was initially considered psychrophilic, a mesophilic subspecies (pectinolytica) and many other mesophilic strains still not attributed to subspecies have been described in the last two decades. These mesophilic strains were sampled from various sources, including humans, and some of them are known to be pathogenic. In this study, we describe a strain, JF2480, which was isolated from the spleen, and also found the kidney and liver of a dead pied avocet (Recurvirostra avosetta), a type of migratory bird inhabiting aquatic environments. A core genome phylogenomic analysis suggests that JF2480 is taxonomically distant from other known A. salmonicida subspecies. The genome sequence confirms that the strain possesses key virulence genes that are present in the typical A. salmonicida psychrophilic subspecies, with the exception of the genes encoding the type three secretion system (T3SS). Bacterial virulence assays conducted on the surrogate host Dictyostelium discoideum amoeba confirmed that the strain is virulent despite the lack of T3SS. Bacterial growth curves showed that strain JF2480 grow well at 40 °C, the body temperature of the pied avocet, and even faster at 41 °C, compared to other mesophilic strains. Discovery of this strain further demonstrates the extent of the phylogenomic tree of this species. This study also suggests that A. salmonicida can infect a wider array of hosts than previously suspected and that we need to rethink the way we perceive A. salmonicida's natural environment.

Comparative genomic analysis of the Hafnia genus reveals an explicit evolutionary relationship between the species alvei and paralvei and provides insights into pathogenicity

BACKGROUND

The Hafnia genus is an opportunistic pathogen that has been implicated in both nosocomial and community-acquired infections. Although Hafnia is fairly often isolated from clinical material, its taxonomy has remained an unsolved riddle, and the involvement and importance of Hafnia in human disease is also uncertain. Here, we used comparative genomic analysis to define the taxonomy of Hafnia, identify species-specific genes that may be the result of ecological and pathogenic specialization, and reveal virulence-related genetic profiles that may contribute to pathogenesis.

RESULTS

One complete genome sequence and 19 draft genome sequences for Hafnia strains were generated and combined with 27 publicly available genomes. We provided high-resolution typing methods by constructing phylogeny and population structure based on single-copy core genes in combination with whole genome average nucleotide identity to identify two distant Hafnia species (alvei and paralvei) and one mislabeled strain. The open pan-genome and the presence of numerous mobile genetic elements reveal that Hafnia has undergone massive gene rearrangements. Presence of species-specific core genomes associated with metabolism and transport suggests the putative niche differentiation between alvei and paralvei. We also identified possession of diverse virulence-related profiles in both Hafnia species., including the macromolecular secretion system, virulence, and antimicrobial resistance. In the macromolecular system, T1SS, Flagellum 1, Tad pilus and T6SS-1 were conserved in Hafnia, whereas T4SS, T5SS, and other T6SSs exhibited the evolution of diversity. The virulence factors in Hafnia are related to adherence, toxin, iron uptake, stress adaptation, and efflux pump. The identified resistance genes are associated with aminoglycoside, beta-lactam, bacitracin, cationic antimicrobial peptide, fluoroquinolone, and rifampin. These virulence-related profiles identified at the genomic level provide insights into Hafnia pathogenesis and the differentiation between alvei and paralvei.

CONCLUSIONS

Our research using core genome phylogeny and comparative genomics analysis of a larger collection of strains provides a comprehensive view of the taxonomy and species-specific traits between Hafnia species. Deciphering the genome of Hafnia strains possessing a reservoir of macromolecular secretion systems, virulence factors, and resistance genes related to pathogenicity may provide insights into addressing its numerous infections and devising strategies to combat the pathogen.

Isolates from Colonic Spirochetosis in Humans Show High Genomic Divergence and Potential Pathogenic Features but Are Not Detected Using Standard Primers for the Human Microbiota

This is the first report of whole-genome analysis of clinical isolates from individuals with colonic spirochetosis. This characterization provides new opportunities in understanding the physiology and potentials of these bacteria that densely colonize the gut in the individuals infected. The observation that standard 16S amplicon primers fail to detect colonic spirochetosis may have major implications for studies searching for associations between members of the microbiota and clinical conditions such as irritable bowel syndrome (IBS) and should be taken into consideration in project design and interpretation of gastrointestinal tract microbiota in population-based and clinical settings.

Colonic spirochetosis, diagnosed based on the striking appearance in histological sections, still has an obscure clinical relevance, and only a few bacterial isolates from this condition have been characterized to date. In a randomized, population-based study in Stockholm, Sweden, 745 healthy individuals underwent colonoscopy with biopsy sampling. Of these individuals, 17 (2.3%) had colonic spirochetosis, which was associated with eosinophilic infiltration and a 3-fold-increased risk for irritable bowel syndrome (IBS). We aimed to culture the bacteria and perform whole-genome sequencing of the isolates from this unique representative population sample. From 14 out of 17 individuals with spirochetosis we successfully isolated, cultured, and performed whole-genome sequencing of in total 17 isolates, including the Brachyspira aalborgi type strain, 513A. Also, 16S analysis of the mucosa-associated microbiota was performed in the cases and nonspirochetosis controls. We found one isolate to be of the species Brachyspira pilosicoli; all remaining isolates were of the species Brachyspira aalborgi. Besides displaying extensive genetic heterogeneity, the isolates harbored several mucin-degrading enzymes and other virulence-associated genes that could confer a pathogenic potential in the human colon. We also showed that 16S amplicon sequencing using standard primers for human microbiota studies failed to detect Brachyspira due to primer incompatibility.

IMPORTANCE This is the first report of whole-genome analysis of clinical isolates from individuals with colonic spirochetosis. This characterization provides new opportunities in understanding the physiology and potentials of these bacteria that densely colonize the gut in the individuals infected. The observation that standard 16S amplicon primers fail to detect colonic spirochetosis may have major implications for studies searching for associations between members of the microbiota and clinical conditions such as irritable bowel syndrome (IBS) and should be taken into consideration in project design and interpretation of gastrointestinal tract microbiota in population-based and clinical settings.