Paeniclostridium sordellii and Clostridioides difficile encode similar and clinically relevant tetracycline resistance loci in diverse genomic locations
BMC Microbiologyvolume 19, Article number: 53 (2019)
With the current rise of antibiotic resistance in bacteria, it is important to monitor the efficacy of antimicrobials in clinical use. Paeniclostridium sordellii (previously Clostridium sordellii) is a bacterial pathogen that causes human uterine infection after spontaneous or medically induced abortion, for which mortality rates approach 100%. Prophylactic antibiotics have been recommended for individuals undergoing medically-induced abortion, one of which is doxycycline, a member of the tetracycline antibiotic family. However, tetracycline resistance had not been well characterized in P. sordellii. This study therefore aimed to determine the levels of tetracycline resistance in P. sordellii isolates, and to identify associated loci and their genomic locations.
Using a MIC assay, five of 24 P. sordellii isolates were found to be resistant to tetracycline, minocycline, and importantly, doxycycline. Analysis of genome sequence data from 46 isolates found that phenotypically resistant isolates encoded a variant of the Clostridium perfringens tetracycline resistance determinant Tet P. Bioinformatic analysis and comparison of the regions surrounding these determinants found variation in the genomic location of Tet P among P. sordellii isolates. The core genome comparison of the 46 isolates revealed genetic diversity and the absence of dominant genetic types among the isolates. There was no strong association between geographic location of isolation, animal host or Tet P carriage with isolate genetic type. Furthermore, the analysis of the Tet P genotype revealed that Tet P is encoded chromosomally, or on one of two, novel, small plasmids, all consistent with multiple acquisition and recombination events. BLAST analysis of Clostridioides difficile draft genome sequences also identified a Tet P locus, the genomic location of which demonstrated an evolutionary relationship with the P. sordellii locus.
The Tet P determinant is found in variable genomic locations within diverse human and animal isolates of P. sordellii and C. difficile, which suggests that it can undergo horizontal transfer, and may disseminate tetracycline resistance between clostridial species. Doxycycline is a suggested prophylactic treatment for P. sordellii infections, however, a small sub-set of the isolates tested are resistant to this antibiotic. Doxycycline may therefore not be an appropriate prophylactic treatment for P. sordellii infections.
With the increase in antibiotic resistance among bacterial populations there is a growing need to develop new antimicrobial drugs, and to investigate the efficacy of those already in use. A family of antimicrobials that have had broad-spectrum applications are the tetracyclines . Tetracyclines are a family of compounds that bind to the 30S ribosomal subunit of the bacterial ribosome and block protein synthesis through inhibition of tRNA binding . However, due to the widespread use of tetracyclines, many resistance determinants are now found in diverse bacterial species . These tetracycline resistance determinants are often associated with mobile genetic elements and act via two main mechanisms [2, 3]. The first mechanism utilizes efflux proteins of the major facilitator superfamily that actively pump tetracyclines out of the bacterial cell . The second mechanism involves ribosomal protection proteins (RPPs) , which bind to the bacterial ribosome, leading to a conformational change and subsequent dissociation of the tetracyclines . While only one type of determinant is generally required to provide tetracycline resistance, phenotypic differences can be observed depending on the type of tetracycline administered. For example, RPPs can also protect against minocycline whereas an efflux protein may not .
Despite the widespread emergence of tetracycline resistance, this antibiotic family is still being used for the treatment of specific conditions such as infections with Brucella or Chlamydia species [7, 8]. They are also recommended for the prophylaxis of infection against organisms known to cause intrauterine infections, including those caused by Paeniclostridium sordellii . P. sordellii, previously known as Clostridium sordellii, is a Gram-positive pathogen that causes severe oedemic, myonecrotic or enterotoxic infections in humans, and farm animals including foals, sheep and cattle . Of particular concern over the last 20 years is an association of P. sordellii-mediated uterine infection post childbirth, spontaneous abortion and medically induced abortion, with the drug misoprostol [10, 11]. While rare, in each of these cases patients have displayed a sudden infection onset, leading to mortality rates approaching 100% [9, 10]. These cases, along with other causes of post-abortive infection, have led to expert bodies recommending antibiotic prophylaxis prior to surgical and medical abortion . Some of these recommendations include the use of the tetracycline family antibiotic doxycycline, despite levels of tetracycline resistance not being well characterized in P. sordellii .
In general, P. sordellii strains are assumed to be tetracycline sensitive [10, 12, 13], despite previous studies showing that tetracycline resistance occurs in this organism [14, 15]. One study in particular identified three P. sordellii isolates from cattle with malignant oedema that were resistant to oxytetracycline, with subsequent PCR analysis indicating that all three isolates carried the tetA(P) and tetB(P) tetracycline resistance genes . These genes are found in Clostridium perfringens within the Tet P determinant (Additional file 1: Figure S1) , although their presence in isolates of Clostridium septicum and Clostridioidies difficile has also been reported [14, 16]. The Tet P determinant is comprised of two genes: tetA(P), which encodes a tetracycline efflux protein, and tetB(P), which encodes a ribosomal protection protein  (Additional file 1: Figure S1). The genes are transcriptionally linked and have an overlap of 17 bp in their protein coding regions. In C. perfringens, the Tet P determinant has only been characterized on large conjugative plasmids related to the plasmid pCW3 [17, 18]. Two genes are found downstream of Tet P on C. perfringens pCW3- like plasmids, specifically, pcw303 that encodes a product of unknown function, and regA, which encodes a putative AraC-like transcriptional regulator [19, 20]. While it appears that Tet P is present in other clostridial members, the genomic context of this locus and any association with mobile genetic elements is unknown.
Here we determined the tetracycline resistance profiles of a diverse group of P. sordellii isolates and identified the genes and genomic location of their tetracycline resistance loci. Using MIC assays, we identified five of 24 isolates that were resistant to three common members of the tetracyclines. BLAST analysis of the genome sequences of these five tetracycline resistant (TetR) isolates identified the Tet P determinant in all of them. Bioinformatic analysis of these genome sequences, along with all publically available P. sordellii genomes, showed that Tet P is not carried on pCW3-like plasmids as seen in C. perfringens but instead is found in highly variable locations on either the chromosome or on small plasmids. Analysis of C. difficile genome sequences also identified Tet P loci and has revealed a link with the Tet P loci observed in P. sordellii.
Diverse P. sordellii strains are resistant to tetracyclines
To test if tetracycline resistance is present in P. sordellii, a resistance profile was determined for 24 geographically diverse human and animal isolates (Additional file 2: Table S2). Resistance to tetracycline, minocycline and doxycycline was examined in these isolates in comparison to positive and negative C. perfringens controls (Table 1). Of the 24 P. sordellii isolates, five displayed much higher minimum inhibitory concentration (MIC) values for all three antibiotics when compared to the other P. sordellii isolates and the C. perfringens negative control (Table 1). The level of resistance to the tetracyclines appeared similar for most of the P. sordellii isolates, however, isolate R28058 displayed MIC values two to four times greater than those of the other four resistant isolates (Table 1). R28058 is the only phenotypically TetR isolate that belongs to the most distant clade of P. sordellii (Fig. 1) , thus the increased resistance to the tetracyclines in this isolate may reflect intrinsic isolate differences that impact on the resistance phenotype.
While obtained from a small sample size, the finding that ~ 20% of the P. sordellii isolates tested are resistant to tetracyclines is of concern because doxycycline has been recommended for prophylaxis against post-abortive uterine infections, including those caused by P. sordellii . A recent screening of women for vaginal and rectal carriage of P. sordellii and C. perfringens concluded that both organisms were transient and occurred at a low rate and therefore screening or prophylactic approaches were not advisable . However, due to the severe nature of infection with these organisms, and the fact that prophylaxis does not solely target the clostridia , it may be that alternatives to doxycycline should be used for prophylaxis prior to abortion, such as metronidazole or metronidazole combined with doxycycline, as has been recommended by other expert agencies . The collection of a larger number of isolates followed by antibiotic susceptibility testing is needed to fully characterise the burden of tetracycline resistance among the P. sordellii population.
Tetracycline resistance in P. sordellii is encoded by a variant of the C. perfringens Tet P determinant
The genome sequencing data for 46 P. sordellii isolates, which included the 24 isolates used in our antibiotic resistance screening, were analysed for the presence of tetracycline resistance genes. Preliminary evaluation of the quality of the publically available data sets revealed two isolates, W2946 and W2947, were not mono-isolates and were removed from the analysis. All five TetR isolates as determined by MIC encoded a variant of the C. perfringens tetracycline resistance determinant, Tet P (Additional file 1: Figure S1). Analysis of the P. sordellii genomic data detected a further eight isolates that carried a Tet P locus, which is therefore present in ~ 30% of isolates tested. This determinant was absent in the genome sequences of all 19 phenotypically tetracycline sensitive (TetS) P. sordellii isolates, with the exception of R15892. PCR analysis of two independent stocks of R15892 could not detect Tet P. We therefore believe that although the published genome sequence of R15892 contains Tet P, the isolate used in the MIC assay has lost this element. This is further validated by our analysis of the genome sequence of this isolate that indicates Tet P is encoded on a plasmid; therefore if the isolate lost the plasmid it may lose its resistance to tetracycline.
In the 13 Tet P-containing P. sordellii isolates, tetA(P) and tetB(P) were the same size and arranged as seen in C. perfringens (Additional file 1: Figure S1). The coding sequences of the genes had ~ 88.3–89.5% and ~ 96.5–97.9% identity to TetA(P) and TetB(P) from C. perfringens, respectively (Additional file 3: Table S1). Sequence comparison of the 13 Tet P loci found in the P. sordellii isolates identifies two distinct locus types (Fig. 1, Additional file 3: Table S1), with five isolates carrying subtype 1 (Fig. 1, red) and the remaining eight isolates carrying subtype 2 (Fig. 1, blue). The upstream regulatory sequences present for Tet P in C. perfringens are observed in the majority of P. sordellii isolates (Additional file 3: Table S1), including the P3 promoter and T1 transcriptional repressor sequence (Additional file 4: Figure S2) . While the sequence of the upstream regions differs between the two P. sordellii Tet P locus subtypes, sequence identity is observed within the subtypes (Additional file 3: Table S1), although isolate W2945 displays a deletion which has resulted in the absence of T1 (Additional file 4: Figure S2). The lack of T1 may cause an upregulation of Tet P transcription , however, further studies are required to examine this hypothesis.
The genomic location of Tet P in P. sordellii isolates is highly variable
In eight of the 13 Tet P+ P. sordellii isolates, the determinant was found on large contigs predicted to be part of the chromosome (via a comparison with the partially closed genome sequence of the P. sordellii type strain ATCC9714 (Fig. 2a) and the presence of chromosomal genes in flanking regions). These isolates carry Tet P in place of two putative genes encoding proteins of unknown function when compared to the genome sequence of ATCC9714 (Fig. 2a), however, the regions directly flanking Tet P in these eight isolates showed significant variation (Fig. 2a). Three isolates (SSCC32135, JGS6961 and R32462) only contain Tet P flanked by the pCW3 homologs regA (purple) and pcw303 (pink) (Additional file 1: Figure S1, Fig. 2a), however, the presence of regA upstream of Tet P in P. sordellii differs to that of C. perfringens. The product of regA is a predicted AraC-like inducible regulator, hypothesized in C. perfringens to play a role in the inducible tetracycline phenotype observed in some isolates , while no function can be predicted for the product of pcw303. While the conservation of regA between Tet P carrying isolates of P. sordellii and C. perfringens suggests that it may play some role in tetracycline resistance, nine of the 13 Tet P+ P. sordellii isolates do not encode RegA anywhere on the genome, suggesting that RegA is not required for tetracycline resistance. The absence of regA in four of the isolates that chromosomally encode Tet P (Fig. 2a) may be explained by the presence of a 19 bp direct repeat that also flanks regA in all remaining isolates and which could serve as a homologous recombination site that leads to the loss of this gene (Additional file 5: Figure S3). Apart from containing pCW3-related genes within it, regions flanking this variable locus among the P. sordellii isolates show no relationship to plasmid sequences. These flanking regions are also highly conserved among P. sordellii isolates and encode tRNAs, chromosomal DNA partitioning proteins and ribosomal components. We therefore believe that this region of the genome is chromosomal and therefore suggest that horizontally acquired DNA has integrated into the same chromosomal hotspot.
Three isolates that encode Tet P chromosomally (W2945, SSCC32135, and W3026) have numerous ORFs positioned either upstream or downstream of Tet P (Fig. 2a), but these flanking regions show little homology between isolates. There are exceptions however, with a region encoding eight ORFs downstream of Tet P in SSCC18392 that has high identity (~ 96%) to a region upstream of Tet P in W3026 (Fig. 2a). The Tet P associated locus of W3026 is also almost entirely conserved with isolate W2922, however the Tet P determinant is absent, instead being located on a separate contig which is flanked by regions of homology to the site of deletion (Fig. 2a). Analysis of the isolate’s pan-genome did not indicate a mixed isolate. We suspect that the two contigs represent a sub-population of the single sequenced isolate in which the Tet P region has potentially looped out or been deleted and subsequently lost from the cell.
The predicted functions of the products encoded in these variable regions are associated with mobile genetic elements or phage, particularly in isolate SSCC18392 (Fig. 2a). Such putative products include a phage DNA methyl transferase, phage primase, a helicase and a type III restriction modification system protein . Despite the differences between the Tet P flanking regions in each isolate, a common feature is that they each encode a predicted large serine-recombinase (Fig. 2a) . No conserved direct or inverted repeat could be identified within the region that could explain the high variation and apparent genetic rearrangement seen among these isolates, however, the conservation of a large serine recombinase gene among these isolates indicates that its product may play a role in this variation. The potential mobility of this chromosomal Tet P locus is further evidenced by their distribution across isolates from distinct clades of P. sordellii, and suggests multiple acquisition and rearrangement events among P. sordellii isolates (Fig. 1).
The remaining five Tet P+ isolates of P. sordellii encode Tet P on one of two small plasmids (Fig. 2b, c). The presence of these plasmids and the chromosomal location of TetP in other isolates of P. sordellii was confirmed using Southern hybridization on five Tet P+ isolates (Additional file 6: Figure S4). The sequence identity between subtypes of Tet P among the population also correlates with genomic location, with sub-type 1 all being plasmid associated and sub-type 2 appearing to be carried on the chromosome (Figs. 1 and 2a, c). Isolates R15892 and R28058 carry Tet P on a novel ~ 9.4 kb plasmid (Fig. 2c), designated pCSTC1 (Fig. 2b). Putative functions could not be assigned to the other four small ORFs on pCSTC1, however pCSTC1_002 contains a predicted C-terminal DNA/RNA-binding domain  that may encode a replication initiation (Rep) protein. The remaining three isolates (SSCC37615, SSCC18838 and R32977) encode Tet P on a novel ~ 9.8 kb plasmid (Fig. 2c), designated pCSTC2 (Fig. 2b). Apart from the high region of identity (~ 99%) of Tet P, the two pCSTC plasmids have no significant homology to one another (Fig. 2b). A putative replication gene, rep, was identified on pCSTC2 , the product of which has 62% identity to the Rep protein from the C. perfringens plasmid pIP404 . pIP404 is used as the base replicon for many C. perfringens shuttle vectors . An origin of replication with significant identity to that of pIP404 is also present on pCSTC2 . Transformation of C. perfringens JIR325 (a TetS isolate), selecting for TetR, using genomic DNA from SSCC37615, SSCC18838 and R28058 resulted in transformation frequencies of 8.9 × 104, 7.2 × 104 and < 10 cfu/ml (below the detection limit), respectively. This finding indicates that pCSTC2, but not pCSTC1, can replicate in C. perfringens, and this was further confirmed using PCR (Additional file 7: Figure S5). These results suggest that an evolutionary relationship exists between pIP404-like plasmids and pCSTC2.
A predicted MobA/MobL family domain (pfam03389) protein is encoded on pCSTC2 (Fig. 2b) . This family of proteins includes MobA from the E. coli plasmid RSF1010, a relaxase that initiates transfer of this plasmid to a recipient cell . The presence of a Mob protein indicates that pCSTC2 may undergo conjugative transfer, thereby transferring tetracycline resistance to sensitive P. sordellii isolates. Bacterial matings were performed using SSCC37651 (pCSTC2), SSCC18838 (pCSTC2) and R28058 (pCSTC1) as donors with TetS P. sordellii and C. perfringens isolates. No transconjugants could be detected under the conditions tested, and therefore the mobility of either pCSTC plasmid could not be confirmed during this study. However, by examining the relatedness of strains that carry these plasmids there is evidence of their mobility. Isolates that carry pCSTC1 or pCSTC2 are not all clonal, but are instead found in distantly related isolates (Fig. 1). The fact that isolates from distinct clades encode identical plasmids suggests a recent acquisition event, thereby indicating that both Tet P-encoding plasmids may be able to undergo mobilisation between P. sordellii strains.
Tet P is found in human clinical isolates of C. difficile
Recent work has detected Tet P in three ribotype 0014 porcine C. difficile isolates . The genomic location of Tet P in all isolates was found to be identical and not associated with plasmid or mobile genetic element-associated genes . We looked for Tet P among draft genome sequences of C. difficile isolates obtained from an Australian hospital. Two isolates, MCD43 and MCD46, contained the Tet P determinant within contigs that were ~ 49.9 kb and ~ 123.5 kb in length, respectively. The nucleotide sequence identity between these contigs was 100%. The predicted coding sequences of TetA(P) and TetB(P) from MCD43/46 were identical to those of the porcine C. difficile isolates, as well as being identical to P. sordellii Tet P subtype 2 (Additional file 3: Table S1). The upstream Tet P regulatory region of MCD43/46 was identical to that of P. sordellii W2945, including the deletion of the T1 transcriptional terminator (Additional file 4: Figure S2). This region in C. difficile porcine isolates, despite the presence of a few SNPs, was highly related to that of P. sordellii Tet P subtype 2-encoding strains and unlike MCD43/46 contained the T1 transcriptional terminator (Additional file 3: Table S1, Additional file 4: Figure S2).
The MCD46 contig was compared to Tet P C. difficile genomes within GenBank, and a subsequent BLAST alignment was performed between the Tet P locus of this isolate, isolate 08ACD0030 (its most similar hit), type strain 630 and isolate P12, a representative of the porcine isolates that contains Tet P (Fig. 3a, b). The majority of the Tet P containing contigs displayed close to 100% identity to a region within the 08ACD0030 draft genome (Fig. 3b), which was predicted to be chromosomal through a comparison with the closed 630 genome . In the porcine Tet P containing isolates an insertion has occurred that contains Tet P, its upstream regulatory region and pcw303 (Fig. 3a, Additional file 1: Figure S1). The same insertion was also found in the human Tet P containing isolates, however these genomes also contain a string of unique ORFs encoded on the same strand directly upstream of Tet P (Fig. 3a). Of particular interest is the presence of an ORF directly upstream of Tet P that encodes a putative large serine recombinase, similar to that found in the P. sordellii Tet P locus (Figs. 2a and 3a). An alignment was constructed between the Tet P locus from MD46 and the P. sordellii isolate W2945. While there were differences in the arrangement of the ORFs flanking Tet P between these two loci, significant levels of nucleotide identity were seen (Fig. 3a).
The variation in the chromosomal location and arrangement of Tet P between C. difficile and P. sordellii (Figs. 2a and 3a) suggests that a number of recombination events have taken place in the evolution of this element. Based on the homology of these sequences between C. difficile and P. sordellii these events may have occurred across species boundaries. Site-specific recombination events usually rely on or are marked by the presence of nucleotide repeats flanking the elements in question . While no conserved direct or indirect repeat could be identified among the P. sordellii sequences, one repeat sequence was found to be present across the C. difficile genomes examined here. An 11 bp sequence was found on the genome of 08ACD0030/630 at the site where the unique regions flanking the Tet P loci are present (green lines, Fig. 3a). This sequence was also found directly flanking Tet P in both the porcine and human isolates, as well as flanking the string of ORFs upstream of Tet P in the human isolates (Fig. 3a). This repeat sequence may represent an insertion site for site specific recombination, being duplicated in the process in Tet P containing isolates, and these events may be catalyzed by the large serine recombinase associated with Tet P. These repeats may also serve as sites for RecA-mediated homologous recombination , which could also explain the insertion, deletion and rearrangement of genes at this locus. The direct repeat can also be found flanking Tet P within P. sordellii W2945 and provides a potential mechanism by which this element could have been excised from the chromosome prior to lateral gene transfer to another species. It should be noted that while not identical, the overall genomic context of the Tet P chromosomal loci appears to be conserved between P. sordellii and C. difficile, particularly in the region upstream of Tet P (Fig. 3b). This region of the chromosome may therefore represent a recombination hotspot within these species.
Although the genome sequences of the three porcine C. difficile isolates contain a Tet P locus, they were reported as being phenotypically sensitive to tetracycline . We determined the MIC for tetracycline and minocycline for the two human Tet P+ C. difficile isolates, in comparison to TetR and TetS C. difficile controls (Table 1). The Clinical and Laboratory Standards Institute (CLSI) guidelines state that 16 mg/mL is the breakpoint MIC for tetracycline resistance in C. difficile , which renders isolates MCD43 and MCD46 as tetracycline sensitive (Table 1). However, the MIC values obtained for both tetracycline and minocycline for these isolates were substantially higher than those for the tetracycline sensitive control (Table 1), which may indicate that Tet P leads to an intermediate level of tetracycline resistance in C. difficile. While the clinical outcomes of such intermediate resistance in C. difficile are currently unknown, similar levels of resistance in other organisms, such as vancomycin intermediate-resistance in Staphylococcus aureus, [33, 34] have resulted in clinical challenges in the treatment of this organism. For S. aureus, both in vitro and clinical evidence indicates that isolates with intermediate-levels of vancomycin resistance may lead to treatment failure, potentially through the development of a higher MIC during treatment . Members of the tetracycline family have been suggested for use in the prevention of C. difficile infection [35, 36]. However, our study and others have shown that determinants that provide full or intermediate resistance to tetracyclines are present in C. difficile strains , which should be considered or assessed before patient treatment commences.
This study has demonstrated that ~ 20–30% of P. sordellii isolates tested in this study are resistant to the tetracycline family of antibiotics, including doxycycline, which has strong clinical relevance for P. sordellii-mediated intrauterine infections . It should be noted that due to the rarity of P. sordellii infections the sample size analysed in this study is relatively small. The collection of further, more recent isolates is required to conclusively determine the levels of tetracycline resistance within P. sordellii. The determinant responsible for this resistance, Tet P, is present in highly variable genomic locations both within P. sordellii isolates and also among other clostridia. This variability in Tet P location and the potential association of this determinant with mobile elements in P. sordellii and C. difficile, and the location of Tet P on the conjugative plasmid pCW3 in C. perfringens [17, 18], suggests that Tet P has the capacity to spread rapidly among the clostridia. The dissemination of Tet P is not just a human clinical issue since we have shown that animal isolates of P. sordellii and C. difficile carry Tet P loci that are closely related to those of human isolates. Three related problems need to be understood when considering the dissemination of the Tet P determinant amongst clostridial strains: human health issues, animal health issues and the factor common to both of these problems, antibiotic use. The findings presented here reinforce the importance of a ‘one health’ approach in considering the use of antibiotics and the importance of genetically and phenotypically understanding the strains that carry resistance determinants.
Additional file 2: Table S2 lists the bacterial isolates used in this study. Unless stated otherwise, bacterial cultures were grown in HIS broth (37 g/L Heart Infusion (Oxoid), 5 g/L yeast extract, 1 g/L L-Cysteine, 0.375% (w/v) glucose) or on HIS agar (HIS broth with 15 g/L agar) at 37 °C in an anaerobic chamber (Coy Laboratory Products, Inc.) in an atmosphere of 10% H2, 10% CO2 and 80% N2. When required, antibiotics were included in the media at the following concentrations unless otherwise specified: tetracycline (Tet) 10 mg/L, streptomycin (Str) 200 mg/L, rifampicin (Rif) 25 mg/L and D-cycloserine (DCy) 250 mg/L.
Minimum inhibitory concentration assays
To determine the tetracycline resistance profile of the P. sordellii isolates, the following assay, based on a previously published method , was performed. Note that C. perfringens isolates were included as controls. Isolates were grown at 37 °C overnight in 20 mL HIS broth with DCy under anaerobic conditions and the turbidity of the cultures determined by optical density at a wavelength of 600 nm (OD600). The cultures were diluted in 20 mL HIS broth with DCy to an OD600 of 0.1 and grown to a late-exponential growth phase, before dilution to an OD600 of 0.3 in fresh HIS broth with DCy. Duplicate 200 μl aliquots were added to 50 μl of HIS broth containing tetracycline (32 to 0.125 mg/L final volume), minocycline or doxycycline (12 to 0.047 mg/L final volume) or no antibiotic, each with an identical concentration of antibiotic solvent and DCy, within wells of a 96-well tissue culture plate. The trays were incubated at 37 °C for ~ 20 h under anaerobic conditions.
To determine a clinically relevant MIC for the C. difficile isolates, a broth microdilution method in standing with the CLSI guidelines and based on a previously published method , was used. An inoculum of 1-2 × 106 cfu/ml from an actively growing HIS broth culture of C. difficile was diluted two fold in 96-well plates containing 100 μl of HIS broth containing tetracycline (32 to 0.125 mg/L final volume), minocycline or doxycycline (12 to 0.047 mg/L final volume) or no antibiotic, each with an identical concentration of antibiotic solvent. Plates were incubated at 37 °C for 24 h under anaerobic conditions. The MIC was determined to be the lowest concentration of the antibiotic for which no growth was observed.
Genomic DNA extractions, PCR and restriction digestions were all performed as previously described . Transformation of C. perfringens cells was carried out as previously described . DNA for Southern hybridization was subjected to electrophoresis on a 0.8% Megabase agarose gel prior to transfer. Southern hybridization was carried out using tetA(P), tetB(P) and sdl specific DIG-labelled probes, as previously described . Primers used for PCR in this study are listed in Additional file 8: Table S3.
Genome sequence analysis
Sequence homology between loci was determined, and resulting graphics were produced, using EasyFig . Multiple sequence alignments were conducted using Clustal Omega from EMBL-EBI tools . Core SNP phylogeny, antibiotic resistance gene identification and pan genome analysis of P. sordellii sequencing data (see Additional file 2: Table S2 for accession numbers) was performed using the Nullarbor pipeline . A contig containing Tet P was obtained from the draft genome sequence of P. sordellii W2922 produced as part of the Nullarbor pipeline . This full contig was not present in the original draft genome of W2922 (CELK00000000) but is present when the original readset (ERR197363) from which the published assembly was derived is assembled with MegaHit version 1.1.3. This contig was therefore uploaded to the Third Party Annotation database under accession BK010701. Genome sequences not acquired as part of this study were obtained from GenBank , with the primary accession numbers for each sequence listed in Additional file 2: Table S2. P. sordellii contigs representing pCSTP1 from R28058 (CEKZ01000007) and pCSTP2 from SSCC37615 (CDNO01000006) were obtained from GenBank , and were trimmed to represent the closed plasmids. These plasmid sequences were re-annotated and deposited into the Third Party Annotation database under accession numbers BK010450 (pCSTC1) and BK010451 (pCSTC2).
C. difficile genomic DNA was sequenced using the Illumina NextSeq500 platform with 2×150 bp PE reads. Megahit v1.1.3 , as part of the Nullarbor pipeline , was used for contig assembly. Abricate v0.2 (Melbourne Bioinformatics) as part of the Nullarbor pipeline , was used for identifying TetR genes among the C. difficile draft genomes. Annotation of sequences was performed using Prokka , followed by further manual annotation using Artemis , based on sequence identity to characterised genes or hits in the Conserved Domain Database or the Phyre2 server [24, 25]. The contigs containing Tet P from C. difficile MCD43 (MH041492) and MCD46 (MH041493) were deposited via BankIt to GenBank .
Clinical and Laboratory Standards Institute
minimum inhibitory concentration
polymerase chain reaction
ribosomal protection protein
- TetR :
- TetS :
Chopra I, Roberts M. Tetracycline antibiotics: mode of action, applications, molecular biology, and epidemiology of bacterial resistance. Microbiol Mol Biol Rev. 2001;65(2):232–60.
Roberts M. Tetracycline resistance determinants: mechanisms of action, regulation of expression, genetic mobility, and distribution. FEMS Microbiol Rev. 1996;19(1):1–24.
Griffin MO, Fricovsky E, Ceballos G, Villarreal F. Tetracyclines: a pleitropic family of compounds with promising therapeutic properties. Review of the literature. Am J Physiol Cell Physiol. 2010;299(3):C539–48.
Thaker M, Spanogiannopoulos P, Wright GD. The tetracycline resistome. Cell Mol Life Sci. 2010;67(3):419–31.
Taylor DE, Chau A. Tetracycline resistance mediated by ribosomal protection. Antimicrob Agents Chemother. 1996;40(1):1–5.
Sloan J, McMurry LM, Lyras D, Levy SB, Rood JI. The Clostridium perfringens Tet P determinant comprises two overlapping genes: tetA(P), which mediates active tetracycline efflux, and tetB(P), which is related to the ribosomal protection family of tetracycline-resistance determinants. Mol Microbiol. 1994;11(2):403–15.
WHO. Critically important antimicrobials for human medicine : ranking of antimicrobial agents for risk management of antimicrobial resistance due to non-human use. In: Advisory Group on Integrated Surveillance of Antimicrobial Resistance; 2017.
Sloan B, Scheinfeld N. The use and safety of doxycycline hyclate and other second-generation tetracyclines. Expert Opin Drug Saf. 2008;7(5):571–7.
Vidor C, Awad M, Lyras D. Antibiotic resistance, virulence factors and genetics of Clostridium sordellii. Res Microbiol. 2015;166(4):368–74.
Aldape MJ, Bryant AE, Stevens DL. Clostridium sordellii infection: epidemiology, clinical findings, and current perspectives on diagnosis and treatment. Clin Infect Dis. 2006;43(11):1436–46.
Aronoff DM, Hao Y, Chung J, Coleman N, Lewis C, Peres CM, Serezani CH, Chen GH, Flamand N, Brock TG, et al. Misoprostol impairs female reproductive tract innate immunity against Clostridium sordellii. J Immunol. 2008;180(12):8222–30.
Fischer M, Bhatnagar J, Guarner J, Reagan S, Hacker JK, Van Meter SH, Poukens V, Whiteman DB, Iton A, Cheung M, et al. Fatal toxic shock syndrome associated with Clostridium sordellii after medical abortion. N Engl J Med. 2005;353(22):2352–60.
Lee E, Sheth C, Estrada O, Eshaghian P. A classic case of toxic shock syndrome due to a not so classic organism, Clostridium sordellii. Chest. 2017;152(4, Supplement):A257.
Sasaki Y, Yamamoto K, Tamura Y, Takahashi T. Tetracycline-resistance genes of Clostridium perfringens, Clostridium septicum and Clostridium sordellii isolated from cattle affected with malignant edema. Vet Microbiol. 2001;83(1):61–9.
Brazier JS, Levett PN, Stannard AJ. Antibiotic susceptibility of clinical isolates of clostridia. J Antimicrob Chemother. 1985;15(2):181–5.
Knight DR, Squire MM, Collins DA, Riley TV. Genome analysis of Clostridium difficile PCR Ribotype 014 lineage in Australian pigs and humans reveals a diverse genetic repertoire and signatures of long-range interspecies transmission. Front Microbiol. 2017;7(1):2138.
Abraham LJ, Rood JI. Molecular analysis of transferable tetracycline resistance plasmids from Clostridium perfringens. J Bacteriol. 1985;161(2):636–40.
Bannam TL, Yan X-X, Harrison PF, Seemann T, Keyburn AL, Stubenrauch C, Weeramantri LH, Cheung JK, McClane BA, Boyce JD, et al. Necrotic enteritis-derived Clostridium perfringens strain with three closely related independently conjugative toxin and antibiotic resistance plasmids. mBio. 2011;2(5):e00190–11.
Johanesen PA, Lyras D, Rood JI. Induction of pCW3-encoded tetracycline resistance in Clostridium perfringens involves a host-encoded factor. Plasmid. 2001;46(3):229–32.
Bannam TL, Teng WL, Bulach D, Lyras D, Rood JI. Functional identification of conjugation and replication regions of the tetracycline resistance plasmid pCW3 from Clostridium perfringens. J Bacteriol. 2006;188(13):4942–51.
Couchman EC, Browne HP, Dunn M, Lawley TD, Songer JG, Hall V, Petrovska L, Vidor C, Awad M, Lyras D, et al. Clostridium sordellii genome analysis reveals plasmid localized toxin genes encoded within pathogenicity loci. BMC Genomics. 2015;16(1):392.
Chong E, Winikoff B, Charles D, Agnew K, Prentice JL, Limbago BM, Platais I, Louie K, Jones HE, Shannon C. Vaginal and rectal Clostridium sordellii and Clostridium perfringens presence among women in the United States. Obstet Gynecol. 2016;127(2):360–8.
Johanesen PA, Lyras D, Bannam TL, Rood JI. Transcriptional analysis of the tet(P) operon from Clostridium perfringens. J Bacteriol. 2001;183(24):7110–9.
Marchler-Bauer A, Derbyshire MK, Gonzales NR, Lu S, Chitsaz F, Geer LY, Geer RC, He J, Gwadz M, Hurwitz DI, et al. CDD: NCBI's conserved domain database. Nucleic Acids Res. 2015;43(Database issue):D222–6.
Kelley LA, Mezulis S, Yates CM, Wass MN, Sternberg MJE. The Phyre2 web portal for protein modeling, prediction and analysis. Nat Protocols. 2015;10(6):845–58.
Garnier T, Cole ST. Identification and molecular genetic analysis of replication functions of the bacteriocinogenic plasmid pIP404 from Clostridium perfringens. Plasmid. 1988;19(2):151–60.
Rood JI, Cole ST. Molecular genetics and pathogenesis of Clostridium perfringens. Microbiol Rev. 1991;55(4):621–48.
Garcillán-Barcia MP, Francia MV, De La Cruz F. The diversity of conjugative relaxases and its application in plasmid classification. FEMS Microbiol Rev. 2009;33(3):657–87.
Riedel T, Bunk B, Thürmer A, Spröer C, Brzuszkiewicz E, Abt B, Gronow S, Liesegang H, Daniel R, Overmann J. Genome resequencing of the virulent and multidrug-resistant reference strain Clostridium difficile 630. Genome Announcements. 2015;3(2):e00276–15.
Frost LS, Leplae R, Summers AO, Toussaint A. Mobile genetic elements: the agents of open source evolution. Nat Rev Microbiol. 2005;3:722.
Rocha EPC, Cornet E, Michel B. Comparative and evolutionary analysis of the bacterial homologous recombination systems. PLoS Genet. 2005;1(2):e15.
CLSI. Methods for antimicrobial susceptibility testing of anaerobic bacteria; approved standard - seventh edition. CLSI document M11-A7. Wayne, PA, USA: Clinical and Laboratory Standards Institute; 2012.
Howden BP, Davies JK, Johnson PDR, Stinear TP, Grayson ML. Reduced vancomycin susceptibility in Staphylococcus aureus, including vancomycin-intermediate and heterogeneous vancomycin-intermediate strains: resistance mechanisms, laboratory detection, and clinical implications. Clin Microbiol Rev. 2010;23(1):99–139.
Wang Y, Li X, Jiang L, Han W, Xie X, Jin Y, He X, Wu R. Novel mutation sites in the development of vancomycin- intermediate resistance in Staphylococcus aureus. Front Microbiol. 2017;7(1):2163.
Khanna S. Do tetracyclines have the potential to reduce the risk of Clostridium difficile infection? Expert Rev Anti Infect Ther. 2018:16(3):183–185.
Doernberg SB, Winston LG, Deck DH, Chambers HF. Does doxycycline protect against development of Clostridium difficile infection? Clin Infect Dis. 2012;55(5):615–20.
Huang H, Weintraub A, Fang H, Nord CE. Antimicrobial resistance in Clostridium difficile. Int J Antimicrob Agents. 2009;34(6):516–22.
Wiegand I, Hilpert K, Hancock REW. Agar and broth dilution methods to determine the minimal inhibitory concentration (MIC) of antimicrobial substances. Nat Protocols. 2008;3(2):163–75.
Vidor CJ, Watts TD, Adams V, Bulach D, Couchman E, Rood JI, Fairweather NF, Awad M, Lyras D. Clostridium sordellii pathogenicity locus plasmid pCS1-1 encodes a novel Clostridial conjugation locus. mBio. 2018;9(1):e01761–17.
Wisniewski JA, Teng WL, Bannam TL, Rood JI. Two novel membrane proteins, TcpD and TcpE, are essential for conjugative transfer of pCW3 in Clostridium perfringens. J Bacteriol. 2015;197(4):774–81.
Rood JI. Transferable tetracycline resistance in Clostridium perfringens strains of porcine origin. Can J Microbiol. 1983;29(10):1241–6.
Rood JI, Maher EA, Somers EB, Campos E, Duncan CL. Isolation and characterization of multiply antibiotic-resistant Clostridium perfringens strains from porcine feces. Antimicrob Agents Chemother. 1978;13(5):871–80.
Sullivan MJ, Petty NK, Beatson SA. Easyfig: a genome comparison visualizer. Bioinformatics. 2011;27(7):1009–10.
McWilliam H, Li W, Uludag M, Squizzato S, Park YM, Buso N, Cowley AP, Lopez R. Analysis tool web services from the EMBL-EBI. Nucleic Acids Res. 2013;41(Web Server issue):W597–600.
Nullarbor [https://github.com/tseemann/nullarbor]. Accessed 25 May 2018.
Clark K, Karsch-Mizrachi I, Lipman DJ, Ostell J, Sayers EW. GenBank. Nucleic Acids Res. 2016;44(D1):D67–72.
Li D, Liu CM, Luo R, Sadakane K, Lam TW. MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics. 2015;31(10):1674–6.
Seemann T. Prokka: rapid prokaryotic genome annotation. Bioinformatics. 2014;30(14):2068–9.
Rutherford K, Parkhill J, Crook J, Horsnell T, Rice P, Rajandream MA, Barrell B. Artemis: sequence visualization and annotation. Bioinformatics. 2000;16(10):944–5.
We thank Neil Fairweather from Imperial College for providing P. sordellii isolates. We thank Allen C. Cheng and Anton Y. Peleg from the Infectious Diseases Unit at the Alfred Hospital for the provision of C. difficile isolates and access to their genome sequence data. We thank Jacob Amy from the Monash Biomedical Discovery Institute (BDI) for providing draft C. difficile genome sequences and his help with the Nullarbor pipeline. We thank David Powell and Adele Barugahare from the Monash Bioinformatics Platform for preliminary sequence data handling. We thank Tom Stent from the Monash Biomedical Discovery Institute (BDI) for his help in preliminary annotation of C. difficile genome sequences.
D.L. was supported by Future Fellowship FT120100779 (Australian Research Council). C.J.V. was supported through an Australian Government Research Training Program Scholarship. Note that funding body played no part in the design of the study; collection, analysis, and interpretation of data; or in writing of the manuscript.
Availability of data and materials
The contig containing Tet P from P. sordellii isolate W2922 was deposited into the Third Party Annotation database under accession number BK010701. The contigs containing Tet P from C. difficile MCD43 and MCD46 were deposited to GenBank under accession numbers MH041492 and MH041493, respectively. P. sordellii sequences of pCSTP1 from R28058 and pCSTP2 from SSCC37615 were deposited into the Third Party Annotation database under accession numbers BK010451 and BK010450 respectively. All other genome sequences or sequencing data used in this analysis was obtained from GenBank , with their respective accession numbers listen in Additional file 2: Table S2.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. The Tet P tetracycline resistance determinant and surrounding regions in C. perfringens. Characterized regulatory elements upstream of the Tet P determinant are annotated and include the promoter P3 and the transcriptional terminator T1. (JPG 607 kb)
Table S2. Bacterial isolates, genomic DNA sequences and sequencing data used in this study. (PDF 274 kb)
Table S1. Matrices displaying the nucleotide identity (%) of the Tet P regulatory region comprising 650 bp directly upstream of the tetA(P) start codon, as well as the amino acid identity (%) of TetA(P) and TetB(P), from isolates of P. sordellii and C. difficile (Cd) compared with that of pCW3 from C. perfringens (Cp) CW92. (PDF 167 kb)
Figure S2. Nucleotide alignment of the upstream regulatory regions of Tet P in P. sordellii and C. difficile isolates compared to that of C. perfringens pCW3 from strain CW92. The sequence from P. sordellii isolates that encode Tet P on small plasmids is represented by R28058 (but includes R15892, SSCC37615, R32977 and SSCC19939, Table S1). The sequence from P. sordellii isolates that appear to encode Tet P on the chromosome, with the exception of W2945, is represented by SSCC18392 (but includes SSCC32135, R32462, JGS6961, W3026, AM370 and W2922, Table S1). The beginning of the tetA(P) gene is colored red (red arrow) with the start codon in bold. The predicted ribosome binding sites (RBS) for these sequences are in black bold font. The region corresponding to the C. perfringens pCW3 Tet P promoter, P3, is highlighted in green with the predicted − 10 and − 35 boxes shown in bold. The inverted repeats that correspond to the T1 transcriptional terminator are highlighted in blue and indicated with inverted black arrows. (PDF 274 kb)
Figure S3. Nucleotide alignment of the regions flanking regA in P. sordellii Tet P+ isolates. Sequences for the regA+ isolates SSCC32125, R32462, JGS6961 and SSCC18392 are compared to the regA− isolate W2945. The direct repeats flanking regA are shown in bold red font and indicated by a black arrow beneath the sequence. The coding region of regA is colored purple. A gap in the alignment representing the majority of regA is indicated by three large dots. The start and stop codons of regA are in bold and annotated on the sequence. The regA− isolates W3026, AM370 and W2922 have not been included in this alignment; however, a single site identical to the repeat is present for this isolate, equidistant to Tet P, as for W2945. (PDF 238 kb)
Figure S4. Comparison of the genomic location of tetA(P) and tetB(P) among P. sordellii isolates using Southern hybridization analysis. Genomic DNA from C. perfringens isolate JIR56 was used as a tetA/B(P) positive control and sdl negative control, and genomic DNA from the P. sordellii isolates was used as sdl positive controls. P. sordellii isolate SSCC33587 was included as a Tet P negative control. Genomic DNA from P. sordellii isolates and C. perfringens strain JIR56 was either undigested (U) or digested with EcoRI (D). Blots were probed with DIG-labelled DNA either specific for sdl (chromosomal P. sordellii marker), tetA(P) or tetB(P). In the chromosomal Tet P isolates, both tetA(P) and tetB(P) probes hybridized to a band which represented undigested chromosomal DNA. This chromosomal location was confirmed using a sdl (chromosomal marker)-specific probe, which produced a fragment of a similar size to that seen with tetA(P) and tetB(P). For the undigested genomic DNA of the remaining 3 TetR isolates, distinct lower-molecular weight bands were observed for both the tetA(P) and tetB(P) blots. When digested with EcoRI and probed with tetA(P) or tetB(P), pCSTC1 and pCSTC2-carrying isolates produced bands of ~ 8.5 kb and ~ 9.8 kb, respectively, which are the expected sizes of these plasmids. (JPG 2793 kb)
Figure S5. PCR analysis of C. perfringens JIR325 transformed with P. sordellii genomic DNA (gDNA) from isolates SSCC37615 or SSCC18838. Three separate isolates from each transformation were analysed (Lanes 1–6). PCR reactions were subjected to agarose gel electrophoresis against a Hyperladder 1 kb marker (Bioline), with relevant sizes indicated in bp. (a) PCR to detect the Tet P determinant using the internal tetA(P) primer DLP104 and the internal tetB(P) primer DLP105. A product of the expected size was observed for all reactions with the exception of non-transformed JIR325 (negative control) and the no DNA control. (b) PCR to detect the C. perfringens chromosomal gene plc using internal primers JRP2873 and JRP2874. A product of the expected size was observed for all reactions with the exception of P. sordellii (Ps) isolate R28058 (negative control) and the no DNA control. (JPG 58 kb)
Table S3. Oligonucleotide primers used in PCR. (+) forward primer, (−) reverse primer. (PDF 89 kb)