Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Microbiol., 17 July 2025

Sec. Antimicrobials, Resistance and Chemotherapy

Volume 16 - 2025 | https://doi.org/10.3389/fmicb.2025.1496275

This article is part of the Research TopicHorizontal Transfer of Antibiotic Resistance Genes in the Environment: Dynamic, Contributing Factors, and ControlView all 9 articles

The restriction impacts of the Type III restriction-modification system on the transmission dynamics of antimicrobial resistance genes in Campylobacter jejuni

Yu QiuYu Qiu1Pengbo GuoPengbo Guo2Hui TianHui Tian3Ye ZhouYe Zhou1Hongling Wen
Hongling Wen1*Hao Liang
Hao Liang1*
  • 1Department of Health Inspection and Quarantine, School of Public Health, Shandong University, Jinan, China
  • 2Department of Pathogen Biology, School of Clinical and Basic Medical Sciences, Shandong First Medical University and Shandong Academy of Medical Sciences, Jinan, China
  • 3Department of Radiation Oncology, Qilu Hospital of Shandong University, Cheeloo College of Medicine, Shandong University, Jinan, China

Introduction: The spread of antibiotic resistance genes among Campylobacter jejuni (C. jejuni) is a serious problem, and the effects of the restriction-modification (R-M) system on the transmission dynamics of these genes in C. jejuni remain poorly understood.

Materials and methods: Complete genome sequences of C. jejuni strains were extracted from the BV-BRC database until March 25, 2024. The phylogenetic and the resistance analysis were used to analyze the distribution of resistance genes in C. jejuni. The impacts of the R-M systems on the AMR genes transmission between C. jejuni strains and the possible mechanisms were explored through recombination, pangenome and mobile genetic elements analysis.

Results: C. jejuni strains carrying the Type III R-M system have a significantly lower number of antimicrobial resistance (AMR) genes compared to strains without this system (p < 0.0001), with covariance value being −0.0526. The recombination analysis also shows that the median number of the number of AMR genes in the strains not possessing the Type III R-M system increases by 19.38% compared to strains carrying that system (p < 0.0001). We also find that the horizontal gene transfer frequency might have limited relationship with the Type III R-M system in C. jejuni through pangenome and mobile genetic elements analysis.

Conclusion: Our research indicates that the Type III R-M system might restrict the transmission of AMR genes potentially by affecting recombination in C. jejuni, which provides a theoretical basis for addressing the drug resistance problem.

1 Introduction

Campylobacter jejuni (C. jejuni) is one of the leading bacterial pathogens causing foodborne diarrhea worldwide (Igwaran and Okoh, 2019; European Food Safety Authority, European Centre for Disease Prevention and Control, 2023). This pathogenic bacterium can also trigger extraintestinal infections such as bacteremia (Graham et al., 2024), myocarditis (Jiffry et al., 2023), in more severe cases, irritable bowel syndrome (Takakura et al., 2022) and Guillain-Barré syndrome (Finsterer, 2022). Antibiotics should be considered when C. jejuni infection is severe. However, with the misuse of antibiotics in human and veterinary medicine, C. jejuni has become increasingly resistant to antimicrobial agents used in animals and clinical settings. Despite several years of active surveillance of C. jejuni and reduced antibiotics usage in livestock production around the world, antimicrobial resistance (AMR) in C. jejuni remains a great public health challenge globally (Devi et al., 2019; van Vliet et al., 2022; Gao et al., 2023). One report has indicated that AMR markers were detected in 53% of C. jejuni isolates in the UK and US from 2001 to 2018 (van Vliet et al., 2022).

Homologous recombination is a genetic process in which DNA molecules exchange genetic material between two similar or identical DNA sequences, playing a significant role in the evolution of antibiotic resistance in bacteria by facilitating the exchange and integration of genetic material that confers resistance traits. And AMR genes mediating AMR carried by mobile genetic elements (MGEs) in C. jejuni could be spread within different bacterial species, exacerbating bacterial resistance to antibiotics (Partridge et al., 2018; Shikov et al., 2023). A mobile florfenicol resistance gene in C. jejuni, fexA, encodes an original florfenicol efflux pump system, conferring a high level of florfenicol resistance. The florfenicol resistance gene is located within a region featuring the tet (L)-fexA-catA-tet (O) arrangement, which has been shown to be transferable among C. jejuni population, thereby worsening the florfenicol resistance issue (Tang et al., 2020).

In addressing the significant issue of AMR in C. jejuni, it is imperative to investigate not only the transmission processes of resistant bacteria but also the dynamics of AMR gene transfer between C. jejuni strains. The entry mechanisms for exogenous DNA into a bacterium encompass transformation, conjugation and transduction, with the restriction-modification (R-M) system functioning as a defense against exogenous genetic material invasion (Ershova et al., 2015; Tao et al., 2022). The R-M system is of great significance for regulating the uptake and integration of exogenous DNA in bacteria. It comprises two active enzymes, restriction endonuclease (REase) which cleaves the target DNA sequences (Pingoud et al., 2014), and methyltransferase (MTase) which methylates the specific DNA sequences to protect them from REase hydrolysis (Oliveira and Fang, 2021; Seong et al., 2021), thus forming a natural barrier against foreign genetic materials invasion. Based on the specific DNA locus recognized, molecular structure and cofactor requirements, the R-M system can be roughly classified into four groups: Type I ~ IV R-M system (Isaev et al., 2021). Type I REases are multisubunit enzymes comprising three functionally distinct components: a sequence recognition subunit, a MTase subunit, and a REase subunit (Loenen et al., 2013). The REases characteristically interact with two cognate recognition sites and induce DNA cleavage at a midpoint between the recognition sites. Type II REases exist as homodimeric or tetrameric complexes and mediate sequence-specific DNA cleavage either within or adjacent to their recognition sites (typically 4–8 bp) (Pingoud et al., 2014). Type III REases consist of two subunits: a DNA recognition/modification subunit and a distinct DNA cleavage subunit. The enzyme interacts with two recognition sites and leads to a break at a location relative to one bound recognition sequence (Rao et al., 2013). Type IV Reases specifically cleave modified DNA (Loenen and Raleigh, 2014). AveC4I, a Type I MTase, might increase drug resistance in Aeromonas veronii (Ma et al., 2023). One research also proved that Cj1051c, a typeIIputative R-M enzyme, could drastically reduce the conjugation efficiency of the standard C. jejuni strain NCTC 11168 (Zeng et al., 2018). Nevertheless, there are few studies about the impact of the R-M system on the distribution of AMR genes in C. jejuni and the underlying mechanisms.

This study aimed to investigate how the R-M system influences the transfer of AMR genes between C. jejuni strains to better control the spread of resistance genes. The findings could offer a theoretical foundation for managing the transmission of AMR genes, which holds significant implications for public health.

2 Materials and methods

2.1 Campylobacter jejuni genome

One thousand and ninety-two C. jejuni strains with high sequencing quality and complete information isolated from Asia, Europe and North America were selected from the BV-BRC database until March 25, 2024 (Olson et al., 2023). Subsequently, their genome files were extracted for our research (Supplementary Table S1).

2.2 Phylogenetic tree

Campylobacter jejuni genome files extracted from the BV-BRC database were used to construct the core gene phylogenetic tree. Through Snippy, C. jejuni genomes were compared with the genome file of C. jejuni strain NCTC11168 to generate a core Single Nucleotide Polymorphism (SNP) alignment. Subsequently, Gubbins was used to remove SNP sites in the recombinant regions to obtain a high-quality core SNP alignment (Croucher et al., 2014). Then, RAxML was used to generate the final phylogenomic tree (Stamatakis, 2014), which could be visualized via the interactive website tvBOT (Xie et al., 2023).

2.3 Prediction of the R-M system in the genomes

To obtain the impacts of the R-M system on the dynamics of the transfer of AMR genes, the REase and MTase genes related to R-M system were extracted from the REBASE database (Roberts et al., 2023). Using these genes collected from the REBASE database as reference, we compared the C. jejuni genome files downloaded with them to obtain the distribution of R-M system. We created a sample pool for local screening of C. jejuni strains with the R-M system. BLAST was performed based on the C. jejuni genome files downloaded from the BV-BRC database (Supplementary Table S2). The interactive website tvBOT was used to visualize the distribution of the R-M system in C. jejuni.

2.4 Prediction of AMR genes in the genomes

Choosing the National Center for Biotechnology Information (NCBI) database as the reference database (Feldgarden et al., 2019), acquired AMR genes in C. jejuni genome files were predicted using the ABRicate tool (Supplementary Table S3). The distribution of AMR genes in C. jejuni was visualized via tvBOT, and to explore the impacts of the R-M system on AMR genes, RStudio packages including ggplot2, RColorBrewer and ggpubr were used to analyze and visualize the differences in the number of AMR genes between C. jejuni with or without the R-M system. The significance was analyzed using the Wilcoxon-Mann–Whitney test.

2.5 Homologous recombination analysis

The recombination_predictions.gff file and per_branch_statistics.csv file were both generated from the Gubbins process. The former one was used to visualize the homologous recombination predictions in C. jejuni via phandango (Hadfield et al., 2017). The latter one contained summary statistics for the number of recombination events reconstructed onto both the internal nodes and the leaves of the phylogenetic tree and along with the.recombination_predictions.embl file which covered detailed information about recombination predictions, we could calculate the total number of recombination events in one C. jejuni strain (Supplementary Tables S4, S5). It was used to analyze the differences of the recombination blocks between C. jejuni with or without the R-M system via RStudio packages including ggplot2, RColorBrewer and ggpubr. The significance was analyzed using the Wilcoxon-Mann–Whitney test.

2.6 Pangenome analysis

Campylobacter jejuni genome files downloaded from the BV-BRC database were initially annotated using Prokka (Seemann, 2014). Subsequently, Pangenome analysis was performed using Roary on the generated GFF files (Page et al., 2015). To access the impacts of the R-M system on the distribution of unique and cloud genes in C. jejuni, using RStudio packages including ggplot2, RColorBrewer and ggpubr, we analyzed and visualized the differences of the number of these genes between C. jejuni with or without R-M system based on the gene_presence_absence.csv file obtained from the Roary process (Supplementary Table S6). The significance was analyzed using the Wilcoxon-Mann–Whitney test.

2.7 Prediction of integrating MGEs

MGEs in C. jejuni genomes were predicted using MobileElementFinder with default parameters (Johansson et al., 2020) (Supplementary Table S7). Then, RStudio packages including ggplot2, RColorBrewer and ggpubr were employed to analyze and visualize the differences in the number of MGEs between C. jejuni with or without the R-M system. The significance was analyzed using the Wilcoxon-Mann–Whitney test.

2.8 Covariance calculation

Covariance is calculated using the following formula: Cov (x, y) in Rstudio. The x is the value of AMR genes or MGEs, and for y, it takes the value 1 when the Type III R-M system is present or 0 otherwise. The covariance value represents the changing trend of two-dimensional random variables. If the value is positive, the changing trend is consistent, or, if not, the changing trend is opposite.

3 Results

3.1 The impacts of the R-M system on the distribution of AMR genes in Campylobacter jejuni

Among C. jejuni strains selected for our study, according to the geographical distribution, 387 (387/1,092, 35.44%), 356 (356/1,092, 32.6%), and 349 (349/1,092, 31.96%) strains were isolated from Europe, North America and Asia, respectively. According to the host origin distribution, 397 (397/1,092, 36.45%), 340 (340/1,092, 31.14%), 167 (167/1,092, 15.29%), 91 (91/1,092, 8.33%), 54 (54/1,092, 4.95%), 35 (35/1,092, 3.21%) and 8 (8/1,092, 0.73%) strains were isolated from patients, poultry, livestock, birds, monkeys, the environment and black bears, respectively.

First, we analyzed the distribution of the R-M system in C. jejuni. Utilizing the REase and MTase genes related to R-M system extracted from the REBASE database, we constructed a sample pool for local screening of C. jejuni strains with the R-M system. BLAST was performed based on the C. jejuni genome files downloaded from the BV-BRC database. We found that 1,043 (1,043/1,092, 95.51%), 1,092 (1,092/1,092, 100%), 334 (334/1,092, 30.59%) and 1,080 (1,080/1,092, 98.90%) strains carry the Type I, Type II, Type III, and Type IV R-M systems, respectively. The Type I and Type III R-M system were distributed unevenly in C. jejuni.

We explored the distribution of AMR genes and the impacts of the R-M systems on the distribution of AMR genes in C. jejuni. Using the ABRicate tool, we identified the acquired AMR genes in C. jejuni and quantified them in the genome of the strains (Figure 1). We divided these strains into groups based on the presence of Type I R-M or Type III R-M systems, as well as whether they are R-M-free. We then analyzed the differences in the number distribution of the AMR genes between these groups (Figure 2). The results revealed that C. jejuni strains carrying the Type III R-M system had a significantly lower number of AMR genes compared to strains without this system (Wilcoxon-Mann–Whitney test, p < 0.0001). To evaluate the directional relationship between the number distribution of AMR genes and the presence of the Type III R-M system, we estimated the covariance value, which shows the relationship between two random variables. The covariance analysis revealed that the covariance value between the Type III R-M group and the Type III R-M-free group was −0.0526, suggesting a negative correlation between the number of AMR genes and the presence of the Type III R-M system. However, there were no significant differences in the distribution of the AMR genes between C. jejuni groups with a Type I R-M system or without it (Wilcoxon-Mann–Whitney test, p = 0.59).

Figure 1
Phylogenetic tree and heatmap depicting the relationship between samples based on antimicrobial resistance (AMR) genes. The tree scale is ten thousand. The heatmap shows variations in AMR genes across samples, grouped by continent and host origin, including livestock, patients, and various animals. Resistance types are indicated by color, with a key provided. Columns represent different gene types like aminoglycoside and tetracycline.

Figure 1. Heatmap of the distribution of the AMR genes in C. jejuni based on the phylogenetic tree. The rows in the heatmap represent different C. jejuni strains, and the columns represent different AMR genes. The color of the blocks represents the number of the specific AMR gene in one C. jejuni strain.

Figure 2
Violin plots comparing the number of antimicrobial genes in R-M-free and R-M types. Panel A shows Type I with no significant difference (ns), and Panel B shows Type III with a significant difference indicated by asterisks. Red indicates R-M-free, blue indicates R-M types.

Figure 2. Differences of AMR genes in C. jejuni with or without (A) the Type I and (B) the Type III R-M system. The thicker lines in the box plot represent the median. The box plot’s upper and lower boundaries indicate the 75th (Q3) and 25th (Q1) percentiles, respectively. According to Wilcoxon-Mann–Whitney test, the number of asterisks indicates the significant differences: * indicates p ≤ 0.05, ** indicates p ≤ 0.01, *** indicates p ≤ 0.001, **** indicates p ≤ 0.0001, and ns indicates no significant difference.

3.2 The existence of the Type III R-M system is responsible for a fewer number of homologous recombination blocks in Campylobacter jejuni

Recombination, as a critical evolutionary force in bacteria, shapes the genomic landscape (Shikov et al., 2023). Using Gubbins, we reconstructed the homologous recombination blocks on the branches of the phylogenetic tree (Figure 3). The total number of the recombination blocks in each C. jejuni strain was counted. Through our research, the median number of homologous recombination blocks in the strains possessing the Type III R-M was 500.5, but in the Type III R-M-free strains, the median number increased by 19.38% to 597.5 (Figure 4, Wilcoxon-Mann–Whitney test, p < 0.0001), and the covariance value between these two groups was −14.1427, implying a negative correlation between them, which is consistent with the conclusion we have gained in our previous study. This suggests that the presence of the Type III R-M system could reduce the homologous recombination frequency in C. jejuni.

Figure 3
Dendrogram on the left classifies samples with linked presence (purple) and absence (white) matrix for Type I to IV R-M systems. Heatmap of varying red intensities represents data distribution, with a black line graph at the bottom showing frequency peaks.

Figure 3. Heatmap of the recombination in C. jejuni. The rows in the heatmap represent different C. jejuni strain genomes. The blue and red blocks represent recombination events reconstructed on the leaves and internal nodes, respectively.

Figure 4
Violin plot comparing the number of recombination events between Type III R-M-free (red) and Type III R-M (blue) systems. The Y-axis ranges from zero to 1200. Statistical significance is indicated with stars above the plot.

Figure 4. Differences of recombination events in C. jejuni with or without the Type III R-M system. The thicker lines in the box plot represent the median. The box plot’s upper and lower boundaries indicate the 75th (Q3) and 25th (Q1) percentiles, respectively. The significance is analyzed using the Wilcoxon-Mann–Whitney test, and the number of asterisks indicates the significant differences: * indicates p ≤ 0.05, ** indicates p ≤ 0.01, *** indicates p ≤ 0.001, **** indicates p ≤ 0.0001, and ns indicates no significant difference.

3.3 The Type III R-M system might have limited effects on the HGT frequency in Campylobacter jejuni

Cloud genes are defined as accessory genes found in 0 ~ 15% ofbacterial species. The HGT frequency of a C. jejuni strain can be evaluated by the number of the cloud and unique genes it carries (Xu et al., 2023). We analyzed the distribution of these genes between C. jejuni groups with or without the Type III R-M system to explore how the HGT frequency was affected by this system in C. jejuni population (Figure 5). Our analysis revealed no significant relationship between the number of cloud and unique genes and the presence of the Type III R-M system (Wilcoxon-Mann–Whitney test, p = 0.061 and p = 0.41, respectively).

Figure 5
Violin plots compare

Figure 5. Differences of (A) the cloud and (B) unique genes in C. jejuni with or without the Type III R-M system. The thicker lines in the box plot represent the median. The box plot’s upper and lower boundaries indicate the 75th (Q3) and 25th (Q1) percentiles, respectively. The significance is analyzed using the Wilcoxon-Mann–Whitney test, and the number of asterisks indicates the significant differences: * indicates p ≤ 0.05, ** indicates p ≤ 0.01, *** indicates p ≤ 0.001, **** indicates p ≤ 0.0001, and ns indicates no significant difference.

AMR genes carried by MGEs can spread within C. jejuni population via HGT (Partridge et al., 2018). We compared the number of the MGEs in C. jejuni strains with or without the Type III R-M system to investigate whether the suppression of the acquired AMR gene distribution by Type III R-M system is related to a decline in HGT frequency (Figure 6). The results showed no significant differences in the number of the MGEs between C. jejuni groups with or without the Type III R-M system (Wilcoxon-Mann–Whitney test, p = 0.13).

Figure 6
Violin plot comparing the number of mobile genetic elements between Type III R-M-free (red) and Type III R-M (blue). Both distributions show a similar spread, with no significant difference noted (ns). The y-axis is labeled

Figure 6. Differences of the MGE in C. jejuni with or without the Type III R-M system. The thicker lines in the box plot represent the median. The box plot’s upper and lower boundaries indicate the 75th (Q3) and 25th (Q1) percentiles, respectively. The significance is analyzed using the Wilcoxon-Mann–Whitney test, and the number of asterisks indicates the significant differences: * indicates p ≤ 0.05, ** indicates p ≤ 0.01, *** indicates p ≤ 0.001, **** indicates p ≤ 0.0001, and ns indicates no significant difference.

In conlusion, our results suggest that the Type III R-M system might have limited effects on HGT frequency in C. jejuni.

4 Discussion

4.1 The possible impact of the R-M system on the transfer of AMR genes in Campylobacter jejuni

Recently, C. jejuni AMR has become a severe global issue. Although the causes of C. jejuni AMR are complex, homologous recombination plays a central in enabling bacteria to acquire AMR elements (Partridge et al., 2018; Shikov et al., 2023). By integrating exogenous genetic materials, bacteria can pre-adapt to changing environment, thereby enhancing their viability. However, since gene acquisition can also lead to gene dysfunction and other side effects, bacteria have evolved ‘defense systems’ to protect their genetic integrity. The restriction-modification (R-M) system is the most prevalent defense mechanism, found in 83% of prokaryotic genomes, which is more than double the incidence of CRISPR-Cas systems (Tesson et al., 2022).

In this study, we first characterized the distribution of R-M systems in C. jejuni, revealing an uneven prevalence for both Type I (1,043/1,092, 95.51%) and Type III (334/1,092, 30.59%) systems. Given this disparity, to analyze the impact of R-M systems on AMR gene transmission in C. jejuni, we grouped strains based on the presence or absence of Type I and Type III R-M systems and compared the number of AMR genes between these groups. Our analysis revealed no significant differences in AMR gene distribution between strains possessing the Type I R-M system or not (Figure 2A; Wilcoxon-Mann–Whitney test, p = 0.59). Conversely, a significant negative correlation was identified between the presence of the Type III R-M system and the number of AMR genes in C. jejuni strains (Figure 2B; Wilcoxon-Mann–Whitney test, p < 0.0001; covariance = −0.0526). This finding suggests that the Type III R-M system may act as a barrier to AMR gene acquisition. The presence of that system may suppress the transmission of AMR genes by cleaving target exogenous DNA sequences, corresponding to the function of the REases.

4.2 The potential underlying mechanisms of the restriction effect on AMR genes transmission by the Type III R-M system in Campylobacter jejuni

AMR genes can spread within bacterial populations through homologous recombination, exacerbating the AMR problem (Partridge et al., 2018; Shikov et al., 2023). The macrolide resistance gene erm (B) could undergo natural transformation via homologous recombination among Campylobacter strains (Wei et al., 2024). Additionally, the horizontally acquired AMR genes could further be transferred via homologous recombination (Gratia, 2017). One study also shows that the horizontal transfer of chromosomal genes via lateral transduction through homologous recombination might provide significant opportunities for the rapid acquisition of virulence factors (Humphrey et al., 2021). To further explore mechanisms underlying the restriction effect on AMR genes transmission by the Type III R-M system, we conducted a series of analyses. Our results show that the presence of the Type III R-M system is associated with a decrease in the number of homologous recombination blocks in C. jejuni, strains lacking this system exhibiting a 19.38% higher median number of homologous recombination events (Figure 4; Wilcoxon-Mann–Whitney test, p < 0.0001; covariance = −14.1427). Thus, we propose that the Type III R-M system restricts the AMR genes transmission within C. jejuni mainly by limiting the frequency of homologous recombination. The homologous recombination frequency of the exogenous AMR genes in C. jejuni might be suppressed by the Type III R-M REases, which can cleave DNA sequences and restrict the transmission of AMR genes among C. jejuni population, resulting in fewer AMR genes in C. jejuni with the Type III R-M system.

Sharing of genes through HGT contributes importantly to the global dissemination of AMR genes, which could be susceptible to bacterial defense systems (Price et al., 2016; Huo et al., 2015). Proposing that the Type III R-M system might influence the acquisition of AMR genes by interfering with the HGT process, we conducted a series of studies. However, results show that there are no significant differences in the number of the cloud or unique genes between C. jejuni groups with or without the Type III R-M system (Figure 5; Wilcoxon-Mann–Whitney test, p = 0.061 and p = 0.41, respectively). MGEs, such as plasmids and the integrative and conjugative elements, disseminating via HGT, have been identified as key vehicles for the dissemination of AMR determinants (Botelho and Schulenburg, 2021). Through the MobileElementFinder software, we found and calculated the number of insertion sequences and unit transposons carried by C. jejuni strains with or without the Type III R-M system. Nevertheless, our research shows no significant difference between the number of MGEs in C. jejuni with or without that system (Figure 6; Wilcoxon-Mann–Whitney test, p = 0.13). Therefore, we suggest that MGEs might have limited correlation with the Type III R-M system in the transfer of AMR genes.

5 Conclusion

The analysis of the number distribution of AMR genes between the two groups of C. jejuni strains with or without type III R-M system shows that this system might have a restriction effect on AMR gene transfer in C. jejuni. Our data indicates that C. jejuni without that system seems to appear a higher recombination frequency, but there is no significant difference in the number of unique or cloud genes and MGEs between the two groups. Hence, the type III R-M system may have limited influence on the dissemination of AMR genes between other bacterial species and C. jejuni. However, it may restrict the exchange and integration of AMR genes within this species probably by limiting the frequency of homologous recombination to affect the distribution of AMR genes.

The transmission of AMR genes among bacterial populations is a complex process influenced by environmental factors, the genetic characteristics and the bacterial defense systems. This study suggests that the Type III R-M system shows a correlation with the dynamics of the AMR genes transmission between C. jejuni strains and explores the possible mechanisms underlying the distribution differences in AMR gene numbers in this species which will be verified through further experiments. Consequently, it provides a theoretical basis for addressing the antibiotic resistance problem in C. jejuni.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding authors.

Author contributions

YQ: Formal analysis, Methodology, Visualization, Writing – original draft, Writing – review & editing. PG: Data curation, Writing – review & editing. HT: Methodology, Writing – review & editing. YZ: Visualization, Writing – review & editing. HW: Validation, Supervision, Writing – review & editing. HL: Conceptualization, Funding acquisition, Writing – review & editing, Project administration.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This work was funded by Shandong Province Natural Science Foundation Youth Branch (ZR2023QC093).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2025.1496275/full#supplementary-material

References

Botelho, J., and Schulenburg, H. (2021). The role of integrative and conjugative elements in antibiotic resistance evolution. Trends Microbiol. 29, 8–18. doi: 10.1016/j.tim.2020.05.011

PubMed Abstract | Crossref Full Text | Google Scholar

Croucher, N. J., Page, A. J., Connor, T. R., Delaney, A. J., Keane, J. A., Bentley, S. D., et al. (2014). Rapid phylogenetic analysis of large samples of recombinant bacterial whole genome sequences using Gubbins. Nucleic Acids Res. 43:e15. doi: 10.1093/nar/gku1196

PubMed Abstract | Crossref Full Text | Google Scholar

Devi, A., Mahony, T. J., Wilkinson, J. M., and Vanniasinkam, T. (2019). Antimicrobial susceptibility of clinical isolates of Campylobacter jejuni from New South Wales, Australia. J. Glob. Antimicrob. Resist. 16, 76–80. doi: 10.1016/j.jgar.2018.09.011

PubMed Abstract | Crossref Full Text | Google Scholar

Ershova, A. S., Rusinov, I. S., Spirin, S. A., Karyagina, A. S., and Alexeevski, A. V. (2015). Role of restriction-modification systems in prokaryotic evolution and ecology. Biochemistry (Mosc) 80, 1373–1386. doi: 10.1134/s0006297915100193

PubMed Abstract | Crossref Full Text | Google Scholar

European Food Safety Authority, European Centre for Disease Prevention and Control (2023). The European Union one health 2022 zoonoses report. EFSA J. 21:e8442. doi: 10.2903/j.efsa.2023.8442

Crossref Full Text | Google Scholar

Feldgarden, M., Brover, V., Haft, D. H., Prasad, A. B., Slotta, D. J., Tolstoy, I., et al. (2019). Validating the AMRFinder tool and resistance gene database by using antimicrobial resistance genotype-phenotype correlations in a collection of isolates. Antimicrob. Agents Chemother. 63:e00483-19. doi: 10.1128/aac.00483-19

PubMed Abstract | Crossref Full Text | Google Scholar

Finsterer, J. (2022). Triggers of Guillain-Barre syndrome: Campylobacter jejuni predominates. Int. J. Mol. Sci. 23:14222. doi: 10.3390/ijms232214222

PubMed Abstract | Crossref Full Text | Google Scholar

Gao, F., Tu, L. H., Chen, M. L., Chen, H. Y., Zhang, X., Zhuang, Y., et al. (2023). Erythromycin resistance of clinical campylobacter jejuni and Campylobacter coli in Shanghai, China. Front. Microbiol. 14:1145581. doi: 10.3389/fmicb.2023.1145581

PubMed Abstract | Crossref Full Text | Google Scholar

Graham, A., Hawkins, L., Balasegaram, S., Narasimhan, S., Wain, J., Clarke, J., et al. (2024). A decade of campylobacter and campylobacter bacteraemias in a district general hospital and the surrounding London and south east region, England. J. Infect. 88, 15–20. doi: 10.1016/j.jinf.2023.11.004

PubMed Abstract | Crossref Full Text | Google Scholar

Gratia, J. P. (2017). Genetic recombinational events in prokaryotes and their viruses: insight into the study of evolution and biodiversity. Anton. Leeuw. Int. J. Gen. Mol. Microbiol. 110, 1493–1514. doi: 10.1007/s10482-017-0916-5

PubMed Abstract | Crossref Full Text | Google Scholar

Hadfield, J., Croucher, N. J., Goater, R. J., Abudahab, K., Aanensen, D. M., and Harris, S. R. (2017). Phandango: an interactive viewer for bacterial population genomics. Bioinformatics 34, 292–293. doi: 10.1093/bioinformatics/btx610

PubMed Abstract | Crossref Full Text | Google Scholar

Humphrey, S., Fillol Salom, A., Quiles Puchalt, N., Ibarra Chávez, R., Haag, A. F., Chen, J., et al. (2021). Bacterial chromosomal mobility via lateral transduction exceeds that of classical mobile genetic elements. Nat. Commun. 12:6509. doi: 10.1038/s41467-021-26004-5

Crossref Full Text | Google Scholar

Huo, W., Adams, H. M., Zhang, M. Q., and Palmer, K. L. (2015). Genome modification in Enterococcus faecalis OG1RF assessed by bisulfite sequencing and single-molecule real-time sequencing. J. Bacteriol. 197, 1939–1951. doi: 10.1128/jb.00130-15

PubMed Abstract | Crossref Full Text | Google Scholar

Igwaran, A., and Okoh, A. I. (2019). Human campylobacteriosis: a public health concern of global importance. Heliyon 5:e02814. doi: 10.1016/j.heliyon.2019.e02814

PubMed Abstract | Crossref Full Text | Google Scholar

Isaev, A. B., Musharova, O. S., and Severinov, K. V. (2021). Microbial arsenal of antiviral defenses part I. Biochemistry (Mosc.) 86, 319–337. doi: 10.1134/s0006297921030081

PubMed Abstract | Crossref Full Text | Google Scholar

Jiffry, M. Z. M., Okam, N. A., Vargas, J., Adekunle, F. A., Pagan, S. C., Khowaja, F., et al. (2023). Myocarditis as a complication of Campylobacter jejuni-associated enterocolitis: a report of two cases. Cureus J. Med. Sci. 15:e36171. doi: 10.7759/cureus.36171

Crossref Full Text | Google Scholar

Johansson, M. H. K., Bortolaia, V., Tansirichaiya, S., Aarestrup, F. M., Roberts, A. P., and Petersen, T. N. (2020). Detection of Mobile genetic elements associated with antibiotic resistance in Salmonella enterica using a newly developed web tool: Mobile element finder. J. Antimicrob. Chemother. 76, 101–109. doi: 10.1093/jac/dkaa390

PubMed Abstract | Crossref Full Text | Google Scholar

Loenen, W. A. M., Dryden, D. T. F., Raleigh, E. A., and Wilson, G. G. (2013). Type I restriction enzymes and their relatives. Nucleic Acids Res. 42, 20–44. doi: 10.1093/nar/gkt847

PubMed Abstract | Crossref Full Text | Google Scholar

Loenen, W. A., and Raleigh, E. A. (2014). The other face of restriction: modification-dependent enzymes. Nucleic Acids Res. 42, 56–69. doi: 10.1093/nar/gkt747

PubMed Abstract | Crossref Full Text | Google Scholar

Ma, J. Y., Zhao, H. H., Mo, S. Y., Li, J. J., Ma, X., Tang, Y. Q., et al. (2023). Acquisition of type I methyltransferase via horizontal gene transfer increases the drug resistance of Aeromonas veronii. Microb. Genom. 9:001107. doi: 10.1099/mgen.0.001107

PubMed Abstract | Crossref Full Text | Google Scholar

Oliveira, P. H., and Fang, G. (2021). Conserved DNA Methyltransferases: a window into fundamental mechanisms of epigenetic regulation in bacteria. Trends Microbiol. 29, 28–40. doi: 10.1016/j.tim.2020.04.007

PubMed Abstract | Crossref Full Text | Google Scholar

Olson, R. D., Assaf, R., Brettin, T., Conrad, N., Cucinell, C., Davis, J. J., et al. (2023). Introducing the bacterial and viral bioinformatics resource center (BV-BRC): a resource combining PATRIC, IRD and ViPR. Nucleic Acids Res. 51, D678–D689. doi: 10.1093/nar/gkac1003

PubMed Abstract | Crossref Full Text | Google Scholar

Page, A. J., Cummins, C. A., Hunt, M., Wong, V. K., Reuter, S., Holden, M. T. G., et al. (2015). Roary: rapid large-scale prokaryote pan genome analysis. Bioinformatics 31, 3691–3693. doi: 10.1093/bioinformatics/btv421

PubMed Abstract | Crossref Full Text | Google Scholar

Partridge, S. R., Kwong, S. M., Firth, N., and Jensen, S. O. (2018). Mobile genetic elements associated with antimicrobial resistance. Clin. Microbiol. Rev. 31:e00088-17. doi: 10.1128/cmr.00088-17

Crossref Full Text | Google Scholar

Pingoud, A., Wilson, G. G., and Wende, W. (2014). Type II restriction endonucleases—a historical perspective and more. Nucleic Acids Res. 42, 7489–7527. doi: 10.1093/nar/gku447

PubMed Abstract | Crossref Full Text | Google Scholar

Price, V. J., Huo, W., Sharifi, A., and Palmer, K. L. (2016). CRISPR-Cas and restriction-modification act additively against conjugative antibiotic resistance plasmid transfer in Enterococcus faecalis. mSphere 1:e00064-16. doi: 10.1128/mSphere.00064-16

Crossref Full Text | Google Scholar

Rao, D. N., Dryden, D. T. F., and Bheemanaik, S. (2013). Type III restriction-modification enzymes: a historical perspective. Nucleic Acids Res. 42, 45–55. doi: 10.1093/nar/gkt616

PubMed Abstract | Crossref Full Text | Google Scholar

Roberts, R. J., Vincze, T., Posfai, J., and Macelis, D. (2023). REBASE: a database for DNA restriction and modification: enzymes, genes and genomes. Nucleic Acids Res. 51, D629–d630. doi: 10.1093/nar/gkac975

PubMed Abstract | Crossref Full Text | Google Scholar

Seemann, T. (2014). Prokka: rapid prokaryotic genome annotation. Bioinformatics 30, 2068–2069. doi: 10.1093/bioinformatics/btu153

PubMed Abstract | Crossref Full Text | Google Scholar

Seong, H. J., Han, S. W., and Sul, W. J. (2021). Prokaryotic DNA methylation and its functional roles. J. Microbiol. 59, 242–248. doi: 10.1007/s12275-021-0674-y

PubMed Abstract | Crossref Full Text | Google Scholar

Shikov, A. E., Savina, I. A., Nizhnikov, A. A., and Antonets, K. S. (2023). Recombination in bacterial genomes: evolutionary trends. Toxins 15:568. doi: 10.3390/toxins15090568

PubMed Abstract | Crossref Full Text | Google Scholar

Stamatakis, A. (2014). RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 30, 1312–1313. doi: 10.1093/bioinformatics/btu033

PubMed Abstract | Crossref Full Text | Google Scholar

Takakura, W., Kudaravalli, P., Chatterjee, C., Pimentel, M., and Riddle, M. S. (2022). Campylobacter infection and the link with irritable bowel syndrome: on the pathway towards a causal association. Pathog. Dis. 80:ftac003. doi: 10.1093/femspd/ftac003

PubMed Abstract | Crossref Full Text | Google Scholar

Tang, B., Tang, Y. Z., Zhang, L., Liu, X., Chang, J., Xia, X. D., et al. (2020). Emergence of fex a in mediating resistance to Florfenicols in campylobacter. Antimicrob. Agents Chemother. 64:e00260-20. doi: 10.1128/aac.00260-20

PubMed Abstract | Crossref Full Text | Google Scholar

Tao, S., Chen, H. M., Li, N., Wang, T., and Liang, W. (2022). The spread of antibiotic resistance genes in vivo model. Can. J. Infect. Dis. Med. Microbiol. 2022:3348695. doi: 10.1155/2022/3348695

PubMed Abstract | Crossref Full Text | Google Scholar

Tesson, F., Hervé, A., Mordret, E., Touchon, M., d'Humières, C., Cury, J., et al. (2022). Systematic and quantitative view of the antiviral arsenal of prokaryotes. Nat. Commun. 13:2561. doi: 10.1038/s41467-022-30269-9

PubMed Abstract | Crossref Full Text | Google Scholar

van Vliet, A. H. M., Thakur, S., Prada, J. M., Mehat, J. W., and La Ragione, R. M. (2022). Genomic screening of antimicrobial resistance markers in UK and US campylobacter isolates highlights stability of resistance over an 18-year period. Antimicrob. Agents Chemother. 66:e0168721. doi: 10.1128/aac.01687-21

PubMed Abstract | Crossref Full Text | Google Scholar

Wei, H. L., Liao, Y. S., Chen, B. H., Teng, R. H., Wang, Y. W., Chang, J. H., et al. (2024). Antimicrobial resistance and genetic relatedness among campylobacter coli and Campylobacter jejuni from humans and retail chicken meat in Taiwan. J. Glob. Antimicrob. Resist. 38, 27–34. doi: 10.1016/j.jgar.2024.05.013

PubMed Abstract | Crossref Full Text | Google Scholar

Xie, J., Chen, Y., Cai, G., Cai, R., Hu, Z., and Wang, H. (2023). Tree visualization by one table (tvBOT): a web application for visualizing, modifying and annotating phylogenetic trees. Nucleic Acids Res. 51, W587–W592. doi: 10.1093/nar/gkad359

PubMed Abstract | Crossref Full Text | Google Scholar

Xu, C., Rao, J., Xie, Y., Lu, J., Li, Z., Dong, C., et al. (2023). The DNA Phosphorothioation restriction-modification system influences the antimicrobial resistance of pathogenic bacteria. Microbiol. Spectr. 11, e03509–e03522. doi: 10.1128/spectrum.03509-22

PubMed Abstract | Crossref Full Text | Google Scholar

Zeng, X., Wu, Z., Zhang, Q., Lin, J., and Elkins, C. A. (2018). A cotransformation method to identify a restriction-modification enzyme that reduces conjugation efficiency in Campylobacter jejuni. Appl. Environ. Microbiol. 84:e02004-18. doi: 10.1128/aem.02004-18

Crossref Full Text | Google Scholar

Keywords: Campylobacter jejuni , restriction-modification system, antibiotic resistance genes, transmission dynamics, drug resistance

Citation: Qiu Y, Guo P, Tian H, Zhou Y, Wen H and Liang H (2025) The restriction impacts of the Type III restriction-modification system on the transmission dynamics of antimicrobial resistance genes in Campylobacter jejuni. Front. Microbiol. 16:1496275. doi: 10.3389/fmicb.2025.1496275

Received: 14 September 2024; Accepted: 07 July 2025;
Published: 17 July 2025.

Edited by:

Bing Li, University of Science and Technology Beijing, China

Reviewed by:

Yuanyue Tang, Yangzhou University, China
Anand Bahadur Karki, Sam Houston State University, United States

Copyright © 2025 Qiu, Guo, Tian, Zhou, Wen and Liang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Hongling Wen, d2VuaG9uZ2xpbmdAc2R1LmVkdS5jbg==; Hao Liang, aGFvbGlhbmdAc2R1LmVkdS5jbg==

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.