Uncovering the boundaries of Campylobacter species through large-scale phylogenetic and nucleotide identity analyses

Campylobacter species are typically helical shaped, Gram-negative, and non-spore-forming bacteria. Species in this genus include established foodborne and animal pathogens as well as emerging pathogens. The accumulation of genomic data from the Campylobacter genus has increased exponentially in recent years, accompanied by the discovery of putative new species. At present, the lack of a standardized species boundary complicates distinguishing established and novel species. We defined the Campylobacter genus core genome (500 loci) using publicly available Campylobacter complete genomes (n = 498) and constructed a core genome phylogeny using 2,193 publicly available Campylobacter genomes to examine inter-species diversity and species boundaries. Utilizing 8,440 Campylobacter genomes representing 33 species and 8 subspecies, we found species delineation based on an average nucleotide identity (ANI) cutoff of 94.2% is consistent with the core genome phylogeny. We identified 60 ANI genomic species that delineated Campylobacter species in concordance with previous comparative genetic studies. All pairwise ANI genomic species pairs had in silico DNA-DNA hybridization scores of less than 70%, supporting their delineation as separate species. We provide the tool Campylobacter Genomic Species typer (CampyGStyper) that assigns ANI genomic species to query genomes based on ANI similarities to medoid genomes from each ANI genomic species with an accuracy of 99.96%. The ANI genomic species definitions proposed here allow consistent species definition in the Campylobacter genus and will facilitate the detection of novel species in the future.IMPORTANCEIn recent years, Campylobacter has gained recognition as the leading cause of bacterial gastroenteritis worldwide, leading to a substantial rise in the collection of genomic data of the Campylobacter genus in public databases. Currently, a standardized Campylobacter species boundary at the genomic level is absent, leading to challenges in detecting emerging pathogens and defining putative novel species within this genus. We used a comprehensive representation of genomes of the Campylobacter genus to construct a core genome phylogenetic tree. Furthermore, we found an average nucleotide identity (ANI) of 94.2% as the optimal cutoff to define the Campylobacter species. Using this cutoff, we identified 60 ANI genomic species which provided a standardized species definition and nomenclature. Importantly, we have developed Campylobacter Genomic Species typer (CampyGStyper), which can robustly and accurately assign these ANI genomic species to Campylobacter genomes, thereby aiding pathogen surveillance and facilitating evolutionary and epidemiological studies of existing and emerging pathogens in the genus Campylobacter.

A core genome multi-locus sequence typing scheme for Streptococcus uberis: an evolution in typing a genetically diverse pathogen

Streptococcus uberis is a globally endemic and poorly controlled cause of bovine mastitis impacting the sustainability of the modern dairy industry. A core genome was derived from 579 newly sequenced S. uberis isolates, along with 305 publicly available genome sequences of S. uberis isolated from 11 countries around the world and used to develop a core genome multi-locus sequence typing (cgMLST) scheme. The S. uberis core genome comprised 1475 genes, and these were used to identify 1447 curated loci that were indexed into the cgMLST scheme. This was able to type 1012 of 1037 (>97  %) isolates used and differentiated the associated sequences into 932 discrete core genome sequence types (cgSTs). Analysis of the phylogenetic relationships of cgSTs revealed no clear clustering of isolates based on metadata such as disease status or year of isolation. Geographical clustering of cgSTs was limited to identification of a UK-centric clade, but cgSTs from UK isolates were also dispersed with those originating from other geographical regions across the entire phylogenetic topology. The cgMLST scheme offers a new tool for the detailed analysis of this globally important pathogen of dairy cattle. Initial analysis has re-emphasized and exemplified the genetically diverse nature of the global population of this opportunistic pathogen.

In-Depth Genome Characterization and Pan-Genome Analysis of Strain KMM 296, a Producer of Highly Active Alkaline Phosphatase; Proposal for the Reclassification of Cobetia litoralis and Cobetia pacifica as the Later Heterotypic Synonyms of Cobetia amphilecti and Cobetia marina, and Emended Description of the Species Cobetia amphilecti and Cobetia marina

A strictly aerobic, Gram-stain-negative, rod-shaped, and motile bacterium, designated strain KMM 296, isolated from the coelomic fluid of the mussel Crenomytilus grayanus, was investigated in detail due to its ability to produce a highly active alkaline phosphatase CmAP of the structural family PhoA. A previous taxonomic study allocated the strain to the species Cobetia marina, a member of the family Halomonadaceae of the class Gammaproteobacteria. However, 16S rRNA gene sequencing showed KMM 296's relatedness to Cobetia amphilecti NRIC 0815T. The isolate grew with 0.5-19% NaCl at 4-42 °C and hydrolyzed Tweens 20 and 40 and L-tyrosine. The DNA G+C content was 62.5 mol%. The prevalent fatty acids were C18:1 ω7c, C12:0 3-OH, C18:1 ω7c, C12:0, and C17:0 cyclo. The polar lipid profile was characterized by the presence of phosphatidylethanolamine, phosphatidylglycerol, phosphatidic acid, and also an unidentified aminolipid, phospholipid, and a few unidentified lipids. The major respiratory quinone was Q-8. According to phylogenomic and chemotaxonomic evidence, and the nearest neighbors, the strain KMM 296 represents a member of the species C. amphilecti. The genome-based analysis of C. amphilecti NRIC 0815T and C. litoralis NRIC 0814T showed their belonging to a single species. In addition, the high similarity between the C. pacifica NRIC 0813T and C. marina LMG 2217T genomes suggests their affiliation to one species. Based on the rules of priority, C. litoralis should be reclassified as a later heterotypic synonym of C. amphilecti, and C. pacifica is a later heterotypic synonym of C. marina. The emended descriptions of the species C. amphilecti and C. marina are also proposed.

Nationwide molecular epidemiology of carbapenemase-producing Citrobacter spp. in France in 2019 and 2020

IMPORTANCE

The emergence of carbapenemase producers in Enterobacterales mostly involves Escherichia coli, Klebsiella pneumoniae, and Enterobacter cloacae complex species. However, in France, we observed the emergence and the rapid dissemination of carbapenemase in Citrobacter spp. In this study, we demonstrated that a wide variety of carbapenemases is produced by many different species of Citrobacter spp. However, we clearly identify three high-risk clones of Citrobacter freundii, ST8, ST22, and ST91 that drive the spread of carbapenemase in France. This epidemiological study paves the way of further analysis that would aim to identify the virulence factors involved in this pellicular ability of these three clones to disseminate at the hospital.

Genomic delineation and description of species and within-species lineages in the genus Pantoea

As the name of the genus Pantoea ("of all sorts and sources") suggests, this genus includes bacteria with a wide range of provenances, including plants, animals, soils, components of the water cycle, and humans. Some members of the genus are pathogenic to plants, and some are suspected to be opportunistic human pathogens; while others are used as microbial pesticides or show promise in biotechnological applications. During its taxonomic history, the genus and its species have seen many revisions. However, evolutionary and comparative genomics studies have started to provide a solid foundation for a more stable taxonomy. To move further toward this goal, we have built a 2,509-gene core genome tree of 437 public genome sequences representing the currently known diversity of the genus Pantoea. Clades were evaluated for being evolutionarily and ecologically significant by determining bootstrap support, gene content differences, and recent recombination events. These results were then integrated with genome metadata, published literature, descriptions of named species with standing in nomenclature, and circumscriptions of yet-unnamed species clusters, 15 of which we assigned names under the nascent SeqCode. Finally, genome-based circumscriptions and descriptions of each species and each significant genetic lineage within species were uploaded to the LINbase Web server so that newly sequenced genomes of isolates belonging to any of these groups could be precisely and accurately identified.

Dynamics of the Streptomyces chromosome: chance and necessity

Streptomyces are prolific producers of specialized metabolites with applications in medicine and agriculture. Remarkably, these bacteria possess a large linear chromosome that is genetically compartmentalized: core genes are grouped in the central part, while the ends are populated by poorly conserved genes including antibiotic biosynthetic gene clusters. The genome is highly unstable and exhibits distinct evolutionary rates along the chromosome. Recent chromosome conformation capture (3C) and comparative genomics studies have shed new light on the interplay between genome dynamics in space and time. Here, we review insights that illustrate how the balance between chance (random genome variations) and necessity (structural and functional constraints) may have led to the emergence of spatial structuring of the Streptomyces chromosome.

Exploring environmental intra-species diversity through non-redundant pangenome assemblies

At the genome level, microorganisms are highly adaptable both in terms of allele and gene composition. Such heritable traits emerge in response to different environmental niches and can have a profound influence on microbial community dynamics. As a consequence, any individual genome or population will contain merely a fraction of the total genetic diversity of any operationally defined "species", whose ecological potential can thus be only fully understood by studying all of their genomes and the genes therein. This concept, known as the pangenome, is valuable for studying microbial ecology and evolution, as it partitions genomes into core (present in all the genomes from a species, and responsible for housekeeping and species-level niche adaptation among others) and accessory regions (present only in some, and responsible for intra-species differentiation). Here we present SuperPang, an algorithm producing pangenome assemblies from a set of input genomes of varying quality, including metagenome-assembled genomes (MAGs). SuperPang runs in linear time and its results are complete, non-redundant, preserve gene ordering and contain both coding and non-coding regions. Our approach provides a modular view of the pangenome, identifying operons and genomic islands, and allowing to track their prevalence in different populations. We illustrate this by analysing intra-species diversity in Polynucleobacter, a bacterial genus ubiquitous in freshwater ecosystems, characterized by their streamlined genomes and their ecological versatility. We show how SuperPang facilitates the simultaneous analysis of allelic and gene content variation under different environmental pressures, allowing us to study the drivers of microbial diversification at unprecedented resolution.

Transporter Proteins as Ecological Assets and Features of Microbial Eukaryotic Pangenomes

Here we review two connected themes in evolutionary microbiology: (a) the nature of gene repertoire variation within species groups (pangenomes) and (b) the concept of metabolite transporters as accessory proteins capable of providing niche-defining "bolt-on" phenotypes. We discuss the need for improved sampling and understanding of pangenome variation in eukaryotic microbes. We then review the factors that shape the repertoire of accessory genes within pangenomes. As part of this discussion, we outline how gene duplication is a key factor in both eukaryotic pangenome variation and transporter gene family evolution. We go on to outline how, through functional characterization of transporter-encoding genes, in combination with analyses of how transporter genes are gained and lost from accessory genomes, we can reveal much about the niche range, the ecology, and the evolution of virulence of microbes. We advocate for the coordinated systematic study of eukaryotic pangenomes through genome sequencing and the functional analysis of genes found within the accessory gene repertoire.

Evolutionary and Pan-genome Analysis of Three Important Black-pigmented Periodontal Pathogens

OBJECTIVE

To analyse the pan-genome of three black-pigmented periodontal pathogens: Porphyromonas gingivalis, Prevotella intermedia and Prevotella nigrescens.

METHODS

Pan-genome analyses of 66, 33 and 5 publicly available whole-genome sequences of P. gingivalis, P. intermedia and P. nigrescens, respectively, were performed using Pan-genome Analysis Pipeline software (version 1.2.1; Beijing Institute of Genomics, Chinese Academy of Sciences, Beijing, PR China). Phylogenetic trees were constructed based on the entire pan-genome and single nucleotide polymorphisms within the core genome. The distribution and abundance of virulence genes in the core and dispensable genomes were also compared in the three species.

RESULTS

All three species possess an open pan-genome. The core genome of P. gingivalis, P. intermedia and P. nigrescens included 1001, 1514 and 1745 orthologous groups, respectively, which were mainly related to basic cellular functions such as metabolism. The dispensable genome of P. gingivalis, P. intermedia and P. nigrescens was composed of 2814, 2689 and 906 orthologous groups, respectively, and it was enriched in genes involved in pathogenicity or with unknown functions. Phylogenetic trees presented a clear separation of P. gingivalis, P. intermedia and P. nigrescens, verifying the reclassification of the black-pigmented species. Furthermore, the three species shared almost the same virulence factors involved in adhesion, proteolysis and evasion of host defences. Some of these virulence genes were conserved across species whereas others belonged to the dispensable genome, which might be acquired through horizontal gene transfer.

CONCLUSION

This study highlighted the usefulness of pan-genome analysis to infer evolutionary cues for black-pigmented species, indicating their homology and phylogenomic diversity.

A panoramic view of the genomic landscape of the genus Streptomyces

We delineate the evolutionary plasticity of the ecologically and biotechnologically important genus Streptomyces, by analysing the genomes of 213 species. Streptomycetes genomes demonstrate high levels of internal homology, whereas the genome of their last common ancestor was already complex. Importantly, we identify the species-specific fingerprint proteins that characterize each species. Even among closely related species, we observed high interspecies variability of chromosomal protein-coding genes, species-level core genes, accessory genes and fingerprints. Notably, secondary metabolite biosynthetic gene clusters (smBGCs), carbohydrate-active enzymes (CAZymes) and protein-coding genes bearing the rare TTA codon demonstrate high intraspecies and interspecies variability, which emphasizes the need for strain-specific genomic mining. Highly conserved genes, such as those specifying genus-level core proteins, tend to occur in the central region of the chromosome, whereas those encoding proteins with evolutionarily volatile species-level fingerprints, smBGCs, CAZymes and TTA-codon-bearing genes are often found towards the ends of the linear chromosome. Thus, the chromosomal arms emerge as the part of the genome that is mainly responsible for rapid adaptation at the species and strain level. Finally, we observed a moderate, but statistically significant, correlation between the total number of CAZymes and three categories of smBGCs (siderophores, e-Polylysin and type III lanthipeptides) that are related to competition among bacteria.

Identifying the core genome of the nucleus-forming bacteriophage family and characterization of Erwinia phage RAY

We recently discovered that some bacteriophages establish a nucleus-like replication compartment (phage nucleus), but the core genes that define nucleus-based phage replication and their phylogenetic distribution were still to be determined. Here, we show that phages encoding the major phage nucleus protein chimallin share 72 conserved genes encoded within seven gene blocks. Of these, 21 core genes are unique to nucleus-forming phage, and all but one of these genes encode proteins of unknown function. We propose that these phages comprise a novel viral family we term Chimalliviridae. Fluorescence microscopy and cryoelectron tomography studies of Erwinia phage vB_EamM_RAY confirm that many of the key steps of nucleus-based replication are conserved among diverse chimalliviruses and reveal variations on this replication mechanism. This work expands our understanding of phage nucleus and PhuZ spindle diversity and function, providing a roadmap for identifying key mechanisms underlying nucleus-based phage replication.

Naegleria genus pangenome reveals new structural and functional insights into the versatility of these free-living amoebae

Introduction

Free-living amoebae of the Naegleria genus belong to the major protist clade Heterolobosea and are ubiquitously distributed in soil and freshwater habitats. Of the 47 Naegleria species described, N. fowleri is the only one being pathogenic to humans, causing a rare but fulminant primary amoebic meningoencephalitis. Some Naegleria genome sequences are publicly available, but the genetic basis for Naegleria diversity and ability to thrive in diverse environments (including human brain) remains unclear.

Methods

Herein, we constructed a high-quality Naegleria genus pangenome to obtain a comprehensive catalog of genes encoded by these amoebae. For this, we first sequenced, assembled, and annotated six new Naegleria genomes.

Results and Discussion

Genome architecture analyses revealed that Naegleria may use genome plasticity features such as ploidy/aneuploidy to modulate their behavior in different environments. When comparing 14 near-to-complete genome sequences, our results estimated the theoretical Naegleria pangenome as a closed genome, with 13,943 genes, including 3,563 core and 10,380 accessory genes. The functional annotations revealed that a large fraction of Naegleria genes show significant sequence similarity with those already described in other kingdoms, namely Animalia and Plantae. Comparative analyses highlighted a remarkable genomic heterogeneity, even for closely related strains and demonstrate that Naegleria harbors extensive genome variability, reflected in different metabolic repertoires. If Naegleria core genome was enriched in conserved genes essential for metabolic, regulatory and survival processes, the accessory genome revealed the presence of genes involved in stress response, macromolecule modifications, cell signaling and immune response. Commonly reported N. fowleri virulence-associated genes were present in both core and accessory genomes, suggesting that N. fowleri's ability to infect human brain could be related to its unique species-specific genes (mostly of unknown function) and/or to differential gene expression. The construction of Naegleria first pangenome allowed us to move away from a single reference genome (that does not necessarily represent each species as a whole) and to identify essential and dispensable genes in Naegleria evolution, diversity and biology, paving the way for further genomic and post-genomic studies.

Phylogenomic diversity of Vibrio species and other Gammaproteobacteria isolated from Pacific oysters (Crassostrea gigas) during a summer mortality outbreak

The Pacific oyster (PO), Crassostrea gigas, is an important commercial marine species but periodically experiences large stock losses due to disease events known as summer mortality. Summer mortality has been linked to environmental perturbations and numerous viral and bacterial agents, indicating this disease is multifactorial in nature. In 2013 and 2014, several summer mortality events occurred within the Port Stephens estuary (NSW, Australia). Extensive culture and molecular-based investigations were undertaken and several potentially pathogenic Vibrio species were identified. To improve species identification and genomically characterise isolates obtained from this outbreak, whole-genome sequencing (WGS) and subsequent genomic analyses were performed on 48 bacterial isolates, as well as a further nine isolates from other summer mortality studies using the same batch of juveniles. Average nucleotide identity (ANI) identified most isolates to the species level and included members of the Photobacterium, Pseudoalteromonas, Shewanella and Vibrio genera, with Vibrio species making up more than two-thirds of all species identified. Construction of a phylogenomic tree, ANI analysis, and pan-genome analysis of the 57 isolates represents the most comprehensive culture-based phylogenomic survey of Vibrios during a PO summer mortality event in Australian waters and revealed large genomic diversity in many of the identified species. Our analysis revealed limited and inconsistent associations between isolate species and their geographical origins, or host health status. Together with ANI and pan-genome results, these inconsistencies suggest that to determine the role that microbes may have in Pacific oyster summer mortality events, isolate identification must be at the taxonomic level of strain. Our WGS data (specifically, the accessory genomes) differentiated bacterial strains, and coupled with associated metadata, highlight the possibility of predicting a strain's environmental niche and level of pathogenicity.

Occurrence and Variability of the Efflux Pump Gene norA across the Staphylococcus Genus

NorA is one of the main native MDR efflux pumps of Staphylococcus aureus, contributing to reduced susceptibility towards fluoroquinolones and biocides, but little is known about its variability within S. aureus or its distribution and conservation among other staphylococci. We screened for sequences homologous to S. aureus norA and found it in 61 out of the 63 Staphylococcus species described. To the best of our knowledge, this is the first study to report the occurrence of norA across the Staphylococcus genus. The norA phylogenetic tree follows the evolutionary relations of staphylococci and the closely related Mammalliicoccus genus. Comparative analyses suggest a conservation of the NorA function in staphylococci. We also analyzed the variability of norA within S. aureus, for which there are several circulating norA alleles, differing up to 10% at the nucleotide level, which may hamper proper norA detection. We demonstrate the applicability of a PCR-based algorithm to detect and differentiate norA alleles in 52 S. aureus representing a wider collection of 89 isolates from different hosts. Our results highlight the prevalence of norAI and norAII in different settings and the association of norA alleles with specific S. aureus clonal lineages. Ultimately, it confirms the applicability of our PCR-based algorithm to rapidly detect and assign the different norA alleles, a trait that may impact antimicrobial efflux capacity and the search for potential NorA inhibitors.

Molecular Epidemiology of Staphylococcus aureus in China Reveals the Key Gene Features Involved in Epidemic Transmission and Adaptive Evolution

Staphylococcus aureus is a Gram-positive pathogen that causes various infections in humans and domestic animals. In China, S. aureus is the most common Gram-positive pathogen that causes clinical infections. However, there are few comprehensive genome-based molecular epidemiology studies to investigate the genotypic background of the major S. aureus clones that are epidemic in China. Here, four S. aureus isolates that were recovered from hospital personnel were sequenced. In combination with whole-genome sequencing (WGS) data of 328 S. aureus strains as references, we performed a comprehensive molecular epidemiology study to reveal the molecular epidemic characterization of S. aureus that is epidemic in China. It was found that 332 S. aureus isolates were phylogenetically categorized into 4 major epidemic groups with different epidemiology phenotypes. Each group has exclusive features in virulence genotypic profiles, antimicrobial resistance genotypic profiles, core and pangenome features representing the differences involved in genetic features, evolutionary processes, and potential future evolutionary directions. Moreover, a comparative core genome analysis of 332 S. aureus isolates indicated several key genes that contributed to differences in molecular epidemic characterization and promoted the adaptive evolutionary process of each group. This study provides a comprehensive understanding of molecular epidemiological characteristics and adaptive evolutionary directions of major S. aureus clones that are epidemic in China. IMPORTANCE Staphylococcus aureus is an important Gram-positive pathogen that is epidemic worldwide and causes various infections in humans and domestic animals. However, there has been relatively little research on comprehensive molecular epidemiology in China. In this research, we reconstructed the phylogenetic relationship based on whole-genome data of strains almost all over China, screened for resistance and virulence genes, and took core and pan genome analysis to perform a comprehensive molecular epidemiology study of S. aureus that is epidemic in China. Our results highlight that there are 4 major epidemic groups with different epidemiology phenotypes after phylogenetic categorization with exclusive genetic features in virulence genotypic profiles, antimicrobial-resistance genotypic profiles, and core and pangenome features, and we found key gene features involved in epidemic transmission and adaptive evolution. Our findings are critical in describing molecular characteristic profiles of S. aureus infection, which could update existing preventive measures and take appropriate strategies.

Comprehensive phylogenomics of Methylobacterium reveals four evolutionary distinct groups and underappreciated phyllosphere diversity

Methylobacterium is a group of methylotrophic microbes associated with soil, fresh water, and particularly the phyllosphere, the aerial part of plants that has been well-studied in terms of physiology but whose evolutionary history and taxonomy are unclear. Recent work has suggested that Methylobacterium is much more diverse than thought previously, questioning its status as an ecologically and phylogenetically coherent taxonomic genus. However, taxonomic and evolutionary studies of Methylobacterium have mostly been restricted to model species, often isolated from habitats other than the phyllosphere, and have yet to utilize comprehensive phylogenomic methods to examine gene trees, gene content, or synteny. By analyzing 189 Methylobacterium genomes from a wide range of habitats, including the phyllosphere, we inferred a robust phylogenetic tree while explicitly accounting for the impact of horizontal gene transfer. We showed that Methylobacterium contains four evolutionarily distinct groups of bacteria (namely A, B, C, D), characterized by different genome size, GC content, gene content and genome architecture, revealing the dynamic nature of Methylobacterium genomes. In addition to recovering 59 described species, we identified 45 candidate species, mostly phyllosphere-associated, stressing the significance of plants as a reservoir of Methylobacterium diversity. We inferred an ancient transition from a free-living lifestyle to association with plant roots in Methylobacteriaceae ancestor, followed by phyllosphere association of three of the major groups (A, B, D), whose early branching in Methylobacterium history has been heavily obscured by HGT. Together, our work lays the foundations for a thorough redefinition of Methylobacterium taxonomy, beginning with the abandonment of Methylorubrum.

The Neisseria gonorrhoeae Accessory Genome and Its Association with the Core Genome and Antimicrobial Resistance

The bacterial accessory genome provides the genetic flexibility needed to facilitate environment and host adaptation. In Neisseria gonorrhoeae, known accessory elements include plasmids which can transfer and mediate antimicrobial resistance (AMR); however, chromosomal accessory genes could also play a role in AMR. Here, the gonococcal accessory genome was characterized using gene-by-gene approaches and its association with the core genome and AMR were assessed. The gonococcal accessory gene pool consisted of 247 genes, which were mainly genes located on large mobile genetic elements, phage associated genes, or genes encoding putative secretion systems. Accessory elements showed similar synteny across genomes, indicating either a predisposition for particular genomic locations or ancestral inheritance that are conserved during strain expansion. Significant associations were found between the prevalence of accessory elements and core genome multi-locus sequence types (cgMLST), consistent with a structured gonococcal population despite frequent horizontal gene transfer (HGT). Increased prevalence of putative DNA exchange regulators was significantly associated with AMR, which included a putative secretion system, methyltransferases and a toxin-antitoxin system. Although frequent HGT results in high genetic diversity in the gonococcus, we found that this is mediated by a small gene pool. In fact, a highly organized genome composition was identified with a strong association between the accessory and core genome. Increased prevalence of DNA exchange regulators in antimicrobial resistant isolates suggests that genetic material exchange plays a role in the development or maintenance of AMR. These findings enhance our understanding of gonococcal genome architecture and have important implications for gonococcal population biology. IMPORTANCE The emergence of antimicrobial resistance (AMR) against third generation cephalosporins in Neisseria gonorrhoeae is a major public health concern, as these are antibiotics of last resort for the effective treatment of gonorrhea. Although the resistance mechanisms against this class of antibiotics have not been entirely resolved, resistance against other classes of antibiotics, such as tetracyclines, is known to be mediated through plasmids, which are known gonococcal extra-chromosomal accessory elements. A complete assessment of the chromosomal accessory genome content and its role in AMR has not yet been undertaken. Here, we comprehensively characterize the gonococcal accessory genome to better understand genome architecture as well as the evolution and mechanisms of AMR in this species.

Genomic Characterization of Bacillus safensis Isolated from Mine Tailings in Peru and Evaluation of Its Cyanide-Degrading Enzyme CynD

Understanding the biochemistry and metabolic pathways of cyanide degradation is necessary to improve the efficacy of cyanide bioremediation processes and industrial requirements. We have isolated and sequenced the genome of a cyanide-degrading Bacillus strain from water in contact with mine tailings from Lima, Peru. This strain was classified as Bacillus safensis based on 16S rRNA gene sequencing and core genome analyses and named B. safensis PER-URP-08. We searched for possible cyanide-degradation enzymes in the genome of this strain and identified a putative cyanide dihydratase (CynD) gene similar to a previously characterized CynD from Bacillus pumilus C1. Sequence analysis of CynD from B. safensis and B. pumilus allow us to identify C-terminal residues that differentiate both CynDs. We then cloned, expressed in Escherichia coli, and purified recombinant CynD from B. safensis PER-URP-08 (CynD_PER-URP-08) and showed that in contrast to CynD from B. pumilus C1, this recombinant CynD remains active at up to pH 9. We also showed that oligomerization of CynD_PER-URP-08 decreases as a function of increased pH. Finally, we demonstrated that transcripts of CynD_PER-URP-08 in B. safensis PER-URP-08 are strongly induced in the presence of cyanide. Our results suggest that the use of B. safensis PER-URP-08 and CynD_PER-URP-08 as potential tool for cyanide bioremediation warrants further investigation. IMPORTANCE Despite being of environmental concern around the world due to its toxicity, cyanide continues to be used in many important industrial processes. Thus, searching for cyanide bioremediation methods is a matter of societal concern and must be present on the political agenda of all governments. Here, we report the isolation, genome sequencing and characterization of cyanide degradation capacity of a bacterial strain isolated from an industrial mining site in Peru. We characterize a cyanide dehydratase (CynD) homolog from one of these bacteria, Bacillus safensis PER-URP-08.

Identification of Putative Vaccine and Drug Targets against the Methicillin-Resistant Staphylococcus aureus by Reverse Vaccinology and Subtractive Genomics Approaches

Methicillin-resistant Staphylococcus aureus (MRSA) is an opportunistic pathogen and responsible for causing life-threatening infections. The emergence of hypervirulent and multidrug-resistant (MDR) S. aureus strains led to challenging issues in antibiotic therapy. Consequently, the morbidity and mortality rates caused by S. aureus infections have a substantial impact on health concerns. The current worldwide prevalence of MRSA infections highlights the need for long-lasting preventive measures and strategies. Unfortunately, effective measures are limited. In this study, we focus on the identification of vaccine candidates and drug target proteins against the 16 strains of MRSA using reverse vaccinology and subtractive genomics approaches. Using the reverse vaccinology approach, 4 putative antigenic proteins were identified; among these, PrsA and EssA proteins were found to be more promising vaccine candidates. We applied a molecular docking approach of selected 8 drug target proteins with the drug-like molecules, revealing that the ZINC4235426 as potential drug molecule with favorable interactions with the target active site residues of 5 drug target proteins viz., biotin protein ligase, HPr kinase/phosphorylase, thymidylate kinase, UDP-N-acetylmuramoyl-L-alanyl-D-glutamate-L-lysine ligase, and pantothenate synthetase. Thus, the identified proteins can be used for further rational drug or vaccine design to identify novel therapeutic agents for the treatment of multidrug-resistant staphylococcal infection.

Pathogenicity and Its Implications in Taxonomy: The Brucella and Ochrobactrum Case

The intracellular pathogens of the genus Brucella are phylogenetically close to Ochrobactrum, a diverse group of free-living bacteria with a few species occasionally infecting medically compromised patients. A group of taxonomists recently included all Ochrobactrum organisms in the genus Brucella based on global genome analyses and alleged equivalences with genera such as Mycobacterium. Here, we demonstrate that such equivalencies are incorrect because they overlook the complexities of pathogenicity. By summarizing Brucella and Ochrobactrum divergences in lifestyle, structure, physiology, population, closed versus open pangenomes, genomic traits, and pathogenicity, we show that when they are adequately understood, they are highly relevant in taxonomy and not unidimensional quantitative characters. Thus, the Ochrobactrum and Brucella differences are not limited to their assignments to different “risk-groups”, a biologically (and hence, taxonomically) oversimplified description that, moreover, does not support ignoring the nomen periculosum rule, as proposed. Since the epidemiology, prophylaxis, diagnosis, and treatment are thoroughly unrelated, merging free-living Ochrobactrum organisms with highly pathogenic Brucella organisms brings evident risks for veterinarians, medical doctors, and public health authorities who confront brucellosis, a significant zoonosis worldwide. Therefore, from taxonomical and practical standpoints, the Brucella and Ochrobactrum genera must be maintained apart. Consequently, we urge researchers, culture collections, and databases to keep their canonical nomenclature.

Meta-analysis of the Ralstonia solanacearum species complex (RSSC) based on comparative evolutionary genomics and reverse ecology

Ralstonia solanacearum species complex (RSSC) strains are bacteria that colonize plant xylem tissue and cause vascular wilt diseases. However, individual strains vary in host range, optimal disease temperatures and physiological traits. To increase our understanding of the evolution, diversity and biology of the RSSC, we performed a meta-analysis of 100 representative RSSC genomes. These 100 RSSC genomes contain 4940 genes on average, and a pangenome analysis found that there are 3262 genes in the core genome (~60 % of the mean RSSC genome) with 13 128 genes in the extensive flexible genome. A core genome phylogenetic tree and a whole-genome similarity matrix aligned with the previously named species (R. solanacearum, R. pseudosolanacearum, R. syzygii) and phylotypes (I–IV). These analyses also highlighted a third unrecognized sub-clade of phylotype II. Additionally, we identified differences between phylotypes with respect to gene content and recombination rate, and we delineated population clusters based on the extent of horizontal gene transfer. Multiple analyses indicate that phylotype II is the most diverse phylotype, and it may thus represent the ancestral group of the RSSC. We also used our genome-based framework to test whether the RSSC sequence variant (sequevar) taxonomy is a robust method to define within-species relationships of strains. The sequevar taxonomy is based on alignments of a single conserved gene (egl). Although sequevars in phylotype II describe monophyletic groups, the sequevar system breaks down in the highly recombinogenic phylotype I, which highlights the need for an improved, cost-effective method for genotyping strains in phylotype I. Finally, we enabled quick and precise genome-based identification of newly sequenced RSSC strains by assigning Life Identification Numbers (LINs) to the 100 strains and by circumscribing the RSSC and its sub-groups in the LINbase Web service.

Strain-level fitness in the gut microbiome is an emergent property of glycans and a single metabolite

The human gut microbiota resides within a diverse chemical environment challenging our ability to understand the forces shaping this ecosystem. Here, we reveal that fitness of the Bacteroidales, the dominant order of bacteria in the human gut, is an emergent property of glycans and one specific metabolite, butyrate. Distinct sugars serve as strain-variable fitness switches activating context-dependent inhibitory functions of butyrate. Differential fitness effects of butyrate within the Bacteroides are mediated by species-level variation in Acyl-CoA thioesterase activity and nucleotide polymorphisms regulating an Acyl-CoA transferase. Using in vivo multi-omic profiles, we demonstrate Bacteroides fitness in the human gut is associated together, but not independently, with Acyl-CoA transferase expression and butyrate. Our data reveal that each strain of the Bacteroides exists within a unique fitness landscape based on the interaction of chemical components unpredictable by the effect of each part alone mediated by flexibility in the core genome.

Pan-genome: setting a new standard for high-quality reference genomes

With the release of high-quality reference genomes assembled by long reads from the third-generation sequencing technology, as well as extensive re-sequencing and population genetic analysis, researchers found that a single reference genome does not represent the diversity within a species. The missing sequences on the reference genome result in an incomplete population genetic polymorphism map. The emergence of pan-genome can well repair the deficiency of single reference genome, which include core genome (responsible for basic biological functions and the main phenotypic characteristics within a species) and the variable genome (related to the genetic diversity or biological characteristics). According to the core and variable genome proportion, the types of pan-genomes can be either open or closed. Here, we review the current exploring of pan-genome for a range of species, to discuss the characteristics of pan-genome in various biological groups. The pan-genome of mammals are more likely closed, while the pan-genomes of microbes, angiosperms, and some invertebrates are likely non-closed. It is possible to complete the reference genome and obtain complete variation information through the pan-genomic study, which will contribute to the study of molecular mechanism for genetic diversity and phenotypic evolution.

Whole Genome Sequence-Based Identification of Clostridium estertheticum Complex Strains Supports the Need for Taxonomic Reclassification Within the Species Clostridium estertheticum

Isolates within the Clostridium estertheticum complex (CEC) have routinely been identified through the 16S rRNA sequence, but the high interspecies sequence similarity reduces the resolution necessary for species level identification and often results in ambiguous taxonomic classification. The current study identified CEC isolates from meat juice (MJS) and bovine fecal samples (BFS) and determined the phylogeny of species within the CEC through whole genome sequence (WGS)-based analyses. About 1,054 MJS were screened for CEC using quantitative real-time PCR (qPCR). Strains were isolated from 33 MJS and 34 BFS qPCR-positive samples, respectively. Pan- and core-genome phylogenomics were used to determine the species identity of the isolates. Average nucleotide identity (ANI) and digital DNA-DNA hybridization (dDDH) were used to validate the species identity. The phylogeny of species within the CEC was determined through a combination of these methods. Twenty-eight clostridia strains were isolated from MJS and BFS samples out of which 13 belonged to CEC. At 95% ANI and 70% dDDH thresholds for speciation, six CEC isolates were identified as genomospecies2 (n=3), Clostridium tagluense (n=2) and genomospecies3 (n=1). Lower thresholds of 94% ANI and 58% dDDH were required for the classification of seven CEC isolates into species C. estertheticum and prevent an overlap between species C. estertheticum and Clostridium frigoriphilum. Combination of the two species and abolishment of current subspecies classification within the species C. estertheticum are proposed. These data demonstrate the suitability of phylogenomics to identify CEC isolates and determine the phylogeny within CEC.

Comparative Genome Analysis and Spore Heat Resistance Assay Reveal a New Component to Population Structure and Genome Epidemiology Within Clostridium perfringens Enterotoxin-Carrying Isolates

Clostridium perfringens causes a variety of human and animal enteric diseases including food poisoning, antibiotic-associated diarrhea, and necrotic enteritis. Yet, the reservoirs of enteropathogenic enterotoxin-producing strains remain unknown. We conducted a genomic comparison of 290 strains and a heat resistance phenotyping of 30 C. perfringens strains to elucidate the population structure and ecology of this pathogen. C. perfringens genomes shared a conserved genetic backbone with more than half of the genes of an average genome conserved in >95% of strains. The cpe-carrying isolates were found to share genetic context: the cpe-carrying plasmids had different distribution patterns within the genetic lineages and the estimated pan genome of cpe-carrying isolates had a larger core genome and a smaller accessory genome compared to that of 290 strains. We characterize cpe-negative strains related to chromosomal cpe-carrying strains elucidating the origin of these strains and disclose two distinct groups of chromosomal cpe-carrying strains with different virulence characteristics, spore heat resistance properties, and, presumably, ecological niche. Finally, an antibiotic-associated diarrhea isolate carrying two copies of the enterotoxin cpe gene and the associated genetic lineage with the potential for the emergence of similar strains are outlined. With C. perfringens as an example, implications of input genome quality for pan genome analysis are discussed. Our study furthers the understanding of genome epidemiology and population structure of enteropathogenic C. perfringens and brings new insight into this important pathogen and its reservoirs.

Re-Emergence and Spread of Haemorrhagic Septicaemia in Germany: The Wolf as a Vector?

Since 2010, outbreaks of haemorrhagic septicaemia (HS) caused by Pasteurella (P.) multocida capsular type B (PmB) emerged in Germany. In 2017, we noticed a close spatiotemporal relationship between HS outbreak sites and wolf (Canis lupus) territories. Thus, the main objectives of our study were to investigate the molecular epidemiology of German PmB-HS-isolates and to assess the role of wolves as putative vectors of this pathogen. We collected 83 PmB isolates from HS outbreaks that occurred between 2010 and 2019 and sampled 150 wolves, which were found dead in the years 2017 to 2019, revealing another three PmB isolates. A maximum-likelihood-based phylogeny of the core genomes of 65 PmB-HS-isolates and the three PmB-wolf-isolates showed high relatedness. Furthermore, all belonged to capsular:LPS:MLST genotype B:L2:ST122^RIRDC and showed highly similar virulence gene profiles, but clustered separately from 35 global ST122^RIRDC strains. Our data revealed that German HS outbreaks were caused by a distinct genomic lineage of PmB-ST122 strains, hinting towards an independent, ongoing epidemiologic event. We demonstrated for the first time, that carnivores, i.e., wolves, might harbour PmB as a part of their oropharyngeal microbiota. Furthermore, the results of our study imply that wolves can carry the pathogen over long distances, indicating a major role of that animal species in the ongoing epidemiological event of HS in Germany.

Pan-Genome of Novel Pantoea stewartii subsp. indologenes Reveals Genes Involved in Onion Pathogenicity and Evidence of Lateral Gene Transfer

Pantoea stewartii subsp. indologenes (Psi) is a causative agent of leafspot on foxtail millet and pearl millet; however, novel strains were recently identified that are pathogenic on onions. Our recent host range evaluation study identified two pathovars; P. stewartii subsp. indologenes pv. cepacicola pv. nov. and P. stewartii subsp. indologenes pv. setariae pv. nov. that are pathogenic on onions and millets or on millets only, respectively. In the current study, we developed a pan-genome using the whole genome sequencing of newly identified/classified Psi strains from both pathovars [pv. cepacicola (n = 4) and pv. setariae (n = 13)]. The full spectrum of the pan-genome contained 7030 genes. Among these, 3546 (present in genomes of all 17 strains) were the core genes that were a subset of 3682 soft-core genes (present in ≥16 strains). The accessory genome included 1308 shell genes and 2040 cloud genes (present in ≤2 strains). The pan-genome showed a clear linear progression with >6000 genes, suggesting that the pan-genome of Psi is open. Comparative phylogenetic analysis showed differences in phylogenetic clustering of Pantoea spp. using PAVs/wgMLST approach in comparison with core genome SNPs-based phylogeny. Further, we conducted a horizontal gene transfer (HGT) study using Psi strains from both pathovars along with strains from other Pantoea species, namely, P. stewartii subsp. stewartii LMG 2715T, P. ananatis LMG 2665T, P. agglomerans LMG L15, and P. allii LMG 24248T. A total of 317 HGT events among four Pantoea species were identified with most gene transfer events occurring between Psi pv. cepacicola and Psi pv. setariae. Pan-GWAS analysis predicted a total of 154 genes, including seven gene-clusters, which were associated with the pathogenicity phenotype (necrosis on seedling) on onions. One of the gene-clusters contained 11 genes with known functions and was found to be chromosomally located.

Single-Nucleotide Polymorphism-Based Genetic Diversity Analysis of Clinical Pseudomonas aeruginosa Isolates

Extensive use of next-generation sequencing has the potential to transform our knowledge on how genomic variation within bacterial species impacts phenotypic versatility. Because different environments have unique selection pressures, they drive divergent evolution. However, there is also parallel or convergent evolution of traits in independent bacterial isolates inhabiting similar environments. The application of tools to describe population-wide genomic diversity provides an opportunity to measure the predictability of genetic changes underlying adaptation. Here, we describe patterns of sequence variations in the core genome among 99 individual Pseudomonas aeruginosa clinical isolates and identified single-nucleotide polymorphisms that are the basis for branching of the phylogenetic tree. We also identified single-nucleotide polymorphisms that were acquired independently, in separate lineages, and not through inheritance from a common ancestor. Although our results demonstrate that the Pseudomonas aeruginosa core genome is highly conserved and in general, not subject to adaptive evolution, instances of parallel evolution will provide an opportunity to uncover genetic changes that underlie phenotypic diversity.

CoreCruncher: Fast and Robust Construction of Core Genomes in Large Prokaryotic Data Sets

The core genome represents the set of genes shared by all, or nearly all, strains of a given population or species of prokaryotes. Inferring the core genome is integral to many genomic analyses, however, most methods rely on the comparison of all the pairs of genomes; a step that is becoming increasingly difficult given the massive accumulation of genomic data. Here, we present CoreCruncher; a program that robustly and rapidly constructs core genomes across hundreds or thousands of genomes. CoreCruncher does not compute all pairwise genome comparisons and uses a heuristic based on the distributions of identity scores to classify sequences as orthologs or paralogs/xenologs. Although it is much faster than current methods, our results indicate that our approach is more conservative than other tools and less sensitive to the presence of paralogs and xenologs. CoreCruncher is freely available from: https://github.com/lbobay/CoreCruncher. CoreCruncher is written in Python 3.7 and can also run on Python 2.7 without modification. It requires the python library Numpy and either Usearch or Blast. Certain options require the programs muscle or mafft.

Phylogenomic Characterization of a Novel Corynebacterium Species Associated with Fatal Diphtheritic Stomatitis in Endangered Yellow-Eyed Penguins

Yellow-eyed penguins, Megadyptes antipodes, are an endangered species that are endemic to New Zealand. Outbreaks of diphtheritic stomatitis have caused significant mortality for this species, especially among young chicks. In this study, we isolated 16 Corynebacterium sp. isolates from the oral cavities of 2- to 14-day-old chicks at a range of infection stages and sequenced the genomes to understand their virulence mechanisms. Phylogenomic and matrix-assisted laser desorption ionization–time of flight (MALDI-TOF) characterization indicate that these strains belong to a novel Corynebacterium species. A simple multiplex PCR-based diagnostic assay has been developed to identify these strains rapidly and reliably. Similar to other corynebacteria, genomic islands and prophages introduced significant diversity among these strains that has potentially led to minor functional variations between the two lineages. Despite the presence of multiple corynebacterial virulence genes and a spaDEF-type pilus gene cluster among these strains, the survival rate was much higher in Galleria mellonella larvae than in those inoculated with Corynebacterium ulcerans NZRM 818 and Corynebacterium pseudotuberculosis NZRM 3004. Therefore, these strains are opportunistic pathogens causing high mortality among young penguin chicks due to a less-developed immune system.

IMPORTANCE Yellow-eyed penguins, Megadyptes antipodes, are endangered species with a sharp decline in the numbers of breeding pairs over the last 2 decades. Diphtheritic stomatitis, characterized by a thick fibrinopurulent exudate in the oral cavities and symptoms, including inanition and significant weight loss, is responsible for significant mortality among the young chicks. These chicks are treated with antibiotics, amoxicillin-clavulanic acid or enrofloxacin, but do not always recover from the infection. The pathogen causing these infections and the mechanism of pathogenesis are unclear. This study has identified a novel Corynebacterium species to be associated with diphtheritic stomatitis in yellow-eyed penguins with potential virulence genes that are likely involved in pathogenesis. Importantly, a gene encoding an exotoxin, phospholipase D, is present among these strains. The inactivated form of this enzyme could potentially be used as an effective vaccine to protect these penguins from infection.

Transition From PCR-Ribotyping to Whole Genome Sequencing Based Typing of Clostridioides difficile

Clostridioides difficile causes nosocomial outbreaks which can lead to severe and even life-threatening colitis. Rapid molecular diagnostic tests allow the identification of toxin-producing, potentially hypervirulent strains, which is critical for patient management and infection control. PCR-ribotyping has been used for decades as the reference standard to investigate transmission in suspected outbreaks. However, the introduction of whole genome sequencing (WGS) for molecular epidemiology provides a realistic alternative to PCR-ribotyping. In this transition phase it is crucial to understand the strengths and weaknesses of the two technologies, and to assess their correlation. We aimed to investigate ribotype prediction from WGS data, and options for analysis at different levels of analytical granularity. Ribotypes cannot be directly determined from short read Illumina sequence data as the rRNA operons including the ribotype-defining ISR fragments collapse in genome assemblies, and comparison with traditional PCR-ribotyping results becomes impossible. Ribotype extraction from long read Oxford nanopore data also requires optimization. We have compared WGS-based typing with PCR-ribotyping in nearly 300 clinical and environmental isolates from Switzerland, and in addition from the Enterobase database (n=1778). Our results show that while multi-locus sequence type (MLST) often correlates with a specific ribotype, the agreement is not complete, and for some ribotypes the resolution is insufficient. Using core genome MLST (cgMLST) analysis, there is an improved resolution and ribotypes can often be predicted within clusters, using cutoffs of 30-50 allele differences. The exceptions are ribotypes within known ribotype complexes such as RT078/RT106, where the genome differences in cgMLST do not reflect the ribotype segregation. We show that different ribotype clusters display different degrees of diversity, which could be important for the definition of ribotype cluster specific cutoffs. WGS-based analysis offers the ultimate resolution to the SNP level, enabling exploration of patient-to-patient transmission. PCR-ribotyping does not sufficiently discriminate to prove nosocomial transmission with certainty. We discuss the associated challenges and opportunities in a switch to WGS from conventional ribotyping for C. difficile.

Pan and Core Genome Analysis of 183 Mycobacterium tuberculosis Strains Revealed a High Inter-Species Diversity among the Human Adapted Strains

Tuberculosis (TB) is an airborne communicable disease with high morbidity and mortality rates, especially in developing countries. The causal agents of TB belong to the complex Mycobacterium tuberculosis (MTBc), which is composed of different human and animal TB associated species. Some animal associated species have zoonotic potential and add to the burden of TB management. The BCG ("Bacillus Calmette-Guérin") vaccine is widely used for the prevention against TB, but its use is limited in immunocompromised patients and animals due to the adverse effects and disseminated life-threatening complications. In this study, we aimed to carry out a comparative genome analysis between the human adapted species including BCG vaccine strains to identify and pinpoint the conserved genes related to the virulence across all the species, which could add a new value for vaccine development. For this purpose, the sequences of 183 Mycobacterium tuberculosis (MTB) strains were retrieved from the freely available WGS dataset at NCBI. The species included: 168 sensu stricto MTB species with other human MTB complex associated strains: M. tuberculosis var. africanum (3), M. tuberculosis var. bovis (2 draft genomes) and 10 BCG species, which enabled the analysis of core genome which contains the conserved genes and some virulence factor determinants. Further, a phylogenetic tree was constructed including the genomes of human (183); animals MTB adapted strains (6) and the environmental Mycobacterium strain "M. canettii". Our results showed that the core genome consists of 1166 conserved genes among these species, which represents a small portion of the pangenome (7036 genes). The remaining genes in the pangenome (5870) are accessory genes, adding a high inter-species diversity. Further, the core genome includes several virulence-associated genes and this could explain the rare infectiousness potential of some attenuated vaccine strains in some patients. This study reveals that low number of conserved genes in human adapted MTBc species and high inter-species diversity of the pan-genome could be considered for vaccine candidate development.

A New Catalog of Structural Variants in 1,301 A. thaliana Lines from Africa, Eurasia, and North America Reveals a Signature of Balancing Selection at Defense Response Genes

Genomic variation in the model plant Arabidopsis thaliana has been extensively used to understand evolutionary processes in natural populations, mainly focusing on single-nucleotide polymorphisms. Conversely, structural variation has been largely ignored in spite of its potential to dramatically affect phenotype. Here, we identify 155,440 indels and structural variants ranging in size from 1 bp to 10 kb, including presence/absence variants (PAVs), inversions, and tandem duplications in 1,301 A. thaliana natural accessions from Morocco, Madeira, Europe, Asia, and North America. We show evidence for strong purifying selection on PAVs in genes, in particular for housekeeping genes and homeobox genes, and we find that PAVs are concentrated in defense-related genes (R-genes, secondary metabolites) and F-box genes. This implies the presence of a "core" genome underlying basic cellular processes and a "flexible" genome that includes genes that may be important in spatially or temporally varying selection. Further, we find an excess of intermediate frequency PAVs in defense response genes in nearly all populations studied, consistent with a history of balancing selection on this class of genes. Finally, we find that PAVs in genes involved in the cold requirement for flowering (vernalization) and drought response are strongly associated with temperature at the sites of origin.

Comparative Genomics of Emerging Lineages and Mobile Resistomes of Contemporary Broiler Strains of Salmonella Infantis and E. coli

Introduction

Commensal and pathogenic strains of multidrug-resistant (MDR) Escherichia coli and non-typhoid strains of Salmonella represent a growing foodborne threat from foods of poultry origin. MDR strains of Salmonella Infantis and E. coli are frequently isolated from broiler chicks and the simultaneous presence of these two enteric bacterial species would potentially allow the exchange of mobile resistance determinants.

Objectives

In order to understand possible genomic relations and to obtain a first insight into the potential interplay of resistance genes between enteric bacteria, we compared genomic diversity and mobile resistomes of S. Infantis and E. coli from broiler sources.

Results

The core genome MLST analysis of 56 S. Infantis and 90 E. coli contemporary strains revealed a high genomic heterogeneity of broiler E. coli. It also allowed the first insight into the genomic diversity of the MDR clone B2 of S. Infantis, which is endemic in Hungary. We also identified new MDR lineages for S. Infantis (ST7081 and ST7082) and for E. coli (ST8702 and ST10088). Comparative analysis of antibiotic resistance genes and plasmid types revealed a relatively narrow interface between the mobile resistomes of E. coli and S. Infantis. The mobile resistance genes tet(A), aadA1, and sul1 were identified at an overall high prevalence in both species. This gene association is characteristic to the plasmid pSI54/04 of the epidemic clone B2 of S. Infantis. Simultaneous presence of these genes and of IncI plasmids of the same subtype in cohabitant caecal strains of E. coli and S. Infantis suggests an important role of these plasmid families in a possible interplay of resistance genes between S. Infantis and E. coli in broilers.

Conclusion

This is the first comparative genomic analysis of contemporary broiler strains of S. Infantis and E. coli. The diversity of mobile resistomes suggests that commensal E. coli could be potential reservoirs of resistance for S. Infantis, but so far only a few plasmid types and mobile resistance genes could be considered as potentially exchangeable between these two species. Among these, IncI1 plasmids could make the greatest contribution to the microevolution and genetic interaction between E. coli and S. Infantis.

Pangenomics of the Symbiotic Rhizobiales. Core and Accessory Functions Across a Group Endowed with High Levels of Genomic Plasticity

Pangenome analyses reveal major clues on evolutionary instances and critical genome core conservation. The order Rhizobiales encompasses several families with rather disparate ecological attitudes. Among them, Rhizobiaceae, Bradyrhizobiaceae, Phyllobacteriacreae and Xanthobacteriaceae, include members proficient in mutualistic symbioses with plants based on the bacterial conversion of N₂ into ammonia (nitrogen-fixation). The pangenome of 12 nitrogen-fixing plant symbionts of the Rhizobiales was analyzed yielding total 37,364 loci, with a core genome constituting 700 genes. The percentage of core genes averaged 10.2% over single genomes, and between 5% to 7% were found to be plasmid-associated. The comparison between a representative reference genome and the core genome subset, showed the core genome highly enriched in genes for macromolecule metabolism, ribosomal constituents and overall translation machinery, while membrane/periplasm-associated genes, and transport domains resulted under-represented. The analysis of protein functions revealed that between 1.7% and 4.9% of core proteins could putatively have different functions.

Taxonomic Evaluation of the Heyndrickxia (Basonym Bacillus) sporothermodurans Group (H. sporothermodurans, H. vini, H. oleronia) Based on Whole Genome Sequences

The genetic heterogeneity of Heyndrickxia sporothermodurans (formerly Bacillussporothermodurans) was evaluated using whole genome sequencing. The genomes of 29 previously identified Heyndrickxiasporothermodurans and two Heyndrickxia vini strains isolated from ultra-high-temperature (UHT)-treated milk were sequenced by short-read (Illumina) sequencing. After sequence analysis, the two H. vini strains could be reclassified as H. sporothermodurans. In addition, the genomes of the H.sporothermodurans type strain (DSM 10599^T) and the closest phylogenetic neighbors Heyndrickxiaoleronia (DSM 9356^T) and Heyndrickxia vini (JCM 19841^T) were also sequenced using both long (MinION) and short-read (Illumina) sequencing. By hybrid sequence assembly, the genome of the H. sporothermodurans type strain was enlarged by 15% relative to the short-read assembly. This noticeable increase was probably due to numerous mobile elements in the genome that are presumptively related to spore heat tolerance. Phylogenetic studies based on 16S rDNA gene sequence, core genome, single-nucleotide polymorphisms and ANI/dDDH, showed that H. vini is highly related to H. sporothermodurans. When examining the genome sequences of all H.sporothermodurans strains from this study, together with 4 H. sporothermodurans genomes available in the GenBank database, the majority of the 36 strains examined occurred in a clonal lineage with less than 100 SNPs. These data substantiate previous reports on the existence and spread of a genetically highly homogenous and heat resistant spore clone, i.e., the HRS-clone.

Core and Accessory Genome Analysis of Vibrio mimicus

Vibrio mimicus is an emerging pathogen, mainly associated with contaminated seafood consumption. However, little is known about its evolution, biodiversity, and pathogenic potential. This study analyzes the pan-, core, and accessory genomes of nine V. mimicus strains. The core genome yielded 2424 genes in chromosome I (ChI) and 822 genes in chromosome II (ChII), with an accessory genome comprising an average of 10.9% of the whole genome for ChI and 29% for ChII. Core genome phylogenetic trees were obtained, and V. mimicus ATCC-33654 strain was the closest to the outgroup in both chromosomes. Additionally, a phylogenetic study of eight conserved genes (ftsZ, gapA, gyrB, topA, rpoA, recA, mreB, and pyrH), including Vibrio cholerae, Vibrio parilis, Vibrio metoecus, and Vibrio caribbenthicus, clearly showed clade differentiation. The main virulence genes found in ChI corresponded with type I secretion proteins, extracellular components, flagellar proteins, and potential regulators, while, in ChII, the main categories were type-I secretion proteins, chemotaxis proteins, and antibiotic resistance proteins. The accessory genome was characterized by the presence of mobile elements and toxin encoding genes in both chromosomes. Based on the genome atlas, it was possible to characterize differential regions between strains. The pan-genome of V. mimicus encompassed 3539 genes for ChI and 2355 genes for ChII. These results give us an insight into the virulence and gene content of V. mimicus, as well as constitute the first approach to its diversity.

OXA-181-Like Carbapenemases in Klebsiella pneumoniae ST14, ST15, ST23, ST48, and ST231 from Septicemic Neonates: Coexistence with NDM-5, Resistome, Transmissibility, and Genome Diversity

Neonatal sepsis is a leading cause of neonatal mortality in low- and middle-income countries (LMICs). Treatment of sepsis in this vulnerable population is dependent on antimicrobials, and resistance to these life-saving antimicrobials is worrisome.

Studies on the epidemiology and genomes of isolates harboring OXA-48-like genes in septicemic neonates are rare. Here, isolates producing these carbapenemases which emerged and persisted in an Indian neonatal unit were characterized in terms of their resistome, transmissibility, and genome diversity. Antibiotic susceptibility and whole-genome sequencing were carried out. The sequence types, resistome, virulome, mobile genetic elements, and transmissibility of carbapenem-resistant plasmids were evaluated. Core genome analysis of isolates was shown in a global context with other OXA-48-like carbapenemase-harboring genomes, including those from neonatal studies. Eleven OXA-48-like carbapenemase-producing Klebsiella pneumoniae (bla_OXA-181, n = 7 and bla_OXA-232, n = 4) isolates belonging to diverse sequence types (ST14, ST15, ST23, ST48, and ST231) were identified. bla_{OXA-181/OXA-232} and bla_NDM-5 were found in a high-risk clone, ST14 (n = 4). bla_{OXA-181/OXA-232} were in small, nonconjugative ColKP3 plasmids located on truncated Tn2013, whereas bla_NDM-5 was in self-transmissible, conjugative IncFII plasmids, within truncated Tn125. Conjugal transfer of bla_{OXA-181/OXA-232} was observed in the presence of bla_NDM-5. The study strains were diverse among themselves and showed various levels of relatedness with non-neonatal strains from different parts of the world and similarity with neonatal strains from Tanzania and Ghana when compared with a representative collection of carbapenemase-positive K. pneumoniae strains. We found that bla_{OXA-181/OXA-232}-harboring isolates from a single neonatal unit had remarkably diverse genomes, ruling out clonal spread and emphasizing the extent of plasmid spreading across different STs. This study is probably the first to report the coexistence of bla_OXA-181/232 and bla_NDM-5 in neonatal isolates.

IMPORTANCE Neonatal sepsis is a leading cause of neonatal mortality in low- and middle-income countries (LMICs). Treatment of sepsis in this vulnerable population is dependent on antimicrobials, and resistance to these life-saving antimicrobials is worrisome. Carbapenemases, enzymes produced by bacteria, can make these antimicrobials useless. Our study describes how OXA-48-like carbapenemases in neonatal septicemic Klebsiella pneumoniae shows remarkable diversity in the genomes of the strains and relatedness with strains from other parts of world and also to some neonatal outbreak strains. It is also the first to describe such resistance due to coproduction of dual carbapenemases, (OXA)-48 and New Delhi metallo-β-lactamase-5, in Klebsiella pneumoniae from neonatal settings. Carbapenemase genes situated on plasmids within high-risk international clones, as seen here, increase the ease and transfer of resistant genetic material. With the WHO treatment protocols not adequately poised to handle such infections, prompt attention to neonatal health care is required.

A Vibrio cholerae Core Genome Multilocus Sequence Typing Scheme To Facilitate the Epidemiological Study of Cholera

Core genome multilocus sequence typing (cgMLST) has gained popularity in recent years in epidemiological research and subspecies-level classification. cgMLST retains the intuitive nature of traditional MLST but offers much greater resolution by utilizing significantly larger portions of the genome. Here, we introduce a cgMLST scheme for Vibrio cholerae, a bacterium abundant in marine and freshwater environments and the etiologic agent of cholera. A set of 2,443 core genes ubiquitous in V. cholerae were used to analyze a comprehensive data set of 1,262 clinical and environmental strains collected from 52 countries, including 65 newly sequenced genomes in this study. We established a sublineage threshold based on 133 allelic differences that creates clusters nearly identical to traditional MLST types, providing backwards compatibility to new cgMLST classifications. We also defined an outbreak threshold based on seven allelic differences that is capable of identifying strains from the same outbreak and closely related isolates that could give clues on outbreak origin. Using cgMLST, we confirmed the South Asian origin of modern epidemics and identified clustering affinity among sublineages of environmental isolates from the same geographic origin. Advantages of this method are highlighted by direct comparison with existing classification methods, such as MLST and single-nucleotide polymorphism-based methods. cgMLST outperforms all existing methods in terms of resolution, standardization, and ease of use. We anticipate this scheme will serve as a basis for a universally applicable and standardized classification system for V. cholerae research and epidemiological surveillance in the future. This cgMLST scheme is publicly available on PubMLST (https://pubmlst.org/vcholerae/).IMPORTANCE Toxigenic Vibrio cholerae isolates of the O1 and O139 serogroups are the causative agents of cholera, an acute diarrheal disease that plagued the world for centuries, if not millennia. Here, we introduce a core genome multilocus sequence typing scheme for V. cholerae Using this scheme, we have standardized the definition for subspecies-level classification, facilitating global collaboration in the surveillance of V. cholerae In addition, this typing scheme allows for quick identification of outbreak-related isolates that can guide subsequent analyses, serving as an important first step in epidemiological research. This scheme is also easily scalable to analyze thousands of isolates at various levels of resolution, making it an invaluable tool for large-scale ecological and evolutionary analyses.

Comparative Genome Analysis of 33 Chlamydia Strains Reveals Characteristic Features of Chlamydia Psittaci and Closely Related Species

To identify genome-based features characteristic of the avian and human pathogen Chlamydia(C.) psittaci and related chlamydiae, we analyzed whole-genome sequences of 33 strains belonging to 12 species. Using a novel genome analysis tool termed Roary ILP Bacterial Annotation Pipeline (RIBAP), this panel of strains was shown to share a large core genome comprising 784 genes and representing approximately 80% of individual genomes. Analyzing the most variable genomic sites, we identified a set of features of C. psittaci that in its entirety is characteristic of this species: (i) a relatively short plasticity zone of less than 30,000 nt without a tryptophan operon (also in C. abortus, C. avium, C. gallinacea, C. pneumoniae), (ii) a characteristic set of of Inc proteins comprising IncA, B, C, V, X, Y (with homologs in C. abortus, C. caviae and C. felis as closest relatives), (iii) a 502-aa SinC protein, the largest among Chlamydia spp., and (iv) an elevated number of Pmp proteins of subtype G (14 in C. psittaci, 14 in Cand. C. ibidis). In combination with future functional studies, the common and distinctive criteria revealed in this study provide important clues for understanding the complexity of host-specific behavior of individual Chlamydia spp.

Acinetobacter portensis sp. nov. and Acinetobacter guerrae sp. nov., isolated from raw meat

The taxonomic status of six strains of Acinetobacter obtained from meat samples, collected from supermarkets in Porto, Portugal, was investigated using polyphasic analysis. Partial rpoB sequence similarities lower than 95 % to other Acinetobacter species with validly published names led to the hypothesis that these strains represented novel species. This was confirmed based on comparative multilocus sequence analysis, which included the gyrB, recA and 16S rRNA genes, revealing that these strains represented two coherent lineages that were distinct from each other and from all known species. The names Acinetobacter portensis sp. nov. (comprising four strains) and Acinetobacter guerrae sp. nov. (comprising two strains) are proposed for these novel species. The species status of these two groups was confirmed by low (below 95 %) whole-genome sequence average nucleotide identity values and low (below 70 %) digital DNA-DNA hybridization similarities between the whole-genome sequences of the proposed type strains of each novel species and the representatives of the known Acinetobacter species. Phylogenomic treeing from core genome analysis supported these results. The coherence of each new species lineage was supported by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry differentiation of the species at the protein level, by cellular fatty acid profiles, and by unique and differential combinations of metabolic and physiological properties shared by each novel species. The type strain of A. portensis sp. nov. is AC 877^T (=CCUG 68672^T=CCM 8789^T) and the type strain of A. guerrae sp. nov. is AC 1271^T (=CCUG 68674^T=CCM 8791^T).

Comparison of Molecular Subtyping and Antimicrobial Resistance Detection Methods Used in a Large Multistate Outbreak of Extensively Drug-Resistant Campylobacter jejuni Infections Linked to Pet Store Puppies

Campylobacter jejuni is a leading cause of enteric bacterial illness in the United States. Traditional molecular subtyping methods, such as pulsed-field gel electrophoresis (PFGE) and 7-gene multilocus sequence typing (MLST), provided limited resolution to adequately identify C. jejuni outbreaks and separate out sporadic isolates during outbreak investigations. Whole-genome sequencing (WGS) has emerged as a powerful tool for C. jejuni outbreak detection. In this investigation, 45 human and 11 puppy isolates obtained during a 2016-2018 outbreak linked to pet store puppies were sequenced. Core genome multilocus sequence typing (cgMLST) and high-quality single nucleotide polymorphism (hqSNP) analysis of the sequence data separated the isolates into the same two clades containing minor within-clade differences; however, cgMLST analysis does not require selection of an appropriate reference genome, making the method preferable to hqSNP analysis for Campylobacter surveillance and cluster detection. The isolates were classified as sequence type 2109 (ST2109)-a rarely seen MLST sequence type. PFGE was performed on 38 human and 10 puppy isolates; PFGE patterns did not reliably predict clustering by cgMLST analysis. Genetic detection of antimicrobial resistance determinants predicted that all outbreak-associated isolates would be resistant to six drug classes. Traditional antimicrobial susceptibility testing (AST) confirmed a high correlation between genotypic and phenotypic antimicrobial resistance determinations. WGS analysis linked C. jejuni isolates in humans and pet store puppies even when canine exposure information was unknown, aiding the epidemiological investigation during the outbreak. WGS data were also used to quickly identify the highly drug-resistant profile of these outbreak-associated C. jejuni isolates.

Rapid Increase of CTX-M-Producing Shigella sonnei Isolates in Switzerland Due to Spread of Common Plasmids and International Clones

The Swiss Centre for Antibiotic Resistance (ANRESIS) has recently noted an increase of extended-spectrum cephalosporin-resistant (ESC-R) Shigella sonnei isolates nationwide (3.8% in 2016 versus 37.5% in 2019). To understand this phenomenon, we analyzed 25 representative isolates (of which 14 were ESC-R) collected in Switzerland during 2016 to 2019. Whole-genome sequencing was achieved using both the Illumina and the Nanopore platforms. Both ESC-R and extended-spectrum cephalosporin-susceptible isolates belonged to sequence type 152 (ST152). The ESC-R isolates carried bla _CTX-M-3 in IncI1-pST57 (n = 5), bla _CTX-M-15 in IncFII (F2:A-:B-) (n = 5), bla _CTX-M-15 in IncI1-pST16, and bla _CTX-M-27, bla _CTX-M-55, or bla _CTX-M-134 in other IncFII plasmids (n = 1 each). Plasmids having the same bla and Inc group exhibited high degrees of genetic identity to each other but also to plasmids previously reported in other Enterobacterales Core-genome analysis showed that there were 4 main clusters, each of which included strains that differed by <58 single nucleotide variants (SNVs) and that consisted of both bla _CTX-M-positive and bla _CTX-M-negative isolates. Moreover, most isolates belonging to the same cluster shared an identical core-genome sequence type (cgST). For instance, cluster 1 included 4 isolates of cgST113036, of which only 3 harbored the IncI1-pST57 bla _CTX-M-3-positive plasmid. The 25 S. sonnei isolates were also subjected to phylogenetic comparison with deposited international strains. As a result, matching isolates (isolates that had the same cgST and that differed by <8 SNVs) have been reported in the United Kingdom, the United States, France, and the Netherlands. Overall, our results suggest that some common S. sonnei clusters can spread between continents and can be imported into other nations after international trips. Such clusters include, in part, isolates that do not possess bla _ESBL-harboring plasmids, indicating their tendency to acquire them from other Enterobacterales.

Reclassification of the Taxonomic Framework of Orders Cellvibrionales, Oceanospirillales, Pseudomonadales, and Alteromonadales in Class Gammaproteobacteria through Phylogenomic Tree Analysis

Orders Oceanospirillales and Pseudomonadales play important roles in various ecosystems as the keystone taxa of microbiomes. However, the two orders present a close evolutionary relationship, which might have caused taxonomic misinterpretation and resulted in an incorrect understanding of their evolutionary history. In this study, first, we used the 16S rRNA gene sequences of 2,049 species of Gammaproteobacteria to build a phylogenetic tree, which demonstrated that reports regarding the evolutionary relationship of orders Cellvibrionales, Oceanospirillales, and Pseudomonadales based on a single conserved gene with a poor resolution have been conflicting; in particular, the major families Moraxellaceae and Pseudomonadaceae of order Pseudomonadales were separated from orders Cellvibrionales and Oceanospirillales Subsequently, we constructed the bac120 trees of all representative reference genomes of class Gammaproteobacteria based on 120 ubiquitous single-copy proteins from bacteria and a phylogenomic tree based on the 119 core genes of 257 reference genomes obtained from orders Cellvibrionales, Oceanospirillales, and Pseudomonadales to cross validate and infer their intrinsic evolutionary relationships. These results indicated that two novel orders, Moraxellales ord. nov. and Kangiellales ord. nov., and three novel families, Marinobacteraceae fam. nov., Perlucidibacaceae fam. nov., and Zooshikellaceae fam. nov., should be proposed. Additionally, orders Cellvibrionales and Oceanospirillales were merged into the order Pseudomonadales except for families Moraxellaceae and Kangiellaceae in class Gammaproteobacteria, which currently includes 18 families. Our work sheds some light on the evolutionary history of class Gammaproteobacteria, which could facilitate the detection and taxonomic analysis of natural communities.IMPORTANCE The orders Cellvibrionales, Oceanospirillales, and Pseudomonadales, as three major orders of the largest bacterial class, Gammaproteobacteria, play important roles in various ecosystems as the keystone taxa of microbiomes, but their evolutionary relationship is currently polyphyletic and chaotic. Here, we constructed a bac120 tree and core-genome tree and calculated the amino acid identity (AAI) value to explore their intrinsic evolutionary history. In this study, we proposed two novel orders and three novel families. This evolution study vastly reconstructed the taxonomic framework of class Gammaproteobacteria and could provide a more distinct perspective on global distribution and evolutionary patterns of these environmental microorganisms.

In silico Prediction of New Drug Candidates Against the Multidrug-Resistant and Potentially Zoonotic Fish Pathogen Serotype III Streptococcus agalactiae

Streptococcus agalactiae is an invasive multi-host pathogen that causes invasive diseases mainly in newborns, elderly, and individuals with underlying health complications. In fish, S. agalactiae causes streptococcosis, which is characterized by septicemia and neurological signs, and leads to great economic losses to the fish farming industry worldwide. These bacteria can be classified into different serotypes based on capsular antigens, and into different sequence types (ST) based on multilocus sequence typing (MLST). In 2015, serotype III ST283 was identified to be associated with a foodborne invasive disease in non-pregnant immunocompetent humans in Singapore, and the infection was related to raw fish consumption. In addition, a serotype III strain isolated from tilapia in Brazil has been reported to be resistant to five antibiotic classes. This specific serotype can serve as a reservoir of resistance genes and pose a serious threat to public health. Thus, new approaches for the control and treatment of S. agalactiae infections are needed. In the present study, 24 S. agalactiae serotype III complete genomes, isolated from human and fish hosts, were compared. The core genome was identified, and, using bioinformatics tools and subtractive criteria, five proteins were identified as potential drug targets. Furthermore, 5,008 drug-like natural compounds were virtually screened against the identified targets. The ligands with the best binding properties are suggested for further in vitro and in vivo analysis.

First Steps in the Analysis of Prokaryotic Pan-Genomes

Pan-genome is defined as the set of orthologous and unique genes of a specific group of organisms. The pan-genome is composed by the core genome, accessory genome, and species- or strain-specific genes. The pan-genome is considered open or closed based on the alpha value of the Heap law. In an open pan-genome, the number of gene families will continuously increase with the addition of new genomes to the analysis, while in a closed pan-genome, the number of gene families will not increase considerably. The first step of a pan-genome analysis is the homogenization of genome annotation. The same software should be used to annotate genomes, such as GeneMark or RAST. Subsequently, several software are used to calculate the pan-genome such as BPGA, GET_HOMOLOGUES, PGAP, among others. This review presents all these initial steps for those who want to perform a pan-genome analysis, explaining key concepts of the area. Furthermore, we present the pan-genomic analysis of 9 bacterial species. These are the species with the highest number of genomes deposited in GenBank. We also show the influence of the identity and coverage parameters on the prediction of orthologous and paralogous genes. Finally, we cite the perspectives of several research areas where pan-genome analysis can be used to answer important issues.

Comparative genome analysis and mining of secondary metabolites of Paenibacillus polymyxa

Paenibacillus polymyxa is a well-known Gram-positive biocontrol bacterium. It has been reported that many P. polymyxa strains can inhibit bacteria, fungi and other plant pathogens. Paenibacillus polymyxa employs a variety of mechanisms to promote plant growth, so it is necessary to understand the biocontrol ability of bacteria at the genome level. In the present study, thanks to the widespread availability of Paenibacillus genome data and the development of bioinformatics tools, we were able to analyze and mine the genomes of 43 P. polymyxa strains. The strain NCTC4744 was determined not to be P. polymyxa according to digital DNA-DNA hybridization and average nucleotide identity. By analysis of the pan-genome and the core genome, we found that the pan-genome of P. polymyxa was open and that there were 3,192 core genes. In a gene cluster analysis of secondary metabolites, 797 secondary metabolite gene clusters were found, of which 343 are not similar to known clusters and are expected to reveal a large number of new secondary metabolites. We also analyzed the plant growth-promoting genes that were mined and found, surpisingly, that these genes are highly conserved. The results of the present study not only reveal a large number of unknown potential secondary metabolite gene clusters in P. polymyxa, but also suggest that plant growth promotion characteristics are evolutionary adaptations of P. polymyxa to plant-related habitats.

Characterization of Burkholderia cepacia Complex Core Genome and the Underlying Recombination and Positive Selection

Recombination and positive selection are two key factors that play a vital role in pathogenic microorganisms’ population adaptation and diversification. The Burkholderia cepacia complex (Bcc) represents bacterial species with high similarity, which can cause severe infections among cases suffering from the chronic granulomatous disorder and cystic fibrosis (CF). At present, no genome-wide study has been carried out focusing on investigating the core genome of Bcc associated with the two evolutionary forces. The general characteristics of the core genome of Bcc species remain scarce as well. In this study, we explored the core orthologous genes of 116 Bcc strains using comparative genomic analysis and studied the two adaptive evolutionary forces: recombination and positive selection. We estimated 1005 orthogroups consisting entirely of single copy genes. These single copy orthologous genes in some Cluster of Orthologous Groups (COG) categories showed significant differences in the comparison of several evolutionary properties, and the encoding proteins were relatively simple and compact. Our findings showed that 5.8% of the core orthologous genes strongly supported recombination; in the meantime, 1.1% supported positive selection. We found that genes involved in protein synthesis as well as material transport and metabolism are favored by selection pressure. More importantly, homologous recombination contributed more genetic variation to a large number of genes and largely maintained the genetic cohesion in Bcc. This high level of recombination between Bcc species blurs their taxonomic boundaries, which leads Bcc species to be difficult or impossible to distinguish phenotypically and genotypically.

Pan-Genome Analyses of Geobacillus spp. Reveal Genetic Characteristics and Composting Potential

The genus Geobacillus is abundant in ecological diversity and is also well-known as an authoritative source for producing various thermostable enzymes. Although it is clear now that Geobacillus evolved from Bacillus, relatively little knowledge has been obtained regarding its evolutionary mechanism, which might also contribute to its ecological diversity and biotechnology potential. Here, a statistical comparison of thirty-two Geobacillus genomes was performed with a specific focus on pan- and core genomes. The pan-genome of this set of Geobacillus strains contained 14,913 genes, and the core genome contained 940 genes. The Clusters of Orthologous Groups (COG) and Carbohydrate-Active Enzymes (CAZymes) analysis revealed that the Geobacillus strains had huge potential industrial application in composting for agricultural waste management. Detailed comparative analyses showed that basic functional classes and housekeeping genes were conserved in the core genome, while genes associated with environmental interaction or energy metabolism were more enriched in the pan-genome. Therefore, the evolution of Geobacillus seems to be guided by environmental parameters. In addition, horizontal gene transfer (HGT) events among different Geobacillus species were detected. Altogether, pan-genome analysis was a useful method for detecting the evolutionary mechanism, and Geobacillus’ evolution was directed by the environment and HGT events.

Whole Genome Sequencing and Characteristics of mcr-1-Harboring Plasmids of Porcine Escherichia coli Isolates Belonging to the High-Risk Clone O25b:H4-ST131 Clade B

Porcine Escherichia coli ST131 isolates are scarcely documented. Here, whole genome sequencing and core genome (CG) and plasmidome analysis of seven isolates collected from diarrheic piglets and four from pork meat were performed. All of the 11 ST131 isolates belonged to serotype O25b:H4 and clade B and showed fimH22 allele or mutational derivatives. The 11 porcine isolates possessed virulence traits that classified the isolates as avian pathogenic, uropathogenic, and extraintestinal pathogenic E. coli–like (APEC-, UPEC-, and ExPEC-like) and constituted virotype D. The CG was performed for all porcine isolates in addition to 73 ST131 reference isolates from different origins. Within clade B, the CG showed nine subclusters, allowing us to describe five new subclades (B6, B6-like, B7, B8, and B9). There was an association between subclade B6, PST43, virotype D2, and food origin, whereas subclade B7 included PST9 isolates with virotype D5 from diarrheic piglets (p = 0.007). The distance between human and porcine isolates from subclades B6 and B7 had an average of 20 and 15 SNP/Mb, respectively. [F2:A-:B1]-IncF, ColE1-like, and IncX plasmids were the most prevalent. Besides, IncF plasmids harbored a ColV region frequent among APEC isolates. Antimicrobial resistance genes conferring resistance to penicillin, tetracycline, quinolones, and colistin were the most common. The mcr-1.1 gene was detected in 5 of 11 porcine isolates, integrated into the chromosome of one isolate and into plasmids in the remainder isolates (two MOB_H11/IncHI2-ST4, one MOB_P3/IncX4, and one MOB_F12/IncF [F2:A-:B1] supposedly cointegrated with an IncHI2). The surrounding environments of the mcr-1 cassette showed variability. However, there were conserved structures within the same plasmid family. In conclusion, CG analysis defined five new subclades. The ST131 porcine isolates belonged to new subclades B6 and B7. Moreover, porcine and clinical human isolates were strongly related. The 11 porcine ST131 isolates harbored a wide variety of plasmids, virulence, and resistance genes. Furthermore, epidemic plasmids IncX4 and IncHI2 are responsible for the acquisition of mcr-1.1 gene. We hypothesize that the APEC-IncF plasmid acquired the mcr-1.1 gene via cointegrating an IncHI2 plasmid, which is worrying due to combination of virulence and resistance attributes in a single mobile genetic element.