1
|
Teichmann L, Wenne M, Luitwieler S, Dugar G, Bengtsson-Palme J, ter Kuile B. Genetic adaptation to amoxicillin in Escherichia coli: The limited role of dinB and katE. PLoS One 2025; 20:e0312223. [PMID: 39970152 PMCID: PMC11838884 DOI: 10.1371/journal.pone.0312223] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2024] [Accepted: 12/27/2024] [Indexed: 02/21/2025] Open
Abstract
Bacteria can quickly adapt to sub-lethal concentrations of antibiotics. Several stress and DNA repair genes contribute to this adaptation process. However, the pathways leading to adaptation by acquisition of de novo mutations remain poorly understood. This study explored the roles of DNA polymerase IV (dinB) and catalase HP2 (katE) in E. coli's adaptation to amoxicillin. These genes are thought to play essential roles in beta-lactam resistance-dinB in increasing mutation rates and katE in managing oxidative stress. By comparing the adaptation rates, transcriptomic profiles, and genetic changes of wild-type and knockout strains, we aimed to clarify the contributions of these genes to beta-lactam resistance. While all strains exhibited similar adaptation rates and mutations in the frdD gene and ampC operon, several unique mutations were acquired in the ΔkatE and ΔdinB strains. Overall, this study distinguishes the contributions of general stress-related genes on the one hand, and dinB, and katE on the other hand, in development of beta-lactam resistance.
Collapse
Affiliation(s)
- Lisa Teichmann
- Swammerdam Institute of Life Sciences, Molecular Biology and Microbial Food Safety, University of Amsterdam, Amsterdam, The Netherlands
| | - Marcus Wenne
- Department of Life Sciences, Division of Systems and Synthetic Biology, SciLifeLab, Chalmers University of Technology, Gothenburg, Sweden
- Centre for Antibiotic Resistance Research (CARe) in Gothenburg, Gothenburg, Sweden
| | - Sam Luitwieler
- Swammerdam Institute of Life Sciences, Molecular Biology and Microbial Food Safety, University of Amsterdam, Amsterdam, The Netherlands
| | - Gaurav Dugar
- Swammerdam Institute of Life Sciences, Molecular Biology and Microbial Food Safety, University of Amsterdam, Amsterdam, The Netherlands
| | - Johan Bengtsson-Palme
- Department of Life Sciences, Division of Systems and Synthetic Biology, SciLifeLab, Chalmers University of Technology, Gothenburg, Sweden
- Centre for Antibiotic Resistance Research (CARe) in Gothenburg, Gothenburg, Sweden
- Department of Infectious Diseases, Institute of Biomedicine, University of Gothenburg, Gothenburg, Sweden
| | - Benno ter Kuile
- Swammerdam Institute of Life Sciences, Molecular Biology and Microbial Food Safety, University of Amsterdam, Amsterdam, The Netherlands
| |
Collapse
|
2
|
Pham NP, Gingras H, Godin C, Feng J, Groppi A, Nikolski M, Leprohon P, Ouellette M. Holistic understanding of trimethoprim resistance in Streptococcus pneumoniae using an integrative approach of genome-wide association study, resistance reconstruction, and machine learning. mBio 2024; 15:e0136024. [PMID: 39120145 PMCID: PMC11389379 DOI: 10.1128/mbio.01360-24] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2024] [Accepted: 07/08/2024] [Indexed: 08/10/2024] Open
Abstract
Antimicrobial resistance (AMR) is a public health threat worldwide. Next-generation sequencing (NGS) has opened unprecedented opportunities to accelerate AMR mechanism discovery and diagnostics. Here, we present an integrative approach to investigate trimethoprim (TMP) resistance in the key pathogen Streptococcus pneumoniae. We explored a collection of 662 S. pneumoniae genomes by conducting a genome-wide association study (GWAS), followed by functional validation using resistance reconstruction experiments, combined with machine learning (ML) approaches to predict TMP minimum inhibitory concentration (MIC). Our study showed that multiple additive mutations in the folA and sulA loci are responsible for TMP non-susceptibility in S. pneumoniae and can be used as key features to build ML models for digital MIC prediction, reaching an average accuracy within ±1 twofold dilution factor of 86.3%. Our roadmap of in silico analysis-wet-lab validation-diagnostic tool building could be adapted to explore AMR in other combinations of bacteria-antibiotic. IMPORTANCE In the age of next-generation sequencing (NGS), while data-driven methods such as genome-wide association study (GWAS) and machine learning (ML) excel at finding patterns, functional validation can be challenging due to the high numbers of candidate variants. We designed an integrative approach combining a GWAS on S. pneumoniae clinical isolates, followed by whole-genome transformation coupled with NGS to functionally characterize a large set of GWAS candidates. Our study validated several phenotypic folA mutations beyond the standard Ile100Leu mutation, and showed that the overexpression of the sulA locus produces trimethoprim (TMP) resistance in Streptococcus pneumoniae. These validated loci, when used to build ML models, were found to be the best inputs for predicting TMP minimal inhibitory concentrations. Integrative approaches can bridge the genotype-phenotype gap by biological insights that can be incorporated in ML models for accurate prediction of drug susceptibility.
Collapse
Affiliation(s)
- Nguyen-Phuong Pham
- Centre de Recherche en Infectiologie du Centre de Recherche du CHU de Québec and Département de Microbiologie, Infectiologie et Immunologie, Faculté de Médecine, Université Laval, Québec City, Québec, Canada
| | - Hélène Gingras
- Centre de Recherche en Infectiologie du Centre de Recherche du CHU de Québec and Département de Microbiologie, Infectiologie et Immunologie, Faculté de Médecine, Université Laval, Québec City, Québec, Canada
| | - Chantal Godin
- Centre de Recherche en Infectiologie du Centre de Recherche du CHU de Québec and Département de Microbiologie, Infectiologie et Immunologie, Faculté de Médecine, Université Laval, Québec City, Québec, Canada
| | - Jie Feng
- State Key Laboratory of Microbial Resources, Institute of Microbiology, Chinese Academy of Sciences, Beijing, China
| | - Alexis Groppi
- Bordeaux Bioinformatics Center and CNRS, Institut de Biochimie et Génétique Cellulaires (IBGC) UMR 5095, Université de Bordeaux, Bordeaux, France
| | - Macha Nikolski
- Bordeaux Bioinformatics Center and CNRS, Institut de Biochimie et Génétique Cellulaires (IBGC) UMR 5095, Université de Bordeaux, Bordeaux, France
| | - Philippe Leprohon
- Centre de Recherche en Infectiologie du Centre de Recherche du CHU de Québec and Département de Microbiologie, Infectiologie et Immunologie, Faculté de Médecine, Université Laval, Québec City, Québec, Canada
| | - Marc Ouellette
- Centre de Recherche en Infectiologie du Centre de Recherche du CHU de Québec and Département de Microbiologie, Infectiologie et Immunologie, Faculté de Médecine, Université Laval, Québec City, Québec, Canada
| |
Collapse
|
3
|
Jugas R, Vitkova H. ProcaryaSV: structural variation detection pipeline for bacterial genomes using short-read sequencing. BMC Bioinformatics 2024; 25:233. [PMID: 38982375 PMCID: PMC11234778 DOI: 10.1186/s12859-024-05843-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2024] [Accepted: 06/13/2024] [Indexed: 07/11/2024] Open
Abstract
BACKGROUND Structural variations play an important role in bacterial genomes. They can mediate genome adaptation quickly in response to the external environment and thus can also play a role in antibiotic resistance. The detection of structural variations in bacteria is challenging, and the recognition of even small rearrangements can be important. Even though most detection tools are aimed at and benchmarked on eukaryotic genomes, they can also be used on prokaryotic genomes. The key features of detection are the ability to detect small rearrangements and support haploid genomes. Because of the limiting performance of a single detection tool, combining the detection abilities of multiple tools can lead to more robust results. There are already available workflows for structural variation detection for long-reads technologies and for the detection of single-nucleotide variation and indels, both aimed at bacteria. Yet we are unaware of structural variations detection workflows for the short-reads sequencing platform. Motivated by this gap we created our workflow. Further, we were interested in increasing the detection performance and providing more robust results. RESULTS We developed an open-source bioinformatics pipeline, ProcaryaSV, for the detection of structural variations in bacterial isolates from paired-end short sequencing reads. Multiple tools, starting with quality control and trimming of sequencing data, alignment to the reference genome, and multiple structural variation detection tools, are integrated. All the partial results are then processed and merged with an in-house merging algorithm. Compared with a single detection approach, ProcaryaSV has improved detection performance and is a reproducible easy-to-use tool. CONCLUSIONS The ProcaryaSV pipeline provides an integrative approach to structural variation detection from paired-end next-generation sequencing of bacterial samples. It can be easily installed and used on Linux machines. It is publicly available on GitHub at https://github.com/robinjugas/ProcaryaSV .
Collapse
Affiliation(s)
- Robin Jugas
- Department of Biomedical Engineering, Brno University of Technology, Brno, Czech Republic
| | - Helena Vitkova
- Department of Biomedical Engineering, Brno University of Technology, Brno, Czech Republic.
| |
Collapse
|
4
|
Wang Z, Gu D, Hong Y, Hu Y, Gu J, Tang Y, Zhou X, Zhang Y, Jiao X, Li Q. Microevolution of Salmonella 4,[5],12:i:- derived from Salmonella enterica serovar Typhimurium through complicated transpositions. Cell Rep 2023; 42:113227. [PMID: 37837619 DOI: 10.1016/j.celrep.2023.113227] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2023] [Revised: 08/28/2023] [Accepted: 09/21/2023] [Indexed: 10/16/2023] Open
Abstract
Salmonella enterica subsp. enterica serovar 4,[5],12:i:- (Salmonella 4,[5],12:i:-), derived from S. Typhimurium, has become the dominant serotype causing human salmonellosis. In this study, we define the genetic mechanism of the generation of Salmonella 4,[5],12:i:- from S. Typhimurium through complicated transpositions and demonstrate that Salmonella 4,[5],12:i:- displays more efficient colonization and survival abilities in mice than its parent S. Typhimurium strain. We identified intermediate strains carrying both resistance regions (RRs) and the fljAB operon for the generation of Salmonella 4,[5],12:i:-. The insertion of RR3 into the chromosomal hin-iroB site of S. Typhimurium produced RR3-S. Typhimurium as a primary intermediate. Salmonella 4,[5],12:i:- was then produced by replacing the fljAB operon and/or its flanking sequences through intramolecular transpositions mediated by IS26 and/or IS1R elements in RR3-S. Typhimurium, which was further confirmed both in vitro and in vivo. Overall, we demonstrate the molecular mechanism underlying the origin, generation, and advantage of RRs-Salmonella 4,[5],12:i:- from S. Typhimurium.
Collapse
Affiliation(s)
- Zhenyu Wang
- Key Laboratory of Prevention and Control of Biological Hazard Factors (Animal Origin) for Agri-food Safety and Quality, Ministry of Agriculture of China, Yangzhou University, 48 Wenhui East Road, Yangzhou 225000, Jiangsu, China; Jiangsu Key Lab of Zoonosis/Jiangsu Co-innovation Center for Prevention and Control of Important Animal Infectious Diseases and Zoonoses, Yangzhou University, 48 Wenhui East Road, Yangzhou 225000, Jiangsu, China; Joint International Research Laboratory of Agriculture and Agri-product Safety, Yangzhou University, 48 Wenhui East Road, Yangzhou 225000, Jiangsu, China
| | - Dan Gu
- Key Laboratory of Prevention and Control of Biological Hazard Factors (Animal Origin) for Agri-food Safety and Quality, Ministry of Agriculture of China, Yangzhou University, 48 Wenhui East Road, Yangzhou 225000, Jiangsu, China; Jiangsu Key Lab of Zoonosis/Jiangsu Co-innovation Center for Prevention and Control of Important Animal Infectious Diseases and Zoonoses, Yangzhou University, 48 Wenhui East Road, Yangzhou 225000, Jiangsu, China; Joint International Research Laboratory of Agriculture and Agri-product Safety, Yangzhou University, 48 Wenhui East Road, Yangzhou 225000, Jiangsu, China
| | - Yaming Hong
- Key Laboratory of Prevention and Control of Biological Hazard Factors (Animal Origin) for Agri-food Safety and Quality, Ministry of Agriculture of China, Yangzhou University, 48 Wenhui East Road, Yangzhou 225000, Jiangsu, China; Jiangsu Key Lab of Zoonosis/Jiangsu Co-innovation Center for Prevention and Control of Important Animal Infectious Diseases and Zoonoses, Yangzhou University, 48 Wenhui East Road, Yangzhou 225000, Jiangsu, China; Joint International Research Laboratory of Agriculture and Agri-product Safety, Yangzhou University, 48 Wenhui East Road, Yangzhou 225000, Jiangsu, China
| | - Yachen Hu
- Key Laboratory of Prevention and Control of Biological Hazard Factors (Animal Origin) for Agri-food Safety and Quality, Ministry of Agriculture of China, Yangzhou University, 48 Wenhui East Road, Yangzhou 225000, Jiangsu, China; Jiangsu Key Lab of Zoonosis/Jiangsu Co-innovation Center for Prevention and Control of Important Animal Infectious Diseases and Zoonoses, Yangzhou University, 48 Wenhui East Road, Yangzhou 225000, Jiangsu, China; Joint International Research Laboratory of Agriculture and Agri-product Safety, Yangzhou University, 48 Wenhui East Road, Yangzhou 225000, Jiangsu, China
| | - Jiaojie Gu
- Key Laboratory of Prevention and Control of Biological Hazard Factors (Animal Origin) for Agri-food Safety and Quality, Ministry of Agriculture of China, Yangzhou University, 48 Wenhui East Road, Yangzhou 225000, Jiangsu, China; Jiangsu Key Lab of Zoonosis/Jiangsu Co-innovation Center for Prevention and Control of Important Animal Infectious Diseases and Zoonoses, Yangzhou University, 48 Wenhui East Road, Yangzhou 225000, Jiangsu, China; Joint International Research Laboratory of Agriculture and Agri-product Safety, Yangzhou University, 48 Wenhui East Road, Yangzhou 225000, Jiangsu, China
| | - Yuanyue Tang
- Key Laboratory of Prevention and Control of Biological Hazard Factors (Animal Origin) for Agri-food Safety and Quality, Ministry of Agriculture of China, Yangzhou University, 48 Wenhui East Road, Yangzhou 225000, Jiangsu, China; Jiangsu Key Lab of Zoonosis/Jiangsu Co-innovation Center for Prevention and Control of Important Animal Infectious Diseases and Zoonoses, Yangzhou University, 48 Wenhui East Road, Yangzhou 225000, Jiangsu, China; Joint International Research Laboratory of Agriculture and Agri-product Safety, Yangzhou University, 48 Wenhui East Road, Yangzhou 225000, Jiangsu, China
| | - Xiaohui Zhou
- Department of Pathobiology and Veterinary Science, University of Connecticut, Storrs, CT 06269, USA; School of Public Health and Emergency Management, Southern University of Science and Technology, 1088 Xueyuan Road, Nanshan District, Shenzhen 518055, Guangdong, China
| | - Yunzeng Zhang
- Key Laboratory of Prevention and Control of Biological Hazard Factors (Animal Origin) for Agri-food Safety and Quality, Ministry of Agriculture of China, Yangzhou University, 48 Wenhui East Road, Yangzhou 225000, Jiangsu, China; Jiangsu Key Lab of Zoonosis/Jiangsu Co-innovation Center for Prevention and Control of Important Animal Infectious Diseases and Zoonoses, Yangzhou University, 48 Wenhui East Road, Yangzhou 225000, Jiangsu, China; Joint International Research Laboratory of Agriculture and Agri-product Safety, Yangzhou University, 48 Wenhui East Road, Yangzhou 225000, Jiangsu, China.
| | - Xinan Jiao
- Key Laboratory of Prevention and Control of Biological Hazard Factors (Animal Origin) for Agri-food Safety and Quality, Ministry of Agriculture of China, Yangzhou University, 48 Wenhui East Road, Yangzhou 225000, Jiangsu, China; Jiangsu Key Lab of Zoonosis/Jiangsu Co-innovation Center for Prevention and Control of Important Animal Infectious Diseases and Zoonoses, Yangzhou University, 48 Wenhui East Road, Yangzhou 225000, Jiangsu, China; Joint International Research Laboratory of Agriculture and Agri-product Safety, Yangzhou University, 48 Wenhui East Road, Yangzhou 225000, Jiangsu, China.
| | - Qiuchun Li
- Key Laboratory of Prevention and Control of Biological Hazard Factors (Animal Origin) for Agri-food Safety and Quality, Ministry of Agriculture of China, Yangzhou University, 48 Wenhui East Road, Yangzhou 225000, Jiangsu, China; Jiangsu Key Lab of Zoonosis/Jiangsu Co-innovation Center for Prevention and Control of Important Animal Infectious Diseases and Zoonoses, Yangzhou University, 48 Wenhui East Road, Yangzhou 225000, Jiangsu, China; Joint International Research Laboratory of Agriculture and Agri-product Safety, Yangzhou University, 48 Wenhui East Road, Yangzhou 225000, Jiangsu, China.
| |
Collapse
|
5
|
Long DR, Penewit K, Lo HY, Almazan J, Holmes EA, Bryan AB, Wolter DJ, Lewis JD, Waalkes A, Salipante SJ. In Vitro Selection Identifies Staphylococcus aureus Genes Influencing Biofilm Formation. Infect Immun 2023; 91:e0053822. [PMID: 36847490 PMCID: PMC10016075 DOI: 10.1128/iai.00538-22] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/27/2022] [Accepted: 02/13/2023] [Indexed: 03/01/2023] Open
Abstract
Staphylococcus aureus generates biofilms during many chronic human infections, which contributes to its growth and persistence in the host. Multiple genes and pathways necessary for S. aureus biofilm production have been identified, but knowledge is incomplete, and little is known about spontaneous mutations that increase biofilm formation as infection progresses. Here, we performed in vitro selection of four S. aureus laboratory strains (ATCC 29213, JE2, N315, and Newman) to identify mutations associated with enhanced biofilm production. Biofilm formation increased in passaged isolates from all strains, exhibiting from 1.2- to 5-fold the capacity of parental lines. Whole-genome sequencing identified nonsynonymous mutations affecting 23 candidate genes and a genomic duplication encompassing sigB. Six candidate genes significantly impacted biofilm formation as isogenic transposon knockouts: three were previously reported to impact S. aureus biofilm formation (icaR, spdC, and codY), while the remaining three (manA, narH, and fruB) were newly implicated by this study. Plasmid-mediated genetic complementation of manA, narH, and fruB transposon mutants corrected biofilm deficiencies, with high-level expression of manA and fruB further enhancing biofilm formation over basal levels. This work recognizes genes not previously identified as contributing to biofilm formation in S. aureus and reveals genetic changes able to augment biofilm production by that organism.
Collapse
Affiliation(s)
- Dustin R. Long
- Division of Critical Care Medicine, Department of Anesthesiology and Pain Medicine, University of Washington School of Medicine, Seattle, Washington, USA
| | - Kelsi Penewit
- Department of Laboratory Medicine and Pathology, University of Washington School of Medicine, Seattle, Washington, USA
| | - Hsin-Yu Lo
- Department of Laboratory Medicine and Pathology, University of Washington School of Medicine, Seattle, Washington, USA
| | - Jared Almazan
- Department of Laboratory Medicine and Pathology, University of Washington School of Medicine, Seattle, Washington, USA
| | - Elizabeth A. Holmes
- Department of Laboratory Medicine and Pathology, University of Washington School of Medicine, Seattle, Washington, USA
| | - Andrew B. Bryan
- Department of Laboratory Medicine and Pathology, University of Washington School of Medicine, Seattle, Washington, USA
| | - Daniel J. Wolter
- Department of Pediatrics, University of Washington School of Medicine, Seattle, Washington, USA
| | - Janessa D. Lewis
- Department of Laboratory Medicine and Pathology, University of Washington School of Medicine, Seattle, Washington, USA
| | - Adam Waalkes
- Department of Laboratory Medicine and Pathology, University of Washington School of Medicine, Seattle, Washington, USA
| | - Stephen J. Salipante
- Department of Laboratory Medicine and Pathology, University of Washington School of Medicine, Seattle, Washington, USA
| |
Collapse
|
6
|
Jiang J, Chen L, Chen X, Li P, Xu X, Fowler VG, van Duin D, Wang M. Carbapenemase-Encoding Gene Copy Number Estimator (CCNE): a Tool for Carbapenemase Gene Copy Number Estimation. Microbiol Spectr 2022; 10:e0100022. [PMID: 35863018 PMCID: PMC9431437 DOI: 10.1128/spectrum.01000-22] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2022] [Accepted: 06/14/2022] [Indexed: 11/20/2022] Open
Abstract
Carbapenemase production is one of the leading mechanisms of carbapenem resistance in Gram-negative bacteria. An increase in carbapenemase gene (blaCarb) copies is an important mechanism of carbapenem resistance. No currently available bioinformatics tools allow for reliable detection and reporting of carbapenemase gene copy numbers. Here, we describe the carbapenemase-encoding gene copy number estimator (CCNE), a ready-to-use bioinformatics tool that was developed to estimate blaCarb copy numbers from whole-genome sequencing data. Its performance on Klebsiella pneumoniae carbapenemase gene (blaKPC) copy number estimation was evaluated by simulation and quantitative PCR (qPCR), and the results were compared with available algorithms. CCNE has two components, CCNE-acc and CCNE-fast. CCNE-acc detects blaCarb copy number in a comprehensive and high-accuracy way, while CCNE-fast rapidly screens blaCarb copy numbers. CCNE-acc achieved the best accuracy (100%) and the lowest root mean squared error (RMSE; 0.07) in simulated noise data sets, compared to the assembly-based method (23.4% accuracy, 1.697 RMSE) and the OrthologsBased method (78.9% accuracy, 0.395 RMSE). In the qPCR validation, a high consistency was observed between the blaKPC copy number determined by qPCR and that determined with CCNE. Reverse transcription-qPCR transcriptional analysis of 40 isolates showed that blaKPC expression was positively correlated with the blaKPC copy numbers detected by CCNE (P < 0.001). An association study of 357 KPC-producing K. pneumoniae isolates and their antimicrobial susceptibility identified a significant association between the estimated blaKPC copy number and MICs of imipenem (P < 0.001) and ceftazidime-avibactam (P < 0.001). Overall, CCNE is a useful genomic tool for the analysis of antimicrobial resistance genes copy number; it is available at https://github.com/biojiang/ccne. IMPORTANCE Globally disseminated carbapenem-resistant Enterobacterales is an urgent threat to public health. The most common carbapenem resistance mechanism is the production of carbapenemases. Carbapenemase-producing isolates often exhibit a wide range of carbapenem MICs. Higher carbapenem MICs have been associated with treatment failure. The increase of carbapenemase gene (blaCarb) copy numbers contributes to increased carbapenem MICs. However, blaCarb gene copy number detection is not routinely conducted during a genomic analysis, in part due to the lack of optimal bioinformatics tools. In this study, we describe a ready-to-use tool we developed and designated the carbapenemase-encoding gene copy number estimator (CCNE) that can be used to estimate the blaCarb copy number directly from whole-genome sequencing data, and we extended the data to support the analysis of all known blaCarb genes and some other antimicrobial resistance genes. Furthermore, CCNE can be used to interrogate the correlations between genotypes and susceptibility phenotypes and to improve our understanding of antimicrobial resistance mechanisms.
Collapse
Affiliation(s)
- Jianping Jiang
- Institute of Antibiotics, Huashan Hospital, Fudan University, Shanghai, China
| | - Liang Chen
- Center for Discovery and Innovation, Hackensack Meridian Health, Nutley, New Jersey, USA
- Department of Medical Sciences, Hackensack Meridian School of Medicine, Nutley, New Jersey, USA
| | - Xin Chen
- Institute of Antibiotics, Huashan Hospital, Fudan University, Shanghai, China
| | - Pei Li
- Institute of Antibiotics, Huashan Hospital, Fudan University, Shanghai, China
| | - Xiaogang Xu
- Institute of Antibiotics, Huashan Hospital, Fudan University, Shanghai, China
| | - Vance G. Fowler
- Division of Infectious Diseases, University of North Carolina, Chapel Hill, North Carolina, USA
| | - David van Duin
- Duke Clinical Research Institute, Duke University Medical Center, Durham, North Carolina, USA
| | - Minggui Wang
- Institute of Antibiotics, Huashan Hospital, Fudan University, Shanghai, China
| |
Collapse
|
7
|
Jangir PK, Yang Q, Shaw LP, Caballero JD, Ogunlana L, Wheatley R, Walsh T, MacLean RC. Pre-existing chromosomal polymorphisms in pathogenic E. coli potentiate the evolution of resistance to a last-resort antibiotic. eLife 2022; 11:e78834. [PMID: 35943060 PMCID: PMC9363117 DOI: 10.7554/elife.78834] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2022] [Accepted: 07/22/2022] [Indexed: 12/17/2022] Open
Abstract
Bacterial pathogens show high levels of chromosomal genetic diversity, but the influence of this diversity on the evolution of antibiotic resistance by plasmid acquisition remains unclear. Here, we address this problem in the context of colistin, a 'last line of defence' antibiotic. Using experimental evolution, we show that a plasmid carrying the MCR-1 colistin resistance gene dramatically increases the ability of Escherichia coli to evolve high-level colistin resistance by acquiring mutations in lpxC, an essential chromosomal gene involved in lipopolysaccharide biosynthesis. Crucially, lpxC mutations increase colistin resistance in the presence of the MCR-1 gene, but decrease the resistance of wild-type cells, revealing positive sign epistasis for antibiotic resistance between the chromosomal mutations and a mobile resistance gene. Analysis of public genomic datasets shows that lpxC polymorphisms are common in pathogenic E. coli, including those carrying MCR-1, highlighting the clinical relevance of this interaction. Importantly, lpxC diversity is high in pathogenic E. coli from regions with no history of MCR-1 acquisition, suggesting that pre-existing lpxC polymorphisms potentiated the evolution of high-level colistin resistance by MCR-1 acquisition. More broadly, these findings highlight the importance of standing genetic variation and plasmid/chromosomal interactions in the evolutionary dynamics of antibiotic resistance.
Collapse
Affiliation(s)
- Pramod K Jangir
- Department of Zoology, University of OxfordOxfordUnited Kingdom
| | - Qiue Yang
- Fujian Provincial Key Laboratory of Soil Environmental Health and Regulation, College of Resources and Environment, Fujian Agriculture and Forestry UniversityFuzhouChina
| | - Liam P Shaw
- Department of Zoology, University of OxfordOxfordUnited Kingdom
| | | | - Lois Ogunlana
- Department of Zoology, University of OxfordOxfordUnited Kingdom
| | - Rachel Wheatley
- Department of Zoology, University of OxfordOxfordUnited Kingdom
| | - Timothy Walsh
- Department of Zoology, University of OxfordOxfordUnited Kingdom
| | - R Craig MacLean
- Department of Zoology, University of OxfordOxfordUnited Kingdom
| |
Collapse
|
8
|
Giulieri SG, Guérillot R, Duchene S, Hachani A, Daniel D, Seemann T, Davis JS, Tong SYC, Young BC, Wilson DJ, Stinear TP, Howden BP. Niche-specific genome degradation and convergent evolution shaping Staphylococcus aureus adaptation during severe infections. eLife 2022; 11:e77195. [PMID: 35699423 PMCID: PMC9270034 DOI: 10.7554/elife.77195] [Citation(s) in RCA: 29] [Impact Index Per Article: 9.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2022] [Accepted: 06/08/2022] [Indexed: 11/13/2022] Open
Abstract
During severe infections, Staphylococcus aureus moves from its colonising sites to blood and tissues and is exposed to new selective pressures, thus, potentially driving adaptive evolution. Previous studies have shown the key role of the agr locus in S. aureus pathoadaptation; however, a more comprehensive characterisation of genetic signatures of bacterial adaptation may enable prediction of clinical outcomes and reveal new targets for treatment and prevention of these infections. Here, we measured adaptation using within-host evolution analysis of 2590 S. aureus genomes from 396 independent episodes of infection. By capturing a comprehensive repertoire of single nucleotide and structural genome variations, we found evidence of a distinctive evolutionary pattern within the infecting populations compared to colonising bacteria. These invasive strains had up to 20-fold enrichments for genome degradation signatures and displayed significantly convergent mutations in a distinctive set of genes, linked to antibiotic response and pathogenesis. In addition to agr-mediated adaptation, we identified non-canonical, genome-wide significant loci including sucA-sucB and stp1. The prevalence of adaptive changes increased with infection extent, emphasising the clinical significance of these signatures. These findings provide a high-resolution picture of the molecular changes when S. aureus transitions from colonisation to severe infection and may inform correlation of infection outcomes with adaptation signatures.
Collapse
Affiliation(s)
- Stefano G Giulieri
- Department of Microbiology and Immunology at the Peter Doherty Institute for Infection and Immunity, University of MelbourneMelbourneAustralia
- Department of Infectious Diseases, Austin HealthHeidelbergAustralia
- Victorian Infectious Diseases Service, Royal Melbourne HospitalMelbourneAustralia
| | - Romain Guérillot
- Department of Microbiology and Immunology at the Peter Doherty Institute for Infection and Immunity, University of MelbourneMelbourneAustralia
| | - Sebastian Duchene
- Department of Microbiology and Immunology at the Peter Doherty Institute for Infection and Immunity, University of MelbourneMelbourneAustralia
| | - Abderrahman Hachani
- Department of Microbiology and Immunology at the Peter Doherty Institute for Infection and Immunity, University of MelbourneMelbourneAustralia
| | - Diane Daniel
- Department of Microbiology and Immunology at the Peter Doherty Institute for Infection and Immunity, University of MelbourneMelbourneAustralia
- Microbiological Diagnostic Unit Public Health Laboratory, The University of Melbourne at the Doherty Institute for Infection and ImmunityMelbourneAustralia
| | - Torsten Seemann
- Microbiological Diagnostic Unit Public Health Laboratory, The University of Melbourne at the Doherty Institute for Infection and ImmunityMelbourneAustralia
| | - Joshua S Davis
- Department of Infectious Diseases, John Hunter HospitalNewcastle, New South WalesAustralia
- Menzies School of Health Research, Charles Darwin UniversityCasuarina, Northern TerritoryAustralia
| | - Steven YC Tong
- Menzies School of Health Research, Charles Darwin UniversityCasuarina, Northern TerritoryAustralia
- Victorian Infectious Disease Service, Royal Melbourne Hospital, and University of Melbourne at the Peter Doherty Institute for Infection and ImmunityMelbourneAustralia
| | | | - Daniel J Wilson
- Big Data Institute, Nuffield Department of Population Health, Li Ka Shing Centre for Health Information and Discovery, Old Road Campus, University of OxfordOxfordUnited Kingdom
| | - Timothy P Stinear
- Department of Microbiology and Immunology at the Peter Doherty Institute for Infection and Immunity, University of MelbourneMelbourneAustralia
| | - Benjamin P Howden
- Department of Microbiology and Immunology at the Peter Doherty Institute for Infection and Immunity, University of MelbourneMelbourneAustralia
- Department of Infectious Diseases, Austin HealthHeidelbergAustralia
- Microbiological Diagnostic Unit Public Health Laboratory, The University of Melbourne at the Doherty Institute for Infection and ImmunityMelbourneAustralia
| |
Collapse
|
9
|
Kapel N, Caballero JD, MacLean RC. Localized pmrB hypermutation drives the evolution of colistin heteroresistance. Cell Rep 2022; 39:110929. [PMID: 35675785 PMCID: PMC9189680 DOI: 10.1016/j.celrep.2022.110929] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2021] [Revised: 12/02/2021] [Accepted: 05/16/2022] [Indexed: 11/16/2022] Open
Abstract
Colistin has emerged as an important last line of defense for the treatment of infections caused by antibiotic-resistant gram-negative pathogens, but colistin resistance remains poorly understood. Here, we investigate the responses of ≈1,000 populations of a multi-drug-resistant (MDR) strain of P. aeruginosa to a high dose of colistin. Colistin exposure causes rapid cell death, but some populations eventually recover due to the growth of sub-populations of heteroresistant cells. Heteroresistance is unstable, and resistance is rapidly lost under culture in colistin-free medium. The evolution of heteroresistance is primarily driven by selection for heteroresistance at two hotspot sites in the PmrAB regulatory system. Localized hypermutation of pmrB generates colistin resistance at 103–104 times the background resistance mutation rate (≈2 × 10-5 per cell division). PmrAB provides resistance to antimicrobial peptides that are involved in host immunity, suggesting that this pathogen may have evolved a highly mutable pmrB as an adaptation to host immunity. Pseudomonas populations recover from colistin due to the growth of heteroresistant cells Heteroresistance is driven by pre-existing mutations in the PmrAB regulatory system pmrB mutations arise at 103–104 times the background mutation rate Heteroresistance is unstable and is rapidly lost in the absence of colistin
Collapse
Affiliation(s)
- Natalia Kapel
- University of Oxford, Department of Zoology, 11a Mansfield Road, Oxford OX1 3SZ, UK
| | - Julio Diaz Caballero
- University of Oxford, Department of Zoology, 11a Mansfield Road, Oxford OX1 3SZ, UK
| | - R Craig MacLean
- University of Oxford, Department of Zoology, 11a Mansfield Road, Oxford OX1 3SZ, UK.
| |
Collapse
|
10
|
Werth BJ, Ashford NK, Penewit K, Waalkes A, Holmes EA, Bryan A, Salipante SJ. OUP accepted manuscript. JAC Antimicrob Resist 2022; 4:dlac011. [PMID: 35156034 PMCID: PMC8827560 DOI: 10.1093/jacamr/dlac011] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2021] [Accepted: 01/19/2022] [Indexed: 11/12/2022] Open
Abstract
Objectives Cefiderocol is a siderophore cephalosporin active against MDR Gram-negatives including Stenotrophomonas maltophilia. Cefiderocol resistance remains uncommon and incompletely understood. We selected for cefiderocol-resistant S. maltophilia in vitro and characterized the genetic mechanisms and potential for cross-resistance to other antimicrobials. Methods We selected cefiderocol resistance in three clinical strains of S. maltophilia by serial passage in escalating concentrations of cefiderocol. Emergent cefiderocol-resistant isolates were subjected to repeat susceptibility testing against a panel of relevant antimicrobials. Isolates with confirmed MIC changes were whole genome sequenced. Results Each parent strain was initially susceptible to cefiderocol (MICs of 0.03125, 0.03125 and 0.125 mg/L), and one initially tested susceptible to ceftazidime/avibactam (MIC 4 mg/L). We recovered evolved isolates achieving cefiderocol resistance at MICs of 8–32 mg/L from each parental strain. Some cefiderocol resistant isolates reverted following one to four drug-free passages. Ceftazidime/avibactam MICs of passaged isolates repeatedly increased to ≥256 mg/L, and while other MICs were largely unchanged, trimethoprim/sulfamethoxazole MICs declined 4-fold in two strains. WGS revealed one evolved isolate carrying six coding mutations, while four were isogenic mutants of tonB, tolQ, smf-1 and the smeT promoter. Mutation of the smeT promoter downregulated the smeDEF efflux pump and reduced susceptibility to penicillins but increased susceptibility to several other classes including sulphonamides. Other mutations occurred in genes putatively involved in iron metabolism including smlt1148 and cirA. Conclusions S. maltophilia strains evolved cefiderocol resistance through different genetic pathways, but often involved iron transport. Future work is required to fully understand the role(s) of other genes in cefiderocol resistance.
Collapse
Affiliation(s)
- Brian J. Werth
- Department of Pharmacy, School of Pharmacy, University of Washington, Seattle, WA, USA
- Corresponding author. E-mail:
| | - Nathaniel K. Ashford
- Department of Pharmacy, School of Pharmacy, University of Washington, Seattle, WA, USA
| | - Kelsi Penewit
- Department of Laboratory Medicine and Pathology, School of Medicine, University of Washington, Seattle, WA, USA
| | - Adam Waalkes
- Department of Laboratory Medicine and Pathology, School of Medicine, University of Washington, Seattle, WA, USA
| | - Elizabeth A. Holmes
- Department of Laboratory Medicine and Pathology, School of Medicine, University of Washington, Seattle, WA, USA
| | - Andrew Bryan
- Department of Laboratory Medicine and Pathology, School of Medicine, University of Washington, Seattle, WA, USA
| | - Stephen J. Salipante
- Department of Laboratory Medicine and Pathology, School of Medicine, University of Washington, Seattle, WA, USA
| |
Collapse
|
11
|
Sicard A, Saponari M, Vanhove M, Castillo AI, Giampetruzzi A, Loconsole G, Saldarelli P, Boscia D, Neema C, Almeida RPP. Introduction and adaptation of an emerging pathogen to olive trees in Italy. Microb Genom 2021; 7. [PMID: 34904938 PMCID: PMC8767334 DOI: 10.1099/mgen.0.000735] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
The invasive plant pathogen Xylella fastidiosa currently threatens European flora through the loss of economically and culturally important host plants. This emerging vector-borne bacterium, native to the Americas, causes several important diseases in a wide range of plants including crops, ornamentals, and trees. Previously absent from Europe, and considered a quarantine pathogen, X. fastidiosa was first detected in Apulia, Italy in 2013 associated with a devastating disease of olive trees (Olive Quick Decline Syndrome, OQDS). OQDS has led to significant economic, environmental, cultural, as well as political crises. Although the biology of X. fastidiosa diseases have been studied for over a century, there is still no information on the determinants of specificity between bacterial genotypes and host plant species, which is particularly relevant today as X. fastidiosa is expanding in the naive European landscape. We analysed the genomes of 79 X. fastidiosa samples from diseased olive trees across the affected area in Italy as well as genomes of the most genetically closely related strains from Central America. We provided insights into the ecological and evolutionary emergence of this pathogen in Italy. We first showed that the outbreak in Apulia is due to a single introduction from Central America that we estimated to have occurred in 2008 [95 % HPD: 1930–2016]. By using a combination of population genomic approaches and evolutionary genomics methods, we further identified a short list of genes that could play a major role in the adaptation of X. fastidiosa to this new environment. We finally provided experimental evidence for the adaptation of the strain to this new environment.
Collapse
Affiliation(s)
- Anne Sicard
- UC Berkeley, Department of Environmental Science, Policy, and Management, Berkeley, CA 94720, U.S.A.,PHIM Plant Health Institute, Univ Montpellier, INRAE, Institut Agro, CIRAD, IRD, Montpellier, France
| | - Maria Saponari
- National Research Council (CNR), Institute for Sustainable Plant Protection, Via Amendola 122/D, 70126 Bari, Italy
| | - Mathieu Vanhove
- UC Berkeley, Department of Environmental Science, Policy, and Management, Berkeley, CA 94720, U.S.A
| | - Andreina I Castillo
- UC Berkeley, Department of Environmental Science, Policy, and Management, Berkeley, CA 94720, U.S.A
| | - Annalisa Giampetruzzi
- University of Bari Aldo Moro, Department of Soil, Plant and Food Sciences, Piazza Umberto I, 70121 Bari, Italy
| | - Giuliana Loconsole
- National Research Council (CNR), Institute for Sustainable Plant Protection, Via Amendola 122/D, 70126 Bari, Italy
| | - Pasquale Saldarelli
- National Research Council (CNR), Institute for Sustainable Plant Protection, Via Amendola 122/D, 70126 Bari, Italy
| | - Donato Boscia
- National Research Council (CNR), Institute for Sustainable Plant Protection, Via Amendola 122/D, 70126 Bari, Italy
| | - Claire Neema
- PHIM Plant Health Institute, Univ Montpellier, INRAE, Institut Agro, CIRAD, IRD, Montpellier, France
| | - Rodrigo P P Almeida
- UC Berkeley, Department of Environmental Science, Policy, and Management, Berkeley, CA 94720, U.S.A
| |
Collapse
|
12
|
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: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/04/2021] [Accepted: 10/01/2021] [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
| |
Collapse
|
13
|
Jugas R, Sedlar K, Vitek M, Nykrynova M, Barton V, Bezdicek M, Lengerova M, Skutkova H. CNproScan: Hybrid CNV detection for bacterial genomes. Genomics 2021; 113:3103-3111. [PMID: 34224809 DOI: 10.1016/j.ygeno.2021.06.040] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2020] [Revised: 06/13/2021] [Accepted: 06/30/2021] [Indexed: 10/20/2022]
Abstract
Discovering copy number variation (CNV) in bacteria is not in the spotlight compared to the attention focused on CNV detection in eukaryotes. However, challenges arising from bacterial drug resistance bring further interest to the topic of CNV and its role in drug resistance. General CNV detection methods do not consider bacteria's features and there is space to improve detection accuracy. Here, we present a CNV detection method called CNproScan focused on bacterial genomes. CNproScan implements a hybrid approach and other bacteria-focused features and depends only on NGS data. We benchmarked our method and compared it to the previously published methods and we can resolve to achieve a higher detection rate together with providing other beneficial features, such as CNV classification. Compared with other methods, CNproScan can detect much shorter CNV events.
Collapse
Affiliation(s)
- Robin Jugas
- Department of Biomedical Engineering, Brno University of Technology, Brno, Czech Republic.
| | - Karel Sedlar
- Department of Biomedical Engineering, Brno University of Technology, Brno, Czech Republic
| | - Martin Vitek
- Department of Biomedical Engineering, Brno University of Technology, Brno, Czech Republic
| | - Marketa Nykrynova
- Department of Biomedical Engineering, Brno University of Technology, Brno, Czech Republic
| | - Vojtech Barton
- Department of Biomedical Engineering, Brno University of Technology, Brno, Czech Republic
| | - Matej Bezdicek
- Department of Internal Medicine-Hematology and Oncology, University Hospital Brno, Brno, Czech Republic
| | - Martina Lengerova
- Department of Internal Medicine-Hematology and Oncology, University Hospital Brno, Brno, Czech Republic
| | - Helena Skutkova
- Department of Biomedical Engineering, Brno University of Technology, Brno, Czech Republic
| |
Collapse
|
14
|
Shared and Unique Evolutionary Trajectories to Ciprofloxacin Resistance in Gram-Negative Bacterial Pathogens. mBio 2021; 12:e0098721. [PMID: 34154405 PMCID: PMC8262867 DOI: 10.1128/mbio.00987-21] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
Abstract
Resistance to the broad-spectrum antibiotic ciprofloxacin is detected at high rates for a wide range of bacterial pathogens. To investigate the dynamics of ciprofloxacin resistance development, we applied a comparative resistomics workflow for three clinically relevant species of Gram-negative bacteria: Escherichia coli, Acinetobacter baumannii, and Pseudomonas aeruginosa. We combined experimental evolution in a morbidostat with deep sequencing of evolving bacterial populations in time series to reveal both shared and unique aspects of evolutionary trajectories. Representative clone characterization by sequencing and MIC measurements enabled direct assessment of the impact of mutations on the extent of acquired drug resistance. In all three species, we observed a two-stage evolution: (i) early ciprofloxacin resistance reaching 4- to 16-fold the MIC for the wild type, commonly as a result of single mutations in DNA gyrase target genes (gyrA or gyrB), and (ii) additional genetic alterations affecting the transcriptional control of the drug efflux machinery or secondary target genes (DNA topoisomerase parC or parE).
Collapse
|
15
|
Shi L, Qin J, Zheng H, Guo Y, Zhang H, Zhong Y, Yang C, Dong S, Yang F, Wu Y, Zhao G, Song Y, Yang R, Wang P, Cui Y. New Genotype of Yersinia pestis Found in Live Rodents in Yunnan Province, China. Front Microbiol 2021; 12:628335. [PMID: 33935990 PMCID: PMC8084289 DOI: 10.3389/fmicb.2021.628335] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2020] [Accepted: 03/24/2021] [Indexed: 12/11/2022] Open
Abstract
Yunnan Province, China is thought to be the original source of biovar Orientalis of Yersinia pestis, the causative agent of the third plague pandemic that has spread globally since the end of the 19th century. Although encompassing a large area of natural plague foci, Y. pestis strains have rarely been found in live rodents during surveillance in Yunnan, and most isolates are from rodent corpses and their fleas. In 2017, 10 Y. pestis strains were isolated from seven live rodents and three fleas in Heqing County of Yunnan. These strains were supposed to have low virulence to local rodents Eothenomys miletus and Apodemus chevrieri because the rodents were healthy and no dead animals were found in surrounding areas, as had occurred in previous epizootic disease. We performed microscopic and biochemical examinations of the isolates, and compared their whole-genome sequences and transcriptome with those of 10 high virulence Y. pestis strains that were isolated from nine rodents and one parasitic flea in adjacent city (Lijiang). We analyzed the phenotypic, genomic, and transcriptomic characteristics of live rodent isolates. The isolates formed a previously undefined monophyletic branch of Y. pestis that was named 1.IN5. Six SNPs, two indels, and one copy number variation were detected between live rodent isolates and the high virulence neighbors. No obvious functional consequence of these variations was found according to the known annotation information. Among genes which expression differential in the live rodent isolates compared to their high virulent neighbors, we found five iron transfer related ones that were significant up-regulated (| log2 (FC) | > 1, p.adjust < 0.05), indicating these genes may be related to the low-virulence phenotype. The novel genotype of Y. pestis reported here provides further insights into the evolution and spread of plague as well as clues that may help to decipher the virulence mechanism of this notorious pathogen.
Collapse
Affiliation(s)
- Liyuan Shi
- Yunnan Institute of Endemic Diseases Control and Prevention, Dali, China
| | - Jingliang Qin
- State Key Laboratory of Pathogen and Biosecurity, Beijing Institute of Microbiology and Epidemiology, Beijing, China.,School of Basic Medical Sciences, Anhui Medical University, Hefei, China
| | - Hongyuan Zheng
- State Key Laboratory of Pathogen and Biosecurity, Beijing Institute of Microbiology and Epidemiology, Beijing, China
| | - Ying Guo
- Yunnan Institute of Endemic Diseases Control and Prevention, Dali, China
| | - Haipeng Zhang
- Yunnan Institute of Endemic Diseases Control and Prevention, Dali, China
| | - Youhong Zhong
- Yunnan Institute of Endemic Diseases Control and Prevention, Dali, China
| | - Chao Yang
- State Key Laboratory of Pathogen and Biosecurity, Beijing Institute of Microbiology and Epidemiology, Beijing, China
| | - Shanshan Dong
- Yunnan Institute of Endemic Diseases Control and Prevention, Dali, China
| | - Fengyi Yang
- Yunnan Institute of Endemic Diseases Control and Prevention, Dali, China
| | - Yarong Wu
- State Key Laboratory of Pathogen and Biosecurity, Beijing Institute of Microbiology and Epidemiology, Beijing, China
| | - Guangyu Zhao
- State Key Laboratory of Pathogen and Biosecurity, Beijing Institute of Microbiology and Epidemiology, Beijing, China.,School of Basic Medical Sciences, Anhui Medical University, Hefei, China
| | - Yajun Song
- State Key Laboratory of Pathogen and Biosecurity, Beijing Institute of Microbiology and Epidemiology, Beijing, China
| | - Ruifu Yang
- State Key Laboratory of Pathogen and Biosecurity, Beijing Institute of Microbiology and Epidemiology, Beijing, China
| | - Peng Wang
- Yunnan Institute of Endemic Diseases Control and Prevention, Dali, China
| | - Yujun Cui
- State Key Laboratory of Pathogen and Biosecurity, Beijing Institute of Microbiology and Epidemiology, Beijing, China.,School of Basic Medical Sciences, Anhui Medical University, Hefei, China
| |
Collapse
|
16
|
Roe C, Stegger M, Lilje B, Johannesen TB, Ng KL, Sieber RN, Driebe E, Engelthaler DM, Andersen PS. Genomic analyses of Staphylococcus aureus clonal complex 45 isolates does not distinguish nasal carriage from bacteraemia. Microb Genom 2020; 6:mgen000403. [PMID: 32667872 PMCID: PMC7641415 DOI: 10.1099/mgen.0.000403] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2019] [Accepted: 06/16/2020] [Indexed: 11/20/2022] Open
Abstract
Staphylococcus aureus is a colonizing opportunistic pathogen and a leading cause of bloodstream infection with high morbidity and mortality. S. aureus carriage frequency is reportedly between 20 and 40 % among healthy adults, with S. aureus colonization considered to be a risk factor for S. aureus bacteraemia. It is unknown whether a genetic component of the bacterium is associated with S. aureus bacteraemia in comparison to nasal carriage strains. Previous association studies primarily focusing on the clinical outcome of an S. aureus infection have produced conflicting results, often limited by study design challenged by sample collections and the clonal diversity of S. aureus. To date, no study has investigated whether genomic features separate nasal carriage isolates from S. aureus bacteraemia isolates within a single clonal lineage. Here we have investigated whether genomic features, including single-nucleotide polymorphisms (SNPs), genes, or kmers, distinguish S. aureus nasal carriage isolates from bacteraemia isolates that all belong to the same clonal lineage [clonal complex 45 (CC45)] using whole-genome sequencing (WGS) and a genome-wide association (GWA) approach. From CC45, 100 isolates (50 bacteraemia and 50 nasal carriage, geographically and temporally matched) from Denmark were whole-genome sequenced and subjected to GWA analyses involving gene copy number variation, SNPs, gene content, kmers and gene combinations, while correcting for lineage effects. No statistically significant association involving SNPs, specific genes, gene variants, gene copy number variation, or a combination of genes was identified that could distinguish bacteraemia isolates from nasal carriage isolates. The presented results suggest that all S. aureus nasal CC45 isolates carry the potential to cause invasive disease, as no core or accessory genome content or variations were statistically associated with invasiveness.
Collapse
Affiliation(s)
- Chandler Roe
- Translational Genomics Research Institute, Flagstaff, AZ, USA
- Northern Arizona University, Flagstaff, AZ, USA
| | - Marc Stegger
- Bacteria, Parasites and Fungi, Statens Serum Institut, Copenhagen, Denmark
| | - Berit Lilje
- Bacteria, Parasites and Fungi, Statens Serum Institut, Copenhagen, Denmark
| | | | - Kim Lee Ng
- Bacteria, Parasites and Fungi, Statens Serum Institut, Copenhagen, Denmark
| | - Raphael N. Sieber
- Bacteria, Parasites and Fungi, Statens Serum Institut, Copenhagen, Denmark
| | | | | | - Paal Skytt Andersen
- Translational Genomics Research Institute, Flagstaff, AZ, USA
- Bacteria, Parasites and Fungi, Statens Serum Institut, Copenhagen, Denmark
- Department of Veterinary and Animal Sciences, University of Copenhagen, Frederiksberg, Denmark
| |
Collapse
|
17
|
Wu L, Wang H, Xia Y, Xi R. CNV-BAC: Copy number Variation Detection in Bacterial Circular Genome. Bioinformatics 2020; 36:3890-3891. [DOI: 10.1093/bioinformatics/btaa208] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2019] [Revised: 03/14/2020] [Accepted: 03/23/2020] [Indexed: 11/14/2022] Open
Abstract
Abstract
Motivation
Whole-genome sequencing (WGS) is widely used for copy number variation (CNV) detection. However, for most bacteria, their circular genome structure and high replication rate make reads more enriched near the replication origin. CNV detection based on read depth could be seriously influenced by such replication bias.
Results
We show that the replication bias is widespread using ∼200 bacterial WGS data. We develop CNV-BAC (CNV-Bacteria) that can properly normalize the replication bias and other known biases in bacterial WGS data and can accurately detect CNVs. Simulation and real data analysis show that CNV-BAC achieves the best performance in CNV detection compared with available algorithms.
Availability and implementation
CNV-BAC is available at https://github.com/XiDsLab/CNV-BAC.
Supplementary information
Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
| | - Han Wang
- School of Mathematical Sciences
- Center for Statistical Science, Peking University, Beijing 100871, China
| | - Yuchao Xia
- School of Applied Science, Beijing Information Science and Technology University, Beijing 100101, China
| | - Ruibin Xi
- School of Mathematical Sciences
- Center for Statistical Science, Peking University, Beijing 100871, China
- Department of Biostatistics, Peking University, Beijing 100871, China
| |
Collapse
|
18
|
San JE, Baichoo S, Kanzi A, Moosa Y, Lessells R, Fonseca V, Mogaka J, Power R, de Oliveira T. Current Affairs of Microbial Genome-Wide Association Studies: Approaches, Bottlenecks and Analytical Pitfalls. Front Microbiol 2020; 10:3119. [PMID: 32082269 PMCID: PMC7002396 DOI: 10.3389/fmicb.2019.03119] [Citation(s) in RCA: 46] [Impact Index Per Article: 9.2] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2019] [Accepted: 12/24/2019] [Indexed: 12/12/2022] Open
Abstract
Microbial genome-wide association studies (mGWAS) are a new and exciting research field that is adapting human GWAS methods to understand how variations in microbial genomes affect host or pathogen phenotypes, such as drug resistance, virulence, host specificity and prognosis. Several computational tools and methods have been developed or adapted from human GWAS to facilitate the discovery of novel mutations and structural variations that are associated with the phenotypes of interest. However, no comprehensive, end-to-end, user-friendly tool is currently available. The development of a broadly applicable pipeline presents a real opportunity among computational biologists. Here, (i) we review the prominent and promising tools, (ii) discuss analytical pitfalls and bottlenecks in mGWAS, (iii) provide insights into the selection of appropriate tools, (iv) highlight the gaps that still need to be filled and how users and developers can work together to overcome these bottlenecks. Use of mGWAS research can inform drug repositioning decisions as well as accelerate the discovery and development of more effective vaccines and antimicrobials for pressing infectious diseases of global health significance, such as HIV, TB, influenza, and malaria.
Collapse
Affiliation(s)
- James Emmanuel San
- Kwazulu-Natal Research and Innovation Sequencing Platform (KRISP), College of Health Sciences, University of KwaZulu-Natal, Durban, South Africa
| | - Shakuntala Baichoo
- Department of Digital Technologies, FoICDT, University of Mauritius, Réduit, Mauritius
| | - Aquillah Kanzi
- Kwazulu-Natal Research and Innovation Sequencing Platform (KRISP), College of Health Sciences, University of KwaZulu-Natal, Durban, South Africa
| | - Yumna Moosa
- Kwazulu-Natal Research and Innovation Sequencing Platform (KRISP), College of Health Sciences, University of KwaZulu-Natal, Durban, South Africa
| | - Richard Lessells
- Kwazulu-Natal Research and Innovation Sequencing Platform (KRISP), College of Health Sciences, University of KwaZulu-Natal, Durban, South Africa
| | - Vagner Fonseca
- Kwazulu-Natal Research and Innovation Sequencing Platform (KRISP), College of Health Sciences, University of KwaZulu-Natal, Durban, South Africa
- Laboratório de Genética Celular e Molecular, ICB, Universidade Federal de Minas Gerais, Belo Horizonte, Brazil
| | - John Mogaka
- Discipline of Public Health, University of Kwazulu-Natal, Durban, South Africa
| | - Robert Power
- St Edmund Hall, Oxford University, Oxford, United Kingdom
| | - Tulio de Oliveira
- Kwazulu-Natal Research and Innovation Sequencing Platform (KRISP), College of Health Sciences, University of KwaZulu-Natal, Durban, South Africa
- Department of Global Health, University of Washington, Seattle, WA, United States
| |
Collapse
|
19
|
The Copy Number of the spoVA 2mob Operon Determines Pressure Resistance of Bacillus Endospores. Appl Environ Microbiol 2019; 85:AEM.01596-19. [PMID: 31375487 DOI: 10.1128/aem.01596-19] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2019] [Accepted: 07/27/2019] [Indexed: 12/27/2022] Open
Abstract
The spoVA 2mob operon confers heat resistance to Bacillus spp., and the resistance correlates to the copy number of the operon. Bacillus endospores also exhibit a strong variation in resistance to pressure, but the underlying mechanisms of endospore resistance to pressure are not fully understood. We determined the effects of multiple spoVA 2mob operons on high-pressure resistance in Bacillus endospores. The copy numbers of the spoVA 2mob operon in 17 strains of Bacillus subtilis, Bacillus amyloliquefaciens, Bacillus cereus, Bacillus velezensis, and Bacillus pumilus were determined via droplet digital PCR (ddPCR) and genome sequencing. These strains contained between 0 and 3 copies of the spoVA 2mob operon; the quantification of the gene copy number by ddPCR was as accurate as whole-genome sequencing. We further tested the pressure resistance of 17 Bacillus endospores at 600 MPa and 80°C. Strains with one or no spoVA 2mob operon had significantly lower pressure resistance than strains with two or three copies of the operons (P < 0.001), indicating that redundant spoVA 2mob operons in Bacillus contributed to higher pressure resistance of endospores. The copy number of the spoVA 2mob operon was not related to the dipicolinic acid (DPA) content of endospores. Overall, the copy number of the spoVA 2mob operon contributes to pressure resistance of Bacillus endospores. This improves our understanding of the pressure resistance mechanisms in Bacillus spp. and may inform the development of high-pressure sterilization in food processing.IMPORTANCE Bacillus spp. are considered pressure-resistant microorganisms, but the resistance mechanisms remain unknown. The spoVA 2mob operon is a mobile genetic element, and it can transfer to pathogenic or spoilage organisms by horizontal gene transfer. Results in this study indicate that multiple copies of the spoVA 2mob operon mediate high-pressure resistance of Bacillus endospores, and it might contribute to the identification of the source of pressure-resistant pathogens and spoilage organisms that may contaminate the food supply. The droplet digital PCR (ddPCR) system is well suited for analysis in some human diseases due to its high efficiency and capability to provide high precision; however, no relevant studies in food microbiology have been reported so far. This study demonstrates a novel application of ddPCR in food microbiology.
Collapse
|
20
|
Wardell SJT, Rehman A, Martin LW, Winstanley C, Patrick WM, Lamont IL. A large-scale whole-genome comparison shows that experimental evolution in response to antibiotics predicts changes in naturally evolved clinical Pseudomonas aeruginosa. Antimicrob Agents Chemother 2019; 63:AAC.01619-19. [PMID: 31570397 PMCID: PMC6879238 DOI: 10.1128/aac.01619-19] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2019] [Accepted: 09/23/2019] [Indexed: 12/13/2022] Open
Abstract
Pseudomonas aeruginosa is an opportunistic pathogen that causes a wide range of acute and chronic infections. An increasing number of isolates have mutations that make them antibiotic resistant, making treatment difficult. To identify resistance-associated mutations we experimentally evolved the antibiotic sensitive strain P. aeruginosa PAO1 to become resistant to three widely used anti-pseudomonal antibiotics, ciprofloxacin, meropenem and tobramycin. Mutants could tolerate up to 2048-fold higher concentrations of antibiotic than strain PAO1. Genome sequences were determined for thirteen mutants for each antibiotic. Each mutant had between 2 and 8 mutations. For each antibiotic at least 8 genes were mutated in multiple mutants, demonstrating the genetic complexity of resistance. For all three antibiotics mutations arose in genes known to be associated with resistance, but also in genes not previously associated with resistance. To determine the clinical relevance of mutations uncovered in this study we analysed the corresponding genes in 558 isolates of P. aeruginosa from patients with chronic lung disease and in 172 isolates from the general environment. Many genes identified through experimental evolution had predicted function-altering changes in clinical isolates but not in environmental isolates, showing that mutated genes in experimentally evolved bacteria can predict those that undergo mutation during infection. Additionally, large deletions of up to 479kb arose in experimentally evolved meropenem resistant mutants and large deletions were present in 87 of the clinical isolates. These findings significantly advance understanding of antibiotic resistance in P. aeruginosa and demonstrate the validity of experimental evolution in identifying clinically-relevant resistance-associated mutations.
Collapse
Affiliation(s)
| | - Attika Rehman
- Department of Biochemistry, University of Otago, Dunedin, New Zealand
| | - Lois W Martin
- Department of Biochemistry, University of Otago, Dunedin, New Zealand
| | - Craig Winstanley
- Institute of Infection and Global Health, University of Liverpool, Liverpool, United Kingdom
| | - Wayne M Patrick
- Department of Biochemistry, University of Otago, Dunedin, New Zealand
- School of Biological Sciences, Victoria University of Wellington, Wellington, New Zealand
| | - Iain L Lamont
- Department of Biochemistry, University of Otago, Dunedin, New Zealand
| |
Collapse
|
21
|
Abstract
Pullorum disease, an acute poultry septicemia caused by Salmonella Gallinarum biovar Pullorum, is fatal for young chickens and is a heavy burden on poultry industry. The pathogen is rare in most developed countries but still extremely difficult to eliminate in China. Efficient epidemiological surveillance necessitates clarifying the origin of the isolates from different regions and their phylogenic relationships. Genomic epidemiological analysis of 97 S. Pullorum strains was carried out to reconstruct the phylogeny and transmission history of S. Pullorum. Further analysis demonstrated that functional gene loss and acquisition occurred simultaneously throughout the evolution of S. Pullorum, both of which reflected adaptation to the changing environment. The result of our study will be helpful in surveillance and prevention of pullorum disease. Salmonella enterica subspecies enterica serovar Gallinarum biovar Pullorum (S. Pullorum) is the etiological agent of pullorum disease, causing white diarrhea with high mortality in chickens. There are many unsolved issues surrounding the epidemiology of S. Pullorum, including its origin and transmission history as well as the discordance between its phenotypic heterogeneity and genetic monomorphism. In this paper, we report the results of whole-genome sequencing of a panel of 97 S. Pullorum strains isolated between 1962 and 2014 from four countries across three continents. We utilized 6,795 core genome single nucleotide polymorphisms (SNPs) to reconstruct a phylogenetic tree within a spatiotemporal Bayesian framework, estimating that the most recent common ancestor of S. Pullorum emerged in ∼914 CE (95% confidence interval [95%CI], 565 to 1273 CE). The extant S. Pullorum strains can be divided into four distinct lineages, each of which is significantly associated with geographical distribution. The intercontinental transmissions of lineages III and IV can be traced to the mid-19th century and are probably related to the “Hen Fever” prevalent at that time. Further genomic analysis indicated that the loss or pseudogenization of functional genes involved in metabolism and virulence in S. Pullorum has been ongoing since before and after divergence from the ancestor. In contrast, multiple prophages and plasmids have been acquired by S. Pullorum, and these have endowed it with new characteristics, especially the multidrug resistance conferred by two large plasmids in lineage I. The results of this study provide insight into the evolution of S. Pullorum and prove the efficiency of whole-genome sequencing in epidemiological surveillance of pullorum disease. IMPORTANCE Pullorum disease, an acute poultry septicemia caused by Salmonella Gallinarum biovar Pullorum, is fatal for young chickens and is a heavy burden on poultry industry. The pathogen is rare in most developed countries but still extremely difficult to eliminate in China. Efficient epidemiological surveillance necessitates clarifying the origin of the isolates from different regions and their phylogenic relationships. Genomic epidemiological analysis of 97 S. Pullorum strains was carried out to reconstruct the phylogeny and transmission history of S. Pullorum. Further analysis demonstrated that functional gene loss and acquisition occurred simultaneously throughout the evolution of S. Pullorum, both of which reflected adaptation to the changing environment. The result of our study will be helpful in surveillance and prevention of pullorum disease.
Collapse
|
22
|
Desvillechabrol D, Bouchier C, Kennedy S, Cokelaer T. Sequana coverage: detection and characterization of genomic variations using running median and mixture models. Gigascience 2018; 7:5091804. [PMID: 30192951 PMCID: PMC6275460 DOI: 10.1093/gigascience/giy110] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/21/2017] [Accepted: 08/23/2018] [Indexed: 11/30/2022] Open
Abstract
Background In addition to mapping quality information, the Genome coverage contains valuable biological information such as the presence of repetitive regions, deleted genes, or copy number variations (CNVs). It is essential to take into consideration atypical regions, trends (e.g., origin of replication), or known and unknown biases that influence coverage. It is also important that reported events have robust statistics (e.g. z-score) associated with their detections as well as precise location. Results We provide a stand-alone application, sequana_coverage, that reports genomic regions of interest (ROIs) that are significantly over- or underrepresented in high-throughput sequencing data. Significance is associated with the events as well as characteristics such as length of the regions. The algorithm first detrends the data using an efficient running median algorithm. It then estimates the distribution of the normalized genome coverage with a Gaussian mixture model. Finally, a z-score statistic is assigned to each base position and used to separate the central distribution from the ROIs (i.e., under- and overcovered regions). A double thresholds mechanism is used to cluster the genomic ROIs. HTML reports provide a summary with interactive visual representations of the genomic ROIs with standard plots and metrics. Genomic variations such as single-nucleotide variants or CNVs can be effectively identified at the same time.
Collapse
Affiliation(s)
| | - Christiane Bouchier
- Institut Pasteur - Pole Biomics - 25-28 Rue du Docteur Roux, 75015 Paris, France
| | - Sean Kennedy
- Institut Pasteur - Pole Biomics - 25-28 Rue du Docteur Roux, 75015 Paris, France
| | - Thomas Cokelaer
- Institut Pasteur - Pole Biomics - 25-28 Rue du Docteur Roux, 75015 Paris, France.,Institut Pasteur - Bioinformatics and Biostatistics Hub - C3BI, USR 3756 IP CNRS - Paris, France
| |
Collapse
|
23
|
Abstract
Whole-genome sequencing with short-read technologies is well suited for calling single nucleotide polymorphisms, but has major problems with the detection of structural variants larger than the read length. One such type of variation is copy number variation (CNV), which entails deletion or duplication of genomic regions, and the expansion or contraction of repeated elements. Duplicated and deleted regions will typically be collapsed during de novo assembly of sequence data, or ignored when mapping reads toward a reference. However, signatures of the copy number variation can be detected in the resultant read depth at each position in the genome. We here provide instructions on how to analyze this read depth signal with the R package CNOGpro, allowing for estimation of copy numbers with uncertainty for each feature in a genome.
Collapse
|
24
|
Parallel Evolution of Group B Streptococcus Hypervirulent Clonal Complex 17 Unveils New Pathoadaptive Mutations. mSystems 2017; 2:mSystems00074-17. [PMID: 28904998 PMCID: PMC5585690 DOI: 10.1128/msystems.00074-17] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2017] [Accepted: 08/10/2017] [Indexed: 12/21/2022] Open
Abstract
The incidence of group B Streptococcus (GBS) neonatal disease continues to be a significant cause of concern worldwide. Strains belonging to clonal complex 17 (CC17) are the most frequently responsible for GBS infections in neonates, especially among late-onset disease cases. Therefore, we undertook the largest genomic study of GBS CC17 strains to date to decipher the genetic bases of their remarkable colonization and infection ability. We show that crucial functions involved in different steps of the colonization or infection process of GBS are distinctly mutated during the adaptation of CC17 to the human host. In particular, our results implicate the CovRS two-component regulator of virulence in the differentiation between carriage- and disease-associated isolates. Not only does this work raise important implications for the ongoing development of a vaccine against GBS but might also drive the discovery of key functions for GBS adaptation and pathogenesis that have been overlooked until now. Group B Streptococcus (GBS) is a commensal of the gastrointestinal and genitourinary tracts, while a prevailing cause of neonatal disease worldwide. Of the various clonal complexes (CCs), CC17 is overrepresented in GBS-infected newborns for reasons that are still largely unknown. Here, we report a comprehensive genomic analysis of 626 CC17 isolates collected worldwide, identifying the genetic traits behind their successful adaptation to humans and the underlying differences between carriage and clinical strains. Comparative analysis with 923 GBS genomes belonging to CC1, CC19, and CC23 revealed that the evolution of CC17 is distinct from that of other human-adapted lineages and recurrently targets functions related to nucleotide and amino acid metabolism, cell adhesion, regulation, and immune evasion. We show that the most distinctive features of disease-specific CC17 isolates were frequent mutations in the virulence-associated CovS and Stk1 kinases, underscoring the crucial role of the entire CovRS regulatory pathway in modulating the pathogenicity of GBS. Importantly, parallel and convergent evolution of major components of the bacterial cell envelope, such as the capsule biosynthesis operon, the pilus, and Rib, reflects adaptation to host immune pressures and should be taken into account in the ongoing development of a GBS vaccine. The presence of recurrent targets of evolution not previously implicated in virulence also opens the way for uncovering new functions involved in host colonization and GBS pathogenesis. IMPORTANCE The incidence of group B Streptococcus (GBS) neonatal disease continues to be a significant cause of concern worldwide. Strains belonging to clonal complex 17 (CC17) are the most frequently responsible for GBS infections in neonates, especially among late-onset disease cases. Therefore, we undertook the largest genomic study of GBS CC17 strains to date to decipher the genetic bases of their remarkable colonization and infection ability. We show that crucial functions involved in different steps of the colonization or infection process of GBS are distinctly mutated during the adaptation of CC17 to the human host. In particular, our results implicate the CovRS two-component regulator of virulence in the differentiation between carriage- and disease-associated isolates. Not only does this work raise important implications for the ongoing development of a vaccine against GBS but might also drive the discovery of key functions for GBS adaptation and pathogenesis that have been overlooked until now. Author Video: An author video summary of this article is available.
Collapse
|
25
|
Brynildsrud O, Gulla S, Feil EJ, Nørstebø SF, Rhodes LD. Identifying copy number variation of the dominant virulence factors msa and p22 within genomes of the fish pathogen Renibacterium salmoninarum. Microb Genom 2016; 2:e000055. [PMID: 28348850 PMCID: PMC5320689 DOI: 10.1099/mgen.0.000055] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2016] [Accepted: 03/17/2016] [Indexed: 11/24/2022] Open
Abstract
Renibacterium salmoninarum is the causative agent of bacterial kidney disease, an important disease of farmed and wild salmonid fish worldwide. Despite the wide spatiotemporal distribution of this disease and habitat pressures ranging from the natural environment to aquaculture and rivers to marine environments, little variation has been observed in the R. salmoninarum genome. Here we use the coverage depth from genomic sequencing corroborated by real-time quantitative PCR to detect copy number variation (CNV) among the genes of R. salmoninarum. CNV was primarily limited to the known dominant virulence factors msa and p22. Among 68 isolates representing the UK, Norway and North America, the msa gene ranged from two to five identical copies and the p22 gene ranged from one to five copies. CNV for these two genes co-occurred, suggesting they may be functionally linked. Isolates carrying CNV were phylogenetically restricted and originated predominantly from sites in North America, rather than the UK or Norway. Although both phylogenetic relationship and geographical origin were found to correlate with CNV status, geographical origin was a much stronger predictor than phylogeny, suggesting a role for local selection pressures in the repeated emergence and maintenance of this trait.
Collapse
Affiliation(s)
- Ola Brynildsrud
- 1Department of Bacteriology and Immunology Lovisenberggata 8, Norwegian Institute of Public Health/Department of Food Safety and Infection Biology, Norwegian University of Life Sciences (NMBU),Oslo,Norway
| | - Snorre Gulla
- 2Department of Bacteriology - Aquatic and Terrestrial Animals, Norwegian Veterinary Institute (NVI), Oslo, Norway
| | - Edward J Feil
- 3Department of Biology and Biochemistry, University of Bath,Claverton Down, Bath,United Kingdom
| | - Simen Foyn Nørstebø
- 4Department of Food Safety and Infection Biology, Norwegian University of Life Sciences (NMBU),Oslo,Norway
| | - Linda D Rhodes
- 5Northwest Fisheries Science Center, National Marine Fisheries Service, NOAA,Seattle, WA,United States
| |
Collapse
|
26
|
Perisin M, Vetter M, Gilbert JA, Bergelson J. 16Stimator: statistical estimation of ribosomal gene copy numbers from draft genome assemblies. ISME JOURNAL 2015; 10:1020-4. [PMID: 26359911 DOI: 10.1038/ismej.2015.161] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/07/2015] [Revised: 06/29/2015] [Accepted: 07/08/2015] [Indexed: 11/09/2022]
Abstract
The 16S rRNA gene (16S) is an accepted marker of bacterial taxonomic diversity, even though differences in copy number obscure the relationship between amplicon and organismal abundances. Ancestral state reconstruction methods can predict 16S copy numbers through comparisons with closely related reference genomes; however, the database of closed genomes is limited. Here, we extend the reference database of 16S copy numbers to de novo assembled draft genomes by developing 16Stimator, a method to estimate 16S copy numbers when these repetitive regions collapse during assembly. Using a read depth approach, we estimate 16S copy numbers for 12 endophytic isolates from Arabidopsis thaliana and confirm estimates by qPCR. We further apply this approach to draft genomes deposited in NCBI and demonstrate accurate copy number estimation regardless of sequencing platform, with an overall median deviation of 14%. The expanded database of isolates with 16S copy number estimates increases the power of phylogenetic correction methods for determining organismal abundances from 16S amplicon surveys.
Collapse
Affiliation(s)
- Matthew Perisin
- Department of Ecology and Evolution, University of Chicago, Chicago, IL, USA.,Committee on Microbiology, University of Chicago, Chicago, IL, USA
| | - Madlen Vetter
- Department of Ecology and Evolution, University of Chicago, Chicago, IL, USA
| | - Jack A Gilbert
- Department of Ecology and Evolution, University of Chicago, Chicago, IL, USA.,Biosciences Department, Institute for Genomic and Systems Biology, Argonne National Laboratory, Argonne, IL, USA.,Marine Biological Laboratory, Woods Hole, MA, USA.,College of Environmental and Resource Sciences, Zhejiang University, Hangzhou, China
| | - Joy Bergelson
- Department of Ecology and Evolution, University of Chicago, Chicago, IL, USA.,Committee on Microbiology, University of Chicago, Chicago, IL, USA
| |
Collapse
|