Evolution and dynamics of megaplasmids with genome sizes larger than 100 kb in the Bacillus cereus group
© Zheng et al.; licensee BioMed Central Ltd. 2013
Received: 21 September 2013
Accepted: 25 November 2013
Published: 2 December 2013
Plasmids play a crucial role in the evolution of bacterial genomes by mediating horizontal gene transfer. However, the origin and evolution of most plasmids remains unclear, especially for megaplasmids. Strains of the Bacillus cereus group contain up to 13 plasmids with genome sizes ranging from 2 kb to 600 kb, and thus can be used to study plasmid dynamics and evolution.
This work studied the origin and evolution of 31 B. cereus group megaplasmids (>100 kb) focusing on the most conserved regions on plasmids, minireplicons. Sixty-five putative minireplicons were identified and classified to six types on the basis of proteins that are essential for replication. Twenty-nine of the 31 megaplasmids contained two or more minireplicons. Phylogenetic analysis of the protein sequences showed that different minireplicons on the same megaplasmid have different evolutionary histories. Therefore, we speculated that these megaplasmids are the results of fusion of smaller plasmids. All plasmids of a bacterial strain must be compatible. In megaplasmids of the B. cereus group, individual minireplicons of different megaplasmids in the same strain belong to different types or subtypes. Thus, the subtypes of each minireplicon they contain may determine the incompatibilities of megaplasmids. A broader analysis of all 1285 bacterial plasmids with putative known minireplicons whose complete genome sequences were available from GenBank revealed that 34% (443 plasmids) of the plasmids have two or more minireplicons. This indicates that plasmid fusion events are general among bacterial plasmids.
Megaplasmids of B. cereus group are fusion of smaller plasmids, and the fusion of plasmids likely occurs frequently in the B. cereus group and in other bacterial taxa. Plasmid fusion may be one of the major mechanisms for formation of novel megaplasmids in the evolution of bacteria.
Horizontal gene transfer (HGT) is the major driving force of bacterial evolution . Plasmids play important roles in this process via their conjugative capability . Additionally, plasmids harbor genes involved in niche specific processes, and are important for bacterial adaptation to changing environmental conditions [3, 4]. As plasmids can transfer frequently among different bacterial strains, they display strain-dependent distributions. Some bacterial strains containing no plasmids, while others have many; sometimes more than 20 . Moreover, the same host can harbor plasmids with a wide size range. For example, B. thuringiensis strain CT-43 has 10 plasmids ranging from 6 kb to 300 kb . However, the origin and evolution of these plasmids remains unclear. To date, studies on the evolution and dynamics have mainly focused on plasmids that have broad host ranges and harbor antibiotic-resistance (AR) genes, for example the plasmids of the incompatibility groups IncW , IncU , IncP [9, 10] and PromA . These plasmids usually have small genome sizes and few of them are larger than 100 kb. Information on the evolution and dynamics of plasmids that have relatively narrow host range is scarce, especially for megaplasmids larger than 100 kb. It was therefore the aim of this study to elucidate the origin, evolution and dynamics of megaplasmids with relatively narrow host range using the Bacillus cereus group as a model.
The B. cereus group includes B. anthracis, the causative agent of anthrax and a potential biological weapon; B. cereus, a ubiquitous soil bacterium and foodborne pathogen; B. thuringiensis, which produces insecticidal crystal proteins; and four additional species, B. cytotoxicus, B. mycoides, B. pseudomycoides, and B. weihenstephanensis[12, 13]. Strains of this group typically contain several plasmids, with some strains containing more than 10 [6, 14, 15]. Plasmids of this group are crucial for the phenotype and virotype of strains. B. anthracis, B. cereus, and B. thuringiensis were defined mainly on the basis of plasmid-encoded phenotypic features [16–18]. Usually, these plasmids are larger than 100 kb. Indeed, most strains of this group contain one or more megaplasmids larger than 100 kb. Only two of these megaplasmids have been studied in depth. One is pXO1 (182 kb) from B. anthracis which harbors two minireplicons that support replication of the plasmid: repX and pXO1-14/pXO1-16. A minireplicon represents the smallest replication region that supports plasmid replication, and contains the origin of replication and genes encoding replication proteins. The origin of replication of plasmid is a particular sequence in a plasmid genome at which replication is initiated. The other well-studied megaplasmid is B. thuringiensis plasmid pBtoxis (128 kb), whose minireplicon consists of two genes: orf156 and orf157. The availability of more than 30 sequences of megaplasmids in genomes of the B. cereus group allows the investigation of their evolution and dynamics.
We collected the genome sequences of plasmids for bioinformatic analyses. First, we studied the distribution of minireplicons for all the megaplasmids. Second, we studied the relationships among different megaplasmids from the same host strain, and from all strains of the B. cereus group. We also studied the distributions of known minireplicons among all plasmids outside of the B. cereus group whose genome sequences were available.
Results and discussion
Six types of minireplicon exist in the megaplasmids of the B. cereus group
The minireplicons are the core part of plasmids and drive plasmid replication and propagation. Their diversity and evolution directly reflects the dynamics and evolution of plasmids [7, 10, 22, 23]. Strains in the B. cereus group are rich in plasmid content, with plasmid numbers ranging from zero to 13 and sizes ranging from 2 kb to 600 kb [6, 15, 24]. Thus, the B. cereus group is an ideal model to study plasmid dynamics and evolution. This study aimed to characterize the origin, evolution and dynamics of megaplasmids with genome sizes larger than 100 kb by studying the distribution and evolution of their minireplicons.
Type and number of minireplicons on plasmids larger than 20 kb
Number of minireplicon
Plasmid size (bp)
rep466, rep26, pXO1-16/pXO1-14
ori44, repA_N, rep466, pXO1-14/pXO1-16(2)
rep466, rep26, pXO1-16/pXO1-14(2)
rep466, rep26, pXO1-16/pXO1-14
rep26 c , pXO1-16/pXO1-14
The other four minireplicons contain one essential protein each. Minireplicon rep26 was found in four megaplasmids larger than 100 kb. In contrast, ori44, ori60 and repA_N were found in both megaplasmids and plasmids smaller than 100 kb and are more widely distributed in smaller plasmids. For example, minireplicon ori44 occurs in eight plasmids, only two of which are larger than 100 kb, while minireplicon ori60 is contained by six plasmids, only one of which is larger than 100 kb. Although the two plasmids containing repA_N in this study are larger than 100 kb, plasmids from other Gram-positive bacteria that contain this type of minireplicon are usually smaller than 100 kb .
Megaplasmids larger than 100 kb contain two or more minireplicons in B. cereusgroup
In plasmids with more than one minireplicon, it is not known whether all the minireplicons are functional for plasmid replication, partitioning and maintenance. For plasmid pXO1, early studies confirmed that two different types of minireplicons, pXO1-14/pXO1-16 and repX, are functional for its replication. Recently, both the repX and pXO1-14/pXO1-16 minireplicons were proven to independently support replication of the full-length pXO1 plasmid, with pXO1-14/pXO1-16 being more effective than repX. Moreover, a 4848-bp DNA fragment within minireplicon pXO1-14/pXO1-16 can be used to deprive plasmid pXO1 from B. anthracis using plasmid incompatibility . This suggests that minireplicon pXO1-14/pXO1-16 is predominantly used for plasmid replication. We therefore speculate that when there is more than one minireplicon on the same plasmid, some of them are more relevant than others. However, how these minireplicons cooperate with each other is not clear.
Minireplicons are conserved during the evolutionary history of a plasmid; however, it would be interesting to determine the evolutionary relationship of multiple minireplicons on the same plasmid. To investigate this, we conducted a comparative analysis between minireplicons tubZ/tubR and pXO1-14/pXO1-16. First, we considered the relative position of the two minireplicons on the same plasmid. When there are only two types of minireplicon on one plasmid, the distance between the minireplicons ranged from 20 to 40 kb, and the distance between two minireplicons on larger plasmids is not larger (Spearman’s r = 0.17, P = 0.41). Multiple minireplicons are frequently clustered in a certain region of the plasmid, which can be recognized as the core region for replication and maintenance. Second, we compared the three phylogenetic trees that were constructed based on protein sequences of TubZ, pXO1-14-like and pXO1-16-like (Figures 1A and 2). The major topologies of these trees were inconsistent. The pXO1-14 or pXO1-16 trees cannot support the classification of four subtypes of TubZ. Plasmids with the same subtype of tubZ/tubR were usually found to have different subtypes of pXO1-14/pXO1-16. This indicates that different minireplicons on the same plasmid evolved independently.
Megaplasmids may be formed by fusion of smaller plasmids in B. cereusgroup
Analysis of the relationships between minireplicon types and plasmid sizes revealed that plasmids with one minireplicon are usually smaller than 100 kb. However, when two or more minireplicons were present on the same plasmid, the genome size could exceed that of either of the presumed original plasmids, usually larger than 100 kb (Table 1). For example, plasmids containing minireplicon ori44 only have genome sizes from 45 to 85 kb. In contrast, other plasmids that combine ori44 and one or more additional minireplicons usually have a genome size larger than 100 kb, even up to 417 kb for the plasmid pBWB401. This indicates that by integrating different minireplicons into a single plasmid, the new plasmid is capable of carrying more genes.
Larger plasmids have lower copy numbers than smaller ones . Formation of larger plasmids by fusion of smaller plasmids thus reduces the amount of DNA that is required for similar plasmid genome sizes. This could provide an evolutionary advantage by reducing the energy requirement for plasmid synthesis and maintenance. Moreover, plasmids with some minireplicons have very low copy numbers, and additional minireplicons are needed to support them to replicate effectively. For example, plasmid pXO1 with only the repX minireplicon was reported to have copy numbers ranging from 0.8 to 1. This indicates that this minireplicon cannot effectively support plasmid replication . If the plasmid contains another minireplicon, pXO1-14/pXO-16, in addition to the repX minireplicon, copy numbers ranging from 3 to 3.6 were observed and the plasmid is stably inherited. For those minireplicons that support effective plasmid replication, there may be dynamic equilibrium between the existence of small plasmids with individual minireplicons and the integration into megaplasmids with multiple minireplicons on plasmids. The selective forces driving plasmid evolution includes factors that determine the fitness of plasmids evolving as autonomous genetic elements as well factors that determine the added ecological fitness of the bacterial host. Ecological determinants that shape maintenance of small plasmids with one replicon or the integration into megaplasmids with multiple replicons are not clear.
Compatibility groups of megaplasmids may depend on each of their minireplicons at the subtype level
Megaplasmids contain more than one minireplicon and many strains contain more than one such megaplasmid, therefore, several minireplicons co-exist in the same host. To determine the compatibility of different minireplicons, we investigated the patterns of coexistence of minireplicons tubZ/tubR and pXO1-14/pXO1-16. As shown in Figure 1 and Additional file 3: Figure S4, each tubZ/tubR from the same strain belongs to one of the four subtypes, with different tubZs, tubRs and putative origins of replication. For example, the two tubZ/tubRs on the two megaplasmids pBC210 and pBCXO1 of B. cereus G9241 belong to the repX and rep228 subtypes, respectively. Many strains have more than one pXO1-14/pXO1-16 minireplicon. In most cases, each of their encoded pXO1-14 or pXO1-16 proteins from a certain strain was found to belong to different subgroups. For example, pXO1-14-like proteins encoded by plasmids pBC210 and pBCXO1 in B. cereus G9241 belong to subgroups I and III, respectively (Figure 2A). Their two corresponding pXO1-16-like proteins are also allocated to the two different subgroups, as shown by the pXO1-16-like protein tree (Figure 2B). In other instances, different plasmids in the same host contain pXO1-14-like or pXO-16-like proteins but only one of the two belongs to the same subgroup. For example, the two pXO1-16-like proteins of pBMB293 and pBMB400 from B. thuringiensis YBT-1520 are located on different branches of the same subgroup (Figure 2B), while the corresponding pXO1-14-like proteins show greater diversity and were allocated into subgroups II and IV, respectively. There was only one instance where different pXO1-14 and pXO1-16-like proteins from different minireplicons in the same strain were grouped together. Plasmid pAH1134_566 contains two pXO-14/16-like minireplicons and both proteins belong to the same subgroup. This may result from gene duplication and indicates that minireplicons of the same subgroups are compatible if they are located on the same plasmid.
All plasmids of a bacterial strain must be compatible. The minireplicon tubZ/tubR has four subtypes; thus, there may be four natural incompatibility groups for tubZ/tubR-containing megaplasmids in the B. cereus group. Different groups have different TubZs, TubRs and putative origins of replication. For pXO1-14/pXO1-16-like minireplicons, as the two essential proteins they encode do not have a concerted evolution, the putative incompatibility groups appear to be determined by the subgroup types of pXO1-14 or/and pXO1-16. Many plasmids contain both of these minireplicons; however, details regarding the coexistence of these plasmids are not clear.
Integrated events are general among plasmids during their evolutionary histories
Plasmids outside of the B. cereus group with more than one minireplicon have been reported, and the most frequently mentioned were plasmids belonging to incompatibility group F (IncF). Most plasmids of this group harbor two or more minireplicons, suggesting that plasmids fusion events occurred in the evolutionary histories of these plasmids . The direct example is plasmid pIP1206, which may have resulted from recombination between pRSB107 and a pAPEC-O1-ColBM-like plasmid. Among its 151 open reading frames, 56 (37%) were also present in pRSB107 and 44 (29%) in pAPEC-O1-ColBM (24) .
In order to reduce the effect of data bias on the results, we used a subset of the plasmid sequence data to repeat the analysis. For each species that has plasmid genome sequences reported, we chose all plasmids from one strain whose plasmid number is the largest in that species. Analysis of this subset of the plasmid genome sequences confirmed the results obtained with the entire data set. Among the 771 plasmids with putative minireplicons (Additional file 5: Table S4), one third of the plasmids have two or more minireplicons (Figure 5C) and plasmids with two or more minireplicons are larger than those with only one minireplicon (Figure 5D, P = 1.035e-06, Mann–Whitney test).
We found that megaplasmids in the B. cereus group larger than 100 kb contain two or more minireplicons. Minireplicons on the same plasmid usually have distinct evolutionary histories. We hypothesize that these megaplasmids are fusions of smaller plasmids. About one third of the plasmids out of the B. cereus group have multiple minireplicons. This indicates that plasmids fusion events occur generally during the plasmids evolutionary histories and plasmids fusion may be an important mechanism for the formation of megaplasmids.
The genome sequences of 45 plasmids of the B. cereus group were retrieved from GenBank (http://www.ncbi.nlm.nih.gov) and those of 11 unpublished plasmids sequenced by our group were used in the analyses. The genome sizes of these 56 plasmids ranged from ≈20 kb to ≈566 kb (Additional file 1: Table S1). All of these data (Dataset 1) were obtained by October 20, 2012.
To study the minireplicons across all of the prokaryotic species, we collected all 3340 plasmid genome sequences from Genbank ftp site (ftp://ftp.ncbi.nlm.nih.gov/genomes/Plasmids/). These data (Dataset 2) were obtained by February 10, 2013.
Replication essential protein sequences and minireplicons prediction
TubZ protein sequences were obtained from Dataset 1 using the hmmsearch command of the hmmer version 3.0 software , with an e-value<0.001, and the model Tubulin/FtsZ family (PF00091) were obtained from the Pfam database . Other types of replicated protein sequences were obtained by BLASTP analysis  using various types of reported replication protein sequences from B. cereus group plasmids as query sequences and the non-redundant protein sequences from Dataset 1 as the database. A minireplicon was approved when all of the essential elements, including one or two genes encoding replication essential proteins and the DNA fragment containing origin of replication, were predicted.
We looked for replicated protein sequences from Dataset 2 by two methods; one was searching the keywords such as “replication protein”, “Rep protein” or “Primase” from the annotation files, and the other one was using hmmsearch command of hmmer software with the models associated with plasmid replication (Additional file 6: Table S2) which were downloaded from Pfam database . Then we combined the results from both of the above methods, and checked these results based on public information. Minireplicons were approved when all the essential elements were predicted. All of the 1285 plasmids with putative minireplicon were showed in Additional file 4: Table S3.
Sequence alignment and phylogenetic analysis
Protein sequences for different minireplicon of B. cereus group were aligned using Muscle . The most disordered regions were eliminated using G-blocks . The evolutionary models that best fit these sequences were determined by ProtTest version 3.0 , and Maximum Likelihood (ML) phylogenetic trees were generated by PhyML software version 3.0 , using the best fitted models (JTT + G + F for pXO1-14 and pXO1-16, LG + G + F for TubZ). Bootstrap supports were calculated as a percent of 1000 replicates. As the identity levels of TubR protein sequences among each type are very high, we collected the DNA sequences from them. Each type of tubR DNA sequences was aligned by Muscle and a ML tree was constructed using PhyML based on the model determined by ModelTest . All the phylogenetic trees were deposited in treeBASE .
All statistical analyses were carried out using in-house Perl scripts and R 2.15.1 .
We want to thank Michael Gänzle from University of Alberta for his critical reading of the manuscript. This work was supported by grants from the National High Technology Research and Development Program (863) of China (2011AA10A203), China 948 Program of Ministry of Agriculture (2011-G25), the National Basic Research Program (973) of China (2009CB118902), the National Natural Science Foundation of China (31170047 and 31000020), the international scientific cooperation of Hubei province (2011BFA019), and the foundmental research fund for the central university (2011PY056).
- Boto L: Horizontal gene transfer in evolution: facts and challenges. Proc Biol Sci. 2010, 277 (1683): 819-827. 10.1098/rspb.2009.1679.PubMed CentralPubMedView ArticleGoogle Scholar
- Norman A, Hansen LH, Sorensen SJ: Conjugative plasmids: vessels of the communal gene pool. Philos Trans R Soc Lond B Biol Sci. 2009, 364 (1527): 2275-2289. 10.1098/rstb.2009.0037.PubMed CentralPubMedView ArticleGoogle Scholar
- Bennett PM: Plasmid encoded antibiotic resistance: acquisition and transfer of antibiotic resistance genes in bacteria. Br J Pharmacol. 2008, 153 (Suppl 1): S347-357.PubMed CentralPubMedGoogle Scholar
- Liu XY, Ruan LF, Hu ZF, Peng DH, Cao SY, Yu ZN, Liu Y, Zheng JS, Sun M: Genome-wide screening reveals the genetic determinants of an antibiotic insecticide in Bacillus thuringiensis. J Biol Chem. 2010, 285 (50): 39191-39200. 10.1074/jbc.M110.148387.PubMed CentralPubMedView ArticleGoogle Scholar
- Schutzer SE, Fraser-Liggett CM, Casjens SR, Qiu WG, Dunn JJ, Mongodin EF, Luft BJ: Whole-genome sequences of thirteen isolates of Borrelia burgdorferi. J Bacteriol. 2011, 193 (4): 1018-1020. 10.1128/JB.01158-10.PubMed CentralPubMedView ArticleGoogle Scholar
- He J, Wang J, Yin W, Shao X, Zheng H, Li M, Zhao Y, Sun M, Wang S, Yu Z: Complete genome sequence of Bacillus thuringiensis subsp. chinensis strain CT-43. J Bacteriol. 2011, 193 (13): 3407-3408. 10.1128/JB.05085-11.PubMed CentralPubMedView ArticleGoogle Scholar
- Fernandez-Lopez R, Garcillan-Barcia MP, Revilla C, Lazaro M, Vielva L, de la Cruz F: Dynamics of the IncW genetic backbone imply general trends in conjugative plasmid evolution. FEMS Microbiol Rev. 2006, 30 (6): 942-966. 10.1111/j.1574-6976.2006.00042.x.PubMedView ArticleGoogle Scholar
- Rhodes G, Parkhill J, Bird C, Ambrose K, Jones MC, Huys G, Swings J, Pickup RW: Complete nucleotide sequence of the conjugative tetracycline resistance plasmid pFBAOT6, a member of a group of IncU plasmids with global ubiquity. Appl Environ Microbiol. 2004, 70 (12): 7497-7510. 10.1128/AEM.70.12.7497-7510.2004.PubMed CentralPubMedView ArticleGoogle Scholar
- Norberg P, Bergstrom M, Jethava V, Dubhashi D, Hermansson M: The IncP-1 plasmid backbone adapts to different host bacterial species and evolves through homologous recombination. Nat Commun. 2011, 2: 268-PubMed CentralPubMedView ArticleGoogle Scholar
- Sen D, Brown CJ, Top EM, Sullivan J: Inferring the Evolutionary History of IncP-1 Plasmids Despite Incongruence among Backbone Gene Trees. Mol Biol Evol. 2013, 30 (1): 154-166. 10.1093/molbev/mss210.PubMed CentralPubMedView ArticleGoogle Scholar
- Van der Auwera GA, Krol JE, Suzuki H, Foster B, Van Houdt R, Brown CJ, Mergeay M, Top EM: Plasmids captured in C. metallidurans CH34: defining the PromA family of broad-host-range plasmids. Antonie Van Leeuwenhoek. 2009, 96 (2): 193-204. 10.1007/s10482-009-9316-9.PubMedView ArticleGoogle Scholar
- Vilas-Boas GT, Peruca AP, Arantes OM: Biology and taxonomy of Bacillus cereus, Bacillus anthracis, and Bacillus thuringiensis. Can J Microbiol. 2007, 53 (6): 673-687. 10.1139/W07-029.PubMedView ArticleGoogle Scholar
- Guinebretiere MH, Auger S, Galleron N, Contzen M, De Sarrau B, De Buyser ML, Lamberet G, Fagerlund A, Granum PE, Lereclus D, et al: Bacillus cytotoxicus sp. nov. is a novel thermotolerant species of the Bacillus cereus Group occasionally associated with food poisoning. Int J Syst Evol Microbiol. 2013, 63 (Pt 1): 31-40. 10.1099/ijs.0.030627-0.PubMedView ArticleGoogle Scholar
- Reyes-Ramirez A, Ibarra JE: Plasmid patterns of Bacillus thuringiensis type strains. Appl Environ Microbiol. 2008, 74 (1): 125-129. 10.1128/AEM.02133-07.PubMed CentralPubMedView ArticleGoogle Scholar
- Zhong C, Peng D, Ye W, Chai L, Qi J, Yu Z, Ruan L, Sun M: Determination of plasmid copy number reveals the total plasmid DNA amount is greater than the chromosomal DNA amount in Bacillus thuringiensis YBT-1520. PLoS One. 2011, 6 (1): e16025-10.1371/journal.pone.0016025.PubMed CentralPubMedView ArticleGoogle Scholar
- Kolsto AB, Tourasse NJ, Okstad OA: What sets Bacillus anthracis apart from other Bacillus species?. Annu Rev Microbiol. 2009, 63: 451-476. 10.1146/annurev.micro.091208.073255.PubMedView ArticleGoogle Scholar
- Ehling-Schulz M, Fricker M, Grallert H, Rieck P, Wagner M, Scherer S: Cereulide synthetase gene cluster from emetic Bacillus cereus: structure and location on a mega virulence plasmid related to Bacillus anthracis toxin plasmid pXO1. BMC Microbiol. 2006, 6: 20-10.1186/1471-2180-6-20.PubMed CentralPubMedView ArticleGoogle Scholar
- Berry C, O'Neil S, Ben-Dov E, Jones AF, Murphy L, Quail MA, Holden MT, Harris D, Zaritsky A, Parkhill J: Complete sequence and organization of pBtoxis, the toxin-coding plasmid of Bacillus thuringiensis subsp. israelensis. Appl Environ Microbiol. 2002, 68 (10): 5082-5095. 10.1128/AEM.68.10.5082-5095.2002.PubMed CentralPubMedView ArticleGoogle Scholar
- Tinsley E, Khan SA: A novel FtsZ-like protein is involved in replication of the anthrax toxin-encoding pXO1 plasmid in Bacillus anthracis. J Bacteriol. 2006, 188 (8): 2829-2835. 10.1128/JB.188.8.2829-2835.2006.PubMed CentralPubMedView ArticleGoogle Scholar
- Pomerantsev AP, Camp A, Leppla SH: A new minimal replicon of Bacillus anthracis plasmid pXO1. J Bacteriol. 2009, 191 (16): 5134-5146. 10.1128/JB.00422-09.PubMed CentralPubMedView ArticleGoogle Scholar
- Tang M, Bideshi DK, Park HW, Federici BA: Minireplicon from pBtoxis of Bacillus thuringiensis subsp. israelensis. Appl Environ Microbiol. 2006, 72 (11): 6948-6954. 10.1128/AEM.00976-06.PubMed CentralPubMedView ArticleGoogle Scholar
- Petersen J, Brinkmann H, Berger M, Brinkhoff T, Pauker O, Pradella S: Origin and Evolution of a Novel DnaA-Like Plasmid Replication Type in Rhodobacterales. Molecular Biology and Evolution. 2010, 28 (3): 1229-1240.PubMedView ArticleGoogle Scholar
- Weaver KE, Kwong SM, Firth N, Francia MV: The RepA_N replicons of Gram-positive bacteria: a family of broadly distributed but narrow host range plasmids. Plasmid. 2009, 61 (2): 94-109. 10.1016/j.plasmid.2008.11.004.PubMed CentralPubMedView ArticleGoogle Scholar
- Challacombe JF, Altherr MR, Xie G, Bhotika SS, Brown N, Bruce D, Campbell CS, Campbell ML, Chen J, Chertkov O, et al: The complete genome sequence of Bacillus thuringiensis Al Hakam. J Bacteriol. 2007, 189 (9): 3680-3681. 10.1128/JB.00241-07.PubMed CentralPubMedView ArticleGoogle Scholar
- Baum JA, Gilbert MP: Characterization and comparative sequence analysis of replication origins from three large Bacillus thuringiensis plasmids. J Bacteriol. 1991, 173 (17): 5280-5289.PubMed CentralPubMedGoogle Scholar
- Akhtar P, Khan SA: Two independent replicons can support replication of the anthrax toxin-encoding plasmid pXO1 of Bacillus anthracis. Plasmid. 2012, 67 (2): 111-117. 10.1016/j.plasmid.2011.12.012.PubMed CentralPubMedView ArticleGoogle Scholar
- Liu X, Wang D, Wang H, Feng E, Zhu L: Curing of plasmid pXO1 from Bacillus anthracis using plasmid incompatibility. PLoS One. 2012, 7 (1): e29875-10.1371/journal.pone.0029875.PubMed CentralPubMedView ArticleGoogle Scholar
- Villa L, Garcia-Fernandez A, Fortini D, Carattoli A: Replicon sequence typing of IncF plasmids carrying virulence and resistance determinants. J Antimicrob Chemother. 2010, 65 (12): 2518-2529. 10.1093/jac/dkq347.PubMedView ArticleGoogle Scholar
- Perichon B, Bogaerts P, Lambert T, Frangeul L, Courvalin P, Galimand M: Sequence of conjugative plasmid pIP1206 mediating resistance to aminoglycosides by 16S rRNA methylation and to hydrophilic fluoroquinolones by efflux. Antimicrob Agents Chemother. 2008, 52 (7): 2581-2592. 10.1128/AAC.01540-07.PubMed CentralPubMedView ArticleGoogle Scholar
- Finn RD, Clements J, Eddy SR: HMMER web server: interactive sequence similarity searching. Nucleic Acids Res. 2011, 39 (Web Server issue): W29-37-PubMedGoogle Scholar
- Punta M, Coggill PC, Eberhardt RY, Mistry J, Tate J, Boursnell C, Pang N, Forslund K, Ceric G, Clements J, et al: The Pfam protein families database. Nucleic Acids Res. 2012, 40 (Database issue): D290-301.PubMed CentralPubMedView ArticleGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25 (17): 3389-3402. 10.1093/nar/25.17.3389.PubMed CentralPubMedView ArticleGoogle Scholar
- Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32 (5): 1792-1797. 10.1093/nar/gkh340.PubMed CentralPubMedView ArticleGoogle Scholar
- Castresana J: Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol Biol Evol. 2000, 17 (4): 540-552. 10.1093/oxfordjournals.molbev.a026334.PubMedView ArticleGoogle Scholar
- Darriba D, Taboada GL, Doallo R, Posada D: ProtTest 3: fast selection of best-fit models of protein evolution. Bioinformatics. 2011, 27 (8): 1164-1165. 10.1093/bioinformatics/btr088.PubMedView ArticleGoogle Scholar
- Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O: New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol. 2010, 59 (3): 307-321. 10.1093/sysbio/syq010.PubMedView ArticleGoogle Scholar
- Posada D: ModelTest Server: a web-based tool for the statistical selection of models of nucleotide substitution online. Nucleic Acids Res. 2006, 34 (Web Server issue): W700-703.PubMed CentralPubMedView ArticleGoogle Scholar
- Anwar N, Hunt E: Improved data retrieval from TreeBASE via taxonomic and linguistic data enrichment. BMC Evol Biol. 2009, 9: 93-10.1186/1471-2148-9-93.PubMed CentralPubMedView ArticleGoogle Scholar
- R Core Team: R: A Language and Environment for Statistical Computing. 2012, Vienna, Austria: R Foundation for Statistical ComputingGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.