Defining the estimated core genome of bacterial populations using a Bayesian decision model

A core genome multi-locus sequence typing scheme for Streptococcus uberis: an evolution in typing a genetically diverse pathogen

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.

Exploring the biosynthetic gene clusters in Brevibacterium: a comparative genomic analysis of diversity and distribution

Exploring Brevibacterium strains from various ecosystems may lead to the discovery of new antibiotic-producing strains. Brevibacterium sp. H-BE7, a strain isolated from marine sediments from Northern Patagonia, Chile, had its genome sequenced to study the biosynthetic potential to produce novel natural products within the Brevibacterium genus. The genome sequences of 98 Brevibacterium strains, including strain H-BE7, were selected for a genomic analysis. A phylogenomic cladogram was generated, which divided the Brevibacterium strains into four major clades. A total of 25 strains are potentially unique new species according to Average Nucleotide Identity (ANIb) values. These strains were isolated from various environments, emphasizing the importance of exploring diverse ecosystems to discover the full diversity of Brevibacterium. Pangenome analysis of Brevibacterium strains revealed that only 2.5% of gene clusters are included within the core genome, and most gene clusters occur either as singletons or as cloud genes present in less than ten strains. Brevibacterium strains from various phylogenomic clades exhibit diverse BGCs. Specific groups of BGCs show clade-specific distribution patterns, such as siderophore BGCs and carotenoid-related BGCs. A group of clade IV-A Brevibacterium strains possess a clade-specific Polyketide synthase (PKS) BGCs that connects with phenazine-related BGCs. Within the PKS BGC, five genes, including the biosynthetic PKS gene, participate in the mevalonate pathway and exhibit similarities with the phenazine A BGC. However, additional core biosynthetic phenazine genes were exclusively discovered in nine Brevibacterium strains, primarily isolated from cheese. Evaluating the antibacterial activity of strain H-BE7, it exhibited antimicrobial activity against Salmonella enterica and Listeria monocytogenes. Chemical dereplication identified bioactive compounds, such as 1-methoxyphenazine in the crude extracts of strain H-BE7, which could be responsible of the observed antibacterial activity. While strain H-BE7 lacks the core phenazine biosynthetic genes, it produces 1-methoxyphenazine, indicating the presence of an unknown biosynthetic pathway for this compound. This suggests the existence of alternative biosynthetic pathways or promiscuous enzymes within H-BE7's genome.

Identifying the core genome of the nucleus-forming bacteriophage family and characterization of Erwinia phage RAY

We recently discovered that some bacteriophages establish a nucleus-like replication compartment (phage nucleus), but the core genes that define nucleus-based phage replication and their phylogenetic distribution were still to be determined. Here, we show that phages encoding the major phage nucleus protein chimallin share 72 conserved genes encoded within seven gene blocks. Of these, 21 core genes are unique to nucleus-forming phage, and all but one of these genes encode proteins of unknown function. We propose that these phages comprise a novel viral family we term Chimalliviridae. Fluorescence microscopy and cryoelectron tomography studies of Erwinia phage vB_EamM_RAY confirm that many of the key steps of nucleus-based replication are conserved among diverse chimalliviruses and reveal variations on this replication mechanism. This work expands our understanding of phage nucleus and PhuZ spindle diversity and function, providing a roadmap for identifying key mechanisms underlying nucleus-based phage replication.

Evolutionary Processes Driving the Rise and Fall of Staphylococcus aureus ST239, a Dominant Hybrid Pathogen

Selection plays a key role in the spread of antibiotic resistance, but the evolutionary drivers of clinically important resistant strains remain poorly understood. Here, we use genomic analyses and competition experiments to study Staphylococcus aureus ST239, a prominent MRSA strain that is thought to have been formed by large-scale recombination between ST8 and ST30. Genomic analyses allowed us to refine the hybrid model for the origin of ST239 and to date the origin of ST239 to 1920 to 1945, which predates the clinical introduction of methicillin in 1959. Although purifying selection has dominated the evolution of ST239, parallel evolution has occurred in genes involved in antibiotic resistance and virulence, suggesting that ST239 has evolved toward an increasingly pathogenic lifestyle. Crucially, ST239 isolates have low competitive fitness relative to both ST8 and ST30 isolates, supporting the idea that fitness costs have driven the demise of this once-dominant pathogen strain.

Helicobacter pylori treatment in the post-antibiotics era-searching for new drug targets

Abstract

Helicobacter pylori, a member of Epsilonproteobacteria, is a Gram-negative microaerophilic bacterium that colonizes gastric mucosa of about 50% of the human population. Although most infections caused by H. pylori are asymptomatic, the microorganism is strongly associated with serious diseases of the upper gastrointestinal tract such as chronic gastritis, peptic ulcer, duodenal ulcer, and gastric cancer, and it is classified as a group I carcinogen. The prevalence of H. pylori infections varies worldwide. The H. pylori genotype, host gene polymorphisms, and environmental factors determine the type of induced disease. Currently, the most common therapy to treat H. pylori is the first line clarithromycin–based triple therapy or a quadruple therapy replacing clarithromycin with new antibiotics. Despite the enormous recent effort to introduce new therapeutic regimens to combat this pathogen, treatment for H. pylori still fails in more than 20% of patients, mainly due to the increased prevalence of antibiotic resistant strains. In this review we present recent progress aimed at designing new anti-H. pylori strategies to combat this pathogen. Some novel therapeutic regimens will potentially be used as an extra constituent of antibiotic therapy, and others may replace current antibiotic treatments.

Key points

• Attempts to improve eradication rate of H. pylori infection.

• Searching for new drug targets in anti-Helicobacter therapies.

A Core Genome Multilocus Sequence Typing Scheme for Pseudomonas aeruginosa

Pseudomonas aeruginosa is a ubiquitous microorganism and an important opportunistic pathogen responsible for a broad spectrum of infections mainly in immunosuppressed and critically ill patients. Molecular investigations traditionally rely on pulsed field gel electrophoresis (PFGE) and multilocus sequence typing (MLST). In this work we propose a core genome multilocus sequence typing (cgMLST) scheme for P. aeruginosa, a methodology that combines traditional MLST principles with whole genome sequencing data. All publicly available complete P. aeruginosa genomes, representing the diversity of this species, were used to establish a cgMLST scheme targeting 2,653 genes. The scheme was then tested using genomes available at contig, chromosome and scaffold levels. The proposed cgMLST scheme for P. aeruginosa typed over 99% (2,314/2,325) of the genomes available for this study considering at least 95% of the cgMLST target genes present. The absence of a certain number gene targets at the threshold considered for both the creation and validation steps due to low genome sequence quality is possibly the main reason for this result. The cgMLST scheme was compared with previously published whole genome single nucleotide polymorphism analysis for the characterization of the population structure of the epidemic clone ST235 and results were highly similar. In order to evaluate the typing resolution of the proposed scheme, collections of isolates belonging to two important STs associated with cystic fibrosis, ST146 and ST274, were typed using this scheme, and ST235 isolates associated with an outbreak were evaluated. Besides confirming the relatedness of all the isolates, earlier determined by MLST, the higher resolution of cgMLST denotes that it may be suitable for surveillance programs, overcoming possible shortcomings of classical MLST. The proposed scheme is publicly available at: https://github.com/BioinformaticsHIAEMolecularMicrobiology/cgMLST-Pseudomonas-aeruginosa.

Genetic, structural, and functional diversity of low and high-affinity siderophores in strains of nitrogen fixing Azotobacter chroococcum

To increase iron (Fe) bioavailability in surface soils, microbes secrete siderophores, chelators with widely varying Fe affinities. Strains of the soil bacterium Azotobacter chroococcum (AC), plant-growth promoting rhizobacteria used as agricultural inoculants, require high Fe concentrations for aerobic respiration and nitrogen fixation. Recently, A. chroococcum str. NCIMB 8003 was shown to synthesize three siderophore classes: (1) vibrioferrin, a low-affinity α-hydroxy carboxylate (pFe = 18.4), (2) amphibactins, high-affinity tris-hydroxamates, and (3) crochelin A, a high-affinity siderophore with mixed Fe-chelating groups (pFe = 23.9). The relevance and specific functions of these siderophores in AC strains remain unclear. We analyzed the genome and siderophores of a second AC strain, A. chroococcum str. B3, and found that it also produces vibrioferrin and amphibactins, but not crochelin A. Genome comparisons indicate that vibrioferrin production is a vertically inherited, conserved strategy for Fe uptake in A. chroococcum and other species of Azotobacter. Amphibactin and crochelin biosynthesis reflects a more complex evolutionary history, shaped by vertical gene transfer, gene gain and loss through recombination at a genomic hotspot. We found conserved patterns of low vs. high-affinity siderophore production across strains: the low-affinity vibrioferrin was produced by mildly Fe limited cultures. As cells became more severely Fe starved, vibrioferrin production decreased in favor of high-affinity amphibactins (str. B3, NCIMB 8003) and crochelin A (str. NCIMB 8003). Our results show the evolution of low and high-affinity siderophore families and conserved patterns for their production in response to Fe bioavailability in a common soil diazotroph.

PPanGGOLiN: Depicting microbial diversity via a partitioned pangenome graph

The use of comparative genomics for functional, evolutionary, and epidemiological studies requires methods to classify gene families in terms of occurrence in a given species. These methods usually lack multivariate statistical models to infer the partitions and the optimal number of classes and don't account for genome organization. We introduce a graph structure to model pangenomes in which nodes represent gene families and edges represent genomic neighborhood. Our method, named PPanGGOLiN, partitions nodes using an Expectation-Maximization algorithm based on multivariate Bernoulli Mixture Model coupled with a Markov Random Field. This approach takes into account the topology of the graph and the presence/absence of genes in pangenomes to classify gene families into persistent, cloud, and one or several shell partitions. By analyzing the partitioned pangenome graphs of isolate genomes from 439 species and metagenome-assembled genomes from 78 species, we demonstrate that our method is effective in estimating the persistent genome. Interestingly, it shows that the shell genome is a key element to understand genome dynamics, presumably because it reflects how genes present at intermediate frequencies drive adaptation of species, and its proportion in genomes is independent of genome size. The graph-based approach proposed by PPanGGOLiN is useful to depict the overall genomic diversity of thousands of strains in a compact structure and provides an effective basis for very large scale comparative genomics. The software is freely available at https://github.com/labgem/PPanGGOLiN.

A proposed core genome scheme for analyses of the Salmonella genus

The salmonellae are found in a wide range of animal hosts and many food products for human consumption. Most cases of human disease are caused by S. enterica subspecies I; however as opportunistic pathogens the other subspecies (II-VI) and S. bongori are capable of causing disease. Loci that were not consistently present in all of the species and subspecies were removed from a previously proposed core genome scheme (EBcgMLSTv2.0), the removal of these 252 loci resulted in a core genus scheme (SalmcgMLSTv1.0). SalmcgMLSTv1.0 clustered isolates from the same subspecies more rapidly and more accurately grouped isolates from different subspecies when compared with EBcgMLSTv2.0. All loci within the EBcgMLSTv2.0 scheme were present in over 98% of S. enterica subspecies I isolates and should, therefore, continue to be used for subspecies I analyses, while the SalmcgMLSTv1.0 scheme is more appropriate for cross genus investigations.

Genomic Analyses of >3,100 Nasopharyngeal Pneumococci Revealed Significant Differences Between Pneumococci Recovered in Four Different Geographical Regions

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.

Identification of Streptococcus pneumoniae by a real-time PCR assay targeting SP2020

Real-time PCR targeting lytA (the major autolysin gene) and piaB (permease gene of the pia ABC transporter) are currently used as the gold-standard culture-independent assays for Streptococcus pneumoniae identification. We evaluated the performance of a new real-time PCR assay – targeting SP2020 (putative transcriptional regulator gene) – and compared its performance with the assays previously described. A collection of 150 pneumococci, 433 non-pneumococci and 240 polymicrobial samples (obtained from nasopharynx, oropharynx, and saliva; 80 from each site) was tested. SP2020 and lytA-CDC assays had the best performance (sensitivity of 100% for each compared to 95.3% for piaB). The specificity for lytA and piaB was 99.5% and for SP2020 was 99.8%. Misidentifications occurred for the three genes: lytA, piaB and SP2020 were found in non-pneumococcal strains; piaB was absent in some pneumococci including a serotype 6B strain. Combining lytA and SP2020 assays resulted in no misidentifications. Most polymicrobial samples (88.8%) yielded concordant results for the three molecular targets. The remaining samples seemed to contain non-typeable pneumococci (0.8%), and non-pneumococci positive for lytA (1.7%) or SP2020 (8.7%). We propose that combined detection of both lytA-CDC and SP2020 is a powerful strategy for the identification of pneumococcus either in pure cultures or in polymicrobial samples.

Genetic Diversity of norA, Coding for a Main Efflux Pump of Staphylococcus aureus

NorA is the best studied efflux system of Staphylococcus aureus and therefore frequently used as a model for investigating efflux-mediated resistance in this pathogen. NorA activity is associated with resistance to fluoroquinolones, several antiseptics and disinfectants and several reports have pointed out the role of efflux systems, including NorA, as a first-line response to antimicrobials in S. aureus. Genetic diversity studies of the gene norA have described three alleles; norAI, norAII and norAIII. However, the epidemiology of these alleles and their impact on NorA activity remains unclear. Additionally, increasing studies do not account for norA variability when establishing relations between resistance phenotypes and norA presence or reported absence, which actually corresponds, as we now demonstrate, to different norA alleles. In the present study we assessed the variability of the norA gene present in the genome of over 1,000 S. aureus isolates, corresponding to 112 S. aureus strains with whole genome sequences publicly available; 917 MRSA strains sourced from a London-based study and nine MRSA isolates collected in a major Hospital in Lisbon, Portugal. Our analyses show that norA is part of the core genome of S. aureus. It also suggests that occurrence of norA variants reflects the population structure of this major pathogen. Overall, this work highlights the ubiquitous nature of norA in S. aureus which must be taken into account when studying the role played by this important determinant on S. aureus resistance to antimicrobials.

Effect of Vaccination on Pneumococci Isolated from the Nasopharynx of Healthy Children and the Middle Ear of Children with Otitis Media in Iceland

Vaccination with pneumococcal conjugate vaccines (PCVs) disrupts the pneumococcal population. Our aim was to determine the impact of the 10-valent PCV on the serotypes, genetic lineages, and antimicrobial susceptibility of pneumococci isolated from children in Iceland.

Vaccination with pneumococcal conjugate vaccines (PCVs) disrupts the pneumococcal population. Our aim was to determine the impact of the 10-valent PCV on the serotypes, genetic lineages, and antimicrobial susceptibility of pneumococci isolated from children in Iceland. Pneumococci were collected between 2009 and 2017 from the nasopharynges of healthy children attending 15 day care centers and from the middle ears (MEs) of children with acute otitis media from the greater Reykjavik capital area. Isolates were serotyped and tested for antimicrobial susceptibility. Whole-genome sequencing (WGS) was performed on alternate isolates from 2009 to 2014, and serotypes and multilocus sequence types (STs) were extracted from the WGS data. Two study periods were defined: 2009 to 2011 (PreVac) and 2012 to 2017 (PostVac). The overall nasopharyngeal carriage rate was similar between the two periods (67.3% PreVac and 61.5% PostVac, P = 0.090). Vaccine-type (VT) pneumococci decreased and nonvaccine-type (NVT) pneumococci (serotypes 6C, 15A, 15B/C, 21, 22F, 23A, 23B, 35F, and 35B) significantly increased in different age strata post-PCV introduction. The total number of pneumococci recovered from ME samples significantly decreased as did the proportion that were VTs, although NVT pneumococci (6C, 15B/C, 23A, and 23B) increased significantly. Most serotype 6C pneumococci were multidrug resistant (MDR). Serotype 19F was the predominant serotype associated with MEs, and it significantly decreased post-PCV introduction: these isolates were predominantly MDR and of the Taiwan^19F-14 PMEN lineage. Overall, the nasopharyngeal carriage rate remained constant and the number of ME-associated pneumococci decreased significantly post-PCV introduction; however, there was a concomitant and statistically significant shift from VTs to NVTs in both collections of pneumococci.

Genome Sequencing Reveals a Large and Diverse Repertoire of Antimicrobial Peptides

Competition among bacterial members of the same ecological niche is mediated by bacteriocins: antimicrobial peptides produced by bacterial species to kill other bacteria. Bacteriocins are also promising candidates for novel antimicrobials. Streptococcus pneumoniae (the “pneumococcus”) is a leading cause of morbidity and mortality worldwide and a frequent colonizer of the human nasopharynx. Here, 14 newly discovered bacteriocin gene clusters were identified among >6,200 pneumococcal genomes. The molecular epidemiology of the bacteriocin clusters was investigated using a large global and historical pneumococcal dataset dating from 1916. These analyses revealed extraordinary bacteriocin diversity among pneumococci and the majority of bacteriocin clusters were also found in other streptococcal species. Genomic hotspots for the integration of different bacteriocin gene clusters were discovered. Experimentally, bacteriocin genes were transcriptionally active when the pneumococcus was under stress and when two strains were co-cultured in broth. These findings reveal much more diversity among bacterial defense mechanisms than previously appreciated, which fundamentally broaden our understanding of bacteriocins relative to intraspecies and interspecies nasopharyngeal competition and bacterial population structure.

Defining and Evaluating a Core Genome Multilocus Sequence Typing Scheme for Genome-Wide Typing of Clostridium difficile

Clostridium difficile, recently renamed Clostridioides difficile, is the most common cause of antibiotic-associated nosocomial gastrointestinal infections worldwide. To differentiate endogenous infections and transmission events, highly discriminatory subtyping is necessary. Today, methods based on whole-genome sequencing data are increasingly used to subtype bacterial pathogens; however, frequently a standardized methodology and typing nomenclature are missing. Here we report a core genome multilocus sequence typing (cgMLST) approach developed for C. difficile Initially, we determined the breadth of the C. difficile population based on all available MLST sequence types with Bayesian inference (BAPS). The resulting BAPS partitions were used in combination with C. difficile clade information to select representative isolates that were subsequently used to define cgMLST target genes. Finally, we evaluated the novel cgMLST scheme with genomes from 3,025 isolates. BAPS grouping (n = 6 groups) together with the clade information led to a total of 11 representative isolates that were included for cgMLST definition and resulted in 2,270 cgMLST genes that were present in all isolates. Overall, 2,184 to 2,268 cgMLST targets were detected in the genome sequences of 70 outbreak-associated and reference strains, and on average 99.3% cgMLST targets (1,116 to 2,270 targets) were present in 2,954 genomes downloaded from the NCBI database, underlining the representativeness of the cgMLST scheme. Moreover, reanalyzing different cluster scenarios with cgMLST were concordant to published single nucleotide variant analyses. In conclusion, the novel cgMLST is representative for the whole C. difficile population, is highly discriminatory in outbreak situations, and provides a unique nomenclature facilitating interlaboratory exchange.

Comparative analysis of core genome MLST and SNP typing within a European Salmonella serovar Enteritidis outbreak

Multi-country outbreaks of foodborne bacterial disease present challenges in their detection, tracking, and notification. As food is increasingly distributed across borders, such outbreaks are becoming more common. This increases the need for high-resolution, accessible, and replicable isolate typing schemes. Here we evaluate a core genome multilocus typing (cgMLST) scheme for the high-resolution reproducible typing of Salmonella enterica (S. enterica) isolates, by its application to a large European outbreak of S. enterica serovar Enteritidis. This outbreak had been extensively characterised using single nucleotide polymorphism (SNP)-based approaches. The cgMLST analysis was congruent with the original SNP-based analysis, the epidemiological data, and whole genome MLST (wgMLST) analysis. Combination of the cgMLST and epidemiological data confirmed that the genetic diversity among the isolates predated the outbreak, and was likely present at the infection source. There was consequently no link between country of isolation and genetic diversity, but the cgMLST clusters were congruent with date of isolation. Furthermore, comparison with publicly available Enteritidis isolate data demonstrated that the cgMLST scheme presented is highly scalable, enabling outbreaks to be contextualised within the Salmonella genus. The cgMLST scheme is therefore shown to be a standardised and scalable typing method, which allows Salmonella outbreaks to be analysed and compared across laboratories and jurisdictions.

•

cgMLST is proposed as a universal typing scheme for Salmonella.

•

cgMLST is congruent with SNP analyses and easier to implement across laboratories.

•

Genomic data are consistent with the epidemiology of the outbreak.

Core Genome Multilocus Sequence Typing Scheme for Stable, Comparative Analyses of Campylobacter jejuni and C. coli Human Disease Isolates

Human campylobacteriosis, caused by Campylobacter jejuni and C. coli, remains a leading cause of bacterial gastroenteritis in many countries, but the epidemiology of campylobacteriosis outbreaks remains poorly defined, largely due to limitations in the resolution and comparability of isolate characterization methods. Whole-genome sequencing (WGS) data enable the improvement of sequence-based typing approaches, such as multilocus sequence typing (MLST), by substantially increasing the number of loci examined. A core genome MLST (cgMLST) scheme defines a comprehensive set of those loci present in most members of a bacterial group, balancing very high resolution with comparability across the diversity of the group. Here we propose a set of 1,343 loci as a human campylobacteriosis cgMLST scheme (v1.0), the allelic profiles of which can be assigned to core genome sequence types. The 1,343 loci chosen were a subset of the 1,643 loci identified in the reannotation of the genome sequence of C. jejuni isolate NCTC 11168, chosen as being present in >95% of draft genomes of 2,472 representative United Kingdom campylobacteriosis isolates, comprising 2,207 (89.3%) C. jejuni isolates and 265 (10.7%) C. coli isolates. Validation of the cgMLST scheme was undertaken with 1,478 further high-quality draft genomes, containing 150 or fewer contiguous sequences, from disease isolate collections: 99.5% of these isolates contained ≥95% of the 1,343 cgMLST loci. In addition to the rapid and effective high-resolution analysis of large numbers of diverse isolates, the cgMLST scheme enabled the efficient identification of very closely related isolates from a well-defined single-source campylobacteriosis outbreak.

Pathogenicity-associated protein domains: The fiercely-conserved evolutionary signatures

Proteins have highly conserved domains that determine their functionality. Out of the thousands of domains discovered so far across all living forms, some of the predominant clinically-relevant domains include IENR1, HNHc, HELICc, Pro-kuma_activ, Tryp_SPc, Lactamase_B, PbH1, ChtBD3, CBM49, acidPPc, G3P_acyltransf, RPOL8c, KbaA, HAMP, HisKA, Hr1, Dak2, APC2, Citrate_ly_lig, DALR, VKc, YARHG, WR1, PWI, ZnF_BED, TUDOR, MHC_II_beta, Integrin_B_tail, Excalibur, DISIN, Cadherin, ACTIN, PROF, Robl_LC7, MIT, Kelch, GAS2, B41, Cyclin_C, Connexin_CCC, OmpH, Bac_rhodopsin, AAA, Knot1, NH, Galanin, IB, Elicitin, ACTH, Cache_2, CHASE, AgrB, PRP, IGR, and Antimicrobial21. These domains are distributed in nucleases/helicases, proteases, esterases, lipases, glycosylase, GTPases, phosphatases, methyltransferases, acyltransferase, acetyltransferase, polymerase, kinase, ligase, synthetase, oxidoreductase, protease inhibitors, nucleic acid binding proteins, adhesion and immunity-related proteins, cytoskeletal component-manipulating proteins, lipid biosynthesis and metabolism proteins, membrane-associated proteins, hormone-like and signaling proteins, etc. These domains are ubiquitous stretches or folds of the proteins in pathogens and allergens. Pathogenesis alleviation efforts can benefit enormously if the characteristics of these domains are known. Hence, this review catalogs and discusses the role of such pivotal domains, suggesting hypotheses for better understanding of pathogenesis at molecular level.

•

Proteins have highly conserved regions or domains across pathogens and allergens.

•

Knowledge on these critical domains can facilitate our understanding of pathogenesis mechanisms.

•

Such immune manipulation-related domains include IENR1, HNHc, HELICc, ACTIN, PROF, Robl_LC7, OmpH etc.

•

These domains are presnt in enzyme, transcription regulators, adhesion proteins, and hormones.

•

This review discusses and hypothesizes on these domains.

Pneumococcal prophages are diverse, but not without structure or history

Bacteriophages (phages) infect many bacterial species, but little is known about the diversity of phages among the pneumococcus, a leading global pathogen. The objectives of this study were to determine the prevalence, diversity and molecular epidemiology of prophages (phage DNA integrated within the bacterial genome) among pneumococci isolated over the past 90 years. Nearly 500 pneumococcal genomes were investigated and RNA sequencing was used to explore prophage gene expression. We revealed that every pneumococcal genome contained prophage DNA. 286 full-length/putatively full-length pneumococcal prophages were identified, of which 163 have not previously been reported. Full-length prophages clustered into four major groups and every group dated from the 1930–40 s onward. There was limited evidence for genes shared between prophage clusters. Prophages typically integrated in one of five different sites within the pneumococcal genome. 72% of prophages possessed the virulence genes pblA and/or pblB. Individual prophages and the host pneumococcal genetic lineage were strongly associated and some prophages persisted for many decades. RNA sequencing provided clear evidence of prophage gene expression. Overall, pneumococcal prophages were highly prevalent, demonstrated a structured population, possessed genes associated with virulence, and were expressed under experimental conditions. Pneumococcal prophages are likely to play a more important role in pneumococcal biology and evolution than previously recognised.

Use of pan-genome analysis for the identification of lineage-specific genes of Helicobacter pylori

The human bacterial pathogen Helicobacter pylori has a highly variable genome, with significant allelic and sequence diversity between isolates and even within well-characterised strains, hampering comparative genomics of H. pylori In this study, pan-genome analysis has been used to identify lineage-specific genes of H. pylori A total of 346 H. pylori genomes spanning the hpAfrica1, hpAfrica2, hpAsia2, hpEurope, hspAmerind and hspEAsia multilocus sequence typing (MLST) lineages were searched for genes specifically overrepresented or underrepresented in MLST lineages or associated with the cag pathogenicity island. The only genes overrepresented in cag-positive genomes were the cag pathogenicity island genes themselves. In contrast, a total of 125 genes were either overrepresented or underrepresented in one or more MLST lineages. Of these 125 genes, alcohol/aldehyde-reducing enzymes linked with acid resistance and production of toxic aldehydes were found to be overrepresented in African lineages. Conversely, the FecA2 ferric citrate receptor was missing from hspAmerind genomes, but present in all other lineages. This work shows the applicability of pan-genome analysis for identification of lineage-specific genes of H. pylori, facilitating further investigation to allow linkage of differential distribution of genes with disease outcome or virulence of H. pylori.

Using bacterial genomes and essential genes for the development of new antibiotics

The shrinking antibiotic development pipeline together with the global increase in antibiotic resistant infections requires that new molecules with antimicrobial activity are developed. Traditional empirical screening approaches of natural and non-natural compounds have identified the majority of antibiotics that are currently available, however this approach has produced relatively few new antibiotics over the last few decades. The vast amount of bacterial genome sequence information that has become available since the sequencing of the first bacterial genome more than 20years ago holds potential for contributing to the discovery of novel antimicrobial compounds. Comparative genomic approaches can identify genes that are highly conserved within and between bacterial species, and thus may represent genes that participate in key bacterial processes. Whole genome mutagenesis studies can also identify genes necessary for bacterial growth and survival under different environmental conditions, making them attractive targets for the development of novel inhibitory compounds. In addition, transcriptomic and proteomic approaches can be used to characterize RNA and protein levels on a cellular scale, providing information on bacterial physiology that can be applied to antibiotic target identification. Finally, bacterial genomes can be mined to identify biosynthetic pathways that produce many intrinsic antimicrobial compounds and peptides. In this review, we provide an overview of past and current efforts aimed at using bacterial genomic data in the discovery and development of novel antibiotics.

Putatively novel serotypes and the potential for reduced vaccine effectiveness: capsular locus diversity revealed among 5405 pneumococcal genomes

The pneumococcus is a leading global pathogen and a key virulence factor possessed by the majority of pneumococci is an antigenic polysaccharide capsule ('serotype'), which is encoded by the capsular (cps) locus. Approximately 100 different serotypes are known, but the extent of sequence diversity within the cps loci of individual serotypes is not well understood. Investigating serotype-specific sequence variation is crucial to the design of sequence-based serotyping methodology, understanding pneumococcal conjugate vaccine (PCV) effectiveness and the design of future PCVs. The availability of large genome datasets makes it possible to assess population-level variation among pneumococcal serotypes and in this study 5405 pneumococcal genomes were used to investigate cps locus diversity among 49 different serotypes. Pneumococci had been recovered between 1916 and 2014 from people of all ages living in 51 countries. Serotypes were deduced bioinformatically, cps locus sequences were extracted and variation was assessed within the cps locus, in the context of pneumococcal genetic lineages. Overall, cps locus sequence diversity varied markedly: low to moderate diversity was revealed among serogroups/types 1, 3, 7, 9, 11 and 22; whereas serogroups/types 6, 19, 23, 14, 15, 18, 33 and 35 displayed high diversity. Putative novel and/or hybrid cps loci were identified among all serogroups/types apart from 1, 3 and 9. This study demonstrated that cps locus sequence diversity varied widely between serogroups/types. Investigation of the biochemical structure of the polysaccharide capsule of major variants, particularly PCV-related serotypes and those that appear to be novel or hybrids, is warranted.

Increasing Metagenomic Resolution of Microbiome Interactions Through Functional Phylogenomics and Bacterial Sub-Communities

The genomic composition of the microbiome and its relationship with the environment is an exciting open question in biology. Metagenomics is a useful tool in the discovery of previously unknown taxa, but its use to understand the functional and ecological capacities of the microbiome is limited until taxonomy and function are understood in the context of the community. We suggest that this can be achieved using a combined functional phylogenomics and co-culture-based experimental strategy that can increase our capacity to measure sub-community interactions. Functional phylogenomics can identify and partition the genome such that hidden gene functions and gene clusters with unique evolutionary signals are revealed. We can test these phylogenomic predictions using an experimental model based on sub-community populations that represent a subset of the diversity directly obtained from environmental samples. These populations increase the detection of mechanisms that drive functional forces in the assembly of the microbiome, in particular the role of metabolites from key taxa in community interactions. Our combined approach leverages the potential of metagenomics to address biological questions from ecological systems.

Identifying targets for antibiotic development using omics technologies

The lack of new compounds in the antibiotic development pipeline together with the increasing incidence of infections caused by antibiotic-resistant bacteria on a global scale represents an alarming public health problem. Advances in genomic, transcriptomic and proteomic technologies permit the characterization of bacterial physiology at an unprecedented scale, and thus can facilitate the identification of bacterial factors that could serve as targets for the development of new antibiotics. Recent studies employing these technologies have permitted the elucidation of key components in multiple bacterial processes such as bacterial survival, persistence in the host and infection. The continued use of these approaches and the incorporation of emerging omics technologies hold great potential in elucidating high value targets for antibiotic development.

Targeted amplification for enhanced detection of biothreat agents by next-generation sequencing

BACKGROUND

Historically, identification of causal agents of disease has relied heavily on the ability to culture the organism in the laboratory and/or the use of pathogen-specific antibodies or sequence-based probes. However, these methods can be limiting: Even highly sensitive PCR-based assays must be continually updated due to signature degradation as new target strains and near neighbors are sequenced. Thus, there has been a need for assays that do not suffer as greatly from these limitations and/or biases. Recent advances in library preparation technologies for Next-Generation Sequencing (NGS) are focusing on the use of targeted amplification and targeted enrichment/capture to ensure that the most highly discriminating regions of the genomes of known targets (organism-unique regions and/or regions containing functionally important genes or phylogenetically-discriminating SNPs) will be sequenced, regardless of the complex sample background.

RESULTS

In the present study, we have assessed the feasibility of targeted sequence enhancement via amplification to facilitate detection of a bacterial pathogen present in low copy numbers in a background of human genomic material. Our results indicate that the targeted amplification of signature regions can effectively identify pathogen genomic material present in as little as 10 copies per ml in a complex sample. Importantly, the correct species and strain calls could be made in amplified samples, while this was not possible in unamplified samples.

CONCLUSIONS

The results presented here demonstrate the efficacy of a targeted amplification approach to biothreat detection, using multiple highly-discriminative amplicons per biothreat organism that provide redundancy in case of variation in some primer regions. Importantly, strain level discrimination was possible at levels of 10 genome equivalents. Similar results could be obtained through use of panels focused on the identification of amplicons targeted for specific genes or SNPs instead of, or in addition to, those targeted for specific organisms (ongoing gene-targeting work to be reported later). Note that without some form of targeted enhancement, the enormous background present in complex clinical and environmental samples makes it highly unlikely that sufficient coverage of key pathogen(s) present in the sample will be achieved with current NGS technology to guarantee that the most highly discriminating regions will be sequenced.

Population genomic datasets describing the post-vaccine evolutionary epidemiology of Streptococcus pneumoniae

Streptococcus pneumoniae is common nasopharyngeal commensal bacterium and important human pathogen. Vaccines against a subset of pneumococcal antigenic diversity have reduced rates of disease, without changing the frequency of asymptomatic carriage, through altering the bacterial population structure. These changes can be studied in detail through using genome sequencing to characterise systematically-sampled collections of carried S. pneumoniae. This dataset consists of 616 annotated draft genomes of isolates collected from children during routine visits to primary care physicians in Massachusetts between 2001, shortly after the seven valent polysaccharide conjugate vaccine was introduced, and 2007. Also made available are a core genome alignment and phylogeny describing the overall population structure, clusters of orthologous protein sequences, software for inferring serotype from Illumina reads, and whole genome alignments for the analysis of closely-related sets of pneumococci. These data can be used to study both bacterial evolution and the epidemiology of a pathogen population under selection from vaccine-induced immunity.

Genomic analyses of pneumococci reveal a wide diversity of bacteriocins - including pneumocyclicin, a novel circular bacteriocin

Background

One of the most important global pathogens infecting all age groups is Streptococcus pneumoniae (the ‘pneumococcus’). Pneumococci reside in the paediatric nasopharynx, where they compete for space and resources, and one competition strategy is to produce a bacteriocin (antimicrobial peptide or protein) to attack other bacteria and an immunity protein to protect against self-destruction. We analysed a collection of 336 diverse pneumococcal genomes dating from 1916 onwards, identified bacteriocin cassettes, detailed their genetic composition and sequence diversity, and evaluated the data in the context of the pneumococcal population structure.

Results

We found that all genomes maintained a blp bacteriocin cassette and we identified several novel blp cassettes and genes. The composition of the ‘bacteriocin/immunity region’ of the blp cassette was highly variable: one cassette possessed six bacteriocin genes and eight putative immunity genes, whereas another cassette had only one of each. Both widely-distributed and highly clonal blp cassettes were identified. Most surprisingly, one-third of pneumococcal genomes also possessed a cassette encoding a novel circular bacteriocin that we called pneumocyclicin, which shared a similar genetic organisation to well-characterised circular bacteriocin cassettes in other bacterial species. Pneumocyclicin cassettes were mainly of one genetic cluster and largely found among seven major pneumococcal clonal complexes.

Conclusions

These detailed genomic analyses revealed a novel pneumocyclicin cassette and a wide variety of blp bacteriocin cassettes, suggesting that competition in the nasopharynx is a complex biological phenomenon.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1729-4) contains supplementary material, which is available to authorized users.

Use of Alignment-Free Phylogenetics for Rapid Genome Sequence-Based Typing of Helicobacter pylori Virulence Markers and Antibiotic Susceptibility

Whole-genome sequencing is becoming a leading technology in the typing and epidemiology of microbial pathogens, but the increase in genomic information necessitates significant investment in bioinformatic resources and expertise, and currently used methodologies struggle with genetically heterogeneous bacteria such as the human gastric pathogen Helicobacter pylori. Here we demonstrate that the alignment-free analysis method feature frequency profiling (FFP) can be used to rapidly construct phylogenetic trees of draft bacterial genome sequences on a standard desktop computer and that coupling with in silico genotyping methods gives useful information for comparative and clinical genomic and molecular epidemiology applications. FFP-based phylogenetic trees of seven gastric Helicobacter species matched those obtained by analysis of 16S rRNA genes and ribosomal proteins, and FFP- and core genome single nucleotide polymorphism-based analysis of 63 H. pylori genomes again showed comparable phylogenetic clustering, consistent with genomotypes assigned by using multilocus sequence typing (MLST). Analysis of 377 H. pylori genomes highlighted the conservation of genomotypes and linkage with phylogeographic characteristics and predicted the presence of an incomplete or nonfunctional cag pathogenicity island in 18/276 genomes. In silico analysis of antibiotic susceptibility markers suggests that most H. pylori hspAmerind and hspEAsia isolates are predicted to carry the T₂₈₁₂C mutation potentially conferring low-level clarithromycin resistance, while levels of metronidazole resistance were similar in all multilocus sequence types. In conclusion, the use of FFP phylogenetic clustering and in silico genotyping allows determination of genome evolution and phylogeographic clustering and can contribute to clinical microbiology by genomotyping for outbreak management and the prediction of pathogenic potential and antibiotic susceptibility.

Defining and Evaluating a Core Genome Multilocus Sequence Typing Scheme for Whole-Genome Sequence-Based Typing of Listeria monocytogenes

Whole-genome sequencing (WGS) has emerged today as an ultimate typing tool to characterize Listeria monocytogenes outbreaks. However, data analysis and interlaboratory comparability of WGS data are still challenging for most public health laboratories. Therefore, we have developed and evaluated a new L. monocytogenes typing scheme based on genome-wide gene-by-gene comparisons (core genome multilocus the sequence typing [cgMLST]) to allow for a unique typing nomenclature. Initially, we determined the breadth of the L. monocytogenes population based on MLST data with a Bayesian approach. Based on the genome sequence data of representative isolates for the whole population, cgMLST target genes were defined and reappraised with 67 L. monocytogenes isolates from two outbreaks and serotype reference strains. The Bayesian population analysis generated five L. monocytogenes groups. Using all available NCBI RefSeq genomes (n = 36) and six additionally sequenced strains, all genetic groups were covered. Pairwise comparisons of these 42 genome sequences resulted in 1,701 cgMLST targets present in all 42 genomes with 100% overlap and ≥90% sequence similarity. Overall, ≥99.1% of the cgMLST targets were present in 67 outbreak and serotype reference strains, underlining the representativeness of the cgMLST scheme. Moreover, cgMLST enabled clustering of outbreak isolates with ≤10 alleles difference and unambiguous separation from unrelated outgroup isolates. In conclusion, the novel cgMLST scheme not only improves outbreak investigations but also enables, due to the availability of the automatically curated cgMLST nomenclature, interlaboratory exchange of data that are crucial, especially for rapid responses during transsectorial outbreaks.

The evolution of Campylobacter jejuni and Campylobacter coli

The global significance of Campylobacter jejuni and Campylobacter coli as gastrointestinal human pathogens has motivated numerous studies to characterize their population biology and evolution. These bacteria are a common component of the intestinal microbiota of numerous bird and mammal species and cause disease in humans, typically via consumption of contaminated meat products, especially poultry meat. Sequence-based molecular typing methods, such as multilocus sequence typing (MLST) and whole genome sequencing (WGS), have been instructive for understanding the epidemiology and evolution of these bacteria and how phenotypic variation relates to the high degree of genetic structuring in C. coli and C. jejuni populations. Here, we describe aspects of the relatively short history of coevolution between humans and pathogenic Campylobacter, by reviewing research investigating how mutation and lateral or horizontal gene transfer (LGT or HGT, respectively) interact to create the observed population structure. These genetic changes occur in a complex fitness landscape with divergent ecologies, including multiple host species, which can lead to rapid adaptation, for example, through frame-shift mutations that alter gene expression or the acquisition of novel genetic elements by HGT. Recombination is a particularly strong evolutionary force in Campylobacter, leading to the emergence of new lineages and even large-scale genome-wide interspecies introgression between C. jejuni and C. coli. The increasing availability of large genome datasets is enhancing understanding of Campylobacter evolution through the application of methods, such as genome-wide association studies, but MLST-derived clonal complex designations remain a useful method for describing population structure.

Tackling Drug Resistant Infection Outbreaks of Global Pandemic Escherichia coli ST131 Using Evolutionary and Epidemiological Genomics

High-throughput molecular screening is required to investigate the origin and diffusion of antimicrobial resistance in pathogen outbreaks. The most frequent cause of human infection is Escherichia coli, which is dominated by sequence type 131 (ST131)-a set of rapidly radiating pandemic clones. The highly infectious clades of ST131 originated firstly by a mutation enhancing conjugation and adhesion. Secondly, single-nucleotide polymorphisms occurred enabling fluoroquinolone-resistance, which is near-fixed in all ST131. Thirdly, broader resistance through beta-lactamases has been gained and lost frequently, symptomatic of conflicting environmental selective effects. This flexible approach to gene exchange is worrying and supports the proposition that ST131 will develop an even wider range of plasmid and chromosomal elements promoting antimicrobial resistance. To stop ST131, deep genome sequencing is required to understand the origin, evolution and spread of antimicrobial resistance genes. Phylogenetic methods that decipher past events can predict future patterns of virulence and transmission based on genetic signatures of adaptation and gene exchange. Both the effect of partial antimicrobial exposure and cell dormancy caused by variation in gene expression may accelerate the development of resistance. High-throughput sequencing can decode measurable evolution of cell populations within patients associated with systems-wide changes in gene expression during treatments. A multi-faceted approach can enhance assessment of antimicrobial resistance in E. coli ST131 by examining transmission dynamics between hosts to achieve a goal of pre-empting resistance before it emerges by optimising antimicrobial treatment protocols.

Genomics Reveals the Worldwide Distribution of Multidrug-Resistant Serotype 6E Pneumococci

The pneumococcus is a leading pathogen infecting children and adults. Safe, effective vaccines exist, and they work by inducing antibodies to the polysaccharide capsule (unique for each serotype) that surrounds the cell; however, current vaccines are limited by the fact that only a few of the nearly 100 antigenically distinct serotypes are included in the formulations. Within the serotypes, serogroup 6 pneumococci are a frequent cause of serious disease and common colonizers of the nasopharynx in children. Serotype 6E was first reported in 2004 but was thought to be rare; however, we and others have detected serotype 6E among recent pneumococcal collections. Therefore, we analyzed a diverse data set of ∼1,000 serogroup 6 genomes, assessed the prevalence and distribution of serotype 6E, analyzed the genetic diversity among serogroup 6 pneumococci, and investigated whether pneumococcal conjugate vaccine-induced serotype 6A and 6B antibodies mediate the killing of serotype 6E pneumococci. We found that 43% of all genomes were of serotype 6E, and they were recovered worldwide from healthy children and patients of all ages with pneumococcal disease. Four genetic lineages, three of which were multidrug resistant, described ∼90% of the serotype 6E pneumococci. Serological assays demonstrated that vaccine-induced serotype 6B antibodies were able to elicit killing of serotype 6E pneumococci. We also revealed three major genetic clusters of serotype 6A capsular sequences, discovered a new hybrid 6C/6E serotype, and identified 44 examples of serotype switching. Therefore, while vaccines appear to offer protection against serotype 6E, genetic variants may reduce vaccine efficacy in the longer term because of the emergence of serotypes that can evade vaccine-induced immunity.

Ten years of pan-genome analyses

Next generation sequencing technologies have engendered a genome sequence data deluge in public databases. Genome analyses have transitioned from single or few genomes to hundreds to thousands of genomes. Pan-genome analyses provide a framework for estimating the genomic diversity of the dataset at hand and predicting the number of additional whole genomes sequences that would be necessary to fully characterize that diversity. We review recent implementations of the pan-genome approach, its impact and limits, and we propose possible extensions, including analyses at the whole genome multiple sequence alignment level.

Comparative genomics of Lactobacillus crispatus suggests novel mechanisms for the competitive exclusion of Gardnerella vaginalis

BACKGROUND

Lactobacillus crispatus is a ubiquitous micro-organism encountered in a wide range of host-associated habitats. It can be recovered from the gastrointestinal tract of animals and it is a common constituent of the vaginal microbiota of humans. Moreover, L. crispatus can contribute to the urogenital health of the host through competitive exclusion and the production of antimicrobial agents. In order to investigate the genetic diversity of this important urogenital species, we performed a comparative genomic analysis of L. crispatus.

RESULTS

Utilizing the completed genome sequence of a strain ST1 and the draft genome sequences of nine other L. crispatus isolates, we defined the scale and scope of the pan- and core genomic potential of L. crispatus. Our comparative analysis identified 1,224 and 2,705 ortholog groups present in all or only some of the ten strains, respectively. Based on mathematical modeling, sequencing of additional L. crispatus isolates would result in the identification of new genes and functions, whereas the conserved core of the ten strains was a good representation of the final L. crispatus core genome, estimated to level at about 1,116 ortholog groups. Importantly, the current core was observed to encode bacterial components potentially promoting urogenital health. Using antibody fragments specific for one of the conserved L. crispatus adhesins, we demonstrated that the L. crispatus core proteins have a potential to reduce the ability of Gardnerella vaginalis to adhere to epithelial cells. These findings thereby suggest that L. crispatus core proteins could protect the vagina from G. vaginalis and bacterial vaginosis.

CONCLUSIONS

Our pan-genome analysis provides insights into the intraspecific genome variability and the collective molecular mechanisms of the species L. crispatus. Using this approach, we described the differences and similarities between the genomes and identified features likely to be important for urogenital health. Notably, the conserved genetic backbone of L. crispatus accounted for close to 60% of the ortholog groups of an average L. crispatus strain and included factors for the competitive exclusion of G. vaginalis, providing an explanation on how this urogenital species could improve vaginal health.