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

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

Non-Structural Proteins (Nsp): A Marker for Detection of Human Coronavirus Families

SARS-CoV-2 was the cause of the global pandemic that caused a total of 14.9 million deaths during the years 2020 and 2021, according to the WHO. The virus presents a mutation rate between 10-5 and 10-3 substitutions per nucleotide site per cell infection (s/n/c). Due to this, studies aimed at knowing the evolution of this virus could help us to foresee (through the future development of new detection strategies and vaccines that prevent the infection of this virus in human hosts) that a pandemic caused by this virus will be generated again. In this research, we performed a functional annotation and identification of changes in Nsp (non-structural proteins) domains in the coronavirus genome. The comparison of the 13 selected coronavirus pangenomes demonstrated a total of 69 protein families and 57 functions associated with the structural domain's differentials between genomes. A marked evolutionary conservation of non-structural proteins was observed. This allowed us to identify and classify highly pathogenic human coronaviruses into alpha, beta, gamma, and delta groups. The designed Nsp cluster provides insight into the trajectory of SARS-CoV-2, demonstrating that it continues to evolve rapidly. An evolutionary marker allows us to discriminate between phylogenetically divergent groups, viral genotypes, and variants between the alpha and betacoronavirus genera. These types of evolutionary studies provide a window of opportunity to use these Nsp as targets of viral therapies.

Comparative genomic analysis of Cohnella hashimotonis sp. nov. isolated from the International Space Station

A single strain from the family Paenibacillaceae was isolated from the wall behind the Waste Hygiene Compartment aboard the International Space Station (ISS) in April 2018, as part of the Microbial Tracking mission series. This strain was identified as a gram-positive, rod-shaped, oxidase-positive, catalase-negative motile bacterium in the genus Cohnella, designated as F6_2S_P_1^T. The 16S sequence of the F6_2S_P_1^T strain places it in a clade with C. rhizosphaerae and C. ginsengisoli, which were originally isolated from plant tissue or rhizosphere environments. The closest 16S and gyrB matches to strain F6_2S_P_1^T are to C. rhizosphaerae with 98.84 and 93.99% sequence similarity, while a core single-copy gene phylogeny from all publicly available Cohnella genomes places it as more closely related to C. ginsengisoli. Average nucleotide identity (ANI) and digital DNA-DNA hybridization (dDDH) values to any described Cohnella species are <89 and <22%, respectively. The major fatty acids for strain F6_2S_P_1^T are anteiso-C_15:0 (51.7%), iso-C_16:0 (23.1%), and iso-C_15:0 (10.5%), and it is able to metabolize a wide range of carbon compounds. Given the results of the ANI and dDDH analyses, this ISS strain is a novel species within the genus Cohnella for which we propose the name Cohnella hashimotonis, with the type strain F6_2S_P_1^T (=NRRL B-65657^T and DSMZ 115098^T). Because no closely related Cohnella genomes were available, this study generated the whole-genome sequences (WGSs) of the type strains for C. rhizosphaerae and C. ginsengisoli. Phylogenetic and pangenomic analysis reveals that F6_2S_P_1^T, C. rhizosphaerae, and C. ginsengisoli, along with two uncharacterized Cohnella strains, possess a shared set of 332 gene clusters which are not shared with any other WGS of Cohnella species, and form a distinct clade branching off from C. nanjingensis. Functional traits were predicted for the genomes of strain F6_2S_P_1^T and other members of this clade.

Emerging Pathogens in Planetary Health and Lessons from Comparative Genome Analyses of Three Clostridia Species

Clostridioides difficile (CD) is a major planetary health burden. A Gram-positive opportunistic pathogen, CD, colonizes the large intestine and is implicated in sepsis, pseudomembranous colitis, and colorectal cancer. C. difficile infection typically following antibiotic exposure results in dysbiosis of the gut microbiome, and is one of the leading causes of diarrhea in the elderly population. While several studies have focused on the toxigenic strains of CD, gut commensals such as Clostridium butyricum (CB) and Clostridium tertium (CT) could harbor toxin/virulence genes, and thus pose a threat to human health. In this study, we sequenced and characterized three isolates, namely, CT (MALS001), CB (MALS002), and CD (MALS003) for their antimicrobial, cytotoxic, antiproliferative, genomic, and proteomic profiles. Although in vitro cytotoxic and antiproliferative potential were observed predominantly in CD MALS003, genome analysis revealed pathogenic potential of CB MALS002 and CT MALS001. Pangenome analysis revealed the presence of several accessory genes typically involved in fitness, virulence, and resistance characteristics in the core genomes of sequenced strains. The presence of an array of virulence and antimicrobial resistance genes in CB MALS002 and CT MALS001 suggests their potential role as emerging pathogens with significant impact on planetary health.

Unravelling the Evolutionary Dynamics of High-Risk Klebsiella pneumoniae ST147 Clones: Insights from Comparative Pangenome Analysis

BACKGROUND

The high prevalence and rapid emergence of antibiotic resistance in high-risk Klebsiella pneumoniae (KP) ST147 clones is a global health concern and warrants molecular surveillance.

METHODS

A pangenome analysis was performed using publicly available ST147 complete genomes. The characteristics and evolutionary relationships among ST147 members were investigated through a Bayesian phylogenetic analysis.

RESULTS

The large number of accessory genes in the pangenome indicates genome plasticity and openness. Seventy-two antibiotic resistance genes were found to be linked with antibiotic inactivation, efflux, and target alteration. The exclusive detection of the bla_OXA-232 gene within the ColKp3 plasmid of KP_SDL79 suggests its acquisition through horizontal gene transfer. The association of seventy-six virulence genes with the acrAB efflux pump, T6SS system and type I secretion system describes its pathogenicity. The presence of Tn6170, a putative Tn7-like transposon in KP_SDL79 with an insertion at the flanking region of the tnsB gene, establishes its transmission ability. The Bayesian phylogenetic analysis estimates ST147's initial divergence in 1951 and the most recent common ancestor for the entire KP population in 1621.

CONCLUSIONS

Present study highlights the genetic diversity and evolutionary dynamics of high-risk clones of K. pneumoniae. Further inter-clonal diversity studies will help us understand its outbreak more precisely and pave the way for therapeutic interventions.

Identification of Subunits for Novel Universal Vaccines against Three Predominant Serogroups and the Emerging O145 among Avian Pathogenic Escherichia coli by Pan-RV Pipeline

Avian pathogenic Escherichia coli, a causative agent of avian colibacillosis, has been causing serious economic losses in the poultry industry. The increase in multidrug-resistant isolates and the complexity of the serotypes of this pathogen, especially the recently reported emergence of a newly predominant serogroup of O145, make the control of this disease difficult. To address this challenge, a high-throughput screening approach, called Pan-RV (Reverse vaccinology based on pangenome analysis), is proposed to search for universal protective antigens against the three traditional serogroups and the newly emerged O145. Using this approach, a total of 61 proteins regarded as probable antigens against the four important serogroups were screened from the core genome of 127 Avian pathogenic Escherichia coli (APEC) genomes, and six were verified by Western blots using antisera. Overall, our research will provide a foundation for the development of an APEC subunit vaccine against avian colibacillosis. Given the exponential growth of whole-genome sequencing (WGS) data, our Pan-RV pipeline will make screening of bacterial vaccine candidates inexpensive, rapid, and efficient. IMPORTANCE With the emergence of drug resistance and the newly predominant serogroup O145, the control of Avian pathogenic Escherichia coli is facing a serious challenge; an efficient immunological method is urgently needed. Here, for the first time, we propose a high-throughput screening approach to search for universal protective antigens against the three traditional serogroups and the newly emerged O145. Importantly, using this approach, a total of 61 proteins regarded as probable antigens against the four important serogroups were screened, and three were shown to be immunoreactive with all antisera (covering the four serogroups), thereby providing a foundation for the development of APEC subunit vaccines against avian colibacillosis. Further, our Pan-RV pipeline will provide immunological control strategies for pathogens with complex and variable genetic backgrounds such as Escherichia coli and will make screening of bacterial vaccine candidates more inexpensive, rapid, and efficient.

Pseudomonas aeruginosa Detection Using Conventional PCR and Quantitative Real-Time PCR Based on Species-Specific Novel Gene Targets Identified by Pangenome Analysis

Mining novel specific molecular targets and establishing efficient identification methods are significant for detecting Pseudomonas aeruginosa, which can enable P. aeruginosa tracing in food and water. Pangenome analysis was used to analyze the whole genomic sequences of 2017 strains (including 1,000 P. aeruginosa strains and 1,017 other common foodborne pathogen strains) downloaded from gene databases to obtain novel species-specific genes, yielding a total of 11 such genes. Four novel target genes, UCBPP-PA14_00095, UCBPP-PA14_03237, UCBPP-PA14_04976, and UCBPP-PA14_03627, were selected for use, which had 100% coverage in the target strain and were not present in nontarget bacteria. PCR primers (PA1, PA2, PA3, and PA4) and qPCR primers (PA12, PA13, PA14, and PA15) were designed based on these target genes to establish detection methods. For the PCR primer set, the minimum detection limit for DNA was 65.4 fg/μl, which was observed for primer set PA2 of the UCBPP-PA14_03237 gene. The detection limit in pure culture without pre-enrichment was 10⁵ colony-forming units (CFU)/ml for primer set PA1, 10³ CFU/ml for primer set PA2, and 10⁴ CFU/ml for primer set PA3 and primer set PA4. Then, qPCR standard curves were established based on the novel species-specific targets. The standard curves showed perfect linear correlations, with R² values of 0.9901 for primer set PA12, 0.9915 for primer set PA13, 0.9924 for primer set PA14, and 0.9935 for primer set PA15. The minimum detection limit of the real-time PCR (qPCR) assay was 10² CFU/ml for pure cultures of P. aeruginosa. Compared with the endpoint PCR and traditional culture methods, the qPCR assay was more sensitive by one or two orders of magnitude. The feasibility of these methods was satisfactory in terms of sensitivity, specificity, and efficiency after evaluating 29 ready-to-eat vegetable samples and was almost consistent with that of the national standard detection method. The developed assays can be applied for rapid screening and detection of pathogenic P. aeruginosa, providing accurate results to inform effective monitoring measures in order to improve microbiological safety.

Rapid Real-Time Polymerase Chain Reaction for Salmonella Serotyping Based on Novel Unique Gene Markers by Pangenome Analysis

An accurate diagnostic method for Salmonella serovars is fundamental to preventing the spread of associated diseases. A diagnostic polymerase chain reaction (PCR)-based method has proven to be an effective tool for detecting pathogenic bacteria. However, the gene markers currently used in real-time PCR to detect Salmonella serovars have low specificity and are developed for only a few serovars. Therefore, in this study, we explored the novel unique gene markers for 60 serovars that share similar antigenic formulas and show high prevalence using pangenome analysis and developed a real-time PCR to detect them. Before exploring gene markers, the 535 Salmonella genomes were evaluated, and some genomes had serovars different from the designated serovar information. Based on these analyses, serovar-specific gene markers were explored. These markers were identified as genes present in all strains of target serovar genomes but absent in strains of other serovar genomes. Serovar-specific primer pairs were designed from the gene markers, and a real-time PCR method that can distinguish between 60 of the most common Salmonella serovars in a single 96-well plate assay was developed. As a result, real-time PCR showed 100% specificity for 199 Salmonella and 29 non-Salmonella strains. Subsequently, the method developed was applied successfully to both strains with identified serovars and an unknown strain, demonstrating that real-time PCR can accurately detect serovars of strains compared with traditional serotyping methods, such as antisera agglutination. Therefore, our method enables rapid and economical Salmonella serotyping compared with the traditional serotyping method.

An object-oriented framework for evolutionary pangenome analysis

Pangenome analysis is fundamental to explore molecular evolution occurring in bacterial populations. Here, we introduce Pagoo, an R framework that enables straightforward handling of pangenome data. The encapsulated nature of Pagoo allows the storage of complex molecular and phenotypic information using an object-oriented approach. This facilitates to go back and forward to the data using a single programming environment and saving any stage of analysis (including the raw data) in a single file, making it sharable and reproducible. Pagoo provides tools to query, subset, compare, visualize, and perform statistical analyses, in concert with other microbial genomics packages available in the R ecosystem. As working examples, we used 1,000 Escherichia coli genomes to show that Pagoo is scalable, and a global dataset of Campylobacter fetus genomes to identify evolutionary patterns and genomic markers of host-adaptation in this pathogen.

Comparative Genomics of Leuconostoc carnosum

Leuconostoc carnosum is a known colonizer of meat-related food matrices. It reaches remarkably high loads during the shelf life in packaged meat products and plays a role in spoilage, although preservative effects have been proposed for some strains. In this study, the draft genomes of 17 strains of L. carnosum (i.e., all the strains that have been sequenced so far) were compared to decipher their metabolic and functional potential and to determine their role in food transformations. Genome comparison and pathway reconstruction indicated that L. carnosum is a compact group of closely related heterofermentative bacteria sharing most of the metabolic features. Adaptation to a nitrogen-rich environment, such as meat, is evidenced by 23 peptidase genes identified in the core genome and by the autotrophy for nitrogen compounds including several amino acids, vitamins, and cofactors. Genes encoding the decarboxylases yielding biogenic amines were not present. All the strains harbored 1–4 of 32 different plasmids, bearing functions associated to proteins hydrolysis, transport of amino acids and oligopeptides, exopolysaccharides, and various resistances (e.g., to environmental stresses, bacteriophages, and heavy metals). Functions associated to bacteriocin synthesis, secretion, and immunity were also found in plasmids. While genes for lactococcin were found in most plasmids, only three harbored the genes for leucocin B, a class IIa antilisterial bacteriocin. Determinants of antibiotic resistances were absent in both plasmids and chromosomes.

Comparative Genomics, Pangenome, and Phylogenomic Analyses of Brenneria spp., and Delineation of Brenneria izadpanahii sp. nov

Brenneria species are bacterial plant pathogens mainly affecting woody plants. Association of all members with devastating disorders (e.g., acute oak decline in Iran and United Kingdom) are due to adaptation and pathogenic behavior in response to host and environmental factors. Some species, including B. goodwinii, B. salicis, and B. nigrifluens, also show endophytic residence. Here we show that all species including novel Brenneria sp. are closely related. Gene-based and genome/pangenome-based phylogeny divide the genus into two distinct lineages, Brenneria clades A and B. The two clades were functionally distinct and were consistent with their common and special potential activities as determined via annotation of functional domains. Pangenome analysis demonstrated that the core pathogenicity factors were highly conserved, an hrp gene cluster encoding a type III secretion system was found in all species except B. corticis. An extensive repertoire of candidate virulence factors was identified. Comparative genomics indicated a repertoire of plant cell wall degrading enzymes, metabolites/antibiotics, and numerous prophages providing new insights into Brenneria-host interactions and appropriate targets for further characterization. This work not only documented the genetic differentiation of Brenneria species but also delineates a more functionally driven understanding of Brenneria by comparison with relevant Pectobacteriaceae thereby substantially enriching the extent of information available for functional genomic investigations.

Not Just a Pathogen? Description of a Plant-Beneficial Pseudomonas syringae Strain

Plants develop in a microbe-rich environment and must interact with a plethora of microorganisms, both pathogenic and beneficial. Indeed, such is the case of Pseudomonas, and its model organisms P. fluorescens and P. syringae, a bacterial genus that has received particular attention because of its beneficial effect on plants and its pathogenic strains. The present study aims to compare plant-beneficial and pathogenic strains belonging to the P. syringae species to get new insights into the distinction between the two types of plant–microbe interactions. In assays carried out under greenhouse conditions, P. syringae pv. syringae strain 260-02 was shown to promote plant-growth and to exert biocontrol of P. syringae pv. tomato strain DC3000, against the Botrytis cinerea fungus and the Cymbidium Ringspot Virus. This P. syringae strain also had a distinct volatile emission profile, as well as a different plant-colonization pattern, visualized by confocal microscopy and gfp labeled strains, compared to strain DC3000. Despite the different behavior, the P. syringae strain 260-02 showed great similarity to pathogenic strains at a genomic level. However, genome analyses highlighted a few differences that form the basis for the following hypotheses regarding strain 260-02. P. syringae strain 260-02: (i) possesses non-functional virulence genes, like the mangotoxin-producing operon Mbo; (ii) has different regulation pathways, suggested by the difference in the autoinducer system and the lack of a virulence activator gene; (iii) has genes encoding DNA methylases different from those found in other P. syringae strains, suggested by the presence of horizontal-gene-transfer-obtained methylases that could affect gene expression.

Comparative Genomics of Rhodococcus equi Virulence Plasmids Indicates Host-Driven Evolution of the vap Pathogenicity Island

The conjugative virulence plasmid is a key component of the Rhodococcus equi accessory genome essential for pathogenesis. Three host-associated virulence plasmid types have been identified: the equine pVAPA and porcine pVAPB circular variants, and the linear pVAPN found in bovine (ruminant) isolates. We recently characterized the R. equi pangenome (Anastasi E, et al. 2016. Pangenome and phylogenomic analysis of the pathogenic actinobacterium Rhodococcus equi. Genome Biol Evol. 8:3140–3148.) and we report here the comparative analysis of the virulence plasmid genomes. Plasmids within each host-associated type were highly similar despite their diverse origins. Variation was accounted for by scattered single nucleotide polymorphisms and short nucleotide indels, while larger indels—mostly in the plasticity region near the vap pathogencity island (PAI)—defined plasmid genomic subtypes. Only one of the plasmids analyzed, of pVAPN type, was exceptionally divergent due to accumulation of indels in the housekeeping backbone. Each host-associated plasmid type carried a unique PAI differing in vap gene complement, suggesting animal host-specific evolution of the vap multigene family. Complete conservation of the vap PAI was observed within each host-associated plasmid type. Both diversity of host-associated plasmid types and clonality of specific chromosomal-plasmid genomic type combinations were observed within the same R. equi phylogenomic subclade. Our data indicate that the overall strong conservation of the R. equi host-associated virulence plasmids is the combined result of host-driven selection, lateral transfer between strains, and geographical spread due to international livestock exchanges.

Pangenome and Phylogenomic Analysis of the Pathogenic Actinobacterium Rhodococcus equi

We report a comparative study of 29 representative genomes of the animal pathogen Rhodococcus equi The analyses showed that R. equi is genetically homogeneous and clonal, with a large core genome accounting for ≈80% of an isolates' gene content. An open pangenome, even distribution of accessory genes among the isolates, and absence of significant core-genome recombination, indicated that gene gain/loss is a main driver of R. equi genome evolution. Traits previously predicted to be important in R. equi physiology, virulence and niche adaptation were part of the core genome. This included the lack of a phosphoenolpyruvate:carbohydrate transport system (PTS), unique among the rhodococci except for the closely related Rhodococcus defluvii, reflecting selective PTS gene loss in the R. equi-R. defluvii sublineage. Thought to be asaccharolytic, rbsCB and glcP non-PTS sugar permease homologues were identified in the core genome and, albeit inefficiently, R. equi utilized their putative substrates, ribose and (irregularly) glucose. There was no correlation between R. equi whole-genome phylogeny and host or geographical source, with evidence of global spread of genomovars. The distribution of host-associated virulence plasmid types was consistent with the exchange of the plasmids (and corresponding host shifts) across the R. equi population, and human infection being zoonotically acquired. Phylogenomic analyses demonstrated that R. equi occupies a central position in the Rhodococcus phylogeny, not supporting the recently proposed transfer of the species to a new genus.