Jejuni Genomes Research Articles

Machine learning to attribute the source of Campylobacter infections in the United States: A retrospective analysis of national surveillance data

Integrating pathogen genomic surveillance with bioinformatics can enhance public health responses by identifying risk and guiding interventions. This study focusses on the two predominant Campylobacter species, which are commonly found in the gut of birds and mammals and often infect humans via contaminated food. Rising incidence and antimicrobial resistance (AMR) are a global concern, and there is an urgent need to quantify the main routes to human infection. During routine US national surveillance (2009-2019), 8856 Campylobacter genomes from human infections and 16,703 from possible sources were sequenced. Using machine learning and probabilistic models, we target genetic variation associated with host adaptation to attribute the source of human infections and estimate the importance of different disease reservoirs. Poultry was identified as the primary source of human infections, responsible for an estimated 68% of cases, followed by cattle (28%), and only a small contribution from wild birds (3%) and pork sources (1%). There was also evidence of an increase in multidrug resistance, particularly among isolates attributed to chickens. National surveillance and source attribution can guide policy, and our study suggests that interventions targeting poultry will yield the greatest reductions in campylobacteriosis and spread of AMR in the US. All sequence reads were uploaded and shared on NCBI's Sequence Read Archive (SRA) associated with BioProjects; PRJNA239251 (CDC / PulseNet surveillance), PRJNA287430 (FSIS surveillance), PRJNA292668 & PRJNA292664 (NARMS) and PRJNA258022 (FDA surveillance). Publicly available genomes, including reference genomes and isolates sampled worldwide from wild birds are associated with BioProject accessions: PRJNA176480, PRJNA177352, PRJNA342755, PRJNA345429, PRJNA312235, PRJNA415188, PRJNA524300, PRJNA528879, PRJNA529798, PRJNA575343, PRJNA524315 and PRJNA689604. Contiguous assemblies of all genome sequences compared are available at Mendeley data (assembled C. coli genomes doi: 10.17632/gxswjvxyh3.1; assembled C. jejuni genomes doi: 10.17632/6ngsz3dtbd.1) and individual project and accession numbers can be found in Supplementary tables S1 and S2, which also includes pubMLST identifiers for assembled genomes. Figshare (10.6084/m9.figshare.20279928). Interactive phylogenies are hosted on microreact separately for C. jejuni (https://microreact.org/project/pascoe-us-cjejuni) and C. coli (https://microreact.org/project/pascoe-us-ccoli).

RefMLST: reference-based multilocus sequence typing enables universal bacterial typing

BackgroundCommonly used approaches for genomic investigation of bacterial outbreaks, including SNP and gene-by-gene approaches, are limited by the requirement for background genomes and curated allele schemes, respectively. As a result, they only work on a select subset of known organisms, and fail on novel or less studied pathogens. We introduce refMLST, a gene-by-gene approach using the reference genome of a bacterium to form a scalable, reproducible and robust method to perform outbreak investigation.ResultsWhen applied to multiple outbreak causing bacteria including 1263 Salmonella enterica, 331 Yersinia enterocolitica and 6526 Campylobacter jejuni genomes, refMLST enabled consistent clustering, improved resolution, and faster processing in comparison to commonly used tools like chewieSnake.ConclusionsrefMLST is a novel multilocus sequence typing approach that is applicable to any bacterial species with a public reference genome, does not require a curated scheme, and automatically accounts for genetic recombination.Availability and implementation: refMLST is freely available for academic use at https://bugseq.com/academic.

Capturing clinically relevant Campylobacter attributes through direct whole genome sequencing of stool.

Campylobacter is the leading bacterial cause of infectious intestinal disease, but the pathogen typically accounts for a very small proportion of the overall stool microbiome in each patient. Diagnosis is even more difficult due to the fastidious nature of Campylobacter in the laboratory setting. This has, in part, driven a change in recent years, from culture-based to rapid PCR-based diagnostic assays which have improved diagnostic detection, whilst creating a knowledge gap in our clinical and epidemiological understanding of Campylobacter genotypes - no isolates to sequence. In this study, direct metagenomic sequencing approaches were used to assess the possibility of replacing genome sequences with metagenome sequences; metagenomic sequencing outputs were used to describe clinically relevant attributes of Campylobacter genotypes. A total of 37 diarrhoeal stool samples with Campylobacter and five samples with an unknown pathogen result were collected and processed with and without filtration, DNA was extracted, and metagenomes were sequenced by short-read sequencing. Culture-based methods were used to validate Campylobacter metagenome-derived genome (MDG) results. Sequence output metrics were assessed for Campylobacter genome quality and accuracy of characterization. Of the 42 samples passing quality checks for analysis, identification of Campylobacter to the genus and species level was dependent on Campylobacter genome read count, coverage and genome completeness. A total of 65% (24/37) of samples were reliably identified to the genus level through Campylobacter MDG, 73% (27/37) by culture and 97% (36/37) by qPCR. The Campylobacter genomes with a genome completeness of over 60% (n=21) were all accurately identified at the species level (100%). Of those, 72% (15/21) were identified to sequence types (STs), and 95% (20/21) accurately identified antimicrobial resistance (AMR) gene determinants. Filtration of stool samples enhanced Campylobacter MDG recovery and genome quality metrics compared to the corresponding unfiltered samples, which improved the identification of STs and AMR profiles. The phylogenetic analysis in this study demonstrated the clustering of the metagenome-derived with culture-derived genomes and revealed the reliability of genomes from direct stool sequencing. Furthermore, Campylobacter genome spiking percentages ranging from 0 to 2% total metagenome abundance in the ONT MinION sequencer, configured to adaptive sequencing, exhibited better assembly quality and accurate identification of STs, particularly in the analysis of metagenomes containing 2 and 1% of Campylobacter jejuni genomes. Direct sequencing of Campylobacter from stool samples provides clinically relevant and epidemiologically important genomic information without the reliance on cultured genomes.

Five centuries of genome evolution and multi-host adaptation of Campylobacter jejuni in Brazil.

Consumption of raw, undercooked or contaminated animal food products is a frequent cause of Campylobacter jejuni infection. Brazil is the world's third largest producer and a major exporter of chicken meat, yet population-level genomic investigations of C. jejuni in the country remain scarce. Analysis of 221 C. jejuni genomes from Brazil shows that the overall core and accessory genomic features of C. jejuni are influenced by the identity of the human or animal source. Of the 60 sequence types detected, ST353 is the most prevalent and consists of samples from chicken and human sources. Notably, we identified the presence of diverse bla genes from the OXA-61 and OXA-184 families that confer beta-lactam resistance as well as the operon cmeABCR related to multidrug efflux pump, which contributes to resistance against tetracyclines, macrolides and quinolones. Based on limited data, we estimated the most recent common ancestor of ST353 to the late 1500s, coinciding with the time the Portuguese first arrived in Brazil and introduced domesticated chickens into the country. We identified at least two instances of ancestral chicken-to-human infections in ST353. The evolution of C. jejuni in Brazil was driven by the confluence of clinically relevant genetic elements, multi-host adaptation and clonal population growth that coincided with major socio-economic changes in poultry farming.

Insights on the genomic diversity, virulence and resistance profile of a Campylobacter jejuni strain isolated from a hospitalized patient in Brazil.

Campylobacteriosis is currently recognized as one of the major causes of foodborne bacterial diseases worldwide. In Brazil, there is insufficient data to estimate the impact of Campylobacter in public health. The aim of this present study was to characterize a C. jejuni CJ-HBSJRP strain isolated from a hospitalized patient in Brazil by its ability to invade human Caco-2 epithelial cells, to survive in U937 human macrophages, and to assess its phenotypic antimicrobial resistance profile. In addition, prophages, virulence and antimicrobial resistance genes were search using whole-genome sequencing data. The genetic relatedness was evaluated by MLST and cgMLST analysis by comparison with 29 other C. jejuni genomes isolated from several countries. The CJ-HBSJRP strain showed an invasion percentage of 50% in Caco-2 polarized cells, 37.5% of survivability in U937 cells and was phenotypically resistant to ampicillin, ciprofloxacin and nalidixic acid. A total of 94 virulence genes related to adherence, biofilm, chemotaxis, immune modulation, invasion process, metabolism, motility and toxin were detected. The resistance genes blaOXA-605 (blaOXA-61), cmeB and mutations in the QRDR region of gyrA were also found and none prophages were detected. The MLST analysis showed 23 different STs among the strains studied. Regarding cgMLST analysis, the CJ-HBSJRP strain was genetically distinct and did not group closely to any other isolate. The results obtained reinforce the pathogenic potential of the CJHBSJRP strain and highlighted the need for more careful attention to Campylobacter spp. infections in Brazil since this pathogen has been the most commonly reported zoonosis in several countries worldwide.

Genomic resistant determinants of multidrug-resistant Campylobacter spp. isolates in Peru

Antimicrobial resistant (AMR) Campylobacter is a global health threat; however, there is limited information on genomic determinants of resistance in low- and middle-income countries. We evaluated genomic determinants of AMR using a collection of whole genome sequenced Campylobacter jejuni and C. coli isolates from Iquitos, Peru. Campylobacter isolates from two paediatric cohort studies enriched with isolates that demonstrated resistance to ciprofloxacin and azithromycin were sequenced and mined for AMR determinants. The gyrA mutation leading to the Thr86Ile amino acid change was the only gyrA mutation associated with fluoroquinolone resistance identified. The A2075G mutation in 23S rRNA was present, but three other 23S rRNA mutations previously associated with macrolide resistance were not identified. A resistant-enhancing variant of the cmeABC efflux pump genotype (RE-cmeABC) was identified in 36.1% (35/97) of C. jejuni genomes and 17.9% (12/67) of C. coli genomes. Mutations identified in the CmeR-binding site, an inverted repeat sequence in the cmeABC promoter region that increases expression of the operon, were identified in 24/97 C. jejuni and 14/67 C. coli genomes. The presence of these variants, in addition to RE-cmeABC, was noted in 18 of the 24 C. jejuni and 9 of the 14 C. coli genomes. Both RE-cmeABC and mutations in the CmeR-binding site were strongly associated with the MDR phenotype in C. jejuni and C. coli. This is the first report of RE-cmeABC in Peru and suggests it is a major driver of resistance to the principal therapies used to treat human campylobacteriosis in this setting.

The Distribution of Campylobacter jejuni Virulence Genes in Genomes Worldwide Derived from the NCBI Pathogen Detection Database.

Campylobacter jejuni (C. jejuni) is responsible for 80% of human campylobacteriosis and is the leading cause of gastroenteritis globally. The relevant public health risks of C. jejuni are caused by particular virulence genes encompassing its virulome. We analyzed 40,371 publicly available genomes of C. jejuni deposited in the NCBI Pathogen Detection Database, combining their epidemiologic metadata with an in silico bioinformatics analysis to increase our current comprehension of their virulome from a global perspective. The collection presented a virulome composed of 126 identified virulence factors that were grouped in three clusters representing the accessory, the softcore, and the essential core genes according to their prevalence within the genomes. The multilocus sequence type distribution in the genomes was also investigated. An unexpected low prevalence of the full-length flagellin flaA and flaB locus of C. jejuni genomes was revealed, and an essential core virulence gene repertoire prevalent in more than 99.99% of genomes was identified. Altogether, this is a pioneer study regarding Campylobacter jejuni that has compiled a significant amount of data about the Multilocus Sequence Type and virulence factors concerning their global prevalence and distribution over this database.

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 blaOXA–61 and tet(O)-tet(O/W/O) genes and the gyrA (T86I) SNP, was found from 2001 to date in Europe.

Host-Dependent Clustering of Campylobacter Strains From Small Mammals in Finland.

Small mammals are known to carry Campylobacter spp.; however, little is known about the genotypes and their role in human infections. We studied intestinal content from small wild mammals collected in their natural habitats in Finland in 2010–2017, and in close proximity to 40 pig or cattle farms in 2017. The animals were trapped using traditional Finnish metal snap traps. Campylobacter spp. were isolated from the intestinal content using direct plating on mCCDA. A total of 19% of the captured wild animals (n = 577) and 41% of the pooled farm samples (n = 227) were positive for C. jejuni, which was the only Campylobacter species identified. The highest prevalence occurred in yellow-necked mice (Apodemus flavicollis) and bank voles (Myodes glareolus) which carried Campylobacter spp. in 66.3 and 63.9% of the farm samples and 41.5 and 24.4% of individual animals trapped from natural habitats, respectively. Interestingly, all house mouse (Mus musculus) and shrew (Sorex spp.) samples were negative for Campylobacter spp. C. jejuni isolates (n = 145) were further characterized by whole-genome sequencing. Core genome multilocus sequence typing (cgMLST) clustering showed that mouse and vole strains were separated from the rest of the C. jejuni population (636 and 671 allelic differences, 94 and 99% of core loci, respectively). Very little or no alleles were shared with C. jejuni genomes described earlier from livestock or human isolates. FastANI results further indicated that C. jejuni strains from voles are likely to represent a new previously undescribed species or subspecies of Campylobacter. Core-genome phylogeny showed that there was no difference between isolates originating from the farm and wild captured animals. Instead, the phylogeny followed the host species-association. There was some evidence (one strain each) of livestock-associated C. jejuni occurring in a farm-caught A. flavicollis and a brown rat (Rattus norvegicus), indicating that although small mammals may not be the original reservoir of Campylobacter colonizing livestock, they may sporadically carry C. jejuni strains occurring mainly in livestock and be associated with disease in humans.

Whole-Genome Sequencing and Bioinformatic Analysis of Environmental, Agricultural, and Human Campylobacter jejuni Isolates From East Tennessee.

As a leading cause of bacterial-derived gastroenteritis worldwide, Campylobacter jejuni has a significant impact on human health in both the developed and developing worlds. Despite its prevalence as a human pathogen, the source of these infections remains poorly understood due to the mutation frequency of the organism and past limitations of whole genome analysis. Recent advances in both whole genome sequencing and computational methods have allowed for the high-resolution analysis of intraspecies diversity, leading multiple groups to postulate that these approaches may be used to identify the sources of Campylobacter jejuni infection. To address this hypothesis, our group conducted a regionally and temporally restricted sampling of agricultural and environmental Campylobacter sources and compared isolated C. jejuni genomes to those that caused human infections in the same region during the same time period. Through a network analysis comparing genomes from various sources, we found that human C. jejuni isolates clustered with those isolated from cattle and chickens, indicating these as potential sources of human infection in the region.

Diversity of transducer-like proteins (Tlps) in Campylobacter.

Campylobacter transducer-like proteins (Tlps), also known as methyl-accepting chemotaxis proteins (MCPs), are associated with virulence as well as niche and host adaptation. While functional attributes of these proteins are being elucidated, little has been published regarding their sequence diversity or chromosomal locations and context, although they appear to define invertible regions within Campylobacter jejuni genomes. Genome assemblies for several species of Campylobacter were obtained from the publicly available NCBI repositories. Genomes from all isolates were obtained from GenBank and assessed for Tlp content, while data from isolates with complete, finished genomes were used to determine the identity of Tlps as well as the gene content of putative invertible elements (IEs) in C. jejuni (Cj) and C. coli (Cc). Tlps from several Campylobacter species were organized into a nomenclature system and novel Tlps were defined and named for Cj and Cc. The content of Tlps appears to be species-specific, though diverse within species. Cj and Cc carried overlapping, related Tlp content, as did the three C. fetus subspecies. Tlp1 was detected in 88% of Cj isolates and approximately 43% of Cc, and was found in a different conserved chromosomal location and genetic context in each species. Tlp1 and Tlp 3 predominated in genomes from Cj whereas other Tlps were detected less frequently. Tlp13 and Tlp20 predominated in genomes from Cc while some Cj/Cc Tlps were not detected at all. Tlps 2–4 and 11–20 were less frequently detected and many showed sequence heterogeneity that could affect substrate binding, signal transduction, or both. Tlps other than Tlp1, 7, and 10 had substantial sequence identity in the C-terminal half of the protein, creating chromosomal repeats potentially capable of mediating the inversion of large chromosomal DNA. Cj and Cc Tlps were both found in association with only 14 different genes, indicating a limited genomic context. In Cj these Tlps defined IEs that were for the most part found at a single chromosomal location and comprised of a conserved set of genes. Cc IEs were situated at very different chromosomal locations, had different structures than Cj IEs, and were occasionally incomplete, therefore not capable of inversion. Tlps may have a role in Campylobacter genome structure and dynamics as well as acting as chemoreceptors mediating chemotactic responses.

Ruminant and chicken: important sources of campylobacteriosis in France despite a variation of source attribution in 2009 and 2015

Pathogen source attribution studies are a useful tool for identifying reservoirs of human infection. Based on Multilocus Sequence Typing (MLST) data, such studies have identified chicken as a major source of C. jejuni human infection. The use of whole genome sequence-based typing methods offers potential to improve the precision of attribution beyond that which is possible from 7 MLST loci. Using published data and 156 novel C. jejuni genomes sequenced in this study, we performed probabilistic host source attribution of clinical C. jejuni isolates from France using three types of genotype data: comparative genomic fingerprints; MLST genes; 15 host segregating genes previously identified by whole genome sequencing. Consistent with previous studies, chicken was an important source of campylobacteriosis in France (31–63% of clinical isolates assigned). There was also evidence that ruminants are a source (22–55% of clinical isolates assigned), suggesting that further investigation of potential transmission routes from ruminants to human would be useful. Additionally, we found evidence of environmental and pet sources. However, the relative importance as sources varied according to the year of isolation and the genotyping technique used. Annual variations in attribution emphasize the dynamic nature of zoonotic transmission and the need to perform source attribution regularly.

Whole Genome Sequencing and Multiplex qPCR Methods to Identify Campylobacter jejuni Encoding cst-II or cst-III Sialyltransferase.

Campylobacter jejuni causes more than 2 million cases of gastroenteritis annually in the United States, and is also linked to the autoimmune sequelae Guillan–Barre syndrome (GBS). GBS often results in flaccid paralysis, as the myelin sheaths of nerve cells are degraded by the adaptive immune response. Certain strains of C. jejuni modify their lipooligosaccharide (LOS) with the addition of neuraminic acid, resulting in LOS moieties that are structurally similar to gangliosides present on nerve cells. This can trigger GBS in a susceptible host, as antibodies generated against C. jejuni can cross-react with gangliosides, leading to demyelination of nerves and a loss of signal transduction. The goal of this study was to develop a quantitative PCR (qPCR) method and use whole genome sequencing data to detect the Campylobacter sialyltransferase (cst) genes responsible for the addition of neuraminic acid to LOS. The qPCR method was used to screen a library of 89 C. jejuni field samples collected by the Food and Drug Administration Pacific Northwest Lab (PNL) as well as clinical isolates transferred to PNL. In silico analysis was used to screen 827 C. jejuni genomes in the FDA GenomeTrakr SRA database. The results indicate that a majority of C. jejuni strains could produce LOS with ganglioside mimicry, as 43.8% of PNL isolates and 46.9% of the GenomeTrakr isolates lacked the cst genes. The methods described in this study can be used by public health laboratories to rapidly determine whether a C. jejuni isolate has the potential to induce GBS. Based on these results, a majority of C. jejuni in the PNL collection and submitted to GenomeTrakr have the potential to produce LOS that mimics human gangliosides.

Virulence and Genomic Feature of Multidrug Resistant Campylobacter jejuni Isolated from Broiler Chicken.

The aim of this study was to reveal the molecular mechanism involved in multidrug resistance and virulence of Campylobacter jejuni isolated from broiler chickens. The virulence of six multidrug resistant C. jejuni was determined by in vitro and in vivo methods. The de novo whole genome sequencing technology and molecular biology methods were used to analyze the genomic features associated with the multidrug resistance and virulence of a selected isolate (C. jejuni 1655). The comparative genomic analyses revealed a large number of single nucleotide polymorphisms, deletions, rearrangements, and inversions in C. jejuni 1655 compared to reference C. jejuni genomes. The co-emergence of Thr-86-Ile mutation in gyrA gene, A2075G mutation in 23S rRNA gene, tetO, aphA and aadE genes and pTet plasmid in C. jejuni 1655 contributed its multidrug resistance to fluoroquinolones, macrolides, tetracycline, and aminoglycosides. The combination of multiple virulence genes may work together to confer the relative higher virulence in C. jejuni 1655. The co-existence of mobile gene elements (e.g., pTet) and CRISPR-Cas system in C. jejuni 1655 may play an important role in the gene transfer and immune defense. The present study provides basic information of phenotypic and genomic features of C. jejuni 1655, a strain recently isolated from a chicken displaying multidrug resistance and relatively high level of virulence.

Association of Campylobacter Jejuni ssp. Jejuni Chemotaxis Receptor Genes with Multilocus Sequence Types and Source of Isolation.

Campylobacter jejuni’s flagellar locomotion is controlled by eleven chemoreceptors. Assessment of the distribution of the relevant chemoreceptor genes in the C. jejuni genomes deposited in the National Center for Biotechnology Information (NCBI) database led to the identification of two previously unknown tlp genes and a tlp5 pseudogene. These two chemoreceptor genes share the same locus in the C. jejuni genome with tlp4 and tlp11, but the gene region encoding the periplasmic ligand binding domain differs significantly from other chemoreceptor genes. Hence, they were named tlp12 and tlp13.Consequently, it was of interest to study their distribution in C. jejuni subpopulations of different clonality, and their cooccurrence with the eleven previously reported chemoreceptor genes. Therefore, the presence of all tlp genes was detected by polymerase chain reaction (PCR) in 292 multilocus sequence typing (MLST)-typed C. jejuni isolates from different hosts.The findings show interesting trends: Tlp4, tlp11, tlp12, and tlp13 appeared to be mutually exclusive and cooccur in a minor subset of isolates. Tlp4 was found to be present in only 33.56% of all tested isolates and was significantly less often detected in turkey isolates. Tlp11 was tested positive in only 17.8% of the isolates, while tlp12 was detected in 29.5% of all isolates, and tlp13 was found to be present in 38.7%.

Genetic import and phenotype specific alleles associated with hyper-invasion in Campylobacter jejuni.

BackgroundCampylobacter jejuni is a major zoonotic pathogen, causing gastroenteritis in humans. Invasion is an important pathogenesis trait by which C. jejuni causes disease. Here we report the genomic analysis of 134 strains to identify traits unique to hyperinvasive isolates.MethodsA total of 134 C. jejuni genomes were used to create a phylogenetic tree to position the hyperinvasive strains. Comparative genomics lead to the identification of mosaic capsule regions. A pan genome approach led to the discovery of unique loci, or loci with unique alleles, to the hyperinvasive strains.ResultsPhylogenetic analysis showed that the hyper-invasive phenotype is a generalist trait. Despite the fact that hyperinvasive strains are only distantly related based on the whole genome phylogeny, they all possess genes within the capsule region with high identity to capsule genes from C. jejuni subsp. doylei and C. lari. In addition there were genes unique to the hyper-invasive strains with identity to non-C. jejuni genes, as well as allelic variants of mainly pathogenesis related genes already known in the other C. jejuni. In particular, the sequence of flagella genes, flgD-E and flgL were highly conserved amongst the hyper-invasive strains and divergent from sequences in other C. jejuni. A novel cytolethal distending toxin (cdt) operon was also identified as present in all hyper-invasive strains in addition to the classic cdt operon present in other C. jejuni.ConclusionsOverall, the hyper-invasive phenotype is strongly linked to the presence of orthologous genes from other Campylobacter species in their genomes, notably within the capsule region, in addition to the observed association with unique allelic variants in flagellar genes and the secondary cdt operon which is unlikely under random sharing of accessory alleles in separate lineages.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-015-2087-y) contains supplementary material, which is available to authorized users.

Differential Distribution of Type II CRISPR-Cas Systems in Agricultural and Nonagricultural Campylobacter coli and Campylobacter jejuni Isolates Correlates with Lack of Shared Environments.

CRISPR (clustered regularly interspaced palindromic repeats)-Cas (CRISPR-associated) systems are sequence-specific adaptive defenses against phages and plasmids which are widespread in prokaryotes. Here we have studied whether phylogenetic relatedness or sharing of environmental niches affects the distribution and dissemination of Type II CRISPR-Cas systems, first in 132 bacterial genomes from 15 phylogenetic classes, ranging from Proteobacteria to Actinobacteria. There was clustering of distinct Type II CRISPR-Cas systems in phylogenetically distinct genera with varying G+C%, which share environmental niches. The distribution of CRISPR-Cas within a genus was studied using a large collection of genome sequences of the closely related Campylobacter species Campylobacter jejuni (N = 3,746) and Campylobacter coli (N = 486). The Cas gene cas9 and CRISPR-repeat are almost universally present in C. jejuni genomes (98.0% positive) but relatively rare in C. coli genomes (9.6% positive). Campylobacter jejuni and agricultural C. coli isolates share the C. jejuni CRISPR-Cas system, which is closely related to, but distinct from the C. coli CRISPR-Cas system found in C. coli isolates from nonagricultural sources. Analysis of the genomic position of CRISPR-Cas insertion suggests that the C. jejuni-type CRISPR-Cas has been transferred to agricultural C. coli. Conversely, the absence of the C. coli-type CRISPR-Cas in agricultural C. coli isolates may be due to these isolates not sharing the same environmental niche, and may be affected by farm hygiene and biosecurity practices in the agricultural sector. Finally, many CRISPR spacer alleles were linked with specific multilocus sequence types, suggesting that these can assist molecular epidemiology applications for C. jejuni and C. coli.

Prevention of biofilm formation and removal of existing biofilms by extracellular DNases of Campylobacter jejuni.

The fastidious nature of the foodborne bacterial pathogen Campylobacter jejuni contrasts with its ability to survive in the food chain. The formation of biofilms, or the integration into existing biofilms by C. jejuni, is thought to contribute to food chain survival. As extracellular DNA (eDNA) has previously been proposed to play a role in C. jejuni biofilms, we have investigated the role of extracellular DNases (eDNases) produced by C. jejuni in biofilm formation. A search of 2791 C. jejuni genomes highlighted that almost half of C. jejuni genomes contains at least one eDNase gene, but only a minority of isolates contains two or three of these eDNase genes, such as C. jejuni strain RM1221 which contains the cje0256, cje0566 and cje1441 eDNase genes. Strain RM1221 did not form biofilms, whereas the eDNase-negative strains NCTC 11168 and 81116 did. Incubation of pre-formed biofilms of NCTC 11168 with live C. jejuni RM1221 or with spent medium from a RM1221 culture resulted in removal of the biofilm. Inactivation of the cje1441 eDNase gene in strain RM1221 restored biofilm formation, and made the mutant unable to degrade biofilms of strain NCTC 11168. Finally, C. jejuni strain RM1221 was able to degrade genomic DNA from C. jejuni NCTC 11168, 81116 and RM1221, whereas strain NCTC 11168 and the RM1221 cje1441 mutant were unable to do so. This was mirrored by an absence of eDNA in overnight cultures of C. jejuni RM1221. This suggests that the activity of eDNases in C. jejuni affects biofilm formation and is not conducive to a biofilm lifestyle. These eDNases do however have a potential role in controlling biofilm formation by C. jejuni strains in food chain relevant environments.

Prevalence of Type VI Secretion System in Spanish Campylobacter jejuni Isolates.

Infections from Campylobacter jejuni pose a serious public health problem and are now considered the leading cause of foodborne bacterial gastroenteritis throughout the world. Sequencing of C. jejuni genomes has previously allowed a number of loci to be identified, which encode virulence factors that aid survival and pathogenicity. Recently, a Type VI secretion system (T6SS) consisting of 13 conserved genes was described in C. jejuni strains and recognised to promote pathogenicity and adaptation to the environment. In this study, we determined the presence of this T6SS in 63 Spanish C. jejuni isolates from the food chain and urban effluents using whole-genome sequencing. Our findings demonstrated that nine (14%) strains harboured the 13 ORFs found in prototype strain C. jejuni 108. Further studies will be necessary to determine the prevalence and importance of T6SS-positive C. jejuni strains.

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.