Various pAQU plasmids possibly contribute to disseminate tetracycline resistance gene tet(M) among marine bacterial community

Emergence of antibiotic-resistant bacteria in the aquaculture environment is a significant problem for disease control of cultured fish as well as in human public health. Conjugative mobile genetic elements (MGEs) are involved in dissemination of antibiotic resistance genes (ARGs) among marine bacteria. In the present study, we first designed a PCR targeting traI gene encoding essential relaxase for conjugation. By this new PCR, we demonstrated that five of 83 strains isolated from a coastal aquaculture site had traI-positive MGEs. While one of the five strains that belonged to Shewanella sp. was shown to have an integrative conjugative element of the SXT/R391 family (ICEVchMex-like), the MGEs of the other four strains of Vibrio spp. were shown to have the backbone structure similar to that of previously described in pAQU1. The backbone structure shared by the pAQU1-like plasmids in the four strains corresponded to a ~100-kbp highly conserved region required for replication, partition and conjugative transfer, suggesting that these plasmids constituted “pAQU group.” The pAQU group plasmids were shown to be capable of conjugative transfer of tet(M) and other ARGs from the Vibrio strains to E. coli. The pAQU group plasmid in one of the examined strains was designated as pAQU2, and its complete nucleotide sequence was determined and compared with that of pAQU1. The results revealed that pAQU2 contained fewer ARGs than pAQU1 did, and most of the ARGs in both of these plasmids were located in the similar region where multiple transposases were found, suggesting that the ARGs were introduced by several events of DNA transposition into an ancestral plasmid followed by drug selection in the aquaculture site. The results of the present study indicate that the “pAQU group” plasmids may play an important role in dissemination of ARGs in the marine environment.


INTRODUCTION
Wide ranges of antimicrobial agents have been used not only for treatment of human infectious diseases but also in the fish farming industry for therapeutic and prophylactic purposes (Cabello et al., 2013). Among these, tetracyclines (TCs) are broadspectrum agents, exhibiting activity against various microorganisms, such as Gram-positive and -negative bacteria, chlamydiae, mycoplasmas, and rickettsiae. In some countries, TCs are added to animal feeds at subtherapeutic levels so that they can act as growth promoters (Chopra and Roberts, 2001). In Japan, oxytetracyline and doxycyline have been included in the drugs permitted to use in the process of aquaculture, and the former is one of the drugs available to use for the largest variety of fish species (Ministry of Agriculture, Forestry and Fisheries, Japan, 2011). When the antibiotics are used in the aquaculture, they are usually mixed with feed and administered to fish, and the extensive use of TCs has led emergence of resistant bacteria (Nonaka et al., 2007). To date, 46 different TC resistance genes have been submitted to the nomenclature center (http://faculty.washington.edu/ marilynr/), of which 30 genes code for efflux proteins, 12 genes for ribosomal protection proteins and four genes for other categories. Emergence and dissemination of drug resistance in the marine environment have a large impact on the fish farming industry. It is also an important issue for clinical medicine because a number of antibiotics resistance genes (ARGs) are shared between marine bacteria and human pathogens (Cabello et al., 2013) and the gene flow between these environments has been suggested. Further, transfer of opportunistic human pathogen from aquatic environments possible to generate antimicrobial resistance infections (Cabello et al., 2013).
A TC resistance gene, tet(M), is among the ARGs shared between marine bacteria and human pathogens and has extensively been studied in terms of the molecular ecology in the environment. The gene is known to be distributed among numerous bacterial species (Chopra and Roberts, 2001;Roberts, 2005) derived from diverse sources, such as human pathogens (Jones et al., 2006), domestic animals Chee-Sanford et al., 2001), food (Gevers et al., 2003;Bertrand et al., 2005), and farmed fish (Kim et al., 2007). In our previous study (Kim et al., 2004;Nonaka et al., 2007), we also found that tet (M) was distributed among several bacterial species isolated from fishes or the marine aquaculture environment. Other recent studies have also reported identification of tet(M) in the bacteria from the aquaculture environments (Nikolakopoulou et al., 2008;Seyfried et al., 2010;Di Cesare et al., 2012. Thus, the marine bacteria have been suggested to be a reservoir of tet (M).
While detailed mechanisms by which tet(M) has propagated in the marine environment are unknown, it has previously been demonstrated that tet(M) is frequently carried by mobile genetic elements (MGEs), such as plasmid and transposon (Chopra and Roberts, 2001). Previously, 73 strains of tet(M)-positive marine bacteria from an aquaculture site were examined, and it was shown that five of them, strains 04Ya001, 04Ya016, 04Ya090, 04Ya108, and 04Ya311, were able to transfer TC resistance to E. coli (Neela et al., 2009). Among these strains, 04Ya311 was characterized in detail, and it was demonstrated that this strain had a novel self-transferable plasmid, which we designated as pAQU1, and that conjugative transfer of pAQU1 conferred TC resistance on E. coli (Nonaka et al., 2012). Determination of its complete nucleotide sequence revealed that pAQU1 encoded tet(M) and other ARGs, such as tet(B), sul2, floR and genes similar to bla CARB−9 , mph(A) and mef (A). Phylogenetic analysis of pAQU1's traI gene encoding relaxase (Nonaka et al., 2012), an only protein component of the conjugative machinery common to all MGEs using a conjugative system (Garcillan-Barcia et al., 2009;Smillie et al., 2010), indicated that the plasmid had a close relationship with IncA/C plasmid and SXT/R391 family of integrative conjugative elements (SXT/R391 ICEs) which belong to the MOB H12 clade in the MOB H family (Nonaka et al., 2012). Since IncA/C plasmid and SXT/R391 ICEs are distributed among many aquatic isolates of genera Enterobacteriaceae and Vibrionaceae (Welch et al., 2007;McIntosh et al., 2008;Fricke et al., 2009), we hypothesized that pAQU1 found in 04Ya311 was a novel plasmid that might play a role in dissemination of tet(M) in the marine environment. In present study, since several pAQU1-related plasmids were found, we aimed to reveal that these pAQU1-related plasmids constituted a novel "pAQU group" characterized by a common backbone structure and conferred to transfer tet(M) and other ARGs to E. coli by conjugation.

BACTERIAL STRAINS
For screening of traI, 83 tetracycline-resistant strains of bacteria isolated from the sediment or seawater of a coastal aquaculture site in Japan were used (Supplementary Table1). At the site, oxytetracycline and flumequine had been used (Nonaka et al., 2007). Theses environmental strains were cultured at 25 • C in the brain heart infusion (BHI) medium (Becton, Dickinson and Company, Sparks, MD, USA) with 0.2% NaCl and 10 μg/ml of TC (Nakalai tesque, Kyoto, Japan). All of the strains were confirmed by PCR to have tet(M) (Nonaka et al., 2007), and strains 04Ya001, 04Ya016, 04Ya090, and 04Ya108 have previously been shown to transfer tet(M) to E. coli (Neela et al., 2009). E. coli W3110 obtained from the National BioResource Project (National Institute of Genetics, Japan) and its rifampicinresistant derivative, W3110Rif r (Nonaka et al., 2012), were cultured at 37 • C in the Luria-Bertani medium (BD) with or without 10 μg/ml (for liquid media) or 20 μg/ml (solid media) of TC and 100 μg/ml of rifampicin (Sigma-Aldrich, Saint-Louis, MO).

DNA EXTRACTION AND PCR
Total bacterial DNA was extracted by using a QuickGene DNA tissue kit (Fujifilm Corporation, Tokyo, Japan) according to the manufacturer's protocol.
For PCR detection of ARGs [tet(B), tet(M), sul2, floR, bla CARB−9 -like, mph(A)-like and mef (A)-like], genes related to the type IV secretion system (repA, traA2, traB, traC, traF, traN, traI, and traL) and ICEVchMex-related genes (int, traF, and traN), respective primers shown in Table 1 were used. PCR amplification was performed with initial denaturation at 94 • C for 1 min followed by 30 cycles of denaturation at 94 • C for 30 s, 1 min of annealing at the temperatures shown in Table 1, and extension at 72 • C for 1 min and a final extension step at 72 • C for 3 min.
BIO INC) and the template DNA (20-50 ng). For amplicons of 2-10 kb, PCR amplification was performed with initial denaturation at 94 • C for 1 min followed by 30 cycles of denaturation at 94 • C for 30 s, annealing at 55 • C for 1 min, and extension at 72 • C for 5 min and a final extension step at 72 • C for 3 min. For amplicons longer than 10 kb, two-step PCR was performed with 30 cycles of denaturation at 96 • C for 20 s and extension at 69 • C for 16 min.  For determination of the phylogenetic position of the strains were identified by the nucleotide sequences of the bacterial 16S rRNA gene which were amplified using primers 9F (5 -GAGTTTGATCCTGGCTCAG-3 ) and 1510R (5 -GGCTACCTTGTTACGA-3 ). PCR reaction was performed as previously described (Nonaka et al., 2007).

FILTER MATING
The donor marine bacteria and recipient (E. coli W3110) were cultured to the mid-log phase (OD600 = 0.5-0.8) and mixed with each other (500 μl each) in triplicates. The mixture was then collected onto a sterile 25 mm in a diameter, 0.2 μm-pore size nitrocellulose membrane filter (Merck millipore, Billerica, MA) by using a KG-25 filter holder (Advantec, Tokyo, Japan) and a vacuum pump at 20 kPa. The filter with the collected cells were rinsed by 10 ml of BHI medium with 2% NaCl or LB and placed on the Marine Broth 2216 plate (Becton, Dickinson and Company) for incubation at 25 • C for 20 h. Each filter was then transferred into 500 μl of phosphate-buffered saline (PBS) and agitated vigorously with a Vortex mixer to detach the bacteria. The suspension was serially diluted by 10-fold in PBS and plated onto the media. To select for the transconjugants, LB agar plates with 20 μg/ml of TC were used and incubated for 20 h at the temperature (42 • C) which inhibits growth of the donors. For counting the total number of donor cells and the transconjugants, the cell suspension was spread on a BHI agar plate with 0.2% NaCl and 10 μg/ml of TC and incubated at 25 • C. When the transconjugants were used as donors, E. coli W3110Rif r was used as a recipient and 2nd transconjugants were counted on LB agar plates with 20 μg/ml of tetracycline and 100 μg/ml of rifampicin after incubation at 37 • C for 20 h. For counting the total number of the donor cells and the transconjugants, LB plates with 20 μg/ml tetracycline were used and incubated at 37 • C for 20 h. The number of the donor cells was obtained by subtracting the number of the transconjugants from the total number of the donor cells and the transconjugants. Transfer rate was calculated as a ratio of the number of transconjugants per donor cell.

PULSED-FIELD GEL ELECTROPHORESIS (PFGE) AND SOUTHERN HYBRIDIZATION
To determine whether tet(M) is carried by the plasmid or integrated in the chromosome, PFGE and Southern hybridization were performed. First, a DNA-containing agar plug of each sample was prepared by the previously reported method (Kijima-Tanaka et al., 2007) with the following modifications. Culture fluid containing about 10 8 cells was put into a 15-ml centrifuge tube and centrifuged at 9,000 rpm for 10 min. The cell pellet was re-suspended in 120 μl of cell suspension buffer [100 mM Tris-HCl, 100 mM EDTA (pH 8.0)] and mixed with 120 μl of melted 2% InCert® Agarose (Lonza, Basel, Switzerland) and 12 μl of lysozyme (20 mg/ml; Sigma-Aldrich) dissolved in TE buffer [10 mM Tris-HCl, 1 mM EDTA (pH 8.0)]. The mixture was quickly dispensed into the plug mold (Bio-Rad, Hercules, CA). After solidification, the plug was transferred into 500 μl of TE buffer containing 1% SDS and 12 μl of lysozyme (20 mg/ml), and incubated at 37 • C for 1 h. Then, the liquid portion in each tube was replaced with 500 μl of lysis buffer [100 mM Tris-HCl, 100 mM EDTA (pH 8.0), 1% N-lauroylsarcosine] and 50 μl of proteinase K (20 mg/ml; TAKARA BIO INC.), and the plug was incubated at 50 • C for overnight. On the next day, the plug was washed three times with 500 μl of TE buffer at room temperature for 15 min, pretreated with 1 ml of S1 nuclease buffer [30 mM sodium acetate (pH 4.6), 1 mM zinc acetate, 5%(v/v) glycerol] at room temperature for 30 min, and treated with 10 U of S1 nuclease (Life Technologies Corporation, Carlsbad, CA) in 200 μl of the S1 nuclease buffer at 37 • C for 15 min. The plugs was then washed twice with 1 ml of ice-cold high concentration EDTA-Tris [10 mM Tris-HCl, 100 mM EDTA (pH 8.0)]. The plug was cut into half and one of the pieces was embedded in a 1% Seakem® Gold Agarose (Lonza) gel and PFGE was performed on a CHFE MAPPER™ (Bio-Rad) under the conditions described previously (Nonaka et al., 2012). The DNA size standard lambda ladder (Bio-Rad) was used as a size marker. The DNA in the gel was stained with SYBR Gold (Life Technologies Corporation) at the concentration recommended by the manufacturer and visualized under UV light. Subsequent hybridization was performed according to the standard method (Sambrook and Russell, 2001). PCR-amplified DNA fragments of tet(M), traI in pAQU1 (traI-pAQU1), traI in ICEVchMex (traI-ICEVchMex) and rep in pAQU1 (rep-pAQU1) were labeled with digoxigenin by using the primers shown in Table 1 and a PCR DIG probe synthesis kit (Roche Diagnostics, Mannheim, Germany). Hybridizations were performed at 41 • C for overnight.

SEQUENCING AND ANNOTATION
The entire sequence of the total DNA of TJ090W1 was determined by the pyrosequencing method with the Titanium chemistry and the Genome Sequence FLX system (Roche Diagnostics), and a set of data consisting of 225,251 shotgun and 271,612 pairedend (average insert size: 7 kb) sequences with an average length of 684 and 398 nucleotides were generated, respectively. The volume of the data was estimated to cover 60 times of the entire DNA sequence. The plasmid-derived sequences were obtained by subtracting the genomic sequence of E. coli W3110 from the entire set of data and then assembled de novo by using the GS assembler (Newbler 2.8) (Roche Diagnostics), generating 1 scaffold and 19 sets of assembled contigs (over 500 bp), of which five selected contigs were used for assembling the whole structure of the plasmid molecule according to the previously described method (Nonaka et al., 2012). Preliminary identification of protein coding regions in the plasmid was carried out by using MetaGeneAnnotator (Noguchi et al., 2008). Then, the results were subjected to BLAST searches and compared with the sequences listed in the NCBI non-redundant protein database for their length, identity and coverage, and final annotation was manually performed by using In Silico Molecular Cloning Genomics Edition® (In Silico Biology, Inc., Yokohama, Japan) with the threshold E-value of 10 −5 .

NUCLEOTIDE SEQUENCE ACCESSION NUMBER
The entire sequence of pAQU2 was deposited in GeneBank/ DDBJ/ EMBL databases under the accession no. AB856327.

SCREENING BY PCR DETECTION AND SEQUENCING OF traI GENE FOR THE MARINE BACTERIAL STRAINS WHICH HARBOR tet (M)-BEARING pAQU1-LIKE PLASMID OR ICEVchMEX-LIKE MGE
For screening for the bacterial strains carrying the pAQU1-like plasmids, we designed degenerate primers, MOBH1_traI062-1F and MOBH1_traI062-2R, based on the traI nucleotide sequences of pAQU1 (GenBank accession number AB571865), five of IncA/C plasmids (Welch et al., 2007;Fricke et al., 2009) and seven of SXT/R391 ICEs (Wozniak et al., 2009) which belong to the MOB H12 clade of the MOB H family ( Table 1). When 83 tet(M)-positive strains of marine bacteria from the aquaculture site were examined, amplified traI DNA of 645 base pairs (bps) was detected in strains 04Ya001, 04Ya016, and 04Ya090 ( Table 3) that had been known to carry transferable tet(M) gene (Neela et al., 2009). The 645-bp traI PCR product was also detected in strains 04Ya007 and 04Ya265 whose ability to transfer tet(M) had been unclear (Table 3). Although strain 04Ya108 had been shown to conjugatively transfer tet(M) to E. coli (Neela et al., 2009), traI was not detected in this strain (Table 3). traI was not detected in 77 other strains, either; 66 strains whose tet(M) was shown not to be transferable (Neela et al., 2009) and 11 strains that had not been examined before. Therefore, five of 83 strains examined appeared to carry MGE related to pAQU1 or IncA/C plasmid or SXT/R391 ICE. Then, nucleotide sequences of the PCR-amplified 645-bp traI DNA of the plasmids in these five strains were determined. The results indicated that the traI sequence of four strains (04Ya001, 04Ya016, 04Ya090, and 04Ya265) were identical to that of pAQU1, whereas the sequence of one strain (04Ya007) was distinct from that of pAQU1, but similar to that of ICEVchMex (97% identity), which belongs to SXT/R391 ICEs found in V. cholerae isolated from sewage in Mexico (Burrus et al., 2006). The ICEVchMex-specific int, traF and traN genes were also detected by PCR in 04Ya007. Therefore, strains 04Ya001, 04Ya016, 04Ya090, and 04Ya265 were shown to have pAQU1 or closely related pAQU1-like plasmids, whereas 04Ya007 likely had another member of the MOB H12 clade, ICEVchMex-like MGE.
To examine whether tet(M) is also plasmid-borne in strains 04Ya001, 04Ya016, 04Ya090, and 04Ya265 carrying the pAQU1like plasmids, 04Ya007 with ICEVchMex-like MGE and 04Ya108 without traI, PFGE followed by Southern hybridization was performed. Hybridization with the tet(M) probe generated extrachromosomal DNA bands ranging from 150 to 350-kbp in all of the strains examined ( Figure 1A, lanes 1-5 and Figure 2A, lane 1). When the traI-pAQU1 or rep-pAQU1 probe was used, bands of the same sizes were detected for 04Ya001, 04Ya016, 04Ya090, and 04Ya265 (Figures 1B,C, lanes 1-3 and 5), but not 04Ya108 (Figures 1B,C, lane 4) compatible with the PCR results. While when a strain 04Ya007 was used, the traI-ICEVchMex probe detected a signal at the position corresponding to the chromosomal DNA ( Figure 2B, lane 1). The results indicated that in four strains 04Ya001, 04Ya016, 04Ya090, and 04Ya265, tet(M) was most likely carried on the pAQU1-like plasmids, whereas in strain 04Ya108, tet(M) appeared to be in the plasmid of non-MOB H12 clade. In strain 04Ya007, tet(M) appeared to be plasmid-borne while ICEVchMex-like MGE was integrated in the chromosome.

CONJUGATIVE TRANSFER OF tet (M) GENE MEDIATED BY THE pAQU1-LIKE PLASMIDS
To confirm and examine the ability of the six strains to transfer TC resistance, filter mating experiments were performed by using E. coli as a recipient. As reported previously (Neela et al., 2009), strains 04Ya001, 04Ya016, and 04Ya090 were able to transfer TC resistance to E. coli (Table 4). Strains 04Ya007 and 04Ya265 were also shown to transfer TC resistance to E. coli with the frequencies of 1.37 × 10 −4 and 1.08 × 10 −4 , respectively, at the levels comparable to the positive control strain 04Ya311 (Table 4). Strain 04Ya108, which did not appear to carry MOB H12 MGE, was also able to transfer TC resistance, but with a much lower frequency (Table 4).
Then, a single clone of TC-resistant transconjugant was selected for each donor strain and analyzed by PCR with MOBH1_traI062-1F and MOBH1_traI062-2R primers. The traI gene was detected in five transconjugants whose donors were 04Ya001, 04Ya007, 04Ya016, 04Ya090, and 04Ya265, but not in the transconjugant whose donor was 04Ya108 ( Table 5). Nucleotide sequence analysis of the PCR products revealed that the traI sequences derived from the transconjugants whose donors were 04Ya001, 04Ya016, 04Ya090, and 04Ya265 were identical to that of pAQU1. On the other hand, the traI sequences of transconjugant TJ007W1 was identical to that of its donor 04Ya007, with 97% identity to ICEVchMex.
These results strongly suggested that the tet(M)-bearing pAQU1-like plasmids could be conjugatively transferred from marine bacteria to E. coli, conferring TC resistance. The mechanism by which tet(M) gene was transferred from 04Ya108 to TJ108W1 is unclear.

TRANSFER OF OTHER ARGs ASSOCIATED WITH CONJUGATIVE TRANSFER OF THE pAQU1-LIKE PLASMIDS
PCR analysis revealed that 04Ya001, 04Ya016, 04Ya090, and 04Ya265 had tet(B) in addition to tet(M) ( Table 5). Strain 04Ya265 was shown by PCR to have other ARGs, such as sul2, bla CARB−9 -like, mph(A)-like and mef (A)-like (Table 5). Generally, all of the ARGs were transferred to the transconjugants when these strains, capable of conjugative transfer of the pAQU1-like plasmids, were used as donors ( Table 5).

SELF-TRANSFERABILITY OF THE pAQU1-LIKE PLASMIDS
To investigate whether the pAQU1-like plasmids are selftransferable, additional filter mating experiments were performed by using the transconjugants as donors and E. coli W3110Rif r as a donor. TJ001W1, TJ090W1 and TJ265W1, which had the tet(M)-bearing pAQU1-like plasmid, transferred TC resistance at the frequencies (3.90 × 10 −5 to 1.08 × 10 −4 ) comparable to that of the positive control TJ311W2 (Table 4). Although TJ016W1 also had a tet(M)-bearing pAQU1-like plasmid, it was unable to transfer TC resistance at a detectable level.  We also tested transconjugants TJ007W1 and TJ108W1 which did not harbor the pAQU1-like plasmid. Ten times higher frequency (1.01 × 10 −3 ) of TC resistance transfer was observed for TJ007W1 which had no detectable tet(M), whereas the transfer frequency of TJ108W1 was below the detection limit (Table 4).

CHARACTERIZATION OF A NOVEL pAQU1-LIKE PLASMID, pAQU2
As shown above, strains 04Ya001, 04Ya016, 04Ya090, and 04Ya265 and their transconjugants, TJ001W1, TJ016W1, TJ090W1, and TJ265W1, had a pAQU1-like plasmid. Among these, we chose the plasmid in TJ090W1, which was shown to be self-transferable, and determined its entire nucleotide sequence. The novel plasmid, which we designated as pAQU2, was 160,406-bp long and contained 195 predicted coding sequences (CDSs). The CDSs were asymmetrically distributed; i.e., 163 CDSs were oriented clockwise, whereas only 32 CDSs were oriented counterclockwise (Figure 3).
Comparison of the pAQU1 and pAQU2 sequences indicated that the regions containing the genes required for basic plasmid functions, such as replication, partition and conjugative transfer, were highly conserved between the two plasmids. Thus, it was suggested that these conserved regions, which totally sum up to ∼100 kbps, might constitute the backbone structure common to the pAQU1-like plasmids. In our previous study (Nonaka et al., 2012), pAQU1 was shown to have a high G+C content region with a cluster of ARGs and transposases. In pAQU2, a region with a cluster of ARGs and transposases were also found at the corresponding position, but with a low G+C content (Figure 3).

pAQU GROUP PLASMIDS CHARACTERIZED BY A COMMON BACKBONE STRUCTURE
To test whether the backbone structure shared by pAQU1 and pAQU2 also exists in the other pAQU1-like plasmids, we characterized the donor strains 04Ya001, 04Ya016, 04Ya090, and 04Ya265 and their transconjugants by PCR with 12 sets of primers covering the entire backbone structure. The results showed that the backbone structure was maintained with no apparent deletion in the pAQU1-like plasmids harbored in all of the donors and the transconjugants examined ( Table 6). Additional PCR analysis with the primers specific to some of the genes (rep, traA2, traB, traC, traF, traN and traL) in the backbone structure revealed that these genes were also retained in all of the donors and the transconjugants examined ( Table 6). These results indicated that the backbone structure common to pAQU1 and pAQU2 existed in all of the other pAQU1-like plasmids examined. Therefore, there may be a novel group (pAQU group) of plasmids characterized by this backbone structure.
Interestingly, strain 04Ya007 in which ICEVchMex-like traI, but not pAQU1-like traI, was shown by PCR to harbor part of the pAQU backbone structure and some backbone-specific genes ( Table 6). In addition, Southern hybridization with the rep-pAQU1 probe generated a band corresponding to the tet(M)bearing plasmid in this strain (Figure 2A, lane 1, Figure 2C, lane 1). Therefore, it is possible that 04Ya007 carry plasmid similar to that of the pAQU group although the plasmid did not appear to be transferred to TJ007W1 (Figures 2A,C lane 2 and Table 6). Neither the pAQU backbone structure nor the backbone-specific genes were detected in 04Ya108 or its transconjugant ( Table 6).

DISCUSSION
Recently, we discovered a novel tet(M)-bearing plasmid pAQU1 (204,052 bp) whose TraI primary structure was similar to that of the MOB H12 MGE of the MOB H family and proposed that pAQU1 might play a role in dissemination of tet(M) and other ARGs in the marine environment (Nonaka et al., 2012). In this study, in order to investigate this possibility, we first designed degenerate PCR primers for detecting traI gene of the MOB H12 MGEs including pAQU1. With these traI primers, we could FIGURE 3 | Circular representation of pAQU1 backbone and specific genes for pAQU2 and pAQU1, respectively. Circles from inside to outside: 1, Putative pAQU back bone regions, 2 and 3, shared CDSs between pAQU1 and pAQU2; 5 and 6, CDSs specific for pAQU1; and 8 and 9, CDSs specific for pAQU2. Each of CDSs inside and outside of the circle coded counter-clockwise and clockwise, respectively. The 4th and 7th circle indicates GC content of pAQU1 and pAQU2 where purple shows upper GC value above the center line and green shows lower GC value below the center line. Putative functions of the products of the CDSs are indicated in color: red, replication, partition and termination; light blue, conjugative transfer; blue, integration; green, transposition; orange, antibiotic resistance; yellow, other functions and gray, unknown function. rapidly identify the bacterial strains which have pAQU1-like plasmids or IncA/C plasmids or SXT/R391 ICEs. Compared with the filter mating usually used for identifying the strains with transferable MGEs, the traI-PCR developed in this study may serve as a specific screening for the MOB H12 MGE-harboring bacterial strains that are frequently found in the marine environment (Welch et al., 2007;McIntosh et al., 2008;Fricke et al., 2009;Nonaka et al., 2012). The filter mating requires an appropriate recipient strain with known markers to distinguish from donor and the experiment is time consuming. On the other hand, the primers for traI-PCR could rapidly identify the sequence of pAQU1-like plasmids or Inc A/C plasmids or SXT/R391 ICEs.
With the traI-PCR method, we examined 83 strains of tet(M)positive marine bacteria and found 4 Vibrio spp. strains (04Ya001, 04Ya016, 04Ya090, and 04Ya265) that harbored the pAQU1-like plasmids. PFGE and following Southern hybridization showed that the plasmids were 150-350 kb long and carried tet(M). It was also shown that the plasmids could conjugatively be transferred to E. coli at the frequencies of 10 −5 to 10 −4 , generating TC-resistant transconjugants. From three (TJ001W1, TJ090W1, and TJ265W1) of the four transconjugants, the pAQU1-like plasmids were further transferred to another strain of E. coli at comparable frequencies, suggesting that the plasmids were selftransferable.
The entire nucleotide sequence of a tet(M)-bearing pAQU1like plasmid, pAQU2, originally possessed by 04Ya090 isolated from marine sediment was determined and compared with that of pAQU1. It was shown that pAQU2 shared conserved regions, which contained the genes required for plasmid replication, partition and transfer, with pAQU1. These regions, which totally sum up to ∼100 kbp, were defined as the pAQU backbone structure. The other pAQU1-like plasmids harbored in 04Ya001, 04Ya016, and 04Ya265 also shared the pAQU backbone, suggesting that these plasmids together with pAQU1 and pAQU2 could be categorized as the "pAQU group." In addition to tet(M), multiple ARGs [tet(B), sul2, floR and genes similar to bla CARB−9 , mph(A) and mef (A)] were accumulated upstream of the Tra2 region in pAQU1 where multiple transposase genes were localized, as previously described (Nonaka et al., 2012). Multiple, but fewer, ARGs were also found in pAQU2 in the corresponding region. Therefore, it was suggested that the ARGs may have been introduced through transposon-mediated DNA insertion events into the backbone structure of a prototypic plasmid. Insertion of ARGs upstream of the Tra2 region    has also been reported for IncA/C plasmids (Welch et al., 2007;Fricke et al., 2009;Fernandez-Alarcon et al., 2011). Also, arrangement of tra genes in the pAQU backbone was quite similar to that in IncA/C plasmids and SXT/R391 ICEs, suggesting that these MGEs shared synteny. These observations are compatible with the previous study showing a close evolutionary relationship among pAQU1, IncA/C and SXT/R391 ICEs (Nonaka et al., 2012). It is possible that these MGEs had a common ancestor and evolved through acquisition of different replication systems. Among 83 bacterial strains examined in this study, six strains were able to conjugatively confer TC resistance upon E. coli, of which four strains were shown to harbor the pAQU group plasmid. Therefore, it is possible that the pAQU group plasmid played a significant role in dissemination of tet(M) in the aquaculture site. In addition to tet(M), most of the ARGs carried on pAQU1 were also found in the four strains, and all of them were transmitted to E. coli in association with transfer of the plasmid. Therefore, it is suggested that all of the pAQU group plasmids found in this study carried not only tet(M) but also other ARGs. The bacterial strains found to harbor the pAQU group plasmids were isolated from the same aquaculture site, but in different times, suggesting that the pAQU group plasmids were distributed and maintained in the aquaculture environment contributing to dissemination of tet(M) and other ARGs among marine bacteria.
Of six bacterial strains that were able to confer TC resistance on E. coli, two strains did not harbor pAQU group plasmid. One of them, marine sediment strain 04Ya007 which belonged to Shewanella sp. and was able to transfer TC resistance at a high frequency, was shown by nucleotide sequencing of the PCR-amplified traI gene and Southern hybridization to have chromosomal ICEVchMex-related MGE instead. In the transconjugant obtained by filter mating by using 04Ya007 as a donor, chromosomal traI, as well as three other essential genes of ICEVchMex, was detected, indicating the transfer of chromosomal ICEVchMex-related MGE of 04Ya007 to the chromosome of the recipient E. coli. Strain 04Ya007 also had partial sequence characteristic to the pAQU backbone and plasmid-borne ARGs including tet(M). However, none of these pAQU group-related sequences were detected in the transconjugant. While the original ICEVchMex was reported to have no ARGs (Burrus et al., 2006), other SXT/R391 ICEs in the same group were shown to have tetA and tetR (Wozniak et al., 2009). Therefore, it is possible that the TC resistance of the transconjugant was due to an ARG(s) other than tet(M) and tet(B) that was transferred in an ICEVchMex-related MGE-dependent manner, while the tet(M)bearing plasmid with the partial pAQU backbone structure had suffered from the mutation which caused loss of the transferability. Further characterization of ICEVchMex-like MGE and the plasmid in this strain is necessary to test this possibility.
The remaining strain, 04Ya108 isolated from marine sediment, which conferred TC resistance on E. coli at a very low frequency, appeared to have tet(M)-bearing plasmid of non-pAQU group, and the transconjugant was shown to have chromosomeencoded tet(M). Although chromosomal integration of conjugatively transferred plasmid has been reported (Jaoua et al., 1990;Ka and Tiedje, 1994;Ghosh et al., 2005), additional studies are required to elucidate the mechanism by which plasmidborn tet(M) of 04Ya108 integrated in the chromosome of its transconjugant.
In conclusion, PCR targeting traI identified four marine bacterial strains that harbor the pAQU group plasmids closely related to previously described pAQU1 (Nonaka et al., 2012). The four strains were shown by filter mating to be able to serve as donors for TC resistance. The results raise a possibility that the pAQU group plasmids might mediate conjugative transfer of tet(M) and other ARGs in the aquaculture site and play a major role in dissemination of drug resistance among marine bacteria.