Agricultural intensification and the evolution of host specialism in the enteric pathogen Campylobacter jejuni.

Modern agriculture has dramatically changed the distribution of animal species on Earth. Changes to host ecology have a major impact on the microbiota, potentially increasing the risk of zoonotic pathogens being transmitted to humans, but the impact of intensive livestock production on host-associated bacteria has rarely been studied. Here, we use large isolate collections and comparative genomics techniques, linked to phenotype studies, to understand the timescale and genomic adaptations associated with the proliferation of the most common food-born bacterial pathogen (Campylobacter jejuni) in the most prolific agricultural mammal (cattle). Our findings reveal the emergence of cattle specialist C. jejuni lineages from a background of host generalist strains that coincided with the dramatic rise in cattle numbers in the 20th century. Cattle adaptation was associated with horizontal gene transfer and significant gene gain and loss. This may be related to differences in host diet, anatomy, and physiology, leading to the proliferation of globally disseminated cattle specialists of major public health importance. This work highlights how genomic plasticity can allow important zoonotic pathogens to exploit altered niches in the face of anthropogenic change and provides information for mitigating some of the risks posed by modern agricultural systems.

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.

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.

Local genes for local bacteria: Evidence of allopatry in the genomes of transatlantic Campylobacter populations.

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.

Core Genome Multilocus Sequence Typing Scheme for Stable, Comparative Analyses of Campylobacter jejuni and C. coli Human Disease Isolates.

Human campylobacteriosis, caused by Campylobacter jejuni and C. coli, remains a leading cause of bacterial gastroenteritis in many countries, but the epidemiology of campylobacteriosis outbreaks remains poorly defined, largely due to limitations in the resolution and comparability of isolate characterization methods. Whole-genome sequencing (WGS) data enable the improvement of sequence-based typing approaches, such as multilocus sequence typing (MLST), by substantially increasing the number of loci examined. A core genome MLST (cgMLST) scheme defines a comprehensive set of those loci present in most members of a bacterial group, balancing very high resolution with comparability across the diversity of the group. Here we propose a set of 1,343 loci as a human campylobacteriosis cgMLST scheme (v1.0), the allelic profiles of which can be assigned to core genome sequence types. The 1,343 loci chosen were a subset of the 1,643 loci identified in the reannotation of the genome sequence of C. jejuni isolate NCTC 11168, chosen as being present in >95% of draft genomes of 2,472 representative United Kingdom campylobacteriosis isolates, comprising 2,207 (89.3%) C. jejuni isolates and 265 (10.7%) C. coli isolates. Validation of the cgMLST scheme was undertaken with 1,478 further high-quality draft genomes, containing 150 or fewer contiguous sequences, from disease isolate collections: 99.5% of these isolates contained ≥95% of the 1,343 cgMLST loci. In addition to the rapid and effective high-resolution analysis of large numbers of diverse isolates, the cgMLST scheme enabled the efficient identification of very closely related isolates from a well-defined single-source campylobacteriosis outbreak.

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

Campylobacter jejuni is a major cause of bacterial gastroenteritis worldwide, primarily associated with the consumption of contaminated poultry. C. jejuni lineages vary in host range and prevalence in human infection, suggesting differences in survival throughout the poultry processing chain. From 7343 MLST-characterised isolates, we sequenced 600 C. jejuni and C. coli isolates from various stages of poultry processing and clinical cases. A genome-wide association study (GWAS) in C. jejuni ST-21 and ST-45 complexes identified genetic elements over-represented in clinical isolates that increased in frequency throughout the poultry processing chain. Disease-associated SNPs were distinct in these complexes, sometimes organised in haplotype blocks. The function of genes containing associated elements was investigated, demonstrating roles for cj1377c in formate metabolism, nuoK in aerobic survival and oxidative respiration, and cj1368-70 in nucleotide salvage. This work demonstrates the utility of GWAS for investigating transmission in natural zoonotic pathogen populations and provides evidence that major C. jejuni lineages have distinct genotypes associated with survival, within the host specific niche, from farm to fork.

Exploiting Bacterial Whole-Genome Sequencing Data for Evaluation of Diagnostic Assays: Campylobacter Species Identification as a Case Study.

The application of whole-genome sequencing (WGS) to problems in clinical microbiology has had a major impact on the field. Clinical laboratories are now using WGS for pathogen identification, antimicrobial susceptibility testing, and epidemiological typing. WGS data also represent a valuable resource for the development and evaluation of molecular diagnostic assays, which continue to play an important role in clinical microbiology. To demonstrate this application of WGS, this study used publicly available genomic data to evaluate a duplex real-time PCR (RT-PCR) assay that targets mapA and ceuE for the detection of Campylobacter jejuni and Campylobacter coli, leading global causes of bacterial gastroenteritis. In silico analyses of mapA and ceuE primer and probe sequences from 1,713 genetically diverse C. jejuni and C. coli genomes, supported by RT-PCR testing, indicated that the assay was robust, with 1,707 (99.7%) isolates correctly identified. The high specificity of the mapA-ceuE assay was the result of interspecies diversity and intraspecies conservation of the target genes in C. jejuni and C. coli Rare instances of a lack of specificity among C. coli isolates were due to introgression in mapA or sequence diversity in ceuE The results of this study illustrate how WGS can be exploited to evaluate molecular diagnostic assays by using publicly available data, online databases, and open-source software.

Rapid host switching in generalist Campylobacter strains erodes the signal for tracing human infections.

Campylobacter jejuni and Campylobacter coli are the biggest causes of bacterial gastroenteritis in the developed world, with human infections typically arising from zoonotic transmission associated with infected meat. Because Campylobacter is not thought to survive well outside the gut, host-associated populations are genetically isolated to varying degrees. Therefore, the likely origin of most strains can be determined by host-associated variation in the genome. This is instructive for characterizing the source of human infection. However, some common strains, notably isolates belonging to the ST-21, ST-45 and ST-828 clonal complexes, appear to have broad host ranges, hindering source attribution. Here whole-genome sequencing has the potential to reveal fine-scale genetic structure associated with host specificity. We found that rates of zoonotic transmission among animal host species in these clonal complexes were so high that the signal of host association is all but obliterated, estimating one zoonotic transmission event every 1.6, 1.8 and 12 years in the ST-21, ST-45 and ST828 complexes, respectively. We attributed 89% of clinical cases to a chicken source, 10% to cattle and 1% to pig. Our results reveal that common strains of C. jejuni and C. coli infectious to humans are adapted to a generalist lifestyle, permitting rapid transmission between different hosts. Furthermore, they show that the weak signal of host association within these complexes presents a challenge for pinpointing the source of clinical infections, underlining the view that whole-genome sequencing, powerful though it is, cannot substitute for intensive sampling of suspected transmission reservoirs.

Wild bird-associated Campylobacter jejuni isolates are a consistent source of human disease, in Oxfordshire, United Kingdom.

The contribution of wild birds as a source of human campylobacteriosis was investigated in Oxfordshire, United Kingdom (UK) over a 10 year period. The probable origin of human Campylobacter jejuni genotypes, as described by multilocus sequence typing, was estimated by comparison with reference populations of isolates from farm animals and five wild bird families, using the STRUCTURE algorithm. Wild bird-attributed isolates accounted for between 476 (2.1%) and 543 (3.5%) cases annually. This proportion did not vary significantly by study year (P = 0.934) but varied seasonally, with wild bird-attributed genotypes comprising a greater proportion of isolates during warmer compared with cooler months (P = 0.003). The highest proportion of wild bird-attributed illness occurred in August (P < 0.001), with a significantly lower proportion in November (P = 0.018). Among genotypes attributed to specific groups of wild birds, seasonality was most apparent for Turdidae-attributed isolates, which were absent during cooler, winter months. This study is consistent with some wild bird species representing a persistent source of campylobacteriosis, and contributing a distinctive seasonal pattern to disease burden. If Oxfordshire is representative of the UK as a whole in this respect, these data suggest that the national burden of wild bird-attributed isolates could be in the order of 10,000 annually.

Defining the estimated core genome of bacterial populations using a Bayesian decision model.

The bacterial core genome is of intense interest and the volume of whole genome sequence data in the public domain available to investigate it has increased dramatically. The aim of our study was to develop a model to estimate the bacterial core genome from next-generation whole genome sequencing data and use this model to identify novel genes associated with important biological functions. Five bacterial datasets were analysed, comprising 2096 genomes in total. We developed a Bayesian decision model to estimate the number of core genes, calculated pairwise evolutionary distances (p-distances) based on nucleotide sequence diversity, and plotted the median p-distance for each core gene relative to its genome location. We designed visually-informative genome diagrams to depict areas of interest in genomes. Case studies demonstrated how the model could identify areas for further study, e.g. 25% of the core genes with higher sequence diversity in the Campylobacter jejuni and Neisseria meningitidis genomes encoded hypothetical proteins. The core gene with the highest p-distance value in C. jejuni was annotated in the reference genome as a putative hydrolase, but further work revealed that it shared sequence homology with beta-lactamase/metallo-beta-lactamases (enzymes that provide resistance to a range of broad-spectrum antibiotics) and thioredoxin reductase genes (which reduce oxidative stress and are essential for DNA replication) in other C. jejuni genomes. Our Bayesian model of estimating the core genome is principled, easy to use and can be applied to large genome datasets. This study also highlighted the lack of knowledge currently available for many core genes in bacterial genomes of significant global public health importance.

Real-time genomic epidemiological evaluation of human Campylobacter isolates by use of whole-genome multilocus sequence typing.

Sequence-based typing is essential for understanding the epidemiology of Campylobacter infections, a major worldwide cause of bacterial gastroenteritis. We demonstrate the practical and rapid exploitation of whole-genome sequencing to provide routine definitive characterization of Campylobacter jejuni and Campylobacter coli for clinical and public health purposes. Short-read data from 384 Campylobacter clinical isolates collected over 4 months in Oxford, United Kingdom, were assembled de novo. Contigs were deposited at the pubMLST.org/campylobacter website and automatically annotated for 1,667 loci. Typing and phylogenetic information was extracted and comparative analyses were performed for various subsets of loci, up to the level of the whole genome, using the Genome Comparator and Neighbor-net algorithms. The assembled sequences (for 379 isolates) were diverse and resembled collections from previous studies of human campylobacteriosis. Small subsets of very closely related isolates originated mainly from repeated sampling from the same patients and, in one case, likely laboratory contamination. Much of the within-patient variation occurred in phase-variable genes. Clinically and epidemiologically informative data can be extracted from whole-genome sequence data in real time with straightforward, publicly available tools. These analyses are highly scalable, are transparent, do not require closely related genome reference sequences, and provide improved resolution (i) among Campylobacter clonal complexes and (ii) between very closely related isolates. Additionally, these analyses rapidly differentiated unrelated isolates, allowing the detection of single-strain clusters. The approach is widely applicable to analyses of human bacterial pathogens in real time in clinical laboratories, with little specialist training required.

A longitudinal 6-year study of the molecular epidemiology of clinical campylobacter isolates in Oxfordshire, United kingdom.

Temporal and seasonal trends in Campylobacter genotypes causing human gastroenteritis were investigated in a 6-year study of 3,300 recent isolates from Oxfordshire, United Kingdom. Genotypes (sequence types [ST]) were defined using multilocus sequence typing and assigned to a clonal complex (a cluster of related strains that share four or more identical alleles with a previously defined central genotype). A previously undescribed clonal complex (ST-464) was identified which, together with ST-42, ST-45, and ST-52 complexes, showed increasing incidence. Concurrently, the incidence of ST-574, ST-607, and ST-658 complexes declined. The relative frequencies of three clonal complexes (ST-45, ST-283, and ST-42) peaked during summer and those of two (ST-353 and ST-403) peaked during winter. Nine clonal complexes (ST-22, ST-45, ST-48, ST-61, ST-257, ST-283, ST-403, ST-658, and ST-677) were significantly associated with ciprofloxacin sensitivity (P < 0.05). Seven clonal complexes (ST-49, ST-206, ST-354, ST-446, ST-460, ST-464, and ST-607) were associated with ciprofloxacin resistance (P < 0.05). Clonal complexes exhibited changing incidence and differences in seasonality and antibiotic resistance phenotype. These data also demonstrated that detailed surveillance at a single site captures information which reflects that observed nationally.

Ribosomal multilocus sequence typing: universal characterization of bacteria from domain to strain.

No single genealogical reconstruction or typing method currently encompasses all levels of bacterial diversity, from domain to strain. We propose ribosomal multilocus sequence typing (rMLST), an approach which indexes variation of the 53 genes encoding the bacterial ribosome protein subunits (rps genes), as a means of integrating microbial genealogy and typing. As with multilocus sequence typing (MLST), rMLST employs curated reference sequences to identify gene variants efficiently and rapidly. The rps loci are ideal targets for a universal characterization scheme as they are: (i) present in all bacteria; (ii) distributed around the chromosome; and (iii) encode proteins which are under stabilizing selection for functional conservation. Collectively, the rps loci exhibit variation that resolves bacteria into groups at all taxonomic and most typing levels, providing significantly more resolution than 16S small subunit rRNA gene phylogenies. A web-accessible expandable database, comprising whole-genome data from more than 1900 bacterial isolates, including 28 draft genomes assembled de novo from the European Bioinformatics Institute (EBI) sequence read archive, has been assembled. The rps gene variation catalogued in this database permits rapid and computationally non-intensive identification of the phylogenetic position of any bacterial sequence at the domain, phylum, class, order, family, genus, species and strain levels. The groupings generated with rMLST data are consistent with current nomenclature schemes and independent of the clustering algorithm used. This approach is applicable to the other domains of life, potentially providing a rational and universal approach to the classification of life that is based on one of its fundamental features, the translation mechanism.

Temporal variation and host association in the Campylobacter population in a longitudinal ruminant farm study.

Campylobacter jejuni and C. coli were quantified and typed, using multilocus sequence typing (MLST), from fecal samples collected from a mixed cattle and sheep farm during summer. Cattle had a significantly higher prevalence than sheep (21.9% [74/338] and 14.0% [30/214], respectively), but both decreased over time. There were no differences in the average Campylobacter concentrations shed by cattle (600 CFU g(-1)) and sheep (820 CFU g(-1)), although sheep did show a significant temporal reduction in the number of Campylobacter organisms shed in their feces. A total of 21 different sequence types (STs) (97.7% C. jejuni, 2.3% C. coli) were isolated from cattle, and 9 different STs were isolated from sheep (40.6% C. jejuni, 59.4% C. coli). The Campylobacter population in cattle was relatively stable, and the frequencies of genotypes isolated showed little temporal variation. However, the composition of subtypes isolated from sheep did show significant temporal differences. The cattle and sheep consistently showed significant differences in their carriage of Campylobacter species, STs, and CCs despite the fact that both were exposed to the same farming environment. This work has highlighted the patterns of a Campylobacter population on a ruminant farm by identifying the existence of both temporal and between-host variations.

Host association of Campylobacter genotypes transcends geographic variation.

Genetic attribution of bacterial genotypes has become a major tool in the investigation of the epidemiology of campylobacteriosis and has implicated retail chicken meat as the major source of human infection in several countries. To investigate the robustness of this approach to the provenance of the reference data sets used, a collection of 742 Campylobacter jejuni and 261 Campylobacter coli isolates obtained from United Kingdom-sourced chicken meat was established and typed by multilocus sequence typing. Comparative analyses of the data with those from other isolates sourced from a variety of host animals and countries were undertaken by genetic attribution, genealogical, and population genetic approaches. The genotypes from the United Kingdom data set were highly diverse, yet structured into sequence types, clonal complexes, and genealogical groups very similar to those seen in chicken isolates from the Netherlands, the United States, and Senegal, but more distinct from isolates obtained from ruminant, swine, and wild bird sources. Assignment analyses consistently grouped isolates from different host animal sources regardless of geographical source; these associations were more robust than geographic associations across isolates from three continents. We conclude that, notwithstanding the high diversity of these pathogens, there is a strong signal of association of multilocus genotypes with particular hosts, which is greater than the geographic signal. These findings are consistent with local and international transmission of host-associated lineages among food animal species and provide a foundation for further improvements in genetic attribution.

Comparative genotyping of Campylobacter jejuni strains from patients with Guillain-Barré syndrome in Bangladesh.

BACKGROUND: Campylobacter jejuni is a common cause of acute gastroenteritis and is associated with post-infectious neuropathies such as the Guillain-Barré syndrome (GBS) and the Miller Fisher syndrome (MFS). We here present comparative genotyping of 49 C. jejuni strains from Bangladesh that were recovered from patients with enteritis or GBS. All strains were serotyped and analyzed by lipo-oligosaccharide (LOS) genotyping, amplified fragment length polymorphism (AFLP) analysis, multilocus sequence typing (MLST), and pulsed-field gel electrophoresis (PFGE). METHODOLOGY/PRINCIPAL FINDINGS: C. jejuni HS:23 was a predominant serotype among GBS patients (50%), and no specific serotype was significantly associated with GBS compared to enteritis. PCR screening showed that 38/49 (78%) of strains could be assigned to LOS classes A, B, C, or E. The class A locus (4/7 vs 3/39; p<0.01) was significantly associated in the GBS-related strains as compared to enteritis strains. All GBS/oculomotor related strains contained the class B locus; which was also detected in 46% of control strains. Overlapping clonal groups were defined by MLST, AFLP and PFGE for strains from patients with gastroenteritis and GBS. MLST defined 22 sequence types (STs) and 7 clonal complexes including 7 STs not previously identified (ST-3742, ST-3741, ST-3743, ST-3748, ST-3968, ST-3969 and ST-3970). C. jejuni HS:23 strains from patients with GBS or enteritis were clonal and all strains belonged to ST-403 complex. Concordance between LOS class B and ST-403 complex was revealed. AFLP defined 25 different types at 90% similarity. The predominant AFLP type AF-20 coincided with the C. jejuni HS:23 and ST-403 complex. CONCLUSION/SIGNIFICANCE: LOS genotyping, MLST, AFLP and PFGE helped to identify the HS:23 strains from GBS or enteritis patients as clonal. Overall, genotypes exclusive for enteritis or for GBS-related strains were not obtained although LOS class A was significantly associated with GBS strains. Particularly, the presence of a clonal and putative neuropathogenic C. jejuni HS:23 serotype may contribute to the high prevalence of C. jejuni related GBS in Bangladesh.

Extended sequence typing of Campylobacter spp., United Kingdom.

Supplementing Campylobacter spp. multilocus sequence typing with nucleotide sequence typing of 3 antigen genes increased the discriminatory index achieved from 0.975 to 0.992 among 620 clinical isolates from Oxfordshire, United Kingdom. This enhanced typing scheme enabled identification of clusters and retained data required for long-range epidemiologic comparisons of isolates.