Population biology of Gram-positive pathogens: high-risk clones for dissemination of antibiotic resistance

Infections caused by multiresistant Gram-positive bacteria represent a major health burden in the community as well as in hospitalized patients. Staphylococcus aureus, Enterococcus faecalis and Enterococcus faecium are well-known pathogens of hospitalized patients, frequently linked with resistance against multiple antibiotics, compromising effective therapy. Streptococcus pneumoniae and Streptococcus pyogenes are important pathogens in the community and S. aureus has recently emerged as an important community-acquired pathogen. Population genetic studies reveal that recombination prevails as a driving force of genetic diversity in E. faecium, E. faecalis, S. pneumoniae and S. pyogenes, and thus, these species are weakly clonal. Although recombination has a relatively modest role driving the genetic variation of the core genome of S. aureus, the horizontal acquisition of resistance and virulence genes plays a key role in the emergence of new clinically relevant clones in this species. In this review, we discuss the population genetics of E. faecium, E. faecalis, S. pneumoniae, S. pyogenes and S. aureus. Knowledge of the population structure of these pathogens is not only highly relevant for (molecular) epidemiological research but also for identifying the genetic variation that underlies changes in clinical behaviour, to improve our understanding of the pathogenic behaviour of particular clones and to identify novel targets for vaccines or immunotherapy.     

Phylogeographic variation in recombination rates within a global clone of methicillin-resistant Staphylococcus aureus

Background

Next-generation sequencing (NGS) is a powerful tool for understanding both patterns of descent over time and space (phylogeography) and the molecular processes underpinning genome divergence in pathogenic bacteria. Here, we describe a synthesis between these perspectives by employing a recently developed Bayesian approach, BRATNextGen, for detecting recombination on an expanded NGS dataset of the globally disseminated methicillin-resistant Staphylococcus aureus (MRSA) clone ST239.

Results

The data confirm strong geographical clustering at continental, national and city scales and demonstrate that the rate of recombination varies significantly between phylogeographic sub-groups representing independent introductions from Europe. These differences are most striking when mobile non-core genes are included, but remain apparent even when only considering the stable core genome. The monophyletic ST239 sub-group corresponding to isolates from South America shows heightened recombination, the sub-group predominantly from Asia shows an intermediate level, and a very low level of recombination is noted in a third sub-group representing a large collection from Turkey.

Conclusions

We show that the rapid global dissemination of a single pathogenic bacterial clone results in local variation in measured recombination rates. Possible explanatory variables include the size and time since emergence of each defined sub-population (as determined by the sampling frame), variation in transmission dynamics due to host movement, and changes in the bacterial genome affecting the propensity for recombination.     

Evidence for clonal evolution among highly polymorphic genes in methicillin-resistant Staphylococcus aureus

The evolution of Staphylococcus aureus has been described as predominantly clonal, based on evidence from seven housekeeping genes. We aimed to test if this was also true for more polymorphic genes. In a collection of 60 isolates including major European epidemic methicillin-resistant S. aureus (MRSA) and sporadic MRSA strains, we compared the partial gene sequences of seven housekeeping genes (arcC, aroE, glpF, gmk, pta, tpi, and yqiL), six core adhesion genes (present in all strains) (clfA, clfB, fnbA, map, sdrC, and spa), and four accessory adhesion genes (not present in all strains) (ebpS, fnbB, sdrD, and sdrE). Nucleotide diversity of adhesion genes was 2- to 10-fold higher than genes used for multilocus sequence typing. All genes showed evidence for purifying selection with a weakly reduced level among accessory adhesion genes. Among these highly variable genes, there was no evidence for a difference in molecular evolution between epidemic and sporadic strains. Gene trees constructed from concatenated sequences of housekeeping, core adhesion, and accessory adhesion genes were highly congruent, indicating clonality, despite some evidence for homologous exchange. Further evidence for clonality was found with an overall positive correlation of allelic and nucleotidic divergence for both seven housekeeping genes and six core adhesion genes. However, for small allelic differences that fit the demarcations of clonal complexes (CCs) there was no such correlation, suggesting that recombination occurred. Therefore, despite an overall clonal population structure, recombination between related isolates within CCs might have contributed to S. aureus evolution.     

Identification of the third type of PVL phage in ST59 methicillin-resistant Staphylococcus aureus (MRSA) strains

The genes lukS-PV and lukF-PV for Panton-Valentine leukocidin (PVL) that confers high virulence to Staphylococcus aureus are located on the prophages (PVL phages) which have been classified into group 1 and 2 sfi21-like Siphoviridae. We report novel PVL phages lysogenized in ST59 methicillin-resistant Staphylococcus aureus (MRSA) strains isolated in Japan (JCSC7247) and Taiwan (JCSC5967). The genomes of φ7247PVL and φ5967PVL showed more than 99% identity, and the regions containing the five genes located at both ends of the prophages, int (integrase), hol (holin), ami (amidase), lukS-PV, and lukF-PV, are highly homologous to extant PVL phages. The genes for the structural module are less homologous to these phages, but are highly homologous to non-PVL phages belonging to group 3 Sfi21-like Siphoviridae, for example φN315. Subsequent PCR identification and nucleotide sequencing of an additional 11 Taiwanese ST59 MRSA isolates suggested they all carry the same phage as φ5967PVL, which differed from φ7247PVL by a single base. This study adds evidence to the notion that novel PVL phages would be generated through illegitimate recombination events by acquiring the region at which hol, ami, luk, and int genes would line up upon lytic growth, and suggests that the PVL-positive MRSA clones that have emerged worldwide may carry distinct phages.     

Single-cell genomics reveal low recombination frequencies in freshwater bacteria of the SAR11 clade

Background

The SAR11 group of Alphaproteobacteria is highly abundant in the oceans. It contains a recently diverged freshwater clade, which offers the opportunity to compare adaptations to salt- and freshwaters in a monophyletic bacterial group. However, there are no cultivated members of the freshwater SAR11 group and no genomes have been sequenced yet.

Results

We isolated ten single SAR11 cells from three freshwater lakes and sequenced and assembled their genomes. A phylogeny based on 57 proteins indicates that the cells are organized into distinct microclusters. We show that the freshwater genomes have evolved primarily by the accumulation of nucleotide substitutions and that they have among the lowest ratio of recombination to mutation estimated for bacteria. In contrast, members of the marine SAR11 clade have one of the highest ratios. Additional metagenome reads from six lakes confirm low recombination frequencies for the genome overall and reveal lake-specific variations in microcluster abundances. We identify hypervariable regions with gene contents broadly similar to those in the hypervariable regions of the marine isolates, containing genes putatively coding for cell surface molecules.

Conclusions

We conclude that recombination rates differ dramatically in phylogenetic sister groups of the SAR11 clade adapted to freshwater and marine ecosystems. The results suggest that the transition from marine to freshwater systems has purged diversity and resulted in reduced opportunities for recombination with divergent members of the clade. The low recombination frequencies of the LD12 clade resemble the low genetic divergence of host-restricted pathogens that have recently shifted to a new host.     

Arginine catabolic mobile element encoded speG abrogates the unique hypersensitivity of Staphylococcus aureus to exogenous polyamines

Polyamines, including spermine (Spm) and spermidine (Spd), are aliphatic cations that are reportedly synthesized by all living organisms. They exert pleiotropic effects on cells and are required for efficient nucleic acid and protein synthesis. Here, we report that the human pathogen Staphylococcus aureus lacks identifiable polyamine biosynthetic genes, and consequently produces no Spm/Spd or their precursor compounds putrescine and agmatine. Moreover, while supplementing defined medium with polyamines generally enhances bacterial growth, Spm and Spd exert bactericidal effects on S. aureus at physiological concentrations. Small colony variants specifically lacking menaquinone biosynthesis arose after prolonged Spm exposure and exhibited reduced polyamine sensitivity. However, other respiratory-defective mutants were no less susceptible to Spm implying menaquinone itself rather than general respiration is required for full Spm toxicity. Polyamine hypersensitivity distinguishes S. aureus from other bacteria and is exhibited by all tested strains save those belonging to the USA-300 group of community-associated methicillin-resistant S. aureus (CA-MRSA). We identified one gene within the USA-300-specific arginine catabolic mobile element (ACME) encoding a Spm/Spd N-acetyltransferase that is necessary and sufficient for polyamine resistance. S. aureus encounters significant polyamine levels during infection; however, the acquisition of ACME encoded speG allows USA-300 clones to circumvent polyamine hypersensitivity, a peculiar trait of S. aureus.     

The effect of sequence divergence on recombination between direct repeats in Arabidopsis

It is well established that sequence divergence has an inhibitory effect on homologous recombination. However, a detailed analysis of this relationship is missing for most higher eukaryotes. We have measured the rate of somatic recombination between direct repeats as a function of the number, type, and position of divergent nucleotides in Arabidopsis. We show that a minor divergence level of 0.16% (one mutation in otherwise identical 618 bp) has a profound effect, decreasing the recombination rate approximately threefold. A further increase in the divergence level affects the recombination rate to a smaller extent until a "divergence saturation" effect is reached at relatively low levels of divergence ( approximately 0.5%). The type of mismatched nucleotide does not affect recombination rates. The decrease in the rate of recombination caused by a single mismatch was not affected by the position of the mismatch along the repeat. This suggests that most recombination intermediate tracts contain a mismatch and thus are as long as the full length of the 618-bp repeats. Finally, we could deduce an antirecombination efficiency of approximately 66% for the first mismatch in the repeat. Altogether, this work shows some degree of conservation across kingdoms when compared to previous reports in yeast; it also provides new insight into the effect of sequence divergence on homologous recombination.     

The ICESt1 element of Streptococcus thermophilus belongs to a large family of integrative and conjugative elements that exchange modules and change their specificity of integration

The 34,734-bp element ICESt1 from Streptococcus thermophilus CNRZ368 is site-specifically integrated into the 3(') end of the gene fda. ICESt1 encodes integrative functions and putative transfer functions. Six proteins of the putative conjugative system of ICESt1 are related to those encoded by the conjugative transposon Tn916 from Enterococcus faecalis. A comparison of these proteins with those encoded by the complete or partial genome sequences of various low G+C bacteria including Bacillus subtilis, Clostridium difficile, E. faecalis, Listeria monocytogenes, Staphylococcus aureus, and Streptococcus mutans revealed the presence of numerous putative site-specific integrative conjugative elements and/or conjugative transposons within these genomes. Sequence comparisons revealed that these elements possess a modular structure and that exchanges of unrelated or distantly related modules and genes have occurred between these elements, and also plasmids and prophages. These exchanges have probably led to modifications in the site specificity of integration of these elements. Therefore, a distinction between low specificity integrative conjugative elements (i.e., conjugative transposons) and site-specific integrative conjugative elements does not appear to be relevant. We propose to call all the conjugative elements that excise by site-specific recombination and integrate by recombination between a specific site of a circular intermediate and another site, "Integrative and Conjugative Elements" (ICEs), irrespective of the integration specificity.     

Nucleotide substitution and recombination at orthologous loci in Staphylococcus aureus

The pattern of nucleotide substitution was examined at 2,129 orthologous loci among five genomes of Staphylococcus aureus, which included two sister pairs of closely related genomes (MW2/MSSA476 and Mu50/N315) and the more distantly related MRSA252. A total of 108 loci were unusual in lacking any synonymous differences among the five genomes; most of these were short genes encoding proteins highly conserved at the amino acid sequence level (including many ribosomal proteins) or unknown predicted genes. In contrast, 45 genes were identified that showed anomalously high divergence at synonymous sites. The latter genes were evidently introduced by homologous recombination from distantly related genomes, and in many cases, the pattern of nucleotide substitution made it possible to reconstruct the most probable recombination event involved. These recombination events introduced genes encoding proteins that differed in amino acid sequence and thus potentially in function. Several of the proteins are known or likely to be involved in pathogenesis (e.g., staphylocoagulase, exotoxin, Ser-Asp fibrinogen-binding bone sialoprotein-binding protein, fibrinogen and keratin-10 binding surface-anchored protein, fibrinogen-binding protein ClfA, and enterotoxin P). Therefore, the results support the hypothesis that exchange of homologous genes among S. aureus genomes can play a role in the evolution of pathogenesis in this species.     

Intraspecies comparison of Streptomyces pratensis genomes reveals high levels of recombination and gene conservation between strains of disparate geographic origin

Background

Streptomyces are widespread bacteria that contribute to the terrestrial carbon cycle and produce the majority of clinically useful antibiotics. While interspecific genomic diversity has been investigated among Streptomyces, information is lacking on intraspecific genomic diversity. Streptomyces pratensis has high rates of homologous recombination but the impact of such gene exchange on genome evolution and the evolution of natural product gene clusters remains uncharacterized.

Results

We report draft genome sequences of four S. pratensis strains and compare to the complete genome of Streptomyces flavogriseus IAF-45-CD (=ATCC 33331), a strain recently reclassified to S. pratensis. Despite disparate geographic origins, the genomes are highly similar with 85.9% of genes present in the core genome and conservation of all natural product gene clusters. Natural products include a novel combination of carbapenem and beta-lactamase inhibitor gene clusters. While high intraspecies recombination rates abolish the phylogenetic signal across the genome, intraspecies recombination is suppressed in two genomic regions. The first region is centered on an insertion/deletion polymorphism and the second on a hybrid NRPS-PKS gene. Finally, two gene families accounted for over 25% of the divergent genes in the core genome. The first includes homologs of bldB (required for spore development and antibiotic production) while the second includes homologs of an uncharacterized protein with a helix-turn-helix motif (hpb). Genes from these families co-occur with fifteen pairs spread across the genome. These genes have evidence for co-evolution of co-localized pairs, supporting previous assertions that these genes may function akin to a toxin-antitoxin system.

Conclusions

S. pratensis genomes are highly similar with exceptional levels of recombination which erase phylogenetic signal among strains of the species. This species has a large core genome and variable terminal regions that are smaller than those found in interspecies comparisons. There is no geographic differentiation between these strains, but there is evidence for local linkage disequilibrium affecting two genomic regions. We have also shown further observational evidence that the DUF397-HTH (bldB and hpb) are a novel toxin-antitoxin pair.     

Molecular evolution of key genes for type II secretion in Legionella pneumophila

Given the role of type II protein secretion system (T2S) in the ecology and pathogenesis of Legionella pneumophila, it is possible that this system is a target for adaptive evolution. The population genetic structure of L.pneumophila was inferred from the partial sequences of rpoB and from the complete sequence of three T2S structural components (lspD, lspE and pilD) and from two T2S effectors critical for intracellular infection of protozoa (proA and srnA) of 37 strains isolated from natural and man-made environments and disease-related from worldwide sources. A phylogenetic analysis was obtained for the concatenated alignment and for each individual locus. Seven main groups were identified containing the same L.pneumophila strains, suggesting an ancient divergence for each cluster and indicating that the operating selective pressures have equally affected the evolution of the five genes. Although linkage disequilibrium analysis indicate a clonal nature for population structure in this sample, our results indicate that recombination is a common phenomenon among T2S-related genes on this species, as 24 of the 37 L.pneumophila isolates contained at least one locus in which recombination was identified. Furthermore, neutral selection acting on the analysed T2S-related genes emerged as a clear result, namely on T2S effectors, ProA and SrnA, indicating that they are probably implicated in conserved virulence mechanisms through legionellae hosts.     

Gene flow-dependent genomic divergence between Anopheles gambiae M and S forms

Anopheles gambiae sensu stricto exists as two often-sympatric races termed the M and S molecular forms, characterized by fixed differences at an X-linked marker. Extreme divergence between M and S forms at pericentromeric "genomic islands" suggested that selection on variants therein could be driving interform divergence in the presence of ongoing gene flow, but recent work has detected much more widespread genomic differentiation. Whether such genomic islands are important in reproductive isolation or represent ancestral differentiation preserved by low recombination is currently unclear. A critical test of these competing hypotheses could be provided by comparing genomic divergence when rates of recent introgression vary. We genotyped 871 single nucleotide polymorphisms (SNPs) in A. gambiae sensu stricto from locations of M and S sympatry and allopatry, encompassing the full range of observed hybridization rates (0-25%). M and S forms were readily partitioned based on genomewide SNP variation in spite of evidence for ongoing introgression that qualitatively reflects hybridization rates. Yet both the level and the heterogeneity of genomic divergence varied markedly in line with levels of introgression. A few genomic regions of differentiation between M and S were common to each sampling location, the most pronounced being two centromere-proximal speciation islands identified previously but with at least one additional region outside of areas expected to exhibit reduced recombination. Our results demonstrate that extreme divergence at genomic islands does not simply represent segregating ancestral polymorphism in regions of low recombination and can be resilient to substantial gene flow. This highlights the potential for islands comprising a relatively small fraction of the genome to play an important role in early-stage speciation when reproductive isolation is limited.     

The Tetracycline Resistance Genetet(M) Exhibits Mosaic Structure

Tetracycline resistance genes of the M class,tet(M), are typically found on mobile genetic elements as the conjugative transposons of gram-positive bacteria. By comparing the sequences of eight differenttet(M) genes (fromEnterococcus faecalis, Streptococcus pneumoniae, Staphylococcus aureus, Ureaplasma urealyticum,andNeisseria), a mosaic structure was detected which could be traced to two distinct alleles. The two alleles displayed a divergence of 8% and a different G/C content. The block structure of these genes provides evidence for the contribution of homologous recombination to the evolution and the heterogeneity of thetet(M) locus. Unlike described cases of chromosomally located mosaic loci,tet(M) is a relatively recently acquired determinant in the species examined and it would appear that mosaic structure withintet(M) has evolved after acquisition of the gene by the mobile genetic elements upon which it is located.     

The Staphylococcus aureus RNome and its commitment to virulence

Staphylococcus aureus is a major human pathogen causing a wide spectrum of nosocomial and community-associated infections with high morbidity and mortality. S. aureus generates a large number of virulence factors whose timing and expression levels are precisely tuned by regulatory proteins and RNAs. The aptitude of bacteria to use RNAs to rapidly modify gene expression, including virulence factors in response to stress or environmental changes, and to survive in a host is an evolving concept. Here, we focus on the recently inventoried S. aureus regulatory RNAs, with emphasis on those with identified functions, two of which are directly involved in pathogenicity.     

Microarrays reveal that each of the ten dominant lineages of Staphylococcus aureus has a unique combination of surface-associated and regulatory genes

Staphylococcus aureus is the most common cause of hospital-acquired infection. In healthy hosts outside of the health care setting, S. aureus is a frequent colonizer of the human nose but rarely causes severe invasive infection such as bacteremia, endocarditis, or osteomyelitis. To identify genes associated with community-acquired invasive isolates, regions of genomic variability, and the S. aureus population structure, we compared 61 community-acquired invasive isolates of S. aureus and 100 nasal carriage isolates from healthy donors using a microarray spotted with PCR products representing every gene from the seven S. aureus sequencing projects. The core genes common to all strains were identified, and 10 dominant lineages of S. aureus were clearly discriminated. Each lineage carried a unique combination of hundreds of "core variable" (CV) genes scattered throughout the chromosome, suggesting a common ancestor but early evolutionary divergence. Many CV genes are regulators of virulence genes or known or predicted to be expressed on the bacterial surface and to interact with the host during nasal colonization and infection. Within each lineage, isolates showed substantial variation in the carriage of mobile genetic elements and their associated virulence and resistance genes, indicating frequent horizontal transfer. However, we were unable to identify any association between lineage or gene and invasive isolates. We suggest that the S. aureus gene combinations necessary for invasive disease may also be necessary for nasal colonization and that community-acquired invasive disease is strongly dependent on host factors.     

Rapid divergence of two classes of Haemophilus ducreyi

Haemophilus ducreyi, the etiologic agent of chancroid, expresses variants of several key virulence factors. While previous reports suggested that H. ducreyi strains formed two clonal populations, the differences between, and diversity within, these populations were unclear. To assess their variability, we examined sequence diversity at 11 H. ducreyi loci, including virulence and housekeeping genes, augmenting published data sets with PCR-amplified genes to acquire data for at least 10 strains at each locus. While sequences from all 11 loci place strains into two distinct groups, there was very little variation within each group. The difference between alleles of the two groups was variable and large at 3 loci encoding surface-exposed proteins (0.4 < K(S) < 1.3, where K(S) is divergence at synonymous sites) but consistently small at genes encoding cytoplasmic or periplasmic proteins (K(S) < 0.09). The data suggest that the two classes have recently diverged, that recombination has introduced variant alleles into at least 3 distinct loci, and that these alleles have been confined to one of the two classes. In addition, recombination is evident among alleles within, but not between, classes. Rather than clones of the same species, these properties indicate that the two classes may form distinct species.