Comparative genomic analyses of seventeen Streptococcus pneumoniae strains: insights into the pneumococcal supragenome

Biological puzzles solved by using Streptococcus pneumoniae: a historical review of the pneumococcal studies that have impacted medicine and shaped molecular bacteriology

The major human pathogen Streptococcus pneumoniae has been the subject of intensive clinical and basic scientific study for over 140 years. In multiple instances, these efforts have resulted in major breakthroughs in our understanding of basic biological principles as well as fundamental tenets of bacterial pathogenesis, immunology, vaccinology, and genetics. Discoveries made with S. pneumoniae have led to multiple major public health victories that have saved the lives of millions. Studies on S. pneumoniae continue today, where this bacterium is being used to dissect the impact of the host on disease processes, as a powerful cell biology model, and to better understand the consequence of human actions on commensal bacteria at the population level. Herein we review the major findings, i.e., puzzle pieces, made with S. pneumoniae and how, over the years, they have come together to shape our understanding of this bacterium's biology and the practice of medicine and modern molecular biology.

Characterization of the family-level Borreliaceae pan-genome and development of an episomal typing protocol

Background

The Borreliaceae family includes many obligate parasitic bacterial species which are etiologically associated with a myriad of zoonotic borrelioses including Lyme disease and vector-borne relapsing fevers. Infections by the Borreliaceae are difficult to detect by both direct and indirect methods, often leading to delayed and missed diagnoses. Efforts to improve diagnoses center around the development of molecular diagnostics (MDx), but due to deep tissue sequestration of the causative spirochaetes and the lack of persistent bacteremias, even MDx assays suffer from a lack of sensitivity. Additionally, the highly extensive genomic heterogeneity among isolates, even within the same species, contributes to the lack of assay sensitivity as single target assays cannot provide universal coverage. This within-species heterogeneity is partly due to differences in replicon repertoires and genomic structures that have likely arisen to support the complex Borreliaceae lifecycle in which these parasites have to survive in multiple hosts each with unique immune responses.

Results

We constructed a Borreliaceae family-level pangenome and characterized the phylogenetic relationships among the constituent taxa which supports the recent taxonomy of splitting the family into at least two genera. Gene content pro les were created for the majority of the Borreliaceae replicons, providing for the first time their unambiguous molecular typing.

Conclusion

Our characterization of the Borreliaceae pan-genome supports the splitting of the former Borrelia genus into two genera and provides for the phylogenetic placement of several non-species designated isolates. Mining this family-level pangenome will enable precision diagnostics corresponding to gene content-driven clinical outcomes while also providing targets for interventions.

RNA cis-regulators are important for Streptococcus pneumoniae in vivo success

Bacteria have evolved complex transcriptional regulatory networks, as well as many diverse regulatory strategies at the RNA level, to enable more efficient use of metabolic resources and a rapid response to changing conditions. However, most RNA-based regulatory mechanisms are not well conserved across different bacterial species despite controlling genes important for virulence or essential biosynthetic processes. Here, we characterize the activity of, and assess the fitness benefit conferred by, twelve cis-acting regulatory RNAs (including several riboswitches and a T-box), in the opportunistic pathogen Streptococcus pneumoniae TIGR4. By evaluating native locus mutants of each regulator that result in constitutively active or repressed expression, we establish that growth defects in planktonic culture are associated with constitutive repression of gene expression, while constitutive activation of gene expression is rarely deleterious. In contrast, in mouse nasal carriage and pneumonia models, strains with either constitutively active and repressed gene expression are significantly less fit than matched control strains. Furthermore, two RNA-regulated pathways, FMN synthesis/transport and pyrimidine synthesis/transport display exceptional sensitivity to mis-regulation or constitutive gene repression in both planktonic culture and in vivo environments. Thus, despite lack of obvious phenotypes associated with constitutive gene expression in vitro, the fitness benefit conferred on bacteria via fine-tuned metabolic regulation through cis-acting regulatory RNAs is substantial in vivo, and therefore easily sufficient to drive the evolution and maintenance of diverse RNA regulatory mechanisms.

Pneumococcal Extracellular Vesicles Mediate Horizontal Gene Transfer via the Transformation Machinery

Bacterial cells secrete extracellular vesicles (EVs), the function of which is a matter of intense investigation. Here, we show that the EVs secreted by the human pathogen Streptococcus pneumoniae (pneumococcus) are associated with bacterial DNA on their surface and can deliver this DNA to the transformation machinery of competent cells. These findings suggest that EVs contribute to gene transfer in Gram-positive bacteria, and in doing so, may promote the spread of drug resistance genes in the population.

Naturally-occurring serotype 3 Streptococcus pneumoniae strains that lack functional pneumolysin and autolysin have attenuated virulence but induce localized protective immune responses

Streptococcus pneumoniae is an important cause of fatal pneumonia in humans. These bacteria express virulence factors, such as the toxins pneumolysin and autolysin, that drive host inflammatory responses. In this study we confirm loss of pneumolysin and autolysin function in a group of clonal pneumococci that have a chromosomal deletion resulting in a pneumolysin-autolysin fusion gene Δ(lytA'-ply')593. The Δ(lytA'-ply')593 pneumococci strains naturally occur in horses and infection is associated with mild clinical signs. Here we use immortalized and primary macrophage in vitro models, which include pattern recognition receptor knock-out cells, and a murine acute pneumonia model to show that a Δ(lytA'-ply')593 strain induces cytokine production by cultured macrophages, however, unlike the serotype-matched ply+lytA+ strain, it induces less tumour necrosis factor α (TNFα) and no interleukin-1β production. The TNFα induced by the Δ(lytA'-ply')593 strain requires MyD88 but, in contrast to the ply+lytA+ strain, is not reduced in cells lacking TLR2, 4 or 9. In comparison to the ply+lytA+ strain in a mouse model of acute pneumonia, infection with the Δ(lytA'-ply')593 strain resulted in less severe lung pathology, comparable levels of interleukin-1α, but minimal release of other pro-inflammatory cytokines, including interferon-γ, interleukin-6 and TNFα. These results suggest a mechanism by which a naturally occurring Δ(lytA'-ply')593 mutant strain of S. pneumoniae that resides in a non-human host has reduced inflammatory and invasive capacity compared to a human S. pneumoniae strain. These data probably explain the relatively mild clinical disease in response to S. pneumoniae infection seen in horses in comparison to humans.

Post-vaccine epidemiology of serotype 3 pneumococci identifies transformation inhibition through prophage-driven alteration of a non-coding RNA

BACKGROUND

The respiratory pathogen Streptococcus pneumoniae (the pneumococcus) is a genetically diverse bacterium associated with over 101 immunologically distinct polysaccharide capsules (serotypes). Polysaccharide conjugate vaccines (PCVs) have successfully eliminated multiple targeted serotypes, yet the mucoid serotype 3 has persisted despite its inclusion in PCV13. This capsule type is predominantly associated with a single globally disseminated strain, GPSC12 (clonal complex 180).

METHODS

A genomic epidemiology study combined previous surveillance datasets of serotype 3 pneumococci to analyse the population structure, dynamics, and differences in rates of diversification within GPSC12 during the period of PCV introductions. Transcriptomic analyses, whole genome sequencing, mutagenesis, and electron microscopy were used to characterise the phenotypic impact of loci hypothesised to affect this strain's evolution.

RESULTS

GPSC12 was split into clades by a genomic analysis. Clade I, the most common, rarely underwent transformation, but was typically infected with the prophage ϕOXC141. Prior to the introduction of PCV13, this clade's composition shifted towards a ϕOXC141-negative subpopulation in a systematically sampled UK collection. In the post-PCV13 era, more rapidly recombining non-Clade I isolates, also ϕOXC141-negative, have risen in prevalence. The low in vitro transformation efficiency of a Clade I isolate could not be fully explained by the ~100-fold reduction attributable to the serotype 3 capsule. Accordingly, prophage ϕOXC141 was found to modify csRNA3, a non-coding RNA that inhibits the induction of transformation. This alteration was identified in ~30% of all pneumococci and was particularly common in the unusually clonal serotype 1 GPSC2 strain. RNA-seq and quantitative reverse transcriptase PCR experiments using a genetically tractable pneumococcus demonstrated the altered csRNA3 was more effective at inhibiting production of the competence-stimulating peptide pheromone. This resulted in a reduction in the induction of competence for transformation.

CONCLUSION

This interference with the quorum sensing needed to induce competence reduces the risk of the prophage being deleted by homologous recombination. Hence the selfish prophage-driven alteration of a regulatory RNA limits cell-cell communication and horizontal gene transfer, complicating the interpretation of post-vaccine population dynamics.

A bacterial pan-genome makes gene essentiality strain-dependent and evolvable

Many bacterial species are represented by a pan-genome, whose genetic repertoire far outstrips that of any single bacterial genome. Here we investigate how a bacterial pan-genome might influence gene essentiality and whether essential genes that are initially critical for the survival of an organism can evolve to become non-essential. By using Transposon insertion sequencing (Tn-seq), whole-genome sequencing and RNA-seq on a set of 36 clinical Streptococcus pneumoniae strains representative of >68% of the species' pan-genome, we identify a species-wide 'essentialome' that can be subdivided into universal, core strain-specific and accessory essential genes. By employing 'forced-evolution experiments', we show that specific genetic changes allow bacteria to bypass essentiality. Moreover, by untangling several genetic mechanisms, we show that gene essentiality can be highly influenced by and/or be dependent on: (1) the composition of the accessory genome, (2) the accumulation of toxic intermediates, (3) functional redundancy, (4) efficient recycling of critical metabolites and (5) pathway rewiring. While this functional characterization underscores the evolvability potential of many essential genes, we also show that genes with differential essentiality remain important antimicrobial drug target candidates, as their inactivation almost always has a severe fitness cost in vivo.

A Comparative Analysis of the Core Proteomes within and among the Bacillus subtilis and Bacillus cereus Evolutionary Groups Reveals the Patterns of Lineage- and Species-Specific Adaptations

By integrating phylogenomic and comparative analyses of 1104 high-quality genome sequences, we identify the core proteins and the lineage-specific fingerprint proteins of the various evolutionary clusters (clades/groups/species) of the Bacillus genus. As fingerprints, we denote those core proteins of a certain lineage that are present only in that particular lineage and absent in any other Bacillus lineage. Thus, these lineage-specific fingerprints are expected to be involved in particular adaptations of that lineage. Intriguingly, with a few notable exceptions, the majority of the Bacillus species demonstrate a rather low number of species-specific fingerprints, with the majority of them being of unknown function. Therefore, species-specific adaptations are mostly attributed to highly unstable (in evolutionary terms) accessory proteomes and possibly to changes at the gene regulation level. A series of comparative analyses consistently demonstrated that the progenitor of the Cereus Clade underwent an extensive genomic expansion of chromosomal protein-coding genes. In addition, the majority (76–82%) of the B. subtilis proteins that are essential or play a significant role in sporulation have close homologs in most species of both the Subtilis and the Cereus Clades. Finally, the identification of lineage-specific fingerprints by this study may allow for the future development of highly specific vaccines, therapeutic molecules, or rapid and low-cost molecular tests for species identification.

Exploration of Streptococcus core genome to reveal druggable targets and novel therapeutics against S. pneumoniae

Streptococcus pneumoniae (S. pneumoniae), the major etiological agent of community-acquired pneumonia (CAP) contributes significantly to the global burden of infectious diseases which is getting resistant day by day. Nearly 30% of the S. pneumoniae genomes encode hypothetical proteins (HPs), and better understandings of these HPs in virulence and pathogenicity plausibly decipher new treatments. Some of the HPs are present across many Streptococcus species, systematic assessment of these unexplored HPs will disclose prospective drug targets. In this study, through a stringent bioinformatics analysis of the core genome and proteome of S. pneumoniae PCS8235, we identified and analyzed 28 HPs that are common in many Streptococcus species and might have a potential role in the virulence or pathogenesis of the bacteria. Functional annotations of the proteins were conducted based on the physicochemical properties, subcellular localization, virulence prediction, protein-protein interactions, and identification of essential genes, to find potentially druggable proteins among 28 HPs. The majority of the HPs are involved in bacterial transcription and translation. Besides, some of them were homologs of enzymes, binding proteins, transporters, and regulators. Protein-protein interactions revealed HP PCS8235_RS05845 made the highest interactions with other HPs and also has TRP structural motif along with virulent and pathogenic properties indicating it has critical cellular functions and might go under unconventional protein secretions. The second highest interacting protein HP PCS8235_RS02595 interacts with the Regulator of chromosomal segregation (RocS) which participates in chromosome segregation and nucleoid protection in S. pneumoniae. In this interacting network, 54% of protein members have virulent properties and 40% contain pathogenic properties. Among them, most of these proteins circulate in the cytoplasmic area and have hydrophilic properties. Finally, molecular docking and dynamics simulation demonstrated that the antimalarial drug Artenimol can act as a drug repurposing candidate against HP PCS8235_RS 04650 of S. pneumoniae. Hence, the present study could aid in drugs against S. pneumoniae.

Comparative Genomic Analysis of Seven Vibrio alginolyticus Strains Isolated From Shrimp Larviculture Water With Emphasis on Chitin Utilization

The opportunistic pathogen Vibrio alginolyticus is gaining attention because of its disease-causing risks to aquatic animals and humans. In this study, seven Vibrio strains isolated from different shrimp hatcheries in Southeast China were subjected to genome sequencing and subsequent comparative analysis to explore their intricate relationships with shrimp aquaculture. The seven isolates had an average nucleotide identity of ≥ 98.3% with other known V. alginolyticus strains. The species V. alginolyticus had an open pan-genome, with the addition of ≥ 161 novel genes following each new genome for seven isolates and 14 publicly available V. alginolyticus strains. The percentages of core genes of the seven strains were up to 83.1–87.5%, indicating highly conserved functions, such as chitin utilization. Further, a total of 14 core genes involved in the chitin degradation pathway were detected on the seven genomes with a single copy, 12 of which had undergone significant purifying selection (dN/dS < 1). Moreover, the seven strains could utilize chitin as the sole carbon-nitrogen source. In contrast, mobile genetic elements (MGEs) were identified in seven strains, including plasmids, prophages, and genomic islands, which mainly encoded accessory genes annotated as hypothetical proteins. The infection experiment showed that four of the seven strains might be pathogenic because the survival rates of Litopenaeus vannamei postlarvae were significantly reduced (P < 0.05) when compared to the control. However, no obvious correlation was noted between the number of putative virulence factors and toxic effects of the seven strains. Collectively, the persistence of V. alginolyticus in various aquatic environments may be attributed to its high genomic plasticity via the acquisition of novel genes by various MGEs. In view of the strong capability of chitin utilization by diverse vibrios, the timely removal of massive chitin-rich materials thoroughly in shrimp culture systems may be a key strategy to inhibit proliferation of vibrios and subsequent infection of shrimp. In addition, transcontinental transfer of potentially pathogenic V. alginolyticus strains should receive great attention to avoid vibriosis.

PclR is a transcriptional activator of the gene that encodes the pneumococcal collagen-like protein PclA

The Gram-positive bacterium Streptococcus pneumoniae is a major human pathogen that shows high levels of genetic variability. The pneumococcal R6 genome harbours several gene clusters that are not present in all strains of the species. One of these clusters contains two divergent genes, pclA, which encodes a putative surface-exposed protein that contains large regions of collagen-like repeats, and spr1404 (here named pclR). PclA was shown to mediate pneumococcal adherence to host cells in vitro. In this work, we demonstrate that PclR (494 amino acids) is a transcriptional activator. It stimulates transcription of the pclA gene by binding to a specific DNA site upstream of the core promoter. In addition, we show that PclR has common features with the MgaSpn transcriptional regulator (493 amino acids), which is also encoded by the R6 genome. These proteins have high sequence similarity (60.3%), share the same organization of predicted functional domains, and generate multimeric complexes on linear double-stranded DNAs. However, on the PpclA promoter region, MgaSpn binds to a site different from the one recognized by PclR. Our results indicate that PclR and MgaSpn have similar DNA-binding properties but different DNA-binding specificities, pointing to a different regulatory role of both proteins.

Inhibition of Streptococcus pneumoniae growth by masarimycin

Despite renewed interest, development of chemical biology methods to study peptidoglycan metabolism has lagged in comparison to the glycobiology field in general. To address this, a panel of diamides were screened against the Gram-positive bacterium Streptococcus pneumoniae to identify inhibitors of bacterial growth. The screen identified the diamide masarimycin as a bacteriostatic inhibitor of S. pneumoniae growth with an MIC of 8 µM. The diamide inhibited detergent-induced autolysis in a concentration-dependent manner, indicating perturbation of peptidoglycan degradation as the mode-of-action. Cell based screening of masarimycin against a panel of autolysin mutants, identified a higher MIC against a ΔlytB strain lacking an endo-N-acetylglucosaminidase involved in cell division. Subsequent biochemical and phenotypic analyses suggested that the higher MIC was due to an indirect interaction with LytB. Further analysis of changes to the cell surface in masarimycin treated cells identified the overexpression of several moonlighting proteins, including elongation factor Tu which is implicated in regulating cell shape. Checkerboard assays using masarimycin in concert with additional antibiotics identified an antagonistic relationship with the cell wall targeting antibiotic fosfomycin, which further supports a cell wall mode-of-action.

A Jack of All Trades: The Role of Pneumococcal Surface Protein A in the Pathogenesis of Streptococcus pneumoniae

Streptococcus pneumoniae (Spn), or the pneumococcus, is a Gram-positive bacterium that colonizes the upper airway. Spn is an opportunistic pathogen capable of life-threatening disease should it become established in the lungs, gain access to the bloodstream, or disseminate to vital organs including the central nervous system. Spn is encapsulated, allowing it to avoid phagocytosis, and current preventative measures against infection include polyvalent vaccines composed of capsular polysaccharide corresponding to its most prevalent serotypes. The pneumococcus also has a plethora of surface components that allow the bacteria to adhere to host cells, facilitate the evasion of the immune system, and obtain vital nutrients; one family of these are the choline-binding proteins (CBPs). Pneumococcal surface protein A (PspA) is one of the most abundant CBPs and confers protection against the host by inhibiting recognition by C-reactive protein and neutralizing the antimicrobial peptide lactoferricin. Recently our group has identified two new roles for PspA: binding to dying host cells via host-cell bound glyceraldehyde 3-phosphate dehydrogenase and co-opting of host lactate dehydrogenase to enhance lactate availability. These properties have been shown to influence Spn localization and enhance virulence in the lower airway, respectively. Herein, we review the impact of CBPs, and in particular PspA, on pneumococcal pathogenesis. We discuss the potential and limitations of using PspA as a conserved vaccine antigen in a conjugate vaccine formulation. PspA is a vital component of the pneumococcal virulence arsenal - therefore, understanding the molecular aspects of this protein is essential in understanding pneumococcal pathogenesis and utilizing PspA as a target for treating or preventing pneumococcal pneumonia.

An Update on "Reverse Vaccinology": The Pathway from Genomes and Epitope Predictions to Tailored, Recombinant Vaccines

In this chapter, we review the computational approaches that have led to a new generation of vaccines in recent years. There are many alternative routes to develop vaccines based on the concept of reverse vaccinology. They all follow the same basic principles-mining available genome and proteome information for antigen candidates, and recombinantly expressing them for vaccine production. Some of the same principles have been used successfully for cancer therapy approaches. In this review, we focus on infectious diseases, describing the general workflow from bioinformatic predictions of antigens and epitopes down to examples where such predictions have been used successfully for vaccine development.

Chromosomal integration of Tn5253 occurs downstream of a conserved 11-bp sequence of the rbgA gene in Streptococcus pneumoniae and in all the other known hosts of this integrative conjugative element (ICE)

Background

Tn5253, a composite Integrative Conjugative Element (ICE) of Streptococcus pneumoniae carrying tet(M) and cat resistance determinants, was found to (i) integrate at specific 83-bp integration site (attB), (ii) produce circular forms joined by a 84-bp sequence (attTn), and (iii) restore the chromosomal integration site. The purpose of this study is to functionally characterize the attB in S. pneumoniae strains with different genetic backgrounds and in other bacterial species, and to investigate the presence of Tn5253 attB site into bacterial genomes.

Results

Analysis of representative Tn5253-carryng transconjugants obtained in S. pneumoniae strains with different genetic backgrounds and in other bacterial species, namely Streptococcus agalactiae, Streptococcus gordonii, Streptococcus pyogenes, and Enterococcus faecalis showed that: (i) Tn5253 integrates in rbgA of S. pneumoniae and in orthologous rbgA genes of other bacterial species, (ii) integration occurs always downstream of a 11-bp sequence conserved among streptococcal and enterococcal hosts, (iii) length of the attB site corresponds to length of the duplication after Tn5253 integration, (iv) attB duplication restores rbgA CDS, (v) Tn5253 produced circular forms containing the attTn site at a concentration ranging between 2.0 × 10⁻⁵ to 1.2 × 10⁻² copies per chromosome depending on bacterial species and strain, (vi) reconstitution of attB sites occurred at 3.7 × 10⁻⁵ to 1.7 × 10⁻² copies per chromosome. A database search of complete microbial genomes using Tn5253 attB as a probe showed that (i) thirteen attB variants were present in the 85 complete pneumococcal genomes, (ii) in 75 pneumococcal genomes (88.3 %), the attB site was 83 or 84 nucleotides in length, while in 10 (11.7 %) it was 41 nucleotides, (iii) in other 19 bacterial species attB was located in orthologous rbgA genes and its size ranged between 17 and 84 nucleotides, (iv) the 11-bp sequence, which correspond to the last 11 nucleotides of attB sites, is conserved among the different bacterial species and can be considered the core of the Tn5253 integration site.

Conclusions

A functional characterization of the Tn5253 attB integration site combined with genome analysis contributed to elucidating the potential of Tn5253 horizontal gene transfer among different bacterial species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13100-021-00253-z.

CrfP, a fratricide protein, contributes to natural transformation in Streptococcus suis

Streptococcussuis (S.suis) is an important zoonotic pathogen that causes septicaemia, meningitis and streptococcal toxic shock-like syndrome in its host, and recent studies have shown that S.suis could be competent for natural genetic transformation. Transformation is an important mechanism for the horizontal transfer of DNA, but some elements that affect the transformation process need to be further explored. Upon entering the competent state, Streptococcus species stimulate the transcription of competence-related genes that are responsible for exogenous DNA binding, uptake and processing. In this study, we performed conserved promoter motif and qRT-PCR analyses and identified CrfP as a novel murein hydrolase that is widespread in S.suis and stimulated with a peptide pheromone in the competent state through a process controlled by ComX. A bioinformatics analysis revealed that CrfP consists of a CHAP hydrolase domain and two bacterial Src homology 3-binding (SH3b) domains. Further characterization showed that CrfP could be exported to extracellular bacterial cells and lytic S.suis strains of different serotypes, and this finding was verified by TEM and a turbidity assay. To investigate the potential effect of CrfP in vivo, a gene-deletion mutant (ΔcrfP) was constructed. Instead of stopping the natural transformation process, the inactivation of CrfP clearly reduced the effective transformation rate. Overall, these findings provide evidence showing that CrfP is important for S.suis serovar 2 competence.

Higher genome variability within metabolism genes associates with recurrent Clostridium difficile infection

BACKGROUND

Clostridium difficile (C. difficile) is a major source of healthcare-associated infection with a high risk of recurrence, attributable to many factors such as usage of antibiotics, older age and immunocompromised status of the patients. C. difficile has also a highly diverse genome, which may contribute to its high virulence. Herein we examined whether the genome conservation, measured as non-synonymous to synonymous mutations ratio (dN/dS) in core genes, presence of single genes, plasmids and prophages increased the risk of reinfection in a subset of 134 C. difficile isolates from our previous study in a singly hemato-oncology ward.

METHODS

C. difficile isolates were subjected to whole-genome sequencing (WGS) on Ion Torrent PGM sequencer. Genomes were assembled with MIRA5 and annotated with prokka and VRprofile. Logistic regression was used to asses the relationship between single gene presence and the odds of infection recurrence. DN/dS ratios were computed with codeml. Functional annotation was conducted with eggNOG-Mapper.

RESULTS

We have found that the presence of certain genes, associated with carbon metabolism and oxidative phosphorylation, increased the odds of infection recurrence. More core genes were under positive selective pressure in recurrent disease isolates - they were mostly associated with the metabolism of aminoacids. Finally, prophage elements were more prevalent in single infection isolates and plasmids did not influence the odds of recurrence.

CONCLUSIONS

Our findings suggest higher genetic plasticity in isolates causing recurrent infection, associated mainly with metabolism. On the other hand, the presence of prophages seems to reduce the isolates' virulence.

Beyond the pan-genome: current perspectives on the functional and practical outcomes of the distributed genome hypothesis

The principle of monoclonality with regard to bacterial infections was considered immutable prior to 30 years ago. This view, espoused by Koch for acute infections, has proven inadequate regarding chronic infections as persistence requires multiple forms of heterogeneity among the bacterial population. This understanding of bacterial plurality emerged from a synthesis of what-were-then novel technologies in molecular biology and imaging science. These technologies demonstrated that bacteria have complex life cycles, polymicrobial ecologies, and evolve in situ via the horizontal exchange of genic characters. Thus, there is an ongoing generation of diversity during infection that results in far more highly complex microbial communities than previously envisioned. This perspective is based on the fundamental tenet that the bacteria within an infecting population display genotypic diversity, including gene possession differences, which result from horizontal gene transfer mechanisms including transformation, conjugation, and transduction. This understanding is embodied in the concepts of the supragenome/pan-genome and the distributed genome hypothesis (DGH). These paradigms have fostered multiple researches in diverse areas of bacterial ecology including host-bacterial interactions covering the gamut of symbiotic relationships including mutualism, commensalism, and parasitism. With regard to the human host, within each of these symbiotic relationships all bacterial species possess attributes that contribute to colonization and persistence; those species/strains that are pathogenic also encode traits for invasion and metastases. Herein we provide an update on our understanding of bacterial plurality and discuss potential applications in diagnostics, therapeutics, and vaccinology based on perspectives provided by the DGH with regard to the evolution of pathogenicity.

Comparative Analysis of the Core Proteomes among the Pseudomonas Major Evolutionary Groups Reveals Species-Specific Adaptations for Pseudomonas aeruginosa and Pseudomonas chlororaphis

The Pseudomonas genus includes many species living in diverse environments and hosts. It is important to understand which are the major evolutionary groups and what are the genomic/proteomic components they have in common or are unique. Towards this goal, we analyzed 494 complete Pseudomonas proteomes and identified 297 core-orthologues. The subsequent phylogenomic analysis revealed two well-defined species (Pseudomonas aeruginosa and Pseudomonas chlororaphis) and four wider phylogenetic groups (Pseudomonas fluorescens, Pseudomonas stutzeri, Pseudomonas syringae, Pseudomonas putida) with a sufficient number of proteomes. As expected, the genus-level core proteome was highly enriched for proteins involved in metabolism, translation, and transcription. In addition, between 39–70% of the core proteins in each group had a significant presence in each of all the other groups. Group-specific core proteins were also identified, with P. aeruginosa having the highest number of these and P. fluorescens having none. We identified several P. aeruginosa-specific core proteins (such as CntL, CntM, PlcB, Acp1, MucE, SrfA, Tse1, Tsi2, Tse3, and EsrC) that are known to play an important role in its pathogenicity. Finally, a holin family bacteriocin and a mitomycin-like biosynthetic protein were found to be core-specific for P. cholororaphis and we hypothesize that these proteins may confer a competitive advantage against other root-colonizers.

Shewanella polaris sp. nov., a psychrotolerant bacterium isolated from Arctic brown algae

A Gram-stain-negative, facultatively anaerobic, flagellated and rod-shaped bacterium, designated strain SM1901T, was isolated from a brown algal sample collected from Kings Bay, Svalbard, Arctic. Strain SM1901T grew at -4‒30 °C and with 0-7.0 % (w/v) NaCl. It reduced nitrate to nitrite and hydrolysed DNA and Tween 80. Results of phylogenetic analyses based on 16S rRNA gene sequences indicated that strain SM1901T was affiliated with the genus Shewanella, showing the highest sequence similarity to the type strain of Shewanella litoralis (97.5%), followed by those of Shewanella vesiculosa, Shewanella livingstonensis and Shewanella saliphila (97.3 % for all three). The major cellular fatty acids were summed feature 3 (C16 : 1 ω7с and/or C16 : 1 ω6с), C16 : 0, C18 : 0, iso-C15 : 0 and C17 : 1 ω8с and the major polar lipids were phosphatidylethanolamine and phosphatidylglycerol. The respiratory quinones were ubiquinones Q-7, Q-8, menaquinones MK-7(H) and MK-8. The genome of strain SM1901T was 4648537 nucleotides long and encoded a variety of cold adaptation related genes, providing clues for better understanding the ecological adaptation mechanisms of polar bacteria. The genomic DNA G+C content of strain SM1901T was 40.5 mol%. Based on the polyphasic evidence presented in this paper, strain SM1901T was considered to represent a novel species, constituting a novel psychrotolerant lineage out of the known SF clade encompassed by polar Shewanella species, within the genus Shewanella, for which the name Shewanella polaris sp. nov. is proposed. The type strain is SM1901T (=KCTC 72047T=MCCC 1K03585T).

Panel 3: Genomics, precision medicine and targeted therapies

OBJECTIVE

To review the most recent advances in human and bacterial genomics as applied to pathogenesis and clinical management of otitis media.

DATA SOURCES

PubMed articles published since the last meeting in June 2015 up to June 2019.

REVIEW METHODS

A panel of experts in human and bacterial genomics of otitis media was formed. Each panel member reviewed the literature in their respective fields and wrote draft reviews. The reviews were shared with all panel members, and a merged draft was created. The panel met at the 20th International Symposium on Recent Advances in Otitis Media in June 2019, discussed the review and refined the content. A final draft was made, circulated, and approved by the panel members.

CONCLUSION

Trans-disciplinary approaches applying pan-omic technologies to identify human susceptibility to otitis media and to understand microbial population dynamics, patho-adaptation and virulence mechanisms are crucial to the development of novel, personalized therapeutics and prevention strategies for otitis media.

IMPLICATIONS FOR PRACTICE

In the future otitis media prevention strategies may be augmented by mucosal immunization, combination vaccines targeting multiple pathogens, and modulation of the middle ear microbiome. Both treatment and vaccination may be tailored to an individual's otitis media phenotype as defined by molecular profiles obtained by using rapidly developing techniques in microbial and host genomics.

Molecular Determinants of Substrate Selectivity of a Pneumococcal Rgg-Regulated Peptidase-Containing ABC Transporter

Peptidase-containing ABC transporters (PCATs) are a widely distributed family of transporters which secrete double-glycine (GG) peptides. In the opportunistic pathogen Streptococcus pneumoniae (pneumococcus), the PCATs ComAB and BlpAB have been shown to secrete quorum-sensing pheromones and bacteriocins related to the competence and pneumocin pathways. Here, we describe another pneumococcal PCAT, RtgAB, encoded by the rtg locus and found intact in 17% of strains. The Rgg/SHP-like quorum-sensing system RtgR/S, which uses a peptide pheromone with a distinctive Trp-X-Trp motif, regulates expression of the rtg locus and provides a competitive fitness advantage in a mouse model of nasopharyngeal colonization. RtgAB secretes a set of coregulated rtg GG peptides. ComAB and BlpAB, which share a substrate pool, do not secrete the rtg GG peptides. Similarly, RtgAB does not efficiently secrete ComAB/BlpAB substrates. We examined the molecular determinants of substrate selectivity between ComAB, BlpAB, and RtgAB and found that the GG peptide signal sequences contain all the information necessary to direct secretion through specific transporters. Secretion through ComAB and BlpAB depends largely on the identity of four conserved hydrophobic signal sequence residues previously implicated in substrate recognition by PCATs. In contrast, a motif situated at the N-terminal end of the signal sequence, found only in rtg GG peptides, directs secretion through RtgAB. These findings illustrate the complexity in predicting substrate-PCAT pairings by demonstrating specificity that is not dictated solely by signal sequence residues previously implicated in substrate recognition.IMPORTANCE The export of peptides from the cell is a fundamental process carried out by all bacteria. One method of bacterial peptide export relies on a family of transporters called peptidase-containing ABC transporters (PCATs). PCATs export so-called GG peptides which carry out diverse functions, including cell-to-cell communication and interbacterial competition. In this work, we describe a PCAT-encoding genetic locus, rtg, in the pathogen Streptococcus pneumoniae (pneumococcus). The rtg locus is linked to increased competitive fitness advantage in a mouse model of nasopharyngeal colonization. We also describe how the rtg PCAT preferentially secretes a set of coregulated GG peptides but not GG peptides secreted by other pneumococcal PCATs. These findings illuminate a relatively understudied part of PCAT biology: how these transporters discriminate between different subsets of GG peptides. Ultimately, expanding our knowledge of PCATs will advance our understanding of the many microbial processes dependent on these transporters.

Synergistic Activity of Mobile Genetic Element Defences in Streptococcus pneumoniae

A diverse set of mobile genetic elements (MGEs) transmit between Streptococcus pneumoniae cells, but many isolates remain uninfected. The best-characterised defences against horizontal transmission of MGEs are restriction-modification systems (RMSs), of which there are two phase-variable examples in S. pneumoniae. Additionally, the transformation machinery has been proposed to limit vertical transmission of chromosomally integrated MGEs. This work describes how these mechanisms can act in concert. Experimental data demonstrate RMS phase variation occurs at a sub-maximal rate. Simulations suggest this may be optimal if MGEs are sometimes vertically inherited, as it reduces the probability that an infected cell will switch between RMS variants while the MGE is invading the population, and thereby undermine the restriction barrier. Such vertically inherited MGEs can be deleted by transformation. The lack of between-strain transformation hotspots at known prophage att sites suggests transformation cannot remove an MGE from a strain in which it is fixed. However, simulations confirmed that transformation was nevertheless effective at preventing the spread of MGEs into a previously uninfected cell population, if a recombination barrier existed between co-colonising strains. Further simulations combining these effects of phase variable RMSs and transformation found they synergistically inhibited MGEs spreading, through limiting both vertical and horizontal transmission.

Genomics and Genetics of Streptococcus pneumoniae

Ninety years after the discovery of pneumococcal Transformation, and 74 years after the work of Avery and colleagues that identified DNA as the genetic material, Streptococcus pneumoniae is still one of the most important model organism to understand Bacterial Genetics and Genomics. In this Chapter special emphasis has been given to Genomics and to Mobile Genetic Elements (the Mobilome) which greatly contribute to the dynamic variation of pneumococcal genomes by horizontal gene transfer. Other topics include molecular mechanisms of Genetic Transformation, Restriction/Modification Systems, Mismatch DNA Repair, and techniques for construction of genetically engineered pneumococcal strains.

Comparative genome-scale modelling of the pathogenic Flavobacteriaceae species Riemerella anatipestifer in China

Riemerella anatipestifer (RA) is a gram-negative bacterium that has a high potential to infect waterfowl. Although more and more genomes of RA have been generated comparaed to genomic analysis of RA still remains at the level of individual species. In this study, we analysed the pan-genome of 27 RA virulent isolates to reveal the intraspecies genomic diversity from various aspects. The multi-locus sequence typing (MLST) analysis suggests that the geographic origin of R. anatipestifer is Guangdong province, China. Results of pan-genome analysis revealed an open pan-genome for all 27 species with the sizes of 2967 genes. We identified 387 genes among 555 unique genes originated by horizontal gene transfer. Further studies showed 204 strain-specific HGT genes were predicted as virulent proteins. Screening the 1113 core genes in RA through subtractive genomic approach, 70 putative vaccine targets out of 125 non-cytoplasmic proteins have been predicted. Further analysis of these non A. platyrhynchos homologous proteins predicted that 56 essential proteins as drug target with more interaction partners were involved in unique metabolic pathways of RA. In conclusion, the present study indicated the essence and the diversity of RA and also provides useful information for identification of vaccine and drugs candidates in future.

Orphan Genes Shared by Pathogenic Genomes Are More Associated with Bacterial Pathogenicity

Recent pangenome analyses of numerous bacterial species have suggested that each genome of a single species may have a significant fraction of its gene content unique or shared by a very few genomes (i.e., ORFans). We selected nine bacterial genera, each containing at least five pathogenic and five nonpathogenic genomes, to compare their ORFans in relation to pathogenicity-related genes. Pathogens in these genera are known to cause a number of common and devastating human diseases such as pneumonia, diphtheria, melioidosis, and tuberculosis. Thus, they are worthy of in-depth systems microbiology investigations, including the comparative study of ORFans between pathogens and nonpathogens. We provide direct evidence to suggest that ORFans shared by more pathogens are more associated with pathogenicity-related genes and thus are more important targets for development of new diagnostic markers or therapeutic drugs for bacterial infectious diseases.

Orphan genes (also known as ORFans [i.e., orphan open reading frames]) are new genes that enable an organism to adapt to its specific living environment. Our focus in this study is to compare ORFans between pathogens (P) and nonpathogens (NP) of the same genus. Using the pangenome idea, we have identified 130,169 ORFans in nine bacterial genera (505 genomes) and classified these ORFans into four groups: (i) SS-ORFans (P), which are only found in a single pathogenic genome; (ii) SS-ORFans (NP), which are only found in a single nonpathogenic genome; (iii) PS-ORFans (P), which are found in multiple pathogenic genomes; and (iv) NS-ORFans (NP), which are found in multiple nonpathogenic genomes. Within the same genus, pathogens do not always have more genes, more ORFans, or more pathogenicity-related genes (PRGs)—including prophages, pathogenicity islands (PAIs), virulence factors (VFs), and horizontal gene transfers (HGTs)—than nonpathogens. Interestingly, in pathogens of the nine genera, the percentages of PS-ORFans are consistently higher than those of SS-ORFans, which is not true in nonpathogens. Similarly, in pathogens of the nine genera, the percentages of PS-ORFans matching the four types of PRGs are also always higher than those of SS-ORFans, but this is not true in nonpathogens. All of these findings suggest the greater importance of PS-ORFans for bacterial pathogenicity.

IMPORTANCE Recent pangenome analyses of numerous bacterial species have suggested that each genome of a single species may have a significant fraction of its gene content unique or shared by a very few genomes (i.e., ORFans). We selected nine bacterial genera, each containing at least five pathogenic and five nonpathogenic genomes, to compare their ORFans in relation to pathogenicity-related genes. Pathogens in these genera are known to cause a number of common and devastating human diseases such as pneumonia, diphtheria, melioidosis, and tuberculosis. Thus, they are worthy of in-depth systems microbiology investigations, including the comparative study of ORFans between pathogens and nonpathogens. We provide direct evidence to suggest that ORFans shared by more pathogens are more associated with pathogenicity-related genes and thus are more important targets for development of new diagnostic markers or therapeutic drugs for bacterial infectious diseases.

Impact of Homologous Recombination on the Evolution of Prokaryotic Core Genomes

Microbial populations exchange genetic material through a process called homologous recombination. Although this process has been studied in particular organisms, we lack an understanding of its differential impact over the genome and across microbes with different life-styles. We used a common analytical framework to assess this process in a representative set of microorganisms. Our results uncovered important trends. First, microbes with different lifestyles are differentially impacted, with endosymbionts and obligate pathogens being those less prone to undergo this process. Second, certain genetic elements such as restriction-modification systems seem to be associated with higher rates of recombination. Most importantly, recombined genomes show the footprints of natural selection in which recombined regions preferentially contain genes that can be related to specific ecological adaptations. Taken together, our results clarify the relative contributions of factors modulating homologous recombination and show evidence for a clear a role of this process in shaping microbial genomes and driving ecological adaptations.

Homologous recombination (HR) enables the exchange of genetic material between and within species. Recent studies suggest that this process plays a major role in the microevolution of microbial genomes, contributing to core genome homogenization and to the maintenance of cohesive population structures. However, we still have a very poor understanding of the possible adaptive roles of intraspecific HR and of the factors that determine its differential impact across clades and lifestyles. Here we used a unified methodological framework to assess HR in 338 complete genomes from 54 phylogenetically diverse and representative prokaryotic species, encompassing different lifestyles and a broad phylogenetic distribution. Our results indicate that lifestyle and presence of restriction-modification (RM) machineries are among the main factors shaping HR patterns, with symbionts and intracellular pathogens having the lowest HR levels. Similarly, the size of exchanged genomic fragments correlated with the presence of RM and competence machineries. Finally, genes exchanged by HR showed functional enrichments which could be related to adaptations to different environments and ecological strategies. Taken together, our results clarify the factors underlying HR impact and suggest important adaptive roles of genes exchanged through this mechanism. Our results also revealed that the extent of genetic exchange correlated with lifestyle and some genomic features. Moreover, the genes in exchanged regions were enriched for functions that reflected specific adaptations, supporting identification of HR as one of the main evolutionary mechanisms shaping prokaryotic core genomes.

Deletion of genes involved in the ketogluconate metabolism, Entner-Doudoroff pathway, and glucose dehydrogenase increase local and invasive virulence phenotypes in Streptococcus pneumoniae

Streptococcus pneumoniae displays increased resistance to antibiotic therapy following biofilm formation. A genome-wide search revealed that SP 0320 and SP 0675 (respectively annotated as 5-keto-D-gluconate-5-reductase and glucose dehydrogenase) contain the highest degree of homology to CsgA of Myxococcus xanthus, a signaling factor that promotes cell aggregation and biofilm formation. Single and double SP 0320 and SP 0675 knockout mutants were created in strain BS72; however, no differences were observed in the biofilm-forming phenotypes of mutants compared to the wild type strain. Using the chinchilla model of otitis media and invasive disease, all three mutants exhibited greatly increased virulence compared to the wild type strain (increased pus formation, tympanic membrane rupture, mortality rates). The SP 0320 gene is located in an operon with SP 0317, SP 0318 and SP 0319, which we bioinformatically annotated as being part of the Entner-Doudoroff pathway. Deletion of SP 0317 also resulted in increased mortality in chinchillas; however, mutations in SP 0318 and SP 0319 did not alter the virulence of bacteria compared to the wild type strain. Complementing the SP 0317, SP 0320 and SP 0675 mutant strains reversed the virulence phenotype. We prepared recombinant SP 0317, SP 0318, SP 0320 and SP 0675 proteins and confirmed their functions. These data reveal that disruption of genes involved in the degradation of ketogluconate, the Entner-Doudoroff pathway, and glucose dehydrogenase significantly increase the virulence of bacteria in vivo; two hypothetical models involving virulence triggered by reduced in carbon-flux through the glycolytic pathways are presented.

Proteotyping bacteria: Characterization, differentiation and identification of pneumococcus and other species within the Mitis Group of the genus Streptococcus by tandem mass spectrometry proteomics

A range of methodologies may be used for analyzing bacteria, depending on the purpose and the level of resolution needed. The capability for recognition of species distinctions within the complex spectrum of bacterial diversity is necessary for progress in microbiological research. In clinical settings, accurate, rapid and cost-effective methods are essential for early and efficient treatment of infections. Characterization and identification of microorganisms, using, bottom-up proteomics, or "proteotyping", relies on recognition of species-unique or associated peptides, by tandem mass spectrometry analyses, dependent upon an accurate and comprehensive foundation of genome sequence data, allowing for differentiation of species, at amino acid-level resolution. In this study, the high resolution and accuracy of MS/MS-based proteotyping was demonstrated, through analyses of the three phylogenetically and taxonomically most closely-related species of the Mitis Group of the genus Streptococcus: i.e., the pathogenic species, Streptococcus pneumoniae (pneumococcus), and the commensal species, Streptococcus pseudopneumoniae and Streptococcus mitis. To achieve high accuracy, a genome sequence database used for matching peptides was created and carefully curated. Here, MS-based, bottom-up proteotyping was observed and confirmed to attain the level of resolution necessary for differentiating and identifying the most-closely related bacterial species, as demonstrated by analyses of species of the Streptococcus Mitis Group, even when S. pneumoniae were mixed with S. pseudopneumoniae and S. mitis, by matching and identifying more than 200 unique peptides for each species.

Puzzling Over the Pneumococcal Pangenome

The Gram positive bacterium Streptococcus pneumoniae (pneumococcus) is a major human pathogen. It is a common colonizer of the human host, and in the nasopharynx, sinus, and middle ear it survives as a biofilm. This mode of growth is optimal for multi-strain colonization and genetic exchange. Over the last decades, the far-reaching use of antibiotics and the widespread implementation of pneumococcal multivalent conjugate vaccines have posed considerable selective pressure on pneumococci. This scenario provides an exceptional opportunity to study the evolution of the pangenome of a clinically important bacterium, and has the potential to serve as a case study for other species. The goal of this review is to highlight key findings in the studies of pneumococcal genomic diversity and plasticity.

Rgg-Shp regulators are important for pneumococcal colonization and invasion through their effect on mannose utilization and capsule synthesis

Microbes communicate with each other by using quorum sensing (QS) systems and modulate their collective 'behavior' for in-host colonization and virulence, biofilm formation, and environmental adaptation. The recent increase in genome data availability reveals the presence of several putative QS sensing circuits in microbial pathogens, but many of these have not been functionally characterized yet, despite their possible utility as drug targets. To increase the repertoire of functionally characterized QS systems in bacteria, we studied Rgg144/Shp144 and Rgg939/Shp939, two putative QS systems in the important human pathogen Streptococcus pneumoniae. We find that both of these QS circuits are induced by short hydrophobic peptides (Shp) upon sensing sugars found in the respiratory tract, such as galactose and mannose. Microarray analyses using cultures grown on mannose and galactose revealed that the expression of a large number of genes is controlled by these QS systems, especially those encoding for essential physiological functions and virulence-related genes such as the capsular locus. Moreover, the array data revealed evidence for cross-talk between these systems. Finally, these Rgg systems play a key role in colonization and virulence, as deletion mutants of these QS systems are attenuated in the mouse models of colonization and pneumonia.

Expression of the lux genes in Streptococcus pneumoniae modulates pilus expression and virulence

Bioluminescence has been harnessed for use in bacterial reporter systems and for in vivo imaging of infection in animal models. Strain Xen35, a bioluminescent derivative of Streptococcus pneumoniae serotype 4 strain TIGR4 was previously constructed for use for in vivo imaging of infections in animal models. We have shown that strain Xen35 is less virulent than its parent TIGR4 and that this is associated with the expression of the genes for bioluminescence. The expression of the luxA-E genes in the pneumococcus reduces virulence and down regulates the expression of the pneumococcal pilus.

Genetic Competence Drives Genome Diversity in Bacillus subtilis

Prokaryote genomes are the result of a dynamic flux of genes, with increases achieved via horizontal gene transfer and reductions occurring through gene loss. The ecological and selective forces that drive this genomic flexibility vary across species. Bacillus subtilis is a naturally competent bacterium that occupies various environments, including plant-associated, soil, and marine niches, and the gut of both invertebrates and vertebrates. Here, we quantify the genomic diversity of B. subtilis and infer the genome dynamics that explain the high genetic and phenotypic diversity observed. Phylogenomic and comparative genomic analyses of 42 B. subtilis genomes uncover a remarkable genome diversity that translates into a core genome of 1,659 genes and an asymptotic pangenome growth rate of 57 new genes per new genome added. This diversity is due to a large proportion of low-frequency genes that are acquired from closely related species. We find no gene-loss bias among wild isolates, which explains why the cloud genome, 43% of the species pangenome, represents only a small proportion of each genome. We show that B. subtilis can acquire xenologous copies of core genes that propagate laterally among strains within a niche. While not excluding the contributions of other mechanisms, our results strongly suggest a process of gene acquisition that is largely driven by competence, where the long-term maintenance of acquired genes depends on local and global fitness effects. This competence-driven genomic diversity provides B. subtilis with its generalist character, enabling it to occupy a wide range of ecological niches and cycle through them.

Comparison of sequential multiplex PCR, sequetyping and whole genome sequencing for serotyping of Streptococcus pneumoniae

Streptococcus pneumoniae is one of the major causes of pneumonia, meningitis and other pneumococcal infections in young children and elders. Determination of circulating S. pneumoniae serotypes is an essential service by public health laboratories for the monitoring of putative serotype replacement following the introduction of pneumococcal conjugate vaccines (PCVs) and of the efficacy of the immunization program. The Quellung method remains the gold standard for typing S. pneumoniae. Although this method is very effective, it is also costly, time consuming and not totally reliable due to its subjective nature. The objectives of this study were to test and evaluate the efficiency of 3 different molecular methods compared to the Quellung method. Sequential multiplex PCR, sequetyping and whole genome sequencing (WGS) were chosen and tested using a set of diverse S. pneumoniae. One-hundred and eighteen isolates covering 83 serotypes were subjected to multiplex PCR and sequetyping while 88 isolates covering 53 serotypes were subjected to WGS. Sequential multiplex PCR allowed the identification of a significant proportion (49%) of serotypes at the serogroup or subset level but only 27% were identified at the serotype level. Using WGS, 55% to 60% of isolates were identified at the serotype level depending on the analysis strategy used. Finally, sequetyping demonstrated the lowest performance, with 17% of misidentified serotypes. The use of Jin cpsB database instead of the GenBank database slightly improved results but did not significantly impact the efficiency of sequetyping. Although none of these molecular methods may currently replace the Quellung method, WGS remains the most promising molecular pneumococcal serotyping method.

Probing Genomic Aspects of the Multi-Host Pathogen Clostridium perfringens Reveals Significant Pangenome Diversity, and a Diverse Array of Virulence Factors

Clostridium perfringens is an important cause of animal and human infections, however information about the genetic makeup of this pathogenic bacterium is currently limited. In this study, we sought to understand and characterise the genomic variation, pangenomic diversity, and key virulence traits of 56 C. perfringens strains which included 51 public, and 5 newly sequenced and annotated genomes using Whole Genome Sequencing. Our investigation revealed that C. perfringens has an "open" pangenome comprising 11667 genes and 12.6% of core genes, identified as the most divergent single-species Gram-positive bacterial pangenome currently reported. Our computational analyses also defined C. perfringens phylogeny (16S rRNA gene) in relation to some 25 Clostridium species, with C. baratii and C. sardiniense determined to be the closest relatives. Profiling virulence-associated factors confirmed presence of well-characterised C. perfringens-associated exotoxins genes including α-toxin (plc), enterotoxin (cpe), and Perfringolysin O (pfo or pfoA), although interestingly there did not appear to be a close correlation with encoded toxin type and disease phenotype. Furthermore, genomic analysis indicated significant horizontal gene transfer events as defined by presence of prophage genomes, and notably absence of CRISPR defence systems in >70% (40/56) of the strains. In relation to antimicrobial resistance mechanisms, tetracycline resistance genes (tet) and anti-defensins genes (mprF) were consistently detected in silico (tet: 75%; mprF: 100%). However, pre-antibiotic era strain genomes did not encode for tet, thus implying antimicrobial selective pressures in C. perfringens evolutionary history over the past 80 years. This study provides new genomic understanding of this genetically divergent multi-host bacterium, and further expands our knowledge on this medically and veterinary important pathogen.

Gene Acquisition by a Distinct Phyletic Group within Streptococcus pneumoniae Promotes Adhesion to the Ocular Epithelium

Streptococcus pneumoniae (pneumococcus) displays broad tissue tropism and infects multiple body sites in the human host. However, infections of the conjunctiva are limited to strains within a distinct phyletic group with multilocus sequence types ST448, ST344, ST1186, ST1270, and ST2315. In this study, we sequenced the genomes of six pneumococcal strains isolated from eye infections. The conjunctivitis isolates are grouped in a distinct phyletic group together with a subset of nasopharyngeal isolates. The keratitis (infection of the cornea) and endophthalmitis (infection of the vitreous body) isolates are grouped with the remainder of pneumococcal strains. Phenotypic characterization is consistent with morphological differences associated with the distinct phyletic group. Specifically, isolates from the distinct phyletic group form aggregates in planktonic cultures and chain-like structures in biofilms grown on abiotic surfaces. To begin to investigate the association between genotype and epidemiology, we focused on a predicted surface-exposed adhesin (SspB) encoded exclusively by this distinct phyletic group. Phylogenetic analysis of the gene encoding SspB in the context of a streptococcal species tree suggests that sspB was acquired by lateral gene transfer from Streptococcus suis. Furthermore, an sspB deletion mutant displays decreased adherence to cultured cells from the ocular epithelium compared to the isogenic wild-type and complemented strains. Together these findings suggest that acquisition of genes from outside the species has contributed to pneumococcal tissue tropism by enhancing the ability of a subset of strains to infect the ocular epithelium causing conjunctivitis. IMPORTANCE Changes in the gene content of pathogens can modify their ability to colonize and/or survive in different body sites in the human host. In this study, we investigate a gene acquisition event and its role in the pathogenesis of Streptococccus pneumoniae (pneumococcus). Our findings suggest that the gene encoding the predicted surface protein SspB has been transferred from Streptococcus suis (a distantly related streptococcal species) into a distinct set of pneumococcal strains. This group of strains distinguishes itself from the remainder of pneumococcal strains by extensive differences in genomic composition and by the ability to cause conjunctivitis. We find that the presence of sspB increases adherence of pneumococcus to the ocular epithelium. Thus, our data support the hypothesis that a subset of pneumococcal strains has gained genes from neighboring species that enhance their ability to colonize the epithelium of the eye, thus expanding into a new niche.

Pan-genome analysis of Bacillus for microbiome profiling

Recent advances in high-throughput sequencing technology allow for in-depth studies on microbial genomes and their communities. While multiple strains of the same species could display genomic variations with different gene contents in diverse habitats and hosts, the essential functions for a specific species are conserved as core genes that are shared among strains. We have comprehensively analyzed 238 strains of five different Bacillus species to identify the properties of core and strain-specific genes. Core and strain-specific genes in each Bacillus species show significant differences in their functions and genomic signatures. Using the core genes defined in this study, we have precisely identified the Bacillus species that exist in food microbiomes. Without resorting to culture-based whole genome sequencing, an unexpectedly large portion of the core genes, 98.22% of core genes in B. amyloliquefaciens and 97.77% of B. subtilis, were reconstructed from the microbiome. We have performed a pan-genome analysis on the core gene data of multiple Bacillus species to investigate the Bacillus species in food microbiome. Our findings provide a comprehensive genetic landscape of the Bacillus species, which is also consistent with previous studies on a limited number of strains and species. Analysis based on comprehensive core genes should thus serve as a powerful profiling tool to better understand major constituents in fermented food microbiomes.

Strain Level Streptococcus Colonization Patterns during the First Year of Life

Pneumococcal pneumonia has decreased significantly since the implementation of the pneumococcal conjugate vaccine (PCV), nevertheless, in many developing countries pneumonia mortality in infants remains high. We have undertaken a study of the nasopharyngeal (NP) microbiome during the first year of life in infants from The Philippines and South Africa. The study entailed the determination of the Streptococcus sp. carriage using a lytA qPCR assay, whole metagenomic sequencing, and in silico serotyping of Streptococcus pneumoniae, as well as 16S rRNA amplicon based community profiling. The lytA carriage in both populations increased with infant age and lytA+ samples ranged from 24 to 85% of the samples at each sampling time point. We next developed informatic tools for determining Streptococcus community composition and pneumococcal serotype from metagenomic sequences derived from a subset of longitudinal lytA-positive Streptococcus enrichment cultures from The Philippines (n = 26 infants, 50% vaccinated) and South African (n = 7 infants, 100% vaccinated). NP samples from infants were passaged in enrichment media, and metagenomic DNA was purified and sequenced. In silico capsular serotyping of these 51 metagenomic assemblies assigned known serotypes in 28 samples, and the co-occurrence of serotypes in 5 samples. Eighteen samples were not typeable using known serotypes but did encode for capsule biosynthetic cluster genes similar to non-encapsulated reference sequences. In addition, we performed metagenomic assembly and 16S rRNA amplicon profiling to understand co-colonization dynamics of Streptococcus sp. and other NP genera, revealing the presence of multiple Streptococcus species as well as potential respiratory pathogens in healthy infants. A range of virulence and drug resistant elements were identified as circulating in the NP microbiomes of these infants. This study revealed the frequent co-occurrence of multiple S. pneumoniae strains along with Streptococcus sp. and other potential pathogens such as S. aureus in the NP microbiome of these infants. In addition, the in silico serotype analysis proved powerful in determining the serotypes in S. pneumoniae carriage, and may lead to developing better targeted vaccines to prevent invasive pneumococcal disease (IPD) in these countries. These findings suggest that NP colonization by S. pneumoniae during the first years of life is a dynamic process involving multiple serotypes and species.

Exploring the Genomic Diversity and Cariogenic Differences of Streptococcus mutans Strains Through Pan-Genome and Comparative Genome Analysis

Pan-genome refers to the sum of genes that can be found in a given bacterial species, including the core-genome and the dispensable genome. In this study, the genomes from 183 Streptococcus mutans (S. mutans) isolates were analyzed from the pan-genome perspective. This analysis revealed that S. mutans has an "open" pan-genome, implying that there are plenty of new genes to be found as more genomes are sequenced. Additionally, S. mutans has a limited core-genome, which is composed of genes related to vital activities within the bacterium, such as metabolism and hereditary information storage or processing, occupying 35.6 and 26.6% of the core genes, respectively. We estimate the theoretical core-genome size to be about 1083 genes, which are fewer than other Streptococcus species. In addition, core genes suffer larger selection pressures in comparison to those that are less widely distributed. Not surprisingly, the distribution of putative virulence genes in S. mutans strains does not correlate with caries status, indicating that other factors are also responsible for cariogenesis. These results contribute to a more understanding of the evolutionary characteristics and dynamic changes within the genome components of the species. This also helps to form a new theoretical foundation for preventing dental caries. Furthermore, this study sets an example for analyzing large genomic datasets of pathogens from the pan-genome perspective.

A novel streptococcal cell-cell communication peptide promotes pneumococcal virulence and biofilm formation

Streptococcus pneumoniae (pneumococcus) is a major human pathogen. It is a common colonizer of the human respiratory track, where it utilizes cell-cell communication systems to coordinate population-level behaviors. We reasoned that secreted peptides that are highly expressed during infection are pivotal for virulence. Thus, we used in silico pattern searches to define a pneumococcal secretome and analyzed the transcriptome of the clinically important PMEN1 lineage to identify which peptide-encoding genes are highly expressed in vivo. In this study, we characterized virulence peptide 1 (vp1), a highly expressed Gly-Gly peptide-encoding gene in chinchilla middle ear effusions. The vp1 gene is widely distributed across pneumococcus as well as encoded in related species. Studies in the chinchilla model of middle ear infection demonstrated that VP1 is a virulence determinant. The vp1 gene is positively regulated by a transcription factor from the Rgg family and its cognate SHP (short hydrophobic peptide). In vitro data indicated that VP1 promotes increased thickness and biomass for biofilms grown on chinchilla middle ear epithelial cells. Furthermore, the wild-type biofilm is restored with the exogenous addition of synthetic VP1. We conclude that VP1 is a novel streptococcal regulatory peptide that controls biofilm development and pneumococcal pathogenesis.

Proteomic variation and diversity in clinical Streptococcus pneumoniae isolates from invasive and non-invasive sites

Streptococcus pneumoniae is responsible for a variety of invasive and non-invasive human infections. There are over 90 serotypes of S. pneumoniae differing in their ability to adapt to the different niches within the host. Two-dimensional gel electrophoresis was used to discriminate clinical S. pneumoniae isolates recovered from either blood cultures (invasive site isolates) or other sites, including sputum, tracheal aspirate, ear, eye and skin swabs (non-invasive site isolates). Global protein expression profiles for five invasive site and six non-invasive site isolates representing five different serotypes (serotypes 4, 6, 9, 14 and 23) were obtained for each isolate and combined into a single data set using Progenesis SameSpots^™ software. One-hundred and eighty six protein spots (39% of the protein spots in the dataset) differed significantly (ANOVA, p<0.05) in abundance between the invasive site (101 upregulated protein spots) and non-invasive site (85 upregulated protein spots) isolates. Correlations between the bacterial proteomes and their sites of isolation were determined by Principal Component Analysis (PCA) using the significantly different protein spots. Out of the 186 variable protein spots, 105 exhibited a serotype-associated pattern of variability. The expression of the remaining 81 protein spots was concluded to be uniquely linked to the site of bacterial isolation. Mass spectrometry was used to identify selected protein spots that showed either constant or differential abundance levels. The identified proteins had a diverse range of functions including, capsule biogenesis, DNA repair, protein deglycation, translation, stress response and virulence as well as amino acid, carbohydrate, lipid and nucleotide metabolism. These findings provide insight on the proteins that contribute towards the adaptation of the bacteria to different sites within the host.

Promiscuous signaling by a regulatory system unique to the pandemic PMEN1 pneumococcal lineage

Streptococcus pneumoniae (pneumococcus) is a leading cause of death and disease in children and elderly. Genetic variability among isolates from this species is high. These differences, often the product of gene loss or gene acquisition via horizontal gene transfer, can endow strains with new molecular pathways, diverse phenotypes, and ecological advantages. PMEN1 is a widespread and multidrug-resistant pneumococcal lineage. Using comparative genomics we have determined that a regulator-peptide signal transduction system, TprA2/PhrA2, was acquired by a PMEN1 ancestor and is encoded by the vast majority of strains in this lineage. We show that TprA2 is a negative regulator of a PMEN1-specific gene encoding a lanthionine-containing peptide (lcpA). The activity of TprA2 is modulated by its cognate peptide, PhrA2. Expression of phrA2 is density-dependent and its C-terminus relieves TprA2-mediated inhibition leading to expression of lcpA. In the pneumococcal mouse model with intranasal inoculation, TprA2 had no effect on nasopharyngeal colonization but was associated with decreased lung disease via its control of lcpA levels. Furthermore, the TprA2/PhrA2 system has integrated into the pneumococcal regulatory circuitry, as PhrA2 activates TprA/PhrA, a second regulator-peptide signal transduction system widespread among pneumococci. Extracellular PhrA2 can release TprA-mediated inhibition, activating expression of TprA-repressed genes in both PMEN1 cells as well as another pneumococcal lineage. Acquisition of TprA2/PhrA2 has provided PMEN1 isolates with a mechanism to promote commensalism over dissemination and control inter-strain gene regulation.

Genome Analysis of Clostridium difficile PCR Ribotype 014 Lineage in Australian Pigs and Humans Reveals a Diverse Genetic Repertoire and Signatures of Long-Range Interspecies Transmission

Clostridium difficile PCR ribotype (RT) 014 is well-established in both human and porcine populations in Australia, raising the possibility that C. difficile infection (CDI) may have a zoonotic or foodborne etiology. Here, whole genome sequencing and high-resolution core genome phylogenetics were performed on a contemporaneous collection of 40 Australian RT014 isolates of human and porcine origin. Phylogenies based on MLST (7 loci, STs 2, 13, and 49) and core orthologous genes (1260 loci) showed clustering of human and porcine strains indicative of very recent shared ancestry. Core genome single nucleotide variant (SNV) analysis found 42% of human strains showed a clonal relationship (separated by ≤2 SNVs in their core genome) with one or more porcine strains, consistent with recent inter-host transmission. Clones were spread over a vast geographic area with 50% of the human cases occurring without recent healthcare exposure. These findings suggest a persistent community reservoir with long-range dissemination, potentially due to agricultural recycling of piggery effluent. We also provide the first pan-genome analysis for this lineage, characterizing its resistome, prophage content, and in silico virulence potential. The RT014 is defined by a large "open" pan-genome (7587 genes) comprising a core genome of 2296 genes (30.3% of the total gene repertoire) and an accessory genome of 5291 genes. Antimicrobial resistance genotypes and phenotypes varied across host populations and ST lineages and were characterized by resistance to tetracycline [tetM, tetA(P), tetB(P) and tetW], clindamycin/erythromycin (ermB), and aminoglycosides (aph3-III-Sat4A-ant6-Ia). Resistance was mediated by clinically important mobile genetic elements, most notably Tn6194 (harboring ermB) and a novel variant of Tn5397 (harboring tetM). Numerous clinically important prophages (Siphoviridae and Myoviridae) were identified as well as an uncommon accessory gene regulator locus (agr3). Conservation in the pathogenicity locus and S-layer correlated with ST affiliation, further extending the concept of clonal C. difficile lineages. This study provides novel insights on the genetic variability and strain relatedness of C. difficile RT014, a lineage of emerging One Health importance. Ongoing molecular and genomic surveillance of strains in humans, animals, food, and the environment is imperative to identify opportunities to reduce the overall CDI burden.

Strain Dependent Genetic Networks for Antibiotic-Sensitivity in a Bacterial Pathogen with a Large Pan-Genome

The interaction between an antibiotic and bacterium is not merely restricted to the drug and its direct target, rather antibiotic induced stress seems to resonate through the bacterium, creating selective pressures that drive the emergence of adaptive mutations not only in the direct target, but in genes involved in many different fundamental processes as well. Surprisingly, it has been shown that adaptive mutations do not necessarily have the same effect in all species, indicating that the genetic background influences how phenotypes are manifested. However, to what extent the genetic background affects the manner in which a bacterium experiences antibiotic stress, and how this stress is processed is unclear. Here we employ the genome-wide tool Tn-Seq to construct daptomycin-sensitivity profiles for two strains of the bacterial pathogen Streptococcus pneumoniae. Remarkably, over half of the genes that are important for dealing with antibiotic-induced stress in one strain are dispensable in another. By confirming over 100 genotype-phenotype relationships, probing potassium-loss, employing genetic interaction mapping as well as temporal gene-expression experiments we reveal genome-wide conditionally important/essential genes, we discover roles for genes with unknown function, and uncover parts of the antibiotic's mode-of-action. Moreover, by mapping the underlying genomic network for two query genes we encounter little conservation in network connectivity between strains as well as profound differences in regulatory relationships. Our approach uniquely enables genome-wide fitness comparisons across strains, facilitating the discovery that antibiotic responses are complex events that can vary widely between strains, which suggests that in some cases the emergence of resistance could be strain specific and at least for species with a large pan-genome less predictable.

Bacterial volatiles and diagnosis of respiratory infections

Breath testing has enormous potential in the medical diagnostic field. The underlying complexity and perceived availability of adequate specimens, combined with a lack of knowledge of the metabolic pathways that give rise to compounds that are sources of analytes detectable in breath, has greatly slowed development. These real obstacles have recently been largely overcome in the use of breath testing to identify patients with cystic fibrosis associated Pseudomonas aeruginosa infection and tuberculosis. This review summarizes progress made in the characterization of microbial volatiles produced by major lower respiratory tract bacterial pathogens, and their potential use as diagnostic markers in patient breath testing.

Comparative Genomics of Carriage and Disease Isolates of Streptococcus pneumoniae Serotype 22F Reveals Lineage-Specific Divergence and Niche Adaptation

Streptococcus pneumoniae is a major cause of meningitis, sepsis, and pneumonia worldwide. Pneumococcal conjugate vaccines have been part of the United Kingdom’s childhood immunization program since 2006 and have significantly reduced the incidence of disease due to vaccine efficacy in reducing carriage in the population. Here we isolated two clones of 22F (an emerging serotype of clinical concern, multilocus sequence types 433 and 698) and conducted comparative genomic analysis on four isolates, paired by Sequence Type (ST) with one of each pair being derived from carriage and the other disease (sepsis). The most compelling observation was of nonsynonymous mutations in pgdA, encoding peptidoglycan N-acetylglucosamine deacetylase A, which was found in the carriage isolates of both ST433 and 698. Deacetylation of pneumococcal peptidoglycan is known to enable resistance to lysozyme upon invasion. Althought no other clear genotypic signatures related to disease or carriage could be determined, additional intriguing comparisons between the two STs were possible. These include the presence of an intact prophage, in addition to numerous additional phage insertions, within the carriage isolate of ST433. Contrasting gene repertoires related to virulence and colonization, including bacteriocins, lantibiotics, and toxin-–antitoxin systems, were also observed.

DiiA is a novel dimorphic cell wall protein of Streptococcus pneumoniae involved in invasive disease

OBJECTIVES

Many outer multidomain proteins play fundamental virulent roles in an allele-dependent manner. We aimed to investigate the influence of the outer SP1992 protein, here renamed DiiA (Dimorphic invasion-involved A), in pneumococcal disease.

METHODS

The presence and type of diiA allele was screened by PCR in 560 clinical isolates. Isogenic mutants carrying progressive diiA deletions were constructed and checked in mouse models of infection. DiiA binding to human molecules was carried out by surface plasmon resonance.

RESULTS

The diiA gene is exclusive of Streptococcus pneumoniae and included in the core genome. DiiA variants contain one or two imperfect repeats (R1 and R2), an unstructured region and a cell-wall anchor domain. Clonal complexes carrying both repeats were associated with invasive disease, while those carrying R2 preferentially caused non-invasive syndromes in patients with underlying risk factors. Mutants lacking both repeats were less efficient in nasopharyngeal colonization and dissemination from lungs. Moreover, the ΔdiiA defective strain suffered a severe impairment in bacterial proliferation in blood. Purified DiiA bound to collagen and lactoferrin with high affinity.

CONCLUSIONS

DiiA is a distinctive pneumococcal virulence factor contributing to colonization and long-term invasion in this pathogen.

BPGA- an ultra-fast pan-genome analysis pipeline

Recent advances in ultra-high-throughput sequencing technology and metagenomics have led to a paradigm shift in microbial genomics from few genome comparisons to large-scale pan-genome studies at different scales of phylogenetic resolution. Pan-genome studies provide a framework for estimating the genomic diversity of the dataset, determining core (conserved), accessory (dispensable) and unique (strain-specific) gene pool of a species, tracing horizontal gene-flux across strains and providing insight into species evolution. The existing pan genome software tools suffer from various limitations like limited datasets, difficult installation/requirements, inadequate functional features etc. Here we present an ultra-fast computational pipeline BPGA (Bacterial Pan Genome Analysis tool) with seven functional modules. In addition to the routine pan genome analyses, BPGA introduces a number of novel features for downstream analyses like core/pan/MLST (Multi Locus Sequence Typing) phylogeny, exclusive presence/absence of genes in specific strains, subset analysis, atypical G + C content analysis and KEGG & COG mapping of core, accessory and unique genes. Other notable features include minimum running prerequisites, freedom to select the gene clustering method, ultra-fast execution, user friendly command line interface and high-quality graphics outputs. The performance of BPGA has been evaluated using a dataset of complete genome sequences of 28 Streptococcus pyogenes strains.