Matthew B. Sullivan

Genome divergence in two Prochlorococcus ecotypes reflects oceanic niche differentiation

The marine unicellular cyanobacterium Prochlorococcus is the smallest-known oxygen-evolving autotroph. It numerically dominates the phytoplankton in the tropical and subtropical oceans, and is responsible for a significant fraction of global photosynthesis. Here we compare the genomes of two Prochlorococcus strains that span the largest evolutionary distance within the Prochlorococcus lineage and that have different minimum, maximum and optimal light intensities for growth. The high-light-adapted ecotype has the smallest genome (1,657,990 base pairs, 1,716 genes) of any known oxygenic phototroph, whereas the genome of its low-light-adapted counterpart is significantly larger, at 2,410,873 base pairs (2,275 genes). The comparative architectures of these two strains reveal dynamic genomes that are constantly changing in response to myriad selection pressures. Although the two strains have 1,350 genes in common, a significant number are not shared, and these have been differentially retained from the common ancestor, or acquired through duplication or lateral transfer. Some of these genes have obvious roles in determining the relative fitness of the ecotypes in response to key environmental variables, and hence in regulating their distribution and abundance in the oceans.

Structure and function of the global ocean microbiome

Microbes are dominant drivers of biogeochemical processes, yet drawing a global picture of functional diversity, microbial community structure, and their ecological determinants remains a grand challenge. We analyzed 7.2 terabases of metagenomic data from 243 Tara Oceans samples from 68 locations in epipelagic and mesopelagic waters across the globe to generate an ocean microbial reference gene catalog with >40 million nonredundant, mostly novel sequences from viruses, prokaryotes, and picoeukaryotes. Using 139 prokaryote-enriched samples, containing >35,000 species, we show vertical stratification with epipelagic community composition mostly driven by temperature rather than other environmental factors or geography. We identify ocean microbial core functionality and reveal that >73% of its abundance is shared with the human gut microbiome despite the physicochemical differences between these two ecosystems.

Genomic islands and the ecology and evolution of Prochlorococcus

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.

Cyanophages infecting the oceanic cyanobacterium Prochlorococcus.

Prochlorococcus is the numerically dominant phototroph in the tropical and subtropical oceans, accounting for half of the photosynthetic biomass in some areas. Here we report the isolation of cyanophages that infect Prochlorococcus, and show that although some are host-strain-specific, others cross-infect with closely related marine Synechococcus as well as between high-light- and low-light-adapted Prochlorococcus isolates, suggesting a mechanism for horizontal gene transfer. High-light-adapted Prochlorococcus hosts yielded Podoviridae exclusively, which were extremely host-specific, whereas low-light-adapted Prochlorococcus and all strains of Synechococcus yielded primarily Myoviridae, which has a broad host range. Finally, both Prochlorococcus and Synechococcus strain-specific cyanophage titres were low (< 103 ml-1) in stratified oligotrophic waters even where total cyanobacterial abundances were high (> 105 cells ml-1). These low titres in areas of high total host cell abundance seem to be a feature of open ocean ecosystems. We hypothesize that gradients in cyanobacterial population diversity, growth rates, and/or the incidence of lysogeny underlie these trends.

Determinants of community structure in the global plankton interactome

Species interaction networks are shaped by abiotic and biotic factors. Here, as part of the Tara Oceans project, we studied the photic zone interactome using environmental factors and organismal abundance profiles and found that environmental factors are incomplete predictors of community structure. We found associations across plankton functional types and phylogenetic groups to be nonrandomly distributed on the network and driven by both local and global patterns. We identified interactions among grazers, primary producers, viruses, and (mainly parasitic) symbionts and validated network-generated hypotheses using microscopy to confirm symbiotic relationships. We have thus provided a resource to support further research on ocean food webs and integrating biological components into ocean models.

Genomic analysis of oceanic cyanobacterial myoviruses compared with T4-like myoviruses from diverse hosts and environments

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.

Patterns and ecological drivers of ocean viral communities

Viruses influence ecosystems by modulating microbial population size, diversity, metabolic outputs, and gene flow. Here, we use quantitative double-stranded DNA (dsDNA) viral-fraction metagenomes (viromes) and whole viral community morphological data sets from 43 Tara Oceans expedition samples to assess viral community patterns and structure in the upper ocean. Protein cluster cataloging defined pelagic upper-ocean viral community pan and core gene sets and suggested that this sequence space is well-sampled. Analyses of viral protein clusters, populations, and morphology revealed biogeographic patterns whereby viral communities were passively transported on oceanic currents and locally structured by environmental conditions that affect host community structure. Together, these investigations establish a global ocean dsDNA viromic data set with analyses supporting the seed-bank hypothesis to explain how oceanic viral communities maintain high local diversity.

The Pacific Ocean Virome (POV): A Marine Viral Metagenomic Dataset and Associated Protein Clusters for Quantitative Viral Ecology

Bacteria and their viruses (phage) are fundamental drivers of many ecosystem processes including global biogeochemistry and horizontal gene transfer. While databases and resources for studying function in uncultured bacterial communities are relatively advanced, many fewer exist for their viral counterparts. The issue is largely technical in that the majority (often 90%) of viral sequences are functionally ‘unknown’ making viruses a virtually untapped resource of functional and physiological information. Here, we provide a community resource that organizes this unknown sequence space into 27 K high confidence protein clusters using 32 viral metagenomes from four biogeographic regions in the Pacific Ocean that vary by season, depth, and proximity to land, and include some of the first deep pelagic ocean viral metagenomes. These protein clusters more than double currently available viral protein clusters, including those from environmental datasets. Further, a protein cluster guided analysis of functional diversity revealed that richness decreased (i) from deep to surface waters, (ii) from winter to summer, (iii) and with distance from shore in surface waters only. These data provide a framework from which to draw on for future metadata-enabled functional inquiries of the vast viral unknown.

Abundant SAR11 viruses in the ocean

Several reports proposed that the extraordinary dominance of the SAR11 bacterial clade in ocean ecosystems could be a consequence of unusual mechanisms of resistance to bacteriophage infection, including ‘cryptic escape’ through reduced cell size and/or K-strategist defence specialism. Alternatively, the evolution of high surface-to-volume ratios coupled with minimal genomes containing high-affinity transporters enables unusually efficient metabolism for oxidizing dissolved organic matter in the world’s oceans that could support vast population sizes despite phage susceptibility. These ideas are important for understanding plankton ecology because they emphasize the potentially important role of top-down mechanisms in predation, thus determining the size of SAR11 populations and their concomitant role in biogeochemical cycling. Here we report the isolation of diverse SAR11 viruses belonging to two virus families in culture, for which we propose the name ‘pelagiphage’, after their host. Notably, the pelagiphage genomes were highly represented in marine viral metagenomes, demonstrating their importance in nature. One of the new phages, HTVC010P, represents a new podovirus subfamily more abundant than any seen previously, in all data sets tested, and may represent one of the most abundant virus subfamilies in the biosphere. This discovery disproves the theory that SAR11 cells are immune to viral predation and is consistent with the interpretation that the success of this highly abundant microbial clade is the result of successfully evolved adaptation to resource competition.

Phage-bacteria infection networks.

Phage and their bacterial hosts are the most abundant and genetically diverse group of organisms on the planet. Given their dominance, it is no wonder that many recent studies have found that phage-bacteria interactions strongly influence global biogeochemical cycles, incidence of human diseases, productivity of industrial microbial commodities, and patterns of microbial genome diversity. Unfortunately, given the extreme diversity and complexity of microbial communities, traditional analyses fail to characterize interaction patterns and underlying processes. Here, we review emerging systems approaches that combine empirical data with rigorous theoretical analysis to study phage-bacterial interactions as networks rather than as coupled interactions in isolation.