Extended-spectrum β-lactamase (ESBL)-producing bacteria causing clinical infections are often also multidrug-resistant (MDR; resistance to ≥3 antimicrobial drug classes), therefore treatment options may be limited. High carriage rates of these potentially zoonotic bacteria have been found in livestock and companion animals. Therefore, people working in veterinary hospitals may be a high-risk population for carriage. This is the first study to determine the prevalence and longitudinal carriage of antimicrobial-resistant (AMR) and ESBL-producing faecal Escherichia coli in veterinary hospital staff and students. Prevalence of faecal AMR and ESBL-producing E coli was determined in 84 staff members and students in three UK veterinary hospitals. Twenty-seven participants were followed for six weeks to investigate longitudinal carriage. Antimicrobial susceptibility and phenotypic ESBL production were determined and selected isolates were whole genome sequenced. ESBL-producing E coli were isolated from five participants (5.95 per cent; 95 per cent CI 0.89 to 11.0 per cent); two participants carried ESBL-producing E coli resistant to all antimicrobials tested. Carriage of MDR E coli was common (32.1 per cent; 95per cent CI 22.2 to 42.1 per cent) and there was a high prevalence of ciprofloxacin resistance (11.9 per cent; 95 per cent CI 4.98 to 18.8 per cent). ESBL-producing E coli were isolated from seven longitudinal participants (25.9 per cent; 95 per cent CI 9.40 to 42.5 per cent); two participants carried ESBL-producing E coli for the entire study period. Twenty-six participants (96.3 per cent; 95 per cent CI 89.2 to 100) carried ≥1 MDR E coli isolate during the six-week period, with seven participants (25.9 per cent) carrying ≥1 MDR isolate for at least five out of six weeks. The prevalence of faecal ESBL-producing E coli in cross-sectional participants is similar to asymptomatic general populations. However, much higher levels of carriage were observed longitudinally in participants. It is vital that veterinary hospitals implement gold-standard biosecurity to prevent transmission of MDR and ESBL-producing bacteria between patients and staff. Healthcare providers should be made aware that people working in veterinary hospitals are a high-risk population for carriage of MDR and ESBL-producing bacteria, and that this poses a risk to the carrier and for transmission of resistance throughout the wider community.

Campylobacter jejuni is the primary cause of bacterial gastroenteritis worldwide, infecting humans mostly through consumption of contaminated poultry. C. jejuni is common in the gut of wild birds, and shows distinct strain-specific association to particular bird species. This contrasts with farm animals, in which several genotypes co-exist. It is unclear if the barriers restricting transmission between host species of such specialist strains are related to environmental factors such as contact between host species, bacterial survival in the environment, etc., or rather to strain specific adaptation to the intestinal environment of specific hosts. We compared colonization dynamics in vivo between two host-specific C. jejuni from a song thrush (ST-1304 complex) and a mallard (ST-995), and a generalist strain from chicken (ST-21 complex) in a wild host, the mallard (Anas platyrhynchos). In 18-days infection experiments, the song thrush strain showed only weak colonization and was cleared from all birds after 10 days, whereas both mallard and chicken strains remained stable. When the chicken strain was given 4 days prior to co-infection of the same birds with a mallard strain, it was rapidly outcompeted by the latter. In contrast, when the mallard strain was given 4 days prior to co-infection with the chicken strain, the mallard strain remained and expansion of the chicken strain was delayed. Our results suggest strain-specific differences in the ability of C. jejuni to colonize mallards, likely associated with host origin. This difference might explain observed host association patterns in C. jejuni from wild birds.

Reference and type strains of well-known bacteria have been a cornerstone of microbiology research for decades. The sharing of well-characterized isolates among laboratories has run in parallel with research efforts and enhanced the reproducibility of experiments, leading to a wealth of knowledge about trait variation in different species and the underlying genetics. Campylobacter jejuni strain NCTC 11168, deposited at the National Collection of Type Cultures in 1977, has been adopted widely as a reference strain by researchers worldwide and was the first Campylobacter for which the complete genome was published (in 2000). In this study, we collected 23 C. jejuni NCTC 11168 reference isolates from laboratories across the UK and compared variation in simple laboratory phenotypes with genetic variation in sequenced genomes. Putatively identical isolates, identified previously to have aberrant phenotypes, varied by up to 281 SNPs (in 15 genes) compared to the most recent reference strain. Isolates also display considerable phenotype variation in motility, morphology, growth at 37 °C, invasion of chicken and human cell lines, and susceptibility to ampicillin. This study provides evidence of ongoing evolutionary change among C. jejuni isolates as they are cultured in different laboratories and highlights the need for careful consideration of genetic variation within laboratory reference strains. This article contains data hosted by Microreact.

The use of antimicrobials in human and veterinary medicine has coincided with a rise in antimicrobial resistance (AMR) in the food-borne pathogens Campylobacter jejuni and Campylobacter coli. Faecal contamination from the main reservoir hosts (livestock, especially poultry) is the principal route of human infection but little is known about the spread of AMR among source and sink populations. In particular, questions remain about how Campylobacter resistomes interact between species and hosts, and the potential role of sewage as a conduit for the spread of AMR. Here, we investigate the genomic variation associated with AMR in 168 C. jejuni and 92 C. coli strains isolated from humans, livestock and urban effluents in Spain. AMR was tested in vitro and isolate genomes were sequenced and screened for putative AMR genes and alleles. Genes associated with resistance to multiple drug classes were observed in both species and were commonly present in multidrug-resistant genomic islands (GIs), often located on plasmids or mobile elements. In many cases, these loci had alleles that were shared among C. jejuni and C. coli consistent with horizontal transfer. Our results suggest that specific antibiotic resistance genes have spread among Campylobacter isolated from humans, animals and the environment.

The genetic structure of bacterial populations can be related to geographical locations of isolation. In some species, there is a strong correlation between geographical distance and genetic distance, which can be caused by different evolutionary mechanisms. Patterns of ancient admixture in Helicobacter pylori can be reconstructed in concordance with past human migration, whereas in Mycobacterium tuberculosis it is the lack of recombination that causes allopatric clusters. In Campylobacter, analyses of genomic data and molecular typing have been successful in determining the reservoir host species, but not geographical origin. We investigated biogeographical variation in highly recombining genes to determine the extent of clustering between genomes from geographically distinct Campylobacter populations. Whole-genome sequences from 294 Campylobacter isolates from North America and the UK were analysed. Isolates from within the same country shared more recently recombined DNA than isolates from different countries. Using 15 UK/American closely matched pairs of isolates that shared ancestors, we identify regions that have frequently and recently recombined to test their correlation with geographical origin. The seven genes that demonstrated the greatest clustering by geography were used in an attribution model to infer geographical origin which was tested using a further 383 UK clinical isolates to detect signatures of recent foreign travel. Patient records indicated that in 46 cases, travel abroad had occurred <2 weeks prior to sampling, and genomic analysis identified that 34 (74%) of these isolates were of a non-UK origin. Identification of biogeographical markers in Campylobacter genomes will contribute to improved source attribution of clinical Campylobacter infection and inform intervention strategies to reduce campylobacteriosis.

Chickens are the most common birds on Earth and colibacillosis is among the most common diseases affecting them. This major threat to animal welfare and safe sustainable food production is difficult to combat because the etiological agent, avian pathogenic Escherichia coli (APEC), emerges from ubiquitous commensal gut bacteria, with no single virulence gene present in all disease-causing isolates. Here, we address the underlying evolutionary mechanisms of extraintestinal spread and systemic infection in poultry. Combining population scale comparative genomics and pangenome-wide association studies, we compare E. coli from commensal carriage and systemic infections. We identify phylogroup-specific and species-wide genetic elements that are enriched in APEC, including pathogenicity-associated variation in 143 genes that have diverse functions, including genes involved in metabolism, lipopolysaccharide synthesis, heat shock response, antimicrobial resistance and toxicity. We find that horizontal gene transfer spreads pathogenicity elements, allowing divergent clones to cause infection. Finally, a Random Forest model prediction of disease status (carriage vs. disease) identifies pathogenic strains in the emergent ST-117 poultry-associated lineage with 73% accuracy, demonstrating the potential for early identification of emergent APEC in healthy flocks.

Campylobacter is the most common cause of bacterial gastroenteritis worldwide, and diarrhoeal disease is a major cause of child morbidity, growth faltering and mortality in low- and middle-income countries. Despite evidence of high incidence and differences in disease epidemiology, there is limited genomic data from studies in developing countries. In this study, we aimed to quantify the extent of gene sharing in local and global populations. We characterized the genetic diversity and accessory-genome content of a collection of Campylobacter isolates from the Cairo metropolitan area, Egypt. In total, 112 Campylobacter isolates were collected from broiler carcasses (n=31), milk and dairy products (n=24), and patients suffering from gastroenteritis (n=57). Among the most common sequence types (STs), we identified the globally disseminated host generalist ST-21 clonal complex (CC21) and the poultry specialists CC206, CC464 and CC48. Notably, CC45 and the cattle-specialist CC42 were under-represented, with a total absence of CC61. Core- and accessory-genome sharing was compared among isolates from Egypt and a comparable collection from the UK (Oxford). Lineage-specific accessory-genome sharing was significantly higher among isolates from the same country, particularly CC21, which demonstrated greater local geographical clustering. In contrast, no geographical clustering was noted in either the core or accessory genome of CC828, suggesting a highly admixed population. A greater proportion of Campylobacter coli isolates were multidrug resistant compared to Campylobacter jejuni . Our results suggest that there is more horizontal transfer of accessory genes between strains in Egypt. This has strong implications for controlling the spread of antimicrobial resistance among this important pathogen.

Listeria monocytogenes is an opportunistic food-borne bacterium that is capable of infecting humans with high rates of hospitalization and mortality. Natural populations are genotypically and phenotypically variable, with some lineages being responsible for most human infections. The success of L. monocytogenes is linked to its capacity to persist on food and in the environment. Biofilms are an important feature that allow these bacteria to persist and infect humans, so understanding the genetic basis of biofilm formation is key to understanding transmission. We sought to investigate the biofilm-forming ability of L. monocytogenes by identifying genetic variation that underlies biofilm formation in natural populations using genome-wide association studies (GWAS). Changes in gene expression of specific strains during biofilm formation were then investigated using RNA sequencing (RNA-seq). Genetic variation associated with enhanced biofilm formation was identified in 273 genes by GWAS and differential expression in 220 genes by RNA-seq. Statistical analyses show that the number of overlapping genes flagged by either type of experiment is less than expected by random sampling. This novel finding is consistent with an evolutionary scenario where rapid adaptation is driven by variation in gene expression of pioneer genes, and this is followed by slower adaptation driven by nucleotide changes within the core genome.

The increasing availability of hundreds of whole bacterial genomes provides opportunities for enhanced understanding of the genes and alleles responsible for clinically important phenotypes and how they evolved. However, it is a significant challenge to develop easy-to-use and scalable methods for characterizing these large and complex data and relating it to disease epidemiology. Existing approaches typically focus on either homologous sequence variation in genes that are shared by all isolates, or non-homologous sequence variation--focusing on genes that are differentially present in the population. Here we present a comparative genomics approach that simultaneously approximates core and accessory genome variation in pathogen populations and apply it to pathogenic species in the genus Campylobacter. A total of 7 published Campylobacter jejuni and Campylobacter coli genomes were selected to represent diversity across these species, and a list of all loci that were present at least once was compiled. After filtering duplicates a 7-isolate reference pan-genome, of 3,933 loci, was defined. A core genome of 1,035 genes was ubiquitous in the sample accounting for 59% of the genes in each isolate (average genome size of 1.68 Mb). The accessory genome contained 2,792 genes. A Campylobacter population sample of 192 genomes was screened for the presence of reference pan-genome loci with gene presence defined as a BLAST match of ≥ 70% identity over ≥ 50% of the locus length--aligned using MUSCLE on a gene-by-gene basis. A total of 21 genes were present only in C. coli and 27 only in C. jejuni, providing information about functional differences associated with species and novel epidemiological markers for population genomic analyses. Homologs of these genes were found in several of the genomes used to define the pan-genome and, therefore, would not have been identified using a single reference strain approach.

Helicobacter pylori is a common component of the human stomach microbiota, possibly dating back to the speciation of Homo sapiens. A history of pathogen evolution in allopatry has led to the development of genetically distinct H. pylori subpopulations, associated with different human populations, and more recent admixture among H. pylori subpopulations can provide information about human migrations. However, little is known about the degree to which some H. pylori genes are conserved in the face of admixture, potentially indicating host adaptation, or how virulence genes spread among different populations. We analyzed H. pylori genomes from 14 countries in the Americas, strains from the Iberian Peninsula, and public genomes from Europe, Africa, and Asia, to investigate how admixture varies across different regions and gene families. Whole-genome analyses of 723 H. pylori strains from around the world showed evidence of frequent admixture in the American strains with a complex mosaic of contributions from H. pylori populations originating in the Americas as well as other continents. Despite the complex admixture, distinctive genomic fingerprints were identified for each region, revealing novel American H. pylori subpopulations. A pan-genome Fst analysis showed that variation in virulence genes had the strongest fixation in America, compared with non-American populations, and that much of the variation constituted non-synonymous substitutions in functional domains. Network analyses suggest that these virulence genes have followed unique evolutionary paths in the American populations, spreading into different genetic backgrounds, potentially contributing to the high risk of gastric cancer in the region.

Helicobacter pylori genetic diversity is known to be influenced by mobile genomic elements. Here we focused on prophages, the least characterized mobile elements of H. pylori. We present the full genomic sequences, insertion sites and phylogenetic analysis of 28 prophages found in H. pylori isolates from patients of distinct disease types, ranging from gastritis to gastric cancer, and geographic origins, covering most continents. The genome sizes of these prophages range from 22.6-33.0 Kbp, consisting of 27-39 open reading frames. A 36.6% GC was found in prophages in contrast to 39% in H. pylori genome. Remarkably a conserved integration site was found in over 50% of the cases. Nearly 40% of the prophages harbored insertion sequences (IS) previously described in H. pylori. Tandem repeats were frequently found in the intergenic region between the prophage at the 3' end and the bacterial gene. Furthermore, prophage genomes present a robust phylogeographic pattern, revealing four distinct clusters: one African, one Asian and two European prophage populations. Evidence of recombination was detected within the genome of some prophages, resulting in genome mosaics composed by different populations, which may yield additional H. pylori phenotypes.

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.

Horizontal gene transfer (HGT) can allow traits that have evolved in one bacterial species to transfer to another. This has potential to rapidly promote new adaptive trajectories such as zoonotic transfer or antimicrobial resistance. However, for this to occur requires gaps to align in barriers to recombination within a given time frame. Chief among these barriers is the physical separation of species with distinct ecologies in separate niches. Within the genus Campylobacter, there are species with divergent ecologies, from rarely isolated single-host specialists to multihost generalist species that are among the most common global causes of human bacterial gastroenteritis. Here, by characterizing these contrasting ecologies, we can quantify HGT among sympatric and allopatric species in natural populations. Analyzing recipient and donor population ancestry among genomes from 30 Campylobacter species, we show that cohabitation in the same host can lead to a six-fold increase in HGT between species. This accounts for up to 30% of all SNPs within a given species and identifies highly recombinogenic genes with functions including host adaptation and antimicrobial resistance. As described in some animal and plant species, ecological factors are a major evolutionary force for speciation in bacteria and changes to the host landscape can promote partial convergence of distinct species through HGT.

Bacterial clades are often ecologically distinct, despite extensive horizontal gene transfer (HGT). How selection works on different parts of bacterial pan-genomes to drive and maintain the emergence of clades is unclear. Focusing on the three largest clades in the diverse and well-studied Bacillus cereus sensu lato group, we identified clade-specific core genes (present in all clade members) and then used clade-specific allelic diversity to identify genes under purifying and diversifying selection. Clade-specific accessory genes (present in a subset of strains within a clade) were characterized as being under selection using presence/absence in specific clades. Gene ontology analyses of genes under selection revealed that different gene functions were enriched in different clades. Furthermore, some gene functions were enriched only amongst clade-specific core or accessory genomes. Genes under purifying selection were often clade-specific, while genes under diversifying selection showed signs of frequent HGT. These patterns are consistent with different selection pressures acting on both the core and the accessory genomes of different clades and can lead to ecological divergence in both cases. Examining variation in allelic diversity allows us to uncover genes under clade-specific selection, allowing ready identification of strains and their ecological niche.

Determining the composition of bacterial communities beyond the level of a genus or species is challenging because of the considerable overlap between genomes representing close relatives. Here, we present the mSWEEP pipeline for identifying and estimating the relative sequence abundances of bacterial lineages from plate sweeps of enrichment cultures. mSWEEP leverages biologically grouped sequence assembly databases, applying probabilistic modelling, and provides controls for false positive results. Using sequencing data from major pathogens, we demonstrate significant improvements in lineage quantification and detection accuracy. Our pipeline facilitates investigating cultures comprising mixtures of bacteria, and opens up a new field of plate sweep metagenomics.

Streptococcus suis is a leading cause of bacterial meningitis in South-East Asia, with frequent zoonotic transfer to humans associated with close contact with pigs. A small number of invasive lineages are responsible for endemic infection in the swine industry, causing considerable global economic losses. A lack of surveillance and a rising trend in clinical treatment failure has raised concerns of growing antimicrobial resistance (AMR) among invasive S. suis . Gene flow between healthy and disease isolates is poorly understood and, in this study, we sample and sequence a collection of isolates predominantly from healthy pigs in Chiang Mai province, Northern Thailand. Pangenome characterization identified extensive genetic diversity and frequent AMR carriage in isolates from healthy pigs. Multiple AMR genes were identified, conferring resistance to aminoglycosides, lincosamides, tetracycline and macrolides. All isolates were non-susceptible to three or more different antimicrobial classes, and 75 % of non-serotype 2 isolates were non-susceptible to six or more classes (compared to 37.5 % of serotype 2 isolates). AMR genes were found on integrative and conjugative elements previously observed in other species, suggesting a mobile gene pool that can be accessed by invasive disease isolates. This article contains data hosted by Microreact.

Campylobacter enterocolitis may lead to post-infection irritable bowel syndrome (PI-IBS) and while some C. jejuni strains are more likely than others to cause human disease, genomic and virulence characteristics promoting PI-IBS development remain uncharacterized. We combined pangenome-wide association studies and phenotypic assays to compare C. jejuni isolates from patients who developed PI-IBS with those who did not. We show that variation in bacterial stress response (Cj0145_phoX), adhesion protein (Cj0628_CapA), and core biosynthetic pathway genes (biotin: Cj0308_bioD; purine: Cj0514_purQ; isoprenoid: Cj0894c_ispH) were associated with PI-IBS development. In vitro assays demonstrated greater adhesion, invasion, IL-8 and TNFα secretion on colonocytes with PI-IBS compared to PI-no-IBS strains. A risk-score for PI-IBS development was generated using 22 genomic markers, four of which were from Cj1631c, a putative heme oxidase gene linked to virulence. Our finding that specific Campylobacter genotypes confer greater in vitro virulence and increased risk of PI-IBS has potential to improve understanding of the complex host-pathogen interactions underlying this condition.

Although molecular genetic approaches have greatly increased our understanding of the evolution and spread of antibiotic resistance genes, there are fewer studies on the dynamics of antibiotic - bacterial (A-B) interactions, especially with respect to stereochemistry. Addressing this knowledge gap requires an interdisciplinary synthesis, and the development of sensitive and selective analytical tools. Here we describe SAM (stereoselective antimicrobial metabolism) workflow, a novel interdisciplinary approach for assessing bacterial resistance mechanisms in the context of A-B interactions that utilise a combination of whole genome sequencing and mass spectrometry. Chloramphenicol was used to provide proof-of-concept to demonstrate the importance of stereoselective metabolism by resistant environmental bacteria. Our data shows that chloramphenicol can be stereoselectively transformed via microbial metabolism with R,R-(-)-CAP being subject to extensive metabolic transformation by an environmental bacterial strain. In contrast S,S-(+)-CAP is not metabolised by this bacterial strain, possibly due to the lack of previous exposure to this isomer in the absence of historical selective pressure to evolve metabolic capacity.

The use of whole-genome phylogenetic analysis has revolutionized our understanding of the evolution and spread of many important bacterial pathogens due to the high resolution view it provides. However, the majority of such analyses do not consider the potential role of accessory genes when inferring evolutionary trajectories. Moreover, the recently discovered importance of the switching of gene regulatory elements suggests that an exhaustive analysis, combining information from core and accessory genes with regulatory elements could provide unparalleled detail of the evolution of a bacterial population. Here we demonstrate this principle by applying it to a worldwide multi-host sample of the important pathogenic E. coli lineage ST131. Our approach reveals the existence of multiple circulating subtypes of the major drug-resistant clade of ST131 and provides the first ever population level evidence of core genome substitutions in gene regulatory regions associated with the acquisition and maintenance of different accessory genome elements.

Campylobacter concisus is an emerging pathogen associated with inflammatory bowel disease (IBD), yet little is known about the genetic diversity of C. concisus in relation to host niches and disease. We isolated 104 C. concisus isolates from saliva, mucosal biopsies and faecal samples from 41 individuals (26 IBD, 3 Gastroenteritis (GE), 12 Healthy controls (HC)). Whole genomes were sequenced and the dataset pan-genome examined, and genomic information was used for typing using multi-locus-sequence typing (MLST). C. concisus isolates clustered into two main groups/genomospecies (GS) with 71 distinct sequence types (STs) represented. Sampling site (p < 0.001), rather than disease phenotype (p = 1.00) was associated with particular GS. We identified 97 candidate genes associated with increase or decrease in prevalence during the anatomical descent from the oral cavity to mucosal biopsies to faeces. Genes related to cell wall/membrane biogenesis were more common in oral isolates, whereas genes involved in cell transport, metabolism and secretory pathways were more prevalent in enteric isolates. Furthermore, there was no correlation between individual genetic diversity and clinical phenotype. This study confirms the genetic heterogeneity of C. concisus and provides evidence that genomic variation is related to the source of isolation, but not clinical phenotype.