BACKGROUND
cyanobacteria are oxygenic phototrophic prokaryotes from which chloroplasts, the light harvesting organelles in plants, evolved.
some are able to convert atmospheric nitrogen into a form usable for plants and animals.
during earth history, cyanobacteria have raised atmospheric oxygen levels starting approximately  <dig>  -  <dig>  billion years ago and provided the basis for the evolution of aerobic respiration  <cit> .
cyanobacteria have also evolved extensive morphological diversity.
various patterns of cell organization exist, ranging from single-celled to differentiated multicellular forms with branching patterns.
species of this phylum occupy various habitats.
they can be found in marine, freshwater or terrestrial environments, ranging from polar to tropical climate zones.
based on their morphology, they have been divided into five sections  <cit>  .
sections i and ii comprise single-celled bacteria, whereas sections iii to v comprise multicellular forms.
the latter sections are distinguished according to their level of organization.
section iii is multicellular and undifferentiated, sections iv and v are multicellular and differentiated.
the latter have the ability to produce heterocysts for nitrogen fixation and akinetes .
in addition, species in section v have the ability to branch in multiple dimensions.
1species used to test substitutional saturation
different interpretations of multicellularity are currently used  <cit> .
for cyanobacteria, characterization of multicellularity has been described in previous studies  <cit> .
cell to cell adhesion, intercellular communication, and for more complex species, terminal cell differentiation seem to be three essential processes that define multicellular, prokaryotic organisms on this planet  <cit> .
some forms of complexity found in several multicellular eukaryotes are not present in prokaryotes, but simple forms of multicellularity can be identified in three sections of the phylum cyanobacteria.
multicellular patterns comprise basic filamentous forms as found for section iii, as well as more complex forms involving terminal differentiation, present in sections iv and v. in eukaryotes, multicellular complexity ranges from what is comparable to cyanobacteria to cases with up to  <dig> cell types as estimated for higher invertebrates such as arthropods or molluscs  <cit> .
considering that cyanobacterial sections iii, iv and v resemble some of the first forms of multicellular filaments on earth, knowing when and how these shapes evolved would further our understanding of complex life forms.
some of the oldest body fossils unambiguously identified as cyanobacteria have been found in the kasegalik and mcleary formations of the belcher subgroup, canada, and are evaluated to be between  <dig>  billion and  <dig>  billion years old  <cit> .
studies from ~  <dig>  billion year old formations  <cit>  contain both unicellular and multicellular morphotypes of cyanobacteria.
cyanobacteria certainly existed as early as  <dig>  billion years ago, if one accepts the assumption that they were responsible for the rapid accumulation of oxygen levels, known as the "great oxygenation event"  <cit> .
multicellular fossils belonging to the cyanobacteria are well known from the late precambrian  <cit>  and possibly already existed  <dig>  billion years ago.
other microbe-like multicellular filaments even older than  <dig>  billion years have been found several times  <cit> .
some of the latter fossils are morphologically similar to species from the cyanobacterial order oscillatoriales  <cit> , but no clear evidence has been adduced yet.
although biogenicity of some of the oldest fossils has been questioned  <cit> , a large variety of bacteria including anoxic phototrophs already existed by the time cyanobacteria evolved oxygenic photosynthesis  <cit> .
though impressive for prokaryotes, the fragmentary fossil record alone is not sufficient to disentangle the origin of cyanobacteria and their morphological phenotypes.
therefore, additional methods such as phylogenetic analysis provide a promising possibility to gather further clues on the evolution of such a complex phylum.
phylogenetic analyses of cyanobacteria have gained in quantity over the past  <dig> years  <cit> .
these studies have shown that morphological characterization does not necessarily reflect true relationships between taxa, and possibly none of the five traditional morphological sections is monophyletic.
similar morphologies must have evolved several times independently, but details on this morphological evolution are scarce.
analyses assessing characteristics of cyanobacterial ancestors  <cit>  provide not only fundamental information on the history of cyanobacteria, but also on the evolution of life forms in the archean eon.
if one studies phylogenetic relationships based on protein coding genes in bacteria, it is possible to encounter the outcome of horizontal gene transfer   <cit> .
this issue is not as problematic for ribosomal dna  <cit> .
nonetheless, the problem could be potentially reduced by analyzing datasets of concatenated conserved genes.
identification of these genes for phylogenetic analyses is not without difficulty, and requires in an ideal case comparison of complete genome data  <cit> .
in cyanobacteria, many phylogenetic studies have concentrated on specific clades or smaller subsets of known species in this diverse phylum  <cit> .
therefore the genomic data presently available are strongly biased towards certain groups.
in particular, genomic studies in cyanobacteria have emphasized marine species from section i. marine microphytoplankton  are a particularly well studied group  <cit> , reflected by  <dig> sequenced genomes out of  <dig> cyanobacterial genomes sequenced to date .
from species belonging to section iii only two genomes  are known.
for sections iv  and v  molecular data are rare or missing.
as genomic data accumulate, promising phylogenomic approaches to cyanobacteria are being established  <cit> .
despite these advances, it is at present difficult to obtain sequences other then 16s rdna to cover a representative sample of species from all five sections.
the aim of this paper is to use molecular phylogenetic methods to address the evolutionary history of cyanobacteria and the evolution of multicellularity.
for this purpose, we established a phylogeny based on 16s rdna sequences belonging to  <dig>  cyanobacterial taxa.
from that phylogeny we sampled  <dig> cyanobacterial taxa that represent all main clades obtained and all five sections described by castenholz et al.
<cit> , and feature a 1: <dig> ratio of unicellular to multicellular species.
we used several methods to reconstruct the morphological evolution of ancestral lineages, and compared our results to known fossil data.
since the fossil record is inconclusive on the timing and taxonomic position of multicellular cyanobacteria, our study provides independent evidence on the first appearance and evolution of multicellularity among the ancestors of living cyanobacteria.
RESULTS
phylogenetic analysis
phylogenetic analyses of all identified cyanobacteria
to infer the evolution of multicellularity in cyanobacteria we carried out several phylogenetic analyses.
to ensure a correct taxon-sampling, a phylogeny containing  <dig>  16s rdna sequences of cyanobacteria obtained from genbank was reconstructed .
cyanobacterial morphotypes were assigned to four groups  which correlate to the five sections described by castenholz et al.
<cit> .
using this nomenclature, sub-groups in the phylogeny were assigned to one of the four different morphological groups  according to their dominant shape.
in total  <dig> sub-groups were identified for the phylogenetic tree.
five sub-groups consist of unicellular species from section i , two sub-groups are composed of single celled section ii bacteria , four sub-groups are made up of multicellular species belonging to section iii  and two sub-groups cover differentiated species from section iv and v .
one sub-group contains both species from section i and iii and is therefore designated as ac <dig>  the phylogeny further contains six chloroplast genomes from the eukaryotic phyla glaucophyta, rhodophyta and the division chlorophyta.
chloroplast sequences branch close to the bases and form a sister group to the cyanobacterial sub-groups mentioned.
furthermore six different eubacteria were included in the phylogeny.
they appear to form a distinct outgroup to the cyanobacteria and chloroplasts.
phylogenetic analyses to identify an outgroup
rooted and unrooted phylogenetic analyses reconstructed with maximum likelihood and bayesian inference and based on 16s rrna gene sequences of  <dig> eubacterial species, including  <dig> cyanobacteria revealed congruent results.
cyanobacteria form a monophyletic group.
figure  <dig> shows the unrooted bayesian consensus tree which supports cyanobacterial monophyly with posterior probabilities /bootstrap values  of  <dig> /100%.
phylogenetic trees constructed with an archaean outgroup support cyanobacterial monophyly with pp/bv of  <dig> /98% .
in both cases, plantomyces brasilienses and chlamydia trachomatis, both gram negative bacteria, form a sister group to the cyanobacteria.
this does not agree with other studies  <cit> , where deinococcus-thermus was suggested to be the closest eubacterial relative to cyanobacteria.
these discrepancies may be due to a lack of information when solely using 16s rrna gene sequences for such distant relations.
furthermore, our results confirm the basal position of gloeobacter violaceus, closest to the rest of the eubacteria, as found elsewhere  <cit> .
this supports previous findings which state that gloeobacter violaceus diverged very early from cyanobacteria living today  <cit> .
gloeobacter shows differences in cell structure and metabolism that clearly distinguish it from the rest of extant cyanobacteria  <cit> .
it lacks thylacoid membranes and many genes from photosystems i and ii.
phylogenetic relations of the other eubacterial species show only weak support and are therefore not discussed further.
we separately tested each of the  <dig> eubacterial species originating from a diverse set of non-cyanobacterial phyla, with a subset of the cyanobacteria .
the latter were chosen from the large dataset containing  <dig>  taxa, and cover all sub-groups of the tree .
this subset was used for all subsequent phylogenetic analyses.
though multicellular species seem to dominate the known cyanobacteria, we chose to sample a taxa set containing unicellular and multicellular morphotypes in a 1: <dig> ratio.
that way biases towards certain character states would be excluded.
furthermore, taxa used in the analyses should represent species from all five sections described by castenholz et al.
<cit> .
given our interest in the base of the phylogeny, a greater number of taxa were sampled from basal sub-groups.
due to a lack of data available on genbank at the present state of research, efforts to build a phylogenetic reconstruction of this size  using additional ribosomal protein sequences failed.
but genomic data are accumulating  and will soon offer possibilities for further extensive analyses.
results of six phylogenetic trees are displayed in figure  <dig> .
the majority of the trees exhibit a topology that agrees with figure  <dig>  with the position of gloeobacter violaceus close to the outgroup.
strong differences are found in group support within the trees.
in  <dig> of the  <dig> trees, three nodes could be identified which lead to three clades, named here e ), ac and c .
gloebacter violaceus and synechococcus p <dig> are found at the base of the cyanobacterial phylogeny in  <dig> trees, from which  <dig> trees exhibit gloeobacter violaceus closest to the eubacterial outgroup.
in total  <dig> trees showed congruent topologies.
from the  <dig> eubacteria which have been used as an outgroup in these trees, we chose beggiatoa sp.
as an outgroup for further analyses because its 16s rrna gene sequence exhibits the shortest distance to the cyanobacteria.
phylogenetic analyses of a cyanobacterial subset
phylogenetic analyses of 16s rrna gene sequences from a subset of  <dig> cyanobacterial taxa were conducted using maximum likelihood  and bayesian inference .
for taxa that diverged a long time ago, there is a possibility of sequence saturation, in which case further mutations would have no effect on the distance between sequences any more.
we could significantly reject the possibility of sequence saturation for our alignment .
a general substitution model  was applied for both analyses.
results of the maximum likelihood and bayesian methods are highly congruent.
result of the bayesian analysis with posterior probabilities  and bootstrap values  displayed at the nodes is pictured in figure  <dig>  posterior probabilities above  <dig>  and bootstrap values over 70% are considered to represent a high phylogenetic support.
bootstrap values between 50% and 70% are considered weak support.
posterior probabilities below  <dig>  and bootstrap values below 50% are not displayed.
at deep nodes, the tree topology is fully resolved with high posterior probabilities.
apart from section v, none of the morphological sections described by castenholz et al.
<cit>  is monophyletic.
compared to the outgroup beggiatoa sp., branch lengths are relatively short, which seems surprising given the old age of the phylum.
rates of evolution in cyanobacteria are extremely slow.
this so called "hypobradytelic" tempo would explain their short evolutionary distances  <cit> .
cyanobacteria form three distinct clades mentioned earlier .
clades e, ac and c exhibit posterior probabilities /bootstrap values  of  <dig> /51%,  <dig> /-, and  <dig> /97% respectively .
clade e comprises all taxa analyzed from section ii, some from section i , some from section iii  and all from sections iv and v. within clade e two subclades, e <dig>  and b , are found.
clade ac contains species from section i and iii .
clade c consists of pseudanabaena species, arthronema gygaxiana and phormidium mucicola belonging to section iii.
gloeobacter violaceus is placed closest to the outgroup.
several phylogenetic studies were conducted showing approximate agreement with the tree topology generated here  <cit> .
to check the consistency of results from the maximum likelihood and bayesian analysis to previous studies, we compare our results to the trees produced by honda et al.
<cit> , turner et al.
<cit>  who used 16s rdna sequences, and swingley et al.
<cit>  who used a genomic approach.
the tree from figure  <dig> in honda et al.
<cit>  shows overall strong congruences with our tree.
the only exception is that in honda et al.
<cit>  "synechococcus elongatus toray" is placed separately between gloeobacter and the rest of the cyanobacteria.
we found that "synechococcus elongatus toray"  is located within clade ac in our study and not next to gloeobacter violaceus.
in turner et al.
<cit> , the major clades are congruent with those inferred in our study, but there are a few differences in the relationships among these clades.
in that study, the analog of clade e <dig> is sister to clade ac, which is not the case in our consensus tree.
furthermore, synechococcus c <dig> is grouped with synechococcus p <dig>  which might be due to long branch attraction.
in our phylogenetic tree, synechococcus c <dig> is grouped within clade ac, a relationship supported by high posterior probabilities and bootstrap values .
clade c in our study is placed in the same position as in the tree from turner et al.
<cit> .
swingley et al.
<cit> , used a phylogenomic approach to investigate cyanobacterial relationships.
due to limited, biased genome data available at present, some clades present in our tree are missing in that study.
even so, the main clades retrieved in that study are mostly congruent with clades in our tree.
monophyly of section v  shown in our tree agrees with turner et al.
<cit>  and other studies  <cit> .
nonetheless it is possible that the monophyly of section v bacteria is due to limited taxon sampling, since polyphyly has been detected for section v in another study  <cit> .
gloeobacter violaceus is placed as the first diverging lineage in the phylogeny after the outgroup, as suggested by previous studies  <cit> .
our phylogenetic reconstruction also confirms the placement of taxa belonging to section i and iii throughout the tree  <cit> .
the finding that possibly none of the traditional morphological sections are monophyletic, clearly indicates that similar morphologies have been gained and lost several times during the evolutionary history of living cyanobacteria.
overall, the strong phylogenetic agreement between this and earlier studies confirms the suitability of the tree presented here for further analyses of morphological evolution.
ancestral character state reconstruction
our analysis indicates that multicellularity is a phylogenetically conservative character .
if the terminal taxa of the bayesian consensus tree are randomly re-shuffled, a count through  <dig>  re-shuffled trees gives an average of  <dig> transition steps.
however an average of only nine parsimonious transitions was observed in a count through  <dig>  randomly sampled trees of our ancestral character state reconstruction.
results of the character state reconstruction using the asymmmk model with transition rates estimated by mesquite  <dig>   <cit>  are displayed in figure  <dig>  using maximum likelihood analysis, average frequencies of the characters were counted across  <dig>  trees randomly sampled from the two metropolis-coupled markov chain monte carlo  searches of the bayesian tree reconstruction.
cyanobacteria share a unicellular ancestor, but multicellularity evolved early in the cyanobacterial lineage.
we identified multicellular character states for three basic ancestors leading to clades e, ac and c in our tree.
together, these clades encompass the entirety of the morphological sections ii, iii, iv and v. additionally character states were reconstructed using maximum likelihood analysis and fixed transition rates to analyze properties of the data set.
transition rates are presented in table  <dig>  probabilities for character states at nodes  <dig>   <dig> and  <dig> were examined in detail .
a multicellular ancestry is very likely for these three nodes.
for node  <dig> the relative probabilities of a multicellular ancestor range from  <dig>  to  <dig> , depending on the probability of the transition rates.
for node  <dig> with varying transition rates, the relative probabilities of a multicellular ancestor range from  <dig>  to  <dig> .
for node  <dig> the probabilities for multicellularity range from  <dig>  to  <dig> .
1asymmetrical markov k-state  <dig> parameter model; rates estimated from the consensus tree
2markov k-state  <dig> parameter model; rates estimated from the consensus tree
3f1-f6: models using different fixed transition rates
4reversible jump for model selection, using a hyper prior
5forward rate describing changes to multicellularity
6backward rate describing changes back to a unicellular state
1maximum likelihood: average frequencies across trees were calculated
2maximum parsimony: uniquely best states across trees were counted
3bayesian analysis: model parameters estimated based on the data
the maximum likelihood analysis is not contradicted by a maximum parsimony optimization .
applying maximum parsimony as a reconstruction method, the uniquely best states were counted across  <dig>  trees randomly sampled from the two  runs of the bayesian tree reconstruction.
the relative probabilities for a multicellular ancestor at nodes  <dig>   <dig> and  <dig> are  <dig> ,  <dig>  and  <dig> , respectively.
in contrast, the relative probabilities for a unicellular ancestor at nodes  <dig>   <dig> and  <dig> under parsimony reconstruction are  <dig> ,  <dig>  and  <dig> , respectively.
using bayesian methods, a similar pattern is observed for these nodes.
as an evolutionary model, bayesfactors revealed that a "hyperprior" approach with exponential prior distributions, whose means were sampled from a uniform distribution between  <dig> and  <dig> gave the best fit.
transition rates were estimated to be almost equal.
figure  <dig> displays the posterior probability distributions of character states at these three nodes as they were estimated over  <dig>  randomly sampled trees.
at nodes  <dig> and  <dig> posterior probabilities of a multicellular character state display values above  <dig>  for most of the trees.
at node  <dig> a multicellular state is more likely as well.
posterior probabilities at node  <dig> are above  <dig>  for most of the trees.
at least five reversals to unicellularity occurred in the tree, three of them within clade ac.
the first transition occurred on a branch which led to a group of thermophilic cyanobacteria: acharyochloris sp., synechococcus lividus c <dig> and thermosynechococcus elongatus.
posterior probabilities  and bootstrap values  for this group are  <dig> /73%, whereas the sister group within ac is supported by  <dig> /66% .
the second transition within clade ac led also to a thermophilic cyanobacterium synechococcus c <dig>  sister relation of this species to a filamentous thermophilic cyanobacterium is supported by  <dig> /99% .
the last transition in clade ac occurred within the group including the marine pico-phytoplankton genera synechococcus and prochlorococcus.
the filamentous prochlorothrix hollandica is supposed to be the closest relative to the group that includes marine pico-phytoplankton, supported by  <dig> /61% .
clade ac has a pp of  <dig> , while its bv is below 50%.
although bootstrap support is below 70% for clade ac and some groups within it, posterior probabilities show a very high support .
simulation studies have shown that posterior probabilities approach the actual probability of a clade  <cit> .
bootstrapping tends to underestimate the actual probability of a true clade.
although, posterior probabilities tend to be erroneous if the model of evolution is underparameterized, overparameterization has only a minor effect on the posterior probabilities.
therefore, using a complex model of evolution, such as the "general time reversible with gamma distributed rate variation", is recommended  <cit> .
we used the gtr+g+i model for our analysis, and assume that nodes with a pp higher than  <dig>  are reliable.
it is very likely that at least one additional reversal to unicellularity occurred in clade e <dig>  but phylogenetic support is not high enough to locate the exact position of this transition.
similarly, support for the nodes where the other transition to multicellularity within clade e occurred is missing.
the exact locations of reversals within clade e therefore are not certain and a scenario where multiple reversals occurred cannot be excluded.
in clade e, there is also a reversal to multicellularity observed in spirulina sp.
pcc  <dig>  the location of this transition is supported by posterior probabilities of  <dig>  at two ancestral nodes.
stucken et al.
<cit>  compared gene sets of multicellular cyanobacteria and found that at least  <dig> genes are essential for the formation of filaments.
besides genes previously thought to be correlated with heterocyst formation  they found seven genes coding for hypothetical proteins.
the species they compare are all located within clade e in our tree, most of them being differentiated.
unfortunately no genomes from multicellular species in more basal clades are available at present.
but genome projects of phormidium sp.
isc  <dig> and plectonema sp.
isc  <dig> are presently being conducted http://www.ncbi.nlm.nih.gov/genomes/lproks.cgi.
if these species turn out to group with phormidium mucicola iam m- <dig> and plectonema sp.
f <dig> from the basal clades c and ac in our study, this could provide important information on the original metabolic pathways in ancient multicellular cyanobacteria and on possible advantages of multicellularity.
the majority of cyanobacteria living today are described as successful ecological generalists growing under diverse conditions  <cit> .
our analysis indicates that this diverse range of cyanobacterial morphotypes found in various habitats today —whether multicellular or unicellular— has evolved from multicellular ancestors.
gaining and losing multicellularity
in eukaryotes, simple multicellular forms build the foundation for the evolution of complex multicellular organisms.
although complex multicellularity exhibiting more than three cell types is presumably missing in prokaryotes, bacteria invented simple multicellular forms possibly more than  <dig>  billion years earlier than eukaryotes  <cit> .
multicellularity has been described as one of several major transitions that occurred in the history of life.
these transitions between different units of selection  <cit>  resulted in changes in the organizational confines of the individual.
maynard smith and szathmary  <cit>  summarize eight major transitions in the evolution of life after which, "entities that were capable of independent replication before the transition can replicate only as part of a larger whole after it".
these transitions can create new units of selection at a higher level of complexity  <cit> .
origin of chromosomes, origin of the eukaryotic cell, origin of multicellular organisms and the origin of eusocial communities are some major transitions that redefine the degree of individuality  <cit> .
some transitions are thought to be unique, such as the evolution of meiosis or the evolution of the genetic code.
other major transitions occurred several times independently, such as the evolution of eusociality  <cit>  and multicellularity  <cit> .
there is a tendency to assume that these transitions occur in a progression that leads to an increase in complexity.
however, it seems that in cyanobacteria this is not the case.
anatomical complexity has been lost during their evolution several times .
in a similar fashion, a complex character such as eusociality has been lost several times in halictid bees  <cit> .
conversely the phylogeny indicates that multicellularity re-evolved in spirulina.
regaining complex characters has been observed in other studies as well  <cit> .
nonetheless, some studies state that re-evolution of a complex character after a previous loss is not possible  <cit> .
such studies argue that according to 'dollo's law', a loss of complexity is irreversible  <cit> , a statement that is not supported in the cyanobacterial case.
repeated transitions in either direction are possible.
prokaryotic fossil record before the "great oxygenation event": evidence for multicellular cyanobacteria?
various claims for life during the early archean eon, more than  <dig>  billion years ago exist.
most of them from two regions: the berberton greenstone belt, south africa  and the pilbara craton, western australia .
for some of these "fossils" a biological origin is questioned  <cit> , but for others biogenicity is very likely  <cit> .
these candidates for early life have clear age constraints and there is no non-biological explanation for these structures.
the ages and possible metabolic features of seven fossils of proposed biological origin are plotted in figure  <dig>   <cit> .
some of these fossils are assumed to have been photosynthetic and mat builders, characteristics that can be identified in cyanobacteria as well.
one of the oldest fossils recorded,  <dig>  billion year old prokaryotic remains found in the panorama formation, east pilbara craton, western australia exhibit a filamentous morphotype and possibly carried out anoxygenic photosynthesis  <cit> .
some late archean fossils show an oscillatorian or chroococcacean morphotype .
<dig>  and  <dig>  billion year old oscillatorian-like fossils  <cit>  could possibly represent close relatives of cyanobacterial ancestors.
<dig>  billion year old filamentous bacteria  <cit>  could potentially represent one of the first multicellular cyanobacteria detected.
for single celled forms,  <dig>  billion year old unicellular fossils  <cit>  could likely represent chroococcacean fossils, relatives of ancestral gloeobacter violaceus or synechococcus sp.
p <dig> .
the first conclusive cyanobacterial fossils from all five sections have been reported from around  <dig>  billion year old rocks.
in  <dig>  hofmann described microfossils from stromatolitic dolomite stones in the kasegalik and mcleary formations of the belcher supergroup in hudson bay, northern canada.
among these fossils are halythrix which seems to belong to the order oscillatoriales , eosynechococcus and entophysalis both presumably order chroococcales  and myxococcoides fossils .
in  <dig> similar fossils were described by amard and bertrand-safarti in paleoproterozoic cherty stromatolites from the "formation c " of the franceville group in gabon, dating back  <dig>  billion years.
they also characterized chroococcalean fossils, particularly eosynechococcus and tetraphycus, and filamentous bacteria  which could likely resemble cyanobacteria and myxococcoides fossils.
furthermore, large microfossils , with akinetes similar to the ones from anabaena-like species were found  <cit> .
akinetes are resting cells which are only present in differentiated cyanobacteria from sections iv and v. as it has been confirmed in several studies, sections iv and v share a most recent common ancestor  <cit> .
therefore these fossil akinetes document the existence of differentiated cyanobacteria  <dig>  billion years ago.
given that differentiation in cyanobacteria is evolutionary stable only in a multicellular setting  <cit> , this again supports the notion that multicellular species belonging to the cyanobacteria must have existed earlier than  <dig>  billion years ago.
several studies have assessed prokaryotic history using phylogenetic dating methods  <cit> .
in these studies the origin of cyanobacteria has been estimated around the time of the "great oxygenation event" of  <dig> - <dig>  billion years ago  <cit> .
other studies have reported elevations of oxygen levels before the great rise of atmospheric oxygen  <cit> .
using small and large ribosomal subunit sequences, blank and sanchez-baracaldo  <cit>  estimated the origin of cyanobacteria between  <dig>  and  <dig>  billion years ago.
they also try to address the evolution of cyanobacterial traits and assess that multicellular cyanobacteria did not originate before  <dig> - <dig>  billion years ago.
in the study of blank and sanchez-baracaldo  <cit> , a smaller set of cyanobacterial taxa was used, with some basal multicellular species that are present in clade c of our analysis missing.
these taxa could have an essential effect on the timing of the first multicellular cyanobacteria.
to resolve this issue further dating analyses would be needed.
clearly, as blank and sanchez-baracaldo point out, for such analyses to ultimately resolve the cyanobacterial history, a larger number of cyanobacterial genome data would be needed to represent all the morphological and genetic diversity within this phylum.
CONCLUSIONS
cyanobacteria, photosynthetic prokaryotes, are one of the oldest phyla still alive on this planet.
approximately  <dig> - <dig>  billion years ago cyanobacteria raised the atmospheric oxygen level and established the basis for the evolution of aerobic respiration  <cit> .
they introduced a dramatic change in the earth's atmosphere, which might have created possibilities for more complex lifeforms to evolve.
considering the importance of cyanobacteria for the evolution of life, it seems unfortunate that data sets for a representative phylogenomic analysis are not yet available.
a coordinated perspective between research groups and a diversified taxon sampling strategy for genome projects would offer the possibility for more comprehensive studies on cyanobacterial evolution.
by presenting results obtained from 16s rdna data analysis here, we hope to boost interest for more extensive genomic studies in this phylum.
phylogenomic approaches would help to further investigate some of the results in the present work.
multicellular prokaryotic fossils from the archean eon are documented  <cit> , and fossil data can support the possibility of multicellular cyanobacteria in the archean eon  <cit> .
furthermore, studies describe smaller accumulations of oxygen levels around  <dig>  to  <dig>  billion years ago  <cit>  and around  <dig>  billion years ago  <cit> .
therefore multicellular cyanobacteria could have evolved before the rise of oxygen in the atmosphere.
the "great oxygenation event", also referred to as "oxygen crisis", could presumably have marked one of the first mass extinction events during earth's history.
new habitats developing around  <dig>  billion years ago, due to a dramatic change of earth's atmosphere could have triggered cyanobacteria to evolve the variety of morphotypes preserved until today.
in terms of cell types, cyanobacteria reached their maximum morphological complexity around  <dig>  billion years ago  <cit> .
by the time eukaryotes evolved, cyanobacteria already exhibited the full range of their morphological diversity.
due to slow evolutionary rates in cyanobacteria, which have been described as "hypobradytelic"  <cit> , extant cyanobacteria that appear to exhibit the same morphotype as in the precambrian eon  <cit>  are reminiscent of the idea of "living fossils".
however, one should consider the possibility that what may appear as morphological stasis may be due to developmental constraints at the phylum level.
cyanobacteria apparently reached their maximum complexity early in earth history, but instead of morphological stasis at the species level, our results suggest that they subsequently changed morphotypes several times during their evolution.
this allowed for the exploration of diverse morphotypes within their developmental constraints, including the loss and regaining of multicellular growth forms.

 cyanobacteria are oxygenic phototrophic prokaryotes from which chloroplasts, the light harvesting organelles in plants, evolved. some are able to convert atmospheric nitrogen into a form usable for plants and animals. during earth history, cyanobacteria have raised atmospheric oxygen levels starting approximately - billion years ago and provided the basis for the evolution of aerobic respiration. cyanobacteria have also evolved extensive morphological diversity. various patterns of cell organization exist, ranging from single-celled to differentiated multicellular forms with branching patterns. species of this phylum occupy various habitats. they can be found in marine, freshwater or terrestrial environments, ranging from polar to tropical climate zones. based on their morphology, they have been divided into five sections. sections i and ii comprise single-celled bacteria, whereas sections iii to v comprise multicellular forms. the latter sections are distinguished according to their level of organization. section iii is multicellular and undifferentiated, sections iv and v are multicellular and differentiated. the latter have the ability to produce heterocysts for nitrogen fixation and akinetes. in addition, species in section v have the ability to branch in multiple dimensions. 1species used to test substitutional saturation different interpretations of multicellularity are currently used. for cyanobacteria, characterization of multicellularity has been described in previous studies. cell to cell adhesion, intercellular communication, and for more complex species, terminal cell differentiation seem to be three essential processes that define multicellular, prokaryotic organisms on this planet. some forms of complexity found in several multicellular eukaryotes are not present in prokaryotes, but simple forms of multicellularity can be identified in three sections of the phylum cyanobacteria. multicellular patterns comprise basic filamentous forms as found for section iii, as well as more complex forms involving terminal differentiation, present in sections iv and v. in eukaryotes, multicellular complexity ranges from what is comparable to cyanobacteria to cases with up to cell types as estimated for higher invertebrates such as arthropods or molluscs. considering that cyanobacterial sections iii, iv and v resemble some of the first forms of multicellular filaments on earth, knowing when and how these shapes evolved would further our understanding of complex life forms. some of the oldest body fossils unambiguously identified as cyanobacteria have been found in the kasegalik and mcleary formations of the belcher subgroup, canada, and are evaluated to be between billion and billion years old. studies from ~ billion year old formations contain both unicellular and multicellular morphotypes of cyanobacteria. cyanobacteria certainly existed as early as billion years ago, if one accepts the assumption that they were responsible for the rapid accumulation of oxygen levels, known as the "great oxygenation event". multicellular fossils belonging to the cyanobacteria are well known from the late precambrian and possibly already existed billion years ago. other microbe-like multicellular filaments even older than billion years have been found several times. some of the latter fossils are morphologically similar to species from the cyanobacterial order oscillatoriales, but no clear evidence has been adduced yet. although biogenicity of some of the oldest fossils has been questioned, a large variety of bacteria including anoxic phototrophs already existed by the time cyanobacteria evolved oxygenic photosynthesis. though impressive for prokaryotes, the fragmentary fossil record alone is not sufficient to disentangle the origin of cyanobacteria and their morphological phenotypes. therefore, additional methods such as phylogenetic analysis provide a promising possibility to gather further clues on the evolution of such a complex phylum. phylogenetic analyses of cyanobacteria have gained in quantity over the past years. these studies have shown that morphological characterization does not necessarily reflect true relationships between taxa, and possibly none of the five traditional morphological sections is monophyletic. similar morphologies must have evolved several times independently, but details on this morphological evolution are scarce. analyses assessing characteristics of cyanobacterial ancestors provide not only fundamental information on the history of cyanobacteria, but also on the evolution of life forms in the archean eon. if one studies phylogenetic relationships based on protein coding genes in bacteria, it is possible to encounter the outcome of horizontal gene transfer. this issue is not as problematic for ribosomal dna. nonetheless, the problem could be potentially reduced by analyzing datasets of concatenated conserved genes. identification of these genes for phylogenetic analyses is not without difficulty, and requires in an ideal case comparison of complete genome data. in cyanobacteria, many phylogenetic studies have concentrated on specific clades or smaller subsets of known species in this diverse phylum. therefore the genomic data presently available are strongly biased towards certain groups. in particular, genomic studies in cyanobacteria have emphasized marine species from section i. marine microphytoplankton are a particularly well studied group, reflected by sequenced genomes out of cyanobacterial genomes sequenced to date. from species belonging to section iii only two genomes are known. for sections iv and v molecular data are rare or missing. as genomic data accumulate, promising phylogenomic approaches to cyanobacteria are being established. despite these advances, it is at present difficult to obtain sequences other then 16s rdna to cover a representative sample of species from all five sections. the aim of this paper is to use molecular phylogenetic methods to address the evolutionary history of cyanobacteria and the evolution of multicellularity. for this purpose, we established a phylogeny based on 16s rdna sequences belonging to cyanobacterial taxa. from that phylogeny we sampled cyanobacterial taxa that represent all main clades obtained and all five sections described by castenholz et al., and feature a 1: ratio of unicellular to multicellular species. we used several methods to reconstruct the morphological evolution of ancestral lineages, and compared our results to known fossil data. since the fossil record is inconclusive on the timing and taxonomic position of multicellular cyanobacteria, our study provides independent evidence on the first appearance and evolution of multicellularity among the ancestors of living cyanobacteria.  phylogenetic analysis phylogenetic analyses of all identified cyanobacteria to infer the evolution of multicellularity in cyanobacteria we carried out several phylogenetic analyses. to ensure a correct taxon-sampling, a phylogeny containing 16s rdna sequences of cyanobacteria obtained from genbank was reconstructed. cyanobacterial morphotypes were assigned to four groups which correlate to the five sections described by castenholz et al.. using this nomenclature, sub-groups in the phylogeny were assigned to one of the four different morphological groups according to their dominant shape. in total sub-groups were identified for the phylogenetic tree. five sub-groups consist of unicellular species from section i, two sub-groups are composed of single celled section ii bacteria, four sub-groups are made up of multicellular species belonging to section iii and two sub-groups cover differentiated species from section iv and v. one sub-group contains both species from section i and iii and is therefore designated as ac the phylogeny further contains six chloroplast genomes from the eukaryotic phyla glaucophyta, rhodophyta and the division chlorophyta. chloroplast sequences branch close to the bases and form a sister group to the cyanobacterial sub-groups mentioned. furthermore six different eubacteria were included in the phylogeny. they appear to form a distinct outgroup to the cyanobacteria and chloroplasts. phylogenetic analyses to identify an outgroup rooted and unrooted phylogenetic analyses reconstructed with maximum likelihood and bayesian inference and based on 16s rrna gene sequences of eubacterial species, including cyanobacteria revealed congruent results. cyanobacteria form a monophyletic group. figure shows the unrooted bayesian consensus tree which supports cyanobacterial monophyly with posterior probabilities /bootstrap values of /100%. phylogenetic trees constructed with an archaean outgroup support cyanobacterial monophyly with pp/bv of /98%. in both cases, plantomyces brasilienses and chlamydia trachomatis, both gram negative bacteria, form a sister group to the cyanobacteria. this does not agree with other studies, where deinococcus-thermus was suggested to be the closest eubacterial relative to cyanobacteria. these discrepancies may be due to a lack of information when solely using 16s rrna gene sequences for such distant relations. furthermore, our results confirm the basal position of gloeobacter violaceus, closest to the rest of the eubacteria, as found elsewhere. this supports previous findings which state that gloeobacter violaceus diverged very early from cyanobacteria living today. gloeobacter shows differences in cell structure and metabolism that clearly distinguish it from the rest of extant cyanobacteria. it lacks thylacoid membranes and many genes from photosystems i and ii. phylogenetic relations of the other eubacterial species show only weak support and are therefore not discussed further. we separately tested each of the eubacterial species originating from a diverse set of non-cyanobacterial phyla, with a subset of the cyanobacteria. the latter were chosen from the large dataset containing taxa, and cover all sub-groups of the tree. this subset was used for all subsequent phylogenetic analyses. though multicellular species seem to dominate the known cyanobacteria, we chose to sample a taxa set containing unicellular and multicellular morphotypes in a 1: ratio. that way biases towards certain character states would be excluded. furthermore, taxa used in the analyses should represent species from all five sections described by castenholz et al.. given our interest in the base of the phylogeny, a greater number of taxa were sampled from basal sub-groups. due to a lack of data available on genbank at the present state of research, efforts to build a phylogenetic reconstruction of this size using additional ribosomal protein sequences failed. but genomic data are accumulating and will soon offer possibilities for further extensive analyses. results of six phylogenetic trees are displayed in figure. the majority of the trees exhibit a topology that agrees with figure with the position of gloeobacter violaceus close to the outgroup. strong differences are found in group support within the trees. in of the trees, three nodes could be identified which lead to three clades, named here e ), ac and c. gloebacter violaceus and synechococcus p are found at the base of the cyanobacterial phylogeny in trees, from which trees exhibit gloeobacter violaceus closest to the eubacterial outgroup. in total trees showed congruent topologies. from the eubacteria which have been used as an outgroup in these trees, we chose beggiatoa sp. as an outgroup for further analyses because its 16s rrna gene sequence exhibits the shortest distance to the cyanobacteria. phylogenetic analyses of a cyanobacterial subset phylogenetic analyses of 16s rrna gene sequences from a subset of cyanobacterial taxa were conducted using maximum likelihood and bayesian inference. for taxa that diverged a long time ago, there is a possibility of sequence saturation, in which case further mutations would have no effect on the distance between sequences any more. we could significantly reject the possibility of sequence saturation for our alignment. a general substitution model was applied for both analyses. results of the maximum likelihood and bayesian methods are highly congruent. result of the bayesian analysis with posterior probabilities and bootstrap values displayed at the nodes is pictured in figure posterior probabilities above and bootstrap values over 70% are considered to represent a high phylogenetic support. bootstrap values between 50% and 70% are considered weak support. posterior probabilities below and bootstrap values below 50% are not displayed. at deep nodes, the tree topology is fully resolved with high posterior probabilities. apart from section v, none of the morphological sections described by castenholz et al. is monophyletic. compared to the outgroup beggiatoa sp., branch lengths are relatively short, which seems surprising given the old age of the phylum. rates of evolution in cyanobacteria are extremely slow. this so called "hypobradytelic" tempo would explain their short evolutionary distances. cyanobacteria form three distinct clades mentioned earlier. clades e, ac and c exhibit posterior probabilities /bootstrap values of /51%, /-, and /97% respectively. clade e comprises all taxa analyzed from section ii, some from section i, some from section iii and all from sections iv and v. within clade e two subclades, e and b, are found. clade ac contains species from section i and iii. clade c consists of pseudanabaena species, arthronema gygaxiana and phormidium mucicola belonging to section iii. gloeobacter violaceus is placed closest to the outgroup. several phylogenetic studies were conducted showing approximate agreement with the tree topology generated here. to check the consistency of results from the maximum likelihood and bayesian analysis to previous studies, we compare our results to the trees produced by honda et al., turner et al. who used 16s rdna sequences, and swingley et al. who used a genomic approach. the tree from figure in honda et al. shows overall strong congruences with our tree. the only exception is that in honda et al. "synechococcus elongatus toray" is placed separately between gloeobacter and the rest of the cyanobacteria. we found that "synechococcus elongatus toray" is located within clade ac in our study and not next to gloeobacter violaceus. in turner et al., the major clades are congruent with those inferred in our study, but there are a few differences in the relationships among these clades. in that study, the analog of clade e is sister to clade ac, which is not the case in our consensus tree. furthermore, synechococcus c is grouped with synechococcus p which might be due to long branch attraction. in our phylogenetic tree, synechococcus c is grouped within clade ac, a relationship supported by high posterior probabilities and bootstrap values. clade c in our study is placed in the same position as in the tree from turner et al.. swingley et al., used a phylogenomic approach to investigate cyanobacterial relationships. due to limited, biased genome data available at present, some clades present in our tree are missing in that study. even so, the main clades retrieved in that study are mostly congruent with clades in our tree. monophyly of section v shown in our tree agrees with turner et al. and other studies. nonetheless it is possible that the monophyly of section v bacteria is due to limited taxon sampling, since polyphyly has been detected for section v in another study. gloeobacter violaceus is placed as the first diverging lineage in the phylogeny after the outgroup, as suggested by previous studies. our phylogenetic reconstruction also confirms the placement of taxa belonging to section i and iii throughout the tree. the finding that possibly none of the traditional morphological sections are monophyletic, clearly indicates that similar morphologies have been gained and lost several times during the evolutionary history of living cyanobacteria. overall, the strong phylogenetic agreement between this and earlier studies confirms the suitability of the tree presented here for further analyses of morphological evolution. ancestral character state reconstruction our analysis indicates that multicellularity is a phylogenetically conservative character. if the terminal taxa of the bayesian consensus tree are randomly re-shuffled, a count through re-shuffled trees gives an average of transition steps. however an average of only nine parsimonious transitions was observed in a count through randomly sampled trees of our ancestral character state reconstruction. results of the character state reconstruction using the asymmmk model with transition rates estimated by mesquite are displayed in figure using maximum likelihood analysis, average frequencies of the characters were counted across trees randomly sampled from the two metropolis-coupled markov chain monte carlo searches of the bayesian tree reconstruction. cyanobacteria share a unicellular ancestor, but multicellularity evolved early in the cyanobacterial lineage. we identified multicellular character states for three basic ancestors leading to clades e, ac and c in our tree. together, these clades encompass the entirety of the morphological sections ii, iii, iv and v. additionally character states were reconstructed using maximum likelihood analysis and fixed transition rates to analyze properties of the data set. transition rates are presented in table probabilities for character states at nodes and were examined in detail. a multicellular ancestry is very likely for these three nodes. for node the relative probabilities of a multicellular ancestor range from to, depending on the probability of the transition rates. for node with varying transition rates, the relative probabilities of a multicellular ancestor range from to. for node the probabilities for multicellularity range from to. 1asymmetrical markov k-state parameter model; rates estimated from the consensus tree 2markov k-state parameter model; rates estimated from the consensus tree 3f1-f6: models using different fixed transition rates 4reversible jump for model selection, using a hyper prior 5forward rate describing changes to multicellularity 6backward rate describing changes back to a unicellular state 1maximum likelihood: average frequencies across trees were calculated 2maximum parsimony: uniquely best states across trees were counted 3bayesian analysis: model parameters estimated based on the data the maximum likelihood analysis is not contradicted by a maximum parsimony optimization. applying maximum parsimony as a reconstruction method, the uniquely best states were counted across trees randomly sampled from the two runs of the bayesian tree reconstruction. the relative probabilities for a multicellular ancestor at nodes and are, and, respectively. in contrast, the relative probabilities for a unicellular ancestor at nodes and under parsimony reconstruction are, and, respectively. using bayesian methods, a similar pattern is observed for these nodes. as an evolutionary model, bayesfactors revealed that a "hyperprior" approach with exponential prior distributions, whose means were sampled from a uniform distribution between and gave the best fit. transition rates were estimated to be almost equal. figure displays the posterior probability distributions of character states at these three nodes as they were estimated over randomly sampled trees. at nodes and posterior probabilities of a multicellular character state display values above for most of the trees. at node a multicellular state is more likely as well. posterior probabilities at node are above for most of the trees. at least five reversals to unicellularity occurred in the tree, three of them within clade ac. the first transition occurred on a branch which led to a group of thermophilic cyanobacteria: acharyochloris sp., synechococcus lividus c and thermosynechococcus elongatus. posterior probabilities and bootstrap values for this group are /73%, whereas the sister group within ac is supported by /66%. the second transition within clade ac led also to a thermophilic cyanobacterium synechococcus c sister relation of this species to a filamentous thermophilic cyanobacterium is supported by /99%. the last transition in clade ac occurred within the group including the marine pico-phytoplankton genera synechococcus and prochlorococcus. the filamentous prochlorothrix hollandica is supposed to be the closest relative to the group that includes marine pico-phytoplankton, supported by /61%. clade ac has a pp of, while its bv is below 50%. although bootstrap support is below 70% for clade ac and some groups within it, posterior probabilities show a very high support. simulation studies have shown that posterior probabilities approach the actual probability of a clade. bootstrapping tends to underestimate the actual probability of a true clade. although, posterior probabilities tend to be erroneous if the model of evolution is underparameterized, overparameterization has only a minor effect on the posterior probabilities. therefore, using a complex model of evolution, such as the "general time reversible with gamma distributed rate variation", is recommended. we used the gtr+g+i model for our analysis, and assume that nodes with a pp higher than are reliable. it is very likely that at least one additional reversal to unicellularity occurred in clade e but phylogenetic support is not high enough to locate the exact position of this transition. similarly, support for the nodes where the other transition to multicellularity within clade e occurred is missing. the exact locations of reversals within clade e therefore are not certain and a scenario where multiple reversals occurred cannot be excluded. in clade e, there is also a reversal to multicellularity observed in spirulina sp. pcc the location of this transition is supported by posterior probabilities of at two ancestral nodes. stucken et al. compared gene sets of multicellular cyanobacteria and found that at least genes are essential for the formation of filaments. besides genes previously thought to be correlated with heterocyst formation they found seven genes coding for hypothetical proteins. the species they compare are all located within clade e in our tree, most of them being differentiated. unfortunately no genomes from multicellular species in more basal clades are available at present. but genome projects of phormidium sp. isc and plectonema sp. isc are presently being conducted http://www.ncbi.nlm.nih.gov/genomes/lproks.cgi. if these species turn out to group with phormidium mucicola iam m- and plectonema sp. f from the basal clades c and ac in our study, this could provide important information on the original metabolic pathways in ancient multicellular cyanobacteria and on possible advantages of multicellularity. the majority of cyanobacteria living today are described as successful ecological generalists growing under diverse conditions. our analysis indicates that this diverse range of cyanobacterial morphotypes found in various habitats today —whether multicellular or unicellular— has evolved from multicellular ancestors. gaining and losing multicellularity in eukaryotes, simple multicellular forms build the foundation for the evolution of complex multicellular organisms. although complex multicellularity exhibiting more than three cell types is presumably missing in prokaryotes, bacteria invented simple multicellular forms possibly more than billion years earlier than eukaryotes. multicellularity has been described as one of several major transitions that occurred in the history of life. these transitions between different units of selection resulted in changes in the organizational confines of the individual. maynard smith and szathmary summarize eight major transitions in the evolution of life after which, "entities that were capable of independent replication before the transition can replicate only as part of a larger whole after it". these transitions can create new units of selection at a higher level of complexity. origin of chromosomes, origin of the eukaryotic cell, origin of multicellular organisms and the origin of eusocial communities are some major transitions that redefine the degree of individuality. some transitions are thought to be unique, such as the evolution of meiosis or the evolution of the genetic code. other major transitions occurred several times independently, such as the evolution of eusociality and multicellularity. there is a tendency to assume that these transitions occur in a progression that leads to an increase in complexity. however, it seems that in cyanobacteria this is not the case. anatomical complexity has been lost during their evolution several times. in a similar fashion, a complex character such as eusociality has been lost several times in halictid bees. conversely the phylogeny indicates that multicellularity re-evolved in spirulina. regaining complex characters has been observed in other studies as well. nonetheless, some studies state that re-evolution of a complex character after a previous loss is not possible. such studies argue that according to'dollo's law', a loss of complexity is irreversible, a statement that is not supported in the cyanobacterial case. repeated transitions in either direction are possible. prokaryotic fossil record before the "great oxygenation event": evidence for multicellular cyanobacteria? various claims for life during the early archean eon, more than billion years ago exist. most of them from two regions: the berberton greenstone belt, south africa and the pilbara craton, western australia. for some of these "fossils" a biological origin is questioned, but for others biogenicity is very likely. these candidates for early life have clear age constraints and there is no non-biological explanation for these structures. the ages and possible metabolic features of seven fossils of proposed biological origin are plotted in figure. some of these fossils are assumed to have been photosynthetic and mat builders, characteristics that can be identified in cyanobacteria as well. one of the oldest fossils recorded, billion year old prokaryotic remains found in the panorama formation, east pilbara craton, western australia exhibit a filamentous morphotype and possibly carried out anoxygenic photosynthesis. some late archean fossils show an oscillatorian or chroococcacean morphotype. and billion year old oscillatorian-like fossils could possibly represent close relatives of cyanobacterial ancestors. billion year old filamentous bacteria could potentially represent one of the first multicellular cyanobacteria detected. for single celled forms, billion year old unicellular fossils could likely represent chroococcacean fossils, relatives of ancestral gloeobacter violaceus or synechococcus sp. p. the first conclusive cyanobacterial fossils from all five sections have been reported from around billion year old rocks. in hofmann described microfossils from stromatolitic dolomite stones in the kasegalik and mcleary formations of the belcher supergroup in hudson bay, northern canada. among these fossils are halythrix which seems to belong to the order oscillatoriales, eosynechococcus and entophysalis both presumably order chroococcales and myxococcoides fossils. in similar fossils were described by amard and bertrand-safarti in paleoproterozoic cherty stromatolites from the "formation c " of the franceville group in gabon, dating back billion years. they also characterized chroococcalean fossils, particularly eosynechococcus and tetraphycus, and filamentous bacteria which could likely resemble cyanobacteria and myxococcoides fossils. furthermore, large microfossils, with akinetes similar to the ones from anabaena-like species were found. akinetes are resting cells which are only present in differentiated cyanobacteria from sections iv and v. as it has been confirmed in several studies, sections iv and v share a most recent common ancestor. therefore these fossil akinetes document the existence of differentiated cyanobacteria billion years ago. given that differentiation in cyanobacteria is evolutionary stable only in a multicellular setting, this again supports the notion that multicellular species belonging to the cyanobacteria must have existed earlier than billion years ago. several studies have assessed prokaryotic history using phylogenetic dating methods. in these studies the origin of cyanobacteria has been estimated around the time of the "great oxygenation event" of - billion years ago. other studies have reported elevations of oxygen levels before the great rise of atmospheric oxygen. using small and large ribosomal subunit sequences, blank and sanchez-baracaldo estimated the origin of cyanobacteria between and billion years ago. they also try to address the evolution of cyanobacterial traits and assess that multicellular cyanobacteria did not originate before - billion years ago. in the study of blank and sanchez-baracaldo, a smaller set of cyanobacterial taxa was used, with some basal multicellular species that are present in clade c of our analysis missing. these taxa could have an essential effect on the timing of the first multicellular cyanobacteria. to resolve this issue further dating analyses would be needed. clearly, as blank and sanchez-baracaldo point out, for such analyses to ultimately resolve the cyanobacterial history, a larger number of cyanobacterial genome data would be needed to represent all the morphological and genetic diversity within this phylum.  cyanobacteria, photosynthetic prokaryotes, are one of the oldest phyla still alive on this planet. approximately - billion years ago cyanobacteria raised the atmospheric oxygen level and established the basis for the evolution of aerobic respiration. they introduced a dramatic change in the earth's atmosphere, which might have created possibilities for more complex lifeforms to evolve. considering the importance of cyanobacteria for the evolution of life, it seems unfortunate that data sets for a representative phylogenomic analysis are not yet available. a coordinated perspective between research groups and a diversified taxon sampling strategy for genome projects would offer the possibility for more comprehensive studies on cyanobacterial evolution. by presenting results obtained from 16s rdna data analysis here, we hope to boost interest for more extensive genomic studies in this phylum. phylogenomic approaches would help to further investigate some of the results in the present work. multicellular prokaryotic fossils from the archean eon are documented, and fossil data can support the possibility of multicellular cyanobacteria in the archean eon. furthermore, studies describe smaller accumulations of oxygen levels around to billion years ago and around billion years ago. therefore multicellular cyanobacteria could have evolved before the rise of oxygen in the atmosphere. the "great oxygenation event", also referred to as "oxygen crisis", could presumably have marked one of the first mass extinction events during earth's history. new habitats developing around billion years ago, due to a dramatic change of earth's atmosphere could have triggered cyanobacteria to evolve the variety of morphotypes preserved until today. in terms of cell types, cyanobacteria reached their maximum morphological complexity around billion years ago. by the time eukaryotes evolved, cyanobacteria already exhibited the full range of their morphological diversity. due to slow evolutionary rates in cyanobacteria, which have been described as "hypobradytelic", extant cyanobacteria that appear to exhibit the same morphotype as in the precambrian eon are reminiscent of the idea of "living fossils". however, one should consider the possibility that what may appear as morphological stasis may be due to developmental constraints at the phylum level. cyanobacteria apparently reached their maximum complexity early in earth history, but instead of morphological stasis at the species level, our results suggest that they subsequently changed morphotypes several times during their evolution. this allowed for the exploration of diverse morphotypes within their developmental constraints, including the loss and regaining of multicellular growth forms. 
the results indicate that a multicellular morphotype evolved early in the cyanobacterial lineage and was regained at least once after a previous loss. most of the morphological diversity exhibited in cyanobacteria today —including the majority of single-celled species— arose from ancient multicellular lineages. we conducted and compared phylogenetic analyses of 16s rdna sequences from a large sample of taxa representing the morphological and genetic diversity of cyanobacteria. comparison to the fossil record supports an early origin of multicellularity, possibly as early as the "great oxygenation event" that occurred - billion years ago. the results suggest that the majority of extant cyanobacteria descend from multicellular ancestors. understanding when and how multicellularity evolved in these ancient organisms would provide fundamental information on the early history of life and further our knowledge of complex life forms.
 cyanobacteria are one of the oldest and morphologically most diverse prokaryotic phyla on our planet. the early development of an oxygen-containing atmosphere approximately - billion years ago is attributed to the photosynthetic activity of cyanobacteria. furthermore, they are one of the few prokaryotic phyla where multicellularity has evolved. understanding when and how multicellularity evolved in these ancient organisms would provide fundamental information on the early history of life and further our knowledge of complex life forms.  we conducted and compared phylogenetic analyses of 16s rdna sequences from a large sample of taxa representing the morphological and genetic diversity of cyanobacteria. we reconstructed ancestral character states on phylogenetic trees. the results suggest that the majority of extant cyanobacteria descend from multicellular ancestors. reversals to unicellularity occurred at least times. multicellularity was established again at least once within a single-celled clade. comparison to the fossil record supports an early origin of multicellularity, possibly as early as the "great oxygenation event" that occurred - billion years ago.  the results indicate that a multicellular morphotype evolved early in the cyanobacterial lineage and was regained at least once after a previous loss. most of the morphological diversity exhibited in cyanobacteria today —including the majority of single-celled species— arose from ancient multicellular lineages. multicellularity could have conferred a considerable advantage for exploring new niches and hence facilitated the diversification of new lineages. 