BACKGROUND
staphylococcus aureus repeat  elements are short gc rich direct repeats found in intergenic regions across the s. aureus genome  <cit> .
star elements consist of 14 bp direct repeats of the consensus sequence ttgttgggccc interspersed with at least 40 bp of recurring sequences  <cit> .
the function, origin and the mechanism by which star elements propagate throughout staphylococcal genomes is unknown.
repetitive dna sequences are ubiquitous in eukaryotic and prokaryotic genomes, and are highly diverse in their structure and function.
while eukaryotic repeat elements often have no clear role within the cell, prokaryotic repeat elements tend to be functionally significant  <cit> .
these roles include transcriptional or translational phase variation of gene expression  <cit> , modulation of mrna transcript stability  <cit>  and in the case of the well characterised crispr elements protecting the genome from invading foreign dna elements  <cit> .
currently no function has been described for star elements.
repetitive elements can evolve rapidly over time.
for simple sequence repeats, such as homopolymeric tracts, slip-strand mispairing during dna replication can result in a change in repeat number after a single generation.
this is the basis of phase variable gene regulation, providing random switching of target genes between on and off states and resulting in bacterial subpopulations that are better adapted to environmental change  <cit> .
mutations in tandem repeats, resulting in changes in repeat number, occur 100– <dig>  times more frequently than point mutations making repeat arrays hotspots for genomic plasticity  <cit> .
interspersed repeats can undergo homologous recombination, resulting in changes in repeat number and the spread of a repeat element throughout the genome  <cit> .
therefore genomic repeats are inherently unstable and can undergo dramatic changes over time, which may or may not be linked to their function.
since their initial discovery over a decade ago there has been little published data regarding star elements, and much of what has been published has focused on their potential as variable number tandem repeats  and their use in s. aureus strain typing  <cit> .
information regarding the abundance and conservation of these repeat elements throughout the staphylococcus genus stem from techniques such as southern blotting that do not provide resolution to the sequence level  <cit> , and a study using comparative genomics which identified some copy number variation in a truncated star element  between a select subset of s. aureus strains  <cit> .
current advances in genome sequencing have meant that there is now a wealth of available staphylococcal genome sequences, allowing us to study the structure and evolution of star elements in much finer detail.
the purpose of this work was to analyse star elements at the molecular level from both a wide variety of s. aureus strains and in other staphylococcal species in order to further our understanding of the origin, propagation and maintenance of this repeat element.
through the use of whole genome pattern searches we have extensively mapped the locations of star elements in 15 s. aureus genomes as well as  <dig> staphylococcal species, alongside a more detailed look at individual star loci from a wider pool of s. aureus strains at the sequence level.
the data show that star elements are associated with distinct flanking genes in each staphylococcal species, suggesting that they are maintained autonomously within each species, and that their positions within each genome are stable over time.
furthermore s. aureus star elements are conserved at the sequence level within ancient evolutionary lineages.
these features point towards an as yet unidentified function for these repeat elements.
RESULTS
star elements are significantly more abundant in both repeat number and genomic location in s. aureus and s. lugdunensis compared with other staphylococci
although star elements have previously been shown to be much more abundant in s. aureus genomes than those of other staphylococcal species  <cit> , the techniques employed only provided semi-quantitative data on the actual numbers of repeat motifs involved and did not give any indication of the exact number of elements present in each genome or into how many distinct loci these fall.
the available s. aureus and staphylococcal genomes were probed in silico for the presence of the degenerate star repeat sequence tntgttgnggccc using the rsa genome-scale pattern-search tool .
the above sequence was chosen to provide enough degeneracy to identify all “true” star elements conforming to the original description from cramton et al.
<cit> , and as used in the mvla schemes  <cit>  while limiting the identification of spurious star elements.
we hypothesised that the abundance of star elements in staphylococcal species other than s. aureus would vary depending on the relatedness of that species to s. aureus, with more closely related species containing similar numbers of elements.
based on 16s rdna sequence comparison  <cit> , s. aureus is most closely related to s. epidermidis, followed by s. haemolyticus, then s. lugdunensis, s. saprophyticus, s. pseudintermedius and finally s. carnosus.
in each s. aureus strain examined, between  <dig> and  <dig> star motifs were found, occurring at  <dig> to  <dig> distinct locations in each genome  .
the number of motifs at a particular locus varied between strains; the majority of loci contain only a single repeat motif however some tracts contain as many as  <dig>  unexpectedly s. lugdunensis contains a similar abundance of star motifs to s. aureus, with  <dig> identified at  <dig> loci, while the more closely related s. epidermidis and s. haemolyticus contain far fewer than s. aureus.
s. epidermidis atcc <dig> contains  <dig> motifs at  <dig> different loci, while s. epidermidis rp62a contains  <dig> motifs at  <dig> different loci and s. haemolyticus contains  <dig> star motifs each at individual loci.
s. pseudintermedius, s. saprophyticus and s. carnosus are all devoid of star motifs.
the prevalence of these repeats is not, therefore, correlated with the phylogenetic relationships of the species, suggesting that the high levels of star motifs found in s. aureus and s. lugdunensis are due to other selective or molecular forces.
star element pattern searches were performed with an increased motif degeneracy of one additional substitution allowed throughout the sequence.
although additional, weaker star motifs were identified in each species tested, the increase in motif number was proportional to the number of “true” elements present so that s. aureus and s. lugdunenesis still showed a higher prevalence of star motifs compared with other staphylococcal species .
star elements locations are conserved within s. aureus, but not between different staphylococcal species
in order to provide insight into the evolution of star elements as species and strains diverged over time, the conservation of the positions of star loci between and within staphylococcal genomes was determined.
a total of  <dig> potential star loci were identified for s. aureus, with each strain containing between  <dig> and  <dig> loci .
strains from the same evolutionary lineage carry the same star loci, and therefore the star elements have not disseminated to new genome positions since the lineages diverged from one another.
this indicates that the elements are stable within the s. aureus genome.
the presence and number of star motifs at each potential star locus identified from each of the following s. aureus genomes:  mssa <dig>   mw <dig>   mrsa <dig>   rf <dig>   jh <dig>   jh <dig>   ed <dig>  , mu <dig>  mu <dig>   n <dig>   col,  nctc <dig>   usa <dig> frp <dig>   usa <dig> tch <dig>   newman.
* indicates only the upstream gene matches.
** indicates only the downstream gene matches.
annotations for unknown genes are taken from mssa <dig>  mrsa  <dig> and rf <dig>  full details of adjacent loci are available in table s <dig> .
the s. aureus star reference set was then used to extend this analysis to the additional staphylococcal genomes, in order to determine whether the star elements are associated with particular genes across different species.
homologues to several of the s. aureus flanking regions in the reference set were identified across the staphylococcal species, but none of these alignments contained star elements.
reference sets for both s. epidermidis and s. haemolyticus, were then used to determine star locus conservation between s. epidermidis, s. haemolyticus and s. lugdunensis.
we did not find a single star associated genomic neighbourhood that was consistent between two species, although the star associated loci were conserved between the two s. epidermidis genomes studied.
these data show that star elements have spread through and been maintained autonomously within each staphylococcal genome.
the gapr star locus differs in structure between strains but contains consistent regions of sequence variability
in order to determine how an individual star locus can alter as isolates diverge from one another, and therefore draw conclusions about how these repeat elements evolve over time, a single star locus was selected and analysed at the sequence level from a diverse pool of s. aureus strains.
the star locus found upstream of the highly conserved s. aureus glycolytic operon, which is essential for glucose metabolism  <cit> , was selected as this star locus showed high variability in the number of motifs between strains in our initial study.
the intergenic region between gapr and the upstream open reading frame was sequenced from a total of 37 s. aureus isolates from a range of sources .
the sequence of this region was also extracted from the  <dig> sequenced s. aureus genomes described above, providing data for a total of 52 s. aureus strains.
comparison of the dna sequence of the gapr star locus between s. aureus strains revealed a large amount of variability in this region, including differences in both repeat number and large scale structural changes .
in the majority of strains  the gapr star locus begins with a “start signature” sequence of gtgggacagaaatgat, which is slightly truncated compared to the sequence initially identified at the hprk star locus  <cit> .
this is followed by between  <dig> and  <dig> conserved star motifs interspersed with 40-44 bp of “spacer” sequence, which shows some variability between strains.
between the star elements and the gapr coding region there is a 380 bp “semi-variable” region, which shares approximately 88% sequence identity between strains.
this is classified as the group  <dig> structure.
in  <dig> of the strains examined  the entire star element locus is missing, as well as the first 39 bp of the 5' end of the semi-variable region .
all of the group  <dig> strains identified share 100% sequence identity across the sequenced region.
an alternative deletion event appears to have resulted in the group  <dig> structure , which retains the star start signature but shows no evidence of any star element repeat sequences .
in addition, the first 70 bp of the semi-variable region in this group shares little similarity with the semi-variable region or the star element sequences identified in any other strains.
the final two structural variants, groups 1b and 2b, appear to be derivatives of group  <dig> and  <dig> respectively.
group 2b is missing the star elements having the same precise deletion site as group  <dig>  both group 1b and group 2b have an identical 77 bp insertion within the semi-variable region , whilst group 2b has a second 37 bp insertion 27 bp upstream of the star element deletion site  .
the 37 bp insertion seen in group 2b does not share any sequence similarity with the 77 bp insertion.
star element structural groups  <dig> and  <dig> are restricted to specific evolutionary lineages
multi locus sequence typing  was used to investigate whether identify the different star element structural groups were associated with particular evolutionary lineages of s. aureus.
st types were derived for each of the strains and then a phylogenetic tree was derived using the neighbour-joining algorithm based on the mlst profiles to determine the evolutionary relationships between these strains .
these st-based phylogenetic trees indicated that the group  <dig> and group  <dig> strains, which do not contain star elements at the gapr locus, fall into distinct evolutionary lineages compared to the group  <dig> strains .
all of the group  <dig> strains , which are 100% conserved across the gapr star locus, fall into clonal complex   <dig> .
as all of the cc <dig> strains examined in this study have a group  <dig> structure, loss of star elements in these strains most probably occurred in a common ancestor and was maintained as the st’s diverged from one another.
interestingly all of the group  <dig> strains, which have a partial loss of the star element locus, belong to st  <dig> .
as the entire sequenced region is 100% conserved between the st <dig> strains, this again suggests that the deletion occurred early in the evolution of this sequence type and has been maintained in subsequent isolates.
surprisingly the group 1b and group 2b strains, which contain the same unique sequence insertion, fall into distinct clonal complexes with very different allelic profiles; the group 1b strains are from st <dig>  and the group 2b strains are from st <dig> and novel st a, which are both in cc <dig>  although it initially appeared that the group 1b and group 2b structures were derived from group  <dig> and group  <dig>  the phylogenetic data indicates that this is not a recent event.
furthermore these structures did not occur due to a recent loss/gain of the star locus between groups 1b and 2b, as these strains are from different cc’s.
taken together these data suggest that the gapr star locus differences occurred in very early lineages of s. aureus and have been maintained at a level equal to that of cc in subsequent strains.
sequence variation of the star element spacers correlates with evolutionary lineage
as the groups  <dig> strains fall into a wide range of st’s and cc’s, it is clear that gapr star locus structure alone does not correlate with any particular lineage.
however analysis of this star locus at the sequence level shows that the sequences of the “spacers”, which occur between star motifs, are strongly conserved within cc’s.
for example, in strains from cc <dig> and cc <dig> the star spacer sequences are 100% identical between isolates even though the number of repeat motifs varies .
interestingly the final spacer sequence  is distinct from the internal spacers, but this “anchor” spacer is again 100% conserved between strains of the same lineage.
in contrast, alignment of the spacers from strains originating from different lineages, even where they contain the same repeat number, detected high levels of variation in these sequences between distinct cc’s .
we have confirmed the conservation of spacer sequences within a cc’s in all strains tested here, with the exception of the two strains representing cc <dig>  the spacer sequences from the cc <dig> strains c <dig> and c <dig> are still highly conserved, but they are not 100% identical.
this is further evidence that the structure and sequence of the gapr star locus is maintained within distinct evolutionary lineages.
star spacer sequences are distinct at different loci within s. aureus strains but still correlate with lineage
two additional star loci were analysed to further investigate the link between star element conservation and evolutionary lineage.
the star loci found upstream of both the hprk gene, encoding a hpr kinase/phosphorylase, and a gene of unknown function sas <dig>  referred to as orf <dig> in this study, were chosen as rsat analysis of these regions shows that they both contained variable numbers of star motifs and are preceded by a start signature.
the star element regions upstream of hprk and orf <dig> were either pcr amplified and sequenced from a selection of s. aureus strains using primer pairs hprk f + hprk r and orf <dig> f + orf <dig> r respectively  or extracted from the  <dig> complete genome sequences.
the strains were chosen to include at least  <dig> examples, where possible, of strains from each lineage identified previously .
interestingly both the hprk and orf <dig> star loci have some key structural differences to that of the gapr star locus.
the star start signature sequence is present at both loci but occurs ~130 bp and 188 bp upstream of the first repeat motif at the hprk and orf <dig> loci respectively, compared to ~70 bp at the gapr star locus .
furthermore there is no evidence for different structural variants in any of the strains examined as both the hprk and orf <dig> star elements only follow the group  <dig> star element structure found at the gapr locus.
there is also less variability in the number of star element repeat motifs at each of these loci, with the hprk locus ranging from 1– <dig> repeats and the orf <dig> locus ranging from 3– <dig>  compared with the 1– <dig> repeats seen at the gapr locus.
sequence analysis of the hprk and orf <dig> star spacers showed that sequence level variation in these repeat regions still strongly correlates with lineage as seen at the gapr locus.
alignments of each individual locus clearly demonstrate high levels of conservation of the spacer sequences within strains from a particular lineage , as shown for cc <dig> .
for strains containing multiple star repeats at locus orf <dig>  we observed two distinct spacer types within the same locus in some strains, as seen in cc <dig> .
however it is important to note that these sequences are still 100% conserved within each lineage and do not occur at either the hprk or gapr loci in any of the strains examined supporting the observation that the spacer sequences are distinct from one another and that there is no frequent transfer of motifs/spacers between the star loci.
discussion
in this study we have taken advantage of the wealth of fully annotated staphylococcal genomes to take a detailed look at star elements.
to our knowledge this is the first in depth study of these interspersed repeats at the sequence level across multiple staphylococcal species, providing a unique insight into their evolution.
star elements are highly abundant in s. aureus and yet we have shown that strain variation in the star element nucleotide sequences strongly correlates with their evolutionary lineage, as derived by mlst.
this is unexpected as intergenic regions such as the star loci, which consist of repetitive elements dispersed throughout the genome, would be expected to show a high level of mutation and hence evolve at a higher rate than the conserved functional mlst loci where mutations are observed at a very low rate  <cit> .
these findings suggest star elements are functional and may be under strong purifying selection.
star elements were sequenced from the gapr, hprk and orf <dig> loci from multiple s. aureus strains.
in the majority of loci where multiple star repeats were present, the spacer sequences were often identical or differed by 1– <dig> nucleotides resulting in tandem repeats of ~ <dig> nucleotides.
these repetitive sequences should be unstable and exhibit frequent alterations in repeat number due to slip-strand mispairing during dna replication.
this process is likely to drive rapid alterations in repeat number, but not sequence, at many of these loci, as found with some other bacterial tandem repeats  <cit> .
congruent with this theory, strains belonging to the same st contain identical or highly conserved spacer sequences between the interspersed star motifs at a specific locus even when repeat numbers varied.
this also suggests that localised expansion and contraction of the repeat region occurs even as the strains diverge from one another.
in contrast, the spacer sequences are distinct at each star locus, even within a particular genome.
due to the repetitive nature of star elements it has previously been suggested that homologous recombination between repeats occurs as a means of large scale genomic rearrangements  <cit> , or could provide a simple means of propagating these repeats at different loci throughout the genome.
as the spacers are distinct between unrelated strains and at different star loci within a strain, homologous recombination is unlikely to be occurring at a high frequency between star loci either intergenomically or intragenomically.
either of these processes would result in gene conversion and the emergence of a dominant spacer sequence variant across multiple loci, a phenomenon we did not identify in this study.
from the evidence presented here we suggest that the process of varying repeat number within a locus is limited to duplication or deletion of motifs from within that locus during dna replication or repair and is not due to recombination with elements present elsewhere in the genome.
we also suggest that the mechanism for dispersal of the star elements to new positions throughout the s. aureus genome may not involve recombination as originally hypothesised.
the gapr star locus was the least structurally stable of the three loci studied.
the loss of the elements in the group  <dig> and 2b structure occurs at the same “deletion” site and the surrounding dna is undisturbed compared to that of the group  <dig> and 1b strains.
this is similar to another class of interspersed bacterial repeats known as enterobacterial repetitive intergenic consensus  sequences, which have been identified across the eubacterial kingdom  <cit> .
the sequence surrounding an inserted eric remains unchanged, indicating a precise insertion or deletion event via a mechanism distinct from classic transposition mechanisms  <cit> .
it is unclear whether a similar conserved mechanism is involved in the total loss or gain of star loci or whether the deletion site is merely acting as a hotspot for star element translocation.
the partial loss of elements seen in strains such as rf <dig>  does not occur at this deletion site, and may represent a different mechanism of repeat propagation or an error in repeat translocation in an ancestral strain that has been maintained in subsequent generations.
there is no evidence of the total loss or gain of the gapr star locus in the recent evolution of s. aureus strains, as both the group  <dig> and group  <dig> isolates fall into distinct evolutionarily lineages.
this strongly implies that the deletion process is infrequent and that the loss or gain of the gapr star locus may have occurred in early ancestors of these lineages and been retained in subsequent isolates.
pourcel et al.
observed a similar complex structure for the star elements in the sa <dig> locus  with restriction of specific structural variants to certain lineages  <cit> .
these findings provide further evidence of the conservation of each of the star loci within a strain and lineage.
our observed correlation between evolutionary lineage and both the structure of the gapr locus and the spacer sequences of the gapr, hprk and orf <dig> loci, suggests that star element loci retain lineage-specific phylogenetic information and may be utilised as major determinants of lineage in typing schemes.
the genome wide mapping of star elements across the 15 s. aureus strains studied here identified  <dig> loci that were present in every genome sequence and a further  <dig> loci that were present in 85% of the genome sequences.
the vast majority of these loci  contain more than one repeat and exhibit variable repeat numbers , making them prime candidates for the development of future typing schemes.
some star loci have already been utilised in typing schemes for s. aureus, first using an rflp typing method  <cit> , and more recently as part of a greater multiple-locus variable-number tandem-repeat analysis  scheme alongside other variable-number tandem repeats  and staphylococcal interspersed repeat units   <cit> .
the recent extended mlva scheme utilised six star element loci of which five were completely conserved in a collection of  <dig> strains  <cit> , although only four are present in up to 85% of the strains studied here.
therefore our highly conserved loci should be examined for their potential value as markers of lineages.
we have found that the star elements are not restricted to specific genomic neighbourhoods across staphylococcal species.
this would suggest that the elements are not simply decaying from some early staphylococcus progenitor as this genus has diverged over time, but rather that each species has acquired star elements as independent events, which have then undergone proliferation to distinct locations in each genome.
furthermore, star elements are maintained at a much higher level in the s. aureus and s. lugdunensis genomes compared to other staphylococcal species.
the higher prevalence of these elements in s. aureus and s. lugdunensis may be due to the presence of a dispersal mechanism  that is absent in the other species studied here, the absence of a mechanism to prevent spread of repetitive elements in these two species or strong selection for the function of these elements.
the highly conserved nature of star elements within a cc suggests a functional role.
unlike eukaryotic genomes which can contain more than 50% repetitive dna  <cit> , prokaryotic genomes are streamlined as the propagation of non-functional “selfish” dna is a burden to the rapidly dividing organisms and selected against  <cit> .
other repeat elements in bacteria have functions in cell physiology, such as transcriptional control  <cit>  and protection of the microbial genome against foreign dna  <cit> .
a functional role for star elements is supported by evidence showing that some star elements are present in the leader regions of mrnas although the significance of this for gene expression has yet to be investigated further  <cit> .
alternatively, these repetitive sequences may have a general function in chromosome structure or stability, as seen with some eukaryotic repeat elements  <cit> , which has led to their maintenance and spread within staphylococcal genomes.
the star repeats are found associated with loci encoding virulence factors, metal transporters and several essential metabolic enzymes.
the significance of the star repeats in the intergenic regions of these particular loci requires further investigation.
interestingly, both s. aureus and s. lugdunensis tend to be much more pathogenic in humans compared to other staphylococcal species  <cit>  with s. lugdunensis n <dig> having several homologues of s. aureus virulence and colonisation factors that are not found in other staphylococcal species  <cit> .
our finding that star elements are present in higher levels in two of the more virulent staphylococcal species may indicate that the star elements play a role in pathogenesis.
with the huge increase in the number of available genome sequences, the occurrence of star repeats in other bacterial species requires further investigation to confirm their existence and function outside of the staphylococcal genus.
CONCLUSIONS
star elements are highly conserved at the sequence level and are maintained at high levels in both s. aureus and s. lugdunensis, but not in the other staphylococcal species studied here.
furthermore star elements are conserved at the sequence level within distinct evolutionary lineages but conversely exhibit localised expansion and contraction of repeats.
this means that these repeat loci retain both ancient and more recent phylogenetic information, making them ideal candidates for strain typing schemes.
the high level of conservation seen in these repeats suggests that star elements may, as with other bacterial repeats, have a functional role in cell physiology and confer fitness advantages on some or all s. aureus lineages.

 staphylococcus aureus repeat elements are short gc rich direct repeats found in intergenic regions across the s. aureus genome. star elements consist of 14 bp direct repeats of the consensus sequence ttgttgggccc interspersed with at least 40 bp of recurring sequences. the function, origin and the mechanism by which star elements propagate throughout staphylococcal genomes is unknown. repetitive dna sequences are ubiquitous in eukaryotic and prokaryotic genomes, and are highly diverse in their structure and function. while eukaryotic repeat elements often have no clear role within the cell, prokaryotic repeat elements tend to be functionally significant. these roles include transcriptional or translational phase variation of gene expression, modulation of mrna transcript stability and in the case of the well characterised crispr elements protecting the genome from invading foreign dna elements. currently no function has been described for star elements. repetitive elements can evolve rapidly over time. for simple sequence repeats, such as homopolymeric tracts, slip-strand mispairing during dna replication can result in a change in repeat number after a single generation. this is the basis of phase variable gene regulation, providing random switching of target genes between on and off states and resulting in bacterial subpopulations that are better adapted to environmental change. mutations in tandem repeats, resulting in changes in repeat number, occur 100– times more frequently than point mutations making repeat arrays hotspots for genomic plasticity. interspersed repeats can undergo homologous recombination, resulting in changes in repeat number and the spread of a repeat element throughout the genome. therefore genomic repeats are inherently unstable and can undergo dramatic changes over time, which may or may not be linked to their function. since their initial discovery over a decade ago there has been little published data regarding star elements, and much of what has been published has focused on their potential as variable number tandem repeats and their use in s. aureus strain typing. information regarding the abundance and conservation of these repeat elements throughout the staphylococcus genus stem from techniques such as southern blotting that do not provide resolution to the sequence level, and a study using comparative genomics which identified some copy number variation in a truncated star element between a select subset of s. aureus strains. current advances in genome sequencing have meant that there is now a wealth of available staphylococcal genome sequences, allowing us to study the structure and evolution of star elements in much finer detail. the purpose of this work was to analyse star elements at the molecular level from both a wide variety of s. aureus strains and in other staphylococcal species in order to further our understanding of the origin, propagation and maintenance of this repeat element. through the use of whole genome pattern searches we have extensively mapped the locations of star elements in 15 s. aureus genomes as well as staphylococcal species, alongside a more detailed look at individual star loci from a wider pool of s. aureus strains at the sequence level. the data show that star elements are associated with distinct flanking genes in each staphylococcal species, suggesting that they are maintained autonomously within each species, and that their positions within each genome are stable over time. furthermore s. aureus star elements are conserved at the sequence level within ancient evolutionary lineages. these features point towards an as yet unidentified function for these repeat elements.  star elements are significantly more abundant in both repeat number and genomic location in s. aureus and s. lugdunensis compared with other staphylococci although star elements have previously been shown to be much more abundant in s. aureus genomes than those of other staphylococcal species, the techniques employed only provided semi-quantitative data on the actual numbers of repeat motifs involved and did not give any indication of the exact number of elements present in each genome or into how many distinct loci these fall. the available s. aureus and staphylococcal genomes were probed in silico for the presence of the degenerate star repeat sequence tntgttgnggccc using the rsa genome-scale pattern-search tool. the above sequence was chosen to provide enough degeneracy to identify all “true” star elements conforming to the original description from cramton et al., and as used in the mvla schemes while limiting the identification of spurious star elements. we hypothesised that the abundance of star elements in staphylococcal species other than s. aureus would vary depending on the relatedness of that species to s. aureus, with more closely related species containing similar numbers of elements. based on 16s rdna sequence comparison, s. aureus is most closely related to s. epidermidis, followed by s. haemolyticus, then s. lugdunensis, s. saprophyticus, s. pseudintermedius and finally s. carnosus. in each s. aureus strain examined, between and star motifs were found, occurring at to distinct locations in each genome. the number of motifs at a particular locus varied between strains; the majority of loci contain only a single repeat motif however some tracts contain as many as unexpectedly s. lugdunensis contains a similar abundance of star motifs to s. aureus, with identified at loci, while the more closely related s. epidermidis and s. haemolyticus contain far fewer than s. aureus. s. epidermidis atcc contains motifs at different loci, while s. epidermidis rp62a contains motifs at different loci and s. haemolyticus contains star motifs each at individual loci. s. pseudintermedius, s. saprophyticus and s. carnosus are all devoid of star motifs. the prevalence of these repeats is not, therefore, correlated with the phylogenetic relationships of the species, suggesting that the high levels of star motifs found in s. aureus and s. lugdunensis are due to other selective or molecular forces. star element pattern searches were performed with an increased motif degeneracy of one additional substitution allowed throughout the sequence. although additional, weaker star motifs were identified in each species tested, the increase in motif number was proportional to the number of “true” elements present so that s. aureus and s. lugdunenesis still showed a higher prevalence of star motifs compared with other staphylococcal species. star elements locations are conserved within s. aureus, but not between different staphylococcal species in order to provide insight into the evolution of star elements as species and strains diverged over time, the conservation of the positions of star loci between and within staphylococcal genomes was determined. a total of potential star loci were identified for s. aureus, with each strain containing between and loci. strains from the same evolutionary lineage carry the same star loci, and therefore the star elements have not disseminated to new genome positions since the lineages diverged from one another. this indicates that the elements are stable within the s. aureus genome. the presence and number of star motifs at each potential star locus identified from each of the following s. aureus genomes: mssa mw mrsa rf jh jh ed, mu mu n col, nctc usa frp usa tch newman. * indicates only the upstream gene matches. ** indicates only the downstream gene matches. annotations for unknown genes are taken from mssa mrsa and rf full details of adjacent loci are available in table s. the s. aureus star reference set was then used to extend this analysis to the additional staphylococcal genomes, in order to determine whether the star elements are associated with particular genes across different species. homologues to several of the s. aureus flanking regions in the reference set were identified across the staphylococcal species, but none of these alignments contained star elements. reference sets for both s. epidermidis and s. haemolyticus, were then used to determine star locus conservation between s. epidermidis, s. haemolyticus and s. lugdunensis. we did not find a single star associated genomic neighbourhood that was consistent between two species, although the star associated loci were conserved between the two s. epidermidis genomes studied. these data show that star elements have spread through and been maintained autonomously within each staphylococcal genome. the gapr star locus differs in structure between strains but contains consistent regions of sequence variability in order to determine how an individual star locus can alter as isolates diverge from one another, and therefore draw conclusions about how these repeat elements evolve over time, a single star locus was selected and analysed at the sequence level from a diverse pool of s. aureus strains. the star locus found upstream of the highly conserved s. aureus glycolytic operon, which is essential for glucose metabolism, was selected as this star locus showed high variability in the number of motifs between strains in our initial study. the intergenic region between gapr and the upstream open reading frame was sequenced from a total of 37 s. aureus isolates from a range of sources. the sequence of this region was also extracted from the sequenced s. aureus genomes described above, providing data for a total of 52 s. aureus strains. comparison of the dna sequence of the gapr star locus between s. aureus strains revealed a large amount of variability in this region, including differences in both repeat number and large scale structural changes. in the majority of strains the gapr star locus begins with a “start signature” sequence of gtgggacagaaatgat, which is slightly truncated compared to the sequence initially identified at the hprk star locus. this is followed by between and conserved star motifs interspersed with 40-44 bp of “spacer” sequence, which shows some variability between strains. between the star elements and the gapr coding region there is a 380 bp “semi-variable” region, which shares approximately 88% sequence identity between strains. this is classified as the group structure. in of the strains examined the entire star element locus is missing, as well as the first 39 bp of the 5' end of the semi-variable region. all of the group strains identified share 100% sequence identity across the sequenced region. an alternative deletion event appears to have resulted in the group structure, which retains the star start signature but shows no evidence of any star element repeat sequences. in addition, the first 70 bp of the semi-variable region in this group shares little similarity with the semi-variable region or the star element sequences identified in any other strains. the final two structural variants, groups 1b and 2b, appear to be derivatives of group and respectively. group 2b is missing the star elements having the same precise deletion site as group both group 1b and group 2b have an identical 77 bp insertion within the semi-variable region, whilst group 2b has a second 37 bp insertion 27 bp upstream of the star element deletion site. the 37 bp insertion seen in group 2b does not share any sequence similarity with the 77 bp insertion. star element structural groups and are restricted to specific evolutionary lineages multi locus sequence typing was used to investigate whether identify the different star element structural groups were associated with particular evolutionary lineages of s. aureus. st types were derived for each of the strains and then a phylogenetic tree was derived using the neighbour-joining algorithm based on the mlst profiles to determine the evolutionary relationships between these strains. these st-based phylogenetic trees indicated that the group and group strains, which do not contain star elements at the gapr locus, fall into distinct evolutionary lineages compared to the group strains. all of the group strains, which are 100% conserved across the gapr star locus, fall into clonal complex. as all of the cc strains examined in this study have a group structure, loss of star elements in these strains most probably occurred in a common ancestor and was maintained as the st’s diverged from one another. interestingly all of the group strains, which have a partial loss of the star element locus, belong to st. as the entire sequenced region is 100% conserved between the st strains, this again suggests that the deletion occurred early in the evolution of this sequence type and has been maintained in subsequent isolates. surprisingly the group 1b and group 2b strains, which contain the same unique sequence insertion, fall into distinct clonal complexes with very different allelic profiles; the group 1b strains are from st and the group 2b strains are from st and novel st a, which are both in cc although it initially appeared that the group 1b and group 2b structures were derived from group and group the phylogenetic data indicates that this is not a recent event. furthermore these structures did not occur due to a recent loss/gain of the star locus between groups 1b and 2b, as these strains are from different cc’s. taken together these data suggest that the gapr star locus differences occurred in very early lineages of s. aureus and have been maintained at a level equal to that of cc in subsequent strains. sequence variation of the star element spacers correlates with evolutionary lineage as the groups strains fall into a wide range of st’s and cc’s, it is clear that gapr star locus structure alone does not correlate with any particular lineage. however analysis of this star locus at the sequence level shows that the sequences of the “spacers”, which occur between star motifs, are strongly conserved within cc’s. for example, in strains from cc and cc the star spacer sequences are 100% identical between isolates even though the number of repeat motifs varies. interestingly the final spacer sequence is distinct from the internal spacers, but this “anchor” spacer is again 100% conserved between strains of the same lineage. in contrast, alignment of the spacers from strains originating from different lineages, even where they contain the same repeat number, detected high levels of variation in these sequences between distinct cc’s. we have confirmed the conservation of spacer sequences within a cc’s in all strains tested here, with the exception of the two strains representing cc the spacer sequences from the cc strains c and c are still highly conserved, but they are not 100% identical. this is further evidence that the structure and sequence of the gapr star locus is maintained within distinct evolutionary lineages. star spacer sequences are distinct at different loci within s. aureus strains but still correlate with lineage two additional star loci were analysed to further investigate the link between star element conservation and evolutionary lineage. the star loci found upstream of both the hprk gene, encoding a hpr kinase/phosphorylase, and a gene of unknown function sas referred to as orf in this study, were chosen as rsat analysis of these regions shows that they both contained variable numbers of star motifs and are preceded by a start signature. the star element regions upstream of hprk and orf were either pcr amplified and sequenced from a selection of s. aureus strains using primer pairs hprk f + hprk r and orf f + orf r respectively or extracted from the complete genome sequences. the strains were chosen to include at least examples, where possible, of strains from each lineage identified previously. interestingly both the hprk and orf star loci have some key structural differences to that of the gapr star locus. the star start signature sequence is present at both loci but occurs ~130 bp and 188 bp upstream of the first repeat motif at the hprk and orf loci respectively, compared to ~70 bp at the gapr star locus. furthermore there is no evidence for different structural variants in any of the strains examined as both the hprk and orf star elements only follow the group star element structure found at the gapr locus. there is also less variability in the number of star element repeat motifs at each of these loci, with the hprk locus ranging from 1– repeats and the orf locus ranging from 3– compared with the 1– repeats seen at the gapr locus. sequence analysis of the hprk and orf star spacers showed that sequence level variation in these repeat regions still strongly correlates with lineage as seen at the gapr locus. alignments of each individual locus clearly demonstrate high levels of conservation of the spacer sequences within strains from a particular lineage, as shown for cc. for strains containing multiple star repeats at locus orf we observed two distinct spacer types within the same locus in some strains, as seen in cc. however it is important to note that these sequences are still 100% conserved within each lineage and do not occur at either the hprk or gapr loci in any of the strains examined supporting the observation that the spacer sequences are distinct from one another and that there is no frequent transfer of motifs/spacers between the star loci. discussion in this study we have taken advantage of the wealth of fully annotated staphylococcal genomes to take a detailed look at star elements. to our knowledge this is the first in depth study of these interspersed repeats at the sequence level across multiple staphylococcal species, providing a unique insight into their evolution. star elements are highly abundant in s. aureus and yet we have shown that strain variation in the star element nucleotide sequences strongly correlates with their evolutionary lineage, as derived by mlst. this is unexpected as intergenic regions such as the star loci, which consist of repetitive elements dispersed throughout the genome, would be expected to show a high level of mutation and hence evolve at a higher rate than the conserved functional mlst loci where mutations are observed at a very low rate. these findings suggest star elements are functional and may be under strong purifying selection. star elements were sequenced from the gapr, hprk and orf loci from multiple s. aureus strains. in the majority of loci where multiple star repeats were present, the spacer sequences were often identical or differed by 1– nucleotides resulting in tandem repeats of ~ nucleotides. these repetitive sequences should be unstable and exhibit frequent alterations in repeat number due to slip-strand mispairing during dna replication. this process is likely to drive rapid alterations in repeat number, but not sequence, at many of these loci, as found with some other bacterial tandem repeats. congruent with this theory, strains belonging to the same st contain identical or highly conserved spacer sequences between the interspersed star motifs at a specific locus even when repeat numbers varied. this also suggests that localised expansion and contraction of the repeat region occurs even as the strains diverge from one another. in contrast, the spacer sequences are distinct at each star locus, even within a particular genome. due to the repetitive nature of star elements it has previously been suggested that homologous recombination between repeats occurs as a means of large scale genomic rearrangements, or could provide a simple means of propagating these repeats at different loci throughout the genome. as the spacers are distinct between unrelated strains and at different star loci within a strain, homologous recombination is unlikely to be occurring at a high frequency between star loci either intergenomically or intragenomically. either of these processes would result in gene conversion and the emergence of a dominant spacer sequence variant across multiple loci, a phenomenon we did not identify in this study. from the evidence presented here we suggest that the process of varying repeat number within a locus is limited to duplication or deletion of motifs from within that locus during dna replication or repair and is not due to recombination with elements present elsewhere in the genome. we also suggest that the mechanism for dispersal of the star elements to new positions throughout the s. aureus genome may not involve recombination as originally hypothesised. the gapr star locus was the least structurally stable of the three loci studied. the loss of the elements in the group and 2b structure occurs at the same “deletion” site and the surrounding dna is undisturbed compared to that of the group and 1b strains. this is similar to another class of interspersed bacterial repeats known as enterobacterial repetitive intergenic consensus sequences, which have been identified across the eubacterial kingdom. the sequence surrounding an inserted eric remains unchanged, indicating a precise insertion or deletion event via a mechanism distinct from classic transposition mechanisms. it is unclear whether a similar conserved mechanism is involved in the total loss or gain of star loci or whether the deletion site is merely acting as a hotspot for star element translocation. the partial loss of elements seen in strains such as rf does not occur at this deletion site, and may represent a different mechanism of repeat propagation or an error in repeat translocation in an ancestral strain that has been maintained in subsequent generations. there is no evidence of the total loss or gain of the gapr star locus in the recent evolution of s. aureus strains, as both the group and group isolates fall into distinct evolutionarily lineages. this strongly implies that the deletion process is infrequent and that the loss or gain of the gapr star locus may have occurred in early ancestors of these lineages and been retained in subsequent isolates. pourcel et al. observed a similar complex structure for the star elements in the sa locus with restriction of specific structural variants to certain lineages. these findings provide further evidence of the conservation of each of the star loci within a strain and lineage. our observed correlation between evolutionary lineage and both the structure of the gapr locus and the spacer sequences of the gapr, hprk and orf loci, suggests that star element loci retain lineage-specific phylogenetic information and may be utilised as major determinants of lineage in typing schemes. the genome wide mapping of star elements across the 15 s. aureus strains studied here identified loci that were present in every genome sequence and a further loci that were present in 85% of the genome sequences. the vast majority of these loci contain more than one repeat and exhibit variable repeat numbers, making them prime candidates for the development of future typing schemes. some star loci have already been utilised in typing schemes for s. aureus, first using an rflp typing method, and more recently as part of a greater multiple-locus variable-number tandem-repeat analysis scheme alongside other variable-number tandem repeats and staphylococcal interspersed repeat units. the recent extended mlva scheme utilised six star element loci of which five were completely conserved in a collection of strains, although only four are present in up to 85% of the strains studied here. therefore our highly conserved loci should be examined for their potential value as markers of lineages. we have found that the star elements are not restricted to specific genomic neighbourhoods across staphylococcal species. this would suggest that the elements are not simply decaying from some early staphylococcus progenitor as this genus has diverged over time, but rather that each species has acquired star elements as independent events, which have then undergone proliferation to distinct locations in each genome. furthermore, star elements are maintained at a much higher level in the s. aureus and s. lugdunensis genomes compared to other staphylococcal species. the higher prevalence of these elements in s. aureus and s. lugdunensis may be due to the presence of a dispersal mechanism that is absent in the other species studied here, the absence of a mechanism to prevent spread of repetitive elements in these two species or strong selection for the function of these elements. the highly conserved nature of star elements within a cc suggests a functional role. unlike eukaryotic genomes which can contain more than 50% repetitive dna, prokaryotic genomes are streamlined as the propagation of non-functional “selfish” dna is a burden to the rapidly dividing organisms and selected against. other repeat elements in bacteria have functions in cell physiology, such as transcriptional control and protection of the microbial genome against foreign dna. a functional role for star elements is supported by evidence showing that some star elements are present in the leader regions of mrnas although the significance of this for gene expression has yet to be investigated further. alternatively, these repetitive sequences may have a general function in chromosome structure or stability, as seen with some eukaryotic repeat elements, which has led to their maintenance and spread within staphylococcal genomes. the star repeats are found associated with loci encoding virulence factors, metal transporters and several essential metabolic enzymes. the significance of the star repeats in the intergenic regions of these particular loci requires further investigation. interestingly, both s. aureus and s. lugdunensis tend to be much more pathogenic in humans compared to other staphylococcal species with s. lugdunensis n having several homologues of s. aureus virulence and colonisation factors that are not found in other staphylococcal species. our finding that star elements are present in higher levels in two of the more virulent staphylococcal species may indicate that the star elements play a role in pathogenesis. with the huge increase in the number of available genome sequences, the occurrence of star repeats in other bacterial species requires further investigation to confirm their existence and function outside of the staphylococcal genus.  star elements are highly conserved at the sequence level and are maintained at high levels in both s. aureus and s. lugdunensis, but not in the other staphylococcal species studied here. furthermore star elements are conserved at the sequence level within distinct evolutionary lineages but conversely exhibit localised expansion and contraction of repeats. this means that these repeat loci retain both ancient and more recent phylogenetic information, making them ideal candidates for strain typing schemes. the high level of conservation seen in these repeats suggests that star elements may, as with other bacterial repeats, have a functional role in cell physiology and confer fitness advantages on some or all s. aureus lineages. 
unexpectedly, sequencing of the internal spacer sequences of individual repeat elements from multiple isolates showed conservation at the sequence level within deep evolutionary lineages of s. aureus. using bioinformatic analyses, we found that the star elements were located in different genomic loci within each staphylococcal species. there was no correlation between the number of star elements in each genome and the evolutionary relatedness of staphylococcal species, however higher levels of repeats were observed in both s. aureus and s. lugdunensis compared to other staphylococcal species. the high degree of lineage and locus-specific conservation of these intergenic repeat regions suggests that star elements are maintained due to selective or molecular forces with some of these elements having an important role in cell physiology. the high prevalence in two of the more virulent staphylococcal species is indicative of a potential role for star elements in pathogenesis. in this study the conservation and variation in these elements was explored by bioinformatic analyses of published staphylococcal genome sequences and through sequencing of specific star element loci from a large set of s. aureus isolates.
 staphylococcus aureus repeat elements are a type of interspersed intergenic direct repeat. in this study the conservation and variation in these elements was explored by bioinformatic analyses of published staphylococcal genome sequences and through sequencing of specific star element loci from a large set of s. aureus isolates.  using bioinformatic analyses, we found that the star elements were located in different genomic loci within each staphylococcal species. there was no correlation between the number of star elements in each genome and the evolutionary relatedness of staphylococcal species, however higher levels of repeats were observed in both s. aureus and s. lugdunensis compared to other staphylococcal species. unexpectedly, sequencing of the internal spacer sequences of individual repeat elements from multiple isolates showed conservation at the sequence level within deep evolutionary lineages of s. aureus. whilst individual star element loci were demonstrated to expand and contract, the sequences associated with each locus were stable and distinct from one another.  the high degree of lineage and locus-specific conservation of these intergenic repeat regions suggests that star elements are maintained due to selective or molecular forces with some of these elements having an important role in cell physiology. the high prevalence in two of the more virulent staphylococcal species is indicative of a potential role for star elements in pathogenesis. 