Chromosome-scale assembly of the streamlined picoeukaryote Picochlorum sp. SENEW3 genome reveals Rabl-like chromatin structure and potential for C4 photosynthesis

Genome sequencing and assembly of the photosynthetic picoeukaryotic Picochlorum sp. SENEW3 revealed a compact genome with a reduced gene set, few repetitive sequences, and an organized Rabl-like chromatin structure. Hi-C chromosome conformation capture revealed evidence of possible chromosomal translocations, as well as putative centromere locations. Maintenance of a relatively few selenoproteins, as compared to similarly sized marine picoprasinophytes Mamiellales, and broad halotolerance compared to others in Trebouxiophyceae, suggests evolutionary adaptation to variable salinity environments. Such adaptation may have driven size and genome minimization and have been enabled by the retention of a high number of membrane transporters. Identification of required pathway genes for both CAM and C4 photosynthetic carbon fixation, known to exist in the marine mamiellale pico-prasinophytes and seaweed Ulva, but few other chlorophyte species, further highlights the unique adaptations of this robust alga. This high-quality assembly provides a significant advance in the resources available for genomic investigations of this and other photosynthetic picoeukaryotes.

Diverse Microbial Hot Spring Mat Communities at Black Canyon of the Colorado River

The thermophilic microbial mat communities at hot springs in the Black Canyon of the Colorado River, thought to harbor the protistan human pathogen Naegleria fowleri, were surveyed using both culture-independent and -dependent methods to further understand the ecology of these hot spring microbiomes. Originating from Lake Mead source water, seven spring sites were sampled, varying in temperature from 25 to 55 °C. Amplicon-based high-throughput sequencing of twelve samples using 16S rRNA primers (hypervariable V4 region) revealed that most mats are dominated by cyanobacterial taxa, some but not all similar to those dominating the mats at other studied hot spring systems. 18S rRNA amplicon sequencing (V9 region) demonstrated a diverse community of protists and other eukaryotes including a highly abundant amoebal sequence related to Echinamoeba thermarum. Additional taxonomic and diversity metric analyses using near full-length 16S and 18S rRNA gene sequencing allowed a higher sequence-based resolution of the community. The mat sequence data suggest a major diversification of the cyanobacterial orders Leptolyngbyales, as well as microdiversity among several cyanobacterial taxa. Cyanobacterial isolates included some representatives of ecologically abundant taxa. A Spearman correlation analysis of short-read amplicon sequencing data supported the co-occurrences of populations of cyanobacteria, chloroflexi, and bacteroidetes providing evidence of common microbial co-occurrences across the Black Canyon hot springs.

Diversity and putative interactions of parasitic alveolates belonging to Syndiniales at a coastal Pacific site

The dinoflagellate lineage Syndiniales currently consists entirely of parasitic species that fall into five well-supported clades. Environmental sequencing studies worldwide have found an abundance of Syndiniales in a variety of marine ecosystems, but very little is known about the majority of Syndiniales species including two entire clades which have only been observed in sequence data. Syndiniales are known to have a wide range of hosts, but only a few dozen interactions have been confirmed through observation of actual infections. Here, we describe the diversity of Syndiniales found at the Scripps Institution of Oceanography pier over the course of a year based on 18S sequencing. We find Syndiniales to be the most species (amplicon sequence variant)-rich taxonomic group and for its members to be present and abundant throughout the year. We used several analytical techniques to identify potential parasite-host interactions which we were then able to visualize over time. Using mock communities and size fractionation of seawater, we suggest that the majority of Syndiniales sequences that are found in environmental studies belong to the free-living dinospore stage rather than representing active infections.

Relating sinking and suspended microbial communities in the California Current Ecosystem: digestion resistance and the contributions of phytoplankton taxa to export

We used 16S, 18S, plastid and internal transcribed spacer (for Synechococcus strains) sequencing to quantify relative microbial abundances in water-column samples and on sediment-trap-collected particles across an environmental gradient in the California Current Ecosystem (CCE) spanning a > 60-fold range of surface chlorophyll. Most mixed-layer dominant eukaryotes and prokaryotes were consistently underrepresented on sinking particles. Diatoms were the only phototrophic taxa consistently overrepresented. Even within this class, however, one genus (Thalassiosira) was a particle-enriched dominant, while a similarly abundant species was poorly represented. Synechococcus was significantly enriched on sinking particles at only one of four sites, but clade I was disproportionately abundant on sinking particles throughout the region compared with clade IV, the euphotic-zone co-dominant. The most abundant microbes on particles across the CCE were organisms with distributional maxima close to the sediment-trap depth (rhizarians), microbes associated with metazoans or sinking particles as a nutritional habitat (certain alveolates, Gammaproteobacteria) and organisms that resist digestive degradation of their DNA (Thalassiosira, Synechococcus). For assessing taxon contributions of phytoplankton to carbon export, our results highlight the need for sequence-based quantitative approaches that can be used to integrate euphotic-zone abundances, compute rates and account for taxon differences in preservation of sequence markers through trophic processing.

Screening and characterization of polyhydroxyalkanoate granules, and phylogenetic analysis of polyhydroxyalkanoate synthase gene PhaC in cyanobacteria

Using Nile Red and BODIPY 493/503 dye-staining and fluorescence microscopy, twenty cyanobacterial strains, including ten commercially available strains and ten environmental isolates from estuaries, freshwater ponds, and lagoons, were screened for the accumulation of ecologically important and potentially biotechnologically significant carbon storage granules such as polyhydroxyalkanoates (PHA). Dye-staining granules were observed in six strains. Three Synechocystis, spp. strains WHSYN, LSNM, and CGF-1, and a Phormidium-like sp. CGFILA were isolated from environmental sources and found to produce granules of polyhydroxyalkanoate (PHA) according to PHA synthase gene (phaC) PCR screening and ¹ H NMR analyses. The environmental isolate, Nodularia sp. Las Olas and commercially available Phormidium cf. iriguum CCALA 759 displayed granules but screened negative for PHA according to phaC PCR and ¹ H NMR analyses. Partial polyhydroxyalkanoate synthase subunit C (phaC) and 16S rRNA gene sequences obtained from the PHA-accumulating strains and analyzed alongside publicly available phaC, phaE, 16S rRNA, and 23S rRNA data help in understanding the distribution and evolutionary history of PHA biosynthesis within the phylum Cyanobacteria. The data show that the presence of phaC is highly conserved within the genus Synechocystis, and present in at least one isolate of Phormidium. Maximum likelihood analyses and cophylogenetic modeling of PHA synthase gene sequences provide evidence of a recent horizontal gene transfer event between distant genera of cyanobacteria related to Pleurocapsa sp. PCC 7327 and Phormidium-like sp. CGFILA. These findings will help guide additional screening for PHA producers, and may explain why some Phormidium species produce PHAs, while others do not.

Spatial and temporal variations in Synechococcus microdiversity in the Southern California coastal ecosystem

The Synechococcus cyanobacterial population at the Scripps Institution of Oceanography pier in La Jolla, CA, shows large increases in abundance, typically in the spring and summer followed, by rapid declines within weeks. Here we used amplicon sequencing of the ribosomal RNA internal transcribed spacer region to examine the microdiversity within this cyanobacterial genus during these blooms as well as further offshore in the Southern California coastal ecosystem (CCE). These analyses revealed numerous Synechococcus amplicon sequence variants (ASVs) and that clade and ASV composition can change over the course of blooms. We also found that a large bloom in August 2016 was highly anomalous both in its overall Synechococcus abundance and in terms of the presence of normally oligotrophic Synechococcus clade II. The dominant ASVs at the pier were found further offshore and in the California Current, but we did observe more oligotrophic ASVs and clades along with depth variation in Synechococcus diversity. We also observed that the dominant sequence variant switched during the peak of multiple Synechococcus blooms, with this switch occurring in multiple clades, but we present initial evidence that this apparent ASV switch is a physiological response rather than a change in the dominant population.

Temporal dynamics of eukaryotic microbial diversity at a coastal Pacific site

High-throughput sequencing of ocean biomes has revealed vast eukaryotic microbial diversity, a significant proportion of which remains uncharacterized. Here we use a temporal approach to understanding eukaryotic diversity at the Scripps Pier, La Jolla, California, USA, via high-throughput amplicon sequencing of the 18S rRNA gene, the abundances of both Synechococcus and Synechococcus grazers, and traditional oceanographic parameters. We also exploit our ability to track operational taxonomic units (OTUs) temporally to evaluate the ability of 18S sequence-based OTU assignments to meaningfully reflect ecological dynamics. The eukaryotic community is highly dynamic in terms of both species richness and composition, although proportional representation of higher-order taxa remains fairly consistent over time. Synechococcus abundance fluctuates throughout the year. OTUs unique to dates of Synechococcus blooms and crashes or enriched in Synechococcus addition incubation experiments suggest that the prasinophyte Tetraselmis sp. and Gymnodinium-like dinoflagellates are likely Synechococcus grazers under certain conditions, and may play an important role in their population fluctuations.

Feeding and grazing impact by the bloom-forming euglenophyte Eutreptiella eupharyngea on marine eubacteria and cyanobacteria

The phototrophic euglenophyte Eutreptiella eupharyngea often causes blooms in the coastal waters of many countries, but its mode of nutrition has not been assessed. This species has previously been considered as exclusively auxotrophic. To explore whether E. eupharyngea is a mixotrophic species, the protoplasm of E. eupharyngea cells were examined using light, epifluorescence, and transmission electron microscopy after eubacteria, the cyanobacterium Synechococcus sp., and diverse algal species were provided as potential prey. Furthermore, the ingestion rates of E. eupharyngea KR on eubacteria or Synechococcus sp. as a function of prey concentration were measured. In addition, grazing by natural populations of euglenophytes on natural populations of eubacteria in Masan Bay was investigated. This study is the first to report that E. eupharyngea is a mixotrophic species. Among the potential prey organisms offered, E. eupharyngea fed only on eubacteria and Synechococcus sp., and the maximum ingestion rates of these two organisms measured in the laboratory were 5.7 and 0.7 cells predator^-1 h^-1, respectively. During the field experiments, the maximum ingestion rates and grazing impacts of euglenophytes, including E. eupharyngea, on natural populations of eubacteria were 11.8 cells predator^-1 h^-1 and 1.228 d^-1, respectively. Therefore, euglenophytes could potentially have a considerable grazing impact on marine bacterial populations.

Use of plankton-derived vitamin B1 precursors, especially thiazole-related precursor, by key marine picoeukaryotic phytoplankton

Several cosmopolitan marine picoeukaryotic phytoplankton are B1 auxotrophs requiring exogenous vitamin B1 or precursor to survive. From genomic evidence, representatives of picoeukaryotic phytoplankton (Ostreococcus and Micromonas spp.) were predicted to use known thiazole and pyrimidine B1 precursors to meet their B1 demands, however, recent culture-based experiments could not confirm this assumption. We hypothesized these phytoplankton strains could grow on precursors alone, but required a thiazole-related precursor other the well-known and extensively tested 4-methyl-5-thiazoleethanol. This hypothesis was tested using bioassays and co-cultures of picoeukaryotic phytoplankton and bacteria. We found that specific B1-synthesizing proteobacteria and phytoplankton are sources of a yet-to-be chemically identified thiazole-related precursor(s) that, along with pyrimidine B1 precursor 4-amino-5-hydroxymethyl-2-methylpyrimidine, can support growth of Ostreococcus spp. (also Micromonas spp.) without B1. We additionally found that the B1-synthesizing plankton do not require contact with picoeukaryotic phytoplankton cells to produce thiazole-related precursor(s). Experiments with wild-type and genetically engineered Ostreococcus lines revealed that the thiazole kinase, ThiM, is required for growth on precursors, and that thiazole-related precursor(s) accumulate to appreciable levels in the euphotic ocean. Overall, our results point to thiazole-related B1 precursors as important micronutrients promoting the survival of abundant phytoplankton influencing surface ocean production and biogeochemical cycling.

Diversity and genome dynamics of marine cyanophages using metagenomic analyses

Cyanophages are abundant in the oceanic environment and directly impact cyanobacterial distributions, physiological processes and evolution. Two samples collected from coastal Maine in July and September 2009 were enriched for Synechococcus cells using flow cytometry and examined through metagenomic sequencing. Homology-based sequence prediction indicated cyanophages, largely myoviruses, accounted for almost half the reads and provided insights into environmental infection events. T4-phage core-gene phylogenetic reconstruction revealed unique diversity among uncultured cyanophages and reference isolates resulting in identification of a new phylogenetic cluster. Genomic comparison of reference cyanophage strains S-SM2 and Syn1 with putative homologous contigs recovered from metagenomes provided evidence that gene insertion, deletion and recombination have occurred among, and are likely important for diversification of, natural populations. Identification of putative genetic exchange between cyanophage and non-cyanophage viruses, i.e. Micromonas virus and Pelagibacter phage, supports hypotheses related to a significant role for viruses in mediating transfer of genetic material between taxonomically diverse organisms with overlapping ecological niches.

Genomes and gene expression across light and productivity gradients in eastern subtropical Pacific microbial communities

Transitions in community genomic features and biogeochemical processes were examined in surface and subsurface chlorophyll maximum (SCM) microbial communities across a trophic gradient from mesotrophic waters near San Diego, California to the oligotrophic Pacific. Transect end points contrasted in thermocline depth, rates of nitrogen and CO2 uptake, new production and SCM light intensity. Relative to surface waters, bacterial SCM communities displayed greater genetic diversity and enrichment in putative sulfur oxidizers, multiple actinomycetes, low-light-adapted Prochlorococcus and cell-associated viruses. Metagenomic coverage was not correlated with transcriptional activity for several key taxa within Bacteria. Low-light-adapted Prochlorococcus, Synechococcus, and low abundance gamma-proteobacteria enriched in the>3.0-μm size fraction contributed disproportionally to global transcription. The abundance of these groups also correlated with community functions, such as primary production or nitrate uptake. In contrast, many of the most abundant bacterioplankton, including SAR11, SAR86, SAR112 and high-light-adapted Prochlorococcus, exhibited low levels of transcriptional activity and were uncorrelated with rate processes. Eukaryotes such as Haptophytes and non-photosynthetic Aveolates were prevalent in surface samples while Mamielles and Pelagophytes dominated the SCM. Metatranscriptomes generated with ribosomal RNA-depleted mRNA (total mRNA) coupled to in vitro polyadenylation compared with polyA-enriched mRNA revealed a trade-off in detection eukaryotic organelle and eukaryotic nuclear origin transcripts, respectively. Gene expression profiles of SCM eukaryote populations, highly similar in sequence identity to the model pelagophyte Pelagomonas sp. CCMP1756, suggest that pelagophytes are responsible for a majority of nitrate assimilation within the SCM.

Molecular mechanisms by which marine phytoplankton respond to their dynamic chemical environment

Marine scientists have long been interested in the interactions of marine phytoplankton with their chemical environments. Nutrient availability clearly controls carbon fixation on a global scale, but the interactions between phytoplankton and nutrients are complex and include both short-term responses (seconds to minutes) and longer-term evolutionary adaptations. This review outlines how genomics and functional genomics approaches are providing a better understanding of these complex interactions, especially for cyanobacteria and diatoms, for which the genome sequences of multiple model organisms are available. Transporters and related genes are emerging as the most likely candidates for biomarkers in stress-specific studies, but other genes are also possible candidates. One surprise has been the important role of horizontal gene transfer in mediating chemical-biological interactions.

Broad-host-range vector system for synthetic biology and biotechnology in cyanobacteria

Inspired by the developments of synthetic biology and the need for improved genetic tools to exploit cyanobacteria for the production of renewable bioproducts, we developed a versatile platform for the construction of broad-host-range vector systems. This platform includes the following features: (i) an efficient assembly strategy in which modules released from 3 to 4 donor plasmids or produced by polymerase chain reaction are assembled by isothermal assembly guided by short GC-rich overlap sequences. (ii) A growing library of molecular devices categorized in three major groups: (a) replication and chromosomal integration; (b) antibiotic resistance; (c) functional modules. These modules can be assembled in different combinations to construct a variety of autonomously replicating plasmids and suicide plasmids for gene knockout and knockin. (iii) A web service, the CYANO-VECTOR assembly portal, which was built to organize the various modules, facilitate the in silico construction of plasmids, and encourage the use of this system. This work also resulted in the construction of an improved broad-host-range replicon derived from RSF1010, which replicates in several phylogenetically distinct strains including a new experimental model strain Synechocystis sp. WHSyn, and the characterization of nine antibiotic cassettes, four reporter genes, four promoters, and a ribozyme-based insulator in several diverse cyanobacterial strains.

Genome of the halotolerant green alga Picochlorum sp. reveals strategies for thriving under fluctuating environmental conditions

An expected outcome of climate change is intensification of the global water cycle, which magnifies surface water fluxes, and consequently alters salinity patterns. It is therefore important to understand the adaptations and limits of microalgae to survive changing salinities. To this end, we sequenced the 13.5 Mbp genome of the halotolerant green alga Picochlorum SENEW3 (SE3) that was isolated from a brackish water pond subject to large seasonal salinity fluctuations. Picochlorum SE3 encodes 7367 genes, making it one of the smallest and most gene dense eukaryotic genomes known. Comparison with the pico-prasinophyte Ostreococcus tauri, a species with a limited range of salt tolerance, reveals the enrichment of transporters putatively involved in the salt stress response in Picochlorum SE3. Analysis of cultures and the protein complement highlight the metabolic flexibility of Picochlorum SE3 that encodes genes involved in urea metabolism, acetate assimilation and fermentation, acetoin production and glucose uptake, many of which form functional gene clusters. Twenty-four cases of horizontal gene transfer from bacterial sources were found in Picochlorum SE3 with these genes involved in stress adaptation including osmolyte production and growth promotion. Our results identify Picochlorum SE3 as a model for understanding microalgal adaptation to stressful, fluctuating environments.

The Marine Microbial Eukaryote Transcriptome Sequencing Project (MMETSP): illuminating the functional diversity of eukaryotic life in the oceans through transcriptome sequencing

Current sampling of genomic sequence data from eukaryotes is relatively poor, biased, and inadequate to address important questions about their biology, evolution, and ecology; this Community Page describes a resource of 700 transcriptomes from marine microbial eukaryotes to help understand their role in the world's oceans.

Microalgal assemblages in a poikilohaline pond

Microalgal strains for algal biofuels production in outdoor ponds will need to have high net growth rates under diverse environmental conditions. A small, variable salinity pond in the San Elijo Lagoon estuary in southern California was chosen to serve as a model pond due to its routinely high chlorophyll content. Profiles of microalgal assemblages from water samples collected from April 2011 to January 2012 were obtained by constructing 18S rDNA environmental clone libraries. Pond assemblages were found to be dominated by green algae Picochlorum sp. and Picocystis sp. throughout the year. Pigment analysis suggested that the two species contributed most of the chlorophyll a of the pond, which ranged from 21.9 to 664.3 μg · L(-1) with the Picocystis contribution increasing at higher salinities. However, changes of temperature, salinity or irradiance may have enabled a bloom of the diatom Chaetoceros sp. in June 2011. Isolates of these microalgae were obtained and their growth rates characterized as a function of temperature and salinity. Chaetoceros sp. had the highest growth rate over the temperature test range while it showed the most sensitivity to high salinity. All three strains showed the presence of lipid bodies during nitrogen starvation, suggesting they have potential as future biofuels strains.

Exposure to bloom-like concentrations of two marine Synechococcus cyanobacteria (strains CC9311 and CC9902) differentially alters fish behaviour

Coastal California experiences large-scale blooms of Synechococcus cyanobacteria, which are predicted to become more prevalent by the end of the 21st century as a result of global climate change. This study investigated whether exposure to bloom-like concentrations of two Synechococcus strains, CC9311 and CC9902, alters fish behaviour. Black perch (Embiotoca jacksoni) were exposed to Synechococcus strain CC9311 or CC9902 (1.5 × 10(6) cells ml(-1)) or to control seawater in experimental aquaria for 3 days. Fish movement inside a testing arena was then recorded and analysed using video camera-based motion-tracking software. Compared with control fish, fish exposed to CC9311 demonstrated a significant preference for the dark zone of the tank in the light-dark test, which is an indication of increased anxiety. Furthermore, fish exposed to CC9311 also had a statistically significant decrease in velocity and increase in immobility and they meandered more in comparison to control fish. There was a similar trend in velocity, immobility and meandering in fish exposed to CC9902, but there were no significant differences in behaviour or locomotion between this group and control fish. Identical results were obtained with a second batch of fish. Additionally, in this second trial we also investigated whether fish would recover after a 3 day period in seawater without cyanobacteria. Indeed, there were no longer any significant differences in behaviour among treatments, demonstrating that the sp. CC9311-induced alteration of behaviour is reversible. These results demonstrate that blooms of specific marine Synechococcus strains can induce differential sublethal effects in fish, namely alterations light-dark preference behaviour and motility.

Impact of DNA damaging agents on genome-wide transcriptional profiles in two marine Synechococcus species

Marine microorganisms, particularly those residing in coastal areas, may come in contact with any number of chemicals of environmental or xenobiotic origin. The sensitivity and response of marine cyanobacteria to such chemicals is, at present, poorly understood. We have looked at the transcriptional response of well characterized Synechococcus open ocean (WH8102) and coastal (CC9311) isolates to two DNA damaging agents, mitomycin C and ethidium bromide, using whole-genome expression microarrays. The coastal strain showed differential regulation of a larger proportion of its genome following “shock” treatment with each agent. Many of the orthologous genes in these strains, including those encoding sensor kinases, showed different transcriptional responses, with the CC9311 genes more likely to show significant changes in both treatments. While the overall response of each strain was considerably different, there were distinct transcriptional responses common to both strains observed for each DNA damaging agent, linked to the mode of action of each chemical. In both CC9311 and WH8102 there was evidence of SOS response induction under mitomycin C treatment, with genes recA, lexA and umuC significantly upregulated in this experiment but not under ethidium bromide treatment. Conversely, ethidium bromide treatment tended to result in upregulation of the DNA-directed RNA polymerase genes, not observed following mitomycin C treatment. Interestingly, a large number of genes residing on putative genomic island regions of each genome also showed significant upregulation under one or both chemical treatments.

Genomic island genes in a coastal marine Synechococcus strain confer enhanced tolerance to copper and oxidative stress

Highly variable regions called genomic islands are found in the genomes of marine picocyanobacteria, and have been predicted to be involved in niche adaptation and the ecological success of these microbes. These picocyanobacteria are typically highly sensitive to copper stress and thus, increased copper tolerance could confer a selective advantage under some conditions seen in the marine environment. Through targeted gene inactivation of genomic island genes that were known to be upregulated in response to copper stress in Synechococcus sp. strain CC9311, we found two genes (sync_1495 and sync_1217) conferred tolerance to both methyl viologen and copper stress in culture. The prevalence of one gene, sync_1495, was then investigated in natural samples, and had a predictable temporal variability in abundance at a coastal monitoring site with higher abundance in winter months. Together, this shows that genomic island genes can confer an adaptive advantage to specific stresses in marine Synechococcus, and may help structure their population diversity.

Selection in coastal Synechococcus (cyanobacteria) populations evaluated from environmental metagenomes

Environmental metagenomics provides snippets of genomic sequences from all organisms in an environmental sample and are an unprecedented resource of information for investigating microbial population genetics. Current analytical methods, however, are poorly equipped to handle metagenomic data, particularly of short, unlinked sequences. A custom analytical pipeline was developed to calculate dN/dS ratios, a common metric to evaluate the role of selection in the evolution of a gene, from environmental metagenomes sequenced using 454 technology of flow-sorted populations of marine Synechococcus, the dominant cyanobacteria in coastal environments. The large majority of genes (98%) have evolved under purifying selection (dN/dS<1). The metagenome sequence coverage of the reference genomes was not uniform and genes that were highly represented in the environment (i.e. high read coverage) tended to be more evolutionarily conserved. Of the genes that may have evolved under positive selection (dN/dS>1), 77 out of 83 (93%) were hypothetical. Notable among annotated genes, ribosomal protein L35 appears to be under positive selection in one Synechococcus population. Other annotated genes, in particular a possible porin, a large-conductance mechanosensitive channel, an ATP binding component of an ABC transporter, and a homologue of a pilus retraction protein had regions of the gene with elevated dN/dS. With the increasing use of next-generation sequencing in metagenomic investigations of microbial diversity and ecology, analytical methods need to accommodate the peculiarities of these data streams. By developing a means to analyze population diversity data from these environmental metagenomes, we have provided the first insight into the role of selection in the evolution of Synechococcus, a globally significant primary producer.

CHARACTERIZATION OF A FUNCTIONAL VANADIUM-DEPENDENT BROMOPEROXIDASE IN THE MARINE CYANOBACTERIUM SYNECHOCOCCUS SP. CC9311(1)

Vanadium-dependent bromoperoxidases (VBPOs) are characterized by the ability to oxidize halides using hydrogen peroxide. These enzymes are well-studied in eukaryotic macroalgae and are known to produce a variety of brominated secondary metabolites. Though genes have been annotated as VBPO in multiple prokaryotic genomes, they remain uncharacterized. The genome of the coastal marine cyanobacterium Synechococcus sp. CC9311 encodes a predicted VBPO (YP_731869.1, sync_2681), and in this study, we show that protein extracts from axenic cultures of Synechococcus possess bromoperoxidase activity, oxidizing bromide and iodide, but not chloride. In-gel activity assays of Synechococcus proteins separated using PAGE reveal a single band having VBPO activity. When sequenced via liquid chromatography/mass spectrometry/mass spectrometry (LC/MS/MS), peptides from the band aligned to the VBPO sequence predicted by the open reading frame (ORF) sync_2681. We show that a VBPO gene is present in a closely related strain, Synechococcus sp. WH8020, but not other clade I Synechococcus strains, consistent with recent horizontal transfer of the gene into Synechococcus. Diverse cyanobacterial-like VBPO genes were detected in a pelagic environment off the California coast using PCR. Investigation of functional VBPOs in unicellular cyanobacteria may lead to discovery of novel halogenated molecules and a better understanding of these organisms' chemical ecology and physiology.

Learning to read the oceans genomics of marine phytoplankton

The phytoplankton are key members of marine ecosystems, generating about half of global primary productivity, supporting valuable fisheries and regulating global biogeochemical cycles. Marine phytoplankton are phylogenetically diverse and are comprised of both prokaryotic and eukaryotic species. In the last decade, new insights have been gained into the ecology and evolution of these important organisms through whole genome sequencing projects and more recently, through both transcriptomics and targeted metagenomics approaches. Sequenced genomes of cyanobacteria are generally small, ranging in size from 1.8 to 9 million base pairs (Mbp). Eukaryotic genomes, in general, have a much larger size range and those that have been sequenced range from 12 to 57 Mbp. Whole genome sequencing projects have revealed key features of the evolutionary history of marine phytoplankton, their varied responses to environmental stress, their ability to scavenge and store nutrients and their unique ability to form elaborate cellular coverings. We have begun to learn how to read the 'language' of marine phytoplankton, as written in their DNA. Here, we review the ecological and evolutionary insights gained from whole genome sequencing projects, illustrate how these genomes are yielding information on marine natural products and informing nanotechnology as well as make suggestions for future directions in the field of marine phytoplankton genomics.

Amine oxidases of marine phytoplankton

Some phytoplankton utilized a novel mechanism for obtaining nitrogen from primary amines. They oxidized the primary amines to produce extracellular hydrogen peroxide and aldehydes and used the third reaction product, ammonium, as a nitrogen source. The specificity, regulation, inhibition by bromoethylamine, and potential dependence on copper of this process are described.

Computational prediction of the osmoregulation network in Synechococcus sp. WH8102

Background

Osmotic stress is caused by sudden changes in the impermeable solute concentration around a cell, which induces instantaneous water flow in or out of the cell to balance the concentration. Very little is known about the detailed response mechanism to osmotic stress in marine Synechococcus, one of the major oxygenic phototrophic cyanobacterial genera that contribute greatly to the global CO₂ fixation.

Results

We present here a computational study of the osmoregulation network in response to hyperosmotic stress of Synechococcus sp strain WH8102 using comparative genome analyses and computational prediction. In this study, we identified the key transporters, synthetases, signal sensor proteins and transcriptional regulator proteins, and found experimentally that of these proteins, 15 genes showed significantly changed expression levels under a mild hyperosmotic stress.

Conclusions

From the predicted network model, we have made a number of interesting observations about WH8102. Specifically, we found that (i) the organism likely uses glycine betaine as the major osmolyte, and others such as glucosylglycerol, glucosylglycerate, trehalose, sucrose and arginine as the minor osmolytes, making it efficient and adaptable to its changing environment; and (ii) σ³⁸, one of the seven types of σ factors, probably serves as a global regulator coordinating the osmoregulation network and the other relevant networks.

PtrA is required for coordinate regulation of gene expression during phosphate stress in a marine Synechococcus

Previous microarray analyses have shown a key role for the two-component system PhoBR (SYNW0947, SYNW0948) in the regulation of P transport and metabolism in the marine cyanobacterium Synechococcus sp. WH8102. However, there is some evidence that another regulator, SYNW1019 (PtrA), probably under the control of PhoBR, is involved in the response to P depletion. PtrA is a member of the cAMP receptor protein transcriptional regulator family that shows homology to NtcA, the global nitrogen regulator in cyanobacteria. To define the role of this regulator, we constructed a mutant by insertional inactivation and compared the physiology of wild-type Synechcococcus sp. WH8102 with the ptrA mutant under P-replete and P-stress conditions. In response to P stress the ptrA mutant failed to upregulate phosphatase activity. Microarrays and quantitative RT-PCR indicate that a subset of the Pho regulon is controlled by PtrA, including two phosphatases, a predicted phytase and a gene of unknown function psip1 (SYNW0165), all of which are highly upregulated during P limitation. Electrophoretic mobility shift assays indicate binding of overexpressed PtrA to promoter sequences upstream of the induced genes. This work suggests a two-tiered response to P depletion in this strain, the first being PhoB-dependent induction of high-affinity PO(4) transporters, and the second the PtrA-dependent induction of phosphatases for scavenging organic P. The levels of numerous other transcripts are also directly or indirectly influenced by PtrA, including those involved in cell-surface modification, metal uptake, photosynthesis, stress responses and other metabolic processes, which may indicate a wider role for PtrA in cellular regulation in marine picocyanobacteria.

Structure of compositionally simple lipopolysaccharide from marine synechococcus

Lipopolysaccharide (LPS) is the first defense against changing environmental factors for many bacteria. Here, we report the first structure of the LPS from cyanobacteria based on two strains of marine Synechococcus, WH8102 and CC9311. While enteric LPS contains some of the most complex carbohydrate residues in nature, the full-length versions of these cyanobacterial LPSs have neither heptose nor 3-deoxy-D-manno-octulosonic acid (Kdo) but instead 4-linked glucose as their main saccharide component, with low levels of glucosamine and galacturonic acid also present. Matrix-assisted laser desorption ionization mass spectrometry of the intact minimal core LPS reveals triacylated and tetraacylated structures having a heterogeneous mix of both hydroxylated and nonhydroxylated fatty acids connected to the diglucosamine backbone and a predominantly glucose outer core-like region for both strains. WH8102 incorporated rhamnose in this region as well, contributing to differences in sugar composition and possibly nutritional differences between the strains. In contrast to enteric lipid A, which can be liberated from LPS by mild acid hydrolysis, lipid A from these organisms could be produced by only two novel procedures: triethylamine-assisted periodate oxidation and acetolysis. The lipid A contains odd-chain hydroxylated fatty acids, lacks phosphate, and contains a single galacturonic acid. The LPS lacks any limulus amoebocyte lysate gelation activity. The highly simplified nature of LPSs from these organisms leads us to believe that they may represent either a primordial structure or an adaptation to the relatively higher salt and potentially growth-limiting phosphate levels in marine environments.

Whole-genome microarray analyses of Synechococcus-Vibrio interactions

Microbes live in diverse communities yet their physiologies are typically studied in axenic culture. To begin to address this dichotomy, whole-genome microarray analyses were used and revealed that several major metabolic pathways were affected in Synechococcus sp. WH8102, a model phototroph, when grown with Vibrio parahaemolyticus, a model heterotroph. In co-cultures with V. parahaemolyticus, although phosphate was not depleted, Synechococcus sp. WH8102 may have experienced phosphate stress since the expression of phosphate acquisition genes increased and alkaline phosphatase activity was higher than in monocultures. Expression of cell wall synthesis genes and the components of a zinc transporter were also upregulated. In contrast, a ferric uptake regulation (Fur) family gene was downregulated as were genes that encode proteins rich in iron or involved in detoxifying oxygen radicals. Nitrogen use may also have been affected in co-cultures as the gene expression changes share similarities with ammonia-grown Synechococcus. This study demonstrates the multiple impacts that interspecific microbial interactions can have on the physiology of a major primary producer and the importance of investigating microbial physiology from a community perspective.

Temporal variation of Synechococcus clades at a coastal Pacific Ocean monitoring site

Marine cyanobacteria from the genus Synechococcus are found throughout the world's oceans and are important contributors to global primary productivity and carbon cycling. Cultured isolates and environmental DNA clone libraries of Synechococcus have demonstrated the diversity of these microbes. However, the natural distribution of this diversity through space and time and the ecological significance of their distribution are still poorly understood. To understand the seasonal dynamics of Synechococcus diversity, we have developed a quantitative PCR strategy using the gene encoding as a subunit of DNA-dependent RNA polymerase (rpoC1) and applied it to a 3-year time series of surface samples from the Scripps Institution of Oceanography pier (La Jolla, CA, USA), a coastal site in the northeastern Pacific Ocean. Synechococcus from clades I and IV were dominant throughout the time series and correlated with total Synechococcus abundance. The relative abundance of these two dominant clades showed evidence of a seasonal cycle. Synechococcus from clade IV were typically more abundant, but those from clade I dominated during periods just before the annual spring bloom of Synechococcus. Synechococcus from clades II and III were absent during spring and early summer, but appeared at low abundances in late summer and winter possibly due to changes in circulation in the Southern California Bight. As the first long-term time series describing Synechococcus population diversity, these temporal dynamics were used to interpret the genetic/genomic diversity observed in the environment and the potential factors regulating their distribution.

Coastal Synechococcus metagenome reveals major roles for horizontal gene transfer and plasmids in population diversity

The extent to which cultured strains represent the genetic diversity of a population of microorganisms is poorly understood. Because they do not require culturing, metagenomic approaches have the potential to reveal the genetic diversity of the microbes actually present in an environment. From coastal California seawater, a complex and diverse environment, the marine cyanobacteria of the genus Synechococcus were enriched by flow cytometry-based sorting and the population metagenome was analysed with 454 sequencing technology. The sequence data were compared with model Synechococcus genomes, including those of two coastal strains, one isolated from the same and one from a very similar environment. The natural population metagenome had high sequence identity to most genes from the coastal model strains but diverged greatly from these genomes in multiple regions of atypical trinucleotide content that encoded diverse functions. These results can be explained by extensive horizontal gene transfer presumably with large differences in horizontally transferred genetic material between different strains. Some assembled contigs showed the presence of novel open reading frames not found in the model genomes, but these could not yet be unambiguously assigned to a Synechococcus clade. At least three distinct mobile DNA elements (plasmids) not found in model strain genomes were detected in the assembled contigs, suggesting for the first time their likely importance in marine cyanobacterial populations and possible role in horizontal gene transfer.

Diversity, function and evolution of genes coding for putative Ni-containing superoxide dismutases

We examined the phylogenetic distribution, functionality and evolution of the sodN gene family, which has been shown to code for a unique Ni-containing isoform of superoxide dismutase (Ni-SOD) in Streptomyces. Many of the putative sodN sequences retrieved from public domain genomic and metagenomic databases are quite divergent from structurally and functionally characterized Ni-SOD. Structural bioinformatics studies verified that the divergent members of the sodN protein family code for similar three-dimensional structures and identified evolutionarily conserved amino acid residues. Structural and biochemical studies of the N-terminus 'Ni-hook' motif coded for by the putative sodN sequences confirmed both Ni (II) ligating and superoxide dismutase activity. Both environmental and organismal genomes expanded the previously noted phylogenetic distribution of sodN, and the sequences form four well-separated clusters, with multiple subclusters. The phylogenetic distribution of sodN suggests that the gene has been acquired via horizontal gene transfer by numerous organisms of diverse phylogenetic background, including both Eukaryotes and Prokaryotes. The presence of sodN correlates with the genomic absence of the gene coding for Fe-SOD, a structurally and evolutionarily distinct isoform of SOD. Given the low levels of Fe found in the marine environment from where many sequences were attained, we suggest that the replacement of Fe-SOD with Ni-SOD may be an evolutionary adaptation to reduce iron requirements.

Characterization of a modular, cell-surface protein and identification of a new gene family in the diatom Thalassiosira pseudonana

We report the characterization of a cell-surface protein isolated from copper-stressed cells of the centric diatom Thalassiosira pseudonana Hasle and Heimdal (CCMP 1335). This protein has an apparent molecular weight of 100kDa and is highly acidic. The 100kDa protein (p100) sequence is comprised almost entirely of a novel domain termed TpRCR for T. pseudonana repetitive cysteine-rich domain, that is repeated 8 times and that contains conserved aromatic, acidic, and potential metal-binding amino acids. The analysis of the T. pseudonana genome suggests that p100 belongs to a large family of modular proteins that consist of a variable number of TpRCR domain repeats. Based on cell surface biotinylation and antibody data, p100 appears to migrate more rapidly with SDS-PAGE when extracted from cells exposed to high levels of copper; however, the discovery of a large family of TpRCR domain-containing proteins leaves open the possibility that the antibody may be cross-reacting with members of this protein family that are responding differently to copper. The response of the gene encoding p100 at the mRNA level during synchronized progression through the normal cell cycle is similar to previously characterized genes in T. pseudonana encoding cell wall proteins called silaffins.

The tiny eukaryote Ostreococcus provides genomic insights into the paradox of plankton speciation

The smallest known eukaryotes, at approximately 1-mum diameter, are Ostreococcus tauri and related species of marine phytoplankton. The genome of Ostreococcus lucimarinus has been completed and compared with that of O. tauri. This comparison reveals surprising differences across orthologous chromosomes in the two species from highly syntenic chromosomes in most cases to chromosomes with almost no similarity. Species divergence in these phytoplankton is occurring through multiple mechanisms acting differently on different chromosomes and likely including acquisition of new genes through horizontal gene transfer. We speculate that this latter process may be involved in altering the cell-surface characteristics of each species. In addition, the genome of O. lucimarinus provides insights into the unique metal metabolism of these organisms, which are predicted to have a large number of selenocysteine-containing proteins. Selenoenzymes are more catalytically active than similar enzymes lacking selenium, and thus the cell may require less of that protein. As reported here, selenoenzymes, novel fusion proteins, and loss of some major protein families including ones associated with chromatin are likely important adaptations for achieving a small cell size.

Modern proteomes contain putative imprints of ancient shifts in trace metal geochemistry

Because of the rise in atmospheric oxygen 2.3 billion years ago (Gya) and the subsequent changes in oceanic redox state over the last 2.3-1 Gya, trace metal bioavailability in marine environments has changed dramatically. Although theorized to have influenced the biological usage of metals leaving discernable genomic signals, a thorough and quantitative test of this hypothesis has been lacking. Using structural bioinformatics and whole-genome sequences, the Fe-, Zn-, Mn-, and Co-binding metallomes of 23 Archaea, 233 Bacteria, and 57 Eukarya were constructed. These metallomes reveal that the overall abundances of these metal-binding structures scale to proteome size as power laws with a unique set of slopes for each Superkingdom of Life. The differences in the power describing the abundances of Fe-, Mn-, Zn-, and Co-binding proteins in the proteomes of Prokaryotes and Eukaryotes are similar to the theorized changes in the abundances of these metals after the oxygenation of oceanic deep waters. This phenomenon suggests that Prokarya and Eukarya evolved in anoxic and oxic environments, respectively, a hypothesis further supported by structures and functions of Fe-binding proteins in each Superkingdom. Also observed is a proliferation in the diversity of Zn-binding protein structures involved in protein-DNA and protein-protein interactions within Eukarya, an event unlikely to occur in either an anoxic or euxinic environment where Zn concentrations would be vanishingly low. We hypothesize that these conserved trends are proteomic imprints of changes in trace metal bioavailability in the ancient ocean that highlight a major evolutionary shift in biological trace metal usage.

Genome sequence of Synechococcus CC9311: Insights into adaptation to a coastal environment

Coastal aquatic environments are typically more highly productive and dynamic than open ocean ones. Despite these differences, cyanobacteria from the genus Synechococcus are important primary producers in both types of ecosystems. We have found that the genome of a coastal cyanobacterium, Synechococcus sp. strain CC9311, has significant differences from an open ocean strain, Synechococcus sp. strain WH8102, and these are consistent with the differences between their respective environments. CC9311 has a greater capacity to sense and respond to changes in its (coastal) environment. It has a much larger capacity to transport, store, use, or export metals, especially iron and copper. In contrast, phosphate acquisition seems less important, consistent with the higher concentration of phosphate in coastal environments. CC9311 is predicted to have differences in its outer membrane lipopolysaccharide, and this may be characteristic of the speciation of some cyanobacterial groups. In addition, the types of potentially horizontally transferred genes are markedly different between the coastal and open ocean genomes and suggest a more prominent role for phages in horizontal gene transfer in oligotrophic environments.

Gene expression induced by copper stress in the diatom Thalassiosira pseudonana

Utilizing a PCR-based subtractive cDNA approach, we demonstrated that the marine diatom Thalassiosira pseudonana exhibits a rapid response at the gene level to elevated concentrations of copper and that this response attenuates over 24 h of continuous exposure. A total of 16 copper-induced genes were identified, 11 of which were completely novel; however, many of the predicted amino acid sequences had characteristics suggestive of roles in ameliorating copper toxicity. Most of the novel genes were not equivalently induced by H2O2- or Cd-induced stress, indicating specificity in response. Two genes that could be assigned functions based on homology were also induced under conditions of general cellular stress. Half of the identified genes were located within two inverted repeats in the genome, and novel genes in one inverted repeat had mRNA levels induced by approximately 500- to 2,000-fold by exposure to copper for 1 h. Additionally, some of the inverted repeat genes demonstrated a dose-dependent response to Cu, but not Cd, and appear to belong to a multigene family. This multigene family may be the diatom functional homolog of metallothioneins.

Computational inference and experimental validation of the nitrogen assimilation regulatory network in cyanobacterium Synechococcus sp. WH 8102

Deciphering the regulatory networks encoded in the genome of an organism represents one of the most interesting and challenging tasks in the post-genome sequencing era. As an example of this problem, we have predicted a detailed model for the nitrogen assimilation network in cyanobacterium Synechococcus sp. WH 8102 (WH8102) using a computational protocol based on comparative genomics analysis and mining experimental data from related organisms that are relatively well studied. This computational model is in excellent agreement with the microarray gene expression data collected under ammonium-rich versus nitrate-rich growth conditions, suggesting that our computational protocol is capable of predicting biological pathways/networks with high accuracy. We then refined the computational model using the microarray data, and proposed a new model for the nitrogen assimilation network in WH8102. An intriguing discovery from this study is that nitrogen assimilation affects the expression of many genes involved in photosynthesis, suggesting a tight coordination between nitrogen assimilation and photosynthesis processes. Moreover, for some of these genes, this coordination is probably mediated by NtcA through the canonical NtcA promoters in their regulatory regions.

The Genome of the Diatom Thalassiosira Pseudonana: Ecology, Evolution, and Metabolism

Diatoms are unicellular algae with plastids acquired by secondary endosymbiosis. They are responsible for approximately 20% of global carbon fixation. We report the 34 million-base pair draft nuclear genome of the marine diatom Thalassiosira pseudonana and its 129 thousand-base pair plastid and 44 thousand-base pair mitochondrial genomes. Sequence and optical restriction mapping revealed 24 diploid nuclear chromosomes. We identified novel genes for silicic acid transport and formation of silica-based cell walls, high-affinity iron uptake, biosynthetic enzymes for several types of polyunsaturated fatty acids, use of a range of nitrogenous compounds, and a complete urea cycle, all attributes that allow diatoms to prosper in aquatic environments.

Operon prediction by comparative genomics: an application to the Synechococcus sp. WH8102 genome

We present a computational method for operon prediction based on a comparative genomics approach. A group of consecutive genes is considered as a candidate operon if both their gene sequences and functions are conserved across several phylogenetically related genomes. In addition, various supporting data for operons are also collected through the application of public domain computer programs, and used in our prediction method. These include the prediction of conserved gene functions, promoter motifs and terminators. An apparent advantage of our approach over other operon prediction methods is that it does not require many experimental data (such as gene expression data and pathway data) as input. This feature makes it applicable to many newly sequenced genomes that do not have extensive experimental information. In order to validate our prediction, we have tested the method on Escherichia coli K12, in which operon structures have been extensively studied, through a comparative analysis against Haemophilus influenzae Rd and Salmonella typhimurium LT2. Our method successfully predicted most of the 237 known operons. After this initial validation, we then applied the method to a newly sequenced and annotated microbial genome, Synechococcus sp. WH8102, through a comparative genome analysis with two other cyanobacterial genomes, Prochlorococcus marinus sp. MED4 and P.marinus sp. MIT9313. Our results are consistent with previously reported results and statistics on operons in the literature.

The genome of a motile marine Synechococcus

Marine unicellular cyanobacteria are responsible for an estimated 20-40% of chlorophyll biomass and carbon fixation in the oceans. Here we have sequenced and analysed the 2.4-megabase genome of Synechococcus sp. strain WH8102, revealing some of the ways that these organisms have adapted to their largely oligotrophic environment. WH8102 uses organic nitrogen and phosphorus sources and more sodium-dependent transporters than a model freshwater cyanobacterium. Furthermore, it seems to have adopted strategies for conserving limited iron stores by using nickel and cobalt in some enzymes, has reduced its regulatory machinery (consistent with the fact that the open ocean constitutes a far more constant and buffered environment than fresh water), and has evolved a unique type of swimming motility. The genome of WH8102 seems to have been greatly influenced by horizontal gene transfer, partially through phages. The genetic material contributed by horizontal gene transfer includes genes involved in the modification of the cell surface and in swimming motility. On the basis of its genome, WH8102 is more of a generalist than two related marine cyanobacteria.

Carbon sequestration in Synechococcus Sp.: from molecular machines to hierarchical modeling

The U.S. Department of Energy recently announced the first five grants for the Genomes to Life (GTL) Program. The goal of this program is to "achieve the most far-reaching of all biological goals: a fundamental, comprehensive, and systematic understanding of life." While more information about the program can be found at the GTL website (www.doegenomestolife.org), this paper provides an overview of one of the five GTL projects funded, "Carbon Sequestration in Synechococcus Sp.: From Molecular Machines to Hierarchical Modeling." This project is a combined experimental and computational effort emphasizing developing, prototyping, and applying new computational tools and methods to elucidate the biochemical mechanisms of the carbon sequestration of Synechococcus Sp., an abundant marine cyanobacteria known to play an important role in the global carbon cycle. Understanding, predicting, and perhaps manipulating carbon fixation in the oceans has long been a major focus of biological oceanography and has more recently been of interest to a broader audience of scientists and policy makers. It is clear that the oceanic sinks and sources of CO(2) are important terms in the global environmental response to anthropogenic atmospheric inputs of CO(2) and that oceanic microorganisms play a key role in this response. However, the relationship between this global phenomenon and the biochemical mechanisms of carbon fixation in these microorganisms is poorly understood. The project includes five subprojects: an experimental investigation, three computational biology efforts, and a fifth which deals with addressing computational infrastructure challenges of relevance to this project and the Genomes to Life program as a whole. Our experimental effort is designed to provide biology and data to drive the computational efforts and includes significant investment in developing new experimental methods for uncovering protein partners, characterizing protein complexes, identifying new binding domains. We will also develop and apply new data measurement and statistical methods for analyzing microarray experiments. Our computational efforts include coupling molecular simulation methods with knowledge discovery from diverse biological data sets for high-throughput discovery and characterization of protein-protein complexes and developing a set of novel capabilities for inference of regulatory pathways in microbial genomes across multiple sources of information through the integration of computational and experimental technologies. These capabilities will be applied to Synechococcus regulatory pathways to characterize their interaction map and identify component proteins in these pathways. We will also investigate methods for combining experimental and computational results with visualization and natural language tools to accelerate discovery of regulatory pathways. Furthermore, given that the ultimate goal of this effort is to develop a systems-level of understanding of how the Synechococcus genome affects carbon fixation at the global scale, we will develop and apply a set of tools for capturing the carbon fixation behavior of complex of Synechococcus at different levels of resolution. Finally, because the explosion of data being produced by high-throughput experiments requires data analysis and models which are more computationally complex, more heterogeneous, and require coupling to ever increasing amounts of experimentally obtained data in varying formats, we have also established a companion computational infrastructure to support this effort as well as the Genomes to Life program as a whole.

Computational inference of regulatory pathways in microbes: an application to phosphorus assimilation pathways in Synechococcus sp. WH8102

We present a computational protocol for inference of regulatory and signaling pathways in a microbial cell, through literature search, mining "high-throughput'' biological data of various types, and computer-assisted human inference. This protocol consists of four key components: (a) construction of template pathways for microbial organisms related to the target genome, which either have been extensively studied and/or have a significant amount of (relevant) experimental data, (b) inference of initial pathway models for the target genome, through combining the template pathway models and target genome-specific information, (c) refinement and expansion of the initial pathway models through applications of various data mining tools, including phylogenetic profile analysis, inference of protein-protein interactions, and prediction of transcription factor binding sites, and (d) validation and refinement of the pathway models using pathway-specific experimental data or other information. To demonstrate the effectiveness of this procedure, we have applied it to the construction of the phosphorus assimilation pathways in cyanobacterium sp. WH8102. We present, in this paper, a model of the core components of this pathway.

Chromatic adaptation in marine Synechococcus strains

Characterization of two genetically distinct groups of marine Synechococcus sp. strains shows that one, but not the other, increases its phycourobilin/phycoerythrobilin chromophore ratio when growing in blue light. This ability of at least some marine Synechococcus strains to chromatically adapt may help explain their greater abundance in particular ocean environments than cyanobacteria of the genus Prochlorococcus.

Swimming marine Synechococcus strains with widely different photosynthetic pigment ratios form a monophyletic group

Unicellular marine cyanobacteria are ubiquitous in both coastal and oligotrophic regimes. The contribution of these organisms to primary production and nutrient cycling is substantial on a global scale. Natural populations of marine Synechococcus strains include multiple genetic lineages, but the link, if any, between unique phenotypic traits and specific genetic groups is still not understood. We studied the genetic diversity (as determined by the DNA-dependent RNA polymerase rpoC1 gene sequence) of a set of marine Synechococcus isolates that are able to swim. Our results show that these isolates form a monophyletic group. This finding represents the first example of correspondence between a physiological trait and a phylogenetic group in marine Synechococcus. In contrast, the phycourobilin (PUB)/phycoerythrobilin (PEB) pigment ratios of members of the motile clade varied considerably. An isolate obtained from the California Current (strain CC9703) displayed a pigment signature identical to that of nonmotile strain WH7803, which is considered a model for low-PUB/PEB-ratio strains, whereas several motile strains had higher PUB/PEB ratios than strain WH8103, which is considered a model for high-PUB/PEB-ratio strains. These findings indicate that the PUB/PEB pigment ratio is not a useful characteristic for defining phylogenetic groups of marine Synechococcus strains.

The marine cyanobacterium Synechococcus sp. WH7805 requires urease (urea amidohydrolase, EC 3.5.1.5) to utilize urea as a nitrogen source: molecular-genetic and biochemical analysis of the enzyme

Cyanobacteria assigned to the genus Synechococcus are an important component of oligotrophic marine ecosystems, where their growth may be constrained by low availability of fixed nitrogen. Urea appears to be a major nitrogen resource in the sea, but little molecular information exists about its utilization by marine organisms, including Synechococcus. Oligonucleotide primers were used to amplify a conserved fragment of the urease (urea amidohydrolase, EC 3.5.1.5) coding region from cyanobacteria. A 5.7 kbp region of the genome of the unicellular marine cyanobacterium Synechococcus sp. strain WH7805 was then cloned, and genes encoding three urease structural subunits and four urease accessory proteins were sequenced and identified by homology. The WH7805 urease had a predicted subunit composition typical of bacterial ureases, but the organization of the WH7805 urease genes was unique. Biochemical characteristics of the WH7805 urease enzyme were consistent with the predictions of the sequence data. Physiological data and sequence analysis both suggested that the urease operon may be nitrogen-regulated by the ntcA system in WH7805. Inactivation of the large subunit of urease, ureC, prevented WH7805 and Synechococcus WH8102 from growing on urea, demonstrating that the urease genes cloned are essential to the ability of these cyanobacteria to utilize urea as a nitrogen source.