Inside the Pan-genome - Methods and Software Overview

Chromosomal gene order defines several structural classes of Staphylococcus epidermidis genomes

The original methodology for describing the pangenome of a prokaryotic species is based on modeling genomes as unordered sets of genes. More recent findings have underlined the importance of considering the ordering of genes along the genetic material as well, when making comparisons among genomes. To further investigate the benefits of gene order when describing genomes of a given species, we applied two distance metrics on a dataset of 84 genomes of Staphylococcus epidermidis. The first metric, GeLev, depends on the order of genes and is a derivative of the Levenshtein distance. The second, the Jaccard distance, depends on gene sets only. The application of these distances reveals information about the global structure of the genomes, and allows clustering of the genomes into classes. The main biological result is that, while genomes within the same class are structurally similar, genomes of different classes have an additional characteristic. Between genomes in different classes we can discover instances where a large segment of the first genome appears in reverse order in the second. This feature suggests that genome rearrangements in S. epidermidis happen on a large scale, while micro-rearrangements of single or a small number of genes are rare. Thus, this paper describes a straight-forward method to classify genomes into structural classes with the same order of genes and makes it possible to visualize reversed segments in pairs of genomes. The method can be readily applied to other species.

Genome and pan-genome analysis of a new exopolysaccharide-producing bacterium Pyschrobacillus sp. isolated from iron ores deposit and insights into iron uptake

Bacterial exopolysaccharides (EPS) have emerged as one of the key players in the field of heavy metal-contaminated environmental bioremediation. This study aimed to characterize and evaluate the metal biosorption potential of EPS produced by a novel Psychrobacillus strain, NEAU-3TGS, isolated from an iron ore deposit at Tamra iron mine, northern Tunisia. Genomic and pan-genomic analysis of NEAU-3TGS bacterium with nine validated published Psychrobacillus species was also performed. The results showed that the NEAU-3TGS genome (4.48 Mb) had a mean GC content of 36%, 4,243 coding sequences and 14 RNA genes. Phylogenomic analysis and calculation of nucleotide identity (ANI) values (less than 95% for new species with all strains) confirmed that NEAU-3TGS represents a potential new species. Pangenomic analysis revealed that Psychrobacillus genomic diversity represents an "open" pangenome model with 33,091 homologous genes, including 65 core, 3,738 shell, and 29,288 cloud genes. Structural EPS characterization by attenuated total reflectance-Fourier transform infrared (ATR-FTIR) spectroscopy showed uronic acid and α-1,4-glycosidic bonds as dominant components of the EPS. X-ray diffraction (XRD) analysis revealed the presence of chitin, chitosan, and calcite CaCO3 and confirmed the amorphous nature of the EPS. Heavy metal bioabsorption assessment showed that iron and lead were more adsorbed than copper and cadmium. Notably, the optimum activity was observed at 37°C, pH=7 and after 3 h contact of EPS with each metal. Genomic insights on iron acquisition and metabolism in Psychrobacillus sp. NEAU-3TGS suggested that no genes involved in siderophore biosynthesis were found, and only the gene cluster FeuABCD and trilactone hydrolase genes involved in the uptake of siderophores, iron transporter and exporter are present. Molecular modelling and docking of FeuA (protein peptidoglycan siderophore-binding protein) and siderophores ferrienterobactine [Fe+3 (ENT)]-3 and ferribacillibactine [Fe+3 (BB)]-3 ligand revealed that [Fe+3 (ENT)]-3 binds to Phe122, Lys127, Ile100, Gln314, Arg215, Arg217, and Gln252. Almost the same for [Fe+3 (ENT)]-3 in addition to Cys222 and Tyr229, but not Ile100.To the best of our knowledge, this is the first report on the characterization of EPS and the adsorption of heavy metals by Psychrobacillus species. The heavy metal removal capabilities may be advantageous for using these organisms in metal remediation.

Investigating Polyextremophilic Bacteria in Al Wahbah Crater, Saudi Arabia: A Terrestrial Model for Life on Saturn's Moon Enceladus

The study of extremophilic microorganisms has sparked interest in understanding extraterrestrial microbial life. Such organisms are fundamental for investigating life forms on Saturn's icy moons, such as Enceladus, which is characterized by potentially habitable saline and alkaline niches. Our study focused on the salt-alkaline soil of the Al Wahbah crater in Saudi Arabia, where we identified microorganisms that could be used as biological models to understand potential life on Enceladus. The search involved isolating 48 bacterial strains, sequencing the genomes of two thermo-haloalkaliphilic strains, and characterizing them for astrobiological application. A deeper understanding of the genetic composition and functional capabilities of the two novel strains of Halalkalibacterium halodurans provided valuable insights into their survival strategies and the presence of coding genes and pathways related to adaptations to environmental stressors. We also used mass spectrometry with a molecular network approach, highlighting various classes of molecules, such as phospholipids and nonproteinogenic amino acids, as potential biosignatures. These are essential features for understanding life's adaptability under extreme conditions and could be used as targets for biosignatures in upcoming missions exploring Enceladus' orbit. Furthermore, our study reinforces the need to look at new extreme environments on Earth that might contribute to the astrobiology field.

Pangenome Analysis Reveals Novel Contact-Dependent Growth Inhibition System and Phenazine Biosynthesis Operons in Proteus mirabilis BL95 That Are Located in An Integrative and Conjugative Element

Proteus mirabilis is a leading cause of urinary tract infections and a common commensal of the gastrointestinal tract. Our recent study (JB) showed that P. mirabilis strain BL95 employs a novel contact-dependent killing system against enteric bacteria in the mouse gut and in vitro. To uncover the genetic determinants of this system, we performed whole-genome sequencing of BL95 and compared it with 98 complete genomes of P. mirabilis. BL95 carries 56 coding sequences (CDSs) not found in other P. mirabilis. Over half of these unique genes are located on a novel integrative conjugative element (ICE) named ICEPm2, inserted in tRNA-Phe and exclusive to BL95. ICEPm2 has integration, conjugation, and DNA replication modules nearly identical to ICEPm1 (common in P. mirabilis), but ICEPm2 of BL95 carries two unique operons for P. mirabilis-a phenazine biosynthesis and a contact-dependent growth inhibition (CDI) system. ICEPm2 is absent in the P. mirabilis (AR_0156) closest to BL95 and it is present in the genomes of several Escherichia coli from mouse intestines, indicating its recent horizontal mobilization. BL95 shares over 100 genes of five different secretion systems with other P. mirabilis, mostly poorly studied, making a large pool of candidate genes for the contact-dependent growth inhibition.

Clinical manifestations, antimicrobial resistance and genomic feature analysis of multidrug-resistant Elizabethkingia strains

BACKGROUND

Elizabethkingia is emerging as an opportunistic pathogen in humans. The aim of this study was to investigate the clinical epidemiology, antimicrobial susceptibility, virulence factors, and genome features of Elizabethkingia spp.

METHODS

Clinical data from 71 patients who were diagnosed with Elizabethkingia-induced pneumonia and bacteremia between August 2019 and September 2021 were analyzed. Whole-genome sequencing was performed on seven isolates, and the results were compared with a dataset of 83 available Elizabethkingia genomes. Genomic features, Kyoto Encyclopedia of Genes and Genomes (KEGG) results and clusters of orthologous groups (COGs) were analyzed.

RESULTS

The mean age of the patients was 56.9 ± 20.7 years, and the in-hospital mortality rate was 29.6% (21/71). Elizabethkingia strains were obtained mainly from intensive care units (36.6%, 26/71) and emergency departments (32.4%, 23/71). The majority of the strains were isolated from respiratory tract specimens (85.9%, 61/71). All patients had a history of broad-spectrum antimicrobial exposure. Hospitalization for invasive mechanical ventilation or catheter insertion was found to be a risk factor for infection. The isolates displayed a high rate of resistance to cephalosporins and carbapenems, but all were susceptible to minocycline and colistin. Genomic analysis identified five β-lactamase genes (blaGOB, blaBlaB, blaCME, blaOXA, and blaTEM) responsible for β-lactam resistance and virulence genes involved in stress adaptation (ureB/G, katA/B, and clpP), adherence (groEL, tufA, and htpB) and immune modulation (gmd, tviB, cps4J, wbtIL, cap8E/D/G, and rfbC). Functional analysis of the COGs revealed that "metabolism" constituted the largest category within the core genome, while "information storage and processing" was predominant in both the accessory and unique genomes. The unique genes in our 7 strains were mostly enriched in KEGG pathways related to microRNAs in cancer, drug resistance (β-lactam and vancomycin), ABC transporters, biological metabolism and biosynthesis, and nucleotide excision repair mechanisms.

CONCLUSION

The Elizabethkingia genus exhibits multidrug resistance and carries carbapenemase genes. This study presents a comparative genomic analysis of Elizabethkingia, providing knowledge that facilitates a better understanding of this microorganism.

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.

Screening Marine Microbial Metabolites as Promising Inhibitors of Borrelia garinii: A Structural Docking Approach towards Developing Novel Lyme Disease Treatment

Lyme disease caused by the Borrelia species is a growing health concern in many parts of the world. Current treatments for the disease may have side effects, and there is also a need for new therapies that can selectively target the bacteria. Pathogens responsible for Lyme disease include B. burgdorferi, B. afzelii, and B. garinii. In this study, we employed structural docking-based screening to identify potential lead-like inhibitors against the bacterium. We first identified the core essential genome fraction of the bacterium, using 37 strains. Later, we screened a library of lead-like marine microbial metabolites (n = 4730) against the arginine deiminase (ADI) protein of Borrelia garinii. This protein plays a crucial role in the survival of the bacteria, and inhibiting it can kill the bacterium. The prioritized lead compounds demonstrating favorable binding energies and interactions with the active site of ADI were then evaluated for their drug-like and pharmacokinetic parameters to assess their suitability for development as drugs. Results from molecular dynamics simulation (100 ns) and other scoring parameters suggest that the compound CMNPD18759 (common name: aureobasidin; IUPAC name: 2-[(4R,6R)-4,6-dihydroxydecanoyl]oxypropan-2-yl (3S,5R)-3,5-dihydroxydecanoate) holds promise as a potential drug candidate for the treatment of Lyme disease, caused by B. garinii. However, further experimental studies are needed to validate the efficacy and safety of this compound in vivo.

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.

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.

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.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Comparative genomic analysis of the Dietzia genus: an insight into genomic diversity, and adaptation

Dietzia strains are widely distributed in the environment, presenting an opportunistic role, and some species have undetermined taxonomic characteristics. Here, we propose the existence of errors in the classification of species in this genus using comparative genomics. We performed ANI, dDDH, pangenome and genomic plasticity analyses better to elucidate the phylogenomic relationships between Dietzia strains. For this, we used 55 genomes of Dietzia downloaded from public databases that were combined with a newly sequenced. Sequence analysis of a phylogenetic tree based on genome similarity comparisons and dDDH, ANI analyses supported grouping different Dietzia species into four distinct groups. The pangenome analysis corroborated the classification of these groups, supporting the idea that some species of Dietzia could be reassigned in a possible classification into three distinct species, each containing less variability than that found within the global pangenome of all strains. Additionally, analysis of genomic plasticity based on groups containing Dietzia strains found differences in the presence and absence of symbiotic Islands and pathogenic islands related to their isolation site. We propose that the comparison of pangenome subsets together with phylogenomic approaches can be used as an alternative for the classification and differentiation of new species of the genus Dietzia.

Dissecting Complex Traits Using Omics Data: A Review on the Linear Mixed Models and Their Application in GWAS

Genome-wide association study (GWAS) is the most popular approach to dissecting complex traits in plants, humans, and animals. Numerous methods and tools have been proposed to discover the causal variants for GWAS data analysis. Among them, linear mixed models (LMMs) are widely used statistical methods for regulating confounding factors, including population structure, resulting in increased computational proficiency and statistical power in GWAS studies. Recently more attention has been paid to pleiotropy, multi-trait, gene-gene interaction, gene-environment interaction, and multi-locus methods with the growing availability of large-scale GWAS data and relevant phenotype samples. In this review, we have demonstrated all possible LMMs-based methods available in the literature for GWAS. We briefly discuss the different LMM methods, software packages, and available open-source applications in GWAS. Then, we include the advantages and weaknesses of the LMMs in GWAS. Finally, we discuss the future perspective and conclusion. The present review paper would be helpful to the researchers for selecting appropriate LMM models and methods quickly for GWAS data analysis and would benefit the scientific society.

Pathogenome comparison and global phylogeny of Escherichia coli ST1485 strains

Escherichia coli ST1485 strains belong to the clinically important phylogroup F and have disseminated worldwide in humans, animals, and the environment. Here, we elucidated the pathogenome of a global collection of E. coli ST1485 isolates from diverse sources retrieved from public databases and a high-quality sequenced complete genome of colistin-resistant E. coli strain CFSAN061771 isolated from raw milk cheese which designated as a reference strain. CFSAN061771 belongs to O83:H42-ST1485 pathotype and carries a conjugative ColV plasmid, pCFSAN061771_01, combining extraintestinal virulence genes (ompt, sitA, iroN, etsC, traT, cvaC, hylF, iss, tsh, mchf, iucC, iutA) with a multidrug resistance island (blaTEM-1, aph(6)-Id, aph(3″)-Ib, sul2, dfrA14). Comparative genomic analysis revealed a high frequency of pCFSAN061771_01-like plasmids in E. coli ST1485. A notable evolutionary genetic event in E. coli ST1485 strains is the acquisition of a pCFSAN061771_02-like plasmid, which confers resistance to several antimicrobials, tellurium, and quaternary ammonium compounds. The identical virulence and antibiotic resistance profiles identified in some human and animal strains are worrisome. This is the first study to emphasize the significance of E. coli ST1485 as a global high-risk virulent and multidrug-resistant clone with zoonotic potential.

Pangenomes-identified singletons for designing specific primers to identify bacterial strains in a plant growth-promoting consortium

BACKGROUND

The use of plant growth-promoting microorganisms represents a sustainable way to increase agricultural yields and plant health. Thus, the identification and tracking of these microorganisms are determinants for validating their positive effects on crops. Pangenomes allow the identification of singletons that can be used to design specific primers for the detection of the studied strains.

OBJECTIVE

This study aimed to establish a strategy based on the use of whole-genome sequencing and pangenomes for designing and validating primer sets for detecting Bacillus cabrialesii TE3^T, Priestia megaterium TRQ8, and Bacillus paralicheniformis TRQ65, a promising beneficial bacterial consortium for wheat.

METHODS AND RESULTS

The identification of singletons of TE3^T, TRQ8, and TRQ65 was performed by pangenomes using the Kbase platform and subsequently analyzed using BLAST®. The identified DNA regions were used for primer design in AlleleID version 7. Primers were validated by multiplex PCR using pure template DNA from each studied strain, combinations of two or three DNA from these strains, and DNA from agricultural soil samples enriched (and not) with the bacterial consortium. Here, we report the first design of primers capable of detecting and identifying the beneficial strains TE3^T, TRQ8, and TRQ65.

CONCLUSIONS

The use of pangenomes allowed the distinction of unique sequences that enables the design of primers for specific identification of the studied bacterial strains. This strategy can be widely used for the design of primer sets to detect other strains of interest for combating biopiracy, and commercial protection of biological products, among other applications.

Insights into genome evolution, pan-genome, and phylogenetic implication through mitochondrial genome sequence of Naegleria fowleri species

In the current study, we have systematically analysed the mitochondrial DNA (mtDNA) sequence of Naegleria fowleri (N. fowleri) isolate AY27, isolated from Karachi, Pakistan. The N. fowleri isolate AY27 has a circular mtDNA (49,541 bp), which harbours 69 genes (46 protein-coding genes, 21 tRNAs and 2 rRNAs). The pan-genome analysis of N. fowleri species showed a B_pan value of 0.137048, which implies that the pan-genome is open. KEGG classified core, accessory and unique gene clusters for human disease, metabolism, environmental information processing, genetic information processing and organismal system. Similarly, COG characterization of protein showed that core and accessory genes are involved in metabolism, information storages and processing, and cellular processes and signaling. The Naegleria species (n = 6) formed a total of 47 gene clusters; 42 single-copy gene clusters and 5 orthologous gene clusters. It was noted that 100% genes of Naegleria species were present in the orthogroups. We identified 44 single nucleotide polymorphisms (SNP) in the N. fowleri isolate AY27 mtDNA using N. fowleri strain V511 as a reference. Whole mtDNA phylogenetic tree analysis showed that N. fowleri isolates AY27 is closely related to N. fowleri (Accession no. JX174181.1). The ANI (Average Nucleotide Identity) values presented a much clear grouping of the Naegleria species compared to the whole mtDNA based phylogenetic analysis. The current study gives a comprehensive understanding of mtDNA architecture as well as a comparison of Naegleria species (N. fowleri and N. gruberi species) at the mitochondrial genome sequence level.

Genomic Characterization of Lactobacillus delbrueckii Strains with Probiotics Properties

Probiotics are health-beneficial microorganisms with mainly immunomodulatory and anti-inflammatory properties. Lactobacillus delbrueckii species is a common bacteria used in the dairy industry, and their benefits to hosting health have been reported. This study analyzed the core genome of nine strains of L. delbrueckii species with documented probiotic properties, focusing on genes related to their host health benefits. For this, a combined methodology including several software and databases (BPGA, SPAAN, BAGEL4, BioCyc, KEEG, and InterSPPI) was used to predict the most important characteristics related to L. delbrueckii strains probiose. Comparative genomics analyses revealed that L. delbrueckii probiotic strains shared essential genes related to acid and bile stress response and antimicrobial activity. Other standard features shared by these strains are surface layer proteins and extracellular proteins-encoding genes, with high adhesion profiles that interacted with human proteins of the inflammatory signaling pathways (TLR2/4-MAPK, TLR2/4-NF-κB, and NOD-like receptors). Among these, the PrtB serine protease appears to be a strong candidate responsible for the anti-inflammatory properties reported for these strains. Furthermore, genes with high proteolytic and metabolic activity able to produce beneficial metabolites, such as acetate, bioactive peptides, and B-complex vitamins were also identified. These findings suggest that these proteins can be essential in biological mechanisms related to probiotics’ beneficial effects of these strains in the host.

Pangenomics in Microbial and Crop Research: Progress, Applications, and Perspectives

Advances in sequencing technologies and bioinformatics tools have fueled a renewed interest in whole genome sequencing efforts in many organisms. The growing availability of multiple genome sequences has advanced our understanding of the within-species diversity, in the form of a pangenome. Pangenomics has opened new avenues for future research such as allowing dissection of complex molecular mechanisms and increased confidence in genome mapping. To comprehensively capture the genetic diversity for improving plant performance, the pangenome concept is further extended from species to genus level by the inclusion of wild species, constituting a super-pangenome. Characterization of pangenome has implications for both basic and applied research. The concept of pangenome has transformed the way biological questions are addressed. From understanding evolution and adaptation to elucidating host–pathogen interactions, finding novel genes or breeding targets to aid crop improvement to design effective vaccines for human prophylaxis, the increasing availability of the pangenome has revolutionized several aspects of biological research. The future availability of high-resolution pangenomes based on reference-level near-complete genome assemblies would greatly improve our ability to address complex biological problems.

Comparative pangenomics: analysis of 12 microbial pathogen pangenomes reveals conserved global structures of genetic and functional diversity

BACKGROUND

With the exponential growth of publicly available genome sequences, pangenome analyses have provided increasingly complete pictures of genetic diversity for many microbial species. However, relatively few studies have scaled beyond single pangenomes to compare global genetic diversity both within and across different species. We present here several methods for "comparative pangenomics" that can be used to contextualize multi-pangenome scale genetic diversity with gene function for multiple species at multiple resolutions: pangenome shape, genes, sequence variants, and positions within variants.

RESULTS

Applied to 12,676 genomes across 12 microbial pathogenic species, we observed several shared resolution-specific patterns of genetic diversity: First, pangenome openness is associated with species' phylogenetic placement. Second, relationships between gene function and frequency are conserved across species, with core genomes enriched for metabolic and ribosomal genes and accessory genomes for trafficking, secretion, and defense-associated genes. Third, genes in core genomes with the highest sequence diversity are functionally diverse. Finally, certain protein domains are consistently mutation enriched across multiple species, especially among aminoacyl-tRNA synthetases where the extent of a domain's mutation enrichment is strongly function-dependent.

CONCLUSIONS

These results illustrate the value of each resolution at uncovering distinct aspects in the relationship between genetic and functional diversity across multiple species. With the continued growth of the number of sequenced genomes, these methods will reveal additional universal patterns of genetic diversity at the pangenome scale.

Optimizing the Parametrization of Homologue Classification in the Pan-Genome Computation for a Bacterial Species: Case Study Streptococcus pyogenes

The paradigm shift associated with the introduction of the pan-genome concept has drawn the attention from singular reference genomes toward the actual sequence diversity within organism populations, strain collections, clades, etc. A single genome is no longer sufficient to describe bacteria of interest, but instead, the genomic repertoire of all existing strains is the key to the metabolic, evolutionary, or pathogenic potential of a species. The classification of orthologous genes derived from a collection of taxonomically related genome sequences is central to bacterial pan-genome computational analysis. In this work, we present a review of methods for computing pan-genome gene clusters including their comparative analysis for the case of Streptococcus pyogenes strain genomes. We exhaustively scanned the parametrization space of the homologue searching procedures and find optimal parameters (sequence identity (60%) and coverage (50-60%) in the pairwise alignment) for the orthologous clustering of gene sequences. We find that the sequence identity threshold influences the number of gene families ~3 times stronger than the sequence coverage threshold.

Comparative Genomic Analysis of Stenotrophomonas maltophilia Strain W18 Reveals Its Adaptative Genomic Features for Degrading Polycyclic Aromatic Hydrocarbons

Polycyclic aromatic hydrocarbons (PAHs) are hazardous pollutants that are ubiquitous in the environment. Numerous bacteria have evolved to have degrading genes or pathways to degrade PAHs. Stenotrophomonas maltophilia strain W18 was found to be able to degrade PAHs. Including 43 other complete genome sequences of S. maltophilia strains, we performed a comparative genomic analysis of 44 S. maltophilia strains by running OrthoFinder. A KEGG pathway enrichment analysis of environmental and clinical isolates of S. maltophilia revealed that environmental isolates tended to enhance gene functions such as "energy metabolism," "amino acid metabolism," "xenobiotic biodegradation and metabolism," and "folding, sorting, and degradation." The pangenome of the 44 S. maltophilia strains was open, while the core genome was estimated to reach a steady plateau. Based on gene annotations, we inferred that most of the degradation potential came from the core genome of S. maltophilia, while character genes and accessory genes also contributed to the degradation ability of S. maltophilia W18. The genes expression level of core genes, character genes and accessory genes were proved by RT-qPCR experiment, and accessory genes encoding alcohol dehydrogenase were upregulated most compared with genes with similar functions. We performed a credible comparative genomic analysis of S. maltophilia strains. S. maltophilia W18 was set as a model PAH-degrading bacterium of this species in this study, which would provide guidance for understanding and predicting the degradation mechanisms of other PAH-degrading S. maltophilia strains lacking complete genome data or waiting to be determined. IMPORTANCE This study provided the latest comparative genomic analysis on Stenotrophomonas maltophilia strains and focused on analyzing their genomic features that allow them to adapt to natural environments. In this study, we set S. maltophilia W18 as a typical PAH-degrading strain of this species. By discussing the genomic adaptative features of degrading PAH, we can predict genomic adaptative features of other S. maltophilia PAH-degrading strains since the core function of this species is stable. The gene functions of how S. maltophilia environmental isolates are enhanced for adaptation to various natural environments compared with clinical isolates have been revealed. Combined with a pangenome analysis and RT-qPCR results, we have proved that core genes, character genes, and accessory genes are all involved in PAH degradation. Accessory genes encoding alcohol dehydrogenase were upregulated most compared with core and character genes with similar functions, which suggests that PAH metabolization potential might be enhanced by horizontal gene transfer.

Comparative Genomics Reveals Insights into Induction of Violacein Biosynthesis and Adaptive Evolution in Janthinobacterium

Violacein has different bioactive properties conferring distinct selective advantages, such as defense from predation and interspecific competition. Adaptation of Janthinobacterium to diverse habitats likely leads to variation in violacein production among phylogenetically closely related species inhabiting different environments, yet genomic mechanisms and the influence of adaptive evolution underpinning violacein biosynthesis in Janthinobacterium are not clear. In this study, we performed genome sequencing, comparative genomic analysis, and phenotypic characterization to investigate genomic factors regulating violacein production in nine Janthinobacterium strains, including a type strain from soil and eight strains we isolated from terrestrial subsurface sediment and groundwater. Results show that although all nine Janthinobacterium strains are phylogenetically closely related and contain genes essential for violacein biosynthesis, they vary in carbon usage and violacein production. Sediment and groundwater strains are weak violacein producers and possess far fewer secondary metabolite biosynthesis genes, indicating genome adaptation compared to soil strains. Further examination suggests that quorum sensing (QS) may play an important role in regulating violacein in Janthinobacterium: the strains exhibiting strong potential in violacein production possess both N-acyl-homoserine lactone (AHL) QS and Janthinobacterium QS (JQS) systems in their genomes, while weaker violacein-producing strains harbor only the JQS system. Preliminary tests of spent media of two Janthinobacterium strains possessing both AHL QS and JQS systems support the potential role of AHLs in inducing violacein production in Janthinobacterium. Overall, results from this study reveal potential genomic mechanisms involved in violacein biosynthesis in Janthinobacterium and provide insights into evolution of Janthinobacterium for adaptation to oligotrophic terrestrial subsurface environment.

IMPORTANCE Phylogenetically closely related bacteria can thrive in diverse environmental habitats due to adaptive evolution. Genomic changes resulting from adaptive evolution lead to variations in cellular function, metabolism, and secondary metabolite biosynthesis. The most well-known secondary metabolite produced by Janthinobacterium is the purple-violet pigment violacein. To date, the mechanisms of induction of violacein biosynthesis in Janthinobacterium is not clear. Comparative genome analysis of closely related Janthinobacterium strains isolated from different environmental habitats not only reveals potential mechanisms involved in induction of violacein production by Janthinobacterium but also provides insights into the survival strategy of Janthinobacterium for adaptation to oligotrophic terrestrial subsurface environment.

Bacillus velezensis: A Treasure House of Bioactive Compounds of Medicinal, Biocontrol and Environmental Importance

Bacillus velezensis gram-positive bacterium, is frequently isolated from diverse niches mainly soil, water, plant roots, and fermented foods. B. velezensis is ubiquitous, non-pathogenic and endospore forming. Being frequently isolated from diverse plant holobionts it is considered host adapted microorganism and recognized of high economic importance given its ability to promote plant growth under diverse biotic and abiotic stress conditions. Additionally, the species suppress many plant diseases, including bacterial, oomycete, and fungal diseases. It is also able after plant host root colonization to induce unique physiological situation of host plant called primed state. Primed host plants are able to respond more rapidly and/or effectively to biotic or abiotic stress. Moreover, B. velezenis have the ability to resist diverse environmental stresses and help host plants to cope with, including metal and xenobiotic stresses. Within species B. velezensis strains have unique abilities allowing them to adopt different life styles. Strain level abilities knowledge is warranted and could be inferred using the ever-expanding new genomes list available in genomes databases. Pangenome analysis and subsequent identification of core, accessory and unique genomes is actually of paramount importance to decipher species full metabolic capacities and fitness across diverse environmental conditions shaping its life style. Despite the crucial importance of the pan genome, its assessment among large number of strains remains sparse and systematic studies still needed. Extensive knowledge of the pan genome is needed to translate genome sequencing efforts into developing more efficient biocontrol agents and bio-fertilizers. In this study, a genome survey of B. velezensis allowed us to (a) highlight B. velezensis species boundaries and show that Bacillus suffers taxonomic imprecision that blurs the debate over species pangenome; (b) identify drivers of their successful acquisition of specific life styles and colonization of new niches; (c) describe strategies they use to promote plant growth and development; (d) reveal the unlocked strain specific orphan secondary metabolite gene clusters (biosynthetic clusters with corresponding metabolites unknown) that product identification is still awaiting to amend our knowledge of their putative role in suppression of pathogens and plant growth promotion, and (e) to describe a dynamic pangenome with a secondary metabolite rich accessory genome.

First report from Bangladesh on genetic diversity of multidrug-resistant Pasteurella multocida type B:2 in fowl cholera

Background and Aim:

Fowl cholera (FC) caused by Pasteurella multocida is a highly contagious bacterial disease of global importance for poultry production. The severity and incidence of FC caused by P. multocida may vary considerably depending on several factors associated with the host (including species and age of infected birds), the environment, and the bacterial strain. This study aimed to investigate the genetic diversity of multidrug-resistant P. multocida strains isolated from FC outbreaks in laying hens from commercial farms of Bangladesh.

Materials and Methods:

We collected 57 samples of suspected FC, including 36 live and 21 dead laying hens. P. multocida isolates were characterized by biochemical and molecular-biological methods.

Results:

Twenty-two strains of P. multocida were isolated from these samples through phenotypic and genotypic characterization. The strains were grouped into two distinct random amplification of polymorphic DNA (RAPD) biotypes harboring a range of pathogenic genes; exbB, ompH, ptfA, nanB, sodC, and hgbA. In this study, 90.90% and 81.82% P. multocida strains were multidrug-resistant and biofilm formers, respectively. Whole-genome sequencing of the two representative RAPD phylotypes confirmed as P. multocida type B: L2:ST122, harboring a number of virulence factors-associated genes (VFGs), and antimicrobial resistance (AMR) genes (ARGs). In addition, pan-genome analysis revealed 90 unique genes in the genomes of P. multocida predicted to be associated with versatile metabolic functions, pathogenicity, virulence, and AMR.

Conclusion:

This is first-ever report on the association of P. multocida genotype B: L2:ST122 and related VFGs and ARGs in the pathogenesis of FC in laying hens. This study also provides a genetic context for future researches on the evolutionary diversity of P. multocida strains and their host adaptation.

Reference-agnostic representation and visualization of pan-genomes

BACKGROUND

The pan-genome of a species is the union of the genes and non-coding sequences present in all individuals (cultivar, accessions, or strains) within that species.

RESULTS

Here we introduce PGV, a reference-agnostic representation of the pan-genome of a species based on the notion of consensus ordering. Our experimental results demonstrate that PGV enables an intuitive, effective and interactive visualization of a pan-genome by providing a genome browser that can elucidate complex structural genomic variations.

CONCLUSIONS

The PGV software can be installed via conda or downloaded from https://github.com/ucrbioinfo/PGV . The companion PGV browser at http://pgv.cs.ucr.edu can be tested using example bed tracks available from the GitHub page.

Comparative genome analysis of Salmonella enterica serovar Gallinarum biovars Pullorum and Gallinarum decodes strain specific genes

Salmonella enterica serovar Gallinarum biovar Pullorum (bvP) and biovar Gallinarum (bvG) are the etiological agents of pullorum disease (PD) and fowl typhoid (FT) respectively, which cause huge economic losses to poultry industry especially in developing countries including India. Vaccination and biosecurity measures are currently being employed to control and reduce the S. Gallinarum infections. High endemicity, poor implementation of hygiene and lack of effective vaccines pose challenges in prevention and control of disease in intensively maintained poultry flocks. Comparative genome analysis unravels similarities and dissimilarities thus facilitating identification of genomic features that aids in pathogenesis, niche adaptation and in tracing of evolutionary history. The present investigation was carried out to assess the genotypic differences amongst S.enterica serovar Gallinarum strains including Indian strain S. Gallinarum Sal40 VTCCBAA614. The comparative genome analysis revealed an open pan-genome consisting of 5091 coding sequence (CDS) with 3270 CDS belonging to core-genome, 1254 CDS to dispensable genome and strain specific genes i.e. singletons ranging from 3 to 102 amongst the analyzed strains. Moreover, the investigated strains exhibited diversity in genomic features such as virulence factors, genomic islands, prophage regions, toxin-antitoxin cassettes, and acquired antimicrobial resistance genes. Core genome identified in the study can give important leads in the direction of design of rapid and reliable diagnostics, and vaccine design for effective infection control as well as eradication. Additionally, the identified genetic differences among the S. enterica serovar Gallinarum strains could be used for bacterial typing, structure based inhibitor development by future experimental investigations on the data generated.

Pan-Resistome Insights into the Multidrug Resistance of Acinetobacter baumannii

Acinetobacter baumannii is an important Gram-negative opportunistic pathogen that is responsible for many nosocomial infections. This etiologic agent has acquired, over the years, multiple mechanisms of resistance to a wide range of antimicrobials and the ability to survive in different environments. In this context, our study aims to elucidate the resistome from the A. baumannii strains based on phylogenetic, phylogenomic, and comparative genomics analyses. In silico analysis of the complete genomes of A. baumannii strains was carried out to identify genes involved in the resistance mechanisms and the phylogenetic relationships and grouping of the strains based on the sequence type. The presence of genomic islands containing most of the resistance gene repertoire indicated high genomic plasticity, which probably enabled the acquisition of resistance genes and the formation of a robust resistome. A. baumannii displayed an open pan-genome and revealed a still constant genetic permutation among their strains. Furthermore, the resistance genes suggest a specific profile within the species throughout its evolutionary history. Moreover, the current study performed screening and characterization of the main genes present in the resistome, which can be used in applied research to develop new therapeutic methods to control this important bacterial pathogen.

The Vaginal Microbial Signatures of Preterm Birth Delivery in Indian Women

Background

The incidence of preterm birth (PTB) in India is around 13%. Specific bacterial communities or individual taxon living in the vaginal milieu of pregnant women is a potential risk factor for PTB and may play an important role in its pathophysiology. Besides, bacterial taxa associated with PTB vary across populations.

Objective

Conduct a comparative analysis of vaginal microbiome composition and microbial genomic repertoires of women who enrolled in the Interdisciplinary Group for Advanced Research on Birth Outcomes – A DBT India Initiative (GARBH-Ini) pregnancy cohort to identify bacterial taxa associated with term birth (TB) and PTB in Indian women.

Methods

Vaginal swabs were collected during all three trimesters from 38 pregnant Indian women who delivered spontaneous term (n=20) and preterm (n=18) neonates. Paired-end sequencing of V3-V4 region of 16S rRNA gene was performed using the metagenomic DNA isolated from vaginal swabs (n=115). Whole genome sequencing of bacterial species associated with birth outcomes was carried out by shotgun method. Lactobacillus species were grown anaerobically in the De Man, Rogosa and Sharpe (MRS) agar culture medium for isolation of genomic DNA and whole genome sequencing.

Results

Vaginal microbiome of both term and preterm samples reveals similar alpha diversity indices. However, significantly higher abundance of Lactobacillus iners (p-value _{All_Trimesters}<0.02), Megasphaera sp (p-value_{1st_Trimester <}0.05), Gardnerella vaginalis (p-value_{2nd_Trimester}= 0.01) and Sneathia sanguinegens (p-value_{2nd_Trimester <}0.0001) were identified in preterm samples whereas higher abundance of L. gasseri (p-value_{3rd_Trimester} =0.010) was observed in term samples by Wilcoxon rank-sum test. The relative abundance of L. iners, and Megasphaera sp. were found to be significantly different over time between term and preterm mothers. Analyses of the representative genomes of L. crispatus and L. gasseri indicate presence of secretory transcriptional regulator and several ribosomally synthesized antimicrobial peptides correlated with anti-inflammatory condition in the vagina. These findings indicate protective role of L. crispatus and L. gasseri in reducing the risk of PTB.

Conclusion

Our findings indicate that the dominance of specific Lactobacillus species and few other facultative anaerobes are associated with birth outcomes.

Genomic signatures of adaptation to natural settings in non-typhoidal Salmonella enterica Serovars Saintpaul, Thompson and Weltevreden

Salmonella enterica is a pathogenic bacterium responsible for intestinal illness and systemic diseases such as typhoid and paratyphoid fevers. Among clinical manifestation classification, non-typhoidal Salmonella is mainly known as foodborne pathogen associated with the consumption of fecal contaminated food and water. Even though Salmonella hosts include humans and warm-blooded animals, it has been found in non-host environments as river water where the bacteria use different strategies to fitness the environment persisting and establishment. Now with the availability of WGS and bioinformatics tools, we can explore bacterial genomes with higher resolution to increase our understanding of specific genetic signatures among environmental and clinical isolates, being the goal of this work. Pangenome construction allowed the detection of specific environmental and clinical gene clusters related to metabolism and secretion systems as the main signature respectively. Specifically, D-galactonate degradation pathway was observed mainly in environmental genomes while T3SS and flagellum genes were detected for all clinical but not for all environmental isolates. Gene duplication and pseudogenes accumulation were detected as the main adaptation strategy for environmental isolates; thus, isolation source may play an important role in genome plasticity, conferring a selective advantage to survive and persist for environmental Salmonella isolates. Intact prophage sequences with cargo genes were observable for both isolation sources playing an important role in virulence contribution.

Biodegradation of phenol by a highly tolerant strain Rhodococcus ruber C1: Biochemical characterization and comparative genome analysis

A novel phenol-degrading strain was isolated and identified as Rhodococcus ruber C1. The degradation analysis shows that 1806 mg/L of phenol can be completely degraded by strain C1 within 38 h, and the maximum specific growth rate (μ_max=1.527 h^-1) and maximum specific phenol degradation rate (q_max=3.674 h^-1) indicate its excellent phenol metabolism capability. More importantly, phenol can be degraded by strain C1 in the temperature range of 20-45 °C within 72 h, and with longer degradation time, phenol can be completely degraded even at 10, 15 and 50 °C. The whole genome of strain C1 was sequenced, and a comparative genome analysis of strain C1 with 36 other genomes of Rhodococcus was performed. A remarkable gene family expansion occurred during the evolution of Rhodococcus, and a comprehensive evolutionary picture of Rhodococcus at genomic level was presented. Moreover, the copy number of genes involved in phenol metabolism was compared among genus Rhodococcus, and the results demonstrate high phenol degradation capability of strain C1 at genomic level. These findings suggest that Rhodococcus ruber C1 is a bacterium capable of degrading phenol efficiently in the temperature range of 10-50 °C.

Pan-genome: A promising resource for noncoding RNA discovery in plants

Plant genomes contain both protein-coding and noncoding sequences including transposable elements (TEs) and noncoding RNAs (ncRNAs). The ncRNAs are recognized as important elements that play fundamental roles in the structural organization and function of plant genomes. Despite various hypotheses, TEs are believed to be a major precursor of ncRNAs. Transposable elements are also prime factors that cause genomic variation among members of a species. Hence, TEs pose a major challenge in the discovery and analysis of ncRNAs. With the increase in the number of sequenced plant genomes, it is now accepted that a single reference genome is insufficient to represent the complete genomic diversity and contents of a species, and exploring the pan-genome of a species is critical. In this review, we summarize the recent progress in the field of plant pan-genomes. We also discuss TEs and their roles in ncRNA biogenesis and present our perspectives on the application of pan-genomes for the discovery of ncRNAs to fully explore and exploit their biological roles in plants.

Genomic Characteristics of Invasive Streptococcus pneumoniae Serotype 1 in New Caledonia Prior to the Introduction of PCV13

Streptococcus pneumoniae serotype 1 is a common cause of global invasive pneumococcal disease. In New Caledonia, serotype 1 is the most prevalent serotype and led to two major outbreaks reported in the 2000s. The pneumococcal conjugate vaccine 13 (PCV13) was introduced into the vaccination routine, intending to prevent the expansion of serotype 1 in New Caledonia. Aiming to provide a baseline for monitoring the post-PCV13 changes, we performed a whole-genome sequence analysis on 67 serotype 1 isolates collected prior to the PCV13 introduction. To highlight the S. pneumoniae serotype 1 population structure, we performed a multilocus sequence typing (MLST) analysis revealing that NC serotype 1 consisted of 2 sequence types: ST3717 and the highly dominant ST306. Both sequence types harbored the same resistance genes to beta-lactams, macrolide, streptogramin B, fluoroquinolone, and lincosamide antibiotics. We have also identified 36 virulence genes that were ubiquitous to all the isolates. Among these virulence genes, the pneumolysin sequence presented an allelic profile associated with disease outbreaks and reduced hemolytic activity. Moreover, recombination hotspots were identified in 4 virulence genes and more notably in the cps locus (cps2L), potentially leading to capsular switching, a major mechanism of the emergence of nonvaccine types. In summary, this study represents the first overview of the genomic characteristics of S. pneumoniae serotype 1 in New Caledonia prior to the introduction of PCV13. This preliminary description represents a baseline to assess the impact of PCV13 on serotype 1 population structure and genomic diversity.

Challenges in gene-oriented approaches for pangenome content discovery

Given a group of genomes, represented as the sets of genes that belong to them, the discovery of the pangenomic content is based on the search of genetic homology among the genes for clustering them into families. Thus, pangenomic analyses investigate the membership of the families to the given genomes. This approach is referred to as the gene-oriented approach in contrast to other definitions of the problem that takes into account different genomic features. In the past years, several tools have been developed to discover and analyse pangenomic contents. Because of the hardness of the problem, each tool applies a different strategy for discovering the pangenomic content. This results in a differentiation of the performance of each tool that depends on the composition of the input genomes. This review reports the main analysis instruments provided by the current state of the art tools for the discovery of pangenomic contents. Moreover, unlike previous works, the presented study compares pangenomic tools from a methodological perspective, analysing the causes that lead a given methodology to outperform other tools. The analysis is performed by taking into account different bacterial populations, which are synthetically generated by changing evolutionary parameters. The benchmarks used to compare the pangenomic tools, in addition to the computational pipeline developed for this purpose, are available at https://github.com/InfOmics/pangenes-review. Contact: V. Bonnici, R. Giugno Supplementary information: Supplementary data are available at Briefings in Bioinformatics online.

Targeting Plasmids to Limit Acquisition and Transmission of Antimicrobial Resistance

Antimicrobial resistance (AMR) is a significant global threat to both public health and the environment. The emergence and expansion of AMR is sustained by the enormous diversity and mobility of antimicrobial resistance genes (ARGs). Different mechanisms of horizontal gene transfer (HGT), including conjugation, transduction, and transformation, have facilitated the accumulation and dissemination of ARGs in Gram-negative and Gram-positive bacteria. This has resulted in the development of multidrug resistance in some bacteria. The most clinically significant ARGs are usually located on different mobile genetic elements (MGEs) that can move intracellularly (between the bacterial chromosome and plasmids) or intercellularly (within the same species or between different species or genera). Resistance plasmids play a central role both in HGT and as support elements for other MGEs, in which ARGs are assembled by transposition and recombination mechanisms. Considering the crucial role of MGEs in the acquisition and transmission of ARGs, a potential strategy to control AMR is to eliminate MGEs. This review discusses current progress on the development of chemical and biological approaches for the elimination of ARG carriers.

The pan-genome of Treponema pallidum reveals differences in genome plasticity between subspecies related to venereal and non-venereal syphilis

BACKGROUND

Spirochetal organisms of the Treponema genus are responsible for causing Treponematoses. Pathogenic treponemes is a Gram-negative, motile, spirochete pathogen that causes syphilis in human. Treponema pallidum subsp. endemicum (TEN) causes endemic syphilis (bejel); T. pallidum subsp. pallidum (TPA) causes venereal syphilis; T. pallidum subsp. pertenue (TPE) causes yaws; and T. pallidum subsp. Ccarateum causes pinta. Out of these four high morbidity diseases, venereal syphilis is mediated by sexual contact; the other three diseases are transmitted by close personal contact. The global distribution of syphilis is alarming and there is an increasing need of proper treatment and preventive measures. Unfortunately, effective measures are limited.

RESULTS

Here, the genome sequences of 53 T. pallidum strains isolated from different parts of the world and a diverse range of hosts were comparatively analysed using pan-genomic strategy. Phylogenomic, pan-genomic, core genomic and singleton analysis disclosed the close connection among all strains of the pathogen T. pallidum, its clonal behaviour and showed increases in the sizes of the pan-genome. Based on the genome plasticity analysis of the subsets containing the subspecies T pallidum subsp. pallidum, T. pallidum subsp. endemicum and T. pallidum subsp. pertenue, we found differences in the presence/absence of pathogenicity islands (PAIs) and genomic islands (GIs) on subsp.-based study.

CONCLUSIONS

In summary, we identified four pathogenicity islands (PAIs), eight genomic islands (GIs) in subsp. pallidum, whereas subsp. endemicum has three PAIs and seven GIs and subsp. pertenue harbours three PAIs and eight GIs. Concerning the presence of genes in PAIs and GIs, we found some genes related to lipid and amino acid biosynthesis that were only present in the subsp. of T. pallidum, compared to T. pallidum subsp. endemicum and T. pallidum subsp. pertenue.

Molecular and Genomic Characterization of PFAB2: A Non-virulent Bacillus anthracis Strain Isolated from an Indian Hot Spring

Background

Thermophilic bacilli in both aerobic or facultative anaerobic forms have been isolated for over a hundred years from different mesophilic or thermophilic environments as they are potential source of bioactive secondary metabolites. But the taxonomic resolution in the Bacillus genus at species or at strain level is very challenging for the insufficient divergence of the 16S rRNA genes. One such recurring problem is among Bacillus anthracis, B. cereus and B. thuringiensis. The disease-causing B. anthracis strains have their characteristic virulence factors coded in two well-known plasmids, namely pXO1 (toxin genes) and pXO2 (capsule genes).

Objective

The present study aimed at the molecular and genomic characterization of a recently reported thermophilic and environmental isolate of B. anthracis, strain PFAB2.

Methods

We performed comparative genomics between the PFAB2 genome and different strains of B. anthracis, along with closely related B. cereus strains.

Results

The pangenomic analysis suggests that the PFAB2 genome harbors no complete prophage genes. Cluster analysis of Bray-Kurtis similarity resemblance matrix revealed that gene content of PFAB2 is more closely related to other environmental strains of B. anthracis. The secretome analysis and the in vitro and in vivo pathogenesis experiments corroborate the avirulent phenotype of this strain. The most probable explanation for this phenotype is the apparent absence of plasmids harboring genes for capsule biosynthesis and toxins secretion in the draft genome. Additional features of PFAB2 are good spore-forming and germinating capabilities and rapid replication ability.

Conclusion

The high replication rate in a wide range of temperatures and culture media, the non-pathogenicity, the good spore forming capability and its genomic similarity to the Ames strain together make PFAB2 an interesting model strain for the study of the pathogenic evolution of B. anthracis.

Whole Genome Sequencing-Based Comparison of Food Isolates of Cronobacter sakazakii

Cronobacter sakazakii is an emerging foodborne pathogen, which is linked to life-threatening infections causing septicemia, meningitis, and necrotizing enterocolitis. These infections have been epidemiologically connected to ingestion of contaminated reconstituted powder infant formula. Even at low water activity C. sakazakii can survive for a long time; it is capable of protective biofilm formation and occasionally shows high virulence and pathogenicity even following stressful environmental conditions. Hence it is a challenging task for the food industry to control contamination of food ingredients and products through the entire production chain, since an increasing number of severe food-related outbreaks of C. sakazakii infections has been observed. The seemingly great capability of C. sakazakii to survive even strict countermeasures combined with its prevalence in many food ingredients requires a greater in depth understanding of its virulence factors to master the food safety issues related to this organism. In this context, we present the whole genome sequence (WGS) of two different C. sakazakii isolated from skimmed milk powder (C7) and ready-to-eat salad mix (C8), respectively. These are compared to other, already sequenced, C. sakazakii genomes. Sequencing of the fusA allele revealed that both isolates were C. sakazakii. We investigated the molecular characteristics of both isolates relevant for genes associated with pathogenesis and virulence factors, resistance to stressful environmental conditions (e.g., osmotic and heat), survival in desiccation as well as conducted a comparative genomic analysis. By using multi-locus sequence typing (MLST), the genetic type of both isolates is assessed and the number of unique genes is determined. DNA of C. sakazakii C8 is shown to hold a novel and unique sequence type; the number of unique genes identified in the genomic sequence of C. sakazakii C7 and C8 were 109 and 188, respectively. Some of the determined unique genes such as the rhs and VgrG genes are linked to the Type VI Secretion System cluster, which is associated with pathogenicity and virulence factors. Moreover, seven genes encoding for multi-drug resistance were found in both isolates. The finding of a number of genes linked to producing capsules and biofilm are likely related to the observed resistance to desiccation.

Proteome and allergenome of the European house dust mite Dermatophagoides pteronyssinus

The European house dust mite Dermatophagoides pteronyssinus is of significant medical importance as it is a major elicitor of allergic illnesses. In this analysis we have undertaken comprehensive bioinformatic and proteomic examination of Dermatophagoides pteronyssinus airmid, identified 12,530 predicted proteins and validated the expression of 4,002 proteins. Examination of homology between predicted proteins and allergens from other species revealed as much as 2.6% of the D. pteronyssinus airmid proteins may cause an allergenic response. Many of the potential allergens have evidence for expression (n = 259) and excretion (n = 161) making them interesting targets for future allergen studies. Comparative proteomic analysis of mite body and spent growth medium facilitated qualitative assessment of mite group allergen localisation. Protein extracts from house dust contain a substantial number of uncharacterised D. pteronyssinus proteins in addition to known and putative allergens. Novel D. pteronyssinus proteins were identified to be highly abundant both in house dust and laboratory cultures and included numerous carbohydrate active enzymes that may be involved in cuticle remodelling, bacteriophagy or mycophagy. These data may have clinical applications in the development of allergen-specific immunotherapy that mimic natural exposure. Using a phylogenomic approach utilising a supermatrix and supertree methodologies we also show that D. pteronyssinus is more closely related to Euroglyphus maynei than Dermatophagoides farinae.

Prediction of new vaccine targets in the core genome of Corynebacterium pseudotuberculosis through omics approaches and reverse vaccinology

Corynebacterium pseudotuberculosis is the etiologic agent of veterinary relevance diseases, such as caseous lymphadenitis, affecting different animal species causing damage to the global agribusiness. So far, there are no completely effective treatment methods to overcome the impacts caused by this pathogen. Several genomes of the species are deposited on public databases, allowing the execution of studies related to the pan-genomic approach. In this study, we used an integrated in silico workflow to prospect novel putative targets using the core genome, a set of shared genes among 65 C. pseudotuberculosis strains. Subsequently, through RNA-Seq data of the same abiotic stresses in two strains, we selected only induced genes to compose the reverse vaccinology workflow based in two different strategies. Our results predicted six probable antigens in both analysis, which indicates that they have a strong potential to be used in further studies as vaccine targets against this bacterium.

Micro-evolution of three Streptococcus species: selection, antigenic variation, and horizontal gene inflow

Background

The genus Streptococcus comprises pathogens that strongly influence the health of humans and animals. Genome sequencing of multiple Streptococcus strains demonstrated high variability in gene content and order even in closely related strains of the same species and created a newly emerged object for genomic analysis, the pan-genome. Here we analysed the genome evolution of 25 strains of Streptococcus suis, 50 strains of Streptococcus pyogenes and 28 strains of Streptococcus pneumoniae.

Results

Fractions of the pan-genome, unique, periphery, and universal genes differ in size, functional composition, the level of nucleotide substitutions, and predisposition to horizontal gene transfer and genomic rearrangements. The density of substitutions in intergenic regions appears to be correlated with selection acting on adjacent genes, implying that more conserved genes tend to have more conserved regulatory regions.

The total pan-genome of the genus is open, but only due to strain-specific genes, whereas other pan-genome fractions reach saturation. We have identified the set of genes with phylogenies inconsistent with species and non-conserved location in the chromosome; these genes are rare in at least one species and have likely experienced recent horizontal transfer between species. The strain-specific fraction is enriched with mobile elements and hypothetical proteins, but also contains a number of candidate virulence-related genes, so it may have a strong impact on adaptability and pathogenicity.

Mapping the rearrangements to the phylogenetic tree revealed large parallel inversions in all species. A parallel inversion of length 15 kB with breakpoints formed by genes encoding surface antigen proteins PhtD and PhtB in S. pneumoniae leads to replacement of gene fragments that likely indicates the action of an antigen variation mechanism.

Conclusions

Members of genus Streptococcus have a highly dynamic, open pan-genome, that potentially confers them with the ability to adapt to changing environmental conditions, i.e. antibiotic resistance or transmission between different hosts. Hence, integrated analysis of all aspects of genome evolution is important for the identification of potential pathogens and design of drugs and vaccines.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1403-6) contains supplementary material, which is available to authorized users.

Bacillus velezensis: phylogeny, useful applications, and avenues for exploitation

Some members of the Bacillus velezensis (Bv) group (e.g., Bv FZB42T and AS3.43) were previously assigned grouping with B. subtilis and B. amyloliquefaciens, based on the fact that they shared a 99% DNA-DNA percentage phylogenetic similarity. However, hinging on current assessments of the pan-genomic reassignments, the differing phylogenomic characteristics of Bv from B. subtilis and B. amyloliquefaciens are now better understood. Within this re-grouping/reassignment, the various strains within the Bv share a close phylogenomic resemblance, and a number of these strains have received a lot of attention in recent years, due to their genomic robustness, and the growing evidence for their possible utilization in the agricultural industry for managing plant diseases. Only a few applications for their use medicinally/pharmaceutically, environmentally, and in the food industry have been reported, and this may be due to the fact that the majority of those strains investigated are those typically occurring in soil. Although the intracellular unique biomolecules of Bv strains have been revealed via in silico genome modeling and investigated using transcriptomics and proteomics, a further inquisition into the Bv metabolome using newer technologies such as metabolomics could elucidate additional applications of this economically relevant Bacillus species, beyond that of primarily the agricultural sector.

Genomic Features, Comparative Genomics, and Antimicrobial Susceptibility Patterns of Elizabethkingia bruuniana

Elizabethkingia bruuniana is a novel species of the Elizabethkingia genus. There is scant information on this microorganism. Here, we report the whole-genome features and antimicrobial susceptibility patterns of E. bruuniana strain EM798-26. Elizabethkingia strain EM798-26 was initially identified as E. miricola. This isolate contained a circular genome of 4,393,011 bp. The whole-genome sequence-based phylogeny revealed that Elizabethkingia strain EM798-26 was in the same group of the type strain E. bruuniana G0146^T. Both in silico DNA-DNA hybridization and average nucleotide identity analysis clearly demonstrated that Elizabethkingia strain EM798-26 was a species of E. bruuniana. The pan-genome analysis identified 2,875 gene families in the core genome and 5,199 gene families in the pan genome of eight publicly available E. bruuniana genome sequences. The unique genes accounted for 0.2–12.1% of the pan genome in each E. bruuniana. A total of 59 potential virulence factor homologs were predicted in the whole-genome of E. bruuniana strain EM798–26. This isolate was nonsusceptible to multiple antibiotics, but susceptible to aminoglycosides, minocycline, and levofloxacin. The whole-genome sequence analysis of E. bruuniana EM798-26 revealed 29 homologs of antibiotic resistance-related genes. This study presents the genomic features of E. bruuniana. Knowledge of the genomic characteristics provides valuable insights into a novel species.

Ochrobactrum quorumnocens sp. nov., a quorum quenching bacterium from the potato rhizosphere, and comparative genome analysis with related type strains

Ochrobactrum spp. are ubiquitous bacteria attracting growing attention as important members of microbiomes of plants and nematodes and as a source of enzymes for biotechnology. Strain Ochrobactrum sp. A44^T was isolated from the rhizosphere of a field-grown potato in Gelderland, the Netherlands. The strain can interfere with quorum sensing (QS) of Gram-negative bacteria through inactivation of N-acyl homoserine lactones (AHLs) and protect plant tissue against soft rot pathogens, the virulence of which is governed by QS. Phylogenetic analysis based on 16S rRNA gene alone and concatenation of 16S rRNA gene and MLSA genes (groEL and gyrB) revealed that the closest relatives of A44^T are O. grignonense OgA9a^T, O. thiophenivorans DSM 7216^T, O. pseudogrignonense CCUG 30717^T, O. pituitosum CCUG 50899^T, and O. rhizosphaerae PR17^T. Genomes of all six type strains were sequenced, significantly expanding the possibility of genome-based analyses in Ochrobactrum spp. Average nucleotide identity (ANIb) and genome-to-genome distance (GGDC) values for A44^T and the related strains were below the single species thresholds (95% and 70%, respectively), with the highest scores obtained for O. pituitosum CCUG 50899^T (87.31%; 35.6%), O. rhizosphaerae PR17^T (86.80%; 34.3%), and O. grignonense OgA9a^T (86.30%; 33.6%). Distinction of A44^T from the related type strains was supported by chemotaxonomic and biochemical analyses. Comparative genomics revealed that the core genome for the newly sequenced strains comprises 2731 genes, constituting 50–66% of each individual genome. Through phenotype-to-genotype study, we found that the non-motile strain O. thiophenivorans DSM 7216^T lacks a cluster of genes related to flagella formation. Moreover, we explored the genetic background of distinct urease activity among the strains. Here, we propose to establish a novel species Ochrobactrum quorumnocens, with A44^T as the type strain (= LMG 30544^T = PCM 2957^T).

Genomic, Phenotypic, and Virulence Analysis of Streptococcus sanguinis Oral and Infective-Endocarditis Isolates

Streptococcus sanguinis, an abundant and benign inhabitant of the oral cavity, is an important etiologic agent of infective endocarditis (IE), particularly in people with predisposing cardiac valvular damage. Although commonly isolated from patients with IE, little is known about the factors that make any particular S. sanguinis isolate more virulent than another or, indeed, whether significant differences in virulence exist among isolates. In this study, we compared the genomes of a collection of S. sanguinis strains comprised of both oral isolates and bloodstream isolates from patients diagnosed with IE. Oral and IE isolates could not be distinguished by phylogenetic analyses, and we did not succeed in identifying virulence genes unique to the IE strains. We then investigated the virulence of these strains in a rabbit model of IE using a variation of the Bar-seq (barcode sequencing) method wherein we pooled the strains and used Illumina sequencing to count unique barcodes that had been inserted into each isolate at a conserved intergenic region. After we determined that several of the genome sequences were misidentified in GenBank, our virulence results were used to inform our bioinformatic analyses, identifying genes that may explain the heterogeneity in virulence. We further characterized these strains by assaying for phenotypes potentially contributing to virulence. Neither strain competition via bacteriocin production nor biofilm formation showed any apparent relationship with virulence. Increased cell-associated manganese was, however, correlated with blood isolates. These results, combined with additional phenotypic assays, suggest that S. sanguinis virulence is highly variable and results from multiple genetic factors.

High genomic variability in the plant pathogenic bacterium Pectobacterium parmentieri deciphered from de novo assembled complete genomes

BACKGROUND

Pectobacterium parmentieri is a newly established species within the plant pathogenic family Pectobacteriaceae. Bacteria belonging to this species are causative agents of diseases in economically important crops (e.g. potato) in a wide range of different environmental conditions, encountered in Europe, North America, Africa, and New Zealand. Severe disease symptoms result from the activity of P. parmentieri virulence factors, such as plant cell wall degrading enzymes. Interestingly, we observe significant phenotypic differences among P. parmentieri isolates regarding virulence factors production and the abilities to macerate plants. To establish the possible genomic basis of these differences, we sequenced 12 genomes of P. parmentieri strains (10 isolated in Poland, 2 in Belgium) with the combined use of Illumina and PacBio approaches. De novo genome assembly was performed with the use of SPAdes software, while annotation was conducted by NCBI Prokaryotic Genome Annotation Pipeline.

RESULTS

The pan-genome study was performed on 15 genomes (12 de novo assembled and three reference strains: P. parmentieri CFBP 8475T, P. parmentieri SCC3193, P. parmentieri WPP163). The pan-genome includes 3706 core genes, a high number of accessory (1468) genes, and numerous unique (1847) genes. We identified the presence of well-known genes encoding virulence factors in the core genome fraction, but some of them were located in the dispensable genome. A significant fraction of horizontally transferred genes, virulence-related gene duplications, as well as different CRISPR arrays were found, which can explain the observed phenotypic differences. Finally, we found also, for the first time, the presence of a plasmid in one of the tested P. parmentieri strains isolated in Poland.

CONCLUSIONS

We can hypothesize that a large number of the genes in the dispensable genome and significant genomic variation among P. parmentieri strains could be the basis of the potential wide host range and widespread diffusion of P. parmentieri. The obtained data on the structure and gene content of P. parmentieri strains enabled us to speculate on the importance of high genomic plasticity for P. parmentieri adaptation to different environments.

Exploration of Survival Traits, Probiotic Determinants, Host Interactions, and Functional Evolution of Bifidobacterial Genomes Using Comparative Genomics

Members of the genus Bifidobacterium are found in a wide-range of habitats and are used as important probiotics. Thus, exploration of their functional traits at the genus level is of utmost significance. Besides, this genus has been demonstrated to exhibit an open pan-genome based on the limited number of genomes used in earlier studies. However, the number of genomes is a crucial factor for pan-genome calculations. We have analyzed the pan-genome of a comparatively larger dataset of 215 members of the genus Bifidobacterium belonging to different habitats, which revealed an open nature. The pan-genome for the 56 probiotic and human-gut strains of this genus, was also found to be open. The accessory- and unique-components of this pan-genome were found to be under the operation of Darwinian selection pressure. Further, their genome-size variation was predicted to be attributed to the abundance of certain functions carried by genomic islands, which are facilitated by insertion elements and prophages. In silico functional and host-microbe interaction analyses of their core-genome revealed significant genomic factors for niche-specific adaptations and probiotic traits. The core survival traits include stress tolerance, biofilm formation, nutrient transport, and Sec-secretion system, whereas the core probiotic traits are imparted by the factors involved in carbohydrate- and protein-metabolism and host-immunomodulations.

Comparative genomic analysis of Staphylococcus lugdunensis shows a closed pan-genome and multiple barriers to horizontal gene transfer

Background

Coagulase negative staphylococci (CoNS) are commensal bacteria on human skin. Staphylococcus lugdunensis is a unique CoNS which produces various virulence factors and may, like S. aureus, cause severe infections, particularly in hospital settings. Unlike other staphylococci, it remains highly susceptible to antimicrobials, and genome-based phylogenetic studies have evidenced a highly conserved genome that distinguishes it from all other staphylococci.

Results

We demonstrate that S. lugdunensis possesses a closed pan-genome with a very limited number of new genes, in contrast to other staphylococci that have an open pan-genome. Whole-genome nucleotide and amino acid identity levels are also higher than in other staphylococci. We identified numerous genetic barriers to horizontal gene transfer that might explain this result. The S. lugdunensis genome has multiple operons encoding for restriction-modification, CRISPR/Cas and toxin/antitoxin systems. We also identified a new PIN-like domain-associated protein that might belong to a larger operon, comprising a metalloprotease, that could function as a new toxin/antitoxin or detoxification system.

Conclusion

We show that S. lugdunensis has a unique genome profile within staphylococci, with a closed pan-genome and several systems to prevent horizontal gene transfer. Its virulence in clinical settings does not rely on its ability to acquire and exchange antibiotic resistance genes or other virulence factors as shown for other staphylococci.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4978-1) contains supplementary material, which is available to authorized users.

Pan4Draft: A Computational Tool to Improve the Accuracy of Pan-Genomic Analysis Using Draft Genomes

High-throughput sequencing technologies are a milestone in molecular biology for facilitating great advances in genomics by enabling the deposit of large volumes of biological data to public databases. The availability of such data has made possible the comparative genomic analysis through pipelines, using the entire gene repertoire of genomes. However, a large number of unfinished genomes exist in public databases; their number is approximately 16-fold higher than the number of complete genomes, which creates bias during comparative analyses. Therefore, the present work proposes a new tool called Pan4Drafts, an automated pipeline for pan-genomic analysis of draft prokaryotic genomes to maximize the representation and accuracy of the gene repertoire of unfinished genomes by using reads from sequencing data. Pan4Draft allows to perform comparative analyses using different methodologies such as combining complete and draft genomes, using only draft genomes or only complete genomes. Pan4Draft is available at http://www.computationalbiology.ufpa.br/pan4drafts and the test dataset is available at https://sourceforge.net/projects/pan4drafts .

Diversity and evolution of the emerging Pandoraviridae family

With DNA genomes reaching 2.5 Mb packed in particles of bacterium-like shape and dimension, the first two Acanthamoeba-infecting pandoraviruses remained up to now the most complex viruses since their discovery in 2013. Our isolation of three new strains from distant locations and environments is now used to perform the first comparative genomics analysis of the emerging worldwide-distributed Pandoraviridae family. Thorough annotation of the genomes combining transcriptomic, proteomic, and bioinformatic analyses reveals many non-coding transcripts and significantly reduces the former set of predicted protein-coding genes. Here we show that the pandoraviruses exhibit an open pan-genome, the enormous size of which is not adequately explained by gene duplications or horizontal transfers. As most of the strain-specific genes have no extant homolog and exhibit statistical features comparable to intergenic regions, we suggest that de novo gene creation could contribute to the evolution of the giant pandoravirus genomes.

Giant viruses are visible by light microscopy and have unusually long genomes. Here, the authors report three new members of the Pandoraviridae family and investigate their evolution and diversity.

Variability of Bacterial Essential Genes Among Closely Related Bacteria: The Case of Escherichia coli

The definition of bacterial essential genes has been widely pursued using different approaches. Their study has impacted several fields of research such as synthetic biology, the construction of bacteria with minimal chromosomes, the search for new antibiotic targets, or the design of strains with biotechnological applications. Bacterial genomes are mosaics that only share a small subset of gene-sequences (core genome) even among members of the same species. It has been reported that the presence of essential genes is highly variable between closely related bacteria and even among members of the same species, due to the phenomenon known as “non-orthologous gene displacement” that refers to the coding for an essential function by genes with no sequence homology due to horizontal gene transfer (HGT). The existence of dormant forms among bacteria and the high incidence of HGT have been proposed to be driving forces of bacterial evolution, and they might have a role in the low level of conservation of essential genes among related bacteria by non-orthologous gene displacement, but this correlation has not been recognized. The aim of this mini-review is to give a brief overview of the approaches that have been taken to define and study essential genes, and the implications of non-orthologous gene displacement in bacterial evolution, focusing mainly in the case of Escherichia coli. To this end, we reviewed the available literature, and we searched for the presence of the essential genes defined by mutagenesis in the genomes of the 63 best-sequenced E. coli genomes that are available in NCBI database. We could not document specific cases of non-orthologous gene displacement among the E. coli strains analyzed, but we found that the quality of the genome-sequences in the database is not enough to make accurate predictions about the conservation of essential-genes among members of this bacterial species.