Genome-wide association of functional traits linked with Campylobacter jejuni survival from farm to fork

Genome-Wide association study to identify genetic markers associated with Campylobacter jejuni motility

The ability of Campylobacter jejuni to survive and persist under harsh conditions is linked to the presence of flagella. This structure promotes the motility of the bacteria towards their optimum environment. The aim of this study was to examine the genetic basis for motility within 136 C. jejuni isolates through two different Genome-Wide Association Studies, gene presence/absence and Single Nucleotide Polymorphisms (SNPs). The motility phenotype was widely distributed across the phylogeny with large intra-lineage swarming performance variabilities. Accessory genes significantly associated with motility were found in four key genomic regions. One of these regions affected the Cj0727-Cj0733 operon, that encodes a putative ABC transporter system for phosphate uptake, while other influenced the capsule biosynthesis locus. Multiple SNPs mostly linked to increased motility were also discovered in clusters of genes, with special relevance to transport and membrane proteins. Therefore, the capsule and membrane composition might influence nutrient transfer, further impacting the protonmotive force that drives flagellar motor rotation in C. jejuni. The study provides novel genetic markers with a potential role in the motility phenotype of the pathogen.

aurora: a machine learning gwas tool for analyzing microbial habitat adaptation

A primary goal of microbial genome-wide association studies is identifying genomic variants associated with a particular habitat. Existing tools fail to identify known causal variants if the analyzed trait shaped the phylogeny. Furthermore, due to inclusion of allochthonous strains or metadata errors, the stated sources of strains in public databases are often incorrect, and strains may not be adapted to the habitat from which they were isolated. We describe a new tool, aurora, that identifies autochthonous strains and the genes associated with habitats while acknowledging the potential role of the habitat adaptation trait in shaping phylogeny.

Characterization of a ST137 multidrug-resistant Campylobacter jejuni strain with a tet(O)-positive genomic island from a bloodstream infection patient

Campylobacter jejuni (C. jejuni) is a major cause of gastroenteritis and rarely cause bloodstream infection. Herein, we characterized a multidrug-resistant C. jejuni strain LZCJ isolated from a tumor patient with bloodstream infection. LZCJ was resistant to norfloxacin, ampicillin, ceftriaxone, ciprofloxacin and tetracycline. It showed high survival rate in serum and acidic environment. Whole genome sequencing (WGS) analysis revealed that strain LZCJ had a single chromosome of 1,629,078 bp (30.6 % G + C content) and belonged to the ST137 lineage. LZCJ shared the highest identity of 99.66 % with the chicken-derived C. jejuni MTVDSCj20. Four antimicrobial resistance genes (ARGs) were detected, blaOXA-61, tet(O), gyrA (T86I), and cmeR (G144D and S207G). In addition, a 12,746 bp genomic island GI_LZCJ carrying 15 open reading frames (ORFs) including the resistance gene tet(O) was identified. Sequence analysis found that the GI_LZCJ was highly similar to the duck-derived C. jejuni ZS004, but with an additional ISChh1-like sequence. 137 non-synonymous mutations in motility related genes (flgF, fapR, flgS), capsular polysaccharide (CPS) coding genes (kpsE, kpsF, kpsM, kpsT), metabolism associated genes (nuoF, nuoG, epsJ, holB), and transporter related genes (comEA, gene0911) were confirmed in LZCJ compared with the best closed chicken-derived strain MTVDSCj20. Our study showed that C. jejuni strain LZCJ was highly similar to the chicken-derived strain MTVDSCj20 but with a lot of SNPs involved in motility, CPS and metabolism coding genes. This strain possessed a tet(O)-positive genomic island GI_LZCJ, which was closed to duck-derived C. jejuni ZS004, but with an additional ISChh1-like sequence. The above data indicated that the LZCJ strain may originate from foodborne bacteria on animals and the importance of continuous surveillance for the spread of foodborne bacteria.

Methods for studying microbial acid stress responses: from molecules to populations

The study of how micro-organisms detect and respond to different stresses has a long history of producing fundamental biological insights while being simultaneously of significance in many applied microbiological fields including infection, food and drink manufacture, and industrial and environmental biotechnology. This is well-illustrated by the large body of work on acid stress. Numerous different methods have been used to understand the impacts of low pH on growth and survival of micro-organisms, ranging from studies of single cells to large and heterogeneous populations, from the molecular or biophysical to the computational, and from well-understood model organisms to poorly defined and complex microbial consortia. Much is to be gained from an increased general awareness of these methods, and so the present review looks at examples of the different methods that have been used to study acid resistance, acid tolerance, and acid stress responses, and the insights they can lead to, as well as some of the problems involved in using them. We hope this will be of interest both within and well beyond the acid stress research community.

Epistasis, core-genome disharmony, and adaptation in recombining bacteria

Recombination of short DNA fragments via horizontal gene transfer (HGT) can introduce beneficial alleles, create genomic disharmony through negative epistasis, and create adaptive gene combinations through positive epistasis. For non-core (accessory) genes, the negative epistatic cost is likely to be minimal because the incoming genes have not co-evolved with the recipient genome and are frequently observed as tightly linked cassettes with major effects. By contrast, interspecific recombination in the core genome is expected to be rare because disruptive allelic replacement is likely to introduce negative epistasis. Why then is homologous recombination common in the core of bacterial genomes? To understand this enigma, we take advantage of an exceptional model system, the common enteric pathogens Campylobacter jejuni and C. coli that are known for very high magnitude interspecies gene flow in the core genome. As expected, HGT does indeed disrupt co-adapted allele pairings, indirect evidence of negative epistasis. However, multiple HGT events enable recovery of the genome's co-adaption between introgressing alleles, even in core metabolism genes (e.g., formate dehydrogenase). These findings demonstrate that, even for complex traits, genetic coalitions can be decoupled, transferred, and independently reinstated in a new genetic background-facilitating transition between fitness peaks. In this example, the two-step recombinational process is associated with C. coli that are adapted to the agricultural niche.IMPORTANCEGenetic exchange among bacteria shapes the microbial world. From the acquisition of antimicrobial resistance genes to fundamental questions about the nature of bacterial species, this powerful evolutionary force has preoccupied scientists for decades. However, the mixing of genes between species rests on a paradox: 0n one hand, promoting adaptation by conferring novel functionality; on the other, potentially introducing disharmonious gene combinations (negative epistasis) that will be selected against. Taking an interdisciplinary approach to analyze natural populations of the enteric bacteria Campylobacter, an ideal example of long-range admixture, we demonstrate that genes can independently transfer across species boundaries and rejoin in functional networks in a recipient genome. The positive impact of two-gene interactions appears to be adaptive by expanding metabolic capacity and facilitating niche shifts through interspecific hybridization. This challenges conventional ideas and highlights the possibility of multiple-step evolution of multi-gene traits by interspecific introgression.

Genomic tailoring of autogenous poultry vaccines to reduce Campylobacter from farm to fork

Campylobacter is a leading cause of food-borne gastroenteritis worldwide, linked to the consumption of contaminated poultry meat. Targeting this pathogen at source, vaccines for poultry can provide short-term caecal reductions in Campylobacter numbers in the chicken intestine. However, this approach is unlikely to reduce Campylobacter in the food chain or human incidence. This is likely as vaccines typically target only a subset of the high genomic strain diversity circulating among chicken flocks, and rapid evolution diminishes vaccine efficacy over time. To address this, we used a genomic approach to develop a whole-cell autogenous vaccine targeting isolates harbouring genes linked to survival outside of the host. We hyper-immunised a whole major UK breeder farm to passively target offspring colonisation using maternally-derived antibody. Monitoring progeny, broiler flocks revealed a near-complete shift in the post-vaccination Campylobacter population with an ~50% reduction in isolates harbouring extra-intestinal survival genes and a significant reduction of Campylobacter cells surviving on the surface of meat. Based on these findings, we developed a logistic regression model that predicted that vaccine efficacy could be extended to target 65% of a population of clinically relevant strains. Immuno-manipulation of poultry microbiomes towards less harmful commensal isolates by competitive exclusion, has major potential for reducing pathogens in the food production chain.

Evidence of potential Campylobacter jejuni zooanthroponosis in captive macaque populations

Non-human primates share recent common ancestry with humans and exhibit comparable disease symptoms. Here, we explored the transmission potential of enteric bacterial pathogens in monkeys exhibiting symptoms of recurrent diarrhoea in a biomedical research facility in China. The common zoonotic bacterium Campylobacter jejuni was isolated from macaques (Macaca mulatta and Macaca fascicularis) and compared to isolates from humans and agricultural animals in Asia. Among the monkeys sampled, 5 % (44/973) tested positive for C. jejuni, 11 % (5/44) of which displayed diarrhoeal symptoms. Genomic analysis of monkey isolates, and 1254 genomes from various sources in Asia, were used to identify the most likely source of human infection. Monkey and human isolates shared high average nucleotide identity, common MLST clonal complexes and clustered together on a phylogeny. Furthermore, the profiles of putative antimicrobial resistance genes were similar between monkeys and humans. Taken together these findings suggest that housed macaques became infected with C. jejuni either directly from humans or via a common contamination source.

Genomic and Antimicrobial Surveillance of Campylobacter Population in Italian Poultry

Campylobacter is one of the most common foodborne diseases worldwide with increasing rates of antibiotic resistance. Most cases of campylobacteriosis can be traced back to the consumption of poultry meat. Despite many efforts to reduce contamination in farms and in slaughterhouses, the persistence of this pathogen in poultry products remains a problem. This study aimed to evaluate the genetic diversity and antibiotic resistance of 542 C. jejuni and C. coli in Italian poultry, in the framework of two National Monitoring Programs. Genomes were screened for antibiotic resistance, virulence determinants and contextualized within a global collection of C. jejuni. ST2116, ST2863 and ST 832 were the most prevalent and significantly associated with Italian poultry. A worrying increase in resistance to quinolones, fluoroquinolones and tetracycline was observed in C. jejuni, while an increased occurrence of multidrug resistant (MDR) strains and strains resistant to macrolides was detected in C. coli. Low resistance rates were found for aminoglycosides. Molecular resistance determinants were consistent with the phenotypic resistance for tetracycline and quinolones. In silico analysis revealed 119 genes associated with virulence factors, with a notably higher prevalence of some genes in ST2863 genomes. This study highlights the increased resistance to macrolides and the emergence of MDR strains for C. coli, the genetic basis of AMR and the predominance of two genotypes among Campylobacter strains isolated from the Italian poultry farms.

Roles of viable but non-culturable state in the survival of Campylobacter jejuni

Despite being considered fragile and fastidious, Campylobacter jejuni is the most prevalent cause of foodborne bacterial gastroenteritis, and chicken meat is considered the main vehicle of transmission to humans. This agent can survive adverse conditions in the form of biofilms, but extreme stress (nutritional, oxidative and thermal) promotes the acquisition of a state called viable but not culturable (VBNC). The emergence of this pathogen worldwide and the recent international requirements in its control instigated us to qualitatively and quantitatively estimate the time required for the acquisition of the VBNC form in 27 strains of C. jejuni, characterize morphological aspects, determine its adaptive and invasive potential and perform comparative metabolomic evaluation. Extreme stress promoted the complete acquisition of the VBNC form in a mean time of 26 days. Starting from an average initial count of 7.8 log CFU/mL, the first four days determined the greatest average reduction of the culturable form of 3.2 log CFU/mL. The scanning and transmission image analyses showed a transition from the typical viable form (VT) to the VBNC form, with initial acquisition of the straight rod shape, followed by loss of the flagella and subdivision into two to 11 imperfect cocci arranged in a chain and rich in cellular content, until their individual release. RT-PCR identified the presence of ciaB and p19 transcripts in the 27 cultivable C. jejuni strains, a character maintained in the VBNC form only for p19 and in 59.3% (16/27) of the VBNC strains for the ciaB gene. The average inoculation of 1.8 log CFU/mL of C. jejuni VBNC into primary chicken embryo hepatocyte cells promoted the occurrence of apoptosis processes significantly after 24 hours of contact by one of the strains tested. In C. jejuni VBNC, we detected higher expression of metabolites linked to protective and adaptation mechanisms and of volatile organic precursor compounds indicative of metabolism interruption. The oscillations in the time of acquisition of the VBNC form together with the presence of transcripts for ciaB and p19, the identification of cell lysis and metabolites that ensure the maintenance of the pathogen alert to the fact that C. jejuni VBNC remains virulent and adapted to stress, which makes evident the potential danger of this latent form, which is not detectable by official methodologies.

The Use of Interdisciplinary Approaches to Understand the Biology of Campylobacter jejuni

Campylobacter jejuni is a bacterial pathogen recognised as a major cause of foodborne illness worldwide. While Campylobacter jejuni generally does not grow outside its host, it can survive outside of the host long enough to pose a health concern. This review presents an up-to-date description and evaluation of biological, mathematical, and statistical approaches used to understand the behaviour of this foodborne pathogen and suggests future avenues which can be explored. Specifically, the incorporation of mathematical modelling may aid the understanding of C. jejuni biofilm formation both outside and inside the host. Predictive studies may be improved by the introduction of more standardised protocols for assessments of disinfection methods and by assessment of novel physical disinfection strategies as well as assessment of the efficiency of plant extracts on C. jejuni eradication. A full description of the metabolic pathways of C. jejuni, which is needed for the successful application of metabolic models, is yet to be achieved. Finally, a shift from animal models (except for those that are a source of human campylobacteriosis) to human-specific data may be made possible due to recent technological advancements, and this may lead to more accurate predictions of human infections.

Genomic insights into zoonotic transmission and antimicrobial resistance in Campylobacter jejuni from farm to fork: a one health perspective

BACKGROUND

Campylobacteriosis represents a global public health threat with various socio-economic impacts. Among different Campylobacter species, Campylobacter jejuni (C. jejuni) is considered to be the foremost Campylobacter species responsible for most of gastrointestinal-related infections. Although these species are reported to primarily inhabit birds, its high genetic and phenotypic diversity allowed their adaptation to other animal reservoirs and to the environment that may impact on human infection.

MAIN BODY

A stringent and consistent surveillance program based on high resolution subtyping is crucial. Recently, different epidemiological investigations have implemented high-throughput sequencing technologies and analytical pipelines for higher resolution subtyping, accurate source attribution, and detection of antimicrobial resistance determinants among these species. In this review, we aim to present a comprehensive overview on the epidemiology, clinical presentation, antibiotic resistance, and transmission dynamics of Campylobacter, with specific focus on C. jejuni. This review also summarizes recent attempts of applying whole-genome sequencing (WGS) coupled with bioinformatic algorithms to identify and provide deeper insights into evolutionary and epidemiological dynamics of C. jejuni precisely along the farm-to-fork continuum.

CONCLUSION

WGS is a valuable addition to traditional surveillance methods for Campylobacter. It enables accurate typing of this pathogen and allows tracking of its transmission sources. It is also advantageous for in silico characterization of antibiotic resistance and virulence determinants, and hence implementation of control measures for containment of infection.

Metaphenotypes associated with recurrent genomic lineages of Campylobacter jejuni responsible for human infections in Luxembourg

Campylobacter jejuni is a leading cause of foodborne illnesses worldwide. Although considered fragile, this microaerophilic bacterium is able to survive in various challenging environments, which subsequently constitutes multiple sources of transmission for human infection. To test the assumption of acquiring specific features for adaptation and survival, we established a workflow of phenotypic tests related to the survival and the persistence of recurrent and sporadic strains. A representative collection of 83 strains isolated over 13 years from human, mammal, poultry, and environmental sources in Luxembourg, representing different spreading patterns (endemic, epidemic, and sporadic), was screened for survival to oxidative stresses, for acclimating to aerobic conditions (AC), and for persistence on abiotic surfaces. Using the cgMLST Oxford typing scheme for WGS data, the collection was classified into genomic lineages corresponding to host-generalist strains (lineages A and D, CC ST-21), host-specific strains (lineage B, CC ST-257 and lineage C, CC ST-464) and sporadic strains. We established that when a strain survives concentrations beyond 0.25 mM superoxide stress, it is six times more likely to survive hyperoxide stress and that a highly adherent strain is 14 times more likely to develop a biofilm. Surprisingly, more than half of the strains could acclimate to AC but this capacity does not explain the difference between recurrent genomic lineages and sporadic strains and the survival to oxidative stresses, while recurrent strains have a significantly higher adhesion/biofilm formation capacity than sporadic ones. From this work, the genomic lineages with more stable genomes could be characterized by a specific combination of phenotypes, called metaphenotypes. From the functional genomic analyses, the presence of a potentially functional T6SS in the strains of lineage D might explain the propensity of these strains to be strong biofilm producers. Our findings support the hypothesis that phenotypical abilities contribute to the spatio-temporal adaptation and survival of stable genomic lineages. It suggests a selection of better-adapted and persistent strains in challenging stress environments, which could explain the prevalence of these lineages in human infections.

Emerging patterns of fluoroquinolone resistance in Campylobacter jejuni in the UK [1998-2018]

Campylobacter jejuni (C.jejuni) is the most common causative agent of bacterial food poisoning worldwide and is known to be genetically highly diverse. C. jejuni is increasingly resistant to fluoroquinolone antibiotics, but very few studies have investigated variant-specific patterns of resistance across time. Here we use statistical modelling and clustering techniques to investigate patterns of fluoroquinolone resistance amongst 10,359 UK isolates from human disease sampled over 20 years. We observed six distinct patterns of fluoroquinolone sensitivity/resistance in C. jejuni across time, grouping by clonal complex (CC). Some CCs were fully resistant, some shifted from susceptible to resistant following a sigmoidal shape, and some remained susceptible over time. Our findings indicate that the fluoroquinolone resistance patterns of C. jejuni are complicated and cannot be analysed as a single species but divided into variant dynamics so that the factors driving resistance can be thoroughly investigated.

Identification of Novel Phage Resistance Mechanisms in Campylobacter jejuni by Comparative Genomics

Phages infecting Campylobacter jejuni are considered a promising intervention strategy at broiler farms, yet phage sensitivity of naturally occurring poultry isolates is not well studied. Here, we investigated phage sensitivity and identified resistance mechanisms of C. jejuni strains originating from Danish broilers belonging to the most prevalent MLST (ST) types. Determining plaque formation of 51 phages belonging to Fletchervirus or Firehammervirus showed that 21 out of 31 C. jejuni strains were susceptible to at least one phage. While C. jejuni ST-21 strains encoded the common phase variable O-methyl phosphoramidate (MeOPN) receptor of the Fletchervirus and were only infected by these phages, ST-45 strains did not encode this receptor and were exclusively infected by Firehammervirus phages. To identify internal phage resistance mechanism in ST-21 strains, we performed comparative genomics of two strains, CAMSA2002 sensitive to almost all Fletchervirus phages and CAMSA2038, resistant to all 51 phages. The strains encoded diverse clustered regularly interspaced short palindromic repeats (CRISPR) spacers but none matched the tested phages. Sequence divergence was also observed in a predicted SspE homolog and putative restriction modification systems including a methyl-specific McrBC endonuclease. Furthermore, when mcrB was deleted, CAMSA2038 became sensitive to 17 out of 43 phages, three being Firehammervirus phages that otherwise did not infect any ST-21 strains. Yet, 16 phages demonstrated significantly lower efficiencies of plating on the mcrB mutant suggesting additional resistance mechanism still restricting phage propagation in CAMSA2038. Thus, our work demonstrates that C. jejuni isolates originating from broilers may have acquired several resistance mechanisms to successfully prevent phage infection in their natural habitat.

Genome-wide association study of gastric cancer- and duodenal ulcer-derived Helicobacter pylori strains reveals discriminatory genetic variations and novel oncoprotein candidates

Genome-wide association studies (GWASs) can reveal genetic variations associated with a phenotype in the absence of any hypothesis of candidate genes. The problem of false-positive sites linked with the responsible site might be bypassed in bacteria with a high homologous recombination rate, such as Helicobacter pylori, which causes gastric cancer. We conducted a small-sample GWAS (125 gastric cancer cases and 115 controls) followed by prediction of gastric cancer and control (duodenal ulcer) H. pylori strains. We identified 11 single nucleotide polymorphisms (eight amino acid changes) and three DNA motifs that, combined, allowed effective disease discrimination. They were often informative of the underlying molecular mechanisms, such as electric charge alteration at the ligand-binding pocket, alteration in subunit interaction, and mode-switching of DNA methylation. We also identified three novel virulence factors/oncoprotein candidates. These results provide both defined targets for further informatic and experimental analyses to gain insights into gastric cancer pathogenesis and a basis for identifying a set of biomarkers for distinguishing these H. pylori-related diseases.

A Campylobacter integrative and conjugative element with a CRISPR-Cas9 system targeting competing plasmids: a history of plasmid warfare?

Microbial genomes are highly adaptable, with mobile genetic elements (MGEs) such as integrative conjugative elements (ICEs) mediating the dissemination of new genetic information throughout bacterial populations. This is countered by defence mechanisms such as CRISPR-Cas systems, which limit invading MGEs by sequence-specific targeting. Here we report the distribution of the pVir, pTet and PCC42 plasmids and a new 70–129 kb ICE (CampyICE1) in the foodborne bacterial pathogens Campylobacter jejuni and Campylobacter coli. CampyICE1 contains a degenerated Type II-C CRISPR system consisting of a sole Cas9 protein, which is distinct from the previously described Cas9 proteins from C. jejuni and C. coli. CampyICE1 is conserved in structure and gene order, containing blocks of genes predicted to be involved in recombination, regulation and conjugation. CampyICE1 was detected in 134/5829 (2.3 %) C. jejuni genomes and 92/1347 (6.8 %) C. coli genomes. Similar ICEs were detected in a number of non-jejuni/coli Campylobacter species, although these lacked a CRISPR-Cas system. CampyICE1 carries three separate short CRISPR spacer arrays containing a combination of 108 unique spacers and 16 spacer-variant families. A total of 69 spacers and 10 spacer-variant families (63.7 %) were predicted to target Campylobacter plasmids. The presence of a functional CampyICE1 Cas9 protein and matching anti-plasmid spacers was associated with the absence of the pVir, pTet and pCC42 plasmids (188/214 genomes, 87.9 %), suggesting that the CampyICE1-encoded CRISPR-Cas has contributed to the exclusion of competing plasmids. In conclusion, the characteristics of the CRISPR-Cas9 system on CampyICE1 suggests a history of plasmid warfare in Campylobacter.

Machine learning to predict the source of campylobacteriosis using whole genome data

Campylobacteriosis is among the world's most common foodborne illnesses, caused predominantly by the bacterium Campylobacter jejuni. Effective interventions require determination of the infection source which is challenging as transmission occurs via multiple sources such as contaminated meat, poultry, and drinking water. Strain variation has allowed source tracking based upon allelic variation in multi-locus sequence typing (MLST) genes allowing isolates from infected individuals to be attributed to specific animal or environmental reservoirs. However, the accuracy of probabilistic attribution models has been limited by the ability to differentiate isolates based upon just 7 MLST genes. Here, we broaden the input data spectrum to include core genome MLST (cgMLST) and whole genome sequences (WGS), and implement multiple machine learning algorithms, allowing more accurate source attribution. We increase attribution accuracy from 64% using the standard iSource population genetic approach to 71% for MLST, 85% for cgMLST and 78% for kmerized WGS data using the classifier we named aiSource. To gain insight beyond the source model prediction, we use Bayesian inference to analyse the relative affinity of C. jejuni strains to infect humans and identified potential differences, in source-human transmission ability among clonally related isolates in the most common disease causing lineage (ST-21 clonal complex). Providing generalizable computationally efficient methods, based upon machine learning and population genetics, we provide a scalable approach to global disease surveillance that can continuously incorporate novel samples for source attribution and identify fine-scale variation in transmission potential.

High-resolution sweep metagenomics using fast probabilistic inference

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.

Heuristic and Hierarchical-Based Population Mining of Salmonella enterica Lineage I Pan-Genomes as a Platform to Enhance Food Safety

The recent incorporation of bacterial whole-genome sequencing (WGS) into Public Health laboratories has enhanced foodborne outbreak detection and source attribution. As a result, large volumes of publicly available datasets can be used to study the biology of foodborne pathogen populations at an unprecedented scale. To demonstrate the application of a heuristic and agnostic hierarchical population structure guided pan-genome enrichment analysis (PANGEA), we used populations of S. enterica lineage I to achieve two main objectives: (i) show how hierarchical population inquiry at different scales of resolution can enhance ecological and epidemiological inquiries; and (ii) identify population-specific inferable traits that could provide selective advantages in food production environments. Publicly available WGS data were obtained from NCBI database for three serovars of Salmonella enterica subsp. enterica lineage I (S. Typhimurium, S. Newport, and S. Infantis). Using the hierarchical genotypic classifications (Serovar, BAPS1, ST, cgMLST), datasets from each of the three serovars showed varying degrees of clonal structuring. When the accessory genome (PANGEA) was mapped onto these hierarchical structures, accessory loci could be linked with specific genotypes. A large heavy-metal resistance mobile element was found in the Monophasic ST34 lineage of S. Typhimurium, and laboratory testing showed that Monophasic isolates have on average a higher degree of copper resistance than the Biphasic ones. In S. Newport, an extra sugE gene copy was found among most isolates of the ST45 lineage, and laboratory testing of multiple isolates confirmed that isolates of S. Newport ST45 were on average less sensitive to the disinfectant cetylpyridimium chloride than non-ST45 isolates. Lastly, data-mining of the accessory genomic content of S. Infantis revealed two cryptic Ecotypes with distinct accessory genomic content and distinct ecological patterns. Poultry appears to be the major reservoir for Ecotype 1, and temporal analysis further suggested a recent ecological succession, with Ecotype 2 apparently being displaced by Ecotype 1. Altogether, the use of a heuristic hierarchical-based population structure analysis that includes bacterial pan-genomes (core and accessory genomes) can (1) improve genomic resolution for mapping populations and accessing epidemiological patterns; and (2) define lineage-specific informative loci that may be associated with survival in the food chain.

Whole Resistome Analysis in Campylobacter jejuni and C. coli Genomes Available in Public Repositories

Campylobacter spp. are the most frequent agent of human gastroenteritis worldwide, and the spread of multidrug-resistant strains makes the clinical treatment difficult. The current study presents the resistome analysis of 39,798 Campylobacter jejuni and 11,920 Campylobacter coli genomes available in public repositories. Determinants of resistance to β-lactams (Be) and tetracyclines (Te) were the most frequent for both species, with resistance to quinolones (Qu) as the third most important on C. jejuni and to aminoglycosides (Am) on C. coli. Moreover, resistance to Te, Qu, and Am was frequently found in co-occurrence with resistance to other antibiotic families. Geographical differences on clonal complexes distribution were found for C. jejuni and on resistome genotypes for both C. jejuni and C. coli species. Attending to the resistome patterns by isolation source, three main clusters of genomes were found on C. jejuni genomes at antimicrobial resistance gene level. The first cluster was formed by genomes from human, food production animals (e.g., sheep, cow, and chicken), and food (e.g., dairy products) isolates. The higher incidence of tet(O), associated with tetracycline resistance, and the gyrA (T86I) single-nucleotide polymorphism (SNP), associated with quinolone resistance, among genomes from this cluster could be due to the intense use of these antibiotics in veterinary and human clinical settings. Similarly, a high incidence of tet(O) genes of C. coli genomes from pig, cow, and turkey was found. Moreover, the cluster based on resistome patterns formed by C. jejuni and C. coli genomes of human, turkey, and chicken origin is in agreement with previous observations reporting chicken or poultry-related environments as the main source of human campylobacteriosis infections. Most clonal complexes (CCs) associated with chicken host specialization (e.g., ST-354, ST-573, ST-464, and ST-446) were the CCs with the highest prevalence of determinants of resistance to Be, Qu, and Te. Finally, a clear trend toward an increase in the occurrence of Te and Qu resistance determinants on C. jejuni, linked to the spread of the co-occurrence of the bla_OXA–61 and tet(O)-tet(O/W/O) genes and the gyrA (T86I) SNP, was found from 2001 to date in Europe.

ProkEvo: an automated, reproducible, and scalable framework for high-throughput bacterial population genomics analyses

Whole Genome Sequence (WGS) data from bacterial species is used for a variety of applications ranging from basic microbiological research, diagnostics, and epidemiological surveillance. The availability of WGS data from hundreds of thousands of individual isolates of individual microbial species poses a tremendous opportunity for discovery and hypothesis-generating research into ecology and evolution of these microorganisms. Flexibility, scalability, and user-friendliness of existing pipelines for population-scale inquiry, however, limit applications of systematic, population-scale approaches. Here, we present ProkEvo, an automated, scalable, reproducible, and open-source framework for bacterial population genomics analyses using WGS data. ProkEvo was specifically developed to achieve the following goals: (1) Automation and scaling of complex combinations of computational analyses for many thousands of bacterial genomes from inputs of raw Illumina paired-end sequence reads; (2) Use of workflow management systems (WMS) such as Pegasus WMS to ensure reproducibility, scalability, modularity, fault-tolerance, and robust file management throughout the process; (3) Use of high-performance and high-throughput computational platforms; (4) Generation of hierarchical-based population structure analysis based on combinations of multi-locus and Bayesian statistical approaches for classification for ecological and epidemiological inquiries; (5) Association of antimicrobial resistance (AMR) genes, putative virulence factors, and plasmids from curated databases with the hierarchically-related genotypic classifications; and (6) Production of pan-genome annotations and data compilation that can be utilized for downstream analysis such as identification of population-specific genomic signatures. The scalability of ProkEvo was measured with two datasets comprising significantly different numbers of input genomes (one with ~2,400 genomes, and the second with ~23,000 genomes). Depending on the dataset and the computational platform used, the running time of ProkEvo varied from ~3-26 days. ProkEvo can be used with virtually any bacterial species, and the Pegasus WMS uniquely facilitates addition or removal of programs from the workflow or modification of options within them. To demonstrate versatility of the ProkEvo platform, we performed a hierarchical-based population structure analyses from available genomes of three distinct pathogenic bacterial species as individual case studies. The specific case studies illustrate how hierarchical analyses of population structures, genotype frequencies, and distribution of specific gene functions can be integrated into an analysis. Collectively, our study shows that ProkEvo presents a practical viable option for scalable, automated analyses of bacterial populations with direct applications for basic microbiology research, clinical microbiological diagnostics, and epidemiological surveillance.

Genome-wide insights into population structure and host specificity of Campylobacter jejuni

The zoonotic pathogen Campylobacter jejuni is among the leading causes of foodborne diseases worldwide. While C. jejuni colonises many wild animals and livestock, persistence mechanisms enabling the bacterium to adapt to host species' guts are not fully understood. In order to identify putative determinants influencing host preferences of distinct lineages, bootstrapping based on stratified random sampling combined with a k-mer-based genome-wide association was conducted on 490 genomes from diverse origins in Germany and Canada. We show a strong association of both the core and the accessory genome characteristics with distinct host animal species, indicating multiple adaptive trajectories defining the evolution of C. jejuni lifestyle preferences in different ecosystems. Here, we demonstrate that adaptation towards a specific host niche ecology is most likely a long evolutionary and multifactorial process, expressed by gene absence or presence and allele variations of core genes. Several host-specific allelic variants from different phylogenetic backgrounds, including dnaE, rpoB, ftsX or pycB play important roles for genome maintenance and metabolic pathways. Thus, variants of genes important for C. jejuni to cope with specific ecological niches or hosts may be useful markers for both surveillance and future pathogen intervention strategies.

Genetic variation and function: revealing potential factors associated with microbial phenotypes

Innovations in sequencing technology have generated voluminous microbial and host genomic data, making it possible to detect these genetic variations and analyze the function influenced by them. Recently, many studies have linked such genetic variations to phenotypes through association or comparative analysis, which have further advanced our understanding of multiple microbial functions. In this review, we summarized the application of association analysis in microbes like Mycobacterium tuberculosis, focusing on screening of microbial genetic variants potentially associated with phenotypes such as drug resistance, pathogenesis and novel drug targets etc.; reviewed the application of additional comparative genomic or transcriptomic methods to identify genetic factors associated with functions in microbes; expanded the scope of our study to focus on host genetic factors associated with certain microbes or microbiome and summarized the recent host genetic variations associated with microbial phenotypes, including susceptibility and load after infection of HIV, presence/absence of different taxa, and quantitative traits of microbiome, and lastly, discussed the challenges that may be encountered and the apparent or potential viable solutions. Gene-function analysis of microbe and microbiome is still in its infancy, and in order to unleash its full potential, it is necessary to understand its history, current status, and the challenges hindering its development.

Next-generation microbiology: from comparative genomics to gene function

Microbiology is at a turning point in its 120-year history. Widespread next-generation sequencing has revealed genetic complexity among bacteria that could hardly have been imagined by pioneers such as Pasteur, Escherich and Koch. This data cascade brings enormous potential to improve our understanding of individual bacterial cells and the genetic basis of phenotype variation. However, this revolution in data science cannot replace established microbiology practices, presenting the challenge of how to integrate these new techniques. Contrasting comparative and functional genomic approaches, we evoke molecular microbiology theory and established practice to present a conceptual framework and practical roadmap for next-generation microbiology.

The flavodoxin FldA activates the class Ia ribonucleotide reductase of Campylobacter jejuni

Campylobacter jejuni is a microaerophilic zoonotic pathogen with an atypical respiratory Complex I that oxidizes a flavodoxin (FldA) instead of NADH. FldA is essential for viability and is reduced via pyruvate and 2-oxoglutarate oxidoreductases (POR/OOR). Here, we show that FldA can also be reduced by FqrB (Cj0559), an NADPH:FldA reductase. An fqrB deletion mutant was viable but displayed a significant growth defect. FqrB is related to flavoprotein reductases from Gram-positive bacteria that can reduce NrdI, a specialized flavodoxin that is needed for tyrosyl radical formation in NrdF, the beta subunit of class 1b-type (Mn) ribonucleotide reductase (RNR). However, C. jejuni possesses a single class Ia-type (Fe) RNR (NrdAB) that would be expected to be ferredoxin dependent. We show that CjFldA is an unusually high potential flavodoxin unrelated to NrdI, yet growth of the fqrB mutant, but not the wild-type or a complemented strain, was stimulated by low deoxyribonucleoside (dRNS) concentrations, suggesting FldA links FqrB and RNR activity. Using purified proteins, we confirmed the NrdB tyrosyl radical could be regenerated in an NADPH, FqrB, and FldA dependent manner, as evidenced by both optical and electron paramagnetic resonance (EPR) spectroscopy. Thus, FldA activates RNR in C. jejuni, partly explaining its essentiality.

Genome Informatics and Machine Learning-Based Identification of Antimicrobial Resistance-Encoding Features and Virulence Attributes in Escherichia coli Genomes Representing Globally Prevalent Lineages, Including High-Risk Clonal Complexes

Escherichia coli, a ubiquitous commensal/pathogenic member from the Enterobacteriaceae family, accounts for high infection burden, morbidity, and mortality throughout the world. With emerging multidrug resistance (MDR) on a massive scale, E. coli has been listed as one of the Global Antimicrobial Resistance and Use Surveillance System (GLASS) priority pathogens. Understanding the resistance mechanisms and underlying genomic features appears to be of utmost importance to tackle further spread of these multidrug-resistant superbugs. While a few of the globally prevalent sequence types (STs) of E. coli, such as ST131, ST69, ST405, and ST648, have been previously reported to be highly virulent and harboring MDR, there is no clarity if certain ST lineages have a greater propensity to acquire MDR. In this study, large-scale comparative genomics of a total of 5,653 E. coli genomes from 19 ST lineages revealed ST-wide prevalence patterns of genomic features, such as antimicrobial resistance (AMR)-encoding genes/mutations, virulence genes, integrons, and transposons. Interpretation of the importance of these features using a Random Forest Classifier trained with 11,988 genomic features from whole-genome sequence data identified ST-specific or phylogroup-specific signature proteins mostly belonging to different protein superfamilies, including the toxin-antitoxin systems. Our study provides a comprehensive understanding of a myriad of genomic features, ST-specific proteins, and resistance mechanisms entailing different lineages of E. coli at the level of genomes; this could be of significant downstream importance in understanding the mechanisms of AMR, in clinical discovery, in epidemiology, and in devising control strategies. IMPORTANCE With the leap in whole-genome data being generated, the application of relevant methods to mine biologically significant information from microbial genomes is of utmost importance to public health genomics. Machine-learning methods have been used not only to mine, curate, or classify the data but also to identify the relevant features that could be linked to a particular class/target. This is perhaps one of the pioneering studies that has attempted to classify a large repertoire of E. coli genome data sets (5,653 genomes) belonging to 19 different STs (including well-studied as well as understudied STs) using machine learning approaches. Important features identified by these approaches have revealed ST-specific signature proteins, which could be further studied to predict possible associations with the phenotypic profiles, thereby providing a better understanding of virulence and the resistance mechanisms among different clonal lineages of E. coli.

Population Biology and Comparative Genomics of Campylobacter Species

The zoonotic pathogen Campylobacter is the leading cause for bacterial foodborne infections in humans. Campylobacters are most commonly transmitted via the consumption of undercooked poultry meat or raw milk products. The decreasing costs of whole genome sequencing enabled large genome-based analyses of the evolution and population structure of this pathogen, as well as the development of novel high-throughput molecular typing methods. Here, we review the evolutionary development and the population diversity of the two most clinically relevant Campylobacter species; C. jejuni and C. coli. The state-of-the-art phylogenetic studies showed clustering of C. jejuni lineages into host specialists and generalists with coexisting lifestyles in chicken and livestock-associated hosts, as well as the separation of C. coli isolates of riparian origin (waterfowl, water) from C. coli isolated from clinical and farm-related samples. We will give an overview of recombination between both species and the potential impact of horizontal gene transfer on host adaptation in Campylobacter. Additionally, this review briefly places the current knowledge of the population structure of other Campylobacter species such as C. lari, C. concisus and C. upsaliensis into perspective. We also provide an overview of how molecular typing methods such as multilocus sequence typing (MLST) and whole genome MLST have been used to detect and trace Campylobacter outbreaks along the food chain.

Stress resistance associated with multi-host transmission and enhanced biofilm formation at 42 °C among hyper-aerotolerant generalist Campylobacter jejuni

One of the emerging conundrums of Campylobacter food-borne illness is the bacterial ability to survive stressful environmental conditions. We evaluated the heterogeneity among 90 C. jejuni and 21 C. coli isolates from different sources in Egypt with respect to biofilm formation capabilities (under microaerobic and aerobic atmosphere) and resistance to a range of stressors encountered along the food chain (aerobic stress, refrigeration, freeze-thaw, heat, peracetic acid, and osmotic stress). High prevalence (63%) of hyper-aerotolerant (HAT) isolates was observed, exhibiting also a significantly high tolerance to heat, osmotic stress, refrigeration, and freeze-thaw stress, coupled with high biofilm formation ability which was clearly enhanced under aerobic conditions, suggesting a potential link between stress adaptation and biofilm formation. Most HAT multi-stress resistant and strong biofilm producing C. jejuni isolates belonged to host generalist clonal complexes (ST-21, ST-45, ST-48 and ST-206). These findings highlight the potential role of oxidative stress response systems in providing cross-protection (resistance to other multiple stress conditions) and enhancing biofilm formation in Campylobacter and suggest that selective pressures encountered in hostile environments have shaped the epidemiology of C. jejuni in Egypt by selecting the transmission of highly adapted isolates, thus promoting the colonization of multiple host species by important disease-causing lineages.

Genome-Wide Identification of Host-Segregating Single-Nucleotide Polymorphisms for Source Attribution of Clinical Campylobacter coli Isolates

Campylobacter is among the most common causes of gastroenteritis worldwide. Campylobacter jejuni and Campylobacter coli are the most common species causing human disease. DNA sequence-based methods for strain characterization have focused largely on C. jejuni, responsible for 80 to 90% of infections, meaning that C. coli epidemiology has lagged behind. Here, we have analyzed the genome of 450 C. coli isolates to determine genetic markers that can discriminate isolates sampled from 3 major reservoir hosts (chickens, cattle, and pigs). These markers then were applied to identify the source of infection of 147 C. coli strains from French clinical cases. Using STRUCTURE software, 259 potential host-segregating markers were revealed by probabilistic characterization of single-nucleotide polymorphism (SNP) frequency variation in strain collections from three different hosts. These SNPs were found in 41 genes or intergenic regions, mostly coding for proteins involved in motility and membrane functions. Source attribution of clinical isolates based on the differential presence of these markers confirmed chickens as the most common source of C. coli infection in France.IMPORTANCE Genome-wide and source attribution studies based on Campylobacter species have shown their importance for the understanding of foodborne infections. Although the use of multilocus sequence typing based on 7 genes from C. jejuni is a powerful method to structure populations, when applied to C. coli, results have not clearly demonstrated its robustness. Therefore, we aim to provide more accurate data based on the identification of single-nucleotide polymorphisms. Results from this study reveal an important number of host-segregating SNPs, found in proteins involved in motility, membrane functions, or DNA repair systems. These findings offer new, interesting opportunities for further study of C. coli adaptation to its environment. Additionally, the results demonstrate that poultry is potentially the main reservoir of C. coli in France.

The impact of key processing stages and flock variables on the prevalence and levels of Campylobacter on broiler carcasses

This study examined the impact of key processing stages and flock variables on the prevalence of Campylobacter on broiler carcasses. Overall, the prevalence of Campylobacter was 62% in caeca, and 68%, 65% and 62% in neck skin samples collected after evisceration, final wash and carcass chilling, respectively. Campylobacter were found in 32% of caeca, and 52%, 40% and 32% of neck skin samples collected after evisceration, final wash and carcass chilling, respectively from first thin broiler batches. Final thin broiler batches were more frequently contaminated with prevalences of 83% found in caeca, 80% in neck skin samples collected after evisceration and 83% found in neck skin samples collected after both final wash and carcass chilling stages (p < 0.05). Thinning status had a significant effect on Campylobacter counts with significantly higher counts observed in samples from final thin batches (p < 0.05). Highest Campylobacter concentrations in neck skin samples were observed at the evisceration stage in both first and final thin samples, with counts ranging from 2.0 to 3.8 log₁₀ CFU/g and 2.3 to 4.8 log₁₀ CFU/g in first and final thin batches, respectively. All first thin samples had counts below the European Union (EU) Process Hygiene Criterion threshold level of 3 log₁₀ CFU/g after chilling while 52% of final thin batches had counts above this limit.

Determining the prevalence, identity and possible origin of bacterial pathogens in soil

Soil biomes are vast, exceptionally diverse and crucial to the health of ecosystems and societies. Soils also contain an appreciable, but understudied, diversity of opportunistic human pathogens. With climate change and other forms of environmental degradation potentially increasing exposure risks to soilborne pathogens, it is necessary to gain a better understanding of their ecological drivers. Here we use the Galleria mellonella insect virulence model to selectively isolate pathogenic bacteria from soils in Cornwall (UK). We find a high prevalence of pathogenic soil bacteria with two genera, Providencia and Serratia, being especially common. Providencia alcalifaciens, P. rustigianii, Serratia liquefaciens and S. plymuthica strains were studied in more detail using phenotypic virulence and antibiotic resistance assays and whole-genome sequencing. Both genera displayed low levels of antibiotic resistance and antibiotic resistance gene carriage. However, Serratia isolates were found to carry the recently characterized metallo-β-lactamase blaSPR-1 that, although not conferring high levels of resistance in these strains, poses a potential risk of horizontal transfer to other pathogens where it could be fully functional. The Galleria assay can be a useful approach to uncover the distribution and identity of pathogenic bacteria in the environment, as well as uncover resistance genes with an environmental origin.

Whole-genome characterisation of multidrug resistant monophasic variants of Salmonella Typhimurium from pig production in Thailand

Background

Monophasic Salmonella Typhimurium or S. enterica 1,4,[5],12:i:- is among the top five serotypes reported in Thailand. In this study, nineteen monophasic S. Typhimurium from the pig production chain in Chiang Mai and Lamphun provinces during 2011–2014 were sequenced and compared to a globally disseminated clone. Isolates were probed in silico for the presence of antimicrobial resistance genes and Salmonella virulence factors, including Pathogenicity Islands.

Results

All isolates were from sequence type 34 (ST-34) and clustered similarly in core and pangenome genealogies. The two closest related isolates showed differences in only eighteen loci from whole-genome multilocus sequence typing analysis. All 19 isolates carried aminoglycoside and beta-lactam class resistance genes and genes for five or more different antibiotic classes. Seven out of 14 known SPIs were detected, including SPI-5, SPI-13 and SPI-14, which were detected in all isolates.

Conclusions

The multi-drug resistant clone, ST-34 was sampled at all stages of pork production. This clone has infiltrated global agricultural processes and poses a significant public health risk. Differences in the core and accessory genomes of the isolates we collected suggest that strains persist though the pork production process, with evidence of mutation within the core-genome and horizontal acquisition of genes, potentially via sharing of pathogenicity islands and plasmids. This highlights the importance of surveillance and targeted intervention measures to successfully control Salmonella contamination.

Genomic epidemiology of Campylobacter jejuni associated with asymptomatic pediatric infection in the Peruvian Amazon

Campylobacter is the leading bacterial cause of gastroenteritis worldwide and its incidence is especially high in low- and middle-income countries (LMIC). Disease epidemiology in LMICs is different compared to high income countries like the USA or in Europe. Children in LMICs commonly have repeated and chronic infections even in the absence of symptoms, which can lead to deficits in early childhood development. In this study, we sequenced and characterized C. jejuni (n = 62) from a longitudinal cohort study of children under the age of 5 with and without diarrheal symptoms, and contextualized them within a global C. jejuni genome collection. Epidemiological differences in disease presentation were reflected in the genomes, specifically by the absence of some of the most common global disease-causing lineages. As in many other countries, poultry-associated strains were likely a major source of human infection but almost half of local disease cases (15 of 31) were attributable to genotypes that are rare outside of Peru. Asymptomatic infection was not limited to a single (or few) human adapted lineages but resulted from phylogenetically divergent strains suggesting an important role for host factors in the cryptic epidemiology of campylobacteriosis in LMICs.

Persistent contamination of raw milk by Campylobacter jejuni ST-883

Campylobacter jejuni has caused several campylobacteriosis outbreaks via raw milk consumption. This study reports follow-up of a milk-borne campylobacteriosis outbreak that revealed persistent C. jejuni contamination of bulk tank milk for seven months or longer. Only the outbreak-causing strain, representing sequence type (ST) 883, was isolated from milk, although other C. jejuni STs were also isolated from the farm. We hypothesized that the outbreak strain harbors features that aid its environmental transmission or survival in milk. To identify such phenotypic features, the outbreak strain was characterized for survival in refrigerated raw milk and in aerobic broth culture by plate counting and for biofilm formation on microplates by crystal violet staining and quantification. Furthermore, whole-genome sequences were studied for such genotypic features. For comparison, we characterized isolates representing other STs from the same farm and an ST-883 isolate that persisted on another dairy farm, but was not isolated from bulk tank milk. With high inocula (105 CFU/ml), ST-883 strains survived in refrigerated raw milk longer (4-6 days) than the other strains (≤3 days), but the outbreak strain showed no outperformance among ST-883 strains. This suggests that ST-883 strains may share features that aid their survival in milk, but other mechanisms are required for persistence in milk. No correlation was observed between survival in refrigerated milk and aerotolerance. The outbreak strain formed a biofilm, offering a potential explanation for persistence in milk. Whether biofilm formation was affected by pTet-like genomic element and phase-variable genes encoding capsular methyltransferase and cytochrome C551 peroxidase warrants further study. This study suggests a phenotypic target candidate for interventions and genetic markers for the phenotype, which should be investigated further with the final aim of developing control strategies against C. jejuni infections.

High-resolution sweep metagenomics using fast probabilistic inference

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.

Gene pool transmission of multidrug resistance among Campylobacter from livestock, sewage and human disease

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.

A systematic review of source attribution of human campylobacteriosis using multilocus sequence typing

BackgroundCampylobacter is a leading global cause of bacterial gastroenteritis, motivating research to identify sources of human infection. Population genetic studies have been increasingly applied to this end, mainly using multilocus sequence typing (MLST) data.ObjectivesThis review aimed to summarise approaches and findings of these studies and identify best practice lessons for this form of genomic epidemiology.MethodsWe systematically reviewed publications using MLST data to attribute human disease isolates to source. Publications were from January 2001, when this type of approach began. Searched databases included Scopus, Web of Science and PubMed. Information on samples and isolate datasets used, as well as MLST schemes and attribution algorithms employed, was obtained. Main findings were extracted, as well as any results' validation with subsequent correction for identified biases. Meta-analysis is not reported given high levels of heterogeneity.ResultsOf 2,109 studies retrieved worldwide, 25 were included, and poultry, specifically chickens, were identified as principal source of human infection. Ruminants (cattle or sheep) were consistently implicated in a substantial proportion of cases. Data sampling and analytical approaches varied, with five different attribution algorithms used. Validation such as self-attribution of isolates from known sources was reported in five publications. No publication reported adjustment for biases identified by validation.ConclusionsCommon gaps in validation and adjustment highlight opportunities to generate improved estimates in future genomic attribution studies. The consistency of chicken as the main source of human infection, across high income countries, and despite methodological variations, highlights the public health importance of this source.

Whole Genome Sequencing for Surveillance of Diphtheria in Low Incidence Settings

Corynebacterium diphtheriae (C. diphtheriae) is a relatively rare pathogen in most Western countries. While toxin producing strains can cause pharyngeal diphtheria with potentially fatal outcomes, the more common presentation is wound infections. The diphtheria toxin is encoded on a prophage and can also be carried by Corynebacterium ulcerans and Corynebacterium pseudotuberculosis. Currently, across Europe, infections are mainly diagnosed in travelers and refugees from regions where diphtheria is more endemic, patients from urban areas with poor hygiene, and intravenous drug users. About half of the cases are non-toxin producing isolates. Rapid identification of the bacterial pathogen and toxin production is a critical element of patient and outbreak management. Beside the immediate clinical management of the patient, public health agencies should be informed of toxigenic C. diphtheriae diagnoses as soon as possible. The collection of case-related epidemiological data from the patient is often challenging due to language barriers and social circumstances. However, information on patient contacts, vaccine status and travel/refugee route, where appropriate, is critical, and should be documented. In addition, isolates should be characterized using high resolution typing, in order to identify transmissions and outbreaks. In recent years, whole genome sequencing (WGS) has become the gold standard of high-resolution typing methods, allowing detailed investigations of pathogen transmissions. De-centralized sequencing strategies with redundancy in sequencing capacities, followed by data exchange may be a valuable future option, especially since WGS becomes more available and portable. In this context, the sharing of sequence data, using public available platforms, is essential. A close interaction between microbiology laboratories, treating physicians, refugee centers, social workers, and public health officials is a key element in successful management of suspected outbreaks. Analyzing bacterial isolates at reference centers may further help to provide more specialized microbiological techniques and to standardize information, but this is also more time consuming during an outbreak. Centralized communication strategies between public health agencies and laboratories helps considerably in establishing and coordinating effective surveillance and infection control. We review the current literature on high-resolution typing of C. diphtheriae and share our own experience with the coordination of a Swiss-German outbreak.

Campylobacter jejuni: collective components promoting a successful enteric lifestyle

Campylobacter jejuni is the leading cause of bacterial diarrhoeal disease in many areas of the world. The high incidence of sporadic cases of disease in humans is largely due to its prevalence as a zoonotic agent in animals, both in agriculture and in the wild. Compared with many other enteric bacterial pathogens, C. jejuni has strict growth and nutritional requirements and lacks many virulence and colonization determinants that are typically used by bacterial pathogens to infect hosts. Instead, C. jejuni has a different collection of factors and pathways not typically associated together in enteric pathogens to establish commensalism in many animal hosts and to promote diarrhoeal disease in the human population. In this Review, we discuss the cellular architecture and structure of C. jejuni, intraspecies genotypic variation, the multiple roles of the flagellum, specific nutritional and environmental growth requirements and how these factors contribute to in vivo growth in human and avian hosts, persistent colonization and pathogenesis of diarrhoeal disease.

Domestication of Campylobacter jejuni NCTC 11168

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 next generation sequencing for improving food safety: Translation into practice

Next Generation Sequencing (NGS) combined with powerful bioinformatic approaches are revolutionising food microbiology. Whole genome sequencing (WGS) of single isolates allows the most detailed comparison possible hitherto of individual strains. The two principle approaches for strain discrimination, single nucleotide polymorphism (SNP) analysis and genomic multi-locus sequence typing (MLST) are showing concordant results for phylogenetic clustering and are complementary to each other. Metabarcoding and metagenomics, applied to total DNA isolated from either food materials or the production environment, allows the identification of complete microbial populations. Metagenomics identifies the entire gene content and when coupled to transcriptomics or proteomics, allows the identification of functional capacity and biochemical activity of microbial populations. The focus of this review is on the recent use and future potential of NGS in food microbiology and on current challenges. Guidance is provided for new users, such as public health departments and the food industry, on the implementation of NGS and how to critically interpret results and place them in a broader context. The review aims to promote the broader application of NGS technologies within the food industry as well as highlight knowledge gaps and novel applications of NGS with the aim of driving future research and increasing food safety outputs from its wider use.

Source attribution of Campylobacter jejuni shows variable importance of chicken and ruminants reservoirs in non-invasive and invasive French clinical isolates

Campylobacter jejuni is the most common cause of bacterial gastroenteritis worldwide. Mainly isolated from stool samples, C. jejuni can also become invasive. C. jejuni belongs to the commensal microbiota of a number of hosts, and infection by this bacterium can sometimes be traced back to exposure to a specific source. Here we genome sequenced 200 clinical isolates (2010–2016) and analyzed them with 701 isolate genomes from human infection, chicken, ruminants and the environment to examine the relative contribution of different reservoirs to non-invasive and invasive infection in France. Host-segregating genetic markers that can discriminate C. jejuni source were used with STRUCTURE software to probabilistically attribute the source of clinical strains. A self-attribution correction step, based upon the accuracy of source apportionment within each potential reservoir, improved attribution accuracy of clinical strains and suggested an important role for ruminant reservoirs in non-invasive infection and a potentially increased contribution of chicken as a source of invasive isolates. Structured sampling of Campylobacter in the clinic and from potential reservoirs provided evidence for variation in the contribution of different infection sources over time and an important role for non-poultry reservoirs in France. This provides a basis for ongoing genomic epidemiology surveillance and targeted interventions.

Integrating Whole-Genome Sequencing Data Into Quantitative Risk Assessment of Foodborne Antimicrobial Resistance: A Review of Opportunities and Challenges

Whole-genome sequencing (WGS) will soon replace traditional phenotypic methods for routine testing of foodborne antimicrobial resistance (AMR). WGS is expected to improve AMR surveillance by providing a greater understanding of the transmission of resistant bacteria and AMR genes throughout the food chain, and therefore support risk assessment activities. At this stage, it is unclear how WGS data can be integrated into quantitative microbial risk assessment (QMRA) models and whether their integration will impact final risk estimates or the assessment of risk mitigation measures. This review explores opportunities and challenges of integrating WGS data into QMRA models that follow the Codex Alimentarius Guidelines for Risk Analysis of Foodborne AMR. We describe how WGS offers an opportunity to enhance the next-generation of foodborne AMR QMRA modeling. Instead of considering all hazard strains as equally likely to cause disease, WGS data can improve hazard identification by focusing on those strains of highest public health relevance. WGS results can be used to stratify hazards into strains with similar genetic profiles that are expected to behave similarly, e.g., in terms of growth, survival, virulence or response to antimicrobial treatment. The QMRA input distributions can be tailored to each strain accordingly, making it possible to capture the variability in the strains of interest while decreasing the uncertainty in the model. WGS also allows for a more meaningful approach to explore genetic similarity among bacterial populations found at successive stages of the food chain, improving the estimation of the probability and magnitude of exposure to AMR hazards at point of consumption. WGS therefore has the potential to substantially improve the utility of foodborne AMR QMRA models. However, some degree of uncertainty remains in relation to the thresholds of genetic similarity to be used, as well as the degree of correlation between genotypic and phenotypic profiles. The latter could be improved using a functional approach based on prediction of microbial behavior from a combination of 'omics' techniques (e.g., transcriptomics, proteomics and metabolomics). We strongly recommend that methodologies to incorporate WGS data in risk assessment be included in any future revision of the Codex Alimentarius Guidelines for Risk Analysis of Foodborne AMR.

Longitudinal Study of Shiga Toxin-Producing Escherichia coli and Campylobacter jejuni on Finnish Dairy Farms and in Raw Milk

Shiga toxin-producing Escherichia coli (STEC) and Campylobacter jejuni are notable health hazards associated with the consumption of raw milk. These bacteria may colonize the intestines of asymptomatic cattle and enter bulk tank milk via fecal contamination during milking. We studied the frequency of STEC O157:H7 and C. jejuni contamination in tank milk (n = 785) and the in-line milk filters of milking machines (n = 631) versus the frequency of isolation from cattle feces (n = 257) on three Finnish dairy farms for 1 year. Despite simultaneous isolation of STEC O157:H7 (17%) or C. jejuni (53%) from cattle, these bacteria were rarely isolated from milk filters (2% or <1%, respectively) and milk (0%). As revealed by phylogenomics, one STEC O157:H7 strain at a time was detected on each farm and persisted for ≤12 months despite rigorous hygienic measures. C. jejuni strains of a generalist sequence type (ST-883 and ST-1080) persisted in the herds for ≥11 months, and several other C. jejuni types were detected sporadically. The stx gene carried by STEC was detected more frequently from milk filters (37%) than from milk (7%), suggesting that milk filters are more suitable sampling targets for monitoring than milk. A questionnaire of on-farm practices suggested lower stx contamination of milk when major cleansing in the barn, culling, or pasturing of dairy cows was applied, while a higher average outdoor temperature was associated with higher stx contamination. Because pathogen contamination occurred despite good hygiene and because pathogen detection from milk and milk filters proved challenging, we recommend heat treatment for raw milk before consumption.IMPORTANCE The increased popularity of raw milk consumption has created demand for relaxing legislation, despite the risk of contamination by pathogenic bacteria, notably STEC and C. jejuni However, the epidemiology of these milk-borne pathogens on the herd level is still poorly understood, and data are lacking on the frequency of milk contamination on farms with cattle shedding these bacteria in their feces. This study suggests (i) that STEC contamination in milk can be reduced, but not prevented, by on-farm hygienic measures while fecal shedding is observable, (ii) that milk filters are more suitable sampling targets for monitoring than milk although pathogen detection from both sample matrices may be challenging, and (iii) that STEC and C. jejuni genotypes may persist in cattle herds for several months. The results can be utilized in developing and targeting pathogen monitoring and risk management on the farm level and contributed to the revision of Finnish legislation in 2017.

Genomic correlates of extraintestinal infection are linked with changes in cell morphology in Campylobacter jejuni

Campylobacter jejuni is the most common cause of bacterial diarrheal disease in the world. Clinical outcomes of infection can range from asymptomatic infection to life-threatening extraintestinal infections. This variability in outcomes for infected patients has raised questions as to whether genetic differences between C. jejuni isolates contribute to their likelihood of causing severe disease. In this study, we compare the genomes of ten C. jejuni isolates that were implicated in extraintestinal infections with reference gastrointestinal isolates, in order to identify unusual patterns of sequence variation associated with infection outcome. We identified a collection of genes that display a higher burden of uncommon mutations in invasive isolates compared with gastrointestinal close relatives, including some that have been previously linked to virulence and invasiveness in C. jejuni. Among the top genes identified were mreB and pgp1, which are both involved in determining cell shape. Electron microscopy confirmed morphological differences in isolates carrying unusual sequence variants of these genes, indicating a possible relationship between extraintestinal infection and changes in cell morphology.

The function, biogenesis and regulation of the electron transport chains in Campylobacter jejuni: New insights into the bioenergetics of a major food-borne pathogen

Campylobacter jejuni is a zoonotic Epsilonproteobacterium that grows in the gastrointestinal tract of birds and mammals, and is the most frequent cause of food-borne bacterial gastroenteritis worldwide. As an oxygen-sensitive microaerophile, C. jejuni has to survive high environmental oxygen tensions, adapt to oxygen limitation in the host intestine and resist host oxidative attack. Despite its small genome size, C. jejuni is a versatile and metabolically active pathogen, with a complex and highly branched set of respiratory chains allowing the use of a wide range of electron donors and alternative electron acceptors in addition to oxygen, including fumarate, nitrate, nitrite, tetrathionate and N- or S-oxides. Several novel enzymes participate in these electron transport chains, including a tungsten containing formate dehydrogenase, a Complex I that uses flavodoxin and not NADH, a periplasmic facing fumarate reductase and a cytochrome c tetrathionate reductase. This review presents an updated description of the composition and bioenergetics of these various respiratory chains as they are currently understood, including recent work that gives new insights into energy conservation during electron transport to various alternative electron acceptors. The regulation of synthesis and assembly of the electron transport chains is also discussed. A deeper appreciation of the unique features of the respiratory systems of C. jejuni may be helpful in informing strategies to control this important pathogen.

PhasomeIt: an 'omics' approach to cataloguing the potential breadth of phase variation in the genus Campylobacter

Hypermutable simple sequence repeats (SSRs) are drivers of phase variation (PV) whose stochastic, high-frequency, reversible switches in gene expression are a common feature of several pathogenic bacterial species, including the human pathogen Campylobacter jejuni. Here we examine the distribution and conservation of known and putative SSR-driven phase variable genes - the phasome - in the genus Campylobacter. PhasomeIt, a new program, was specifically designed for rapid identification of SSR-mediated PV. This program detects the location, type and repeat number of every SSR. Each SSR is linked to a specific gene and its putative expression state. Other outputs include conservation of SSR-driven phase-variable genes and the 'core phasome' - the minimal set of PV genes in a phylogenetic grouping. Analysis of 77 complete Campylobacter genome sequences detected a 'core phasome' of conserved PV genes in each species and a large number of rare PV genes with few, or no, homologues in other genome sequences. Analysis of a set of partial genome sequences, with food-chain-associated metadata, detected evidence of a weak link between phasome and source host for disease-causing isolates of sequence type (ST)-828 but not the ST-21 or ST-45 complexes. Investigation of the phasomes in the genus Campylobacter provided evidence of overlapping but distinctive mechanisms of PV-mediated adaptation to specific niches. This suggests that the phasome could be involved in host adaptation and spread of campylobacters. Finally, this tool is malleable and will have utility for studying the distribution and genic effects of other repetitive elements in diverse bacterial species.