Sensitive, specific polymorphism discovery in bacteria using massively parallel sequencing

Validation and Development of an Escherichia coli Riboflavin Pathway Phenotypic Screen Hit as a Small-Molecule Ligand of the Flavin Mononucleotide Riboswitch

A riboflavin biosynthesis pathway-specific phenotypic screen using a library of compounds, all with unspecified antibiotic activity, identified one small molecule later named ribocil, for which intrinsic antibacterial activity against Escherichia coli was completely suppressed by addition of exogenous riboflavin to the bacterial growth medium. The ability of riboflavin to suppress the activity of ribocil, and further demonstration that ribocil inhibited riboflavin synthesis (IC₅₀ = 0.3 μM), supported that a component of the riboflavin synthesis pathway was the molecular target. Remarkably, resistance mutation selection and whole-genome sequencing showed that the target of ribocil was not an enzyme in the riboflavin biosynthesis pathway, but instead the flavin mononucleotide riboswitch, a noncoding structural RNA element in the ribB gene that encodes a key riboflavin synthesis enzyme. Although ribocil is structurally distinct from the natural riboswitch regulatory ligand flavin mononucleotide, ribocil binding to the riboswitch results in efficient repression of ribB expression and inhibition of riboflavin biosynthesis and bacterial growth. A cell-based riboswitch regulated gene reporter assay as well as an in vitro riboswitch RNA aptamer-binding assay, both of which are described in detail here along with the riboflavin pathway-specific screen, were developed to further validate the mechanism of action of ribocil and to facilitate the discovery of more potent analogues.

Canada-Wide Epidemic of emm74 Group A Streptococcus Invasive Disease

Background

The number of invasive group A Streptococcus (iGAS) infections due to hitherto extremely rare type emm74 strains has increased in several Canadian provinces since late 2015. We hypothesized that the cases recorded in the different provinces are linked and caused by strains of an emm74 clone that recently emerged and expanded explosively.

Methods

We analyzed both active and passive surveillance data for iGAS infections and used whole-genome sequencing to investigate the phylogenetic relationships of the emm74 strains responsible for these invasive infections country-wide.

Results

Genome analysis showed that highly clonal emm74 strains, genetically different from emm74 organisms previously circulating in Canada, were responsible for a country-wide epidemic of >160 invasive disease cases. The emerging clone belonged to multilocus sequence typing ST120. The analysis also revealed dissemination patterns of emm74 subclonal lineages across Canadian provinces. Clinical data analysis indicated that the emm74 epidemic disproportionally affected middle-aged or older male individuals. Homelessness, alcohol abuse, and intravenous drug usage were significantly associated with invasive emm74 infections.

Conclusions

In a period of 20 months, an emm74 GAS clone emerged and rapidly spread across several Canadian provinces located more than 4500 km apart, causing invasive infections primarily among disadvantaged persons.

Capsular Switching and Other Large-Scale Recombination Events in Invasive Sequence Type 1 Group B Streptococcus

We report several cases of recombination events leading to capsular switching among sequence type (ST) 1 group B Streptococcus strains. These strains otherwise shared a common genome backbone with serotype V ST1 strains. However, the genomes of ST1 serotype V strains and those of serotypes VI, VII, and VIII strains differed substantially.

Differential daptomycin resistance development in Staphylococcus aureus strains with active and mutated gra regulatory systems

The first-in-class lipopeptide antibiotic daptomycin (DAP) is highly active against Gram-positive pathogens including ß-lactam and glycopeptide resistant strains. Its molecular mode of action remains enigmatic, since a defined target has not been identified so far and multiple effects, primarily on the cell envelope have been observed. Reduced DAP susceptibility has been described in S. aureus and enterococci after prolonged treatment courses. In line with its pleiotropic antibiotic activities, a unique, defined molecular mechanism of resistance has not emerged, instead non-susceptibility appears often accompanied by alterations in membrane composition and changes in cell wall homeostasis. We compared S. aureus strains HG001 and SG511, which differ primarily in the functionality of the histidine kinase GraS, to evaluate the impact of the GraRS regulatory system on the development of DAP non-susceptibility. After extensive serial passing, both DAP^R variants reached a minimal inhibitory concentration of 31 μg/ml and shared some phenotypic characteristics (e.g. thicker cell wall, reduced autolysis). However, based on comprehensive analysis of the underlying genetic, transcriptomic and proteomic changes, we found that both strains took different routes to achieve DAP resistance. Our study highlights the impressive genetic and physiological capacity of S. aureus to counteract pleiotropic activities of cell wall- and membrane-active compounds even when a major cell wall regulatory system is dysfunctional.

Whole-genome sequencing of mutants with increased resistance against the two-peptide bacteriocin plantaricin JK reveals a putative receptor and potential docking site

By whole-genome sequencing of resistant mutants, a putative receptor for plantaricin JK, a two-peptide bacteriocin produced by some Lactobacillus plantarum strains, was identified in Lactobacillus plantarum NCFB 965 and Weissella viridescens NCFB 1655. The receptors of the two species had 66% identical amino acid sequences and belong to the amino acid-polyamine-organocation (APC) transporter protein family. The resistant mutants contained point mutations in the protein-encoding gene resulting in either premature stop codons, leading to truncated versions of the protein, or single amino acid substitutions. The secondary structure of the W. viridescens protein was predicted to contain 12 transmembrane (TM) helices, a core structure shared by most members of the APC protein family. The single amino acid substitutions that resulted in resistant strains were located in a confined region of the protein that consists of TM helix 10, which is predicted to be part of an inner membrane pore, and an extracellular loop between TM helix 11 and 12. By use of template-based modeling a 3D structure model of the protein was obtained, which visualizes this mutational hotspot region and further strengthen the hypothesis that it represents a docking site for plantaricin JK.

Comparative pathogenomic characterization of a non-invasive serotype M71 strain Streptococcus pyogenes NS53 reveals incongruent phenotypic implications from distinct genotypic markers

The strains serotyped as M71 from group A Streptococcus are common causes of pharyngeal and skin diseases worldwide. Here we characterize the genome of a unique non-invasive M71 human isolate, NS53. The genome does not contain structural rearrangements or large-scale gene gains/losses, but encodes a full set of non-truncated known virulence factors, thus providing an ideal reference for comparative studies. However, the NS53 genome showed incongruent phenotypic implications from distinct genotypic markers. NS53 is characterized as an emm pattern D and FCT (fibronectin-collagen-T antigen) type-3 strain, typical of skin tropic strains, but is phylogenetically close to emm pattern E strains with preference for both skin and pharyngeal infections. We propose that this incongruence could result from recombination within the emm gene locus, or, alternatively, selection has been against those genetic alterations. Combined with the inability to select for CovS switching, a process is indicated whereby NS53 has been pre-adapted to specific host niches selecting against variations in CovS and many other genes. This may allow the strain to attain successful colonization and long-term survival. A balance between genetic variations and fitness may exist for this bacterium to form a stabilized genome optimized for survival in specific host environments.

Antibacterial small molecules targeting the conserved TOPRIM domain of DNA gyrase

To combat the threat of antibiotic-resistant Gram-negative bacteria, novel agents that circumvent established resistance mechanisms are urgently needed. Our approach was to focus first on identifying bioactive small molecules followed by chemical lead prioritization and target identification. Within this annotated library of bioactives, we identified a small molecule with activity against efflux-deficient Escherichia coli and other sensitized Gram-negatives. Further studies suggested that this compound inhibited DNA replication and selection for resistance identified mutations in a subunit of E. coli DNA gyrase, a type II topoisomerase. Our initial compound demonstrated weak inhibition of DNA gyrase activity while optimized compounds demonstrated significantly improved inhibition of E. coli and Pseudomonas aeruginosa DNA gyrase and caused cleaved complex stabilization, a hallmark of certain bactericidal DNA gyrase inhibitors. Amino acid substitutions conferring resistance to this new class of DNA gyrase inhibitors reside exclusively in the TOPRIM domain of GyrB and are not associated with resistance to the fluoroquinolones, suggesting a novel binding site for a gyrase inhibitor.

Affinity Selection-Mass Spectrometry Identifies a Novel Antibacterial RNA Polymerase Inhibitor

The growing prevalence of drug resistant bacteria is a significant global threat to human health. The antibacterial drug rifampin, which functions by inhibiting bacterial RNA polymerase (RNAP), is an important part of the antibacterial armamentarium. Here, in order to identify novel inhibitors of bacterial RNAP, we used affinity-selection mass spectrometry to screen a chemical library for compounds that bind to Escherichia coli RNAP. We identified a novel small molecule, MRL-436, that binds to RNAP, inhibits RNAP, and exhibits antibacterial activity. MRL-436 binds to RNAP through a binding site that differs from the rifampin binding site, inhibits rifampin-resistant RNAP derivatives, and exhibits antibacterial activity against rifampin-resistant strains. Isolation of mutants resistant to the antibacterial activity of MRL-436 yields a missense mutation in codon 622 of the rpoC gene encoding the RNAP β' subunit or a null mutation in the rpoZ gene encoding the RNAP ω subunit, confirming that RNAP is the functional cellular target for the antibacterial activity of MRL-436, and indicating that RNAP β' subunit residue 622 and the RNAP ω subunit are required for the antibacterial activity of MRL-436. Similarity between the resistance determinant for MRL-436 and the resistance determinant for the cellular alarmone ppGpp suggests a possible similarity in binding site and/or induced conformational state for MRL-436 and ppGpp.

Serotype Distribution, Population Structure, and Antimicrobial Resistance of Group B Streptococcus Strains Recovered from Colonized Pregnant Women

Using serotyping, multilocus sequence typing, and whole-genome sequencing (WGS) of selected strains, we studied the population structure of 102 group B Streptococcus (GBS) isolates prospectively sampled in 2014 from vaginal/rectal swabs of healthy pregnant women in metropolitan Toronto, Canada. We also determined the susceptibilities of each of the colonizing isolates to penicillin, erythromycin, clindamycin, tetracycline, and other antimicrobial agents. Overall, we observed a high rate of tetracycline resistance (89%) among colonizing GBS isolates. We found resistance to erythromycin in 36% of the strains, and 33% were constitutively or inducibly resistant to clindamycin. The most frequently identified serotypes were III (25%), Ia (23%), and V (19%). Serotype IV accounted for 6% of the colonizing isolates, a rate consistent with that observed among patients with invasive GBS infections in metropolitan Toronto. The majority of serotype IV isolates belonged to sequence type (ST)459, a tetracycline-, erythromycin-, and clindamycin-resistant ST first identified in Minnesota, which is considered to be the main driver of serotype IV GBS expansion in North America. WGS revealed that ST459 isolates from Canada are clonally related to colonizing and invasive ST459 organisms circulating in regions of the United States. We also used WGS to study recombination in selected colonizing strains from metropolitan Toronto, which revealed multiple episodes of capsular switching. Present and future circulating GBS organisms and their genetic diversity may influence GBS vaccine development.

Phenotypic differentiation of Streptococcus pyogenes populations is induced by recombination-driven gene-specific sweeps

Genomic recombination plays an important role in driving adaptive evolution and population differentiation in bacteria. However, controversy exists as to the effects of recombination on population diversity and differentiation, i.e., recombination is frequent enough to sweep through the population at selected gene loci (gene-specific sweeps), or the recombination rate is low without interfering genome-wide selective sweeps. Observations supporting either view are sparse. Pathogenic bacteria causing infectious diseases are promising candidates to provide observations of recombination. However, phenotype-associated differentiations are usually vague among them due to diverse disease manifestations. Here we report a population genomic study of the group A Streptococcus pyogenes (GAS), a human pathogen with highly recombining genomes. By employing a genome-wide association study on single nucleotide polymorphisms (SNPs), we demonstrate a phenotypic differentiation of GAS, represented by separate clustering of two sublineages associated with niche-specific infections, i.e., skin infection and pharyngitis-induced acute rheumatic fever. By quantifying SNPs associated with the differentiation in a statistical and phylogenetic context, we propose that the phenotype-associated differentiation arose through recombination-driven gene-specific sweeps, rather than genome-wide sweeps. Our work provides a novel paradigm of phenotype-associated differentiation induced by gene-specific sweeps in a human pathogen and has implications for understanding of driving forces of bacterial evolution.

Capsular Switching and Other Large-Scale Recombination Events in Invasive Sequence Type 1 Group B Streptococcus

We report several cases of recombination events leading to capsular switching among sequence type (ST) 1 group B Streptococcus strains. These strains otherwise shared a common genome backbone with serotype V ST1 strains. However, the genomes of ST1 serotype V strains and those of serotypes VI, VII, and VIII strains differed substantially.

Genomic Recombination Leading to Decreased Virulence of Group B Streptococcus in a Mouse Model of Adult Invasive Disease

Adult invasive disease caused by Group B Streptococcus (GBS) is increasing worldwide. Whole-genome sequencing (WGS) now permits rapid identification of recombination events, a phenomenon that occurs frequently in GBS. Using WGS, we described that strain NGBS375, a capsular serotype V GBS isolate of sequence type (ST)297, has an ST1 genomic background but has acquired approximately 300 kbp of genetic material likely from an ST17 strain. Here, we examined the virulence of this strain in an in vivo model of GBS adult invasive infection. The mosaic ST297 strain showed intermediate virulence, causing significantly less systemic infection and reduced mortality than a more virulent, serotype V ST1 isolate. Bacteremia induced by the ST297 strain was similar to that induced by a serotype III ST17 strain, which was the least virulent under the conditions tested. Yet, under normalized bacteremia levels, the in vivo intrinsic capacity to induce the production of pro-inflammatory cytokines was similar between the ST297 strain and the virulent ST1 strain. Thus, the diminished virulence of the mosaic strain may be due to reduced capacity to disseminate or multiply in blood during a systemic infection which could be mediated by regulatory factors contained in the recombined region.

Fluoroquinolone Resistance among Clonal Complex 1 Group B Streptococcus Strains

Fluoroquinolone resistance in group B Streptococcus is increasingly being reported worldwide. Here, we correlated fluoroquinolone resistance with mutations in gyrA, gyrB, parC, and parE genes, identified by mining whole-genome sequencing (WGS) data of 190 clonal complex 1 group B Streptococcus strains recovered from patients with invasive diseases in North America. We report a high prevalence of fluoroquinolone resistance (12%) among GBS strains in our collection. Our approach is the first step towards accurate prediction of fluoroquinolone resistance from WGS data in this opportunistic pathogen.

Genomic Characterization of a Pattern D Streptococcus pyogenes emm53 Isolate Reveals a Genetic Rationale for Invasive Skin Tropicity

The genome of an invasive skin-tropic strain (AP53) of serotype M53 group A Streptococcus pyogenes (GAS) is composed of a circular chromosome of 1,860,554 bp and carries genetic markers for infection at skin locales, viz, emm gene family pattern D and FCT type 3. Through genome-scale comparisons of AP53 with other GAS genomes, we identified 596 candidate single-nucleotide polymorphisms (SNPs) that reveal a potential genetic basis for skin tropism. The genome of AP53 differed by ∼30 point mutations from a noninvasive pattern D serotype M53 strain (Alab49), 4 of which are located in virulence genes. One pseudogene, yielding an inactive sensor kinase (CovS(-)) of the two-component transcriptional regulator CovRS, a major determinant for invasiveness, severely attenuated the expression of the secreted cysteine protease SpeB and enhanced the expression of the hyaluronic acid capsule compared to the isogenic noninvasive AP53/CovS(+) strain. The collagen-binding protein transcript sclB differed in the number of 5'-pentanucleotide repeats in the signal peptides of AP53 and Alab49 (9 versus 15), translating into different lengths of their signal peptides, which nonetheless maintained a full-length translatable coding frame. Furthermore, GAS strain AP53 acquired two phages that are absent in Alab49. One such phage (ΦAP53.2) contains the known virulence factor superantigen exotoxin gene tandem speK-slaA Overall, we conclude that this bacterium has evolved in multiple ways, including mutational variations of regulatory genes, short-tandem-repeat polymorphisms, large-scale genomic alterations, and acquisition of phages, all of which may be involved in shaping the adaptation of GAS in specific infectious environments and contribute to its enhanced virulence.

IMPORTANCE

Infectious strains of S. pyogenes (GAS) are classified by their serotypes, relating to the surface M protein, the emm-like subfamily pattern, and their tropicity toward the nasopharynx and/or skin. It is generally agreed that M proteins from pattern D strains, which also directly bind human host plasminogen, are skin tropic. We have sequenced and characterized the genome of an invasive pattern D GAS strain (AP53) in comparison to a very similar strain (Alab49) that is noninvasive and developed a genomic rationale as to possible reasons for the skin tropicity of these two strains and the greater invasiveness of AP53.

A putative amino acid transporter determines sensitivity to the two-peptide bacteriocin plantaricin JK

Lactobacillus plantarum produces a number of antimicrobial peptides (bacteriocins) that mostly target closely related bacteria. Although bacteriocins are important for the ecology of these bacteria, very little is known about how the peptides target sensitive cells. In this work, a putative membrane protein receptor of the two-peptide bacteriocin plantaricin JK was identified by comparing Illumina sequence reads from plantaricin JK-resistant mutants to a crude assembly of the sensitive wild-type Weissella viridescens genome using the polymorphism discovery tool VAAL. Ten resistant mutants harbored altogether seven independent mutations in a gene encoding an APC superfamily protein with 12 transmembrane helices. The APC superfamily transporter thus is likely to serve as a target for plantaricin JK on sensitive cells.

Population structure and drug resistance patterns of emerging non-PCV-13 Streptococcus pneumoniae serotypes 22F, 15A, and 8 isolated from adults in Ontario, Canada

The introduction of pneumococcal conjugate vaccines has led to the emergence of non-vaccine serotypes, which contributed to invasive pneumococcal disease in Canada and worldwide. A significant increase in the prevalence of non-13-valent pneumococcal conjugate vaccine (PCV-13)-included serotypes 22F, 15A, and 8 was observed from 2009 to 2013 in Ontario (all p values<0.01). In this study, whole genome sequencing was conducted on the 25 isolates of serotype 22F, seven of 15A and 10 of 8 to investigate the population structure and antibiotic resistance. All seven serotype 15A isolates were found to be multidrug resistant. From whole genome analysis, we observed recombination events among serotypes 22F, 15A and 8 populations. Serotype 22F (ST433) has emerged into two sub-populations, with 28% (7/25) exhibiting recombination events, and five also acquiring macrolide resistance as a result of recombination. This study enhances the knowledge on the molecular evolution of emerging non-PCV-13 vaccine serotype 22F, including acquisition of resistance genes through recombination events. It underpins the importance of whole genome sequencing in studying Streptococcus pneumoniae population structures and dynamics, and its utility in molecular surveillance.

Population Structure and Antimicrobial Resistance Profiles of Streptococcus suis Serotype 2 Sequence Type 25 Strains

Strains of serotype 2 Streptococcus suis are responsible for swine and human infections. Different serotype 2 genetic backgrounds have been defined using multilocus sequence typing (MLST). However, little is known about the genetic diversity within each MLST sequence type (ST). Here, we used whole-genome sequencing to test the hypothesis that S. suis serotype 2 strains of the ST25 lineage are genetically heterogeneous. We evaluated 51 serotype 2 ST25 S. suis strains isolated from diseased pigs and humans in Canada, the United States of America, and Thailand. Whole-genome sequencing revealed numerous large-scale rearrangements in the ST25 genome, compared to the genomes of ST1 and ST28 S. suis strains, which result, among other changes, in disruption of a pilus island locus. We report that recombination and lateral gene transfer contribute to ST25 genetic diversity. Phylogenetic analysis identified two main and distinct Thai and North American clades grouping most strains investigated. These clades also possessed distinct patterns of antimicrobial resistance genes, which correlated with acquisition of different integrative and conjugative elements (ICEs). Some of these ICEs were found to be integrated at a recombination hot spot, previously identified as the site of integration of the 89K pathogenicity island in serotype 2 ST7 S. suis strains. Our results highlight the limitations of MLST for phylogenetic analysis of S. suis, and the importance of lateral gene transfer and recombination as drivers of diversity in this swine pathogen and zoonotic agent.

Clonal Complex 17 Group B Streptococcus strains causing invasive disease in neonates and adults originate from the same genetic pool

A significant proportion of group B Streptococcus (GBS) neonatal disease, particularly late-onset disease, is associated with strains of serotype III, clonal complex (CC) 17. CC17 strains also cause invasive infections in adults. Little is known about the phylogenetic relationships of isolates recovered from neonatal and adult CC17 invasive infections. We performed whole-genome-based phylogenetic analysis of 93 temporally and geographically matched CC17 strains isolated from both neonatal and adult invasive infections in the metropolitan region of Toronto/Peel, Canada. We also mined the whole-genome data to reveal mobile genetic elements carrying antimicrobial resistance genes. We discovered that CC17 GBS strains causing neonatal and adult invasive disease are interspersed and cluster tightly in a phylogenetic tree, signifying that they are derived from the same genetic pool. We identified limited variation due to recombination in the core CC17 genome. We describe that loss of Pilus Island 1 and acquisition of different mobile genetic elements carrying determinants of antimicrobial resistance contribute to CC17 genetic diversity. Acquisition of some of these mobile genetic elements appears to correlate with clonal expansion of the strains that possess them. Our results provide a genome-wide portrait of the population structure and evolution of a major disease-causing clone of an opportunistic pathogen.

Population structure and antimicrobial resistance of invasive serotype IV group B Streptococcus, Toronto, Ontario, Canada

Conjugate vaccines should include polysaccharide or virulence proteins of this serotype to provide complete protection.

We recently showed that 37/600 (6.2%) invasive infections with group B Streptococcus (GBS) in Toronto, Ontario, Canada, were caused by serotype IV strains. We report a relatively high level of genetic diversity in 37 invasive strains of this emerging GBS serotype. Multilocus sequence typing identified 6 sequence types (STs) that belonged to 3 clonal complexes. Most isolates were ST-459 (19/37, 51%) and ST-452 (11/37, 30%), but we also identified ST-291, ST-3, ST-196, and a novel ST-682. We detected further diversity by performing whole-genome single-nucleotide polymorphism analysis and found evidence of recombination events contributing to variation in some serotype IV GBS strains. We also evaluated antimicrobial drug resistance and found that ST-459 strains were resistant to clindamycin and erythromycin, whereas strains of other STs were, for the most part, susceptible to these antimicrobial drugs.

High Incidence of Invasive Group A Streptococcus Disease Caused by Strains of Uncommon emm Types in Thunder Bay, Ontario, Canada

An outbreak of type emm59 invasive group A Streptococcus (iGAS) disease was declared in 2008 in Thunder Bay District, Northwestern Ontario, 2 years after a countrywide emm59 epidemic was recognized in Canada. Despite a declining number of emm59 infections since 2010, numerous cases of iGAS disease continue to be reported in the area. We collected clinical information on all iGAS cases recorded in Thunder Bay District from 2008 to 2013. We also emm typed and sequenced the genomes of all available strains isolated from 2011 to 2013 from iGAS infections and from severe cases of soft tissue infections. We used whole-genome sequencing data to investigate the population structure of GAS strains of the most frequently isolated emm types. We report an increased incidence of iGAS in Thunder Bay compared to the metropolitan area of Toronto/Peel and the province of Ontario. Illicit drug use, alcohol abuse, homelessness, and hepatitis C infection were underlying diseases or conditions that might have predisposed patients to iGAS disease. Most cases were caused by clonal strains of skin or generalist emm types (i.e., emm82, emm87, emm101, emm4, emm83, and emm114) uncommonly seen in other areas of the province. We observed rapid waxing and waning of emm types causing disease and their replacement by other emm types associated with the same tissue tropisms. Thus, iGAS disease in Thunder Bay District predominantly affects a select population of disadvantaged persons and is caused by clonally related strains of a few skin and generalist emm types less commonly associated with iGAS in other areas of Ontario.

The majority of 9,729 group A streptococcus strains causing disease secrete SpeB cysteine protease: pathogenesis implications

Group A streptococcus (GAS), the causative agent of pharyngitis and necrotizing fasciitis, secretes the potent cysteine protease SpeB. Several lines of evidence suggest that SpeB is an important virulence factor. SpeB is expressed in human infections, protects mice from lethal challenge when used as a vaccine, and contributes significantly to tissue destruction and dissemination in animal models. However, recent descriptions of mutations in genes implicated in SpeB production have led to the idea that GAS may be under selective pressure to decrease secreted SpeB protease activity during infection. Thus, two divergent hypotheses have been proposed. One postulates that SpeB is a key contributor to pathogenesis; the other, that GAS is under selection to decrease SpeB during infection. In order to distinguish between these alternative hypotheses, we performed casein hydrolysis assays to measure the SpeB protease activity secreted by 6,775 GAS strains recovered from infected humans. The results demonstrated that 84.3% of the strains have a wild-type SpeB protease phenotype. The availability of whole-genome sequence data allowed us to determine the relative frequencies of mutations in genes implicated in SpeB production. The most abundantly mutated genes were direct transcription regulators. We also sequenced the genomes of 2,954 GAS isolates recovered from nonhuman primates with experimental necrotizing fasciitis. No mutations that would result in a SpeB-deficient phenotype were identified. Taken together, these data unambiguously demonstrate that the great majority of GAS strains recovered from infected humans secrete wild-type levels of SpeB protease activity. Our data confirm the important role of SpeB in GAS pathogenesis and help end a long-standing controversy.

Complex Population Structure and Virulence Differences among Serotype 2 Streptococcus suis Strains Belonging to Sequence Type 28

Streptococcus suis is a major swine pathogen and a zoonotic agent. Serotype 2 strains are the most frequently associated with disease. However, not all serotype 2 lineages are considered virulent. Indeed, sequence type (ST) 28 serotype 2 S. suis strains have been described as a homogeneous group of low virulence. However, ST28 strains are often isolated from diseased swine in some countries, and at least four human ST28 cases have been reported. Here, we used whole-genome sequencing and animal infection models to test the hypothesis that the ST28 lineage comprises strains of different genetic backgrounds and different virulence. We used 50 S. suis ST28 strains isolated in Canada, the United States and Japan from diseased pigs, and one ST28 strain from a human case isolated in Thailand. We report a complex population structure among the 51 ST28 strains. Diversity resulted from variable gene content, recombination events and numerous genome-wide polymorphisms not attributable to recombination. Phylogenetic analysis using core genome single-nucleotide polymorphisms revealed four discrete clades with strong geographic structure, and a fifth clade formed by US, Thai and Japanese strains. When tested in experimental animal models, strains from this latter clade were significantly more virulent than a Canadian ST28 reference strain, and a closely related Canadian strain. Our results highlight the limitations of MLST for both phylogenetic analysis and virulence prediction and raise concerns about the possible emergence of ST28 strains in human clinical cases.

Emergence of Serotype IV Group B Streptococcus Adult Invasive Disease in Manitoba and Saskatchewan, Canada, Is Driven by Clonal Sequence Type 459 Strains

Serotype IV group B Streptococcus (GBS) is emerging in Canada and the United States with rates as high as 5% of the total burden of adult invasive GBS disease. To understand this emergence, we studied the population structure and assessed the antimicrobial susceptibility of serotype IV isolates causing adult invasive infection in Manitoba and Saskatchewan, Canada, between 2010 and 2014. Whole-genome sequencing was used to determine multilocus sequence typing information and identify genes encoding antimicrobial resistance in 85 invasive serotype IV GBS strains. Antimicrobial susceptibility testing was performed by standard methods. Strain divergence was assessed using genome-wide single-nucleotide polymorphism analysis. Serotype IV strains were responsible for 16.9% of adult invasive GBS infections in Manitoba and Saskatchewan during the period. The majority of serotype IV isolates (89%) were clonally related, tetracycline-, erythromycin-, and clindamycin-resistant sequence type 459 (ST459) strains that possessed genes tetM and ermTR. Genome comparisons between ST459 and serotype V ST1 GBS identified several areas of recombination in an overall similar genomic background. Serotype IV ST459 GBS strains are expanding and causing a substantial percentage of adult invasive GBS disease. This emergence may be linked to the acquisition of resistance to tetracycline, macrolides, and lincosamides.

Sequence type 1 group B Streptococcus, an emerging cause of invasive disease in adults, evolves by small genetic changes

The molecular mechanisms underlying pathogen emergence in humans is a critical but poorly understood area of microbiologic investigation. Serotype V group B Streptococcus (GBS) was first isolated from humans in 1975, and rates of invasive serotype V GBS disease significantly increased starting in the early 1990s. We found that 210 of 229 serotype V GBS strains (92%) isolated from the bloodstream of nonpregnant adults in the United States and Canada between 1992 and 2013 were multilocus sequence type (ST) 1. Elucidation of the complete genome of a 1992 ST-1 strain revealed that this strain had the highest homology with a GBS strain causing cow mastitis and that the 1992 ST-1 strain differed from serotype V strains isolated in the late 1970s by acquisition of cell surface proteins and antimicrobial resistance determinants. Whole-genome comparison of 202 invasive ST-1 strains detected significant recombination in only eight strains. The remaining 194 strains differed by an average of 97 SNPs. Phylogenetic analysis revealed a temporally dependent mode of genetic diversification consistent with the emergence in the 1990s of ST-1 GBS as major agents of human disease. Thirty-one loci were identified as being under positive selective pressure, and mutations at loci encoding polysaccharide capsule production proteins, regulators of pilus expression, and two-component gene regulatory systems were shown to affect the bacterial phenotype. These data reveal that phenotypic diversity among ST-1 GBS is mainly driven by small genetic changes rather than extensive recombination, thereby extending knowledge into how pathogens adapt to humans.

PBP2a mutations causing high-level Ceftaroline resistance in clinical methicillin-resistant Staphylococcus aureus isolates

Ceftaroline is the first member of a novel class of cephalosporins approved for use in the United States. Although prior studies have identified eight ceftaroline-resistant methicillin-resistant Staphylococcus aureus (MRSA) isolates in Europe and Asia with MICs ranging from 4 to 8 mg/liter, high-level resistance to ceftaroline (>32 mg/liter) has not been described in MRSA strains isolated in the United States. We isolated a ceftaroline-resistant (MIC > 32 mg/liter) MRSA strain from the blood of a cystic fibrosis patient and five MRSA strains from the respiratory tract of this patient. Whole-genome sequencing identified two amino acid-altering mutations uniquely present in the ceftaroline-binding pocket of the transpeptidase region of penicillin-binding protein 2a (PBP2a) in ceftaroline-resistant isolates. Biochemical analyses and the study of isogenic mutant strains confirmed that these changes caused ceftaroline resistance. Thus, we identified the molecular mechanism of ceftaroline resistance in the first MRSA strain with high-level ceftaroline resistance isolated in the United States.

Absence of patient-to-patient intrahospital transmission of Staphylococcus aureus as determined by whole-genome sequencing

Nosocomial transmission of pathogens is a major health care challenge. The increasing spread of antibiotic-resistant strains represents an ongoing threat to public health. Previous Staphylococcus aureus transmission studies have focused on transmission of S. aureus between asymptomatic carriers or used low-resolution typing methods such as multilocus sequence typing (MLST) or spa typing. To identify patient-to-patient intrahospital transmission using high-resolution genetic analysis, we sequenced the genomes of a consecutive set of 398 S. aureus isolates from sterile-site infections. The S. aureus strains were collected from four hospitals in the Houston Methodist Hospital System over a 6-month period. Importantly, we discovered no evidence of transmission of S. aureus between patients with sterile-site infections. The lack of intrahospital transmission may reflect a fundamental difference between day-to-day transmission events in the hospital setting and the more frequently studied outbreak scenarios.

Previous studies have suggested that nosocomial transmission of S. aureus is common. Our data revealed an unexpected lack of evidence for intrahospital transmission of S. aureus between patients with invasive infections. This finding has important implications for hospital infection control and public health efforts. In addition, our data demonstrate that highly related pools of S. aureus strains exist in the community which may complicate outbreak investigations.

Isolation of Salmonella enterica serovar Kentucky strain ST 198 and its H2S-negative variant from a patient: implications for diagnosis

H2S-producing multiresistant Salmonella enterica serovar Kentucky strain sequence type (ST) 198 and its non-H2S-producing variant were isolated from a patient. Whole-genome comparison showed a base addition in the gene encoding molybdenum cofactor biosynthesis protein C, which could affect H2S production in the variant. Lack of H2S production has implications for diagnosis of salmonella.

Multihost experimental evolution of the pathogen Ralstonia solanacearum unveils genes involved in adaptation to plants

Ralstonia solanacearum, the causal agent of a lethal bacterial wilt plant disease, infects an unusually wide range of hosts. These hosts can further be split into plants where R. solanacearum is known to cause disease (original hosts) and those where this bacterium can grow asymptomatically (distant hosts). Moreover, this pathogen is able to adapt to many plants as supported by field observations reporting emergence of strains with enlarged pathogenic properties. To investigate the genetic bases of host adaptation, we conducted evolution experiments by serial passages of a single clone of the pathogen on three original and two distant hosts over 300 bacterial generations and then analyzed the whole-genome of nine evolved clones. Phenotypic analysis of the evolved clones showed that the pathogen can increase its fitness on both original and distant hosts although the magnitude of fitness increase was greater on distant hosts. Only few genomic modifications were detected in evolved clones compared with the ancestor but parallel evolutionary changes in two genes were observed in independent evolved populations. Independent mutations in the regulatory gene efpR were selected for in three populations evolved on beans, a distant host. Reverse genetic approaches confirmed that these mutations were associated with fitness gain on bean plants. This work provides a first step toward understanding the within-host evolutionary dynamics of R. solanacearum during infection and identifying bacterial genes subjected to in planta selection. The discovery of EfpR as a determinant conditioning host adaptation of the pathogen illustrates how experimental evolution coupled with whole-genome sequencing is a potent tool to identify novel molecular players involved in central life-history traits.

Whole-genome sequencing of the Mycobacterium tuberculosis Manila sublineage results in less clustering and better resolution than mycobacterial interspersed repetitive-unit-variable-number tandem-repeat (MIRU-VNTR) typing and spoligotyping

Mycobacterium tuberculosis isolates of the Manila sublineage are genetically homogeneous. In this study, we used whole-genome sequencing (WGS) to type a collection of 36 M. tuberculosis isolates of the Manila family. WGS enabled the subtyping of these 36 isolates into at least 10 distinct clusters. Our results indicate that WGS is a powerful approach to determining the relatedness of Manila family M. tuberculosis isolates.

Sensitivity to the two-peptide bacteriocin lactococcin G is dependent on UppP, an enzyme involved in cell-wall synthesis

Most bacterially produced antimicrobial peptides (bacteriocins) are thought to kill target cells by a receptor-mediated mechanism. However, for most bacteriocins the receptor is unknown. For instance, no target receptor has been identified for the two-peptide bacteriocins (class IIb), whose activity requires the combined action of two individual peptides. To identify the receptor for the class IIb bacteriocin lactococcin G, which targets strains of Lactococcus lactis, we generated 12 lactococcin G-resistant mutants and performed whole-genome sequencing to identify mutations causing the resistant phenotype. Remarkably, all had a mutation in or near the gene uppP (bacA), encoding an undecaprenyl pyrophosphate phosphatase; a membrane protein involved in peptidoglycan synthesis. Nine mutants had stop codons or frameshifts in the uppP gene, two had point mutations in putative regulatory regions and one caused an amino acid substitution in UppP. To verify the receptor function of UppP, it was shown that growth of non-sensitive Streptococcus pneumoniae could be inhibited by lactococcin G when L. lactis uppP was expressed in this bacterium. Furthermore, we show that the related class IIb bacteriocin enterocin 1071 also uses UppP as receptor. The approach used here should be broadly applicable to identify receptors for other bacteriocins as well.

Characterization of invasive group B streptococcus strains from the greater Toronto area, Canada

We determined the capsular polysaccharide (CPS) type of 600 group B Streptococcus (GBS) (also known as Streptococcus agalactiae) strains recovered from patients with invasive infections in the greater Toronto area, Canada, between 2009 and 2012. GBS strains of CPS type III were the most prevalent among infants (44% in those with early-onset disease, 75% in those with late-onset disease), while type V strains were most frequently isolated from adult patients (26% in patients≥19 years old). We next investigated the presence in our collection of GBS strains belonging to the hypervirulent multilocus sequence typing clonal complex 17 (CC17). We used a PCR test described as specific for the detection of CC17 strains, which targets the gene encoding the major virulence factor HvgA. We identified 91 hvgA-positive strains; of these, 88 were CPS type III, 2 were CPS type IV, and 1 was CPS type V. Using whole-genome sequencing, we showed that the two hvgA-positive CPS type IV strains are CC17 strains which underwent capsular switching. However, sequence analysis revealed that the hvgA-positive CPS type V strain does not belong to CC17 but instead is a bona fide CC1 strain which acquired hvgA, probably by recombination from a CC17 donor. Our findings underline the importance of recombination in GBS pathogenesis and caution against the use of single-gene-based PCR tests to detect CC17 GBS strains.

A Zn-dependent metallopeptidase is responsible for sensitivity to LsbB, a class II leaderless bacteriocin of Lactococcus lactis subsp. lactis BGMN1-5

Lactococcus lactis subsp. lactis BGMN1-5 produces a leaderless class II bacteriocin called LsbB. To identify the receptor for LsbB, a cosmid library of the LsbB-sensitive strain BGMN1-596 was constructed. About 150 cosmid clones were individually isolated and transferred to LsbB-resistant mutants of BGMN1-596. Cosmid pAZILcos/MN2, carrying a 40-kb insert, was found to restore LsbB sensitivity in LsbB-resistant mutants. Further subcloning revealed that a 1.9-kb fragment, containing only one open reading frame, was sufficient to restore sensitivity. The fragment contains the gene yvjB coding for a Zn-dependent membrane-bound metallopeptidase, suggesting that this gene may serve as the receptor for LsbB. Further support for this notion derives from several independent experiments: (i) whole-genome sequencing confirmed that all LsbB-resistant mutants contain mutations in yvjB; (ii) disruption of yvjB by direct gene knockout rendered sensitive strains BGMN1-596 and IL1403 resistant to LsbB; and (iii) most compellingly, heterologous expression of yvjB in naturally resistant strains of other species, such as Lactobacillus paracasei and Enterococcus faecalis, also rendered them sensitive to the bacteriocin. To our knowledge, this is the first time a membrane-bound peptidase gene has been shown to be involved in bacteriocin sensitivity in target cells. We also demonstrated a novel successful approach for identifying bacteriocin receptors.

Bi-modal distribution of the second messenger c-di-GMP controls cell fate and asymmetry during the caulobacter cell cycle

Many bacteria mediate important life-style decisions by varying levels of the second messenger c-di-GMP. Behavioral transitions result from the coordination of complex cellular processes such as motility, surface adherence or the production of virulence factors and toxins. While the regulatory mechanisms responsible for these processes have been elucidated in some cases, the global pleiotropic effects of c-di-GMP are poorly understood, primarily because c-di-GMP networks are inherently complex in most bacteria. Moreover, the quantitative relationships between cellular c-di-GMP levels and c-di-GMP dependent phenotypes are largely unknown. Here, we dissect the c-di-GMP network of Caulobacter crescentus to establish a global and quantitative view of c-di-GMP dependent processes in this organism. A genetic approach that gradually reduced the number of diguanylate cyclases identified novel c-di-GMP dependent cellular processes and unraveled c-di-GMP as an essential component of C. crescentus cell polarity and its bimodal life cycle. By varying cellular c-di-GMP concentrations, we determined dose response curves for individual c-di-GMP-dependent processes. Relating these values to c-di-GMP levels modeled for single cells progressing through the cell cycle sets a quantitative frame for the successive activation of c-di-GMP dependent processes during the C. crescentus life cycle. By reconstructing a simplified c-di-GMP network in a strain devoid of c-di-GMP we defined the minimal requirements for the oscillation of c-di-GMP levels during the C. crescentus cell cycle. Finally, we show that although all c-di-GMP dependent cellular processes were qualitatively restored by artificially adjusting c-di-GMP levels with a heterologous diguanylate cyclase, much higher levels of the second messenger are required under these conditions as compared to the contribution of homologous c-di-GMP metabolizing enzymes. These experiments suggest that a common c-di-GMP pool cannot fully explain spatiotemporal regulation by c-di-GMP in C. crescentus and that individual enzymes preferentially regulate specific phenotypes during the cell cycle.

Bacterial processes like virulence, motility or biofilm formation are governed by the second messenger c-di-GMP. In most bacteria, c-di-GMP is produced and degraded by a complex network comprising dozens of enzymes. This has hindered a comprehensive analysis of the cellular role of c-di-GMP. Here we mutate the entire c-di-GMP network in Caulobacter crescentus, a model organism with inherent cell polarity and bimodal life-cycle. We find that a c-di-GMP free strain (cdG⁰) shows severe developmental defects, a loss of cell polarity and defective cell division. By determining c-di-GMP dose-response curves for individual processes and relating these to c-di-GMP levels, which were modeled for single cells progressing through the cell cycle, we define a quantitative frame for the c-di-GMP dependent program during the C. crescentus life cycle. We then show that the defects of cdG⁰ can be largely rescued by restoring c-di-GMP levels with a single heterologous enzyme producing c-di-GMP. However, much higher levels of the second messenger are required under these conditions as compared to the contribution of homologous enzymes. Our data argue for specific regulatory fine-tuning of the enzymes mediating c-di-GMP oscillation during the cell cycle and provide evidence for both global as well as insulated c-di-GMP pools.

Bacterial genome evolution within a clonal population: from in vitro investigations to in vivo observations

Bacteria are faced with a diversity of environmental stresses that include high salt concentrations, heavy metals and pH fluctuations. Adaptation to resist such stresses is a complex phenomenon that involves global pathways and simultaneous acquisition of multiple unrelated properties. During the last 3 years, the development of new technologies in the field of molecular biology has led to numerous fundamental and quantitative in vitro and in vivo evolutionary studies that have improved our understanding of the principles underlying bacterial adaptations, and helped us develop strategies to cope with the health burden of bacterial virulence. In this review, the authors discuss the evolution of bacteria in the laboratory and in human patients.

Integrated whole-genome sequencing and temporospatial analysis of a continuing Group A Streptococcus epidemic

Analysis of microbial epidemics has been revolutionized by whole-genome sequencing. We recently sequenced the genomes of 601 type emm59 Group A Streptococcus (GAS) organisms responsible for an ongoing epidemic of invasive infections in Canada and some of the United States. The epidemic has been caused by the emergence of a genetically distinct, hypervirulent clone that has genetically diversified. The ease of obtaining genomic data contrasts with the relatively difficult task of translating them into insightful epidemiological information. Here, we sequenced the genomes of 90 additional invasive Canadian emm59 GAS organisms, including 80 isolated recently in 2010–2011. We used an improved bioinformatics pipeline designed to rapidly process and analyze whole-genome data and integrate strain metadata. We discovered that emm59 GAS organisms are undergoing continued multiclonal evolutionary expansion. Previously identified geographic patterns of strain dissemination are being diluted as mixing of subclones over time and space occurs. Our integrated data analysis strategy permits prompt and accurate mapping of the dissemination of bacterial organisms in an epidemic wave, permitting rapid generation of hypotheses that inform public health and virulence studies.

Aminoglycoside cross-resistance in Mycobacterium tuberculosis due to mutations in the 5' untranslated region of whiB7

Since the discovery of streptomycin's bactericidal activity against Mycobacterium tuberculosis, aminoglycosides have been utilized to treat tuberculosis (TB). Today, the aminoglycosides kanamycin and amikacin are used to treat multidrug-resistant (MDR) TB, and resistance to any of the second-line injectable antibiotics, including kanamycin, amikacin, or capreomycin, is a defining characteristic of extensively drug-resistant (XDR) TB. Resistance to kanamycin and streptomycin is thought to be due to the acquisition of unlinked chromosomal mutations. However, we identified eight independent mutations in the 5' untranslated region of the transcriptional activator whiB7 that confer low-level resistance to both aminoglycosides. The mutations lead to 23- to 145-fold increases in whiB7 transcripts and subsequent increased expression of both eis (Rv2416c) and tap (Rv1258c). Increased expression of eis confers kanamycin resistance in these mutants, while increased expression of tap, which encodes an efflux pump, is a previously uncharacterized mechanism of low-level streptomycin resistance. Additionally, high-level resistance to streptomycin arose at a much higher frequency in whiB7 mutants than in a wild-type (WT) strain. Although whiB7 is typically associated with intrinsic antibiotic resistance in M. tuberculosis, these data suggest that mutations in an uncharacterized regulatory region of whiB7 contribute to cross-resistance against clinically used second-line antibiotics. As drug resistance continues to develop and spread, understanding the mechanisms and molecular basis of antibiotic resistance is critical for the development of rapid molecular tests to diagnose drug-resistant TB strains and ultimately for designing regimens to treat drug-resistant cases of TB.

The role of clonal interference in the evolutionary dynamics of plasmid-host adaptation

Promiscuous plasmids replicate in a wide range of bacteria and therefore play a key role in the dissemination of various host-beneficial traits, including antibiotic resistance. Despite the medical relevance, little is known about the evolutionary dynamics through which drug resistance plasmids adapt to new hosts and thereby persist in the absence of antibiotics. We previously showed that the incompatibility group P-1 (IncP-1) minireplicon pMS0506 drastically improved its stability in novel host Shewanella oneidensis MR-1 after 1,000 generations under antibiotic selection for the plasmid. The only mutations found were those affecting the N terminus of the plasmid replication initiation protein TrfA1. Our aim in this study was to gain insight into the dynamics of plasmid evolution. Changes in stability and genotype frequencies of pMS0506 were monitored in evolving populations of MR-1 (pMS0506). Genotypes were determined by sequencing trfA1 amplicons from individual clones and by 454 pyrosequencing of whole plasmids from entire populations. Stability of pMS0506 drastically improved by generation 200. Many evolved plasmid genotypes with point mutations as well as in-frame and frameshift deletions and duplications in trfA1 were observed in all lineages with both sequencing methods. Strikingly, multiple genotypes were simultaneously present at high frequencies (>10%) in each population. Their relative abundances changed over time, but after 1,000 generations only one or two genotypes dominated the populations. This suggests that hosts with different plasmid genotypes were competing with each other, thus affecting the evolutionary trajectory. Plasmids can thus rapidly improve their stability, and clonal interference plays a significant role in plasmid-host adaptation dynamics.

Promiscuous plasmids play an important role in the spread of antibiotic resistance and many other traits between closely and distantly related bacteria. However, little is known about the dynamics by which these broad-host-range antibiotic resistance plasmids adapt to novel bacteria and thereby become more persistent, even in the absence of antibiotics. In this study, we show that after no more than 200 generations of growth in the presence of antibiotics, a plasmid that was initially poorly maintained in a novel bacterial host evolved to become drastically more persistent in the absence of antibiotics. In each of the evolving populations, an unexpectedly large number of bacterial variants arose with distinct mutations in the plasmid’s replication initiation protein. Our results suggest that clonal interference, characterized by competition between variant clones in a population, plays a major role in the evolution of the persistence of drug resistance.

Serine/threonine phosphatase Stp1 contributes to reduced susceptibility to vancomycin and virulence in Staphylococcus aureus

The genetic mechanisms that contribute to reduced susceptibility to vancomycin in Staphylococcus aureus are complex and heterogeneous. In addition, debate is emerging as to the true effect of reduced susceptibility to vancomycin on staphylococcal virulence. To investigate this, comparative genomics was performed on a collection of vancomycin-exposed isogenic S. aureus pairs (14 strains in total). Previously described mutations were observed in genes such as vraG, agrA, yvqF, and rpoB; however, a new mechanism was identified involving a serine/threonine phosphatase, Stp1. After constructing an stp1 deletion mutant, we showed that stp1 is important in vancomycin susceptibility and cell wall biosynthesis. Gene expression studies showed that stp1 also regulates virulence genes, including a hemolysin, superantigen-like protein, and phenol-soluble modulin, and that the deletion mutant is attenuated in virulence in vivo. Stp1 provides a new link between vancomycin susceptibility and virulence in S. aureus.

Full-genome dissection of an epidemic of severe invasive disease caused by a hypervirulent, recently emerged clone of group A Streptococcus

Group A Streptococcus (GAS) causes an exceptionally broad range of infections in humans, from relatively mild pharyngitis and skin infections to life-threatening necrotizing fasciitis and toxic shock syndrome. An epidemic of severe invasive human infections caused by type emm59 GAS, heretofore an exceedingly rare cause of disease, spread west to east across Canada over a 3-year period (2006 to 2008). By sequencing the genomes of 601 epidemic, historic, and other emm59 organisms, we discovered that a recently emerged, genetically distinct emm59 clone is responsible for the Canadian epidemic. Using near-real-time genome sequencing, we were able to show spread of the Canadian epidemic clone into the United States. The extensive genome data permitted us to identify patterns of geographic dissemination as well as links between emm59 subclonal lineages that cause infections. Mouse and nonhuman primate models of infection demonstrated that the emerged clone is unusually virulent. Transmission of epidemic emm59 strains may have occurred primarily by skin contact, as suggested by an experimental model of skin transmission. In addition, the emm59 strains had a significantly impaired ability to persist in human saliva and to colonize the oropharynx of mice, and seldom caused human pharyngitis. Our study contributes new information to the rapidly emerging field of molecular pathogenomics of bacterial epidemics and illustrates how full-genome data can be used to precisely illuminate the landscape of strain dissemination during a bacterial epidemic.

Whole genome characterization of the mechanisms of daptomycin resistance in clinical and laboratory derived isolates of Staphylococcus aureus

Background

Daptomycin remains one of our last-line anti-staphylococcal agents. This study aims to characterize the genetic evolution to daptomycin resistance in S. aureus.

Methods

Whole genome sequencing was performed on a unique collection of isogenic, clinical (21 strains) and laboratory (12 strains) derived strains that had been exposed to daptomycin and developed daptomycin-nonsusceptibility. Electron microscopy (EM) and lipid membrane studies were performed on selected isolates.

Results

On average, six coding region mutations were observed across the genome in the clinical daptomycin exposed strains, whereas only two mutations on average were seen in the laboratory exposed pairs. All daptomycin-nonsusceptible strains had a mutation in a phospholipid biosynthesis gene. This included mutations in the previously described mprF gene, but also in other phospholipid biosynthesis genes, including cardiolipin synthase (cls2) and CDP-diacylglycerol-glycerol-3-phosphate 3-phosphatidyltransferase (pgsA). EM and lipid membrane composition analyses on two clinical pairs showed that the daptomycin-nonsusceptible strains had a thicker cell wall and an increase in membrane lysyl-phosphatidylglycerol.

Conclusion

Point mutations in genes coding for membrane phospholipids are associated with the development of reduced susceptibility to daptomycin in S. aureus. Mutations in cls2 and pgsA appear to be new genetic mechanisms affecting daptomycin susceptibility in S. aureus.

LOCAS--a low coverage assembly tool for resequencing projects

MOTIVATION

Next Generation Sequencing (NGS) is a frequently applied approach to detect sequence variations between highly related genomes. Recent large-scale re-sequencing studies as the Human 1000 Genomes Project utilize NGS data of low coverage to afford sequencing of hundreds of individuals. Here, SNPs and micro-indels can be detected by applying an alignment-consensus approach. However, computational methods capable of discovering other variations such as novel insertions or highly diverged sequence from low coverage NGS data are still lacking.

RESULTS

We present LOCAS, a new NGS assembler particularly designed for low coverage assembly of eukaryotic genomes using a mismatch sensitive overlap-layout-consensus approach. LOCAS assembles homologous regions in a homology-guided manner while it performs de novo assemblies of insertions and highly polymorphic target regions subsequently to an alignment-consensus approach. LOCAS has been evaluated in homology-guided assembly scenarios with low sequence coverage of Arabidopsis thaliana strains sequenced as part of the Arabidopsis 1001 Genomes Project. While assembling the same amount of long insertions as state-of-the-art NGS assemblers, LOCAS showed best results regarding contig size, error rate and runtime.

CONCLUSION

LOCAS produces excellent results for homology-guided assembly of eukaryotic genomes with short reads and low sequencing depth, and therefore appears to be the assembly tool of choice for the detection of novel sequence variations in this scenario.

Site-specific mobilization of vinyl chloride respiration islands by a mechanism common in Dehalococcoides

Background

Vinyl chloride is a widespread groundwater pollutant and Group 1 carcinogen. A previous comparative genomic analysis revealed that the vinyl chloride reductase operon, vcrABC, of Dehalococcoides sp. strain VS is embedded in a horizontally-acquired genomic island that integrated at the single-copy tmRNA gene, ssrA.

Results

We targeted conserved positions in available genomic islands to amplify and sequence four additional vcrABC -containing genomic islands from previously-unsequenced vinyl chloride respiring Dehalococcoides enrichments. We identified a total of 31 ssrA-specific genomic islands from Dehalococcoides genomic data, accounting for 47 reductive dehalogenase homologous genes and many other non-core genes. Sixteen of these genomic islands contain a syntenic module of integration-associated genes located adjacent to the predicted site of integration, and among these islands, eight contain vcrABC as genetic 'cargo'. These eight vcrABC -containing genomic islands are syntenic across their ~12 kbp length, but have two phylogenetically discordant segments that unambiguously differentiate the integration module from the vcrABC cargo. Using available Dehalococcoides phylogenomic data we estimate that these ssrA-specific genomic islands are at least as old as the Dehalococcoides group itself, which in turn is much older than human civilization.

Conclusions

The vcrABC -containing genomic islands are a recently-acquired subset of a diverse collection of ssrA-specific mobile elements that are a major contributor to strain-level diversity in Dehalococcoides, and may have been throughout its evolution. The high similarity between vcrABC sequences is quantitatively consistent with recent horizontal acquisition driven by ~100 years of industrial pollution with chlorinated ethenes.

Inviability of a DNA2 deletion mutant is due to the DNA damage checkpoint

Dna2 is a dual polarity exo/endonuclease, and 5' to 3' DNA helicase involved in Okazaki Fragment Processing (OFP) and Double-Strand Break (DSB) Repair. In yeast, DNA2 is an essential gene, as expected for a DNA replication protein. Suppression of the lethality of dna2Δ mutants has been found to occur by two mechanisms: overexpression of RAD27 (scFEN1) , encoding a 5' to 3' exo/endo nuclease that processes Okazaki fragments (OFs) for ligation, or deletion of PIF1, a 5' to 3' helicase involved in mitochondrial recombination, telomerase inhibition and OFP. Mapping of a novel, spontaneously arising suppressor of dna2Δ now reveals that mutation of rad9 and double mutation of rad9 mrc1 can also suppress the lethality of dna2Δ mutants. Interaction of dna2Δ and DNA damage checkpoint mutations provides insight as to why dna2Δ is lethal but rad27Δ is not, even though evidence shows that Rad27 (ScFEN1) processes most of the Okazaki fragments, while Dna2 processes only a subset.

Probing bacterial pathogenesis with genetics, genomics, and chemical biology: past, present, and future approaches

Classical genetic approaches for studying bacterial pathogenesis have provided a solid foundation for our current understanding of microbial physiology and the interactions between pathogen and host. During the past decade however, advances in several arenas have expanded the ways in which the biology of pathogens can be studied. This review discussed the impact of these advances on bacterial genetics, including the application of genomics and chemical biology to the study of pathogenesis.

Distinct signatures of diversifying selection revealed by genome analysis of respiratory tract and invasive bacterial populations

Many pathogens colonize different anatomical sites, but the selective pressures contributing to survival in the diverse niches are poorly understood. Group A Streptococcus (GAS) is a human-adapted bacterium that causes a range of infections. Much effort has been expended to dissect the molecular basis of invasive (sterile-site) infections, but little is known about the genomes of strains causing pharyngitis (streptococcal "sore throat"). Additionally, there is essentially nothing known about the genetic relationships between populations of invasive and pharyngitis strains. In particular, it is unclear if invasive strains represent a distinct genetic subpopulation of strains that cause pharyngitis. We compared the genomes of 86 serotype M3 GAS pharyngitis strains with those of 215 invasive M3 strains from the same geographical location. The pharyngitis and invasive groups were highly related to each other and had virtually identical phylogenetic structures, indicating they belong to the same genetic pool. Despite the overall high degree of genetic similarity, we discovered that strains from different host environments (i.e., throat, normally sterile sites) have distinct patterns of diversifying selection at the nucleotide level. In particular, the pattern of polymorphisms in the hyaluronic acid capsule synthesis operon was especially different between the two strain populations. This finding was mirrored by data obtained from full-genome analysis of strains sequentially cultured from nonhuman primates. Our results answer the long-standing question of the genetic relationship between GAS pharyngitis and invasive strains. The data provide previously undescribed information about the evolutionary history of pathogenic microbes that cause disease in different anatomical sites.

Analyzing and minimizing PCR amplification bias in Illumina sequencing libraries

Despite the ever-increasing output of Illumina sequencing data, loci with extreme base compositions are often under-represented or absent. To evaluate sources of base-composition bias, we traced genomic sequences ranging from 6% to 90% GC through the process by quantitative PCR. We identified PCR during library preparation as a principal source of bias and optimized the conditions. Our improved protocol significantly reduces amplification bias and minimizes the previously severe effects of PCR instrument and temperature ramp rate.

Analyzing and minimizing PCR amplification bias in Illumina sequencing libraries

Despite the ever-increasing output of Illumina sequencing data, loci with extreme base compositions are often under-represented or absent. To evaluate sources of base-composition bias, we traced genomic sequences ranging from 6% to 90% GC through the process by quantitative PCR. We identified PCR during library preparation as a principal source of bias and optimized the conditions. Our improved protocol significantly reduces amplification bias and minimizes the previously severe effects of PCR instrument and temperature ramp rate.