- Research article
Highly plastic genome of Microcystis aeruginosa PCC 7806, a ubiquitous toxic freshwater cyanobacterium
BMC Genomicsvolume 9, Article number: 274 (2008)
The colonial cyanobacterium Microcystis proliferates in a wide range of freshwater ecosystems and is exposed to changing environmental factors during its life cycle. Microcystis blooms are often toxic, potentially fatal to animals and humans, and may cause environmental problems. There has been little investigation of the genomics of these cyanobacteria.
Deciphering the 5,172,804 bp sequence of Microcystis aeruginosa PCC 7806 has revealed the high plasticity of its genome: 11.7% DNA repeats containing more than 1,000 bases, 6.8% putative transposases and 21 putative restriction enzymes. Compared to the genomes of other cyanobacterial lineages, strain PCC 7806 contains a large number of atypical genes that may have been acquired by lateral transfers. Metabolic pathways, such as fermentation and a methionine salvage pathway, have been identified, as have genes for programmed cell death that may be related to the rapid disappearance of Microcystis blooms in nature. Analysis of the PCC 7806 genome also reveals striking novel biosynthetic features that might help to elucidate the ecological impact of secondary metabolites and lead to the discovery of novel metabolites for new biotechnological applications. M. aeruginosa and other large cyanobacterial genomes exhibit a rapid loss of synteny in contrast to other microbial genomes.
Microcystis aeruginosa PCC 7806 appears to have adopted an evolutionary strategy relying on unusual genome plasticity to adapt to eutrophic freshwater ecosystems, a property shared by another strain of M. aeruginosa (NIES-843). Comparisons of the genomes of PCC 7806 and other cyanobacterial strains indicate that a similar strategy may have also been used by the marine strain Crocosphaera watsonii WH8501 to adapt to other ecological niches, such as oligotrophic open oceans.
Dated approximately 3 billion years old by fossil records, cyanobacteria were the first oxyphototrophic prokaryotes present on Earth . As architects of the Earth's atmosphere they had a major impact on the evolution of aerobic metabolism and the evolution of life . Cyanobacteria still play a fundamental role in the functioning of global ecosystems by significantly contributing to carbon fluxes [3, 4] and by providing nitrogen used for primary production . On the other hand, cyanobacterial blooms may lead to a loss of biodiversity in the phytoplanktonic communities and, by generating very high quantities of organic matter used by anoxygenic bacteria in the bottom layers of water resources, can cause massive death of fish by asphyxia . The financial costs resulting from cyanobacterial proliferations are considerable (e.g. 200 million Australian dollars/year in Australia) .
Freshwater cyanobacteria of the genus Microcystis are distributed worldwide, and are involved in numerous proliferation events in stratified lakes . In their natural environment, Microcystis cells are organized in large colonies of various sizes and shapes, which were used to define various morphospecies. Five of these have recently been reunified as a single species, Microcystis aeruginosa . The determinism of the morphogical variations within this polymorphic cyanobacterial species is currently under debate.
The ecology of M. aeruginosa is characterized by an annual life cycle comprising a spring and summer pelagic phase, and an overwintering benthic phase . During the pelagic phase, M. aeruginosa colonies migrate daily in the water column  and may accumulate to form blooms or scums on the surface of the water. Thus, on a daily basis, as well as during the benthic and pelagic phases, colonies are exposed to changing environmental conditions of light, temperature and oxygen concentrations.
In the last decade, cyanobacterial blooms have been involved in numerous cases of animal  and human  poisonings, mainly due to the ability of Microcystis cells to synthesize toxins, in particular variants of microcystin . Many other oligopeptides, such as cyanopeptolins, aeruginosins, microginins, microviridins and cyclamides may also be produced . Other peptides and congeners doubtless remain to be discovered, as do their respective biosynthesis pathways.
To gain further insight into the ecophysiology of Microcystis aeruginosa, we deciphered the genome sequence of the toxic strain PCC 7806. The results presented here associate descriptive genomics and comparisons with the genomes of other cyanobacteria isolated from freshwater and marine ecosystems to highlight the ecophysiological peculiarities of this strain, and put its particularly high genome plasticity into a cyanobacterial context.
Results and discussion
General features of the M. aeruginosa PCC 7806 genome
The 12× shotgun sequencing project produced 90,000 sequence reads, and their assembly resulted in more than 500 contigs. After the first steps of a long finishing process performed using CAAT-Box  and Consed  software, the number of contigs was reduced to 328 (N50 = 100kb), 116 of which were more than 3,000 bases in length (up to 533,374 bases). The genome contains an unusually high number of long DNA repeats. Most of the extremities of these contigs consist of DNA repeated sequences including gene coding for transposases (see below). The 116 contigs were deposited in the EMBL database (AM778843–AM778958). The genome sequence of M. aeruginosa PCC 7806 (Mic-PCC7806), represented by these contigs, consists of 5,172,804 bases, with an average G+C content of 42%. These values are consistent with those previously determined using thermally denatured DNA . The contigs were annotated using CAAT-Box software and a total of 5,292 predicted protein-coding sequences (CDSs) were validated manually. These CDSs were compared to several protein (Uniprot, COG and 45 cyanobacterial proteomes) and motif databases (Prosite and Pfam).
All the genomes used for the comparative studies described below are listed in the Methods section.
Comparison with other cyanobacterial genomes
A concatenated dataset of large and small subunit rRNA sequences (23S and 16S rRNA) was used to construct a phylogenetic tree including Mic-PCC7806 and 37 other cyanobacterial strains (Figure 1). The tree is congruent with previously published ones based on 16S rRNA sequences [19, 20], but shows higher statistical support at most nodes (especially internal ones), probably due to the larger number of positions used. The strains of the genus Microcystis form a well-supported group (BV of 853‰) with Synechocystis sp. (Syn-PCC6803), Crocosphaera watsonii (Cwa-WH8501) and Cyanothece sp. (Cth-CCY0110 and Cth-ATCC51142). Within this group, Microcystis is most closely related to Syn-PCC6803 (BV of 990‰).
The Mic-PCC7806 genome was compared to the recently publicly available genome of Microcystis aeruginosa strain NIES-843 (Mic-NIES843) . Although the average similarity between the orthologous genes is 94%, their comparison emphasizes that the two genomes largely differ both in length and gene composition (Table 1). Indeed, the Mic-NIES843 genome is 0.6 Mb longer than that of Mic-PCC7806. Moreover, the two genomes display a high number of strain-specific genes (838 for Mic-PCC7806 and 1760 for Mic-NIES843). Interestingly, most of these genes are absent from 44 other cyanobacterial complete genomes suggesting that they have recently been acquired in each of the two Microcystis strains independently. Although the two genomes contain the same proportions of large DNA repeats (~12%, see below), their distribution and size partly differ since Mic-PCC7806 contains 48 repeats longer than 3,000 bases for only 11 in Mic-NIES843. The comparison of the location of similar genes in the largest contig of the Mic-PCC7806 assembly (contig328) and in the Mic-NIES843 genome shows numerous genomic rearrangements (see Additional file 1). These rearrangements, probably facilitated by the presence of large repeats, render the Mic-NIES843 genome of little help for the finishing of the assembly process of the Mic-PCC7806 genome sequence.
The 5292 CDSs of the Mic-PCC7806 genome were also compared to the proteomes of 44 strains representing the diversity of the cyanobacterial lineages (all publicly available genomes excluding Mic-NIES843). The distribution of the best High Scoring Pairs (HSPs) found using Blastall software indicates a high similarity between the proteome of Mic-PCC7806 and a group of three strains Cth-ATCC51142, Cth-CCY0110 and Cwa-WH8501 (Table 2). This is puzzling, since Mic-PCC7806 is closer to Syn-PCC6803 than to this group in the 23S-16S phylogeny (Figure 1). In order to exclude possible bias introduced by uneven distribution of CDSs in these genomes, we analyzed only the orthologs shared by three of these genomes, Mic-PCC7806, Syn-PCC6803 and Cwa-WH8501. Based on BiDirectional Best Hit (BDBH) analyses, 1789 CDSs of the Mic-PCC7806 genome were found to correspond to putative orthologs in Cwa-WH8501 and Syn-PCC6803. The mean Blast score of these CDSs was 381 for the comparison between Mic-PCC7806 and Cwa-WH8501, and only 366 for Mic-PCC7806 versus Syn-PCC6803. The distribution curve of all the Blast scores (see Additional file 2) showed that the Mic-PCC7806 genome was more closely related to Cwa-WH8501 than to Syn-PCC6803 for all score values considered. The absence of congruence between the results obtained with rDNA sequences and the core proteins means that additional data sets for other members of these three cyanobacterial genera are required. Nevertheless, the results obtained by comparing all the orthologous genes shared by Mic-PCC7806 (freshwater strain) and Cwa-WH8501 (marine strain) are consistent with the fact that freshwater and marine cyanobacteria are interspersed in global 16S rDNA phylogenetic trees .
Three distinct groups of proteins were identified on the basis of Blastp analyses of the 5,292 CDSs of Mic-PCC7806, with a selection of 15 other cyanobacterial genomes displaying at least 1% of best Blastp hits with Mic-PCC7806 (Table 2). The composition of these groups largely depends on the threshold chosen to consider that two proteins are similar. Without an obvious breakpoint in the distribution of protein similarities between different genomes (see Additional file 2), we arbitrarily chose a threshold of 40% of similarity, considering that below this value two proteins do not share the same function. The three groups are as follows:
- The "maeru40" group included 764 CDSs (14.4%) specific to the Mic-PCC7806 genome and not found in the 15 selected genomes; 438 (8.3%) of them have no homolog in the uniprot database;
- The "core40" group comprised 652 proteins (12.3%) sharing significant Blastp scores with at least one CDS in each of the 15 other genomes tested;
- The last group, designated "other40", consisted of 3,876 CDSs (73%) sharing significant Blastp scores with CDSs in only some of the other 15 genomes tested.
The small percentage of CDSs in the core40 group reflects the wide diversity of the cyanobacterial genomes analyzed. In the other40 group, the distribution of the Mic-PCC7806 CDSs among the tested genomes matches their phylogenetic distances based on 23S-16S rDNA sequences. For example, in this group, 10% of the CDSs were present in all the genomes, apart from that of Gvi-PCC7421, which is the most distant phylogenetically (Figure 1). Moreover, the four closest genomes to Mic-PCC7806 (Syn-PCC6803 and the group including Cwa-WH8501, Cth-CCY0110 and Cth-ATCC51142) appear to have the same percentage (2%) of CDSs, shared only with Mic-PCC7806.
Plasticity of the genome of M. aeruginosa PCC 7806
Large number of long repeated sequences
The Mic-PCC7806 genome includes a very large number of DNA sequences containing more than 1000 bases that are repeated at least twice in the genome with more than 90% identity. A comparative analysis of all the cyanobacterial genome sequences available in databases showed that Mic-PCC7806, Mic-NIES843 and Cwa-WH8501 are particularly rich in such DNA repeats. Indeed, they account for 11.7%, 11.7% and 19.8% of the total DNA length, respectively (Figure 2). The cumulative size of the DNA repeated sequences is not strictly a function of genome length as Mic-PCC7806 and Cwa-WH8501 genomes have the highest percentage of DNA repeats, but are of intermediate size relative to the other cyanobacterial genomes (see Additional file 3). In the Mic-PCC7806 genome, 1346 CDSs (25%) are located within these DNA repeats. Among these CDSs, only 256 and 92 belong to the maeru40 and core40 groups, respectively. Most of the CDSs of the core40 group correspond to orthologs that are not located within DNA repeats in other cyanobacterial genomes. This implies that over the course of evolution, resident genes were probably captured by genetic mobile elements. A large number of CDSs (362) are very similar to transposases from the COG database, and 93% of them are located within long DNA repeated sequences. At least 46 transposases correspond to ISMae1A/2/3/4 that had previously been characterized in strain PCC7806 , but a large majority of the other transposases cannot be clearly associated with any known insertion sequence (only 17 are associated to IS30, 7 to IS1 and 3 to IS5). The genome of Cwa-WH8501 also contains numerous putative transposases. One third of them are associated to IS5, but none to IS30; the DNA repeated sequences are therefore different in each genome, and cannot account for the close phylogenetic relationship between these two strains.
Synteny of cyanobacterial genomes
Although Mic-PCC7806 and Mic-NIES843 are very closely related strains (Figure 1), their genomes contain a high number of rearrangements. Moreover, an unexpectedly low level of synteny was also observed between the Microcystis strains and two close relatives, Cwa-WH8501 and Syn-PCC6803 (68% mean CDS similarity). Since the same observation was made for all the cyanobacterial genomes tested, we compared the dynamics of these genomes using a large set of other bacterial genomes chosen on the basis of their sizes and phylogenetic distances. To this end, a synteny score was calculated for a number of genome pairs (see Methods), and then compared to their evolutionary distance based on the 23S-16S rDNA tree. This analysis showed that the synteny scores for cyanobacterial genomes were significantly lower than those obtained for pairs of non-cyanobacterial genomes with similar genome lengths and 23S-16S phylogenetic distances (Figure 3). Similar results were obtained for all the cyanobacterial genomes tested. This means that the low synteny scores observed cannot be related to the long DNA repeated sequences, which occur only in the Mic-PCC7806 and Cwa-WH8501 genomes. These results are in agreement with those of Fang et al. , who showed that both persistent and rare genes are significantly clustered in most of the 169 bacterial genomes analyzed. However, in a minority subset of bacterial genomes that includes the cyanobacteria, persistent genes were found to be fairly uniformly distributed throughout the genome.
Interestingly, only 8 clusters with at least 4 CDSs remain syntenic in the genomes of Mic-PCC7806, Cwa-WH8501 and Syn-PCC6803. Four of these clusters correspond to ribosomal proteins. The other clusters are shown in Table 3. Considering the very low level of synteny between cyanobacterial genomes, it is likely that these specific clusters have been subjected to strong positive selection pressure and may play essential roles in these cyanobacteria. Some of these clusters are clearly linked to a specific biological function, such as the transport of phosphate (see Additional file 4) , while others consist of conserved proteins with unknown functions. One can thus speculate that these proteins may be involved in the same biological pathway as their close neighbors.
Four groups can clearly be identified among the cyanobacterial genomes studied on the basis of their intergenic distances (Figure 4). The first consists solely of the genome of Ter-IMS101, which harbors exceptionally long intergenic regions. To the best of our knowledge, no data has been published on this genome, which makes it impossible to rule out the possibility that these regions result from the poor quality of the sequence or the syntaxic annotation. The second group includes the genome of Mic-PCC7806 and, among others, those of Cwa-WH8501 and Syn-PCC6803 which have a high proportion of intergenic sequences around 300 bases long; in the case of the Mic-PCC7806 genome, less than 35% of intergenic sequences are shorter than 100 bases. The third group comprises the genomes of Syn-PCC7942, Tel-BP1 and Gvi-PCC7421, which have short intergenic regions, similar in size to those found in a number of other bacterial genomes (see Additional file 5). The fourth group includes some members of the Prochlorococcus genus that have very small genomes with short or no intergenic regions.
The mean length of the intergenic sequences seems to be linked to the genome size of the cyanobacterium, except for the genome of Syn-PCC6803, which is smaller (3.6 Mb) than that of Gvi-PCC7421 (4.6 Mb), but harbors longer intergenic sequences. Although the role of long intergenic sequences in most cyanobacterial genomes remains unclear, we can surmise that they might be involved in the modulation of gene expression, which would allow cells to acclimate to rapid environmental changes.
Cluster of atypical genes
In order to explore the plasticity of the Mic-PCC7806 genome further, the number of CDSs with an atypical dinucleotide composition was determined using a one-order Markov chain-based methodology . This method can identify genes that may have been acquired recently by lateral transfers. In the Mic-PCC7806 genome, a total of 1971 atypical genes were found, including 1402 within 159 clusters of atypical genes (CAGs) that probably correspond to recently acquired foreign genomic elements (Table 4). As expected, more than 98% of Mic-PCC7806 genes belonging to the core40 group were not in CAGs, and 31% of the atypical genes were in the maeru40 group. Moreover, a high percentage (80%) of the transposase genes were in CAGs (16% of the genes present in CAGs encode putative transposases). Compared to seven other cyanobacterial genomes, those of Mic-PCC7806 and Mic-NIES843 harbor the highest percentages of atypical genes (37%) and CAGs (34% and 36%, respectively). These findings may indicate that the Microcystis genomes contain a higher proportion of genes recently acquired by lateral transfers than the other genomes studied.
Putative restriction and modification systems
Blast searches for restriction enzymes and examination of genes surrounding DNA methylases, identified 21 potential restriction enzymes (see Additional file 6), seven of which were found to be co-localized with putative methylases (see Additional file 7) in the Mic-PCC7806 genome. The Mic-NIES843 genome also contains a high number (at least 17) of putative restriction enzymes . Blast searches revealed that 14 restriction enzymes are common to both genomes. In contrast, seven and eight restriction enzymes seem specific to Mic-PCC7806 and Mic-NIES843, respectively. The Microcystis aeruginosa strains might thus constitute a rich source of novel restriction enzymes potentially useful in biotechnology. According to Zhao et al. , filamentous cyanobacteria (Anabaena, Spirulina and Nostoc strains) contain more restriction and modification genes than unicellular cyanobacteria (Synechocystis, Synechococcus and Prochlorococcus strains). Based on COG annotations, at least as many restriction-modification genes were found in Mic-PCC7806, Mic-NIES843 and Cwa-WH8501 as in filamentous cyanobacteria. Thus, rather than corresponding to a difference between filamentous and unicellular cyanobacteria, the restriction-modification gene content of Microcystis aeruginosa may reflect the potential exposure of the cells to high concentrations of foreign DNA due to the presence of numerous other bacterial cells or viruses associated with Microcystis colonies . This exposure to foreign DNA is also consistent with the high number of CAGs putatively acquired by lateral transfers. Whether such a hypothesis might also hold true for planktonic cyanobacteria of the genus Crocosphaera remains an open question.
In bacterial genomes containing a high number of genes for restriction enzymes, short palindromic sequences corresponding to the target sites of these enzymes may be under-represented . Since the genomes of Microcystis aeruginosa and Cwa-WH8501 contain a very high number of putative restriction enzymes, there should be a number of under-represented short sequences that correspond to restriction sites. To test this hypothesis, the number of occurrences of each 6-mer was counted, and a frequency distribution calculated for Mic-PCC7806, Mic-NIES843, Cwa-WH8501 and Syn-PCC6803 (Table 5). The under-represented sites in the three first genomes were not found in Syn-PCC6803, a genome devoid of restriction enzymes , supporting the idea that these rare 6-mers could indeed correspond to restriction enzyme sites. In total, there are 4096 possible 6-mers, 1.5% of which are palindromes. Fifty-one percent of the rarest 1% of 6-mers in the Mic-PCC7806 genome are palindromes (see Additional file 8). Palindromes are thus over-represented among the rarest 6-mers, further supporting the hypothesis that they could correspond to sites cut by restriction enzymes. The identity of the rarest 1% of 6-mers in the Mic-PCC7806 genome was compared to known restriction sites in other organisms as identified by New England Biolabs . We found that 20 of the 41 sites corresponded to sites cut by restriction enzymes in other organisms.
A novel DNA modification system was discovered recently in the Gram-positive bacterium Streptomyces lividans 66 . This system results in the degradation of DNA in vitro by oxidative, double-stranded, site-specific cleavage during electrophoresis, and is determined by a cluster of five genes (dndA-B-C-D-E). The dnd gene products incorporate sulfur into the DNA backbone as a sequence-selective, stereospecific phosphorothioate modification . According to He et al. , the resistance of phosphorothiate linkages to a variety of nuclease activities, and the site specific nature of such a modification suggest that phosphorothioates could have a role comparable to that of DNA methylation in protection against nucleases. Although the presence of dndB homologs is not clear in the genomes of cyanobacteria, the rest of the cluster was found in several of them including Mic-PCC7806 (see Additional file 9). Despite the low level of synteny in cyanobacterial genomes (see above), the dndC-D-E genes are still clustered.
Unraveling genetic features related to the ecophysiology of M. aeruginosa PCC 7806
Life cycle, colony formation and floatation
During the overwintering benthic phase of their life cycle, Microcystis colonies withstand long periods of darkness. A fermentation pathway has been proposed based on biochemical data . All the genes coding for the enzymes required for the various steps in this pathway have been identified in the genome sequence (see Additional file 10). During the benthic phase, Microcystis colonies are exposed to lower temperature and higher pressure. In this respect, it is interesting to note the presence of a gene (mic5251) coding for a protein similar to Hik33 that perceives osmotic stress and cold stress in Syn-PCC6803 . Another gene, mic5237, is similar to the Ana-PCC7120 orrA gene whose product is involved in osmoregulation . A genomic island carrying actM and pfnM, two genes that encode eukaryotic-like proteins, actin and profilin (an actin cognate binding partner), respectively, have been discovered in the Mic-PCC7806 genome. As shown by Guljamow et al. , this eukaryotic-like actin forms a shell-like structure that could strengthen cell resistance to hydrostatic and osmotic pressures. Interestingly, these genes are only present in Microcystis cells that inhabit the Braakman water reservoir (The Netherlands), which was cut off from the sea in the 20th century, and from which the Mic-PCC7806 strain was originally isolated.
Although several different M. aeruginosa morphotypes have been described , little is known about their colony formation. The genome sequence of strain Mic-PCC7806 revealed a gene coding for a lectin (mvn; mic3128), which binds specifically to a sugar moiety present on the surface of Mic-PCC7806 cells, and a binding partner has been identified in the lipoplysaccharide fraction . A functional correlation between the potent toxin microcystin and this lectin has been demonstrated, with possible implications for the formation of colonial aggregates that are characteristic of different Microcystis morphotypes. Another protein, MrpC (m icrocystin-r elated p rotein C), has been shown to be a potential target of an O-glycosyltransferase of the SPINDLY family . In situ, this protein accumulates at the cell surface, and is involved in cellular interactions. Microcystins may therefore have an impact on the aggregation of Microcystis cells, which is very important for the competitive advantage of these organisms over other phytoplankton species. Mvn and MrpC are predominantly encoded in toxic strains [ and E. Dittmann, unpublished data], but not in the genome of Mic-NIES843. The latter strain may thus represent an ecotype that differs from Mic-PCC7806 in the characteristics of the cell surface. Genes coding for a Ser/Thr kinase (mic0129) and a Ser/Thr phosphatase of the PPP family (mic4622) are found within two clusters that may be involved in cell wall synthesis. Mic-PCC7806 also has two genes that encode Wzc-like protein Tyr kinases (mic2086 and mic1089) and three genes coding for Wzb-like protein Tyr phosphatases (mic3515, mic3588 and mic6566). In E. coli, the function of these systems is known to be related to the synthesis of the cell wall and polysaccharides . These kinases/phosphatases could potentially be involved in colony formation. Colony migration depends not only on the cell ballast resulting from the accumulation of photosynthates and the size of the colonies, but also on the synthesis of gas vesicles (GV), intracellular structures providing cells with buoyancy . The Mic-PCC7806 genome carries a cluster of 12 genes required for GV synthesis, two of which, gvp V and gvp W, are novel . The mic1271 and mic1270 genes are highly similar to the genes coding for a light-regulated two-component system in Syn-PCC6803. This system, which consists of a cyanobacterial phytochrome (Cph1) and its response regulator (Rcp1), has been proposed to play a role in the control of processes required for the adaptation from light to dark conditions and vice-versa . Moreover, all the genes involved in circadian rhythm  are present in Mic-PCC7806 (see Additional file 11). Whether day-night cycles and the timing of vertical migration of Microcystis colonies in the water column are controlled by this phytochrome and by the circadian clock mechanism would be worth being tested.
In natural populations of Microcystis, oxidative stress was shown to induce programmed cell death (PCD) . Accordingly, 5 putative eukaryotic caspase-like genes were identified by PSI-Blast in the genome of strain Mic-PCC7806. Three of them (Mic0980, Mic3930 and Mic4051) showed best similarity with Mic-NIES843 proteins that lack caspase-like motifs. Consequently, these three proteins are likely involved in other functions than PCD. In contrast, the Mic1068 protein showed similarity in the caspase-like region with one protein of Mic-NIES843 (MAE24870). The last caspase-like protein of Mic-PCC7806 (Mic5406) is strain-specific. Both mic1068 and mic5406 are expressed, and a cross-reaction with human caspase-3 polyclonal antisera was observed indicating that the proteins are synthesized (data not shown). Alignment of the regions containing the conserved caspase domains of Mic1068, Mic5406, MAE24870 and a yeast metacaspase shows that the Histidine-Cysteine catalytic diad of the key functionnal regions of the capases is conserved (see Additional file 12). PCD might thus be triggered when Microcystis cells are exposed to severe environmental stress conditions, leading to the rapid decline of blooms, as has been suggested by Berman-Frank et al. in the case of Ter-IMS101 . Mic-PCC7806 and Mic-NIES843 are the only unicellular cyanobacteria known to have genes coding for HstK-like kinases (mic1879 and mic1015), proteins characterized by the presence of both His and Ser/Thr kinase domains [48, 49]. Some of these kinases are implicated in either the iron homeostasis/oxidative stress response or in the differentiation of N2-fixing cells in filamentous cyanobacteria [[48, 49] and C-C Zhang, unpublished data]. Cell differentiation does not occur in M. aeruginosa, but it would be interesting to test whether these HstK-like protein kinases are involved in iron homeostasis and/or in the control of programmed cell death in response to oxidative stress. It has been proposed that the methionine recycling pathway may contribute to preventing oxidative stress in Bacillus subtilis [50, 51]. Interestingly, all the genes involved in this pathway are present in the Mic-PCC7806 genome (see Additional file 13). One of these genes, mtnW (rbcLIV), encodes a 2,3-diketo-5-methylthiopentyl-1-phosphate enolase that has been identified in all the Microcystis strains tested including Mic-NIES843 [21, 52], but not in other cyanobacteria for which the genome sequences are available, except Lae-PCC8106 (accession n° ZP_01618990) and Cth-PCC8801 (accession n° ZP_02940034). The putative methionine recycling pathway may thus have a specific role related to the lifestyle or ecological niches inhabited by members of the genera Microcystis, Lyngbya and Cyanothece.
Genetic potential for the production of secondary metabolites
Cyanobacteria are known as prolific producers of natural products, in particular of the nonribosomal peptide and polyketide classes [15, 53]. However, the potential to produce complex secondary metabolites largely varies among the cyanobacterial genera and species, and even among individual strains. Remarkably, the genomes of Mic-PCC7806, Mic-NIES843 and Cwa-WH8501 differ from unicellular cyanobacteria of other genera in that they contain a large number of genes that encode nonribosomal peptide synthetases (NRPS) and polyketide synthases (PKS). Interestingly, such genes in Mic-PCC7806 outnumber those found in Mic-NIES843 and Cwa-WH8501 (Table 6). Apart from the terrestrial filamentous strain Npu-PCC73102, Mic-PCC7806 devotes the largest percentage of its genome (~3.5%) to secondary metabolite production (Table 6) .
The strain Mic-PCC7806 is known to produce two isoforms of microcystin . The corresponding genes in the bi-directional mcyA-J gene cluster encoding NRPS, PKS and tailoring enzymes [56, 57] could be re-assigned during the genome sequencing project (Figure 5). Genes for cyanopeptolin biosynthesis (mcn cluster) could be assigned based on the amino acid specificities of the substrate-activating domains of a second NRPS gene cluster that was congruent with the amino acid moieties contained in the cyanopeptolin structure  (Figure 5). The mcn genes of Mic-PCC7806 display some similarity to the anabaenopeptilide genes of Anabaena strain 90  and to the cyanopeptolin genes of Microcystis wesenbergii . In addition, the genome of Mic-PCC7806 harbors three NRPS and PKS gene clusters (Figure 5). One of the clusters displays some similarity to the cluster involved in the production of the protease inhibitor aeruginoside in Planktothrix agardhii Cya 126 . The genomic data therefore clearly indicate that strain Mic-PCC7806 might be capable of producing a variant of aeruginosin (Figure 5).
The two remaining PKS I gene clusters do not show significant similarity to any known cyanobacterial biosynthetic gene clusters, and may be involved in the production of hitherto unknown compounds (Figure 5 and Table 6). The first gene cluster encodes an iterative PKS I that is similar in both architecture and sequence to the PksE of various actinobacteria, and is accompanied by several tailoring enzymes including three halogenases. The actinobacterial enzyme is involved in the biosynthesis of enedyine type antitumor antibiotics . The second PKS gene cluster encodes a modular PKS I complex accompanied by several putative tailoring enzymes, and a PKS III type enzyme that is capable of synthesizing compounds of the chalcone/stilbene family. These biosynthetic enzymes are widespread in plants but have only recently been discovered in bacteria . A comparison of the biosynthetic potential of Mic-PCC7806 and Mic-NIES843 reveals that three of the large NRPS/PKS complexes, namely those dedicated to microcystin, cyanopeptolin and aeruginosin production, are encoded on both genomes, whereas some other gene clusters are not shared by both genomes. The biosynthetic versatility of members of the genus Microcystis may thus be larger than expected, since the two strains selected for genome sequencing have similar chemotypes. Beside the NRPS and PKS encoding genes, the genome of Mic-PCC7806 contains a gene cluster similar to the patellamide genes that were recently detected in symbiotic cyanobacterial strains of ascidians . Patellamides are a family of cyclic peptides generated from a ribosomally-synthesized precursor. Mic-PCC7806 is the first freshwater cyanobacterium showing the capability to produce patellamide-like peptides. A peptide with striking similarity to the patellamides, microcyclamide, has been reported in M. aeruginosa strain NIES-298 . Chemical analyses have revealed that the gene cluster discovered in Mic-PCC7806 is indeed dedicated to the production of a microcyclamide-type compound . The genome of Mic-PCC7806 could attract further attention, as it also contains gene clusters comprising unique features that have yet to be characterized and which may well produce so-far unidentified natural substances.
Transporter genes are commonly found in the immediate vicinity of the secondary metabolite biosynthetic genes. These secondary metabolites may therefore at least partly function at the surface of Microcystis cells, in the colony-surrounding sheath or in their planktonic environment. Gene clusters involved in the synthesis of secondary metabolites are frequently associated with genes that confer resistance to these metabolites, which would otherwise be toxic to the cells producing them. In Mic-PCC7806, only the transport system associated with the uncharacterized PKS I/PKS III hybrid compound (Figure 5) shows any similarity to typical efflux transporters that potentially confer self-resistance. The compound produced could therefore have an allelopathic or antibacterial role in the environment .
Among bacteria, members of the genus Microcystis have a particularly high potential for the production of complex secondary metabolites, although this is lower than that of some actinobacterial and myxobacterial genomes that have been shown to devote up to 10% of their coding capacity to the production of secondary metabolites . Genomics has already been useful to the study of secondary metabolites, and has restored natural product research as a major field of pharmaceutical research . Analysis of the Mic-PCC7806 genome has revealed striking novel biosynthetic features that might help to explain the ecological impact of these compounds, as well as guide the search for novel metabolites of biotechnological importance.
Data mining of the genome sequence of Mic-PCC7806 has also shed light on genes that are of importance for the colonial life style and survival of this cyanobacterium in its natural habitat, either during the benthic phase or when it forms blooms on the surface of the water. One of the most intriguing features of this genome is its exceptional plasticity, characterized by a very large number of long repeated sequences, and genes encoding transposases and putative restriction enzymes. These biological entities may generate deletions, duplications, conversions, and rearrangements in the chromosome . One illustration of these changes is the marked loss of synteny between this genome and other cyanobacterial genomes. In addition, the presence of a large number of clustered atypical genes in the genome of Mic-PCC7806 suggests that frequent gene acquisition events by lateral transfers have occurred.
Genome plasticity in prokaryotes is often considered to be an adaptive strategy allowing microorganisms to promote diversification in a way similar to sexual reproduction in eukaryotic organisms. However, genomic rearrangements can also impede the co-expression of genes  and disrupt gene dosage effects . The resulting trade-off between gene conservation and rearrangement in the chromosome depends on various factors and processes linked to the ecophysiology of the microorganisms. The cost of chromosome rearrangements may be greater for fast-growing bacteria, than for slow-growing ones such as cyanobacteria . The relative importance of the process of gene co-expression in cyanobacteria is more difficult to evaluate. However, it is worth noting that some of the eight syntenic clusters found in Mic-PCC7806 concern transport systems for nutrients, such as phosphate, which is often the limiting factor in marine and freshwater ecosystems.
Although Syn-PCC6803, Cwa-WH8501 and Mic-PCC7806 are phylogenetically closely related, only the last two strains have highly plastic genomes containing high proportions of long DNA repeats and transposase genes. No obvious explanation can be deduced from the ecophysiological features of these two strains. Indeed, members of the genus Microcystis are freshwater colonial cyanobacteria that proliferate in eutrophic ecosystems (e.g. ≤ 2.107cells/ml in ) while the Crocosphaera are marine nitrogen-fixing cyanobacteria living in oligotrophic open oceans (≤ 103 cells/ml ). Microcystis colonies may display chaotic population dynamics, with alternating explosion and crash phases , but to the best of our knowledge, no such data are available for Crocosphaera. Such chaotic population dynamics could explain the widespread occurrence of rearrangements in the Mic-PCC7806 genome, if, as proposed by Helm et al.  for Salmonella serovars, bottlenecks and genetic drifts generally promote the fixation of mildly harmful rearrangements.
More genome sequences of members of the Microcystis and Crocosphaera genera are required to clarify the molecular basis of their genome plasticity, at both the intergeneric and intraspecies levels. This will also provide a deeper understanding of the evolutionary significance of this mode of adaptation to the environment. The ongoing sequencing of such genomes should make it possible to reach this goal in the near future. More generally, large cyanobacterial genomes constitute excellent model systems for studying genome dynamics and the mechanism(s) by which some gene clusters may escape rearrangement and retain the same physical organization in several different lineages.
Pairs of cyanobacterial genomes used in Figure 3
Other bacterial strains used in Figure 3 (genome accession number)
Shigella dysenteriae, serovar 1, strain Sd97/Sd197 (CP000034_GR)
Acidovorax avenae subsp. citrulli AAC00-1 (NC_008752)
Agrobacterium tumefaciens str. C58 (NC_003062)
Bacillus subtilis subsp. subtilis str. 168 (NC_000964)
Bordetella parapertussis 12822 (NC_002928)
Escherichia coli APEC O1 (NC_008563)
Enterobacter sp. 638 (NC_009436)
Janthinobacterium sp. Marseille (NC_009659)
Klebsiella pneumoniae subsp. pneumoniae MGH 78578 (CP000647)
Listeria monocytogenes EGD-e (NC_003210)
Methylococcus capsulatus str. Bath (NC_002977)
Ochrobactrum anthropi ATCC 49188 chromosome 1 (NC_009667)
Polaromonas naphthalenivorans CJ2 (NC_008781)
Pseudomonas aeruginosa PA7 (NC_009656)
Pseudomonas fluorescens PfO-1 (NC_007492)
Rhizobium etli CFN 42 (NC_007761)
Rhizobium leguminosarum bv. viciae 3841 (NC_008380)
Rhodobacter sphaeroides ATCC 17025 (NC_009428)
Rhodoferax ferrireducens T118 (NC_007908)
Shewanella loihica PV-4 (NC_009092)
Shewanella oneidensis MR-1 (NC_004347)
Shewanella sp. W3-18-1 (NC_008750)
Shigella boydii Sb227 (NC_007613)
Silicibacter sp. TM1040 (NC_008044)
Yersinia enterocolitica subsp. enterocolitica 8081 (NC_008800)
Yersinia pestis CO92 (NC_003143)
Photorhabdus luminescens subsp. laumondii TTO1 (NC_005126)
DNA preparation and sequencing
The strain Microcystis aeruginosa PCC 7806 (kept in constant culture since its isolation in 1978; Pasteur Culture Collection, Paris, France ) was grown as described . The genome sequence of Mic-PCC7806 was determined by a whole-genome shotgun strategy. Two libraries were generated using genomic DNA extracted with the kit Nucleobond AGX500 (Macherey-Nagel, Hoerdt, France) and shared by nebulization. The first library contained inserts from 1 to 4 kb cloned in pcDNA2.1 (Invitrogen Life Technologies, Carlsbad, CA, USA) and the second included inserts from 5 to 8 kb cloned in the low-copy vector pSYX34 (gift of F. Kunst, Institut Pasteur, Paris, France). A BAC library was constructed into the vector pBeloBAC11 (inserts ≤ 20 kb) (Epicentre, Madison, USA) using spooled DNA extracted as described  and partially hydrolyzed with Hin dIII.
Plasmid DNA purification was performed using the Montage Plasmid Miniprep96 Kit (Millipore, Molsheim, France) or the TempliPhi DNA sequencing template amplification kit (GE Healthcare, Uppsala, Sweden). BAC Miniprep96 Kit (Millipore, Molsheim, France) was used for BAC templates. Sequencing reactions were done, from both ends of DNA inserts, using ABI PRISM BigDye Terminator cycle sequencing ready reactions kit and run on a 3700 Genetic Analyzer (Applied Biosystems, Foster City, CA, USA). The trace file was used with the Phred-Phrap-Consed package to perform the assembly . Sequencing reactions were performed to close gaps, improve coverage and resolve sequence ambiguities using PCR products amplified from genomic DNA or DNA plasmid templates.
A dataset containing a concatenation of the 16S and 23S sequences was aligned by Muscle , and the alignment was manually edited to remove ambiguously aligned positions, giving a final dataset of 4195 nucleotide positions for phylogenetic analysis. From this dataset, a maximum likelihood tree was calculated by Phyml , using the HKY model of nucleotide evolution with an estimation of the transition/transversion ratio, including 4 rates of site heterogeneity, an estimated number of invariable positions, and an estimated alpha shape parameter. The numbers at the nodes correspond to the bootstrap values calculated on 1000 resampled datasets by Phyml.
Syntenic score computation
Ten orthologs located on either side of one pair of putatively orthologous CDS (linked by BDBH) were analyzed. For each pair of orthologous genes located in the proximity of the tested gene and of its ortholog, the synteny score was incremented by 1. Using this method of calculation, two totally syntenic genomes will have a score of 20 attributed to each of their orthologs, whereas two-non syntenic genomes will have a score of 0.
Putative restriction enzymes were identified by Blast searching of known type I and II restriction enzymes against the Mic-PCC7806 genome. Because DNA methylases are more reliably identified by Blast than restriction enzymes, we also identified all methylases, and examined the surrounding genes for potential restriction enzymes.
Detection of atypical CDSs
A first-order Markov model was built based on the dinucleotide composition of the core genes of a group of 8 selected cyanobacterial genomes (Table 4), identified by bi-directional best hits using BLASTp (bitscore of 30% against itself). This Markov model takes into account the Markov probability matrix of the core genes to analyse whether the composition of the CDS under study is "atypical", using the formula described in . For each CDS, the model calculates an index that represents the likelihood that CDS will have a dinucleotide composition compatible with that of the core genes. In order to assess significance cutoffs, we applied the following statistics : for each gene analyzed, one million random sequences were generated based on the Markov model probability matrix of the core genes, and the Markov index was calculated for each of these random sequences. The results were then analyzed by a one-tailed test with cut-offs of 0.1%. The cut-off was defined after several in silico horizontal gene transfer simulations, during which random genes from different genomes were introduced artificially into the genome sequences under study. The optimal threshold (0.1%) was defined for all the genomes of the group as the value at which the model had the highest detection of the in silico introduced genes (true positives), and the lowest detection of core genes (false positives).
Clustering of atypical genes
We defined an initial cluster of at least 4 neighboring atypical genes which was allowed to grow (in both directions) searching for other nearby atypical genes, until regions containing 4 or more non-atypical genes appeared. By this process, a reduced number of less-atypical genes and of normal genes could be included in a larger CAG.
Acaryochloris marina MBIC11017 (embl: CP000828)
Anabaena/Nostoc sp. PCC 7120 (embl: BA000019)
Anabaena variabilis ATCC 29413 (embl: CP000117)
Cyanobium sp. PCC 7001 (gb: 1106012173546)
Cyanothece sp. ATCC 51142 (embl: CP000806)
Cyanothece sp. CCY0110 (gb: 1101676644636–1101676644658)
Crocosphaera watsonii WH8501 (embl: AADV02000100)
Cyanobacteria Yellowstone JA-3-3Ab (embl: CP000239)
Cyanobacteria Yellowstone JA-2-3B'a (embl: CP000240)
Gloeobacter violaceus PCC 7421 (embl: BA000045)
Lyngbya aestuari PCC 8106 (gb: 1099428180563–1099428180584)
Microcoleus chthonoplastes PCC 7420 (gb:1103659003780–1103659003836)
Microcystis aeruginosa NIES-843 (embl: AP009552)
Microcystis aeruginosa PCC 7806 (embl: AM778843–AM778958)
Nodularia spumigena CCY9414 (gb:1099428179735–1099428179797)
Nostoc punctiforme PCC 73102 (kindly provided by J. C. Meeks) 
Prochlorococcus marinus SS120 (embl: AE017126)
Prochlorococcus marinus AS9601 (embl: CP000551)
Prochlorococcus marinus MED4 (embl: BX548174)
Prochlorococcus marinus MIT9211 (embl: AALP01000001)
Prochlorococcus marinus MIT9215 (embl: CP000825)
Prochlorococcus marinus MIT9301 (embl: CP000576)
Prochlorococcus marinus MIT9303 (embl: CP000554)
Prochlorococcus marinus MIT9312 (embl: CP000111)
Prochlorococcus marinus MIT9313 (embl: BX572095)
Prochlorococcus marinus MIT9515 (embl: CP000552)
Prochlorococcus marinus NATL1A (embl: CP000553)
Prochlorococcus marinus NATL2A (embl: CP000095)
Synechococcus sp. BL107 (gb: 1099739244347)
Synechococcus sp. CC9311 (embl: CP000435)
Synechococcus sp. CC9605 (embl: CP000110)
Synechococcus sp. CC9902 (embl: CP000097)
Synechococcus elongatus PCC 6301 (embl: AP008231)
Synechococcus sp. PCC 7002 (embl: CP000951)
Synechococcus sp. PCC 7335 (gb: 1103496006889–1103496006899)
Synechococcus elongatus PCC 7942 (embl: CP000100)
Synechococcus sp. RCC307 (embl: CT978603)
Synechococcus sp. RS9916 (gb: 1100013018508)
Synechococcus sp. RS9917 (gb: 1099465004208)
Synechococcus sp. WH5701 (gb: 1099465003749–1099465003864)
Synechococcus sp. WH7803 (embl: CT971583)
Synechococcus sp. WH7805 (gb: 1099646010155–1099646010157)
Synechococcus sp. WH8102 (gb: BX548020)
Synechocystis sp. PCC 6803 (embl: BA000022)
Thermosynechococcus elongatus BP-1 (embl: BA000039)
Trichodesmium erythreum IMS101 (embl: CP000393)
high scoring segment pair
bidirectional best hit
cluster of atypical gene
bootstrap value. NRPS: nonribosomal peptide synthetase
contig size such that all the larger contigs contain 50% of the bases of the assembly.
Awramic SM: The oldest records of photosynthesis. Photosynth Res. 1992, 33: 75-89. 10.1007/BF00039172.
Dismukes GC, Klimov VV, Baranov SV, Kozlov YN, DasGupta J, Tyryshkin A: The origin of atmospheric oxygen on Earth: The innovation of oxygenic photosynthesis. Proc Natl Acad Sci USA. 2001, 98: 2170-2175. 10.1073/pnas.061514798.
Morán XAG: Annual cycle of picophytoplankton photosynthesis and growth rates in a temperate coastal ecosystem: a major contribution to carbon fluxes. Aquat Microb Ecol. 2007, 49: 267-279. 10.3354/ame01151. [http://www.int-res.com/abstracts/ame/v49/n3/p267-279]
Goericke R, Welschmeyer NA: The marine prochlorophyte Prochlorococcus contributes significantly to phytoplankton biomass and primary production in the Sargasso Sea. Deep Sea Res. 1993, 40: 2283-2294. 10.1016/0967-0637(93)90104-B.
Zehr JP, Waterbury JB, Turner PJ, Montoya JP, Omoregle E, Steward GF, Hansen A, Karl DM: Unicellular cyanobacteria fix N2 in the subtropical North Pacific Ocean. Nature. 2001, 412: 635-638. 10.1038/35088063.
Paerl HW, Fulton RS, Moisander PH, Dyble J: Harmful freshwater algal blooms, with an emphasis on cyanobacteria. ScientificWorldJournal. 2001, 1: 76-113.
Steffensen DA: Economic costs of cyanobacterial blooms. Interagency, International Symposium on Cyanobacterial Harmful Algal Blooms (ISOC-HAB), Advances in Experimental Medicine and Biology. 2008, New York: Springer Verlag, 619: 843-853.
Mur LR: Some aspects of the ecophysiology of cyanobacteria. Ann Microbiol (Paris). 1983, 134B (1): 61-72.
Otsuka S, Suda S, Shibata S, Oyaizu H, Matsumoto S, Watanabe MM: A proposal for the unification of five species of the cyanobacterial genus Microcystis Kützing ex Lemmermann 1907 under the rules of the bacteriological code. Int J Syst Evol Microbiol. 2001, 51: 873-879.
Reynolds CS, Jaworski GHM, Cmiech HA, Leedale GF: On the annual cycle of the blue-green alga Microcystis aeruginosa Kütz. Emend. Elenkin. Philos Trans R Soc Lond B Biol Sci. 1981, 293: 419-477. 10.1098/rstb.1981.0081.
Thomas RH, Walsby AE: Buoyancy regulation in a strain of Microcystis. J Gen Microbiol. 1985, 131: 799-809.
Briand JF, Jacquet S, Bernard C, Humbert JF: Health hazards for terrestrial vertebrates from toxic cyanobacteria in surface water ecosystems. Vet Res. 2003, 34: 361-378. 10.1051/vetres:2003019.
Falconer IR, Humpage AR: Health risk assessment of cyanobacterial (blue-green algal) toxins in drinking water. Int J Environ Res Public Health. 2005, 2: 43-50.
Soares RM, Yuan M, Servaites JC, Delgado A, Magalhaes VF, Hilborn ED, Carmichael WW, Azevedo SLE: Sublethal exposure from microcystins to renal insufficiency patients in Rio de Janeiro, Brazil. Environ Toxicol. 2006, 21: 95-103. 10.1002/tox.20160.
Welker M, von Döhren H: Cyanobacterial peptides – Nature's own combinatorial biosynthesis. FEMS Microbiol Rev. 2006, 30: 530-563. 10.1111/j.1574-6976.2006.00022.x.
Frangeul L, Glaser P, Rusniok C, Buchrieser C, Duchaud E, Dehoux P, Kunst F: CAAT-Box, contigs-assembly and annotation tool-box for genome sequencing projects. Bioinformatics. 2004, 20: 790-797. 10.1093/bioinformatics/btg490.
Gordon D, Desmarais C, Green P: Automated finishing with autofinish. Genome Res. 2001, 11: 614-625. 10.1101/gr.171401.
Rippka R, Herdman M: Catalogue of Strains. Pasteur Culture Collection of Cyanobacterial Strains in Axenic Culture: Catalogue & Taxonomic Handbook. 1992, Paris: Institut Pasteur, I:
Voss B, Gierga G, Axmann IM, Hess WR: A motif-based search in bacterial genomes identifies the ortholog of the small RNA Yfr1 in all lineages of cyanobacteria. BMC Genomics. 2007, 8: 375-385. 10.1186/1471-2164-8-375.
Castenholz RW: Phylum BX. Cyanobacteria – Oxygenic Photosynthetic Bacteria. Bergey's Manual of Systematic Bacteriology. Volume One – The Archaea and the Deeply Branching and Phototrophic Bacteria. Edited by: Boone DR, Castenholz RW, Garrity GM. 2001, Springer Verlag, New York, 473-487. [http://www.springer.com/life+sci/microbiology/book/978-0-387-98771-2]2
Kaneko T, Narajima N, Okamoto S, Suzuki I, Tanabe Y, Tamaoki M, Nakamura Y, Kasai F, Watanabe A, Kawashima K: Complete genomic structure of the bloom-forming toxic cyanobacterium Microcystis aeruginosa NIES-843. DNA Res. 2007, 14: 247-256. 10.1093/dnares/dsm002. [http://dnaresearch.oxfordjournals.org/cgi/content/full/14/6/247]
Mlouka A, Comte K, Tandeau de Marsac N: Mobile DNA elements in the gas vesicle gene cluster of the planktonic cyanobacteria Microcystis aeruginosa. FEMS Microbiol Lett. 2004, 237: 27-34. 10.1111/j.1574-6968.2004.tb09674.x.
Fang G, Rocha EP, Danchin A: Persistence drives gene clustering in bacterial genomes. BMC Genomics. 2008, 9: 4-10.1186/1471-2164-9-4.
Dyhrman ST, Haley ST: Phosphorus scavenging in the unicellular marine diazotroph Crocosphaera watsonii. Appl Environ Microbiol. 2006, 72: 1452-1458. 10.1128/AEM.72.2.1452-1458.2006.
Cortez DQ, Lazcano A, Becerra A: Comparative analysis of methodologies for the detection of horizontally transferred genes: a reassessment of first-order Markov models. In Silico Biol. 2005, 5: 581-592.
Zhao F, Zhang X, Liang C, Wu J, Bao Q, Qin S: Genome-wide analysis of restriction-modification system in unicellular and filamentous cyanobacteria. Physiol Genomics. 2006, 24: 181-190.
Maruyama T, Kato K, Yokoyama A, Tanaka T, Hiraishi A, Park HD: Dynamics of microcystin-degrading bacteria in mucilage of Microcystis. Microb Ecol. 2003, 46: 279-288. 10.1007/s00248-002-3007-7.
Gelfand MS, Koonin EV: Avoidance of palindromic words in bacterial and archaeal genomes: a close connection with restriction enzymes. Nucleic Acids Res. 1997, 25: 2430-2439. 10.1093/nar/25.12.2430.
Scharnagl M, Richter S, Hagemann M: The cyanobacterium Synechocystis sp. strain PCC 6803 expresses a DNA methyltransferase specific for the recognition sequence of the restriction endonuclease Pvu I. J Bacteriol. 1998, 180: 4116-4122.
New England Biolabs, inc. [http://www.neb.com]
Zhou X, He X, Liang J, Li A, Xu T, T K, Helmann JD, Deng Z: A novel DNA modification by sulphur. Mol Microbiol. 2005, 57: 1428-1438. 10.1111/j.1365-2958.2005.04764.x.
Wang L, Chen S, Xu T, Taghizadeh K, Wishnok JS, Zhou X, You D, Deng Z, Dedon P: Phosphorothioation of DNA in bacteria by dnd genes. Nat Chem Biol. 2007, 3: 709-710. 10.1038/nchembio.2007.39.
He X, Ou HY, Yu Q, Zhou X, Wu J, Liang J, Zhang W, Rajakumar K, Deng Z: Analysis of a genomic island housing genes for DNA S-modification system in Streptomyces lividans 66 and its counterparts in other distantly related bacteria. Mol Microbiol. 2007, 65: 1034-1048. 10.1111/j.1365-2958.2007.05846.x.
Moezelaar R, Stal LJ: A comparison of fermentation in the cyanobacterium Microcystis PCC7806 grown under a light/dark cycle and continuous light. Eur J Phycol. 1997, 32: 373-378. [http://www.informaworld.com/smpp/content~content=a714029679~db=all~order=page]
Mikami K, Kanesaki Y, Suzuki I, Murata N: The histidine kinase Hik33 perceives osmotic stress and cold stress in Synechocystis sp. PCC 6803. Mol Microbiol. 2002, 46: 905-915. 10.1046/j.1365-2958.2002.03202.x.
Schwartz SH, Black TA, Jäger K, Panoff J-M, Wolk CP: Regulation of an osmoticum-responsive gene in Anabaena sp. strain PCC 7120. J Bacteriol. 1998, 180: 6332-6337.
Guljamow A, Jenke-Kodama H, Saumweber H, Quillardet P, Frangeul L, Castets AM, Bouchier C, Tandeau de Marsac N, Dittmann E: Horizontal gene transfer of two cytoskeletal elements from a eukaryote to a cyanobacterium. Curr Biol. 2007, 17: R757-R759. 10.1016/j.cub.2007.06.063.
Via-Ordorika L, Fastner J, Kurmayer R, Hisbergues M, Dittmann E, Komarek J, Erhard M, Chorus I: Distribution of microcystin-producing and non-microcystin-producing Microcystis sp. in European freshwater bodies: detection of microcystins and microcystin genes in individual colonies. Syst Appl Microbiol. 2004, 27: 592-602. 10.1078/0723202041748163.
Kehr J-C, Zilliges Y, Springer A, Disney MD, Ratner DD, Bouchier C, Seeberger PH, Tandeau de Marsac N, Dittmann E: A mannan binding lectin is involved in cell-cell attachment in a toxic strain of Microcystis aeruginosa. Mol Microbiol. 2006, 59: 893-906. 10.1111/j.1365-2958.2005.05001.x.
Zilliges Y, Kehr JC, Mikkat S, Bouchier C, Tandeau de Marsac N, Börner T, Dittmann E: An extracellular glycoprotein is implicated in cell-cell contacts in the toxic cyanobacterium Microcystis aeruginosa PCC 7806. J Bacteriol. 2008, 190: 2871-2879. 10.1128/JB.01867-07.
Whitfield C, Paiment A: Biosynthesis and assembly of Group 1 capsular polysaccharides in Escherichia coli and related extracellular polysaccharides in other bacteria. Carbohydr Res. 2003, 338: 2491-2502. 10.1016/j.carres.2003.08.010.
Walsby AE: Gas Vesicles. Microbiol Rev. 1994, 58: 94-144.
Mlouka A, Comte K, Castets AM, Bouchier C, Tandeau de Marsac N: The gas vesicle gene cluster from Microcystis aeruginosa and DNA rearrangements that lead to loss of cell buoyancy. J Bacteriol. 2004, 186: 2355-2365. 10.1128/JB.186.8.2355-2365.2004.
García-Domínguez M, Muro-Pastor MI, Reyes JC, Florencio FJ: Light-dependent regulation of cyanobacterial phytochrome expression. J Bacteriol. 2000, 182: 38-44.
Mackey SR, Golden SS: Winding up the cyanobacterial circadian clock. Trends Microbiol. 2007, 15: 381-388. 10.1016/j.tim.2007.08.005.
Ross C, Santiago-Vazquez L, Paul V: Toxin release in response to oxidative stress and programmed cell death in the cyanobacterium Microcystis aeruginosa. Aquat Toxicol. 2006, 78: 66-73. 10.1016/j.aquatox.2006.02.007.
Berman-Frank I, Bidle K, Haramaty L, Falkowski P: The demise of the marine cyanobacterium, Trichodesmium spp. via an autocatalyzed cell death pathway. Limnol Oceanogr. 2004, 49: 997-1005. [http://aslo.org/lo/toc/vol_49/issue_4/0997.pdf]
Cheng Y, Li JH, Shi L, Wang L, Latifi A, Zhang C-C: A pair of iron-responsive genes encoding protein kinases with a Ser/Thr kinase domain and a His kinase domain are regulated by NtcA in the Cyanobacterium Anabaena sp. strain PCC 7120. J Bacteriol. 2006, 188: 4822-4829. 10.1128/JB.00258-06.
Shi L, Li JH, Cheng Y, Wang L, Chen WL, Zhang C-C: Two genes encoding protein kinases of the HstK family are involved in synthesis of the minor heterocyst-specific glycolipid in the cyanobacterium Anabaena sp. strain PCC 7120. J Bacteriol. 2007, 189: 5075-5081. 10.1128/JB.00323-07.
Sekowska A, Danchin A: The methionine salvage pathway in Bacillus subtilis. BMC Microbiol. 2002, 2: 8-10.1186/1471-2180-2-8.
Sekowska A, Dénervaud V, Ashida H, Michoud K, Haas D, Yokota A, Danchin A: Bacterial variations on the methionine salvage pathway. BMC Microbiol. 2004, 4: 9-10.1186/1471-2180-4-9.
Carré-Mlouka A, Méjean A, Quillardet P, Ashida H, Saito Y, Yokita A, Callebaut I, Sekowska A, Dittmann E, Bouchier C, Tandeau de Marsac N: A new RuBisCO-like protein coexists with a photosynthetic RuBisCO in the planktonic cyanobacteria Microcystis. J Biol Chem. 2006, 281: 24462-24471. 10.1074/jbc.M602973200.
Ehrenreich IM, Waterbury JB, Webb EA: Distribution and diversity of natural product genes in marine and freshwater cyanobacterial cultures and genomes. Appl Environ Microbiol. 2005, 71: 7401-7413. 10.1128/AEM.71.11.7401-7413.2005.
Meeks JC, Elhai J, Thiel T, Potts M, Larimer F, Lamerdin J, Predki P, Atlas R: An overview of the genome of Nostoc punctiforme, a multicellular, symbiotic cyanobacterium. Photosynth Res. 2001, 70: 85-106. 10.1023/A:1013840025518.
Dierstein R, Kaiser I, Weckesser J, Matern U, König WA, Krebber R: Two closely related peptide toxins in axenically grown Microcystis aeruginosa PCC 7806. Syst Appl Microbiol. 1990, 13: 86-91.
Nishizawa T, Ueda A, Asayama M, Fujii K, Harada K-I, Ochi K, Shirai M: Polyketide synthase gene coupled to the peptide synthetase module involved in the biosynthesis of the cyclic heptapeptide microcystin. J Biochem. 2000, 127: 779-789.
Tillett D, Dittmann E, Erhard M, von Döhren H, Börner T, Neilan BA: Structural organization of microcystin biosynthesis in Microcystis aeruginosa PCC7806: an integrated peptide-polyketide synthetase system. Chem Biol. 2000, 7: 753-764. 10.1016/S1074-5521(00)00021-1.
Martin C, Oberer L, Ino T, König WA, Busch M, Weckesser J: Cyanopeptolins, new depsispetides from the cyanobacterium Microcystis sp. PCC 7806. J Antibiot (Tokyo). 1993, 46: 1550-1556.
Rouhiainen L, Paulin L, Suomalainen S, Hyytiäinen H, Buikema W, Haselkorn R, Sivonen K: Genes encoding synthetases of cyclic depsipeptides, anabaenopeptilides, in Anabaena strain 90. Mol Microbiol. 2000, 37: 156-157. 10.1046/j.1365-2958.2000.01982.x.
Tooming-Klunderud A, Rohrlack T, Shalchian-Tabrizi K, Kristensen T, Jakobsen KS: Structural analysis of a non-ribosomal halogenated cyclic peptide and its putative operon from Microcystis : implications for evolution of cyanopeptolins. Microbiology. 2007, 153: 1382-1393. 10.1099/mic.0.2006/001123-0.
Ishida K, Christiansen G, Yoshida WY, Kurmayer R, Welker M, Valls N, Bonjoch J, Hertweck C, Börner T, Hemscheidt T, Dittmann E: Biosynthesis and structure of aeruginoside 126A and 126B, cyanobacterial peptide glycosides bearing a 2-carboxy-6-hydroxyoctahydroindole moiety. Chem Biol. 2007, 14: 565-576. 10.1016/j.chembiol.2007.04.006.
Zazopoulos E, Huang K, Staffa A, Liu W, Bachmann BO, Nonaka K, Ahlert J, Thorson JS, Shen B, Farnet CM: A genomics-guided approach for discovering and expressing cryptic metabolic pathways. Nat Biotechnol. 2003, 21: 187-190. 10.1038/nbt784.
Gross F, Luniak N, Pervola O, Gaitatzis N, Jenke-Kodama H, Gerth K, Gottschalk D, Dittmann E, Muller R: Bacterial type III polyketide synthases: phylogenetic analysis and potential for the production of novel secondary metabolites by heterologous expression in pseudomonads. Arch Microbiol. 2006, 185: 28-38. 10.1007/s00203-005-0059-3.
Schmidt EW, Nelson JT, Rasko DA, Sudek S, Eisen JA, Haygood MG, Ravel J: Patellamide A and C biosynthesis by a microcin-like pathway in Prochloron didemni, the cyanobacterial symbiont of Lissoclinum patella. Proc Natl Acad Sci USA. 2005, 102: 7315-7320. 10.1073/pnas.0501424102.
Ishida K, Nakagawa H, Murakami M: Microcyclamide, a cytotoxic cyclic hexapeptide from the cyanobacterium Microcystis aeruginosa. J Nat Prod. 2000, 63: 1315-1317. 10.1021/np000159p.
Ziemert N, Ishida K, Quillardet P, Bouchier C, Hertweck C, Tandeau de Marsac N, Dittmann E: Microcyclamide biosynthesis in two strains of Microcystis aeruginosa : from structure to genes and vice versa. Appl Environ Microbiol. 2008, 74: 1791-1797. 10.1128/AEM.02392-07.
Schatz D, Keren Y, Vardi A, Sukenik A, Carmeli S, B-örner T, Dittmann E, Kaplan A: Towards clarification of the biological role of microcystins, a family of cyanobacterial toxins. Environ Microbiol. 2007, 9: 965-970. 10.1111/j.1462-2920.2006.01218.x.
Udwary DW, Zeigler L, Asolkar RN, Singan V, Lapidus A, Fenical W, Jensen PR, Moore BS: Genome sequencing reveals complex secondary metabolome in the marine actinomycete Salinispora tropica. Proc Natl Acad Sci USA. 2007, 104: 10376-10381. 10.1073/pnas.0700962104.
Bode HB, Müller R: The impact of bacterial genomics on natural product research. Angew Chem Int Ed Engl. 2005, 44: 6828-6846. 10.1002/anie.200501080.
Rocha EPC: Order and disorder in bacterial genomes. Curr Opin Microbiol. 2004, 7: 519-527. 10.1016/j.mib.2004.08.006.
Korbel JO, Jensen LJ, von Mering C, Bork P: Analysis of genomic context: prediction of functional associations from conserved bidirectionally transcribed gene pairs. Nat Biotechnol. 2004, 22: 911-917. 10.1038/nbt988.
Rocha EPC: The replication-related organization of bacterial genomes. Microbiology. 2004, 150: 1609-1627. 10.1099/mic.0.26974-0.
Lehman PW, Boyer G, Stachwell M, Waller S: The influence of environmental conditions on the seasonnal variations of Microcystis cell density and microcystins concentration in San Fransisco estuary. Hydrobiologia. 2008, 600: 187-204. 10.1007/s10750-007-9231-x.
Montoya JP, Holl CM, Zehr JP, Hansen A, Villareal TA, Capone DG: High rates of N2 fixation by unicellular diazotrophs in the oligotrophic Pacific Ocean. Nature. 2004, 430: 1027-1032. 10.1038/nature02824.
Manage PM, Kawabata Z, Nakano S: Seasonal changes in densities of cyanophage infections to Microcystis aeruginosa in a hypereutrophic pond. Hydrobiologia. 2001, 411: 211-216. 10.1023/A:1003868803832.
Helm RA, Lee AG, Christman HD, Maloy S: Genomic rearrangements at rrn operons in Salmonella. Genetics. 2003, 165: 951-959.
Nostoc punctiforme ATCC 29133 (PCC 73102). [http://genome.jgi-psf.org/finished_microbes/nospu/nospu.home.html]
Sakamoto T, Shirai M, Asayama M, Aida T, Sato A, Tanaka K, Takahashi H, Nakano M: Characteristics of DNA and multiple rpoD homologs of Microcystis (Synechocystis) strains. Int J Syst Bacteriol. 1993, 43: 844-847.
Phred, Phrap, Consed. [http://www.phrap.org]
Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32: 1792-1797. 10.1093/nar/gkh340.
Guindon S, Gascuel O: A simple, fast and accurate algorithm to estimate large phylogenies by maximum-likelihood. Syst Biol. 2003, 52: 696-704. 10.1080/10635150390235520.
Nakamura Y, Itoh T, Matsuda H, Gojobori T: Biased biological functions of horizontally transferred genes in prokaryotic genomes. Nat Genet. 2004, 36: 760-766. 10.1038/ng1381.
We are grateful to S. Cole, P. Glaser and F. Kunst (Pasteur Genopole®), who were involved in the genome sequencing project financed by the Institut Pasteur, the Ministère de l'Education Nationale, de la Recherche et de la Technologie (MENRT), the Consortium national de la recherche en génomique and the Centre National de la Recherche Scientifique (URA 2172). We acknowledge support from the Gordon and Betty Moore Foundation, as part of its Marine Microbial Genome Sequencing Project. We also thank the JCVI software team (leader, S.A. Kravitz) and the JCVI Joint Technology Center (leader, Y.-H. Rogers and sequencing production manager, S. Ferriera). We are grateful to A. Marcel and S. Bun for their contribution to sofware development, and to L. Ma and S. Ferris, for their technical assistance. We would also like to thank A. Danchin and R. Rippka for helpful discussions. M. Ghosh is acknowledged for revising the English version of the manuscript.
LF carried out the bioinformatics studies. PQ and A–MC carried out the molecular genetic studies. AB and A–MC constructed the DNA libraries. CB, AL, PQ and A–MC carried out the sequence of the genome. LF, PQ, A–MC, NTM, ED, J–FH, J–CK, YZ, and NZ annotated the genome. HCPM, SB and NTM analyzed the metabolic pathways. DC carried out the CAG analyses. AT, PQ and LF carried out the enzyme and 6-mer analyses. SG and LF performed the phylogenetic analyses. CCZ, ET and AL analyzed the sensor and regulatory systems. NTM, LF and PQ designed the research. NTM coordinated the study. C–CZ, SG, DC, HCPM and AT drafted the manuscript. LF, NTM, J–FH, PQ and ED wrote the manuscript.
All authors read and approved the final manuscript.
Lionel Frangeul, Philippe Quillardet contributed equally to this work.