LESSONS FROM SEQUENCING OF THE GENOME OF A UNICELLULAR CYANOBACTERIUM, SYNECHOCYSTIS SP. PCC6803

Analysis of correlation structures in the Synechocystis PCC6803 genome

Transfer of nucleotide strings in the Synechocystis sp. PCC6803 genome is investigated to exhibit periodic and non-periodic correlation structures by using the recurrence plot method and the phase space reconstruction technique. The periodic correlation structures are generated by periodic transfer of several substrings in long periodic or non-periodic nucleotide strings embedded in the coding regions of genes. The non-periodic correlation structures are generated by non-periodic transfer of several substrings covering or overlapping with the coding regions of genes. In the periodic and non-periodic transfer, some gaps divide the long nucleotide strings into the substrings and prevent their global transfer. Most of the gaps are either the replacement of one base or the insertion/reduction of one base. In the reconstructed phase space, the points generated from two or three steps for the continuous iterative transfer via the second maximal distance can be fitted by two lines. It partly reveals an intrinsic dynamics in the transfer of nucleotide strings. Due to the comparison of the relative positions and lengths, the substrings concerned with the non-periodic correlation structures are almost identical to the mobile elements annotated in the genome. The mobile elements are thus endowed with the basic results on the correlation structures.

The gene ssl3076 encodes a protein mediating the salt-induced expression of ggpS for the biosynthesis of the compatible solute glucosylglycerol in Synechocystis sp. strain PCC 6803

Acclimation to high salt concentrations involves concerted changes in gene expression. For the majority of salt-regulated genes, the mechanism underlying the induction process is not known. The gene ggpS (sll1566), which encodes the glucosylglycerol-phosphate synthase responsible for the synthesis of the compatible solute glucosylglycerol (GG), is specifically induced by salt in the cyanobacterial model strain Synechocystis sp. strain PCC 6803. To identify mechanisms mediating this salt-specific gene regulation, the ggpS promoter was analyzed in more detail. 5' rapid amplification of cDNA ends (5'-RACE) experiments revealed that the adjacent open reading frame (ORF), which is annotated as unknown protein Ssl3076, overlaps with the transcriptional start site of the ggpS gene. Reporter gene expression analyses indicated an essential role for the intact ssl3076 gene in the salt-regulated transcription of a gfp reporter gene. Promoter fragments containing a mutated ssl3076 lost the salt regulation; similarly, a frameshift mutation in ssl3076 resulted in a high level of ggpS expression under low-salt conditions, thereby establishing this small ORF, named ggpR, as a negative regulator of ggpS. Interestingly, small ORFs were also found adjacent to ggpS genes in the genomes of other GG-accumulating cyanobacteria. These results suggest that the GgpR protein represses ggpS expression under low-salt conditions, whereas in salt-shocked and salt-acclimated cells a stress-proportional ggpS expression occurs, leading to GG accumulation.

Structural characterization of NDH-1 complexes of Thermosynechococcus elongatus by single particle electron microscopy

The structure of the multifunctional NAD(P)H dehydrogenase type 1 (NDH-1) complexes from cyanobacteria was investigated by growing the wild type and specific ndh His-tag mutants of Thermosynechococcus elongatus BP-1 under different CO(2) conditions, followed by an electron microscopy (EM) analysis of their purified membrane protein complexes. Single particle averaging showed that the complete NDH-1 complex (NDH-1L) is L-shaped, with a relatively short hydrophilic arm. Two smaller complexes were observed, differing only at the tip of the membrane-embedded arm. The smallest one is considered to be similar to NDH-1M, lacking the NdhD1 and NdhF1 subunits. The other fragment, named NDH-1I, is intermediate between NDH-1L and NDH-1M and only lacks a mass compatible with the size of the NdhF1 subunit. Both smaller complexes were observed under low- and high-CO(2) growth conditions, but were much more abundant under the latter conditions. EM characterization of cyanobacterial NDH-1 further showed small numbers of NDH-1 complexes with additional masses. One type of particle has a much longer peripheral arm, similar to the one of NADH: ubiquinone oxidoreductase (complex I) in E. coli and other organisms. This indicates that Thermosynechococcus elongatus must have protein(s) which are structurally homologous to the E. coli NuoE, -F, and -G subunits. Another low-abundance type of particle (NDH-1U) has a second labile hydrophilic arm at the tip of the membrane-embedded arm. This U-shaped particle has not been observed before by EM in a NDH-I preparation.

Genome comparison using Gene Ontology (GO) with statistical testing

BACKGROUND

Automated comparison of complete sets of genes encoded in two genomes can provide insight on the genetic basis of differences in biological traits between species. Gene ontology (GO) is used as a common vocabulary to annotate genes for comparison. Current approaches calculate the fold of unweighted or weighted differences between two species at the high-level GO functional categories. However, to ensure the reliability of the differences detected, it is important to evaluate their statistical significance. It is also useful to search for differences at all levels of GO.

RESULTS

We propose a statistical approach to find reliable differences between the complete sets of genes encoded in two genomes at all levels of GO. The genes are first assigned GO terms from BLAST searches against genes with known GO assignments, and for each GO term the abundance of genes in the two genomes is compared using a chi-squared test followed by false discovery rate (FDR) correction. We applied this method to find statistically significant differences between two cyanobacteria, Synechocystis sp. PCC6803 and Anabaena sp. PCC7120. We then studied how the set of identified differences vary when different BLAST cutoffs are used. We also studied how the results vary when only subsets of the genes were used in the comparison of human vs. mouse and that of Saccharomyces cerevisiae vs. Schizosaccharomyces pombe.

CONCLUSION

There is a surprising lack of statistical approaches for comparing complete genomes at all levels of GO. With the rapid increase of the number of sequenced genomes, we hope that the approach we proposed and tested can make valuable contribution to comparative genomics.

Cyanobacterial two-component proteins: structure, diversity, distribution, and evolution

A survey of the already characterized and potential two-component protein sequences that exist in the nine complete and seven partially annotated cyanobacterial genome sequences available (as of May 2005) showed that the cyanobacteria possess a much larger repertoire of such proteins than most other bacteria. By analysis of the domain structure of the 1,171 potential histidine kinases, response regulators, and hybrid kinases, many various arrangements of about thirty different modules could be distinguished. The number of two-component proteins is related in part to genome size but also to the variety of physiological properties and ecophysiologies of the different strains. Groups of orthologues were defined, only a few of which have representatives with known physiological functions. Based on comparisons with the proposed phylogenetic relationships between the strains, the orthology groups show that (i) a few genes, some of them clustered on the genome, have been conserved by all species, suggesting their very ancient origin and an essential role for the corresponding proteins, and (ii) duplications, fusions, gene losses, insertions, and deletions, as well as domain shuffling, occurred during evolution, leading to the extant repertoire. These mechanisms are put in perspective with the different genetic properties that cyanobacteria have to achieve genome plasticity. This review is designed to serve as a basis for orienting further research aimed at defining the most ancient regulatory mechanisms and understanding how evolution worked to select and keep the most appropriate systems for cyanobacteria to develop in the quite different environments that they have successfully colonized.

"Evolution Canyon": A Microcosm of Life's Evolution Focusing on Adaptation and Speciation

Local microcosmic natural laboratories, dubbed "Evolution Canyon" (EC) models, reinforce studies of regional and global macrocosmic ecological theaters across life and unravelevolution in action.The EC model laboratories permit genomic, proteomic, and phenomic studies highlighting speciation and adaptation at a microscale. Critical transplant experiment tests can evaluate interslope differential fitness. Novel techniques of genetic mapping, sequence nucleotide polymorphisms (SNPs), and wide genome coding and noncoding expressions can unravel evolutionary dynamics. Finally, fundamental problems such as stress effects on nonrandom mutations, lateral transfers, splicing variations, sex, and social evolutions, and adaptive strategies of prokaryotes and eukaryotes are testable. We are studying four "Evolution Canyons" (EC I-IV) in the Carmel, Galilee, Negev, and Golan mountains. We've identified 2,500 species in EC I (Carmel) from bacteria to mammals in an area of 7,000 m. Higher terrestrial species richness was found on the more stressful tropical "African" slope (AS). Aquatic species richness was higher on the milder, temperate "European" slope (ES). In soil fungi we found interslope local and regional adaptive divergence in sex, melanism, and conidia. In nine out of 14 (64%) model organisms across life we identified largely higher genetic polymorphisms on the more stressful "African" slope. In some model species we found higher mutation rates, gene conversion, recombination, DNA repair, and larger genome size as well as interslope divergent micro-satellites, molecular polymorphisms, retrotransposons, and wide-genome gene expression on the more stressful AS. Remarkably, interslope incipient sympatric speciation was found across life. The "Evolution Canyon" model represents the Israeli ecological equivalent of the Galapagos Islands. Micro-climatic selection overrides drift and drives both interslope adaptive radiation and incipient sympatric speciation. The EC model could potentially highlight many mysteries of evolutionary biology.

Adaptation of the photosynthetic electron transport chain in cyanobacteria to iron deficiency: The function of IdiA and IsiA

In this review we give an overview on the adaptational responses of photosystem (PS) II and PSI in cyanobacteria to iron starvation, mainly summarizing our results with the mesophilic Synechococcus elongatus PCC 7942. We also discuss this process with respect to the strong interrelationship between iron limitation and oxidative stress that exists in cyanobacteria as oxygenic photosynthetic organisms. The adaptation of the multiprotein complexes PSII and PSI to iron starvation is a sequential process, which is characterized by the enhanced expression of two major iron-regulated proteins, IdiA (iron deficiency induced protein A) and IsiA (iron stress induced protein A). Our results suggest that IdiA protects the acceptor side of PSII against oxidative stress under conditions of mild iron limitation in a currently unclear way, whereas prolonged iron deficiency leads to the synthesis of a chlorophyll a antenna around PSI-trimers consisting of IsiA molecules. The physiological consequences of these alterations under prolonged iron starvation, as shown by acridine yellow fluorescence measurements, are a reduction of linear electron transport activity through PSII and an increase of cyclic electron flow around PSI as well as an increase in respiratory activity. IdiA and IsiA expression are mediated by two distinct helix-turn-helix transcriptional regulators of the Crp/Fnr family. IdiB positively regulates expression of idiA under iron starvation, and Fur represses transcription of isiA under iron-sufficient conditions. Although both transcriptional regulators seem to operate independently of each other, our results indicate that a cross-talk between the signal transduction pathways exists. Moreover, IdiA as well as IsiA expression are affected by hydrogen peroxide. We suggest that due to the interdependence of iron limitation and the formation of reactive oxygen species, peroxide stress might be the superior trigger that leads to expression of these proteins under iron starvation. The modifications of PSII and PSI under iron starvation influence the redox state of redox-sensitive components of the electron transport chain, and thus the activity of metabolic pathways being regulated in dependence of the redox state of these components.

Plastid protein synthesis is required for plant development in tobacco

Chloroplasts fulfill important functions in cellular metabolism. The majority of plastid genome-encoded genes is involved in either photosynthesis or chloroplast gene expression. Whether or not plastid genes also can determine extraplastidic functions has remained controversial. We demonstrate here an essential role of plastid protein synthesis in tobacco leaf development. By using chloroplast transformation, we have developed an experimental system that produces recombination-based knockouts of chloroplast translation in a cell-line-specific manner. The resulting plants are chimeric and, in the presence of translational inhibitors, exhibit severe developmental abnormalities. In the absence of active plastid protein synthesis, leaf blade development is abolished because of an apparent arrest of cell division. This effect appears to be cell-autonomous in that adjacent sectors of cells with translating plastids are phenotypically normal but cannot complement for the absence of plastid translation in mutant sectors. Developmental abnormalities also are seen in flower morphology, indicating that the defects are not caused by inhibited expression of plastid photosynthesis genes. Taken together, our data point to an unexpected essential role of plastid genes and gene expression in plant development and cell division.

Hydrogenases and hydrogen metabolism of cyanobacteria

Cyanobacteria may possess several enzymes that are directly involved in dihydrogen metabolism: nitrogenase(s) catalyzing the production of hydrogen concomitantly with the reduction of dinitrogen to ammonia, an uptake hydrogenase (encoded by hupSL) catalyzing the consumption of hydrogen produced by the nitrogenase, and a bidirectional hydrogenase (encoded by hoxFUYH) which has the capacity to both take up and produce hydrogen. This review summarizes our knowledge about cyanobacterial hydrogenases, focusing on recent progress since the first molecular information was published in 1995. It presents the molecular knowledge about cyanobacterial hupSL and hoxFUYH, their corresponding gene products, and their accessory genes before finishing with an applied aspect--the use of cyanobacteria in a biological, renewable production of the future energy carrier molecular hydrogen. In addition to scientific publications, information from three cyanobacterial genomes, the unicellular Synechocystis strain PCC 6803 and the filamentous heterocystous Anabaena strain PCC 7120 and Nostoc punctiforme (PCC 73102/ATCC 29133) is included.

How do cyanobacteria sense and respond to light?

Cyanobacteria exhibit numerous responses to changes in the intensity and spectral quality of light. What sensors do cyanobacteria use to detect light and what are the mechanisms of signal transduction? The publication in 1996 of the complete genome sequence of the cyanobacterium Synechocystis 6803 provided a tremendous stimulus for research in this field, and many light-sensors and signal transducers have now been identified. However, our knowledge of cyanobacterial light-signal transduction remains fragmentary. This review summarizes what we know about the ways in which cyanobacteria perceive light, some of the ways which they respond to light signals and some recent achievements in elucidating the signal transduction mechanisms. Some problems in characterizing cyanobacterial signal transduction pathways are outlined and alternative experimental strategies are discussed.

The complete DNA sequence of the Ectocarpus siliculosus Virus EsV-1 genome

The Ectocarpus siliculosus Virus-1, EsV-1, is the type-species of a genus of Phycodnaviridae, the phaeoviruses, infecting marine filamentous brown algae. The EsV-1 genome of 335,593 bp contains tandem and dispersed repetitive elements in addition to a large number of open reading frames of which 231 are currently counted as genes. Many genes can be assigned to functional groups involved in DNA synthesis, DNA integration, transposition, and polysaccharide metabolism. Furthermore, EsV-1 contains components of a surprisingly complex signal transduction system with six different hybrid histidine protein kinases and four putative serine/threonine protein kinases. Several other genes encode polypeptides with protein-protein interaction domains. However, 50% of the predicted genes have no counterparts in data banks. Only 28 of the 231 identified genes have significant sequence similarities to genes of the Chlorella virus PBCV-1, another phycodnavirus. To our knowledge, the EsV-1 genome is the largest viral DNA sequenced to date.

Evolution of genome-phenome diversity under environmental stress

The genomic era revolutionized evolutionary biology. The enigma of genotypic-phenotypic diversity and biodiversity evolution of genes, genomes, phenomes, and biomes, reviewed here, was central in the research program of the Institute of Evolution, University of Haifa, since 1975. We explored the following questions. (i) How much of the genomic and phenomic diversity in nature is adaptive and processed by natural selection? (ii) What is the origin and evolution of adaptation and speciation processes under spatiotemporal variables and stressful macrogeographic and microgeographic environments? We advanced ecological genetics into ecological genomics and analyzed globally ecological, demographic, and life history variables in 1,200 diverse species across life, thousands of populations, and tens of thousands of individuals tested mostly for allozyme and partly for DNA diversity. Likewise, we tested thermal, chemical, climatic, and biotic stresses in several model organisms. Recently, we introduced genetic maps and quantitative trait loci to elucidate the genetic basis of adaptation and speciation. The genome-phenome holistic model was deciphered by the global regressive, progressive, and convergent evolution of subterranean mammals. Our results indicate abundant genotypic and phenotypic diversity in nature. The organization and evolution of molecular and organismal diversity in nature at global, regional, and local scales are nonrandom and structured; display regularities across life; and are positively correlated with, and partly predictable by, abiotic and biotic environmental heterogeneity and stress. Biodiversity evolution, even in small isolated populations, is primarily driven by natural selection, including diversifying, balancing, cyclical, and purifying selective regimes, interacting with, but ultimately overriding, the effects of mutation, migration, and stochasticity.

Sequence and analysis of the Arabidopsis genome

The comprehensive analysis of the genome sequence of the plant Arabidopsis thaliana has been completed recently. The genome sequence and associated analyses provide the foundations for rapid progress in many fields of plant research, such as the exploitation of genetic variation in Arabidopsis ecotypes, the assessment of the transcriptome and proteome, and the association of genome changes at the sequence level with evolutionary processes. Nevertheless, genome sequencing and analysis are only the first steps towards a new plant biology. Much remains to be done to refine the analysis of encoded genes, to define the functions of encoded proteins systematically, and to establish new generations of databases to capture and relate diverse data sets generated in widely distributed laboratories.

Analysis of the genome sequence of the flowering plant Arabidopsis thaliana

The flowering plant Arabidopsis thaliana is an important model system for identifying genes and determining their functions. Here we report the analysis of the genomic sequence of Arabidopsis. The sequenced regions cover 115.4 megabases of the 125-megabase genome and extend into centromeric regions. The evolution of Arabidopsis involved a whole-genome duplication, followed by subsequent gene loss and extensive local gene duplications, giving rise to a dynamic genome enriched by lateral gene transfer from a cyanobacterial-like ancestor of the plastid. The genome contains 25,498 genes encoding proteins from 11,000 families, similar to the functional diversity of Drosophila and Caenorhabditis elegans--the other sequenced multicellular eukaryotes. Arabidopsis has many families of new proteins but also lacks several common protein families, indicating that the sets of common proteins have undergone differential expansion and contraction in the three multicellular eukaryotes. This is the first complete genome sequence of a plant and provides the foundations for more comprehensive comparison of conserved processes in all eukaryotes, identifying a wide range of plant-specific gene functions and establishing rapid systematic ways to identify genes for crop improvement.

Genomic survey of cAMP and cGMP signalling components in the cyanobacterium Synechocystis PCC 6803

Cyanobacteria modulate intracellular levels of cAMP and cGMP in response to environmental conditions (light, nutrients and pH). In an attempt to identify components of the cAMP and cGMP signalling pathways in Synechocystis PCC 6803, the authors screened its complete genome sequence by using bioinformatic tools and data from sequence-function studies performed on both eukaryotic and prokaryotic cAMP/cGMP-dependent proteins. Sll1624 and Slr2100 were tentatively assigned as being two putative cyclic nucleotide phosphodiesterases. Five proteins were identified as having all the determinants required to be cyclic nucleotide receptors, two of them being probably more specific for cGMP (an element of two-component regulatory systems - Slr2104 - and a putative cyclic-nucleotide-gated cation channel - Slr1575), the three others being probably more specific for cAMP: (i) a protein of unidentified function (Slr0842); (ii) a putative cyclic-nucleotide-modulated permease (Slr0593), previously annotated as a kinase A regulatory subunit; and (iii) a putative transcription factor (CRP-SYN: =Sll1371), which possesses cAMP- and DNA-binding determinants homologous to those of the cAMP receptor protein of Escherichia coli (CRP-EC:). This homology, together with the presence in Synechocystis of CRP-EC:-like binding sites upstream of crp, cya1, slr1575, and several genes encoding enzymes involved in transport and metabolism, strongly suggests that CRP-SYN: is a global regulator.

Isolation of regulated genes of the cyanobacterium Synechocystis sp. strain PCC 6803 by differential display

Global identification of differentially regulated genes in prokaryotes is constrained because the mRNA does not have a 3' polyadenylation extension; this precludes specific separation of mRNA from rRNA and tRNA and synthesis of cDNAs from the entire mRNA population. Knowledge of the entire genome sequence of Synechocystis sp. strain PCC 6803 has enabled us to develop a differential display procedure that takes advantage of a short palindromic sequence that is dispersed throughout the Synechocystis sp. strain PCC 6803 genome. This sequence, designated the HIP (highly iterated palindrome) element, occurs in approximately half of the Synechocystis sp. strain PCC 6803 genes but is absent in rRNA and tRNA genes. To determine the feasibility of exploiting the HIP element, alone or in combination with specific primer subsets, for analyzing differential gene expression, we used HIP-based primers to identify light intensity-regulated genes. Several gene fragments, including those encoding ribosomal proteins and phycobiliprotein subunits, were differentially amplified from RNA templates derived from cells grown in low light or exposed to high light for 3 h. One novel finding was that expression of certain genes of the pho regulon, which are under the control of environmental phosphate levels, were markedly elevated in high light. High-light activation of pho regulon genes correlated with elevated growth rates that occur when the cells are transferred from low to high light. These results suggest that in high light, the rate of growth of Synechocystis sp. strain PCC 6803 exceeds its capacity to assimilate phosphate, which, in turn, may trigger a phosphate starvation response and activation of the pho regulon.

Multiplication of antenna genes as a major adaptation to low light in a marine prokaryote

Two ecotypes of the prokaryote Prochlorococcus adapted to distinct light niches in the ocean have been described recently. These ecotypes are characterized by their different (divinyl-) chlorophyll (Chl) a to Chl b ratios and 16S rRNA gene signatures, as well as by their significantly distinct irradiance optima for growth and photosynthesis [Moore, L. R., Rocap, G. & Chisholm, S. W. (1998) Nature (London) 393, 464-467]. However, the molecular basis of their physiological differences remained, so far, unexplained. In this paper, we show that the low-light-adapted Prochlorococcus strain SS120 possesses a gene family of seven transcribed genes encoding different Chl a/b-binding proteins (Pcbs). In contrast, Prochlorococcus sp. MED4, a high-light-adapted ecotype, possesses a single pcb gene. The presence of multiple antenna genes in another low-light ecotype (NATL2a), but not in another high-light ecotype (TAK9803-2), is demonstrated. Thus, the multiplication of pcb genes appears as a key factor in the capacity of deep Prochlorococcus populations to survive at extremely low photon fluxes.

Physiologic processes affecting the content and distribution of phytonutrients in plants

Plant physiologists and biochemists are unraveling the transport mechanisms and biosynthetic pathways that determine each plant's unique phytonutrient composition so that crop plants can be modified to improve their nutrition quality. However, before this goal can be achieved, more information is needed on various plant and human phytonutrient-related processes, and some new technical capabilities are required. The current status of our knowledge base and recommendations for technical improvements and research priorities will be discussed.

A gene required for the regulation of photosynthetic light harvesting in the cyanobacterium Synechocystis 6803

A gene required for the short-term regulation of photosynthetic light harvesting (the state transition) has been identified in the cyanobacterium Synechocystis sp. PCC6803. The open reading frame is designated sll1926 in the complete Synechocystis gene sequence. The deduced amino acid sequence has no homologues in current sequence databases and no recognizable sequence motifs. It encodes a putative integral membrane protein of 16 kDa, which we have designated RpaC (regulator of phycobilisome association C). Fluorescence measurements of an insertional inactivation mutant of rpaC (Deltasll1926) show that it is specifically unable to perform state transitions. Deltasll1926 has approximately wild-type levels of PS1, PS2 and phycobilisomes. Measurements of oxygen evolution and uptake show Deltasll1926 to have no deficiency in electron transport rates. In vitro [gamma-32P]-ATP labelling experiments suggest that RpaC is not the 15 kDa membrane phosphoprotein previously implicated in state transitions. Deltasll1926 grows more slowly than the wild type only at very low light intensities.