Identification of a chimeric emm gene and novel emm pattern in currently circulating strains of emm4 Group A Streptococcus

GAS-J, a User-Friendly Browser Application for Genome Assembly, emm-Typing, MLST Typing, and Virulence Factor Gene Detection of Streptococcus pyogenes

Clinical isolates of Streptococcus pyogenes are usually classified using emm and multilocus sequence typing (MLST). Recently, whole genome sequencing (WGS) has been employed for emm typing and MLST analysis. WGS data provides additional information on the presence of virulence factor genes. To enable researchers unfamiliar with bioinformatics to analyze WGS data of S. pyogenes, we opened an online tool named GAS-J, which automatically outputs emm types, sequence types (STs), carriage of virulence factor genes, and phylogenetic trees. The tool accepts raw short-read data as inputs, since it includes the velvet assembler. In all isolates, the emm typing results from this tool were identical to those obtained by conventional PCR and Sanger sequencing, even in cases where isolates had pseudo-emm (emm-like) genes. STs are determined by performing a BLAST search using seven alleles as references. To detect S. pyogenes virulence factor genes, we prepared a new data set containing 620 related proteins. Users may choose which isolates to include in SNP-based phylogenetic tree from a pool of 406 isolates with epidemiological data. The data set includes isolates whose symptoms (STSS or non-STSS) were diagnosed based on the STSS criteria of the Japan Communicable Disease Prevention Law. GAS-J application is available at http://gasj.ncgm.go.jp. The data of isolates are going to be updated in the future.

Emergent emm4 group A Streptococcus evidences a survival strategy during interaction with immune effector cells

The major gram-positive pathogen group A Streptococcus (GAS) is a model organism for studying microbial epidemics as it causes waves of infections. Since 1980, several GAS epidemics have been ascribed to the emergence of clones producing increased amounts of key virulence factors such as streptolysin O (SLO). Herein, we sought to identify mechanisms underlying our recently identified temporal clonal emergence among emm4 GAS, given that emergent strains did not produce augmented levels of virulence factors relative to historic isolates. By creating and analyzing isoallelic strains, we determined that a conserved mutation in a previously undescribed gene encoding a putative carbonic anhydrase was responsible for the defective in vitro growth observed in the emergent strains. We also identified that the emergent strains survived better inside macrophages and killed macrophages at lower rates than the historic strains. Via the creation of isogenic mutant strains, we linked the emergent strain "survival" phenotype to the downregulation of the SLO encoding gene and upregulation of the msrAB operon which encodes proteins involved in defense against extracellular oxidative stress. Our findings are in accord with recent surveillance studies which found a high ratio of mucosal (i.e., pharyngeal) relative to invasive infections among emm4 GAS. Since ever-increasing virulence is unlikely to be evolutionarily advantageous for a microbial pathogen, our data further understanding of the well-described oscillating patterns of virulent GAS infections by demonstrating mechanisms by which emergent strains adapt a "survival" strategy to outcompete previously circulating isolates.

Presence of Group A streptococcus frequently assayed virulence genes in invasive disease: a systematic review and meta-analysis

Introduction

It is currently unclear what the role of Group A streptococcus (GAS) virulence factors (VFs) is in contributing to the invasive potential of GAS. This work investigated the evidence for the association of GAS VFs with invasive disease.

Methods

We employed a broad search strategy for studies reporting the presence of GAS VFs in invasive and non-invasive GAS disease. Data were independently extracted by two reviewers, quality assessed, and meta-analyzed using Stata®.

Results

A total of 32 studies reported on 45 putative virulence factors [invasive (n = 3,236); non-invasive (n = 5,218)], characterized by polymerase chain reaction (PCR) (n = 30) and whole-genome sequencing (WGS) (n = 2). The risk of bias was rated as low and moderate, in 23 and 9 studies, respectively. Meta-,analyses of high-quality studies (n = 23) revealed a significant association of speM [OR, 1.64 (95%CI, 1.06; 2.52)] with invasive infection. Meta-analysis of WGS studies demonstrated a significant association of hasA [OR, 1.91 (95%CI, 1.36; 2.67)] and speG [OR, 2.83 (95%CI, 1.63; 4.92)] with invasive GAS (iGAS). Meta-analysis of PCR studies indicated a significant association of speA [OR, 1.59 (95%CI, 1.10; 2.30)] and speK [OR, 2.95 (95%CI, 1.81; 4.80)] with invasive infection. A significant inverse association was observed between prtf1 [OR, 0.42 (95%CI, 0.20; 0.87)] and invasive infection.

Conclusion

This systematic review and genomic meta-analysis provides evidence of a statistically significant association with invasive infection for the hasA gene, while smeZ, ssa, pnga3, sda1, sic, and NaDase show statistically significantly inverse associations with invasive infection. SpeA, speK, and speG are associated with GAS virulence; however, it is unclear if they are markers of invasive infection. This work could possibly aid in developing preventative strategies.

Emergent emm4 group A Streptococcus evidences a survival strategy during interaction with immune effector cells

The major gram-positive pathogen group A Streptococcus (GAS) is a model organism for studying microbial epidemics as it causes waves of infections. Since 1980, several GAS epidemics have been ascribed to the emergence of clones producing increased amounts of key virulence factors such as streptolysin O (SLO). Herein, we sought to identify mechanisms underlying our recently identified temporal clonal emergence amongst emm4 GAS, given that emergent strains did not produce augmented levels of virulence factors relative to historic isolates. Through the creation and analysis of isoallelic strains, we determined that a conserved mutation in a previously undescribed gene encoding a putative carbonic anhydrase was responsible for the defective in vitro growth observed in the emergent strains. We also identified that the emergent strains survived better inside macrophages and killed macrophages at lower rates relative to the historic strains. Via creation of isogenic mutant strains, we linked the emergent strain "survival" phenotype to the downregulation of the SLO encoding gene and upregulation of the msrAB operon which encodes proteins involved in defense against extracellular oxidative stress. Our findings are in accord with recent surveillance studies which found high ratio of mucosal (i.e., pharyngeal) relative to invasive infections amongst emm4 GAS. Inasmuch as ever-increasing virulence is unlikely to be evolutionary advantageous for a microbial pathogen, our data furthers understanding of the well described oscillating patterns of virulent GAS infections by demonstrating mechanisms by which emergent strains adapt a "survival" strategy to outcompete previously circulating isolates.

Genomic Characterization of Skin and Soft Tissue Streptococcus pyogenes Isolates from a Low-Income and a High-Income Setting

Streptococcus pyogenes is a leading cause of human morbidity and mortality, especially in resource-limited settings. The development of a vaccine against S. pyogenes is a global health priority to reduce the burden of postinfection rheumatic heart disease. To support this, molecular characterization of circulating S. pyogenes isolates is needed. We performed whole-genome analyses of S. pyogenes isolates from skin and soft tissue infections in Sukuta, The Gambia, a low-income country (LIC) in West Africa where there is a high burden of such infections. To act as a comparator to these LIC isolates, skin infection isolates from Sheffield, United Kingdom (a high-income country [HIC]), were also sequenced. The LIC isolates from The Gambia were genetically more diverse (46 emm types in 107 isolates) than the HIC isolates from Sheffield (23 emm types in 142 isolates), with only 7 overlapping emm types. Other molecular markers were shared, including a high prevalence of the skin infection-associated emm pattern D and the variable fibronectin-collagen-T antigen (FCT) types FCT-3 and FCT-4. Fewer of the Gambian LIC isolates carried prophage-associated superantigens (64%) and DNases (26%) than did the Sheffield HIC isolates (99% and 95%, respectively). We also identified streptococcin genes unique to 36% of the Gambian LIC isolates and a higher prevalence (48%) of glucuronic acid utilization pathway genes in the Gambian LIC isolates than in the Sheffield HIC isolates (26%). Comparison to a wider collection of HIC and LIC isolate genomes supported our findings of differing emm diversity and prevalence of bacterial factors. Our study provides insight into the genetics of LIC isolates and how they compare to HIC isolates. IMPORTANCE The global burden of rheumatic heart disease (RHD) has triggered a World Health Organization response to drive forward development of a vaccine against the causative human pathogen Streptococcus pyogenes. This burden stems primarily from low- and middle-income settings where there are high levels of S. pyogenes skin and soft tissue infections, which can lead to RHD. Our study provides much needed whole-genome-based molecular characterization of isolates causing skin infections in Sukuta, The Gambia, a low-income country (LIC) in West Africa where infection and RHD rates are high. Although we identified a greater level of diversity in these LIC isolates than in isolates from Sheffield, United Kingdom (a high-income country), there were some shared features. There were also some features that differed by geographical region, warranting further investigation into their contribution to infection. Our study has also contributed data essential for the development of a vaccine that would target geographically relevant strains.

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.

Population Genomics of emm4 Group A Streptococcus Reveals Progressive Replacement with a Hypervirulent Clone in North America

Clonal replacement is a major driver for changes in bacterial disease epidemiology. Recently, it has been proposed that episodic emergence of novel, hypervirulent clones of group A Streptococcus (GAS) results from acquisition of a 36-kb DNA region leading to increased expression of the cytotoxins Nga (NADase) and SLO (streptolysin O). We previously described a gene fusion event involving the gene encoding the GAS M protein (emm) and an adjacent M-like protein (enn) in the emm4 GAS population, a GAS emm type that lacks the hyaluronic acid capsule. Using whole-genome sequencing of a temporally and geographically diverse set of 1,126 isolates, we discovered that the North American emm4 GAS population has undergone clonal replacement with emergent GAS strains completely replacing historical isolates by 2017. Emergent emm4 GAS strains contained a handful of small genetic variations, including the emm-enn gene fusion, and showed a marked in vitro growth defect compared to historical strains. In contrast to other previously described GAS clonal replacement events, emergent emm4 GAS strains were not defined by acquisition of exogenous DNA and had no significant increase in transcript levels of nga and slo toxin genes via RNA sequencing and quantitative real-time PCR analysis relative to historic strains. Despite the in vitro growth differences, emergent emm4 GAS strains were hypervirulent in mice and ex vivo growth in human blood compared to historical strains. Thus, these data detail the emergence and dissemination of a hypervirulent acapsular GAS clone defined by small, endogenous genetic variation, thereby defining a novel model for GAS strain replacement.

IMPORTANCE Severe invasive infections caused by group A Streptococcus (GAS) result in substantial morbidity and mortality in children and adults worldwide. Previously, GAS clonal strain replacement has been attributed to acquisition of exogenous DNA leading to novel virulence gene acquisition or increased virulence gene expression. Our study of type emm4 GAS identified emergence of a hypervirulent GAS clade defined by variation in endogenous DNA content and lacking augmented toxin gene expression relative to replaced strains. These findings expand our understanding of the molecular mechanisms underlying bacterial clonal emergence.

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.

Global genomic epidemiology of Streptococcus pyogenes

Streptococcus pyogenes is one of the Top 10 human infectious disease killers worldwide causing a range of clinical manifestations in humans. Colonizing a range of ecological niches within its sole host, the human, is key to the ability of this opportunistic pathogen to cause direct and post-infectious manifestations. The expansion of genome sequencing capabilities and data availability over the last decade has led to an improved understanding of the evolutionary dynamics of this pathogen within a global framework where epidemiological relationships and evolutionary mechanisms may not be universal. This review uses the recent publication by Davies et al., 2019 as an updated global framework to address S. pyogenes population genomics, highlighting how genomics is being used to gain new insights into evolutionary processes, transmission pathways, and vaccine design.

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.