v
Search
Advanced Search

Publications > Journals > Gene Expression > Article Full Text

  • OPEN ACCESS

Decoding the Complex Genetic Network of Antimicrobial Resistance in Campylobacter jejuni Using Advanced Gene Network Analysis

  • Prasanna Kumar Selvam1,
  • Santhosh Mudipalli Elavarasu1,
  • Hrituraj Dey1,
  • Karthick Vasudevan1,*  and
  • George Priya Doss2
 Author information
Gene Expression   2024;23(2):106-115

doi: 10.14218/GE.2023.00107

Abstract

Background and objectives

Antimicrobial resistance (AMR) poses a significant threat to public health in the 21st century, with bacteria such as Campylobacter jejuni (C. jejuni) exhibiting multidrug resistance due to the presence of AMR genes. Understanding the evolutionary patterns and functional relationships of these genes is crucial for addressing this issue effectively.

Methods

We conducted phylogenetic analysis to examine the evolution of AMR genes in C. jejuni. Additionally, we constructed and analyzed a gene interaction network comprising 39 functional relationships. Clustering analysis was employed to identify interconnected clusters associated with AMR processes. Functional enrichment analysis was performed to explore the involvement of cellular components, molecular functions, and biological processes.

Results

Our analysis revealed two interconnected clusters (C1 and C2) closely associated with AMR processes. Furthermore, genes encoding ribosomal proteins (rplE, rplV, rplG, rplK, rplA, rplJ, rpsE, rplB, rpsL, and rpmA) were identified as hub genes within the gene interaction network. These genes interact frequently with their functional counterparts, indicating their significance in AMR mechanisms. Enriched Kyoto Encyclopedia of Genes and Genomes pathway analysis highlighted the importance of the ribosome pathway in understanding antibiotic resistance mechanisms in C. jejuni.

Conclusions

The findings of this study enhance our understanding of the molecular mechanisms underlying AMR in C. jejuni. By elucidating the evolutionary patterns, gene interactions, and pathway enrichment, our study provides valuable insights that may contribute to the development of novel treatments for illnesses caused by this pathogen.

Graphical Abstract

Keywords

Antimicrobial resistance, Gene interaction network, Campylobacter jejuni, AMR, rpl, Cytoscape, ClueGo

Introduction

One of the greatest threats to public health in the 21st century is antimicrobial resistance (AMR) in bacteria. AMR occurs when bacteria undergo genetic changes that reduce the effectiveness of antibiotics used to treat infections. According to the UK Government-commissioned Study on Antimicrobial Resistance, AMR might result in the yearly death of 10 million people by 2050.1 The World Health Organization, as well as numerous other organizations and researchers, concur that the development of AMR is a pressing issue that must be addressed through a global, coordinated action plan.2,3

Understanding the full cost of resistance is a significant obstacle in the fight against AMR, especially in areas with scant surveillance and limited available information.1 Many studies have reported the impact of AMR for particular pathogen-drug combinations, particularly on incidence, fatalities, hospital length of stay, and healthcare costs.4–7 The occurrence of AMR and the proliferation of antibiotic-resistant microorganisms encompass a range of significant human diseases. The spread of AMR from hospital environments, which are generally closed communities, is seen as a threat to public health.

Campylobacter species are gram-negative, spiral-shaped, and nonspore-forming bacteria that thrive best in microaerophilic environments. The first Campylobacter may have been discovered as early as 1913; however, it was not until 1963 that the genus Campylobacter was formally recognized, having been categorized as Vibrio spp. Currently, the family Campylobacteriaceae includes the genera Campylobacter and Arcobacter.8Campylobacter comprises 14 Campylobacter species, and Campylobacter jejuni (C. jejuni) is frequently linked to human gastroenteritis.

Most human infections caused by Campylobacter usually involve C. jejuni, one of several species within the Campylobacter spp.9 In underdeveloped nations, Campylobacter infections are most commonly recorded in young people. Both children under one year of age and those under five in Southeast Asia have shown peak infection rates.10 Between 2.9% and 15% of children in Southeast Asia were previously found to have Campylobacter spp.11 According to the World Health Organization, AMR in Campylobacter spp. is a growing global concern and poses a significant public health threat.12 The Centers for Disease Control and Prevention estimates that each year in the United States, Campylobacter causes 1.3 million cases of human illness.13

Understanding the molecular factors underlying the AMR pattern is crucial to address this challenge and propose effective solutions. To accomplish this, we developed a gene interaction network to identify highly interacting genes.14 Our study highlights the critical importance of identifying hub genes within these networks to elucidate the mechanisms of AMR in C. jejuni. By focusing on these hub genes, researchers have gained valuable insights into the functional aspects of AMR, facilitating the identification of potential therapeutic targets. Furthermore, our approach highlights the significance of assessing biological pathways at the gene level, providing a dynamic framework for understanding AMR mechanisms comprehensively.15,16 This study not only contributes to the scientific understanding of AMR but also provides essential insights necessary for the development of effective treatments for infections caused by C. jejuni.

Materials and methods

NDARO and ABRicate

NDARO (National Database of Antibiotic Resistance Organisms) is an online repository run by the National Center for Biotechnology Information that is used to store information on the AMR genes of disease-causing bacteria. It was developed to provide background information and host-specific data on bacterial AMR genes. Each of the AMR databases in ABRicate v. 0.8 (https://github.com/tseemann/abricate/ ), including NCBI AMRFinderPlus, CARD, and ResFinder, contains data on thousands of AMRs.17Figure 1 visually outlines the sequential steps and methodology employed in the study, offering a concise overview of the research process.

Visual representation of the methodology involved in the study.
Fig. 1  Visual representation of the methodology involved in the study.

iTOL-a tool for phylogenetic tree development

iTOL is a web-based tool used for displaying, manipulating, and annotating phylogenetic trees. It facilitates interactive rooting and pruning of trees. The tool also enables the mapping of various data types onto trees, including genome sizes and protein domain repertoires. Additionally, iTOL supports exporting images in a variety of bitmap and vector graphic formats.18

The STRING database for network analysis

The STRING database (Search Tool for the Retrieval of Interaction Genes/Proteins) comprises several online platforms devoted to organism-wide protein interaction networks. The entire STRING database is precomputed, stored in a relational database, and independently accessible for download. Based on probabilistic confidence scores, the functional partners engaged in these interactions are described. STRING features a proprietary scoring method based on a reference set and several relationships.19

Cytoscape tool for network construction

The open-source network analysis and visualization software Cytoscape were used to construct and analyze networks of social and biological interactions. Cytoscape includes essential tools for network analysis, data integration, and visualization.20,21 By implementing Cytoscape plugins, additional functions can be introduced. Installing integrated database plugins provides access to data from other databases. In biological networks, each component is referred to as a node, and the connections between nodes are referred to as edges.22

MCODE tool for cluster analysis

Molecular Complex Detection (MCODE) is a Cytoscape tool that locates cluster regions with many interconnections in a network. It is a clustering technique that is reasonably quick.23,24 This approach is appropriate for academics who are focused on computation and biological research because of its simple interface. MCODE rates each cluster according to its size and density by giving it a score. Finally, the data are graphically represented as clusters.

NetworkAnalyzer

NetworkAnalyzer is a flexible and intuitive tool for examining biological and other networks. This plugin seamlessly integrates with Cytoscape and utilizes effective graph algorithms to compute a comprehensive list of simple and sophisticated topology parameters.25,26 It adds node properties for the outcomes and provides useful visualization settings to display and export the generated distributions.

ClueGO tool for gene ontology interpretation

ClueGO offers preconfigured functional analysis settings that range from broad to highly detailed. Additionally, Gene Ontology (GO) is a standardized system for annotating genes and their products across different species. It provides a controlled vocabulary to describe gene and gene product attributes in any organism, consisting of three structured networks: Biological Process, Molecular Function, and Cellular Component. Users can adjust analysis parameters to focus on terms within a specific range of GO, enabling a more precise exploration of gene functions and biological processes. GO levels have a specific evidence code, or have a specific number and percentage of linked genes. ClueGO first generates a binary gene-term matrix consisting of the selected terms and their related genes. Using a term-term similarity matrix built on this matrix and chance-corrected kappa statistics, the association strength between the words is calculated. Similar to BiNGO, ClueGO may be used in conjunction with GOlorize for the functional analysis of a Cytoscape gene network.27

Results

Data compilation from databases

For this study, 25 AMR genes associated with C. jejuni were extracted from the ABRicate and NDARO databases. To ensure the exclusion of duplicate entries, a process of elimination was implemented. Consequently, a set of 25 distinct AMR genes was identified for analysis.

Phylogenetic tree construction

The TREEFILE was aligned using MAFFT (Multiple Alignment using Fast Fourier Transform), a tool specifically designed for multiple sequence alignment. Following this alignment, iTOL, a program tailored for visualizing phylogenetic trees in a rooted format, was utilized to present the tree data (for details on network generation, refer to Appendix A in Supplementary file). Phylogenetic analysis was then carried out to examine the evolutionary connections among the C. jejuni strains.28Figure 2 displays the phylogenetic tree.

Schematic depiction of the rooted phylogenetic tree for the 25 genomic sequences for <italic>C. jejuni</italic> strains using the MAFFT (Multiple Alignment using Fast Fourier Transform).
Fig. 2  Schematic depiction of the rooted phylogenetic tree for the 25 genomic sequences for C. jejuni strains using the MAFFT (Multiple Alignment using Fast Fourier Transform).

STRING analysis of gene interactions

Seven specific genes (bla, fusA, gyrA, lepA, cjaA, sodB, and peb1A) associated with the 25 collected AMR genes were extracted from the STRING database. With a set minimum interaction score of 0.4 for medium confidence, our goal is to expand the network to include 39 nodes and 438 edges. This expansion aimed to maximize the representation of interactions among AMR genes, thereby enhancing our insight into their relationships. (For details on network generation, refer to Appendix B in supplementary file.)

Network analysis and hub gene identification

NetworkAnalyzer, a tool in Cytoscape, was used to conduct the analysis. The results showed a total number of 39 nodes and 438 edges in the network. A strong association between genes was revealed by the data summary, which comprised the highest degree count, the clustering factor, and the shortest path length at the minimum (Table 1). Among the genes, rplE, rplV, rpsG, rplK, rplA, rplJ, rpsE, rplB, rpsL, and rpmA were identified with the greatest number of connections (Table 1).

Table 1

List of the top 20 genes analyzed using NetworkAnalyzer considering various parameters, such as degree, average shortest path length, and betweenness centrality

Sl.GeneDegreeAverage shortest path lengthBetweenness centralityCloseness centralityClustering coefficient
1rplE321.2432432430.0213549840.8043478260.796370968
2rplV311.270270270.0125461760.7872340430.838709677
3rpsG301.2972972970.0142542610.7708333330.855172414
4rplK301.2972972970.006954870.7708333330.882758621
5rplA301.2972972970.0110562930.7708333330.857471264
6rplJ301.1891891890.248269940.8409090910.816091954
7rpsE301.2972972970.006954870.7708333330.882758621
8rplB301.2972972970.006954870.7708333330.882758621
9rpsL291.3243243240.0058403630.7551020410.903940887
10rpmA291.3243243240.0036705250.7551020410.918719212
11rpsA291.3243243240.0139221990.7551020410.849753695
12rpmF291.3243243240.0089132050.7551020410.901477833
13rpsF291.3243243240.0089132050.7551020410.901477833
14rpsD291.3243243240.0050153680.7551020410.913793103
15rpsP281.3513513510.0017310230.740.955026455
16rplS281.3513513510.0017310230.740.955026455
17rpsO281.3513513510.0017310230.740.955026455
18rpmB271.3783783780.0014954630.7254901960.96011396
19rpll271.3783783780.0006170.7254901960.985754986
20rplU271.3783783780.0006170.7254901960.985754986

Clustering analysis of gene networks

The application of the MCODE clustering technique led to the identification of two distinct groups, In Figure 3, labelled as C1 and C2, there are densely interconnected regions. Out of the total 39 genes analyzed, 33 were organized into these clusters, while the remaining 6 genes did not form part of any specific cluster (refer to Appendix C in supplementary file for a detailed description of the network generation process).29,30

Clustering analysis of the gene interaction network using MCODE tool resulted in two clusters, C1 (orange) and C2 (blue) where C1 had the highest level of clustering.
Fig. 3  Clustering analysis of the gene interaction network using MCODE tool resulted in two clusters, C1 (orange) and C2 (blue) where C1 had the highest level of clustering.

Nodes highlighted in yellow represents the no zero degree of the clustering of the genes. MCODE, Molecular Complex Detection.

Functional enrichment analysis of gene clusters

ClueGo was used to analyze overrepresented sequences and their functional components. The contributions of the numerous processes and pathways are clarified by this enrichment research. A moderate level of network specificity was used to evaluate the words and annotations from the Gene Ontology (GO) analysis. Moreover, pathway information was extracted from the Kyoto Encyclopedia of Genes and Genomes (KEGG) and UniProt databases using the STRING database, providing scientifically valuable pathway data for different genes. The analysis revealed that many pathways and activities, including biological processes (BPs), molecular functions (MFs), and cellular components (CCs), were associated with the AMR genes and their interactors. The network genes were shown to be enriched in 23 BPs, 21 CCs, and 16 MFs according to the enrichment analysis. The metabolic pathways involved in one KEGG pathway, 2 UniProt keywords along with gene expression (GO: 10467), cellular biosynthetic processes (GO: 44249), macromolecule biosynthetic processes (GO: 9059), organic substance biosynthetic processes (GO:1901576), and translation were among the highly enriched BPs (GO:6412) (Table 2).31Figure 4 depicts the biosynthetic process facilitated by genes like fusA, lepA, nusG, rplV, rplE, among others.

Table 2

Gene Ontology terms significantly enriched in biological process’s associated with genes

SIIDDescriptionGenes
1GO:10467gene expressionfusA, lepA, nusG, rplA, rplE, rplJ, rplK, rplL, rplS, rplT, rplU, rplV, rpmA, rpmB, rpmF, rpmG, rpsD, rpsE, rpsF, rpsG, rpsL, rpsO, rpsP.
2GO:44249cellular biosynthetic processfusA, lepA, nusG, rplA, rplE, rplJ, rplK, rplL, rplS, rplT, rplU, rplV, rpmA, rpmB, rpmF, rpmG, rpsD, rpsE, rpsF, rpsG, rpsL, rpsO, rpsP.
3GO:9059macromolecule biosynthetic processfusA, lepA, nusG, rplA, rplE, rplJ, rplK, rplL, rplS, rplT, rplU, rplV, rpmA, rpmB, rpmF, rpmG, rpsD, rpsE, rpsF, rpsG, rpsL, rpsO, rpsP.
4GO:1901576organic substance biosynthetic processfusA, lepA, nusG, rplA, rplE, rplJ, rplK, rplL, rplS, rplT, rplU, rplV, rpmA, rpmB, rpmF, rpmG, rpsD, rpsE, rpsF, rpsG, rpsL, rpsO, rpsP.
5GO:6412translationfusA, lepA, rplA, rplE, rplJ, rplK, rplL, rplS, rplT, rplU, rplV, rpmA, rpmB, rpmF, rpmG, rpsD, rpsE, rpsF, rpsG, rpsL, rpsO, rpsP.
Enrichment of the AMR (antimicrobial resistance) genes using ClueGO involved in BPs.
Fig. 4  Enrichment of the AMR (antimicrobial resistance) genes using ClueGO involved in BPs.

The abundant genes in distinct BPs are highlighted in yellow, while the resistance-related rpl genes are highlighted in red. BPs, biological processes.

Discussion

Understanding the whole expense of resistance poses a significant challenge in the battle against AMR, particularly in regions with limited monitoring and information accessibility. To address this requirement and offer a workable remedy, one must comprehend the molecular mechanisms underlying AMR. Our analysis revealed two phylogenetic clusters for the OXA genes, with the second cluster further subdivided into two subclusters. Additionally, the phylogenetic clusters of tetM and tetO appeared to be closely related. Tetracyclines, a group of drugs that encompasses tigecycline, minocycline, doxycycline, and tetracycline, are utilized for the management and treatment of various bacterial infections. The evolution of tetracycline resistance in Campylobacter, a gram-negative bacterium, is supported by the similarities between tetM and tetO.32 This finding reveals the evolution and family links of these gene variations. Genes within a network exhibiting the highest number of interactions are commonly known as hub genes. Due to their involvement in crucial BPs, these genes are considered essential.33 Understanding these genes allows us to assess the molecular mechanisms and processes that underlie an organism’s antibiotic resistance. In alignment with the results from the same study, the genes with the highest numbers of interactions included rplE, rplV, rpsG, rplK, rplA, rplJ, rpsE, rplB, rpsL, and rpmA (Table 1). The rpl gene is responsible for the assembly of the 60 s subunit of the eukaryotic ribosome, and plays a vital role in protein translation and cellular functions. Disruptions in ribosome assembly typically trigger a cellular stress response. Ribosomal proteins are essential for ribosome biogenesis and protein synthesis, and are crucial in diverse developmental processes. Furthermore, DNA replication, transcription, strand separation, repair, and DNA topoisomerase type II all require negative supercoiling, which is accomplished by the gyrA gene.14,34

Within the identified clusters, C1 emerged as the region with the highest connectivity, boasting 28 nodes and a notable score of 27.63. This suggests a robust interplay among the genes within C1, suggesting potential functional relationships. On the other hand, C2 comprised five nodes, each with a score of 5 (Table 3). This information provides insights into the structural organization of the network and highlights specific gene clusters that may play pivotal roles in antibiotic resistance.35,36 BPs provide essential knowledge about cellular repair and interactions with cells, facilitating the interaction between molecular machinery and catalytic processes in MFs. CCs, such as cell structures and complexes, are essential for DNA storage and detoxification. Bacteria must rewire their cellular metabolic pathways to survive in the host and respond to antibiotic exposure, necessitating ATP synthesis.37 Cellular metabolic processes are potential drug targets that are crucial for bacterial growth and survival.38 Several genes, including rpsL, are responsible for resistance to aminoglycosides, likely due to efflux pump mechanisms or other unidentified resistance mechanisms.39 The CmeABC multidrug resistance efflux pump is vital for C. jejuni colonization, regulating resistance to bile salts in the intestinal tract. The resistance-nodulation-division superfamily of bacterial transporters includes the tripartite efflux mechanism known as the multidrug resistance pump. The CmeABC complex allows C. jejuni to develop intrinsic resistance to a variety of antibiotics and other antimicrobial substances.40,41

Table 3

Clustering analysis of the gene interaction network

Number of nodesNumber of edgesCluster score = (density × no of nodes)Gene
Cluster 12837327.63rplJ, rpml, rpsD, rplV, rpmB, rpsA, rpsG, rpsF, rplT, rplY, fusA, rpsL, rpsO, rpmA, rpsE, rpmG, rplA, rpll, lepA, rplL, rplE, rplb, rplK, rpmF, rpsP, rplU, nusG, rplS.
Cluster 25105pebC, peb1A, Cj0919c, Cj0920c, cjaA.

The study also revealed that several ribosomal proteins are involved in the resistance mechanism. Clinical resistance to macrolides, lincosamides, streptogramins, and ketolides is related to modifications in ribosomal proteins L4 and L22 in various bacteria.42 Ribosomal proteins L4 and L22 interact with the CmeABC efflux pump to confer macrolide resistance. Understanding of the processes underlying Campylobacter resistance to macrolides is based on alterations in ribosomal proteins.43 The research also indicated that a few GO terms were associated with mutations in rpl genes linked to azithromycin resistance. Fusidic acid (FA) resistance develops through either horizontal acquisition of the fusA gene, which encodes a binding protein shielding the translation apparatus from FA inhibition, or genetic anomalies in the EF-G (fusA) gene.44 The antibiotic FA targets ribosome-bound EF-G in both translocation and ribosome recycling, thus inhibiting protein synthesis.45 FA has been used to treat infections caused by gram-positive bacteria since it was first identified in the early 1960s. However, the increase in bacterial pathogen resistance has become a growing clinical concern, with FA resistance often resulting from mutations in the EF-G-encoding fusA gene or the reduction/loss of ribosomal protein L6. To control both intrinsic termination and global gene expression, NusG and NusA collaborate. Loss of NusG leads to an altered pattern in fla/che operon expression, causing a reduced motility phenotype.46

In this study, a few of the enriched CC terms were related to intracellular organelles (GO: 43229), nonmembrane-bound organelles (GO: 43228), intracellular nonmembrane-bound organelles (GO: 43232), and ribosomes (GO: 5840) (Table 4). These terms were associated with a set of genes, Genes such as rplV, rplK, rplE, and rpsO are shown in Figure 5. One of the mechanisms by which bacteria can resist the effects of antibiotics is through drug efflux pumps, primarily found in gram-negative bacteria, which maintain the internal environment. Among the significant genes identified, rpsO encodes one of the main rRNA binding proteins that directly binds to 16S rRNA, aiding in the assembly of the 30S subunit platform. The rplV gene encodes the ribosomal 50S subunit protein L22, which is essential for early ribosomal 50S subunit synthesis.47,48 The nascent peptide exit tunnel of the growing ribosome needs to form.49rplK is a component of the ribosomal stalk, facilitating the interaction of the ribosome with guanosine Triphosphate-bound translation factors. The enrichment of molecular activities, including ribosome structural elements and nucleic acid binding, was also a significant discovery from the study (GO:3676, GO:3723, and GO: 19843) (Table 5). These functions are linked to a cluster of genes, including fusA, lepA, rplE, rplK, rplV, and rpsO, represented in Figure 6. Erythromycin, spectinomycin, and streptomycin resistance have all been linked to mutations in ribosomal proteins L22, S5, and S12.50 When bacteria are exposed to these medications, alterations in ribosomal proteins are integrated into the bacterium. In numerous bacterial species, including the intestinal pathogen C. jejuni, RNA-binding regulators have been shown to control posttranscriptional protein expression.51,52

Table 4

Significantly enriched Gene Ontology terms in cellular component’s with the associated genes

SIIDDescriptionGenes
1GO:43229intracellular organellerplA, rplE, rplJ, rplK, rplL, rplS, rplT, rplU, rplV, rpmA, rpmB, rpmF, rpmG, rpsA, rpsD, rpsE, rpsF, rpsG, rpsL, rpsO, rpsP.
2GO:43228non-membrane-bounded organellerplA, rplE, rplJ, rplK, rplL, rplS, rplT, rplU, rplV, rpmA, rpmB, rpmF, rpmG, rpsA, rpsD, rpsE, rpsF, rpsG, rpsL, rpsO, rpsP.
3GO:43232intracellular non-membrane-bounded organellerplA, rplE, rplJ, rplK, rplL, rplS, rplT, rplU, rplV, rpmA, rpmB, rpmF, rpmG, rpsA, rpsD, rpsE, rpsF, rpsG, rpsL, rpsO, rpsP.
4GO:5840ribosomerplA, rplE, rplJ, rplK, rplL, rplS, rplT, rplU, rplV, rpmA, rpmB, rpmF, rpmG, rpsA, rpsD, rpsE, rpsF, rpsG, rpsL, rpsO, rpsP.
Enrichment of the AMR (antimicrobial resistance) genes using ClueGO involved in CCs.
Fig. 5  Enrichment of the AMR (antimicrobial resistance) genes using ClueGO involved in CCs.

The genes enriched in CCs are highlighted in yellow, whereas the top resistance genes are shown in red. CCs, cellular components.

Table 5

Gene Ontology terms significantly enriched in molecular function’s associated with genes

SIIDDescriptionGenes
1GO:3676nucleic acid bindingfusA, lepA, rplA, rplE, rplJ, rplK, rplT, rplU, rplV, rpsA, rpsD, rpsE, rpsF, rpsG, rpsL, rpsO.
2GO:3723RNA bindingfusA, lepA, rplA, rplE, rplJ, rplK, rplT, rplU, rplV, rpsD, rpsE, rpsF, rpsG, rpsL, rpsO.
3GO:19843rRNA bindingrplA, rplE, rplJ, rplK, rplT, rplU, rplV, rpsD, rpsE, rpsF, rpsG, rpsL, rpsO.
Enrichment of the AMR (antimicrobial resistance) genes using ClueGO involved in MFs.
Fig. 6  Enrichment of the AMR (antimicrobial resistance) genes using ClueGO involved in MFs.

The genes enriched in MFs were highlighted in yellow, whereas the top resistance genes are shown in red. MFs, molecular functions.

Transformations in certain ribosomal proteins are known to be connected with antibiotic resistance, affecting the stability or translation of the mRNAs they bind to.53 According to several studies, antibiotic exposure alters the metabolic state of bacteria. The activation of efflux pumps alters bacterial susceptibility to antibiotics, leading to the development of resistance mechanisms and influencing biofilm formation.54,55 As a result, the enrichment of the ribosome pathway in the KEGG pathway is essential for understanding antibiotic resistance in bacteria.56

Through our studies, we have gained a better understanding of the intricate functional relationships between genes and their variations. Our results may have significant implications for the development of innovative therapies and diagnostic tools for C. jejuni infections because they illuminate the complex interplay of genetic variables in a range of BPs.

Conclusions

C. jejuni stands out as a predominant pathogen in global food-borne outbreaks, notably amid increasing concerns about AMR. A recent study focused on tetracycline resistance genes tetO and tetM. By employing phylogenetic tree analysis, this research has provided valuable insights into the genetic landscape and variants associated with C. jejuni. The investigation highlighted the key hub genes such as rplE, rplV, rplG, and others, revealing their integral roles in AMR through GO keywords such as gene expression, cellular biosynthetic processes, and RNA binding. Crucially, this study highlighted the significance of the rpl gene in driving the AMR phase of C. jejuni. These hub genes, exhibiting a high degree of clustering with their functional partners, have emerged as potential drug targets. The study’s findings raise hope that targeting these genes could pave the way for innovative treatments combating AMR in C. jejuni infections. This comprehensive exploration of genetic and functional aspects offers valuable insights into the complex dynamics of AMR, providing a foundation for future therapeutic interventions and strategies in the ongoing battle against antibiotic resistance in C. jejuni.

Abbreviations

AMR: 

antimicrobial resistance

BPs: 

biological processes

C. jejuni

Campylobacter jejuni

CCs: 

cellular components

FA: 

Fusidic acid

GO: 

gene ontology

KEGG: 

Kyoto Encyclopedia of Genes and Genomes

MCODE: 

Molecular Complex Detection

MFs: 

molecular functions

STRING: 

Search Tool for the Retrieval of Interaction Genes/Proteins

Declarations

Acknowledgement

The authors express deep gratitude to the management of REVA University and the Vellore Institute of Technology for providing all the support, necessary facilities, assistance, and constant encouragement to carry out this work.

Data sharing statement

The data will be made available upon request.

Funding

No funding was received for this work.

Conflict of interest

The authors declare no conflict of interests.

Authors’ contributions

Writing-review & editing (PKS, SME, KV, GPDC), Writing-original draft (PKS), Visualization (PKS), Validation (PKS), Methodology (PKS), Formal analysis (PKS, SME, HD), Data curation (PKS), Conceptualization (PKS, SME, KV), Project administration (KV), Supervision (GPDC).

References

  1. Antimicrobial Resistance Collaborators. Global burden of bacterial antimicrobial resistance in 2019: A systematic analysis. Lancet 2022;399(10325):629-655 View Article PubMed/NCBI
  2. de Kraker ME, Stewardson AJ, Harbarth S. Will 10 Million People Die a Year due to Antimicrobial Resistance by 2050?. PLoS Med 2016;13(11):e1002184 View Article PubMed/NCBI
  3. Kadri SS. Key Takeaways From the U.S. CDC’s 2019 Antibiotic Resistance Threats Report for Frontline Providers. Critical care medicine 2020;48(7):939-945 View Article
  4. Naylor NR, Atun R, Zhu N, Kulasabanathan K, Silva S, Chatterjee A, et al. Estimating the burden of antimicrobial resistance: a systematic literature review. Antimicrob Resist Infect Control 2018;7:58 View Article PubMed/NCBI
  5. Cassini A, Högberg LD, Plachouras D, Quattrocchi A, Hoxha A, Simonsen GS, et al. Attributable deaths and disability-adjusted life-years caused by infections with antibiotic-resistant bacteria in the EU and the European Economic Area in 2015: a population-level modelling analysis. Lancet Infect Dis 2019;19(1):56-66 View Article PubMed/NCBI
  6. Lim C, Takahashi E, Hongsuwan M, Wuthiekanun V, Thamlikitkul V, Hinjoy S, et al. Epidemiology and burden of multidrug-resistant bacterial infection in a developing country. Elife 2016;5:e18082 View Article PubMed/NCBI
  7. Temkin E, Fallach N, Almagor J, Gladstone BP, Tacconelli E, Carmeli Y. Estimating the number of infections caused by antibiotic-resistant Escherichia coli and Klebsiella pneumoniae in 2014: a modelling study. Lancet Glob Health 2018;6(9):e969-e979 View Article PubMed/NCBI
  8. Debruyne L, Gevers D, Vandamme P. Campylobacter. New Jersey: Wiley; 2008 View Article
  9. Tay ST, Puthucheary SD, Devi S, Kautner I. Characterisation of Campylobacters from Malaysia. Singapore Med J 1995;36(3):282-284 PubMed/NCBI
  10. Friedman CR, Neimann J, Wegener HC, Tauxe RV. Epidemiology of Campylobacter jejuni infections in the United States and other industrialized nations. Novelty (OH): ASM International; 2000
  11. Padungton P, Kaneene JB. Campylobacter spp in human, chickens, pigs and their antimicrobial resistance. J Vet Med Sci 2003;65(2):161-70 View Article PubMed/NCBI
  12. Moore JE, Barton MD, Blair IS, Corcoran D, Dooley JSG, Fanning S, et al. The epidemiology of antibiotic resistance in Campylobacter. Microbes Infect 2006;8(7):1955-1966 View Article PubMed/NCBI
  13. Fischer GH, Hashmi MF, Paterek E. StatPearls. StatPearls Publishing; 2024
  14. Dey H, Vasudevan K, Dasegowda KR, Rambabu M, Cn P, C GPD. An integrated gene network analysis to decode the multi-drug resistance mechanism in Klebsiella pneumoniae. Microb Pathog 2022;173(Pt A):105878 View Article PubMed/NCBI
  15. Anitha P, Anbarasu A, Ramaiah S. Gene network analysis reveals the association of important functional partners involved in antibiotic resistance: A report on an important pathogenic bacterium Staphylococcus aureus. Gene 2016;575(2 Pt 1):253-263 View Article PubMed/NCBI
  16. Ashok G, Miryala SK, Anbarasu A, Ramaiah S. Integrated systems biology approach using gene network analysis to identify the important pathways and new potential drug targets for Neuroblastoma. Gene Rep 2021;23:1101 View Article
  17. Xihui Z, Yanlan L, Zhiwei W, Zheyu P, Zhenshu S, Cheng L, et al. Antibiotic resistance of Riemerella anatipestifer and comparative analysis of antibiotic-resistance gene detection methods. Poult Sci 2023;102(3):102405 View Article PubMed/NCBI
  18. Letunic I, Bork P. Interactive Tree Of Life (iTOL) v5: an online tool for phylogenetic tree display and annotation. Nucleic Acids Res 2021;49(W1):W293-W296 View Article PubMed/NCBI
  19. Naha A, Kumar Miryala S, Debroy R, Ramaiah S, Anbarasu A. Elucidating the multi-drug resistance mechanism of Enterococcus faecalis V583: A gene interaction network analysis. Gene 2020;748:4704 View Article PubMed/NCBI
  20. Otasek D, Morris JH, Bouças J, Pico AR, Demchak B. Cytoscape Automation: empowering workflow-based network analysis. Genome Biol 2019;20(1):185 View Article PubMed/NCBI
  21. Mousavian Z, Khodabandeh M, Sharifi-Zarchi A, Nadafian A, Mahmoudi A. StrongestPath: a Cytoscape application for protein-protein interaction analysis. BMC Bioinformatics 2021;22(1):352 View Article PubMed/NCBI
  22. Miryala SK, Anbarasu A, Ramaiah S. Gene interaction network approach to elucidate the multidrug resistance mechanisms in the pathogenic bacterial strain Proteus mirabilis. J Cell Physiol 2021;236(1):468-479 View Article PubMed/NCBI
  23. Menon SMP, Elengoe A. Evaluation of the role of kras gene in colon cancer pathway using string and Cytoscape software. Biomed Res Ther 2020;7(6):3835-3842 View Article
  24. Majeed A, Mukhtar S. Protein-Protein Interaction Network Exploration Using Cytoscape. Methods Mol Biol 2023;2690:419-427 View Article PubMed/NCBI
  25. Yang Y, Xu X. Bioinformatic identification of hub genes and related transcription factors in low shear stress treated endothelial cells. BMC Med Genomics 2021;14(1):120 View Article PubMed/NCBI
  26. Díaz–Montaña JJ, Díaz–Díaz N, Barranco CD, Ponzoni I. Development and use of a Cytoscape app for GRNCOP2. Comput Methods Programs Biomed 2019;177:211-218 View Article PubMed/NCBI
  27. Bindea G, Mlecnik B, Hackl H, Charoentong P, Tosolini M, Kirilovsky A, et al. ClueGO: a Cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks. Bioinformatics 2009;25(8):1091-1093 View Article PubMed/NCBI
  28. Zhou T, Xu K, Zhao F, Liu W, Li L, Hua Z, et al. itol.toolkit accelerates working with iTOL (Interactive Tree of Life) by an automated generation of annotation files. Bioinformatics 2023;39(6):btad339 View Article PubMed/NCBI
  29. Abdelfattah EM, Ekong PS, Okello E, Chamchoy T, Karle BM, Black RA, et al. Epidemiology of antimicrobial resistance (AMR) on California dairies: descriptive and cluster analyses of AMR phenotype of fecal commensal bacteria isolated from adult cows. PeerJ 2021;9:e11108 View Article PubMed/NCBI
  30. Suwono B, Eckmanns T, Kaspar H, Merle R, Zacher B, Kollas C, et al. Cluster analysis of resistance combinations in Escherichia coli from different human and animal populations in Germany 2014-2017. PLoS One 2021;16(1):e0244413 View Article PubMed/NCBI
  31. Udaondo Z, Matilla MA. Mining for novel antibiotics in the age of antimicrobial resistance. Microb Biotechnol 2020;13(6):1702-1704 View Article PubMed/NCBI
  32. Sougakoff W, Papadopoulou B, Nordmann P, Courvalin P. Nucleotide Sequence and Distribution of Gene TetO Encoding Tetracycline Resistance in Campylobacter Coli. FEMS Microbiol Let 1987;44:153-159 View Article
  33. Anusha M, Tejaswini V, Udhaya Kumar S, Prashantha CN, Vasudevan K, George Priya Doss C. Gene network interaction analysis to elucidate the antimicrobial resistance mechanisms in the Clostridium difficile. Microb Pathog 2023;178:106083 View Article
  34. Miryala SK, Anbarasu A, Ramaiah S. Role of SHV-11, a Class A β-Lactamase, Gene in Multidrug Resistance Among Klebsiella pneumoniae Strains and Understanding Its Mechanism by Gene Network Analysis. Microbial Drug Resistance 2020;26(8):900-908 View Article PubMed/NCBI
  35. Feldgarden M, Brover V, Haft DH, Prasad AB, Slotta DJ, Tolstoy I, et al. Validating the AMRFinder Tool and Resistance Gene Database by Using Antimicrobial Resistance Genotype-Phenotype Correlations in a Collection of Isolates. Antimicrob Agents Chemother 2019;63(11):e00483-19 View Article PubMed/NCBI
  36. Liu Y, Cheng D, Xue J, Weaver L, Wakelin SA, Feng Y, et al. Changes in microbial community structure during pig manure composting and its relationship to the fate of antibiotics and antibiotic resistance genes. J Hazard Mater 2020;389:122082 View Article PubMed/NCBI
  37. Baquero F, Martínez JL, F Lanza V, Rodríguez-Beltrán J, Galán JC, San Millán A, et al. Evolutionary Pathways and Trajectories in Antibiotic Resistance. Clin Microbiol Rev 2021;34(4):e0005019 View Article PubMed/NCBI
  38. Moo CL, Yang SK, Yusoff K, Ajat M, Thomas W, Abushelaibi A, et al. Mechanisms of Antimicrobial Resistance (AMR) and Alternative Approaches to Overcome AMR. Curr Drug Discov Technol 2020;17(4):430-447 View Article PubMed/NCBI
  39. Marotta F, Di Marcantonio L, Janowicz A, Pedonese F, Di Donato G, Ardelean A, et al. Genotyping and antibiotic resistance traits in campylobacter jejuni and coli from pigs and wild boars in italy. Front Cell Infect Microbiol 2020;10:592512 View Article PubMed/NCBI
  40. Vieira A, Ramesh A, Seddon AM, Karlyshev AV. CmeABC multidrug efflux pump contributes to antibiotic resistance and promotes campylobacter jejuni survival and multiplication in acanthamoeba polyphaga. Appl Environ Microbiol 2017;83(22):e01600-17 View Article PubMed/NCBI
  41. Lin J, Michel LO, Zhang Q. CmeABC functions as a multidrug efflux system in Campylobacter jejuni. Antimicrob Agents Chemother 2002;46(7):2124-2131 View Article PubMed/NCBI
  42. Roberts MC. Resistance to macrolide, lincosamide, streptogramin, ketolide, and oxazolidinone antibiotics. Mol Biotechnol 2004;28(1):47-62 View Article PubMed/NCBI
  43. Cagliero C, Mouline C, Cloeckaert A, Payot S. Synergy between efflux pump CmeABC and modifications in ribosomal proteins L4 and L22 in conferring macrolide resistance in Campylobacter jejuni and Campylobacter coli. Antimicrob Agents Chemother 2006;50(11):3893-3896 View Article PubMed/NCBI
  44. NOAH. NOAH response to final O Neill AMR review. London: NOAH; 2016
  45. Borg A, Pavlov M, Ehrenberg M. Mechanism of fusidic acid inhibition of RRF- and EF-G-dependent splitting of the bacterial post-termination ribosome. Nucleic Acids Res 2016;44(7):3264-3275 View Article PubMed/NCBI
  46. Mandell ZF, Oshiro RT, Yakhnin AV, Vishwakarma R, Kashlev M, Kearns DB, et al. NusG is an intrinsic transcription termination factor that stimulates motility and coordinates gene expression with NusA. Elife 2021;10:e61880 View Article PubMed/NCBI
  47. Koutsoumanis K, Allende A, Álvarez-Ordóñez A, Bolton D, Bover-Cid S, Chemaly M, et al. Role played by the environment in the emergence and spread of antimicrobial resistance (AMR) through the food chain. EFSA J 2021;19(6):e06651 View Article PubMed/NCBI
  48. Hiller CX, Hübner U, Fajnorova S, Schwartz T, Drewes JE. Antibiotic microbial resistance (AMR) removal efficiencies by conventional and advanced wastewater treatment processes: A review. Sci Total Environ 2019;685:596-608 View Article PubMed/NCBI
  49. Han D, Liu Y, Li J, Liu C, Gao Y, Feng J, et al. Twenty-seven-nucleotide repeat insertion in the rplV gene confers specific resistance to macrolide antibiotics in Staphylococcus aureus. Oncotarget 2018;9(40):26086-26095 View Article PubMed/NCBI
  50. Wilcox SK, Cavey GS, Pearson JD. Single ribosomal protein mutations in antibiotic-resistant bacteria analyzed by mass spectrometry. Antimicrob Agents Chemother 2001;45(11):3046-3055 View Article PubMed/NCBI
  51. El Abbar FM, Li J, Owen HC, Daugherty CL, Fulmer CA, Bogacz M, et al. RNA Binding by the Campylobacter jejuni Post-transcriptional Regulator CsrA. Front Microbiol 2019;10:1776 View Article PubMed/NCBI
  52. Bairán G, Rebollar-Pérez G, Chávez-Bravo E, Torres E. Treatment Processes for Microbial Resistance Mitigation: The Technological Contribution to Tackle the Problem of Antibiotic Resistance. Int J Environ Res Public Health 2020;17(23):8866 View Article PubMed/NCBI
  53. Han D, Liu Y, Li J, Liu C, Gao Y, Feng J, et al. Twenty-seven-nucleotide repeat insertion in the rplV gene confers specific resistance to macrolide antibiotics in Staphylococcus aureus. Oncotarget 2018;9(40):26086-26095 View Article PubMed/NCBI
  54. Soto SM. Role of efflux pumps in the antibiotic resistance of bacteria embedded in a biofilm. Virulence 2013;4(3):223-229 View Article PubMed/NCBI
  55. Kumar S, Varela MF. Biochemistry of bacterial multidrug efflux pumps. Int J Mol Sci 2012;13(4):4484-4495 View Article PubMed/NCBI
  56. Anusha M, Tejaswini V, Udhaya Kumar S, Prashantha CN, Vasudevan K, George Priya Doss C. Gene network interaction analysis to elucidate the antimicrobial resistance mechanisms in the Clostridiumdifficile. Microb Pathog 2023;178:106083 View Article PubMed/NCBI
  • Gene Expression
  • pISSN 1052-2166
  • eISSN 1555-3884
Back to Top

Decoding the Complex Genetic Network of Antimicrobial Resistance in Campylobacter jejuni Using Advanced Gene Network Analysis

Prasanna Kumar Selvam, Santhosh Mudipalli Elavarasu, Hrituraj Dey, Karthick Vasudevan, George Priya Doss
  • Reset Zoom
  • Download TIFF