Within-Genome Evolution of REPINs: a New Family of Miniature Mobile DNA in Bacteria
[['Short repetitive sequences are a feature of most genomes and have consequences for genome function and evolution [1] , [2] .', 'Often attributable to the proliferation of selfish elements [3] , [4] , short repeats also arise from amplification processes , such as replication slippage [5] and via selection on genome architecture [6]–[8] .', 'Repetitive DNA in bacterial genomes is less prominent than in eukaryotes , nonetheless , an over abundance of short oligomers is a hallmark of almost every microbial genome [9] .', 'Known generically as interspersed repetitive sequences , these elements have a history of exploitation as signatures of genetic diversity ( e . g . , [10]–[12] ) , but their evolution , maintenance and mechanism of within- and between-genome dissemination are poorly understood [9] , [13]–[16] .', 'Interspersed repetitive sequences fall into several broad groups each sharing short length ( individual units range from ∼20 to ∼130 bp ) , extragenic placement , and palindromic structure [9] , [17] .', 'REPs ( repetitive extragenic palindromic sequences ) – also known as PUs ( palindromic units ) – range from ∼20 to ∼60 bp in length , possess an imperfect palindromic core , are widespread among bacteria , and occur hundreds of times per genome [13] , [18]–[23] .', 'While often existing as singlets , REPs also form a range of complex higher order structures termed BIMEs ( bacterial interspersed mosaic elements ) [14] .', 'CRISPRs ( clustered regularly interspaced short palindromic repeats ) are a further , higher order composite of REP-like sequences that are formed from direct repeats of short ( ∼30 bp ) palindromic sequences interspersed by similar size unique non-repeated DNA ( [24]; reviewed in [25] ) .', 'Recent work shows that the unique sequences are often phage derived and that CRISPRs , along with associated proteins , confer resistance to phage by targeting viral DNA [25] , [26] .', 'Non-autonomous DNA transposons form a more distinct family of repetitive sequences defined by their size ( ∼100 to ∼400 bp ) and presence of terminal inverted repeats .', 'Also known generically as MITEs ( miniature inverted repeat transposable elements ) , non-autonomous transposons depend on transposase activity encoded by co-existing autonomous transposons for dissemination [4] .', 'Identified initially in plants [27] , where evidence of active transposition has been obtained [28] , recent bioinformatic analyses suggest that they also occur in bacteria [29] , [30] .', 'For example , ERICs ( enterobacterial repetitive intergenic consensus ) – found in a range of enteric bacteria including Escherichia coli , Salmonella and Yersinia [31] – and NEMISs ( Neisseria miniature insertion sequences ) in pathogenic neisseriae [32] are thought to be non-autonomous transposons ( MITEs ) .', 'Scenarios for the origins and functional significance of non-autonomous elements , and to a lesser extent CRISPRs , can be envisaged , but this is not so for the majority of short interspersed repetitive sequences .', 'Nonetheless , studies of specific elements in particular genetic contexts have uncovered evidence of functional roles ranging from transcription termination and control of mRNA stability , to binding sites for DNA polymerase I ( reviewed in [9] ) .', 'However , the fact that the distribution and abundance of elements show substantial among-strain diversity [16] , [22] suggests that the range of functional roles is incidental , arising from , for example , co-option or genetic accommodation [31] .', 'Differences in the distribution and abundance of repetitive elements among closely related strains carries additional significance in that it suggests that the evolution of these elements is independent of the core genome .', 'This is particularly apparent from comparisons of closely related strains .', 'For example , Pseudomonas fluorescens isolates SBW25 and Pf0-1 are closely related and yet highly dissimilar in terms of the nature , abundance and distribution of interspersed repetitive elements [22] , even , as we show here , at the level of REPs . While this may reflect unequal rates of element loss , an alternative possibility is independent acquisition .', 'Implicit in this suggestion is the notion that repetitive elements are genetic parasites [13] , [31] , [33] .', 'The idea that REPs are selfish elements is not new [13] , [31] , [33]; however , there is little evidence – either direct or indirect – to support such an assertion .', 'Indeed , the small size of REPs makes a mechanism for autonomous replication difficult to envision , however , the recent discovery of a proximal association between REPs and IS200-like elements , termed RAYTs ( REP-associated tyrosine transposases ) [23] , raises interesting possibilities and suggests shared ancestry between RAYTs and certain REP families .', 'Evolutionary approaches to the analysis of sequence motifs can be highly informative [34] .', 'While there is a ready tendency to assume that motifs recognized by search algorithms have functional significance , this need not be so .', 'Neutral evolutionary processes alone ( nothing more than random chance ) ensure that short sequences will occur multiple times within any given genome .', 'Thus , before concluding functional significance , it is necessary to test the null hypothesis of chance .', 'Should this hypothesis be rejected , then the conclusion that over-abundance of short sequences is attributable – at least in part – to natural selection is sound .', 'Moreover , evidence for selection justifies the assumption of functional significance .', 'A key issue , however , is the level of biological organization at which functionality has been selected .', 'There are two distinct possibilities: short repeats may have evolved because of selective benefits conferred on the cell , but alternatively , they may deliver benefits at the level of the gene – more specifically , at the level of a genetic element , of which the repeat sequence is a component .', 'Distinguishing between these two alternatives is possible , although not necessarily straightforward .', 'Indeed , whereas on initial emergence , selection is likely to operate exclusively at one level , over time , it is likely to shift to encompass multiple levels [4] , [16] .', 'Here , we take a fresh and unbiased look at bacterial genome sequences in order to analyze the frequency and nature of short sequence repeats .', 'Our approach is informed by evolutionary theory and begins free of assumptions regarding functional significance .', 'Accordingly , the null hypothesis that short sequence repeats are no more frequent than expected by chance is the initial focus .', 'We begin by interrogating the P . fluorescens SBW25 genome .', 'Using suitable null models we show that over-abundant oligomers – which cannot be accounted for by chance alone – fall into three separate groups , each with characteristics typical of REPs . Highly significant differences in patterns of REP abundance and diversity between SBW25 and a second closely related P . fluorescens strain led us to question the hypothesis that the causes of REP diversity are linked to cellular function .', 'This prompted a search for a replicative unit , which , based on patterns of REP distribution , we argue is a REP doublet .', 'We refer to these entities as REPINs ( REP doublets forming hairpins ) and provide evidence from population sequencing for the existence of a putative transposition intermediate .', 'Finally , extension to a range of RAYT-containing bacterial genomes including E . coli K-12 and Nostoc punctiforme indicate that REP sequences , organized as REPINs , define a class of hitherto unrecognized miniature non-autonomous mobile DNA .'], ['Defining repetitive DNA on the basis of short sequences ranging from 10–20 nucleotides is simple and can be done logically without invoking heuristics and approximations ( for longer sequences exact repetitions are rare ) .', 'Figure 1 shows that the P . fluorescens SBW25 genome harbors numerous repetitive sequences: the most common 10-mer occurs 832 times; the most common 20-mer occurs 427 times .', 'While these numbers appear significant , it is possible that they are no more than expected by random chance .', 'To test this hypothesis , 100 random genomes were generated , with the same dinucleotide content , replication bias and length , as the SBW25 genome .', 'The frequency of the most abundant oligonucleotides was determined from both leading and lagging strands .', 'Figure 1 shows that the most abundant 10-mer from the randomly generated genomes occurs 304 times .', 'For longer sequence lengths this number rapidly decreases ( four instances in the case of 20-mers ) : the number of repeats expected by chance alone is thus much lower than observed .', 'In total , there are 108 different 10-mers and 14 , 351 different 20-mers that occur significantly more often in the P . fluorescens genome than the most abundant oligonucleotides from randomly generated genomes ( P<0 . 01 , Figure S1 ) .', 'While compelling evidence for the existence of over-representation of short sequences , gene duplications could in part account for these findings [35] .', 'We therefore sought an alternative null model .', 'P . fluorescens Pf0-1 , one of the closest relatives of SBW25 , shares the same GC-content and has a highly similar dinucleotide content ( Table S1 ) ; coding density differs by 1 . 7% and the genome length differs by 4% ( 6 , 722 , 539 bp for SBW25 and 6 , 438 , 405 bp for Pf0-1 , [22] ) .', 'The close similarity means that any bias in the representation of short sequences due to duplicative evolutionary processes , or other selective mechanisms , should be similar in both genomes .', 'As in SBW25 , over-represented short sequences in Pf0-1 are more frequent than expected by chance ( Figure 1 ) , however , a considerable difference in short sequence frequency is apparent .', 'The difference between SBW25 and Pf0-1 is greatest at a sequence length of 16 , where the most abundant sequence in SBW25 occurs 618 times – over 11 times more frequently than the most abundant 16-mer in Pf0-1 ( Figure S2 ) .', 'On the basis of comparisons to both the random null model and the Pf0-1 genome we deemed all SBW25 16-mers occurring more than 55 times ( the frequency of the most abundant 16-mer in Pf0-1 ) to be over-represented .', 'This led us to reject the null hypothesis that chance alone explains the occurrence of short repetitive sequences in the SBW25 genome .', 'Accordingly , we attribute over-representation of oligonucleotides to selective processes .', 'The collection of over-represented 16-mers together encompasses 96 different sequences; however , a cursory glance suggested that many share similarity .', 'Using a grouping method designed to detect overlapping subsets of sequences ( Methods and Figure S3 ) , the 96 sequences were found to be members of just three separate sequence groups ( GI , GII and GIII ( Figure S4 ) ) , each containing an imperfect palindrome ( the palindrome overlaps the most abundant 16-mer in GI and GII , but is part of the most abundant 16-mer in GIII ( Table 1 ) ) .', 'The most abundant 16-mers of each group together occur 1 , 067 times .', 'The majority of these sequences are extragenic; only 14 16-mers overlap with genes .', 'Together these data show that the three groups of 16-mers are over-represented in the SBW25 genome , contain an imperfect palindromic core and are primarily extragenic .', 'Possessing the hallmarks of repetitive extragenic palindromic ( REP ) sequences , we conclude that the three groups of 16-mers are , for all intents and purposes , synonymous with REPs . In order to accommodate the possibility of related family members , we generated a pool of sequences that differed to GI , GII and GIII sequences by up to four bases .', 'This generated 488 , 373 different 16-mers of which 1 , 861 were located in extragenic space .', 'To define the proportion of false positives the search was repeated by interrogating randomly generated extragenic space ( with the same dinucleotide content and length of each individual extragenic space ) for matches to the 488 , 373 different 16-mers .', 'This showed that 12% of all sequences with up to four substitutions are false positives ( sequences unrelated to GI , GII or GIII ) .', 'Repeating the analysis with the subset of sequences , which differ firstly by three and subsequently , two substitutions showed that 2% and 0 . 2% of matches are false positive , respectively .', 'For two substitutions the false positive rate is low enough to conclude that the described repetitive sequence families consist of at least 1 , 422 members ( Table 2 ) .', 'The precise number of members belonging to each of the GI , GII and GIII groups cannot be determined because with a degeneracy of two , some sequences fall into more than one group .', 'The selective causes for the prevalence of GI , GII and GIII sequences in the SBW25 genome are of considerable interest .', 'Although implicit in many studies is the notion that REP-like sequences have evolved because of their selective benefit to the cell ( as transcription binding sites , termination signals and the like [20] , [36] , [37] ) , it is also possible that selection has favored their evolution as a consequence of benefits delivered to a genetic ( parasitic ) element , of which the repeat sequence is a component .', 'The highly significant differences in the frequency , nature and genomic location of short repetitive sequences in SBW25 , compared to Pf0-1 make a compelling case for the latter .', 'If the prevalence of GI , GII and GIII sequences is a consequence of gene-level selection , then this implies the existence of a replicative entity – a genetic element that has the capacity to reproduce within the genome .', 'The distribution of REP sequences is likely to provide some information .', 'One way to quantify the distribution is to measure distances between neighboring REP sequences and compare these to distances between REPs generated by a null ( random ) model .', 'If individual REPs are randomly distributed then this would suggest the individual REP as replicative unit .', 'If the distance between adjacent REPs is non-random , then this may suggest the evolving entity is some higher order arrangement of REPs . To construct the null model , 1 , 053 ( the number of invariant GI , GII and GIII sequences in extragenic space ) non-overlapping 16 bp segments were positioned at random within the extragenic space of the SBW25 genome .', 'This process was repeated 10 , 000 times and the average occurrence of the distance between neighboring elements calculated .', 'Equivalent data for the 1 , 053 over-represented REPs is shown in Figure 2 .', 'A comparison between the two histograms reveals marked differences in the distributions of distances between next-neighbors .', 'Most striking is the strong bias toward specific inter-element distances .', 'This marked skew shows that REPs are not independently distributed and is suggestive of an underlying copying mechanism involving at least two REP sequences .', 'Of note is the fact that doublets typically comprise pairs of identical GI , GII or GIII sequences and are rarely mixed ( although some exceptions are discussed below ) ( Figure 2 ) .', 'To explore the possibility that the replicative unit is an entity comprised of two REP elements ( a REP doublet ) we determined the number of singlets , doublets , triplets and higher order arrangements of REPs ( REP clusters ) by examining the 400 bp flanking either side of each REP for the presence of REP sequences ( Figure S5 ) .', 'Once again , the results of this analysis were compared to the null ( random ) model used above .', 'According to the random model , 58% of all REP sequences are expected to occur as singlets , whereas data from SBW25 shows that just 18% are singlets .', 'In contrast , 61% of all REPs are organized as doublets , which is significantly greater than the 17% expected by chance ( Table 3 ) .', 'Interestingly , REP triplets are rarer than expected , whereas several higher order arrangements of REPs , including two sets of twelve ( see below ) , are more frequent than expected ( Table 3 ) .', 'The highly significant over-representation of REP doublets suggests that the doublet defines an appropriate replicative unit .', 'If true , then the distribution of doublets across extragenic space should be unaffected by neighboring REP elements and should thus conform approximately to a null ( random ) model .', 'To test this hypothesis , random distributions of REP doublets over extragenic space were compared to actual REP clusters found in SBW25 ( Table 4 ) .', 'However , because the distance between REPs ( in the doublet conformation ) varies ( Figure 2 ) , two random models were generated based on the two most common inter-REP spacings: 71 bp ( a doublet of GI REPs ) and 110 bp ( a doublet of GII REPs ) .', 'Simulations were based on the random assignment of 560 REP doublets ( corresponding to the sum of REP clusters ( of two or more ) in Table 3 ) to extragenic space and were repeated 10 , 000 times .', 'Although the two segments differ significantly in size , simulations for each family gave remarkably similar results ( Table 4 ) .', 'Together these data show that the observed number resembles that predicted if the doublets are randomly distributed .', 'A further prediction concerns evolutionary processes affecting doublets vs . singlets .', 'If REP doublets are the replicative unit , then singlets are likely to derive from doublets , either by decay ( divergence ) of the neighboring element , or by destruction of the doublet through insertion or deletion .', 'In either case the REP singlet is expected to be non-functional ( immobile ) and thus subject to random genetic drift .', 'REP doublets on the other hand – being ( according to our hypothesis ) functional and potentially mobile – are expected to be shaped by selection: genetic diversity of REP singlets should thus be greater than doublets .', 'To test this hypothesis we extracted GI , GII and GIII sequences from the SBW25 genome plus all related sequences that varied by up to two positions .', 'Since only two nucleotide differences distinguish GII and GIII sequences from a GI sequence , GII and GIII sequences were defined by two fixed ( invariant ) positions ( GII: 2T , 6C; GIII: 6A , 13T ) .', 'After extraction , sequences from each group were divided into a set of 16-mers obtained from singlets , a set of 16-mers from doublets and a set of 16-mers obtained from clusters ( where a cluster contains three or more REPs ) .', 'For all nine sequence groups ( three from each GI , GII and GIII group ) the pairwise identity was calculated ( Figure 3 , see Methods for details ) .', 'The average pairwise identity of 16-mers obtained from REP doublets is significantly greater than the average pairwise identity of 16-mers obtained from REP singlets: this is true for comparisons within each of the REP groups ( P<1e-10 for GI; P<1e-8 for GII and GIII ) .', 'Analysis of the organization of REP doublets shows that in the majority of cases , pairs of REPs ( 93% of all 430 REP doublets ) – of either the GI , GII , or GIII types – are organized as two inverted REP sequences that overlap the most abundant 16-mer ( Figure 4A and 4B ) .', 'While the spacer region between REPs shows less conservation than evident in the REPs themselves , secondary structure predictions for ssDNA shows that the conserved bases on each side pair resulting in a hairpin ( Figure 4E ) .', 'Thus , while selection appears to favor highly conserved nucleotide arrangements for REP and adjacent sequences , the critical features of the intervening sequence would appear to be length , and capacity to form a hairpin .', 'Indeed , compensatory changes on either side of the predicted hairpin are common ( Figure 4A ) .', 'Finally , if our assertion that the doublet defines a replicative entity is correct , then evidence of movement could in principle come from population sequencing .', 'To this end we interrogated 55 , 768 , 706 paired-end Illumina reads ( 36–76 bp long ) obtained from sequencing DNA extracted from 5×10e9 SBW25 cells , for evidence of insertion and excision events .', 'A total of 18 putative insertions were detected , however , the possibility of false positives could not be discounted .', 'A similar search for excision events proved more profitable: three single reads were identified which mapped to three different locations on the genome , each corresponding to unique sequences flanking a GI REP doublet ( Figure 4C and Figure S6 ) .', 'However , the expected doublet was absent from all sequence reads leading us to conclude that these sequences were from DNA molecules from which the doublet had excised .', 'Additionally , we observed 200 individual sequence reads spanning a GII REP doublet indicating its excision from the entire population ( Figure S6 ) .', 'That these events could result from machine and/or chemistry error is improbably low .', 'Furthermore , a search for evidence of REP singlet deletions from the ∼56 million Illumina reads failed to find evidence of a single such event ( see Methods ) .', 'Details of the three excised GI doublets are shown in Figure 4C and 4D .', 'Of particular interest is the asymmetrical nature of the deleted sequence: in all instances it begins ( in the left-hand ( 5′ ) end ( Figure 4B ) ) at the start of the invariant sequence defined by the most conserved 16-mer and extends through the spacer region into the second REP sequence .', 'However , rather than finish at the end of the conserved 16-mer , the deletion truncates at the 3′-end of the right-hand REP sequence , leaving the last ∼6 bp of invariant sequence intact ( Figure 4C ) .', 'Secondary structure predictions show a hairpin structure with a 5′-single strand tail .', 'Although the structures of the hairpins are not identical ( due to differences in the sequence of the space region ) the 5′-tail is a feature of the excised entity in all instances ( Figure 4E ) .', 'It is possible that the excised sequences define a putative transposition intermediate .', 'Together the above analyses implicate REP doublets as a unit of selection: a family of mobile DNA that has , until now , eluded recognition .', 'Although REP doublets have previously been noted as one of many different higher order arrangements of REPs , they have not before been implicated as replicative entities [16]–[20] .', 'Furthermore , in previous discussions of higher order arrangements it has been assumed that the singlet is the basic building block .', 'In contrast , our data supports the view that REP singlets are defunct remnants of once functional REPINs .', 'Because of their likely evolutionary relevance , a label that defines the replicative entity appears warranted .', 'Henceforth we refer to REP doublets forming hairpins as REPINs .', 'While the majority of REPINs exist as singlets , some higher order arrangements are apparent ( above and Table 4 ) .', 'These are of two main types: those showing a distinctive ordering and those with no apparent structure .', 'REPINs occurring in ordered clusters are typically arranged as tandem repeats of nearly identical REPINs – including the flanking sequences ( Figure S7 ) .', 'With 16 such clusters distributed throughout the genome , these arrays are the most common higher order arrangement of REPINs in SBW25 .', 'The largest cluster consists of four REPINs ( plus an additional REP sequence ) with a total length of over 700 bp .', 'Three higher order REPIN clusters are of particular note: one from each of the three distinctive REPIN groups ( GI , GII and GIII ) each located adjacent to one of the three recently identified REP-associate tyrosine transposases ( RAYTs , [23] ) ( pflu3939 , pflu4255 and pflu2165 ) .', 'The fact that a different REPIN cluster is located beside each of the RAYTs , combined with the fact that REPINs ( and REPs ) in SBW25 come in three distinct flavors , raises the possibility that RAYTs are intimately linked to REPIN mobilization ( Figure 5 ) .', 'REPINs in clusters lacking obvious organization are found in five regions of the genome and typically consist of two unrelated REPINs .', 'Close inspection suggests that these clusters are formed by insertion of REPINs into , or next to , existing REPINs .', 'REPs also form higher order arrangements .', 'These are of two distinct types: the first involves highly organized tandem arrays of GI and GIII REP sequences: GI REPs are separated from GIII REPs by 112 bp; GIII REPs are separated from GI REPs by 72 bp .', 'Five such tandem arrays are located at ∼2 Mbp all of which are found in forward orientation , six are found ∼4 Mbp in reverse orientation ( at a distance of ∼2 Mbp from the origin of replication ) .', 'The two largest tandem arrays both contain 12 GI and GIII sequences , one found at ∼4 . 1 Mbp the other at ∼2 . 5 Mbp ( Figure S8 ) .', 'These two arrays are almost identical copies of each other , but found in opposite orientations on opposite sides of the genome .', 'The second type of tandemly organized REP sequences consists solely of evenly spaced GI sequences found at two positions in the genome .', 'Similar to the GI–GIII tandem arrays one GI tandem array is found in forward and the other one in reverse orientation .', 'REPIN dissemination could occur via the exploitation of a functional transposase encoded separately within the genome .', 'Non-autonomous DNA transposons ( MITEs ) do precisely this and typically consist of two inverted repeats .', 'REPINs also consist of two inverted repeats ( REP sequences ) and , as mentioned above , may exploit the putative transposase encoded by RAYTs .', 'If REP sequences in other genomes are components of REPINs – and disseminate via RAYT-encoded transposase activity – then , given the broad distributions of RAYTs [23] , REPINs are likely to be a common feature of bacterial genomes; they are also likely to share common ancestry .', 'Although a fully comprehensive among-genome analysis is beyond the scope of this paper we nonetheless analyzed REP sequence clusters in a variety of genomes containing RAYTs .', 'To this end REP sequences were selected from 18 different bacterial strains including all fully sequenced Pseudomonas genomes , the genomes of E . coli K-12 DH10B and Salmonella enterica serovar Paratyphi A AKU 12601 ( chosen because of their significance for REP research ) and the genomes of Thioalkalivibrio HL-EbGR7 and N . punctiforme PCC73102 ( chosen because of their distant relation to Pseudomonas ) .', 'A phylogenetic analysis of the RAYTs was firstly undertaken ( Figure S9 ) .', 'Notably , RAYTs from these strains form two distinct evolutionary lineages with evidence of multiple independent introductions .', 'For example , the genus Pseudomonas is separated into two sets of species defined by the presence of either ‘clade I’ or ‘clade II’ RAYTs .', 'The genome of Thioalkalivibrio contains one clade I and one clade II RAYT .', 'Several other genomes , in addition to SBW25 , contain more than a single RAYT , but these almost never cluster .', 'In fact the most closely related RAYTs are found in different genomes .', 'Overall the distribution of RAYTs among distantly related organisms shows evidence of lateral gene transfer; however , at the species level , lateral gene transfer does not seem to occur frequently as evident by the fact that RAYT phylogeny is largely congruent with the relationship among species ( Figure S9 ) .', 'Since REP sequences have been shown to be associated with RAYT genes ( this work and [23] ) , we interrogated non-coding DNA flanking each RAYT for 16-mers that were repetitive , extragenic and palindromic , that is , are REPs . In each instance a REP was identified ( Table S2 ) .', 'To test the hypothesis that REPs are organized as REPINs an analysis of the distribution of REPs was performed on each genome as described above ( also see Methods ) and included all REP sequences that differed from the consensus by up to two nucleotides .', 'Results were expressed as the ratio of REP singlets to doublets , where ratios greater than two indicate that REPs occur predominantly as singlets .', 'Ratios less than two mean that REPs occur predominantly as doublets .', 'Figure 6 shows a histogram of singlet to doublet ratios for REP sequences associated with clade I RAYTs .', 'Of the 20 REP sequence classes ( one associated with each RAYT , some genomes contain more than one RAYT e . g . , SBW25 ) 17 had singlet to doublet ratios of less than two , indicating that most REPs occur as doublets .', 'The majority of doublets contained REPs as inverted pairs ( Table S3 ) as expected of REPINs .', 'Our simple search method did not return conclusive results for clade II REP sequences .', 'One possibility is that the REPIN structure in these genomes is less conserved .', 'To this end we performed a secondary structure prediction on a sample of REP sequences .', 'In all instances we found the general REPIN composition to hold ( two inverted REP sequences separated by a short stretch of DNA and forming a hairpin , Figure S10 ) , with the exception of REP sequences found in P . stutzeri: interestingly no REPINs were identified in this genome .', 'We also analyzed higher order arrangements for clade I REP sequences , but these were not present in all analyzed genomes .', 'They were predominantly found in P . syringae and P . fluorescens , although two REP sequence classes were also detected in P . putida ( Table S3 ) .', 'No correlation was found between the singlet to doublet ratio and cluster formation .', 'Taken together , the systematic cluster analysis of clade I REP sequences and secondary structure prediction of a selection of clade II REP sequences suggest that the organization of REP sequences into REPINs is a necessary condition for REP sequence distribution .'], ['Short interspersed repetitive sequences are widely distributed in bacteria , but past studies have shed little light on their evolutionary origins .', 'We began by examining the abundance and distribution of short sequences in P . fluorescens SBW25 and showed , by comparison against a random ( null ) model , and subsequently against Pf0-1 , that short sequences are over-represented .', 'Moreover , we found that short repetitive sequences fall into three distinct groups ( GI , GII and GIII ) , each bearing characteristics typical of REP sequences , that is , they are repetitive , extragenic and palindromic .', 'In order to discount the possibility that REP sequences are the product of mutation pressure ( a possibility already called into doubt by comparison to the random model ) we took advantage of the closely related Pf0-1 genome .', 'Comparisons using this null model – based upon a genome likely to have been shaped by similar underlying evolutionary processes – allowed us to emphatically reject the possibility that REP evolution can be explained by drift .', 'Our data thus indicate natural selection as the primary driver of REP sequence evolution .', 'A critical issue is the nature of the entity upon which selection acts .', 'Evidence that this entity comprises a doublet of REP sequences – a REP doublet forming a hairpin structure ( REPIN ) – came firstly from analysis of the distribution of REPs in extragenic space .', 'The striking departure from a random model shown in Figure 2 , along with clear bias toward specific distances between REPs , pointed to the REPIN as the replicative entity .', 'The hypothesis was further tested by examining the distribution of REP doublets in extragenic space , by measuring nucleotide diversity in singlets versus doublets , and by analysis of the conserved features of REPINs .', 'Finally , the existence of REPINs as actively mobile entities was bolstered through the discovery of four deletion events that may define putative transposition intermediates ( Figure 4 ) .', 'A previous analysis of the SBW25 genome using various repetitive DNA finding algorithms [22] revealed numerous repeat families .', 'Two of these , the so named R0 and R2 repeats have characteristics similar to REPINs; indeed , a comparison ( Table S4 ) shows a correspondence between REPINs and the R0 and R2 repeats .', 'In general R0 repeats map to GI REPINs , while R2 repeats correspond to a mixture of both GII and GIII REPINs .', 'The mechanism by which REPINs are disseminated is a central , but unresolved issue .', 'Recently , a hypothesis concerning REP seq
['Introduction', 'Results', 'Discussion', 'Methods']
['Repetitive sequences are a conserved feature of many bacterial genomes .', 'While first reported almost thirty years ago , and frequently exploited for genotyping purposes , little is known about their origin , maintenance , or processes affecting the dynamics of within-genome evolution .', 'Here , beginning with analysis of the diversity and abundance of short oligonucleotide sequences in the genome of Pseudomonas fluorescens SBW25 , we show that over-represented short sequences define three distinct groups ( GI , GII , and GIII ) of repetitive extragenic palindromic ( REP ) sequences .', 'Patterns of REP distribution suggest that closely linked REP sequences form a functional replicative unit: REP doublets are over-represented , randomly distributed in extragenic space , and more highly conserved than singlets .', 'In addition , doublets are organized as inverted repeats , which together with intervening spacer sequences are predicted to form hairpin structures in ssDNA or mRNA .', 'We refer to these newly defined entities as REPINs ( REP doublets forming hairpins ) and identify short reads from population sequencing that reveal putative transposition intermediates .', 'The proximal relationship between GI , GII , and GIII REPINs and specific REP-associated tyrosine transposases ( RAYTs ) , combined with features of the putative transposition intermediate , suggests a mechanism for within-genome dissemination .', 'Analysis of the distribution of REPs in a range of RAYT–containing bacterial genomes , including Escherichia coli K-12 and Nostoc punctiforme , show that REPINs are a widely distributed , but hitherto unrecognized , family of miniature non-autonomous mobile DNA .']
['DNA sequences that copy themselves throughout genomes , and make no specific contribution to reproductive success , are by definition “selfish . ”', 'Such DNA is a feature of the genomes of all organisms and evident by virtue of its repetitive nature .', 'In bacteria the predominant repetitive sequences are short ( ∼20 bp ) , extragenic , and palindromic .', 'These so-called REP sequences may occur many hundreds of times per genome , but their origins and means of dissemination have been a longstanding mystery .', 'We show that REPs are components of higher-order replicative entities termed REPINs , which are themselves thought to be derived from REP sequences that flanked an ancestral autonomous selfish element .', 'In this ancestral state the REP sequences were likely to have been critical for the movement of the selfish element , but were devoid of any capacity to replicate independently .', 'REPINs , on the other hand , have evolved to have a life of their own , albeit one that exploits—even enslaves—a genetic element upon which their existence depends .', 'REPINs are the ultimate non-autonomous , super-streamlined , selfish element and are widespread among bacteria .']
['genome evolution', 'genome sequencing', 'coevolution', 'genome complexity', 'comparative genomics', 'biology', 'evolutionary genetics', 'genetic drift', 'adaptation', 'natural selection', 'genetics', 'genomics', 'evolutionary biology', 'evolutionary processes', 'genetics and genomics']