Waleed Nasser

Evolutionary pathway to increased virulence and epidemic group A Streptococcus disease derived from 3,615 genome sequences

Significance Epidemics of microbial infections are a considerable threat to human and animal health. Analysis of 3,615 genome sequences, coupled with virulence studies in animals, permitted us to delineate the nature and timing of molecular events that contributed to an ongoing global human epidemic of infections caused by group A Streptococcus, the “flesh-eating” pathogen. We clarified decades-long uncertainty about the timing and sequence of genomic alterations that underpinned the global epidemic. Analyses of this type are crucial for developing better strategies to predict and monitor strain emergence and epidemics, formulate effective protective public health maneuvers, and develop or modify vaccines. We sequenced the genomes of 3,615 strains of serotype Emm protein 1 (M1) group A Streptococcus to unravel the nature and timing of molecular events contributing to the emergence, dissemination, and genetic diversification of an unusually virulent clone that now causes epidemic human infections worldwide. We discovered that the contemporary epidemic clone emerged in stepwise fashion from a precursor cell that first contained the phage encoding an extracellular DNase virulence factor (streptococcal DNase D2, SdaD2) and subsequently acquired the phage encoding the SpeA1 variant of the streptococcal pyrogenic exotoxin A superantigen. The SpeA2 toxin variant evolved from SpeA1 by a single-nucleotide change in the M1 progenitor strain before acquisition by horizontal gene transfer of a large chromosomal region encoding secreted toxins NAD+-glycohydrolase and streptolysin O. Acquisition of this 36-kb region in the early 1980s into just one cell containing the phage-encoded sdaD2 and speA2 genes was the final major molecular event preceding the emergence and rapid intercontinental spread of the contemporary epidemic clone. Thus, we resolve a decades-old controversy about the type and sequence of genomic alterations that produced this explosive epidemic. Analysis of comprehensive, population-based contemporary invasive strains from seven countries identified strong patterns of temporal population structure. Compared with a preepidemic reference strain, the contemporary clone is significantly more virulent in nonhuman primate models of pharyngitis and necrotizing fasciitis. A key finding is that the molecular evolutionary events transpiring in just one bacterial cell ultimately have produced millions of human infections worldwide.

Transcriptome Remodeling Contributes to Epidemic Disease Caused by the Human Pathogen Streptococcus pyogenes

ABSTRACT For over a century, a fundamental objective in infection biology research has been to understand the molecular processes contributing to the origin and perpetuation of epidemics. Divergent hypotheses have emerged concerning the extent to which environmental events or pathogen evolution dominates in these processes. Remarkably few studies bear on this important issue. Based on population pathogenomic analysis of 1,200 Streptococcus pyogenes type emm89 infection isolates, we report that a series of horizontal gene transfer events produced a new pathogenic genotype with increased ability to cause infection, leading to an epidemic wave of disease on at least two continents. In the aggregate, these and other genetic changes substantially remodeled the transcriptomes of the evolved progeny, causing extensive differential expression of virulence genes and altered pathogen-host interaction, including enhanced immune evasion. Our findings delineate the precise molecular genetic changes that occurred and enhance our understanding of the evolutionary processes that contribute to the emergence and persistence of epidemically successful pathogen clones. The data have significant implications for understanding bacterial epidemics and for translational research efforts to blunt their detrimental effects. IMPORTANCE The confluence of studies of molecular events underlying pathogen strain emergence, evolutionary genetic processes mediating altered virulence, and epidemics is in its infancy. Although understanding these events is necessary to develop new or improved strategies to protect health, surprisingly few studies have addressed this issue, in particular, at the comprehensive population genomic level. Herein we establish that substantial remodeling of the transcriptome of the human-specific pathogen Streptococcus pyogenes by horizontal gene flow and other evolutionary genetic changes is a central factor in precipitating and perpetuating epidemic disease. The data unambiguously show that the key outcome of these molecular events is evolution of a new, more virulent pathogenic genotype. Our findings provide new understanding of epidemic disease. The confluence of studies of molecular events underlying pathogen strain emergence, evolutionary genetic processes mediating altered virulence, and epidemics is in its infancy. Although understanding these events is necessary to develop new or improved strategies to protect health, surprisingly few studies have addressed this issue, in particular, at the comprehensive population genomic level. Herein we establish that substantial remodeling of the transcriptome of the human-specific pathogen Streptococcus pyogenes by horizontal gene flow and other evolutionary genetic changes is a central factor in precipitating and perpetuating epidemic disease. The data unambiguously show that the key outcome of these molecular events is evolution of a new, more virulent pathogenic genotype. Our findings provide new understanding of epidemic disease.

Trading Capsule for Increased Cytotoxin Production: Contribution to Virulence of a Newly Emerged Clade of emm89 Streptococcus pyogenes

ABSTRACT Strains of emm89 Streptococcus pyogenes have become one of the major causes of invasive infections worldwide in the last 10 years. We recently sequenced the genome of 1,125 emm89 strains and identified three major phylogenetic groups, designated clade 1, clade 2, and the epidemic clade 3. Epidemic clade 3 strains, which now cause the great majority of infections, have two distinct genetic features compared to clade 1 and clade 2 strains. First, all clade 3 organisms have a variant 3 nga promoter region pattern, which is associated with increased production of secreted cytolytic toxins SPN (S. pyogenes NADase) and SLO (streptolysin O). Second, all clade 3 strains lack the hasABC locus mediating hyaluronic acid capsule synthesis, whereas this locus is intact in clade 1 and clade 2 strains. We constructed isogenic mutant strains that produce different levels of SPN and SLO toxins and capsule (none, low, or high). Here we report that emm89 strains with elevated toxin production are significantly more virulent than low-toxin producers. Importantly, we also show that capsule production is dispensable for virulence in strains that already produce high levels of SPN and SLO. Our results provide new understanding about the molecular mechanisms contributing to the rapid emergence and molecular pathogenesis of epidemic clade 3 emm89 S. pyogenes. IMPORTANCE S. pyogenes (group A streptococcus [GAS]) causes pharyngitis (“strep throat”), necrotizing fasciitis, and other human infections. Serious infections caused by emm89 S. pyogenes strains have recently increased in frequency in many countries. Based on whole-genome sequence analysis of 1,125 strains recovered from patients on two continents, we discovered that a new emm89 clone, termed clade 3, has two distinct genetic features compared to its predecessors: (i) absence of the genes encoding antiphagocytic hyaluronic acid capsule virulence factor and (ii) increased production of the secreted cytolytic toxins SPN and SLO. emm89 S. pyogenes strains with the clade 3 phenotype (absence of capsule and high expression of SPN and SLO) are highly virulent in mice. These findings provide new understanding of how new virulent clones emerge and cause severe infections worldwide. This newfound knowledge of S. pyogenes virulence can be used to help understand future epidemics and conduct new translational research. S. pyogenes (group A streptococcus [GAS]) causes pharyngitis (“strep throat”), necrotizing fasciitis, and other human infections. Serious infections caused by emm89 S. pyogenes strains have recently increased in frequency in many countries. Based on whole-genome sequence analysis of 1,125 strains recovered from patients on two continents, we discovered that a new emm89 clone, termed clade 3, has two distinct genetic features compared to its predecessors: (i) absence of the genes encoding antiphagocytic hyaluronic acid capsule virulence factor and (ii) increased production of the secreted cytolytic toxins SPN and SLO. emm89 S. pyogenes strains with the clade 3 phenotype (absence of capsule and high expression of SPN and SLO) are highly virulent in mice. These findings provide new understanding of how new virulent clones emerge and cause severe infections worldwide. This newfound knowledge of S. pyogenes virulence can be used to help understand future epidemics and conduct new translational research.

Bacterial Discrimination by Dictyostelid Amoebae Reveals the Complexity of Ancient Interspecies Interactions

BACKGROUND Amoebae and bacteria interact within predator-prey and host-pathogen relationships, but the general response of amoeba to bacteria is not well understood. The amoeba Dictyostelium discoideum feeds on, and is colonized by, diverse bacterial species, including Gram-positive [Gram(+)] and Gram-negative [Gram(-)] bacteria, two major groups of bacteria that differ in structure and macromolecular composition. RESULTS Transcriptional profiling of D. discoideum revealed sets of genes whose expression is enriched in amoebae interacting with different species of bacteria, including sets that appear specific to amoebae interacting with Gram(+) or with Gram(-) bacteria. In a genetic screen utilizing the growth of mutant amoebae on a variety of bacteria as a phenotypic readout, we identified amoebal genes that are only required for growth on Gram(+) bacteria, including one that encodes the cell-surface protein gp130, as well as several genes that are only required for growth on Gram(-) bacteria, including one that encodes a putative lysozyme, AlyL. These genes are required for parts of the transcriptional response of wild-type amoebae, and this allowed their classification into potential response pathways. CONCLUSIONS We have defined genes that are critical for amoebal survival during feeding on Gram(+), or Gram(-), bacteria that we propose form part of a regulatory network that allows D. discoideum to elicit specific cellular responses to different species of bacteria in order to optimize survival.

The Majority of 9,729 Group A Streptococcus Strains Causing Disease Secrete SpeB Cysteine Protease: Pathogenesis Implications

ABSTRACT 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.

Genomic Characteristics Behind the Spread of Bacteremic Group A Streptococcus Type emm89 in Finland, 2004–2014

Background. Many countries worldwide have reported increasing numbers of emm89 group A Streptococcus (GAS) infections during last decade. Pathogen genetic factors linked to this increase need assessment. Methods. We investigated epidemiological characteristics of emm89 GAS bacteremic infections, including 7-day and 30-day case-fatality rates, in Finland during 2004–2014 and linked them to whole-genome sequencing data obtained from corresponding strains. The Fisher exact test and exact logistic regression were used to compare differences between bacteremic infections due to emm89 GAS belonging to different genetic clades and subclades. Results. Out of 1928 cases of GAS bacteremic infection, 278 were caused by emm89 GAS. We identified 2 genetically distinct clades, arbitrarily designated clade 2 and clade 3. Both clades were present during 2004–2008, but clade 3 increased rapidly from 2009 onward. Six subclades (designated subclades A–F) were identified within clade 3, based on phylogenetic core genome analysis. The case-fatality rate differed significantly between subclades (P < .05), with subclade D having the highest 30-day estimated case-fatality rate (19% vs 3%–14%). Conclusions. A new emm89 clone, clade 3, emerged in 2009 and spread rapidly in Finland. Patients infected with certain subclades of clade 3 were significantly more likely to die. A specific polymerase chain reaction assay was developed to follow the spread of subclade D in 2015.

Genomic landscape of intrahost variation in Group A Streptococcus: repeated and abundant mutational inactivation of the fabT gene encoding a regulator of fatty acid synthesis.

ABSTRACT To obtain new information about Streptococcus pyogenes intrahost genetic variation during invasive infection, we sequenced the genomes of 2,954 serotype M1 strains recovered from a nonhuman primate experimental model of necrotizing fasciitis. A total of 644 strains (21.8%) acquired polymorphisms relative to the input parental strain. The fabT gene, encoding a transcriptional regulator of fatty acid biosynthesis genes, contained 54.5% of these changes. The great majority of polymorphisms were predicted to deleteriously alter FabT function. Transcriptome-sequencing (RNA-seq) analysis of a wild-type strain and an isogenic fabT deletion mutant strain found that between 3.7 and 28.5% of the S. pyogenes transcripts were differentially expressed, depending on the growth temperature (35°C or 40°C) and growth phase (mid-exponential or stationary phase). Genes implicated in fatty acid synthesis and lipid metabolism were significantly upregulated in the fabT deletion mutant strain. FabT also directly or indirectly regulated central carbon metabolism genes, including pyruvate hub enzymes and fermentation pathways and virulence genes. Deletion of fabT decreased virulence in a nonhuman primate model of necrotizing fasciitis. In addition, the fabT deletion strain had significantly decreased survival in human whole blood and during phagocytic interaction with polymorphonuclear leukocytes ex vivo. We conclude that FabT mutant progeny arise during infection, constitute a metabolically distinct subpopulation, and are less virulent in the experimental models used here.

Musser et al. Reply to "Emergence of the Same Successful Clade among Distinct Populations of emm89 Streptococcus pyogenes in Multiple Geographic Regions"

The gist of the letter by Friaes et al. (1) seems to be that an emerged clone of emm 89 Streptococcus pyogenes documented in three peer-reviewed publications (2–4) to be present in the United Kingdom, United States, Iceland, and Finland also may be present in Portugal. This conclusion was reached by analysis of publicly available full-genome sequence data for >900 emm 89 strains and comparisons performed with very limited genetic information (data from multilocus sequence type characterization, PCR analysis of the hasABC region, and DNA sequencing of the nga gene and its promoter region) for 95 Portuguese emm 89 strains. These Portuguese strains appear to be a convenience sample, rather than population based, an investigative approach that has led to erroneous conclusions, as described in previous studies (4). Unfortunately, the validity of the speculative claim by Friaes et al. (1) cannot be appropriately evaluated for at least three reasons. First, full-genome sequencing of the strains recovered from Portuguese patients was not conducted, which means that genomic relationships among the Portuguese and previously studied strains …

Whole exome sequencing study identifies novel rare and common Alzheimer’s-Associated variants involved in immune response and transcriptional regulation

The Alzheimer’s Disease Sequencing Project (ADSP) undertook whole exome sequencing in 5,740 late-onset Alzheimer disease (AD) cases and 5,096 cognitively normal controls primarily of European ancestry (EA), among whom 218 cases and 177 controls were Caribbean Hispanic (CH). An age-, sex- and APOE based risk score and family history were used to select cases most likely to harbor novel AD risk variants and controls least likely to develop AD by age 85 years. We tested ~1.5 million single nucleotide variants (SNVs) and 50,000 insertion-deletion polymorphisms (indels) for association to AD, using multiple models considering individual variants as well as gene-based tests aggregating rare, predicted functional, and loss of function variants. Sixteen single variants and 19 genes that met criteria for significant or suggestive associations after multiple-testing correction were evaluated for replication in four independent samples; three with whole exome sequencing (2,778 cases, 7,262 controls) and one with genome-wide genotyping imputed to the Haplotype Reference Consortium panel (9,343 cases, 11,527 controls). The top findings in the discovery sample were also followed-up in the ADSP whole-genome sequenced family-based dataset (197 members of 42 EA families and 501 members of 157 CH families). We identified novel and predicted functional genetic variants in genes previously associated with AD. We also detected associations in three novel genes: IGHG3 (p = 9.8 × 10−7), an immunoglobulin gene whose antibodies interact with β-amyloid, a long non-coding RNA AC099552.4 (p = 1.2 × 10−7), and a zinc-finger protein ZNF655 (gene-based p = 5.0 × 10−6). The latter two suggest an important role for transcriptional regulation in AD pathogenesis.

xAtlas: Scalable small variant calling across heterogeneous next-generation sequencing experiments

Motivation The rapid development of next-generation sequencing (NGS) technologies has lowered the barriers to genomic data generation, resulting in millions of samples sequenced across diverse experimental designs. The growing volume and heterogeneity of these sequencing data complicate the further optimization of methods for identifying DNA variation, especially considering that curated highconfidence variant call sets commonly used to evaluate these methods are generally developed by reference to results from the analysis of comparatively small and homogeneous sample sets. Results We have developed xAtlas, an application for the identification of single nucleotide variants (SNV) and small insertions and deletions (indels) in NGS data. xAtlas is easily scalable and enables execution and retraining with rapid development cycles. Generation of variant calls in VCF or gVCF format from BAM or CRAM alignments is accomplished in less than one CPU-hour per 30× short-read human whole-genome. The retraining capabilities of xAtlas allow its core variant evaluation models to be optimized on new sample data and user-defined truth sets. Obtaining SNV and indels calls from xAtlas can be achieved more than 40 times faster than established methods while retaining the same accuracy. Availability Freely available under a BSD 3-clause license at https://github.com/jfarek/xatlas. Contact farek@bcm.edu Supplementary information Supplementary data are available at Bioinformatics online.