Search

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.

Although gene regulation can occur at both transcriptional and epigenetic (microRNA) levels, combined transcriptomic and microRNAomic responses to environmental stress are still largely unexplored for marine plankton. Here, we conducted transcriptome and microRNAome sequencing for Prorocentrum donghaiense to understand the molecular mechanisms by which this dinoflagellate copes with phosphorus (P) deficiency. Under P-depleted conditions, G1/S specific cyclin gene was markedly downregulated, consistent with growth inhibition, and genes related to dissolved organic phosphorus (DOP) hydrolysis, carbon fixation, nitrate assimilation, glycolysis, and cellular motility were upregulated. The elevated expression of ATP-generating genes (for example, rhodopsin) and ATP-consuming genes suggests some metabolic reconfiguration towards accelerated ATP recycling under P deficiency. MicroRNAome sequencing revealed 17 microRNAs, potentially regulating 3268 protein-coding genes. Functional enrichment analysis of these microRNA-targeted genes predicted decreases in sulfatide (sulfolipid) catabolism under P deficiency. Strikingly, we detected a significant increase in sulfolipid sulfatide content (but not in sulphoquinovosyldiacylglycerol content) and its biosynthesis gene expression, indicating a different sulfolipid-substituting-phospholipid mechanism in this dinoflagellate than other phytoplankters studied previously. Taken together, our integrative transcriptomic and microRNAomic analyses show that enhanced DOP utilization, accelerated ATP cycling and repressed sulfolipid degradation constitute a comprehensive strategy to cope with P deficiency in a model dinoflagellate.

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.

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.

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., (© The Author(s) 2014. Published by Oxford University Press on behalf of Nucleic Acids Research.)

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.

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.

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.

Primary productivity of open ocean environments, such as those inhabited by marine picocyanobacteria, is often limited by low inorganic phosphate (P). To observe how these organisms cope with P starvation, we constructed a full genome microarray for Synechococcus sp. WH8102 and compared differences in gene expression under P-replete and P-limited growth conditions, including both early P stress, during extracellular alkaline phosphatase induction, and late P stress. A total of 36 genes showed significant upregulation (>log(2) fold) whereas 23 genes were highly downregulated at the early time point; however, these changes in expression were maintained during late P stress for only 5 of the upregulated genes. Knockout mutants were constructed for genes SYNW0947 and SYNW0948, comprising a two-component regulator hypothesized to have a key function in regulating P metabolism. A high degree of overlap in the sets of genes affected by P stress conditions and in the knockout mutants supports this hypothesis; however, there is some indication that other regulators may be involved in this response in Synechococcus sp. WH8102. Consistent with what has been observed in many other cyanobacteria, the Pho regulon of this strain is comprised largely of genes for alkaline phosphatases, P transport or P metabolism. Interestingly, however, the exact composition and arrangement of the Pho regulon appears highly variable in marine cyanobacteria.

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.

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.