Search

Phytoplankton and associated microbial communities provide organic carbon to oceanic food webs and drive ecosystem dynamics. However, capturing those dynamics is challenging. Here, an in situ, semi-Lagrangian, robotic sampler profiled pelagic microbes at 4 h intervals over ~2.6 days in North Pacific high-nutrient, low-chlorophyll waters. We report on the community structure and transcriptional dynamics of microbes in an operationally large size class (>5 μm) predominantly populated by dinoflagellates, ciliates, haptophytes, pelagophytes, diatoms, cyanobacteria (chiefly Synechococcus), prasinophytes (chiefly Ostreococcus), fungi, archaea, and proteobacteria. Apart from fungi and archaea, all groups exhibited 24-h periodicity in some transcripts, but larger portions of the transcriptome oscillated in phototrophs. Periodic photosynthesis-related transcripts exhibited a temporal cascade across the morning hours, conserved across diverse phototrophic lineages. Pronounced silica:nitrate drawdown, a high flavodoxin to ferredoxin transcript ratio, and elevated expression of other Fe-stress markers indicated Fe-limitation. Fe-stress markers peaked during a photoperiodically adaptive time window that could modulate phytoplankton response to seasonal Fe-limitation. Remarkably, we observed viruses that infect the majority of abundant taxa, often with total transcriptional activity synchronized with putative hosts. Taken together, these data reveal a microbial plankton community that is shaped by recycled production and tightly controlled by Fe-limitation and viral activity.

Oil spills are a frequent perturbation to the marine environment that has rapid and significant impacts on the local microbiome. Previous studies have shown that exposure to synthetic dispersant alone did not enhance heterotrophic microbial activity or oxidation rates of specific hydrocarbon components but increased the abundance of some taxa (e.g., Colwellia ). In contrast, exposure to oil, but not dispersants, increased the abundance of other taxa (e.g., Marinobacter ) and stimulated hydrocarbon oxidation rates. Here, we advance these findings by interpreting metatranscriptomic data from this experiment to explore how and why specific components of the microbial community responded to distinct organic carbon exposure regimes. Dispersant alone was selected for a unique community and for dominant organisms that reflected treatment- and time-dependent responses. Dispersant amendment also led to diverging functional profiles among the different treatments. Similarly, oil alone was selected for a community that was distinct from treatments amended with dispersants. The presence of oil and dispersants with added nutrients led to substantial differences in microbial responses, likely suggesting increased fitness driven by the presence of additional inorganic nutrients. The oil-only additions led to a marked increase in the expression of phages, prophages, transposable elements, and plasmids (PPTEPs), suggesting that aspects of microbial community response to oil are driven by the "mobilome," potentially through viral-associated regulation of metabolic pathways in ciliates and flagellates that would otherwise throttle the microbial community through grazing.IMPORTANCEMicrocosm experiments simulated the April 2010 Deepwater Horizon oil spill by applying oil and synthetic dispersants (Corexit EC9500A and EC9527A) to deep ocean water samples. The exposure regime revealed severe negative alterations in the treatments' heterotrophic microbial activity and hydrocarbon oxidation rates. We expanded these findings by exploring metatranscriptomic signatures of the microbial communities during the chemical amendments in the microcosm experiments. Here we report how dominant organisms were uniquely associated with treatment- and time-dependent trajectories during the exposure regimes; nutrient availability was a significant factor in driving changes in metatranscriptomic responses. Remarkable signals associated with PPTEPs showed the potential role of mobilome and viral-associated survival responses. These insights underscore the time-dependent environmental perturbations of fragile marine environments under oil and anthropogenic stress., Competing Interests: The authors declare no conflict of interest.

Environmental perturbations shape the structure and function of microbial communities. Oil spills are a major perturbation and resolving spills often requires active measures like dispersant application that can exacerbate the initial disturbance. Species-specific responses of microorganisms to oil and dispersant exposure during such perturbations remain largely unknown. We merged metatranscriptomic libraries with pangenomes to generate Core-Accessory Metatranscriptomes (CA-Metatranscriptomes) for two microbial hydrocarbon degraders that played important roles in the aftermath of the Deepwater Horizon oil spill. The Colwellia CA-Metatranscriptome illustrated pronounced dispersant-driven acceleration of core (~41%) and accessory gene (~59%) transcription, suggesting an opportunistic strategy. Marinobacter responded to oil exposure by expressing mainly accessory genes (~93%), suggesting an effective hydrocarbon-degrading lifestyle. The CA-Metatranscriptome approach offers a robust way to identify the underlying mechanisms of key microbial functions and highlights differences of specialist-vs-opportunistic responses to environmental disturbance., (© 2023. ISME Publications B.V.)

Nearly half of carbon fixation and primary production originates from marine phytoplankton, and much of it occurs in episodic blooms in upwelling regimes. Here, we simulated blooms limited by nitrogen and iron by incubating Monterey Bay surface waters with subnutricline waters and inorganic nutrients and measured the whole-community transcriptomic response during mid- and late-bloom conditions. Cell counts revealed that centric and pennate diatoms (largely Pseudo-nitzschia and Chaetoceros spp.) were the major blooming taxa, but dinoflagellates, prasinophytes, and prymnesiophytes also increased. Viral mRNA significantly increased in late bloom and likely played a role in the bloom's demise. We observed conserved shifts in the genetic similarity of phytoplankton populations to cultivated strains, indicating adaptive population-level changes in community composition. Additionally, the density of single nucleotide variants (SNVs) declined in late-bloom samples for most taxa, indicating a loss of intraspecific diversity as a result of competition and a selective sweep of adaptive alleles. We noted differences between mid- and late-bloom metabolism and differential regulation of light-harvesting complexes (LHCs) under nutrient stress. While most LHCs are diminished under nutrient stress, we showed that diverse taxa upregulated specialized, energy-dissipating LHCs in low iron. We also suggest the relative expression of NRT2 compared to the expression of GSII as a marker of cellular nitrogen status and the relative expression of iron starvation-induced protein genes ( ISIP1 , ISIP2 , and ISIP3 ) compared to the expression of the thiamine biosynthesis gene ( thiC ) as a marker of iron status in natural diatom communities. IMPORTANCE Iron and nitrogen are the nutrients that most commonly limit phytoplankton growth in the world's oceans. The utilization of these resources by phytoplankton sets the biomass available to marine systems and is of particular interest in high-nutrient, low-chlorophyll (HNLC) coastal fisheries. Previous research has described the biogeography of phytoplankton in HNLC regions and the transcriptional responses of representative taxa to nutrient limitation. However, the differential transcriptional responses of whole phytoplankton communities to iron and nitrogen limitation has not been previously described, nor has the selective pressure that these competitive bloom environments exert on major players. In addition to describing changes in the physiology of diverse phytoplankton, we suggest practical indicators of cellular nitrogen and iron status for future monitoring.

Diatoms are one of the most successful phytoplankton groups in our oceans, being responsible for over 20% of the Earth's photosynthetic productivity. Their chimeric genomes have genes derived from red algae, green algae, bacteria, and heterotrophs, resulting in multiple isoenzymes targeted to different cellular compartments with the potential for differential regulation under nutrient limitation. The resulting interactions between metabolic pathways are not yet fully understood. We previously showed how acclimation to Cu limitation enhanced susceptibility to overreduction of the photosynthetic electron transport chain and its reorganization to favor photoprotection over light harvesting in the oceanic diatom Thalassiosira oceanica (Hippmann et al., 2017, 10.1371/journal.pone.0181753). In order to gain a better understanding of the overall metabolic changes that help alleviate the stress of Cu limitation, we have further analyzed the comprehensive proteomic datasets generated in that study to identify differentially expressed proteins involved in carbon, nitrogen, and oxidative stress-related metabolic pathways. Metabolic pathway analysis showed integrated responses to Cu limitation. The upregulation of ferredoxin (Fdx) was correlated with upregulation of plastidial Fdx-dependent isoenzymes involved in nitrogen assimilation as well as enzymes involved in glutathione synthesis, thus suggesting an integration of nitrogen uptake and metabolism with photosynthesis and oxidative stress resistance. The differential expression of glycolytic isoenzymes located in the chloroplast and mitochondria may enable them to channel both excess electrons and/or ATP between these compartments. An additional support for chloroplast-mitochondrial cross-talk is the increased expression of chloroplast and mitochondrial proteins involved in the proposed malate shunt under Cu limitation., Competing Interests: The Authors did not report any conflict of interest., (© 2022 The Authors. Plant Direct published by American Society of Plant Biologists and the Society for Experimental Biology and John Wiley & Sons Ltd.)

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

The Southern Ocean (SO) harbors some of the most intense phytoplankton blooms on Earth. Changes in temperature and iron availability are expected to alter the intensity of SO phytoplankton blooms, but little is known about how these changes will influence community composition and downstream biogeochemical processes. We performed light-saturated experimental manipulations on surface ocean microbial communities from McMurdo Sound in the Ross Sea to examine the effects of increased iron availability (+2 nM) and warming (+3 and +6 °C) on nutrient uptake, as well as the growth and transcriptional responses of two dominant diatoms, Fragilariopsis and Pseudo-nitzschia We found that community nutrient uptake and primary productivity were elevated under both warming conditions without iron addition (relative to ambient -0.5 °C). This effect was greater than additive under concurrent iron addition and warming. Pseudo-nitzschia became more abundant under warming without added iron (especially at 6 °C), while Fragilariopsis only became more abundant under warming in the iron-added treatments. We attribute the apparent advantage Pseudo-nitzschia shows under warming to up-regulation of iron-conserving photosynthetic processes, utilization of iron-economic nitrogen assimilation mechanisms, and increased iron uptake and storage. These data identify important molecular and physiological differences between dominant diatom groups and add to the growing body of evidence for Pseudo-nitzschia 's increasingly important role in warming SO ecosystems. This study also suggests that temperature-driven shifts in SO phytoplankton assemblages may increase utilization of the vast pool of excess nutrients in iron-limited SO surface waters and thereby influence global nutrient distribution and carbon cycling., Competing Interests: The authors declare no competing interest., (Copyright   2021 the Author(s). Published by PNAS.)

Marine microeukaryotes play a fundamental role in biogeochemical cycling through the transfer of energy to higher trophic levels and vertical carbon transport. Despite their global importance, microeukaryote physiology, nutrient metabolism and contributions to carbon cycling across offshore ecosystems are poorly characterized. Here, we observed the prevalence of dinoflagellates along a 4,600-km meridional transect extending across the central Pacific Ocean, where oligotrophic gyres meet equatorial upwelling waters rich in macronutrients yet low in dissolved iron. A combined multi-omics and geochemical analysis provided a window into dinoflagellate metabolism across the transect, indicating a continuous taxonomic dinoflagellate community that shifted its functional transcriptome and proteome as it extended from the euphotic to the mesopelagic zone. In euphotic waters, multi-omics data suggested that a combination of trophic modes were utilized, while mesopelagic metabolism was marked by cytoskeletal investments and nutrient recycling. Rearrangement in nutrient metabolism was evident in response to variable nitrogen and iron regimes across the gradient, with no associated change in community assemblage. Total dinoflagellate proteins scaled with particulate carbon export, with both elevated in equatorial waters, suggesting a link between dinoflagellate abundance and total carbon flux. Dinoflagellates employ numerous metabolic strategies that enable broad occupation of central Pacific ecosystems and play a dual role in carbon transformation through both photosynthetic fixation in the euphotic zone and remineralization in the mesopelagic zone.

Warming, eutrophication (nutrient fertilization) and brownification (increased loading of allochthonous organic matter) are three global trends impacting lake ecosystems. However, the independent and synergistic effects of resource addition and warming on autotrophic and heterotrophic microorganisms are largely unknown. In this study, we investigate the independent and interactive effects of temperature, dissolved organic carbon (DOC, both allochthonous and autochthonous) and nitrogen (N) supply, in addition to the effect of spatial variables, on the composition, richness, and evenness of prokaryotic and eukaryotic microbial communities in lakes across elevation and N deposition gradients in the Sierra Nevada mountains of California, USA. We found that both prokaryotic and eukaryotic communities are structured by temperature, terrestrial (allochthonous) DOC and latitude. Prokaryotic communities are also influenced by total and aquatic (autochthonous) DOC, while eukaryotic communities are also structured by nitrate. Additionally, increasing N availability was associated with reduced richness of prokaryotic communities, and both lower richness and evenness of eukaryotes. We did not detect any synergistic or antagonistic effects as there were no interactions among temperature and resource variables. Together, our results suggest that (a) organic and inorganic resources, temperature, and geographic location (based on latitude and longitude) independently influence lake microbial communities; and (b) increasing N supply due to atmospheric N deposition may reduce richness of both prokaryotic and eukaryotic microbes, probably by reducing niche dimensionality. Our study provides insight into abiotic processes structuring microbial communities across environmental gradients and their potential roles in material and energy fluxes within and between ecosystems., (© 2020 John Wiley & Sons Ltd.)

Diatoms are among the most globally distributed and ecologically successful organisms in the modern ocean, contributing upwards of 40% of total marine primary productivity 1,2 . By converting dissolved silicon into biogenic silica, and photosynthetically fixing carbon dioxide into particulate organic carbon, diatoms effectively couple the silicon (Si) and carbon cycles and ballast substantial vertical flux of carbon out of the euphotic zone into the mesopelagic and deep ocean 3-5 . Viruses are key players in ocean biogeochemical cycles 6,7 , yet little is known about how viral infection specifically impacts diatom populations. Here, we show that Si limitation facilitates virus infection and mortality in diatoms in the highly productive coastal waters of the California Current Ecosystem. Using metatranscriptomic analysis of cell-associated diatom viruses and targeted quantification of extracellular viruses, we found a link between Si stress and the early, active and lytic stages of viral infection. This relationship was also observed in cultures of the bloom-forming diatom Chaetoceros tenuissimus, where Si stress accelerated virus-induced mortality. Together, these findings contextualize viruses within the ecophysiological framework of Si availability and diatom-mediated biogeochemical cycling.

Diatoms outcompete other phytoplankton for nitrate, yet little is known about the mechanisms underpinning this ability. Genomes and genome-enabled studies have shown that diatoms possess unique features of nitrogen metabolism however, the implications for nutrient utilization and growth are poorly understood. Using a combination of transcriptomics, proteomics, metabolomics, fluxomics, and flux balance analysis to examine short-term shifts in nitrogen utilization in the model pennate diatom in Phaeodactylum tricornutum, we obtained a systems-level understanding of assimilation and intracellular distribution of nitrogen. Chloroplasts and mitochondria are energetically integrated at the critical intersection of carbon and nitrogen metabolism in diatoms. Pathways involved in this integration are organelle-localized GS-GOGAT cycles, aspartate and alanine systems for amino moiety exchange, and a split-organelle arginine biosynthesis pathway that clarifies the role of the diatom urea cycle. This unique configuration allows diatoms to efficiently adjust to changing nitrogen status, conferring an ecological advantage over other phytoplankton taxa.

Oceanic harmful algal blooms of Pseudo-nitzschia diatoms produce the potent mammalian neurotoxin domoic acid (DA). Despite decades of research, the molecular basis for its biosynthesis is not known. By using growth conditions known to induce DA production in Pseudo-nitzschia multiseries , we implemented transcriptome sequencing in order to identify DA biosynthesis genes that colocalize in a genomic four-gene cluster. We biochemically investigated the recombinant DA biosynthetic enzymes and linked their mechanisms to the construction of DA's diagnostic pyrrolidine skeleton, establishing a model for DA biosynthesis. Knowledge of the genetic basis for toxin production provides an orthogonal approach to bloom monitoring and enables study of environmental factors that drive oceanic DA production., (Copyright © 2018 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

In vast areas of the ocean, the scarcity of iron controls the growth and productivity of phytoplankton. Although most dissolved iron in the marine environment is complexed with organic molecules, picomolar amounts of labile inorganic iron species (labile iron) are maintained within the euphotic zone and serve as an important source of iron for eukaryotic phytoplankton and particularly for diatoms. Genome-enabled studies of labile iron utilization by diatoms have previously revealed novel iron-responsive transcripts, including the ferric iron-concentrating protein ISIP2A, but the mechanism behind the acquisition of picomolar labile iron remains unknown. Here we show that ISIP2A is a phytotransferrin that independently and convergently evolved carbonate ion-coordinated ferric iron binding. Deletion of ISIP2A disrupts high-affinity iron uptake in the diatom Phaeodactylum tricornutum, and uptake is restored by complementation with human transferrin. ISIP2A is internalized by endocytosis, and manipulation of the seawater carbonic acid system reveals a second-order dependence on the concentrations of labile iron and carbonate ions. In P. tricornutum, the synergistic interaction of labile iron and carbonate ions occurs at environmentally relevant concentrations, revealing that carbonate availability co-limits iron uptake. Phytotransferrin sequences have a broad taxonomic distribution and are abundant in marine environmental genomic datasets, suggesting that acidification-driven declines in the concentration of seawater carbonate ions will have a negative effect on this globally important eukaryotic iron acquisition mechanism.

There is an intricate interaction between iron (Fe) and copper (Cu) physiology in diatoms. However, strategies to cope with low Cu are largely unknown. This study unveils the comprehensive restructuring of the photosynthetic apparatus in the diatom Thalassiosira oceanica (CCMP1003) in response to low Cu, at the physiological and proteomic level. The restructuring results in a shift from light harvesting for photochemistry-and ultimately for carbon fixation-to photoprotection, reducing carbon fixation and oxygen evolution. The observed decreases in the physiological parameters Fv/Fm, carbon fixation, and oxygen evolution, concomitant with increases in the antennae absorption cross section (σPSII), non-photochemical quenching (NPQ) and the conversion factor (φe:C/ηPSII) are in agreement with well documented cellular responses to low Fe. However, the underlying proteomic changes due to low Cu are very different from those elicited by low Fe. Low Cu induces a significant four-fold reduction in the Cu-containing photosynthetic electron carrier plastocyanin. The decrease in plastocyanin causes a bottleneck within the photosynthetic electron transport chain (ETC), ultimately leading to substantial stoichiometric changes. Namely, 2-fold reduction in both cytochrome b6f complex (cytb6f) and photosystem II (PSII), no change in the Fe-rich PSI and a 40- and 2-fold increase in proteins potentially involved in detoxification of reactive oxygen species (ferredoxin and ferredoxin:NADP+ reductase, respectively). Furthermore, we identify 48 light harvesting complex (LHC) proteins in the publicly available genome of T. oceanica and provide proteomic evidence for 33 of these. The change in the LHC composition within the antennae in response to low Cu underlines the shift from photochemistry to photoprotection in T. oceanica (CCMP1003). Interestingly, we also reveal very significant intra-specific strain differences. Another strain of T. oceanica (CCMP 1005) requires significantly higher Cu concentrations to sustain both its maximal and minimal growth rate compared to CCMP 1003. Under low Cu, CCMP 1005 decreases its growth rate, cell size, Chla and total protein per cell. We argue that the reduction in protein per cell is the main strategy to decrease its cellular Cu requirement, as none of the other parameters tested are affected. Differences between the two strains, as well as differences between the well documented responses to low Fe and those presented here in response to low Cu are discussed.

The ecological prominence of diatoms in the ocean environment largely results from their superior competitive ability for dissolved nitrate (NO 3 - ). To investigate the cellular and genetic basis of diatom NO 3 - assimilation, we generated a knockout in the nitrate reductase gene ( NR -KO) of the model pennate diatom Phaeodactylum tricornutum In NR -KO cells, N-assimilation was abolished although NO 3 - transport remained intact. Unassimilated NO 3 - accumulated in NR -KO cells, resulting in swelling and associated changes in biochemical composition and physiology. Elevated expression of genes encoding putative vacuolar NO 3 - chloride channel transporters plus electron micrographs indicating enlarged vacuoles suggested vacuolar storage of NO 3 - Triacylglycerol concentrations in the NR -KO cells increased immediately following the addition of NO 3 - , and these increases coincided with elevated gene expression of key triacylglycerol biosynthesis components. Simultaneously, induction of transcripts encoding proteins involved in thylakoid membrane lipid recycling suggested more abrupt repartitioning of carbon resources in NR -KO cells compared with the wild type. Conversely, ribosomal structure and photosystem genes were immediately deactivated in NR -KO cells following NO 3 - addition, followed within hours by deactivation of genes encoding enzymes for chlorophyll biosynthesis and carbon fixation and metabolism. N-assimilation pathway genes respond uniquely, apparently induced simultaneously by both NO 3 - replete and deplete conditions., (© 2017 American Society of Plant Biologists. All rights reserved.)

[This corrects the article DOI: 10.1371/journal.pgen.1006490.].

Metagenomic and metatranscriptomic data were generated from size-fractionated samples from 11 sites within the Baltic Sea and adjacent marine waters of Kattegat and freshwater Lake Torneträsk in order to investigate the diversity, distribution, and transcriptional activity of virioplankton. Such a transect, spanning a salinity gradient from freshwater to the open sea, facilitated a broad genome-enabled investigation of natural as well as impacted aspects of Baltic Sea viral communities. Taxonomic signatures representative of phages within the widely distributed order Caudovirales were identified with enrichments in lesser-known families such as Podoviridae and Siphoviridae . The distribution of phage reported to infect diverse and ubiquitous heterotrophic bacteria (SAR11 clades) and cyanobacteria ( Synechococcus sp.) displayed population-level shifts in diversity. Samples from higher-salinity conditions (>14 practical salinity units [PSU]) had increased abundances of viruses for picoeukaryotes, i.e., Ostreococcus . These data, combined with host diversity estimates, suggest viral modulation of diversity on the whole-community scale, as well as in specific prokaryotic and eukaryotic lineages. RNA libraries revealed single-stranded DNA (ssDNA) and RNA viral populations throughout the Baltic Sea, with ssDNA phage highly represented in Lake Torneträsk. Further, our data suggest relatively high transcriptional activity of fish viruses within diverse families known to have broad host ranges, such as Nodoviridae (RNA), Iridoviridae (DNA), and predicted zoonotic viruses that can cause ecological and economic damage as well as impact human health. IMPORTANCE Inferred virus-host relationships, community structures of ubiquitous ecologically relevant groups, and identification of transcriptionally active populations have been achieved with our Baltic Sea study. Further, these data, highlighting the transcriptional activity of viruses, represent one of the more powerful uses of omics concerning ecosystem health. The use of omics-related data to assess ecosystem health holds great promise for rapid and relatively inexpensive determination of perturbations and risk, explicitly with regard to viral assemblages, as no single marker gene is suitable for widespread taxonomic coverage.

Environmental fluctuations affect distribution, growth and abundance of diatoms in nature, with iron (Fe) availability playing a central role. Studies on the response of diatoms to low Fe have either utilized continuous (24 hr) illumination or sampled a single time of day, missing any temporal dynamics. We profiled the physiology, metabolite composition, and global transcripts of the pennate diatom Phaeodactylum tricornutum during steady-state growth at low, intermediate, and high levels of dissolved Fe over light:dark cycles, to better understand fundamental aspects of genetic control of physiological acclimation to growth under Fe-limitation. We greatly expand the catalog of genes involved in the low Fe response, highlighting the importance of intracellular trafficking in Fe-limited diatoms. P. tricornutum exhibited transcriptomic hallmarks of slowed growth leading to prolonged periods of cell division/silica deposition, which could impact biogeochemical carbon sequestration in Fe-limited regions. Light harvesting and ribosome biogenesis transcripts were generally reduced under low Fe while transcript levels for genes putatively involved in the acquisition and recycling of Fe were increased. We also noted shifts in expression towards increased synthesis and catabolism of branched chain amino acids in P. tricornutum grown at low Fe whereas expression of genes involved in central core metabolism were relatively unaffected, indicating that essential cellular function is protected. Beyond the response of P. tricornutum to low Fe, we observed major coordinated shifts in transcript control of primary and intermediate metabolism over light:dark cycles which contribute to a new view of the significance of distinctive diatom pathways, such as mitochondrial glycolysis and the ornithine-urea cycle. This study provides new insight into transcriptional modulation of diatom physiology and metabolism across light:dark cycles in response to Fe availability, providing mechanistic understanding for the ability of diatoms to remain metabolically poised to respond quickly to Fe input and revealing strategies underlying their ecological success., Competing Interests: The authors have declared that no competing interests exist.

Background: Improvement in the performance of eukaryotic microalgae for biofuel and bioproduct production is largely dependent on characterization of metabolic mechanisms within the cell. The marine diatom Cyclotella cryptica, which was originally identified in the Aquatic Species Program, is a promising strain of microalgae for large-scale production of biofuel and bioproducts, such as omega-3 fatty acids., Results: We sequenced the nuclear genome and methylome of this oleaginous diatom to identify the genetic traits that enable substantial accumulation of triacylglycerol. The genome is comprised of highly methylated repetitive sequence, which does not significantly change under silicon starved lipid induction, and data further suggests the primary role of DNA methylation is to suppress DNA transposition. Annotation of pivotal glycolytic, lipid metabolism, and carbohydrate degradation processes reveal an expanded enzyme repertoire in  C. cryptica that would allow for an increased metabolic capacity toward triacylglycerol production. Identification of previously unidentified genes, including those involved in carbon transport and chitin metabolism, provide potential targets for genetic manipulation of carbon flux to further increase its lipid phenotype. New genetic tools were developed, bringing this organism on a par with other microalgae in terms of genetic manipulation and characterization approaches., Conclusions: Functional annotation and detailed cross-species comparison of key carbon rich processes in C. cryptica highlights the importance of enzymatic subcellular compartmentation for regulation of carbon flux, which is often overlooked in photosynthetic microeukaryotes. The availability of the genome sequence, as well as advanced genetic manipulation tools enable further development of this organism for deployment in large-scale production systems.

The Amazon River has the largest discharge of all rivers on Earth, and its complex plume system fuels a wide array of biogeochemical processes, across a large area of the western tropical North Atlantic. The plume thus stimulates microbial processes affecting carbon sequestration and nutrient cycles at a global scale. Chromosomal gene expression patterns of the 2.0 to 156 μm size-fraction eukaryotic microbial community were investigated in the Amazon River Plume, generating a robust dataset (more than 100 million mRNA sequences) that depicts the metabolic capabilities and interactions among the eukaryotic microbes. Combining classical oceanographic field measurements with metatranscriptomics yielded characterization of the hydrographic conditions simultaneous with a quantification of transcriptional activity and identity of the community. We highlight the patterns of eukaryotic gene expression for 31 biogeochemically significant gene targets hypothesized to be valuable within forecasting models. An advantage to this targeted approach is that the database of reference sequences used to identify the target genes was selectively constructed and highly curated optimizing taxonomic coverage, throughput, and the accuracy of annotations. A coastal diatom bloom highly expressed nitrate transporters and carbonic anhydrase presumably to support high growth rates and enhance uptake of low levels of dissolved nitrate and CO2. Diatom-diazotroph association (DDA: diatoms with nitrogen fixing symbionts) blooms were common when surface salinity was mesohaline and dissolved nitrate concentrations were below detection, and hence did not show evidence of nitrate utilization, suggesting they relied on ammonium transporters to aquire recently fixed nitrogen. These DDA blooms in the outer plume had rapid turnover of the photosystem D1 protein presumably caused by photodegradation under increased light penetration in clearer waters, and increased expression of silicon transporters as silicon became limiting. Expression of these genes, including carbonic anhydrase and transporters for nitrate and phosphate, were found to reflect the physiological status and biogeochemistry of river plume environments. These relatively stable patterns of eukaryotic transcript abundance occurred over modest spatiotemporal scales, with similarity observed in sample duplicates collected up to 2.45 km in space and 120 minutes in time. These results confirm the use of metatranscriptomics as a valuable tool to understand and predict microbial community function., Competing Interests: The authors have declared that no competing interests exist.

Metacaspases are distant homologs of metazoan caspase proteases, implicated in stress response, and programmed cell death (PCD) in bacteria and phytoplankton. While the few previous studies on metacaspases have relied on cultured organisms and sequenced genomes, no studies have focused on metacaspases in a natural setting. We here present data from the first microbial community-wide metacaspase survey; performed by querying metagenomic and metatranscriptomic datasets from the brackish Baltic Sea, a water body characterized by pronounced environmental gradients and periods of massive cyanobacterial blooms. Metacaspase genes were restricted to ~4% of the bacteria, taxonomically affiliated mainly to Bacteroidetes, Alpha- and Betaproteobacteria and Cyanobacteria. The gene abundance was significantly higher in larger or particle-associated bacteria (>0.8 μm), and filamentous Cyanobacteria dominated metacaspase gene expression throughout the bloom season. Distinct seasonal expression patterns were detected for the three metacaspase genes in Nodularia spumigena, one of the main bloom-formers. Clustering of normalized gene expression in combination with analyses of genomic and assembly data suggest functional diversification of these genes, and possible roles of the metacaspase genes related to stress responses, i.e., sulfur metabolism in connection to oxidative stress, and nutrient stress induced cellular differentiation. Co-expression of genes encoding metacaspases and nodularin toxin synthesis enzymes was also observed in Nodularia spumigena. The study shows that metacaspases represent an adaptation of potentially high importance for several key organisms in the Baltic Sea, most prominently Cyanobacteria, and open up for further exploration of their physiological roles in microbes and assessment of their ecological impact in aquatic habitats.

Diatoms are a dominant group of eukaryotic phytoplankton that contribute substantially to global primary production and the cycling of important elements such as carbon and nitrogen. Heterotrophic bacteria, including members of the gammaproteobacteria, are commonly associated with diatom populations and may rely on them for organic carbon while potentially competing with them for other essential nutrients. Considering that bacterioplankton drive oceanic release of CO 2 (i.e., bacterial respiration) while diatoms drive ocean carbon sequestration vial the biological pump, the outcome of such competition could influence the direction and magnitude of carbon flux in the upper ocean. Nitrate availability is commonly a determining factor for the growth of diatom populations, particularly in coastal and upwelling regions. Diatoms as well as many bacterial species can utilize nitrate, however the ability of bacteria to compete for nitrate may be hindered by carbon limitation. Here we have developed a genetically tractable model system using the pennate diatom Phaeodactylum tricornutum and the widespread heterotrophic bacteria Alteromonas macleodii to examine carbon-nitrogen dynamics. While subsisting solely on P. tricornutum derived carbon, A. macleodii does not appear to be an effective competitor for nitrate, and may in fact benefit the diatom; particularly in stationary phase. However, allochthonous dissolved organic carbon addition in the form of pyruvate triggers A. macleodii proliferation and nitrate uptake, leading to reduced P. tricornutum growth. Nitrate reductase deficient mutants of A. macleodii (ΔnasA) do not exhibit such explosive growth and associated competitive ability in response to allochthonous carbon when nitrate is the sole nitrogen source, but could survive by utilizing solely P. tricornutum-derived nitrogen. Furthermore, allocthonous carbon addition enables wild-type A. macleodii to rescue nitrate reductase deficient P. tricornutum populations from nitrogen starvation, and RNA-seq transcriptomic evidence supports nitrogen-based interactions between diatoms and bacteria at the molecular level. This study provides key insights into the roles of carbon and nitrogen in phytoplankton-bacteria dynamics and lays the foundation for developing a mechanistic understanding of these interactions using co-culturing and genetic manipulation.

Background: Prostate cancer incidence and mortality rates are significantly increased in African-American men, but limited studies have been performed within Sub-Saharan African populations. As mitochondria control energy metabolism and apoptosis we speculate that somatic mutations within mitochondrial genomes are candidate drivers of aggressive prostate carcinogenesis., Methods: We used matched blood and prostate tissue samples from 87 South African men (77 with African ancestry) to perform deep sequencing of complete mitochondrial genomes. Clinical presentation was biased toward aggressive disease (Gleason score >7, 64%), and compared with men without prostate cancer either with or without benign prostatic hyperplasia., Results: We identified 144 somatic mtDNA single nucleotide variants (SNVs), of which 80 were observed in 39 men presenting with aggressive disease. Both the number and frequency of somatic mtDNA SNVs were associated with higher pathological stage., Conclusions: Besides doubling the total number of somatic PCa-associated mitochondrial genome mutations identified to date, we associate mutational load with aggressive prostate cancer status in men of African ancestry., (© 2015 The Authors. The Prostate published by Wiley Periodicals, Inc.)

Southern Ocean primary productivity plays a key role in global ocean biogeochemistry and climate. At the Southern Ocean sea ice edge in coastal McMurdo Sound, we observed simultaneous cobalamin and iron limitation of surface water phytoplankton communities in late Austral summer. Cobalamin is produced only by bacteria and archaea, suggesting phytoplankton-bacterial interactions must play a role in this limitation. To characterize these interactions and investigate the molecular basis of multiple nutrient limitation, we examined transitions in global gene expression over short time scales, induced by shifts in micronutrient availability. Diatoms, the dominant primary producers, exhibited transcriptional patterns indicative of co-occurring iron and cobalamin deprivation. The major contributor to cobalamin biosynthesis gene expression was a gammaproteobacterial population, Oceanospirillaceae ASP10-02a. This group also contributed significantly to metagenomic cobalamin biosynthesis gene abundance throughout Southern Ocean surface waters. Oceanospirillaceae ASP10-02a displayed elevated expression of organic matter acquisition and cell surface attachment-related genes, consistent with a mutualistic relationship in which they are dependent on phytoplankton growth to fuel cobalamin production. Separate bacterial groups, including Methylophaga, appeared to rely on phytoplankton for carbon and energy sources, but displayed gene expression patterns consistent with iron and cobalamin deprivation. This suggests they also compete with phytoplankton and are important cobalamin consumers. Expression patterns of siderophore- related genes offer evidence for bacterial influences on iron availability as well. The nature and degree of this episodic colimitation appear to be mediated by a series of phytoplankton-bacterial interactions in both positive and negative feedback loops.

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.

Bacterial community composition and functional potential change subtly across gradients in the surface ocean. In contrast, while there are significant phylogenetic divergences between communities from freshwater and marine habitats, the underlying mechanisms to this phylogenetic structuring yet remain unknown. We hypothesized that the functional potential of natural bacterial communities is linked to this striking divide between microbiomes. To test this hypothesis, metagenomic sequencing of microbial communities along a 1,800 km transect in the Baltic Sea area, encompassing a continuous natural salinity gradient from limnic to fully marine conditions, was explored. Multivariate statistical analyses showed that salinity is the main determinant of dramatic changes in microbial community composition, but also of large scale changes in core metabolic functions of bacteria. Strikingly, genetically and metabolically different pathways for key metabolic processes, such as respiration, biosynthesis of quinones and isoprenoids, glycolysis and osmolyte transport, were differentially abundant at high and low salinities. These shifts in functional capacities were observed at multiple taxonomic levels and within dominant bacterial phyla, while bacteria, such as SAR11, were able to adapt to the entire salinity gradient. We propose that the large differences in central metabolism required at high and low salinities dictate the striking divide between freshwater and marine microbiomes, and that the ability to inhabit different salinity regimes evolved early during bacterial phylogenetic differentiation. These findings significantly advance our understanding of microbial distributions and stress the need to incorporate salinity in future climate change models that predict increased levels of precipitation and a reduction in salinity.

Recent advances in metagenomics have generated numerous insights related to the cell biology, ecology, evolution, and biogeochemistry of microbial life in the ocean. Notably, advances in single cell genomics (SCG), fluorescence activated cell sorting (FACS) capture of specific populations and metagenomic assembly are providing valuable information related to within and between lineage genomic content. Improved reference genome databases have assisted biogeographic studies of particular taxa, including in some cases different ecotypes. Several studies targeting picoplankton, associated with various particle size classes, have begun to define contrasting trends in gene family evolution between free-living compared to particle associated microbes. Also viruses and eukaryotic microbes are increasingly considered in metagenomic studies and specific associations between viruses, bacteria, archaea, and eukaryotic microbes are emerging.

Metagenomic data sets were generated from samples collected along a coastal to open ocean transect between Southern California Bight and California Current waters during a seasonal upwelling event, providing an opportunity to examine the impact of episodic pulses of cold nutrient-rich water into surface ocean microbial communities. The data set consists of ~5.8 million predicted proteins across seven sites, from three different size classes: 0.1-0.8, 0.8-3.0 and 3.0-200.0 μm. Taxonomic and metabolic analyses suggest that sequences from the 0.1-0.8 μm size class correlated with their position along the upwelling mosaic. However, taxonomic profiles of bacteria from the larger size classes (0.8-200 μm) were less constrained by habitat and characterized by an increase in Cyanobacteria, Bacteroidetes, Flavobacteria and double-stranded DNA viral sequences. Functional annotation of transmembrane proteins indicate that sites comprised of organisms with small genomes have an enrichment of transporters with substrate specificities for amino acids, iron and cadmium, whereas organisms with larger genomes have a higher percentage of transporters for ammonium and potassium. Eukaryotic-type glutamine synthetase (GS) II proteins were identified and taxonomically classified as viral, most closely related to the GSII in Mimivirus, suggesting that marine Mimivirus-like particles may have played a role in the transfer of GSII gene functions. Additionally, a Planctomycete bloom was sampled from one upwelling site providing a rare opportunity to assess the genomic composition of a marine Planctomycete population. The significant correlations observed between genomic properties, community structure and nutrient availability provide insights into habitat-driven dynamics among oligotrophic versus upwelled marine waters adjoining each other spatially.

The characterization of global marine microbial taxonomic and functional diversity is a primary goal of the Global Ocean Sampling Expedition. As part of this study, 19 water samples were collected aboard the Sorcerer II sailing vessel from the southern Indian Ocean in an effort to more thoroughly understand the lifestyle strategies of the microbial inhabitants of this ultra-oligotrophic region. No investigations of whole virioplankton assemblages have been conducted on waters collected from the Indian Ocean or across multiple size fractions thus far. Therefore, the goals of this study were to examine the effect of size fractionation on viral consortia structure and function and understand the diversity and functional potential of the Indian Ocean virome. Five samples were selected for comprehensive metagenomic exploration; and sequencing was performed on the microbes captured on 3.0-, 0.8- and 0.1 µm membrane filters as well as the viral fraction (<0.1 µm). Phylogenetic approaches were also used to identify predicted proteins of viral origin in the larger fractions of data from all Indian Ocean samples, which were included in subsequent metagenomic analyses. Taxonomic profiling of viral sequences suggested that size fractionation of marine microbial communities enriches for specific groups of viruses within the different size classes and functional characterization further substantiated this observation. Functional analyses also revealed a relative enrichment for metabolic proteins of viral origin that potentially reflect the physiological condition of host cells in the Indian Ocean including those involved in nitrogen metabolism and oxidative phosphorylation. A novel classification method, MGTAXA, was used to assess virus-host relationships in the Indian Ocean by predicting the taxonomy of putative host genera, with Prochlorococcus, Acanthochlois and members of the SAR86 cluster comprising the most abundant predictions. This is the first study to holistically explore virioplankton dynamics across multiple size classes and provides unprecedented insight into virus diversity, metabolic potential and virus-host interactions.

One of the big operational problems facing laboratories today is the ability to rapidly distinguish between strains of bacteria that, while physiologically distinct, are nearly impossible to separate based on 16S rRNA gene sequence differences. Here we demonstrate that ITS-DGGE provides a convenient approach to distinguishing between closely related strains of Shewanella, some of which were impossible to separate and identify by 16 rRNA gene sequence alone. Examined Shewanella genomes contain 8-11 copies of rrn (ribosomal RNA gene) operons, and variable size and sequence of 16S-23S ITS (intergenic transcribed spacer) regions which result in distinct ITS-DGGE profiles. Phylogenetic constructions based on ITS are congruent with the genomic trees generated from concatenated core genes as well as with those based on conserved indels, suggesting that ITS patterns appear to be linked with evolutionary lineages and physiology. In addition, three new Shewanella strains (MFC 2, MFC 6, and MFC 14) were isolated from microbial fuel cells enriched from wastewater sludge and identified by ITS-DGGE. Subsequent physiological and electrochemical studies of the three isolates confirmed that each strain is phenotypically/genotypically distinct. Thus, this study validates ITS-DGGE as a quick fingerprint approach to identifying and distinguishing between closely related but novel Shewanella ecotypes., (Published by Elsevier B.V.)

Diatoms dominate the biomass of phytoplankton in nutrient-rich conditions and form the basis of some of the world's most productive marine food webs. The diatom nuclear genome contains genes with bacterial and plastid origins as well as genes of the secondary endosymbiotic host (the exosymbiont), yet little is known about the relative contribution of each gene group to diatom metabolism. Here we show that the exosymbiont-derived ornithine-urea cycle, which is similar to that of metazoans but is absent in green algae and plants, facilitates rapid recovery from prolonged nitrogen limitation. RNA-interference-mediated knockdown of a mitochondrial carbamoyl phosphate synthase impairs the response of nitrogen-limited diatoms to nitrogen addition. Metabolomic analyses indicate that intermediates in the ornithine-urea cycle are particularly depleted and that both the tricarboxylic acid cycle and the glutamine synthetase/glutamate synthase cycles are linked directly with the ornithine-urea cycle. Several other depleted metabolites are generated from ornithine-urea cycle intermediates by the products of genes laterally acquired from bacteria. This metabolic coupling of bacterial- and exosymbiont-derived proteins seems to be fundamental to diatom physiology because the compounds affected include the major diatom osmolyte proline and the precursors for long-chain polyamines required for silica precipitation during cell wall formation. So far, the ornithine-urea cycle is only known for its essential role in the removal of fixed nitrogen in metazoans. In diatoms, this cycle serves as a distribution and repackaging hub for inorganic carbon and nitrogen and contributes significantly to the metabolic response of diatoms to episodic nitrogen availability. The diatom ornithine-urea cycle therefore represents a key pathway for anaplerotic carbon fixation into nitrogenous compounds that are essential for diatom growth and for the contribution of diatoms to marine productivity.

The understanding of marine microbial ecology and metabolism has been hampered by the paucity of sequenced reference genomes. To this end, we report the sequencing of 137 diverse marine isolates collected from around the world. We analysed these sequences, along with previously published marine prokaryotic genomes, in the context of marine metagenomic data, to gain insights into the ecology of the surface ocean prokaryotic picoplankton (0.1-3.0 μm size range). The results suggest that the sequenced genomes define two microbial groups: one composed of only a few taxa that are nearly always abundant in picoplanktonic communities, and the other consisting of many microbial taxa that are rarely abundant. The genomic content of the second group suggests that these microbes are capable of slow growth and survival in energy-limited environments, and rapid growth in energy-rich environments. By contrast, the abundant and cosmopolitan picoplanktonic prokaryotes for which there is genomic representation have smaller genomes, are probably capable of only slow growth and seem to be relatively unable to sense or rapidly acclimate to energy-rich conditions. Their genomic features also lead us to propose that one method used to avoid predation by viruses and/or bacterivores is by means of slow growth and the maintenance of low biomass.

Among eukaryotes, four major phytoplankton lineages are responsible for marine photosynthesis; prymnesiophytes, alveolates, stramenopiles, and prasinophytes. Contributions by individual taxa, however, are not well known, and genomes have been analyzed from only the latter two lineages. Tiny "picoplanktonic" members of the prymnesiophyte lineage have long been inferred to be ecologically important but remain poorly characterized. Here, we examine pico-prymnesiophyte evolutionary history and ecology using cultivation-independent methods. 18S rRNA gene analysis showed pico-prymnesiophytes belonged to broadly distributed uncultivated taxa. Therefore, we used targeted metagenomics to analyze uncultured pico-prymnesiophytes sorted by flow cytometry from subtropical North Atlantic waters. The data reveal a composite nuclear-encoded gene repertoire with strong green-lineage affiliations, which contrasts with the evolutionary history indicated by the plastid genome. Measured pico-prymnesiophyte growth rates were rapid in this region, resulting in primary production contributions similar to the cyanobacterium Prochlorococcus. On average, pico-prymnesiophytes formed 25% of global picophytoplankton biomass, with differing contributions in five biogeographical provinces spanning tropical to subpolar systems. Elements likely contributing to success include high gene density and genes potentially involved in defense and nutrient uptake. Our findings have implications reaching beyond pico-prymnesiophytes, to the prasinophytes and stramenopiles. For example, prevalence of putative Ni-containing superoxide dismutases (SODs), instead of Fe-containing SODs, seems to be a common adaptation among eukaryotic phytoplankton for reducing Fe quotas in low-Fe modern oceans. Moreover, highly mosaic gene repertoires, although compositionally distinct for each major eukaryotic lineage, now seem to be an underlying facet of successful marine phytoplankton.

High levels of alignment errors associated with gaps have generally meant their exclusion from phylogenetic analysis. Conserved inserts and deletions (indels) may in some cases be less subject to errors than amino acid substitutions for inferring the history of genomes and identifying recently laterally transferred genes, but alignment error near gaps must be evaluated prior to using indels as phylogenetic characters. A method is presented for evaluating the phylogenetic unambiguity of gaps in multiple sequence alignments by allowing a defined amount of pairwise alignment ambiguity. This work considers the bacterial genus Shewanella, which is of particular interest for applications of bioremediation and environmental engineering. Understanding the genetic history of these species is vital for these applications. A set of pairwise dynamic programming alignments is constructed to test positions in multiple alignments for phylogenetic unambiguity, and a whole genome scan is done on protein sequences from 11 sequenced species of the bacterial genus Shewanella. The splits defined by phylogenetically unambiguous indels are then used as characters for phylogenetic analysis, and results are compared to whole genome Maximum Likelihood phylogeny. A comparable description of the history of the species is found, as well as a set of lateral gene transfer candidates undetectable by traditional analysis of amino acid substitutions. This analysis is applicable to other taxonomic units at all levels and has the potential to allow cataloging of clear genome-wide phylogenetic markers for taxonomic profiling down to the species level.