OASIS: an automated program for global investigation of bacterial and archaeal insertion sequences

Insertion sequences (ISs) are simple transposable elements present in most bacterial and archaeal genomes and play an important role in genomic evolution. The recent expansion of sequenced genomes offers the opportunity to study ISs comprehensively, but this requires efficient and accurate tools for IS annotation. We have developed an open-source program called OASIS, or Optimized Annotation System for Insertion Sequences, which automatically annotates ISs within sequenced genomes. OASIS annotations of 1737 bacterial and archaeal genomes offered an unprecedented opportunity to examine IS evolution. At a broad scale, we found that most IS families are quite widespread; however, they are not present randomly across taxa. This may indicate differential loss, barriers to exchange and/or insufficient time to equilibrate across clades. The number of ISs increases with genome length, but there is both tremendous variation and no increase in IS density for genomes >2 Mb. At the finer scale of recently diverged genomes, the proportion of shared IS content falls sharply, suggesting loss and/or emergence of barriers to successful cross-infection occurs rapidly. Surprisingly, even after controlling for 16S rRNA sequence divergence, the same ISs were more likely to be shared between genomes labeled as the same species rather than as different species.     

A survey of bacterial insertion sequences using IScan

Bacterial insertion sequences (ISs) are the simplest kinds of bacterial mobile DNA. Evolutionary studies need consistent IS annotation across many different genomes. We have developed an open-source software package, IScan, to identify bacterial ISs and their sequence elements—inverted and target direct repeats—in multiple genomes using multiple flexible search parameters. We applied IScan to 438 completely sequenced bacterial genomes and 20 IS families. The resulting data show that ISs within a genome are extremely similar, with a mean synonymous divergence of K_s = 0.033. Our analysis substantially extends previously available information, and suggests that most ISs have entered bacterial genomes recently. By implication, their population persistence may depend on horizontal transfer. We also used IScan's ability to analyze the statistical significance of sequence similarity among many IS inverted repeats. Although the inverted repeats of insertion sequences are evolutionarily highly flexible parts of ISs, we show that this ability can be used to enrich a dataset for ISs that are likely to be functional. Applied to the thousands of genomes that will soon be available, IScan could be used for many purposes, such as mapping the evolutionary history and horizontal transfer patterns of different ISs.     

Identification and distribution of new insertion sequences in the genome of the extremely halotolerant and alkaliphilic Oceanobacillus iheyensis HTE831.

Six kinds of new insertion sequences (ISs), IS667 to IS672, a group II intron (Oi.Int), and an incomplete transposon (Tn852loi) were identified in the 3,630,528-bp genome of the extremely halotolerant and alkaliphilic Oceanobacillus iheyensis HTE831. Of 19 ISs identified in the HTE831 genome, 7 were truncated, indicating the occurrence of internal rearrangement of the genome. All ISs except IS669 generated a 4- to 8-bp duplication of the target site sequence, and these ISs carried 23- to 28-bp inverted repeats (IRs). Sequence analysis revealed that four ISs (IS669, IS670, IS671, and IS672) were newly identified as belonging to separate IS families (IS200/IS605, IS30, IS5, and IS3, respectively). IS667 and IS668 were also characterized as new members of the ISL3 family. Tn8521oi, which belongs to the Tn3 family as a new member, generated a 5-bp duplication of the target site sequence and carried complete 38-bp IRs. Of the eight protein-coding sequences (CDSs) identified in Tn8521oi, three CDSs (OB481, OB482, and OB483) formed a ger gene cluster, and two other paralogous gene clusters were found in the HTE831 genome. Most of the ISs and the group II intron widely distributed throughout the genome were inserted in noncoding regions, while two ISs (IS667-08 and IS668-02) and Oi.Int-04 were inserted in the coding regions.     

Causes of insertion sequences abundance in prokaryotic genomes

Insertion sequences (ISs) are the smallest and most frequent transposable elements in prokaryotes where they play an important evolutionary role by promoting gene inactivation and genome plasticity. Their genomic abundance varies by several orders of magnitude for reasons largely unknown and widely speculated. The current availability of hundreds of genomes renders testable many of these hypotheses, notably that IS abundance correlates positively with the frequency of horizontal gene transfer (HGT), genome size, pathogenicity, nonobligatory ecological associations, and human association. We thus reannotated ISs in 262 prokaryotic genomes and tested these hypotheses showing that when using appropriate controls, there is no empirical basis for IS family specificity, pathogenicity, or human association to influence IS abundance or density. HGT seems necessary for the presence of ISs, but cannot alone explain the absence of ISs in more than 20% of the organisms, some of which showing high rates of HGT. Gene transfer is also not a significant determinant of the abundance of IS elements in genomes, suggesting that IS abundance is controlled at the level of transposition and ensuing natural selection and not at the level of infection. Prokaryotes engaging in obligatory associations have fewer ISs when controlled for genome size, but this may be caused by some being sexually isolated. Surprisingly, genome size is the only significant predictor of IS numbers and density. Alone, it explains over 40% of the variance of IS abundance. Because we find that genome size and IS abundance correlate negatively with minimal doubling times, we conclude that selection for rapid replication cannot account for the few ISs found in small genomes. Instead, we show evidence that IS numbers are controlled by the frequency of highly deleterious insertion targets. Indeed, IS abundance increases quickly with genome size, which is the exact inverse trend found for the density of genes under strong selection such as essential genes. Hence, for ISs, the bigger the genome the better.     

Identification and distribution of new insertion sequences in the genome of alkaliphilic Bacillus halodurans C-125.

Fifteen kinds of new insertion sequences (ISs), IS641 to IS643, IS650 to IS658, IS660, IS662, and IS663, and a group II intron (Bh.Int) were identified in the 4,202,352-bp genome of alkaliphilic Bacillus halodurans C-125. Out of 120 ISs identified in the C-125 genome, 29 were truncated, indicating the occurrence of internal rearrangements of the genome. The ISs other than IS650, IS653, IS660, and IS663 generated a 2- to 9-bp duplication of the target site sequence, and the ISs other than IS650, IS653, and IS657 carry 14- to 64-bp inverted repeats. Sequence analysis revealed that six kinds of ISs (IS642, IS643, IS654, IS655, IS657, and IS658) belong to a separate IS family (IS630, IS21, IS256, IS3, IS200/IS605, and IS30, respectively) as a new member. Also, IS651 and IS652 were characterized as new members of the ISL3 family. Significant similarity was found between the transposase (Tpase) sequences between IS650 and IS653 (78.2%), IS651 and IS652 (56.3%), IS656 and IS662 (71.0%), and IS660 and IS663 (44.5%), but the others showed no similarity to one another. Tpases in 28 members of IS651 in the C-125 genome were found to have become diversified. Most of the IS elements widely distributed throughout the genome were inserted in noncoding regions, although some genes, such as those coding for an ATP-binding cassette transporter/permease, a response regulator, and L-indole 2-dehydrogenase, have been mutated through the insertion of IS elements. It is evident, however, that not all IS elements have transposed and caused rearrangements of the genome in the past 17 years during which strain C-125 was subcultured under neutral and alkaline conditions.     

Experimental analysis of recently transposed insertion sequences in the cyanobacterium Synechocystis sp. PCC 6803.

The genome DNA of the cyanobacterium Synechocystis sp. PCC 6803 carries a number of insertion sequences (Kaneko, T. et al. 1996, DNA Res., 3, 109-136). We analyzed one of the abundant ISs (ISY203 group of IS4 family) in the common three substrains of Synechocystis and found that the four ISs with identical nucleotide sequences were present only in the "Kazusa" strain, whose complete genome sequence had been determined, while absent in ancestral strains (the original strain from Pasteur Culture Collection and its glucose-tolerant derivative). Three of these ISs were found in the genomic sequence as transposase genes of sll1474, sll1780 and slr1635. The fourth was on the plasmid, pSYSM. On the other hand, all three strains had a novel IS (denoted ISY203x), of which the nucleotide sequence was totally identical to the four ISs found only in the Kazusa strain. Since the flanking regions of ISY203x did not match any part of the genome or of the known plasmids of Synechocystis, it is presumably located on a yet uncharacterized plasmid. These suggest that the four ISs in Kazusa strain were recently transposed from ISY203x. Apparently, the transposition inactivated four preexisting genes, of which modified forms are presented as putative genes (sll1473, sll1475, slr1862, slr1863, slr1635 and ssl2982) in the list of the complete genome (CyanoBase: http://www.kazusa.or.jp/cyano/cyano.html). The possible effects of transposition of ISs in Synechocystis are discussed in relation to phenotypic mutations and microevolution.     

Analysis of insertion sequences in thermophilic cyanobacteria: exploring the mechanisms of establishing, maintaining, and withstanding high insertion sequence abundance

Insertion sequences (ISs) are simple mobile genetic elements capable of relocating within a genome. Through this transposition activity, they are known to create mutations which are mostly deleterious to the cell, although occasionally they are beneficial. Two closely related isolates of thermophilic Synechococcus species from hot spring microbial mats are known to harbor a large number of diverse ISs. To explore the mechanism of IS acquisition within natural populations and survival in the face of high IS abundance, we examined IS content and location in natural populations of Synechococcus by comparing metagenomic data to the genomes of fully sequenced cultured isolates. The observed IS distribution in the metagenome was equivalent to the distribution in the isolates, indicating that the cultured isolates are appropriate models for the environmental population. High sequence conservation between IS families shared between the two isolates suggests that ISs are able to move between individuals within populations and between species via lateral gene transfer, consistent with models for IS family accumulation. Most IS families show evidence of recent activity, and interruption of critical genes in some individuals was observed, demonstrating that transposition is an ongoing mutational force in the populations.     

Local hopping mobile DNA implicated in pseudogene formation and reductive evolution in an obligate cyanobacteria-plant symbiosis

Background

Insertion sequences (ISs) are approximately 1 kbp long “jumping” genes found in prokaryotes. ISs encode the protein Transposase, which facilitates the excision and reinsertion of ISs in genomes, making these sequences a type of class I (“cut-and-paste”) Mobile Genetic Elements. ISs are proposed to be involved in the reductive evolution of symbiotic prokaryotes. Our previous sequencing of the genome of the cyanobacterium ‘Nostoc azollae’ 0708, living in a tight perpetual symbiotic association with a plant (the water fern Azolla), revealed the presence of an eroding genome, with a high number of insertion sequences (ISs) together with an unprecedented large proportion of pseudogenes. To investigate the role of ISs in the reductive evolution of ‘Nostoc azollae’ 0708, and potentially in the formation of pseudogenes, a bioinformatic investigation of the IS identities and positions in 47 cyanobacterial genomes was conducted. To widen the scope, the IS contents were analysed qualitatively and quantitatively in 20 other genomes representing both free-living and symbiotic bacteria.

Results

Insertion Sequences were not randomly distributed in the bacterial genomes and were found to transpose short distances from their original location (“local hopping”) and pseudogenes were enriched in the vicinity of IS elements. In general, symbiotic organisms showed higher densities of IS elements and pseudogenes than non-symbiotic bacteria. A total of 1108 distinct repeated sequences over 500 bp were identified in the 67 genomes investigated. In the genome of ‘Nostoc azollae’ 0708, IS elements were apparent at 970 locations (14.3%), with 428 being full-length. Morphologically complex cyanobacteria with large genomes showed higher frequencies of IS elements, irrespective of life style.

Conclusions

The apparent co-location of IS elements and pseudogenes found in prokaryotic genomes implies earlier IS transpositions into genes. As transpositions tend to be local rather than genome wide this likely explains the proximity between IS elements and pseudogenes. These findings suggest that ISs facilitate the reductive evolution in for instance in the symbiotic cyanobacterium ‘Nostoc azollae’ 0708 and in other obligate prokaryotic symbionts.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1386-7) contains supplementary material, which is available to authorized users.     

Insertion sequences in prokaryotic genomes

Insertion sequences (ISs) are small DNA segments that are often capable of moving neighbouring genes. Over 1500 different ISs have been identified to date. They can have large and spectacular effects in shaping and reshuffling the bacterial genome. Recent studies have provided dramatic examples of such IS activity, including massive IS expansion during the emergence of some pathogenic bacterial species and the intimate involvement of ISs in assembling genes into complex plasmid structures. However, a global understanding of their impact on bacterial genomes requires detailed knowledge of their distribution across the eubacterial and archaeal kingdoms, understanding their partition between chromosomes and extra-chromosomal elements (e.g. plasmids and viruses) and the factors which influence this, and appreciation of the different transposition mechanisms in action, the target preferences and the host factors that influence transposition. In addition, defective (non- autonomous) elements, which can be complemented by related active elements in the same cell, are often overlooked in genome annotations but also contribute to the evolution of genome organisation.     

Molecular features of Mycobacterium avium human isolates carrying a single copy of IS1245 and IS1311 per genome

Human clinical isolates of the Mycobacterium avium complex, from hospitals in Bogotá, were studied using a wide range of molecular tests including PCR restriction-enzyme analysis (PRA) of the hsp65 gene. Up to 21 of the isolates were identified as M. avium PRA variant III (Mav III), a variant obtained only from isolates on the American continent. In contrast to previous reports, restriction fragment length polymorphism analysis using IS1245 and IS1311 showed a single copy for each insertion sequence (IS) in the majority (19/21) of the Colombian Mav III isolates under study. In order to analyse whether the ISs were inserted in a relevant genomic region, experimental conditions were established to determine the insertion loci of each single copy of both ISs in the genome. Analysis of genomic insertion loci indicated that both IS1245 and IS1311 were present in areas containing putatively truncated integrases and/or transposases, which may have an influence on the mobility of the inserted IS. In addition, a conserved genomic region was identified for the insertion of IS1311; this region could be part of the IS1311 itself.     

IS<Emphasis Type="Italic">4</Emphasis> family goes genomic

Background  

Insertion sequences (ISs) are small, mobile DNA entities able to expand in prokaryotic genomes and trigger important rearrangements. To understand their role in evolution, accurate IS taxonomy is essential. The IS 4 family is composed of ~70 elements and, like some other families, displays extremely elevated levels of internal divergence impeding its classification. The increasing availability of complete genome sequences provides a valuable source for the discovery of additional IS 4 elements. In this study, this genomic database was used to update the structural and functional definition of the IS 4 family.     

Insertion sequence diversity in archaea.

Insertion sequences (ISs) can constitute an important component of prokaryotic (bacterial and archaeal) genomes. Over 1,500 individual ISs are included at present in the ISfinder database (www-is.biotoul.fr), and these represent only a small portion of those in the available prokaryotic genome sequences and those that are being discovered in ongoing sequencing projects. In spite of this diversity, the transposition mechanisms of only a few of these ubiquitous mobile genetic elements are known, and these are all restricted to those present in bacteria. This review presents an overview of ISs within the archaeal kingdom. We first provide a general historical summary of the known properties and behaviors of archaeal ISs. We then consider how transposition might be regulated in some cases by small antisense RNAs and by termination codon readthrough. This is followed by an extensive analysis of the IS content in the sequenced archaeal genomes present in the public databases as of June 2006, which provides an overview of their distribution among the major archaeal classes and species. We show that the diversity of archaeal ISs is very great and comparable to that of bacteria. We compare archaeal ISs to known bacterial ISs and find that most are clearly members of families first described for bacteria. Several cases of lateral gene transfer between bacteria and archaea are clearly documented, notably for methanogenic archaea. However, several archaeal ISs do not have bacterial equivalents but can be grouped into Archaea-specific groups or families. In addition to ISs, we identify and list nonautonomous IS-derived elements, such as miniature inverted-repeat transposable elements. Finally, we present a possible scenario for the evolutionary history of ISs in the Archaea.     

An insertion sequence-dependent plasmid rearrangement in Aeromonas salmonicida causes the loss of the type three secretion system

Aeromonas salmonicida, a bacterial fish pathogen, possesses a functional Type Three Secretion System (TTSS), which is essential for its virulence. The genes for this system are mainly located in a single region of the large pAsa5 plasmid. Bacteria lose the TTSS region from this plasmid through rearrangements when grown in stressful growth conditions. The A. salmonicida genome is rich in insertion sequences (ISs), which are mobile DNA elements that can cause DNA rearrangements in other bacterial species. pAsa5 possesses numerous ISs. Three IS11s from the IS256 family encircle the rearranged regions. To confirm that these IS11s are involved in pAsa5 rearrangements, 26 strains derived from strain A449 and two Canadian isolates (01-B526 and 01-B516) with a pAsa5 rearrangement were tested using a PCR approach to determine whether the rearrangements were the result of an IS11-dependent process. Nine out of the 26 strains had a positive PCR result, suggesting that the rearrangement in these strains were IS-dependent. The PCR analysis showed that all the rearrangements in the A449-derived strains were IS11-dependent process while the rearrangements in 01-B526 and 01-B516 could only be partially coupled to the action of IS11. Unidentified elements that affect IS-dependent rearrangements may be present in 01-B526 and 01-B516. Our results suggested that pAsa5 rearrangements involve IS11. This is the first study showing that ISs are involved in plasmid instability in A. salmonicida.     

Positive selection on transposase genes of insertion sequences in the Crocosphaera watsonii genome

Insertion sequences (ISs) are mobile elements that are commonly found in bacterial genomes. Here, the structural and functional diversity of these mobile elements in the genome of the cyanobacterium Crocosphaera watsonii WH8501 is analyzed. The number, distribution, and diversity of nucleotide and amino acid stretches with similarity to the transposase gene of this IS family suggested that this genome harbors many functional as well as truncated IS fragments. The selection pressure acting on full-length transposase open reading frames of these ISs suggested (i) the occurrence of positive selection and (ii) the presence of one or more positively selected codons. These results were obtained using three data sets of transposase genes from the same IS family that were collected based on the level of amino acid similarity, the presence of an inverted repeat, and the number of sequences in the data sets. Neither recombination nor ribosomal frameshifting, which may interfere with the selection analyses, appeared to be important forces in the transposase gene family. Some positively selected codons were located in a conserved domain, suggesting that these residues are functionally important. The finding that this type of selection acts on IS-carried genes is intriguing, because although ISs have been associated with the adaptation of the bacterial host to new environments, this has typically been attributed to transposition or transformation, thus involving different genomic locations. Intragenic adaptation of IS-carried genes identified here may constitute a novel mechanism associated with bacterial diversification and adaptation.     

The distribution of insertion sequences in the genome of Shigella flexneri strain 2457T

Shigella flexneri, which causes shigellosis in humans, evolved from Escherichia coli. The sequencing of Shigella genomes has revealed that a large number of insertion sequence (IS) elements (over 200 elements) reside in the genome. Although the presence of these elements has been noted previously and summarized, more detailed analyses are required to understand their evolutionary significance. Here, the genome of S. flexneri strain 2457T is used to investigate the spatial distribution of IS copies around the chromosome and the location of elements with respect to genes. It is found that most IS isoforms occur essentially randomly around the genome. Two exceptions are IS91 and IS911, which appear to cluster due to local hopping. The location of IS elements with respect to genes is biased, however, revealing the action of natural selection. The non-coding regions of the genome (no more than 21%) carry disproportionally more IS elements (at least 28%) than the coding regions, implying that selection acts against insertion into genes. Of the genes disrupted by ISs, those involved in signal transduction, intracellular trafficking, and cell motility are most commonly targeted, suggesting selection against genes in these categories.     

Accessing the SEED Genome Databases via Web Services API: Tools for Programmers

Background  

The SEED integrates many publicly available genome sequences into a single resource. The database contains accurate and up-to-date annotations based on the subsystems concept that leverages clustering between genomes and other clues to accurately and efficiently annotate microbial genomes. The backend is used as the foundation for many genome annotation tools, such as the Rapid Annotation using Subsystems Technology (RAST) server for whole genome annotation, the metagenomics RAST server for random community genome annotations, and the annotation clearinghouse for exchanging annotations from different resources. In addition to a web user interface, the SEED also provides Web services based API for programmatic access to the data in the SEED, allowing the development of third-party tools and mash-ups.     

New insertion sequences of Sulfolobus: functional properties and implications for genome evolution in hyperthermophilic archaea

Analyses of complete genomes indicate that insertion sequences (ISs) are abundant and widespread in hyperthermophilic archaea, but few experimental studies have measured their activities in these hosts. As a way to investigate the impact of ISs on Sulfolobus genomes, we identified seven transpositionally active ISs in a widely distributed Sulfolobus species, and measured their functional properties. Six of the seven were found to be distinct from previously described ISs of Sulfolobus, and one of the six could not be assigned to any known IS family. A type II 'Miniature Inverted-repeat Transposable Element' (MITE) related to one of the ISs was also recovered. Rates of transposition of the different ISs into the pyrEF region of their host strains varied over a 250-fold range. The Sulfolobus ISs also differed with respect to target-site selectivity, although several shared an apparent preference for the pyrEF promoter region. Despite the number of distinct ISs assayed and their molecular diversity, only one demonstrated precise excision from the chromosomal target region. The fact that this IS is the only one lacking inverted repeats and target-site duplication suggests that the observed precise excision may be promoted by the IS itself. Sequence searches revealed previously unidentified partial copies of the newly identified ISs in the Sulfolobus tokodaii and Sulfolobus solfataricus genomes. The structures of these fragmentary copies suggest several distinct molecular mechanisms which, in the absence of precise excision, inactivate ISs and gradually eliminate the defective copies from Sulfolobus genomes.     

Anvaya: a workflows environment for automated genome analysis

Anvaya is a workflow environment for automated genome analysis that provides an interface for several bioinformatics tools and databases, loosely coupled together in a coordinated system, enabling the execution of a set of analyses tools in series or in parallel. It is a client-server workflow environment that has an advantage over existing software as it enables extensive pre & post processing of biological data in an efficient manner. "Anvaya" offers the user, novel functionalities to carry out exhaustive comparative analysis via "custom tools," which are tools with new functionality not available in standard tools, and "built-in PERL parsers," which automate data-flow between tools that hitherto, required manual intervention. It also provides a set of 11 pre-defined workflows for frequently used pipelines in genome annotation and comparative genomics ranging from EST assembly and annotation to phylogenetic reconstruction and microarray analysis. It provides a platform that serves as a single-stop solution for biologists to carry out hassle-free and comprehensive analysis, without being bothered about the nuances involved in tool installation, command line parameters, format conversions required to connect tools and manage/process multiple data sets at a single instance.     

Identification and distribution of insertion sequences of Paracoccus solventivorans

Three novel insertion sequences (ISs) (ISPso1, ISPso2, and ISPso3) of the soil bacterium Paracoccus solventivorans DSM 11592 were identified by transposition into entrapment vector pMEC1. ISPso1 (1,400 bp) carries one large open reading frame (ORF) encoding a putative basic protein (with a DDE motif conserved among transposases [Tnps] of elements belonging to the IS256 family) with the highest levels of similarity with the hypothetical Tnps of Rhodospirillum rubrum and Sphingopyxis macrogoltabida. ISPso2 (832 bp) appeared to be closely related to ISPpa2 of Paracoccus pantotrophus DSM 11072 and IS1248 of Paracoccus denitrificans PdX22, both of which belong to the IS427 group (IS5 family). These elements contain two overlapping ORFs and a putative frameshift motif (AAAAG) responsible for production of a putative transframe Tnp. ISPso3 (1,286 bp) contains a single ORF, whose putative product showed homology with Tnps of ISs classified as members of a distinct subgroup of the IS5 group of the IS5 family. The highest levels of similarity were observed with ISSsp126 of Sphingomonas sp. and IS1169 of Bacteroides fragilis. Analysis of the distribution of ISs of P. solventivorans revealed that ISPso2-like elements are the most widely spread of the elements in nine species of the genus PARACOCCUS: ISPso1 and ISPso3 are present in only a few paracoccal strains, which suggests that they were acquired by lateral transfer. Phylogenetic analysis of Tnps of the novel ISs and their closest relatives showed their evolutionary relationships and possible directions of lateral transfer between various bacterial hosts.