1
|
Vieira A, Wan Y, Ryan Y, Li HK, Guy RL, Papangeli M, Huse KK, Reeves LC, Soo VWC, Daniel R, Harley A, Broughton K, Dhami C, Ganner M, Ganner MA, Mumin Z, Razaei M, Rundberg E, Mammadov R, Mills EA, Sgro V, Mok KY, Didelot X, Croucher NJ, Jauneikaite E, Lamagni T, Brown CS, Coelho J, Sriskandan S. Rapid expansion and international spread of M1 UK in the post-pandemic UK upsurge of Streptococcus pyogenes. Nat Commun 2024; 15:3916. [PMID: 38729927 PMCID: PMC11087535 DOI: 10.1038/s41467-024-47929-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2024] [Accepted: 04/15/2024] [Indexed: 05/12/2024] Open
Abstract
The UK observed a marked increase in scarlet fever and invasive group A streptococcal infection in 2022 with severe outcomes in children and similar trends worldwide. Here we report lineage M1UK to be the dominant source of invasive infections in this upsurge. Compared with ancestral M1global strains, invasive M1UK strains exhibit reduced genomic diversity and fewer mutations in two-component regulator genes covRS. The emergence of M1UK is dated to 2008. Following a bottleneck coinciding with the COVID-19 pandemic, three emergent M1UK clades underwent rapid nationwide expansion, despite lack of detection in previous years. All M1UK isolates thus-far sequenced globally have a phylogenetic origin in the UK, with dispersal of the new clades in Europe. While waning immunity may promote streptococcal epidemics, the genetic features of M1UK point to a fitness advantage in pathogenicity, and a striking ability to persist through population bottlenecks.
Collapse
Affiliation(s)
- Ana Vieira
- Department of Infectious Disease, Imperial College London, London, UK
- Centre for Bacterial Resistance Biology, Imperial College London, London, UK
- NIHR Health Protection Research Unit in Healthcare-associated Infections and AMR, Imperial College London, London, UK
| | - Yu Wan
- Department of Infectious Disease, Imperial College London, London, UK
- NIHR Health Protection Research Unit in Healthcare-associated Infections and AMR, Imperial College London, London, UK
- Healthcare-Associated Infections, Fungal, AMR, AMU, and Sepsis Division, UK Health Security Agency, London, UK
| | - Yan Ryan
- Reference Services Division, UK Health Security Agency, London, UK
| | - Ho Kwong Li
- Department of Infectious Disease, Imperial College London, London, UK
- Centre for Bacterial Resistance Biology, Imperial College London, London, UK
| | - Rebecca L Guy
- Healthcare-Associated Infections, Fungal, AMR, AMU, and Sepsis Division, UK Health Security Agency, London, UK
| | - Maria Papangeli
- Department of Infectious Disease, Imperial College London, London, UK
- Centre for Bacterial Resistance Biology, Imperial College London, London, UK
| | - Kristin K Huse
- Department of Infectious Disease, Imperial College London, London, UK
- Centre for Bacterial Resistance Biology, Imperial College London, London, UK
| | - Lucy C Reeves
- Department of Infectious Disease, Imperial College London, London, UK
- Centre for Bacterial Resistance Biology, Imperial College London, London, UK
| | - Valerie W C Soo
- Department of Infectious Disease, Imperial College London, London, UK
- Centre for Bacterial Resistance Biology, Imperial College London, London, UK
| | - Roger Daniel
- Reference Services Division, UK Health Security Agency, London, UK
| | | | - Karen Broughton
- Reference Services Division, UK Health Security Agency, London, UK
| | - Chenchal Dhami
- Reference Services Division, UK Health Security Agency, London, UK
| | - Mark Ganner
- Reference Services Division, UK Health Security Agency, London, UK
| | | | - Zaynab Mumin
- Reference Services Division, UK Health Security Agency, London, UK
| | - Maryam Razaei
- Reference Services Division, UK Health Security Agency, London, UK
| | - Emma Rundberg
- Reference Services Division, UK Health Security Agency, London, UK
| | - Rufat Mammadov
- Reference Services Division, UK Health Security Agency, London, UK
| | - Ewurabena A Mills
- Department of Infectious Disease, Imperial College London, London, UK
- Centre for Bacterial Resistance Biology, Imperial College London, London, UK
| | - Vincenzo Sgro
- Department of Infectious Disease, Imperial College London, London, UK
| | - Kai Yi Mok
- Department of Infectious Disease, Imperial College London, London, UK
| | - Xavier Didelot
- School of Life Sciences and Department of Statistics, University of Warwick, Coventry, UK
| | - Nicholas J Croucher
- School of Public Health, Imperial College London, London, UK
- MRC Centre for Global Infectious Disease Analysis, Imperial College London, London, UK
| | - Elita Jauneikaite
- NIHR Health Protection Research Unit in Healthcare-associated Infections and AMR, Imperial College London, London, UK
- School of Public Health, Imperial College London, London, UK
- MRC Centre for Global Infectious Disease Analysis, Imperial College London, London, UK
| | - Theresa Lamagni
- NIHR Health Protection Research Unit in Healthcare-associated Infections and AMR, Imperial College London, London, UK
- Healthcare-Associated Infections, Fungal, AMR, AMU, and Sepsis Division, UK Health Security Agency, London, UK
| | - Colin S Brown
- NIHR Health Protection Research Unit in Healthcare-associated Infections and AMR, Imperial College London, London, UK
- Healthcare-Associated Infections, Fungal, AMR, AMU, and Sepsis Division, UK Health Security Agency, London, UK
| | - Juliana Coelho
- NIHR Health Protection Research Unit in Healthcare-associated Infections and AMR, Imperial College London, London, UK.
- Healthcare-Associated Infections, Fungal, AMR, AMU, and Sepsis Division, UK Health Security Agency, London, UK.
- Reference Services Division, UK Health Security Agency, London, UK.
| | - Shiranee Sriskandan
- Department of Infectious Disease, Imperial College London, London, UK.
- Centre for Bacterial Resistance Biology, Imperial College London, London, UK.
- NIHR Health Protection Research Unit in Healthcare-associated Infections and AMR, Imperial College London, London, UK.
| |
Collapse
|
2
|
Lux J, Portmann H, Sánchez García L, Erhardt M, Holivololona L, Laloli L, Licheri MF, Gallay C, Hoepner R, Croucher NJ, Straume D, Veening JW, Dijkman R, Heller M, Grandgirard D, Leib SL, Hathaway LJ. Klebsiella pneumoniae peptide hijacks a Streptococcus pneumoniae permease to subvert pneumococcal growth and colonization. Commun Biol 2024; 7:425. [PMID: 38589539 PMCID: PMC11001997 DOI: 10.1038/s42003-024-06113-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/22/2023] [Accepted: 03/26/2024] [Indexed: 04/10/2024] Open
Abstract
Treatment of pneumococcal infections is limited by antibiotic resistance and exacerbation of disease by bacterial lysis releasing pneumolysin toxin and other inflammatory factors. We identified a previously uncharacterized peptide in the Klebsiella pneumoniae secretome, which enters Streptococcus pneumoniae via its AmiA-AliA/AliB permease. Subsequent downregulation of genes for amino acid biosynthesis and peptide uptake was associated with reduction of pneumococcal growth in defined medium and human cerebrospinal fluid, irregular cell shape, decreased chain length and decreased genetic transformation. The bacteriostatic effect was specific to S. pneumoniae and Streptococcus pseudopneumoniae with no effect on Streptococcus mitis, Haemophilus influenzae, Staphylococcus aureus or K. pneumoniae. Peptide sequence and length were crucial to growth suppression. The peptide reduced pneumococcal adherence to primary human airway epithelial cell cultures and colonization of rat nasopharynx, without toxicity. We identified a peptide with potential as a therapeutic for pneumococcal diseases suppressing growth of multiple clinical isolates, including antibiotic resistant strains, while avoiding bacterial lysis and dysbiosis.
Collapse
Affiliation(s)
- Janine Lux
- Faculty of Medicine, Institute for Infectious Diseases, University of Bern, Bern, Switzerland
- Graduate School for Cellular and Biomedical Sciences, University of Bern, Bern, Switzerland
| | - Hannah Portmann
- Faculty of Medicine, Institute for Infectious Diseases, University of Bern, Bern, Switzerland
| | - Lucía Sánchez García
- Faculty of Medicine, Institute for Infectious Diseases, University of Bern, Bern, Switzerland
| | - Maria Erhardt
- Faculty of Medicine, Institute for Infectious Diseases, University of Bern, Bern, Switzerland
| | - Lalaina Holivololona
- Faculty of Medicine, Institute for Infectious Diseases, University of Bern, Bern, Switzerland
| | - Laura Laloli
- Faculty of Medicine, Institute for Infectious Diseases, University of Bern, Bern, Switzerland
| | - Manon F Licheri
- Faculty of Medicine, Institute for Infectious Diseases, University of Bern, Bern, Switzerland
| | - Clement Gallay
- Department of Fundamental Microbiology, University of Lausanne, Lausanne, Switzerland
| | - Robert Hoepner
- Department of Neurology, Bern University Hospital and University of Bern, Bern, Switzerland
| | - Nicholas J Croucher
- MRC Centre for Global Infectious Disease Analysis, Sir Michael Uren Hub, White City Campus, Imperial College London, London, UK
| | - Daniel Straume
- Faculty of Chemistry, Biotechnology and Food Science, Norwegian University of Life Sciences, 1430, Ås, Norway
| | - Jan-Willem Veening
- Department of Fundamental Microbiology, University of Lausanne, Lausanne, Switzerland
| | - Ronald Dijkman
- Faculty of Medicine, Institute for Infectious Diseases, University of Bern, Bern, Switzerland
| | - Manfred Heller
- Proteomics and Mass Spectrometry Core Facility, Department for BioMedical Research (DBMR), University of Bern, Bern, Switzerland
| | - Denis Grandgirard
- Faculty of Medicine, Institute for Infectious Diseases, University of Bern, Bern, Switzerland
| | - Stephen L Leib
- Faculty of Medicine, Institute for Infectious Diseases, University of Bern, Bern, Switzerland
| | - Lucy J Hathaway
- Faculty of Medicine, Institute for Infectious Diseases, University of Bern, Bern, Switzerland.
| |
Collapse
|
3
|
Belman S, Pesonen H, Croucher NJ, Bentley SD, Corander J. Estimating Between Country Migration in Pneumococcal Populations. G3 (Bethesda) 2024:jkae058. [PMID: 38507601 DOI: 10.1093/g3journal/jkae058] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/22/2023] [Revised: 02/29/2024] [Accepted: 03/11/2024] [Indexed: 03/22/2024]
Abstract
Streptococcus pneumoniae (the pneumococcus) is a globally distributed, human obligate opportunistic bacterial pathogen which, although often carried commensally, is also a significant cause of invasive disease. Apart from multi-drug resistant and virulent clones, the rate and direction of pneumococcal dissemination between different countries remains largely unknown. The ability for the pneumococcus to take a foothold in a country depends on existing population configuration, the extent of vaccine implementation, as well as human mobility since it is a human obligate bacterium. To shed light on its international movement, we used extensive genome data from the Global Pneumococcal Sequencing (GPS) project and estimated migration parameters between multiple countries in Africa. Data on allele frequencies of polymorphisms at housekeeping-like loci for multiple different lineages circulating in the populations of South Africa, Malawi, Kenya, and The Gambia were used to calculate the fixation index (Fst) between countries. We then further used these summaries to fit migration coalescent models with the likelihood-free inference algorithms available in the ELFI software package. Synthetic data were additionally used to validate the inference approach. Our results demonstrate country-pair specific migration patterns and heterogeneity in the extent of migration between different lineages. Our approach demonstrates that coalescent models can be effectively used for inferring migration rates for bacterial species and lineages provided sufficiently granular population genomics surveillance data. Further it can demonstrate the connectivity of respiratory disease agents between countries to inform intervention policy in the longer term.
Collapse
Affiliation(s)
- Sophie Belman
- Parasites and Microbes, Wellcome Sanger Institute, Hinxton, Cambridgeshire, UK, CB10 1SA
| | - Henri Pesonen
- Oslo Centre for Biostatistics and Epidemiology, Oslo University Hospital, Oslo, Norway, 0372
| | - Nicholas J Croucher
- MRC Centre for Global Infectious Disease Analysis, Department of Infectious Disease Epidemiology, School of Public Health, White City Campus, Imperial College London, London, UK, W12 0BZ
| | - Stephen D Bentley
- Parasites and Microbes, Wellcome Sanger Institute, Hinxton, Cambridgeshire, UK, CB10 1SA
| | - Jukka Corander
- Parasites and Microbes, Wellcome Sanger Institute, Hinxton, Cambridgeshire, UK, CB10 1SA
- Department of Biostatistics, University of Oslo, Oslo, Norway, 0371
- Helsinki Institute for Information Technology HIIT, Department of Mathematics and Statistics, University of Helsinki, Helsinki, Finland, 02150
| |
Collapse
|
4
|
Pöntinen AK, Gladstone RA, Pesonen H, Pesonen M, Cléon F, Parcell BJ, Kallonen T, Simonsen GS, Croucher NJ, McNally A, Parkhill J, Johnsen PJ, Samuelsen Ø, Corander J. Modulation of multidrug-resistant clone success in Escherichia coli populations: a longitudinal, multi-country, genomic and antibiotic usage cohort study. Lancet Microbe 2024; 5:e142-e150. [PMID: 38219757 DOI: 10.1016/s2666-5247(23)00292-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/17/2023] [Revised: 08/24/2023] [Accepted: 09/01/2023] [Indexed: 01/16/2024]
Abstract
BACKGROUND The effect of antibiotic usage on the success of multidrug-resistant (MDR) clones in a population remains unclear. With this genomics-based molecular epidemiology study, we aimed to investigate the contribution of antibiotic use to Escherichia coli clone success, relative to intra-strain competition for colonisation and infection. METHODS We sequenced all the available E coli bloodstream infection isolates provided by the British Society for Antimicrobial Chemotherapy (BSAC) from 2012 to 2017 (n=718) and combined these with published data from the UK (2001-11; n=1090) and Norway (2002-17; n=3254). Defined daily dose (DDD) data from the European Centre for Disease Prevention and Control (retrieved on Sept 21, 2021) for major antibiotic classes (β-lactam, tetracycline, macrolide, sulfonamide, quinolone, and non-penicillin β-lactam) were used together with sequence typing, resistance profiling, regression analysis, and non-neutral Wright-Fisher simulation-based modelling to enable systematic comparison of resistance levels, clone success, and antibiotic usage between the UK and Norway. FINDINGS Sequence type (ST)73, ST131, ST95, and ST69 accounted for 892 (49·3%) of 1808 isolates in the BSAC collection. In the UK, the proportion of ST69 increased between 2001-10 and 2011-17 (p=0·0004), whereas the proportions of ST73 and ST95 did not vary between periods. ST131 expanded quickly after its emergence in 2003 and its prevalence remained consistent throughout the study period (apart from a brief decrease in 2009-10). The extended-spectrum β-lactamase (ESBL)-carrying, globally disseminated MDR clone ST131-C2 showed overall greater success in the UK (154 [56·8%] of 271 isolates in 2003-17) compared with Norway (51 [18·3%] of 278 isolates in 2002-17; p<0·0001). DDD data indicated higher total use of antimicrobials in the UK, driven mainly by the class of non-penicillin β-lactams, which were used between 2·7-times and 5·1-times more in the UK per annum (ratio mean 3·7 [SD 0·8]). This difference was associated with the higher success of the MDR clone ST131-C2 (pseudo-R2 69·1%). A non-neutral Wright-Fisher model replicated the observed expansion of non-MDR and MDR sequence types under higher DDD regimes. INTERPRETATION Our study indicates that resistance profiles of contemporaneously successful clones can vary substantially, warranting caution in the interpretation of correlations between aggregate measures of resistance and antibiotic usage. Our study further suggests that in countries with low-to-moderate use of antibiotics, such as the UK and Norway, the extent of non-penicillin β-lactam use modulates rather than determines the success of widely disseminated MDR ESBL-carrying E coli clones. Detailed understanding of underlying causal drivers of success is important for improved control of resistant pathogens. FUNDING Trond Mohn Foundation, Marie Skłodowska-Curie Actions, European Research Council, Royal Society, and Wellcome Trust.
Collapse
Affiliation(s)
- Anna K Pöntinen
- Department of Biostatistics, Faculty of Medicine, University of Oslo, Oslo, Norway; Norwegian National Advisory Unit on Detection of Antimicrobial Resistance, Department of Microbiology and Infection Control, University Hospital of North Norway, Tromsø, Norway.
| | - Rebecca A Gladstone
- Department of Biostatistics, Faculty of Medicine, University of Oslo, Oslo, Norway
| | - Henri Pesonen
- Department of Biostatistics, Faculty of Medicine, University of Oslo, Oslo, Norway
| | - Maiju Pesonen
- Department of Biostatistics, Faculty of Medicine, University of Oslo, Oslo, Norway; Oslo Centre for Biostatistics and Epidemiology, Oslo University Hospital Research Support Services, Oslo, Norway
| | - François Cléon
- Department of Pharmacy, Faculty of Health Sciences, UiT The Arctic University of Norway, Tromsø, Norway
| | | | - Teemu Kallonen
- Department of Clinical Microbiology, Turku University Hospital, Turku, Finland
| | - Gunnar Skov Simonsen
- Research Group for Host-Microbe Interaction, Faculty of Health Sciences, UiT The Arctic University of Norway, Tromsø, Norway; Department of Microbiology and Infection Control, University Hospital of North Norway, Tromsø, Norway
| | - Nicholas J Croucher
- MRC Centre for Global Infectious Disease Analysis, Department of Infectious Disease Epidemiology, Imperial College London, London, UK
| | - Alan McNally
- Institute of Microbiology and Infection, University of Birmingham, Birmingham, UK
| | - Julian Parkhill
- Department of Veterinary Medicine, University of Cambridge, Cambridge, UK
| | - Pål J Johnsen
- Department of Pharmacy, Faculty of Health Sciences, UiT The Arctic University of Norway, Tromsø, Norway
| | - Ørjan Samuelsen
- Norwegian National Advisory Unit on Detection of Antimicrobial Resistance, Department of Microbiology and Infection Control, University Hospital of North Norway, Tromsø, Norway; Department of Pharmacy, Faculty of Health Sciences, UiT The Arctic University of Norway, Tromsø, Norway
| | - Jukka Corander
- Department of Biostatistics, Faculty of Medicine, University of Oslo, Oslo, Norway; Parasites and Microbes, Wellcome Sanger Institute, Cambridge, UK; Helsinki Institute of Information Technology, Department of Mathematics and Statistics, University of Helsinki, Helsinki, Finland.
| |
Collapse
|
5
|
Croucher NJ, Campo JJ, Le TQ, Pablo JV, Hung C, Teng AA, Turner C, Nosten F, Bentley SD, Liang X, Turner P, Goldblatt D. Genomic and panproteomic analysis of the development of infant immune responses to antigenically-diverse pneumococci. Nat Commun 2024; 15:355. [PMID: 38191887 PMCID: PMC10774285 DOI: 10.1038/s41467-023-44584-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2023] [Accepted: 12/19/2023] [Indexed: 01/10/2024] Open
Abstract
Streptococcus pneumoniae (pneumococcus) is a nasopharyngeal commensal and respiratory pathogen. This study characterises the immunoglobulin G (IgG) repertoire recognising pneumococci from birth to 24 months old (mo) in a prospectively-sampled cohort of 63 children using a panproteome array. IgG levels are highest at birth, due to transplacental transmission of maternal antibodies. The subsequent emergence of responses to individual antigens exhibit distinct kinetics across the cohort. Stable differences in the strength of individuals' responses, correlating with maternal IgG concentrations, are established by 6 mo. By 12 mo, children develop unique antibody profiles that are boosted by re-exposure. However, some proteins only stimulate substantial responses in adults. Integrating genomic data on nasopharyngeal colonisation demonstrates rare pneumococcal antigens can elicit strong IgG levels post-exposure. Quantifying such responses to the diverse core loci (DCL) proteins is complicated by cross-immunity between variants. In particular, the conserved N terminus of DCL protein zinc metalloprotease B provokes the strongest early IgG responses. DCL proteins' ability to inhibit mucosal immunity likely explains continued pneumococcal carriage despite hosts' polyvalent antibody repertoire. Yet higher IgG levels are associated with reduced incidence, and severity, of pneumonia, demonstrating the importance of the heterogeneity in response strength and kinetics across antigens and individuals.
Collapse
Affiliation(s)
- Nicholas J Croucher
- MRC Centre for Global Infectious Disease Analysis, Department of Infectious Disease Epidemiology, School of Public Health, Imperial College London, London, W12 0BZ, UK.
| | - Joseph J Campo
- Antigen Discovery Inc, 1 Technology Drive, Irvine, CA, 92618, USA
| | - Timothy Q Le
- Antigen Discovery Inc, 1 Technology Drive, Irvine, CA, 92618, USA
| | - Jozelyn V Pablo
- Antigen Discovery Inc, 1 Technology Drive, Irvine, CA, 92618, USA
| | - Christopher Hung
- Antigen Discovery Inc, 1 Technology Drive, Irvine, CA, 92618, USA
| | - Andy A Teng
- Antigen Discovery Inc, 1 Technology Drive, Irvine, CA, 92618, USA
| | - Claudia Turner
- Cambodia Oxford Medical Research Unit, Angkor Hospital for Children, Siem Reap, 9V54+8FQ, Cambodia
- Centre for Tropical Medicine and Global Health, Nuffield Department of Medicine, University of Oxford, Oxford, OX3 7BN, UK
| | - François Nosten
- Centre for Tropical Medicine and Global Health, Nuffield Department of Medicine, University of Oxford, Oxford, OX3 7BN, UK
- Shoklo Malaria Research Unit, Mahidol-Oxford Tropical Medicine Research Unit, Faculty of Tropical Medicine, Mahidol University, Mae Sot, 63110, Thailand
| | - Stephen D Bentley
- Parasites & Microbes, Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton, Cambridge, CB10 1SA, UK
| | - Xiaowu Liang
- Antigen Discovery Inc, 1 Technology Drive, Irvine, CA, 92618, USA
| | - Paul Turner
- Cambodia Oxford Medical Research Unit, Angkor Hospital for Children, Siem Reap, 9V54+8FQ, Cambodia
- Centre for Tropical Medicine and Global Health, Nuffield Department of Medicine, University of Oxford, Oxford, OX3 7BN, UK
| | - David Goldblatt
- Great Ormond Street Institute of Child Health, University College London, London, WC1N 1EH, UK
| |
Collapse
|
6
|
Horsfield ST, Tonkin-Hill G, Croucher NJ, Lees JA. Accurate and fast graph-based pangenome annotation and clustering with ggCaller. Genome Res 2023; 33:1622-1637. [PMID: 37620118 PMCID: PMC10620059 DOI: 10.1101/gr.277733.123] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2023] [Accepted: 08/18/2023] [Indexed: 08/26/2023]
Abstract
Bacterial genomes differ in both gene content and sequence mutations, which underlie extensive phenotypic diversity, including variation in susceptibility to antimicrobials or vaccine-induced immunity. To identify and quantify important variants, all genes within a population must be predicted, functionally annotated, and clustered, representing the "pangenome." Despite the volume of genome data available, gene prediction and annotation are currently conducted in isolation on individual genomes, which is computationally inefficient and frequently inconsistent across genomes. Here, we introduce the open-source software graph-gene-caller (ggCaller). ggCaller combines gene prediction, functional annotation, and clustering into a single workflow using population-wide de Bruijn graphs, removing redundancy in gene annotation and resulting in more accurate gene predictions and orthologue clustering. We applied ggCaller to simulated and real-world bacterial data sets containing hundreds or thousands of genomes, comparing it to current state-of-the-art tools. ggCaller has considerable speed-ups with equivalent or greater accuracy, particularly with data sets containing complex sources of error, such as assembly contamination or fragmentation. ggCaller is also an important extension to bacterial genome-wide association studies, enabling querying of annotated graphs for functional analyses. We highlight this application by functionally annotating DNA sequences with significant associations to tetracycline and macrolide resistance in Streptococcus pneumoniae, identifying key resistance determinants that were missed when using only a single reference genome. ggCaller is a novel bacterial genome analysis tool with applications in bacterial evolution and epidemiology.
Collapse
Affiliation(s)
- Samuel T Horsfield
- MRC Centre for Global Infectious Disease Analysis, Department of Infectious Disease Epidemiology, Imperial College London, London W12 0BZ, United Kingdom;
- European Molecular Biology Laboratory, European Bioinformatics Institute, Wellcome Genome Campus, Hinxton CB10 1SD, United Kingdom
| | - Gerry Tonkin-Hill
- Department of Biostatistics, University of Oslo, Blindern, 0372 Oslo, Norway
| | - Nicholas J Croucher
- MRC Centre for Global Infectious Disease Analysis, Department of Infectious Disease Epidemiology, Imperial College London, London W12 0BZ, United Kingdom
| | - John A Lees
- MRC Centre for Global Infectious Disease Analysis, Department of Infectious Disease Epidemiology, Imperial College London, London W12 0BZ, United Kingdom
- European Molecular Biology Laboratory, European Bioinformatics Institute, Wellcome Genome Campus, Hinxton CB10 1SD, United Kingdom
| |
Collapse
|
7
|
Glanville DG, Gazioglu O, Marra M, Tokars VL, Kushnir T, Habtom M, Croucher NJ, Nebenzahl YM, Mondragón A, Yesilkaya H, Ulijasz AT. Pneumococcal capsule expression is controlled through a conserved, distal cis-regulatory element during infection. PLoS Pathog 2023; 19:e1011035. [PMID: 36719895 PMCID: PMC9888711 DOI: 10.1371/journal.ppat.1011035] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2022] [Accepted: 11/29/2022] [Indexed: 02/01/2023] Open
Abstract
Streptococcus pneumoniae (the pneumococcus) is the major cause of bacterial pneumonia in the US and worldwide. Studies have shown that the differing chemical make-up between serotypes of its most important virulence factor, the capsule, can dictate disease severity. Here we demonstrate that control of capsule synthesis is also critical for infection and facilitated by two broadly conserved transcription factors, SpxR and CpsR, through a distal cis-regulatory element we name the 37-CE. Strikingly, changing only three nucleotides within this sequence is sufficient to render pneumococcus avirulent. Using in vivo and in vitro approaches, we present a model where SpxR interacts as a unique trimeric quaternary structure with the 37-CE to enable capsule repression in the airways. Considering its dramatic effect on infection, variation of the 37-CE between serotypes suggests this molecular switch could be a critical contributing factor to this pathogen's serotype-specific disease outcomes.
Collapse
Affiliation(s)
- David G. Glanville
- Department of Microbiology and Immunology, Loyola University Chicago, Maywood, Illinois, United States of America
| | - Ozcan Gazioglu
- Department of Respiratory Sciences, University of Leicester, University Road, Leicester, United Kingdom
| | - Michela Marra
- Department of Microbiology and Immunology, Loyola University Chicago, Maywood, Illinois, United States of America
| | - Valerie L. Tokars
- Department of Pharmacology, Feinberg School of Medicine, Northwestern University, Chicago, Illinois, United States of America
| | - Tatyana Kushnir
- The Shraga Segal Department of Microbiology, Immunology and Genetics, Faculty of Health Sciences, Ben-Gurion University of The Negev, Beer-Sheva, Israel
| | - Medhanie Habtom
- Department of Respiratory Sciences, University of Leicester, University Road, Leicester, United Kingdom
| | - Nicholas J. Croucher
- MRC Centre for Global Infectious Disease Analysis, Department of Infectious Disease Epidemiology, Sir Michael Uren Hub, Imperial College London, London, United Kingdom
| | - Yaffa Mizrachi Nebenzahl
- The Shraga Segal Department of Microbiology, Immunology and Genetics, Faculty of Health Sciences, Ben-Gurion University of The Negev, Beer-Sheva, Israel
| | - Alfonso Mondragón
- Department of Molecular Biosciences, Northwestern University, Evanston, Illinois, United States of America
| | - Hasan Yesilkaya
- Department of Respiratory Sciences, University of Leicester, University Road, Leicester, United Kingdom
| | - Andrew T. Ulijasz
- Department of Microbiology and Immunology, Loyola University Chicago, Maywood, Illinois, United States of America
| |
Collapse
|
8
|
Kwun MJ, Ion AV, Cheng HC, D’Aeth JC, Dougan S, Oggioni MR, Goulding DA, Bentley SD, Croucher NJ. Post-vaccine epidemiology of serotype 3 pneumococci identifies transformation inhibition through prophage-driven alteration of a non-coding RNA. Genome Med 2022; 14:144. [PMID: 36539881 PMCID: PMC9764711 DOI: 10.1186/s13073-022-01147-2] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/09/2022] [Accepted: 11/29/2022] [Indexed: 12/24/2022] Open
Abstract
BACKGROUND The respiratory pathogen Streptococcus pneumoniae (the pneumococcus) is a genetically diverse bacterium associated with over 101 immunologically distinct polysaccharide capsules (serotypes). Polysaccharide conjugate vaccines (PCVs) have successfully eliminated multiple targeted serotypes, yet the mucoid serotype 3 has persisted despite its inclusion in PCV13. This capsule type is predominantly associated with a single globally disseminated strain, GPSC12 (clonal complex 180). METHODS A genomic epidemiology study combined previous surveillance datasets of serotype 3 pneumococci to analyse the population structure, dynamics, and differences in rates of diversification within GPSC12 during the period of PCV introductions. Transcriptomic analyses, whole genome sequencing, mutagenesis, and electron microscopy were used to characterise the phenotypic impact of loci hypothesised to affect this strain's evolution. RESULTS GPSC12 was split into clades by a genomic analysis. Clade I, the most common, rarely underwent transformation, but was typically infected with the prophage ϕOXC141. Prior to the introduction of PCV13, this clade's composition shifted towards a ϕOXC141-negative subpopulation in a systematically sampled UK collection. In the post-PCV13 era, more rapidly recombining non-Clade I isolates, also ϕOXC141-negative, have risen in prevalence. The low in vitro transformation efficiency of a Clade I isolate could not be fully explained by the ~100-fold reduction attributable to the serotype 3 capsule. Accordingly, prophage ϕOXC141 was found to modify csRNA3, a non-coding RNA that inhibits the induction of transformation. This alteration was identified in ~30% of all pneumococci and was particularly common in the unusually clonal serotype 1 GPSC2 strain. RNA-seq and quantitative reverse transcriptase PCR experiments using a genetically tractable pneumococcus demonstrated the altered csRNA3 was more effective at inhibiting production of the competence-stimulating peptide pheromone. This resulted in a reduction in the induction of competence for transformation. CONCLUSION This interference with the quorum sensing needed to induce competence reduces the risk of the prophage being deleted by homologous recombination. Hence the selfish prophage-driven alteration of a regulatory RNA limits cell-cell communication and horizontal gene transfer, complicating the interpretation of post-vaccine population dynamics.
Collapse
Affiliation(s)
- Min Jung Kwun
- grid.7445.20000 0001 2113 8111MRC Centre for Global Infectious Disease Analysis, Department of Infectious Disease Epidemiology, School of Public Health, White City Campus, Imperial College London, London, W12 0BZ UK
| | - Alexandru V. Ion
- grid.7445.20000 0001 2113 8111MRC Centre for Global Infectious Disease Analysis, Department of Infectious Disease Epidemiology, School of Public Health, White City Campus, Imperial College London, London, W12 0BZ UK
| | - Hsueh-Chien Cheng
- grid.10306.340000 0004 0606 5382Parasites & Microbes, Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton, Cambridge, CB10 1SA UK
| | - Joshua C. D’Aeth
- grid.7445.20000 0001 2113 8111MRC Centre for Global Infectious Disease Analysis, Department of Infectious Disease Epidemiology, School of Public Health, White City Campus, Imperial College London, London, W12 0BZ UK
| | - Sam Dougan
- grid.10306.340000 0004 0606 5382Parasites & Microbes, Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton, Cambridge, CB10 1SA UK
| | - Marco R. Oggioni
- grid.9918.90000 0004 1936 8411Department of Genetics, University of Leicester, University Road, Leicester, LE1 7RH UK ,grid.6292.f0000 0004 1757 1758Dipartimento di Farmacia e Biotecnologie, Università di Bologna, Via Irnerio 42, 40126 Bologna, Italy
| | - David A. Goulding
- grid.10306.340000 0004 0606 5382Parasites & Microbes, Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton, Cambridge, CB10 1SA UK
| | - Stephen D. Bentley
- grid.10306.340000 0004 0606 5382Parasites & Microbes, Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton, Cambridge, CB10 1SA UK
| | - Nicholas J. Croucher
- grid.7445.20000 0001 2113 8111MRC Centre for Global Infectious Disease Analysis, Department of Infectious Disease Epidemiology, School of Public Health, White City Campus, Imperial College London, London, W12 0BZ UK
| |
Collapse
|
9
|
Lehtinen S, Croucher NJ, Blanquart F, Fraser C. Epidemiological dynamics of bacteriocin competition and antibiotic resistance. Proc Biol Sci 2022; 289:20221197. [PMID: 36196547 PMCID: PMC9532987 DOI: 10.1098/rspb.2022.1197] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
Bacteriocins, toxic peptides involved in the competition between bacterial strains, are extremely diverse. Previous work on bacteriocin dynamics has highlighted the role of non-transitive 'rock-paper-scissors' competition in maintaining the coexistence of different bacteriocin profiles. The focus to date has primarily been on bacteriocin interactions at the within-host scale (i.e. within a single bacterial population). Yet in species such as Streptococcus pneumoniae, with relatively short periods of colonization and limited within-host diversity, ecological outcomes are also shaped by processes at the epidemiological (between-host) scale. Here, we first investigate bacteriocin dynamics and diversity in epidemiological models. We find that in these models, bacteriocin diversity is more readily maintained than in within-host models, and with more possible combinations of coexisting bacteriocin profiles. Indeed, maintenance of diversity in epidemiological models does not require rock-paper-scissors dynamics; it can also occur through a competition-colonization trade-off. Second, we investigate the link between bacteriocin diversity and diversity at antibiotic resistance loci. Previous work has proposed that bacterial duration of colonization modulates the fitness of antibiotic resistance. Due to their inhibitory effects, bacteriocins are a plausible candidate for playing a role in the duration of colonization episodes. We extend the epidemiological model of bacteriocin dynamics to incorporate an antibiotic resistance locus and demonstrate that bacteriocin diversity can indeed maintain the coexistence of antibiotic-sensitive and -resistant strains.
Collapse
Affiliation(s)
- Sonja Lehtinen
- Department of Environmental System Science, Institute for Integrative Biology, ETH Zürich, Zürich, Switzerland
| | - Nicholas J Croucher
- MRC Centre for Global Infectious Disease Epidemiology, Imperial College London, London, UK
| | - François Blanquart
- Centre for Interdisciplinary Research in Biology (CIRB), Collège de France, CNRS, INSERM, PSL Research University, Paris, France.,Infection Antimicrobials Modelling Evolution, UMR, 1137, INSERM, Université de Paris, Paris, France
| | - Christophe Fraser
- Big Data Institute, Nuffield Department of Medicine, University of Oxford, Oxford, UK
| |
Collapse
|
10
|
Chaguza C, Ebruke C, Senghore M, Lo SW, Tientcheu PE, Gladstone RA, Tonkin-Hill G, Cornick JE, Yang M, Worwui A, McGee L, Breiman RF, Klugman KP, Kadioglu A, Everett DB, Mackenzie G, Croucher NJ, Roca A, Kwambana-Adams BA, Antonio M, Bentley SD. Comparative Genomics of Disease and Carriage Serotype 1 Pneumococci. Genome Biol Evol 2022; 14:evac052. [PMID: 35439297 PMCID: PMC9048925 DOI: 10.1093/gbe/evac052] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 04/12/2022] [Indexed: 11/14/2022] Open
Abstract
The isolation of Streptococcus pneumoniae serotypes in systemic tissues of patients with invasive disease versus the nasopharynx of healthy individuals with asymptomatic carriage varies widely. Some serotypes are hyper-invasive, particularly serotype 1, but the underlying genetics remain poorly understood due to the rarity of carriage isolates, reducing the power of comparison with invasive isolates. Here, we use a well-controlled genome-wide association study to search for genetic variation associated with invasiveness of serotype 1 pneumococci from a serotype 1 endemic setting in Africa. We found no consensus evidence that certain genomic variation is overrepresented among isolates from patients with invasive disease than asymptomatic carriage. Overall, the genomic variation explained negligible phenotypic variability, suggesting a minimal effect on the disease status. Furthermore, changes in lineage distribution were seen with lineages replacing each other over time, highlighting the importance of continued pathogen surveillance. Our findings suggest that the hyper-invasiveness is an intrinsic property of the serotype 1 strains, not specific for a "disease-associated" subpopulation disproportionately harboring unique genomic variation.
Collapse
Affiliation(s)
- Chrispin Chaguza
- Parasites and Microbes Programme, Wellcome Sanger Institute, Wellcome Genome Campus, Cambridge, UK
- Darwin College, University of Cambridge, Silver Street, Cambridge, UK
- Department of Clinical Infection, Microbiology and Immunology, Institute of Infection, Veterinary and Ecological Sciences, University of Liverpool, Liverpool, UK
| | - Chinelo Ebruke
- Medical Research Council (MRC) Unit The Gambia at the London School of Hygiene and Tropical Medicine, Fajara, The Gambia
| | - Madikay Senghore
- Medical Research Council (MRC) Unit The Gambia at the London School of Hygiene and Tropical Medicine, Fajara, The Gambia
- Department of Epidemiology, Center for Communicable Disease Dynamics, Harvard T.H. Chan School of Public Health, Boston, MA, USA
| | - Stephanie W. Lo
- Parasites and Microbes Programme, Wellcome Sanger Institute, Wellcome Genome Campus, Cambridge, UK
| | - Peggy-Estelle Tientcheu
- Medical Research Council (MRC) Unit The Gambia at the London School of Hygiene and Tropical Medicine, Fajara, The Gambia
| | - Rebecca A. Gladstone
- Parasites and Microbes Programme, Wellcome Sanger Institute, Wellcome Genome Campus, Cambridge, UK
- Department of Biostatistics, University of Oslo, Oslo, Norway
| | - Gerry Tonkin-Hill
- Parasites and Microbes Programme, Wellcome Sanger Institute, Wellcome Genome Campus, Cambridge, UK
| | - Jennifer E. Cornick
- Department of Clinical Infection, Microbiology and Immunology, Institute of Infection, Veterinary and Ecological Sciences, University of Liverpool, Liverpool, UK
- Malawi-Liverpool-Wellcome Trust Clinical Research Programme, Blantyre, Malawi
| | - Marie Yang
- Department of Clinical Infection, Microbiology and Immunology, Institute of Infection, Veterinary and Ecological Sciences, University of Liverpool, Liverpool, UK
| | - Archibald Worwui
- Medical Research Council (MRC) Unit The Gambia at the London School of Hygiene and Tropical Medicine, Fajara, The Gambia
| | - Lesley McGee
- Respiratory Diseases Branch, Centers for Disease Control and Prevention, Atlanta, GA, USA
| | - Robert F. Breiman
- Hubert Department of Global Health, Rollins School of Public Health, Emory University, Atlanta, GA, USA
| | - Keith P. Klugman
- Hubert Department of Global Health, Rollins School of Public Health, Emory University, Atlanta, GA, USA
| | - Aras Kadioglu
- Department of Clinical Infection, Microbiology and Immunology, Institute of Infection, Veterinary and Ecological Sciences, University of Liverpool, Liverpool, UK
| | - Dean B. Everett
- College of Medicine and Health Sciences, Khalifa University, Abu Dhabi, UAE
| | - Grant Mackenzie
- Medical Research Council (MRC) Unit The Gambia at the London School of Hygiene and Tropical Medicine, Fajara, The Gambia
- Murdoch Children’s Research Institute, Parkville, Melbourne, VIC, Australia
- London School of Hygiene & Tropical Medicine, London, UK
| | - Nicholas J. Croucher
- MRC Centre for Global Infectious Disease Analysis, Department of Infectious Disease Epidemiology, School of Public Health, Imperial College London, London, UK
| | - Anna Roca
- Medical Research Council (MRC) Unit The Gambia at the London School of Hygiene and Tropical Medicine, Fajara, The Gambia
- London School of Hygiene & Tropical Medicine, London, UK
| | - Brenda A. Kwambana-Adams
- Medical Research Council (MRC) Unit The Gambia at the London School of Hygiene and Tropical Medicine, Fajara, The Gambia
- NIHR Global Health Research Unit on Mucosal Pathogens, Division of Infection and Immunity, University College London, London, UK
| | - Martin Antonio
- Medical Research Council (MRC) Unit The Gambia at the London School of Hygiene and Tropical Medicine, Fajara, The Gambia
- London School of Hygiene & Tropical Medicine, London, UK
- Warwick Medical School, University of Warwick, Coventry, UK
| | - Stephen D. Bentley
- Parasites and Microbes Programme, Wellcome Sanger Institute, Wellcome Genome Campus, Cambridge, UK
| |
Collapse
|
11
|
Stevens EJ, Morse DJ, Bonini D, Duggan S, Brignoli T, Recker M, Lees JA, Croucher NJ, Bentley S, Wilson DJ, Earle SG, Dixon R, Nobbs A, Jenkinson H, van Opijnen T, Thibault D, Wilkinson OJ, Dillingham MS, Carlile S, McLoughlin RM, Massey RC. Targeted control of pneumolysin production by a mobile genetic element in Streptococcus pneumoniae. Microb Genom 2022; 8:000784. [PMID: 35416147 PMCID: PMC9453066 DOI: 10.1099/mgen.0.000784] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Streptococcus pneumoniae is a major human pathogen that can cause severe invasive diseases such as pneumonia, septicaemia and meningitis. Young children are at a particularly high risk, with an estimated 3-4 million cases of severe disease and between 300 000 and 500 000 deaths attributable to pneumococcal disease each year. The haemolytic toxin pneumolysin (Ply) is a primary virulence factor for this bacterium, yet despite its key role in pathogenesis, immune evasion and transmission, the regulation of Ply production is not well defined. Using a genome-wide association approach, we identified a large number of potential affectors of Ply activity, including a gene acquired horizontally on the antibiotic resistance-conferring Integrative and Conjugative Element (ICE) ICESp23FST81. This gene encodes a novel modular protein, ZomB, which has an N-terminal UvrD-like helicase domain followed by two Cas4-like domains with potent ATP-dependent nuclease activity. We found the regulatory effect of ZomB to be specific for the ply operon, potentially mediated by its high affinity for the BOX repeats encoded therein. Using a murine model of pneumococcal colonization, we further demonstrate that a ZomB mutant strain colonizes both the upper respiratory tract and lungs at higher levels when compared to the wild-type strain. While the antibiotic resistance-conferring aspects of ICESp23FST81 are often credited with contributing to the success of the S. pneumoniae lineages that acquire it, its ability to control the expression of a major virulence factor implicated in bacterial transmission is also likely to have played an important role.
Collapse
Affiliation(s)
- Emily J Stevens
- School of Cellular and Molecular Medicine, University of Bristol, Bristol, BS8 1TD, UK
| | - Daniel J Morse
- School of Cellular and Molecular Medicine, University of Bristol, Bristol, BS8 1TD, UK
| | - Dora Bonini
- School of Cellular and Molecular Medicine, University of Bristol, Bristol, BS8 1TD, UK
| | - Seána Duggan
- School of Cellular and Molecular Medicine, University of Bristol, Bristol, BS8 1TD, UK
| | - Tarcisio Brignoli
- School of Cellular and Molecular Medicine, University of Bristol, Bristol, BS8 1TD, UK
| | - Mario Recker
- Centre for Ecology and Conservation, University of Exeter, Penryn Campus, Exeter, TR10 9FE, UK.,Institute of Tropical Medicine, University of Tübingen, Tübingen, Germany
| | - John A Lees
- MRC Centre for Global Infectious Disease Analysis, Department of Infectious Disease Epidemiology, St. Mary's Campus, Imperial College London, London, W2 1PG, UK
| | - Nicholas J Croucher
- MRC Centre for Global Infectious Disease Analysis, Department of Infectious Disease Epidemiology, St. Mary's Campus, Imperial College London, London, W2 1PG, UK
| | - Stephen Bentley
- The Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge, CB10 1SA, UK
| | - Daniel J Wilson
- Big Data Institute, Nuffield Department of Population Health, University of Oxford, Oxford, OX3 7LF, UK
| | - Sarah G Earle
- Big Data Institute, Nuffield Department of Population Health, University of Oxford, Oxford, OX3 7LF, UK
| | - Robert Dixon
- Big Data Institute, Nuffield Department of Population Health, University of Oxford, Oxford, OX3 7LF, UK
| | - Angela Nobbs
- Bristol Dental School, University of Bristol, Bristol, BS1 2LY, UK
| | - Howard Jenkinson
- Bristol Dental School, University of Bristol, Bristol, BS1 2LY, UK
| | | | - Derek Thibault
- Biology Department, Boston College, Chestnut Hill, MA, USA
| | - Oliver J Wilkinson
- DNA-Protein Interactions Unit, School of Biochemistry, University of Bristol, Bristol, BS8 1TD, UK
| | - Mark S Dillingham
- DNA-Protein Interactions Unit, School of Biochemistry, University of Bristol, Bristol, BS8 1TD, UK
| | - Simon Carlile
- Host Pathogen Interactions Group, School of Biochemistry and Immunology, Trinity Biomedical Sciences Institute, Trinity College Dublin, Dublin 2, Ireland
| | - Rachel M McLoughlin
- Host Pathogen Interactions Group, School of Biochemistry and Immunology, Trinity Biomedical Sciences Institute, Trinity College Dublin, Dublin 2, Ireland
| | - Ruth C Massey
- School of Cellular and Molecular Medicine, University of Bristol, Bristol, BS8 1TD, UK.,Schools of Microbiology and Medicine and APC Microbiome Ireland, University College Cork, Cork, Ireland
| |
Collapse
|
12
|
Mallawaarachchi S, Tonkin-Hill G, Croucher NJ, Turner P, Speed D, Corander J, Balding D. Genome-wide association, prediction and heritability in bacteria with application to Streptococcus pneumoniae. NAR Genom Bioinform 2022; 4:lqac011. [PMID: 35211669 PMCID: PMC8862724 DOI: 10.1093/nargab/lqac011] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2021] [Revised: 01/06/2022] [Accepted: 02/01/2022] [Indexed: 11/14/2022] Open
Abstract
Whole-genome sequencing has facilitated genome-wide analyses of association, prediction and heritability in many organisms. However, such analyses in bacteria are still in their infancy, being limited by difficulties including genome plasticity and strong population structure. Here we propose a suite of methods including linear mixed models, elastic net and LD-score regression, adapted to bacterial traits using innovations such as frequency-based allele coding, both insertion/deletion and nucleotide testing and heritability partitioning. We compare and validate our methods against the current state-of-art using simulations, and analyse three phenotypes of the major human pathogen Streptococcus pneumoniae, including the first analyses of minimum inhibitory concentrations (MIC) for penicillin and ceftriaxone. We show that the MIC traits are highly heritable with high prediction accuracy, explained by many genetic associations under good population structure control. In ceftriaxone MIC, this is surprising because none of the isolates are resistant as per the inhibition zone criteria. We estimate that half of the heritability of penicillin MIC is explained by a known drug-resistance region, which also contributes a quarter of the ceftriaxone MIC heritability. For the within-host carriage duration phenotype, no associations were observed, but the moderate heritability and prediction accuracy indicate a moderately polygenic trait.
Collapse
Affiliation(s)
| | - Gerry Tonkin-Hill
- Parasites and Microbes, Wellcome Sanger Institute, Cambridge CB10 1SA, UK
| | - Nicholas J Croucher
- Faculty of Medicine, School of Public Health, Imperial College, London SW7 2AZ, UK
| | - Paul Turner
- Cambodia-Oxford Medical Research Unit, Angkor Hospital for Children, Siem Reap 1710, Cambodia,Centre for Tropical Medicine and Global Health, Nuffield Department of Medicine, University of Oxford, Oxford OX3 7LG, UK
| | - Doug Speed
- Aarhus Institute of Advanced Studies (AIAS), Aarhus University, 8000 Aarhus, Denmark,Bioinformatics Research Centre, Aarhus University, 8000 Aarhus, Denmark,UCL Genetics Institute, University College London, London WC1E 6BT, United Kingdom
| | - Jukka Corander
- Parasites and Microbes, Wellcome Sanger Institute, Cambridge CB10 1SA, UK,Department of Biostatistics, Faculty of Medicine, University of Oslo, 0372 Oslo, Norway,Helsinki Institute of Information Technology, Department of Mathematics and Statistics, University of Helsinki, Helsinki 00014, Finland
| | - David Balding
- Correspondence may also be addressed to David Balding.
| |
Collapse
|
13
|
Tonkin-Hill G, Ling C, Chaguza C, Salter SJ, Hinfonthong P, Nikolaou E, Tate N, Pastusiak A, Turner C, Chewapreecha C, Frost SDW, Corander J, Croucher NJ, Turner P, Bentley SD. Pneumococcal within-host diversity during colonization, transmission and treatment. Nat Microbiol 2022; 7:1791-1804. [PMID: 36216891 PMCID: PMC9613479 DOI: 10.1038/s41564-022-01238-1] [Citation(s) in RCA: 17] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2022] [Accepted: 07/18/2022] [Indexed: 11/13/2022]
Abstract
Characterizing the genetic diversity of pathogens within the host promises to greatly improve surveillance and reconstruction of transmission chains. For bacteria, it also informs our understanding of inter-strain competition and how this shapes the distribution of resistant and sensitive bacteria. Here we study the genetic diversity of Streptococcus pneumoniae within 468 infants and 145 of their mothers by deep sequencing whole pneumococcal populations from 3,761 longitudinal nasopharyngeal samples. We demonstrate that deep sequencing has unsurpassed sensitivity for detecting multiple colonization, doubling the rate at which highly invasive serotype 1 bacteria were detected in carriage compared with gold-standard methods. The greater resolution identified an elevated rate of transmission from mothers to their children in the first year of the child's life. Comprehensive treatment data demonstrated that infants were at an elevated risk of both the acquisition and persistent colonization of a multidrug-resistant bacterium following antimicrobial treatment. Some alleles were enriched after antimicrobial treatment, suggesting that they aided persistence, but generally purifying selection dominated within-host evolution. Rates of co-colonization imply that in the absence of treatment, susceptible lineages outcompeted resistant lineages within the host. These results demonstrate the many benefits of deep sequencing for the genomic surveillance of bacterial pathogens.
Collapse
Affiliation(s)
- Gerry Tonkin-Hill
- grid.10306.340000 0004 0606 5382Parasites and Microbes, Wellcome Sanger Institute, Cambridge, UK ,grid.5510.10000 0004 1936 8921Department of Biostatistics, University of Oslo, Blindern, Norway
| | - Clare Ling
- grid.10223.320000 0004 1937 0490Shoklo Malaria Research Unit, Mahidol-Oxford Tropical Medicine Research Unit, Faculty of Tropical Medicine, Mahidol University, Mae Sot, Thailand ,grid.4991.50000 0004 1936 8948Centre for Tropical Medicine and Global Health, Nuffield Department of Medicine, University of Oxford, Oxford, UK
| | - Chrispin Chaguza
- grid.10306.340000 0004 0606 5382Parasites and Microbes, Wellcome Sanger Institute, Cambridge, UK ,grid.47100.320000000419368710Department of Epidemiology of Microbial Diseases, Yale School of Public Health, Yale University, New Haven, CT USA
| | - Susannah J. Salter
- grid.5335.00000000121885934Department of Veterinary Medicine, University of Cambridge, Cambridge, UK
| | - Pattaraporn Hinfonthong
- grid.10223.320000 0004 1937 0490Shoklo Malaria Research Unit, Mahidol-Oxford Tropical Medicine Research Unit, Faculty of Tropical Medicine, Mahidol University, Mae Sot, Thailand
| | - Elissavet Nikolaou
- grid.48004.380000 0004 1936 9764Department of Clinical Sciences, Liverpool School of Tropical Medicine, Liverpool, UK ,grid.1058.c0000 0000 9442 535XInfection and Immunity, Murdoch Children’s Research Institute, Melbourne, Victoria Australia ,grid.1008.90000 0001 2179 088XDepartment of Microbiology and Immunology, Peter Doherty Institute for Infection and Immunity, University of Melbourne, Melbourne, Victoria Australia
| | - Natalie Tate
- grid.48004.380000 0004 1936 9764Department of Clinical Sciences, Liverpool School of Tropical Medicine, Liverpool, UK
| | | | - Claudia Turner
- grid.4991.50000 0004 1936 8948Centre for Tropical Medicine and Global Health, Nuffield Department of Medicine, University of Oxford, Oxford, UK ,grid.459332.a0000 0004 0418 5364Cambodia-Oxford Medical Research Unit, Angkor Hospital for Children, Siem Reap, Cambodia
| | - Claire Chewapreecha
- grid.10306.340000 0004 0606 5382Parasites and Microbes, Wellcome Sanger Institute, Cambridge, UK ,grid.10223.320000 0004 1937 0490Mahidol-Oxford Tropical Medicine Research Unit, Faculty of Tropical Medicine, Mahidol University, Bangkok, Thailand
| | - Simon D. W. Frost
- grid.419815.00000 0001 2181 3404Microsoft Research, Redmond, WA USA ,grid.8991.90000 0004 0425 469XLondon School of Hygiene and Tropical Medicine, London, UK
| | - Jukka Corander
- grid.10306.340000 0004 0606 5382Parasites and Microbes, Wellcome Sanger Institute, Cambridge, UK ,grid.5510.10000 0004 1936 8921Department of Biostatistics, University of Oslo, Blindern, Norway ,grid.7737.40000 0004 0410 2071Helsinki Institute for Information Technology HIIT, Department of Mathematics and Statistics, University of Helsinki, Helsinki, Finland
| | - Nicholas J. Croucher
- grid.7445.20000 0001 2113 8111MRC Centre for Global Infectious Disease Analysis, Department of Infectious Disease Epidemiology, Imperial College London, London, UK
| | - Paul Turner
- grid.4991.50000 0004 1936 8948Centre for Tropical Medicine and Global Health, Nuffield Department of Medicine, University of Oxford, Oxford, UK ,grid.459332.a0000 0004 0418 5364Cambodia-Oxford Medical Research Unit, Angkor Hospital for Children, Siem Reap, Cambodia
| | - Stephen D. Bentley
- grid.10306.340000 0004 0606 5382Parasites and Microbes, Wellcome Sanger Institute, Cambridge, UK
| |
Collapse
|
14
|
Wan Y, Mills E, Leung RC, Vieira A, Zhi X, Croucher NJ, Woodford N, Jauneikaite E, Ellington MJ, Sriskandan S. Alterations in chromosomal genes nfsA, nfsB, and ribE are associated with nitrofurantoin resistance in Escherichia coli from the United Kingdom. Microb Genom 2021; 7:000702. [PMID: 34860151 PMCID: PMC8767348 DOI: 10.1099/mgen.0.000702] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023] Open
Abstract
Antimicrobial resistance in enteric or urinary Escherichia coli is a risk factor for invasive E. coli infections. Due to widespread trimethoprim resistance amongst urinary E. coli and increased bacteraemia incidence, a national recommendation to prescribe nitrofurantoin for uncomplicated urinary tract infection was made in 2014. Nitrofurantoin resistance is reported in <6% urinary E. coli isolates in the UK, however, mechanisms underpinning nitrofurantoin resistance in these isolates remain unknown. This study aimed to identify the genetic basis of nitrofurantoin resistance in urinary E. coli isolates collected from north west London and then elucidate resistance-associated genetic alterations in available UK E. coli genomes. As a result, an algorithm was developed to predict nitrofurantoin susceptibility. Deleterious mutations and gene-inactivating insertion sequences in chromosomal nitroreductase genes nfsA and/or nfsB were identified in genomes of nine confirmed nitrofurantoin-resistant urinary E. coli isolates and additional 11 E. coli isolates that were highlighted by the prediction algorithm and subsequently validated to be nitrofurantoin-resistant. Eight categories of allelic changes in nfsA, nfsB, and the associated gene ribE were detected in 12412 E. coli genomes from the UK. Evolutionary analysis of these three genes revealed homoplasic mutations and explained the previously reported order of stepwise mutations. The mobile gene complex oqxAB, which is associated with reduced nitrofurantoin susceptibility, was identified in only one of the 12412 genomes. In conclusion, mutations and insertion sequences in nfsA and nfsB were leading causes of nitrofurantoin resistance in UK E. coli. As nitrofurantoin exposure increases in human populations, the prevalence of nitrofurantoin resistance in carriage E. coli isolates and those from urinary and bloodstream infections should be monitored.
Collapse
Affiliation(s)
- Yu Wan
- NIHR Health Protection Research Unit in Healthcare Associated Infections and Antimicrobial Resistance, Department of Infectious Disease, Imperial College London, London, United Kingdom
| | - Ewurabena Mills
- NIHR Health Protection Research Unit in Healthcare Associated Infections and Antimicrobial Resistance, Department of Infectious Disease, Imperial College London, London, United Kingdom
- MRC Centre for Molecular Bacteriology and Infection, Imperial College London, London, United Kingdom
| | - Rhoda C.Y. Leung
- NIHR Health Protection Research Unit in Healthcare Associated Infections and Antimicrobial Resistance, Department of Infectious Disease, Imperial College London, London, United Kingdom
- Present address: Department of Microbiology, Queen Mary Hospital, Hong Kong S.A.R., PR China
| | - Ana Vieira
- NIHR Health Protection Research Unit in Healthcare Associated Infections and Antimicrobial Resistance, Department of Infectious Disease, Imperial College London, London, United Kingdom
- MRC Centre for Molecular Bacteriology and Infection, Imperial College London, London, United Kingdom
| | - Xiangyun Zhi
- NIHR Health Protection Research Unit in Healthcare Associated Infections and Antimicrobial Resistance, Department of Infectious Disease, Imperial College London, London, United Kingdom
- MRC Centre for Molecular Bacteriology and Infection, Imperial College London, London, United Kingdom
| | - Nicholas J. Croucher
- Department of Infectious Disease Epidemiology, School of Public Health, Imperial College London, London, United Kingdom
- MRC Centre for Global Infectious Disease Analysis, School of Public Health, Imperial College London, London, United Kingdom
| | - Neil Woodford
- Antimicrobial Resistance and Healthcare Associated Infections Reference Unit, National Infection Service, Public Health England, Colindale, London, United Kingdom
| | - Elita Jauneikaite
- NIHR Health Protection Research Unit in Healthcare Associated Infections and Antimicrobial Resistance, Department of Infectious Disease, Imperial College London, London, United Kingdom
- Department of Infectious Disease Epidemiology, School of Public Health, Imperial College London, London, United Kingdom
- MRC Centre for Global Infectious Disease Analysis, School of Public Health, Imperial College London, London, United Kingdom
| | - Matthew J. Ellington
- NIHR Health Protection Research Unit in Healthcare Associated Infections and Antimicrobial Resistance, Department of Infectious Disease, Imperial College London, London, United Kingdom
- Antimicrobial Resistance and Healthcare Associated Infections Reference Unit, National Infection Service, Public Health England, Colindale, London, United Kingdom
| | - Shiranee Sriskandan
- NIHR Health Protection Research Unit in Healthcare Associated Infections and Antimicrobial Resistance, Department of Infectious Disease, Imperial College London, London, United Kingdom
- MRC Centre for Molecular Bacteriology and Infection, Imperial College London, London, United Kingdom
- *Correspondence: Shiranee Sriskandan,
| |
Collapse
|
15
|
Chaguza C, Tonkin-Hill G, Lo SW, Hadfield J, Croucher NJ, Harris SR, Bentley SD. RCandy: an R package for visualizing homologous recombinations in bacterial genomes. Bioinformatics 2021; 38:1450-1451. [PMID: 34864895 PMCID: PMC8826011 DOI: 10.1093/bioinformatics/btab814] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2021] [Revised: 11/25/2021] [Accepted: 11/29/2021] [Indexed: 01/05/2023] Open
Abstract
SUMMARY Homologous recombination is an important evolutionary process in bacteria and other prokaryotes, which increases genomic sequence diversity and can facilitate adaptation. Several methods and tools have been developed to detect genomic regions recently affected by recombination. Exploration and visualization of such recombination events can reveal valuable biological insights, but it remains challenging. Here, we present RCandy, a platform-independent R package for rapid, simple and flexible visualization of recombination events in bacterial genomes. AVAILABILITY AND IMPLEMENTATION RCandy is an R package freely available for use under the MIT license. It is platform-independent and has been tested on Windows, Linux and MacOSX. The source code comes together with a detailed vignette available on GitHub at https://github.com/ChrispinChaguza/RCandy. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
| | - Gerry Tonkin-Hill
- Parasites and Microbes Programme, Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton, Cambridge, UK
| | - Stephanie W Lo
- Parasites and Microbes Programme, Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton, Cambridge, UK
| | - James Hadfield
- Vaccine and Infectious Disease Division, Fred Hutchinson Cancer Research Center, Seattle, WA, USA
| | - Nicholas J Croucher
- Department of Infectious Disease Epidemiology, Imperial College London, London, UK
| | - Simon R Harris
- Microbiotica Ltd, Biodata Innovation Centre, Wellcome Genome Campus, Hinxton, UK
| | - Stephen D Bentley
- Parasites and Microbes Programme, Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton, Cambridge, UK
| |
Collapse
|
16
|
D'Aeth JC, van der Linden MPG, McGee L, de Lencastre H, Turner P, Song JH, Lo SW, Gladstone RA, Sá-Leão R, Ko KS, Hanage WP, Breiman RF, Beall B, Bentley SD, Croucher NJ. The role of interspecies recombination in the evolution of antibiotic-resistant pneumococci. eLife 2021; 10:e67113. [PMID: 34259624 PMCID: PMC8321556 DOI: 10.7554/elife.67113] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2021] [Accepted: 04/16/2021] [Indexed: 12/27/2022] Open
Abstract
Multidrug-resistant Streptococcus pneumoniae emerge through the modification of core genome loci by interspecies homologous recombinations, and acquisition of gene cassettes. Both occurred in the otherwise contrasting histories of the antibiotic-resistant S. pneumoniae lineages PMEN3 and PMEN9. A single PMEN3 clade spread globally, evading vaccine-induced immunity through frequent serotype switching, whereas locally circulating PMEN9 clades independently gained resistance. Both lineages repeatedly integrated Tn916-type and Tn1207.1-type elements, conferring tetracycline and macrolide resistance, respectively, through homologous recombination importing sequences originating in other species. A species-wide dataset found over 100 instances of such interspecific acquisitions of resistance cassettes and flanking homologous arms. Phylodynamic analysis of the most commonly sampled Tn1207.1-type insertion in PMEN9, originating from a commensal and disrupting a competence gene, suggested its expansion across Germany was driven by a high ratio of macrolide-to-β-lactam consumption. Hence, selection from antibiotic consumption was sufficient for these atypically large recombinations to overcome species boundaries across the pneumococcal chromosome.
Collapse
Affiliation(s)
- Joshua C D'Aeth
- MRC Centre for Global Infectious Disease Analysis, Department of Infectious Disease Epidemiology, Imperial College LondonLondonUnited Kingdom
| | - Mark PG van der Linden
- Institute for Medical Microbiology, National Reference Center for Streptococci, University Hospital RWTH AachenAachenGermany
| | - Lesley McGee
- Respiratory Diseases Branch, Centers for Disease Control and PreventionAtlantaUnited States
| | - Herminia de Lencastre
- Laboratory of Molecular Genetics, Instituto de Tecnologia Química e Biológica, Universidade Nova de LisboaOeirasPortugal
- Laboratory of Microbiology and Infectious Diseases, The Rockefeller UniversityNew YorkUnited States
| | - Paul Turner
- Cambodia Oxford Medical Research Unit, Angkor Hospital for ChildrenSiem ReapCambodia
- Centre for Tropical Medicine and Global Health, Nuffield Department of Medicine, University of OxfordOxfordUnited Kingdom
| | - Jae-Hoon Song
- Department of Molecular Cell Biology, Sungkyunkwan University School of MedicineSuwonRepublic of Korea
| | - Stephanie W Lo
- Parasites & Microbes, Wellcome Sanger Institute, Wellcome Genome CampusHinxtonUnited Kingdom
| | - Rebecca A Gladstone
- Parasites & Microbes, Wellcome Sanger Institute, Wellcome Genome CampusHinxtonUnited Kingdom
| | - Raquel Sá-Leão
- Laboratory of Molecular Microbiology of Human Pathogens, Instituto de Tecnologia Química e Biológica, Universidade Nova de LisboaOeirasPortugal
| | - Kwan Soo Ko
- Department of Molecular Cell Biology, Sungkyunkwan University School of MedicineSuwonRepublic of Korea
| | - William P Hanage
- Center for Communicable Disease Dynamics, Harvard T.H. Chan School of Public HealthBostonUnited States
| | - Robert F Breiman
- Department of Global Health, Rollins School of Public Health, Emory UniversityAtlantaUnited States
| | - Bernard Beall
- Respiratory Diseases Branch, Centers for Disease Control and PreventionAtlantaUnited States
| | - Stephen D Bentley
- Parasites & Microbes, Wellcome Sanger Institute, Wellcome Genome CampusHinxtonUnited Kingdom
| | - Nicholas J Croucher
- MRC Centre for Global Infectious Disease Analysis, Department of Infectious Disease Epidemiology, Imperial College LondonLondonUnited Kingdom
| |
Collapse
|
17
|
Gladstone RA, McNally A, Pöntinen AK, Tonkin-Hill G, Lees JA, Skytén K, Cléon F, Christensen MOK, Haldorsen BC, Bye KK, Gammelsrud KW, Hjetland R, Kümmel A, Larsen HE, Lindemann PC, Löhr IH, Marvik Å, Nilsen E, Noer MT, Simonsen GS, Steinbakk M, Tofteland S, Vattøy M, Bentley SD, Croucher NJ, Parkhill J, Johnsen PJ, Samuelsen Ø, Corander J. Emergence and dissemination of antimicrobial resistance in Escherichia coli causing bloodstream infections in Norway in 2002-17: a nationwide, longitudinal, microbial population genomic study. Lancet Microbe 2021; 2:e331-e341. [PMID: 35544167 PMCID: PMC7614948 DOI: 10.1016/s2666-5247(21)00031-8] [Citation(s) in RCA: 30] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/13/2020] [Revised: 02/01/2021] [Accepted: 02/03/2021] [Indexed: 01/23/2023]
Abstract
BACKGROUND The clonal diversity underpinning trends in multidrug resistant Escherichia coli causing bloodstream infections remains uncertain. We aimed to determine the contribution of individual clones to resistance over time, using large-scale genomics-based molecular epidemiology. METHODS This was a longitudinal, E coli population, genomic, cohort study that sampled isolates from 22 512 E coli bloodstream infections included in the Norwegian surveillance programme on resistant microbes (NORM) from 2002 to 2017. 15 of 22 laboratories were able to share their isolates, and the first 22·5% of isolates from each year were requested. We used whole genome sequencing to infer the population structure (PopPUNK), and we investigated the clade composition of the dominant multidrug resistant clonal complex (CC)131 using genetic markers previously reported for sequence type (ST)131, effective population size (BEAST), and presence of determinants of antimicrobial resistance (ARIBA, PointFinder, and ResFinder databases) over time. We compared these features between the 2002-10 and 2011-17 time periods. We also compared our results with those of a longitudinal study from the UK done between 2001 and 2011. FINDINGS Of the 3500 isolates requested from the participating laboratories, 3397 (97·1%) were received, of which 3254 (95·8%) were successfully sequenced and included in the analysis. A significant increase in the number of multidrug resistant CC131 isolates from 71 (5·6%) of 1277 in 2002-10 to 207 (10·5%) of 1977 in 2011-17 (p<0·0001), was the largest clonal expansion. CC131 was the most common clone in extended-spectrum β-lactamase (ESBL)-positive isolates (75 [58·6%] of 128) and fluoroquinolone non-susceptible isolates (148 [39·2%] of 378). Within CC131, clade A increased in prevalence from 2002, whereas the global multidrug resistant clade C2 was not observed until 2007. Multiple de-novo acquisitions of both blaCTX-M ESBL-encoding genes in clades A and C1 and gain of phenotypic fluoroquinolone non-susceptibility across the clade A phylogeny were observed. We estimated that exponential increases in the effective population sizes of clades A, C1, and C2 occurred in the mid-2000s, and in clade B a decade earlier. The rate of increase in the estimated effective population size of clade A (Ne=3147) was nearly ten-times that of C2 (Ne=345), with clade A over-represented in Norwegian CC131 isolates (75 [27·0%] of 278) compared with the UK study (8 [5·4%] of 147 isolates). INTERPRETATION The early and sustained establishment of predominantly antimicrobial susceptible CC131 clade A isolates, relative to multidrug resistant clade C2 isolates, suggests that resistance is not necessary for clonal success. However, even in the low antibiotic use setting of Norway, resistance to important antimicrobial classes has rapidly been selected for in CC131 clade A isolates. This study shows the importance of genomic surveillance in uncovering the complex ecology underlying multidrug resistance dissemination and competition, which have implications for the design of strategies and interventions to control the spread of high-risk multidrug resistant clones. FUNDING Trond Mohn Foundation, European Research Council, Marie Skłodowska-Curie Actions, and the Wellcome Trust.
Collapse
Affiliation(s)
| | - Alan McNally
- Institute of Microbiology and Infection, University of Birmingham, Birmingham, UK
| | - Anna K Pöntinen
- Department of Biostatistics, University of Oslo, Oslo, Norway
| | | | - John A Lees
- Faculty of Medicine, School of Public Health, Imperial College, London, UK
| | - Kusti Skytén
- Department of Biostatistics, University of Oslo, Oslo, Norway
| | - François Cléon
- Department of Pharmacy, Faculty of Health Sciences UiT The Arctic University of Norway, Tromsø, Norway
| | - Martin O K Christensen
- Norwegian National Advisory Unit on Detection of Antimicrobial Resistance, Department of Microbiology and Infection Control, University Hospital of North Norway, Tromsø, Norway
| | - Bjørg C Haldorsen
- Norwegian National Advisory Unit on Detection of Antimicrobial Resistance, Department of Microbiology and Infection Control, University Hospital of North Norway, Tromsø, Norway
| | - Kristina K Bye
- Laboratory of Microbiology, Department of Medical Biochemistry, Oslo University Hospital Radiumhospitalet, Oslo, Norway
| | - Karianne W Gammelsrud
- Department of Microbiology, Division of Laboratory Medicine, Oslo University Hospital Ullevål, Oslo, Norway
| | - Reidar Hjetland
- Department of Microbiology, Førde General Hospital, Førde Health Trust, Førde, Norway
| | - Angela Kümmel
- Department of Laboratory Medicine, Levanger Hospital, Nord-Trøndelag Hospital Trust, Levanger, Norway
| | - Hege E Larsen
- Department of Microbiology, Nordland Hospital, Bodø, Norway
| | | | - Iren H Löhr
- Department of Medical Microbiology, Stavanger University Hospital, Stavanger, Norway
| | - Åshild Marvik
- Department of Microbiology, Vestfold Hospital, Tønsberg, Norway
| | - Einar Nilsen
- Department of Microbiology, Moere and Romsdal Hospital Trust, Molde, Norway
| | - Marie T Noer
- Institute of Clinical Medicine, Faculty of Medicine, University of Oslo, Oslo, Norway; Institute of Medical Microbiology, Oslo University Hospital Rikshospitalet, Oslo, Norway
| | - Gunnar S Simonsen
- Department of Medical Biology, Faculty of Health Sciences UiT The Arctic University of Norway, Tromsø, Norway; Norwegian National Advisory Unit on Detection of Antimicrobial Resistance, Department of Microbiology and Infection Control, University Hospital of North Norway, Tromsø, Norway; Norwegian Institute of Public Health, Oslo, Norway
| | - Martin Steinbakk
- Centre for Laboratory Medicine, Sections for Microbiology, Østfold Hospital, Kalnes, Norway
| | - Ståle Tofteland
- Department of Medical Microbiology, Sørlandet Hospital, Kristiansand, Norway
| | - Marit Vattøy
- Department of Microbiology, Akershus University Hospital, Lørenskog, Norway
| | | | | | - Julian Parkhill
- Department of Veterinary Medicine, University of Cambridge, Cambridge, UK
| | - Pål J Johnsen
- Department of Pharmacy, Faculty of Health Sciences UiT The Arctic University of Norway, Tromsø, Norway
| | - Ørjan Samuelsen
- Department of Pharmacy, Faculty of Health Sciences UiT The Arctic University of Norway, Tromsø, Norway; Norwegian National Advisory Unit on Detection of Antimicrobial Resistance, Department of Microbiology and Infection Control, University Hospital of North Norway, Tromsø, Norway
| | - Jukka Corander
- Department of Biostatistics, University of Oslo, Oslo, Norway; Parasites and Microbes, Wellcome Sanger Institute, Cambridge, UK
| |
Collapse
|
18
|
Harrow GL, Lees JA, Hanage WP, Lipsitch M, Corander J, Colijn C, Croucher NJ. Negative frequency-dependent selection and asymmetrical transformation stabilise multi-strain bacterial population structures. ISME J 2021; 15:1523-1538. [PMID: 33408365 PMCID: PMC8115253 DOI: 10.1038/s41396-020-00867-w] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/11/2020] [Revised: 12/01/2020] [Accepted: 12/03/2020] [Indexed: 02/06/2023]
Abstract
Streptococcus pneumoniae can be divided into many strains, each a distinct set of isolates sharing similar core and accessory genomes, which co-circulate within the same hosts. Previous analyses suggested the short-term vaccine-associated dynamics of S. pneumoniae strains may be mediated through multi-locus negative frequency-dependent selection (NFDS), which maintains accessory loci at equilibrium frequencies. Long-term simulations demonstrated NFDS stabilised clonally-evolving multi-strain populations through preventing the loss of variation through drift, based on polymorphism frequencies, pairwise genetic distances and phylogenies. However, allowing symmetrical recombination between isolates evolving under multi-locus NFDS generated unstructured populations of diverse genotypes. Replication of the observed data improved when multi-locus NFDS was combined with recombination that was instead asymmetrical, favouring deletion of accessory loci over insertion. This combination separated populations into strains through outbreeding depression, resulting from recombinants with reduced accessory genomes having lower fitness than their parental genotypes. Although simplistic modelling of recombination likely limited these simulations' ability to maintain some properties of genomic data as accurately as those lacking recombination, the combination of asymmetrical recombination and multi-locus NFDS could restore multi-strain population structures from randomised initial populations. As many bacteria inhibit insertions into their chromosomes, this combination may commonly underlie the co-existence of strains within a niche.
Collapse
Affiliation(s)
- Gabrielle L Harrow
- MRC Centre for Global Infectious Disease Analysis, Department of Infectious Disease Epidemiology, Imperial College London, Norfolk Place, London, W2 1PG, UK
| | - John A Lees
- MRC Centre for Global Infectious Disease Analysis, Department of Infectious Disease Epidemiology, Imperial College London, Norfolk Place, London, W2 1PG, UK
| | - William P Hanage
- Center for Communicable Disease Dynamics, Department of Epidemiology, Harvard T.H. Chan School of Public Health, 677 Huntington Avenue, Boston, MA, 02115, USA
| | - Marc Lipsitch
- Center for Communicable Disease Dynamics, Department of Epidemiology, Harvard T.H. Chan School of Public Health, 677 Huntington Avenue, Boston, MA, 02115, USA
| | - Jukka Corander
- Department of Biostatistics, University of Oslo, Oslo, Norway
- Helsinki Institute of Information Technology, Department of Mathematics and Statistics, University of Helsinki, Helsinki, Finland
- Parasites & Microbes Programme, Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton, Cambridge, UK
| | - Caroline Colijn
- Parasites & Microbes Programme, Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton, Cambridge, UK
- Department of Mathematics, Simon Fraser University, Burnaby, BC, Canada
| | - Nicholas J Croucher
- MRC Centre for Global Infectious Disease Analysis, Department of Infectious Disease Epidemiology, Imperial College London, Norfolk Place, London, W2 1PG, UK.
| |
Collapse
|
19
|
Zamudio R, Haigh RD, Ralph JD, De Ste Croix M, Tasara T, Zurfluh K, Kwun MJ, Millard AD, Bentley SD, Croucher NJ, Stephan R, Oggioni MR. Lineage-specific evolution and gene flow in Listeria monocytogenes are independent of bacteriophages. Environ Microbiol 2020; 22:5058-5072. [PMID: 32483914 PMCID: PMC7614921 DOI: 10.1111/1462-2920.15111] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2020] [Revised: 05/26/2020] [Accepted: 05/30/2020] [Indexed: 11/29/2022]
Abstract
Listeria monocytogenes is a foodborne pathogen causing systemic infection with high mortality. To allow efficient tracking of outbreaks a clear definition of the genomic signature of a cluster of related isolates is required, but lineage-specific characteristics call for a more detailed understanding of evolution. In our work, we used core genome MLST (cgMLST) to identify new outbreaks combined to core genome SNP analysis to characterize the population structure and gene flow between lineages. Whilst analysing differences between the four lineages of L. monocytogenes we have detected differences in the recombination rate, and interestingly also divergence in the SNP differences between sub-lineages. In addition, the exchange of core genome variation between the lineages exhibited a distinct pattern, with lineage III being the best donor for horizontal gene transfer. Whilst attempting to link bacteriophage-mediated transduction to observed gene transfer, we found an inverse correlation between phage presence in a lineage and the extent of recombination. Irrespective of the profound differences in recombination rates observed between sub-lineages and lineages, we found that the previously proposed cut-off of 10 allelic differences in cgMLST can be still considered valid for the definition of a foodborne outbreak cluster of L. monocytogenes.
Collapse
Affiliation(s)
- Roxana Zamudio
- Department of Genetics and Genome Biology, University of Leicester, Leicester, UK
| | - Richard D Haigh
- Department of Respiratory Sciences, University of Leicester, Leicester, UK
| | - Joseph D Ralph
- Department of Genetics and Genome Biology, University of Leicester, Leicester, UK
| | - Megan De Ste Croix
- Department of Genetics and Genome Biology, University of Leicester, Leicester, UK
| | - Taurai Tasara
- Institute for Food Safety and Hygiene, Vetsuisse Faculty, University of Zurich, Zurich, Switzerland
| | - Katrin Zurfluh
- Institute for Food Safety and Hygiene, Vetsuisse Faculty, University of Zurich, Zurich, Switzerland
| | - Min Jung Kwun
- MRC Centre for Global Infectious Disease Analysis, Department of Infectious Disease Epidemiology, Imperial College London, London, UK
| | - Andrew D Millard
- Department of Genetics and Genome Biology, University of Leicester, Leicester, UK
| | - Stephen D Bentley
- Parasites and Microbes, Wellcome Trust Sanger Institute, Hinxton, UK
| | - Nicholas J Croucher
- MRC Centre for Global Infectious Disease Analysis, Department of Infectious Disease Epidemiology, Imperial College London, London, UK
| | - Roger Stephan
- Institute for Food Safety and Hygiene, Vetsuisse Faculty, University of Zurich, Zurich, Switzerland
| | - Marco R Oggioni
- Department of Genetics and Genome Biology, University of Leicester, Leicester, UK
| |
Collapse
|
20
|
Løchen A, Croucher NJ, Anderson RM. Divergent serotype replacement trends and increasing diversity in pneumococcal disease in high income settings reduce the benefit of expanding vaccine valency. Sci Rep 2020; 10:18977. [PMID: 33149149 PMCID: PMC7643077 DOI: 10.1038/s41598-020-75691-5] [Citation(s) in RCA: 61] [Impact Index Per Article: 15.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2020] [Accepted: 10/08/2020] [Indexed: 11/25/2022] Open
Abstract
Streptococcus pneumoniae is a significant cause of otitis media, pneumonia, and meningitis. Only seven of the approximately 100 serotypes were initially included in the pneumococcal polysaccharide conjugate vaccine (PCV) in 2000 before it was expanded in subsequent years. Although the invasive pneumococcal disease (IPD) incidence due to vaccine serotypes (VT) has declined, partial replacement by non-vaccine serotypes (NVT) was observed following widespread vaccine uptake. We conducted a trend analysis assembling the available evidence for PCV impact on European, North American and Australian national IPD. Significant effectiveness against VT IPD in infants was observed, although the impact on national IPD incidence varied internationally due to serotype replacement. Currently, NVT serotypes 8, 9N, 15A and 23B are increasing in the countries assessed, although a variety of other NVTs are affecting each country and age group. Despite these common emerging serotypes, there has not been a dominant IPD serotype post-vaccination as there was pre-vaccination (serotype 14) or post-PCV7 (serotype 19A), suggesting that future vaccines with additional serotypes will be less effective at targeting and reducing IPD in global populations than previous PCVs. The rise of diverse NVTs in all settings’ top-ranked IPD-causing serotypes emphasizes the urgent need for surveillance data on serotype distribution and serotype-specific invasiveness post-vaccination to facilitate decision making concerning both expanding current vaccination programmes and increasing vaccine valency.
Collapse
Affiliation(s)
- Alessandra Løchen
- Department of Infectious Disease Epidemiology, School of Public Health, Imperial College, London, UK.,MRC Centre for Global Infectious Disease Analysis, School of Public Health, Imperial College London, Norfolk Place, London, W2 1PG, UK
| | - Nicholas J Croucher
- Department of Infectious Disease Epidemiology, School of Public Health, Imperial College, London, UK. .,MRC Centre for Global Infectious Disease Analysis, School of Public Health, Imperial College London, Norfolk Place, London, W2 1PG, UK.
| | - Roy M Anderson
- Department of Infectious Disease Epidemiology, School of Public Health, Imperial College, London, UK.,MRC Centre for Global Infectious Disease Analysis, School of Public Health, Imperial College London, Norfolk Place, London, W2 1PG, UK
| |
Collapse
|
21
|
Azarian T, Martinez PP, Arnold BJ, Qiu X, Grant LR, Corander J, Fraser C, Croucher NJ, Hammitt LL, Reid R, Santosham M, Weatherholtz RC, Bentley SD, O’Brien KL, Lipsitch M, Hanage WP. Frequency-dependent selection can forecast evolution in Streptococcus pneumoniae. PLoS Biol 2020; 18:e3000878. [PMID: 33091022 PMCID: PMC7580979 DOI: 10.1371/journal.pbio.3000878] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2020] [Accepted: 09/18/2020] [Indexed: 11/30/2022] Open
Abstract
Predicting how pathogen populations will change over time is challenging. Such has been the case with Streptococcus pneumoniae, an important human pathogen, and the pneumococcal conjugate vaccines (PCVs), which target only a fraction of the strains in the population. Here, we use the frequencies of accessory genes to predict changes in the pneumococcal population after vaccination, hypothesizing that these frequencies reflect negative frequency-dependent selection (NFDS) on the gene products. We find that the standardized predicted fitness of a strain, estimated by an NFDS-based model at the time the vaccine is introduced, enables us to predict whether the strain increases or decreases in prevalence following vaccination. Further, we are able to forecast the equilibrium post-vaccine population composition and assess the invasion capacity of emerging lineages. Overall, we provide a method for predicting the impact of an intervention on pneumococcal populations with potential application to other bacterial pathogens in which NFDS is a driving force.
Collapse
Affiliation(s)
- Taj Azarian
- Burnett School of Biomedical Sciences, University of Central Florida, Orlando, Florida, United States of America
- Center for Communicable Disease Dynamics, Department of Epidemiology, T.H. Chan School of Public Health, Harvard University, Boston, Massachusetts, United States of America
| | - Pamela P. Martinez
- Center for Communicable Disease Dynamics, Department of Epidemiology, T.H. Chan School of Public Health, Harvard University, Boston, Massachusetts, United States of America
| | - Brian J. Arnold
- Center for Communicable Disease Dynamics, Department of Epidemiology, T.H. Chan School of Public Health, Harvard University, Boston, Massachusetts, United States of America
| | - Xueting Qiu
- Center for Communicable Disease Dynamics, Department of Epidemiology, T.H. Chan School of Public Health, Harvard University, Boston, Massachusetts, United States of America
| | - Lindsay R. Grant
- Center for American Indian Health, Johns Hopkins Bloomberg School of Public Health, Baltimore, Maryland, United States of America
| | - Jukka Corander
- Helsinki Institute for Information Technology, Department of Mathematics and Statistics, University of Helsinki, Helsinki, Finland
- Department of Biostatistics, University of Oslo, Oslo, Norway
- Infection Genomics, The Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, United Kingdom
| | - Christophe Fraser
- Big Data Institute, Nuffield Department of Medicine, University of Oxford, Oxford, United Kingdom
| | - Nicholas J. Croucher
- MRC Centre for Global Infectious Disease Analysis, Department of Infectious Disease Epidemiology, Imperial College London, London, United Kingdom
| | - Laura L. Hammitt
- Center for American Indian Health, Johns Hopkins Bloomberg School of Public Health, Baltimore, Maryland, United States of America
| | - Raymond Reid
- Center for American Indian Health, Johns Hopkins Bloomberg School of Public Health, Baltimore, Maryland, United States of America
| | - Mathuram Santosham
- Center for American Indian Health, Johns Hopkins Bloomberg School of Public Health, Baltimore, Maryland, United States of America
| | - Robert C. Weatherholtz
- Center for American Indian Health, Johns Hopkins Bloomberg School of Public Health, Baltimore, Maryland, United States of America
| | - Stephen D. Bentley
- Infection Genomics, The Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, United Kingdom
| | | | - Marc Lipsitch
- Center for Communicable Disease Dynamics, Department of Epidemiology, T.H. Chan School of Public Health, Harvard University, Boston, Massachusetts, United States of America
- Department of Immunology and Infectious Diseases, T.H. Chan School of Public Health, Harvard University, Boston, Massachusetts, United States of America
| | - William P. Hanage
- Center for Communicable Disease Dynamics, Department of Epidemiology, T.H. Chan School of Public Health, Harvard University, Boston, Massachusetts, United States of America
| |
Collapse
|
22
|
Jeffrey B, Aanensen DM, Croucher NJ, Bhatt S. Predicting the future distribution of antibiotic resistance using time series forecasting and geospatial modelling. Wellcome Open Res 2020. [DOI: 10.12688/wellcomeopenres.16153.1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022] Open
Abstract
Background: Increasing antibiotic resistance in a location may be mitigated by changes in treatment policy, or interventions to limit transmission of resistant bacteria. Therefore, accurate forecasting of the distribution of antibiotic resistance could be advantageous. Two previously published studies addressed this, but neither study compared alternative forecasting algorithms or considered spatial patterns of resistance spread. Methods: We analysed data describing the annual prevalence of antibiotic resistance per country in Europe from 2012 – 2016, and the quarterly prevalence of antibiotic resistance per clinical commissioning group in England from 2015 – 2018. We combined these with data on rates of possible covariates of resistance. These data were used to compare the previously published forecasting models, with other commonly used forecasting models, including one geospatial model. Covariates were incorporated into the geospatial model to assess their relationship with antibiotic resistance. Results: For the European data, which was recorded on a coarse spatiotemporal scale, a naïve forecasting model was consistently the most accurate of any of the forecasting models tested. The geospatial model did not improve on this accuracy. However, it did provide some evidence that antibiotic consumption can partially explain the distribution of resistance. The English data were aggregated at a finer scale, and expected-trend-seasonal (ETS) forecasts were the most accurate. The geospatial model did not significantly improve upon the median accuracy of the ETS model, but it appeared to be less sensitive to noise in the data, and provided evidence that rates of antibiotic prescription and bacteraemia are correlated with resistance. Conclusion: Annual, national-level surveillance data appears to be insufficient for fitting accurate antibiotic resistance forecasting models, but there is evidence that data collected at a finer spatiotemporal scale could be used to improve forecast accuracy. Additionally, incorporating antibiotic prescription or consumption data into the model could improve the predictive accuracy.
Collapse
|
23
|
Lehtinen S, Chewapreecha C, Lees J, Hanage WP, Lipsitch M, Croucher NJ, Bentley SD, Turner P, Fraser C, Mostowy RJ. Horizontal gene transfer rate is not the primary determinant of observed antibiotic resistance frequencies in Streptococcus pneumoniae. Sci Adv 2020; 6:eaaz6137. [PMID: 32671212 PMCID: PMC7314567 DOI: 10.1126/sciadv.aaz6137] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/24/2019] [Accepted: 03/10/2020] [Indexed: 06/11/2023]
Abstract
The extent to which evolution is constrained by the rate at which horizontal gene transfer (HGT) allows DNA to move between genetic lineages is an open question, which we address in the context of antibiotic resistance in Streptococcus pneumoniae. We analyze microbiological, genomic, and epidemiological data from the largest-to-date sequenced pneumococcal carriage study in 955 infants from a refugee camp on the Thailand-Myanmar border. Using a unified framework, we simultaneously test prior hypotheses on rates of HGT and a key evolutionary covariate (duration of carriage) as determinants of resistance frequencies. We conclude that in this setting, there is little evidence of HGT playing a major role in determining resistance frequencies. Instead, observed resistance frequencies are best explained as the outcome of selection acting on a pool of variants, irrespective of the rate at which resistance determinants move between genetic lineages.
Collapse
Affiliation(s)
- Sonja Lehtinen
- Big Data Institute, University of Oxford, Oxford, UK
- Institute of Integrative Biology, Department of Environmental Systems Science, ETH Zurich, Zurich, Switzerland
| | - Claire Chewapreecha
- Wellcome Sanger Institute, Hinxton, UK
- Mahidol-Oxford Tropical Medicine Research Unit, Bangkok, Thailand
- Bioinformatics and Systems Biology Program, School of Bioresource and Technology, King Mongkut’s University of Technology Thonburi, Bangkok, Thailand
| | - John Lees
- Department of Infectious Disease Epidemiology, School of Public Health, Imperial College London, London, UK
| | - William P. Hanage
- Center for Communicable Disease Dynamics, Department of Epidemiology, Harvard T. H. Chan School of Public Health, Harvard University, Boston, MA, USA
| | - Marc Lipsitch
- Center for Communicable Disease Dynamics, Department of Epidemiology, Harvard T. H. Chan School of Public Health, Harvard University, Boston, MA, USA
| | - Nicholas J. Croucher
- Center for Communicable Disease Dynamics, Department of Epidemiology, Harvard T. H. Chan School of Public Health, Harvard University, Boston, MA, USA
| | | | - Paul Turner
- Cambodia Oxford Medical Research Unit, Angkor Hospital for Children, Siem Reap, Cambodia
- Centre for Tropical Medicine and Global Health, University of Oxford, Oxford, UK
| | | | - Rafał J. Mostowy
- Big Data Institute, University of Oxford, Oxford, UK
- Department of Infectious Disease Epidemiology, School of Public Health, Imperial College London, London, UK
- Malopolska Centre of Biotechnology, Jagiellonian University, Krakow, Poland
| |
Collapse
|
24
|
Gladstone RA, Lo SW, Goater R, Yeats C, Taylor B, Hadfield J, Lees JA, Croucher NJ, van Tonder AJ, Bentley LJ, Quah FX, Blaschke AJ, Pershing NL, Byington CL, Balaji V, Hryniewicz W, Sigauque B, Ravikumar K, Almeida SCG, Ochoa TJ, Ho PL, du Plessis M, Ndlangisa KM, Cornick JE, Kwambana-Adams B, Benisty R, Nzenze SA, Madhi SA, Hawkins PA, Pollard AJ, Everett DB, Antonio M, Dagan R, Klugman KP, von Gottberg A, Metcalf BJ, Li Y, Beall BW, McGee L, Breiman RF, Aanensen DM, Bentley SD. Visualizing variation within Global Pneumococcal Sequence Clusters (GPSCs) and country population snapshots to contextualize pneumococcal isolates. Microb Genom 2020; 6:e000357. [PMID: 32375991 PMCID: PMC7371119 DOI: 10.1099/mgen.0.000357] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2019] [Accepted: 03/03/2020] [Indexed: 11/21/2022] Open
Abstract
Knowledge of pneumococcal lineages, their geographic distribution and antibiotic resistance patterns, can give insights into global pneumococcal disease. We provide interactive bioinformatic outputs to explore such topics, aiming to increase dissemination of genomic insights to the wider community, without the need for specialist training. We prepared 12 country-specific phylogenetic snapshots, and international phylogenetic snapshots of 73 common Global Pneumococcal Sequence Clusters (GPSCs) previously defined using PopPUNK, and present them in Microreact. Gene presence and absence defined using Roary, and recombination profiles derived from Gubbins are presented in Phandango for each GPSC. Temporal phylogenetic signal was assessed for each GPSC using BactDating. We provide examples of how such resources can be used. In our example use of a country-specific phylogenetic snapshot we determined that serotype 14 was observed in nine unrelated genetic backgrounds in South Africa. The international phylogenetic snapshot of GPSC9, in which most serotype 14 isolates from South Africa were observed, highlights that there were three independent sub-clusters represented by South African serotype 14 isolates. We estimated from the GPSC9-dated tree that the sub-clusters were each established in South Africa during the 1980s. We show how recombination plots allowed the identification of a 20 kb recombination spanning the capsular polysaccharide locus within GPSC97. This was consistent with a switch from serotype 6A to 19A estimated to have occured in the 1990s from the GPSC97-dated tree. Plots of gene presence/absence of resistance genes (tet, erm, cat) across the GPSC23 phylogeny were consistent with acquisition of a composite transposon. We estimated from the GPSC23-dated tree that the acquisition occurred between 1953 and 1975. Finally, we demonstrate the assignment of GPSC31 to 17 externally generated pneumococcal serotype 1 assemblies from Utah via Pathogenwatch. Most of the Utah isolates clustered within GPSC31 in a USA-specific clade with the most recent common ancestor estimated between 1958 and 1981. The resources we have provided can be used to explore to data, test hypothesis and generate new hypotheses. The accessible assignment of GPSCs allows others to contextualize their own collections beyond the data presented here.
Collapse
Affiliation(s)
| | - Stephanie W. Lo
- Parasites and microbes, Wellcome Sanger InstituteHinxton, UK
| | - Richard Goater
- Centre for Genomic Pathogen Surveillance, Wellcome Genome CampusHinxton, UK
- Big Data Institute, Li Ka Shing Centre for Health Information and Discovery, University of Oxford, Oxford, UK
| | - Corin Yeats
- Centre for Genomic Pathogen Surveillance, Wellcome Genome CampusHinxton, UK
- Big Data Institute, Li Ka Shing Centre for Health Information and Discovery, University of Oxford, Oxford, UK
| | - Ben Taylor
- Centre for Genomic Pathogen Surveillance, Wellcome Genome CampusHinxton, UK
- Big Data Institute, Li Ka Shing Centre for Health Information and Discovery, University of Oxford, Oxford, UK
| | - James Hadfield
- Vaccine and Infectious Disease Division, Fred Hutchinson Cancer Research Center, Seattle, WA, USA
| | - John A. Lees
- Faculty of Medicine, School of Public Health, Imperial College London, UK
| | | | - Andries J. van Tonder
- Parasites and microbes, Wellcome Sanger InstituteHinxton, UK
- Department of Veterinary Medicine, University of Cambridge, Cambridge, UK
| | - Leon J. Bentley
- Parasites and microbes, Wellcome Sanger InstituteHinxton, UK
| | - Fu Xiang Quah
- Parasites and microbes, Wellcome Sanger InstituteHinxton, UK
| | - Anne J. Blaschke
- Division of Pediatric Infectious Diseases, Department of Pediatrics, School of Medicine, University of Utah, 295 Chipeta Way, Salt Lake City, UT, 84108, USA
| | - Nicole L. Pershing
- Division of Pediatric Infectious Diseases, Department of Pediatrics, School of Medicine, University of Utah, 295 Chipeta Way, Salt Lake City, UT, 84108, USA
| | | | | | - Waleria Hryniewicz
- National Medicines Institute, Division of Clinical Microbiology and Infection Prevention, Warsaw, Poland
| | - Betuel Sigauque
- Fundação Manhiça / Centro de Investigação em Saúde da Manhiça (CISM), Maputo Mozambique, Instituto Nacional de Saúde, inistério de Saúde, Maputo, Mozambique
| | - K.L. Ravikumar
- Central Research Laboratory, Department of Microbiology, Kempegowda Institute of Medical Sciences Hospital & Research Center, Bangalore, India
| | | | - Theresa J. Ochoa
- Instituto de Medicina Tropical, Universidad Peruana Cayetano Heredia, Lima, Peru
| | - Pak Leung Ho
- Department of Microbiology and Carol Yu Centre for Infection, The University of Hong Kong, Queen Mary Hospital, Hong Kong, PR China
| | - Mignon du Plessis
- Centre for Respiratory Diseases and Meningitis, National Institute for Communicable Diseases, Johannesburg, South Africa
| | - Kedibone M. Ndlangisa
- Centre for Respiratory Diseases and Meningitis, National Institute for Communicable Diseases, Johannesburg, South Africa
| | | | - Brenda Kwambana-Adams
- NIHR Global Health Research Unit on Mucosal Pathogens, Division of Infection and Immunity, University College London, London, UK
- WHO Collaborating Centre for New Vaccines Surveillance, Medical Research Council Unit The Gambia at The London School of Hygiene & Tropical Medicine, Fajara, The Gambia
| | - Rachel Benisty
- The Faculty of Health Sciences, Ben-Gurion University of the NegevBeer-Sheva, Israel
| | - Susan A. Nzenze
- Medical Research Council: Respiratory and Meningeal Pathogens Research Unit, University of the Witwatersrand, Johannesburg, South Africa
- Department of Science and Technology/National Research Foundation: Vaccine Preventable Diseases, University of the Witwatersrand, Johannesburg, South Africa
| | - Shabir A. Madhi
- Medical Research Council: Respiratory and Meningeal Pathogens Research Unit, University of the Witwatersrand, Johannesburg, South Africa
- Department of Science and Technology/National Research Foundation: Vaccine Preventable Diseases, University of the Witwatersrand, Johannesburg, South Africa
| | | | - Andrew J. Pollard
- Oxford Vaccine Group, Department of Paediatrics, University of Oxford, and the NIHR Oxford Biomedical Research Centre, Oxford, UK
| | | | - Martin Antonio
- WHO Collaborating Centre for New Vaccines Surveillance, Medical Research Council Unit The Gambia at The London School of Hygiene & Tropical Medicine, Fajara, The Gambia
| | - Ron Dagan
- The Faculty of Health Sciences, Ben-Gurion University of the NegevBeer-Sheva, Israel
| | | | - Anne von Gottberg
- Department of Microbiology and Carol Yu Centre for Infection, The University of Hong Kong, Queen Mary Hospital, Hong Kong, PR China
| | | | - Yuan Li
- Centers for Disease Control and Prevention, Atlanta, GA, USA
| | | | - Lesley McGee
- Centers for Disease Control and Prevention, Atlanta, GA, USA
| | - Robert F. Breiman
- Rollins School Public Health, Emory University, GA, USA
- Emory Global Health Institute, Atlanta, GA, USA
| | - David M. Aanensen
- Centre for Genomic Pathogen Surveillance, Wellcome Genome CampusHinxton, UK
- Big Data Institute, Li Ka Shing Centre for Health Information and Discovery, University of Oxford, Oxford, UK
| | | |
Collapse
|
25
|
Colijn C, Corander J, Croucher NJ. Designing ecologically optimized pneumococcal vaccines using population genomics. Nat Microbiol 2020; 5:473-485. [PMID: 32015499 PMCID: PMC7614922 DOI: 10.1038/s41564-019-0651-y] [Citation(s) in RCA: 27] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2019] [Accepted: 12/03/2019] [Indexed: 12/14/2022]
Abstract
Streptococcus pneumoniae (the pneumococcus) is a common nasopharyngeal commensal that can cause invasive pneumococcal disease (IPD). Each component of current protein-polysaccharide conjugate vaccines (PCVs) generally induces immunity specific to one of the approximately 100 pneumococcal serotypes, and typically eliminates it from carriage and IPD through herd immunity. Overall carriage rates remain stable owing to replacement by non-PCV serotypes. Consequently, the net change in IPD incidence is determined by the relative invasiveness of the pre- and post-PCV-carried pneumococcal populations. In the present study, we identified PCVs expected to minimize the post-vaccine IPD burden by applying Bayesian optimization to an ecological model of serotype replacement that integrated epidemiological and genomic data. We compared optimal formulations for reducing infant-only or population-wide IPD, and identified potential benefits to including non-conserved pneumococcal carrier proteins. Vaccines were also devised to minimize IPD resistant to antibiotic treatment, despite the ecological model assuming that resistance levels in the carried population would be preserved. We found that expanding infant-administered PCV valency is likely to result in diminishing returns, and that complementary pairs of infant- and adult-administered vaccines could be a superior strategy. PCV performance was highly dependent on the circulating pneumococcal population, further highlighting the advantages of a diversity of anti-pneumococcal vaccination strategies.
Collapse
Affiliation(s)
- Caroline Colijn
- Department of Mathematics, Simon Fraser University, Burnaby, BC, Canada.
- Department of Mathematics, Imperial College London, London, UK.
| | - Jukka Corander
- Department of Biostatistics, University of Oslo, Oslo, Norway
- Helsinki Institute of Information Technology, Department of Mathematics and Statistics, University of Helsinki, Helsinki, Finland
- Parasites & Microbes, Wellcome Sanger Institute, Wellcome Genome Campus, Cambridge, UK
| | - Nicholas J Croucher
- MRC Centre for Global Infectious Disease Analysis, Department of Infectious Disease Epidemiology, Imperial College London, London, UK
| |
Collapse
|
26
|
Cremers AJH, Mobegi FM, van der Gaast-de Jongh C, van Weert M, van Opzeeland FJ, Vehkala M, Knol MJ, Bootsma HJ, Välimäki N, Croucher NJ, Meis JF, Bentley S, van Hijum SAFT, Corander J, Zomer AL, Ferwerda G, de Jonge MI. The Contribution of Genetic Variation of Streptococcus pneumoniae to the Clinical Manifestation of Invasive Pneumococcal Disease. Clin Infect Dis 2020; 68:61-69. [PMID: 29788414 DOI: 10.1093/cid/ciy417] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2017] [Accepted: 05/10/2018] [Indexed: 01/02/2023] Open
Abstract
Background Different clinical manifestations of invasive pneumococcal disease (IPD) have thus far mainly been explained by patient characteristics. Here we studied the contribution of pneumococcal genetic variation to IPD phenotype. Methods The index cohort consisted of 349 patients admitted to 2 Dutch hospitals between 2000-2011 with pneumococcal bacteremia. We performed genome-wide association studies to identify pneumococcal lineages, genes, and allelic variants associated with 23 clinical IPD phenotypes. The identified associations were validated in a nationwide (n = 482) and a post-pneumococcal vaccination cohort (n = 121). The contribution of confirmed pneumococcal genotypes to the clinical IPD phenotype, relative to known clinical predictors, was tested by regression analysis. Results Among IPD patients, the presence of pneumococcal gene slaA was a nationwide confirmed independent predictor of meningitis (odds ratio [OR], 10.5; P = .001), as was sequence cluster 9 (serotype 7F: OR, 3.68; P = .057). A set of 4 pneumococcal genes co-located on a prophage was a confirmed independent predictor of 30-day mortality (OR, 3.4; P = .003). We could detect the pneumococcal variants of concern in these patients' blood samples. Conclusions In this study, knowledge of pneumococcal genotypic variants improved the clinical risk assessment for detrimental manifestations of IPD. This provides us with novel opportunities to target, anticipate, or avert the pathogenic effects related to particular pneumococcal variants, and indicates that information on pneumococcal genotype is important for the diagnostic and treatment strategy in IPD. Ongoing surveillance is warranted to monitor the clinical value of information on pneumococcal variants in dynamic microbial and susceptible host populations.
Collapse
Affiliation(s)
- Amelieke J H Cremers
- Section of Pediatric Infectious Diseases, Laboratory of Medical Immunology, Radboud Institute for Molecular Life Sciences, Nijmegen, The Netherlands.,Radboud Center for Infectious Diseases, Center for Molecular and Biomolecular Informatics, Radboudumc, Nijmegen, The Netherlands.,Department of Medical Microbiology, Center for Molecular and Biomolecular Informatics, Radboudumc, Nijmegen, The Netherlands
| | - Fredrick M Mobegi
- Section of Pediatric Infectious Diseases, Laboratory of Medical Immunology, Radboud Institute for Molecular Life Sciences, Nijmegen, The Netherlands.,Radboud Center for Infectious Diseases, Center for Molecular and Biomolecular Informatics, Radboudumc, Nijmegen, The Netherlands.,Bacterial Genomics Group, Center for Molecular and Biomolecular Informatics, Radboudumc, Nijmegen, The Netherlands
| | - Christa van der Gaast-de Jongh
- Section of Pediatric Infectious Diseases, Laboratory of Medical Immunology, Radboud Institute for Molecular Life Sciences, Nijmegen, The Netherlands.,Radboud Center for Infectious Diseases, Center for Molecular and Biomolecular Informatics, Radboudumc, Nijmegen, The Netherlands
| | - Michelle van Weert
- Section of Pediatric Infectious Diseases, Laboratory of Medical Immunology, Radboud Institute for Molecular Life Sciences, Nijmegen, The Netherlands.,Radboud Center for Infectious Diseases, Center for Molecular and Biomolecular Informatics, Radboudumc, Nijmegen, The Netherlands
| | - Fred J van Opzeeland
- Section of Pediatric Infectious Diseases, Laboratory of Medical Immunology, Radboud Institute for Molecular Life Sciences, Nijmegen, The Netherlands.,Radboud Center for Infectious Diseases, Center for Molecular and Biomolecular Informatics, Radboudumc, Nijmegen, The Netherlands
| | - Minna Vehkala
- Department of Mathematics and Statistics, University of Helsinki, Finland
| | - Mirjam J Knol
- Centre for Infectious Disease Control, National Institute for Public Health and the Environment (RIVM), Bilthoven, The Netherlands
| | - Hester J Bootsma
- Centre for Infectious Disease Control, National Institute for Public Health and the Environment (RIVM), Bilthoven, The Netherlands
| | - Niko Välimäki
- Department of Mathematics and Statistics, University of Helsinki, Finland
| | - Nicholas J Croucher
- Medical Research Council Centre for Outbreak Analysis and Modelling, Department of Infectious Disease Epidemiology, Imperial College London, United Kingdom
| | - Jacques F Meis
- Department of Medical Microbiology and Infectious Diseases, Canisius-Wilhelmina Hospital, Nijmegen, The Netherlands
| | - Stephen Bentley
- Wellcome Trust Sanger Institute, Pathogen Genomics Group, Hinxton, Cambridge, United Kingdom
| | - Sacha A F T van Hijum
- Radboud Center for Infectious Diseases, Center for Molecular and Biomolecular Informatics, Radboudumc, Nijmegen, The Netherlands.,Bacterial Genomics Group, Center for Molecular and Biomolecular Informatics, Radboudumc, Nijmegen, The Netherlands.,NIZO, Ede, The Netherlands
| | - Jukka Corander
- Department of Mathematics and Statistics, University of Helsinki, Finland.,Wellcome Trust Sanger Institute, Pathogen Genomics Group, Hinxton, Cambridge, United Kingdom.,Department of Biostatistics, University of Oslo, Norway
| | - Aldert L Zomer
- Department of Infectious Diseases and Immunology, Faculty of Veterinary Medicine, Utrecht University, The Netherlands
| | - Gerben Ferwerda
- Section of Pediatric Infectious Diseases, Laboratory of Medical Immunology, Radboud Institute for Molecular Life Sciences, Nijmegen, The Netherlands.,Radboud Center for Infectious Diseases, Center for Molecular and Biomolecular Informatics, Radboudumc, Nijmegen, The Netherlands
| | - Marien I de Jonge
- Section of Pediatric Infectious Diseases, Laboratory of Medical Immunology, Radboud Institute for Molecular Life Sciences, Nijmegen, The Netherlands.,Radboud Center for Infectious Diseases, Center for Molecular and Biomolecular Informatics, Radboudumc, Nijmegen, The Netherlands
| |
Collapse
|
27
|
Pensar J, Puranen S, Arnold B, MacAlasdair N, Kuronen J, Tonkin-Hill G, Pesonen M, Xu Y, Sipola A, Sánchez-Busó L, Lees JA, Chewapreecha C, Bentley SD, Harris SR, Parkhill J, Croucher NJ, Corander J. Genome-wide epistasis and co-selection study using mutual information. Nucleic Acids Res 2019; 47:e112. [PMID: 31361894 PMCID: PMC6765119 DOI: 10.1093/nar/gkz656] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/04/2019] [Revised: 07/09/2019] [Accepted: 07/19/2019] [Indexed: 01/19/2023] Open
Abstract
Covariance-based discovery of polymorphisms under co-selective pressure or epistasis has received considerable recent attention in population genomics. Both statistical modeling of the population level covariation of alleles across the chromosome and model-free testing of dependencies between pairs of polymorphisms have been shown to successfully uncover patterns of selection in bacterial populations. Here we introduce a model-free method, SpydrPick, whose computational efficiency enables analysis at the scale of pan-genomes of many bacteria. SpydrPick incorporates an efficient correction for population structure, which adjusts for the phylogenetic signal in the data without requiring an explicit phylogenetic tree. We also introduce a new type of visualization of the results similar to the Manhattan plots used in genome-wide association studies, which enables rapid exploration of the identified signals of co-evolution. Simulations demonstrate the usefulness of our method and give some insight to when this type of analysis is most likely to be successful. Application of the method to large population genomic datasets of two major human pathogens, Streptococcus pneumoniae and Neisseria meningitidis, revealed both previously identified and novel putative targets of co-selection related to virulence and antibiotic resistance, highlighting the potential of this approach to drive molecular discoveries, even in the absence of phenotypic data.
Collapse
Affiliation(s)
- Johan Pensar
- Department of Mathematics and Statistics, Helsinki Institute for Information Technology (HIIT), Faculty of Science, University of Helsinki, FI-00014 Helsinki, Finland
| | - Santeri Puranen
- Department of Mathematics and Statistics, Helsinki Institute for Information Technology (HIIT), Faculty of Science, University of Helsinki, FI-00014 Helsinki, Finland.,Department of Computer Science, Aalto University, Espoo, FI-00014, Finland
| | - Brian Arnold
- Division of Informatics, Faculty of Arts and Sciences, Harvard University, Cambridge, MA 02138, USA
| | - Neil MacAlasdair
- Parasites and Microbes, Wellcome Sanger Institute, Cambridge, CB10 1SA, UK
| | - Juri Kuronen
- Department of Biostatistics, University of Oslo, Oslo, 0317, Norway
| | - Gerry Tonkin-Hill
- Parasites and Microbes, Wellcome Sanger Institute, Cambridge, CB10 1SA, UK
| | - Maiju Pesonen
- Department of Mathematics and Statistics, Helsinki Institute for Information Technology (HIIT), Faculty of Science, University of Helsinki, FI-00014 Helsinki, Finland.,Department of Computer Science, Aalto University, Espoo, FI-00014, Finland
| | - Yingying Xu
- Department of Mathematics and Statistics, Helsinki Institute for Information Technology (HIIT), Faculty of Science, University of Helsinki, FI-00014 Helsinki, Finland.,Department of Computer Science, Aalto University, Espoo, FI-00014, Finland
| | - Aleksi Sipola
- Department of Mathematics and Statistics, Helsinki Institute for Information Technology (HIIT), Faculty of Science, University of Helsinki, FI-00014 Helsinki, Finland
| | | | - John A Lees
- Department of Microbiology, New York University School of Medicine, New York, NY 10016, USA
| | - Claire Chewapreecha
- Department of Medicine, University of Cambridge, Cambridge CB2 0QQ, UK.,Bioinformatics & Systems Biology program, King Mongkut's University of Technology Thonburi, Bangkok 10150, Thailand
| | - Stephen D Bentley
- Parasites and Microbes, Wellcome Sanger Institute, Cambridge, CB10 1SA, UK
| | - Simon R Harris
- Parasites and Microbes, Wellcome Sanger Institute, Cambridge, CB10 1SA, UK
| | - Julian Parkhill
- Department of Veterinary Medicine, University of Cambridge, Madingley Road, Cambridge, CB3 0ES, UK
| | - Nicholas J Croucher
- MRC Centre for Global Infectious Disease Analysis, Department of Infectious Disease Epidemiology, St. Mary's Campus, Imperial College London, London, W2 1PG, UK
| | - Jukka Corander
- Department of Mathematics and Statistics, Helsinki Institute for Information Technology (HIIT), Faculty of Science, University of Helsinki, FI-00014 Helsinki, Finland.,Parasites and Microbes, Wellcome Sanger Institute, Cambridge, CB10 1SA, UK.,Department of Biostatistics, University of Oslo, Oslo, 0317, Norway
| |
Collapse
|
28
|
Abstract
Penicillin-non-susceptible Streptococcus pneumoniae (PNSP) were first detected in the 1960s, and are now common worldwide, predominantly through the international spread of a limited number of strains. Extant PNSP are characterized by mosaic pbp2x, pbp2b and pbp1a genes generated by interspecies recombinations, with the extent of these alterations determining the range and concentrations of β-lactams to which the genotype is non-susceptible. The complexity of the genetics underlying these phenotypes has been the subject of both molecular microbiology and genome-wide association and epistasis analyses. Such studies can aid our understanding of PNSP evolution and help improve the already highly-performing bioinformatic methods capable of identifying PNSP from genomic surveillance data.
Collapse
Affiliation(s)
- Tamsin C M Dewé
- MRC Centre for Global Infectious Disease Analysis, Department of Infectious Disease Epidemiology, St. Mary's Campus, Imperial College London, London, W2 1PG, UK
| | - Joshua C D'Aeth
- MRC Centre for Global Infectious Disease Analysis, Department of Infectious Disease Epidemiology, St. Mary's Campus, Imperial College London, London, W2 1PG, UK
| | - Nicholas J Croucher
- MRC Centre for Global Infectious Disease Analysis, Department of Infectious Disease Epidemiology, St. Mary's Campus, Imperial College London, London, W2 1PG, UK
| |
Collapse
|
29
|
Kwun MJ, Oggioni MR, Bentley SD, Fraser C, Croucher NJ. Synergistic Activity of Mobile Genetic Element Defences in Streptococcus pneumoniae. Genes (Basel) 2019; 10:genes10090707. [PMID: 31540216 PMCID: PMC6771155 DOI: 10.3390/genes10090707] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2019] [Revised: 09/01/2019] [Accepted: 09/06/2019] [Indexed: 01/02/2023] Open
Abstract
A diverse set of mobile genetic elements (MGEs) transmit between Streptococcus pneumoniae cells, but many isolates remain uninfected. The best-characterised defences against horizontal transmission of MGEs are restriction-modification systems (RMSs), of which there are two phase-variable examples in S. pneumoniae. Additionally, the transformation machinery has been proposed to limit vertical transmission of chromosomally integrated MGEs. This work describes how these mechanisms can act in concert. Experimental data demonstrate RMS phase variation occurs at a sub-maximal rate. Simulations suggest this may be optimal if MGEs are sometimes vertically inherited, as it reduces the probability that an infected cell will switch between RMS variants while the MGE is invading the population, and thereby undermine the restriction barrier. Such vertically inherited MGEs can be deleted by transformation. The lack of between-strain transformation hotspots at known prophage att sites suggests transformation cannot remove an MGE from a strain in which it is fixed. However, simulations confirmed that transformation was nevertheless effective at preventing the spread of MGEs into a previously uninfected cell population, if a recombination barrier existed between co-colonising strains. Further simulations combining these effects of phase variable RMSs and transformation found they synergistically inhibited MGEs spreading, through limiting both vertical and horizontal transmission.
Collapse
Affiliation(s)
- Min Jung Kwun
- MRC Centre for Global Infectious Disease Analysis, Department of Infectious Disease Epidemiology, St. Mary's Campus, Imperial College London, London W2 1PG, UK.
| | - Marco R Oggioni
- Department of Genetics and Genome Biology, University of Leicester, Leicester LE1 7RH, UK.
| | - Stephen D Bentley
- Pathogens and Microbes, Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton, Cambridge CB10 1SA, UK.
| | - Christophe Fraser
- Big Data Institute, Nuffield Department of Medicine, Old Road Campus, University of Oxford, Oxford OX3 7LF, UK.
| | - Nicholas J Croucher
- MRC Centre for Global Infectious Disease Analysis, Department of Infectious Disease Epidemiology, St. Mary's Campus, Imperial College London, London W2 1PG, UK.
| |
Collapse
|
30
|
Nasher F, Kwun MJ, Croucher NJ, Heller M, Hathaway LJ. Peptide Occurring in Enterobacteriaceae Triggers Streptococcus pneumoniae Cell Death. Front Cell Infect Microbiol 2019; 9:320. [PMID: 31552200 PMCID: PMC6748166 DOI: 10.3389/fcimb.2019.00320] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2019] [Accepted: 08/26/2019] [Indexed: 11/15/2022] Open
Abstract
Non-encapsulated Streptococcus pneumoniae often possess two genes, aliB-like ORF 1 and aliB-like ORF 2, in place of capsule genes. AliB-like ORF 1 is thought to encode a substrate binding protein of an ABC transporter which binds peptide SETTFGRDFN, found in 50S ribosomal subunit protein L4 of Enterobacteriaceae. Here, we investigated the effect of binding of AliB-like ORF 1 peptide on the transcriptome and proteome of non-encapsulated pneumococci. We found upregulation of gene expression of a metacaspase and a gene encoding N-acetylmuramoyl-L-alanine amidase, both of which are proposed to be involved in programmed cell death in prokaryotic cells. Proteome profiling indicated upregulation of transcriptional regulators and downregulation of metabolism-associated genes. Exposure to the peptide specifically triggered death in pneumococci which express AliB-like ORF 1, with the bacteria having an apoptotic appearance by electron microscopy. We propose that binding of the AliB-like ORF 1 peptide ligand by the pneumococcus signals a challenging environment with hostile bacterial species leading to death of a proportion of the pneumococcal population.
Collapse
Affiliation(s)
- Fauzy Nasher
- Institute for Infectious Diseases, Faculty of Medicine, University of Bern, Bern, Switzerland.,Graduate School for Cellular and Biomedical Sciences, University of Bern, Bern, Switzerland
| | - Min Jung Kwun
- Department of Infectious Disease Epidemiology, MRC Centre for Global Infectious Disease Analysis, Imperial College London, London, United Kingdom
| | - Nicholas J Croucher
- Department of Infectious Disease Epidemiology, MRC Centre for Global Infectious Disease Analysis, Imperial College London, London, United Kingdom
| | - Manfred Heller
- Proteomics and Mass Spectrometry Core Facility, Department for BioMedical Research, University of Bern, Bern, Switzerland
| | - Lucy J Hathaway
- Institute for Infectious Diseases, Faculty of Medicine, University of Bern, Bern, Switzerland
| |
Collapse
|
31
|
Kwun MJ, Oggioni MR, De Ste Croix M, Bentley SD, Croucher NJ. Excision-reintegration at a pneumococcal phase-variable restriction-modification locus drives within- and between-strain epigenetic differentiation and inhibits gene acquisition. Nucleic Acids Res 2019; 46:11438-11453. [PMID: 30321375 PMCID: PMC6265443 DOI: 10.1093/nar/gky906] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2018] [Accepted: 10/08/2018] [Indexed: 12/21/2022] Open
Abstract
Phase-variation of Type I restriction-modification systems can rapidly alter the sequence motifs they target, diversifying both the epigenetic patterns and endonuclease activity within clonally descended populations. Here, we characterize the Streptococcus pneumoniae SpnIV phase-variable Type I RMS, encoded by the translocating variable restriction (tvr) locus, to identify its target motifs, mechanism and regulation of phase variation, and effects on exchange of sequence through transformation. The specificity-determining hsdS genes were shuffled through a recombinase-mediated excision-reintegration mechanism involving circular intermediate molecules, guided by two types of direct repeat. The rate of rearrangements was limited by an attenuator and toxin-antitoxin system homologs that inhibited recombinase gene transcription. Target motifs for both the SpnIV, and multiple Type II, MTases were identified through methylation-sensitive sequencing of a panel of recombinase-null mutants. This demonstrated the species-wide diversity observed at the tvr locus can likely specify nine different methylation patterns. This will reduce sequence exchange in this diverse species, as the native form of the SpnIV RMS was demonstrated to inhibit the acquisition of genomic islands by transformation. Hence the tvr locus can drive variation in genome methylation both within and between strains, and limits the genomic plasticity of S. pneumoniae.
Collapse
Affiliation(s)
- Min Jung Kwun
- MRC Centre for Global Infectious Disease Analysis, Department of Infectious Disease Epidemiology, Imperial College London, London W2 1PG, UK
| | - Marco R Oggioni
- Department of Genetics, University of Leicester, Leicester LE1 7RH, UK
| | | | - Stephen D Bentley
- Parasites and Microbes, Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton, Cambridge CB10 1SA, UK
| | - Nicholas J Croucher
- MRC Centre for Global Infectious Disease Analysis, Department of Infectious Disease Epidemiology, Imperial College London, London W2 1PG, UK
| |
Collapse
|
32
|
Didelot X, Croucher NJ, Bentley SD, Harris SR, Wilson DJ. Bayesian inference of ancestral dates on bacterial phylogenetic trees. Nucleic Acids Res 2019; 46:e134. [PMID: 30184106 PMCID: PMC6294524 DOI: 10.1093/nar/gky783] [Citation(s) in RCA: 116] [Impact Index Per Article: 23.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2018] [Accepted: 08/21/2018] [Indexed: 12/15/2022] Open
Abstract
The sequencing and comparative analysis of a collection of bacterial genomes from a single species or lineage of interest can lead to key insights into its evolution, ecology or epidemiology. The tool of choice for such a study is often to build a phylogenetic tree, and more specifically when possible a dated phylogeny, in which the dates of all common ancestors are estimated. Here, we propose a new Bayesian methodology to construct dated phylogenies which is specifically designed for bacterial genomics. Unlike previous Bayesian methods aimed at building dated phylogenies, we consider that the phylogenetic relationships between the genomes have been previously evaluated using a standard phylogenetic method, which makes our methodology much faster and scalable. This two-step approach also allows us to directly exploit existing phylogenetic methods that detect bacterial recombination, and therefore to account for the effect of recombination in the construction of a dated phylogeny. We analysed many simulated datasets in order to benchmark the performance of our approach in a wide range of situations. Furthermore, we present applications to three different real datasets from recent bacterial genomic studies. Our methodology is implemented in a R package called BactDating which is freely available for download at https://github.com/xavierdidelot/BactDating.
Collapse
Affiliation(s)
- Xavier Didelot
- Department of Infectious Disease Epidemiology, School of Public Health, Imperial College London, London, UK
| | - Nicholas J Croucher
- Department of Infectious Disease Epidemiology, School of Public Health, Imperial College London, London, UK
| | - Stephen D Bentley
- The Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton, Cambridge, UK
| | - Simon R Harris
- The Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton, Cambridge, UK
| | - Daniel J Wilson
- Big Data Institute, Nuffield Department of Population Health, University of Oxford, Oxford, UK
| |
Collapse
|
33
|
Lees JA, Ferwerda B, Kremer PHC, Wheeler NE, Serón MV, Croucher NJ, Gladstone RA, Bootsma HJ, Rots NY, Wijmega-Monsuur AJ, Sanders EAM, Trzciński K, Wyllie AL, Zwinderman AH, van den Berg LH, van Rheenen W, Veldink JH, Harboe ZB, Lundbo LF, de Groot LCPGM, van Schoor NM, van der Velde N, Ängquist LH, Sørensen TIA, Nohr EA, Mentzer AJ, Mills TC, Knight JC, du Plessis M, Nzenze S, Weiser JN, Parkhill J, Madhi S, Benfield T, von Gottberg A, van der Ende A, Brouwer MC, Barrett JC, Bentley SD, van de Beek D. Joint sequencing of human and pathogen genomes reveals the genetics of pneumococcal meningitis. Nat Commun 2019; 10:2176. [PMID: 31092817 PMCID: PMC6520353 DOI: 10.1038/s41467-019-09976-3] [Citation(s) in RCA: 58] [Impact Index Per Article: 11.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2018] [Accepted: 04/11/2019] [Indexed: 12/21/2022] Open
Abstract
Streptococcus pneumoniae is a common nasopharyngeal colonizer, but can also cause life-threatening invasive diseases such as empyema, bacteremia and meningitis. Genetic variation of host and pathogen is known to play a role in invasive pneumococcal disease, though to what extent is unknown. In a genome-wide association study of human and pathogen we show that human variation explains almost half of variation in susceptibility to pneumococcal meningitis and one-third of variation in severity, identifying variants in CCDC33 associated with susceptibility. Pneumococcal genetic variation explains a large amount of invasive potential (70%), but has no effect on severity. Serotype alone is insufficient to explain invasiveness, suggesting other pneumococcal factors are involved in progression to invasive disease. We identify pneumococcal genes involved in invasiveness including pspC and zmpD, and perform a human-bacteria interaction analysis. These genes are potential candidates for the development of more broadly-acting pneumococcal vaccines.
Collapse
Affiliation(s)
- John A Lees
- Department of Microbiology, New York University School of Medicine, New York, NY, 10016, USA
- Parasites and Microbes, Wellcome Sanger Institute, Cambridge, CB10 1SA, UK
| | - Bart Ferwerda
- Amsterdam UMC, University of Amsterdam, Department of Neurology, Amsterdam Neuroscience, Meibergdreef 9, Amsterdam, 1105 AZ, The Netherlands
| | - Philip H C Kremer
- Amsterdam UMC, University of Amsterdam, Department of Neurology, Amsterdam Neuroscience, Meibergdreef 9, Amsterdam, 1105 AZ, The Netherlands
| | - Nicole E Wheeler
- Parasites and Microbes, Wellcome Sanger Institute, Cambridge, CB10 1SA, UK
- The Centre for Genomic Pathogen Surveillance, Wellcome Sanger Institute, Cambridge, CB10 1SA, UK
| | - Mercedes Valls Serón
- Amsterdam UMC, University of Amsterdam, Department of Neurology, Amsterdam Neuroscience, Meibergdreef 9, Amsterdam, 1105 AZ, The Netherlands
| | - Nicholas J Croucher
- MRC Centre for Global Infectious Disease Analysis, Department of Infectious Disease Epidemiology, Imperial College London, London, W2 1PG, UK
| | | | - Hester J Bootsma
- Centre for Infectious Disease Control, National Institute for Public Health and the Environment, Bilthoven, 3721 MA, The Netherlands
| | - Nynke Y Rots
- Centre for Infectious Disease Control, National Institute for Public Health and the Environment, Bilthoven, 3721 MA, The Netherlands
| | - Alienke J Wijmega-Monsuur
- Centre for Infectious Disease Control, National Institute for Public Health and the Environment, Bilthoven, 3721 MA, The Netherlands
| | - Elisabeth A M Sanders
- Centre for Infectious Disease Control, National Institute for Public Health and the Environment, Bilthoven, 3721 MA, The Netherlands
- Department of Pediatric Immunology and Infectious Diseases, Wilhelmina Children's Hospital, University Medical Centre Utrecht, Utrecht, 3508 AB, The Netherlands
| | - Krzysztof Trzciński
- Department of Pediatric Immunology and Infectious Diseases, Wilhelmina Children's Hospital, University Medical Centre Utrecht, Utrecht, 3508 AB, The Netherlands
| | - Anne L Wyllie
- Department of Pediatric Immunology and Infectious Diseases, Wilhelmina Children's Hospital, University Medical Centre Utrecht, Utrecht, 3508 AB, The Netherlands
- Epidemiology of Microbial Diseases, Yale School of Public Health, New Haven, CT, 06520, USA
| | - Aeilko H Zwinderman
- Amsterdam UMC, University of Amsterdam, Department of Clinical Epidemiology, Biostatistics and Bioinformatics, Amsterdam Public Health, Meibergdreef 9, Amsterdam, 1105 AZ, The Netherlands
| | - Leonard H van den Berg
- Department of Neurology, Brain Center Rudolf Magnus, University Medical Center Utrecht, Utrecht, 3584 CG, The Netherlands
| | - Wouter van Rheenen
- Department of Neurology, Brain Center Rudolf Magnus, University Medical Center Utrecht, Utrecht, 3584 CG, The Netherlands
| | - Jan H Veldink
- Department of Neurology, Brain Center Rudolf Magnus, University Medical Center Utrecht, Utrecht, 3584 CG, The Netherlands
| | - Zitta B Harboe
- Department of Microbiological Surveillance and Research, Statens Serum Institut, Copenhagen, DK-2300, Denmark
| | - Lene F Lundbo
- Department of Infectious Diseases, Hvidovre Hospital, University of Copenhagen, Hvidovre, 2650, Denmark
| | - Lisette C P G M de Groot
- Department of Human Nutrition, Wageningen University, P.O. Box 17, 6700 AA, Wageningen, The Netherlands
| | - Natasja M van Schoor
- Amsterdam UMC, VU University, Department of Epidemiology and Biostatistics, Amsterdam Public Health, Van der Boechorststraat 7, Amsterdam, 1007 MB, The Netherlands
| | - Nathalie van der Velde
- Amsterdam UMC, University of Amsterdam, Department of Internal Medicine, Geriatrics, Amsterdam Public Health, Meibergdreef 9, Amsterdam, 1105 AZ, The Netherlands
- Department of Internal Medicine, Erasmus MC, University Medical Centre Rotterdam, P.O. Box 2040, 3000 CA, Rotterdam, The Netherlands
| | - Lars H Ängquist
- Center for Clinical Research and Disease Prevention, Bispebjerg and Frederiksberg Hospitals, The Capital Region, Copenhagen, DK-2000, Denmark
| | - Thorkild I A Sørensen
- The Novo Nordisk Foundation Center for Basic Metabolic Research, Section of Metabolic Genetics, Copenhagen, DK-2200, Denmark
- The Department of Public Health, Section of Epidemiology, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, DK-1014, Denmark
| | - Ellen A Nohr
- Institute of Clinical Research, University of Southern Denmark, Odense, DK-5000, Denmark
| | - Alexander J Mentzer
- Wellcome Centre for Human Genetics, Nuffield Department of Medicine, University of Oxford, Oxford, OX3 7BN, UK
| | - Tara C Mills
- Wellcome Centre for Human Genetics, Nuffield Department of Medicine, University of Oxford, Oxford, OX3 7BN, UK
| | - Julian C Knight
- Wellcome Centre for Human Genetics, Nuffield Department of Medicine, University of Oxford, Oxford, OX3 7BN, UK
| | - Mignon du Plessis
- School of Pathology, Faculty of Health Sciences, University of Witwatersrand, Johannesburg, 2000, South Africa
| | - Susan Nzenze
- School of Pathology, Faculty of Health Sciences, University of Witwatersrand, Johannesburg, 2000, South Africa
| | - Jeffrey N Weiser
- Department of Microbiology, New York University School of Medicine, New York, NY, 10016, USA
| | - Julian Parkhill
- Parasites and Microbes, Wellcome Sanger Institute, Cambridge, CB10 1SA, UK
| | - Shabir Madhi
- National Institute for Communicable Diseases, Johannesburg, 2192, South Africa
| | - Thomas Benfield
- Department of Infectious Diseases, Hvidovre Hospital, University of Copenhagen, Hvidovre, 2650, Denmark
| | - Anne von Gottberg
- School of Pathology, Faculty of Health Sciences, University of Witwatersrand, Johannesburg, 2000, South Africa
- National Institute for Communicable Diseases, Johannesburg, 2192, South Africa
| | - Arie van der Ende
- Amsterdam UMC, University of Amsterdam, Department of Medical Microbiology, Amsterdam Infection and Immunity, Meibergdreef 9, Amsterdam, 1105 AZ, The Netherlands
- Netherlands Reference Laboratory for Bacterial Meningitis, Amsterdam UMC/RIVM, University of Amsterdam, Meibergdreef 9, Amsterdam, 1105 AZ, The Netherlands
| | - Matthijs C Brouwer
- Amsterdam UMC, University of Amsterdam, Department of Neurology, Amsterdam Neuroscience, Meibergdreef 9, Amsterdam, 1105 AZ, The Netherlands
| | - Jeffrey C Barrett
- Parasites and Microbes, Wellcome Sanger Institute, Cambridge, CB10 1SA, UK
- Genomics Plc, East Road, Cambridge, CB1 1BH, UK
| | - Stephen D Bentley
- Parasites and Microbes, Wellcome Sanger Institute, Cambridge, CB10 1SA, UK.
| | - Diederik van de Beek
- Amsterdam UMC, University of Amsterdam, Department of Neurology, Amsterdam Neuroscience, Meibergdreef 9, Amsterdam, 1105 AZ, The Netherlands.
| |
Collapse
|
34
|
Gladstone RA, Lo SW, Lees JA, Croucher NJ, van Tonder AJ, Corander J, Page AJ, Marttinen P, Bentley LJ, Ochoa TJ, Ho PL, du Plessis M, Cornick JE, Kwambana-Adams B, Benisty R, Nzenze SA, Madhi SA, Hawkins PA, Everett DB, Antonio M, Dagan R, Klugman KP, von Gottberg A, McGee L, Breiman RF, Bentley SD. International genomic definition of pneumococcal lineages, to contextualise disease, antibiotic resistance and vaccine impact. EBioMedicine 2019; 43:338-346. [PMID: 31003929 PMCID: PMC6557916 DOI: 10.1016/j.ebiom.2019.04.021] [Citation(s) in RCA: 121] [Impact Index Per Article: 24.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2019] [Revised: 04/05/2019] [Accepted: 04/09/2019] [Indexed: 11/25/2022] Open
Abstract
BACKGROUND Pneumococcal conjugate vaccines have reduced the incidence of invasive pneumococcal disease, caused by vaccine serotypes, but non-vaccine-serotypes remain a concern. We used whole genome sequencing to study pneumococcal serotype, antibiotic resistance and invasiveness, in the context of genetic background. METHODS Our dataset of 13,454 genomes, combined with four published genomic datasets, represented Africa (40%), Asia (25%), Europe (19%), North America (12%), and South America (5%). These 20,027 pneumococcal genomes were clustered into lineages using PopPUNK, and named Global Pneumococcal Sequence Clusters (GPSCs). From our dataset, we additionally derived serotype and sequence type, and predicted antibiotic sensitivity. We then measured invasiveness using odds ratios that relating prevalence in invasive pneumococcal disease to carriage. FINDINGS The combined collections (n = 20,027) were clustered into 621 GPSCs. Thirty-five GPSCs observed in our dataset were represented by >100 isolates, and subsequently classed as dominant-GPSCs. In 22/35 (63%) of dominant-GPSCs both non-vaccine serotypes and vaccine serotypes were observed in the years up until, and including, the first year of pneumococcal conjugate vaccine introduction. Penicillin and multidrug resistance were higher (p < .05) in a subset dominant-GPSCs (14/35, 9/35 respectively), and resistance to an increasing number of antibiotic classes was associated with increased recombination (R2 = 0.27 p < .0001). In 28/35 dominant-GPSCs, the country of isolation was a significant predictor (p < .05) of its antibiogram (mean misclassification error 0.28, SD ± 0.13). We detected increased invasiveness of six genetic backgrounds, when compared to other genetic backgrounds expressing the same serotype. Up to 1.6-fold changes in invasiveness odds ratio were observed. INTERPRETATION We define GPSCs that can be assigned to any pneumococcal genomic dataset, to aid international comparisons. Existing non-vaccine-serotypes in most GPSCs preclude the removal of these lineages by pneumococcal conjugate vaccines; leaving potential for serotype replacement. A subset of GPSCs have increased resistance, and/or serotype-independent invasiveness.
Collapse
Affiliation(s)
| | - Stephanie W Lo
- Parasites and microbes, Wellcome Sanger Institute, Hinxton, UK
| | - John A Lees
- New York University School of Medicine, New York, NY, USA
| | | | | | - Jukka Corander
- Parasites and microbes, Wellcome Sanger Institute, Hinxton, UK; Department of Biostatistics, University of Oslo, 0317 Oslo, Norway
| | - Andrew J Page
- Parasites and microbes, Wellcome Sanger Institute, Hinxton, UK
| | - Pekka Marttinen
- Department of Computer Science, Helsinki Institute for Information Technology HIIT, Espoo, Finland
| | - Leon J Bentley
- Parasites and microbes, Wellcome Sanger Institute, Hinxton, UK
| | - Theresa J Ochoa
- Instituto de Medicina Tropical, Universidad Peruana Cayetano Heredia, Lima, Peru
| | - Pak Leung Ho
- Department of Microbiology, Carol Yu Centre for Infection, Queen Mary Hospital, The University of Hong Kong, Hong Kong, China
| | - Mignon du Plessis
- Centre for Respiratory Diseases and Meningitis, National Institute for Communicable Diseases, Johannesburg, South Africa
| | - Jennifer E Cornick
- Malawi-Liverpool-Wellcome-Trust Clinical Research Programme, Blantyre, Malawi
| | - Brenda Kwambana-Adams
- NIHR Global Health Research Unit on Mucosal Pathogens, Division of Infection and Immunity, University College London, London, UK; WHO Collaborating Centre for New Vaccines Surveillance, Medical Research Council Unit The Gambia at London School of Hygiene & Tropical Medicine, Atlantic Boulevard, Fajara, PO Box 273 Banjul, the Gambia
| | - Rachel Benisty
- The Faculty of Health Sciences, Ben-Gurion University of the Negev, Beer-Sheva, Israel
| | - Susan A Nzenze
- Medical Research Council: Respiratory and Meningeal Pathogens Research Unit, University of the Witwatersrand, South Africa; Department of Science and Technology, National Research Foundation: Vaccine Preventable Diseases, University of the Witwatersrand, South Africa
| | - Shabir A Madhi
- Medical Research Council: Respiratory and Meningeal Pathogens Research Unit, University of the Witwatersrand, South Africa; Department of Science and Technology, National Research Foundation: Vaccine Preventable Diseases, University of the Witwatersrand, South Africa
| | | | | | - Martin Antonio
- WHO Collaborating Centre for New Vaccines Surveillance, Medical Research Council Unit The Gambia at London School of Hygiene & Tropical Medicine, Atlantic Boulevard, Fajara, PO Box 273 Banjul, the Gambia; Division of Microbiology & Immunity, Warwick Medical School, University of Warwick, Coventry, CV4 7AL, UK
| | - Ron Dagan
- The Faculty of Health Sciences, Ben-Gurion University of the Negev, Beer-Sheva, Israel
| | | | - Anne von Gottberg
- Centre for Respiratory Diseases and Meningitis, National Institute for Communicable Diseases, Johannesburg, South Africa
| | - Lesley McGee
- Centers for Disease Control and Prevention, Atlanta, USA
| | - Robert F Breiman
- Rollins School Public Health, Emory University, USA; Emory Global Health Institute, Atlanta, USA
| | | |
Collapse
|
35
|
McNally A, Kallonen T, Connor C, Abudahab K, Aanensen DM, Horner C, Peacock SJ, Parkhill J, Croucher NJ, Corander J. Diversification of Colonization Factors in a Multidrug-Resistant Escherichia coli Lineage Evolving under Negative Frequency-Dependent Selection. mBio 2019; 10:e00644-19. [PMID: 31015329 PMCID: PMC6479005 DOI: 10.1128/mbio.00644-19] [Citation(s) in RCA: 71] [Impact Index Per Article: 14.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2019] [Accepted: 03/21/2019] [Indexed: 12/18/2022] Open
Abstract
Escherichia coli is a major cause of bloodstream and urinary tract infections globally. The wide dissemination of multidrug-resistant (MDR) strains of extraintestinal pathogenic E. coli (ExPEC) poses a rapidly increasing public health burden due to narrowed treatment options and increased risk of failure to clear an infection. Here, we present a detailed population genomic analysis of the ExPEC ST131 clone, in which we seek explanations for its success as an emerging pathogenic strain beyond the acquisition of antimicrobial resistance (AMR) genes. We show evidence for evolution toward separate ecological niches for the main clades of ST131 and differential evolution of anaerobic metabolism, key colonization, and virulence factors. We further demonstrate that negative frequency-dependent selection acting across accessory loci is a major mechanism that has shaped the population evolution of this pathogen.IMPORTANCE Infections with multidrug-resistant (MDR) strains of Escherichia coli are a significant global public health concern. To combat these pathogens, we need a deeper understanding of how they evolved from their background populations. By understanding the processes that underpin their emergence, we can design new strategies to limit evolution of new clones and combat existing clones. By combining population genomics with modelling approaches, we show that dominant MDR clones of E. coli are under the influence of negative frequency-dependent selection, preventing them from rising to fixation in a population. Furthermore, we show that this selection acts on genes involved in anaerobic metabolism, suggesting that this key trait, and the ability to colonize human intestinal tracts, is a key step in the evolution of MDR clones of E. coli.
Collapse
Affiliation(s)
- Alan McNally
- Institute of Microbiology and Infection, University of Birmingham, Birmingham, United Kingdom
| | - Teemu Kallonen
- Infection Genomics, Wellcome Sanger Institute, Cambridge, United Kingdom
- Department of Biostatistics, University of Oslo, Oslo, Norway
| | - Christopher Connor
- Institute of Microbiology and Infection, University of Birmingham, Birmingham, United Kingdom
| | - Khalil Abudahab
- Infection Genomics, Wellcome Sanger Institute, Cambridge, United Kingdom
| | - David M Aanensen
- Infection Genomics, Wellcome Sanger Institute, Cambridge, United Kingdom
| | - Carolyne Horner
- British Society of Antimicrobial Chemotherapy, Birmingham, United Kingdom
| | - Sharon J Peacock
- Infection Genomics, Wellcome Sanger Institute, Cambridge, United Kingdom
- Department of Medicine, University of Cambridge, Cambridge, United Kingdom
- London School of Hygiene and Tropical Medicine, London, United Kingdom
| | - Julian Parkhill
- Infection Genomics, Wellcome Sanger Institute, Cambridge, United Kingdom
| | - Nicholas J Croucher
- Faculty of Medicine, School of Public Health, Imperial College, London, United Kingdom
| | - Jukka Corander
- Infection Genomics, Wellcome Sanger Institute, Cambridge, United Kingdom
- Department of Biostatistics, University of Oslo, Oslo, Norway
- Department of Mathematics and Statistics, University of Helsinki, Helsinki, Finland
| |
Collapse
|
36
|
Mitchell PK, Azarian T, Croucher NJ, Callendrello A, Thompson CM, Pelton SI, Lipsitch M, Hanage WP. Population genomics of pneumococcal carriage in Massachusetts children following introduction of PCV-13. Microb Genom 2019; 5. [PMID: 30777813 PMCID: PMC6421351 DOI: 10.1099/mgen.0.000252] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022] Open
Abstract
The 13-valent pneumococcal conjugate vaccine (PCV-13) was introduced in the United States in 2010. Using a large paediatric carriage sample collected from shortly after the introduction of PCV-7 to several years after the introduction of PCV-13, we investigate alterations in the composition of the pneumococcal population following the introduction of PCV-13, evaluating the extent to which the post-vaccination non-vaccine type (NVT) population mirrors that from prior to vaccine introduction and the effect of PCV-13 on vaccine type lineages. Draft genome assemblies from 736 newly sequenced and 616 previously published pneumococcal carriage isolates from children in Massachusetts between 2001 and 2014 were analysed. Isolates were classified into one of 22 sequence clusters (SCs) on the basis of their core genome sequence. We calculated the SC diversity for each sampling period as the probability that any two randomly drawn isolates from that period belong to different SCs. The sampling period immediately after the introduction of PCV-13 (2011) was found to have higher diversity than preceding (2007) or subsequent (2014) sampling periods {Simpson’s D 2007: 0.915 [95 % confidence interval (CI) 0.901, 0.929]; 2011: 0.935 [0.927, 0.942]; 2014 : 0.912 [0.901, 0.923]}. Amongst NVT isolates, we found the distribution of SCs in 2011 to be significantly different from that in 2007 or 2014 (Fisher’s exact test P=0.018, 0.0078), but did not find a difference comparing 2007 to 2014 (Fisher’s exact test P=0.24), indicating greater similarity between samples separated by a longer time period than between samples from closer time periods. We also found changes in the accessory gene content of the NVT population between 2007 and 2011 to have been reduced by 2014. Amongst the new serotypes targeted by PCV-13, four were present in our sample. The proportion of our sample composed of PCV-13-only vaccine serotypes 19A, 6C and 7F decreased between 2007 and 2014, but no such reduction was seen for serotype 3. We did, however, observe differences in the genetic composition of the pre- and post-PCV-13 serotype 3 population. Our isolates were collected during discrete sampling periods from a small geographical area, which may limit the generalizability of our findings. Pneumococcal diversity increased immediately following the introduction of PCV-13, but subsequently returned to pre-vaccination levels. This is reflected in the distribution of NVT lineages, and, to a lesser extent, their accessory gene frequencies. As such, there may be a period during which the population is particularly disrupted by vaccination before returning to a more stable distribution. The persistence and shifting genetic composition of serotype 3 is a concern and warrants further investigation.
Collapse
Affiliation(s)
- Patrick K Mitchell
- 1Center for Communicable Disease Dynamics, Department of Epidemiology, T.H. Chan School of Public Health, Harvard University, Boston, MA, USA
| | - Taj Azarian
- 1Center for Communicable Disease Dynamics, Department of Epidemiology, T.H. Chan School of Public Health, Harvard University, Boston, MA, USA
| | - Nicholas J Croucher
- 2MRC Centre for Global Infectious Disease Analysis, Department of Infectious Disease Epidemiology, Imperial College London, London W2 1PG, UK
| | - Alanna Callendrello
- 1Center for Communicable Disease Dynamics, Department of Epidemiology, T.H. Chan School of Public Health, Harvard University, Boston, MA, USA
| | - Claudette M Thompson
- 1Center for Communicable Disease Dynamics, Department of Epidemiology, T.H. Chan School of Public Health, Harvard University, Boston, MA, USA
| | - Stephen I Pelton
- 3Division of Pediatric Infectious Diseases, Maxwell Finland Laboratory for Infectious Diseases, Boston Medical Center, Boston, MA, USA
| | - Marc Lipsitch
- 1Center for Communicable Disease Dynamics, Department of Epidemiology, T.H. Chan School of Public Health, Harvard University, Boston, MA, USA
| | - William P Hanage
- 1Center for Communicable Disease Dynamics, Department of Epidemiology, T.H. Chan School of Public Health, Harvard University, Boston, MA, USA
| |
Collapse
|
37
|
Lees JA, Harris SR, Tonkin-Hill G, Gladstone RA, Lo SW, Weiser JN, Corander J, Bentley SD, Croucher NJ. Fast and flexible bacterial genomic epidemiology with PopPUNK. Genome Res 2019; 29:304-316. [PMID: 30679308 PMCID: PMC6360808 DOI: 10.1101/gr.241455.118] [Citation(s) in RCA: 178] [Impact Index Per Article: 35.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/05/2018] [Accepted: 12/10/2018] [Indexed: 12/02/2022]
Abstract
The routine use of genomics for disease surveillance provides the opportunity for high-resolution bacterial epidemiology. Current whole-genome clustering and multilocus typing approaches do not fully exploit core and accessory genomic variation, and they cannot both automatically identify, and subsequently expand, clusters of significantly similar isolates in large data sets spanning entire species. Here, we describe PopPUNK (Population Partitioning Using Nucleotide K-mers), a software implementing scalable and expandable annotation- and alignment-free methods for population analysis and clustering. Variable-length k-mer comparisons are used to distinguish isolates’ divergence in shared sequence and gene content, which we demonstrate to be accurate over multiple orders of magnitude using data from both simulations and genomic collections representing 10 taxonomically widespread species. Connections between closely related isolates of the same strain are robustly identified, despite interspecies variation in the pairwise distance distributions that reflects species’ diverse evolutionary patterns. PopPUNK can process 103–104 genomes in a single batch, with minimal memory use and runtimes up to 200-fold faster than existing model-based methods. Clusters of strains remain consistent as new batches of genomes are added, which is achieved without needing to reanalyze all genomes de novo. This facilitates real-time surveillance with consistent cluster naming between studies and allows for outbreak detection using hundreds of genomes in minutes. Interactive visualization and online publication is streamlined through the automatic output of results to multiple platforms. PopPUNK has been designed as a flexible platform that addresses important issues with currently used whole-genome clustering and typing methods, and has potential uses across bacterial genetics and public health research.
Collapse
Affiliation(s)
- John A Lees
- Department of Microbiology, New York University School of Medicine, New York, New York 10016, USA
| | - Simon R Harris
- Parasites and Microbes, Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton CB10 1SA, United Kingdom
| | - Gerry Tonkin-Hill
- Parasites and Microbes, Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton CB10 1SA, United Kingdom
| | - Rebecca A Gladstone
- Parasites and Microbes, Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton CB10 1SA, United Kingdom
| | - Stephanie W Lo
- Parasites and Microbes, Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton CB10 1SA, United Kingdom
| | - Jeffrey N Weiser
- Department of Microbiology, New York University School of Medicine, New York, New York 10016, USA
| | - Jukka Corander
- Parasites and Microbes, Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton CB10 1SA, United Kingdom.,Department of Biostatistics, University of Oslo, 0372 Oslo, Norway.,Helsinki Institute of Information Technology, Department of Mathematics and Statistics, University of Helsinki, 00014 Helsinki, Finland
| | - Stephen D Bentley
- Parasites and Microbes, Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton CB10 1SA, United Kingdom.,Institute of Infection and Global Health, University of Liverpool, Liverpool L7 3EA, United Kingdom.,Department of Pathology, University of Cambridge, Cambridge CB2 1QP, United Kingdom
| | - Nicholas J Croucher
- MRC Centre for Global Infectious Disease Analysis, Department of Infectious Disease Epidemiology, Imperial College London, London W2 1PG, United Kingdom
| |
Collapse
|
38
|
Campo JJ, Le TQ, Pablo JV, Hung C, Teng AA, Tettelin H, Tate A, Hanage WP, Alderson MR, Liang X, Malley R, Lipsitch M, Croucher NJ. Panproteome-wide analysis of antibody responses to whole cell pneumococcal vaccination. eLife 2018; 7:37015. [PMID: 30592459 PMCID: PMC6344088 DOI: 10.7554/elife.37015] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2018] [Accepted: 12/25/2018] [Indexed: 11/13/2022] Open
Abstract
Pneumococcal whole cell vaccines (WCVs) could cost-effectively protect against a greater strain diversity than current capsule-based vaccines. Immunoglobulin G (IgG) responses to a WCV were characterised by applying longitudinally-sampled sera, available from 35 adult placebo-controlled phase I trial participants, to a panproteome microarray. Despite individuals maintaining distinctive antibody ‘fingerprints’, responses were consistent across vaccinated cohorts. Seventy-two functionally distinct proteins were associated with WCV-induced increases in IgG binding. These shared characteristics with naturally immunogenic proteins, being enriched for transporters and cell wall metabolism enzymes, likely unusually exposed on the unencapsulated WCV’s surface. Vaccine-induced responses were specific to variants of the diverse PclA, PspC and ZmpB proteins, whereas PspA- and ZmpA-induced antibodies recognised a broader set of alleles. Temporal variation in IgG levels suggested a mixture of anamnestic and novel responses. These reproducible increases in IgG binding to a limited, but functionally diverse, set of conserved proteins indicate WCV could provide species-wide immunity. Clinical trial registration: The trial was registered with ClinicalTrials.gov with Identifier NCT01537185; the results are available from https://clinicaltrials.gov/ct2/show/results/NCT01537185. Streptococcus pneumoniae is a bug that causes pneumonia and meningitis, killing around a million people each year. Vaccines now exist to protect young children against these diseases, but they are expensive and do not work against all the strains of the bacteria. This is because these shots train the body’s immune system to recognize and attack the bacterium’s capsule, a layer of sugars that surrounds the microbe and is often different between strains. One possible solution could be a cheap, whole cell vaccine. These injections expose the body to genetically modified S. pneumoniae that do not carry the capsule. Such treatment has now been tested in a small number of people during a clinical trial. Here, Campo et al. use a technique known as panproteome array to scan samples collected during this trial, and identify which elements the body learns to recognize when it is exposed to the genetically manipulated strain of S. pneumoniae. The results show that when volunteers receive this vaccine, their body targets proteins that the capsule normally shields from the immune system. Many of these proteins are very similar across all strains of S. pneumoniae, which means that the whole cell vaccine could potentially better protect against a broad spectrum of bacteria. However, further studies are needed to assess whether this is the case, especially in infants.
Collapse
Affiliation(s)
| | - Timothy Q Le
- Antigen Discovery Inc, California, United States
| | | | | | - Andy A Teng
- Antigen Discovery Inc, California, United States
| | - Hervé Tettelin
- Institute for Genome Sciences, School of Medicine, University of Maryland, Baltimore, United States
| | | | - William P Hanage
- Center for Communicable Disease Dynamics, Department of Epidemiology, Harvard TH Chan School of Public Health, Boston, United States
| | | | - Xiaowu Liang
- Antigen Discovery Inc, California, United States
| | - Richard Malley
- Division of Infectious Diseases, Department of Medicine, Boston Children's Hospital and Harvard Medical School, Boston, United States
| | - Marc Lipsitch
- Center for Communicable Disease Dynamics, Department of Epidemiology, Harvard TH Chan School of Public Health, Boston, United States.,Department of Immunology and Infectious Diseases, Harvard TH Chan School of Public Health, Boston, United States
| | - Nicholas J Croucher
- MRC Centre for Global Infectious Disease Analysis, Department of Infectious Disease Epidemiology, Imperial College London, London, United Kingdom
| |
Collapse
|
39
|
Abstract
The standard workhorse for genomic analysis of the evolution of bacterial populations is phylogenetic modelling of mutations in the core genome. However, a notable amount of information about evolutionary and transmission processes in diverse populations can be lost unless the accessory genome is also taken into consideration. Here, we introduce panini (Pangenome Neighbour Identification for Bacterial Populations), a computationally scalable method for identifying the neighbours for each isolate in a data set using unsupervised machine learning with stochastic neighbour embedding based on the t-SNE (t-distributed stochastic neighbour embedding) algorithm. panini is browser-based and integrates with the Microreact platform for rapid online visualization and exploration of both core and accessory genome evolutionary signals, together with relevant epidemiological, geographical, temporal and other metadata. Several case studies with single- and multi-clone pneumococcal populations are presented to demonstrate the ability to identify biologically important signals from gene content data. panini is available at http://panini.pathogen.watch and code at http://gitlab.com/cgps/panini.
Collapse
Affiliation(s)
- Khalil Abudahab
- 1Centre for Genomic Pathogen Surveillance, Wellcome Genome Campus, Hinxton, UK
| | - Joaquín M Prada
- 2School of Veterinary Medicine, University of Surrey, Guildford, UK
| | - Zhirong Yang
- 3Department of Mathematics and Statistics, Helsinki Institute of Information Technology, University of Helsinki, FI-00014 Helsinki, Finland
| | | | - Nicholas J Croucher
- 5Department of Infectious Disease Epidemiology, Imperial College London, London, UK
| | - Jukka Corander
- 3Department of Mathematics and Statistics, Helsinki Institute of Information Technology, University of Helsinki, FI-00014 Helsinki, Finland.,6Department of Biostatistics, Institute of Basic Medical Sciences, University of Oslo, N-0317 Oslo, Norway
| | - David M Aanensen
- 1Centre for Genomic Pathogen Surveillance, Wellcome Genome Campus, Hinxton, UK.,7Big Data Institute, Li Ka Shing Centre for Health Informatics, University of Oxford, Oxford, UK
| |
Collapse
|
40
|
Abstract
Streptococcus pneumoniae (the pneumococcus) is a nasopharyngeal commensal and respiratory pathogen. Most isolates express a capsule, the species-wide diversity of which has been immunologically classified into ∼100 serotypes. Capsule polysaccharides have been combined into multivalent vaccines widely used in adults, but the T cell independence of the antibody response means they are not protective in infants. Polysaccharide conjugate vaccines (PCVs) trigger a T cell–dependent response through attaching a carrier protein to capsular polysaccharides. The immune response stimulated by PCVs in infants inhibits carriage of vaccine serotypes (VTs), resulting in population-wide herd immunity. These were replaced in carriage by non-VTs. Nevertheless, PCVs drove reductions in infant pneumococcal disease, due to the lower mean invasiveness of the postvaccination bacterial population; age-varying serotype invasiveness resulted in a smaller reduction in adult disease. Alternative vaccines being tested in trials are designed to provide species-wide protection through stimulating innate and cellular immune responses, alongside antibodies to conserved antigens.
Collapse
Affiliation(s)
- Nicholas J. Croucher
- Department of Infectious Disease Epidemiology, Imperial College London, London W2 1PG, United Kingdom
| | - Alessandra Løchen
- Department of Infectious Disease Epidemiology, Imperial College London, London W2 1PG, United Kingdom
| | - Stephen D. Bentley
- Infection Genomics Programme, Wellcome Sanger Institute, Hinxton, Cambridge CB10 1SA, United Kingdom
| |
Collapse
|
41
|
Abstract
The potential for genome-wide modelling of epistasis has recently surfaced given the possibility of sequencing densely sampled populations and the emerging families of statistical interaction models. Direct coupling analysis (DCA) has previously been shown to yield valuable predictions for single protein structures, and has recently been extended to genome-wide analysis of bacteria, identifying novel interactions in the co-evolution between resistance, virulence and core genome elements. However, earlier computational DCA methods have not been scalable to enable model fitting simultaneously to 104-105 polymorphisms, representing the amount of core genomic variation observed in analyses of many bacterial species. Here, we introduce a novel inference method (SuperDCA) that employs a new scoring principle, efficient parallelization, optimization and filtering on phylogenetic information to achieve scalability for up to 105 polymorphisms. Using two large population samples of Streptococcus pneumoniae, we demonstrate the ability of SuperDCA to make additional significant biological findings about this major human pathogen. We also show that our method can uncover signals of selection that are not detectable by genome-wide association analysis, even though our analysis does not require phenotypic measurements. SuperDCA, thus, holds considerable potential in building understanding about numerous organisms at a systems biological level.
Collapse
Affiliation(s)
- Santeri Puranen
- 2Department of Mathematics and Statistics, Helsinki Institute of Information Technology (HIIT), FI-00014 University of Helsinki, Finland.,1Department of Computer Science, Aalto University, FI-00076 Espoo, Finland
| | - Maiju Pesonen
- 1Department of Computer Science, Aalto University, FI-00076 Espoo, Finland.,2Department of Mathematics and Statistics, Helsinki Institute of Information Technology (HIIT), FI-00014 University of Helsinki, Finland
| | - Johan Pensar
- 2Department of Mathematics and Statistics, Helsinki Institute of Information Technology (HIIT), FI-00014 University of Helsinki, Finland
| | - Ying Ying Xu
- 1Department of Computer Science, Aalto University, FI-00076 Espoo, Finland.,2Department of Mathematics and Statistics, Helsinki Institute of Information Technology (HIIT), FI-00014 University of Helsinki, Finland
| | - John A Lees
- 3Pathogen Genomics, Wellcome Trust Sanger Institute, Cambridge CB10 1SA, UK
| | - Stephen D Bentley
- 3Pathogen Genomics, Wellcome Trust Sanger Institute, Cambridge CB10 1SA, UK
| | - Nicholas J Croucher
- 4Department of Infectious Disease Epidemiology, Imperial College London, London W2 1PG, UK
| | - Jukka Corander
- 5Department of Biostatistics, University of Oslo, 0317 Oslo, Norway.,2Department of Mathematics and Statistics, Helsinki Institute of Information Technology (HIIT), FI-00014 University of Helsinki, Finland.,3Pathogen Genomics, Wellcome Trust Sanger Institute, Cambridge CB10 1SA, UK
| |
Collapse
|
42
|
Wymant C, Blanquart F, Golubchik T, Gall A, Bakker M, Bezemer D, Croucher NJ, Hall M, Hillebregt M, Ong SH, Ratmann O, Albert J, Bannert N, Fellay J, Fransen K, Gourlay A, Grabowski MK, Gunsenheimer-Bartmeyer B, Günthard HF, Kivelä P, Kouyos R, Laeyendecker O, Liitsola K, Meyer L, Porter K, Ristola M, van Sighem A, Berkhout B, Cornelissen M, Kellam P, Reiss P, Fraser C. Easy and accurate reconstruction of whole HIV genomes from short-read sequence data with shiver. Virus Evol 2018; 4:vey007. [PMID: 29876136 PMCID: PMC5961307 DOI: 10.1093/ve/vey007] [Citation(s) in RCA: 44] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023] Open
Abstract
Studying the evolution of viruses and their molecular epidemiology relies on accurate viral sequence data, so that small differences between similar viruses can be meaningfully interpreted. Despite its higher throughput and more detailed minority variant data, next-generation sequencing has yet to be widely adopted for HIV. The difficulty of accurately reconstructing the consensus sequence of a quasispecies from reads (short fragments of DNA) in the presence of large between- and within-host diversity, including frequent indels, may have presented a barrier. In particular, mapping (aligning) reads to a reference sequence leads to biased loss of information; this bias can distort epidemiological and evolutionary conclusions. De novo assembly avoids this bias by aligning the reads to themselves, producing a set of sequences called contigs. However contigs provide only a partial summary of the reads, misassembly may result in their having an incorrect structure, and no information is available at parts of the genome where contigs could not be assembled. To address these problems we developed the tool shiver to pre-process reads for quality and contamination, then map them to a reference tailored to the sample using corrected contigs supplemented with the user's choice of existing reference sequences. Run with two commands per sample, it can easily be used for large heterogeneous data sets. We used shiver to reconstruct the consensus sequence and minority variant information from paired-end short-read whole-genome data produced with the Illumina platform, for sixty-five existing publicly available samples and fifty new samples. We show the systematic superiority of mapping to shiver's constructed reference compared with mapping the same reads to the closest of 3,249 real references: median values of 13 bases called differently and more accurately, 0 bases called differently and less accurately, and 205 bases of missing sequence recovered. We also successfully applied shiver to whole-genome samples of Hepatitis C Virus and Respiratory Syncytial Virus. shiver is publicly available from https://github.com/ChrisHIV/shiver.
Collapse
Affiliation(s)
- Chris Wymant
- Big Data Institute, Li Ka Shing Centre for Health Information and Discovery, Nuffield Department of Medicine, University of Oxford, Oxford, UK.,Medical Research Council Centre for Outbreak Analysis and Modelling, Department of Infectious Disease Epidemiology, Imperial College London, London, UK
| | - François Blanquart
- Medical Research Council Centre for Outbreak Analysis and Modelling, Department of Infectious Disease Epidemiology, Imperial College London, London, UK
| | - Tanya Golubchik
- Big Data Institute, Li Ka Shing Centre for Health Information and Discovery, Nuffield Department of Medicine, University of Oxford, Oxford, UK.,Wellcome Centre for Human Genetics, Nuffield Department of Medicine, University of Oxford, Oxford, UK
| | - Astrid Gall
- Department of Veterinary Medicine, University of Cambridge, Cambridge, UK.,Virus Genomics, Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton, Cambridge, UK
| | - Margreet Bakker
- Laboratory of Experimental Virology, Department of Medical Microbiology, Center for Infection and Immunity Amsterdam (CINIMA), Academic Medical Center of the University of Amsterdam, Amsterdam, The Netherlands
| | | | - Nicholas J Croucher
- Medical Research Council Centre for Outbreak Analysis and Modelling, Department of Infectious Disease Epidemiology, Imperial College London, London, UK
| | - Matthew Hall
- Big Data Institute, Li Ka Shing Centre for Health Information and Discovery, Nuffield Department of Medicine, University of Oxford, Oxford, UK.,Medical Research Council Centre for Outbreak Analysis and Modelling, Department of Infectious Disease Epidemiology, Imperial College London, London, UK
| | | | - Swee Hoe Ong
- Virus Genomics, Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton, Cambridge, UK
| | - Oliver Ratmann
- Medical Research Council Centre for Outbreak Analysis and Modelling, Department of Infectious Disease Epidemiology, Imperial College London, London, UK.,Department of Mathematics, Imperial College London, London, UK
| | - Jan Albert
- Department of Microbiology, Tumor and Cell Biology, Karolinska Institutet, Stockholm, Sweden.,Department of Clinical Microbiology, Karolinska University Hospital, Stockholm, Sweden
| | - Norbert Bannert
- Division for HIV and Other Retroviruses, Robert Koch Institute, Berlin, Germany
| | - Jacques Fellay
- School of Life Sciences, Ecole Polytechnique Fédérale de Lausanne, Lausanne, Switzerland.,Swiss Institute of Bioinformatics, Lausanne, Switzerland
| | - Katrien Fransen
- HIV/STI Reference Laboratory, Department of Clinical Science, WHO Collaborating Centre, Institute of Tropical Medicine, Antwerpen, Belgium
| | - Annabelle Gourlay
- Institute for Global Health, University College London, London, UK.,Department of Population Health, Faculty of Epidemiology and Population Health, London School of Hygiene and Tropical Medicine, London, UK
| | - M Kate Grabowski
- Department of Pathology, John Hopkins University, Baltimore, MD, USA
| | | | - Huldrych F Günthard
- Division of Infectious Diseases and Hospital Epidemiology, University Hospital Zurich, Zurich, Switzerland.,Institute of Medical Virology, University of Zurich, Zurich, Switzerland
| | - Pia Kivelä
- Department of Infectious Diseases, Helsinki University Hospital, Helsinki, Finland
| | - Roger Kouyos
- Division of Infectious Diseases and Hospital Epidemiology, University Hospital Zurich, Zurich, Switzerland.,Institute of Medical Virology, University of Zurich, Zurich, Switzerland
| | | | - Kirsi Liitsola
- Department of Infectious Diseases, Helsinki University Hospital, Helsinki, Finland
| | - Laurence Meyer
- INSERM CESP U1018, Université Paris Sud, Université Paris Saclay, APHP, Service de Santé Publique, Hôpital de Bicêtre, Le Kremlin-Bicêtre, France
| | - Kholoud Porter
- Institute for Global Health, University College London, London, UK
| | - Matti Ristola
- Department of Infectious Diseases, Helsinki University Hospital, Helsinki, Finland
| | | | - Ben Berkhout
- Laboratory of Experimental Virology, Department of Medical Microbiology, Center for Infection and Immunity Amsterdam (CINIMA), Academic Medical Center of the University of Amsterdam, Amsterdam, The Netherlands
| | - Marion Cornelissen
- Laboratory of Experimental Virology, Department of Medical Microbiology, Center for Infection and Immunity Amsterdam (CINIMA), Academic Medical Center of the University of Amsterdam, Amsterdam, The Netherlands
| | - Paul Kellam
- Kymab Ltd, Cambridge, UK.,Division of Infectious Diseases, Department of Medicine, Imperial College London, London, UK
| | - Peter Reiss
- Stichting HIV Monitoring, Amsterdam, The Netherlands.,Department of Global Health, Academic Medical Center and Amsterdam Institute for Global Health and Development, Amsterdam, The Netherlands
| | - Christophe Fraser
- Big Data Institute, Li Ka Shing Centre for Health Information and Discovery, Nuffield Department of Medicine, University of Oxford, Oxford, UK.,Medical Research Council Centre for Outbreak Analysis and Modelling, Department of Infectious Disease Epidemiology, Imperial College London, London, UK
| | | |
Collapse
|
43
|
De Ste Croix M, Vacca I, Kwun MJ, Ralph JD, Bentley SD, Haigh R, Croucher NJ, Oggioni MR. Phase-variable methylation and epigenetic regulation by type I restriction-modification systems. FEMS Microbiol Rev 2018; 41:S3-S15. [PMID: 28830092 DOI: 10.1093/femsre/fux025] [Citation(s) in RCA: 61] [Impact Index Per Article: 10.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/18/2017] [Accepted: 05/09/2017] [Indexed: 12/26/2022] Open
Abstract
Epigenetic modifications in bacteria, such as DNA methylation, have been shown to affect gene regulation, thereby generating cells that are isogenic but with distinctly different phenotypes. Restriction-modification (RM) systems contain prototypic methylases that are responsible for much of bacterial DNA methylation. This review focuses on a distinctive group of type I RM loci that , through phase variation, can modify their methylation target specificity and can thereby switch bacteria between alternative patterns of DNA methylation. Phase variation occurs at the level of the target recognition domains of the hsdS (specificity) gene via reversible recombination processes acting upon multiple hsdS alleles. We describe the global distribution of such loci throughout the prokaryotic kingdom and highlight the differences in loci structure across the various bacterial species. Although RM systems are often considered simply as an evolutionary response to bacteriophages, these multi-hsdS type I systems have also shown the capacity to change bacterial phenotypes. The ability of these RM systems to allow bacteria to reversibly switch between different physiological states, combined with the existence of such loci across many species of medical and industrial importance, highlights the potential of phase-variable DNA methylation to act as a global regulatory mechanism in bacteria.
Collapse
Affiliation(s)
| | - Irene Vacca
- Department of Genetics, University of Leicester, Leicester LE1 7RH, UK
| | - Min Jung Kwun
- Department of Infectious Disease Epidemiology, Imperial College London, London W2 1PG, UK
| | - Joseph D Ralph
- Department of Genetics, University of Leicester, Leicester LE1 7RH, UK
| | - Stephen D Bentley
- Infection Genomics, Wellcome Trust Sanger Institute, Hinxton CB10 1SA, UK
| | - Richard Haigh
- Department of Genetics, University of Leicester, Leicester LE1 7RH, UK
| | - Nicholas J Croucher
- Department of Infectious Disease Epidemiology, Imperial College London, London W2 1PG, UK
| | - Marco R Oggioni
- Department of Genetics, University of Leicester, Leicester LE1 7RH, UK
| |
Collapse
|
44
|
Corander J, Fraser C, Gutmann MU, Arnold B, Hanage WP, Bentley SD, Lipsitch M, Croucher NJ. Frequency-dependent selection in vaccine-associated pneumococcal population dynamics. Nat Ecol Evol 2017; 1:1950-1960. [PMID: 29038424 PMCID: PMC5708525 DOI: 10.1038/s41559-017-0337-x] [Citation(s) in RCA: 79] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2017] [Accepted: 09/01/2017] [Indexed: 12/21/2022]
Abstract
Many bacterial species are composed of multiple lineages distinguished by extensive variation in gene content. These often cocirculate in the same habitat, but the evolutionary and ecological processes that shape these complex populations are poorly understood. Addressing these questions is particularly important for Streptococcus pneumoniae, a nasopharyngeal commensal and respiratory pathogen, because the changes in population structure associated with the recent introduction of partial-coverage vaccines have substantially reduced pneumococcal disease. Here we show that pneumococcal lineages from multiple populations each have a distinct combination of intermediate-frequency genes. Functional analysis suggested that these loci may be subject to negative frequency-dependent selection (NFDS) through interactions with other bacteria, hosts or mobile elements. Correspondingly, these genes had similar frequencies in four populations with dissimilar lineage compositions. These frequencies were maintained following substantial alterations in lineage prevalences once vaccination programmes began. Fitting a multilocus NFDS model of post-vaccine population dynamics to three genomic datasets using Approximate Bayesian Computation generated reproducible estimates of the influence of NFDS on pneumococcal evolution, the strength of which varied between loci. Simulations replicated the stable frequency of lineages unperturbed by vaccination, patterns of serotype switching and clonal replacement. This framework highlights how bacterial ecology affects the impact of clinical interventions.
Collapse
Affiliation(s)
- Jukka Corander
- Helsinki Institute for Information Technology, Department of Mathematics and Statistics, University of Helsinki, 00014, Helsinki, Finland
- Department of Biostatistics, University of Oslo, 0317, Oslo, Norway
- Infection Genomics, The Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge, CB10 1SA, UK
| | - Christophe Fraser
- Big Data Institute, Nuffield Department of Medicine, University of Oxford, Oxford, OX3 7LF, UK
| | - Michael U Gutmann
- School of Informatics, University of Edinburgh, Edinburgh, EH8 9AB, UK
| | - Brian Arnold
- Center for Communicable Disease Dynamics, Harvard T. H. Chan School of Public Health, 677 Huntington Avenue, Boston, MA, 02115, USA
| | - William P Hanage
- Center for Communicable Disease Dynamics, Harvard T. H. Chan School of Public Health, 677 Huntington Avenue, Boston, MA, 02115, USA
| | - Stephen D Bentley
- Infection Genomics, The Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge, CB10 1SA, UK
| | - Marc Lipsitch
- Center for Communicable Disease Dynamics, Harvard T. H. Chan School of Public Health, 677 Huntington Avenue, Boston, MA, 02115, USA
- Departments of Epidemiology and Immunology and Infectious Diseases, Harvard T. H. Chan School of Public Health, 677 Huntington Avenue, Boston, MA, 02115, USA
| | - Nicholas J Croucher
- MRC Centre for Outbreak Analysis and Modelling, Department of Infectious Disease Epidemiology, Imperial College London, London, W2 1PG, UK.
| |
Collapse
|
45
|
Abstract
Bacterial transformation can insert or delete genomic islands (GIs), depending on the donor and recipient genotypes, if an homologous recombination spans the GI’s integration site and includes sufficiently long flanking homologous arms. Combining mathematical models of recombination with experiments using pneumococci found GI insertion rates declined geometrically with the GI’s size. The decrease in acquisition frequency with length (1.08×10−3 bp−1) was higher than a previous estimate of the analogous rate at which core genome recombinations terminated. Although most efficient for shorter GIs, transformation-mediated deletion frequencies did not vary consistently with GI length, with removal of 10-kb GIs ∼50% as efficient as acquisition of base substitutions. Fragments of 2 kb, typical of transformation event sizes, could drive all these deletions independent of island length. The strong asymmetry of transformation, and its capacity to efficiently remove GIs, suggests nonmobile accessory loci will decline in frequency without preservation by selection.
Collapse
Affiliation(s)
- Katinka J Apagyi
- MRC Centre for Outbreak Analysis and Modelling, Department of Infectious Disease Epidemiology, Imperial College London, London, United Kingdom
| | - Christophe Fraser
- Big Data Institute, Nuffield Department of Medicine, University of Oxford, Oxford, United Kingdom
| | - Nicholas J Croucher
- MRC Centre for Outbreak Analysis and Modelling, Department of Infectious Disease Epidemiology, Imperial College London, London, United Kingdom
| |
Collapse
|
46
|
Mostowy RJ, Croucher NJ, De Maio N, Chewapreecha C, Salter SJ, Turner P, Aanensen DM, Bentley SD, Didelot X, Fraser C. Pneumococcal Capsule Synthesis Locus cps as Evolutionary Hotspot with Potential to Generate Novel Serotypes by Recombination. Mol Biol Evol 2017; 34:2537-2554. [PMID: 28595308 PMCID: PMC5850285 DOI: 10.1093/molbev/msx173] [Citation(s) in RCA: 50] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Abstract
Diversity of the polysaccharide capsule in Streptococcus pneumoniae-main surface antigen and the target of the currently used pneumococcal vaccines-constitutes a major obstacle in eliminating pneumococcal disease. Such diversity is genetically encoded by almost 100 variants of the capsule biosynthesis locus, cps. However, the evolutionary dynamics of the capsule remains not fully understood. Here, using genetic data from 4,519 bacterial isolates, we found cps to be an evolutionary hotspot with elevated substitution and recombination rates. These rates were a consequence of relaxed purifying selection and positive, diversifying selection acting at this locus, supporting the hypothesis that the capsule has an increased potential to generate novel diversity compared with the rest of the genome. Diversifying selection was particularly evident in the region of wzd/wze genes, which are known to regulate capsule expression and hence the bacterium's ability to cause disease. Using a novel, capsule-centered approach, we analyzed the evolutionary history of 12 major serogroups. Such analysis revealed their complex diversification scenarios, which were principally driven by recombination with other serogroups and other streptococci. Patterns of recombinational exchanges between serogroups could not be explained by serotype frequency alone, thus pointing to nonrandom associations between co-colonizing serotypes. Finally, we discovered a previously unobserved mosaic serotype 39X, which was confirmed to carry a viable and structurally novel capsule. Adding to previous discoveries of other mosaic capsules in densely sampled collections, these results emphasize the strong adaptive potential of the bacterium by its ability to generate novel antigenic diversity by recombination.
Collapse
Affiliation(s)
- Rafał J. Mostowy
- Department of Infectious Disease Epidemiology, School of Public Health, Imperial College London, London, United Kingdom
| | - Nicholas J. Croucher
- Department of Infectious Disease Epidemiology, School of Public Health, Imperial College London, London, United Kingdom
| | - Nicola De Maio
- Nuffield Department of Medicine, University of Oxford, Oxford, United Kingdom
- Institute for Emerging Infections, Oxford Martin School, Oxford, United Kingdom
| | - Claire Chewapreecha
- Department of Medicine, University of Cambridge, Cambridge, United Kingdom
- Bioinformatics and Systems Biology Program, School of Bioresources and Technology, King Mongkut’s University of Technology Thonburi, Bangkok, Thailand
| | - Susannah J. Salter
- The Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge, United Kingdom
| | - Paul Turner
- Nuffield Department of Medicine, Centre for Tropical Medicine and Global Health, University of Oxford, Oxford, United Kingdom
- Cambodia-Oxford Medical Research Unit, Angkor Hospital for Children, Siem Reap, Cambodia
| | - David M. Aanensen
- Department of Infectious Disease Epidemiology, School of Public Health, Imperial College London, London, United Kingdom
- Centre for Genomic Pathogen Surveillance, Wellcome Genome Campus, Hinxton, Cambridge, United Kingdom
| | - Stephen D. Bentley
- The Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge, United Kingdom
| | - Xavier Didelot
- Department of Infectious Disease Epidemiology, School of Public Health, Imperial College London, London, United Kingdom
| | - Christophe Fraser
- Department of Infectious Disease Epidemiology, School of Public Health, Imperial College London, London, United Kingdom
- Nuffield Department of Medicine, Li Ka Shing Centre for Health Information and Discovery, Oxford Big Data Institute, University of Oxford, Oxford, United Kingdom
| |
Collapse
|
47
|
Hadfield J, Croucher NJ, Goater RJ, Abudahab K, Aanensen DM, Harris SR. Phandango: an interactive viewer for bacterial population genomics. Bioinformatics 2017; 34:292-293. [PMID: 29028899 PMCID: PMC5860215 DOI: 10.1093/bioinformatics/btx610] [Citation(s) in RCA: 343] [Impact Index Per Article: 49.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/24/2017] [Revised: 09/06/2017] [Accepted: 09/22/2017] [Indexed: 12/13/2022] Open
Abstract
Summary Fully exploiting the wealth of data in current bacterial population genomics datasets requires synthesizing and integrating different types of analysis across millions of base pairs in hundreds or thousands of isolates. Current approaches often use static representations of phylogenetic, epidemiological, statistical and evolutionary analysis results that are difficult to relate to one another. Phandango is an interactive application running in a web browser allowing fast exploration of large-scale population genomics datasets combining the output from multiple genomic analysis methods in an intuitive and interactive manner. Availability and implementation Phandango is a web application freely available for use at www.phandango.net and includes a diverse collection of datasets as examples. Source code together with a detailed wiki page is available on GitHub at https://github.com/jameshadfield/phandango.
Collapse
Affiliation(s)
- James Hadfield
- Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Cambridge, UK
| | - Nicholas J Croucher
- Department of Infectious Disease Epidemiology, Imperial College London, London, UK
| | - Richard J Goater
- Centre for Genomic Pathogen Surveillance, Wellcome Trust Genome Campus, Cambridge, UK
| | - Khalil Abudahab
- Centre for Genomic Pathogen Surveillance, Wellcome Trust Genome Campus, Cambridge, UK
| | - David M Aanensen
- Department of Infectious Disease Epidemiology, Imperial College London, London, UK.,Centre for Genomic Pathogen Surveillance, Wellcome Trust Genome Campus, Cambridge, UK
| | - Simon R Harris
- Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Cambridge, UK
| |
Collapse
|
48
|
Mostowy R, Croucher NJ, Andam CP, Corander J, Hanage WP, Marttinen P. Efficient Inference of Recent and Ancestral Recombination within Bacterial Populations. Mol Biol Evol 2017; 34:1167-1182. [PMID: 28199698 PMCID: PMC5400400 DOI: 10.1093/molbev/msx066] [Citation(s) in RCA: 105] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023] Open
Abstract
Prokaryotic evolution is affected by horizontal transfer of genetic material through recombination. Inference of an evolutionary tree of bacteria thus relies on accurate identification of the population genetic structure and recombination-derived mosaicism. Rapidly growing databases represent a challenge for computational methods to detect recombinations in bacterial genomes. We introduce a novel algorithm called fastGEAR which identifies lineages in diverse microbial alignments, and recombinations between them and from external origins. The algorithm detects both recent recombinations (affecting a few isolates) and ancestral recombinations between detected lineages (affecting entire lineages), thus providing insight into recombinations affecting deep branches of the phylogenetic tree. In simulations, fastGEAR had comparable power to detect recent recombinations and outstanding power to detect the ancestral ones, compared with state-of-the-art methods, often with a fraction of computational cost. We demonstrate the utility of the method by analyzing a collection of 616 whole-genomes of a recombinogenic pathogen Streptococcus pneumoniae, for which the method provided a high-resolution view of recombination across the genome. We examined in detail the penicillin-binding genes across the Streptococcus genus, demonstrating previously undetected genetic exchanges between different species at these three loci. Hence, fastGEAR can be readily applied to investigate mosaicism in bacterial genes across multiple species. Finally, fastGEAR correctly identified many known recombination hotspots and pointed to potential new ones. Matlab code and Linux/Windows executables are available at https://users.ics.aalto.fi/~pemartti/fastGEAR/ (last accessed February 6, 2017).
Collapse
Affiliation(s)
- Rafal Mostowy
- Department of Infectious Disease Epidemiology, St. Mary's Campus, Imperial College London, London, United Kingdom
| | - Nicholas J Croucher
- Department of Infectious Disease Epidemiology, St. Mary's Campus, Imperial College London, London, United Kingdom
| | - Cheryl P Andam
- Department of Epidemiology, Harvard TH Chan School of Public Health, Center for Communicable Disease Dynamics, Boston, MA
| | - Jukka Corander
- Department of Mathematics and Statistics, Helsinki Institute for Information Technology HIIT, University of Helsinki, Helsinki, Finland.,Department of Biostatistics, University of Oslo, Oslo, Norway
| | - William P Hanage
- Department of Epidemiology, Harvard TH Chan School of Public Health, Center for Communicable Disease Dynamics, Boston, MA
| | - Pekka Marttinen
- Department of Computer Science, Helsinki Institute for Information Technology HIIT, Aalto University, Espoo, Finland
| |
Collapse
|
49
|
Lees JA, Croucher NJ, Goldblatt D, Nosten F, Parkhill J, Turner C, Turner P, Bentley SD. Genome-wide identification of lineage and locus specific variation associated with pneumococcal carriage duration. eLife 2017; 6:e26255. [PMID: 28742023 PMCID: PMC5576492 DOI: 10.7554/elife.26255] [Citation(s) in RCA: 58] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2017] [Accepted: 07/21/2017] [Indexed: 01/04/2023] Open
Abstract
Streptococcus pneumoniae is a leading cause of invasive disease in infants, especially in low-income settings. Asymptomatic carriage in the nasopharynx is a prerequisite for disease, but variability in its duration is currently only understood at the serotype level. Here we developed a model to calculate the duration of carriage episodes from longitudinal swab data, and combined these results with whole genome sequence data. We estimated that pneumococcal genomic variation accounted for 63% of the phenotype variation, whereas the host traits considered here (age and previous carriage) accounted for less than 5%. We further partitioned this heritability into both lineage and locus effects, and quantified the amount attributable to the largest sources of variation in carriage duration: serotype (17%), drug-resistance (9%) and other significant locus effects (7%). A pan-genome-wide association study identified prophage sequences as being associated with decreased carriage duration independent of serotype, potentially by disruption of the competence mechanism. These findings support theoretical models of pneumococcal competition and antibiotic resistance.
Collapse
Affiliation(s)
- John A Lees
- Infection GenomicsWellcome Trust Sanger InstituteHinxtonUnited Kingdom
| | - Nicholas J Croucher
- Department of Infectious Disease EpidemiologySt. Mary’s Campus, Imperial College LondonLondonUnited Kingdom
| | - David Goldblatt
- Institute of Child HealthUniversity College LondonLondonUnited Kingdom
| | - François Nosten
- Shoklo Malaria Research Unit, Mahidol-Oxford Tropical Medicine Research Unit, Faculty of Tropical MedicineMahidol UniversityMae SotThailand
- Centre for Tropical Medicine and Global Health, Nuffield Department of MedicineUniversity of OxfordOxfordUnited Kingdom
| | - Julian Parkhill
- Infection GenomicsWellcome Trust Sanger InstituteHinxtonUnited Kingdom
| | - Claudia Turner
- Shoklo Malaria Research Unit, Mahidol-Oxford Tropical Medicine Research Unit, Faculty of Tropical MedicineMahidol UniversityMae SotThailand
- Centre for Tropical Medicine and Global Health, Nuffield Department of MedicineUniversity of OxfordOxfordUnited Kingdom
| | - Paul Turner
- Shoklo Malaria Research Unit, Mahidol-Oxford Tropical Medicine Research Unit, Faculty of Tropical MedicineMahidol UniversityMae SotThailand
- Centre for Tropical Medicine and Global Health, Nuffield Department of MedicineUniversity of OxfordOxfordUnited Kingdom
| | - Stephen D Bentley
- Infection GenomicsWellcome Trust Sanger InstituteHinxtonUnited Kingdom
| |
Collapse
|
50
|
Skwark MJ, Croucher NJ, Puranen S, Chewapreecha C, Pesonen M, Xu YY, Turner P, Harris SR, Beres SB, Musser JM, Parkhill J, Bentley SD, Aurell E, Corander J. Interacting networks of resistance, virulence and core machinery genes identified by genome-wide epistasis analysis. PLoS Genet 2017; 13:e1006508. [PMID: 28207813 PMCID: PMC5312804 DOI: 10.1371/journal.pgen.1006508] [Citation(s) in RCA: 74] [Impact Index Per Article: 10.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2016] [Accepted: 11/24/2016] [Indexed: 12/05/2022] Open
Abstract
Recent advances in the scale and diversity of population genomic datasets for bacteria now provide the potential for genome-wide patterns of co-evolution to be studied at the resolution of individual bases. Here we describe a new statistical method, genomeDCA, which uses recent advances in computational structural biology to identify the polymorphic loci under the strongest co-evolutionary pressures. We apply genomeDCA to two large population data sets representing the major human pathogens Streptococcus pneumoniae (pneumococcus) and Streptococcus pyogenes (group A Streptococcus). For pneumococcus we identified 5,199 putative epistatic interactions between 1,936 sites. Over three-quarters of the links were between sites within the pbp2x, pbp1a and pbp2b genes, the sequences of which are critical in determining non-susceptibility to beta-lactam antibiotics. A network-based analysis found these genes were also coupled to that encoding dihydrofolate reductase, changes to which underlie trimethoprim resistance. Distinct from these antibiotic resistance genes, a large network component of 384 protein coding sequences encompassed many genes critical in basic cellular functions, while another distinct component included genes associated with virulence. The group A Streptococcus (GAS) data set population represents a clonal population with relatively little genetic variation and a high level of linkage disequilibrium across the genome. Despite this, we were able to pinpoint two RNA pseudouridine synthases, which were each strongly linked to a separate set of loci across the chromosome, representing biologically plausible targets of co-selection. The population genomic analysis method applied here identifies statistically significantly co-evolving locus pairs, potentially arising from fitness selection interdependence reflecting underlying protein-protein interactions, or genes whose product activities contribute to the same phenotype. This discovery approach greatly enhances the future potential of epistasis analysis for systems biology, and can complement genome-wide association studies as a means of formulating hypotheses for targeted experimental work. Epistatic interactions between polymorphisms in DNA are recognized as important drivers of evolution in numerous organisms. Study of epistasis in bacteria has been hampered by the lack of densely sampled population genomic data, suitable statistical models and inference algorithms sufficiently powered for extremely high-dimensional parameter spaces. We introduce the first model-based method for genome-wide epistasis analysis and use two of the largest available bacterial population genome data sets on Streptococcus pneumoniae (the pneumococcus) and Streptococcus pyogenes (group A Streptococcus) to demonstrate its potential for biological discovery. Our approach reveals interacting networks of resistance, virulence and core machinery genes in the pneumococcus, which highlights putative candidates for novel drug targets. We also discover a number of plausible targets of co-selection in S. pyogenes linked to RNA pseudouridine synthases. Our method significantly enhances the future potential of epistasis analysis for systems biology, and can complement genome-wide association studies as a means of formulating hypotheses for targeted experimental work.
Collapse
Affiliation(s)
- Marcin J Skwark
- Department of Chemistry, Vanderbilt University, Nashville, TN, United States of America
| | - Nicholas J Croucher
- Department of Infectious Disease Epidemiology, Imperial College London, London, United Kingdom
| | - Santeri Puranen
- Department of Computer Science, Aalto University, Espoo, Finland
| | | | - Maiju Pesonen
- Department of Computer Science, Aalto University, Espoo, Finland
| | - Ying Ying Xu
- Department of Computer Science, Aalto University, Espoo, Finland
| | - Paul Turner
- Shoklo Malaria Research Unit, Mahidol-Oxford Tropical Medicine Research Unit, Faculty of Tropical Medicine, Mahidol University, Mae Sot, Thailand.,Centre for Tropical Medicine, Nuffield Department of Medicine, University of Oxford, Oxford, United Kingdom
| | - Simon R Harris
- Pathogen Genomics, Wellcome Trust Sanger Institute, Cambridge, United Kingdom
| | - Stephen B Beres
- Center for Molecular and Translational Human Infectious Diseases Research, Department of Pathology and Genomic Medicine, Houston Methodist Research Institute, and Houston Methodist Hospital, Houston, Texas, United States of America
| | - James M Musser
- Center for Molecular and Translational Human Infectious Diseases Research, Department of Pathology and Genomic Medicine, Houston Methodist Research Institute, and Houston Methodist Hospital, Houston, Texas, United States of America.,Departments of Pathology and Laboratory Medicine and Microbiology and Immunology, Weill Cornell Medical College, New York, New York, United States of America
| | - Julian Parkhill
- Pathogen Genomics, Wellcome Trust Sanger Institute, Cambridge, United Kingdom
| | - Stephen D Bentley
- Pathogen Genomics, Wellcome Trust Sanger Institute, Cambridge, United Kingdom
| | - Erik Aurell
- Department of Computational Biology, KTH-Royal Institute of Technology, Stockholm, Sweden.,Departments of Applied Physics and Computer Science, Aalto University, Espoo, Finland.,Institute of Theoretical Physics, Chinese Academy of Sciences, Beijing, China
| | - Jukka Corander
- Pathogen Genomics, Wellcome Trust Sanger Institute, Cambridge, United Kingdom.,Department of Mathematics and Statistics, University of Helsinki, Helsinki, Finland.,Department of Biostatistics, University of Oslo, Oslo, Norway.,Department of Veterinary Medicine, University of Cambridge, Cambridge, United Kingdom
| |
Collapse
|