{'fraction_non_alphanumeric': 0.02994144906623748, 'fraction_numerical': 0.01631656353192552, 'mean_word_length': 5.133083730428863, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 1, 'https://': 0, 'lorem ipsum': 0, 'www.': 1, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 5, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Characterization of Species-Specific Repeats in 613 Prokaryotic Species

Prokaryotes are in general believed to possess small, compactly organized genomes, with repetitive sequences forming only a small part of them. Nonetheless, many prokaryotic genomes in fact contain species-specific repeats (>85 bp long genomic sequences with less than 60% identity to other species) as we have previously demonstrated. However, it is not known at present how frequent such species-specific repeats are and what their functional roles in bacterial genomes may be. Therefore, we have conducted a comprehensive survey of prokaryotic species-specific repeats and characterized them to examine as to whether there are functional classes among different repeats or not and how they are mutually related to each other. Of the 613 distinct prokaryotic species analyzed, 97% were found to contain at least one species-specific repeats. It seems interesting to note that the species-specific repeats thus identified appear to be functionally variable in different genomes: in some genomes, they are mostly associated with duplicated protein-coding genes, whereas in some other genomes with rRNA and tRNA genes. Contrary to what may be expected, only one-fourth of the species-specific repeats were found to be associated with mobile genetic elements.


Introduction
To date, more than 1000 prokaryotic species have been completely sequenced. These genome data provide an opportunity to conduct large-scale genome-wide studies in comparative genomics. Many analyses are carried out to explore eukaryotic genomes side by side, but not many comparative analyses are performed with prokaryotes. Even though prokaryotic organisms demonstrate recognizable, common architectural principles, they have a great variability in ecology and in metabolic and genomic complexity. 1 For example, the lengths of sequenced bacterial genomes range from 180 kb (Carsonella ruddii) 2 to 13 Mb (Sorangium cellulosum). 3 Although a lot of this variability can be attributed to varying metabolic complexity, it also suggests that some species have more compact genomes than others.
Large differences exist in the fraction of intergenic regions. 1 The percentage of intergenic regions in the genomic sequence varies from 5% (Thermotoga neapolitana, NC_011978) to 50% (Sodalis glossinidius, NC_007712). The median value of the fraction of intergenic regions per genome size is 12%.
Though prokaryotic genomes are known to be compactly organized, there is still room for different repetitive sequences in them. 4,5 There are very few comprehensive studies of repetitive sequences; most of the studies focus on specific types of repeats or on a limited number of bacterial species. 4 -8 The processes generating repeats include duplication, horizontal gene transfer, transposition, and replicon fusion. 1,5,9 -17 These processes allow bacteria to adapt to environmental changes or to evolve into pathogens (owing to the distribution of genes encoding toxins, effector proteins, cell wall modification enzymes, fitness factors, and antibiotic and heavy metal resistance determinants). 10 A large fraction of known types of repetitive sequences belong to integrative and conjugative elements, also known as mobilome. 18,9 The mobilome consists of bacteriophages, plasmids, and transposable elements. It mediates the movement of DNA within and between genomes. Integrative elements play a key role in the emergence of infectious diseases, antibiotic resistance and other problems. 1,19 -21 Thus, different mechanisms give rise to repeats and also different roles are suggested for repeats. Though many papers about completely sequenced prokaryotic genomes also try to analyze repetitive sequences, the characterization is frequently limited to some types of repeats (e.g. CRISPR, MITE, IS) or repeats as whole without classification. 22,23 No exhaustive research about the characterization of prokaryotic speciesspecific repeats is published. However, as further shown, different types of species-specific repeats with different functions are present in prokaryotic genomes. At least some of them are the characteristic of a particular species (e.g. specific protein-biosynthesis-coding genes or specific bacterial toxins).
We have previously reported the existence of species-specific repeats in several bacterial species of medical importance and we have developed a computational method for de novo detection of such repeats previously. 24 However, two questions remain unanswered: do all species have species-specific repeats and can we associate these repeats with sequences other than mobile genetic elements?
The first topic we cover in this paper is the general description of prokaryotic repeats. It is important for species characterization and identification to know what fraction of a prokaryotic genome contains repetitive DNA regions, how many of them are species-specific, and what is the common size range of the repetitive sequences. Our main focus is on species-specific repeats, with a few comparisons with universal repeats (definitions of these repeat types are given below, at the beginning of the 'Results' section).
Thereafter, we focus on the functional characterization of species-specific repeats. The main interest here is whether species-specific repeats are reliable for use as target sequences in identifying a species or whether they are really only selfish elements or sequences with unknown function according to the common view. For this, different functional classes of repeats are distinguished. In this paper we use the following classes to describe repeats: (I) mobile genetic elements, (II) repeats associated with RNA genes, (III) repeats associated with protein-coding genes, and (IV) non-coding repeats.
Finally, we discuss the different practical aspects of detecting prokaryotic repetitive sequences. The Supplementary material contains a list of all of the detected sequence-specific repeats.

Genomic sequences used
To identify prokaryotic repeats, DNA sequences for all completely sequenced prokaryotic chromosomes and plasmids were retrieved from the FTP site of the NCBI RefSeq (ftp://ftp.ncbi.nih.gov/genomes/Bacteria). The human genome (version NCBI36) was downloaded from Ensembl (ftp://ensembl.org/pub/). A full list of analyzed genomes and their basic features (genome length, superkingdom, number of chromosomes, number of plasmids, etc.) is provided in the Supplementary Table S3. In total, 876 chromosomal genomes (63 archaeal and 819 eubacterial) and 613 prokaryotic species (containing 54 archaeal and 561 eubacterial genomes) were analyzed. Of the 613 species, 95 are represented with more than one strain of chromosomal genome sequence. For 232 prokaryotic species at least one plasmid genome was available, 733 prokaryotic plasmid genomes were downloaded in total. The RefSeq accession numbers of the virus genomes used in the current work are provided in the additional files (Additional File 1: Supplementary  Table S3). To find strain-specific repeats online, the NCBI BLAST database for microbes was used on genomic sequences with sequencing status 'Complete, in Progress or Assembly'.

Annotation data
To explore the positions of prokaryotic repeats in relation to genes, gene coordinates were retrieved from the FTP site of the NCBI RefSeq (ftp://ftp.ncbi.nih.gov/ genomes/Bacteria) for all analyzed prokaryotic chromosomal genomes.

Method for identifying microbial repetitive sequences
The method for finding microbial repeats has been described earlier. 24 The availability of complete genome sequences for the species of interest is required in order to find repetitive regions. The species of interest is called the 'target genome' in this article. In the context of species-specific repetitive sequences, the 'non-target genome' is defined as genomic DNA sequence(s) of all other organisms not in the target genome group. In the case of bacterial species-specific repeats, all sequenced chromosomal genomes of prokaryotic species and their plasmid genomes were used as non-target genomes. In this study, the human genome was always used by default as one of the non-target genomes.

220
Characterization of Species-Specific Repeats in 613 Prokaryotic Species [Vol. 19, To identify repeats, the genomic sequence of the target genome was segmented into consecutive windows. The length of a window was 100 bp and the length of overlap between consecutive windows was 50 bp. We also tested different length for the overlap (75 bp) to find species-specific repeats (50 randomly chosen species was used). However, no additional species-specific repeats could be found. Thus, to achieve smaller search time, we decided to use 50 bp for overlap. Matches of sequence windows inside the target genome were identified using the similarity search software BLAST. A sequence window was classified as a candidate repeat if the following two criteria were met: the length of the BLAST match in target genome was between 85 and 115 bp and the BLAST identity between the matching region and the query window was .80%. Thereafter, in the case of finding species-specific repeats, the specificity of candidate repeats was checked. A candidate repeat was defined as a non-species-specific repeat if the length of any BLAST match in any non-target genome had a length of at least 50 bp and the identity between the match of background sequence and the query window was .60%. Finally, overlapping candidate repeats were joined to form complete (species-specific) repeats. All scripts were written in the Perl programming language.
2.4. Definition of length and number of copies of the repeat The length of different copies of the same repeat is, to some extent, variable between species or between strains of the same species. In this analysis, the median length of all copies of a given repeat is shown as the length of the repeat. To define the copy number of a given species-specific repeat, the median value of all copy numbers of that repeat over the analyzed strains of that species was used. Similarly, the median copy number of particular repeat over different species where this repeat was present was used as the copy number of a universal repeat.

Repetitive sequence functionality analysis
First, we identified prophage regions using the Phage_finder program. 25 After that, IS elements and transposons were identified from the ISfinder database 26 and from RefSeq annotations. Next, a BLAST search against plasmid sequences (2473 completely sequenced plasmid genomes downloaded from the NCBI FTP site ftp://ftp.ncbi.nih.gov/genomes/ Plasmids) was conducted to find repeats of possible plasmidic origin. Clustered regularly interspaced short palindromic repeat (CRISPR) sequences were identified within repetitive sequences with CRISPRFinder. 27 Intergenic repeat units (IRUs; ERIC/ IRU sequences) were retrieved from the collection of short repeated palindromes (http://www.pasteur.fr/ recherche/unites/pmtg/repet/intro.IRU.html). We used RefSeq gene annotation data for calculating how many species-specific and universal repetitive sequences cover genes (RNA genes and proteincoding genes) in more than 50% of the repeat length. Even if only one copy of a repeat contained a gene annotation then the repeat was considered a repeat associated with a gene (we assumed annotation faults are responsible for the missing annotation of other copies). The same logic was used to identify other features of the repeats except intergenic regions. None of the copies of intergenic repeats were allowed to contain any of the aforementioned features by more than 5% of the length of the repeated region.
Finally, Pfam database was used to predict the functions of NCBI RefSeq hypothetical proteins. 28

Definition of repeats
Throughout this paper, 'repeats' are defined as DNA sequences at least 85 bp long that are present in at least two copies in the genome of each strain of the species. Repeat copies within the same species are required to have at least 80% identity with each other. 'Species-specific repeat' is defined as a repeat that is present in only one species and there are no similar sequences (more than 60% identity over 50 bp or more) in chromosomes of other fully sequenced species. If more than one strain of a species has been sequenced, the species-specific repeat must be present in all strains, with at least two copies. The minimum length and homology cutoffs for repeat definition in this paper were chosen for practical reasons-to allow the design of species-specific polymerase chain reaction (PCR) primers for each identified repeat. 24 Our methodology is also able to identify 'strain-specific repeats', but we deliberately chose not to analyze these in order to maintain the focus of this article on species-specific repeats.
Some repeats are detected in more than one species. These can be further divided into 'universal repeats' and 'intermediate repeats'. Universal repeat is a repeat that is present in at least two species (in both as a repeat) which must belong to different genera. Thus, universal repeats are repeats that are shared between some species of separate genera. For example, widely known universal repeats are ribosomal RNA genes and elongation factor Tu. All the rest of the repeats are defined as intermediate repeats. An additional difficulty in repeat definition comes from the existence of bacterial plasmids. Typically, people search for repeats in bacterial chromosome(s), but some repeat copies can also be located in plasmids. In our data set, 232 of 613 species contained plasmids. Some of these plasmids also contained species-specific repeats. However, the analysis of species-specific plasmidic repeats is unequivocal because plasmid transfer between different bacterial species is a common phenomenon. 29,30 It has been shown that at least 51% of sequenced proteobacterial plasmids are transmissible. 9 Thus, the plasmid-based repeats are not reliable targets for characterization of bacterial species. For this reason, we have omitted plasmid sequences from our analysis and identified only the repeats that have multiple copies on chromosomes (no repetitive sequences were searched from species plasmid sequences).
As a first step, we analyzed the frequencies and lengths of all types of repeats in prokaryotic genomes. In the following analysis, we show the characteristic features of the species-specific repeats.

Characterization of repeats
3.2.1. Almost all prokaryotes contain repetitive sequences In total, 613 microbial species comprising 876 different strains were analyzed.
The number of different repeats per genome was highly variable (Fig. 1), ranging from 1 to 690 with the median 32 repeats per species. 95% of species contained fewer than 149 repeats. The highest number of repeats was found in the genome of the myxobacterium S. cellulosum' So ce 56 0 , which also has the largest size (13 Mb) of all the analyzed genomes.

The fraction of genome sequence covered by
repeats The median value of the repeat coverage was 1.8% of the bacterial genome. In only 12 bacterial species, repetitive sequences made up 7% of the total size of their genome (Fig. 2). Most of these bacteria are pathogens (insect pathogen Wolbachia sp. wRi NC_012416, plant pathogens Phytoplasma mali NC_011047 and Hamiltonella defensa NC_012751, fish pathogen Aliivibrio salmonicida NC_011312 and NC_011313, and mammalian pathogens Mycoplasma mycoides NC_005364 and NC_015431, Bartonella tribocorum NC_010161, Anaplasma phagocytophilum NC_007797, Bartonella grahamii NC_012846, and Bartonella henselae NC_005956). In general, larger genomes tend to contain larger number of repeats and also a higher fraction of repetitive sequences from the genome sequence.
The fraction of pathogens among species with elevated percentages of repetitive sequences [9 of 12  Characterization of Species-Specific Repeats in 613 Prokaryotic Species [Vol. 19, species (75%) are pathogenic] is somewhat higher than one would expect. The fraction of pathogenic species among all analyzed genomes (all that are known to be pathogenic or non-pathogenic) is 37.7%. We suggest that pathogenic bacteria have a pressure to adapt to the environment more quickly and repetitive sequences may aid this or, alternatively, be the result of this. In many pathogenic bacteria, a fraction of repetitive sequences are associated with pathogenicity genes. For example, repetitive sequences in A. phagocytophilum ( 9% of the genome comprises repeats) contain genes functionally involved with the type IV secretion system. Also, repetitive sequences in B. tribocorum and B. henselae ( 11 and 7% of the genomes are repeats, respectively) contain several genes associated with pathogenicity-the Trw-conjugation system, filamentous hemagglutinin, specific adhesions, and components of the type IV secretion system. As many pathogenic processes remain uncharacterized, the fraction of genes involved in pathogenesis may be larger than we could interpret. An example of a nonpathogenic genome containing a high fraction of repeats is Dehalococcoides ethenogenes. It contains large duplications that include genes involved in the ability of dechlorinate groundwater pollutants. Supplementary Table S1 gives an overview of the numbers of identified repetitive sequences from all analyzed prokaryotic species.

3.2.3.
Most of the repeats are species-specific From a total of 33 921 different repeats found in this work, 29 771 (88%) were species-specific (Fig. 3). The total number of universal repeats (those that are present in two or more different genera) found in this work was 554; approximately half of them were present in only two species. After analysis of the universal repeats, we estimated that only 0.25% (84 repeats) of all repeats were truly universal, present in more than six species. Measured in nucleotides, species-specific repeats constituted more than half of the total length of all repetitive sequences. Perhaps more informative is the number of nucleotides included in either universal or species-specific repeats. In Fig. 3 (second set of three columns), one can see that species-specific repeats constitute ca. 58% of the total nucleotides included in all repeats. In all analyzed genomes, there are three times more nucleotides in species-specific repeats than in universal repeats (from Fig. 3, 58% divided to 19%).
Of the analyzed species, 97% contained at least one species-specific repeat. Species-specific repetitive DNA was not detected in 20 of 613 species (listed in Supplementary Table S1). There are several possible explanations as to why we were not able to detect species-specific repeats in these 20 species. For example, B. aphidicola (six sequenced strains in our data set) is very widely defined and contains many variable genomic sequences, not sharing sufficient similarity over their repeats. Two species (Candidatus Carsonella ruddii and Candidatus Hodgkinia cicadicola) contain no repeats at all, probably because of their small genomes. In the other species (e.g. Brucella ovis, Brucella canis, and Brucella abortus), the existence of closely related fully sequenced species makes it impossible to find truly species-specific repeats.  Fig. 4, the length distribution for species-specific repeats is shown. We found that 95% of species-specific repeats were shorter than 1000 bp. Universal repeats had a similar repeat length distribution.
Although most repeats are very short, many long species-specific repeats also exist in prokaryotic genomes. The longest species-specific repeat (36 kb) belongs to Alkaliphilus metalliredigens (NC_009633), which has a genome size of 4.9 Mb. The repeat has three copies overlapping 178 gene sequences with different RefSeq gene functions, mostly annotated as hypothetical proteins. Part of this repeat (23 kb of 36 kb) contains a prophage sequence identified by the computer software Phage_finder. 25 The remaining sequences of this repeat may also have a phage origin as the genes present in this region are associated with regulation and initiation of transcription, cell wall component breakdown, and regeneration.

Most of the species-specific repeats have two
copies As the copy number of a speciesspecific repeat might differ between different strains of the same species, the median value of copies of particular repeat in the analyzed strains of particular species was used as the copy number of a given repeat. Similarly, the median copy number of particular repeat in different species was used as a copy number of a universal repeat.
Approximately 77% of species-specific repeats had two copies per repeat (Fig. 5). Only a small fraction had a large number of copies per repeat. For example, only 0.4% of species-specific repeats had more than 20 copies per repeat. The maximum number of copies of a species-specific repeat was identified in the marine bacterium Hahella chejuensis KCTC 2396 (NC_007645, genome length 7.2 Mb) with 111 copies. The median length of the copies of this repeat was 141 bp. The copies are dispersed over the genome and the origin of this sequence is unclear. The fact that several copies overlap by a few   In general, shorter repeats tend to have more copies (Fig. 6). However, Fig. 6 also shows somewhat higher copy numbers for species-specific repeats with lengths between 1000 and 2000 bp. This is due to species-specific transposons that contribute to the higher copy numbers in the length range between 900 and 1900 nucleotides.
From the repeats analyzed in this study, 11% of species-specific repeats were located in the genome in tandem. Interestingly, species-specific repeats were much more frequently located in tandem than universal repeats, 2% of which were located in tandem.

Functional analysis of species-specific repeats
3.3.1. A large fraction of the species-specific repeats are associated with protein-coding genes We analyzed 29 771 species-specific repeats in order to understand the possible functions of repetitive sequences. We were interested in the general functional categories encoded within speciesspecific and universal repeats.
Here, we have used the following classification of bacterial repetitive elements: (I) mobile genetic elements such as insertion sequences (ISs), transposons, phage sequences, and plasmid sequences; (II) RNA genes; (III) protein-coding genes; (IV) non-coding short interspersed repeats (ERIC/IRU, CRISPR, intergenic regions). All these repeat classes were represented among the species-specific repeats (Fig. 7). The protein-coding genes (class III) were the most abundant class (64% of all repeats). Mobile genetic elements (class I) and non-coding short interspersed repeats (class IV) both constituted 13 -14% of all species-specific repeats. The smallest class (class II) of repeats included tRNA and rRNA genes. Approximately 6% of species-specific repeats could not be classified with our methodology.
For comparison, the universal repeats ( Supplementary Fig. S1) contained mainly mobilome elements (class I), some long regions of universally repeated rRNA genes (class II), and 40% were protein-coding genes (class III).  We identified that 19% of species-specific repeats belong to the class of non-coding short intergenic repeats and repeats with unknown origin. A fraction of these repeats may comprise a certain conserved regulatory elements. For example, IRUs and specific intergenic repeat motifs (CIR motifs, .100 bp) are supposedly associated with gene conversion or with regulatory or structural requirements of the bacterial genomes. 31 -33 In our data set, we found that more than half of the species (365 species) contain at least one species-specific repeat which is located prior to, following to, or between RNA gene sequences (e.g. between 23S rRNA and 5S rRNA gene sequences, between different tRNA sequences, etc). The overlap between the gene sequence and the repeat constitutes a small extent from repeat/gene sequence. Thus, in many cases, sequences between RNA genes may characterize a species uniquely. Still, many intergenic repeats have no associated function and are considered to be predominantly genomic parasites. 34,35 An overview of the most frequent functions of species-specific repeats associated with genes (class III) is shown in Table 1. A notable group of species-specific repeated genes contains membrane proteins, transport proteins, and repeat-containing proteins. High sequence variation of these genes between species may be related to the need to adapt to specific environments or the need of pathogenic microbes to vary their immunogenic properties. The largest group, 31% of all repeat-associated genes, is annotated as 'hypothetical protein' by NCBI RefSeq. We analyzed the regions of hypothetical proteins containing repeats against Pfam database and found Pfam entries for 28% (2372 repeats) of them. Over the half of the repetitive sequences searched for are found from automatically generated Pfam-B families (866 families) or from the curated Pfam-A collection (471 families) predicted as proteins of unknown function. The rest of the matches represent variable functions of proteins. For example, the highest number of matches (for 48 repeats) was gained from the tetratricopeptide repeat superfamily, which are found from numerous and diverse proteins involved in such functions as cell The number shows how many repeats are overlapping the gene by at least 50% of the length of a given repeat. The total number of repeats overlapping gene sequences is 18 961, of which 15 742 repeats are shown in this table.

The function of repetitive full-length speciesspecific genes remains unknown
We analyzed repeats that contain full-length protein-coding gene sequences that are represented in at least two copies of a particular repeat. We found that 1034 repetitive sequences contain full-length species-specific protein-coding genes in at least two copies of that repeat, resulting in 4001 different genes in total. Most of these genes (56.7 versus the 30.6% expected on the basis of RefSeq annotations of prokaryotic protein-coding genes) are hypothetical proteins (statistically over-represented P , 0.0001).
However, several interesting details emerged. For example, we found that pseudogenes are statistically significantly over-represented (P , 0.0001). The RefSeq annotations of prokaryotic protein-coding genes imply that on average 1.68% of randomly chosen protein-coding genes are pseudogenes. In contrast, we found that 6.0% of the analyzed genes are pseudogenes. We suggest that this can be explained by gene duplication. It is known that gene duplication is a significant contributor to the evolution of genomes. In some cases, one duplicated copy may become non-functionalized by randomly accumulating degenerative mutations in the absence of selective advantage. 36 Another notable group of genes, encoding transport proteins, is under-represented (P , 0.0001) among the analyzed genes (12.6 and 5.8% are the percentages of expected and observed transport proteins, respectively). It is shown that at least some transport proteins are conserved among a wide range of species. 37,38 This is also supported by our analysis of universal repeats where transport proteins are over-represented.

Example of species-specific repeats of
H. defensa Hamiltonella defensa (NC_012751) is an endosymbiont found in sapfeeding insects, including aphids, psyllids, and whiteflies. Hamiltonella defensa is an example of species containing relatively various repeat classes, whereas many species contain only one or two different classes of repeats.
The complete list of species-specific repeats in H. defensa contains 95 repeats ( plus three universal repeats that contain rRNA and/or tRNA sequences and four repeats that were present in other species but not as repeats). In total, 55 species-specific repetitive sequences are classified as mobile genetic elements; 4 are classified as ISs; 34 are classified as phage-related sequences; and 17 are classified as plasmidic repeats.
The remaining 40 species-specific repeats are associated with protein-coding genes. Sixteen of them are related to proteins of the RTX family, which contains putative virulence factors, RTX toxins (repeats in toxins), which are exported proteins. 39 This could explain why there are 10 species-specific repeats associated with different inner and outer membrane proteins, some of which are suggested to be autotransporters. Also, one repeat contains genes of which the functions involve different type IV pilus biosynthesis proteins (type IV pilus proteins are also essential for virulence).
Furthermore, only two repeats contain hypothetical protein-coding genes and only four species-specific repeats contain pseudogenes.
The final six repetitive sequences are associated with genes with various functions (methionine sulfoxide reductase, phosphomannomutase, heat shock protein, amidase and lipoprotein, serine endoprotease, transcriptional regulator, and GTP cyclohydrolase). Some of these genes could also be involved with pathogenicity but no direct reference can be indicated.

Firmness of species-specific repeats
One of the major questions that arose during this work was related to the definition of species and the genetic variation within species. Does a speciesspecific repeat exist as a repeat in every newly sequenced strain? Or, if a new, closely related species is sequenced, does the repeat remain species-specific? In reality, of all the prokaryotic genomes existing on the Earth, only a small part of each has been sequenced; this part may represent distantly related genomes. Because of these questions, we have validated species-specific repeats against a prokaryotic database that consists of 1.42 times more species than the database used for finding those repeats (data not shown). Approximately 88% of detected species-specific repeats remained species-specific in this larger data set. Interestingly, closest species where the species-specific repeats could still be found are 99% identical in terms of 16S rRNA (e.g. Shewanella pealeana NC_009901 and Shewanella halifaxensis NC_010334).
Nevertheless, some repeats turned out to be non-species-specific in the larger data set, because they existed in another newly sequenced species or the repeat was not present in the genome of a newly sequenced strain of the same species. This raises the question of whether the species-specific repeats we have identified are real species-specific repeats or just data analysis artifacts, dependent on the number of available genomes and the choice of species for analysis.

No. 3]
T. Koressaar and M. Remm 227 The possibility of finding, or not finding, a speciesspecific repeat is partly determined by the controversial definition of bacterial species. Currently, the definition comprises 70% DNA -DNA re-association as a standard. 40 As this definition is somewhat arbitrary, the definition of species-specific repeats is also arbitrary. However, the problem of the species concept is not the topic of this article. We use the definition of species as it comes from the NCBI genome database, together with the genomic DNA sequence. In most cases, the repeats we define as species-specific repeats will remain specific to a group of phylogenetically related strains, even if the concept of the species will be changed in future. Yet, we have identified strain-specific repetitive sequences (data not shown) in 50 of 91 species for which 2 -26 completely sequenced strains per species were available. Strainspecific repetitive sequences add some confidence to the existence of repeats specific to a small phylogenetically related group.

Discussion
Although repetitive sequences are less common in prokaryotes than in eukaryotes, almost all prokaryotic species still have some identifiable repeats. In this work, we have characterized mainly the species-specific repetitive sequences-repeats common to all strains of a species, but not to other species.
In accordance with widely accepted understanding that prokaryotic genomes are compactly organized, we have demonstrated that the median value of the repeat coverage of bacterial genome is 1.8% (species-specific repeats constitute 1% from this). Despite of a small fraction of repetitive sequences in prokaryotic genomes, many novel findings have arisen. We have found that not all repetitive sequences in prokaryotic genomes belong to mobile genetic elements. Our results indicate that only 14% of the species-specific repetitive sequences and 53% of universal repeats are associated with mobile genetic elements. Further, the functions of many intergenic repeats and repeats from the class of unknown origin (in total 17% from the speciesspecific repeats) are not understood. However, a fraction of intergenic repeats may comprise a certain conserved regulatory elements. For example, specific intergenic repeats (CIR motifs, .100 bp) in Caulobacter crescentus may be associated with gene conversion and with gene regulation. 33 A large number of various functions of protein-coding genes, reflecting the diversity of sequenced prokaryotic genomes, have emerged from the functional analysis. In different species, genes associated with repeats have variable roles and no common phenomenon can be identified for all species or for groups of particular species. As shown, the function of many protein-coding genes associated with species-specific repeats remains unknown. We rather think that this interesting question cannot be answered by automatic analysis but must be left to biologists who investigate a particular genome in more detail. The functional analysis also revealed that different repeat classes appear in different species. In some species, phage-related sequences are in prevalence; in the other, only rDNA-related repeats appear; or in the third, one large duplication containing different protein-coding genes is found. Every species seems to be somewhat unique regarding to its repeats.
One important application of species-specific repeats is their usage in the identification of prokaryotic species. The repetitive sequences increase the sensitivity of the PCR test simply because of the higher concentration (double or more) of the initial sequence. 24 In this work, we have shown that many of the species-specific repetitive sequences are associated with protein-coding genes, so at least these repeats can be reliably used as target sequences in PCR tests. Species-specific or a group of phylogenetically related strains-specific repeats may sometimes not be the preferred target in medical diagnostics. This is in cases where pathogenic genes are known and the objective of the diagnostic test is to detect such a gene rather than to detect the species. 41 However, we believe that this does not diminish the value of species-specific repeats as a PCR target region for the detection of species. Pathogenic gene-specific PCR primers can be designed with similar methodology and used in addition to species-specific PCR primers in medical tests.
In conclusion, the current work enhances knowledge about prokaryotic species-specific repeats in many novel aspects. Most prokaryotic species contain several short species-specific repeats, which are represented in the genome with a low copy number. In a typical prokaryotic species, the speciesspecific repeats cover 1% of the genome. Over half of these repeats contain protein-coding genes, approximately one-eighth are associated with mobile genetic elements and one-eighth are associated with short non-coding interspersed elements. Thus, most of the species-specific repeats are associated with protein-coding genes and probably related to the biochemical processes required for adaptation and survival of any given prokaryote. Unfortunately, the functions of many species-specific repeats related to protein-coding genes currently remain unexplained.