Bioinformatic genome comparisons for taxonomic and phylogenetic assignments using Aeromonas as a test case

Genome sequencing and comparative genomic analysis of bovine mastitis-associated non-aureus staphylococci and mammaliicocci (NASM) strains from India

We describe the whole-genome sequencing and comparative genomic analysis of 22 mastitis-associated NASM strains isolated from India. The mean genome size of the strains was 2.55 Mbp, with an average GC content of 32.2%. We identified 14 different sequence types (STs) among the 22 NASM strains. Of these, ST1 and ST6 of S. chromogenes were exclusively associated with bovine mastitis. Genome-wide SNP-based minimum spanning tree revealed the intricate phylogenetic relationships among NASM strains from India, categorizing them into five major clades. Interestingly, mastitis-associated strains formed separate subclades in all the NASM species studied, indicating distinct host-specific co-evolution. The study identified 32 antimicrobial resistance (AMR) genes and 53 virulence-associated genes, providing insights into the genetic factors that could contribute to the pathogenicity of NASM species. Some virulence and AMR genes were found in the predicted genomic islands, suggesting possible horizontal transfer events.

Investigating Additive and Replacing Horizontal Gene Transfers Using Phylogenies and Whole Genomes

Horizontal gene transfer (HGT) is fundamental to microbial evolution and adaptation. When a gene is horizontally transferred, it may either add itself as a new gene to the recipient genome (possibly displacing nonhomologous genes) or replace an existing homologous gene. Currently, studies do not usually distinguish between "additive" and "replacing" HGTs, and their relative frequencies, integration mechanisms, and specific roles in microbial evolution are poorly understood. In this work, we develop a novel computational framework for large-scale classification of HGTs as either additive or replacing. Our framework leverages recently developed phylogenetic approaches for HGT detection and classifies HGTs inferred between terminal edges based on gene orderings along genomes and phylogenetic relationships between the microbial species under consideration. The resulting method, called DART, is highly customizable and scalable and can classify a large fraction of inferred HGTs with high confidence and statistical support. Our application of DART to a large dataset of thousands of gene families from 103 Aeromonas genomes provides insights into the relative frequencies, functional biases, and integration mechanisms of additive and replacing HGTs. Among other results, we find that (i) the relative frequency of additive HGT increases with increasing phylogenetic distance, (ii) replacing HGT dominates at shorter phylogenetic distances, (iii) additive and replacing HGTs have strikingly different functional profiles, (iv) homologous recombination in flanking regions of a novel gene may be a frequent integration mechanism for additive HGT, and (v) phages and mobile genetic elements likely play an important role in facilitating additive HGT.

A coordinated attack by a bacterial secretion system and a small molecule drives prey specificity

Vibrio species are recognized for their role in food- and water-borne diseases in humans, fish, and aquatic invertebrates. We screened bacterial strains isolated from raw food shrimp for those that are bactericidal to Vibrio strains. Here we identify and characterize Aeromonas dhakensis strain A603 which shows robust bactericidal activity specifically towards Vibrio and related taxa but less potency toward other Gram-negative species. Using the A603 genome and genetic analysis, we show that two antibacterial mechanisms account for its vibriocidal activity -- a highly potent Type Six Secretion System (T6SS) and biosynthesis of a vibriocidal phenazine-like small molecule, named here as Ad-Phen. Further analysis indicates coregulation between Ad-Phen and a pore-forming T6SS effector TseC, which potentiates V. cholerae to killing by Ad-Phen.

In silico analysis of Ffp1, an ancestral Porphyromonas spp. fimbrillin, shows differences with Fim and Mfa

Background. Scant information is available regarding fimbrillins within the genus Porphyromonas, with the notable exception of those belonging to Porphyromonas gingivalis, which have been extensively researched for several years. Besides fim and mfa, a third P. gingivalis adhesin called filament-forming protein 1 (Ffp1) has recently been described and seems to be pivotal for outer membrane vesicle (OMV) production. Objective. We aimed to investigate the distribution and diversity of type V fimbrillin, particularly Ffp1, in the genus Porphyromonas. Methods. A bioinformatics phylogenomic analysis was conducted using all accessible Porphyromonas genomes to generate a domain search for fimbriae, using hidden Markov model profiles. Results. Ffp1 was identified as the sole fimbrillin present in all analysed genomes. After manual verification (i.e. biocuration) of both structural and functional annotations and 3D modelling, this protein was determined to be a type V fimbrillin, with a closer structural resemblance to a Bacteroides ovatus fimbrillin than to FimA or Mfa1 from P. gingivalis. Conclusion. It appears that Ffp1 is an ancestral fimbria, transmitted through vertical inheritance and present across all Porphyromonas species. Additional investigations are necessary to elucidate the biogenesis of Ffp1 fimbriae and their potential role in OMV production and niche adaptation.

Updates and Current Challenges in Reproductive Microbiome: A Comparative Analysis between Cows and Women

The microbiota plays an important role in numerous physiological processes, pathogenesis, development, and metabolism in different animal species. In humans, several studies have demonstrated an association between the vaginal microbiota and fertility rates, and even success in assisted reproduction techniques. In the context of cattle reproduction, although few studies have addressed the microbiota in a healthy state (which is not associated with diseases that affect the reproductive tract of cows), changes in its composition also seem to influence fertility. This review aims to explain the importance of the reproductive microbiota in female bovines and what is available in the literature regarding its possible role in increasing fertility. What are the challenges involved in this process? Future perspectives on its use and manipulation as a selection or intervention tool. Will it be possible to one day extrapolate the findings to reality and apply them in the field? In short, understanding the role of the reproductive microbiota of female bovines can signal the prospect of increasing production, whether of milk or meat, from the same number of animals, as it can optimize reproductive efficiency and perhaps become an allied tool for the economic profitability and sustainability of livestock farming.

Aeromonas dhakensis: A Zoonotic Bacterium of Increasing Importance in Aquaculture

Aeromonas dhakensis is increasingly recognised to be an important pathogen responsible for disease losses in warm-water aquaculture and, similar to several other Aeromonas species, it can infect humans. Knowledge of A. dhakensis is accumulating, but this species remains relatively under-investigated compared to its close relative, Aeromonas hydrophila. The significance of A. dhakensis may have been overlooked in disease events of aquatic animals due to issues with reliable identification. Critical to appreciating the importance of this pathogen is the application of dependable molecular tools that enable accurate identification and discrimination from A. hydrophila and other motile aeromonads. This review aims to synthesise the key literature on A. dhakensis, particularly with relevance to aquaculture, including knowledge of the bacterium derived from disease case studies in aquatic hosts. Identification methods and strain phylogeny are discussed, with accurate detection important for prompt diagnosis and for distinguishing strains with heightened virulence. Increasing evidence suggests that A. dhakensis may be more virulent than A. hydrophila and correct identification is required to determine the zoonotic risks posed, which includes concerns for antibiotic-resistant strains. This review provides an impetus to improve species identification in the future and screen strain collections of presumptive Aeromonas spp. retrospectively to reveal the true prevalence and impact of A. dhakensis in aquaculture, the environment, and healthcare settings.

Analysis of global Aeromonas caviae genomes revealed that strains carrying T6SS are more common in human gastroenteritis than in environmental sources and are often phylogenetically related

Aeromonas caviae is an emerging human enteric pathogen. However, the genomic features and virulence genes of A. caviae strains from human gastroenteritis and other sources have not been fully elucidated. Here, we conducted a genomic analysis of 565 global A. caviae strains isolated from different sources, including 261 strains isolated from faecal samples of gastroenteritis patients, of which 18 genomes were sequenced in this study. The presence of bacterial virulence genes and secretion systems in A. caviae strains from different sources was compared, and the phylogenetic relationship of A. caviae strains was assessed based on the core genome. The complete genome of A. caviae strain A20-9 isolated from a gastroenteritis patient was obtained in this study, from which 300 putative virulence factors and a T4SS-encoding plasmid, pAC, were identified. Genes encoding T4SS were also identified in a novel genomic island, ACI-1, from other T4SS-positive strains. The prevalence of T4SS was significantly lower in A. caviae strains from gastroenteritis patients than in environmental strains (3 %, P<0.0001 vs 14 %, P<0.01). Conversely, the prevalence of T6SS was significantly higher in A. caviae strains isolated from gastroenteritis patients than in environmental strains (25 %, P<0.05 vs 13  %, P<0.01). Four phylogenetic clusters were formed based on the core genome of 565 A. caviae strains, and strains carrying T6SS often showed close phylogenetic relationships. T3SS, aerolysin and thermostable cytotonic enterotoxin were absent in all 565 A. caviae strains. Our findings provide novel information on the genomic features of A. caviae and suggest that T6SS may play a role in A. caviae-induced human gastroenteritis.

Methylomonas defluvii sp. nov., a type I methane-oxidizing bacterium from a secondary sedimentation tank of a wastewater treatment plant

An aerobic methanotroph was isolated from a secondary sedimentation tank of a wastewater treatment plant and designated strain OY6T. Cells of OY6T were Gram-stain-negative, pink-pigmented, motile rods and contained an intracytoplasmic membrane structure typical of type I methanotrophs. OY6T could grow at a pH range of 4.5-7.5 (optimum pH 6.5) and at temperatures ranging from 20 °C to 37 °C (optimum 30 °C). The major cellular fatty acids were C14 : 0, C16 : 1ω7c/C16 : 1ω6c and C16 : 1ω5c; the predominant respiratory quinone was MQ-8. The genome size was 5.41 Mbp with a DNA G+C content of 51.7 mol%. OY6T represents a member of the family Methylococcaceae of the class Gammaproteobacteria and displayed 95.74-99.64 % 16S rRNA gene sequence similarity to the type strains of species of the genus Methylomonas. Whole-genome comparisons based on average nucleotide identity (ANI) and digital DNA-DNA hybridisation (dDDH) confirmed that OY6T should be classified as representing a novel species. The most closely related type strain was Methylomonas fluvii EbBT, with 16S rRNA gene sequence similarity, ANI by blast (ANIb), ANI by MUMmer (ANIm) and dDDH values of 99.64, 90.46, 91.92 and 44.5 %, respectively. OY6T possessed genes encoding both the particulate methane monooxygenase enzyme and the soluble methane monooxygenase enzyme. It grew only on methane or methanol as carbon sources. On the basis of phenotypic, genetic and phylogenetic data, strain OY6T represents a novel species within the genus Methylomonas for which the name Methylomonas defluvii sp. nov. is proposed, with strain OY6T (=GDMCC 1.4114T=KCTC 8159T=LMG 33371T) as the type strain.

Proteus appendicitidis sp. nov., isolated from the appendiceal pus of an appendicitis patient in Yongzhou, China

A Gram-negative, facultatively anaerobic, short rod-shaped bacterium, designated as strain HZ0627T, was isolated from the appendiceal pus of a patient with appendicitis in Yongzhou, Hunan, China. This strain was subjected to comprehensive phenotypic, phylogenetic, and genomic analyses using polyphasic taxonomic methods. Phylogenetic analysis of the 16S rRNA gene sequence revealed that this strain belonged to the genus Proteus and the family Morganellaceae, whereas that based on the rpoB gene sequence and phylogenomic analysis demonstrated that this strain was distinctly separated from other type strains of Proteus species. Moreover, whole-genome-based analyses, including in silico DNA-DNA hybridization (isDDH) and average nucleotide identity (ANI), revealed that strain HZ0627T had much lower isDDH rates (24.5-55.6%) and ANI (82.04-93.90%) than those of the thresholds (i.e., 70% and 95%, respectively) for species delineation, when compared to the type strains of other Proteus species. The cellular fatty acid profile of strain HZ0627T was dominated by C16:0 (34.5%), cyclo C17:0 (25.8%), C14:0 (12.6%), C16:1 iso I/14:0 3-OH (7.7%), C18:1ω7c/18:1ω6c (6.5%), and C16:1ω7c/16:1ω6c (4.9%), which clearly differentiated it from the documented type strains of Proteus species. In addition, several specific physiological traits, including optimal growth temperature, tolerance to sodium chloride, and carbon source utilization, differed from those of other Proteus species. Therefore, we propose the name Proteus appendicitidis sp. nov. for strain HZ0627T (= CCTCC AB 2022380T = KCTC 92986T), which represents the type strain of this novel Proteus species.

Genomic analysis of clinical Aeromonas isolates reveals genetic diversity but little evidence of genetic determinants for diarrhoeal disease

Aeromonas spp. are associated with a number of infectious syndromes in humans including gastroenteritis and dysentery. Our understanding of the genetic diversity, population structure, virulence determinants and antimicrobial resistance of the genus has been limited by a lack of sequenced genomes linked to metadata. We performed a comprehensive analysis of the whole genome sequences of 447 Aeromonas isolates from children in Karachi, Pakistan, with moderate-to-severe diarrhoea (MSD) and from matched controls without diarrhoea that were collected as part of the Global Enteric Multicenter Study (GEMS). Human-associated Aeromonas isolates exhibited high species diversity and extensive antimicrobial and virulence gene content. Aeromonas caviae, A. dhankensis, A. veronii and A. enteropelogenes were all significantly associated with MSD in at least one cohort group. The maf2 and lafT genes that encode components of polar and lateral flagella, respectively, exhibited a weak association with isolates originating from cases of gastroenteritis.

Taxonomic classification of genus Aeromonas using open reading frame-based binarized structure network analysis

Objectives

Taxonomic assignment based on whole-genome sequencing data facilitates clear demarcation of species within a complex genus. Here, we applied a unique pan-genome phylogenetic method, open reading frame (ORF)-based binarized structure network analysis (OSNA), for taxonomic inference of Aeromonas spp., a complex taxonomic group consisting of 30 species.

Methods

Data from 335 publicly available Aeromonas genomes, including the reference genomes of 30 species, were used to build a phylogenetic tree using OSNA. In OSNA, whole-genome structures are expressed as binary sequences based on the presence or absence of ORFs, and a tree is generated using neighbor-net, a distance-based method for constructing phylogenetic networks from binary sequences. The tree built by OSNA was compared to that constructed by a core-genome single-nucleotide polymorphism (SNP)-based analysis. Furthermore, the orthologous average nucleotide identity (OrthoANI) values of the sequences that clustered in a single clade in the OSNA-based tree were calculated.

Results

The phylogenetic tree constructed with OSNA successfully delineated the majority of species of the genus Aeromonas forming conspecific clades for individual species, which was corroborated by OrthoANI values. Moreover, the OSNA-based phylogenetic tree demonstrated high compositional similarity to the core-genome SNP-based phylogenetic tree, supported by the Fowlkes-Mallows index.

Conclusions

We propose that OSNA is a useful tool in predicting the taxonomic classification of complex bacterial genera.

Taxonomic resolution of different 16S rRNA variable regions varies strongly across plant-associated bacteria

Plant-microbiome research plays a pivotal role in understanding the relationships between plants and their associated microbial communities, with implications for agriculture and ecosystem dynamics. Metabarcoding analysis on variable regions of the 16S ribosomal RNA (rRNA) gene remains the dominant technology to study microbiome diversity in this field. However, the choice of the targeted variable region might affect the outcome of the microbiome studies. In our in silico analysis, we have evaluated whether the targeted variable region has an impact on taxonomic resolution in 16 plant-related microbial genera. Through a comparison of 16S rRNA gene variable regions with whole-genome data, our findings suggest that the V1-V3 region is generally a more suitable option than the widely used V3-V4 region for targeting microbiome analysis in plant-related genera. However, sole reliance on one region could introduce detection biases for specific genera. Thus, we are suggesting that while transitioning to full-length 16S rRNA gene and whole-genome sequencing for plant-microbiome analysis, the usage of genus-specific variable regions can achieve more precise taxonomic assignments. More broadly, our approach provides a blueprint to identify the most discriminating variable regions of the 16S rRNA gene for genera of interest.

A Class 4-like Chromosomal Integron Found in Aeromonas sp. Genomospecies paramedia Isolated from Human Feces

Integrons are genetic elements that store, express and exchange gene cassettes. These elements are characterized by containing a gene that codes for an integrase (intI), a cassette integration site (attI) and a variable region holding the cassettes. Using bioinformatics and molecular biology methods, a functional integron found in Aeromonas sp. 3925, a strain isolated from diarrheal stools, is described. To confirm the integron class, a phylogenetic analysis with amino acid sequences was conducted. The integrase was associated to class 4 integrases; however, it is clearly different from them. Thus, we classified the associated element as a class 4-like integron. We found that the integrase activity is not under the control of the SOS or catabolic repression, since the expression was not increased in the presence of mitomycin or arabinose. The class-4-like integron is located on the chromosome and contains two well-defined gene cassettes: aadA1 that confers resistance to streptomycin and lpt coding for a lipoprotein. It also includes eight Open Reading frames (ORFs) with unknown functions. The strain was characterized through a Multilocus Phylogenetic Analyses (MLPA) of the gyrB, gyrA, rpoD, recA, dnaJ and dnaX genes. The phylogenetic results grouped it into a different clade from the species already reported, making it impossible to assign a species. We resorted to undertaking complete genome sequencing and a phylogenomic analysis. Aeromonas sp. 3925 is related to A. media and A. rivipollensis clusters, but it is clearly different from these species. In silico DNA-DNA hybridization (isDDH) and Average Nucleotide Identity (ANI) analyses suggested that this isolate belongs to the genomospecies paramedia. This paper describes the first class 4-like integron in Aeromonas and contributes to the establishment of genomospecies paramedia.

Pontibacter kalidii sp. nov., isolated from rhizosphere soil of Kalidium foliatum

A Gram-negative, aerobic, gliding motile, rod-shaped bacterium, designated XAAS-72T, was isolated from the rhizosphere soil of Kalidium foliatum sampled in the Xinjiang Uyghur Autonomous Region, PR China. Cells grew at 4-45 °C, pH 5.0-8.0 and 0-8% NaCl, with optimal growth at 20-30 °C, pH 6.0-7.0 and 1-2 % NaCl. Strain XAAS-72T is closely related to members of the genus Pontibacter, namely Pontibacter korlensis CCTCC AB 206081T (97.6%) and Pontibacter flavimaris ACCC 19859T (97.2 %), and <94.6 % related to other currently described Pontibacter strains. The average nucleotide identity values between XAAS-72T and P. korlensis CCTCC AB 206081T and P. flavimaris ACCC 19859T were 77.9 and 86.9 %, respectively; the corresponding digital DNA-DNA hybridization values were 21.7 and 31.8 %. Menaquinone-7 was the predominant respiratory menaquinone. The polar lipids consisted of phosphatidylethanolamine, two unidentified aminophospholipids, two unidentified glycolipids and five unidentified lipids. The major cellular fatty acids were summed feature 4 (containing iso-C17 : 1 I/anteiso-C17 : 1 B), summed feature 3 (containing C16 : 1  ω7c/C16 : 1  ω6c) and iso-C15 : 0. The genome length of strain XAAS-72T was 5 054 860 bp with a genomic DNA G+C content of 54.5 mol%. The phenotypic and genotypic data suggest that strain XAAS-72T represents a novel species of the genus Pontibacter, for which the name Pontibacter kalidii sp. nov. is proposed. The strain is XAAS-72T (CGMCC 16594T=KCTC 72095T).

Rhodoalgimonas zhirmunskyi gen. nov., sp. nov., a Marine Alphaproteobacterium Isolated from the Pacific Red Alga Ahnfeltia tobuchiensis: Phenotypic Characterization and Pan-Genome Analysis

A novel Gram-staining negative, strictly aerobic, rod-shaped, and non-motile bacterium, designated strain 10Alg 79T, was isolated from the red alga Ahnfeltia tobuchiensis. A phylogenetic analysis based on 16S rRNA gene sequences placed the novel strain within the family Roseobacteraceae, class Alphaproteobacteria, phylum Pseudomonadota, where the nearest neighbor was Shimia sediminis ZQ172T (97.33% of identity). However, a phylogenomic study clearly showed that strain 10Alg 79T forms a distinct evolutionary lineage at the genus level within the family Roseobacteraceae combining with strains Aquicoccus porphyridii L1 8-17T, Marimonas arenosa KCTC 52189T, and Lentibacter algarum DSM 24677T. The ANI, AAI, and dDDH values between them were 75.63-78.15%, 67.41-73.08%, and 18.8-19.8%, respectively. The genome comprises 3,754,741 bp with a DNA GC content of 62.1 mol%. The prevalent fatty acids of strain 10Alg 79T were C18:1 ω7c and C16:0. The polar lipid profile consisted of phosphatidylethanolamine, phosphatidylglycerol, phosphatidylcholine, an unidentified aminolipid, an unidentified phospholipid and an unidentified lipid. A pan-genome analysis showed that the unique part of the 10Alg 79T genome consists of 13 genus-specific clusters and 413 singletons. The annotated singletons were more often related to transport protein systems, transcriptional regulators, and enzymes. A functional annotation of the draft genome sequence revealed that this bacterium could be a source of a new phosphorylase, which may be used for phosphoglycoside synthesis. A combination of the genotypic and phenotypic data showed that the bacterial isolate represents a novel species and a novel genus, for which the name Rhodoalgimonas zhirmunskyi gen. nov., sp. nov. is proposed. The type strain is 10Alg 79T (=KCTC 72611T = KMM 6723T).

Alteromonas salexigens sp. nov., isolated from coastal seawater

A Gram-negative, aerobic, short rod-shaped bacterium, designated ASW11-19T, was isolated from a coastal seawater sample of the Yellow Sea, PR China. Strain ASW11-19T grew optimally at 37 °C, 3.0-5.0% (w/v) NaCl and pH 7.5. Phylogenetic analysis based on the 16S rRNA gene sequences revealed that strain ASW11-19T belonged to the genus Alteromonas and most closely related to Alteromonas profundi 345S023T and Alteromonas fortis 1T (98.4%, both). The draft genome was 3.55 Mb with 3150 protein-coding genes, 18 contigs, and a DNA G+C content was 44.4%. The digital DNA-DNA hybridization and average nucleotide identity values were below the species-delineating thresholds. The major fatty acids were summed feature 3 (C16:1ω7c/C16:1ω6c), summed feature 8 (C18:1ω7c/C18:1ω6c), and C16:0. The sole respiratory quinone was ubiquinone 8. The polar lipids were phosphatidylethanolamine, phosphatidylglycerol, phospholipid, and two unidentified lipids. Based on these genomic data, phenotypic and chemotaxonomic properties, strain ASW11-19T is considered to represent a novel species of the genus Alteromonas. The name Alteromonas salexigens sp.nov. is proposed for ASW11-19T (=MCCC 1K07239T=KCTC 92247T).

Whole genome sequence analysis of Aeromonas spp. isolated from ready-to-eat seafood: antimicrobial resistance and virulence factors

Aeromonas are widespread in aquatic environments and are considered emerging pathogens in humans and animals. Multidrug resistant (MDR) Aeromonas circulating in the aquatic environment and food production chain can potentially disseminate antimicrobial resistance (AMR) to humans via the foodborne route. In this study, we aimed to investigate AMR and virulence factors of 22 Aeromonas strains isolated from ready-to-eat (RTE) seafood. A multilocus phylogenetic analysis (MLPA) using the concatenated sequences of six housekeeping genes (gyrB, rpoD, gyrA, recA, dnaJ, and dnaX) in the 22 Aeromonas genomes and average nucleotide identity (ANI) analysis revealed eight different species; A. caviae, A. dhakensis, A. hydrophila, A. media, A. rivipollensis, A. salmonicida, A. bestiarum, and A. piscicola. The presence of virulence genes, AMR genes and mobile genetic elements (MGEs) in the Aeromonas genomes was predicted using different databases. Our data showed that the genes responsible for adherence and motility (Msh type IV pili, tap type IV pili, polar flagella), type II secretion system (T2SS) and hemolysins were present in all strains, while the genes encoding enterotoxins and type VI secretion system (T6SS) including major effectors were highly prevalent. Multiple AMR genes encoding β-lactamases such as cphA and bla_OXA were detected, and the distribution of those genes was species-specific. In addition, the quinolone resistance gene, qnrS2 was found in a IncQ type plasmid of the A. rivopollensis strain A539. Furthermore, we observed the co-localization of a class I integron (intl1) with two AMR genes (sul1 and aadA1), and a Tn521 transposon carrying a mercury operon in A. caviae strain SU4-2. Various MGEs including other transposons and insertion sequence (IS) elements were identified without strongly associating with detected AMR genes or virulence genes. In conclusion, Aeromonas strains in RTE seafood were potentially pathogenic, carrying several virulence-related genes. Aeromonas carrying multiple AMR genes and MGEs could potentially be involved in the dissemination and spread of AMR genes to other bacterial species residing in the same environment and possibly to humans. Considering a One-Health approach, we highlight the significance of monitoring AMR caused by Aeromonas circulating in the food chain.

ComQXPA quorum-sensing systems contribute to enhancing the protease activity of Bacillus velezensis DMB05 from fermented soybeans

Bacillus velezensis DMB05, isolated from traditionally fermented soybean, meju, exhibited no protease activity on a TSA plate containing skim milk. To shed light on the genetic background behind this phenotypic non-protease activity, we analyzed the complete genome sequence of strain DMB05 and compared it with those of two B. velezensis strains which did exhibit protease activity. Comparative genome analyses showed no significant difference in the kind or number of proteases between the genomes of the three strains and that all strains possessed the degSU two-component system involved in the gene regulation of protease. However, strain DMB05 possessed a truncated comP which is part of the comQXPA operon that regulates the expression of degQ involved in the activation of DegSU. When the entire comQXPA operon derived from DMB06 was introduced into DMB05, the recombinant expressed proteolytic activity. The results of this experimental study provide evidence for the presence of regulatory genes involved in protease activity, one of several important factors involved in fermentation.

Genomic Characterization of Aeromonas veronii Provides Insights into Taxonomic Assignment and Reveals Widespread Virulence and Resistance Genes throughout the World

Aeromonas veronii is a Gram-negative bacterial species that causes disease in fish and is nowadays increasingly recurrent in enteric infections of humans. This study was performed to characterize newly sequenced isolates by comparing them with complete genomes deposited at the NCBI (National Center for Biotechnology Information). Nine isolates from fish, environments, and humans from the São Francisco Valley (Petrolina, Pernambuco, Brazil) were sequenced and compared with complete genomes available in public databases to gain insight into taxonomic assignment and to better understand virulence and resistance profiles of this species within the One Health context. One local genome and four NCBI genomes were misidentified as A. veronii. A total of 239 virulence genes were identified in the local genomes, with most encoding adhesion, motility, and secretion systems. In total, 60 genes involved with resistance to 22 classes of antibiotics were identified in the genomes, including mcr-7 and cphA. The results suggest that the use of methods such as ANI is essential to avoid misclassification of the genomes. The virulence content of A. veronii from local isolates is similar to those complete genomes deposited at the NCBI. Genes encoding colistin resistance are widespread in the species, requiring greater attention for surveillance systems.

Pathogenicity and Genomic Characterization of a Novel Genospecies, Bacillus shihchuchen, of the Bacillus cereus Group Isolated from Chinese Softshell Turtle (Pelodiscus sinensis)

The Chinese softshell turtle (CST; Pelodiscus sinensis) is a freshwater aquaculture species of substantial economic importance that is commercially farmed across Asia, particularly in Taiwan. Although diseases caused by the Bacillus cereus group (Bcg) pose a major threat to commercial CST farming systems, information regarding its pathogenicity and genome remains limited. Here, we investigated the pathogenicity of Bcg strains isolated in a previous study and performed whole-genome sequencing. Pathogenicity analysis indicated that QF108-045 isolated from CSTs caused the highest mortality rate, and whole-genome sequencing revealed that it was an independent group distinct from other known Bcg genospecies. The average nucleotide identity compared to other known Bcg genospecies was below 95%, suggesting that QF108-045 belongs to a new genospecies, which we named Bacillus shihchuchen. Furthermore, genes annotation revealed the presence of anthrax toxins, such as edema factor and protective antigen, in QF108-045. Therefore, the biovar anthracis was assigned, and the full name of QF108-045 was Bacillus shihchuchen biovar anthracis. In addition to possessing multiple drug-resistant genes, QF108-045 demonstrated resistance to various types of antibiotics, including penicillins (amoxicillin and ampicillin), cephalosporins (ceftifour, cephalexin, and cephazolin), and polypeptides, such as vancomycin.

Comparative genomic analysis provides insights into taxonomy and temperature adaption of Aeromonas salmonicida

Aeromonas salmonicida has long been known as psychrophiles since it is mainly isolated from cold water fish, and recent reports have revealed the existence of mesophilic strains isolated from warm sources. However, the genetic differences between mesophilic and psychrophilic strains remain unclear due to few complete genomes of mesophilic strain are available. In this study, six A. salmonicida (2 mesophilic and 4 psychrophilic) were genome-sequenced, and comparative analyses of 25 A. salmonicida complete genomes were conducted. The ANI values and phylogenetic analysis revealed that 25 strains formed three independent clades, which were referred as typical psychrophilic, atypical psychrophilic and mesophilic groups. Comparative genomic analysis showed that two chromosomal gene clusters, related to lateral flagella and outer membrane proteins (A-layer and T2SS proteins), and insertion sequences (ISAs4, ISAs7 and ISAs29) were unique to the psychrophilic groups, while the complete MSH type IV pili were unique to the mesophilic group, all of which may be considered as lifestyle-related factors. The results of this study not only provide new insights into the classification, lifestyle adaption and pathogenic mechanism of different strains of A. salmonicida, but also contributes to the prevention and control of disease caused by psychrophilic and mesophilic A. salmonicida.

A Global Survey of Hypervirulent Aeromonas hydrophila (vAh) Identified vAh Strains in the Lower Mekong River Basin and Diverse Opportunistic Pathogens from Farmed Fish and Other Environmental Sources

Hypervirulent Aeromonas hydrophila (vAh) has emerged as the etiologic agent of epidemic outbreaks of motile Aeromonas septicemia (MAS) in high-density aquaculture of farmed carp in China and catfish in the United States, which has caused millions of tons of lost fish. We conducted a global survey to better understand the evolution, geographical distribution, and phylogeny of vAh. Aeromonas isolates were isolated from fish that showed clinical symptoms of MAS, and pure cultures were screened for the ability to utilize myo-inositol as the sole carbon source. A total of 113 myo-inositol-utilizing bacterial strains were included in this study, including additional strains obtained from previously published culture collections. Based on a gyrB phylogeny, this collection included 66 A. hydrophila isolates, 48 of which were vAh. This collection also included five new vAh isolates from diseased Pangas catfish (Pangasius pangasius) and striped catfish (Pangasianodon hypophthalmus) obtained in Cambodia and Vietnam, respectively. Genome sequences were generated from representative vAh and non-vAh isolates to evaluate the potential for lateral genetic transfer of the myo-inositol catabolism pathway. Phylogenetic analyses of each of the nine genes required for myo-inositol utilization revealed the close affiliation of vAh strains regardless of geographic origin and suggested lateral genetic transfer of this catabolic pathway from an Enterobacter species. Prediction of virulence factors was conducted to determine differences between vAh and non-vAh strains in terms of virulence and secretion systems. Core genome phylogenetic analyses on vAh isolates and Aeromonas spp. disease isolates (55 in total) were conducted to evaluate the evolutionary relationships among vAh and other Aeromonas sp. isolates, which supported the clonal nature of vAh isolates. IMPORTANCE This global survey of vAh brought together scientists that study fish disease to evaluate the evolution, geographical distribution, phylogeny, and hosts of vAh and other Aeromonas sp. isolates. In addition to vAh isolates from China and the United States, four new vAh isolates were isolated from the lower Mekong River basin in Cambodia and Vietnam, indicating the significant threat of vAh to modern aquaculture and the need for improved biosecurity to prevent vAh spread.

Comparative genomics and informational content analysis uncovered internal regions of the core genes rpoD, pepN and gltX for an MLSA with genome-level resolving power within the genus Pseudomonas

In the field of prokaryotic taxonomy, there has been a recent transition towards phylogenomics as the gold standard approach. However, genome-based phylogenetics is still restrictive for its cost when managing large amounts of isolates. Fast, cheap, and taxonomically competent alternatives, like multilocus sequence analysis (MLSA) are thus recommendable. Nevertheless, the criteria for selecting the conserved genes for MLSA have not been explicit for different bacterial taxa, including the broadly diverse Pseudomonas genus. Here, we have carried out an unbiased and rational workflow to select internal sequence regions of Pseudomonas core genes (CG) for a MLSA with the best phylogenetic power, and with a resolution comparable to the genome-based ANI approach. A computational workflow was established to inspect 126 complete genomes of representatives from over 60 Pseudomonas species and subspecies, in order to identify the most informative CG internal regions and determine which combinations in sets of three partial CG sequences have comparable phylogenetic resolution to that of the current ANI standard. We found that the rpoD^346-1196-pepN^1711-2571-gltX^86-909 concatenated sequences were the best performing in terms of phylogenetic robustness and resulted highly sensitive and specific when contrasted with ANI. The rpoD-pepN-gltX MLSA was validated in silico and in vitro. Altogether, the results presented here supports the proposal of the rpoD-pepN-gltX MLSA as a fast, affordable, and robust phylogenetic tool for members of the Pseudomonas genus.

First Record of the Rare Species Aeromonas lusitana from Rainbow Trout (Oncorhynchus mykiss, Walbaum): Comparative Analysis with the Existing Strains

The species Aeromonas lusitana was first described in 2016 with five strains recovered from untreated water and vegetables from Portugal. Since then, no further records exist of this species. During a surveillance study on the presence of Aeromonas in fish farms in Mexico, a new strain (ESV-351) of the mentioned species isolated from a rainbow trout was recovered. It was identified because it clustered phylogenetically with the type strain of A. lusitana based on the analysis of the rpoD gene sequences. In the present study, phenotypic characteristics, antimicrobial resistance profiles, and the presence of putative virulence genes of this novel strain (ESV-351) were determined in parallel to the five isolates from the original species description. Phenotypic differential characteristics exhibited by A. lusitana ESV-351 depicted an evident similarity to the characteristics exhibited by the other evaluated strains. However, the novel strain was positive for the production of indole using conventional methods, while the rest of the strains, including the type strain, were negative for its production. Furthermore, intermediate resistance to ampicillin, amoxicillin-clavulanic acid and cephalothin was detected in both the novel and the type strain. Five different virulence-related genes were detected in the novel strain and in the previously described strains, with the type strain exhibiting the highest number of virulence-related genes. In addition to this, the genome of the novel strain (ESV-351) was sequenced and compared with the genomes from the type strain (A. lusitana CECT 7828T) and other Aeromonas spp. The genomic analysis defined Aeromonas tecta as the closest species to A. lusitana with a highly similar number of predicted proteins. The genomic size, the number of protein-encoding genes and the number of different tRNAs, among other characteristics, make it possible to propose that the ESV-351 strain could potentially have the capacity to adapt to different environments. Genome comparison of the ESV-351 strain with the type strain revealed that both possess a similar sequence of the citrate synthase gene. In addition to this finding, the chromosomal region containing the citrate synthase locus of the novel strain exhibits some similarity to the chromosomal region in the genome of the A. hydrophila type strain and other known human pathogens, such as Vibrio cholerae. This could suggest a possible virulence role for the citrate synthase gene in A. lusitana (ESV-351).

Draft genome sequence of a new carotenoid-producing strain Brevibacterium sp. XU54, isolated from radioactive soil in Xinjiang, China

Some species of the genus Brevibacterium are orange bacteria involved in cheese ripening, synthesis of odoriferous compounds, and carotenoids with aromatic end groups. Here, we report the genome sequence of Brevibacterium sp. XU54, isolated from radioactive soil in Xinjiang, China. The genome of XU54 consists of 4,899,099 base pairs with a GC content of 62.2%. The genome sequence was annotated with 4453 genes, encoding 4260 proteins, 13 rRNAs, and 49 tRNAs. 16S rRNA BLAST and comparative genomic analysis both indicated that XU54 may be a new species of Brevibacterium. In addition, compared to the type strains, some enzymes related to sulfur metabolism showed a low similarity of 66.85, 79.53 and 14.61%, respectively. The carotenoids biosynthesis gene cluster was identified and analyzed according to the genomic data, which revealed relatively low identity (5-85%) with existing strains. The optimum conditions for its growth and carotenoid production were then discussed. The whole-genome sequence of Brevibacterium sp. XU54 will be beneficial for utilizing these newly identified genes in carotenoid biosynthesis and regulation of sulfur metabolism pathway to promote the production of novel carotenoids and other structurally diverse compounds through combinatorial biosynthesis, which facilitates cheese ripening and coloration.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03366-1.

Aeromonas allosaccharophila Strain AE59-TE2 Is Highly Antagonistic towards Multidrug-Resistant Human Pathogens, What Does Its Genome Tell Us?

Multidrug-resistant bacteria are of critical importance and a problem for human health and food preservation; the discovery of new antimicrobial substances to control their proliferation is part of the solution. This work reports on 57 antagonistic Aeromonas strains, of which 38 strains were antagonistic towards problematic human pathogens. The genome of the most antagonistic strain was sequenced and identified as Aeromonas allosaccharophila. Its genome was fully annotated and mined for genes that might explain that activity. Strain AE59-TE was antagonistic toward clinically relevant gram-negative and gram-positive multidrug-resistant bacteria, including Klebsiella pneumoniae KPC, Escherichia coli ESBL, Salmonella typhimurium, and Staphylococcus aureus MRSA. Strain AE59-TE2 was identified by multilocus sequence analysis. Genome mining identified four genes homologous to the bacteriocin, zoocin A from Streptococcus equi and a gene 98% similar to cvpA linked to colicin V production. A. allosaccharophila strain AE59-TE2 produced antimicrobial activity against a broad range of bacteria, including important gram-negative bacteria, not typically targeted by bacteriocins. Herewere described novel zoocin genes that are promising for industrial applications in the food and health sectors. Interesting and important antagonistic activity is described combined with the first detailed genomic analysis of the species Aeromonas allosaccharophila.

Genome-based analyses reveal heterotypic synonyms of Streptomyces species and associated subspecies

In the genus Streptomyces, several validly described species have been reduced to synonyms of earlier described species though additional synonyms remain to be detected given the previous dependence on traditional phenotypic methods. In this study, genome-based procedures, including DNA-DNA hybridization analyses, overall genome-related indices, such as ANI, dDDH and AAI, revealed that certain strains recorded genomic indices above the threshold values used to define species boundaries. The results of phylogenetic and phylogenomic trees based on concatenated and phylogenomic analyses showed that 33 out of 364 tested species could be assigned to 15 species groups and that 18 Streptomyces species names be reclassified as later heterotypic synonyms of earlier validly published species. Consequently, it is proposed that S. albaduncus is a later heterotypic synonym of S. griseoloalbus; S. bellus is a synonym of S. coeruleorubidus; S. gancidicus and S. rubiginosus are synonyms of S. pseudogriseolus; S. niveoruber is a synonym of S. griseoviridis; S. griseomycini is a synonym of S. griseostramineus; S. jietaisiensis is a synonym of S. griseoaurantiacus; S. pluricolorescens is a synonym of S. rubiginosohelvolus; S. nashvillensis is a synonym of S. tanashiensis; S. yerevanensis is a synonym of S. flaveus; S. durhamensis is a synonym of S. filipinensis; S. recifensis is a synonym of S. griseoluteus; S. canaries and S. olivaceoviridis are synonyms of S. corchorusii; S. melanosporofaciens is a synonym of S. antimycoticus; S. albulus is a synonym of S. noursei; and S. janthinus and S. violarus are synonyms of S. violaceus. Additionally, seven of these 18 Streptomyces species have been designated subspecies.

Algicella marina gen. nov., sp. nov., a novel marine bacterium isolated from a Pacific red alga

A novel Gram-staining negative, strictly aerobic, rod-shaped, and non-motile bacterium, designated strain 9Alg 56^T, was isolated from the red alga Tichocarpus crinitus. The phylogenetic analysis based on 16S rRNA gene sequences placed the novel strain within the family Rhodobacteraceae, the order Rhodobacterales, the class Alphaproteobacteria, the phylum Pseudomonadota. The nearest neighbors of the new strain were Pontivivens insulae KCTC 42458^T, Oceanibium sediminis KCTC 62076^T, Halovulum dunhuangense YYQ-30^T and Monaibacterium marinum C7^T with 16S rRNA gene sequence similarity of 94.7, 94.4%, 93.1 and 92.7%, respectively. The AAI/ANI/dDDH values between 9Alg 56^T and the five species of the closest genera (Pontivivens, Oceanibium, Halovulum, Monaibacterium, and 'Oceanomicrobium') were 58.63-63.91%/ 75.91-77.37%/ 19.3-20.4%. The prevalent fatty acids of strain 9Alg 56^T were C_18:1 ω7c, C_18:0 and C_14:0 3-OH. The polar lipid profile consisted of phosphatidylethanolamine, phosphatidylglycerol, phosphatidic acid, phosphatidylcholine, and two unidentified lipids. The DNA G+C content of strain 9Alg 56^T was 61.5 mol%. A combination of the genotypic and phenotypic data showed that the algal isolate represents a novel genus and species, for which the name Algicella marina gen. nov., sp. nov. is proposed. The type strain is 9Alg 56^T (= KCTC 72005^T = KMM 6775^T).

Comparative analysis of genome-based CAZyme cassette in Antarctic Microbacterium sp. PAMC28756 with 31 other Microbacterium species

BACKGROUND

The genus Microbacterium belongs to the family Microbacteriaceae and phylum Actinobacteria. A detailed study on the complete genome and systematic comparative analysis of carbohydrate-active enzyme (CAZyme) among the Microbacterium species would add knowledge on metabolic and environmental adaptation. Here we present the comparative genomic analysis of CAZyme using the complete genome of Antarctic Microbacterium sp. PAMC28756 with other complete genomes of 31 Microbacterium species available.

OBJECTIVE

The genomic and CAZyme comparison of Microbacterium species and to rule out the specific features of CAZyme for the environmental and metabolic adaptation.

METHODS

Bacterial source were collected from NCBI database, CAZyme annotation of Microbacterium species was analyzed using dbCAN2 Meta server. Cluster of orthologous groups (COGs) analysis was performed using the eggNOG4.5 database. Whereas, KEGG database was used to compare and obtained the functional genome annotation information in carbohydrate metabolism and glyoxylate cycle.

RESULTS

Out of 32 complete genomes of Microbacterium species, strain No. 7 isolated from Activated Sludge showed the largest genomic size at 4.83 Mb. The genomic size of PAMC28756 isolated from Antarctic lichen species Stereocaulons was 3.54 Mb, the G + C content was 70.4% with 3,407 predicted genes, of which 3.36% were predicted CAZyme. In addition, while comparing the Glyoxylate cycle among 32 bacteria, except 10 strains, all other, including our strain have Glyoxylate pathway. PAMC28756 contained the genes that degrade cellulose, hemicellulose, amylase, pectinase, chitins and other exo-and endo glycosidases. Utilizing these polysaccharides can provides source of energy in an extreme environment. In addition, PAMC28756 assigned the (10.15%) genes in the carbohydrate transport and metabolism functional group closely related to the CAZyme for polysaccharides degradation.

CONCLUSIONS

The genomic content and CAZymes distribution was varied in Microbacterium species. There was the presence of more than 10% genes in the carbohydrate transport and metabolism functional group closely related to the CAZyme for polysaccharides degradation. In addition, occurrence of glyoxylate cycle for alternative utilization of carbon sources suggest the adaptation of PAMC28756 in the harsh microenvironment.

Establishment of Genome Based Criteria for Classification of the Family Desulfovibrionaceae and Proposal of Two Novel Genera, Alkalidesulfovibrio gen. nov. and Salidesulfovibrio gen. nov

Bacteria in the Desulfovibrionaceae family, which contribute to S element turnover as sulfate-reducing bacteria (SRB) and disproportionation of partially oxidized sulfoxy anions, have been extensively investigated since the importance of the sulfur cycle emerged. Novel species belonging to this taxon are frequently reported, because they exist in various environments and are easy to culture using established methods. Due to the rapid expansion of the taxon, correction and reclassification have been conducted. The development of high-throughput sequencing facilitated rapid expansion of genome sequence database. Genome-based criteria, based on these databases, proved to be potential classification standard by overcoming the limitations of 16S rRNA-based phylogeny. Although standards methods for taxogenomics are being established, the addition of a novel genus requires extensive calculations with taxa, including many species, such as Desulfovibrionaceae. Thus, the genome-based criteria for classification of Desulfovibrionaceae were established and validated in this study. The average amino-acid identity (AAI) cut-off value, 63.43 ± 0.01, was calculated to be an appropriate criterion for genus delineation of the family Desulfovibrionaceae. By applying the AAI cut-off value, 88 genomes of the Desulfovibrionaceae were divided into 27 genera, which follows the core gene phylogeny results. In this process, two novel genera (Alkalidesulfovibrio and Salidesulfovibrio) and one former invalid genus (“Psychrodesulfovibrio”) were officially proposed. Further, by applying the 95–96% average nucleotide identity (ANI) standard and the 70% digital DNA–DNA hybridization standard values for species delineation of strains that were classified as the same species, five strains have the potential to be newly classified. After verifying that the classification was appropriately performed through relative synonymous codon usage analysis, common characteristics were listed by group. In addition, by detecting metal resistance related genes via in silico analysis, it was confirmed that most strains display metal tolerance.

Comparative Genomic Analysis Reveals Potential Pathogenicity and Slow-Growth Characteristics of Genus Brevundimonas and Description of Brevundimonas pishanensis sp. nov

The genus Brevundimonas consists of Gram-negative bacteria widely distributed in environment and can cause human infections. However, the genomic characteristics and pathogenicity of Brevundimonas remain poorly studied. Here, the whole-genome features of 24 Brevundimonas type strains were described. Brevundimonas spp. had relatively small genomes (3.13 ± 0.29 Mb) within the family Caulobacteraceae but high G+C contents (67.01 ± 2.19 mol%). Two-dimensional hierarchical clustering divided those genomes into 5 major clades, in which clades II and V contained nine and five species, respectively. Interestingly, phylogenetic analysis showed a one-to-one match between core and accessory genomes, which suggested coevolution of species within the genus Brevundimonas. The unique genes were annotated to biological functions like catalytic activity, signaling and cellular processes, multisubstance metabolism, etc. The majority of Brevundimonas spp. harbored virulence-associated genes icl, tufA, kdsA, htpB, and acpXL, which encoded isocitrate lyase, elongation factor, 2-dehydro-3-deoxyphosphooctonate aldolase, heat shock protein, and acyl carrier protein, respectively. In addition, genomic islands (GIs) and phages/prophages were identified within the Brevundimonas genus. Importantly, a novel Brevundimonas species was identified from the feces of a patient (suffering from diarrhea) by the analyses of biochemical characteristics, phylogenetic tree of 16S rRNA gene, multilocus sequence analysis (MLSA) sequences, and genomic data. The name Brevundimonas pishanensis sp. nov. was proposed, with type strain CHPC 1.3453 (= GDMCC 1.2503^T= KCTC 82824^T). Brevundimonas spp. also showed obvious slow growth compared with that of Escherichia coli. Our study reveals insights into genomic characteristics and potential virulence-associated genes of Brevundimonas spp., and provides a basis for further intensive study of the pathogenicity of Brevundimonas.

IMPORTANCEBrevundimonas spp., a group of bacteria from the family Caulobacteraceae, is associated with nosocomial infections, deserve widespread attention. Our study elucidated genes potentially associated with the pathogenicity of the Brevundimonas genus. We also described some new characteristics of Brevundimonas spp., such as small chromosome size, high G+C content, and slow-growth phenotypes, which made the Brevundimonas genus a good model organism for in-depth studies of growth rate traits. Apart from the comparative analysis of the genomic features of the Brevundimonas genus, we also reported a novel Brevundimonas species, Brevundimonas pishanensis, from the feces of a patient with diarrhea. Our study promotes the understanding of the pathogenicity characteristics of Brevundimonas species bacteria.

Improving Phylogenies Based on Average Nucleotide Identity, Incorporating Saturation Correction and Nonparametric Bootstrap Support

Whole-genome comparisons based on average nucleotide identities (ANI) and the genome-to-genome distance calculator have risen to prominence in rapidly classifying prokaryotic taxa using whole-genome sequences. Some implementations have even been proposed as a new standard in species classification and have become a common technique for papers describing newly sequenced genomes. However, attempts to apply whole-genome divergence data to the delineation of higher taxonomic units and to phylogenetic inference have had difficulty matching those produced by more complex phylogenetic methods. We present a novel method for generating statistically supported phylogenies of archaeal and bacterial groups using a combined ANI and alignment fraction-based metric. For the test cases to which we applied the developed approach, we obtained results comparable with other methodologies up to at least the family level. The developed method uses nonparametric bootstrapping to gauge support for inferred groups. This method offers the opportunity to make use of whole-genome comparison data, that is already being generated, to quickly produce phylogenies including support for inferred groups. Additionally, the developed ANI methodology can assist the classification of higher taxonomic groups.[Average nucleotide identity (ANI); genome evolution; prokaryotic species delineation; taxonomy.].

To Be or Not to Be Mesophilic, That Is the Question for Aeromonas salmonicida

The bacterium Aeromonas salmonicida has long been known to be one of the most feared pathogens in fish farming. However, the more we discover about this bacterial species, the more we question whether it is really exclusively an aquatic pathogen. In recent years, it has become obvious that this bacterial species includes a myriad of strains with various lifestyle and ecological niches, including the well-known strict psychrophiles, the first bacteria known of the species, and the newly described mesophilic strains. The mesophiles are able to grow at low temperatures, but even better at temperatures of approximately 37 °C, which strict psychrophiles cannot do. In this perspective article, we address some aspects surrounding this dual lifestyle in A. salmonicida, including the impact of mobile genetic elements, and how future research around this bacterial species may focus on the psychrophilic/mesophilic dichotomy, which makes A. salmonicida an increasingly interesting and relevant model for the study of speciation.

Comparative Genomics of Typical and Atypical Aeromonas salmonicida Complete Genomes Revealed New Insights into Pathogenesis Evolution

Aeromonas salmonicida is a global distributed Gram-negative teleost pathogen, affecting mainly salmonids in fresh and marine environments. A. salmonicida strains are classified as typical or atypical depending on their origin of isolation and phenotype. Five subspecies have been described, where A. salmonicida subsp. salmonicida is the only typical subspecies, and the subsp. achromogenes, masoucida, smithia, and pectinolytica are considered atypical. Genomic differences between A. salmonicida subsp. salmonicida isolates and their relationship with the current classification have not been explored. Here, we sequenced and compared the complete closed genomes of four virulent strains to elucidate their molecular diversity and pathogenic evolution using the more accurate genomic information so far. Phenotypes, biochemical, and enzymatic profiles were determined. PacBio and MiSeq sequencing platforms were utilized for genome sequencing. Comparative genomics showed that atypical strains belong to the subsp. salmonicida, with 99.55% ± 0.25% identity with each other, and are closely related to typical strains. The typical strain A. salmonicida J223 is closely related to typical strains, with 99.17% identity with the A. salmonicida A449. Genomic differences between atypical and typical strains are strictly related to insertion sequences (ISs) activity. The absence and presence of genes encoding for virulence factors, transcriptional regulators, and non-coding RNAs are the most significant differences between typical and atypical strains that affect their phenotypes. Plasmidome plays an important role in A. salmonicida virulence and genome plasticity. Here, we determined that typical strains harbor a larger number of plasmids and virulence-related genes that contribute to its acute virulence. In contrast, atypical strains harbor a single, large plasmid and a smaller number of virulence genes, reflected by their less acute virulence and chronic infection. The relationship between phenotype and A. salmonicida subspecies’ taxonomy is not evident. Comparative genomic analysis based on completed genomes revealed that the subspecies classification is more of a reflection of the ecological niche occupied by bacteria than their divergences at the genomic level except for their accessory genome.

Comparative Genomics of Pseudomonas stutzeri Complex: Taxonomic Assignments and Genetic Diversity

Pseudomonas stutzeri is a species complex with extremely broad phenotypic and genotypic diversity. However, very little is known about its diversity, taxonomy and phylogeny at the genomic scale. To address these issues, we systematically and comprehensively defined the taxonomy and nomenclature for this species complex and explored its genetic diversity using hundreds of sequenced genomes. By combining average nucleotide identity (ANI) evaluation and phylogenetic inference approaches, we identified 123 P. stutzeri complex genomes covering at least six well-defined species among all sequenced Pseudomonas genomes; of these, 25 genomes represented novel members of this species complex. ANI values of ≥∼95% and digital DNA-DNA hybridization (dDDH) values of ≥∼60% in combination with phylogenomic analysis consistently and robustly supported the division of these strains into 27 genomovars (most likely species to some extent), comprising 16 known and 11 unknown genomovars. We revealed that 12 strains had mistaken taxonomic assignments, while 16 strains without species names can be assigned to the species level within the species complex. We observed an open pan-genome of the P. stutzeri complex comprising 13,261 gene families, among which approximately 45% gene families do not match any sequence present in the COG database, and a large proportion of accessory genes. The genome contents experienced extensive genetic gain and loss events, which may be one of the major mechanisms driving diversification within this species complex. Surprisingly, we found that the ectoine biosynthesis gene cluster (ect) was present in all genomes of P. stutzeri species complex strains but distributed at very low frequency (43 out of 9548) in other Pseudomonas genomes, suggesting a possible origin of the ancestors of P. stutzeri species complex in high-osmolarity environments. Collectively, our study highlights the potential of using whole-genome sequences to re-evaluate the current definition of the P. stutzeri complex, shedding new light on its genomic diversity and evolutionary history.

Pancreaticobiliary Cancers and Aeromonas Isolates Carrying Type Ⅲ Secretion System Genes ascF-ascG Are Associated With Increased Mortality: An Analysis of 164 Aeromonas Infection Episodes in Southern Taiwan

This prospective study aimed to investigate the clinical and microbiological characteristics of different Aeromonas species. Clinical isolates of Aeromonas species between 2016 to 2018 were collected in a university hospital in southern Taiwan. The species was determined by rpoD or gyrB sequencing. A total of 222 Aeromonas isolates from 160 patients in 164 episodes were identified. The crude in-hospital mortality was 17.2%. The most frequently isolated species was Aeromonas veronii (30.6%), followed by A. caviae (24.8%), A. hydrophila (23%), and A. dhakensis (16.7%). The major clinical manifestations were primary bacteremia (31.1%), skin and soft tissue infection (22.6%), and biliary tract infection (18.3%). The most common underlying diseases were malignancy (45.1%), diabetes mellitus (27.4%), and liver cirrhosis or chronic hepatitis (26.2%). A. hydrophila and A. dhakensis predominated in the skin and soft tissue infection (p<0.0001), whereas A. vernoii and A. caviae prevailed in primary bacteremia and biliary tract infections (p=0.012). Pneumonia, malignancy, and ascF-ascG genotype were independent factors associated with mortality. Ertapenem susceptibility was decreased in A. sobria (42.9%), A. veronii (66.7%), A. dhakensis (73%), and A. hydrophila (84.3%). Cefotaxime resistance was found in 30.9% of A. caviae and 18.9% of A. dhakensis isolates, much more prevalent than the other species. The metallo-β-lactamase blaCphA was almost invariably present in A. dhakensis, A. hydrophila, and A. veronii (100%, 100% and 89.9%, respectively). Amp-C β-lactamases such as blaMOX and blaAQU-1 were identified in all A. caviae and 91.9% of A. dhakensis isolates. Cefepime, fluoroquinolones and tigecycline showed good in vitro activity against aeromonads.

Novosphingobium decolorationis sp. nov., an aniline blue-decolourizing bacterium isolated from East Pacific sediment

Aniline blue-decolourizing bacterial strain 502str22^T, isolated from sediment collected in the East Pacific, was subjected to characterization by a polyphasic taxonomic approach. Phylogenetic analyses based on 16S rRNA gene sequences revealed that strain 502str22^T belongs to the genus Novosphingobium, with closely related type strains 'Novosphingobium profundi' F72^T (97.6%), N. mathurense SM117^T (97.1%) and N. arvoryzae Jyi-02^T (97.0%). Digital DNA-DNA hybridization and average nucleotide identity values between strain 502str22^T and closely related type strains were 20.3-24.8% and 74.1-81.9%, respectively. The major cellular fatty acid (>10%) was C_18:1 ω7c. The polar lipid profile consisted of a mixture of phosphatidylcholine, one sphingoglycolipid, diphosphatidylglycerol, phosphatidylglycerol, phosphatidylethanolamine and phosphatidylmonomethylethanolamine. The DNA G+C content of strain 502str22^T was 65.5 mol%. The polyphasic taxonomic results indicated that strain 502str22^T represents a novel species of the genus Novosphingobium, for which the name Novosphingobium decolorationis sp. nov is proposed. The type strain is 502str22^T (=KCTC 82134^T= MCCC 1K04799 ^T).

Comparative Genome Analysis Provides Molecular Evidence for Reclassification of the Photosynthetic Bacterium Rhodobacter sphaeroides EBL0706 as a Strain of Luteovulum azotoformans

In this study, we conducted a genome-wide comparative analysis of a former Rhodobacter sphaeroides strain EBL0706, which is now recorded as Luteovulum sphaeroides EBL0706. The genome of EBL0706 was compared with that of Luteovulum azotoformans ATCC 17025, Luteovulum azotoformans KA25, and Luteovulum sphaeroides 2.4.1. The average nucleotide identity (ANI), tetra nucleotide signatures (Tetra), digital DNA–DNA hybridization (dDDH) values, comparative genome, and phylogenetic analysis proposed that EBL0706 is a strain of Luteovulum azotoformans. Functional annotations identified a total of 4034 protein-coding genes in the genome of EBL0706, including a complete photosynthetic gene cluster. This study provides genomic molecular verification for the strain EBL0706 to be reclassified to Luteovulum azotoformans.

Systematic Detection of Large-Scale Multigene Horizontal Transfer in Prokaryotes

Horizontal gene transfer (HGT) is central to prokaryotic evolution. However, little is known about the “scale” of individual HGT events. In this work, we introduce the first computational framework to help answer the following fundamental question: How often does more than one gene get horizontally transferred in a single HGT event? Our method, called HoMer, uses phylogenetic reconciliation to infer single-gene HGT events across a given set of species/strains, employs several techniques to account for inference error and uncertainty, combines that information with gene order information from extant genomes, and uses statistical analysis to identify candidate horizontal multigene transfers (HMGTs) in both extant and ancestral species/strains. HoMer is highly scalable and can be easily used to infer HMGTs across hundreds of genomes. We apply HoMer to a genome-scale data set of over 22,000 gene families from 103 Aeromonas genomes and identify a large number of plausible HMGTs of various scales at both small and large phylogenetic distances. Analysis of these HMGTs reveals interesting relationships between gene function, phylogenetic distance, and frequency of multigene transfer. Among other insights, we find that 1) the observed relative frequency of HMGT increases as divergence between genomes increases, 2) HMGTs often have conserved gene functions, and 3) rare genes are frequently acquired through HMGT. We also analyze in detail HMGTs involving the zonula occludens toxin and type III secretion systems. By enabling the systematic inference of HMGTs on a large scale, HoMer will facilitate a more accurate and more complete understanding of HGT and microbial evolution.

Comparative and Evolutionary Genomics of Isolates Provide Insight into the Pathoadaptation of Aeromonas

Aeromonads are ubiquitous aquatic bacteria that cause opportunistic infections in humans, but their pathogenesis remains poorly understood. A pathogenomic approach was undertaken to provide insights into the emergence and evolution of pathogenic traits in aeromonads. The genomes of 64 Aeromonas strains representative of the whole genus were analyzed to study the distribution, phylogeny, and synteny of the flanking sequences of 13 virulence-associated genes. The reconstructed evolutionary histories varied markedly depending on the gene analyzed and ranged from vertical evolution, which followed the core genome evolution (alt and colAh), to complex evolution, involving gene loss by insertion sequence-driven gene disruption, horizontal gene transfer, and paraphyly with some virulence genes associated with a phylogroup (aer, ser, and type 3 secretion system components) or no phylogroup (type 3 secretion system effectors, Ast, ExoA, and RtxA toxins). The general pathogenomic overview of aeromonads showed great complexity with diverse evolution modes and gene organization and uneven distribution of virulence genes in the genus; the results provided insights into aeromonad pathoadaptation or the ability of members of this group to emerge as pathogens. Finally, these findings suggest that aeromonad virulence-associated genes should be examined at the population level and that studies performed on type or model strains at the species level cannot be generalized to the whole species.

Comparative Genomics Reveals Potential Mechanisms of Plant Beneficial Effects of a Novel Bamboo-Endophytic Bacterial Isolate Paraburkholderia sacchari Suichang626

Plant-beneficial microbes have drawn wide attention due to their potential application as bio-control agents and bio-fertilizers. Moso bamboo, which is among the monocots with the highest growth rate, lives perennially with abundant microbes that may benefit annually growing crops. Genome information of moso bamboo associated bacteria remains underexplored. We isolated and identified a novel Paraburkholderia strain Suichang626 from moso bamboo roots. Growth promoting effects of Suichang626 on both moso bamboo and seedlings of the model dicot Arabidopsis thaliana were documented in laboratory conditions. To gain insight into the genetic basis of this growth promotion effect, we sequenced the genome of Suichang626. Evidenced by genome-wide phylogeny data, we propose that Suichang626 is a novel strain of Paraburkholderia sacchari. Gene homologs encoding biosynthesis of the plant growth-promoting chemicals, acetoin and 2,3-butanediol, were identified in the genome of Suichang626. Comparative genomics was further performed with plant-beneficial and plant/animal pathogenic species of Paraburkholderia and Burkholderia. Genes related to volatile organic compounds, nitrogen fixation, and auxin biosynthesis were discovered specifically in the plant growth-promoting species of both genera.

Comparative Genomic Analyses Reveal Functional Insights Into Key Determinants of the Pathogenesis of Pectobacterium actinidiae in Kiwifruit

The Gram-negative bacterial species Pectobacterium actinidiae causes summer canker in kiwifruit plants. However, little is known about its virulence factors and mechanisms of genetic adaptation. We aimed to identify the key determinants that control the virulence of P. actinidiae in kiwifruit by genomic and functional analyses. Analysis of four P. actinidiae isolates indicated low genetic variability with an average of 98.7% genome-level sequence similarity and 82% shared protein-coding gene content. Phylogenetic analysis, based on both bulk single nucleotide polymorphisms (SNPs) and single-copy genes, revealed that P. actinidiae strains cluster into a single clade, which is closely related to the clades of P. odoriferum (species with a completely different host range). Through comparison between these two clades of strains, 746 unique core orthologs/genes were clustered in the clades of P. actinidiae, especially key virulence determinants involved in the biosynthesis of secretion systems (type III, IV, and VI), iron, flagellar structure, and the quorum-sensing system. Our results provide insights into the pathogenomics underlying the genetic diversification and evolution of pathogenicity in P. actinidiae species.

Muricauda sediminis sp. nov., isolated from western Pacific Ocean sediment

A Gram-stain-negative bacterium, designated strain 40Bstr401^T, was isolated from a sediment sample collected from the western Pacific Ocean. Analysis of its 16S rRNA gene sequence revealed that strain 40Bstr401^T belongs to the genus Muricauda and is closely related to type strains Muricauda antarctica Ar-22^T (98.2 %), Muricauda taeanensis 105^T (98.2 %) and Muricauda beolgyonensis BB-My12^T (97.4 %). The average nucleotide identity values for 40Bstr401^T with M. antarctica Ar-22^T and M. taeanensis 105^T are 79.3 % and 78.8 %, respectively. The in silico DNA-DNA hybridization values between strain 40Bstr401^T and M. antarctica Ar-22^T and M. taeanensis 105^T are 26.7 and 26.6 %, respectively. The major isoprenoid quinone of 40Bstr401^T is MK-6, and iso-C_17 : 0 3-OH and iso-C_15 : 0 are the dominant cellular fatty acids. The major polar lipids are phosphatidylethanolamine, four unidentified amino lipids and two unidentified lipids. The G+C content of the genomic DNA is 42.9 mol%. Its phylogenetic distinctiveness and chemotaxonomic differences, together with the phenotypic properties observed in this study, indicate that strain 40Bstr401^T can be differentiated from closely related species. Therefore, we propose strain 40Bstr401^T represents a novel species in the genus Muricauda, for which the name Muricauda sediminis sp. nov. is suggested. The type strain is 40Bstr401^T (=MCCC 1K04568^T=KCTC 82139^T).

Pelistega ratti sp. nov. from Rattus norvegicus of Hainan island

Two strains (NLN63^T and NLN82) of Gram-stain-negative, oxidase- and catalase-positive, bacilli-shaped organisms were isolated from the faecal samples of two separate Rattus norvegicus in Baisha county of Hainan Province, Southern PR China. Phylogenetic analysis based on the near full-length 16S rRNA sequences revealed that strain NLN63^T belongs to the genus Pelistega, having maximum similarity to Pelistega suis CCUG 64465^T (97.1 %), Pelistega europaea CCUG 39967^T (96.2 %) and Pelistega indica DSM 27484^T (96.2 %), respectively. The phylogenomic tree built on 553 core genes from genomes of 20 species in the genus Pelistega and other adjacent genera further confirmed that strains NLN63^T and NLN82 form a distinct subline and exhibit specific phylogenetic affinity with P. europaea CCUG 39967^T. In digital DNA-DNA hybridization analyses, strain NLN63^T showed low estimated DNA reassociation values (21.4-22.6 %) with the type strains of the species in the genus Pelistega. The DNA G+C contents of strains NLN63^T and NLN82 were 37.3 and 37.1 mol%, respectively. Strain NLN63^T had a unique MALDI-TOF MS profile, contained Q-8 as the major quinone and C_16 : 0, summed feature 8 (C_18 : 1  ω7c/C_18 : 1  ω6c or both) and summed feature 3 (C_16 : 1  ω7c/C_16 : 1  ω6c or both) as the dominant fatty acids. Based upon these polyphasic characterization data obtained from the present study, a novel species of the genus Pelistega, Pelistega ratti sp. nov., is proposed with NLN63^T (=GDMCC 1.1697^T=JCM 33788^T) as the type strain.

High-Resolution Differentiation of Enteric Bacteria in Premature Infant Fecal Microbiomes Using a Novel rRNA Amplicon

Identifying and tracking microbial strains as microbiomes evolve are major challenges in the field of microbiome research. We utilized a new sequencing kit that combines DNA extraction with PCR amplification of a large region of the rRNA operon and downstream bioinformatic data analysis. Longitudinal microbiome samples of coadmitted twins from two different neonatal intensive care units (NICUs) were analyzed using an ∼2,500-base amplicon that spans the 16S and 23S rRNA genes and mapped to a new, custom 16S-23S rRNA database. Amplicon sequence variants (ASVs) inferred using DADA2 provided sufficient resolution for the differentiation of rRNA variants from closely related but not previously sequenced Klebsiella, Escherichia coli, and Enterobacter strains, among the first bacteria colonizing the gut of these infants after admission to the NICU. Distinct ASV groups (fingerprints) were monitored between coadmitted twins over time, demonstrating the potential to track the source and spread of both commensals and pathogens. The high-resolution taxonomy obtained from long amplicon sequencing enables the tracking of strains temporally and spatially as microbiomes are established in infants in the hospital environment.

Comparative Genomics Suggests a Taxonomic Revision of the Staphylococcus cohnii Species Complex

Staphylococcus cohnii (SC), a coagulase-negative bacterium, was first isolated in 1975 from human skin. Early phenotypic analyses led to the delineation of two subspecies (subsp.), Staphylococcus cohnii subsp. cohnii (SCC) and Staphylococcus cohnii subsp. urealyticus (SCU). SCC was considered to be specific to humans, whereas SCU apparently demonstrated a wider host range, from lower primates to humans. The type strains ATCC 29974 and ATCC 49330 have been designated for SCC and SCU, respectively. Comparative analysis of 66 complete genome sequences-including a novel SC isolate-revealed unexpected patterns within the SC complex, both in terms of genomic sequence identity and gene content, highlighting the presence of 3 phylogenetically distinct groups. Based on our observations, and on the current guidelines for taxonomic classification for bacterial species, we propose a revision of the SC species complex. We suggest that SCC and SCU should be regarded as two distinct species: SC and SU (Staphylococcus urealyticus), and that two distinct subspecies, SCC and SCB (SC subsp. barensis, represented by the novel strain isolated in Bari) should be recognized within SC. Furthermore, since large-scale comparative genomics studies recurrently suggest inconsistencies or conflicts in taxonomic assignments of bacterial species, we believe that the approach proposed here might be considered for more general application.

On the transformation of MinHash-based uncorrected distances into proper evolutionary distances for phylogenetic inference

Recently developed MinHash-based techniques were proven successful in quickly estimating the level of similarity between large nucleotide sequences. This article discusses their usage and limitations in practice to approximating uncorrected distances between genomes, and transforming these pairwise dissimilarities into proper evolutionary distances. It is notably shown that complex distance measures can be easily approximated using simple transformation formulae based on few parameters. MinHash-based techniques can therefore be very useful for implementing fast yet accurate alignment-free phylogenetic reconstruction procedures from large sets of genomes. This last point of view is assessed with a simulation study using a dedicated bioinformatics tool.

Croceivirga litoralis sp. nov., isolated from coastal surface water, and reclassification of Muricauda lutea as Croceivirga lutea comb. nov

A novel Gram-stain-negative, aerobic, yellow-pigmented bacterium was isolated from seawater of Aoshan Bay, and designated as strain ASW18^T. Strain ASW18^T was a long-rod-shaped bacterium without flagellum and lacked gliding ability. Based on 16S rRNA gene phylogeny, strain ASW18^T showed the closest relationship to Croceivirga radicis MCCC 1A06690^T, with a sequence similarity of 97.0 %. Strain ASW18^T was able to grow at 25-40 °C, at pH 5.5-9.5 and with 0.5-9 % (w/v) NaCl. The genomic DNA G+C content of strain ASW18^T was 37.3 %. The predominant cellular fatty acids of strain ASW18^T were iso-C_15 : 0, iso-C_17 : 0 3-OH and iso-C_15 : 1 G. The major polar lipids were phosphatidylethanolamine, phosphatidyldimethylethanolamine, an aminolipid and three unidentified lipids. The respiratory quinone of strain ASW18^T was menaquinone with six isoprene units (MK-6). Based on the present polyphasic analysis, strain ASW18^T represents a novel species of the genus Croceivirga, for which the name Croceivirga litoralis sp. nov. is proposed; the type strain is ASW18^T (=MCCC 1K04203^T=KCTC 72852^T). In addition, it is also proposed that Muricauda lutea should be reclassified as Croceivirga lutea comb. nov.; the type strain is CSW06^T (=CGMCC 1.15761^T=JCM 31455^T=KCTC 52375^T=MCCC 1K03195^T).

Genetic relatedness and novel sequence types of clinical Aeromonas dhakensis from Malaysia

Aeromonas dhakensis is an emergent human pathogen with medical importance. This study was aimed to determine the sequence types (STs), genetic diversity, and phylogenetic relationships of different clinical sources of 47 A. dhakensis from Malaysia using multilocus sequence typing (MLST), goeBURST, and phylogenetic analyses. The analysis of a concatenated six-gene tree with a nucleotide length of 2994 bp based on six housekeeping genes (gyrB, groL, gltA, metG, ppsA, and recA) and independent analyses of single gene fragments was performed. MLST was able to group 47 A. dhakensis from our collection into 36 STs in which 34 STs are novel STs. The most abundant ST521 consisted of five strains from peritoneal fluid and two strains from stools. Comparison of 62 global A. dhakensis was carried out via goeBURST; 94.4% (34/36) of the identified STs are novel and unique in Malaysia. Two STs (111 and 541) were grouped into clonal complexes among our strains and 32 STs occurred as singletons. Single-gene phylogenetic trees showed varying topologies; groL and rpoD grouped all A. dhakensis into a tight-cluster with bootstrap values of 100% and 99%, respectively. A poor phylogenetic resolution encountered in single-gene analyses was buffered by the multilocus phylogenetic tree that offered high discriminatory power (bootstrap value = 100%) in resolving all A. dhakensis from A. hydrophila and delineating the relationship among other taxa. Genetic diversity analysis showed groL as the most conserved gene and ppsA as the most variable gene. This study revealed novel STs and high genetic diversity among clinical A. dhakensis from Malaysia.

Novel xylanase producing Bacillus strain X2: molecular phylogenetic analysis and its application for production of xylooligosaccharides

A Bacillus strain X2 that produced extracellular endo-xylanase (GH 11) (EC: 3.2.1.8) was isolated from the soil of the Northeast India region. This aerobic culture was Gram positive and endospore forming. Chemotaxonomic characterization showed variance with the fatty acid profile of related species in the Bacillus subtilis group. In Bacillus strain X2, distinct occurrence of iso-C14:0 lipids is absent in other related species. The 16S rRNA gene sequence homology showed 99% similarity with Bacillus subtilis subsp. inaquosorum. The phylogenetic analysis by the multilocus sequence analysis (MLSA) of the nucleotide sequence of six concatenated genes (16S rRNA, groEL, gyrA, polC, purH and rpoB) resolved the taxonomic position of the Bacillus strain X2 in the Bacillus subtilis subsp. group. The MLSA showed that it is a member of a clade that includes Bacillus subtilis subsp. stercoris. In in silico DNA-DNA hybridization (DDH), the highest matching score was obtained with Bacillus subtilis subsp. stercoris (87%). The in silico DDH of the genome (G + C 43.7 mol %) shared 48.5%, with Bacillus subtilis subsp. inaquosorum. The MLSA phylogenetic tree and the highest degree of DNA hybridization, indicating that it belongs to the Bacillus subtilis subspecies stercoris.