Natural genetic transformation generates a population of merodiploids in Streptococcus pneumoniae

Negative regulation of MurZ and MurA underlies the essentiality of GpsB- and StkP-mediated protein phosphorylation in Streptococcus pneumoniae D39

GpsB links peptidoglycan synthases to other proteins that determine the shape of the respiratory pathogen Streptococcus pneumoniae (pneumococcus; Spn) and other low-GC Gram-positive bacteria. GpsB is also required for phosphorylation of proteins by the essential StkP(Spn) Ser/Thr protein kinase. Here we report three classes of frequently arising chromosomal duplications (≈21-176 genes) containing murZ (MurZ-family homolog of MurA) or murA that suppress ΔgpsB or ΔstkP. These duplications arose from three different repeated sequences and demonstrate the facility of pneumococcus to modulate gene dosage of numerous genes. Overproduction of MurZ or MurA alone or overproduction of MurZ caused by ΔkhpAB mutations suppressed ΔgpsB or ΔstkP phenotypes to varying extents. ΔgpsB and ΔstkP were also suppressed by MurZ amino-acid changes distant from the active site, including one in commonly studied laboratory strains, and by truncation or deletion of the homolog of IreB(ReoM). Unlike in other Gram-positive bacteria, MurZ is predominant to MurA in pneumococcal cells. However, ΔgpsB and ΔstkP were not suppressed by ΔclpCP, which did not alter MurZ or MurA amounts. These results support a model in which regulation of MurZ and MurA activity, likely by IreB(Spn), is the only essential requirement for StkP-mediated protein phosphorylation in exponentially growing D39 pneumococcal cells.

Chromosomal Duplications of MurZ (MurA2) or MurA (MurA1), Amino Acid Substitutions in MurZ (MurA2), and Absence of KhpAB Obviate the Requirement for Protein Phosphorylation in Streptococcus pneumoniae D39

GpsB links peptidoglycan synthases to other proteins that determine the shape of the respiratory pathogen Streptococcus pneumoniae (pneumococcus; Spn ) and other low-GC Gram-positive bacteria. GpsB is also required for phosphorylation of proteins by the essential StkP( Spn ) Ser/Thr protein kinase. Here we report three classes of frequently arising chromosomal duplications (≈21-176 genes) containing murZ (MurZ-family homolog of MurA) or murA that suppress Δ gpsB or Δ stkP . These duplications arose from three different repeated sequences and demonstrate the facility of pneumococcus to modulate gene dosage of numerous genes. Overproduction of MurZ or MurA alone or overexpression of MurZ caused by Δ khpAB mutations suppressed Δ gpsB or Δ stkP phenotypes to varying extents. Δ gpsB and Δ stkP were also suppressed by MurZ amino-acid changes distant from the active site, including one in commonly studied laboratory strains, and by truncation or deletion of the homolog of IreB(ReoM). Unlike in other Gram-positive bacteria, MurZ is predominant to MurA in pneumococcal cells. However, Δ gpsB and Δ stkP were not suppressed by Δ clpCP , which did not alter MurZ or MurA amounts. These results support a model in which regulation of MurZ and MurA activity, likely by IreB( Spn ), is the only essential requirement for protein phosphorylation in exponentially growing D39 pneumococcal cells.

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.

Reconstructing Gene Gains and Losses with BadiRate

Estimating gene gain and losses is paramount to understand the molecular mechanisms underlying adaptive evolution. Despite the advent of high-throughput sequencing, such analyses have been so far hampered by the poor contiguity of genome assemblies. The increasing affordability of long-read sequencing technologies will however revolutionize our capacity to identify gene gains and losses at an unprecedented resolution, even in non-model organisms. To thoroughly exploit all such multigene family variation, the software BadiRate implements a collection of birth-and-death stochastic models, aiming at estimating by maximum likelihood the gene turnover rates along the internal and external branches of a given phylogenetic species tree. Its statistical framework also provides versatility for inferring the gene family content at the internal phylogenetic nodes (and to estimate the minimum number of gene gains and losses in each branch), for statistically contrasting competing hypotheses (e.g., accelerations of the gene turnover rates at pre-defined clades), and for pinpointing gene family expansions or contractions likely driven by natural selection. In this chapter we review the theoretical models implemented in BadiRate and illustrate their applicability by analyzing a hypothetical data set of 14 microbial species.

Tight Interplay between Replication Stress and Competence Induction in Streptococcus pneumoniae

Cells respond to genome damage by inducing restorative programs, typified by the SOS response of Escherichia coli. Streptococcus pneumoniae (the pneumococcus), with no equivalent to the SOS system, induces the genetic program of competence in response to many types of stress, including genotoxic drugs. The pneumococcal competence regulon is controlled by the origin-proximal, auto-inducible comCDE operon. It was previously proposed that replication stress induces competence through continued initiation of replication in cells with arrested forks, thereby increasing the relative comCDE gene dosage and expression and accelerating the onset of competence. We have further investigated competence induction by genome stress. We find that absence of RecA recombinase stimulates competence induction, in contrast to SOS response, and that double-strand break repair (RexB) and gap repair (RecO, RecR) initiation effectors confer a similar effect, implying that recombinational repair removes competence induction signals. Failure of replication forks provoked by titrating PolC polymerase with the base analogue HPUra, over-supplying DnaA initiator, or under-supplying DnaE polymerase or DnaC helicase stimulated competence induction. This induction was not correlated with concurrent changes in origin-proximal gene dosage. Our results point to arrested and unrepaired replication forks, rather than increased comCDE dosage, as a basic trigger of pneumococcal competence.

Reconstructing the Evolutionary History of a Highly Conserved Operon Cluster in Gammaproteobacteria and Bacilli

The evolution of gene order rearrangements within bacterial chromosomes is a fast process. Closely related species can have almost no conservation in long-range gene order. A prominent exception to this rule is a >40 kb long cluster of five core operons (secE-rpoBC-str-S10-spc-alpha) and three variable adjacent operons (cysS, tufB, and ecf) that together contain 57 genes of the transcriptional and translational machinery. Previous studies have indicated that at least part of this operon cluster might have been present in the last common ancestor of bacteria and archaea. Using 204 whole genome sequences, ∼2 Gy of evolution of the operon cluster were reconstructed back to the last common ancestors of the Gammaproteobacteria and of the Bacilli. A total of 163 independent evolutionary events were identified in which the operon cluster was altered. Further examination showed that the process of disconnecting two operons generally follows the same pattern. Initially, a small number of genes is inserted between the operons breaking the concatenation followed by a second event that fully disconnects the operons. While there is a general trend for loss of gene synteny over time, there are examples of increased alteration rates at specific branch points or within specific bacterial orders. This indicates the recurrence of relaxed selection on the gene order within bacterial chromosomes. The analysis of the alternation events indicates that segmental genome duplications and/or transposon-directed recombination play a crucial role in rearrangements of the operon cluster.

A positive perspective on DNA methylation: regulatory functions of DNA methylation outside of host defense in Gram-positive bacteria

The presence of post-replicative DNA methylation is pervasive among both prokaryotic and eukaryotic organisms. In bacteria, the study of DNA methylation has largely been in the context of restriction-modification systems, where DNA methylation serves to safeguard the chromosome against restriction endonuclease cleavage intended for invading DNA. There has been a growing recognition that the methyltransferase component of restriction-modification systems can also regulate gene expression, with important contributions to virulence factor gene expression in bacterial pathogens. Outside of restriction-modification systems, DNA methylation from orphan methyltransferases, which lack cognate restriction endonucleases, has been shown to regulate important processes, including DNA replication, DNA mismatch repair, and the regulation of gene expression. The majority of research and review articles have been focused on DNA methylation in the context of Gram-negative bacteria, with emphasis toward Escherichia coli, Caulobacter crescentus, and related Proteobacteria. Here we summarize the epigenetic functions of DNA methylation outside of host defense in Gram-positive bacteria, with a focus on the regulatory effects of both phase variable methyltransferases and DNA methyltransferases from traditional restriction-modification systems.

Bacterial Transformation Buffers Environmental Fluctuations through the Reversible Integration of Mobile Genetic Elements

Horizontal gene transfer (HGT) promotes the spread of genes within bacterial communities. Among the HGT mechanisms, natural transformation stands out as being encoded by the bacterial core genome. Natural transformation is often viewed as a way to acquire new genes and to generate genetic mixing within bacterial populations. Another recently proposed function is the curing of bacterial genomes of their infectious parasitic mobile genetic elements (MGEs). Here, we propose that these seemingly opposing theoretical points of view can be unified. Although costly for bacterial cells, MGEs can carry functions that are at points in time beneficial to bacteria under stressful conditions (e.g., antibiotic resistance genes). Using computational modeling, we show that, in stochastic environments, an intermediate transformation rate maximizes bacterial fitness by allowing the reversible integration of MGEs carrying resistance genes, although these MGEs are costly for host cell replication. Based on this dual function (MGE acquisition and removal), transformation would be a key mechanism for stabilizing the bacterial genome in the long term, and this would explain its striking conservation.IMPORTANCE Natural transformation is the acquisition, controlled by bacteria, of extracellular DNA and is one of the most common mechanisms of horizontal gene transfer, promoting the spread of resistance genes. However, its evolutionary function remains elusive, and two main roles have been proposed: (i) the new gene acquisition and genetic mixing within bacterial populations and (ii) the removal of infectious parasitic mobile genetic elements (MGEs). While the first one promotes genetic diversification, the other one promotes the removal of foreign DNA and thus genome stability, making these two functions apparently antagonistic. Using a computational model, we show that intermediate transformation rates, commonly observed in bacteria, allow the acquisition then removal of MGEs. The transient acquisition of costly MGEs with resistance genes maximizes bacterial fitness in environments with stochastic stress exposure. Thus, transformation would ensure both a strong dynamic of the bacterial genome in the short term and its long-term stabilization.

PGAP-X: extension on pan-genome analysis pipeline

Background

Since PGAP (pan-genome analysis pipeline) was published in 2012, it has been widely employed in bacterial genomics research. Though PGAP has integrated several modules for pan-genomics analysis, how to properly and effectively interpret and visualize the results data is still a challenge.

Result

To well present bacterial genomic characteristics, a novel cross-platform software was developed, named PGAP-X. Four kinds of data analysis modules were developed and integrated: whole genome sequences alignment, orthologous genes clustering, pan-genome profile analysis, and genetic variants analysis. The results from these analyses can be directly visualized in PGAP-X. The modules for data visualization in PGAP-X include: comparison of genome structure, gene distribution by conservation, pan-genome profile curve and variation on genic and genomic region. Meanwhile, result data produced by other programs with similar function can be imported to be further analyzed and visualized in PGAP-X. To test the performance of PGAP-X, we comprehensively analyzed 14 Streptococcus pneumonia strains and 14 Chlamydia trachomatis. The results show that, S. pneumonia strains have higher diversity on genome structure and gene contents than C. trachomatis strains. In addition, S. pneumonia strains might have suffered many evolutionary events, such genomic rearrangements, frequent horizontal gene transfer, homologous recombination, and other evolutionary process.

Conclusion

Briefly, PGAP-X directly presents the characteristics of bacterial genomic diversity with different visualization methods, which could help us to intuitively understand dynamics and evolution in bacterial genomes. The source code and the pre-complied executable programs are freely available from http://pgapx.ybzhao.com.

Electronic supplementary material

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

A programmed cell division delay preserves genome integrity during natural genetic transformation in Streptococcus pneumoniae

Competence for genetic transformation is a differentiation program during which exogenous DNA is imported into the cell and integrated into the chromosome. In Streptococcus pneumoniae, competence develops transiently and synchronously in all cells during exponential phase, and is accompanied by a pause in growth. Here, we reveal that this pause is linked to the cell cycle. At least two parallel pathways impair peptidoglycan synthesis in competent cells. Single-cell analyses demonstrate that ComM, a membrane protein induced during competence, inhibits both initiation of cell division and final constriction of the cytokinetic ring. Competence also interferes with the activity of the serine/threonine kinase StkP, the central regulator of pneumococcal cell division. We further present evidence that the ComM-mediated delay in division preserves genomic integrity during transformation. We propose that cell division arrest is programmed in competent pneumococcal cells to ensure that transformation is complete before resumption of cell division, to provide this pathogen with the maximum potential for genetic diversity and adaptation.

In Streptococcus pneumoniae, competence for genetic transformation is accompanied by a pause in growth. Here, Bergé et al. show that this pause is linked to the cell cycle via at least two pathways that impair peptidoglycan synthesis and preserve genomic integrity during transformation.

Genome-wide fitness analyses of the foodborne pathogen Campylobacter jejuni in in vitro and in vivo models

Campylobacter is the most common cause of foodborne bacterial illness worldwide. Faecal contamination of meat, especially chicken, during processing represents a key route of transmission to humans. There is a lack of insight into the mechanisms driving C. jejuni growth and survival within hosts and the environment. Here, we report a detailed analysis of C. jejuni fitness across models reflecting stages in its life cycle. Transposon (Tn) gene-inactivation libraries were generated in three C. jejuni strains and the impact on fitness during chicken colonisation, survival in houseflies and under nutrient-rich and -poor conditions at 4 °C and infection of human gut epithelial cells was assessed by Tn-insertion site sequencing (Tn-seq). A total of 331 homologous gene clusters were essential for fitness during in vitro growth in three C. jejuni strains, revealing that a large part of its genome is dedicated to growth. We report novel C. jejuni factors essential throughout its life cycle. Importantly, we identified genes that fulfil important roles across multiple conditions. Our comprehensive screens showed which flagella elements are essential for growth and which are vital to the interaction with host organisms. Future efforts should focus on how to exploit this knowledge to effectively control infections caused by C. jejuni.

Overcoming the Barrier of Low Efficiency during Genetic Transformation of Streptococcus mitis

Objective:Streptococcus mitis is a predominant oral colonizer, but difficulties in genetic manipulation of this species have hampered our understanding of the mechanisms it uses for colonization of oral surfaces. The aim of this study was to reveal optimal conditions for natural genetic transformation in S. mitis and illustrate its application in direct genome editing.

Methods: Luciferase reporter assays were used to assess gene expression of the alternative sigma factor (σ^X) in combination with natural transformation experiments to evaluate the efficiency by which S. mitis activates the competence system and incorporates exogenous DNA. Optimal amounts and sources of donor DNA (chromosomal, amplicon, or replicative plasmid), concentrations of synthetic competence-stimulating peptide, and transformation media were assessed.

Results: A semi-defined medium showed much improved results for response to the competence stimulating peptide when compared to rich media. The use of a donor amplicon with large homology flanking regions also provided higher transformation rates. Overall, an increase of transformation efficiencies from 0.001% or less to over 30% was achieved with the developed protocol. We further describe the construction of a markerless mutant based on this high efficiency strategy.

Conclusion: We optimized competence development in S. mitis, by use of semi-defined medium and appropriate concentrations of synthetic competence factor. Combined with the use of a large amplicon of donor DNA, this method allowed easy and direct editing of the S. mitis genome, broadening the spectrum of possible downstream applications of natural transformation in this species.

Sequence diversity of NanA manifests in distinct enzyme kinetics and inhibitor susceptibility

Streptococcus pneumoniae is the leading pathogen causing bacterial pneumonia and meningitis. Its surface-associated virulence factor neuraminidase A (NanA) promotes the bacterial colonization by removing the terminal sialyl residues from glycoconjugates on eukaryotic cell surface. The predominant role of NanA in the pathogenesis of pneumococci renders it an attractive target for therapeutic intervention. Despite the highly conserved activity of NanA, our alignment of the 11 NanAs revealed the evolutionary diversity of this enzyme. The amino acid substitutions we identified, particularly those in the lectin domain and in the insertion domain next to the catalytic centre triggered our special interest. We synthesised the representative NanAs and the mutagenized derivatives from E. coli for enzyme kinetics study and neuraminidase inhibitor susceptibility test. Via molecular docking we got a deeper insight into the differences between the two major variants of NanA and their influence on the ligand-target interactions. In addition, our molecular dynamics simulations revealed a prominent intrinsic flexibility of the linker between the active site and the insertion domain, which influences the inhibitor binding. Our findings for the first time associated the primary sequence diversity of NanA with the biochemical properties of the enzyme and with the inhibitory efficiency of neuraminidase inhibitors.

RecFOR is not required for pneumococcal transformation but together with XerS for resolution of chromosome dimers frequently formed in the process

Homologous recombination (HR) is required for both genome maintenance and generation of diversity in eukaryotes and prokaryotes. This process initiates from single-stranded (ss) DNA and is driven by a universal recombinase, which promotes strand exchange between homologous sequences. The bacterial recombinase, RecA, is loaded onto ssDNA by recombinase loaders, RecBCD and RecFOR for genome maintenance. DprA was recently proposed as a third loader dedicated to genetic transformation. Here we assessed the role of RecFOR in transformation of the human pathogen Streptococcus pneumoniae. We firstly established that RecFOR proteins are not required for plasmid transformation, strongly suggesting that DprA ensures annealing of plasmid single-strands internalized in the process. We then observed no reduction in chromosomal transformation using a PCR fragment as donor, contrasting with the 10,000-fold drop in dprA^- cells and demonstrating that RecFOR play no role in transformation. However, a ∼1.45-fold drop in transformation was observed with total chromosomal DNA in recFOR mutants. To account for this limited deficit, we hypothesized that transformation with chromosomal DNA stimulated unexpectedly high frequency (>30% of cells) formation of chromosome dimers as an intermediate in the generation of tandem duplications, and that RecFOR were crucial for dimer resolution. We validated this hypothesis, showing that the site-specific recombinase XerS was also crucial for dimer resolution. An even higher frequency of dimer formation (>80% of cells) was promoted by interspecies transformation with Streptococcus mitis chromosomal DNA, which contains numerous inversions compared to pneumococcal chromosome, each potentially promoting dimerization. In the absence of RecFOR and XerS, dimers persist, as confirmed by DAPI staining, and can limit the efficiency of transformation, since resulting in loss of transformant chromosome. These findings strengthen the view that different HR machineries exist for genome maintenance and transformation in pneumococci. These observations presumably apply to most naturally transformable species.

Homologous recombination (HR) is a widespread process which maintains genome integrity and promotes diversity. In bacteria, HR mends damaged DNA to ensure genome integrity and is also involved in transformation, a mechanism of horizontal gene transfer allowing acquisition of new genetic traits. HR is driven by recombinases, which are loaded onto single-stranded DNA by the recombinase loaders RecBCD and RecFOR for genome maintenance. DprA was recently proposed as another loader dedicated to transformation. During transformation, foreign DNA is taken up as single strands and integrated into the chromosome by HR. In this study, we show that RecFOR is not involved in transformation in Streptococcus pneumoniae. These results provide further support to the existence of different HR machineries dedicated to genetic transformation and genome maintenance in this pathogen. In addition, we show that transformation with chromosomal DNA generates chromosome dimers with unexpectedly high frequency, and that their resolution requires RecFOR and the site-specific recombinase XerS. In cells lacking these proteins, dimers persist and have a detrimental effect on the efficiency of transformation. Since the HR mechanisms leading to dimer formation are most likely conserved, this effect is presumably general to naturally transformable species.

Parallel evolution of Streptococcus pneumoniae and Streptococcus mitis to pathogenic and mutualistic lifestyles

The bacterium Streptococcus pneumoniae is one of the leading causes of fatal infections affecting humans. Intriguingly, phylogenetic analysis shows that the species constitutes one evolutionary lineage in a cluster of the otherwise commensal Streptococcus mitis strains, with which humans live in harmony. In a comparative analysis of 35 genomes, including phylogenetic analyses of all predicted genes, we have shown that the pathogenic pneumococcus has evolved into a master of genomic flexibility while lineages that evolved into the nonpathogenic S. mitis secured harmonious coexistence with their host by stabilizing an approximately 15%-reduced genome devoid of many virulence genes. Our data further provide evidence that interspecies gene transfer between S. pneumoniae and S. mitis occurs in a unidirectional manner, i.e., from S. mitis to S. pneumoniae. Import of genes from S. mitis and other mitis, anginosus, and salivarius group streptococci ensured allelic replacements and antigenic diversification and has been driving the evolution of the remarkable structural diversity of capsular polysaccharides of S. pneumoniae. Our study explains how the unique structural diversity of the pneumococcal capsule emerged and conceivably will continue to increase and reveals a striking example of the fragile border between the commensal and pathogenic lifestyles. While genomic plasticity enabling quick adaptation to environmental stress is a necessity for the pathogenic streptococci, the commensal lifestyle benefits from stability. Importance: One of the leading causes of fatal infections affecting humans, Streptococcus pneumoniae, and the commensal Streptococcus mitis are closely related obligate symbionts associated with hominids. Faced with a shortage of accessible hosts, the two opposing lifestyles evolved in parallel. We have shown that the nonpathogenic S. mitis secured harmonious coexistence with its host by stabilizing a reduced genome devoid of many virulence genes. Meanwhile, the pathogenic pneumococcus evolved into a master of genomic flexibility and imports genes from S. mitis and other related streptococci. This process ensured antigenic diversification and has been driving the evolution of the remarkable structural diversity of capsular polysaccharides of S. pneumoniae, which conceivably will continue to increase and present a challenge to disease prevention.

Streptococcus pneumoniae, le transformiste

Streptococcus pneumoniae (the pneumococcus) is an important human pathogen. Natural genetic transformation, which was discovered in this species, involves internalization of exogenous single-stranded DNA and its incorporation into the chromosome. It allows acquisition of pathogenicity islands and antibiotic resistance and promotes vaccine escape via capsule switching. This opinion article discusses how recent advances regarding several facets of pneumococcal transformation support the view that the process has evolved to maximize plasticity potential in this species, making the pneumococcus le transformiste of the bacterial kingdom and providing an advantage in the constant struggle between this pathogen and its host.