Search

Many heterotrophic bacteria are known to release extracellular vesicles, facilitating interactions between cells and their environment from a distance. Vesicle production has not been described in photoautotrophs, however, and the prevalence and characteristics of vesicles in natural ecosystems is unknown. Here, we report that cultures of Prochlorococcus, a numerically dominant marine cyanobacterium, continuously release lipid vesicles containing proteins, DNA, and RNA. We also show that vesicles carrying DNA from diverse bacteria are abundant in coastal and open-ocean seawater samples. Prochlorococcus vesicles can support the growth of heterotrophic bacterial cultures, which implicates these structures in marine carbon flux. The ability of vesicles to deliver diverse compounds in discrete packages adds another layer of complexity to the flow of information, energy, and biomolecules in marine microbial communities.

Viruses that infect marine cyanobacteria-cyanophages-often carry genes with orthologs in their cyanobacterial hosts, and the frequency of these genes can vary with habitat. To explore habitat-influenced genomic diversity more deeply, we used the genomes of 28 cultured cyanomyoviruses as references to identify phage genes in three ocean habitats. Only about 6-11% of genes were consistently observed in the wild, revealing high gene-content variability in these populations. Numerous shared phage/host genes differed in relative frequency between environments, including genes related to phosphorous acquisition, photorespiration, photosynthesis and the pentose phosphate pathway, possibly reflecting environmental selection for these genes in cyanomyovirus genomes. The strongest emergent signal was related to phosphorous availability; a higher fraction of genomes from relatively low-phosphorus environments-the Sargasso and Mediterranean Sea-contained host-like phosphorus assimilation genes compared with those from the N. Pacific Gyre. These genes are known to be upregulated when the host is phosphorous starved, a response mediated by pho box motifs in phage genomes that bind a host regulatory protein. Eleven cyanomyoviruses have predicted pho boxes upstream of the phosphate-acquisition genes pstS and phoA; eight of these have a conserved cyanophage-specific gene (PhCOG173) between the pho box and pstS. PhCOG173 is also found upstream of other shared phage/host genes, suggesting a unique regulatory role. Pho boxes are found upstream of high light-inducible (hli) genes in cyanomyoviruses, suggesting that this motif may have a broader role than regulating phosphorous-stress responses in infected hosts or that these hlis are involved in the phosphorous-stress response.

The marine cyanobacteria Prochlorococcus and Synechococcus are highly abundant in the global oceans, as are the cyanophage with which they co-evolve. While genomic analyses have been relatively extensive for cyanomyoviruses, only three cyanopodoviruses isolated on marine cyanobacteria have been sequenced. Here we present nine new cyanopodovirus genomes, and analyse them in the context of the broader group. The genomes range from 42.2 to 47.7 kb, with G+C contents consistent with those of their hosts. They share 12 core genes, and the pan-genome is not close to being fully sampled. The genomes contain three variable island regions, with the most hypervariable genes concentrated at one end of the genome. Concatenated core-gene phylogeny clusters all but one of the phage into three distinct groups (MPP-A and two discrete clades within MPP-B). The outlier, P-RSP2, has the smallest genome and lacks RNA polymerase, a hallmark of the Autographivirinae subfamily. The phage in group MPP-B contain photosynthesis and carbon metabolism associated genes, while group MPP-A and the outlier P-RSP2 do not, suggesting different constraints on their lytic cycles. Four of the phage encode integrases and three have a host integration signature. Metagenomic analyses reveal that cyanopodoviruses may be more abundant in the oceans than previously thought., (© 2012 Society for Applied Microbiology and Blackwell Publishing Ltd.)

Prochlorococcus is the numerically dominant photosynthetic organism throughout much of the world's oceans, yet little is known about the ecology and genetic diversity of populations inhabiting tropical waters. To help close this gap, we examined natural Prochlorococcus communities in the tropical Pacific Ocean using a single-cell whole-genome amplification and sequencing. Analysis of the gene content of just 10 single cells from these waters added 394 new genes to the Prochlorococcus pan-genome--that is, genes never before seen in a Prochlorococcus cell. Analysis of marker genes, including the ribosomal internal transcribed sequence, from dozens of individual cells revealed several representatives from two uncultivated clades of Prochlorococcus previously identified as HNLC1 and HNLC2. While the HNLC clades can dominate Prochlorococcus communities under certain conditions, their overall geographic distribution was highly restricted compared with other clades of Prochlorococcus. In the Atlantic and Pacific oceans, these clades were only found in warm waters with low Fe and high inorganic P levels. Genomic analysis suggests that at least one of these clades thrives in low Fe environments by scavenging organic-bound Fe, a process previously unknown in Prochlorococcus. Furthermore, the capacity to utilize organic-bound Fe appears to have been acquired horizontally and may be exchanged among other clades of Prochlorococcus. Finally, one of the single Prochlorococcus cells sequenced contained a partial genome of what appears to be a prophage integrated into the genome.

Primary productivity in the ocean's oligotrophic regions is often limited by phosphorus (P) availability. In low phosphate environments, the prevalence of many genes involved in P acquisition is elevated, suggesting that the ability to effectively access diverse P sources is advantageous for organisms inhabiting these regions. Prochlorococcus, the numerically dominant primary producer in the oligotrophic ocean, encodes high-affinity P transporters, P regulatory proteins and enzymes for organic phosphate utilization, but its ability to use reduced P compounds has not been previously demonstrated. Because Prochlorococcus strain MIT9301 encodes genes similar to phnY and phnZ, which constitute a novel marine bacterial 2-aminoethylphosphonate (2-AEPn) utilization pathway, it has been suggested that this organism might use 2-AEPn as an alternative P source. We show here that although MIT9301 was unable to use 2-AEPn as a sole P source under standard culture conditions, it was able to use phosphite. Phosphite utilization by MIT9301 appears to be mediated by an NAD-dependent phosphite dehydrogenase encoded by ptxD. We show that phosphite utilization genes are present in diverse marine microbes and that their abundance is higher in low-P waters. These results strongly suggest that phosphite represents a previously unrecognized component of the marine P cycle., (© 2011 Society for Applied Microbiology and Blackwell Publishing Ltd.)

Phosphorus (P) availability, which often limits productivity in marine ecosystems, shapes the P-acquisition gene content of the marine cyanobacteria Prochlorococcus [1-4] and its viruses (cyanophages). As in other bacteria, in Prochlorococcus these genes are regulated by the PhoR/PhoB two-component regulatory system that is used to sense and respond to P availability and is typical of signal transduction systems found in diverse organisms. Replication of cyanophage genomes requires a significant amount of P, and therefore these phages could gain a fitness advantage by influencing host P acquisition in P-limited environments. Here we show that the transcription of a phage-encoded high-affinity phosphate-binding protein gene (pstS) and alkaline phosphatase gene (phoA)-both of which have host orthologs-is elevated when the phages are infecting host cells that are P starved, relative to P-replete control cells. We further show that the phage versions of these genes are regulated by the host's PhoR/PhoB system. This not only extends this fundamental signaling mechanism to viruses but is also the first example of regulation of lytic phage genes by nutrient limitation in the host. As such, it reveals an important new dimension of the intimate coevolution of phage, host, and environment in the world's oceans., (Copyright © 2012 Elsevier Ltd. All rights reserved.)

Background: Growth of the ocean's most abundant primary producer, the cyanobacterium Prochlorococcus, is tightly synchronized to the natural 24-hour light-dark cycle. We sought to quantify the relationship between transcriptome and proteome dynamics that underlie this obligate photoautotroph's highly choreographed response to the daily oscillation in energy supply., Methodology/principal Findings: Using RNA-sequencing transcriptomics and mass spectrometry-based quantitative proteomics, we measured timecourses of paired mRNA-protein abundances for 312 genes every 2 hours over a light-dark cycle. These temporal expression patterns reveal strong oscillations in transcript abundance that are broadly damped at the protein level, with mRNA levels varying on average 2.3 times more than the corresponding protein. The single strongest observed protein-level oscillation is in a ribonucleotide reductase, which may reflect a defense strategy against phage infection. The peak in abundance of most proteins also lags that of their transcript by 2-8 hours, and the two are completely antiphase for some genes. While abundant antisense RNA was detected, it apparently does not account for the observed divergences between expression levels. The redirection of flux through central carbon metabolism from daytime carbon fixation to nighttime respiration is associated with quite small changes in relative enzyme abundances., Conclusions/significance: Our results indicate that expression responses to periodic stimuli that are common in natural ecosystems (such as the diel cycle) can diverge significantly between the mRNA and protein levels. Protein expression patterns that are distinct from those of cognate mRNA have implications for the interpretation of transcriptome and metatranscriptome data in terms of cellular metabolism and its biogeochemical impact.

ProPortal (http://proportal.mit.edu/) is a database containing genomic, metagenomic, transcriptomic and field data for the marine cyanobacterium Prochlorococcus. Our goal is to provide a source of cross-referenced data across multiple scales of biological organization--from the genome to the ecosystem--embracing the full diversity of ecotypic variation within this microbial taxon, its sister group, Synechococcus and phage that infect them. The site currently contains the genomes of 13 Prochlorococcus strains, 11 Synechococcus strains and 28 cyanophage strains that infect one or both groups. Cyanobacterial and cyanophage genes are clustered into orthologous groups that can be accessed by keyword search or through a genome browser. Users can also identify orthologous gene clusters shared by cyanobacterial and cyanophage genomes. Gene expression data for Prochlorococcus ecotypes MED4 and MIT9313 allow users to identify genes that are up or downregulated in response to environmental stressors. In addition, the transcriptome in synchronized cells grown on a 24-h light-dark cycle reveals the choreography of gene expression in cells in a 'natural' state. Metagenomic sequences from the Global Ocean Survey from Prochlorococcus, Synechococcus and phage genomes are archived so users can examine the differences between populations from diverse habitats. Finally, an example of cyanobacterial population data from the field is included.

The marine cyanobacterium Prochlorococcus, the smallest and most abundant oxygenic phototroph, has an extremely streamlined genome and a high rate of protein evolution. High-light adapted strains of Prochlorococcus in particular have seemingly inadequate DNA repair systems, raising the possibility that inadequate repair may lead to high mutation rates. Prochlorococcus mutation rates have been difficult to determine, in part because traditional methods involving quantifying colonies on solid selective media are not straightforward for this organism. Here we used a liquid dilution method to measure the approximate number of antibiotic-resistant mutants in liquid cultures of Prochlorococcus strains previously unexposed to antibiotic selection. Several antibiotics for which resistance in other bacteria is known to result from a single base pair change were used. The resulting frequencies of antibiotic resistance in Prochlorococcus cultures allowed us to then estimate maximum spontaneous mutation rates, which were similar to those in organisms such as E. coli (∼5.4 × 10(-7) per gene per generation). Therefore, despite the lack of some DNA repair genes, it appears unlikely that the Prochlorcoccus genomes studied here are currently being shaped by unusually high mutation rates., (© 2011 Society for Applied Microbiology and Blackwell Publishing Ltd.)

Prochlorococcus contributes significantly to ocean primary productivity. The link between primary productivity and iron in specific ocean regions is well established and iron limitation of Prochlorococcus cell division rates in these regions has been shown. However, the extent of ecotypic variation in iron metabolism among Prochlorococcus and the molecular basis for differences is not understood. Here, we examine the growth and transcriptional response of Prochlorococcus strains, MED4 and MIT9313, to changing iron concentrations. During steady state, MIT9313 sustains growth at an order-of-magnitude lower iron concentration than MED4. To explore this difference, we measured the whole-genome transcriptional response of each strain to abrupt iron starvation and rescue. Only four of the 1159 orthologs of MED4 and MIT9313 were differentially expressed in response to iron in both strains. However, in each strain, the expression of over a hundred additional genes changed, many of which are in labile genomic regions, suggesting a role for lateral gene transfer in establishing diversity of iron metabolism among Prochlorococcus. Furthermore, we found that MED4 lacks three genes near the iron-deficiency-induced gene (idiA) that are present and induced by iron stress in MIT9313. These genes are interesting targets for studying the adaptation of natural Prochlorococcus assemblages to local iron conditions as they show more diversity than other genomic regions in environmental metagenomic databases.

Cyanophages infecting the marine cyanobacteria Prochlorococcus and Synechococcus encode and express genes for the photosynthetic light reactions. Sequenced cyanophage genomes lack Calvin cycle genes, however, suggesting that photosynthetic energy harvested via phage proteins is not used for carbon fixation. We report here that cyanophages carry and express a Calvin cycle inhibitor, CP12, whose host homologue directs carbon flux from the Calvin cycle to the pentose phosphate pathway (PPP). Phage CP12 was coexpressed with phage genes involved in the light reactions, deoxynucleotide biosynthesis, and the PPP, including a transaldolase gene that is the most prevalent PPP gene in cyanophages. Phage transaldolase was purified to homogeneity from several strains and shown to be functional in vitro, suggesting that it might facilitate increased flux through this key reaction in the host PPP, augmenting production of NADPH and ribose 5-phosphate. Kinetic measurements of phage and host transaldolases revealed that the phage enzymes have k(cat)/K(m) values only approximately one third of the corresponding host enzymes. The lower efficiency of phage transaldolase may be a tradeoff for other selective advantages such as reduced gene size: we show that more than half of host-like cyanophage genes are significantly shorter than their host homologues. Consistent with decreased Calvin cycle activity and increased PPP and light reaction activity under infection, the host NADPH/NADP ratio increased two-fold in infected cells. We propose that phage-augmented NADPH production fuels deoxynucleotide biosynthesis for phage replication, and that the selection pressures molding phage genomes involve fitness advantages conferred through mobilization of host energy stores.

Interactions between microorganisms shape microbial ecosystems. Systematic studies of mixed microbes in co-culture have revealed widespread potential for growth inhibition among marine heterotrophic bacteria, but similar synoptic studies have not been done with autotroph/heterotroph pairs, nor have precise descriptions of the temporal evolution of interactions been attempted in a high-throughput system. Here, we describe patterns in the outcome of pair-wise co-cultures between two ecologically distinct, yet closely related, strains of the marine cyanobacterium Prochlorococcus and hundreds of heterotrophic marine bacteria. Co-culture with the collection of heterotrophic strains influenced the growth of Prochlorococcus strain MIT9313 much more than that of strain MED4, reflected both in the number of different types of interactions and in the magnitude of the effect of co-culture on various culture parameters. Enhancing interactions, where the presence of heterotrophic bacteria caused Prochlorococcus to grow faster and reach a higher final culture chlorophyll fluorescence, were much more common than antagonistic ones, and for a selected number of cases were shown to be mediated by diffusible compounds. In contrast, for one case at least, temporary inhibition of Prochlorococcus MIT9313 appeared to require close cellular proximity. Bacterial strains whose 16S gene sequences differed by 1-2% tended to have similar effects on MIT9313, suggesting that the patterns of inhibition and enhancement in co-culture observed here are due to phylogenetically cohesive traits of these heterotrophs.

T4-like myoviruses are ubiquitous, and their genes are among the most abundant documented in ocean systems. Here we compare 26 T4-like genomes, including 10 from non-cyanobacterial myoviruses, and 16 from marine cyanobacterial myoviruses (cyanophages) isolated on diverse Prochlorococcus or Synechococcus hosts. A core genome of 38 virion construction and DNA replication genes was observed in all 26 genomes, with 32 and 25 additional genes shared among the non-cyanophage and cyanophage subsets, respectively. These hierarchical cores are highly syntenic across the genomes, and sampled to saturation. The 25 cyanophage core genes include six previously described genes with putative functions (psbA, mazG, phoH, hsp20, hli03, cobS), a hypothetical protein with a potential phytanoyl-CoA dioxygenase domain, two virion structural genes, and 16 hypothetical genes. Beyond previously described cyanophage-encoded photosynthesis and phosphate stress genes, we observed core genes that may play a role in nitrogen metabolism during infection through modulation of 2-oxoglutarate. Patterns among non-core genes that may drive niche diversification revealed that phosphorus-related gene content reflects source waters rather than host strain used for isolation, and that carbon metabolism genes appear associated with putative mobile elements. As well, phages isolated on Synechococcus had higher genome-wide %G+C and often contained different gene subsets (e.g. petE, zwf, gnd, prnA, cpeT) than those isolated on Prochlorococcus. However, no clear diagnostic genes emerged to distinguish these phage groups, suggesting blurred boundaries possibly due to cross-infection. Finally, genome-wide comparisons of both diverse and closely related, co-isolated genomes provide a locus-to-locus variability metric that will prove valuable for interpreting metagenomic data sets., (© 2010 Society for Applied Microbiology and Blackwell Publishing Ltd.)

Bacterial populations harbor vast genetic diversity that is continually shaped by abiotic and biotic selective pressures, as well as by neutral processes. Individuals coexisting in the same geographically defined population often have significantly different gene content, but whether this variation is largely adaptive or neutral remains poorly understood. Here we quantify heterogeneity in gene content for two model marine microbes, Prochlorococcus and Pelagibacter, within and between populations in the Atlantic and Pacific Oceans, to begin to understand the selective pressures that are shaping these "population genomes." We discovered a large fraction of genes that are rare in each population, reflecting continual gene transfer and loss. Despite this high variation within each population, only a few genes significantly differ in abundance between the two biogeochemically distinct environments; nearly all of these are related to phosphorus acquisition and are enriched in the Atlantic relative to the Pacific. Moreover, P-related genes from the two sites form phylogenetically distinct clusters, whereas housekeeping genes do not, consistent with a recent spread of adaptive P-related genes in the Atlantic populations. These findings implicate phosphorus availability as the dominant selective force driving divergence between these populations, and demonstrate the promise of this approach for revealing selective agents in more complex microbial systems.

To better understand the temporal and spatial dynamics of Prochlorococcus populations, and how these populations co-vary with the physical environment, we followed monthly changes in the abundance of five ecotypes-two high-light adapted and three low-light adapted-over a 5-year period in coordination with the Bermuda Atlantic Time Series (BATS) and Hawaii Ocean Time-series (HOT) programs. Ecotype abundance displayed weak seasonal fluctuations at HOT and strong seasonal fluctuations at BATS. Furthermore, stable 'layered' depth distributions, where different Prochlorococcus ecotypes reached maximum abundance at different depths, were maintained consistently for 5 years at HOT. Layered distributions were also observed at BATS, although winter deep mixing events disrupted these patterns each year and produced large variations in ecotype abundance. Interestingly, the layered ecotype distributions were regularly reestablished each year after deep mixing subsided at BATS. In addition, Prochlorococcus ecotypes each responded differently to the strong seasonal changes in light, temperature and mixing at BATS, resulting in a reproducible annual succession of ecotype blooms. Patterns of ecotype abundance, in combination with physiological assays of cultured isolates, confirmed that the low-light adapted eNATL could be distinguished from other low-light adapted ecotypes based on its ability to withstand temporary exposure to high-intensity light, a characteristic stress of the surface mixed layer. Finally, total Prochlorococcus and Synechococcus dynamics were compared with similar time series data collected a decade earlier at each location. The two data sets were remarkably similar-testimony to the resilience of these complex dynamic systems on decadal time scales.

Marine dissolved organic matter (DOM) contains as much carbon as the Earth's atmosphere, and represents a critical component of the global carbon cycle. To better define microbial processes and activities associated with marine DOM cycling, we analyzed genomic and transcriptional responses of microbial communities to high-molecular-weight DOM (HMWDOM) addition. The cell density in the unamended control remained constant, with very few transcript categories exhibiting significant differences over time. In contrast, the DOM-amended microcosm doubled in cell numbers over 27 h, and a variety of HMWDOM-stimulated transcripts from different taxa were observed at all time points measured relative to the control. Transcripts significantly enriched in the HMWDOM treatment included those associated with two-component sensor systems, phosphate and nitrogen assimilation, chemotaxis, and motility. Transcripts from Idiomarina and Alteromonas spp., the most highly represented taxa at the early time points, included those encoding TonB-associated transporters, nitrogen assimilation genes, fatty acid catabolism genes, and TCA cycle enzymes. At the final time point, Methylophaga rRNA and non-rRNA transcripts dominated the HMWDOM-amended microcosm, and included gene transcripts associated with both assimilatory and dissimilatory single-carbon compound utilization. The data indicated specific resource partitioning of DOM by different bacterial species, which results in a temporal succession of taxa, metabolic pathways, and chemical transformations associated with HMWDOM turnover. These findings suggest that coordinated, cooperative activities of a variety of bacterial "specialists" may be critical in the cycling of marine DOM, emphasizing the importance of microbial community dynamics in the global carbon cycle.

Background: Different high-throughput nucleic acid sequencing platforms are currently available but a trade-off currently exists between the cost and number of reads that can be generated versus the read length that can be achieved., Methodology/principal Findings: We describe an experimental and computational pipeline yielding millions of reads that can exceed 200 bp with quality scores approaching that of traditional Sanger sequencing. The method combines an automatable gel-less library construction step with paired-end sequencing on a short-read instrument. With appropriately sized library inserts, mate-pair sequences can overlap, and we describe the SHERA software package that joins them to form a longer composite read., Conclusions/significance: This strategy is broadly applicable to sequencing applications that benefit from low-cost high-throughput sequencing, but require longer read lengths. We demonstrate that our approach enables metagenomic analyses using the Illumina Genome Analyzer, with low error rates, and at a fraction of the cost of pyrosequencing.

Exposure to solar radiation can cause mortality in natural communities of pico-phytoplankton, both at the surface and to a depth of at least 30 m. DNA damage is a significant cause of death, mainly due to cyclobutane pyrimidine dimer formation, which can be lethal if not repaired. While developing a UV mutagenesis protocol for the marine cyanobacterium Prochlorococcus, we isolated a UV-hyper-resistant variant of high light-adapted strain MED4. The hyper-resistant strain was constitutively upregulated for expression of the mutT-phrB operon, encoding nudix hydrolase and photolyase, both of which are involved in repair of DNA damage that can be caused by UV light. Photolyase (PhrB) breaks pyrimidine dimers typically caused by UV exposure, using energy from visible light in the process known as photoreactivation. Nudix hydrolase (MutT) hydrolyses 8-oxo-dGTP, an aberrant form of GTP that results from oxidizing conditions, including UV radiation, thus impeding mispairing and mutagenesis by preventing incorporation of the aberrant form into DNA. These processes are error-free, in contrast to error-prone SOS dark repair systems that are widespread in bacteria. The UV-hyper-resistant strain contained only a single mutation: a 1 bp deletion in the intergenic region directly upstream of the mutT-phrB operon. Two subsequent enrichments for MED4 UV-hyper-resistant strains from MED4 wild-type cultures gave rise to strains containing this same 1 bp deletion, affirming its connection to the hyper-resistant phenotype. These results have implications for Prochlorococcus DNA repair mechanisms, genome stability and possibly lysogeny.

Podovirus P-SSP7 infects Prochlorococcus marinus, the most abundant oceanic photosynthetic microorganism. Single-particle cryo-electron microscopy yields icosahedral and asymmetrical structures of infectious P-SSP7 with 4.6-A and 9-A resolution, respectively. The asymmetric reconstruction reveals how symmetry mismatches are accommodated among five of the gene products at the portal vertex. Reconstructions of infectious and empty particles show a conformational change of the 'valve' density in the nozzle, an orientation difference in the tail fibers, a disordering of the C terminus of the portal protein and the disappearance of the core proteins. In addition, cryo-electron tomography of P-SSP7 infecting Prochlorococcus showed the same tail-fiber conformation as that in empty particles. Our observations suggest a mechanism whereby, upon binding to the host cell, the tail fibers induce a cascade of structural alterations of the portal vertex complex that triggers DNA release.

Our understanding of secondary metabolite production in bacteria has been shaped primarily by studies of attached varieties such as symbionts, pathogens, and soil bacteria. Here we show that a strain of the single-celled, planktonic marine cyanobacterium Prochlorococcus-which conducts a sizable fraction of photosynthesis in the oceans-produces many cyclic, lanthionine-containing peptides (lantipeptides). Remarkably, in Prochlorococcus MIT9313 a single promiscuous enzyme transforms up to 29 different linear ribosomally synthesized peptides into a library of polycyclic, conformationally constrained products with highly diverse ring topologies. Genes encoding this system are found in variable abundances across the oceans-with a hot spot in a Galapagos hypersaline lagoon-suggesting they play a habitat- and/or community-specific role. The extraordinarily efficient pathway for generating structural diversity enables these cyanobacteria to produce as many secondary metabolites as model antibiotic-producing bacteria, but with much smaller genomes.

Our view of marine microbes is transforming, as culture-independent methods facilitate rapid characterization of microbial diversity. It is difficult to assimilate this information into our understanding of marine microbe ecology and evolution, because their distributions, traits, and genomes are shaped by forces that are complex and dynamic. Here we incorporate diverse forces--physical, biogeochemical, ecological, and mutational--into a global ocean model to study selective pressures on a simple trait in a widely distributed lineage of picophytoplankton: the nitrogen use abilities of Synechococcus and Prochlorococcus cyanobacteria. Some Prochlorococcus ecotypes have lost the ability to use nitrate, whereas their close relatives, marine Synechococcus, typically retain it. We impose mutations for the loss of nitrogen use abilities in modeled picophytoplankton, and ask: in which parts of the ocean are mutants most disadvantaged by losing the ability to use nitrate, and in which parts are they least disadvantaged? Our model predicts that this selective disadvantage is smallest for picophytoplankton that live in tropical regions where Prochlorococcus are abundant in the real ocean. Conversely, the selective disadvantage of losing the ability to use nitrate is larger for modeled picophytoplankton that live at higher latitudes, where Synechococcus are abundant. In regions where we expect Prochlorococcus and Synechococcus populations to cycle seasonally in the real ocean, we find that model ecotypes with seasonal population dynamics similar to Prochlorococcus are less disadvantaged by losing the ability to use nitrate than model ecotypes with seasonal population dynamics similar to Synechococcus. The model predictions for the selective advantage associated with nitrate use are broadly consistent with the distribution of this ability among marine picocyanobacteria, and at finer scales, can provide insights into interactions between temporally varying ocean processes and selective pressures that may be difficult or impossible to study by other means. More generally, and perhaps more importantly, this study introduces an approach for testing hypotheses about the processes that underlie genetic variation among marine microbes, embedded in the dynamic physical, chemical, and biological forces that generate and shape this diversity.

Background: Bacterial viruses (phages) play a critical role in shaping microbial populations as they influence both host mortality and horizontal gene transfer. As such, they have a significant impact on local and global ecosystem function and human health. Despite their importance, little is known about the genomic diversity harbored in phages, as methods to capture complete phage genomes have been hampered by the lack of knowledge about the target genomes, and difficulties in generating sufficient quantities of genomic DNA for sequencing. Of the approximately 550 phage genomes currently available in the public domain, fewer than 5% are marine phage., Methodology/principal Findings: To advance the study of phage biology through comparative genomic approaches we used marine cyanophage as a model system. We compared DNA preparation methodologies (DNA extraction directly from either phage lysates or CsCl purified phage particles), and sequencing strategies that utilize either Sanger sequencing of a linker amplification shotgun library (LASL) or of a whole genome shotgun library (WGSL), or 454 pyrosequencing methods. We demonstrate that genomic DNA sample preparation directly from a phage lysate, combined with 454 pyrosequencing, is best suited for phage genome sequencing at scale, as this method is capable of capturing complete continuous genomes with high accuracy. In addition, we describe an automated annotation informatics pipeline that delivers high-quality annotation and yields few false positives and negatives in ORF calling., Conclusions/significance: These DNA preparation, sequencing and annotation strategies enable a high-throughput approach to the burgeoning field of phage genomics.

Background: RNA turnover plays an important role in the gene regulation of microorganisms and influences their speed of acclimation to environmental changes. We investigated whole-genome RNA stability of Prochlorococcus, a relatively slow-growing marine cyanobacterium doubling approximately once a day, which is extremely abundant in the oceans., Results: Using a combination of microarrays, quantitative RT-PCR and a new fitting method for determining RNA decay rates, we found a median half-life of 2.4 minutes and a median decay rate of 2.6 minutes for expressed genes - twofold faster than that reported for any organism. The shortest transcript half-life (33 seconds) was for a gene of unknown function, while some of the longest (approximately 18 minutes) were for genes with high transcript levels. Genes organized in operons displayed intriguing mRNA decay patterns, such as increased stability, and delayed onset of decay with greater distance from the transcriptional start site. The same phenomenon was observed on a single probe resolution for genes greater than 2 kb., Conclusions: We hypothesize that the fast turnover relative to the slow generation time in Prochlorococcus may enable a swift response to environmental changes through rapid recycling of nucleotides, which could be advantageous in nutrient poor oceans. Our growing understanding of RNA half-lives will help us interpret the growing bank of metatranscriptomic studies of wild populations of Prochlorococcus. The surprisingly complex decay patterns of large transcripts reported here, and the method developed to describe them, will open new avenues for the investigation and understanding of RNA decay for all organisms.

Prochlorococcus, an abundant phototroph in the oceans, are infected by members of three families of viruses: myo-, podo- and siphoviruses. Genomes of myo- and podoviruses isolated on Prochlorococcus contain DNA replication machinery and virion structural genes homologous to those from coliphages T4 and T7 respectively. They also contain a suite of genes of cyanobacterial origin, most notably photosynthesis genes, which are expressed during infection and appear integral to the evolutionary trajectory of both host and phage. Here we present the first genome of a cyanobacterial siphovirus, P-SS2, which was isolated from Atlantic slope waters using a Prochlorococcus host (MIT9313). The P-SS2 genome is larger than, and considerably divergent from, previously sequenced siphoviruses. It appears most closely related to lambdoid siphoviruses, with which it shares 13 functional homologues. The approximately 108 kb P-SS2 genome encodes 131 predicted proteins and notably lacks photosynthesis genes which have consistently been found in other marine cyanophage, but does contain 14 other cyanobacterial homologues. While only six structural proteins were identified from the genome sequence, 35 proteins were detected experimentally; these mapped onto capsid and tail structural modules in the genome. P-SS2 is potentially capable of integration into its host as inferred from bioinformatically identified genetic machinery int, bet, exo and a 53 bp attachment site. The host attachment site appears to be a genomic island that is tied to insertion sequence (IS) activity that could facilitate mobility of a gene involved in the nitrogen-stress response. The homologous region and a secondary IS-element hot-spot in Synechococcus RS9917 are further evidence of IS-mediated genome evolution coincident with a probable relic prophage integration event. This siphovirus genome provides a glimpse into the biology of a deep-photic zone phage as well as the ocean cyanobacterial prophage and IS element 'mobilome'.

Bacterial microcompartments (BMCs) are polyhedral bodies, composed entirely of proteins, that function as organelles in bacteria; they promote subcellular processes by encapsulating and co-localizing targeted enzymes with their substrates. The best-characterized BMC is the carboxysome, a central part of the carbon-concentrating mechanism that greatly enhances carbon fixation in cyanobacteria and some chemoautotrophs. Here we report the first structural insights into the carboxysome of Prochlorococcus, the numerically dominant cyanobacterium in the world's oligotrophic oceans. Bioinformatic methods, substantiated by analysis of gene expression data, were used to identify a new carboxysome shell component, CsoS1D, in the genome of Prochlorococcus strain MED4; orthologs were subsequently found in all cyanobacteria. Two independent crystal structures of Prochlorococcus MED4 CsoS1D reveal three features not seen in any BMC-domain protein structure solved to date. First, CsoS1D is composed of a fused pair of BMC domains. Second, this double-domain protein trimerizes to form a novel pseudohexameric building block for incorporation into the carboxysome shell, and the trimers further dimerize, forming a two-tiered shell building block. Third, and most strikingly, the large pore formed at the 3-fold axis of symmetry appears to be gated. Each dimer of trimers contains one trimer with an open pore and one whose pore is obstructed due to side-chain conformations of two residues that are invariant among all CsoS1D orthologs. This is the first evidence of the potential for gated transport across the carboxysome shell and reveals a new type of building block for BMC shells.

Background: Single-cell genome sequencing has the potential to allow the in-depth exploration of the vast genetic diversity found in uncultured microbes. We used the marine cyanobacterium Prochlorococcus as a model system for addressing important challenges facing high-throughput whole genome amplification (WGA) and complete genome sequencing of individual cells., Methodology/principal Findings: We describe a pipeline that enables single-cell WGA on hundreds of cells at a time while virtually eliminating non-target DNA from the reactions. We further developed a post-amplification normalization procedure that mitigates extreme variations in sequencing coverage associated with multiple displacement amplification (MDA), and demonstrated that the procedure increased sequencing efficiency and facilitated genome assembly. We report genome recovery as high as 99.6% with reference-guided assembly, and 95% with de novo assembly starting from a single cell. We also analyzed the impact of chimera formation during MDA on de novo assembly, and discuss strategies to minimize the presence of incorrectly joined regions in contigs., Conclusions/significance: The methods describe in this paper will be useful for sequencing genomes of individual cells from a variety of samples.

In order to expand our understanding of the diversity and biogeography of Prochlorococcus ribotypes, we PCR-amplified, cloned and sequenced the 16S/23S rRNA ITS region from sites in the Atlantic and Pacific oceans. Ninety-three per cent of the ITS sequences could be assigned to existing Prochlorococcus clades, although many novel subclades were detected. We assigned the sequences to operational taxonomic units using a graduated scale of sequence identity from 80% to 99.5% and correlated Prochlorococcus diversity with respect to environmental variables and dispersal time between the sites. Dispersal time was estimated using a global ocean circulation model. The significance of specific environmental variables was dependent on the degree of sequence identity used to define a taxon: light correlates with broad-scale diversity (90% cut-off), temperature with intermediate scale (95%) whereas no correlation with phosphate was observed. Community structure was correlated with dispersal time between sample sites only when taxa were defined using the finest sequence similarity cut-off. Surprisingly, the concentration of nitrate, which cannot be used as N source by the Prochlorococcus strains in culture, explains some variation in community structure for some definitions of taxa. This study suggests that the spatial distribution of Prochlorococcus ecotypes is shaped by a hierarchy of environmental factors as well dispersal limitation.

Prochlorococcus and Synechococcus are the two most abundant marine cyanobacteria. They represent a significant fraction of the total primary production of the world oceans and comprise a major fraction of the prey biomass available to phagotrophic protists. Despite relatively rapid growth rates, picocyanobacterial cell densities in open-ocean surface waters remain fairly constant, implying steady mortality due to viral infection and consumption by predators. There have been several studies on grazing by specific protists on Prochlorococcus and Synechococcus in culture, and of cell loss rates due to overall grazing in the field. However, the specific sources of mortality of these primary producers in the wild remain unknown. Here, we use a modification of the RNA stable isotope probing technique (RNA-SIP), which involves adding labelled cells to natural seawater, to identify active predators that are specifically consuming Prochlorococcus and Synechococcus in the surface waters of the Pacific Ocean. Four major groups were identified as having their 18S rRNA highly labelled: Prymnesiophyceae (Haptophyta), Dictyochophyceae (Stramenopiles), Bolidomonas (Stramenopiles) and Dinoflagellata (Alveolata). For the first three of these, the closest relative of the sequences identified was a photosynthetic organism, indicating the presence of mixotrophs among picocyanobacterial predators. We conclude that the use of RNA-SIP is a useful method to identity specific predators for picocyanobacteria in situ, and that the method could possibly be used to identify other bacterial predators important in the microbial food-web.

The marine cyanobacterium Prochlorococcus MED4 has the smallest genome and cell size of all known photosynthetic organisms. Like all phototrophs at temperate latitudes, it experiences predictable daily variation in available light energy which leads to temporal regulation and partitioning of key cellular processes. To better understand the tempo and choreography of this minimal phototroph, we studied the entire transcriptome of the cell over a simulated daily light-dark cycle, and placed it in the context of diagnostic physiological and cell cycle parameters. All cells in the culture progressed through their cell cycles in synchrony, thus ensuring that our measurements reflected the behavior of individual cells. Ninety percent of the annotated genes were expressed, and 80% had cyclic expression over the diel cycle. For most genes, expression peaked near sunrise or sunset, although more subtle phasing of gene expression was also evident. Periodicities of the transcripts of genes involved in physiological processes such as in cell cycle progression, photosynthesis, and phosphorus metabolism tracked the timing of these activities relative to the light-dark cycle. Furthermore, the transitions between photosynthesis during the day and catabolic consumption of energy reserves at night- metabolic processes that share some of the same enzymes--appear to be tightly choreographed at the level of RNA expression. In-depth investigation of these patterns identified potential regulatory proteins involved in balancing these opposing pathways. Finally, while this analysis has not helped resolve how a cell with so little regulatory capacity, and a 'deficient' circadian mechanism, aligns its cell cycle and metabolism so tightly to a light-dark cycle, it does provide us with a valuable framework upon which to build when the Prochlorococcus proteome and metabolome become available.

Oceanic phages are critical components of the global ecosystem, where they play a role in microbial mortality and evolution. Our understanding of phage diversity is greatly limited by the lack of useful genetic diversity measures. Previous studies, focusing on myophages that infect the marine cyanobacterium Synechococcus, have used the coliphage T4 portal-protein-encoding homologue, gene 20 (g20), as a diversity marker. These studies revealed 10 sequence clusters, 9 oceanic and 1 freshwater, where only 3 contained cultured representatives. We sequenced g20 from 38 marine myophages isolated using a diversity of Synechococcus and Prochlorococcus hosts to see if any would fall into the clusters that lacked cultured representatives. On the contrary, all fell into the three clusters that already contained sequences from cultured phages. Further, there was no obvious relationship between host of isolation, or host range, and g20 sequence similarity. We next expanded our analyses to all available g20 sequences (769 sequences), which include PCR amplicons from wild uncultured phages, non-PCR amplified sequences identified in the Global Ocean Survey (GOS) metagenomic database, as well as sequences from cultured phages, to evaluate the relationship between g20 sequence clusters and habitat features from which the phage sequences were isolated. Even in this meta-data set, very few sequences fell into the sequence clusters without cultured representatives, suggesting that the latter are very rare, or sequencing artefacts. In contrast, sequences most similar to the culture-containing clusters, the freshwater cluster and two novel clusters, were more highly represented, with one particular culture-containing cluster representing the dominant g20 genotype in the unamplified GOS sequence data. Finally, while some g20 sequences were non-randomly distributed with respect to habitat, there were always numerous exceptions to general patterns, indicating that phage portal proteins are not good predictors of a phage's host or the habitat in which a particular phage may thrive.

Prochlorococcus, an extremely small cyanobacterium that is very abundant in the world's oceans, has a very streamlined genome. On average, these cells have about 2,000 genes and very few regulatory proteins. The limited capability of regulation is thought to be a result of selection imposed by a relatively stable environment in combination with a very small genome. Furthermore, only ten non-coding RNAs (ncRNAs), which play crucial regulatory roles in all forms of life, have been described in Prochlorococcus. Most strains also lack the RNA chaperone Hfq, raising the question of how important this mode of regulation is for these cells. To explore this question, we examined the transcription of intergenic regions of Prochlorococcus MED4 cells subjected to a number of different stress conditions: changes in light qualities and quantities, phage infection, or phosphorus starvation. Analysis of Affymetrix microarray expression data from intergenic regions revealed 276 novel transcriptional units. Among these were 12 new ncRNAs, 24 antisense RNAs (asRNAs), as well as 113 short mRNAs. Two additional ncRNAs were identified by homology, and all 14 new ncRNAs were independently verified by Northern hybridization and 5'RACE. Unlike its reduced suite of regulatory proteins, the number of ncRNAs relative to genome size in Prochlorococcus is comparable to that found in other bacteria, suggesting that RNA regulators likely play a major role in regulation in this group. Moreover, the ncRNAs are concentrated in previously identified genomic islands, which carry genes of significance to the ecology of this organism, many of which are not of cyanobacterial origin. Expression profiles of some of these ncRNAs suggest involvement in light stress adaptation and/or the response to phage infection consistent with their location in the hypervariable genomic islands., Competing Interests: The authors have declared that no competing interests exist.

Although the oceanic cyanobacterium Prochlorococcus harvests light with a chlorophyll antenna [1-3] rather than with the phycobilisomes that are typical of cyanobacteria, some strains express genes that are remnants of the ancestral Synechococcus phycobilisomes [4]. Similarly, some Prochlorococcus cyanophages, which often harbor photosynthesis-related genes [5], also carry homologs of phycobilisome pigment biosynthesis genes [6, 7]. Here, we investigate four such genes in two cyanophages that both infect abundant Prochlorococcus strains [8]: homologs of heme oxygenase (ho1), 15,16-dihydrobiliverdin:ferredoxin oxidoreductase (pebA), ferredoxin (petF) in the myovirus P-SSM2, and a phycocyanobilin:ferredoxin oxidoreductase (pcyA) homolog in the myovirus P-SSM4. We demonstrate that the phage homologs mimic the respective host activities, with the exception of the divergent phage PebA homolog. In this case, the phage PebA single-handedly catalyzes a reaction for which uninfected host cells require two consecutive enzymes, PebA and PebB. We thus renamed the phage enzyme phycoerythrobilin synthase (PebS). This gene, and other pigment biosynthesis genes encoded by P-SSM2 (petF and ho1), are transcribed during infection, suggesting that they can improve phage fitness. Analyses of global ocean metagenomes show that PcyA and Ho1 occur in both cyanobacteria and their phages, whereas the novel PebS-encoding gene is exclusive to phages.

Metagenomics is expanding our knowledge of the gene content, functional significance, and genetic variability in natural microbial communities. Still, there exists limited information concerning the regulation and dynamics of genes in the environment. We report here global analysis of expressed genes in a naturally occurring microbial community. We first adapted RNA amplification technologies to produce large amounts of cDNA from small quantities of total microbial community RNA. The fidelity of the RNA amplification procedure was validated with Prochlorococcus cultures and then applied to a microbial assemblage collected in the oligotrophic Pacific Ocean. Microbial community cDNAs were analyzed by pyrosequencing and compared with microbial community genomic DNA sequences determined from the same sample. Pyrosequencing-based estimates of microbial community gene expression compared favorably to independent assessments of individual gene expression using quantitative PCR. Genes associated with key metabolic pathways in open ocean microbial species-including genes involved in photosynthesis, carbon fixation, and nitrogen acquisition-and a number of genes encoding hypothetical proteins were highly represented in the cDNA pool. Genes present in the variable regions of Prochlorococcus genomes were among the most highly expressed, suggesting these encode proteins central to cellular processes in specific genotypes. Although many transcripts detected were highly similar to genes previously detected in ocean metagenomic surveys, a significant fraction ( approximately 50%) were unique. Thus, microbial community transcriptomic analyses revealed not only indigenous gene- and taxon-specific expression patterns but also gene categories undetected in previous DNA-based metagenomic surveys.

Background: Phages infecting marine picocyanobacteria often carry a psbA gene, which encodes a homolog to the photosynthetic reaction center protein, D1. Host encoded D1 decays during phage infection in the light. Phage encoded D1 may help to maintain photosynthesis during the lytic cycle, which in turn could bolster the production of deoxynucleoside triphosphates (dNTPs) for phage genome replication., Methodology/principal Findings: To explore the consequences to a phage of encoding and expressing psbA, we derive a simple model of infection for a cyanophage/host pair--cyanophage P-SSP7 and Prochlorococcus MED4--for which pertinent laboratory data are available. We first use the model to describe phage genome replication and the kinetics of psbA expression by host and phage. We then examine the contribution of phage psbA expression to phage genome replication under constant low irradiance (25 microE m(-2) s(-1)). We predict that while phage psbA expression could lead to an increase in the number of phage genomes produced during a lytic cycle of between 2.5 and 4.5% (depending on parameter values), this advantage can be nearly negated by the cost of psbA in elongating the phage genome. Under higher irradiance conditions that promote D1 degradation, however, phage psbA confers a greater advantage to phage genome replication., Conclusions/significance: These analyses illustrate how psbA may benefit phage in the dynamic ocean surface mixed layer.

Prochlorococcus is a marine cyanobacterium that numerically dominates the mid-latitude oceans and is the smallest known oxygenic phototroph. Numerous isolates from diverse areas of the world's oceans have been studied and shown to be physiologically and genetically distinct. All isolates described thus far can be assigned to either a tightly clustered high-light (HL)-adapted clade, or a more divergent low-light (LL)-adapted group. The 16S rRNA sequences of the entire Prochlorococcus group differ by at most 3%, and the four initially published genomes revealed patterns of genetic differentiation that help explain physiological differences among the isolates. Here we describe the genomes of eight newly sequenced isolates and combine them with the first four genomes for a comprehensive analysis of the core (shared by all isolates) and flexible genes of the Prochlorococcus group, and the patterns of loss and gain of the flexible genes over the course of evolution. There are 1,273 genes that represent the core shared by all 12 genomes. They are apparently sufficient, according to metabolic reconstruction, to encode a functional cell. We describe a phylogeny for all 12 isolates by subjecting their complete proteomes to three different phylogenetic analyses. For each non-core gene, we used a maximum parsimony method to estimate which ancestor likely first acquired or lost each gene. Many of the genetic differences among isolates, especially for genes involved in outer membrane synthesis and nutrient transport, are found within the same clade. Nevertheless, we identified some genes defining HL and LL ecotypes, and clades within these broad ecotypes, helping to demonstrate the basis of HL and LL adaptations in Prochlorococcus. Furthermore, our estimates of gene gain events allow us to identify highly variable genomic islands that are not apparent through simple pairwise comparisons. These results emphasize the functional roles, especially those connected to outer membrane synthesis and transport that dominate the flexible genome and set it apart from the core. Besides identifying islands and demonstrating their role throughout the history of Prochlorococcus, reconstruction of past gene gains and losses shows that much of the variability exists at the "leaves of the tree," between the most closely related strains. Finally, the identification of core and flexible genes from this 12-genome comparison is largely consistent with the relative frequency of Prochlorococcus genes found in global ocean metagenomic databases, further closing the gap between our understanding of these organisms in the lab and the wild., Competing Interests: Competing interests. The authors have declared that no competing interests exist.

Interactions between bacterial hosts and their viruses (phages) lead to reciprocal genome evolution through a dynamic co-evolutionary process. Phage-mediated transfer of host genes--often located in genome islands--has had a major impact on microbial evolution. Furthermore, phage genomes have clearly been shaped by the acquisition of genes from their hosts. Here we investigate whole-genome expression of a host and phage, the marine cyanobacterium Prochlorococcus MED4 and the T7-like cyanophage P-SSP7, during lytic infection, to gain insight into these co-evolutionary processes. Although most of the phage genome was linearly transcribed over the course of infection, four phage-encoded bacterial metabolism genes formed part of the same expression cluster, even though they are physically separated on the genome. These genes--encoding photosystem II D1 (psbA), high-light inducible protein (hli), transaldolase (talC) and ribonucleotide reductase (nrd)--are transcribed together with phage DNA replication genes and seem to make up a functional unit involved in energy and deoxynucleotide production for phage replication in resource-poor oceans. Also unique to this system was the upregulation of numerous genes in the host during infection. These may be host stress response genes and/or genes induced by the phage. Many of these host genes are located in genome islands and have homologues in cyanophage genomes. We hypothesize that phage have evolved to use upregulated host genes, leading to their stable incorporation into phage genomes and their subsequent transfer back to hosts in genome islands. Thus activation of host genes during infection may be directing the co-evolution of gene content in both host and phage genomes.

Prochlorococcus is a simple cyanobacterium that is abundant throughout large regions of the oceans, and has become a useful model for studying the nature and regulation of biological diversity across all scales of complexity. Recent work has revealed that environmental factors such as light, nutrients and predation influence diversity in different ways, changing our image of the structure and dynamics of the global Prochlorococcus population. Advances in metagenomics, transcription profiling and global ecosystem modeling promise to deliver an even greater understanding of this system and further demonstrate the power of cross-scale systems biology.

A marine ecosystem model seeded with many phytoplankton types, whose physiological traits were randomly assigned from ranges defined by field and laboratory data, generated an emergent community structure and biogeography consistent with observed global phytoplankton distributions. The modeled organisms included types analogous to the marine cyanobacterium Prochlorococcus. Their emergent global distributions and physiological properties simultaneously correspond to observations. This flexible representation of community structure can be used to explore relations between ecosystems, biogeochemical cycles, and climate change.

Prochlorococcus MED4 has, with a total of only 1,716 annotated protein-coding genes, the most compact genome of a free-living photoautotroph. Although light quality and quantity play an important role in regulating the growth rate of this organism in its natural habitat, the majority of known light-sensing proteins are absent from its genome. To explore the potential for light sensing in this phototroph, we measured its global gene expression pattern in response to different light qualities and quantities by using high-density Affymetrix microarrays. Though seven different conditions were tested, only blue light elicited a strong response. In addition, hierarchical clustering revealed that the responses to high white light and blue light were very similar and different from that of the lower-intensity white light, suggesting that the actual sensing of high light is mediated via a blue-light receptor. Bacterial cryptochromes seem to be good candidates for the blue-light sensors. The existence of a signaling pathway for the redox state of the photosynthetic electron transport chain was suggested by the presence of genes that responded similarly to red and blue light as well as genes that responded to the addition of DCMU [3-(3,4-dichlorophenyl)-1,1-N-N'-dimethylurea], a specific inhibitor of photosystem II-mediated electron transport.

The cyanobacterium Prochlorococcus is the numerically dominant phototroph in the oligotrophic oceans. This group consists of multiple ecotypes that are physiologically and phylogenetically distinct and occur in different abundances along environmental gradients. Here we examine adaptations to phosphate (P) limitation among ecotypes. First, we used DNA microarrays to identify genes involved in the P-starvation response in two strains belonging to different ecotypes, MED4 (high-light-adapted) and MIT9313 (low-light-adapted). Most of the up-regulated genes under P starvation were unique to one strain. In MIT9313, many ribosomal genes were down-regulated, suggesting a general stress response in this strain. We also observed major differences in regulation. The P-starvation-induced genes comprise two clusters on the chromosome, the first containing the P master regulator phoB and most known P-acquisition genes and the second, absent in MIT9313, containing genes of unknown function. We examined the organization of the phoB gene cluster in 11 Prochlorococcus strains belonging to diverse ecotypes and found high variability in gene content that was not congruent with rRNA phylogeny. We hypothesize that this genome variability is related to differences in P availability in the oceans from which the strains were isolated. Analysis of a metagenomic library from the Sargasso Sea supports this hypothesis; most Prochlorococcus cells in this low-P environment contain the P-acquisition genes seen in MED4, although a number of previously undescribed gene combinations were observed.

Cyanophages (cyanobacterial viruses) are important agents of horizontal gene transfer among marine cyanobacteria, the numerically dominant photosynthetic organisms in the oceans. Some cyanophage genomes carry and express host-like photosynthesis genes, presumably to augment the host photosynthetic machinery during infection. To study the prevalence and evolutionary dynamics of this phenomenon, 33 cultured cyanophages of known family and host range and viral DNA from field samples were screened for the presence of two core photosystem reaction center genes, psbA and psbD. Combining this expanded dataset with published data for nine other cyanophages, we found that 88% of the phage genomes contain psbA, and 50% contain both psbA and psbD. The psbA gene was found in all myoviruses and Prochlorococcus podoviruses, but could not be amplified from Prochlorococcus siphoviruses or Synechococcus podoviruses. Nearly all of the phages that encoded both psbA and psbD had broad host ranges. We speculate that the presence or absence of psbA in a phage genome may be determined by the length of the latent period of infection. Whether it also carries psbD may reflect constraints on coupling of viral- and host-encoded PsbA-PsbD in the photosynthetic reaction center across divergent hosts. Phylogenetic clustering patterns of these genes from cultured phages suggest that whole genes have been transferred from host to phage in a discrete number of events over the course of evolution (four for psbA, and two for psbD), followed by horizontal and vertical transfer between cyanophages. Clustering patterns of psbA and psbD from Synechococcus cells were inconsistent with other molecular phylogenetic markers, suggesting genetic exchanges involving Synechococcus lineages. Signatures of intragenic recombination, detected within the cyanophage gene pool as well as between hosts and phages in both directions, support this hypothesis. The analysis of cyanophage psbA and psbD genes from field populations revealed significant sequence diversity, much of which is represented in our cultured isolates. Collectively, these findings show that photosynthesis genes are common in cyanophages and that significant genetic exchanges occur from host to phage, phage to host, and within the phage gene pool. This generates genetic diversity among the phage, which serves as a reservoir for their hosts, and in turn influences photosystem evolution.

Genome sequencing currently requires DNA from pools of numerous nearly identical cells (clones), leaving the genome sequences of many difficult-to-culture microorganisms unattainable. We report a sequencing strategy that eliminates culturing of microorganisms by using real-time isothermal amplification to form polymerase clones (plones) from the DNA of single cells. Two Escherichia coli plones, analyzed by Affymetrix chip hybridization, demonstrate that plonal amplification is specific and the bias is randomly distributed. Whole-genome shotgun sequencing of Prochlorococcus MIT9312 plones showed 62% coverage of the genome from one plone at a sequencing depth of 3.5x, and 66% coverage from a second plone at a depth of 4.7x. Genomic regions not revealed in the initial round of sequencing are recovered by sequencing PCR amplicons derived from plonal DNA. The mutation rate in single-cell amplification is <2 x 10(5), better than that of current genome sequencing standards. Polymerase cloning should provide a critical tool for systematic characterization of genome diversity in the biosphere.

Prochlorococcus is the numerically dominant phytoplankter in the oligotrophic oceans, accounting for up to half of the photosynthetic biomass and production in some regions. Here, we describe how the abundance of six known ecotypes, which have small subunit ribosomal RNA sequences that differ by less than 3%, changed along local and basin-wide environmental gradients in the Atlantic Ocean. Temperature was significantly correlated with shifts in ecotype abundance, and laboratory experiments confirmed different temperature optima and tolerance ranges for cultured strains. Light, nutrients, and competitor abundances also appeared to play a role in shaping different distributions.

Prochlorococcus ecotypes are a useful system for exploring the origin and function of diversity among closely related microbes. The genetic variability between phenotypically distinct strains that differ by less that 1% in 16S ribosomal RNA sequences occurs mostly in genomic islands. Island genes appear to have been acquired in part by phage-mediated lateral gene transfer, and some are differentially expressed under light and nutrient stress. Furthermore, genome fragments directly recovered from ocean ecosystems indicate that these islands are variable among cooccurring Prochlorococcus cells. Genomic islands in this free-living photoautotroph share features with pathogenicity islands of parasitic bacteria, suggesting a general mechanism for niche differentiation in microbial species.

Prochlorococcus is a marine cyanobacterium which is found at high abundances in world's tropical and subtropical oligotrophic oceans. The genus Prochlorococcus can be divided into two major groups based on light physiology. Both of these groups can be further subdivided into genetically distinct lineages, or ecotypes. Real-time polymerase chain reaction (PCR) assays based on sequence differences in the 16S-23S rDNA internal transcribed spacer or the 23S rDNA were developed to examine the distribution of each ecotype in the field. The real-time PCR assays enabled linear quantification of concentrations ranging from 10 to 4 x 10(5) cells ml(-1). These assays were applied to a stratified water column in the Sargasso Sea. The majority of Prochlorococcus cells above 110 m belonged to the one of the low chlorophyll b/a ratio (high-light adapted) ecotypes, while two types of high chlorophyll b/a ratio (low-light adapted) cells dominated below 110 m. The other three types were found at significantly lower numbers or not detected at all. Differences in the abundance of ecotypes within the major light physiology groupings suggest that other factors, such as nutrient utilization and differential mortality, are driving their relative distributions. Real-time PCR assays will enable further exploration of these factors and temporal and geographic variability in ecotype abundance.

Microbial life predominates in the ocean, yet little is known about its genomic variability, especially along the depth continuum. We report here genomic analyses of planktonic microbial communities in the North Pacific Subtropical Gyre, from the ocean's surface to near-sea floor depths. Sequence variation in microbial community genes reflected vertical zonation of taxonomic groups, functional gene repertoires, and metabolic potential. The distributional patterns of microbial genes suggested depth-variable community trends in carbon and energy metabolism, attachment and motility, gene mobility, and host-viral interactions. Comparative genomic analyses of stratified microbial communities have the potential to provide significant insight into higher-order community organization and dynamics.

The cyanobacterium Prochlorococcus numerically dominates the photosynthetic community in the tropical and subtropical regions of the world's oceans. Six evolutionary lineages of Prochlorococcus have been described, and their distinctive physiologies and genomes indicate that these lineages are "ecotypes" and should have different oceanic distributions. Two methods recently developed to quantify these ecotypes in the field, probe hybridization and quantitative PCR (QPCR), have shown that this is indeed the case. To facilitate a global investigation of these ecotypes, we modified our QPCR protocol to significantly increase its speed, sensitivity, and accessibility and validated the method in the western and eastern North Atlantic Ocean. We showed that all six ecotypes had distinct distributions that varied with depth and location, and, with the exception of the deeper waters at the western North Atlantic site, the total Prochlorococcus counts determined by QPCR matched the total counts measured by flow cytometry. Clone library analyses of the deeper western North Atlantic waters revealed ecotypes that are not represented in the culture collections with which the QPCR primers were designed, explaining this discrepancy. Finally, similar patterns of relative ecotype abundance were obtained in QPCR and probe hybridization analyses of the same field samples, which could allow comparisons between studies.

Nitrogen (N) often limits biological productivity in the oceanic gyres where Prochlorococcus is the most abundant photosynthetic organism. The Prochlorococcus community is composed of strains, such as MED4 and MIT9313, that have different N utilization capabilities and that belong to ecotypes with different depth distributions. An interstrain comparison of how Prochlorococcus responds to changes in ambient nitrogen is thus central to understanding its ecology. We quantified changes in MED4 and MIT9313 global mRNA expression, chlorophyll fluorescence, and photosystem II photochemical efficiency (Fv/Fm) along a time series of increasing N starvation. In addition, the global expression of both strains growing in ammonium-replete medium was compared to expression during growth on alternative N sources. There were interstrain similarities in N regulation such as the activation of a putative NtcA regulon during N stress. There were also important differences between the strains such as in the expression patterns of carbon metabolism genes, suggesting that the two strains integrate N and C metabolism in fundamentally different ways.

Cyanobacteria, and the viruses (phages) that infect them, are significant contributors to the oceanic 'gene pool'. This pool is dynamic, and the transfer of genetic material between hosts and their phages probably influences the genetic and functional diversity of both. For example, photosynthesis genes of cyanobacterial origin have been found in phages that infect Prochlorococcus and Synechococcus, the numerically dominant phototrophs in ocean ecosystems. These genes include psbA, which encodes the photosystem II core reaction centre protein D1, and high-light-inducible (hli) genes. Here we show that phage psbA and hli genes are expressed during infection of Prochlorococcus and are co-transcribed with essential phage capsid genes, and that the amount of phage D1 protein increases steadily over the infective period. We also show that the expression of host photosynthesis genes declines over the course of infection and that replication of the phage genome is a function of photosynthesis. We thus propose that the phage genes are functional in photosynthesis and that they may be increasing phage fitness by supplementing the host production of these proteins.

The oceanic cyanobacteria Prochlorococcus are globally important, ecologically diverse primary producers. It is thought that their viruses (phages) mediate population sizes and affect the evolutionary trajectories of their hosts. Here we present an analysis of genomes from three Prochlorococcus phages: a podovirus and two myoviruses. The morphology, overall genome features, and gene content of these phages suggest that they are quite similar to T7-like (P-SSP7) and T4-like (P-SSM2 and P-SSM4) phages. Using the existing phage taxonomic framework as a guideline, we examined genome sequences to establish "core" genes for each phage group. We found the podovirus contained 15 of 26 core T7-like genes and the two myoviruses contained 43 and 42 of 75 core T4-like genes. In addition to these core genes, each genome contains a significant number of "cyanobacterial" genes, i.e., genes with significant best BLAST hits to genes found in cyanobacteria. Some of these, we speculate, represent "signature" cyanophage genes. For example, all three phage genomes contain photosynthetic genes (psbA, hliP) that are thought to help maintain host photosynthetic activity during infection, as well as an aldolase family gene (talC) that could facilitate alternative routes of carbon metabolism during infection. The podovirus genome also contains an integrase gene (int) and other features that suggest it is capable of integrating into its host. If indeed it is, this would be unprecedented among cultured T7-like phages or marine cyanophages and would have significant evolutionary and ecological implications for phage and host. Further, both myoviruses contain phosphate-inducible genes (phoH and pstS) that are likely to be important for phage and host responses to phosphate stress, a commonly limiting nutrient in marine systems. Thus, these marine cyanophages appear to be variations of two well-known phages-T7 and T4-but contain genes that, if functional, reflect adaptations for infection of photosynthetic hosts in low-nutrient oceanic environments.