Next Article in Journal
The Significance of Exosomal RNAs in the Development, Diagnosis, and Treatment of Gastric Cancer
Next Article in Special Issue
Genomic Investigation of Carbapenem-Resistant Klebsiella pneumonia Colonization in an Intensive Care Unit in South Africa
Previous Article in Journal
Comparative Genomics Provides Insights into the Taxonomy of Azoarcus and Reveals Separate Origins of Nif Genes in the Proposed Azoarcus and Aromatoleum Genera
Previous Article in Special Issue
Long-Term Impact of Suppressive Antibiotic Therapy on Intestinal Microbiota
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

In Silico Characterization of Toxin-Antitoxin Systems in Campylobacter Isolates Recovered from Food Sources and Sporadic Human Illness

by
Bishoy Wadie
1,
Mohamed A. Abdel-Fattah
2,
Alshymaa Yousef
1,
Shaimaa F. Mouftah
1,
Mohamed Elhadidy
1,3,* and
Tamer Z. Salem
1,4,*
1
Biomedical Sciences Program, University of Science and Technology, Zewail City of Science and Technology, October Gardens, 6th of October City, Giza 12578, Egypt
2
Department of Agricultural Biotechnology, Faculty of Agriculture, Ain Shams University, Cairo 11566, Egypt
3
Department of Bacteriology, Mycology, and Immunology, Faculty of Veterinary Medicine, Mansoura University, Mansoura 35516, Egypt
4
Department of Microbial Genetics, AGERI, ARC, Giza 12619, Egypt
*
Authors to whom correspondence should be addressed.
Genes 2021, 12(1), 72; https://doi.org/10.3390/genes12010072
Submission received: 30 October 2020 / Revised: 30 November 2020 / Accepted: 1 December 2020 / Published: 7 January 2021
(This article belongs to the Special Issue Genetics of Antimicrobial Resistance)

Abstract

:
Campylobacter spp. represents the most common cause of gastroenteritis worldwide with the potential to cause serious sequelae. The ability of Campylobacter to survive stressful environmental conditions has been directly linked with food-borne illness. Toxin-antitoxin (TA) modules play an important role as defense systems against antimicrobial agents and are considered an invaluable strategy harnessed by bacterial pathogens to survive in stressful environments. Although TA modules have been extensively studied in model organisms such as Escherichia coli K12, the TA landscape in Campylobacter remains largely unexplored. Therefore, in this study, a comprehensive in silico screen of 111 Campylobacter (90 C. jejuni and 21 C. coli) isolates recovered from different food and clinical sources was performed. We identified 10 type II TA systems belonging to four TA families predicted in Campylobacter genomes. Furthermore, there was a significant association between the clonal population structure and distribution of TA modules; more specifically, most (12/13) of the Campylobacter isolates belonging to ST-21 isolates possess HicB-HicA TA modules. Finally, we observed a high degree of shared synteny among isolates bearing certain TA systems or even coexisting pairs of TA systems. Collectively, these findings provide useful insights about the distribution of TA modules in a heterogeneous pool of Campylobacter isolates from different sources, thus developing a better understanding regarding the mechanisms by which these pathogens survive stressful environmental conditions, which will further aid in the future designing of more targeted antimicrobials.

1. Introduction

Campylobacter is considered the most common bacterial pathogen responsible for human gastroenteritis worldwide [1]. Food-producing animals represent the main reservoir for Campylobacter infection, and transmission was previously reported to occur through the consumption of food of animal origins, including undercooked chicken, unpasteurized/incomplete milk, and dairy products, and through contaminated drinking water [2].
Several studies have focused on the identification and characterization of virulence factors in Campylobacter through which infection is mediated. These factors are primarily associated with the expression of genes involved in colonization, cell invasion, motility, and toxin production [3,4]. Such investigations have provided the research community with a better understanding to identify potential drug targets and to develop therapeutic interventions. However, the prevalence of cases with campylobacteriosis is still on the rise, albeit with the availability of a wide variety of antimicrobial agents. Despite the considerable importance of dissecting the mechanisms by which the virulence factors drive the infectivity of Campylobacter, our understanding of how Campylobacter is defending itself against antimicrobials is still lacking.
Bacterial cells are continuously exposed to exogenous genetic elements by means of horizontal gene transfer and phage infections that put them at risk of cell death. Generally, plasmids are transferred to other bacterial cells either by transformation or conjugation, and while it is commonly known that plasmids confer a selective advantage to the recipient bacterial cells, such as antibiotic resistance, other plasmids (e.g., Col plasmids) might contain bacteriocins that are proteins that kill other bacterial cells [5]. On the other hand, phages are well-known for their ability to infect and eventually lyse the bacterial cells. To be protected against such hazards, bacterial cells utilize three general defense strategies: (1) preventing pathogen entry, (2) adopting abortive infection (suicide/dormancy when infected), and (3) utilizing cellular immunity, which encompasses innate mechanisms (such as restriction modification systems) and adaptive mechanisms (such as CRISPR-Cas systems) [6]. In response to such defense systems, phages have, in turn, developed ways to counteract these systems, leading to an evolutionary arms race influenced by the extensive coevolution of both phage and host [7]. Since some defense systems are antagonized by phages, bacterial genomes encode multiple defense systems in discrete genomic loci called “defense islands”; these sites are frequently found in sites undergoing horizontal gene transfer, which explains the rapid diversification of these systems [6]. One of these defense systems is the toxin-antitoxin (TA) system, which is the primary focus of this study.
Toxin-antitoxin (TA) modules are comprised of a stable toxin that inhibits cellular growth and an unstable antitoxin that counteracts its action. TA modules can be generally classified into four major types (types I, II, III, and IV) based on the mechanism through which the antitoxin neutralizes the activity of the toxin protein, in addition to two minor types (types V and VI). Type I and type III antitoxins are RNA molecules that regulate the activity of their toxin counterparts through inhibiting the translation of toxin mRNA or directly inhibiting the toxin protein, respectively. Type II antitoxins are proteins that directly bind and inhibit toxin proteins, unlike type IV antitoxin proteins that neutralize the toxin’s activity without direct interaction [8]. Unlike types I–IV, types V and VI antitoxins contribute to the degradation of their toxin counterparts, either by acting as an RNase leading to the cleavage of toxin mRNA or forming a complex with a protease to cleave the toxin protein, respectively [9,10]. TA systems were first discovered on plasmids, and since then, they have been labeled as mobile genetic elements that are frequently transferred horizontally [8,11]. All types of TA modules are also enriched on bacterial chromosomes; however, the composition and quantities of chromosomal TA modules are highly variable among different bacterial species and even between strains [12,13]. However, no clear distinctions can be drawn between TA modules and the biology of a given species. Nonetheless, these TA systems are more frequently found in organisms living in hostile environments with recurrent events of horizontal gene transfer [12,13]. Toxins have been shown to utilize a myriad of molecular mechanisms to inhibit cell growth, most of which are well-characterized [8,14,15]. On the contrary, the biological functions of most TA modules remain elusive. Some functions have been proposed, however; the three most common ones include postsegregational killing (PSK) aimed at the stabilization of mobile genetic elements, abortive infection through altruistic suicide, and dormancy through which bacterial cells remains dormant and become antibiotic-tolerant [8,16].
As previously mentioned, the composition and quantities of TA modules are highly variable among organisms. Most studies are biased towards the characterization of TA modules in model organisms such as Escherichia coli K12. Nonetheless, little is known about the composition and numbers of TA modules in Campylobacter. In this study, we performed a comprehensive in silico screen for TA modules on 111 Campylobacter isolates sampled from different sources and assessed the distribution and conservation of the predicted TA modules among those isolates.

2. Materials and Methods

2.1. Bacterial Isolates and Culture Conditions

A total of 111 Campylobacter isolates (90 C. jejuni and 21 C. coli), collected between 2017 and 2018 in Cairo, Egypt, were isolated from foods of animal origin (commercial broiler carcasses (n = 30), milk and unpasteurized/incompletely pasteurized dairy products (n = 24)) and fecal clinical samples from patients suffering from diarrheal illness (n = 57) admitted to two different hospitals in Cairo, Egypt (Supplementary File S1). A stratified randomized sampling was conducted to collect food samples from different retail stores located around the study region. Isolation and enumeration of Campylobacter isolates from various food matrices was achieved according to the ISO 10272-1 (enrichment method; detection of Campylobacter spp. after selective enrichment). Stool samples from human patients were grown on modified charcoal cefoperazone deoxycholate agar (mCCDA). Plates were then incubated for 48 h at 42 °C under anaerobic conditions using AnaeroGen™ 2.5-L sachets (Oxoid, Basingstoke, UK). Confirmation of Campylobacter at the genus level was performed using PCR identification of the 16S rRNA gene [17]. Confirmation of C. jejuni and C. coli at the species level was performed by PCR detection of the mapA gene and ceuE gene [18,19], respectively. All isolates were sub-cultured from −80 °C frozen stocks onto Mueller-Hinton (MH) agar (Oxoid, Basingstoke, UK).

2.2. Multilocus Sequence Typing (MLST)

Sequence types (ST) and clonal complexes (CC) of all isolates were determined from Whole Genome Sequencing (WGS) data and queried against the sequences in the BIGSdb database [20] which automatically identifies alleles and assigns sequence types (ST) and a clonal complex (CC) for each isolate sequence.

2.3. Data Acquisition

Campylobacter genome assemblies (FASTA) for each corresponding isolate record were retrieved from the BIGSdb database (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA576513). BIGSdb is an open-source web-based platform that integrates phenotype and sequence data for a multitude of bacterial species building on multilocus sequence typing (MLST) data [20].

2.4. Genome Annotation

Since the bacterial species for all the isolates identified are either C. jejuni or C. coli, genome annotation was performed against different strains for the above species. Contigs from BIGSdb (FASTA files) for each isolate were annotated using prokka v1.14.5 against a custom database [21]. The custom database was built using annotated genomes for different C. jejuni and C. coli based on the TADB (Toxin-Antitoxin database) and Microbesonline list of available strains [22,23]. A list of Campylobacter reference genomes can be found in Supplementary File S2. Annotated genomes (GBK files) were downloaded from Microbesonline (http://www.microbesonline.org/) and NCBI (https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3570591/) and were converted to protein format (FAA files) using in-house python scripts. Then, all protein sequences were concatenated, clustered, and added to prokka’s genus directory to be used as a reference for annotation. The remaining parameters were set to default. For each isolate sequence, GFF and FNA output files were stored and used for downstream analyses.

2.5. In Silico Screening for TA Systems

To identify toxin-antitoxin genes in each isolate, TADB and PADS Arsenal (prokaryotic defense system-related genes) databases were used, which contain both predicted and experimentally validated TA sequences. TADB provides ample data about bacterial type II TA systems covering 6193 loci, 105 of which are experimentally validated [22]. PADS Arsenal is a public database integrating 18 different categories of defense systems (including TA systems), with a total of 6,600,624 annotated genes from 63,701 genomes covering 33,390 species of bacteria and archaea [24]. FASTA files of toxin and antitoxin genes were downloaded from TADB and PADS Arsenal databases. FNA output files for each isolate were aligned against TA genes from both TADB and PADS Arsenal using the BLSATN command line tool to identify TA genes present in each isolate [25]. All csv outputs for each database were merged and records having e-values > 10−5 were filtered out. In addition, only reference genome assemblies with a level of assembly labeled as “Complete” or “Chromosome” were included, as opposed to “Contig” or “Scaffold” from PADS Arsenal hits. Finally, in order to account for gene coverage, the percentage identity was multiplied by the alignment length followed by scaling values from 0 to 1 using “min–max normalization”. The scaled score represents the quality of the alignment given that multiple e-values were equal to zero because of the relatively narrow search space. It is worth noting that only four out of 111 Campylobacter isolates harbored type I antitoxin genes with no toxin counterparts, 3 of which are Campylobacter jejuni (7672, 7673, and 7690) and 1 Campylobacter coli (7700). Interestingly, these genes represent experimentally validated antitoxin genes according to TADB. An annotated file for all TA modules from TADB and PADS Arsenal can be found in Supplementary File S3.

2.6. Protein Domain Analysis

The protein domain for each toxin-antitoxin gene was identified using the NCBI conserved domain database (CDD) [26] and InterProScan web-based search tool [27]. These tools were used to confirm that the identified TA genes bear a TA-related functional domain.

2.7. Sequence Alignment and Clustering

In order to evaluate the sequence conservation of TA-related domains among hits from different genome assemblies, multiple sequence alignment was implemented for each locus tag possessing a TA-related domain. Sequence alignment was performed using the msa R package using default parameters [28]. In addition, hierarchical sequence clustering was performed using the seqinr R package [29] for all locus tags regardless of having a domain for visualization purposes.

2.8. Synteny Analysis

Syntenic regions between isolate genomes were identified to better assess the degree of conservation of homologous genes and gene order among isolates containing TA systems. Accordingly, for each isolate, only the contigs containing at least one TA system were selected for synteny analysis. Matching syntenic regions between two contigs are accomplished by reciprocal blast implemented in the DECIPHER R package using k-mer exact matching, which is identified as a syntenic hit [30]. Pairwise analysis was conducted on 64 contigs from 54 isolates containing at least one TA system; for each contig pair, the number of syntenic hits, number of syntenic blocks, and percentage coverage of contig length were recorded. Hits coverage was calculated as:
k n W k min ( l i , l j )
where Wk corresponds to width of syntenic hit k, n is the total number of syntenic hits, and li and lj correspond to the length of contigs i and j, respectively. The hits coverage was used to draw a clustered heatmap to represent the degree of overlap between the corresponding contigs.

2.9. Protein–Protein Interaction (PPI) and Functional Enrichment Analysis

A protein–protein interaction network was constructed using STRINGdb with an interaction score cutoff of 0.9 (highest confidence), and a second shell of no more than five interactors was added [31]. For each cluster of the PPI network, a GO (gene ontology) enrichment analysis was performed on the genes comprising each cluster for “Biological Process” and “Molecular Function” ontologies using the enrichGO function with default parameter, which is implemented in the R package clusterProfiler [32].

2.10. Statistical Analysis

All statistical tests were performed using R version 4.0.0 (2020-04-24) [33]. p-value < 0.05 were considered significant.

3. Results

3.1. Distribution of TA Systems across Campylobacter Isolates

The diversity of toxin-antitoxin within and between isolates allowed us to better dissect their distribution based on their structural domain and reference genome. Therefore, we performed extensive screening using TADB and PADS Arsenal databases as references and only considered aligned sequences with >60% coverage. After removing false positives due to multiple matches with the same sequence from different genome assemblies or same genomic regions in each isolate mapping to different sequences, we obtained a total of 12 toxin genes and 14 antitoxin genes (Figure S1) across 66 out of 111 isolates. From these isolates, 54 isolates possess at least one TA system, and 27 isolates possess structurally validated TA systems. According to TADB, the predicted TA systems belong to four type II TA families (Figure 1). Furthermore, the alignment quality for each toxin-antitoxin across the isolates was assessed based on a scaled score that takes into account the alignment coverage (see Methods). Most of the toxin-antitoxin genes (grouped by domain) had alignment scores > 95% (Figure S2) reflecting a high similarity between the isolate sequences and their corresponding genes from various genome assemblies. All of the toxin/antitoxin genes discussed hereafter belong to type II TA systems; however, we also obtained hits for cjrA-RNA (a type I antitoxin found in pVir plasmid identified in C. jejuni 81-176) in four isolates (Supplementary File S4) from the TADB database. Interestingly, this antitoxin was the only experimentally validated hit from TADB [34], but since we could not find its associated toxin, we decided not to perform further downstream analyses on cjrA-RNA.
Next, the distribution of toxin-antitoxin genes across isolates as a function of the isolate source and clonal complex was investigated. We wanted to check if isolates that are sampled from the same animal/clinical source or share the same clonal complex also share the same toxin-antitoxin genes. The distribution of TA genes was found to be highly sparse as a function of the isolate source (Figure 2A). On the other hand, the same distribution appears to be less sparse as a function of the clonal complex (Figure 2B). To further investigate whether the isolate source or clonal complex are dependent on the presence of different TA families, a Pearson chi-square test of independence between each clonal complex/isolate source and family pair were performed. According to both the chi-square and Fisher test, no significant association between the isolate source and family type was found, suggesting that the isolate source and TA family type are independent of each other. In contrast, there was a strong association between four CCs and three TA family types confirmed through Fisher’s exact test (Supplementary File S5), suggesting a dependent relationship between some TA family types and certain clonal complexes (Figure 2D).
Since most of the identified TA systems are based on BLAST alignment against the TADB and PADS databases, and since contigs were assembled in isolates from different veterinary/clinical sources, it is highly plausible that the matched TA systems are found in different strains. Accordingly, we wanted to connect each TA system identified in each isolate to the genome where it was identified (Figures S3 and S4). In addition, we also checked the association between TA family types and the genomic strains where they were identified (Figure 3). Interestingly, we observed that most of HicB-HicA TA systems were identified from the C. jejuni 81-176 strain, while almost all the other families were identified from C. jejuni RM1221 (Figure S3). It is also worth noting that the HicB-HicA family was the only family identified that harbors a TA-related domain. (Figure 1). On the other hand, either the toxin or antitoxin in the other families contain a TA-related domain. This pattern can be seen in Figure S4, where only the TA systems that contain both the toxin and antitoxin domain are considered and showed the absence of connections from C. jejuni RM1221.

3.2. Sequence Alignment and Clustering

Given the variety of genome assemblies from which the TA genes were identified and the heterogeneity of the assembled isolate genomes, we investigated the degree of similarity between the toxin and antitoxin genes from different genome assemblies representing the bacterial strains. Multiple sequence alignments were performed between toxin/antitoxin genes for each given domain, as well as those with no reported domain (Figure 4). Some of the domains (HicA_SF, YAFQ, and S24-LexA-like) showed >95% similarity between their constituent genes, while other domains such as HicA and HicB showed relatively low similarity. Despite the low similarity, those genes had highly conserved sites, as shown in HicB (Figure 4B), and their e-values, as reported by the CDD, with at least 10−3 and 10−8 for genes designated to harbor the HicA and HicB domains, respectively, which explains why these genes were assigned to their respective domains. Finally, we clustered all the locus tags by their sequence similarity to show the hierarchical relationship between the locus tags from different genome assemblies (Figure 4 C,D).

3.3. Synteny Analysis

In this study, we sought to characterize the TA systems in a diversified set of Campylobacter isolates with a sparse genomic organization. Hence, it is crucial to assess the homology and the degree of conservation between the isolates in terms of their constituent TA systems. Additionally, despite the fact that most TA systems are scattered in different regions across the genome, the synteny between those regions or so-called blocks between different isolates should be conserved. Therefore, we quantified the syntenic blocks between all pairs of contigs that are known to contain a TA system; we restricted the analysis only to contigs containing TA systems to increase the signal-to-noise ratio compared to a whole-genome comparison. The shared synteny between two respective contigs was quantified using the hits coverage, as described above. The contigs were then hierarchically clustered based on the coverage, and the TA systems found in each contig pair were overlaid to show which cluster of contig pairs shared a specific TA system (Figure 5). Some clusters showed high synteny and were significantly enriched in only one TA system, such as those contig pairs containing TADB|4371 (HicB-HicA) and TADB|1511 (Xre-HipA). Interestingly, there were some clusters containing contig pairs that shared two TA systems in each contig, such as (phd-doc and relBE) or (relBE and HicB-HicA), so, in those isolates, those two TA systems could coexist together.

3.4. Functional Enrichment and PPI Network

In order to comprehensively characterize the identified TA systems, the protein–protein interaction/association between the identified structural domains found in multiple toxin-antitoxin genes and their underlying functional insights were investigated. The resulting network was divided into three main clusters, including the following seed domains (i.e., the ones present in our dataset): (1) yafQ-hicA-hicB interactors, (2) lexA interactors, and (3) Fic interactors (Figure 6). Only yhfG was reported to interact with fic, and although it is an uncharacterized protein, it is believed that it might be the putative antitoxin for Fic [35]. Next, a GO enrichment analysis was performed on the remaining two clusters as a means of exploring the molecular function and biological pathways through which these genes act. Both clusters were enriched for nuclease activities, with the yafQ cluster mainly enriched for endonucleases while the LexA cluster was enriched in exonucleases and nucleotidyltransferases (Figure 7 A,E). On the other hand, they were not part of the same biological processes, as LexA was enriched in DNA repair and the SOS response in contrast to yafQ, which was enriched in mRNA metabolism and negative regulation of the gene expression.

4. Discussion

In this study, we set out to identify and characterize toxin-antitoxin modules in Campylobacter isolates sampled from different clinical/animal sources using in silico screening. We first performed a comprehensive BLAST search against both TADB and PADS arsenal databases to identify previously predicted toxin-antitoxin genes; we then performed a domain analysis to separate between TA modules where both genes contain a structurally validated TA domain and those with only one component possessing a TA-related domain (Figure S1). Most of the TA families identified were mapped to TA modules where either the toxin or antitoxin possess a TA-related domain (relBE, phd-doc, and Xre-HipA) compared to HicB-HicA as the only family where both components are structurally validated (Figure 1). All of the identified families are type II TA modules, which contain antitoxins that usually contain an N-terminal DNA-binding domain for transcriptional autoregulation and a C-terminal domain that binds directly to its toxin counterpart [37,38]. RelBE is one of the examples where the direct inactivation mediated by the antitoxin occur through the inhibition of catalysis at the toxin’s active site [39]. On the other hand, a defining feature of many toxins is their ability to bind and degrade/cleave their targets enzymatically and, thus, impede the physiological cascades, leading to bacterial growth inhibition. HicB-HicA and RelBE are both mRNA endoribonucleases that bind to and cleave cellular mRNAs, with the latter being ribosome-dependent and the former ribosome-independent [8]. HicB-HicA and RelBE modules are also functionally similar to the well-characterized MazEF TA modules, especially in E. coli (strain K-12). The toxin MAZF is a sequence-specific ribosome-independent mRNA endoribonuclease that recognizes and cleaves a ~7-nt region with a central ACA sequence at the 5′ end of the underlying mRNA, whose secondary structure conformation also influences the recognition and cleavage site [40,41]. The impact of this TA module on the physiology of E. coli cells is still debatable; however, it has been suggested that programmed cell death (PCD) is highly dependent on the presence of MazEF modules, in addition to quorum-sensing pentapeptide and the extracellular death factor (EDF) [42,43]. However, it has also been reported that MazEF-mediated PCD was not reproducible, where the E. coli strains (MC4100 and its MazEF derivative) failed to show a RelA phenotype [44] which demands further scrutiny and skepticism regarding the MazEF-mediated PCD hypothesis in future studies. It has also been suggested that both MazEF, along with DinJ-YafQ TA systems, contribute to biofilm formation [45]. Despite the strong association between HicB-HicA and RelBE functions to that of MazEF TA modules (Figure 6), we cannot conclude that this is indeed the case, especially in Campylobacter, because the stress response mediated by MazEF has been shown to be strain-dependent in E. coli [46], which makes it even more challenging to make inferences on the species level. Unfortunately, there are no/few functional studies concerned with the characterization of the stress response mediated by HicB-HicA in Campylobacter to reinforce the findings mentioned in this study.
As for the other identified TA families, like Phd-Doc (as mentioned in TADB), we observed Fic domains in the toxin component belonging to this TA module, and unlike HicB-HicA, Fic toxins were previously observed in Campylobacter fetus [47], and a Phd-Doc TA module was recently identified in C. jejuni YH002 [48]. Fic domains function through post-translational modification of their target proteins, specifically by the addition of AMP on key functional residues, thus regulating metabolic functions [49]. Accordingly, activation of the Fic toxin can slow cell growth and allow bacterial cells to enter a dormant state [47,50]. These functional cues provide insights about possible mechanisms through which Campylobacter species utilize TA modules to survive under different stresses in different environments. The last TA family identified was Xre/HipA (as mentioned in TADB), which was observed in seven Campylobacter isolates. The antitoxin of this family contained a S24-LexA-like domain, which corresponds to the well-known LexA transcriptional regulator of the SOS response to DNA damage [51] (Figure 6 and Figure 7). Although it is not considered a putative antitoxin, it was previously reported that some TA systems were regulated by the SOS response, including the yafNO TA module in E. coli [52]. This was partially explained, because this module lies downstream of the dinB gene, which encodes the translesion DNA polymerase IV that is associated with adaptive mutation; the dinB gene is also repressed by LexA and is part of the yafNO operon, which explains why this TA module is regulated by the SOS response [53]. Furthermore, there are also multiple domain architectures where S24-LexA-like domains lie downstream of the HTH (helix-turn-helix) Xre family, which constitutes a HipB antitoxin that lies in the same operon with the HipA toxin [54,55]. This might explain why TADB has attributed the S24-lexA-like domain to the Xre-HipA module. However, it is worth noting that this is just a prediction and that the role of transcriptional regulators like LexA in the regulating and determining of the functionality of TA systems is unknown [53].
Next, the overrepresentation of TA systems and families was assessed in relation to the isolate sources and MLST clonal complex (Figure 4). No significant differences between the distributions of TA families were detected among different isolate sources. However, some MLST clonal complexes overrepresented TA families, according to the Fisher exact test. MLST provided invaluable insights into the population genetics of Campylobacter species, especially in discerning bacterial relationships, albeit not as discriminatory as pulse-field gel electrophoresis (PFGE) [56]. In addition, it was proven useful to identify sequence types (ST) that fall into clusters called clonal complexes associated with specific niches, which might help in delineating transmission routes for human infection [56]. A significant association between the ST-21 complex and HicB-HicA was observed compared to other complexes; more specifically, 12 out of 13 ST-21 isolates harbor a HicB-HicA TA module. The ST-21 complex is considered one of the largest clonal complexes found in Campylobacter isolates, constituting 26% of all submitted isolates to the pubmlst database [57]. It has also been previously reported that Campylobacter isolates belonging to the ST-21 complex are resistant to ciprofloxacin [58], which would be an interesting association to investigate with respect to the presence of HicB-HicA TA modules. As per our knowledge, this is the first study reporting the association between TA modules and Campylobacter clonal complex, and further investigations are required to validate this premise.
Given the dynamic nature and high mobility of TA modules through recurring events of horizontal gene transfer, it was imperative to evaluate the shared synteny between isolates in terms of the TA modules they harbor, thus quantifying the relatedness of the isolates to each other as a function of the TA systems they share (Figure 5). Interestingly, we observed defined clusters of isolates sharing a certain TA system and even coexisting pairs of TA systems belonging to different families. This type of analysis offers a new perspective to evaluate the relationship between isolates in terms of the TA systems that reside in shared syntenic blocks.
In this study, we only scratched the surface to identify and characterize different TA modules in a fairly large set of Campylobacter isolates. However, there are still so many unanswered questions and future interrogations that need to take place. Nathan Fraikin et al. proposed some interesting challenges facing the field of TA systems that currently remain untackled. For example, how genes that are so unpredictably distributed and impermeable to fixations can be central in essential processes such as the stress response. In addition, the mechanisms by which TA systems are transferred between hosts and integrate between genomes remain unresolved [53].
More fundamentally, we still lack basic functional understandings regarding wild-type TA modules. For instance, mRNA targets of TA modules such as HicB-HicA remain largely unknown, as well as the biological processes through which these modules function. Additionally, given the heterogeneous pool of TA modules and their diversity among different bacterial species, our basic understanding about the context-specific cellular response to stress is limited. Indeed, more attention is required to properly characterize and gain mechanistic insights about how these small modules regulate cellular physiology in different environments.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4425/12/1/72/s1: Figure S1: Isolate count per toxin-antitoxin gene. Figure S2: Distribution of alignment score for each toxin-antitoxin domain. Figure S3: Circos plot showing the distribution of TA systems from different strains. Figure S4: Circos plot showing the distribution of structurally validated TA systems from different strains. Supplementary File S1: Annotation of isolates and additional metadata. Supplementary File S2: List of Campylobacter reference genomes used for genome annotation. Supplementary File S3: An annotated file for all predicted TA modules from TADB and PADS Arsenal in this study. Supplementary File S4: Predicted type I antitoxins from TADB. Supplementary File S5: Statistical association between the clonal complex and TA family type.

Author Contributions

Conceptualization, T.Z.S., and B.W.; data curation, B.W., M.A.A.-F., A.Y., and S.F.M.; formal analysis, B.W. and M.A.A.-F.; investigation, B.W., M.A.A.-F., and A.Y.; methodology, B.W. and M.A.A.; resources, M.E. and S.F.M.; visualization, B.W. and M.A.A.-F.; writing—original draft, B.W.; writing—review and editing, T.Z.S., M.E., and B.W.; project administration, T.Z.S.; supervision, T.Z.S. and M.E.; and funding acquisition, T.Z.S. and M.E. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded with ZC internal funds 004-2019 and 016-2019.

Acknowledgments

The authors are grateful to Samuel Sheppard, Ben Pascoe, and Jessica Calland at the Milner Centre of Evolution, University of Bath, UK for their excellent technical assistance in genome sequencing and assemblies.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Silva, J.; Leite, D.; Fernandes, M.; Mena, C.; Gibbs, P.A.; Teixeira, P. Campylobacter spp. as a Foodborne Pathogen: A Review. Front. Microbiol. 2011, 2. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Newell, D.G.; Mughini-Gras, L.; Kalupahana, R.S.; Wagenaar, J.A. Campylobacter epidemiology-sources and routes of transmission for human infection. Campylobacter Features Detect. Prev. Foodborne Dis. 2017, 85–110. [Google Scholar] [CrossRef]
  3. Wieczorek, K.; Wołkowicz, T.; Osek, J. Antimicrobial Resistance and Virulence-Associated Traits of Campylobacter jejuni Isolated from Poultry Food Chain and Humans with Diarrhea. Front. Microbiol. 2018, 9. [Google Scholar] [CrossRef] [PubMed]
  4. Liu, B.; Zheng, D.; Jin, Q.; Chen, L.; Yang, J. VFDB 2019: A comparative pathogenomic platform with an interactive web interface. Nucleic Acids Res. 2019, 47, D687–D692. [Google Scholar] [CrossRef] [PubMed]
  5. Inglis, R.F.; Bayramoglu, B.; Gillor, O.; Ackermann, M. The role of bacteriocins as selfish genetic elements. Biol. Lett. 2013, 9. [Google Scholar] [CrossRef] [Green Version]
  6. Forsberg, K.J.; Malik, H.S. Microbial Genomics: The Expanding Universe of Bacterial Defense Systems. Curr. Biol. 2018, 28, R361–R364. [Google Scholar] [CrossRef] [Green Version]
  7. Stern, A.; Sorek, R. The phage-host arms race: Shaping the evolution of microbes. Bioessays 2011, 33, 43–51. [Google Scholar] [CrossRef] [Green Version]
  8. Harms, A.; Brodersen, D.E.; Mitarai, N.; Gerdes, K. Toxins, Targets, and Triggers: An Overview of Toxin-Antitoxin Biology. Mol. Cell 2018, 70, 768–784. [Google Scholar] [CrossRef] [Green Version]
  9. Wang, X.; Lord, D.M.; Cheng, H.-Y.; Osbourne, D.O.; Hong, S.H.; Sanchez-Torres, V.; Quiroga, C.; Zheng, K.; Herrmann, T.; Peti, W.; et al. A new type V toxin-antitoxin system where mRNA for toxin GhoT is cleaved by antitoxin GhoS. Nat. Chem. Biol. 2012, 8, 855–861. [Google Scholar] [CrossRef] [Green Version]
  10. Aakre, C.D.; Phung, T.N.; Huang, D.; Laub, M.T. A Bacterial Toxin Inhibits DNA Replication Elongation through a Direct Interaction with the β Sliding Clamp. Mol. Cell 2013, 52, 617–628. [Google Scholar] [CrossRef] [Green Version]
  11. Gerdes, K.; Rasmussen, P.B.; Molin, S. Unique type of plasmid maintenance function: Postsegregational killing of plasmid-free cells. Proc. Natl. Acad. Sci. USA 1986, 83, 3116–3120. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Leplae, R.; Geeraerts, D.; Hallez, R.; Guglielmini, J.; Drèze, P.; Van Melderen, L. Diversity of bacterial type II toxin–antitoxin systems: A comprehensive search and functional analysis of novel families. Nucleic Acids Res. 2011, 39, 5513–5525. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Pandey, D.P.; Gerdes, K. Toxin–antitoxin loci are highly abundant in free-living but lost from host-associated prokaryotes. Nucleic Acids Res. 2005, 33, 966–976. [Google Scholar] [CrossRef] [PubMed]
  14. Masuda, H.; Inouye, M. Toxins of Prokaryotic Toxin-Antitoxin Systems with Sequence-Specific Endoribonuclease Activity. Toxins 2017, 9, 140. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Wen, J.; Fozo, E.M. sRNA Antitoxins: More than One Way to Repress a Toxin. Toxins 2014, 6, 2310–2335. [Google Scholar] [CrossRef]
  16. Van Melderen, L. Toxin–antitoxin systems: Why so many, what for? Curr. Opin. Microbiol. 2010, 13, 781–785. [Google Scholar] [CrossRef]
  17. Linton, D.; Lawson, A.J.; Owen, R.J.; Stanley, J. PCR detection, identification to species level, and fingerprinting of Campylobacter jejuni and Campylobacter coli direct from diarrheic samples. J. Clin. Microbiol. 1997, 35, 2568–2572. [Google Scholar] [CrossRef] [Green Version]
  18. Stucki, U.; Frey, J.; Nicolet, J.; Burnens, A.P. Identification of Campylobacter jejuni on the basis of a species-specific gene that encodes a membrane protein. J. Clin. Microbiol. 1995, 33, 855–859. [Google Scholar] [CrossRef] [Green Version]
  19. Gonzalez, I.; Grant, K.A.; Richardson, P.T.; Park, S.F.; Collins, M.D. Specific identification of the enteropathogens Campylobacter jejuni and Campylobacter coli by using a PCR test based on the ceuE gene encoding a putative virulence determinant. J. Clin. Microbiol. 1997, 35, 759–763. [Google Scholar] [CrossRef] [Green Version]
  20. Jolley, K.A.; Maiden, M.C. BIGSdb: Scalable analysis of bacterial genome variation at the population level. BMC Bioinform. 2010, 11, 595. [Google Scholar] [CrossRef] [Green Version]
  21. Seemann, T. Prokka: Rapid prokaryotic genome annotation. Bioinformatics 2014, 30, 2068–2069. [Google Scholar] [CrossRef] [PubMed]
  22. Xie, Y.; Wei, Y.; Shen, Y.; Li, X.; Zhou, H.; Tai, C.; Deng, Z.; Ou, H.-Y. TADB 2.0: An updated database of bacterial type II toxin–antitoxin loci. Nucleic Acids Res. 2018, 46, D749–D753. [Google Scholar] [CrossRef] [PubMed]
  23. Dehal, P.S.; Joachimiak, M.P.; Price, M.N.; Bates, J.T.; Baumohl, J.K.; Chivian, D.; Friedland, G.D.; Huang, K.H.; Keller, K.; Novichkov, P.S.; et al. MicrobesOnline: An integrated portal for comparative and functional genomics. Nucleic Acids Res. 2010, 38, D396–D400. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Zhang, Y.; Zhang, Z.; Zhang, H.; Zhao, Y.; Zhang, Z.; Xiao, J. PADS Arsenal: A database of prokaryotic defense systems related genes. Nucleic Acids Res. 2020, 48, D590–D598. [Google Scholar] [CrossRef] [PubMed]
  25. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic local alignment search tool. J. Mol. Biol. 1990, 215, 403–410. [Google Scholar] [CrossRef]
  26. Marchler-Bauer, A.; Derbyshire, M.K.; Gonzales, N.R.; Lu, S.; Chitsaz, F.; Geer, L.Y.; Geer, R.C.; He, J.; Gwadz, M.; Hurwitz, D.I.; et al. CDD: NCBI’s conserved domain database. Nucleic Acids Res. 2015, 43, D222–D226. [Google Scholar] [CrossRef] [Green Version]
  27. Quevillon, E.; Silventoinen, V.; Pillai, S.; Harte, N.; Mulder, N.; Apweiler, R.; Lopez, R. InterProScan: Protein domains identifier. Nucleic Acids Res. 2005, 33, W116–W120. [Google Scholar] [CrossRef] [Green Version]
  28. Bodenhofer, U.; Bonatesta, E.; Horejš-Kainrath, C.; Hochreiter, S. msa: An R package for multiple sequence alignment. Bioinformatics 2015, 31, 3997–3999. [Google Scholar] [CrossRef]
  29. Charif, D.; Lobry, J.R. SeqinR 1.0-2: A Contributed Package to the R Project for Statistical Computing Devoted to Biological Sequences Retrieval and Analysis. In Structural Approaches to Sequence Evolution: Molecules, Networks, Populations; Biological and Medical Physics, Biomedical Engineering; Bastolla, U., Porto, M., Roman, H.E., Vendruscolo, M., Eds.; Springer: Berlin/Heidelberg, Germany, 2007; pp. 207–232. ISBN 978-3-540-35306-5. [Google Scholar]
  30. Wright, E.S. Using DECIPHER v2.0 to Analyze Big Biological Sequence Data in R. R J. 2016, 8, 352–359. [Google Scholar] [CrossRef] [Green Version]
  31. Szklarczyk, D.; Gable, A.L.; Lyon, D.; Junge, A.; Wyder, S.; Huerta-Cepas, J.; Simonovic, M.; Doncheva, N.T.; Morris, J.H.; Bork, P.; et al. STRING v11: Protein–protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019, 47, D607–D613. [Google Scholar] [CrossRef] [Green Version]
  32. Yu, G.; Wang, L.-G.; Han, Y.; He, Q.-Y. clusterProfiler: An R Package for Comparing Biological Themes Among Gene Clusters. OMICS 2012, 16, 284–287. [Google Scholar] [CrossRef] [PubMed]
  33. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2020. [Google Scholar]
  34. Shen, Z.; Patil, R.D.; Sahin, O.; Wu, Z.; Pu, X.-Y.; Dai, L.; Plummer, P.J.; Yaeger, M.J.; Zhang, Q. Identification and functional analysis of two toxin-antitoxin systems in Campylobacter jejuni. Mol. Microbiol. 2016, 101, 909–923. [Google Scholar] [CrossRef] [PubMed]
  35. Rajagopala, S.V.; Sikorski, P.; Kumar, A.; Mosca, R.; Vlasblom, J.; Arnold, R.; Franca-Koh, J.; Pakala, S.B.; Phanse, S.; Ceol, A.; et al. The binary protein-protein interaction landscape of Escherichia coli. Nat. Biotechnol. 2014, 32, 285–290. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Tandon, H.; Sharma, A.; Wadhwa, S.; Varadarajan, R.; Singh, R.; Srinivasan, D.N.; Sandhya, S. Bioinformatic and mutational studies of related toxin–antitoxin pairs in M. tuberculosis predict and identify key functional residues. J. Biol. Chem. 2019. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Goeders, N.; Van Melderen, L. Toxin-Antitoxin Systems as Multilevel Interaction Systems. Toxins 2014, 6, 304–324. [Google Scholar] [CrossRef] [Green Version]
  38. Chan, W.T.; Espinosa, M.; Yeo, C.C. Keeping the Wolves at Bay: Antitoxins of Prokaryotic Type II Toxin-Antitoxin Systems. Front. Mol. Biosci. 2016, 3, 9. [Google Scholar] [CrossRef] [Green Version]
  39. Bøggild, A.; Sofos, N.; Andersen, K.R.; Feddersen, A.; Easter, A.D.; Passmore, L.A.; Brodersen, D.E. The Crystal Structure of the Intact E. coli RelBE Toxin-Antitoxin Complex Provides the Structural Basis for Conditional Cooperativity. Structure 2012, 20, 1641–1648. [Google Scholar] [CrossRef] [Green Version]
  40. Zhang, Y.; Zhang, J.; Hara, H.; Kato, I.; Inouye, M. Insights into the mRNA cleavage mechanism by MazF, an mRNA interferase. J. Biol. Chem. 2005, 280, 3143–3150. [Google Scholar] [CrossRef] [Green Version]
  41. Culviner, P.H.; Laub, M.T. Global Analysis of the E. coli Toxin MazF Reveals Widespread Cleavage of mRNA and the Inhibition of rRNA Maturation and Ribosome Biogenesis. Mol. Cell 2018, 70, 868–880.e10. [Google Scholar] [CrossRef] [Green Version]
  42. Kolodkin-Gal, I.; Hazan, R.; Gaathon, A.; Carmeli, S.; Engelberg-Kulka, H. A linear pentapeptide is a quorum-sensing factor required for mazEF-mediated cell death in Escherichia coli. Science 2007, 318, 652–655. [Google Scholar] [CrossRef] [Green Version]
  43. Kolodkin-Gal, I.; Engelberg-Kulka, H. The extracellular death factor: Physiological and genetic factors influencing its production and response in Escherichia coli. J. Bacteriol. 2008, 190, 3169–3175. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Ramisetty, B.C.M.; Raj, S.; Ghosh, D. Escherichia coli MazEF toxin-antitoxin system does not mediate programmed cell death. J. Basic Microbiol. 2016, 56, 1398–1402. [Google Scholar] [CrossRef] [PubMed]
  45. Kolodkin-Gal, I.; Verdiger, R.; Shlosberg-Fedida, A.; Engelberg-Kulka, H. A differential effect of E. coli toxin-antitoxin systems on cell death in liquid media and biofilm formation. PLoS ONE 2009, 4, e6785. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Hazan, R.; Engelberg-Kulka, H. Escherichia coli mazEF-mediated cell death as a defense mechanism that inhibits the spread of phage P1. Mol. Genet. Genom. 2004, 272, 227–234. [Google Scholar] [CrossRef]
  47. Sprenger, H.; Kienesberger, S.; Pertschy, B.; Pöltl, L.; Konrad, B.; Bhutada, P.; Vorkapic, D.; Atzmüller, D.; Feist, F.; Högenauer, C.; et al. Fic Proteins of Campylobacter fetus subsp. venerealis Form a Network of Functional Toxin–Antitoxin Systems. Front. Microbiol. 2017, 8. [Google Scholar] [CrossRef] [Green Version]
  48. Ghatak, S.; He, Y.; Reed, S.; Irwin, P. Comparative Genomic Analysis of a Multidrug-Resistant Campylobacter jejuni Strain YH002 Isolated from Retail Beef Liver. Foodborne Pathog. Dis. 2020. [Google Scholar] [CrossRef]
  49. Garcia-Pino, A.; Zenkin, N.; Loris, R. The many faces of Fic: Structural and functional aspects of Fic enzymes. Trends Biochem. Sci. 2014, 39, 121–129. [Google Scholar] [CrossRef]
  50. Castro-Roa, D.; Garcia-Pino, A.; De Gieter, S.; van Nuland, N.A.J.; Loris, R.; Zenkin, N. The Fic protein Doc uses an inverted substrate to phosphorylate and inactivate EF-Tu. Nat. Chem. Biol. 2013, 9, 811–817. [Google Scholar] [CrossRef] [Green Version]
  51. Maslowska, K.H.; Makiela-Dzbenska, K.; Fijalkowska, I.J. The SOS system: A complex and tightly regulated response to DNA damage. Environ. Mol. Mutagenes. 2019, 60, 368–384. [Google Scholar] [CrossRef] [Green Version]
  52. Singletary, L.A.; Gibson, J.L.; Tanner, E.J.; McKenzie, G.J.; Lee, P.L.; Gonzalez, C.; Rosenberg, S.M. An SOS-Regulated Type 2 Toxin-Antitoxin System. J. Bacteriol. 2009, 191, 7456–7465. [Google Scholar] [CrossRef] [Green Version]
  53. Fraikin, N.; Goormaghtigh, F.; Melderen, L.V. Type II Toxin-Antitoxin Systems: Evolution and Revolutions. J. Bacteriol. 2020, 202. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Geer, L.Y.; Domrachev, M.; Lipman, D.J.; Bryant, S.H. CDART: Protein homology by domain architecture. Genome Res. 2002, 12, 1619–1623. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Hansen, S.; Vulić, M.; Min, J.; Yen, T.-J.; Schumacher, M.A.; Brennan, R.G.; Lewis, K. Regulation of the Escherichia coli HipBA Toxin-Antitoxin System by Proteolysis. PLoS ONE 2012, 7, e39185. [Google Scholar] [CrossRef]
  56. McTavish, S.M.; Pope, C.E.; Nicol, C.; Campbell, D.; French, N.; Carter, P.E. Multilocus sequence typing of Campylobacter jejuni, and the correlation between clonal complex and pulsed-field gel electrophoresis macrorestriction profile. FEMS Microbiol. Lett. 2009, 298, 149–156. [Google Scholar] [CrossRef]
  57. Best, E.L.; Fox, A.J.; Frost, J.A.; Bolton, F.J. Identification of Campylobacter jejuni Multilocus Sequence Type ST-21 Clonal Complex by Single-Nucleotide Polymorphism Analysis. J. Clin. Microbiol. 2004, 42, 2836–2839. [Google Scholar] [CrossRef] [Green Version]
  58. Cody, A.J.; McCarthy, N.M.; Wimalarathna, H.L.; Colles, F.M.; Clark, L.; Bowler, I.C.J.W.; Maiden, M.C.J.; Dingle, K.E. A Longitudinal 6-Year Study of the Molecular Epidemiology of Clinical Campylobacter Isolates in Oxfordshire, United Kingdom. J. Clin. Microbiol. 2012, 50, 3193–3201. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Isolate count per toxin-antitoxin (TA) family type. Bar lengths correspond to the number of isolates containing TA systems belonging to a given TA family. Blue bars represent TA systems where both toxin and antitoxin genes harbor a TA-related domain, while red bars represent TA systems where either the toxin or antitoxin gene harbor a TA-related domain.
Figure 1. Isolate count per toxin-antitoxin (TA) family type. Bar lengths correspond to the number of isolates containing TA systems belonging to a given TA family. Blue bars represent TA systems where both toxin and antitoxin genes harbor a TA-related domain, while red bars represent TA systems where either the toxin or antitoxin gene harbor a TA-related domain.
Genes 12 00072 g001
Figure 2. Heatmap showing the distribution of TA genes and families across isolates, isolate sources, and clonal complexes. In (A,B), each cell represents the occurrence of a TA gene (represented by the locus tag shown on the x-axis) for each isolate (shown on the y-axis). For the identification of a hit, a cut-off threshold of alignment length > 60% was applied, as proposed in [36]. The color conventions show the scaled score representing the alignment quality (see Methods). TA genes are grouped by their protein domain (shown as a colored track on top). In (C,D), each cell is colored by a normalized score, where the number of isolates in a given isolate source/clonal complex and a given family type is divided by the total number of isolates in a given isolate source/clonal complex, regardless of family type. In addition, cells are labeled by their standardized residuals as calculated from Pearson’s chi-square test. Residuals higher or lower than 2 or −2, respectively, show a statistically significant association (p < 0.05). (A) Isolates are grouped based on the isolate source (represented by a colored track on the left). (B) Isolates are grouped based on their clonal complex (represented by a colored track on the left). (C) Hierarchically clustered heat map between the TA family type and clonal complex. (D) Hierarchically clustered heat map between the TA family type and isolate source.
Figure 2. Heatmap showing the distribution of TA genes and families across isolates, isolate sources, and clonal complexes. In (A,B), each cell represents the occurrence of a TA gene (represented by the locus tag shown on the x-axis) for each isolate (shown on the y-axis). For the identification of a hit, a cut-off threshold of alignment length > 60% was applied, as proposed in [36]. The color conventions show the scaled score representing the alignment quality (see Methods). TA genes are grouped by their protein domain (shown as a colored track on top). In (C,D), each cell is colored by a normalized score, where the number of isolates in a given isolate source/clonal complex and a given family type is divided by the total number of isolates in a given isolate source/clonal complex, regardless of family type. In addition, cells are labeled by their standardized residuals as calculated from Pearson’s chi-square test. Residuals higher or lower than 2 or −2, respectively, show a statistically significant association (p < 0.05). (A) Isolates are grouped based on the isolate source (represented by a colored track on the left). (B) Isolates are grouped based on their clonal complex (represented by a colored track on the left). (C) Hierarchically clustered heat map between the TA family type and clonal complex. (D) Hierarchically clustered heat map between the TA family type and isolate source.
Genes 12 00072 g002
Figure 3. Circos plot showing the distribution of type II TA families from different strains. The lower nodes correspond to different genomic strains from which arcs are drawn towards the corresponding TA families in the upper nodes. The thickness of each arc corresponds to the number of isolates containing TA systems that belong to a given TA family and are matched from a given strain. C.j, Campylobacter jejuni and C.c, Campylobacter coli.
Figure 3. Circos plot showing the distribution of type II TA families from different strains. The lower nodes correspond to different genomic strains from which arcs are drawn towards the corresponding TA families in the upper nodes. The thickness of each arc corresponds to the number of isolates containing TA systems that belong to a given TA family and are matched from a given strain. C.j, Campylobacter jejuni and C.c, Campylobacter coli.
Genes 12 00072 g003
Figure 4. Multiple sequence alignment of locus tags for each domain. (A) Sequence alignment of toxin genes for each identified domain. (B) Sequence alignment of antitoxin genes for each identified domain. (C) Dendogram representing sequence clustering of all toxin genes, locus tags highlighted in red contain toxin domains. (D) Dendogram representing sequence clustering of all antitoxin genes, locus tags highlighted in red contain antitoxin domains. Residues highlighted in blue are ≥50% conserved, and those highlighted in pink show similar residues.
Figure 4. Multiple sequence alignment of locus tags for each domain. (A) Sequence alignment of toxin genes for each identified domain. (B) Sequence alignment of antitoxin genes for each identified domain. (C) Dendogram representing sequence clustering of all toxin genes, locus tags highlighted in red contain toxin domains. (D) Dendogram representing sequence clustering of all antitoxin genes, locus tags highlighted in red contain antitoxin domains. Residues highlighted in blue are ≥50% conserved, and those highlighted in pink show similar residues.
Genes 12 00072 g004
Figure 5. Clustered heatmap mapping the coverage of syntentic hits between each contig. Each cell is labeled by an index indicating the TA systems shared by each contig pair (table above). The color conventions show the coverage proportion between each contig pair; 1 means 100% coverage, and 0 indicates no synteny between corresponding contigs. The columns and row labels represent the isolate id, along with the contig number.
Figure 5. Clustered heatmap mapping the coverage of syntentic hits between each contig. Each cell is labeled by an index indicating the TA systems shared by each contig pair (table above). The color conventions show the coverage proportion between each contig pair; 1 means 100% coverage, and 0 indicates no synteny between corresponding contigs. The columns and row labels represent the isolate id, along with the contig number.
Genes 12 00072 g005aGenes 12 00072 g005b
Figure 6. Protein–protein interaction (PPI) network. Each node represents a given protein, and the edges represent the degree of association between the given nodes. The network is divided into 3 independent clusters; the top-left cluster shows the PPIs between the interactors with yafQ-hicA-hicB as the seed/query nodes. The bottom-right cluster shows the PPIs between the interactors with lexA as the seed/query node. The top-right cluster shows a single interaction involving fic as the seed node.
Figure 6. Protein–protein interaction (PPI) network. Each node represents a given protein, and the edges represent the degree of association between the given nodes. The network is divided into 3 independent clusters; the top-left cluster shows the PPIs between the interactors with yafQ-hicA-hicB as the seed/query nodes. The bottom-right cluster shows the PPIs between the interactors with lexA as the seed/query node. The top-right cluster shows a single interaction involving fic as the seed node.
Genes 12 00072 g006
Figure 7. Functional enrichment analysis. (AD) Show the enrichment analysis of the cluster involving lexA as the seed node, and (EH) show the enrichment results of the cluster involving yafQ-hicA-hicB as seed nodes (Figure 6). In panels (A,B,E,F), upset plots are drawn to visualize the intersection between the gene ontology (GO) terms in a given ontology (horizontal labels). Every possible intersection is represented by the bottom plot, and the frequency (i.e., number of genes) is shown on the top vertical bar plot. In panels (C,D,G,H), circular gene-concept networks are drawn to visualize the interaction between each gene and its associated GO term. The edges are colored according to the GO term, and the node size reflects the number of genes associated with each category. Panels (A,C,E,G) correspond to enrichment results for GO “molecular function” ontology, while panels (B,D,F,H) correspond to enrichment results for GO “biological process” ontology.
Figure 7. Functional enrichment analysis. (AD) Show the enrichment analysis of the cluster involving lexA as the seed node, and (EH) show the enrichment results of the cluster involving yafQ-hicA-hicB as seed nodes (Figure 6). In panels (A,B,E,F), upset plots are drawn to visualize the intersection between the gene ontology (GO) terms in a given ontology (horizontal labels). Every possible intersection is represented by the bottom plot, and the frequency (i.e., number of genes) is shown on the top vertical bar plot. In panels (C,D,G,H), circular gene-concept networks are drawn to visualize the interaction between each gene and its associated GO term. The edges are colored according to the GO term, and the node size reflects the number of genes associated with each category. Panels (A,C,E,G) correspond to enrichment results for GO “molecular function” ontology, while panels (B,D,F,H) correspond to enrichment results for GO “biological process” ontology.
Genes 12 00072 g007
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wadie, B.; Abdel-Fattah, M.A.; Yousef, A.; Mouftah, S.F.; Elhadidy, M.; Salem, T.Z. In Silico Characterization of Toxin-Antitoxin Systems in Campylobacter Isolates Recovered from Food Sources and Sporadic Human Illness. Genes 2021, 12, 72. https://doi.org/10.3390/genes12010072

AMA Style

Wadie B, Abdel-Fattah MA, Yousef A, Mouftah SF, Elhadidy M, Salem TZ. In Silico Characterization of Toxin-Antitoxin Systems in Campylobacter Isolates Recovered from Food Sources and Sporadic Human Illness. Genes. 2021; 12(1):72. https://doi.org/10.3390/genes12010072

Chicago/Turabian Style

Wadie, Bishoy, Mohamed A. Abdel-Fattah, Alshymaa Yousef, Shaimaa F. Mouftah, Mohamed Elhadidy, and Tamer Z. Salem. 2021. "In Silico Characterization of Toxin-Antitoxin Systems in Campylobacter Isolates Recovered from Food Sources and Sporadic Human Illness" Genes 12, no. 1: 72. https://doi.org/10.3390/genes12010072

APA Style

Wadie, B., Abdel-Fattah, M. A., Yousef, A., Mouftah, S. F., Elhadidy, M., & Salem, T. Z. (2021). In Silico Characterization of Toxin-Antitoxin Systems in Campylobacter Isolates Recovered from Food Sources and Sporadic Human Illness. Genes, 12(1), 72. https://doi.org/10.3390/genes12010072

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop