The bacterial pangenome as a new tool for analysing pathogenic bacteria

Investigating Anthrax-Associated Virulence Genes among Archival and Contemporary Bacillus cereus Group Genomes

Bacillus anthracis causes anthrax through virulence factors encoded on two plasmids. However, non-B. anthracis organisms within the closely related, environmentally ubiquitous Bacillus cereus group (BCG) may cause an anthrax-like disease in humans through the partial adoption of anthrax-associated virulence genes, challenging the definition of anthrax disease. To elucidate these phenomena and their evolutionary past, we performed whole-genome sequencing on non-anthracis BCG isolates, including 93 archival (1967-2003) and 5 contemporary isolates (2019-2023). We produced annotated genomic assemblies and performed a pan-genome analysis to identify evidence of virulence gene homology and virulence gene acquisition by linear inheritance or horizontal gene transfer. At least one anthrax-associated virulence gene was annotated in ten isolates. Most homologous sequences in archival isolates showed evidence of pseudogenization and subsequent gene loss. The presence or absence of accessory genes, including anthrax-associated virulence genes, aligned with the phylogenetic structure of the BCG core genome. These findings support the hypothesis that anthrax-associated virulence genes were inherited from a common ancestor in the BCG and were retained or lost across different lineages, and contribute to a growing body of work informing public health strategies related to anthrax surveillance and identification.

Whole genome sequencing and characterization of Corynebacterium isolated from the healthy and dry eye ocular surface

BACKGROUND

The purpose of this study was to characterize Corynebacterium isolated from the ocular surface of dry eye disease patients and healthy controls. We aimed to investigate the pathogenic potential of these isolates in relation to ocular surface health. To this end, we performed whole genome sequencing in combination with biochemical, enzymatic, and antibiotic susceptibility tests. In addition, we employed deferred growth inhibition assays to examine how Corynebacterium isolates may impact the growth of potentially competing microorganisms including the ocular pathogens Pseudomonas aeruginosa and Staphylococcus aureus, as well as other Corynebacterium present on the eye.

RESULTS

The 23 isolates were found to belong to 8 different species of Corynebacterium with genomes ranging from 2.12 mega base pairs in a novel Corynebacterium sp. to 2.65 mega base pairs in C. bovis. Whole genome sequencing revealed the presence of a range of antimicrobial targets present in all isolates. Pangenome analysis showed the presence of 516 core genes and that the pangenome is open. Phenotypic characterization showed variously urease, lipase, mucinase, protease and DNase activity in some isolates. Attention was particularly drawn to a potentially new or novel Corynebacterium species which had the smallest genome, and which produced a range of hydrolytic enzymes. Strikingly the isolate inhibited in vitro the growth of a range of possible pathogenic bacteria as well as other Corynebacterium isolates. The majority of Corynebacterium species included in this study did not seem to possess canonical pathogenic activity.

CONCLUSIONS

This study is the first reported genomic and biochemical characterization of ocular Corynebacterium. A number of potential virulence factors were identified which may have direct relevance for ocular health and contribute to the finding of our previous report on the ocular microbiome, where it was shown that DNA libraries were often dominated by members of this genus. Particularly interesting in this regard was the observation that some Corynebacterium, particularly new or novel Corynebacterium sp. can inhibit the growth of other ocular Corynebacterium as well as known pathogens of the eye.

"Phylogenomic insights into brucellaceae: The Pseudochrobactrum algeriensis case"

The genus Pseudochrobactrum encompasses free-living bacteria phylogenetically close to Ochrobactrum opportunistic pathogens and to Brucella, facultative intracellular parasites causing brucellosis, a worldwide-extended and grave zoonosis. Recently, Pseudochrobactrum strains were isolated from Brucella natural hosts on Brucella selective media, potentially causing diagnostic confusions. Strikingly, P. algeriensis was isolated from cattle lymph nodes, organs that are inimical to bacteria. Here, we analyse P. algeriensis potential virulence factors in comparison with Ochrobactrum and Brucella. Consistent with genomic analyses, Western-Blot analyses confirmed that P. algeriensis lacks the ability to synthesize the N-formylperosamine O-polysaccharide characteristic of the lipopolysaccharide (LPS) of smooth Brucella core species. However, unlike other Pseudochrobactrum but similar to some early diverging brucellae, P. algeriensis carries genes potentially synthetizing a rhamnose-based O-polysaccharide LPS. Lipid A analysis by MALDI-TOF demonstrated that P. algeriensis LPS bears a lipid A with a reduced pathogen-associated molecular pattern, a trait shared with Ochrobactrum and Brucella that is essential to generate a highly stable outer membrane and to delay immune activation. Also, although not able to multiply intracellularly in macrophages, the analysis of P. algeriensis cell lipid envelope revealed the presence of large amounts of cationic aminolipids, which may account for the extremely high resistance of P. algeriensis to bactericidal peptides and could favor colonization of mucosae and transient survival in Brucella hosts. However, two traits critical in Brucella pathogenicity are either significantly different (T4SS [VirB]) or absent (erythritol catabolic pathway) in P. algeriensis. This work shows that, while diverging in other characteristics, lipidic envelope features relevant in Brucella pathogenicity are conserved in Brucellaceae. The constant presence of these features strongly suggests that reinforcement of the envelope integrity as an adaptive advantage in soil was maintained in Brucella because of the similarity of some environmental challenges, such as the action of cationic peptide antibiotics and host defense peptides. This information adds knowledge about the evolution of Brucellaceae, and also underlines the taxonomical differences of the three genera compared.

RIBAP: a comprehensive bacterial core genome annotation pipeline for pangenome calculation beyond the species level

Microbial pangenome analysis identifies present or absent genes in prokaryotic genomes. However, current tools are limited when analyzing species with higher sequence diversity or higher taxonomic orders such as genera or families. The Roary ILP Bacterial core Annotation Pipeline (RIBAP) uses an integer linear programming approach to refine gene clusters predicted by Roary for identifying core genes. RIBAP successfully handles the complexity and diversity of Chlamydia, Klebsiella, Brucella, and Enterococcus genomes, outperforming other established and recent pangenome tools for identifying all-encompassing core genes at the genus level. RIBAP is a freely available Nextflow pipeline at github.com/hoelzer-lab/ribap and zenodo.org/doi/10.5281/zenodo.10890871.

Escherichia coli non-coding regulatory regions are highly conserved

Microbial genome sequences are rapidly accumulating, enabling large-scale studies of sequence variation. Existing studies primarily focus on coding regions to study amino acid substitution patterns in proteins. However, non-coding regulatory regions also play a distinct role in determining physiologic responses. To investigate intergenic sequence variation on a large-scale, we identified non-coding regulatory region alleles across 2350 Escherichia coli strains. This 'alleleome' consists of 117 781 unique alleles for 1169 reference regulatory regions (transcribing 1975 genes) at single base-pair resolution. We find that 64% of nucleotide positions are invariant, and variant positions vary in a median of just 0.6% of strains. Additionally, non-coding alleles are sufficient to recover E. coli phylogroups. We find that core promoter elements and transcription factor binding sites are significantly conserved, especially those located upstream of essential or highly-expressed genes. However, variability in conservation of transcription factor binding sites is significant both within and across regulons. Finally, we contrast mutations acquired during adaptive laboratory evolution with wild-type variation, finding that the former preferentially alter positions that the latter conserves. Overall, this analysis elucidates the wealth of information found in E. coli non-coding sequence variation and expands pangenomic studies to non-coding regulatory regions at single-nucleotide resolution.

Genomic, Phylogenetic and Physiological Characterization of the PAH-Degrading Strain Gordonia polyisoprenivorans 135

The strain Gordonia polyisoprenivorans 135 is able to utilize a wide range of aromatic compounds. The aim of this work was to study the features of genetic organization and biotechnological potential of the strain G. polyisoprenivorans 135 as a degrader of aromatic compounds. The study of the genome of the strain 135 and the pangenome of the G. polyisoprenivorans species revealed that some genes, presumably involved in PAH catabolism, are atypical for Gordonia and belong to the pangenome of Actinobacteria. Analyzing the intergenic regions of strain 135 alongside the "panIGRome" of G. polyisoprenivorans showed that some intergenic regions in strain 135 also differ from those located between the same pairs of genes in related strains. The strain G. polyisoprenivorans 135 in our work utilized naphthalene (degradation degree 39.43%) and grew actively on salicylate. At present, this is the only known strain of G. polyisoprenivorans with experimentally confirmed ability to utilize these compounds.

Analysis of the limited M. tuberculosis accessory genome reveals potential pitfalls of pan-genome analysis approaches

Pan-genome analysis is a fundamental tool for studying bacterial genome evolution; however, the variety of methods used to define and measure the pan-genome poses challenges to the interpretation and reliability of results. To quantify sources of bias and error related to common pan-genome analysis approaches, we evaluated different approaches applied to curated collection of 151 Mycobacterium tuberculosis ( Mtb ) isolates. Mtb is characterized by its clonal evolution, absence of horizontal gene transfer, and limited accessory genome, making it an ideal test case for this study. Using a state-of-the-art graph-genome approach, we found that a majority of the structural variation observed in Mtb originates from rearrangement, deletion, and duplication of redundant nucleotide sequences. In contrast, we found that pan-genome analyses that focus on comparison of coding sequences (at the amino acid level) can yield surprisingly variable results, driven by differences in assembly quality and the softwares used. Upon closer inspection, we found that coding sequence annotation discrepancies were a major contributor to inflated Mtb accessory genome estimates. To address this, we developed panqc, a software that detects annotation discrepancies and collapses nucleotide redundancy in pan-genome estimates. When applied to Mtb and E. coli pan-genomes, panqc exposed distinct biases influenced by the genomic diversity of the population studied. Our findings underscore the need for careful methodological selection and quality control to accurately map the evolutionary dynamics of a bacterial species.

Putative genome contamination has minimal impact on the GTDB taxonomy

The Genome Taxonomy Database (GTDB) provides a species to domain classification of publicly available genomes based on average nucleotide identity (ANI) (for species) and a concatenated gene phylogeny normalized by evolutionary rates (for genus to phylum), which has been widely adopted by the scientific community. Here, we use the Genome UNClutterer (GUNC) software to identify putatively contaminated genomes in GTDB release 07-RS207. We found that GUNC reported 35,723 genomes as putatively contaminated, comprising 11.25 % of the 317,542 genomes in GTDB release 07-RS207. To assess the impact of this high level of inferred contamination on the delineation of taxa, we created 'clean' versions of the 34,846 putatively contaminated bacterial genomes by removing the most contaminated half. For each clean half, we re-calculated the ANI and concatenated gene phylogeny and found that only 77 (0.22 %) of the genomes were not consistent with their original classification. We conclude that the delineation of taxa in GTDB is robust to the putative contamination detected by GUNC.

Genetic barriers more than environmental associations explain Serratia marcescens population structure

Bacterial species often comprise well-separated lineages, likely emerged and maintained by genetic isolation and/or ecological divergence. How these two evolutionary actors interact in the shaping of bacterial population structure is currently not fully understood. In this study, we investigate the genetic and ecological drivers underlying the evolution of Serratia marcescens, an opportunistic pathogen with high genomic flexibility and able to colonise diverse environments. Comparative genomic analyses reveal a population structure composed of five deeply-demarcated genetic clusters with open pan-genome but limited inter-cluster gene flow, partially explained by Restriction-Modification (R-M) systems incompatibility. Furthermore, a large-scale research on hundred-thousands metagenomic datasets reveals only a partial habitat separation of the clusters. Globally, two clusters only show a separate gene composition coherent with ecological adaptations. These results suggest that genetic isolation has preceded ecological adaptations in the shaping of the species diversity, an evolutionary scenario coherent with the Evolutionary Extended Synthesis.

Analysis of twelve genomes of the bacterium Kerstersia gyiorum from brown-throated sloths (Bradypus variegatus), the first from a non-human host

Kerstersia gyiorum is a Gram-negative bacterium found in various animals, including humans, where it has been associated with various infections. Knowledge of the basic biology of K. gyiorum is essential to understand the evolutionary strategies of niche adaptation and how this organism contributes to infectious diseases; however, genomic data about K. gyiorum is very limited, especially from non-human hosts. In this work, we sequenced 12 K. gyiorum genomes isolated from healthy free-living brown-throated sloths (Bradypus variegatus) in the Parque Estadual das Fontes do Ipiranga (São Paulo, Brazil), and compared them with genomes from isolates of human origin, in order to gain insights into genomic diversity, phylogeny, and host specialization of this species. Phylogenetic analysis revealed that these K. gyiorum strains are structured according to host. Despite the fact that sloth isolates were sampled from a single geographic location, the intra-sloth K. gyiorum diversity was divided into three clusters, with differences of more than 1,000 single nucleotide polymorphisms between them, suggesting the circulation of various K. gyiorum lineages in sloths. Genes involved in mobilome and defense mechanisms against mobile genetic elements were the main source of gene content variation between isolates from different hosts. Sloth-specific K. gyiorum genome features include an IncN2 plasmid, a phage sequence, and a CRISPR-Cas system. The broad diversity of defense elements in K. gyiorum (14 systems) may prevent further mobile element flow and explain the low amount of mobile genetic elements in K. gyiorum genomes. Gene content variation may be important for the adaptation of K. gyiorum to different host niches. This study furthers our understanding of diversity, host adaptation, and evolution of K. gyiorum, by presenting and analyzing the first genomes of non-human isolates.

Sequence-subtype association of multi-drug-resistant diarrheagenic Escherichia coli

Background and Objectives

Multi-drug-resistant pathogens pose a significant threat as they can rapidly spread, leading to severe healthcare-associated invasive infections. In developing countries, diarrheagenic Escherichia coli (DEC) is a major bacterial pathogen responsible for causing diarrhea. However, the outbreak of resistant strains has made the treatment of DEC infections much more challenging. This study aimed to investigate the relationship between antibiotic resistance genes and other virulence categories in E. coli strains that cause diarrhea, particularly DEC.

Materials and Methods

The phylogenetic grouping was defined using PCR and multi-locus sequence type (MLST) methods.

Results

Among the isolates analyzed, 14 were identified as resistant and were classified into eight distinct sequence types: ST3, ST53, ST77, ST483, ST512, ST636, ST833, and ST774, indicating genetic diversity among the resistant strains. Certain sequence types, notably ST512 and ST636, were found to be associated with multiple antibiotic resistance in DEC. Regarding antibiotic susceptibility, strains showed the highest resistance to amoxicillin, suggesting that this antibiotic may not be effective in treating DEC infections. On the other hand, the isolates demonstrated susceptibility to amikacin and chloramphenicol, implying that these antibiotics could be more suitable treatment options for DEC infections.

Conclusion

The findings underscore the importance of promptly identifying antibiotic resistance patterns and their correlation with specific pathogenic virulence categories, as this knowledge can aid in selecting the most appropriate antibiotics for treating DEC infections. Considering the antibiotic resistance profiles and associated resistance genes is crucial in managing and containing diarrheal outbreaks and in selecting effective antibiotic therapies for DEC infections.

Pan-Genome Plasticity and Virulence Factors: A Natural Treasure Trove for Acinetobacter baumannii

Acinetobacter baumannii is a Gram-negative pathogen responsible for a variety of community- and hospital-acquired infections. It is recognized as a life-threatening pathogen among hospitalized individuals and, in particular, immunocompromised patients in many countries. A. baumannii, as a member of the ESKAPE group, encompasses high genomic plasticity and simultaneously is predisposed to receive and exchange the mobile genetic elements (MGEs) through horizontal genetic transfer (HGT). Indeed, A. baumannii is a treasure trove that contains a high number of virulence factors. In accordance with these unique pathogenic characteristics of A. baumannii, the authors aim to discuss the natural treasure trove of pan-genome and virulence factors pertaining to this bacterial monster and try to highlight the reasons why this bacterium is a great concern in the global public health system.

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.

Advancement of Autogenous Vaccines in the Poultry Industry

Autogenous vaccines, also known as "custom" vaccines, have become an essential instrument in the production veterinarian's toolbox for the control of emerging and evolving diseases. Autogenous vaccines require a reduced burden of U.S. Department of Agriculture licensing, making them rapidly accessible. Autogenous vaccines have made significant advancements in the ability to reduce disease within the poultry industry from a combination of several different advancements in regulation requirements, rapid and accurate diagnostic assessments, and improvements in manufacturing. The use of autogenous vaccines by poultry health professionals has also increased, and these custom-made products have been instrumental in combating diseases resulting from antigenic variants such as salmonellosis, colibacillosis, infectious coryza, infectious bursal disease, inclusion body hepatitis, viral enteritis, and viral arthritis and tenosynovitis.

Ancient Roman bacterium against current issues: strain Aquil_B6, Paenisporosarcina quisquiliarum, or Psychrobacillus psychrodurans?

IMPORTANCE

Since 1988, through the United States government's founding, the National Center for Biotechnology Information (NCBI) has provided an invaluable service to scientific advancement. The universality and total freedom of use if on the one hand allow the use of this database on a global level by all researchers for their valuable work, on the other hand, it has the disadvantage of making it difficult to check the correctness of all the materials present. It is, therefore, of fundamental importance for the correctness and ethics of research to improve the databases at our disposal, identifying and amending the critical issues. This work aims to provide the scientific community with a new sequence for the type strain Paenisporosarcina quisquiliarum SK 55 and broaden the knowledge of the Psychrobacillus psychrodurans species, in particular, considering the ancient strain Aquil_B6 found in an ancient Roman amphora.

PanDelos-frags: A methodology for discovering pangenomic content of incomplete microbial assemblies

Pangenomics was originally defined as the problem of comparing the composition of genes into gene families within a set of bacterial isolates belonging to the same species. The problem requires the calculation of sequence homology among such genes. When combined with metagenomics, namely for human microbiome composition analysis, gene-oriented pangenome detection becomes a promising method to decipher ecosystem functions and population-level evolution. Established computational tools are able to investigate the genetic content of isolates for which a complete genomic sequence is available. However, there is a plethora of incomplete genomes that are available on public resources, which only a few tools may analyze. Incomplete means that the process for reconstructing their genomic sequence is not complete, and only fragments of their sequence are currently available. However, the information contained in these fragments may play an essential role in the analyses. Here, we present PanDelos-frags, a computational tool which exploits and extends previous results in analyzing complete genomes. It provides a new methodology for inferring missing genetic information and thus for managing incomplete genomes. PanDelos-frags outperforms state-of-the-art approaches in reconstructing gene families in synthetic benchmarks and in a real use case of metagenomics. PanDelos-frags is publicly available at https://github.com/InfOmics/PanDelos-frags.

Comparative Genome analysis of the Genus Curvibacter and the Description of Curvibacter microcysteis sp. nov. and Curvibacter cyanobacteriorum sp. nov., Isolated from Fresh Water during the Cyanobacterial Bloom Period

The three Gram-negative, catalase- and oxidase-positive bacterial strains RS43T, HBC28, and HBC61T, were isolated from fresh water and subjected to a polyphasic study. Comparison of 16S rRNA gene sequence initially indicated that strains RS43T, HBC28, and HBC61T were closely related to species of genus Curvibacter and shared the highest sequence similarity of 98.14%, 98.21%, and 98.76%, respectively, with Curvibacter gracilis 7-1T. Phylogenetic analysis based on genome sequences placed all strains within the genus Curvibacter. The average nucleotide identity (ANI) and digital DNA-DNA hybridization (dDDH) values between the three strains and related type strains supported their recognition as two novel genospecies in the genus Curvibacter. Comparative genomic analysis revealed that the genus possessed an open pangenome. Based on KEGG BlastKOALA analyses, Curvibacter species have the potential to metabolize benzoate, phenylacetate, catechol, and salicylate, indicating their potential use in the elimination of these compounds from the water systems. The results of polyphasic characterization indicated that strain RS43T and HBC61T represent two novel species, for which the name Curvibacter microcysteis sp. nov. (type strain RS43T =KCTC 92793T=LMG 32714T) and Curvibacter cyanobacteriorum sp. nov. (type strain HBC61T =KCTC 92794T =LMG 32713T) are proposed.

Identification of novel Carnobacterium maltaromaticum strains in bone marrow samples of patients with acute myeloid leukemia using a metagenomic binning approach

The aim of this study aimed to examine the existence of a bacterial metagenome in the bone marrow of patients with acute myeloid leukemia (AML). We re-examined whole-genome sequencing data from the bone marrow samples of seven patients with AML, four of whom were remitted after treatment, for metagenomic analysis. After the removal of human reads, unmapped reads were used to profile the species-level composition. We used the metagenomic binning approach to confirm whether the identified taxon was a complete genome of known or novel strains. We observed a unique and novel microbial signature in which Carnobacterium maltaromaticum was the most abundant species in five patients with AML or remission. The complete genome of C. maltaromaticum "BMAML_KR01," which was observed in all samples, was 100% complete with 8.5% contamination and closely clustered with C. maltaromaticum strains DSM20730 and SF668 based on single nucleotide polymorphism variations. We identified five unique proteins that could contribute to cancer progression and 104 virulent factor proteins in the BMAML_KR01 genome. To our knowledge, this is the first report of a new strain of C. maltaromaticum in patients with AML. The presence of C. maltaromaticum and its new strain in patients indicates an urgent need to validate the existence of this bacterium and evaluate its pathophysiological role.

One Earth: The Equilibrium between the Human and the Bacterial Worlds

Misuse and abuse of antibiotics on humans, cattle, and crops have led to the selection of multi-resistant pathogenic bacteria, the most feared 'superbugs'. Infections caused by superbugs are progressively difficult to treat, with a subsequent increase in lethality: the toll on human lives is predicted to reach 10 million by 2050. Here we review three concepts linked to the growing resistance to antibiotics, namely (i) the Resistome, which refers to the collection of bacterial genes that confer resistance to antibiotics, (ii) the Mobilome, which includes all the mobile genetic elements that participate in the spreading of antibiotic resistance among bacteria by horizontal gene transfer processes, and (iii) the Nichome, which refers to the set of genes that are expressed when bacteria try to colonize new niches. We also discuss the strategies that can be used to tackle bacterial infections and propose an entente cordiale with the bacterial world so that instead of war and destruction of the 'fierce enemy' we can achieve a peaceful coexistence (the One Earth concept) between the human and the bacterial worlds. This, in turn, will contribute to microbial biodiversity, which is crucial in a globally changing climate due to anthropogenic activities.

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.

happi: a hierarchical approach to pangenomics inference

Recovering metagenome-assembled genomes (MAGs) from shotgun sequencing data is an increasingly common task in microbiome studies, as MAGs provide deeper insight into the functional potential of both culturable and non-culturable microorganisms. However, metagenome-assembled genomes vary in quality and may contain omissions and contamination. These errors present challenges for detecting genes and comparing gene enrichment across sample types. To address this, we propose happi, an approach to testing hypotheses about gene enrichment that accounts for genome quality. We illustrate the advantages of happi over existing approaches using published Saccharibacteria MAGs, Streptococcus thermophilus MAGs, and via simulation.

Pangenome analysis of Corynebacterium striatum: insights into a neglected multidrug-resistant pathogen

BACKGROUND

Over the past two decades, Corynebacterium striatum has been increasingly isolated from clinical cultures with most isolates showing increased antimicrobial resistance (AMR) to last resort agents. Advances in the field of pan genomics would facilitate the understanding of the clinical significance of such bacterial species previously thought to be among commensals paving the way for identifying new drug targets and control strategies.

METHODS

We constructed a pan-genome using 310 genome sequences of C. striatum. Pan-genome analysis was performed using three tools including Roary, PIRATE, and PEPPAN. AMR genes and virulence factors have been studied in relation to core genome phylogeny. Genomic Islands (GIs), Integrons, and Prophage regions have been explored in detail.

RESULTS

The pan-genome ranges between a total of 5253-5857 genes with 2070 - 1899 core gene clusters. Some antimicrobial resistance genes have been identified in the core genome portion, but most of them were located in the dispensable genome. In addition, some well-known virulence factors described in pathogenic Corynebacterium species were located in the dispensable genome. A total of 115 phage species have been identified with only 44 intact prophage regions.

CONCLUSION

This study presents a detailed comparative pangenome report of C. striatum. The species show a very slowly growing pangenome with relatively high number of genes in the core genome contributing to lower genomic variation. Prophage elements carrying AMR and virulence elements appear to be infrequent in the species. GIs appear to offer a prominent role in mobilizing antibiotic resistance genes in the species and integrons occur at a frequency of 50% in the species. Control strategies should be directed against virulence and resistance determinants carried on the core genome and those frequently occurring in the accessory genome.

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.

Phylogenomic analysis of the Porphyromonas gingivalis - Porphyromonas gulae duo: approaches to the origin of periodontitis

Porphyromonas gingivalis is an oral human pathogen associated with the onset and progression of periodontitis, a chronic immune-inflammatory disease characterized by the destruction of the teeth-supporting tissue. P. gingivalis belongs to the genus Porphyromonas, which is characterized by being composed of Gram-negative, asaccharolytic, non-spore-forming, non-motile, obligatory anaerobic species, inhabiting niches such as the oral cavity, urogenital tract, gastrointestinal tract and infected wound from different mammals including humans. Among the Porphyromonas genus, P. gingivalis stands out for its specificity in colonizing the human oral cavity and its keystone pathogen role in periodontitis pathogenesis. To understand the evolutionary process behind P. gingivalis in the context of the Pophyoromonas genus, in this study, we performed a comparative genomics study with publicly available Porphyromonas genomes, focused on four main objectives: (A) to confirm the phylogenetic position of P. gingivalis in the Porphyromonas genus by phylogenomic analysis; (B) the definition and comparison of the pangenomes of P. gingivalis and its relative P. gulae; and (C) the evaluation of the gene family gain/loss events during the divergence of P. gingivalis and P. gulae; (D) the evaluation of the evolutionary pressure (represented by the calculation of Tajima-D values and dN/dS ratios) comparing gene families of P. gingivalis and P. gulae. Our analysis found 84 high-quality assemblies representing P. gingivalis and 14 P. gulae strains (from a total of 233 Porphyromonas genomes). Phylogenomic analysis confirmed that P. gingivalis and P. gulae are highly related lineages, close to P. loveana. Both organisms harbored open pangenomes, with a strong core-to-accessory ratio for housekeeping genes and a negative ratio for unknown function genes. Our analyses also characterized the gene set differentiating P. gulae from P. gingivalis, mainly associated with unknown functions. Relevant virulence factors, such as the FimA, Mfa1, and the hemagglutinins, are conserved in P. gulae, P. gingivalis, and P. loveana, suggesting that the origin of those factors occurred previous to the P. gulae - P. gingivalis divergence. These results suggest an unexpected evolutionary relationship between the P. gulae - P. gingivalis duo and P. loveana, showing more clues about the origin of the role of those organisms in periodontitis.

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

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

Advanced Methods for Natural Products Discovery: Bioactivity Screening, Dereplication, Metabolomics Profiling, Genomic Sequencing, Databases and Informatic Tools, and Structure Elucidation

Natural Products (NP) are essential for the discovery of novel drugs and products for numerous biotechnological applications. The NP discovery process is expensive and time-consuming, having as major hurdles dereplication (early identification of known compounds) and structure elucidation, particularly the determination of the absolute configuration of metabolites with stereogenic centers. This review comprehensively focuses on recent technological and instrumental advances, highlighting the development of methods that alleviate these obstacles, paving the way for accelerating NP discovery towards biotechnological applications. Herein, we emphasize the most innovative high-throughput tools and methods for advancing bioactivity screening, NP chemical analysis, dereplication, metabolite profiling, metabolomics, genome sequencing and/or genomics approaches, databases, bioinformatics, chemoinformatics, and three-dimensional NP structure elucidation.

Comparative genomics reveals the acquisition of mobile genetic elements by the plant growth-promoting Pantoea eucrina OB49 in polluted environments

Heavy metal-tolerant plant growth-promoting bacteria (PGPB) have gained popularity in bioremediation in recent years. A genome-assisted study of a heavy metal-tolerant PGPB Pantoea eucrina OB49 isolated from the rhizosphere of wheat grown on a heavy metal-contaminated site is presented. Comparative pan-genome analysis indicated that OB49 acquired heavy metal resistance genes through horizontal gene transfer. On contigs S10 and S12, OB49 has two arsRBCH operons that give arsenic resistance. On the S12 contig, an arsRBCH operon was discovered in conjunction with the merRTPCADE operon, which provides mercury resistance. P. eucrina OB49 may be involved in an ecological alternative for heavy metal remediation and growth promotion of wheat grown in metal-polluted soils. Our results suggested the detection of mobile genetic elements that harbour the ars operon and the fluoride resistance genes adjacent to the mer operon.

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.

Pan-genome association study of Mycobacterium tuberculosis lineage-4 revealed specific genes related to the high and low prevalence of the disease in patients from the North-Eastern area of Medellín, Colombia

Mycobacterium tuberculosis (Mtb) lineage 4 is responsible for the highest burden of tuberculosis (TB) worldwide. This lineage has been the most prevalent lineage in Colombia, especially in the North-Eastern (NE) area of Medellin, where it has been shown to have a high prevalence of LAM9 SIT42 and Haarlem1 SIT62 sublineages. There is evidence that regardless of environmental factors and host genetics, differences among sublineages of Mtb strains play an important role in the course of infection and disease. Nevertheless, the genetic basis of the success of a sublineage in a specific geographic area remains uncertain. We used a pan-genome-wide association study (pan-GWAS) of 47 Mtb strains isolated from NE Medellin between 2005 and 2008 to identify the genes responsible for the phenotypic differences among high and low prevalence sublineages. Our results allowed the identification of 12 variants in 11 genes, of which 4 genes showed the strongest association to low prevalence (mmpL12, PPE29, Rv1419, and Rv1762c). The first three have been described as necessary for invasion and intracellular survival. Polymorphisms identified in low prevalence isolates may suggest related to a fitness cost of Mtb, which might reflect a decrease in their capacity to be transmitted or to cause an active infection. These results contribute to understanding the success of some sublineages of lineage-4 in a specific geographical area.

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.

Genomic analysis of Paenibacillus sp. MDMC362 from the Merzouga desert leads to the identification of a potentially thermostable catalase

Microorganisms in hot deserts face heat and other environmental conditions, such as desiccation, UV radiation, or low nutrient availability. Therefore, this hostile environment harbour microorganisms with acquired characteristics related to survival in their habitat, which can be exploited in biotechnology. In this work, the genome of Paenibacillus sp. MDMC362 isolated from the Merzouga desert in Morocco was sequenced to understand its survival strategy's genetic basis; and to evaluate the thermostability of a catalase extracted from genomic annotation files using molecular dynamics. Paenibacillus sp. MDMC362 genome was rich in genetic elements involved in the fight against different stresses, notably temperature stress, UV radiations, osmotic stress, carbon starvation, and oxidative stress. Indeed, we could identify genes of the operons groES-groEL and hrcA-grpE-dnaK and those involved in the different stages of sporulation, which can help the bacteria to survive the high temperatures imposed by a desertic environment. We also observed the genetic components of the UvrABC system and additional mechanisms involved in DNA repair, which help overcome UV radiation damage. Other genes have been identified in the genome, like those coding for ectoine and proline, that aids fight osmotic stress and desiccation. Catalase thermostability investigation using molecular dynamics showed that the protein reached stability and conserved its compactness at temperatures up to 373.15 K. These results suggest a potential thermostability of the enzyme. Since the studied protein is a core protein, thermostability could be conserved among Paenibacillus sp. MDMC362 closely related strains; however, bacteria from harsh environments may have a slight advantage regarding protein stability.

Contribution of -Omics Technologies in the Study of Porphyromonas gingivalis during Periodontitis Pathogenesis: A Minireview

Periodontitis is a non-communicable chronic inflammatory disease characterized by the progressive and irreversible breakdown of the soft periodontal tissues and resorption of teeth-supporting alveolar bone. The etiology of periodontitis involves dysbiotic shifts in the diversity of microbial communities inhabiting the subgingival crevice, which is dominated by anaerobic Gram-negative bacteria, including Porphyromonas gingivalis. Indeed, P. gingivalis is a keystone pathogen with a repertoire of attributes that allow it to colonize periodontal tissues and influence the metabolism, growth rate, and virulence of other periodontal bacteria. The pathogenic potential of P. gingivalis has been traditionally analyzed using classical biochemical and molecular approaches. However, the arrival of new techniques, such as whole-genome sequencing, metagenomics, metatranscriptomics, proteomics, and metabolomics, allowed the generation of high-throughput data, offering a suitable option for bacterial analysis, allowing a deeper understanding of the pathogenic properties of P. gingivalis and its interaction with the host. In the present review, we revise the use of the different -omics technologies and techniques used to analyze bacteria and discuss their potential in studying the pathogenic potential of P. gingivalis.

Phylogenomic characterization and pangenomic insights into the surfactin-producing bacteria Bacillus subtilis strain RI4914

Bacillus subtilis is a versatile bacterial species able to produce surfactin, a lipopeptide biosurfactant. We carried out the phylogenomic characterization and pangenomic analyses using available B. subtilis complete genomes. Also, we report the whole genome of the biosurfactant-producing B. subtilis strain RI4914 that was isolated from effluent water from an oil exploration field. We applied a hybrid sequencing approach using both long- and short-read sequencing technologies to generate a highly accurate, single-chromosome genome. The pangenomics analysis of 153 complete genomes classified as B. subtilis retrieved from the NCBI shows an open pangenome composed of 28,511 accessory genes, which agrees with the high genetic plasticity of the species. Also, this analysis suggests that surfactin production is a common trait shared by members of this species since the srfA operon is highly conserved among the B. subtilis strains found in most of the assemblies available. Finally, increased surfactin production corroborates the higher srfAA gene expression in B. subtilis strain RI4914.

Genomics and pathotypes of the many faces of Escherichia coli

Escherichia coli is the most researched microbial organism in the world. Its varied impact on human health, consisting of commensalism, gastrointestinal disease, or extraintestinal pathologies, has generated a separation of the species into at least eleven pathotypes (also known as pathovars). These are broadly split into two groups, intestinal pathogenic E. coli (InPEC) and extraintestinal pathogenic E. coli (ExPEC). However, components of E. coli's infinite open accessory genome are horizontally transferred with substantial frequency, creating pathogenic hybrid strains that defy a clear pathotype designation. Here, we take a birds-eye view of the E. coli species, characterizing it from historical, clinical, and genetic perspectives. We examine the wide spectrum of human disease caused by E. coli, the genome content of the bacterium, and its propensity to acquire, exchange, and maintain antibiotic resistance genes and virulence traits. Our portrayal of the species also discusses elements that have shaped its overall population structure and summarizes the current state of vaccine development targeted at the most frequent E. coli pathovars. In our conclusions, we advocate streamlining efforts for clinical reporting of ExPEC, and emphasize the pathogenic potential that exists throughout the entire species.

Whole-Genome Sequence Comparisons of Listeria monocytogenes Isolated from Meat and Fish Reveal High Inter- and Intra-Sample Diversity

Interpretation of whole-genome sequencing (WGS) data for foodborne outbreak investigations is complex, as the genetic diversity within processing plants and transmission events need to be considered. In this study, we analyzed 92 food-associated Listeria monocytogenes isolates by WGS-based methods. We aimed to examine the genetic diversity within meat and fish production chains and to assess the applicability of suggested thresholds for clustering of potentially related isolates. Therefore, meat-associated isolates originating from the same samples or processing plants as well as fish-associated isolates were analyzed as distinct sets. In silico serogrouping, multilocus sequence typing (MLST), core genome MLST (cgMLST), and pangenome analysis were combined with screenings for prophages and genetic traits. Isolates of the same subtypes (cgMLST types (CTs) or MLST sequence types (STs)) were additionally compared by SNP calling. This revealed the occurrence of more than one CT within all three investigated plants and within two samples. Analysis of the fish set resulted in predominant assignment of isolates from pangasius catfish and salmon to ST2 and ST121, respectively, potentially indicating persistence within the respective production chains. The approach not only allowed the detection of distinct subtypes but also the determination of differences between closely related isolates, which need to be considered when interpreting WGS data for surveillance.

Pan-Genome Analysis Reveals Functional Divergences in Gut-Restricted Gilliamella and Snodgrassella

Gilliamella and Snodgrassella, members of core gut microbiota in corbiculate bees, have high species diversity and adaptability to a wide range of hosts. In this study, we performed species taxonomy and phylogenetic analysis for Gilliamella and Snodgrassella strains that we isolated in our laboratory, in combination with published whole-genome. Functional effects of accessory and unique genes were investigated by KEGG category and pathway annotation in pan-genome analysis. Consequently, in Gilliamella, we inferred the importance of carbohydrate metabolism, amino acid metabolism, membrane transport, energy metabolism, and metabolism of cofactors and vitamins in accessory or unique genes. The pathway mentioned above, plus infectious disease, lipid metabolism, nucleotide metabolism as well as replication and repair exert a pivotal role in accessory or unique genes of Snodgrassella. Further analysis revealed the existence of functional differentiation of accessory and unique genes among Apis-derived genomes and Bombus-derived genomes. We also identified eight and four biosynthetic gene clusters in all Gilliamella and Snodgrassella genomes, respectively. Our study provides a good insight to better understand how host heterogeneity influences the bacterial speciation and affects the versatility of the genome of the gut bacteria.

Comparative genomics reveals the evolution of antimicrobial resistance in Bacteroides nordii

Bacteroides nordii, is an understudied member of the pathogenic B. fragilis group which comprises several multidrug-resistant (MDR) strains. Thus, it is of great interest to study the genome biology of Bacteroides nordii. However, no detailed study is available that characterized B. nordii at the genetic level and explored its role as a potential pathogen. We isolated an MDR strain viz., B. nordii PGMM4098 from the pus sample and subjected it to whole genome sequencing using Illumina technology. The draft genome was de-novo assembled and annotated, followed by comprehensive comparative genomics analyses using the publicly available genome dataset of B. nordii. The pan-genome analysis revealed the open nature of B. nordii, indicating the continuous accumulation of novel genes in non-core components leading to the emergence of new strains of this species. The thirteen antimicrobial resistance (AMR) genes identified in the genomes of all B. nordii strains were part of the non-core component of the pan-genome. Of these, four AMR genes, nimE, aadS, mef(En2), and ermB/F/G were found to be acquired via the process of horizontal gene transfer (HGT) from anaerobic Bacteroidetes. Importantly, the nimE gene conferring metronidazole resistance was found to be present only in B. nordii PGMM4098, which harbors five other AMR genes encoded in its genome. Of these, nimE (metronidazole resistance), ermB/F/G (macrolide-lincosamide-streptogramin B resistance), and cfxA2/A3 (class A β-lactam resistance) genes were further validated using targeted polymerase chain reaction assay. Notably, these three genes were also found to be under the operation of positive selective pressure suggesting the diversification of these genes, which might lead to the emergence of new MDR strains of B. nordii in the near future. Our study reported and characterized the genome of the first MDR strain of B. nordii and revealed the AMR evolution in this species using a comprehensive comparative genomics approach.

ProPan: a comprehensive database for profiling prokaryotic pan-genome dynamics

Compared with conventional comparative genomics, the recent studies in pan-genomics have provided further insights into species genomic dynamics, taxonomy and identification, pathogenicity and environmental adaptation. To better understand genome characteristics of species of interest and to fully excavate key metabolic and resistant genes and their conservations and variations, here we present ProPan (https://ngdc.cncb.ac.cn/propan), a public database covering 23 archaeal species and 1,481 bacterial species (in a total of 51,882 strains) for comprehensively profiling prokaryotic pan-genome dynamics. By analyzing and integrating these massive datasets, ProPan offers three major aspects for the pan-genome dynamics of the species of interest: 1) the evaluations of various species' characteristics and composition in pan-genome dynamics; 2) the visualization of map association, the functional annotation and presence/absence variation for all contained species' gene clusters; 3) the typical characteristics of the environmental adaptation, including resistance genes prediction of 126 substances (biocide, antimicrobial drug and metal) and evaluation of 31 metabolic cycle processes. Besides, ProPan develops a very user-friendly interface, flexible retrieval and multi-level real-time statistical visualization. Taken together, ProPan will serve as a weighty resource for the studies of prokaryotic pan-genome dynamics, taxonomy and identification as well as environmental adaptation.

How GBS Got Its Hump: Genomic Analysis of Group B Streptococcus from Camels Identifies Host Restriction as well as Mobile Genetic Elements Shared across Hosts and Pathogens

Group B Streptococcus (GBS) literature largely focuses on humans and neonatal disease, but GBS also affects numerous animals, with significant impacts on health and productivity. Spill-over events occur between humans and animals and may be followed by amplification and evolutionary adaptation in the new niche, including changes in the core or accessory genome content. Here, we describe GBS from one-humped camels (Camelus dromedarius), a relatively poorly studied GBS host of increasing importance for food security in arid regions. Genomic analysis shows that virtually all GBS from camels in East Africa belong to a monophyletic clade, sublineage (SL)609. Capsular types IV and VI, including a new variant of type IV, were over-represented compared to other host species. Two genomic islands with signatures of mobile elements contained most camel-associated genes, including genes for metal and carbohydrate utilisation. Lactose fermentation genes were associated with milk isolates, albeit at lower prevalence in camel than bovine GBS. The presence of a phage with high identity to Streptococcus pneumoniae and Streptococcus suis suggests lateral gene transfer between GBS and bacterial species that have not been described in camels. The evolution of camel GBS appears to combine host restriction with the sharing of accessory genome content across pathogen and host species.

A genome-wide association study of the occurrence of genetic variations in Edwardsiella piscicida, Vibrio harveyi, and Streptococcus parauberis under stressed environments

Bacterial mutation and genetic diversity in aquaculture have led to increasing phenotypic variances, which can weaken or invalidate strategies for controlling diseases. However, few studies have monitored the degree of mutation in fish bacterial pathogens caused by environmental pressure within a short period. In this study, transcriptomic sequences from Edwardsiella piscicida, Vibrio harveyi and Streptococcus parauberis under stressed environments were used for investigating the emergence of variants. In detail, a sub-inhibitory concentration of formalin and phenol for E. piscicida, sea water at 30°C for V. harveyi and flounder serum for S. parauberis were used as stressed environments, and significant single-nucleotide polymorphisms (SNPs) and/or mutation sites were investigated after culture in the ordinary liquid media (control) and the stressed environment through a genome-wide association study. As results, several SNPs or mutations during incubation were observed under different environments in E. piscicida and/or V. harveyi in the genes relevant to flagella, fimbria type 3 secretion systems, and outer and inner membranes that have been directly exposed to external environments. In particular, given that flagella and fimbriae are considered important factors in differentiating the serotypes in some bacterial pathogens, it can be speculated that different environmental pressures are the source of phenotypic or serotypic differentiation from the same origin. On the other hands, S. parauberis did not exhibit notable changes for 4 h when inoculated in the serum from olive flounder. The results presented in this study provide examples of possible molecular evolution in pathogens relevant to the aquaculture industry as a response to different environmental pressure.

Analysis of the Taxonomy, Synteny, and Virulence Factors for Soft Rot Pathogen Pectobacterium aroidearum in Amorphophallus konjac Using Comparative Genomics

Bacterial soft rot is a devastating disease for a wide range of crops, vegetables, and ornamental plants including konjac (Amorphophallus konjac). However, the pangenome and genomic plasticity of the konjac soft rot pathogens is little explored. In this study, we reported the complete genome sequences of 11 bacterial isolates that can cause typical soft rot symptoms in konjac by in vitro and in vivo pathogenicity tests. Based on in silico DNA-DNA hybridization, average nucleotide identity and phylogenomic analysis, all 11 isolates were determined to be Pectobacterium aroidearum. In addition, synteny analysis of these genomes revealed considerable chromosomal inversions, one of which is triggered by homologous recombination of ribose operon. Pangenome analysis and COG enrichment analysis showed that the pangenome of P. aroidearum is open and that accessory genes are enriched in replication, recombination, and repair. Variations in type IV secretion system and type VI secretion system were found, while plant cell wall degrading enzymes were conserved. Furthermore, sequence analyses also provided evidence for the presence of a type V secretion system in Pectobacterium. These findings advance our understanding of the pathogenicity determinants, genomic plasticity, and evolution of P. aroidearum.

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.

Comparative Genomics Provides Insights Into Genetic Diversity of Clostridium tyrobutyricum and Potential Implications for Late Blowing Defects in Cheese

Clostridium tyrobutyricum has been recognized as the main cause of late blowing defects (LBD) in cheese leading to considerable economic losses for the dairy industry. Although differences in spoilage ability among strains of this species have been acknowledged, potential links to the genetic diversity and functional traits remain unknown. In the present study, we aimed to investigate and characterize genomic variation, pan-genomic diversity and key traits of C. tyrobutyricum by comparing the genomes of 28 strains. A comparative genomics analysis revealed an “open” pangenome comprising 9,748 genes and a core genome of 1,179 genes shared by all test strains. Among those core genes, the majority of genes encode proteins related to translation, ribosomal structure and biogenesis, energy production and conversion, and amino acid metabolism. A large part of the accessory genome is composed of sets of unique, strain-specific genes ranging from about 5 to more than 980 genes. Furthermore, functional analysis revealed several strain-specific genes related to replication, recombination and repair, cell wall, membrane and envelope biogenesis, and defense mechanisms that might facilitate survival under stressful environmental conditions. Phylogenomic analysis divided strains into two clades: clade I contained human, mud, and silage isolates, whereas clade II comprised cheese and milk isolates. Notably, these two groups of isolates showed differences in certain hypothetical proteins, transcriptional regulators and ABC transporters involved in resistance to oxidative stress. To the best of our knowledge, this is the first study to provide comparative genomics of C. tyrobutyricum strains related to LBD. Importantly, the findings presented in this study highlight the broad genetic diversity of C. tyrobutyricum, which might help us understand the diversity in spoilage potential of C. tyrobutyricum in cheese and provide some clues for further exploring the gene modules responsible for the spoilage ability of this species.

Sequence-based pangenomic core detection

One of the most basic kinds of analysis to be performed on a pangenome is the detection of its core, i.e., the information shared among all members. Pangenomic core detection is classically done on the gene level and many tools focus exclusively on core detection in prokaryotes. Here, we present a new method for sequence-based pangenomic core detection. Our model generalizes from a strict core definition allowing us to flexibly determine suitable core properties depending on the research question and the dataset under consideration. We propose an algorithm based on a colored de Bruijn graph that runs in linear time with respect to the number of k-mers in the graph. An implementation of our method is called Corer. Because of the usage of a colored de Bruijn graph, it works alignment-free, is provided with a small memory footprint, and accepts as input assembled genomes as well as sequencing reads.

•

Pangenomic core detection for large collections of prokaryotes or higher eukaryotes

•

Whole-genome analysis with assemblies or even read data as input

•

Alignment-free, linear time algorithm with small memory footprint

•

Variation tolerance and quorum for flexible core detection

Genome analysis uncovers the prolific antagonistic and plant growth-promoting potential of endophyte Bacillus velezensis K1

Insights into the application of endophytic bacilli in sustainable agricultural practices have opened up new avenues for the inhibition of soil-borne pathogens and the improvement of plant health. Bacillus subtilis K1, an endophytic bacterium originally isolated from aerial roots of Ficus benghalensis is a potential biocontrol agent secreting a mixture of surfactins, iturins and fengycins. The current study extends the characterization of this bacterium through genomic and comparative genomics approaches. The sequencing of the bacterial genome at Illumina MiSeq platform revealed that it possessed a 4,103,502-bp circular chromosome with 45.98% GC content and 4325 predicted protein-coding sequences. Based on phylogenomics and whole-genome average nucleotide identity, the B. subtilis K1 was taxonomically classified as Bacillus velezensis. The formerly evaluated phenotypic traits viz. C-source utilization and lipopeptide-mediated fungal antagonism were correlated to their molecular determinants. The genome also harbored several genes associated with induced systemic resistance and plant growth promotion i.e, phytohormone production, nitrogen assimilation and reduction, siderophore production, phosphate solubilization, biofilm formation, swarming motility, acetoin and butanediol synthesis. The production of antifungal volatile organic compounds and plant growth promotion was experimentally demonstrated by volatile compound assay and seed germination assay on cumin and groundnut. The isolate also holds great prospects for application as a soil inoculant as indicated by enhancement in the growth of groundnut via in planta pot studies. Bacterial pan-genome analysis based on a comparison of whole genomes with eighteen other Bacillus strains was also conducted. Comparative examination of biosynthetic gene clusters across all genomes indicated that the largest number of gene clusters were harbored by the K1 genome. Based on the findings, we propose K1 as a model for scrutinizing non-ribosomally synthesized peptide synthetase and polyketide synthetase derived molecules.

GAAP: A GUI-based Genome Assembly and Annotation Package

Background: Next-generation sequencing (NGS) technologies are being continuously used for high-throughput sequencing data generation that requires easy-to-use GUI-based data analysis software. These kinds of software could be used in-parallel with sequencing for the automatic data analysis. At present, very few software are available for use and most of them are commercial, thus creating a gap between data generation and data analysis. Methods: GAAP is developed on the NodeJS platform that uses HTML, JavaScript as the front-end for communication with users. We have implemented FastQC and trimmomatic tool for quality checking and control. Velvet and Prodigal are integrated for genome assembly and gene prediction. The annotation will be done with the help of remote NCBI Blast and IPR-Scan. In the back- end, we have used PERL and JavaScript for the processing of data. To evaluate the performance of GAAP, we have assembled a viral (SRR11621811), bacterial (SRR17153353) and human genome (SRR16845439). Results: We have used GAAP software to assemble, and annotate a COVID-19 genome on a desktop computer that resulted in a single contig of 27994bp with 99.57% reference genome coverage. This assembly predicted 11 genes, of which 10 were annotated using annotation module of GAAP. We have also assembled a bacterial and human genome 138 and 194281 contigs with N50 value 100399 and 610, respectively. Conclusion: In this study, we have developed freely available, platform-independent genome assembly and annotation (GAAP) software (www.deepaklab.com/gaap). The software itself acts as a complete data analysis package with quality check, quality control, de-novo genome assembly, gene prediction and annotation (Blast, PFAM, GO-Term, pathway and enzyme mapping) modules.

Comparative genome analysis of Pasteurella multocida serogroup B:2 strains causing haemorrhagic septicaemia (HS) in bovines

Pasteurella multocida, a Gram-negative bacterium with ubiquitous nature, is known to affect wide range of host species worldwide with varied clinical manifestations including haemorrhagic septicaemia (HS) in bovines. Although, HS causing P. multocida strains were identified and characterized by conventional tools and PCR assays, diverse strains are indistinguishable by these tools in the face of disease outbreaks. In this study, draft genomes of three virulent P. multocida serotype B:2 strains (NIVEDIPm32, NIVEDIPm34 and NIVEDIPm35) were analyzed following whole genome sequencing, assembly, annotation and compared them with existing global genomes (n = 43) of bovine origin in the database. Three draft genomes of NIVEDIPm strains consisted of 40-52 contigs with GC content of ∼40.4%. The genome size and predicted genes content was ∼2.3 Mb and 2181-2189, respectively. Besides, the presence of various mobile genetic elements, antimicrobial resistance genes and biofilm related genes suggested their vital roles in virulence; further, adaptation to the host immune system as well as host pathogen interaction. Multi locus sequence analysis based on RIRDC scheme showed the presence of ST122 in all the three strains. wgMLST based phylogenic analysis suggested that HS causing Indian virulent field strains differed geographically and showed diversity from existing HS vaccine strain P52. The phylogenetic tree revealed that North Indian strains share high similarity with strains of Pakistan than South Indian Strain. Notably, a high divergence of SNPs between the HS causing circulating virulent strains of India and current HS vaccine strain P52 suggested an imminent need for relook in to HS vaccination strategy for livestock in India.

Pangenome analysis of Enterobacteria reveals richness of secondary metabolite gene clusters and their associated gene sets

In silico genome mining provides easy access to secondary metabolite biosynthetic gene clusters (BGCs) encoding the biosynthesis of many bioactive compounds, which are the basis for many important drugs used in human medicine. However, the association between BGCs and other functions encoded in the genomes of producers have remained elusive. Here, we present a systems biology workflow that integrates genome mining with a detailed pangenome analysis for detecting genes associated with a particular BGC. We analyzed 3,889 enterobacterial genomes and found 13,266 BGCs, represented by 252 distinct BGC families and 347 additional singletons. A pangenome analysis revealed 88 genes putatively associated with a specific BGC coding for the colon cancer-related colibactin that code for diverse metabolic and regulatory functions. The presented workflow opens up the possibility to discover novel secondary metabolites, better understand their physiological roles, and provides a guide to identify and analyze BGC associated gene sets.