Phylogenomics picks out the par excellence markers for species phylogeny in the genus Staphylococcus

Genus-targeted markers for the taxonomic identification and monitoring of coagulase-positive and coagulase-negative Staphylococcus species

The Staphylococcus genus comprises multiple pathogenic and opportunistic species that represent a risk to public health. Epidemiological studies require accurate taxonomic classification of isolates with enough resolution to distinguish clonal complexes. Unfortunately, 16 S rRNA molecular analysis and phenotypic characterization cannot distinguish all species and do not offer enough resolution to assess intraspecific diversity. Other approaches, such as Multilocus Sequence Tagging, provide higher resolution; however, they have been developed for Staphylococcus aureus and a few other species. Here, we developed a set of genus-targeted primers using five orthologous genes (pta, tuf, tpi, groEs, and sarA) to identify all Staphylococcus species within the genus. The primers were initially evaluated using 20 strains from the Collection of Microorganisms of Interest in Animal Health from AGROSAVIA (CMISA), and their amplified sequences were compared to a set of 33 Staphylococcus species. This allowed the taxonomic identification of the strains even on close species and the establishment of intraspecies diversity. To enhance the scope and cost-effectiveness of the proposed strategy, we customized the primer sets for an Illumina paired-end amplicon protocol, enabling gene multiplexing. We assessed five genes across 177 strains, generating 880 paired-end libraries from the CMISA. This approach significantly reduced sequencing costs, as all libraries can be efficiently sequenced in a single MiSeq run at a fraction (one-fourth or less) of the cost associated with Sanger sequencing. In summary, this method can be used for precise identification and diversity analysis of Staphylococcus species, offering an advancement over traditional techniques in both resolution and cost-effectiveness.

Phylogenomic Analyses of Three Distinct Lineages Uniting Staphylococcus cohnii and Staphylococcus urealyticus from Diverse Hosts

We sequenced and assembled genomes for 17 isolates of Staphylococcus cohnii isolated from osteomyelitis lesions in young broilers from two separate experiments where we induced lameness using a hybrid wire-litter flooring system. Whole genome comparisons using three different methods support a close relationship of genomes from both S. cohnii and Staphylococcus urealyticus. The data support three different lineages, which we designated as Lineage 1, Lineage 2, and Lineage 3, uniting these two species within an evolving complex. We present evidence for horizontal transfer between lineages of genomic regions from 50-440 kbp. The transfer of a 186 kbp region from Lineage 1 to Lineage 2 appears to have generated Lineage 3. Human-associated isolates appear to be limited to Lineages 2 and 3 but Lineage 2 appears to contain a higher number of human pathogenic isolates. The chicken isolates from our lameness trials included genomically diverse isolates from both Lineage 1 and 2, and isolates from both lineages were obtained from osteomyelitis lesions of individual birds. Our results expand the diversity of Staphylococci associated with osteomyelitis in poultry and suggest a high diversity in the microbiome of day-old chicks. Our data also support a reevaluation and unification of the taxonomic classifications of S. cohnii and S. urealyticus.

Discovery of novel treponemes associated with pododermatitis in elk (Cervus canadensis)

Pododermatitis, also known as treponeme-associated hoof disease (TAHD), presents a significant challenge to elk (Cervus canadensis) populations in the northwestern USA, with Treponema spp. consistently implicated in the lesion development. However, identifying species-specific Treponema strains from these lesions is hindered by its culture recalcitrance and limited genomic information. This study utilized shotgun sequencing, in silico genome reconstruction, and comparative genomics as a culture-independent approach to identify metagenome-assembled Treponema genomes (MATGs) from skin scraping samples collected from captive elk experimentally challenged with TAHD. The genomic analysis revealed 10 new MATGs, with 6 representing novel genomospecies associated with pododermatitis in elk and 4 corresponding to previously identified species-Treponema pedis and Treponema phagedenis. Importantly, genomic signatures of novel genomospecies identified in this study were consistently detected in biopsy samples of free-ranging elk diagnosed with TAHD, indicating a potential etiologic association. Comparative metabolic profiling of the MATGs against other Treponema genomes showed a distinct metabolic profile, suggesting potential host adaptation or geographic uniqueness of these newly identified genomospecies. The discovery of novel Treponema genomospecies enhances our understanding of the pathogenesis of pododermatitis and lays the foundation for the development of improved molecular surveillance tools to monitor and manage the disease in free-ranging elk.IMPORTANCETreponema spp. play an important role in the development of pododermatitis in free-ranging elk; however, the species-specific detection of Treponema from pododermatitis lesions is challenging due to culture recalcitrance and limited genomic information. The study utilized shotgun sequencing and in silico genome reconstruction to identify novel Treponema genomospecies from elk with pododermatitis. The discovery of the novel Treponema species opens new avenues to develop molecular diagnostic and epidemiologic tools for the surveillance of pododermatitis in elk. These findings significantly enhance our understanding of the genomic landscape of the Treponemataceae consortium while offering valuable insights into the etiology and pathogenesis of emerging pododermatitis in elk populations.

First report on in-depth genome and comparative genome analysis of a metal-resistant bacterium Acinetobacter pittii S-30, isolated from environmental sample

A newly isolated bacterium Acinetobacter pittii S-30 was recovered from waste-contaminated soil in Ranchi, India. The isolated bacterium belongs to the ESKAPE organisms which represent the major nosocomial pathogens that exhibit high antibiotic resistance. Furthermore, average nucleotide identity (ANI) analysis also showed its closest match (>95%) to other A. pittii genomes. The isolate showed metal-resistant behavior and was able to survive up to 5 mM of ZnSO4. Whole genome sequencing and annotations revealed the occurrence of various genes involved in stress protection, motility, and metabolism of aromatic compounds. Moreover, genome annotation identified the gene clusters involved in secondary metabolite production (biosynthetic gene clusters) such as arylpolyene, acinetobactin like NRP-metallophore, betalactone, and hserlactone-NRPS cluster. The metabolic potential of A. pittii S-30 based on cluster of orthologous, and Kyoto Encyclopedia of Genes and Genomes indicated a high number of genes related to stress protection, metal resistance, and multiple drug-efflux systems etc., which is relatively rare in A. pittii strains. Additionally, the presence of various carbohydrate-active enzymes such as glycoside hydrolases (GHs), glycosyltransferases (GTs), and other genes associated with lignocellulose breakdown suggests that strain S-30 has strong biomass degradation potential. Furthermore, an analysis of genetic diversity and recombination in A. pittii strains was performed to understand the population expansion hypothesis of A. pittii strains. To our knowledge, this is the first report demonstrating the detailed genomic characterization of a heavy metal-resistant bacterium belonging to A. pittii. Therefore, the A. pittii S-30 could be a good candidate for the promotion of plant growth and other biotechnological applications.

In-depth genome and pan-genome analysis of a metal-resistant bacterium Pseudomonas parafulva OS-1

A metal-resistant bacterium Pseudomonas parafulva OS-1 was isolated from waste-contaminated soil in Ranchi City, India. The isolated strain OS-1 showed its growth at 25-45°C, pH 5.0-9.0, and in the presence of ZnSO₄ (upto 5 mM). Phylogenetic analysis based on 16S rRNA gene sequences revealed that strain OS-1 belonged to the genus Pseudomonas and was most closely related to parafulva species. To unravel the genomic features, we sequenced the complete genome of P. parafulva OS-1 using Illumina HiSeq 4,000 sequencing platform. The results of average nucleotide identity (ANI) analysis indicated the closest similarity of OS-1 to P. parafulva PRS09-11288 and P. parafulva DTSP2. The metabolic potential of P. parafulva OS-1 based on Clusters of Othologous Genes (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) indicated a high number of genes related to stress protection, metal resistance, and multiple drug-efflux, etc., which is relatively rare in P. parafulva strains. Compared with other parafulva strains, P. parafulva OS-1 was found to have the unique β-lactam resistance and type VI secretion system (T6SS) gene. Additionally, its genomes encode various CAZymes such as glycoside hydrolases and other genes associated with lignocellulose breakdown, suggesting that strain OS-1 have strong biomass degradation potential. The presence of genomic complexity in the OS-1 genome indicates that horizontal gene transfer (HGT) might happen during evolution. Therefore, genomic and comparative genome analysis of parafulva strains is valuable for further understanding the mechanism of resistance to metal stress and opens a perspective to exploit a newly isolated bacterium for biotechnological applications.

Pan-genome evolution and its association with divergence of metabolic functions in Bifidobacterium genus

Previous studies were mainly focused on genomic evolution and diversity of type species of Bifidobacterium genus due to their health-promoting effect on host. However, those studies were mainly based on species-level taxonomic resolution, adaptation, and characterization of carbohydrate metabolic features of the bifidobacterial species. Here, a comprehensive analysis of the type strain genome unveils the association of pan-genome evolution with the divergence of metabolic function of the Bifidobacterium genus. This study has also demonstrated that horizontal gene transfer, as well as genome expansion and reduction events, leads to the divergence of metabolic functions in Bifidobacterium genus. Furthermore, the genome-based search of probiotic traits among all the available bifidobacterial type strains gives hints on type species, that could confer health benefits to nutrient-deficient individuals. Altogether, the present study provides insight into the developments of genomic evolution, functional divergence, and potential probiotic type species of the Bifidobacterium genus.

A High-Throughput Short Sequence Typing Scheme for Serratia marcescens Pure Culture and Environmental DNA

Molecular typing methods are used to characterize the relatedness between bacterial isolates involved in infections. These approaches rely mostly on discrete loci or whole-genome sequencing (WGS) analyses of pure cultures. On the other hand, their application to environmental DNA profiling to evaluate epidemiological relatedness among patients and environments has received less attention. We developed a specific, high-throughput short sequence typing (HiSST) method for the opportunistic human pathogen Serratia marcescens. Genes displaying the highest polymorphism were retrieved from the core genome of 60 S. marcescens strains. Bioinformatics analyses showed that use of only three loci (within bssA, gabR, and dhaM) distinguished strains with a high level of efficiency. This HiSST scheme was applied to an epidemiological survey of S. marcescens in a neonatal intensive care unit (NICU). In a first case study, a strain responsible for an outbreak in the NICU was found in a sink drain of this unit, by using HiSST scheme and confirmed by WGS. The HiSST scheme was also applied to environmental DNA extracted from sink-environment samples. Diversity of S. marcescens was modest, with 11, 6, and 4 different sequence types (ST) of gabR, bssA, and dhaM loci among 19 sink drains, respectively. Epidemiological relationships among sinks were inferred on the basis of pairwise comparisons of ST profiles. Further research aimed at relating ST distribution patterns to environmental features encompassing sink location, utilization, and microbial diversity is needed to improve the surveillance and management of opportunistic pathogens. IMPORTANCE Serratia marcescens is an important opportunistic human pathogen, often multidrug resistant and involved in outbreaks of nosocomial infections in neonatal intensive care units. Here, we propose a quick and user-friendly method to select the best typing scheme for nosocomial outbreaks in relating environmental and clinical sources. This method, named high-throughput short sequence typing (HiSST), allows to distinguish strains and to explore the diversity profile of nonculturable S. marcescens. The application of HiSST profile analysis for environmental DNA offers new possibilities to track opportunistic pathogens, identify their origin, and relate their distribution pattern with environmental features encompassing sink location, utilization, and microbial diversity. Adaptation of the method to other opportunistic pathogens is expected to improve knowledge regarding their ecology, which is of significant interest for epidemiological risk assessment and elaborate outbreak mitigation strategies.

The mutL Gene as a Genome-Wide Taxonomic Marker for High Resolution Discrimination of Lactiplantibacillus plantarum and Its Closely Related Taxa

The current taxonomy of the Lactiplantibacillus plantarum group comprises of 17 closely related species that are indistinguishable from each other by using commonly used 16S rRNA gene sequencing. In this study, a whole-genome-based analysis was carried out for exploring the highly distinguished target genes whose interspecific sequence identity is significantly less than those of 16S rRNA or conventional housekeeping genes. In silico analyses of 774 core genes by the cano-wgMLST_BacCompare analytics platform indicated that csbB, morA, murI, mutL, ntpJ, rutB, trmK, ydaF, and yhhX genes were the most promising candidates. Subsequently, the mutL gene was selected, and the discrimination power was further evaluated using Sanger sequencing. Among the type strains, mutL exhibited a clearly superior sequence identity (61.6–85.6%; average: 66.6%) to the 16S rRNA gene (96.7–100%; average: 98.4%) and the conventional phylogenetic marker genes (e.g., dnaJ, dnaK, pheS, recA, and rpoA), respectively, which could be used to separat tested strains into various species clusters. Consequently, species-specific primers were developed for fast and accurate identification of L. pentosus, L. argentoratensis, L. plantarum, and L. paraplantarum. During this study, one strain (BCRC 06B0048, L. pentosus) exhibited not only relatively low mutL sequence identities (97.0%) but also a low digital DNA–DNA hybridization value (78.1%) with the type strain DSM 20314^T, signifying that it exhibits potential for reclassification as a novel subspecies. Our data demonstrate that mutL can be a genome-wide target for identifying and classifying the L. plantarum group species and for differentiating novel taxa from known species.

Whole genome sequence and comparative genome analyses of multi-resistant Staphylococcus warneri GD01 isolated from a diseased pig in China

Staphylococcus warneri is a coagulase-negative staphylococcus that is a normal inhabitant of the skin. It is also considered to be an opportunistic etiological agent causing significant infections in human and animals. Currently, relatively little attention has been paid to the genome biology of S. warneri pathogenicity and antibiotic resistance, which are emerging issues for this etiological agent with considerably clinical significance. In this study, we determined the complete genome sequence of S. warneri strain GD01 recovered from the sampled muscle abscess tissue of a diseased pig in South China. The genome of S. warneri is composed of a circular chromosome of 2,473,911 base pairs as well as eight plasmid sequences. Genome-wide metabolic reconstruction revealed 82 intact functional modules driving the catabolism of respiration and fermentation for energy production, uptake of distinct sugars as well as two-component regulatory systems. The evidence uncovered herein enables better understanding for metabolic potential and physiological traits of this etiological agent. The antibiotic susceptibility test demonstrated that S. warneri GD01 was resistant to penicillin, amoxicillin, ampicillin, cefalexin, vancomycin, and sulfisoxazole. The associations between antibiotic phenotypes and the related genotypes were identified to reveal the molecular basis conferring resistance to this pathogen. A number of genes coding for potential virulence factors were firstly depicted in the genome of S. warneri GD01, including adhesins, exoenzymes, capsule, and iron acquisition proteins. Our study provides a valuable genomic context of the genes/modules devoting to metabolism, antibiotic resistance, and virulence of S. warneri.

Phylogenomic Analyses of Members of the Widespread Marine Heterotrophic Genus Pseudovibrio Suggest Distinct Evolutionary Trajectories and a Novel Genus, Polycladidibacter gen. nov

Bacteria belonging to the Pseudovibrio genus are widespread, metabolically versatile, and able to thrive as both free-living and host-associated organisms. Although more than 50 genomes are available, a comprehensive comparative genomics study to resolve taxonomic inconsistencies is currently missing. We analyzed all available genomes and used 552 core genes to perform a robust phylogenomic reconstruction. This in-depth analysis revealed the divergence of two monophyletic basal lineages of strains isolated from polyclad flatworm hosts, namely, Pseudovibrio hongkongensis and Pseudovibrio stylochi These strains have reduced genomes and lack sulfur-related metabolisms and major biosynthetic gene clusters, and their environmental distribution appears to be tightly associated with invertebrate hosts. We showed experimentally that the divergent strains are unable to utilize various sulfur compounds that, in contrast, can be utilized by the type strain Pseudovibrio denitrificans Our analyses suggest that the lineage leading to these two strains has been subject to relaxed purifying selection resulting in great gene loss. Overall genome relatedness indices (OGRI) indicate substantial differences between the divergent strains and the rest of the genus. While 16S rRNA gene analyses do not support the establishment of a different genus for the divergent strains, their substantial genomic, phylogenomic, and physiological differences strongly suggest a divergent evolutionary trajectory and the need for their reclassification. Therefore, we propose the novel genus Polycladidibacter gen. nov.IMPORTANCE The genus Pseudovibrio is commonly associated with marine invertebrates, which are essential for ocean health and marine nutrient cycling. Traditionally, the phylogeny of the genus has been based on 16S rRNA gene analysis. The use of the 16S rRNA gene or any other single marker gene for robust phylogenetic placement has recently been questioned. We used a large set of marker genes from all available Pseudovibrio genomes for in-depth phylogenomic analyses. We identified divergent monophyletic basal lineages within the Pseudovibrio genus, including two strains isolated from polyclad flatworms. These strains showed reduced sulfur metabolism and biosynthesis capacities. The phylogenomic analyses revealed distinct evolutionary trajectories and ecological adaptations that differentiate the divergent strains from the other Pseudovibrio members and suggest that they fall into a novel genus. Our data show the importance of widening the use of phylogenomics for better understanding bacterial physiology, phylogeny, and evolution.

Phylogenomics Reveals Clear Cases of Misclassification and Genus-Wide Phylogenetic Markers for Acinetobacter

The Gram-negative Acinetobacter genus has several species of clear medical relevance. Many fully sequenced genomes belonging to the genus have been published in recent years; however, there has not been a recent attempt to infer the evolutionary history of Acinetobacter with that vast amount of information. Here, through a phylogenomic approach, we established the most up-to-date view of the evolutionary relationships within this genus and highlighted several cases of poor classification, especially for the very closely related species within the Acinetobacter calcoaceticus-Acinetobacter baumannii complex (Acb complex). Furthermore, we determined appropriate phylogenetic markers for this genus and showed that concatenation of the top 13 gives a very decent reflection of the evolutionary relationships for the genus Acinetobacter. The intersection between our top markers and previously defined universal markers is very small. In general, our study shows that, although there seems to be hardly any universal markers, bespoke phylogenomic approaches can be used to infer the phylogeny of different bacterial genera. We expect that ad hoc phylogenomic approaches will be the standard in the years to come and will provide enough information to resolve intricate evolutionary relationships like those observed in the Acb complex.