Pan-Genome Analysis of Campylobacter: Insights on the Genomic Diversity and Virulence Profile

Identification of Novel Gene-Specific Markers for Differentiating Various Pathogenic Campylobacter Species Using a Pangenome Analysis Approach

Campylobacter spp. are the causative agents of campylobacteriosis, a major foodborne illness globally, with millions of cases reported annually. These pathogens pose significant risks to both human and animal health. Conventional culture-based diagnostic methods are labor-intensive and time-consuming, underscoring the need for more efficient molecular detection strategies. This study employed a pangenomic analysis to identify novel gene-specific markers for pathogenic Campylobacter species and subspecies, laying the groundwork for their application in diverse diagnostic assays. A curated dataset of 105 high-quality genomes, representing 33 species and 9 subspecies, was analyzed using the Roary ILP Bacterial Annotation Pipeline. The results revealed substantial genomic diversity within the genus, with core gene counts varying across different nucleotide identity thresholds. Ribosomal genes such as rpsL, rpsJ, rpsS, rpmA, rpsK, rpsU, rpsG, rpmH, and rpsZ were consistently identified in the core genome, whereas accessory genes exhibited marked variability. This study uncovered novel and highly specific genetic markers for various Campylobacter species, including petB, clpX, and carB for C. coli; hypothetical proteins for C. jejuni and C. fetus; porA2 for C. lari; and mdtJ for C. upsaliensis. These markers demonstrated a specificity of at least 90% with minimal cross-reactivity with non-target organisms. The findings underscore the genomic heterogeneity within Campylobacter and provide essential genetic targets for the enhanced molecular detection of its pathogenic species, subspecies, and biovars.

Multi-country and intersectoral assessment of cluster congruence between pipelines for genomics surveillance of foodborne pathogens

Different laboratories employ different Whole-Genome Sequencing (WGS) pipelines for Food and Waterborne disease (FWD) surveillance, casting doubt on the comparability of their results and hindering optimal communication at intersectoral and international levels. Through a collaborative effort involving eleven European institutes spanning the food, animal, and human health sectors, we aimed to assess the inter-pipeline clustering congruence across all resolution levels and perform an in-depth comparative analysis of cluster composition at outbreak level for four important foodborne pathogens: Listeria monocytogenes, Salmonella enterica, Escherichia coli, and Campylobacter jejuni. We found a general concordance between allele-based pipelines for all species, except for C. jejuni, where the different resolution power of allele-based schemas led to marked discrepancies. Still, we identified non-negligible differences in outbreak detection and demonstrated how a threshold flexibilization favors the detection of similar outbreak signals by different laboratories. These results, together with the observation that different traditional typing groups (e.g., serotypes) exhibit a remarkably different genetic diversity, represent valuable information for future outbreak case-definitions and WGS-based nomenclature design. This study reinforces the need, while demonstrating the feasibility, of conducting continuous pipeline comparability assessments, and opens good perspectives for a smoother international and intersectoral cooperation towards an efficient One Health FWD surveillance.

PanKB: An interactive microbial pangenome knowledgebase for research, biotechnological innovation, and knowledge mining

The exponential growth of microbial genome data presents unprecedented opportunities for unlocking the potential of microorganisms. The burgeoning field of pangenomics offers a framework for extracting insights from this big biological data. Recent advances in microbial pangenomic research have generated substantial data and literature, yielding valuable knowledge across diverse microbial species. PanKB (pankb.org), a knowledgebase designed for microbial pangenomics research and biotechnological applications, was built to capitalize on this wealth of information. PanKB currently includes 51 pangenomes from 8 industrially relevant microbial families, comprising 8402 genomes, over 500 000 genes and over 7M mutations. To describe this data, PanKB implements four main components: (1) Interactive pangenomic analytics to facilitate exploration, intuition, and potential discoveries; (2) Alleleomic analytics, a pangenomic-scale analysis of variants, providing insights into intra-species sequence variation and potential mutations for applications; (3) A global search function enabling broad and deep investigations across pangenomes to power research and bioengineering workflows; (4) A bibliome of 833 open-access pangenomic papers and an interface with an LLM that can answer in-depth questions using its knowledge. PanKB empowers researchers and bioengineers to harness the potential of microbial pangenomics and serves as a valuable resource bridging the gap between pangenomic data and practical applications.

Genome characterisation and comparative analysis of Schaalia dentiphila sp. nov. and its subspecies, S. dentiphila subsp. denticola subsp. nov., from the human oral cavity

BACKGROUND

Schaalia species are primarily found among the oral microbiota of humans and other animals. They have been associated with various infections through their involvement in biofilm formation, modulation of host responses, and interaction with other microorganisms. In this study, two strains previously indicated as Actinomyces spp. were found to be novel members of the genus Schaalia based on their whole genome sequences.

RESULTS

Whole-genome sequencing revealed both strains with a genome size of 2.3 Mbp and GC contents of 65.5%. Phylogenetics analysis for taxonomic placement revealed strains NCTC 9931 and C24 as distinct species within the genus Schaalia. Overall genome-relatedness indices including digital DNA-DNA hybridization (dDDH), and average nucleotide/amino acid identity (ANI/AAI) confirmed both strains as distinct species, with values below the species boundary thresholds (dDDH < 70%, and ANI and AAI < 95%) when compared to nearest type strain Schaalia odontolytica NCTC 9935 T. Pangenome and orthologous analyses highlighted their differences in gene properties and biological functions compared to existing type strains. Additionally, the identification of genomic islands (GIs) and virulence-associated factors indicated their genetic diversity and potential adaptive capabilities, as well as potential implications for human health. Notably, CRISPR-Cas systems in strain NCTC 9931 underscore its adaptive immune mechanisms compared to strain C24.

CONCLUSIONS

Based on these findings, strain NCTC 9931T (= ATCC 17982T = DSM 43331T = CIP 104728T = CCUG 18309T = NCTC 14978T = CGMCC 1.90328T) represents a novel species, for which the name Schaalia dentiphila subsp. dentiphila sp. nov. subsp. nov. is proposed, while strain C24T (= NCTC 14980T = CGMCC 1.90329T) represents a distinct novel subspecies, for which the name Schaalia dentiphila subsp. denticola. subsp. nov. is proposed. This study enriches our understanding of the genomic diversity of Schaalia species and paves the way for further investigations into their roles in oral health.

SIGNIFICANCE

This research reveals two Schaalia strains, NCTC 9931 T and C24T, as novel entities with distinct genomic features. Expanding the taxonomic framework of the genus Schaalia, this study offers a critical resource for probing the metabolic intricacies and resistance patterns of these bacteria. This work stands as a cornerstone for microbial taxonomy, paving the way for significant advances in clinical diagnostics.

Uncovering the boundaries of Campylobacter species through large-scale phylogenetic and nucleotide identity analyses

Campylobacter species are typically helical shaped, Gram-negative, and non-spore-forming bacteria. Species in this genus include established foodborne and animal pathogens as well as emerging pathogens. The accumulation of genomic data from the Campylobacter genus has increased exponentially in recent years, accompanied by the discovery of putative new species. At present, the lack of a standardized species boundary complicates distinguishing established and novel species. We defined the Campylobacter genus core genome (500 loci) using publicly available Campylobacter complete genomes (n = 498) and constructed a core genome phylogeny using 2,193 publicly available Campylobacter genomes to examine inter-species diversity and species boundaries. Utilizing 8,440 Campylobacter genomes representing 33 species and 8 subspecies, we found species delineation based on an average nucleotide identity (ANI) cutoff of 94.2% is consistent with the core genome phylogeny. We identified 60 ANI genomic species that delineated Campylobacter species in concordance with previous comparative genetic studies. All pairwise ANI genomic species pairs had in silico DNA-DNA hybridization scores of less than 70%, supporting their delineation as separate species. We provide the tool Campylobacter Genomic Species typer (CampyGStyper) that assigns ANI genomic species to query genomes based on ANI similarities to medoid genomes from each ANI genomic species with an accuracy of 99.96%. The ANI genomic species definitions proposed here allow consistent species definition in the Campylobacter genus and will facilitate the detection of novel species in the future.IMPORTANCEIn recent years, Campylobacter has gained recognition as the leading cause of bacterial gastroenteritis worldwide, leading to a substantial rise in the collection of genomic data of the Campylobacter genus in public databases. Currently, a standardized Campylobacter species boundary at the genomic level is absent, leading to challenges in detecting emerging pathogens and defining putative novel species within this genus. We used a comprehensive representation of genomes of the Campylobacter genus to construct a core genome phylogenetic tree. Furthermore, we found an average nucleotide identity (ANI) of 94.2% as the optimal cutoff to define the Campylobacter species. Using this cutoff, we identified 60 ANI genomic species which provided a standardized species definition and nomenclature. Importantly, we have developed Campylobacter Genomic Species typer (CampyGStyper), which can robustly and accurately assign these ANI genomic species to Campylobacter genomes, thereby aiding pathogen surveillance and facilitating evolutionary and epidemiological studies of existing and emerging pathogens in the genus Campylobacter.

(Pan)genomic analysis of two Rhodococcus isolates and their role in phenolic compound degradation

The genus Rhodococcus is recognized for its potential to degrade a large range of aromatic substances, including plant-derived phenolic compounds. We used comparative genomics in the context of the broader Rhodococcus pan-genome to study genomic traits of two newly described Rhodococcus strains (type-strain Rhodococcus pseudokoreensis R79T and Rhodococcus koreensis R85) isolated from apple rhizosphere. Of particular interest was their ability to degrade phenolic compounds as part of an integrated approach to treat apple replant disease (ARD) syndrome. The pan-genome of the genus Rhodococcus based on 109 high-quality genomes was open with a small core (1.3%) consisting of genes assigned to basic cell functioning. The range of genome sizes in Rhodococcus was high, from 3.7 to 10.9 Mbp. Genomes from host-associated strains were generally smaller compared to environmental isolates which were characterized by exceptionally large genome sizes. Due to large genomic differences, we propose the reclassification of distinct groups of rhodococci like the Rhodococcus equi cluster to new genera. Taxonomic species affiliation was the most important factor in predicting genetic content and clustering of the genomes. Additionally, we found genes that discriminated between the strains based on habitat. All members of the genus Rhodococcus had at least one gene involved in the pathway for the degradation of benzoate, while biphenyl degradation was mainly restricted to strains in close phylogenetic relationships with our isolates. The ~40% of genes still unclassified in larger Rhodococcus genomes, particularly those of environmental isolates, need more research to explore the metabolic potential of this genus.IMPORTANCERhodococcus is a diverse, metabolically powerful genus, with high potential to adapt to different habitats due to the linear plasmids and large genome sizes. The analysis of its pan-genome allowed us to separate host-associated from environmental strains, supporting taxonomic reclassification. It was shown which genes contribute to the differentiation of the genomes based on habitat, which can possibly be used for targeted isolation and screening for desired traits. With respect to apple replant disease (ARD), our isolates showed genome traits that suggest potential for application in reducing plant-derived phenolic substances in soil, which makes them good candidates for further testing against ARD.

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.

Antibiotic resistance, plasmids, and virulence-associated markers in human strains of Campylobacter jejuni and Campylobacter coli isolated in Italy

Campylobacteriosis, a prevalent foodborne gastrointestinal infection in Europe, is primarily caused by Campylobacter jejuni and Campylobacter coli, with rising global concerns over antimicrobial resistance in these species. This study comprehensively investigates 133 human-origin Campylobacter spp. strains (102 C. jejuni and 31 C. coli) collected in Italy from 2013 to 2021. The predominant Multilocus Sequence Typing Clonal complexes (CCs) were ST-21 CC and ST-206 CC in C. jejuni and ST-828 CC in C. coli. Ciprofloxacin and tetracycline resistance, mainly attributed to GyrA (T86I) mutation and tet(O) presence, were prevalent, while erythromycin resistance was associated with 23S rRNA gene mutation (A2075G), particularly in C. coli exhibiting multidrug-resistant pattern CipTE. Notable disparities in virulence factors among strains were observed, with C. jejuni exhibiting a higher abundance compared to C. coli. Notably, specific C. jejuni sequence types, including ST-21, ST-5018, and ST-1263, demonstrated significantly elevated counts of virulence genes. This finding underscores the significance of considering both the species and strain-level variations in virulence factor profiles, shedding light on potential differences in the pathogenicity and clinical outcomes associated with distinct C. jejuni lineages. Campylobacter spp. plasmids were classified into three groups comprising pVir-like and pTet-like plasmids families, exhibiting diversity among Campylobacter spp. The study underscores the importance of early detection through Whole Genome Sequencing to identify potential emergent virulence, resistance/virulence plasmids, and new antimicrobial resistance markers. This approach provides actionable public health data, supporting the development of robust surveillance programs in Italy.

Characterising the gut microbiome of stranded harbour seals (Phoca vitulina) in rehabilitation

Animal rehabilitation centres provide a unique opportunity to study the microbiome of wild animals because subjects will be handled for their treatment and can therefore be sampled longitudinally. However, rehabilitation may have unintended consequences on the animals' microbiome because of a less varied and suboptimal diet, possible medical treatment and exposure to a different environment and human handlers. Our study describes the gut microbiome of two large seal cohorts, 50 pups (0-30 days old at arrival) and 23 weaners (more than 60 days old at arrival) of stranded harbour seals admitted for rehabilitation at the Sealcentre Pieterburen in the Netherlands, and the effect of rehabilitation on it. Faecal samples were collected from all seals at arrival, two times during rehabilitation and before release. Only seals that did not receive antimicrobial treatment were included in the study. The average time in rehabilitation was 95 days for the pups and 63 days for the weaners. We observed that during rehabilitation, there was an increase in the relative abundance of some of the Campylobacterota spp and Actinobacteriota spp. The alpha diversity of the pups' microbiome increased significantly during their rehabilitation (p-value <0.05), while there were no significant changes in alpha diversity over time for weaners. We hypothesize that aging is the main reason for the observed changes in the pups' microbiome. At release, the sex of a seal pup was significantly associated with the microbiome's alpha (i.e., Shannon diversity was higher for male pups, p-value <0.001) and beta diversity (p-value 0.001). For weaners, variation in the microbiome composition (beta diversity) at release was partly explained by sex and age of the seal (p-values 0.002 and 0.003 respectively). We mainly observed variables known to change the gut microbiome composition (e.g., age and sex) and conclude that rehabilitation in itself had only minor effects on the gut microbiome of seal pups and seal weaners.

Contraction and expansion dynamics: deciphering genomic underpinnings of growth rate and pathogenicity in Mycobacterium

BACKGROUND

Mycobacterium bacteria, encompassing both slow growth (SGM) and rapid growth mycobacteria (RGM), along with true pathogenic (TP), opportunistic pathogenic (OP), and non-pathogenic (NP) types, exhibit diverse phenotypes. Yet, the genetic underpinnings of these variations remain elusive.

METHODS

Here, We conducted a comprehensive comparative genomics study involving 53 Mycobacterium species to unveil the genomic drivers behind growth rate and pathogenicity disparities.

RESULTS

Our core/pan-genome analysis highlighted 1,307 shared gene families, revealing an open pan-genome structure. A phylogenetic tree highlighted clear boundaries between SGM and RGM, as well as TP and other species. Gene family contraction emerged as the primary alteration associated with growth and pathogenicity transitions. Specifically, ABC transporters for amino acids and inorganic ions, along with quorum sensing genes, exhibited significant contractions in SGM species, potentially influencing their distinct traits. Conversely, TP strains displayed contraction in lipid and secondary metabolite biosynthesis and metabolism-related genes. Across the 53 species, we identified 26 core and 64 accessory virulence factors. Remarkably, TP and OP strains stood out for their expanded mycobactin biosynthesis and type VII secretion system gene families, pivotal for their pathogenicity.

CONCLUSION

Our findings underscore the importance of gene family contraction in nucleic acids, ions, and substance metabolism for host adaptation, while emphasizing the significance of virulence gene family expansion, including type VII secretion systems and mycobactin biosynthesis, in driving mycobacterial pathogenicity.

Genetic characteristics, antimicrobial susceptibility, and virulence genes distribution of Campylobacter isolated from local dual-purpose chickens in central China

Food-borne antibiotic-resistant Campylobacter poses a serious threat to public health. To understand the prevalence and genetic characteristics of Campylobacter in Chinese local dual-purpose (meat and eggs) chickens, the genomes of 30 Campylobacter isolates, including 13 C. jejuni and 17 C. coli from Jianghan-chickens in central China, were sequenced and tested for antibiotic susceptibility. The results showed that CC-354 and CC-828 were the dominant clonal complexes of C. jejuni and C. coli, respectively, and a phylogenetic analysis showed that three unclassified multilocus sequence types of C. coli were more closely genetically related to C. jejuni than to other C. coli in this study. Of the six antibiotics tested, the highest resistance rates were to ciprofloxacin and tetracycline (100%), followed by lincomycin (63.3%), erythromycin (30.0%), amikacin (26.7%), and cefotaxime (20.0%). The antibiotic resistance rate of C. coli was higher than that of C. jejuni. The GyrA T86I mutation and 15 acquired resistance genes were detected with whole-genome sequencing (WGS). Among those, the GyrA T86I mutation and tet(O) were most prevalent (both 96.7%), followed by the blaOXA-type gene (90.0%), ant(6)-Ia (26.7%), aac(6')-aph(3'') (23.3%), erm(B) (13.3%), and other genes (3.3%). The ciprofloxacin and tetracycline resistance phenotypes correlated strongly with the GyrA T86I mutation and tet(O)/tet(L), respectively, but for other antibiotics, the correlation between genes and resistance phenotypes were weak, indicating that there may be resistance mechanisms other than the resistance genes detected in this study. Virulence gene analysis showed that several genes related to adhesion, colonization, and invasion (including cadF, porA, ciaB, and jlpA) and cytolethal distending toxin (cdtABC) were only present in C. jejuni. Overall, this study extends our knowledge of the epidemiology and antibiotic resistance of Campylobacter in local Chinese dual-purpose chickens.

Application of TraDIS to define the core essential genome of Campylobacter jejuni and Campylobacter coli

Campylobacter species are the major cause of bacterial gastroenteritis. As there is no effective vaccine, combined with the rapid increase in antimicrobial resistant strains, there is a need to identify new targets for intervention. Essential genes are those that are necessary for growth and/or survival, making these attractive targets. In this study, comprehensive transposon mutant libraries were created in six C. jejuni strains, four C. coli strains and one C. lari and C. hyointestinalis strain, allowing for those genes that cannot tolerate a transposon insertion being called as essential. Comparison of essential gene lists using core genome analysis can highlight those genes which are common across multiple strains and/or species. Comparison of C. jejuni and C. coli, the two species that cause the most disease, identified 316 essential genes. Genes of interest highlighted members of the purine pathway being essential for C. jejuni whilst also finding that a functional potassium uptake system is essential. Protein-protein interaction networks using these essential gene lists also highlighted proteins in the purine pathway being major 'hub' proteins which have a large number of interactors across the network. When adding in two more species (C. lari and C. hyointestinalis) the essential gene list reduces to 261. Within these 261 essential genes, there are many genes that have been found to be essential in other bacteria. These include htrB and PEB4, which have previously been found as core virulence genes across Campylobacter species in other studies. There were 21 genes which have no known function with eight of these being associated with the membrane. These surface-associated essential genes may provide attractive targets. The essential gene lists presented will help to prioritise targets for the development of novel therapeutic and preventative interventions.