BACKGROUND
two-component systems  are the dominant phosphorylation based signalling pathways in prokaryotes.
the average eubacterial genome encodes ~ <dig> tcs proteins  <cit> , and they regulate diverse and important behaviours, including pathcomparative genomics and experimentalogenesis, cell cycle progression, sporulation and drug resistance  <cit> .
typical tcss consist of two proteins – a histidine kinase , which can sense an environmental parameter through an ‘input’ domain, and a response regulator , which pocomparative genomics and experimentalssesses an effector function  via an ‘output’ domain.
in response to an environmental signal the hk autophosphorylates within its ‘transmitter’ domain and the phosphoryl group is then transferred onto the ‘receiver’ domain of the partner rr, modulating its output domain effector activity.
common modifications to the basic scheme include hybrid kinases – proteins which contain both transmitter and receiver domains, and phosphotransfer proteins  - which can shuttle phosphoryl groups between receiver domains forming a phosphorelay.
for comprehensive reviews of tcs biology, refer to the recent book by gross and beier  <cit>  and the special issue of current opinion in microbiology edited by bourret and silversmith  <cit> .
comparative genomics and experimental approaches have shed light on the general principles of tcs evolution  <cit> .
for instance tcss have evolved primarily by duplication and divergence of ancestral systems, with the contemporary diversity of tcs architecture generated by domain shuffling.
this is important, as an understanding of natural variation due to evolution can inform the rational engineering of tcss for biotechnological and research applications  <cit> .
nevertheless there are important questions left unanswered regarding tcs evolution and most studies have been forced to make several assumptions in their analyses.
for instance, inferences regarding tcs evolution have been based on the comparison of various extant genomes.
common approaches are to use all available genomes, or to use multiple organisms from a single taxonomic unit.
in both cases this typically results in comparisons between genomes which exhibit different relatedness to one another.
thus relatively distant and close genomes are considered equally, potentially distorting the conclusions drawn.
this is particularly problematic for studies of tcs evolution, as it is thought that different evolutionary processes work at different timescales.
for example, shuffling of hk input domains seems to occur rapidly post gene duplication, and gene gain/loss apparently dominates over gene fusion/fission at long timescales  <cit> .
in addition, studies are also often limited in power as there is sequencing bias within/between taxa, most organisms contain just dozens of tcs proteins, and at most taxonomic levels there is a relative shortfall in the number of available genomes, all of which can distort conclusions .
myxobacteria  are well known for their large numbers of tcss, many of which have been characterised experimentally .
they exhibit complex developmental and predatory behaviours and are of ongoing importance to the biotechnology industry .
a major focus of myxobacterial research has been elucidation of the molecular mechanisms underpinning multicellular development, which are surprisingly plastic  <cit> , and involve large numbers of tcs proteins  <cit> .
at the time of analysis, genome sequences were available for  <dig> myxobacteria .
for ease of reference, genomes are abbreviated in the remaining text as shown in table  <dig>  two of these genome sequences  are incomplete, but were included in the study anyway, as they were nearly complete, being composed of relatively small numbers of contigs .
bdellovibrio bacteriovorus hd <dig>  was also included in the study as the bdellovibrionales are the closest known relatives of the myxobacteria.table  <dig> order myxococcales – genome properties, taxonomy and ecology
myxococcus fulvus hw- <dig>
ncbi.nlm.nih.gov
ncbi.nlm.nih.gov
ncbi.nlm.nih.gov
myxobacterial genomes are generally large, with a characteristically high %gc.
bdellovibrio bacteriovorus hd <dig> is included as a representative of the bdellovibrios – the closest known relatives of the myxobacteria
athe m. fulvus genome includes a plasmid, which encodes no tcs proteins
this study sought to investigate how the apparent evolutionary mechanisms affecting myxobacterial tcss change with timescale and/or phylogenetic distance.
the analysis was restricted to myxobacterial genomes, as they possess large numbers of tcss, and are available across a broad taxonomic range.
thus comparisons could be performed between organisms within the same species, genus, family, sub-order or order, with representative genomes available for all three sub-orders, and every family of the myxobacteria.
secondary aims of the study were to define myxobacterial tcs gene orthology relationships, and to identify core and accessory tcss.
RESULTS
relatedness of sequenced genomes
to obtain a measure of the divergence of the sequenced myxobacterial genomes, the 16s rrna genes of each genome were aligned and distances  between each pair of sequences calculated, allowing construction of a phylogenetic tree .
there was strong bootstrap support for the tree , which also agreed well with that of garcia et al.,  <cit> .
distances clustered into four distinct groups, which largely shared particular taxonomic ranks .
the ‘genus’ group is composed of distances between strains within the same genus .
comparisons between genomes in the same formally-defined family  were mainly distributed into two groups of genetic distances, here denoted ‘family’ and ‘sub-order’ in fig.
<dig>  the final ‘order’ group was dominated by distances between genomes from the same order, but different sub-orders .
a redrawn phylogenetic tree with distances discretised according to the groupings defined above is shown in fig.
1b.fig.
<dig>
a correlation between 16s rrna gene sequence conservation and taxonomic classification.
pairwise genetic distances between myxococcales 16s rrna gene sequences , plotted in rank order.
distances cluster into four distinct groups:  <dig>  ho/pp vs. sc/ca vs. mx/mf/cc/sa/ak/af/ad1/adc and ho vs. pp ,  <dig>  ak/af/ad1/adc vs. mx/mf/cc/sa ,  <dig>  sc vs. ca, and ak/ad1/adc vs. af, and mx/mf vs. cc vs. sa , and  <dig>  ak vs. ad <dig> vs. adc, and mx vs. mf .
see text for abbreviations.
b phylogenetic tree with distances discretised into four vertical bands as grouped in part a
in several instances, distance-based groupings of genome pairs were not consistent with the formal taxonomy of the myxobacteria.
for instance:  <dig>  bb is only as distant from the myxococcales as pp is from myxococcales of different sub-orders,  <dig>  ho and pp should be assigned to separate sub-orders,  <dig>  af should be in a separate genus to the other anaeromyxobacter,  <dig>  sa belongs to the myxococcaceae family , and  <dig>  the anaeromyxobacter should be a separate family from the other myxococcaceae.
for the rest of this study we have classified inter-genome relationships according to fig.
<dig> rather than using the current taxonomic classification.
thus for our purposes af is considered to be in a different genus  to ak, ad <dig> and adc, while ho and pp are in separate sub-orders.
ak, af, ad <dig> and adc are considered to be in the same sub-order, but in a different family from mf, mx, cc and sa, a suggestion previously made by garcia et al.
<cit> .
tcs protein sets
for each genome, a list of tcs proteins was inferred as described in the methods section.
table  <dig> presents summary statistics for each tcs gene set, alongside myxobacterial mean values.
williams and whitworth  <cit>  provided a characterisation of the tcs genes from  <dig> replicons as a pangenomic yardstick against which to compare tcs gene sets and these are included in table  <dig> as bacterial mean values.table  <dig> two-component gene sets of the myxobacteria
genomea
mx
mf
cc
sa
adc
ad1
af
ak
sc
ca
ho
pp
bn
bacterial mean:b
for each genome, the numbers of tcs genes are presented, broken down by type.
hks are characterised according to the percentages that are hybrid proteins and that contain transmembrane helices.
genomic tcs gene organisation is also presented as the percentage of tcs genes which are orphaned, paired or in complex gene clusters.
values are presented for bdellovibrio bacteriovorus , and the myxobacterial and bacterial means are also shown for comparison
a
ho = haliangium ochraceum dsm  <dig>  pp = plesiocystis pacifica sir- <dig>  sc = sorangium cellulosum so ce  <dig>  ca = chondromyces apiculatus dsm  <dig>  mx = myxococcus xanthus dk <dig>  mf = myxococcus fulvus hw- <dig>  cc = corallococcus coralloides dsm  <dig>  sa = stigmatella aurantiaca dw3/4- <dig>  ak = anaeromyxobacter sp.
k, af = anaeromyxobacter sp.
fw109- <dig>  ad1 = anaeromyxobacter dehalogenans 2cp- <dig>  adc = anaeromyxobacter dehalogenans 2cp-c, and bb = bdellovibrio bacteriovorus hd100
bcalculated from williams and whitworth  <cit>
differences between the myxobacterial and bacterial means are apparent for some metrics.
myxobacteria are characterised by large numbers of tcs genes, greater than expected even for organisms with such large genomes.
figure 2a shows the number of tcs as a function of genome size for  <dig> bacteria and the  <dig> myxobacteria.
myxobacteria also tend to have a high proportion of hybrid hks and a high proportion of tcs genes in complex foci .
they also exhibit exceptionally low proportions of tm hks , implying they sense their internal state to an unusual degree  <cit> .
bb exhibits values intermediate between the myxobacterial and bacterial means, except that  it has a substantially greater proportion of paired tcs genes and tm hks than the average bacterium .fig.
<dig>
a numbers of tcs genes as a function of genome size.
myxobacterial genomes  tend to have larger than expected numbers of tcs genes for genomes of their size.
bdellovibrio bacteriovorus hd <dig>  is also highlighted.
a trendline is shown for all bacteria.
b the proportion of transmembrane histidine kinases , as a function of total hk genes in myxobacterial  and bacterial  genomes.
bdellovibrio bacteriovorus hd <dig>  is also highlighted.
trendlines and their equations are shown for bacteria  and myxobacteria .
updated from  <cit>
on inspection of table  <dig>  some myxobacteria stand out as unusual.
the anaeromyxobacter strains generally exhibit percentages of hybrid kinases more similar to bacteria than myxobacteria, and also possess relatively high percentages of tm hks.
sc, ho and pp have relatively small numbers of tcs genes for myxobacterial genomes of their size.
ho and pp are marine organisms, and marine organisms generally exhibit reduced environmental responsiveness compared to their terrestrial counterparts  <cit> .
sc could be thought of as having an exceptionally large genome for a myxobacterium, rather than having relatively low numbers of tcss for a myxobacterial genome of its size.
orthology relationships between tcs proteins
tcs gene sets were clustered into groups of orthologous and paralagous proteins as described in the methods section.
even at the lowest similarity cut-off used , proteins clustered into pairs from mx and mf were also reciprocal highest-scoring blast hits, supporting the validity of this approach to define orthologues.
the percentage of proteins successfully grouped into clusters of orthologues decreased exponentially with increasing genetic distance .
a notable departure from the trend was the ca and sc comparison, which gave a remarkably low proportion of tcs genes being grouped as orthologues .
this was the only pair-wise comparison which departed substantially from the trend, and the identified orthologues exhibited surprisingly low similarity.
the cause of this extra variability is unclear.fig.
<dig> the percentage of orthologous tcs proteins varies exponentially with genetic distance.
as genetic distance between two strains increases , the percentage of proteins from either genome that are successfully clustered decreases exponentially.
each black diamond represents one pairwise comparison.
comparisons shown include all possible intra-genus and intra-family comparisons , and a representative sample of within sub-order and order comparisons .
a notable departure from the trend is the comparison between ca and sc
within-genus orthologues
we began a detailed analysis of orthologues by comparing the gene set of mf to that of mx, the best characterised myxobacterium.
mx and mf had  <dig> and  <dig> tcs genes respectively which were unique to that organism .
the remaining tcs genes  all clustered into mx/mf pairs except for a single cluster of three proteins , which seems to have arisen by a local duplication in the mf lineage.
a dot-plot of the relative location of tcs orthologues in the two genomes  recapitulated the global synteny observed between mx and mf, which is dominated by a single inversion  <cit> .
there were no observable changes in relative genomic location of orthologues, suggesting that the relocation of individual tcs gene foci around a genome is a rare event, compared to re-location via large-scale chromosomal recombination.
five singleton tcs genes were found close to the putative inversion breakpoints , indicating that gain/loss of those genes was potentially associated with the chromosomal inversion.fig.
<dig> dot-plots of tcs gene orthologues.
orthologues are plotted according to their genomic location  in both genomes.
each genome starts at the origin, with the dnaa gene being the first gene.
a central line of orthologues extending from top left to bottom right indicates an inversion of a large section of the chromosome, centered around the origin/terminus.
a dot-plot of mx/mf tcs gene orthologues.
b dot-plot of adc/ak orthologues
within mx/mf paired clusters, tcs character  was conserved in the vast majority  of cases.
in four pairs gene organisation was altered /lilab_ <dig> and mxan_5147/lilab_33705), in two cases domain architecture/family was altered , and in one case both gene organisation and family membership changed between the two organisms .
of the five observed changes in organisation, two were caused by gain/loss of intervening genes at a single focus, while the remaining three were caused by gain/loss of tcs genes.
all three examples of family membership changes involved gain/loss of rr effector domains, and were associated with frameshift mutations.
intriguingly all three pairs included genes which were either mispredicted or annotated as pseudogenes, so it isn’t clear whether the observed differences represent true changes in family membership, or gene loss in progress.
there were seven cases of altered gene length  between mx and mf pairs, four of which involved gain/loss of effector/input domains.
two were changes of family membership mentioned previously , one was gain/loss of a duplicated pas domain , and there was one instance of expansion/contraction of a poly-hamp domain tract .
the remaining three changes in gene length involved gain/loss of a c-terminal region /lilab_30440), extension/contraction of an n-terminal region , and gain/loss of an internal region , all of which may represent the gain/loss of currently uncharacterised domains.
across the whole mx/mf genome, 23 % of protein-encoding genes were unique to either genome, whereas only 10–11 % of tcs genes were singletons, implying that tcs genes are a relatively conserved gene set.
mx tcs genes with names have published phenotypes, and historically mx has been subjected to exhaustive genetic screens for developmental and motility mutations.
the set of named tcs genes is therefore enriched for genes with significant phenotypes .
mx/mf paired tcs genes were enriched for named genes –24 % of all m. xanthus tcs genes are named, however only two of the  <dig> singletons were named , and those two genes  have minor phenotypes  <cit> , suggesting that singletons are not subject to a strong retentive selective pressure.
the mx/mf singletons were very different in character - mx encoded eight hybrid kinases,  <dig> tcs genes with complex gene organisations and  <dig> paired tcs genes, whereas the mf singletons included no hybrid kinases, only three complex tcs genes, but  <dig> paired genes.
to see whether mx/mf differences were typical of species-level comparisons, anaeromyxobacter strains ak and adc, were also compared.
the adc/ak genetic distance  is identical to that between mx and mf, however the anaeromyxobacter species are anaerobic, and non-fruiting, so any commonality between the ak/adc and mx/mf comparisons is likely to be due to genetic similarity and independent of lifestyle.
no duplications were observed, while 9 % of genes were singletons .
six changes in gene organisation were observed , one change in size and domain architecture  and one other gene size change .
a dot-plot of ak/adc orthologues  shows global synteny, including across a large chromosomal inversion, albeit with three pairs of orthologues demonstrating re-location of tcs genes within the chromosome , including one causing a change in gene organisation .
insertions/deletions  of ≥ <dig> gene at three foci , resulted in four singletons, and the remaining five changes in gene organisation.
one small  indel  resulted in one singleton, and a domain architecture change and gene size change, as the consequences of changing between a tr,t system and a tt system .
finally, an indel resulted in an apparent gain/loss of two input domains from one hk .
unlike with the mx/mf singletons, there were no obvious differences in character between the adc and ak singletons.
in summary , these two within-genus comparisons , suggest that inter-species differences are dominated by tcs gene gain/loss, occasionally causing changes in gene organisation.
this also leads to ‘individuality’ within each genome, with organism-specific ‘accessory’ genes.
relatively infrequently, changes in domain architecture occur as a consequence of small indels or frameshift mutations.
as expected, such alterations most frequently affect the output domains of rrs  and the sensory domains of hks.table  <dig> mutational events and their consequences for tcs gene sets
percentages of orthologous clusters which show evidence of evolutionary changes when comparing genomes from the same genus  and the same family .
percentages obtained from pairwise comparisons between cystobacterineae members
within-family orthologues
to assess tcs gene changes between genera, three ‘within-family’ groups of comparisons were made; cc/sa/mx/mf, ad1/adc/ak/af, and sc/ca.
in cc/sa/mx/mf comparisons, some clusters of homologues were observed with > <dig> members.
one cluster of six proteins was observed due to an apparent double-duplication in sa , while single duplications in cc and mf gave rise to two clusters of five proteins .
a final cluster of five proteins was composed entirely of sa proteins .
<dig> proteins were found in clusters of four , all but two of which included a single representative from each genome .
another  <dig> proteins were found in groups of three, with members from different genomes, however in two cases there was an apparent duplication in one of the genomes .
<dig> proteins were paired, including five pairs of duplicated proteins from sa .
this left  <dig> singletons, of which the majority belonged to sa , with large numbers also in cc , and relatively few in mx  or mf .
the pattern of protein orthology can be seen in additional file 2: figure s <dig> as a venn diagram, and numbers are consistent with genetic distance, with closer genomes sharing more orthologues and possessing fewer singletons.
the same trends were observed for ad1/adc/ak/af orthologues .
across the  <dig> cc/sa/mx/mf clusters with two or more members, in addition to the mx/mf changes described earlier, there were  <dig> cases of altered gene organisation , and three changes of domain architecture, all three of which involved hybrid kinases.
a small indel caused a rrt/rt change , a trt/trttrt change appears associated with a tandem duplication and fusion , while a tr/t change is associated with a large indel, also involving gain/loss of an adjacent gene .
in addition to these changes in domain architecture,  <dig> cases of genes changing size by > <dig> aa were observed, of which the majority  affected hks, mainly in their n-terminal regions .
for ad1/adc/ak/af orthologues ,  <dig> instances of changes in gene organization were observed .
these were mainly  associated with large  indels, with five cases also involving gene rearrangements , and two of the  <dig> were due to frameshift mutations generating pseudogenes .
four changes in domain architecture were observed, each of which was associated with a change in gene size.
in one case an n-terminal pilz domain was gained/lost from a rr , while the other changes were due to gene fusion/fission events .
seven gene duplications at five loci were apparent , with an additional duplication also associated with a change in gene size .
finally, two additional changes in gene size were seen, both of which were due to gain/loss of c-terminal regions .
the final within-family comparison possible is between ca and sc, however that comparison is characterised by an unusually low proportion of orthologous proteins  and low similarity between orthologues.
of the  <dig> clusters of orthologues which included members from both organisms, there were no apparent changes in domain architecture, gene organisation or gene size.
further comparisons between more distant genomes were also uninformative due to the small number of orthologous protein clusters identified by cd-hit or by bidirectional best-scoring blast hits.
thus, in summary, comparisons between same-family members show similar trends to those from same-genus comparisons albeit at higher frequencies .
complexity
another notable feature of myxobacterial tcs proteins is their complexity – both in terms of their domain architectures, and their gene organisation  <cit> .
we define here the most complex systems as ‘intricate’, where at least two receiver and at least two transmitter domains are found in the same protein and/or genetic focus.
trtr proteins are encoded as singletons in four genomes , while sa also encodes a trttrt protein .
according to our analysis, the trtr proteins are not orthologous to each other or other myxobacterial tcs proteins, while the trttrt protein appears to be a tandem duplication and fusion of two adjacent trt proteins.
in addition, myxobacterial genomes possessed between two and nine intricate multi-gene tcs foci , considerably more than the bacterial mean.
the intricate foci of ho, sc and af were largely unique to those genomes, otherwise orthologous relationships between the intricate foci were relatively clear at the genus and family level .
the foci highlighted exhibit several interesting properties.
there seems to be a great deal of plasticity in these foci with many exhibiting orthology with non-intricate tcs foci as a consequence of gene gain/loss and/or recombination .
in addition, a large proportion of the tcs domains encoded at intricate foci are found in hybrid kinases .
the proportion of foci containing large numbers of tcs domains are similar between bacteria and myxobacteria .
nevertheless, when considering all tcs proteins , those of myxobacteria are generally more complex than other bacteria.
table  <dig> shows the relative frequencies of tcs proteins with different domain compositions for myxobacteria and bacteria.
myxobacteria possess far greater proportions of tcs proteins with two or more transmitter domains, and hybrid kinases with three or more receiver domains than other bacteria.table  <dig> transmitter and receiver domain composition of tcs proteins
values are the % of myxobacterial tcs proteins exhibiting the number of receiver  and transmitter  domains in question divided by the % of all such bacterial tcs proteins  <cit>
athere are no such proteins in bacteria
bno such proteins in myxobacteria .
numbers of such myxobacterial proteins are shown in parentheses
tcs gene expansion
nearly two thirds  of all mx/mf singletons were found to have orthologues in cc and/or sa, implying that the formation of the majority of singletons was through progressive gene loss.
however, genes can also be gained by duplication and divergence.
yet duplication is responsible for a relatively small amount of differences between tcs gene sets.
for instance, clustering reveals only  <dig> lineage-specific duplications in sa, four cases in cc and two in mf.
eight duplications were apparent within the anaeromyxobacter strains .
nevertheless, it appears that certain tcs families have been expanded hugely in the myxobacteria and this occurred relatively recently.
for instance, ntrc family members are unusually numerous  <cit> .
indeed, eight of the  <dig> mx ntrc family members exhibited more similarity to other mx ntrc paralogues than to homologues in other organisms, suggesting origin by duplication after the divergence of contemporary genera of myxobacteria.
another mechanism that can increase gene number is horizontal gene transfer .
to test the frequency of horizontal acquisition, the mx/mf singletons lacking homologues in cc/sa/mx/mf were queried against the nr protein database using blast.
eight of the ten mf singletons and six of the  <dig> mx singletons gave blast hits across the myxococcales, implying linear descent.
however nine of the  <dig> mx/mf genes gave best non-self blast hits to non-myxococcales genomes, implying horizontal acquisition.
such evidence of hgt is also abundant in the singletons of other myxobacteria ).
methods
phylogenetic analysis
16s rrna gene sequences were aligned using clustalw  <cit> , using default parameter settings, and neighbour-joining trees were created using phylip  <cit>  with  <dig> bootstraps, excluding columns containing gaps.
inference of tcs proteins sets
eight myxobacterial genomes  and that of bb are included within the latest version of the p2cs  database, available at www.p2cs.org  <cit> , which provides an automated and consistent pipeline for the identification and annotation of tcs genes.
the same categorisation pipeline is available for unpublished sequences at www.p2rp.org  <cit>  and this service was used to characterise the tcs gene sets for sa, cc, pp and ca nucleotide sequence files from www.ncbi.nlm.nih.gov.
tcs gene sets were compared with those previously published for mx, sa, adc and sc  <cit> , and discrepancies manually curated.
for the tcs gene sets available through p2cs, potentially ‘mispredicted’ tcss were identified, and these were validated manually and added to the tcs datasets .
these two processes resulted in the addition of five, six and  <dig> tcs proteins to the gene sets for mx, sc and adc respectively, as originally described by whitworth and cock  <cit> .
three, one, one and three ‘mispredicted’ tcs genes were also validated for mf, ho, ad <dig> and af respectively.
tcs datasets for each genome are available as additional data .
characterisation of tcs proteins
domain architectures and the presence of transmembrane  helices within tcs proteins were determined as described by barakat et al.
<cit> .
tcs gene organisation was defined using a classification scheme based on tcs domains  encoded unidirectionally at a ‘focus’, as proposed by williams and whitworth  <cit> .
if a transmitter and a receiver domain were found adjacent in the genome , and they had at least two non-tcs genes on each side, then the relevant genes were described as ‘paired’ .
if a gene contained a single transmitter, receiver or phosphotransfer domain and was found more than two genes from another tcs gene, it was described as ‘orphan’.
tcs gene foci with any other domain/gene organisation were described as ‘complex’.
the genomes of pp and ca are incomplete, so gene organisations were not defined for those organisms.
focus tcs domain organisation is described in the text using r, t and h to represent receiver, transmitter and hiska/hpt domains respectively, with commas separating genes.
thus a ‘rt,r’ organisation represents a hybrid kinase gene encoding an n-terminal receiver domain and a c-terminal transmitter domain, adjacent to a response regulator gene encoding a single receiver domain.
sequence-based clustering
to define orthologues and paralogues, sets of tcs genes were clustered using cd-hit  <cit> .
a high similarity cut-off was used  to give high confidence clusters, and then a lower similarity  cut-off was used to identify more divergent clusters.
clusters obtained were then subjected to curation, whereby any large clusters containing multiple paralogues obtained using the 30 % cut-off were manually split into the high-confidence clusters obtained with a 90 % cut-off.
clusters which included proteins of significantly different length , or of varying domain architecture, family membership and/or gene organisation, or with varying numbers of members per genome were earmarked for further analysis.
discussion
genetic distance and inferred evolutionary events
comparing the tcss from organisms with differing degrees of relatedness allowed inference of the contribution of different types of mutation to the evolution of contemporary sets of tcss.
it appears that the relative impact of different mutational events are not affected by the relatedness of the organisms being compared, however, the apparent frequency of mutations is dependent on the relatedness of compared genomes , echoing the findings of he et al.
<cit>  in their study of clostridia.
this means that when comparing the gene sets of multiple genomes, it is important to characterize the similarity of those genomes and compensate for differing degree of relatedness.
it is even worth including organisms of diverse relatedness, to increase the power of analysis, as suggested by eisen et al.,  <cit> .
in this study we used 16s rrna gene sequence similarity as a straightforward indicator of genomic relatedness/distance and found that the ability to identify tcs orthologues between genomes was dependent on relatedness .
at low levels of relatedness, difficulties in the identification of orthologues reduced the power of comparative analysis making it difficult to infer evolutionary events.
this difficulty was likely exacerbated by the large numbers of tcs in myxobacterial genomes, and should be less of a problem for the majority of prokaryotes.
there have been several studies comparing the tcs gene sets of related organisms, and although not noted, the phenomena we describe can often be seen in the tcs gene sets of diverse taxa.
for instance, tcs gene complements can vary by 10–30 % when comparing organisms within the same family/genus, and the frequency of within-species variation is about half that observed when comparing organisms within the same genus .
evolution of myxobacterial tcs gene sets
our data suggest that most changes in myxobacterial tcs gene sets are due to large indels, some of which are associated with chromosomal recombination.
there is evidence of frequent gene loss , but also gene gain via hgt and occasional gene duplications.
some  of these indels result in changes in gene organization, when the indel only affects part of a tcs locus.
small  indels are less frequent and often result in domain gain/loss with associated changes in tcs family membership.
other studies have reported broadly consistent results.
for example, in a study of recombination rates in salmonella, sun et al.,  <cit>  found that deletions are significantly more frequent than duplications.
both mx and bb exhibit genome-wide evidence of hgt  <cit> , and we find that hk and rr genes tend to be acquired/lost as pairs, which has been observed in several other taxa, including mutans streptococci and xanthomonas spp.
<cit> .
domain shuffling appears to have happened occasionally during tcs evolution, giving rise to hks with a multiplicity of input domains, and rrs with diverse output domains  <cit> .
our methods of defining orthologues used whole gene sequences and may not have been capable of identifying shuffled domains, nevertheless we did observe several instances of input/output domain gain/loss.
with domain shuffling being a rare evolutionary event, it seems probable that apparent domain shuffling is a consequence of domain loss followed by domain gain rather than direct shuffling between proteins.
we also did not observe significant relocation of tcs within genomes.
the few cases observed  were associated with breakpoints of chromosomal recombination.
it is possible that those chromosomal relocations occurred due to recombination between homologous tcs genes at those points, however there is too little conservation of synteny in the break point regions to provide any evidence of this happening in myxobacterial genomes.
myxobacteria were chosen as the subject of this study due, in part, to the large number of tcs genes they possess.
potentially, this could limit extrapolation of our conclusions to other organisms, as increased numbers of tcs genes may have promoted some mutagenic events , and the large numbers of paralogous genes makes robust orthologue identification difficult.
in  <dig> we described the tcs of four sequenced myxobacteria and used domain-based orthology searches to identify likely evolutionary events and peculiarities of myxobacteial tcs gene sets  <cit> .
for this study we added a further eight myxobacterial genomes, dramatically increasing the number of inter-genome comparisons possible, and defined orthologues using whole gene sequence-based methods.
the results obtained here agree well with that previous work, and provide a less anecdotal analysis of evolutionary events.
benchmarking our genome comparisons by relatedness also allows straightforward comparison with tcs gene sets from other genomes in the future.
the tcs genes of myxobacteria seem to be relatively conserved compared to non-tcs genes, and we also found that the mx tcs genes with mf orthologues were enriched for named genes, which tend to have observable phenotypes.
this suggests that there is a ‘core’ set of tcs genes in myxobacteria responsible for regulating important functions, with other ‘accessory’ tcs genes being more variable between genomes, and with minor  observable phenotypes.
the ad1/adc/af/ak tcs gene sets exhibited much less variability than those of mx/mf/cc/sa , however this effect was proportional with the reduced numbers of tcs genes in the anaeromyxobacteria.
it is possible that the fruiting myxobacteria need enhanced ‘individuality’ to allow easy discrimination of competing genotypes, but whether this is a selected trait, or purely a consequence of their increased numbers of tcs is uncertain.
among strain-specific genes, there is particularly good evidence of acquisition by hgt, which is known to have had a major effect on the shaping of contemporary myxobacterial genomes, alongside lineage-specific expansion  <cit> .
myxobacterial genomics
the field of myxobacterial genomics is currently hampered by the large number of genes encoded in each genome , the lack of easily assayable behavioural outputs , and pleiotropy of many genes.
this means that phenotypic properties cannot easily be correlated with particular genomic differences.
however, further myxobacterial genomes are becoming available at an accelerating rate, and it is possible that in the near future, physiological characterisation of the new genome host organisms will enable correlations with genomic features.
to that end, it is important to develop benchmarked catalogs of genomic differences.
taking into consideration the relative distance of compared genomes is an important feature of benchmarking, and this study provides a framework to develop such benchmarks using 16s ribosomal rna gene sequence differences as indicators of genetic distance, which can be related to formal taxonomy and other systematic methods  <cit> .
during the course of this study draft genomes became available for strains of cystobacter fuscus, cystobacter violaceus, enhygromyxa salina, hyalangium minutum, sandaracinus amylolyticus, and for two further strains of myxococcus xanthus, including dz <dig> and dzf <dig>  <cit> , however these genomes are not complete or scheduled for finishing.
thankfully, two complete myxobacterial genomes have also became available, myxococcus stipitatus dsm  <dig> and sorangium cellulosum so0157- <dig>  <cit> .
both of these genomes are now characterized in p2cs and seem typical members of their genera, albeit with slightly more tcs genes , except that the lyttr family of rrs is expanded substantially in ms .
as more genomes are completed it will become easier to identify more distant orthologues via ‘gap-bridging’ genomes.
this will in turn improve deduction of the directionality of evolutionary changes such as indels, and allow us to create a more finely resolved timeline of myxobacterial tcs gene evolution.
CONCLUSIONS
focusing on the tcs genes of myxobacteria, this study investigated the relationship between apparent evolutionary events, and genomic distance or ‘relatedness’.
as genomic distance increases it becomes more difficult to identify tcs orthologues, and the relative frequency of evolutionary events increases.
however, the relative frequency of evolutionary changes remain unaffected by scale , and it seems that in the myxobacteria, as for other organisms, tcs gene evolution is dominated by large indels , with smaller indels being less frequent .
tcs genes are ‘born’ by duplication and acquisition through hgt, but the individuality of tcs gene sets appears to be dominated by lineage-specific gene loss.
more generally, our results suggest that researchers should be mindful that the numbers of genomes compared can have a substantial effect on the apparent frequency of observed events.
this is particularly true when comparing organisms that share different levels of ‘genetic relatedness’, which should not be assumed to always be consistent with the taxonomic assignments of compared organisms.
additional files
additional file 1: figure s <dig>  genetic distance  and neighbour-joining tree inferred from myxobacterial 16s rrna gene sequences.
bootstrap support for each clade was > 95 %  in all cases except for the clade indicated by a filled circle .
the tree is rooted with b. bacteriovorus.
the bar represents  <dig>  substitutions per nucleotide position.
additional file 2: figure s <dig>  orthology relationships within-family.
venn diagrams showing the distribution of orthologous genes between organisms.
some orthologous clusters exhibited multiple members from the same organism, suggestive of lineage-specific gene duplications.
for the purpose of these diagrams, duplicated proteins were treated as singletons.
for instance, if a cluster contained members from sa, cc, mx and mf, but the sa protein was duplicated, then  <dig> was entered in the sa/cc/mx/mf cell, and  <dig> into the sa only cell.
for a pair of sa proteins,  <dig> was entered into the sa cell.
a: sa/cc/mx/mf, b: af/ak/ad1/adc.
additional file 3:
the sequence clusters of orthologues of every myxobacterial tcs protein.
additional file 4: table s <dig>  intricate multi-gene tcs foci.
each row shows all myxobacterial orthologues of intricate multi-gene foci, including non-intricate orthologues, arranged by organism.
foci in particular genomes are denoted by the number from the locus tag of the first protein within the focus.
‘//’ refers to an unspecified number  of non-tcs genes.
additional file 5: table s <dig>  evidence of horizontal gene transfer  amongst cystobacterial tcs proteins.
singletons  from comparisons between ad1/adc/ak/af, and between mx/mf/cc/sa, are presented alongside the number of singletons acquired by hgt .
also tabulated are the number of tcs genes found in foci where every tcs gene was acquired by putative hgt and the number of those which constitute an ‘intact’ tcs
additional file 6:
datasets of the annotated tcs proteins from each genome.
abbreviations
tcstwo-component system
hkhistidine kinase
rrresponse regulator
tmtransmembrane
competing interests
the author declares that he has no competing interests.

 two-component systems are the dominant phosphorylation based signalling pathways in prokaryotes. the average eubacterial genome encodes ~ tcs proteins, and they regulate diverse and important behaviours, including pathcomparative genomics and experimentalogenesis, cell cycle progression, sporulation and drug resistance. typical tcss consist of two proteins – a histidine kinase, which can sense an environmental parameter through an ‘input’ domain, and a response regulator, which pocomparative genomics and experimentalssesses an effector function via an ‘output’ domain. in response to an environmental signal the hk autophosphorylates within its ‘transmitter’ domain and the phosphoryl group is then transferred onto the ‘receiver’ domain of the partner rr, modulating its output domain effector activity. common modifications to the basic scheme include hybrid kinases – proteins which contain both transmitter and receiver domains, and phosphotransfer proteins - which can shuttle phosphoryl groups between receiver domains forming a phosphorelay. for comprehensive reviews of tcs biology, refer to the recent book by gross and beier and the special issue of current opinion in microbiology edited by bourret and silversmith. comparative genomics and experimental approaches have shed light on the general principles of tcs evolution. for instance tcss have evolved primarily by duplication and divergence of ancestral systems, with the contemporary diversity of tcs architecture generated by domain shuffling. this is important, as an understanding of natural variation due to evolution can inform the rational engineering of tcss for biotechnological and research applications. nevertheless there are important questions left unanswered regarding tcs evolution and most studies have been forced to make several assumptions in their analyses. for instance, inferences regarding tcs evolution have been based on the comparison of various extant genomes. common approaches are to use all available genomes, or to use multiple organisms from a single taxonomic unit. in both cases this typically results in comparisons between genomes which exhibit different relatedness to one another. thus relatively distant and close genomes are considered equally, potentially distorting the conclusions drawn. this is particularly problematic for studies of tcs evolution, as it is thought that different evolutionary processes work at different timescales. for example, shuffling of hk input domains seems to occur rapidly post gene duplication, and gene gain/loss apparently dominates over gene fusion/fission at long timescales. in addition, studies are also often limited in power as there is sequencing bias within/between taxa, most organisms contain just dozens of tcs proteins, and at most taxonomic levels there is a relative shortfall in the number of available genomes, all of which can distort conclusions. myxobacteria are well known for their large numbers of tcss, many of which have been characterised experimentally. they exhibit complex developmental and predatory behaviours and are of ongoing importance to the biotechnology industry. a major focus of myxobacterial research has been elucidation of the molecular mechanisms underpinning multicellular development, which are surprisingly plastic, and involve large numbers of tcs proteins. at the time of analysis, genome sequences were available for myxobacteria. for ease of reference, genomes are abbreviated in the remaining text as shown in table two of these genome sequences are incomplete, but were included in the study anyway, as they were nearly complete, being composed of relatively small numbers of contigs. bdellovibrio bacteriovorus hd was also included in the study as the bdellovibrionales are the closest known relatives of the myxobacteria.table order myxococcales – genome properties, taxonomy and ecology myxococcus fulvus hw- ncbi.nlm.nih.gov ncbi.nlm.nih.gov ncbi.nlm.nih.gov myxobacterial genomes are generally large, with a characteristically high %gc. bdellovibrio bacteriovorus hd is included as a representative of the bdellovibrios – the closest known relatives of the myxobacteria athe m. fulvus genome includes a plasmid, which encodes no tcs proteins this study sought to investigate how the apparent evolutionary mechanisms affecting myxobacterial tcss change with timescale and/or phylogenetic distance. the analysis was restricted to myxobacterial genomes, as they possess large numbers of tcss, and are available across a broad taxonomic range. thus comparisons could be performed between organisms within the same species, genus, family, sub-order or order, with representative genomes available for all three sub-orders, and every family of the myxobacteria. secondary aims of the study were to define myxobacterial tcs gene orthology relationships, and to identify core and accessory tcss.  relatedness of sequenced genomes to obtain a measure of the divergence of the sequenced myxobacterial genomes, the 16s rrna genes of each genome were aligned and distances between each pair of sequences calculated, allowing construction of a phylogenetic tree. there was strong bootstrap support for the tree, which also agreed well with that of garcia et al.,. distances clustered into four distinct groups, which largely shared particular taxonomic ranks. the ‘genus’ group is composed of distances between strains within the same genus. comparisons between genomes in the same formally-defined family were mainly distributed into two groups of genetic distances, here denoted ‘family’ and ‘sub-order’ in fig.  the final ‘order’ group was dominated by distances between genomes from the same order, but different sub-orders. a redrawn phylogenetic tree with distances discretised according to the groupings defined above is shown in fig. 1b.fig. a correlation between 16s rrna gene sequence conservation and taxonomic classification. pairwise genetic distances between myxococcales 16s rrna gene sequences, plotted in rank order. distances cluster into four distinct groups: ho/pp vs. sc/ca vs. mx/mf/cc/sa/ak/af/ad1/adc and ho vs. pp, ak/af/ad1/adc vs. mx/mf/cc/sa, sc vs. ca, and ak/ad1/adc vs. af, and mx/mf vs. cc vs. sa, and ak vs. ad vs. adc, and mx vs. mf. see text for abbreviations. b phylogenetic tree with distances discretised into four vertical bands as grouped in part a in several instances, distance-based groupings of genome pairs were not consistent with the formal taxonomy of the myxobacteria. for instance: bb is only as distant from the myxococcales as pp is from myxococcales of different sub-orders, ho and pp should be assigned to separate sub-orders, af should be in a separate genus to the other anaeromyxobacter, sa belongs to the myxococcaceae family, and the anaeromyxobacter should be a separate family from the other myxococcaceae. for the rest of this study we have classified inter-genome relationships according to fig.  rather than using the current taxonomic classification. thus for our purposes af is considered to be in a different genus to ak, ad and adc, while ho and pp are in separate sub-orders. ak, af, ad and adc are considered to be in the same sub-order, but in a different family from mf, mx, cc and sa, a suggestion previously made by garcia et al.. tcs protein sets for each genome, a list of tcs proteins was inferred as described in the methods section. table  presents summary statistics for each tcs gene set, alongside myxobacterial mean values. williams and whitworth provided a characterisation of the tcs genes from replicons as a pangenomic yardstick against which to compare tcs gene sets and these are included in table  as bacterial mean values.table two-component gene sets of the myxobacteria genomea mx mf cc sa adc ad1 af ak sc ca ho pp bn bacterial mean:b for each genome, the numbers of tcs genes are presented, broken down by type. hks are characterised according to the percentages that are hybrid proteins and that contain transmembrane helices. genomic tcs gene organisation is also presented as the percentage of tcs genes which are orphaned, paired or in complex gene clusters. values are presented for bdellovibrio bacteriovorus, and the myxobacterial and bacterial means are also shown for comparison a ho = haliangium ochraceum dsm pp = plesiocystis pacifica sir- sc = sorangium cellulosum so ce ca = chondromyces apiculatus dsm mx = myxococcus xanthus dk mf = myxococcus fulvus hw- cc = corallococcus coralloides dsm sa = stigmatella aurantiaca dw3/4- ak = anaeromyxobacter sp. k, af = anaeromyxobacter sp. fw109- ad1 = anaeromyxobacter dehalogenans 2cp- adc = anaeromyxobacter dehalogenans 2cp-c, and bb = bdellovibrio bacteriovorus hd100 bcalculated from williams and whitworth differences between the myxobacterial and bacterial means are apparent for some metrics. myxobacteria are characterised by large numbers of tcs genes, greater than expected even for organisms with such large genomes. figure 2a shows the number of tcs as a function of genome size for bacteria and the myxobacteria. myxobacteria also tend to have a high proportion of hybrid hks and a high proportion of tcs genes in complex foci. they also exhibit exceptionally low proportions of tm hks, implying they sense their internal state to an unusual degree. bb exhibits values intermediate between the myxobacterial and bacterial means, except that it has a substantially greater proportion of paired tcs genes and tm hks than the average bacterium.fig. a numbers of tcs genes as a function of genome size. myxobacterial genomes tend to have larger than expected numbers of tcs genes for genomes of their size. bdellovibrio bacteriovorus hd is also highlighted. a trendline is shown for all bacteria. b the proportion of transmembrane histidine kinases, as a function of total hk genes in myxobacterial and bacterial genomes. bdellovibrio bacteriovorus hd is also highlighted. trendlines and their equations are shown for bacteria and myxobacteria. updated from on inspection of table  some myxobacteria stand out as unusual. the anaeromyxobacter strains generally exhibit percentages of hybrid kinases more similar to bacteria than myxobacteria, and also possess relatively high percentages of tm hks. sc, ho and pp have relatively small numbers of tcs genes for myxobacterial genomes of their size. ho and pp are marine organisms, and marine organisms generally exhibit reduced environmental responsiveness compared to their terrestrial counterparts. sc could be thought of as having an exceptionally large genome for a myxobacterium, rather than having relatively low numbers of tcss for a myxobacterial genome of its size. orthology relationships between tcs proteins tcs gene sets were clustered into groups of orthologous and paralagous proteins as described in the methods section. even at the lowest similarity cut-off used, proteins clustered into pairs from mx and mf were also reciprocal highest-scoring blast hits, supporting the validity of this approach to define orthologues. the percentage of proteins successfully grouped into clusters of orthologues decreased exponentially with increasing genetic distance. a notable departure from the trend was the ca and sc comparison, which gave a remarkably low proportion of tcs genes being grouped as orthologues. this was the only pair-wise comparison which departed substantially from the trend, and the identified orthologues exhibited surprisingly low similarity. the cause of this extra variability is unclear.fig. the percentage of orthologous tcs proteins varies exponentially with genetic distance. as genetic distance between two strains increases, the percentage of proteins from either genome that are successfully clustered decreases exponentially. each black diamond represents one pairwise comparison. comparisons shown include all possible intra-genus and intra-family comparisons, and a representative sample of within sub-order and order comparisons. a notable departure from the trend is the comparison between ca and sc within-genus orthologues we began a detailed analysis of orthologues by comparing the gene set of mf to that of mx, the best characterised myxobacterium. mx and mf had and tcs genes respectively which were unique to that organism. the remaining tcs genes all clustered into mx/mf pairs except for a single cluster of three proteins, which seems to have arisen by a local duplication in the mf lineage. a dot-plot of the relative location of tcs orthologues in the two genomes recapitulated the global synteny observed between mx and mf, which is dominated by a single inversion. there were no observable changes in relative genomic location of orthologues, suggesting that the relocation of individual tcs gene foci around a genome is a rare event, compared to re-location via large-scale chromosomal recombination. five singleton tcs genes were found close to the putative inversion breakpoints, indicating that gain/loss of those genes was potentially associated with the chromosomal inversion.fig. dot-plots of tcs gene orthologues. orthologues are plotted according to their genomic location in both genomes. each genome starts at the origin, with the dnaa gene being the first gene. a central line of orthologues extending from top left to bottom right indicates an inversion of a large section of the chromosome, centered around the origin/terminus. a dot-plot of mx/mf tcs gene orthologues. b dot-plot of adc/ak orthologues within mx/mf paired clusters, tcs character was conserved in the vast majority of cases. in four pairs gene organisation was altered /lilab_ and mxan_5147/lilab_33705), in two cases domain architecture/family was altered, and in one case both gene organisation and family membership changed between the two organisms. of the five observed changes in organisation, two were caused by gain/loss of intervening genes at a single focus, while the remaining three were caused by gain/loss of tcs genes. all three examples of family membership changes involved gain/loss of rr effector domains, and were associated with frameshift mutations. intriguingly all three pairs included genes which were either mispredicted or annotated as pseudogenes, so it isn’t clear whether the observed differences represent true changes in family membership, or gene loss in progress. there were seven cases of altered gene length between mx and mf pairs, four of which involved gain/loss of effector/input domains. two were changes of family membership mentioned previously, one was gain/loss of a duplicated pas domain, and there was one instance of expansion/contraction of a poly-hamp domain tract. the remaining three changes in gene length involved gain/loss of a c-terminal region /lilab_30440), extension/contraction of an n-terminal region, and gain/loss of an internal region, all of which may represent the gain/loss of currently uncharacterised domains. across the whole mx/mf genome, 23 % of protein-encoding genes were unique to either genome, whereas only 10–11 % of tcs genes were singletons, implying that tcs genes are a relatively conserved gene set. mx tcs genes with names have published phenotypes, and historically mx has been subjected to exhaustive genetic screens for developmental and motility mutations. the set of named tcs genes is therefore enriched for genes with significant phenotypes. mx/mf paired tcs genes were enriched for named genes –24 % of all m. xanthus tcs genes are named, however only two of the singletons were named, and those two genes have minor phenotypes, suggesting that singletons are not subject to a strong retentive selective pressure. the mx/mf singletons were very different in character - mx encoded eight hybrid kinases, tcs genes with complex gene organisations and paired tcs genes, whereas the mf singletons included no hybrid kinases, only three complex tcs genes, but paired genes. to see whether mx/mf differences were typical of species-level comparisons, anaeromyxobacter strains ak and adc, were also compared. the adc/ak genetic distance is identical to that between mx and mf, however the anaeromyxobacter species are anaerobic, and non-fruiting, so any commonality between the ak/adc and mx/mf comparisons is likely to be due to genetic similarity and independent of lifestyle. no duplications were observed, while 9 % of genes were singletons. six changes in gene organisation were observed, one change in size and domain architecture and one other gene size change. a dot-plot of ak/adc orthologues shows global synteny, including across a large chromosomal inversion, albeit with three pairs of orthologues demonstrating re-location of tcs genes within the chromosome, including one causing a change in gene organisation. insertions/deletions of ≥ gene at three foci, resulted in four singletons, and the remaining five changes in gene organisation. one small indel resulted in one singleton, and a domain architecture change and gene size change, as the consequences of changing between a tr,t system and a tt system. finally, an indel resulted in an apparent gain/loss of two input domains from one hk. unlike with the mx/mf singletons, there were no obvious differences in character between the adc and ak singletons. in summary, these two within-genus comparisons, suggest that inter-species differences are dominated by tcs gene gain/loss, occasionally causing changes in gene organisation. this also leads to ‘individuality’ within each genome, with organism-specific ‘accessory’ genes. relatively infrequently, changes in domain architecture occur as a consequence of small indels or frameshift mutations. as expected, such alterations most frequently affect the output domains of rrs and the sensory domains of hks.table mutational events and their consequences for tcs gene sets percentages of orthologous clusters which show evidence of evolutionary changes when comparing genomes from the same genus and the same family. percentages obtained from pairwise comparisons between cystobacterineae members within-family orthologues to assess tcs gene changes between genera, three ‘within-family’ groups of comparisons were made; cc/sa/mx/mf, ad1/adc/ak/af, and sc/ca. in cc/sa/mx/mf comparisons, some clusters of homologues were observed with > members. one cluster of six proteins was observed due to an apparent double-duplication in sa, while single duplications in cc and mf gave rise to two clusters of five proteins. a final cluster of five proteins was composed entirely of sa proteins. proteins were found in clusters of four, all but two of which included a single representative from each genome. another proteins were found in groups of three, with members from different genomes, however in two cases there was an apparent duplication in one of the genomes. proteins were paired, including five pairs of duplicated proteins from sa. this left singletons, of which the majority belonged to sa, with large numbers also in cc, and relatively few in mx or mf. the pattern of protein orthology can be seen in additional file 2: figure s as a venn diagram, and numbers are consistent with genetic distance, with closer genomes sharing more orthologues and possessing fewer singletons. the same trends were observed for ad1/adc/ak/af orthologues. across the cc/sa/mx/mf clusters with two or more members, in addition to the mx/mf changes described earlier, there were cases of altered gene organisation, and three changes of domain architecture, all three of which involved hybrid kinases. a small indel caused a rrt/rt change, a trt/trttrt change appears associated with a tandem duplication and fusion, while a tr/t change is associated with a large indel, also involving gain/loss of an adjacent gene. in addition to these changes in domain architecture, cases of genes changing size by > aa were observed, of which the majority affected hks, mainly in their n-terminal regions. for ad1/adc/ak/af orthologues, instances of changes in gene organization were observed. these were mainly associated with large indels, with five cases also involving gene rearrangements, and two of the were due to frameshift mutations generating pseudogenes. four changes in domain architecture were observed, each of which was associated with a change in gene size. in one case an n-terminal pilz domain was gained/lost from a rr, while the other changes were due to gene fusion/fission events. seven gene duplications at five loci were apparent, with an additional duplication also associated with a change in gene size. finally, two additional changes in gene size were seen, both of which were due to gain/loss of c-terminal regions. the final within-family comparison possible is between ca and sc, however that comparison is characterised by an unusually low proportion of orthologous proteins and low similarity between orthologues. of the clusters of orthologues which included members from both organisms, there were no apparent changes in domain architecture, gene organisation or gene size. further comparisons between more distant genomes were also uninformative due to the small number of orthologous protein clusters identified by cd-hit or by bidirectional best-scoring blast hits. thus, in summary, comparisons between same-family members show similar trends to those from same-genus comparisons albeit at higher frequencies. complexity another notable feature of myxobacterial tcs proteins is their complexity – both in terms of their domain architectures, and their gene organisation. we define here the most complex systems as ‘intricate’, where at least two receiver and at least two transmitter domains are found in the same protein and/or genetic focus. trtr proteins are encoded as singletons in four genomes, while sa also encodes a trttrt protein. according to our analysis, the trtr proteins are not orthologous to each other or other myxobacterial tcs proteins, while the trttrt protein appears to be a tandem duplication and fusion of two adjacent trt proteins. in addition, myxobacterial genomes possessed between two and nine intricate multi-gene tcs foci, considerably more than the bacterial mean. the intricate foci of ho, sc and af were largely unique to those genomes, otherwise orthologous relationships between the intricate foci were relatively clear at the genus and family level. the foci highlighted exhibit several interesting properties. there seems to be a great deal of plasticity in these foci with many exhibiting orthology with non-intricate tcs foci as a consequence of gene gain/loss and/or recombination. in addition, a large proportion of the tcs domains encoded at intricate foci are found in hybrid kinases. the proportion of foci containing large numbers of tcs domains are similar between bacteria and myxobacteria. nevertheless, when considering all tcs proteins, those of myxobacteria are generally more complex than other bacteria. table  shows the relative frequencies of tcs proteins with different domain compositions for myxobacteria and bacteria. myxobacteria possess far greater proportions of tcs proteins with two or more transmitter domains, and hybrid kinases with three or more receiver domains than other bacteria.table transmitter and receiver domain composition of tcs proteins values are the % of myxobacterial tcs proteins exhibiting the number of receiver and transmitter domains in question divided by the % of all such bacterial tcs proteins athere are no such proteins in bacteria bno such proteins in myxobacteria. numbers of such myxobacterial proteins are shown in parentheses tcs gene expansion nearly two thirds of all mx/mf singletons were found to have orthologues in cc and/or sa, implying that the formation of the majority of singletons was through progressive gene loss. however, genes can also be gained by duplication and divergence. yet duplication is responsible for a relatively small amount of differences between tcs gene sets. for instance, clustering reveals only lineage-specific duplications in sa, four cases in cc and two in mf. eight duplications were apparent within the anaeromyxobacter strains. nevertheless, it appears that certain tcs families have been expanded hugely in the myxobacteria and this occurred relatively recently. for instance, ntrc family members are unusually numerous. indeed, eight of the mx ntrc family members exhibited more similarity to other mx ntrc paralogues than to homologues in other organisms, suggesting origin by duplication after the divergence of contemporary genera of myxobacteria. another mechanism that can increase gene number is horizontal gene transfer. to test the frequency of horizontal acquisition, the mx/mf singletons lacking homologues in cc/sa/mx/mf were queried against the nr protein database using blast. eight of the ten mf singletons and six of the mx singletons gave blast hits across the myxococcales, implying linear descent. however nine of the mx/mf genes gave best non-self blast hits to non-myxococcales genomes, implying horizontal acquisition. such evidence of hgt is also abundant in the singletons of other myxobacteria ). methods phylogenetic analysis 16s rrna gene sequences were aligned using clustalw, using default parameter settings, and neighbour-joining trees were created using phylip with bootstraps, excluding columns containing gaps. inference of tcs proteins sets eight myxobacterial genomes and that of bb are included within the latest version of the p2cs database, available at www.p2cs.org, which provides an automated and consistent pipeline for the identification and annotation of tcs genes. the same categorisation pipeline is available for unpublished sequences at www.p2rp.org and this service was used to characterise the tcs gene sets for sa, cc, pp and ca nucleotide sequence files from www.ncbi.nlm.nih.gov. tcs gene sets were compared with those previously published for mx, sa, adc and sc, and discrepancies manually curated. for the tcs gene sets available through p2cs, potentially ‘mispredicted’ tcss were identified, and these were validated manually and added to the tcs datasets. these two processes resulted in the addition of five, six and tcs proteins to the gene sets for mx, sc and adc respectively, as originally described by whitworth and cock. three, one, one and three ‘mispredicted’ tcs genes were also validated for mf, ho, ad and af respectively. tcs datasets for each genome are available as additional data. characterisation of tcs proteins domain architectures and the presence of transmembrane helices within tcs proteins were determined as described by barakat et al.. tcs gene organisation was defined using a classification scheme based on tcs domains encoded unidirectionally at a ‘focus’, as proposed by williams and whitworth. if a transmitter and a receiver domain were found adjacent in the genome, and they had at least two non-tcs genes on each side, then the relevant genes were described as ‘paired’. if a gene contained a single transmitter, receiver or phosphotransfer domain and was found more than two genes from another tcs gene, it was described as ‘orphan’. tcs gene foci with any other domain/gene organisation were described as ‘complex’. the genomes of pp and ca are incomplete, so gene organisations were not defined for those organisms. focus tcs domain organisation is described in the text using r, t and h to represent receiver, transmitter and hiska/hpt domains respectively, with commas separating genes. thus a ‘rt,r’ organisation represents a hybrid kinase gene encoding an n-terminal receiver domain and a c-terminal transmitter domain, adjacent to a response regulator gene encoding a single receiver domain. sequence-based clustering to define orthologues and paralogues, sets of tcs genes were clustered using cd-hit. a high similarity cut-off was used to give high confidence clusters, and then a lower similarity cut-off was used to identify more divergent clusters. clusters obtained were then subjected to curation, whereby any large clusters containing multiple paralogues obtained using the 30 % cut-off were manually split into the high-confidence clusters obtained with a 90 % cut-off. clusters which included proteins of significantly different length, or of varying domain architecture, family membership and/or gene organisation, or with varying numbers of members per genome were earmarked for further analysis. discussion genetic distance and inferred evolutionary events comparing the tcss from organisms with differing degrees of relatedness allowed inference of the contribution of different types of mutation to the evolution of contemporary sets of tcss. it appears that the relative impact of different mutational events are not affected by the relatedness of the organisms being compared, however, the apparent frequency of mutations is dependent on the relatedness of compared genomes, echoing the findings of he et al. in their study of clostridia. this means that when comparing the gene sets of multiple genomes, it is important to characterize the similarity of those genomes and compensate for differing degree of relatedness. it is even worth including organisms of diverse relatedness, to increase the power of analysis, as suggested by eisen et al.,. in this study we used 16s rrna gene sequence similarity as a straightforward indicator of genomic relatedness/distance and found that the ability to identify tcs orthologues between genomes was dependent on relatedness. at low levels of relatedness, difficulties in the identification of orthologues reduced the power of comparative analysis making it difficult to infer evolutionary events. this difficulty was likely exacerbated by the large numbers of tcs in myxobacterial genomes, and should be less of a problem for the majority of prokaryotes. there have been several studies comparing the tcs gene sets of related organisms, and although not noted, the phenomena we describe can often be seen in the tcs gene sets of diverse taxa. for instance, tcs gene complements can vary by 10–30 % when comparing organisms within the same family/genus, and the frequency of within-species variation is about half that observed when comparing organisms within the same genus. evolution of myxobacterial tcs gene sets our data suggest that most changes in myxobacterial tcs gene sets are due to large indels, some of which are associated with chromosomal recombination. there is evidence of frequent gene loss, but also gene gain via hgt and occasional gene duplications. some of these indels result in changes in gene organization, when the indel only affects part of a tcs locus. small indels are less frequent and often result in domain gain/loss with associated changes in tcs family membership. other studies have reported broadly consistent results. for example, in a study of recombination rates in salmonella, sun et al., found that deletions are significantly more frequent than duplications. both mx and bb exhibit genome-wide evidence of hgt, and we find that hk and rr genes tend to be acquired/lost as pairs, which has been observed in several other taxa, including mutans streptococci and xanthomonas spp.. domain shuffling appears to have happened occasionally during tcs evolution, giving rise to hks with a multiplicity of input domains, and rrs with diverse output domains. our methods of defining orthologues used whole gene sequences and may not have been capable of identifying shuffled domains, nevertheless we did observe several instances of input/output domain gain/loss. with domain shuffling being a rare evolutionary event, it seems probable that apparent domain shuffling is a consequence of domain loss followed by domain gain rather than direct shuffling between proteins. we also did not observe significant relocation of tcs within genomes. the few cases observed were associated with breakpoints of chromosomal recombination. it is possible that those chromosomal relocations occurred due to recombination between homologous tcs genes at those points, however there is too little conservation of synteny in the break point regions to provide any evidence of this happening in myxobacterial genomes. myxobacteria were chosen as the subject of this study due, in part, to the large number of tcs genes they possess. potentially, this could limit extrapolation of our conclusions to other organisms, as increased numbers of tcs genes may have promoted some mutagenic events, and the large numbers of paralogous genes makes robust orthologue identification difficult. in we described the tcs of four sequenced myxobacteria and used domain-based orthology searches to identify likely evolutionary events and peculiarities of myxobacteial tcs gene sets. for this study we added a further eight myxobacterial genomes, dramatically increasing the number of inter-genome comparisons possible, and defined orthologues using whole gene sequence-based methods. the results obtained here agree well with that previous work, and provide a less anecdotal analysis of evolutionary events. benchmarking our genome comparisons by relatedness also allows straightforward comparison with tcs gene sets from other genomes in the future. the tcs genes of myxobacteria seem to be relatively conserved compared to non-tcs genes, and we also found that the mx tcs genes with mf orthologues were enriched for named genes, which tend to have observable phenotypes. this suggests that there is a ‘core’ set of tcs genes in myxobacteria responsible for regulating important functions, with other ‘accessory’ tcs genes being more variable between genomes, and with minor observable phenotypes. the ad1/adc/af/ak tcs gene sets exhibited much less variability than those of mx/mf/cc/sa, however this effect was proportional with the reduced numbers of tcs genes in the anaeromyxobacteria. it is possible that the fruiting myxobacteria need enhanced ‘individuality’ to allow easy discrimination of competing genotypes, but whether this is a selected trait, or purely a consequence of their increased numbers of tcs is uncertain. among strain-specific genes, there is particularly good evidence of acquisition by hgt, which is known to have had a major effect on the shaping of contemporary myxobacterial genomes, alongside lineage-specific expansion. myxobacterial genomics the field of myxobacterial genomics is currently hampered by the large number of genes encoded in each genome, the lack of easily assayable behavioural outputs, and pleiotropy of many genes. this means that phenotypic properties cannot easily be correlated with particular genomic differences. however, further myxobacterial genomes are becoming available at an accelerating rate, and it is possible that in the near future, physiological characterisation of the new genome host organisms will enable correlations with genomic features. to that end, it is important to develop benchmarked catalogs of genomic differences. taking into consideration the relative distance of compared genomes is an important feature of benchmarking, and this study provides a framework to develop such benchmarks using 16s ribosomal rna gene sequence differences as indicators of genetic distance, which can be related to formal taxonomy and other systematic methods. during the course of this study draft genomes became available for strains of cystobacter fuscus, cystobacter violaceus, enhygromyxa salina, hyalangium minutum, sandaracinus amylolyticus, and for two further strains of myxococcus xanthus, including dz and dzf, however these genomes are not complete or scheduled for finishing. thankfully, two complete myxobacterial genomes have also became available, myxococcus stipitatus dsm and sorangium cellulosum so0157-. both of these genomes are now characterized in p2cs and seem typical members of their genera, albeit with slightly more tcs genes, except that the lyttr family of rrs is expanded substantially in ms. as more genomes are completed it will become easier to identify more distant orthologues via ‘gap-bridging’ genomes. this will in turn improve deduction of the directionality of evolutionary changes such as indels, and allow us to create a more finely resolved timeline of myxobacterial tcs gene evolution.  focusing on the tcs genes of myxobacteria, this study investigated the relationship between apparent evolutionary events, and genomic distance or ‘relatedness’. as genomic distance increases it becomes more difficult to identify tcs orthologues, and the relative frequency of evolutionary events increases. however, the relative frequency of evolutionary changes remain unaffected by scale, and it seems that in the myxobacteria, as for other organisms, tcs gene evolution is dominated by large indels, with smaller indels being less frequent. tcs genes are ‘born’ by duplication and acquisition through hgt, but the individuality of tcs gene sets appears to be dominated by lineage-specific gene loss. more generally, our results suggest that researchers should be mindful that the numbers of genomes compared can have a substantial effect on the apparent frequency of observed events. this is particularly true when comparing organisms that share different levels of ‘genetic relatedness’, which should not be assumed to always be consistent with the taxonomic assignments of compared organisms. additional files additional file 1: figure s genetic distance and neighbour-joining tree inferred from myxobacterial 16s rrna gene sequences. bootstrap support for each clade was > 95 % in all cases except for the clade indicated by a filled circle. the tree is rooted with b. bacteriovorus. the bar represents substitutions per nucleotide position. additional file 2: figure s orthology relationships within-family. venn diagrams showing the distribution of orthologous genes between organisms. some orthologous clusters exhibited multiple members from the same organism, suggestive of lineage-specific gene duplications. for the purpose of these diagrams, duplicated proteins were treated as singletons. for instance, if a cluster contained members from sa, cc, mx and mf, but the sa protein was duplicated, then was entered in the sa/cc/mx/mf cell, and into the sa only cell. for a pair of sa proteins, was entered into the sa cell. a: sa/cc/mx/mf, b: af/ak/ad1/adc. additional file 3: the sequence clusters of orthologues of every myxobacterial tcs protein. additional file 4: table s intricate multi-gene tcs foci. each row shows all myxobacterial orthologues of intricate multi-gene foci, including non-intricate orthologues, arranged by organism. foci in particular genomes are denoted by the number from the locus tag of the first protein within the focus. ‘//’ refers to an unspecified number of non-tcs genes. additional file 5: table s evidence of horizontal gene transfer amongst cystobacterial tcs proteins. singletons from comparisons between ad1/adc/ak/af, and between mx/mf/cc/sa, are presented alongside the number of singletons acquired by hgt. also tabulated are the number of tcs genes found in foci where every tcs gene was acquired by putative hgt and the number of those which constitute an ‘intact’ tcs additional file 6: datasets of the annotated tcs proteins from each genome. abbreviations tcstwo-component system hkhistidine kinase rrresponse regulator tmtransmembrane competing interests the author declares that he has no competing interests.
genome relatedness/evolutionary distance was found to have a large effect on the apparent frequency of evolutionary events affecting tcs genes, but not on the relative dominance of different types of mutations. however, there is also evidence of extensive acquisition of genes by lateral transfer, with gene duplication also creating new tcs genes. comparative genomics can provide important insights into the evolution of these genes, but inferences are dependent on the relatedness of the compared genomes. myxobacteriamyxococcus xanthushistidine kinaseresponse regulatormutationcomparative genomicsissue-copyright-statement© the author 2015 myxobacterial tcs gene sets were compared, orthologues defined, and changes in tcs properties such as gene organisation, domain architecture and size identified. large indels were the most common changes, often giving rise to gene organisation changes. two-component systems are abundant prokaryotic signaling pathways, whose evolution is of particular importance because of their role in bacterial pathogenicity. this study investigated the relationship between evolutionary distance and tcs evolution in myxobacterial genomes, of which there are several sequenced examples, of varying relatedness, and which encode large numbers of tcss. it also provides insights into the relationship between evolutionary distance and the inference of tcs estudies of tcs evolution beyond the myxobacteriavolution, which may be important for studies of tcs evolutithe online version of this articleon beyond the myxobacteria. smaller indels were also common, sometimes changing domain architecture, and/or leading to pseudogene formation. individuality of myxobacterial tcs gene sets seems primarily due to lineage specific gene loss. this study provides catalogues of myxobacterial tcs gene sets and their orthology relationships, benchmarked against genome relatedness.
 two-component systems are abundant prokaryotic signaling pathways, whose evolution is of particular importance because of their role in bacterial pathogenicity. comparative genomics can provide important insights into the evolution of these genes, but inferences are dependent on the relatedness of the compared genomes. this study investigated the relationship between evolutionary distance and tcs evolution in myxobacterial genomes, of which there are several sequenced examples, of varying relatedness, and which encode large numbers of tcss. methods myxobacterial tcs gene sets were compared, orthologues defined, and changes in tcs properties such as gene organisation, domain architecture and size identified.  genome relatedness/evolutionary distance was found to have a large effect on the apparent frequency of evolutionary events affecting tcs genes, but not on the relative dominance of different types of mutations. large indels were the most common changes, often giving rise to gene organisation changes. smaller indels were also common, sometimes changing domain architecture, and/or leading to pseudogene formation. individuality of myxobacterial tcs gene sets seems primarily due to lineage specific gene loss. however, there is also evidence of extensive acquisition of genes by lateral transfer, with gene duplication also creating new tcs genes.  this study provides catalogues of myxobacterial tcs gene sets and their orthology relationships, benchmarked against genome relatedness. it also provides insights into the relationship between evolutionary distance and the inference of tcs estudies of tcs evolution beyond the myxobacteriavolution, which may be important for studies of tcs evolutithe online version of this articleon beyond the myxobacteria. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords myxobacteriamyxococcus xanthushistidine kinaseresponse regulatormutationcomparative genomicsissue-copyright-statement© the author 2015 