- Volume 6, Issue 5, 2020
Volume 6, Issue 5, 2020
- Review
-
- Microbial Communities
- Environmental
-
-
Microbial genomics amidst the Arctic crisis
The Arctic is warming – fast. Microbes in the Arctic play pivotal roles in feedbacks that magnify the impacts of Arctic change. Understanding the genome evolution, diversity and dynamics of Arctic microbes can provide insights relevant for both fundamental microbiology and interdisciplinary Arctic science. Within this synthesis, we highlight four key areas where genomic insights to the microbial dimensions of Arctic change are urgently required: the changing Arctic Ocean, greenhouse gas release from the thawing permafrost, 'biological darkening' of glacial surfaces, and human activities within the Arctic. Furthermore, we identify four principal challenges that provide opportunities for timely innovation in Arctic microbial genomics. These range from insufficient genomic data to develop unifying concepts or model organisms for Arctic microbiology to challenges in gaining authentic insights to the structure and function of low-biomass microbiota and integration of data on the causes and consequences of microbial feedbacks across scales. We contend that our insights to date on the genomics of Arctic microbes are limited in these key areas, and we identify priorities and new ways of working to help ensure microbial genomics is in the vanguard of the scientific response to the Arctic crisis.
-
- Mini Review
-
- Responses to Human Interventions
- Antibiotics
-
-
Antibiotic resistance in Pseudomonas aeruginosa and adaptation to complex dynamic environments
More LessAntibiotic resistance has become a serious threat to human health (WHO Antibacterial Agents in Clinical Development: an Analysis of the Antibacterial Clinical Development Pipeline, Including Tuberculosis. Geneva: World Health Organization; 2017), and the ability to predict antibiotic resistance from genome sequencing has become a focal point for the medical community. With this genocentric prediction in mind, we were intrigued about two particular findings for a collection of clinical Pseudomonas aeruginosa isolates (Marvig et al. Nature Genetics 2015;47:57–64; Frimodt-Møller et al. Scientific Reports 2018;8:12512; Bartell et al. Nature Communications 2019;10:629): (i) 15 out of 52 genes found to be frequently targeted by adaptive mutations during the initial infection stage of cystic fibrosis airways (‘candidate pathoadaptive genes’) (Marvig et al. Nature Genetics 2015;47:57–64) were associated with antibiotic resistance (López-Causapé et al. Fronters in Microbiology 2018;9:685; López-Causapé et al. Antimicrobal Agents and Chemotherapy 2018;62:e02583-17); (ii) there was a parallel lack of resistance development and linkage to the genetic changes in these antibiotic-resistance-associated genes (Frimodt-Møller et al. Scientific Reports 2018;8:12512; Bartell et al. Nature Communications 2019;10:629). In this review, we highlight alternative selective forces that potentially enhance the infection success of P. aeruginosa and focus on the linkage to the 15 pathoadaptive antibiotic-resistance-associated genes, thereby showing the problems we may face when using only genomic information to predict and inform about relevant antibiotic treatment.
-
- BioResource
-
- Microbial Evolution and Epidemiology
- Population Genomics
-
-
Visualizing variation within Global Pneumococcal Sequence Clusters (GPSCs) and country population snapshots to contextualize pneumococcal isolates
Rebecca A. Gladstone, Stephanie W. Lo, Richard Goater, Corin Yeats, Ben Taylor, James Hadfield, John A. Lees, Nicholas J. Croucher, Andries J. van Tonder, Leon J. Bentley, Fu Xiang Quah, Anne J. Blaschke, Nicole L. Pershing, Carrie L. Byington, Veeraraghavan Balaji, Waleria Hryniewicz, Betuel Sigauque, K.L. Ravikumar, Samanta Cristine Grassi Almeida, Theresa J. Ochoa, Pak Leung Ho, Mignon du Plessis, Kedibone M. Ndlangisa, Jennifer E. Cornick, Brenda Kwambana-Adams, Rachel Benisty, Susan A. Nzenze, Shabir A. Madhi, Paulina A. Hawkins, Andrew J. Pollard, Dean B. Everett, Martin Antonio, Ron Dagan, Keith P. Klugman, Anne von Gottberg, Benjamin J. Metcalf, Yuan Li, Bernard W. Beall, Lesley McGee, Robert F. Breiman, David M. Aanensen, Stephen D. Bentley and The Global Pneumococcal Sequencing ConsortiumKnowledge 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.
-
- Research Article
-
- Microbial Evolution and Epidemiology
- Population Genomics
-
-
Genomic analysis of 40 prophages located in the genomes of 16 carbapenemase-producing clinical strains of Klebsiella pneumoniae
Klebsiella pneumoniae is the clinically most important species within the genus Klebsiella and, as a result of the continuous emergence of multi-drug resistant (MDR) strains, the cause of severe nosocomial infections. The decline in the effectiveness of antibiotic treatments for infections caused by MDR bacteria has generated particular interest in the study of bacteriophages. In this study, we characterized a total of 40 temperate bacteriophages (prophages) with a genome range of 11.454–84.199 kb, predicted from 16 carbapenemase-producing clinical strains of K. pneumoniae belonging to different sequence types, previously identified by multilocus sequence typing. These prophages were grouped into the three families in the order Caudovirales (27 prophages belonging to the family Myoviridae, 10 prophages belonging to the family Siphoviridae and 3 prophages belonging to the family Podoviridae). Genomic comparison of the 40 prophage genomes led to the identification of four prophages isolated from different strains and of genome sizes of around 33.3, 36.1, 39.6 and 42.6 kb. These prophages showed sequence similarities (query cover >90 %, identity >99.9 %) with international Microbe Versus Phage (MVP) (http://mvp.medgenius.info/home) clusters 4762, 4901, 3499 and 4280, respectively. Phylogenetic analysis revealed the evolutionary proximity among the members of the four groups of the most frequently identified prophages in the bacterial genomes studied (33.3, 36.1, 39.6 and 42.6 kb), with bootstrap values of 100 %. This allowed the prophages to be classified into three clusters: A, B and C. Interestingly, these temperate bacteriophages did not infect the highest number of strains as indicated by a host-range assay, these results could be explained by the development of superinfection exclusion mechanisms. In addition, bioinformatic analysis of the 40 identified prophages revealed the presence of 2363 proteins. In total, 59.7 % of the proteins identified had a predicted function, mainly involving viral structure, transcription, replication and regulation (lysogenic/lysis). Interestingly, some proteins had putative functions associated with bacterial virulence (toxin expression and efflux pump regulators), phage defence profiles such as toxin–antitoxin modules, an anti-CRISPR/Cas9 protein, TerB protein (from terZABCDE operon) and methyltransferase proteins.
-
-
-
Genomic profiling of Escherichia coli isolates from bacteraemia patients: a 3-year cohort study of isolates collected at a Sydney teaching hospital
This study sought to assess the genetic variability of Escherichia coli isolated from bloodstream infections (BSIs) presenting at Concord Hospital, Sydney during 2013–2016. Whole-genome sequencing was used to characterize 81 E. coli isolates sourced from community-onset (CO) and hospital-onset (HO) BSIs. The cohort comprised 64 CO and 17 HO isolates, including 35 multidrug-resistant (MDR) isolates exhibiting phenotypic resistance to three or more antibiotic classes. Phylogenetic analysis identified two major ancestral clades. One was genetically diverse with 25 isolates distributed in 16 different sequence types (STs) representing phylogroups A, B1, B2, C and F, while the other comprised phylogroup B2 isolates in subclades representing the ST131, ST73 and ST95 lineages. Forty-seven isolates contained a class 1 integron, of which 14 carried bla CTX -M-gene. Isolates with a class 1 integron carried more antibiotic resistance genes than isolates without an integron and, in most instances, resistance genes were localized within complex resistance loci (CRL). Resistance to fluoroquinolones could be attributed to point mutations in chromosomal parC and gyrB genes and, in addition, two isolates carried a plasmid-associated qnrB4 gene. Co-resistance to fluoroquinolone and broad-spectrum beta-lactam antibiotics was associated with ST131 (HO and CO), ST38 (HO), ST393 (CO), ST2003 (CO) and ST8196 (CO and HO), a novel ST identified in this study. Notably, 10/81 (12.3 %) isolates with ST95 (5 isolates), ST131 (2 isolates), ST88 (2 isolates) and a ST540 likely carry IncFII–IncFIB plasmid replicons with a full spectrum of virulence genes consistent with the carriage of ColV-like plasmids. Our data indicate that IncF plasmids play an important role in shaping virulence and resistance gene carriage in BSI E. coli in Australia.
-
- Phylogeography
-
-
Origin, genomic diversity and microevolution of the Clostridium difficile B1/NAP1/RT027/ST01 strain in Costa Rica, Chile, Honduras and Mexico
Clostridium difficile B1/NAP1/RT027/ST01 has been responsible for outbreaks of antibiotic-associated diarrhoea in clinical settings worldwide and is associated with severe disease presentations and increased mortality rates. Two fluoroquinolone-resistant (FQR) lineages of the epidemic B1/NAP1/RT027/ST01 strain emerged in the USA in the early 1990s and disseminated trans continentally (FQR1 and FQR2). However, it is unclear when and from where they entered Latin America (LA) and whether isolates from LA exhibit unique genomic features when compared to B1/NAP1/RT027/ST01 isolates from other regions of the world. To answer the first issue we compared whole-genome sequences (WGS) of 25 clinical isolates typed as NAP1, RT027 or ST01 in Costa Rica (n=16), Chile (n=5), Honduras (n=3) and Mexico (n=1) to WGS of 129 global isolates from the same genotype using Bayesian phylogenomics. The second question was addressed through a detailed analysis of the number and type of mutations of the LA isolates and their mobile resistome. All but two B1/NAP1/RT027/ST01 isolates from LA belong to the FQR2 lineage (n=23, 92 %), confirming its widespread distribution. As indicated by analysis of a dataset composed of 154 WGS, the B1/NAP1/RT027/ST01 strain was introduced into the four LA countries analysed between 1998 and 2005 from North America (twice) and Europe (at least four times). These events occurred soon after the emergence of the FQR lineages and more than one decade before the first report of the detection of the B1/NAP1/RT027/ST01 in LA. A total of 552 SNPs were identified across all genomes examined (3.8–4.3 Mb) in pairwise comparisons to the R20291 reference genome. Moreover, pairwise SNP distances were among the smallest distances determined in this species so far (0 to 55). Despite this high level of genomic conservation, 39 unique SNPs (7 %) in genes that play roles in the infection process (i.e. slpA) or antibiotic resistance (i.e. rpoB, fusA) distinguished the LA isolates. In addition, isolates from Chile, Honduras and Mexico had twice as many antibiotic resistance genes (ARGs, n=4) than related isolates from other regions. Their unique set of ARGs includes a cfr-like gene and tetM, which were found as part of putative mobile genetic elements whose sequences resemble undescribed integrative and conjugative elements. These results show multiple, independent introductions of B1/NAP1/RT027/ST01 isolates from the FQR1 and FQR2 lineages from different geographical sources into LA and a rather rapid accumulation of distinct mutations and acquired ARG by the LA isolates.
-
-
-
Molecular epidemiology of Salmonella Infantis in Europe: insights into the success of the bacterial host and its parasitic pESI-like megaplasmid
Salmonella Infantis is one of the five serovars most frequently causing human salmonellosis in Europe, mainly associated with poultry. A clone harbouring a conjugative plasmid of emerging S. Infantis (pESI)-like megaplasmid, carrying multidrug resistant (MDR) and extended-spectrum beta-lactamases (ESBL) genes, has spread in the Italian broiler chicken industry also causing human illness. This work is aimed at elucidating the molecular epidemiology of S. Infantis and pESI-like in Europe using whole-genome sequencing and bioinformatics analysis, and to investigate the genetic relatedness of S. Infantis clones and pESI-like from animals, meat, feed and humans provided by institutions of nine European countries. Two genotyping approaches were used: chromosome or plasmid SNP-based analysis and the minimum spanning tree (MST) algorithm based on core-genome multilocus sequence typing (cgMLST). The European S. Infantis population appeared heterogeneous, with different genetic clusters defined at core-genome level. However, pESI-like variants present in 64.1 % of the isolates were more genetically homogeneous and capable of infecting different clonal lineages in most of the countries. Two different pESI-like with ESBL genes (n=82) were observed: bla CTX-M-1-positive in European isolates and bla CTX-M-65-positive in American isolates (study outgroup). Both variants had toxin-antitoxin systems, resistance genes towards tetracyclines, trimethoprim, sulphonamides and aminoglycosides, heavy metals (merA) and disinfectants (qacEΔ). Worryingly, 66 % of the total isolates studied presented different gyrA chromosomal point mutations associated with (fluoro)quinolone resistance (MIC range 0.125–0.5 mg/L), while 18 % displayed transferable macrolide resistance mediated by mph, mef and erm(B) genes. Proper intervention strategies are needed to prevent further dissemination/transmission of MDR S. Infantis and pESI-like along the food chain in Europe.
-
- Mechanisms of Evolution
-
-
Transposable elements contribute to the genome plasticity of Ralstonia solanacearum species complex
The extensive genetic diversity of Ralstonia solanacearum , a serious soil-borne phytopathogen, has led to the concept that R. solanacearum encompasses a species complex [ R. solanacearum species complex (RSSC)]. Insertion sequences (ISs) are suggested to play an important role in the genome evolution of this pathogen. Here, we identified and analysed transposable elements (TEs), ISs and transposons, in 106 RSSC genomes and 15 Ralstonia spp. We mapped 10 259 IS elements in the complete genome of 62 representative RSSC strains and closely related Ralstonia spp. A unique set of 20 IS families was widespread across the strains, IS5 and IS3 being the most abundant. Our results showed six novel transposon sequences belonging to the Tn3 family carrying passenger genes encoding antibiotic resistance and avirulence proteins. In addition, internal rearrangement events associated with ISs were demonstrated in Ralstonia pseudosolanacearum strains. We also mapped IS elements interrupting avirulence genes, which provided evidence that ISs plays an important role in virulence evolution of RSSC. Additionally, the activity of ISs was demonstrated by transcriptome analysis and DNA hybridization in R. solanacearum isolates. Altogether, we have provided collective data of TEs in RSSC genomes, opening a new path for understanding their evolutionary impact on the genome evolution and diversity of this important plant pathogen.
-
- Microbe-Niche Interactions
- Host Adaptation
-
-
Meningococcal core and accessory phasomes vary by clonal complex
More LessNeisseria meningitidis is a Gram-negative human commensal pathogen, with extensive phenotypic plasticity afforded by phase-variable (PV) gene expression. Phase variation is a stochastic switch in gene expression from an ON to an OFF state, mediated by localized hypermutation of simple sequence repeats (SSRs). Circulating N. meningitidis clones vary in propensity to cause disease, with some clonal complexes (ccs) classified as hypervirulent and others as carriage-associated. We examined the PV gene repertoires, or phasome, of these lineages in order to determine whether phase variation contributes to disease propensity. We analysed 3328 genomes representative of nine circulating meningococcal ccs with PhasomeIt, a tool that identifies PV genes by the presence of SSRs and homologous gene clusters. The presence, absence and functions of all identified PV gene clusters were confirmed by annotation or blast searches within the Neisseria PubMLST database. While no significant differences were detected in the number of PV genes or the core, conserved phasome content between hypervirulent and carriage lineages, individual ccs exhibited major variations in PV gene numbers. Phylogenetic clusters produced by phasome or core genome analyses were similar, indicating co-evolution of PV genes with the core genome. While conservation of PV clusters is high, with 76 % present in all meningococcal isolates, maintenance of an SSR is variable, ranging from conserved in all isolates to present only in a single cc, indicating differing evolutionary trajectories for each lineage. Diverse functional groups of PV genes were present across the meningococcal lineages; however, the majority directly or indirectly influence bacterial surface antigens and could impact on future vaccine development. Finally, we observe that meningococci have open pan phasomes, indicating ongoing evolution of PV gene content and a significant potential for adaptive changes in this clinically relevant genus.
-
- Genomic Methodologies
- Genome Variation Detection
-
-
Identification of integrative and conjugative elements in pathogenic and commensal Neisseriaceae species via genomic distributions of DNA uptake sequence dialects
More LessMobile genetic elements (MGEs) are key factors responsible for dissemination of virulence determinants and antimicrobial-resistance genes amongst pathogenic bacteria. Conjugative MGEs are notable for their high gene loads donated per transfer event, broad host ranges and phylogenetic ubiquity amongst prokaryotes, with the subclass of chromosomally inserted integrative and conjugative elements (ICEs) being particularly abundant. The focus on a small number of model systems has biased the study of ICEs towards those conferring readily selectable phenotypes to host cells, whereas the identification and characterization of integrated cryptic elements remains challenging. Even though antimicrobial resistance and horizontally acquired virulence genes are major factors aggravating neisserial infection, conjugative MGEs of Neisseria gonorrhoeae and Neisseria meningitidis remain poorly characterized. Using a phenotype-independent approach based on atypical distributions of DNA uptake sequences (DUSs) in MGEs relative to the chromosomal background, we have identified two groups of chromosomally integrated conjugative elements in Neisseria : one found almost exclusively in pathogenic species possibly deriving from the genus Kingella , the other belonging to a group of Neisseria mucosa -like commensals. The former element appears to enable transfer of traditionally gonococcal-specific loci such as the virulence-associated toxin–antitoxin system fitAB to N. meningitidis chromosomes, whilst the circular form of the latter possesses a unique attachment site (attP) sequence seemingly adapted to exploit DUS motifs as chromosomal integration sites. In addition to validating the use of DUS distributions in Neisseriaceae MGE identification, the >170 identified ICE sequences provide a valuable resource for future studies of ICE evolution and host adaptation.
-
- Systems Microbiology
- Large-scale Comparative Genomics
-
-
Computational genomic discovery of diverse gene clusters harbouring Fe-S flavoenzymes in anaerobic gut microbiota
More LessThe gut contains an enormous diversity of simple as well as complex molecules from highly diverse food sources, together with host-secreted molecules. This presents a large metabolic opportunity for the gut microbiota, but little is known about how gut microbes are able to catabolize this large chemical diversity. Recently, Fe-S flavoenzymes were found to be key in the transformation of bile acids, catalysing the key step in the 7α-dehydroxylation pathway that allows gut bacteria to transform cholic acid into deoxycholic acid, an exclusively microbe-derived molecule with major implications for human health. While this enzyme family has also been implicated in a limited number of other catalytic transformations, little is known about the extent to which it is of more global importance in gut microbial metabolism. Here, we perform a large-scale computational genomic analysis to show that this enzyme superfamily has undergone a remarkable expansion in Clostridiales, and occurs throughout a diverse array of >1000 different families of putative metabolic gene clusters. Analysis of the enzyme content of these gene clusters suggests that they encode pathways with a wide range of predicted substrate classes, including saccharides, amino acids/peptides and lipids. Altogether, these results indicate a potentially important role of this protein superfamily in the human gut, and our dataset provides significant opportunities for the discovery of novel pathways that may have significant effects on human health.
-
- Short Communication
-
- Systems Microbiology
- Pangenome Analysis, Big-data Approaches
-
-
Accelerating surveillance and research of antimicrobial resistance – an online repository for sharing of antimicrobial susceptibility data associated with whole-genome sequences
Antimicrobial resistance (AMR) is an emerging threat to modern medicine. Improved diagnostics and surveillance of resistant bacteria require the development of next-generation analysis tools and collaboration between international partners. Here, we present the ‘AMR Data Hub’, an online infrastructure for storage and sharing of structured phenotypic AMR data linked to bacterial whole-genome sequences. Leveraging infrastructure built by the European COMPARE Consortium and structured around the European Nucleotide Archive (ENA), the AMR Data Hub already provides an extensive data collection of more than 2500 isolates with linked genome and AMR data. Representing these data in standardized formats, we provide tools for the validation and submission of new data and services supporting search, browse and retrieval. The current collection was created through a collaboration by several partners from the European COMPARE Consortium, demonstrating the capacities and utility of the AMR Data Hub and its associated tools. We anticipate growth of content and offer the hub as a basis for future research into methods to explore and predict AMR.
-
- Method
-
- Genomic Methodologies
- Genome-phenotype Association
-
-
prewas: data pre-processing for more informative bacterial GWAS
More LessWhile variant identification pipelines are becoming increasingly standardized, less attention has been paid to the pre-processing of variants prior to their use in bacterial genome-wide association studies (bGWAS). Three nuances of variant pre-processing that impact downstream identification of genetic associations include the separation of variants at multiallelic sites, separation of variants in overlapping genes, and referencing of variants relative to ancestral alleles. Here we demonstrate the importance of these variant pre-processing steps on diverse bacterial genomic datasets and present prewas, an R package, that standardizes the pre-processing of multiallelic sites, overlapping genes, and reference alleles before bGWAS. This package facilitates improved reproducibility and interpretability of bGWAS results. prewas enables users to extract maximal information from bGWAS by implementing multi-line representation for multiallelic sites and variants in overlapping genes. prewas outputs a binary SNP matrix that can be used for SNP-based bGWAS and will prevent the masking of minor alleles during bGWAS analysis. The optional binary gene matrix output can be used for gene-based bGWAS, which will enable users to maximize the power and evolutionary interpretability of their bGWAS studies. prewas is available for download from GitHub.
-