Search

Bacillus subtilis currently encompasses four subspecies, Bacillus subtilis subsp. subtilis, Bacillus subtilis subsp. inaquosorum, Bacillus subtilis subsp. spizizenii and Bacillus subtilis subsp. stercoris. Several studies based on genomic comparisons have suggested these subspecies should be promoted to species status. Previously, one of the main reasons for leaving them as subspecies was the lack of distinguishing phenotypes. In this study, we used comparative genomics to determine the genes unique to each subspecies and used these to lead us to the unique phenotypes. The results show that one difference among the subspecies is they produce different bioactive secondary metabolites. B. subtilis subsp. spizizenii is shown conserve the genes to produce mycosubtilin, bacillaene and 3,3′-neotrehalosadiamine. B. subtilis subsp. inaquosorum is shown conserve the genes to produce bacillomycin F, fengycin and an unknown PKS/NRPS cluster. B. subtilis subsp. stercoris is shown conserve the genes to produce fengycin and an unknown PKS/NRPS cluster. While B. subtilis subsp. subtilis is shown to conserve the genes to produce 3,3′-neotrehalosadiamine. In addition, we update the chemotaxonomy and phenotyping to support their promotion to species status. [ABSTRACT FROM AUTHOR]

Campylobacteriosis, primarily caused by Campylobacter jejuni and C. coli , is the main bacterial zoonosis worldwide. While poultry is recognized as the main reservoir, bovines are considered another important reservoir for Campylobacter spp. found in human infections. In contrast to chicken, retail beef is seldom contaminated by Campylobacter species. However, beef liver is recognized to be frequently contaminated and has been linked to human infections via epidemiological investigations. Our aims were to evaluate the prevalence of Campylobacter spp. inside and on the surface of beef liver pieces at retail in Finland and to analyse the population in more detail using whole genome sequencing (WGS) to assess the public health relevance. A total of 90 retail beef livers were studied using both enrichment of the external peptone-saline rinse of the liver piece and direct culture from the inside after surface sterilization. Furthermore, 46 of the livers were also studied using direct culture of retail beef liver juice collected from the bottom of the consumer package to estimate the concentration of Campylobacter species. Overall, 44 (49 %) of the samples were positive for Campylobacter species, C. jejuni , C. fetus and C. lari being identified in 42 %, 8.9 % and 1.1 % of the samples, respectively. Direct culture of retail liver juice was a sensitive and convenient method for Campylobacter spp. detection, resulting in 48 % prevalence and a mean concentration of 49 cfu/ml (maximum 335 cfu/ml). Two samples (2.2 %), containing large hepatic ducts, were positive for C. jejuni internally, representing multilocus sequence typing (MLST) sequence type ST-19 and ST-21. WGS, core genome phylogeny and core genome MLST revealed that in most cases only one clearly distinct clone of clinically relevant C. jejuni or C. fetus was isolated from a single lot of samples. However, in some cases several distinct clones were identified simultaneously even from a single liver piece. In epidemiological investigations, it is thus highly advisable to genotype multiple isolates to capture the whole diversity of Campylobacter spp. from suspected food sources. Good kitchen hygiene, avoidance of cross-contamination and thorough cooking are important for limiting the transmission of campylobacteriosis. • Retail liver also internally contaminated with Campylobacter spp. • Campylobacter spp. concentration in retail liver juice correlated with shelf-life. • Epidemiological implications with multiple C. jejuni clones in a single lot • All C. fetus subsp. fetus isolates predicted susceptibility to common antimicrobials. [ABSTRACT FROM AUTHOR]

Since 2010, outbreaks of haemorrhagic septicaemia (HS) caused by Pasteurella (P.) multocida capsular type B (PmB) emerged in Germany. In 2017, we noticed a close spatiotemporal relationship between HS outbreak sites and wolf (Canis lupus) territories. Thus, the main objectives of our study were to investigate the molecular epidemiology of German PmB-HS-isolates and to assess the role of wolves as putative vectors of this pathogen. We collected 83 PmB isolates from HS outbreaks that occurred between 2010 and 2019 and sampled 150 wolves, which were found dead in the years 2017 to 2019, revealing another three PmB isolates. A maximum-likelihood-based phylogeny of the core genomes of 65 PmB-HS-isolates and the three PmB-wolf-isolates showed high relatedness. Furthermore, all belonged to capsular:LPS:MLST genotype B:L2:ST122RIRDC and showed highly similar virulence gene profiles, but clustered separately from 35 global ST122RIRDC strains. Our data revealed that German HS outbreaks were caused by a distinct genomic lineage of PmB-ST122 strains, hinting towards an independent, ongoing epidemiologic event. We demonstrated for the first time, that carnivores, i.e., wolves, might harbour PmB as a part of their oropharyngeal microbiota. Furthermore, the results of our study imply that wolves can carry the pathogen over long distances, indicating a major role of that animal species in the ongoing epidemiological event of HS in Germany.

Pantoea stewartii subsp. indologenes (Psi) is a causative agent of leafspot on foxtail millet and pearl millet; however, novel strains were recently identified that are pathogenic on onions. Our recent host range evaluation study identified two pathovars; P. stewartii subsp. indologenes pv. cepacicola pv. nov. and P. stewartii subsp. indologenes pv. setariae pv. nov. that are pathogenic on onions and millets or on millets only, respectively. In the current study, we developed a pan-genome using the whole genome sequencing of newly identified/classified Psi strains from both pathovars [pv. cepacicola (n = 4) and pv. setariae (n = 13)]. The full spectrum of the pan-genome contained 7030 genes. Among these, 3546 (present in genomes of all 17 strains) were the core genes that were a subset of 3682 soft-core genes (present in ≥16 strains). The accessory genome included 1308 shell genes and 2040 cloud genes (present in ≤2 strains). The pan-genome showed a clear linear progression with >6000 genes, suggesting that the pan-genome of Psi is open. Comparative phylogenetic analysis showed differences in phylogenetic clustering of Pantoea spp. using PAVs/wgMLST approach in comparison with core genome SNPs-based phylogeny. Further, we conducted a horizontal gene transfer (HGT) study using Psi strains from both pathovars along with strains from other Pantoea species, namely, P. stewartii subsp. stewartii LMG 2715T, P. ananatis LMG 2665T, P. agglomerans LMG L15, and P. allii LMG 24248T. A total of 317 HGT events among four Pantoea species were identified with most gene transfer events occurring between Psi pv. cepacicola and Psi pv. setariae. Pan-GWAS analysis predicted a total of 154 genes, including seven gene-clusters, which were associated with the pathogenicity phenotype (necrosis on seedling) on onions. One of the gene-clusters contained 11 genes with known functions and was found to be chromosomally located.

Introduction. The invasive meningococcal disease (meningitis and/or septicemia) is actual problem of public health in Russia. Neisseria meningitidis isolates are classified into serogroups, PorA/FetA VRs, sequence types and clonal complexes. The growth of the invasive forms of meningococcal infection caused by isolates with «W: P1.5,2: F1-1: ST-11 (cc11)» profile requires attention for extended genotyping because the discriminating ability of classical MLST and antigens typing does not allow to answer the question about genetic and antigenic features of the pathogens and their epidemic potential. Materials and Methods. Four N. meningitidis serogroup W isolates associated with invasive meningococcal disease in Moscow (Russia) were characterized by next-generation sequencing. Three isolates were sequenced directly from cerebrospinal fluid samples and one -as a bacterial culture. All isolates were characterized earlier and the data were published in the PubMLST data base (id38565, id38573, id50225 and id50241). Genomic DNA was sequenced on Illumina MiSeq instrument. Results and Discussion. Obtained sequences allowed us to characterize four meningococci isolates for more than 1400 loci from the core genome MLST scheme. We have analyzed the core genome MLST scheme information about surface-antigen coding sequences. Housekeeping genes sequences were used to determine eMLST profile, ribosomal protein genes and some antibiotic resistance associated genes. We have characterized some ribosomal protein genes and antibiotic resistance associated genes. Based on eMLST profiles we noticed that there are at list two clones of N. meningitidis serogroup W inside complex ST-11/ET-37 clonal complex circulating in Moscow during 2016. An eMLST profile of isolates id50225 and id50241 differs in 3 loci out of 20. Application of the approach based on next-generation sequencing in routine epidemiological surveillance dramatically increases the amount of data and genotyping discriminating ability.

Abstract Background Meningococcal colonization is a prerequisite for transmission and disease, but the bacterium only very infrequently causes disease while asymptomatic carriage is common. Carriage is highly dynamic, showing a great variety across time and space within and across populations, but also within individuals. The understanding of genetic changes in the meningococcus during carriage, when the bacteria resides in its natural niche, is important for understanding not only the carriage state, but the dynamics of the entire meningococcal population. Results Paired meningococcal isolates, obtained from 50 asymptomatic carriers about 2 months apart were analyzed with whole genome sequencing (WGS). Phylogenetic analysis revealed that most paired isolates from the same individual were closely related, and the average and median number of allelic differences between paired isolates defined as the same strain was 35. About twice as many differences were seen between isolates from different individuals within the same sequence type (ST). In 8%, different strains were detected at different time points. A difference in ST was observed in 6%, including an individual who was found to carry three different STs over the course of 9 weeks. One individual carried different strains from the same ST. In total, 566 of 1605 cgMLST genes had undergone within-host genetic changes in one or more pairs. The most frequently changed cgMLST gene was relA that was changed in 47% of pairs. Across the whole genome, pilE, differed mostly, in 85% of the pairs. The most frequent mechanisms of genetic difference between paired isolates were phase variation and recombination, including gene conversion. Different STs showed variation with regard to which genes that were most frequently changed, mostly due to absence/presence of phase variation. Conclusions This study revealed within-host genetic differences in meningococcal isolates during short-term asymptomatic carriage. The most frequently changed genes were genes belonging to the pilin family, the restriction/modification system, opacity proteins and genes involved in glycosylation. Higher resolution genome-wide sequence typing is necessary to resolve the diversity of isolates and reveals genetic differences not discovered by traditional typing schemes, and would be the preferred choice of technology.

ABSTRACT There are more than 200 published species within the Lactobacillus genus complex (LGC), the majority of which have sequenced type strain genomes available. Although genome-based species delimitation cutoffs are accepted as the gold standard by the community, these are seldom actually checked for new or already published species. In addition, the availability of genome data is revealing inconsistencies in the species-level classification of many strains. We constructed a de novo species taxonomy for the LGC based on 2,459 publicly available genomes, using a 94% core nucleotide identity cutoff. We reconciled these de novo species with published species and subspecies names by (i) identifying genomes of type strains and (ii) comparing 16S rRNA genes of the genomes with 16S rRNA genes of type strains. We found that genomes within the LGC could be divided into 239 de novo species that were discontinuous and exclusive. Comparison of these de novo species to published species led to the identification of nine sets of published species that can be merged and one species that can be split. Further, we found at least eight de novo species that constitute new, unpublished species. Finally, we reclassified 74 genomes on the species level and identified for the first time the species of 98 genomes. Overall, the current state of LGC species taxonomy is largely consistent with genome-based species delimitation cutoffs. There are, however, exceptions that should be resolved to evolve toward a taxonomy where species share a consistent diversity in terms of sequence divergence. IMPORTANCE The Lactobacillus genus complex is a group of bacteria that constitutes an important source of strains with medical and food applications. The number of bacterial whole-genome sequences available for this taxon has been increasing rapidly in recent years. Despite this wealth of information, the species within this group are still largely defined by older techniques. Here, we constructed a completely new species-level taxonomy for the Lactobacillus genus complex based on ∼2,500 whole-genome sequences. As a result of this effort, we found that many genomes are not classified to their correct species, and we were able to correct these. In addition, we found that some published species are abnormally large, while others are too small. Finally, we discovered at least eight completely novel species that have not been published before. Our work will help the field to evolve toward a more meaningful and complete taxonomy, based on whole-genome sequences.

Streptomyces rimosus is best known as the primary source of the tetracycline class of antibiotics, most notably oxytetracycline, which have been widely used against many gram-positive and gram-negative pathogens and protozoan parasites. However, despite the medical and agricultural importance of S. rimosus, little is known of its evolutionary history and genome dynamics. In this study, we aim to elucidate the pan-genome characteristics and phylogenetic relationships of 32 S. rimosus genomes. The S. rimosus pan-genome contains more than 22,000 orthologous gene clusters, and approximately 8.8% of these genes constitutes the core genome. A large part of the accessory genome is composed of 9,646 strain-specific genes. S. rimosus exhibits an open pan-genome (decay parameter α = 0.83) and high gene diversity between strains (genomic fluidity φ = 0.12). We also observed strain-level variation in the distribution and abundance of biosynthetic gene clusters (BGCs) and that each individual S. rimosus genome has a unique repertoire of BGCs. Lastly, we observed variation in recombination, with some strains donating or receiving DNA more often than others, strains that tend to frequently recombine with specific partners, genes that often experience recombination more than others, and variable sizes of recombined DNA sequences. We conclude that the high levels of inter-strain genomic variation in S. rimosus is partly explained by differences in recombination among strains. These results have important implications on current efforts for natural drug discovery, the ecological role of strain-level variation in microbial populations, and addressing the fundamental question of why microbes have pan-genomes.

Understanding the structure of a bacterial population is essential in order to understand bacterial evolution. Estimating the core genome (those genes common to all, or nearly all, strains of a species) is a key component of such analyses. The size and composition of the core genome varies by dataset, but we hypothesized that the variation between different collections of the same bacterial species would be minimal. To investigate this, we analyzed the genome sequences of 3,118 pneumococci recovered from healthy individuals in Reykjavik (Iceland), Southampton (United Kingdom), Boston (United States), and Maela (Thailand). The analyses revealed a “supercore” genome (genes shared by all 3,118 pneumococci) of 558 genes, although an additional 354 core genes were shared by pneumococci from Reykjavik, Southampton, and Boston. Overall, the size and composition of the core and pan-genomes among pneumococci recovered in Reykjavik, Southampton, and Boston were similar. Maela pneumococci were distinctly different in that they had a smaller core genome and larger pan-genome. The pan-genome of Maela pneumococci contained several >25 Kb sequence regions (flanked by pneumococcal genes) that were homologous to genomic regions found in other bacterial species. Overall, our work revealed that some subsets of the global pneumococcal population are highly heterogeneous, and our hypothesis was rejected. This is an important finding in terms of understanding genetic variation among pneumococci and is also an essential point of consideration before generalizing the findings from a single dataset to the wider pneumococcal population.

Background The Streptococcus genus is relevant to both public health and food safety because of its ability to cause pathogenic infections. It is well-represented (>100 genomes) in publicly available databases. Streptococci are ubiquitous, with multiple sources of isolation, from human pathogens to dairy products. The Streptococcus genus has traditionally been classified by morphology, serum types, the 16S ribosomal RNA (rRNA) gene, and multi-locus sequence types subject to in-depth comparative genomic analysis. Methods Core and pan-genomes described the genomic diversity of 108 strains belonging to 16 Streptococcus species. The core genome nucleotide diversity was calculated and compared to phylogenomic distances within the genus Streptococcus. The core genome was also used as a resource to recruit metagenomic fragment reads from streptococci dominated environments. A conventional 16S rRNA gene phylogeny reconstruction was used as a reference to compare the resulting dendrograms of average nucleotide identity (ANI) and genome similarity score (GSS) dendrograms. Results The core genome, in this work, consists of 404 proteins that are shared by all 108 Streptococcus. The average identity of the pairwise compared core proteins decreases proportionally to GSS lower scores, across species. The GSS dendrogram recovers most of the clades in the 16S rRNA gene phylogeny while distinguishing between 16S polytomies (unresolved nodes). The GSS is a distance metric that can reflect evolutionary history comparing orthologous proteins. Additionally, GSS resulted in the most useful metric for genus and species comparisons, where ANI metrics failed due to false positives when comparing different species. Discussion Understanding of genomic variability and species relatedness is the goal of tools like GSS, which makes use of the maximum pairwise shared orthologous sequences for its calculation. It allows for long evolutionary distances (above species) to be included because of the use of amino acid alignment scores, rather than nucleotides, and normalizing by positive matches. Newly sequenced species and strains could be easily placed into GSS dendrograms to infer overall genomic relatedness. The GSS is not restricted to ubiquitous conservancy of gene features; thus, it reflects the mosaic-structure and dynamism of gene acquisition and loss in bacterial genomes.

Streptococcus uberis is a globally endemic and poorly controlled cause of bovine mastitis impacting the sustainability of the modern dairy industry. A core genome was derived from 579 newly sequenced S. uberis isolates, along with 305 publicly available genome sequences of S. uberis isolated from 11 countries around the world and used to develop a core genome multi-locus sequence typing (cgMLST) scheme. The S. uberis core genome comprised 1475 genes, and these were used to identify 1447 curated loci that were indexed into the cgMLST scheme. This was able to type 1012 of 1037 (>97  %) isolates used and differentiated the associated sequences into 932 discrete core genome sequence types (cgSTs). Analysis of the phylogenetic relationships of cgSTs revealed no clear clustering of isolates based on metadata such as disease status or year of isolation. Geographical clustering of cgSTs was limited to identification of a UK-centric clade, but cgSTs from UK isolates were also dispersed with those originating from other geographical regions across the entire phylogenetic topology. The cgMLST scheme offers a new tool for the detailed analysis of this globally important pathogen of dairy cattle. Initial analysis has re-emphasized and exemplified the genetically diverse nature of the global population of this opportunistic pathogen.

Vibrio parahaemolyticus is a significant cause of foodborne illness, and its incidence worldwide is on the rise. It is thus imperative to develop a straightforward and efficient method for typing strains of this pathogen. In this study, we conducted a pangenome analysis of 75 complete genomes of V. parahaemolyticus and identified the core gene mtlA with the highest degree of variation, which distinguished 44 strains and outperformed traditional seven-gene-based MLST when combined with aer , another core gene with high degree of variation. The mtlA gene had higher resolution to type strains with a close relationship compared to the traditional MLST genes in the phylogenetic tree built by core genomes. Strong positive selection was also detected in the gene mtlA (ω > 1), representing adaptive and evolution in response to the environment. Therefore, the panel of gene mtlA and aer may serve as a tool for the typing of V. parahaemolyticus , potentially contributing to the prevention and control of this foodborne disease., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2024 The Authors.)

Prokaryotic genomes constantly undergo gene flux via lateral gene transfer, generating a pangenome structure consisting of a conserved core genome surrounded by a more variable accessory genome shell. Over time, flux generates change in genome content. Here, we measure and compare the rate of genome flux for 5655 prokaryotic genomes as a function of amino acid sequence divergence in 36 universally distributed proteins of the informational core (IC). We find a clock of gene content change. The long-term average rate of gene content flux is remarkably constant across all higher prokaryotic taxa sampled, whereby the size of the accessory genome-the proportion of the genome harboring gene content difference for genome pairs-varies across taxa. The proportion of species-level accessory genes per genome, varies from 0% (Chlamydia) to 30%-33% (Alphaproteobacteria, Gammaproteobacteria, and Clostridia). A clock-like rate of gene content change across all prokaryotic taxa sampled suggest that pangenome structure is a general feature of prokaryotic genomes and that it has been in existence since the divergence of bacteria and archaea., (© The Author(s) 2024. Published by Oxford University Press on behalf of FEMS.)

Tuberculosis (TB) is an airborne communicable disease with high morbidity and mortality rates, especially in developing countries. The causal agents of TB belong to the complex Mycobacterium tuberculosis (MTBc), which is composed of different human and animal TB associated species. Some animal associated species have zoonotic potential and add to the burden of TB management. The BCG (“Bacillus Calmette-Guérin”) vaccine is widely used for the prevention against TB, but its use is limited in immunocompromised patients and animals due to the adverse effects and disseminated life-threatening complications. In this study, we aimed to carry out a comparative genome analysis between the human adapted species including BCG vaccine strains to identify and pinpoint the conserved genes related to the virulence across all the species, which could add a new value for vaccine development. For this purpose, the sequences of 183 Mycobacterium tuberculosis (MTB) strains were retrieved from the freely available WGS dataset at NCBI. The species included: 168 sensu stricto MTB species with other human MTB complex associated strains: M. tuberculosis var. africanum (3), M. tuberculosis var. bovis (2 draft genomes) and 10 BCG species, which enabled the analysis of core genome which contains the conserved genes and some virulence factor determinants. Further, a phylogenetic tree was constructed including the genomes of human (183); animals MTB adapted strains (6) and the environmental Mycobacterium strain “M. canettii”. Our results showed that the core genome consists of 1166 conserved genes among these species, which represents a small portion of the pangenome (7036 genes). The remaining genes in the pangenome (5870) are accessory genes, adding a high inter-species diversity. Further, the core genome includes several virulence-associated genes and this could explain the rare infectiousness potential of some attenuated vaccine strains in some patients. This study reveals that low number of conserved genes in human adapted MTBc species and high inter-species diversity of the pan-genome could be considered for vaccine candidate development.

Background: Cronobacter sakazakii is an emerging opportunistic bacterial pathogen known to cause neonatal and pediatric infections, including meningitis, necrotizing enterocolitis, and bacteremia. Multiple disease outbreaks of C. sakazakii have been documented in the past few decades, yet little is known of its genomic diversity, adaptation, and evolution. Here, we analyzed the pan-genome characteristics and phylogenetic relationships of 237 genomes of C. sakazakii and 48 genomes of related Cronobacter species isolated from diverse sources. Results: The C. sakazakii pan-genome contains 17,158 orthologous gene clusters, and approximately 19.5% of these constitute the core genome. Phylogenetic analyses reveal the presence of at least ten deep branching monophyletic lineages indicative of ancestral diversification. We detected enrichment of functions involved in proton transport and rotational mechanism in accessory genes exclusively found in human-derived strains. In environment-exclusive accessory genes, we detected enrichment for those involved in tryptophan biosynthesis and indole metabolism. However, we did not find significantly enriched gene functions for those genes exclusively found in food strains. The most frequently detected virulence genes are those that encode proteins associated with chemotaxis, enterobactin synthesis, ferrienterobactin transporter, type VI secretion system, galactose metabolism, and mannose metabolism. The genes fos which encodes resistance against fosfomycin, a broad-spectrum cell wall synthesis inhibitor, and mdf(A) which encodes a multidrug efflux transporter were found in nearly all genomes. We found that a total of 2991 genes in the pan-genome have had a history of recombination. Many of the most frequently recombined genes are associated with nutrient acquisition, metabolism and toxin production. Conclusions: Overall, our results indicate that the presence of a large accessory gene pool, ability to switch between ecological niches, a diverse suite of antibiotic resistance, virulence and niche-specific genes, and frequent recombination partly explain the remarkable adaptability of C. sakazakii within and outside the human host. These findings provide critical insights that can help define the development of effective disease surveillance and control strategies for Cronobacter-related diseases. [ABSTRACT FROM AUTHOR]

Background: Reverse vaccinology accelerates the discovery of potential vaccine candidates (PVCs) prior to experimental validation. Current programs typically use one bacterial proteome to identify PVCs through a filtering architecture using feature prediction programs or a machine learning approach. Filtering approaches may eliminate potential antigens based on limitations in the accuracy of prediction tools used. Machine learning approaches are heavily dependent on the selection of training datasets with experimentally validated antigens (positive control) and non-protective-antigens (negative control). The use of one or few bacterial proteomes does not assess PVC conservation among strains, an important feature of vaccine antigens. Results: We present ReVac, which implements both a panoply of feature prediction programs without filtering out proteins, and scoring of candidates based on predictions made on curated positive and negative control PVCs datasets. ReVac surveys several genomes assessing protein conservation, as well as DNA and protein repeats, which may result in variable expression of PVCs. ReVac's orthologous clustering of conserved genes, identifies core and dispensable genome components. This is useful for determining the degree of conservation of PVCs among the population of isolates for a given pathogen. Potential vaccine candidates are then prioritized based on conservation and overall feature-based scoring. We present the application of ReVac, applied to 69 Moraxella catarrhalis and 270 non-typeable Haemophilus influenzae genomes, prioritizing 64 and 29 proteins as PVCs, respectively. Conclusion: ReVac's use of a scoring scheme ranks PVCs for subsequent experimental testing. It employs a redundancy-based approach in its predictions of features using several prediction tools. The protein's features are collated, and each protein is ranked based on the scoring scheme. Multi-genome analyses performed in ReVac allow for a comprehensive overview of PVCs from a pan-genome perspective, as an essential pre-requisite for any bacterial subunit vaccine design. ReVac prioritized PVCs of two human respiratory pathogens, identifying both novel and previously validated PVCs. [ABSTRACT FROM AUTHOR]

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. • The core genome Corynebacterium pseudotuberculosis is comprised by 768 genes. • According to RNA-Seq data, 230 genes are induced in at least one stress condition. • Six genes were predicted as candidates in different reverse vaccinology pipelines. • These genes may contribute in further studies to develop new vaccines. [ABSTRACT FROM AUTHOR]

Taxonomic and phylogenetic relationships of Streptococcus mitis and Streptococcus oralis have been difficult to establish biochemically and genetically. We used core-genome analyses of S. mitis and S. oralis, as well as the closely related species Streptococcus pneumoniae and Streptococcus parasanguinis, to clarify the phylogenetic relationships between S. mitis and S. oralis, as well as within subclades of S. oralis. All S. mitis (n  = 67), S. oralis (n  = 89), S. parasanguinis (n  = 27), and 27 S. pneumoniae genome assemblies were downloaded from NCBI and reannotated. All genes were delineated into homologous clusters and maximum-likelihood phylogenies built from putatively nonrecombinant core gene sets. Population structure was determined using Bayesian genome clustering, and patristic distance was calculated between populations. Population-specific gene content was assessed using a phylogenetic-based genome-wide association approach. Streptococcus mitis and S. oralis formed distinct clades, but species mixing suggests taxonomic misassignment. Patristic distance between populations suggests that S. oralis subsp. dentisani is a distinct species, whereas S. oralis subsp. tigurinus and subsp. oralis are supported as subspecies, and that S. mitis comprises two subspecies. None of the genes within the pan-genomes of S. mitis and S. oralis could be statistically correlated with either, and the dispensable genomes showed extensive variation among isolates. These are likely important factors contributing to established overlap in biochemical characteristics for these taxa. Based on core-genome analysis, the substructure of S. oralis and S. mitis should be redefined, and species assignments within S. oralis and S. mitis should be made based on whole-genome analysis to be robust to misassignment. [ABSTRACT FROM AUTHOR]

Streptomyces rimosus is best known as the primary source of the tetracycline class of antibiotics, most notably oxytetracycline, which have been widely used against many gram-positive and gram-negative pathogens and protozoan parasites. However, despite the medical and agricultural importance of S. rimosus , little is known of its evolutionary history and genome dynamics. In this study, we aim to elucidate the pan-genome characteristics and phylogenetic relationships of 32 S. rimosus genomes. The S. rimosus pan-genome contains more than 22,000 orthologous gene clusters, and approximately 8.8% of these genes constitutes the core genome. A large part of the accessory genome is composed of 9,646 strain-specific genes. S. rimosus exhibits an open pan-genome (decay parameter α = 0.83) and high gene diversity between strains (genomic fluidity φ = 0.12). We also observed strain-level variation in the distribution and abundance of biosynthetic gene clusters (BGCs) and that each individual S. rimosus genome has a unique repertoire of BGCs. Lastly, we observed variation in recombination, with some strains donating or receiving DNA more often than others, strains that tend to frequently recombine with specific partners, genes that often experience recombination more than others, and variable sizes of recombined DNA sequences. We conclude that the high levels of inter-strain genomic variation in S. rimosus is partly explained by differences in recombination among strains. These results have important implications on current efforts for natural drug discovery, the ecological role of strain-level variation in microbial populations, and addressing the fundamental question of why microbes have pan-genomes. [ABSTRACT FROM AUTHOR]

Background: A revolutionary diversion from classical vaccinology to reverse vaccinology approach has been observed in the last decade. The ever-increasing genomic and proteomic data has greatly facilitated the vaccine designing and development process. Reverse vaccinology is considered as a cost-effective and proficient approach to screen the entire pathogen genome. To look for broad-spectrum immunogenic targets and analysis of closelyrelated bacterial species, the assimilation of pangenome concept into reverse vaccinology approach is essential. The categories of species pangenome such as core, accessory, and unique genes sets can be analyzed for the identification of vaccine candidates through reverse vaccinology. Results: We have designed an integrative computational pipeline term as “PanRV” that employs both the pangenome and reverse vaccinology approaches. PanRV comprises of four functional modules including i) Pangenome Estimation Module (PGM) ii) Reverse Vaccinology Module (RVM) iii) Functional Annotation Module (FAM) and iv) Antibiotic Resistance Association Module (ARM). The pipeline is tested by using genomic data from 301 genomes of Staphylococcus aureus and the results are verified by experimentally known antigenic data. Conclusion: The proposed pipeline has proved to be the first comprehensive automated pipeline that can precisely identify putative vaccine candidates exploiting the microbial pangenome. PanRV is a Linux based package developed in JAVA language. [ABSTRACT FROM AUTHOR]

Abstract In 1986, Bouvet and Grimont delineated two related taxa of the genus Acinetobacter termed genospecies (GS) 8 and 9. They proposed the name Acinetobacter lwoffii for GS8, which included the supposed type strain (CIP 64.10). As the authenticity of CIP 64.10 was later questioned, this study aimed at reassessing the taxonomy of these genospecies. We investigated 52 strains of GS8 or GS9, including CIP 64.10 and the genuine type strain of A. lwoffii (NCTC 5866T). All strains were subjected to the genus-wide comparative analyses of MALDI-TOF whole-cell mass spectra, rpoB gene sequences and metabolic traits while whole-genome sequences were analysed for 16 strains. The strains were classified into two distinct groups corresponding to GS8 (n = 15) and GS9 (n = 37). CIP 64.10 fell within GS8 whereas NCTC 5866T belonged to GS9. Intraspecies ANIb values for the genomes of GS8 (n = 6) and GS9 (n = 10) were ≥96.1% and ≥95.4%, respectively, whereas the ANIb values between them were 86.8–88.6%. Based on core genome phylogeny, GS8 and GS9 formed a distinct clade within the genus, with two respective, strongly supported subclades. GS8 and GS9 were similar in physiological and catabolic properties but were separable by MALDI-TOF MS. We conclude that the name A. lwoffii pertains to GS9 and not to GS8 as originally assumed and that these groups represent two species. We propose the name Acinetobacter pseudolwoffii sp. nov. for GS8, with ANC 5044T (= CCM 8638T = CCUG 67963T = CIP 111642T) as the type strain, and provide the emended description of A. lwoffii. [ABSTRACT FROM AUTHOR]

Understanding the structure of a bacterial population is essential in order to understand bacterial evolution. Estimating the core genome (those genes common to all, or nearly all, strains of a species) is a key component of such analyses. The size and composition of the core genome varies by dataset, but we hypothesized that the variation between different collections of the same bacterial species would be minimal. To investigate this, we analyzed the genome sequences of 3,118 pneumococci recovered from healthy individuals in Reykjavik (Iceland), Southampton (United Kingdom), Boston (United States), and Maela (Thailand). The analyses revealed a "supercore" genome (genes shared by all 3,118 pneumococci) of 558 genes, although an additional 354 core genes were shared by pneumococci from Reykjavik, Southampton, and Boston. Overall, the size and composition of the core and pan-genomes among pneumococci recovered in Reykjavik, Southampton, and Boston were similar. Maela pneumococci were distinctly different in that they had a smaller core genome and larger pan-genome. The pan-genome of Maela pneumococci contained several >25 Kb sequence regions (flanked by pneumococcal genes) that were homologous to genomic regions found in other bacterial species. Overall, our work revealed that some subsets of the global pneumococcal population are highly heterogeneous, and our hypothesis was rejected. This is an important finding in terms of understanding genetic variation among pneumococci and is also an essential point of consideration before generalizing the findings from a single dataset to the wider pneumococcal population. [ABSTRACT FROM AUTHOR]

The huge increase in the availability of bacterial genomes led us to a point in which we can investigate and query pan-genomes, for example, the full set of genes of a given bacterial species or clade. Here, we used a data set of 1,311 high-quality genomes from the human pathogen Pseudomonas aeruginosa, 619 of which were newly sequenced, to show that a pan-genomic approach can greatly refine the population structure of bacterial species, provide new insights to define species boundaries, and generate hypotheses on the evolution of pathogenicity. The 665-gene P. aeruginosa core genome presented here, which constitutes only 1% of the entire pan-genome, is the first to be in the same order of magnitude as the minimal bacterial genome and represents a conservative estimate of the actual core genome. Moreover, the phylogeny based on this core genome provides strong evidence for a five-group population structure that includes two previously undescribed groups of isolates. Comparative genomics focusing on antimicrobial resistance and virulence genes showed that variation among isolates was partly linked to this population structure. Finally, we hypothesized that horizontal gene transfer had an important role in this respect, and found a total of 3,010 putative complete and fragmented plasmids, 5% and 12% of which contained resistance or virulence genes, respectively. This work provides data and strategies to study the evolutionary trajectories of resistance and virulence in P. aeruginosa. [ABSTRACT FROM AUTHOR]

Background: The Streptococcus genus is relevant to both public health and food safety because of its ability to cause pathogenic infections. It is well-represented (>100 genomes) in publicly available databases. Streptococci are ubiquitous, with multiple sources of isolation, from human pathogens to dairy products. The Streptococcus genus has traditionally been classified by morphology, serum types, the 16S ribosomal RNA (rRNA) gene, and multi-locus sequence types subject to in-depth comparative genomic analysis. Methods: Core and pan-genomes described the genomic diversity of 108 strains belonging to 16 Streptococcus species. The core genome nucleotide diversity was calculated and compared to phylogenomic distances within the genus Streptococcus. The core genome was also used as a resource to recruit metagenomic fragment reads from streptococci dominated environments. A conventional 16S rRNA gene phylogeny reconstruction was used as a reference to compare the resulting dendrograms of average nucleotide identity (ANI) and genome similarity score (GSS) dendrograms. Results: The core genome, in this work, consists of 404 proteins that are shared by all 108 Streptococcus. The average identity of the pairwise compared core proteins decreases proportionally to GSS lower scores, across species. The GSS dendrogram recovers most of the clades in the 16S rRNA gene phylogeny while distinguishing between 16S polytomies (unresolved nodes). The GSS is a distance metric that can reflect evolutionary history comparing orthologous proteins. Additionally, GSS resulted in the most useful metric for genus and species comparisons, where ANI metrics failed due to false positives when comparing different species. Discussion: Understanding of genomic variability and species relatedness is the goal of tools like GSS, which makes use of the maximum pairwise shared orthologous sequences for its calculation. It allows for long evolutionary distances (above species) to be included because of the use of amino acid alignment scores, rather than nucleotides, and normalizing by positive matches. Newly sequenced species and strains could be easily placed into GSS dendrograms to infer overall genomic relatedness. The GSS is not restricted to ubiquitous conservancy of gene features; thus, it reflects the mosaic-structure and dynamism of gene acquisition and loss in bacterial genomes. [ABSTRACT FROM AUTHOR]

Pangenome analyses reveal major clues on evolutionary instances and critical genome core conservation. The order Rhizobiales encompasses several families with rather disparate ecological attitudes. Among them, Rhizobiaceae, Bradyrhizobiaceae, Phyllobacteriacreae and Xanthobacteriaceae, include members proficient in mutualistic symbioses with plants based on the bacterial conversion of N2 into ammonia (nitrogen-fixation). The pangenome of 12 nitrogen-fixing plant symbionts of the Rhizobiales was analyzed yielding total 37,364 loci, with a core genome constituting 700 genes. The percentage of core genes averaged 10.2% over single genomes, and between 5% to 7% were found to be plasmid-associated. The comparison between a representative reference genome and the core genome subset, showed the core genome highly enriched in genes for macromolecule metabolism, ribosomal constituents and overall translation machinery, while membrane/periplasm-associated genes, and transport domains resulted under-represented. The analysis of protein functions revealed that between 1.7% and 4.9% of core proteins could putatively have different functions.

Vibrio mimicus is an emerging pathogen, mainly associated with contaminated seafood consumption. However, little is known about its evolution, biodiversity, and pathogenic potential. This study analyzes the pan-, core, and accessory genomes of nine V. mimicus strains. The core genome yielded 2424 genes in chromosome I (ChI) and 822 genes in chromosome II (ChII), with an accessory genome comprising an average of 10.9% of the whole genome for ChI and 29% for ChII. Core genome phylogenetic trees were obtained, and V. mimicus ATCC-33654 strain was the closest to the outgroup in both chromosomes. Additionally, a phylogenetic study of eight conserved genes (ftsZ, gapA, gyrB, topA, rpoA, recA, mreB, and pyrH), including Vibrio cholerae, Vibrio parilis, Vibrio metoecus, and Vibrio caribbenthicus, clearly showed clade differentiation. The main virulence genes found in ChI corresponded with type I secretion proteins, extracellular components, flagellar proteins, and potential regulators, while, in ChII, the main categories were type-I secretion proteins, chemotaxis proteins, and antibiotic resistance proteins. The accessory genome was characterized by the presence of mobile elements and toxin encoding genes in both chromosomes. Based on the genome atlas, it was possible to characterize differential regions between strains. The pan-genome of V. mimicus encompassed 3539 genes for ChI and 2355 genes for ChII. These results give us an insight into the virulence and gene content of V. mimicus, as well as constitute the first approach to its diversity.

To identify genome-based features characteristic of the avian and human pathogen Chlamydia (C.) psittaci and related chlamydiae, we analyzed whole-genome sequences of 33 strains belonging to 12 species. Using a novel genome analysis tool termed Roary ILP Bacterial Annotation Pipeline (RIBAP), this panel of strains was shown to share a large core genome comprising 784 genes and representing approximately 80% of individual genomes. Analyzing the most variable genomic sites, we identified a set of features of C. psittaci that in its entirety is characteristic of this species: (i) a relatively short plasticity zone of less than 30,000 nt without a tryptophan operon (also in C. abortus, C. avium, C. gallinacea, C. pneumoniae), (ii) a characteristic set of of Inc proteins comprising IncA, B, C, V, X, Y (with homologs in C. abortus, C. caviae and C. felis as closest relatives), (iii) a 502-aa SinC protein, the largest among Chlamydia spp., and (iv) an elevated number of Pmp proteins of subtype G (14 in C. psittaci, 14 in Cand. C. ibidis). In combination with future functional studies, the common and distinctive criteria revealed in this study provide important clues for understanding the complexity of host-specific behavior of individual Chlamydia spp.

Mycobacterium abscessus is a multidrug resistant bacterium that causes pulmonary and extrapulmonary disease. The reported prevalence of pulmonary M. abscessus infections appears to be increasing in the United States (US) and around the world. In the last five years, multiple studies have utilized whole genome sequencing to investigate the genetic epidemiology of two clinically relevant subspecies, M. abscessus subsp. abscessus (MAB) and M. abscessus subsp. massiliense (MMAS). Phylogenomic comparisons of clinical isolates revealed that substantial proportions of patients have MAB and MMAS isolates that belong to genetically similar clusters also known as ‘dominant clones’. Unlike the genetic lineages of Mycobacterium tuberculosis that tend to be geographically clustered, the MAB and MMAS clones have been found in clinical populations from the US, Europe, Australia and South America. Moreover, the clones have been associated with worse clinical outcomes and show increased pathogenicity in macrophage and mouse models. While some have suggested that they may have spread locally and then globally through ‘indirect transmission’ within cystic fibrosis (CF) clinics, isolates of these clones have also been associated with sporadic pulmonary infections in non-CF patients and unrelated hospital-acquired soft tissue infections. M. abscessus has long been thought to be acquired from the environment, but the prevalence, exposure risk and environmental reservoirs of the dominant clones are currently not known. This review summarizes the genomic studies of M. abscessus and synthesizes the current knowledge surrounding the geographically diverse dominant clones identified from patient samples. Furthermore, it discusses the limitations of core genome comparisons for studying these genetically similar isolates and explores the breadth of accessory genome variation that has been observed to date. The combination of both core and accessory genome variation among these isolates may be the key to elucidating the origin, spread and evolution of these frequent genotypes.

Milk and dairy products are rich in nutrients and are therefore habitats for various microbiomes. However, the composition of nutrients can be quite diverse, in particular among the sulfur containing amino acids. In milk, methionine is present in a 25-fold higher abundance than cysteine. Interestingly, a fraction of strains of the species L. paracasei – a flavor-enhancing adjunct culture species – can grow in medium with methionine as the sole sulfur source. In this study, we focus on genomic and evolutionary aspects of sulfur dependence in L. paracasei strains. From 24 selected L. paracasei strains, 16 strains can grow in medium with methionine as sole sulfur source. We sequenced these strains to perform gene-trait matching. We found that one gene cluster – consisting of a cysteine synthase, a cystathionine lyase, and a serine acetyltransferase – is present in all strains that grow in medium with methionine as sole sulfur source. In contrast, strains that depend on other sulfur sources do not have this gene cluster. We expanded the study and searched for this gene cluster in other species and detected it in the genomes of many bacteria species used in the food production. The comparison to these species showed that two different versions of the gene cluster exist in L. paracasei which were likely gained in two distinct events of horizontal gene transfer. Additionally, the comparison of 62 L. paracasei genomes and the two versions of the gene cluster revealed that this gene cluster is mobile within the species.

Strains belonging to the genus Amycolatopsis are well known for the production of a number of important antimicrobials and other bioactive molecules. In this study, we have sequenced the genomes of five Amycolatopsis strains including Amycolatopsis circi DSM 45561T, Amycolatopsis palatopharyngis DSM 44832T and Amycolatopsis thermalba NRRL B-24845T. The genome sequences were analyzed with 52 other publically available Amycolatopsis genomes, representing 34 species, and 12 representatives from related genera including Saccharomonospora, Saccharopolyspora, Saccharothrix, Pseudonocardia and Thermobispora. Based on the core genome phylogeny, Amycolatopsis strains were subdivided into four major clades and several singletons. The genus Amycolatopsis is homogeneous with only three strains noted to group with other genera. Amycolatopsis halophila YIM93223T is quite distinct from other Amycolatopsis strains, both phylogenetically and taxonomically, and belongs to a distinct genus. In addition, Amycolatopsis palatopharyngis DSM 44832T and Amycolatopsis marina CGMCC4 3568T grouped in a clade with Saccharomonospora strains and showed similar taxogenomic differences to this genus as well as other Amycolatopsis strains. The study found a number of strains, particularly those identified as Amycolatopsis orientalis, whose incorrect identification could be resolved by taxogenomic analyses. Similarly, some unclassified strains could be assigned with species designations. The genome sequences of some strains that were independently sequenced by different laboratories were almost identical (99–100% average nucleotide and amino acid identities) consistent with them being the same strain, and confirming the reproducibility and robustness of genomic data. These analyses further demonstrate that whole genome sequencing can reliably resolve intra- and, inter-generic structures and should be incorporated into prokaryotic systematics.

The Lactobacillus genus encompasses a genetically and functionally diverse group of species, and contains many strains widely formulated in the human food supply chain as probiotics and starter cultures. Within this genetically expansive group, there are several distinct clades that have high levels of homology, one of which is the Lactobacillus acidophilus group. Of the uniting features, small genomes, low GC content, adaptation to dairy environments, and fastidious growth requirements, are some of the most defining characteristics of this group. To better understand what truly links and defines this clade, we sought to characterize the genomic organization and content of the genomes of several members of this group. Through core genome analysis we explored the synteny and intrinsic genetic underpinnings of the L. acidophilus clade, and observed key features related to the evolution and adaptation of these organisms. While genetic content is able to provide a large map of the potential of each organism, it does not always reflect their functionality. Through transcriptomic data we inferred the core transcriptome of the L. acidophilus complex to better define the true metabolic capabilities that unite this clade. Using this approach we have identified seven small ORFs that are both highly conserved and transcribed in diverse members of this clade and could be potential novel small peptide or untranslated RNA regulators. Overall, our results reveal the core features of the L. acidophilus complex and open new avenues for the enhancement and formulation and of next generation probiotics and starter cultures.