Deriving group A Streptococcus typing information from short-read whole-genome sequencing data

The clinical and epidemiological impacts of whole genomic sequencing on bacterial and virological agents

Whole Genome Sequencing (WGS) is a molecular biology tool consisting in the sequencing of the entire genome of a given organism. Due to its ability to provide the finest available resolution of bacterial and virological genetics, it is used at several levels in the field of infectiology. On an individual scale and through application of a single technique, it enables the typological identification and characterization of strains, the characterization of plasmids, and enhanced search for resistance genes and virulence factors. On a collective scale, it enables the characterization of strains and the determination of phylogenetic links between different microorganisms during community outbreaks and healthcare-associated epidemics. The information provided by WGS enables real-time monitoring of strain-level epidemiology on a worldwide scale, and facilitates surveillance of the resistance dissemination and the introduction or emergence of pathogenic variants in humans or their environment. There are several possible approaches to completion of an entire genome. The choice of one method rather than another is essentially dictated by the matrix, either a clinical sample or a culture isolate, and the clinical objective. WGS is an advanced technology that remains costly despite a gradual decrease in its expenses, potentially hindering its implementation in certain laboratories and thus its use in routine microbiology. Even though WGS is making steady inroads as a reference method, efforts remain needed in view of so harmonizing its interpretations and decreasing the time to generation of conclusive results.

Promiscuous evolution of Group A Streptococcal M and M-like proteins

Group A Streptococcus (GAS) M and M-like proteins are essential virulence factors and represent the primary epidemiological marker of this pathogen. Protein sequences encoding 1054 M, Mrp and Enn proteins, from 1668 GAS genomes, were analysed by SplitsTree4, partitioning around medoids and co-occurrence. The splits network and groups-based analysis of all M and M-like proteins revealed four large protein groupings, with multiple evolutionary histories as represented by multiple edges for most splits, leading to 'M-family-groups' (FG) of protein sequences: FG I, Mrp; FG II, M protein and Protein H; FG III, Enn; and FG IV, M protein. M and Enn proteins formed two groups with nine sub-groups and Mrp proteins formed four groups with ten sub-groups. Discrete co-occurrence of M and M-like proteins were identified suggesting that while dynamic, evolution may be constrained by a combination of functional and virulence attributes. At a granular level, four distinct family-groups of M, Enn and Mrp proteins are observable, with Mrp representing the most genetically distinct of the family-group of proteins. While M and Enn protein families generally group into three distinct family-groups, horizontal and vertical gene flow between distinct GAS strains is ongoing.

Distribution of emm, superantigen and other virulence genotypes and detection of phylogenetic relationships in group A streptococcal isolates

Group A streptococci are important pathogens with various virulence factors, such as M protein, superantigens, hemolysins, deoxyribonuclease, and proteases. The aims of this study are to investigate the detection of emm genotypes and other virulence genes, such as SAgs, DNase, protease, antibiotic resistance, and phylogenetic relationships in GAS strains isolated from clinical samples.Test strains were obtained from Çukurova University Balcalı Hospital and regional hospitals in Adana province. The M proteins were detected by sequence analysis of emm genes. SAgs and other virulence gene profiles were determined using the Multiplex-PCR method. The antibiotic susceptibility of the isolates was performed by the disc diffusion method and evaluated according to CLSI criteria. The PFGE method was used to determine the clonal relationship between the strains.The emm gene was positive in 86 isolates. The most common emm genotypes were emm28 (22%), emm1 (18.6%), emm12 (13.9%), and emm3 (11.6%). Also, the most common virulence genes were speG (58.1%), speC (56.9%), sdaB (53.4%), and mac (53.4%). The rates of resistance to erythromycin, clindamycin, levofloxacin, ciprofloxacin and telithromycin were 19.8%, 16.3%, 4.7%, 3.5%, and 3.5%, respectively.As a result, additional regional studies on the detection and prevalence of GAS virulence factors in Turkey are required. We believe that this study will provide valuable information for epidemiological studies on emm sequences, Sags, and other virulence factors of Streptococcus pyogenes in Turkey.

Deciphering mobile genetic elements disseminating macrolide resistance in Streptococcus pyogenes over a 21 year period in Barcelona, Spain

OBJECTIVES

To phenotypically and genetically characterize the antibiotic resistance determinants and associated mobile genetic elements (MGEs) among macrolide-resistant (MR) Streptococcus pyogenes [Group A streptococci (GAS)] clinical isolates collected in Barcelona, Spain.

METHODS

Antibiotic susceptibility testing was performed by microdilution. Isolates were emm and MLST typed and 55 were whole-genome sequenced to determine the nature of the macrolide resistance (MR) determinants and their larger MGE and chromosomal context.

RESULTS

Between 1998 and 2018, 142 of 1028 GAS (13.8%) were MR. Among 108 isolates available for molecular characterization, 41.7% had cMLSB, 30.5% iMLSB and 27.8% M phenotype. Eight erm(B)-containing strains were notable in having an MDR phenotype conferred by an MGE encoding several antibiotic resistance genes. MR isolates were comprised of several distinct genetic lineages as defined by the combination of emm and ST. Although most lineages were only transiently present, the emm11/ST403 clone persisted throughout the period. Two lineages, emm9/ST75 with erm(B) and emm77/ST63 with erm(TR), emerged in 2016-18. The erm(B) was predominantly encoded on the Tn916 family of transposons (21/31) with different genetic contexts, and in other MGEs (Tn6263, ICESpHKU372 and one harbouring an MDR cluster called ICESp1070HUB). The erm(TR) was found in ICESp2905 (8/17), ICESp1108-like (4/17), ICESpHKU165 (3/17) and two structures described in this study (IMESp316HUB and ICESp3729HUB). The M phenotype [mef(A)-msr(D)] was linked to phage φ1207.3. Eight integrative conjugative element/integrative mobilizable element (ICE/IME) cluster groups were classified on the basis of gene content within conjugation modules. These groups were found among MGEs, which corresponded with the MR-containing element or the site of integration.

CONCLUSIONS

We detected several different MGEs harbouring erm(B) or erm(TR). This is the first known description of Tn6263 in GAS and three MGEs [IMESp316HUB, ICESp3729HUB and ICESp1070HUB] associated with MR. Periods of high MR rates in our area were mainly associated with the expansion of certain predominant lineages, while in low MR periods different sporadic and low prevalence lineages were more frequent.

Application of the random forest algorithm to Streptococcus pyogenes response regulator allele variation: from machine learning to evolutionary models

Group A Streptococcus (GAS) is a globally significant bacterial pathogen. The GAS genotyping gold standard characterises the nucleotide variation of emm, which encodes a surface-exposed protein that is recombinogenic and under immune-based selection pressure. Within a supervised learning methodology, we tested three random forest (RF) algorithms (Guided, Ordinary, and Regularized) and 53 GAS response regulator (RR) allele types to infer six genomic traits (emm-type, emm-subtype, tissue and country of sample, clinical outcomes, and isolate invasiveness). The Guided, Ordinary, and Regularized RF classifiers inferred the emm-type with accuracies of 96.7%, 95.7%, and 95.2%, using ten, three, and four RR alleles in the feature set, respectively. Notably, we inferred the emm-type with 93.7% accuracy using only mga2 and lrp. We demonstrated a utility for inferring emm-subtype (89.9%), country (88.6%), invasiveness (84.7%), but not clinical (56.9%), or tissue (56.4%), which is consistent with the complexity of GAS pathophysiology. We identified a novel cell wall-spanning domain (SF5), and proposed evolutionary pathways depicting the 'contrariwise' and 'likewise' chimeric deletion-fusion of emm and enn. We identified an intermediate strain, which provides evidence of the time-dependent excision of mga regulon genes. Overall, our workflow advances the understanding of the GAS mga regulon and its plasticity.

Cryptic prophages within a Streptococcus pyogenes genotype emm4 lineage

The major human pathogen Streptococcus pyogenes shares an intimate evolutionary history with mobile genetic elements, which in many cases carry genes encoding bacterial virulence factors. During recent whole-genome sequencing of a longitudinal sample of S. pyogenes isolates in England, we identified a lineage within emm4 that clustered with the reference genome MEW427. Like MEW427, this lineage was characterized by substantial gene loss within all three prophage regions, compared to MGAS10750 and isolates outside of the MEW427-like lineage. Gene loss primarily affected lysogeny, replicative and regulatory modules, and to a lesser and more variable extent, structural genes. Importantly, prophage-encoded superantigen and DNase genes were retained in all isolates. In isolates where the prophage elements were complete, like MGAS10750, they could be induced experimentally, but not in MEW427-like isolates with degraded prophages. We also found gene loss within the chromosomal island SpyCIM4 of MEW427-like isolates, although surprisingly, the SpyCIM4 element could not be experimentally induced in either MGAS10750-like or MEW427-like isolates. This did not, however, appear to abolish expression of the mismatch repair operon, within which this element resides. The inclusion of further emm4 genomes in our analyses ratified our observations and revealed an international emm4 lineage characterized by prophage degradation. Intriguingly, the USA population of emm4 S. pyogenes appeared to constitute predominantly MEW427-like isolates, whereas the UK population comprised both MEW427-like and MGAS10750-like isolates. The degraded and cryptic nature of these elements may have important phenotypic and fitness ramifications for emm4 S. pyogenes, and the geographical distribution of this lineage raises interesting questions on the population dynamics of the genotype.

In silico characterisation of stand-alone response regulators of Streptococcus pyogenes

Bacterial “stand-alone” response regulators (RRs) are pivotal to the control of gene transcription in response to changing cytosolic and extracellular microenvironments during infection. The genome of group A Streptococcus (GAS) encodes more than 30 stand-alone RRs that orchestrate the expression of virulence factors involved in infecting multiple tissues, so causing an array of potentially lethal human diseases. Here, we analysed the molecular epidemiology and biological associations in the coding sequences (CDSs) and upstream intergenic regions (IGRs) of 35 stand-alone RRs from a collection of global GAS genomes. Of the 944 genomes analysed, 97% encoded 32 or more of the 35 tested RRs. The length of RR CDSs ranged from 297 to 1587 nucleotides with an average nucleotide diversity (π) of 0.012, while the IGRs ranged from 51 to 666 nucleotides with average π of 0.017. We present new evidence of recombination in multiple RRs including mga, leading to mga-2 switching, emm-switching and emm-like gene chimerization, and the first instance of an isolate that encodes both mga-1 and mga-2. Recombination was also evident in rofA/nra and msmR loci with 15 emm-types represented in multiple FCT (fibronectin-binding, collagen-binding, T-antigen)-types, including novel emm-type/FCT-type pairings. Strong associations were observed between concatenated RR allele types, and emm-type, MLST-type, core genome phylogroup, and country of sampling. No strong associations were observed between individual loci and disease outcome. We propose that 11 RRs may form part of future refinement of GAS typing systems that reflect core genome evolutionary associations. This subgenomic analysis revealed allelic traits that were informative to the biological function, GAS strain definition, and regional outbreak detection.

Analysis of Global Collection of Group A Streptococcus Genomes Reveals that the Majority Encode a Trio of M and M-Like Proteins

While the GAS M protein has been the leading vaccine target for decades, the bacteria encode many other virulence factors of interest for vaccine development. In this work, we show that emm-like genes are encoded in a remarkable majority of GAS genomes and expressed at a level similar to that for the emm gene. In collaboration with the U.S. Centers for Disease Control, we developed molecular definitions of the different emm and emm-like gene families. This clarification should abrogate mistyping of strains, especially in the area of whole-genome typing. We have also updated the emm-typing collection by removing emm-like gene sequences and provided in-depth analysis of Mrp and Enn protein sequence structure and diversity.

The core Mga (multiple gene activator) regulon of group A Streptococcus (GAS) contains genes encoding proteins involved in adhesion and immune evasion. While all GAS genomes contain genes for Mga and C5a peptidase, the intervening genes encoding M and M-like proteins vary between strains. The genetic make-up of the Mga regulon of GAS was characterized by utilizing a collection of 1,688 GAS genomes that are representative of the global GAS population. Sequence variations were examined with multiple alignments, and the expression of all core Mga regulon genes was examined by quantitative reverse transcription-PCR in a representative strain collection. In 85.2% of the sampled genomes, the Mga locus contained genes encoding Mga, Mrp, M, Enn, and C5a peptidase proteins. These isolates account for 53% of global infections. Only 9.1% of genomes did not contain either an mrp or an enn gene. The pairwise identity within Enn (68.6%) and Mrp (83.2%) protein sequences was higher than within M proteins (44.7%). Gene expression varied between strains tested, but high expression was recorded for all genes in at least one strain. Previous nomenclature issues were clarified with molecular gene definitions. Our findings support a shift in focus in the GAS research field to further consider the role of Mrp and Enn in virulence and vaccine development.

IMPORTANCE While the GAS M protein has been the leading vaccine target for decades, the bacteria encode many other virulence factors of interest for vaccine development. In this work, we show that emm-like genes are encoded in a remarkable majority of GAS genomes and expressed at a level similar to that for the emm gene. In collaboration with the U.S. Centers for Disease Control, we developed molecular definitions of the different emm and emm-like gene families. This clarification should abrogate mistyping of strains, especially in the area of whole-genome typing. We have also updated the emm-typing collection by removing emm-like gene sequences and provided in-depth analysis of Mrp and Enn protein sequence structure and diversity.

The Emergence of Successful Streptococcus pyogenes Lineages through Convergent Pathways of Capsule Loss and Recombination Directing High Toxin Expression

Gene transfer and homologous recombination in Streptococcus pyogenes has the potential to trigger the emergence of pandemic lineages, as exemplified by lineages of emm1 and emm89 that emerged in the 1980s and 2000s, respectively. Although near-identical replacement gene transfer events in the nga (NADase) and slo (streptolysin O) loci conferring high expression of these toxins underpinned the success of these lineages, extension to other emm genotype lineages is unreported. The emergent emm89 lineage was characterized by five regions of homologous recombination additional to nga-slo, including complete loss of the hyaluronic acid capsule synthesis locus hasABC, a genetic trait replicated in two other leading emm types and recapitulated by other emm types by inactivating mutations. We hypothesized that other leading genotypes may have undergone similar recombination events. We analyzed a longitudinal data set of genomes from 344 clinical invasive disease isolates representative of locations across England, dating from 2001 to 2011, and an international collection of S. pyogenes genomes representing 54 different genotypes and found frequent evidence of recombination events at the nga-slo locus predicted to confer higher toxin genotype. We identified multiple associations between recombination at this locus and inactivating mutations within hasAB, suggesting convergent evolutionary pathways in successful genotypes. This included common genotypes emm28 and emm87. The combination of no or low capsule and high expression of nga and slo may underpin the success of many emergent S. pyogenes lineages of different genotypes, triggering new pandemics, and could change the way S. pyogenes causes disease.IMPORTANCEStreptococcus pyogenes is a genetically diverse pathogen, with over 200 different genotypes defined by emm typing, but only a minority of these genotypes are responsible for the majority of human infection in high-income countries. Two prevalent genotypes associated with disease rose to international dominance following recombination of a toxin locus that conferred increased expression. Here, we found that recombination of this locus and promoter has occurred in other diverse genotypes, events that may allow these genotypes to expand in the population. We identified an association between the loss of hyaluronic acid capsule synthesis and high toxin expression, which we propose may be associated with an adaptive advantage. As S. pyogenes pathogenesis depends both on capsule and toxin production, new variants with altered expression may result in abrupt changes in the molecular epidemiology of this pathogen in the human population over time.

Atlas of group A streptococcal vaccine candidates compiled using large-scale comparative genomics

Group A Streptococcus (GAS; Streptococcus pyogenes) is a bacterial pathogen for which a commercial vaccine for humans is not available. Employing the advantages of high-throughput DNA sequencing technology to vaccine design, we have analyzed 2,083 globally sampled GAS genomes. The global GAS population structure reveals extensive genomic heterogeneity driven by homologous recombination and overlaid with high levels of accessory gene plasticity. We identified the existence of more than 290 clinically associated genomic phylogroups across 22 countries, highlighting challenges in designing vaccines of global utility. To determine vaccine candidate coverage, we investigated all of the previously described GAS candidate antigens for gene carriage and gene sequence heterogeneity. Only 15 of 28 vaccine antigen candidates were found to have both low naturally occurring sequence variation and high (>99%) coverage across this diverse GAS population. This technological platform for vaccine coverage determination is equally applicable to prospective GAS vaccine antigens identified in future studies.

Molecular Epidemiology, Ecology, and Evolution of Group A Streptococci

The clinico-epidemiological features of diseases caused by group A streptococci (GAS) is presented through the lens of the ecology, population genetics, and evolution of the organism. The serological targets of three typing schemes (M, T, SOF) are themselves GAS cell surface proteins that have a myriad of virulence functions and a diverse array of structural forms. Horizontal gene transfer expands the GAS antigenic cell surface repertoire by generating numerous combinations of M, T, and SOF antigens. However, horizontal gene transfer of the serotype determinant genes is not unconstrained, and therein lies a genetic organization that may signify adaptations to a narrow ecological niche, such as the primary tissue reservoirs of the human host. Adaptations may be further shaped by selection pressures such as herd immunity. Understanding the molecular evolution of GAS on multiple levels-short, intermediate, and long term-sheds insight on mechanisms of host-pathogen interactions, the emergence and spread of new clones, rational vaccine design, and public health interventions.

Methods for the detection and characterization of Streptococcus suis: from conventional bacterial culture methods to immunosensors

One of the most important zoonotic pathogens worldwide, Streptococcus suis is a swine pathogen that is responsible for meningitis, toxic shock and even death in humans. S. suis infection develops rapidly with nonspecific clinical symptoms in the early stages and a high fatality rate. Recently, much attention has been paid to the high prevalence of S. suis as well as the increasing incidence and its epidemic characteristics. As laboratory-acquired infections of S. suis can occur and it is dangerous to public health security, timely and early diagnosis has become key to controlling S. suis prevalence. Here, the techniques that have been used for the detection, typing and characterization of S. suis are reviewed and the prospects for future detection methods for this bacterium are also discussed.

In silico characterisation of the two-component system regulators of Streptococcus pyogenes

Bacteria respond to environmental changes through the co-ordinated regulation of gene expression, often mediated by two-component regulatory systems (TCS). Group A Streptococcus (GAS), a bacterium which infects multiple human body sites and causes multiple diseases, possesses up to 14 TCS. In this study we examined genetic variation in the coding sequences and non-coding DNA upstream of these TCS as a method for evaluating relationships between different GAS emm-types, and potential associations with GAS disease. Twelve of the 14 TCS were present in 90% of the genomes examined. The length of the intergenic regions (IGRs) upstream of TCS coding regions varied from 39 to 345 nucleotides, with an average nucleotide diversity of 0.0064. Overall, IGR allelic variation was generally conserved with an emm-type. Subsequent phylogenetic analysis of concatenated sequences based on all TCS IGR sequences grouped genomes of the same emm-type together. However grouping with emm-pattern and emm-cluster-types was much weaker, suggesting epidemiological and functional properties associated with the latter are not due to evolutionary relatedness of emm-types. All emm5, emm6 and most of the emm18 genomes, all historically considered rheumatogenic emm-types clustered together, suggesting a shared evolutionary history. However emm1, emm3 and several emm18 genomes did not cluster within this group. These latter emm18 isolates were epidemiologically distinct from other emm18 genomes in study, providing evidence for local variation. emm-types associated with invasive disease or nephritogenicity also did not cluster together. Considering the TCS coding sequences (cds), correlation with emm-type was weaker than for the IGRs, and no strong correlation with disease was observed. Deletion of the malate transporter, maeP, was identified that serves as a putative marker for the emm89.0 subtype, which has been implicated in invasive outbreaks. A recombination-related, subclade-forming DNA motif was identified in the putative receiver domain of the Spy1556 response regulator that correlated with throat-associated emm-pattern-type A-C strains.

An Outbreak of Streptococcus pyogenes in a Mental Health Facility: Advantage of Well-Timed Whole-Genome Sequencing Over emm Typing

OBJECTIVEWe report the utility of whole-genome sequencing (WGS) conducted in a clinically relevant time frame (ie, sufficient for guiding management decision), in managing a Streptococcus pyogenes outbreak, and present a comparison of its performance with emm typing.SETTINGA 2,000-bed tertiary-care psychiatric hospital.METHODSActive surveillance was conducted to identify new cases of S. pyogenes. WGS guided targeted epidemiological investigations, and infection control measures were implemented. Single-nucleotide polymorphism (SNP)-based genome phylogeny, emm typing, and multilocus sequence typing (MLST) were performed. We compared the ability of WGS and emm typing to correctly identify person-to-person transmission and to guide the management of the outbreak.RESULTSThe study included 204 patients and 152 staff. We identified 35 patients and 2 staff members with S. pyogenes. WGS revealed polyclonal S. pyogenes infections with 3 genetically distinct phylogenetic clusters (C1-C3). Cluster C1 isolates were all emm type 4, sequence type 915 and had pairwise SNP differences of 0-5, which suggested recent person-to-person transmissions. Epidemiological investigation revealed that cluster C1 was mediated by dermal colonization and transmission of S. pyogenes in a male residential ward. Clusters C2 and C3 were genomically diverse, with pairwise SNP differences of 21-45 and 26-58, and emm 11 and mostly emm120, respectively. Clusters C2 and C3, which may have been considered person-to-person transmissions by emm typing, were shown by WGS to be unlikely by integrating pairwise SNP differences with epidemiology.CONCLUSIONSWGS had higher resolution than emm typing in identifying clusters with recent and ongoing person-to-person transmissions, which allowed implementation of targeted intervention to control the outbreak.Infect Control Hosp Epidemiol 2018;852-860.

Prevalent emm Types among Invasive GAS in Europe and North America since Year 2000

BACKGROUND

Streptococcus pyogenes or group A streptococcus (GAS) is an important human pathogen responsible for a broad range of infections, from uncomplicated to more severe and invasive diseases with high mortality and morbidity. Epidemiological surveillance has been crucial to detect changes in the geographical and temporal variation of the disease pattern; for this purpose the M protein gene (emm) gene typing is the most widely used genotyping method, with more than 200 emm types recognized. Molecular epidemiological data have been also used for the development of GAS M protein-based vaccines.

METHODS

The aim of this paper was to provide an updated scenario of the most prevalent GAS emm types responsible for invasive infections in developed countries as Europe and North America (US and Canada), from 1st January 2000 to 31st May 2017. The search, performed in PubMed by the combined use of the terms ("emm") and ("invasive") retrieved 264 articles, of which 38 articles (31 from Europe and 7 from North America) met the inclusion criteria and were selected for this study. Additional five papers cited in the European articles but not retrieved by the search were included.

RESULTS

emm1 represented the dominant type in both Europe and North America, replaced by other emm types in only few occasions. The seven major emm types identified (emm1, emm28, emm89, emm3, emm12, emm4, and emm6) accounted for approximately 50-70% of the total isolates; less common emm types accounted for the remaining 30-50% of the cases. Most of the common emm types are included in either one or both the 26-valent and 30-valent vaccines, though some well-represented emm types found in Europe are not.

CONCLUSION

This study provided a picture of the prevalent emm types among invasive GAS (iGAS) in Europe and North America since the year 2000 onward. Continuous surveillance on the emm-type distribution among iGAS infections is strongly encouraged also to determine the potential coverage of the developing multivalent vaccines.

Virulence Role of the GlcNAc Side Chain of the Lancefield Cell Wall Carbohydrate Antigen in Non-M1-Serotype Group A Streptococcus

Classification of streptococci is based upon expression of unique cell wall carbohydrate antigens. All serotypes of group A Streptococcus (GAS; Streptococcus pyogenes), a leading cause of infection-related mortality worldwide, express the group A carbohydrate (GAC). GAC, the classical Lancefield antigen, is comprised of a polyrhamnose backbone with N-acetylglucosamine (GlcNAc) side chains. The immunodominant GlcNAc epitope of GAC is the basis of all rapid diagnostic testing for GAS infection. We previously identified the 12-gene GAC biosynthesis gene cluster and determined that the glycosyltransferase GacI was required for addition of the GlcNAc side chain to the polyrhamnose core. Loss of the GAC GlcNAc epitope in serotype M1 GAS resulted in attenuated virulence in two animal infection models and increased GAS sensitivity to killing by whole human blood, serum, neutrophils, and antimicrobial peptides. Here, we report that the GAC biosynthesis gene cluster is ubiquitous among 520 GAS isolates from global sources, representing 105 GAS emm serotypes. Isogenic ΔgacI mutants were constructed in M2, M3, M4, M28, and M89 backgrounds and displayed an array of phenotypes in susceptibility to killing by whole human blood, baby rabbit serum, human platelet releasate, human neutrophils, and antimicrobial peptide LL-37. The contribution of the GlcNAc side chain to GAS survival in vivo also varied by strain, demonstrating that it is not a prerequisite for virulence in the murine infection model. Thus, the relative contribution of GAC to virulence in non-M1 serotypes appears to depend on the quorum of other virulence factors that each strain possesses.IMPORTANCE The Lancefield group A carbohydrate (GAC) is the species-defining antigen for group A Streptococcus (GAS), comprising ~50% of the cell wall of this major human pathogen. We previously showed that the GlcNAc side chain of GAC contributes to the innate immune resistance and animal virulence phenotypes of the globally disseminated strain of serotype M1 GAS. Here, we use isogenic mutagenesis to examine the role of GAC GlcNAc in five additional medically relevant GAS serotypes. Overall, the GlcNAc side chain of GAC contributes to the innate immune resistance of GAS, but the relative contribution varies among individual strains. Moreover, the GAC GlcNAc side chain is not a universal prerequisite for GAS virulence in the animal model.

Integration of Genomic and Other Epidemiologic Data to Investigate and Control a Cross-Institutional Outbreak of Streptococcus pyogenes

Single-strain outbreaks of Streptococcus pyogenes infections are common and often go undetected. In 2013, two clusters of invasive group A Streptococcus (iGAS) infection were identified in independent but closely located care homes in Oxfordshire, United Kingdom. Investigation included visits to each home, chart review, staff survey, microbiologic sampling, and genome sequencing. S. pyogenes emm type 1.0, the most common circulating type nationally, was identified from all cases yielding GAS isolates. A tailored whole-genome reference population comprising epidemiologically relevant contemporaneous isolates and published isolates was assembled. Data were analyzed independently using whole-genome multilocus sequencing and single-nucleotide polymorphism analyses. Six isolates from staff and residents of the homes formed a single cluster that was separated from the reference population by both analytical approaches. No further cases occurred after mass chemoprophylaxis and enhanced infection control. Our findings demonstrate the ability of 2 independent analytical approaches to enable robust conclusions from nonstandardized whole-genome analysis to support public health practice.

Management of children with persistent group A streptococcal carriage

INTRODUCTION

Chronic GAS carrier state is best defined as the prolonged presence of group A β-haemolytic Streptococcus (GAS) in the pharynx without evidence of infection or inflammation. Chronic GAS carriers have a low risk of immune mediated complications. Persistent pharyngeal carriage often raises management issues. In this study, we review the evidence on the management of persistent GAS carriage in children and propose a management algorithm. Areas covered: Chronic GAS pharyngeal carriage is quite common affecting 10-20% of school-aged children. Pathogenesis of carriage has been related to the pharynx microflora and to special properties of GAS, but several aspects are yet to be elucidated. Management greatly depends on whether the individual child belongs to a 'high-risk' group and might benefit from eradication regimens or not, when observation-only and reassurance are enough. Penicillin plus rifampin and clindamycin monotherapy have been recommended for eradication; limited evidence of effectiveness of azithromycin has been reported. Surgical intervention is not indicated. Expert commentary: GAS infection is a common reason for antibiotic use and abuse in children and asymptomatic carriers constitute the major reservoir of GAS in the community. Several aspects are yet to be elucidated and well-designed studies are needed for firm conclusions to be drawn.

Whole genome sequencing of group A Streptococcus: development and evaluation of an automated pipeline for emmgene typing

Streptococcus pyogenes group A Streptococcus (GAS) is the most common cause of bacterial throat infections, and can cause mild to severe skin and soft tissue infections, including impetigo, erysipelas, necrotizing fasciitis, as well as systemic and fatal infections including septicaemia and meningitis. Estimated annual incidence for invasive group A streptococcal infection (iGAS) in industrialised countries is approximately three per 100,000 per year. Typing is currently used in England and Wales to monitor bacterial strains of S. pyogenes causing invasive infections and those isolated from patients and healthcare/care workers in cluster and outbreak situations. Sequence analysis of the emm gene is the currently accepted gold standard methodology for GAS typing. A comprehensive database of emm types observed from superficial and invasive GAS strains from England and Wales informs outbreak control teams during investigations. Each year the Bacterial Reference Department, Public Health England (PHE) receives approximately 3,000 GAS isolates from England and Wales. In April 2014 the Bacterial Reference Department, PHE began genomic sequencing of referred S. pyogenes isolates and those pertaining to selected elderly/nursing care or maternity clusters from 2010 to inform future reference services and outbreak analysis (n = 3, 047). In line with the modernizing strategy of PHE, we developed a novel bioinformatics pipeline that can predict emmtypes using whole genome sequence (WGS) data. The efficiency of this method was measured by comparing the emmtype assigned by this method against the result from the current gold standard methodology; concordance to emmsubtype level was observed in 93.8% (2,852/3,040) of our cases, whereas in 2.4% (n = 72) of our cases concordance was observed to emm type level. The remaining 3.8% (n = 117) of our cases corresponded to novel types/subtypes, contamination, laboratory sample transcription errors or problems arising from high sequence similarity of the allele sequence or low mapping coverage. De novo assembly analysis was performed in the two latter groups (n = 72 + 117) and was able to diagnose the problem and where possible resolve the discordance (60/72 and 20/117, respectively). Overall, we have demonstrated that our WGS emm-typing pipeline is a reliable and robust system that can be implemented to determine emm type for the routine service.

Epidemiological and Molecular Characterization of an Invasive Group A Streptococcus emm32.2 Outbreak

An emm32.2 invasive group A streptococcus (iGAS) outbreak occurred in Liverpool from January 2010 to September 2012. This genotype had not previously been identified in Liverpool, but was responsible for 32% (14/44) of all iGAS cases reported during this time period. We performed a case-case comparison of emm32.2 iGAS cases with non-emm32.2 control iGAS cases identified in the Liverpool population over the same time period to assess patient risk factors for emm32.2 iGAS infection. The emm32.2 iGAS cases were confined to the adult population. We show that homelessness, intravenous drug use, and alcohol abuse predisposed patients to emm32.2 iGAS disease; however, no obvious epidemiological linkage between the patients with emm32.2 iGAS could be identified. Comparative whole-genome sequencing analysis of emm32.2 iGAS and non-emm32.2 control isolates was also performed to identify pathogen factors which might have driven the outbreak. We identified 19 genes, five of which had previously been implicated in virulence, which were present in all of the emm32.2 iGAS isolates but not present in any of the non-emm32.2 control isolates. We report that a novel emm32.2 genotype emerged in Liverpool in 2010 and identified a specific subset of genes, which could have allowed this novel emm32.2 genotype to persist in a disadvantaged population in the region over a 3-year period.

Whole-genome Sequencing for Surveillance of Invasive Pneumococcal Diseases in Ontario, Canada: Rapid Prediction of Genotype, Antibiotic Resistance and Characterization of Emerging Serotype 22F

Background: Molecular typing is essential for inferring genetic relatedness between bacterial pathogens. In this study, we applied whole genome sequencing (WGS) for rapid prediction of sequence type and antibiotic resistance for invasive pneumococcal isolates. Methods: 240 isolates from adults (≥50 years old) in Ontario, Canada during 2009 to 2013 were subjected to WGS. Sequence type, antibiotic susceptibility and resistance were predicted directly from short reads. Emerging non-vaccine serotype 22F was further characterized by WGS. Results: Sequence type was successfully determined for 98.3% of isolates. The overall sensitivity and specificity for antibiotic resistance prediction were 95 and 100% respectively, compared to standard susceptibility testing methods. WGS-based phylogeny divided emerging 22F (ST433) strains into two distinct clades: clade A harboring a 23 kb-prophage and anti-phage PhD/Doc system and clade B with virulence-related proteases. Five isolates in clade A developed macrolide resistance via 5.1 kb mega element recombination (encoding mefE and msrD), while one isolate in clade B displayed quinolone resistance via a gyrA mutation. Conclusions: WGS is valuable for routine surveillance of pneumococcal clinical isolates and facilitates prediction of genotype and antibiotic resistance. The emergence of 22F in Ontario in the post-vaccine era and evidence of evolution and divergence of the 22F population warrants heightened pneumococcal molecular surveillance.

International circumpolar surveillance interlaboratory quality control program for emm typing of Streptococcus pyogenes, 2011-2015

In 2011, an interlaboratory quality control (QC) program for emm typing group A streptococci (GAS) was incorporated into existing international circumpolar surveillance QC programs. From 2011 - 2015, 35 GAS isolates were distributed to three laboratories; emm type-level concordance was 100%, while the overall sub-type level concordance was 83%.

Determining Streptococcus suis serotype from short-read whole-genome sequencing data

Background

Streptococcus suis is divided into 29 serotypes based on a serological reaction against the capsular polysaccharide (CPS). Multiplex PCR tests targeting the cps locus are also used to determine S. suis serotypes, but they cannot differentiate between serotypes 1 and 14, and between serotypes 2 and 1/2. Here, we developed a pipeline permitting in silico serotype determination from whole-genome sequencing (WGS) short-read data that can readily identify all 29 S. suis serotypes.

Results

We sequenced the genomes of 121 strains representing all 29 known S. suis serotypes. We next combined available software into an automated pipeline permitting in silico serotyping of strains by differential alignment of short-read sequencing data to a custom S. suis cps loci database. Strains of serotype pairs 1 and 14, and 2 and 1/2 could be differentiated by a missense mutation in the cpsK gene. We report a 99 % match between coagglutination- and pipeline-determined serotypes for strains in our collection. We used 375 additional S. suis genomes downloaded from the NCBI’s Sequence Read Archive (SRA) to validate the pipeline. Validation with SRA WGS data resulted in a 92 % match. Included pipeline subroutines permitted us to assess strain virulence marker content and obtain multilocus sequence typing directly from WGS data.

Conclusions

Our pipeline permits rapid and accurate determination of S. suis serotype, and other lineage information, directly from WGS data. By discriminating between serotypes 1 and 14, and between serotypes 2 and 1/2, our approach solves a three-decade longstanding S. suis typing issue.

Electronic supplementary material

The online version of this article (doi:10.1186/s12866-016-0782-8) contains supplementary material, which is available to authorized users.

Whole genome sequencing in clinical and public health microbiology

SummaryGenomics and whole genome sequencing (WGS) have the capacity to greatly enhance knowledge and understanding of infectious diseases and clinical microbiology.

The growth and availability of bench-top WGS analysers has facilitated the feasibility of genomics in clinical and public health microbiology.

Given current resource and infrastructure limitations, WGS is most applicable to use in public health laboratories, reference laboratories, and hospital infection control-affiliated laboratories.

As WGS represents the pinnacle for strain characterisation and epidemiological analyses, it is likely to replace traditional typing methods, resistance gene detection and other sequence-based investigations (e.g., 16S rDNA PCR) in the near future.

Although genomic technologies are rapidly evolving, widespread implementation in clinical and public health microbiology laboratories is limited by the need for effective semi-automated pipelines, standardised quality control and data interpretation, bioinformatics expertise, and infrastructure.

Epidemiology of Invasive Group A Streptococcal Disease in Alaska, 2001 to 2013

The Arctic Investigations Program (AIP) began surveillance for invasive group A streptococcal (GAS) infections in Alaska in 2000 as part of the invasive bacterial diseases population-based laboratory surveillance program. Between 2001 and 2013, there were 516 cases of GAS infection reported, for an overall annual incidence of 5.8 cases per 100,000 persons with 56 deaths (case fatality rate, 10.7%). Of the 516 confirmed cases of invasive GAS infection, 422 (82%) had isolates available for laboratory analysis. All isolates were susceptible to penicillin, cefotaxime, and levofloxacin. Resistance to tetracycline, erythromycin, and clindamycin was seen in 11% (n = 8), 5.8% (n = 20), and 1.2% (n = 4) of the isolates, respectively. A total of 51 emm types were identified, of which emm1 (11.1%) was the most prevalent, followed by emm82 (8.8%), emm49 (7.8%), emm12 and emm3 (6.6% each), emm89 (6.2%), emm108 (5.5%), emm28 (4.7%), emm92 (4%), and emm41 (3.8%). The five most common emm types accounted for 41% of isolates. The emm types in the proposed 26-valent and 30-valent vaccines accounted for 56% and 78% of all cases, respectively. GAS remains an important cause of invasive bacterial disease in Alaska. Continued surveillance of GAS infections will help improve understanding of the epidemiology of invasive disease, with an impact on disease control, notification of outbreaks, and vaccine development.

SRST2: Rapid genomic surveillance for public health and hospital microbiology labs

Rapid molecular typing of bacterial pathogens is critical for public health epidemiology, surveillance and infection control, yet routine use of whole genome sequencing (WGS) for these purposes poses significant challenges. Here we present SRST2, a read mapping-based tool for fast and accurate detection of genes, alleles and multi-locus sequence types (MLST) from WGS data. Using >900 genomes from common pathogens, we show SRST2 is highly accurate and outperforms assembly-based methods in terms of both gene detection and allele assignment. We include validation of SRST2 within a public health laboratory, and demonstrate its use for microbial genome surveillance in the hospital setting. In the face of rising threats of antimicrobial resistance and emerging virulence among bacterial pathogens, SRST2 represents a powerful tool for rapidly extracting clinically useful information from raw WGS data. Source code is available from http://katholt.github.io/srst2/.

SRST2: Rapid genomic surveillance for public health and hospital microbiology labs

Rapid molecular typing of bacterial pathogens is critical for public health epidemiology, surveillance and infection control, yet routine use of whole genome sequencing (WGS) for these purposes poses significant challenges. Here we present SRST2, a read mapping-based tool for fast and accurate detection of genes, alleles and multi-locus sequence types (MLST) from WGS data. Using >900 genomes from common pathogens, we show SRST2 is highly accurate and outperforms assembly-based methods in terms of both gene detection and allele assignment. We include validation of SRST2 within a public health laboratory, and demonstrate its use for microbial genome surveillance in the hospital setting. In the face of rising threats of antimicrobial resistance and emerging virulence among bacterial pathogens, SRST2 represents a powerful tool for rapidly extracting clinically useful information from raw WGS data. Source code is available from http://katholt.github.io/srst2/.

Clinical laboratory response to a mock outbreak of invasive bacterial infections: a preparedness study

Large hospital-based clinical laboratories must be prepared to rapidly investigate potential infectious disease outbreaks. To challenge the ability of our molecular diagnostics laboratory to use whole-genome sequencing in a potential outbreak scenario and identify impediments to these efforts, we studied 84 invasive serotype emm59 group A streptococcus (GAS) strains collected in the United States. We performed a rapid-response exercise to the mock outbreak scenario using whole-genome sequencing, genome-wide transcript analysis, and mouse virulence studies. The protocol changes installed in response to the lessons learned were tested in a second iteration. The initial investigation was completed in 9 days. Whole-genome sequencing showed that the invasive infections were caused by multiple subclones of epidemic emm59 GAS strains likely spread to the United States from Canada. The phylogenetic tree showed a strong temporal-spatial structure with diversity in mobile genetic element content, features that are useful for identifying closely related strains and possible transmission events. The genome data informed the epidemiology, identifying multiple patients who likely acquired the organisms through direct person-to-person transmission. Transcriptome analysis unexpectedly revealed significantly altered expression of genes encoding a two-component regulator and the hyaluronic acid capsule virulence factor. Mouse infection studies confirmed a high-virulence capacity of these emm59 organisms. Whole-genome sequencing, coupled with transcriptome analysis and animal virulence studies, can be rapidly performed in a clinical environment to effectively contribute to patient care decisions and public health maneuvers.