Search

Funder: Wellcome Trust (Wellcome); doi: https://doi.org/10.13039/100004440, Funder: Cystic Fibrosis Trust (CF); doi: https://doi.org/10.13039/501100000292, Funder: Fondation Botnar; doi: https://doi.org/10.13039/501100011318, Funder: European Molecular Biology Organization (EMBO); doi: https://doi.org/10.13039/100004410, Funder: Wellcome Trust (Wellcome); doi: https://doi.org/10.13039/100010269, The medical and scientific response to emerging and established pathogens is often severely hampered by ignorance of the genetic determinants of virulence, drug resistance and clinical outcomes that could be used to identify therapeutic drug targets and forecast patient trajectories. Taking the newly emergent multidrug-resistant bacteria Mycobacterium abscessus as an example, we show that combining high-dimensional phenotyping with whole-genome sequencing in a phenogenomic analysis can rapidly reveal actionable systems-level insights into bacterial pathobiology. Through phenotyping of 331 clinical isolates, we discovered three distinct clusters of isolates, each with different virulence traits and associated with a different clinical outcome. We combined genome-wide association studies with proteome-wide computational structural modelling to define likely causal variants, and employed direct coupling analysis to identify co-evolving, and therefore potentially epistatic, gene networks. We then used in vivo CRISPR-based silencing to validate our findings and discover clinically relevant M. abscessus virulence factors including a secretion system, thus illustrating how phenogenomics can reveal critical pathways within emerging pathogenic bacteria.

Mycobacterium abscessus, a multidrug-resistant nontuberculous mycobacterium, has emerged as a major pathogen affecting people with cystic fibrosis (CF). Although originally thought to be acquired independently from the environment, most individuals are infected with one of several dominant circulating clones (DCCs), indicating the presence of global transmission networks of M. abscessus. How and when these clones emerged and spread globally is unclear. Here, we use evolutionary analyses of isolates from individuals both with and without CF to reconstruct the population history, spatiotemporal spread and recent transmission networks of the DCCs. We demonstrate synchronous expansion of six unrelated DCCs in the 1960s, a period associated with major changes in CF care and survival. Each of these clones has spread globally as a result of rare intercontinental transmission events. We show that the DCCs, but not environmentally acquired isolates, exhibit a specific smoking-associated mutational signature and that current transmission networks include individuals both with and without CF. We therefore propose that the DCCs initially emerged in non-CF populations but were then amplified and spread through the CF community. While individuals with CF are probably the most permissive host, non-CF individuals continue to play a key role in transmission networks and may facilitate long-distance transmission., In this Article, the authors perform evolutionary analyses of M. abscessus clinical isolates and report the emergence of dominant circulating clones (DCCs) in non-cystic fibrosis (CF) individuals followed by amplification in the CF community.

Funder: Department of Health {UK), Genome-wide association studies (GWAS) are increasingly being applied to investigate the genetic basis of bacterial traits. However, approaches to perform power calculations for bacterial GWAS are limited. Here we implemented two alternative approaches to conduct power calculations using existing collections of bacterial genomes. First, a sub-sampling approach was undertaken to reduce the allele frequency and effect size of a known and detectable genotype-phenotype relationship by modifying phenotype labels. Second, a phenotype-simulation approach was conducted to simulate phenotypes from existing genetic variants. We implemented both approaches into a computational pipeline (PowerBacGWAS) that supports power calculations for burden testing, pan-genome and variant GWAS; and applied it to collections of Enterococcus faecium, Klebsiella pneumoniae and Mycobacterium tuberculosis. We used this pipeline to determine sample sizes required to detect causal variants of different minor allele frequencies (MAF), effect sizes and phenotype heritability, and studied the effect of homoplasy and population diversity on the power to detect causal variants. Our pipeline and user documentation are made available and can be applied to other bacterial populations. PowerBacGWAS can be used to determine sample sizes required to find statistically significant associations, or the associations detectable with a given sample size. We recommend to perform power calculations using existing genomes of the bacterial species and population of study.

Escherichia coli is the leading cause of urinary tract infection, one of the most common bacterial infections in humans. Despite this, a genomic perspective is lacking regarding the phylogenetic distribution of isolates associated with different clinical syndromes. Here, we present a large-scale phylogenomic analysis of a spatiotemporally and clinically diverse set of 907 E. coli isolates, including 722 uropathogenic E. coli (UPEC) isolates. A genome-wide association approach identifies the (P-fimbriae-encoding) papGII locus as the key feature distinguishing invasive UPEC, defined as isolates associated with severe UTI, i.e., kidney infection (pyelonephritis) or urinary-source bacteremia, from non-invasive UPEC, defined as isolates associated with asymptomatic bacteriuria or bladder infection (cystitis). Within the E. coli population, distinct invasive UPEC lineages emerged through repeated horizontal acquisition of diverse papGII-containing pathogenicity islands. Our findings elucidate the molecular determinants of severe UTI and have implications for the early detection of this pathogen.

Funder: KEH is supported by a Senior Medical Research Fellowship from the Viertel Foundation of Australia RJB through a UKRI MRC New Investigator Research Grant (held by KSB, MR/R020787/1), Shigella sonnei is the most common agent of shigellosis in high-income countries, and causes a significant disease burden in low- and middle-income countries. Antimicrobial resistance is increasingly common in all settings. Whole genome sequencing (WGS) is increasingly utilised for S. sonnei outbreak investigation and surveillance, but comparison of data between studies and labs is challenging. Here, we present a genomic framework and genotyping scheme for S. sonnei to efficiently identify genotype and resistance determinants from WGS data. The scheme is implemented in the software package Mykrobe and tested on thousands of genomes. Applying this approach to analyse >4,000 S. sonnei isolates sequenced in public health labs in three countries identified several common genotypes associated with increased rates of ciprofloxacin resistance and azithromycin resistance, confirming intercontinental spread of highly-resistant S. sonnei clones and demonstrating the genomic framework can facilitate monitoring the spread of resistant clones, including those that have recently emerged, at local and global scales.

In order to control and eradicate epidemic cholera, we need to understand how epidemics begin, how they spread, and how they decline and eventually end. This requires extensive sampling of epidemic disease over time, alongside the background of endemic disease that may exist concurrently with the epidemic. The unique circumstances surrounding the Argentinian cholera epidemic of 1992–1998 presented an opportunity to do this. Here, we use 490 Argentinian V. cholerae genome sequences to characterise the variation within, and between, epidemic and endemic V. cholerae. We show that, during the 1992–1998 cholera epidemic, the invariant epidemic clone co-existed alongside highly diverse members of the Vibrio cholerae species in Argentina, and we contrast the clonality of epidemic V. cholerae with the background diversity of local endemic bacteria. Our findings refine and add nuance to our genomic definitions of epidemic and endemic cholera, and are of direct relevance to controlling current and future cholera epidemics., Pandemic cholera was reintroduced to Argentina in 1992, leading to epidemic spread. Here, the authors use whole genome sequencing to show how, over 6 years, epidemic cholera was caused by invariant 7PET lineage Vibrio cholerae, against a background of sporadic disease caused by diverse local strains.

Restriction-Modification systems (RMS) are one of the main mechanisms of defence against foreign DNA invasion and can have an important role in the regulation of gene expression. The obligate human pathogen Neisseria gonorrhoeae carries one of the highest loads of RMS in its genome; between 13 to 15 of the three main types. Previous work has described their organization in the reference genome FA1090 and has experimentally inferred the associated methylated motifs. Here, we studied the structure of RMS and target methylated motifs in 25 gonococcal strains sequenced with Single Molecule Real-Time (SMRT) technology, which provides data on DNA modification. The results showed a variable picture of active RMS in different strains, with phase variation switching the activity of Type III RMS, and both the activity and specificity of a Type I RMS. Interestingly, the Dam methylase was found in place of the NgoAXI endonuclease in two of the strains, despite being previously thought to be absent in the gonococcus. We also identified the real methylation target of NgoAX as 5’-GCAGA-3’, different from that previously described. Results from this work give further insights into the diversity and dynamics of RMS and methylation patterns in N. gonorrhoeae.

Of the ten human-restricted Neisseria species two, Neisseria meningitidis, and Neisseria gonorrhoeae, cause invasive disease: the other eight are carried asymptomatically in the pharynx, possibly modulating meningococcal and gonococcal infections. Consequently, characterizing their diversity is important for understanding the microbiome in health and disease. Whole genome sequences from 181 Neisseria isolates were examined, including those of three well-defined species (N. meningitidis; N. gonorrhoeae; and Neisseria polysaccharea) and genomes of isolates unassigned to any species (Nspp). Sequence analysis of ribosomal genes, and a set of core (cgMLST) genes were used to infer phylogenetic relationships. Average Nucleotide Identity (ANI) and phenotypic data were used to define species clusters, and morphological and metabolic differences among them. Phylogenetic analyses identified two polyphyletic clusters (N. polysaccharea and Nspp.), while, cgMLST data grouped Nspp isolates into nine clusters and identified at least three N. polysaccharea clusters. ANI results classified Nspp into seven putative species, and also indicated at least three putative N. polysaccharea species. Electron microscopy identified morphological differences among these species. This genomic approach provided a consistent methodology for species characterization using distinct phylogenetic clusters. Seven putative novel Neisseria species were identified, confirming the importance of genomic studies in the characterization of the genus Neisseria.

Bloodstream infections by Salmonella enterica serovar Typhimurium constitute a major health burden in sub-Saharan Africa (SSA). These invasive non-typhoidal (iNTS) infections are dominated by isolates of the antibiotic resistance-associated sequence type (ST) 313. Here, we report emergence of ST313 sublineage II.1 in the Democratic Republic of the Congo. Sublineage II.1 exhibits extensive drug resistance, involving a combination of multidrug resistance, extended spectrum β-lactamase production and azithromycin resistance. ST313 lineage II.1 isolates harbour an IncHI2 plasmid we name pSTm-ST313-II.1, with one isolate also exhibiting decreased ciprofloxacin susceptibility. Whole genome sequencing reveals that ST313 II.1 isolates have accumulated genetic signatures potentially associated with altered pathogenicity and host adaptation, related to changes observed in biofilm formation and metabolic capacity. Sublineage II.1 emerged at the beginning of the 21st century and is involved in on-going outbreaks. Our data provide evidence of further evolution within the ST313 clade associated with iNTS in SSA., Invasive non-typhoidal Salmonella (iNTS) infections are dominated by antibiotic resistant isolates of the sequence type (ST) 313. Here, the authors identify the ST313 sublineage II.1 in the Democratic Republic of the Congo exhibiting extensive drug resistance and genetic signatures potentially associated with host adaptation.

The environmental bacterium Burkholderia pseudomallei causes melioidosis, an important endemic human disease in tropical and sub-tropical countries. This bacterium occupies broad ecological niches including soil, contaminated water, single-cell microbes, plants and infection in a range of animal species. Here, we performed genome-wide association studies for genetic determinants of environmental and human adaptation using a combined dataset of 1,010 whole genome sequences of B. pseudomallei from Northeast Thailand and Australia, representing two major disease hotspots. With these data, we identified 47 genes from 26 distinct loci associated with clinical or environmental isolates from Thailand and replicated 12 genes in an independent Australian cohort. We next outlined the selective pressures on the genetic loci (dN/dS) and the frequency at which they had been gained or lost throughout their evolutionary history, reflecting the bacterial adaptability to a wide range of ecological niches. Finally, we highlighted loci likely implicated in human disease., Claire Chewapreecha et al. combine 753 newly sequenced Thai Burkholderia pseudomallei isolates with 258 Australian isolates to identify genes associated with either clinical or environmental strains. They find 47 genes, 12 of which replicate in both locations, that may provide clues to the strategy used by this microbe to adapt to survive in wide range of ecological niches, including human hosts.

Bloodstream infections by Salmonella enterica serovar Typhimurium constitute a major health burden in sub-Saharan Africa (SSA). These invasive non-typhoidal (iNTS) infections are dominated by isolates of the antibiotic resistance-associated sequence type (ST) 313. Here, we report emergence of ST313 sublineage II.1 in the Democratic Republic of the Congo. Sublineage II.1 exhibits extensive drug resistance, involving a combination of multidrug resistance, extended spectrum β-lactamase production and azithromycin resistance. ST313 lineage II.1 isolates harbour an IncHI2 plasmid we name pSTm-ST313-II.1, with one isolate also exhibiting decreased ciprofloxacin susceptibility. Whole genome sequencing reveals that ST313 II.1 isolates have accumulated genetic signatures potentially associated with altered pathogenicity and host adaptation, related to changes observed in biofilm formation and metabolic capacity. Sublineage II.1 emerged at the beginning of the 21st century and is involved in on-going outbreaks. Our data provide evidence of further evolution within the ST313 clade associated with iNTS in SSA.

Of the ten human-restricted Neisseria species two, Neisseria meningitidis, and Neisseria gonorrhoeae, cause invasive disease: the other eight are carried asymptomatically in the pharynx, possibly modulating meningococcal and gonococcal infections. Consequently, characterizing their diversity is important for understanding the microbiome in health and disease. Whole genome sequences from 181 Neisseria isolates were examined, including those of three well-defined species (N. meningitidis; N. gonorrhoeae; and Neisseria polysaccharea) and genomes of isolates unassigned to any species (Nspp). Sequence analysis of ribosomal genes, and a set of core (cgMLST) genes were used to infer phylogenetic relationships. Average Nucleotide Identity (ANI) and phenotypic data were used to define species clusters, and morphological and metabolic differences among them. Phylogenetic analyses identified two polyphyletic clusters (N. polysaccharea and Nspp.), while, cgMLST data grouped Nspp isolates into nine clusters and identified at least three N. polysaccharea clusters. ANI results classified Nspp into seven putative species, and also indicated at least three putative N. polysaccharea species. Electron microscopy identified morphological differences among these species. This genomic approach provided a consistent methodology for species characterization using distinct phylogenetic clusters. Seven putative novel Neisseria species were identified, confirming the importance of genomic studies in the characterization of the genus Neisseria.

Funder: U.S. Department of Health & Human Services | National Institutes of Health (NIH), In order to control and eradicate epidemic cholera, we need to understand how epidemics begin, how they spread, and how they decline and eventually end. This requires extensive sampling of epidemic disease over time, alongside the background of endemic disease that may exist concurrently with the epidemic. The unique circumstances surrounding the Argentinian cholera epidemic of 1992–1998 presented an opportunity to do this. Here, we use 490 Argentinian V. cholerae genome sequences to characterise the variation within, and between, epidemic and endemic V. cholerae. We show that, during the 1992–1998 cholera epidemic, the invariant epidemic clone co-existed alongside highly diverse members of the Vibrio cholerae species in Argentina, and we contrast the clonality of epidemic V. cholerae with the background diversity of local endemic bacteria. Our findings refine and add nuance to our genomic definitions of epidemic and endemic cholera, and are of direct relevance to controlling current and future cholera epidemics.

Funder: U.S. Department of Health & Human Services | National Institutes of Health (NIH), In order to control and eradicate epidemic cholera, we need to understand how epidemics begin, how they spread, and how they decline and eventually end. This requires extensive sampling of epidemic disease over time, alongside the background of endemic disease that may exist concurrently with the epidemic. The unique circumstances surrounding the Argentinian cholera epidemic of 1992–1998 presented an opportunity to do this. Here, we use 490 Argentinian V. cholerae genome sequences to characterise the variation within, and between, epidemic and endemic V. cholerae. We show that, during the 1992–1998 cholera epidemic, the invariant epidemic clone co-existed alongside highly diverse members of the Vibrio cholerae species in Argentina, and we contrast the clonality of epidemic V. cholerae with the background diversity of local endemic bacteria. Our findings refine and add nuance to our genomic definitions of epidemic and endemic cholera, and are of direct relevance to controlling current and future cholera epidemics.

Funder: Thailand National Science and Technology Development Agency (FDA-CO-2562-8764-TH) and Thailand Science Research and Innovation fund (MRG6280226), Funder: Food Standards Agency Fellow and is supported by the BBSRC Institute Strategic Programme Microbes in the Food Chain BB/R012504/1 and its constituent projects BBS/E/F/000PR10348 (Theme 1, Epidemiology and Evolution of Pathogens in the Food Chain) and BBS/E/F/000PR10351 (Theme 3, Microbial Communities in the Food Chain), The environmental bacterium Burkholderia pseudomallei causes melioidosis, an important endemic human disease in tropical and sub-tropical countries. This bacterium occupies broad ecological niches including soil, contaminated water, single-cell microbes, plants and infection in a range of animal species. Here, we performed genome-wide association studies for genetic determinants of environmental and human adaptation using a combined dataset of 1,010 whole genome sequences of B. pseudomallei from Northeast Thailand and Australia, representing two major disease hotspots. With these data, we identified 47 genes from 26 distinct loci associated with clinical or environmental isolates from Thailand and replicated 12 genes in an independent Australian cohort. We next outlined the selective pressures on the genetic loci (dN/dS) and the frequency at which they had been gained or lost throughout their evolutionary history, reflecting the bacterial adaptability to a wide range of ecological niches. Finally, we highlighted loci likely implicated in human disease.