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

Uncharacterized and lineage-specific accessory genes within the Proteus mirabilis pan-genome landscape

Proteus mirabilis is a Gram-negative bacterium recognized for its unique swarming motility and urease activity. A previous proteomic report on four strains hypothesized that, unlike other Gram-negative bacteria, P. mirabilis may not exhibit significant intraspecies variation in gene content. However, there has not been a comprehensive analysis of large numbers of P. mirabilis genomes from various sources to support or refute this hypothesis. We performed comparative genomic analysis on 2,060 Proteus genomes. We sequenced the genomes of 893 isolates recovered from clinical specimens from three large US academic medical centers, combined with 1,006 genomes from NCBI Assembly and 161 genomes assembled from Illumina reads in the public domain. We used average nucleotide identity (ANI) to delineate species and subspecies, core genome phylogenetic analysis to identify clusters of highly related P. mirabilis genomes, and pan-genome annotation to identify genes of interest not present in the model P. mirabilis strain HI4320. Within our cohort, Proteus is composed of 10 named species and 5 uncharacterized genomospecies. P. mirabilis can be subdivided into three subspecies; subspecies 1 represented 96.7% (1,822/1,883) of all genomes. The P. mirabilis pan-genome includes 15,399 genes outside of HI4320, and 34.3% (5,282/15,399) of these genes have no putative assigned function. Subspecies 1 is composed of several highly related clonal groups. Prophages and gene clusters encoding putatively extracellular-facing proteins are associated with clonal groups. Uncharacterized genes not present in the model strain P. mirabilis HI4320 but with homology to known virulence-associated operons can be identified within the pan-genome. IMPORTANCE Gram-negative bacteria use a variety of extracellular facing factors to interact with eukaryotic hosts. Due to intraspecies genetic variability, these factors may not be present in the model strain for a given organism, potentially providing incomplete understanding of host-microbial interactions. In contrast to previous reports on P. mirabilis, but similar to other Gram-negative bacteria, P. mirabilis has a mosaic genome with a linkage between phylogenetic position and accessory genome content. P. mirabilis encodes a variety of genes that may impact host-microbe dynamics beyond what is represented in the model strain HI4320. The diverse, whole-genome characterized strain bank from this work can be used in conjunction with reverse genetic and infection models to better understand the impact of accessory genome content on bacterial physiology and pathogenesis of infection.

Portrait of a generalist bacterium: pathoadaptation, metabolic specialization and extreme environments shape diversity of Staphylococcus saprophyticus

Staphylococcus saprophyticus is a Gram-positive, coagulase-negative staphylococcus found in diverse environments including soil and freshwater, meat, and dairy foods. S. saprophyticus is also an important cause of urinary tract infections (UTIs) in humans, and mastitis in cattle. However, the genetic determinants of virulence have not yet been identified, and it remains unclear whether there are distinct sub-populations adapted to human and animal hosts. Using a diverse sample of S. saprophyticus isolates from food, animals, environmental sources, and human infections, we characterized the population structure and diversity of global populations of S. saprophyticus . We found that divergence of the two major clades of S. saprophyticus is likely facilitated by barriers to horizontal gene transfer (HGT) and differences in metabolism. Using genome-wide association study (GWAS) tools we identified the first Type VII secretion system (T7SS) described in S. saprophyticus and its association with bovine mastitis. Finally, we found that in general, strains of S. saprophyticus from different niches are genetically similar with the exception of built environments, which function as a 'sink' for S. saprophyticus populations. This work increases our understanding of the ecology of S. saprophyticus and of the genomics of bacterial generalists.

Data summary

Raw sequencing data for newly sequenced S. saprophyticus isolates have been deposited to the NCBI SRA under the project accession PRJNA928770. A list of all genomes used in this work and their associated metadata are available in the supplementary material. Custom scripts used in the comparative genomics and GWAS analyses are available here: https://github.com/myoungblom/sapro_genomics .

Impact statement

It is not known whether human and cattle diseases caused by S. saprophyticus represent spillover events from a generalist adapted to survive in a range of environments, or whether the capacity to cause disease represents a specific adaptation. Seasonal cycles of S. saprophyticus UTIs and molecular epidemiological evidence suggest that these infections may be environmentally-acquired rather than via transmission from person to person. Using comparative genomics and genome wide association study tools, we found that S. saprophyticus appears adapted to inhabit a wide range of environments (generalist), with isolates from animals, food, natural environments and human infections being closely related. Bacteria that routinely switch environments, particularly between humans and animals, are of particular concern when it comes to the spread of antibiotic resistance from farm environments into human populations. This work provides a framework for comparative genomic analyses of bacterial generalists and furthers our understanding of how bacterial populations move between humans, animals, and the environment.

Monitoring within-farm transmission dynamics of antimicrobial-resistant Campylobacter in dairy cattle using broth microdilution and long-read whole genome sequencing

Campylobacter jejuni and Campylobacter coli are important foodborne zoonotic pathogens and cause for concern due to the increasing trend in antimicrobial resistance. A long-run surveillance study was conducted in animals from different age groups in five dairy cattle farms to investigate the within-farm diversity and transmission dynamics of resistant Campylobacter throughout time. The resistance phenotype of the circulating isolates (170 C. jejuni and 37 C. coli) was determined by broth microdilution and a selection of 56 isolates were whole genome sequenced using the Oxford-Nanopore long-fragment sequencing technology resulting in completely resolved and circularized genomes (both chromosomes and plasmids). C. jejuni was isolated from all farms while C. coli was isolated from only two farms, but resistance rates were higher in C. coli than in C. jejuni and in calves than in adult animals. Some genotypes (e.g. ST-48, gyrA_T86I/tet(O)/blaOXA-61 in farm F1; ST-12000, aadE-Cc/tet(O)/blaOXA-489 in F4) persisted throughout the study while others were only sporadically detected. Acquisition of extracellular genes from other isolates and intracellular mutational events were identified as the processes that led to the emergence of the resistant genotypes that spread within the herds. Monitoring with Oxford Nanopore Technologies sequencing helped to decipher the complex molecular epidemiology underlying the within-farm dissemination of resistant Campylobacter.

Whole-genome sequencing reveals changes in genomic diversity and distinctive repertoires of T3SS and T6SS effector candidates in Chilean clinical Campylobacter strains

Campylobacter is the leading cause of bacterial gastroenteritis worldwide and an emerging and neglected pathogen in South America. This zoonotic pathogen colonizes the gastrointestinal tract of a wide range of mammals and birds, with poultry as the most important reservoir for human infections. Apart from its high morbidity rates, the emergence of resistant strains is of global concern. The aims of this work were to determine genetic diversity, presence of antimicrobial resistance determinants and virulence potential of Campylobacter spp. isolated from patients with acute gastrointestinal disease at 'Clinica Alemana', Santiago de Chile. The study considered the isolation of Campylobacter spp., from stool samples during a 20-month period (January 2020 to September 2021). We sequenced (NextSeq, Illumina) and performed an in-depth analysis of the genome sequences of 88 Campylobacter jejuni and 2 Campylobacter coli strains isolated from clinical samples in Chile. We identified a high genetic diversity among C. jejuni strains and the emergence of prevalent clonal complexes, which were not identified in our previous reports. While ~40% of strains harbored a mutation in the gyrA gene associated with fluoroquinolone resistance, no macrolide-resistance determinants were detected. Interestingly, gene clusters encoding virulence factors such as the T6SS or genes associated with long-term sequelae such as Guillain-Barré syndrome showed lineage-relatedness. In addition, our analysis revealed a high degree of variability regarding the presence of fT3SS and T6SS effector proteins in comparison to type strains 81-176, F38011, and NCTC 11168 and 488. Our study provides important insights into the molecular epidemiology of this emerging foodborne pathogen. In addition, the differences observed regarding the repertoire of fT3SS and T6SS effector proteins could have an impact on the pathogenic potential and transmissibility of these Latin American isolates, posing another challenge in characterizing the infection dynamics of this emergent and neglected bacterial pathogen.

Oncogenic potential of Campylobacter infection in the gastrointestinal tract: narrative review

BACKGROUND

Campylobacter jejuni is the leading cause of zoonotic gastroenteritis. The other emerging group of Campylobacters spp. are part of human oral commensal, represented by C. concisus (CC), which has been recently linked to non-oral conditions. Although long-term gastrointestinal (GI) complications from these two groups of Campylobacters have been previously reviewed individually, overall impact of Campylobacter infection on GI carcinogenesis and their inflammatory precursor lesions has not been assessed collectively.

AIMS

To evaluate the available evidence concerning the association between Campylobacter infection/colonization and inflammatory bowel disease (IBD), reflux esophagitis/metaplasia colorectal cancer (CRC) and esophageal cancer (EC).

METHODS

We performed a comprehensive literature search of PubMed for relevant original publications and systematic reviews/meta-analyses of epidemiological and clinical studies. In addition, we gathered additional information concerning microbiological data, animal models and mechanistic data from in vitro studies.

RESULTS

Both retrospective and prospective studies on IBD showed relatively consistent increased risk associated with Campylobacter infection. Despite lack of supporting prospective studies, retrospective studies based on tissue/fecal microbiome revealed consistent enrichment of Campylobacter in CRC samples. Studies on EC precursor lesions (esophagitis and metaplasia) were generally supportive for the association with Campylobacter, while inconsistent observations on EC. Studies on both IBD and EC precursors suggested the predominant role of CC, but studies on CRC were not informative of species.

CONCLUSIONS

There is sufficient evidence calling for concerted effort in unveiling direct and indirect connection of this organism to colorectal and esophageal cancer in humans.

Genome-wide association reveals host-specific genomic traits in Escherichia coli

BACKGROUND

Escherichia coli is an opportunistic pathogen which colonizes various host species. However, to what extent genetic lineages of E. coli are adapted or restricted to specific hosts and the genomic determinants of such adaptation or restriction is poorly understood.

RESULTS

We randomly sampled E. coli isolates from four countries (Germany, UK, Spain, and Vietnam), obtained from five host species (human, pig, cattle, chicken, and wild boar) over 16 years, from both healthy and diseased hosts, to construct a collection of 1198 whole-genome sequenced E. coli isolates. We identified associations between specific E. coli lineages and the host from which they were isolated. A genome-wide association study (GWAS) identified several E. coli genes that were associated with human, cattle, or chicken hosts, whereas no genes associated with the pig host could be found. In silico characterization of nine contiguous genes (collectively designated as nan-9) associated with the human host indicated that these genes are involved in the metabolism of sialic acids (Sia). In contrast, the previously described sialic acid regulon known as sialoregulon (i.e. nanRATEK-yhcH, nanXY, and nanCMS) was not associated with any host species. In vitro growth experiments with a Δnan-9 E. coli mutant strain, using the sialic acids 5-N-acetylneuraminic acid (Neu5Ac) and N-glycolylneuraminic acid (Neu5Gc) as sole carbon source, showed impaired growth behaviour compared to the wild-type.

CONCLUSIONS

This study provides an extensive analysis of genetic determinants which may contribute to host specificity in E. coli. Our findings should inform risk analysis and epidemiological monitoring of (antimicrobial resistant) E. coli.

Whole-genome-based characterization of Campylobacter jejuni from human patients with gastroenteritis collected over an 18 year period reveals increasing prevalence of antimicrobial resistance

Campylobacteriosis is the most common cause of acute gastrointestinal bacterial infection in Europe, with most infections linked to the consumption of contaminated food. While previous studies found an increasing rate of antimicrobial resistance (AMR) in Campylobacter spp. over the past decades, the investigation of additional clinical isolates is likely to provide novel insights into the population structure and mechanisms of virulence and drug resistance of this important human pathogen. Therefore, we combined whole-genome sequencing and antimicrobial-susceptibility testing of 340 randomly selected Campylobacter jejuni isolates from humans with gastroenteritis, collected in Switzerland over an 18 year period. In our collection, the most common multilocus sequence types (STs) were ST-257 (n=44), ST-21 (n=36) and ST-50 (n=35); the most common clonal complexes (CCs) were CC-21 (n=102), CC-257 (n=49) and CC-48 (n=33). High heterogeneity was observed among STs, with the most abundant STs recurring over the entire study period, while others were observed only sporadically. Source attribution based on ST assigned more than half of the strains to the 'generalist' category (n=188), 25  % as 'poultry specialist' (n=83), and only a few to 'ruminant specialist' (n=11) or 'wild bird' origin (n=9). The isolates displayed an increased frequency of AMR from 2003 to 2020, with the highest rates of resistance observed for ciprofloxacin and nalidixic acid (49.8 %), followed by tetracycline (36.9 %). Quinolone-resistant isolates carried chromosomal gyrA mutations T86I (99.4 %) and T86A (0.6 %), whereas tetracycline-resistant isolates carried tet(O) (79.8 %) or mosaic tetO/32/O (20.2 %) genes. A novel chromosomal cassette carrying several resistance genes, including aph(3')-III, satA and aad(6), and flanked by insertion sequence elements was detected in one isolate. Collectively, our data revealed an increasing prevalence of resistance to quinolones and tetracycline in C. jejuni isolates from Swiss patients over time, linked to clonal expansion of gyrA mutants and acquisition of the tet(O) gene. Investigation of source attribution suggests that infections are most likely related to isolates from poultry or generalist backgrounds. These findings are relevant to guide future infection prevention and control strategies.

In silico prediction and expression analysis of vaccine candidate genes of Campylobacter jejuni

Campylobacter jejuni (C. jejuni) is the most common food-borne pathogen that causes human gastroenteritis in the United States. Consumption of contaminated poultry products is considered as the major source of human Campylobacter infection. An effective vaccine would be a promising alternative to antibiotic supplements to curb C. jejuni colonization in poultry gastrointestinal (GI) tract. However, the genetic diversity among the C. jejuni isolates makes vaccine production more challenging. Despite many attempts, an effective Campylobacter vaccine is not yet available. This study aimed to identify suitable candidates to develop a subunit vaccine against C. jejuni, which could reduce colonization in the GI tract of the poultry. In the current study, 4 C. jejuni strains were isolated from retail chicken meat and poultry litter samples and their genomes were sequenced utilizing next-generation sequencing technology. The genomic sequences of C. jejuni strains were screened to identify potential antigens utilizing the reverse vaccinology approach. In silico genome analysis predicted 3 conserved potential vaccine candidates (phospholipase A [PldA], TonB dependent vitamin B12 transporter [BtuB], and cytolethal distending toxin subunit B [CdtB]) suitable for the development of a vaccine. Furthermore, the expression of predicted genes during host-pathogen interaction was analyzed by an infection study using an avian macrophage-like immortalized cell line (HD11). The HD11 was infected with C. jejuni strains, and the RT-qPCR assay was performed to determine the expression of the predicted genes. The expression difference was analyzed using ΔΔCt methods. The results indicate that all 3 predicted genes, PldA, BtuB, and CdtB, were upregulated in 4 tested C. jejuni strains irrespective of their sources of isolation. In conclusion, in silico prediction and gene expression analysis during host-pathogen interactions identified 3 potential vaccine candidates for C. jejuni.

Phylogenetic Association and Genetic Factors in Cold Stress Tolerance in Campylobacter jejuni

Campylobacter jejuni is a major foodborne pathogen transmitted to humans primarily via contaminated poultry meat. Since poultry meat is generally processed, distributed, and stored in the cold chain, the survival of C. jejuni at refrigeration temperatures crucially affects human exposure to C. jejuni. Here, we investigated genetic factors associated with cold stress tolerance in C. jejuni. Seventy-nine C. jejuni strains isolated from retail raw chicken exhibited different survival levels at 4°C for 21 days. Multilocus sequence typing (MLST) clonal complex 21 (CC-21) and CC-443 were dominant among cold stress-tolerant strains, whereas CC-45 was common among cold stress-sensitive strains. Genome-wide average nucleotide identity (ANI) analysis identified a phylogenetic cluster associated with cold stress tolerance. Moreover, a pangenome analysis revealed 58 genes distinctively present in the cold stress-tolerant phylogenetic cluster. Among these 58 genes, cfrA, encoding the ferric enterobactin receptor involved in ion transport and metabolism, was selected for further analysis. Remarkably, the viability of a ΔcfrA mutant at 4°C was significantly decreased, while the levels of total reactive oxygen species and intracellular iron exceeded those of the wild type. Additionally, a knockout mutation of cfrA also significantly decreased the viability of three cold stress-tolerant isolates at 4°C, confirming the role of cfrA in cold stress tolerance. The results of this study demonstrate that unique phylogenetic clusters of C. jejuni associated with cold stress tolerance exist and that cfrA is a genetic factor contributing to cold stress tolerance in C. jejuni. IMPORTANCE The tolerance of foodborne pathogens to environmental stresses significantly affects food safety. Several studies have demonstrated that C. jejuni survives extended exposures to low temperatures, but the mechanisms of cold stress tolerance are not fully understood. Here, we demonstrate that C. jejuni strains in certain phylogenetic groups exhibit increased tolerance to cold stress. Notably, cfrA is present in the phylogenetic cluster associated with cold stress tolerance and plays a role in the survival of C. jejuni at low temperatures by alleviating oxidative stress. This is the first study to discover phylogenetic associations involving cold stress tolerance and to identify genetic elements conferring cold stress tolerance to C. jejuni.

Local accessory gene sharing among Egyptian Campylobacter potentially promotes the spread of antimicrobial resistance

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

A Complex Competitive Exclusion Culture Reduces Campylobacter jejuni Colonization in Broiler Chickens at Slaughter Age In Vivo

Diminishing Campylobacter prevalence in poultry flocks has proven to be extremely challenging. To date, efficacious control measures to reduce Campylobacter prevalence are still missing. A potential approach to control Campylobacter in modern poultry productions is to occupy its niche in the mucosal layer by administering live intestinal microbiota from adult chickens to dayold-chicks (competitive exclusion (CE)). Therefore, this in vivo study investigates the efficacy of a complex CE culture to reduce Campylobacter (C.) jejuni colonization in broiler chickens. For this purpose, the complex CE culture was applied twice: once by spray application to day-old chicks immediately after hatching (on the 1st day of life) and subsequently by an additional application via drinking water on the 25th day of life. We observed a consistent and statistically significant reduction of C. jejuni counts in cloacal swabs throughout the entire fattening period. At the end of the trial after necropsy (at 33 days of age), C. jejuni cecal counts also showed a statistically significant decrease of 1 log₁₀ MPN/g compared to the control group. Likewise, colon counts were reduced by 2.0 log₁₀ MPN/g. These results suggest that CE cultures can be considered a practically relevant control strategy to reduce C. jejuni colonization in broiler chickens on poultry farms.

Genomic Analysis and Antimicrobial Resistance of Campylobacter jejuni and Campylobacter coli in Peru

Campylobacter is the leading cause of human bacterial gastroenteritis worldwide and has a major impact on global public health. Whole Genome Sequencing (WGS) is a powerful tool applied in the study of foodborne pathogens. The objective of the present study was to apply WGS to determine the genetic diversity, virulence factors and determinants of antimicrobial resistance of the populations of C. jejuni and C. coli in Peru. A total of 129 Campylobacter strains (108 C. jejuni and 21 C. coli) were sequenced using Illumina Miseq platform. In silico MLST analysis identified a high genetic diversity among those strains with 30 sequence types (STs), several of them within 11 clonal complexes (CC) for C. jejuni, while the strains of C. coli belonged to a single CC with 8 different STs. Phylogeny analysis showed that Peruvian C. jejuni strains were divided into 2 clades with 5 populations, while C. coli formed a single clade with 4 populations. Furthermore, in silico analyses showed the presence of several genes associated with adherence, colonization and invasion among both species: cadF (83.7%), jlpA (81.4%), racR (100%), dnaJ (83.7%), pebA (83.7%), pldA (82.1%), porA (84.5%), ceuE (82.9%), ciaB (78.3%), iamB (86.8%), and flaC (100%). The majority (82.9%) of the Campylobacter strains carried the cdtABC operon which code for cytolethal distending toxin (CDT). Half of them (50.4%) carried genes associated with the presence of T6SS, while the frequency of genes associated with T4SS were relatively low (11.6%). Genetic markers associated with resistance to quinolones, tetracycline (tetO, tetW/N/W), beta-lactamases (bla_oxa–61), macrolides (A2075G in 23S rRNA) were found in 94.5, 21.7, 66.7, 6.2, 69.8, and 18.6% of strains, respectively. The cmeABC multidrug efflux operon was present in 78.3% of strains. This study highlights the importance of using WGS in the surveillance of emerging pathogens associated with foodborne diseases, providing genomic information on genetic diversity, virulence mechanisms and determinants of antimicrobial resistance. The description of several Campylobacter genotypes having many virulence factors and resistance to quinolones and tetracyclines circulating in Peru provides important information which helps in the monitoring, control and prevention strategies of this emerging pathogen in our country.

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.