Detecting co-selection through excess linkage disequilibrium in bacterial genomes.

Population genomics has revolutionized our ability to study bacterial evolution by enabling data-driven discovery of the genetic architecture of trait variation. Genome-wide association studies (GWAS) have more recently become accompanied by genome-wide epistasis and co-selection (GWES) analysis, which offers a phenotype-free approach to generating hypotheses about selective processes that simultaneously impact multiple loci across the genome. However, existing GWES methods only consider associations between distant pairs of loci within the genome due to the strong impact of linkage-disequilibrium (LD) over short distances. Based on the general functional organisation of genomes it is nevertheless expected that majority of co-selection and epistasis will act within relatively short genomic proximity, on co-variation occurring within genes and their promoter regions, and within operons. Here, we introduce LDWeaver, which enables an exhaustive GWES across both short- and long-range LD, to disentangle likely neutral co-variation from selection. We demonstrate the ability of LDWeaver to efficiently generate hypotheses about co-selection using large genomic surveys of multiple major human bacterial pathogen species and validate several findings using functional annotation and phenotypic measurements. Our approach will facilitate the study of bacterial evolution in the light of rapidly expanding population genomic data.

Population structure and antimicrobial resistance among Klebsiella isolates sampled from human, animal, and environmental sources in Ghana: a cross-sectional genomic One Health study.

BACKGROUND: One Health approaches to address the increasing threat of antimicrobial resistance (AMR) are gaining attention. However, data on the distribution and movement of bacteria and their AMR-associated genes between clinical and non-clinical sources are scarce, especially from low-income and middle-income countries. We aimed to analyse Klebsiella isolates from various sources in Ghana and compare the prevalence of AMR with datasets from two other countries. METHODS: We conducted a cross-sectional genomic One Health study. Multiple clinical, environmental, and animal sources were sampled from 78 locations (eg, hospitals, residential areas, and farms) in and around Tamale, Ghana. Clinical samples were collected through routine screening and in cases of suspected infection between March 15 and Sept 15, 2019, and samples from the wider environment were collected during a dedicated sampling effort between the dates of Aug 19, 2018, and Sept 26, 2019. Sampling locations were approximately evenly distributed from the centre of the city and steadily outwards to capture both rural and urban locations. Samples with positive growth for Klebsiella were included. Isolates of Klebsiella were obtained from the samples using Simmons citrate agar medium and characterised by antimicrobial susceptibility testing and whole-genome sequencing. A comparative analysis with Klebsiella population surveys from Pavia, Italy, and Tromsø, Norway, was performed. AMR-associated and virulence genes were detected, and the population distribution of these genes was studied. FINDINGS: Of 957 samples collected around Tamale, Ghana, 620 were positive for Klebsiella spp. 573 Klebsiella isolates were successfully sequenced, of which 370 were Klebsiella pneumoniae. Only two hospital isolates were carbapenem-resistant. Extended-spectrum ß-lactamase (ESBL) genes were relatively common among the Ghanaian clinical isolates but rare in the environmental samples. Prevalence of ESBL genes in human-hospital disease samples was 64% (14 of 22 isolates) in Ghana and 44% (four of nine isolates) in Italy, and prevalence in human-hospital carriage samples was 7% (eight of 107) in Ghana and 13% (54 of 428) in Italy; the prevalence was higher in human-hospital disease samples than in human-hospital carriage samples in both countries, and prevalence across both samples in both countries was higher than in Norway. Ghanaian isolates showed evidence of high recombination rates (recombination events compared with point mutations [r/m] 9·455) and a considerable accessory gene overlap with isolates from Italy and Norway. INTERPRETATION: Although several AMR-associated gene classes were observed relatively frequently in non-clinical sources, ESBL, carbapenemase, and virulence genes were predominantly present only in hospital samples. These results suggest that interventions should be focused on clinical settings to have the greatest effect on the prevalence and dissemination of AMR-associated genes. FUNDING: European Research Council (742158), Academy of Finland EuroHPC grant, Trond Mohn Foundation (BATTALION grant), and Wellcome Trust.

Strong pathogen competition in neonatal gut colonisation.

Opportunistic bacterial pathogen species and their strains that colonise the human gut are generally understood to compete against both each other and the commensal species colonising this ecosystem. Currently we are lacking a population-wide quantification of strain-level colonisation dynamics and the relationship of colonisation potential to prevalence in disease, and how ecological factors might be modulating these. Here, using a combination of latest high-resolution metagenomics and strain-level genomic epidemiology methods we performed a characterisation of the competition and colonisation dynamics for a longitudinal cohort of neonatal gut microbiomes. We found strong inter- and intra-species competition dynamics in the gut colonisation process, but also a number of synergistic relationships among several species belonging to genus Klebsiella, which includes the prominent human pathogen Klebsiella pneumoniae. No evidence of preferential colonisation by hospital-adapted pathogen lineages in either vaginal or caesarean section birth groups was detected. Our analysis further enabled unbiased assessment of strain-level colonisation potential of extra-intestinal pathogenic Escherichia coli (ExPEC) in comparison with their propensity to cause bloodstream infections. Our study highlights the importance of systematic surveillance of bacterial gut pathogens, not only from disease but also from carriage state, to better inform therapies and preventive medicine in the future.

Repeated out-of-Africa expansions of Helicobacter pylori driven by replacement of deleterious mutations.

Helicobacter pylori lives in the human stomach and has a population structure resembling that of its host. However, H. pylori from Europe and the Middle East trace substantially more ancestry from modern African populations than the humans that carry them. Here, we use a collection of Afro-Eurasian H. pylori genomes to show that this African ancestry is due to at least three distinct admixture events. H. pylori from East Asia, which have undergone little admixture, have accumulated many more non-synonymous mutations than African strains. European and Middle Eastern bacteria have elevated African ancestry at the sites of these mutations, implying selection to remove them during admixture. Simulations show that population fitness can be restored after bottlenecks by migration and subsequent admixture of small numbers of bacteria from non-bottlenecked populations. We conclude that recent spread of African DNA has been driven by deleterious mutations accumulated during the original out-of-Africa bottleneck.

A large-scale genomic snapshot of Klebsiella spp. isolates in Northern Italy reveals limited transmission between clinical and non-clinical settings.

The Klebsiella group, found in humans, livestock, plants, soil, water and wild animals, is genetically and ecologically diverse. Many species are opportunistic pathogens and can harbour diverse classes of antimicrobial resistance genes. Healthcare-associated Klebsiella pneumoniae clones that are non-susceptible to carbapenems can spread rapidly, representing a high public health burden. Here we report an analysis of 3,482 genome sequences representing 15 Klebsiella species sampled over a 17-month period from a wide range of clinical, community, animal and environmental settings in and around the Italian city of Pavia. Northern Italy is a hotspot for hospital-acquired carbapenem non-susceptible Klebsiella and thus a pertinent setting to examine the overlap between isolates in clinical and non-clinical settings. We found no genotypic or phenotypic evidence for non-susceptibility to carbapenems outside the clinical environment. Although we noted occasional transmission between clinical and non-clinical settings, our data point to a limited role of animal and environmental reservoirs in the human acquisition of Klebsiella spp. We also provide a detailed genus-wide view of genomic diversity and population structure, including the identification of new groups.

Quantifying bacterial evolution in the wild: A birthday problem for Campylobacter lineages.

Measuring molecular evolution in bacteria typically requires estimation of the rate at which nucleotide changes accumulate in strains sampled at different times that share a common ancestor. This approach has been useful for dating ecological and evolutionary events that coincide with the emergence of important lineages, such as outbreak strains and obligate human pathogens. However, in multi-host (niche) transmission scenarios, where the pathogen is essentially an opportunistic environmental organism, sampling is often sporadic and rarely reflects the overall population, particularly when concentrated on clinical isolates. This means that approaches that assume recent common ancestry are not applicable. Here we present a new approach to estimate the molecular clock rate in Campylobacter that draws on the popular probability conundrum known as the 'birthday problem'. Using large genomic datasets and comparative genomic approaches, we use isolate pairs that share recent common ancestry to estimate the rate of nucleotide change for the population. Identifying synonymous and non-synonymous nucleotide changes, both within and outside of recombined regions of the genome, we quantify clock-like diversification to estimate synonymous rates of nucleotide change for the common pathogenic bacteria Campylobacter coli (2.4 x 10-6 s/s/y) and Campylobacter jejuni (3.4 x 10-6 s/s/y). Finally, using estimated total rates of nucleotide change, we infer the number of effective lineages within the sample time frame-analogous to a shared birthday-and assess the rate of turnover of lineages in our sample set over short evolutionary timescales. This provides a generalizable approach to calibrating rates in populations of environmental bacteria and shows that multiple lineages are maintained, implying that large-scale clonal sweeps may take hundreds of years or more in these species.

Genome of Superficieibacter maynardsmithii, a novel, antibiotic susceptible representative of Enterobacteriaceae.

During a citywide microbiological screening project in Pavia (Italy) a bacterial strain isolated from the surface of an Automated Teller Machine was classified as a Klebsiella sp. by MALDI-TOF spectrometry, and shown to be susceptible to the most antimicrobial classes by phenotypic testing. After Illumina genome sequencing and subsequent assembly, a high-quality draft genome was obtained (size = 5,051,593 bp, N50 = 615,571 bp, largest contig = 1,328,029 bp, N_contig = 17, GC content = 51.58%, coverage = 141.42), absence of antimicrobial resistance genes was confirmed, but the strain resulted to be highly divergent from all Klebsiella, and more related to other Enterobacteriaceae. The higher values of 16S rRNA identity were with members of the genera Citrobacter, Salmonella, and "Superficieibacter." An ortholog-based phylogenomic analysis indicated a sister group relationship with "Superficieibacter electus," in a distinct clade from other members of the Enterobacteriaceae family. In order to evaluate whether the novel genome represents a new species of "Superficiebacter," average nucleotide identity (ANI) and Hadamard analysis were performed on a dataset of 78 Enterobacteriaceae. The novel genome showed an ANI of 87.51% with S. electus, which compared on identity values between other members of the family, clearly indicates that the genome represents a new species within the genus "Superficieibacter." We propose for the new species the name "Superficieibacter maynardsmithii."

A high prevalence of blaOXA-48 in Klebsiella (Raoultella) ornithinolytica and related species in hospital wastewater in South West England.

Klebsiella species occupy a wide range of environmental and animal niches, and occasionally cause opportunistic infections that are resistant to multiple antibiotics. In particular, Klebsiella pneumoniae (Kpne) has gained notoriety as a major nosocomial pathogen, due principally to the rise in non-susceptibility to carbapenems and other beta-lactam antibiotics. Whilst it has been proposed that the urban water cycle facilitates transmission of pathogens between clinical settings and the environment, the level of risk posed by resistant Klebsiella strains in hospital wastewater remains unclear. We used whole genome sequencing (WGS) to compare Klebsiella species in contemporaneous samples of wastewater from an English hospital and influent to the associated wastewater treatment plant (WWTP). As we aimed to characterize representative samples of Klebsiella communities, we did not actively select for antibiotic resistance (other than for ampicillin), nor for specific Klebsiella species. Two species, Kpne and K. (Raoultella) ornithinolytica (Korn), were of equal dominance in the hospital wastewater, and four other Klebsiella species were present in low abundance in this sample. In contrast, despite being the species most closely associated with healthcare settings, Kpne was the dominant species within the WWTP influent. In total, 29 % of all isolates harboured the blaOXA-48 gene on a pOXA-48-like plasmid, and these isolates were almost exclusively recovered from the hospital wastewater. This gene was far more common in Korn (68 % of isolates) than in Kpne (3.4 % of isolates). In general plasmid-borne, but not chromosomal, resistance genes were significantly enriched in the hospital wastewater sample. These data implicate hospital wastewater as an important reservoir for antibiotic-resistant Klebsiella, and point to an unsuspected role of species within the Raoultella group in the maintenance and dissemination of plasmid-borne blaOXA-48. This article contains data hosted by Microreact.

Description of Klebsiella spallanzanii sp. nov. and of Klebsiella pasteurii sp. nov.

Klebsiella oxytoca causes opportunistic human infections and post-antibiotic haemorrhagic diarrhea. This Enterobacteriaceae species is genetically heterogeneous and is currently subdivided into seven phylogroups (Ko1 to Ko4 and Ko6 to Ko8). Here we investigated the taxonomic status of phylogroups Ko3 and Ko4. Genomic sequence-based phylogenetic analyses demonstrate that Ko3 and Ko4 formed well-defined sequence clusters related to, but distinct from, Klebsiella michiganensis (Ko1), K. oxytoca (Ko2), K. huaxiensis (Ko8), and K. grimontii (Ko6). The average nucleotide identity (ANI) of Ko3 and Ko4 were 90.7% with K. huaxiensis and 95.5% with K. grimontii, respectively. In addition, three strains of K. huaxiensis, a species so far described based on a single strain from a urinary tract infection patient in China, were isolated from cattle and human feces. Biochemical and MALDI-ToF mass spectrometry analysis allowed differentiating Ko3, Ko4, and Ko8 from the other K. oxytoca species. Based on these results, we propose the names Klebsiella spallanzanii for the Ko3 phylogroup, with SPARK_775_C1T (CIP 111695T and DSM 109531T) as type strain, and Klebsiella pasteurii for Ko4, with SPARK_836_C1T (CIP 111696T and DSM 109530T) as type strain. Strains of K. spallanzanii were isolated from human urine, cow feces, and farm surfaces, while strains of K. pasteurii were found in fecal carriage from humans, cows, and turtles.

PIRATE: A fast and scalable pangenomics toolbox for clustering diverged orthologues in bacteria.

BACKGROUND: Cataloguing the distribution of genes within natural bacterial populations is essential for understanding evolutionary processes and the genetic basis of adaptation. Advances in whole genome sequencing technologies have led to a vast expansion in the amount of bacterial genomes deposited in public databases. There is a pressing need for software solutions which are able to cluster, catalogue and characterise genes, or other features, in increasingly large genomic datasets. RESULTS: Here we present a pangenomics toolbox, PIRATE (Pangenome Iterative Refinement and Threshold Evaluation), which identifies and classifies orthologous gene families in bacterial pangenomes over a wide range of sequence similarity thresholds. PIRATE builds upon recent scalable software developments to allow for the rapid interrogation of thousands of isolates. PIRATE clusters genes (or other annotated features) over a wide range of amino acid or nucleotide identity thresholds and uses the clustering information to rapidly identify paralogous gene families and putative fission/fusion events. Furthermore, PIRATE orders the pangenome using a directed graph, provides a measure of allelic variation, and estimates sequence divergence for each gene family. CONCLUSIONS: We demonstrate that PIRATE scales linearly with both number of samples and computation resources, allowing for analysis of large genomic datasets, and compares favorably to other popular tools. PIRATE provides a robust framework for analysing bacterial pangenomes, from largely clonal to panmictic species.

Whole genome sequencing and phylogenetic analysis of strains of the agent of Lyme disease Borrelia burgdorferi from Canadian emergence zones.

Lyme disease is emerging in southern Canada due to range expansion of the tick vector, followed by invasion of the agent of Lyme disease Borrelia burgdorferi sensu stricto. Strain diversity, as determined by Multi Locus Sequence Typing, occurs in this zone of emergence, and this may have its origins in adaptation to ecological niches, and have phenotypic consequences for pathogenicity and serological test performance. Sixty-four unique strains were cultured from ticks collected in southern Canada and the genomes sequenced using the Illumina MiSeq platform. A maximum likelihood phylogenetic tree of the chromosome revealed two large clades with multiple subclades. Consistent with previous studies on this species, the clades were not geographically defined, and some Canadian strains were highly divergent from previously sequenced US strains. There was evidence for recombination in the chromosome but this did not affect the phylogeny. Analysis of chromosomal genes indicated that these are under intense purifying selection. Phylogenies of the accessory genome and chromosome were congruent. Therefore strain differences identified in the phylogeny of chromosomal genes likely act as a proxy for genetic determinants of phenotypic differences amongst strains that are harboured in the accessory genome. Further studies on health implications of strain diversity are needed.

Piggy: a rapid, large-scale pan-genome analysis tool for intergenic regions in bacteria.

Background: The concept of the "pan-genome," which refers to the total complement of genes within a given sample or species, is well established in bacterial genomics. Rapid and scalable pipelines are available for managing and interpreting pan-genomes from large batches of annotated assemblies. However, despite overwhelming evidence that variation in intergenic regions in bacteria can directly influence phenotypes, most current approaches for analyzing pan-genomes focus exclusively on protein-coding sequences. Findings: To address this we present Piggy, a novel pipeline that emulates Roary except that it is based only on intergenic regions. A key utility provided by Piggy is the detection of highly divergent ("switched") intergenic regions (IGRs) upstream of genes. We demonstrate the use of Piggy on large datasets of clinically important lineages of Staphylococcus aureus and Escherichia coli. Conclusions: For S. aureus, we show that highly divergent (switched) IGRs are associated with differences in gene expression and we establish a multilocus reference database of IGR alleles (igMLST; implemented in BIGSdb).

The use of Oxford Nanopore native barcoding for complete genome assembly.

Background: The Oxford Nanopore Technologies MinION(TM) is a mobile DNA sequencer that can produce long read sequences with a short turn-around time. Here we report the first demonstration of single contig genome assembly using Oxford Nanopore native barcoding when applied to a multiplexed library of 12 samples and combined with existing Illumina short read data. This paves the way for the closure of multiple bacterial genomes from a single MinION(TM) sequencing run, given the availability of existing short read data. The strain we used, MHO_001, represents the important community-acquired methicillin-resistant Staphylococcus aureus lineage USA300. Findings: Using a hybrid assembly of existing short read and barcoded long read sequences from multiplexed data, we completed a genome of the S. aureus USA300 strain MHO_001. The long read data represented only ∼5% to 10% of an average MinION(TM) run (∼7x genomic coverage), but, using standard tools, this was sufficient to complete the circular chromosome of S. aureus strain MHO_001 (2.86 Mb) and two complete plasmids (27 Kb and 3 Kb). Minor differences were noted when compared to USA300 reference genome, USA300_FPR3757, including the translocation, loss, and gain of mobile genetic elements. Conclusion: Here we demonstrate that MinION(TM) reads, multiplexed using native barcoding, can be used in combination with short read data to fully complete a bacterial genome. The ability to complete multiple genomes, for which short read data is already available, from a single MinION(TM) run is set to impact our understanding of accessory genome content, plasmid diversity, and genome rearrangements.

Comparative Analyses of Selection Operating on Nontranslated Intergenic Regions of Diverse Bacterial Species.

Nontranslated intergenic regions (IGRs) compose 10-15% of bacterial genomes, and contain many regulatory elements with key functions. Despite this, there are few systematic studies on the strength and direction of selection operating on IGRs in bacteria using whole-genome sequence data sets. Here we exploit representative whole-genome data sets from six diverse bacterial species: Staphylococcus aureus, Streptococcus pneumoniae, Mycobacterium tuberculosis, Salmonella enterica, Klebsiella pneumoniae, and Escherichia coli We compare patterns of selection operating on IGRs using two independent methods: the proportion of singleton mutations and the dI/dS ratio, where dI is the number of intergenic SNPs per intergenic site. We find that the strength of purifying selection operating over all intergenic sites is consistently intermediate between that operating on synonymous and nonsynonymous sites. Ribosome binding sites and noncoding RNAs tend to be under stronger selective constraint than promoters and Rho-independent terminators. Strikingly, a clear signal of purifying selection remains even when all these major categories of regulatory elements are excluded, and this constraint is highest immediately upstream of genes. While a paucity of variation means that the data for M. tuberculosis are more equivocal than for the other species, we find strong evidence for positive selection within promoters of this species. This points to a key adaptive role for regulatory changes in this important pathogen. Our study underlines the feasibility and utility of gauging the selective forces operating on bacterial IGRs from whole-genome sequence data, and suggests that our current understanding of the functionality of these sequences is far from complete.

Transmission of Staphylococcus aureus from Humans to Green Monkeys in The Gambia as Revealed by Whole-Genome Sequencing.

UNLABELLED: Staphylococcus aureus is an important pathogen of humans and animals. We genome sequenced 90 S. aureus isolates from The Gambia: 46 isolates from invasive disease in humans, 13 human carriage isolates, and 31 monkey carriage isolates. We inferred multiple anthroponotic transmissions of S. aureus from humans to green monkeys (Chlorocebus sabaeus) in The Gambia over different time scales. We report a novel monkey-associated clade of S. aureus that emerged from a human-to-monkey switch estimated to have occurred 2,700 years ago. Adaptation of this lineage to the monkey host is accompanied by the loss of phage-carrying genes that are known to play an important role in human colonization. We also report recent anthroponotic transmission of the well-characterized human lineages sequence type 6 (ST6) and ST15 to monkeys, probably because of steadily increasing encroachment of humans into the monkeys' habitat. Although we have found no evidence of transmission of S. aureus from monkeys to humans, as the two species come into ever-closer contact, there might be an increased risk of additional interspecies exchanges of potential pathogens. IMPORTANCE: The population structures of Staphylococcus aureus in humans and monkeys in sub-Saharan Africa have been previously described using multilocus sequence typing (MLST). However, these data lack the power to accurately infer details regarding the origin and maintenance of new adaptive lineages. Here, we describe the use of whole-genome sequencing to detect transmission of S. aureus between humans and nonhuman primates and to document the genetic changes accompanying host adaptation. We note that human-to-monkey switches tend to be more common than the reverse and that a novel monkey-associated clade is likely to have emerged from such a switch approximately 2,700 years ago. Moreover, analysis of the accessory genome provides important clues as to the genetic changes underpinning host adaptation and, in particular, shows that human-to-monkey switches tend to be associated with the loss of genes known to confer adaptation to the human host.

The Stealthy Superbug: the Role of Asymptomatic Enteric Carriage in Maintaining a Long-Term Hospital Outbreak of ST228 Methicillin-Resistant Staphylococcus aureus.

UNLABELLED: Whole-genome sequencing (WGS) of 228 isolates was used to elucidate the origin and dynamics of a long-term outbreak of methicillin-resistant Staphylococcus aureus (MRSA) sequence type 228 (ST228) SCCmec I that involved 1,600 patients in a tertiary care hospital between 2008 and 2012. Combining of the sequence data with detailed metadata on patient admission and movement confirmed that the outbreak was due to the transmission of a single clonal variant of ST228, rather than repeated introductions of this clone into the hospital. We note that this clone is significantly more frequently recovered from groin and rectal swabs than other clones (P < 0.0001) and is also significantly more transmissible between roommates (P < 0.01). Unrecognized MRSA carriers, together with movements of patients within the hospital, also seem to have played a major role. These atypical colonization and transmission dynamics can help explain how the outbreak was maintained over the long term. This "stealthy" asymptomatic colonization of the gut, combined with heightened transmissibility (potentially reflecting a role for environmental reservoirs), means the dynamics of this outbreak share some properties with enteric pathogens such as vancomycin-resistant enterococci or Clostridium difficile. IMPORTANCE: Using whole-genome sequencing, we showed that a large and prolonged outbreak of methicillin-resistant Staphylococcus aureus was due to the clonal spread of a specific strain with genetic elements adapted to the hospital environment. Unrecognized MRSA carriers, the movement of patients within the hospital, and the low detection with clinical specimens were also factors that played a role in this occurrence. The atypical colonization of the gut means the dynamics of this outbreak may share some properties with enteric pathogens.

Evolutionary Trade-Offs Underlie the Multi-faceted Virulence of Staphylococcus aureus.

Bacterial virulence is a multifaceted trait where the interactions between pathogen and host factors affect the severity and outcome of the infection. Toxin secretion is central to the biology of many bacterial pathogens and is widely accepted as playing a crucial role in disease pathology. To understand the relationship between toxicity and bacterial virulence in greater depth, we studied two sequenced collections of the major human pathogen Staphylococcus aureus and found an unexpected inverse correlation between bacterial toxicity and disease severity. By applying a functional genomics approach, we identified several novel toxicity-affecting loci responsible for the wide range in toxic phenotypes observed within these collections. To understand the apparent higher propensity of low toxicity isolates to cause bacteraemia, we performed several functional assays, and our findings suggest that within-host fitness differences between high- and low-toxicity isolates in human serum is a contributing factor. As invasive infections, such as bacteraemia, limit the opportunities for onward transmission, highly toxic strains could gain an additional between-host fitness advantage, potentially contributing to the maintenance of toxicity at the population level. Our results clearly demonstrate how evolutionary trade-offs between toxicity, relative fitness, and transmissibility are critical for understanding the multifaceted nature of bacterial virulence.