- Research article
- Open Access
Plasmid composition in Aeromonas salmonicida subsp. salmonicida 01-B526 unravels unsuspected type three secretion system loss patterns
BMC Genomicsvolume 18, Article number: 528 (2017)
Aeromonas salmonicida subsp. salmonicida is a ubiquitous psychrophilic waterborne bacterium and a fish pathogen. The numerous mobile elements, especially insertion sequences (IS), in its genome promote rearrangements that impact its phenotype. One of the main virulence factors of this bacterium, its type three secretion system (TTSS), is affected by these rearrangements. In Aeromonas salmonicida subsp. salmonicida most of the TTSS genes are encoded in a single locus on a large plasmid called pAsa5, and may be lost when the bacterium is cultivated at a higher temperature (25 °C), producing non-virulent mutants. In a previous study, pAsa5-rearranged strains that lacked the TTSS locus on pAsa5 were produced using parental strains, including 01-B526. Some of the generated deletions were explained by homologous recombination between ISs found on pAsa5, whereas the others remained unresolved. To investigate those rearrangements, short- and long-read high-throughput sequencing technologies were used on the A. salmonicida subsp. salmonicida 01-B526 whole genome.
Whole genome sequencing of the 01-B526 strain revealed that its pAsa5 has an additional IS copy, an ISAS5, compared to the reference strain (A449) sequence, which allowed for a previously unknown rearrangement to occur. It also appeared that 01-B526 bears a second large plasmid, named pAsa9, which shares 40 kbp of highly similar sequences with pAsa5. Following these discoveries, previously unexplained deletions were elucidated by genotyping. Furthermore, in one of the derived strains a fusion of pAsa5 and pAsa9, involving the newly discovered ISAS5 copy, was observed.
The loss of TTSS and hence virulence is explained by one consistent mechanism: IS-driven homologous recombination. The similarities between pAsa9 and pAsa5 also provide another example of genetic diversity driven by ISs.
The bacterium Aeromonas salmonicida subsp. salmonicida is the causative agent of furunculosis, a disease that affects salmonids worldwide. It has a significant economic impact on the fish farming industry . Vaccination and antibiotherapy are the available treatments for furunculosis. However, resistant A. salmonicida subsp. salmonicida strains have emerged [2,3,4,5,6,7] and vaccination is not always effective . Alternative treatment options against this pathogen would thus be beneficial [8,9,10] but developing new treatments requires a better understanding of its underlying mechanisms, such as pathogenicity [4, 11].
This ubiquitous waterborne bacterium shows evidence of lateral gene transfer. A. salmonicida subsp. salmonicida strain genomes have been shown to bear a wide array of mobile elements [4, 12], genomic islands [13,14,15], and transferable plasmids [2, 3, 5, 16]. The close relationship between the mobile elements from this species and those of Salmonella enterica has raised concerns about its ability to act as a reservoir of antibiotic resistance genes, which could transfer to human pathogens [5, 7]. Aeromonas salmonicida subsp. salmonicida is also known for its genomic variability. Different kinds of recombination have formed plasmid and genomic island variants with altered gene contents that result in phenotypic changes [5, 6, 15,16,17,18,19]. This bacterium, beyond mediating gene transfer, could thus provide a new combination of those genes.
In many cases, insertion sequences (ISs) and transposons are responsible for these variations. They were found to participate in plasmid rearrangements, virulence factor inactivation, and the exchange of antibiotic resistance genes [16,17,18,19,20]. It was even suggested that the abundance of ISs (88 copies divided into 10 different types) in the reference genome A449 could maintain the bacterium in a psychrophilic lifestyle to preserve its genomic integrity . Unfortunately, transposable elements also hinder genome sequencing and assembly, which are used to study structural variations [22,23,24,25].
One structural variation mediated by ISs that has been studied in A. salmonicida subsp. salmonicida is the loss of its type three secretion system (TTSS). This needle-like apparatus injects proteins named effectors into the host cells . An interesting point for furunculosis management is the fact that the TTSS is essential for the pathogenicity of this bacterium [11, 27]. In A. salmonicida subsp. salmonicida, the TTSS structural components and many effectors are encoded on a plasmid called pAsa5 (also known as pASvirA) [12, 28], which has been found in many strains [12, 27, 28]. This 155-kilobase-pair (kbp) replicon harbours all TTSS structural genes in a single locus, four TTSS effector genes (aopH, ati2, aopN, and aopO), conjugative transfer genes, many uncharacterized open reading frames, and ISs [11, 12]. Interestingly, pAsa5 is a thermolabile plasmid that is known to lose segments, including the TTSS locus, when the cells are exposed to 25 °C and above [28, 29]. The resulting mutants are non-virulent in fish and in model hosts, making the TTSS loss mechanism an interesting target for avirulence-producing treatment development [27, 29,30,31,32,33,34].
In a previous study, a collection of TTSS-negative mutants were produced by prolonged cultivation at 25 °C of three A. salmonicida subsp. salmonicida parental strains: A449, 01-B526, and 01-B516 . These plasmid-rearranged strains were shown to have two different types of truncated pAsa5.
Plasmids displaying a type 1 deletion profile had lost their TTSS locus, which includes the essential structural proteins, effectors ati2 and aopN, and their respective chaperones (Fig. 1a). Plasmids with a type 2 deletion profile lost their TTSS locus as well as 40 kbp upstream of it. Both deletion types led to complete virulence loss against an alternative host .
Further work revealed that all type 2 and some type 1 deletions could be explained by the same mechanism based on IS-mediated homologous recombination driven by the copies of ISAS11 found on pAsa5 . However, many strains that show a type 1 deletion profile could not be explained by the known IS recombination . It was thus hypothesized that their plasmids had undergone incomplete or more complex recombinations.
As mentioned before, a better understanding of the genomic content and rearrangements in A. salmonicida subsp. salmonicida would be beneficial for virulence-targeted treatment development and for assessing the contributions of this species to gene flow by horizontal transfers. Different TTSS-loss mechanisms seem to occur in strains 01-B526 and A449. The genome of A449 has been already sequenced and fully assembled . In this study, we combined Illumina and single molecule real-time (SMRT) sequencing for the 01-B526 genome to elucidate its plasmid composition. This helped us to identify another TTSS-loss scenario involved in unexplained pAsa5-rearranged strains.
Characterization of pAsa9
A new plasmid sequence was obtained from the 01-B526 strain by combining short and long reads in a hybrid de novo assembly. The new plasmid, named pAsa9, is 76,724 base pairs with a G + C content of 52.76%, and 90 annotated open reading frames (ORF). It can be divided into sections according to its gene content and nucleotide similarities (Fig. 2). The first segment has replication-associated repBA genes (Fig. 2, 0 to 5 kbp) and has high nucleotide identity with two pAsa5 regions at 98.13% and 93.12%, respectively (Fig. 1a). The second region has many conjugation genes (tra, blue rectangles) and again was over 83% similar to pAsa5 (Fig. 1a). The third region, directly following tra genes, had many ORFs coding for hypothetical proteins (cyan, Fig. 2) or with less described functions (black, Fig. 2), but was still 84% similar (Fig. 1a) to pAsa5. Finally, pAsa9, the last region, was mostly filled with fragmented and complete transposase genes (Fig. 2, red rectangles).
Since ISs are a prominent component of A. salmonicida subsp. salmonicida genome , we investigated those of pAsa9. ISs and IS-related sequences were found solely in the last quarter of the sequence (Fig. 2, 57 kbp to end). More partial than complete ISs could be identified, suggesting ancestral transposition events . Of the three complete ISs found on pAsa9, two, ISAS10 and ISAS5, have already been identified and described in A. salmonicida subsp. salmonicida notably on the reference pAsa5 . ISAS35, an IS that had not been described previously, was identified based on its transposase gene and inverted repeats homology. It belongs to the ISBst12 group, of which members have been found in Klebsiella pneumoniae (ISKpn15, best score 248) and A. salmonicida subsp. salmonicida (ISAs21, best score 167) based on the ISFinder database . A blast analysis against deposited Aeromonas genomes (taxid:642, October 2016) showed ISAS35 is only found in three A. salmonicida subsp. salmonicida genomes: 01-B526 (accession number AGVO01000000), J231 (LSGY01000000), and J227 (LSGX01000000).
Comparison between pAsa9 and pAsa5
pAsa9 was sequenced from the strain 01-B526, which is already known to contain other plasmids: pAsa1, pAsa2, pAsa3, pAsal1 and pAsa5 . Interestingly, pAsa9 shared long (up to 24 kbp), similar regions with pAsa5 (Fig. 1a). Shared similarities made it impossible to detect, assemble, and characterize pAsa9 in the previous studies focusing on 01-B526 genomics that were based on short-read sequencing only [22, 23]. High-throughput sequencing was sufficient to cover its specific regions since some contigs of A. salmonicida subsp. salmonicida J231 and J227 strains found in GenBank matched against pAsa9. However, the repeated regions between pAsa5 and pAsa9 did not allow for a complete assembly with previous 01-B526 datasets [22, 23]. SMRT and Illumina sequencing resulted in the assembly of both pAsa5 and pAsa9 as single contigs, which allowed us to distinguish them for the first time.
Apart from ISs, both plasmids share regions that encode conjugative transfer and plasmid replication genes. pAsa5 and pAsa9 have similar and contiguous set of tra genes (between 83% and 91%) and similar repBA genes (between 93% and 98%). Some regions, such as the rep surroundings and the traI gene, are repeated thrice between the plasmids: once on pAsa9 and twice on pAsa5. The distinct regions on pAsa9 compared to pAsa5 contained mostly hypothetical protein-coding sequences. On the other hand, pAsa5 unique regions carried hypothetical protein-coding sequences, but also ISs and, more importantly, TTSS-related genes. pAsa9 ISs and IS-associated genes are concentrated in the same region, while pAsa5 ISs are scattered throughout its length, and some are proximal to TTSS genes. For instance, the virulence-essential TTSS locus is framed by two ISAS11, as described earlier . This region also disrupts the G + C skew in pAsa5 compared to pAsa9, suggesting a recent acquisition (Additional file 1). A 10 kbp region comprising another TTSS effector gene, aopO, is also framed by ISs, two ISAS10s. The third TTSS effector gene, aopH, is framed by non-identical ISs (ISAS11 and ISAS4).
pAsa5 also shares similarities with pAsa6 (18.5 kbp), another A. salmonicida subsp. salmonicida plasmid that also has the aopH gene . Sequence alignment showed that the genetic content and architecture of pAsa6 are closer to pAsa5 than pAsa9, since it carries the aopH section (Additional file 2). Furthermore, its similarities are concentrated around pAsa5’s first rep region and insertion sequences. The significant similarities between pAsa5 and pAsa9 suggest that pAsa6 has probably derived from pAsa5 through successive deletions instead of pAsa5 evolved from the fusion of pAsa6 with other plasmid(s).
Genotyping unexplained type 1 loss profile and pAsa9 occurrence
01-B526 pAsa5 was aligned against the reference plasmid from the strain A449. The two plasmids had a high nucleotide identity (more than 99.5%) across their length. The only structural variation is an additional ISAS5 inserted between the first ISAS4 and the second rep region on 01-B526 pAsa5 (Fig. 1a). Since this IS copy is upstream from ISAS5A from the previous annotation , it was designated as ISAS5Z for the purpose of this paper. ISAS5Z (upstream from tra genes) and ISAS5A (downstream from the TTSS locus) are in the same orientation, hence they form a valid pair for IS-mediated loss by recombination.
This deletion would encompass both conjugation and TTSS loci, and should result in another loss profile for pAsa5 genotyping, a feature that had not been observed previously (Fig. 1b) [17, 29]. However, due to its similarity with pAsa9, primers hybridizing in the conjugation locus (traE, traC, and traD) can also hybridize on pAsa9, leading to false positives for this region (Fig. 1). ISAS5Z-ISAS5A recombination thus leads to a type 1 deletion profile in pAsa9-bearing strains when the previous genotyping primer set is used (Fig. 1b). The unexplained type 1 deletion profiles found in pAsa5-rearranged strains could thus be explained by ISAS5Z-A recombination, but only if the strains bore pAsa9 and an ISAS5Z.
All parental and derived strains were directly tested with primers flanking ISAS5Z 3′ and ISAS5A 5′ to assess the new ZA recombination. All parental strains and daughter strains for which a recombination pattern had already been found were negative to that new primer pair.
In the unexplained derived strains, all but one (01-B526-R19) were positive for the new ZA amplification, which gave a 2888 bp amplification product (Table 1). Sanger sequencing for that amplicon showed that it encloses an ISAS5, ISAS5Z flanking 3′, and ISAS5A flanking 5′, confirming the recombination. Those strains had thus experienced a homologous recombination between ISAS5Z and ISAS5A, leading to the loss of both TTSS and conjugation loci. All derived strains were then tested for the presence of pAsa9 using a set of genotyping primers that targeted its unique region (Additional file 3). A449-derived strains were negative for pAsa9, which was expected, since A449 does not carry this plasmid. All the other pAsa5-rearranged strains were positive for pAsa9. Since all newly assigned ZA pAsa5-bearing strains also carried pAsa9, it explained why the conjugation locus loss had not been detected in those strains previously .
01-B516 and its assigned daughter strains
The parental strains, A449, 01-B526 and 01-B516 were then tested for the presence of pAsa9. It was found that neither A449 (as expected) nor 01-B516 had the plasmid. However, all three 01-B516-derived strains (01-B516–3, 01-B516–11, and 01-B516–30) had pAsa9, creating a disparity between parental and offspring strains. It was thus hypothesized that during mutant production, some 01-B526-derived strains had been mistakenly labelled as 01-B516.
To verify that 01-B516 did not have any form of pAsa9, plasmid profiles of all parental strains, one BC pAsa5-bearing strain (01-B526-R2), and one possibly mislabelled derived strain (01-B516–11), were visualized in agarose gel electrophoresis (Fig. 3a). As expected, A449 showed one band corresponding to pAsa5 and pAsa4 (155 kbp and 167 kbp, respectively). 01-B526 showed two bands, one for pAsa5 and one for pAsa9. 01-B516 showed a single large band whose weight matched pAsa5 (Fig. 3a). 01-B526-R2 showed the pAsa9 band as well as a BC-rearranged pAsa5 band. However, 01-B516–11, which was found to bear a ZA-rearranged pAsa5 according to our previous genotyping, gave two bands: one for the ZA pAsa5 and one whose size matched pAsa9. This result suggested this strain has a 01-B526-derived background and was thus mislabelled.
01-B516 genomic DNA was then sequenced using short-read, high-throughput sequencing (Illumina). When the resulting contigs were aligned on pAsa5 and pAsa9, it showed 01-B516 did not have sequences unique to pAsa9. Moreover, 01-B516 did not contain ISAS5Z, making an ISAS5Z-ISAS5A recombination impossible in this strain, nor in any daughter strain.
Derived strains 01-B516–3, 01-B516–11 and 01-B516–30 were also genotyped for AsaGEI1a (accession number KJ626178) presence and were found positive for that genomic island, whereas 01-B516 carries AsaGEI2a . AsaGEIs (for Aeromonas salmonicida GEnomic Island) are found in some A. salmonicida subsp. salmonicida genomes. Each carrying strain only has one copy of one type of AsaGEIs. They have variable insertion sites (indicated by the number in their name) and gene content (indicated by their letter). They contain some phage-related genes and many ORFs coding for hypothetical protein. Consequently, their role remains to be determined [13,14,15].
Moreover, AsaGEI types are linked with their geographic origin and strains bearing AsaGEI1a are found in North America . Since all 01-B516 ‘derived’ strains genotyped positive for pAsa9 and AsaGEI1a, and that two of them have a ZA-rearranged pAsa5, we concluded these three strains were mislabelled when they were obtained, and are instead 01-B526-derived strains.
The prevalence of pAsa9 in a collection of A. salmonicida subsp. salmonicida isolates was then determined by genotyping with primers covering its specific region (Additional file 3). Among the 154 genotyped strains, 35 had the plasmid and two generated amplicons using some primers only, suggesting another plasmid (Additional file 4). Furthermore, in our strain set, pAsa9 was always found with another element, the genomic island AsaGEI1a.
Plasmid fusion in 01-B526-R19
One type 1 deletion profile strain, 01-B526-R19, could not be explained by any IS recombination, even with the ZA recombination discovery. 01-B526-R19 genomic DNA was thus sequenced by Illumina technology to detect potential new genomic content. Contig alignment to 01-B526 pAsa5 and pAsa9 showed it lacked the tra and TTSS regions on pAsa5, suggesting an ISAS5Z-A recombination, a fact which had been suspected but not proven by PCR. It also had a complete pAsa9, as shown by contig alignment.
The 01-B526-R19 recombination was deemed to be complex, so its plasmids were visualized on gel (Fig. 3b). 01-B526-R19 displayed a single large plasmid band at approximately 150 kbp, higher than a ZA-rearranged pAsa5 band, and no pAsa9 band was detected, even if its presence had been assessed by PCR and sequencing (Fig. 3b).
However, since the weight of the 01-B526-R19 band matched the sum of ZA-rearranged pAsa5 and pAsa9, and since both plasmids had one ISAS5 copy, we hypothesized that both large plasmids had merged together by ISAS5 recombination (Fig. 4a). Primers flanking ISAS5ZA in pAsa5 and pAsa9 ISAS5 were used in a long amplification PCR (Additional file 4). Both sides of the merge were successfully amplified and confirmed by Sanger sequencing, indicating 01-B526-R19’s rearranged pAsa5 and its pAsa9 had indeed recombined into a single plasmid, and shedding light on the last unexplained strain recombination (Fig. 4b).
In this study, we used long-read sequencing (PacBio SMRT technology) to uncover a new A. salmonicida subsp. salmonicida plasmid, pAsa9. It would have been impossible to assemble pAsa9 with short-read sequencing information alone, due to its similarity with another well-studied, virulence-essential plasmid called pAsa5 [12, 29, 30]. Even with mate-pair sequencing, which allowed us, with a certain configuration, to assemble the A. salmonicida subsp. salmonicida chromosome in one scaffold , the pAsa9 assembly was complicated by its 24 kbp-long repeated sequence with pAsa5.
SMRT sequencing can be used to reduce contig numbers and manual gap filling for most microbial genomes, especially those that have many long repeats; this feature applies to A. salmonicida subsp. salmonicida . Single-read sequencing was previously used to identify and assemble plasmids that shared a common backbone in the same strain . PacBio or other long-read sequencing technologies should be used in future studies to facilitate the discovery of new plasmids in A. salmonicida subsp. salmonicida, especially if they are presumed to share repeats with the rest of the genome. Also, some species in the Aeromonas genus, such as Aeromonas media , have complex genomes and would thus benefit from that technology .
The newly discovered plasmid, pAsa9, was found in A. salmonicida subsp. salmonicida strains that also bear the genomic island AsaGEI1a. Thus, its presence is limited to North America so far based on the strains analysed. Interestingly, pAsa9 is found concurrently with pAsa5, even though their replication genes (repBA) are similar. Although the discovered plasmid, pAsa9, shares highly similar sequences with pAsa5, its unique content and architecture show that it is not a mere copy.
In fact, we propose that pAsa9 provides information on pAsa5 evolution and gene acquisition. The two plasmids are similar and have syntenic sequences, which suggests they share a common ancestor whose backbone should be closer to pAsa9 than pAsa5. Apart from its last quarter (Fig. 2), which encloses all mobile element-related sequences, pAsa9 is free from IS or sequence duplication. In contrast, pAsa5 has sequence duplications in core genes (replication and transfer) , and many ISs. Furthermore, some IS localizations follow TTSS genes.
We thus suggest that ISs were key elements in TTSS-related gene acquisition in pAsa5, which arose on a plasmid backbone close to pAsa9. The discovery of pAsa9 also places pAsa6, another pAsa5-related plasmid, as a pAsa5-like descendant rather than an ancestor, since its structure is closer to the latter (Additional file 2). Together, those results further prove that insertion sequences have a role in A. salmonicida subsp. salmonicida genetic diversity, as proposed earlier .
PacBio sequencing also provided new architectural details that allowed us to revisit previous experimental results on pAsa5 recombination and the associated loss of TTSS and virulence. With this information, all pAsa5 deletions that were generated in the study of Daher et al. can be regrouped under a consistent and simple mechanism, which is IS-mediated homologous recombination . These results also allowed us to explain previous observations of pAsa5-rearranged strains that had lost their TTSS locus .
For instance, while attempting to produce TTSS-negative mutants, strain 01-B526 was more prone to producing plasmid-rearranged strains than strain A449 . This can be explained by the fact that in 01-B526 pAsa5, another IS (ISAS5Z) can be a target for recombination.
With 2.6 kbp, ISAS5 offers a large span of homologous sequences for recombination to occur, which could explain why ISAS5 rearrangements were observed more frequently in 01-B526-derived strains than ISAS11 rearrangements. ISAS5 was already known as an active insertion sequence in A. salmonicida subsp. salmonicida. It was found to have transposed at multiple locations on pAsal1 [18, 19]. The results thus suggest that ISAS5 is both an insertion sequence capable of transposition, and a target for host-mediated recombination. 01-B526 also gave only one ISAS11A-C recombinant with an atypical genotyping profile, 01-B526-R4, which can now be explained by the presence of pAsa9 and the false positive tra targets it creates (Fig. 1). Thus, PacBio sequencing both provided an explanation for shaded areas left by the previous study, and simplified some cases.
One exception is 01-B526-R19, which has a fusion between its ZA-rearranged pAsa5 and its pAsa9. This merged plasmid is a rarity since it resulted from two sequential recombinations, one between ISAS5Z and 5A, and one with ISAS5ZA and pAsa9 ISAS5. Such double recombinations between ISAS11s were tested in 2012 with no success . This may be due to the fact, as mentioned earlier, ISAS11 is half the size of ISAS5.
The reassignment of derived strains 01-B516–3, 01-B516–11 and 01-B516–30 to 01-B526 leaves 01-B516 without pAsa5-rearranged descendants that would prove the strain is capable of IS-targeted recombination. Difficulties in producing 01-B516 pAsa5-rearranged strains were observed in the original experiment . This problem could be caused either by the structure of 01-B516 pAsa5 or by the lack of suitable recombinase. The Illumina sequencing showed that ISAS11B and ISAS11C are in the correct position relative to each other to mediate the deletion. BC recombination and TTSS loss should thus have been observed. Consequently, the lack of rearrangement in 01-B516 may be due to a missing or inactive recombinase.
Furunculosis management needs new alternative treatments to antibiotherapy and vaccination. Avirulent treatment suppresses essential pathogenicity mechanisms without killing the bacteria, and could be an option against A. salmonicida subsp. salmonicida. However, full understanding of the TTSS loss and the underlying IS-mediated recombination is a prerequisite to finding new ways to attenuate furunculosis-causing strains. Now, factors promoting these events can be investigated in the hopes of finding a specific, aquaculture-compatible condition that would trigger recombination of pAsa5.
Our results showed that PacBio sequencing could provide crucial information for plasmid architecture studies by allowing for proper separation of long repeated regions. In this case, it led to the discovery of a new plasmid, pAsa9, which in turn allowed refinement and completion of the recombination model described for pAsa5 in 2012 .
Regrouping all TTSS-loss scenarios under a simple rule was an important step in investigating the mechanism for avirulence-producing treatments. The relationship between the pAsa9 and pAsa5 architectures demonstrates the importance of ISs in producing biologically-relevant genomic diversity in A. salmonicida subsp. salmonicida.
Bacterial isolates and growth conditions
The A. salmonicida subsp. salmonicida strains that were used in this study and are listed in Table 1 (strains listed in Tanaka et al.) and Additional file 4 (all other strains, used for genotyping). In Table 1, pAsa5-rearranged strains (that display type 1 or type 2 loss profiles) were obtained by culturing virulent parental strains (A449, 01-B516 and 01-B526) for 2 weeks at 25 °C . All A. salmonicida strains were grown on furunculosis agar (10 g of tryptone, 5 g of yeast extract, 1 g of L-tyrosine, 2.5 g of NaCl, and 15 g of agar per litre of distilled water) for 2 or 3 days at 18 °C . Paracoccus aminophilus JCM 7686 was used as a DNA size marker in plasmid visualization . It was grown on tryptic soy agar (EMD Millipore, Ontario, Canada) for 3 days at 30 °C.
DNA extraction and sequencing
Total genomic DNA was extracted for strains 01-B526 and 01-B516, and derived strain 01-B526-R19 (Table 1) using DNeasy Blood and Tissue kits (Qiagen, Canada). Libraries were prepared with a KAPA Hyper Prep Kit and sequencing was done at the Plateforme d’Analyse Génomique of the Institut de Biologie Intégrative et des Systèmes (IBIS, Université Laval, Quebec City, Canada) using the Illumina MiSeq system. Strain 01-B526 total genomic DNA was also extracted by phenol/chloroform by following the protocol Extracting DNA Using Phenol-Chloroform provided by Pacific Biosciences (http://www.pacb.com) and was sequenced at the Génome Québec Innovation Centre (McGill University, Montreal, Canada) using a Pacific BioScience RS II system. A hybrid de novo assembly was used for 01-B526 reads using Spades version 3.6.0 (kmer lengths of 21, 33, 55, 77, 99, 127) . 01-B516 and 01-B526-R19 reads were assembled using A5-miseq version 20150522 .
The 01-B526 assembly allowed us to obtain the large plasmids as complete closed sequences. pAsa5 and pAsa9 were annotated using the webserver RAST (February 2016) . Manual verification was done using Artemis version 16.0.0 using pAsa5 reference sequence (NC_009350) for comparison and blastp (non-redundant protein sequences database, May 2016) [47, 48]. Manual verification of the mobile genetic elements was done with ISFinder assistance (June 2016) . For 01-B516 and 01-B526-R19, contigs were aligned on 01-B526 pAsa5 and pAsa9 sequences using CONTIGuator version 2.7.4 . Sequence alignment visualization was done with ACT version 13.0.0, DNA Plotter version 10.2 and EasyFig version 2.1 [50,51,52].
For PCR, amplifications up to 3 kbp were performed as previously described , with the exception of the elongation time, which was 1 min per 1 kbp of amplicon. For amplicons longer than 3 kbp, each reaction was built as follows: 3.75 μL of 2 mM dNTPs, 1 μL of each 100 ng/μL forward and reverse primers, 1 μL of 100 ng/μL template, 12.25 μL of H2O, 5 μL of 5X LongAmp Taq buffer (New England BioLabs, Massachusetts, US), and 1 μL LongAmp Taq (New England BioLabs). For those amplicons, the elongation time was 50 s per 1 kb of amplicon. New PCR primers were designed manually and verified with Oligoanalyser 3.1 from Integrated DNA Technologies (IDT, http://www.idtdna.com/calc/analyzer). The PCR assays were performed at least twice, and appropriate positive and negative controls were included with each assay. The PCR primers are listed in Additional file 3.
Plasmid visualization by agarose gel electrophoresis
A. salmonicida subsp. salmonicida plasmids were visualized by agarose gel electrophoresis using a protocol adapted from Wheatcroft et al. . 160 μL of overnight P. aminophilus JCM 7686 in tryptic soy broth (EMD Millipore) or A. salmonicida subsp. salmonicida in lysogeny broth (EMD Millipore) were mixed with 1 mL ice-cold 2% N-lauroyl sarcosine and centrifuged at 17200 x g for 10 min at 4 °C. Supernatants were discarded and pellets were suspended in 40 μL 10 mM Tris, 10 mM EDTA, 4 mg/mL RNAse, 1 mg/mL lysozyme, 1 mg/mL bromophenol blue and 15.2% Ficoll solution. Samples were incubated on ice for 45 min. A 0.75% agarose gel in Tris-borate-EDTA buffer (TBE 1X) was loaded with 40 μL of SDS 10%, 0.5% xylene cyanol solution, which was migrated backwards at 100 V for 15 min. DNA samples were then centrifuged at 376 x g for 20 s, then 20 μL were loaded on the gel. The polarity was inverted back and samples were migrated 30 min at 40 V, and then 5–7 h at 100 V. Gels were stained in 1 μg/mL ethidium bromide for 40 min, rinsed in water for 1 h to remove the excess of ethidium bromide, and photographed using short-length UV transilluminator.
Open reading frame
Single molecule real-time
Type three secretion system
Derome N, Boutin S, Llewellyn M, Gauthier J. Opportunistic pathogens of fish. In: Hurst CJ, editor. The Rasputin effect: when commensals and symbionts become parasitic. New York: Springer; 2016. p. 649. doi:10.1007/978-3-319-28170-4.
Sørum H, L'Abée-Lund T, Solberg A, Wold A. Integron-containing IncU R plasmids pRAS1 and pAr-32 from the fish pathogen Aeromonas salmonicida. Antimicrob Agents Chemother. 2003;47(4):1285–90.
McIntosh D, Cunningham M, Ji B, Fekete FA, Parry EM, Clark SE, et al. Transferable, multiple antibiotic and mercury resistance in Atlantic Canadian isolates of Aeromonas salmonicida subsp. salmonicida is associated with carriage of an IncA/C plasmid similar to the Salmonella enterica plasmid pSN254. J Antimicrob Chemother. 2008;61(6):1221–8.
Piotrowska M, Popowska M. Insight into the mobilome of Aeromonas strains. Front Microbiol. 2015;6:494.
Vincent AT, Trudel MV, Paquet VE, Boyle B, Tanaka KH, Dallaire-Dufresne S, et al. Detection of variants of the pRAS3, pAB5S9, and pSN254 plasmids in Aeromonas salmonicida subsp. salmonicida: multidrug resistance, interspecies exchanges, and plasmid reshaping. Antimicrob Agents Chemother. 2014;58(12):7367–74.
Vincent AT, Emond-Rheault JG, Barbeau X, Attéré SA, Frenette M, Lagüe P, et al. Antibiotic resistance due to an unusual ColE1-type replicon plasmid in Aeromonas salmonicida. Microbiology. 2016;162(6):942–53.
Trudel MV, Vincent AT, Attéré SA, Labbé M, Derome N, Culley AI, et al. Diversity of antibiotic-resistance genes in Canadian isolates of Aeromonas salmonicida subsp. salmonicida: dominance of pSN254b and discovery of pAsa8. Sci Rep. 2016;6:35617.
Dallaire-Dufresne S, Tanaka KH, Trudel MV, Lafaille A, Charette SJ. Virulence, genomic features, and plasticity of Aeromonas salmonicida subsp. salmonicida, the causative agent of fish furunculosis. Vet Microbiol. 2014;169(1–2):1–7.
Imbeault S, Parent S, Lagacé M, Uhland CF, Blais J-F. Using Bacteriophages to prevent furunculosis caused by Aeromonas salmonicida in farmed brook trout. J Aquat Anim Health. 2006;18(3):203–14.
Tanaka KH, Frenette M, Charette SJ. IS-mediated loss of virulence by Aeromonas salmonicida: a tangible piece of an evolutionary puzzle. Mobile genetic elements. 2013;3(1):e23498.
Frey J, Origgi FC. Type III secretion system of Aeromonas salmonicida undermining the Host's immune response. Front Mar Sci. 2016;3:130.
Reith ME, Singh RK, Curtis B, Boyd JM, Bouevitch A, Kimball J, et al. The genome of Aeromonas salmonicida subsp. salmonicida A449: insights into the evolution of a fish pathogen. BMC Genomics. 2008;9:427.
Emond-Rheault JG, Vincent AT, Trudel MV, Brochu F, Boyle B, Tanaka KH, et al. Variants of a genomic island in Aeromonas salmonicida subsp. salmonicida link isolates with their geographical origins. Vet Microbiol. 2015;175(1):68–76.
Emond-Rheault JG, Vincent AT, Trudel MV, Frey J, Frenette M, Charette SJ. AsaGEI2b: a new variant of a genomic island identified in the Aeromonas salmonicida subsp. salmonicida JF3224 strain isolated from a wild fish in Switzerland. FEMS Microbiol Lett. 2015;362(13):fnv093.
Long M, Nielsen TK, Leisner JJ, Hansen LH, Shen ZX, Zhang QQ, et al. Aeromonas salmonicida subsp. salmonicida Strains isolated from Chinese freshwater fish contain a novel genomic island and possible regional-specific mobile genetic elements profiles. FEMS Microbiol Lett. 2016;363(17):fnw190.
Tanaka KH, Vincent AT, Trudel MV, Paquet VE, Frenette M, Charette SJ. The mosaic architecture of Aeromonas salmonicida subsp. salmonicida pAsa4 plasmid and its consequences on antibiotic resistance. PeerJ. 2016;4:e2595.
Tanaka KH, Dallaire-Dufresne S, Daher RK, Frenette M, Charette SJ. An insertion sequence-dependent plasmid rearrangement in Aeromonas salmonicida causes the loss of the type three secretion system. PLoS One. 2012;7(3):e33725.
Trudel MV, Tanaka KH, Filion G, Daher RK, Frenette M, Charette SJ. Insertion sequence AS5 (ISAS5) is involved in the genomic plasticity of Aeromonas salmonicida. Mobile genetic elements. 2013;3(3):e25640.
Attéré SA, Vincent AT, Trudel MV, Chanut R, Charette SJ. Diversity and homogeneity among small plasmids of Aeromonas salmonicida subsp. salmonicida linked with geographical origin. Front Microbiol. 2015;6:1274.
Gustafson CE, Chu SJ, Trust TJ. Mutagenesis of the paracrystalline surface protein array of Aeromonas salmonicida by endogenous insertion elements. J Mol Biol. 1994;237(4):452–63.
Vincent AT, Trudel MV, Freschi L, Nagar V, Gagné-Thivierge C, Levesque RC, et al. Increasing genomic diversity and evidence of constrained lifestyle evolution due to insertion sequences in Aeromonas salmonicida. BMC Genomics. 2016;17(1):44.
Charette SJ, Brochu F, Boyle B, Filion G, Tanaka KH, Derome N. Draft genome sequence of the virulent strain 01-B526 of the fish pathogen Aeromonas salmonicida. J Bacteriol. 2012;194(3):722–3.
Vincent AT, Boyle B, Derome N, Charette SJ. Improvement in the DNA sequencing of genomes bearing long repeated elements. J Microbiol Methods. 2014;107:186–8.
Beatson SA, Walker MJ. Tracking antibiotic resistance. Science. 2014;345(6203):1454–5.
Vincent AT, Tanaka KH, Trudel MV, Frenette M, Derome N, Charette SJ. Draft genome sequences of two Aeromonas salmonicida subsp. salmonicida isolates harboring plasmids conferring antibiotic resistance. FEMS Microbiol Lett. 2015;362(4):1–4.
Buttner D. Protein export according to schedule: architecture, assembly, and regulation of type III secretion systems from plant- and animal-pathogenic bacteria. Microbiol Mol Biol Rev. 2012;76(2):262–310.
Burr SE, Pugovkin D, Wahli T, Segner H, Frey J. Attenuated virulence of an Aeromonas salmonicida subsp. salmonicida type III secretion mutant in a rainbow trout model. Microbiology. 2005;151(Pt 6):2111–8.
Stuber K, Burr SE, Braun M, Wahli T, Frey J. Type III secretion genes in Aeromonas salmonicida subsp. salmonicida are located on a large thermolabile virulence plasmid. J Clin Microbiol. 2003;41(8):3854–6.
Daher RK, Filion G, Tan SG, Dallaire-Dufresne S, Paquet VE, Charette SJ. Alteration of virulence factors and rearrangement of pAsa5 plasmid caused by the growth of Aeromonas salmonicida in stressful conditions. Vet Microbiol. 2011;152(3–4):353–60.
Burr SE, Stuber K, Wahli T, Frey J. Evidence for a type III secretion system in Aeromonas salmonicida subsp. salmonicida. J Bacteriol. 2002;184(21):5966–70.
Burr SE, Wahli T, Segner H, Pugovkin D, Frey J. Association of Type III secretion genes with virulence of Aeromonas salmonicida subsp. salmonicida. Dis Aquat Org. 2003;57(1–2):167–71.
Dacanay A, Knickle L, Solanky KS, Boyd JM, Walter JA, Brown LL, et al. Contribution of the type III secretion system (TTSS) to virulence of Aeromonas salmonicida subsp. salmonicida. Microbiology. 2006;152(Pt 6):1847–56.
Dacanay A, Boyd JM, Fast MD, Knickle LC, Reith ME. Aeromonas salmonicida Type I pilus system contributes to host colonization but not invasion. Dis Aquat Org. 2010;88(3):199–206.
Froquet R, Cherix N, Burr SE, Frey J, Vilches S, Tomas JM, et al. Alternative host model to evaluate Aeromonas virulence. Appl Environ Microbiol. 2007;73(17):5657–9.
Siguier P, Gourbeyre E, Chandler M. Bacterial insertion sequences: their genomic impact and diversity. FEMS Microbiol Rev. 2014;38(5):865–91.
Siguier P, Perochon J, Lestrade L, Mahillon J, Chandler M. ISfinder: the reference centre for bacterial insertion sequences. Nucleic Acids Res. 2006;34(Database issue):D32–6.
Najimi M, Balado M, Lemos ML, Osorio CR. Genetic characterization of pAsa6, a new plasmid from Aeromonas salmonicida subsp. salmonicida that encodes a type III effector protein AopH homolog. Plasmid. 2009;61(3):176–81.
Koren S, Harhay GP, Smith TP, Bono JL, Harhay DM, McVey SD, et al. Reducing assembly complexity of microbial genomes with single-molecule sequencing. Genome Biol. 2013;14(9):R101.
Mehdizadeh Gohari I, Kropinski AM, Weese SJ, Parreira VR, Whitehead AE, Boerlin P, et al. Plasmid characterization and chromosome analysis of two netF+ Clostridium perfringens isolates associated with foal and canine necrotizing enteritis. PLoS One. 2016;11(2):e0148344.
Chai B, Wang H, Chen X. Draft genome sequence of high-melanin-yielding Aeromonas media strain WS. J Bacteriol. 2012;194(23):6693–4.
Garcillan-Barcia MP, Alvarado A, de la Cruz F. Identification of bacterial plasmids based on mobility and plasmid population biology. FEMS Microbiol Rev. 2011;35(5):936–56.
Hanninen ML, Hirvela-Koski V. Molecular and phenotypic methods for the characterization of atypical Aeromonas salmonicida. Vet Microbiol. 1997;56(1–2):147–58.
Dziewit L, Czarnecki J, Wibberg D, Radlinska M, Mrozek P, Szymczak M, et al. Architecture and functions of a multipartite genome of the methylotrophic bacterium Paracoccus aminophilus JCM 7686, containing primary and secondary chromids. BMC Genomics. 2014;15
Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol. 2012;19(5):455–77.
Coil D, Jospin G, Darling AE. A5-miseq: an updated pipeline to assemble microbial genomes from Illumina MiSeq data. Bioinformatics. 2015;31(4):587–9.
Overbeek R, Olson R, Pusch GD, Olsen GJ, Davis JJ, Disz T, et al. The SEED and the rapid annotation of microbial genomes using subsystems technology (RAST). Nucleic Acids Res. 2014;42(Database issue):D206–14.
Carver T, Harris SR, Berriman M, Parkhill J, McQuillan JA. Artemis: an integrated platform for visualization and analysis of high-throughput sequence-based experimental data. Bioinformatics. 2012;28(4):464–9.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215(3):403–10.
Galardini M, Biondi EG, Bazzicalupo M, Mengoni A. CONTIGuator: a bacterial genomes finishing tool for structural insights on draft genomes. Source Code Biol Med. 2011;6:11.
Carver TJ, Rutherford KM, Berriman M, Rajandream MA, Barrell BG, Parkhill J. ACT: the Artemis comparison tool. Bioinformatics. 2005;21(16):3422–3.
Carver T, Thomson N, Bleasby A, Berriman M, Parkhill J. DNAPlotter: circular and linear interactive genome visualization. Bioinformatics. 2009;25(1):119–20.
Sullivan MJ, Petty NK, Beatson SA. Easyfig: a genome comparison visualizer. Bioinformatics. 2011;27(7):1009–10.
Wheatcroft R, McRae DG, Miller RW. Changes in the Rhizobium meliloti genome and the ability to detect Supercoiled plasmids during Bacteroid development. Mol Plant-Microbe Interact. 1990;3(1):9–17.
Rutherford K, Parkhill J, Crook J, Horsnell T, Rice P, Rajandream MA, et al. Artemis: sequence visualization and annotation. Bioinformatics. 2000;16(10):944–5.
Dautremepuits C, Fortier M, Croisetiere S, Belhumeur P, Fournier M. Modulation of juvenile brook trout (Salvelinus fontinalis) cellular immune system after Aeromonas salmonicida challenge. Vet Immunol Immunopathol. 2006;110(1–2):27–36.
We thank Mégane Paccaud, Mathilde Goldschmitt, and Romain Chanut for technical support. We also thank the Laboratoire de bactériologie clinique, Faculté de médecine vétérinaire, Université de Montréal (Montreal, QC, Canada), the Félix d’Hérelle Reference Center (Département de biochimie, de microbiologie et de bio-informatique, Université Laval, Quebec City, QC, Canada), M. Reith (National Research Council laboratories, Halifax, Nova Scotia), J. Frey (University of Bern), and R.M.W. Stevenson (University of Guelph) for the A. salmonicida subsp. salmonicida isolates.
KHT and ATV received scholarships from the Fonds de recherche québécois – Nature et Technologies (FRQNT) and the Natural Sciences and Engineering Research Council of Canada (NSERC). JGER received a scholarship from Ressources Aquatiques Québec. This work was funded by a NSERC grant to SJC (RGPIN-2014-04595) and the Innovamer Program of the Ministère de l’Agriculture, des Pêcheries et de l’Alimentation du Québec. The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Availability of data and materials
Aeromonas salmonicida subsp. salmonicida strain 01-B526 plasmids pAsa5 and pAsa9 complete sequences were deposited in GenBank with the accession numbers KY555069 (http://www.ncbi.nlm.nih.gov/nuccore/KY555069) and KY555070 (http://www.ncbi.nlm.nih.gov/nuccore/KY555070), respectively. The authors declare that all other data supporting the findings of this study are available within the article and its supplementary information files.
Ethics approval and consent to participate
No animal material was used in this study. All fish-related information given in this study was provided by each strain’s donors (see below) as additional data.
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.
Circular map of pAsa5 of 01-B526 strain with G + C skew. pAsa5 sequence and features were visualized with Artemis and DNAPlotter . From the outermost ring moving inwards, the first two circles shows open reading frames (forward, then reverse strand) in the colours described in Fig. 1. The third circle shows mobile elements. The fourth circle shows the G + C skew, using the following colours: purple: below average; yellow: above average. (PDF 171 kb)
Alignment between pAsa5 and pAsa9 from 01-B526 with pAsa6. Methods and colour shown are the same as Fig. 1. (PDF 165 kb)
Primers used in this study. This table lists the primers used for pAsa5, pAsa9 and AsaGEI1a genotyping, IS recombination detection, plasmid fusion detection and positive control. (DOCX 66 kb)
Strains genotyped in this study. These strains were genotyped for the presence of pAsa9 and AsaGEI1a. (DOCX 255 kb)