Campylobacter fetus subspecies specific PCR assays inferred from comparative genomic analysis for accurate subspecies identification

Bovine Genital Campylobacteriosis (BGC) is caused by Campylobacter fetus subsp. venerealis and is a notifiable disease to the WOAH (World Organisation for Animal Health). For an effective BGC control program, the reliable differentiation of Campylobacter fetus subsp. venerealis (Cfv) from the closely related Campylobacter fetus subsp. fetus (Cff) is required. However, the available molecular C. fetus subspecies identification assays lack sensitivity and specificity to differentiate C. fetus isolates based on their phenotypic or genotypic differences. Furthermore, the current biochemical subspecies identification is not fully congruent with the genomic differentiation of C. fetus strains. In this study, the genome sequences of 41C. fetus strains with well identified subspecies, were analyzed with the large-scale BLAST score ratio (LS-BSR) pipeline to identify Cff and Cfv specific sequences. With this analysis, the asd gene encoding an aspartate-semialdehyde dehydrogenase was identified, which contained a 6-bp Cff-specific sequence, and this 6-bp sequence was absent in the asd gene of Cfv strains. This sequence was used for the development of PCR assays to differentiate Cff and Cfv strains. The C. fetus subspecies identification of the developed asd PCR assays was in full congruence with the genomic classification of strains and are recommended for molecular identification of C. fetus subspecies in BGC control programs. The asd PCR can be assessed on sequenced genomes using a web interface containing the Cfvcatch tool, which includes placement of the tested genome in a phylogenetic tree with reference C. fetus genomes to distinguish the two subspecies and to detect antimicrobial resistance genes.