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

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

Genomic diversity and antimicrobial resistance of Campylobacter spp. from humans and livestock in Nigeria.

BACKGROUND: Campylobacter spp. are zoonotic pathogens, ubiquitous and are found naturally as commensals in livestock from where they can be transmitted to humans directly or through animal products. The genomic diversity and antimicrobial resistance profile of Campylobacter was investigated with a focus on C. jejuni and C. coli in humans and livestock (poultry and cattle) from Nigeria. METHODS: 586 human stool samples and 472 faecal samples from livestock were cultured for thermophilic Campylobacter species on modified charcoal cefoperazone deoxycholate agar (mCCDA). Culture in combination with whole genome sequencing identified and confirmed the presence of Campylobacter in humans and animals from the study area. Further analysis of the sequences was performed to determine multilocus sequence types and genetic determinants of antimicrobial resistance to fluoroquinolone, betalactam, tetracycline and macrolide classes of antimicrobials. RESULTS: From the human stool samples tested, 50 (9%) were positive of which 33 (66%) were C. jejuni, 14 (28%) were C. coli while 3 (6%) were C. hyointestinalis. In livestock, 132 (28%) were positive. Thirty one (7%) were C. jejuni while 101 (21%) were C. coli. Whole genome sequencing and MLST of the isolates revealed a total of 32 sequence types (STs) identified from 47 human isolates while 48 STs were identified in 124 isolates from livestock indicating a population which was overall, genetically diverse with a few more dominant strains. The antimicrobial resistance profiles of the isolates indicated a higher prevalence of resistance in Campylobacter isolated from livestock than in humans. Generally, resistance was greatest for betalactams (42%) closely followed by fluoroquinolones (41%), tetracyclines (15%) and lastly macrolides (2%). Multidrug resistance to three or more antimicrobials was observed in 24 (13%) isolates from humans (n = 1, 4%) and chicken (n = 23, 96%). CONCLUSIONS: This study has further contributed information about the epidemiology, genetic diversity and antimicrobial resistance profile of thermophilic Campylobacter in Nigeria.

Importance of untested infectious individuals for interventions to suppress COVID-19.

The impact of the extent of testing infectious individuals on suppression of COVID-19 is illustrated from the early stages of outbreaks in Germany, the Hubei province of China, Italy, Spain and the UK. The predicted percentage of untested infected individuals depends on the specific outbreak but we found that they typically represent 60-80% of all infected individuals during the early stages of the outbreaks. We propose that reducing the underlying transmission from untested cases is crucial to suppress the virus. This can be achieved through enhanced testing in combination with social distancing and other interventions that reduce transmission such as wearing face masks. Once transmission from silent carriers is kept under control by these means, the virus could have been fully suppressed through fast isolation and contact tracing of tested cases.

Estimated Dissemination Ratio-A Practical Alternative to the Reproduction Number for Infectious Diseases.

Policymakers require consistent and accessible tools to monitor the progress of an epidemic and the impact of control measures in real time. One such measure is the Estimated Dissemination Ratio (EDR), a straightforward, easily replicable, and robust measure of the trajectory of an outbreak that has been used for many years in the control of infectious disease in livestock. It is simple to calculate and explain. Its calculation and use are discussed below together with examples from the current COVID-19 outbreak in the UK. These applications illustrate that EDR can demonstrate changes in transmission rate before they may be clear from the epidemic curve. Thus, EDR can provide an early warning that an epidemic is resuming growth, allowing earlier intervention. A conceptual comparison between EDR and the commonly used reproduction number is also provided.

Mining whole genome sequence data to efficiently attribute individuals to source populations.

Whole genome sequence (WGS) data could transform our ability to attribute individuals to source populations. However, methods that efficiently mine these data are yet to be developed. We present a minimal multilocus distance (MMD) method which rapidly deals with these large data sets as well as methods for optimally selecting loci. This was applied on WGS data to determine the source of human campylobacteriosis, the geographical origin of diverse biological species including humans and proteomic data to classify breast cancer tumours. The MMD method provides a highly accurate attribution which is computationally efficient for extended genotypes. These methods are generic, easy to implement for WGS and proteomic data and have wide application.

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.

Variability in growth responses of non-O157 EHEC isolates in leafy vegetables, sprouted seeds and soil extracts occurs at the isolate level.

Foods of plant origin are recognised as a major source of foodborne pathogens, in particular for Shigatoxigenic Escherichia coli (STEC). Most work for STEC and plant-based fresh produce has focused on the most prevalent outbreak serogroup, O157. However, non-O157 STEC is an emerging hazard, and as such it is important to characterise aspects within this group that reflect their ability to colonise alternative hosts and habitats relevant to horticultural production. Growth kinetics were quantified for a diverse set of clinical enterohaemorrhagic E. coli isolates in extracts made from different tissues of spinach, lettuce or sprouted seeds, or from soil, to represent association with ready-to-eat fresh produce production. For leafy vegetables, spinach apoplast supported the fastest rates of growth and lettuce root extracts generated the slowest growth rates. Growth rates were similar for the majority of isolates in fenugreek or alfalfa sprouted seed extracts. Monosaccharides were the major driver of bacterial growth. No correlations were found for growth rates between different serotypes or for Shigatoxin gene carriage. Thus, growth rates varied in a plant-dependent and isolate-dependent manner, for all plant or soil extracts tested, indicative of isolate-specific differences in metabolic flexibility. These findings are relevant for risk assessment of non-O157 STEC.

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.

Nationwide Stepwise Emergence and Evolution of Multidrug-Resistant Campylobacter jejuni Sequence Type 5136, United Kingdom.

We examined whole-genome-sequenced Campylobacter jejuni and C. coli from 2012-2015 isolated from birds and human stool samples in North East Scotland for the presence of antimicrobial resistance genes. We found that sequence type (ST) 5136 (clonal complex 464) was the most prevalent multidrug-resistant strain of C. jejuni exclusively associated with poultry host reservoirs and recovered from human cases of campylobacteriosis. Tetracycline resistance in ST5136 isolates was due to a tet(O/32/O) mosaic gene, ampicillin resistance was conferred by G → T transversion in the -10 promoter region of blaOXA-193, fluoroquinolone resistance was due to C257T change in gyrA, and aminoglycoside resistance was conferred by aac. Whole-genome analysis showed that the strain ST5136 evolved from ST464. The nationwide emergence of ST5136 was probably due to stepwise acquisition of antimicrobial resistance genes selected by high use of ß-lactam, tetracycline, fluoroquinolone, and aminoglycoside classes of drugs in the poultry industry.

Challenges of biofilm control and utilization: lessons from mathematical modelling.

This article reviews modern applications of mathematical descriptions of biofilm formation. The focus is on theoretically obtained results which have implications for areas including the medical sector, food industry and wastewater treatment. Examples are given as to how models have contributed to the overall knowledge on biofilms and how they are used to predict biofilm behaviour. We conclude that the use of mathematical models of biofilms has demonstrated over the years the ability to significantly contribute to the vast field of biofilm research. Among other things, they have been used to test various hypotheses on the nature of interspecies interactions, viability of biofilm treatment methods or forces behind observed biofilm pattern formations. Mathematical models can also play a key role in future biofilm research. Many models nowadays are analysed through computer simulations and continue to improve along with computational capabilities. We predict that models will keep on providing answers to important challenges involving biofilm formation. However, further strengthening of the ties between various disciplines is necessary to fully use the tools of collective knowledge in tackling the biofilm phenomenon.

Influence of Plant Species, Tissue Type, and Temperature on the Capacity of Shiga-Toxigenic Escherichia coli To Colonize, Grow, and Be Internalized by Plants.

Contamination of fresh produce with pathogenic Escherichia coli, including Shiga-toxigenic E. coli (STEC), represents a serious risk to human health. Colonization is governed by multiple bacterial and plant factors that can impact the probability and suitability of bacterial growth. Thus, we aimed to determine whether the growth potential of STEC for plants associated with foodborne outbreaks (two leafy vegetables and two sprouted seed species) is predictive of the colonization of living plants, as assessed from growth kinetics and biofilm formation in plant extracts. The fitness of STEC isolates was compared to that of environmental E. coli isolates at temperatures relevant to plant growth. Growth kinetics in plant extracts varied in a plant-dependent and isolate-dependent manner for all isolates, with spinach leaf lysates supporting the highest rates of growth. Spinach extracts also supported the highest levels of biofilm formation. Saccharides were identified to be the major driver of bacterial growth, although no single metabolite could be correlated with growth kinetics. The highest level of in planta colonization occurred on alfalfa sprouts, though internalization was 10 times more prevalent in the leafy vegetables than in sprouted seeds. Marked differences in in planta growth meant that the growth potential of STEC could be inferred only for sprouted seeds. In contrast, biofilm formation in extracts related to spinach colonization. Overall, the capacity of E. coli to colonize, grow, and be internalized within plants or plant-derived matrices was influenced by the isolate type, plant species, plant tissue type, and temperature, complicating any straightforward relationship between in vitro and in planta behaviors.IMPORTANCE Fresh produce is an important vehicle for STEC transmission, and experimental evidence shows that STEC can colonize plants as secondary hosts, but differences in the capacity to colonize occur between different plant species and tissues. Therefore, an understanding of the impact that these plant factors have on the ability of STEC to grow and establish is required for food safety considerations and risk assessment. Here, we determined whether growth and the ability of STEC to form biofilms in plant extracts could be related to specific plant metabolites or could predict the ability of the bacteria to colonize living plants. Growth rates for sprouted seeds (alfalfa and fenugreek) but not those for leafy vegetables (lettuce and spinach) exhibited a positive relationship between plant extracts and living plants. Therefore, the detailed variations at the level of the bacterial isolate, plant species, and tissue type all need to be considered in risk assessment.

Climate, human behaviour or environment: individual-based modelling of Campylobacter seasonality and strategies to reduce disease burden.

BACKGROUND: With over 800 million cases globally, campylobacteriosis is a major cause of food borne disease. In temperate climates incidence is highly seasonal but the underlying mechanisms are poorly understood, making human disease control difficult. We hypothesised that observed disease patterns reflect complex interactions between weather, patterns of human risk behaviour, immune status and level of food contamination. Only by understanding these can we find effective interventions. METHODS: We analysed trends in human Campylobacter cases in NE England from 2004 to 2009, investigating the associations between different risk factors and disease using time-series models. We then developed an individual-based (IB) model of risk behaviour, human immunological responses to infection and environmental contamination driven by weather and land use. We parameterised the IB model for NE England and compared outputs to observed numbers of reported cases each month in the population in 2004-2009. Finally, we used it to investigate different community level disease reduction strategies. RESULTS: Risk behaviours like countryside visits (t = 3.665, P < 0.001 and t = - 2.187, P = 0.029 for temperature and rainfall respectively), and consumption of barbecued food were strongly associated with weather, (t = 3.219, P = 0.002 and t = 2.015, P = 0.045 for weekly average temperature and average maximum temperature respectively) and also rain (t = 2.254, P = 0.02527). This suggests that the effect of weather was indirect, acting through changes in risk behaviour. The seasonal pattern of cases predicted by the IB model was significantly related to observed patterns (r = 0.72, P < 0.001) indicating that simulating risk behaviour could produce the observed seasonal patterns of cases. A vaccination strategy providing short-term immunity was more effective than educational interventions to modify human risk behaviour. Extending immunity to 1 year from 20 days reduced disease burden by an order of magnitude (from 2412-2414 to 203-309 cases per 50,000 person-years). CONCLUSIONS: This is the first interdisciplinary study to integrate environment, risk behaviour, socio-demographics and immunology to model Campylobacter infection, including pathways to mitigation. We conclude that vaccination is likely to be the best route for intervening against campylobacteriosis despite the technical problems associated with understanding both the underlying human immunology and genetic variation in the pathogen, and the likely cost of vaccine development.

Phylogeographic Analysis Reveals Multiple International transmission Events Have Driven the Global Emergence of Escherichia coli O157:H7.

BACKGROUND: Shiga toxin-producing Escherchia coli (STEC) O157:H7 is a zoonotic pathogen that causes numerous food and waterborne disease outbreaks. It is globally distributed, but its origin and the temporal sequence of its geographical spread are unknown. METHODS: We analyzed whole-genome sequencing data of 757 isolates from 4 continents, and performed a pan-genome analysis to identify the core genome and, from this, extracted single-nucleotide polymorphisms. A timed phylogeographic analysis was performed on a subset of the isolates to investigate its worldwide spread. RESULTS: The common ancestor of this set of isolates occurred around 1890 (1845-1925) and originated from the Netherlands. Phylogeographic analysis identified 34 major transmission events. The earliest were predominantly intercontinental, moving from Europe to Australia around 1937 (1909-1958), to the United States in 1941 (1921-1962), to Canada in 1960 (1943-1979), and from Australia to New Zealand in 1966 (1943-1982). This pre-dates the first reported human case of E. coli O157:H7, which was in 1975 from the United States. CONCLUSIONS: Inter- and intra-continental transmission events have resulted in the current international distribution of E. coli O157:H7, and it is likely that these events were facilitated by animal movements (eg, Holstein Friesian cattle). These findings will inform policy on action that is crucial to reduce the further spread of E. coli O157:H7 and other (emerging) STEC strains globally.

Comparative Genomic Analysis Identifies a Campylobacter Clade Deficient in Selenium Metabolism.

The nonthermotolerant Campylobacter species C. fetus, C. hyointestinalis, C. iguaniorum, and C. lanienae form a distinct phylogenetic cluster within the genus. These species are primarily isolated from foraging (swine) or grazing (e.g., cattle, sheep) animals and cause sporadic and infrequent human illness. Previous typing studies identified three putative novel C. lanienae-related taxa, based on either MLST or atpA sequence data. To further characterize these putative novel taxa and the C. fetus group as a whole, 76 genomes were sequenced, either to completion or to draft level. These genomes represent 26 C. lanienae strains and 50 strains of the three novel taxa. C. fetus, C. hyointestinalis and C. iguaniorum genomes were previously sequenced to completion; therefore, a comparative genomic analysis across the entire C. fetus group was conducted (including average nucleotide identity analysis) that supports the initial identification of these three novel Campylobacter species. Furthermore, C. lanienae and the three putative novel species form a discrete clade within the C. fetus group, which we have termed the C. lanienae clade. This clade is distinguished from other members of the C. fetus group by a reduced genome size and distinct CRISPR/Cas systems. Moreover, there are two signature characteristics of the C. lanienae clade. C. lanienae clade genomes carry four to ten unlinked and similar, but nonidentical, flagellin genes. Additionally, all 76 C. lanienae clade genomes sequenced demonstrate a complete absence of genes related to selenium metabolism, including genes encoding the selenocysteine insertion machinery, selenoproteins, and the selenocysteinyl tRNA.

Genomic plasticity and rapid host switching can promote the evolution of generalism: a case study in the zoonotic pathogen Campylobacter.

Horizontal gene transfer accelerates bacterial adaptation to novel environments, allowing selection to act on genes that have evolved in multiple genetic backgrounds. This can lead to ecological specialization. However, little is known about how zoonotic bacteria maintain the ability to colonize multiple hosts whilst competing with specialists in the same niche. Here we develop a stochastic evolutionary model and show how genetic transfer of host segregating alleles, distributed as predicted for niche specifying genes, and the opportunity for host transition could interact to promote the emergence of host generalist lineages of the zoonotic bacterium Campylobacter. Using a modelling approach we show that increasing levels of homologous recombination enhance the efficiency with which selection can fix combinations of beneficial alleles, speeding adaptation. We then show how these predictions change in a multi-host system, with low levels of recombination, consistent with real r/m estimates, increasing the standing variation in the population, allowing a more effective response to changes in the selective landscape. Our analysis explains how observed gradients of host specialism and generalism can evolve in a multihost system through the transfer of ecologically important loci among coexisting strains.

Novel Sampling Method for Assessing Human-Pathogen Interactions in the Natural Environment Using Boot Socks and Citizen Scientists, with Application to Campylobacter Seasonality.

This paper introduces a novel method for sampling pathogens in natural environments. It uses fabric boot socks worn over walkers' shoes to allow the collection of composite samples over large areas. Wide-area sampling is better suited to studies focusing on human exposure to pathogens (e.g., recreational walking). This sampling method is implemented using a citizen science approach: groups of three walkers wearing boot socks undertook one of six routes, 40 times over 16 months in the North West (NW) and East Anglian (EA) regions of England. To validate this methodology, we report the successful implementation of this citizen science approach, the observation that Campylobacter bacteria were detected on 47% of boot socks, and the observation that multiple boot socks from individual walks produced consistent results. The findings indicate higher Campylobacter levels in the livestock-dominated NW than in EA (55.8% versus 38.6%). Seasonal differences in the presence of Campylobacter bacteria were found between the regions, with indications of winter peaks in both regions but a spring peak in the NW. The presence of Campylobacter bacteria on boot socks was negatively associated with ambient temperature (P = 0.011) and positively associated with precipitation (P < 0.001), results consistent with our understanding of Campylobacter survival and the probability of material adhering to boot socks. Campylobacter jejuni was the predominant species found; Campylobacter coli was largely restricted to the livestock-dominated NW. Source attribution analysis indicated that the potential source of C. jejuni was predominantly sheep in the NW and wild birds in EA but did not differ between peak and nonpeak periods of human incidence.IMPORTANCE There is debate in the literature on the pathways through which pathogens are transferred from the environment to humans. We report on the success of a novel method for sampling human-pathogen interactions using boot socks and citizen science techniques, which enable us to sample human-pathogen interactions that may occur through visits to natural environments. This contrasts with traditional environmental sampling, which is based on spot sampling techniques and does not sample human-pathogen interactions. Our methods are of practical value to scientists trying to understand the transmission of pathogens from the environment to people. Our findings provide insight into the risk of Campylobacter exposure from recreational visits and an understanding of seasonal differences in risk and the factors behind these patterns. We highlight the Campylobacter species predominantly encountered and the potential sources of C. jejuni.

Whole genome sequence analysis indicates recent diversification of mammal-associated Campylobacter fetus and implicates a genetic factor associated with H2S production.

BACKGROUND: Campylobacter fetus (C. fetus) can cause disease in both humans and animals. C. fetus has been divided into three subspecies: C. fetus subsp. fetus (Cff), C. fetus subsp. venerealis (Cfv) and C. fetus subsp. testudinum (Cft). Subspecies identification of mammal-associated C. fetus strains is crucial in the control of Bovine Genital Campylobacteriosis (BGC), a syndrome associated with Cfv. The prescribed methods for subspecies identification of the Cff and Cfv isolates are: tolerance to 1 % glycine and H2S production. RESULTS: In this study, we observed the deletion of a putative cysteine transporter in the Cfv strains, which are not able to produce H2S from L-cysteine. Phylogenetic reconstruction of the core genome single nucleotide polymorphisms (SNPs) within Cff and Cfv strains divided these strains into five different clades and showed that the Cfv clade and a Cff clade evolved from a single Cff ancestor. CONCLUSIONS: Multiple C. fetus clades were observed, which were not consistent with the biochemical differentiation of the strains. This suggests the need for a closer evaluation of the current C. fetus subspecies differentiation, considering that the phenotypic differentiation is still applied in BGC control programs.

Comparative Genomics of Campylobacter fetus from Reptiles and Mammals Reveals Divergent Evolution in Host-Associated Lineages.

Campylobacter fetus currently comprises three recognized subspecies, which display distinct host association. Campylobacter fetus subsp. fetus and C fetus subsp. venerealis are both associated with endothermic mammals, primarily ruminants, whereas C fetus subsp. testudinum is primarily associated with ectothermic reptiles. Both C. fetus subsp. testudinum and C. fetus subsp. fetus have been associated with severe infections, often with a systemic component, in immunocompromised humans. To study the genetic factors associated with the distinct host dichotomy in C. fetus, whole-genome sequencing and comparison of mammal- and reptile-associated C fetus was performed. The genomes of C fetus subsp. testudinum isolated from either reptiles or humans were compared with elucidate the genetic factors associated with pathogenicity in humans. Genomic comparisons showed conservation of gene content and organization among C fetus subspecies, but a clear distinction between mammal- and reptile-associated C fetus was observed. Several genomic regions appeared to be subspecies specific, including a putative tricarballylate catabolism pathway, exclusively present in C fetus subsp. testudinum strains. Within C fetus subsp. testudinum, sapA, sapB, and sapAB type strains were observed. The recombinant locus iamABC (mlaFED) was exclusively associated with invasive C fetus subsp. testudinum strains isolated from humans. A phylogenetic reconstruction was consistent with divergent evolution in host-associated strains and the existence of a barrier to lateral gene transfer between mammal- and reptile-associated C fetus Overall, this study shows that reptile-associated C fetus subsp. testudinum is genetically divergent from mammal-associated C fetus subspecies.

Campylobacter fetus Subspecies Contain Conserved Type IV Secretion Systems on Multiple Genomic Islands and Plasmids.

The features contributing to differences in pathogenicity of the Campylobacter fetus subspecies are unknown. Putative factors involved in pathogenesis are located in genomic islands that encode a type IV secretion system (T4SS) and fic domain (filamentation induced by cyclic AMP) proteins, which may disrupt host cell processes. In the genomes of 27 C. fetus strains, three phylogenetically-different T4SS-encoding regions (T4SSs) were identified: one was located in both the chromosome and in extra-chromosomal plasmids; one was located exclusively in the chromosome; and one exclusively in extra-chromosomal plasmids. We observed that C. fetus strains can contain multiple T4SSs and that homologous T4SSs can be present both in chromosomal genomic islands (GI) and on plasmids in the C. fetus strains. The GIs of the chromosomally located T4SS differed mainly by the presence of fic genes, insertion sequence elements and phage-related or hypothetical proteins. Comparative analysis showed that T4SS sequences, inserted in the same locations, were conserved in the studied C. fetus genomes. Using phylogenetic analysis of the T4SSs, it was shown that C. fetus may have acquired the T4SS regions from other Campylobacter species by horizontal gene transfer. The identified T4SSs and fic genes were found in Cff and Cfv strains, although the presence of T4SSs and fic genes were significantly associated with Cfv strains. The T4SSs and fic genes could not be associated with S-layer serotypes or geographical origin of the strains.