DNA methyltransferases of the cyanobacterium Anabaena PCC 7120

Characterization of Limnospira platensis PCC 9108 R-M and CRISPR-Cas systems

The filamentous cyanobacterium Limnospira platensis, formerly known as Arthrospira platensis or spirulina, is one of the most commercially important species of microalgae. Due to its high nutritional value, pharmacological and industrial applications it is extensively cultivated on a large commercial scale. Despite its widespread use, its precise manipulation is still under development due to the lack of effective genetic protocols. Genetic transformation of Limnospira has been attempted but the methods reported have not been generally reproducible in other laboratories. Knowledge of the transformation defense mechanisms is essential for understanding its physiology and for broadening their applications. With the aim to understand more about the genetic defenses of L. platensis, in this work we have identified the restriction-modification and CRISPR-Cas systems and we have cloned and characterized thirteen methylases. In parallel, we have also characterized the methylome and orphan methyltransferases using genome-wide analysis of DNA methylation patterns and RNA-seq. The identification and characterization of these enzymes will be a valuable resource to know how this strain avoids being genetically manipulated and for further genomics studies.

Effects of a Type I RM System on Gene Expression and Glycogen Catabolism in Synechocystis sp. PCC 6803

Increasing evidence has shown that DNA methylation is involved in gene regulation in prokaryotes. However, there have been very limited reports about the role of DNA methylation in regulation of gene expression and physiological functions in cyanobacteria. In Synechocystis sp. PCC 6803, four genes on the plasmid pSYSX are predicted to encode the type I restriction-methylation system, slr6095 and slr6096 for the M subunit, slr6097 for the S subunit and slr6102 for the R subunit. Compared to the wild type, slr6095, slr6096, and slr6097 mutants lacked the GG^m6AN₇TTGG/CCA^m6AN₇TCC methylation in genomic DNA. Transcriptomic analysis indicated that 171 genes were reproducibly up- or down-regulated in all three mutants relative to the wild type. The changed expression of some genes, including sll1356 for glycogen phosphorylase (GlgP), was associated with the loss of GG^m6AN₇TTGG/CCA^m6AN₇TCC methylation in the coding regions or the upstream non-coding sequences. Inactivation of slr6095, slr6096, or slr6097 increased the expression of sll1356 and the GlgP activity but lowered the glycogen content. These results indicated that the DNA methylation by a type I RM system could alter the expression of certain genes and physiological functions in Synechocystis sp. PCC 6803.

Cyanobacterial multi-copy chromosomes and their replication

While the model bacteria Escherichia coli and Bacillus subtilis harbor single chromosomes, which is known as monoploidy, some freshwater cyanobacteria contain multiple chromosome copies per cell throughout their cell cycle, which is known as polyploidy. In the model cyanobacteria Synechococcus elongatus PCC 7942 and Synechocystis sp. PCC 6803, chromosome copy number (ploidy) is regulated in response to growth phase and environmental factors. In S. elongatus 7942, chromosome replication is asynchronous both among cells and chromosomes. Comparative analysis of S. elongatus 7942 and S. sp. 6803 revealed a variety of DNA replication mechanisms. In this review, the current knowledge of ploidy and DNA replication mechanisms in cyanobacteria is summarized together with information on the features common with plant chloroplasts. It is worth noting that the occurrence of polyploidy and its regulation are correlated with certain cyanobacterial lifestyles and are shared between some cyanobacteria and chloroplasts.

ABBREVIATIONS

NGS: next-generation sequencing; Repli-seq: replication sequencing; BrdU: 5-bromo-2'-deoxyuridine; TK: thymidine kinase; GCSI: GC skew index; PET: photosynthetic electron transport; RET: respiration electron transport; Cyt b₆f complex: cytochrome b₆f complex; PQ: plastoquinone; PC: plastocyanin.

Divergent methylation of CRISPR repeats and cas genes in a subtype I-D CRISPR-Cas-system

Background

The presence and activity of CRISPR-Cas defense systems is a hallmark of many prokaryotic microorganisms. Here, the distribution of sequences related to the highly iterated palindrome 1 (HIP1) element and the DNA methylation of CGATCG motifs embedded within HIP1 as a vital part of the CRISPR1 repeat sequence was analyzed in the cyanobacterium Synechocystis sp. PCC 6803. Previously suggested functions of HIP1 include organization of chromosomal structure, DNA recombination or gene regulation, all of which could be relevant in CRISPR-Cas functionality.

Results

The CRISPR1 repeat-spacer array contains more than 50 CGATCG elements that are double-methylated (^5mCG^6mATCG) by the enzymes M.Ssp6803I and M.Ssp6803III. Hence, more than 200 possible methylation events cluster over a stretch of 3600 bp of double-stranded DNA. Bisulfite sequencing showed that these motifs were highly methylated at the ^m5CGATCG positions whereas specific motifs within the CRISPR1 cas genes were hypomethylated suggesting a lowered accessibility for the DNA methylase to these regions. Assays for conjugation and CRISPR1-mediated DNA interference revealed a 50% drop in conjugation efficiency in the mutant lacking the ^5mC methylation of CGATCG motifs, while the highly efficient DNA interference activity was not affected by the lack of ^m5CGATCG DNA-methylation, nor was the capability to differentiate between self and non-self targets based on the protospacer adjacent motifs (PAMs) GTA and GTC versus the non-PAM AGC. A third DNA methylation mediated by M.Ssp6803II modifies the first cytosine in the motif GGCC yielding GG^m4CC. We found a remarkable absence of GGCC motifs and hence the corresponding methylation over an 11 kb stretch encompassing all the cas genes involved in interference and crRNA maturation but not adaptation of the CRISPR1 system.

Conclusions

The lack of GGCC tetranucleotides along the CRISPR1 interference and maturation genes supports the reported hybrid character of subtype I-D CRISPR-Cas systems. We report tight and very high ^5mC methylation of the CRISPR1 repeat sequences. Nevertheless, cells lacking the ^5mC methylation activity were unaffected in their CRISPR1-mediated interference response but the efficiency of conjugation was reduced by 50%. These results point to an unknown role of ^m5CGATCG DNA-methylation marks in conjugation and DNA transformation.

Cytosine N4-Methylation via M.Ssp6803II Is Involved in the Regulation of Transcription, Fine- Tuning of DNA Replication and DNA Repair in the Cyanobacterium Synechocystis sp. PCC 6803

DNA methylation plays a crucial role for gene regulation among eukaryotes, but its regulatory function is less documented in bacteria. In the cyanobacterium Synechocystis sp. PCC 6803 five DNA methyltransferases have been identified. Among them, M.Ssp6803II is responsible for the specific methylation of the first cytosine in the frequently occurring motif GGCC, leading to N4-methylcytosine (GGm4CC). The mutation of the corresponding gene sll0729 led to lowered chlorophyll/phycocyanin ratio and slower growth. Transcriptomics only showed altered expression of sll0470 and sll1526, two genes encoding hypothetical proteins. Moreover, prolonged cultivation revealed instability of the initially obtained phenotype. Colonies with normal pigmentation and wild-type-like growth regularly appeared on agar plates. These colonies represent suppressor mutants, because the sll0729 gene was still completely inactivated and the GGCC sites remained unmethylated. The suppressor strains showed smaller cell size, lowered DNA content per cell, and decreased tolerance against UV compared to wild type. Promoter assays revealed that the transcription of the sll0470 gene was still stimulated in the suppressor clones. Proteomics identified decreased levels of DNA topoisomerase 4 subunit A in suppressor cells. Collectively, these results indicate that GGm4CC methylation is involved in the regulation of gene expression, in the fine-tuning of DNA replication, and DNA repair mechanisms.

Interaction between cyanophage MaMV-DC and eight Microcystis strains, revealed by genetic defense systems

Cyanophage MaMV-DC is a member of Myoviridae that was reported to specifically infect and lyse Microcystis aeruginosa FACHB-524 among 21 selected cyanobacterial strains. We reidentified the infection specificity of MaMV-DC among seven other Microcystis strains of different species. In our experiments, MaMV-DC infected three Microcystis strains but did not form plaque in Microcystis lawns. This indicated that MaMV-DC is at least a genus- rather than strain-specific virus. Cyanophage MaMV-DC genes were transcribed in M. aeruginosa FACHB-524, M. flos-aquae TF09, M. aeruginosa TA09 and M. wesenbergii DW09, and the growth of these Microcystis strains was inhibited by the addition of MaMV-DC. The predicted defense of eight Microcystis strains by CRISPR-Cas systems has shown mixed consistency with the infection experiment results, suggesting other defense or anti-defense systems play roles during infection process. Restriction-modification (RM) system analysis revealed an abundance of four types of RM proteins that may play roles in defense against cyanophages.

Identification of the DNA methyltransferases establishing the methylome of the cyanobacterium Synechocystis sp. PCC 6803

DNA methylation in bacteria is important for defense against foreign DNA, but is also involved in DNA repair, replication, chromosome partitioning, and regulatory processes. Thus, characterization of the underlying DNA methyltransferases in genetically tractable bacteria is of paramount importance. Here, we characterized the methylome and orphan methyltransferases in the model cyanobacterium Synechocystis sp. PCC 6803. Single molecule real-time (SMRT) sequencing revealed four DNA methylation recognition sequences in addition to the previously known motif m5CGATCG, which is recognized by M.Ssp6803I. For three of the new recognition sequences, we identified the responsible methyltransferases. M.Ssp6803II, encoded by the sll0729 gene, modifies GGm4CC, M.Ssp6803III, encoded by slr1803, represents the cyanobacterial dam-like methyltransferase modifying Gm6ATC, and M.Ssp6803V, encoded by slr6095 on plasmid pSYSX, transfers methyl groups to the bipartite motif GGm6AN7TTGG/CCAm6AN7TCC. The remaining methylation recognition sequence GAm6AGGC is probably recognized by methyltransferase M.Ssp6803IV encoded by slr6050. M.Ssp6803III and M.Ssp6803IV were essential for the viability of Synechocystis, while the strains lacking M.Ssp6803I and M.Ssp6803V showed growth similar to the wild type. In contrast, growth was strongly diminished of the Δsll0729 mutant lacking M.Ssp6803II. These data provide the basis for systematic studies on the molecular mechanisms impacted by these methyltransferases.

Transgenerational Epigenetic Inheritance Under Environmental Stress by Genome-Wide DNA Methylation Profiling in Cyanobacterium

Epigenetic modifications such as DNA methylation are well known as connected with many important biological processes. Rapid accumulating evidence shows environmental stress can generate particular defense epigenetic changes across generations in eukaryotes. This transgenerational epigenetic inheritance in animals and plants has gained interest over the last years. Cyanobacteria play very crucial role in the earth, and as the primary producer they can adapt to nearly all diverse environments. However, few knowledge about the genome wide epigenetic information such as methylome information in cyanobacteria, especially under any environment stress, was reported so far. In this study we profiled the genome-wide cytosine methylation from a model cyanobacterium Synechocystis sp. PCC 6803, and explored the possibility of transgenerational epigenetic process in this ancient single-celled prokaryote by comparing the DNA methylomes among normal nitrogen medium cultivation, nitrogen starvation for 72 h and nitrogen recovery for about 12 generations. Our results shows that DNA methylation patterns in nitrogen starvation and nitrogen recovery are much more similar with each other, significantly different from that of the normal nitrogen. This study reveals the difference in global DNA methylation pattern of cyanobacteria between normal and nutrient stress conditions and reports the evidence of transgenerational epigenetic process in cyanobacteria. The results of this study may contribute to a better understanding of epigenetic regulation in prokaryotic adaptation to and survive in the ever changing environment.

A closely-related clade of globally distributed bloom-forming cyanobacteria within the Nostocales

In order to better understand the relationships among current Nostocales cyanobacterial blooms, eight genomes were sequenced from cultured isolates or from environmental metagenomes of recent planktonic Nostocales blooms. Phylogenomic analysis of publicly available sequences placed the new genomes among a group of 15 genomes from four continents in a distinct ADA clade (Anabaena/Dolichospermum/Aphanizomenon) within the Nostocales. This clade contains four species-level groups, two of which include members with both Anabaena-like and Aphanizomenon flos-aquae-like morphology. The genomes contain many repetitive genetic elements and a sizable pangenome, in which ABC-type transporters are highly represented. Alongside common core genes for photosynthesis, the differentiation of N₂-fixing heterocysts, and the uptake and incorporation of the major nutrients P, N and S, we identified several gene pathways in the pangenome that may contribute to niche partitioning. Genes for problematic secondary metabolites-cyanotoxins and taste-and-odor compounds-were sporadically present, as were other polyketide synthase (PKS) and nonribosomal peptide synthetase (NRPS) gene clusters. By contrast, genes predicted to encode the ribosomally generated bacteriocin peptides were found in all genomes.

The highly heterogeneous methylated genomes and diverse restriction-modification systems of bloom-forming Microcystis

The occurrence of harmful Microcystis blooms is increasing in frequency in a myriad of freshwater ecosystems. Despite considerable research pertaining to the cause and nature of these blooms, the molecular mechanisms behind the cosmopolitan distribution and phenotypic diversity in Microcystis are still unclear. We compared the patterns and extent of DNA methylation in three strains of Microcystis, PCC 7806SL, NIES-2549 and FACHB-1757, using Single Molecule Real-Time (SMRT) sequencing technology. Intact restriction-modification (R-M) systems were identified from the genomes of these strains, and from two previously sequenced strains of Microcystis, NIES-843 and TAIHU98. A large number of methylation motifs and R-M genes were identified in these strains, which differ substantially among different strains. Of the 35 motifs identified, eighteen had not previously been reported. Strain NIES-843 contains a larger number of total putative methyltransferase genes than have been reported previously from any bacterial genome. Genomic comparisons reveal that methyltransferases (some partial) may have been acquired from the environment through horizontal gene transfer.

Biosynthesis and Function of Modified Bases in Bacteria and Their Viruses

Naturally occurring modification of the canonical A, G, C, and T bases can be found in the DNA of cellular organisms and viruses from all domains of life. Bacterial viruses (bacteriophages) are a particularly rich but still underexploited source of such modified variant nucleotides. The modifications conserve the coding and base-pairing functions of DNA, but add regulatory and protective functions. In prokaryotes, modified bases appear primarily to be part of an arms race between bacteriophages (and other genomic parasites) and their hosts, although, as in eukaryotes, some modifications have been adapted to convey epigenetic information. The first half of this review catalogs the identification and diversity of DNA modifications found in bacteria and bacteriophages. What is known about the biogenesis, context, and function of these modifications are also described. The second part of the review places these DNA modifications in the context of the arms race between bacteria and bacteriophages. It focuses particularly on the defense and counter-defense strategies that turn on direct recognition of the presence of a modified base. Where modification has been shown to affect other DNA transactions, such as expression and chromosome segregation, that is summarized, with reference to recent reviews.

Premethylation of foreign DNA improves integrative transformation efficiency in Synechocystis sp. strain PCC 6803

Restriction digestion of foreign DNA is one of the key biological barriers against genetic transformation in microorganisms. To establish a high-efficiency transformation protocol in the model cyanobacterium, Synechocystis sp. strain PCC 6803 (Synechocystis 6803), we investigated the effects of premethylation of foreign DNA on the integrative transformation of this strain. In this study, two type II methyltransferase-encoding genes, i.e., sll0729 (gene M) and slr0214 (gene C), were cloned from the chromosome of Synechocystis 6803 and expressed in Escherichia coli harboring an integration plasmid. After premethylation treatment in E. coli, the integration plasmid was extracted and used for transformation of Synechocystis 6803. The results showed that although expression of methyltransferase M had little impact on the transformation of Synechocystis 6803, expression of methyltransferase C resulted in 11- to 161-fold-higher efficiency in the subsequent integrative transformation of Synechocystis 6803. Effective expression of methyltransferase C, which could be achieved by optimizing the 5' untranslated region, was critical to efficient premethylation of the donor DNA and thus high transformation efficiency in Synechocystis 6803. Since premethylating foreign DNA prior to transforming Synechocystis avoids changing the host genetic background, the study thus provides an improved method for high-efficiency integrative transformation of Synechocystis 6803.

Highly Iterated Palindromic Sequences (HIPs) and Their Relationship to DNA Methyltransferases

The sequence GCGATCGC (Highly Iterated Palindrome, HIP1) is commonly found in high frequency in cyanobacterial genomes. An important clue to its function may be the presence of two orphan DNA methyltransferases that recognize internal sequences GATC and CGATCG. An examination of genomes from 97 cyanobacteria, both free-living and obligate symbionts, showed that there are exceptional cases in which HIP1 is at a low frequency or nearly absent. In some of these cases, it appears to have been replaced by a different GC-rich palindromic sequence, alternate HIPs. When HIP1 is at a high frequency, GATC- and CGATCG-specific methyltransferases are generally present in the genome. When an alternate HIP is at high frequency, a methyltransferase specific for that sequence is present. The pattern of 1-nt deviations from HIP1 sequences is biased towards the first and last nucleotides, i.e., those distinguish CGATCG from HIP1. Taken together, the results point to a role of DNA methylation in the creation or functioning of HIP sites. A model is presented that postulates the existence of a GmeC-dependent mismatch repair system whose activity creates and maintains HIP sequences.

Methylation by a unique α-class N4-cytosine methyltransferase is required for DNA transformation of Caldicellulosiruptor bescii DSM6725

Thermophilic microorganisms capable of using complex substrates offer special advantages for the conversion of lignocellulosic biomass to biofuels and bioproducts. Members of the gram-positive bacterial genus Caldicellulosiruptor are anaerobic thermophiles with optimum growth temperatures between 65°C and 78°C and are the most thermophilic cellulolytic organisms known. In fact, they efficiently use biomass non-pretreated as their sole carbon source and in successive rounds of application digest 70% of total switchgrass substrate. The ability to genetically manipulate these organisms is a prerequisite to engineering them for use in conversion of these complex substrates to products of interest as well as identifying gene products critical for their ability to utilize non-pretreated biomass. Here, we report the first example of DNA transformation of a member of this genus, C. bescii. We show that restriction of DNA is a major barrier to transformation (in this case apparently absolute) and that methylation with an endogenous unique α-class N4-Cytosine methyltransferase is required for transformation of DNA isolated from E. coli. The use of modified DNA leads to the development of an efficient and reproducible method for DNA transformation and the combined frequencies of transformation and recombination allow marker replacement between non-replicating plasmids and chromosomal genes providing the basis for rapid and efficient methods of genetic manipulation.

Characterization of type II and III restriction-modification systems from Bacillus cereus strains ATCC 10987 and ATCC 14579

The genomes of two Bacillus cereus strains (ATCC 10987 and ATCC 14579) have been sequenced. Here, we report the specificities of type II/III restriction (R) and modification (M) enzymes. Found in the ATCC 10987 strain, BceSI is a restriction endonuclease (REase) with the recognition and cut site CGAAG 24-25/27-28. BceSII is an isoschizomer of AvaII (G/GWCC). BceSIII cleaves at ACGGC 12/14. The BceSIII C terminus resembles the catalytic domains of AlwI, MlyI, and Nt.BstNBI. BceSIV is composed of two subunits and cleaves on both sides of GCWGC. BceSIV activity is strongly stimulated by the addition of cofactor ATP or GTP. The large subunit (R1) of BceSIV contains conserved motifs of NTPases and DNA helicases. The R1 subunit has no endonuclease activity by itself; it strongly stimulates REase activity when in complex with the R2 subunit. BceSIV was demonstrated to hydrolyze GTP and ATP in vitro. BceSIV is similar to CglI (GCSGC), and homologs of R1 are found in 11 sequenced bacterial genomes, where they are paired with specificity subunits. In addition, homologs of the BceSIV R1-R2 fusion are found in many sequenced microbial genomes. An orphan methylase, M.BceSV, was found to modify GCNGC, GGCC, CCGG, GGNNCC, and GCGC sites. A ParB-methylase fusion protein appears to nick DNA nonspecifically. The ATCC 14579 genome encodes an active enzyme Bce14579I (GCWGC). BceSIV and Bce14579I belong to the phospholipase D (PLD) family of endonucleases that are widely distributed among Bacteria and Archaea. A survey of type II and III restriction-modification (R-M) system genes is presented from sequenced B. cereus, Bacillus anthracis, and Bacillus thuringiensis strains.

Cloning and characterization of the TneDI restriction: modification system of Thermotoga neapolitana

A putative Type II restriction-modification system of Thermotoga neapolitana, TneDI, was cloned into Escherichia coli XL1-Blue MRF' and characterized. Gene CTN_0339 specifies the endonuclease R.TneDI, while CTN_0340 encodes the cognate DNA methyltransferase M.TneDI. Both enzymes were purified simply by heating the cell lysates of E. coli followed by centrifugation. The enzymes were active over a broad range of temperatures, from 42°C to at least 77°C, with the highest activities observed at 77°C. R.TneDI cleaved at the center of the recognition sequence (CG↓CG) and generated blunt-end cuts. Overexpression of R.TneDI in BL21(DE3) was confirmed by both SDS-PAGE and Western blotting.

The insertion sequences of Anabaena sp. strain PCC 7120 and their effects on its open reading frames

Anabaena sp. strain PCC 7120, widely studied, has 145 annotated transposase genes that are part of transposable elements called insertion sequences (ISs). To determine the entirety of the ISs, we aligned transposase genes and their flanking regions; identified the ISs' possible terminal inverted repeats, usually flanked by direct repeats; and compared IS-interrupted sequences with homologous sequences. We thereby determined both ends of 87 ISs bearing 110 transposase genes in eight IS families (http://www-is.biotoul.fr/) and in a cluster of unclassified ISs, and of hitherto unknown miniature inverted-repeat transposable elements. Open reading frames were then identified to which ISs contributed and others--some encoding proteins of predictable function, including protein kinases, and restriction endonucleases--that were interrupted by ISs. Anabaena sp. ISs were often more closely related to exogenous than to other endogenous ISs, suggesting that numerous variant ISs were not degraded within PCC 7120 but transferred from without. This observation leads to the expectation that further sequencing projects will extend this and similar analyses. We also propose an adaptive role for poly(A) sequences in ISs.

Predicted glycosyl transferase genes located outside the HEP island are required for formation of heterocyst envelope polysaccharide in Anabaena sp. strain PCC 7120

During maturation, heterocysts form an envelope layer of polysaccharide, called heterocyst envelope polysaccharide (HEP), whose synthesis depends on a cluster of genes, the HEP island, and on an additional, distant gene, hepB, or a gene immediately downstream from hepB. We show that HEP formation depends upon the predicted glycosyl transferase genes all4160 at a third locus and alr3699, which is adjacent to hepB and is cotranscribed with it. Mutations in the histidine kinase genes hepN and hepK appear to silence the promoter of hepB and incompletely down-regulate all4160.

Identification of a new subfamily of HNH nucleases and experimental characterization of a representative member, HphI restriction endonuclease

The restriction endonuclease (REase) R. HphI is a Type IIS enzyme that recognizes the asymmetric target DNA sequence 5'-GGTGA-3' and in the presence of Mg(2+) hydrolyzes phosphodiester bonds in both strands of the DNA at a distance of 8 nucleotides towards the 3' side of the target, producing a 1 nucleotide 3'-staggered cut in an unspecified sequence at this position. REases are typically ORFans that exhibit little similarity to each other and to any proteins in the database. However, bioinformatics analyses revealed that R.HphI is a member of a relatively big sequence family with a conserved C-terminal domain and a variable N-terminal domain. We predict that the C-terminal domains of proteins from this family correspond to the nuclease domain of the HNH superfamily rather than to the most common PD-(D/E)XK superfamily of nucleases. We constructed a three-dimensional model of the R.HphI catalytic domain and validated our predictions by site-directed mutagenesis and studies of DNA-binding and catalytic activities of the mutant proteins. We also analyzed the genomic neighborhood of R.HphI homologs and found that putative nucleases accompanied by a DNA methyltransferase (i.e. predicted REases) do not form a single group on a phylogenetic tree, but are dispersed among free-standing putative nucleases. This suggests that nucleases from the HNH superfamily were independently recruited to become REases in the context of RM systems multiple times in the evolution and that members of the HNH superfamily may be much more frequent among the so far unassigned REase sequences than previously thought.

Genome-wide analysis of restriction-modification system in unicellular and filamentous cyanobacteria

Cyanobacteria are an ancient group of gram-negative bacteria with strong genome size variation ranging from 1.6 to 9.1 Mb. Here, we first retrieved all the putative restriction-modification (RM) genes in the draft genome of Spirulina and then performed a range of comparative and bioinformatic analyses on RM genes from unicellular and filamentous cyanobacterial genomes. We have identified 6 gene clusters containing putative Type I RMs and 11 putative Type II RMs or the solitary methyltransferases (MTases). RT-PCR analysis reveals that 6 of 18 MTases are not expressed in Spirulina, whereas one hsdM gene, with a mutated cognate hsdS, was detected to be expressed. Our results indicate that the number of RM genes in filamentous cyanobacteria is significantly higher than in unicellular species, and this expansion of RM systems in filamentous cyanobacteria may be related to their wide range of ecological tolerance. Furthermore, a coevolutionary pattern is found between hsdM and hsdR, with a large number of site pairs positively or negatively correlated, indicating the functional importance of these pairing interactions between their tertiary structures. No evidence for positive selection is found for the majority of RMs, e.g., hsdM, hsdS, hsdR, and Type II restriction endonuclease gene families, while a group of MTases exhibit a remarkable signature of adaptive evolution. Sites and genes identified here to have been under positive selection would provide targets for further research on their structural and functional evaluations.

Identification of a new family of putative PD-(D/E)XK nucleases with unusual phylogenomic distribution and a new type of the active site

Background

Prediction of structure and function for uncharacterized protein families by identification of evolutionary links to characterized families and known structures is one of the cornerstones of genomics. Theoretical assignment of three-dimensional folds and prediction of protein function even at a very general level can facilitate the experimental determination of the molecular mechanism of action and the role that members of a given protein family fulfill in the cell. Here, we predict the three-dimensional fold and study the phylogenomic distribution of members of a large family of uncharacterized proteins classified in the Clusters of Orthologous Groups database as COG4636.

Results

Using protein fold-recognition we found that members of COG4636 are remotely related to Holliday junction resolvases and other nucleases from the PD-(D/E)XK superfamily. Structure modeling and sequence analyses suggest that most members of COG4636 exhibit a new, unusual variant of the putative active site, in which the catalytic Lys residue migrated in the sequence, but retained similar spatial position with respect to other functionally important residues. Sequence analyses revealed that members of COG4636 and their homologs are found mainly in Cyanobacteria, but also in other bacterial phyla. They undergo horizontal transfer and extensive proliferation in the colonized genomes; for instance in Gloeobacter violaceus PCC 7421 they comprise over 2% of all protein-encoding genes. Thus, members of COG4636 appear to be a new type of selfish genetic elements, which may fulfill an important role in the genome dynamics of Cyanobacteria and other species they invaded. Our analyses provide a platform for experimental determination of the molecular and cellular function of members of this large protein family.

Conclusion

After submission of this manuscript, a crystal structure of one of the COG4636 members was released in the Protein Data Bank (code 1wdj; Idaka, M., Wada, T., Murayama, K., Terada, T., Kuramitsu, S., Shirouzu, M., Yokoyama, S.: Crystal structure of Tt1808 from Thermus thermophilus Hb8, to be published). Our analysis of the Tt1808 structure reveals that we correctly predicted all functionally important features of the COG4636 family, including the membership in the PD-(D/E)xK superfamily of nucleases, the three-dimensional fold, the putative catalytic residues, and the unusual configuration of the active site.

DNA methyltransferases from Neisseria meningitidis and Neisseria gonorrhoeae FA1090 associated with mismatch nicking endonucleases

The genes encoding the DNA methyltransferases M.NmeDI and M.NmeAI from Neisseria meningitidis associated with the genes encoding putative Vsr endonucleases were overexpressed in Escherichia coli. The enzymes were purified to apparent homogeneity on Ni-NTA agarose columns, yielding proteins of 49+/-1 kDa and 39.6+/-1 kDa, respectively, under denaturing conditions. M.NmeDI recognizes the degenerate sequence 5'-RCCGGB-3'. It methylates the first 5' cytosine residue on both strands within the core sequence CCGG. The enzyme shows higher affinity with the hemimethylated degenerate sequence than with the unmethylated degenerate sequence. Comparison of the amino acid sequence of the target-recognizing domain of M.NmeDI with the closest neighbours recognizing the sequence 5'-RCCGGY-3' showed the presence of the homologous domain and an additional domain that may be responsible for recognizing the degenerate sequence. M.NmeAI recognizes the sequence 5'-CCGG-3' and methylates the second 5' cytosine residue on both DNA strands. In Neisseria gonorrhoeae strain FA1090 the homologues of these ORFs are truncated due to a variety of mutations.

Improved Genetic Transformation of the Thermophilic Cyanobacterium, Thermosynechococcus elongatus BP-1

We improved genetic transformation of the thermophilic cyanobacterium, Thermosynechococcus elongatus BP-1, by combining electroporation with a top agar method. Transformation was also improved when a disruptant of a putative type I restriction endonuclease (tll2230) was used as recipient cells. In particular, some constructs, with which wild type has never been transformed, were successfully integrated into the tll2230-disruptant. Single-crossover recombination was detected more frequently than the double-crossover recombination. In accordance with the presence of all the homologs of pil genes in Synechocystis sp. PCC 6803, we found that T. elongatus is naturally transformable with exogenous DNA.

An overview of the genome of Nostoc punctiforme, a multicellular, symbiotic cyanobacterium

Nostoc punctiforme is a filamentous cyanobacterium with extensive phenotypic characteristics and a relatively large genome, approaching 10 Mb. The phenotypic characteristics include a photoautotrophic, diazotrophic mode of growth, but N. punctiforme is also facultatively heterotrophic; its vegetative cells have multiple developmental alternatives, including terminal differentiation into nitrogen-fixing heterocysts and transient differentiation into spore-like akinetes or motile filaments called hormogonia; and N. punctiforme has broad symbiotic competence with fungi and terrestrial plants, including bryophytes, gymnosperms and an angiosperm. The shotgun-sequencing phase of the N. punctiforme strain ATCC 29133 genome has been completed by the Joint Genome Institute. Annotation of an 8.9 Mb database yielded 7432 open reading frames, 45% of which encode proteins with known or probable known function and 29% of which are unique to N. punctiforme. Comparative analysis of the sequence indicates a genome that is highly plastic and in a state of flux, with numerous insertion sequences and multilocus repeats, as well as genes encoding transposases and DNA modification enzymes. The sequence also reveals the presence of genes encoding putative proteins that collectively define almost all characteristics of cyanobacteria as a group. N. punctiforme has an extensive potential to sense and respond to environmental signals as reflected by the presence of more than 400 genes encoding sensor protein kinases, response regulators and other transcriptional factors. The signal transduction systems and any of the large number of unique genes may play essential roles in the cell differentiation and symbiotic interaction properties of N. punctiforme.