First Steps in the Analysis of Prokaryotic Pan-Genomes

Ecogenomic insights into the resilience of keystone Blastococcus Species in extreme environments: a comprehensive analysis

BACKGROUND

The stone-dwelling genus Blastococcus plays a key role in ecosystems facing extreme conditions such as drought, salinity, alkalinity, and heavy metal contamination. Despite its ecological significance, little is known about the genomic factors underpinning its adaptability and resilience in such harsh environments. This study investigates the genomic basis of Blastococcus's adaptability within its specific microniches, offering insights into its potential for biotechnological applications.

RESULTS

Comprehensive pangenome analysis revealed that Blastococcus possesses a highly dynamic genetic composition, characterized by a small core genome and a large accessory genome, indicating significant genomic plasticity. Ecogenomic assessments highlighted the genus's capabilities in substrate degradation, nutrient transport, and stress tolerance, particularly on stone surfaces and archaeological sites. The strains also exhibited plant growth-promoting traits, enhanced heavy metal resistance, and the ability to degrade environmental pollutants, positioning Blastococcus as a candidate for sustainable agriculture and bioremediation. Interestingly, no correlation was found between the ecological or plant growth-promoting traits (PGPR) of the strains and their isolation source, suggesting that these traits are not linked to their specific environments.

CONCLUSIONS

This research highlights the ecological and biotechnological potential of Blastococcus species in ecosystem health, soil fertility improvement, and stress mitigation strategies. It calls for further studies on the adaptation mechanisms of the genus, emphasizing the need to validate these findings through wet lab experiments. This study enhances our understanding of microbial ecology in extreme environments and supports the use of Blastococcus in environmental management, particularly in soil remediation and sustainable agricultural practices.

From the depths of the Java Trench: genomic analysis of Priestia flexa JT4 reveals bioprospecting and lycopene production potential

BACKGROUND

The marine environment boasts distinctive physical, chemical, and biological characteristics. While numerous studies have delved into the microbial ecology and biological potential of the marine environment, exploration of genetically encoded, deep-sea sourced secondary metabolites remains scarce. This study endeavors to investigate marine bioproducts derived from deep-sea water samples at a depth of 1,000 m in the Java Trench, Indonesia, utilizing both culture-dependent and whole-genome sequencing methods.

RESULTS

Our efforts led to the successful isolation and cultivation of a bacterium Priestia flexa JT4 from the water samples, followed by comprehensive genome sequencing. The resultant high-quality draft genome, approximately 4 Mb, harbored 5185 coding sequences (CDSs). Notably, 61.97% of these CDSs were inadequately characterized, presenting potential novel CDSs. This study is the first to identify the "open-type" (α < 1) pangenome within the genus Priestia. Moreover, our analysis uncovered eight biosynthetic gene clusters (BGCs) using the common genome mining pipeline, antiSMASH. Two non-ribosomal peptide synthetase (NRPS) BGCs within these clusters exhibited the potential to generate novel biological compounds. Noteworthy is the confirmation that the terpene BGC in P. flexa JT4 can produce lycopene, a compound in substantial industrial demand. The presence of lycopene in the P. flexa JT4 cells was verified using Ultra-performance liquid chromatography-mass spectrometry (UPLC-MS/MS) in multiple reaction modes.

CONCLUSIONS

This study highlights the bioprospecting opportunity to explore novel bioproducts and lycopene compounds from P. flexa JT4. It marks the pioneering exploration of deep-sea bacterium bioprospecting in Indonesia, seeking to unveil novel bioproducts and lycopene compounds through a genome mining approach.

Comparative genomics of Fervidobacterium: a new phylogenomic landscape of these wide-spread thermophilic anaerobes

BACKGROUND

Fervidobacterium is a genus of thermophilic anaerobic Gram-negative rod-shaped bacteria belonging to the phylum Thermotogota. They can grow through fermentation on a wide range of sugars and protein-rich substrates. Some can also break down feather keratin, which has significant biotechnological potential. Fervidobacteria genomes have undergone several horizontal gene transfer events, sharing DNA with unrelated microbial taxa. Despite increasing biotechnological and evolutionary interest in this genus, only seven species have been described to date. Here, we present and describe six new and complete Fervidobacterium genomes, including the type strains Fervidobacterium gondwanense CBS-1 T, F. islandicum H-21 T and F. thailandense FC2004T, one novel isolate from Georgia (strain GSH) and two strains (DSM 21710 and DSM 13770) that have not been previously described along with an evolutionary and phylogenomic analysis of the genus.

RESULTS

The complete genomes were around 2 Mb with approximately 2,000 CDS identified and annotated in each of them and a G + C content ranging from 38.9 mol% to 45.8 mol%. Phylogenomic comparisons of all currently available Fervidobacterium genomes, including OrthoANI and TYGS analyses, as well as a phylogenetic analysis based on the 16S rRNA gene, identified six species and nine subspecies clusters across the genus, with a consistent topology and a distant and separately branching species, Fervidobacterium thailandense. F. thailandense harbored the highest number of transposases, CRISPR clusters, pseudo genes and horizontally transferred regions The pan genome of the genus showed that 44% of the genes belong to the cloud pangenome, with most of the singletons found also in F. thailandense.

CONCLUSIONS

The additional genome sequences described in this work and the comparison with all available Fervidobacterium genome sequences provided new insights into the evolutionary history of this genus and supported a phylogenetic reclassification. The phylogenomic results from OrthoANI and TYGS analyses revealed that F. riparium and F. gondwanense belong to the same genome species, and includes Fervidobacterium sp. 13770, while "F. pennivorans" strain DYC belongs to a separate genome species, whereas Fervidobacterium sp. 21710 and Fervidobacterium sp. GSH within the Fervidobacterium pennivorans clade represent two subspecies. F. changbaicum is reclassified as F. islandicum.

Unitig-centered pan-genome machine learning approach for predicting antibiotic resistance and discovering novel resistance genes in bacterial strains

In current genomic research, the widely used methods for predicting antimicrobial resistance (AMR) often rely on prior knowledge of known AMR genes or reference genomes. However, these methods have limitations, potentially resulting in imprecise predictions owing to incomplete coverage of AMR mechanisms and genetic variations. To overcome these limitations, we propose a pan-genome-based machine learning approach to advance our understanding of AMR gene repertoires and uncover possible feature sets for precise AMR classification. By building compacted de Brujin graphs (cDBGs) from thousands of genomes and collecting the presence/absence patterns of unique sequences (unitigs) for Pseudomonas aeruginosa, we determined that using machine learning models on unitig-centered pan-genomes showed significant promise for accurately predicting the antibiotic resistance or susceptibility of microbial strains. Applying a feature-selection-based machine learning algorithm led to satisfactory predictive performance for the training dataset (with an area under the receiver operating characteristic curve (AUC) of > 0.929) and an independent validation dataset (AUC, approximately 0.77). Furthermore, the selected unitigs revealed previously unidentified resistance genes, allowing for the expansion of the resistance gene repertoire to those that have not previously been described in the literature on antibiotic resistance. These results demonstrate that our proposed unitig-based pan-genome feature set was effective in constructing machine learning predictors that could accurately identify AMR pathogens. Gene sets extracted using this approach may offer valuable insights into expanding known AMR genes and forming new hypotheses to uncover the underlying mechanisms of bacterial AMR.

Unlocking the Mycobacteroides abscessus pan-genome using computational tools: insights into evolutionary dynamics and lifestyle

Mycobacteroides abscessus is a non-tuberculous mycobacteria implicated in causing lung infections. It is difficult to control owing to resistance to antibiotics and disinfectants. This work was aimed at comprehending: the pan-genome architecture, evolutionary dynamics, and functionalities of pan-genome components linked to COGs and KEGG. Around 2802 core genes were present in each strain of the M. abscessus genome. The number of accessory genes ranged from 1615 to 2481. The open pan-genome of M. abscessus was attributed to the accessory genes underlining its adaptability in the host. Phylogenetic analysis revealed cluster-based relationships and highlighted factors shaping variability and adaptive capabilities. Transcription, metabolism, and pathogenic genes were vital for M. abscessus lifestyle. The accessory genes contributed to the diverse metabolic capability. The incidence of a significant portion of secondary metabolite biosynthesis genes provided insights for investigating their biosynthetic gene clusters. Additionally, a high proportion of xenobiotic biodegradation genes highlighted potential metabolic capabilities. In silico screening identified a potential vaccine candidate among hypothetical proteins in COGs. Functional analysis of M. abscessus pan-genome components unveiled factors associated with virulence, pathogenicity, infection establishment, persistence, and resistance. Notable amongst them were: MMPL family transporters, PE-PPE domain-containing proteins, TetR family transcriptional regulators, ABC transporters, Type-I, II, III, VII secretion proteins, DUF domain-containing proteins, cytochrome P450, VapC family toxin, virulence factor Mce family protein, type II toxin-antitoxin system. Overall, these results enhanced understanding of the metabolism, host-pathogen dynamics, pathogenic lifestyle, and adaptations. This will facilitate further investigations for combating infections and designing suitable therapies.

Pangenome analysis of Paenibacillus polymyxa strains reveals the existence of multiple and functionally distinct Paenibacillus species

Paenibacillus polymyxa, a Gram-positive bacterium commonly found in soil and plant roots, plays an important role in the environment due to its nitrogen-fixing ability and is renowned for producing antibiotics like polymyxin. In this study, we present a robust framework for investigating the evolutionary and taxonomic connections of strains belonging to P. polymyxa available at the National Center for Biotechnology Information, as well as five new additional strains isolated at the University of Camerino (Italy), through pangenome analysis. These strains can produce secondary metabolites active against Staphylococcus aureus and Klebsiella pneumoniae. Employing techniques such as digital DNA-DNA hybridization (dDDH), average nucleotide identity (ANI) estimation, OrthoFinder, and ribosomal multilocus sequence typing, we consistently divided these P. polymyxa strains into four clusters, which differ significantly in terms of ANI and dDDH percentages, both considered as reference indices for separating bacterial species. Moreover, the strains of Cluster 2 were re-classified as belonging to the Paenibacillus ottowii species. By comparing the pangenomes, we identified the core genes of each cluster and analyzed them to recognize distinctive features in terms of biosynthetic/metabolic potential. The comparison of pangenomes also allowed us to pinpoint differences between clusters in terms of genetic variability and the percentage of the genome dedicated to core and accessory genes. In conclusion, the data obtained from our analyses of strains belonging to the P. polymyxa species converge toward a necessary reclassification, which will require a fundamental contribution from microbiologists in the near future.

IMPORTANCE

The development of sequencing technologies has led to an exponential increase in microbial sequencing data. Accurately identifying bacterial species remains a challenge because of extensive intra-species variability, the need for multiple identification methods, and the rapid rate of taxonomic changes. A substantial contribution to elucidating the relationships among related bacterial strains comes from comparing their genomic sequences. This comparison also allows for the identification of the "pangenome," which is the set of genes shared by all individuals of a species, as well as the set of genes that are unique to subpopulations. Here, we applied this approach to Paenibacillus polymyxa, a species studied for its potential as a biofertilizer and biocontrol agent and known as an antibiotic producer. Our work highlights the need for a more efficient classification of this bacterial species and provides a better delineation of strains with different properties.

Chromosome architecture as a determinant for biosynthetic diversity in Micromonospora

Natural products - small molecules generated by organisms to facilitate ecological interactions - are of great importance to society and are used as antibacterial, antiviral, antifungal and anticancer drugs. However, the role and evolution of these molecules and the fitness benefits they provide to their hosts in their natural habitat remain an outstanding question. In bacteria, the genes that encode the biosynthetic proteins that generate these molecules are organised into discrete loci termed biosynthetic gene clusters (BGCs). In this work, we asked the following question: How are biosynthetic gene clusters organised at the chromosomal level? We sought to answer this using publicly available high-quality assemblies of Micromonospora, an actinomycete genus with members responsible for biosynthesizing notable natural products, such as gentamicin and calicheamicin. By orienting the Micromonospora chromosome around the origin of replication, we demonstrated that Micromonospora has a conserved origin-proximal region, which becomes progressively more disordered towards the antipodes of the origin. We then demonstrated through genome mining of these organisms that the conserved origin-proximal region and the origin-distal region of Micromonospora have distinct populations of BGCs and, in this regard, parallel the organization of Streptomyces, which possesses linear chromosomes. Specifically, the origin-proximal region contains highly syntenous, conserved BGCs predicted to biosynthesize terpenes and a type III polyketide synthase. In contrast, the ori-distal region contains a highly diverse population of BGCs, with many BGCs belonging to unique gene cluster families. These data highlight that genomic plasticity in Micromonospora is locus-specific, and highlight the importance of using high-quality genome assemblies for natural product discovery and guide future natural product discovery by highlighting that biosynthetic novelty may be enriched in specific chromosomal neighbourhoods.

Pangenomes of human gut microbiota uncover links between genetic diversity and stress response

The genetic diversity of the gut microbiota has a central role in host health. Here, we created pangenomes for 728 human gut prokaryotic species, quadrupling the genes of strain-specific genomes. Each of these species has a core set of a thousand genes, differing even between closely related species, and an accessory set of genes unique to the different strains. Functional analysis shows high strain variability associates with sporulation, whereas low variability is linked with antibiotic resistance. We further map the antibiotic resistome across the human gut population and find 237 cases of extreme resistance even to last-resort antibiotics, with a predominance among Enterobacteriaceae. Lastly, the presence of specific genes in the microbiota relates to host age and sex. Our study underscores the genetic complexity of the human gut microbiota, emphasizing its significant implications for host health. The pangenomes and antibiotic resistance map constitute a valuable resource for further research.

A gentle introduction to pangenomics

Pangenomes have emerged in response to limitations associated with traditional linear reference genomes. In contrast to a traditional reference that is (usually) assembled from a single individual, pangenomes aim to represent all of the genomic variation found in a group of organisms. The term 'pangenome' is currently used to describe multiple different types of genomic information, and limited language is available to differentiate between them. This is frustrating for researchers working in the field and confusing for researchers new to the field. Here, we provide an introduction to pangenomics relevant to both prokaryotic and eukaryotic organisms and propose a formalization of the language used to describe pangenomes (see the Glossary) to improve the specificity of discussion in the field.

Improving Beneficial Traits in Bacillus cabrialesii subsp. cabrialesii TE3T through UV-Induced Genomic Changes

It is essential to hunt for new technologies that promote sustainable practices for agroecosystems; thus, the bioprospecting of beneficial microorganisms complementing with mutation induction techniques to improve their genomic, metabolic, and functional traits is a promising strategy for the development of sustainable microbial inoculants. Bacillus cabrialesii subsp. cabrialesii strain TE3T, a previously recognized plant growth-promoting and biological control agent, was subjected to UV mutation induction to improve these agro-biotechnological traits. Dilutions were made which were spread on Petri dishes and placed under a 20 W UV lamp at 10-min intervals for 60 min. After the UV-induced mutation of this strain, 27 bacterial colonies showed morphological differences compared to the wild-type strain; however, only a strain named TE3T-UV25 showed an improvement in 53.6% of the biocontrol against Bipolaris sorokiniana vs. the wild-type strain, by competition of nutrient and space (only detected in the mutant strain), as well as diffusible metabolites. Furthermore, the ability to promote wheat growth was evaluated by carrying out experiments under specific greenhouse conditions, considering un-inoculated, strain TE3T, and strain TE3T-UV25 treatments. Thus, after 120 days, biometric traits in seedlings were quantified and statistical analyses were performed, which showed that strain TE3T-UV25 maintained its ability to promote wheat growth in comparison with the wild-type strain. On the other hand, using bioinformatics tools such as ANI, GGDC, and TYGS, the Overall Genome Relatedness Index (OGRI) and phylogenomic relationship of mutant strain TE3T-UV25 were performed, confirming that it changed its taxonomic affiliation from B. cabrialesii subsp. cabrialesii to Bacillus subtilis. In addition, genome analysis showed that the mutant, wild-type, and B. subtilis strains shared 3654 orthologous genes; however, a higher number of shared genes (3954) was found between the TE3T-UV25 mutant strain and B. subtilis 168, while the mutant strain shared 3703 genes with the wild-type strain. Genome mining was carried out using the AntiSMASH v7.0 web server and showed that mutant and wild-type strains shared six biosynthetic gene clusters associated with biocontrol but additionally, pulcherriminic acid cluster only was detected in the genome of the mutant strain and Rhizocticin A was exclusively detected in the genome of the wild-type strain. Finally, using the PlaBase tool, differences in the number of genes (17) associated with beneficial functions in agroecosystems were detected in the genome of the mutant vs. wild-type strain, such as biofertilization, bioremediation, colonizing plant system, competitive exclusion, phytohormone, plant immune response stimulation, putative functions, stress control, and biocontrol. Thus, the UV-induced mutation was a successful strategy to improve the bioactivity of B. cabrialesii subsp. cabrialesii TE3T related to the agro-biotecnology applications. The obtained mutant strain, B. subtilis TE3T-UV25, is a promising strain to be further studied as an active ingredient for the bioformulation of bacterial inoculants to migrate sustainable agriculture.

Seven quick tips for gene-focused computational pangenomic analysis

Pangenomics is a relatively new scientific field which investigates the union of all the genomes of a clade. The word pan means everything in ancient Greek; the term pangenomics originally regarded genomes of bacteria and was later intended to refer to human genomes as well. Modern bioinformatics offers several tools to analyze pangenomics data, paving the way to an emerging field that we can call computational pangenomics. Current computational power available for the bioinformatics community has made computational pangenomic analyses easy to perform, but this higher accessibility to pangenomics analysis also increases the chances to make mistakes and to produce misleading or inflated results, especially by beginners. To handle this problem, we present here a few quick tips for efficient and correct computational pangenomic analyses with a focus on bacterial pangenomics, by describing common mistakes to avoid and experienced best practices to follow in this field. We believe our recommendations can help the readers perform more robust and sound pangenomic analyses and to generate more reliable results.

Genomic Functional Analysis of Novel Radiation-Resistant Species of Knollia sp. nov. S7-12T from the North Slope of Mount Everest

Radiation protection is an important field of study, as it relates to human health and environmental safety. Radiation-resistance mechanisms in extremophiles are a research hotspot, as this knowledge has great application value in bioremediation and development of anti-radiation drugs. Mount Everest, an extreme environment of high radiation exposure, harbors many bacterial strains resistant to radiation. However, owing to the difficulties in studying them because of the extreme terrain, many remain unexplored. In this study, a novel species (herein, S7-12T) was isolated from the moraine of Mount Everest, and its morphology and functional and genomic characteristics were analyzed. The strain S7-12T is white in color, smooth and rounded, non-spore-forming, and non-motile and can survive at a UV intensity of 1000 J/m2, showing that it is twice as resistant to radiation as Deinococcus radiodurans. Radiation-resistance genes, including IbpA and those from the rec and CspA gene families, were identified. The polyphasic taxonomic approach revealed that the strain S7-12T (=KCTC 59114T =GDMCC 1.3458T) is a new species of the genus Knoellia and is thus proposed to be named glaciei. The in-depth study of the genome of strain S7-12T will enable us to gain further insights into its potential use in radiation resistance. Understanding how microorganisms resist radiation damage could reveal potential biomarkers and therapeutic targets, leading to the discovery of potent anti-radiation compounds, thereby improving human resistance to the threat of radiation.

In vitro and genomic mining studies of anti-Clostridium perfringens Compounds Derived from Bacillus amyloliquefaciens

Clostridium perfringens is an important opportunistic microorganism in commercial poultry production that is implicated in necrotic enteritis (NE) outbreaks. This disease poses a severe financial burden on the global poultry industry, causing estimated annual losses of $6 billion globally. The ban on in-feed antibiotic growth promoters has spurred investigations into approaches of alternatives to antibiotics, among which Bacillus probiotics have demonstrated varying degrees of effectiveness against NE. However, the precise mechanisms underlying Bacillus-mediated beneficial effects on host responses in NE remain to be further elucidated. In this manuscript, we conducted in vitro and genomic mining analysis to investigate anti-C. perfringens activity observed in the supernatants derived from 2 Bacillus amyloliquefaciens strains (FS1092 and BaD747). Both strains demonstrated potent anti-C. perfringens activities in in vitro studies. An analysis of genomes from 15 B. amyloliquefaciens, 11 B. velezensis, and 2 B. subtilis strains has revealed an intriguing clustering pattern among strains known to possess anti-C. perfringens activities. Furthermore, our investigation has identified 7 potential antimicrobial compounds, predicted as secondary metabolites through antiSMASH genomic mining within the published genomes of B. amyloliquefaciens species. Based on in vitro analysis, BaD747 may have the potential as a probiotic in the control of NE. These findings not only enhance our understanding of B. amyloliquefaciens's action against C. perfringens but also provide a scientific rationale for the development of novel antimicrobial therapeutic agents against NE.

MAGs-centric crack: how long will, spore-positive Frankia and most Protofrankia, microsymbionts remain recalcitrant to axenic growth?

Nearly 50 years after the ground-breaking isolation of the primary Comptonia peregrina microsymbiont under axenic conditions, efforts to isolate a substantial number of Protofrankia and Frankia strains continue with enduring challenges and complexities. This study aimed to streamline genomic insights through comparative and predictive tools to extract traits crucial for isolating specific Frankia in axenic conditions. Pangenome analysis unveiled significant genetic diversity, suggesting untapped potential for cultivation strategies. Shared metabolic strategies in cellular components, central metabolic pathways, and resource acquisition traits offered promising avenues for cultivation. Ecological trait extraction indicated that most uncultured strains exhibit no apparent barriers to axenic growth. Despite ongoing challenges, potential caveats, and errors that could bias predictive analyses, this study provides a nuanced perspective. It highlights potential breakthroughs and guides refined cultivation strategies for these yet-uncultured strains. We advocate for tailored media formulations enriched with simple carbon sources in aerobic environments, with atmospheric nitrogen optionally sufficient to minimize contamination risks. Temperature adjustments should align with strain preferences-28-29°C for Frankia and 32-35°C for Protofrankia-while maintaining an alkaline pH. Given potential extended incubation periods (predicted doubling times ranging from 3.26 to 9.60 days, possibly up to 21.98 days), patience and rigorous contamination monitoring are crucial for optimizing cultivation conditions.

Morphological, Metabolomic and Genomic Evidences on Drought Stress Protective Functioning of the Endophyte Bacillus safensis Ni7

The metabolomic and genomic characterization of an endophytic Bacillus safensis Ni7 was carried out in this study. This strain has previously been isolated from the xerophytic plant Nerium indicum L. and reported to enhance the drought tolerance in Capsicum annuum L. seedlings. The effects of drought stress on the morphology, biofilm production, and metabolite production of B. safensis Ni7 are analyzed in the current study. From the results obtained, the organism was found to have multiple strategies such as aggregation and clumping, robust biofilm production, and increased production of surfactin homologues under the drought induced condition when compared to non-stressed condition. Further the whole genome sequencing (WGS) based analysis has demonstrated B. safensis Ni7 to have a genome size of 3,671,999 bp, N50 value of 3,527,239, and a mean G+C content of 41.58%. Interestingly the organism was observed to have the presence of various stress-responsive genes (13, 20U, 16U,160, 39, 17M, 18, 26, and ctc) and genes responsible for surfactin production (srfAA, srfAB, srfAC, and srfAD), biofilm production (epsD, epsE, epsF, epsG, epsH, epsI, epsK, epsL, epsM, epsN, and pel), chemotaxis (cheB_1, cheB_2, cheB_3, cheW_1, cheW_2 cheR, cheD, cheC, cheA, cheY, cheV, and cheB_4), flagella synthesis (flgG_1, flgG_2, flgG_3, flgC, and flgB) as supportive to the drought tolerance. Besides these, the genes responsible for plant growth promotion (PGP), including the genes for nitrogen (nasA, nasB, nasC, nasD, and nasE) and sulfur assimilation (cysL_1&L_2, cysI) and genes for phosphate solubilization (phoA, phoP_1& phoP_2, and phoR) could also be predicted. Along with the same, the genes for catalase, superoxide dismutase, protein homeostasis, cellular fitness, osmoprotectants production, and protein folding could also be predicted from its WGS data. Further pan-genome analysis with plant associated B. safensis strains available in the public databases revealed B. safensis Ni7 to have the presence of a total of 5391 gene clusters. Among these, 3207 genes were identified as core genes, 954 as shell genes and 1230 as cloud genes. This variation in gene content could be taken as an indication of evolution of strains of Bacillus safensis as per specific conditions and hence in the case of B. safensis Ni7 its role in habitat adaptation of plant is well expected. This diversity in endophytic bacterial genes may attribute its role to support the plant system to cope up with stress conditions. Overall, the study provides genomic evidence on Bacillus safensis Ni7 as a stress alleviating microbial partner in plants.

Genome plasticity shapes the ecology and evolution of Phocaeicola dorei and Phocaeicola vulgatus

Phocaeicola dorei and Phocaeicola vulgatus are very common and abundant members of the human gut microbiome and play an important role in the infant gut microbiome. These species are closely related and often confused for one another; yet, their genome comparison, interspecific diversity, and evolutionary relationships have not been studied in detail so far. Here, we perform phylogenetic analysis and comparative genomic analyses of these two Phocaeicola species. We report that P. dorei has a larger genome yet a smaller pan-genome than P. vulgatus. We found that this is likely because P. vulgatus is more plastic than P. dorei, with a larger repertoire of genetic mobile elements and fewer anti-phage defense systems. We also found that P. dorei directly descends from a clade of P. vulgatus¸ and experienced genome expansion through genetic drift and horizontal gene transfer. Overall, P. dorei and P. vulgatus have very different functional and carbohydrate utilisation profiles, hinting at different ecological strategies, yet they present similar antimicrobial resistance profiles.

Ovine and Caprine Strains of Corynebacterium pseudotuberculosis on Czech Farms-A Comparative Study

Caseous lymphadenitis (CLA) is a worldwide disease of small ruminants caused by Corynebacterium pseudotuberculosis, a facultative intracellular pathogen that is able to survive and multiply in certain white blood cells of the host. In this study, 33 strains of C. pseudotuberculosis were isolated from sheep and goats suffering from CLA on nine farms in the Czech Republic. All these strains were tested for their antibiotic susceptibility, ability to form a biofilm and resistance to the effects of commonly used disinfectant agents. To better understand the virulence of C. pseudotuberculosis, the genomes of strains were sequenced and comparative genomic analysis was performed with another 123 genomes of the same species, including ovis and equi biovars, downloaded from the NCBI. The genetic determinants for the virulence factors responsible for adherence and virulence factors specialized for iron uptake and exotoxin phospholipase D were revealed in every analyzed genome. Carbohydrate-Active Enzymes were compared, revealing the presence of genetic determinants encoding exo-α-sialidase (GH33) and the CP40 protein in most of the analyzed genomes. Thirty-three Czech strains of C. pseudotuberculosis were identified as the biovar ovis on the basis of comparative genome analysis. All the compared genomes of the biovar ovis strains were highly similar regardless of their country of origin or host, reflecting their clonal behavior.

First report of coexistence of blaKPC-2 and blaNDM-1 in carbapenem-resistant clinical isolates of Klebsiella aerogenes in Brazil

Klebsiella aerogenes is an important opportunistic pathogen with the potential to develop resistance against last-line antibiotics, such as carbapenems, limiting the treatment options. Here, we investigated the antibiotic resistance profiles of 10 K. aerogenes strains isolated from patient samples in the intensive-care unit of a Brazilian tertiary hospital using conventional PCR and a comprehensive genomic characterization of a specific K. aerogenes strain (CRK317) carrying both the blaKPC-2 and blaNDM-1 genes simultaneously. All isolates were completely resistant to β-lactam antibiotics, including ertapenem, imipenem, and meropenem with differencing levels of resistance to aminoglycosides, quinolones, and tigecycline also observed. Half of the strains studied were classified as multidrug-resistant. The carbapenemase-producing isolates carried many genes of interest including: β-lactams (blaNDM-1, blaKPC-2, blaTEM-1, blaCTX-M-1 group, blaOXA-1 group and blaSHVvariants in 20-80% of the strains), aminoglycoside resistance genes [aac(6')-Ib and aph(3')-VI, 70 and 80%], a fluoroquinolone resistance gene (qnrS, 80%), a sulfonamide resistance gene (sul-2, 80%) and a multidrug efflux system transporter (mdtK, 70%) while all strains carried the efflux pumps Acr (subunit A) and tolC. Moreover, we performed a comprehensive genomic characterization of a specific K. aerogenes strain (CRK317) carrying both the blaKPC-2 and blaNDM-1 genes simultaneously. The draft genome assembly of the CRK317 had a total length of 5,462,831 bp and a GC content of 54.8%. The chromosome was found to contain many essential genes. In silico analysis identified many genes associated with resistance phenotypes, including β-lactamases (blaOXA-9, blaTEM-1, blaNDM-1, blaCTX-M-15, blaAmpC-1, blaAmpC-2), the bleomycin resistance gene (bleMBL), an erythromycin resistance methylase (ermC), aminoglycoside-modifying enzymes [aac(6')-Ib, aadA/ant(3")-Ia, aph(3')-VI], a sulfonamide resistance enzyme (sul-2), a chloramphenicol acetyltransferase (catA-like), a plasmid-mediated quinolone resistance protein (qnrS1), a glutathione transferase (fosA), PEtN transferases (eptA, eptB) and a glycosyltransferase (arnT). We also detected 22 genomic islands, eight families of insertion sequences, two putative integrative and conjugative elements with a type IV secretion system, and eight prophage regions. This suggests the significant involvement of these genetic structures in the dissemination of antibiotic resistance. The results of our study show that the emergence of carbapenemase-producing K. aerogenes, co-harboring blaKPC-2 and blaNDM-1, is a worrying phenomenon which highlights the importance of developing strategies to detect, prevent, and control the spread of these microorganisms.

In Silico Safety Assessment of Bacillus Isolated from Polish Bee Pollen and Bee Bread as Novel Probiotic Candidates

Bacillus species isolated from Polish bee pollen (BP) and bee bread (BB) were characterized for in silico probiotic and safety attributes. A probiogenomics approach was used, and in-depth genomic analysis was performed using a wide array of bioinformatics tools to investigate the presence of virulence and antibiotic resistance properties, mobile genetic elements, and secondary metabolites. Functional annotation and Carbohydrate-Active enZYmes (CAZYme) profiling revealed the presence of genes and a repertoire of probiotics properties promoting enzymes. The isolates BB10.1, BP20.15 (isolated from bee bread), and PY2.3 (isolated from bee pollen) genome mining revealed the presence of several genes encoding acid, heat, cold, and other stress tolerance mechanisms, adhesion proteins required to survive and colonize harsh gastrointestinal environments, enzymes involved in the metabolism of dietary molecules, antioxidant activity, and genes associated with the synthesis of vitamins. In addition, genes responsible for the production of biogenic amines (BAs) and D-/L-lactate, hemolytic activity, and other toxic compounds were also analyzed. Pan-genome analyses were performed with 180 Bacillus subtilis and 204 Bacillus velezensis genomes to mine for any novel genes present in the genomes of our isolates. Moreover, all three isolates also consisted of gene clusters encoding secondary metabolites.

The zoonotic pathogen Wohlfahrtiimonas chitiniclastica - current findings from a clinical and genomic perspective

The zoonotic pathogen Wohlfahrtiimonas chitiniclastica can cause several diseases in humans, including sepsis and bacteremia. Although the pathogenesis is not fully understood, the bacterium is thought to enter traumatic skin lesions via fly larvae, resulting in severe myiasis and/or wound contamination. Infections are typically associated with, but not limited to, infestation of an open wound by fly larvae, poor sanitary conditions, cardiovascular disease, substance abuse, and osteomyelitis. W. chitiniclastica is generally sensitive to a broad spectrum of antibiotics with the exception of fosfomycin. However, increasing drug resistance has been observed and its development should be monitored with caution. In this review, we summarize the currently available knowledge and evaluate it from both a clinical and a genomic perspective.

Biosynthetic gene cluster profiling from North Java Sea Virgibacillus salarius reveals hidden potential metabolites

Virgibacillus salarius 19.PP.SC1.6 is a coral symbiont isolated from Indonesia's North Java Sea; it has the ability to produce secondary metabolites that provide survival advantages and biological functions, such as ectoine, which is synthesized by an ectoine gene cluster. Apart from being an osmoprotectant for bacteria, ectoine is also known as a chemical chaperone with numerous biological activities such as maintaining protein stability, which makes ectoine in high demand in the market industry and makes it beneficial to investigate V. salarius ectoine. However, there has been no research on genome-based secondary metabolite and ectoine gene cluster characterization from Indonesian marine V. salarius. In this study, we performed a genomic analysis and ectoine identification of V. salarius. A high-quality draft genome with total size of 4.45 Mb and 4426 coding sequence (CDS) was characterized and then mapped into the Cluster of Orthologous Groups (COG) category. The genus Virgibacillus has an "open" pangenome type with total of 18 genomic islands inside the V. salarius 19.PP.SC1.6 genome. There were seven clusters of secondary metabolite-producing genes found, with a total of 80 genes classified as NRPS, PKS (type III), terpenes, and ectoine biosynthetic related genes. The ectoine gene cluster forms one operon consists of ectABC gene with 2190 bp gene cluster length, and is successfully characterized. The presence of ectoine in V. salarius was confirmed using UPLC-MS/MS operated in Multiple Reaction Monitoring (MRM) mode, which indicates that V. salarius has an intact ectoine gene clusters and is capable of producing ectoine as compatible solutes.

Identification of potential novel therapeutic drug target against Elizabethkingia anophelis by integrative pan and subtractive genomic analysis: An in silico approach

Elizabethkingia anophelis is a human pathogen responsible for severe nosocomial infections in neonates and immunocompromised patients. The significantly higher mortality rate from E. anophelis infections and the lack of available regimens highlight the critical need to explore novel drug targets. The current study investigated effective novel drug targets by employing a comprehensive in silico subtractive genomic approach integrated with pangenomic analysis of E. anophelis strains. A total of 2809 core genomic proteins were found by pangenomic analysis of non-paralogous proteins. Subsequently, 156 pathogen-specific, 442 choke point, 202 virulence factor, 53 antibiotic resistant and 119 host-pathogen interacting proteins were identified in E. anophelis. By subtractive genomic approach, at first 791 proteins were found to be indispensable for the survival of E. anophelis. 558 and 315 proteins were detected as non-homologous to human and gut microflora respectively. Following that 245 cytoplasmic, 245 novel, and 23 broad-spectrum targets were selected and finally four proteins were considered as potential therapeutic targets of E. anophelis based on highest degree score in PPI network. Among those, three proteins were subjected to molecular docking and subsequent MD simulation as one protein did not contain a plausible binding pocket with sufficient surface area and volume. All the complexes were found to be stable and compact in 100 ns molecular dynamics simulation studies as measured by RMSD, RMSF, and Rg. These three short-listed targets identified in this study may lead to the development of novel antimicrobials capable of curing infections and pave the way to prevent and control the disease progression caused by the deadly agent E. anophelis.

Genomic evidence of genetic diversity and functional evolution in Flavobacterium columnare

Flavobacterium columnare is the causative agent of columnaris disease in freshwater fish. Columnaris disease can cause heavy economic losses in aquaculture. In this study, whole-genome sequencing was used to characterize this pathogen. F. columnare isolate AH-01 had a circular chromosome and plasmid that encoded a total of 3,022 genes. Isolate GX-01 only had a circular chromosome and encoded 2,965 genes. Genomic islands, prophage regions, and CRISPR/Cas systems were identified in both genomes. Both genomes presented evidence of gene variation and horizontal transfer, both of which are the essential components of genetic diversity, genome plasticity, and functional evolution. Single-gene phylogeny and comparative genome analyses were performed to investigate the variation and evolution of this pathogen. Genetic analysis of 16S rRNA and housekeeping gene sequences significantly clustered 55 F. columnare isolates into four clades. The intragroup identity of the 16S rRNA gene exceeded 99%, while the intergroup identity was below the species delineation threshold. We discovered significant translocation, inversion, and rearrangement events that influenced local synteny within each group. Notably, the observed alignments varied considerably among all the studied groups. The core genomes of all strains with available sequences comprised 747 genes, corresponding to approximately 25% of the genome. Core genome multilocus sequence typing, genome-wide orthology and phylogenetic analyses, and average nucleotide identity suggested that the currently existing F. columnare was an assemblage of several distinct species, with levels of divergence at least equivalent to those between recognized bacterial species. The present investigation provided genomic evidence of gene variation and horizontal transfer, which were the basis of genetic diversity, genome plasticity, and functional evolution. The findings supported a proposed new taxonomic perspective on F. columnare.

Insights into early evolutionary adaptations of the Akkermansia genus to the vertebrate gut

Akkermansia, a relevant mucin degrader from the vertebrate gut microbiota, is a member of the deeply branched Verrucomicrobiota, as well as the only known member of this phylum to be described as inhabitants of the gut. Only a few Akkermansia species have been officially described so far, although there is genomic evidence addressing the existence of more species-level variants for this genus. This niche specialization makes Akkermansia an interesting model for studying the evolution of microorganisms to their adaptation to the gastrointestinal tract environment, including which kind of functions were gained when the Akkermansia genus originated or how the evolutionary pressure functions over those genes. In order to gain more insight into Akkermansia adaptations to the gastrointestinal tract niche, we performed a phylogenomic analysis of 367 high-quality Akkermansia isolates and metagenome-assembled genomes, in addition to other members of Verrucomicrobiota. This work was focused on three aspects: the definition of Akkermansia genomic species clusters and the calculation and functional characterization of the pangenome for the most represented species; the evolutionary relationship between Akkermansia and their closest relatives from Verrucomicrobiota, defining the gene families which were gained or lost during the emergence of the last Akkermansia common ancestor (LAkkCA) and; the evaluation of the evolutionary pressure metrics for each relevant gene family of main Akkermansia species. This analysis found 25 Akkermansia genomic species clusters distributed in two main clades, divergent from their non-Akkermansia relatives. Pangenome analyses suggest that Akkermansia species have open pangenomes, and the gene gain/loss model indicates that genes associated with mucin degradation (both glycoside hydrolases and peptidases), (micro)aerobic metabolism, surface interaction, and adhesion were part of LAkkCA. Specifically, mucin degradation is a very ancestral innovation involved in the origin of Akkermansia. Horizontal gene transfer detection suggests that Akkermansia could receive genes mostly from unknown sources or from other Gram-negative gut bacteria. Evolutionary metrics suggest that Akkemansia species evolved differently, and even some conserved genes suffered different evolutionary pressures among clades. These results suggest a complex evolutionary landscape of the genus and indicate that mucin degradation could be an essential feature in Akkermansia evolution as a symbiotic species.

Acquisition of Type I methyltransferase via horizontal gene transfer increases the drug resistance of Aeromonas veronii

Aeromonas veronii is an opportunistic pathogen that affects both fish and mammals, including humans, leading to bacteraemia, sepsis, meningitis and even death. The increasing virulence and drug resistance of A. veronii are of significant concern and pose a severe risk to public safety. The Type I restriction-modification (RM) system, which functions as a bacterial defence mechanism, can influence gene expression through DNA methylation. However, little research has been conducted to explore its origin, evolutionary path, and relationship to virulence and drug resistance in A. veronii. In this study, we analysed the pan-genome of 233 A. veronii strains, and the results indicated that it was 'open', meaning that A. veronii has acquired additional genes from other species. This suggested that A. veronii had the potential to adapt and evolve rapidly, which might have contributed to its drug resistance. One Type I methyltransferase (MTase) and two complete Type I RM systems were identified, namely AveC4I, AveC4II and AveC4III in A. veronii strain C4, respectively. Notably, AveC4I was exclusive to A. veronii C4. Phylogenetic analysis revealed that AveC4I was derived from horizontal gene transfer from Thiocystis violascens and exchanged genes with the human pathogen Comamonas kerstersii. Single molecule real-time sequencing was applied to identify the motif methylated by AveC4I, which was unique and not recognized by any reported MTases in the REBASE database. We also annotated the functions and pathways of the genes containing the motif, revealing that AveC4I may control drug resistance in A. veronii C4. Our findings provide new insight on the mechanisms underlying drug resistance in pathogenic bacteria. By identifying the specific genes and pathways affected by AveC4I, this study may aid in the development of new therapeutic approaches to combat A. veronii infections.

Comparative Genomic Analysis Reveals Genetic Diversity and Pathogenic Potential of Haemophilus seminalis and Emended Description of Haemophilus seminalis

Haemophilus seminalis is a newly proposed species that is phylogenetically related to Haemophilus haemolyticus. The distribution of H. seminalis in the human population, its genomic diversity, and its pathogenic potential are still unclear. This study reports the finding of our comparative genomic analyses of four newly isolated Haemophilus strains (SZY H8, SZY H35, SZY H36, and SZY H68) from human sputum specimens (Guangzhou, China) along with the publicly available genomes of other phylogenetically related Haemophilus species. Based on pairwise comparisons of the 16S rRNA gene sequences, the four isolates showed <98.65% sequence identity to the type strains of all known Haemophilus species but were identified as belonging to H. seminalis, based on comparable phenotypic and genotypic features. Additionally, the four isolates showed high genome-genome relatedness indices (>95% ANI values) with 17 strains that were previously identified as either "Haemophilus intermedius" or hemin (X-factor)-independent H. haemolyticus and therefore required a more detailed classification study. Phylogenetically, these isolates, along with the two previously described H. seminalis isolates (a total of 23 isolates), shared a highly homologous lineage that is distinct from the clades of the main H. haemolyticus and Haemophilus influenzae strains. These isolates present an open pangenome with multiple virulence genes. Notably, all 23 isolates have a functional heme biosynthesis pathway that is similar to that of Haemophilus parainfluenzae. The phenotype of hemin (X-factor) independence and the analysis of the ispD, pepG, and moeA genes can be used to distinguish these isolates from H. haemolyticus and H. influenzae. Based on the above findings, we propose a reclassification for all "H. intermedius" and two H. haemolyticus isolates belonging to H. seminalis with an emended description of H. seminalis. This study provides a more accurate identification of Haemophilus isolates for use in the clinical laboratory and a better understanding of the clinical significance and genetic diversity in human environments. IMPORTANCE As a versatile opportunistic pathogen, the accurate identification of Haemophilus species is a challenge in clinical practice. In this study, we characterized the phenotypic and genotypic features of four H. seminalis strains that were isolated from human sputum specimens and propose the "H. intermedius" and hemin (X-factor)-independent H. haemolyticus isolates as belonging to H. seminalis. The prediction of virulence-related genes indicates that H. seminalis isolates carry several virulence genes that are likely to play an important role in its pathogenicity. In addition, we depict that the genes ispD, pepG, and moeA can be used as biomarkers for distinguishing H. seminalis from H. haemolyticus and H. influenzae. Our findings provide some insights into the identification, epidemiology, genetic diversity, pathogenic potential, and antimicrobial resistance of the newly proposed H. seminalis.

Pangenome-driven insights into nitrogen metabolic characteristics of Citrobacter portucalensis strain AAK_AS5 associated with wastewater nitrogen removal

Nitrogen metabolism in the genus Citrobacter is very poorly studied despite its several implications in wastewater treatment. In the current study, Citrobacter portucalensis strain AAK_AS5 was assessed for remediation of simulated wastewater supplemented with different inorganic nitrogen sources. Combination of (NH₄)₂SO₄ with KNO₃ was the most preferred for achieving high growth density followed by (NH₄)₂SO₄ and KNO₃ alone. This was in agreement with highest ammonical nitrogen removal of 92.9% in the presence of combined nitrogen sources and the corresponding nitrate nitrogen removal of 93% in the presence of KNO₃. Furthermore, these removal capacities were validated by investigating the uniqueness and the spread of metabolic features through pan-genomic approach that revealed the largest number of unique genes (2097) and accessory genes (705) in strain AAK_AS5. Of the total 44 different types of nitrogen metabolism-related genes, 39 genes were associated with the core genome, while 5 genes such as gltI, nasA, nasR, nrtA, and ntrC uniquely belonged to the accessory genome. Strain AAK_AS5 possessed three major nitrate removal pathways viz., assimilatory and dissimilatory nitrate reduction to ammonia (ANRA & DNRA), and denitrification; however, the absence of nitrification was compensated by ammonia assimilation catalyzed by gene products of the GDH and GS-GOGAT pathways. narGHIJ encoding the respiratory nitrate reductase was commonly identified in all the studied genomes, while genes such as nirK, norB, and nosZ were uniquely present in the strain AAK_AS5 only. A markedly different genetic content and metabolic diversity between the strains reflected their adaptive evolution in the environment thus highlighting the significance of C. portucalensis AAK_AS5 for potential application in nitrogen removal from wastewater.

Genome Organization of Four Brazilian Xanthomonas albilineans Strains Does Not Correlate with Aggressiveness

An integrative approach combining genomics, transcriptomics, and cell biology is presented to address leaf scald disease, a major problem for the sugarcane industry. To gain insight into the biology of the causal agent, the complete genome sequences of four Brazilian Xanthomonas albilineans strains with differing virulence capabilities are presented and compared to the GPEPC73 reference strain and FJ1. Based on the aggressiveness index, different strains were compared: Xa04 and Xa11 are highly aggressive, Xa26 is intermediate, and Xa21 is the least, while, based on genome structure, Xa04 shares most of its genomic features with Xa26, and Xa11 share most of its genomic features with Xa21. In addition to presenting more clustered regularly interspaced short palindromic repeats (CRISPR) clusters, four more novel prophage insertions are present than the previously sequenced GPEPC73 and FJ1 strains. Incorporating the aggressiveness index and in vitro cell biology into these genome features indicates that disease establishment is not a result of a single determinant factor, as in most other Xanthomonas species. The Brazilian strains lack the previously described plasmids but present more prophage regions. In pairs, the most virulent and the least virulent share unique prophages. In vitro transcriptomics shed light on the 54 most highly expressed genes among the 4 strains compared to ribosomal proteins (RPs), of these, 3 outer membrane proteins. Finally, comparative albicidin inhibition rings and in vitro growth curves of the four strains also do not correlate with pathogenicity. In conclusion, the results disclose that leaf scald disease is not associated with a single shared characteristic between the most or the least pathogenic strains. IMPORTANCE An integrative approach is presented which combines genomics, transcriptomics, and cell biology to address leaf scald disease. The results presented here disclose that the disease is not associated with a single shared characteristic between the most pathogenic strains or a unique genomic pattern. Sequence data from four Brazilian strains are presented that differ in pathogenicity index: Xa04 and Xa11 are highly virulent, Xa26 is intermediate, and Xa21 is the least pathogenic strain, while, based on genome structure, Xa04 shares with Xa26, and Xa11 shares with X21 most of the genome features. Other than presenting more CRISPR clusters and prophages than the previously sequenced strains, the integration of aggressiveness and cell biology points out that disease establishment is not a result of a single determinant factor as in other xanthomonads.

Brevundimonas brasiliensis sp. nov.: a New Multidrug-Resistant Species Isolated from a Patient in Brazil

To increase knowledge on Brevundimonas pathogens, we conducted in-depth genomic and phenotypic characterization of a Brevundimonas strain isolated from the cerebrospinal fluid of a patient admitted in a neonatal intensive care unit. The strain was identified as a member of the genus Brevundimonas based on Vitek 2 system results and 16S rRNA gene sequencing and presented a multidrug resistance profile (MDR). Several molecular and biochemical tests were used to characterize and identify the species for in-depth results. The draft genome assembly of the isolate has a total length of 3,261,074 bp and a G+C of 66.86%, similar to other species of the genus. Multilocus sequence analysis, Type (Strain) Genome Server, digital DNA-DNA hybridization, and average nucleotide identity confirmed that the Brevundimonas sp. studied represents a distinct species, for which we propose the name Brevundimonas brasiliensis sp. nov. In silico analysis detected antimicrobial resistance genes (AMRGs) mediating resistance to β-lactams (penP, blaTEM-16, and blaBKC-1) and aminoglycosides [strA, strB, aac(6')-Ib, and aac(6')-Il]. We also found AMRGs encoding the AcrAB efflux pump that confers resistance to a broad spectrum of antibiotics. Colistin and quinolone resistance can be attributed to mutation in qseC and/or phoP and GyrA/GyrB, respectively. The Brevundimonas brasiliensis sp. nov. genome contained copies of type IV secretion system (T4SS)-type integrative and conjugative elements (ICEs); integrative mobilizable elements (IME); and Tn3-type and IS3, IS6, IS5, and IS1380 families, suggesting an important role in the development and dissemination of antibiotic resistance. The isolate presented a range of virulence-associated genes related to biofilm formation, adhesion, and invasion that can be relevant for its pathogenicity. Our findings provide a wealth of data to hinder the transmission of MDR Brevundimonas and highlight the need for monitoring and identifying new bacterial species in hospital environments. IMPORTANCE Brevundimonas species is considered an opportunistic human pathogen that can cause multiple types of invasive and severe infections in patients with underlying pathologies. Treatment of these pathogens has become a major challenge because many isolates are resistant to most antibiotics used in clinical practice. Furthermore, there are no consistent therapeutic results demonstrating the efficacy of antibacterial agents. Although considered a rare pathogen, recent studies have provided evidence of the emergence of Brevundimonas in clinical settings. Hence, we identified a novel pathogenic bacterium, Brevundimonas brasiliensis sp. nov., that presented a multidrug resistance (MDR) profile and carried diverse genes related to drug resistance, virulence, and mobile genetic elements. Such data can serve as a baseline for understanding the genomic diversity, adaptation, evolution, and pathogenicity of MDR Brevundimonas.

Using bacterial pan-genome-based feature selection approach to improve the prediction of minimum inhibitory concentration (MIC)

Background: Predicting the resistance profiles of antimicrobial resistance (AMR) pathogens is becoming more and more important in treating infectious diseases. Various attempts have been made to build machine learning models to classify resistant or susceptible pathogens based on either known antimicrobial resistance genes or the entire gene set. However, the phenotypic annotations are translated from minimum inhibitory concentration (MIC), which is the lowest concentration of antibiotic drugs in inhibiting certain pathogenic strains. Since the MIC breakpoints that classify a strain to be resistant or susceptible to specific antibiotic drug may be revised by governing institutes, we refrained from translating these MIC values into the categories "susceptible" or "resistant" but instead attempted to predict the MIC values using machine learning approaches. Results: By applying a machine learning feature selection approach on a Salmonella enterica pan-genome, in which the protein sequences were clustered to identify highly similar gene families, we showed that the selected features (genes) performed better than known AMR genes, and that models built on the selected genes achieved very accurate MIC prediction. Functional analysis revealed that about half of the selected genes were annotated as hypothetical proteins (i.e., with unknown functional roles), and that only a small portion of known AMR genes were among the selected genes, indicating that applying feature selection on the entire gene set has the potential of uncovering novel genes that may be associated with and may contribute to pathogenic antimicrobial resistances. Conclusion: The application of the pan-genome-based machine learning approach was indeed capable of predicting MIC values with very high accuracy. The feature selection process may also identify novel AMR genes for inferring bacterial antimicrobial resistance phenotypes.

Comparative Genomics of Pseudomonas aeruginosa Strains Isolated from Different Ecological Niches

The Pseudomonas aeruginosa genome can change to adapt to different ecological niches. We compared four genomes from a Mexican hospital and 59 genomes from GenBank from different niches, such as urine, sputum, and environmental. The ST analysis showed that high-risk STs (ST235, ST773, and ST27) were present in the genomes of the three niches from GenBank, and the STs of Mexican genomes (ST167, ST2731, and ST549) differed from the GenBank genomes. Phylogenetic analysis showed that the genomes were clustering according to their ST and not their niche. When analyzing the genomic content, we observed that environmental genomes had genes involved in adapting to the environment not found in the clinics and that their mechanisms of resistance were mutations in antibiotic resistance-related genes. In contrast, clinical genomes from GenBank had resistance genes, in mobile/mobilizable genetic elements in the chromosome, except for the Mexican genomes that carried them mostly in plasmids. This was related to the presence of CRISPR-Cas and anti-CRISPR; however, Mexican strains only had plasmids and CRISPR-Cas. bla_OXA-488 (a variant of bla_OXA50) with higher activity against carbapenems was more prevalent in sputum genomes. The virulome analysis showed that exoS was most prevalent in the genomes of urinary samples and exoU and pldA in sputum samples. This study provides evidence regarding the genetic variability among P. aeruginosa isolated from different niches.

Characterization and Association of Rips Repertoire to Host Range of Novel Ralstonia solanacearum Strains by In Silico Approaches

Ralstonia solanacearum species complex (RSSC) cause several phytobacteriosis in many economically important crops around the globe, especially in the tropics. In Brazil, phylotypes I and II cause bacterial wilt (BW) and are indistinguishable by classical microbiological and phytopathological methods, while Moko disease is caused only by phylotype II strains. Type III effectors of RSSC (Rips) are key molecular actors regarding pathogenesis and are associated with specificity to some hosts. In this study, we sequenced and characterized 14 newly RSSC isolates from Brazil's Northern and Northeastern regions, including BW and Moko ecotypes. Virulence and resistance sequences were annotated, and the Rips repertoire was predicted. Confirming previous studies, RSSC pangenome is open as α≅0.77. Genomic information regarding these isolates matches those for R. solanacearum in NCBI. All of them fit in phylotype II with a similarity above 96%, with five isolates in phylotype IIB and nine in phylotype IIA. Almost all R. solanacearum genomes in NCBI are actually from other species in RSSC. Rips repertoire of Moko IIB was more homogeneous, except for isolate B4, which presented ten non-shared Rips. Rips repertoire of phylotype IIA was more diverse in both Moko and BW, with 43 common shared Rips among all 14 isolates. New BW isolates shared more Rips with Moko IIA and Moko IIB than with other public BW genome isolates from Brazil. Rips not shared with other isolates might contribute to individual virulence, but commonly shared Rips are good avirulence candidates. The high number of Rips shared by new Moko and BW isolates suggests they are actually Moko isolates infecting solanaceous hosts. Finally, infection assays and Rips expression on different hosts are needed to better elucidate the association between Rips repertoire and host specificities.

Comparative genomic analysis of Stenotrophomonas maltophilia unravels their genetic variations and versatility trait

Stenotrophomonas maltophilia is a species with immensely broad phenotypic and genotypic diversity that could widely distribute in natural and clinical environments. However, little attention has been paid to reveal their genome plasticity to diverse environments. In the present study, a comparative genomic analysis of S. maltophilia isolated from clinical and natural sources was systematically explored its genetic diversity of 42 sequenced genomes. The results showed that S. maltophilia owned an open pan-genome and had strong adaptability to different environments. A total of 1612 core genes were existed with an average of 39.43% of each genome, and the shared core genes might be necessary to maintain the basic characteristics of those S. maltophilia strains. Based on the results of the phylogenetic tree, the ANI value, and the distribution of accessory genes, genes associated with the fundamental process of those strains from the same habitat were found to be mostly conserved in evolution. Isolates from the same habitat had a high degree of similarity in COG category, and the most significant KEGG pathways were mainly involved in carbohydrate and amino acid metabolism, indicating that genes related to essential processes were mostly conserved in evolution for the clinical and environmental settings. Meanwhile, the number of resistance and efflux pump gene was significantly higher in the clinical setting than that of in the environmental setting. Collectively, this study highlights the evolutionary relationships of S. maltophilia isolated from clinical and environmental sources, shedding new light on its genomic diversity.

Bipartite Genomes in Enterobacterales: Independent Origins of Chromids, Elevated Openness and Donors of Horizontally Transferred Genes

Multipartite bacteria have one chromosome and one or more chromid. Chromids are believed to have properties that enhance genomic flexibility, making them a favored integration site for new genes. However, the mechanism by which chromosomes and chromids jointly contribute to this flexibility is not clear. To shed light on this, we analyzed the openness of chromosomes and chromids of the two bacteria, Vibrio and Pseudoalteromonas, both which belong to the Enterobacterales order of Gammaproteobacteria, and compared the genomic openness with that of monopartite genomes in the same order. We applied pangenome analysis, codon usage analysis and the HGTector software to detect horizontally transferred genes. Our findings suggest that the chromids of Vibrio and Pseudoalteromonas originated from two separate plasmid acquisition events. Bipartite genomes were found to be more open compared to monopartite. We found that the shell and cloud pangene categories drive the openness of bipartite genomes in Vibrio and Pseudoalteromonas. Based on this and our two recent studies, we propose a hypothesis that explains how chromids and the chromosome terminus region contribute to the genomic plasticity of bipartite genomes.

The Emerging Nosocomial Pathogen Klebsiella michiganensis: Genetic Analysis of a KPC-3 Producing Strain Isolated from Venus Clam

The recovery and characterization of a multidrug-resistant, KPC-3-producing Klebsiella michiganensis that was obtained from Venus clam samples is reported in this study. A whole-genome sequencing (WGS) analysis using Illumina and Nanopore technologies of the K. michiganensis 23999A2 isolate revealed that the strain belonged to the new sequence type 382 (ST382) and carried seven plasmid replicon sequences, including four IncF type plasmids (FII, FIIY, FIIk, and FIB), one IncHI1 plasmid, and two Col plasmids. The FIB and FIIk plasmids showed high homology to each other and to multireplicon pKpQIL-like plasmids that are found in epidemic KPC-K. pneumoniae clones worldwide. The strain carried multiple β-lactamase genes on the IncF plasmids: blaOXA-9 and blaTEM-1A on FIB, blaKPC-3 inserted in a Tn4401a on FIIK, and blaSHV-12 on FIIY. The IncHI1-ST11 harbored no resistance gene. The curing of the strain caused the loss of all of the bla genes and a rearrangement of the IncF plasmids. Conjugal transfer of the blaOXA-9, blaTEM-1A and blaKPC-3 genes occurred at a frequency of 5 × 10-7, using K. quasipneumoniae as a recipient, and all of the bla genes were transferred through a pKpQIL that originated from the recombination of the FIB and FIIk plasmids of the donor. A comparison with 31 K. michiganensis genomes that are available in the NCBI database showed that the closest phylogenetic relatives of K. michiganensis 23999A2 are an environmental isolate from soil in South Korea and a clinical isolate from human sputum in Japan. Finally, a pan-genome analysis showed a large accessory genome of the strain as well as the great genomic plasticity of the K. michiganensis species. IMPORTANCE Klebsiella michiganensis is an emerging nosocomial pathogen, and, so far, few studies describe isolates of clinical origin in the environment. This study contributes to the understanding of how the dissemination of carbapenem-resistance outside the hospital setting may be related to the circulation of pKpQIL-like plasmids that are derived from epidemic Klebsiella pneumoniae strains. The recovery of a carbapenem-resistant isolate in clams is of great concern, as bivalves could represent vehicles of transmission of pathogens and resistance genes to humans via the food chain. The study demonstrates the plasticity of K. michiganensis genome, which is probably useful to multiple environment adaptation and to the evolution of the species.

Integration of text mining and biological network analysis: Identification of essential genes in sulfate-reducing bacteria

The growth and survival of an organism in a particular environment is highly depends on the certain indispensable genes, termed as essential genes. Sulfate-reducing bacteria (SRB) are obligate anaerobes which thrives on sulfate reduction for its energy requirements. The present study used Oleidesulfovibrio alaskensis G20 (OA G20) as a model SRB to categorize the essential genes based on their key metabolic pathways. Herein, we reported a feedback loop framework for gene of interest discovery, from bio-problem to gene set of interest, leveraging expert annotation with computational prediction. Defined bio-problem was applied to retrieve the genes of SRB from literature databases (PubMed, and PubMed Central) and annotated them to the genome of OA G20. Retrieved gene list was further used to enrich protein-protein interaction and was corroborated to the pangenome analysis, to categorize the enriched gene sets and the respective pathways under essential and non-essential. Interestingly, the sat gene (dde_2265) from the sulfur metabolism was the bridging gene between all the enriched pathways. Gene clusters involved in essential pathways were linked with the genes from seleno-compound metabolism, amino acid metabolism, secondary metabolite synthesis, and cofactor biosynthesis. Furthermore, pangenome analysis demonstrated the gene distribution, where 69.83% of the 116 enriched genes were mapped under "persistent," inferring the essentiality of these genes. Likewise, 21.55% of the enriched genes, which involves specially the formate dehydrogenases and metallic hydrogenases, appeared under "shell." Our methodology suggested that semi-automated text mining and network analysis may play a crucial role in deciphering the previously unexplored genes and key mechanisms which can help to generate a baseline prior to perform any experimental studies.

Comparative Genome Analysis of 19 Trueperella pyogenes Strains Originating from Different Animal Species Reveal a Genetically Diverse Open Pan-Genome

Trueperella pyogenes is a Gram-positive opportunistic pathogen that causes severe cases of mastitis, metritis, and pneumonia in a wide range of animals, resulting in significant economic losses. Although little is known about the virulence factors involved in the disease pathogenesis, a comprehensive comparative genome analysis of T. pyogenes genomes has not been performed till date. Hence, present investigation was carried out to characterize and compare 19 T. pyogenes genomes originating in different geographical origins including the draftgenome of the first Indian origin strain T. pyogenes Bu5. Additionally, candidate virulence determinants that could be crucial for their pathogenesis were also detected and analyzed by using various bioinformatics tools. The pan-genome calculations revealed an open pan-genome of T. pyogenes. In addition, an inventory of virulence related genes, 190 genomic islands, 31 prophage sequences, and 40 antibiotic resistance genes that could play a significant role in organism's pathogenicity were detected. The core-genome based phylogeny of T. pyogenes demonstrates a polyphyletic, host-associated group with a high degree of genomic diversity. The identified core-genome can be further used for screening of drug and vaccine targets. The investigation has provided unique insights into pan-genome, virulome, mobiliome, and resistome of T. pyogenes genomes and laid the foundation for future investigations.

Genetic Features of Mycobacterium avium subsp. paratuberculosis Strains Circulating in the West of France Deciphered by Whole-Genome Sequencing

Paratuberculosis is a chronic infection of the intestine, mainly the ileum, caused by Mycobacterium avium subsp. paratuberculosis in cattle and other ruminants. This enzootic disease is present worldwide and has a negative impact on the dairy cattle industry. For this subspecies, the current genotyping tools do not provide the needed resolution to investigate the genetic diversity of closely related strains. These limitations can be overcome by the application of whole-genome sequencing (WGS), particularly for clonal populations such as M. avium subsp. paratuberculosis. The purpose of the present study was to undertake a WGS analysis with a panel of 200 animal field M. avium subsp. paratuberculosis strains selected based on a previous large-scale longitudinal study of Prim'Holstein and Normande dairy breeds naturally infected with M. avium subsp. paratuberculosis in the West of France. The pangenome analysis revealed that M. avium subsp. paratuberculosis has a closed pangenome. The phylogeny, based on alignment of 2,786 nonhomoplasic single nucleotide polymorphisms (SNPs), showed that the strain population is structured into three clades independently of the cattle breed or geographic distribution. The increased resolution of phylogeny obtained by WGS confirmed the homoplasic nature of the markers variable-number tandem repeat (VNTR) and short sequence repeat (SSR) used for M. avium subsp. paratuberculosis genotyping. These phylogenetic data also revealed independent introductions of the different genotypes in two main waves since at least 2003. WGS applied to this sampling demonstrated the presence of mixed infections in herds and at the individual animal level. Collectively, the phylogeny results inferred with French isolates compared to M. avium subsp. paratuberculosis isolates from around the world suggest introductions of M. avium subsp. paratuberculosis genotypes through the animal trade. Relationships between genetic traits and epidemiological data can now be investigated to better understand transmission dynamics of the disease. IMPORTANCE Mycobacterium avium subsp. paratuberculosis causes Johne's disease in ruminants, which is present worldwide and has significant negative impacts on the dairy cattle industry and animal welfare. Prevention and control of M. avium subsp. paratuberculosis infection are hampered by knowledge gaps in strain virulence, genotype distribution, and transmission dynamics. This work has revealed new insights into M. avium subsp. paratuberculosis strains currently circulating in western France and how they are related to strains circulating globally. We applied whole-genome sequencing (WGS) to obtain comprehensive information on genome evolution and discrimination of closely related strains. This approach revealed the history of M. avium subsp. paratuberculosis infection in France, refined the pangenomic characteristics of M. avium subsp. paratuberculosis, and demonstrated the existence of mixed infection in animals. Finally, this study identified predominant genotypes, which allow a better understanding of disease transmission dynamics. This information will facilitate tracking of this pathogen on farms and across agricultural regions, thus informing transmission pathways and disease control points.

Navigating and expanding the roadmap of natural product genome mining tools

Natural products are structurally highly diverse and exhibit a wide array of biological activities. As a result, they serve as an important source of new drug leads. Traditionally, natural products have been discovered by bioactivity-guided fractionation. The advent of genome sequencing technology has resulted in the introduction of an alternative approach towards novel natural product scaffolds: Genome mining. Genome mining is an in-silico natural product discovery strategy in which sequenced genomes are analyzed for the potential of the associated organism to produce natural products. Seemingly universal biosynthetic principles have been deciphered for most natural product classes that are used to detect natural product biosynthetic gene clusters using pathway-encoded conserved key enzymes, domains, or motifs as bait. Several generations of highly sophisticated tools have been developed for the biosynthetic rule-based identification of natural product gene clusters. Apart from these hard-coded algorithms, multiple tools that use machine learning-based approaches have been designed to complement the existing genome mining tool set and focus on natural product gene clusters that lack genes with conserved signature sequences. In this perspective, we take a closer look at state-of-the-art genome mining tools that are based on either hard-coded rules or machine learning algorithms, with an emphasis on the confidence of their predictions and potential to identify non-canonical natural product biosynthetic gene clusters. We highlight the genome mining pipelines' current strengths and limitations by contrasting their advantages and disadvantages. Moreover, we introduce two indirect biosynthetic gene cluster identification strategies that complement current workflows. The combination of all genome mining approaches will pave the way towards a more comprehensive understanding of the full biosynthetic repertoire encoded in microbial genome sequences.

Nitrosophilus kaiyonis sp. nov., a hydrogen-, sulfur- and thiosulfate-oxidizing chemolithoautotroph within "Campylobacteria" isolated from a deep-sea hydrothermal vent in the Mid-Okinawa Trough

A novel bacterium, strain MOT50^T, was isolated from the chimney structure at the Iheya North field in the Mid-Okinawa Trough. The cells were motile short rods with a single polar flagellum. Growth was observed between 40 and 65 ℃ (optimum, 52 ℃), at pH values between 5.0 and 7.1 (optimum, pH 6.1) and in the presence of 2.0-4.0% NaCl (optimum, 2.5%). The isolates utilized molecular hydrogen, thiosulfate, or elemental sulfur as the sole electron donor. Thiosulfate, elemental sulfur, nitrate, and molecular oxygen are utilized as the sole electron acceptor. Ammonium is required as a nitrogen source. Thiosulfate, elemental sulfur, sulfate, or sulfite serves as a sulfur source for growth. The G + C content of the genomic DNA was 28.9%. Phylogenetic analysis based on the 16S rRNA gene sequences indicated that strain MOT50^T belonged to the genus Nitrosophilus of the class "Campylobacteria", and its closest relative was Nitrosophilus labii HRV44^T (97.20%). On the basis of the phylogenetic, physiological, and molecular characteristics, it is proposed that the organism represents a novel species within the genus Nitrosophilus, Nitrosophilus kaiyonis sp. nov. The type strain is MOT50^T (= JCM 39187^T = KCTC 25251^T).

Characterization and Identification of Probiotic Features in Lacticaseibacillus Paracasei Using a Comparative Genomic Analysis Approach

Lacticaseibacillus paracasei species are widely used for their health-promoting properties in food and agricultural applications. These bacteria have been isolated from various habitats such as the oral cavity, cereals, vegetables, meats, and dairy products conferring them the ability to consume different carbohydrates. Two subspecies are recognized, Lacticaseibacillus paracasei subsp. paracasei and Lacticaseibacillus paracasei subsp. tolerans according to their acid production from carbohydrates. Some strains are currently used as probiotics. In this study, we performed a comparative genomic analysis of 181 genomes of the Lacticaseibacillus paracasei species to reveal genomic differences at the subspecies level and to reveal adaptive and probiotic features, and special emphasis is given to inulin consumption. No clear distinction at the subspecies level for L. paracasei was shown using a phylogenetic tree with orthologous genes from the core-genome set. In general, a good correlation was observed between genomic distance and isolation origin, suggesting that L. paracasei strains are adapted to their natural habitat, giving rise to genetic differences at the genomic level. A low frequency of undesirable characteristics such as plasmids, prophages, antibiotic resistance genes, absence of virulence factors, and frequent bacteriocin production supports these species being good candidates for use as probiotics. Lastly, we found that the inulin gene cluster in L. paracasei strains seems to differ slightly in the presence or absence of some genes but maintains a core defined by at least three fructose-PTS proteins, one hypothetical protein, and extracellular β-fructosidase. Finally, we conclude that further work has to be done for L. paracasei subspecies classification. Improving outgroup selection criteria is a key factor for their correct subspecies assignation.

Specialized Metabolism of Gordonia Genus: An Integrated Survey on Chemodiversity Combined with a Comparative Genomics-Based Analysis

Members of the phylum Actinomycetota (formerly Actinobacteria) have historically been the most prolific providers of small bioactive molecules. Although the genus Streptomyces is the best-known member for this issue, other genera, such as Gordonia, have shown interesting potential in their specialized metabolism. Thus, we combined herein the result of a comprehensive literature survey on metabolites derived from Gordonia strains with a comparative genomic analysis to examine the potential of the specialized metabolism of the genus Gordonia. Thirty Gordonia-derived compounds of different classes were gathered (i.e., alkaloids, amides, phenylpropanoids, and terpenoids), exhibiting antimicrobial and cytotoxic activities, and several were also isolated from Streptomyces (e.g., actinomycin, nocardamin, diolmycin A1). With the genome data, we estimated an open pan-genome of 57,901 genes, most of them being part of the cloud genome. Regarding the BGCs content, 531 clusters were found, including Terpenes, RiPP-like, and NRPS clusters as the most frequent clusters. Our findings demonstrated that Gordonia is a poorly studied genus in terms of its specialized metabolism production and potential applications. Nevertheless, given their BGCs content, Gordonia spp. are a valuable biological resource that could expand the chemical spectrum of the phylum Actinomycetota, involving novel BGCs for inspiring innovative outlines for synthetic biology and further use in biotechnological initiatives. Therefore, further studies and more efforts should be made to explore different environments and evaluate other bioactivities.

Impact of key parameters involved with plant-microbe interaction in context to global climate change

Anthropogenic activities have a critical influence on climate change that directly or indirectly impacts plant and microbial diversity on our planet. Due to climate change, there is an increase in the intensity and frequency of extreme environmental events such as temperature rise, drought, and precipitation. The increase in greenhouse gas emissions such as CO2, CH4, NOx, water vapor, increase in global temperature, and change in rainfall patterns have impacted soil–plant-microbe interactions, which poses a serious threat to food security. Microbes in the soil play an essential role in plants’ resilience to abiotic and biotic stressors. The soil microbial communities are sensitive and responsive to these stressors. Therefore, a systemic approach to climate adaptation will be needed which acknowledges the multidimensional nature of plant-microbe-environment interactions. In the last two scores of years, there has been an enhancement in the understanding of plant’s response to microbes at physiological, biochemical, and molecular levels due to the availability of techniques and tools. This review highlights some of the critical factors influencing plant-microbe interactions under stress. The association and response of microbe and plants as a result of several stresses such as temperature, salinity, metal toxicity, and greenhouse gases are also depicted. New tools to study the molecular complexity of these interactions, such as genomic and sequencing approaches, which provide researchers greater accuracy, reproducibility, and flexibility for exploring plant-microbe–environment interactions under a changing climate, are also discussed in the review, which will be helpful in the development of resistant crops/plants in present and future.

Comparative pan-genomic analysis of 51 Renibacterium salmoninarum indicates heterogeneity in the principal virulence factor, the 57 kDa protein

Renibacterium salmoninarum, a Gram-positive intracellular pathogen, is the causative agent of bacterial kidney disease (BKD), the impacts of which are high mortalities and economic losses for the salmon industry. This study provides novel analyses for the whole-genome sequences of 50 R. salmoninarum isolates and the reference strain ATCC 33209 using a pan-genomic approach to elucidate phylogenomic relationships and identify unique and shared genes associated with pathogenicity and infection mechanisms. Genome size varied from 3,061,638 to 3,155,332 bp; gene count from 3452 to 3580; and predicted coding sequences from 3402 to 3527. Comparative analyses revealed an open, but approaching closed, pan-genome. The pan-genome analysis recovered 4064 genes, with a core genome containing 3306 genes. Phylogenetic analysis of R. salmoninarum showed high genomic homogeneity, apart from one isolate obtained from Salmo trutta in Norway. All genomes presented the 57-kDa protein (p57). Strain ATCC 33209 and the Chilean isolates H-2 and DJ2R presented two copies of the msa gene, while the remaining isolates had one copy. The pan-genome analysis further identified differences in the number of copies and length of the signalling peptide for p57, the principal virulence factor reported for this bacterium. This heterogeneity could be associated with the secretion levels of p57, potentially influencing virulence. Additionally identified were numerous common genes related to iron uptake, the stress response and regulation, and cell signalling-all of which constitute the pathogenic repertoire of R. salmoninarum. This investigation provides information that is applicable in future studies for identifying therapeutic targets and/or for designing new strategies (e.g., vaccines) to prevent BKD infections in salmon farming.

Current Techniques to Study Beneficial Plant-Microbe Interactions

Many different experimental approaches have been applied to elaborate and study the beneficial interactions between soil bacteria and plants. Some of these methods focus on changes to the plant and others are directed towards assessing the physiology and biochemistry of the beneficial plant growth-promoting bacteria (PGPB). Here, we provide an overview of some of the current techniques that have been employed to study the interaction of plants with PGPB. These techniques include the study of plant microbiomes; the use of DNA genome sequencing to understand the genes encoded by PGPB; the use of transcriptomics, proteomics, and metabolomics to study PGPB and plant gene expression; genome editing of PGPB; encapsulation of PGPB inoculants prior to their use to treat plants; imaging of plants and PGPB; PGPB nitrogenase assays; and the use of specialized growth chambers for growing and monitoring bacterially treated plants.

Comparative Genomic Analysis of the Human Pathogen Wohlfahrtiimonas Chitiniclastica Provides Insight Into the Identification of Antimicrobial Resistance Genotypes and Potential Virulence Traits

Recent studies suggest that Wohlfahrtiimonas chitiniclastica may be the cause of several diseases in humans including sepsis and bacteremia making the bacterium as a previously underappreciated human pathogen. However, very little is known about the pathogenicity and genetic potential of W. chitiniclastica; therefore, it is necessary to conduct systematic studies to gain a deeper understanding of its virulence characteristics and treatment options. In this study, the entire genetic repertoire of all publicly available W. chitiniclastica genomes was examined including in silico characterization of bacteriophage content, antibiotic resistome, and putative virulence profile. The pan-genome of W. chitiniclastica comprises 3819 genes with 1622 core genes (43%) indicating a putative metabolic conserved species. Furthermore, in silico analysis indicated presumed resistome expansion as defined by the presence of genome-encoded transposons and bacteriophages. While macrolide resistance genes macA and macB are located within the core genome, additional antimicrobial resistance genotypes for tetracycline (tetH, tetB, and tetD), aminoglycosides (ant(2'')-Ia, aac(6')-Ia,aph(3'')-Ib, aph(3')-Ia, and aph(6)-Id)), sulfonamide (sul2), streptomycin (strA), chloramphenicol (cat3), and beta-lactamase (blaVEB) are distributed among the accessory genome. Notably, our data indicate that the type strain DSM 18708T does not encode any additional clinically relevant antibiotic resistance genes, whereas drug resistance is increasing within the W. chitiniclastica clade. This trend should be monitored with caution. To the best of our knowledge, this is the first comprehensive genome analysis of this species, providing new insights into the genome of this opportunistic human pathogen.

Genomic Characterization of Enterococcus hirae From Beef Cattle Feedlots and Associated Environmental Continuum

Enterococci are commensal bacteria of the gastrointestinal tract of humans, animals, and insects. They are also found in soil, water, and plant ecosystems. The presence of enterococci in human, animal, and environmental settings makes these bacteria ideal candidates to study antimicrobial resistance in the One-Health continuum. This study focused on Enterococcus hirae isolates (n = 4,601) predominantly isolated from beef production systems including bovine feces (n = 4,117, 89.5%), catch-basin water (n = 306, 66.5%), stockpiled bovine manure (n = 24, 0.5%), and natural water sources near feedlots (n = 145, 32%), and a few isolates from urban wastewater (n = 9, 0.2%) denoted as human-associated environmental samples. Antimicrobial susceptibility profiling of a subset (n = 1,319) of E. hirae isolates originating from beef production systems (n = 1,308) showed high resistance to tetracycline (65%) and erythromycin (57%) with 50.4% isolates harboring multi-drug resistance, whereas urban wastewater isolates (n = 9) were resistant to nitrofurantoin (44.5%) and tigecycline (44.5%) followed by linezolid (33.3%). Genes for tetracycline (tetL, M, S/M, and O/32/O) and macrolide resistance erm(B) were frequently found in beef production isolates. Antimicrobial resistance profiles of E. hirae isolates recovered from different environmental settings appeared to reflect the kind of antimicrobial usage in beef and human sectors. Comparative genomic analysis of E. hirae isolates showed an open pan-genome that consisted of 1,427 core genes, 358 soft core genes, 1701 shell genes, and 7,969 cloud genes. Across species comparative genomic analysis conducted on E. hirae, Enterococcus faecalis and Enterococcus faecium genomes revealed that E. hirae had unique genes associated with vitamin production, cellulose, and pectin degradation, traits which may support its adaptation to the bovine digestive tract. E. faecium and E. faecalis more frequently harbored virulence genes associated with biofilm formation, iron transport, and cell adhesion, suggesting niche specificity within these species.

Comparative genomic analysis revealed genetic divergence between Bifidobacterium catenulatum subspecies present in infant versus adult guts

BACKGROUND

The two subspecies of Bifidobacterium catenulatum, B. catenulatum subsp. kashiwanohense and B. catenulatum subsp. catenulatum, are usually from the infant and adult gut, respectively. However, the genomic analysis of their functional difference and genetic divergence has been rare. Here, 16 B. catenulatum strains, including 2 newly sequenced strains, were analysed through comparative genomics.

RESULTS

A phylogenetic tree based on 785 core genes indicated that the two subspecies of B. catenulatum were significantly separated. The comparison of genomic characteristics revealed that the two subspecies had significantly different genomic sizes (p < 0.05) but similar GC contents. The functional comparison revealed the most significant difference in genes of carbohydrate utilisation. Carbohydrate-active enzymes (CAZyme) present two clustering patterns in B. catenulatum. The B. catenulatum subsp. kashiwanohense specially including the glycoside hydrolases 95 (GH95) and carbohydrate-binding modules 51 (CBM51) families involved in the metabolism of human milk oligosaccharides (HMO) common in infants, also, the corresponding fucosylated HMO gene clusters were detected. Meanwhile, B. catenulatum subsp. catenulatum rich in GH3 may metabolise more plant-derived glycan in the adult intestine.

CONCLUSIONS

These findings provide genomic evidence of carbohydrate utilisation bias, which may be a key cause of the genetic divergence of two B. catenulatum subspecies.

To kill or to be killed: pangenome analysis of Escherichia coli strains reveals a tailocin specific for pandemic ST131

BACKGROUND

Escherichia coli (E. coli) has been one of the most studied model organisms in the history of life sciences. Initially thought just to be commensal bacteria, E. coli has shown wide phenotypic diversity including pathogenic isolates with great relevance to public health. Though pangenome analysis has been attempted several times, there is no systematic functional characterization of the E. coli subgroups according to the gene profile.

RESULTS

Systematically scanning for optimal parametrization, we have built the E. coli pangenome from 1324 complete genomes. The pangenome size is estimated to be ~25,000 gene families (GFs). Whereas the core genome diminishes as more genomes are added, the softcore genome (≥95% of strains) is stable with ~3000 GFs regardless of the total number of genomes. Apparently, the softcore genome (with a 92% or 95% generation threshold) can define the genome of a bacterial species listing the critically relevant, evolutionarily most conserved or important classes of GFs. Unsupervised clustering of common E. coli sequence types using the presence/absence GF matrix reveals distinct characteristics of E. coli phylogroups B1, B2, and E. We highlight the bi-lineage nature of B1, the variation of the secretion and of the iron acquisition systems in ST11 (E), and the incorporation of a highly conserved prophage into the genome of ST131 (B2). The tail structure of the prophage is evolutionarily related to R2-pyocin (a tailocin) from Pseudomonas aeruginosa PAO1. We hypothesize that this molecular machinery is highly likely to play an important role in protecting its own colonies; thus, contributing towards the rapid rise of pandemic E. coli ST131.

CONCLUSIONS

This study has explored the optimized pangenome development in E. coli. We provide complete GF lists and the pangenome matrix as supplementary data for further studies. We identified biological characteristics of different E. coli subtypes, specifically for phylogroups B1, B2, and E. We found an operon-like genome region coding for a tailocin specific for ST131 strains. The latter is a potential killer weapon providing pandemic E. coli ST131 with an advantage in inter-bacterial competition and, suggestively, explains their dominance as human pathogen among E. coli strains.

Metabolism of the Genus Guyparkeria Revealed by Pangenome Analysis

Halophilic sulfur-oxidizing bacteria belonging to the genus Guyparkeria occur at both marine and terrestrial habitats. Common physiological characteristics displayed by Guyparkeria isolates have not yet been linked to the metabolic potential encoded in their genetic inventory. To provide a genetic basis for understanding the metabolism of Guyparkeria, nine genomes were compared to reveal the metabolic capabilities and adaptations. A detailed account is given on Guyparkeria’s ability to assimilate carbon by fixation, to oxidize reduced sulfur, to oxidize thiocyanate, and to cope with salinity stress.