Genomic and Phenotypic Characterization of Campylobacter fetus subsp. venerealis Strains

Campylobacter fetus isolates from both human patients and healthy cattle carry three distinct cytolethal distending toxin (cdt) gene clusters

Campylobacter fetus is a zoonotic pathogen. Although the precise virulence mechanisms have not yet been fully elucidated, cytolethal distending toxin (CDT) is considered as one of the well-characterized virulence factors in Campylobacter. In silico analysis of the genome of C. fetus type strain ATCC27374T indicates that there are three cdt gene clusters, Cfcdt-I, Cfcdt-II and Cfcdt-III. However, it is not clear whether these clusters are ubiquitously present in C. fetus and their association with diseases in humans and animals. In this study, we have analyzed the distribution and nucleotide sequences of these cdt gene clusters in 137 C. fetus strains isolated from human patients and healthy cattle. MLST and PFGE were also applied to determine clonal relationship between C. fetus strains isolated from patients and cattle. We found all C. fetus strains carry three Cfcdt gene clusters by colony hybridization assay and the strains belonged to 38 different pulsotypes. Whole genome sequencing of 38 C. fetus strains was carried out to determine the entire cdt gene cluster sequences and their sequence type (ST). Among 38 strains, six STs were identified, and each cdt gene cluster showed high similarity (>99%). Interestingly, some of these Cfcdt genes are more similar to the cdt genes of other Campylobacter species than other Cfcdt gene types. Altogether, the results suggest that three Cfcdt gene clusters are highly conserved in C. fetus and the strains belonging to ST-6 may be more pathogenic to human.

Whole-genome comparison using complete genomes from Campylobacter fetus strains revealed single nucleotide polymorphisms on non-genomic islands for subspecies differentiation

Introduction

Bovine Genital Campylobacteriosis (BGC), caused by Campylobacter fetus subsp. venerealis, is a sexually transmitted bacterium that significantly impacts cattle reproductive performance. However, current detection methods lack consistency and reliability due to the close genetic similarity between C. fetus subsp. venerealis and C. fetus subsp. fetus. Therefore, this study aimed to utilize complete genome analysis to distinguish genetic features between C. fetus subsp. venerealis and other subspecies, thereby enhancing BGC detection for routine screening and epidemiological studies.

Methods and results

This study reported the complete genomes of four C. fetus subsp. fetus and five C. fetus subsp. venerealis, sequenced using long-read sequencing technologies. Comparative whole-genome analyses (n = 25) were conducted, incorporating an additional 16 complete C. fetus genomes from the NCBI database, to investigate the genomic differences between these two closely related C. fetus subspecies. Pan-genomic analyses revealed a core genome consisting of 1,561 genes and an accessory pangenome of 1,064 genes between the two C. fetus subspecies. However, no unique predicted genes were identified in either subspecies. Nonetheless, whole-genome single nucleotide polymorphisms (SNPs) analysis identified 289 SNPs unique to one or the C. fetus subspecies. After the removal of SNPs located on putative genomic islands, recombination sites, and those causing synonymous amino acid changes, the remaining 184 SNPs were functionally annotated. Candidate SNPs that were annotated with the KEGG "Peptidoglycan Biosynthesis" pathway were recruited for further analysis due to their potential association with the glycine intolerance characteristic of C. fetus subsp. venerealis and its biovar variant. Verification with 58 annotated C. fetus genomes, both complete and incomplete, from RefSeq, successfully classified these seven SNPs into two groups, aligning with their phenotypic identification as CFF (Campylobacter fetus subsp. fetus) or CFV/CFVi (Campylobacter fetus subsp. venerealis and its biovar variant). Furthermore, we demonstrated the application of mraY SNPs for detecting C. fetus subspecies using a quantitative PCR assay.

Discussion

Our results highlighted the high genetic stability of C. fetus subspecies. Nevertheless, Campylobacter fetus subsp. venerealis and its biovar variants encoded common SNPs in genes related to glycine intolerance, which differentiates them from C. fetus subsp. fetus. This discovery highlights the potential of employing a multiple-SNP assay for the precise differentiation of C. fetus subspecies.

Comparative Genomic Analysis of Campylobacter rectus and Closely Related Species

Campylobacter rectus is a gram-negative, anaerobic bacterium strongly associated with periodontitis. It also causes various extraoral infections and is linked to adverse pregnancy outcomes in humans and murine models. C. rectus and related oral Campylobacters have been termed "emerging Campylobacter species" because infections by these organisms are likely underreported. Previously, no comparative methods have been used to analyze more than single C. rectus strains and until recently, very few C. rectus genomes have been publicly available. More sequenced genomes and comparative analyses are needed to study the genomic features and pathogenicity of this species. We sequenced eight new C. rectus strains and used comparative methods to identify regions of interest. An emphasis was put on the type III flagellar secretion system (T3SS), type IV secretion system (T4SS), and type VI secretion system (T6SS) because these protein complexes are important for pathogenesis in other Campylobacter species. RAST, BV-BRC, and other bioinformatics tools were used to assemble, annotate, and compare these regions in the genomes. The pan-genome of C. rectus consists of 2670 genes with core and accessory genomes of 1429 and 1241 genes, respectively. All isolates analyzed in this study have T3SS and T6SS hallmark proteins, while five of the isolates are missing a T4SS system. Twenty-one prophage clusters were identified across the panel of isolates, including four that appear intact. Overall, significant genomic islands were found, suggesting regions in the genomes that underwent horizontal gene transfer. Additionally, the high frequency of CRISPR arrays and other repetitive elements has led to genome rearrangements across the strains, including in areas adjacent to secretion system gene clusters. This study describes the substantial diversity present among C. rectus isolates and highlights tools/assays that have been developed to permit functional genomic studies. Additionally, we have expanded the studies on C. showae T4SS since we have two new C. showae genomes to report. We also demonstrate that unlike C. rectus , C showae does not demonstrate evidence of intact T6SS except for the strain CAM. The only strain of sequenced C. massilensis has neither T4SS or T6SS.

Design of a multi-epitope-based vaccine candidate against Bovine Genital Campylobacteriosis using a reverse vaccinology approach

BACKGROUND

Bovine Genital Campylobacteriosis (BGC), a worldwide distributed venereal disease caused by Campylobacter fetus subsp. venerealis (Cfv), has a relevant negative economic impact in cattle herds. The control of BGC is hampered by the inexistence of globally available effective vaccines. The present in silico study aimed to develop a multi-epitope vaccine candidate against Cfv through reverse vaccinology.

RESULTS

The analysis of Cfv strain NCTC 10354 proteome allowed the identification of 9 proteins suitable for vaccine development. From these, an outer membrane protein, OmpA, and a flagellar protein, FliK, were selected for prediction of B-cell and T-cell epitopes. The top-ranked epitopes conservancy was assessed in 31 Cfv strains. The selected epitopes were integrated to form a multi-epitope fragment of 241 amino acids, which included 2 epitopes from OmpA and 13 epitopes from FliK linked by GPGPG linkers and connected to the cholera toxin subunit B by an EAAAK linker. The vaccine candidate was predicted to be antigenic, non-toxic, non-allergenic, and soluble upon overexpression. The protein structure was predicted and optimized, and the sequence was successfully cloned in silico into a plasmid vector. Additionally, immunological simulations demonstrated the vaccine candidate's ability to stimulate an immune response.

CONCLUSIONS

This study developed a novel vaccine candidate suitable for further in vitro and in vivo experimental validation, which may become a useful tool for the control of BGC.

Biochemical and molecular characterization of Campylobacter fetus isolates from bulls subjected to bovine genital campylobacteriosis diagnosis in Spain

BACKGROUND

Bovine genital campylobacteriosis (BGC) is caused by Campylobacter fetus subsp. venerealis (Cfv) including its biovar intermedius (Cfvi). This sexually transmitted disease induces early reproductive failure causing considerable economic losses in the cattle industry. Using a collection of well-characterized isolates (n = 13), C. fetus field isolates (n = 64) and saprophytic isolates resembling Campylobacter (n = 75) obtained from smegma samples of breeding bulls, this study evaluated the concordance of the most used phenotypic (H2S production in cysteine medium and 1% glycine tolerance) and molecular (PCR) methods for the diagnosis of BGC and assessed possible cross-reactions in the molecular diagnostic methods.

RESULTS

Characterization at the subspecies level (fetus vs. venerealis) of C. fetus isolated from bull preputial samples using phenotypic and molecular (PCR targeting nahE and ISCfe1) methods showed moderate concordance (κ = 0.462; CI: 0.256-0.669). No cross-reactions were observed with other saprophytic microaerophilic species or with other Campylobacter species that can be present in preputial samples. Whole genome sequencing (WGS) of discrepant isolates showed 100% agreement with PCR identification. For the differentiation of Cfv biovars, comparison of the H2S test (at 72 h and 5 days of incubation) and a PCR targeting the L-cysteine transporter genes showed higher concordance when H2S production was assessed after 5 days (72 h; κ = 0.553, 0.329-0.778 CI vs. 5 days; κ = 0.881, 0.631-1 CI), evidencing the efficacy of a longer incubation time.

CONCLUSIONS

This study confirmed the limitations of biochemical tests to correctly identify C. fetus subspecies and biovars. However, in the case of biovars, when extended incubation times for the H2S test (5 days) were used, phenotypic identification results were significantly improved, although PCR-based methods produced more accurate results. Perfect agreement of WGS with the PCR results and absence of cross-reactions with non-C. fetus saprophytic bacteria from the smegma demonstrated the usefulness of these methods. Nevertheless, the identification of new C. fetus subspecies-specific genes would help to improve BGC diagnosis.

Comparative pangenomic analysis of Campylobacter fetus isolated from Spanish bulls and other mammalian species

Campylobacter fetus comprises two closely related mammal-associated subspecies: Campylobacter fetus subsp. fetus (Cff) and Campylobacter fetus subsp. venerealis (Cfv). The latter causes bovine genital campylobacteriosis, a sexually-transmitted disease endemic in Spain that results in significant economic losses in the cattle industry. Here, 33 C. fetus Spanish isolates were whole-genome sequenced and compared with 62 publicly available C. fetus genomes from other countries. Genome-based taxonomic identification revealed high concordance with in silico PCR, confirming Spanish isolates as Cff (n = 4), Cfv (n = 9) and Cfv biovar intermedius (Cfvi, n = 20). MLST analysis assigned the Spanish isolates to 6 STs, including three novel: ST-76 and ST-77 for Cfv and ST-78 for Cff. Core genome SNP phylogenetic analysis of the 95 genomes identified multiple clusters, revealing associations at subspecies and biovar level between genomes with the same ST and separating the Cfvi genomes from Spain and other countries. A genome-wide association study identified pqqL as a Cfv-specific gene and a potential candidate for more accurate identification methods. Functionality analysis revealed variations in the accessory genome of C. fetus subspecies and biovars that deserve further studies. These results provide valuable information about the regional variants of C. fetus present in Spain and the genetic diversity and predicted functionality of the different subspecies.

Bovine campylobacteriosis in heifer: pathogenesis study and insights in the conventional and molecular diagnosis in an experimental bovine model and field cases

Campylobacter fetus spp. is a bacterium associated to reproductive losses in cattle worldwide. It is a venereal infectious disease known as bovine campilobacteriosis, with high impact mainly in countries with extensive production systems. Here, we show pathogenesis and diagnostic methods for Campylobacter fetus detection in cervico-vaginal mucus (CVM) samples from heifers experimentally infected and field cases from herds with low reproductive performance by campylobacteriosis infection. Bacterial culture, direct immunofluorescence test and qPCR were used as diagnostic methods to evaluate detection of C. fetus. In the experimental model 30 Aberdeen Angus and crossbred heifers and 4 Aberdeen Angus bulls for natural mating were assigned to 3 groups experimentally challenged with C. fetus subsp. fetus (Cff), C. fetus subsps venerealis (Cfv) and C. fetus subsp venerealis biovar intermedius (Cfvi), respectively, and a negative control group, all followed for 9 months. Also, field samples of CVM and aborted fetuses were recollected from seven beef cattle farms. Bacteriological culture had the higher C. fetus detection rate in CVM being the most appropriate, followed by qPCR (with commercial extraction DNA kit), direct immunofluorescence test and qPCR (with in-house extraction DNA method), in both, experimental model and field cases. From experimental model after natural mating, 62.5% and 25% heifers got pregnant from Cff and Cfvi groups, respectively, while from Cfv no pregnancy was detected. The strain more frequently detected was Cfvi, followed by Cff and Cfv. Colonization of Cff in female genital tract with high number of carriers and presence in aborted fetuses was evidenced, suggesting a high risk to bovine reproductive health. Bacteriemia was not detected after genital infection. Given the low detection rate of either test, we suggest the use of both, PCR based methods and bacterial culture could result in higher detection rate in farms with endemic campylobacteriosis.

Whole genome characterization of thermophilic Campylobacter species isolated from dairy manure in small specialty crop farms of Northeast Ohio

Introduction

With more public interest in consuming locally grown produce, small specialty crop farms (SSCF) are a viable and growing segment of the food production chain in the United States.

Methods

The goal of this study was to investigate the genomic diversity of Campylobacter isolated from dairy manure (n = 69) collected from 10 SSCF in Northeast Ohio between 2018 and 2020.

Results

A total of 56 C. jejuni and 13 C. coli isolates were sequenced. Multi-locus sequence typing (MLST) identified 22 sequence types (STs), with ST-922 (18%) and ST-61 (13%) predominant in C. jejuni and ST-829 (62%) and ST-1068 (38%) predominant in C. coli. Interestingly, isolates with similar genomic and gene contents were detected within and between SSCF over time, suggesting that Campylobacter could be transmitted between farms and may persist in a given SSCF over time. Virulence-associated genes (n = 35) involved in the uptake and utilization of potassium and organic compounds (succinate, gluconate, oxoglutarate, and malate) were detected only in the C. jejuni isolates, while 45 genes associated with increased resistance to environmental stresses (capsule production, cell envelope integrity, and iron uptake) were detected only in the C. coli isolates. Campylobacter coli isolates were also sub-divided into two distinct clusters based on the presence of unique prophages (n = 21) or IncQ conjugative plasmid/type-IV secretion system genes (n = 15). Campylobacter coli isolates harbored genes associated with resistance to streptomycin (aadE-Cc; 54%) and quinolone (gyrA-T86I; 77%), while C. jejuni had resistance genes for kanamycin (aph3'-IIIa; 20%). Both species harbored resistance genes associated with β-lactam (especially, blaOXA-193; up to 100%) and tetracycline (tetO; up to 59%).

Discussion/Conclusion

Our study demonstrated that Campylobacter genome plasticity associated with conjugative transfer might provide resistance to certain antimicrobials and viral infections via the acquisition of protein-encoding genes involved in mechanisms such as ribosomal protection and capsule modification.

Dysgonomonas mossii Strain Shenzhen WH 0221, a New Member of the Genus Dysgonomonas Isolated from the Blood of a Patient with Diabetic Nephropathy, Exhibits Multiple Antibiotic Resistance

Herein, we present a new bacterial strain isolated from infected blood of a patient with diabetic nephropathy. Matrix-assisted laser desorption/ionization time-of-flight mass spectrometry failed to identify the strain. 16S rRNA gene sequencing showed the highest similarity (>99.5%) with genus Dysgonomonas, but the strain could not be distinguished from Dysgonomonas oryzarvi and Dysgonomonas mossii. Whole genome sequencing, followed by phylogenetic analysis and average nucleotide identity (>95%) analysis, confirmed that the new strain represented Dysgonomonas mossii, leading it to be named Dysgonomonas mossii strain Shenzhen WH 0221. Shenzhen WH 0221 was 3.60 Mb with 37.4% GC content. It was Gram-stain negative, facultatively anaerobic, and grown on Columbia agar supplemented with 5% of sheep blood, exhibiting a smooth surface and pinpoint morphology. The morphological characteristics of this strain include a short rod shape without flagella and a size of 0.45-0.55 × 0.95-1.52 μm observed under transmission electron microscopy. The physiological and biochemical features and major cellular fatty acids (characterized by C_14:0 3-OH, C_14:0 9-CH₃, and C_16:0) differed from D. mossii CCUG 43457^T and other members of the genus Dysgonomonas. The isolate was found resistant to most cephalosporins, penicillin, norfloxacin, vancomycin, and chloramphenicol, but was susceptible to meropenem, imipenem, tetracycline, clindamycin, and amoxicillin-clavulanic acid. Genes kdpE, ykkD, cmeB, TLA-3, and vanRM found in its genome are probably associated with multiple antibiotic resistance. Lipopolysaccharides, capsules, and cytolysin may also help to illuminate its potential pathogenicity. This is the first report of a case of sepsis caused by Dysgonomonas mossii, and its pathogenic system was analyzed by whole genome sequencing. IMPORTANCE This study identified a new strain, Dysgonomonas mossii strain Shenzhen WH 0221, which has been first reported to cause sepsis isolated from infected blood of a patient with diabetic nephropathy. Physiological and biochemical characterizations, as well as overall fatty acid profile, distinguish Shenzhen WH 0221 from other species of the same genus. However, limited antibiotics were researched for Dysgonomonas mossii. Seventeen antibiotics spanning at least 6 classes were studied, providing a valuable guide to the clinical usage of drugs to treat Dysgonomonas mossii infection. For the first time, we report genome-based functional predictions for Dysgonomonas mossii. Five antibiotic resistance ontologies and more than 200 virulence factors likely underlie the multidrug resistance of Shenzhen WH 0221 and its potential pathogenicity.

Evaluation of PCR assays for Campylobacter fetus detection and discrimination between C. fetus subspecies in bovine preputial wash samples

Campylobacter fetus is a zoonotic pathogen found in cattle, in which it is one of the main causes of infectious infertility. Most diagnostic laboratories use PCR as quick easy tool for C. fetus identification. However, there is no standardized PCR assay for C. fetus detection and subspecies differentiation, hindering the comparison of results. In this study, we evaluated selected PCR assays targeting the 16S rRNA, gyrB, cpn60, cstA, cdtB and nahE genes for C. fetus identification and ISCfe1, sapB2, parA and virB11 for subspecies differentiation. Analytical sensitivity and specificity were assessed for each PCR assay, and the assays were then tested on 289 bull preputial samples that had also been analysed by 16S rRNA barcode metagenomics. In total, 41 C. fetus-positive samples were included. The P12 PCR assay targeting the gyrB gene performed best, detecting the pathogen in 95.1% of positive samples. For the discrimination of C. fetus subspecies, we were able to identify a proportion (85.4%) of the C. fetus-positive samples correctly as C. fetus venerealis with at least one subspecies-specific PCR, but C. fetus fetus was not detected in any of the samples tested. Remarkably, C. fetus subspecies amplification was observed following PCR on some samples (33.1%) considered C. fetus-negative, highlighting the need for rigorous criteria for discriminating between C. fetus subspecies, to improve understanding of the role of the two C. fetus subspecies in the epidemiology and pathogenesis of bovine infectious infertility.